QTL affecting stress response to crowding in a rainbow trout broodstock population

Background Genomic analyses have the potential to impact selective breeding programs by identifying markers that serve as proxies for traits which are expensive or difficult to measure. Also, identifying genes affecting traits of interest enhances our understanding of their underlying biochemical pathways. To this end we conducted genome scans of seven rainbow trout families from a single broodstock population to identify quantitative trait loci (QTL) having an effect on stress response to crowding as measured by plasma cortisol concentration. Our goal was to estimate the number of major genes having large effects on this trait in our broodstock population through the identification of QTL. Results A genome scan including 380 microsatellite markers representing 29 chromosomes resulted in the de novo construction of genetic maps which were in good agreement with the NCCCWA genetic map. Unique sets of QTL were detected for two traits which were defined after observing a low correlation between repeated measurements of plasma cortisol concentration in response to stress. A highly significant QTL was detected in three independent analyses on Omy16, many additional suggestive and significant QTL were also identified. With linkage-based methods of QTL analysis such as half-sib regression interval mapping and a variance component method, we determined that the significant and suggestive QTL explain about 40-43% and 13-27% of the phenotypic trait variation, respectively. Conclusions The cortisol response to crowding stress is a complex trait controlled in a sub-sample of our broodstock population by multiple QTL on at least 8 chromosomes. These QTL are largely different from others previously identified for a similar trait, documenting that population specific genetic variants independently affect cortisol response in ways that may result in different impacts on growth. Also, mapping QTL for multiple traits associated with stress response detected trait specific QTL which indicate the significance of the first plasma cortisol measurement in defining the trait. Fine mapping these QTL can lead towards the identification of genes affecting stress response and may influence approaches to selection for this economically important stress response trait.


Background
Employing molecular approaches to understanding complex genetic traits which are difficult or expensive to measure has the potential to facilitate genetic improvement through selective breeding and the identification of population specific genetic variants which affect important phenotypes [1]. Although some populations may exhibit similar phenotypes, the genetic architectures underlying a trait may be vastly different. Such traits include those associated with the physiological responses of rainbow trout to the stressors of aquaculture production environments. Typical stressors can be categorized under handling and manipulation, overcrowding, suboptimal water quality parameters, and social interactions, all of which have been shown to negatively affect production traits such as growth, feed intake, feed efficiency, disease resistance, flesh quality, and reproductive performance [2][3][4][5][6][7]. Fish respond to these stressors in similar patterns to those exhibited by terrestrial vertebrates, by stimulating the hypothalamic-sympatheticchromaffin cell axis which releases catecholamines to increase oxygen uptake and energy mobilization and the hypothalamic-pituitary-interrenal axis, which produces cortisol to effect carbohydrate, protein, and lipid metabolism [8]. As in mammals, an appropriate hormonal response enables the fish to endure a stress; however, excessive or prolonged exposure to these hormones, particularly cortisol, can also have deleterious effects on traits associated with aquaculture production [3].
Despite clear evidence that these hormonal responses are central to the overall stress response, attempts to use measures of the hormonal responses as predictors of performance in rainbow trout, have been equivocal at best. A good example is the response to crowding stress in rainbow trout. In an European population, response to crowding stress was observed to be heritable, with fish responding with elevated plasma cortisol concentrations having inferior growth when reared in co-culture with low responders [9,10]. Similarly studies on fish in the US have also shown that this trait is heritable; however, in contrast, fish responding to stress with elevated concentrations of plasma cortisol were observed to be superior growers [11,12]. These differences in association between the physiological marker for stress response and the production trait indicates genetic variation in multiple biochemical mechanisms controlling growth and/or stress responses in these populations. Therefore our goal was to survey multiple families from a single broodstock population to characterize the diversity of loci potentially affecting stress response to crowding as a first step in using allelic variation for parsing out the crude phenotype of post-stressor plasma cortisol levels into more refined phenotypes that will likely better associate with production traits as predicted.
In 2002 the National Center for Cool and Cold Water Aquaculture (NCCCWA) initiated a selective breeding program with genetically diverse strains of rainbow trout with the aim of creating stocks that are improved for aquaculture production efficiency [13]. To date selection has focused on characteristics associated with growth and disease resistance [14][15][16]. Stress has been shown to affect both of these traits [8]. For initial evaluation seven multi-generation families from the growth line were selected for identifying quantitative trait loci (QTL) for stress response based on their plasma cortisol response to crowding [12]. Previously, complex segregation analysis (CSA) has suggested that one or more major genes is affecting stress response to crowding in this population, and that a dominant gene has a large negative effect on plasma cortisol concentration in response to crowding stress [17]. Therefore our goals in this study were to: 1) identify the number of loci affecting stress response in our broodstock population; and 2) validate the findings of the previous CSA. To this end we evaluated seven families containing a total of 222 offspring that were measured for their stress response according to the protocol of Pottinger and Carrick [18] as modified by Weber and Silverstein [12], and conducted a genome scan with 380 microsatellite markers selected from the NCCCWA genetic map [19]. We defined two traits for QTL mapping. One phenotype was based on estimated breeding values (EBV) using four plasma cortisol repeated measurements per animal. The second phenotype was based on best linear unbiased predictors of the last three repeated measurements (BLUP3) as we observed that the initial cortisol measurement for each fish did not correlate well with the last three measurements.

Methods
The stress challenge and collection of blood for measuring plasma cortisol concentration was conducted with the approval of the IUCAC of the USDA/ARS National Center for Cool and Cold Water Aquaculture, protocol #50.

Mapping population
The QTL mapping population was identified from a broodstock population at the NCCCWA in Leetown, West Virginia, USA. This selective breeding program was initiated in 2002 and designed to select for growth in even year classes and disease resistance in odd year classes [13]. Stress responsiveness was initially evaluated to determine whether or not it is associated with performance traits affecting aquaculture production efficiency. To this end a total of 584 fish representing 64 families and 9 replicated families (8 fish per family) from the 2002 year class were evaluated for stress response to crowding according to the protocol of Pottinger and Carrick [9] to identify phenotypic variation of 11.6-93.9 ng cortisol/mL and a significant positive association between plasma cortisol response and growth performance [12]. In the 2004 year class, 8 fish from each of the top 15% of the growth selected families were similarly evaluated for stress response to identify a heritability of 46%; mean family values were used to develop a high and low responding (P 1 ) generation [20]. In 2006, high and low responding P 1 fish were crossed to create a F1 generation consisting of 7 full-sib (FS) families including 222 offspring from 12 P 1 parents which were evaluated for their potential for use in understanding the genetic control of this trait [21].

Stress challenge
The stress challenge method for crowding was modified from Pottinger and Carrick [9] as described by Weber and Silverstein [12]. The challenge for the parental generation is described in Weber et al. [20] and for the offspring in Vallejo et al. [17] Briefly, fish were reared as individual families and were exposed to an artificial ambient photoperiod and reared in continuous flow spring water with temperatures between~11.5 and 13.5°C and dissolved oxygen near saturation. Fish were placed in 120 l blue polypropylene tanks at approximately one month of age. Seven weeks prior to the start of the experiment, fish were tagged with passive integrated transponder (PIT) tags (Avid Identification Systems Inc., Norco CA) and split by family into tanks containing eight fish each. Fish in the stress study were fed Zeigler Gold (Zeigler Bros. Inc., Gardners PA) at 2% body weight/ day. The eight fish from each tank were sampled four times at 4-week intervals when the fish were approximately 160 g body weight. Fish were not fed the day of sampling or the afternoon preceding the sampling. For the crowding stress challenge, fish from a single family were netted and transferred from a 120 l tank to a 6 l or 15 l tank and left undisturbed for three hours. The 6 l or 15 l tanks were used in an effort to keep fish densities in the challenge consistent as fish continue to grow throughout the experiment. After three hours of crowding, the fish were then netted and transferred into an anesthesia bath followed by blood collection. Plasma cortisol was measured (ng/ml) by tritium radioimmunoassay following procedures described by Redding et al. [22].

Stress response derived traits
The parents and offspring fish used in this study (sires = 5; dams = 7; offspring = 222) had four repeated measurements of plasma cortisol recorded at about 4-week intervals as described elsewhere [17]. These plasma cortisol measurements are highly variable and heavily impacted by environmental effects with a polygenic model heritability of~0.26 [17]. So, in order to minimize the plasma cortisol variation due to non-genetic factors, and consequently increase the statistical power of QTL detection, and reduce the false positive rate of detected QTL, we derived two stress response phenotypes to use in QTL analyses.
First, we estimated animal breeding value (EBV) using four repeated measurements of plasma cortisol and fitting a mixed inheritance linear model under a Bayesian framework with software iBay version 1.46 [23]. We decided to fit a permanent environmental effect in the repeated measures mixed model analysis to account for the covariance between the records of an individual (i.e., repeated measurements), and capture individual variation across measurements in the EBV computation. The covariates body weight, body length and sexual maturity which had significant effect on the predictive power of the response variable plasma cortisol [17] were also included in the mixed model to minimize the variance in the sampled population.
Second, we also estimated an animal effect using best linear unbiased prediction (BLUP3) and fitting three repeated measurements of plasma cortisol in a multivariate mixed model with the software ASReml version 2.0 [24]. The mixed model also included covariates indicated above that had significant effect on the predictive power of the response variable plasma cortisol. Here, we estimated heritability of each repeated measurement and genetic correlations between repeated measurements of plasma cortisol with ASReml version 2.0 [24], and noticed that measurement at time 1 was different than measurements at time 2, 3 and 4 in heritability and genetic correlations. So, we decided to calculate an index BLUP3 that was weighted by their relative heritability using measurements at time points 2, 3 and 4, respectively.
In order to determine the effect of using adjusted animal effects for mid-parent genetic effect in the QTL analysis, we performed QTL analysis using HS regression interval mapping as outlined above using offspring animal effects that were adjusted for mid-parent genetic effects. The adjustment of offspring animal effects for mid-parent genetic effect allowed accounting for the effects of relatives, and the use of a pure measure of offspring individual's genetic value in the QTL analysis [25,26].
As expected, we determined that the estimated animal effects EBV and BLUP3 had a Pearson's correlation of r = 0.84 (P <0.0001) with SAS Procedure REG [27]. In this study, we used the estimated animal effects EBV and BLUP3 in the QTL analysis. The use of estimated animal effects such as EBV and daughter-yield deviation (DYD) in whole genome QTL scans has been well documented in livestock species [25,[28][29][30][31].

Genotyping and linkage analysis
A panel of 412 microsatellite loci identified from the NCCCWA Genetic Map [19] were used to genotype the parents and offspring of seven FS QTL mapping families (sires = 5; dams = 7; offspring = 222; total marker genotyped fish = 234). Markers were either genotyped using the tailed protocol [32] or by direct fluorescent labelling (with FAM, HEX, or NED) of the forward primer. Primer pairs were obtained from commercial sources (forward primers labelled with FAM or HEX from Alpha DNA, Montreal, Quebec, Canada, or NED from ABI, Foster City, CA, USA). PCR reactions consisted of 12 μl reaction volumes containing 12.5 ng DNA, 1.5-2.5 mM MgCl 2 , 1.0 μM of each primer, 200 μM of dNTPs, 1X manufacturer's reaction buffer and 0.5 units Taq DNA polymerase. Thermal cycling consisted of an initial denaturation at 95°C for 15 min followed by 30 cycles of 95°C for 1 min, annealing temperature for 45 s, 72°C extension for 45 s and a final extension at 72°C for 10 min. PCR products were visualized on agarose gels after staining with ethidium bromide. Markers were grouped in combinations of two or three markers based on differences in fluorescent dye color and amplicon size. Three μl of each PCR product was diluted with 20 μl of water, 1 μl of the diluted sample was added to 12.5 μl of loading mixture made up with 12 μl of HiDi formamide and 0.5 of Genscan 400 ROX internal size standard. Samples were denatured at 95°C for 5 min and kept on ice until loading on an automated DNA sequencer ABI 3730 DNA Analyzer (ABI, Foster City, CA, USA). Output files were analyzed using GeneMapper version 3.7 (ABI, Foster City, CA, USA), formatted using Microsoft Excel and stored in Microsoft Access database. Linkage maps were constructed for each chromosome using marker genotypes from all seven mapping families with the software MULTIMAP version 2.0 [33] using the approach of Rexroad et al [19].

Testing loci for Mendelian segregation distortion
Before QTL analysis, all STR loci were tested for Mendelian segregation distortion (MSD). In outbreed populations, the progeny of an informative QTL mapping family can have any of these marker genotype proportions: 1:1; 1:2:1 and 1:1:1:1. Within individual FS families, the marker genotype counts were performed using a Perl script (Written by G. GAO, unpublished). Then, Chi-Square goodness-of-fit test of marker genotype counts to expected proportions under Mendelian segregation was performed with SAS Procedure FREQ (SAS, 2007) using a default significance level of α = 0.01. Loci with significant MSD were not used in the QTL analysis unless these loci had high quality marker genotypes (i.e., minimum genotyping errors). The rationale to use STR loci that had significant MSD in QTL genome scan is that loci with significant MSD can be linked to viability/survivability and stress response QTL.

QTL analysis using HS regression interval mapping
For each stress response trait, we performed combined sire-family (five HS families) and dam-family (seven HS families) HS regression analysis, separately, using the web-based software GridQTL [34]. This software implements a multi-marker approach of interval mapping in HS families as described by Knott et al. [35]. This method of QTL analysis does not assume the parents had fixed QTL alleles, instead it relaxes the assumption of fixed QTL allele [35]. A QTL with a gene substitution effect is fitted at 1-cM intervals along the chromosome using this one-QTL model, y ij = a i + b i x ij + e ij where y ij is trait score of individual j from sire or dam i; a i is average effect for HS family i; b i is regression coefficient within HS family i (substitution effect of QTL); x ij is conditional probability for individual j within HS family i of inheriting allele 1 (or 2); and e ij is the residual error.
In the HS regression analysis, the likelihood ratio (LR) test statistic was defined as where ln stands for natural logarithm, Î(z) is the likelihood function evaluated at the maximum likelihood estimate (MLE) for the full model that includes polygenic and QTL effects, and Î r (z) is the MLE for the restricted model under which r parameters of the full model are assigned fixed values [36]. The P-value was calculated assuming an F-value distributed with numerator DF equal to the number of sires or dams, and denominator DF equal to the total number of offspring minus twice the number of sires or dams [35]. The chromosome-wide F-value (F ChromWide P=0.05 ) and experiment-wide F-value (F ExperWide P=0.05 ) were estimated using 10,000 permutations with software GridQTL [34]. The genome-wide significance level (P GenomeWide ) for detected QTL was estimated as P GenomeWide = 1 − (1 − P) g [37] where P is the nominal P-value, and g = 380 STR loci used in the HS regression analysis. The QTL with F-value ≥ F ChromWide P=0.05 was defined as suggestive QTL (*); and QTL with F-value ≥ F ExperWide P=0.05 or P GenomeWide ≤ 0.05 was defined as significant QTL (**). The proportion of phenotypic variance explained by the QTL was calculated as h 2 where MSE full and MSE reduced are the mean squared error of the full and reduced model, respectively [35]. The 95% QTL confidence intervals were estimated using 10,000 bootstraps with re-sampling with software GridQTL [34].
Briefly, the HS regression analysis was performed following these steps: First, all chromosomes were genome scanned for EBV/BLUP3 QTL performing sire-family and dam-family regression analysis, separately, using one-QTL model. At this stage: (a) the chromosome-wide significance threshold level (F chromWideP=0.05 ) to declare suggestive QTL was determined using 10,000 permutations; and (b) the 95% confidence interval (CI95) for each suggestive QTL was determined using 10,000 bootstraps with resampling. Second, chromosomes with suggestive QTL were re-scanned for additional QTL using one-QTL model that accounted for the effect of already detected QTL until not detecting more QTL in a chromosome. Third, the chromosomes with two detected QTL were rescanned using two-QTL models. In the chromosomes with two detected QTL: (a) the (F chromWideP=0.05 ) significance level for one QTL at a time was estimated while accounting for the effect of the other detected QTL using 10,000 permutations; and (b) the CI95 for each QTL was determined by fixing alternatively the effect of each detected QTL using 10,000 bootstraps with re-sampling [38]. Fourth, the experiment-wide significance threshold level (F ExperWideP=0.05 ) to declare significant QTL was determined by analyzing altogether marker genotype data from 29 chromosomes using one-QTL models and 10,000 permutations [25,26].

QTL analysis using variance components approach
We also performed a genome scan for genetic loci linked to stress response phenotypes (animal BLUP3 index and EBV for plasma cortisol) using a variance components (VARCOMP) method of QTL analysis with software SOLAR version 4.0 [39]. The aim on using a second linkage-based method of QTL analysis was to reduce the rate of false positive claims of QTL. This robust VAR-COMP method of QTL analysis was performed using the seven FS families altogether (total parents and offspring n =234). This combined family QTL analysis enables accounting for the impact of common environmental effects on QTL mapping which increases the resemblance between full-sibs as they share the same environment during early stages of rainbow trout rearing.
In order to perform multipoint QTL analysis, the multipoint identical by descent (MIBD) relationship matrix for each pair of individuals was estimated with the program LOKI version 2.4.5 [40]. We used LOKI because SOLAR cannot estimate MIBDs in complex pedigrees with double grandparent-grandchild relationships as those used in this study. The VARCOMP approach performs multipoint genome scan at 1-cM intervals. At each testing interval, SOLAR calculates the residual genetic variance or proportion of the total variance due to the polygenic component (h u 2 ); the heritability associated with the QTL or proportion of the total variance due to the QTL (h q 2 ); and the logarithm of odds (LOD) score as LOD = log 10 [L(QTL)/L(polygenic)], where L(QTL) stands for the likelihood of the full model that includes the polygenic and QTL effect. The genome-wide significance level of detected QTL was estimated using this expression P GenomeWide = 1 − (1 − P) g [37] which assumes the use of sparse genetic maps; here P is the nominal Pvalue, and g = 365 STR loci used with VARCOMP analysis (average marker distance = 6.6 cM). The QTL with LOD ≥ 2 was defined as suggestive QTL (*); and the QTL with LOD ≥ 3 or P GenomeWide ≤ 0.05 was defined as significant QTL (**).

QTL analysis using LDLA method
We performed a genome scan for stress response QTL using an LDLA method with the module from the webbased software GridQTL [34,41]. We reasoned that the rate of false positive QTL may be reduced by following up the linkage-based QTL analyses with an LDLA-based method of QTL analysis because the QTL signals have to conform with both LD and LA assumptions [42].
The LR test statistic is defined as indicated above in the section of QTL mapping using HS regression analysis. The LDLA method performs LR tests at 1-cM intervals along the chromosome. Here, the LR test statistic is defined as LR = L(QTL)/L(polygenic) where L (QTL) is the likelihood of the full model that includes the polygenic and QTL effects. The nominal P-value was estimated assuming the LR test statistic follows a chisquare distribution with two degrees of freedom [43]. The genome-wide significance level of detected QTL was also estimated as P GenomeWide = 1 − (1 − P) g [37] where P is the nominal P-value, and g = 365 STR loci used with LDLA analysis. The QTL with LR ≥ 13.82 or nominal P ≤ 0.001 was declared as suggestive QTL (*), and the QTL with P GenomeWide ≤ 0.05 was declared as significant QTL (**). The variance components due to polygenic (σ u 2 ), additive QTL (σ a 2 ), dominant QTL (σ d 2 ) and residual error (σ e 2 ) were used to estimate the proportion of total variance due to the polygenic component [h u , the proportion of total variance due to additive effect QTL [h addQTL , and the proportion of total variance due to dominance effect QTL [h domQTL , respectively.

Correlations among repeated measurements
The heritability, genetic and phenotypic correlations between the four plasma cortisol measurements are presented in Table 1. Heritability estimates for the first measurements were significantly less than the second through fourth measurements, with the heritability increasing over time 1 to 4. The ranges of genetic and phenotypic correlations were much higher when the first measurement is not included.

Genotyping and linkage analysis
The genetic linkage analyses resulted in the construction of genetic maps representing 29 chromosomes and including 380 microsatellite markers whose orders were found to be in good agreement with markers on the NCCCWA genetic map [19]. Marker densities ranged from 4.5 to 10.9 cM per chromosome with a genome average of 6.4 cM spanning 2426 cM. A summary of genetic maps for each chromosome including marker information content across sires and dams is presented in Table 4.

Identification of QTL using half-sib regression interval mapping
Results of QTL analyses for EBV and BLUP3 using halfsib regression interval mapping (HS) are presented in 0QTL; DF=5 =3.58; LR = 32.50) when fitting two-QTL models with sire HS families ( Table 2). The estimated sire QTL effects with their corresponding significance level suggest that Sire 5 and 3 were likely segregating the QTL detected on Omy16 and Sex, respectively (Additional file 2). The measurement unit is expressed in ng/mL of plasma cortisol. 2 The heritability estimates with their corresponding S.E. are bold text highlighted and shown in the main diagonal. 3 The genetic correlation estimates with their corresponding S.E. are shown in the lower diagonal. 4 The phenotypic correlation estimates with their corresponding S.E. are shown in the upper diagonal. The likelihood ratio (LR) test statistic. 2 Asymptotically F-test statistic with the degrees of freedom(DF) being the number of sires or dams included for the numerator, and the total number of offspring minus twice the number of sires or dams for the denominator. The QTL with F-value ≥ F ChromWide P=0.05 was defined as suggestive QTL (*); and QTL with F-value ≥ F ExperWide P=0.05 or P GenomeWide ≤ 0.05 was defined as significant QTL (**). 3 Chromosome-wide F-value at P=0.05 was estimated using 10,000 permutations with software GridQTL. 4 Experiment-wide F-value at P=0.05 was estimated using 10,000 bootstraps with re-sampling with software GridQTL. 5 The nominal P-value was calculated assuming an F-value distributed with numerator DF equal to the number of sires or dams, and denominator DF equal to the total number of offspring minus twice the number of sires or dams. 6 The genome-wide significance level for detected QTL. 7 The proportion of phenotypic variance explained by the QTL. 8 The QTL flanking markers OMM1772, OMM1130, BX913059 and OMM5153 had significant Mendelian segregation distortion (P <0.01) in families 2, 4, 5 and 6, respectively. 9 The Omy16 QTL had significant statistical support for two QTL segregating (F 2QTL vs. 0QTL; DF=10 =3.58; F 1QTL vs. 0QTL; DF=5 =3.60; LR 2QTL vs. 0QTL = 33.62) for BLUP3 when fitting a two-QTL model with sire HS families. 10 The Sex chromosome QTL had significant statistical support for two QTL segregating (F 2QTL vs. 0QTL; DF=10 =3.45; F 1QTL vs. 0QTL; DF=5 =3.10; LR 2QTL vs. 0QTL = 32.50) for EBV when fitting a two-QTL model with sire HS families.

Identification of QTL using a variance components approach
Results of the QTL analyses using the variance component (VC) approach are presented in Table 3. Two suggestive QTL explaining 24 and 43% of the phenotypic variance were detected for EBV (Omy12 and Sex) and a significant QTL explaining 52% of the phenotypic variation was detected for BLUP3BLUP33 (Omy16). Both of the QTL for EBV overlap with similar QTL detected with HS regression interval mapping ( Table 2). The significant QTL for BLUP3 on Omy16 was also determined to be significant by HS regression analyses ( Table 2).

Discussion
Identifying genetic variation in mechanisms underlying complex traits such as response to stress may improve our ability to mitigate the negative effects of stressors on aquaculture production through selective breeding or management practices. To this end we sought to identify QTL affecting response to crowding in seven families from a broodstock population under selection for increased growth rate for a single generation. The families evaluated in this study were not bred specifically for QTL analyses; therefore multiple approaches to QTL discovery were performed as a means of validation and to identify an inclusive list of chromosomes with genes affecting stress phenotypes. Although the 7 families including 5 sires, 7 dams and a total of 222 offspring is low for estimating heritability and genetic correlations, we were able to detect a single significant and 7 suggestive QTL. Although we have detected locations of QTL at the sub-chromosomal level on the genetic maps, it is possible that multiple QTL on a single chromosome are actually the same QTL identified in different families and map positions differ because of imbalances in female:male recombination ratios [44]. Therefore, we take a conservative approach and focus our discussion at the level of the chromosome. These results are summarized in Table 4, where suggestive and significant QTL are tabulated according to each analysis and trait.

Comparisons of EBV and BLUP3
In estimating the heritability and genetic correlation between repeated plasma cortisol measurements, we Logarithm of odds (LOD) score was calculated as LOD = log 10 [L(QTL)/L(polygenic)] where L = likelihood of the model. QTL with LOD ≥ 2 was defined as suggestive QTL (*); and QTL with LOD ≥ 3 or P GenomeWide ≤ 0.05 was defined as significant QTL (**). 2 The genome-wide significance level for detected QTL was estimated as P GenomeWide = 1 − (1 − P) g , where P is the nominal P-value, and g = 365 STR loci used with variance components method of QTL analysis. 3 h u 2 is the residual genetic variance or proportion of the total variance due to the polygenic component. 4 h q 2 is the heritability associated with the QTL or proportion of the total variance due to the QTL. 5 The QTL flanking markers OMM5106a, OMM1408, OMM1130, BX913059 and OMM5143 had significant Mendelian segregation distortion (P <0.01) in families 2, 3, 4, 5 and 6, respectively.
observed that measurement at time point 1 was markedly different in terms of heritability and genetic correlations than those measurements at times 2, 3 and 4 which were all about a month apart. This may indicate that physiological response to stress is altered after the initial challenge and genes participating in a first response to a stressor might differ in part from that of subsequent exposures. Overall our interest is in the responses of rainbow trout to the stressors of the aquaculture production environment throughout the entire grow-out cycle, therefore we sought to identify a robust trait which characterizes the chronic response as closely as possible. Given our observation that the measurements taken at the first time point was somewhat of an outlier; we developed BLUP3 using only measurements at time points 2, 3 and 4 weighted by their corresponding relative heritability. Our thought was there might be something unique and biologically interesting about the first measurement such that not including it would result in the detection of different sets of QTL than when all four measurements are Marker information content combines the multilocus probability of individuals inheriting allele 1 or 2 from the common parent with marker segregation distortion [35]. [2][3][4] The HS, VC, and LDLA columns show chromosomes containing suggestive (*) and significant (**) QTL for each trait (EBV, BLUP3) identified by half-sib regression interval mapping, variance component analysis, or linkage disequilibrium and linkage analysis, respectively. included (EBV). In fact, the HS analysis resulted in the detection of QTL for both traits on Omy19, with five chromosomes having QTL specific for EBV and two specific for BLUP3. Three QTL detected by VC analysis were trait specific. Eventually we hope to realize the impacts these QTL have on production traits and determine if one or the other is a more accurate predictor of long term response to stress.

Comparisons of QTL detected by HS, LDLA, and VC
The three methods of QTL analysis (HS, LDLA and VC) used in this study detected unique sets of QTL along the scanned rainbow trout genome due to several reasons. First, each of these methods uses a different algorithm of QTL analysis. The HS and VC methods are linkagebased methods, and the LDLA method exploits simultaneously linkage and LD information from the sampled where ln stands for natural logarithm, Î(z) is the likelihood function evaluated at the maximum likelihood estimate (MLE) for the full model that includes polygenic and QTL effects, and Î r (z) is the MLE for the restricted model under which r parameters of the full model are assigned fixed values. The QTL was declared as suggestive if it had a LR ≥ 13.82 or nominal P ≤ 0.001 (*), and as significant if it had a P GenomeWide ≤0.05 (**). 2 The nominal P-value was estimated assuming the LR test statistic follows a chi-square distribution with two degrees of freedom. 3 The genome-wide significance level for detected QTL was estimated as P GenomeWide = 1 − (1 − P) g [37] where P is the nominal P-value, and g = 365 STR loci used with LDLA method of QTL analysis. 4 Logarithm of odds (LOD) score was calculated as LOD = log 10 [L(QTL)/L(polygenic)] where L = likelihood of the model. 5 h u 2 is the residual genetic variance or proportion of the total variance due to the polygenic component. families. The HS method uses a least-square based regression interval mapping approach, and the VC method is a variance components based approach. From others and our extensive testing of these computer packages, we noticed that the VC algorithm is more conservative and robust to violations of assumptions on trait multivariate normal distribution with low rate of type I errors in the analysis of complex traits [39,45]. In contrast, the HS and LDLA methods seem to be more liberal and powerful than the VC method with a trade-off of more type I errors in the HS and LDLA engines. Clearly, the test-statistic values were generally larger for the HS and LDLA methods than for those of the VC method. Comparing the test-statistics of these three methods of QTL analysis, we noticed that they display a similar trend along the genome although at different scale of magnitude: higher scale of magnitude for the HS and LDLA method, and lower scale of magnitude for the VC method; consequently reaching most often statistical significance with the HS and LDLA method, and less often with the VC method. Coupling with complementing use of EBV and BLUP3 to estimate the effects of initial vs. repeated stressors, we believe that this strategy of using three different algorithms of QTL mapping allowed us to reduce the false positive rate of detected QTL in this genome scan. To this end across the three analyses significant QTL were detected on Omy16 for BLUP3; and at least suggestive QTL for EBV were detected on Omy12 and Sex. However, we acknowledge that there is also a likelihood of false negatives as these QTL detection methods are not optimized for our mapping families, and that lack of validation across analyses does not necessarily invalidate QTL detected in a single analyses. LDLA detected QTL for EBV on 22 chromosomes and for BLUP3 on 5 chromosomes, identifying a high degree of false positives. The robustness of LDLA software has not yet been extensively tested, and the LDLA developers warned about the potential pitfalls [41]. They indicated that a large significant result obtained with LDLA and a very insignificant result obtained with LA method may indicate a false positive result. So, here the LDLA results were handled with caution and compared with HS and VC results.
Drew et al. [11] compared response to netting stress in three clonal lines of rainbow trout differing in their histories of domestication and observed significant variation in plasma cortisol concentrations. Genotyping of doubled haploid offspring created from crossing two of those strains identified two QTL with significant and antagonistic additive effects that explained 43% of the phenotypic variance, one of which overlapped with a QTL for juvenile body mass. Our results from HS and VC do not overlap with the findings of that study. However, given that the stressor was not exactly the same in the two experiments and the differences in populations, family structures, and analytical methods; it is possible there exists some overlap in QTL affecting stress response in the populations from both studies that were just not detected due to differences in experimental designs. The number of QTL identified across studies highlights the large degree of genetic variation affecting stress responses in rainbow trout.

Comparing results with CSA
Performing CSA for stress response traits, we found that more than one major gene with dominant cortisoldecreasing alleles and also major genes with additive effect of −42 ng/mL were segregating in the NCCCWA rainbow trout broodstock population [17]. The results found in this QTL genome scan partially validate the previous CSA report. Here, using linkage-based methods of QTL analysis (HS and VC), we identified eight suggestive/significant QTL that were validated with LDLA method. These detected QTL explained a wide range of phenotypic trait variation (17-52%). The additive effect of suggestive/significant QTL detected with HS regression analysis ranged from -17.98 to 17.0 ng/mL (Additional file 2). However, the additive effect of significant QTL detected on Omy16 was -17.9 ng/mL which even as a cumulative effect was remarkably smaller than the predicted effect of the CSA QTL.
By design, the HS regression analysis performed with GridQTL [34] can make inferences only about additive genetic effects. Although we performed combined FS family analysis with the VC method, the SOLAR software was implemented to fit only additive effect QTL models. In contrast, the combined FS family analysis performed with the LDLA method allowed testing for additive and dominance QTL effects. From the QTL detected with the LDLA method (Table 5) about 1/3 and 2/3 of detected QTL had additive and dominance effects, respectively.

Conclusions
We report the detection of multiple QTL for stress response as defined by Pottinger and Carrick [18] and modified by Weber and Silverstein [12] in a mapping population sampled from broodstock under selection for growth. Suggestive and significant QTL affecting stress responses were detected on eight chromosomes in a survey of seven families from a single broodstock population. These QTL are largely different from others previously identified for plasma cortisol level in response to stress in rainbow trout, implying that many genetic variants may affect this trait which may independently affect other traits. Also, mapping QTL for EBV and BLUP3 produced trait specific QTL which indicate the significance of the first measurement in defining traits associated with stress response. Validation of these QTL in other families and populations followed by fine mapping will lead towards the identification of genes affecting stress response and may influence approaches to selection regarding these economically important traits.