The Alveolate Perkinsus marinus: Biological Insights from EST Gene Discovery
- Sandeep J Joseph†1, 2, 5,
- José A Fernández-Robledo†3,
- Malcolm J Gardner4, 6,
- Najib M El-Sayed4, 7,
- Chih-Horng Kuo2, 8,
- Eric J Schott3, 9,
- Haiming Wang1,
- Jessica C Kissinger1, 2Email author and
- Gerardo R Vasta3Email author
© Joseph et al; licensee BioMed Central Ltd. 2010
Received: 21 August 2009
Accepted: 7 April 2010
Published: 7 April 2010
Perkinsus marinus, a protozoan parasite of the eastern oyster Crassostrea virginica, has devastated natural and farmed oyster populations along the Atlantic and Gulf coasts of the United States. It is classified as a member of the Perkinsozoa, a recently established phylum considered close to the ancestor of ciliates, dinoflagellates, and apicomplexans, and a key taxon for understanding unique adaptations (e.g. parasitism) within the Alveolata. Despite intense parasite pressure, no disease-resistant oysters have been identified and no effective therapies have been developed to date.
To gain insight into the biological basis of the parasite's virulence and pathogenesis mechanisms, and to identify genes encoding potential targets for intervention, we generated >31,000 5' expressed sequence tags (ESTs) derived from four trophozoite libraries generated from two P. marinus strains. Trimming and clustering of the sequence tags yielded 7,863 unique sequences, some of which carry a spliced leader. Similarity searches revealed that 55% of these had hits in protein sequence databases, of which 1,729 had their best hit with proteins from the chromalveolates (E-value ≤ 1e-5). Some sequences are similar to those proven to be targets for effective intervention in other protozoan parasites, and include not only proteases, antioxidant enzymes, and heat shock proteins, but also those associated with relict plastids, such as acetyl-CoA carboxylase and methyl erythrithol phosphate pathway components, and those involved in glycan assembly, protein folding/secretion, and parasite-host interactions.
Our transcriptome analysis of P. marinus, the first for any member of the Perkinsozoa, contributes new insight into its biology and taxonomic position. It provides a very informative, albeit preliminary, glimpse into the expression of genes encoding functionally relevant proteins as potential targets for chemotherapy, and evidence for the presence of a relict plastid. Further, although P. marinus sequences display significant similarity to those from both apicomplexans and dinoflagellates, the presence of trans-spliced transcripts confirms the previously established affinities with the latter. The EST analysis reported herein, together with the recently completed sequence of the P. marinus genome and the development of transfection methodology, should result in improved intervention strategies against dermo disease.
The protozoan parasite Perkinsus marinus is a facultative intracellular parasite of mollusks, and the causative agent of "dermo" disease in both wild and farmed Eastern oyster (Crassostrea virginica) populations along the Atlantic and Gulf coasts of the USA [1–4]. Since its initial description, P. marinus has generated substantial controversy with regards to its taxonomic placement [5, 6]. A close relationship with the Apicomplexa was initially proposed based on the ultrastructural analysis of the zoospore, which revealed the presence of organelles resembling an apical complex . Molecular evidence gathered over the following years revealed affinities with the Dinozoa , and even suggested this group to be its closest extant taxon [9, 10]. More recently, however, ultrastructural similarities and molecular phylogenetic affinities to Parvilucifera sp., a parasite of microeukaryotes, led to the establishment of the phylum Perkinsozoa, which like the Apicomplexa, is only comprised of endoparasites . This new phylum is considered to be one of the earliest diverging groups from the lineage leading to dinoflagellates, albeit close to the ancestor from which the ciliates, dinoflagellates, and apicomplexans originated [12–14].
Together with the various emerging diseases in the estuarine and marine environment , infections by P. marinus and other Perkinsus spp. are responsible for devastating losses in shellfisheries of economically relevant mollusk species worldwide, including oysters, clams, and abalone . Further, given the critical role that oysters and other filter-feeding bivalves play in maintaining environmental water quality, the dramatic declines in bivalve populations, caused by dermo disease, have had correspondingly detrimental impacts on the estuarine environment. Although several intervention strategies have been implemented to control dermo disease, they have had little or no success . During the early 1990s, however, the development of in vitro culture methods for P. marinus[18–20] provided a key resource that enabled studies leading to the identification of new targets for intervention . These studies ranged from fundamental cellular, biochemical, molecular and genetic studies of P. marinus biology to the direct in vitro testing of potential chemotherapeutic drugs that might suppress its proliferation [21–27].
The life cycle of P. marinus includes a free-living motile stage (zoospore) and a non-motile vegetative stage (trophozoite). The parasite is ingested by the host during filter-feeding, and is then phagocytosed by hemocytes present in the alimentary canal. Although the infection mechanism has not been fully elucidated, once inside the host, P. marinus trophozoites are recognized via a galectin present on the surface of phagocytic hemocytes [28, 29], internalized, and localize inside phagosome-like structures. The phagocytosed trophozoites remain viable by abrogating the host's respiratory burst through their effective antioxidative machinery [30, 31], and retain their proliferative capacity . Hemocyte migration throughout host tissues leads to systemic infection and eventually death [6, 33]. Therefore, in addition to those genes that mediate intrahemocytic survival, the trophozoite is likely to express additional genes that are involved in nutrient acquisition, proliferation, and pathogenesis.
Expressed sequence tag (EST) surveys  have proven to be a viable approach for gene discovery and therapeutic target identification in a variety of microbial pathogens and parasites [35–38]. EST analysis offers a rapid and valuable first glimpse of gene expression at a particular life cycle stage or under certain environmental conditions. The number of gene sequences published or available in GenBank for the genus Perkinsus is very limited, and mostly consist of ribosomal RNA (rRNA) sequences. In this study, we analyzed 31,727 ESTs generated from trophozoites from two different P. marinus strains, one of which was studied in both the presence and absence of oyster serum. Grouping of the ESTs into clusters and singletons resulted in 7,863 unique sequences. Together, they provide the first broad-based molecular view into the basic biology and cellular metabolism of this protozoan parasite of unique phylogenetic position.
All P. marinus strains were maintained in DME: Ham's F12 (1:2) supplemented with 5% fetal bovine serum (FBS) as reported elsewhere . Two cultured strains of P. marinus trophozoites, CB5D4-ATCC PRA-240 (the strain used for the sequencing the P. marinus genome; ) grown in standard medium, and TXsc-ATCC 50983  grown in both standard medium and medium supplemented with C. virginica serum (50%), were used for RNA isolation and subsequent cDNA library and EST generation.
RNA extraction and library construction
A total of four (non-normalized) libraries were constructed using P. marinus trophozoites. Three libraries were constructed from strains propagated in standard culture medium (TXsc and CB5D4), and one from the TXsc strain propagated in medium supplemented with C. virginica oyster serum (50%). Oyster serum was prepared as reported earlier . Briefly, shells were notched at the posterior end and dorsal side of the shells, cleaned with ethanol and hemolymph drawn from the adductor muscle (approximately 2 ml hemolymph per individual oyster) with a sterile syringe fitted with a 19-gauge needle. Hemolymph samples were centrifuged at 800 × g at 4°C for 10 min, the supernatant serum separated from the cell pellet, mixed at equal parts with P. marinus culture medium and used for the experimental cultures described above.
Perkinsus marinus TXsc cultures were centrifuged for 5 min at 490 × g, and RNA extracted from the pellets (1.5-2.0 ml packed P. marinus trophozoites) by the guanidine isothiocyanate method . One mg of total RNA was used for Poly A+ isolation and 5 μg of Poly A+ enriched P. marinus RNA was used to construct each P. marinus TXsc Lambda ZAP library following the manufacturer's instructions (Stratagene, La Jolla, CA). Each library was amplified once through Escherichia coli XL1-Blue MRF and stored at -80°C. Total RNA from the strain CB5D4 was extracted as above. A commercial service (Express Genomics, Frederick, MD) was used to construct the cDNA libraries. Since the libraries were not normalized and to avoid redundancy of the ESTs, a large-insert library (average insert size 2.4 kb) and a small-insert library (average insert size 1.4 kb) were constructed in pExpress 1 (Not I-Eco RV cloned into T1 phage-resistant DH10B E. coli).
P. marinus TXsc ESTs were obtained using two methods. First, 308 individual recombinant phage were selected randomly and cored out from the LB agar plates into micro-centrifuge tubes containing 400 μl of TMG (Tris-magnesium-gelatin) buffer and two drops (~50 μl) of chloroform. cDNA inserts were amplified with T3 and T7 primers, resolved in a 1.5% agarose gel (w/v), recovered using the QIAex II kit (QIAGEN, Valencia, CA), and used for direct sequencing. Second, pBluescriptKS(-) phagemid from the Lambda ZAP vectors was mass excised from both libraries and individual colonies grown for 24 h in a volume of 3 ml for plasmid preparations. In both cases, single-pass sequencing of the 5' end of the cDNA clone was carried out to generate the ESTs. P. marinus CB5D4 ESTs were sequenced using the M13 primers.
Clustering and assembly of EST/mRNA sequences
EST sequences were pre-processed to determine the sequence quality and to remove cloning vector sequences from the reads using the Phred and Cross-match software http://www.phrap.org. Poly-A/T tail trimming was done using the 'EMBOSS' Trimmest program  before submission to GenBank (accessions #EH059339 - EH090757; GR954914-GR955219; GR955352-GR955353). Sequence assembly was performed by first clustering the ESTs into groups of similar sequences, using TIGR's TGICL . Subsequently, each cluster was separately assembled into consensus sequences consisting of the longest non-redundant stretch of multiple aligned ESTs, using the CAP4 algorithm (Parcel Inc.; http://www.paracel.com). The sequences that did not cluster were treated as singletons. The cluster consensus and singleton sequences are named Pm00001-Pm07863 and are available in Additional File 1.
Annotation of the P. marinus EST sequences
Consensus sequences of the EST assemblies and singletons were compared with the NCBI non-redundant (nr) protein database (May 2009) using the BLASTX algorithm and the GenBank dbEST database (May 2009) using the TBLASTX algorithm  with the default parameters and with a cut-off E-value ≤ 1e-5. P. marinus ESTs were removed from dbEST to avoid self hits during screening. Determination of the taxonomic affinities of hits was based upon an NCBI taxonomic trace-back of best hits. For ease of presentation we have grouped the red and green algae together with the plants, and the brown algae with the Stramenopiles.
The PLAN web system (Personal BLAST Navigator, Noble Foundation) was used to assign functional annotation based on the top BLASTX hit and the gene ontology (GO) sequence database . To identify poorly conserved, or short fragments of genes contained in the ESTs, six-frame translations of the sequences were generated. This resulted in 23,888 open reading frames (ORFs) that are ≥ 75 amino acids. We searched the ORFs with Pfam (Protein Families Database) (ver. 22.0) with an E-value threshold of 0.1 to identify protein family domains. Putative signal peptides were identified using SignalP 3.0 server http://www.cbs.dtu.dk/services/SignalP/ and SecretomeP 2.0 server http://www.cbs.dtu.dk/services/SecretomeP/ using default parameters.
Analysis of orthology
We used the annotated proteins from 21 genomes [Additional file 2: Supplemental Table S1] from diverse organisms across the tree of life together with the P. marinus proteome (ORF translations of ESTs) to identify orthologs with the OrthoMCL algorithm . The OrthoMCL parameters used for the analysis were: BLASTP E-value > 1e-25 and an inflation parameter of 1.5. Multiple sequence alignment (MSA) was performed on ortholog groups that are shared by ≥ 4 taxa including P. marinus. MSA was carried out using ClustalW  enabling the 'TOSSGAPS' option and using the default values for all other parameters. The regions that contained gaps or were highly divergent were removed from the MSA by GBLOCKS  using its default settings. Phylogenetic analysis was performed on the filtered multiple sequence alignments using Seqboot, Protdist, Neighbor and Consense Tree (Phylip) , and the nearest neighbor of each P. marinus sequence was determined.
Searches for spliced leaders
The first 22 bp from the 5' end of each of the consensus ESTs were extracted and analyzed with the de novo pattern-finding algorithm MEME to identify over-represented patterns. The width for patterns was set at 22 nt and the zoops (zero or one per sequence) mode of occurrence was specified. Patterns that contained similarity or partial similarity to previously identified spliced leader (SL) sequences in dinoflagellates and Perkinsus marinus were examined in further detail. ESTs that contain potential partial SL sequences at the 5' end were enumerated, and the draft sequence of the P. marinus genome (GenBank Project ID: 12736) was searched using regular expressions to find exact SL consensus sequences. Genomic sequence (+/- 200 nt) flanking each of the putative SL loci was obtained, compared, and used to search the P. marinus ESTs.
Results and Discussion
EST sequence analysis
EST sequence numbers by library, Perkinsus strain and sequence source.
P. marinus trophozoite cultures and EST source
No. of ESTs
P. marinus CB5D4* (TIGR-JCVI)
P. marinus TXsc (without serum) (SOM/IMET)
P. marinus TXsc (oyster serum-supplemented) (SOM/IMET)
Comparison of ESTs obtained in the presence and absence of oyster serum
Previously, we have shown that supplementation of the Perkinsus culture medium with plasma from heavily infected C. virginica oysters resulted in a 32% decrease in proliferation of P. marinus; the inhibitory effect is less pronounced with plasma from moderately infected and uninfected oysters . C. virginica plasma can also specifically inhibit protease activity, including from P. marinus. Analysis of the EST sequences from the serum-supplemented library revealed: aldolase, proteases, histone-specific proteins, serine/threonin-protein kinase, deoxycytidylate deaminase and oxidoreductase enzymes, and various transporters [Additional file 2: Supplemental Table S2]. However, a comparison between the ESTs obtained from the standard and oyster serum - supplemented culture (Table 1) is not possible. There are insufficient ESTs for conducting statistically-significant comparisons under the different conditions due to the very small sample size from the serum-supplemented library. Combined sequences from both culture conditions were used for all subsequent analyses.
EST gene discovery projects for several members of the dinoflagellates including Alexandrium tamarense, Alexandrium fundyense, Karenia brevis, Lingulodinium polyedrum[60, 61], Amphidinium carterae, and Symbiodinium sp.  have been completed. TBLASTX comparisons of all P. marinus sequences to the NCBI dbEST database were performed. Figure 1 shows the top hits of the Perkinsus sequences to ESTs from various taxonomic groups. TBLASTX analysis of 7,863 sequences against dbEST resulted in 3,698 top hits (47%), within which the Alveolata were the most well-represented. Among the chromalveolates, dinoflagellates had the maximum number of top hits (690 hits; Figure 2), an observation that supports the existing phylogenetic placement of P. marinus closer to the dinoflagellates than to the apicomplexans [12, 13].
Molecular phylogenetic data have shown that P. marinus is a basal alveolate derived from ancestral dinoflagellates just after the split from apicomplexans . Further evidence indicating a relationship of P. marinus to the dinoflagellates was revealed with the discovery of SL trans-splicing in P. marinus and P. chesapeaki cDNAs , an mRNA processing phenomenon found in several dinoflagellate species . Trans-splicing in P. marinus has been reported, but the SL sequence bears a single nucleotide change with respect to the canonical sequence in dinoflagellates . A considerable number (14.7%), of the P. marinus ESTs contain partial SL sequences identical to the reported single-nucleotide non-canonical SL sequence. The longest SL identified was 13 nts and was present in 40 different sequences. Partial SL sequences of ≥ 5 nt were observed in 1,156 consensus sequences. Incomplete SL sequences in cDNA are common  and believed to be the result of inhibitors, often 5' message modifications, that interfere with reverse transcription. We have identified a putative full-length 22 nt SL in three P. marinus Nramp divalent cation transporters (Lin, Z., Fernández-Robledo, J.A., Cellier, M.F.M. Vasta, G.R., unpublished results). Numerous 22 nt SL sequences are encoded in the P. marinus genome (Joseph, S.J., and Kissinger, J.C., unpublished data), and although no full-length spliced-leader sequences were observed in the ESTs presented here, putative ESTs corresponding to portions of genomic SL loci were identified.
Functional categorization of Perkinsus sequences
Functional analysis of protein sequences
The 20 most common Pfam domains
Description of Pfam domain
No. of P. marinus ORFs
Transmembrane amino acid transporter protein
ADP-ribosylation factor family
Helicase conserved C-terminal domain
Mitochondrial carrier protein
Papain family cysteine peptidase
Protein kinase domain
Proteasome A-type and B-type
RNA recognition motif 1
RNA recognition motif 2
Sugar (and other) transporter
Tetra tricopeptide repeat
WD domain, G-beta repeat
Zinc finger C3HC4 type
Zinc finger C8C5C3H type
ATPase family associated with various cellular activities
DEAD/DEAH box helicase
Dna J domain: associated with hsp 70 heat-shock protein
Cytochrome b5-like Heme/Steroid binding domain
Pfam domains of interest identified in P. marinus EST cluster open reading frames.
Pfam family domain
Description of the domain
P. marinus Cluster IDs
Papain family cysteine protease
Pm02598, Pm05550, Pm02598, Pm01528, Pm00748, Pm00174, Pm00175, Pm00188, Pm01093, Pm02540, Pm00187
Pm03998, Pm00516, Pm01387
UIp1 protease family, C-terminal catalytic domain
Peptidase family C54
Peptidase family C69
Pm00327, Pm01089, Pm00100, Pm00136
Peptidase family M1
Zinc carboxy peptidase
Pm05750, Pm04910, Pm01740, Pm00611, Pm01740, Pm07103
Cytosol amino peptidase
Aminopeptidase I zinc metalloprotease
Peptidase family M20
Metallopeptidase family M24
Pm00721, Pm05108, Pm01923, Pm01484, Pm03328, Pm05608, Pm01369
Peptidase family M3
Peptidase family M41
Peptidase family M48
Serine carboxy peptidase
Pm04952, Pm03877, Pm07791, Pm01720, Pm07092, Pm05386, Pm02482, Pm07066, Pm02294, Pm00207, Pm00346, Pm06037
X-Pro dipeptidyl peptidase
Peptidase S24 like
Pm05197, Pm00155, Pm01368, Pm00344, Pm00265, Pm02193, Pm03660, Pm00235, Pm00154
Prolyl oligo peptidase family
Iron/manganese superoxide dismutase, alpha-hairpin
Pm00420, Pm05057, Pm02278
Iron/manganese superoxide dismutase, C-terminal domain
Pm00420, Pm05057, Pm02278
Pm07679, Pm07604, Pm04790, Pm04765
Pm01390, Pm02808, Pm00345, Pm02770, Pm07686, Pm03543, Pm01102, Pm00367, Pm03381, Pm01163, Pm03631, Pm07411, Pm02203, Pm03850, Pm05817
Fatty Acid Synthesis
Fatty acid desaturase
Pm02182, Pm03006, Pm06024, Pm07340, Pm02313, Pm02648, Pm03710, Pm03609
MEP cytidylyl transferase
Glycosyl transferase Family
Alg9-like mannosyl transferase
Lipopolysaccharide galactosyl tranferase
Oligosaccharide biosynthesis protein Alg14-like
Heat Shock Proteins
Heat shock protein 90
Pm00142, Pm00170, Pm00836, Pm03154
Heat shock protein 70
Pm00172, Pm05287, Pm00172, Pm03153, Pm05233, Pm06933, Pm00065, Pm00137, Pm01460, Pm00363, Pm01884
HSp60 chaperonin family and the TCP-1 family
Pm03177, Pm02627, Pm01137, Pm02465, Pm00825, Pm03422, Pm04033, Pm02165, Pm01657, Pm05905
Ricin-type beta-trefoil domain
Pm01439, Pm02264, Pm03133
Legume -like lectin family
Phosphoribosyl transferase adenylosuccinate synthetase
Pm04986, Pm07718, Pm03005
To determine the number of potential orthologs between P. marinus and proteins from 21 organisms representative of the major divisions of the tree of life, [Additional file 2: Supplemental Table S1] we performed a clustering analysis on the non-redundant 21 genome protein sets along with the P. marinus EST ORF proteome using the OrthoMCL program . In brief, OrthoMCL defines clusters based on reciprocal best BLASTP hits and sorts proteins without a reciprocal best BLASTP hit into the best matching cluster. The input dataset for OrthoMCL consisted of 325,188 protein sequences from 22 proteomes. P. marinus sequences (5,661 ORFs) were found in 2,878 ortholog groups, of which 1,715 are unique to P. marinus i.e., no sequences from other taxa are present in this group. The extent to which P. marinus shares orthologous genes with the 21 other taxa examined is listed in Additional file 2: Supplemental Table S4. As the number of taxa increases, the number of shared genes decreases.
Automated phylogenetic analyses of the orthologous clusters was performed to determine the nearest taxonomic neighbor for each sequence. This approach is useful for establishing gene origins, identifying genes with restricted phylogenetic distribution and to identify putative horizontal gene transfer (HGT) events in the P. marinus genome. Phylogenetic trees were constructed from the 1,042, ortholog groups that shared ≥ 4 taxa including P. marinus [Additional file 2: Supplemental Table S4]. Multiple alignments were created with ClustalW and filtered using GBLOCKS to remove ambiguous regions. Only ortholog groups with ≥ 50 aligned aa were used for phylogenetic analysis, and 291 alignments met this criterion. Neighbor-joining trees were constructed with bootstrap support to determine the P. marinus nearest neighbor. Not surprisingly, 54% of the orthologous groups have alveolates as the closest neighbor to P. marinus (48% apicomplexans and 6% ciliates (dinoflagellates are not represented)). Toxoplasma gondii was the most highly represented species. P. marinus sequences also show nearest neighbors of kinetoplastids, bacteria, archaea and red algae. For the gene encoding 2-C-methyl-D-erythritol 4-phosphate cytidyly transferase (IspD), one of the seven enzymes involved in plastid metabolism recently discovered in P. marinus[75–77] its nearest neighbor is that from the red algae Cyanidioschyzon merolae. Other important taxonomic groups that were found to be closest neighbors to P. marinus include the Heterokontophyta [Thalassiosira pseudonana (17 genes), and Phytopthora ramorum (36 genes)], Plantae [Arabidopsis thaliana (26 genes)], Animalia [Homo sapiens (24 genes), Drosophila melanogaster (7 genes)], and Fungi [Saccharomyces cerevisiae (11 genes)]. It should be noted that the nearest neighbors are expected to change as more taxa are sequenced and added to the analysis.
Sequences of particular interest
P. marinus consensus sequences with similarity to known proteases.
P. marinus EST cluster ID
GenBank best hit accession
Species with best hit
Subtilisin-like serine protease
Cathepsin-like cysteine protease
Papain cysteine protease
Rhomboid-like protease 1
Calcium-dependent cysteine protease, putative
ATP-dependent protease La
Serine protease, putative
Metacaspase 1 precursor, putative
Putative LON protease
Peptidase family M48 domain
P. marinus sequences with similarity to known oxidative enzymes.
P. marinus EST cluster ID
GenBank best hit accession
Species with best blast hit
Superoxide dismutase 1
Superoxide dismutase 2
Pediculus humanus corporis
Fatty acid synthesis
Lipid analysis of 7 day-old in vitro cultured P. marinus trophozoites indicated that triacylglycerol represents 48.7% of the total lipids . P. marinus trophozoites utilize 13C-acetate to synthesize a range of saturated and unsaturated fatty acids and the parasite's ability to synthesize 20:4(n-6) de novo is unique within parasitic protozoa . Eukaryotes employ either the delta-6 or delta-8 desaturase pathway, or both, to synthesize arachidonic acid, an essential fatty acid. The meront stage of P. marinus synthesizes arachidonic acid through the delta-8 pathway . In addition, it has been suggested that P. marinus cannot synthesize sterols and must sequester them from its host. Perkinsus cells are able to proliferate in complete lipid supplement medium (cod liver oil, cholesterol and alpha tocopherol acetate in detergent) and media containing cholesterol or cholesterol+alpha tocopherol acetate, but fail to proliferate in control medium and medium containing just alpha tocopherol acetate . However, the genome of P. marinus encodes 6 out of 7 methylerythrithol phosphate (MEP) pathway genes [75–77] indicating that Perkinsus is able to synthesize de novo sterols (see below: isoprenoid metabolism). ESTs matching enzymes involved in sterol metabolism also include sterol glucosyltransferases (Pm04113, E = 1e-47; Pm01717, E = 4e-42; Pm04859, E = 6e-44), sterol C-24 reductase (Pm01729, E = 3e-50), and sterol desaturase (Pm00771, E = 9e-23). Three P. marinus genes encoding the enzymes responsible for arachidonic acid biosynthesis (C18 delta-9-elongating activity, C20 delta-8 desaturase, C20 delta-5 desaturase) are clustered and co-transcribed as an operon . Sequences highly similar to delta-9 desaturase (Pm03609, E-value = 1e-66), delta-8 fatty acid desaturase (Pm07340, E-value = 1e-161), and delta-5 desaturase (Pm03710, E-value = 0.00) were present, while no sequences similar to delta-6 fatty acid desaturase were identified, in agreement with the abovementioned observations. Pfam analysis identified 8 fatty acid desaturase domains including the above sequences (Table 3).
Sequences similar to the plastid-localized enzyme, acetyl-CoA carboxylase (Pm00609 (E-value = 2e-41) and Pm05907 (E-value = 3e-16)) involved in fatty acid biosynthesis were also identified in this study. Indeed, it has been shown that Perkinsus proliferation is inhibited by Triclosan and cerulenin, which has been interpreted as evidence for the presence of a plastidic FAS II pathway [24, 25]. However, when considering the effect of Triclosan as indicative of the relevance of the apicoplast FASII biosynthesis, results should be interpreted with caution, since in Plasmodium, the antimalarial activity of Triclosan is not targeted to FabI . Further, although Theileria is also susceptible to this drug, genes coding for FASII are lacking in the Theileria genome .
Heat shock proteins
Expression of heat shock protein 70 (HSP70) in P. marinus suggests that the parasite might use general stress response genes to overcome the stress imposed by the host environment. In Toxoplasma, the HSP70 gene is expressed during the transition from the active to latent form . In virulent Toxoplasma strains, HSP70 contains seven amino acids not present in the HSP70 from non-virulent strains, and HSP70 expression is elevated 2-fold in virulent versus non-virulent strains . A highly similar consensus sequence to the T. gondii HSP70 protein (Pm00065, Evalue = 0.00) was identified in the P. marinus EST analysis. Further, two sequences similar to HSP90 were identified, Pm00171 (E-value = 1e-176) and Pm05367 (E-value = 3e-60). Many heat shock proteins are chaperonins, which function in concert with proteins similar to the delta subunit of the t-complex family of chaperonins, and are found in virtually all organisms including Leishmania. In different species, t-complex chaperonins are involved in protein folding after stress-related denaturation. Two chaperonin homologs, encoded by Pm03177 and Pm01657, are similar to proteins with roles in thermo-tolerance, cell-cycle progression and hematopoeisis. Pfam analysis detected the same HSP90 and HSP70 proteins as well as and Cpn_TCP1 (HSP60 chaperonin family-TCP-1 family) domains (Table 3).
P. marinus consensus sequences with similarity to enzymes involved in plastid isoprenoid metabolism and other plant-type plastid proteins.
P. marinus EST Cluster ID
GenBank best hit accession
Species with best hit
2C-methyl-D-erythritol 2,4-cyclodiphosphate synthase
4-hydroxy-3-methyl-2-butenyl 1-diphosphate synthase
Triose phosphate/phosphate translocator, non-green plastid, chloroplast precursor, putative
Cytidylyltransferase family protein
2-C-methyl-D-erythritol 4-phosphate cytidylyltransferase
Candidatus Protochlamydia amoebophila
Glycan assembly and carbohydrate-binding proteins
Oligosaccharides on cell surface and intracellular glycoconjugates are assembled by the combined activity of glycosyltransferases and glycosidases, and may interact with carbohydrate-binding proteins (lectins). Further, protein-carbohydrate recognition plays a key role in intracellular processes such as protein folding and transport, as well as interactions between cells or cells and the extracellular matrix in functions related to development, immune responses, and host-parasite interactions. It is widely accepted that recognition of parasitic surface glycoconjugates by host humoral or cell-associated lectins are frequently involved in cellular recognition and colonization of the host. Several studies have partially characterized multiple lectins in the oyster plasma and hemocytes, which recognize glycoproteins that display terminal galactose and N-acetylated sugars [96, 97]. Flow cytometry analysis of P. marinus trophozoites labeled with commercial lectins  enabled a tentative identification of some of the sugars present on the surface of the P. marinus trophozoites. Recently, we demonstrated that an oyster galectin present on the hemocyte surface functions as a parasite receptor and facilitates host entry . Once inside the host, trophozoites survive intra-hemocytic killing mechanisms and proliferate . The P. marinus ESTs include several with high similarity to glycosyltransferases and other members of the 'sugar moiety' transferase enzyme family. Four sequences (Pm01326 (E-value = 2e-49), Pm02039 (E-value = 9e-31), Pm07403 (E-value = 1e-08), Pm02460, E-value = 3e-09) show high similarity to the genes that encode glycosyltransferase in Cryptosporidium hominis and Verrucomicrobium spinosum. BlastX hits to N-acetylglucosaminyltransferase (Pm07672, E-value = 4e-18), mannosyltransferase (Pm03101, E-value = 6e-36), sterol-glucosyltransferase (Pm04859, E-value = 6e-44), prolipoprotein diacylglyceryl transferase (Pm04168, E-value = 2e-23) and alpha glucosidase (Pm03682 (E-value = 2e-44), Pm01066 (E-value = 2e-50) were also observed. In addition, Pfam domains for the glycosyltransferase family, Glyco_transf_20 (trehalose-6-phosphate transferase), Glyco_transf_22 (Alg9-like mannosyl transferase), Glyco_transf_28 (UDP-N-acetyl glucosaminine transferase and Glyco_transf_8 (lipopolysaccaride galactosyl transferase) were observed (Table 3). A single EST ORF encoding the domain called Alg14 (oligosaccharide biosynthesis protein Alg14 like), which represents an important protein in the synthesis of glycoconjugates, was also identified (Table 3). Prior studies have revealed glycosidase activity, including β-D-glucosidase, β-D-xylosidase, N-acetyl β-D-glucosaminidase and N-acetyl β-D-galactosaminidase in Perkinsus trophozoites, and N-acetyl β-D-glucosaminidase in spent culture medium (Ahmed, H., Fernández-Robledo, J.A., Vasta, G.R., unpublished results).
Three unique EST clusters (Pm03024 (E-value = 1e-36), Pm03327 (E-value = 1e-30) and Pm02250 (E-value = 7e-26)) show similarity to the ERGIC-53-like mannose binding lectin also present in Cryptosporidium, Toxoplasma and Plasmodium species. ERGIC-53 is a type 1 transmembrane L-type lectin present in the endoplasmic reticulum (ER) that captures correctly folded glycoproteins, and mediates their transport along the secretory pathway. The yeast L-lectins Emp46p and EMP47p, homologues of ERGIC-53, have been proposed to be transport receptors that facilitate recruitment of glycoproteins into vesicles budding from the ER . One unique sequence (Pm04230, E-value = 7e-12) shows similarity to the mannose-binding lectin derived from Crinum asiaticum, a plant lectin with homologies to known monocot mannose-binding lectins from Amaryllidaceae, Orchidaceae, Alliaceae and Liliaceae, high similarity to the gastrodianin-type antifungal proteins and a predicted structure similar to that of the Galanthus nivalis agglutinin . The sequence Pm07056 shows similarity to a galactose-specific lectin. Pfam analysis also detected similarities with plant lectins, such as domains for the ricin-type beta-trefoil lectin and legume-like-lectin (Table 3).
Several nucleotide salvage enzymes were identified. Hypoxanthine-guanine phosphoribosyl transferase, HGPRT (Pm04986, E-value = 2e-57), and hypoxanthine-xanthine-guanine phosphoribosyl transferase, HXGPRT (Pm07697, E-value = 6e-53); uridine kinase - uracil phosphoribosyltransferase, UK-UPRT (Pm01230, E-value = 1e-40) and uracil phosphoribosyltransferase, UPRT (Pm07824, E-value = 2e-43) were detected. Pfam analysis revealed the presence of phosphoribosyl transferase in 3 ESTs including Pm04986 (Table 3). Evidence for a few de novo biosynthetic enzymes was detected. There is a hit to dihydroorotate oxidase (Pm05343, E-value = 7e-30), uridylate kinase (Pm05827, E-value = 3e-19) and adenylosuccinate synthetase (Pm02370, E-value = 1e-128) all involved in pyrimidine biosynthesis. Pfam analysis also revealed an adenylosuccinate synthetase domain in Pm02370.
Potential targets for chemotherapy
Several potential candidates for chemotherapy in P. marinus were identified in this EST study. Like the proteases from other parasites, known to degrade host proteins for acquisition of the nutrients, proteases produced by P. marinus could be targeted for chemotherapy. Further, the plasma of the eastern oyster C. virginica contains a serine protease inhibitor that binds tightly to subtilisin and perkinsin, and is potentially involved in blocking the parasite's proteolytic activity . Therefore, dermo disease-resistant oyster populations could be established by introducing through selective breeding selected gene variants that exhibit enhanced expression of protease inhibitors. Further, based on the observed in vitro susceptibility of the parasite to selected ROIs  and hemocyte-based defense against the parasite in vivo, therapeutic agents could be applied to infected oysters to enhance their respiratory burst in response to P. marinus. Similarly, and based on the observation that parasites are usually more susceptible to ROIs than are their hosts, components of the P. marinus anti-oxidative stress pathway could be the basis for development of therapeutic drugs. Conventional molecular modeling approaches such as the Drug Design by Receptor Fit (DDRF) methods can be applied to the parasite enzymes such as glutathione reductase, peroxiredoxin, thioredoxin, and glutaredoxins, whose sequences were identified in this study. The P. marinus EST analysis also identified sequences similar to aldolase, an enzyme that has been used as a target for therapy in several Plasmodium species because of its central role in energy metabolism (glycolysis) . Aldolases have been also used for intervention in Trypanosoma brucei infection . Therefore, the potential for aldolase inhibitors, already in use for other parasites, to inhibit P. marinus growth in vitro warrants further investigation. The recent characterization of the isoprenoid pathway in P. marinus strongly suggests the presence of a cryptic plastid in the parasite, which has been identified as an excellent target for drug therapy in apicomplexan parasites [104–107]. Analytical methods  as well as EST evidence from the current study point towards the presence of a DOXP/ME pathway (Table 5) to produce isopentenyl diphosphate. The enzymatic machinery of the DOXP/MEP pathway in Plasmodium falciparum has been fully characterized  and fosmidomycin, a specific inhibitor of the DOXP reductoisomerase, is very effective against malaria . Interestingly, fosmidomycin has also been tested against Perkinsus trophozoites and it appears that up to 5 days concentrations up to 1 mM show no effect on growth inhibition . The genes for the DOXP pathway are present in Toxoplasma and it has been proposed  that toxoplasmosis could be treated by targeting a downstream pathway enzyme, farnesyl diphosphate synthase (FPPS), using bisphosphonates, which are specific FPPS inhibitors. A similar strategy could also be applied to dermo disease if bisphosphonates do not affect the oyster or if fosmidomycin does not show differential effectiveness against the parasite DOXP pathway. Therapeutic strategies using existing drugs such as bisphosphonates and fosmidomycin have the advantage of avoiding costs of de novo drug design and development. Moreover, virtual screening initiatives could provide new avenues for drug development against numerous protozoan parasites . Conversely, these observations also highlight the potential for using P. marinus as a readily-cultured, non-pathogenic model for early screening of potential drugs against a variety of protistan human parasites.
P. marinus sequences display the greatest similarity to EST sequences from dinoflagellates. No significant differences were observed between EST populations obtained from parasites propagated under standard conditions and those exposed to oyster serum. This finding is a consequence of the small number of ESTs sampled, and a more rigorous analysis should be carried out by expanding both sample size and incorporating additional experimental strategies. Concerning the latter, to gain a better understanding P. marinus virulence, approaches such as subtractive techniques [112, 113] and microarray analysis will be very useful. Although our results are a first step in that direction, the application of subtractive techniques should result in ESTs overlooked in our study due to the limitations of non-normalized libraries. The libraries and ESTs generated here, however, may find further use in the production of microarrays to visualize changes in gene expression, such as expression of parasite genes related to defense against the oyster's immune system.
As indicated above, P. marinus expresses sequences with significant similarity of dinoflagellates followed by the apicomplexans. Although the fraction of P. marinus transcripts that are trans-spliced is still unknown, the identification of a SL in >14% of the ESTs confirms the previously established affinity with dinoflagellates. It also suggests that PCR amplification based on the Perkinsus variant of the SL would provide a rapid and efficient method of amplifying and cloning full-length transcripts in the future.
The EST analysis reported herein, together with the recently completed P. marinus genome sequencing project (GenBank Project ID: 12736), and the development of a transfection system for Perkinsus trophozoites  will enhance the community's ability to improve the status of both natural and farmed oyster stocks by identifying gene products suitable for drug targeting, which will lead to therapeutic applications that may be effective in closed (especially hatchery) systems. The genes relevant to host-parasite interactions, particularly those involved in host-cell entry and/or pathway signaling may lead to genetically-selected or -engineered oysters that either block the entry of the parasite or enhance the response of the oyster defense against P. marinus. Production of seed oysters that remain disease-free and reach marketable size will be critical for full recovery of wild eastern oyster populations, which provide irreplaceable environmental services. Disease-resistant oysters would also form the basis of a viable shellfishery, as well as sustainable production of farmed oysters.
Funding for this work was provided by NOAA - MD Sea Grant SA7528068-I to G.R.V and J.C.K, and in part by the NSF/USDA-CSREES Microbial Sequencing Program (Award NSF/USDA 0333240) to M.J.G and N.M.E.S (JCVI), and G.R.V and J.A.F.R. (SOM/IMET) and by NSF Award IOS-0618409 to G.R.V. We thank Adam G. Marsh for making the initial P. marinus Lambda ZAP libraries, Dr. Jeremy DeBarry for useful discussions and Ganesh Srinivasamoorthy for help with PERL scripts and data parsing.
- Ford SE: Range extension by the oyster parasite Perkinsus marinus into the northeastern United States: Response to climate change?. J Shellfish Res. 1996, 15 (1): 45-56.Google Scholar
- Ford SE, Chintala MM: Northward expansion of a marine parasite: Testing the role of temperature adaptation. J Exp Mar Biol Ecol. 2006, 339 (2): 226-235. 10.1016/j.jembe.2006.08.004.Google Scholar
- Ford SE, Smolowitz R: Infection dynamics of an oyster parasite in its newly expanded range. Mar Biol. 2007, 151 (1): 119-133. 10.1007/s00227-006-0454-6.Google Scholar
- Soniat TM: Epizootiology of Perkinsus marinus disease of eastern oysters in the Gulf of México. J Shellfish Res. 1996, 15 (1): 35-43.Google Scholar
- Reece KS, Siddall ME, Burreson EM, Graves JE: Phylogenetic analysis of Perkinsus based on actin gene sequences. J Parasitol. 1997, 83 (3): 417-423. 10.2307/3284403.PubMedGoogle Scholar
- Perkins FO: The structure of Perkinsus marinus (Mackin, Owen and Collier, 1950) Levine, 1978 with comments on taxonomy and phylogeny of Perkinsus spp. J Shellfish Res. 1996, 15 (1): 67-87.Google Scholar
- Levine ND: Perkinsus genus and other new taxa in the protozoan phylum Apicomplexa. J Parasitol. 1978, 64 (3): 10.2307/3279807.Google Scholar
- Goggin CL, Barker SC: Phylogenetic position of the genus Perkinsus (Protista, Apicomplexa) based on small subunit ribosomal RNA. Mol Biochem Parasitol. 1993, 60 (1): 65-70. 10.1016/0166-6851(93)90029-W.PubMedGoogle Scholar
- Siddall ME, Reece KS, Graves JE, Burreson EM: 'Total evidence' refutes the inclusion of Perkinsus species in the phylum Apicomplexa. Parasitology. 1997, 115 (Pt 2): 165-176. 10.1017/S0031182097001157.PubMedGoogle Scholar
- Kuvardina ON, Leander BS, Aleshin VV, Myl'nikov AP, Keeling PJ, Simdyanov TG: The phylogeny of colpodellids (Alveolata) using small subunit rRNA gene sequences suggests they are the free-living sister group to apicomplexans. J Eukaryot Microbiol. 2002, 49 (6): 498-504. 10.1111/j.1550-7408.2002.tb00235.x.PubMedGoogle Scholar
- Norén F, Moestrup Ø, Rehnstam-Holm AS: Parvilucifera infectans Norén et Moestrup gen. et sp. nov. (Perkinsozoa phylum nov.): A parasitic flagellate capable of killing toxic microalgae. Eur J Protist. 1999, 35 (3): 233-254.Google Scholar
- Saldarriaga JF, McEan ML, Fast NM, Taylor FJR, Keeling PJ: Multiple protein phylogenies show that Oxyrrhis marina and Perkinsus marinus are early branches of the dinoflagellate lineage. Int J Syst Evol Microbiol. 2003, 53: 355-365. 10.1099/ijs.0.02328-0.PubMedGoogle Scholar
- Gile GH, Patron NJ, Keeling PJ: EFL GTPase in cryptomonads and the distribution of EFL and EF-1alpha in chromalveolates. Protist. 2006, 157 (4): 435-444. 10.1016/j.protis.2006.06.002.PubMedGoogle Scholar
- Hoppenrath M, Leander BS: Molecular phylogeny of Parvilucifera prorocentri (Alveolata, Myzozoa): Insights into perkinsid character evolution. J Eukaryot Microbiol. 2009, 56 (3): 251-256. 10.1111/j.1550-7408.2009.00395.x.PubMedGoogle Scholar
- Harvell CD, Kim K, Burkholder JM, Colwell RR, Epstein PR, Grimes DJ, Hofmann EE, Lipp EK, Osterhaus AD, Overstreet RM: Emerging marine diseases-climate links and anthropogenic factors. Science. 1999, 285 (5433): 1505-1510. 10.1126/science.285.5433.1505.PubMedGoogle Scholar
- Villalba A, Reece KS, Ordas A, Casas SM, Figueras A: Perkinsosis in molluscs: A review. Aquat Living Res. 2004, 17: 411-432. 10.1051/alr:2004050.Google Scholar
- Krantz GE, Jordan SJ: Management alternatives for protecting Crassostrea virginica fisheries in Perkinsus marinus enzootic and epizootic areas. J Shellfish Res. 1996, 15 (1): 167-176.Google Scholar
- Gauthier JD, Vasta GR: Continuous in vitro culture of the eastern oyster parasite Perkinsus marinus. J Invert Pathol. 1993, 62 (3): 321-323. 10.1006/jipa.1993.1122.Google Scholar
- Kleinschuster SJ, Swink SL: A simple method for the in vitro culture of Perkinsus marinus. Nautilus. 1993, 107 (2): 76-78.Google Scholar
- La Peyre JF, Faisal M, Burreson EM: In vitro propagation of the protozoan Perkinsus marinus, a pathogen of the eastern oyster, Crassostrea virginica. J Eukaryot Microbiol. 1993, 40 (3): 304-310. 10.1111/j.1550-7408.1993.tb04920.x.Google Scholar
- Faisal M, Schafhauser DY, Garreis KA, Elsayed E, La Peyre JF: Isolation and characterization of Perkinsus marinus proteases using bacitracin-sepharose affinity chromatography. Comp Biochem Physiol, B. 1999, 123B (4): 417-426. 10.1016/S0305-0491(99)00088-7.Google Scholar
- Delaney MA, Brady YJ, Worley SD, Huels KL: The effectiveness of N-halamine disinfectant compounds on Perkinsus marinus, a parasite of the eastern oyster Crassostrea virginica. J Shellfish Res. 2003, 22 (1): 91-94.Google Scholar
- Elandalloussi LM, Leite RB, Rodrigues PM, Afonso R, Nunes PA, Cancela ML: Effect of antiprotozoal drugs on the proliferation of the bivalve parasite Perkinsus olseni. Aquaculture. 2005, 243 (1-4): 9-17. 10.1016/j.aquaculture.2004.09.038.Google Scholar
- Lund ED, Soudant P, Chu F-LE, Harvey E, Bolton S, Flowers A: Effects of Triclosan on growth, viability and fatty acid synthesis of the oyster protozoan parasite Perkinsus marinus. Dis Aquat Org. 2005, 67 (3): 217-224. 10.3354/dao067217.PubMedGoogle Scholar
- Stelter K, El-Sayed NM, Seeber F: The expression of a plant-type ferredoxin redox system provides molecular evidence for a plastid in the early dinoflagellate Perkinsus marinus. Protist. 2007, 158 (1): 119-130. 10.1016/j.protis.2006.09.003.PubMedGoogle Scholar
- Teles-Grilo ML, Tato-Costa J, Duarte SM, Maia A, Casal G, Azevedo C: Is there a plastid in Perkinsus atlanticus (Phylum Perkinsozoa)?. Eur J Protistol. 2007, 43 (2): 163-167. 10.1016/j.ejop.2007.02.002.PubMedGoogle Scholar
- Venegas-Caleron M, Beaudoin F, Sayanova O, Napier JA: Co-transcribed genes for long chain polyunsaturated fatty acid biosynthesis in the protozoon Perkinsus marinus include a plant-like FAE1 3-ketoacyl coenzyme A synthase. J Biol Chem. 2007, 282 (5): 2996-3003. 10.1074/jbc.M607051200.PubMedGoogle Scholar
- Tasumi S, Vasta GR: A galectin of unique domain organization from hemocytes of the eastern oyster (Crassostrea virginica) Is a receptor for the protistan parasite Perkinsus marinus. J Immunol. 2007, 179 (5): 3086-3098.PubMedGoogle Scholar
- Vasta GR: Roles of galectins in infection. Nat Rev Microbiol. 2009, 7: 424-438. 10.1038/nrmicro2146.PubMed CentralPubMedGoogle Scholar
- Volety AK, Chu FL: Suppression of chemiluminescence of eastern oyster (Crassostrea virginica) hemocytes by the protozoan parasite Perkinsus marinus. Dev Comp Immunol. 1995, 19 (2): 135-142. 10.1016/0145-305X(94)00059-O.PubMedGoogle Scholar
- Anderson RS: Perkinsus marinus secretory products modulate superoxide anion production by oyster (Crassostrea virginica) haemocytes. Fish Shellfish Immunol. 1999, 9 (1): 51-60. 10.1006/fsim.1998.0174.Google Scholar
- Alavi MR, Fernández-Robledo JA, Vasta GR: In vitro intracellular survival of Perkinsus marinus trophozoites upon phagocytosis by oyster (Crassostrea virginica and Crassostrea ariakensis) hemocytes. J Parasitol. 2009, 1-Google Scholar
- Chu F-LE: Laboratory investigations of susceptibility, infectivity, and transmission of Perkinsus marinus in oysters. J Shellfish Res. 1996, 15 (1): 57-66.Google Scholar
- Adams MD, Kerlavage AR, Fleischmann RD, Fuldner RA, Bult CJ, Lee NH, Kirkness EF, Weinstock KG, Gocayne JD, White O: Initial assessment of human gene diversity and expression patterns based upon 83 million nucleotides of cDNA sequence. Nature. 1995, 377 (6547 Suppl): 3-174.PubMedGoogle Scholar
- Ajioka JW, Boothroyd JC, Brunk BP, Hehl A, Hillier L, Manger ID, Marra M, Overton GC, Roos DS, Wan KL: Gene discovery by EST sequencing in Toxoplasma gondii reveals sequences restricted to the Apicomplexa. Genome Res. 1998, 8 (1): 18-28.PubMedGoogle Scholar
- Prichard R, Tait A: The role of molecular biology in veterinary parasitology. Vet Parasitol. 2001, 98 (1-3): 169-194. 10.1016/S0304-4017(01)00429-0.PubMedGoogle Scholar
- Tarleton RL, Kissinger J: Parasite genomics: current status and future prospects. Curr Opin Immunol. 2001, 13 (4): 395-402. 10.1016/S0952-7915(00)00233-8.PubMedGoogle Scholar
- Pappas GJ, Benabdellah K, Zingales B, Gonzalez A: Expressed sequence tags from the plant trypanosomatid Phytomonas serpens. Mol Biochem Parasitol. 2005, 142 (2): 149-157. 10.1016/j.molbiopara.2005.03.017.PubMedGoogle Scholar
- Gauthier JD, Vasta GR: In vitro culture of the eastern oyster parasite Perkinsus marinus: Optimization of the methodology. J Invert Pathol. 1995, 66 (2): 156-168. 10.1006/jipa.1995.1079.Google Scholar
- Fernández-Robledo JA, Lin Z, Vasta GR: Transfection of the protozoan parasite Perkinsus marinus. Mol Biochem Parasitol. 2008, 157 (1): 44-53. 10.1016/j.molbiopara.2007.09.007.PubMedGoogle Scholar
- Gauthier JD, Vasta GR: Effects of plasma from bivalve mollusk species on the in vitro proliferation of the protistan parasite Perkinsus marinus. J Exp Zool. 2002, 292 (3): 221-230. 10.1002/jez.10013.PubMedGoogle Scholar
- Sambrook J, Russell DW: Molecular Cloning A Laboratory Manual. 2001, New York: CHSL Press, ThirdGoogle Scholar
- Rice P, Longden I, Bleasby A: EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000, 16 (6): 276-277. 10.1016/S0168-9525(00)02024-2.PubMedGoogle Scholar
- Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B: TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003, 19 (5): 651-652. 10.1093/bioinformatics/btg034.PubMedGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.PubMedGoogle Scholar
- He J, Dai X, Zhao X: PLAN: a web platform for automating high-throughput BLAST searches and for managing and mining results. BMC Bioinformatics. 2007, 8: 53-10.1186/1471-2105-8-53.PubMed CentralPubMedGoogle Scholar
- Li L, Stoeckert CJ, Roos DS: OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003, 13 (9): 2178-2189. 10.1101/gr.1224503.PubMed CentralPubMedGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: Clustal-W - Improving the Sensitivity of Progressive Multiple Sequence Alignment through Sequence Weighting, Position-Specific Gap Penalties and Weight Matrix Choice. Nucl Acids Res. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralPubMedGoogle Scholar
- Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000, 17 (4): 540-552.PubMedGoogle Scholar
- Felsenstein J: Phylogeny Inference Package (PHYLIP). Version 3.5c edn. 1993, Seattle: University of Washington, Seattle, WashingtonGoogle Scholar
- Li L, Brunk BP, Kissinger JC, Pape D, Tang K, Cole RH, Martin J, Wylie T, Dante M, Fogarty SJ: Gene discovery in the apicomplexa as revealed by EST sequencing and assembly of a comparative gene database. Genome Res. 2003, 13 (3): 443-454. 10.1101/gr.693203.PubMed CentralPubMedGoogle Scholar
- Leggat W, Hoegh-Guldberg O, Dove S, Yellowlees D: Analysis of an EST library from the dinoflagellate (Symbiodinium sp.) symbiont of reef-building corals. J Phycol. 2007, 43 (5): 1010-1021. 10.1111/j.1529-8817.2007.00387.x.Google Scholar
- Coss CA, Robledo JA, Ruíz GM, Vasta GR: Description of Perkinsus andrewsi n. sp. isolated from the Baltic clam (Macoma balthica) by characterization of the ribosomal RNA locus, and development of a species-specific PCR-based diagnostic assay. J Eukaryot Microbiol. 2001, 48 (1): 52-61. 10.1111/j.1550-7408.2001.tb00415.x.PubMedGoogle Scholar
- Coss CA, Robledo JA, Vasta GR: Fine structure of clonally propagated in vitro life stages of a Perkinsus sp. isolated from the Baltic clam Macoma balthica. J Eukaryot Microbiol. 2001, 48 (1): 38-51. 10.1111/j.1550-7408.2001.tb00414.x.PubMedGoogle Scholar
- Ellis JT, Morrison DA, Jeffries AC: The phylum Apicomplexa: an update on the molecular phylogeny. 1998, Boston: Kluwer Academis PublishersGoogle Scholar
- Leander BS, Clopton RE, Keeling PJ: Phylogeny of gregarines (Apicomplexa) as inferred from small-subunit rDNA and beta-tubulin. Int J Syst Evol Microbiol. 2003, 53 (Pt 1): 345-354. 10.1099/ijs.0.02284-0.PubMedGoogle Scholar
- Hackett JD, Scheetz TE, Yoon HS, Soares MB, Bonaldo MF, Casavant TL, Bhattacharya D: Insights into a dinoflagellate genome through expressed sequence tag analysis. BMC Genomics. 2005, 6 (1): 80-10.1186/1471-2164-6-80.PubMed CentralPubMedGoogle Scholar
- Erdner DL, Anderson DM: Global transcriptional profiling of the toxic dinoflagellate Alexandrium fundyense using Massively Parallel Signature Sequencing. BMC Genomics. 2006, 7: 88-10.1186/1471-2164-7-88.PubMed CentralPubMedGoogle Scholar
- Lidie KB, Ryan JC, Barbier M, Van Dolah FM: Gene expression in Florida red tide dinoflagellate Karenia brevis: analysis of an expressed sequence tag library and development of DNA microarray. Mar Biotechnol (NY). 2005, 7 (5): 481-493. 10.1007/s10126-004-4110-6.Google Scholar
- Bachvaroff TR, Concepcion GT, Rogers CR, Herman EM, Delwiche CF: Dinoflagellate expressed sequence tag data indicate massive transfer of chloroplast genes to the nuclear genome. Protist. 2004, 155 (1): 65-78. 10.1078/1434461000165.PubMedGoogle Scholar
- Tanikawa N, Akimoto H, Ogoh K, Chun W, Ohmiya Y: Expressed sequence tag analysis of the dinoflagellate Lingulodinium polyedrum during dark phase. Photochem Photobiol. 2004, 80: 31-35. 10.1562/2004-03-12-RA-110.1.PubMedGoogle Scholar
- Adl SM, Simpson AG, Farmer MA, Andersen RA, Anderson OR, Barta JR, Bowser SS, Brugerolle G, Fensome RA, Fredericq S: The new higher level classification of eukaryotes with emphasis on the taxonomy of protists. J Eukaryot Microbiol. 2005, 52 (5): 399-451. 10.1111/j.1550-7408.2005.00053.x.PubMedGoogle Scholar
- Zhang H, Hou Y, Miranda L, Campbell DA, Sturm NR, Gaasterland T, Lin S: Spliced leader RNA trans-splicing in dinoflagellates. Proc Natl Acad Sci USA. 2007, 104 (11): 4618-4623. 10.1073/pnas.0700258104.PubMed CentralPubMedGoogle Scholar
- Slamovits CH, Keeling PJ: Widespread recycling of processed cDNAs in dinoflagellates. Curr Biol. 2008, 18 (13): R550-552. 10.1016/j.cub.2008.04.054.PubMedGoogle Scholar
- Zhang H, Lin S: Retrieval of missing spliced leader in dinoflagellates. PLoS One. 2009, 4 (1): e4129-10.1371/journal.pone.0004129.PubMed CentralPubMedGoogle Scholar
- Wan KL, Blackwell JM, Ajioka JW: Toxoplasma gondii expressed sequence tags: insight into tachyzoite gene expression. Mol Biochem Parasitol. 1996, 75 (2): 179-186. 10.1016/0166-6851(95)02524-3.PubMedGoogle Scholar
- Strong WB, Nelson RG: Preliminary profile of the Cryptosporidium parvum genome: an expressed sequence tag and genome survey sequence analysis. Mol Biochem Parasitol. 2000, 107 (1): 1-32. 10.1016/S0166-6851(99)00225-X.PubMedGoogle Scholar
- Wan KL, Chong SP, Ng ST, Shirley MW, Tomley FM, Jangi MS: A survey of genes in Eimeria tenella merozoites by EST sequencing. Int J Parasitol. 1999, 29 (12): 1885-1892. 10.1016/S0020-7519(99)00160-5.PubMedGoogle Scholar
- Asojo OA, Schott EJ, Vasta GR, Silva AM: Structures of PmSOD1 and PmSOD2, two superoxide dismutases from the protozoan parasite Perkinsus marinus. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2006, 62 (Pt 11): 1072-1075. 10.1107/S1744309106040425.PubMed CentralPubMedGoogle Scholar
- Ahmed H, Schott EJ, Gauthier JD, Vasta GR: Superoxide dismutases from the oyster parasite Perkinsus marinus: purification, biochemical characterization, and development of a plate microassay for activity. Anal Biochem. 2003, 318 (1): 132-141. 10.1016/S0003-2697(03)00192-1.PubMedGoogle Scholar
- Schott EJ, Vasta GR: The PmSOD1 Gene of the protistan parasite Perkinsus marinus complements the sod2 Δ mutant of Saccharomyces cerevisiae, and directs an iron superoxide dismutase to mitochondria. Mol Biochem Parasitol. 2003, 126 (1): 81-92. 10.1016/S0166-6851(02)00271-2.PubMedGoogle Scholar
- Schott EJ, Robledo JAF, Wright AC, Silva AM, Vasta GR: Gene organization and homology modeling of two iron superoxide dismutases of the early branching protist Perkinsus marinus. Gene. 2003, 309 (1): 1-9. 10.1016/S0378-1119(03)00469-4.PubMedGoogle Scholar
- Schott EJ, Pecher WT, Okafor F, Vasta GR: The protistan parasite Perkinsus marinus is resistant to selected reactive oxygen species. Exp Parasitol. 2003, 105 (3-4): 232-240. 10.1016/j.exppara.2003.12.012.PubMedGoogle Scholar
- Wright A, Ahmed H, Gauthier J, Silva A, Vasta G: cDNA cloning and characterization of two iron superoxide dismutases from the oyster parasite Perkinsus marinus. Mol Biochem Parasitol. 2002, 123 (1): 73-10.1016/S0166-6851(02)00090-7.PubMedGoogle Scholar
- Matsuzaki M, Kita K, Nozaki H: Orthologs of plastid isoplenoids biosynthesis pathway genes from Perkinsus. J Plant Res. 2006, 119 (suppl): 122-Google Scholar
- Grauvogel C, Reece KS, Brinkmann H, Petersen J: Plastid isoprenoid metabolism in the oyster parasite Perkinsus marinus connects dinoflagellates and malaria pathogens-new impetus for studying alveolates. J Mol Evol. 2007, 65 (6): 725-729. 10.1007/s00239-007-9053-5.PubMedGoogle Scholar
- Matsuzaki M, Kuroiwa H, Kuroiwa T, Kita K, Nozaki H: A cryptic algal group unveiled: a plastid biosynthesis pathway in the oyster parasite Perkinsus marinus. Mol Biol Evol. 2008, 25 (6): 1167-1179. 10.1093/molbev/msn064.PubMedGoogle Scholar
- Earnhart CG, Kaattari SL: Potential novel epitopes in the extracellular products of oyster homogenate-supplemented Perkinsus marinus cells are not detected by subtractive immunization. J Parasitol. 2005, 91 (3): 689-691. 10.1645/GE-3427RN.PubMedGoogle Scholar
- MacIntyre EA, Earnhart CG, Kaattari SL: Host oyster tissue extracts modulate in vitro protease expression and cellular differentiation in the protozoan parasite, Perkinsus marinus. Parasitology. 2003, 126 (Pt 4): 293-302. 10.1017/S003118200200286X.PubMedGoogle Scholar
- Romestand B, Corbier F, Roch P: Protease inhibitors and haemagglutinins associated with resistance to the protozoan parasite, Perkinsus marinus, in the Pacific oyster, Crassostrea gigas. Parasitology. 2002, 125 (Pt 4): 323-329.PubMedGoogle Scholar
- Faisal M, MacIntyre EA, Adham KG, Tall BD, Kothary MH, La Peyre JF: Evidence for the presence of protease inhibitors in eastern (Crassostrea virginica) and Pacific (Crassostrea gigas) oysters. Comp Biochem Physiol, B. 1998, 121B (2): 161-168. 10.1016/S0305-0491(98)10084-6.Google Scholar
- Oliver JL, Lewis TD, Faisal M, Kaattari SL: Analysis of the effects of Perkinsus marinus proteases on plasma proteins of the eastern oyster (Crassostrea virginica) and the Pacific oyster (Crassostrea gigas). J Invertebr Pathol. 1999, 74 (2): 173-183. 10.1006/jipa.1999.4879.PubMedGoogle Scholar
- La Peyre JF, Schafhauser DY, Rizkalla EH, Faisal M: Production of serine proteases by the oyster pathogen Perkinsus marinus (Apicomplexa) in vitro. J Eukaryot Microbiol. 1995, 42 (5): 544-551. 10.1111/j.1550-7408.1995.tb05903.x.Google Scholar
- Xue QG, Waldrop GL, Schey KL, Itoh N, Ogawa M, Cooper RK, Losso JN, La Peyre JF: A novel slow-tight binding serine protease inhibitor from eastern oyster (Crassostrea virginica) plasma inhibits perkinsin, the major extracellular protease of the oyster protozoan parasite Perkinsus marinus. Comp Biochem Physiol B Biochem Mol Biol. 2006, 145 (1): 16-26. 10.1016/j.cbpb.2006.05.010.PubMedGoogle Scholar
- Fernández-Robledo JA, Schott EJ, Vasta GR: Perkinsus marinus superoxide dismutase 2 (PmSOD2) localizes to single-membrane subcellular compartments. Biochem Biophys Res Commun. 2008, 375 (2): 215-219. 10.1016/j.bbrc.2008.07.162.PubMedGoogle Scholar
- Soudant P, Chu FL, Marty Y: Lipid class composition of the protozoan Perkinsus marinus, an oyster parasite, and its metabolism of a fluorescent phosphatidylcholine analog. Lipids. 2000, 35 (12): 1387-1395. 10.1007/s11745-000-0656-1.PubMedGoogle Scholar
- Chu FL, Lund E, Soudant P, Harvey E: De novo arachidonic acid synthesis in Perkinsus marinus, a protozoan parasite of the eastern oyster Crassostrea virginica. Mol Biochem Parasitol. 2002, 119 (2): 179-190. 10.1016/S0166-6851(01)00413-3.PubMedGoogle Scholar
- Chu FL, Lund ED, Harvey E, Adlof R: Arachidonic acid synthetic pathways of the oyster protozoan parasite, Perkinsus marinus: evidence for usage of a delta-8 pathway. Mol Biochem Parasitol. 2004, 133 (1): 45-51. 10.1016/j.molbiopara.2003.08.012.PubMedGoogle Scholar
- Lund ED, Chu FL, Soudant P, Harvey E: Perkinsus marinus, a protozoan parasite of the eastern oyster, has a requirement for dietary sterols. Comp Biochem Physiol A Mol Integr Physiol. 2007, 146 (1): 141-147. 10.1016/j.cbpa.2006.10.004.PubMedGoogle Scholar
- Yu M, Kumar T, Nkrumah L, Coppi A, Retzlaff S, Li C, Kelly B, Moura P, Lakshmanan V, Freundlich J: The Fatty Acid Biosynthesis Enzyme FabI Plays a Key Role in the Development of Liver-Stage Malarial Parasites. Cell host & microbe. 2008, 4 (6): 567-578.Google Scholar
- Gardner MJ, Bishop R, Shah T, Villiers EP, Carlton JM, Hall N, Ren Q, Paulsen LT, Pain A, Berriman M: Genome sequences of Theileria parva, a bovine pathogen that transforms lymphocytes. Science. 2005, 309 (5731): 134-137. 10.1126/science.1110439.PubMedGoogle Scholar
- Weiss LM, Ma YF, Takvorian PM, Tanowitz HB, Wittner M: Bradyzoite development in Toxoplasma gondii and the hsp70 stress response. Infect Immun. 1998, 66 (7): 3295-3302.PubMed CentralPubMedGoogle Scholar
- Lyons RE, Johnson AM: Gene sequence and transcription differences in 70 kDa heat shock protein correlate with murine virulence of Toxoplasma gondii. Int J Parasitol. 1998, 28 (7): 1041-1051. 10.1016/S0020-7519(98)00074-5.PubMedGoogle Scholar
- Ivens AC, Lewis SM, Bagherzadeh A, Zhang L, Chan HM, Smith DF: A physical map of the Leishmania major Friedlin genome. Genome Res. 1998, 8 (2): 135-145.PubMed CentralPubMedGoogle Scholar
- Sunila I, Hamilton RM, Dungan CF: Ultrastructural characteristics of the in vitro cell cycle of the protozoan pathogen of oysters, Perkinsus marinus. J Eukaryot Microbiol. 2001, 48 (3): 348-361. 10.1111/j.1550-7408.2001.tb00324.x.PubMedGoogle Scholar
- Vasta GR, Sullivan JT, Cheng TC, Marchalonis JJ, Warr GW: A cell membrane-associated lectin of the oyster hemocyte. J Invertbr Pathol. 1982, 40: 367-377. 10.1016/0022-2011(82)90175-6.Google Scholar
- Vasta GR, Cheng TC, Marchalonis JJ: A lectin on the hemocyte membrane of the oyster (Crassostrea virginica). Cell Immunol. 1984, 88 (2): 475-488. 10.1016/0008-8749(84)90179-5.PubMedGoogle Scholar
- Gauthier JD, Jenkins JA, La Peyre JF: Flow cytometric analysis of lectin binding to in vitro -cultured Perkinsus marinus surface carbohydrates. J Parasitol. 2004, 90 (3): 446-454. 10.1645/GE-3269.PubMedGoogle Scholar
- Sato K, Nakano A: Emp47p and its close homolog Emp46p have a tyrosine-containing endoplasmic reticulum exit signal and function in glycoprotein secretion in Saccharomyces cerevisiae. Mol Biol Cell. 2002, 13 (7): 2518-2532. 10.1091/mbc.E02-01-0027.PubMed CentralPubMedGoogle Scholar
- Chai Y, Pang Y, Liao Z, Zhang L, Sun X, Lu Y, Wang S, Tang K: Molecular cloning and characterization of a mannose-binding lectin gene from Crinum asiaticum. J Plant Physiol. 2003, 160 (8): 913-920. 10.1078/0176-1617-01115.PubMedGoogle Scholar
- Schirmer RH, Schollhammer T, Eisenbrand G, Krauth-Siegel RL: Oxidative stress as a defense mechanism against parasitic infections. Free Radic Res Commun. 1987, 3 (1-5): 3-12. 10.3109/10715768709069763.PubMedGoogle Scholar
- Wanidworanun C, Nagel RL, Shear HL: Antisense oligonucleotides targeting malarial aldolase inhibit the asexual erythrocytic stages of Plasmodium falciparum. Mol Biochem Parasitol. 1999, 102 (1): 91-101. 10.1016/S0166-6851(99)00087-0.PubMedGoogle Scholar
- Misset O, Bos OJ, Opperdoes FR: Glycolytic enzymes of Trypanosoma brucei. Simultaneous purification, intraglycosomal concentrations and physical properties. Eur J Biochem. 1986, 157 (2): 441-453. 10.1111/j.1432-1033.1986.tb09687.x.PubMedGoogle Scholar
- Fichera ME, Roos DS: A plastid organelle as a drug target in apicomplexan parasites. Nature. 1997, 390 (6658): 407-409. 10.1038/37132.PubMedGoogle Scholar
- McFadden GI, Roos DS: Apicomplexan plastids as drug targets. Trends Microbiol. 1999, 7 (8): 328-333. 10.1016/S0966-842X(99)01547-4.PubMedGoogle Scholar
- Soldati D: The apicoplast as a potential therapeutic target in and other apicomplexan parasites. Parasitol Today. 1999, 15 (1): 5-7. 10.1016/S0169-4758(98)01363-5.PubMedGoogle Scholar
- Wiesner J, Reichenberg A, Heinrich S, Schlitzer M, Jomaa H: The plastid-like organelle of apicomplexan parasites as drug target. Curr Pharm Des. 2008, 14 (9): 855-871. 10.2174/138161208784041105.PubMedGoogle Scholar
- Wiesner J, Jomaa H: Isoprenoid biosynthesis of the apicoplast as drug target. Curr Drug Targets. 2007, 8 (1): 3-13. 10.2174/138945007779315551.PubMedGoogle Scholar
- Jomaa H, Wiesner J, Sanderbrand S, Altincicek B, Weidemeyer C, Hintz M, Turbachova I, Eberl M, Zeidler J, Lichtenthaler HK: Inhibitors of the nonmevalonate pathway of isoprenoid biosynthesis as antimalarial drugs. Science. 1999, 285 (5433): 1573-1576. 10.1126/science.285.5433.1573.PubMedGoogle Scholar
- Moreno SN, Li ZH: Anti-infectives targeting the isoprenoid pathway of Toxoplasma gondii. Expert Opin Ther Targets. 2008, 12 (3): 253-263. 10.1517/14728220.127.116.11.PubMedGoogle Scholar
- de Azevedo WF, Soares MB: Selection of targets for drug development against protozoan parasites. Curr Drug Targets. 2009, 10 (3): 193-201. 10.2174/138945009787581186.PubMedGoogle Scholar
- Yamagishi H, Garg V, Matsuoka R, Thomas T, Srivastava D: A molecular pathway revealing a genetic basis for human cardiac and craniofacial defects. Science. 1999, 283 (5405): 1158-1161. 10.1126/science.283.5405.1158.PubMedGoogle Scholar
- Dessens JT, Margos G, Rodriguez MC, Sinden RE: Identification of differentially regulated genes of Plasmodium by suppression subtractive hybridization. Parasitol Today. 2000, 16 (8): 354-356. 10.1016/S0169-4758(00)01710-5.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.