- Research article
- Open Access
Comprehensive EST analysis of Atlantic halibut (Hippoglossus hippoglossus), a commercially relevant aquaculture species
BMC Genomics volume 8, Article number: 144 (2007)
An essential first step in the genomic characterisation of a new species, in this case Atlantic halibut (Hippoglossus hippoglossus), is the generation of EST information. This forms the basis for subsequent microarray design, SNP detection and the placement of novel markers on genetic linkage maps.
Normalised directional cDNA libraries were constructed from five different larval stages (hatching, mouth-opening, midway to metamorphosis, premetamorphosis, and post-metamorphosis) and eight different adult tissues (testis, ovary, liver, head kidney, spleen, skin, gill, and intestine). Recombination efficiency of the libraries ranged from 91–98% and insert size averaged 1.4 kb. Approximately 1000 clones were sequenced from the 5'-end of each library and after trimming, 12675 good sequences were obtained. Redundancy within each library was very low and assembly of the entire EST collection into contigs resulted in 7738 unique sequences of which 6722 (87%) had matches in Genbank. Removal of ESTs and contigs that originated from bacteria or food organisms resulted in a total of 7710 unique halibut sequences.
A Unigene collection of 7710 functionally annotated ESTs has been assembled from Atlantic halibut. These have been incorporated into a publicly available, searchable database and form the basis for an oligonucleotide microarray that can be used as a tool to study gene expression in this economically important aquacultured fish.
Atlantic halibut is a cold-water flatfish native to the North Atlantic that shows excellent potential for production in aquaculture due to its highly prized white flesh. Flatfish have long been a choice food fish, with many members of the group e.g., halibuts, flounders, soles, turbot, and plaice, having great commercial value especially in Asia. With the general worldwide decline in the wild fishery, and the predicted global collapse of all currently fished species by the year 2048 , it is crucial that alternatives such as aquaculture be pursued. Investigations into producing flatfish by aquaculture have been underway for the last fifteen to twenty years. Aquaculture production of Japanese flounder, turbot, Atlantic halibut and others has now been successfully achieved, although improvements in efficiency are still clearly required.
Production of Atlantic halibut is relatively recent and is currently underway in Norway, Iceland, Scotland and Canada. Significant hurdles must still be overcome, particularly with regard to judging when to spawn females, reproduction and sex determination , nutrition [3, 4], and enhancing disease resistance. The application of genomics technologies to thoroughly characterize the biological processes of reproduction, development, nutrition, and immunity promises to improve our knowledge of this poorly understood fish and provide for long-term enhancements in aquaculture production.
Flatfish, members of the order Pleuronectiformes, comprise a biologically interesting group of fish. During development, in a process known as metamorphosis, these fish reorient themselves to lie on one side, the body flattens, and the eye migrates to the other side of the body. This settling of the fish on the side vacated by the migrating eye requires a complex reorganization of skeletal, nervous and muscle tissues . Significant losses due to mortality in the early larval stages, as well as developmental abnormalities such as malpigmentation [6, 7], bone deformities , and incomplete eye migration , have hampered the successful production of flatfish. A better understanding of these processes at the molecular level and the impact that rearing conditions have on survival, metamorphosis, and growth will improve the commercial feasibility of flatfish aquaculture.
Expressed sequence tag (EST) surveys of species of interest provide a great deal of information-rich data on the expressed portion of an organism's genome  and are invaluable in comparative genomics . In addition, they can be an important source of microsatellites and single nucleotide polymorphisms (SNPs) that can be used for genetic mapping . EST surveys in teleosts have been performed from both non-normalised cDNA libraries and libraries generated by suppression subtractive hybridisation  in order to elucidate genes involved in immunity [13–28], muscle formation , endocrinology [30–32], and toxin production . The most highly represented teleost in dbEST is the zebrafish, Danio rerio, with over 1.3 million ESTs. Salmonids are also well-represented with over half a millions ESTs. Economically important fish species such as catfish, cod, Japanese flounder, sea bream and sea bass have increasingly been the subject of genomic studies and these species are now represented by several thousand or even tens of thousands of ESTs.
Among the flatfish, considerable effort has been made to determine ESTs for Japanese flounder [13, 17, 22]. There are currently 8842 ESTs for this species, although most are unannotated. Furthermore, the majority of these sequences arise from non-normalised cDNA libraries or those constructed from immune-stimulated fish; few ESTs have been sequenced from non-immune tissues or different developmental stages of flatfish. For turbot, the other main commercially relevant flatfish, there are no published EST studies and most of the sequences from this species in GenBank are microsatellites or rRNA.
As a first step towards developing genomics tools for Atlantic halibut, a large-scale EST survey was performed, annotations undertaken where possible, and a searchable database set up. Considerable effort has been made to annotate these ESTs and associate them with Gene Ontology (GO) terms to facilitate subsequent microarray analyses. The species-independent structured GO vocabulary  is widely accepted and used in most large scale genome annotation projects.
cDNA libraries and ESTs
Normalized cDNA libraries were constructed from eight Atlantic halibut tissues and five developmental stages. For all libraries, the percentage of clones with inserts was high, ranging from 91 to 98%, with average insert sizes between 1.3 and 1.6 kb (Table 1). After trimming and vector removal, 12675 good sequences were obtained. Average read length after trimming for each library was approximately 600–700 bp with the majority of reads compiled from all 13 libraries being between 700 and 800 bp (Figure 1). Clustering of the sequences using Paracel transcript assembler yielded 7738 unique sequences (redundancy factor of 1.6), demonstrating the excellent normalization achieved by this approach. Twenty-eight ESTs or EST clusters, mostly from later larval stages, had matches to food organisms such as Artemia and were removed from the database. The 7710 remaining unique ESTs consisted of 2556 unique clusters and 5154 singletons. The normalization of the libraries precludes an accurate analysis of transcript abundance, although sequences present multiple times in the normalized libraries presumably correspond to high abundance transcripts. Of the clones that were represented 5 times or more in a given library, the most common were a lectin-like protein, a 14 Kda apolipoprotein, tropomyosin, parvalbumin, and apolipoprotein A1 (Table 2). Parvalbumin, lectin-like protein and the 14 Kda apolipoprotein were each over-represented in two libraries.
Screening of EST sequences for short tandem repeats (2 – 5 bp) identified 129 that contained microsatellite sequences. Of these, 60 had 2 bp repeats (mostly GT or GA), 58 had 3 bp repeats, 7 had 4 bp repeats and 4 had 5 bp repeats. Sixty of these loci were polymorphic and were incorporated into our halibut linkage map (D. Reid, C. Smith, D. Martin-Robichaud, M. Reith, unpublished). All EST data have been deposited in GenBank (accession numbers EB029285–EB041700 &EB080851–EB080975), and preliminary annotations are available on the Pleurogene website .
All ESTs and clusters were initially annotated with AutoFACT, which carries out a series of Blast analyses to identify the EST and attempts to assign function and/or pathway information. AutoFACT provided product or gene names for 4786 sequences (62%) and 254 additional sequences had domain names associated with them (Table 3). Of the remaining sequences, 28 matched ribosomal RNAs, 824 matched unassigned proteins, 802 matched unknown ESTs, and 1016 had no database hits ("unclassified").
Of the 5040 functionally annotated and domain name-containing sequences, Gene Ontology (GO) classifications were found for 3878 using AutoFACT, Goblet and InterPro. The distribution of sequences among GO classifications are represented in Figures 2, 3, 4. Sequences with GO terms corresponding to cellular component fell into 27 categories (Figure 2), molecular function into 34 categories (Figure 3) and biological process into 46 categories (Figure 4). Clusters of Orthologous Groups (COG; ) categories were determined for 1191 sequences, and Kyoto Encyclopedia of Genes and Genomes (KEGG; ) categories were found for 578 sequences. Sequences assigned to COG fell into 22 categories and those assigned to KEGG fell into 185 pathways (Table 3).
We also performed a similar analysis on the 4110 partially annotated Atlantic halibut sequences deposited in GenBank by other groups. After assembly, a unigene set of 2337 sequences was obtained of which 40 were rRNA, 781 were unclassified, 80 matched unassigned proteins and 1436 were informative annotations. However, of the informative annotations, 531 fell into five clusters: nuclease diphosphate kinase B (73), cytochrome c oxidase subunit III (109), cytochrome c oxidase subunit II (182), cytochrome c oxidase subunit I (62) and cytochrome b (105). Similarly, of the 788 sequences that received COG annotations, over half (487) were associated with energy production and conversion, and were predominantly of mitochondrial origin.
When the complete EST set of 12675 sequences was searched for GO terms using the GOSLIM classification, two libraries were substantially more enriched in GOSLIM terms than others: the ovary library had 4496, and the testis library had 3371 hits to GOSLIM terms, respectively. Some of this enrichment can be explained by the increased number of ESTs sequenced from these libraries (50%), but even taking this into account, there are still many more GOSLIM terms in each of these libraries than the other tissue-specific libraries, which have between 1017 (gill) and 1871 (liver) terms. Of the larval libraries, that from the mouth-opening stage had 2415 hits to GOSLIM terms. Interestingly, the library constructed from larvae midway to metamorphosis had only 443 hits to GOSLIM terms. The remainder of the libraries had between 1444 and 1788 GOSLIM terms.
Functional categories of cDNAs
The functional categories included in the COG and GO classification schemes are not entirely congruent and thus not directly comparable. Broadly speaking, the COG classification gave a good overview of the cellular processes represented in the EST library. As can be seen from Table 4, approximately one third of the ESTs are associated with a COG category related to transport, metabolism and energy, 44% of the ESTs are associated with nucleic acid processes, and 4% of the ESTs are associated with cell structure. The remainder of the classified ESTs represent genes involved in signal transduction, defense and intracellular trafficking. KEGG annotations were found for 578 different sequences, the most common pathways represented are shown in Table 5.
To improve our understanding of flatfish biology and the problems associated with their development and rearing, a comparative genomics program focusing on Atlantic halibut and Senegal sole (Solea senegalensis) has been initiated (see ). As a prelude to construction of a DNA microarray, the EST survey reported here has been carried out.
Two previous EST surveys have been conducted in Atlantic halibut: one from a study of the effect of vaccination , resulting in approximately 1000 sequences, and a second from an investigation of muscle somite formation , resulting in approximately 4250 sequences. The study reported here greatly enriches the genomic resources for this commercially important flatfish by adding more than 12,000 ESTs to the partially annotated sequences that had already been deposited.
Over 5000 of the 7710 unique transcripts represented by the ESTs have been functionally annotated. These annotations and the development of a searchable database containing all of the information associated with each EST add enormous value to such a study. The main categories of genes represented in the ESTs are involved in binding, catalytic activity, transport, metabolism, response to stimuli, signal transduction, nucleic acid processes, and cellular biogenesis. Again, this adds substantially to the Atlantic halibut sequences from other research groups, many of which were of mitochondrial origin, and which contained slightly over 900 informative AutoFACT annotations.
With only modest resources for EST sequencing available, we choose to normalize our cDNA libraries so as to maximize the number of different ESTs sequenced. The normalization method used (Evrogen Trimmer kit) was very effective at reducing the number of highly expressed cDNAs. Since the libraries are well-normalised (redundancy factor of only 1.5), it is not possible to gain insights into the actual abundance of different types of transcripts; however, the enrichment in GO terms in the reproductive tissues and one of the larval libraries indicates that these libraries represent a broad diversity of transcripts, indicative of the high metabolic and proliferative characteristics of ovary and testis tissues. Larvae at the early mouth-opening stage of development are also undergoing tremendous metabolic changes as they transition from the yolk-sac to first-feeding stages. On the other hand, the library made from larvae midway to metamorphosis has very few GO terms associated with it, possibly because a large number of genes associated with the unusual metamorphic process in flatfish have not yet been described.
A number of ESTs were restricted to only a single tissue library and as such, may be good tissue-specific markers for in situ hybridisation and aid in tracking the appearance of different tissues during development [39, 40]. For example, ESTs for a renal organic ion transporter and nephrosin are only found in the head kidney library. Several transporters (for amino acids and various solutes), binding proteins (for lipid, sterol and lipoproteins) as well as digestive enzymes (elastases, peptidase, aminopeptidase N, carboxypeptidase B, triglyceride lipase) are only found in the intestine library. Various complement components, apolipoproteins, fatty acid binding protein, alpha-2-macroglobulin and biliverdin reductase A are only present in the liver library. A single EST unique to spleen was identified – metaxin 2, similar to the von Willebrand clotting factor. Unique to the ovary library were a number of ESTs specific to reproduction – zona pellucida protein, vitelline envelope protein, chorion protein, choriolytic enzyme, aquaporin, alveolin, estrogen receptor binding protein and luteinizing hormone beta. Similarly, the testis library uniquely contained ESTs specific to reproduction – spindlin protein C, testis intracellular mediator protein, a cysteine and glycine-rich protein, and periostin. The gill library uniquely contained two ESTs involved in chloride transport and the skin uniquely contained keratin, epithelial membrane protein-3, epiplakin and dermatopontin.
The liver, head kidney and spleen libraries proved to be an excellent source of genetic information concerning hematopoiesis and immune function in this fish. The head kidney is the major site of hematopoiesis in fish and an EST survey of zebrafish kidney revealed many insights into this process in fish . From our Atlantic halibut EST survey, many complement components, immune type receptors, lectins, defense proteins, MHC I and II components, cytokines, chemokines as well as signal transduction molecules and transcription factors involved in expression of immune genes were identified. It should be noted that ESTs for components of the immune system were also found in other tissue libraries, particularly those exposed to the environment such as skin and gill; these arose from circulating or resident immune cells in these tissues. This new sequence information will greatly enhance our understanding of the immune system of flatfish and provide molecular tools for further studying disease resistance.
The identification of microsatellite sequences in the Atlantic halibut ESTs will aid in the completion of a genetic linkage map of Atlantic halibut that is currently being constructed. Since these microsatellites are linked to genes, they are useful as Type I markers.
The addition of over 7700 ESTs, of which 5040 are functionally annotated, significantly enhances the genomic tools available for non-model fish species. Given the high degree of sequence similarity between flatfish species, the Atlantic halibut ESTs will be of great interest to the flatfish researchers in general, as well as the halibut aquaculture research community. The publicly accessible, searchable database also adds substantial value to the genomic data generated in this study. This EST survey has provided a number of microsatellite markers that have been placed on the Atlantic halibut genetic linkage map (Reith, pers. comm.) as well as probes for cellular localisation studies by in situ hybridisation. It has also laid the groundwork for the design and construction of a microarray for studying gene expression under different environmental conditions to better understand the impact of nutrition, stress, and environmental conditions on aquaculture production.
Fish rearing and sampling
Larvae were reared at Scotian Halibut Limited (Clarks Harbour, NS, Canada) in constant light (approximately 1000 lux at the surface) in 7 m3 tanks with flow-through salt water (32 ppt) maintained at 11 ± 0.2°C using a heat exchanger. Larvae were fed Artemia until weaning onto artificial feed at 65 DPH. The ages and sizes of the larvae at the different stages were as follows: hatching (1 dph; 10 mm), mouth-opening (21 dph;15 mm), midway to metamorphosis (64 dph; 20 mm), premetamorphosis (91 dph; 25 mm), and post-metamorphosis (104 dph; 30–35 mm).
Livers were sampled from one male (104 cm; 18 kg) and one female (130 cm; 31 kg) adult broodstock maintained at St. Andrews Biological Station, St. Andrews, NB, Canada. All other tissues were sampled from 3 male and 3 female immature fish (672–945 g) reared at the Institute for Marine Biosciences Marine Research Station (Sambro, NS, Canada). Immediately before sampling, fish were transferred to a bucket containing an overdose of TMS-Aqua MS-222 (Syndel, Vancouver, BC, Canada). Ovaries from 3 females were pooled, as were testes from 3 males. Gill, head kidney, intestine, skin, and spleen samples from each group of 6 fish were pooled, and all tissue samples preserved in RNALater (Ambion, Austin, TX, USA), and stored at -80°C until use. Larvae were pooled and preserved in RNALater, and stored at -80°C until use. All animal procedures were approved by the NRC Institute for Marine Biosciences Animal Care Committee.
cDNA library construction
For liver, mRNA was extracted from one female and one male halibut using the FastTrack kit (Invitrogen, Burlington, ON, Canada) and equal amounts of mRNA were combined. For spleen, total RNA was extracted from pooled tissues using Trizol Reagent (Invitrogen). For other tissues, mRNA was extracted from pooled tissues using the Micro-FastTrack kit (Invitrogen). For larval libraries, RNA was extracted from pooled samples (~20 larvae each for post-hatch and mouth-opening stages, 15 larvae for midway to metamorphosis, and 5 larvae each for pre- and post-metamorphosis) using Trizol Reagent. All RNA isolation kits were used according to the manufacturer's protocols.
First strand cDNA was prepared from 0.25–0.4 μg mRNA or 2 μg total RNA using the Creator SMART cDNA method (Clontech, Palo Alto, CA, USA) and PowerScript reverse transcriptase (Clontech). The CDS-3M adaptor, included in the TRIMMER-DIRECT kit (Evrogen, Moscow, Russia), was used instead of the SMART CDSIII primer. cDNA was amplified by LD-PCR according to the Creator SMART cDNA method (Clontech) using the 5' PCR primer as the forward and reverse primer. The optimal number of cycles to yield sufficient cDNA for normalisation, but remain within the exponential phase of amplification, was determined by analysing aliquots of the PCR reaction after every second cycle on agarose gels. In all cases, sufficient cDNA was obtained in 18 cycles or less, ensuring even the rarest messages were represented. Amplified cDNA was purified using the QIAquick PCR purification kit (Qiagen, Valencia, CA, USA), quantitated using a NanoDrop® ND-1000 spectrophotometer (NanoDrop® Technologies Inc, Wilmington, DE, USA) and normalized using the TRIMMER-DIRECT protocol (Evrogen). After digestion with SfiI, products smaller than 500 bp were removed using the Chroma Spin-400 column column as described in the Creator SMART protocol.
The resulting cDNAs were directionally cloned into the SfiI sites of pDNR-LIB (Clontech) and transformed into ElectroMAX DH10B T1 phage-resistant cells (Invitrogen) by electroporation using the Cell Porator and Voltage Booster system (Gibco BRL). The Cell Porator settings were 400 V, 330 μF capacitance, low Ω impedance and fast charge rate, and the Voltage Booster was set at 4 kΩ. For each library, 106 primary transformants were amplified by the semi-solid amplification method described in Stratagene's pBluescript XR cDNA library construction kit manual. Randomly picked clones (96 from each library) were screened for insert size by protoplasting  or by PCR using the M13 forward and reverse primers flanking the multiple cloning site of the vector.
Individual bacterial colonies were picked into 96- or 384-well plates containing LB/glycerol using the QPix colony picker (Genetix Ltd., New Milton, Hampshire, UK). A 96-well test plate was prepared from each library for sequencing and if the quality of the library was good, clones were sequenced from two additional 96-well and two 384-well plates, giving a total of 1056 reads. Two additional 384-well plates were sequenced from each of the ovary and testis libraries. Plates were incubated overnight at 37°C. The resulting bacterial suspensions were inoculated into lysis buffer and denatured at 95°C for 5 minutes. DNA from each clone was amplified using TempliPhi™ DNA polymerase (GE Healthcare, Baie d'Urfe, QC, Canada) according to manufacturer's instructions. DNA sequencing was performed using ET terminator chemistry (GE Healthcare) in the 5' direction (primer sequence GGCCGCATAACTTCGTATAGC). Reactions were processed using Sera-Mag™ magnetic carboxylate-modified microparticles (Seradyn™, Indianapolis, IN, USA) to remove excess fluorescent terminators before loading onto GE Healthcare MegaBACE 4000 capillary DNA sequencers. Clones from each library were replicated into glycerol stocks and stored at -80°C.
ESTs and annotation
ESTs were clustered using Paracel Transcript Assembler 3.0 (Paracel Inc., Pasadena, CA), which is based on the CAP4 clustering algorithm . Annotation was performed using AutoFACT  and the default parameters with UniProt's UniRef90, NCBI's nr, KEGG, COG, PFAM, LSU, SSU and, for contigs, est_others. AutoFACT summary results are stored in our database . Each sequence has one or more AutoFACT results associated with it. Each AutoFACT result is related to the most informative BLAST hit from each of the queried set of databases. Those hits are also stored in the Pleurogene database. GO annotations associated with either the contributing hits or the AutoFACT results are stored in a separate table.
Due to the low level of GO annotation obtained with AutoFACT (1640 sequences out of 7710), we chose to run the unannotated sequences through Goblet  using the vertebrate database and a BLAST cutoff of e-10. This increased the annotation level by 1736 sequences. Because the default criteria for a match in Goblet is lower than in AutoFACT and these sequences had already passed through AutoFACT, these GO annotations should be regarded as less reliable. GO terms were also identified for an additional 502 sequences by searching InterPro. Functionally annotated ESTs with GO annotations were classified using GOSlim and each category with more than 5 hits was plotted. Individual ESTs in each library (13,000 total) were also searched for GOSLIM terms and compiled by library.
Identification of microsatellite sequences
Short tandem repeats (2 – 5 bp) were detected using Tandem Repeats Finder .
Worm B, Barbier EB, Beaumont N, Duffy JE, Folke C, Halpern BS, Jackson JBC, Lotze HK, Micheli F, Palumbi SR, Sala E, Selkoe KA, Stachowicz JJ, Watson R: Impacts of biodiversity loss on ocean ecosystem services. Science. 2006, 314: 787-790. 10.1126/science.1132294.
Hendry CI, Martin-Robichaud DJ, Benfey TJ: Hormonal sex reversal of Atlantic halibut (Hippoglossus hippoglossus L.). Aquaculture. 2003, 219: 769-781. 10.1016/S0044-8486(02)00344-7.
Gara B, Shields RJ, McEvoy L: Feeding strategies to achieve correct metamorphosis of Atlantic halibut, Hippoglossus hippoglossus L., using enriched Artemia. Aquaculture Research. 1998, 29: 935-948. 10.1046/j.1365-2109.1998.29120935.x.
Hamre K, Opstad I, Espe M, Solbakken J, Hemre G-I, Pittman K: Nutrient composition and metamorphosis success of Atlantic halibut (Hippoglossus hippoglossus, L.) larvae fed natural zooplankton or Artemia. Aquaculture Nutrition. 2002, 8: 139-148. 10.1046/j.1365-2095.2002.00201.x.
Okada N, Takagi Y, Seikai T, Tanaka M, Tagawa M: Asymmetrical development of bones and soft tissues during eye migration of metamorphosing Japanese flounder, Paralichthys olivaceus. Cell Tissue Research. 2001, 304: 59-66. 10.1007/s004410100353.
Naess T, Lie O: A sensitive period during first feeding for the determination of pigmentation pattern in Atlantic halibut, Hippoglossus hippoglossus L., juveniles: the role of diet. Aquaculture Research. 1998, 29: 925-934. 10.1046/j.1365-2109.1998.29120925.x.
Venizelos A, Benetti D: Pigment abnormalities in flatfish. Aquaculture. 1999, 176: 181-188. 10.1016/S0044-8486(99)00060-5.
Saele O, Solbakken JS, Watanabe K, Hamre K, Pittman K: The effect of diet on ossification and eye migration in Atlantic halibut larvae (Hippoglossus hippoglossus L.). Aquaculture. 2003, 220: 683-696. 10.1016/S0044-8486(02)00584-7.
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.
Marra MA, Hillier L, Waterston RH: Expressed sequence tags – ESTablishing bridges between genomes. Trends in Genetics. 1998, 14: 4-7. 10.1016/S0168-9525(97)01355-3.
Serapion J, Kucuktas H, Feng J, Liu Z: Bioinformatic mining of type I microsatellites from expressed sequence tags of channel catfish (Ictalurus punctatus). Marine Biotechnology. 2004, 6: 364-377. 10.1007/s10126-003-0039-z.
Diatchenko L, Lau Y-FC, Campbell AP, Chenchik A, Moqadam F, Huang B, Lukyanov S, Lukyanov K, Gurskaya N, Sverdlov ED, Siebert PD: Suppression subtractive hybridization: A method for generating differentially regulated or tissue-specific cDNA probes and libraries. Proceedings of the National Academy of Science of the USA. 1996, 93: 6025-6030. 10.1073/pnas.93.12.6025.
Nam B-H, Yamamoto E, Hirono I, Aoki T: A survey of expressed genes in the leukocytes of Japanese flounder, Paralichthys olivaceus, infected with Hirame rhabdovirus. Developmental and Comparative Immunology. 2000, 24: 13-24. 10.1016/S0145-305X(99)00058-0.
Kono T, Sakai M: The analysis of expressed genes in the kidney of Japanese flounder, Paralichthys olivaceus, injected with the immunostimulant peptidoglycan. Fish and Shellfish Immunology. 2001, 11: 357-366. 10.1006/fsim.2000.0314.
Davey GC, Caplice NC, Martin SA, Powell R: A survey of genes in the Atlantic salmon (Salmo salar) as identified by expressed sequence tags. Gene. 2001, 263: 121-130. 10.1016/S0378-1119(00)00587-4.
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 and Shellfish Immunology. 2000, 10: 623-630. 10.1006/fsim.2000.0279.
Aoki T, Nam B-H, Hirono I, Yamamota E: Sequences of 596 cDNA clones (565,977 bp) of Japanese flounder (Paralichthys olivaceus) leukocytes infected with Hirame rhabdovirus. Marine Biotechnology. 1999, 1: 477-488. 10.1007/PL00011804.
Bayne CJ, Gerwick L, Fujiki K, Nakao M, Yano T: Immune-relevant (including acute phase) genes identified in the livers of rainbow trout, Oncorhynchus mykiss, by means of suppression subtractive hybridization. Developmental and Comparative Immunology. 2001, 25: 205-217. 10.1016/S0145-305X(00)00057-4.
Douglas SE, Belanger JC, Tsoi SCM, Williams J, Ewart KV: Identification of differentially expressed genes involved in Atlantic salmon (Salmo salar) host-pathogen interaction using microarray technology. Sixth International Marine Biotechnology Conference – Aquatic Genomics Symposium; Sept. 20–27, 2003. 2003, Tokyo, Japan. Marine Biotechnology, S335-339.
Gerwick L, Steinhauer R, LaPatra SE, Sandell T, Ortuno J, Hajiseyedjavadi N, Bayne CJ: The acute phase response of rainbow trout (Oncorhynchus mykiss) plasma proteins to viral, bacterial and fungal inflammatory agents. Fish and Shellfish Immunology. 2002, 12: 229-242. 10.1006/fsim.2001.0367.
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. Marine Biotechnology. 2002, 4: 526-536. 10.1007/s10126-002-0067-0.
Kono T, Sakai M, LaPatra SE: Expressed sequence tag analysis of kidney and gill tissues from rainbow trout (Oncorhynchus mykiss) infected with Infectious Hematopoietic Necrosis virus. Marine Biotechnology. 2000, 2: 493-498.
Martin SA, Blaney SC, Houlihan DF, Secombes CJ: Transcriptome response following administration of a live bacterial vaccine in Atlantic salmon (Salmo salar). Molecular Immunology. 2006, 43: 1900-1911. 10.1016/j.molimm.2005.10.007.
Savan R, Sakai M: Analysis of expressed sequence tags (EST) obtained from common carp, Cyprinus carpio L., head kidney cells after stimulation by two mitogens, lipopolysaccharide and concanavalin-A. Comparative Biochemistry and Physiology Part B. 2002, 131: 71-82. 10.1016/S1096-4959(01)00488-2.
Suzuki T, Shin IT, Kohara Y, Kasahara M: Transcriptome analysis of hagfish leukocytes: a framework for understanding the immune system of jawless fishes. Developmental and Comparative Immunology. 2004, 28: 993-1003. 10.1016/j.dci.2004.04.005.
Tsoi SC, Ewart KV, Penny S, Melville K, Liebscher RS, Brown LL, Douglas SE: Identification of immune-relevant genes from Atlantic salmon using suppression subtractive hybridization. Marine Biotechnology. 2004, 6: 199-214. 10.1007/s10126-002-0101-2.
Chini V, Rimoldi S, Terova G, Saroglia M, Rossi F, Bernardini G, Gornati R: EST-based identification of genes expressed in the liver of adult seabass (Dicentrarchus labrax, L.). Gene. 2006, 376: 102-106. 10.1016/j.gene.2006.02.010.
Park KC, Osborne JA, Tsoi SC, Brown LL, Johnson SC: Expressed sequence tags analysis of Atlantic halibut (Hippoglossus hippoglossus) liver, kidney and spleen tissues following vaccination against Vibrio anguillarum and Aeromonas salmonicida. Fish and Shellfish Immunology. 2005, 18: 393-415. 10.1016/j.fsi.2004.10.003.
Fernandes JM, Mackenzie MG, Elgar G, Suzuki Y, Watabe S, Kinghorn JR, Johnston IA: A genomic approach to reveal novel genes associated with myotube formation in the model teleost, Takifugu rubripes. Physiological Genomics. 2005, 22: 327-338. 10.1152/physiolgenomics.00087.2005.
Ju Z, Karsi A, Kocobas 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: 373-382. 10.1016/S0378-1119(00)00491-1.
Karsi A, Li P, Dunham R, Liu ZJ: Transcriptional activities in the pituitaries of channel catfish before and after induced ovulation by injection of carp pituitary extract as revealed by expressed sequence tag analysis. Journal of Molecular Endocrinology. 1998, 21: 121-129. 10.1677/jme.0.0210121.
Zhou L, Yao B, Xia W, Li CJ, Wang Y, Shi YH, Gui JF: EST-based identification of genes expressed in the hypothalamus of male orange-spotted grouper (Epinephelus coioides). Aquaculture. 256: 129-139. 10.1016/j.aquaculture.2006.01.022.
Magalhaes GS, Junqueira-de-Azevedo IL, Lopes-Ferreira M, Lorenzini DM, Ho PL, Moura-da-Silva AM: Transcriptome analysis of expressed sequence tags from the venom glands of the fish Thalassophryne nattereri. Biochimie. 2006, 88: 693-699. 10.1016/j.biochi.2005.12.008.
Gene Ontology. [http://www.geneontology.org/]
Clusters of Orthologous Genes. [http://www.ncbi.nlm.nih.gov/COG/]
Kyoto Encyclopedia of Genes and Genomes. [http://www.genome.jp/kegg/]
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. Journal of Experimental Biology. 2006, 209: 2432-2441. 10.1242/jeb.02269.
Murray HM, Perez-Casanova JC, Gallant JW, Johnson SC, Douglas SE: Trypsinogen expression during the development of the exocrine pancreas in winter flounder (Pleuronectes americanus). Comparative Biochemistry and Physiology Part A. 2004, 138: 53-59. 10.1016/j.cbpb.2004.02.020.
Murray HM, Gallant JW, Johnson SC, Douglas SE: Cloning and expression analysis of three digestive enzymes from Atlantic halibut during early development: predicting gastrointestinal functionality. Aquaculture. 2006, 252: 394-408. 10.1016/j.aquaculture.2005.03.030.
Song HD, Sun XJ, Deng M, Zhang GW, Zhou Y, Wu XY, Sheng Y, Chen Y, Ruan Z, Jiang CL, Fan HY, Zon LI, Kanki JP, Liu TX, Look AT, Chen Z: Hematopoietic gene expression profile in zebrafish kidney marrow. Proceedings of the National Academy of Science of the USA. 2004, 101: 16240-16245. 10.1073/pnas.0407241101.
Sekar V: A rapid screening procedure for the identification of recombinant bacterial clones. BioTechniques. 5: 11-13.
Huang X, Madan A: CAP3: A DNA Sequence Assembly Program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.
Koski LB, Gray MW, Lang BF, Burger G: AutoFACT: An Auto matic F unctional A nnotation and C lassification T ool. BMC Bioinformatics. 2005, 6: 151-161. 10.1186/1471-2105-6-151.
Groth D, Lehrach H, Hennig S: GOblet: a platform for Gene Ontology annotation of anonymous sequence data. Nucleic Acids Research. 2004, 32: W313-317. 10.1093/nar/gkh406.
Benson G: Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Research. 1999, 27: 573-580. 10.1093/nar/27.2.573.
This project entitled: PLEUROGENE :Flatfish genomics – Enhancing commercial culture of Atlantic halibut and Senegal sole, was funded by the Genome Canada-Genome España joint program. We thank Debbie Martin-Robichaud (St. Andrews Biological Station) and Scotian Halibut Ltd., Clark's Harbour, Nova Scotia for tissue and larval samples of Atlantic halibut and Jeff Gallant and Harry Murray for storage and inventory of samples. Sequencing by The Atlantic Genome Centre, Halifax, Nova Scotia, a partnership between Genome Atlantic and the Canadian National Research Council Institute for Marine Biosciences, is gratefully acknowledged. This is NRC publication number 2007-42694.
SED and MER conceived of and designed the project. LCK constructed cDNA libraries. JK and MER carried out bioinformatic analyses. SED wrote the manuscript. MER and SED edited the manuscript. All authors read and approved the manuscript.
Susan E Douglas, Leah C Knickle, Jennifer Kimball and Michael E Reith contributed equally to this work.