# Early Salmonella Typhimurium infection in pigs disrupts Microbiome composition and functionality principally at the ileum mucosa

Héctor Argüello, Jordi Estellé, Sara Zaldívar-López, Ana Carvajal, Mª

Asunción López-Bascón, Fiona Crispie, Orla O'sullivan, Paul D. Cotter, Feliciano Priego-Capote, Luis Morera, et al.

► To cite this version:

Héctor Argüello, Jordi Estellé, Sara Zaldívar-López, Ana Carvajal, Mª Asunción López-Bascón, et al. Early Salmonella Typhimurium infection in pigs disrupts Microbiome composition and functionality principally at the ileum mucosa. Scientific Reports, 2018, 8 (1), 10.1038/s41598-018-26083-3. hal-02621111

HAL Id: hal-02621111

https://hal.inrae.fr/hal-02621111

Submitted on 26 May 2020

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

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

OPEN

# Early *Salmonella* Typhimurium infection in pigs disrupts Microbiome composition and functionality principally at the ileum mucosa

Received: 5 October 2017

Accepted: 2 May 2018

Published online: 17 May 2018

Héctor Argüello<fnref n="1" />, Jordi Estellé², Sara Zaldívar-López<fnref n="1" />, Ángeles Jiménez-Marín<fnref n="1" />, Ana Carvajal³, Mª Asunción López-Bascón⁴, Fiona Crispie⁵,⁶, Orla O'Sullivan⁵,⁶, Paul D. Cotter⁵,⁶, Feliciano Priego-Capote⁴, Luis Morera<fnref n="1" /> & Juan J. Garrido<fnref n="1" />

*Salmonella* is a major foodborne pathogen which successfully infects animal species for human consumption such as swine. The pathogen has a battery of virulence factors which it uses to colonise and persist within the host. The host microbiota may play a role in resistance to, and may also be indirectly responsible from some of the consequences of, *Salmonella* infection. To investigate this, we used 16S rRNA metagenomic sequencing to determine the changes in the gut microbiota of pigs in response to infection by *Salmonella* Typhimurium at three locations: ileum mucosa, ileum content and faeces. Early infection (2 days post-infection) impacted on the microbiome diversity at the mucosa, reflected in a decrease in representatives of the generally regarded as desirable genera (i.e., *Bifidobacterium* and *Lactobacillus*). Severe damage in the epithelium of the ileum mucosa correlated with an increase in synergistic (with respect to *Salmonella* infection; *Akkermansia*) or opportunistically pathogenic bacteria (*Citrobacter*) and a depletion in anaerobic bacteria (*Clostridium spp.*, *Ruminococcus*, *or Dialliser*). Predictive functional analysis, together with metabolomic analysis revealed changes in glucose and lipid metabolism in infected pigs. The observed changes in commensal healthy microbiota, including the growth of synergistic or potentially pathogenic bacteria and depletion of beneficial or competing bacteria, could contribute to the pathogen's ability to colonize the gut successfully. The findings from this study could be used to form the basis for further research aimed at creating intervention strategies to mitigate the effects of *Salmonella* infection.

*Salmonella* is a major foodborne pathogen with a significant number of cases and outbreaks every year<fnref n="1" />. Poultry meat, eggs and pork are the main sources of human salmonellosis. As a consequence of the successful establishment of control programmes in poultry, the relative proportion of human cases attributed to pork consumption has increased².

Swine are a natural host for *Salmonella* spp. Apart from infections caused by S. Cholerasuis, the host adapted serovar to swine, *Salmonella* infections in pigs usually vary from subclinical to a self-limiting diarrhoea³. Among serovars affecting swine, *Salmonella* Typhimurium stands out due to its zoonotic importance in human salmonellosis<fnref n="1" /> and its widespread distribution in pork production⁴,⁵. The relevance of S. Typhimurium infections in swine is reflected by the vast literature published with respect to the pathogenesis⁶, the host immune response⁷,⁸,

<footnote n="1">¹Grupo de Genómica y Mejora Animal, Departamento de Genética, Facultad de Veterinaria, Universidad de Córdoba, 14047, Córdoba, Spain. ²GABI, INRA, AgroParisTech, Université Paris-Saclay, 78350, Jouy-en-Josas, France. ³Departamento de Sanidad Animal, Facultad de Veterinaria, Universidad de León, 24071, León, Spain. ⁴Departamento de Química Analítica Universidad de Córdoba, Córdoba, CeiA3 Campus de Excelencia Agroalimentaria, Universidad de Córdoba, 14047, Córdoba, Spain. ⁵Teagasc, Food Research Centre, Moorepark, Fermoy, Co., Cork, Ireland. ⁶APC Microbiome Institute, Cork, Ireland. Correspondence and requests for materials should be addressed to H.A. (email: ge2arroh@uco.es)</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Estimation of evenness and richness (Shannon and Chao1) in faeces, ileum content and ileum mucosa samples during the time course of S. Typhimurium infection. Significant differences were observed regardless the estimator used when the type of sample factor was analysed (*p* < 0.001).</figcaption></figure>

epidemiology of infections³ and control⁹,¹⁰ of this serotype in pigs. Notably, it has been reported that within a few hours of ingestion, the pathogen is detected in high concentrations in faeces¹¹. Early after infection, epithelial damage and inflammation can be observed, principally at the ileum mucosa¹², together with a migration of the pathogen to the gut associated lymphoid tissue via monocyte-derived cells¹³.

High throughput sequencing methods have provided an opportunity to study the porcine microbiota with unprecedented detail through sequencing of the 16S rRNA gene as a phylogenetic marker¹⁴⁻¹⁶. Indeed, these high throughput sequencing methods offer the chance to understand how infections such as swine salmonellosis impact the gut microbiota. For instance, this approach has been employed to identify differences in the microbiota among high and low *Salmonella* shedders¹⁷. In their study, Bearson and colleagues revealed a higher abundance of representatives of the family *Ruminococcaceae* in low-shedder pigs and also depletion of *Prevotella* in high shedder pigs early after infection. Another recent study, where pigs were challenged with S. Typhimurium, showed that strain virulence and the resulting intensity of the inflammatory response following infection by *Salmonella* caused a reduction in the abundance of short-chain fatty acids-producing bacteria¹⁸.

Understanding the changes occurring in the host microbiome due to *Salmonella* infection is of paramount importance with respect to developing a better understanding of *Salmonella* pathogenesis and disease. Here we characterise for the first time how the gut microbiome is modified after a *Salmonella* Typhimurium challenge in recently weaned pigs within three target locations of the gastrointestinal tract, i.e., faeces, ileum content and ileum mucosa. By studying changes in the diversity, taxonomy and abundance, and by combining this information with the predicted functionality of the microbiome and metabolite analysis, we aimed to expand on previous studies¹⁷⁻¹⁹, and provide a greater insight into the mechanisms underlying successful colonization of the gut by *Salmonella*¹⁷⁻¹⁹. Ultimately, it is hoped that identification of the groups of bacteria affected by disruption of gut homeostasis following *Salmonella* infection, could serve as a preliminary step to design strategies to limit the alteration of the microbiota by *Salmonella* infection and to increase resistance to gut colonization by this pathogen.

## Results

58014 operational taxonomic units (OTUs) were detected among the samples following matching of the 16S rRNA sequencing data to the Greengenes database. After filtering low abundance OTUs using a 0.005% threshold, a final number of 1192 taxa across seven taxonomic ranks were retained for further analysis.

**Factors affecting microbiome richness.** Sample richness was estimated by Chao1 and Shannon indices while Simpson index was used to estimate sample evenness (Fig. 1; Supplementary file 1). By repeated measure of analysis of variance (ANOVA), it was apparent that there were no significant differences in richness or evenness between infected and control pigs (*p* = 0.977; *p* = 0.75) or across sampling time points, (*p* = 0.96; *p* = 0.67). However, when ileum mucosa samples were specifically analysed, differences in richness among sampling time points were close to significance (Shannon index; *p* = 0.061). Across sample type, significant differences were observed in the richness estimators between faeces, ileum content and ileum mucosa (*p* < 0.001). Faecal samples showed higher richness than ileum content or ileum mucosa regardless of the estimator used (Fig. 1).

**Ordination of the pig gut microbiota varies according to the location in the gut from which samples are collected and between Day 0 and Day 2 samples.** Two factors, i.e., type of sample (*p* < 0.001) and individual pig (*p* < 0.01), influenced the ordination of samples (Fig. 2a). The effect was maintained even after splitting samples by sampling day (Fig. 2b). Considering this output, the next ordination analyses were split by type of sample and the analysis of faeces, ileum mucosa and ileum content samples separately revealed that both infection (*p* < 0.01) and sampling day (*p* < 0.05) had a significant influence (Fig. 2). When samples from Day 0 (control group) and 2 days post-infection (d.p.i.) were specifically compared, significant
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Ordination analysis of the samples obtained through the time-course of a S. Typhimurium infection in pigs, using Bray-Curtis dissimilarity index and plot by Non-metric Multidimensional Scaling (NMDS). (a) Ordination of samples by sample type factor (ileum content, ileum mucosa and faeces). (b) Splits the NMDS ordination by sampling day. (c,d) Include the ordination of samples from Day 0 (control pigs) and Day 2 (infected pigs). In (c) dots are coloured by infection factor, non-infected control pigs (red) and infected pigs at 2 days post-infection (dpi) (blue). (d) Represents the same ordination colouring the type of sample factor. See that ileum mucosa samples from 2 dpi pigs are positioned far from the rest of the samples green circle.</figcaption></figure>

differences in Bray-Curtis non-linear multi-dimensional scaling (NMDS) ordination were detected at within the ileum mucosa microbiota (*p* < 0.001) and close to significant differences (*p* = 0.053) were apparent in faecal samples (Fig. 2c and d).

**Taxonomic differences in faeces, ileum content and mucosa microbiomes at the time points evaluated.** We compared the relative abundances of bacteria present in faeces, ileum content and ileum mucosa at different time points, including late infection sampling (14 d.p.i. and 30 d.p.i.; Fig. 3). OTUs were ranked by relative abundance, selecting up to the 10th most abundant taxa (families and genera) in each sample type (faeces, ileum mucosa and ileum content) at each sampling point (Day 0, Day 2, Day 14 and Day 30). Overall, we observed large differences among samples across a variety of variables. On Day 0, *Lactobacillaceae* (*Lactobacillus*) dominated the bacteria present within the ileum contents. At the ileum mucosa *Lactobacillus*, *Helicobacteraceae* (*Flexispira*) and *Ruminococcaceae* (*Ruminococcus*) were the most abundant bacteria. In faeces, proportions of *Lactobacillaceae* were lower and the dominant family was *Ruminococcaceae*. Clear changes in dominant families and genera were observed in infected pigs at Days 2, 14 and 30 across intestinal locations. These differences were evident within ileum content and ileum mucosa samples, and are described in greater depth below. In contrast, the dominant families *Prevotellaceae* (*Prevotella*), *Lachnospiraceae* (*Rosseburia* and *Coprococcus*) and *Ruminococcaceae* (*Ruminococcus*) remained relatively stable in faeces at different sampling points.

**Changes in the abundance of OTUs as a consequence of the S. Typhimurium infection.** The analysis of changes in microbiota abundance at 2 d.p.i. revealed 60 differentially abundant OTUs (FDR-value < 0.05) in the ileum mucosa (Supplementary file 2), when compared with the corresponding samples from the control group (0 d.p.i.). Similarly, 25 and 10 OTUs were found to be differentially abundant (FDR-value < 0.05) across the ileum content (Supplementary file 3) and faecal samples (Supplementary file 4), respectively. Particular OTUs differed significantly in abundance across the gut regions analysed, as depicted in Fig. 4a. Table 1 summarises the principal increased and decreased OTUs in the ileum mucosa, ileum content and in faeces after two days of infection. Notably, a decrease in abundance of OTUs belonging to the class *Clostridia* was apparent across all three locations analysed. Specifically, genera such as *Clostridium* (FDR = 0.037) and *Dialister* (FDR = 0.002) in the ileum mucosa, *Megasphaera* (FDR = 0.017) and *Mitsuokella* (FDR = 0.035) in the ileum content and two
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Relative percentage abundance of the main taxa (Phyla, Familiae and Genera) detected on each sampling day in each type of sample analysed (faeces, ileum content and ileum mucosa).</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Differences in OTUs between control pigs (Day 0) and infected pigs (2 days post-infection) at the ileum mucosa, ileum content and faeces. (a) Venn diagram showing the overlapping of differentially expressed OTUs in ileum mucosa, ileum content and faeces between Day 0 and Day 2 p.i., pigs. (b) Heatmap illustrating the mean relative abundance of multiple OTUs at the ileum mucosa of non-infected control pigs (green) and Day 2 p.i. infected pigs (red). White colour indicates low abundance while dark blue high values of abundance. Dendrogram was built using hierarchical cluster analysis with Bray-Curtis dissimilarity indices. OTUs in the area included in the red frame belong to the family *Veillonellaceae* (*Mitsuokella*, *Megasphaera*, and *Dialister* genera).</figcaption></figure>

OTUs from families *Lachnospiraeae* (FDR = 0.033) and *Ruminococcaceae* (FDR = 0.027) in faeces, were all present in decreased proportions after two days of infection. In addition, OTUs from the genus *Lactobacillus* were also decreased in the ileum mucosa and content samples (FDR = 0.007, FDR = 0.042; respectively) at this time point.
<table><caption>**Table 1.** Principal OTUs changed (FDR < 0.05) in faeces, ileum content and ileum mucosa samples between control pigs (Day 0) and *Salmonella* infected pigs (Day 2 post-infection). <sup>a</sup>Normalised counts of OTUs. <sup>b</sup>Adjusted p-value by false discovery rate.</caption><thead><tr><th></th><th>Phylum</th><th>Class</th><th>Family</th><th>Genus</th><th>Counts Control group<sup>a</sup></th><th>Counts Infected group<sup>a</sup></th><th>Scaling factor</th><th>Adjusted p-value<sup>b</sup></th></tr></thead><tbody><tr><td rowspan="10">FAECES</td><td>Firmicutes</td><td>Clostridiales</td><td>Lachnospiraceae</td><td>Unclassified</td><td>1055</td><td>153</td><td>5.59</td><td>0.027</td></tr><tr><td>Firmicutes</td><td>Clostridiales</td><td>Ruminococcaceae</td><td>Unclassified</td><td>8601</td><td>821</td><td>4.26</td><td>0.033</td></tr><tr><td>Firmicutes</td><td>Erysipelotrichales</td><td>Erysipelotrichaceae</td><td>p-75-a5</td><td>3874</td><td>141</td><td>2.28</td><td>0.018</td></tr><tr><td>Chlamydiae</td><td>Chlamydiales</td><td>Chlamydiaceae</td><td>Chlamydia</td><td>1</td><td>1580</td><td>2.23</td><td>0.008</td></tr><tr><td>Firmicutes</td><td>Clostridiales</td><td>Ruminococcaceae</td><td>Unclassified</td><td>3273</td><td>22</td><td>1.69</td><td>0.008</td></tr><tr><td>Firmicutes</td><td>Clostridiales</td><td>Ruminococcaceae</td><td>Unclassified</td><td>22</td><td>2529</td><td>0.39</td><td>0.008</td></tr><tr><td>Firmicutes</td><td>Clostridiales</td><td>Mogibacteriaceae</td><td>Unclassified</td><td>0</td><td>344</td><td>0.25</td><td>0.025</td></tr><tr><td>Proteobacteria</td><td>Tremblayales</td><td>Unclassified</td><td>Unclassified</td><td>9</td><td>570</td><td>-0.52</td><td>0.040</td></tr><tr><td>Firmicutes</td><td>Lactobacillales</td><td>Lactobacillaceae</td><td>Lactobacillus mucosae</td><td>474</td><td>24</td><td>9.46</td><td>0.042</td></tr><tr><td>Proteobacteria</td><td>Campylobacterales</td><td>Campylobacteraceae</td><td>Campylobacter</td><td>267</td><td>2</td><td>9.14</td><td>0.036</td></tr><tr><td rowspan="6">ILEUM CONTENT</td><td>Firmicutes</td><td>Lactobacillales</td><td>Streptococcaceae</td><td>Streptococcus</td><td>707</td><td>3</td><td>8.53</td><td>0.020</td></tr><tr><td>Firmicutes</td><td>Clostridiales</td><td>Veillonellaceae</td><td>Megasphaera</td><td>312</td><td>0</td><td>7.35</td><td>0.017</td></tr><tr><td>Firmicutes</td><td>Clostridiales</td><td>Veillonellaceae</td><td>Mitsuokella</td><td>1210</td><td>66</td><td>5.96</td><td>0.035</td></tr><tr><td>Proteobacteria</td><td>Enterobacteriales</td><td>Enterobacteriaceae</td><td>Unclassified</td><td>0</td><td>690</td><td>-2.93</td><td>0.027</td></tr><tr><td>Firmicutes</td><td>Clostridiales</td><td>Lachnospiraceae</td><td>Clostridium</td><td>688</td><td>74</td><td>3.93</td><td>0.041</td></tr><tr><td>Firmicutes</td><td>Clostridiales</td><td>Lachnospiraceae</td><td>Clostridium</td><td>802</td><td>83</td><td>3.71</td><td>0.037</td></tr><tr><td rowspan="15">ILEUM MUCOSA</td><td>Firmicutes</td><td>Clostridiales</td><td>Veillonellaceae</td><td>Dialister</td><td>764</td><td>0</td><td>3.35</td><td>0.002</td></tr><tr><td>Actinobacteria</td><td>Bifidobacteriales</td><td>Bifidobacteriaceae</td><td>Bifidobacterium</td><td>6456</td><td>4</td><td>3.05</td><td>1.081E-06</td></tr><tr><td>Actinobacteria</td><td>Bifidobacteriales</td><td>Bifidobacteriaceae</td><td>Bifidobacterium</td><td>3132</td><td>1</td><td>3.04</td><td>1.022E-05</td></tr><tr><td>Actinobacteria</td><td>Bifidobacteriales</td><td>Bifidobacteriaceae</td><td>Bifidobacterium</td><td>1767</td><td>3</td><td>3.02</td><td>1.565E-04</td></tr><tr><td>Actinobacteria</td><td>Bifidobacteriales</td><td>Bifidobacteriaceae</td><td>Bifidobacterium</td><td>669</td><td>2</td><td>2.86</td><td>0.002</td></tr><tr><td>Actinobacteria</td><td>Bifidobacteriales</td><td>Bifidobacteriaceae</td><td>Bifidobacterium</td><td>45161</td><td>38</td><td>2.67</td><td>1.695E-09</td></tr><tr><td>Firmicutes</td><td>Lactobacillales</td><td>Lactobacillaceae</td><td>Lactobacillus</td><td>445</td><td>0</td><td>2.59</td><td>0.007</td></tr><tr><td>Verrucomicrobia</td><td>Verrucomicrobiales</td><td>Verrucomicrobiaceae</td><td>Akkermansia</td><td>103</td><td>1691</td><td>-2.39</td><td>0.016</td></tr><tr><td>Firmicutes</td><td>Clostridiales</td><td>Ruminococcaceae</td><td>Oscillospira</td><td>13</td><td>300</td><td>-2.84</td><td>0.032</td></tr><tr><td>Bacteroidetes</td><td>Bacteroidales</td><td>Odoribacteriaceae</td><td>Odoribacter</td><td>30</td><td>367</td><td>-4.51</td><td>0.035</td></tr><tr><td>Actinobacteria</td><td>Actinomycetales</td><td>Corynebacteriaceae</td><td>Corynebacterium</td><td>4</td><td>537</td><td>-5.00</td><td>0.002</td></tr><tr><td>Firmicutes</td><td>Lactobacillales</td><td>Streptococcaceae</td><td>Streptococcus</td><td>154</td><td>22415</td><td>-6.43</td><td>0.011</td></tr><tr><td>Firmicutes</td><td>Clostridiales</td><td>Clostridiaceae</td><td>Candidatus Arthromitus</td><td>1</td><td>62633</td><td>-11.07</td><td>2.991E-08</td></tr></tbody></table>

The proportions of a number of taxa also increased at this time point (i.e., 2 dpi pigs). Genera such as *Candidatus arthromitus* (FDR < 0.001), *Akkermansia* (FDR = 0.016) and *Odoribacter* (FDR = 0.035) were increased within the mucosa of infected pigs, one OTU belonging to the *Enterobacteriaceae* family (genera unclassified) (FDR = 0.027) was significantly increased in ileum content samples and OTUs that were significantly increased in faeces were assigned to the families *Mogibacteriaceae* (FDR = 0.025) and *Ruminococcaceae* (FDR = 0.08).

Hierarchical clustering of pig samples according to the relative count of differentially abundant OTUs in the ileum mucosa, clustered pigs on the basis of infection (Fig. 4b). Samples collected from the four necropsied pigs at 2 d.p.i. (infected pigs) exhibited a similar pattern of OTUs abundance that diverged from the control pigs necropsied on Day 0. This clustering of pigs by the relative abundance of OTUs was not observed either in the ileum content or faeces for Day 0 and Day 2 pigs (Supplementary Figure S1a and b respectively).

**Correlation between histopathology scoring and mucosa-associated microbiota.** Histological lesions in ileum mucosa samples from pigs necropsied on Day 0 and Day 2 were scored and correlated with OTU proportions at the mucosa (Fig. 5; Supplementary file 5). Severe epithelial damage was positively correlated to taxa from the phyla *Proteobacteria*, such as *Citrobacter* (Spearman's ρ = 0.97 and *p* < 0.001), *Verrucomicrobia* such as *Akkermansia muciniphila* (Spearman's ρ = 0.93 and *p* < 0.001), *Firmicutes*, such as *Allobaculum* (Spearman's ρ = 0.92 and *p* < 0.001) or *Roseburia* (Spearman's ρ = 0.88 and *p* < 0.01), and *Bacteroidetes*, such as *Bacteroides* (Spearman's ρ = 0.87 and *p* < 0.01).

In contrast, the integrity of the intestinal epithelial cell wall was positively associated with different *Firmicutes*, such as *Lactobacillus* (Spearman's ρ = -0.90 and *p* < 0.01), *Clostridium* (Spearman's ρ = -0.88 and *p* < 0.01), *Bulleidia* (Spearman's ρ = -0.87 and *p* < 0.01) and *Ruminococcus* (Spearman's ρ = -0.87 and *p* < 0.01), *Proteobacteria*, such as *Campylobacter* (Spearman's ρ = -0.88 and *p* < 0.01) and *Flexispira* (Spearman's ρ = -0.87 and *p* < 0.01), and the genus *Bifidobacterium* (phylum *Actinobacteria*) (Spearman's ρ = -0.88 and *p* < 0.01).
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Genus correlated to integrity and damage of the epithelial of the mucosa. (A) Non-damaged mucosa from a control pig (score = 0), (B) damaged mucosa with villi shortened and *Salmonella* infiltration (score = 3) and (C) damaged mucosa with complete erosion of villi and high concentration of *Salmonella* (score = 5). (D) Relative counts of taxa correlated to epithelium damage at the ileum mucosa (*p* < 0.01).</figcaption></figure>

**Analysis of modulation of predicted microbiome functionality by *Salmonella* infection.** The PICRUSt tool was used to predict Kegg Orthology group (KOs) abundances based on 16S rRNA sequences. The analysis yielded a prediction of 5441 KOs. As for OTU-based analysis, Bray-Curtis compositional dissimilarity highlighted significant differences in predicted functions (*p* < 0.001) by type of sample (faeces, ileum content or ileum mucosa (Supplementary Figure S2)). Further analyses focused on differences between ileum mucosa samples from Day 0 and Day 2 pigs. No significant differences in metagenome functions were observed at KO level 2 functions although a tendency (*p* < 0.1) was observed in functions such as membrane transport or biosynthesis of secondary metabolites, all with lower abundance in infected pigs (Supplementary Figure S3). Analyses of functions at L3 in these samples identified 1138 differentially abundant KOs among control pigs (Day 0) and infected pigs (2 d.p.i) (Supplementary file 6). The inclusion of differentially expressed KOs in the software Interactive pathway explorer (Ipath2) showed an increased abundance of genes linked with the metabolism of lipids (biosynthesis of fatty acids) and nucleotide metabolism, and a decreased abundance of genes associated with bacterial chemotaxis, glycolysis and gluconeogenesis pathways in the ileum mucosa of Day 2 infected pigs.

Ileum content samples from control pigs (Day 0) and infected pigs (Day 2) were also analysed by liquid and gas chromatography. The analysis targeted fatty acids, fatty acyls and short-chain fatty acids compounds, carboxylic acid compounds and sugar derivatives. Table 2 summarises the 19 compounds that were significantly changed between control and infected pigs. Among fatty acids, there was a remarkable increase in concentration of dihydroxybutanoic (*p* = 0.026) acid, octanoic acid (*p* = 0.026) and acetoacetic acid (*p* = 0.002) in 2 dpi infected pigs. Among carboxylic acids, the levels of pyruvic (*p* = 0.002) and aconitic (*p* = 0.015) acids were significantly decreased in the ileum contents of infected pigs. Finally, glycerol, a sugar derivate, was also significantly decreased (*p* = 0.002) in the ileum content of Day 2 infected pigs.

Discussion

This study involved an in-depth investigation of the microbial changes that occur in the gut environment at target locations, ileum content and mucosa as well as faeces, after a S. Typhimurium infection in swine. S. Typhimurium, either under controlled or field conditions, is able to successfully infect pigs in a short time-frame, with relatively high concentrations of the pathogen being present in faeces within a few hours of exposure¹¹,²⁰. A number of studies agree that the local inflammation, resulting from the immune response elicited by *Salmonella*, triggers changes in the gut environment which favour the survival of *Salmonella* and prompt changes in the gut microbiota
<table><caption>**Table 2.** Chromatography table (mean normalised value) of differentiated metabolites between control pigs (Day 0) and *Salmonella* infected pigs (Day 2 post-infection) in ileum content samples.</caption><thead><tr><th>Compound</th><th>Formula</th><th>Control group</th><th>Infected group (2 dpi)</th><th>p-value</th></tr></thead><tbody><tr><td colspan="5"><strong>Fatty acids, fatty acids and short chain fatty acids derivatives</strong></td></tr><tr><td>Methylvaleric acid</td><td>C6H12O2</td><td>6765.50</td><td>10608.50</td><td>0.026</td></tr><tr><td>Oxovaleric acid</td><td>C5H8O3</td><td>1043565.50</td><td>6997420.00</td><td>0.0021</td></tr><tr><td>Acetoacetic acid</td><td>C4H6O3</td><td>56857.00</td><td>438702.50</td><td>0.0022</td></tr><tr><td>Dihydroxyisovaleric acid</td><td>C5HO4</td><td>412918.00</td><td>216793.50</td><td>0.0022</td></tr><tr><td>Dihydroxybutanoic acid</td><td>C4H8O2</td><td>182645.00</td><td>1595772.50</td><td>0.026</td></tr><tr><td>Octanoic acid</td><td>C8H16O2</td><td>196746.50</td><td>563696.00</td><td>0.026</td></tr><tr><td>Dihydroxy-docosanoic acid</td><td>C22H44O4</td><td>30227.16</td><td>5058.92</td><td>0.0021</td></tr><tr><td>Dihydroxystearic acid</td><td>C18H36O4</td><td>3918.66</td><td>314.47</td><td>0.0021</td></tr><tr><td>HODE</td><td>C18H32O3</td><td>19330.85</td><td>1830.28</td><td>0.0021</td></tr><tr><td>Hydroxyhexadecatetraenoic acid</td><td>C16H24O2</td><td>5224.16</td><td>5210.28</td><td>0.0086</td></tr><tr><td>Hydroxyoctadecenoic acid</td><td>C18H36O3</td><td>13607.09</td><td>17868.10</td><td>0.0259</td></tr><tr><td>Hydroxy tetracosanoic acid</td><td>C24H48O3</td><td>22198.33</td><td>4904.09</td><td>0.0259</td></tr><tr><td>Oxohexadecanoic acid</td><td>C16H30O3</td><td>280040.30</td><td>183366.84</td><td>0.0021</td></tr><tr><td colspan="5"><strong>Carboxilic acids</strong></td></tr><tr><td>Aconitic acid</td><td>C6H6O6</td><td>54789.00</td><td>20940.50</td><td>0.0152</td></tr><tr><td>Pyruvic acid</td><td>C3H4O3</td><td>7913.00</td><td>2186.00</td><td>0.0022</td></tr><tr><td>Butenedioic acid</td><td>C4H4O4</td><td>95914.50</td><td>21361.00</td><td>0.0022</td></tr><tr><td>Pantothenic acid</td><td>C9H17NO5</td><td>1741.02</td><td>7536.10</td><td>0.0086</td></tr><tr><td colspan="5"><strong>Sugar derivatives</strong></td></tr><tr><td>Glycerol</td><td>C3H8O3</td><td>10788.00</td><td>3411.00</td><td>0.0022</td></tr><tr><td>Mannitol</td><td>C6H14O6</td><td>71854.00</td><td>505248.50</td><td>0.0022</td></tr></tbody></table>

composition<sup>18,21,22</sup>. This inflammatory response is thought to enable *Salmonella* to compete with the resident microbiota through the use of tetrathionate produced by the oxidation of thiosulfate by reactive oxygen species (ROS)<sup>23</sup> or due to its resistance to defensive molecules such as ROS or lipocalin-2, which is secreted by neutrophils and limits the acquisition of iron<sup>24</sup>. As a result, changes in the composition of the microbiota and gut functionality were anticipated. Our study included a limited number of pigs per sampling day. Although, ideally, a larger sample size would have been employed, our investigation with the limited number of pigs available, together with the appropriate statistical analysis, allow us to report, for the first time to our knowledge, that *Salmonella*’s biggest impact on microbiota composition takes places at the ileum mucosa as early as 2 d.p.i.

Methods of evaluating microbial diversity include α-diversity indexes such as richness (the total number of unique species within a community) and the relative abundance of different species that make up the richness in that community, and as well as β-diversity, dissimilarity-based ordinations. These metrics were clearly influenced by the origin of the sample (ileum content, ileum mucosa or faeces) in accordance with previous results<sup>25,26</sup>. When samples were split by sample type, the influence of the infection on the ordination of the microbiota at the ileum mucosa was particularly apparent. In contrast, it was surprising to note that although challenged pigs exhibited clinical signs, such as mild to watery diarrhoea and low-grade fever, as well as the steady shedding of *Salmonella* during the first week (data not shown), no clear significant differences in diversity were observed in faeces (only a trend), or ileum content samples. Although the sample size limits the power to detect differences in diversity in faeces, our observation is consistent with the results obtained by Bearson and colleagues<sup>17</sup>, who observed that neither alpha nor beta diversity were affected when challenged pigs shed *Salmonella* at low concentrations in faeces. Furthermore, Drumo *et al.*<sup>18</sup> showed that *Salmonella* strain virulence influences the observed changes of diversity results in faeces. In our results, we establish for the first time that the major diversity differences occur at the ileum mucosa and that this effect is less apparent within the ileum content or the downstream segments of the intestine (i.e., faeces).

In agreement with previous studies in pigs<sup>17-19</sup>, we found that *Salmonella* infection triggered changes in the relative abundance of certain bacteria at 2 d.p.i. The shift in the microbiota was apparent at the ileum mucosa (see dendogram in Fig. 4b), where samples from 2 d.p.i. infected pigs clustered together in the same clade, influenced by the drop in abundance of genera from the family *Veillonellaceae* (*Mitsuokella*, *Dialister* and *Megasphaera*). Since the ileum constitutes the principal target for *Salmonella* invasion in the gut<sup>27</sup>, it can be inferred from our results that the aforementioned impacts on epithelial cell damage, inflammation, defensive molecules and metabolites principally affect the microbiota associated to the ileum mucosa, while the effect on ileum content or faeces is more limited.

To evaluate the effect of the infection on the relative abundance of the microbiota, relative changes in taxa abundance between Day 0 and Day 2 were evaluated. In agreement with previous studies, we observed dissimilarities in OTU changes at different locations after *Salmonella* infection<sup>19</sup>. Although the microbiota is not stable and evolves through the pigs’ life<sup>26</sup>, the short time-frame between challenge and infection (2 days) ensures that the changes observed cannot be attributed to age-related influences. Notably, our results demonstrated a decrease
in abundance of beneficial gut bacteria such as *Clostridium* spp., *Lactobacillus* spp., *Prevotella*, *Bifidobacterium* or *Megaesphaera* in the ileum (mucosa and content) after infection. The decrease in *Lactobacillus* was one of the most robust results with respect to the microbiota of the ileum content and ileum mucosa, with a number of different OTUs of these species ranking among the most decreased. S. Typhimurium infection in pigs prompts an early inflammatory response which is associated with a decreased expression of genes coding for apical sodium-dependent bile acid transporter (ASBT) and the fatty and bile acid transporters FABP2 and FABP6, impairing normal bile absorption in the ileum⁸. As a result, higher concentrations of bile salts are present in the lumen of the small intestine. As some beneficial species such as *Lactobacillus* and *Bifidobacterium*, are susceptible to bile salts²⁸,²⁹, the marked drop of abundance of these bacteria in the ileum could be the consequence of the disruption of the entero-hepatic circulation of the biliary acids. However, it remains to be seen if these changes are transitory, i.e., as a consequence of the challenge with a high dose of *Salmonella*, or if they persist for longer periods. It is notable that previous work does show that the changes in microbiome community structure persist over the first week post-challenge¹⁷,¹⁹.

We also observed taxa changes specific to each location, e.g., *Dialister* only changed at the mucosa, *Veillonella* at the ileum content and *Treponema* in faeces. Even if divergences in changes in the microbiota at different locations could be explained by the differences of the background microbiota²⁵, we propose that a diluted effect in faeces compared to the ileum may also be a contributory factor. Through the analysis of the most abundant taxa in samples Day 14 and Day 30, it was apparent that the faecal microbiome was more stable throughout the course of the study than was the case for the ileal microbiomes¹⁴. This result highlights that the microbiome is more susceptible to changes in targeted locations (e.g. ileum mucosa) than in faeces (which can be considered a global picture of the gut microbiome), a fact that is particularly notable in light of the fact that faeces is regularly used in microbiome studies due to the fact that sample procurement is non-invasive³⁰. Although our results may be limited by the study sample size, we illustrate that the interpretation of the results should be validated in the context of the anatomical location of the studied effect when faeces are used as proxy samples.

Damage of ileal mucosa was correlated to an increase of bacteria that are synergistic to the infectious process, i.e., bacteria that boost the damage caused by the pathogen, and/or opportunistic pathogens, such as *Akkermansia muciniphila* and *Citrobacter*, respectively. Previous studies have reported an increase in levels of other *Enterobacteriaceae*, such as *Citrobacter*, arising from *Salmonella* infection²²,³¹. Our study confirms that the modification of the intestinal environment at the ileum favours the expansion of closely related microorganisms²²,³²,³³. *Akkermansia* is an opportunistic bacterium which develops synergistic relationships with *Salmonella* by disturbing mucus secretion and siltation¹⁹. Similarly to our study, Borton *et al*.²² detected an increase in the proportions of the genus *Akkermansia* genera in *Salmonella* infection with inflamed murine gut samples. In contrast, typically regarded as beneficial, anaerobic, Gram-positive gut bacteria such as *Lactobacillus*, *Ruminococcus*, *Clostridium* or *Allobaculum* abundance decreased in severe damaged ileum mucosa samples. The release of RNS and ROS during intestinal non-typhoidal *Salmonella* infection creates a highly oxidative environment which is not permissive for the growth of anaerobic bacteria³⁴, while transient aerotolerant bacteria such as *Bacteroides*³⁵ might take an advantage of the disruption in the gut environmental conditions.

Finally, based on PICRUSt predictions and metabolomic analysis, we inferred changes in functional pathways encoded within the ileum mucosa microbiome associated with infection. Ordination of predicted functions by Bray-Curtis dissimilarity clustered samples by tissue, similarly to 16S rRNA OTU results. At KEGG level 2, there were no clear differences in microbiome functions among control and early infected pigs at the ileum mucosa. The sample size together with the overlapping of functions among different bacteria are two possibly limiting factors in finding potential differences³⁰. Nevertheless, some functions related to basal metabolism (secondary metabolism, metabolism of enzymes, amino acids and biosynthesis of lipids) demonstrated a tendency to be reduced at 2 d.p.i., an observation that suggests a general downshift in microbial intestinal activity. These results are in accordance with a study of induced colitis in pigs³⁶, which described a reduced capacity for energy acquisition and dysregulated microbial signaling and repair pathways after inflammation. Further analysis of functions at KEGG level 3 interestingly revealed that chemotaxis and glucose metabolism (glycolysis and gluconeogenesis) functions were more abundant at the ileum mucosa prior to *Salmonella* infection (control pigs). The activation of the glucose metabolism could be associated with the degradation of oligo- and polysacharides by lactic acid bacteria. Levels of glycerol, a sugar derivative, and pyruvate, a carboxylic acid, were decreased in infected pigs, according to the predictions in the metabolism of glucose and its derivatives. In contrast, levels of pantothenic acid were increased in the ileum contents of infected pigs at 2 d.p.i. Analysis of predicted metabolic routes in Interactive Pathways Explorer (Ipath2) showed an activation of valine, leucine and isoleucine metabolism, all of which result in the production of pantothenic acid. At 2 d.p.i., an increase of fatty acid metabolism KOs was observed, with particular activation of fatty acid biosynthesis routes. Levels of methylvaleric, oxovaleric, acetoacetic, butenedioic and octanoic acids were increased in ileum contents from 2 d.p.i. pigs, confirming the activation of biosynthesis of fatty acids by the gut microbiota early after infection. The link between the microbiota and fatty acid metabolism could be attributed to changes in the nutrients available in the gut environment. *Salmonella* infection disrupts the entero-hepatic circulation of the biliary acids as mentioned above⁸,³⁰, which in turn alters the lipid metabolism. As a result, we can expect higher concentrations of available fatty acids in the gut lumen, favouring the growth of bacteria that use them in their metabolism. Overall, these results illustrate that S. Typhimurium infection triggers impairment of metabolic functions, which in turn can affect the physiology of digestive processes in the gut³⁰.

## Conclusions

The present study demonstrates that *Salmonella* infection modifies the microbiota in the gut early after infection, and shows for the first time that the major changes occur in the ileum mucosa and are less evident in the faeces and even the ileum lumen, highlighting the merits of study host mucosal samples, where practical. Ironically, host response to infection (immune response and metabolic changes) could be a major contributor to the
depletion of commensal/beneficial inhabitants of the intestinal tract (*Lactobacillus*, *Bifidobacterium*, *Prevotella* or *Megasphaera*), and the increase of synergists of infection such as *Akkermansia* or *Citrobacter*. Indeed, the relative abundance of these synergists in the ileal microbial ecosystem is positively correlated to the degree of epitelium damage in the ileum. Finally, we illustrate that the disruption of microbiota composition triggered by infection also has apparent consequences for the microbiome function, with a general downshift of activities early after infection and additional changes in relevant metabolic routes. Future research should focus on identifying the key host response-related factors that prompt changes in the microbiome structure, which might help to establish novel microbiome restoration strategies to limit the impact of *Salmonella* infection in pigs and, ultimately, zoonotic infections in humans.

## Methods

**Salmonella challenge outline.** Sixteen crossbred commercial weaned female pigs of approximately four weeks of age from a *Salmonella*-free herd were moved into the isolation facilities at the University of León (Spain). After an adaptation period of seven days, four pigs were euthanized and necropsied (control group). The remaining twelve pigs were orally challenged with 5 ml of a broth containing $10^8$ cfu/ml of S. Typhimurium phage type DT104 and subsequently necropsied in groups of four at 2 days post-infection (d.p.i.), 14 d.p.i., and 30 d.p.i.

Faeces and ileum (approx., 20 cm backwards from ileocaecal junction) were collected from each animal. Ileum luminal content (approx. 5 g) was poured in a plastic tube and the ileum wall was preserved in another tube. All samples were collected aseptically and frozen in liquid nitrogen immediately. Frozen samples were kept at -80 °C until further processing. A subsample of faeces was cultured by the current ISO:6579/2007 to detect the presence/absence of *Salmonella* in all pigs included in the study. Another sub-sample of ileum wall was fixed in formalin and embedded in paraffin wax for further processing.

All procedures involving animals were performed in accordance with the European regulations regarding the protection of animals used for experimental and other scientific purposes, under the supervision of the Ethical and Animal Welfare Committee of the University of León (Spain).

**Histopathology and immunohistochemistry.** Two paraffin sections (5 µm) from each formalin fixed samples were routinely processed and stained with hematoxylin and eosin (H&E) to evaluate tissue morphology. For immunohistochemistry assays, a standard avidin-biotin peroxidase method was performed as described elsewhere³⁷ employing 1F12 monoclonal antibody and a biotinylated anti-mouse IgG (Dako, Barcelona, Spain) as a secondary antibody. Sections were evaluated by microscopy and scored (from 0 to 5) by adding up the following parameters: no lesions = 0; villi shorten = 1; complete villi erosion = 1; presence of *Salmonella* = 1; high concentration of *Salmonella* = 1; inflammatory cells infiltration = 1.

**DNA extraction.** Frozen ileum wall (5–10 cm of ileum) was opened using a scalpel, mucosa wall was rinsed carefully with sterile Phosphate Saline Buffer to remove any faecal matter and its surface was scraped without excessive pressure using a sterile slide.

Purelink Genomic DNA mini-kit (Invitrogen) was used to extract DNA from 200 mg ileum mucosa (IM) while Purelink Microbiome DNA Extraction kit (Invitrogen) was used for Ileum content (200 mg) and faeces (200 mg) samples. Both kits were used in accordance with manufacturer's instructions. After extraction, DNA was quantified (Nanodrop®, ThermoFisher) and the presence of microbial DNA was confirmed by 16S rRNA amplification by conventional PCR³⁸.

**16S rRNA amplicon sequencing.** All samples were prepared for MiSeq compositional sequencing using the specifications outlined by Illumina Inc. (Illumina Inc., Cambridge, UK). The V3–V4 region of the 16S rRNA gene was amplified and Illumina index primers attached in two separate PCR reactions³⁹. All PCR reaction conditions and clean up procedures using AMPure XP (Labplan, Kildare, Ireland) were followed as outlined by Illumina Inc. Quantified samples were then sequenced using an Illumina MiSeq system at the Teagasc Sequencing Centre (Moorepark, Fermoy, Ireland).

**Bioinformatic processing and analysis.** Raw sequence reads generated by MiSeq sequencing were processed using the version 1.9.1 of the Quantitative Insights Into Microbial Ecology (QIIME) pipeline⁴⁰ by using the subsampled open-reference OTU calling approach⁴¹. Demultiplexing and trimming of Miseq reads was done with the default QIIME parameters⁴². After trimming, the reads were merged into a single fasta file and clustered into OTU against the GreenGenes database⁴³ (release 2013-08: gg_13_8_otus) by using the parallel *uclust_ref* method. Reads that failed this step were clustered to *de novo* OTU with the *uclust*⁴⁴ method. The filtering of chimeric OTU was performed by using ChimeraSlayer⁴⁵ against the GreenGenes reference alignment. After removing singletons and doubletons OTUs, only those OTUs representing more than 0.005% of the total filtered were retained as suggested by Bokulich *et al.*⁴². For the genera analyses, OTU were collapsed into genera taxonomic level by using the *tax_glam* function in PhyloSeq⁴⁶. The full data sets have been submitted to NCBI Sequence Read Archive (SRA) under accession: SRP111505, (BioProject:PRJNA393762).

**Functional metagenomic predictions.** A functional prediction based on 16S rRNA marker gene sequences was performed using PICRUSt⁴⁷. After excluding the unknown OTUs from the GreenGenes 13.5 reference database and normalizing by 16S rRNA gene copy number, functional metagenomes for each sample were predicted from the Kyoto Encyclopedia of Genes and Genomes (KEGG) catalogue and categorised to a specified KEGG level. In addition, Kegg orthology groups (KOs) were mapped to KEGG and visualized using the Interactive Pathway Explorer (iPath2.0) web-based tool⁴⁸.
**Detection of metabolites by metabolomics.** Detection of metabolites was performed from 160 mg of ileum content samples from Day 0 and Day 2 pigs by liquid chromatography time-of-flight mass spectrometry (LC-QTOF MS/MS) and Gas chromatography time-of-flight mass spectrometry (GC-QTOF MS). Prior to the analysis, an extraction of the metabolites was carried out with methanol and dichloromethane.

LC-QTOF MS/MS was run in Agilent 1200 Series LC system coupled to an Agilent 6540 UHD Accurate-Mass QTOF hybrid mass spectrometer (Santa Clara, CA, USA). Chromatographic separation was performed by using a Poroshell 120 EC-C18 column (50 mm × 2.1 mm i.d., 2.7 µm particle size) from Agilent Technologies. The mobile phases were 95:5 water:ACN (phase A) and 95:5 ACN:water (phase B), both containing 0.1% (v/v) formic acid and 5 mM ammonium acetate as ionization agents. The LC pump was programmed at a flow rate of 0.25 mL min⁻¹ and the elution gradient was as follows: from min 0 to 30, the percentage of phase B was modified from 50% to 100%, and then the final percentage was held for 20 min. Thus, the total analysis time per sample was 61 min (including postprocessing). Accurate mass spectra in MS scan were acquired in the m/z range 40–1200. MassHunter Workstation software (version B7.00 Qualitative Analysis, Agilent Technologies, Santa Clara, CA, USA) was used to process all data obtained by LC-QTOF in auto MS/MS mode. Identification of the metabolites was supported on MS and MS/MS information and search in the METLIN MS and MS/MS databases (http://metlin.scripps.edu), the Human Metabolome Database (HMDB, 3.6 version) and the LIPID MAPS website (http://www.lipidmaps.org), using in all cases the MFs obtained in the previous step. A table with the peak area of all identified compounds in the different samples injected was obtained as a result.

GC-QTOF/MS analysis was performed by Agilent 7890 A Series GC system coupled to an Agilent 7200 UHD Accurate-Mass QTOF hybrid mass spectrometer equipped with an electron impact (EI) source (Santa Clara, CA, USA). The samples were analyzed by EI ionization mode at 70 eV. Chromatographic separation was carried out with a fused silica DB-5MS-UI 30 m × 0.25 mm i.d., × 0.25 µm film thickness capillary column. The GC oven temperature program started at 60 °C (1 min held), followed by a temperature ramp of 10 °C min⁻¹ to final 300 °C (2 min held). Post-run time was programmed for 4 min up to 310 °C to assure complete elution of the injected sample. MassHunter Workstation software (version B7.00 Qualitative Analysis, Agilent Technologies, Santa Clara, CA, USA) was used to process all data obtained by GC-TOF/MS in full scan mode. The list of MFs obtained for each analysis was exported as data files in compound exchange format (.cef files). Tentative identification of compounds was performed by searching each mass spectrum in the NIST 11 and Fiehn databases using the retention index or retention time value, respectively.

**Statistical analysis.** Further analyses were performed in R v3.1.1. Microbiota and study variables (type of sample, sampling, infection) were included in the estimation of alpha diversity (Shannon and Chao1 indexes) and beta diversity (Simpson index) analysis by Vegan R package⁴⁹. Potential differences in richness of factors included in the study were estimated by analysis of variance (ANOVA). Dissimilarities between pairs of samples were estimated with the Bray-Curtis dissimilarity index⁵⁰ and weighted unifrac index⁵¹ and analysed with non-linear multi-dimensional scaling (NMDS) in Vegan. The Vegan envfit function was used to evaluate if the factors of study (i.e. sampling day and infection status) where associated to the NMDS ordinations; the significance of the fitted factors was estimated by using 999 permutations.

Normalisation of OTUs counts is recommended for further statistical processing⁵². The metagenomeSeq R package⁵³ was used for data normalisation and further statistical analysis of differences in OTUs abundance. By using the *fitZig* function, a zero-inflated Gaussian mixture model, OTUs counts were normalised by calculating the scaling factors equal to the sum of counts. The parameters for the mixture model were estimated with an expectation-maximization algorithm, coupled with a moderated *t* statistic. Individual pig was used as a covariate in the analyses. Differences in OTUs threshold was fixed at a false discovery rate (FDR) value of 0.05. Finally, KO abundance obtained from PICRUSt was introduced into R using PhyloSeq and statistical differences were estimated using Deseq²⁵⁴ using the Wald-test and a FDR threshold of 0.05. The relationship between histopathological lesions and abundance of bacteria in the ileum mucosa was explored by non-parametric Spearman's rank correlation coefficient test implemented in Qiime.

Metabolites abundance were measured by the normalised area under the curve. For each metabolite, comparisons were made by the non-parametric Wilcoxon-test in R with a significant threshold of 0.05.

**Ethical approval.** All procedures involving animals were approved by the institutional bioethical committee to the University of Leon (license number ULE_003_2015) and performed according to European regulations regarding animal welfare and protection of animals used for experimental and other scientific purposes.

**References**

1. EFSA (European Food Safety Authority) and ECDC (European Centre for Disease Prevention and Control) The European Union summary report on trends and sources of zoonoses, zoonotic agents and food-borne outbreaks in 2015. *EFSA J* (2016).

2. Pires, S. M. et al. Source attribution of human salmonellosis: an overview of methods and estimates. *Foodborne Pathog Dis* **11**, 667–676 (2014).

3. Boyen, F. et al. Non-typhoidal Salmonella infections in pigs: a closer look at epidemiology, pathogenesis and control. *Vet Microbiol* **130**, 1–19 (2008).

4. Arguello, H. et al. Prevalence, serotypes and resistance patterns of Salmonella in Danish pig production. *Res Vet Sci* **95**, 334–342 (2013).

5. García-Feliz, C. et al. Salmonella enterica infections in Spanish swine fattening units. *Zoonoses Public Health* **54**, 294–300 (2007).

6. Fedorka-Cray, P. J. et al. Transmission of Salmonella typhimurium to swine. *Vet Microbiol* **41**, 333–344 (1994).

7. Collado-Romero, M. et al. Quantitative analysis of the immune response upon Salmonella typhimurium infection along the porcine intestinal gut. *Vet Res* **41**, 23 (2010).

8. Uribe, J. H. et al. Transcriptional analysis of porcine intestinal mucosa infected with Salmonella Typhimurium revealed a massive inflammatory response and disruption of bile acid absorption in ileum. *Vet Res* **47**, 11 (2016).

9. Argüello, H. et al. Effect of the addition of organic acids in drinking water or feed during part of the finishing period on the prevalence of Salmonella in finishing pigs. *Foodborne Pathog Dis* **10**, 842–849 (2016).
10. Arguello, H. et al. Evaluation of protection conferred by a Salmonella Typhimurium inactivated vaccine in Salmonella-infected finishing pig farms. *Comp Immunol Microbiol Infect Dis* **36**, 489–498 (2013).

11. Hurd, H. S. et al. Rapid infection in market-weight swine following exposure to a Salmonella typhimurium-contaminated environment. *Am J Vet Res* **62**, 1194–1197 (2008).

12. Collado-Romero, M. et al. An in vivo proteomic study of the interaction between Salmonella Typhimurium and porcine ileum mucosa. *J Proteomics* **75**, 2015–2026 (2012).

13. Martins, R. P. et al. Exploring the immune response of porcine mesenteric lymph nodes to Salmonella enterica serovar Typhimurium: an analysis of transcriptional changes, morphological alterations and pathogen burden. *Comp Immunol Microbiol Infect Dis* **36**, 149–160 (2013).

14. Mach, N. et al. Early-life establishment of the swine gut microbiome and impact on host phenotypes. *Environ Microbiol Rep* **7**, 554–569 (2015).

15. Ramayo-Caldas, Y. et al. Phylogenetic network analysis applied to pig gut microbiota identifies an ecosystem structure linked with growth traits. *ISME J* **10**, 2973–2977 (2016).

16. Xiao, L. et al. A reference gene catalogue of the pig gut microbiome. *Nat Microbiol*, **16**161 (2016).

17. Bearson, S. M. et al. Profiling the gastrointestinal microbiota in response to Salmonella: low versus high Salmonella shedding in the natural porcine host. *Infect Genet Evol* **16**, 330–340 (2013).

18. Drumo, R. et al. Salmonella enterica Serovar Typhimurium Exploits Inflammation to Modify Swine Intestinal Microbiota. *Front Cell Infect Microbiol* **5**, 106 (2015).

19. Borewicz, K. A. et al. Changes in the Porcine Intestinal Microbiome in Response to Infection with Salmonella enterica and Lawsonia intracellularis. *PLoS One* **10**, e0139106 (2015).

20. Lynch, H. et al. Evaluation of an Alternative Experimental Infection Method, Which Closely Mimics the Natural Route of Transmission of Monophasic Salmonella Typhimurium in Pigs. *Foodborne Pathog Dis* **14**, 23–28 (2017).

21. Gart, E. V. et al. Salmonella Typhimurium and Multidirectional Communication in the Gut. *Front Microbiol* **7**, 1827 (2016).

22. Borton, M. A. et al. Chemical and pathogen-induced inflammation disrupt the murine intestinal microbiome. *Microbiome*. **5**, 47 (2017).

23. Thiennimitr, P. et al. Intestinal inflammation allows Salmonella to use ethanolamine to compete with the microbiota. *Proc Natl Acad Sci USA* **108**, 17480–17485 (2011).

24. Raffatellu, M. et al. Lipocalin-2 resistance confers an advantage to Salmonella enterica serotype Typhimurium for growth and survival in the inflamed intestine. *Cell Host Microbe* **5**, 476–486 (2009).

25. Holman, D. B. et al. Meta-analysis To Define a Core Microbiota in the Swine Gut. *mSystems*, **2** (2017).

26. Zhao, W. et al. The dynamic distribution of porcine microbiota across different ages and gastrointestinal tract segments. *PLoS One* **10**, e0117441 (2015).

27. Alvarez-Ordóñez, A. et al. Salmonella spp. survival strategies within the host gastrointestinal tract. *Microbiology* **157**, 3268–3281 (2011).

28. Kurdi, P. et al. Mechanism of growth inhibition by free bile acids in lactobacilli and bifidobacteria. *J Bacteriol* **188**, 1979–1986 (2006).

29. Ruiz, L., Margolles, A. & Sánchez, B. Bile resistance mechanisms in Lactobacillus and Bifidobacterium. *Front Microbiol* **4**, 396 (2013).

30. Derakhshani, H. et al. The Features of Fecal and Ileal Mucosa-Associated Microbiota in Dairy Calves during Early Infection with Mycobacterium avium Subspecies paratuberculosis. *Front Microbiol* **7**, 426 (2013).

31. Singh, P. et al. Intestinal microbial communities associated with acute enteric infections and disease recovery. *Microbiome*. **3**, 45 (2015).

32. Stecher, B. et al. Salmonella enterica serovar typhimurium exploits inflammation to compete with the intestinal microbiota. *PLoS Biol* **5**, 2177–2189 (2007).

33. Stecher, B. et al. Like will to like: abundances of closely related species can predict susceptibility to intestinal colonization by pathogenic and commensal bacteria. *PLoS Pathog* **6**, e1000711 (2010).

34. Behnsen, J. et al. Exploiting host immunity: the Salmonella paradigm. *Trends Immunol* **36**, 112–120 (2015).

35. Baughn, A. D. & Malamy, M. H. The strict anaerobe Bacteroides fragilis grows in and benefits from nanomolar concentrations of oxygen. *Nature* **427**(6973), 441–444 (2004).

36. Munyaka, P. M. et al. Carrageenan Gum and Adherent Invasive Escherichia coli in a Piglet Model of Inflammatory Bowel Disease: Impact on Intestinal Mucosa-associated Microbiota. *Front Microbiol* **7**, 462 (2016).

37. Martins, R. P. et al. Pyroptosis and adaptive immunity mechanisms are promptly engendered in mesenteric lymph-nodes during pig infections with Salmonella enterica serovar Typhimurium. *Vet Res* **44**, 120 (2013).

38. Martin-Peláez, S. et al. Evaluation of the inclusion of a mixture of organic acids or lactulose into the feed of pigs experimentally challenged with Salmonella Typhimurium. *Vet Microbiol* **142**, 337–345 (2010).

39. Klindworth, A. et al. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. *Nucleic Acids Res* **41**, e1 (2013).

40. Caporaso, J. et al. QIIME allows analysis of high-throughput community sequencing data. *Nat Methods* **7**, 335–336 (2010).

41. Rideout, J. R. et al. Subsampled open-reference clustering creates consistent, comprehensive OTU definitions and scales to billions of sequences. *Peer J* **2**, e545 (2014).

42. Bokulich, N. A. et al. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. *Nat Methods* **10**, 57–59 (2013).

43. DeSantis, T. Z. et al. Greengenes, a Chimera-Checked 16S rRNA Gene Database and Workbench Compatible with ARB. *Appl Env Microb* **72**, 5069–5072 (2006).

44. Edgar, R. C. Search and clustering orders of magnitude faster than BLAST. *Bioinformatics* **26**, 2460–2461 (2010).

45. Haas, B. J. et al. Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. *Genome Res* **21**, 494–504 (2011).

46. McMurdie, P. J. & Holmes, S. Phyloseq An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. *PLoS ONE* **8**, e61217 (2013).

47. Langille, M. G. I. et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. *Nat Biotech* **31**, 814–821 (2013).

48. Yamada, T. et al. iPath2.0: interactive pathway explorer. *Nucleic Acids Research* **39**, W412–W415 (2011).

49. Oksanen, J. et al. vegan: Community Ecology Package. *R package version* **2**, 0–9 (2009).

50. Bray, J. R. & Curtis, J. T. An ordination of upland forest communities of southern Wisconsin. *Ecological Monographs*. 325–349 (1957).

51. Lozupone, C. A. et al. Quantitative and qualitative beta diversity measures lead to different insights into factors that structure microbial communities. *Appl Environ Microbiol* **73**, 1576–1585 (2007).

52. Weiss, S. et al. Normalization and microbial differential abundance strategies depend upon data characteristics. *Microbiome* **5**, 27 (2017).

53. Paulson, J. N. et al. Differential abundance analysis for microbial marker-gene surveys. *Nat Meth* **10**, 1200–1202 (2013).

54. Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq 2. *Genome Biol* **15**, 550 (2014).
## Acknowledgements

We want to acknowledge the staff from the Genomics and Animal breeding group at the University of Córdoba for their technical support to carry out this study. This article is based upon work from COST Action FA1401 (PiGutNet), supported by COST (European Cooperation in Science and Technology. This work was supported by the Spanish Ministry of Economy and Competitiveness (AGL2014-54089-R/AGL2017-87415-R). HA was funded by the PiGutNet COST action (FA1401) for a Short-Term Scientific Mission at INRA's GABI laboratory (Jouy-en-Josas, France). SZL is a postdoctoral researcher supported by the Postdoctoral Trainee Program of the Spanish Ministry of Economy and Competitiveness (FPDI-2013-15619). HA is a postdoctoral researcher supported by the Juan de la Cierva Postdoctoral Trainee Program of the Spanish Ministry of Economy and Competitiveness (FJCI-2014-22877).

## Author Contributions

H.A. performed the sample processing, statistical analysis and manuscript writing. J.E. supervised and performed the statistical analysis and revised the manuscript. S.Z.L. helped in the study design, sequences upload and manuscript writing. A.J.M. helped in the sample processing and manuscript writing. A.C. was responsible for the challenge study with pigs. F.C. performed sequencing analysis. O.O. performed the bioinformatic analysis of the sequences. F.P.C. and A.L.B. performed the metabolomics and the statistical analysis of the metabolome compounds. P.D.C. participated in the revision of the manuscript. L.M. helped in the statistical analysis. J.J.G. designed the study, supervised the challenge study, writing and revised the final manuscript.

## Additional Information

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-018-26083-3.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018