Smelt was the likely beneficiary of an antifreeze gene laterally transferred between fishes
© Graham et al.; licensee BioMed Central Ltd. 2012
Received: 8 May 2012
Accepted: 20 August 2012
Published: 25 September 2012
Type II antifreeze protein (AFP) from the rainbow smelt, Osmerus mordax, is a calcium-dependent C-type lectin homolog, similar to the AFPs from herring and sea raven. While C-type lectins are ubiquitous, type II AFPs are only found in a few species in three widely separated branches of teleost fishes. Furthermore, several other non-homologous AFPs are found in intervening species. We have previously postulated that this sporadic distribution has resulted from lateral gene transfer. The alternative hypothesis, that the AFP evolved from a lectin present in a shared ancestor and that this gene was lost in most species, is not favored because both the exon and intron sequences are highly conserved.
Here we have sequenced and annotated a 160 kb smelt BAC clone containing a centrally-located AFP gene along with 14 other genes. Quantitative PCR indicates that there is but a single copy of this gene within the smelt genome, which is atypical for fish AFP genes. The corresponding syntenic region has been identified and searched in a number of other species and found to be devoid of lectin or AFP sequences. Unlike the introns of the AFP gene, the intronic sequences of the flanking genes are not conserved between species. As well, the rate and pattern of mutation in the AFP gene are radically different from those seen in other smelt and herring genes.
These results provide stand-alone support for an example of lateral gene transfer between vertebrate species. They should further inform the debate about genetically modified organisms by showing that gene transfer between ‘higher’ eukaryotes can occur naturally. Analysis of the syntenic regions from several fishes strongly suggests that the smelt acquired the AFP gene from the herring.
KeywordsIce-binding Horizontal gene transfer Antifreeze protein Osmerus mordax
Some organisms that are unable to avoid subzero temperatures use antifreeze proteins (AFPs) that bind to ice and lower the non-equilibrium freezing point by a non-colligative mechanism [1–3]. This allows the organism to survive in a supercooled state in the presence of ice. AFPs have now been reported from disparate branches of the tree of life, including bacteria, fungi, diatoms, plants, insects and fish [3–8]. The structures of these proteins (even within one kingdom) are proving to be as diverse as their taxonomic distribution, but most have repetitive sequences and structures with a flat surface that binds to ice. In fish, four different classes of AFPs provide freezing resistance to their hosts . Two are linear and repetitive; the alpha-helical type I AFPs and the antifreeze glycoproteins (AFGPs), whereas the type II and type III AFPs are globular. The scattered distribution of these AFP types, which often defies taxonomic conventions, has been noted on many occasions and Figure 1 of the 1998 review by Cheng nicely demonstrates this conundrum .
The type II AFPs, which are the focus of this study, are derived from the carbohydrate-recognition domain of calcium-dependent (C-type) lectins . Those from herring (Clupea harengus) and rainbow smelt (Osmerus mordax) retain this calcium dependence, whereas those of the sea raven (Hemitripterus americanus) and longsnout poacher (Brachyopsis rostratus) do not [11–13]. The structures of herring and poacher AFPs were determined by X-ray crystallography [13, 14] and that of the sea raven AFP by NMR . They all contain five disulfide bridges, one to three more than found in other members of this lectin family, and are much more similar to each other than to any other family members [14, 16]. This led to the proposal that an ancestral AFP had arisen prior to the separation of the herring and smelt lineages , estimated to have taken place well over 200 million years ago [17–19]. The absence of this AFP within the genomes of intervening species was attributed to gene loss. However, several lines of evidence, discussed later, indicated to us a different scenario in which the gene was acquired via lateral gene transfer .
Type II AFP transcripts were cloned from both closely related and distantly related fish, leading to the curious observation that the AFPs of the rainbow smelt and the Japanese smelt (Hypomesus nipponensis, same family) are no more similar to each other than they are to the AFP of the herring (different superorder) [11, 16, 20]. Previously, we cloned the AFP genes of rainbow smelt and herring, and found that the high similarity observed in the coding regions extended to the intronic sequences, which were up to 97% identical. The rate of non-synonymous substitution was not significantly lower than for synonymous substitution, suggesting that strong selection could not explain this similarity. In addition, hybridization bands were only detected for herrings, smelt and sea raven when the DNA of a variety of species was subjected to Southern blot analysis. Together, this evidence suggested that the AFP gene was laterally transferred between fish species .
The amplification of advantageous genes is thought to be a common adaptive response to rapid environmental changes, such as those that occur at the onset of an ice age . This process has been clearly demonstrated in both bacteria and yeast subjected to other stresses [22, 23]. Gene duplication is many orders of magnitude more frequent than advantageous point mutations, although reversion to a single copy may occur once the genes are optimized for their new function . Not surprisingly, the vast majority of AFP-producing organisms contain multiple copies of their AFP genes, usually encoding a number of different isoforms. In this, fish are no exception. There are more than 40 gene copies in both winter flounder (type I AFP)  and different pout species living near either pole (type III AFP) [26, 27]. In the case of AFGPs, they have arisen at least twice via convergent evolution and in both cases, not only are there multiple copies of the genes, but they encode polyprotein precursors that are processed to yield up to 46 mature AFGPs . The type II gene is also amplified in herring, as multiple bands were evident when the DNA from a single individual was subjected to Southern blotting. However, only a single band was evident with rainbow smelt DNA .
Lateral (or horizontal) gene transfer (LGT) is a well-accepted phenomenon in unicellular organisms and is particularly evident when bacteria acquire genes for antibiotic resistance . For multicellular organisms the transfer must be to the germ line and the mechanisms (transposons/viruses excepted) are not as obvious as for prokaryotes. We previously postulated that gene transfer may occur between fish during spawning , as fish release huge quantities of sperm into the water at this time . Indeed, sperm-mediated transformation is an established laboratory procedure for making transgenic organisms, including fish . Recently, other studies have supported LGT to metazoans , including the transfer of carotenoid biosynthetic genes from fungi to aphids  and cellulase genes from bacteria to nematodes .
The AFP genes in herring are likely dispersed and polymorphic as multiple bands were previously observed on a genomic Southern blot, whereas the single band observed in rainbow smelt suggested either a single-copy or a tandemly-duplicated gene . To explore this locus within the genome, we screened a rainbow smelt bacterial artificial chromosome (BAC) library (CHORI-74, ) for AFP genes. One of 18 overlapping BAC isolates was sequenced in its entirety and the single copy status of the AFP gene in the smelt genome was confirmed by quantitative PCR. Syntenic comparisons of the genes flanking the smelt AFP were made to sequenced fish genomes. The results remain consistent with and strongly reinforce the hypothesis that the smelt AFP gene was acquired via LGT.
BAC screening, sequencing and annotation
Probe and primer sequences
Type II AFP probe
Type II AFP R
AFP (80) qF2
AFP (80) qR2
PBAC-GMR primer F2
PBAC-GMR primer R2
O0139C19 SP6 F
O0139C19 SP6 R
O0139C19 T7 F
O0139C19 T7 R
RBP-3-2 5′ inner
RBP-3-2 3′ inner
RBP-3-2 5′ outer
RBP-3-2 3′ outer
RBP-3-2 5′ herring
RBP-3-2 3′ herring
A minimum tiling path was built for the 18 BACs using the BAC end sequences of O0139C19 as sequence tag sites (Additional file 1: Figure S1A). All but one was PCR positive using primers specific for the 3′ of this BAC (Table 1, O0139C19_T7_F & R), indicating that they contained the same AFP gene and came from the same region of the rainbow smelt genome. The remaining BAC (O0059H07) was positive for the 5′ end. As the AFP gene was centrally located in the sequenced BAC, the biased tiling path was a concern. So, the restriction fragment distributions of the three clones (O0068I08, O0119M24, O0139C19) were analyzed. Their BglII restriction fragment patterns were highly similar with 25 or more discernible bands (Additional file 2: Table S1). The in silico BglII digest of the assembled sequence was in complete agreement with the observed restriction pattern of the sequenced clone O0139C19 (Additional file 2: Table S1). Clones O0119M24 and O0068I08 shared all but the fragments found in the first ~10 and ~40 kb, respectively. End sequencing revealed that O0119M24 begins at the EcoRI site at base 6274 (EcoRI was used to generate these BAC clones via a partial digest of smelt genomic DNA), whereas O0068I08 began at the site at base 42446. Clone O0059H07, which overlapped at the 5′ end of the sequenced clone, had the expected unique 5′ end and terminated at the EcoRI site at base 93176. The dearth of clones extending upstream of the sequence insert is likely a result of the unequal distribution of EcoRI sites as 11 of the 28 sites are found within the first 25 kb of this 162 kb clone (Additional file 1: Figure S1B), which would decrease the odds of recovering clones spanning this section.
BAC clone O0139C19 contains fifteen genes
Annotation of the smelt BAC insert was aided by the large number of fish ESTs in the database, and the only gene that could not be accurately predicted was FRMPD2, just upstream of the AFP gene. A total of 33 exons were identified but several of those near the 3′ end of the gene are similar to earlier exons and are annotated separately as FRMPD2-like. This may represent a local duplication event but this is difficult to determine as this gene is also poorly annotated in other fish genomes. The other genes contain 1 to 31 exons and encode a variety of proteins with generally unrelated and poorly characterized functions as described in the UniProt database . For example, little is known about UBTD1 other than it contains a single ubiquitin-like domain. The best characterized is RBP3, which is essential for vision as it shuttles retinoids between the retinal pigment epithelium and the photoreceptor cells.
Quantitative PCR indicates that AFP is a single-copy gene in smelt
The similarity in the restriction fragment patterns of the three BACs mentioned above, as well as the presence of only a single copy of the AFP gene on the sequenced BAC, suggested that there is but a single copy of the AFP gene in smelt. To test this, qPCR was performed by amplifying portions of the AFP gene, LDH-A (which is a single copy gene in smelt ), LUC7La (which is present as a single copy near the AFP in BAC O0139C19) and the BAC vector (pBAC-GMR) using the primers listed in Table 1. When DNA from BAC O0139C19 was used as the template, the ΔCT values for each of the serial dilutions were the same for the BAC vector, AFP and LUC7La, indicating that each was present in a 1:1:1 ratio in the BAC clone. Similar results were obtained when smelt genomic DNA was used as template with a 1:1:1 ratio of AFP: LUC7La: LDH-A (Additional file 3: Table S2). Thus the AFP gene of smelt is single copy.
The dearth of polymorphisms in smelt ESTs is consistent with a single-copy AFP gene
BLAST searches were performed on both the non-redundant and EST databases of NCBI. In addition to one AFP gene and two AFP mRNA sequences, 406 ESTs out of almost 37,000 smelt ESTs in the database coded for the AFP. This represents just over 1% of the total. However, only six of those AFP sequences are found in the close to 33,000 ESTs derived from mixed tissues (brain, kidney and spleen ), so 10% of the ~4,000 liver transcripts encode the AFP. The liver libraries were derived either from an unspecified number of individuals  or from six and nine individuals respectively ( and GenBank Accession GR557461, for example)). For Atlantic herring, there are 397 ESTs in the database, derived from the livers of nine individuals . Of these, 34 (~9%) encode the AFP, along with one gene and four mRNA sequences. Transcriptome sequencing of the Pacific herring has generated a large number of additional AFP sequences  and the assembled AFP sequence from liver library was over 99% identical to the Atlantic herring genomic sequence (Genbank accession DQ003023.1).
The smelt AFP gene has had little time to diverge
The homologs used above were also used to determine overall codon usage within each species (data not shown). There were no striking or consistent differences between herring and smelt (data not shown), nor between these two species and fugu (Takifugu rubripes) . Therefore, the codon usage pattern of the AFP gene was of no use in determining the source species.
The AFP gene is absent from conserved syntenic regions of other fishes
Phylogenetic comparisons suggest the neighboring genes to AFP were not transferred
Genes flanking the AFP gene were evaluated for their usefulness in determining the relatedness of the corresponding smelt and herring sequences. RBP3-2 and MMS19 were chosen as they are quite dissimilar to their paralogs. MMS19 was recovered from the Pacific herring liver and testes transcriptome assemblies  but RBP3-2 was absent as it is only expressed in the eye , so it was amplified from Atlantic herring genomic DNA using degenerate primers.
The non-coding regions of syntenic genes and the AFP genes differ in their conservation
The intergenic regions flanking the AFP gene are devoid of viral and transposon sequences
Vista comparisons between the 24.5 kb region between the FRMPD2 and GDF10 genes in smelt, which contains the AFP gene (Figure 1), and the corresponding intergenic regions of the five other fish showed that there are short regions of similarity close to the two flanking genes but no similarities anywhere near the AFP gene (Figure 6D). Therefore, there is no evidence for AFP pseudogenes in these other fish, which does not support the gene loss hypothesis, but is entirely consistent with LGT.
This intergenic region was split into overlapping segments of 2 kb that were used as queries in both BLASTN (nr and EST databases) and BLASTX (nr database) searches and known transposable element or viral sequences were not detected. The sequence was also scanned using Censor  which searches for similarity to repetitive elements catalogued in the repbase database. The 90 putative matches throughout the smelt insert were quite short (<300 bp). Some of these, including the twelve hits within the intergenic regions flanking the AFP gene, may represent residual transposon or retrotransposon sequences, but their level of degradation is extensive compared to the AFP gene. Therefore, there is no evidence to suggest that this LGT was mediated by a transposon or viral vector.
Lateral gene transfer (LGT) is a well-recognized phenomenon in prokaryotes, particularly where antibiotic resistance is acquired by pathogenic strains of bacteria, conferring a selective advantage that has tremendous consequences for human health . Although LGT is far less frequent in eukaryotes, it has been observed in single-celled eukaryotes [49, 50] and the review by Zhaxybayeva and Doolittle provides several salient examples from across the ‘tree’ of life . Despite this, vertical descent is usually taken as the null hypothesis, even in the face of evidence to the contrary, and Keeling and Palmer state that “the erroneous report of substantial HGT in the human genome has probably chilled the field, tarnished more credible claims for eukaryotic HGT” . There may be other examples of lateral AFP gene transfer as surprising similarities have been observed between the AFPs of fungi, diatoms and some bacteria, while similar sequences are absent from most other species throughout these three biological domains . Unicellular eukaryotes, including fish parasites  and plant parasites  have been the recipients of intradomain LGT from other eukaryotes, but given the relative paucity of complete eukaryotic genome sequences, the true extent of such transfers is currently unknown .
In higher organisms, for transferred genes to persist, they must be incorporated into the germ line. So, it is not surprising that extensive LGT has occurred between the intracellular bacterial parasite, Wolbachia, which infects gonadal tissues, and a number of its host insects and nematodes . Reports of transfers to and between multicellular organisms are on the rise , including the previously mentioned transfer of carotenoid biosynthesis genes between fungi and aphids . However, reports of LGT in vertebrates are rare and limited to “selfish” DNAs such as the SPIN transposon that has been transferred between 14 families of lizards .
There are three main criteria of use in establishing LGT in eukaryotes; unexpected phylogenetic relationships, a patchy phyletic distribution and an overlap of life habits , We previously reported the first instance of LGT of an antifreeze gene between two vertebrates. Southern blotting indicated that many intervening species lacked this gene and the sequence of the rainbow smelt AFP differed about as much from that of the Japanese smelt (same family) as it did from the herring (different superorder). We also discussed the possibility of sperm-mediated transformation as a mechanism, whereby sperm absorb foreign DNA fragments, given that these fishes employ external fertilization and have overlapping ranges . A true hybridization, in which an egg was fertilized by the sperm of the other species, is not possible given that herring and smelt belong to different superorders that diverged over 200 million years ago [17–19]. The incompatibility is so severe that even whole nuclei, transplanted into the enucleated cytoplasm of fish eggs from different orders, die early in development .
Here we report several new lines of evidence to support the finding of LGT in fishes. A BAC clone, containing the single copy of the AFP gene found in a smelt, was sequenced (Figure 1). The corresponding microsyntenic region was identified in the five publicly available genome sequences from fish and was found to be devoid of any AFP-like sequences (Figure 4). Not even so much as a pseudogene was found here (Figure 6D), or in any other location within these five genomes. The twelve smelt genes flanking this AFP were unremarkable in that they did not appear unexpectedly similar to their homologs in any other species and their intron sequences were unique (Figure 6A). As four of these fish are more closely related to smelt than herring, [17–19] this suggests that the high conservation of the non-coding sequences of the AFP between herring and smelt is atypical (Figure 6C) . Phylogenetic gene comparisons indicated that it is only the AFP that defies taxonomic conventions (Figure 5, Additional file 5: Figure S2C, [10, 16]), suggesting that the AFP is the only gene within this region that was laterally transferred.
The recent deposition of numerous EST sequences for Atlantic herring and smelt has allowed us to demonstrate the unusual mutation pattern within the AFP gene (Figure 3). All other genes showed similar rates of silent site mutation that were 5- to14-fold higher than that of the AFP gene, while the rate of missense mutation in the AFP suggests the protein sequence is not under strong selection. This is not likely due to codon selection as several of these other genes also appear to be highly expressed as evidenced by their abundance in the EST database. As well, the codon usage of the AFPs does not differ significantly from the other genes.
An unusual feature of the AFP of smelt is that it is encoded by a single-copy gene, whereas the AFPs of herring  the AFGPs of numerous notothenioid fishes , the type I AFPs of flounders  and the type III AFP from ocean pout  are all encoded by multigene families in which most of the copies appear to be clustered, sometimes in tandem repeats. The dearth of polymorphisms in the smelt ESTs is consistent with this difference. Gene amplification is a more efficient method of responding to rapid environmental change, as duplication of a gene encoding a weakly active protein is several orders of magnitude more rapid than the acquisition of point mutations that improve activity. It is thought that gene copy number reduction follows once the genes have improved under selection , and as this has not yet occurred in other AFP-producing fishes, it would be puzzling for the smelt to be an exception. Smelt is also different than most fish in that it does not rely upon AFP alone for freeze resistance but synthesizes high levels (~0.4 M) of the small molecule antifreeze, glycerol. This drops the serum freezing point close to that of seawater, but this is not the case for herring, which has much lower levels of solutes in wintertime [62, 63]. Glycerol is an energetically costly solution to freeze resistance as it is readily lost to the environment, likely through the skin and gills . Although the directionality of the transfer cannot be conclusively determined at this time and the independent evolution of two parallel strategies is not impossible, it seems more plausible that the AFP of smelt was acquired secondarily and is an energetically favorable adjunct to an existing antifreeze strategy.
Phylogenetic comparisons indicate that the type II AFP gene is more similar to fish lectins (<40% identity) than to those of any other group of organisms . A newly acquired AFP gene would confer an immediate selective advantage to fish at risk of freezing in ice-laden waters and as such, it is a likely candidate for LGT since it is not part of a complex of interacting proteins or processes . Since most fish spawn externally in water, their germ line is somewhat vulnerable to DNA uptake. When the evidence is taken together, it all supports the hypothesis that the patchy distribution of highly similar type II AFP genes in herring and smelt is a result of LGT between fishes.
Here we have documented five pieces of new evidence to support the case for the lateral transfer of a lectin-like AFP gene between fish species. i) The arrangement of genes in the region where the AFP gene is located in smelt is shared by other fishes, but only the smelt has the AFP gene. ii) There is also no evidence of AFP pseudogenes in these other fish. iii) Genes flanking the smelt AFP gene are smelt-like and not herring-like, indicating that the sequenced region is normal except for the AFP gene. iv) The introns of genes adjacent to the AFP gene are highly divergent between species, which strongly contrasts with the similarity of the introns of the herring and smelt AFP genes. v) Lastly, the rate of silent vs non-silent mutation of the AFP gene is inconsistent with descent from a common ancestor. As well, smelt is unique in that other species of AFP-producing fishes (including the herring) have multiple copies of their AFP genes.
Screening of rainbow smelt BAC library
The rainbow smelt BAC library filters numbers one to three (CHORI-74, ) were screened with a 40-mer probe (Table 1) designed from GenBank accession DQ004949. The 40-mer was end-labeled with 32PγATP using T4 polynucleotide kinase. Prehybridization, hybridization, imaging and growth of positive clones were performed as in Johnstone et al. (2008 ). A small amount of DNA from each positive clone was amplified using the 40-mer probe and Type II AFP R primer (Table 1) as previously described , except that denaturation was done for 45 s, and the annealing temperature was 65°C.
Isolation and sequencing of BAC clone
BAC DNA was isolated using the Qiagen® Large-Construct Kit (Qiagen, Mississauga, Ontario) according to the manufacturer’s instructions with the addition of a proteinase K digestion (100 μg/mL) for 1 h at 50°C immediately after the nuclease digestion step, plus two phenol/chloroform extractions and one chloroform extraction once the final pellet was resuspended in 1 ml of 10 mM Tris–HCl, 1 mM EDTA (pH 8). The DNA was then extensively dialyzed against the same buffer. The shearing, shotgun cloning and majority of the sequencing, using an Applied Biosystems 3730xl DNA Analyzer and standard ABI dye-terminator technology, were done by the McGill University and Génome Québec Innovation Centre (Montreal, Canada). Problematic regions were sequenced at Robarts Research Institute (London, Canada) using custom primers flanking the region, in combination with the dGTP BigDye terminator v3.0 kit and/or betaine as needed. The sequence has been deposited under GenBank Accession JQ514278.
Assembly and gene annotation
The BAC insert sequence was assembled using the CAP3 program  and genes identified using overlapping 10 kb segments as the query in a variety of BLAST searches  against the non-redundant DNA, protein and EST GenBank databases. Genes were annotated using a wide variety of techniques. This included splice site prediction at the Galaxy website  using SplicePredict trained on spotted pufferfish (Tetraodon nigroviridis) sequences, comparison with ESTs using GeneSeqer  once the ESTs were assembled using the CAP EST Assembler (http://host9.bioinfo3.ifom-ieo-campus.it/cap/), comparisons with protein sequences using Wise2 , as well as a variety of gene prediction programs such as Augustus .
Comparisons with other species
The deduced smelt protein sequences were used as queries for BLAST searches on the Ensembl website (http://www.ensembl.org/index.html) , against the genomes of other fishes and humans; Zv9 Danio rerio (zebrafish) (The Danio rerio Sequencing Project (http://www.sanger.ac.uk/Projects/D_rerio/) Wellcome Trust Sanger Institute), BROADS1 Gasterosteus aculeatus (stickleback) (The Broad Institute at MIT and Harvard), MEDAKA1 Oryzias latipes (medaka) , FUGU4 T. rubripes (fugu) , TETRADON8 T. nigroviridis (spotted pufferfish)  and human GRCh37 . The corresponding syntenic regions from these species were downloaded and viewed using the Apollo Genome Annotation Tool v1.11.6 . Genes were reannotated, when necessary, using the techniques mentioned above. mVISTA  was used for between species comparisons using LAGAN alignments, unless otherwise mentioned.
Determination of dN/dSratios
All available Atlantic herring sequences were downloaded from the non-redundant and EST databases and mitochondrial sequences, short sequences as well as non-protein coding sequences, such as satellite DNA and rDNA were removed. Following contig assembly using DNAMAN v6 (Lynnon Corporation, Pointe-Claire, Canada) and the CAP3 program , homologous smelt sequences were identified using a TBLASTX search. These were categorized as orthologs if both the smelt and herring sequence showed the highest similarity to the same stickleback or zebrafish sequence using TBLASTN and/or BLASTP (http://www.ensembl.org/index.html). Codon usage and percent identity were determined using DNAMAN. Codon aligned orthologous sequence pairs were generated by removing gaps, ambiguous regions including runs of single amino acids with gaps, misidentified introns or likely sequencing errors within poorly represented ESTs. Only alignments longer than 150 bp in length were used to calculate the ratio of non-synonymous changes per non-synonymous site to synonymous substitutions per synonymous site (d N /d S ) using SNAP (Synonymous Non-synonymous Analysis Program) .
Portions of the Altantic herring gene encoding the RBP-3-2 sequence were amplified from genomic DNA by PCR as above but with an annealing temperature of 50°C. Two fully degenerate primers (Table 1, named inner), corresponding to amino acids 414–421 and 442–449 of zebrafish NP_571526.2, were used to obtain a 61 bp segment. Specific primers were then designed and used in combination with degenerate primers corresponding to amino acids 235–242 and 472–479 to obtain two overlapping segments spanning three exons. Fragments were cloned using the TOPO® PCR Cloning Kit (Life Technologies, Burlington, Canada) and sequenced at Robarts Research Institute (London, Canada). MMS19 sequences were obtained as described above, except for the Pacific herring sequence that was deduced from the transcriptome data . Phylogenetic analysis of the aligned sequences was done using MEGA 5.01 . The WAG + I and JTT + G substitution models were the best fit for the RBP-3-2 and MMS-19 alignments respectively. Maximum likelihood analysis with gap deletion was performed using 500 bootstrap replicates.
Quantitative PCR (qPCR)
qPCR reactions were performed with forward and reverse primers for AFP, LDH-A, LUC7La and pBAC-GMR sequences using both smelt genomic DNA and BAC DNA as templates. The primer sequences are listed in Table 1. Five serial dilutions of the DNA were used for each qPCR reaction. The qPCR assay mixture comprised: 12.5 μL of Quanta SYBR, 2 μL of 5 μM forward primer, 2 μL of 5 μM reverse primer and 7.5 μL of nuclease-free water per reaction. The qPCR protocol was: 95°C for 3 min, followed by 40 cycles of 95°C for 15 s, 65°C for 30 s, and 72°C for 15 s. The reactions were recorded and analyzed by the program RQ manager 1.2.1. The efficiencies of qPCR reactions for each set of primers were calculated as: E = (10(−1/slope)-1) × 100 where the slope is derived from the linear regression line of the plot log (DNA concentration) vs. average CT of the 5 serial dilutions for each set of primers using either BAC DNA or genomic DNA as template.
Bacterial artificial chromosome
- (d N /d S ):
Ratio of non-synonymous changes per non-synonymous site to synonymous substitutions per synonymous site
Lateral gene transfer
- C-type lectin:
Calcium-dependent type lectin.
This work was supported by the Genome Canada and Genome BC funded Consortium for Genomic Research on All Salmonids Project (cGRASP) (WSD) and by CIHR (Grant 106612) to PLD, who holds a Canada Research Chair in Protein Engineering.
- DeVries AL, Wohlschlag DE: Freezing resistance in some Antarctic fishes. Science. 1969, 163 (871): 1073-1075. 10.1126/science.163.3871.1073.PubMedView ArticleGoogle Scholar
- Raymond JA, DeVries AL: Adsorption inhibition as a mechanism of freezing resistance in polar fishes. Proc Natl Acad Sci USA. 1977, 74 (6): 2589-2593. 10.1073/pnas.74.6.2589.PubMedPubMed CentralView ArticleGoogle Scholar
- Duman JG: Antifreeze and ice nucleator proteins in terrestrial arthropods. Annu Rev Physiol. 2001, 63: 327-357. 10.1146/annurev.physiol.63.1.327.PubMedView ArticleGoogle Scholar
- Duman JG, Olsen TM: Thermal hysteresis protein activity in bacteria, fungi and phylogenetically diverse plants. Cryobiology. 1993, 30: 322-328. 10.1006/cryo.1993.1031.View ArticleGoogle Scholar
- Fletcher GL, Hew CL, Davies PL: Antifreeze proteins of teleost fishes. Annu Rev Physiol. 2001, 63: 359-390. 10.1146/annurev.physiol.63.1.359.PubMedView ArticleGoogle Scholar
- Hoshino T, Kiriaki M, Nakajima T: Novel thermal hysteresis proteins from low temperature basidiomycete, Coprinus psychromorbidus. Cryo-Letters. 2003, 24 (3): 135-142.PubMedGoogle Scholar
- Raymond JA, Knight CA: Ice binding, recrystallization inhibition, and cryoprotective properties of ice-active substances associated with Antarctic sea ice diatoms. Cryobiology. 2003, 46 (2): 174-181. 10.1016/S0011-2240(03)00023-3.PubMedView ArticleGoogle Scholar
- Gilbert JA, Davies PL, Laybourn-Parry J: A hyperactive, Ca2+-dependent antifreeze protein in an Antarctic bacterium. FEMS Microbiol Letts. 2005, 245 (1): 67-72. 10.1016/j.femsle.2005.02.022.View ArticleGoogle Scholar
- Jia Z, Davies PL: Antifreeze proteins: an unusual receptor-ligand interaction. Trends Biochem Sci. 2002, 27 (2): 101-106. 10.1016/S0968-0004(01)02028-X.PubMedView ArticleGoogle Scholar
- Cheng CH: Evolution of the diverse antifreeze proteins. Curr Opin Genet Dev. 1998, 8 (6): 715-720. 10.1016/S0959-437X(98)80042-7.PubMedView ArticleGoogle Scholar
- Ewart KV, Rubinsky B, Fletcher GL: Structural and functional similarity between fish antifreeze proteins and calcium-dependent lectins. Biochem Biophys Res Commun. 1992, 185 (1): 335-340.PubMedView ArticleGoogle Scholar
- Ewart KV, Fletcher GL: Herring antifreeze protein: primary structure and evidence for a C-type lectin evolutionary origin. Mol Mar Biol Biotechnol. 1993, 2 (1): 20-27.PubMedGoogle Scholar
- Nishimiya Y, Kondo H, Takamichi M, Sugimoto H, Suzuki M, Miura A, Tsuda S: Crystal structure and mutational analysis of Ca2 + −independent type II antifreeze protein from longsnout poacher, Brachyopsis rostratus. J Mol Biol. 2008, 382 (3): 734-746. 10.1016/j.jmb.2008.07.042.PubMedView ArticleGoogle Scholar
- Liu Y, Li Z, Lin Q, Kosinski J, Seetharaman J, Bujnicki JM, Sivaraman J, Hew CL: Structure and evolutionary origin of Ca2 + −dependent herring type II antifreeze protein. PLoS One. 2007, 2 (6): e548-10.1371/journal.pone.0000548.PubMedPubMed CentralView ArticleGoogle Scholar
- Loewen MC, Gronwald W, Sonnichsen FD, Sykes BD, Davies PL: The ice-binding site of sea raven antifreeze protein is distinct from the carbohydrate-binding site of the homologous C-type lectin. Biochemistry. 1998, 37 (51): 17745-17753. 10.1021/bi9820513.PubMedView ArticleGoogle Scholar
- Graham LA, Lougheed SC, Ewart KV, Davies PL: Lateral transfer of a lectin-like antifreeze protein gene in fishes. PLoS One. 2008, 3 (7): e2616-10.1371/journal.pone.0002616.PubMedPubMed CentralView ArticleGoogle Scholar
- Steinke D, Salzburger W, Meyer A: Novel relationships among ten fish model species revealed based on a phylogenomic analysis using ESTs. J Mol Evol. 2006, 62 (6): 772-784. 10.1007/s00239-005-0170-8.PubMedView ArticleGoogle Scholar
- Peng Z, He S, Wang J, Wang W, Diogo R: Mitochondrial molecular clocks and the origin of the major Otocephalan clades (Pisces: Teleostei): A new insight. Gene. 2006, 370: 113-124.PubMedView ArticleGoogle Scholar
- Azuma Y, Kumazawa Y, Miya M, Mabuchi K, Nishida M: Mitogenomic evaluation of the historical biogeography of cichlids toward reliable dating of teleostean divergences. BMC Evol Biol. 2008, 8: 215-10.1186/1471-2148-8-215.PubMedPubMed CentralView ArticleGoogle Scholar
- Yamashita Y, Miura R, Takemoto Y, Tsuda S, Kawahara H, Obata H: Type II antifreeze protein from a mid-latitude freshwater fish, Japanese smelt (Hypomesus nipponensis). Biosci Biotechnol Biochem. 2003, 67 (3): 461-466. 10.1271/bbb.67.461.PubMedView ArticleGoogle Scholar
- Scott GK, Hew CL, Davies PL: Antifreeze protein genes are tandemly linked and clustered in the genome of the winter flounder. Proc Natl Acad Sci USA. 1985, 82 (9): 2613-2617. 10.1073/pnas.82.9.2613.PubMedPubMed CentralView ArticleGoogle Scholar
- Libuda DE, Winston F: Alterations in DNA replication and histone levels promote histone gene amplification in Saccharomyces cerevisiae. Genetics. 2010, 184 (4): 985-997. 10.1534/genetics.109.113662.PubMedPubMed CentralView ArticleGoogle Scholar
- Hendrickson H, Slechta ES, Bergthorsson U, Andersson DI, Roth JR: Amplification-mutagenesis: evidence that "directed" adaptive mutation and general hypermutability result from growth with a selected gene amplification. Proc Natl Acad Sci USA. 2002, 99 (4): 2164-2169. 10.1073/pnas.032680899.PubMedPubMed CentralView ArticleGoogle Scholar
- Bergthorsson U, Andersson DI, Roth JR: Ohno's dilemma: evolution of new genes under continuous selection. Proc Natl Acad Sci USA. 2007, 104 (43): 17004-17009. 10.1073/pnas.0707158104.PubMedPubMed CentralView ArticleGoogle Scholar
- Davies PL, Hough C, Scott GK, Ng N, White BN, Hew CL: Antifreeze protein genes of the winter flounder. J Biol Chem. 1984, 259 (14): 9241-9247.PubMedGoogle Scholar
- Hew CL, Wang NC, Joshi S, Fletcher GL, Scott GK, Hayes PH, Buettner B, Davies PL: Multiple genes provide the basis for antifreeze protein diversity and dosage in the ocean pout, Macrozoarces americanus. J Biol Chem. 1988, 263 (24): 12049-12055.PubMedGoogle Scholar
- Wang X, DeVries AL, Cheng CH: Genomic basis for antifreeze peptide heterogeneity and abundance in an Antarctic eel pout: gene structures and organization. Mol Mar Biol Biotechnol. 1995, 4 (2): 135-147.PubMedGoogle Scholar
- Chen L, DeVries AL, Cheng CH: Convergent evolution of antifreeze glycoproteins in Antarctic notothenioid fish and Arctic cod. Proc Natl Acad Sci USA. 1997, 94 (8): 3817-3822. 10.1073/pnas.94.8.3817.PubMedPubMed CentralView ArticleGoogle Scholar
- Tenover FC: Mechanisms of antimicrobial resistance in bacteria. Am J Infect Control. 2006, 34 (5 Suppl 1): S3-S10. discussion S64-73PubMedView ArticleGoogle Scholar
- Hourston AS, Rosenthal H: Sperm density during active spawning of pacific herring Clupea harengus pallasi. J Fish Res Board Can. 1976, 33: 1788-1790. 10.1139/f76-226.View ArticleGoogle Scholar
- Smith K, Spadafora C: Sperm-mediated gene transfer: applications and implications. Bioessays. 2005, 27 (5): 551-562. 10.1002/bies.20211.PubMedView ArticleGoogle Scholar
- Keeling PJ, Palmer JD: Lateral transfer at the gene and subgenic levels in the evolution of eukaryotic enolase. Proc Natl Acad Sci USA. 2001, 98 (19): 10745-10750. 10.1073/pnas.191337098.PubMedPubMed CentralView ArticleGoogle Scholar
- Moran NA, Jarvik T: Lateral transfer of genes from fungi underlies carotenoid production in aphids. Science. 2010, 328 (5978): 624-627. 10.1126/science.1187113.PubMedView ArticleGoogle Scholar
- Mayer WE, Schuster LN, Bartelmes G, Dieterich C, Sommer RJ: Horizontal gene transfer of microbial cellulases into nematode genomes is associated with functional assimilation and gene turnover. BMC Evol Biol. 2011, 11 (1): 13-10.1186/1471-2148-11-13.PubMedPubMed CentralView ArticleGoogle Scholar
- von Schalburg KR, Leong J, Cooper GA, Robb A, Beetz-Sargent MR, Lieph R, Holt RA, Moore R, Ewart KV, Driedzic WR, et al: Rainbow smelt (Osmerus mordax) genomic library and EST resources. Mar Biotechnol (NY). 2008, 10 (5): 487-491. 10.1007/s10126-008-9089-6.View ArticleGoogle Scholar
- Consortium TU: Ongoing and future developments at the Universal Protein Resource. Nucleic Acids Res. 2011, 39: D214-D219. Database issueView ArticleGoogle Scholar
- Ma X: The characterization of lactate dehydrogenase genes in rainbow smelt (Osmerus mordax). 2009, Canada: Simon Fraser University,Google Scholar
- Richards RC, Achenbach JC, Short CE, Kimball J, Reith ME, Driedzic WR, Ewart KV: Seasonal expressed sequence tags of rainbow smelt (Osmerus mordax) revealed by subtractive hybridization and the identification of two genes up-regulated during winter. Gene. 2008, 424 (1–2): 56-62.PubMedView ArticleGoogle Scholar
- Olsvik PA, Waagbo R, Pedersen SA, Meier S: Transcriptional effects of dietary exposure of oil-contaminated Calanus finmarchicus in Atlantic herring (Clupea harengus). J Toxicol Environ Health A. 2011, 74 (7–9): 508-528.PubMedView ArticleGoogle Scholar
- Roberts SB, Hauser L, Seeb LW, Seeb JE: Development of genomic resources for Pacific Herring through targeted transcriptome pyrosequencing. PLoS One. 2012, 7 (2): e30908-10.1371/journal.pone.0030908.PubMedPubMed CentralView ArticleGoogle Scholar
- Nakamura Y, Gojobori T, Ikemura T: Codon usage tabulated from international DNA sequence databases: status for the year 2000. Nucleic Acids Res. 2000, 28 (1): 292-10.1093/nar/28.1.292.PubMedPubMed CentralView ArticleGoogle Scholar
- van Veen T, Katial A, Shinohara T, Barrett DJ, Wiggert B, Chader GJ, Nickerson JM: Retinal photoreceptor neurons and pinealocytes accumulate mRNA for interphotoreceptor retinoid-binding protein (IRBP). FEBS Lett. 1986, 208 (1): 133-137. 10.1016/0014-5793(86)81547-2.PubMedView ArticleGoogle Scholar
- Miya M, Takeshima H, Endo H, Ishiguro NB, Inoue JG, Mukai T, Satoh TP, Yamaguchi M, Kawaguchi A, Mabuchi K, et al: Major patterns of higher teleostean phylogenies: a new perspective based on 100 complete mitochondrial DNA sequences. Mol Phylogenet Evol. 2003, 26 (1): 121-138. 10.1016/S1055-7903(02)00332-9.PubMedView ArticleGoogle Scholar
- Ishiguro NB, Miya M, Nishida M: Basal euteleostean relationships: a mitogenomic perspective on the phylogenetic reality of the "Protacanthopterygii". Mol Phylogenet Evol. 2003, 27 (3): 476-488. 10.1016/S1055-7903(02)00418-9.PubMedView ArticleGoogle Scholar
- Lavoue S, Miya M, Saitoh K, Ishiguro NB, Nishida M: Phylogenetic relationships among anchovies, sardines, herrings and their relatives (Clupeiformes), inferred from whole mitogenome sequences. Mol Phylogenet Evol. 2007, 43 (3): 1096-1105. 10.1016/j.ympev.2006.09.018.PubMedView ArticleGoogle Scholar
- Yamanoue Y, Miya M, Doi H, Mabuchi K, Sakai H, Nishida M: Multiple invasions into freshwater by pufferfishes (teleostei: tetraodontidae): a mitogenomic perspective. PLoS One. 2011, 6 (2): e17410-10.1371/journal.pone.0017410.PubMedPubMed CentralView ArticleGoogle Scholar
- Kohany O, Gentles AJ, Hankus L, Jurka J: Annotation, submission and screening of repetitive elements in Repbase: RepbaseSubmitter and Censor. BMC Bioinforma. 2006, 7: 474-10.1186/1471-2105-7-474.View ArticleGoogle Scholar
- Barlow M: What antimicrobial resistance has taught us about horizontal gene transfer. Meth Mol Biol. 2009, 532: 397-411. 10.1007/978-1-60327-853-9_23.View ArticleGoogle Scholar
- Andersson JO: Lateral gene transfer in eukaryotes. Cell Mol Life Sci. 2005, 62: 1-16. 10.1007/s00018-004-4437-9.View ArticleGoogle Scholar
- Keeling PJ, Palmer JD: Horizontal gene transfer in eukaryotic evolution. Nat Rev Genet. 2008, 9 (8): 605-618. 10.1038/nrg2386.PubMedView ArticleGoogle Scholar
- Zhaxybayeva O, Doolittle WF: Lateral gene transfer. Curr Biol. 2011, 21 (7): R242-R246. 10.1016/j.cub.2011.01.045.PubMedView ArticleGoogle Scholar
- Raymond JA, Janech MG: Ice-binding proteins from enoki and shiitake mushrooms. Cryobiology. 2009, 58 (2): 151-156. 10.1016/j.cryobiol.2008.11.009.PubMedView ArticleGoogle Scholar
- Andersson JO, Sjogren AM, Horner DS, Murphy CA, Dyal PL, Svard SG, Logsdon JM, Ragan MA, Hirt RP, Roger AJ: A genomic survey of the fish parasite Spironucleus salmonicida indicates genomic plasticity among diplomonads and significant lateral gene transfer in eukaryote genome evolution. BMC Genomics. 2007, 8: 51-10.1186/1471-2164-8-51.PubMedPubMed CentralView ArticleGoogle Scholar
- Richards TA, Soanes DM, Jones MD, Vasieva O, Leonard G, Paszkiewicz K, Foster PG, Hall N, Talbot NJ: Horizontal gene transfer facilitated the evolution of plant parasitic mechanisms in the oomycetes. Proc Natl Acad Sci USA. 2011, 108 (37): 15258-15263. 10.1073/pnas.1105100108.PubMedPubMed CentralView ArticleGoogle Scholar
- Andersson JO: Gene transfer and diversification of microbial eukaryotes. Annu Rev Microbiol. 2009, 63: 177-193. 10.1146/annurev.micro.091208.073203.PubMedView ArticleGoogle Scholar
- Dunning Hotopp JC, Clark ME, Oliveira DC, Foster JM, Fischer P, Munoz Torres MC, Giebel JD, Kumar N, Ishmael N, Wang S, et al: Widespread lateral gene transfer from intracellular bacteria to multicellular eukaryotes. Science. 2007, 317 (5845): 1753-1756. 10.1126/science.1142490.PubMedView ArticleGoogle Scholar
- Sinkovics JG: Horizontal gene transfers with or without cell fusions in all categories of the living matter. Adv Exp Med Biol. 2011, 714: 5-89.PubMedView ArticleGoogle Scholar
- Gilbert C, Hernandez SS, Flores-Benabib J, Smith EN, Feschotte C: Rampant Horizontal Transfer of SPIN Transposons in Squamate Reptiles. Mol Biol Evol. 2011, 29 (2): 503-515.PubMedPubMed CentralView ArticleGoogle Scholar
- Yan SY, Tu M, Yang HY, Mao ZG, Zhao ZY, Fu LJ, Li GS, Huang GP, Li SH, Jin GQ, et al: Developmental incompatibility between cell nucleus and cytoplasm as revealed by nuclear transplantation experiments in teleost of different families and orders. Int J Dev Biol. 1990, 34 (2): 255-266.PubMedGoogle Scholar
- Cheng CH, Detrich HW: Molecular ecophysiology of Antarctic notothenioid fishes. Philos Trans R Soc Lond B Biol Sci. 2007, 362 (1488): 2215-2232. 10.1098/rstb.2006.1946.PubMedPubMed CentralView ArticleGoogle Scholar
- Davies PL: Conservation of antifreeze protein-encoding genes in tandem repeats. Gene. 1992, 112 (2): 163-170. 10.1016/0378-1119(92)90372-V.PubMedView ArticleGoogle Scholar
- Raymond JA: Glycerol is a colligative antifreeze in some northern fishes. J Exp Zool. 1992, 262 (3): 347-352. 10.1002/jez.1402620316.View ArticleGoogle Scholar
- Driedzic WR, Ewart KV: Control of glycerol production by rainbow smelt (Osmerus mordax) to provide freeze resistance and allow foraging at low winter temperatures. Comp Biochem Physiol B Biochem Mol Biol. 2004, 139 (3): 347-357. 10.1016/j.cbpc.2004.04.007.PubMedView ArticleGoogle Scholar
- Jain R, Rivera MC, Lake JA: Horizontal gene transfer among genomes: the complexity hypothesis. Proc Natl Acad Sci USA. 1999, 96 (7): 3801-3806. 10.1073/pnas.96.7.3801.PubMedPubMed CentralView ArticleGoogle Scholar
- Johnstone KA, Lubieniecki KP, Chow W, Phillips RB, Koop BF, Davidson WS: Genomic organization and characterization of two vomeronasal 1 receptor-like genes (ora1 and ora2) in Atlantic salmon Salmo salar. Mar Genomics. 2008, 1 (1): 23-31. 10.1016/j.margen.2008.04.003.PubMedView ArticleGoogle Scholar
- Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9 (9): 868-877. 10.1101/gr.9.9.868.PubMedPubMed CentralView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.PubMedView ArticleGoogle Scholar
- Blankenberg D, Von Kuster G, Coraor N, Ananda G, Lazarus R, Mangan M, Nekrutenko A, Taylor J: Galaxy: a web-based genome analysis tool for experimentalists. Curr Protoc Mol Biol. 2010, Chapter 19 (Unit 19 10): 11-21.Google Scholar
- Usuka J, Zhu W, Brendel V: Optimal spliced alignment of homologous cDNA to a genomic DNA template. Bioinformatics. 2000, 16 (3): 203-211. 10.1093/bioinformatics/16.3.203.PubMedView ArticleGoogle Scholar
- Birney E, Clamp M, Durbin R: GeneWise and Genomewise. Genome Res. 2004, 14 (5): 988-995. 10.1101/gr.1865504.PubMedPubMed CentralView ArticleGoogle Scholar
- Stanke M, Steinkamp R, Waack S, Morgenstern B: AUGUSTUS: a web server for gene finding in eukaryotes. Nucleic Acids Res. 2004, 32: W309-W312. 10.1093/nar/gkh379. (Web Server issue)PubMedPubMed CentralView ArticleGoogle Scholar
- Flicek P, Amode MR, Barrell D, Beal K, Brent S, Chen Y, Clapham P, Coates G, Fairley S, Fitzgerald S, et al: Ensembl 2011. Nucleic Acids Res. 2011, 39: D800-D806. 10.1093/nar/gkq1064. (Database issue)PubMedPubMed CentralView ArticleGoogle Scholar
- Kasahara M, Naruse K, Sasaki S, Nakatani Y, Qu W, Ahsan B, Yamada T, Nagayasu Y, Doi K, Kasai Y, et al: The medaka draft genome and insights into vertebrate genome evolution. Nature. 2007, 447 (7145): 714-719. 10.1038/nature05846.PubMedView ArticleGoogle Scholar
- Aparicio S, Chapman J, Stupka E, Putnam N, Chia JM, Dehal P, Christoffels A, Rash S, Hoon S, Smit A, et al: Whole-genome shotgun assembly and analysis of the genome of Fugu rubripes. Science. 2002, 297 (5585): 1301-1310. 10.1126/science.1072104.PubMedView ArticleGoogle Scholar
- Jaillon O, Aury JM, Brunet F, Petit JL, Stange-Thomann N, Mauceli E, Bouneau L, Fischer C, Ozouf-Costaz C, Bernot A, et al: Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype. Nature. 2004, 431 (7011): 946-957. 10.1038/nature03025.PubMedView ArticleGoogle Scholar
- Consortium IHGS: Finishing the euchromatic sequence of the human genome. Nature. 2004, 431 (7011): 931-945. 10.1038/nature03001.View ArticleGoogle Scholar
- Lewis SE, Searle SM, Harris N, Gibson M, Lyer V, Richter J, Wiel C, Bayraktaroglir L, Birney E, Crosby MA, et al: Apollo: a sequence annotation editor. Genome Biol. 2002, 3 (12): RESEARCH0082-PubMedPubMed CentralView ArticleGoogle Scholar
- Frazer KA, Pachter L, Poliakov A, Rubin EM, Dubchak I: VISTA: computational tools for comparative genomics. Nucleic Acids Res. 2004, 32: W273-W279. 10.1093/nar/gkh458. Web Server issuePubMedPubMed CentralView ArticleGoogle Scholar
- Korber B: HIV sequence signatures and similarities. Computational and evolutionary analysis of HIV molecular sequences. Edited by: Rodrigo AG, Learn GH. 2000, Dordrecht, Netherlands: Kluwer Academic Publishers, 55-72.Google Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol. 2011, 10: 2731-2739.View ArticleGoogle Scholar