- Research article
- Open Access
Generation and analysis of expressed sequence tags from the ciliate protozoan parasite Ichthyophthirius multifiliis
© Abernathy et al; licensee BioMed Central Ltd. 2007
- Received: 21 November 2006
- Accepted: 18 June 2007
- Published: 18 June 2007
The ciliate protozoan Ichthyophthirius multifiliis (Ich) is an important parasite of freshwater fish that causes 'white spot disease' leading to significant losses. A genomic resource for large-scale studies of this parasite has been lacking. To study gene expression involved in Ich pathogenesis and virulence, our goal was to generate expressed sequence tags (ESTs) for the development of a powerful microarray platform for the analysis of global gene expression in this species. Here, we initiated a project to sequence and analyze over 10,000 ESTs.
We sequenced 10,368 EST clones using a normalized cDNA library made from pooled samples of the trophont, tomont, and theront life-cycle stages, and generated 9,769 sequences (94.2% success rate). Post-sequencing processing led to 8,432 high quality sequences. Clustering analysis of these ESTs allowed identification of 4,706 unique sequences containing 976 contigs and 3,730 singletons. These unique sequences represent over two million base pairs (~10% of Plasmodium falciparum genome, a phylogenetically related protozoan). BLASTX searches produced 2,518 significant (E-value < 10-5) hits and further Gene Ontology (GO) analysis annotated 1,008 of these genes. The ESTs were analyzed comparatively against the genomes of the related protozoa Tetrahymena thermophila and P. falciparum, allowing putative identification of additional genes. All the EST sequences were deposited by dbEST in GenBank (GenBank: EG957858–EG966289). Gene discovery and annotations are presented and discussed.
This set of ESTs represents a significant proportion of the Ich transcriptome, and provides a material basis for the development of microarrays useful for gene expression studies concerning Ich development, pathogenesis, and virulence.
- Gene Ontology
- Enzyme Commission
- Infected Fish
- BLASTX Search
- Normalize cDNA Library
The ciliate protozoan Ichthyophthirius multifiliis (Ich)is one of the most devastating pathogens. It infects fish skin and gills, and causes white spot diseases in many species of freshwater fish worldwide, which leads to significant losses in the aquaculture industry. The ciliate parasite has three main life-cycle stages: the reproductive tomont, the infective theront, and a parasitic trophont [1–3]. The mature trophont drops off the host to become the tomont where it attaches to a substrate, and undergoes multiple divisions to produce hundreds to thousands of tomites within a cyst. Tomites bore their way through the cyst into water, and differentiate into theronts that infect fish. Once they burrow into the fish epithelium, theronts become trophonts that feed and mature in the host.
In spite of great losses caused by Ich to the aquaculture industry, molecular studies of the parasite have been scarce [see a recent review ]. Limited studies have concentrated on immune responses of the host and factors affecting them [5–11]. One of the difficulties for the studies of Ich is the problem involved in long-term maintenance of Ich isolates. Ich isolates appear to lose infectivity or become senescent after a certain number of passages [12–15]. Most often a significant decrease in infectivity is observed after about 50 passages . Not only the infectivity decreases with higher numbers of passages, but also the development of the parasite as measured by the period required for trophonts to emerge from fish .
The Ich senescence phenomenon is interesting not only as a developmental biology issue, but also as a potential research system to study the virulence factors involved in the parasite pathogenesis. Assuming the life cycles of Ich and its infectivity are controlled by gene products, then it would be of great interest to learn what genes are involved in the loss of infectivity, and in the slowing down of its development. However, as very limited molecular information is available from Ich, in-depth research is limited by the lack of information and the lack of genomic resources.
EST analysis is one of the most effective means for gene discoveries, gene expression profiling, and functional genome studies [16–23]. It is also one of the most efficient ways for the identification of differentially expressed genes [24–28]. In order to provide genomic resources for the analysis of differentially expressed genes at different developmental stages of the Ich parasite, and for the analysis of genes differentially expressed when infectivity is being lost, the objectives of this study were to create cDNA libraries suitable for the analysis of expressed sequence tags (ESTs) and to generate an EST resource for Ich to allow cDNA-based design of microarrays for the study of gene expression in relation to the passages and development of the parasite. Before this work, there were only 511 Ich sequences in the GenBank dbEST (release 100606) . A brief examination of these existing EST sequences indicated that a large proportion of them were trophont only reads, histones, ribosomal proteins, and immobilization antigen-related sequences. Here we report sequencing of 10,368 Ich EST clones, and generation of 8,432 high quality EST sequences. This EST resource should provide the material basis for the development of microarrays for Ich, and serve as a platform for its functional genomic studies including the development and pathogenesis of Ich, and the host-parasite interactions.
Generation of the Ich ESTs
A summary of the EST analysis.
Total number of clones sequenced
Total number of successful sequences
Number of high quality sequences
Number of contigs
Number of clones included in the contigs
Average clones per contig
Number of singletons
Number of known genes
Unique unknown genes
The processed sequences were subjected to cluster analysis using Vector NTI to evaluate sequence redundancies. Of the 8,432 sequences, 4,702 sequences fell within 976 contigs while 3,730 sequences were singletons. On average, each contig contained 4.8 sequences. Taken together, the 976 contigs and the 3,730 singletons made up 4,706 unique sequences (Table 1).
The most abundant ESTs detected from the EST sequencing
# of Sequences
% of Total
Hypothetical protein TTHERM_02141640 from Tetrahymena thermophila
Hypothetical protein TTHERM_02641280 from Tetrahymena thermophila
Heat shock protein 90
Zinc finger ZZ type family protein
Heat shock protein 90
Heat shock protein 70 (dnaK)
Conserved hypothetical protein from Paracoccus denitrificans
Hypothetical protein PY05925 from Plasmodium yoelii
ER type HSP70
Unknown protein from Oryza sativa
Dynein heavy chain protein
Outer surface protein from Rickettsia typhi
This work demonstrated that pooling of samples from all three stages of Ich life cycle followed by normalization was an effective way to reduce common messages across all three life stages. As one would expect, many structural genes would be expressed highly abundantly in all stages of the life cycle. In addition to making savings economically, pooling of samples allowed very effective normalization of these common transcripts without going through three rounds of normalization. This is consistent with our previous experience for the generation of a large number of catfish and oyster ESTs [32–36]. It is obvious that the pooling of samples from three developmental stages made it impossible to provide information concerning expression profiling in relation to developmental stages. However, such information would not be highly meaningful in normalized cDNA libraries where the major focus was to develop EST resources, rather than expression profiling. The other limitation caused by construction of a pooled cDNA library is the loss of sequencing flexibility as to the number of clones to be sequenced from each developmental stage library if they had been separately constructed.
A summary of simple sequence repeats identified from the Ich ESTs. Percentages indicated in the parentheses are percentage of each type of repeat among all repeats
Total number of sequences analyzed
Number of dinucleotide repeats
Number of AC repeats
Number of AG repeats
Number of AT repeats
Number of CT repeats
Number of GT repeats
Number of GC repeats
Number of trinucleotide repeats
Number of tetranucleotide repeats
Total number simple sequence repeats
The putative identities of the sequenced ESTs were assessed using BLASTX searches against the non-redundant (NR) database in GenBank . All the search results are summarized in supplemental Table 1. Of the 4,706 unique ESTs, 2,518 (53.5%) had significant (E-value < 10-5) hits. The remaining 2,188 (46.5%) EST sequences were not similar to any known sequences. Additional searches using the Swiss-Prot database resulted in putative identities for six additional unknown ESTs (Supplemental Table 1).
Identification of putative secretory proteins
Secretory proteins have been shown to be an important component in many biological processes, including pathogenesis of parasites [40–42]. We therefore searched for transcripts with putative signal peptides (suggestive of peptides of secretory proteins) within the EST set using the program SignalP 3.0 . We found 314 ESTs with signal peptides, representing 6.7% of the unique sequences. Of these, 180 (3.8%) were from ESTs with no significant (E-value < 10-5) BLASTX hit to the NR database in GenBank (Supplemental Table 1).
Comparative analysis to related taxa
The parasite Ich is phylogenetically placed between the protozoan's Plasmodium falciparum and Tetrahymena thermophila. Previous studies using 18S rDNA, histone genes, and I-antigens [44–46] suggested that Ich was more related to T. thermophila than to P. falciparum. Furthermore, T. thermophila and Ich share the ciliate nuclear genetic code, while P. falciparum uses the standard genetic code for translation. As the entire genome sequence of P. falciparum is available and the macronuclear sequencing project of T. thermophila was just recently completed, we made comparative BLAST analyses against both genome sequences.
Of the 1,759 significant hits against the T. thermophila genome, 1,673 had been identified with a putative identity using BLASTX searches against the NR database, while the tBLASTx searches against the T. thermophila genome allowed identification of putative identities for additional 86 unique ESTs. Similarly, BLASTX searches against the P. falciparum genome allowed identification of 9 additional ESTs. Taken together, the BLAST searches against these two genomes allowed putative identities of 95 additional unique ESTs, bringing the total number of ESTs with significant similarities to known genes to 2,613.
Such genome searches also revealed that of the 2,518 ESTs that had significant hits in BLASTX searches against the NR database, 845 had no significant hits to the Tetrahymena genome. Clearly, these ESTs were similar to sequences of organisms other than the ciliate Tetrahymena.
These results clearly suggest conservation of a large fraction of gene sequences among the three protozoa parasites, with a higher level of conservation between Ich and the T. thermophila genome than between the Ich genome and the P. falciparum genome; although a significant fraction of gene sequences are also shared between the genomes of T. thermophila and P. falciparum. The results of this comparative analysis are compatible with existing phylogenetic analyses using several molecular markers such as the 18S rDNA, histone genes, and the I-antigens. Obviously, use of a large set of sequences should provide a greater confidence concerning genome evolution. The comparative analysis suggested that the EST resource generated from this study should be useful for phylogenetic analysis and studies concerning genome evolution.
The unique Ich sequences were compared to annotations through the Gene Ontology Consortium  using the automated software Blast2GO . We were able to obtain GO terms for 1,008 unique sequences using this method. Of these, 304 were contigs and 704 were singletons. Sequence descriptions, gene ontology (GO) and enzyme commission (EC) numbers are summarized in Supplemental Table 3. There were 258 sequences with both GO terms and EC numbers.
We have produced 8,432 high quality I. multifiliis EST sequences. Sequence analysis indicated the presence of 4,706 unique sequences in the EST set. This should represent a significant fraction of the Ich genes, although the exact gene number of the parasite is unknown at present. The majority of the unique EST sequences had similarities to known genes, making them more amenable to functional analysis. The EST sequences should enhance the effectiveness of molecular studies, especially for gene expression profiling and the analysis of genes involved in virulence and infectivity. Microarrays can now be designed using either cDNA microarray or oligo-based platforms using the EST information. Additionally, the cluster and redundancy information should be useful for further subtraction of the most abundant transcripts included in the cDNA library, making further EST analysis in the parasite more effective.
The source of mRNA for this analysis was derived and expanded from a single parasite cloned from the infected fish. The source of the original Ichthyophthirius multifiliis was isolated from an infected fish obtained from a local pet shop and the parasite was transmitted to channel catfish held in a 50-l glass aquarium at the USDA-ARS Aquatic Animal Health Research Laboratory, Auburn, AL. The transmission of I. multifiliis was achieved through co-habitation of the infected fish with two fingerlings of channel catfish (3 inches in size). When the two catfish were infected, the symptoms of Ich, white spots, started to emerge when trophonts were collected by scraping with a glass slide. Channel catfish infected with maturing trophonts were rinsed in dechlorinated water and the skin was gently scraped to dislodge the parasites. Trophonts were harvested by filtering through a 0.22 μ m filter to remove fish skin. The trophonts were placed into a Petri dish to allow them to develop into theronts that were used to infect 8 fish each for the collection of trophonts, toments, and theronts, respectively. Trophonts were directly collected from the skin surface of the 8 infected fish. To collect tomonts and theronts, trophonts isolated from fish were placed in Petri dishes and allowed to attach. After replacing the water in the Petri dishes with fresh dechlorinated water to remove contaminating mucus, the trophonts were incubated at 24°C for 8 h to harvest tomonts (32–128 cells/cyst) or 24 h to harvest theronts. Trophonts, tomonts and theronts were washed with PBS (pH 7.4), concentrated with a centrifuge (Beckman Coulter, Inc., Miami, FL) at 228 × g for 5 min and discarded supernatant. After washing 3 times with PBS, parasite samples from the three life stages were stored in liquid nitrogen and used for the isolation of RNA for the construction of normalized cDNA library.
Total RNA was isolated from the samples using the TRIzol reagent method from Invitrogen (Carlsbad, CA, USA) according to manufacturer's instructions. Briefly, samples of tomont, theront, and trophont were resuspended after thawing on ice, and 100 μl each were combined in a sterile tube to provide a total of 300 μl of Ich samples with equal fractions from each of its three life stages. As the major objective of this study was to generate EST resources with maximal efficiency of gene discovery, a pooled sample followed with normalization would allow inclusion of all transcripts in the library while reducing cost for library construction and increasing gene discovery rate. Three milliliters of TRIzol reagent was added to the sample tube. Cells were lysed by repetitively pipetting up and down. RNA was isolated following the manufacturer's protocol. The RNA pellet was resuspended in 100 μl of RNase-free double distilled water and divided into 25 μl aliquots. RNA aliquots were checked for quality using agarose gel electrophoreses containing formaldehyde.
Normalized library construction
The Creator Smart cDNA Library Construction Kit from Clontech (Mountain View, CA) and components from the TRIMMER-DIRECT Kit from Evrogen (Moscow, Russia) were used for the construction of the normalized cDNA library. Total RNA concentration was checked on a spectrophotometer and 1 μg RNA was combined with 1 μl of SMART IV oligonucleotide (Clontech) and 1 μl CDS-3M adapter (Evrogen) for first-strand cDNA synthesis. The reaction was incubated at 72°C for 2 min followed by immediate cooling on ice for 2 min. Next, 2 μl of 5× first strand buffer, 1 μl of DTT (20 mM), 1 μl of dNTP mix (10 mM), and 1 μl of PowerScript reverse transcriptase were added to the tube and incubated at 42°C for 1 h in a thermal cycler (PTC-100, Bio-Rad, Hercules, CA) then placed on ice. The SMART cDNA cloning system allows the enrichment of full-length cDNA through the use of a 5'-linker with 3'-GGG tails. Reverse transcriptase has terminal transferase activity that preferentially adds three additional Cs at the end of first strand cDNA. As a result, the first strand cDNA is able to base pair with the 5'-linker with 3'-GGG tails. Once base paired, the reverse transcriptase would switch the template and extend into the linker sequences allowing PCR amplification of full-length cDNA using a single primer (the 5'-linker has the same sequences as the linker containing poly T used for the synthesis of the first strand cDNA). Truncated cDNAs are not able to base pair with the 5'-linker and, therefore, get lost in the PCR amplification of the full-length cDNA.
The first strand cDNA was initially amplified by long-distance PCR (LD-PCR) using hot-start amplification. For the reaction, the following were combined in a reaction tube: 1.5 μl of the first strand cDNA, 60 μl of sterile deionized water, 7.5 μl of 10× Advantage 2 PCR buffer, 1.5 μl of 50× dNTP mix, 3 μl of 5' PCR primer and 1.5 μl of 50× Advantage 2 polymerase mix. The tube was mixed and briefly centrifuged and added to a pre-heated (95°C) thermal cycler. Cycle settings were 95°C for 1 min followed by 19 cycles of 95°C for 7 s, 66°C for 20 s, and 72°C for 5.5 min. The product was analyzed on a 1.1% agarose gel to determine the sizes and amount of the cDNA products before proceeding to the next step. The LD-PCR reaction was purified and eluted in 30 μl of sterile Nanopure water using the QIAquick PCR Purification Kit (Qiagen, Valencia, CA).
For the normalization procedure, the TRIMMER-DIRECT Kit from Evrogen (Moscow, Russia) was used. This system is specially developed to normalize cDNA enriched with full length sequences [49, 50]. The cDNA from the LD-PCR was quantified [~100 ng/μl] and 1 μl was mixed with 1 μl of the 4× hybridization buffer and 2 μl of sterile water. The mix was overlaid with mineral oil and incubated for 3 min at 98°C followed by 4 h at 70°C. Then, 5 μl of 2× DSN buffer (preheated to 70°C) and 0.25 Kunitz units of DSN enzyme were added and incubated at 70°C for 20 min. The DSN enzyme specifically degrades double-stranded molecules. The reaction was inactivated by adding 10 μl of DSN stop solution, and sterile water added to a final volume of 40 μl.
Following normalization, two rounds of PCR were performed using 1 μl of the normalization reaction as template. A shorter primer M1 (first 23 bases of the SMART IV oligonucleotide) was used in the first round of PCR with 15 amplification cycles using the same thermal cycling parameters as above; and an even shorter primer M2 (first 20 bases of the SMART IV oligonucleotide) was used in the second round of PCR for 15 amplification cycles of 95°C for 7 s, 64°C for 20 s, and 72°C for 5.5 min. Products were checked on a 1.1% agarose gel. The PCR products were quantified and 3 μg were used for treatment with proteinase K. All the subsequent procedures including proteinase K treatment, restriction digestion with Sfi I, size fractionation, and ligation followed the manufacture's instructions (Clontech). The cDNA was ligated to the pDNR-LIB vector. Electroporation (MicroPulser, Bio-Rad, Hercules, CA) was performed using DH12S electrocompetent cells following supplier's instructions (Invitrogen). A total of approximately 700,000 primary recombinant clones were obtained, and the library was amplified, titered, and stored in glycerol stocks in a -80°C freezer.
Plasmid isolation and EST sequencing
Independent colonies were picked and grown for 20 h at 37°C in 1.2 ml LB broth containing 30 μg/ml chloramphenicol. Plasmid DNA was isolated using the Perfectprep Plasmid 96 Vacuum Direct Bind Kit from Eppendorf (Westbury, NY). Plasmids were stored at -20°C until usage. The cDNA inserts were directionally sequenced from the 5'-end of the cDNAs using universal M13(-21) primer and the BigDye terminator sequencing kit version 3.1 from Applied Biosystems (Foster City, CA) on a 3130XL DNA analyzer (Applied Biosystems).
Base calling was performed using the Phred program [30, 31] at quality cut-off set at Phred 20. Raw sequences were then imported into the Vector NTI Advance 10 software (Invitrogen) and were subjected to trimming of vector sequences and 5'adapter sequences using default settings. Afterwards, poly (A) tails were trimmed where necessary and sequences less than 100 bases were removed. Contigs were built in Vector NTI ContigExpress using default settings. All unique sequences were compared to the GenBank database using BLASTX in the non-redundant (NR), Swiss-Prot, and Plasmodium falciparum 3D7 genome database. For comparison to the Tetrahymena thermophila SB210 genome, tBLASTx was used. The cut-off for sequence similarity used was E-value < 10-5 for all analyses. Ciliate nuclear translation code was used in the BLAST searches. Search results from genome comparisons were summarized using a Venn diagram.
Gene ontology (GO) annotations were assigned using the program Blast2GO . BLASTX results were loaded into the program and the default settings were used to assign GO terms to all unique sequences. From these annotations, pie charts were made using 2nd level GO terms based on biological process, molecular function, and cellular component.
Putative secretory proteins and signal peptides were identified using both neural networks and hidden Markov model methods in SignalP 3.0 . All 4,706 unique ESTs were used as the tester sequences. Open reading frames were predicted using both OrfPredictor  and BLASTX, with ciliate nuclear genetic code for ESTs of known genes and just OrfPredictor with unknown ESTs. The resulting deduced protein sequences from the ORFs were uploaded into SignalP 3.0. Sequences were identified as putatively secretory, predicted with signal peptides if both D-score in the neural network model and prediction probability in the hidden Markov model were significant.
Total lengths of all ESTs, G+C% content and simple repetitive elements were estimated using the Repeatmasker program .
All Ich EST sequences were submitted to the dbEST database of NCBI. Continuous accession numbers are from EG957858–EG966289.
This project was supported by a Specific Cooperative Agreement with USDA ARS Aquatic Animal Health Laboratory under the Contract Number 58-6420-5-030, and in part by CSREES from a grant of USDA NRI Animal Genome Tools and Resources Program (award # 2006-35616-16685) to support the first year stipend of JWA. We are grateful for an equipment grant from the National Research Initiative Competitive Grant no. 2005-35206-15274 from the USDA Cooperative State Research, Education, and Extension Service.
- MacLennan RF: Observations on the life cycle of Ichthyophthirius, a ciliate parasitic on fish. Northwestern Scientist. 1935, 9: 12-14.Google Scholar
- Hines RS, Spira DT: Ichthyophthiriasis in the mirror carp Cyprinus carpio (L.) V. Acquired immunity. J Fish Biology. 1974, 6: 373-378. 10.1111/j.1095-8649.1974.tb04554.x.View ArticleGoogle Scholar
- Nigrelli RF, Pokorny KS, Ruggieri GD: Notes on Ichthyophthirius multifilis, a ciliate parasitic on fresh-water fishes, with some remarks on possible physiological races and species. Trans Ame Microscopical Soc. 1976, 95 (4): 607-613. 10.2307/3225383.View ArticleGoogle Scholar
- Matthews RA: Ichthyophthirius multifiliis Fouquet and Ichthyophthiriosis in Freshwater Teleosts. Adv parasitol. 2005, 59: 159-241.PubMedView ArticleGoogle Scholar
- Xu DH, Klesius PH: Protective effect of cutaneous antibodyproduced by channel catfish, Ictalurus punctatus (Rafinesque), immune to Ichthyophthirius multifiliis Fouquet on cohabited non-immune catfish. J fish dis. 2003, 26 (5): 287-291. 10.1046/j.1365-2761.2003.00463.x.PubMedView ArticleGoogle Scholar
- Xu DH, Klesius PH, Panangala VS: Induced cross-protectionin channel catfish, Ictalurus punctatus (Rafinesque), against different immobilization serotypes of Ichthyophthirius multifiliis. J fish dis. 2006, 29 (3): 131-138. 10.1111/j.1365-2761.2006.00700.x.PubMedView ArticleGoogle Scholar
- Xu DH, Klesius PH, Shelby RA: Immune responses and hostprotection of channel catfish, Ictalurus punctatus (Rafinesque), against Ichthyophthirius multifiliis after immunization with live theronts and sonicated trophonts. J fish dis. 2004, 27 (3): 135-141. 10.1111/j.1365-2761.2004.00523.x.PubMedView ArticleGoogle Scholar
- Xu DH, Klesius PH, Shoemaker CA: Effect of lectins on theinvasion of Ichthyophthirius theront to channel catfish tissue. Dis Aquat organisms. 2001, 45 (2): 115-120.View ArticleGoogle Scholar
- Xu DH, Klesius PH, Shoemaker CA: Cutaneous antibodies fromchannel catfish, Ictalurus punctatus (Rafinesque), immune to Ichthyophthirius multifiliis (Ich) may induce apoptosis of Ich theronts. J fish dis. 2005, 28 (4): 213-220. 10.1111/j.1365-2761.2005.00622.x.PubMedView ArticleGoogle Scholar
- Dickerson HW, Clark TG, Findly RC: Ichthyophthirius multifiliis has membrane-associated immobilization antigens. J Protozool. 1989, 36 (2): 159-164.PubMedView ArticleGoogle Scholar
- Dickerson HW, Clark TG, Leff AA: Serotypic variationamong isolates of Ichthyophthirius multifiliis based on immobilization. J Eeukaryotic Microbiol. 1993, 40 (6): 816-820. 10.1111/j.1550-7408.1993.tb04480.x.View ArticleGoogle Scholar
- Ekless L, Matthews R: Ichthyophthirius multifiliis: axenic isolation and short-term maintenance in selected monophasic media. J fish dis. 1993, 16: 437-447. 10.1111/j.1365-2761.1993.tb00877.x.View ArticleGoogle Scholar
- Burkart M, Clark T, Dickerson H: Immunization of channelcatfish, Ictalurus punctatus Rafinesque, against Ichthyophthirius multifiliis (Fouquet): killed versus live vaccines. J Fish Biol. 1990, 13: 401-410.Google Scholar
- Houghton G, Matthews RA: Immunosuppression of carp (Cyprinus carpio L.) to ichthyophthiriasis using the corticosteroid triamcinolone acetonide. Vet Immunol Immunopathol. 1986, 12 (1–4): 413-419. 10.1016/0165-2427(86)90148-0.PubMedView ArticleGoogle Scholar
- Xu DH, Klesius PH: Two year study on the infectivity of Ichthyophthirius multifiliis in channel catfish Ictalurus punctatus. Dis Aquat Organisms. 2004, 59 (2): 131-134.View ArticleGoogle Scholar
- Deng Y, Dong Y, Thodima V, Clem RJ, Passarelli AL: Analysis and functional annotation of expressed sequence tags from the fall armyworm Spodoptera frugiperda. BMC genomics. 2006, 7 (1): 264-10.1186/1471-2164-7-264.PubMed CentralPubMedView ArticleGoogle Scholar
- Govoroun M, Le Gac F, Guiguen Y: Generation of a large scale repertoire of Expressed Sequence Tags (ESTs) from normalised rainbow trout cDNA libraries. BMC genomics. 2006, 7: 196-10.1186/1471-2164-7-196.PubMed CentralPubMedView ArticleGoogle Scholar
- Hagen-Larsen H, Laerdahl JK, Panitz F, Adzhubei A, Hoyheim B: An EST-based approach for identifying genes expressed in the intestine and gills of pre-smolt Atlantic salmon (Salmo salar). BMC genomics. 2005, 6: 171-10.1186/1471-2164-6-171.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim TH, Kim NS, Lim D, Lee KT, Oh JH, Park HS, Jang GW, Kim HY, Jeon M, Choi BH: Generation and analysis of large-scale expressed sequence tags (ESTs) from a full-length enriched cDNA library of porcine backfat tissue. BMC genomics. 2006, 7: 36-10.1186/1471-2164-7-36.PubMed CentralPubMedView ArticleGoogle Scholar
- Vizcaino JA, Gonzalez FJ, Suarez MB, Redondo J, Heinrich J, Delgado-Jarana J, Hermosa R, Gutierrez S, Monte E, Llobell A: Generation, annotation and analysis of ESTs from Trichoderma harzianum CECT 2413. BMC genomics. 2006, 7: 193-10.1186/1471-2164-7-193.PubMed CentralPubMedView ArticleGoogle Scholar
- Gonzalez SF, Chatziandreou N, Nielsen ME, Li W, Rogers J, Taylor R, Santos Y, Cossins A: Cutaneous immune responses in the common carp detected using transcript analysis. Mol Immunol. 2007, 44 (7): 1675-1690. 10.1016/j.molimm.2006.08.016.View ArticleGoogle Scholar
- La Claire JW: Analysis of expressed sequence tagsfrom the harmful alga, Prymnesium parvum (Prymnesiophyceae, haptophyta). Mar Biotechnol (NY). 2006, 8 (5): 534-546. 10.1007/s10126-005-5182-2.View ArticleGoogle Scholar
- Perez F, Ortiz J, Zhinaula M, Gonzabay C, Calderon J, Volckaert FA: Development of EST-SSR markers by data mining in three species of shrimp: Litopenaeus vannamei, Litopenaeus stylirostris, and Trachypenaeus birdy. Mar Biotechnol (NY). 2005, 7 (5): 554-569. 10.1007/s10126-004-5099-1.View ArticleGoogle Scholar
- Brenner ED, Katari MS, Stevenson DW, Rudd SA, Douglas AW, Moss WN, Twigg RW, Runko SJ, Stellari GM, McCombie WR: EST analysis in Ginkgo biloba: an assessment of conserved developmental regulators and gymnosperm specific genes. BMC genomics. 2005, 6: 143-10.1186/1471-2164-6-143.PubMed CentralPubMedView ArticleGoogle Scholar
- Gandhe AS, Arunkumar KP, John SH, Nagaraju J: Analysis ofbacteria-challenged wild silkmoth, Antheraea mylitta (lepidoptera) transcriptome reveals potential immune genes. BMC genomics. 2006, 7: 184-10.1186/1471-2164-7-184.PubMed CentralPubMedView ArticleGoogle Scholar
- Hecht J, Kuhl H, Haas SA, Bauer S, Poustka AJ, Lienau J, Schell H, Stiege AC, Seitz V, Reinhardt R: Gene identification and analysis of transcripts differentially regulated in fracture healing by EST sequencing in the domestic sheep. BMC genomics. 2006, 7: 172-10.1186/1471-2164-7-172.PubMed CentralPubMedView ArticleGoogle Scholar
- Nelson RT, Shoemaker R: Identification and analysis of gene families from the duplicated genome of soybean using EST sequences. BMC genomics. 2006, 7: 204-10.1186/1471-2164-7-204.PubMed CentralPubMedView ArticleGoogle Scholar
- Ribichich KF, Georg RC, Gomes SL: Comparative EST analysis provides insights into the basal aquatic fungus Blastocladiella emersonii. BMC genomics. 2006, 7: 177-10.1186/1471-2164-7-177.PubMed CentralPubMedView ArticleGoogle Scholar
- Boguski MS, Lowe TM, Tolstoshev CM: dbEST – database for "expressed sequence tags". Nat Genet. 1993, 4 (4): 332-333. 10.1038/ng0893-332.PubMedView ArticleGoogle Scholar
- Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8 (3): 186-194.PubMedView ArticleGoogle Scholar
- Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8 (3): 175-185.PubMedView ArticleGoogle Scholar
- Kocabas AM, Li P, Cao D, Karsi A, He C, Patterson A, Ju Z, Dunham RA, Liu Z: Expression profile of the channel catfish spleen: analysis of genes involved in immune functions. Mar Biotechnol (NY). 2002, 4 (6): 526-536. 10.1007/s10126-002-0067-0.View ArticleGoogle Scholar
- Karsi A, Cao D, Li P, Patterson A, Kocabas A, Feng J, Ju Z, Mickett KD, Liu Z: Transcriptome analysis of channel catfish (Ictalurus punctatus): initial analysis of gene expression and microsatellite-containing cDNAs in the skin. Gene. 2002, 285 (1–2): 157-168. 10.1016/S0378-1119(02)00414-6.PubMedView ArticleGoogle Scholar
- Cao D, Kocabas A, Ju Z, Karsi A, Li P, Patterson A, Liu Z: Transcriptome of channel catfish (Ictalurus punctatus): initial analysis of genes and expression profiles of the head kidney. Anim Genet. 2001, 32 (4): 169-188. 10.1046/j.1365-2052.2001.00753.x.PubMedView ArticleGoogle Scholar
- Ju Z, Karsi A, Kocabas A, Patterson A, Li P, Cao D, Dunham R, Liu Z: Transcriptome analysis of channel catfish (Ictalurus punctatus): genes and expression profile from the brain. Gene. 2000, 261 (2): 373-382. 10.1016/S0378-1119(00)00491-1.PubMedView ArticleGoogle Scholar
- Quilang J, Wang S, Li P, Abernathy J, Peatman E, Wang Y, Wang L, Shi Y, Wallace R, Guo X, Liu Z: Generation and analysis ofESTs from the eastern oyster, Crassostrea virginica Gmelin and identification of microsatellite and SNP markers. BMC Genomics. 2007,Google Scholar
- Wuitschick JD, Karrer KM: Analysis of genomic G + C content, codon usage, initiator codon context and translation termination sites in Tetrahymena thermophila. J Eukaryotic Microbiol. 1999, 46 (3): 239-247. 10.1111/j.1550-7408.1999.tb05120.x.View ArticleGoogle Scholar
- Eisen JA, Coyne RS, Wu M, Wu D, Thiagarajan M, Wortman JR, Badger JH, Ren Q, Amedeo P, Jones KM: MacronuclearGenome Sequence of the Ciliate Tetrahymena thermophila, a Model Eukaryote. PLoS Biol. 2006, 4 (9): e286-10.1371/journal.pbio.0040286.PubMed CentralPubMedView ArticleGoogle Scholar
- The National Center for Biotechnology Information. [http://www.ncbi.nlm.nih.gov/]
- O'Donnell RA, Blackman MJ: The role of malaria merozoite proteases in red blood cell invasion. Curr Opinion Microbiol. 2005, 8 (4): 422-427. 10.1016/j.mib.2005.06.018.View ArticleGoogle Scholar
- Prato M, Giribaldi G, Polimeni M, Gallo V, Arese P: Phagocytosis of hemozoin enhances matrix metalloproteinase-9 activity and TNF-alpha production in human monocytes: role of matrix metalloproteinases in the pathogenesis of falciparum malaria. J Immunol. 2005, 175 (10): 6436-6442.PubMedView ArticleGoogle Scholar
- Tosini F, Trasarti E, Pozio E: Apicomplexa genes involvedin the host cell invasion: the Cpa135 protein family. Parassitologia. 2006, 48 (1–2): 105-107.PubMedGoogle Scholar
- Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 340 (4): 783-795. 10.1016/j.jmb.2004.05.028.PubMedView ArticleGoogle Scholar
- Wright AD, Lynn DH: Phylogeny of the fish parasite Ichthyophthirius and its relatives Ophryoglena and Tetrahymena (Ciliophora, Hymenostomatia) inferred from 18S ribosomal RNA sequences. Mol Biol Evol. 1995, 12 (2): 285-290.PubMedGoogle Scholar
- Van Den Bussche RA, Hoofer SR, Drew CP, Ewing MS: Characterization of histone H3/H4 gene region and phylogenetic affinity of Ichthyophthirius multifiliis based on H4 DNA sequence variation. Mol Phylogenet Evol. 2000, 14 (3): 461-468. 10.1006/mpev.1999.0700.PubMedView ArticleGoogle Scholar
- Clark TG, McGraw RA, Dickerson HW: Developmental expression of surface antigen genes in the parasitic ciliate Ichthyophthirius multifiliis. Proc Natl Acad Sci USA. 1992, 89 (14): 6363-6367. 10.1073/pnas.89.14.6363.PubMed CentralPubMedView ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.PubMed CentralPubMedView ArticleGoogle Scholar
- Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics (Oxford, England). 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.View ArticleGoogle Scholar
- Zhulidov PA, Bogdanova EA, Shcheglov AS, Vagner LL, Khaspekov GL, Kozhemyako VB, Matz MV, Meleshkevitch E, Moroz LL, Lukyanov SA: Simple cDNA normalization using kamchatka crab duplex-specific nuclease. Nucl Acids Res. 2004, 32 (3): e37-10.1093/nar/gnh031.PubMed CentralPubMedView ArticleGoogle Scholar
- Liu Z: Transcriptome characterization through the generation and analysis of expressed sequence tags: Factors to consider for a successful EST project. Israel J Aquaculture – Bamidgeh. 2006, 58 (4): 328-340.Google Scholar
- Min XJ, Butler G, Storms R, Tsang A: OrfPredictor:predicting protein-coding regions in EST-derived sequences. Nucl Acids Res. 2005, W677-680. 10.1093/nar/gki394. 33 Web ServerGoogle Scholar
- Smit A, Hubley R, Green P: RepeatMasker Open – 3.0. 1996, [http://www.repeatmasker.org]Google 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.