Expressed sequence tags (ESTs) from immune tissues of turbot (Scophthalmus maximus) challenged with pathogens
© Pardo et al; licensee BioMed Central Ltd. 2008
Received: 23 June 2008
Accepted: 25 September 2008
Published: 25 September 2008
The turbot (Scophthalmus maximus; Scophthalmidae; Pleuronectiformes) is a flatfish species of great relevance for marine aquaculture in Europe. In contrast to other cultured flatfish, very few genomic resources are available in this species. Aeromonas salmonicida and Philasterides dicentrarchi are two pathogens that affect turbot culture causing serious economic losses to the turbot industry. Little is known about the molecular mechanisms for disease resistance and host-pathogen interactions in this species. In this work, thousands of ESTs for functional genomic studies and potential markers linked to ESTs for mapping (microsatellites and single nucleotide polymorphisms (SNPs)) are provided. This information enabled us to obtain a preliminary view of regulated genes in response to these pathogens and it constitutes the basis for subsequent and more accurate microarray analysis.
A total of 12584 cDNAs partially sequenced from three different cDNA libraries of turbot (Scophthalmus maximus) infected with Aeromonas salmonicida, Philasterides dicentrarchi and from healthy fish were analyzed. Three immune-relevant tissues (liver, spleen and head kidney) were sampled at several time points in the infection process for library construction. The sequences were processed into 9256 high-quality sequences, which constituted the source for the turbot EST database. Clustering and assembly of these sequences, revealed 3482 different putative transcripts, 1073 contigs and 2409 singletons. BLAST searches with public databases detected significant similarity (e-value ≤ 1e-5) in 1766 (50.7%) sequences and 816 of them (23.4%) could be functionally annotated. Two hundred three of these genes (24.9%), encoding for defence/immune-related proteins, were mostly identified for the first time in turbot. Some ESTs showed significant differences in the number of transcripts when comparing the three libraries, suggesting regulation in response to these pathogens. A total of 191 microsatellites, with 104 having sufficient flanking sequences for primer design, and 1158 putative SNPs were identified from these EST resources in turbot.
A collection of 9256 high-quality ESTs was generated representing 3482 unique turbot sequences. A large proportion of defence/immune-related genes were identified, many of them regulated in response to specific pathogens. Putative microsatellites and SNPs were identified. These genome resources constitute the basis to develop a microarray for functional genomics studies and marker validation for genetic linkage and QTL analysis in turbot.
The turbot (Scophthalmus maximus; Scophthalmidae; Pleuronectiformes) is a commercially valuable flatfish that has been intensively cultured for the last decade. Its production has steadily increased to 7120 tonnes in 2006 (80% from Spain; FEAP, 2006) and represents one of the most promising aquaculture species in Europe. However, disease outbreaks in turbot have occurred frequently and losses due to infections constitute a serious problem for its culture [1, 2]. The use of antibiotics, vaccines and fish health management practices has partially solved the problem, but the achievement of large-scale production in the highly competitive world market requires enhancing resistance of cultured fish to diseases. Information on the immune response of turbot is still limited, and little is known about host-pathogen interactions in fish species. The screening and identification of immune-relevant genes is essential to analyze the genetic basis for infection, immunity and resistance to pathogens of economic relevance in aquaculture. Expressed sequence tag (EST) analysis is a powerful approach to provide a rapid and efficient method to go from expressed sequences to genes. ESTs are essential for studies of gene function [3, 4], but are also useful to identify polymorphic gene markers, such as microsatellites and single nucleotide polymorphisms (SNPs) [5–8]. These markers are the basis for genetic and physical mapping, and for comparative genome analysis [9–11]. From a practical perspective, maps can be applied for assisted selection programmes (MAS) and eventually for identification of genes related with quantitative traits (QTL) [12, 13]. In addition, ESTs constitute the basic resources to develop microarrays for functional genomics studies .
EST sequence resources are rapidly growing in molecular databases. However, the number of ESTs in fish is generally scarce, excluding some model species and Atlantic salmon among cultured fish [15–19]. The scarcity of EST resources in cultured fish limits the use of modern functional genomic approaches for selective breeding purposes . Among flatfish, aquaculture production has been successfully achieved in turbot, Japanese flounder and Atlantic halibut. Compared to the very large efforts for the development of EST resources in Japanese flounder (8856 ESTs) and Atlantic halibut (17659 ESTs) [21–28], EST resources in turbot are scarce (3171 ESTs). Less than 800 sequences have been deposited to date in the NCBI nucleotide database, most of them from anonymous microsatellite searching .
With the aim of increasing the genomic resources in turbot and identifying relevant genes for immunity, three cDNA libraries were constructed from mRNA isolated from immunity-related tissues of turbot (liver, spleen and head kidney) at different times after infection with Aeromonas salmonicida and the scuticociliate parasite Philasterides dicentrarchi. These pathogens are responsible for important disease outbreaks in turbot, as well as in other culture fish species [30–33]. Our main goal was to obtain the most accurate information possible to address functional genomic studies on disease resistance. However, the use of non-normalized cDNA libraries made it feasible to get a preliminary picture of the turbot genetic response to pathogens through analyzing transcript distribution among infected vs. control libraries. A total of 12584 ESTs were sequenced and compared to GenBank database and a large array of defence or immune-related genes was identified. Also, this large scale EST study increased the number of putative markers for mapping. A total of 191 microsatellites, of which 104 exhibited sufficient flanking sequences for primer design, and 2197 good quality SNPs were identified for the first time in turbot. The cDNA sequences generated will serve as a basis for microarray construction. This first EST study in turbot will provide the support for further research into the genetics, genomics and even proteomics of this important aquaculture species.
Results and Discussion
cDNA libraries and ESTs
EST analysis is an efficient and fast method for gene discovery [15, 17]. In Pleuronectiformes, this approach has been recently applied in Japanese flounder (Paralichthys olivaceus; [21, 22, 24, 25]), winter flounder (Pseudopleuronectes americanus; ), flounder (Platichthys flesus; [35, 36]) and Atlantic halibut (Hippoglossus hippoglossus; [26–28]. However, this order comprises around 600 species, many of them of great economic value both for fisheries and farming. In this study, we have addressed the construction of an EST database for the identification of genes related with immunity and defence in turbot.
Summary statistics of ESTs from turbot libraries
Unique sequences with no BLAST hits
Unique sequences with BLAST hits
Unique sequences with functional annotation
ESTs were identified by BLAST searches against nucleotide database at GenBank. Due to the low representation of fish genes, a protein-based homology strategy was also used in the international database searches. Protein sequences have been demonstrated to be more suitable to detect homology over long periods of evolutionary time . Our EST database pprovided a graphical view of all contigs, their PROSITE/PRINTS protein patterns and a search interface by keywords, microsatellite, gene, and UniGene/GO/KEGG information. Tools to search for sequences and markers based on annotations, to perform local BLAST searches, and to select sequences for a prospective microarray were also included. These tools used RepeatMasker http://www.repeatmasker.org/ for masking low-complexity sequences and OligoArray 2.1 http://berry.engin.umich.edu/oligoarray2_1/ to predict secondary structures and potential cross-hybridization.
As shown in Table 1, 1716 (49.3%) unique sequences displayed no significant similarity to known sequences or ESTs in the public databases, whereas the remaining 1766 (50.7%) showed significant matches with e-values ≤ 1e-5. Among these, BLAST database searches allowed assignment of putative function to 816 sequences (23.4%). In spite of their lower frequency among unique sequences, contigs were annotated more frequently (59.9%) than singletons (40.1%). As in other ESTs fish studies, the lower percentage of annotated singletons suggests that these are either novel fish-specific or rapidly evolving genes . Also, it is possible that bioinformatic errors could have a greater impact on singletons, since they are unique sequences whose information cannot be contrasted with other sequences in the database. The lower annotation success regarding similar genomic projects [28, 39] was probably related to the read length (around 500 bp) and specially the direction (from the 3'end) of sequencing. The 3' untranslated region (UTR) is approximately double the length of the 5'UTR according to GenBank fish entries. So, what is gained in specificity for microarray oligo-design is lost for gene annotation.
All unique sequences were annotated based on similarity using BLASTX or BLASTN  in the public databases GenBank NR and Unigene. The multiple annotations provided greater assurance about gene description and frequency of annotation than in a single database. The use of consensus sequences allowed sequences without significant similarity regions with a known protein (e.g., 5' or 3' noncoding regions) to be annotated if they were members of an annotated contig. All hits with e-value ≤ 1e-5 and their associated alignments were stored in the database and tracked with any associated functional annotation. We also ran AutoFACT  on all sequences in the database. It is interesting to note that while AutoFACT was able to come up with more function-specific information than our online tool (provided external database hits were found), it was not able to annotate as many sequences as our online custom tool.
Defence and immune-related annotated ESTs from turbot libraries
No. unique sequences
Major histocompatibility complex
Like in mammals, the immune system of fish is composed of non-specific and specific defence. The innate immune response is an important and highly developed defence mechanism against pathogens in fish . Examples of innate immunity include anatomic barriers, mechanical removal of pathogens, bacterial antagonism, pattern-recognition receptors, antigen-nonspecific defence compounds, complement pathways, phagocytosis, and inflammation . In Table 2, a list of defence/immune-related genes that were found three or more times in our libraries is shown. Complement related genes were predominant (7.9%) followed by apoptosis (4.9%) and immunoglobulin (Ig)-related (3.9%) genes. Glutathione S-transferase, heat shock proteins and cytochrome P450, elastases, major histocompatibility complex (MHC) and coagulation factors involved in innate immunity were also present, as well as others like interferon, perforin, hepcidin, nephrosin or alpha-2-macroglobulin. Interestingly, the majority of our immune-related cDNAs (75.0%) were reported for the first time in turbot, even though some of them have important roles in the immune response like B-cell linker, chemotaxin, complement components, IgD, IgM, interferon stimulated gene 12 (b2), lipopolysaccharide (LPS)-binding protein, natural killer (NK)-lysine type 1, peptidoglycan recognition protein, skin mucus lectin and tumour necrosis factor (TNF) receptor associated factor2, among others. Previously, only a hepcidin  and a natural killer cell enhancing factor  had been characterized in turbot.
The availability of a large number of sequences from immune-related cDNA libraries, both from non-infected fish (controls) as well as from fish challenged with specific pathogens, suggested a comparison of transcript profiles to identify genes regulated in response to these pathogens. Four or more sequences per gene (contig) are the minimum necessary to get statistical power to check the null hypothesis of even distribution of sequences among the three libraries. Ninety six genes presented four or more sequences in our libraries. However, taking into account the large number of tests performed, we decided to use a more conservative set of genes (with 6 or more sequences) to avoid type I errors. This new set comprised 72 genes, whose distribution among libraries is shown in Additional file 1 together with their probability of departure from the null hypothesis using a chi-square test. To use all information available and to get higher statistical power, we decided to perform an additional analysis by comparing the profiles of all genes grouped according to GO categories (see Additional file 2). Though this approach could be mixing genes with opposite expression patterns, it could provide significant trends associated to specific functions for complementing the "individual-level gene" analysis. Bonferroni correction for multiple tests was considered in this case, because of the higher statistical power achieved with this approach.
As expected, defence/immune-related genes appeared overrepresented in the set of 72 genes (37.5%) with regard to the total number of annotated genes in our libraries (24.9%). Some genes apparently responded in a similar way to both pathogens, being down- (myosin, nephrosin and several peptidases) or up-regulated (actin and lysozyme) regarding control. However, most genes responded to only one pathogen. The different infection profile of the pathogens used in this study, a bacterium (A. salmonicida) and a parasite (P. dicentrarchi), is expected to stimulate a specific set of genes in the host. So, protein biosynthesis was up-regulated in response to Aeromonas salmonicida, as suggested the increased amount or ribosome-related and elongation factor proteins, and their correspondent GO terms (see Additional file 2). Also, a general down-regulation of genes involved in proteolytic activity in response to this bacterium (chymotrypsin B precursor, trypsinogen1, chymotrypsinogen 2, trypsynogen-like serin protease and elastase precursor) was observed, in turn correlated to the lower presence of peptidase and proteolysis GO terms in this library. The possibility that A. salmonicida could be blocking these genes cannot be ruled out, since it would facilitate the infection process. Additionally, it was noteworthy that several iron metabolism-related genes turned out to be up-regulated after infection with A. salmonicida (haptoglobin fragment 1, globin-related proteins, hepcidin precursor). The high relevance of iron for bacteria infection could explain this observation. The up-regulation of hepcidin after bacteria infection has been recently described in turbot and gilthead seabream [44, 46]. Finally, some other immune relevant genes like MHC II alpha and beta antigens, complement C9, thrombin, chemotaxin, bactericidal permeability-increasing protein, apolipoprotein A-IV3, heat shock protein 90 beta and chemotaxin showed significant up-regulation in response to A. salmonicida.
However, the strongest genetic signature was observed in response to P. dicentrarchi infection. Specific genes were exclusively detected in this library and sometimes at high frequencies, like those homologous to Oryzias latipes and Paralichthys olivaceous cDNA clones and to mitochondrial ATP synthase alpha-subunit gene. Functional annotation of these unknown genes now appears relevant to understand the response of turbot to this pathogen. The same was reflected when GO categories were analyzed (see Additional file 2). A large number of differentially regulated gene categories (P = 0) appeared associated to P. dicentrarchi library according to the different GO criteria. Cytoskeleton organization and biogenesis, as well as carbohydrate metabolic process appeared up-regulated categories among biological processes. Transporter and hydrolase activities, as well as protein, ion and heme binding also were overexpressed regarding molecular function classification. Finally, cytoskeleton and cytosol appeared as the highest active cellular components. Looking at specific categories, a significant increase in the expression of antioxidant genes, like glutathione peroxidase and glutathione-S-transferase, was observed. These genes play a pivotal role to prevent cellular damage due to the increased reactive oxygen species (ROS) during infection . Paramá et al.  have recently demonstrated the increase in intracellular ROS by proteases of turbot kidney in response to this pathogen. Tissue trauma or invasion by pathogens induces changes in the quantities of several macromolecules in animal body fluids, which comprise one aspect of the acute phase response (APR). In fish, APR proteins include pentraxins, serum amyloid P, several components of the complement system, transferrin and thrombin. Up-regulation of pentraxin and serotransferrin was observed in response to P. dicentrarchi in our study. Transferrin up-regulation could be related to the increase in inflammation and enhanced oxidative stress typical of infections with this parasite . Finally, up-regulation of other important components of the immune response such the profilin and lysozyme was observed. Profilin-like protein has been reported as a toll-like receptor (TLR) 11 ligand in some parasites . TLRs are evolutionary conserved transmembrane proteins that recognize a unique pattern of molecules derived from pathogens or damaged cells, triggering robust but defined innate immune responses .
EST studies can also provide resources for the identification of polymorphic DNA markers such as microsatellites and SNPs. In our study, screening of EST sequences for short tandem repeats (2–6 bp) identified 191 microsatellites using a conservative criterion (≥ 6 and ≥ 8 repeats for tri/tetra/penta and dinucleotide motifs, respectively). Of these, 120 had significant hits in BLAST with e-value cut off ≤ 1e-5 and 71 were annotated. Most microsatellites were dinucleotide (128) and trinucleotide (56), while only 6 tetranucleotide and 1 pentanucleotide were found. Among these 191 ESTs sequences, 104 contained sufficient flanking sequences length for primer design. Fifty of these 104 microsatellites-containing ESTs were contigs, therefore the in silico comparison of the sequences included in these contigs allowed us to identify 11 putative polymorphic microsatellites.
Summary statistics of SNP identification from turbot EST resources
True SNPs andSNPs
Total sequences analysed
Number of contigs
Total SNPs detected
Total number of transitions
Total number of transversions
Total number of indels
Real and true quality SNP distribution in contigs with 4 or more ESTs
Number of contigs
with 4 sequences
with 5–10 sequences
with 11–20 sequences
with 21–30 sequences
with 31–50 sequences
with > 50 sequences
To our knowledge, this is the first report on a large transcriptional analysis in turbot providing new genomic resources in this important European aquaculture species. This study describes a collection of 9256 ESTs representing 3482 unique sequences obtained from three directionally cloned cDNA libraries from S. maximus, all of these being novel ESTs for this species. Therefore, this is a valuable EST collection, which increases genomic resources of turbot and enhances the genomic tools available for non-model fish species. The transcript profile comparison among the three libraries allowed the identification of putative genes generally or specifically related with infections in turbot. In addition, a high number of putative microsatellite and SNP EST-markers are now available for turbot map and highly useful for comparative mapping. These ESTs will be the basis for the development of a turbot microarray, focused on the characterization of the transcriptional response to pathogen exposure.
Tissue source and challenge
Two batches of 40 individuals (20–30 g each) obtained from a mixture of heterogeneous genetic families were collected at a specialized turbot fish farm. Fish from each batch were challenged intraperitoneally with A. salmonicida subsp. salmonicida and P. dicentrarchi, respectively [53, 54]. The dose was adjusted to obtain around 50% survival (LD50). The challenges were performed at the CETGA facilities in quarantine tanks. Fish were sacrificed prior to organ extraction using a lethal dose of MS222 anesthetic. In order to obtain mRNA representative of both innate and adaptive immune systems across the infection process, liver, spleen and head kidney tissues were collected from 5 sacrificed fish at five sampling points along the infective process for A. salmonicida (12 h, 1 day, 3 days, 7 days and 21 days post-inoculation) and at four sampling points for P. dicentrarchi (1 day, 3 days, 7 days and 15 days post-inoculation). Analogous batches of control fish were injected with saline serum and sampled at the same days of challenged fish. For each sample time, equal amounts of tissue from liver, kidney and spleen of each fish were pooled and immediately frozen in liquid nitrogen, constituting 15 pools for A. salmonicida (3 tissues × 5 sampling points), 12 for P. dicentrarchi (3 tissues × 4 sampling points) and the same pools for their respective controls. Each sample of pooled tissues was ground to a fine powder in a mortar and pestled with liquid nitrogen and stored at -80°C until being used for RNA extraction. The use of pools of individuals at each sampling point was planned, both for identifying putative SNPs, as well as for averaging individual effects in future microarray analysis on turbot immune response.
RNA isolation and cDNA library construction
Total RNA was extracted from pooled tissues of control and infected fish using TRIZOL Reagent (Life Technologies) according to manufacturer's recommendations. RNA quality was assessed in a Bioanalyzer (Bonsai Technologies). RNA was quantified using gel and NanoDrop® ND-1000 spectrophotometer (NanoDrop® Technologies Inc) estimations. Poly-A mRNA was isolated using the Dynabeads® mRNA Purification Kit (INVITROGEN). Three cDNA libraries (from A. salmonicida and P. dicentrarchi infected fish and control) were directionally constructed (5' EcoRI, 3' XhoI), with equal amounts of RNA from each tissue at each sampling time, using the ZAP-cDNA Library Construction Kit (STRATAGENE) following manufacturer's instructions except for size fractioning. This was performed on cDNAs prior to ligation into vector and carried out with the SizeSep 400 Spun Columns (GE HEALTHCARE). To allow characterization of the inserts in a plasmid system the library was mass excised according to the manufacturer's recommendations, converting it from plaque forming units to a phagemid in E. coli SOLR strain. Manually isolated colonies were randomly picked and arrayed into 96-well microtiter plates containing liquid selective medium. Glycerol stocks of overnight cultures were prepared in 96 well plates and stored at -80°C.
Sequencing and bioinformatics
Plasmid DNA was isolated from around 4000 clones from each library using the DirectPrep® 96 Miniprep kit (QIAGEN) and the Plasmid Miniprep96 kit (MILLIPORE) with a robotic platform (BIOMEK 3000), and checked by electrophoresis in 1% agarose gels. One part of the purified DNA was sequenced following the ABI Prism BigDye™ Terminator v3.1 Cycle Sequencing Kit protocol on an ABI 3100 DNA sequencer (Applied Biosystems) and the other with the DTCS kit protocol on a CEQ2000 DNA sequencer (BECKMAN COULTER). All clones were sequenced from their 3' ends using a standard T7 primer to obtain the highest specific sequences of genes for oligo-microarray design. Those clones which suffered a systematic drop on sequencing signal after poly-A tails were sequenced from 5' end.
A bioinformatic tool was developed in order to process all data. Basecalling from chromatogram traces was performed by using PHRED [55, 56]. Vector, poly-A tails and low-quality regions were trimmed from EST sequences using a custom Perl script, a local BLAST search engine and the trimmest utility from the EMBOSS suite. High quality ESTs (at least 100 bp and PHRED score ≥ 20 after removal of vector sequence, adapter, and poly-A tail) of both normal and infected libraries were combined and assembled to form clusters using CAP3 http://seq.cs.iastate.edu./.  with the overlapping identity percentage and minimum overlapping length parameters set to > 85.0% and 50 bp, respectively, in order to obtain highly reliable contig sequences. These contigs were manually revised to detect possible errors all along the bioinformatics process. ESTs that did not form contigs (singletons) and the contigs resultant of assemblage of multiple sequences were referred to as unique sequences. Singletons and consensus sequences of each contig (unique sequences) were compared against public databases. Unique sequences were searched by both Blastn and Blastx, and the corresponding outputs subsequently parsed using BioPerl for fish-relevant hits and significant UniGene information. GO, KEGG and COG terms were extracted from AutoFACT output, which in turn feeds from Tblastx and RPS-blast output. All results were stored in a mySQL database to be consulted and searched using a custom-designed, friendly AJAX web interface. The e-value cut off was ≤ 1e-5. All annotations thus obtained were complemented with AutoFACT output, whenever available.
Bioinformatic mining of microsatellites and SNPs
The set of 3482 unique sequences was searched for microsatellites using the program SPUTNIK http://espressosoftware.com/pages/sputnik.jsp. The minimum repeat number used for this search was 8 for dinucleotide and 6 for tri-, tetra-and pentanucleotide microsatellites. Microsatellite-containing ESTs were identified as candidates for marker development if they presented enough flanking sequences on either side of the repeats for primer design. Detection of single nucleotide polymorphisms (SNPs) was done using the pipeline QualitySNP with default settings . This program uses three filters for the identification of reliable SNPs: filter 1 screens for all potential SNPs; filter 2 uses a haplotype-based strategy to detect reliable SNPs; and filter 3 screens SNPs by calculating a confidence score based on sequence redundancy and quality. SNPs that pass filters 1 and 2 are called real SNPs, and those that pass all of the three filters are called true SNPs.
This study was supported by a Consellería de Pesca e Asuntos Marítimos and the Dirección Xeral de I+D – Xunta de Galicia project (2004/CP480). Authors wish to thank to Lucía Insua, Sonia Gómez, Ana Teijido and Susana Sánchez for technical assistance.
- Toranzo AE, Novoa B, Romalde JL, Núñez S, Devesa S, Mariño E, Silva R, Martínez E, Figueras A, Barja JL: Microflora associated with healthy and diseased turbot (Scophthalmus maximus) from three farms in northwest Spain. Aquaculture. 1993, 114: 189-202. 10.1016/0044-8486(93)90295-A.View ArticleGoogle Scholar
- Iglesias R, Paramá A, Álvarez MF, Leiro J, Ubeira FM, Sanmartín ML: Philasterides dicentrarchi (Ciliophora: Scuticociliatida) expresses surface immobilization antigens that probably induce protective immune responses in turbot. Parasitology. 2003, 126: 125-134. 10.1017/S0031182002002688.PubMedView ArticleGoogle Scholar
- Adams MD, Kelley JM, Gocayne JD, Dubnick M, Polymeropoulos MH, Xiao H, Merril CR, Wu A, Olde B, Moreno RF: Complementary DNA sequencing: expressed sequence tags and human genome project. Science. 1991, 252: 1651-1656. 10.1126/science.2047873.PubMedView ArticleGoogle Scholar
- Marra MA, Hillier L, Waterston RH: Expressed sequence tags-ESTablishing bridges between genomes. Trends Genet. 1998, 14: 4-7. 10.1016/S0168-9525(97)01355-3.PubMedView ArticleGoogle Scholar
- Liu ZJ, Karsi A, Dunham RA: Development of polymorphic EST markers suitable for genetic linkage mapping of catfish. Mar Biotechnol. 1999, 1: 437-447. 10.1007/PL00011800.PubMedView ArticleGoogle Scholar
- Liu ZJ, Li P, Kocabas A, Ju Z, Karsi A, Cao D, Patterson A: Microsatellite-containing genes from the channel catfish brain: evidence of trinucleotide repeat expansion in the coding region of nucleotide excision repair gene RAD23B. Biochem Biophys Res Commun. 2001, 289: 317-324. 10.1006/bbrc.2001.5978.PubMedView ArticleGoogle Scholar
- He C, Chen L, Simmons M, Li P, Kim S, Liu ZJ: Putative SNP discovery in interspecific hybrids of catfish by comparative EST analysis. Anim Genet. 2003, 34: 445-448. 10.1046/j.0268-9146.2003.01054.x.PubMedView ArticleGoogle Scholar
- Serapion J, Kucuktas H, Feng J, Liu Z: Bioinformatic mining of type I microsatellites from expressed sequence tags of channel catfish (Ictalurus punctatus). Mar Biotechnol. 2004, 6: 364-377. 10.1007/s10126-003-0039-z.PubMedView ArticleGoogle Scholar
- Murakawa K, Matsubara K, Fukushima A, Yoshii J, Okubo K: Chromosomal assignments of 3Vdirected partial cDNA sequences representing novel genes expressed in granulocytoid cells. Genomics. 1994, 23: 379-389. 10.1006/geno.1994.1514.PubMedView ArticleGoogle Scholar
- Bouza C, Hermida M, Pardo BG, Fernández C, Castro J, Fortes G, Sánchez L, Presa P, Pérez M, Sanjuán A, Comesaña S, Álvarez JA, Calaza M, Cal R, Piferrer F, Martínez P: A microsatellite genetic map of the turbot (Scophthalmus maximus). Genetics. 2007, 177: 2457-2467. 10.1534/genetics.107.075416.PubMed CentralPubMedView ArticleGoogle Scholar
- Moen T, Hayes B, Baranski M, Berg PR, Kjøglum S, Koop BF, Davidson WS, Omholt SW, Lien S: A linkage map of the Atlantic salmon (Salmo salar) based on EST-derived SNP markers. BMC Genomics. 2008, 9: 223-10.1186/1471-2164-9-223.PubMed CentralPubMedView ArticleGoogle Scholar
- Poompuang S, Hallerman EM: Toward detection of quantitative trait loci and marker-assisted selection in fish. Reviews Fish Sci. 1997, 5: 253-277.View ArticleGoogle Scholar
- Liu ZJ, Cordes JF: DNA marker technologies and their applications in aquaculture genetics. Aquaculture. 2004, 238: 1-37. 10.1016/j.aquaculture.2004.05.027.View ArticleGoogle Scholar
- Jørgensen SM, Afanasyev S, Krasnov A: Gene expression analyses in Atlantic salmon challenged with infectious salmon anemia virus reveal differences between individuals with early, intermediate and late mortality. BMC Genomics. 2008, 9: 179-10.1186/1471-2164-9-179.PubMed CentralPubMedView ArticleGoogle Scholar
- Hirono I, Aoki T: Expressed sequence tags of medaka (Oryzias latipes) liver mRNA. Mol Mar Biol Biotechnol. 1997, 6: 345-350.PubMedGoogle Scholar
- Lo J, Lee S, Xu M, Liu F, Ruan H, Eun A, He Y, Ma W, Wang W, Wen Z, Peng J: 15,000 Unique zebrafish EST clusters and their future use in microarray for profiling gene expression patterns during embryogenesis. Genome Res. 2003, 13: 455-466. 10.1101/gr.885403.PubMed CentralPubMedView ArticleGoogle Scholar
- Clark MS, Edwards YJ, Peterson D, Clifton SW, Thompson AJ, Sasaki M, Suzuki Y, Kikuchi K, Watabe S, Kawakami K, Sugano S, Elgar G, Johnson SL: Fugu ESTs: new resources for transcription analysis and genome annotation. Genome Res. 2003, 13 (12): 2747-2753. 10.1101/gr.1691503.PubMed CentralPubMedView ArticleGoogle Scholar
- Rise ML, von Schalburg KR, Brown GD, Mawer MA, Devlin RH, Kuipers N, Busby M, Beetz-Sargent M, Alberto R, Gibbs AR, Hunt P, Shukin R, Zeznik JA, Nelson C, Jones SR, Smailus DE, Jones SJ, Schein JE, Marra MA, Butterfield YS, Stott JM, Ng SH, Davidson WS, Koop BF: Development and application of a salmonid EST database and cDNA microarray: data mining and interspecific hybridization characteristics. Genome Res. 2004, 14: 478-490. 10.1101/gr.1687304.PubMed CentralPubMedView ArticleGoogle Scholar
- Jaillon O, Aury JM, Brunet F, Petit JL, Stange-Thomann N, Mauceli E, Bouneau L, Fischer C, Ozouf-Costaz C, Bernot A, Nicaud S, Jaffe D, Fisher S, Lutfalla G, Dossat C, Segurens B, Dasilva C, Salanoub M, Levy M, Boudet N, Castellano S, Anthouard V, Jubin C, Castelli V, Katinka M, Vacherie B, Biemont C, Skalli Z, Cattolico L, Poulain J, De Berardinis V, Cruaud C, Duprat S, Brottier P, Coutanceau JP, Gouzy J, Parra G, Lardier G, Chapple C, McKernan KJ, McEwan P, Bosak S, Kellis M, Volff JN, Guigo R, Zody MC, Mesirov J, Lindblad-Toh K, Birren B, Nusbaum C, Kahn D, Robinson-Rechav M, Laudet V, Schachter V, Quetier F, Saurin W, Scarpelli C, Wincker P, Lander ES, Weissenbach J, Roest Crollius H: Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype. Nature. 2004, 431: 946-957. 10.1038/nature03025.PubMedView ArticleGoogle Scholar
- Chen SL, Xu MY, Hu SN, Li L: Analysis of immune relevant genes expressed in red sea bream (Chrysophrys major) spleen. Aquaculture. 2004, 240: 115-130. 10.1016/j.aquaculture.2004.07.008.View ArticleGoogle Scholar
- Inoue S, Nam BH, Hirono I, Aoki T: A survey of expressed genes in Japanese flounder (Paralichthys olivaceus) liver and spleen. Mol Mar Biol Biotechnol. 1997, 6: 376-380.PubMedGoogle Scholar
- Nam BH, Yamamoto E, Hirono I, Aoki T: A survey of expressed genes in the leukocytes of Japanese flounder, Paralichthys olivaceus, infected with Hirame rhabdovirus. Dev Comp Immunol. 2000, 24: 13-24. 10.1016/S0145-305X(99)00058-0.PubMedView ArticleGoogle Scholar
- Aoki T, Hirono I, Kim M-G, Katagiri T, Tokuda Y, Toyohara H, Yamamota E: Identification of viral induced genes in Ig+ leucocytes of Japanese flounder Paralichthys olivaceus, by differential hybridisation with subtracted and un-subtracted cDNA probes. Fish Shellfish Immunol. 2000, 10: 623-630. 10.1006/fsim.2000.0279.PubMedView ArticleGoogle Scholar
- Arma NR, Hirono I, Aoki T: Genes expressed in Japanese flounder Paralichthys olivaceus spleen: analysis of genes involved in immune function. Fisheries Sci. 2005, 71: 1304-1323. 10.1111/j.1444-2906.2005.01097.x.View ArticleGoogle Scholar
- Aoki T, Hirono I: Immune relevant genes of Japanese flounder, Paralichthys olivaceus. Comp Biochem Phys D. 2006, 1: 115-121.Google Scholar
- Park KC, Osborne JA, Tsoi SCM, Brown LL, Johnson SC: Expressed sequence tags analysis of Atlantic halibut (Hippoglossus hippoglossus) liver, kidney and spleen tissues following vaccination against Vibrio anguillarun and Aeromonas salmonicida. Fish Shellfish Immunol. 2005, 18: 393-415. 10.1016/j.fsi.2004.10.003.PubMedView ArticleGoogle Scholar
- Galloway TF, Bardal T, Kvam SN, Dahle SW, Nesse G, Randol M, Kjorsvik E, Andersen O: Somite formation and expression of MyoD, myogenin and myosin in Atlantic halibut (Hippoglossus hippoglossus L.) embryos incubated at different temperatures: transient asymmetric expression of MyoD. J Exp Biol. 2006, 209: 2432-2441. 10.1242/jeb.02269.PubMedView ArticleGoogle Scholar
- Douglas SE, Knickle LC, Kimball J, Reith ME: Comprehensive EST analysis of Atlantic halibut (Hippoglossus hippoglossus), a commercially relevant aquaculture species. BMC Genomics. 2007, 8: 144-10.1186/1471-2164-8-144.PubMed CentralPubMedView ArticleGoogle Scholar
- Pardo BG, Fernández C, Hermida M, Vázquez A, Pérez M, Presa P, Calaza M, Alvarez-Dios JA, Comesaña AS, Raposo-Guillán J, Bouza C, Martínez P: Development and characterization of 248 novel microsatellite markers in turbot (Scophthalmus maximus). Genome. 2007, 50: 329-332. 10.1139/G06-154.PubMedView ArticleGoogle Scholar
- Toranzo AE, Barja JL: First report of furunculosis in turbot reared in floating cages in North-West of Spain. Bull Eur Ass Fish Pathol. 1992, 12: 147-149.Google Scholar
- Björnsdóttir B, Gudmundsdóttir S, Bambir SH, Gudmundsdóttir BK: Experimental infection of turbot, Scophthalmus maximus (L.), by Aeromonas salmonicida subsp. achromogenes and evaluation of cross protection induced by a furunculosis vaccine. J Fish Dis. 2005, 28: 181-188. 10.1111/j.1365-2761.2005.00617.x.PubMedView ArticleGoogle Scholar
- Iglesias R, Paramá A, Alvarez MF, Leiro J, Fernandez J, Sanmartin ML: Philasterides dicentrarchi (Ciliophora, Scuticociliatida) as the causative agent of scuticociliatosis in farmed turbot Scophthalmus maximus in Galicia (NW Spain). Dis Aquat Organ. 2001, 46: 47-55. 10.3354/dao046047.PubMedView ArticleGoogle Scholar
- Paramá A, Arranz JA, Alvarez MF, Sanmartin ML, Leiro J: Ultrastructure and phylogeny of Philasterides dicentrarchi (Ciliophora, Scuticociliatia) from farmed turbot in NW Spain. Parasitology. 2006, 132: 555-564. 10.1017/S0031182005009534.PubMedView ArticleGoogle Scholar
- Douglas SE, Gallant JW, Bullerwell CE, Wolff C, Munholland J, Reith ME: Winter flounder expressed sequence tags: establishment of an EST database and identification of novel fish genes. Mar Biotechnol. 1999, 1: 458-464. 10.1007/PL00011802.PubMedView ArticleGoogle Scholar
- Williams TD, Gensberg K, Minchin SD, Chipman JK: A DNA expression array to detect toxic stress response in European flounder (Platichthys flesus). Aquat Toxicol. 2003, 65: 141-157. 10.1016/S0166-445X(03)00119-X.PubMedView ArticleGoogle Scholar
- Sheader DL, Gensberg K, Lyons BP, Chipman K: Isolation of differentially expressed genes from contaminant exposed European flounder by suppressive, subtractive hybridisation. Mar Environ Res. 2004, 58: 553-557. 10.1016/j.marenvres.2004.03.044.PubMedView ArticleGoogle Scholar
- Pearson WR: Identifying distantly related protein sequences. Comput Appl Biosci. 1997, 13: 325-332.PubMedGoogle Scholar
- Paschall JE, Oleksiak MF, VanWye JD, Roach JL, Whitehead JA, Wyckoff GJ, Kolell KJ, Crawford DL: FunnyBase: a systems level functional annotation of Fundulus ESTs for the analysis of gene expression. BMC Genomics. 2004, 5: 96-10.1186/1471-2164-5-96.PubMed CentralPubMedView 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 of ESTs from the eastern oyster, Crassostrea virginica Gmelin and identification of microsatellite and SNP markers. BMC Genomics. 2007, 8: 157-10.1186/1471-2164-8-157.PubMed CentralPubMedView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.PubMedView ArticleGoogle Scholar
- Koski LB, Gray MW, Lang BF, Burger G: AutoFACT: An Automatic Functional Annotation and Classification Tool. BMC Bioinformatics. 6: 151-10.1186/1471-2105-6-151.Google Scholar
- Gene Ontology: tool for the unification of biology. The Gene Ontology Consortium. Nature Genet. 2000, 25: 25-29. 10.1038/75556.View ArticleGoogle Scholar
- Ellis AE: Innate host defence mechanisms of fish against viruses and bacteria. Dev Comp Immunol. 2001, 25: 827-839. 10.1016/S0145-305X(01)00038-6.PubMedView ArticleGoogle Scholar
- Chen SL, Li W, Meng L, Sha ZX, Wang ZJ, Ren GC: Molecular cloning and expression analysis of a hepcidin antimicrobial peptide gene from turbot (Scophthalmus maximus). Fish Shellfish Immunol. 2007, 22: 172-181. 10.1016/j.fsi.2006.04.004.PubMedView ArticleGoogle Scholar
- Chen Y, Zhang YX, Fan TJ, Meng L, Ren GC, Chen SL: Molecular identification and expression analysis of the natural killer cell enhancing factor (NKEF) gene from turbot (Scophthalmus maximus). Aquaculture. 2006, 261: 1186-1193. 10.1016/j.aquaculture.2006.09.034.View ArticleGoogle Scholar
- Cuesta A, Meseguer J, Esteban MA: The antimicrobial peptide hepcidin exerts an important role in the innate immunity against bacteria in the bony fish gilthead seabream. Mol Immunol. 2008, 45: 2333-2342.PubMedView ArticleGoogle Scholar
- Valko M, Rhodes CJ, Moncol J, Izakovic M, Mazur M: Free radicals, metals and antioxidants in oxidative stress-induced cancer. Chem Biol Interact. 2006, 160: 1-40. 10.1016/j.cbi.2005.12.009.PubMedView ArticleGoogle Scholar
- Paramá A, Castro R, Arranz JA, Sanmartín ML, Lamas J, Leiro J: Scuticociliate cysteine proteinases modulate turbot leucocyte functions. Fish Shellfish Immunol. 2007, 23: 945-956. 10.1016/j.fsi.2007.02.001.PubMedView ArticleGoogle Scholar
- Ishii KJ, Uematsu S, Akira S: 'Toll' gates for future immunotherapy. Curr Pharm Des. 2006, 12: 4135-4142. 10.2174/138161206778743484.PubMedView ArticleGoogle Scholar
- Picoult-Newberg L, Ideker TE, Pohl MG, Taylor SL, Donaldson MA, Nickerson DA, Boyce-Jacino M: Mining SNPs from EST databases. Genome Res. 1999, 9: 167-174.PubMed CentralPubMedGoogle Scholar
- Syvanen AC: Accessing genetic variation: genotyping single nucleotide polymorphisms. Nat Rev Genet. 2001, 2: 930-942. 10.1038/35103535.PubMedView ArticleGoogle Scholar
- Tang J, Vosman B, Voorrips RE, Linden van der CG, Leunissen JAM: QualitySNP: a pipeline for detecting single nucleotide polymorphisms and insertions/deletions in EST data from diploid and polyploid species. BMC Bioinformatics. 2006, 7: 438-10.1186/1471-2105-7-438.PubMed CentralPubMedView ArticleGoogle Scholar
- Magariños B, Santos Y, Romalde JL, Barja JL, Toranzo AE: Pathogenic activities of live cells and extracellular products of the fish pathogen Pasteurella piscicida. J Gen Microbiol. 1992, 138: 2491-2498.PubMedView ArticleGoogle Scholar
- Paramá A, Iglesias R, Álvarez MF, Leiro J, Aja C, Sanmartín ML: Philasterides dicentrarchi (Ciliophora, Scuticociliatida): experimental infection and possible routes of entry in farmed turbot (Scophthalmus maximus). Aquaculture. 2003, 217: 73-80. 10.1016/S0044-8486(02)00523-9.View 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: 175-185.PubMedView ArticleGoogle Scholar
- Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8: 186-194.PubMedView ArticleGoogle Scholar
- Huang X, Madan A: CAP3: a DNA sequence assembly program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.PubMed CentralPubMedView ArticleGoogle 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.