An EST-based approach for identifying genes expressed in the intestine and gills of pre-smolt Atlantic salmon (Salmo salar)
© Hagen-Larsen et al; licensee BioMed Central Ltd. 2005
Received: 15 June 2005
Accepted: 01 December 2005
Published: 01 December 2005
The Atlantic salmon is an important aquaculture species and a very interesting species biologically, since it spawns in fresh water and develops through several stages before becoming a smolt, the stage at which it migrates to the sea to feed. The dramatic change of habitat requires physiological, morphological and behavioural changes to prepare the salmon for its new environment. These changes are called the parr-smolt transformation or smoltification, and pre-adapt the salmon for survival and growth in the marine environment. The development of hypo-osmotic regulatory ability plays an important part in facilitating the transition from rivers to the sea. The physiological mechanisms behind the developmental changes are largely unknown. An understanding of the transformation process will be vital to the future of the aquaculture industry. A knowledge of which genes are expressed prior to the smoltification process is an important basis for further studies.
In all, 2974 unique sequences, consisting of 779 contigs and 2195 singlets, were generated for Atlantic salmon from two cDNA libraries constructed from the gills and the intestine, accession numbers [Genbank: CK877169-CK879929, CK884015-CK886537 and CN181112-CN181464]. Nearly 50% of the sequences were assigned putative functions because they showed similarity to known genes, mostly from other species, in one or more of the databases used. The Swiss-Prot database returned significant hits for 1005 sequences. These could be assigned predicted gene products, and 967 were annotated using Gene Ontology (GO) terms for molecular function, biological process and/or cellular component, employing an annotation transfer procedure.
This paper describes the construction of two cDNA libraries from pre-smolt Atlantic salmon (Salmo salar) and the subsequent EST sequencing, clustering and assigning of putative function to 1005 genes expressed in the gills and/or intestine.
The number of known genes in Atlantic salmon (Salmo salar) is limited compared with the very large numbers known in mammals. However, the salmon is so important as an aquaculture species that a lot of effort has been put into acquiring knowledge of its genome compared to many of the other aquatic species. Although only 14 full-length genes are listed for Atlantic salmon at NCBI , there are 111 457-nucleotide sequences (September 2005), many containing complete cDNA sequences. The genome and biology of the Atlantic salmon are complex and the genome has not been well characterised. The salmonidae species underwent a genome duplication event about 25–100 million years ago [2, 3]. This event resulted in a tetraploid genome, but since then the salmon has been gradually returning to the diploid condition. This process is still continuing, and the salmon now has a partially tetraploid genome. Although many genes have been lost after the duplication event , the fact that regions of the genome are duplicated and almost identical does complicate the hunt for genes in this species. It means that there may be four almost identical copies of a gene in the genome, but it is an open question whether all four copies are active. Furthermore, the Atlantic salmon goes through a process known as smoltification as a step in its maturation. The salmon spends the first part of its life in fresh water before migrating to the ocean, where it lives until returning to rivers to spawn. Smoltification  involves synchronised morphological, physiological and behavioural changes that enable the young salmon (parr) to survive in the ocean and to grow and migrate normally. Since smoltification occurs while the fish are still in fresh water, they are pre-adapted to the marine environment. Given this enormous change in both fish biology and environment, it seems reasonable to suggest that the pre-smolt expresses a different set of genes from the smolt. It is therefore important to identify the genes that are expressed both before and after the smoltification process and at each of the different developmental stages in order to obtain a complete set of genes expressed in the salmon genome. One strategy for identifying the genes expressed in specific life stages and tissues is to use expressed sequence tags or ESTs. These are short stretches of single pass sequences obtained from sequencing cDNA [6, 7], and are widely used for gene discovery, mapping, polymorphism analysis, expression studies and gene prediction. Gene discovery methods using ESTs include hunting for new members of gene families in the same species (paralogues), for functionally equivalent genes in other species (orthologues), or even for alternatively spliced forms of known genes. ESTs are also used to predict or refine computational predictions of the location of genes in genomic DNA. Recently, a lot of effort has been put into acquiring ESTs from Atlantic salmon, and over 108,242 EST sequences (September 2005), mainly from smoltified fish, have been deposited in GenBank [8–11]The aim of the present study was to increase the number of ESTs available from pre-smolts in order to identify genes that are expressed in the early life stages of Atlantic salmon. We have focused on the gills and the intestine.
Fraction sizes and number of sequences for the gills and intestine cDNA libraries
Mean insert sizeb
Number of hits observed in the different databases
Gene Ontologya(Molecular function) of the sequences with a significant blastx hit in SwissProt
Nr of Hits
%glutathione-disulfide reductase activity
%metal ion binding
%nucleic acid binding
%pyridoxal phosphate binding
%small protein activating enzyme activity
%small protein conjugating enzyme activity
%heat shock protein activity
%enzyme regulator activity
%caspase regulator activity
%enzyme activator activity
%enzyme inhibitor activity
%GTPase regulator activity
%kinase regulator activity
%nitric-oxide synthase regulator activity
%microtubule motor activity
%signal transducer activity
%receptor signaling protein activity
%structural molecule activity
%extracellular matrix structural constituent
%structural constituent of bone
%structural constituent of cytoskeleton
%structural constituent of eye lens
%structural constituent of muscle
%structural constituent of ribosome
%transcription regulator activity
%RNA polymerase I transcription factor activity
%RNA polymerase II transcription factor activity
%transcription cofactor activity
%transcription factor activity
%transcriptional repressor activity
%translation regulator activity
%translation factor activity, nucleic acid binding
%amine/polyamine transporter activity
%auxiliary transport protein activity
%carbohydrate transporter activity
%channel/pore class transporter activity
%drug transporter activity
%electron transporter activity
%intracellular transporter activity
%ion transporter activity
%lipid transporter activity
%neurotransmitter transporter activity
%organic acid transporter activity
%oxygen transporter activity
%peptide transporter activity
%protein transporter activity
Table 3 shows that the sequences can be classified in 12 of the major categories for molecular function, the two largest groups being catalytic activity (26%) and binding (36%).
The total number of gill and intestine clones available was 8256. Table 1 shows the sampled means of the insert sizes after fractionation, which varied between 700 and 1500 base pairs. cDNA cloning rarely produces full-length gene products but rather products containing varying proportions of the whole sequence. 5'sequencing may overcome this problem since clustering can be used to cover larger parts of the transcripts. In this study, after trimming and quality control followed by clustering, approximately 65% of the trimmed sequences from both libraries ended up in contigs.
Sequence annotation was carried out using GO terms. A fairly strict criterion of E-value below 1.0 × 10-15 for a significant hit was used both for Blastx run against the Swiss-Prot and nr protein databases, and blastn against the nr nucleotide database. This rather low value was chosen to reduce the number of hits on very remote homologous genes with functions unrelated or only remotely related to the sequences in the current libraries. For blastx against the PDB database, the significance level chosen was an E-value below 1.0 × 10-10, and 489 or 16% of the sequences gave significant hits in this database. Almost all of these showed sequence similarity above 30%, indicating that they might be suitable targets for homology modelling of protein 3D structure if the hit produces a biologically meaningful alignment.
In the two libraries, we discovered several genes that may be involved in the smoltification process, genes involved in cell homeostasis and genes coding for hormones that may influence salmon maturation. Certain genes were of particular interest, for example the glucagon-family neuropeptide precursor, tyrosine 3/ tryptophan 5-monooxygenase activities and thyroid receptor interactin protein 6 (TRIP6). In addition, five different genes related to growth factor activity were also present: granulin, neuregulin 1, syntenin, Bmp2 and pleiotrohin 1. All of the above-mentioned genes are thought to play a part in the smoltification process in salmon, along with the osmoregulatory genes. However, it is possible that some of the more interesting genes involved in the smoltification process have not yet been annotated and characterised and are therefore listed as unknown. Approximately 50% of the sequences gave no significant hits in the blastx-blastn pipeline. A large proportion of these gave hits with higher E-values to more remote homologous sequences, but a significant number, at least 500–1000 sequences, are of completely unknown function at present. However, the 5' UTR of fish may differ significantly from that of mammals and terrestrial organisms. This means that the failure to match 500–1000 sequences does not confirm that the genes are unique to fish or markedly different from the coding region of transcripts of other organisms. A number of factors, including the 5'UTR, may help to explain why no similarity was found for these sequences, but the unusual biology of salmon and the fact that the sequences are from pre-smolt fish may be contributory factors. Pre-smolt salmon live in fresh water, and presumably different sets of genes are turned on and off in response to the two different habitats, the freshwater and the marine environment. In addition, the lack of similarity to already known sequences in other vertebrates could influence these results. Although large amounts of sequence data from other species have been generated, a large fraction of the genes in higher organisms is still uncharacterised. It therefore seems likely that the lack of similarity of these ESTs to known genes is partly explained by the fact that many genes are still uncharacterised.
This work has made available approximately 3000 sequences from pre-smolt salmon, one third of which have been annotated for function, biological process and cellular component. Several hundred sequences code for proteins that show enough sequence-similarity to proteins of known structure to be suitable candidates for homology modelling. The results of this study will be a valuable resource for future studies of Atlantic salmon biology since it is now possible to search for genes sequenced in salmon rather than having to use comparative data from other species.
1. Construction of cDNA libraries
Pre-smolt Atlantic salmon (Salmo salar) from a commercial aquaculture-bred stock (National Norwegian Breeding Programme for salmon, Norway), weighing 50–100 g, were collected in late September for the construction of two cDNA libraries from the intestine and gills. The entire intestine and the whole of the gills were used to construct the libraries. Tissue samples of 400 mg were used for total RNA extraction, which was performed using the guanidine thiocyanate-phenol-chloroform extraction method . mRNA was further isolated using the Poly (A) Quick® mRNA isolation kit (Stratagene Cloning Systems, California, USA). Three μg of total RNA and 0.1 μg of mRNA were run on an agarose gel to investigate the integrity/quality of the samples. cDNA synthesis was performed using the pBluescript® II XR cDNA library construction kit (Stratagene Cloning Systems, California, USA), which generates directional cDNA libraries using an oligo (dT)18 primer. For each library the cDNA was size fractionated using a gel filtration column, and three fractions containing most of the cDNA were ligated into the plasmid vector and subsequently transformed into XL10 Gold ultra competent cells.
2. Pre-screening of the cDNA libraries
Pre-screening was performed to separate abundant genes from rare ones. From each of the three cDNA fractions from both tissues, approximately 50 000 colonies were plated out on agar plates (140 mm in diameter), with an average of 2000 colonies per plate. The colonies were transferred to 132 mm BIOTRANS nylon membranes, pore size 1.2 Micron, before hybridisation to total cDNA. Total cDNA from the intestine was used to screen the intestine filters and total cDNA from the gills to screen the gill filters. The cDNA was labelled using α32P-dCTP in a PCR assay. In order not to label the vector, the plasmids containing the cDNA were digested with XhoI/EcoRI to open up the circle. Then M13–20 primer was used in a 'one way' PCR to label the insert. The colonies that hybridised were presumed to be the more abundant ones in the libraries. A total of 3840 non-hybridising clones (the presumed rare genes) and 288 hybridising clones were picked from each of the gill and intestine libraries and gridded in duplicate into 384 well microtiterplates. The plates were grown overnight in LB, then glycerol was added and the plates were sealed and stored at -80°C.
3. Enrichment of pBluescript cDNA libraries
The remaining cDNA was amplified in semi-solid fashion, as described by the manufacturer (Stratagene Cloning Systems, California, USA). Amplification in suspension allowed three-dimensional, uniform growth, which reduced the likelihood of under-representation of particular clones. Up to 5 × 105 cfu/bottle of primary library was added to each bottle of the semi-solid suspension, before incubation on an ice-water bath for an hour to achieve the required semi-solid state of the medium, followed by incubation for 40–45 hours at 30°C. Afterwards the contents of the bottles were spun down and resuspended in 50 ml of 2xLB-glycerol (12.5%) for the gill library and 100 ml for the intestine library, ending up with 8 × 109 clones from the gills and 3.1 × 1010 clones from the intestine. Total volumes of 1 ml were aliquoted into tubes and subsequently frozen at -80°C for future use.
4. DNA isolation procedures
A few microlitres of the pre-screened cDNA library was cultured overnight in LB before the DNA was isolated using the QIAprep 96 Turbo Miniprep kit from QIAGEN in a QIAvac96 according to the manufacturer's manual.
5. Size determination of the library fractions
To determine the insert sizes of the clones from the different size fractions, a rough estimate was made for 192 randomly picked clones from each fraction. The insert size was determined by digesting the insert out of the vector using the restriction enzymes EcoRI and XhoI, followed by separation on an agarose gel. The gels were scanned on a Typhoon 9410 Imager and analysed using the ImageQuant TL software, both from Amersham Biosciences.
6. DNA sequencing
Sequencing was done from the 5'-end using T3 as sequencing primer. 5'-sequencing was chosen in order to assign functional annotation to as many transcripts as possible. The sequencing reactions were performed using the ABI PRISM® BigDye™ Terminators Cycle Sequencing Kit (Applied Biosystems), and run on the ABI 377 (Applied Biosystems) or on the Megabace 1000 (Amersham Pharmacia).
7. Clustering, sequence comparison and annotation transfer
Phred [13, 14] and cross_match (Green P: unpublished) were used for basecalling and trimming of vector, respectively. The sequences were masked for repeats against a Danio rerio repeat library , and sequences containing contaminants or that were of poor quality were removed. The remaining high quality sequences were assembled and clustered with the Sequencher 4.1 package. In order to examine sequence similarity to known genes, NCBI blast sequence alignment  was performed against the NCBI databases  (on 08.03.2004). For each sequence, blastx was run against the PDB, Swiss-Prot and nr protein sequence databases, and blastn was run against the nr nucleotide sequence database. We defined a significant database hit as having an expectation value (E-value) below 1.0 × 10-15 for all sequence alignments, except for blastx against PDB where 1.0 × 10-10 was used. All sequences that gave a significant blastx hit in Swiss-Prot were annotated by annotation transfer, applying the Gene Ontology (GO) [18, 19] assignments for the UniProt database produced by the GOA project of the European Bioinformatics Institute. The gene_association.goa_uniprot database of 26.04.2004 was used  together with GO terms from the GO release of 11.05.2004 . The sequences were annotated on the basis of the single best hit in the Swiss-Prot database.
The authors would like to thank Agate Noer and Kristin Vekterud for their assistance with library construction and sequencing. This work was supported by the EU FAIR Programme SALGENE project (FAIR CT98-4314); and grant 130162/130 as part of a project entitled "Strategic QTL research plan for disease resistance in Atlantic salmon and cattle" and grant 139617/140 as part of a project entitled "Salmon Genome Project (SGP)", both supported by the Research Council of Norway.
- NCBI - National Center for Biotechnology Information [http://www.ncbi.nlm.nih.gov/entrez].
- Allendorf FWTGH: Tetraploidy and the evolution of salmonid fishes. Evolutionary genetcs of fishes. Edited by: Turner BJ. 1984, New York , Plenum Press, 1-53.View ArticleGoogle Scholar
- Ohno S, Wolf U, Atkin NB: Evolution from fish to mammals by gene duplication. Hereditas. 1968, 59 (1): 169-187.PubMedView ArticleGoogle Scholar
- Wolfe KH: Yesterday's polyploids and the mystery of diploidization. Nat Rev Genet. 2001, 2 (5): 333-341. 10.1038/35072009.PubMedView ArticleGoogle Scholar
- Boeuf G: Salmonid smolting: a pre-adaptation to the oceanic environment. Fish Ecophysiology. Edited by: Rankin JCJFB. 1993, London , Chapman & Hall, pp 105-135.View 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 (5013): 1651-1656.PubMedView ArticleGoogle Scholar
- Marra MA, Hillier L, Waterston RH: Expressed sequence tags--ESTablishing bridges between genomes. Trends Genet. 1998, 14 (1): 4-7. 10.1016/S0168-9525(97)01355-3.PubMedView ArticleGoogle Scholar
- 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 (1-2): 121-130. 10.1016/S0378-1119(00)00587-4.PubMedView ArticleGoogle Scholar
- Martin SA, Caplice NC, Davey GC, Powell R: EST-based identification of genes expressed in the liver of adult Atlantic salmon (Salmo salar). Biochem Biophys Res Commun. 2002, 293 (1): 578-585. 10.1016/S0006-291X(02)00263-2.PubMedView 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 (3): 478-490. 10.1101/gr.1687304.PubMedPubMed CentralView ArticleGoogle Scholar
- SGP - Salmon genome project [http://www.salmongenome.no].
- Chomczynski P, Sacchi N: Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987, 162 (1): 156-159. 10.1016/0003-2697(87)90021-2.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
- 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
- Repeat Masker [http://ftp.genomewashington.edu/RM/RepeatMasker.html].
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410. 10.1006/jmbi.1990.9999.PubMedView ArticleGoogle Scholar
- NCBI - National Center for Biotechnology Information [ftp://ftp.ncbi.nih.gov/blast/db].
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.PubMedPubMed CentralView ArticleGoogle Scholar
- Ashburner M, Lewis S: On ontologies for biologists: the Gene Ontology--untangling the web. Novartis Found Symp. 2002, 247: 66-80; discussion 80-3, 84-90, 244-52.PubMedView ArticleGoogle Scholar
- European Bioinformatics Institute [http://www.ebi.ac.uk/GOA].
- European Bioinformatics Institute [ftp://ftp.ebi.ac.uk/pub/databases/GO/goa/UNIPROT].
- The Gene Ontology [http://www.geneontology.org].
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.