- Open Access
Heritability of cortisol response to confinement stress in European sea bass dicentrarchus labrax
Genetics Selection Evolution volume 44, Article number: 15 (2012)
In fish, the most studied production traits in terms of heritability are body weight or growth, stress or disease resistance, while heritability of cortisol levels, widely used as a measure of response to stress, is less studied. In this study, we have estimated heritabilities of two growth traits (body weight and length) and of cortisol response to confinement stress in the European sea bass.
The F1 progeny analysed (n = 922) belonged to a small effective breeding population with contributions from an unbalanced family structure of just 10 males and 2 females. Heritability values ranged from 0.54 (±0.21) for body weight to 0.65 (±0.22) for standard body length and were low for cortisol response i.e. 0.08 (±0.06). Genetic correlations were positive (0.94) between standard body length and body weight and negative between cortisol and body weight and between cortisol and standard body length (−0.60 and −0.55, respectively).
This study confirms that in European sea bass, heritability of growth-related traits is high and that selection on such traits has potential. However, heritability of cortisol response to stress is low in European sea bass and since it is known to vary greatly among species, further studies are necessary to understand the reasons for these differences.
Farming of European sea bass (Dicentrarchus labrax, Moronidae, Teleostei), represents about 100 000 tons produced per year  and attracts extensive interest as a major fish species for establishing breeding programmes to improve production traits. In fish, the most studied production traits in terms of heritability are body weight or growth, stress or disease resistance [2–6].
In this work, our aim was to set up European sea bass families by assigning parentage and heritability for three traits i.e. cortisol response to stress, body weight and standard body length to the progeny derived from the batch of a single spawning day.
The methodology used for producing, phenotyping and genotyping the F1 population has been described by Massault et al . In summary, the broodstock consisted of 34 females, 22 males and one individual of undetermined sex originating from wild and caged fish. From this broodstock, 2000 offspring were raised for 254 days under standard farm conditions and then distributed into four tanks of 45 m3, each with a net covering the inner surface. After this period of acclimatization, a confinement stress was applied, which consisted in slowly pulling the inner net of each tank so that the fish were confined in a volume of approximately 0.2 m3. After 4 h of confinement, the net was lifted and emptied into a tank of icy water, a process which stuns the fish within 3 min. Each group of 500 fish was numbered serially, bled, weighed and digitally photographed within 140 min after stunning, either in the morning (11 am) or afternoon (3 pm). Blood plasma was stored at −20°C for cortisol analysis and red blood cells stored in absolute ethanol for genetic analysis. Cortisol was measured by radioimmunoassay and microsatellite genotyping was carried out with three multiplex PCR (polymerase chain reaction) [See Additional file 1. All pedigree genotypes from the 11 larger families (n = 922) were checked for Mendelian errors before estimating heritabilities and correlations. Parentage assignment was implemented using three software packages [i.e., CERVUS v.3.0; , PAPA v.2.0; , VITASSIGN v.1.0; ] in order to constitute families with the highest possible certainty. The genotyping error rate was set to 1%. The assignment was tested for power and performance and locus-specific polymorphism information content (PIC) values were calculated.
Heritabilities and phenotypic correlations were calculated using phenotypic data collected on 930 animals. Eight animals were removed because of missing phenotypes. Thus, the dataset used to estimate heritability values comprised 922 animals, with missing values in some variables (see Table 1). Heritabilities were estimated using ASReml fitting an animal model. Several fixed effects were tested (sample set, day, tube number and assay number in cortisol assays) to check if they influenced the trait in question. With the exception of cortisol, the model was
TRAIT represents the phenotypic trait, μ the trait mean, A the additive genetic effect and E the environmental effect. Only sample set was found to have an effect on cortisol and therefore sample set was added as fixed effect to the model.
Phenotypic correlations were calculated with the software genstat v.10.
Parentage assignment and contributions
The 1151 progeny and 56 parents were genotyped at 29 microsatellite loci. The number of alleles per locus varied between two and 10 and PIC values varied between 0.124 and 0.767 [see Additional file 2]. A power analysis was conducted with the rarefaction method and showed that 10 loci were sufficient for a reliable assignment (details not shown). A core group of five families contributed 71.5% to the progeny, six families made a measurable contribution and 80 families only a very small contribution. This is a highly skewed family representation of at least 748 dam x sire combinations with a low effective breeding size, which might affect the estimates through the unwanted genetic drift (or limited Mendelian sampling) caused by the skewed representation. For the heritability and correlation analyses, 922 offspring were used from which two females (5.9%) and 10 males (45.5%) contributed the most [see Additional file 3]. Our study shows that within a single breeding day, the majority of the progeny can be produced with the contribution of just two females and 10 males, which amounts to an effective population size (Ne) of 6.7 while there were 43 participating breeders and a total number of initial individuals of 57 (Ne/N = 0.12). Thus, artificial insemination provides the best guarantee to set up experimental crosses since during natural spawning the number of families with significant contribution can be small.
Basic descriptive statistics for the phenotypes are shown in Table 1. Mean cortisol levels were constant over the time of blood collection as indicated by the horizontal regression lines in Figure 1 (range of coefficients of linear regression per tank −0.153 to 0.109; r2 = 2.10-3).
Heritabilities and correlations
Heritability values for growth traits ranged from 0.54 ± 0.21 for BW to 0.65 ± 0.22 for SL (Table 2). Such heritability estimates support the large proportion of phenotypic variation explained by the QTL detected in Massault et al.  and the reasons are discussed therein. For response to stress, the heritability of CORT was 0.08 ± 0.06. However, as noted above, the family structure is clearly suboptimal to estimate heritabilities as evidenced by the high standard errors of the estimates . As expected, the phenotypic correlation between BW and SL was high (0.94), whereas CORT was not phenotypically correlated to either growth trait. The genetic correlation between SL and BW was high (0.94) and that between CORT and SL or BW was negative (−0.55 and −0.60, respectively). These estimates confirm data from previous studies in European seabass , rainbow trout  and Nile tilapia . The value for body weight heritability obtained in our study also agrees with estimates previously reported in European sea bass, which range from medium heritabilities (0.2 in , 0.38-0.44 in  and 0.39 in [16, 17]) to high heritabilities when taking into account different environments (0.31-0.60 in ). Generally, body length and body weight have moderate to high heritability values in teleost fishes: 0.6 in Coho salmon, Oncorhynchus kisutch, 0.12-0.47 in brown trout, Salmo trutta, 0.09-0.44 in carp, Cyprinus carpio[21, 22], 0.38-0.79 in Nile tilapia, Oreochromis niloticus[23, 24], 0.64 (± 0.12) in cod, Gadus morhua, 0.38 ± 0.07 in gilthead seabream, Sparus aurata.
To our knowledge, studies on the heritability of cortisol response to stress in fish have been limited to salmonids [27–29] and cyprinids , which limits generalizations. In addition, lines with high and low cortisol response to stress have been selected in rainbow trout . However, the heritability of cortisol response to stress appears to be variable even among related species: 0.27-0.50 in rainbow trout, Oncorhynchus mykiss[27, 29, 30] and 0.60 in carp , but only 0.05 in Atlantic salmon Salmo salar[27, 31]. These discrepancies can be partly explained by the differences in species and methodologies used to determine the cortisol response. It should be noted that the time-dependence of cortisol response to stress is a key element and a potential source of error. However, the methodology used in our study seems reliable since no apparent drift in cortisol levels with time was observed after applying the confinement stress (Figure 1). A comparative study on stunning methods used in different fish farms for European sea bass reported mean levels of cortisol response similar to that obtained here, corresponding to a 5-fold increase in cortisol compared to resting values when using ice . In a pilot study on gilthead seabream using the same methodology, a significant change in cortisol levels of a control group not subjected to stress was observed, but there was no significant additional effect of ice-water on a group subjected to confinement stress (A. Canario, unpublished observations). Thus, the methodology adopted here for European sea bass is appropriate and even 2 h after stunning, the levels of cortisol obtained are due to the response to confinement stress and should reflect individual variation. In conclusion, in European sea bass, the growth traits measured have a moderate to high heritability but the cortisol level, as an indicator of response to stress, has a low heritability. Whether this low heritability derives from an artefact or an unbalanced family structure or whether it has a true biological base needs further clarification.
FAO: The state of the world fisheries and aquaculture. 2010, Food and Agriculture Organization of the United Nations, Rome, 197. NationsRome197FAO: The state of the world fisheries and aquaculture. Rome: Food and Agriculture Organization of the United Nations; 2010:197.
Henryon M, Jokumsen A, Berg P, Lund I, Pedersen PB, Olesen NJ, Slierendrecht WJ: Genetic variation for growth rate, feed conversion efficiency, and disease resistance exists within a farmed population of rainbow trout. Aquaculture. 2002, 209: 59-76. 10.1016/S0044-8486(01)00729-3.
Fevolden SE, Roed KH, Fjalestad KT, Stien J: Poststress levels of lysozyme and cortisol in adult rainbow trout: heritabilities and genetic correlations. J Fish Biol. 1999, 54: 900-910. 10.1111/j.1095-8649.1999.tb02040.x.
del Valle G, Taniguchi N, Tsujimura A: Genetic differences in some haematological traits of communally reared clonal ayu, Plecoglossus altivelis Temminck & Schlegel, under stressed and non-stressed conditions. Aquaculture Res. 1996, 27: 787-793.
Beacham TD, Evelyn TPT: Genetic variation in disease resistance and growth of chinook, coho, and chum salmon with respect to vibriosis, furunculosis, and bacterial kidney disease. Trans Am Fish Soc. 1992, 121: 456-485. 10.1577/1548-8659(1992)121<0456:GVIDRA>2.3.CO;2.
Gjedrem T, Salte R, Gjoen HM: Genetic variation in susceptibility of Atlantic salmon to furunculosis. Aquaculture. 1991, 97: 1-6. 10.1016/0044-8486(91)90274-B.
Massault C, Hellemans B, Louro B, Batargias C, Van Houdt JK, Canario A, Volckaert FA, Bovenhuis H, Haley C, de Koning DJ: QTL for body weight, morphometric traits and stress response in European sea bass Dicentrarchus labrax. Anim Genet. 2010, 41: 337-345.
Marshall TC, Slate J, Kruuk LEB, Pemberton JM: Statistical confidence for likelihood-based paternity inference in natural populations. Mol Ecol. 1998, 7: 639-655. 10.1046/j.1365-294x.1998.00374.x.
Duchesne P, Godbout M-H, Bernatchez : PAPA (package for the analysis of parental allocation): a computer program for simulated and real parental allocation. Mol Ecol Notes. 2002, 2: 191-193. 10.1046/j.1471-8286.2002.00164.x.
Vandeputte M, Mauger S, Dupont-Nivet M: An evaluation of allowing for mismatches as a way to manage genotyping errors in parentage assignment by exclusion. Mol Ecol Notes. 2006, 6: 265-267. 10.1111/j.1471-8286.2005.01167.x.
Dupont-Nivet M, Vandeputte M, Chevassus B: Optimization of factorial mating designs for inference on heritability in fish species. Aquaculture. 2002, 204: 361-370. 10.1016/S0044-8486(01)00839-0.
Chatziplis D, Batargias C, Tsigenopoulos CS, Magoulas A, Kollias S, Kotoulas G, Volckaert FAM, Haley CS: Mapping quantitative trait loci in European sea bass (Dicentrarchus labrax): The BASSMAP pilot study. Aquaculture. 2007, 272: S172-S182.
Gjerde B, Schaeffer LR: Body traits in rainbow trout. 2. Estimates of heritabilities and of phenotypic and genetic correlations. Aquaculture. 1989, 80: 25-44. 10.1016/0044-8486(89)90271-8.
Rutten MJM, Bovenhuis H, Komen H: Genetic parameters for fillet traits and body measurements in Nile tilapia (Oreochromis niloticus L.). Aquaculture. 2005, 246: 125-132. 10.1016/j.aquaculture.2005.01.006.
Dupont-Nivet M, Vandeputte M, Vergnet A, Merdy O, Haffray P, Chavanne H, Chatain B: Heritabilities and GxE interactions for growth in the European sea bass (Dicentrarchus labrax L.) using a marker-based pedigree. Aquaculture. 2008, 275: 81-87. 10.1016/j.aquaculture.2007.12.032.
Grima L, Chatain B, Ruelle F, Vergnet A, Launay A, Mambrini M, Vandeputte M: In search for indirect criteria to improve feed utilization efficiency in sea bass (Dicentrarchus labrax): Part II: Heritability of weight loss during feed deprivation and weight gain during re-feeding periods. Aquaculture. 2010, 302: 169-174. 10.1016/j.aquaculture.2010.02.016.
Costa C, Vandeputte M, Antonucci F, Boglione C, Menesatti P, Cenadelli S, Parati K, Chavanne H, Chatain B: Genetic and environmental influences on shape variation in the European sea bass (Dicentrarchus labrax). Biol J Linn Soc. 2010, 101: 427-436. 10.1111/j.1095-8312.2010.01512.x.
Saillant E, Dupont-Nivet M, Haffray P, Chatain B: Estimates of heritability and genotype-environment interactions for body weight in sea bass (Dicentrarchus labrax L.) raised under communal rearing conditions. Aquaculture. 2006, 254: 139-147. 10.1016/j.aquaculture.2005.10.018.
Withler RE, Beacham TD: Genetic variation in body weight and flesh color of the coho salmon (Oncorhynchus kisutch) in British Columbia. Aquaculture. 1994, 119: 135-148. 10.1016/0044-8486(94)90170-8.
Blanc JM, Poisson H: Genetic variation of body size, condition and pyloric caeca number in juvenile brown trout, Salmo trutta L. Aquaculture Res. 2006, 37: 637-642. 10.1111/j.1365-2109.2006.01474.x.
Tanck MWT, Vermeulen KJ, Bovenhuis H, Komen H: Heredity of stress-related cortisol response in androgenetic common carp (Cyprinus carpio L.). Aquaculture. 2001, 199: 283-294. 10.1016/S0044-8486(01)00518-X.
Vandeputte M, Kocour M, Mauger S, Rodina M, Launay A, Gela D, Dupont-Nivet M, Hulak M, Linhart O: Genetic variation for growth at one and two summers of age in the common carp (Cyprinus carpio L.): Heritability estimates and response to selection. Aquaculture. 2008, 277: 7-13. 10.1016/j.aquaculture.2008.02.009.
Charo-Karisa H, Bovenhuis H, Rezk MA, Ponzoni RW, van Arendonk JAM, Komen H: Phenotypic and genetic parameters for body measurements, reproductive traits and gut length of Nile tilapia (Oreochromis niloticus) selected for growth in low-input earthen ponds. Aquaculture. 2007, 273: 15-23. 10.1016/j.aquaculture.2007.09.011.
Charo-Karisa H, Komen H, Rezk MA, Ponzoni RW, van Arendonk JAM, Bovenhuis H: Heritability estimates and response to selection for growth of Nile tilapia (Oreochromis niloticus) in low-input earthen ponds. Aquaculture. 2006, 261: 479-486. 10.1016/j.aquaculture.2006.07.007.
Kolstad K, Thorland I, Refstie T, Gjerde B: Genetic variation and genotype by location interaction in body weight, spinal deformity and sexual maturity in Atlantic cod (Gadus morhua) reared at different locations off Norway. Aquaculture. 2006, 259: 66-73. 10.1016/j.aquaculture.2005.12.022.
Antonello J, Massault C, Franch R, Haley C, Pellizzari C, Bovo G, Patarnello T, de Koning DJ, Bargelloni L: Estimates of heritability and genetic correlation for body length and resistance to fish pasteurellosis in the gilthead sea bream (Sparus aurata L.). Aquaculture. 2009, 298: 29-35. 10.1016/j.aquaculture.2009.10.022.
Fevolden SE, Refstie T, Gjerde B: Genetic and phenotypic parameters for cortisol and glucose stress-response in Atlantic salmon and rainbow trout. Aquaculture. 1993, 118: 205-216. 10.1016/0044-8486(93)90457-A.
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. N Am J Aqualcult. 2008, 70: 425-433. 10.1577/A07-059.1.
Pottinger TG, Carrick TR: Modification of the plasma cortisol response to stress in rainbow trout by selective breeding. Gen Comp Endocrinol. 1999, 116: 122-132. 10.1006/gcen.1999.7355.
Fevolden SE, Roed KH, Fjalestad KT: Selection response of cortisol and lysozyme in rainbow trout and correlation to growth. Aquaculture. 2002, 205: 61-75. 10.1016/S0044-8486(01)00660-3.
Fevolden SE, Roed KH, Gjerde B: Genetic components of post-stress cortisol and lysozyme activity in Atlantic salmon: Correlations to disease resistance. Fish Shellfish Immun. 1994, 4: 507-519. 10.1006/fsim.1994.1045.
Acerete L, Reig L, Alvarez D, Flos R, Tort L: Comparison of two stunning/slaughtering methods on stress response and quality indicators of European sea bass (Dicentrarchus labrax). Aquaculture. 2009, 287: 139-144. 10.1016/j.aquaculture.2008.10.012.
This study was funded by the European Commission (project AQUAFIRST contract number FP6-STREP-2004-513692). DJK and CSH acknowledge financial support from the BBSRC. The authors are indebted to E Couto, S Darivianaku, J Fuentes, P Guerreiro, M Kampakli, G Kotoulas and DM Power for their invaluable help with sample collection.
The authors declare no competing interests.
FAMV planned the genotyping and wrote the manuscript, BH performed the genotyping, CB participated in the experimental planning, organized and carried out the experiment, participated in the sampling and wrote the paper, BL participated in the sampling and genotyping, CM performed the heritability and correlation analyses, JKJVH participated in the experimental planning and genetic analysis, CH participated in the experimental planning and genetic analysis, D-JK participated in the experimental planning, genetic analysis and wrote the manuscript, AVMC participated in the experimental planning, participated in the sampling, performed the cortisol analyses and wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Loci included in each multiplex set and assigned to D. labrax linkage groups, primer concentrations and comments on the 98 microsatellite markers used to scan the genome of European sea bass. (DOCX 15 KB)
Additional file 2: Parentage assignment statistics. Locus-specific PIC (Polymorphic Information Content) values, test of Hardy-Weinberg equilibrium (**: P < 0.05), null alleles, cumulative parentage assignment (in percent) with one and no parent known at high (95%) and low (80%) stringency. (DOCX 16 KB)
Additional file 3: resspecies ID code of male and female parent, and number of offspring of the 11 largest families (FS01 to FS11) of the experimental population of European sea bass. The pedigree is based on assignments obtained with the software packages cervus, papa, and vitassign, and submission to the resspecies database (n = 922). (DOCX 14 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Volckaert, F.A., Hellemans, B., Batargias, C. et al. Heritability of cortisol response to confinement stress in European sea bass dicentrarchus labrax. Genet Sel Evol 44, 15 (2012). https://doi.org/10.1186/1297-9686-44-15
- Rainbow Trout
- Cortisol Level
- Polymorphism Information Content
- Growth Trait