BACKGROUND
water molds such as saprolegnia and aphanomyces species are responsible for devastating infections on fish in aquaculture, fish farms and hobby fish tanks  <cit> .
members of the genus saprolegnia cause saprolegniosis, a disease that is characterized by visible white or grey patches of filamentous mycelium on the body or fins of freshwater fish  <cit> .
the oomycete saprolegnia parasitica is economically one of the most important fish pathogens, especially on salmon and trout species.
it causes tens of million dollar losses to aquaculture business worldwide, notably in scotland, scandinavia, chile, japan, canada, and the usa  <cit> .
s. parasitica infections are second only to bacterial diseases.
in japan, there is an annual mortality rate of 50% in coho salmon and elver due to s. parasitica infections  <cit> .
in the united states, "winter kill" in catfish caused by saprolegnia results in financial loses of up to 50%, which represents an economic loss of $ <dig> million  <cit> .
in scotland, saprolegniosis also causes significant losses with the main problem occurring in salmon hatcheries.
previously, saprolegnia infections were kept under control with malachite green, an organic dye that is very efficient at killing the pathogen.
however, since  <dig> the use of malachite green has been banned around the world, due to its carcinogenic and toxicological effects.
this has resulted in dramatic recrudescence of saprolegnia infections.
therefore, there is an urgent need for novel alternative methods of management of saprolegniosis.
saprolegnia is often considered an opportunistic pathogen that is saprotrophic and necrotrophic  <cit> .
however, it has become apparent that some s. parasitica strains are highly virulent and able to cause primary infections on salmon  <cit> .
infections occur on both eggs and fish.
on eggs the disease is manifested by profuse mycelial growth on the egg surface resulting in rapid death.
on fish, saprolegnia invades epidermal tissues and can infect the entire surface of the body  <cit> .
it causes cellular necrosis as well as dermal and epidermal damage, which ultimately leads to death by heamodilution  <cit> .
severe saprolegnia infections result in lethargic behaviour, loss of equilibrium and commonly death of the fish  <cit> .
oomycete species can be pathogenic on plants, insects, crustaceans, fish, vertebrate animals, and various microorganisms  <cit> .
oomycetes, including saprolegnia, have many fungus-like characteristics, but are not true fungi.
a number of studies have indicated that they should be classified with the golden-brown algae and diatoms as stramenopiles  <cit> .
this implies that oomycetes evolved genetic and biochemical mechanisms for interaction with animals and plants that are different from those of true fungi  <cit> .
indeed, oomycetes have several clearly defined developmental stages that are not found in fungal pathogens.
for example, saprolegnia species have a complex life cycle that includes both sexual and asexual reproduction.
the asexual spore or sporangium is formed at the end of hyphal cells and can release many motile primary zoospores  <cit> .
the primary zoospores swim only for a short time before they encyst and release a secondary zoospore.
secondary zoospores are motile for a longer period and are the main infection spore  <cit> .
secondary zoospores are able to encyst and release new zoospores several times.
this process is called "polyplanetism"  <cit> , and may have evolved to allow the zoospores to have several attempts to locate and infect a host  <cit> .
uniquely within the class of oomycetes, secondary zoospores of saprolegnia can possess hairs that are thought to be required for attachment to the host  <cit> .
for example zoospores of s. parasitica have long hooked hairs that are believed to increase efficiency of the attachment to the fish hosts  <cit> .
although different in their selection of host organisms, plant and fish pathogenic oomycetes have many features in common.
evidently, the formation of specialised spore structures including zoospores, sporangia and oospores are similar.
also, infection strategies are comparable to some extent, involving encystment and attachment of zoospores on host surfaces, and penetration of host tissues.
furthermore, it is hypothesized that similar to biotrophic plant pathogenic oomycetes, such as peronospora and several phytophthora species, suppression of host defenses is likely to play a critical role in saprolegnia pathogenesis.
host defense suppression by oomycetes remains poorly understood and only a few pathogen molecules that suppress host defenses have been identified in pathogenic oomycetes  <cit> .
there is intriguing evidence that saprolegnia-infected fish appear to be immuno-compromised  <cit> .
possibly, virulence factors secreted by the pathogen might account for the immuno-suppression and the lack of an effective response to pathogen infection.
despite the huge economic importance of animal pathogenic oomycetes, such as s. parasitica, very little is known about the fundamental molecular mechanisms underlying development, pathogenicity and host specificity  <cit> .
a thorough understanding of the basic molecular processes in saprolegnia, the nature of the interactions with its hosts, and the identification of genes and proteins involved in these processes, could lead to novel control strategies that increase fish health, reduce disease losses and increase profits.
in this study we used a genomics approach to gain a first insight into the transcriptome of s. parasitica.
we generated random cdna sequences  from a cdna library of s. parasitica to identify genes that inform us about the biology of this organism and that could be involved in pathogenicity.
we provide an overview of the identified sequences as well as a detailed description of a number of notable cdnas.
the data is available through a publicly accessible website as part of the oomycete genomics database   <cit> .
RESULTS
cdna library and sequencing
we constructed a unidirectional cdna library using mrna isolated from mycelium of atcc <dig>  a s. parasitica strain isolated from diseased salmon  <cit> .
mycelium was obtained from nutrient deprived 29-day-old in vitro culture.
in other oomycetes, such as p. infestans, similar treatments promote the expression of stress-related genes and possibly mimic infection conditions  <cit> .
a total of  <dig> sequencing reactions corresponding to the 5' end of the cdna insert were performed.
of these,  <dig> gave readable sequences.
the sequences and the quality  scores were fed into ncgr's x genome initiative   <cit>  annotation pipeline and subjected to further quality controls  <cit> .
<dig> ests remained after vector and low quality sequences were removed.
of these, 5% were assessed to be in reverse orientation based on the occurrence of at least eight consecutive a residues within the first  <dig> bp.
following additional quality screening and assembly,  <dig> consensus sequences  were obtained consisting of  <dig> singletons and  <dig> consensus with two or more ests.
in total,  <dig>  bp of assembled sequences were obtained corresponding to an average consensus sequence length of  <dig> bp.
at 61%, the gc content of the assembled sequences was relatively high and similar to the gc content reported for phytophthora spp.
<cit> .
sequence annotation
the  <dig> consensus sequences were annotated using the methods implemented in the xgi pipeline .
a total of  <dig> sequences  showed significant similarities to known protein sequences  based on blastx searches,  <dig>  gave significant hits to protein motifs in the blocks+ database, and  <dig>  gave hits to the interpro protein motif database with at least one of the  <dig> algorithms implemented in interproscan.
among these, interpro database searches with hmmpfam revealed  <dig> hits.
in total,  <dig> consensus sequences  could be assigned identities based on the gene ontology  consortium .
the differences between the different analyses are expected for such bioinformatics annotations and reflect, among other things, differences in sensitivity between the various programs.
a total of  <dig> sequences were positive with the pexfinder algorithm and are candidate for carrying signal peptides and encoding extracellular proteins.
taxonomic identity of the homologs of s. parasitica cdnas
considering the classification of oomycetes as stramenopiles, it was interesting to systematically examine the taxonomic identity of the homologs of s. parasitica cdnas.
to this end, we took advantage of the availability of several eukaryotic genomes to compile a data set of  <dig> proteins covering six major phyla of eukaryotes: fungi, animals, plants, alveolates, discicristates, and heterokonts.
the data included the complete proteomes of at least one species for all phyla except for discicristates .
we used blastx to compare the  <dig> s. parasitica unigenes to these eukaryotic proteins.
in total,  <dig> sequences showed no significant hits .
of the  <dig> sequences that showed significant hits, 32%  had the top hit to a diatom protein confirming the affinity between diatoms and oomycetes as stramenopiles.
in contrast, only  <dig>  sequences had a fungal protein as a top hit.
phylogenetic analyses
we also exploited the sequence data to examine phylogenetic affinities between s. parasitica, three phytophthora spp., and the diatom thalassiosira pseudonana.
we performed reciprocal blast searches to identify a common set of conserved protein sequences between the five species.
multiple alignments of the conserved portions of  <dig> different proteins covering  <dig> amino acids were concatenated and used in phylogenetic analyses.
the obtained tree clearly supported a monophyletic relationship between the four oomycetes, and consistent with published phylogenies of phytophthora  <cit>  suggested that p. sojae and p. ramorum are more closely related .
average amino acid identity among sequences of the phytophthora spp.
was 93% .
in contrast, average amino acid identity between s. parasitica and the three phytophthora spp.
was 77% .
average amino acid identity between the diatom and the four oomycetes was  <dig> % .
most represented protein domains
as described above, searches of the interpro database with hmmpfam revealed  <dig> hits among the  <dig> consensus sequences.
we classified the interpro domains based on their abundance .
the  <dig> hits corresponded to  <dig> different interpro domains.
of these,  <dig> domains were represented three times or more .
by far the most abundant domain was ipr <dig> for eukaryotic protein kinase, which was represented  <dig> times.
other domains that function in signal transduction, such as ankyrin-repeat, g-protein, ras gtpase, myb dna binding, were also relatively abundant.
other well-represented domains comprised fungal type i cellulose-binding domain , pan domain , and papain cysteine protease family .
fungal type i cellulose binding domain
four sequences showed similarities to fungal-type i cellulose binding domain  .
three of these occurred in cdnas predicted to encode extracellular proteins.
two cdnas, sp_002_ <dig> and sp_001_ <dig>  encoded putative proteins with two cbds.
the full length sequence of cdna spm5f <dig>  was obtained .
this cdna contained an orf of  <dig> bp corresponding to a protein of  <dig> amino acids.
signalp  <cit>  analysis of the predicted protein identified a 20-amino acid signal peptide with a significant mean s value of  <dig> .
domain ipb <dig> has been mainly reported in fungi  <cit> , but also occurs in one protein from ectocarpus siliculosus virus esv- <dig> , a viral pathogen of brown algae  <cit> .
to determine the extent to which the cbd occurs in oomycetes, we performed iterative blast searches of all publicly available phytophthora sequences using the s. parasitica cbd sequences.
in total,  <dig> different sequences similar to the s. parasitica cbds were recovered from p. infestans , p. sojae , and p. ramorum .
the  <dig> oomycete cbd sequences aligned perfectly over a  <dig> amino acid region.
multiple alignments of the oomycete cbds with the well-studied cbd of cellobiohydrolase i  of trichoderma resei  <cit>  suggested that the major features of the domain are conserved in oomycetes .
amino acid residues defining the cel6a domain including the four cysteine backbone, as well as a glutamine  and the three tyrosines  that are important for binding to cellulose, were frequently conserved.
nonetheless, tyrosines, particularly tyr <dig> and tyr <dig>  were often replaced by other aromatic residues, such as tryptophane and phenylalanine, as observed for various fungal cbds, such as endoglucanase i of t. resei  <cit> .
cbel-like sequences
seven distinct cdna sequences showed similarity to the cellulose binding, elicitor and lectin-like protein , a 34-kda extracellular glycoprotein that was first isolated from phytophthora parasitica  <cit> .
five of these cbel-like cdnas were also predicted to encode extracellular proteins based on pexfinder analyses.
cbel has a dual function:  it is required for attachment to plant surfaces,  it elicits necrosis and defense gene expression in tobacco plants  <cit> .
cbel contains two regions with similarity to the pan module , a conserved domain that includes the apple domain and functions in protein-protein or protein-carbohydrate interactions  <cit> .
the similarity to cbel in the identified s. parasitica cdnas centered on the pan module regions.
three of the seven s. parasitica cdnas contained two pan-like domains while the other four had a single domain.
we used these sequences to survey the available phytophthora sequences for additional pan/cbel-like domains using iterative blast searches.
in total,  <dig> pan/cbel-like domains were identified in  <dig> putative proteins of p. infestans, p. sojae, p. ramorum, and p. parasitica suggesting that this domain is widely distributed in oomycetes .
multiple alignment of the  <dig> oomycete cbel-like domains revealed a conserved pattern centered around a conserved core of six cysteines .
glycosyl hydrolases
six cdnas showed similarity to various classes of glucanases.
one of these, sp_001_ <dig>  showed significant similarity to microbial endo- <dig> -β-glucanases , as well as high similarity to the recently described gene, piendo <dig>  from p. infestans  <cit> .
the  <dig> bp insert of a full-length cdna  corresponding to this endo- <dig> -β-glucanase was fully sequenced .
the sequence revealed a single orf of  <dig> bp encoding a predicted protein of  <dig> amino acids with 38% identity to piendo <dig>  signalp  <cit>  analysis of the predicted protein identified a 19-amino acid signal peptide with a significant mean s value of  <dig> .
previously, phylogenetic analyses of several phytophthora hydrolytic enzymes, including piendo <dig>  revealed unexpected affinity to fungal proteins  <cit> .
indeed, blastp searches of the protein encoded by spm16a <dig> showed significant similarities to fungal and bacterial proteins  but none to other eukaryotic proteins.
proteases
we found a set of  <dig> cdnas with similarity to aspartyl , serine , and cysteine  proteases among the annotated sequences of s. parasitica.
the sequence of spm3b <dig>  a full length cdna corresponding to unigene sp_004_ <dig> was obtained .
this cdna encoded a putative protein of  <dig> amino acids.
blastp searches of the merops database  <cit>  revealed significant similarity to pepsin aspartic proteases such as cathepsin d .
signalp  <cit>  analysis of the predicted protein identified a 17-amino acid signal peptide with a significant mean s value of  <dig> .
we also determined the full length sequence of cdna spm9f <dig>  .
this cdna encoded a putative protein of  <dig> amino acids with significant similarity to papain cysteine proteases .
signalp  <cit>  analysis of the predicted protein identified a 22-amino acid signal peptide with a significant mean s value of  <dig> .
blastp searches against genbank nr and the phytophthora data sets revealed that both proteases are widely distributed among eukaryotes and oomycetes.
protease inhibitors
we have also identified two s. parasitica cdnas with similarity to protease inhibitor domains of two structural classes:  kazal-like serine protease inhibitor ,  cysteine protease inhibitor .
we further analyzed these two cdnas by aligning their putative inhibitor domains to those of known protease inhibitors .
sp_001_ <dig> showed significant similarity to the kazal-like inhibitors recently described by tian et al.
<cit>  from p. infestans and other plant pathogenic oomycetes.
amino acid residues defining the kazal motif, including the six cysteine backbone, tyrosine and asparagine residues, were conserved in sp_001_ <dig> .
the predicted active site p <dig>  which is central to the specificity of kazal inhibitors  <cit> , consisted of a proline, and therefore differed from all previously reported oomycete kazal domains  <cit> .
sp_001_ <dig> is predicted to encode a secreted protein that bears the hallmark of the cystatin class of cysteine protease inhibitors including the highly conserved qxvxg motif in the first binding loop   <cit>  .
these findings suggest that secretion of protease inhibitors is a common feature of oomycetes.
thiamine biosynthetic enzyme
we identified one sequence  with significant similarity to a thiamine biosynthetic enzyme from plants  and fungi .
unlike s. parasitica and other oomycetes, members of the genus phytophthora are thiamine auxotrophs, they require exogenous sources of thiamine for growth  <cit> .
interestingly, blast searches of the genome sequence reads of p. sojae and p. ramorum, as well as all available sequences of p. infestans, failed to reveal sequences with similarity to sp_001_ <dig> or to the plant and fungal enzymes.
these findings suggest that this thiamine biosynthetic enzyme may have been lost in the phytophthora lineage and could be related to thiamine auxotrophy in this genus.
discussion
in this study we generated  <dig> high quality ests from s. parasitica, an economically important and reemerging oomycete pathogen that causes multimillion dollar losses in the aquaculture industry.
the ests were generated from a cdna library constructed from one-month old nutrient deprived mycelium cultures.
so far significant data sets of oomycete ests have been described for three plant pathogenic species, p. infestans  <cit> , p. sojae  <cit>  and p. parasitica  <cit> .
therefore, the s. parasitica ests offer some insights into the transcriptome of animal pathogenic oomycetes, which have been extremely understudied.
prior to this work, only  <dig> nucleotide and  <dig> protein sequences of s. parasitica could be retrieved from genbank .
the sequence data and the corresponding annotations described in this study are accessible through an interactive public resource, the oomycete genomics database .
we hope that this pilot genomics project will accelerate research on this important pathogen and lays the foundation for more significant genome and cdna sequencing initiatives of animal pathogenic oomycetes.
we used the s. parasitica ests to confirm the phylogenetic affinities between oomycetes and diatoms  <cit> .
about 32% of the s. parasitica sequences that showed significant similarities to eukaryotic proteins matched a protein of the diatom t. pseudonana as a top hit.
within the oomycetes, s. parasitica is classified with other water molds, such as achlya and aphanomyces, in the order saprolegniales  <cit> .
these species are morphologically very distinct from the great majority of plant pathogens, such as the peronosporales phytophthora and downy mildews, or the pythiales pythium  <cit> .
indeed, the s. parasitica sequences tended to be relatively divergent from phytophthora sequences.
for example, based on the sequence alignments of  <dig> different conserved proteins, the average amino acid identity between s. parasitica and three phytophthora spp., p. infestans, p. sojae, and p. ramorum, was 77% compared to 93% within phytophthora.
differences in transcript content were also noted.
cdnas with similarity to elicitins, a group of 10-kda proteins that occur in all phytophthora and some pythium species  <cit>  and form 1–2% of mycelial ests in phytophthora  <cit> , were not found in the s. parasitica dataset.
although, elicitin-like genes could very well occur in the genome of s. parasitica, they do not seem to be abundantly expressed in mycelium.
elicitins were shown to function as sterol carriers  <cit> .
among the oomycetes, members of the saprolegniales are able to synthesize sterols de novo whereas phytophthora and pythium spp.
are sterol auxotrophs  <cit> .
possibly, s. parasitica may not require sterol carriers, such as elicitins, for optimal hyphal growth.
another difference between phytophthora and saprolegnia involves thiamine metabolism.
members of the genus phytophthora are thiamine auxotrophs and require exogenous sources of thiamine for growth  <cit> .
we identified one sequence in s. parasitica  that shows significant similarity to a thiamine biosynthetic enzyme from plants and fungi but that is absent in the draft genome sequences of p. sojae and p. ramorum.
this finding suggests that this thiamine biosynthetic enzyme may have been lost in the phytophthora lineage and could be related to thiamine auxotrophy in this genus.
we searched the annotated data set for s. parasitica sequences that show similarities to known proteins and protein motifs that could inform us about the biology and pathology of this microbe.
about half of the unigenes showed similarities to known protein sequences and could be assigned a putative function.
a number of sequences showed particularly interesting similarities.
cdnas with similarity to signal transduction proteins, such as kinases and transcription factors, were particularly abundant.
in total,  <dig> cdnas encoded proteins with a putative signal peptide that are potentially secreted to the extracellular space.
secretion is an essential mechanism for delivery of virulence factors by eukaryotic pathogens to their appropriate site in infected host tissue.
therefore, several putative secreted proteins of s. parasitica, such as cbd proteins, cbel-like proteins, glycosyl hydrolases, proteases, and protease inhibitors could function in virulence and will be worthy of additional studies.
phylogenetic analyses indicated that several phytophthora proteins, particularly hydrolytic enzymes such endopolygalacturonases, pectate lyases, exo- <dig> -beta-glucanases, and an endo- <dig> -beta-glucanase, are more similar to fungal proteins than to their counterparts in other eukaryotes  <cit> .
these observations are in sharp contrast with phylogenies constructed from ribosomal sequences or compiled protein sequences from mitochondrial and housekeeping chromosomal genes, which indicate considerable evolutionary distance between oomycetes and fungi  <cit> .
the apparent discrepancies between these phylogenies could reflect convergent evolution in the arsenal of hydrolytic enzymes between these pathogens, perhaps as a result of common mechanisms of infection among filamentous microbes  <cit> .
the s. parasitica sequences allowed us to evaluate whether the similarity to fungal proteins extends to oomycetes other than phytophthora and to animal pathogenic oomycetes.
although no s. parasitica sequences similar to the cell wall degrading enzymes endopolygalacturonases and pectate lyases were found, we identified a cdna, spm9f <dig>  that encodes a  <dig> amino acid protein with high similarity to endo- <dig> -β-glucanases, including the recently described piendo <dig> of p. infestans  <cit> .
similar to piendo <dig>  spm9f <dig> was most similar to fungal glucanases and no significant blastp hits  were observed to non-fungal eukaryotic proteins.
therefore, conservation in the arsenal of hydrolytic enzymes appears to extend beyond phytophthora spp.
to the saprolegniales and animal pathogenic oomycetes.
domain annotation of the s. parasitica sequences revealed the occurrence of a protein domain typically associated with fungi.
type i cbds  are thought to be unique to fungi  <cit> , although a related domain also occurs in the brown algae viral pathogen ectocarpus siliculosus virus esv- <dig>   <cit> .
in this study, we found that this domain is widespread and diverse in s. parasitica and other oomycetes.
a total of  <dig> domains from four oomycete species were found to share a  <dig> amino acid region that aligns perfectly with the canonical t. resei cel6a cbd highlighting a core of conserved four cysteines and aromatic residues known to bind the cellulosic substrate  <cit> .
interestingly, the occurrence of this cbd in a virus of brown algae, which are related to oomycetes, suggests that type i cbds might be more widespread in stramenopiles although we did not detect them in the draft genome sequence of the diatom t. pseudonana.
in fungi, the cbds are usually located in the n- or c-terminal regions of hydrolytic enzymes, such as cellulases and xylanases, and function by concentrating the catalytic domains on the surface of the insoluble cellulose substrate  <cit> .
one of the s. parasitica cdnas, spm5f <dig>  encodes a small  <dig> amino acid protein with two cbds.
such a protein could function as a scaffolding component of the multienzyme complex known as cellulosome  <cit> .
the function of this and other cbd proteins in s. parasitica may relate to attachment to organic debris on the host surface or during saprophytic growth.
alternatively, since cellulose is a major component of the cell wall of oomycetes, these proteins may play endogenous function in cell wall biogenesis.
seven s. parasitica cdnas showed similarity to cbel, a 34-kda cell wall glycoprotein of p. parasitica that binds to cellulose and host surfaces, functions in the agglutination of red blood cells, and elicits necrosis and defense gene expression in tobacco  <cit> .
the similarity centered mainly on two regions of cbel that match the pan module/apple domain .
the cbel-like pan module, which is thought to function in protein-protein or protein-carbohydrate interactions  <cit> , appeared to be particularly diverse in oomycetes with  <dig> different sequences identified in five species.
the pan module was found in proteins with diverse functions, such as the blood coagulation factor xi and the plasma protein pre-kallikrein  <cit> .
recently, several secreted proteins from apicomplexan mammalian parasites were found to contain apple-like domains and are thought to play a role during parasite attachment and invasion of host cells  <cit> .
for example, mic <dig>  an adhesin secreted by the apicomplexan toxoplasma gondi, contains six apple domains  <cit> .
it remains to be determined the extent to which the secreted pan/cbel-like proteins of s. parasitica play a role in attachment and invasion during interaction with the fish host.
nonetheless, it appears that in oomycetes, similar to the apicomplexan parasites, some adhesins are secreted pan module proteins.
proteolytic enzymes are considered important virulence factors that aid in host colonization and release of nutrients by animal pathogenic microbes.
it has long been known that the saprolegnia spp.
pathogenic on fish exhibit significant extracellular protease activity and it was postulated that this enzymatic activity contributes to pathogenesis  <cit> .
a serine protease gene, aasp <dig> from the related crayfish pathogen aphanomyces astacus, was recently characterized and shown to be highly expressed during in vivo growth  <cit> .
however, besides aasp <dig>  genes for secreted proteases of animal pathogenic oomycetes have not been reported.
in this study, we identified a diverse set of  <dig> cdnas of s. parasitica with similarity to the major catalytic classes of proteases.
a number of the identified proteases had a signal peptide that would predict them to be localised at the interface between pathogen and host and suggests that they are candidate virulence factors.
tian et al.
<cit>  recently reported that plant pathogenic oomycetes secrete a diverse family of kazal-like serine protease inhibitors with at least  <dig> members identified from p. infestans, p. sojae, p. ramorum, p. brassicae, and the downy mildew plasmopara halstedii.
among these, the two-domain epi <dig> protein and the three domain epi <dig> of p. infestans were found to inhibit and interact with p69b, a defense subtilase of tomato, and were suggested to play a role in counterdefense  <cit> .
inhibitors of serine protease might be ubiquitous among eukaryotic parasites.
for instance, the apicomplexan obligate parasite toxoplasma gondii secretes tgpi- <dig> and tgpi- <dig>  four-domain serine protease inhibitors of the kazal family  <cit> , and the intestinal hookworm ancylostoma ceylanicum secretes an 8-kda broad spectrum serine protease inhibitor of the kunitz family  <cit> .
here we found that kazal-like motifs also occur in saprolegniales proteins.
in addition to kazal-like motifs, we also discovered a cdna that encodes a secreted protein with similarity to the cystatin class of cysteine protease inhibitors  <cit> .
cysteine protease inhibitors, such as chagasin, have been reported in animal parasites, mainly trypanosomids, and are thought to target proteases of the insect vector or the mammalian host  <cit> .
perhaps, inhibition of host proteases is a widespread counterdefense strategy in animal and plant pathogenic eukaryotes.
future studies will help to address whether the discovered protease inhibitors play a role in s. parasitica-fish interactions.
CONCLUSIONS
this pilot cdna sequencing project provides a first look into the gene content of s. parasitica and sets the basis for genomics research in this reemerging animal pathogen.
annotation of the ests revealed a number of genes that could function in virulence.
future work will focus on developing molecular tools for functional analysis of s. parasitica genes.
in this regards, stable transformation of saprolegnia monoica has been reported  <cit> , and the rnai protocol recently developed for p. infestans  <cit>  should be adaptable to s. parasitica.
gene expression profiling will also be applied to investigate transcriptome changes during s. parasitica-fish interactions.
overall, these resources will greatly accelerate research on this important pathogen and could lead to novel perspectives for controlling saprolegniosis.

 water molds such as saprolegnia and aphanomyces species are responsible for devastating infections on fish in aquaculture, fish farms and hobby fish tanks. members of the genus saprolegnia cause saprolegniosis, a disease that is characterized by visible white or grey patches of filamentous mycelium on the body or fins of freshwater fish. the oomycete saprolegnia parasitica is economically one of the most important fish pathogens, especially on salmon and trout species. it causes tens of million dollar losses to aquaculture business worldwide, notably in scotland, scandinavia, chile, japan, canada, and the usa. s. parasitica infections are second only to bacterial diseases. in japan, there is an annual mortality rate of 50% in coho salmon and elver due to s. parasitica infections. in the united states, "winter kill" in catfish caused by saprolegnia results in financial loses of up to 50%, which represents an economic loss of $ million. in scotland, saprolegniosis also causes significant losses with the main problem occurring in salmon hatcheries. previously, saprolegnia infections were kept under control with malachite green, an organic dye that is very efficient at killing the pathogen. however, since the use of malachite green has been banned around the world, due to its carcinogenic and toxicological effects. this has resulted in dramatic recrudescence of saprolegnia infections. therefore, there is an urgent need for novel alternative methods of management of saprolegniosis. saprolegnia is often considered an opportunistic pathogen that is saprotrophic and necrotrophic. however, it has become apparent that some s. parasitica strains are highly virulent and able to cause primary infections on salmon. infections occur on both eggs and fish. on eggs the disease is manifested by profuse mycelial growth on the egg surface resulting in rapid death. on fish, saprolegnia invades epidermal tissues and can infect the entire surface of the body. it causes cellular necrosis as well as dermal and epidermal damage, which ultimately leads to death by heamodilution. severe saprolegnia infections result in lethargic behaviour, loss of equilibrium and commonly death of the fish. oomycete species can be pathogenic on plants, insects, crustaceans, fish, vertebrate animals, and various microorganisms. oomycetes, including saprolegnia, have many fungus-like characteristics, but are not true fungi. a number of studies have indicated that they should be classified with the golden-brown algae and diatoms as stramenopiles. this implies that oomycetes evolved genetic and biochemical mechanisms for interaction with animals and plants that are different from those of true fungi. indeed, oomycetes have several clearly defined developmental stages that are not found in fungal pathogens. for example, saprolegnia species have a complex life cycle that includes both sexual and asexual reproduction. the asexual spore or sporangium is formed at the end of hyphal cells and can release many motile primary zoospores. the primary zoospores swim only for a short time before they encyst and release a secondary zoospore. secondary zoospores are motile for a longer period and are the main infection spore. secondary zoospores are able to encyst and release new zoospores several times. this process is called "polyplanetism", and may have evolved to allow the zoospores to have several attempts to locate and infect a host. uniquely within the class of oomycetes, secondary zoospores of saprolegnia can possess hairs that are thought to be required for attachment to the host. for example zoospores of s. parasitica have long hooked hairs that are believed to increase efficiency of the attachment to the fish hosts. although different in their selection of host organisms, plant and fish pathogenic oomycetes have many features in common. evidently, the formation of specialised spore structures including zoospores, sporangia and oospores are similar. also, infection strategies are comparable to some extent, involving encystment and attachment of zoospores on host surfaces, and penetration of host tissues. furthermore, it is hypothesized that similar to biotrophic plant pathogenic oomycetes, such as peronospora and several phytophthora species, suppression of host defenses is likely to play a critical role in saprolegnia pathogenesis. host defense suppression by oomycetes remains poorly understood and only a few pathogen molecules that suppress host defenses have been identified in pathogenic oomycetes. there is intriguing evidence that saprolegnia-infected fish appear to be immuno-compromised. possibly, virulence factors secreted by the pathogen might account for the immuno-suppression and the lack of an effective response to pathogen infection. despite the huge economic importance of animal pathogenic oomycetes, such as s. parasitica, very little is known about the fundamental molecular mechanisms underlying development, pathogenicity and host specificity. a thorough understanding of the basic molecular processes in saprolegnia, the nature of the interactions with its hosts, and the identification of genes and proteins involved in these processes, could lead to novel control strategies that increase fish health, reduce disease losses and increase profits. in this study we used a genomics approach to gain a first insight into the transcriptome of s. parasitica. we generated random cdna sequences from a cdna library of s. parasitica to identify genes that inform us about the biology of this organism and that could be involved in pathogenicity. we provide an overview of the identified sequences as well as a detailed description of a number of notable cdnas. the data is available through a publicly accessible website as part of the oomycete genomics database.  cdna library and sequencing we constructed a unidirectional cdna library using mrna isolated from mycelium of atcc a s. parasitica strain isolated from diseased salmon. mycelium was obtained from nutrient deprived 29-day-old in vitro culture. in other oomycetes, such as p. infestans, similar treatments promote the expression of stress-related genes and possibly mimic infection conditions. a total of sequencing reactions corresponding to the 5' end of the cdna insert were performed. of these, gave readable sequences. the sequences and the quality scores were fed into ncgr's x genome initiative annotation pipeline and subjected to further quality controls. ests remained after vector and low quality sequences were removed. of these, 5% were assessed to be in reverse orientation based on the occurrence of at least eight consecutive a residues within the first bp. following additional quality screening and assembly, consensus sequences were obtained consisting of singletons and consensus with two or more ests. in total, bp of assembled sequences were obtained corresponding to an average consensus sequence length of bp. at 61%, the gc content of the assembled sequences was relatively high and similar to the gc content reported for phytophthora spp.. sequence annotation the consensus sequences were annotated using the methods implemented in the xgi pipeline. a total of sequences showed significant similarities to known protein sequences based on blastx searches, gave significant hits to protein motifs in the blocks+ database, and gave hits to the interpro protein motif database with at least one of the algorithms implemented in interproscan. among these, interpro database searches with hmmpfam revealed hits. in total, consensus sequences could be assigned identities based on the gene ontology consortium. the differences between the different analyses are expected for such bioinformatics annotations and reflect, among other things, differences in sensitivity between the various programs. a total of sequences were positive with the pexfinder algorithm and are candidate for carrying signal peptides and encoding extracellular proteins. taxonomic identity of the homologs of s. parasitica cdnas considering the classification of oomycetes as stramenopiles, it was interesting to systematically examine the taxonomic identity of the homologs of s. parasitica cdnas. to this end, we took advantage of the availability of several eukaryotic genomes to compile a data set of proteins covering six major phyla of eukaryotes: fungi, animals, plants, alveolates, discicristates, and heterokonts. the data included the complete proteomes of at least one species for all phyla except for discicristates. we used blastx to compare the s. parasitica unigenes to these eukaryotic proteins. in total, sequences showed no significant hits. of the sequences that showed significant hits, 32% had the top hit to a diatom protein confirming the affinity between diatoms and oomycetes as stramenopiles. in contrast, only sequences had a fungal protein as a top hit. phylogenetic analyses we also exploited the sequence data to examine phylogenetic affinities between s. parasitica, three phytophthora spp., and the diatom thalassiosira pseudonana. we performed reciprocal blast searches to identify a common set of conserved protein sequences between the five species. multiple alignments of the conserved portions of different proteins covering amino acids were concatenated and used in phylogenetic analyses. the obtained tree clearly supported a monophyletic relationship between the four oomycetes, and consistent with published phylogenies of phytophthora suggested that p. sojae and p. ramorum are more closely related. average amino acid identity among sequences of the phytophthora spp. was 93%. in contrast, average amino acid identity between s. parasitica and the three phytophthora spp. was 77%. average amino acid identity between the diatom and the four oomycetes was %. most represented protein domains as described above, searches of the interpro database with hmmpfam revealed hits among the consensus sequences. we classified the interpro domains based on their abundance. the hits corresponded to different interpro domains. of these, domains were represented three times or more. by far the most abundant domain was ipr for eukaryotic protein kinase, which was represented times. other domains that function in signal transduction, such as ankyrin-repeat, g-protein, ras gtpase, myb dna binding, were also relatively abundant. other well-represented domains comprised fungal type i cellulose-binding domain, pan domain, and papain cysteine protease family. fungal type i cellulose binding domain four sequences showed similarities to fungal-type i cellulose binding domain. three of these occurred in cdnas predicted to encode extracellular proteins. two cdnas, sp_002_ and sp_001_ encoded putative proteins with two cbds. the full length sequence of cdna spm5f was obtained. this cdna contained an orf of bp corresponding to a protein of amino acids. signalp analysis of the predicted protein identified a 20-amino acid signal peptide with a significant mean s value of. domain ipb has been mainly reported in fungi, but also occurs in one protein from ectocarpus siliculosus virus esv-, a viral pathogen of brown algae. to determine the extent to which the cbd occurs in oomycetes, we performed iterative blast searches of all publicly available phytophthora sequences using the s. parasitica cbd sequences. in total, different sequences similar to the s. parasitica cbds were recovered from p. infestans, p. sojae, and p. ramorum. the oomycete cbd sequences aligned perfectly over a amino acid region. multiple alignments of the oomycete cbds with the well-studied cbd of cellobiohydrolase i of trichoderma resei suggested that the major features of the domain are conserved in oomycetes. amino acid residues defining the cel6a domain including the four cysteine backbone, as well as a glutamine and the three tyrosines that are important for binding to cellulose, were frequently conserved. nonetheless, tyrosines, particularly tyr and tyr were often replaced by other aromatic residues, such as tryptophane and phenylalanine, as observed for various fungal cbds, such as endoglucanase i of t. resei. cbel-like sequences seven distinct cdna sequences showed similarity to the cellulose binding, elicitor and lectin-like protein, a 34-kda extracellular glycoprotein that was first isolated from phytophthora parasitica. five of these cbel-like cdnas were also predicted to encode extracellular proteins based on pexfinder analyses. cbel has a dual function: it is required for attachment to plant surfaces, it elicits necrosis and defense gene expression in tobacco plants. cbel contains two regions with similarity to the pan module, a conserved domain that includes the apple domain and functions in protein-protein or protein-carbohydrate interactions. the similarity to cbel in the identified s. parasitica cdnas centered on the pan module regions. three of the seven s. parasitica cdnas contained two pan-like domains while the other four had a single domain. we used these sequences to survey the available phytophthora sequences for additional pan/cbel-like domains using iterative blast searches. in total, pan/cbel-like domains were identified in putative proteins of p. infestans, p. sojae, p. ramorum, and p. parasitica suggesting that this domain is widely distributed in oomycetes. multiple alignment of the oomycete cbel-like domains revealed a conserved pattern centered around a conserved core of six cysteines. glycosyl hydrolases six cdnas showed similarity to various classes of glucanases. one of these, sp_001_ showed significant similarity to microbial endo- -β-glucanases, as well as high similarity to the recently described gene, piendo from p. infestans. the bp insert of a full-length cdna corresponding to this endo- -β-glucanase was fully sequenced. the sequence revealed a single orf of bp encoding a predicted protein of amino acids with 38% identity to piendo signalp analysis of the predicted protein identified a 19-amino acid signal peptide with a significant mean s value of. previously, phylogenetic analyses of several phytophthora hydrolytic enzymes, including piendo revealed unexpected affinity to fungal proteins. indeed, blastp searches of the protein encoded by spm16a showed significant similarities to fungal and bacterial proteins but none to other eukaryotic proteins. proteases we found a set of cdnas with similarity to aspartyl, serine, and cysteine proteases among the annotated sequences of s. parasitica. the sequence of spm3b a full length cdna corresponding to unigene sp_004_ was obtained. this cdna encoded a putative protein of amino acids. blastp searches of the merops database revealed significant similarity to pepsin aspartic proteases such as cathepsin d. signalp analysis of the predicted protein identified a 17-amino acid signal peptide with a significant mean s value of. we also determined the full length sequence of cdna spm9f. this cdna encoded a putative protein of amino acids with significant similarity to papain cysteine proteases. signalp analysis of the predicted protein identified a 22-amino acid signal peptide with a significant mean s value of. blastp searches against genbank nr and the phytophthora data sets revealed that both proteases are widely distributed among eukaryotes and oomycetes. protease inhibitors we have also identified two s. parasitica cdnas with similarity to protease inhibitor domains of two structural classes: kazal-like serine protease inhibitor, cysteine protease inhibitor. we further analyzed these two cdnas by aligning their putative inhibitor domains to those of known protease inhibitors. sp_001_ showed significant similarity to the kazal-like inhibitors recently described by tian et al. from p. infestans and other plant pathogenic oomycetes. amino acid residues defining the kazal motif, including the six cysteine backbone, tyrosine and asparagine residues, were conserved in sp_001_. the predicted active site p which is central to the specificity of kazal inhibitors, consisted of a proline, and therefore differed from all previously reported oomycete kazal domains. sp_001_ is predicted to encode a secreted protein that bears the hallmark of the cystatin class of cysteine protease inhibitors including the highly conserved qxvxg motif in the first binding loop. these findings suggest that secretion of protease inhibitors is a common feature of oomycetes. thiamine biosynthetic enzyme we identified one sequence with significant similarity to a thiamine biosynthetic enzyme from plants and fungi. unlike s. parasitica and other oomycetes, members of the genus phytophthora are thiamine auxotrophs, they require exogenous sources of thiamine for growth. interestingly, blast searches of the genome sequence reads of p. sojae and p. ramorum, as well as all available sequences of p. infestans, failed to reveal sequences with similarity to sp_001_ or to the plant and fungal enzymes. these findings suggest that this thiamine biosynthetic enzyme may have been lost in the phytophthora lineage and could be related to thiamine auxotrophy in this genus. discussion in this study we generated high quality ests from s. parasitica, an economically important and reemerging oomycete pathogen that causes multimillion dollar losses in the aquaculture industry. the ests were generated from a cdna library constructed from one-month old nutrient deprived mycelium cultures. so far significant data sets of oomycete ests have been described for three plant pathogenic species, p. infestans, p. sojae and p. parasitica. therefore, the s. parasitica ests offer some insights into the transcriptome of animal pathogenic oomycetes, which have been extremely understudied. prior to this work, only nucleotide and protein sequences of s. parasitica could be retrieved from genbank. the sequence data and the corresponding annotations described in this study are accessible through an interactive public resource, the oomycete genomics database. we hope that this pilot genomics project will accelerate research on this important pathogen and lays the foundation for more significant genome and cdna sequencing initiatives of animal pathogenic oomycetes. we used the s. parasitica ests to confirm the phylogenetic affinities between oomycetes and diatoms. about 32% of the s. parasitica sequences that showed significant similarities to eukaryotic proteins matched a protein of the diatom t. pseudonana as a top hit. within the oomycetes, s. parasitica is classified with other water molds, such as achlya and aphanomyces, in the order saprolegniales. these species are morphologically very distinct from the great majority of plant pathogens, such as the peronosporales phytophthora and downy mildews, or the pythiales pythium. indeed, the s. parasitica sequences tended to be relatively divergent from phytophthora sequences. for example, based on the sequence alignments of different conserved proteins, the average amino acid identity between s. parasitica and three phytophthora spp., p. infestans, p. sojae, and p. ramorum, was 77% compared to 93% within phytophthora. differences in transcript content were also noted. cdnas with similarity to elicitins, a group of 10-kda proteins that occur in all phytophthora and some pythium species and form 1–2% of mycelial ests in phytophthora, were not found in the s. parasitica dataset. although, elicitin-like genes could very well occur in the genome of s. parasitica, they do not seem to be abundantly expressed in mycelium. elicitins were shown to function as sterol carriers. among the oomycetes, members of the saprolegniales are able to synthesize sterols de novo whereas phytophthora and pythium spp. are sterol auxotrophs. possibly, s. parasitica may not require sterol carriers, such as elicitins, for optimal hyphal growth. another difference between phytophthora and saprolegnia involves thiamine metabolism. members of the genus phytophthora are thiamine auxotrophs and require exogenous sources of thiamine for growth. we identified one sequence in s. parasitica that shows significant similarity to a thiamine biosynthetic enzyme from plants and fungi but that is absent in the draft genome sequences of p. sojae and p. ramorum. this finding suggests that this thiamine biosynthetic enzyme may have been lost in the phytophthora lineage and could be related to thiamine auxotrophy in this genus. we searched the annotated data set for s. parasitica sequences that show similarities to known proteins and protein motifs that could inform us about the biology and pathology of this microbe. about half of the unigenes showed similarities to known protein sequences and could be assigned a putative function. a number of sequences showed particularly interesting similarities. cdnas with similarity to signal transduction proteins, such as kinases and transcription factors, were particularly abundant. in total, cdnas encoded proteins with a putative signal peptide that are potentially secreted to the extracellular space. secretion is an essential mechanism for delivery of virulence factors by eukaryotic pathogens to their appropriate site in infected host tissue. therefore, several putative secreted proteins of s. parasitica, such as cbd proteins, cbel-like proteins, glycosyl hydrolases, proteases, and protease inhibitors could function in virulence and will be worthy of additional studies. phylogenetic analyses indicated that several phytophthora proteins, particularly hydrolytic enzymes such endopolygalacturonases, pectate lyases, exo- -beta-glucanases, and an endo- -beta-glucanase, are more similar to fungal proteins than to their counterparts in other eukaryotes. these observations are in sharp contrast with phylogenies constructed from ribosomal sequences or compiled protein sequences from mitochondrial and housekeeping chromosomal genes, which indicate considerable evolutionary distance between oomycetes and fungi. the apparent discrepancies between these phylogenies could reflect convergent evolution in the arsenal of hydrolytic enzymes between these pathogens, perhaps as a result of common mechanisms of infection among filamentous microbes. the s. parasitica sequences allowed us to evaluate whether the similarity to fungal proteins extends to oomycetes other than phytophthora and to animal pathogenic oomycetes. although no s. parasitica sequences similar to the cell wall degrading enzymes endopolygalacturonases and pectate lyases were found, we identified a cdna, spm9f that encodes a amino acid protein with high similarity to endo- -β-glucanases, including the recently described piendo of p. infestans. similar to piendo spm9f was most similar to fungal glucanases and no significant blastp hits were observed to non-fungal eukaryotic proteins. therefore, conservation in the arsenal of hydrolytic enzymes appears to extend beyond phytophthora spp. to the saprolegniales and animal pathogenic oomycetes. domain annotation of the s. parasitica sequences revealed the occurrence of a protein domain typically associated with fungi. type i cbds are thought to be unique to fungi, although a related domain also occurs in the brown algae viral pathogen ectocarpus siliculosus virus esv-. in this study, we found that this domain is widespread and diverse in s. parasitica and other oomycetes. a total of domains from four oomycete species were found to share a amino acid region that aligns perfectly with the canonical t. resei cel6a cbd highlighting a core of conserved four cysteines and aromatic residues known to bind the cellulosic substrate. interestingly, the occurrence of this cbd in a virus of brown algae, which are related to oomycetes, suggests that type i cbds might be more widespread in stramenopiles although we did not detect them in the draft genome sequence of the diatom t. pseudonana. in fungi, the cbds are usually located in the n- or c-terminal regions of hydrolytic enzymes, such as cellulases and xylanases, and function by concentrating the catalytic domains on the surface of the insoluble cellulose substrate. one of the s. parasitica cdnas, spm5f encodes a small amino acid protein with two cbds. such a protein could function as a scaffolding component of the multienzyme complex known as cellulosome. the function of this and other cbd proteins in s. parasitica may relate to attachment to organic debris on the host surface or during saprophytic growth. alternatively, since cellulose is a major component of the cell wall of oomycetes, these proteins may play endogenous function in cell wall biogenesis. seven s. parasitica cdnas showed similarity to cbel, a 34-kda cell wall glycoprotein of p. parasitica that binds to cellulose and host surfaces, functions in the agglutination of red blood cells, and elicits necrosis and defense gene expression in tobacco. the similarity centered mainly on two regions of cbel that match the pan module/apple domain. the cbel-like pan module, which is thought to function in protein-protein or protein-carbohydrate interactions, appeared to be particularly diverse in oomycetes with different sequences identified in five species. the pan module was found in proteins with diverse functions, such as the blood coagulation factor xi and the plasma protein pre-kallikrein. recently, several secreted proteins from apicomplexan mammalian parasites were found to contain apple-like domains and are thought to play a role during parasite attachment and invasion of host cells. for example, mic an adhesin secreted by the apicomplexan toxoplasma gondi, contains six apple domains. it remains to be determined the extent to which the secreted pan/cbel-like proteins of s. parasitica play a role in attachment and invasion during interaction with the fish host. nonetheless, it appears that in oomycetes, similar to the apicomplexan parasites, some adhesins are secreted pan module proteins. proteolytic enzymes are considered important virulence factors that aid in host colonization and release of nutrients by animal pathogenic microbes. it has long been known that the saprolegnia spp. pathogenic on fish exhibit significant extracellular protease activity and it was postulated that this enzymatic activity contributes to pathogenesis. a serine protease gene, aasp from the related crayfish pathogen aphanomyces astacus, was recently characterized and shown to be highly expressed during in vivo growth. however, besides aasp genes for secreted proteases of animal pathogenic oomycetes have not been reported. in this study, we identified a diverse set of cdnas of s. parasitica with similarity to the major catalytic classes of proteases. a number of the identified proteases had a signal peptide that would predict them to be localised at the interface between pathogen and host and suggests that they are candidate virulence factors. tian et al. recently reported that plant pathogenic oomycetes secrete a diverse family of kazal-like serine protease inhibitors with at least members identified from p. infestans, p. sojae, p. ramorum, p. brassicae, and the downy mildew plasmopara halstedii. among these, the two-domain epi protein and the three domain epi of p. infestans were found to inhibit and interact with p69b, a defense subtilase of tomato, and were suggested to play a role in counterdefense. inhibitors of serine protease might be ubiquitous among eukaryotic parasites. for instance, the apicomplexan obligate parasite toxoplasma gondii secretes tgpi- and tgpi- four-domain serine protease inhibitors of the kazal family, and the intestinal hookworm ancylostoma ceylanicum secretes an 8-kda broad spectrum serine protease inhibitor of the kunitz family. here we found that kazal-like motifs also occur in saprolegniales proteins. in addition to kazal-like motifs, we also discovered a cdna that encodes a secreted protein with similarity to the cystatin class of cysteine protease inhibitors. cysteine protease inhibitors, such as chagasin, have been reported in animal parasites, mainly trypanosomids, and are thought to target proteases of the insect vector or the mammalian host. perhaps, inhibition of host proteases is a widespread counterdefense strategy in animal and plant pathogenic eukaryotes. future studies will help to address whether the discovered protease inhibitors play a role in s. parasitica-fish interactions.  this pilot cdna sequencing project provides a first look into the gene content of s. parasitica and sets the basis for genomics research in this reemerging animal pathogen. annotation of the ests revealed a number of genes that could function in virulence. future work will focus on developing molecular tools for functional analysis of s. parasitica genes. in this regards, stable transformation of saprolegnia monoica has been reported, and the rnai protocol recently developed for p. infestans should be adaptable to s. parasitica. gene expression profiling will also be applied to investigate transcriptome changes during s. parasitica-fish interactions. overall, these resources will greatly accelerate research on this important pathogen and could lead to novel perspectives for controlling saprolegniosis. 
there is a dramatic recrudescence of saprolegnia infections in aquaculture since the use of the toxic organic dye malachite green was banned in little is known about the molecular mechanisms underlying pathogenicity in s. parasitica and other animal pathogenic oomycetes. these resources will greatly accelerate research on this important pathogen. about half of the unigenes showed similarities to known protein sequences or motifs. several s. parasitica cdnas, such as those with similarity to fungal type i cellulose binding domain proteins, pan/apple module proteins, glycosyl hydrolases, proteases, as well as serine and cysteine protease inhibitors, were predicted to encode secreted proteins that could function in virulence. based on the sequence alignments of conserved proteins, the average amino acid identity between s. parasitica and three phytophthora species was 77% compared to 93% within phytophthora. some of these cdnas were more similar to fungal proteins than to other eukaryotic proteins confirming that oomycetes and fungi share some virulence components despite their evolutionary distance in this study we used a genomics approach to gain a first insight into the transcriptome of s. parasitica. we generated expressed sequence tags from a mycelial cdna library of s. parasitica. the s. parasitica sequences tended to be relatively divergent from phytophthora sequences. we provide a first glimpse into the gene content of s. parasitica, a reemerging oomycete fish pathogen. the oomycete saprolegnia parasitica is one of the most economically important fish pathogens.
 the oomycete saprolegnia parasitica is one of the most economically important fish pathogens. there is a dramatic recrudescence of saprolegnia infections in aquaculture since the use of the toxic organic dye malachite green was banned in little is known about the molecular mechanisms underlying pathogenicity in s. parasitica and other animal pathogenic oomycetes. in this study we used a genomics approach to gain a first insight into the transcriptome of s. parasitica.  we generated expressed sequence tags from a mycelial cdna library of s. parasitica. a total of consensus sequences corresponding to base pairs were assembled. about half of the unigenes showed similarities to known protein sequences or motifs. the s. parasitica sequences tended to be relatively divergent from phytophthora sequences. based on the sequence alignments of conserved proteins, the average amino acid identity between s. parasitica and three phytophthora species was 77% compared to 93% within phytophthora. several s. parasitica cdnas, such as those with similarity to fungal type i cellulose binding domain proteins, pan/apple module proteins, glycosyl hydrolases, proteases, as well as serine and cysteine protease inhibitors, were predicted to encode secreted proteins that could function in virulence. some of these cdnas were more similar to fungal proteins than to other eukaryotic proteins confirming that oomycetes and fungi share some virulence components despite their evolutionary distance  we provide a first glimpse into the gene content of s. parasitica, a reemerging oomycete fish pathogen. these resources will greatly accelerate research on this important pathogen. the data is available online through the oomycete genomics database. 