An evolutionary perspective on Elovl5 fatty acid elongase: comparison of Northern pike and duplicated paralogs from Atlantic salmon
© Carmona-Antoñanzas et al.; licensee BioMed Central Ltd. 2013
Received: 23 November 2012
Accepted: 11 April 2013
Published: 19 April 2013
The ability to produce physiologically critical LC-PUFA from dietary fatty acids differs greatly among teleost species, and is dependent on the possession and expression of fatty acyl desaturase and elongase genes. Atlantic salmon, as a result of a recently duplicated genome, have more of these enzymes than other fish. Recent phylogenetic studies show that Northern pike represents the closest extant relative of the preduplicated ancestral salmonid. Here we characterise a pike fatty acyl elongase, elovl5, and compare it to Atlantic salmon elovl5a and elovl5b duplicates.
Phylogenetic analyses show that Atlantic salmon paralogs are evolving symmetrically, and they have been retained in the genome by purifying selection. Heterologous expression in yeast showed that Northern pike Elovl5 activity is indistinguishable from that of the salmon paralogs, efficiently elongating C18 and C20 substrates. However, in contrast to salmon, pike elovl5 was predominantly expressed in brain with negligible expression in liver and intestine.
We suggest that the predominant expression of Elovl5b in salmon liver and Elovl5a in salmon intestine is an adaptation, enabled by genome duplication, to a diet rich in terrestrial invertebrates which are relatively poor in LC-PUFA. Pike have retained an ancestral expression profile which supports the maintenance of PUFA in the brain but, due to a highly piscivorous LC-PUFA-rich diet, is not required in liver and intestine. Thus, the characterisation of elovl5 in Northern pike provides insights into the evolutionary divergence of duplicated genes, and the ecological adaptations of salmonids which have enabled colonisation of nutrient poor freshwaters.
KeywordsAtlantic salmon Elongase of very long-chain fatty acids Northern pike Paralogous genes Whole-genome duplication
Atlantic salmon (Salmo salar) have been the focus of considerable research effort as a result of their widespread environmental and economic importance as a sporting and cultured species. In addition, in common with all other salmonids, they possess a comparatively recently duplicated genome, believed to have arisen as a result of a relatively recent autotetraploidisation event between 25 and 100 mya [1, 2]. Whole-genome duplication (WGD) has been argued as a powerful evolutionary force creating new raw material for evolution to act upon , thus enabling the divergence and neo- or subfunctionalisation of duplicated loci promoting adaptation and speciation. The imprints of three or four ancient duplications can be detected in vertebrate genomes, including a specific event early in teleost evolution and the recent one in salmonids [4, 5]. Esocids (members of the pike family) are regarded as having the closest extant preduplicated (diploid) genomes to salmonids, based on molecular phylogenetic studies [6, 7], karyotype data , and comparative analyses of expressed gene sequences [9, 10]. Therefore, Northern pike (Esox lucius) is representative of a sister-group to salmonids, and can be viewed as an appropriate species to study the consequences of genome duplication in salmonids. Despite their shared ancestry and overlapping habitats, Atlantic salmon and Northern pike have differing life histories and feeding behaviours where, in freshwaters, pike have a largely piscivorous diet, and salmon a diet rich in terrestrial insects . These differences may be reflected in differing nutritional physiology and, in particular, lipid biochemistry.
Teleosts, like all vertebrates, are unable to synthesise polyunsaturated fatty acids (PUFA) de novo, and so they are essential and required in the diet [12, 13]. However, which PUFA can satisfy the dietary requirement for essential fatty acids (EFA) varies with species. The long-chain PUFA (LC-PUFA), arachidonic acid (20:4n-6, ARA), eicosapentaenoic acid (20:5n-3, EPA), and docosahexaenoic acid (22:6n-3, DHA), which have essential functions in vertebrate immune defense systems and neuronal membranes [14–16] can be produced endogenously, in some but not all vertebrates, from the base EFA, α-linolenic acid (18:3n-3, ALA) and linoleic acid (18:2n-6, LA), by desaturation and elongation . The capability to produce LC-PUFA from EFA varies between fish species [13, 18], and salmonids, including Atlantic salmon, brown trout (Salmo trutta), and Arctic charr (Salvelinus alpinus) have substantially higher LC-PUFA biosynthetic efficiency in comparison with other freshwater species, including zebrafish (Danio rerio), Nile tilapia (Oreochromis niloticus) [11, 19, 20], and Northern Pike [21–24]. The ability to produce LC-PUFA is dependent on the possession and expression of fatty acyl desaturase (fad) and fatty acid elongase (elovl) genes, and salmonids, in contrast to many other fish species examined, have a complete set of genes and expressed enzymes required for the production of ARA from LA, and EPA and DHA from ALA [25, 26]. In Atlantic salmon some of these LC-PUFA biosynthetic enzymes appear to have arisen from duplicated genes, and the subsequent neo- or subfunctionalisation has been hypothesised as an enabling adaptation for salmonids to thrive in relatively nutrient-poor freshwater environments .
The aim of the present study was to characterise a critical gene and enzyme of LC-PUFA biosynthesis, Elovl5, in Northern pike and to compare with previously identified, duplicated elovl5 paralogs in Atlantic salmon. Elovl5 catalyses the first and second elongations of LA and ALA and is therefore essential for the production of LC-PUFA. Thus, comparison between the sequence, activity and expression of pike and salmon Elovl5 genes may provide insights into mechanisms that have driven the evolution and ecological adaptations of salmonids.
Northern pike Elovl5 sequence and phylogenetics
Purifying selection on salmonid Elovl5 paralogs
Functional characterisation of pike Elovl5
Functional characterisation of Northern pike Elovl5 elongase, and its role in the biosynthesis of long-chain polyunsaturated fatty acids (LC-PUFA)
C18 → 20
C20 → 22
C18 → 20
C20 → 22
C20 → 22
C22 → 24
C22 → 24
C18 → 20
C20 → 22
C18 → 20
C20 → 22
C20 → 22
C22 → 24
C22 → 24
Tissue distribution of pike Elovl5
Segregation analysis and mapping of salmon elovl5duplicated loci
The amplicons derived from the elovl5-linked microsatellite primer sets resolved clearly and consistently. In both cases allelic size variants, consistent with amplification of a single discrete locus, were detected. Among the pedigree panel parental samples, 12 different alleles were observed for MS_ elovl5a and four alleles for MS_elovl5b. Joint segregation analysis of the two panels (Br5 and Br5/2), informative for sire based linkage only, did not detect linkage between the two loci (Ho = independent assortment; P = 0.78 and 0.11 for Br5 and Br5/2, respectively). Genetic mapping of the two elovl5 loci in the SalMap family (Br5) confirmed the Atlantic salmon paralogs to be located in distinct linkage groups (LG): elovl5b on LG 5, and elovl5a on LG 33 (LOD scores > 3.5).
The primary aim of this study was to characterise Northern pike Elovl5, a critical enzyme of LC-PUFA biosynthesis in vertebrates. The precise reasons for undertaking the work were ultimately to gain understanding of the evolutionary and ecological adaptations of salmonids. Phylogenetic evidence indicates that esocids are the nearest living relatives of salmonids, having diverged at some point prior to a WGD event in the common ancestor of all salmonids between 25 and 100 mya [1, 2]. As WGD has been widely suggested as a major enabling event in evolutionary innovation , comparison of single preduplicated genes in pike with their duplicated paralogs in Atlantic salmon has the potential to shed light on evolutionary mechanisms and adaptation in salmonids. The genes of the LC-PUFA biosynthetic pathway are interesting candidates for studies of this type because both genetic and biochemical evidence suggests that salmonids have a higher LC-PUFA biosynthetic capacity than many other fish species [11, 20]. Atlantic salmon have duplicated genes for both desaturases and elongases of fatty acids and, in the case of desaturases, duplicates appear to have diverged and neo- or subfunctionalised to provide enzymatic activities for the entire LC-PUFA pathway [26, 27]. In contrast, other fish species, particularly marine carnivorous species, are incapable of endogenous production of significant amounts of LC-PUFA because they lack critical genes of the biosynthetic pathway [17, 25, 35]. It has been suggested that this enhancement of LC-PUFA capacity in salmonids has been a factor in their success in colonising relatively nutrient poor freshwater environments . Salmonids can thrive, and are often the only fish group, in flowing bodies of water, which are reliant on allochthonous terrestrial inputs as a major source of energy. Terrestrial inputs whether directly from leaf litter, or from insect species are poor sources of LC-PUFA, especially DHA, in contrast with food sources in marine or eutrophic freshwater environments which are underpinned by blooms of phytoplankton rich in LC-PUFA [17, 20, 35, 36]. Despite the fact that the tissues of freshwater and marine fish are generally rich in C20 and C22 LC-PUFA, the strategies they utilise to fulfil such requirements vary depending on the species and the dietary input. Northern pike is a strictly freshwater species, whose distribution overlaps that of salmonids, and which shares a relatively recent common ancestor with salmonids. However, pike differ from salmonids in exhibiting a far more piscivorous feeding behaviour. Hence, our particular interest in LC-PUFA biosynthetic genes in this species in comparison to those in Atlantic salmon.
Genetic linkage analyses established that the Atlantic salmon elovl5 duplicates are located on different linkage groups: elovl5a on LG 33, and elovl5b on LG 5. Cytogenetic mapping using fluorescence in situ hybridisation has assigned salmon LG 33 and LG 5 to chromosomes 28 and 13, respectively [37, 38], both of which are acrocentric. While our data clearly show that the two elovl5 loci are not physically linked, in a recent analysis  the Atlantic salmon chromosomes ssa28 and ssa13 were not homeologous and it is therefore not possible to conclude that the salmon elovl5 paralogous genes are the result of a WGD. It is also possible that this duplication is unique to Atlantic salmon, since no clear evidence of duplicated elovl5 genes in other salmonids exist in the current sequence databases. However, compared to all other salmonids which have c. 100 chromosome arms, Atlantic salmon is unique in possessing only 72–74 chromosome arms, believed to be the result of species-specific tandem fusions and other rearrangements . Furthermore, linkage maps show  that salmon chromosome ssa13 has homeologous regions on at least three other salmon chromosomes, and shares syntenic blocks with at least four separate chromosomes from the diploid stickleback (Gasterosteus aculeatus). Thus, given the complexity of the Atlantic salmon genome, it would be premature to reject a WGD origin for the two salmon elovl5 loci.
Phylogenetic analysis confirmed the basal nature of Northern pike within the protacanthopterigyans. Analysis of the rates of nucleotide substitution in the Atlantic salmon Elovl5 paralogous genes indicated that they are under an evolutionary regime of purifying selective pressure (ω < 1), and are currently evolving at comparable evolutionary rates at the protein level, thus supporting the idea that both duplicates are physiologically required, and have the same biochemical activity as the pike Elovl5. In a larger phylogenetic study,  performed pairwise dN/dS analyses on 408 sets of duplicated salmon genes using a preduplicated set of Northern pike orthologs as outgroups, and similarly concluded that salmon paralogs were predominantly exposed to purifying selection, although some loci may be showing a relaxation of selective pressure suggesting that evolution was acting asymmetrically on some paralogs due to reduced constraints. A closer inspection of the dN/dS along the coding sequences of vertebrate elovl5 suggested the accommodation of a major catalytic site where stronger functional constraints seemed to have acted against the retention of nonsynonymous mutations that would compromise the performance of the enzyme, or impair its activity [40, 41].
Here we also tested for functional similarity of pike Elovl5 to salmon paralogs by heterologous expression of the pike enzyme in yeast. Supporting the phylogenetic results, the activity of the pike Elovl5 was indistinguishable from previous assays of the salmon paralogs . Pike Elovl5 was able to lengthen PUFA substrates with chain lengths from C18 to C22. The 18:4n-3, 18:3n-6 and C20 specificity of pike Elovl5 and both salmon Elovl5 paralogs was very similar, and similar to that in other vertebrates . Only very low, residual activity for the production of 24:5n-3 was detected in yeast transformed with the pike elongase when incubated with 22:5n-3, as previously observed in other species, including salmonids [25, 32, 42]. As 24:5n-3 is an important intermediate in most vertebrates for the biosynthesis of DHA, and in salmonids C22 to C24 activity has been demonstrated in Elovl2 and Elovl4 [25, 28], it would be interesting to look for and study these elongase genes in pike.
Until recently the production of LC-PUFA from C18 PUFA was thought to proceed via an alternating desaturation/elongation cycle with the initial step being Δ6 desaturation. However, recent evidence suggests that an alternative pathway is also possible, with the initial step being C18 to C20 elongation, based upon the ability of Δ6 desaturases to also catalyse Δ8 desaturation. For such an alternative loop to exist, elongase activity on LA and ALA is required, and this has been demonstrated in other species [42–44] and, according to results presented here, is also shown by pike Elovl5. Although not yet demonstrated in salmon, we would expect that, given the Δ8 activity of salmon Δ6 desaturases , salmon Elovl5 enzymes would also possess this activity.
Although the phylogenetic and functional analyses point towards the maintenance of ancestral activities, the expression profiles indicate functional partitioning in the salmon elovl5 paralogs. Salmon Elovl5a and Elovl5b have different tissue expression profiles, with Elovl5a being expressed at higher level in intestine and Elovl5b in liver . In addition, the nutritional regulation of mRNA transcription of these genes differs in tissues from fish fed diets containing low levels LC-PUFA . The tissue distribution of pike Elovl5 transcripts showed that the highest expression across the tissues tested was in brain, whereas most other tissues, including liver and intestine, showed very low expression in comparison. In contrast, salmon Elovl5 transcripts were predominantly expressed in liver and intestine, with much lower expression in brain, and even less in other tissues . The pattern of pike Elovl5 tissue expression closely resembles the pattern of LC-PUFA biosynthetic gene expression in carnivorous marine fish, where expression and activity is low in liver and intestine, but high in the brain . Brain, as with all neural tissues, has a high LC-PUFA level [17, 46] but, as yet, mechanisms for the biosynthesis, or the transport of fatty acids to the brain are not fully understood . Mammalian brain is only capable of biosynthesising a restricted set of fatty acids  and it is clear that fatty acid uptake in the brain is different to uptake in most other tissues, probably due the requirement for passage across the blood–brain barrier. Radiolabeled PUFA injected intraperitoneally in Northern pike could be detected in high concentrations throughout the body, with the exception of brain, which consistently contained the lowest amounts of injected PUFA . Currently, most studies support an energy-dependant mechanism that facilitates and regulates fatty acid uptake influenced by chain length and degree of unsaturation [49, 50]. The gene expression results from pike and other carnivorous fish suggest that brain may have endogenous biosynthetic machinery for LC-PUFA to supply the high requirements of this tissue, whereas low expression in liver and intestine indicate a reduced requirement for LC-PUFA, or at least DHA, in these tissues. Although low levels of LC-PUFA biosynthesis have been detected in isolated pike hepatocytes , experiments in vivo failed to find evidence of significant conversion in liver . In salmon, although Elovl5 and other LC-PUFA biosynthetic genes are expressed in brain, the highest expression levels are in liver and intestine, the main tissues responsible for dietary fatty acid uptake, biosynthesis and distribution [17, 25]. Salmonid liver has a comparatively high capacity for biosynthesising LC-PUFA such as EPA and DHA , whereas marine carnivorous fish, such as sea bass (Dicentrarchus labrax) have negligible capacity depending on dietary LC-PUFA .
This above discussion should be qualified by noting that studies on various salmonids and other freshwater species fed with artificial diets with varying LC-PUFA contents and compositions have shown that hepatic desaturase and elongase enzymes exhibited higher expression when the amount of dietary LC-PUFA decreased and LA and/or ALA increased [11, 20, 25, 51]. Although, it is possible that LC-PUFA biosynthetic enzymes in liver are also under nutritional regulation in pike, the fish in the present study were obtained from wild stock and thus essentially piscivorous and so would be equivalent to salmon fed diets containing fish oil (i.e. high in LC-PUFA) .
Northern pike possess an active Elovl5 gene, which is homologous to the duplicated salmon elovl5a and elovl5b genes. There is no evidence of positive (diversifying) selection acting on the salmon genes, or on the pike gene. On the contrary there is evidence of purifying selection maintaining the activity of all of these genes under functional constraints at the protein level. Supporting this, the enzymatic activities of the pike and the two salmon enzymes are indistinguishable. However, a sharply contrasting expression profile was noted between pike Elovl5 and salmon Elovl5 paralogs. Pike elovl5 expression was restricted to brain, in common with other carnivorous fish, whereas salmon have highest expression in liver and intestine. This may be the result of adaptations; salmon to a diet in freshwater relatively poor in LC-PUFA, while pike have a highly piscivorous diet containing higher levels of LC-PUFA. Future studies on the possible diversification and/or neo-or subfunctionalisation of the duplicated salmon genes should focus on mechanisms of gene expression and transcriptional regulation.
Wild Northern pike were caught in Airthrey Loch, Stirling (Scotland). The individuals were euthanised humanely, immediately and rapidly dissected and the tissue samples stored in RNAlater® (Invitrogen, Paisley, UK) at room temperature until tissue homogenisation was performed 24 h later. Fin clips were obtained and stored in 70% ethanol for genomic DNA extraction. Pike were obtained conforming to local ethical regulations regarding field studies on wild animals, and with the permission of the local fisheries managers.
Molecular cloning of pike Elovl5 cDNA
Details of primer pairs (restriction sites for Xba I and Bam HI, and fluorescent CAG and M13R binding sites underlined) used for the cloning of pike Elovl5 ORF in pYES2, qPCR analysis of tissue expression and for the genotyping of Atlantic salmon paralogous genes
elovl5 cDNA cloning
5’- ATGGATGGGTCCCAGAGA -3’
5’- AGTTCATAACGAACCACCAGAT -3’
5’- TGGGGCCCAAGTACATGA -3’
5’- TGGACGAAGCTGTTAAGGG -3’
5’- GATGGCAGAGAGCCCATAGT -3’
5’- CCACACAGCAGACACCATCT -3’
5’- AGATGGTGTCTGCTGTGTGG -3’
5’- ATGCTCAACATCTGGTGGTTT -3’
5’- GCCCAGGTTCGCATCACCCAG -3’
5’- ATTCCGGGGGTCATTTGAGATAGACG -3’
5’- CCCGAGCTCGGATCCAAATGGAGACTTTTAATCAGAGACTTAACACC -3’
5’- GGGTCTAGACTCGAGCTTCAGTCCGCCCTCACTTTCCT -3’
5’- TTCGAATGTCTGCCCTATCAAC -3’
5’- CCTTCCTTGGATGTGGTAGC -3’
5’- CCTTTGCACTGACCGTGATA -3’
5’- GCGTGTCCTGGCAGTAGAA -3’
5’- AAGATCGACCATCGTTCTGG -3’
5’- CTGGCAGCCTTCTTATCGAC -3’
5’- ACAATTGCCATTTTTGCAGATAGC -3’
5’- GGATAACAATTTCACACAGGAGCCATTCTTGATCCGCTTAT -3’
5’- CAGTCGGGCGTCATCAAGCCCGATATGATATTACCGTATT -3’
5’- GTAAAATGGTGACTTTGGGTTCAG -3’
Sequence and phylogenetic analyses
The amino acid (AA) sequence deduced from pike elovl5 cDNA was aligned with salmonid orthologs, including Atlantic salmon, rainbow trout (Oncorhynchus mykiss, [GenBank:NM_001124636.1]), and masu salmon (Oncorhynchus masou, [GenBank:DQ067616.1]) using the BLOSUM62 matrix (BioEdit 7.1.3, Tom Hall, Ibis Biosciences, Abbott Laboratories). Further global pairwise alignments of corresponding DNA sequences and those of other vertebrate elovl5 genes were performed using the Needleman-Wunsch algorithm [http://www.ebi.ac.uk/Tools/psa/emboss_needle/]. A phylogenetic tree was constructed on the basis of the nucleotide sequence alignment between pike elovl5 coding sequence, and other vertebrate elovl5-like ORF in order to take account of single nucleotide substitutions. The 16 vertebrate elovl5 sequences used to construct the phylogenetic tree constitute the species that were considered for further analyses. The tree was rooted to Mortierella alpina PUFA elongase [GenBank:AF206662] as the outgroup species applying the Neighbour Joining algorithm , and the maximum composite likelihood (ML) model using MEGA4 . The bootstrap test was conducted by burning 10,000 trees. Northern pike Elovl5 membrane topology was predicted using PredictProtein tool [http://www.predictprotein.org], which applies the hydrophobicity PHD method (HeiDelberg) with an expected average accuracy of > 98%.
Tests of selection
The influence of selection on Elovl5 genes in pike and salmon was assessed by determining the rates of nonsynonymous substitutions (dN) versus synonymous substitutions (dS) per site in the coding sequence of pike Elovl5, and vertebrate orthologs. For a nucleotide change in a single codon, a synonymous substitution does not result in a change of the amino acid that is specified, whereas a nonsynonymous substitution does. The “acceptance rate” (ω = dN/dS) termed by  is related to evolutionary constraints at the protein level. A value of ω > 1 indicates evidence of Darwinian positive selection whereas ω < 1 suggests negative (purifying) selection, and ω = 1 neutral selection .
For greater confidence, a variety of stochastic algorithms were used to assess the dN/dS ratio based on codon-aligned nucleotide sequences. SNAP  employs the Nei and Gojobori method  incorporating a statistic that allows a more accurate estimation of the variances and covariances of dS and dN . This method allows the use of a best estimate of an ancestral sequence to perform pairwise comparisons between the predicted ancestor and the rest of the sequences, providing an overall ω for the sequences analysed. Similarly, Datamonkey [http://www.datamonkey.org/] is a publically available set of tools implemented by HyPhy package that include robust maximum likelihood approaches to determine the total species dN/dS ratio. The GA-branch analysis (HKY85 model)  is based on a genetic algorithm approach that assigns rate classes of ω to lineages in a phylogenetic context. Finally, Tajima’s test implemented by MEGA4 was used to determine the relative rate of nucleotide and amino acid substitution of the salmon Elovl5 paralogs over evolutionary time using pike Elovl5 as an outgroup. This test performs a chi-square (χ2) to determine whether the rates of evolution of the two compared sequences are comparable (Ho), or not (HA), thus disrupting the rate of the molecular clock .
Substrate specificity of pike Elovl5
The pike Elovl5 coding sequence was obtained by PCR from full-length cDNA with primers containing restriction sites, PIKE5ORF_F (Bam HI) and PIKE5ORF_R (Xho I) (Table 2). This fragment and the yeast episomal plasmid pYES2 (Invitrogen) were digested with the corresponding restriction endonucleases (New England BioLabs, Herts, UK), and ligated using T4 DNA Ligase (Bioline, London, UK). S. cerevisiae (strain INVSc1) transformation, selection and growth were as described in . The substrate specificity of pike Elovl5 was determined by incubating transformed yeast with one of the following fatty acid substrates: ALA (18:3n-3), LA (18:2n-6), stearidonic acid (18:4n-3, SDA), γ-linolenic acid (18:3n-6, GLA), EPA (20:5n-3), and ARA (20:4n-6). The fatty acid substrates (> 99% pure) and chemicals used to prepare the S. cerevisiae minimal medium-uracil were purchased from Sigma Chemical Co. Ltd. Yeast transformed with pYES2 containing no insert were cultured under the same conditions described above and used as control treatment. Yeast fatty acid analysis was as described in .
Tissue distribution of pike Elovl5
Tissue expression of Elovl5 was determined by quantitative real-time PCR (qPCR). Total RNA was extracted from brain, intestine, liver, adipose tissue, white muscle, kidney, spleen, heart and gill from juvenile Northern pike (n = 3) using organic solvent Tri Reagent (Sigma). RNA quality and quantity were assessed by electrophoresis and spectrophotometry (Nanodrop 1000, Thermo Scientific, Wilmington, USA). One microgram of DNase (DNA-free kit, Ambion, Applied Biosystems, Warrington, UK) treated total RNA was reverse transcribed into cDNA using SuperScriptTM III Reverse Transcriptase (Invitrogen) and primed with 5 μM of oligo(dT) and 15 μM of random hexamers (AB Applied Biosystems). Following primer annealing at 25°C for 5 min, and cDNA synthesis at 55°C for 1 h, reactions were stopped by heating at 70°C for 15 min, and cDNA was diluted 5-fold with nuclease-free water. The elovl5 qPCR primers were designed to anneal to each of two predicted exons spanning an intron, whereas the ef1-α1 (reference gene) primer pair was designed in a region corresponding to the coding sequence. Table 2 shows the sequence of the primers used, their specific annealing temperatures, the size of fragments produced, and the reference sequences used for primer design, using the Primer3 software . qPCR analyses were performed using a Mastercycler RealPlex2 (Eppendorf, Stevenage, UK) in a final volume of 20 μl containing 5 μl of diluted cDNA (1/20 for ef1-α1 and elovl5, and 1/5,000 for 18S rRNA), 0.5 μM of each primer, and 10 μl of SensiMix SYBR No-ROX (Bioline). Amplifications were carried out with a systematic negative control (NTC, no template control) containing no cDNA. Standard curves with dilutions: 1/5, 1/10, 1/20, 1/50, 1/100, 1/200, and 1/500 were prepared for ef1-α1 and elovl5, whereas dilutions ranging from 1/10 to 1/1,000,000 were made up for 18S rRNA given its predominant abundance in the total RNA. Thermal cycling was initiated at 95°C for 10 min, followed by 40 cycles: 15 s denaturing step at 95°C, 15 s at the specific primer pair annealing Tm (Table 2), and 15 s extension at 72°C. After the amplification phase, a melting curve from 60°C to 95°C was performed, enabling confirmation of amplification of a single product in each reaction. The qPCR product sizes were determined by agarose gel electrophoresis and their identity confirmed by sequencing. No primer-dimer formation occurred in the NTC.
Normalised expression values were generated by the ΔCt method  using the geometric mean of the Ct values of the two reference genes (ef1-α1 and 18S rRNA). Results were expressed as the percentage of normalised expression relative to the sum of expression across all tissues tested. Normalised expression values of salmon elovl5a and elovl5b reported in  were treated in the same way. Differences in the expression of pike Elovl5 among different tissues were determined by one-way analysis of variance (ANOVA) followed by Tukey’s multiple comparison test at the significance level of P ≤ 0.05 (PASWS 18.0, SPSS Inc., Chicago, USA).
Linkage analysis of duplicated loci in Atlantic salmon
BlastN searches of elovl5a and elovl5b cDNA sequences ([GenBank:AY170327] and [GenBank:FJ237531], respectively) identified larger contig fragments for both genes [GenBank:AGKD01037727, GenBank: AGKD01030045] within the preliminary Atlantic salmon genome assembly [GenBank:GCA_000233375.1]. RepeatMasker  was used to identify potentially informative microsatellite markers within both sequences. A large dinucleotide repeat (CA)36 was located 0.5 kb downstream of elovl5a, while the closest marker identified for elovl5b was a trinucleotide repeat, (TTA)10, situated 13 kb downstream. PrimerSelect software (DNASTAR Inc., USA) was used to design flanking PCR primers for allele detection (Table 2). Fluorescently-labelled tailed primers  were used (detailed in Table 2), the allelic products being detected and sized by capillary electrophoresis (CEQ-8800 Genetic Analysis System, Beckman Coulter Inc., USA). The loci were amplified separately, with 0.5 μL of both reactions being pooled later and electrophoresed in duplex. Each reaction (6 μL final volume) comprised 50 ng template DNA, 20 nM specific forward primer (tailed), 300 nM each of specific reverse primer and fluorescently labelled tail primer and 3 μL 2× MyTaq HS mix (a Taq polymerase mastermix from Bioline Reagents Ltd, UK). Hot-start PCR was carried out on a T-gradient thermocycler, (Biometra GmbH, Germany) using the following profiles 95°C for 1 min (for Taq activation and initial denaturation), followed by 34 cycles consisting of 95°C for 15 s, annealing at 58°C (elovl5b), or 56°C (elovl5a) for 30 s, and extension at 72°C for 45 s.
The parents of six available Atlantic salmon pedigree panels were screened for both loci, with two of these (Br5 and Br5/2) being identified as informative for linkage analysis (i.e. at least one parent being heterozygous for both loci). A total of 46 progeny from family Br5 were screened (more than 400 loci having been previously mapped for this panel (SalMap project, ), while 31 progeny available for family Br5/2 were genotyped. Joint segregation analysis , as implemented in  was carried out on both panels, while linkage mapping was performed on the SalMap reference family (Br5) using JoinMap .
We thank S Carboni for the sample collection, S Morais who kindly provided salmonid expression data for analysis and Bjørn Høyheim for assisting in linkage assignments. This research was supported by a Horizon Scholarship awarded by the University of Stirling (Scotland).
- Allendorf FW, Thorgaard GH: Tetraploidy and the evolution of salmonid fishes. Evolutionary Genetics of Fishes. Edited by: Turner BJ. 1984, New York: Plenum PressGoogle Scholar
- Koop BF, Davidson WS: Genomics and the genome duplication in salmonids. Fisheries for Global Welfare and Environment 5th World Fisheries Congress. Edited by: Tsukamoto K, Kawamura T, Takeuchi T, Beard JTD, Kaiser MJ. 2008, Tokyo: Terrapub, 77-86.Google Scholar
- Ohno S: Evolution by gene duplication. 1970, Berlin: Springer-VerlagView ArticleGoogle Scholar
- Roth C, Rastogi S, Arvestad L, Dittmar K, Light S, Ekman D, Liberles DA: Evolution after gene duplication: models, mechanisms, sequences, systems, and organisms. J Exp Zool. 2007, 308 (1): 58-73.View ArticleGoogle Scholar
- Sémon M, Wolfe KH: Consequences of genome duplication. Curr Opin Genetics Dev. 2007, 17 (6): 505-512. 10.1016/j.gde.2007.09.007.View ArticleGoogle Scholar
- López JA, Chen W, Ortí G: Esociform phylogeny. Copeia. 2004, 3: 449-464.View ArticleGoogle Scholar
- Li C, Lu G, Orti G: Optimal data partitioning and a test case for ray-finned fishes (Actinopterygii) based on ten nuclear loci. Syst Biol. 2008, 57 (4): 519-539. 10.1080/10635150802206883.PubMedView ArticleGoogle Scholar
- Animal genome size database. http://www.genomesize.com,
- Ishiguro NB, Miya M, Nishida M: Basal euteleostean relationships: a mitogenomic perspective on the phylogenetic reality of the “Protacanthopterygii”. Mol Phylogenet Evol. 2003, 27: 476-488. 10.1016/S1055-7903(02)00418-9.PubMedView ArticleGoogle Scholar
- Leong JS, Jantzen SG, von Schalburg KR, Cooper GA, Messmer AM, Liao NY, Munro S, Moore R, Holt RA, Jones SJ: Salmo salar and Esox lucius full-length cDNA sequences reveal changes in evolutionary pressures on a post-tetraploidization genome. BMC Genomics. 2010, 11: 279-10.1186/1471-2164-11-279.PubMed CentralPubMedView ArticleGoogle Scholar
- Tocher DR, Bell JG, MacGlaughlin P, McGhee F, Dick JR: Hepatocyte fatty acid desaturation and polyunsaturated fatty acid composition of liver in salmonids: effects of dietary vegetable oil. Comp Biochem Physiol B Comp Biochem. 2001, 130 (2): 257-270. 10.1016/S1096-4959(01)00429-8.View ArticleGoogle Scholar
- Holman RT: Control of polyunsaturated acids in tissue lipids. J Am Coll Nutr. 1986, 5 (2): 183-211.PubMedView ArticleGoogle Scholar
- Sargent JR, Tocher DR, Bell JG: The lipids. Fish nutrition. Edited by: Halver JE. 2002, San Diego, California: Academic Press, 181-257. 3Google Scholar
- Wainwright PE: Dietary essential fatty acids and brain function: a developmental perspective on mechanisms. P Nutr Soc. 2002, 61 (1): 61-69. 10.1079/PNS2001130.View ArticleGoogle Scholar
- Dangour AD, Uauy R: N-3 long-chain polyunsaturated fatty acids for optimal function during brain development and ageing. Asia Pac J Clin Nutr. 2008, 1: 185-188.Google Scholar
- Calder PC, Albers R, Antoine JM, Blum S, Bourdet-Sicard R, Ferns GA, Folkerts G, Friedmann PS, Frost GS, Guarner F: Inflammatory disease processes and interactions with nutrition. Br J Nutr. 2009, 1: 1-45.Google Scholar
- Tocher DR: Metabolism and functions of lipids and fatty acids in teleost fish. Rev Fish Sci. 2003, 11 (2): 107-184. 10.1080/713610925.View ArticleGoogle Scholar
- Sargent JR, Bell JG, Bell MV, Henderson RJ, Tocher DR: Requirement criteria for essential fatty acids. J Appl Ichthyol. 1995, 11 (3–4): 183-198.View ArticleGoogle Scholar
- Bell JG, McEvoy J, Tocher DR, McGhee F, Campbell PJ, Sargent JR: Replacement of fish oil with rapeseed oil in diets of Atlantic salmon (Salmo salar) affects tissue lipid compositions and hepatocyte fatty acid metabolism. J Nutr. 2001, 131: 1535-1543.PubMedGoogle Scholar
- Tocher DR, Agaba M, Hastings N, Bell JG, Dick JR, Teale AJ: Nutritional regulation of hepatocyte fatty acid desaturation and polyunsaturated fatty acid composition in zebrafish (Danio rerio) and tilapia (Oreochromis niloticus). Fish Physiol Biochem. 2002, 24 (4): 309-320.View ArticleGoogle Scholar
- Kluytmans JHFM, Zandee DI: Lipid metabolism in the northern pike (Esox lucius L.) in vivo incorporation of 1-14C-acetate in the lipids. Comp Biochem Physiol B Comp Biochem. 1974, 48 (4): 641-649. 10.1016/0305-0491(74)90144-8.View ArticleGoogle Scholar
- Schwalme K: A quantitative comparison between diet and body fatty acid composition in wild northern pike (Esox lucius L.). Fish Physiol Biochem. 1992, 10 (2): 91-98. 10.1007/BF00004520.PubMedView ArticleGoogle Scholar
- Henderson RJ, Park MT, Sargent JR: The desaturation and elongation of 14C-labelled polyunsaturated fatty acids by pike (Esox lucius L.) in vivo. Fish Physiol Biochem. 1995, 14 (3): 223-235. 10.1007/BF00004313.PubMedView ArticleGoogle Scholar
- Buzzi M, Henderson RJ, Sargent JR: The biosynthesis of docosahexaenoic acid [22:6(n-3)] from linolenic acid in primary hepatocytes isolated from wild northern pike. J Fish Biol. 1997, 51 (6): 1197-1208. 10.1111/j.1095-8649.1997.tb01136.x.View ArticleGoogle Scholar
- Morais S, Monroig O, Zheng XZ, Leaver MJ, Tocher DR: Highly unsaturated fatty acid synthesis in Atlantic salmon: Characterization of ELOVL5-and ELOVL2-like elongases. Mar Biotechnol. 2009, 11: 627-639. 10.1007/s10126-009-9179-0.PubMedView ArticleGoogle Scholar
- Monroig Ó, Zheng X, Morais S, Leaver MJ, Taggart JB, Tocher DR: Multiple genes for functional 6 fatty acyl desaturases (Fad) in Atlantic salmon (Salmo salar L.): gene and cDNA characterization, functional expression, tissue distribution and nutritional regulation. Biochim Biophys Acta. 2010, 1801: 1072-1081. 10.1016/j.bbalip.2010.04.007.PubMedView ArticleGoogle Scholar
- Leaver MJ, Bautista JM, Björnsson BT, Jönsson E, Krey G, Tocher DR, Torstensen BE: Towards fish lipid nutrigenomics: current state and prospects for fin-fish aquaculture. Rev Fish Sci. 2008, 16 (1): 73-94. 10.1080/10641260802325278.View ArticleGoogle Scholar
- Carmona-Antoñanzas G, Monroig Ó, Dick JR, Davie A, Tocher DR: Biosynthesis of very long-chain fatty acids (C > 24) in Atlantic salmon: Cloning, functional characterisation, and tissue distribution of an Elovl4 elongase. Comp Biochem Physiol B Comp Biochem. 2011, 159: 122-129. 10.1016/j.cbpb.2011.02.007.View ArticleGoogle Scholar
- Engelman DM, Steitz TA, Goldman A: Identifying nonpolar transbilayer helices in amino acid sequences of membrane proteins. Annu Rev Biophys Biophys Chem. 1986, 15: 321-353. 10.1146/annurev.bb.15.060186.001541.PubMedView ArticleGoogle Scholar
- Tvrdik P, Westerberg R, Silve S, Asadi A, Jakobsson A, Cannon B, Loison G, Jacobsson A: Role of a new mammalian gene family in the biosynthesis of very long chain fatty acids and sphingolipids. J Cell Biol. 2000, 149: 707-717. 10.1083/jcb.149.3.707.PubMed CentralPubMedView ArticleGoogle Scholar
- Leonard AE, Bobik EG, Dorado J, Kroeger PE, Chuang LT, Thurmond JM, Parker-Barnes JM, Das T, Huang YS, Mukerji P: Cloning of a human cDNA encoding a novel enzyme involved in the elongation of long-chain polyunsaturated fatty acids. Biochem J. 2000, 350: 765-770. 10.1042/0264-6021:3500765.PubMed CentralPubMedView ArticleGoogle Scholar
- Agaba KM, Tocher DR, Zheng X, Dickson CA, Dick JR, Teale AJ: Cloning and functional characterisation of polyunsaturated fatty acid elongases of marine and freshwater teleost fish. Comp Biochem Physiol. 2005, 142: 342-352. 10.1016/j.cbpb.2005.08.005.View ArticleGoogle Scholar
- Monroig Ó, Li Y, Tocher DR: Delta-8 desaturation activity varies among fatty acyl desaturases of teleost fish: High activity in delta-6 desaturases of marine species. Comp Biochem Physiol B Comp Biochem. 2011, 159 (4): 206-213. 10.1016/j.cbpb.2011.04.007.View ArticleGoogle Scholar
- Zhang J: Evolution by gene duplication: an update. Trends Ecol Evol. 2003, 18 (6): 292-298. 10.1016/S0169-5347(03)00033-8.View ArticleGoogle Scholar
- Ghioni C, Tocher DR, Bell MV, Dick JR, Sargent JR: Low C18 to C20 fatty acid elongase activity and limited conversion of stearidonic acid, 18:4(n-3), to eicosapentaenoic acid, 20:5(n-3), in a cell line from the turbot, Scophthalmus maximus. Biochim Biophys Acta. 1999, 1437: 170-181. 10.1016/S1388-1981(99)00010-4.PubMedView ArticleGoogle Scholar
- Ramsden SD, Brinkmann H, Hawryshyn CW, Taylor JS: Mitogenomics and the sister of Salmonidae. Trends Ecol Evolut. 2003, 18 (12): 607-610. 10.1016/j.tree.2003.09.020.View ArticleGoogle Scholar
- Phillips BR, Keatley AK, Morasch RM, Ventura BA, Lubieniecki PK, Koop FB, Danzmann GR, Davidson SW: Assignment of Atlantic salmon (Salmo salar) linkage groups to specific chromosomes: Conservation of large syntenic blocks corresponding to whole chromosome arms in rainbow trout (Oncorhynchus mykiss). BMC Genet. 2009, 10: 46-PubMed CentralPubMedView ArticleGoogle Scholar
- ASalBase. http://www.asalbase.org,
- Lien S, Gidskehaug L, Moen T, Hayes BJ, Berg PR, Davidson WS, Omholt SW, Kent MP: A dense SNP-based linkage map for Atlantic salmon (Salmo salar) reveals extended chromosome homeologies and striking differences in sex-specific recombination patterns. BMC Genomics. 2011, 12 (1): 615-10.1186/1471-2164-12-615.PubMed CentralPubMedView ArticleGoogle Scholar
- Newcomer ME, Liljas A, Sundelin J, Rask L, Peterson PA: Crystallization of and preliminary X-ray data for the plasma retinol binding protein. J Biol Chem. 1984, 259 (8): 5230-5232.PubMedGoogle Scholar
- Fernandes JMO, Ruangsri J, Kiron V: Atlantic cod piscidin and its diversification through positive selection. PLoS One. 2010, 5: e9501-10.1371/journal.pone.0009501.PubMed CentralPubMedView ArticleGoogle Scholar
- Inagaki K, Aki T, Fukuda Y, Kawamoto S, Shigeta S, Ono K, Suzuki O: Identification and expression of a rat fatty acid elongase involved in the biosynthesis of C18 fatty acids. Biosci Biotechnol Biochem. 2002, 66 (3): 613-621. 10.1271/bbb.66.613.PubMedView ArticleGoogle Scholar
- Gregory MK, See VHL, Gibson RA, Schuller KA: Cloning and functional characterisation of a fatty acyl elongase from southern bluefin tuna (Thunnus maccoyii). Comp Biochem Physiol B Comp Biochem. 2010, 155 (2): 178-185. 10.1016/j.cbpb.2009.11.002.View ArticleGoogle Scholar
- Tu WC, Muhlhausler BS, James MJ, Stone DAJ, Gibson RA: An alternative n-3 fatty acid elongation pathway utilising 18:3n-3 in barramundi (Lates calcarifer). Biochem Biophys Res Commun. 2012, 423 (1): 176-182. 10.1016/j.bbrc.2012.05.110.PubMedView ArticleGoogle Scholar
- Tocher DR, Zheng X, Schlechtriem C, Hastings N, Dick JR, Teale AJ: Highly unsaturated fatty acid synthesis in marine fish: cloning, functional characterization, and nutritional regulation of fatty acyl delta 6 desaturase of Atlantic cod (Gadus morhua L.). Lipids. 2006, 41 (11): 1003-1016. 10.1007/s11745-006-5051-4.PubMedView ArticleGoogle Scholar
- Lauritzen L, Hansen HS, Jorgensen MH, Michaelsen KF: The essentiality of long chain n-3 fatty acids in relation to development and function of the brain and retina. Prog Lipid Res. 2001, 40: 1-94. 10.1016/S0163-7827(00)00017-5.PubMedView ArticleGoogle Scholar
- Hamilton J, Brunaldi K: A model for fatty acid transport into the brain. J Mol Neurosci. 2007, 33 (1): 12-17. 10.1007/s12031-007-0050-3.PubMedView ArticleGoogle Scholar
- Moore SA: Polyunsaturated fatty acid synthesis and release by brain-derived cells in vitro. J Mol Neurosci. 2001, 16 (2–3): 195-200.PubMedView ArticleGoogle Scholar
- Schwenk RW, Holloway GP, Luiken JJFP, Bonen A, Glatz JFC: Fatty acid transport across the cell membrane: regulation by fatty acid transporters. Prostaglandins Leukot Essent Fatty Acids. 2010, 82 (4–6): 149-154.PubMedView ArticleGoogle Scholar
- Mitchell RW, On NH, Del Bigio MR, Miller DW, Hatch GM: Fatty acid transport protein expression in human brain and potential role in fatty acid transport across human brain microvessel endothelial cells. J Neurochem. 2011, 117 (4): 735-746.PubMedGoogle Scholar
- Zheng X, Tocher DR, Dickson CA, Bell JG, Teale AJ: Effects of diets containing vegetable oil on expression of genes involved in highly unsaturated fatty acid biosynthesis in liver of Atlantic salmon (Salmo salar). Aquaculture. 2004, 236: 467-483. 10.1016/j.aquaculture.2004.02.003.View ArticleGoogle Scholar
- Mourente G, Bell JG: Partial replacement of dietary fish oil with blends of vegetable oils (rapeseed, linseed and palm oils) in diets for European sea bass (Dicentrarchus labrax L.) over a long term growth study: effects on muscle and liver fatty acid composition and effectiveness of a fish oil finishing diet. Comp Biochem Physiol B Comp Biochem. 2006, 145 (3–4): 389-399.View ArticleGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4 (4): 406-425.PubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.PubMedView ArticleGoogle Scholar
- Miyata T, Yasunaga T: Molecular evolution of mRNA: a method for estimating evolutionary rates of synonymous and amino acid substitutions from homologous nucleotide sequences and its application. J Mol Evol. 1980, 16 (1): 23-36. 10.1007/BF01732067.PubMedView ArticleGoogle Scholar
- Yang Z, Nielsen R, Goldman N, Pedersen AM: Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics. 2000, 155 (1): 431-449.PubMed CentralPubMedGoogle Scholar
- Korber B: HIV signature and sequence variation analysis. Computational and evolutionary analysis of HIV molecular sequences. Edited by: Rodrigo AG, Learn GH. 2000, Dordrecht, Netherlands: Kluwer Academic Publishers, 55-72.Google Scholar
- Nei M, Gojobori T: Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol. 1986, 3 (5): 418-426.PubMedGoogle Scholar
- Ota T, Nei M: Variance and covariances of the numbers of synonymous and nonsynonymous substitutions per site. Mol Biol Evol. 1994, 11 (4): 613-619.PubMedGoogle Scholar
- Kosakovsky Pond SL, Frost SDW: A genetic algorithm approach to detecting lineage-specific variation in selection pressure. Mol Biol Evol. 2004, 22 (3): 478-485. 10.1093/molbev/msi031.View ArticleGoogle Scholar
- Pevsner J: Bioinformatics and Functional Genomics. 2009, Hoboken: John Wiley & Sons, 2View ArticleGoogle Scholar
- Primer3. http://biotools.umassmed.edu/bioapps/primer3_www.cgi,
- Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29 (9): e45-10.1093/nar/29.9.e45.PubMed CentralPubMedView ArticleGoogle Scholar
- RepeatMasker. http://www.repeatmasker.org/cgi-bin/WEBRepeatMasker,
- Boutin-Ganache I, Raposo M, Raymond M, Deschepper CF: M13-tailed primers improved the readability and usability of microsatellites analyses performed with two different allele-sizing methods. Biotechniques. 2001, 31: 24-28.PubMedGoogle Scholar
- Danzmann RG, Cairney M, Davidson WS, Ferguson MM, Gharbi K, Guyomard R, Holm L-E, Leder E, Okamoto N, Ozaki A: A comparative analysis of the rainbow trout genome with 2 other species of fish (Arctic charr and Atlantic salmon) within the tetraploid derivative Salmonidae family (subfamily: Salmoninae). Genome. 2005, 48: 1037-1051. 10.1139/g05-067.PubMedView ArticleGoogle Scholar
- Mather K: The measurement of linkage in heredity. 1951, London: MethuenGoogle Scholar
- Taggart JB, Ferguson A: Allozyme variation in the brown trout (Salmo trutta L.): Single locus and joint segregation inheritance studies. Heredity. 1984, 53: 339-359. 10.1038/hdy.1984.92.View ArticleGoogle Scholar
- Stam P: Construction of integrated genetic linkage maps by means of a new computer package: JoinMap. Plant J. 1993, 3 (5): 739-744. 10.1111/j.1365-313X.1993.00739.x.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.