<title>Introduction</title> <p>Osmeriformes are close relatives of the Salmoniformes. The Osmeroidei and Salmoniformes clades are separated by at least 200 My with the Salmonidae having undergone at least one genome duplication event since their divergence (Ohno et al. ##REF##5662632##1968##; Allendorf and Thorgaard ##UREF##0##1984##; Ishiguro et al. ##REF##12742752##2003##). Osmerids, such as the rainbow smelt, have less than half the amount of genomic deoxyribonucleic acid (DNA) as salmonids and are considered to represent the ancestral state prior to the salmonid genome duplication (Ohno ##UREF##1##1970##). The bacterial artificial chromosome (BAC) resources provide a unique opportunity to study differences between orthologs (and paralog numbers), as well as chromosome alterations (through syntenic BAC comparisons), between species.</p> <p>Rainbow smelt and salmon are closely related and have similar life histories; however, they represent different scales of cold adaptation. Smelt, unlike salmonids, are completely cold adapted, fully freeze resistant, active, and feed voraciously at low temperature (reviewed by Driedzic and Ewart ##REF##15544960##2004##). Smelt have adapted to these conditions by producing and accumulating an antifreeze protein (AFP), glycerol, trimethylamine <italic>N</italic>-oxide, and urea that each contribute to lowering the freezing point of the fish (Driedzic and Ewart ##REF##15544960##2004##). Glycerol can be synthesized from glucose or amino acid precursors in smelt (Walter et al. ##REF##16506225##2006##). It is interesting to note that the abbreviated pathway by which glycerol is produced from amino acids is well known in mammals and termed glyceroneogenesis (Hanson and Reshef ##REF##14739071##2003##). The seasonal accumulation of glycerol and AFP do not appear to be linked transcriptionally or metabolically (Liebscher et al. ##REF##16555199##2006##).</p> <p>To isolate and identify genes involved in cold adaptation and other physiological functions, we have constructed a large BAC clone and BAC library and generated a large expressed sequence tag (EST) clone and sequence database. Our large smelt EST resource facilitates further gene discovery and determination of how genes (proteins) evolve new functions and processes between species and provide an opportunity for future microarray and microsatellite studies.</p>
<title>Materials and Methods</title> <title>BAC Resources</title> <p>To provide a genomic clone resource, a BAC library, CHORI-74, was prepared following Osoegawa et al. ##REF##9740665##(1998##; Children’s Hospital Oakland Research Institute [CHORI], Oakland, CA, USA). High-molecular-weight DNA was isolated from blood cells from a female individual, ID number 4, partially digested with a combination of <italic>Eco</italic>RI restriction and <italic>Eco</italic>RI methylase enzymes and then size fractionated by pulsed-field gel electrophoresis. DNA fragments were cloned into the pBAC-GMR vector. The library was arrayed into 144 384-well microtiter plates and gridded onto three 22 × 22-cm nylon high-density filters. Each hybridization membrane represents more than 18,000 distinct BAC clones, stamped in duplicate.</p> <title>EST Resources</title> <p>To identify genes in <italic>Osmerus mordax</italic>, complementary DNA (cDNA) libraries were constructed from ribonucleic acid (RNA) isolated from samples obtained from the Memorial University of Newfoundland Ocean Sciences Center, Logy Bay, NL, Canada. Smelt were collected in October 2002 in Long Harbour, Placentia Bay, Newfoundland, then transferred to the Ocean Sciences Centre, held under a natural photoperiod, and fed chopped herring twice per week. Fish were maintained in seawater at ambient temperature, which followed a profile similar to that presented in Lewis et al. (##REF##15286915##2004##). Fish were sampled in January and April 2003. Brain, liver, head kidney, and spleen tissues were flash frozen and stored at −80°C until RNA extraction. Total RNA (Trizol reagent; Invitrogen, Carlsbad, CA, USA) or poly(A)+ RNA (FastTrack MAG kit; Invitrogen) was extracted from the flash-frozen tissues. Conventional libraries of low- and high-molecular-weight smelt brain, liver, kidney, and spleen cDNAs were individually constructed using pBluescript II XR cDNA library construction kits (Stratagene, La Jolla, CA, USA). Mixed tissue libraries were normalized by either the negative subtraction-based normalization method (Invitrogen; Research Genetics, California) or the duplex-specific nuclease normalization method (Evrogen, Moscow, Russia). The normalized libraries were directionally constructed in pCMV-Sport6.1 (Invitrogen) or pAL-17.3 (Evrogen) vectors.</p> <title>Bioinformatic Resources</title> <p>Plasmid DNAs were extracted and BigDye™ Terminator (ABI, Foster City, CA, USA) cycle sequenced on ABI 3730 sequencers using conventional procedures and the following primers: 5′-T<sub>18</sub>-3′, M13 forward (5′-GTAAAACGACGGCCAGT-3′), and M13 reverse (5′-AACAGCTATGACCATG-3′ or 5′-CAGGAAACAGCTATGAC-3′). Base calling and trimming of vector, poly-A tails, and low-quality regions were addressed as described by Rise et al. (##REF##14962987##2004##). Initial assembly of ESTs into contigs used PHRAP (<ext-link ext-link-type="uri" xlink:href="http://bozeman.mbt.washington.edu">http://bozeman.mbt.washington.edu</ext-link>), under stringent clustering parameters (minimum score = 100; repeat stringency = 0.99). A second-stage assembly used the consensus sequences (with quality scores) from the first stage and parameters of 96% repeat frequency and 300 minscore to build final contigs and consensus sequences. Assemblies using CAP3 (Huang and Madan ##REF##10508846##1999##) using default parameters of 75% identity over an area of 30 bp resulted in similar contigs. Contig consensus sequences and singleton sequences were aligned with nonredundant GenBank nucleotide and several amino acid sequence databases (Gene Ontology [GO], swissprot, Conserved Domain Database [CDD], and Uniref90) using BLASTN and BLASTX, respectively (Boguski et al. ##REF##8401577##1993##; Altschul et al. ##REF##9254694##1997##; Schwede et al. ##REF##12824332##2003##; Camon et al. ##REF##14681408##2004##; Harris et al. ##REF##14681407##2004##; Marchler-Bauer et al. ##REF##15608175##2005##; Kopp and Schwede ##REF##16381875##2006##). Using the swissprot database cross-reference, alignments of the second-stage contigs with entries in the database were used to assign GO terms to the contigs.</p> <p>The EST resources have been submitted to GenBank with the following accession numbers: for the normalized libraries, EL518196 to EL551831, and for the non-normalized libraries, CB484654 to CB484815, CN442489 to CN442491, CX349771 to CX351193, and EL517809 to EL518195. Sequence databases, assemblies, consensus sequences, tools such as BLAST, and sequence and consensus annotations are available at the Genomics Research on Atlantic Salmon Project website (<ext-link ext-link-type="uri" xlink:href="http://www.uvic.ca/cbr/grasp">http://www.uvic.ca/cbr/grasp</ext-link>).</p>
<title>Results and Discussion</title> <p>In the study, 52,410 BAC clones with an average insert size of 146 kb were obtained. Determination of the average insert size was calculated by taking one sample from each plate and, following minipreping and <italic>Not</italic>1 digestion, sizing by contour-clamped homogeneous electric field electrophoresis (CHORI). The insert size distribution is shown in Fig. ##FIG##0##1##. Given an estimated genome size of 0.69 pg (Hardie and Hebert ##REF##12897876##2003##), the BAC clone library represents approximately 11-fold genome coverage. These BAC clones enable us to isolate and characterize gene regions of interest and are available through CHORI BAC resources (<ext-link ext-link-type="uri" xlink:href="http://bacpac.chori.org/library.php?id=421">http://bacpac.chori.org/library.php?id=421</ext-link>).

</p> <p>At least 33,636 sequences were sequenced from the two normalized, mixed tissue libraries, and these were combined with 1,975 ESTs from the non-normalized libraries for a total of 35,611 EST sequences submitted to GenBank. Of the 12,159 second-stage contigs or transcripts assembled from 36,758 EST sequences, 6,139 had a BLASTX hit with <italic>E</italic> values less than 1e−10 to a well-annotated protein entry in the swissprot, CDD, or Uniref90 database (Table ##TAB##0##1##).

</p> <p>Alignments of the 6,139 contigs to entries in the swissprot database led to a total of 9,921 GO annotations to 2,500 different terms. The three ontologies comprising GO—molecular function, biological process, and cellular component—provided annotations for 3,534, 3,846, and 2,506 contigs, respectively. A further breakdown of the annotations is provided in Table ##TAB##1##2##. The complete GO hierarchy and the annotations corresponding to the contigs are available at <ext-link ext-link-type="uri" xlink:href="http://www.uvic.ca/cbr/grasp">http://www.uvic.ca/cbr/grasp</ext-link>.

</p> <p>For molecular function, 1,078, 1,635, and 497 contigs have been ascribed catalytic, nucleotide- and protein-binding, or regulator and transducer activities, respectively. The cellular component presents contigs that comprise various cellular regions, partitioning representatives to extra- or intracellular regions, as well as to mitochondrial, endoplasmic reticular, or nuclear regions. For the 3,846 contigs assigned a biological process, 1,370 represented metabolism of macromolecules, proteins, and lipids, and 2,869 represented cellular processes, such as reproductive, immune system, cell communication, cell cycle, proliferation, and development (including morphogenesis, differentiation and localization) processes (Table ##TAB##1##2##).</p> <p>When the 12,159 contigs were compared (BLASTN) to EST sequences in GenBank, 4,697 rainbow smelt transcripts aligned with an Atlantic salmon EST (<italic>E</italic> value less than 1e−25 over more than 200 bp) with an average identity of 86.2% (over an average of 431 bp), and 4,347 transcripts aligned with rainbow trout ESTs with an average identity of 86.1% (over 419 bp). These comparisons provide only a very general indication of the similarity between transcriptomes of rainbow smelt and salmonids, as assemblies contain both 5′- (generally genic regions) and 3′- (generally 3′-untranslated regions) transcript reads. However, these DNA sequence similarity values corroborate a more ancient separation of rainbow smelt and salmonid species than duplicated salmonid major histocompatibility complex class IA and B genes (Lukacs et al. ##REF##17651474##2007##) or growth hormone genes (McKay et al. ##REF##15284876##2004##). Comparisons of sequence identity between the Atlantic salmon gene duplicates are closer to one another (88% to 95%) than to any of the aligned smelt EST sequences (86%), consistent with an ancestral salmonid genome duplication hypothesis. Moreover, the high level of similarity between rainbow smelt ESTs and salmonid ESTs (86% identity) explains the observed high level of rainbow smelt cDNA hybridization to salmonid cDNA microarrays (Rise et al. ##REF##14962987##2004##; von Schalburg et al. ##REF##16164747##2005##).</p> <p>The primary function of the AFP in smelt tissues is likely to be freezing point depression, although roles for AFPs in low-temperature tolerance have also been suggested (reviewed by Inglis et al. ##REF##17168784##2006##). Seasonal expression of smelt AFP has been shown (Liebscher et al. ##REF##16555199##2006##). However, the tissue distribution of expression was unknown. Our liver libraries predominately contained type II AFP transcripts. In fact, sequences representing AFP clustered to one contig with the highest frequency of all genes in the smelt database. The AFP does not appear to be expressed in the brain, head kidney, or spleen libraries, suggesting that the liver is exclusive or predominant among these tissues in expressing AFP in smelt. Further insight into the evolution, diversity, and structure/function of the smelt AFP may arise from studies using the resources developed here.</p> <p>Cold adaptation is normally multifactorial, and it is likely that smelt have adaptations in addition to the known glycerol and AFP. Studies to identify other adaptations will draw largely on the resources presented here. The opportunity to further study low-temperature adaptation in this thoroughly cold adapted vertebrate may present unique opportunities for new applications in animal biology and in medicine.</p>
<title>Results and Discussion</title> <p>In the study, 52,410 BAC clones with an average insert size of 146 kb were obtained. Determination of the average insert size was calculated by taking one sample from each plate and, following minipreping and <italic>Not</italic>1 digestion, sizing by contour-clamped homogeneous electric field electrophoresis (CHORI). The insert size distribution is shown in Fig. ##FIG##0##1##. Given an estimated genome size of 0.69 pg (Hardie and Hebert ##REF##12897876##2003##), the BAC clone library represents approximately 11-fold genome coverage. These BAC clones enable us to isolate and characterize gene regions of interest and are available through CHORI BAC resources (<ext-link ext-link-type="uri" xlink:href="http://bacpac.chori.org/library.php?id=421">http://bacpac.chori.org/library.php?id=421</ext-link>).

</p> <p>At least 33,636 sequences were sequenced from the two normalized, mixed tissue libraries, and these were combined with 1,975 ESTs from the non-normalized libraries for a total of 35,611 EST sequences submitted to GenBank. Of the 12,159 second-stage contigs or transcripts assembled from 36,758 EST sequences, 6,139 had a BLASTX hit with <italic>E</italic> values less than 1e−10 to a well-annotated protein entry in the swissprot, CDD, or Uniref90 database (Table ##TAB##0##1##).

</p> <p>Alignments of the 6,139 contigs to entries in the swissprot database led to a total of 9,921 GO annotations to 2,500 different terms. The three ontologies comprising GO—molecular function, biological process, and cellular component—provided annotations for 3,534, 3,846, and 2,506 contigs, respectively. A further breakdown of the annotations is provided in Table ##TAB##1##2##. The complete GO hierarchy and the annotations corresponding to the contigs are available at <ext-link ext-link-type="uri" xlink:href="http://www.uvic.ca/cbr/grasp">http://www.uvic.ca/cbr/grasp</ext-link>.

</p> <p>For molecular function, 1,078, 1,635, and 497 contigs have been ascribed catalytic, nucleotide- and protein-binding, or regulator and transducer activities, respectively. The cellular component presents contigs that comprise various cellular regions, partitioning representatives to extra- or intracellular regions, as well as to mitochondrial, endoplasmic reticular, or nuclear regions. For the 3,846 contigs assigned a biological process, 1,370 represented metabolism of macromolecules, proteins, and lipids, and 2,869 represented cellular processes, such as reproductive, immune system, cell communication, cell cycle, proliferation, and development (including morphogenesis, differentiation and localization) processes (Table ##TAB##1##2##).</p> <p>When the 12,159 contigs were compared (BLASTN) to EST sequences in GenBank, 4,697 rainbow smelt transcripts aligned with an Atlantic salmon EST (<italic>E</italic> value less than 1e−25 over more than 200 bp) with an average identity of 86.2% (over an average of 431 bp), and 4,347 transcripts aligned with rainbow trout ESTs with an average identity of 86.1% (over 419 bp). These comparisons provide only a very general indication of the similarity between transcriptomes of rainbow smelt and salmonids, as assemblies contain both 5′- (generally genic regions) and 3′- (generally 3′-untranslated regions) transcript reads. However, these DNA sequence similarity values corroborate a more ancient separation of rainbow smelt and salmonid species than duplicated salmonid major histocompatibility complex class IA and B genes (Lukacs et al. ##REF##17651474##2007##) or growth hormone genes (McKay et al. ##REF##15284876##2004##). Comparisons of sequence identity between the Atlantic salmon gene duplicates are closer to one another (88% to 95%) than to any of the aligned smelt EST sequences (86%), consistent with an ancestral salmonid genome duplication hypothesis. Moreover, the high level of similarity between rainbow smelt ESTs and salmonid ESTs (86% identity) explains the observed high level of rainbow smelt cDNA hybridization to salmonid cDNA microarrays (Rise et al. ##REF##14962987##2004##; von Schalburg et al. ##REF##16164747##2005##).</p> <p>The primary function of the AFP in smelt tissues is likely to be freezing point depression, although roles for AFPs in low-temperature tolerance have also been suggested (reviewed by Inglis et al. ##REF##17168784##2006##). Seasonal expression of smelt AFP has been shown (Liebscher et al. ##REF##16555199##2006##). However, the tissue distribution of expression was unknown. Our liver libraries predominately contained type II AFP transcripts. In fact, sequences representing AFP clustered to one contig with the highest frequency of all genes in the smelt database. The AFP does not appear to be expressed in the brain, head kidney, or spleen libraries, suggesting that the liver is exclusive or predominant among these tissues in expressing AFP in smelt. Further insight into the evolution, diversity, and structure/function of the smelt AFP may arise from studies using the resources developed here.</p> <p>Cold adaptation is normally multifactorial, and it is likely that smelt have adaptations in addition to the known glycerol and AFP. Studies to identify other adaptations will draw largely on the resources presented here. The opportunity to further study low-temperature adaptation in this thoroughly cold adapted vertebrate may present unique opportunities for new applications in animal biology and in medicine.</p>
<p>Genomic resources in rainbow smelt (<italic>Osmerus mordax</italic>) enable us to examine the genome duplication process in salmonids and test hypotheses relating to the fate of duplicated genes. They further enable us to pursue physiological and ecological studies in smelt. A bacterial artificial chromosome library containing 52,410 clones with an average insert size of 146 kb was constructed. This library represents an 11-fold average coverage of the rainbow smelt (<italic>O. mordax</italic>) genome. In addition, several complementary deoxyribonucleic acid libraries were constructed, and 36,758 sequences were obtained and combined into 12,159 transcripts. Over half of these transcripts have been identified, several of which have been associated with cold adaptation. These basic resources show high levels of similarity (86%) to salmonid genes and provide initial support for genome duplication in the salmonid ancestor. They also facilitate identification of genes important to fish and direct us toward new technologies for other studies in fish biology.</p> <title>Keywords</title>
<p>This research was supported by Genome Canada and Genome British Columbia. We thank Stephen O’Leary (NRC IMB) for critical review of the manuscript. The authors thank Connie Short for animal husbandry and tissue collection. This is NRC publication number 2007-42762.</p> <title>Open Access</title> <p>This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.</p>
<fig id="Fig1"><label>Fig. 1</label><caption><p>The insert size distribution</p></caption></fig>
<table-wrap id="Tab1"><label>Table 1</label><caption><p>Rainbow smelt EST project summary</p></caption><table frame="hsides" rules="groups"><thead><tr><th> </th><th>Rainbow smelt</th></tr></thead><tbody><tr><td>Number of EST sequences<sup>a</sup>
</td><td char="." align="char">36,758</td></tr><tr><td>Assembly stage 1<sup>b</sup>
</td><td char="." align="char"> </td></tr><tr><td> Number of contigs (2+ seq)<sup>c</sup>
</td><td char="." align="char">9,044</td></tr><tr><td> Number of singletons<sup>d</sup>
</td><td char="." align="char">7,019</td></tr><tr><td> Number of transcripts<sup>e</sup>
</td><td char="." align="char">16,063</td></tr><tr><td> Max. contig size<sup>f</sup>
</td><td char="." align="char">185</td></tr><tr><td> Ave. contig size<sup>g</sup>
</td><td char="." align="char">2.29</td></tr><tr><td> Number of transcripts with BLASTX hits<sup>h</sup>
</td><td char="." align="char">8,457</td></tr><tr><td>Assembly stage 2</td><td char="." align="char"> </td></tr><tr><td> Number of transcripts<sup>i</sup>
</td><td char="." align="char">12,159</td></tr><tr><td> Number with BLASTX hits<sup>h</sup>
</td><td char="." align="char">6,139</td></tr><tr><td> Percent with hits<sup>j</sup>
</td><td char="." align="char">50</td></tr></tbody></table></table-wrap> <table-wrap id="Tab2"><label>Table 2</label><caption><p>GO annotation of contigs</p></caption><table frame="hsides" rules="groups"><thead><tr><th>GO accession</th><th>GO term name</th><th>Number of contigs</th></tr></thead><tbody><tr><td colspan="3">Biological process ontology</td></tr><tr><td>GO:0000003</td><td>Reproduction</td><td char="." align="char">88</td></tr><tr><td>GO:0002376</td><td>Immune system process</td><td char="." align="char">97</td></tr><tr><td>GO:0008152</td><td>Metabolic process</td><td char="." align="char">1,370</td></tr><tr><td>GO:0009987</td><td>Cellular process</td><td char="." align="char">2,869</td></tr><tr><td>GO:0016032</td><td>Viral reproduction</td><td char="." align="char">6</td></tr><tr><td>GO:0022414</td><td>Reproductive process</td><td char="." align="char">45</td></tr><tr><td>GO:0022610</td><td>Biological adhesion</td><td char="." align="char">85</td></tr><tr><td>GO:0032501</td><td>Multicellular organismal process</td><td char="." align="char">581</td></tr><tr><td>GO:0032502</td><td>Developmental process</td><td char="." align="char">671</td></tr><tr><td>GO:0040007</td><td>Growth</td><td char="." align="char">34</td></tr><tr><td>GO:0040011</td><td>Locomotion</td><td char="." align="char">15</td></tr><tr><td>GO:0043473</td><td>Pigmentation</td><td char="." align="char">5</td></tr><tr><td>GO:0048511</td><td>Rhythmic process</td><td char="." align="char">11</td></tr><tr><td>GO:0050896</td><td>Response to stimulus</td><td char="." align="char">373</td></tr><tr><td>GO:0051179</td><td>Localization</td><td char="." align="char">459</td></tr><tr><td>GO:0051234</td><td>Establishment of localization</td><td char="." align="char">373</td></tr><tr><td>GO:0051235</td><td>Maintenance of localization</td><td char="." align="char">8</td></tr><tr><td>GO:0051704</td><td>Multiorganism process</td><td char="." align="char">37</td></tr><tr><td>GO:0065007</td><td>Biological regulation</td><td char="." align="char">911</td></tr><tr><td colspan="3">Cellular component ontology</td></tr><tr><td>GO:0005576</td><td>Extracellular region</td><td char="." align="char">91</td></tr><tr><td>GO:0005623</td><td>Cell</td><td char="." align="char">2,401</td></tr><tr><td>GO:0031012</td><td>Extracellular matrix</td><td char="." align="char">25</td></tr><tr><td>GO:0031974</td><td>Membrane-enclosed lumen</td><td char="." align="char">162</td></tr><tr><td>GO:0031975</td><td>Envelope</td><td char="." align="char">105</td></tr><tr><td>GO:0032991</td><td>Macromolecular complex</td><td char="." align="char">337</td></tr><tr><td>GO:0043226</td><td>Organelle</td><td char="." align="char">1,272</td></tr><tr><td>GO:0044420</td><td>Extracellular matrix part</td><td char="." align="char">11</td></tr><tr><td>GO:0044421</td><td>Extracellular region part</td><td char="." align="char">66</td></tr><tr><td>GO:0044422</td><td>Organelle part</td><td char="." align="char">592</td></tr><tr><td>GO:0044456</td><td>Synapse part</td><td char="." align="char">6</td></tr><tr><td>GO:0045202</td><td>Synapse</td><td char="." align="char">11</td></tr><tr><td colspan="3">Molecular function ontology</td></tr><tr><td>GO:0003774</td><td>Motor activity</td><td char="." align="char">16</td></tr><tr><td>GO:0003824</td><td>Catalytic activity</td><td char="." align="char">1,078</td></tr><tr><td>GO:0005198</td><td>Structural molecule activity</td><td char="." align="char">53</td></tr><tr><td>GO:0005215</td><td>Transporter activity</td><td char="." align="char">140</td></tr><tr><td>GO:0005488</td><td>Binding</td><td char="." align="char">1,963</td></tr><tr><td>GO:0015457</td><td>Auxiliary transport protein activity</td><td char="." align="char">11</td></tr><tr><td>GO:0016209</td><td>Antioxidant activity</td><td char="." align="char">13</td></tr><tr><td>GO:0030188</td><td>Chaperone regulator activity</td><td char="." align="char">7</td></tr><tr><td>GO:0030234</td><td>Enzyme regulator activity</td><td char="." align="char">119</td></tr><tr><td>GO:0030528</td><td>Transcription regulator activity</td><td char="." align="char">210</td></tr><tr><td>GO:0031386</td><td>Protein tag</td><td char="." align="char">1</td></tr><tr><td>GO:0045182</td><td>Translation regulator activity</td><td char="." align="char">21</td></tr><tr><td>GO:0060089</td><td>Molecular transducer activity</td><td char="." align="char">140</td></tr></tbody></table></table-wrap>
<table-wrap-foot><p>
<sup>a</sup>Number of EST sequences for all of the species including those in GenBank as of September 2007</p><p>
<sup>b</sup>Assembly stage 1 refers to PHRAP assembly using parameters 99 repeat_frequency and 100 minscore</p><p>
<sup>c</sup>Number of contigs with two or more sequences</p><p>
<sup>d</sup>Number of contigs with one sequence</p><p>
<sup>e</sup>Total number of transcripts including singletons</p><p>
<sup>f</sup>The largest number of sequences that are contained within any single contig</p><p>
<sup>g</sup>The average number of sequences within each contig (includes singletons)</p><p>
<sup>h</sup>Number of transcripts that have a BLASTX hit of less than 1e−10 to swissprot databases</p><p>
<sup>i</sup>The number of transcripts that result from a reassembly of all stage 1 transcripts using PHRAP parameters 96 repeat_frequency and 300 minscore</p><p>
<sup>j</sup>Percent of stage 2 assembled transcripts that have a BLASTX hit</p></table-wrap-foot>
<graphic xlink:href="10126_2008_9089_Fig1_HTML" id="MO1"/>
[{"surname": ["Allendorf", "Thorgaard", "Turner"], "given-names": ["FW", "GH", "BJ"], "article-title": ["Tetraploidy and the evolution of salmonid fishes"], "source": ["Evolutionary genetics of fishes"], "year": ["1984"], "publisher-loc": ["New York"], "publisher-name": ["Plenum"], "fpage": ["1"], "lpage": ["53"]}, {"surname": ["Ohno"], "given-names": ["S"], "source": ["Evolution by gene duplication"], "year": ["1970"], "publisher-loc": ["Heidelberg, Germany"], "publisher-name": ["Springer"]}]
Mar Biotechnol (NY). 2008 Apr 2; 10(5):487