Skip to main content

Effects of divergent selection upon adrenocortical activity on immune traits in pig



The sustainability of farming and animal welfare requires the reconsideration of current selection schemes. In particular, implementation of new selection criteria related to animal health and welfare should help to produce more robust animals and to reduce anti-microbial use. The hypothalamo-pituitary-adrenocortical (HPA) axis plays a major role in metabolic regulation and adaptation processes and its activity is strongly influenced by genetic factors. A positive association between HPA axis activity and robustness was recently described. To explore whether selecting pigs upon HPA axis activity could increase their robustness, a divergent selection experiment was carried out in the Large White pig breed. This allowed the generation of low (HPAlo) and high (HPAhi) responders to adrenocorticotropic hormone administration.


In this study, we compared 23 hematologic and immune parameters of 6-week-old, HPAlo and HPAhi piglets and analysed their response to a low dose of lipopolysaccharide (LPS) two weeks later. At six weeks of age, HPAhi piglets displayed greater red blood cell and leucocyte number including CD8α+ γδ cells, cytotoxic T lymphocytes, naive T helper (Th) cells and B lymphocytes as compared to HPAlo individuals. The ability of blood cells to secrete TNFα in response to LPS ex vivo was higher for HPAhi pigs. At week eight, the inflammatory response to the LPS in vivo challenge was poorly affected by the HPA axis activity.


Divergent selection upon HPA axis activity modulated hematologic and immune parameters in 6-week-old pigs, which may confer an advantage to HPAhi pigs at weaning. However, HPAlo and HPAhi piglets did not exhibit major differences in the parameters analysed two weeks later, i. e. in 8-week-old pigs. In conclusion, chronic exposure to high cortisol levels in HPAhi pigs does not negatively impact immunity.


For decades, genetic selection in the pig industry focused on production criteria with deleterious consequences on health and welfare (reviewed in [1]). The sustainability of farming, including animal welfare, requires reconsideration of selection schemes. Implementation of new genetic selection criteria in relation to pig health could help reducing antibiotic use, which is urgently needed to address public health concerns [2, 3] and is strongly encouraged by the European Commission (for details, see [4]). In particular, selecting pigs for robustness could increase animal resilience to pathogens.

A positive association between HPA axis activity and robustness was previously described [5]. The HPA axis plays a major role in metabolic regulation and adaptation processes and its activity is strongly influenced by genetic factors [6, 7]. Together with the autonomous nervous system, it helps the organism to cope with threats and infections. Cortisol is the main terminal effector of the HPA axis in humans and pigs. Its secretion by adrenal glands is controlled by the release of the corticotropin-releasing factor by the hypothalamus, which induces subsequent adrenocorticotropic hormone (ACTH) secretion by the anterior pituitary gland. Adrenal gland sensitivity to ACTH was demonstrated to differ largely among individuals while being a highly heritable trait [8]. We hypothesized that selecting pigs upon the adrenal response to ACTH could increase their robustness. To explore this strategy, a divergent selection experiment was carried out in the Large White pig breed based on plasma cortisol levels measured one hour after ACTH administration [8, 9]. In second generation offsprings, we showed that high- (HPAhi) and low- (HPAlo) responders to ACTH stimulation differ in their ability to cope with an acute social stress, with HPAhi piglets appearing more resilient [10]. We also highlighted some differences in immune traits between the two populations of pigs. Indeed, HPAhi pigs displayed greater CD4+ T cell counts and their blood cells secreted higher amounts of TNFα after LPS stimulation.

In third generation offsprings, the post-ACTH plasma cortisol level was found to be 2.16 times higher in the HPAhi line than in the HPAlo line [9]. In the current study, we aimed to compare hematologic and immune parameters of 16 HPAhi and 16 HPAlo piglets belonging to this third generation of divergent selection at six (experiment 1) and eight (experiment 2) weeks of age.


Effects of HPA axis activity on immune traits in pigs

In experiment 1, hematologic and immune parameters of 6-week-old HPAlo (n = 16) and HPAhi (n = 16) pigs were assessed. As expected, pigs exhibited significantly different steady-state cortisolemia with HPAhi piglets secreting twice as much cortisol than HPAlo pigs (Fig. 1).

Fig. 1

Steady-state cortisolemia in HPAlo and HPAhi piglets. Plasma cortisol levels were measured in 6-week-old HPAlo (empty squares) and HPAhi (black squares) piglets (n = 16 per group). Each point represents an individual and medians are shown. Stars denote statistical significance (p < 0.001)

To get an overview of our data, we first performed a Principal Component Analysis (PCA) on raw data. Individuals are mapped along the first two dimensions of the PCA (Fig. 2), while contributions of individual variables to the first two dimensions are presented in Additional file 1: Figure S1 and Additional file 2: Table S1. The first axis was mainly determined by white blood cell parameters (CD4+ and CD8αhi T cells), while the second axis mainly resulted from red blood cell parameters (RBC, HGB and HCT). HPAlo and HPAhi piglet populations appeared clearly separated, as shown by their confidence ellipses, especially along the “white blood cell” axis.

Fig. 2

Principal component analysis: individuals factor map. PCA was performed using values obtained for 20 variables measured in 6-week-old HPAlo (black dots) and HPAhi (red dots) pigs (n = 16 per group). Each point represents a single individual in the first two dimensions of the PCA

While divergent selection did not modify platelet counts (data not shown), some red and white blood cell parameters appeared to be affected by selection of pigs upon HPA axis activity. As shown in Fig. 3, among seven red blood cell parameters measured, five significantly differed between HPAlo and HPAhi pigs, namely red blood cell count (RBC), mean corpuscular volume (MCV), haemoglobin (HGB), reticulocyte count (Ret) and immature reticulocyte fraction (IRF). While RBC, HGB and IRF were higher in HPAhi pigs, MCV and Ret were lower as compared to HPAlo individuals.

Fig. 3

Red blood cell parameters in HPAlo and HPAhi pigs. Red blood cell parameters (RBC, HCT, MCV, HGB, MCHC, Ret and IRF) were measured in EDTA-treated blood of HPAlo (empty squares) and HPAhi (black squares) piglets (n = 16 per group) using the Procyte™ hematimeter. Each point represents a single piglet. Stars denote statistical significance (*p < 0.05; ** p < 0.01)

Concerning white blood cell parameters, leucocyte count was significantly higher in HPAhi as compared to HPAlo pigs at six weeks of age (p < 0.01, data not shown). This difference in leucocyte counts essentially relied on lymphocyte counts since neutrophil and monocyte counts did not vary significantly between the two groups (Fig. 4). To better characterise the difference in lymphocyte numbers, we used flow cytometry, which allowed the discrimination of seven lymphocyte subpopulations (Additional file 3: Figure S2 and not shown for B cells). Four lymphocyte subpopulations appeared higher in HPAhi as compared to HPAlo individuals: CD8 α+ γδ T cells, CD8αhi cytotoxic T lymphocytes, naive CD4+ Th cells and CD21+ B lymphocytes (Fig. 5). Of note, CD8α γδ cell, NK cell and antigen-experienced Th lymphocyte numbers were similar between the two groups.

Fig. 4

White blood cell counts of both lines of pigs. Neutrophil, monocyte and lymphocyte cell counts were measured in EDTA-treated blood of HPAlo (empty squares) and HPAhi (black squares) piglets (n = 16 per group) using the Procyte™ hematimeter. Each point represents a single piglet. Stars denote statistical significance (** p < 0.01)

Fig. 5

Specific variations in lymphocyte subpopulations of HPAlo and HPAhi piglets. Lymphocyte subpopulations were characterized in EDTA-treated blood of 6-week old HPAlo (empty squares) and HPAhi (black squares) piglets (n = 16 per group) by multicolor staining using flow cytometry. Seven populations were identified: B lymphocytes (CD3 CD21+), naive CD4+ Th cells (CD3+ CD4+), antigen-experienced Th cells (CD3+ CD4+ CD8α+), CD8α γδ T cells (CD3+ CD4 CD8α), CD8α+ γδ T cells (CD3+ CD4 CD8αmed), CD8αhi cytotoxic T lymphocytes (CD3+ CD4 CD8αhi) and NK cells (CD3 CD4 CD8αmed); gating strategy is shown in Additional file 3: Figure S2. Each point represents a single piglet. Stars denote statistical significance (*p < 0.05; ** p < 0.01; ***p < 0.001)

To investigate whether selection upon HPA axis strength also modulated leucocyte functions in 6-week-old pigs, we performed ex-vivo phagocytosis and whole blood assays (Fig. 6). Mononuclear and polynuclear cell phagocytosis were similar for HPAlo and HPAhi pigs (Fig. 6a). Concerning LPS-induced cytokine secretion by blood cells, while IL-8 and IL-10 were not different between the two groups, TNFα secretion from cells derived from pigs with strong HPA axis activity was significantly higher (Fig. 6b).

Fig. 6

Functional ability of circulating immune cells in both lines of pigs. Phagocytosis (a) and cytokine secretion (b) were evaluated on heparin-treated blood of 6-week old HPAlo (empty squares) and HPAhi (black squares) pigs (n = 16 per group). Each point represents a single piglet. a, Phagocytosis of fluorescent E. coli. Bacteria by mononuclear and polynuclear cells was analysed by flow cytometry. b, Cytokine secretion by whole blood cells was evaluated after an ex vivo stimulation for 18-h with 10 ng/ml LPS. The amount of TNFα, IL-8 and IL-10 in the supernatants was measured by ELISA. Stars denote statistical significance (** p < 0.01)

Effects of HPA axis activity on acute response to LPS in 8-week-old pigs

Experiment 2 was performed on the same animals two weeks later. We administered a low dose of LPS (15 μg/kg) and investigated its consequences on several parameters focusing on short-term effects, i.e. within 24 h after intramuscular LPS injection. While this protocol is not known to induce sickness behaviour (vomiting, diarrhoea and/or somnolence), one pig from the HPAhi line died between 4 and 24 h after LPS injection and was hence excluded from the analysis. In pigs from both populations, LPS induced a rise in tympanic temperature, although less intense in HPAhi pigs (Fig. 7, Table 1).

Fig. 7

Effect of LPS injection on body temperature, cortisol, TNFα and IL-6 plasma levels in 8-week-old pigs. Blood was collected in heparin-coated tubes from 8-week old HPAlo (empty squares, n = 16) and HPAhi (black squares, n = 15) pigs just before (t0) and 1, 4 and 24 h following the i.m. injection of LPS. At each time point, temperature was measured into the ear. Cortisol, TNFα and IL-6 levels were measured in plasma samples. Each point represents a single piglet. See Table 1 for statistical analysis

Table 1 Results of the statistical analysis of the in vivo experiment (p-values)

Also, cortisolemia rapidly increased after LPS injection with a peak at 4 h (Fig. 7). As expected, HPAhi pigs released significantly greater amounts of cortisol at every time point after LPS challenge than HPAlo pigs did (Table 1). In nearly all HPAlo and HPAhi individuals, stress response to LPS was accompanied by a decrease in glycaemia and an increase in free-fatty acid release and bilirubinaemia. The profiles of white blood cell count and lymphocyte-to-granulocyte ratio in response to LPS were similar between pigs from both populations (Table 1, p = 0.052 and 0.92 respectively).

A transcriptomic analysis was also conducted to monitor divergent selection effects on the expression of 34 genes modulated by LPS. Among them, 23 genes were up-regulated at t1 h, 13 genes were down-regulated at t1 h, 1 gene was up-regulated at t4 h, 21 genes were down-regulated at t4 h and 4 genes were down-regulated at t24 h (Additional file 4: Table S2). No differences in LPS-modulated transcript levels were observed between the two populations of pigs at any time point (data not shown). Finally, plasma TNFα and IL-6 levels after LPS administration were measured. Typical secretion curves were obtained, peaking at t1 h and t4 h, respectively (Fig. 7). While TNFα levels did not differ between individuals from the two populations, plasma IL-6 levels were higher in HPAlo as compared to HPAhi pigs (Table 1, p = 0.013). Strikingly, two full-sib individuals of the HPAhi group behaved differently and responded poorly to the LPS challenge.


In this study, we aimed to document on the effects of divergent genetic selection upon HPA axis activity on blood cell counts and functions by comparing data obtained from 16 HPAlo and 16 HPAhi pigs in basal conditions (experiment 1) and after LPS challenge (experiment 2), respectively at 6 and 8 weeks of age.

At 6 weeks of age, HPAhi pigs displayed a higher red blood cell count and haemoglobin content, and a lower mean corpuscular volume. HPAhi individuals also exhibited a lower reticulocyte number associated with a higher immature reticulocyte fraction, which reflects the proportion of younger reticulocytes [11]. The higher red blood cell count in association with a greater immature reticulocyte fraction suggests a stronger erythropoietin activity in HPAhi pigs, which could increase the oxygen carrying capacity of blood (reviewed in [12]). Interestingly, the plasma cortisol level was recently shown to be positively correlated to the erythropoietin (EPO) level in healthy human volunteers [13]. EPO is a glycoprotein controlling early events of erythropoiesis in the bone marrow thereby increasing the number of red blood cells [14]. In this context, it could be interesting to measure serum EPO levels in piglets from both populations. Whether the higher oxygen carrying capacity of blood could boost HPAhi piglet growth and performance is currently under investigation.

While moderate and short-term elevations of circulating glucocorticoids (GCs) may promote immunity, high and long-term elevations of circulating GCs were rather shown to be immunosuppressive (reviewed in [15]). In particular, GCs induce lymphocyte sequestration into tissues thereby reducing their number in the circulation [16, 17]. While studying diurnal rhythm effects on blood cell numbers in pigs, Engert et al. recently underscored negative associations between plasma cortisol levels and most lymphocyte subtypes [18]. However, in our study, despite a higher cortisolemia, 6-week-old HPAhi piglets displayed a higher leucocyte count, which mainly relied on B lymphocytes, naive Th lymphocytes, cytotoxic T lymphocytes and CD8α+ γδ T cells. Whether this increase in circulating lymphocytes may foster the immune capacity of HPA pigs is still unknown. We next evaluated the “innate immune response signature” of HPAlo and HPAhi piglets by assessing LPS-induced TNFα, IL-8 and IL-10 secretion using whole blood assays (WBA) and phagocytosis of opsonized E. coli bacteria. While most functions were preserved in HPAhi piglets, TNFα secretion was significantly higher in 6-week-old HPAhi pigs as compared to HPAlo individuals, confirming our previous results obtained in animals from the second generation of selection [10]. In this previous experiment, HPAhi piglets exhibited a higher naive Th lymphocyte number and greater LPS-induced TNFα secretion than did HPAlo pigs did. Interestingly, after another round of divergent genetic selection, we also evidenced an increase in three other lymphocyte subpopulations, namely cytotoxic T lymphocytes, CD8α+ γδ T cells and B lymphocytes, in 6-week-old, HPAhi pigs. Given their immaturity, young piglets are more fragile, weaning being the most critical period of their lives. From the first experiment, we presume that HPAhi piglets might be better prepared to cope with stress including exposure to pathogens in the first months of their life.

We next evaluated metabolic and immune responses to in vivo LPS challenge at week 8. In both lines, LPS induced activation of the HPA axis, as shown by the elevation of plasma cortisol levels following injection, confirming previous observations in unselected pigs [19]. The LPS-induced metabolic changes also consisted in an increase in free fatty acid release and bilirubinaemia associated to a decrease in glycaemia, probably due to increased glycolysis and decreased neoglucogenesis. No difference between HPAlo and HPAhi pigs was evidenced for most immune parameters observed, except for plasma IL-6 levels, which were higher in the HPAlo line. Thus, the divergent genetic selection program poorly affected the LPS response of 8-week-old pigs, which is in agreement with our previous work. Indeed, in unselected Large White pigs, responses to LPS were not influenced by adrenal gland reactivity as assessed by the cortisol response to ACTH.

Interestingly, in 8-week-old pigs, the difference in baseline blood cell counts previously observed between the two lines of pigs had disappeared (data not shown). Chronic exposure to high cortisol levels is known to induce GC resistance as an adaptive process, notably through down-regulation of glucocorticoid receptors (GR) expression. Age-related effects on lymphocyte cortisol sensitivity were also documented in pigs [20].

In addition, in the case of this divergent selection experiment, exposure to differential GC levels might have led to select individuals with different intrinsic sensitivity to GC. A Genome Wide Association Study conducted on siblings actually revealed that the selection process segregated two alleles of the NR3C1 gene, which encodes GR [9]. HPAlo pigs, in contrast to HPAhi pigs, carry the previously described mutated form of swine NR3C1 gene [21]. This allele encodes an A610V substitution in the ligand-binding domain of GR, which induces a gain of function of the receptor leading to GC hypersensitivity. Pigs carrying this mutation did not exhibit major phenotypic differences in growth and body composition, but they displayed a transient lower leucocyte number compared to the carriers of the wild type gene [22]. In this context, we can wonder whether the overall exposure of HPAlo and HPAhi pigs to GC is different.


In this study, we confirmed that young HPAhi piglets exhibit greater leucocyte counts and TNFα secretion in whole-blood assay than HPAlo pigs, but these differences decreased with age. We also demonstrated that chronic exposure to high cortisol levels in HPAhi pigs does not negatively impact immunity.


Animals and experimental design

All experimental procedures were carried out in accordance with the EU directive 2010/63 for animal experiments and submitted to and approved by the C2EA 84 Poitou-Charentes Region ethics committee (decision #CE2013–1, 2013.01.21).

A divergent genetic selection experiment was conducted in purebred Large White pigs by selecting, at each generation, founders with extreme cortisolemia following ACTH stimulation as an indicator of HPA axis activity. This allowed the generation of pigs low (HPAlo) and high (HPAhi) responders to ACTH injection [9]], which were all born and raised at the INRA experimental unit of Le Magneraud (UE1372 GenESI) until euthanasia by lethal electronarcosis.

This study involved blood samples collected from 32 healthy female pigs belonging to the third generation of divergent selection (16 HPAlo and 16 HPAhi pigs). In each population, the 16 pigs were offspring of 8 sows mated with 8 different barrows, constituting 8 pairs of full-sibs.

After weaning at 4 weeks, pigs were housed by groups of 20 to 22 individuals per pen in controlled ambient conditions (27 °C, 12 h/12 h light/dark cycle) with unrestricted access to food and water. Two experiments were conducted. At week 6, hematologic and immune parameters were analysed on blood samples collected from 9 to 10 a.m (Experiment 1). Two weeks later, a kinetic response to intra-muscular LPS injection (15 μg/kg in the neck muscles) was performed on the same animals with blood samples collected just before (t0, between 9 and 10 a.m.) and 1, 4 and 24 h after LPS administration as previously described [19] (Experiment 2).

In both experiments, blood was collected from the jugular vein of pigs slightly maintained in a supine position by experienced staff, in less than 30 s per pig as previously described, with a total volume collected never exceeding 25 ml per pig [23]. All blood samples were processed by blinded investigators.

Complete blood cell counts

Complete blood cell counts (CBC) were determined from EDTA-treated blood samples using the clinical-grade Procyte Dx haematology analyser (IDEXX, Saint-Denis, France). Red blood cell parameters included red blood cell count (RBC), haematocrit (HCT), haemoglobin (HGB), mean corpuscular volume (MCV), mean corpuscular haemoglobin concentration (MCHC), reticulocyte count (Ret) and immature reticulocyte fraction (IRF). CBC also included platelet count as well as white blood cell parameters: white blood cell count (WBC), neutrophil, lymphocyte and monocyte counts.

Analysis of lymphocyte subpopulations

Lymphocyte subpopulations were identified from EDTA-treated blood samples by two types of staining using monoclonal antibodies: (1) FITC-conjugated anti-pig CD4 (Acrys Antibodies, Herford, Germany), PE-conjugated anti-pig CD8α (BD Biosciences) and PerCP-conjugated anti-pig CD3 (Bio-Techne, Lille, France) antibodies, or (2) PE-conjugated anti-pig CD21 (Acrys Antibodies) and PerCP-conjugated anti-pig CD3 in the dark for 15 min at room temperature. Erythrocytes were removed using a lysis solution (BD Biosciences, Le Pont de Claix, France). The DNA marker DRAQ5™ (1 μM final; Biostatus, Leicestershire, UK) was added to identify nucleated cells, then samples were analysed with a FACS Aria flow cytometer (BD Biosciences) and data computed using the FlowJo software (FlowJo, Ashland, Oregon). Staining allowed the identification of seven different lymphocyte subpopulations as previously described [18]: B lymphocytes (CD3 CD21+), naive CD4+ Th cells (CD3+ CD4+), antigen-experienced Th cells (CD3+ CD4+ CD8α+), CD8α γδ T cells (CD3+ CD4 CD8α), CD8α+ γδ T cells (CD3+ CD4 CD8αmed), CD8αhi cytotoxic T lymphocytes (CD3+ CD4 CD8αhi) and NK cells (CD3 CD4 CD8αmed). The cell count for each specific subpopulation was then determined as the product of the percentage of this specific subpopulation in the DRAQ5+ FSC/SSC white blood cell gate by the absolute number of white blood cells obtained from the Procyte analyser.


Ex vivo phagocytosis was assessed using the Phagotest™ Kit (BD Biosciences) according to manufacturer’s instructions. Briefly, heparinized blood samples were incubated for 10 min at 37 °C with opsonized FITC-labelled E. coli bacteria. For each analysis, ice-incubated negative controls were included. Samples were immediately analysed by flow cytometry using FSC/SSC dot plot to discriminate granulocytes from mononuclear cells.

Whole blood assay (WBA)

Heparinized blood samples were five-fold diluted in RPMI 1640 medium supplemented with 2 mM L-glutamine, 100 IU/ml penicillin and 100 mg/ml streptomycin and stimulated in duplicate using 10 ng/ml Escherichia coli O111:B4 LPS (Sigma-Aldrich, St-Quentin-Fallavier, France). After an 18-h incubation at 37 °C, in 5% CO2, 95% relative humidity, supernatants were collected, centrifuged at 500 g for 5 min at 4 °C and stored at − 80 °C before cytokine concentration measurements. Porcine TNFα, IL-8 and IL-10 were quantified by ELISA with respective detection limits of 15, 62 and 12 pg/ml (Bio-Techne). For statistical analyses, we attributed to undetected samples half of the corresponding threshold value.

Kinetic response to LPS

Cortisol, glucose and free fatty acid plasma levels were assessed after LPS administration as previously described [23].

Plasma cytokine levels in response to LPS injection were quantified in heparin-treated samples by ELISA (porcine IL-6 and TNFα, Bio-Techne).

In order to investigate divergent selection effects on gene expression response to LPS, we tested expression of 34 genes modulated by LPS. These genes were provided by the literature and from our own data ([23] and unpublished data). Gene expression analysis was performed as previously described [23]. Total RNA extraction was done using the Nucleospin RNA Blood kit (Macherey-Nagel, Hoerdt, France) followed by DNase treatment. The quality of each RNA preparation was verified through the Bioanalyser Agilent 2100 (Agilent Technologies, Massy, France) and low-quality RNA samples were discarded (RIN < 8). One μg of total RNA was reverse-transcribed. The transferrin receptor and glutamyl-prolyl-tRNA synthetase genes were used as housekeeping genes, according to previous data. Pre-amplified samples were analysed with a 96 × 96 Dynamic Array™ IFC (Fluidigm) following the protocol defined by Spurgeon et al. [24]. Each gene was tested twice for each sample. Four dilution points containing a pool of all samples were used to determine PCR efficiency. Data were analysed using BioMark Gene Expression Data Analysis software (Fluidigm) to obtain Ct values. The Pfaffl method was applied to compute the relative expression of each gene [25]. A linear model was performed to reveal the genes differentially expressed according to the population (high, low) and the time point (0, 1, 4, 24 h):

$$ {expr}_i={\beta}_0+{\beta}_1t+{\beta}_2{population}_i+{\beta}_3{t}^{\ast }{population}_i+{\varepsilon}_i $$

Statistical analysis

A PCA was performed on raw data obtained on 6-week-old piglets using the FactoMineR package in the R software (version 3.2.1, R Core Team, 2016) [26]. Redundant parameters were excluded from the analysis, i.e. MCV, MCHC, % Ret, WBC and lymphocyte counts. The first two components explained more than 40% of the total inertia and each contributed to more than 15% of the total variance.

Univariate statistical analyses for experiment 1 were made using Graph Pad Prism (version 6.0). For each parameter measured, 16 HPAlo and 16 HPAhi pigs were compared using non-parametric Mann-Whitney tests; statistical significance appears as stars on charts (*: p < 0.05; **: p < 0.01; ***: p < 0.001); each point represents a single animal and lines denote medians.

For experiment 2, cortisol, TNFα, IL-6, white blood cells and platelets data were transformed into their logarithmic scores for normalisation. The square root transformation was used for NEFA (non-esterified fatty acids) levels [19]. For each parameter, the effect of HPA activity, time, and their interaction was tested in a two-way analysis of variance with repeated measures for the time factor. The effects of LPS on tympanic temperature, cortisolemia, plasma TNFα and IL-6 levels were graphed, each point representing a single pig.



Adrenocorticotropic hormone


Complete blood cell counts






Glucocorticoid receptors








Immature reticulocyte fraction




Mean corpuscular hemoglobin concentration


Mean corpuscular volume


Non-esterified fatty acids


Principal component analysis


Red blood cell count


Reticulocyte count


Helper T cell


Whole-blood assay


White blood cell count


  1. 1.

    Prunier A, Heinonen M, Quesnel H. High physiological demands in intensively raised pigs: impact on health and welfare. Animal. 2010;4(6):886–98.

    CAS  PubMed  Google Scholar 

  2. 2.

    WHO Global Action Plan on Antimicrobial Resistance; accessed on 18/10/18;

  3. 3.

    United Nations: Draft political declaration of the high-level meeting of the General Assembly on antimicrobial resistance; accessed on 18/10/18;

  4. 4.

    The new EU One Health Action Plan against Antimicrobial Resistance ; accessed on 18/10/18;

  5. 5.

    Mormede P, Foury A, Terenina E, Knap PW. Breeding for robustness: the role of cortisol. Animal. 2011;5(5):651–7.

    CAS  PubMed  Google Scholar 

  6. 6.

    Mormede P, Foury A, Barat P, Corcuff J-B, Terenina E, Marissal-Arvy N, Moisan M-P. Molecular genetics of hypothalamic-pituitary-adrenal axis activity and function. Annals of the New-York Academy of Sciences. 2011;1220:127–36.

    CAS  PubMed  Google Scholar 

  7. 7.

    Mormede P, Terenina E. Molecular genetics of the adrenocortical axis and breeding for robustness. Domest Anim Endocrinol. 2012;43(2):116–31.

    CAS  PubMed  Google Scholar 

  8. 8.

    Larzul C, Terenina E, Foury A, Billon Y, Louveau I, Merlot E, Mormede P. The cortisol response to ACTH in pigs, heritability and influence of corticosteroid-binding globulin. Animal. 2015;9(12):1929–34.

    CAS  PubMed  Google Scholar 

  9. 9.

    Larzul C, Terenina E, Billon Y, Gress L, Mormede P. Genetic determinism of cortisol levels in pigs. In: World Congress on Genetics Applied to Livestock Production, vol. 2. Auckland, New Zealand: Volume species- Porcine; 2018. p. 214.

    Google Scholar 

  10. 10.

    Bacou E, Haurogné K, Mignot G, Allard M, De Beaurepaire L, Marchand J, Terenina E, Billon Y, Jacques J, Bach JM, et al. Acute social stress-induced immunomodulation in pigs high and low responders to ACTH. Physiol Behav. 2017;169:1–8.

    CAS  PubMed  Google Scholar 

  11. 11.

    Piva E, Brugnara C, Spolaore F, Plebani M. Clinical utility of reticulocyte parameters. Clin Lab Med. 2015;35(1):133.

    PubMed  Google Scholar 

  12. 12.

    Buttarello M. Laboratory diagnosis of anemia: are the old and new red cell parameters useful in classification and treatment, how? Int J Lab Hematol. 2016;38:123–32.

    PubMed  Google Scholar 

  13. 13.

    Montero D, Breenfeldt-Andersen A, Oberholzer L, Haider T, Goetze JP, Meinild-Lundby AK, Lundby C. Erythropoiesis with endurance training: dynamics and mechanisms. Am J Phys Regul Integr Comp Phys. 2017;312(6):R894–902.

    Google Scholar 

  14. 14.

    Citartan M, Gopinath SCB, Chen Y, Lakshmipriya T, Tang TH. Monitoring recombinant human erythropoietin abuse among athletes. Biosens Bioelectron. 2015;63:86–98.

    CAS  PubMed  Google Scholar 

  15. 15.

    Martin LB. Stress and immunity in wild vertebrates: timing is everything. Gen Comp Endocrinol. 2009;163(1–2):70–6.

    CAS  PubMed  Google Scholar 

  16. 16.

    Davis AK, Maney DL, Maerz JC. The use of leukocyte profiles to measure stress in vertebrates: a review for ecologists. Funct Ecol. 2008;22(5):760–72.

    Google Scholar 

  17. 17.

    Dhabhar FS, Malarkey WB, Neri E, McEwen BS. Stress-induced redistribution of immune cells-from barracks to boulevards to battlefields: a tale of three hormones - Curt Richter award winner. Psychoneuroendocrinology. 2012;37(9):1345–68.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Engert LC, Weiler U, Pfaffinger B, Stefanski V, Schmucker SS. Diurnal rhythms in peripheral blood immune cell numbers of domestic pigs. Dev Comp Immunol. 2018;79:11–20.

    CAS  PubMed  Google Scholar 

  19. 19.

    Terenina E, Sautron V, Ydier C, Bazovkina D, Sevin-Pujol A, Gress L, Lippi Y, Naylies C, Billon Y, Liaubet L, et al. Time course study of the response to LPS targeting the pig immune gene networks. BMC Genomics. 2017;18:988.

    PubMed  PubMed Central  Google Scholar 

  20. 20.

    Tuchscherer M, Kanitz E, Tuchscherer A, Puppe B. Effects of social support on glucocorticoid sensitivity of lymphocytes in socially deprived piglets. Stress-the International Journal on the Biology of Stress. 2016;19(3):325–32.

    CAS  Google Scholar 

  21. 21.

    Murani E, Reyer H, Ponsuksili S, Fritschka S, Wimmers K. A substitution in the ligand binding domain of the porcine glucocorticoid receptor affects activity of the adrenal gland. PLoS One. 2012;7(9):e45518.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Murani E, Ponsuksili S, Jaeger A, Gorres A, Tuchscherer A, Wimmers K. A naturally hypersensitive glucocorticoid receptor elicits a compensatory reduction of hypothalamus-pituitary-adrenal axis activity early in ontogeny. Open Biol. 2016;6(7):150193.

  23. 23.

    Sautron V, Terenina E, Gress L, Lippi Y, Billon Y, Larzul C, Liaubet L, Villa-Vialaneix N, Mormede P. Time course of the response to ACTH in pig: biological and transcriptomic study. BMC Genomics. 2015;16:961.

    PubMed  PubMed Central  Google Scholar 

  24. 24.

    Spurgeon SL, Jones RC, Ramakrishnan R. High throughput gene expression measurement with real time PCR in a microfluidic dynamic Array. PLoS One. 2008;3(2):e1662.

    PubMed  PubMed Central  Google Scholar 

  25. 25.

    Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29(9):e45.

    CAS  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Le S, Josse J, Husson F. FactoMineR: an R package for multivariate analysis. J Stat Softw. 2008;25(1):1–18.

    Google Scholar 

Download references


We acknowledge the staff of the INRA Experimental Station of “le Magneraud” for pig care, handling and blood sampling.


This work was partly funded by the French National Research Agency (SuSoSTRESS project, #ANR-12-ADAP-0008) to PM. The funding agency had no role in the design and conduct of the study, the collection, management, analysis, and interpretation of the data; or the preparation, review, or approval of the manuscript.

Availability of data and materials

The materials and datasets used and/or analysed during the current study are available from the corresponding authors on reasonable request.

Author information




JH, ET, PM and BL conceived and designed the experiments. JH, ET, KH, EB, MA, EK performed the experiments. YB was in charge of animal management. JH, ET, PM and BL analysed the data. JH, ET, JMB, PM and BL prepared the manuscript. All authors read, commented on and approved the manuscript.

Corresponding authors

Correspondence to Julie Hervé or Elena Terenina.

Ethics declarations

Ethics approval

All experimental procedures were performed in accordance with the EU directive 2010/63 for animal experiments and submitted to and approved by the C2EA 84 Poitou-Charentes Region ethics committee (decision #CE2013–1, 2013.01.21).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S1. PCA variables’ correlogram. Correlations of each variable to the first two dimensions of the PCA are shown. CTL, CD8αhi cytotoxic T lymphocytes. (PDF 129 kb)

Additional file 2:

Table S1. Variable’s contributions to the first two dimensions of the PCA. Variables that contributed to more than 10% to one axis are underlined. (PDF 35 kb)

Additional file 3:

Figure S2. Gating strategy for non B-lymphocyte subsets analysis. Representative flow cytometry profile is shown. Nucleated single cells were identified as DRAQ5-positive cells. Among SSC-Alo cells, CD4/CD8α/CD3 co-staining allowed the determination of 6 lymphocyte subsets as mentioned on the plots. (PDF 935 kb)

Additional file 4:

Table S2. List of genes whose expression was significantly modulated after LPS injection compared to their t0 level. (PDF 43 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Hervé, J., Terenina, E., Haurogné, K. et al. Effects of divergent selection upon adrenocortical activity on immune traits in pig. BMC Vet Res 15, 71 (2019).

Download citation


  • Stress
  • Hypothalamic-pituitary-adrenocortical axis
  • Cortisol
  • LPS challenge
  • Genetics
  • Robustness
  • Immunity