Isolation, characterization and comparison of Atlantic and Chinook salmon growth hormone 1 and 2

Background Growth hormone (GH) is an important regulator of skeletal growth, as well as other adapted processes in salmonids. The GH gene (gh) in salmonids is represented by duplicated, non-allelic isoforms designated as gh1 and gh2. We have isolated and characterized gh-containing bacterial artificial chromosomes (BACs) of both Atlantic and Chinook salmon (Salmo salar and Oncorhynchus tshawytscha) in order to further elucidate our understanding of the conservation and regulation of these loci. Results BACs containing gh1 and gh2 from both Atlantic and Chinook salmon were assembled, annotated, and compared to each other in their coding, intronic, regulatory, and flanking regions. These BACs also contain the genes for skeletal muscle sodium channel oriented in the same direction. The sequences of the genes for interferon alpha-1, myosin alkali light chain and microtubule associated protein Tau were also identified, and found in opposite orientations relative to gh1 and gh2. Viability of each of these genes was examined by PCR. We show that transposon insertions have occurred differently in the promoters of gh, within and between each species. Other differences within the promoters and intronic and 3'-flanking regions of the four gh genes provide evidence that they have distinct regulatory modes and possibly act to function differently and/or during different times of salmonid development. Conclusion A core proximal promoter for transcription of both gh1 and gh2 is conserved between the two species of salmon. Nevertheless, transposon integration and regulatory element differences do exist between the promoters of gh1 and gh2. Additionally, organization of transposon families into the BACs containing gh1 and for the BACs containing gh2, are very similar within orthologous regions, but much less clear conservation is apparent in comparisons between the gh1- and gh2-containing paralogous BACs for the two fish species. This is consistent with the hypothesis that a burst of transposition activity occurred during the speciation events which led to Atlantic and Pacific salmon. The Chinook and other Oncorhynchus GH1s are strikingly different in comparison to the other GHs and this change is not apparent in the surrounding non-coding sequences.


Background
Salmonids are used as models for studies in environmental toxicology, physiology, comparative immunology, growth, gametogenesis, olfaction and osmoregulation [1]. Although considerable knowledge of the basic biology of salmonids exists [1][2][3], the characterization of salmonid genomes will better enable scientific decisions on the conservation and enhancement of wild stocks, improve knowledge of fish health and increase the commercial viability of aquaculture. Sequencing of salmonid genomes also permits investigators to study fundamental questions concerning genome evolution, genome duplication and rearrangements, repeat-rich structures, transposon activity, gene silencing, and the re-establishment of a more stable diploid genome from a pseudotetraploid state.
Growth hormone (GH) plays a very important role in many regulatory, metabolic and developmental processes in various vertebrate tissues [4]. In salmonids, GH is the principle stimulator of skeletal growth and plays a key role in lipid mobilization, protein synthesis and feeding behaviour [5]. In various fishes, GH also manifests some functions not found in all vertebrates, such as activities that influence sexual maturation and saltwater adaptation [6,7]. The GH gene (gh) in salmonids is represented by duplicated, non-allelic isoforms designated as gh1 and gh2 [8], which diverged at least 30 million years (MY) ago [9].
In vertebrates, the major source of GH is the somatotrophic cells in the pituitary, from which it is secreted into and borne by the plasma to act on receptors throughout the organism [5,10]. Splice variants have also been isolated from extrapituitary tissue in humans (placenta, testis, blood mononuclear cells) and chicken (eye, heart) [ [11] and refs. therein]. In trout, gh transcripts have been detected in the pituitary, liver, head kidney, spleen, thymus, intestine and leukocytes and gut [12,13]. Improved understanding of the chromosomal environment in which gh genes reside in salmonids will assist our understanding of the transcriptional controls regulating these genes. By analyzing two paralogous genes (e.g. gh1 and gh2) and doing so in two species (which have been separated for approximately 20 MY; [9]), identification of common conserved regions with presumed broad functional importance can be achieved. Further, regions which are found to be conserved only between species in one paralogue type may be important for differential regulatory control between the paralogues, while differences between species within a single paralogue type may help identify regions not important to regulation. Regions which are conserved between paralogues in only one species are candidates for regions which have undergone gene conversion subsequent to divergence between the two species. While gene conversion does not appear to have occurred at salmonid GH loci when examined at the gene level [9], an analysis of gene conversion has not been examined at larger scales now possible with genomic analyses. Understanding these sequence relationships has important practical ramifications since complete genomic resources are being developed for Atlantic salmon. It is important to determine the degree of conservation of this information with other important salmonid species such as Chinook salmon for extrapolation of emerging genomic information from one species to another.
Bacterial artificial chromosomes (BACs) containing gh1 and gh2 from both Atlantic and Chinook salmon were assembled, annotated, and compared to each other in their coding, intronic, regulatory, and flanking regions. A core proximal promoter for transcription of both gh1 and gh2 is conserved between the two species of salmon. A 1600 bp insertion of a Tc1-like DNA transposon sequence was found within the promoter region of both gh2 genes, but not in the promoter of the gh1 genes. Furthermore, a Polinton-1 transposon is inserted in only the promoter for Chinook salmon GH1. Other differences within the promoters and intronic and 3'-flanking regions of the four genes provide evidence that supports the notion that they are regulated differently and thus may possess different functions. Intriguingly, Chinook salmon GH1 has undergone more than twice as many changes than any of the other GHs; changes not reflected in the surrounding noncoding DNA.

Results
Confirmation that the Atlantic salmon BACs we isolated were gh-containing BACs was performed by comparing HindIII-digested BAC fragment profiles to profiles on the internet Contig Explorer version 3.4 (iCE 3.4) database [14]. The gene for the skeletal muscle sodium channel (scn) was also identified by PCR for each isolated BAC, suggesting that each BAC contained the 5'-region upstream of each gh gene. gh type (gh1 vs. gh2) was determined for BACs in both Atlantic and Chinook salmon by paralogue-specific PCR.

BAC and GH comparisons
Atlantic salmon (AS) and Chinook salmon (CS)gh loci were analyzed using DIGIT [15], which identified the presence, location, and direction of putative genes on each BAC. Each of these putative genes were assessed by BLASTX [16] to protein databases. BACs containing the gh genes also bore the genes for scn oriented in the same direction ( Figure 1A). The sequences of the genes for interferon alpha-1 (ifna1), myosin alkali light chain (mlc) and microtubule associated protein Tau (mapt) were also identified, and found in opposite orientations relative to gh1 and gh2 ( Figure 1A). The Chinook salmon gh1 and gh2containing BACs are similarly organized.
An overall comparison of Atlantic and Chinook salmon gh1 and gh2 loci are shown in VISTA plots [17] ( Figure  1B). This figure shows many regions of high similarity separated by gaps that represent repeat sequences, deletions and insertions. Visually orthologous comparisons (AS gh1 and CS gh1, and AS gh2 and CS gh2) appear more similar than paralogous comparisons (gh1 and gh2).
Using RepeatMasker [18] and a new salmonid repeat database [19], 30% of the CS gh1, 28% of the CS gh2, 33% of AS gh1 and 35% of AS gh2 are repeated elements or simple repeats (approximately 50% of the total repeats are 1.6 kb Tc1-like transposons). With repeats removed from the genomic sequences, paralogous comparisons showed that Atlantic salmon gh1 and gh2 have 89.5% identity over 46,425 bp, and Chinook salmon gh1 and gh2 showed 87.9% identity over 52,782 bp. These alignments also revealed numerous insertions and deletions among these loci.
In comparisons of Atlantic and Chinook salmon BACs, gh1 showed 92.3% identity over 85,987 bp and gh2 showed 93.4% identity over 78,106 bp. Estimates of 93% identity between Atlantic salmon and Chinook salmon are similar to the 94% identity found between Atlantic A) The genomic organization of growth hormone 1 and 2 gene loci in Atlantic and Chinook salmon  salmon and rainbow trout (Oncorhynchus mykiss) over 125 kb [20].
Comparisons of the 210 amino acid residues of GH show very strong similarity ( Figure 2). Atlantic salmon GH1 and GH2 shows 97.1% identity, while Chinook salmon GH1 and GH2 shows 93.1% identity, and Atlantic salmon and Chinook salmon GH1 and GH2 shows 94.1% identity and 97.1% identity, respectively. An analysis of these relationships indicates that the CS GH1 has undergone more than twice as many changes than for any of the other GHs. This change in CS GH1 is not apparent in the surrounding non-coding sequences and indicates there may have been a shift of natural selection pressures on the coding region of the CS gh1.

Gene annotation
The gh paralogs were fully annotated and the genes consist of 6 exons and 5 introns. All upstream positions we document are presented in relationship to the translation start site (in negative numbers) and, due to variations in insertions and deletions between the gh genes within their intronic regions, all downstream information is based on nucleotide (nt) numbering as presented in Additional file 1. Intron 3 of the gh2 genes are much shorter than the corresponding intron for the gh1 paralogs, partly accounted for by insertion of SSsp2201 microsatellites within them. Intron 5 of CS gh1 is at least 360 nts longer than for any of the other genes. Also, numerous polymorphic loci exist

GH promoter and intron analysis
Amino acid residue alignments of Atlantic and Chinook salmon GH1 and 2 Figure 2 Amino acid residue alignments of Atlantic and Chinook salmon GH1 and 2. Identical residues designated by exclamation marks; asterisks denote strong residue similarity.
in intron 4 of the gh genes that have been previously characterized for salmonids, including AS and CS (coordinates: 2895/2917; 3060/3083; 3635/3673 and 3822/ 3865) [8]. The sequences from position 3060 to 3083 match identically [see Additional file 1], but we show that there is one nucleotide difference in the AS gh1 paralogue at position 3069. All of the other polymorphic loci match identically.
There are three polyadenylation (poly(A)) tail signal sequences (AATAAA) beginning at positions 5264, 5617 and 5717 [see Additional file 1]. All four gh genes have at least one poly(A) signal sequence, but only the gh1 genes potentially possess three canonical signal sequences.  (Figure 3). This fragment is truncated at the 3'-end, consistent with partial retrotranscription. Furthermore, an unusual insertion 10 kb upstream of only the CS gh1 gene was found ( Figure 3). This insertion is approximately 20 kb in length containing a 10 kb palindromic sequence which is part of a Polinton-1 transposon ( Figure 4) [21].

Tissue expression PCRs
To examine the functionality of the genes identified on each gh-containing BAC, primers specific for each gene were designed and used in PCRs on a panel of tissues for each species [see Additional files 2 and 3]. This was not done for genes that were clearly pseudogenic or for genes for which information was incomplete, such as the genes that encode dynein heavy chain (data not shown) and microtubule associated protein Tau. Microtubule associated protein Tau is pseudogenic due to frameshifts and stop codons within the gene in each species (data not shown).
We conducted tissue expression studies of the 7 genes we had characterized in the BACs to determine if similarities existed in expression pattern between the two species; if expression was confined to specific tissues and whether indeed each of these genes were functional. Of the 7 genes examined, the expression of the gh paralogues was the least diverse. The expression of gh1 and gh2 was restricted to only muscle, brain, gill or eye (pituitary gland was not examined). Expression of the scn genes was not restricted to just skeletal muscle, but indeed was found in every tissue except the skin. It is also evident that ifna1 has fairly broad expression in the tissues examined here.
Myosin alkali light chain 1 also was expressed in a large range of different tissues. However, myosin alkali light chain 2 (mlc2) appears to not be expressed at detectable levels in any of the tissues we examined in either species, despite in silico analysis indicating the viability of the mlc2 transcript. Each gene that encodes mlc2 appears intact, with no obvious problems for transcription or translation.

GH transcription
Transactivation and footprinting studies have been conducted on the gh promoters of both Chinook salmon and rainbow trout [22,23]. Assembly of at least part of the transcriptional machinery was delineated to binding elements between positions -300 to -30. It was demonstrated that cyclic-AMP/PKA regulation of these promoters required nucleation and recruitment of both pituitary transcriptional activator-1 (Pit-1) and cAMP-response element (CRE) binding proteins (CREB) [22,23].
We note strong sequence identity of potential Pit-1, CREB, and retinoic acid receptor (RAR)/retinoid × receptor (RXR) binding elements that form the core of the proximal gh promoters of both Atlantic salmon (AS) and Chinook salmon (CS). An asymmetric CRE (TGACG) at -217/ -213 is flanked by at least four Pit-1 binding elements [see Additional file 1]. The more distal Pit-1 binding elements follow more strictly the core consensus sequence for Pit-1 (W 3 CAT) that is conserved in fish and mammals [24]. The potential RAR/RXR response element present at position -190/-171 could hypothetically permit interaction of chicken ovulbumin upstream promoter transcription factor (COUP-TF), or retinoic acid, estrogen or thyroid hormone receptors with an assembled Pit-1/CREB complex. Different variations of GTCT-rich blocks immediately upstream of each proximal promoter, as well as within the introns of each gene, may make regulatory or structural contributions to gh transcription. For example, some of the CTCT blocks found in these regions could serve as binding elements for GAGA-binding factors [25]. However, differences do exist between the AS and CS gh1 and gh2 promoters. For example, there is the potential for Pit-1 to bind a reverse-oriented response element roughly 250 bp further upstream in each gh2 promoter that is not present in the gh1 promoters. No palindromic CREs are found within 15 kb upstream of the gh genes, but one is present in intron 4 of the gh1 genes [see Additional file 1]. The different CRE, ERE and Pit-1 motif locations suggest that the gh1 and gh2 paralogues may have distinct regulatory modes and possibly act to function differently and/or at different times of salmonid development.

GH expression and function
The more traditional models that have regarded divergence among duplicate genes as a means for acquisition of new functions have met with recent challenges. For example, it has been demonstrated that mutation in many duplicated genes of rice is much slower than in 'singletons' and less drastic in terms of amino acid change [26]. Conservation of the apparent crucial functionality of the duplicates outweighs processes that lead to neofunctionalizaton. Indeed, it may be more common for gene nov- At the amino acid level, CS GH1 appears to be undergoing change much faster than any of the other GHs ( Figure 2). In terms of retention of amino acid residue identity, the CS GH1 ORF has more changes than any of the other translated regions (nine differences overall compared to no more than three for any other gene). These changes do alter the polarity and charge of the proteins and therefore may subtly effect the kinetics and affinity of hormone binding to their receptor(s). We did find colocalized expression of both cDNAs in the CS gill and brain [see Additional files 2 and 3]. Interestingly, both GH receptors have been found in various regions of the brain and in the gill filaments of trout [10]. Furthermore, although the very limited amount of change between all four genes indicates that their functions are conserved, we do observe that GH1 and GH2 may potentially have different subfunctions. Our RT-PCR work demonstrates that some tissue expression segregation has occurred, whereby only gh2 appears to be expressed in muscle or eye [see Additional files 2 and 3].

Transposon integration in Atlantic and Chinook salmon gh promoter regions
To determine the function of locally expressed gh2 in the gill, it will be necessary to elucidate the specific cells GH synthesis is associated with and to determine whether local production is a component of known osmoregulatory processes or some other activity (see [7]). It is also possible that the gh2 expression observed in the gill is from leukocyte accumulation as leukocytes have been The distribution of DNA transposons and LINE sequence elements in gh1-and gh2-containing BACs Figure 4 The distribution of DNA transposons and LINE sequence elements in gh1-and gh2-containing BACs. The growth hormone coding sequence (consisting of 6 exons) is aligned for the four BACs and connected by a light blue line. The conservation of the order of transposon families is indicated by lines, connecting corresponding sequences.

Chinook salmon Atlantic salmon
Chinook salmon Atlantic salmon Growth Hormone 2 Growth Hormone 1 shown to produce both gh1 and gh2, with gh2 the predominant transcript [28].

Expression of genes linked with GH
Much less is known about the expression of the skeletal muscle sodium channel (scn), microtubule associated protein Tau (mapt), myosin alkali light chain (mlc) and interferon alpha-1 (ifna1) genes in salmonids. To the best of our knowledge, for the other genes examined in this study, only the promoters of two ifn genes have been characterized in salmon. Two different AS type 1 ifn genes have been demonstrated to each possess two promoters that respond to different phases of infection or are tissue-specific in their transcriptional contexts [29].
The scn gene has been characterized in zebrafish [30] and fugu [31]. Some subfunctional partitioning of the four sets of duplicated zebrafish scn genes does exist, but there is still expression overlap across neuronal and muscle tissues [32]. The expression of the scn genes in the salmonids examined here also are not restricted to skeletal muscle, but their activities are found in a large number of different tissues [see Additional files 2 and 3].
There appears to be a number of ifn genes in both the Atlantic and Chinook salmon GH1 BACs analyzed (data not shown). At least one functional ifna1 gene and two to five different ifna genes (that might be pseudo or partial genes) exist at these loci for each species. It should be noted, however, that this region was difficult to sequence and assemble due to repetitive regions, therefore we cannot conclude unequivocally on the number of ifn genes present at these loci. A genomic BAC clone containing the ifna1 gene and seven other ifn genes has been characterized in Atlantic salmon (Robertsen, B: 7 th International Symposium on Fish Immunology. Abstract, 2007). However, to date, we find no other ifna genes present in our database other than ifna1 to enable us to conclude that expression of these other ifn genes does occur. Zebrafish also has two to three ifn genes in the 3'-flanking region and the scn gene in the 5'-flanking region of the gh gene on chromosome 3 [33].

Genome rearrangements
The gh2 genes have undergone some rearrangements in comparison to the gh1 paralogues (differences in promoter sequences, partly due to variable transposon integrations; deletions or lack of microsatellite insertions in intron 3; intronic microsatellite variations in intron 4). Part of this may be due to chromosome structure or positioning [8]. Also, some of the genomic variations observed between the gh1 and gh2 paralogues might be due to vestigial effects caused by a rearrangement that led to the rise of the Y-linked gh2 pseudogene [34]. Whether the additional transposon insertions found in the AS gh2 promoter can account for why we do see gh-Y in CS, but not in AS [35], will require characterization of other Pacific salmon gh regions. It should be noted, as an example, that dramatic differences have also been found between the promoters of the salmon form of gonadotropin-releasing hormone in AS and sockeye salmon [36].
An indication of the changes that can occur when sequences of various transposon families integrate into the promoters of genes are shown in Figure 3. Insertions such as these may help to shape the dynamism of genomes by influencing modes of transcription through the introduction of new regulatory elements, by changing the distances between functional binding elements or by disabling them. It can be envisaged how alterations to transcription through these integration events could lead to changes in the function of their translated products over time. For example, note there are at least two different transposon sequences found within the AS gh2 proximal regulatory sequences that are not present in any of the other gh promoters (Figure 3). These integrated sequences might impact the observed specific expression of gh2 in AS compared to CS [see Additional files 2 and 3]. Also, the presence of a Tss transposon in the 5' sequence of each gh2 and an HpaI short interspersed nuclear element (SINE) in intron 3 of each gh1 is consistent with speciation events occurring after the gene or genome duplication events (Figures 3 and 4).
Vertebrate genomes contain large numbers of transposon sequences. The mouse and human genome contain approximately 3 million such sequences and zebrafish approximately 1.4 million. However, transposons occur in certain genome regions more frequently than in socalled transposon-free regions (which range from 5 to 66 kb in length), frequently associated with developmental genes [37,38]. This association with such genes suggests that transposons are deleterious to gene integrity possibly through increased mutagenesis or recombination, or by changing regulation. The presence of the Polinton sequence found upstream of the CS GH1 gene might therefore affect its genomic stability.

Conclusion
This is one of the first in a series of studies that are needed to document coding and non-coding changes that have occurred subsequent to a whole genome duplication. As GH has been the focus of biotechnological advances in aquaculture, it is important to investigate expression changes and genomic organizational changes in important economic traits such as growth. In this paper, we report on the impact of repeat elements and transposon integrations and show that in Oncorhynchus species the GH1 duplicate has undergone a higher rate of change.
Genomes, particularly salmonid genomes, are dynamic. Part of the dynamism is the result of the purported whole genome duplication and the past integration of SINEs and LINEs into them. LINEs and SINEs introduce repetitive elements that may misalign causing unequal recombination to occur, and thereby introduce deletions and insertions among duplicated genomes. Coupled with this is the integration of a wide assortment of transposon sequences as described here and elsewhere [19]. Over time, these genomic rearrangements have led to the creation of pseudogenes, structural differences between duplicated genes and to differential regulation of paralogues.
In AS and CS, we show that transposon insertions have occurred differently in the promoters of gh, within and between each species. The organization of transposon families in the BACs containing gh1 and in the BACs containing gh2, is very similar within orthologous regions ( Figure 4). However, much less conservation is apparent in comparisons between the gh1-and gh2-containing paralogous BACs for the two fish species (Figure 4). The appearance of repeated elements and differential rates of change in the gh1 and gh2 regions is consistent with the hypothesis that a burst of transposition activity occurred during the speciation events which led to Atlantic and Pacific salmon [19]. Genome and gene duplication (30 to 100 MYA; [39,40]) has taken place much earlier than the speciation (14 to 23 MYA; [9,41,42]) and these data suggest that after gene duplication, gh1 in Oncorhynchus has evolved much faster than gh2, possibly because the genomic region for gh1 has undergone more reorganization compared to the region containing gh2.

Methods
Atlantic salmon (AS) CHORI-214 [43] and Chinook salmon (CS) CHORI-217 [44] bacterial artificial chromosome (BAC) libraries were obtained from BACPAC Resources, Children's Hospital Oakland Research Institute (CHORI) [45]. AS BAC library filters were hybridized with an oligonucleotide probe (5'-TCCCAAACAAACAG-CAACATACTCAACCGACCACCGCACT-3') designed from the AS gh1 EST [GenBank:X61938] that had been end-labeled with gamma-32 P-ATP (Amersham). The CS BAC library filters were screened using a gh cDNA probe (GH2-8) generated from sockeye salmon genomic sequence by a PCR-based intron deletion methodology [35]. The CS GH2-8 probe was labelled with alpha-32 P dCTP using the RediprimeII Random Prime Labelling system, and purified through a Probequant G-50 micro column (Amersham). Filter hybridizations were conducted as described by CHORI [45]. Probed BAC library filters were visualized using a Molecular Dynamics Storm Phos-phorImaging system.

BAC DNA confirmation
Confirmation of AS gh-containing BACs was performed by comparisons of HindIII restriction digests of the isolated clones to in silico digests for each BAC. The entire Salmo salar genome BAC library has been digested by Hin-dIII and fragment profiles are available on the internet Contig Explorer version 3.4 (iCE 3.4) database [14]. Fingerprinting of the CS BACs was carried out by SnaPshot (ABI) labeling of restriction-digested fragments and samples were analyzed on an ABI 3130xl genetic analyzer. Data was processed utilizing Genemapper and fingerprinted contigs assembled as described by Luo et al [46].
PCRs were performed using 200 ng of BAC DNA utilizing hot start PCRs: reaction mixtures without any enzyme were heated to 94°C for 3 min, cooled to 80°C and then 0.625 U of Taq DNA polymerase (Invitrogen) was added. Each PCR then underwent 35 cycles of the following parameters: 94°C for 20 sec, 50°C for 30 sec and 72°C for 45 sec. AS gh1 was contained on BACs 11-I-04 and 73-D-15 and gh2 in BAC 63-I-10. The CS gh1 was localized to BAC 108-O-24 and gh2 to BAC 206-E-17.

BAC preparation and library construction
BAC DNA was isolated by an alkaline lysis procedure using Nucleobond columns (Clontech) following the manufacturer's protocol. The isolated BAC DNA was nebulized and the DNA ends were made blunt by filling with T4 polymerase. The blunt-ended, repaired DNA was size fractioned by electrophoresis and the gel region corresponding to 1.6 to 4.0 kb was excised and gel purified (Qiagen). The fragments were blunt-end ligated into pUC19 plasmid cut with HincII (NEB) and transformed into electrocompetent DH5α E. coli cells using a Bio-Rad Gene Pulser system. Library quality was evaluated and high redundancy plating was followed by large-scale colony picking (Genetix). Extracted recombinant plasmid templates were sequenced on an ABI 3730 DNA sequencer.

BAC contig assembly
Bases were called using PHRED [47,48]. High quality sequence reads were assembled using PHRAP [49] and then viewed and edited using Consed [50]. Some gaps in BAC assembly were filled by designing primers to the contiguous sequence ends, followed by amplification of the BAC region by PCR and subsequent cloning and sequencing of the fragments. Each BAC has been deposited in GenBank as follows: Dotter [51] and PipMaker [52] were used to compare each BAC sequence to itself and to identify duplicated and repeated regions. Identification of other repeat elements was done with RepeatMasker [18] using repeat library 4.01 from Repbase [53], as well as a salmonid repeat database [19]. Gene location and direction on each BAC was determined using Digit Integrated Gene Identification Tools (DIGIT) [15].

GH alignment, annotation, and comparison
AS gh1 [GenBank:AY614010] and gh2 [GenBank: M21573] gene sequences were used to align the gh1 and gh2 paralogues with BioEdit [54]. BioEdit was also used to annotate the two paralogues and calculate percent similarities between the two genes, their coding sequences, and their amino acid sequences. Searches and comparisons of regulatory elements for the two paralogues were also performed using BioEdit. Identification of transposon insertions was performed using Dotter plots, comparing 15 kb of each promoter region with salmon transposon sequences [19].

Reverse transcription and cDNA amplification
Total RNA was extracted in TRIzol reagent (Invitrogen) from flash-frozen, adult AS (Mowi stock, DFO, West Vancouver) and CS (Chehalis River Hatchery) kidney, muscle, skin, gut, gill, spleen, brain, heart, testis, liver, eye and pyloric caecum tissues. The extracted total RNAs were cleaned using MEGAclear (Ambion) and then quantified and quality-checked by spectrophotometer and agarose gel, respectively.
The cDNAs were synthesized in 25-μL reactions from 1.0 μg total RNA using oligo(dT) 15 (Promega) and Supercript II RNase Hreverse transcriptase according to the manufacturer's instructions (Invitrogen). The reactions were incubated at 37°C for 90 min and the transcriptase heatinactivated at 70°C for 30 min. Approximately 200 ng of cDNA was used in each 25-μL PCR reaction containing 1.25 U Taq polymerase, 1 × Taq buffer, 1.25 mM MgCl 2 , 10 mM dNTPs (Invitrogen) and 15 pmol of each genespecific 5' and 3' primer [see Additional file 4]. Each PCR was carried out under the following cycling parameters: 94°C for 2 min, then 35 cycles of 94°C for 30 sec, 55°C for 30 sec, and 72°C for 1 min using a Perkin Elmer 9600. The AS scn1 and scn2 PCRs were similarly amplified but the anneal temperature used was 57.5°C. The integrity of each cDNA used was confirmed by control PCRs using ubiquitin primers.
The PCR products were separated by electrophoresis on 1.25 to 1.30% agarose gels and photographs were stored using an UVP GelDoc-It documentation system (UVP). Representative products were isolated and cloned into pCR2.1-TOPO vector (Invitrogen) and sequenced to confirm gene identities.