Repeatability of cortisol stress response in the European sea bass (Dicentrarchus labrax) and transcription differences between individuals with divergent responses

Understanding the stress responses of organisms is of importance in the performance and welfare of farmed animals, including fish. Especially fish in aquaculture commonly face stressors, and better knowledge of their responses may assist in proper husbandry and selection of breeding stocks. European sea bass (Dicentrarchus labrax), a species with high cortisol concentrations, is of major importance in this respect. The main objectives of the present study were to assess the repeatability and consistency of cortisol stress response and to identify differences in liver transcription profiles of European sea bass individuals, showing a consistent low (LR) or high (HR) cortisol response. The progeny of six full sib families was used, and sampled for plasma cortisol after an acute stress challenge once per month, for four consecutive months. Results suggest that cortisol responsiveness was a repeatable trait with LR and HR fish showing low or high resting, free and post-stress cortisol concentrations respectively. Finally, the liver transcription profiles of LR and HR fish showed some important differences, indicating differential hepatic regulation between these divergent phenotypes. These transcription differences were related to various metabolic and immunological processes, with 169 transcripts being transcribed exclusively in LR fish and 161 exclusively in HR fish.

Quantitative trait loci for magnitude of the plasma cortisol response to confinement in rainbow trout

Better understanding of the mechanisms underlying interindividual variation in stress responses and their links with production traits is a key issue for sustainable animal breeding. In this study, we searched for quantitative trait loci (QTL) controlling the magnitude of the plasma cortisol stress response and compared them to body size traits in five F2 full-sib families issued from two rainbow trout lines divergently selected for high or low post-confinement plasma cortisol level. Approximately 1000 F2 individuals were individually tagged and exposed to two successive acute confinement challenges (1 month interval). Post-stress plasma cortisol concentrations were determined for each fish. A medium density genome scan was carried out (268 markers, overall marker spacing less than 10 cM). QTL detection was performed using qtlmap software, based on an interval mapping method (http://www.inra.fr/qtlmap). Overall, QTL of medium individual effects on cortisol responsiveness (<10% of phenotypic variance) were detected on 18 chromosomes, strongly supporting the hypothesis that control of the trait is polygenic. Although a core array of QTL controlled cortisol concentrations at both challenges, several QTL seemed challenge specific, suggesting that responses to the first and to a subsequent exposure to the confinement stressor are distinct traits sharing only part of their genetic control. Chromosomal location of the steroidogenic acute regulatory protein (STAR) makes it a good potential candidate gene for one of the QTL. Finally, comparison of body size traits QTL (weight, length and body conformation) with cortisol-associated QTL did not support evidence for negative genetic relationships between the two types of traits.

Assessment of genetic correlation between bacterial cold water disease resistance and spleen index in a domesticated population of rainbow trout: identification of QTL on chromosome Omy19

Selective breeding of animals for increased disease resistance is an effective strategy to reduce mortality in aquaculture. However, implementation of selective breeding programs is limited by an incomplete understanding of host resistance traits. We previously reported results of a rainbow trout selection program that demonstrated increased survival following challenge with Flavobacterium psychrophilum, the causative agent of bacterial cold water disease (BCWD). Mechanistic study of disease resistance identified a positive phenotypic correlation between post-challenge survival and spleen somatic-index (SI). Herein, we investigated the hypothesis of a genetic correlation between the two traits influenced by colocalizing QTL. We evaluated the inheritance and calculated the genetic correlation in five year-classes of odd- and even-year breeding lines. A total of 322 pedigreed families (n = 25,369 fish) were measured for disease resistance, and 251 families (n = 5,645 fish) were evaluated for SI. Spleen index was moderately heritable in both even-year (h(2)  = 0.56±0.18) and odd-year (h(2)  = 0.60±0.15) lines. A significant genetic correlation between SI and BCWD resistance was observed in the even-year line (rg  = 0.45±0.20, P = 0.03) but not in the odd-year line (rg  = 0.16±0.12, P = 0.19). Complex segregation analyses of the even-year line provided evidence of genes with major effect on SI, and a genome scan of a single family, 2008132, detected three significant QTL on chromosomes Omy19, 16 and 5, in addition to ten suggestive QTL. A separate chromosome scan for disease resistance in family 2008132 identified a significant BCWD QTL on Omy19 that was associated with time to death and percent survival. In family 2008132, Omy19 microsatellite alleles that associated with higher disease resistance also associated with increased spleen size raising the hypothesis that closely linked QTL contribute to the correlation between these traits. To our knowledge, this is the first estimation of spleen size heritability and evidence for genetic linkage with specific disease resistance in a teleost fish.

Quantitative trait loci affecting response to crowding stress in an F(2) generation of rainbow trout produced through phenotypic selection

Selective breeding programs for salmonids typically aim to improve traits associated with growth and disease resistance. It has been established that stressors common to production environments can adversely affect these and other traits which are important to producers and consumers. Previously, we employed phenotypic selection to create families that exhibit high or low plasma cortisol concentrations in response to crowding stress. Subsequent crosses of high × low phenotypes founded a multigenerational breeding scheme with the aim of dissecting the genetic basis for variation underlying stress response through the identification of quantitative trait loci (QTL). Multiple methods of QTL analyses differing in their assumptions of homozygosity of the causal alleles in the grandparental generation yielded similar results in the F1 generation, and the analysis of two stress response phenotype measurement indexes were highly correlated. In the current study, we conducted a genome scan with microsatellites to detect QTL in the F2 generation of two families created through phenotypic selection and having larger numbers of offspring than families screened in the previous generation. Seven suggestive and three significant QTL were detected, seven of which were not previously detected in the National Center for Cool and Cold Water Aquaculture germplasm, bringing the total number of chromosomes containing significant and suggestive stress response QTL to 4 and 15, respectively. One significant QTL which peaks at 7 cM on chromosome Omy12 spans 12 cM and explains 25 % of the phenotypic variance in family 2008052 particularly warrants further investigation. Five QTL with significant parent-of-origin effects were detected in family 2008052, including two QTL on Omy12. The 95 % confidence intervals for the remaining QTL we detected were broad, requiring validation and fine mapping with other genotyping approaches and mapping strategies. These results will facilitate identification of potential casual alleles that can be employed in strategies aimed at better understanding the genetic and physiological basis of stress responses to crowding in rainbow trout aquaculture production.

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.

Genetics of behavioural adaptation of livestock to farming conditions

Behavioural adaptation of farm animals to environmental changes contributes to high levels of production under a wide range of farming conditions, from highly controlled indoor systems to harsh outdoor systems. The genetic variation in livestock behaviour is considerable. Animals and genotypes with a larger behavioural capacity for adaptation may cope more readily with varying farming conditions than those with a lower capacity for adaptation. This capacity should be exploited when the aim is to use a limited number of species extensively across the world. The genetics of behavioural traits is understood to some extent, but it is seldom accounted for in breeding programmes. This review summarizes the estimates of genetic parameters for behavioural traits in cattle, pigs, poultry and fish. On the basis of the major studies performed in the last two decades, we focus the review on traits of common interest in the four species. These concern the behavioural responses to both acute and chronic stressors in the physical environment (feed, temperature, etc.) and those in the social environment (other group members, progeny, humans). The genetic strategies used to improve the behavioural capacity for adaptation of animals differ between species. There is a greater emphasis on responses to acute environmental stress in fish and birds, and on responses to chronic social stress in mammals.

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

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.

The mRNA expression of cortisol axis related genes differs in Atlantic cod (Gadus morhua) categorized as high or low responders

Cortisol is a major stress hormone in fish and is known, under normal or stressful conditions, to affect several physiological processes including growth and immunity. Thus, efforts have been made for several cultured finfish species, including the Atlantic cod, to determine whether fish with a high or low cortisol response to stress can be identified and selected. However, we have a limited understanding of the mechanisms that determine these two phenotypes. Thus, we measured total and free plasma cortisol levels in high and low responding cod when subjected to a 30 s handling stress, and the mRNA expression of four key genes in the glucocorticoid (i.e. cortisol) stress axis both pre- and post-stress. The cortisol data is consistent with our previous findings for cod, with high responding (HR) fish having ∼3-fold higher total and free plasma cortisol levels when compared to low responding (LR) fish. Three of the transcripts studied encode key proteins involved in steroidogenesis (StAR, P450scc and 3βHSD), and the constitutive mRNA expression of all three genes was significantly higher (∼2-fold) in the head kidney of HR fish when compared to LR cod. The other gene of interest was the glucocorticoid receptor (GR). We partly cloned and characterized a cDNA from Atlantic cod likely to be this fish's ortholog of the teleost GR1, and showed that while there was no difference in hepatic constitutive GR mRNA expression between groups, HR fish had liver GR mRNA levels that were significantly (1.8-fold) higher at 3 h post-stress as compared to LR fish. Our results suggest that the different magnitude of cortisol response between LR and HR fish is at least partially determined by the capacity of the interrenal tissue to produce steroids.

Detection of QTL with effects on osmoregulation capacities in the rainbow trout (Oncorhynchus mykiss)

Background

There is increasing evidence that the ability to adapt to seawater in teleost fish is modulated by genetic factors. Most studies have involved the comparison of species or strains and little is known about the genetic architecture of the trait. To address this question, we searched for QTL affecting osmoregulation capacities after transfer to saline water in a nonmigratory captive-bred population of rainbow trout.

Results

A QTL design (5 full-sib families, about 200 F2 progeny each) was produced from a cross between F0 grand-parents previously selected during two generations for a high or a low cortisol response after a standardized confinement stress. When fish were about 18 months old (near 204 g body weight), individual progeny were submitted to two successive hyper-osmotic challenges (30 ppt salinity) 14 days apart. Plasma chloride and sodium concentrations were recorded 24 h after each transfer. After the second challenge, fish were sacrificed and a gill index (weight of total gill arches corrected for body weight) was recorded. The genome scan was performed with 196 microsatellites and 85 SNP markers. Unitrait and multiple-trait QTL analyses were carried out on the whole dataset (5 families) through interval mapping methods with the QTLMap software. For post-challenge plasma ion concentrations, significant QTL (P < 0.05) were found on six different linkage groups and highly suggestive ones (P < 0.10) on two additional linkage groups. Most QTL affected concentrations of both chloride and sodium during both challenges, but some were specific to either chloride (2 QTL) or sodium (1 QTL) concentrations. Six QTL (4 significant, 2 suggestive) affecting gill index were discovered. Two were specific to the trait, while the others were also identified as QTL for post-challenge ion concentrations. Altogether, allelic effects were consistent for QTL affecting chloride and sodium concentrations but inconsistent for QTL affecting ion concentrations and gill morphology. There was no systematic lineage effect (grand-parental origin of QTL alleles) on the recorded traits.

Conclusions

For the first time, genomic loci associated with effects on major physiological components of osmotic adaptation to seawater in a nonmigratory fish were revealed. The results pave the way for further deciphering of the complex regulatory mechanisms underlying seawater adaptation and genes involved in osmoregulatory physiology in rainbow trout and other euryhaline fishes.

Identification of genes associated with heat tolerance in Arctic charr exposed to acute thermal stress

Arctic charr is an especially attractive aquaculture species given that it features the desirable tissue traits of other salmonids and is bred and grown at inland freshwater tank farms year round. It is of interest to develop upper temperature tolerant (UTT) strains of Arctic charr to increase the robustness of the species in the face of climate change and to enable production in more southern regions. We used a genomics approach that takes advantage of the well-studied Atlantic salmon genome to identify genes that are associated with UTT in Arctic charr. Specifically, we conducted an acute temperature trial to identify temperature tolerant and intolerant Arctic charr individuals, which were subject to microarray and qPCR analysis to identify candidate UTT genes. These were compared with genes annotated in a quantitative trait locus (QTL) region that was previously identified as associated with UTT in rainbow trout and Arctic charr and that we sequenced in Atlantic salmon. Our results suggest that small heat shock proteins as well as HSP-90 genes are associated with UTT. Furthermore, hemoglobin expression was significantly downregulated in tolerant compared with intolerant fish. Finally, QTL analysis and expression profiling identified COUP-TFII as a candidate UTT gene, although its specific role is unclear given the identification of two transcripts, which appear to have different expression patterns. Our results highlight the importance of using more than one approach to identify candidate genes, particularly when examining a complicated trait such as UTT in a highly complex genome for which there is no reference genome.

Estimates of linkage disequilibrium and effective population size in rainbow trout

Background

The use of molecular genetic technologies for broodstock management and selective breeding of aquaculture species is becoming increasingly more common with the continued development of genome tools and reagents. Several laboratories have produced genetic maps for rainbow trout to aid in the identification of loci affecting phenotypes of interest. These maps have resulted in the identification of many quantitative/qualitative trait loci affecting phenotypic variation in traits associated with albinism, disease resistance, temperature tolerance, sex determination, embryonic development rate, spawning date, condition factor and growth. Unfortunately, the elucidation of the precise allelic variation and/or genes underlying phenotypic diversity has yet to be achieved in this species having low marker densities and lacking a whole genome reference sequence. Experimental designs which integrate segregation analyses with linkage disequilibrium (LD) approaches facilitate the discovery of genes affecting important traits. To date the extent of LD has been characterized for humans and several agriculturally important livestock species but not for rainbow trout.

Results

We observed that the level of LD between syntenic loci decayed rapidly at distances greater than 2 cM which is similar to observations of LD in other agriculturally important species including cattle, sheep, pigs and chickens. However, in some cases significant LD was also observed up to 50 cM. Our estimate of effective population size based on genome wide estimates of LD for the NCCCWA broodstock population was 145, indicating that this population will respond well to high selection intensity. However, the range of effective population size based on individual chromosomes was 75.51 - 203.35, possibly indicating that suites of genes on each chromosome are disproportionately under selection pressures.

Conclusions

Our results indicate that large numbers of markers, more than are currently available for this species, will be required to enable the use of genome-wide integrated mapping approaches aimed at identifying genes of interest in rainbow trout.