- Research article
- Open Access
QTL affecting stress response to crowding in a rainbow trout broodstock population
BMC Genetics volume 13, Article number: 97 (2012)
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.
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.
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.
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. 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, sub-optimal 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–7]. Fish respond to these stressors in similar patterns to those exhibited by terrestrial vertebrates, by stimulating the hypothalamic-sympathetic-chromaffin 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. 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.
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. To date selection has focused on characteristics associated with growth and disease resistance[14–16]. Stress has been shown to affect both of these traits. 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. 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. 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 as modified by Weber and Silverstein, and conducted a genome scan with 380 microsatellite markers selected from the NCCCWA genetic map. 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.
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.
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. 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 to identify phenotypic variation of 11.6-93.9 ng cortisol/mL and a significant positive association between plasma cortisol response and growth performance. 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 (P1) generation. In 2006, high and low responding P1 fish were crossed to create a F1 generation consisting of 7 full-sib (FS) families including 222 offspring from 12 P1 parents which were evaluated for their potential for use in understanding the genetic control of this trait.
The stress challenge method for crowding was modified from Pottinger and Carrick as described by Weber and Silverstein. The challenge for the parental generation is described in Weber et al. and for the offspring in Vallejo et al. 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..
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. These plasma cortisol measurements are highly variable and heavily impacted by environmental effects with a polygenic model heritability of ~0.26. 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. 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 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. 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, 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. 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–31].
Genotyping and linkage analysis
A panel of 412 microsatellite loci identified from the NCCCWA Genetic Map 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 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 MgCl2, 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 using the approach of Rexroad et al.
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. This software implements a multi-marker approach of interval mapping in HS families as described by Knott et al.. This method of QTL analysis does not assume the parents had fixed QTL alleles, instead it relaxes the assumption of fixed QTL allele. 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. 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. The chromosome-wide F-value (F ChromWide P=0.05 ) and experiment-wide F-value (FExperWide P=0.05) were estimated using 10,000 permutations with software GridQTL. The genome-wide significance level (P GenomeWide ) for detected QTL was estimated as P GenomeWide = 1 − (1 − P)g 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 where MSE full and MSE reduced are the mean squared error of the full and reduced model, respectively. The 95% QTL confidence intervals were estimated using 10,000 bootstraps with re-sampling with software GridQTL.
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 (FchromWideP=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 re-sampling. 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 re-scanned using two-QTL models. In the chromosomes with two detected QTL: (a) the (FchromWideP=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. Fourth, the experiment-wide significance threshold level (FExperWideP=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. 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 VARCOMP 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. 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 = log10L(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 which assumes the use of sparse genetic maps; here P is the nominal P-value, 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 web-based 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.
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 chi-square distribution with two degrees of freedom. The genome-wide significance level of detected QTL was also estimated as P GenomeWide = 1 − (1 − P)g 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 2 = σ u 2/(σ u 2 + σ a 2 + σ d 2 + σ e 2)], the proportion of total variance due to additive effect QTL h addQTL 2 = σ a 2/(σ u 2 + σ a 2 + σ d 2 + σ e 2)], and the proportion of total variance due to dominance effect QTL h domQTL 2 = σ d 2/(σ u 2 + σ a 2 + σ d 2 + σ e 2)], respectively.
Correlations among repeated measurements
The heritability, genetic and phenotypic correlations between the four plasma cortisol measurements are presented in Table1. 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.
Loci testing for Mendelian segregation distortion
Testing 346 short tandem repeat (STR) loci for Mendelian segregation distortion (MSD) resulted in about 1%, 6%, 4%, 2%, 2%, 3% and 2% of markers having significant MSD (P ≤ 0.01) in families 1, 2, 3, 4, 5, 6 and 7, respectively (Additional file1). From all markers with significant MSD, four of these markers were flanking suggestive QTL for stress response: OMM1130 (Omy12); OMM5153 (Omy14); BX913059 (Omy22); and OMM1772 (OmySex) (Tables2 and3).
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. 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 Table4.
Identification of QTL using half-sib regression interval mapping
Results of QTL analyses for EBV and BLUP3 using half-sib regression interval mapping (HS) are presented in Table2. Six suggestive QTL were detected for EBV (Omy10, 12 14, 19, 22 and Sex) which each explain 17 - 30% of the phenotypic variance. For BLUP3, suggestive (Omy6 and 19) and one significant (Omy16) QTL were detected that explain 18-32% of the phenotypic variance. The only overlapping QTL between traits are suggestive QTLs on Omy19. Multiple QTL were detected on Omy12, 16, 19 and Sex.
On chromosomes with two detected QTL using one-QTL models, we found significant statistical support for two QTL segregating on Omy16 (BLUP3BLUP33; F2QTL vs. 0QTL; DF=10=3.58; F1QTL vs. 0QTL; DF=5=3.60; LR = 33.62) and Sex chromosome (EBV; F2QTL vs. 0QTL; DF=10=3.45; F1QTL vs. 0QTL; DF=5=3.58; LR = 32.50) when fitting two-QTL models with sire HS families (Table2). 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 file2).
Identification of QTL using a variance components approach
Results of the QTL analyses using the variance component (VC) approach are presented in Table3. 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 (Table2). The significant QTL for BLUP3 on Omy16 was also determined to be significant by HS regression analyses (Table2).
Identification of QTL using linkage disequilibrium and linkage analysis
Results of efforts to detect QTL for EBV and BLUP3 using linkage disequilibrium and linkage analysis (LDLA) are presented in Table5. Three suggestive (Omy14, 18, and 21) and 19 significant (Omy3-13, 16, 17, 19, and 25-Sex) QTL were detected for EBV; three suggestive (Omy3, 13, and Sex) and two significant (Omy5 and 16) QTL were detected for BLUP3.
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. Therefore, we take a conservative approach and focus our discussion at the level of the chromosome. These results are summarized in Table4, 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 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 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 linkage-based methods, and the LDLA method exploits simultaneously linkage and LD information from the sampled 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. 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. 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 cortisol-decreasing alleles and also major genes with additive effect of −42 ng/mL were segregating in the NCCCWA rainbow trout broodstock population. 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 file2). 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 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 (Table5) about 1/3 and 2/3 of detected QTL had additive and dominance effects, respectively.
We report the detection of multiple QTL for stress response as defined by Pottinger and Carrick and modified by Weber and Silverstein 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.
Liu ZJ, Cordes JF: DNA marker technologies and their applications in aquaculture genetics. Aquaculture. 2004, 238 (1–4): 1-37.
Øverli Ø, Sørensen C, Kiessling A, Pottinger TG, Gjøen HM: Selection for improved stress tolerance in rainbow trout (Oncorhynchus mykiss) leads to reduced feed waste. Aquaculture. 2006, 261 (2): 776-781. 10.1016/j.aquaculture.2006.08.049.
Pickering AD: Rainbow trout husbandry: management of the stress response. Aquaculture. 1992, 100 (1–3): 125-139.
Wendelaar Bonga SE: The stress response in fish. Physiol Rev. 1997, 77 (3): 591-625.
Lefèvre F, Bugeon J, Aupérin B, Aubin J: Rearing oxygen level and slaughter stress effects on rainbow trout flesh quality. Aquaculture. 2008, 284 (1–4): 81-89.
Petri Hoskonen JP: Effects of repeated handling, with or without anaesthesia, on feed intake and growth in juvenile rainbow trout, Oncorhynchus mykiss (Walbaum). Aquacult Res. 2006, 37 (4): 409-415. 10.1111/j.1365-2109.2005.01448.x.
Olsen RE, Sundell K, Mayhew TM, Myklebust R, Ringø E: Acute stress alters intestinal function of rainbow trout, Oncorhynchus mykiss (Walbaum). Aquaculture. 2005, 250 (1–2): 480-495.
Bonga SEW: The stress response in fish. Physiol Rev. 1997, 77 (3): 591-625.
Pottinger TG, Carrick TR: Modification of the plasma cortisol response to stress in rainbow trout by selective breeding. Gen Comp Endocrinol. 1999, 116 (1): 122-132. 10.1006/gcen.1999.7355.
Fevolden S-E, Røed KH, Fjalestad KT: Selection response of cortisol and lysozyme in rainbow trout and correlation to growth. Aquaculture. 2002, 205 (1–2): 61-75.
Drew RE, Schwabl H, Wheeler PA, Thorgaard GH: Detection of QTL influencing cortisol levels in rainbow trout (Oncorhynchus mykiss). Aquaculture. 2007, 272 (Supplement 1): S183-S194.
Weber GM, Silverstein JT: Evaluation of a Stress Response for Use in a Selective Breeding Program for Improved Growth and Disease Resistance in Rainbow Trout. North Am J Aquaculture. 2007, 69 (1): 69-79. 10.1577/A05-103.1.
Silverstein JT, King T, Rexroad CE: Genetic Variation Measured by Microsatellites Among Three Strains of Domesticated Rainbow Trout. Aquacult Res. 2004, 35: 40-48. 10.1111/j.1365-2109.2004.00979.x.
Silverstein JT, Vallejo RL, Palti Y, Leeds TD, Rexroad CE, Welch TJ, Wiens GD, Ducrocq V: Rainbow trout resistance to bacterial cold-water disease is moderately heritable and is not adversely correlated with growth. J Anim Sci. 2009, 87 (3): 860-867.
Leeds T, Silverstein J, Weber G, Vallejo R, Palti Y, Rexroad CE, Evenhuis J, Hadidi S, Welch T, Wiens G: Response to selection for bacterial cold water disease resistance in rainbow trout. J Anim Sci. 2010, 88: 1936-1946. 10.2527/jas.2009-2538.
Silverstein JT, Hostuttler M, Blemings KP: Strain differences in feed efficiency measured as residual feed intake in individually reared rainbow trout, Oncorhynchus mykiss (Walbaum). Aquacult Res. 2005, 36: 704-711. 10.1111/j.1365-2109.2005.01278.x.
Vallejo RL, Rexroad CE: Silverstein JT, Janss LL, Weber GM: Evidence of major genes affecting stress response in rainbow trout using Bayesian methods of complex segregation analysis. J Anim Sci. 2009, 87 (11): 3490-3505. 10.2527/jas.2008-1616.
Pottinger TG, Carrick TR: Stress responsiveness affects dominant-subordinate relationships in rainbow trout. Horm Behav. 2001, 40 (3): 419-427. 10.1006/hbeh.2001.1707.
Rexroad CE, Palti Y, Gahr SA, Vallejo RL: A second generation genetic map for rainbow trout (Oncorhynchus mykiss). BMC Genet. 2008, 9: 74-
Weber GM, Vallejo RL, Lankford SE, Silverstein JT, Welch TJ: Cortisol Response to a Crowding Stress: Heritability and Association with Disease Resistance to Yersinia ruckeri in Rainbow Trout. North Am J Aquaculture. 2008, 70: 425-433. 10.1577/A07-059.1.
Vallejo RL CER, Silverstein JT, Weber GM: Multivariate Complex Segregation Analysis Provides Compelling Evidence of Major Genes Underlying Genetic Variation for Stress Response in Rainbow Trout. J Anim Sci. 2009, 87: 3490-3505. 10.2527/jas.2008-1616.
Redding JM, Schreck CB, Birks EK, Ewing RD: Cortisol and its effects on plasma thyroid hormone and electrolyte concentrations in fresh water and during seawater acclimation in yearling coho salmon, Oncorhynchus kisutch. Gen Comp Endocrinol. 1984, 56 (1): 146-155. 10.1016/0016-6480(84)90071-6.
Janss LLG: iBay manual version 1.46. 2008, Leiden, Netherlands: Janss Biostatistics, P.O. Box 535
Gilmour AR, Gogel BJ, Cullis BR, Thompson R: ASReml User Guide Release 2.0. 2006, Hemel Hempstead: VSN International Ltd
de Koning DJ, Windsor D, Hocking PM, Burt DW, Law A, Haley CS, Morris A, Vincent J, Griffin H: Quantitative trait locus detection in commercial broiler lines using candidate regions. J Anim Sci. 2003, 81 (5): 1158-1165.
Van Kaam JBCMH, van Arendonk JAM, Groenen MAM, Bovenhuis H, Vereijken ALJ, Crooijmans R, van der Poel JJ, Veenendaal A: Whole genome scan for quantitative trait loci affecting body weight in chickens using a three generation design. Livest Prod Sci. 1998, 54: 133-150. 10.1016/S0301-6226(97)00171-1.
SAS: SAS 9.1.3 Help and Documentation. 2007, Cary, NC: SAS Institute Inc
Weller JI: Quantitative trait analysis in animals. 2009, Cambridge, MA: CAB, 2
Freyer G, Kuhn C, Weikard R, Zhang Q, Mayer M, Hoeschele I: Multiple QTL on chromosome six in dairy cattle affecting yield and content traits. J Anim Breed Genet. 2002, 119 (2): 69-82. 10.1046/j.1439-0388.2002.00322.x.
Freyer G, Stricker C, Kuhn C: Comparison of estimated breeding values and daughter yield deviations used in segregation and linkage analyses. Czech J Anim Sci. 2002, 47 (6): 247-252.
Thomsen H, Reinsch N, Xu N, Looft C, Grupe S, Kuhn C, Brockmann GA, Schwerin M, Leyhe-Horn B, Hiendleder S: Comparison of estimated breeding values, daughter yield deviations and de-regressed proofs within a whole genome scan for QTL. J Anim Breed Genet. 2001, 118 (6): 357-370. 10.1046/j.1439-0388.2001.00302.x.
Boutin-Ganache I, Raposo M, Raymond M, Deschepper CF: M13-tailed primers improve the readability and usability of microsatellite analyses performed with two different allele-sizing methods. Biotechniques. 2001, 31 (1): 24-26, 28.
Matise TC, Perlin M, Chakravarti A: Automated construction of genetic linkage maps using an expert system (MultiMap): a human genome linkage map. Nat Genet. 1994, 6 (4): 384-390. 10.1038/ng0494-384.
Seaton G, Haley CS, Knott SA, Kearsey M, Visscher PM: QTL Express: mapping quantitative trait loci in of simple and complex pedigrees. Bioinformatics. 2002, 18 (2): 339-340. 10.1093/bioinformatics/18.2.339.
Knott SA, Elsen JM, Haley CS: Methods for multiple-marker mapping of quantitative trait loci in half-sib populations. Theor Appl Genet. 1996, 93 (1–2): 71-80.
Lynch M, Walsh B: Genetics and analysis of quantitative traits. 1997, Sunderland, MA: Sinauer Associates Inc, 1
Ott J: Analysis of Human Genetic Linkage. 1999, Baltimore: The John Hopkins Unversity Press, 3
Leach RJ, Craigmile SC, Knott SA, Williams JL, Glass EJ: Quantitative trait loci for variation in immune response to a Foot-and-Mouth Disease virus peptide. BMC Genet. 2010, 11: 107-(7 December 2010)
Almasy L, Blangero J: Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998, 62 (5): 1198-1211. 10.1086/301844.
Heath SC: Markov chain Monte Carlo segregation and linkage analysis for oligogenic models. Am J Hum Genet. 1997, 61 (3): 748-760. 10.1086/515506.
Hernandez-Sanchez J, Grunchec JA, Knott S: A web application to perform linkage disequilibrium and linkage analyses on a computational grid. Bioinformatics. 2009, 25 (11): 1377-1383. 10.1093/bioinformatics/btp171.
Meuwissen TH, Karlsen A, Lien S, Olsaker I, Goddard ME: Fine mapping of a quantitative trait locus for twinning rate using combined linkage and linkage disequilibrium mapping. Genetics. 2002, 161 (1): 373-379.
Gong YF, Lu X, Wang ZP, Hu F, Luo YR, Cai SQ, Qi CM, Li S, Niu XY, Qiu XT: Detection of quantitative trait loci affecting haematological traits in swine via genome scanning. BMC Genet. 2010, 11 (1): 56-
Sakamoto T, Danzmann RG, Gharbi K, Howard P, Ozaki A, Khoo SK, Woram RA, Okamoto N, Ferguson MM, Holm LE: A microsatellite linkage map of rainbow trout (Oncorhynchus mykiss) characterized by large sex-specific differences in recombination rates. Genetics. 2000, 155 (3): 1331-1345.
Williams JT, Blangero J: Power of variance component linkage analysis to detect quantitative trait loci. Ann Hum Genet. 1999, 63: 545-563. 10.1046/j.1469-1809.1999.6360545.x.
We would like to acknowledge the following people for providing excellent technical assistance including Roseanna Long, Kristy Shewbridge, and Brian Smith for genotyping; Lisa Radler, Jill Birkett, and David Payne for phenotyping fish for stress response; and Josh Kretzer, Jim Everson, Kyle Jenkins, Jenea McGowan, Kevin Melody for animal care. We are very grateful to Guangtu Gao for writing several PERL scripts used in genetic data analysis of this study. We are also very grateful to Dirk-Jan de Koning for generously sharing his expertise in methods of QTL analysis. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture. USDA is an equal opportunity provider and employer.
CER and GMW designed overall study including the breeding scheme, GMW executed the stress challenge and phenotyping, RLV performed QTL analyses, and SL and YP participated in genotyping and linkage analysis. All authors read and approved the final manuscript.
Caird E Rexroad III, Roger L Vallejo contributed equally to this work.
Electronic supplementary material
Additional file 1:STR loci testing for Mendelian segregation distortion (MSD). In outbred 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 each of seven FS families, the marker genotype counts were performed using a Perl script (Written by G. GAO, unpublished). Then , the STR loci were tested for MSD using Chi-Square goodness-of-fit test of marker genotype counts to expected proportions under Mendelian segregation with SAS Procedure FREQ (SAS, 2007) using a default significance level of α = 0.01. (XLSX 132 KB)
Additional file 2: Summary of QTL location and effect for stress response traits using combined sire (and dam) half-sib family regression analysis performed with software GridQTL (Seaton et al. 2006). The QTL effect is expressed in ng/mL of plasma cortisol, and the allele substitution effect for each parent was tested using a one-sided t-test (testing absolute t-values) with one DF. Within each trait QTL group, the sire or dam parent that is most likely segregating the QTL allele is indicated with an asterisk (*), and the parent with the lowest P-value is bold text highlighted. The average QTL location was determined using 10000 bootstraps with re-sampling. (DOCX 32 KB)
About this article
Cite this article
Rexroad, C.E., Vallejo, R.L., Liu, S. et al. QTL affecting stress response to crowding in a rainbow trout broodstock population. BMC Genet 13, 97 (2012). https://doi.org/10.1186/1471-2156-13-97
- Quantitative Trait Locus
- Plasma Cortisol
- Quantitative Trait Locus Analysis
- Quantitative Trait Locus Effect
- Significant Quantitative Trait Locus