BACKGROUND
sows have commonly been housed under field conditions in individual stalls throughout pregnancy because it eases animal handling, reduces social stress and allows appropriate feeding.
this individual housing system has been considered to be stressful and harmful for animals by animal welfare experts  <cit>  and consequently, this practice has been banned by the european union .
in gestational stalls, sows have very limited space for moving or laying down  <cit> .
as a consequence, skin abrasions, locomotion difficulties, and loss of muscle mass and bone resistance have been observed .
moreover, abnormal behaviours  <cit>  and reproductive problems  <cit>  are more common in individual-housed sows than those housed in groups.
despite the importance of the objective evaluation of stress, information on laboratorial biomarkers is scarce.
some articles have assessed the validity of phenotypic animal-based measures such as lameness or oral stereotypies to monitor the welfare of sows  <cit> .
several behavioral measures such as sham chewing or excessive drinking suggest that sows in stalls found the conditions less comfortable and reflected a heightened arousal that prevented them from lying down due to stress and frustration  <cit> .
nevertheless, contradictory results have been also observed, since others have found that the drinking frequency was higher in group-housed sows, although an agonistic behavior and sham chewing was higher in stalls  <cit> .
thus, despite the utility of behavioral measures, the availability of a set of objective and easily measurable laboratory parameters that may be combined as a “welfare score” would be an outstanding element to monitor welfare in farm animals.
the complementarity of animal behavior, biochemistry, physiology and immunology together with the use of modern, high-throughput technologies should help to create a more comprehensive scientific basis for animal care and management  <cit> .
previous studies of our group have defined potential laboratory markers to evaluate stress in male pigs housed at different stocking densities  <cit>  that are monitored in gilts housed under stress conditions in the present study.
our first objective was to evaluate the metabolic and hormonal markers under a housing switch in females.
a selection of proteins related to the acute phase response  and oxidative stress pathways were also included in the analysis.
on the other hand, a second objective was to extend the application of gel-based  and gel-free  proteomic technologies to search for potential novel serum markers for stress.
methods
experimental design
six-month old gilts  from picber sa  were used.
the study began in the quarantine facility ), where animals  had been housed for a month, in groups of  <dig>  the following day, all the gilts were moved simultaneously to individual stalls accomplishing all-in and all-out management by productive phase at the seponts s.l.
piglet production farm .
gilts were fed ad libitum in the quarantine facility  or received approximately  <dig>  kg of feed in the individual stalls .
animals were vaccinated against skin diamond disease, swine influenza virus, ausjezky disease and mycoplasma hyopneumoniae with the following commercial vaccines: parvosuin mr, gripork and auskipra , and stellamune mycoplasma , respectively.
the last vaccine was injected a week before the study.
moreover, gilts received ivermectin to control the parasite load  three weeks before the study.table  <dig> calculated composition of the diets for gilts allocated in pen or individual stalls based on its feed composition and its average daily feed in each case
two groups of animals from the original quarantine group  were considered:group ’hs’: a group of animals  were sampled at d <dig>  d <dig>  d <dig> and d <dig>  these animals were subjected to two types of stress: individual housing  and repeated sampling .
group ’h’: three different subgroups of gilts were formed and sampled at d <dig> , d <dig>  or d <dig> .
these animals were only subjected to the stress induced by individual housing  and not to repeated sampling.
this design allowed us to distinguish the effect unleashed by change of housing from the stress induced by repeated sampling, without disturbing the all-in and all-out management system of the farm.
in this experimental design, parameters measured at d <dig> corresponded to basal levels.
the experiment was carried out in a piglet production farm under field conditions before the actual eu legislation against separate housing came into force.
treatment, housing and husbandry conditions conformed to the european union guidelines .
the experiment received approval from the ethical committee for animal experimentation from the universitat de lleida .
sample collection and preparation
gilts were sampled in the morning, from  <dig> to 11 a.m. blood was collected by jugular venipuncture in tubes without anticoagulant  or with edta-k <dig> .
serum was obtained by centrifugation at 2000 g for 10 min.
erythrocyte lysate was prepared after blood centrifugation, as described for serum, and the cellular pellet was washed with  <dig> %  nacl and lysed with cold desionized water.
deproteinized blood was obtained by mixing blood with 5%  metaphosphoric acid followed by centrifugation at 3000 g for 5 min.
saliva samples were collected on salivette™ tubes  by allowing the gilts to chew spontaneously the cotton swabs for 30–60 s without immobilization.
cotton swabs were centrifuged at 4000 g for 10 min.
all centrifugations were at room temperature.
supernatants were frozen at -80 °c until analysis.
analytical parameters
serum non-esterified fatty acids , triglycerides and cholesterol were determined as markers of lipid metabolism.
serum total protein, urea and creatinine were determined as indicators of protein metabolism.
the measurement of metabolic parameters was carried out with the analyzer olympus au <dig> .
technical details are given in table  <dig>  all parameters were analysed in duplicate.table  <dig> list of analytical methods
c-reactive protein a
avalidated for porcine samples  <cit>
immunoblotting
the carbonyl protein content was detected by slot blot and serum apolipoprotein a-i  by immunoblotting, as described by marco-ramell et al.
<cit> , respectively.
a control sample was used in both immunoassays to compare results from different membranes.
ratios between samples and control protein bands were calculated.
proteomic analyses
proteome changes were studied in two parallel pipelines.
discovery proteome studies were performed on pooled serum samples representing d <dig>  d <dig> and d <dig> by 2d-dige and itraq.
2d-differential gel electrophoresis  and spot identification by mass spectrometry
serum samples corresponding to d <dig>  and to d <dig> and d <dig> for both groups  were pooled, obtaining five pools in total , and then desalted and quantified.
fifty μg protein was labelled with cy <dig>  cy <dig> or cy <dig> dyes  and subjected to a bidimensional electrophoresis on four ipg strips  followed by  <dig> % sds-page, as described in  <cit> .
the pooled d <dig> sample was labelled with cy <dig> and used as internal standard in all four gels.
dye-swap was performed and samples were randomly paired, as described in additional file  <dig>  image analysis and statistical quantification of relative protein abundance was performed with progenesis samespots  .
the spot intensity at d <dig> was compared with d <dig> intensity for both animal groups.
spots that met the criteria p <  <dig>  and fold-change  ≥ ± <dig>   were selected for further mass spectrometry  identification.
spots located near the gel borders and small or faint spots were excluded from protein identification.
protein identification and confirmation was performed on an ultraflex tof-tof instrument .
the resulting final peak list was used for protein identification by peptide mass fingerprint.
mascot  <dig>   was used to search the swiss-prot  database.
criterion for positive identification was a significant mascot probability .
when protein identification was not achieved, spots were analyzed by ion trap ms on an esquire hct ultra it mass spectrometer , coupled to a nano-hplc system .
ms/ms fragmentation  was performed on the most intense ions.
ms/ms spectra were searched as:  <dig>  da precursor mass tolerance,  <dig>  da fragment tolerance,  <dig> trypsin missed cleavage, cys carbamidomethylation as fixed modification and met oxidation as variable modification.
a positive identification criterion was set as an individual mascot score for each peptide ms/ms spectrum higher than the corresponding homology threshold score.
full details about the procedures are described in  <cit> .
isobaric tag for relative and absolute quantitation
pooled samples at d <dig>  and both groups  at d <dig> and d <dig>  were used.
pooled samples were enriched with proteominer™  as described in  <cit>  and then desalted and quantified.
proteominer was used to decrease the complexity of serum proteome since itraq requires relatively simple protein mixtures  <cit> .
proteins were reduced, alkylated and digested.
100 μg peptides were differentially labeled with itraq® reagents .
an internal pool, formed by all the samples, was also labeled for further normalization.
peptides were combined and fractionated on a  <dig> offgel fractionator  prior to ms analysis.
<dig>  μg peptides were separated on a nano-uplc acquity system  and analyzed on a ltq-orbitrap velos .
peptides  were analyzed in data dependent mode .
the ten most abundant peptides  from each scan were chosen and fragmented using high energy collision dissociation in a c-trap .
fragments were analyzed at  <dig>  fwhm at 400 m/z.
full ms scan was 250 ms  and ms/ms 300 ms .
raw data was analyzed with xcalibur  <dig> .
<dig>  and protein identification and relative quantification with proteome discoverer .
search parameters were: 10 ppm precursor mass tolerance,  <dig>  da fragment mass tolerance, two missed trypsin cleavages, cys carbamidomethylation and lys and n-term itraq labelling as fixed modifications, met oxidation, n-term acetylation and tyr itraq labelling as variable modifications, mascot score >  <dig> and ‘other mammalia’ taxonomy.
the correction factor for each isobaric label, provided by the manufacturer, was applied and the intensity of each isobaric tag was normalized with the intensity of the internal pool.
the average intensity of all the peptides from the same protein was calculated and normalized by mean-centering, and the ratio between the intensity of each isobaric tag  versus d <dig>  was calculated for each protein.
two criteria were then used to select differential proteins.
first, a fc ≥ ± <dig>  at d <dig> versus d <dig> was required, independently of the variation at d <dig>  then fc were normalized due to the presence of intracellular proteins released to the bloodstream that might affect the overall extracellular proteins levels: only proteins exhibiting a fc greater than the median fc for increased and decreased proteins were selected.
identifications with less than two unique peptides were excluded from the selection.
statistical analysis of biochemical data
all statistical analyses were carried out using the sas system v. <dig> .
<dig> .
for all biochemical analyses, the individual pig was used as the experimental unit.
when not specified, the significance level  was set at  <dig> .
a test to discard outliers was not performed.
the variables included in the statistical analyses were continuous.
shapiro-wilk and levene tests were used to evaluate the normality of the distribution of the variables and the homogeneity of variances, respectively.
a non-parametric  was chosen to compare the results obtained for all the measured parameters between d <dig>  and the rest of the days for the h and hs groups.
the effect of repeated sampling on the parameters was carried out using a non-parametric  between h and hs groups.
RESULTS
general biochemical profiling
the metabolic and the stress response after transition to individual stalls was evaluated by measuring several blood parameters.
results at d <dig>  d <dig> and d <dig> were always compared with basal levels  in h and hs groups.
differences between h and hs groups are also described.
all statistical differences are reported in table  <dig> and figs.
<dig> and  <dig> table  <dig> serum concentration of compounds related to lipid and protein metabolism, enzyme activities and cortisol in serum and saliva in the hs and the h group throughout the experimental period
<dig>  ±  <dig> **
138 ± 46*
<dig>  ±  <dig> *
<dig>  ±  <dig> *
<dig>  ±  <dig> *a
<dig>  ±  <dig> *
a,b within columns, means with a different superscript differed significantly  between h and hs group
asterisks show significant differences versus d1: *p <  <dig> , **p <  <dig> , ***p <  <dig>
fig.
<dig> average serum concentrations of acute phase proteins in h and hs groups throughout the trial.
asterisks represent statistically significant differences between d <dig>  d <dig> or d <dig> and basal levels : *p <  <dig> , **p <  <dig> , ***p <  <dig> .
different letter means significant differences between h and hs groups
fig.
<dig> average serum concentrations of oxidative stress markers in h and hs groups throughout the trial.
asterisks represent statistically significant differences between d <dig>  d <dig> or d <dig> and basal levels : *p <  <dig> , **p <  <dig> , ***p <  <dig> .
different letter means significant differences between h and hs groups
a transient increase in non-esterified fatty acids  was observed at d <dig> in both groups that was rapidly normalized at d <dig>  triglycerides  significantly increased at the end of the study  in both groups.
no changes were observed in total cholesterol except an increase at d <dig> in the hs group.
no differences were observed between hs and h groups for any lipid marker with the exception of nefa at d <dig>  which were lower in the h group.
urea decreased at d <dig> and d <dig> in both groups.
creatinine increased at d <dig> in both groups, and at d <dig> and d <dig> only in h group.
serum total protein concentration increased after the transition to individual stalls at d <dig> and d <dig>  no differences between hs and h groups were observed, except for creatinine and total protein at d <dig>
creatine kinase , alanine aminotransferase  and alkaline phosphatase  enzymatic activities were measured in serum as markers of muscular  or hepatic damage .
a high inter-individual variability in ck and alt activities was found, especially in the hs group at d <dig>  an increase in ck at d <dig> was observed in the hs group.
cortisol was determined in serum and saliva as stress marker.
a great inter-individual variability in this parameter was observed, especially in saliva.
no significant effect was observed due either to the change of housing or to the repeated sampling procedure.
innate immune system and oxidative stress defences monitoring
porcine app levels showed changes throughout the study and different kinetics between h and hs groups .
haptoglobin  slightly increased at d <dig> and normalized afterwards.
c-reactive protein  increased at d <dig> and onwards only in the hs group, whereas apo a-i only decreased at d <dig> in the same animal group.
pig-map decreased at d <dig> and onwards only in the h group.
differences between both groups were observed in hp and apo a-i at d <dig>  and throughout all the study  in crp and pig-map levels.
oxidative stress markers also showed different kinetics after the housing change and between h and hs groups .
in the h group, glutathione peroxidase activity  increased at d <dig>  superoxide dismutase  increased at d <dig>  and slight increases were observed in tgsh and protein-carbonyl groups.
significant differences between hs and h groups were only observed for gpx at d <dig> and d <dig>  and sod at d <dig>
proteomic analysis
proteome changes were analyzed by 2d-dige and itraq.
table  <dig> shows differentially expressed proteins at d <dig> and d <dig> versus d <dig> in the h and hs groups, identified by dige or itraq.
the complete list of identified proteins is shown in additional file  <dig>  and additional file  <dig> .table  <dig> selected differentially expressed proteins at d <dig> and d <dig> versus d <dig>  in the h and hs groups identified by proteomics
increased proteins at d <dig> are represented as positive fold-changes  and decreased proteins as negative fc.
proteins are classified according to their biological function.
the complete list of differentially expressed proteins identified by ms is shown in additional file  <dig>  and additional file  <dig>
2d-dige
a subset of  <dig> spots was found in all four bidimensional gels.
of these,  <dig> spots were differentially expressed.
ten of these were selected, as described in the methods section, but only five of them could be identified by ms. these identifications include the app hp , apo a-i and α1-antichymotrypsin  <dig> and the antioxidant enzyme peroxiredoxin  <dig>  .
hp, α1-antichymotrypsin  <dig> and prdx <dig> significantly increased whereas apo a-i significantly decreased at d <dig> in h and hs groups.
all four proteins tended to reach basal levels at d <dig> in both groups .fig.
<dig> serum proteome analysis by dige.
representative image of one of the four 2d gels containing 150 μg protein  and stained with flamingo® .
differentially expressed spots identified by ms are marked
itraq
the 8-plex itraq® analysis yielded the identification of  <dig> proteins.
the overall analysis of the data reflects a slight increase in intracellular proteins released to the bloodstream and a relative decrease in many extracellular proteins .
to counteract this protein bias, itraq data was normalized as described in the methods section.
sixty-seven decreased and  <dig> increased proteins at d <dig> met these two criteria .
most proteins showed larger fold-change variations over time in the hs than the h group, and larger fold-change variations at d <dig> than at d <dig>
amongst the differentially expressed proteins, two of them were common with the 2d-dige analysis: apoa-i was found to be decreased and prdx <dig> increased at d <dig> and d <dig>  other functionally-related proteins such as apolipoproteins c-iii, e, m and r  and glutathione peroxidase  isoform  <dig>  were also identified by itraq.
another group of differentially expressed proteins with a relevant function was related to the immune system, including lactotransferrin, protegrin  <dig> and galectin  <dig>  which increased at d <dig> in both h and hs animals, and tended to basal levels at d <dig>  finally, another relevant protein group were intracellular proteins such as actin-binding proteins, tubulins, enzymes involved in glycolysis and several proteasome subunits, which were found to be increased in serum .
some of these findings were already observed in the biochemical analysis of all the samples .
discussion
the allocation of sows in isolated farrowing stalls is believed to be a stressful procedure and the change of environmental conditions is a challenge to the individual, that must adapt to the new situation.
several problems and abnormal behaviours have been reported in individual-housed sows, including locomotion difficulties, skin lesions and loss of muscle mass and bone resistance, together with behavioural symptoms .
in this study, hs gilts were subjected to two parallel stressors, namely individual housing/diet  and repeated blood sampling  the day before  and at the three following days post stressor , whereas h gilts only experienced a switch of housing/diet.
a detailed biochemical profile was performed in all the samples  and two proteomics approaches were performed for the identification of new potential markers of this stress condition.
all the biochemical parameters analyzed were within the reference range indicating that gilts were healthy during all the days of the experiment.
nefa levels indicated an increase of lipid mobilization from fat stores at d <dig>  whereas the increase of tg at the end of the study in both animal groups indicated that gilts adapted to the new environment and normalized their food intake.
cholesterol was elevated at short-time in the hs group, similarly to high-density housed pigs in our previous study  <cit> .
variations in lipid metabolism were also observed in the proteomic analyses, as several apolipoproteins involved in lipid transport were found decreased both by itraq and 2d-dige studies.
these observations confirm that altered lipid metabolism is associated to physiological stress, likely as a consequence of the lipolytic activity of epinephrine and cortisol  <cit> .
stress, change of diet, and low food intake might also have an accelerating effect on tissue protein catabolism, as reflected by urea and creatinine variation.
ck, a marker for skeletal muscle, increased in hs gilts probably as a consequence of tissue damage during blood extraction  <cit> .
a strong correlation between ck and alt was observed , being both enzymes highly abundant in skeletal muscle.
the large individual variation in ck and alt levels suggests that blood sampling was not equally harmful for all the animals.
the measurement of salivary cortisol is recommended due to the simplicity of sampling and direct biological relevance  <cit>  and salivary cortisol levels have previously been found to increase in group-housed sows compared with stalled sows  <cit> .
in the present study, serum and salivary cortisol were unaffected, but this latter presented high individual variation, suggesting that the perception of stress was variable between animals  <cit> .
our previous studies in pigs identified pig-map increased by high-density housing, whereas other app were not altered  <cit> .
in the present situation, hp increased and apo a-i decreased confirming the role of app as stress markers .
pig-map concentrations were unexpectedly high at the beginning of the experiment, which could have masked the expected increase of this protein.
the high mean value was due to five out of the  <dig> individuals having a pig-map concentration between  <dig>  and  <dig>  mg/ml but with normal values for the other app, ruling out the possibility of these individuals being ill. other authors have reported that hp and pig-map behave similar in sows kept in stalls or in group but in that case animals were maintained in the same housing system all time, without changing conditions  <cit> .
on the other hand, crp profile was remarkably different between the hs and h groups, suggesting that it may be a good marker of inflammation caused by repeated blood sampling and not affected by housing.
regarding oxidative stress, gilts of the h group had higher activity of the antioxidant defence pathways than gilts of the hs group.
repetitive blood sampling might have increased the degree of oxidative stress in the hs gilts by reducing the antioxidant defences or raising the reactive oxygen species  production .
these results were similar to our previous studies  <cit> .
two different proteomic approaches, 2d-dige and itraq, revealed larger variations in protein composition in the hs group than in the h group, confirming that gilts repeatedly sampled and subjected to the change of housing/diet were exposed to a stronger challenge than gilts only subjected to this latter situation.
working with pooled samples has some limitations as it can reduce the effects of biological variation.
nevertheless, pooling does not mask the biological variation itself and it has been shown to be beneficial when working with a high number of individuals  <cit> .
only two proteins, apo a-i and prdx <dig>  were identified as differential proteins in both procedures, probably as a consequence of the complementary nature of the two techniques  <cit> , the lower fold-changes achieved by itraq  <cit>  and the use of different serum proteomes .
nevertheless, both methods confirmed the innate immune system and the redox defence pathways as the two main homeostatic mechanisms involved in the stress response.
thus, hp increased after the housing switch, confirming previous results.
other altered proteins related to the immune defences were α1-antichymotrypsin  <dig>  lactotransferrin, ficolin- <dig>  protegrin  <dig> and galectin  <dig>  protegrins are short cathelicidins exclusively produced by leukocytes  <cit> , which are upregulated by proinflammatory cytokines and apps, and may be proposed as indicators of the innate immune response in the pig  <cit> .
the immune profile of sows housed in stalls is different from those housed in a greater floor space  <cit> , and it has been suggested that it may reflect a physiological stress response that enabled them to adapt to their environment.
all together, these results suggest that an immunological response is required to reestablish homeostasis.
the involvement of oxidative processes was also confirmed by proteomics, since several antioxidant enzymes were found differentially expressed i.e.
prdx <dig> increased, whereas gpx isoform  <dig>  decreased.
gpx activity was found increased in erythrocyte lysates, but the involvement of different gpx isoforms could explain the apparent contradiction between itraq and enzymatic results.
proteomics also indicated that a certain degree of cellular damage has occurred, entailing a relative increase of cellular proteins in plasma, such as actin-binding proteins, tubulins, enzymes involved in glycolysis and several proteasome subunits.
some of these proteins are potentially good markers, since their increase in serum is large.
this is the case of tropomyosin-beta, glycogen phosphorylase and glyceraldehyde-3-phosphate dehydrogenase .
cellular leakage could be increased due to physical injuries linked to the reduced space of the individual stall and to sample extraction.
it is interesting to note that the cytoplasmatic protein β-actin has been previously validated as a stress biomarker in pigs housed at a high density  <cit> .
nevertheless, a thorough validation of these proteins in a higher number of individuals and in different stress conditions should be performed before they could be considered adequate biomarkers.
in all the cases  the serum analytical variations observed in the h and hs group took place in the same direction.
global changes were always larger in the hs than in the h group, reinforcing the idea that damage was more important in gilts subjected to both types of stressors.
further research is needed to confirm the potential role of the differentially expressed proteins as stress biomarkers and to characterize their specificity, sensitivity and kinetics.
CONCLUSIONS
the confinement of gilts in individual stalls promotes relevant changes in serum composition including an increase of lipid, protein catabolism and oxidative stress, the activation of the innate immune system and may cause tissue damage.
therefore we propose that a combination of plasma lipid markers, app or other innate immune indicators, redox components, intracellular proteins and salivary cortisol may be proposed as a welfare laboratory profile for its application in the pig routine.
additional files
additional file 1: experimental design of the 2d-dige experiment .
the pooled d <dig> sample was labelled with cy <dig> and used as internal standard in all four gels.
dye-swap was performed and samples were randomly paired.
additional file 2: complete list of differentially expressed proteins identified by dige at d <dig> and d <dig> versus d <dig>  in h and hs groups .
increased proteins at d <dig> are represented as positive fold-changes  and decreased proteins as negative fc.
proteins are ordered by number of spot.
additional file 3: list of differentially expressed proteins identified and quantified by itraq at d <dig> and d <dig> versus d <dig>  in h and hs groups .
increased proteins at d <dig> are represented as positive fold-changes  and decreased proteins as negative fc.
proteins are classified according to their biological function and ordered by mascot score.
abbreviations
2dbidimensional electrophoresis
alpalkaline phosphatase
altalanine aminotransferase
apo a-iapolipoprotein a-i
appacute phase protein
ckcreatine kinase
crpc-reactive protein
digedifferential gel electrophoresis
fcfold change
gpxglutathione peroxidase
hphaptoglobin
itraqisobaric tag for relative and absolute quantitation
msmass spectrometry
nefanon-esterified fatty acids
pig-mappig major acute phase protein
prdx2peroxiredoxin 2
rosreactive oxygen species
sodsuperoxide dismutase
tgtriglycerides
tgshtotal glutathione

 sows have commonly been housed under field conditions in individual stalls throughout pregnancy because it eases animal handling, reduces social stress and allows appropriate feeding. this individual housing system has been considered to be stressful and harmful for animals by animal welfare experts and consequently, this practice has been banned by the european union. in gestational stalls, sows have very limited space for moving or laying down. as a consequence, skin abrasions, locomotion difficulties, and loss of muscle mass and bone resistance have been observed. moreover, abnormal behaviours and reproductive problems are more common in individual-housed sows than those housed in groups. despite the importance of the objective evaluation of stress, information on laboratorial biomarkers is scarce. some articles have assessed the validity of phenotypic animal-based measures such as lameness or oral stereotypies to monitor the welfare of sows. several behavioral measures such as sham chewing or excessive drinking suggest that sows in stalls found the conditions less comfortable and reflected a heightened arousal that prevented them from lying down due to stress and frustration. nevertheless, contradictory results have been also observed, since others have found that the drinking frequency was higher in group-housed sows, although an agonistic behavior and sham chewing was higher in stalls. thus, despite the utility of behavioral measures, the availability of a set of objective and easily measurable laboratory parameters that may be combined as a “welfare score” would be an outstanding element to monitor welfare in farm animals. the complementarity of animal behavior, biochemistry, physiology and immunology together with the use of modern, high-throughput technologies should help to create a more comprehensive scientific basis for animal care and management. previous studies of our group have defined potential laboratory markers to evaluate stress in male pigs housed at different stocking densities that are monitored in gilts housed under stress conditions in the present study. our first objective was to evaluate the metabolic and hormonal markers under a housing switch in females. a selection of proteins related to the acute phase response and oxidative stress pathways were also included in the analysis. on the other hand, a second objective was to extend the application of gel-based and gel-free proteomic technologies to search for potential novel serum markers for stress. methods experimental design six-month old gilts from picber sa were used. the study began in the quarantine facility ), where animals had been housed for a month, in groups of the following day, all the gilts were moved simultaneously to individual stalls accomplishing all-in and all-out management by productive phase at the seponts s.l. piglet production farm. gilts were fed ad libitum in the quarantine facility or received approximately  kg of feed in the individual stalls. animals were vaccinated against skin diamond disease, swine influenza virus, ausjezky disease and mycoplasma hyopneumoniae with the following commercial vaccines: parvosuin mr, gripork and auskipra, and stellamune mycoplasma, respectively. the last vaccine was injected a week before the study. moreover, gilts received ivermectin to control the parasite load three weeks before the study.table calculated composition of the diets for gilts allocated in pen or individual stalls based on its feed composition and its average daily feed in each case two groups of animals from the original quarantine group were considered:group ’hs’: a group of animals were sampled at d d d and d these animals were subjected to two types of stress: individual housing and repeated sampling. group ’h’: three different subgroups of gilts were formed and sampled at d, d or d. these animals were only subjected to the stress induced by individual housing and not to repeated sampling. this design allowed us to distinguish the effect unleashed by change of housing from the stress induced by repeated sampling, without disturbing the all-in and all-out management system of the farm. in this experimental design, parameters measured at d corresponded to basal levels. the experiment was carried out in a piglet production farm under field conditions before the actual eu legislation against separate housing came into force. treatment, housing and husbandry conditions conformed to the european union guidelines. the experiment received approval from the ethical committee for animal experimentation from the universitat de lleida. sample collection and preparation gilts were sampled in the morning, from to 11 a.m. blood was collected by jugular venipuncture in tubes without anticoagulant or with edta-k. serum was obtained by centrifugation at 2000 g for 10 min. erythrocyte lysate was prepared after blood centrifugation, as described for serum, and the cellular pellet was washed with % nacl and lysed with cold desionized water. deproteinized blood was obtained by mixing blood with 5% metaphosphoric acid followed by centrifugation at 3000 g for 5 min. saliva samples were collected on salivette™ tubes by allowing the gilts to chew spontaneously the cotton swabs for 30–60 s without immobilization. cotton swabs were centrifuged at 4000 g for 10 min. all centrifugations were at room temperature. supernatants were frozen at -80 °c until analysis. analytical parameters serum non-esterified fatty acids, triglycerides and cholesterol were determined as markers of lipid metabolism. serum total protein, urea and creatinine were determined as indicators of protein metabolism. the measurement of metabolic parameters was carried out with the analyzer olympus au. technical details are given in table  all parameters were analysed in duplicate.table list of analytical methods c-reactive protein a avalidated for porcine samples immunoblotting the carbonyl protein content was detected by slot blot and serum apolipoprotein a-i by immunoblotting, as described by marco-ramell et al., respectively. a control sample was used in both immunoassays to compare results from different membranes. ratios between samples and control protein bands were calculated. proteomic analyses proteome changes were studied in two parallel pipelines. discovery proteome studies were performed on pooled serum samples representing d d and d by 2d-dige and itraq. 2d-differential gel electrophoresis and spot identification by mass spectrometry serum samples corresponding to d and to d and d for both groups were pooled, obtaining five pools in total, and then desalted and quantified. fifty μg protein was labelled with cy cy or cy dyes and subjected to a bidimensional electrophoresis on four ipg strips followed by % sds-page, as described in. the pooled d sample was labelled with cy and used as internal standard in all four gels. dye-swap was performed and samples were randomly paired, as described in additional file image analysis and statistical quantification of relative protein abundance was performed with progenesis samespots. the spot intensity at d was compared with d intensity for both animal groups. spots that met the criteria p <  and fold-change  ≥ ± were selected for further mass spectrometry identification. spots located near the gel borders and small or faint spots were excluded from protein identification. protein identification and confirmation was performed on an ultraflex tof-tof instrument. the resulting final peak list was used for protein identification by peptide mass fingerprint. mascot was used to search the swiss-prot database. criterion for positive identification was a significant mascot probability. when protein identification was not achieved, spots were analyzed by ion trap ms on an esquire hct ultra it mass spectrometer, coupled to a nano-hplc system. ms/ms fragmentation was performed on the most intense ions. ms/ms spectra were searched as:  da precursor mass tolerance,  da fragment tolerance, trypsin missed cleavage, cys carbamidomethylation as fixed modification and met oxidation as variable modification. a positive identification criterion was set as an individual mascot score for each peptide ms/ms spectrum higher than the corresponding homology threshold score. full details about the procedures are described in. isobaric tag for relative and absolute quantitation pooled samples at d and both groups at d and d were used. pooled samples were enriched with proteominer™ as described in and then desalted and quantified. proteominer was used to decrease the complexity of serum proteome since itraq requires relatively simple protein mixtures. proteins were reduced, alkylated and digested. 100 μg peptides were differentially labeled with itraq® reagents. an internal pool, formed by all the samples, was also labeled for further normalization. peptides were combined and fractionated on a offgel fractionator prior to ms analysis.  μg peptides were separated on a nano-uplc acquity system and analyzed on a ltq-orbitrap velos. peptides were analyzed in data dependent mode. the ten most abundant peptides from each scan were chosen and fragmented using high energy collision dissociation in a c-trap. fragments were analyzed at fwhm at 400 m/z. full ms scan was 250 ms and ms/ms 300 ms. raw data was analyzed with xcalibur. and protein identification and relative quantification with proteome discoverer. search parameters were: 10 ppm precursor mass tolerance,  da fragment mass tolerance, two missed trypsin cleavages, cys carbamidomethylation and lys and n-term itraq labelling as fixed modifications, met oxidation, n-term acetylation and tyr itraq labelling as variable modifications, mascot score >  and ‘other mammalia’ taxonomy. the correction factor for each isobaric label, provided by the manufacturer, was applied and the intensity of each isobaric tag was normalized with the intensity of the internal pool. the average intensity of all the peptides from the same protein was calculated and normalized by mean-centering, and the ratio between the intensity of each isobaric tag versus d was calculated for each protein. two criteria were then used to select differential proteins. first, a fc ≥ ± at d versus d was required, independently of the variation at d then fc were normalized due to the presence of intracellular proteins released to the bloodstream that might affect the overall extracellular proteins levels: only proteins exhibiting a fc greater than the median fc for increased and decreased proteins were selected. identifications with less than two unique peptides were excluded from the selection. statistical analysis of biochemical data all statistical analyses were carried out using the sas system v... for all biochemical analyses, the individual pig was used as the experimental unit. when not specified, the significance level was set at. a test to discard outliers was not performed. the variables included in the statistical analyses were continuous. shapiro-wilk and levene tests were used to evaluate the normality of the distribution of the variables and the homogeneity of variances, respectively. a non-parametric was chosen to compare the results obtained for all the measured parameters between d and the rest of the days for the h and hs groups. the effect of repeated sampling on the parameters was carried out using a non-parametric between h and hs groups.  general biochemical profiling the metabolic and the stress response after transition to individual stalls was evaluated by measuring several blood parameters. results at d d and d were always compared with basal levels in h and hs groups. differences between h and hs groups are also described. all statistical differences are reported in table  and figs.  and table serum concentration of compounds related to lipid and protein metabolism, enzyme activities and cortisol in serum and saliva in the hs and the h group throughout the experimental period  ±  ** 138 ± 46*  ±  *  ±  *  ±  *a  ±  * a,b within columns, means with a different superscript differed significantly between h and hs group asterisks show significant differences versus d1: *p < , **p < , ***p <  fig. average serum concentrations of acute phase proteins in h and hs groups throughout the trial. asterisks represent statistically significant differences between d d or d and basal levels: *p < , **p < , ***p < . different letter means significant differences between h and hs groups fig. average serum concentrations of oxidative stress markers in h and hs groups throughout the trial. asterisks represent statistically significant differences between d d or d and basal levels: *p < , **p < , ***p < . different letter means significant differences between h and hs groups a transient increase in non-esterified fatty acids was observed at d in both groups that was rapidly normalized at d triglycerides significantly increased at the end of the study in both groups. no changes were observed in total cholesterol except an increase at d in the hs group. no differences were observed between hs and h groups for any lipid marker with the exception of nefa at d which were lower in the h group. urea decreased at d and d in both groups. creatinine increased at d in both groups, and at d and d only in h group. serum total protein concentration increased after the transition to individual stalls at d and d no differences between hs and h groups were observed, except for creatinine and total protein at d creatine kinase, alanine aminotransferase and alkaline phosphatase enzymatic activities were measured in serum as markers of muscular or hepatic damage. a high inter-individual variability in ck and alt activities was found, especially in the hs group at d an increase in ck at d was observed in the hs group. cortisol was determined in serum and saliva as stress marker. a great inter-individual variability in this parameter was observed, especially in saliva. no significant effect was observed due either to the change of housing or to the repeated sampling procedure. innate immune system and oxidative stress defences monitoring porcine app levels showed changes throughout the study and different kinetics between h and hs groups. haptoglobin slightly increased at d and normalized afterwards. c-reactive protein increased at d and onwards only in the hs group, whereas apo a-i only decreased at d in the same animal group. pig-map decreased at d and onwards only in the h group. differences between both groups were observed in hp and apo a-i at d and throughout all the study in crp and pig-map levels. oxidative stress markers also showed different kinetics after the housing change and between h and hs groups. in the h group, glutathione peroxidase activity increased at d superoxide dismutase increased at d and slight increases were observed in tgsh and protein-carbonyl groups. significant differences between hs and h groups were only observed for gpx at d and d and sod at d proteomic analysis proteome changes were analyzed by 2d-dige and itraq. table  shows differentially expressed proteins at d and d versus d in the h and hs groups, identified by dige or itraq. the complete list of identified proteins is shown in additional file and additional file.table selected differentially expressed proteins at d and d versus d in the h and hs groups identified by proteomics increased proteins at d are represented as positive fold-changes and decreased proteins as negative fc. proteins are classified according to their biological function. the complete list of differentially expressed proteins identified by ms is shown in additional file and additional file 2d-dige a subset of spots was found in all four bidimensional gels. of these, spots were differentially expressed. ten of these were selected, as described in the methods section, but only five of them could be identified by ms. these identifications include the app hp, apo a-i and α1-antichymotrypsin and the antioxidant enzyme peroxiredoxin. hp, α1-antichymotrypsin and prdx significantly increased whereas apo a-i significantly decreased at d in h and hs groups. all four proteins tended to reach basal levels at d in both groups.fig. serum proteome analysis by dige. representative image of one of the four 2d gels containing 150 μg protein and stained with flamingo®. differentially expressed spots identified by ms are marked itraq the 8-plex itraq® analysis yielded the identification of proteins. the overall analysis of the data reflects a slight increase in intracellular proteins released to the bloodstream and a relative decrease in many extracellular proteins. to counteract this protein bias, itraq data was normalized as described in the methods section. sixty-seven decreased and increased proteins at d met these two criteria. most proteins showed larger fold-change variations over time in the hs than the h group, and larger fold-change variations at d than at d amongst the differentially expressed proteins, two of them were common with the 2d-dige analysis: apoa-i was found to be decreased and prdx increased at d and d other functionally-related proteins such as apolipoproteins c-iii, e, m and r and glutathione peroxidase isoform were also identified by itraq. another group of differentially expressed proteins with a relevant function was related to the immune system, including lactotransferrin, protegrin and galectin which increased at d in both h and hs animals, and tended to basal levels at d finally, another relevant protein group were intracellular proteins such as actin-binding proteins, tubulins, enzymes involved in glycolysis and several proteasome subunits, which were found to be increased in serum. some of these findings were already observed in the biochemical analysis of all the samples. discussion the allocation of sows in isolated farrowing stalls is believed to be a stressful procedure and the change of environmental conditions is a challenge to the individual, that must adapt to the new situation. several problems and abnormal behaviours have been reported in individual-housed sows, including locomotion difficulties, skin lesions and loss of muscle mass and bone resistance, together with behavioural symptoms. in this study, hs gilts were subjected to two parallel stressors, namely individual housing/diet and repeated blood sampling the day before and at the three following days post stressor, whereas h gilts only experienced a switch of housing/diet. a detailed biochemical profile was performed in all the samples and two proteomics approaches were performed for the identification of new potential markers of this stress condition. all the biochemical parameters analyzed were within the reference range indicating that gilts were healthy during all the days of the experiment. nefa levels indicated an increase of lipid mobilization from fat stores at d whereas the increase of tg at the end of the study in both animal groups indicated that gilts adapted to the new environment and normalized their food intake. cholesterol was elevated at short-time in the hs group, similarly to high-density housed pigs in our previous study. variations in lipid metabolism were also observed in the proteomic analyses, as several apolipoproteins involved in lipid transport were found decreased both by itraq and 2d-dige studies. these observations confirm that altered lipid metabolism is associated to physiological stress, likely as a consequence of the lipolytic activity of epinephrine and cortisol. stress, change of diet, and low food intake might also have an accelerating effect on tissue protein catabolism, as reflected by urea and creatinine variation. ck, a marker for skeletal muscle, increased in hs gilts probably as a consequence of tissue damage during blood extraction. a strong correlation between ck and alt was observed, being both enzymes highly abundant in skeletal muscle. the large individual variation in ck and alt levels suggests that blood sampling was not equally harmful for all the animals. the measurement of salivary cortisol is recommended due to the simplicity of sampling and direct biological relevance and salivary cortisol levels have previously been found to increase in group-housed sows compared with stalled sows. in the present study, serum and salivary cortisol were unaffected, but this latter presented high individual variation, suggesting that the perception of stress was variable between animals. our previous studies in pigs identified pig-map increased by high-density housing, whereas other app were not altered. in the present situation, hp increased and apo a-i decreased confirming the role of app as stress markers. pig-map concentrations were unexpectedly high at the beginning of the experiment, which could have masked the expected increase of this protein. the high mean value was due to five out of the individuals having a pig-map concentration between and  mg/ml but with normal values for the other app, ruling out the possibility of these individuals being ill. other authors have reported that hp and pig-map behave similar in sows kept in stalls or in group but in that case animals were maintained in the same housing system all time, without changing conditions. on the other hand, crp profile was remarkably different between the hs and h groups, suggesting that it may be a good marker of inflammation caused by repeated blood sampling and not affected by housing. regarding oxidative stress, gilts of the h group had higher activity of the antioxidant defence pathways than gilts of the hs group. repetitive blood sampling might have increased the degree of oxidative stress in the hs gilts by reducing the antioxidant defences or raising the reactive oxygen species production. these results were similar to our previous studies. two different proteomic approaches, 2d-dige and itraq, revealed larger variations in protein composition in the hs group than in the h group, confirming that gilts repeatedly sampled and subjected to the change of housing/diet were exposed to a stronger challenge than gilts only subjected to this latter situation. working with pooled samples has some limitations as it can reduce the effects of biological variation. nevertheless, pooling does not mask the biological variation itself and it has been shown to be beneficial when working with a high number of individuals. only two proteins, apo a-i and prdx were identified as differential proteins in both procedures, probably as a consequence of the complementary nature of the two techniques, the lower fold-changes achieved by itraq and the use of different serum proteomes. nevertheless, both methods confirmed the innate immune system and the redox defence pathways as the two main homeostatic mechanisms involved in the stress response. thus, hp increased after the housing switch, confirming previous results. other altered proteins related to the immune defences were α1-antichymotrypsin lactotransferrin, ficolin- protegrin and galectin protegrins are short cathelicidins exclusively produced by leukocytes, which are upregulated by proinflammatory cytokines and apps, and may be proposed as indicators of the innate immune response in the pig. the immune profile of sows housed in stalls is different from those housed in a greater floor space, and it has been suggested that it may reflect a physiological stress response that enabled them to adapt to their environment. all together, these results suggest that an immunological response is required to reestablish homeostasis. the involvement of oxidative processes was also confirmed by proteomics, since several antioxidant enzymes were found differentially expressed i.e. prdx increased, whereas gpx isoform decreased. gpx activity was found increased in erythrocyte lysates, but the involvement of different gpx isoforms could explain the apparent contradiction between itraq and enzymatic results. proteomics also indicated that a certain degree of cellular damage has occurred, entailing a relative increase of cellular proteins in plasma, such as actin-binding proteins, tubulins, enzymes involved in glycolysis and several proteasome subunits. some of these proteins are potentially good markers, since their increase in serum is large. this is the case of tropomyosin-beta, glycogen phosphorylase and glyceraldehyde-3-phosphate dehydrogenase. cellular leakage could be increased due to physical injuries linked to the reduced space of the individual stall and to sample extraction. it is interesting to note that the cytoplasmatic protein β-actin has been previously validated as a stress biomarker in pigs housed at a high density. nevertheless, a thorough validation of these proteins in a higher number of individuals and in different stress conditions should be performed before they could be considered adequate biomarkers. in all the cases the serum analytical variations observed in the h and hs group took place in the same direction. global changes were always larger in the hs than in the h group, reinforcing the idea that damage was more important in gilts subjected to both types of stressors. further research is needed to confirm the potential role of the differentially expressed proteins as stress biomarkers and to characterize their specificity, sensitivity and kinetics.  the confinement of gilts in individual stalls promotes relevant changes in serum composition including an increase of lipid, protein catabolism and oxidative stress, the activation of the innate immune system and may cause tissue damage. therefore we propose that a combination of plasma lipid markers, app or other innate immune indicators, redox components, intracellular proteins and salivary cortisol may be proposed as a welfare laboratory profile for its application in the pig routine. additional files additional file 1: experimental design of the 2d-dige experiment. the pooled d sample was labelled with cy and used as internal standard in all four gels. dye-swap was performed and samples were randomly paired. additional file 2: complete list of differentially expressed proteins identified by dige at d and d versus d in h and hs groups. increased proteins at d are represented as positive fold-changes and decreased proteins as negative fc. proteins are ordered by number of spot. additional file 3: list of differentially expressed proteins identified and quantified by itraq at d and d versus d in h and hs groups. increased proteins at d are represented as positive fold-changes and decreased proteins as negative fc. proteins are classified according to their biological function and ordered by mascot score. abbreviations 2dbidimensional electrophoresis alpalkaline phosphatase altalanine aminotransferase apo a-iapolipoprotein a-i appacute phase protein ckcreatine kinase crpc-reactive protein digedifferential gel electrophoresis fcfold change gpxglutathione peroxidase hphaptoglobin itraqisobaric tag for relative and absolute quantitation msmass spectrometry nefanon-esterified fatty acids pig-mappig major acute phase protein prdx2peroxiredoxin 2 rosreactive oxygen species sodsuperoxide dismutase tgtriglycerides tgshtotal glutathione 
biochemical analyses showed adaptive changes in lipid and protein metabolism after the housing switch, whereas cortisol and muscular markers showed a large variability between animals. housing of gilts in individual stalls and diet change increase lipid and protein catabolism, oxidative stress, activate the innate immune system and cause a certain degree of tissue damage. both techniques showed alterations in two main homeostatic mechanisms: the innate immune and redox systems. blood and saliva samples were collected a day before and up to four days after changing the housing system. other proteins related to the innate immune system, including lactotransferrin, protegrin and galectin were also identified by itraq, as well as oxidative stress enzymes such as peroxiredoxin and glutathione peroxidase proteomics also revealed the decrease of apolipoproteins, and the presence of intracellular proteins in serum, which may indicate physical injury to tissues. we propose that valuable assays for stress assessment in gilts may be based on a score composed by a combination of salivary cortisol, lipid metabolites, innate immunity and oxidative stress markers and intracellular proteins. acute phase proteinsbiomarkerindividual confinementinnate immunitylipidsoxidative stresspigproteomicsstresshttp://dx.doi.org/ /501100003030agència de gestió d’ajuts universitaris i de recerca2014sgr-4bassols anna http://dx.doi.org/ /501100003329ministerio de economía y competitividadagl2011-30598-c03-02bassols anna issue-copyright-statement© the author 2016 in this work, we have analyzed the changes in serum composition in gilts after switching their housing, from pen to individual stalls, which is generally accepted to cause animal discomfort.
 the objective assessment of animal stress and welfare requires proper laboratory biomarkers. in this work, we have analyzed the changes in serum composition in gilts after switching their housing, from pen to individual stalls, which is generally accepted to cause animal discomfort.  blood and saliva samples were collected a day before and up to four days after changing the housing system. biochemical analyses showed adaptive changes in lipid and protein metabolism after the housing switch, whereas cortisol and muscular markers showed a large variability between animals. 2d-dige and itraq proteomic approaches revealed variations in serum protein composition after changing housing and diet of gilts. both techniques showed alterations in two main homeostatic mechanisms: the innate immune and redox systems. the acute phase proteins haptoglobin, apolipoprotein a-i and α1-antichymotrypsin and the antioxidant enzyme peroxiredoxin were found differentially expressed by 2d-dige. other proteins related to the innate immune system, including lactotransferrin, protegrin and galectin were also identified by itraq, as well as oxidative stress enzymes such as peroxiredoxin and glutathione peroxidase proteomics also revealed the decrease of apolipoproteins, and the presence of intracellular proteins in serum, which may indicate physical injury to tissues.  housing of gilts in individual stalls and diet change increase lipid and protein catabolism, oxidative stress, activate the innate immune system and cause a certain degree of tissue damage. we propose that valuable assays for stress assessment in gilts may be based on a score composed by a combination of salivary cortisol, lipid metabolites, innate immunity and oxidative stress markers and intracellular proteins. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords acute phase proteinsbiomarkerindividual confinementinnate immunitylipidsoxidative stresspigproteomicsstresshttp://dx.doi.org/ /501100003030agència de gestió d’ajuts universitaris i de recerca2014sgr-4bassols anna http://dx.doi.org/ /501100003329ministerio de economía y competitividadagl2011-30598-c03-02bassols anna issue-copyright-statement© the author 2016 