# Coding Gene SNP Mapping Reveals QTL Linked to Growth and Stress Response in Brook Charr (Salvelinus fontinalis)

Christopher Sauvage, Marie Vagner, Nicolas Derôme, Céline Audet, Louis Bernatchez

► To cite this version:

Christopher Sauvage, Marie Vagner, Nicolas Derôme, Céline Audet, Louis Bernatchez. Coding Gene SNP Mapping Reveals QTL Linked to Growth and Stress Response in Brook Charr (Salvelinus fontinalis). G3, 2012, 2 (6), pp.707 - 720. 10.1534/g3.112.001990. hal-01251103

HAL Id: hal-01251103

https://hal.science/hal-01251103

Submitted on 3 Feb 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Coding Gene SNP Mapping Reveals QTL Linked to Growth and Stress Response in Brook Charr (*Salvelinus fontinalis*)

Christopher Sauvage,*,†,1,2 Marie Vagner,†,8,1 Nicolas Derôme,* Céline Audet,† and Louis Bernatchez*

*Département de Biologie, Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Québec (Québec)
Canada, G1V 0A6, †INRA, UR1052, Unité de Génétique et d'Amélioration des Fruits et Légumes, 84143 Montfavet,
France, ‡Institut des sciences de la mer de Rimouski (ISMER), Université du Québec à Rimouski, Rimouski (Québec)
Canada, G5L 3A1, and §Institut du Littoral et de l'Environnement, LIENS, UMR7266, 17000 La Rochelle, France

**ABSTRACT** Growth performance and reduced stress response are traits of major interest in fish production. Growth and stress-related quantitative trait loci (QTL) have been already identified in several salmonid species, but little effort has been devoted to charrs (genus *Salvelinus*). Moreover, most QTL studies to date focused on one or very few traits, and little investigation has been devoted to QTL identification for gene expression. Here, our objective was to identify QTL for 27 phenotypes related to growth and stress responses in brook charr (*Salvelinus fontinalis*), which is one of the most economically important freshwater aquaculture species in Canada. Phenotypes included 12 growth parameters, six blood and plasma variables, three hepatic variables, and one plasma hormone level as well as the relative expression measurements of five genes of interest linked to growth regulation. QTL analysis relied on a linkage map recently built from *S. fontinalis* consisting of both single-nucleotide polymorphism (SNP, *n* = 266) and microsatellite (*n* =81) markers in an F₂ interstrain hybrid population (*n* = 171). We identified 63 growth-related QTL and four stress-related QTL across 18 of the 40 linkage groups of the brook charr linkage map. Percent variance explained, confidence interval, and allelic QTL effects also were investigated to provide insight into the genetic architecture of growth- and stress-related QTL. QTL related to growth performance and stress response that were identified could be classified into two groups: (1) a group composed of the numerous, small-effect QTL associated with some traits related to growth (i.e., weight) that may be under the control of a large number of genes or pleiotropic genes, and (2) a group of less numerous QTL associated with growth (i.e., gene expression) and with stress-related QTL that display a larger effect, suggesting that these QTL are under the control of a limited number of genes of major effect. This study represents a first step toward the identification of genes potentially linked to phenotypic variation of growth and stress response in brook charr. The ultimate goal is to provide new tools for developing Molecular Assisted Selection for this species.

**KEYWORDS**

linkage mapping
QTL detection
single-nucleotide polymorphism
growth
stress response
*Salvelinus*
*fontinalis*

Quantitative trait loci (QTL), the portions of a species' genome that affect the variation of heritable phenotypic traits, are revealed by the

association of phenotypes with molecular markers. They provide insight into the number of loci affecting a trait and on the distribution of the QTL effects of each locus (Lynch and Walsh 1998). QTL analysis has become a popular method for studying the genetic basis of continuous variation in a variety of systems and is now an integral tool in medical genetics, livestock production, plant breeding, and population genetics of model organisms (Slate 2005). Thus, this approach is of great interest when studying experimental populations obtained from livestock species. Knowledge gained from these analyses can be used to improve traits of economic importance such as growth, resistance to pathogens, sexual maturation, or stress response through the clarification of their underlying genetic basis. Such information helps the development of selection programs aiming to improve efficiency,

Copyright © 2012 Sauvage et al.

doi: 10.1534/g3.112.001990

Manuscript received January 11, 2012; accepted for publication April 16, 2012

This is an open-access article distributed under the terms of the Creative Commons Attribution Unported License (http://creativecommons.org/licenses/by/3.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

<figure><img src="image_2.png" /></figure>

<footnote n="1">¹These authors contributed equally to this work.</footnote>

<footnote n="2">²Corresponding author: INRA, UR1052, Unité de Génétique et d'Amélioration des Fruits et Légumes, Domaine Saint-Maurice BP94, Montfavet 84143, France. E-mail: christopher.sauvage@avignon.inra.fr</footnote>
yield, and production sustainability (Haley and Koning 2006). Indeed, the use of QTL in marker-assisted selection for phenotypic traits of interests offers great potential and provides the basis for cloning genes underlying the genetic architecture of these traits (Mackay 2001; Remington et al. 2001). This approach has been applied in many livestock species but was introduced relatively recently (a decade ago) in farmed aquatic species such as rainbow trout *Oncorhynchus mykiss* (Ozaki et al. 2001), Atlantic salmon *Salmo salar* (Houston et al. 2008), and Pacific oyster *Crassostrea gigas* (Sauvage et al. 2010a) to improve disease resistance or growth. QTL studies provide a framework for the identification of genes and genetic architecture underlying heritable variation within populations and divergence among them. However, this has not proven true from QTL studies alone, which need to be supported by candidate genes approach to fully detect and understand the complex traits architecture (for review, see Rockman 2011).

Growth is one of the most important fitness traits targeted toward a more efficient production of livestock species. The variation of this complex trait relies on a network of genes (e.g., pleiotropy) and on many surrounding environmental inputs (Wang et al. 2002), such as seasonal variations of environmental conditions (Makinen and Ruohonen 1992), food availability (Ali et al. 2003; Bureau et al. 2006), competition (Metcalfe 1986; Blanchet et al. 2007), and other biotic and abiotic factors (e.g., thermal tolerance) (Jackson et al. 1998). Moreover, growth is known to be correlated with variations of other life-history traits, such as gonad maturation processes and reproductive timing (Schaffer 1979; Thorpe 1994; Devlin and Nagahama 2002). Despite the numerous factors influencing growth, in most studies in which investigated its heritability revealed moderate-to-high levels of heritability throughout a wide range of taxa (Wringe et al. 2010).

Stress response, which has been defined as a “diversion of metabolic energy from animal’s normal activities” (Barton and Schreck 1987), is another important fitness-related trait in aquaculture production. In aquaculture facilities, fish are submitted to many stressful manipulations (handling, sorting, transportation, vaccination). All these have the potential to initiate a severe stress response (Barton and Iwama 1991; Portz et al. 2006), which can affect other relevant production traits, including growth performance, feed conversion, immunocompetence, reproductive performance, and disease resistance (Pickering 1981; Adams 1990; Pottinger and Pickering 1997; Wendelaar Bonga 1997; Iversen et al. 1998; Barton 2002).

Salmonids are the most important farmed fish group in Canada. As is the case for other livestock, their growth performance and stress response are of particular economical interest. The mapping of QTL associated with growth traits has been extensively documented in several salmonid species, including rainbow trout (Martyniuk et al. 2003; O’Malley et al. 2003; Perry et al. 2005; Drew et al. 2007; Moghadam et al. 2007a; Wringe et al. 2010), coho salmon *Oncorhynchus kisutch* (McClelland and Naish 2010), Arctic charr *Salvelinus alpinus* (Moghadam et al. 2007b), Atlantic salmon (Reid et al. 2005), and chinook salmon *Oncorhynchus tshawyscha* (Du et al. 1993). The results of these studies have provided insight into the genomic architecture of growth-regulating regions within the salmonid genome. For example, homologous linkage groups with similar QTL effects on fork length and body weight have been observed among different species (O’Malley et al. 2003; Drew et al. 2007; Moghadam et al. 2007b; Wringe et al. 2010). It has also been demonstrated that duplicate copies of growth hormone coding sequences are located in the homologous linkage groups RT-2/9 and that genetic markers close to these regions have been identified as body weight QTL regions in both rainbow trout and Arctic charr (Moghadam et al. 2007b). In addition, recent studies

have reported the identification of QTL and candidate genes related to plasma cortisol concentration in rainbow trout (Drew et al. 2007; Vallejo et al. 2009) as well as three potential QTL related to stress response in sea bass *Dicentrarchus labrax* (Massault et al. 2010). Despite these studies, QTL related to stress response remain poorly studied in fish.

Using brook charr (*Salvelinus fontinalis*), one of the most economically important freshwater aquaculture species in Canada, we aimed to extend the work on salmonids by the identification of QTL underlying two phenotypic traits highly relevant to aquaculture production, i.e., growth performance and stress response. Our analyses were based on a single-nucleotide polymorphism (SNP)-based consensus linkage map (Sauvage et al. 2012) identified by RNA-seq and thus all located in coding genes and a set of 27 traits related to growth and stress response that were phenotyped in 171 F₂ full-sib individuals. These phenotypes included measurements on 12 growth parameters, six blood and plasma variables, three hepatic variables, one stress hormone plasma level, and the expression of five genes of interest related to growth. This study represents a first step toward the identification of genes potentially linked to phenotypic variation of growth and stress response in brook charr. The ultimate goal is to provide new tools for developing molecular-assisted selection for this species.

## MATERIALS AND METHODS

### Biological material and fish crosses

The F₂ population used in this study was obtained from a cross between a domestic population (D) that has been used in aquaculture in Québec (Canada) for more than 100 years and another one (L) that was derived from an anadromous population originating from the Laval River near Forestville (north of the St. Lawrence River, QC, Canada; see Castric and Bernatchez 2003). In previous studies investigators showed that these two strains are highly genetically distinct on the basis of both on gene expression analyses (Bougas et al. 2010) and Fst (The fixation index, Fst is a measure of population differentiation) estimate of 0.187 (± 0.009) on the basis of microsatellite data (Martin et al. 1997). Breeders from the L population were kept in captivity for three generations at the Station aquicole de l’Institut des Sciences de la Mer (ISMER, Pointe-au-Père, QC, Canada, 48°31’N, 68°28’W), whereas those from the D population were obtained from Pisciculture de la Jacques Cartier (Cap-Santé, QC, Canada). In 2005, 10 sires of each population (L and D; F₀ generation) were crossed with 10 dams (L and D) to generate 10 full-sib outbred hybrid (LD) crosses (F₁ generation). The biparental cross of six F₁ individuals then resulted in three F₂ families. At each stage, families were kept separately at the Laboratoire Régional des Sciences Aquatiques (LARSA, Université Laval, QC, Canada) under identical controlled conditions of temperature and photoperiod. Fertilized eggs were incubated at 6°. After hatching, progeny were kept at 8° with a 12L:12D photoperiod. From the three F₂ families, a single F₂ family was used in the present study and chosen according to its lower mortality rate to minimize bias (e.g., segregation distortion) in the subsequent analyses.

### Fish-rearing conditions

Fish from the F₂ hybrid population hatched in January 2008 and were reared in a single indoor flow-through circular tank in fresh water at the Station aquicole de l'ISMER under natural photoperiod and temperature conditions. All fish were tagged with electronic PIT-tags at 5 months for individual identification at all sampling points. Density was maintained below 35 kg of fish per m³. Fish were fed daily with commercial pellets according to fish age and water temperature (from 1.5 to 5.3% of body weight).
## Genetic linkage map

Marker development and construction of the genetic map used for the QTL analysis of this study have been presented in details in a companion study (Sauvage et al. 2012). In brief, a normalized cDNA library was sequenced on a GS-FLX 454 Titanium sequencer. Assembled contigs were screened for SNPs using the implemented tool in CLC Genomic Workbench v3.7 based on Brockman et al. (2008). A subset of 300 SNP markers in the mapping population was validated from the whole set of potential SNP detected in the contigs using a four step validation approach detailed in Sauvage et al. (2012). In addition to SNP markers that were genotyped on F2 progeny, a total of 81 microsatellites markers available in the literature were also used to build the map.

## Data collection

**Sampling and phenotyping for growth:** Samplings were performed on 171 F₂ progeny of age 1+ on four occasions: T1: May 2009 (n = 171; 72.6 ± 1.5 g), T2: July 2009 (n = 171; 145.8 ± 2.5 g), T3: August 2009 (n = 85; 223.9 ± 5.5 g), and T4: November 2009 (n = 86; 273.7 ± 6.4 g). Fish had fasted for 12 hr before each sampling. Fish were captured and immediately placed in an anesthetic solution (3-amino-benzoic-ethyl-ester-acid MS-222, 0.16 g L⁻¹) with constant aeration. They were identified using a PIT-tag reader, measured (± 0.1 cm), and weighed (± 0.1 g).

Two growth-related parameters were calculated for each period:

$$ \text{Specific growth rate (SGR)} = 100 \times (\ln w_{\text{final}} - \ln w_{\text{initial}}) \text{day}^{-1} \quad (1) $$

where SGR is % day⁻¹ and wₗinitial and wₗfinal are the initial and final mean body weights for each period, respectively; and

$$ \text{Fulton index} = 100 \times [(L/w)^3] \quad (2) $$

where L is fish length and w is fish weight.

Blood from the 85 anesthetized fish randomly sampled in August 2009 was sampled by caudal puncture using cooled 1-mL heparinized syringes. Blood samples were used for hematocrit, plasma cortisol, plasma osmolality, and plasma chloride measurements. Fish were killed by decapitation immediately after caudal puncture; this was done according to Canadian Council of Animal Protection recommendations and protocols approved by the University Animal Care Committee. All manipulations were performed quickly so that blood samples were obtained within 2 to 3 min of transfer into the anesthetic solution. Blood samples were immediately centrifuged at 5000 rpm (8500g) for 5 min and collected plasma was immediately frozen at -80° until analyses. Hematocrit (percentage of red blood cells in the centrifuged blood volume) was measured in duplicate in capillary tubes centrifuged for 3 min at 4000 rpm (6800g). Fish were dissected on ice (4°) to collect the whole liver without the gall bladder because it affects RNAs integrity. Livers were weighed, immediately frozen in liquid nitrogen, and stored at -80° for subsequent gene expression measurements. The hepatosomatic index (HSI) was determined for all fish by use of the following relationship:

$$ \text{HSI} = (w_{\text{liver}} \times 100)/(w_{\text{fish}} - w_{\text{liver}}) \quad (3) $$

where wₗliver is the liver weight and wₗfish is the fish weight.

**Quantitative analyses of gene expression:** The relative expression of five genes involved in growth metabolism were measured in

<table><caption>Table 1 Primers used for gene expression analysis by RT-PCR</caption><thead><tr><th>Target Gene</th><th>Primer Set (5' → 3')</th></tr></thead><tbody><tr><td rowspan="2">ghr</td><td>Forward: CCCACTGCCCCCTGTATCT</td></tr><tr><td>Reverse: CTTCAGAAGGAGGCTGTTTTGC</td></tr><tr><td rowspan="2">β-actin</td><td>Forward: GCTGTCTTCCCCTCCATCGT</td></tr><tr><td>Reverse: TCTCCCACGTAGCTGTCTTTCTG</td></tr><tr><td rowspan="2">igf1</td><td>Forward: CAGGCATCCAGATTGTGCAA</td></tr><tr><td>Reverse: ACCATGTTCTGAGAATTCCCTGTGTT</td></tr><tr><td rowspan="2">igf1r</td><td>Forward: AGACCCAGTTTCTGAATTTCACC</td></tr><tr><td>Reverse: GTTCTTATAAGGCGCCTCTTTGT</td></tr><tr><td rowspan="2">ef1</td><td>Forward: GCCCCTCCAGGATGTCTACA</td></tr><tr><td>Reverse: ACGGCCCACGGGTACTG</td></tr><tr><td rowspan="2">18s</td><td>Forward: CCCCGTAATTGGAATGAGTACACTTT</td></tr><tr><td>Reverse: ACGCTATTGGAGCTGGAATTACC</td></tr></tbody></table>

For details, see *Materials and Methods*. RT-PCR, reverse transcription polymerase chain reaction; ghr, growth hormone receptor; β-actin, beta actin; igf1, insulin growth factor-1; igf1r, insulin growth factor-1 receptor; ef1, elongation factor-1; 18s, 18s ribosomal subunit.

liver tissue from the 85 fish sampled in August 2009; these genes were growth hormone receptor (ghr), insulin growth factor-1 (igf1), insulin growth factor-1 receptor (igf1r), elongation factor-1 (ef1), and β-actin. The gene encoding the 18s ribosomal unit was used as a reference gene because its expression was constant between samples. Total RNA was extracted from 30 mg of liver from each fish using RNeasy Plus Mini Kit (QIAGEN, Inc., ON, Canada). Total RNA integrity and quantity were determined using a Nanodrop spectrophotometer (Nanodrop ND-1000 v3.3.0, NanoDrop Technologies, Inc., Wilmington, DE) and a 1.2% agarose gel. To obtain cDNAs, reverse transcription was performed on 1 µg of total RNA in duplicate for each sample using a Quantitect Reverse Transcription kit with integrated removal of genomic DNA contamination (QIAGEN, Inc.). cDNA integrity and quantity were also checked using Nanodrop spectrophotometer. Duplicate cDNAs were then pooled for each sample and real-time polymerase chain reaction (PCR) was performed using the iCycler iQ (Bio-Rad Laboratories Inc., ON, Canada).

The igf1 and ghr primers were designed specifically for S. fontinalis by Côté et al. (2007) using sequence information from *Oncorhynchus kisutch* [GenBank:AF403539] and O. keta [GenBank:AF063216]. The mRNA sequences for igf1r, ef1, β-actin, and 18s were not available for S. fontinalis in the GenBank databases. Consequently, we designed primers from closely related salmonids species by using Primer 3 software (Rozen and Skaletsky 2000) to obtain PCR products ranging from 90 to 150 bp. We used mRNA igf1r sequences from O. mykiss [GenBank: AF062496]; alignment between mRNA β-actin sequences from O. mykiss [GenBank:AF157514] and S. salar [GenBank:NM_001123525]; mRNA ef1 sequences from O. tshawytscha [GenBank:FJ890356]; and mRNA 18s sequence from S. salar [GenBank:AJ427629]. Sequences of primers used are summarized in Table 1.

Amplicons were sequenced to check the specificity of forward and reverse primers: ligation was performed with the TOPO TA Cloning Kit for Sequencing (Invitrogen Inc., ON, Canada) and then transformation was performed using One Shot Chemically Competent E. coli (Invitrogen Inc.). Bacterial cDNA was extracted using EZNA Plasmid Mini Kit I (Omega Bio-Tek, Norcross, GA). Nucleotides were isolated with the Ultra-Step Dye Terminator Removal Kit (Eazy Nucleic Isolation; Ezna, Omega Bio-Tek) and sequenced in forward and reverse sense using the Big Dye Terminator v3 chemistry (Applied Biosystems, Foster City, CA). Alignment between the sequence obtained and the sequence used for primer design was performed for each gene; the similarity percentages obtained were 100% for 18s, 98% for β-actin, 94% for ef1, and 97% for igf1r.
Real-time PCR analyses for each gene were performed in duplicate for each pool of cDNA in a total volume of 15 µL containing 5 µL of cDNA (dilution: 10⁻²), 0.5 µL of primers (10 µmol L⁻¹), and 7.5 µL of 2X iQ SYBR Green Supermix (Bio-Rad Laboratories, Inc.). Thermal cycling of real-time PCR was initiated with incubation at 95° for 13.5 min for activation of the hot-start enzyme, iTaq DNA polymerase. After this initial step, 45 cycles of PCR were performed. Each PCR cycle consisted of heating at 95° for 30 s for denaturing, at 60° for 1 min for annealing, and at 72° for 30 s for extension. Cycle threshold (CT) values corresponded to the number of cycles at which the fluorescence emission monitored in real time exceeded the threshold limit. CT values were automatically calculated on the log curve for each gene. For each plate, a melt curve was established to ensure the presence of a unique amplicon. Thus, the 45 cycles for cDNA amplification were followed by one cycle at 95° for 1 min, one cycle at 55° for 1 min, and 80 cycles at 55° for 10 sec. Standard curves were established in triplicate for each gene by plotting the CT values against log₁₀ of five different dilutions (in triplicate) of a pool of all cDNA sample solutions. Real-time PCR efficiency (E) was determined for each gene from the slope of the mean standard curve according to the equation (4):

$$E = 10^{[-1/\text{slope}]} \quad (4)$$

To determine the relative quantity of target gene-specific transcripts present in the different samples, relative expression ratios were calculated according to equation (5). The relative expression ratio for a considered gene is based on the PCR efficiency (E) and the CT of a sample vs. the control (standard group); it is expressed in comparison to the reference gene (18s) according to Pfaffl's equation (Pfaffl 2001):

$$\text{Ratio} = \left[ \left( E_{\text{gene}} \right)^{\Delta \text{CT}_{\text{gene}(\text{control-sample})}} \right] / \left[ \left( E_{18s} \right)^{\Delta \text{CT}_{18s}(\text{control-sample})} \right] \quad (5)$$

The mean CT of samples from standard curves and diluted at 10⁻² was used as the standard group because it was the most representative cDNA of the population tested. Normalization was used to correct for intensity distortions as well as spatial variation in signal level across the different real-time PCR runs.

**Physiological measurements:** Plasma glucose concentrations and hepatic glycogen levels were measured in all 85 fish sampled in August 2009. Plasma glucose concentrations (mg·mL plasma⁻¹) were measured using a commercial kit (QuantiChrom Glucose Assay kit; BioAssay Systems, Hayward, CA). Hepatic glycogen levels (mg·g of liver⁻¹) were assessed using the amyloglucosidase digestion method (Carr and Neff 1984) followed by glucose concentration determinations (Quanti-Chrom Glucose Assay kit; BioAssay Systems). All measurements were done in duplicates, which all varied by less than 5% (data not shown).

**Sampling and phenotyping for stress:** In November 2009, the baseline of the stress response for 86 F₂ progeny was determined. Water level was gradually lowered in the single circular tank, and fish were individually captured without being pursued. They were individually anesthetized (MS-222; 0.16 g L⁻¹), identified using a PIT-tag reader, and then blood was sampled (0.3 mL) by caudal puncture using cooled heparinized syringes. After a full recovery, fish were returned to their initial tank.

One week later, the stress experiment was performed on the 86 fish. Water level was again gradually lowered in the single circular tank, and fish were individually captured without being pursued. After being captured, each fish was individually stressed by one minute of

handling out of water in a small net. After the handling stress, fish were placed in groups of six in smaller opaque tanks (volume: 80 L; fish density: 1.95 10⁻³ kg L⁻¹) until blood sampling. Three hours after the handling stress (Bastien 2010), fish were anesthetized (MS-222; 0.16 g L⁻¹). Fish were identified using a PIT-tag reader, weighed (± 0.1 g), and measured (± 0.1 cm) before caudal puncture. Fish were killed by decapitation immediately after caudal puncture; this was done according to Canadian Council of Animal Protection recommendations and protocols approved by the University Animal Care Committee. Sex and maturation status were visually determined by inspection of the gonads. Fish were classified as mature (when eggs and sperm could be collected by stripping) or immature.

All samplings (before and after acute stress) were done between 1:00 PM and 2:00 PM to avoid circadian variation effects on plasma cortisol concentration (Audet and Claireaux 1992). All manipulations were done quickly so that blood was obtained within 2 to 3 min after transfer to the anesthetic solution. Plasma was obtained by centrifugation (5 min, 5000 rpm; 8500g) and then stored at -80° until analysis.

**Physiological measurements of stress response:** Plasma cortisol concentrations (µg dL⁻¹ of plasma) were measured for each fish using a commercial radioimmunoassay kit (ImmuChem Cortisol ¹²⁵I RIA kit; MP Biomedicals, Cleveland, OH) validated in fish (Vijayan and Moon 1992). Radioactivity of the [¹²⁵I]-labeled cortisol tracer was quantified using the automatic CliniGamma 1272 gamma counter (LKB-Wallac, Wallac, Finland). For each assay, a standard curve was constructed and cortisol levels were within the linear range of the assay. All measurements were done in duplicate, and replicate measurements varied by less than 5%.

Plasma osmolality (mmol kg⁻¹) was measured using a Vapro Vapor Pressure Osmometer 5520 (Wescor Inc., Logan, UT) and plasma chloride (mmol L⁻¹) was measured using the Chloride Analyzer 925 (Corning Medical and Scientific, England). All measurements were performed in duplicate, and replicate measurements varied by less than 5%.

**Phenotype data analyses:** Because some traits were measured at several sampling times (i.e., size and weight were measured at T1, T2, and T3 on 85 of the F₂ progeny and again at T4 on the remaining 86 fish), we decided to perform the QTL analyses on the largest time interval (from T1 to T3 for the first 85 sampled fish and from T1 to T4 for the last 86 sampled fish) to maximize QTL detection and variations in phenotype measurements. QTL detection related to size, weight, SGR, and Fulton index were then performed with the data set obtained at T4 (November), whereas the QTL detection for traits related to blood parameters (hematocrit, plasma chloride, plasma osmolality, and plasma glucose), gene expression (ghr, igf1, igf1R, efl, and β-actin), and liver variables (liver fresh weight, hepato-somatic index, and hepatic glycogen) were performed on the data set obtained at T3 (August).

Correlations within phenotypic traits related to growth performance (collected at T1, T2, T3, and T4) and stress response (collected at T4) were determined by correlation matrix using Statistica 6 for Windows (StatSoft Inc., Tulsa, OK). Normality was tested on residuals using Kolmogorov-Smirnov test. Phenotypic traits that were not normally distributed were log₁₀ transformed. Another assumption for correlation matrix is the absence of outliers. Thus, outliers (n = 10) were identified and removed from the database using the scatterplot box-plot function (Statistica). Differences were considered significant at *P* < 0.05. The graphical correlation matrix displaying the positive and negative correlations between the phenotypes related to growth and the phenotypes related to the stress response were obtained using the R package "Corrplot" (Friendly 2002).
**QTL detection:** QTL analyses for the aforementioned 27 phenotypic traits were performed using the [R] package R/qlt (v. 1.18-7, August 2010, http://www.rqtl.org/) (Broman *et al.* 2003) only on the sex-averaged (consensus) linkage map. The following approach was used: (1) A single QTL analysis was performed using the Haley-Knott (HK) regression method (10 000 permutations) (Haley and Knott 1992) to reveal which linkage groups (LGs) were carrying QTL. The most probable position of the QTL was defined at the position giving the largest log₁₀ of the odd ratio (LOD) score; this QTL was fixed. (2) Then, a full model that incorporates all those loci identified in the single QTL scan was used to refine the positions and to estimate effects and PVE across the genome with a resolution of 5 cM. (3) Finally, the model best fitting our data were used to compute the percent variance explained (PVE) associated with the QTL. The chromosome-wide and genome-wide thresholds were calculated for each LG using 10,000 permutations. The 1.5 LOD confidence intervals were determined for all analyses following the Bayesian method implemented in the "bayesint" function in R/qlt. The bayesint function calculates an approximate interval (end points around the maximum LOD) for a given chromosome using the genome scan output. Allele effects were determined using the effect plot function in R/qlt with the QTL peak marker or the marker nearest to the peak as the reference marker. The additive effect was estimated as one-half of the difference between the two homozygous genotype values. The dominance effect was estimated as the deviation of the heterozygous genotype values from the average of the two homozygous genotypic values (Lynch and Walsh 1998). The resulting value indicated that the progeny should closely resemble one of the two parental lines rather than having an intermediate phenotype (Bennewitz and Meuwissen 2010).

# RESULTS

## Genetic linkage map

The detailed building of the linkage map is described in Sauvage *et al.* (2012). Briefly, the dataset used to build the linkage map comprised 81 microsatellite markers and 256 SNP located in coding gene regions, for a total of 337 markers. Forty LGs were generated, which is close to the haploid number of chromosomes in brook charr (2n = 84). The consensus (sex-averaged) map contained 266 markers (191 SNPs and 75 microsatellites) displayed in the 40 LGs (see Table 2 in Sauvage *et al.* 2012 for details). The LG length ranged from 1.4 cM to 132 cM, for a total map length of 2047.5 cM. The average marker spacing per LG ranged from 0.7 to 21.3 cM and was estimated at 8.3 cM over the whole genome. The genome coverage was estimated at 89% as following: a complete female map is usually expected to be roughly 25 Morgans, assuming approximately 50 cM per chromosome arm pair (50 chromosome arm sets in the brook charr). In the present study, the female map that covers a total of 22.48 Morgans represents about 89% of the genome. The exact position and order of the 266 markers among the 40 LGs are given in supplementary Table 2 in Sauvage *et al.* (2012).

## Phenotyping

The detailed statistical descriptions of the 27 phenotypes analyzed are reported in Table 2 and Table 3. Among the phenotypes measured for growth related trait, standard deviations (SD) of the mean varied highly, from 0.07 for the Fulton's condition factor at T3 up to 14.53 for hepatic glycogen. SD for SGR T2-T3 (0.24) was almost 2.5-fold greater than that for SGR T1-T3 (0.10). SD for stress-related phenotypes also varied highly, from 5.81 for the plasma chloride difference between before and after acute stress to 17.26 for plasma osmolality

difference between, before, and after acute stress. The hierarchical clustering of phenotypes on the basis of correlation revealed two main clusters. The first is composed of six closely related phenotypes; size T4, weight T4, SGR T1-T4, Fulton's condition factor T4, HSI T3, and liver weight T3 (Figure 1). The second cluster is composed of gene expression (mRNA levels of *ef1*, *igf1*, *igf1r*, *ghr*, and *β-actin*; Figure 1). Phenotypes related to stress response and measured at T4 were expressed as the difference in levels before and after acute stress (response intensity). A negative correlation (r = -0.25; P = 0.021) was found between plasma chloride and plasma cortisol response intensity, while a positive correlation (r = +0.34; P = 0.001) was measured between plasma cortisol and plasma osmolality response intensities.

# QTL detection

**Traits related to growth:** A total of 18 of the 40 LGs comprised QTL related to one of the phenotypic traits. For growth, 63 QTL were identified across the 22 growth-related phenotypes, with 8 QTL at the genome-wide level of significance and 55 at the chromosome-wide level. The greatest number of QTLs identified for a single trait was estimated at 12 for individual body weight (T4) and specific growth rate (SGR T1-T4), whereas the lowest number was estimated at one for blood hematocrit, plasma osmolality, and plasma glucose phenotypes as well as for the expression of *igf1*, *igf1r*, and *ghr* genes. LG22 carried the greatest number of QTL, with a total of six QTL related to six different phenotypes (size T4, weight T4, SGR T1-T3, SGR T1-T4, HSI T3, and liver weight T3). The minimal, maximal, and average LOD scores associated with the 63 QTL were estimated at 3.01 (QTL for Fulton index T4), 15.99 (QTL for HSI T3), and 8.29, respectively (Table 4). The respective minimal, maximal, and average PVE by the identified QTL related to growth were estimated at 3.08% (QTL linked to SGR T1-T4 on LG 3 at 44.8 cM), 17.27% (QTL linked to hematocrit on LG 4), and 6.49%, respectively (Table 4). Some QTL linked to different traits were colocalized on a same linkage group. Thus, QTL related to size, weight, and SGR were colocalized on LGs 3, 12, 20, 22, 24, and 32 and QTL linked to HSI and liver weight traits were co-localized on LGs 20 and 22.

For QTL related to growth, the distribution of additive and dominance effects estimated for each of the QTL was relatively large, ranging from -16.17 ± 0.45 (QTL for mRNA *ghr* relative expression) to 44.32 ± 8.36 (g) (mean = 5.58 ± 0.88) for additive effects and from -40.24 ± 0.01 g to 34.58 ± 8.27 g (mean = 3.35 ± 0.66) for dominance effects. Moreover, only one QTL showed an additive effect below 10, and 13 QTL showed additive effects above 10. Two QTL showed dominance effects below 10, and 10 QTL showed dominance effects above 10. For the 63 QTL identified for growth related traits, 30 were associated with an additive gene effect and 20 of these 30 displayed a positive effect. The 33 remaining QTL were associated with a dominance gene effect, and 11 of these 33 QTL displayed a positive effect. Only 12 QTL were linked to SNP markers in contrast with 51 for microsatellite markers. Within these 12 SNP markers, five were identified in coding regions, and two of these were characterized as nonsynonymous changes (Table 6). Seven of these 12 SNP were characterized as transitions, and five were annotated (Table 6). Finally, no significant QTL were detected for five of the 24 (20.8%; 5/24) growth-related phenotypes. These were related to specific growth rate (SGR T1-T2), Fulton's condition factor measured in May (T1) and July (T2), and to level of expression of the *ef1* and *β-actin* genes.

**Traits related to the stress response:** Four QTL were identified for the three traits related to the stress response: two QTL were related

<figure><img src="image_3.png" /></figure>
<table><caption>Table 2 Descriptive statistics of phenotypic traits related to growth measured in sampled fish</caption><thead><tr><th rowspan="2"></th><th colspan="4">Specific Growth Rate, % d<sup>-1</sup></th><th colspan="4">Length, cm</th></tr><tr><th>T1-T2</th><th>T2-T3</th><th>T1-T3</th><th>T1-T4</th><th>T1</th><th>T2</th><th>T3</th><th>T4</th></tr></thead><tbody><tr><td>N</td><td>171</td><td>171</td><td>85</td><td>86</td><td>171</td><td>171</td><td>85</td><td>86</td></tr><tr><td>Min</td><td>-0.12</td><td>0.27</td><td>0.78</td><td>0.51</td><td>15.60</td><td>18.00</td><td>22.00</td><td>24.60</td></tr><tr><td>Max</td><td>1.53</td><td>2.16</td><td>1.20</td><td>0.89</td><td>23.40</td><td>31.60</td><td>31.60</td><td>32.90</td></tr><tr><td>Mean ± SD</td><td>1.16 ± 0.17</td><td>0.63 ± 0.24</td><td>0.98 ± 0.10</td><td>0.68 ± 0.08</td><td>19.14 ± 1.61</td><td>23.09 ± 1.92</td><td>27.04 ± 1.86</td><td>28.75 ± 1.75</td></tr><tr><th rowspan="2"></th><th colspan="4">Fulton's Condition Factor</th><th colspan="4">Blood and Plasma Variables</th></tr><tr><th>T1</th><th>T2</th><th>T3</th><th>T4</th><th>Hematocrit, % of red cells</th><th>Chloride, mmol L<sup>-1</sup></th><th>Osmolality, mmol kg<sup>-1</sup></th><th>Glucose, mg mL<sup>-1</sup></th></tr><tr><td>N</td><td>171</td><td>171</td><td>85</td><td>86</td><td>81</td><td>85</td><td>85</td><td>85</td></tr><tr><td>Min</td><td>0.62</td><td>0.37</td><td>0.97</td><td>0.82</td><td>28.00</td><td>111.93</td><td>301.50</td><td>0.35</td></tr><tr><td>Max</td><td>1.36</td><td>1.49</td><td>1.53</td><td>1.50</td><td>46.00</td><td>155.00</td><td>350.50</td><td>1.22</td></tr><tr><td>Mean ± SD</td><td>1.02 ± 0.12</td><td>1.18 ± 0.13</td><td>1.11 ± 0.07</td><td>1.14 ± 0.1</td><td>36.60 ± 3.72</td><td>133.58 ± 5.09</td><td>319.87 ± 10.21</td><td>0.67 ± 0.18</td></tr><tr><th rowspan="2"></th><th>Liver Fresh Weight, g</th><th>Hepato-somatic Index, HSI</th><th>Hepatic glycogen, mg g<sup>-1</sup></th><th colspan="5">mRNA Level in Liver (With Respect to 18s mRNA Level)</th></tr><tr><th></th><th></th><th></th><th>ghr</th><th>igf1</th><th>igf1r</th><th>ef1</th><th>β-actin</th></tr><tr><td>N</td><td>171</td><td>171</td><td>85</td><td>83</td><td>83</td><td>83</td><td>82</td><td>83</td></tr><tr><td>Min</td><td>0.95</td><td>0.35</td><td>65.39</td><td>0.003</td><td>0.002</td><td>0.001</td><td>0.0003</td><td>0.002</td></tr><tr><td>Max</td><td>10.35</td><td>3.34</td><td>146.32</td><td>1.41</td><td>1.13</td><td>0.57</td><td>0.73</td><td>1.04</td></tr><tr><td>Mean ± SD</td><td>3.65 ± 1.86</td><td>1.57 ± 0.82</td><td>109.42 ± 14.53</td><td>0.16 ± 0.24</td><td>0.09 ± 0.17</td><td>0.07 ± 0.10</td><td>0.07 ± 0.12</td><td>0.09 ± 0.15</td></tr></tbody></table>

N, number of observations, Min, minimum value obtained; Max, maximal value obtained, Mean ± SD: mean ± SD of the values measured. T1, May; T2, July; T3, August; T4, November; 18s, 18s ribosomal subunit; ghr, growth hormone receptor; igf1, insulin growth factor-1; igf1r, insulin growth factor-1 receptor; ef1, elongation factor-1; β-actin, beta-actin.

to plasma cortisol, whereas one significant QTL was found for plasma chloride and another for osmolality. The QTL for cortisol level on LG 23 was significant at the genome-wide level, whereas the three others were significant at the chromosome-wide level (QTL for plasma chloride on LG7, QTL for osmolality on LG10, and the other QTL for plasma cortisol on LG 14). The minimal, maximal, and average LOD scores associated with the four QTL were estimated at 3.21, 8.16, and 4.64, respectively (Table 5). The minimal, maximal, and average PVE were estimated at 3.85% (QTL for plasma cortisol on LG7), 31.35% (QTL for the level of cortisol on LG23), and 13.79%, respectively (Table 5). The distribution of additive (from -2.12 µg dL⁻¹ ± 0.59 for the plasma osmolality QTL to 0.63 mmol L⁻¹ ± 0.49 for the plasma chloride QTL) and dominance (from -3.98 mmol kg⁻¹ ± 0.53 for the plasma osmolality QTL to 8.87 µg dL⁻¹ ± 2.53 for the plasma cortisol QTL) effects estimated for each of the four QTL related to stress response was relatively weak (Table 5). Of these four QTL, only one was associated with a positive additive gene effect while the three others were associated with a dominance effect (positive for two and negative for one). Among the four QTL, two were linked to SNP markers (sf003382, sf004319) and two to microsatellite markers (SAL5UoG and OMM-5312 were linked to the QTL associated with plasma osmolality and plasma chloride, respectively). Both SNP markers were identified in coding regions; one was characterized as

a transition (sf003382), the other as a transversion (sf004319) annotated as a hydroxymethylglutaryl-coA lyase (HMG-CoA, sf004319) (Table 6). Figure 2 displays the QTL identified over all LG for the traits associated with growth and stress response and their respective significance levels.

## DISCUSSION

### Growth-related QTL

This study allowed the identification of QTL for a total of 27 phenotypes related to two physiological traits of economical importance in *S. fontinalis*, one of the most economically important salmonid species for freshwater aquaculture. First, this allowed the identification of 63 growth-related QTL. The low number of QTL identified genome-wide (only eight) compared with the high number identified chromosome-wide (55) is in accordance with previous studies in salmonids in which authors identified QTL related to body weight and condition factor in Arctic charr *S. alpinus* (Kuttner *et al.* 2011) or QTL related to growth in domestic and wild rainbow trout *Onchorynchus mykiss* (Wringe *et al.* 2010). The limited number of molecular markers and sample size used in our study may partly explain these observations because these factors tend to overestimate the magnitude of large-effect QTL and hide the numerous small-effect

<table><caption>Table 3 Descriptive statistics of phenotypic traits related to stress measured in the plasma of sampled fish in November 2009 (T4)</caption><thead><tr><th></th><th>Cortisol, µg dL<sup>-1</sup><br>Difference between before<br>and after acute stress</th><th>Osmolality, mmol kg<sup>-1</sup><br>Difference between before<br>and after acute stress</th><th>Chloride, mmol L<sup>-1</sup><br>Difference between before<br>and after acute stress</th></tr></thead><tbody><tr><td>N</td><td>86</td><td>86</td><td>86</td></tr><tr><td>Min</td><td>-8.39</td><td>-68.00</td><td>-15.00</td></tr><tr><td>Max</td><td>29.44</td><td>50.00</td><td>21.00</td></tr><tr><td>Mean ± SD</td><td>5.02 ± 6.63</td><td>3.16 ± 17.26</td><td>-1.47 ± 5.81</td></tr></tbody></table>

N, number of fish tested; Min, minimum value obtained; Max, maximal value obtained; Mean ± SD, mean ± SD of the values measured.
<figure><img src="image_4.png" /><figcaption>**Figure 1** Correlation matrix displaying the positive (blue) and negative (red) correlations between the phenotypes related to growth measured at T3. Phenotypes have been clustered according to their correlation factor. The graphical representation was given by the R package "Corrplot" (Friendly 2002). SGR, specific growth rate; HSI, hepato-somatic index; *igf1*, insulin growth factor-1; *ef1*, elongation factor-1; *β-actin*, *beta-actin*; *ghr*, growth hormone receptor; *igf1r*, insulin growth factor-1 receptor. Color and size variation of circles in the figure are proportional level and direction (positive or negative) of the correlation between traits.</figcaption></figure>

QTL (Mackay *et al.* 2009). Note, however, Lander and Kruglyak (1995) strongly recommend the use of a genome-wide significance threshold, regardless of the density of the map used, to further deploy many additional markers in regions that show signs of a segregating QTL after an initial sparse map search (Van Ooijen 1999). Similarly, although the low number of QTL identified per trait may indicate that the expression of these traits is under the control of one major or a limited number of genes, we cannot refute the possibility that it may also reflect the limited lack of power of the study. Namely, the limited number of progeny may influence QTL detection experiments by underestimating the QTL number and overestimating their respective effects (Beavis 1998; Xu 2003). Moreover, the fact that the number of detected QTL varied substantially depending on traits suggests that power limitation is not the only factor explaining the detection of one or very few QTL for some of these, especially for polygenic traits (i.e., genetic architecture) for which, individual locus contribution (i.e., distributions of effects) is relatively low and remain difficult to detect.

The 63 QTL associated with the growth-related traits detected in this study is comparable with the 64 and 31 growth-related QTL identified by Kuttner *et al.* (2011) and Wringe *et al.* (2010) in Arctic charr and rainbow trout, respectively, on the basis of maps comprising 106 and 137 markers, respectively. Admittedly, some assumptions made for the analysis may have contributed to overestimates in the number of detected QTL. Namely, we assumed that each QTL corresponded to an independent repeated observation. This assumption was a simplification because some QTL were found at very similar

chromosomal positions (i.e., colocalization) and thus might be the same (i.e., growth-related QTL). In addition, some traits are correlated and hence QTL for two correlated traits at a similar chromosomal position might in fact be a single pleiotropic QTL (i.e., because of the effect of a single gene on multiple phenotypic traits). Jansen *et al.* (1995) proposed statistical methods to test for the pleiotropic effect of genes. However, the application of these methods remains challenging when using low- to medium-density genetic linkage map, as was the case here and in the majority of QTL analysis in salmonids thus far. This is also particularly problematic for salmonids in general, given the lower recombination rate in males that leads to a restricted resolution of the male genetic map (Moghadam *et al.* 2007b). Nevertheless, estimations of QTL positions at the genome scale show relatively large 95% confidence interval (95% CI) (Visscher *et al.* 1996) and hence are not estimated with enough precision to claim that the QTL positions are a single pleiotropic QTL (Bennewitz and Meuwissen 2010). The most probable locations of the 63 growth-related QTL were identified with varying degrees of accuracy, with minimal, maximal, average, and median 95% CI of 2.4, 66, 19.7, and 12 cM, respectively. Again, these values are comparable with those previously reported in Arctic charr (Kuttner *et al.* 2011) and rainbow trout (Wringe *et al.* 2010).

The number of QTL detected and their associated PVE for each trait also provide further insight into the genetic architecture underlying the variance of the phenotypic traits studied. Phenotypes such as weight, length, SGR, HSI, and liver fresh weight typically characterizing growth traits generally showed numerous QTL (ranging from four to nine) associated with small PVE (ranging from 3.07 to 9.37). This finding suggests that these phenotypes are under the control of a large number of genes of relatively small effects. On the contrary, hematocrit, plasma chloride, plasma glucose, plasma osmolality, hepatic glycogen, and mRNA levels of *ghr*, *igf1*, and *igf1r* showed few QTL (<2) with generally larger PVE (up to 17.27%). This finding suggests that these phenotypes are under the control of one or two major genes of larger effects. In total, however, only six QTL (linked to hematocrit, plasma glucose, hepatic glycogen, and mRNA levels of *ghr*, *igf1*, and *igf1r*) showed PVE over the 10% threshold and can be characterized as large-effect QTL.

Although correlations between microsatellite markers and the variation of traits related to growth and reproduction were also investigated in other salmonids, such as rainbow trout (Wringe *et al.* 2010) and Arctic charr (Kuttner *et al.* 2011), it was not possible to compare the colocalization of these previously identified QTL with those identified in the present study because of the low number of microsatellite marker pairs shared between these studies.

## Identification of SNP markers linked to growth-related QTL

For each QTL position, the closest underlying molecular marker was identified. Of the 63 growth-related QTL, only 12 were linked to a SNP marker, whereas the 51 others were linked to a microsatellite marker. Moreover, only five of these 12 were identified in coding regions; two of these five were characterized as nonsynonymous changes, and five of the 12 (41.6%; 5/12) had a significant annotation. The seven remaining SNP were identified in UTR or noncoding regions. Thus, direct assumptions of the role played by the genes linked to QTL in the genetic architecture of the phenotypic traits remain difficult, especially since the molecular function of these five genes was not directly related to growth, considering their annotation. We thus suggest that these SNP markers are in linkage disequilibrium
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /><figcaption>Table 4 Descriptive statistics of LG, including the position, 95% CI, LOD score, PVE, associated P value, and specific additive and dominance effects of each QTL linked to every phenotype related to growth</figcaption></figure>

<table><thead><tr><th>Phenotype<sup>a</sup></th><th>LG</th><th>Position, cM</th><th>95% CI, cM</th><th>LOD</th><th>% PVE</th><th>P Value, χ<sup>2</sup></th><th>P Value, F</th><th>Additive Effect</th><th>SE</th><th>Dominance Effect</th><th>SE</th><th>Nearest Marker</th></tr></thead><tbody><tr><td rowspan="12">Weight T4, g</td><td>1</td><td>70.5</td><td>70-83</td><td>4.742</td><td>4.477</td><td>0.032</td><td>0.105</td><td>19.696</td><td>0.2342</td><td>29.897</td><td>0.1300</td><td>OMM-1195</td></tr><tr><td>3</td><td>44.8</td><td>20-48</td><td>4.945</td><td>3.255</td><td>0.009</td><td>0.046</td><td>22.488</td><td>0.2365</td><td>25.514</td><td>0.0845</td><td>OMM-5007</td></tr><tr><td>4</td><td>43.0</td><td>39-47</td><td>7.438</td><td>5.386</td><td>0</td><td>0.005</td><td>28.711</td><td>0.8577</td><td>12.158</td><td>0.4404</td><td>OMM5155i</td></tr><tr><td>5</td><td>26.7</td><td>21-28</td><td>5.267</td><td>8.826</td><td>0.197</td><td>0.352</td><td>10.377</td><td>0.1783</td><td>14.236</td><td>0.0952</td><td>CA-378164</td></tr><tr><td>9</td><td>40.2</td><td>39-43</td><td>5.137</td><td>4.009</td><td>0.303</td><td>0.468</td><td>23.987</td><td>0.1054</td><td>34.588</td><td>0.2607</td><td>SSA-0072BSFU</td></tr><tr><td>12</td><td>78.4</td><td>73-78.4</td><td>8.803</td><td>6.783</td><td>0.203</td><td>0.359</td><td>3.5396</td><td>0.0340</td><td>-40.241</td><td>0.0346</td><td>OMM-5060</td></tr><tr><td>22</td><td>17.1</td><td>14-19</td><td>5.137</td><td>5.841</td><td>0.087</td><td>0.204</td><td>44.324</td><td>0.2342</td><td>5.8255</td><td>0.1300</td><td>OMM-3015i</td></tr><tr><td>24</td><td>11.9</td><td>41-94</td><td>8.583</td><td>6.433</td><td>0.162</td><td>0.309</td><td>8.9338</td><td>0.2365</td><td>31.112</td><td>0.0845</td><td>OMM-1201</td></tr><tr><td>27</td><td>0.0</td><td>0-6</td><td>5.858</td><td>6.027</td><td>0.003</td><td>0.02</td><td>20.090</td><td>0.8577</td><td>19.243</td><td>0.4404</td><td>OMM-5102ii</td></tr><tr><td>31</td><td>0.0</td><td>0-5</td><td>7.293</td><td>4.866</td><td>0.083</td><td>0.199</td><td>-4.786</td><td>0.1783</td><td>12.887</td><td>0.0952</td><td>C129</td></tr><tr><td>34</td><td>0.0</td><td>0-12</td><td>5.185</td><td>5.521</td><td>0.142</td><td>0.283</td><td>9.2407</td><td>0.1054</td><td>23.328</td><td>0.2607</td><td>BX319411i</td></tr><tr><td>40</td><td>52.4</td><td>46-90</td><td>5.633</td><td>4.745</td><td>0.021</td><td>0.079</td><td>17.955</td><td>0.0340</td><td>-0.681</td><td>0.0346</td><td>BX-079862</td></tr><tr><td>12</td><td>78.4</td><td>76-78.4</td><td>4.276</td><td>4.669</td><td>0.036</td><td>0.051</td><td>-0.016</td><td>0.2365</td><td>0.2152</td><td>0.0845</td><td>OMM-5060</td></tr><tr><td>20</td><td>0.0</td><td>0-37</td><td>4.405</td><td>3.984</td><td>0.041</td><td>0.057</td><td>0.0519</td><td>0.8577</td><td>0.1643</td><td>0.4404</td><td>OMM1210</td></tr><tr><td>22</td><td>17.1</td><td>0-17</td><td>3.227</td><td>5.515</td><td>0.048</td><td>0.056</td><td>-0.306</td><td>0.1783</td><td>-0.089</td><td>0.0952</td><td>OMM-3015i</td></tr><tr><td>34</td><td>0.0</td><td>0-5</td><td>3.605</td><td>4.474</td><td>0.042</td><td>0.051</td><td>-0.166</td><td>0.1054</td><td>-0.151</td><td>0.2607</td><td>BX319411i</td></tr><tr><td rowspan="12">Specific growth rate, T1-T3, % day<sup>-1</sup></td><td>1</td><td>70.5</td><td>70-74.5</td><td>14.275</td><td>4.587</td><td>0</td><td>0.009</td><td>-0.171</td><td>0.0340</td><td>-0.158</td><td>0.0346</td><td>OMM-1195</td></tr><tr><td>2</td><td>78.2</td><td>78-115</td><td>8.016</td><td>4.119</td><td>0</td><td>0.022</td><td>0.1342</td><td>0.2342</td><td>0.0340</td><td>0.1300</td><td>OMY-RGT2TUfii</td></tr><tr><td>3</td><td>44.8</td><td>10-46</td><td>8.509</td><td>3.083</td><td>0</td><td>0.003</td><td>-0.121</td><td>0.0371</td><td>-0.112</td><td>0.0407</td><td>OMM-5007</td></tr><tr><td>9</td><td>40.2</td><td>39-42</td><td>10.781</td><td>5.022</td><td>0.028</td><td>0.019</td><td>-0.175</td><td>0.0371</td><td>-0.128</td><td>0.0370</td><td>SSA-0072BSFU</td></tr><tr><td>12</td><td>78.4</td><td>76-78.4</td><td>11.862</td><td>4.628</td><td>0.003</td><td>0.049</td><td>-0.039</td><td>0.1488</td><td>0.1737</td><td>0.0532</td><td>OMM-5060</td></tr><tr><td>14</td><td>45.9</td><td>36-47</td><td>6.128</td><td>9.376</td><td>0.039</td><td>0.225</td><td>0.1685</td><td>0.1574</td><td>0.0419</td><td>0.0829</td><td>OMI-30TUfii</td></tr><tr><td>20</td><td>0.0</td><td>0-38</td><td>10.361</td><td>3.449</td><td>0</td><td>0.012</td><td>0.0416</td><td>0.5659</td><td>0.1424</td><td>0.2905</td><td>OMM1210</td></tr><tr><td>22</td><td>17.1</td><td>0-17.1</td><td>9.466</td><td>3.468</td><td>0.005</td><td>0.08</td><td>-0.216</td><td>0.1228</td><td>-0.079</td><td>0.0656</td><td>OMM-3015i</td></tr><tr><td>23</td><td>13.1</td><td>9-25</td><td>8.034</td><td>5.03</td><td>0.027</td><td>0.018</td><td>0.1479</td><td>0.1771</td><td>0.1234</td><td>0.0891</td><td>BHMS-238</td></tr><tr><td>24</td><td>11.9</td><td>9-13.5</td><td>13.409</td><td>4.148</td><td>0</td><td>0.021</td><td>-0.016</td><td>0.0343</td><td>-0.175</td><td>0.0334</td><td>OMM-1201</td></tr><tr><td>31</td><td>42.4</td><td>36-42.4</td><td>13.392</td><td>4.122</td><td>0.019</td><td>0.015</td><td>0.1125</td><td>0.3519</td><td>0.0904</td><td>0.4061</td><td>sf004732_01CG</td></tr><tr><td>34</td><td>0.0</td><td>0-5</td><td>12.315</td><td>4.175</td><td>0</td><td>0.02</td><td>-0.109</td><td>0.0753</td><td>-0.085</td><td>0.1861</td><td>BX319411i</td></tr><tr><td rowspan="12">Size T4, mm</td><td>3</td><td>49.2</td><td>45-49.2</td><td>11.915</td><td>5.067</td><td>0.027</td><td>0.035</td><td>8.904</td><td>1.7452</td><td>4.320</td><td>0.8313</td><td>sf00071_02CT</td></tr><tr><td>4</td><td>32.7</td><td>34-50</td><td>13.65</td><td>4.618</td><td>0.036</td><td>0.048</td><td>4.422</td><td>0.5432</td><td>-0.093</td><td>0.0803</td><td>sf003084_01AG</td></tr><tr><td>12</td><td>78.4</td><td>52-78.4</td><td>14.042</td><td>6.148</td><td>0.024</td><td>0.041</td><td>0.3066</td><td>1.7646</td><td>-1.215</td><td>0.6311</td><td>OMM-5060</td></tr><tr><td>20</td><td>0.0</td><td>0-66</td><td>12.296</td><td>5.012</td><td>0.049</td><td>0.042</td><td>-1.032</td><td>0.8970</td><td>-0.89</td><td>0.2640</td><td>OMM1210</td></tr><tr><td>22</td><td>20.6</td><td>0-20.6</td><td>12.262</td><td>6.302</td><td>0.045</td><td>0.038</td><td>10.499</td><td>4.1333</td><td>3.6166</td><td>1.4460</td><td>sf004055_02CT</td></tr><tr><td>24</td><td>1.5</td><td>0-6</td><td>15.655</td><td>4.703</td><td>0.047</td><td>0.046</td><td>6.2293</td><td>2.3095</td><td>1.7749</td><td>1.2464</td><td>sf003334_01CT</td></tr><tr><td>26</td><td>11.4</td><td>8-14</td><td>12.8</td><td>3.657</td><td>0.045</td><td>0.043</td><td>0.0000</td><td>0.094</td><td>8.8897</td><td>3.7901</td><td>sf003018_09CT</td></tr><tr><td>38</td><td>28.5</td><td>25-29.5</td><td>13.278</td><td>6.273</td><td>0.023</td><td>0.031</td><td>-1.209</td><td>0.8342</td><td>-0.576</td><td>0.5330</td><td>OMM1345</td></tr><tr><td>24</td><td>49.3</td><td>21-49.3</td><td>3.018</td><td>4.359</td><td>0.036</td><td>0.053</td><td>3.1857</td><td>8.3655</td><td>2.0000</td><td>2.1284</td><td>sf003442_01CG</td></tr><tr><td>38</td><td>0.0</td><td>0-11</td><td>3.57</td><td>8.214</td><td>0.027</td><td>0.035</td><td>0.1998</td><td>0.4860</td><td>0.3300</td><td>0.7019</td><td>OMY21INRAiii</td></tr><tr><td>Fulton's condition factor T4</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

(continued)
<table><caption>Table 4, continued</caption><thead><tr><th>Phenotype<sup>a</sup></th><th>LG</th><th>Position, cM</th><th>95% CI, cM</th><th>LOD</th><th>% PVE</th><th>P Value, χ<sup>2</sup></th><th>P Value, F</th><th>Additive Effect</th><th>SE</th><th>Dominance Effect</th><th>SE</th><th>Nearest Marker</th></tr></thead><tbody><tr><td>Liver fresh weight T3, mg</td><td>11</td><td>0.0</td><td>0-11</td><td>6.752</td><td>9.037</td><td>0.006</td><td>0.036</td><td>-0.893</td><td>0.5955</td><td>-0.803</td><td>0.3896</td><td>BX-870052i</td></tr><tr><td></td><td>15</td><td>52.6</td><td>48-52.6</td><td>5.698</td><td>5.24</td><td>0</td><td>0.002</td><td>-1.024</td><td>4.0172</td><td>-0.419</td><td>2.0618</td><td>sf004811_AT</td></tr><tr><td></td><td>20</td><td>58.1</td><td>30-59</td><td>9.649</td><td>5.99</td><td>0</td><td>0.001</td><td>0.9517</td><td>0.3583</td><td>1.2504</td><td>0.3946</td><td>CA-376300i</td></tr><tr><td></td><td>22</td><td>16.3</td><td>12-17</td><td>7.517</td><td>9.788</td><td>0.047</td><td>0.053</td><td>1.6254</td><td>1.2706</td><td>0.0679</td><td>0.8132</td><td>OMM-5147</td></tr><tr><td></td><td>24</td><td>11.9</td><td>8-16</td><td>9.996</td><td>4.01</td><td>0.006</td><td>0.037</td><td>0.4826</td><td>0.0568</td><td>1.3452</td><td>0.0945</td><td>BHMS-238</td></tr><tr><td></td><td>34</td><td>0.0</td><td>0-9</td><td>5.595</td><td>6.965</td><td>0.006</td><td>0.039</td><td>0.3635</td><td>0.3583</td><td>0.4824</td><td>0.1946</td><td>BX319411i</td></tr><tr><td></td><td>40</td><td>52.4</td><td>49-105</td><td>9.808</td><td>7.119</td><td>0.031</td><td>0.051</td><td>-0.595</td><td>1.2706</td><td>0.2837</td><td>0.0132</td><td>SFO-226</td></tr><tr><td>Hematocrit T3, % of red blood cells</td><td>4</td><td>10</td><td>5-45</td><td>3.337</td><td>17.279</td><td>0.002</td><td>0.002</td><td>15.473</td><td>4.4098</td><td>4.3119</td><td>2.2291</td><td>OMM-1220</td></tr><tr><td>Plasma chloride T3, mmol L<sup>-1</sup></td><td>7</td><td>36.6</td><td>19-36.6</td><td>4.428</td><td>5.066</td><td>0.045</td><td>0.051</td><td>0.5712</td><td>0.6385</td><td>0.2201</td><td>0.9368</td><td>OMM-5312i</td></tr><tr><td>Plasma osmolality T3, mmol kg<sup>-1</sup></td><td>10</td><td>42</td><td>35-46</td><td>3.814</td><td>9.257</td><td>0.039</td><td>0.046</td><td>-3.019</td><td>0.6014</td><td>-4.342</td><td>0.4977</td><td>SAL5UoG</td></tr><tr><td>Plasma glucose T3, mg mL<sup>-1</sup></td><td>10</td><td>26</td><td>12-62</td><td>3.241</td><td>11.432</td><td>0.016</td><td>0.02</td><td>14.754</td><td>2.3124</td><td>-1.468</td><td>1.2598</td><td>sf004715_AC</td></tr><tr><td>Hepatic glycogen T3, mg g liver<sup>-1</sup></td><td>2</td><td>40.2</td><td>20-77</td><td>9.48</td><td>4.495</td><td>0.029</td><td>0.04</td><td>9.0189</td><td>4.3652</td><td>4.8747</td><td>2.4390</td><td>sf004818_06CT</td></tr><tr><td></td><td>36</td><td>0.0</td><td>0-21.5</td><td>5.149</td><td>10.576</td><td>0.017</td><td>0.022</td><td>8.9957</td><td>1.7665</td><td>4.3860</td><td>1.7135</td><td>sf004038_AG</td></tr><tr><td>ghr mRNA level T3</td><td>40</td><td>66</td><td>30-95</td><td>4.401</td><td>12.479</td><td>0.032</td><td>0.043</td><td>-16.17</td><td>0.4591</td><td>-13.82</td><td>0.9399</td><td>BX-079862</td></tr><tr><td>igf1 mRNA level T3</td><td>3</td><td>8.1</td><td>5-18</td><td>3.391</td><td>15.612</td><td>0.001</td><td>0.002</td><td>38.289</td><td>1.6790</td><td>14.576</td><td>0.8526</td><td>sf004560_GT</td></tr><tr><td>igf1r mRNA level T3</td><td>3</td><td>10.1</td><td>5-68</td><td>3.178</td><td>13.362</td><td>0.035</td><td>0.047</td><td>41.485</td><td>0.5037</td><td>9.7899</td><td>0.3448</td><td>OMM5008</td></tr><tr><td>Hepato-somatic Index T3</td><td>2</td><td>112.4</td><td>105-117</td><td>8.503</td><td>5.88</td><td>0.019</td><td>0.044</td><td>0.9172</td><td>0.2425</td><td>0.4344</td><td>0.3797</td><td>BX-318599</td></tr><tr><td></td><td>3</td><td>44.8</td><td>12-48</td><td>6.31</td><td>4.443</td><td>0.007</td><td>0.027</td><td>1.2470</td><td>0.1765</td><td>0.8431</td><td>0.2290</td><td>OMM-5007</td></tr><tr><td></td><td>4</td><td>43.0</td><td>32-44</td><td>13.266</td><td>4.454</td><td>0.005</td><td>0.045</td><td>1.1688</td><td>0.2144</td><td>1.0091</td><td>8.2750</td><td>OMM5155i</td></tr><tr><td></td><td>5</td><td>26.7</td><td>23-38</td><td>11.077</td><td>4.866</td><td>0.004</td><td>0.09</td><td>-0.275</td><td>0.1639</td><td>0.6595</td><td>0.1752</td><td>CA-378164</td></tr><tr><td></td><td>7</td><td>13.7</td><td>10-46</td><td>4.624</td><td>4.387</td><td>0</td><td>0.003</td><td>0.9594</td><td>0.2867</td><td>-0.025</td><td>0.5421</td><td>OMM-1263</td></tr><tr><td></td><td>20</td><td>58.1</td><td>20-58.1</td><td>15.991</td><td>4.948</td><td>0.004</td><td>0.049</td><td>0.9517</td><td>0.1440</td><td>1.2504</td><td>0.1586</td><td>CA-376300i</td></tr><tr><td></td><td>22</td><td>17.1</td><td>14-17.1</td><td>13.496</td><td>6.774</td><td>0.002</td><td>0.063</td><td>2.1143</td><td>0.6137</td><td>0.2365</td><td>0.3277</td><td>OMM-3015i</td></tr><tr><td></td><td>23</td><td>13.1</td><td>2-13.1</td><td>10.931</td><td>4.985</td><td>0</td><td>0.03</td><td>-1.252</td><td>0.7850</td><td>-0.535</td><td>0.3950</td><td>BHMS-238</td></tr><tr><td></td><td>38</td><td>28.5</td><td>24-29</td><td>14.168</td><td>4.007</td><td>0.026</td><td>0.021</td><td>0.3635</td><td>0.3590</td><td>0.4824</td><td>0.2296</td><td>OMM1345</td></tr></tbody></table>

Sign (+/-) of the dominance and additive effect indicates to the parent whose allele increases the phenotypic values of the studied trait in the progeny. LG, linkage age; CI, confidence interval; LOD, log<sub>10</sub> of the odd ratio; PVE, percent variance explained; QTL, quantitative trait loci; SE: standard error of the mean; T4, November; T1, May; T3, August; ghr, growth hormone receptor; igf1, insulin growth factor-1; receptor.

a For each phenotype, the sampling period was added (T1, T3, T4) and the unit of measure is given in parentheses.
<table><caption>■ Table 5 Descriptive statistics of LG including the position, 95% CI, LOD score, PVE, associated *P* value, and specific additive, dominance, and interaction effects of each QTL linked to every phenotype related to stress response</caption><thead><tr><th>Phenotype<sup>a</sup></th><th>LG</th><th>Position, cM</th><th>95% CI, cM</th><th>LOD</th><th>% PVE</th><th>P Value, χ<sup>2</sup></th><th>P Value, F</th><th>Additive Effect</th><th>SE</th><th>Dominance Effect</th><th>SE</th><th>Nearest Marker</th></tr></thead><tbody><tr><td>Plasma cortisol, µg dL<sup>-1</sup></td><td>14</td><td>68.5</td><td>27-68.5</td><td>3.217</td><td>3.854</td><td>0.013</td><td>0.015</td><td>0.0056</td><td>0.0091</td><td>8.8723</td><td>2.5333</td><td>sf003382_AG</td></tr><tr><td>Plasma osmolality, mmol kg<sup>-1</sup></td><td>23</td><td>0.1</td><td>0-16</td><td>8.161</td><td>31.35</td><td>0</td><td>0</td><td>0.0000</td><td>0.0098</td><td>0.1119</td><td>0.0168</td><td>sf004319_GT</td></tr><tr><td>Plasma chloride, mmol L<sup>-1</sup></td><td>10</td><td>42</td><td>37-47</td><td>3.667</td><td>11.232</td><td>0.039</td><td>0.046</td><td>-2.123</td><td>0.5945</td><td>-3.9858</td><td>0.5363</td><td>SAL5UoG</td></tr><tr><td></td><td>7</td><td>36.6</td><td>17-38</td><td>3.532</td><td>8.734</td><td>0.045</td><td>0.049</td><td>0.6324</td><td>0.4987</td><td>0.2201</td><td>0.9856</td><td>OMM-5312i</td></tr></tbody></table>

LG, linkage age; CI, confidence interval; LOD, log₁₀ of the odd ratio; PVE, percent variance explained; QTL, quantitative trait loci; SE: standard error of the mean.

with the actual “causative locus.” For instance, previous results obtained for another salmonid, the lake whitefish *Coregonus clupea-formis*, showed that expression QTL were located within a hotspot of recombination where a putative *trans*-regulating element was identified (Zinc finger protein 35; reviewed in Bernatchez *et al.* 2010).

Of the 63 growth-related QTL in this study, similar proportions were associated with dominance (52.4%; 33/63) or additive (47.6%; 30/63) effect. This measure gives insight about the genetic architecture of the trait, i.e., about the gene action associated with the QTL and the effect of the underlying alleles. Sign (−/+), of the dominance or additive effect, at each locus, corresponds to the parent whose allele increases the phenotypic values of the studied trait in the progeny (Lynch and Walsh 1998).

## Correlations between phenotypes related to growth performance

A positive correlation was observed between glucose plasma level and the expression of genes related to growth performance (*ghr*, *igf1*, *igf1r*). It is known that glucose availability and utilization are regulated by endocrine factors involved in appetite, growth, and energy homeostasis, and that glucose exerts a feedback loop on these endocrine pathways. In particular, glucose plays a role in the GH/IGF1 growth axis in all vertebrates,

including fish (Pérez-Sánchez *et al.* 1994; Banos *et al.* 1998; Riley *et al.* 2009). In particular, our results are in accordance with those of Riley *et al.* (2009), who demonstrated that a dietary glucose load significantly increased the mRNA level of two genes coding for two growth hormone receptors (*ghr1* and *ghr2*) in tilapia *Oreochromis mossambicus*. Moreover, Banos *et al.* (1998) showed that rainbow trout fed high-carbohydrate diets had an increased liver *Igf1r* number and specific binding. Contrary to our results, however, Riley *et al.* (2009) reported a stable liver *igf1* mRNA level as well as a reduced plasma *igf1* level in tilapia with a plasma glucose load. This finding could indicate that the relationship between the plasma glucose level and the GH/IGF1 axis may differ according to the species considered. We also measured a positive correlation between glucose plasma level and mRNA levels of *ef1* and *β-actin*. *β-actin* is involved in cell motility, structure, and integrity, whereas *ef1* has been shown to play a central role in protein synthesis in teleosts (Jurss *et al.* 1992). Our results thus suggest the essential role of glucose as energy supplier to the cell to perform these vital functions in *S. fontinalis*.

## QTL related to stress response

The present study identified four QTL associated with the three phenotypes related to stress response and localized on four different

<table><caption>■ Table 6 Detailed information related to the SNP markers linked to QTL associated with growth and stress response</caption><thead><tr><th>Phenotype</th><th>SNP Name</th><th>Accession Number</th><th>GI</th><th>Variation</th><th>Ts/Tv</th><th>C/NC</th><th>S/NS</th><th>LG</th><th>Available Annotation</th></tr></thead><tbody><tr><td>Growth</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Size</td><td>sf000071</td><td></td><td></td><td>A/C</td><td>Tv</td><td>NC</td><td>—</td><td>3</td><td>—</td></tr><tr><td>Size</td><td>sf003018</td><td></td><td></td><td>A/G</td><td>Ts</td><td>NC</td><td>—</td><td>26</td><td>—</td></tr><tr><td>Size</td><td>sf003084</td><td></td><td></td><td>A/G</td><td>Ts</td><td>C</td><td>—</td><td>4</td><td>—</td></tr><tr><td>Size</td><td>sf003334</td><td></td><td></td><td>A/G</td><td>Ts</td><td>NC</td><td>S</td><td>24</td><td>—</td></tr><tr><td>Fulton index</td><td>sf003442</td><td></td><td></td><td>A/G</td><td>Ts</td><td>C</td><td>S</td><td>24</td><td>—</td></tr><tr><td>Hepatic glycogen</td><td>sf004038</td><td>NM_001165151</td><td>GI:259089170</td><td>A/G</td><td>Ts</td><td>NC</td><td>—</td><td>34</td><td>Ubiquitin-conjugating enzyme E2W</td></tr><tr><td>Size</td><td>sf004055</td><td>XP_001921123</td><td>GI:189530039</td><td>A/T</td><td>Tv</td><td>C</td><td>NS</td><td>22</td><td>Neural-cadherin-like</td></tr><tr><td>igf1 mRNA level</td><td>sf004560</td><td></td><td></td><td>G/T</td><td>Tv</td><td>NC</td><td>—</td><td>3</td><td>—</td></tr><tr><td>Plasma glucose</td><td>sf004715</td><td></td><td></td><td>C/T</td><td>Ts</td><td>C</td><td>S</td><td>10</td><td>—</td></tr><tr><td>SGR</td><td>sf004732</td><td>AB258536</td><td>GI:118596560</td><td>C/T</td><td>Ts</td><td>NC</td><td>—</td><td>29</td><td>Onmy-LDA gene for MHC class I antigen</td></tr><tr><td>Liver fresh weight</td><td>sf004811</td><td>EU025709</td><td>GI:158702304</td><td>A/T</td><td>Tv</td><td>C</td><td>NS</td><td>15</td><td>—</td></tr><tr><td>Hepatic glycogen</td><td>sf004818</td><td>EU481821</td><td>GI:171474994</td><td>A/C</td><td>Tv</td><td>NC</td><td>—</td><td>2</td><td>Formin-binding protein 1</td></tr><tr><td>Stress response</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Plasma cortisol</td><td>sf003382</td><td></td><td></td><td>A/G</td><td>Ts</td><td>C</td><td>NS</td><td>14</td><td>—</td></tr><tr><td>Plasma cortisol</td><td>sf004319</td><td>BT058994</td><td>GI:223647897</td><td>G/T</td><td>Tv</td><td>C</td><td>S</td><td>23</td><td>Hydroxymethylglutaryl-CoA lyase, mitochondrial precursor putative mRNA</td></tr></tbody></table>

SNP, single nucleotide polymorphism; QTL, quantitative trait loci; GI, GenInfo identifiers; Ts, transition; Tv, transversion; C, coding region; NC, noncoding region; S, synonymous; NS, nonsynonymous; LG, linkage group.

<footnote>ª All the phenotypes related to the stress response were measured at T4 (November) and the unit of measure is given in parentheses.</footnote>
<figure><img src="image_12.png" /><figcaption>**Figure 2** Distribution of QTL related to growth and stress response across the 40 LG identifying the chromosome-wide (gray blocks) and genome-wide significance levels (black blocks). SGR, specific growth rate; HSI, hepatosomatic index; T1, May; T2, July; T3, August; T4, November; *ghr*, *growth hormone receptor*; *igf1*, *insulin growth factor-1*; *igf1r*, *insulin growth factor-1 receptor*.</figcaption></figure>

LGs. Two QTL were identified for the level of plasma cortisol (LG 14 and 23; chromosome- and genome-wide level of significance, respectively), whereas only one QTL was identified for plasma osmolality (LG 10; chromosome-wide level of significance) and one for plasma chloride (LG 7; chromosome-wide level of significance). The PVE estimation of 31.35% revealed a major effect for the QTL associated with plasma cortisol (LG 23). The PVE was weaker for another QTL associated with plasma cortisol (LG 14; PVE = 3.85%), as well as QTL associated with plasma osmolality (LG 10; PVE = 11.23%), and plasma chloride (LG 7; PVE = 8.73%). The small number of QTL identified for these traits combined with general high PVE values suggests that the genetic architecture of these four phenotypes is under the control of a limited number of genes or a major gene effect. Low number of QTL related to the cortisol stress response has previously been reported in other vertebrates (pigs: Ousova *et al.* 2004; rodents: Garlow *et al.* 2005 and Llamas *et al.* 2005; and poultry: Beaumont *et al.* 2005). Moreover, Vallejo *et al.* (2009) and Quillet *et al.* (2010) reported a limited number of major genes affecting plasma cortisol levels after a crowding stress in rainbow trout. Again, because QTL analyses are biased toward the identification of large QTL effects, the power of our analysis (n = 171 individuals) may be insufficient to detect low to moderate QTL effects influencing these traits. However, this cannot in itself explain the contrast observed between these QTL and those for phenotypic traits related to growth for which a larger number of small-effect QTL was revealed.

In the present study, both QTL associated with plasma cortisol were colocalized with two growth-related QTL (SGR located on LG 14 and 23), suggesting a genetic correlation between stress response and growth performance. This finding is in accordance with Drew *et al.* (2007), who detected one QTL associated with juvenile rainbow trout body mass overlapping with a QTL for cortisol levels, indicating a putative genetic link between these two traits. Moreover, these authors reported additive effects of these two overlapping QTL, suggesting

a positive relationship between cortisol levels and growth rate. However, this finding is in contrast with Fevolden *et al.* (2002), who previously reported a negative correlation between stress response and growth performance in rainbow trout. However, this statement should be interpreted cautiously when comparing both studies, as fish were younger in Fevolden *et al.* (2002) and were not subjected to any domestication process that may alter stress response and growth. Taken together, our results and previous ones suggest that the cortisol response after a stress exposure may be under variable genetic control depending on species. This also demonstrates that stress response measured through plasma cortisol is a highly complex trait and may influence a range of underlying physiological mechanisms. In the case of brook charr however, it seems clear that stress responses are under strong genetic control given recently published high heritability estimates ($h^2$ = 0.60 ± 0.20; $h^2$ = 0.61 ± 0.20) for plasma cortisol and plasma glucose, respectively (Crespel *et al.* 2011). The identification of QTL related to the same phenotypic traits supports the strong potential role of these QTL for genetic improvement of the stress response in brook charr and possibly other salmonid species.

## Identification of SNP markers linked to identified QTL related to stress response

For each QTL position, the closest underlying molecular marker has been identified. Within the four stress-related QTL, two were linked to an SNP marker (sf003382 and sf004319 linked to plasma cortisol on LG 14 and LG 23, respectively), whereas the two others were linked to SSR markers (SAL5UoG and OMM-5312i linked to plasma osmolality and plasma chloride on LG 10 and LG 7, respectively). For the two QTL related to plasma cortisol, the closest linked SNP markers were both identified in a coding region but were characterized as a nonsynonymous transition in one case (sf003382) and as a silent transversion in the other (sf004319). SNP sf004319 was significantly
annotated as a precursor of the mitochondrial HMG-CoA (accession number [GenBank:BT058994]). This gene catalyses the transformation of HMG-CoA into acetyl-CoA and acetoacetate. In vertebrates, HMG-CoA is a mitochondrial and paroxysmal enzyme that is involved in ketogenesis and in leucine catabolism (Mitchell *et al.* 1993). As was the case for growth-related QTL that were linked to SNP markers, direct assumptions of the role played by these two SNPs linked to stress-related QTL remain difficult, especially since the molecular function of the only annotated marker was not directly related to stress response. Here again, we suggest that these SNP markers are in linkage disequilibrium with the causative locus (or transcription factors). As such, our results should be perceived as representing a necessary first step toward the identification of the genes underlying the genetic architecture of stress response in fish and of growth performance in *S. fontinalis*. The identification of actual causative locus will require the availability of intensive genetic information and genomic tools, such as dense, gene-rich linkage maps and whole-genome sequences, that are currently in development in salmonids (Davidson *et al.* 2010).

## Conclusion

For the first time in brook charr, and to our knowledge only the second time in the genus *Salvelinus*, this study examined the association between molecular markers and the variance of numerous traits related to growth performance and stress response. Our results can be divided into two groups. The first group may be composed of growth performance-related QTL that are numerous per trait (>2; size, weight, SGR, HSI, liver fresh weight, Fulton index) and display small-to-moderate effects. This suggests that these growth-related traits, as expected for such a complex trait, are under the control of a large number of genes (pleiotropic genes). The second group may be composed of QTL that are less numerous per trait (< 2) and display larger effects, suggesting that these QTL are under the control of a limited number of genes of more important effect. These less numerous QTL were mainly identified for some growth-related traits (hematocrit, plasma chloride, plasma glucose, plasma osmolality, hepatic glycogen and mRNA levels of *igf1*, *igf1r*, *gh*, *ef1*, and *β-actine*) and for all the traits related to stress response (plasma chloride, plasma cortisol, and plasma osmolality). Colocalization of QTL associated with plasma cortisol and with SGR has been identified for the first time in *S. fontinalis*. Admittedly, the identification of genes underlying the genetic architecture of traits of interest, such as growth or stress response, will require improvement in a subsequent map generation. This will require an ever-increasing amount of genetic information and research tools to increase the density of linkage maps. Linkage map improvement will also benefit from the increasing availability of whole-genome sequences in fish that are currently in development in salmonids (Davidson *et al.* 2010) and the more widespread use of the latest NGS techniques for development in nonmodel species in general (Amores *et al.* 2011; Rowe *et al.* 2011).

## ACKNOWLEDGMENTS

We thank S. Granier, R. Gagné, H. Segouani, and G. Côté for their invaluable help during the phenotyping experiments and sample collection; A. Lemieux for her help with qPCR experiments; E. Parent for his help in sequencing; E. Normandeau for his help in bioinformatics; L. Papillon and N. Brodeur for their help with microsatellite genotyping; and A. Montpetit and A. Belisle from Genome Québec Innovation Centre (McGill University, Montréal, QC, Canada) for their help in cDNA library sequencing and SNP genotyping. We also thank Associate editor Andrew Paterson and an anonymous referee for their useful comments and suggestions. This work was founded by a Part-

nership grant (Strategic program) of the Natural Sciences and Engineering Research Council of Canada (NSERC) to L.B., N.D., and C.A., and by an FQRNT (Fonds Québécois pour la Recherche, Nature et Technologies) post-doctoral fellowship for M.V. This is a contribution to the research program of Réseau Aquaculture Québec (RAQ). The authors have no conflict of interest to declare and note that the sponsors of the issue had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

## LITERATURE CITED

Adams, P., 1990 *Saltmarsh Ecology*, Cambridge University Press, Cambridge, UK.

Ali, M., G. Nicieza, and R. Wootton, 2003 Compensatory growth in fishes: a response to growth depression. Fish Fish. 4: 147-190.

Amores, A., J. Catchen, A. Ferrara, Q. Fontenot, and J. H. Postlethwait, 2011 Genome evolution and meiotic maps by massively parallel DNA sequencing: spotted gar, an outgroup for the teleost genome duplication. Genetics 188: 799-808.

Audet, C., and G. Claireaux, 1992 Diel and seasonal changes in resting levels of various blood parameters in brook trout (*Salvelinus fontinalis*). Can. J. Fish. Aquat. Sci. 49: 870-877.

Banos, N., J. Baro, C. Castejon, I. Navarro, and J. Gutierrez, 1998 Influence of high-carbohydrate enriched diets on plasma insulin levels and insulin and IGF-I receptors in trout. Regul. Pept. 77: 55-62.

Barton, B. A., 2002 Stress in fishes: a diversity of responses with particular reference to changes in circulating corticosteroids. Integr. Comp. Biol. 42: 517-525.

Barton, B. A., and G. K. Iwama, 1991 Physiological changes in fish from stress in aquaculture with emphasis on the response and effects of corticosteroids. Annu. Rev. Fish Dis. 1: 26.

Barton, B. A., and C. B. Schreck, 1987 Influence of acclimation temperature on interrenal and carbohydrate stress responses in juvenile chinook salmon (*Oncorhynchus tshawytscha*). Aquaculture 62: 299-310.

Bastien, A., 2010 Évaluation d'un programme de sélection et identification des traits physiologiques liés à l'anadromie chez l'omble de fontaine (*Salvelinus fontinalis*). Ph.D. Thesis, Université du Québec à Rimouski, Rimouski (Québec) Canada.

Beaumont, C., O. Roussot, K. Feve, F. Vignoles, S. Leroux *et al.*, 2005 A genome scan with AFLP markers to detect fearfulness-related QTL in Japanese quail. Anim. Genet. 36: 401-407.

Beavis, W. D., 1998 QTL analyses: Power, precision and accuracy, pp. 145-161 in *Molecular Dissection of Complex Traits*, edited by A. H. Paterson. CRC Press, Boca Raton, FL.

Bennewitz, J., and T. H. Meuwissen, 2010 The distribution of QTL additive and dominance effects in porcine F2 crosses. J. Anim. Breed. Genet. 127: 171-179.

Bernatchez, L., S. Renaut, A. R. Whiteley, N. Derome, J. Jeukens *et al.*, 2010 On the origin of species: insights from the ecological genomics of lake whitefish. Philos. Trans. R. Soc. Lond. B Biol. Sci. 365: 1783-1800.

Blanchet, S., G. Loot, G. Grenouillet, and S. Brosse, 2007 Competitive interactions between native and exotic salmonids: a combined field and laboratory demonstration. Ecol. Freshwat. Fish 16: 133-143.

Bougas, B., S. Granier, C. Audet, and L. Bernatchez, 2010 The transcriptional landscape of cross-specific hybrids and its possible link with growth in brook charr (*Salvelinus fontinalis* Mitchell). Genetics 186: 97-107.

Brockman, W., P. Alvarez, S. Young, M. Garber, G. Giannoukos *et al.*, 2008 Quality scores and SNP detection in sequencing-by-synthesis systems. Genet. Res. 18: 763-770.

Broman, K. W., H. Wu, S. Sen, and G. A. Churchill, 2003 *R/qlt*: QTL mapping in experimental crosses. Bioinformatics 19: 889-890.

Bureau, D., K. Hua, and C. Cho, 2006 Effect of feeding level on growth and nutrient deposition in rainbow trout (*Oncorhynchus mykiss* Walbaum) growing from 150 to 600 g. Aquacult. Res. 37: 1090-1098.

Carr, R. S., and J. M. Neff, 1984 Quantitative semi-automated enzymatic assay for tissue glycogen. Comp. Biochem. Physiol. 77: 447-449.

Castric, V., and L. Bernatchez, 2003 The rise and fall of isolation by distance in the anadromous brook charr (*Salvelinus fontinalis*, Mitchell). Genetics 163: 983-996.
Côté, G., G. Perry, P. Blier, and L. Bernatchez, 2007 The influence of gene-environment interactions on GHR and IGF-1 expression and their association with growth in brook charr, *Salvelinus fontinalis* (Mitchill). BMC Genet. 8: 87.

Crespel, A., L. Bernatchez, D. Garant, and C. Audet, 2011 Quantitative genetic analysis of the physiological stress response in three strains of brook charr, *Salvelinus fontinalis* (Mitchill), and their hybrids. J. Fish Biol. (in press).

Davidson, W. S., B. F. Koop, S. J. Jones, P. Iturra, R. Vidal et al., 2010 Sequencing the genome of the Atlantic salmon (*Salmo salar*). Genome Biol. 11: 403.

Devlin, R., and Y. Nagahama, 2002 Sex determination and sex differentiation in fish: an overview of genetic, physiological, and environmental influences. Aquaculture 208: 191-364.

Drew, R., H. Schwabl, P. Wheeler, and G. Thorgaard, 2007 Detection of QTL influencing cortisol levels in rainbow trout (*Oncorhynchus mykiss*). Aquaculture 272: 183-194.

Du, S., R. Devlin, and C. Hew, 1993 Genomic structure of growth hormone genes in Chinook salmon (*Oncorhynchus tshawyscha*): presence of two functional genes, GH-I and GH-II, and a male specific pseudogene, GH-psi. DNA Cell Biol. 12: 739-751.

Fevolden, S. E., K. H. Roed, and K. T. Fjalestad, 2002 Selection response of cortisol and lysozyme in rainbow trout and correlation to growth. Aquaculture 205: 61-75.

Friendly, M., 2002 Corrgrams: exploratory displays for correlation matrices. Am. Stat. 56: 316-324.

Garlow, S. J., E. Boone, W. Li, M. J. Owens, and C. B. Nemeroff, 2005 Genetic analysis of the hypothalamic corticotropin-releasing factor system. Endocrinology 146: 2362-2368.

Haley, C., and D. J. D. Koning, 2006 Genetical genomics in livestock: potentials and pitfalls. Anim. Genet. 37: 10-12.

Haley, C. S., and S. A. Knott, 1992 A simple regression method for mapping quantitative trait loci in line crosses using flanking markers. Heredity 69: 315-324.

Houston, R. D., C. S. Haley, A. Hamilton, D. R. Guy, A. E. Tinch et al., 2008 Major quantitative trait loci affect resistance to infectious pancreatic necrosis in Atlantic salmon (*Salmo salar*). Genetics 178: 1109-1115.

Iversen, M., B. Finstad, and K. Nilssen, 1998 Recovery from loading and transport stress in Atlantic salmon (*Salmo salar* L.) smolts. Aquaculture 168: 387-394.

Jackson, T. R., M. M. Ferguson, R. G. Danzmann, A. G. Fishback, P. E. Ihssen et al., 1998 Identification of two QTL influencing upper temperature tolerance in three rainbow trout (*Oncorhynchus mykiss*) half-sib families. Heredity 80: 143-151.

Jansen, R. C., J. W. Ooijen, P. Stam, C. Lister, and C. Dean, 1995 Genotype-by-environment interaction in genetic mapping of multiple quantitative trait loci. TAG 91: 33-37.

Jurss, K., I. Junghahn, and R. Bastrop, 1992 The role of elongation factors in protein synthesis rate variation in white teleost muscle. J. Comp. Physiol. B 162: 345-350.

Kuttner, E., H. K. Moghadam, S. Skulason, R. G. Danzmann, and M. M. Ferguson, 2011 Genetic architecture of body weight, condition factor and age of sexual maturation in Icelandic Arctic charr (*Salvelinus alpinus*). Mol. Genet. Genomics 286: 67-79.

Lander, E., and L. Kruglyak, 1995 Genetic dissection of complex traits: guidelines for interpreting and reporting linkage results. Nat. Genet. 11: 241-247.

Llamas, B., V. Contesse, V. Guyonnet-Duperat, H. Vaudry, P. Mormede et al., 2005 QTL mapping for traits associated with stress neuroendocrine reactivity in rats. Mamm. Genome 16: 505-515.

Lynch, M., and B. Walsh, 1998 Genetics and Analysis of Quantitative Traits. Sinauer Associates, Sunderland, MA.

Mackay, T. F., 2001 The genetic architecture of quantitative traits. Annu. Rev. Genet. 35: 303-339.

Mackay, T. F., E. A. Stone, and J. F. Ayroles, 2009 The genetics of quantitative traits: challenges and prospects. Nat. Rev. Genet. 10: 565-577.

Makinen, T., and K. Ruohonen, 1992 Effect of delayed photoperiod on the growth of a Finnish rainbow trout (*Oncorhynchus mykiss* Walbaum) stock. J. Appl. Ichthyology 8: 40-50.

Martin, S., J.-Y. Savaria, C. Audet, and L. Bernatchez, 1997 Microsatellites reveal no evidence for inbreeding effects but low-inter-stock genetic diversity among brook charr stocks used for production in Québec. Bull. Aquacult. Assoc. Can. 97: 21-23.

Martyniuk, C., G. Perry, H. Moghadam, M. Ferguson, and R. Danzmann, 2003 The genetic architecture of correlations among growth-related traits and male age at maturation in rainbow trout. J. Fish Biol. 63: 746-764.

Massault, C., B. Hellemans, B. Louro, C. Batargias, J. K. Van Houdt et al., 2010 QTL for body weight, morphometric traits and stress response in European sea bass *Dicentrarchus labrax*. Anim. Genet. 41: 337-345.

McClelland, E. K., and K. Naish, 2010 Quantitative trait locus analysis of hatch timing, weight, length and growth rate in coho salmon. Heredity 105: 562-573.

Metcalfe, N., 1986 Intraspecific variation in competitive ability and food intake in salmonids: consequences for energy budgets and growth rates. J. Fish Biol. 28: 525-531.

Mitchell, G. A., M. F. Robert, P. W. Hruz, S. Wang, G. Fontaine et al., 1993 3-Hydroxy-3-methylglutaryl coenzyme A lyase (HL). Cloning of human and chicken liver HL cDNAs and characterization of a mutation causing human HL deficiency. J. Biol. Chem. 268: 4376-4381.

Moghadam, H., M. Ferguson, C. Rexroad, I. Coulibaly, and R. Danzmann, 2007a Genomic organization of the IGF1, IGF2, MYF5, MYF6, and GRF/PACAP genes across *Salmoninae* genera. Anim. Genet. 38: 527-532.

Moghadam, H., J. Poissant, H. Fotherby, L. Haidle, M. Ferguson et al., 2007b Quantitative trait loci for body weight, condition factor and age at sexual maturation in Arctic charr (*Salvelinus alpinus*): comparative analysis with rainbow trout (*Oncorhynchus mykiss*) and Atlantic salmon (*Salmo salar*). Mol. Genet. Genomics 277: 647-661.

O'Malley, K., T. Sakamoto, R. Danzmann, and M. Ferguson, 2003 Quantitative trait loci for spawning date and body weight in rainbow trout: testing for conserved effects across ancestrally duplicated chromosomes. J. Hered. 94: 273-284.

Ousova, O., V. Guyonnet-Duperat, N. Iannuccelli, J. P. Bidanel, D. Milan et al., 2004 Corticosteroid binding globulin: a new target for cortisol-driven obesity. Mol. Endocrinol. 18: 1687-1696.

Ozaki, A., T. Sakamoto, S. Khoo, K. Nakamura, M. Coimbra et al., 2001 Quantitative trait loci (QTL) associated with resistance/susceptibility to infectious pancreatic necrosis virus (IPNV) in rainbow trout (*Oncorhynchus mykiss*). Mol. Genet. Genomics 265: 23-31.

Perry, G. M. L., M. M. Ferguson, T. Sakamoto, and R. G. Danzmann, 2005 Sex-linked quantitative trait loci for thermotolerance and length in the rainbow trout. J. Hered. 96: 97-107.

Pfaffl, M. W., 2001 A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 29: e45.

Pickering, A. D., 1981 Stress and Fish, Academic Press, New York.

Portz, D. E., C. M. Woodley, and J. J. Cech, 2006 Stress-associated impacts of short-term holding on fishes. Rev. Fish Biol. Fish. 16: 125-170.

Pottinger, T. G., and A. D. Pickering, 1997 Genetic basis to the stress response: selective breeding for stress-tolerant fish, pp. 171-193 in *Fish stress and health in aquaculture*, edited by Iwama, G. K., A. D. Pickering, J. P. Sumpter, and C. B. Schreck. Cambridge University Press, Cambridge, UK.

Quillet, E., F. Krieg, N. Dechamp, C. Hervet, P. Le Roy et al., 2010 Mapping QTL affecting cortisol response to confinement stress in rainbow trout. 9th World Congress on Genetics Applied to Livestock Production, Leipzig, Germany, August 1-6, 2010.

Reid, D., A. Szanto, B. Glebe, R. Danzmann, and M. Ferguson, 2005 QTL for body weight and condition factor in Atlantic salmon (*Salmo salar*): comparative analysis with rainbow trout (*Oncorhynchus mykiss*) and Arctic charr (*Salvelinus alpinus*). Heredity 94: 166-172.

Remington, D. L., J. M. Thornsberry, Y. Matsuoka, L. M. Wilson, S. R. Whitt et al., 2001 Structure of linkage disequilibrium and phenotypic associations in the maize genome. Proc. Natl. Acad. Sci. USA 98: 11479-11484.

Riley, L. G. Jr., A. P. Walker, C. P. Dborough, S. E. Schwandt, and E. G. Grau, 2009 Glucose regulates ghrelin, neuropeptide Y, and the GH/IGF-I axis in the tilapia, *Oreochromis mossambicus*. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 154: 541-546.

Rockman, M. V., 2011 The QTN program and alleles that matter for evolution: all that's gold does not glitter. Evolution 66: 1-17.
Rowe, H. C., S. Renaut, and A. Guggisberg, 2011 RAD in the realm of next-generation sequencing technologies. Mol. Ecol. 20: 3499–3502.

Rozen, S., and H. Skaletsky, 2000 Primer3 on the WWW for general users and for biologist programmers, pp. 365–386 in *Bioinformatics Methods and Protocols: Methods in Molecular Biology*, edited by S. Krawetz, and S. Misener. Humana Press, Totowa, NJ.

Sauvage, C., P. Boudry, D. J. de Koning, C. S. Haley, S. Heurtebise et al., 2010a QTL for resistance to summer mortality and OsHV-1 load in the Pacific oyster (*Crassostrea gigas*). Anim. Genet. 41: 390–399.

Sauvage, C., M. Vagner, N. Derome, C. Audet, and L. Bernatchez 2012 Genetic mapping of SNP markers and its application in QTL detection for reproductive traits in brook charr, *Salvelinus fontinalis*. G3 (Bethesda) 2: 379–392.

Schaffer, W., 1979 Equivalence of maximizing reproductive value and fitness in the case of reproductive strategies. Proc. Natl. Acad. Sci. USA 76: 3567–3569.

Slate, J. O. N., 2005 Quantitative trait locus mapping in natural populations: progress, caveats and future directions. Mol. Ecol. 14: 363–379.

Thorpe, J., 1994 Performance thresholds and life-history flexibility in salmonids. Cons Bio 8: 877–879.

Vallejo, R. L., C. E. Rexroad III, J. T. Silverstein, L. L. G. Janss, and G. M. Weber, 2009 Evidence of major genes affecting stress response in

rainbow trout using Bayesian methods of complex segregation analysis. J. Anim. Sci. 87: 3490–3505.

Van Ooijen, J. W., 1999 LOD significance thresholds for QTL analysis in experimental populations of diploid species. Heredity 83: 613–624.

Vijayan, M. M., and T. W. Moon, 1992 Acute handling stress alters hepatic glycogen metabolism in food-deprived rainbow trout (*Oncorhynchus mykiss*). Can. J. Fish. Aquat. Sci. 49: 2260–2266.

Visscher, P. M., C. Haley, and S. A. Knott, 1996 Mapping QTL for binary traits in backcross and F2 populations. Genet. Res. 68: 55–63.

Wang, Y., D. Zhu, X. Wu, Y. W. Terry Healy, and H. Judy-Ann, 2002 Tidal flats and associated muddy coast of China, pp. 319–345 in *Proceedings in Marine Science*. Elsevier, Amsterdam.

Wendelaar Bonga, S. E., 1997 The stress response in fish. Physiol. Rev. 77: 591–625.

Wringe, B. F., R. H. Devlin, M. M. Ferguson, H. K. Moghadam, D. Sakhrani et al., 2010 Growth-related quantitative trait loci in domestic and wild rainbow trout (*Oncorhynchus mykiss*). BMC Genet. 11: 63.

Xu, S., 2003 Theoretical basis of the Beavis effect. Genetics 165: 2259–2268.

Communicating editor: A. H. Paterson