Association of HPA axis-related genetic variation with stress reactivity and aggressive behaviour in pigs
BMC Genetics volume 11, Article number: 74 (2010)
Stress, elicited for example by aggressive interactions, has negative effects on various biological functions including immune defence, reproduction, growth, and, in livestock, on product quality. Stress response and aggressiveness are mutually interrelated and show large interindividual variation, partly attributable to genetic factors. In the pig little is known about the molecular-genetic background of the variation in stress responsiveness and aggressiveness. To identify candidate genes we analyzed association of DNA markers in each of ten genes (CRH g.233C>T, CRHR1 c.*866_867insA, CRHBP c.51G>A, POMC c.293_298del, MC2R c.306T>G, NR3C1 c.*2122A>G, AVP c.207A>G, AVPR1B c.1084A>G, UCN g.1329T>C, CRHR2 c.*13T>C) related to the hypothalamic-pituitary-adrenocortical (HPA) axis, one of the main stress-response systems, with various stress- and aggression-related parameters at slaughter. These parameters were: physiological measures of the stress response (plasma concentrations of cortisol, creatine kinase, glucose, and lactate), adrenal weight (which is a parameter reflecting activity of the central branch of the HPA axis over time) and aggressive behaviour (measured by means of lesion scoring) in the context of psychosocial stress of mixing individuals with different aggressive temperament.
The SNP NR3C1 c.*2122A>G showed association with cortisol concentration (p = 0.024), adrenal weight (p = 0.003) and aggressive behaviour (front lesion score, p = 0.012; total lesion score p = 0.045). The SNP AVPR1B c.1084A>G showed a highly significant association with aggressive behaviour (middle lesion score, p = 0.007; total lesion score p = 0.003). The SNP UCN g.1329T>C showed association with adrenal weight (p = 0.019) and aggressive behaviour (front lesion score, p = 0.029). The SNP CRH g.233C>T showed a significant association with glucose concentration (p = 0.002), and the polymorphisms POMC c.293_298del and MC2R c.306T>G with adrenal weight (p = 0.027 and p < 0.0001 respectively).
The multiple and consistent associations shown by SNP in NR3C1 and AVPR1B provide convincing evidence for genuine effects of their DNA sequence variation on stress responsiveness and aggressive behaviour. Identification of the causal functional molecular polymorphisms would not only provide markers useful for pig breeding but also insight into the molecular bases of the stress response and aggressive behaviour in general.
Stress responses promote the maintenance of homeostasis and adaptation to physiological and psychosocial challenges of a changing environment. This complex process involves coordinated activation of behavioural, autonomic, and neuroendocrine reactions. Concomitantly, pathways that promote vegetative functions such as growth, reproduction, and feeding are inhibited, the extent to which is dependent upon the duration and intensity of the stressor ("biological cost" of the stress response ).
Aggression is a powerful stressor and has been shown to activate the hypothalamic-pituitary-adrenocortical (HPA) axis as well as the sympatho-adrenomedullar (SAM) system in various species including pigs [2–4]. In the pig, aggression commonly occurs when mixing unfamiliar individuals which disturbs the social dominance order . Besides negative effects on animal welfare, aggression has also been shown to have a negative impact on immune response , growth performance , and product quality in pigs [4, 8]. Aggressive behaviour in turn is affected by the functional properties of the HPA axis and of the SAM system (reviewed in ). Baseline levels of glucocorticoids are inversely related with aggressiveness in various species including pigs (reviewed in [9, 6, 10]). In contrast, as reported in rodents, an acute increase in glucocorticoid levels may promote aggressive behaviour by a fast feed forward mechanism .
Functional properties of the stress response systems and aggressive behaviour show large interindividual variation, depending upon a variety of factors including genetic predisposition [12–14]. In humans and in model animals several genomic regions and gene variants associated with variation in stress responsiveness and aggressiveness have been identified using quantitative trait loci (QTL) mapping and candidate gene approaches (reviewed in [15–18]). In the pig, such studies are scarce, in spite of the widely recognized impact of aggressive behaviour and stress on the expression of meat and carcass quality traits . Désautés et al. mapped QTL for behavioural and neuroendocrine stress responses in a Meishan × Large White intercross leading to the identification of the corticosteroid binding globulin encoding gene as a major QTL for plasma cortisol levels in this population (, reviewed in ). In addition Geldermann et al. mapped QTL for plasma creatine kinase levels after pharmacological challenge. Using the candidate gene approach Fujii et al. identified a mutation (SNP c.1843C>T) in the ryanodine receptor 1 (RYR1) gene responsible for malignant hyperthermia and porcine stress syndrome. The mutation has been shown to affect the basal functioning of the HPA axis in vivo and in vitro, and analysis of this important mutation was included in the present study.
The aim of the present study was to expand the current knowledge of the molecular-genetic basis of the variation in stress responsiveness and aggressive behaviour in the pig. To this end we analyzed the association of candidate gene DNA markers with physiological parameters of the stress response (cortisol, creatine kinase, glucose, and lactate concentration in plasma) and aggressive behaviour (measured by means of lesion scoring) in the context of the psychosocial stress of mixing individuals with different aggressive temperaments in a commercial pig herd . In addition, we also analyzed the association with adrenal weight, a parameter reflecting activity of the central branch of the HPA axis (i.e. the release of corticotropin-releasing hormone and ACTH) over time, in a different herd of commercial crossbred pigs. To relate position of the candidate genes with known porcine QTL we physically mapped those candidates whose position had not previously been determined.
The candidates represented genes encoding members of the HPA axis pathway (corticotropin-releasing hormone, CRH; CRH type 1 receptor, CRHR1; CRH binding protein, CRHBP; vasopressin, AVP; vasopressin V1B receptor, AVPR1B; proopiomelanocortin, POMC; melanocortin type 2 (ACTH) receptor, MC2R; glucocorticoid receptor, NR3C1) and genes encoding members of the corticotropin-releasing hormone (CRH) system (urocortin, UCN; CRH type 2 receptor, CRHR2). CRH and vasopressin act synergistically to activate expression of proopiomelanocortin, the precursor of adrenocorticotropin (ACTH), in the pituitary through activation of CRH type 1 and vasopressin V1B receptors respectively. CRH is the dominant trigger for HPA axis activation during acute stress while vasopressin, which itself is a weak ACTH secretagogue, is important in mediating response to chronic stress. Vasopressin also has an important function in the control of aggressive behaviour in various species  including the pig . CRH binding protein functions as a buffer for CRH and related peptides and plays an inhibitory role in the modulation of CRH activity. ACTH stimulates synthesis and secretion of glucocorticoids from the adrenal cortex via melanocortin type 2 receptor. The action of glucocorticoids on target tissues is mediated by glucocorticoid receptor, which in the hypothalamus, in the pituitary, and in the hippocampus acts to terminate the stress response (reviewed in ). The urocortins, although not directly involved in the HPA axis action, may play a modulatory role via stimulation of the CRH receptor type 2 .
Results and Discussion
Allele distribution of the candidate genes
The candidate gene DNA markers are located mainly in the transcribed region (Table 1) with the exception of the SNPs CRH g.233C>T and UCN g.1329T>C, where the DNA markers are located in the promoter region in evolutionarily conserved segments, ~350bp and ~90bp upstream of the transcription start site respectively ( and unpublished data respectively). Two DNA markers lead to amino acid exchanges (CRHBP c.51G>A, AVPR1B c.1084A>G), and one to an amino acid deletion (POMC c.293_298del), respectively. The remaining markers are either silent SNPs (MC2R c.306T>G, AVP c.207A>G) or polymorphisms located in the 3' untranslated region (CRHR1 c.*866_867insA, CRHR2 c.*13T>C, NR3C1 c.*2122A>G). The SNP c.306T>G in MC2R was reported previously by others .
Allele frequencies of the candidate gene DNA markers in the two phenotyped commercial crossbred herds (SYN and PiF1a), along with allele frequencies in sets of pigs from three commercial pure breeds German Landrace (LR), German Large White (LW), and Pietrain (Pi), in a set from Vietnamese local breed Muong Khuong (MK), and in a set of European Wild Boars (WB) are summarized in Table 1. Contrary to the common assumption that selective breeding reduces genetic variability, the commercial breeds and crosses showed the highest variability, whereas Wild Boar the least. Six out of the ten candidate gene DNA markers tested, namely CRHR1 c.*866_867insA, CRHBP c.51G>A, AVPR1B c.1084A>G, POMC c.293_298del, MC2R c.306T>G, and UCN g.1329T>C were fixed in Wild Boar. The increased variability of the commercial breeds compared to Wild Boar may be a result of the introgression of Asian genetics into European breeds in 18th and 19th centuries. Ramirez et al. showed that the proportion of Asian genetics in the genetic pool of European commercial breeds is still significant at around 47-61%. The majority of the candidate gene DNA markers, as for example the SNP AVPR1B c.1084A>G, showed large differences in allele distribution between Wild Boar and Muong Khuong, indicating that the alleles missing in Wild Boar but segregating in commercial breeds may originate from the introgression of the Asian genetics.
In view of the intriguing differences in allele frequency of the DNA markers of the HPA axis-related genes between commercial, local, and wild pigs it is interesting to note that domestication of several species and selection for tameness (i.e. reduced interspecific aggression) in silver foxes and rats, which is an important aspect of domestication, was coupled with decreased activity of the HPA axis (, reviewed in ). Likewise, domestic pig breeds show lower cortisol levels compared to their wild ancestor . Although the reduced cortisol level in modern pigs may result from selection for leaner body composition , it is conceivable that it might be related to domestication-related behaviour as well.
A test of Hardy-Weinberg equilibrium revealed significant deviation in both crossbred herds for the SNPs RYR1 c.1843C>T and AVPR1B c.1084A>G, for the SNPs AVP c.207A>G and UCN g.1329T>C in SYN, and for the SNP NR3C1 c.*2122A>G in PiF1a. In contrast in the pure breeds, only the polymorphisms CRHR1 c.*866_867insA in Pietrain and CRHR2 c.*13T>C in Muong Khuong showed significant deviation, i.e. in two tests out of fifty performed, which is a proportion that would be expected to occur at the significance level of 5% simply by chance. The more frequent deviation from Hardy-Weinberg equilibrium in SYN and PiF1a is thus likely a result of crossbreeding, rather than an indication of selection.
Association with stress responsiveness and aggressive behaviour
The physiological parameters of the stress response included in this study were shown previously to be affected by aggressive interactions [2, 36] and also in the present experiment as reported earlier . In this context creatine kinase (CK), lactate and glucose concentration reflect the amount of physical activity and mobilization of energy sources, and indicate activation of the SAM system, whilst cortisol level indicates activation of the HPA axis. For evaluating aggressive behaviour at mixing we used the number of skin lesions, which is an established indicator trait. Turner et al.[14, 37, 38] showed that front lesion score is the best indicator of engagement in reciprocal aggression at mixing, whereas rear lesion score is the best indicator of receipt of nonreciprocal aggression at mixing.
Effects with p ≤ 0.10 are summarized in Tables 2 and 3 for stress parameters and in Table 4 for lesion scores. We refer to results as significant when p < 0.05 and as showing a tendency when 0.05 ≤ p ≤ 0.10. To advise the reader of the increased risk of type I error due to multiple testing we also provide the corresponding false discovery rate (q-value; Tables 2, 3, 4). The q value is similar to the well known p value, except it is a measure of significance in terms of the false discovery rate rather than the false positive rate .
All physiological stress parameters were significantly affected by the RYR1 SNP c.1843C>T (Table 2). In accordance with the known effect of the RYR1 SNP on metabolism of the skeletal muscle [40, 41] the mutated T allele highly significantly (p < 0.0001) increased plasma CK and lactate concentration. Plasma glucose concentration, in turn, was reduced. Consistent with Weaver et al., heterozygous carriers of the T allele showed significantly lower cortisol concentration compared to CC homozygous individuals. Furthermore the adrenal weight tended to be lower in the heterozygous carriers of the T allele in the PiF1a line, providing independent supporting evidence that the RYR1 SNP affects activity of the HPA axis (Table 3). With increased sample size this effect became significant (n = 316; data not shown). Concerning aggressive behaviour the T allele of the RYR1 SNP c.1843C>T significantly increased rear lesion score (Table 4). The effect of RYR1 on aggressive behaviour might be related to its effect on HPA axis activity. Guárdia et al. also reported an effect of RYR1 on skin lesions, but the direction of the effects was opposite to our findings, perhaps due to differences in recording of the skin lesions (in the study of Guárdia et al. as a whole carcass score on a 5-point scale).
The SNP CRH g.233C>T showed a significant association with plasma glucose concentration and a tendency to affect plasma CK concentration; however the effect on glucose showed neither an additive nor dominance direction (Table 2). We mapped CRH previously on SSC4 in the marker interval SW724-S0107 . So far no QTL for physiological stress parameters have been reported in this genomic region. The effect on CK concentration is possibly related to aggressive behaviour, because carriers of the T allele, which tended to have higher CK concentration, also tended to have a higher rear lesion score (Table 4). It could be speculated that the association of the SNP CRH g.233C>T with aggressive behaviour might be related to the anxiogenic effect of CRH [43, 44] because the SNP showed no association with activity of the HPA axis. The anxiogenic effect of CRH is, at least partly, independent of its action on the HPA axis . Anxiety is related to aggressive behaviour in a complex manner, depending on the model used. Veenema and Neumann  for example, found an inverse relationship between anxiety and offensive aggression in rats divergently selected for anxiety-related behaviour.
CRHR1 maps on SSC12 in the marker interval SW957-SW943 (Additional File 1, Table S6). No effects, on either the stress parameters or on aggressive behaviour were detected for the polymorphism CRHR1 c.*866_867insA. CRHBP maps on SSC2 in the marker interval SW1602-SW1320 (Additional File 1, Table S6). The SNP CRHBP c.51G>A showed a tendency to affect middle lesion score (Table 4), but showed no association with physiological stress parameters or adrenal weight.
The polymorphisms POMC c.293_298del and MC2R c.306T>G both showed association only with adrenal weight in the PiF1a line (Table 3). This effect is in line with the established positive effect of POMC-derived peptides, in particular ACTH, on adrenal growth . We mapped POMC on SSC3, in the marker interval SW314-S0002 (Additional File 1, Table S6), close to a QTL for basal glucose level . POMC is involved in the control of glucose homeostasis via the HPA axis and other pathways , however in the present study we found no association of the polymorphism POMC c.293_298del with glucose concentration. MC2R was mapped on SSC6 close to marker SW2173 by Jacobs et al. and close to QTL for plasma creatine kinase levels in the Meishan × Pietrain and Wild Boar × Pietrain families . However these QTL are most likely caused by the SNP c.1843C>T in RYR1, which is also located in the same QTL region.
The SNP NR3C1 c.*2122A>G showed a significant association with plasma cortisol concentration in SYN and a significant association with adrenal weight in PiF1a (Tables 2 and 3). The allele A, which is associated with lower cortisol concentration, is also associated with lower adrenal weight (Tables 2 and 3), suggesting that it is associated with an enhanced negative feedback effect on the activity of the HPA axis. A similar phenotype, including decreased basal plasma corticosterone and ACTH levels, reduced adrenal weight and adrenocortical size, has been reported in a knock-in mouse line showing increased functioning of a modified glucocorticoid receptor . A clue about possible molecular mechanisms underlying the enhanced negative feedback effect might be obtained from the report of Perreau et al. who observed a higher density of glucocorticoid receptors in the pituitary in Large White pigs which exhibit lower activity of the HPA axis compared to Meishan pigs, suggesting that genetic factors cause variation in glucocorticoid receptor density in the pituitary of the pig. Furthermore, the SNP NR3C1 c.*2122A>G showed a significant association with the total and front lesion score, but not with the middle or rear lesion score (Table 4). This indicates that the animals accumulated lesions mainly through reciprocal fighting. Because, as mentioned in the introduction, baseline levels of glucocorticoids are inversely related with aggressiveness, the decreased activity of the HPA axis in the carriers of the allele A provides a possible explanation for their enhanced aggressive behaviour. NR3C1 maps on SSC2 in the marker interval SW1879-SWR308 (Additional File 1, Table S6). So far no QTL for physiological stress parameters have been reported in this genomic region.
AVP maps on SSC17 in a QTL region for post-stress ACTH level . However, in the present study the SNP AVP c.207A>G showed associations with neither the physiological stress parameters, nor with aggressive behaviour. The SNP AVPR1B c.1084A>G in turn showed significant association only with aggressive behaviour (Table 4). The allele G consistently decreased lesion score. However, the effect was most pronounced for the middle lesion score, less pronounced for the front lesion score, and did not reach the 0.1 level for the rear lesion score (Table 4). The vasopressin pathway plays a prominent role in the regulation of social behaviour, including aggression [25, 26]. Vasopressin V1B receptor knockout mice display reduced levels of social forms of aggressive behaviour. While male vasopressin V1B receptor knockout mice demonstrate deficits in offensive and defensive aggression, female knockout mice have deficits in maternal aggression. However, there is no global deficit in aggressive behaviour as vasopressin V1B receptor knockout mice show normal predatory aggression . Hence, there is strong functional evidence supporting the identified association of the SNP AVPR1B c.1084A>G with aggressive behaviour in the context of social stress of mixing. AVPR1B maps on SSC9 in the marker interval SW1879-SWR308 (Additional File 1, Table S6).
The SNP UCN g.1329T>C showed no association with physiological stress parameters but a significant association with adrenal weight in PiF1a (Table 3) indicating that it might affect activity of the HPA axis. Urocortin is thought to be involved in the autonomic stress response  but so far studies with knockout mouse models have revealed no consistent evidence for an involvement of urocortin in the HPA axis response to acute stress (reviewed in ). The study of Zalutskaya et al. on the response of urocortin knockout mice to repeated restraint indicates that urocortin may play a role in the adaptation of the HPA axis to chronic stress. Little is known about the function of urocortin in the pig. Parrott et al. showed that intracerebroventricular injection of urocortin increases cortisol release in the pig. Besides association with adrenal weight the SNP UCN g.1329T>C also showed significant association with front lesion score (Table 4). Urocortin, similar to CRH, possess an anxiogenic effect [43, 56]. As discussed above for the SNP CRH g.233C>T, this might also underlie the association of the SNP UCN g.1329T>C with aggressive behaviour. UCN maps on SSC3 in the marker interval SW730-SWR201 (Additional File 1, Table S6), close to the position of POMC and to a QTL for basal glucose level . Urocortins are involved in the regulation of glucose homeostasis [57, 58]. However in the present study we found no association of the SNP UCN g.1329T>C with the glucose concentration.
CRHR2 maps on SSC18 in the marker interval SW787-SW1682 (Additional File 1, Table S6). The only effect we found for the SNP CRHR2 c.*13T>C was a tendency for adrenal weight in PiF1a (Table 3). On SSC18 QTL for basal and stress-induced increase in cortisol level were detected , however these map distal to CRHR2.
In the present study to examine phenotypic effects of the HPA axis-related genes we used a single DNA marker per gene. However, the HPA axis-related genes are usually highly polymorphic with several polymorphisms affecting gene expression and/or function of the encoded protein [29, 59, 60], with the consequence that a single DNA marker most likely does not capture all of the genetic information. Furthermore, for several polymorphisms, e.g. for POMC c.293_298del or CRHBP c.51G>A, the power to detect an association was limited by the low frequency of the minor allele. Therefore, genes that showed no significant associations here could still harbour functional DNA sequence variation with phenotypic effects. On the other hand the SNP in NR3C1 and AVPR1B showed multiple consistent effects, partly significant even after correction for multiple testing (the SNP NR3C1 c.*2122A>G on adrenal weight and the SNP AVPR1B c.1084A>G on middle lesion score respectively), providing convincing evidence for a genuine effect of the DNA sequence variation of these two genes on stress responsiveness and aggressive behaviour. Consequently, the SNP used here are either directly involved or are in linkage disequilibrium with the causal variants.
In the pig knowledge about the molecular basis of the stress response, aggressive behaviour and their interindividual variation is very limited. In the present study we analyzed the association between DNA markers of ten HPA axis related genes with stress reactivity and aggressive behaviour in the context of psychosocial stress of mixing individuals with different aggressive temperaments. From this we obtained convincing evidence for an effect of two genes: NR3C1 on HPA axis activity and aggressive behaviour, and AVPR1B on aggressive behaviour. Our results provide a foundation for future studies directed at the identification of the causal functional DNA sequence variation which would not only provide markers useful for pig breeding but also insight into the molecular basis of the stress response and aggressive behaviour in general.
The structure and phenotyping of the pig line designated SYN (synthetic) was described in detail by D'Eath et al.. Briefly, pigs of the SYN line were progeny of Pietrain sires and crossbred dams from six different commercial parent lines of the Pig Improvement Company (PIC). The RYR1 SNP c.1843C>T was segregating among the Pietrain sires. Aggressive temperament was measured by counting skin lesions (lesion scoring) immediately before and 24 h after mixing pigs into new groups at approximately 10-11 weeks of age. The increase in number of skin lesions has been shown by Turner et al. to be positively associated with the duration of involvement in reciprocal fighting behaviour. Pigs in each mixed group were ordered by the change in total skin lesions: half of the pigs in each group (those with the most lesions) were designated as high aggressiveness and the remaining half as low aggressiveness. In four slaughter batches pigs were assigned to one of four mixing treatments based on their aggressiveness (number of animals per treatment included in this study: high with high n = 63, high with low n = 61, low with low n = 65, and unmixed n = 64). Pigs were mixed into their treatment groups as they were loaded onto a vehicle for transport to the abattoir. In a fifth treatment, experienced by another four batches, mixing of pigs occurred at loading onto the truck and at lairage in an uncontrolled way, typical of commercial practice (n = 165). Skin lesions were counted before mixing and after slaughter on the carcass, dividing the body into front (head, neck, shoulders and front legs), middle (flanks and back) and rear (rump, hind legs and tail) sections. The difference in lesion number was again taken as the lesion score.
The mixing treatment had significant effect on aggressive behaviour as reported by D'Eath et al.. The animal experiment received approval from the Scottish Agricultural College Animal Experiments Committee.
Pigs were stunned by means of CO2 gas and at exsanguination, a 50 ml sample of trunk blood was collected from each pig in a plastic tube containing 1 ml of 0.5 M EDTA and was stored on ice until plasma preparation, after which they were stored at -80°C. Glucose, lactate and creatine kinase activity were measured with a clinical biochemistry automate (COBAS-MIRA Plus, Roche). Cortisol concentration was measured with the automated analyzer Centaur (Siemens) using a kit designed for human serum and that we validated for pig serum.
The PiF1a line consisted of performance tested pigs (n = 208) of the German commercial cross Pietrain × (German Large White × German Landrace). At slaughter in the FBN experimental slaughterhouse the left adrenal gland was dissected, trimmed of visceral fat and weighed.
Detection of DNA sequence variation and genotyping
Genomic DNA was isolated from skin or liver samples according to the standard phenol-chlorophorm extraction protocol.
DNA sequence variation was detected either in silico (POMC, NR3C1, AVP, CRHR2) by alignment of available porcine sequences and confirmed by direct sequencing or de novo by comparative sequencing of each two individuals of the breeds Pietrain, German Large White, German Landrace and Wild Boar (CRHR1, CRHBP, AVPR1B, UCN). The SNPs MC2R c.306T>G and CRH g.233C>T were published previously by Jacobs et al. and Murani et al. respectively.
Genotyping of the SNPs RYR 1 c.1843C>T and CRH g.233C>T was performed by PCR-RFLP and SSCP respectively as described previously (D'Eath et al. and Murani et al. respectively). The polymorphisms CRHR1 c.*866_867insA, CRHBP c.51G>A, MC2R c.306T>G, NR3C1 c.*2122A>G, AVPR1B c.1084A>G and UCN g.1329T>C were genotyped by PCR-RFLP. Briefly, the polymorphic sites were amplified in 20 μl PCR reactions containing 100 ng genomic DNA, 0.2 mM dNTP, primer as listed in Additional File 1 (Table S1), and 0.5 U SupraTherm Taq-polymerase (Ares Biosciences, Köln, Germany). The temperature profile included initial denaturation at 95°C for 3 min, followed by 40 cycles of denaturation at 95°C for 15 s, annealing at the specific temperature (Additional File 1, Table S1) for 60 s, extension at 72°C for 60 s, and one cycle of final extension at 72°C for 5 min. To detect the polymorphisms 10 μl of the amplified DNA were digested using the appropriate enzyme (Additional File 1, Table S1) overnight according to manufacturer recommendations (Fermentas, St. Leon-Rot, Germany), and the resulting RFLP was analyzed on 2% ethidium bromide stained agarose gel.
The SNP AVP c.207A>G was genotyped by SSCP. The PCR was performed as described above, only scaled down to 10 μl (Additional File 1, Table S2). PCR products were separated on a 12% (49:1 AA:Bis) PAA gel containing 10% urea at 450 V for 4 hours at room temperature and subsequently visualized by silver staining.
The SNP CRHR2 c.*13T>C was genotyped by pyrosequencing. The PCR was performed as described above, only scaled up to 25 μl, using a step-down temperature profile (Additional File 1, Table S3). The subsequent pyrosequencing reaction was performed as described by Srikanchai et al..
The insertion/deletion c.293_298del in POMC was analyzed as a length polymorphism on a MegaBACE 750 capillary sequencer using the MegaBACE Fragment Profiler v1.2 software (GE Healthcare, Munich, Germany). The PCR was performed as described above, only scaled down to 10 μl, using a step-down temperature profile (Additional File 1, Table S4).
Population genetic analyses were performed using the PowerMarker V3.25 software . Association between the candidate gene DNA markers and phenotypic variation was analyzed using general linear model (Proc GLM, SAS V9.1, SAS Institute, Cary, NC, USA). Genotypes of the RYR1 SNP c.1843C>T were considered as a fixed effect in the model for analysis of every other DNA marker because it has been shown to affect stress responsiveness and aggressiveness in previous and also in the present study (see the Results and Discussion section). For the behaviour traits the model included fixed effects of the marker genotype, RYR1 c.1843C>T genotype, slaughter batch, sex and treatment and for the physiological stress parameters also the fixed effect of the dam line. Skin lesion scores and creatine kinase concentrations were log10-transformed before analysis. For adrenal weight the model included fixed effects of the marker genotype, RYR1 c.1843C>T genotype, farm, sex, and body weight as a covariate. Least square means for marker genotypes were compared by a t-test and the p-values were adjusted by Tukey-Kramer correction. False discovery rate (q-value, ) was computed using the JMP Genomics 3 software (SAS Institute, Cary, NC, USA).
Rydiation hybrid physical Mapping
Physical mapping was performed using the INRA-University of Minnesota porcine Radiation Hybrid (IMpRH) panel. The panel was screened by a standard PCR and the products were resolved on 2% agarose gels. The primers and PCR conditions used are detailed in Additional File 1 (Table S5). Regional assignment was obtained using the multipoint analysis option of the IMpRH mapping tool at the IMpRH server http://www.toulouse.inra.fr.
Moberg GP: Biological Response to Stress: Implications for Animal Welfare. The Biology of Animal Stress: Basic Principles and Implications for Animal Welfare. 2000, CABI Publishing; Wallingford; UK, 1-21. full_text.
Fernandez X, Meunier-Salaun MC, Mormede P: Agonistic behavior, plasma stress hormones, and metabolites in response to dyadic encounters in domestic pigs: interrelationships and effect of dominance status. Physiol Behav. 1994, 56: 841-847. 10.1016/0031-9384(94)90313-1.
Otten W, Puppe B, Kanitz E, Schon PC, Stabenow B: Effects of dominance and familiarity on behaviour and plasma stress hormones in growing pigs during social confrontation. Zentralbl Veterinarmed A. 1999, 46: 277-292.
D'Eath RB, Turner SP, Kurt E, Evans G, Tholking L, Looft H, et al: Pigs' aggressive temperament affects pre-slaughter mixing aggression, stress and meat quality. animal. 2010, 4: 604-616. 10.1017/S1751731109991406.
Meese GB, Ewbank R: Exploratory behaviour and leadership in the domesticated pig. Br Vet J. 1973, 129: 251-259.
Tuchscherer M, Puppe B, Tuchscherer A, Kanitz E: Effects of social status after mixing on immune, metabolic, and endocrine responses in pigs. Physiol Behav. 1998, 64: 353-360. 10.1016/S0031-9384(98)00084-5.
Sherritt GW, Graves HB, Gobble JL, Hazlett VE: Effects of Mixing Pigs during the Growing-Finishing Period. J Anim Sci. 1974, 39: 834-837.
Terlouw EMC, Arnould C, Auperin B, Berri C, Le Bihan-Duval E, Deiss V, et al: Pre-slaughter conditions, animal stress and welfare: current status and possible future research. animal. 2008, 2: 1501-1517. 10.1017/S1751731108002723.
Haller J, Kruk MR: Neuroendocrine stress responses and aggression. 2003, Humana Press Inc
Touma C, Bunck M, Glasl L, Nussbaumer M, Palme R, Stein H, et al: Mice selected for high versus low stress reactivity: a new animal model for affective disorders. Psychoneuroendocrinology. 2008, 33: 839-862. 10.1016/j.psyneuen.2008.03.013.
Kruk MR, Halasz J, Meelis W, Haller J: Fast positive feedback between the adrenocortical stress response and a brain mechanism involved in aggressive behavior. Behav Neurosci. 2004, 118: 1062-1070. 10.1037/0735-7044.118.5.1062.
Desautes C, Sarrieau A, Caritez JC, Mormede P: Behavior and pituitary-adrenal function in large white and Meishan pigs. Domest Anim Endocrinol. 1999, 16: 193-205. 10.1016/S0739-7240(99)00014-4.
Mormede P, Courvoisier H, Ramos A, Marissal-Arvy N, Ousova O, Desautes C, et al: Molecular genetic approaches to investigate individual variations in behavioral and neuroendocrine stress responses. Psychoneuroendocrinology. 2002, 27: 563-583. 10.1016/S0306-4530(01)00093-2.
Turner SP, Roehe R, D'Eath RB, Ison SH, Farish M, Jack MC, et al: Genetic validation of postmixing skin injuries in pigs as an indicator of aggressiveness and the relationship with injuries under more stable social conditions. J Anim Sci. 2009, 87: 3076-3082. 10.2527/jas.2008-1558.
Nelson RJ, Trainor BC: Neural mechanisms of aggression. Nat Rev Neurosci. 2007, 8: 536-546. 10.1038/nrn2174.
Derijk RH: Single nucleotide polymorphisms related to HPA axis reactivity. Neuroimmunomodulation. 2009, 16: 340-352. 10.1159/000216192.
Feder A, Nestler EJ, Charney DS: Psychobiology and molecular genetics of resilience. Nat Rev Neurosci. 2009, 10: 446-457. 10.1038/nrn2649.
Moisan MP: Genotype-phenotype associations in understanding the role of corticosteroid-binding globulin in health and disease animal models. Mol Cell Endocrinol. 2010, 316: 35-41. 10.1016/j.mce.2009.07.017.
Desautes C, Bidanelt JP, Milant D, Iannuccelli N, Amigues Y, Bourgeois F, et al: Genetic linkage mapping of quantitative trait loci for behavioral and neuroendocrine stress response traits in pigs. J Anim Sci. 2002, 80: 2276-2285.
Ousova O, Guyonnet-Duperat V, Iannuccelli N, Bidanel JP, Milan D, Genet C, et al: Corticosteroid binding globulin: a new target for cortisol-driven obesity. Mol Endocrinol. 2004, 18: 1687-1696. 10.1210/me.2004-0005.
Geldermann H, Muller E, Moser G, Reiner G, Bartenschlager H, Cepica S, et al: Genome-wide linkage and QTL mapping in porcine F-2 families generated from Pietrain, Meishan and Wild Boar crosses. Journal of Animal Breeding and Genetics. 2003, 120: 363-393. 10.1046/j.0931-2668.2003.00408.x.
Fujii J, Otsu K, Zorzato F, de LS, Khanna VK, Weiler JE, et al: Identification of a mutation in porcine ryanodine receptor associated with malignant hyperthermia. Science. 1991, 253: 448-451. 10.1126/science.1862346.
Weaver SA, Dixon WT, Schaefer AL: The effects of mutated skeletal ryanodine receptors on hypothalamic-pituitary-adrenal axis function in boars. J Anim Sci. 2000, 78: 1319-1330.
Li LA, Xia D, Wei S, Li X, Parvizi N, Zhao RQ: Diminished expression of ACTH signaling proteins and steroidogenic limiting factors in adrenocortical cells isolated from halothane nn pigs. Domest Anim Endocrinol. 2008, 35: 1-7. 10.1016/j.domaniend.2007.12.005.
Caldwell HK, Lee HJ, Macbeth AH, Young WS: Vasopressin: behavioral roles of an "original" neuropeptide. Prog Neurobiol. 2008, 84: 1-24. 10.1016/j.pneurobio.2007.10.007.
D'Eath RB, Ormandy E, Lawrence AB, Sumner BE, Meddle SL: Resident-intruder trait aggression is associated with differences in lysine vasopressin and serotonin receptor 1A (5-HT1A) mRNA expression in the brain of pre-pubertal female domestic pigs (Sus scrofa). J Neuroendocrinol. 2005, 17: 679-686. 10.1111/j.1365-2826.2005.01359.x.
Steckler T: The molecular neurobiology of stress--evidence from genetic and epigenetic models. Behav Pharmacol. 2001, 12: 381-427.
Bale TL: Sensitivity to stress: dysregulation of CRF pathways and disease development. Horm Behav. 2005, 48: 1-10. 10.1016/j.yhbeh.2005.01.009.
Murani E, Muraniova M, Ponsuksili S, Schellander K, Wimmers K: Molecular characterization and evidencing of the porcine CRH gene as a functional-positional candidate for growth and body composition. Biochem Biophys Res Commun. 2006, 342: 394-405. 10.1016/j.bbrc.2006.01.143.
Jacobs K, Van Poucke PM, Mattheeuws M, Chardon P, Yerle M, Rohrer G, et al: Characterization of the porcine melanocortin 2 receptor gene (MC2R). Anim Genet. 2002, 33: 415-421. 10.1046/j.1365-2052.2002.00899.x.
Ramirez O, Ojeda A, Tomas A, Gallardo D, Huang LS, Folch JM, et al: Integrating Y-chromosome, mitochondrial, and autosomal data to analyze the origin of pig breeds. Mol Biol Evol. 2009, 26: 2061-2072. 10.1093/molbev/msp118.
Albert FW, Shchepina O, Winter C, Rompler H, Teupser D, Palme R, et al: Phenotypic differences in behavior, physiology and neurochemistry between rats selected for tameness and for defensive aggression towards humans. Horm Behav. 2008, 53: 413-421. 10.1016/j.yhbeh.2007.11.010.
Trut L, Oskina I, Kharlamova A: Animal evolution during domestication: the domesticated fox as a model. Bioessays. 2009, 31: 349-360. 10.1002/bies.200800070.
Weiler U, Claus R, Schnoebelen-Combes S, Louveau I: Influence of age and genotype on endocrine parameters and growth performance: a comparative study in Wild boars, Meishan and Large White boars. Livestock Production Science. 1998, 54: 21-31. 10.1016/S0301-6226(97)00165-6.
Foury A, Tribout T, Bazin C, Billon Y, Bouffaud M, Gogue JM, et al: Estimation of genetic trends from 1977 to 2000 for stress-responsive systems in French Large White and Landrace pig populations using frozen semen. animal. 2009, 3: 1681-1687. 10.1017/S1751731109990504.
Warriss PD, Brown SN: The Physiological-Responses to Fighting in Pigs and the Consequences for Meat Quality. Journal of the Science of Food and Agriculture. 1985, 36: 87-92. 10.1002/jsfa.2740360207.
Turner SP, Farnworth MJ, White IMS, Brotherstone S, Mendl M, Knap P, et al: The accumulation of skin lesions and their use as a predictor of individual aggressiveness in pigs. Applied Animal Behaviour Science. 2006, 96: 245-259. 10.1016/j.applanim.2005.06.009.
Turner SP, Roehe R, Mekkawy W, Farnworth MJ, Knap PW, Lawrence AB: Bayesian analysis of genetic associations of skin lesions and behavioural traits to identify genetic components of individual aggressiveness in pigs. Behavior Genetics. 2008, 38: 67-75. 10.1007/s10519-007-9171-2.
Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci USA. 2003, 100: 9440-9445. 10.1073/pnas.1530509100.
Kuchenmeister U, Nurnberg K, Fiedler I, Kuhn G, Nurnberg G, Ender K: Cell injury and meat duality of pig in the time period post mortem from two genotypes susceptible or resistant to malignant hyperthermia. European Food Research and Technology. 1999, 209: 97-103. 10.1007/s002170050465.
Shen QW, Underwood KR, Means WJ, McCormick RJ, Du M: The halothane gene, energy metabolism, adenosine monophosphate-activated protein kinase, and glycolysis in postmortem pig longissimus dorsi muscle. J Anim Sci. 2007, 85: 1054-1061. 10.2527/jas.2006-114.
Guárdia MD, Estany J, Balasch S, Oliver MA, Gispert M, Diestre A: Risk assessment of skin damage due to pre-slaughter conditions and RYR1 gene in pigs. Meat Science. 2009, 81: 745-751. 10.1016/j.meatsci.2008.11.020.
Parrott RF, Vellucci SV, Goode JA: Behavioral and hormonal effects of centrally injected "anxiogenic" neuropeptides in growing pigs. Pharmacol Biochem Behav. 2000, 65: 123-129. 10.1016/S0091-3057(99)00134-3.
Smagin GN, Heinrichs SC, Dunn AJ: The role of CRH in behavioral responses to stress. Peptides. 2001, 22: 713-724. 10.1016/S0196-9781(01)00384-9.
Deussing JM, Wurst W: Dissecting the genetic effect of the CRH system on anxiety and stress-related behaviour. C R Biol. 2005, 328: 199-212. 10.1016/j.crvi.2005.01.001.
Veenema AH, Neumann ID: Neurobiological mechanisms of aggression and stress coping: a comparative study in mouse and rat selection lines. Brain Behav Evol. 2007, 70: 274-285. 10.1159/000105491.
Karpac J, Ostwald D, Bui S, Hunnewell P, Shankar M, Hochgeschwender U: Development, maintenance, and function of the adrenal gland in early postnatal proopiomelanocortin-null mutant mice. Endocrinology. 2005, 146: 2555-2562. 10.1210/en.2004-1290.
Parton LE, Ye CP, Coppari R, Enriori PJ, Choi B, Zhang CY, et al: Glucose sensing by POMC neurons regulates glucose homeostasis and is impaired in obesity. Nature. 2007, 449: 228-232. 10.1038/nature06098.
Yue G, Stratil A, Kopecny M, Schroffelova D, Schroffel J, Hojny J, et al: Linkage and QTL mapping for Sus scrofa chromosome 6. Journal of Animal Breeding and Genetics. 2003, 120: 45-55. 10.1046/j.0931-2668.2003.00423.x.
Zhang J, Ge R, Matte-Martone C, Goodwin J, Shlomchik WD, Mamula MJ, et al: Characterization of a novel gain of function glucocorticoid receptor knock-in mouse. J Biol Chem. 2009, 284: 6249-6259. 10.1074/jbc.M807997200.
Perreau V, Sarrieau A, Mormede P: Characterization of mineralocorticoid and glucocorticoid receptors in pigs: comparison of Meishan and Large White breeds. Life Sci. 1999, 64: 1501-1515. 10.1016/S0024-3205(99)00087-9.
Wersinger SR, Caldwell HK, Christiansen M, Young WS: Disruption of the vasopressin 1b receptor gene impairs the attack component of aggressive behavior in mice. Genes Brain Behav. 2007, 6: 653-660. 10.1111/j.1601-183X.2006.00294.x.
Weninger SC, Peters LL, Majzoub JA: Urocortin expression in the Edinger-Westphal nucleus is up-regulated by stress and corticotropin-releasing hormone deficiency. Endocrinology. 2000, 141: 256-263. 10.1210/en.141.1.256.
Bale TL, Vale WW: CRF and CRF receptors: role in stress responsivity and other behaviors. Annu Rev Pharmacol Toxicol. 2004, 44: 525-557. 10.1146/annurev.pharmtox.44.101802.121410.
Zalutskaya AA, Arai M, Bounoutas GS, Abou-Samra AB: Impaired adaptation to repeated restraint and decreased response to cold in urocortin 1 knockout mice. Am J Physiol Endocrinol Metab. 2007, 293: E259-E263. 10.1152/ajpendo.00616.2006.
Moreau JL, Kilpatrick G, Jenck F: Urocortin, a novel neuropeptide with anxiogenic-like properties. Neuroreport. 1997, 8: 1697-1701. 10.1097/00001756-199705060-00027.
Tanaka C, Asakawa A, Ushikai M, Sakoguchi T, Amitani H, Terashi M, et al: Comparison of the anorexigenic activity of CRF family peptides. Biochem Biophys Res Commun. 2009, 390: 887-891. 10.1016/j.bbrc.2009.10.069.
Kuperman Y, Chen A: Urocortins: emerging metabolic and energy homeostasis perspectives. Trends Endocrinol Metab. 2008, 19: 122-129. 10.1016/j.tem.2007.12.002.
Guyonnet-Duperat V, Geverink N, Plastow GS, Evans G, Ousova O, Croisetiere C, et al: Functional Implication of an Arg307Gly Substitution in Corticosteroid-Binding Globulin, a Candidate Gene for a Quantitative Trait Locus Associated With Cortisol Variability and Obesity in Pig. Genetics. 2006, 173: 2143-2149. 10.1534/genetics.105.053983.
Niu N, Manickam V, Kalari KR, Moon I, Pelleymounter LL, Eckloff BW, et al: Human glucocorticoid receptor alpha gene (NR3C1) pharmacogenomics: gene resequencing and functional genomics. J Clin Endocrinol Metab. 2009, 94: 3072-3084. 10.1210/jc.2008-2109.
Srikanchai T, Murani E, Phatsara C, Schwerin M, Schellander K, Ponsuksili S, et al: Association of ZYX polymorphisms with carcass and meat quality traits in commercial pigs. Meat Science. 2010, 84: 159-164. 10.1016/j.meatsci.2009.08.042.
Liu K, Muse SV: PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics. 2005, 21: 2128-2129. 10.1093/bioinformatics/bti282.
The authors thank Angela Garve, Marlies Fuchs, and Sarah H. Ison for excellent technical help. The Muong Khuong samples were kindly provided by Dr. Do Vo Anh Khoa. We are grateful to Peter Andresen for allowing his pigs and farm to be used for the study of SYN pigs.
This work was part-funded by 6th Research Framework Programme of the European Union European Commission as part of the SABRE [Cutting-Edge Genomics for Sustainable Animal Breeding] Integrated Project. FBN, PIC and INRA provided their own matched funding. SAC matched funding was provided by the Rural and Environment Research and Analysis Directorate of the Scottish Government.
EM performed the statistic analysis and drafted the manuscript. SP assisted the statistic analysis, participated in sampling and assisted drafting the manuscript and revising it critically for scientific content. RBD, SPT and EK designed the mixing experiment and performed lesion scoring. GE, LT and RK provided input to the design of the mixing experiment, and participated in sampling and phenotyping. AF and PM performed measurements of the physiological stress parameters. KW significantly contributed to the concept, design and coordination of the study and in drafting the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
About this article
Cite this article
Muráni, E., Ponsuksili, S., D'Eath, R.B. et al. Association of HPA axis-related genetic variation with stress reactivity and aggressive behaviour in pigs. BMC Genet 11, 74 (2010). https://doi.org/10.1186/1471-2156-11-74