# Hair cortisol concentration in finishing pigs on commercial farms: variability between pigs, batches, and farms

Pierre Levallois, Mily Leblanc-Maridor, Anne Lehébel, Solenn Gavaud, Blandine Lieubeau, Julie Hervé, Christine Fourichon, Catherine Belloc

► To cite this version:

Pierre Levallois, Mily Leblanc-Maridor, Anne Lehébel, Solenn Gavaud, Blandine Lieubeau, et al. Hair cortisol concentration in finishing pigs on commercial farms: variability between pigs, batches, and farms. Frontiers in Veterinary Science, 2024, 10, pp.1298756. 10.3389/fvets.2023.1298756. hal-04410821

HAL Id: hal-04410821

https://oniris.hal.science/hal-04410821v1

Submitted on 9 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN ACCESS

EDITED BY
Flaviana Gottardo,
University of Padua, Italy

REVIEWED BY
Anna Valros,
University of Helsinki, Finland
Ulrike Gimsa,
Research Institute for Farm Animal Biology,
Germany

*CORRESPONDENCE

Catherine Belloc

catherine.belloc@oniris-nantes.fr

RECEIVED 22 September 2023

ACCEPTED 26 December 2023

PUBLISHED 22 January 2024

CITATION

Levallois P, Leblanc-Maridor M, Lehébel A,
Gavaud S, Lieubeau B, Hervé J,
Fourichon C and Belloc C (2024) Hair cortisol
concentration in finishing pigs on commercial
farms: variability between pigs, batches, and
farms.

Front. Vet. Sci. 10:1298756.
doi: 10.3389/fvets.2023.1298756

COPYRIGHT

© 2024 Levallois, Leblanc-Maridor, Lehébel, Gavaud, Lieubeau, Hervé, Fourichon and Belloc. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

# Hair cortisol concentration in finishing pigs on commercial farms: variability between pigs, batches, and farms

Pierre Levallois¹, Mily Leblanc-Maridor¹, Anne Lehébel¹, Solenn Gavaud², Blandine Lieubeau², Julie Hervé², Christine Fourichon¹ and Catherine Belloc¹*

¹Oniris, INRAE, BIOEPAR, Nantes, France, ²Oniris, INRAE, IECM, Nantes, France

Hair cortisol is a stress indicator and could be used to assess the pigs' exposure to stressors in the weeks/months prior to non-invasive hair sampling. The main aim of this study was to describe the hair cortisol concentration (HCC) variability between individuals within a batch, between farms and between batches within a farm. The secondary aim was to determine how the number of sampled pigs influences the characterization of HCC within a batch. Twenty farrow-to-finish pig farms were recruited considering the diversity of their management practices and health status (data collected). Hair was sampled in two separate batches, 8 months apart. The necks of 24 finishing pigs were clipped per batch the week prior to slaughter. To describe the variability in HCC, an analysis of the variance model was run with three explanatory variables (batch, farm and their interaction). To identify farm clusters, a principal component analysis followed by a hierarchical clustering was carried out with four active variables (means and standard deviations of the two batches per farm) and 17 supplementary variables (management practices, herd health data). We determined how the number of sampled pigs influenced the characterization of HCC within a batch by selecting subsamples of the results. HCC ranged from 0.4 to 121.6 pg/mg, with a mean of 25.9 ± 16.2 pg/mg. The variability in HCC was mainly explained by differences between pigs (57%), then between farms (24%), between batches within the same farm (16%) and between batches (3%). Three clusters of farms were identified: low homogeneous concentrations (n = 3 farms), heterogeneous concentrations with either higher (n = 7) or lower (n = 10) HCC in batch 2 than in batch 1. The diversity of management practices and health statuses allowed to discuss hypotheses explaining the HCC variations observed. We highlighted the need to sample more than 24 pigs to characterize HCC in a pig batch. HCC differences between batches on six farms suggest sampling pigs in more than one batch to describe the HCC at the farm level. HCC variations described here confirm the need to study its links with exposure of pigs to stressors.

KEYWORDS

hair cortisol, stress, biomarker, pig, welfare, health, field study
Introduction

Minimizing distress is essential to ensure the welfare of pigs on farms (1, 2). Distress is an aversive, negative state in which coping and adaptation processes fail to return an organism to physiological homeostasis (3). Distress is caused by the exposure of an individual to stressors. Stressors can be physiological or psychological threats, either real or perceived (4, 5).

On pig farms, stressors can be multiple and occur at various stages of life. It can be assumed that the exposure of pigs to stressors varies from farm to farm and from batch to batch within the same farm. For instance, the presence or absence of castration, tail docking or ear tagging from one farm to another could lead to different exposure of piglets to stressors (6).

The intensity of the stress response can vary from one animal to another when facing with the same stressor, depending in particular on their past experiences (7). Cortisol is a widely used indicator to assess the stress response in most mammals (8). Cortisol can be quantified in various matrices, including blood, saliva, urine, feces, and more recently in hair (8–10). Depending on the matrix, quantified cortisol can provide information about stress levels in different time windows prior to sampling. Serum or salivary cortisol provides information on stress levels at a specific moment, corresponding to a few minutes before sampling. Urinary or fecal cortisol (or its metabolites) reflects the stress level experienced during a period of accumulation, corresponding to several hours until urination or defecation occurs (10). Since cortisol accumulates during hair growth and remains a priori stored in the hair for several weeks or even months, hair cortisol is considered as a biomarker of chronic stress in humans (11, 12). In pigs, hair cortisol could potentially be used to assess stress throughout a long time period prior to sampling on commercial farms. Consequently, hair cortisol could possibly be a measure of on-farm pig welfare assessment. However, this possibility still needs to be assessed.

Measuring hair cortisol offers several advantages. Hair sampling is non-invasive, painless, and often does not require the animal to be restrained. Animal handling during hair sampling does not distort the measured concentration, as demonstrated by studies on dairy cows and reindeer (13, 14). Cortisol secreted during handling is unlikely to be present in hair samples, unlike blood samples. Hair samples can be stored at room temperature and protected from light for at least 1 year without significantly altering cortisol concentration, as studies on dairy cows, bears, and chimpanzees have shown (13, 15, 16).

Little is known about the mechanism of incorporation of cortisol into hair. Therefore, the definition of the time window considered by hair cortisol is uncertain (a few weeks or several months) (12). Cortisol could incorporate into hair through passive diffusion from blood (17). In that case, cortisol would incorporate into hair only during hair growth when hair follicle is linked to blood capillaries. Knowledge about the pig hair cycle on commercial farm is scarce. Watson and Moore (18) showed that hair growth is continuous but varies through time in piglets. Evidences showed that cortisol seems to accumulate in hair of growing pigs both before weaning and between weaning and slaughter (19). Therefore, hair cortisol measured at the end of the finishing period could potentially provide information on the exposure of pigs to stressors throughout their life on farm.

To our knowledge, studies providing results on hair cortisol in finishing pigs on commercial farms were carried out on a single farm (19–21). There is a lack of data to determine how hair cortisol concentration at the end of the finishing period varies in different non-experimental farming contexts. In absence of variation at the end of the finishing period, the hair cortisol concentration would provide only limited information regarding long-term stress. Therefore, it is important to describe the variability of hair cortisol concentration in finishing pigs before carrying out further research. It is assumed that hair cortisol concentration of finishing pigs may vary between farms and between batches within the same farm due to different exposure to stressors throughout their lives.

The number of animals sampled to describe hair cortisol concentration within a batch must consider the interindividual variability in stress response. Although little is known about interindividual variability in the literature, the study by Gavaud et al. (19) reported hair cortisol concentrations ranging from 30.0 to 177.6 pg/mg in 67 21-week-old pigs from three batches on the same conventional farm. The sample size must therefore be sufficient to take account of these interindividual differences.

Nevertheless, the feasibility of a hair cortisol measurement protocol is a limit to the determination of sample size, particularly when sampling is carried out on several farms. Total sample size in previous studies of hair cortisol at the end of the finishing period has ranged from 42 to 97 pigs, distributed over one to three batches within a single farm (19–23). Sampling several dozens of pigs per batch across several batches and farms may not be feasible given the constraints of laboratory protocol (e.g., labor time and cost of analysis). Consequently, to characterize the variability in hair cortisol concentration in finishing pigs, the sampling design must strike a balance between the number of farms and the number of pigs sampled per farm.

The aim of this exploratory study was to describe the variability in hair cortisol concentration in finishing pigs. Firstly, we wanted to characterize interindividual differences within a given batch. Secondly, we wanted to identify potential differences between farms (reflecting different levels of exposure to stressors). Thirdly, we wanted to identify potential differences between batches on the same farm. Finally, a secondary objective was to determine how the number of pigs sampled influences the characterization of hair cortisol concentration within a batch.

Methods

Description and recruitment of farms

This study was carried out as part of the European research project HealthyLivestock. A cohort of 20 farrow-to-finish pig farms in western France was monitored longitudinally. Farms were recruited based on criteria relating to the diversity in management practices and health status. This diversity was assumed to introduce variability in pig exposure to stressors across farms (in terms of nature, number, and intensity). The farms housed from 70 to 800 sows, with an average of 244 sows. Within the cohort, two farms were organic and 18 were conventional. Seven farms out of the 18 conventional had a specification: four farms were Label Rouge [product quality specification with minimal age at slaughter of 182 days and conditions
for husbandry including lower animal densities, for details see (24)), two farms were antibiotic-free from birth, one farm was antibiotic-free from 42-day-old. Data were collected during two visits. During the initial visit, tailor-made health plans were formulated by veterinarians. The second visit was carried out on average 8 months later to assess the herd health evolution (25).

## Selection of sampled pigs

Pigs were sampled at the end of the finishing period at an average age of 165 ± 10 days. Sampling was carried out in two batches (one batch per visit), 8 months apart. Therefore, sampled pigs were different animals in the two batches on the same farm. The first sampling was performed in winter/early spring then the second one in autumn/early winter. Twenty-four pigs were sampled in each batch. Sample size was determined by reconciling the need for a sufficient sample size to describe the interindividual variability in hair cortisol concentration within a batch and the feasibility of the protocol. Sampled pigs were randomly selected within three or four pens, except on two farms where pigs were sampled in, respectively, one and two pens due to housing facilities. Pens were randomly selected on farm with more than four pens. Sex and breed of sampled pigs were not recorded. However, no variation of the hair cortisol concentration should be induced due to a sex effect for pigs aged of 165 ± 10 days (26). Most farms housed only castrated boars and females (Supplementary Table S1). Regarding breeds, all sows were Large White x Landrace and boars were predominantly Pietrain or sometimes Duroc. Our sample did not allow us to test a putative breed effect.

## Hair sampling procedure

Sampling was carried out without restraint by clipping the hair as close as possible to the skin of the dorsal area of the neck region, to sample the total length of the hair. Moreover, this method avoided to sample hair follicles likely to contain cortisol (via blood or their potential endocrine activity) (23, 27). The sampling region remained unchanged to avoid potential concentration variations between body regions (15, 16, 23). The neck was chosen as the sampling region to minimize potential contamination of hair with urinary or fecal cortisol (28, 29). At least 40 mg of hair were sampled per pig. The samples were placed in envelopes and stored in the dark at room temperature for a period ranging from 2 to 10 months until the measurements of hair cortisol concentration.

## Measurement of hair cortisol concentration

The protocol used to measure hair cortisol concentration was based on Gavaud et al. (19), adapted from Heimbürge et al. (26). Each hair sample was placed into a 15 mL polypropylene tube for two wash cycles. The first cycle consisted of two 3-min washes with water (5 mL of water per wash) to remove dust particles. The second cycle involved two 3-min washes with isopropanol (HPLC grade, Sigma-Aldrich) to remove potential lipophilic contaminants. Hair samples were then placed in Petri

dishes and dried in the dark at room temperature for at least 48 h. Dried hairs were finely cut into 2–3 mm portions, and around 35 mg to 40 mg of each sample were weighed and placed in microtubes. To grind the hair samples, three stainless steel grinding balls were added to each tube, and the samples were briefly frozen by immersing the tubes in liquid nitrogen. The microtubes were then placed directly into a Retsch ball mill (MM 400, Verder Scientific, Eragny-sur-Oise, France) for two consecutive grinding cycles at 30 Hz for 10 min. A brief centrifugation was performed to gather the hair powder at the bottom of the tubes. To solubilize cortisol from the hair samples, 1.5 mL of methanol (HPLC grade, Sigma-Aldrich) was added to each tube. Extraction of cortisol from hair samples was carried out at room temperature for 24 h, with the tubes placed on a plate shaker set at 800 rpm. Then, the tubes were centrifuged at 1,500×g for 5 min, and a volume of 750 µL of supernatant was taken from each tube and transferred to new tubes. The new tubes were placed in a vacuum concentrator (5,303 vacuum concentrator, Eppendorf, Montesson, France) at 45°C for approximately 1 h to evaporate the methanol. Finally, tubes containing cortisol in solid form were stored at –20°C until further analysis. Cortisol was re-solubilized in 140 µL of PBS solution and measured in duplicate using an ELISA kit for free cortisol in saliva (Demeditec Diagnostics GmbH, Kiel, Germany) following the manufacturer's instructions. Concentrations were converted from ng/mL to pg/mg of hair for data analysis. The detection limit of the ELISA kit was 0.019 ng/mL (i.e., 0.13–0.15 pg/mg for hair samples weighing between 35 and 40 mg).

## Recording of management practices and herd health data

Management practices and herd health data were recorded at the first visit. The same recording was repeated 8 months later, at the second visit, to identify any changes in management practices or the occurrence of health events between the two visits. Data requested from farmers included: (1) the procedures applied to piglets in the farrowing unit (castration, tail docking, teeth grinding, iron administration), (2) the age at which weaning and other social regrouping took place (i.e., when pigs were moved from one pen to another, resulting in a new social structure and potentially the establishment of a new social hierarchy) and the criteria used to distribute pigs in the new pens (i.e., litter, sow parity, weight or pig sex), (3) the minimal floor space allowance per pig in the post-weaning and finishing units. Housing facilities were described at each visit, indicating whether they were indoors, indoors with outdoors access, or outdoors. With regard to herd health, the herd veterinarian was asked to provide information on: (1) any health disorder requiring the formulation of a health plan on the farms recruited, (2) the physiological stages affected by these health disorders, and (3) whether there had been any change in herd health between the first and second visits. Change in herd health statuses was identified by veterinarians through clinical observations of animals. All data are displayed in Supplementary Tables S1, S2.

## Data analyses

Statistical analyses were carried out using R version 4.0.5. No transformation or adjustment were made to the raw data. No outlier was identified in the raw data.
Firstly, the distribution of cortisol concentration was described for the entire population of sampled pigs and separately for each farm's batch.

Secondly, to characterize the variability in hair cortisol concentration based on farm and batch, an analysis of variance model was implemented using the aov function (stats package). The model was validated by an analysis of residuals. Normality was assessed using a graphical analysis (QQ-plot) and a Shapiro-Wilk test, while homoscedasticity was assessed using a graphical examination of residual linearity and a Levene test. The analysis of variance model included the effects of batch, farm, and their interaction. Variables were considered significant if p-values were < 0.05 (anova test). To determine the proportion of variance explained by each predictor variable (batch, farm and their interaction), the effect size η² was calculated (ratio of the sum of squares of a predictor variable to the total sum of squares of the model). The η² indicator describes the proportion of concentration variability explained by each variable, after adjustment for other predictor variables. A pairwise post-hoc test (Tukey's test) was performed to compare mean hair cortisol concentrations between batches from the same farm. The absolute value of Cohen's d was calculated for comparisons of means between the two batches on all farms (ratio between the difference of the means of compared batches and the standard deviation) (30). Cohen's d is a complementary indicator to the p-value, providing an estimate of the magnitude of the difference between two means.

Thirdly, to determine the presence of farm profiles, a Principal Component Analysis (PCA) followed by a Hierarchical Clustering on Principal Components (HCPC) was carried out. Farms were distinguished using four active variables (i.e., variables contributing to the distinction of farms). The active variables represented the means and standard deviations of hair cortisol concentration obtained per farm, considering both batches. Seventeen supplementary variables were considered (i.e., they do not contribute to the distinction of farms but provide information on the characteristics of the farm clusters identified). The supplementary variables included management practices and pig health data, as reported by farmers and veterinarians. They included: (1) castration, (2) tail docking, (3) teeth grinding, (4) iron administration, age at potential social regroupings, (5) during the suckling period (i.e., early socialization), (6) at weaning, (7) during fattening, and their (8–10) respective modalities (i.e., regrouping based on litter, sow parity, weight or pig sex), type of housing facilities (i.e., indoors, indoors with outdoors access, outdoors) during the (11) suckling, (12) post-weaning and (13) fattening periods, the minimal floor space allowance per pig in the (14) post-weaning and (15) fattening units, (16–17) the presence of a health disorder for each batch.

Finally, a sampling procedure was implemented to determine how the mean and standard deviation of hair cortisol concentration in a batch varied with decreasing sample size. For all batches, all the possible samples of 5, 10, 15, and 20 pigs were drawn from the initial sample of 24 animals. Two calculations were performed for each sample: (1) the relative difference between the means of the drawn samples and the initial sample (i.e., with the 24 sampled animals; Equation 1), and (2) the relative difference between the standard deviations of the drawn samples and the initial sample (Equation 2). The distributions of the relative differences in means and standard

deviations were described. Summary statistics (5th, 50th, and 95th percentiles) were computed for the relative differences in means and standard deviations for all sample sizes (5, 10, 15, and 20 pigs) in all batches studied (40 batches). Forty percentile values were obtained for each sample size. The means of the summary statistics (5th, 50th, and 95th percentiles) were computed for each sample size (5, 10, 15, and 20 pigs) in all batches.

$$ \text{Relative difference between means (\%): } \Delta_{i,j,n} = \frac{\mu_{i,j} - \mu_n}{\mu_n} \quad (1) $$

Relative difference between standard deviations (%):

$$ \Delta'_{i,j,n} = \frac{\sigma_{i,j} - \sigma_n}{\sigma_n} \quad (2) $$

Where μ=mean; σ=standard deviation; i=number of the drawn sample, i ∈ [1; N]; j=size of the drawn sample, j={5; 10; 15; 20}; n=identification number of a given farm, n ∈ [1; 20].

# Results

## Interindividual variability in hair cortisol concentration at the sampled population level

The distribution of hair cortisol concentration of 950 pigs raised on 20 farms (two batches per farm) is shown in Figure 1. Of note, 10 pigs from five batches on five farms were not sampled due to premature departure for slaughter. The distribution appears to be asymmetrical. Hair cortisol concentrations ranged from 0.4 to 121.6 pg/mg, with 50% of the values <25 pg/mg. The mean hair cortisol concentration was 25.9 pg/mg (standard deviation = 16.2 pg/mg).

## Interindividual variability in hair cortisol concentration within a batch

The hair cortisol concentrations in pigs sampled in the two batches 8 months apart are shown farm by farm in Figure 2. The coefficients of variation for hair cortisol concentration within a batch ranged from 0.17 to 0.71. The standard deviation of hair cortisol concentration within a batch increased as the mean hair cortisol concentration of the batch increased (correlation coefficient: 0.83; p < 0.001; Figure 3). No batch had a high mean and a low standard deviation of hair cortisol concentration. For mean values ranging from 20 to 50 pg/mg, some batches showed higher homogeneity while others displayed higher heterogeneity (coefficient of variation ranging from 0.18 to 0.64).

The proportion of variability (η²) explained by the three predictor variables in the analysis of the variance model (batch, farm and their interaction) is presented in Table 1. The residuals explained 57% of the total variability in hair cortisol concentration (Table 1). Residuals represent the portion of the variability not explained by the predictor variables in the analysis of the variance model (batch, farm and their interaction), thus corresponding to interindividual differences.
<figure><img src="image_3.png" /><figcaption>FIGURE 1
Distribution of hair cortisol concentrations in 950 finishing pigs sampled in two batches per farm on 20 farms.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>FIGURE 2
Hair cortisol concentration in finishing pigs sampled in two batches per farm (1 and 2) on 20 farms (A to T). Each dot corresponds to the concentration of one pig and the red lines indicate the median concentration of a batch.</figcaption></figure>

## Variability in hair cortisol concentration between farms

The farm effect explained 24% of the total variability in hair cortisol concentration (without considering its interaction with the batch effect, Table 1). Farms A, B, and G had the lowest hair cortisol concentration (Table 2).

## Variability in hair cortisol concentration between batches

The mean hair cortisol concentration ranged from 4.8 to 54.5 pg/mg across batches. The batch effect, without considering its interaction with the farm effect, explained the lowest proportion of the variability in hair cortisol concentration (η²=3%, Table 1). The batch:farm
<figure><img src="image_5.png" /><figcaption>FIGURE 3
Linear regression of standard deviation against mean hair cortisol concentration in finishing pigs sampled in two batches per farm on 20 farms (one dot corresponds to one batch).</figcaption></figure>

<table><caption>TABLE 1 Proportion of the variability in hair cortisol concentration of finishing pigs sampled from two batches per farm on 20 farms, as a function of the predictor variables included in the analysis of variance model.</caption><thead><tr><td>Predictor variable</td><td>Percentage of variance explained (η²)</td></tr></thead><tbody><tr><td>Batch</td><td>3.1</td></tr><tr><td>Farm</td><td>24.1</td></tr><tr><td>Batch: Farm</td><td>16.3</td></tr><tr><td>Residuals</td><td>56.5</td></tr><tr><td>Total</td><td>100</td></tr></tbody></table>

interaction effect explained 16% of the total variability in hair cortisol concentration (Table 1).

Within a given farm, the mean hair cortisol concentration in batch 1 could be significantly lower (on 4 farms), higher (on 2 farms), or not significantly different (on 14 farms) from that in batch 2 (Table 2). All Cohen's *d* values were > 1 on farms where the mean hair cortisol concentration differed significantly between the two batches, except for farm O (indicating a difference of more than one standard deviation between the two compared means). Farm A was the only one with a Cohen's *d* > 1.0 and a non-significant difference in mean concentrations between the two batches (*p* > 0.05, Tukey's test).

On the two farms where the mean hair cortisol concentration in batch 1 was significantly higher than that in batch 2, veterinarians reported an improvement in health (farms J, Q). Also, castration was ceased on farm Q, while no change in management practices was noted on farm J between the two samplings. Among the four farms where the mean hair cortisol concentration in batch 1 was significantly lower than that in batch 2, veterinarians observed: an improvement in health between the two batches sampled 8 months apart on farm D; the occurrence of a

health disorder in the second batch on farm O; the presence of an unchanged severity of a health disorder on farm K; and no health disorder on farm E. No change in management practices was reported on these four farms between the two batches (farms D, E, K, O).

## Clustering of farms based on hair cortisol concentration

The first two principal dimensions of the PCA explained 91.5% of the variance (Supplementary Figure S1). HCPC allowed the identification of three clusters of farms, denoted as C1 (three farms), C2 (seven farms), and C3 (10 farms). The *v.test* value is provided in Table 3 when *P* was <0.05.

### Cluster 1 (A, B, G)

This cluster was defined by lower means of hair cortisol concentration than those of the cohort for both batches, and a lower intra-batch variability for batch 2 only. This cluster was characterized by seven supplementary variables, with the following modalities (Table 3): no tail docking, early socialization in the farrowing unit, weaning age of 42-day-old, post-weaning and finishing housing facilities providing outdoors access, a minimal floor space allowance per pig, respectively, not <0.3 m² and > 1.5 m² in the post-weaning and the finishing pens. These characteristics corresponded to farms A and B, which were the only two organic farms in the cohort. Farm G, a conventional farm without specification, had very different management practices. Indeed, pigs on farm G were tail docked, not early socialized in the farrowing unit, weaned at 28-day-old, raised in post-weaning and finishing facilities without outdoors access, with a minimal floor space allowance per pig <0.3 m²/pig in post-weaning and 0.7–0.9 m²/pig in finishing (Supplementary Tables S1, S2).

### Cluster 2 (D, E, F, H, I, K, O)

This cluster was defined by higher means and intra-batch variabilities of hair cortisol concentration in batch 2 compared to the cohort. The mean hair cortisol concentration was significantly different between the two batches on four farms out of seven (farms D, E, K, O). This cluster was characterized by one supplementary variable with the following modality (Table 3): the composition of pig pens in the finishing stage was based on the weight of pigs at the end of the post-weaning period, resulting in new groups of pigs and novel social interactions among pigs.

### Cluster 3 (C, J, L, M, N, P, Q, R, S, T)

This cluster was defined by higher means and intra-batch variabilities of hair cortisol concentration in batch 1 compared to the cohort. The mean hair cortisol concentration was similar and high for the two batches in this cluster, except for farms J and Q. No supplementary variable was associated with this cluster.

## Sample size and variability in hair cortisol concentration

The percentage of pigs sampled per batch ranged from about 5 to 17% of the total number of finishing pigs, with a median
<table><caption>TABLE 2 Hair cortisol concentrations (pg/mg) in two batches of finishing pigs sampled 8 months apart on 20 farms (24 pigs were sampled per batch regardless of the total number of pigs in the batch).</caption><thead><tr><th rowspan="2">Farm</th><th colspan="2">First batch</th><th colspan="2">Second batch</th><th rowspan="2">Cohen's d<sup>1</sup><br>(absolute value)</th><th rowspan="2">p-value<br>(Tukey's test)</th></tr><tr><th>Mean (μ1)</th><th>Standard<br>deviation</th><th>Mean (μ2)</th><th>Standard<br>deviation</th></tr></thead><tbody><tr><td>A</td><td>4.9</td><td>1.6</td><td>13.7</td><td>6.1</td><td>1.4</td><td>ns</td></tr><tr><td>B</td><td>4.8</td><td>3.1</td><td>6.5</td><td>1.7</td><td>0.7</td><td>ns</td></tr><tr><td>C</td><td>35.9</td><td>25.7</td><td>42.5</td><td>15.2</td><td>0.3</td><td>ns</td></tr><tr><td>D</td><td>13.6</td><td>8.3</td><td>34.0</td><td>18.0</td><td>1.2</td><td>&lt;0.001</td></tr><tr><td>E</td><td>13.7</td><td>3.7</td><td>54.5</td><td>24.1</td><td>1.5</td><td>&lt;0.001</td></tr><tr><td>F</td><td>24.5</td><td>6.0</td><td>32.0</td><td>13.2</td><td>0.7</td><td>ns</td></tr><tr><td>G</td><td>15.6</td><td>4.7</td><td>13.8</td><td>3.7</td><td>0.4</td><td>ns</td></tr><tr><td>H</td><td>17.6</td><td>5.5</td><td>25.0</td><td>14.3</td><td>0.7</td><td>ns</td></tr><tr><td>I</td><td>20.3</td><td>3.6</td><td>32.6</td><td>21.1</td><td>0.8</td><td>ns</td></tr><tr><td>J</td><td>36.5</td><td>9.7</td><td>18.4</td><td>7.5</td><td>1.4</td><td>&lt;0.001</td></tr><tr><td>K</td><td>11.3</td><td>3.6</td><td>29.6</td><td>11.8</td><td>1.4</td><td>&lt;0.001</td></tr><tr><td>L</td><td>30.2</td><td>9.1</td><td>26.5</td><td>8.8</td><td>0.4</td><td>ns</td></tr><tr><td>M</td><td>29.7</td><td>7.1</td><td>24.6</td><td>10.4</td><td>0.6</td><td>ns</td></tr><tr><td>N</td><td>26.9</td><td>12.6</td><td>28.6</td><td>6.0</td><td>0.2</td><td>ns</td></tr><tr><td>O</td><td>29.5</td><td>8.8</td><td>43.8</td><td>24.5</td><td>0.7</td><td>&lt;0.05</td></tr><tr><td>P</td><td>31.8</td><td>13.7</td><td>29.2</td><td>15.1</td><td>0.2</td><td>ns</td></tr><tr><td>Q</td><td>37.1</td><td>21.2</td><td>16.6</td><td>4.8</td><td>1.1</td><td>&lt;0.001</td></tr><tr><td>R</td><td>21.3</td><td>12.4</td><td>29.4</td><td>9.9</td><td>0.7</td><td>ns</td></tr><tr><td>S</td><td>27.2</td><td>11.8</td><td>37.4</td><td>15.0</td><td>0.7</td><td>ns</td></tr><tr><td>T</td><td>25.0</td><td>9.0</td><td>33.9</td><td>12.9</td><td>0.6</td><td>ns</td></tr></tbody></table>

<table><caption>TABLE 3 Characterization of farm clusters identified by hierarchical clustering on the four active variables and the 17 supplementary variables.</caption><thead><tr><th rowspan="2"></th><th colspan="3">Cluster</th></tr><tr><th>C1</th><th>C2</th><th>C3</th></tr></thead><tbody><tr><td>Number of farms in the cluster</td><td>3</td><td>7</td><td>10</td></tr><tr><td>Value of <i>v.test</i> of active variables</td><td></td><td></td><td></td></tr><tr><td>Mean of the hair cortisol concentration in pigs in batch 1</td><td>-2.73</td><td></td><td>3.29</td></tr><tr><td>Mean of the hair cortisol concentration in pigs in batch 2</td><td>-2.85</td><td>2.10</td><td></td></tr><tr><td>Standard deviation of the hair cortisol concentration in pigs in batch 1</td><td></td><td></td><td>3.07</td></tr><tr><td>Standard deviation of the hair cortisol concentration in pigs in batch 2</td><td>-2.43</td><td>3.02</td><td></td></tr><tr><td>Value of <i>v.test</i> of supplementary variables</td><td></td><td></td><td></td></tr><tr><td>Tail docking (no)</td><td>2.41</td><td></td><td></td></tr><tr><td>Early socialization in the suckling period (no)</td><td>-1.99</td><td></td><td></td></tr><tr><td>Age at weaning (42-day-old)</td><td>2.41</td><td></td><td></td></tr><tr><td>Regrouping in fattening units (weight)</td><td></td><td>2.68</td><td></td></tr><tr><td>Housing facilities during the post-weaning period (closed with outdoor facilities)</td><td>2.41</td><td></td><td></td></tr><tr><td>Housing facilities during the fattening period (closed with outdoor facilities)</td><td>2.41</td><td></td><td></td></tr><tr><td>Minimal floor space allowance per pig during post-weaning period (&lt;0.3 m<sup>2</sup>/porc)</td><td>-1.99</td><td></td><td></td></tr><tr><td>Minimal floor space allowance per pig during fattening period (&gt;1.5 m<sup>2</sup>/porc)</td><td>2.41</td><td></td><td></td></tr></tbody></table>

The *v.test* value is indicated when *P* < 0.05. A positive (or negative) *v.test* value for quantitative variables indicates a higher (or lower) mean or standard deviation in the cluster compared to the cohort. A positive (or negative) *v.test* value for qualitative variables indicates overrepresentation (or underrepresentation) of a modality in the cluster compared to the cohort. The higher the absolute value of the *v.test*, the more the group is characterized by this variable (or modality). Only supplementary variables for which at least one modality had a *p* < 0.05 are shown in the table (modality indicated in brackets). C1, Group of farms with lower means of hair cortisol concentration for both batches and lower interindividual differences only for batch 2, compared to the cohort. C2, Group of farms with higher mean and interindividual differences of hair cortisol concentration in batch 2, compared to the cohort. C3, Group of farms with higher mean and interindividual differences of hair cortisol concentration in batch 1, compared to the cohort.

<footnote n="1">¹Cohen's *d* is an indicator to estimate the magnitude of the difference between two means (expressed as a ratio between the difference of the two means and the standard deviation).</footnote>
<figure><img src="image_6.png" /><figcaption>**FIGURE 4**
Linear regression of the coefficient of variation of hair cortisol concentrations against the percentage of finishing pigs sampled in two batches per farm on 20 farms (one dot corresponds to one batch).</figcaption></figure>

percentage around 8%. The coefficient of variation of hair cortisol concentration per batch was not correlated with the percentage of pigs sampled per batch (correlation coefficient: 0.15; *p* = 0.33; **Figure 4**). For example, when considering a similar proportion of pigs sampled per batch (about 8%), the coefficient of variation ranged from 0.18 to 0.71. The farms with the lowest and highest proportion of pigs sampled per batch in the cohort showed a high and a very low intra-batch variability, respectively. The farm with the lowest percentage of pigs sampled per batch in the cohort (farm H, 5%) had a standard deviation of hair cortisol concentration higher than the third quartile in batch 2 compared to the cohort (14.3 pg/mg). The farm with the highest percentage of pigs sampled per batch in the cohort (farm A, 17%) had the lowest standard deviation of hair cortisol concentration (1.6 pg/mg, batch 1).

**Figure 5** shows the relative differences from the initial means (Figure 5A) and initial standard deviations (Figure 5B), according to the size of the drawn samples (5, 10, 15, or 20 pigs from the initial batch of 24 pigs). The relative differences were lower in means than in standard deviations, all sizes of the random samples considered. The

<figure><img src="image_7.png" /><figcaption>**FIGURE 5**
Percentiles (5th, 50th, 95th) of relative differences between the means (A) or standard deviations (B) in hair cortisol concentration in the initial sample and drawn sample (DS) of n pigs (DS5: n=5 pigs; DS10: n=10 pigs; DS15: n=15 pigs; DS20: n=20 pigs). In all batches of finishing pigs (two batches per farm on 20 farms), all possible samples of 5, 10, 15, and 20 pigs were drawn from the initial sample of 24 pigs. For each drawn sample, relative differences were calculated between the means and standard deviations of the initial sample and the random sample. For each batch (n=40), summary statistics (5th, 50th, 95th percentiles) of these relative differences were computed for all the drawn samples (one dot corresponds to the percentile value of one batch, for a given sample size). The mean of a summary statistic for one size of drawn sample is indicated by a red line.</figcaption></figure>
relative difference in standard deviations was higher in batches where interindividual differences in hair cortisol concentration were large, even with drawn samples of 20 pigs.

## Discussion

This exploratory study is to our knowledge the first one to describe the variability in hair cortisol concentration in finishing pigs on several commercial farms. We aimed to characterize differences between pigs within a batch, between farms (which could reflect different levels of exposure to stressors), and between batches within the same farm. The secondary aim was to determine how the number of sampled pigs influences the characterization of hair cortisol concentration within a batch.

## Variability in hair cortisol concentration explained by interindividual differences in stress response

Interindividual differences explained more than half of the variability in hair cortisol concentration. The range of reported hair cortisol concentration for a total of 421 finishing pigs in six previous studies (ranging from 4.5 to 177.6 pg/mg) is similar to the one reported in this study (with 2.3 times more sampled pigs; ranging from 0.4 to 121.6 pg/mg), except for the lowest values (19–23, 26). The lowest hair cortisol concentration measured in this study had not yet been described in finishing pigs but had been reported in sows (31). Of note, no study comparing inter-laboratory assay methods to measure hair cortisol concentration was, to our knowledge, published. To validate the comparability between hair cortisol concentration measured on different laboratories, further comparison should be carried out.

Differences between pigs of a same batch in hair cortisol concentration could potentially be explained by different stress responses to the exposures of same stressors. For instance, Casal et al. (23) demonstrated that hair cortisol concentration in 21-week-old pigs, subjected to repeated regroupings once a week for the 4 weeks preceding sampling, could differ among individuals (with minimum and maximum concentrations of 6.4 and 40.0 pg/mg). Differences in hair cortisol concentration between pigs after repeated regroupings could be explained by differences in their individual characteristics, such as differences in their dominance ranks (32). Nonetheless, this study was not designed to describe different stress response between pigs for the same exposure to a stressor. Reasons explaining the concentration variability between finishing pigs within the same batch should be explored to better understand the information provided by hair cortisol at the pig level.

## Variability in hair cortisol concentration explained by differences between farms

Differences between farms explained nearly a quarter of the variability in hair cortisol concentration. To our knowledge, this is the

first description of differences in hair cortisol concentration at the end of the finishing period between commercial farms. It is likely that the diverse management practices and health disorders in recruited pig farms contribute to the differences between farms. These farm-level differences suggest the need for further investigations into the factors influencing hair cortisol concentration throughout the lifetime of pigs in farming systems.

The clustering revealed that a minority of farms had low means and differences between individuals in hair cortisol concentration whatever the batch. The two organic farms of the cohort were included in this cluster. These results are interesting because organic management practices, such as early socialization, weaning at 6 weeks of age, and outdoors access, are generally perceived to provide better animal welfare than conventional management practices (33, 34). Moreover, all growing pigs within batches of these two organic farms were housed in one or two pens since weaning. Social stress in these farms could be limited since social hierarchy was early established and not further modified. However, it is not possible to conclude from this study that organic management practices are less stressful for pigs than the other types of management practices. Indeed, only two organic farms were monitored. Moreover, the study design was not suitable for assessing the association between hair cortisol concentration and management practices. Furthermore, the other farm included in the same cluster than the two organic farms was conventional and without specification. We can only suppose that the pigs having the lower hair cortisol concentrations were less exposed to stressors than the others. We hypothesize that a low exposure to stressors limits the observations of large interindividual differences within a batch.

## Variability in hair cortisol concentration explained by differences between batches

This study is, to our knowledge, the first to evidence differences in hair cortisol concentration between two batches within the same farm. These differences imply that hair should be sampled in more than one batch on farm, if the perspective of using hair cortisol is to assess long-term stress of pigs at the farm level. The previously reported mean or median values of a batch (ranging from 4.5 to 54.3 pg/mg; 19–23.26), for pigs sampled at the end of the finishing period (between 20 and 40 weeks of age), were comparable to the mean batch values obtained in this field study (ranging from 4.8 to 54.5 pg/mg). Differences in hair cortisol concentration between two batches within the same farm could potentially be attributed to changes in the nature or intensity of stressors. Hypotheses to further investigate could include changes in management practices, social stress due to pig behavior, thermal stress or health events (22, 35–39).

## Sample size for observing variability in hair cortisol concentration

The low interindividual differences in a few batches do not appear to be due to sampling bias resulting from a small number of sampled pigs. Indeed, the farm with the lowest percentage of pigs sampled per batch had high interindividual differences. Inversely,
the farm with the highest percentage of pigs sampled per batch showed very low interindividual differences. So, the variability in concentration within a batch seems to be a characteristic of this batch.

Characterizing the hair cortisol concentration of a batch requires an adequate sample size to calculate both the mean and the standard deviation. While a sample size of 15 pigs per batch can be sufficient to estimate the mean hair cortisol concentration in a batch, more than 24 pigs seem to be needed to better describe the variability between pigs. Indeed, the mean value of hair cortisol concentration in a drawn sample of 15 pigs differed relatively little from the initial sample of 24 pigs, resulting in a similar classification as low, intermediate, or high mean concentration. However, the differences in described variability between samples of 20 and 24 pigs suggest that there is still room for more precision in describing the variability with a larger sample size. That is why more than 24 pigs seem to be needed to better describe the variability between pigs.

## Conclusion

The hair cortisol concentration in finishing pigs on commercial farms showed large differences between pigs within a batch, between farms and to a lesser extent between batches of the same farm. The described variability suggests that there is an interest in further investigating if hair cortisol can be a biomarker to assess long-term stress in pigs on commercial farms. The differences between farms and batches confirm the need to assess whether it can be linked to long-term exposure of pigs to stressors. Insights about the sampling of finishing pigs were obtained in this study. More than 24 finishing pigs per batch are needed to describe the variability between pigs within a batch whereas 15 pigs could be sufficient to estimate the mean hair cortisol concentration of the batch. More than one batch is needed to describe the hair cortisol concentration at the farm level.

## Data availability statement

The datasets presented in this study can be found in an online repository. The names of the repository and accession number can be found at: https://zenodo.org/record/8359352 (40).

## Ethics statement

Ethical approval was not required for the study involving animals in accordance with the local legislation and institutional requirements because no ethical approval was required for this study, since the sampling of pig hair was carried out on commercial farms. According to French regulation (Rural and Sea fishing Code, book II, first title, chapter IV, section 6), ethic committees only provide authorizations for research conducted under experimental conditions. Hair sampling was non-invasive and was carried out with the approval of farmers.

## Author contributions

PL: Data curation, Formal analysis, Investigation, Validation, Visualization, Writing – original draft. ML-M: Conceptualization, Investigation, Supervision, Validation, Writing – review & editing. AL: Data curation, Formal analysis, Validation, Writing – review & editing. SG: Resources, Validation, Writing – review & editing. BL: Validation, Writing – review & editing. JH: Validation, Writing – review & editing. CF: Conceptualization, Funding acquisition, Supervision, Validation, Writing – review & editing. CB: Conceptualization, Funding acquisition, Investigation, Supervision, Validation, Writing – review & editing.

## Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was co-funded by the EU part of the HealthyLivestock project, funded by the EU Horizon 2020 Research and Innovation Program under grant agreement number 773436 and the Région Pays de la Loire (number of contract 2020_10347).

## Acknowledgments

The authors would like to thank all the farmers and herd veterinarians who contributed to this study for the time they dedicated and the talks we shared. Thanks to Tracy Delon who contributed to data collection. Thanks to Morgane Rémond and Alexandre Popiolek who contributed to farm visits. Thanks to Guillaume Morgant who contributed to measure hair cortisol concentration in the laboratory. PL used ChatGPT-3.5 (May 24 2023 Version) to translate into English the very first version of the manuscript then all the authors revised this draft.

## Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

## Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2023.1298756/full#supplementary-material
References

1. Farm Animal Welfare Council. Press statement. *The national archives* (1979) Available at: https://webarchive.nationalarchives.gov.uk/ukgwa/20121007104210/http://www.fawc.org.uk/pdf/fivefreedoms1979.pdf.

2. Fraser D, Weary DM, Pajor EA, Milligan BN. A scientific conception of animal welfare that reflects ethical concerns. *Anim Welf.* (1997) 6:187–205. doi: 10.1017/S0962728600019795

3. National Research Council. Stress and distress: definitions In: *Recognition and alleviation of distress in laboratory animals*. Washington, DC: National Academies Press (2008). 13–24. Available at: http://www.nap.edu/catalog/11931

4. Selye H. Stress and the general adaptation syndrome. *Br Med J.* (1950) 1:1383–92.

5. Mason JW. A re-evaluation of the concept of “non-specificity” in stress theory. *J Psychiatr Res.* (1971) 8:323–33.

6. Numberger J, Ritzmann M, Übel N, Eddicks M, Reese S, Zöls S. Ear tagging in piglets: the cortisol response with and without analgesia in comparison with castration and tail docking. *Animal*. (2016) 10:1864–70. doi: 10.1017/S1751731116000811

7. Koolhaas JM, Korte SM, De Boer SF, Van Der Vegt BJ, Van Reenen CG, Hopster H, et al. Coping styles in animals: current status in behavior and stress-physiology. *Neurosci Biobehav Rev.* (1999) 23:925–35. doi: 10.1016/S0149-7634(99)00026-3

8. Mormède P, Andanson S, Aupérin B, Beerda B, Guémené D, Malmkvist J, et al. Exploration of the hypothalamic-pituitary-adrenal function as a tool to evaluate animal welfare. *Physiol Behav.* (2007) 92:317–39. doi: 10.1016/j.physbeh.2006.12.003

9. Cirimele V, Kintz P, Dumestre V, Goullé JP, Ludes B. Identification of ten corticosteroids in human hair by liquid chromatography–ionspray mass spectrometry. *Forensic Sci Int.* (2000) 107:381–8. doi: 10.1016/S0379-0738(99)00180-2

10. Palme R, Rettenbacher S, Touma C, El-Bahr SM, Möstl E. Stress hormones in mammals and birds: comparative aspects regarding metabolism, excretion, and noninvasive measurement in Fecal samples. *Ann N Y Acad Sci.* (2005) 1040:162–71. doi: 10.1196/annals.1327.021

11. Meyer JS, Novak MA. Minireview: hair cortisol: a novel biomarker of hypothalamic-pituitary-adrenocortical activity. *Endocrinology*. (2012) 153:4120–7. doi: 10.1210/en.2012-1226

12. Kalliokoski O, Jellestad FK, Murison R. A systematic review of studies utilizing hair glucocorticoids as a measure of stress suggests the marker is more appropriate for quantifying short-term stressors. *Sci Rep.* (2019) 9:1997. doi: 10.1038/s41598-019-48517-2

13. González-De-La-Vara Del R, Valdez RA, Lemos-Ramirez V, Vázquez-Chagoyán JC, Villa-Godoy A, Romano MC. Effects of adrenocorticotropic hormone challenge and age on hair cortisol concentrations in dairy cattle. *Can J Vet Res.* (2011) 75:216–21.

14. Ashley NT, Barboza PS, Macbeth BJ, Janz DM, Cattet MRL, Booth RK, et al. Glucocorticosteroid concentrations in feces and hair of captive caribou and reindeer following adrenocorticotropic hormone challenge. *Gen Comp Endocrinol.* (2011) 172:382–91. doi: 10.1016/j.ygcen.2011.03.029

15. Macbeth BJ, Cattet MRL, Stenhouse GB, Gibeau ML, Janz DM. Hair cortisol concentration as a noninvasive measure of long-term stress in free-ranging grizzly bears (*Ursus arctos*): considerations with implications for other wildlife. *Can J Zool.* (2010) 88:935–49. doi: 10.1139/Z10-057

16. Yamanashi Y, Teramoto M, Morimura N, Hirata S, Suzuki J, Hayashi M, et al. Analysis of hair cortisol levels in captive chimpanzees: effect of various methods on cortisol stability and variability. *MethodsX*. (2016) 3:110–7. doi: 10.1016/j.mex.2016.01.004

17. Kapoor A, Schultz-Darken N, Ziegler TE. Radiolabel validation of cortisol in the hair of rhesus monkeys. *Psychoneuroendocrinology*. (2018) 97:190–5. doi: 10.1016/j.psyneuen.2018.07.022

18. Watson SA, Moore GP. Postnatal development of the hair cycle in the domestic pig. *J Anat.* (1990) 70:1–9.

19. Gavaud S, Hauagné K, Buchet A, Garcia Vinado I, Allard M, Lehébel A, et al. Effects of improved early-life conditions on health, welfare, and performance of pigs raised on a conventional farm. *Animal*. (2023) 17:100810. doi: 10.1016/j.animal.2023.100810

20. Bergamin C, Comin A, Corazzin M, Faustini M, Peric T, Scollo A, et al. Cortisol, DHEA, and sexual steroid concentrations in fattening pigs' hair. *Animals*. (2019) 9:345. doi: 10.3390/ani9060345

21. Montillo M, Rota Nodari S, Peric T, Polloni A, Corazzin M, Bergamin C, et al. Steroids in pig hair and welfare evaluation systems: combined approaches to improve management in pig breeding? *Vet Ital.* (2020) 56:177–84. doi: 10.12834/VetIt.1974.11885.1

22. Carroll GA, Boyle LA, Hanlon A, Palmer MA, Collins L, Griffin K, et al. Identifying physiological measures of lifetime welfare status in pigs: exploring the usefulness of haptoglobin, C-reactive protein and hair cortisol sampled at the time of slaughter. *Ir Vet J.* (2018) 71:8. doi: 10.1186/s13620-018-0118-0

23. Casal N, Manteca X, Peña LR, Bassols A, Fàbrega E. Analysis of cortisol in hair samples as an indicator of stress in pigs. *J Vet Behav.* (2017) 19:1–6. doi: 10.1016/j.jveb.2017.01.002

24. République Française. Cahier des charges du label rouge n° LA 02/05 homologué par l'arrêté du 11 août 2017. *JORF du 23 août 2017-Bulletin officiel du Ministère de l'agriculture et de l'alimentation n°2017-36.* (2017).

25. Levallois P, Leblanc-Maridor M, Scollo A, Ferrari P, Belloc C, Fourichon C. Combining several indicators to assess the effectiveness of tailor-made health plans in pig farms. *Peer Community J.* (2023). 3:e86. doi: 10.24072/pcjournal.318

26. Heimbürge S, Kanitz E, Tuchscherer A, Otten W. Within a hair's breadth – factors influencing hair cortisol levels in pigs and cattle. *Gen Comp Endocrinol.* (2020) 288:113359. doi: 10.1016/j.ygcen.2019.113359

27. Ito N, Ito T, Kromminga A, Bettermann A, Takigawa M, Kees F, et al. Human hair follicles display a functional equivalent of the hypothalamic-pituitary-adrenal (HPA) axis and synthesize cortisol. *FASEB J.* (2005) 19:1332–4. doi: 10.1096/fj.04-1968fje

28. Otten W, Heimbürge S, Kanitz E, Tuchscherer A. It's getting hairy – external contamination may affect the validity of hair cortisol as an indicator of stress in pigs and cattle. *Gen Comp Endocrinol.* (2020) 295:113531. doi: 10.1016/j.ygcen.2020.113531

29. Otten W, Heimbürge S, Tuchscherer A, Kanitz E. The age of hair matters – the incorporation of cortisol by external contamination is enhanced in distal hair segments of pigs and cattle. *Animal*. (2022) 16:100495. doi: 10.1016/j.animal.2022.100495

30. Cohen J. Statistical power analysis for the behavioral sciences. 2nd ed. Hillsdale, NJ: Lawrence Erlbaum Associates (1988).

31. Wiechers DH, Brunner S, Herbrandt S, Kemper N, Fels M. Analysis of hair cortisol as an indicator of chronic stress in pigs in two different farrowing systems. *Front Vet Sci.* (2021) 8:605078. doi: 10.3389/fvets.2021.605078

32. Otten W, Puppe B, Kanitz E, Schön PC, Stabenow B. Effects of dominance and familiarity on behaviour and plasma stress hormones in growing pigs during social confrontation. *J Vet Med Ser A.* (1999) 46:277–92. doi: 10.1046/j.1439-0442.1999.00216.x

33. Harper GC, Makatouni A. Consumer perception of organic food production and farm animal welfare. *Br Food J.* (2002) 104:287–99. doi: 10.1108/00070700210425723

34. Duval E, von Keyserlingk MAG, Lecorps B. Organic dairy cattle: do European Union regulations promote animal welfare? *Animals*. (2020) 10:1786. doi: 10.3390/ani10101786

35. Bacci ML, Nannoni E, Govoni N, Scorrano F, Zannoni A, Forni M, et al. Hair cortisol determination in sows in two consecutive reproductive cycles. *Reprod Biol.* (2014) 14:218–23. doi: 10.1016/j.repbio.2014.06.001

36. Pathak PK, Roychoudhury R, Saharia J, Borah MC, Dutta DJ, Bhuyan R, et al. Impact of seasonal thermal stress on physiological and blood biochemical parameters in pigs under different dietary energy levels. *Trop Anim Health Prod.* (2018) 50:1025–32. doi: 10.1007/s11250-018-1526-6

37. Comin A, Peric T, Corazzin M, Veronesi MC, Meloni T, Zufferli V, et al. Hair cortisol as a marker of hypothalamic-pituitary-adrenal axis activation in Friesian dairy cows clinically or physiologically compromised. *Livest Sci.* (2013) 152:36–41. doi: 10.1016/j.livsci.2012.11.021

38. Burnett TA, Madureira AML, Silper BF, Tahmasbi A, Nadalin A, Veira DM, et al. Relationship of concentrations of cortisol in hair with health, biomarkers in blood, and reproductive status in dairy cows. *J Dairy Sci.* (2015) 98:4414–26. doi: 10.3168/jds.2014-8871

39. Morgan L, Itin-Shwartz B, Koren L, Meyer JS, Matas D, Younis A, et al. Physiological and economic benefits of abandoning invasive surgical procedures and enhancing animal welfare in swine production. *Sci Rep.* (2019) 9:16093. doi: 10.1038/s41598-019-52677-6

40. Levallois P, Leblanc-Maridor M, Lehebel A, Gavaud S, Lieubeau B, Herve J, et al. Dataset of hair cortisol concentration in 950 finishing pigs on 20 commercial farms. *Zenodo* (2023). Available at: https://zenodo.org/record/8359352