Intertwining of the IGF system and animal welfare

TL;DR

This study explores how different housing conditions affect the IGF system in sows, suggesting it could help objectively assess animal welfare.

Contribution

The study introduces a novel bioassay to evaluate IGF system activity as a potential indicator of animal welfare.

Findings

Differential housing conditions altered serum IGF bioactivities in sows.

IGF system factors varied depending on housing conditions during pregnancy and lactation.

Principal component analysis suggests IGF factors may distinguish sows by housing environment.

Abstract

To date, husbandry categories are defined based on the human concept of animal needs, however, to objectively assess animal welfare in defined husbandry systems, animal-related indicators must be evaluated. Here, we studied the impact of differential housing conditions on factors of the insulin-like growth factor system (IGF-system) in mother sows either housed in a conventional husbandry system corresponding to the regulations of the German husbandry form 1, or in the associated ecological facility, corresponding to the German husbandry form 4/5. Utilizing a proprietary bioassay that determines the activity of the IGF system at the level of intracellular signal transduction, we revealed altered serum IGF bioactivities in differentially housed sows. Moreover, we present a comprehensive analysis of the individual factors of the IGF system in serum of pregnant and lactating sows to elucidate which of these are modified depending on the housing conditions. An unsupervised principal component analysis suggests that IGF system factors might distinguish sows based on their housing environment. These results imply that the IGF system provides a solid basis for the documentation and objective assessment of animal health and welfare. However, a robust biomarker system for animal welfare assessment will likely require a multifactorial and integrative approach.

The online version contains supplementary material available at 10.1038/s41598-026-42315-3.

Figures

Graphical summary of factors affecting IGF bioactivity.

Experimental design to study effects of the housing environment on factors of the somatotropic axis. Blood and saliva samples were obtained from differentially housed sows 30 and 105 days post coitum (dpc) as well as 8 days post partum (dpp) to analyze abundance and biological activity of candidate factors (CORT: cortisol).

Effects of the husbandry system on cortisol and IGF bioactivity. In differentially housed sows, cortisol concentrations were measured (A) in serum (*n* = 9–10) and (B) in saliva, for which values were normalized to protein content to account for the varying water content (*n* = 9–10). (C) Normalized IGF-dependent bioactivity was determined based on recombinant human IGF1 standards as reference (*n* = 20). Statistical comparisons are based on estimated marginal means. Brackets indicate statistically significant differences with the corresponding p-values shown above.

Husbandry dependent IGF concentrations. Concentrations of total IGF1 (A) and IGF2 (B) were determined in serum of differentially housed pregnant and lactating sows (*n* = 10). Statistical comparisons are based on estimated marginal means. Brackets indicate statistically significant differences with the corresponding p-values shown above.

Husbandry effects on regulators of IGF bioactivity (A) IGFBP2 concentrations were increased in conventional housing resulting in an altered IGFBP2/IGFBP3 ratio (*n* = 18). (B) STC1 concentrations in serum of pregnant and lactating sows (*n* = 10). Statistical comparisons are based on estimated marginal means. Brackets indicate statistically significant differences with the corresponding p-values shown above.

PCA results for candidate parameters determined in the study (A) with a reduced dataset containing only complete cases and (B) with a greater dataset containing imputed data for missing values in the original dataset.

Overview on the differential impact of the husbandry system on the abundance and bioactivity of factors of the IGF system.

Full text

Introduction

The increasing interest of consumers in the farming conditions of animals, motivates considerations about animal welfare labeling as discussed in the framework of the EU Farm to Fork Strategy. In Germany, voluntary labeling, based on different husbandry forms concerning the type of shelter, available space, and opportunities to express natural behaviors, has already been introduced (husbandry forms 1–5)^1^. However, these categories are defined based on the human concept of animal needs; to assess whether animals actually benefit from the superior husbandry forms, animal-related indicators must be recorded and evaluated. The Welfare Quality^®^ Assessment Protocol^2^ is one attempt to assess animal welfare but it only records externally recognizable characteristics such as behavioral disorders, diseases or injuries and does not include any laboratory tests in stages before adverse manifestations. However, in order to be able to take timely action, it must be achieved to detect discomfort and stress before the onset of mental and physical damage.

To date, hormones of the pituitary-adrenal axis, such as cortisol, are commonly used as biomarkers, particularly for acute stress. However, these hormones are less suitable as indicators of long-term stress and well-being since they are characterized by relatively short half-lives and a pulsatile secretion usually following circadian and seasonal rhythms, which generally necessitates repeated measures^3^. More importantly, several lines of evidence suggest that chronic stress in pigs might not even be effective in deflecting basal cortisol levels^4^. Even if the triggering stimulus is maintained, plasma cortisol levels usually decline after the acute response. For instance, pigs receiving unpredictable and inescapable electric shocks for a month have plasma cortisol levels similar to those of control animals not receiving any shock. Even though, they still react behaviorally to the situation (agitation followed by inhibition), suggesting that they are continuously stressed^5^. While hair cortisol is suggested to better represent chronic stress, there are still several practical limitations such as local and external sources of cortisol “contaminating” the hair and differences in hair growth rate within the same body region that could affect hair cortisol concentrations^6^. Taken together, hormones of the pituitary-adrenal axis, such as cortisol, are not well suited to assess animal welfare in dependence on the husbandry conditions, particularly in pigs.

The pituitary-liver axis (somatotropic axis), which initially induces growth hormone (GH) and subsequently controls the insulin-like growth factor (IGF) system, is also regulated by stress^7^. Here, we propose that the activity of this system is informative not only for short- and medium-term stress, as we have already shown^8,9^, but also for longer-term stress and well-being on the other hand.

The IGF system includes two forms of insulin-like growth factors (IGF1 and IGF2), their receptors, and IGF-binding proteins (IGFBPs) that regulate the bioavailability of IGFs^10^. Additionally, its activity is modulated by several other factors. For example, the binding of IGFBPs is regulated by proteases such as PAPP-A and PAPP-A2, which are in turn modulated by inhibitors such as stanniocalcin 1 and 2 (STC1 and STC2)^11,12^. Recently, it was shown that also calcium affects IGF signaling^13^. Figure 1 provides an illustration of these several layers of IGF bioactivity modulation. The IGF1 receptor (IGF1R) belongs to the receptor tyrosine kinase family and induces pathways such as PI3K/AKT/mTOR and MAPK/ERK. Our proprietary BIRA assay is based on the detection of AKT phosphorylation to assess IGF-related bioactivity downstream of the IGF1R, thereby accounting for the integration of multiple signals exerted by free IGFs, IGFBPs, IGFBP-proteases, and protease inhibitors^14^.

Here, we consider both, the overall IGF-related bioactivity as well as quantities of individual factors to investigate the impact of the environment, and welfare aspects in particular, on the IGF system. More generally, we aim to clarify if the husbandry categories for pigs, defined by humans, represent distinct levels of welfare that can be measured objectively. To approach this question, we tested the feasibility of discriminating between husbandry systems based on selected candidate parameters utilizing unsupervised principal component analysis (PCA).

Fig. 1. Graphical summary of factors affecting IGF bioactivity.

Methods

Housing conditions in conventional and ecological facility

For the study, German landrace sows were housed in the experimental animal facility at FBN Dummerstorf in compliance with European animal welfare regulations for appropriate husbandry and experimental conditions as well as the national Animal Welfare Act and Animal Protection Farm-Animal-Husbandry-Ordinance (German TierSchNutztV). The animals were divided in two experimental groups, one being housed in the conventional husbandry system corresponding to the definition of the German husbandry form 1, and the other being housed in the associated ecological (eco) facility, corresponding to the definition of the German husbandry form 4/5 (Fig. 2). The housing experiment involved the impregnation of the sows and consecutive rearing of up to three litters. For each parity, the sows age at insemination is given in Supplementary Table S1 and S3. The sows were mated with 10 different boars, covering individuals from both husbandry systems. Mating assignments were determined by the institute’s pig herd manager, without external influence. On the 105th day of gestation, the sows were moved from group pens into the farrowing pens.

For the conventional husbandry, pregnant sows were housed in three-area pens (walking area - concrete slatted floor, lying area - flat paved with litter, feeding stalls - self-catching stalls with partially slatted floor) for 4 sows each with an area of 22 m² for old sows and 15.9 m² for gilts (corresponds to a pen area of 5.5 m² per old sow and 4.0 m² per gilts). Animals were fed twice a day with a complete feed for pregnant sows (Sauen Struktur tragend, Trede and von Pein, Itzehoe, Germany), according to the recommendations of the Society for Nutritional Physiology for the energy and nutrient supply of pigs^15^. The animals were provided with enrichment material in the form of long straw in racks, shavings and spelt litter as well as chains with wood. Water was accessible at all times via a drinking trough. The farrowing pen offered the sow and suckling piglets an area of 6.5 m², while allowing for the restrainment of the sows from the 114th day of pregnancy until the piglets were four days old. The pens are equipped with plastic slatted floors, cast-iron elements in the sow’s lying area and a piglet nest (0.75 m²) with a climate cover and a heat lamp. The sows were fed twice daily as required (gradually increasing the daily feed quantity from 1 kg after birth to 7–8 kg 14 days post partum) with a complete feed for lactating sows (Sauen Struktur säugend, Trede and von Pein, Itzehoe, Germany). In addition to the bedding mix, wood on chains was provided as enrichment material. The birth was induced by administering 2 ml PGF Veyx (Cloprostenol, Prostaglandin F2α) one day before the calculated day of birth and 1 ml Longacton (Carbetocin, Oxytocin) on the day of birth.

In the eco facility, the group pen offers a stall area of 66 m² and an exercise area of 45 m² housing 8–16 sows during the experiment. This corresponds to a space of minimally 4.1 m² stall area plus 2.8 m² outside exercise area per sow. The animals are fed a high-fiber, GMO-free, organic feed (B Sauenfutter NT Pell., Cera-Green) via a call-off station. Fixed amounts of feed (depending on their condition and stage of pregnancy) can be called off throughout the day and water is available ad libitum via two trough drinkers. The pens are littered with straw in all areas and an automatic cattle brush is provided as an additional activity material. The farrowing pens offer approx. 8.4 m² of littered inside pen space and 5 m² of outside area. A piglet nest (0.75 m²) with underfloor heating and a heat lamp is integrated into the pen. The sows and piglets can move freely in the pen, but the sows can be restrained for short periods. The sows are fed a high-fiber, nutritionally balanced organic lactation feed (B Sauenfutter LAC Pell., CeraGreeen) several times a day following the recommendations of the Society for Nutritional Physiology for the energy and nutrient supply of pigs^15^. Composition and nutritional contents of all animal feed is declared in Supplemental Table S2. Water is available to the animals at all times via a drinking trough and a mother-child drinking trough ad libitum. There is no routine induction of labor with prostaglandin and/or oxytocin.

Monitoring of general health parameters of mother sows and offspring

Generally, all animals are cared for and looked after by experienced staff of the experimental animal facility at FBN recording any abnormalities that indicate illness. For the litters in particular, a number of parameters are routinely recorded, such as the total number and gender of the piglets and their birth weight. Moreover, any abnormalities such as the occurrence of still born piglets, piglets with defective leg muscles (splay leg cases), and low weight and feeble animals are documented (see Supplementary Table S1). Of note, while all piglets per litter were considered for the ratio of males and the ratio of still born piglets and splay leg cases, only piglets born alive were considered for the calculation of the mean birth weight and ratio of low-birth-weight piglets in this study.

Sampling of blood and saliva

At day 30 and 105 of gestation as well as day 8 after farrowing, approximately 10 ml of blood were collected from each sow in 50 ml tubes by puncture of the cranial vena cava under fixation. For serum preparation, whole blood was centrifuged at 1.500 g for 10 min. To improve the purity of the samples, the supernatant was transferred to new tubes and centrifuged again at 1.500 g for 5 min before the serum was transferred to cryo-tubes and stored at −80 °C until usage.

Additionally, non-invasive sampling of saliva using swabs (5001.06, Salimetrics, Carlsbad, USA) was conducted before blood sampling. For this, the swabs were kept in the sows’ mouths for at least one minute until satisfactory moistening was achieved. To collect the saliva, the swabs were placed in syringe without cannula and put in 15 ml tubes for a centrifugation at 2.000 g for 10 min. Saliva was then transferred in cryo-tubes and stored at −80 °C until usage. Samplings were generally conducted in the morning (7:00–7:30). Serum and saliva samples were used for a number of assays as outlined in the following. An overview of the experimental design is illustrated in Fig. 2.

Fig. 2. Experimental design to study effects of the housing environment on factors of the somatotropic axis. Blood and saliva samples were obtained from differentially housed sows 30 and 105 days post coitum (dpc) as well as 8 days post partum (dpp) to analyze abundance and biological activity of candidate factors (CORT: cortisol).

Quantification of cortisol concentrations

For the determination of cortisol concentrations in serum and saliva, respectively, commercial ELISA Kits (DEH3388 and DES6611, demeditec, Kiel, Germany) were used according to the instructions of the manufacturer. Briefly, serum samples were homogenized by vortexing, while saliva samples were diluted 1:2 in sample buffer and centrifuged at 2.000 g for 3 min to pellet any contaminants. 10 µl (for serum) and 50 µl (for saliva) of each sample, control and calibrator were dispensed in duplicates in the anti-cortisol antibody coated 96-well plate. 200 µl and 50 µl enzyme conjugate, respectively, were added to the wells and incubated for 1 h at room temperature under gentle shaking followed by several washing steps. Subsequently, 200 µl of substrate buffer was pipetted and incubated for 30 min at room temperature in the dark for both analyzed matrices. Finally, 50 µl of stop solution were added to each well and absorbance was measured at 450 nm. Cortisol concentrations were calculated based on a four-parameter logistics curve fit of the calibrator values. For saliva samples, protein concentrations were determined using BCA (Bicinchoninic acid) tests to normalize the obtained results regarding potential dilution effects due to variations in salivation and drinking behavior of the sows.

Determination of IGF-related bioactivity - BIRA assay

For the BIRA (IGFBP2-enhanced IGF-related AKT activation) assay we used a transfected human embryonic kidney fibroblast cell line (HEK293), obtained from the American Type Culture Collection (ATCC, Rockville, MD, USA). HEK293 cells were transfected with 8 µg pCMV-int-mIGFBP2 (XhoI linearized) and 0.8 µg EcoRI linearized neomycin resistance plasmid pSV2neo (Clontech, Heidelberg, Germany)^16^, which renders the cells more sensitive to bioactive IGF^14^. HEK293 clones were cultured in Eagle’s minimal essential medium (EMEM) supplemented with 10% fetal bovine serum (FBS; 10500-064; Gibco, Carlsbad, CA, USA) and 1% 100× penicillin/streptomycin/amphotericin B mixture (882087; Lonza, Gampel, Switzerland) at 37 °C and 5% CO_2_.

To prepare the bioassay, cells were detached using 0.25% trypsin in phosphate-buffered saline (PBS, 137 mM NaCl, 2.7 mM KCl, 10 mM Na_2_HPO_4_, 1.8 mM KH_2_PO_4_) and seeded in 24-well-plates with 1–2 × 10^5^ cells per well in 500 µL cell culture medium. After reaching 90% confluence, medium was exchanged to EMEM with 0.5% FBS for 16–20 h to reduce triggering stimuli originating from the serum. After the cells were washed with PBS, 240 µl of serum samples (diluted 1:10 in PBS) and serial dilutions from 6.25ng/ml to 100 ng/ml of human recombinant IGF1 (Cat. #100 − 11, PeproTech, Inc., Rocky Hill, NJ, USA), were added for 20 min at 37 °C. Subsequently, the test media were discarded and cells were washed twice with PBS to eliminate any protein residues, before the cells were lyzed in 100 µL lysis buffer (150mM NaCl, 50mM Tris/Cl, 1% Tween 20, Complete Mini protease inhibitor cocktail dissolved in VE water).

For sodium dodecylsulfate (SDS)-polyacrylamide gel electrophoresis (PAGE), 100 µl 2x Laemmli (125mM tris(hydroxymethyl)-aminomethane, 4% SDS, 20% glycerine) containing 0.01% bromophenol blue and 4% beta-mercaptoethanol were added before lysates were homogenized for 30 min at 45 °C in an ultrasonic bath. Samples were denatured at 95 °C for 5 min prior to PAGE utilizing the Bio-Rad TGX Stain-Free FastCast Acrylamide-Kit (Bio-Rad Laboratories GmbH, Munich, Germany). All gels were blotted onto polyvinylidene fluoride membranes (PVDF, pore size 0.45 μm, Carl Roth GmbH + Co. KG, Karlsruhe, Germany), blocked for 1 h with 3% milk powder in Tris-buffered saline with Tween20 (TBST) and incubated overnight with a primary antibody for AKT phosphorylated at serine 473 (CST #4058, dilution 1:1000, Cell Signaling Technologies). After washing, the membranes were incubated for 2 h with a secondary antibody (antirabbit IgG HRP, CST #7074, dilution 1:2000, Cell Signaling Technologies). Bands were visualized using Lumigen ECL Ultra (Lumigen Inc., Southfield, MI, USA) in a Bio-Rad Chemi-Doc MP system (Bio-Rad Laboratories GmbH). Images were analyzed using Image Lab Software Ver. 6.0.1 (Bio-Rad Laboratories GmbH, Hercules, CA, USA) and signal intensities normalized to total protein concentration per lane. Finally, values of the IGF1 standards were used to determine IGF1 equivalents in ng/ml for all samples. However, since for complex matrices like serum, phosphorylation of AKT cannot be translated into IGF1 concentrations in ng/ml directly, we used the mean values for standard housing at 30dpc as reference value to calculate a relative IGF-dependent bioactivity for the investigated samples from differentially housed pregnant (105 dpc) and lactating (8dpp) sows.

Quantification of IGF1 and IGF2 concentrations

For the determination of total IGF1 and IGF2 concentrations in serum commercial ELISA Kits (E20 and E30, Mediagnost, Reutlingen, Germany) were used according to the manufacturers protocols. In short, serum samples were briefly vortexed for sample homogenization and diluted in sample buffer (1:21 for IGF1, 1:404 for IGF2). Two serum samples with known IGF concentrations each were included for assay quality control and prepared similarly. The microtiter plate was prepared with 80 µl/50 µl antibody conjugate before 20 µl/50 µl of samples, controls, blanks and calibrators were added to each well and incubated for 1 h/2 h at room temperature under gentle shaking. After multiple washing steps, 100 µl of enzyme conjugate was added for 30 min at room temperature under gentle shaking. Following further washing steps, 100 µl of substrate buffer was pipetted and incubated for 15 min at room temperature in the dark, before 100 µl of stop solution was added to each well. Absorbance was measured at 450 nm and IGF concentrations were calculated based on a third-degree polynomial curve fit of the calibrator values.

Determination of IGFBP concentrations

Serum samples and serial dilutions of recombinant human IGFBP standards (R&D Systems, Wiesbaden, Germany) were diluted 1:20 and boiled in sample buffer [312.5 mM Tris (pH 6.8), 50% (w/v) glycerol, 5mM EDTA, 1% (w/v) SDS, and 0.02% bromophenol blue] for 5 min. Proteins were separated by 12% SDS-PAGE using TGX Stain-Free™ FastCast™ acrylamide kit (Bio-Rad Laboratories GmbH, Feldkirchen, Germany) and transferred to a PVDF membrane (Millipore, Bedford, USA). The blots were blocked with 1% fish gelatine in TBST and then incubated with biotin labeled human IGF2 (1:500; BioIGF2-10; ibt-systems, Binzwangen, Germany), followed by streptavidin conjugate incubation and a detection step with Lumigen ECL Ultra (TMA-6) (Lumigen, Southfield, USA). Bands were visualized on ChemiDoc XRS + gel imaging system (Bio-Rad Laboratories GmbH, Feldkirchen, Germany) and quantified using Image Lab software Ver. 6.0.1 (Bio-Rad Laboratories GmbH, Hercules, CA, USA). Signal intensities were normalized using stain-free gel images and IGFBP concentrations were quantified using the standard serial dilutions. For each sample two technical replicates were generated on independent gels. If the results of technical replicates differed considerably, two additional technical replicates were generated and average values were used for downstream analyses.

Detection of protease activity

To test for proteolytic cleavage of IGFBPs in serum, samples were diluted 1:5 in PBS and incubated in a thermal shaker at 39 °C and 40 rpm for up to 44 h. At different time points (4, 20, 28 & 44 h), samples were mixed and centrifuged, an aliqot was taken and thoroughly mixed with an equal amount of 2x sample buffer (312.5 mM Tris (pH 6.8), 50% (w/v) glycerol, 5mM EDTA, 1% (w/v) SDS, and 0.02% bromophenol blue). Samples without incubation served as positive controls. Detection of intact IGFBPs followed the protocol for SDS-PAGE and Western ligand Blotting as described above for the “Determination of IGFBP concentrations”.

Quantification of STC1 concentrations

For the quantitative detection of STC1 levels in serum a commercial ELISA kit (EK1404, Boster Biological Technology, Pleasanton, USA) was used according to the manufacturer´s protocol with some adjustments accounting for the highly variable STC1 concentrations recognized in a pilot experiment. For a wide concentration range of standards, a 10ng/ml stock solution was diluted to 5000, 2000, 1000, 500, 250, 125, 62.5 and 0 pg/ml. 100 µl of standards and serum samples were pipetted in duplicates in a 96-well plate and incubated for 1.5 h at 37 °C. After discarding the fluids, 100 µl of a biotinylated polyclonal STC1 antibody was added per well and incubated for 1 h at 37 °C. After washing, 100 µl per well of an avidin-biotin-peroxidase complex were pipetted and incubated for 30 min at 37 °C. Following further thorough washing steps, 90 µl color developing reagent per well was pipetted and kept for a reduced incubation time of 20 min at room temperature in the dark, before 100 µl stop solution were added. Absorbance was measured at 450 nm and STC1 concentrations were calculated based on a four-parameter logistics curve fit of the standard concentrations.

Determination of calcium concentrations

To determine calcium concentrations in serum, a commercial kit was utilized (Calcium oCPC, LT-CA0100, Labor+Technik Eberhard Lehmann GmbH, Germany) with some adaptations to allow for a higher throughput using a 96-well format instead of single cuvettes. For this, serum samples were first vortexed for homogenization before 2 µl per sample were transferred into a 96-well plate in duplicates. Similarly, 2 µl of the provided calcium standard was pipetted in duplicates. The cresolphthalein-complexon and an alkaline buffer, both included in the kit, were mixed 1:1 and 200 µl of the mix were added to each well. Wells containing neither a standard nor a sample were used as blank. After an incubation for 10 min at room temperature extinction was measured at 546 nm in a plate reader. Calcium concentrations were calculated based on the provided standard according to the manufacturers’ protocol.

Statistical analysis

In this study, we detected eight serum parameters, with cortisol additionally measured in saliva, in 15 mother sows housed in two different systems: a conventional (n = 7) and an ecological facility (n = 8). Additionally, phenotypic data of their litters were collected (Supplementary Table S1). Each pregnancy was considered an experimental unit, and data were collected across up to three pregnancies per sow, resulting in a total of 20 litters per husbandry system (Supplementary Table S3). For each parameter, we performed an outlier test (ROUT, Q = 1%), to curate the data before further statistical tests. Data processing and statistical analyses were conducted in R^17^ (version 4.4.2). A power analysis using the pwr package^18^ determined that a sample size of 20 per group could detect effect sizes of Cohen’s d = 0.9 with a power of 0.80, classified as a large effect. To analyze the effects of housing and time on the detected neuroendocrine parameters, we fitted linear mixed-effects models (LMMs) using the lmerTest package^19^, which allows for hypothesis testing of fixed effects based on Satterthwaite’s degrees of freedom approximation. The response variables were log-transformed to meet model assumptions. The model included husbandry conditions (housing: conventional or eco), sampling time points (time: 105dpc or 8 dpp), and their interaction (housing × time) as fixed effects.

To account for the hierarchical structure of the data, we included a random intercept for sows and an additional random intercept for litters nested within sows (1 | sow/litter). This structure reflects that each sow gave birth to up to three litters. The model was fitted using restricted maximum likelihood (REML), and parameter estimation was optimized using the bobyqa optimizer. To assess the significance of fixed effects, we conducted type II Wald chi-square tests using the car package^20^. We also performed pairwise comparisons of housing conditions within each time point (and vice versa) using estimated marginal means (EMMs) via the emmeans package^21^. Pairwise contrasts were computed using Tukey’s adjustment for multiple comparisons to control for Type I error. The degrees of freedom were approximated using the Satterthwaite method and results were presented on the response scale. Standardized effect sizes (Cohen’s d) were calculated from the EMMs, accounting for uncertainty in both the estimated effects and the population standard deviation.

Principal component analysis and data imputation

Since the long-term ambition is to develop an objective biomarker for animal welfare, the application of unbiased methods such as unsupervised machine learning algorithms is an adequate method to test for the eligibility of candidate parameters. To investigate whether selected parameters could differentiate between the husbandry systems, we applied an unsupervised principal component analysis (PCA). For this, we condensed the data table to include only chosen parameters found to be deflected in dependence on the husbandry system or considered as the gold standard for stress detection (cortisol concentration). Since PCA requires complete cases, the analysis was initially limited to eight conventional and nine ecological litters, reducing sample size by over half and consequently lowering statistical power. Noteworthy, this reduction could also introduce bias.

To mitigate this issue, we applied an imputation method to replace missing values, ensuring that all cases were preserved. Notably, Missing values were imputed separately for each husbandry group to preserve group-specific characteristics. The K-Nearest Neighbors (KNN) imputer, accessible through the scikit-learn package in Python^22^, was used to identify k-nearest neighbors for incomplete cases based on Euclidean distance that is suitable for continuous variables. Missing values were replaced with the mean of the nearest neighbors, applying uniform weighting. Neighbors adaptively adjusted when multiple missing values were present, ensuring context-sensitive imputation. To test for the feasibility of discriminating between the two husbandry systems based on the collected parameters, we conducted PCA using the mixOmics package^23^ for both the original complete cases and the imputed dataset.

Results

The husbandry system does not affect quantity and weight of piglets but potentially the incidence of abnormalities

Given the involvement of the IGF system in growth and development, we examined the characteristics of litters born within our study cohort in the conventional and eco housing facilities. We could analyze 20 litters per husbandry system (Supplementary Table S1). During the study, sows gave birth to maximally three litters following inseminations at on average 310, 470 and 630 days of age for the standard husbandry and 270, 430 and 600 days for the ecological husbandry in our tested cohort. In both groups, data were collected relatively evenly across all parities (n = 7–8 for 1 st, n = 6–7 for 2nd and n = 6 for 3rd parity). The data revealed no significant differences in total litter size or sex distribution between conventional and ecological husbandry (Table 1). Also the proportion of low-birth-weight piglets (less than 1.1 kg at birth) was comparable between the two systems, with the most underweight piglets occurring in the first litters (conv: 56.3% (IQR 41.2–61.6), eco: 22.7% (IQR 11.8–37.9)), which is in line with data of other studies^24^.

However, the incidence of litters with stillborn piglets was higher under conventional conditions. Furthermore, the occurrence of splay leg cases (piglets with defective leg muscles) was twice as high in the conventional system compared to the eco facility. To account for the potential confounding effects of the boars, we tested their influence on litter parameters. Except for the proportion of splay leg cases, boars had no significant impact. Only 9 out of 40 litters exhibited splay leg cases at all, with one boar showing relatively high occurrences (10.5% and 15.0% of piglets in a litter), leading to a significant confounding effect (Supplementary Table S4). However, this trait was not included in further analyses.

Table 1. Comparison of basic features of litters depending on the husbandry system.Litter featureConventional (n = 20)Eco (n = 20)Total number of piglets per litter17.9 ± 3.116.7 ± 3.9Males per litter50.1% ± 13.7%47.7% ± 10.6%Average birth weight1.31 ± 0.24 kg1.37 ± 0.23 kgLow birth weight piglets per litter29.8 ± 24.9%24.5 ± 21.9%Litters with still born70.0%50.0%Litters with splay leg cases30.0%15.0%

Cortisol levels in sows are similar in conventional and eco husbandry

Total cortisol concentrations were determined in the serum of pregnant and lactating sows across the two husbandry systems. As it represents a much less invasive approach, saliva cortisol concentrations are now more commonly analyzed to assess stress levels in pigs. Thus, we additionally detected cortisol in saliva. Since the housing conditions did not allow for controlling the hydration status of the sows, saliva cortisol concentrations were normalized to saliva protein content to account for potential dilution effects due to variations in salivary flow and drinking behavior. Of note, saliva sampling was generally performed before blood sampling.

For serum, cortisol concentrations did not differ significantly between husbandry systems or sampling time points (Table 2; Fig. 3A). Values ranged between 17.1 and 54.7 ng/ml. For saliva, slightly different results were obtained. Pregnant sows tended to have higher salivary cortisol concentrations in the conventional system compared to the eco husbandry (10.6 ng/mg (IQR 6.7–17.0) vs. 5.3 ng/mg (IQR 4.3–9.5)), while in lactating sows levels were similar in both husbandry systems (3.8 ng/mg (IQR 2.8–5.1) vs. 5.3 ng/mg (IQR 3.2–11.2)) as the values obtained of conventional husbandry have converged with those of the eco husbandry (Fig. 3B). Despite these numerical trends, none of the differences reached statistical significance (Table 2).

IGF bioactivity is affected by the housing environment

The proprietary BIRA assay was used to specifically detect IGF1R-mediated bioactivity, quantified based on AKT phosphorylation levels^14^. Since serum is a complex biological matrix potentially containing various compounds that could induce AKT phosphorylation, we verified the assay’s specificity using the IGF1R inhibitor BMS-754,807. As expected, PBS with 50ng/ml IGF1 yielded a strong signal, while adding 0.1µM BMS reduced signal intensity to levels comparable to the negative control (PBS) (Supplementary Fig. S1, see also^14^. Similarly, in serum diluted 1:5 in PBS, the signal intensity dropped to negative control levels upon BMS supplementation, confirming that serum-induced AKT phosphorylation was specifically mediated by IGF1R. Even though IGF1 standards were utilized to assess IGF-dependent bioactivities, we provide data on IGF1R-mediated bioactivity as relative quantification to avoid potentially misleading concentration values (Fig. 3C).

The comparison of serum IGF bioactivities of differentially housed sows revealed significantly increased values under standard conditions (Table 2; Fig. 3C) in pregnant sows (86.9 (IQR 53.1–108.5) vs. 51.24 (IQR 35.4–80.7), p = 0.018) and a trend to increased values under standard conditions in lactating sows (91.3 (IQR 67.6–118.2) vs. 70.5 (IQR 45.5–86.1), p = 0.109). Calculated effect sizes for the husbandry conditions were 1.11 ± 0.46 at 105dpc, rendering this observation a large effect (Supplementary Table S5). In contrast, the time of sampling did not significantly affect IGF bioactivity, a finding further supported by the additional measurements at 30dpc (Supplementary Table S3).

Fig. 3. Effects of the husbandry system on cortisol and IGF bioactivity. In differentially housed sows, cortisol concentrations were measured (A) in serum (n = 9–10) and (B) in saliva, for which values were normalized to protein content to account for the varying water content (n = 9–10). (C) Normalized IGF-dependent bioactivity was determined based on recombinant human IGF1 standards as reference (n = 20). Statistical comparisons are based on estimated marginal means. Brackets indicate statistically significant differences with the corresponding p-values shown above.

Both IGF1 and IGF2 concentrations were measured in a subset of serum samples. At 105dpc, IGF1 levels showed a trend toward higher concentrations in conventional husbandry compared to the eco system (154.7 ng/ml (IQR 59.3–288.4) vs. 72.3 ng/ml (IQR 59.2–123.0), p = 0.064); however, this difference was no longer evident at 8dpp (222.7 ng/ml (IQR 189.0 −325.2) vs. 222.8 ng/ml (IQR 149.4–265.8), p = 0.779). Since IGF2 also activates IGF1R, thereby contributing to IGF-related bioactivity, IGF2 concentrations were measured in the same serum samples used for IGF1 analysis. Unlike IGF1, IGF2 levels were significantly elevated in conventional husbandry compared to eco husbandry at 105dpc (72.3 ng/ml (IQR 55.5–86.9) vs. 52.3 ng/ml (IQR 49.9–58.7), p = 0.008), with an effect size of 1.38 ± 0.52. At 8dpp, this difference was no longer observed (111.4 ng/ml (IQR 101.2–116.6) vs. 123.1 ng/ml (IQR 102.1–144.9), p = 0.381). Both growth factors exhibited a significant increase in concentration from 105dpc to 8dpp in both housing systems (Table 2, p < 0.001). However, this increase was more pronounced in sows housed in the eco facility, resulting in a convergence of the values at 8dpp (Fig. 4).

Fig. 4. Husbandry dependent IGF concentrations. Concentrations of total IGF1 (A) and IGF2 (B) were determined in serum of differentially housed pregnant and lactating sows (n = 10). Statistical comparisons are based on estimated marginal means. Brackets indicate statistically significant differences with the corresponding p-values shown above.

IGFBP2 but not IGFBP3 is affected by the housing environment

To address the question, if IGF-binding capacities, which control the activity of the IGFs, were affected by the husbandry system, Western ligand blotting was applied in the serum of our mother sow cohort. As expected IGFBP2 and IGFBP3 were the dominant binding proteins in serum, while IGFBP4 reached the detection level only in single replicates, which is why it was omitted for the comparative analysis.

For IGFBP2, housing conditions significantly influenced serum levels (Table 2, p = 0.011). Similar to IGF1 and IGF2 concentrations this effect was more pronounced at 105dpc (conventional: 668.5 ng/ml (IQR 555.8–818.3) vs. ecological: 452.4 ng/ml (IQR 356.2–593.7), p = 0.013), while at 8dpp IGFBP2 serum concentrations measured in sows in the eco facility were nearing those values observed for conventional husbandry (conventional: 593.3 ng/ml (IQR 501.3–784.0) vs. ecological: 573.9 ng/ml (IQR 325.3–619.1)) rendering the numerical difference no longer significant (p = 0.067). Yet, the time of sampling was not found to have a statistically significant effect on serum IGFBP2 levels (Table 2, p = 0.577).

In contrast, IGFBP3 levels were significantly increased in lactating sows compared to pregnant sows (Table 2, p < 0.001), while housing conditions had no significant effect on this parameter (Table 2, p = 0.806). In pregnant sows, IGFBP3 concentrations were 769.4 ng/ml (IQR 478.0–1572.0) under standard conditions and 840.8 ng/ml (IQR 533.9–1408.0) in eco husbandry. In lactating sows, concentrations were markedly higher, measuring 1573.6 ng/ml (IQR 751.8–3335.2) in conventional husbandry and 1954.0 ng/ml (630.1–3682.1) in eco husbandry. Due to this increase in IGFBP3 levels the IGFBP2/IGFBP3 ratio is significantly decreased in lactating sows (Fig. 5A).

Proteolytic activity is minor in serum

Since IGF bioactivity also depends on the proteolytic activity of PAPPA-A1/2, which cleaves binding proteins and thereby releases bioactive IGF, we also assessed proteolytic activity in the serum of the differentially housed sows. For this, we employed a Western ligand blot, detecting the amount of binding protein standards after increasing incubation times with the serum samples. Surprisingly, we observed extremely low proteolysis, which we then quantified by calculating the percentage decrease in detectable binding proteins after 24 h. Even after this extended incubation period, the observed values ranged only from 0% to 21%, with considerable variation between individuals and across different time points within the same individual (Supplementary Table S3). There were no significant differences based on housing conditions or between pregnant and lactating sows (Table 2). The mean values for all conditions fell within a narrow range of 6–7%. Estimates and standard errors of all measured parameters together with P-values for pairwise comparisons within time point and husbandry system are summarized in Supplementary Table S5.

STC1 levels are highly variable but overall lower under conventional husbandry conditions

As another crucial modulator of the overall IGF bioactivity, we investigated STC1 levels in serum of the differentially housed mother sows. For the conventional husbandry system, most measured values were below 200 pg/ml in both pregnant and lactating sows, with only a few samples exceeding 2000 pg/ml (Fig. 5B). A similar bimodal distribution was observed in sows from the eco husbandry system, but at a generally higher level. While some sows had STC1 concentrations below 2000 pg/ml, others exhibited values as high as 8000 pg/ml. Due to the unusual data distribution statistical models cannot be fitted simply. However, model diagnostics demonstrated a decent quality of a linear mixed model, especially regarding the normality of residues. Applying this model, sampling time point (pregnant vs. lactating) as well as husbandry form (ecological vs. conventional) were shown to demonstrate significant effects on STC1 concentrations (Fig. 5B and Supplemental Table S5).

Fig. 5. Husbandry effects on regulators of IGF bioactivity (A) IGFBP2 concentrations were increased in conventional housing resulting in an altered IGFBP2/IGFBP3 ratio (n = 18). (B) STC1 concentrations in serum of pregnant and lactating sows (n = 10). Statistical comparisons are based on estimated marginal means. Brackets indicate statistically significant differences with the corresponding p-values shown above.

Calcium levels are unaffected by the housing environment

Neither housing nor time point of sampling had a significant effect on calcium levels (Table 2), for which we measured values ranging from 2.5 to 26.5 mg/dl. Estimates and standard errors for both husbandry systems and measuring time points are presented in Supplementary Table S5.

Table 2. Analysis of deviance tables (type II Wald chisquare tests) of 10 analyzed parameters. For each factor (housing, time, and their interaction) chi-square values, degree of freedoms (df) and p values are given. Statistically significant factors are in bold.FactorChi-squaredfp-valueCortisol (serum) Housing0.04110.840 Time0.39810.528 Housing×time0.62010.431Cortisol (saliva) Housing0.00510.944 Time0.89710.344 Housing×time2.54610.111IGF bioactivity Housing 5.456 1 0.020  Time1.11010.295 Housing×time0.68410.408IGF1 Housing1.54510.214 Time 37.204 1 < 0.001  Housing×time3.17010.075IGF2 Housing1.555110.212 Time 72.720 1 < 0.001  Housing×time 8.298 1 0.004 IGFBP2 Housing 6.451 1 0.011  Time0.31110.577 Housing×time1.43410.231IGFBP3 Housing0.06010.806 Time 23.658 1 < 0.001  Housing×Time0.03010.863IGFBP2/IGFBP3 Housing3.22010.073 Time 26.410 1 < 0.001  Housing×time0.84110.359Proteolytic activity Housing0.01910.892 Time0.00110.974 Housing×time0.41610.519STC1 Housing 6.124 1 0.013  Time 12.473 1 < 0.001  Housing×time0.30210.583Calcium Housing0.09010.764 Time0.53310.465 Housing×time1.53610.215

Factors of the IGF system partially distinguish husbandry systems

To identify underlying patterns and unsupervised groupings within the data, a PCA was performed based on 14 parameters, including serum IGF bioactivities, saliva cortisol concentrations, total IGF concentrations, binding protein concentrations (IGFBP2 and IGFBP3) and STC1 concentrations for 105 and 8dpp.

As PCA requires a complete dataset without missing values, the initial analysis was conducted on complete cases only, reducing the dataset to 10 cases from the conventional husbandry and 9 cases from the eco facility. While this analysis did not result in completely distinct data clusters, labeling the data points according to the husbandry system revealed that subgroups could be at least partially differentiated (Fig. 6A).

To minimize data loss and improve analytical accuracy, missing values were imputed using KNN imputation. This approach allowed for a more extensive dataset, increasing the sample size to 20 cases per husbandry system. PCA results for the imputed dataset reaffirmed the distinction between subgroups based on the candidate parameters (Fig. 6B). Importantly, variance patterns in the imputed dataset closely resembled those of the original non-imputed data, suggesting that the imputation algorithm preserved subgroup characteristics without introducing bias (Supplementary Fig. S2).

Analysis of feature weights for PC1 and PC2 indicated that serum IGF1 and STC1 concentrations, as well as IGFBP3 abundance accounted for the most variance in the reduced dataset (Fig. 6C). As expected, the principal components of the imputed dataset differed slightly from those of the reduced dataset due to the increased data input, reflected in differences in feature weights for PC1 and PC2 (Fig. 6D). Nevertheless, IGF1 and STC1 levels remained primary drivers of variation, while serum IGF bioactivities and IGFBP2 abundance also substantially contributed to differentiation between the two husbandry systems.

Fig. 6PCA results for candidate parameters determined in the study (A) with a reduced dataset containing only complete cases and (B) with a greater dataset containing imputed data for missing values in the original dataset.

Discussion

Modulations of neuroendocrine factors are part of the animal’s response to its environment. As such external conditions can affect the somatotropic axis, which controls the IGF system. Notably, activity of the IGF system cannot be derived from isolated hormone concentrations, because a number of mechanisms are known to impact on bioavailability and bioactivity of the IGFs (Fig. 1)^25^. In this study, we considered both, the overall IGF-related bioactivity and quantities of individual factors to investigate the impact of the housing environment, and welfare aspects in particular, on the downstream targets of the somatotropic axis in farmed pigs. More generally, we aimed to clarify if the husbandry categories, defined by humans, represent distinct levels of welfare that can be measured objectively. Indeed, we demonstrated a significantly increased serum IGF-dependent bioactivity under conventional husbandry conditions (German husbandry form 1), compared to the eco husbandry (German husbandry form 4/5) in pregnant sows. Moreover, we identified several individual factors of the IGF system to be significantly dependent on the husbandry system such as IGF2, IGFBP2 and STC1 (Table 2).

In contrast, the classical stress marker cortisol failed to reveal significant differences between the husbandry systems when single samplings were conducted at the measuring days. Even though samplings were conducted at similar time points for all animals, physiological fluctuations in cortisol secretion throughout the day and the seasons might contribute to the inconclusive results obtained for this hormone^3^. Similarly inconclusive results were obtained by Meunier-Salaun et al., who report lower cortisol levels in growing pigs both in groups with high and low area allocations compared to medium area allocations^26^. Moreover, our results align with a study of Pearce and Paterson, who conducted repeated blood samplings in male pigs and demonstrated that basal cortisol concentrations are neither elevated by chronic space restriction nor reduced by access to toys^27^. Intriguingly, growth rate was depressed under space restriction implying effects on the GH-IGF-axis.

However, another study in group-housed pregnant gilts by Barnett et al. reports higher plasma cortisol concentrations under space restriction^28^. Of note, the offered space in this study from 1992 was generally extremely restricted compared to our study (0.98 m^2^ or 1.97 m^2^ per gilt vs. 4.0 m^2^ or minimally 6.9 m^2^ per gilt), which might explain these more pronounced differences in cortisol secretion. On the other hand, measured average cortisol concentrations where higher in our study (28.2–31.5 ng/ml) compared to the study of Barnett et al. (converted from nmol/l: 16.5–21.0 ng/ml), in which an implanted cannula was utilized for repeated measures. These elevated serum cortisol concentrations can potentially be attributed to the invasive blood sampling itself, which poses an acute stressor that affects cortisol secretion^29^. For comparison, at 105 dpc normalized cortisol concentrations measured in saliva of the same sows ranged from 11.8 ± 8.4 ng/mg in conventional husbandry to 6.4 ± 4.0 ng/mg in eco husbandry, demonstrating a trend that is in accordance with the results of Barnett et al.

Overall, IGF bioactivities were not significantly affected by the measuring time point (pregnant vs. lactating sows), but concentrations of both IGF1 and IGF2 exhibited a significant increase in concentration from 105 dpc to 8 dpp in both husbandry systems. Increased IGF1 concentrations in lactating sows compared to pregnant sows have been observed by others before^30,31^ and earlier results demonstrating increasing IGF1 and IGF2 serum concentrations throughout lactation^32^ imply a similar course. A study in rats suggests an effect of suckling on IGF1 release^33^, which might be one aspect underlying the findings in our and other studies.

Interestingly, the increase in IGF1 and IGF2 concentrations is more pronounced in the eco facility, which results in diminished differences in serum concentrations between the husbandry systems at 8 dpp. A similar convergence was observed for IGFBP2 even though differences in IGFBP2 levels before and after farrowing did not reach statistical significance. In contrast, we found serum IGFBP3 levels to be significantly increased in lactating sows compared to pregnant sows. To our knowledge, there are no other studies investigating serum IGFBP3 concentrations in pregnant and lactating sows, but several lines of evidence suggest that IGFBP3 might act locally in brain regions responsible for maternal responsiveness and prolactin release following suckling^34^.

Despite the presence of binding proteins in the serum, it might be noted that in some cases (4 out of 40) the IGF bioactivity, determined as IGF1-equivalents in ng/ml based on IGF1 standards, was higher than the sum of the measured total IGF1 and IGF2 concentrations. Several reasons can underly this non-intuitive result. First, due to its structural similarities, insulin can bind the IGF1R with low affinity and thereby activate the downstream signaling cascade, which can, to a limited extent, contribute to the final detected bioactivity^35,36^. Of note, we did not analyze IGF1R levels in target tissues, as the mother sows were not killed for tissue sampling during this experiment. Second, even though IGFBPs are generally known as carriers that inhibit IGF bioactivity until its release, it is now known that IGFBP2 can enhance IGF1 signaling by vimentin mediated binding of receptor protein tyrosine phosphatase β, which eventually results in inhibiting the dephosphorylation of AKT by Phosphatase and Tensin homolog as demonstrated in vascular smooth muscle cells^37^. We confirmed this promoting effect of IGFBP2 on IGF induced AKT phosphorylation in human embryonic kidney cells and tumorigenic epithelial-like cells^14^. In fact, we employ this effect to increase the sensitivity of our BIRA-Assay. In view of this, the elevated IGFBP2 abundance that we observed in the conventional husbandry might contribute to the increased IGF bioactivity that was demonstrated here.

For another crucial modulator of the overall IGF bioactivity, namely STC1, levels seemed to be decreased under standard conditions both in pregnant and lactating sows. Since STC1 inhibits the proteolytic activity of PAPP-A, thereby limiting the release of IGF from its binding proteins (Fig. 1)^12^, this is in line with the increased bioactivity that was observed under standard conditions. However, surprisingly, the proteolytic activity was extremely low in all samples and did not reveal any significant differences between the husbandry systems. Accordingly, there might be another mechanism by which differences in STC1 concentrations translate into differential bioactivities. Literature on the linkages between STC1 and AKT phosphorylation (readout of our bioactivity assay) are inconsistent. While a study by Yang et al. reports that STC1 promotes the phosphorylation of AKT^38^, a recent study of Kim et al. states that STC1 strongly inhibited AKT phosphorylation^39^. Further studies imply that STC1 expression itself is positively^40^ or negatively^41^ regulated by AKT phosphorylation. Apparently, the function and regulation of STC1 is highly tissue and species dependent, rendering further studies in pig crucial to elucidate the relevance of differential serum STC1 concentrations. However, studies involving mammalian STC1 are complicated by the fact that circulating serum STC1 is rapidly eliminated^42^. In fact, STC1 expression predominantly occurs in reproductive tissues of pregnant females, while circulating STC1 is reported to be hardly detectable in nonpregnant individuals^43^. For mice before and after parturition maximum STC1 levels of 1,600 pg/ml were reported^44^, while for cows values ranged between 5,000 and 8,000 pg/ml^45^. Mean values of 8,400 pg/ml were detected in porcine follicular fluids^46^, which is similar to the peak values measured in serum of our mother sow cohort.

Fig. 7. Overview on the differential impact of the husbandry system on the abundance and bioactivity of factors of the IGF system.

While we identified several individual factors of the IGF system to be dependent on the husbandry system (Table 2; Fig. 7), there are still some limitations to their applicability as biomarkers. For example, STC1 concentrations demonstrated very high inter-individual variances, while effects of the husbandry system on IGF bioactivity as well as IGF2, and IGFBP2 levels were only significant in pregnant but not in lactating sows. Nevertheless, our unsupervised approach suggests that factors of the IGF system can be potentially utilized to identify groups based on their housing environment.

The imputation of missing data by bioinformatics methods can further support data analysis by unsupervised approaches. Among available methods, we selected KNN imputation, which is particularly well-suited for small datasets with continuous data and complex biological variation. By this means, we identified STC1 levels and IGF1 concentration as the main contributors to the observed variance in our dataset. This leads us to the conclusion that the IGF system might be well suited to assess differences in general welfare, even if the development of a suitable biomarker could still be a long way off. Moreover, our results imply that husbandry categories for pigs indeed represent distinct levels of welfare that can be measured objectively. However, it needs to be noted that we compared husbandry form 1 and 4, while we believe that we would not have obtained similarly clear results for the comparisons of husbandry form 1 vs. 2 or husbandry form 3 vs. 4.

Surprisingly, IGFBP3 abundance contributed significantly to the principal components that allowed for a distinction of the husbandry systems, despite the fact that differences in IGFBP3 concentrations between the husbandry systems were not found to be statistically significant neither for the reduced nor for the full dataset. Similarly, serum IGF 1 and saliva cortisol concentrations at 105dpc were involved in the observed group separation while not demonstrating statistically significant differences between the two groups. This underpins the notion that measuring individual molecule concentrations might not be sufficient to assess chronic stress and distinguish between husbandry forms. While considering the interplay of multiple factors is one approach to tackle this issue^8,9^, we showcased here that detecting downstream signaling might be another well suited approach. Ultimately, however, only a multifactorial and integrative approach is likely to provide the necessary robustness for the development of biomarkers to objectively evaluate animal welfare.

Concerning the physiological relevance of our findings, a number of studies suggest that deflections of the IGF system are not only a by-product of stress responses that can be exploited for biomarker development but that the system functionally contributes to the animals´ stress coping capacity. Findings across different species indicate that circulating IGF1 is involved in mood homeostasis^47,48^. For example, IGF1 administration was repeatedly shown to initiate cascades of neurochemical effects resulting in antidepressant activity as reviewed by Szczęsny et al.^49^. Similar effects are observed under treatment with the nonspecific IGFBP inhibitor NBI-31,772^50^. Suggested underlying mechanisms include the promotion of other neurotrophic factors such as Brain derived neurotrophic factor^51^, upregulation of serotonin levels in the hippocampus^52^, sensitivity regulation of neuronal glucocorticoid receptors^47^, and modification of GABA/Glutamate synaptic structures in orexin neurons^48^. Accordingly, the increased IGF1 bioactivity that we found in conventional husbandry might reflect the sows’ compensative answer to the more stressful environment.

Moreover, it was reported that overexpression of mutant IGFBP2 reduced anxiety in mice indicating central functions of the IGF system^53^. Intriguingly, we found increased serum IGFBP2 levels in sows housed under conventional husbandry conditions, which is also in line with the hypothesis that factors of the IGF system play a relevant role for stress coping mechanisms.

It needs to be mentioned that the fact that the test animals were only sampled alive and not slaughtered during our longitudinal study leads to a certain limitation. While we monitored all relevant parameters upstream of the IGF1R as well as downstream bioactivity, we could not detect IGF1R expression levels in target tissues. Therefore, we could not asses any regulatory mechanisms on the level of the IGF1 receptor, which are known to be crucial for the translation of signals into pleiotropic biological functions^54^, while differential sensitivity of receptor variants can even add another layer of complexity^55^. For the same reason, it was not possible to detect STC1 in reproductive tissues of investigated sows.

Another potential limitation was introduced by the situation that successfully inseminated sows happened to be younger on average in the ecological facility. This led us to analyze the potential of age as confounding factor. In fact, we observed a significant effect of the animals age on serum concentrations of IGF1, IGFBP2 and STC1, which might imply a habituation process in the sows. Nevertheless, IGF bioactivity under eco conditions was found to be generally lower and these effects were extremely small (Supplementary Table S4). Moreover, this effect is unlikely to be relevant in industrial pig farming, as the time to slaughter is much shorter than the lifetime of the sows examined here. However, further studies will be necessary to assess the compliance of the German husbandry form 4/5 with the animals needs and to elucidate options to continuously improve animal welfare in agriculture.

Conclusions

Our results imply that factors of the somatotropic axis might provide a solid basis for the documentation and objective assessment of animal health and welfare. Potentially, alterations in the IGF system might even indicate illness or chronic stress before symptoms manifest themselves externally in the animal thus allowing for early interventions. However, while we identified several individual factors of the IGF system to be dependent on the husbandry system, there are limitations to their applicability as biomarkers such as high inter-individual variances and changing effect sizes. Ultimately, only a multifactorial and integrative approach is likely to provide the necessary robustness for developing biomarkers to objectively evaluate animal welfare.

Supplementary Information

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Supplementary Material 3

Supplementary Material 4

Supplementary Material 5

Supplementary Material 6

Supplementary Material 7