Skip to main content

Pathogen exposure influences immune parameters around weaning in pigs reared in commercial farms



Multiple antigenic stimulations are crucial to immune system training during early post-natal life. These stimulations can be either due to commensals, which accounts for the acquisition and maintenance of tolerance, or to pathogens, which triggers immunity. In pig, only few works previously explored the influence of natural exposition to pathogens upon immune competence. We propose herein the results of a multicentric, field study, conducted on 265 piglets exposed to contrasted pathogen levels in their living environment. Piglets were housed in 15 different commercial farms, sorted in two groups, low (HSLOW)- and high (HSHIGH)-health status farms, depending on their recurrent exposition to five common swine pathogens.


Using animal-based measures, we compared the immune competence and growth performances of HSLOW and HSHIGH pigs around weaning. As expected, we observed a rise in the number of circulating leucocytes with age, which affected different cell populations. Monocyte, antigen-experienced and cytotoxic lymphocyte subpopulation counts were higher in piglets reared in HSLOW farms as compared to their HSHIGH homologs. Also, the age-dependent evolution in γδ T cell and neutrophil counts was significantly affected by the health status. With age, circulating IFNα level decreased and IgM level increased while being greater in HSLOW piglets at any time. After weaning, LPS-stimulated blood cells derived from HSLOW piglets were more prone to secrete IL-8 than those derived from HSHIGH pigs did. Monocytes and granulocytes issued from HSLOW pigs also exhibited comparable phagocytosis capacity. Altogether our data emphasize the more robust immunophenotype of HSLOW piglets. Finally, piglets raised under higher pathogen pressure grew less than HSHIGH piglets did and exhibited a different metabolic profile. The higher cost of the immune responses associated with the low farm health status may account for lower HSLOW piglet performances.


Altogether, our data, obtained in field conditions, provide evidence that early exposure to pathogens shapes the immune competence of piglets. They also document the negative impact of an overstimulation of the immune system on piglets’ growth.

Peer Review reports


If the development of the immune system is initiated during the foetal life, it is still largely incomplete at birth. This accounts for the increased susceptibility of neonates to clinically relevant infections [1]. During the neonatal period, the microbial colonisation of mucosal tissues plays a key role in the shaping of immunity, especially facilitating tolerance to commensal micro-organisms [2]. Post-natal maturation of the immune system also requires training through multiple non-commensal antigenic stimulations. Thus, the overall immune competence of an individual depends on the diversity of antigens he had been confronted to, particularly early in life. Early pathogen exposure obviously influences the immune landscape of an individual, through multiple mechanisms. For example, it can imprint the human T cell repertoire with potential consequences for pathogen defences [3].

Classical animal models are poorly relevant to elucidate the training effect of pathogen exposure on individuals. It is obviously the case for mice housed in standardised, specific pathogen-free (SPF) conditions, which are immunologically inexperienced. To overcome this limitation, Beura et al. worked on mice trapped from free-living feral populations [4]. They demonstrated that these “dirty” mice exhibited an experienced immune phenotype, close to the human one, in contrast with the one of age-matched SPF mice. Pig exhibits strong anatomical, physiological and immunological similarities with human, while being housed in diverse environments [5]. As such, it appears as a valuable model to address this question.

In pig, to our knowledge, only few works purposefully explored how natural exposition to pathogens triggers heterologous immunity. These studies were essentially conducted on adult pigs, reared in experimental facilities. For example, Clapperton et al. analysed the effects of the farm health status, defined by the presence or not of specific pathogens, on some immune parameters. In pigs older than 3 months, they reported no major effect on blood monocyte and lymphocyte sub-populations counts [6]. Other studies compared the effects of hygiene conditions on immune parameters, using “dirty” vs “clean” environments, to mirror different levels of antigenic stimulation. For example, van der Meer et al. evidenced a poor effect of hygiene conditions on white blood cell counts during the fattening period [7]. In other studies, bad hygiene conditions were rather shown to induce greater circulating neutrophil numbers in pigs [8, 9].

In commercial farms, pigs are naturally exposed to variable pathogen pressures. As a consequence, field studies are particularly relevant to document the link between pathogen exposure and immune competence. Along pig’s life, weaning is the most critical step. Piglets are abruptly separated from their mother and mixed with non-littermates, while being exposed to a novel biotic environment. This stressful transition, usually characterised by a decreased food intake and a subsequent growth slowdown, can also triggers long-term consequences for pigs’ health [10].

In this context, we report herein the results of a field study conducted on 265 pigs naturally exposed to contrasted pathogen levels. For this purpose, we selected 15 commercial farms sorted in two groups, low (HSLOW) and high (HSHIGH) health status farms, depending on their recurrent exposition to five common swine pathogens. Using animal-based measures, we compared the immune competence and growth performances of HSLOW and HSHIGH pigs around weaning.


We followed the evolution of 17 immune and 4 metabolic parameters in 142 HSLOW and 123 HSHIGH piglets, sampled and weighed two days before and seven days after weaning. Individual raw data are available in Additional file 6 while a summary is presented in Table 1.

Table 1 Raw data of immune markers, metabolic parameters and performances in HSLOW and HSHIGH piglets before (t0) and after (t1) weaning

White blood cell counts

Neutrophil counts were not different between HSLOW and HSHIGH piglets both before (p = 0.82) and after (p = 0.14) weaning (Fig. 1a). They increased with age (p < 0.0001) but to a larger extent in piglets reared under high pathogen pressure (HS × Age, p < 0.05). Monocyte (Fig. 1b) and lymphocyte (Fig. 1c) numbers were affected by both the farm health status (p < 0.05) and the age (p < 0.0001). At any time, both were significantly greater in piglets reared in the worst environmental conditions. The age-dependent increase in monocyte and lymphocyte counts was not affected by pathogen exposure (HS × Age, p = 0.26 and 0.30 respectively). Eosinophils were undetectable or present in very low numbers (detection threshold 104 cells/ml) in piglets before weaning under any condition. Their count dramatically increased after weaning, reaching comparable levels in both groups (Additional file 1). Of note, basophil counts were undetectable-to-low under any condition (data not shown).

Fig. 1
figure 1

White blood cell counts in HSLOW versus HSHIGH piglets, before (grey) and after (black) weaning. Neutrophil (A), monocyte (B) and lymphocyte (C) numbers were obtained from the hematology analyser. Emmeans with lower and upper confidence intervals’ limits are shown on graphs. Insets, results of the linear mixed-effects model analyses

Lymphocyte subpopulations analyses

We next assessed whether the different lymphocyte subpopulations behaved as the whole lymphocyte population did. The gating strategies used for T and B cells’ analyses are described in Additional files 2 and 3 respectively. For Ag-experienced Th (A) and cytotoxic T (B) lymphocytes, cell counts were affected by both the farm health status (p < 0.01) and the age (p < 0.0001), without interaction (Additional file 4). Cell numbers increased with age in any of the environmental conditions and were higher in HSLOW piglets at any time. After weaning, the proportion of antigen-experienced cells among CD4+ T lymphocytes was 25.8% in piglets reared under low farm health status versus 21.7% in HSHIGH pigs.

As shown in Fig. 2, three lymphocyte subpopulations, CD8α γδ T (A), naive Th (B) and B (C) lymphocytes were affected by the age (p < 0.0001) but not by the farm health status (p = 0.87, 0.11 and 0.31 respectively). Interestingly, the age-dependent increase in the CD8α γδ T lymphocyte count was greater in HSHIGH piglets which were reared under low pathogen pressure (HS × Age, p < 0.0001).

Fig. 2
figure 2

γδ T, naive CD4+ T and B lymphocyte counts in HSLOW versus HSHIGH piglets before (grey) and after (black) weaning. yδ T lymphocytes (A CD3+ CD4 CD8−/med), naive CD4+ Th lymphocytes (B CD3+ CD4+) and B lymphocytes (C CD3 CD21+) were discriminated using flow cytometry and their numbers calculated. Emmeans with lower and upper confidence intervals’ limits are shown on graphs. Insets, results of the linear mixed-effects model analyses

Plasma levels of IgM and IFNα

Both the farm health status (p < 0.05) and the age (p < 0.0001) affected circulating IgM concentration, without interaction (HS × Age, p = 0.60) (Fig. 3a). Indeed, piglets reared under HSLOW conditions displayed greater IgM levels, both before and after weaning. With age, we observed a significant rise in IgM concentration (p < 0.0001). We also evidenced a decrease in IFNα concentration with age (Fig. 3b; p < 0.01), independently of pathogen pressure.

Fig. 3
figure 3

Plasmatic IgM and IFNα in HSLOW versus HSHIGH piglets before (grey) and after (black) weaning. IgM (A) and IFNα (B) levels were measured in blood samples. Emmeans with lower and upper confidence intervals’ limits are shown on graphs. Insets, results of the linear mixed-effects model analyses

Assessment of immune competence after weaning

To further evaluate the ability of piglets to cope with pathogen challenge, functional tests were performed seven days after weaning and raw data are shown on Fig. 4.

Fig. 4
figure 4

Cytokine secretion and phagocytic capacity of leucocytes issued from HSLOW versus HSHIGH piglets after weaning. A Secretion of IL-10, TNFα and IL-8 was measured in response to the ex vivo exposure of blood cells to LPS in WBA. B Fluorescent bacteria phagocytosis by blood cells was analysed by flow cytometry in HSLOW versus HSHIGH piglets. Medians are shown on graphs

We first measured the ability of blood cells to secrete IL-10, TNFα and IL-8 in response to an ex vivo stimulation by LPS (Fig. 4a). IL-10 secretion by blood cells was not significantly modulated by the health status of piglets (p = 0.5). Since TNFα and IL-8 were not detectable in 49 and 94 upon 264 samples respectively, we analysed these parameters as qualitative variables. TNFα detection was not modulated by the health status (p = 0.35). In contrast, blood cells derived from HSLOW piglets were more likely to secrete detectable levels of IL-8 than HSHIGH piglets did (Odd ratio = 0.41 [0.23–0.67], p < 0.001). Finally, the phagocytosis ability of blood neutrophils and monocytes was not affected by the farm health status (Fig. 4b).

Plasma metabolite concentrations after weaning

We measured circulating glucose, free fatty acid, urea and creatinine levels in fasted piglets one week after weaning and synthesised them using a Principal Component Analysis (PCA) (Additional file 5). The first three components explained more than 80% of the total inertia and each contributed to more than 20% of the total variance. The first dimension was related to “alternative energy supply” while the second one represented “major energy stores”, finally, the third dimension appeared more integrative with equivalent contributions of glucose, free-fatty acids and urea. Dimensions 1 and 2 were not modulated by the farm health status. In contrast, the third dimension was significantly affected by the health status (p < 0.001) with the metabolic profile of HSLOW piglets (emmeans = − 0.37 [− 0.70 to − 0.03]) shifted towards urea instead of glucose and free-fatty acids such as the one of HSHIGH piglets (emmeans = 0.43 [0.07–0.79]).

Individual growth performances and clinical signs

Finally, piglets were weighed at the two time points (see Table 1 for raw data). Before weaning, the weight of HSLOW and HSHIGH piglets was not different with estimated marginal means of 7.82 kg [7.31–8.32] for HSLOW piglets, and 8.20 kg [7.66–8.75] for HSHIGH piglets (p = 0.22). In contrast, along the period we studied, HSLOW piglets grew significantly less than HSHIGH piglets did, with estimated marginal means of 158 [126–190] and 219 [185–253] g/day respectively (p < 0.05).

Along the period we studied, we detected no major health problems in any of the groups. After weaning, at day 35, only fifteen individuals upon 265 exhibited diarrhoea and the diarrhoeic piglets were equally distributed between HSLOW and HSHIGH farms (Additional file 6).


In this field study, we demonstrated the impact of pathogen exposure on the immune system maturation and growth performances of piglets around weaning. For this purpose, we took advantage of the natural exposure of piglets to common swine pathogens in commercial farms. Using individual blood indicators, we compared the immune competence and the growth of male piglets two days before and one week after weaning.

As already described, we observed a rise in the number of circulating leucocytes with age, which affected different cell populations. At weaning, the immune system is still developing, circulating cell numbers reaching stable levels after the age of 10 weeks [11]. The neutrophil count increased less with age than the lymphocyte count did, which resulted in a decrease in the neutrophil-to-lymphocyte ratio (data not shown). According to Stull et al., this decrease persists until the acquisition of the immune competence [12].

As we expected, the farm health status significantly affected several immune parameters. Plasma IgM levels were greater in piglets grown under HSLOW conditions compared to HSHIGH piglets. After birth, piglets acquire maternal immunoglobulins, mainly IgG from colostrum, before producing their own, starting with IgM [13]. Thus, the greater IgM levels appear as the consequence of a stronger antigenic stimulation due to the higher microbial pressure in HSLOW farms.

To characterise the effects of the farm health status on lymphocytes, we focused on subpopulations that were robustly quantified by flow cytometry on blood samples. Indeed, while it could have been interesting to quantify regulatory T cells and Th17 cells in our samples, intracellular staining is poorly adapted to whole-blood samples. We observed significantly greater numbers of antigen-experienced and cytotoxic T lymphocytes in piglets raised in HSLOW farms. In swine, CD3+ γδ T lymphocytes can represent up to 50% of total circulating lymphocytes [14]. This lymphocyte subtype is crucial to young animals, before the complete maturation of their immune system [15]. In particular, the CD8α fraction of γδ T cells was shown to dominate within the peripheral blood lymphocyte compartment in pigs [16]. In this context, the evolution of CD8α γδ T cell counts with age in piglets raised under low or high pathogen pressure is particularly interesting. Indeed, the increase in CD8α γδ T cell counts around weaning was lower in HSLOW piglets raised under high pathogen pressure. In contrast, in these piglets, the age-dependent rise in the neutrophil count was more pronounced. Altogether, these characteristics may be considered as hallmarks of a more mature immunophenotype in HSLOW piglets.

The concept of trained immunity was initially described in plants and invertebrates. It refers to the ability of innate immune cells to secrete higher amounts of pro-inflammatory cytokines when exposed to a second, unrelated, damage-associated molecular pattern [17]. This innate immune memory was demonstrated to be crucial for juveniles that still lack efficient adaptative immunity [18]. After weaning, while we did not evidence any difference in LPS-induced IL-10 and TNFα secretion, we observed that leucocytes derived from HSLOW piglets were more prone to secrete the pro-inflammatory cytokine IL-8. Although we cannot exclude a bias related to the greater monocyte count evidenced in HSLOW pigs, it might reveal different trained immunity patterns. Indeed, the stronger capacity of their leucocytes to respond to LPS could confer to HSLOW pigs a better ability to face Gram-negative bacterial infections as previously suggested for humans [19]. HSLOW pigs also exhibited greater numbers of monocytes with comparable phagocytosis capacity, emphasising their more robust immunophenotype.

We also confirmed herein the extra cost of the immune training through antigenic stimulation in HSLOW pigs. It is well known that individuals have to balance the fitness benefits of a robust immune response against the costs of drawing resources [20]. In agreement with this concept, piglets born and reared in HSHIGH farms grew better than their homologs reared under low health status. Since we did not measure the feed consumption of piglets, we do not know whether HSLOW piglets ate less than HSHIGH piglets did. Nevertheless, it has already been documented that poor sanitary conditions after weaning can deteriorate growth performances despite a limited negative effect on feed intake [21]. We presume that the decreased performances of HSLOW piglets were likely due to the nutritional cost of adaptation to the higher pathogen pressure as previously described [22] and suggested by the different metabolic profiles of the HSLOW and HSHIGH piglets one week after weaning.

One of the limitations of this study is that it was performed on males only, with the aim to reduce heterogeneity due to the well-known influences of the gender on immune responses [23]. Early in life, in humans, circulating immune cell counts are generally greater in males than females. In young pigs also, some immune parameters were shown to differ between males and females [24]. In this context, part of the differences we observed might not have been seen in females. The season might constitute another bias of this study. Indeed, piglets were either born in winter or early spring, which may have influenced the circulation of pathogens. However, samplings in HSLOW farms were interspersed by experiments in HSHIGH farms. Finally, in this multicentric, on-field study, although we could not control the whole environmental conditions, we selected the farms as relatively homogenous in terms of genetics and batch management production system. In order to reduce the bias related to other parameters such as feed or ambient conditions (air quality including relative humidity or dust and ammonia concentration…), we took into account the farm effect in the statistical models we used.

In this study, the two defined health statuses did not represent highly contrasted situations since they were representative of field conditions found in commercial farms. As a consequence, at any time point, the large majority of piglets was healthy. After weaning, only few piglets expressed digestive symptoms, which is in accordance with the farmers’ health recordings made between day 26 and day 48 in each post-weaning room. Indeed, as previously reported, the total number of days with recorded diarrhoea was 26 in the HSLOW farms compared to 22 in the HSHIGH farms [25]. Nevertheless, the impact of pathogen exposure on immune parameters was not insignificant which confirms that early exposure to pathogens shapes the immune competence of piglets. Recent findings reported that housing conditions in early life influenced specific antibodies and leucocyte populations in older piglets, suggesting quite long-term effects [26]. Whether the differences we observed in some immune parameters persist or alleviate all along pigs’ life will have to be explored.


Our data obtained in field conditions describe the immune and metabolic consequences of variable levels of pathogen exposure on piglets, echoing previous findings connecting pathogen exposure to immune system shaping. It also dovetails with previous literature confirming the nutritional cost associated with the overstimulation of the immune system.


Ethical statement

All experimental procedures were carried out in accordance with the European directive 2010/63/EU for animal experiments. It was submitted to and approved by the Oniris’ ethics committee for veterinary clinical research (decision CERVO-2016-6-V).

Selection of the farms

This study was conducted in 15 commercial farms in the West of France during a 5 month-period, from January to May. Farms were selected a priori based on the analysis of the annual farm health records made by the attending veterinarians. Recurrent exposure to five major porcine pathogens, namely Influenza virus H1N1, Porcine Reproductive and Respiratory Syndrome virus, Streptococcus suis, Actinobacillus pleuropneumoniae and Lawsonia intracellularis, was used to determine farms’ health status (HS). Farms were further classified as HSLOW, when more than 2 pathogens were detected, or conversely HSHIGH. Each farm health status was confirmed at the end of the study using the record of the year. The overall characteristics of the farms (batch management system, animals’ genetics, nutrition, post-weaning facilities…) are provided in Additional file 7. Of note, farms were located in the periphery of the city of Nantes to allow rapid processing of the samples (transport time < 2 h).

Animals, experimental design and sample collection

We studied healthy entire male piglets [Nucleus germ line: Pietrain x (Landrace x Large White)]. Indeed, in the selected farms, castration was no longer practiced. In each farm, 18 individuals, issued from 9 sows (2 piglets/sow) of different parities (representative of the demographic pyramid of the farm), were initially included. We selected middleweight piglets upon visual inspection of the litters. Unless specified in the Table 1, analyses were performed on 142 HSLOW versus 123 HSHIGH piglets (respectively from 8 and 7 farms), due to the exclusion of five deceased piglets equally distributed between both groups.

Blood sampling was performed in the morning, two days before and one week after weaning (that occurred at day 28), when piglets were in average 26 (t0) and 35 (t1) day-old. For the second sampling, piglets had been fasted overnight to allow biochemical analyses (i.e. glycemia, creatininemia, uremia and free-fatty acidemia). Blood was collected from the jugular vein of pigs slightly maintained in a supine position, by experienced staff, in less than 30 s per pig. The total volume of blood collected never exceeded 25 ml per pig. Immediately after blood sampling, piglets were individually weighed using a portable electronic scale. Careful examination of the animal was performed to detect any clinical sign. At t1, rectal temperature was measured and the occurrence of diarrhoea was estimated upon faeces consistency (liquid or not).

Complete blood cell counts

Complete blood cell counts were determined from EDTA-treated blood samples using the clinical-grade Procyte Dx haematology analyser (IDEXX, Saint-Denis, France). White blood cell parameters included total white blood cell (WBC), neutrophil, eosinophil, lymphocyte and monocyte counts.

Analyses of lymphocyte subpopulations

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


Ex vivo phagocytosis was assessed using the Phagotest™ Kit (BD Biosciences). Briefly, 50 µl of heparinised blood were incubated for 10 min at 37 °C with opsonised FITC-labelled E. coli bacteria. For each analysis, ice-incubated negative controls were included. Samples were immediately analysed by flow cytometry (20 000 cells/sample analysed) using FSC/SSC dot plot to discriminate granulocytes (FSChigh) from mononuclear cells (FSClow). Results are expressed as percent of phagocytes among mono- and poly-nuclear cells.

Biochemical analyses

Glycemia, uremia and creatininemia were analysed on heparinised plasma samples collected from fasted piglets, using an automated analyser (Konelab, Thermo Scientific, Courtaboeuf, France) and commercial kits (Thermo Scientific). Free fatty acids levels were determined on EDTA-treated plasma samples with a commercial kit (Wako-Chemicals, Neuss, Deutschland). IgM concentrations were determined by a homemade quantitative sandwich ELISA on heparinised plasma samples. Briefly, serially diluted plasmas were incubated for 1 h in 96-well plates coated with an excess of goat antibody against swine IgM Fc. After washing, a peroxidase-conjugated goat anti-swine IgM Fc was added. Capture and conjugated antibodies were purchased from MyBiosource (Clinisciences, Nanterre, France).

Plasma and whole blood assay (WBA) cytokine measurements

Porcine IFNα was analysed in heparinised plasma samples using a homemade sandwich ELISA test. The test used K9 anti-IFNα and biotinylated-F17 anti-IFNα antibodies, and recombinant swine IFNα as standard (PBL Assay Science, Piscataway, USA). The intra-assay coefficient of variation (CV) was 3% while inter-assay CV was 18%.

For WBA, heparinised blood samples were diluted five-fold in RPMI 1640 medium and stimulated for 18 h in duplicate using 10 ng/ml Escherichia coli O111:B4 LPS (Sigma-Aldrich, St-Quentin-Fallavier, France). Porcine TNFα, IL-8 and IL-10 were then quantified by ELISA with respective detection limits of 15.6, 62.5 and 23 pg/ml (Bio-Techne, Noyal-Chatillon sur Seiche, France).

Statistical analyses

Statistical analyses were performed using R software (version 4.1.2, R Core Team, 2016) while charts were graphed using Graph Pad Prism (version 9.3.1).

To describe the consequences of the farm health status and/or the age upon immune and endocrine systems in piglets, linear mixed-effects models were performed using the lme4 package [27]. Whenever required, data were initially transformed. We used logarithm (neutrophil, γδ T lymphocyte, Ag-experienced Th cell and cytotoxic lymphocyte counts, blood IgM, IFNα, urea and free fatty acid concentrations) or square root (monocyte, total lymphocyte, naive CD4+ T and B lymphocyte counts). The goodness-of-fit of all models (normality and homoscedasticity of residuals and distribution of random effects) was checked by graphic procedures and by Shapiro–Wilk tests.

The main model included, as fixed effects, the farm health status (HS: low and high), the age (Age: t0 and t1), the sow parity (1 versus 2 and above) and the interaction between HS and Age, and, as random effects, the farm (1 to 15) and the individual (Id), as described below.

$$\begin{aligned} & Y_{ijt} = \beta_{0} + Parity_{ijt} \beta_{1} + Score_{ijt} \beta_{2} + Age_{ijt} \beta_{3} + Score_{ijt} Age_{ijt} \beta_{4} + \eta_{ij} + \nu_{j} + \varepsilon_{ijt} \\ & \eta_{ij} \sim Normal\left( {0,\sigma_{\eta }^{2} } \right) \\ & \nu_{j} \sim Normal \left( {0,\sigma_{\nu }^{2} } \right) \\ & \varepsilon_{ijt} \sim Normal\left( {0,\sigma_{\varepsilon }^{2} } \right) \\ \end{aligned}$$

where \(Y_{ijt}\) is the quantitative variable of interest of the animal i in the farm j at age t, \(Parity_{ijt}\) is the parity in 2 classes (1 vs 2 and +), \(Score_{ijt}\) is the Health Status (LOW vs HIGH), \(Age_{ijt}\) is the sampling date (t0 vs t1), \(Score_{ijt} Age_{ijt}\) is the interaction between age and score, \(\eta_{ij}\) is the random effect for animal i in the farm j. Animal random effect follows a normal distribution with mean 0 and variance \(\sigma_{\eta }^{2}\). \(\nu_{j}\) is the random effect for farm j. Farm random effect follows a normal distribution with mean 0 and variance \(\sigma_{\nu }^{2}\). \(\varepsilon_{ijt}\) is the residual for the animal i in the farm j at age t which follows a normal distribution with mean 0 and variance \(\sigma_{\varepsilon }^{2}\).

We presented only the effects we were interested in, i.e. the ones of the HS, the age and their interaction, and all results were adjusted on sow parity. Of note, 7 parameters were also influenced by sow parity, namely neutrophil, lymphocyte and monocyte counts, experienced Th cell numbers, cytotoxic T cell counts and IFNα level. Also, as expected, piglets issued from parity 1 sows exhibited a lower initial weight (7.47 kg [7.03–7.91]) than the other ones (8.55 kg [8.17–8.93]).

When the interaction effect was significant, a Tukey’s Honest Significant Difference post-hoc test was performed to search for relevant differences and results are shown beside the graphs.

For quantitative variables measured only at t1, only the HS and the sow parity were included in the model as fixed effects, and the farm as a random one. For this model, metabolic parameters (glycemia, free-fatty acidaemia, uraemia and creatininemia) were firstly processed using a Principal Component Analysis (PCA) [28] with the FactoMineR package in the R software [29] in order to reduce the number of variables while minimizing information loss.

TNFα and IL-8 variables were analysed as qualitative variables (< or > 15.6 for TNFα and (< or > 62.5 for IL-8) using mixed logistic regression models included the HS, the farm and the sow parity as fixed effects and the farm as a random effect.

To assess the effect of the HS on growth (ADG), the mixed linear model included the health status, the sow parity and the t0 weight as fixed effects, and the farm as a random one.

For most response variables, estimated marginal means (emmeans) and lower and upper limits of the 95% confidence interval [IC95] were calculated using the emmeans package [30] and graphed on charts.

Availability of data and materials

The datasets used and/or analysed during the current study are provided in the supplementary file section.



Average Daily Gain


Estimated marginal mean


Free-Fatty Acids


Health status


Confidence Interval


Principal Component Analysis


Specific Pathogen Free


Whole Blood Assay


White Blood Cell


  1. Ygberg S, Nilsson A. The developing immune system—from foetus to toddler. Acta Paediatr Oslo Nor. 1992;2012(101):120–7.

    Article  CAS  Google Scholar 

  2. Gensollen T, Iyer SS, Kasper DL, Blumberg RS. How colonization by microbiota in early life shapes the immune system. Science. 2016;352:539–44.

    Article  CAS  Google Scholar 

  3. Gaimann MU, Nguyen M, Desponds J, Mayer A. Early life imprints the hierarchy of T cell clone sizes. Elife. 2020;9:e61639.

    Article  CAS  Google Scholar 

  4. Beura LK, Hamilton SE, Bi K, Schenkel JM, Odumade OA, Casey KA, et al. Normalizing the environment recapitulates adult human immune traits in laboratory mice. Nature. 2016;532:512–6.

    Article  CAS  Google Scholar 

  5. Pabst R. The pig as a model for immunology research. Cell Tissue Res. 2020;380:287–304.

    Article  Google Scholar 

  6. Clapperton M, Diack AB, Matika O, Glass EJ, Gladney CD, Mellencamp MA, et al. Traits associated with innate and adaptive immunity in pigs: heritability and associations with performance under different health status conditions. Genet Sel Evol GSE. 2009;41:54.

    Article  CAS  Google Scholar 

  7. van der Meer Y, Lammers A, Jansman AJM, Rijnen MMJA, Hendriks WH, Gerrits WJJ. Performance of pigs kept under different sanitary conditions affected by protein intake and amino acid supplementation. J Anim Sci. 2016;94:4704–19.

    Article  Google Scholar 

  8. Le Floc’h N, Knudsen C, Gidenne T, Montagne L, Merlot E, Zemb O. Impact of feed restriction on health, digestion and faecal microbiota of growing pigs housed in good or poor hygiene conditions. Anim Int J Anim Biosci. 2014;8:1632–42.

    Article  Google Scholar 

  9. Chatelet A, Gondret F, Merlot E, Gilbert H, Friggens NC, Le Floc’h N. Impact of hygiene of housing conditions on performance and health of two pig genetic lines divergent for residual feed intake. Anim Int J Anim Biosci. 2018;12:350–8.

    Article  CAS  Google Scholar 

  10. Campbell JM, Crenshaw JD, Polo J. The biological stress of early weaned piglets. J Anim Sci Biotechnol. 2013;4:19.

    Article  Google Scholar 

  11. Juul-Madsen HR, Jensen KH, Nielsen J, Damgaard BM. Ontogeny and characterization of blood leukocyte subsets and serum proteins in piglets before and after weaning. Vet Immunol Immunopathol. 2010;133:95–108.

    Article  CAS  Google Scholar 

  12. Stull CL, Kachulis CJ, Farley JL, Koenig GJ. The effect of age and teat order on alpha1-acid glycoprotein, neutrophil-to-lymphocyte ratio, cortisol, and average daily gain in commercial growing pigs. J Anim Sci. 1999;77:70–4.

    Article  CAS  Google Scholar 

  13. Butler JE, Wertz N, Sinkora M. Antibody repertoire development in swine. Annu Rev Anim Biosci. 2017;5:255–79.

    Article  CAS  Google Scholar 

  14. Talker SC, Käser T, Reutner K, Sedlak C, Mair KH, Koinig H, et al. Phenotypic maturation of porcine NK- and T-cell subsets. Dev Comp Immunol. 2013;40:51–68.

    Article  CAS  Google Scholar 

  15. Le Page L, Baldwin CL, Telfer JC. γδ T cells in artiodactyls: Focus on swine. Dev Comp Immunol. 2022;128:104334.

    Article  CAS  Google Scholar 

  16. Gerner W, Käser T, Saalmüller A. Porcine T lymphocytes and NK cells—an update. Dev Comp Immunol. 2009;33:310–20.

    Article  CAS  Google Scholar 

  17. Netea MG, Joosten LAB, Latz E, Mills KHG, Natoli G, Stunnenberg HG, et al. Trained immunity: a program of innate immune memory in health and disease. Science. 2016;352:aaf1098.

    Article  CAS  Google Scholar 

  18. Levy O, Wynn JL. A prime time for trained immunity: innate immune memory in newborns and infants. Neonatology. 2014;105:136–41.

    Article  CAS  Google Scholar 

  19. Wurfel MM, Park WY, Radella F, Ruzinski J, Sandstrom A, Strout J, et al. Identification of high and low responders to lipopolysaccharide in normal subjects: an unbiased approach to identify modulators of innate immunity. J Immunol Baltim Md. 1950;2005(175):2570–8.

    Article  Google Scholar 

  20. Lochmiller RL, Deerenberg C. Trade-offs in evolutionary immunology: just what is the cost of immunity? Oikos. 2000;88:87–98.

    Article  Google Scholar 

  21. Pastorelli H, LeFloch N, Merlot E, Meunier-Salaün MC, van Milgen J, Montagne L. Sanitary housing conditions modify the performance and behavioural response of weaned pigs to feed- and housing-related stressors. Anim Int J Anim Biosci. 2012;6:1811–20.

    Article  CAS  Google Scholar 

  22. Le Floch N, Jondreville C, Matte JJ, Seve B. Importance of sanitary environment for growth performance and plasma nutrient homeostasis during the post-weaning period in piglets. Arch Anim Nutr. 2006;60:23–34.

    Article  CAS  Google Scholar 

  23. Klein SL, Flanagan KL. Sex differences in immune responses. Nat Rev Immunol. 2016;16:626–38.

    Article  CAS  Google Scholar 

  24. Christoforidou Z, MoraOrtiz M, Poveda C, Abbas M, Walton G, Bailey M, et al. Sexual dimorphism in immune development and in response to nutritional intervention in neonatal piglets. Front Immunol. 2019;10:5.

    Article  CAS  Google Scholar 

  25. Luise D, Sciellour ML, Buchet A, Resmond R, Clement C, Rossignol M-N, et al. The fecal microbiota of piglets during weaning transition and its association with piglet growth across various farm environments. PLoS ONE. 2021;16:e0250655.

    Article  CAS  Google Scholar 

  26. Luo L, Jansen CA, Bolhuis JE, Arts JAJ, Kemp B, Parmentier HK. Early and later life environmental enrichment affect specific antibody responses and blood leukocyte subpopulations in pigs. Physiol Behav. 2020;217:112799.

    Article  CAS  Google Scholar 

  27. Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.

    Article  Google Scholar 

  28. Jolliffe IT, Cadima J. Principal component analysis: a review and recent development. Philos Trans R Soc A. 2016;374:20150202.

    Article  Google Scholar 

  29. Lê S, Josse J, Husson F. FactoMineR: an R package for multivariate analysis. J Stat Softw. 2008;25:1–18.

    Article  Google Scholar 

  30. Lenth RV. emmeans: estimated marginal means, aka least-squares means. 2021.

Download references


The authors thank pig farmers that were involved in this field study. The authors are grateful to JY Audiard, M Lefebvre and G Martin for on-farm technical assistance as well as to F. Thomas for laboratory assistance.


This study was funded by the INRAE Integrated Management of Animal Health metaprogram as part of the “GISA-SEVROBUST” project and Cooperl Arc Atlantique. EB and SG received an INRAE/Pays de la Loire Region PhD funding.

Author information

Authors and Affiliations



JH, AB, ET, PM, EM, CB and BL conceived and designed experiments; JH, KH, AB, EB, GM, MA, MLM, EM and BL performed experiments; JH, KH, SG, EM, AL, JMB and BL analysed results; JH and BL prepared the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Julie Hervé.

Ethics declarations

Ethics approval and consent to participate

All experimental procedures were carried out in accordance with the European directive 2010/63/EU for animal experiments. They were submitted to and approved by the Oniris’ ethics committee for veterinary clinical research (decision CERVO-2016-6-V).

Consent for publication

Not applicable.

Competing interests

PhD funding of AB was obtained from Cooperl Innovation, but Cooperl Innovation did not have any additional role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Other authors declare no competing interest.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1

. Eosinophil cell counts in HSLOW and HSHIGH piglets after weaning.

Additional file 2

. Gating strategy for the identification of T lymphocyte sub-populations.

Additional file 3

. Gating strategy for the identification of B lymphocytes.

Additional file 4

. Ag-experienced Th (A) and cytotoxic T lymphocyte (B) counts inHSLOW versus HSHIGH piglets before (grey) and after (black) weaning.

Additional file 5

. Principal Component Analysis on metabolic data.

Additional file 6

. Individual raw data.

Additional file 7

. Characterisation of the farms.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hervé, J., Haurogné, K., Buchet, A. et al. Pathogen exposure influences immune parameters around weaning in pigs reared in commercial farms. BMC Immunol 23, 61 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: