The genetic basis of salinity tolerance traits in Arctic charr (Salvelinus alpinus)
© Norman et al; licensee BioMed Central Ltd. 2011
Received: 3 June 2011
Accepted: 21 September 2011
Published: 21 September 2011
Skip to main content
© Norman et al; licensee BioMed Central Ltd. 2011
Received: 3 June 2011
Accepted: 21 September 2011
Published: 21 September 2011
The capacity to maintain internal ion homeostasis amidst changing conditions is particularly important for teleost fishes whose reproductive cycle is dependent upon movement from freshwater to seawater. Although the physiology of seawater osmoregulation in mitochondria-rich cells of fish gill epithelium is well understood, less is known about the underlying causes of inter- and intraspecific variation in salinity tolerance. We used a genome-scan approach in Arctic charr (Salvelinus alpinus) to map quantitative trait loci (QTL) correlated with variation in four salinity tolerance performance traits and six body size traits. Comparative genomics approaches allowed us to infer whether allelic variation at candidate gene loci (e.g., ATP1α1b, NKCC1, CFTR, and cldn10e) could have underlain observed variation.
Combined parental analyses yielded genome-wide significant QTL on linkage groups 8, 14 and 20 for salinity tolerance performance traits, and on 1, 19, 20 and 28 for body size traits. Several QTL exhibited chromosome-wide significance. Among the salinity tolerance performance QTL, trait co-localizations occurred on chromosomes 1, 4, 7, 18 and 20, while the greatest experimental variation was explained by QTL on chromosomes 20 (19.9%), 19 (14.2%), 4 (14.1%) and 12 (13.1%). Several QTL localized to linkage groups exhibiting homeologous affinities, and multiple QTL mapped to regions homologous with the positions of candidate gene loci in other teleosts. There was no gene × environment interaction among body size QTL and ambient salinity.
Variation in salinity tolerance capacity can be mapped to a subset of Arctic charr genomic regions that significantly influence performance in a seawater environment. The detection of QTL on linkage group 12 was consistent with the hypothesis that variation in salinity tolerance may be affected by allelic variation at the ATP1α1b locus. IGF2 may also affect salinity tolerance capacity as suggested by a genome-wide QTL on linkage group 19. The detection of salinity tolerance QTL in homeologous regions suggests that candidate loci duplicated from the salmonid-specific whole-genome duplication may have retained their function on both sets of homeologous chromosomes. Homologous affinities suggest that loci affecting salinity tolerance in Arctic charr may coincide with QTL for smoltification and salinity tolerance traits in rainbow trout. The effects of body size QTL appear to be independent of changes in ambient salinity.
The life history of anadromous salmonids entails migration between freshwater and seawater environments. To ensure that internal ion concentrations remain homeostatic in the face of abrupt changes in ambient salinity, an individual's osmoregulatory mechanisms must switch between states of ion absorption (i.e., hyper-osmoregulation) and ion excretion (i.e., hypo-osmoregulation). These changes are particularly important in salmonids whose reproductive cycle is dependent upon anadromous behaviour. As described by the seawater mitochondria-rich cell model, hypo-osmoregulation is primarily achieved by mechanisms associated with mitochondria-rich cells in gill tissue . In conjunction with pavement cells and accessory cells, mitochondria-rich cells form the epithelial layer of gill tissue, where the removal of Cl- is facilitated by three interdependent membrane-bound ion transporters. An electrochemical gradient across the basolateral membrane is actively maintained by Na+/K+-ATPase pumps that exchange intracellular Na+ for extracellular K+, thereby driving Na+/K+/2Cl- (NKCC) cotransporters to move Cl- from blood plasma into the cell. Once the intracellular electrochemical equilibrium is reached, Cl- exits the cell through cystic fibrosis transmembrane conductance-like regulator anion channels (CFTR) embedded in the apical membrane. Na+ secretion is thought to occur passively through leaky cation-selective paracellular pores between accessory cells and mitochondria-rich cells [1, 2]. A current hypothesis suggests that alternate claudin isoforms may confer differential permeability characteristics at these junctions [3, 4], as evident in kidney tissue .
Research at the molecular level has largely focused on Na+/K+-ATPase genes, where two isoforms in particular have important implications for hyper- and hypo-osmoregulation [6, 7]. Seawater immersion is correlated with transcriptional up-regulation of ATP1α1b, and down-regulation of ATP1α1a, whereas the reciprocal pattern is evident in freshwater. This phenomenon has been observed in rainbow trout (Oncorhynchus mykiss) [6, 7], Atlantic salmon (Salmo salar) [7–9], Arctic charr (Salvelinus alpinus) , and brown trout (Salmo trutta) .
Variation in salinity tolerance is well documented within and among the Salmoninae (Oncorhynchus, Salmo and Salvelinus) [7, 8, 11–14]. Salvelinus species are typically far less anadromous than Salmo species  and are considered to be less efficient osmoregulators [7, 16]. The evidence that variation in salinity tolerance capacity of salmonids is attributed in part to underlying genetic variation is mostly indirect. For example, different strains of Arctic charr show wide variation in seawater-induced mortality, where the salinity tolerance capacity of some individuals is commensurate with that of Atlantic salmon . Other experiments with salmonids conclude that the basis for such variation likely has a genetic component [18–20]. The most convincing evidence comes from a study on rainbow trout, where variation in genomic regions was associated with seawater-induced fluctuations in Na+ and Cl- ion concentrations in blood plasma .
The rich evolutionary history of salmonids has provided multiple opportunities for the pseudogenization, neofunctionalization and subfunctionalization of genes involved in osmoregulation. The 1-2-4 model of vertebrate evolution demonstrates that the vertebrate ancestor underwent two genome duplication events (i.e., 1R, 2R) prior to the divergence of gnathostomes . Support for a third, fish-specific genome duplication event (3R; the 1-2-4-8 model) is evident in the striking pattern of doubly conserved synteny blocks among the genomes of Tetraodon nigroviridis and Homo sapiens . It is now generally accepted that a fourth, salmonid-specific genome duplication (4R), in the form of an autotetraploidization event, occurred in the salmonid ancestor some 25-100 million years ago . Although disomic inheritance has been largely restored in modern salmonids, segregation patterns characteristic of tetrasomic inheritance and the occurrence of multivalent formations during meiosis continue to be observed [24, 25]. Multiple Na+/K+-ATPase genes have been detected and mapped within the genomes of Atlantic salmon and rainbow trout, where some duplicated loci reside on homeologous linkage groups [26, 27]. Although unmapped, transcription patterns of ATP1α1a and ATP1α1b isoforms in gill tissue provide evidence that they are conserved in Arctic charr . Furthermore, comparative analysis among homologous chromosome blocks within the salmonids can facilitate a more direct comparison of the Na+/K+-ATPase genes with model species such as rainbow trout [28, 29].
A complete physiological model inclusive of all the relevant mechanisms that affect variation in salinity tolerance capacity would be ideal prior to characterizing their underlying genetic bases. Nonetheless, the mitochondria-rich cell model remains useful, since genome scans coupled with quantitative trait locus (QTL) analyses can be utilized to identify chromosomal regions correlated with variation in salinity tolerance performance traits. QTL analyses are necessary as the traits by which seawater performance is assessed can not be objectively divided into discrete groups, for their distributions are typically continuous, and as such are not amenable to study using classical Mendelian genetics. Once QTL are identified, comparative genomics methods can determine if they occupy regions sharing homology with chromosomes holding salinity tolerance candidate genes (e.g., ATP1α1a, NKCC1, CFTR, cldn10e) in the sequenced genomes of 3R teleosts, such as zebrafish (Danio rerio) and medaka (Oryzias latipes). Confirmation of homology would advance hypotheses that variation in salinity tolerance capacity is a function of allelic variation at those loci.
In this study we employ a genome-scan approach to identify the genomic regions in Arctic char that correlate with variation in four salinity tolerance performance traits: Na+/K+-ATPase activity, blood plasma osmolality, and specific growth rates from two distinct time intervals. We describe multiple QTL for each trait among several linkage groups, and find that some QTL co-localize with the putative locations of candidate genes predicted by comparative genomics with zebrafish and medaka. We also describe putative homologous QTL from comparisons with smoltification and salinity tolerance QTL locations in the rainbow trout genome [21, 30].
Four sex-specific genetic maps were generated from the parents of two full-sib families (i.e., 10, 12; see Additional files 1, 2, 3, 4). Each map was comprised of over 100 markers arranged into a minimum of 27 linkage groups. Unlinked markers represented three to four linkage groups per map. Linkage group names were ascribed following designations from existing Arctic charr maps [28, 31]. A total of 35 linkage groups from a potential 39 in Arctic charr were represented by at least one informative marker among all families.
We found several QTL for seawater traits (i.e., Na+/K+-ATPase activity, blood plasma osmolality, specific growth rate 1, and 2) over 26 linkage groups. Among these, genome-wide significant QTL were restricted to AC-8, -14 (specific growth rate 2) and -19 (specific growth rate 1). Co-localization of QTL for multiple traits occurred on nine linkage groups (i.e., AC-1, -4, -5, -7, -18, -20, -22, -26, and 32). We also detected several QTL on linkage groups exhibiting homeologous affinities [28, 29] (see Additional file 8) suggesting the possibility that salinity tolerance may be affected by functional gene duplicates derived from the 4R genome duplication.
The QTL effect on AC-12 may also be associated with variation in claudins. It has been suggested that permeability characteristics of Na+ ions through paracellular junctions between mitochondria-rich cells and accessory cells in gill tissue may be correlated with claudin isoform identity [3, 4]. In fact, among the claudin 26 isoforms identified in Atlantic salmon, qPCR analysis shows that mRNA of cldn10e is significantly elevated after seawater exposure . Interspecific homologies suggest that the genomic positions of ATP1α1b and claudin isoform 10e (cldn10e) may have converged to a single linkage group in salmonids despite the fact that their apparent homologues are located on separate linkage groups in both zebrafish (Zv9 database; ATP1α1b on Dr-9, GenBank:NP_571765; cldn10e on Dr-6, GenBank:XM_678711) and medaka (MEDAKA1 database; ATP1α1 on Ol-2, ENSEMBL:ENSORLG00000002047; cldn10 on Ol-21, ENSEMBL:ENSORLG00000017717). Further, Dr-6 and Ol-2 show extensive synteny (supported by ≥ 3 markers) with RT-12 and AS-22 where the Na+/K+-ATPase α-isoform clusters are located [26, 27], which suggests that cldn10e and ATP1α1b may reside on the same linkage group in salmonids. Chromosomes Dr-6/9 and Ol-2/21 are also derived from the C ancestral chromosome grouping in teleost fishes, and this ancestral lineage is the most conserved syntenic block on the rainbow trout linkage groups RT-12p/16p and Atlantic salmon AS-22qb (homologous with AC-12/27) chromosome blocks . This further supports a likely conserved evolutionary origin for the ATP1α1b and cldn10e genes.
Genomic homologies among Arctic charr linkage groups containing seawater QTL with the genomes of zebrafish and Atlantic salmon suggest that the QTL on AC-4 may stem from the combined effects of multiple candidate loci. Two isoforms of NKCC1 (a and b) are known to exist , one of which has been mapped to AS-17 in Atlantic salmon (unpublished data) . Among the Arctic charr linkage group that share homology with AS-17 (i.e, AC-4, -9, -21, and -28) , AC-4 (i.e., Na+/K+-ATPase activity, specific growth rate 1 and 2), -21 (i.e., specific growth rate 2) and -28 (i.e., Na+/K+-ATPase activity) contain salinity tolerance QTL. Considering the role of NKCC1 in ion excretion, putative effects related to allelic variation at this locus would more likely be reflected by variation in Na+/K+-ATPase activity levels rather than growth measurements. Thus if variation at NKCC1 affected salinity tolerance capacity in the present study, it may not be unexpected to find QTL on AC-4 and/or -28 arising from allelic variation at this locus. Presently, a single NKCC1 locus is annotated in the zebrafish genome, whereas two NKCC1 loci have been characterized in European eel (Anguilla anguilla) . This suggests that the tetraploid salmonid ancestor might have possessed up to four functional NKCC1 genes. As such it is possible that the Arctic charr genome could contain multiple functionally active NKCC1 isoforms, which could be associated with QTL on AC-4 and -28. This is also supported by the fact that a putative homeologous affinity has been detected between AC-4/28 .
Comparative genomics suggests that a CFTR locus may also reside on AC-4. BLASTN searches of the CFTR gene in zebrafish (ENSEMBL v.59; http://www.ensembl.org/) identified significant homology with complete cDNA sequences of CFTR-I (GenBank:AF155237) and CFTR-II (GenBank:AF161070) from Atlantic salmon . Both copies coalesced to the same region on Dr-18 (CFTR-I, E-value = 5.1E-95; CFTR-II, E-value = 1.2E-106; Zv9; ENSEMBL:ENSDARG00000041107), supporting the annotation of only a single CFTR gene in zebrafish. Similarly, in medaka and stickleback single gene copies for CFTR have been localized to an unassigned scaffold (E-value = 0) and to chromosome Ga-XIX at position 10.186Mb (E-value = 0), respectively. This also suggests that duplicates of CFTR in salmonids may be derived from the 4R whole-genome duplication. The known locations of the CFTR gene in zebrafish and stickleback suggest an origin of this gene from the K ancestral lineage of teleost fishes which may thus share homology to linkage groups RT-7, -15, and 27q and possibly more extensive regions on RT-6p . Linkage group RT-27q shares some homology with AC-19a and/or AC-19b and is syntenic with AC-4a [28, 29], suggesting that these linkage groups may house CFTR duplicates. Incongruence of trait-specific QTL among AC-4 (i.e., Na+/K+-ATPase activity, blood plasma osmolality, specific growth rate 1and 2) and -19 (i.e, specific growth rate 1) suggests that the QTL on AC-4 confers a stronger effect on mitochondria-rich cell physiology than that of AC-19, indicating that AC-4 is a better candidate for the location of CFTR. However, it is also possible that each of AC-4 and -19 may contain a CFTR isoform, for as mentioned previously, two CFTR genes have been detected in the Atlantic salmon genome . Moreover, over a small region AC-4 and -19 appear to share a homeologous affinity . As postulated for ATP1α1b and cldn10e, copies of NKCC1 and CFTR may have converged to the same linkage group in salmonids (e.g., AC-4 in Arctic charr), or retained expression of one of the duplicate copies of these genes (i.e, on AC-4), if derived from AC-4/19 duplicates, despite being on separate linkage groups in zebrafish (i.e., Dr-8 and -18).
Comparison of salinity tolerance and body size trait QTL (P ≤ 0.05) in Arctic charr (Salvelinus alpinus) with growth-related genes mapped by Moghadam et al. (2007)
Salinity Tolerance QTL
Body Size QTL
BW3, K1, K2
NKAA, OSMO, SGR1, SGR2
BW1, BW2, BW3
OSMO, SGR1, SGR2
BW1, BW2, BW3, K1, K2, K3
BW1, BW2, K2
Allelic diversity at IGF2 loci may also contribute to variation in salinity tolerance capacity in non-salmonids. Transcriptional responses in multiple teleost species demonstrate that mRNA production levels of several genes not included in the mitochondria-rich cell model are contingent upon seawater exposure [40–43]. IGF2 has been implicated in salinity tolerance capacity in black-chinned tilapia (Sarotherodon melanotheron), where it was reported that gill tissue IGF2 mRNA levels were elevated in response to seawater .
Some of the salinity tolerance QTL reported here may co-localize to homologous linkage groups containing smoltification QTL in rainbow trout . This is evident for the RT-16p and AC-12/27 [28, 29] linkage groups, which share homologous affinities. The genetic marker associated with variation in skin reflectance (i.e., Omy77DU) on RT-16  maps to within 0.8 cM of the EST linked to a QTL for blood plasma osmolality on AC-27 (i.e., CA345149) . These QTL may derive from a single locus and thereby be pleiotropic, or alternatively, may be part of a salinity tolerance or smoltification-related gene-cluster. Intriguingly, thyroid hormone levels have been associated both with purine deposition (i.e., skin reflectance) and changes in Na+/K+-ATPase activity [15, 38], suggesting that, irrespective of genomic position, these loci may be regulated by thyroid hormone.
Certain salinity tolerance QTL shared tentative homologies with smoltification QTL in rainbow trout. The rainbow trout reference map reveals that a major QTL associated with multiple smoltification-related traits was linked to markers on RT-20q (i.e., OMM5108, OMM5017) [30, 32], which has homologous affinities with AC-13b (a putative metacentric) and AC-34 [28, 29]. AC-13 contains a QTL for specific growth rate 2 (i.e., OMM1174-OMM5180) that spans both the a- and b-arms while a QTL for specific growth rate 1 has been localized (i.e., OMM1657-OMM5312) on AC-34.
Other genome-wide significant QTL showed little to no homology with smoltification QTL in rainbow trout. Although the marker linked to a specific growth rate 2 QTL on AC-8b (i.e., OmyRGT6TUF) is found on RT-17 (centromeric regions) and RT-22p , where QTL for reflectance (i.e, OMM1117) and body shape (i.e., Ssa289) were localized , the respective markers are each > 30 cM away from OmyRGT6TUF , indicating that these QTL are likely not homologous. Also, a QTL on AC-14 was not homologous with any smoltification QTL in rainbow trout. These discrepancies could be related to differences in experimental design, for Nichols et al.  exposed their fish only to freshwater. This suggests that the effects conferred by the aforementioned QTL may be specifically induced by seawater exposure, and furthermore, that the genes involved in preparing a fish for migration to seawater (i.e., smoltification) are not necessarily the same as those involved in hypo-osmoregulation. Differences could also be due to species, ontological status, environmental conditions, or any combination thereof.
For the salinity tolerance QTL in rainbow trout described by Le Bras et al. , putative QTL homologies were apparent with the chromosome-wide and genome-wide QTL detected from the combined analysis in Arctic charr within the following regions: RT-4 (AC-4); RT-10 (AC-15); RT-25 (AC-4); and RT-26 (AC-4, -25) for physiological traits. Homologies with QTL for gill weight on RT-9q (AC-20), RT-9p (AC-23) and RT-26 (AC-4, -25 and possibly -22) were also evident. Also, the strong QTL region detected on AC-19 may share homology to the QTL regions on RT-19 although the Arctic charr markers appear to lie more in the central cluster of RT-19 linkage group, while those reported span the RT-19q arm . Similarly a homeologous region to AC-19 (i.e., AC-32) overlaps the QTL region reported on RT-23q . Multiple homology assignments for AC-8 exist within the rainbow trout genome (i.e., RT-2p; RT-3q; RT-7q and RT-10q) with the largest tract homologous to RT-3q. Nonetheless, homologies to the region reported on RT-10  cannot be excluded at present. Finally, AC-14 shares homology with RT-24q and also shares marker synteny with two locations around the centromere on RT-19 and the p-arm. This region may be related to the QTL localized to RT-19 .
We identified genome-wide significant body size QTL on AC-1 (body weight 1), -19 (body weight 1 and 2), -20 (condition factor 3), and -28 (condition factor 1 and 2). A marker linked to a QTL for body weight 1 on AC-1 (i.e., BX311884) was previously associated with body weight QTL in Arctic charr . Further, the association of AC-1 with QTL for Na+/K+-ATPase activity, specific growth rate 1 and 2 is suggestive of its relevance to salinity tolerance. The close proximity of QTL for body weight 1 (BX311884) and Na+/K+-ATPase activity (OMM1330) suggests they may be controlled by a pleiotropic locus, for their respective markers are separated by only 2.1 cM [28, 29]. However, at this time alternative explanations such as tight linkage among loci cannot be ruled out, since each marker interval will likely contain clusters of several genes .
Our genome-wide significant results for condition factor support the findings of Moghadam et al. , who described a locus on AC-28 that affected this trait. In fact, the same genetic marker that we found to be associated with the QTL for condition factors 1 and 2 on AC-28, was identified by Moghadam et al.  (i.e., Omi34TUF) as localizing the QTL. We also identified a QTL for body weight 1 on AC-28, however according to reference maps [28, 29] the associated marker (i.e., OMM1459) is 49 cM away from Omi34TUF, suggesting that separate QTL regions on this linkage group affect either trait. These regions are clearly autonomous in rainbow trout, for OMM1459 (RT-23, p-arm) and Omi34TUF (RT-1) are located on distinct linkage groups .
QTL detected within the first sampling time period remained constant temporally, indicating that there were no gene × time interaction effects. Further, the QTL effects appeared to be largely unaffected by changes in ambient salinity, given that some QTL were associated with both freshwater (body weight 1) and seawater (body weight 2 and 3) traits, suggesting that gene × environment interaction effects were also minimal. In fact, this pattern was evident among several linkage groups affiliated with body size QTL (i.e., AC-1, -7, -19, -24, -25, -26, -27). This result is somewhat unexpected given that gene × environment interactions have been observed for growth and gene expression in salmonid fishes .
Although uncommon, the effects of some QTL were restricted to either seawater or freshwater. Though tentative observations were made, the discernment of which loci, if any, that could have exerted such effects was difficult given that our experiment was not designed to make such assessments. The only QTL with putative effects specific to seawater exposure was located on AC-9 (body weight 3 and condition factor 3). A QTL on AC-4 also appeared to be seawater-specific (i.e., QTL for only body weight 3 and condition factor 3). However, AC-4 has also been associated with body weight QTL in freshwater Arctic charr , and thus separate loci on AC-4 may influence growth differentially in either rearing environment. Identifying QTL with potential freshwater -specific effects proved more difficult, for only a single freshwater sample was collected. Regardless, our data tentatively suggest that freshwater QTL reside on AC-10 and -17, both of which contain QTL for body weight 1 and condition factor 1. Moghadam et al.  did not find QTL on either linkage group, despite using freshwater Arctic charr in their study.
We were able to confirm several QTL previously detected by other researchers. In our evaluation of body weight and condition factor of individuals in freshwater, we found three body weight 1 QTL (i.e., AC-1, -13, and -25), and four condition factor QTL (i.e., AC-18, -20, -28, -36) that coincided with those of Moghadam et al. . Approximately 60% of the growth-related QTL detected by Moghadam et al.  were also detected in the current study (see Figure 2). This agreement is not surprising, for both studies used the same Arctic charr strain (i.e., Fraser River) and performed genome scans using genetic markers derived from the same reference mapping panel. Though the fish in the present study were 1.5 years old, those used by Moghadam et al.  were only 1 year-old, which could explain some of the discrepancies, as QTL regions are differentially expressed as salmonids age .
This study was the first to use a genome-scan approach to assess the genetic basis of salinity tolerance in Arctic charr. However, it involved a relatively low density genome scan which may encompass 10's to 100's of differentially interacting genes spanning the existing interval regions mapped. Under these constraints, it is perhaps more likely to observe little if any overlap in reported QTL positions rather than the relative congruence that has been reported in QTL studies to date. There is a need for multiple replicated studies examining the association between putative QTL regions and trait expression before a final consensus can be reached as to which regions tend to have the strongest influences upon any given trait. While family based and age-related differences are likely to exist, major strain-specific differences may also exist [10, 14], highlighting the need for additional studies of this important physiological trait.
We identified several genomic regions associated with seawater performance traits. The detection of QTL on AC-12 provided further support for the hypothesis that genetic variation at the ATP1α1b locus may confer an effect on salinity tolerance. In addition, the detection of strong QTL on AC-19 led us to propose that genetic variation at the IGF2 locus may also affect hypo-osmoregulation. Multiple instances of apparent conserved effects among homeologous linkage groups within Arctic charr were observed, and some salinity tolerance QTL appeared to be in regions homologous with salinity tolerance and smoltification QTL in rainbow trout [21, 30]. Body size QTL were not affected by changes in ambient salinity, and coincidence with the findings of others further confirmed multiple freshwater -based body size QTL in Arctic charr. Though these findings are only suggestive, they provide support for the justification of future work in elucidating the genetic basis of salinity tolerance, for it is clear that several regions in the Arctic charr genome affect the salinity tolerance capacity of the individual, thereby providing a foundation for more detailed candidate gene-based experiments.
Six families (denoted 9, 10, 11, 12, 18, and 28) were produced in November, 2006, at the Coastal Zones Research Institute (Shippigan, New Brunswick, Canada) using full-sib crosses from an F3 generation originally derived from an anadromous Arctic charr population from the Fraser River, Labrador, Canada. In the summer of 2007, approximately 150 progeny from each family were PIT (passive integrated transponder) tagged and transferred to St. Andrews Biological Station (St. Andrews, New Brunswick, Canada) where they were reared in two cylindrical 1 m3 freshwater tanks, under controlled simulated-natural photoperiod and water temperature regimes. Tissue was collected from the adipose fin from all individuals for genetic analysis in March, 2008. On June 12, 2008, body weight and fork length measurements were obtained for all individuals (= body weight 1), that were then randomly sorted among six cylindrical 1 m3 tanks to facilitate growth and as a prelude to future experiments. The tanks were equally stocked (~120 fish·tank-1) and families equally represented (~20 fish·family-1tank-1). Each tank was supplied with filtered, aerated freshwater (9.9-10.7°C, flow rate 18 L·min-1, dissolved O2 10.0-10.6 mg·L-1), and covered with a clear Plexiglas lid. To minimize the confounding effects associated with a naturally changing photoperiod, a 16 h-light/8 h-dark photoperiod regime was maintained for the duration of the sampling period (June 9, 2008 to July 6, 2008), after which a simulated-natural photoperiod was restored. Light was provided by 30 W incandescent bulbs, centred approximately 0.75 m above each tank. Feeding occurred daily to satiation with Skretting Optiline salmonid feed pellets (Skretting, Bayside, NB, Canada). Fish were held in these conditions for 7-14 days prior to seawater introduction, dependent upon the tank sampling order.
All tanks were converted to seawater in the period of June 19, 2008 to June 26, 2008, at a rate of one tank·day-1. For each tank, freshwater was replaced with filtered seawater over a 24-hour period: at 6-h intervals freshwater and seawater inputs were proportionately decreased and increased, respectively, such that after 24-h the tank consisted of 100% seawater (31-33‰, 10.5-11.9°C, flow rate 18 L·min-1, dissolved O2 8.1-11.4 mg·L-1). Fish were fasted for 24-h prior to sampling. All sampling was preceded by anaesthetization with tricaine methanesulfonate (MS 222; 150 mg·L-1), and concluded with replacement of the individual into a seawater recovery tank. Ninety-six hours post-full seawater exposure, blood was collected by caudal puncture with a heparinized syringe (500 U mL-1 heparin) and placed on ice. Within 15 minutes of sample acquisition, blood was centrifuged at 13500 g for 4 minutes at 4°C, after which plasma was immediately removed and frozen in liquid nitrogen at -80°C for future analysis. Ten days post full seawater exposure non-lethal gill biopsies  were collected from the same individuals. Upon excision, tissue samples were immersed in 500 μL of ice-cold SEI buffer (250 mM sucrose, 10 mM EDTA, 50 mM imidazole, pH 7.3) in 2 mL cryovials and were frozen within 0.5 h in liquid nitrogen at -80°C for future analysis. We sampled ten days post-seawater exposure as historically this is when differences in Na+/K+-ATPase activity in Arctic charr have been observed . Individuals remained in seawater from June 19-26, 2008, to November 14, 2008. Body weight and fork length measurements were taken again on August 28, 2008 and November 14, 2008, to facilitate the calculation of specific growth rates. At each sampling time, fish were monitored for signs of early maturity, as indicated by secondary sexual characteristics (darkening body colour, kype development) or gamete extrusion. Mortality was monitored daily.
Blood plasma osmolality (mOsmol·kg-1) was measured using a vapour pressure osmometer (Wescor model 5520; Wescor Inc., Utah, USA). Na+/K+-ATPase activity (μmol ADP·mg protein-1·h-1) was determined spectrophotometrically following the methods of McCormick . Gill filaments were homogenized on ice in SEI buffer for 30 s using a disposable pestle grinder system (Fisher Scientific). Homogenates were centrifuged at 5000 g for 30 s at 4°C to separate insoluble material from the supernatant, which was used directly in the assay mixture (189 mM NaCl, 42 mM KCl, 10.5 mM MgCl2, 50 mM imidazole, 0.7 mM ATP, 2.8 mM phospho(enol)pyruvate, 0.22 mM NADH, 4.0 U mL-1 lactic dehydrogenase, 5.0 U mL-1 pyruvate kinase, pH 7.5) or the assay mixture plus ouabain (0.7 mM), a Na+/K+-ATPase enzyme inhibitor. Na+/K+-ATPase activity was measured in triplicate at 340 nm for 10 minutes using a SpectraMax 190 microplate reader (Molecular Devices, Sunnyvale, CA, USA) maintained at 25°C. Protein concentration was determined with a commercial bicinchoninic acid (BCA) protein assay kit (Pierce, Rockford, Illinois, USA). Specific growth rates were calculated for two periods: June 12 to August 28, 2008 (= Specific Growth Rate 1), and August 29, to November 14, 2008 (= Specific Growth Rate 2), with the formula, G = [Ln(Wt) - Ln(Wi)]/t, where Wt was weight at time t, and Wi was the initial weight . Fulton's condition factor was calculated for June 12, 2008 (= Condition Factor 1), August 28, 2008 (= Condition Factor 2), and November 14, 2008 (= Condition Factor 3), with the formula, K = [100 × BW × FL-3], where BW was body weight (g) and FL was fork length (cm). Body weight measurements obtained on August 28, 2008, and November 14, 2008 were designated as body weight 2 and body weight 3, respectively, in addition to the body weight 1 measurement made on June 12, 2008 (as mentioned above).
Families 11, 18, and 28 exhibited high rates of early maturation (46-57%) and thus were not considered for genome scans. Variation coefficients for blood plasma osmolality (Cv = 0.056-0.064) and Na+/K+-ATPase activity (Cv = 0.33-0.40) were similar among the remaining three families (i.e., 9, 10, 12). Genetic maps were created for families 10 (n = 116) and 12 (n = 118), given that they exhibited lower early maturity than family 9, and thus provided a larger number of progeny for QTL analyses. Genetic markers were selected based on pre-existing Arctic charr linkage group assignments [28, 31]. Where possible, markers were chosen at 20 cM intervals to ensure adequate QTL detection power  and comprehensive genome coverage. A standard phenol chloroform protocol was used for genomic DNA . Forward or reverse marker primers were 5'-flourescently end-labeled with tetrachloro-6-carboxy-flourescent (TET) or 6-carboxy-floursecein (FAM). Polymerase chain reaction (PCR) mixtures were made in 7 μL volumes (2.6 ng genomic DNA·μL-1, 1× PCR buffer, 0.125 mM dNTP, 1.5-2.0 mM MgCl2, 0.1 mg·ml-1 BSA, 0.3-0.6 μM of each primer, 0.021 U μL-1 Taq DNA polymerase). PCR conditions began with initial denaturation (95°C for 10 min), followed by 35 cycles of denaturation (95°C for 1 min), annealing (50-58°C for 30 s), and extension (72°C for 30 s), and concluded with final extension (72°C for 5 min). Amplified PCR products were detected using polyacrylamide gel electrophoresis (6% polyacrylamide gel, 19:1 ratio of acrylamide to bisacrylamide, 8 M urea, 0.5× TBE buffer). Prior to loading, PCR products were mixed with 10 μl loading dye (95% formamide, 10 mM NaOH, 0.25% bromophenol blue) and denatured at 95°C for 10 min. Electrophoresis occurred under denaturing conditions in 1× TBE running buffer for 1.5-2.5 h at 1600 V. Gels were scanned using an FMBIO III Fluorescence scanner (MiraiBio Inc., Alameda, CA, USA).
Independent linkage maps were created for males and females due to large differences in recombination rates between the sexes [28, 31, 53]. Linkage of genetic markers and their relative order within linkage groups was established using several modules within the LINKMFEX software package (v2.3; LINKFMEX, LINKGRP, MAPORD, MAPDIS) . Linkage was assigned based on a minimum logarithm of odds (LOD) score of 3.0.
Prior to QTL analysis all salinity tolerance traits were tested for deviation from normality (Kolmogorov-Smirnov and Lilliefors tests). Residual trait values (ANOVA) were used in QTL analyses as tank effects were detected independently in both families for all traits measured (i.e., blood plasma osmolality, Na+/K+-ATPase activity, specific growth rate 1, specific growth rate 2, body weight 1, body weight 2, body weight 3, condition factor 1, condition factor 2, and condition factor 3). Subsequent to the removal of tank effects a secondary effect of body weight 1 on Na+/K+-ATPase activity was detected in each family, therefore the residuals from a linear regression of Na+/K+-ATPase activity with body weight 1 were used in the QTL analysis for Na+/K+-ATPase activity. Statistics were performed with SYSTAT 12 for Windows (SYSTAT Software, Inc., 2007).
QTL analyses were based on 68 individuals per family for physiological QTL (i.e., Na+/K+-ATPase activity and blood plasma osmolality), and 112 individuals per family for body size and growth QTL. Linear regression-based interval analyses were conducted for each trait and parent independently using MultiQTL software (v2.5) . Single-marker analysis was performed across all parents and families combined. Since interval distances among parental maps were often quite variable (i.e., large sex-specific differences in salmonid recombination rates), with random marker positions missing within the different parents used due to chance homozygous genotypes present in these outbred parents, it was not possible to perform the combined family analysis using interval analysis. For all QTL analyses, chromosome-wide LOD thresholds were determined empirically with 1000 permutations of the trait data against the genotypes . Chromosome-wide significant QTL were assigned at a threshold of P ≤ 0.05 (though QTL slightly above this threshold were also acknowledged), and then further tested for genome-wide significance using a B-H False Discovery Rate (FDR) test (α = 0.05).
The authors wish to thank Claude Pelletier of the Coastal Zones Research Institute for generating the family crosses. We also thank Anne Easton and Evan Timusk for their assistance in conducting the salinity tolerance trials, and Steven Leadbeater and Wilfred Young-Lai of the St. Andrews Biological Station for rearing and maintaining the fish. Special thanks to Xia Yue for her assistance in the lab throughout the study. We acknowledge NSERC (Canada) for their funding support.
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.