# Transient effect of single or repeated acute deoxynivalenol and zearalenone dietary challenge on fecal microbiota composition in female finishing pigs

M. Le Sciellour, Olivier Zemb, Aira Maye Serviento, David Renaudeau

► To cite this version:

M. Le Sciellour, Olivier Zemb, Aira Maye Serviento, David Renaudeau. Transient effect of single or repeated acute deoxynivalenol and zearalenone dietary challenge on fecal microbiota composition in female finishing pigs. Animal, 2020, 14 (11), pp.2277-2287. 10.1017/S1751731120001299. hal-02929238

HAL Id: hal-02929238

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

Submitted on 3 Sep 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
<figure><img src="image_2.png" /></figure>

Transient effect of single or repeated acute deoxynivalenol and
zearalenone dietary challenge on fecal microbiota composition in
female finishing pigs

M. Le Sciellour¹, O. Zemb², A.-M. Serviento¹ and D. Renaveau¹†

¹PEGASE, INRAE, Agrocampus Ouest, 16 Le Clos Domaine de la Prise, Saint-Gilles 35590, France; ²GenPhySE, Université de Toulouse, INRA, INPT, INP-ENVT, 24, chemin de Borde-Rouge - Auzeville Tolosane, Castanet Tolosan 31320, France

(Received 18 September 2019; Accepted 22 May 2020; First published online 18 June 2020)

Mycotoxins are a major contaminant of pig feed and have negative effects on health and performance. The present study investigated the impact of single or repeated acute challenges with a diet naturally contaminated with deoxynivalenol (DON) and zearalenone (ZEN) on growth performances of finishing pigs and their fecal microbiota composition. A total of 160 pigs (castrated males and females) in two successive batches were randomly divided into four experimental groups of 40 pigs each. The control group received a control finisher diet from 99 to 154 days of age. Challenged groups were subjected to a 7-day acute challenge by being fed a DON- and ZEN-contaminated diet (3.02 mg DON/kg feed and 0.76 mg ZEN/kg feed) at 113 days (group DC), 134 days (group CD) or both 113 and 134 days (group DD). Microbiota composition was analyzed via 16S rRNA sequencing from fecal samples collected from the 80 females at 99, 119, 140 and 154 days. Challenged pigs (i.e. groups DC, CD and DD) reduced their average daily feed intake by 25% and 27% (P < 0.001) and feed efficiency by 34% and 28% (P < 0.05) during the first and second mycotoxin exposure, respectively. Microbiota composition was affected by mycotoxin exposure (P = 0.07 during the first exposure and P = 0.01 during the second exposure). At the family level, mycotoxin exposure significantly (P < 0.05) decreased the relative abundances of Ruminococcaceae, Streptococcaceae and Veillonellaceae and increased that of Erysipelotrichaceae at both 119 and 140 days of age. After the 7-day DON/ZEN challenge, the relative abundance of 6 to 148 operational taxonomic units (OTUs) differed among the treatment groups. However, none of these OTUs changed in all treatment groups. Using 27 functional pathways, pigs exposed to DON/ZEN challenges could be distinguished from control pigs using sparse partial least squares discriminant analysis, with a 15% misclassification rate. Regarding the functionality of these predictors, two pathways were involved in detoxifying mycotoxins: drug metabolism and xenobiotic metabolism by cytochrome P450. In challenged pigs, microbiota composition returned to the initial state within 3 weeks after the end of a single or repeated DON/ZEN challenge, highlighting the resilience of the gut microbiome. The feeding and growth performances of the pigs during challenge periods were significantly correlated with biological pathways related to health problems and modifications in host metabolism. To conclude, short-term DON/ZEN challenges resulted in transient modifications in the composition and functions of fecal microbiota.

**Keywords:** swine, microbiome, mycotoxin, performance, acute challenge

Implications

The present study aimed to investigate short- and long-term
effects of a mycotoxin challenge (based on deoxynivalenol
and zearalenone) on microbiota composition in pigs and
to look at the relationships between these changes and host
performance during and after these challenges. The study
also demonstrated the transient impact of an acute challenge
on fecal microbiota composition, suggesting that pigs in a

commercial context may recover quickly from this type of
disturbance.

Introduction

Deoxynivalenol (DON) and zearalenone (ZEN) are secondary
metabolites produced by *Fusarium* spp. of fungi. These
mycotoxins are known in livestock production since they
impair animal health and growth when ingested via
contaminated cereals, such as wheat and maize, in feed
(Bouhet and Oswald, 2005; Antonissen *et al.*, 2014). In

<footnote>† E-mail: david.renaveau@inrae.fr</footnote>
particular, exposure of pigs to DON injures the epithelium, which compromises their immune functions, results in an inflammatory response and may cause diarrhea (Bouhet and Oswald, 2005; Sobrova et al., 2010). This results in reduced feed intake and growth rate (Andretta et al., 2012; Serviento et al., 2018). Modifications in pig microbiota composition have been observed after chronic ingestion of the mycotoxins fumonisin (Mateos et al., 2018) and DON (Waché et al., 2008). In the latter study, chronic exposure to DON resulted in transient modifications in the abundance of fecal aerobic mesophilic bacteria during the first week of the DON challenge (Waché et al., 2008). However, short- and long-term impacts of single or repeated acute DON/ZEN challenges on pig fecal microbiota composition have not yet been described. *In vitro* studies demonstrated that the microbiota could detoxify DON (He, 1992; Kollarczik et al., 1994; Eriksen et al., 2002 and 2003; Young et al., 2007). Thus, microbiota could help hosts to continue functioning despite mycotoxin exposure via feed. Consequently, we assume that animals' ability to cope with a mycotoxin challenge may be related to microbiota composition and modifications of it. The present study first evaluated the impact of single or repeated acute DON/ZEN challenges on short- and long-term changes in fecal microbiota. It then explored the relationship between microbiota functions and variations in pig performance due to the DON/ZEN challenge.

## Materials and methods

### Animals and experimental design

The experiment was conducted on Pietrain × (Large White × Landrace) pigs 99 to 154 days old. The 160 pigs (80 castrated males and 80 females) were raised in two batches from April to October 2017 in experimental facilities at INRAE's Unité Expérimentale Porcs de Rennes, located in Saint Gilles, France. Within each batch, pigs were randomly allotted to four experimental groups in a randomized complete block design. Within each batch, all the pigs were housed in one single pen. Each litter and sex was equally represented in each treatment. The experimental room was equipped with a weighing machine and the type of automatic precision feeders described by Pomar et al.

(2011). Serviento et al. (2018) describe the experimental room in detail. During the challenge periods, the living areas of control and challenged pigs were separated to avoid cross-contamination via feces. Pigs had *ad libitum* access to water and feed.

A naturally DON- and ZEN-contaminated feed was used to challenge the pigs. Both the control finisher diet and contaminated diet were based on maize and soybean meal, with 0.8 g standardized ileal digestible lysine per MJ of net energy, and met the ideal profile for essential amino acids. Serviento et al. (2018) describe the diet compositions in detail (online Supplementary Table S1). The control diet was based on maize, with trace DON concentration. The contaminated diet was formulated with naturally contaminated maize containing 4.8 mg DON/kg. According to analyses of mycotoxin composition, the control diet contained 0.14 mg DON/kg and 0.10 mg ZEN/kg, while the contaminated diet contained 3.02 mg DON/kg and 0.76 mg ZEN/kg. In the experimental design (Figure 1), the control group (CC) received the control finisher feed throughout the experiment (from 99 to 154 days of age). The challenged groups were subjected to the 7-day DON/ZEN challenge on day 113 (group DC), day 134 (group CD) or days 113 and 134 (group DD). Before and after the challenge periods, challenged pigs received the control finisher feed.

### Measurements

Live BW was recorded each time a pig passed through the weighing machine and was averaged daily into an average daily BW. The automatic feeders were used to record individual daily feed intake. Feces were collected from females at 99, 119, 140 and 154 days of age. Each fecal sample was taken directly from the rectum, immediately stored in barcode tubes and snap-frozen in liquid N₂ within 10 min of sampling. All samples were stored at -80°C until DNA extraction.

### DNA sequencing

DNA was extracted from 50 mg samples of feces. Bead-beating lysis and DNA purification using the ZR-96 Soil Microbe DNA kit (Zymo Research, Irvine, CA, USA) were conducted according to the manufacturer instructions (Frese et al., 2015). The V3-V4 region of the 16S rRNA gene was

<figure><img src="image_3.png" /><figcaption>**Figure 1** (colour online) Experimental design. The pigs were given a control diet all over the experimental period (group CC) except for DC, CD and DD groups fed a deoxynivalenol- and zearalenone-contaminated diet (DON) during P2 (113 to 119 days), P4 (134 to 140 days) and both P2 and P4, respectively.</figcaption></figure>
amplified through PCR using F460 and R460 primers (F460: CTTTCCCTACACGACGCTCTTCCGATCTACGGRAGGCAGCAG, R460: GGAGTTCAGACGTGTGCTCTTCCGATCTTACCAGGGT-ATCTAATCCT) and 30 cycles of PCR with an annealing temperature of 65°C. After sequencing on the Illumina MiSeq platform, the 2 × 250 bp paired-end sequences were cleaned internally for length, homopolymers and undetermined nucleotides. Sequencing and sequence cleaning were performed in the Get-PlaGE platform (Toulouse, France).

Chimeras were removed from the sequences using VSEARCH (Rognes et al., 2016). The resulting sequences were clustered de novo into operational taxonomic units (OTUs) using a similarity threshold of 0.97, and the abundance of each OTU in each sample was recorded. Each OTU was assigned to a taxon based on the SILVA database using USEARCH (9.2.64_i86linux32 version). The PICRUSt pipeline was used to categorize KEGG (Kyoto Encyclopedia of Genes and Genomes) orthologs resulting from the 16S rRNA sequencing data into biological pathways.

## Calculations and statistical analysis

Statistical analyses were conducted using R 3.5.1. (R Core Team, 2018). Average daily gain (ADG), average daily feed intake (ADFI) and feed efficiency (FE = ADG/ADFI) of the female pigs were compared using ANOVA models with the following effects: period, treatment and their interaction. Each ANOVA model also included the fixed effect of the batch and the initial BW as a covariable.

Microbiota diversity was assessed based on the number of OTU and Shannon index using the phyloseq package of R (McMurdie and Holmes, 2013). The structure of the bacterial community was analyzed by calculating a Bray–Curtis distance matrix. The homogeneity of the dispersion within each group was checked using a multivariate analogue of Levene's test for homogeneity of variances. A pairwise multi-variate ANOVA (MANOVA) was applied to the distance matrix to compare the control and challenged groups at 99, 119, 140 and 154 days of age. Differences in groups' bacterial communities were considered significant at an adjusted false discovery rate (*P*<sub>FDR</sub) < 0.05 after 999 permutations based on the Benjamini–Hochberg procedure. In a second step, the impacts of the DON/ZEN challenges on microbiota composition at the family and genus levels were assessed using the non-parametric Wilcoxon test by comparing the relative abundance of family and genus of control and challenged groups at the end of each period (i.e. 99, 119, 140 and 154 days of age). Finally, the relative abundance of OTU of control and challenged pigs was compared using generalized linear models (GLMs) developed with the edgeR package of R (Robinson et al., 2010). The GLMs included a batch effect within each treatment. The GLM analyses were conducted as pairwise comparisons between the control group and each treatment group within each period: at 99 days to assess differences before the challenges, at 119 and 140 days to assess short-term effects of DON/ZEN mycotoxins on microbiota composition and at 154 days to assess long-term effects of DON/ZEN mycotoxins. In the

GLM analyses, abundances of an OTU were considered to differ significantly at *P*<sub>FDR</sub> < 0.05.

To extract the biological functions that best distinguished control pigs from challenged pigs at 119 and/or 140 days of age, multivariate integrative sparse partial least squares discriminant analysis (MINT sPLS-DA) was used (Lê Cao et al., 2011). This approach can aggregate functions even when they are provided by different microbial species. Using microbiota information collected at 119 and 140 days as two different datasets, this analysis was applied to the pathway abundance tables, which had been normalized using total-sum scaling and centered log-ratio transformation. The error rates of the resulting model were validated from 10 replicates of 10-fold cross-validation.

The ADG, ADFI and FE were averaged during the DON/ZEN challenges (113 to 119 and 134 to 140 days of age) and corrected for the batch effect and initial BW at 99 days of age using a linear model. The relationship between each pair of performance traits (ADG, ADFI and FE) was analyzed via Pearson correlations, and the correlations were considered significant at *P*<sub>FDR</sub> < 0.05. Regularized canonical correlation analysis (rCCA) was used to analyze correlations between animal performance (ADG, ADFI and FE) and the relative abundance of biological pathways of the microbiota. As a multivariate approach, rCCA calculates correlations between two datasets acquired from the same experimental units when the datasets have more variables than samples (González et al., 2008). To determine a significant threshold for the correlation of each performance-pathway pair, the matrix of the relative abundance of the pathways was randomized by the pathway. The correlation of each performance-pathway pair was then calculated and recorded. This sequence of randomization and correlation was repeated 100 times. For each pair, the threshold of a significant positive or negative correlation was the maximum or minimum correlation, respectively, obtained during 100 randomizations.

## Results

### Animal performance

Impacts of the DON/ZEN challenges on the performances of all pigs (males and females combined) were previously published (Serviento et al., 2018). Thus, we provide results only for the performance of females because microbiota data were collected only from them. Mean (±1 SD) BW was 54.2 ± 5.0 kg at 99 days of age and 106.3 ± 10.0 kg at 154 days. Briefly, when exposed to a DON/ZEN challenge, mean ADFI and ADG decreased significantly. For instance, ADFI decreased from 2.39 to 1.74 kg/day from days 113 to 119, respectively, and 2.91 to 2.17 kg/day from days 134 to 140, respectively (*P* < 0.001) (itemized values are presented in Table 1). However, mean ADFI and ADG did not differ between the control and challenged pigs over the entire period (ADG: 0.98 and 0.94 kg/day, respectively; *P* > 0.05). Mean FE decreased during DON/ZEN challenge periods
<table><caption>**Table 1** Effect of the mycotoxin challenge on the pigs performance</caption><thead><tr><th rowspan="2"></th><th colspan="4">Treatments¹</th><th rowspan="2">RSD²</th><th rowspan="2">Statistics</th></tr><tr><th>CC</th><th>DC</th><th>CD</th><th>DD</th></tr></thead><tbody><tr><td>No. of pigs</td><td>15</td><td>16</td><td>15</td><td>18</td><td></td><td></td></tr><tr><td>ADFI (kg/day)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>P1 (99 to 112 days)</td><td>2.19<sup>u</sup></td><td>2.26<sup>uv</sup></td><td>2.24<sup>u</sup></td><td>2.13<sup>u</sup></td><td>0.37</td><td>P***, P<sup>x</sup>T***</td></tr><tr><td>P2 (113 to 119 days)</td><td>2.43<sup>au</sup></td><td>1.83<sup>bu</sup></td><td>2.35<sup>au</sup></td><td>1.64<sup>bv</sup></td><td></td><td></td></tr><tr><td>P3 (120 to 133 days)</td><td>2.58<sup>uv</sup></td><td>2.72<sup>vw</sup></td><td>2.56<sup>u</sup></td><td>2.55<sup>w</sup></td><td></td><td></td></tr><tr><td>P4 (134 to 140 days)</td><td>2.92<sup>av</sup></td><td>3.10<sup>aw</sup></td><td>2.28<sup>bu</sup></td><td>2.07<sup>buv</sup></td><td></td><td></td></tr><tr><td>P5 (141 to 154 days)</td><td>3.01<sup>v</sup></td><td>3.10<sup>w</sup></td><td>3.13<sup>v</sup></td><td>3.18<sup>x</sup></td><td></td><td></td></tr><tr><td>P1 to P5</td><td>2.62</td><td>2.65</td><td>2.57</td><td>2.43</td><td>0.30</td><td>ns</td></tr><tr><td>ADG (kg/day)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>P1</td><td>0.77<sup>u</sup></td><td>0.85<sup>v</sup></td><td>0.78<sup>u</sup></td><td>0.76<sup>u</sup></td><td>0.21</td><td>P***, P<sup>x</sup>T***, P<sup>x</sup>B**</td></tr><tr><td>P2</td><td>1.10<sup>av</sup></td><td>0.52<sup>bu</sup></td><td>1.01<sup>avw</sup></td><td>0.49<sup>bv</sup></td><td></td><td></td></tr><tr><td>P3</td><td>0.97<sup>u</sup></td><td>1.11<sup>vw</sup></td><td>1.02<sup>vw</sup></td><td>1.08<sup>w</sup></td><td></td><td></td></tr><tr><td>P4</td><td>1.16<sup>av</sup></td><td>1.19<sup>aw</sup></td><td>0.65<sup>bu</sup></td><td>0.65<sup>buv</sup></td><td></td><td></td></tr><tr><td>P5</td><td>1.01<sup>u</sup></td><td>1.07<sup>vw</sup></td><td>1.15<sup>w</sup></td><td>1.15<sup>w</sup></td><td></td><td></td></tr><tr><td>P1 to P5</td><td>0.98</td><td>0.97</td><td>0.95</td><td>0.89</td><td>0.12</td><td>ns</td></tr><tr><td>FE</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>P1</td><td>0.36<sup>u</sup></td><td>0.38<sup>u</sup></td><td>0.35<sup>uv</sup></td><td>0.36<sup>uv</sup></td><td>0.08</td><td>P<sup>x</sup>T***, P<sup>x</sup>B**</td></tr><tr><td>P2</td><td>0.46<sup>av</sup></td><td>0.27<sup>bv</sup></td><td>0.42<sup>au</sup></td><td>0.30<sup>bu</sup></td><td></td><td></td></tr><tr><td>P3</td><td>0.38<sup>uv</sup></td><td>0.41<sup>u</sup></td><td>0.40<sup>u</sup></td><td>0.42<sup>v</sup></td><td></td><td></td></tr><tr><td>P4</td><td>0.40<sup>auv</sup></td><td>0.39<sup>au</sup></td><td>0.28<sup>bv</sup></td><td>0.31<sup>abu</sup></td><td></td><td></td></tr><tr><td>P5</td><td>0.33<sup>u</sup></td><td>0.34<sup>uv</sup></td><td>0.37<sup>uv</sup></td><td>0.36<sup>uv</sup></td><td></td><td></td></tr><tr><td>P1 to P5</td><td>0.37</td><td>0.37</td><td>0.37</td><td>0.37</td><td>0.03</td><td>ns</td></tr></tbody></table>

ADFI = average daily feed intake; ADG = average daily gain; FE = feed efficiency.

a,bLeast square means within a row with different superscript differ according to the treatment group.

u,v,w,x¹Least square means within a column with different superscript differ according to the period.

¹Control pigs (CC) were fed a control diet. Pigs in the challenged groups were fed a deoxynivalenol- and zearalenone-contaminated diet between 113 and 119 days (DC), between 134 and 140 days (CD) and both between 113 and 119 days and between 134 and 140 days (DD).

²Residual SD from an ANOVA model accounting for the period (P), the treatment (T), the batch (B) and their interactions. The initial BW was included as a covariate in the models. **P < 0.01, ***P < 0.001, ns = not significant.

(from 0.44 to 0.29 from days 113 to 119, respectively, and from 0.40 to 0.29 from days 134 to 140, respectively; P < 0.05). However, mean FE did not differ between control and challenged pigs during recovery periods (i.e. 0.39 v. 0.41, respectively, from 119 to 134 days of age; 0.34 and 0.37, respectively, from 140 to 154 days of age; P > 0.05).

### Microbiota composition

A total of 277 fecal samples were collected from 74 female pigs. Six females were removed from the experiment due to health problems or their inability to adapt to the sorting machine or automatic feeders. In addition, we were not able to collect fecal samples from some female pigs. The number of sequences in the samples ranged from 15 924 to 76 582 (mean ± 1 SD = 29 482 ± 11 282). No sample was discarded due to having too few sequences. After filtering out rare OTUs (<0.01% of all sequences in the dataset), 1556 OTUs remained in the dataset.

At 99 days of age, mean ± 1 SD of OTU abundance of the fecal samples was 3649 ± 3642, while the mean ± 1 SD Shannon diversity index was 7.13 ± 0.16 (Table 2). Based on the Wilcoxon tests applied to the richness and the Shannon index of the control and challenged pigs at 119 or 140 days of age, the mycotoxin exposure did not impact

microbial diversity (P > 0.05) (Table 2). *Firmicutes* was the most abundant phylum (80.10% of the affiliated sequences), followed by *Bacteroidetes* (15.37%). *Actinobacteria*, *Proteobacteria* and *Spirochaetes* were represented by 0.69%, 0.49% and 0.18% of the sequences, respectively. *Tenericutes* and *Fibrobacteres* each represented less than 0.01% of the sequences. The remaining sequences were unclassified.

According to the Wilcoxon tests, eight families were impacted by the DON/ZEN challenges (Table 2). Compared to the control pigs, mycotoxin exposure from days 113 to 119 decreased the relative abundances of *Ruminococcaceae* (P < 0.05), *Streptococcaceae* (P < 0.001) and *Veillonellaceae* (P < 0.001), and increased those of *Erysipelotrichaceae* (P < 0.001) and *Peptostreptococcaceae* (P < 0.05). Mycotoxin exposure from days 134 to 140 decreased the relative abundances of *Lachnospiraceae* (P < 0.05), *Prevotellaceae* (P < 0.01), *Ruminococcaceae* (P < 0.01), *Streptococcaceae* (P < 0.001) and *Veillonellaceae* (P < 0.01), and increased those of *Desulfovibrionaceae* (P < 0.01) and *Erysipelotrichaceae* (P < 0.05). At the genus level, mycotoxin exposure from days 113 to 119 decreased the relative abundances of *Blautia* (P < 0.05), *Dialister* (P < 0.001), *Fusicatenibacter* (P < 0.05), *Mitsuokella* (P < 0.001), *Olsenella* (P < 0.01),
<table><caption>**Table 2** Diversity indexes (mean ± SD) and relative abundance (mean ± SD) of families and genera in pigs' fecal samples collected in the four different experimental groups<fnref n="1" /> at 99, 119, 140 and 154 days of age</caption><thead><tr><th rowspan="2"></th><th>d99</th><th colspan="3">d119</th><th colspan="3">d140</th><th>d154</th></tr><tr><th>All<sup>2</sup></th><th>Ctrl<sup>3</sup></th><th>DON/ZEN<sup>3</sup></th><th>P<sup>4</sup></th><th>Ctrl<sup>5</sup></th><th>DON/ZEN<sup>5</sup></th><th>P<sup>6</sup></th><th>All<sup>2</sup></th></tr></thead><tbody><tr><td>Number of pigs</td><td>71</td><td>31</td><td>35</td><td></td><td>34</td><td>37</td><td></td><td>72</td></tr><tr><td>Diversity indexes</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Richness (number of operational taxonomic units)</td><td>3649 ± 364</td><td>3653 ± 279</td><td>3507 ± 385</td><td>0.10</td><td>3764 ± 261</td><td>3626 ± 347</td><td>0.10</td><td>3746 ± 295</td></tr><tr><td>Shannon index</td><td>7.13 ± 0.16</td><td>7.14 ± 0.14</td><td>7.13 ± 0.16</td><td>0.91</td><td>7.21 ± 0.07</td><td>7.16 ± 0.13</td><td>0.05</td><td>7.18 ± 0.11</td></tr><tr><td>Relative abundance (%)</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Families</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Desulfovibrionaceae</td><td>0.05 ± 0.06</td><td>0.05 ± 0.04</td><td>0.08 ± 0.06</td><td>0.11</td><td>0.07 ± 0.04</td><td>0.10 ± 0.06</td><td>&lt;0.01</td><td>0.08 ± 0.06</td></tr><tr><td>Erysipelotrichaceae</td><td>0.57 ± 0.40</td><td>0.98 ± 0.54</td><td>1.61 ± 0.81</td><td>&lt;0.001</td><td>1.23 ± 0.67</td><td>1.90 ± 1.22</td><td>&lt;0.05</td><td>1.45 ± 0.94</td></tr><tr><td>Lachnospiraceae</td><td>9.28 ± 2.60</td><td>5.70 ± 2.03</td><td>5.35 ± 2.60</td><td>0.36</td><td>5.45 ± 1.77</td><td>4.35 ± 1.97</td><td>&lt;0.05</td><td>4.24 ± 2.14</td></tr><tr><td>Peptostreptococcaceae</td><td>0.04 ± 0.03</td><td>0.14 ± 0.06</td><td>0.19 ± 0.09</td><td>&lt;0.05</td><td>0.55 ± 0.21</td><td>0.68 ± 0.35</td><td>0.18</td><td>0.62 ± 0.35</td></tr><tr><td>Prevotellaceae</td><td>7.41 ± 4.26</td><td>5.97 ± 4.07</td><td>5.02 ± 4.31</td><td>0.21</td><td>5.49 ± 4.66</td><td>2.74 ± 2.59</td><td>&lt;0.01</td><td>2.61 ± 2.59</td></tr><tr><td>Ruminococcaceae</td><td>4.42 ± 1.54</td><td>2.78 ± 0.91</td><td>2.36 ± 1.45</td><td>&lt;0.05</td><td>2.53 ± 1.16</td><td>1.76 ± 0.84</td><td>&lt;0.01</td><td>1.80 ± 0.96</td></tr><tr><td>Streptococcaceae</td><td>13.76 ± 6.07</td><td>9.86 ± 5.17</td><td>5.29 ± 4.19</td><td>&lt;0.001</td><td>9.60 ± 3.78</td><td>5.49 ± 3.67</td><td>&lt;0.001</td><td>9.21 ± 4.62</td></tr><tr><td>Veillonellaceae</td><td>1.79 ± 1.75</td><td>0.88 ± 1.19</td><td>0.27 ± 0.37</td><td>&lt;0.001</td><td>0.36 ± 0.57</td><td>0.13 ± 0.19</td><td>&lt;0.01</td><td>0.12 ± 0.18</td></tr><tr><td>Genera</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Bacteroides</td><td>0.00 ± 0.00</td><td>0.00 ± 0.00</td><td>0.01 ± 0.05</td><td>0.96</td><td>0.00 ± 0.01</td><td>0.00 ± 0.00</td><td>&lt;0.05</td><td>0.05 ± 0.20</td></tr><tr><td>Blautia</td><td>1.96 ± 0.81</td><td>0.80 ± 0.32</td><td>0.70 ± 0.55</td><td>&lt;0.05</td><td>0.54 ± 0.25</td><td>0.62 ± 0.41</td><td>0.74</td><td>0.51 ± 0.33</td></tr><tr><td>Clostridium IV</td><td>0.02 ± 0.02</td><td>0.04 ± 0.03</td><td>0.06 ± 0.05</td><td>&lt;0.01</td><td>0.06 ± 0.04</td><td>0.06 ± 0.03</td><td>0.61</td><td>0.07 ± 0.04</td></tr><tr><td>Clostridium XI</td><td>0.02 ± 0.02</td><td>0.09 ± 0.05</td><td>0.14 ± 0.08</td><td>&lt;0.01</td><td>0.57 ± 0.31</td><td>0.55 ± 0.25</td><td>0.92</td><td>0.57 ± 0.32</td></tr><tr><td>Dialister</td><td>1.32 ± 1.21</td><td>0.63 ± 0.89</td><td>0.12 ± 0.26</td><td>&lt;0.001</td><td>0.14 ± 0.36</td><td>0.16 ± 0.33</td><td>0.73</td><td>0.08 ± 0.16</td></tr><tr><td>Fusicatenibacter</td><td>0.37 ± 0.20</td><td>0.21 ± 0.13</td><td>0.14 ± 0.14</td><td>&lt;0.05</td><td>0.12 ± 0.12</td><td>0.15 ± 0.12</td><td>0.20</td><td>0.11 ± 0.11</td></tr><tr><td>Megasphera</td><td>0.02 ± 0.02</td><td>0.00 ± 0.01</td><td>0.00 ± 0.00</td><td>&lt;0.05</td><td>0.00 ± 0.00</td><td>0.00 ± 0.00</td><td>0.59</td><td>0.00 ± 0.01</td></tr><tr><td>Mitsuokella</td><td>0.35 ± 0.48</td><td>0.15 ± 0.27</td><td>0.03 ± 0.07</td><td>&lt;0.001</td><td>0.03 ± 0.08</td><td>0.02 ± 0.05</td><td>0.32</td><td>0.01 ± 0.02</td></tr><tr><td>Olsenella</td><td>0.47 ± 1.02</td><td>0.12 ± 0.24</td><td>0.03 ± 0.06</td><td>&lt;0.01</td><td>0.01 ± 0.02</td><td>0.00 ± 0.00</td><td>0.07</td><td>0.00 ± 0.01</td></tr><tr><td>Ruminococcus</td><td>0.23 ± 0.14</td><td>0.36 ± 0.24</td><td>0.23 ± 0.13</td><td>&lt;0.05</td><td>0.37 ± 0.24</td><td>0.34 ± 0.25</td><td>0.59</td><td>0.31 ± 0.28</td></tr><tr><td>Sarcina</td><td>0.19 ± 0.27</td><td>0.21 ± 0.20</td><td>0.56 ± 0.64</td><td>&lt;0.05</td><td>0.23 ± 0.35</td><td>0.29 ± 0.35</td><td>0.38</td><td>0.14 ± 0.26</td></tr><tr><td>Streptococcus</td><td>13.76 ± 6.07</td><td>9.86 ± 5.17</td><td>5.29 ± 4.19</td><td>&lt;0.001</td><td>7.47 ± 4.37</td><td>7.27 ± 4.14</td><td>0.90</td><td>9.21 ± 4.62</td></tr><tr><td>Turicibacter</td><td>0.45 ± 0.40</td><td>0.90 ± 0.56</td><td>1.55 ± 0.83</td><td>&lt;0.001</td><td>1.61 ± 1.24</td><td>1.49 ± 0.90</td><td>0.87</td><td>1.41 ± 0.96</td></tr></tbody></table>

<footnote n="1">¹The treatment groups are: the CC group fed a control diet during the whole experiment, the DC group fed a deoxynivalenol (DON)- and zearalenone (ZEN)-contaminated diet between 113 and 119 days, the CD group fed a DON- and ZEN-contaminated diet between 134 and 140 days, and the DD group fed a DON- and ZEN-contaminated diet between 113 and 119 days and between 134 and 140 days.</footnote>

<footnote n="2">²All pigs were in control conditions at 99 and 154 days of age.</footnote>

<footnote n="3">³At 119 days of age, Ctrl referred to the pigs in the CC and CD groups, and DON referred to the pigs in the DC and DD groups.</footnote>

<footnote n="4">⁴The *P*-value resulted from a Wilcoxon test between Ctrl and challenged pigs at 119 days of age.</footnote>

<footnote n="5">⁵At 140 days of age, Ctrl referred to the pigs in the CC and DC groups, and DON/ZEN referred to the pigs in the CD and DD groups.</footnote>

<footnote n="6">⁶The *P*-value resulted from a Wilcoxon test between Ctrl and challenged pigs at 140 days of age.</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 2** (colour online) Bray–Curtis distance between treatment groups and controls at 99, 119, 140 and 154 days of age represented in a non-metric multidimensional scaling (nMDS). The control group (CC) was compared with treatment groups: group of pigs fed a deoxynivalenol (DON)- and zearalenone (ZEN)-contaminated diet between 113 and 119 days (DC), group fed a DON- and ZEN-contaminated diet between 134 and 140 days (CD), group fed a DON- and ZEN-contaminated diet between 113 and 119 days and between 134 and 140 days (DD). The *P*-values are presented in Table 1.</figcaption></figure>

*Ruminococcus* (*P* < 0.05) and *Streptococcus* (*P* < 0.001), and increased those of *Clostridium* IV (*P* < 0.01), *Clostridium* XI (*P* < 0.01), *Sarcina* (*P* < 0.05) and *Turicibacter* (*P* < 0.001). Mycotoxin exposure from days 134 to 140 did not significantly influence the relative abundance of genera. However, the relative abundance of two genera (*Lactobacillus* and *Collinsella*) and three families (*Lactobacillaceae*, *Coriobacteriaceae* and *Spirochaetaceae*) differed significantly between groups CD and DD at 140 days (online Supplementary Table S2).

According to the MANOVA, microbiota composition did not differ among control or treatment groups (*P*<sub>FDR</sub> > 0.52) at the beginning of the experiment (99 days of age) (Figure 2, Table 3). At the end of the first challenge period, microbiota composition of challenged pigs tended to differ from that of control pigs (*P* = 0.07). At the end of the second challenge period, microbiota composition of challenged pigs differed significantly from that of control pigs (*P* = 0.01). After 7 days of DON/ZEN challenge, microbiota composition did not differ between groups DC and DD at 119 days (*P*<sub>FDR</sub> = 0.61) or between groups CD and DD at 140 days (*P*<sub>FDR</sub> = 0.43). After a 2-week recovery period, at

154 days, microbiota composition did not differ among control or treatment groups (*P*<sub>FDR</sub> = 0.99).

According to the MANOVA, the global microbiota composition recovers within 2 weeks after the end of a DON/ZEN challenge. Then, we considered that the GLM analyses could be applied at each age without introducing unwanted biases. In the GLM analysis, relative OTU abundance did not differ between control and challenged pigs at 99 days (Figure 3). At 119 days, relative abundance of 148 of the 1556 OTUs in the dataset differed between groups CC and DC and 6 between groups CC and DD. At 140 days, relative abundance of 15 OTUs differed between groups CC and CD, 32 between groups CC and DD and 2 between groups CD and DD. None of the OTUs that differed in relative abundance after the first challenge period (day 119) also differed after the second challenge period (day 140), and vice versa. After the recovery period, at 154 days, relative OTU abundance did not differ between control and challenged pigs.

Based on pathway abundance, the MINT sPLS-DA distinguished the pigs that received the contaminated diet from control pigs using 27 pathways, with a 15% misclassification

<table><caption>**Table 3** Adjusted false discovery rate *P*-values (*P*<sub>FDR</sub>) using Benjamini–Hochberg method resulting from 999 iterations of a pairwise multivariate ANOVA between control (CC) and treatment groups of pigs exposed to a deoxynivalenol- and zearalenone-contaminated diet</caption><thead><tr><th rowspan="2">Days of age</th><th rowspan="2">Number of pigs</th><th colspan="5">P<sub>FDR</sub></th></tr><tr><th>DC v. CC</th><th>CD v. CC</th><th>DD v. CC</th><th>DC v. DD</th><th>CD v. DD</th></tr></thead><tbody><tr><td>d99</td><td>71</td><td>0.52</td><td>0.52</td><td>0.68</td><td>0.72</td><td>0.58</td></tr><tr><td>d119</td><td>66</td><td>0.07</td><td>0.82</td><td>0.07</td><td>0.61</td><td>0.06</td></tr><tr><td>d140</td><td>68</td><td>0.49</td><td>0.01</td><td>0.01</td><td>0.01</td><td>0.43</td></tr><tr><td>d154</td><td>72</td><td>0.99</td><td>0.99</td><td>0.99</td><td>0.99</td><td>0.99</td></tr></tbody></table>

Pigs in the treatment groups were fed a deoxynivalenol- and zearalenone-contaminated diet between 113 and 119 days (DC), between 134 and 140 days (CD) and both between 113 and 119 days and between 134 and 140 days (DD).
<figure><img src="image_5.png" /><figcaption>**Figure 3** (colour online) Operational taxonomic units differentially abundant between control group (CC) and treatment groups: group of pigs fed a deoxynivalenol (DON)- and zearalenone (ZEN)-contaminated diet between 113 and 119 days (DC), group fed a DON- and ZEN-contaminated diet between 134 and 140 days (CD), group fed a DON- and ZEN-contaminated diet between 113 and 119 days and between 134 and 140 days (DD). Log2 fold changes (logFC) resulted from a 2 × 2 comparison using a generalized linear model analysis (*P* < 0.05).</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 4** (colour online) Multivariate integrative sparse partial least squares discriminant analysis based on pathways abundance at 119 days (○) and 140 days (Δ). Pigs receiving a deoxynivalenol (DON)- and zearalenone (ZEN)-contaminated diet (in red) were discriminated from the control pigs (in black) with 27 pathways and 15% error rate of misclassification. Expl. var. refers to *t* amount of variation explained per component.</figcaption></figure>

rate (Figure 4). Based on loading values, the main pathways that contributed to the first component of the sPLS-DA were lipid metabolism, amyotrophic lateral sclerosis, amoebiasis, RIG-I-like receptor signaling, and carbohydrate digestion and absorption (online Supplementary Table S3). The main pathways that contributed to the second component were retinol metabolism, drug metabolism by cytochrome P450, xenobiotic metabolism by cytochrome P450, biosynthesis

of siderophore group non-ribosomal peptides and inositol phosphate metabolism (online Supplementary Table S3).

*Relationships between performance and microbial functioning*

For correlations among performance traits, ADG correlated positively with ADFI from days 113 to 119 (P2) for control pigs (*r* = 0.63, *P* < 0.001) and challenged pigs (*r* = 0.49,
$P<0.05$), and from days 134 to 140 (P4) for control pigs ($r=0.46, P<0.01$) and challenged pigs ($r=0.68, P<0.001$).

Based on the rCCA, correlation coefficients between performance traits of control or challenged pigs during P2 and P4 and the relative abundance of microbiota pathways depended on the trait and the period. Microbiota pathways were weakly correlated with ADG or FE of the control and challenged pigs (Table 4). Compared to P2, fewer correlations were significant during P4. During P2, ADFI was significantly correlated with 11 pathways in the challenged pigs, but with no pathways in the control pigs (Table 4). The 11 pathways corresponded to functions related to environmental information processing, the organismal system, human diseases and metabolism (terpenoid and polyketide metabolism, carbohydrate and glycan metabolism, and xenobiotic biodegradation and metabolism). No significant correlation with ADFI was found during P4 during the second challenge. ADG and FE were not correlated with the relative abundance of pathways during the challenge periods. Feed efficiency of control pigs was correlated with two pathways during P2, but ADFI and ADG showed no significant correlations. However, for control pigs during P4, ADFI was correlated with three pathways and ADG with one pathway. No pathway was correlated with any performance indicator of either control or challenged pigs (Table 4).

## Discussion

Because the impact of DON/ZEN challenges on pig performance was previously described (Serviento et al., 2018), the present study focused mainly on the challenges' impacts on gut microbiota and the latter's relationship with pig performance. Regardless of age at exposure to DON and ZEN, the relative abundances of the families *Streptococcaceae*, *Ruminococcaceae* and *Veillonellaceae* decreased and that of *Erysipelotrichaceae* increased. These families were consistently impacted by the DON/ZEN challenge. In contrast, Reddy et al. (2018a) observed no difference in these families in 8-week-old castrated pigs subjected to a 4-week DON challenge. This difference from the present study could be related to the difference in sex (castrated males *v*. females), age of the pigs (8 *v*. 17 weeks old), breed or environmental factors such as diet or location. At the sub-family level, the OTUs or genera in the present study that were impacted by the DON/ZEN challenge at 119 days of age differed from those impacted at 140 days of age. This confirmed that the effects of acute DON/ZEN exposure on microbiota composition are time-dependent.

The present study was designed to evaluate the effect of repeated DON/ZEN challenges on microbiota composition. According to its results, the relative abundance of only two OTUs differed between groups CD and DD at 140 days of age. In addition, the relative abundance of only two genera (*Lactobacillus* and *Collinsella*) and their corresponding families differed significantly between groups CD and DD

at 140 days of age. This result was confirmed by the MANOVA, which suggests that no habituation mechanisms were observed in microbiota composition when pigs were challenged twice. Accordingly, and based on the ADG and ADFI observed, pigs that had experienced a previous DON/ZEN challenge (group DD) performed similar to those challenged only once (group CD) (Serviento et al., 2018). At a higher DON concentration (>5 mg/kg BW) (Flannery et al., 2011), however, higher ADFI was observed in mice challenged twice. Despite differences in species, diet and other raising conditions, these studies could suggest a dose-dependent effect of DON on host responses and habituation.

In contrast to the effects of chronic ingestion of DON (Robert et al., 2017), our results indicate that acute exposure to a DON- and ZEN-contaminated diet has no detrimental effects on overall pig ADG and FE after a 2-week recovery period. Based on patterns in ADFI, pigs given the DON- and ZEN-contaminated diet started to consume similar amounts than the pigs given the control diet after 3 days of challenge (Serviento et al., 2018). This resulted in a full performance recovery at the end of the trial. With a bit less naturally contaminated diet (2.8 *v*. 3.02 mg DON/kg; 0.28 *v* 0.76 mg ZEN/kg), Waché et al. (2008) observed that pigs recovered full performance after 1 week exposure. Then, we assume that consuming a diet naturally contaminated with DON and ZEN at the concentrations under 3.02 and 0.76 mg/kg, respectively, during 1 week will result in a full performance recovery. Looking at the microbiota composition between control and challenged pigs 1 week after the beginning of the challenge, Waché et al. (2008) did not highlight any difference. Based on our results obtained with a bit more contaminated diet, the number of OTUs that differed in relative abundance between the control pigs and pigs after 1 week of exposure (119 or 140 days of age) ranged from 6 to 148 out of 1556 OTUs in the dataset. This clearly suggests that microbiota composition in the present study also changed moderately after a 7-day exposure to DON and ZEN.

Depending on the level of contamination, 45% to 95% of the DON ingested by pigs end up in the urine without being transformed by the microbiota (Prelusky et al., 1988; Eriksen et al., 2003; Dänicke et al., 2004; Nagl et al., 2014). In other words, gut microbiota have the opportunity to detoxify 5% to 55% of the DON ingested. Similarly, ca. 26% of ZEN end up in the urine (Binder et al., 2017), which could explain in part the moderate changes in fecal microbiota composition observed in the present study during acute mycotoxin exposure. This agrees with previous studies in which fecal microbiota composition changed moderately during chronic DON exposure in pigs (Waché et al., 2008), or not at all in rodents (Saint-Cyr et al., 2013; Payros et al., 2017).

As mentioned, relative bacterial abundance changed moderately by the end of the DON/ZEN challenges, but these differences were no longer present 2 weeks later. This clearly indicates that short-term exposure to DON and ZEN had a short and transient impact on microbiota composition in pigs. The lack of a long-term effect of DON/ZEN contamination
<table><caption>**Table 4** Regularized canonical correlations (rCCA)<fnref n="1" /> between functional pathways relative abundance and the average daily feed intake (ADFI)<fnref n="2" />, the average daily gain (ADG) and the feed efficiency (FE) of pigs</caption><thead><tr><th rowspan="3">Function level A</th><th rowspan="3">Function level B</th><th rowspan="3">Pathway</th><th colspan="4">ADFI</th><th colspan="4">ADG</th><th colspan="4">FE</th></tr><tr><th colspan="2">Control<sup>3</sup></th><th colspan="2">Challenged<sup>4</sup></th><th colspan="2">Control</th><th colspan="2">Challenged</th><th colspan="2">Control</th><th colspan="2">Challenged</th></tr><tr><th>P2<sup>5</sup></th><th>P4<sup>6</sup></th><th>P2</th><th>P4</th><th>P2</th><th>P4</th><th>P2</th><th>P4</th><th>P2</th><th>P4</th><th>P2</th><th>P4</th></tr></thead><tbody><tr><td>Environmental information processing</td><td>Membrane transport</td><td>Secretion system</td><td>—</td><td>—</td><td>-0.52</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Organismal system</td><td>Circulatory system</td><td>Cardiac muscle contraction</td><td>—</td><td>—</td><td>0.43</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td rowspan="3">Human diseases</td><td>Cancers</td><td>Bladder cancer</td><td>—</td><td>—</td><td>-0.47</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Immune diseases</td><td>Primary immunodeficiency</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>0.47</td><td>—</td><td>—</td></tr><tr><td>Neurodegenerative diseases</td><td>Huntington's disease</td><td>—</td><td>—</td><td>-0.44</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td rowspan="6">Metabolism</td><td>Carbohydrate metabolism</td><td>Butanoate metabolism</td><td>—</td><td>—</td><td>-0.48</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Glycan biosynthesis and metabolism</td><td>Other glycan degradation</td><td>—</td><td>—</td><td>0.46</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Lipid metabolism</td><td>Glycosaminoglycan degradation</td><td>—</td><td>—</td><td>—</td><td>—</td><td>0.46</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Metabolism of other amino acids</td><td>Linoleic acid metabolism</td><td>—</td><td>-0.44</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Metabolism of terpenoids and polyketides</td><td>Phosphonate and phosphinate metabolism</td><td>—</td><td>—</td><td>-0.43</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>0.55</td><td>—</td><td>—</td></tr><tr><td>Xenobiotics biodegradation and metabolism</td><td>Limonene and pinene degradation</td><td>—</td><td>—</td><td>-0.49</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td></td><td></td><td>Benzoate degradation</td><td>—</td><td>—</td><td>-0.48</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td></td><td></td><td>Tetracycline degradation</td><td>—</td><td>—</td><td>-0.41</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td></td><td></td><td>Ethylbenzene degradation</td><td>—</td><td>—</td><td>-0.47</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td></td><td></td><td>Bisphenol degradation</td><td>—</td><td>-0.43</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td></td><td></td><td>Toluene degradation</td><td>—</td><td>0.43</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr></tbody></table>

<footnote n="1">¹To determine a significant threshold to each performance-pathway pair, the matrix of the relative abundance of the pathways was randomized by pathway. A rCCA was then applied, and each correlation by pair was kept. This randomization followed by rCCA was repeated 1000 times. For each pair, the thresholds for a significant positive or negative correlation were the maximum or minimum correlation obtained during the 1000 randomizations, respectively. This table presents all and only the significant correlations.</footnote>

<footnote n="2">²The performances (ADFI, ADG and FE) were averaged during the deoxynivalenol (DON)/zearalenone (ZEN) challenge periods (P2 and P4) and corrected for the batch effect and the initial BW at 99 days of age using a linear model. The square mean values of ADFI, ADG and FE were used for the rCCA.</footnote>

<footnote n="3">³Pathways significantly correlated with performance (ADFI, ADG and FE) in control pigs but not in challenged pigs.</footnote>

<footnote n="4">⁴Pathways significantly correlated with performance (ADFI, ADG and FE) in challenged pigs but not in control pigs.</footnote>

<footnote n="5">⁵The treatment groups are: the CC group fed a control diet during the whole experiment, the DC group fed a DON- and ZEN-contaminated diet between 113 and 119 days, the CD group fed a DON- and ZEN-contaminated diet between 134 and 140 days, and the DD group fed a DON- and ZEN-contaminated diet between 113 and 119 days and between 134 and 140 days. In P2 (between 113 and 119 days of age), the control pigs belong to the CC and CD experimental groups and the challenged pigs belonged to the DC and DD groups. The performances were correlated with the microbial data at 119 days of age.</footnote>

<footnote n="6">⁶In P4 (between 134 and 140 days of age), the control pigs belong to the CC and DC experimental groups and the challenged pigs belonged to the CD and DD groups. The performances were correlated with the microbial data at 140 days of age.</footnote>
could be explained in part by pigs' rapid elimination of myco-
toxins. Half of the ingested DON is excreted in the urine
within 5.8 h (Dänicke et al., 2004), and all is excreted within
24 h (Prelusky et al., 1988).

The pathways extracted from the sPLS-DA differed from
those that were significantly correlated with the perfor-
mances. This lack of common pathways highlights limits of
a functional interpretation of the variables extracted from
the discriminant analysis. The information used to distinguish
groups does not always have biological significance or utility
and should be interpreted with caution. Nevertheless, on a
short-term basis, changes in relative OTU abundance due
to DON/ZEN exposure were due to a slight increase in biologi-
cal functions related to drug and xenobiotic metabolism by
cytochrome P450 (pathways used in the discriminant
analysis). This function restores gut health by degrading toxic
compounds such as xenobiotic mycotoxins (Sobrova et al.,
2010; Wilson and Nicholson, 2017). In the present study,
some Clostridium bacteria were impacted by the DON
challenge. In chickens, the Clostridiales can de-epoxify up
to 100% of DON into a much less toxic form: DOM-1 (Yu
et al., 2010). This ability of microbiota to detoxify DON
was previously observed in chickens (Yu et al., 2010;
Pierron et al., 2016), humans (Gratz et al., 2013) and pigs
(Eriksen et al., 2002; Gratz et al., 2018). In an in vitro study,
Chlebicz and Ślizewska (2019) identified Lactobacillus strains
that can detoxify DON, ZEN and other mycotoxins. These
results confirm that the microbiota seems to play a key role
in detoxifying mycotoxins. In livestock species, growing
evidence suggests a connection between microbiota compo-
sition and host performance. However, it is unclear whether a
change in the functioning of gut microbiota is the cause or
the consequence of variations in host performance.

As mentioned, microbiota could help hosts to maintain
good health when exposed to mycotoxins, with subsequent
positive effects on performance. Some microbial functions
were significantly correlated with ADFI only during challenge
periods. We hypothesize that mycotoxin challenge exacer-
bated differences between the pigs and emphasized signifi-
cant correlations between performance and the microbial
functions. It can be hypothesized that changes in animal per-
formance, especially ADFI, due to mycotoxin exposure would
also affect the composition of gut microbiota. In the present
study, one microbiota functional pathway related to xenobi-
otic biodegradation and metabolism and three pathways
related to the metabolism of terpenoids and polyketides
(from which the mycotoxins originated (Huffman et al.,
2010)) were negatively correlated with ADFI during the first
challenge period. This suggests that the pigs with ADFI
decreased by the DON/ZEN challenge had increased abun-
dance of these microbial functions. Specifically, some
changes in the functioning of gut microbiota were likely
related indirectly to the effect of DON/ZEN mycotoxins on
ADFI, rather than related directly to the effect of the myco-
toxins on microbiota composition.

Two pathways related to human diseases were signifi-
cantly correlated with ADFI of challenged pigs. Although

the meaning of these pathways in pigs is questionable, we
assume that these modifications were related to disturbance
in the health status of hosts that ingested DON/ZEN mycotox-
ins. Chronic ingestion of mycotoxins has impacted the health
of human and pig hosts: intestinal pathologies, disturbance
of intestinal barrier function and changes in the immune
response (Pinton and Oswald, 2014; Robert et al., 2017).
Effects are particularly strong for pigs: 4-week chronic expo-
sure of 8 mg DON/kg or 0.8 mg ZEN/kg in the feed changed
the immune response and resulted in a decrease in ADFI and
BW (Reddy et al., 2018b). The additional energy cost of the
inflammation response reduced the amount of energy avail-
able for growth.

Conclusion

Single or repeated short-term DON/ZEN challenges resulted
in transient modifications in the composition of fecal micro-
biota and microbial biological functions related to mycotoxin
detoxification. The challenges had no effect on long-term
changes in microbiota composition. Functions related to
metabolism and health status were correlated with ADFI
when pigs were exposed to a mycotoxin challenge, which
suggests that the microbiota helped their host manage the
mycotoxin disturbance.

Acknowledgements

The authors thank B. Duteil, P. Roger, N. Muller, J. Delamarre and H. Demay from the UEPR for animal care and sample collection and C. Barilly from UMR GenPhySE, INRAE, 31320 Castanet-Tolosan, France for lab analysis.

This study is part of the Feed-a-Gene Project funded by the
European Union's H2020 Program (grant 633531).

D. Renaureau 0000-0002-9306-2109

Declaration of interest

The authors declare that there is no conflict of interest.

Ethics statement

The present study was conducted in accordance with the French legislation on animal experimentation and ethics (authorization: 2016022415253973).

Software and data repository resources

None of the data were deposited in an official repository.

Supplementary material

To view supplementary material for this article, please visit
https://doi.org/10.1017/S1751731120001299

References

Andretta I, Kipper M, Lehnen CR, Hauschild L, Vale MM and Lovatto PA 2012. Meta-analytical study of productive and nutritional interactions of mycotoxins in growing pigs. Animal 6, 1476–1482.
Impact of mycotoxin challenge on pig microbiota

Antonissen G, Martel A, Pasmans F, Ducatelle R, Verbrugghe E, Vandenbroucke V, Li S, Haesebrouck F, Van Immerseel F and Croubels S 2014. The impact of Fusarium mycotoxins on human and animal host susceptibility to infectious diseases. Toxins 6, 430–452.

Binder SB, Schwartz-Zimmermann H, Varga E, Bichl G, Michlmayr H, Adam G and Berthiller F 2017. Metabolism of zearalenone and its major modified forms in pigs. Toxins 9, 56.

Bouhet S and Oswald IP 2005. The effects of mycotoxins, fungal food contaminants, on the intestinal epithelial cell-derived innate immune response. Veterinary Immunology and Immunopathology 108, 199–209.

Chlebicz A and Ślizewska K 2019. In vitro detoxification of Aflatoxin B1, Deoxynivalenol, Fumonisins, T-2 toxin and Zearalenone by probiotic bacteria from genus *Lactobacillus* and *Saccharomyces cerevisiae* yeast. Probiotics and Antimicrobial Proteins 12, 289–301.

Dänicke S, Valenta H and Döll S 2004. On the toxicokinetics and the metabolism of deoxynivalenol (don) in the pig. Archives of Animal Nutrition 58, 169–180.

Eriksen GS, Pettersson H, Johnsen K and Lindberg JE 2002. Transformation of trichothecenes in ileal digesta and faeces from pigs. Archiv für Tierernaehrung 56, 263–274.

Eriksen GS, Pettersson H and Lindberg JE 2003. Absorption, metabolism and excretion of 3-acetyl DON in pigs. Archiv für Tierernaehrung 57, 335–345.

Flannery BM, Wu W and Pestka JJ 2011. Characterization of deoxynivalenol-induced anorexia using mouse bioassay. Food and Chemical Toxicology 49, 1863–1869.

Frese SA, Parker K, Calvert CC and Mills DA 2015. Diet shapes the gut microbiome of pigs during nursing and weaning. Microbiome 3, 28.

González I, Déjean S, Martin PG and Baccini A 2008. CCA: an R package to extend canonical correlation analysis. Journal of Statistical Software 23, 1–14.

Gratz SW, Currie V, Richardson AJ, Duncan G, Holtrop G, Farquharson F, Louis P, Pinton P and Oswald IP 2018. Porcine small and large intestinal microbiota rapidly hydrolyze the masked mycotoxin deoxynivalenol-3-glucoside and release deoxynivalenol in spiked batch cultures in vitro. Applied and Environmental Microbiology 84, e02106–17.

Gratz SW, Duncan G and Richardson AJ 2013. The human fecal microbiota metabolizes deoxynivalenol and deoxynivalenol-3-glucoside and may be responsible for urinary deepoxy-deoxynivalenol. Applied and Environmental Microbiology 79, 1821–1825.

He P 1992. Microbial transformation of deoxynivalenol (vomitoxin). Applied and Environmental Microbiology 58, 3857–3863.

Huffman J, Gerber R and Du L 2010. Recent advancements in the biosynthetic mechanisms for polyketide-derived mycotoxins. Biopolymers 93, 764–776.

Kollarczik B, Gareis M and Hanelt M 1994. In vitro transformation of the fusarium mycotoxins deoxynivalenol and zearalenone by the normal gut microflora of pigs. Natural Toxins 2, 105–110.

Lê Cao K-A, Boitard S and Besse P 2011. Sparse PLS discriminant analysis: biologically relevant feature selection and graphical displays for multiclass problems. BMC Bioinformatics 12, 253.

Mateos I, Combes S, Pascal G, Cauquil L, Barilly C, Cossalter A-M, Laffitte J, Botti S, Pinton P and Oswald IP 2018. Fumonisin-exposure impairs age-related ecological succession of bacterial species in weaned pig gut microbiota. Toxins 10, 230.

McMurdie PJ and Holmes S 2013. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE 8, e61217.

Nagl V, Woechtl B, Schwartz-Zimmermann HE, Hennig-Pauka I, Moll W-D, Adam G and Berthiller F 2014. Metabolism of the masked mycotoxin deoxynivalenol-3-glucoside in pigs. Toxicology Letters 229, 190–197.

Payros D, Dobrindt U, Martin P, Secher T, Bracarense APFL, Boury M, Laffitte J, Pinton P, Oswald E and Oswald IP 2017. The food contaminant deoxynivalenol exacerbates the genotoxicity of gut microbiota. mBio 8, e00007–17.

Pierron A, Mimoun S, Murate LS, Loiseau N, Lippi Y, Bracarense A-PFL, Schatzmayr G, He JW, Zhou T, Moll W-D and Oswald IP 2016. Microbial biotransformation of DON: molecular basis for reduced toxicity. Scientific Reports 6, 29105.

Pinton P and Oswald IP 2014. Effect of deoxynivalenol and other type B trichothecenes on the intestine: a review. Toxins 6, 1615–1643.

Pomar J, López V and Pomar C 2011. Agent-based simulation framework for virtual prototyping of advanced livestock precision feeding systems. Computers and Electronics in Agriculture 78, 88–97.

Prelusky DB, Hartin KE, Trenholm HL and Miller JD 1988. Pharmacokinetic fate of 14C-labeled deoxynivalenol in swine. Fundamental and Applied Toxicology 10, 276–286.

R Core Team 2018. R: a language and environment for statistical computing. R foundation for statistical computing, Vienna, Austria.

Reddy K, Song J, Lee H-J, Kim M, Kim D-W, Jung H, Kim B, Lee Y, Yu D, Kim D-W, Oh Y and Lee S 2018b. Effects of high levels of deoxynivalenol and zearalenone on growth performance, and hematological and immunological parameters in pigs. Toxins 10, 114.

Reddy KE, Jeong JY, Song J, Lee Y, Lee H-J, Kim D-W, Jung HJ, Kim KH, Kim M, Oh YK, Lee SD and Kim M 2018a. Colon microbiome of pigs fed diet contaminated with commercial purified deoxynivalenol and zearalenone. Toxins 10, 347.

Robert H, Payros D, Pinton P, Théodorou V, Mercier-Bonin M and Oswald IP 2017. Impact of mycotoxins on the intestine: are mucus and microbiota new targets? Journal of Toxicology and Environmental Health, Part B 20, 249–275.

Robinson MD, McCarthy DJ and Smyth GK 2010. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140.

Rognes T, Flouri T, Nichols B, Quince C and Mahé F 2016. VSEARCH: a versatile open source tool for metagenomics. PeerJ 4, e2584.

Saint-Cyr MJ, Perrin-Guyomard A, Houée P, Rolland J-G and Laurentie M 2013. Evaluation of an oral subchronic exposure of deoxynivalenol on the composition of human gut microbiota in a model of human microbiota-associated rats. PLoS ONE 8, e80578.

Serviento AM, Brossard L and Renaudeau D 2018. An acute challenge with a deoxynivalenol-contaminated diet has short- and long-term effects on performance and feeding behavior in finishing pigs. Journal of Animal Science 96, 5209–5221.

Sobrova P, Adam V, Vasatkova A, Beklova M, Zeman L and Kizek R 2010. Deoxynivalenol and its toxicity. Interdisciplinary Toxicology 3, 94–99.

Waché YJ, Valat C, Postollec G, Bougeard S, Burel C, Oswald IP and Fravalo P 2008. Impact of deoxynivalenol on the intestinal microflora of pigs. International Journal of Molecular Sciences 10, 1–17.

Wilson ID and Nicholson JK 2017. Gut microbiome interactions with drug metabolism, efficacy and toxicity. Translational Research 179, 204–222.

Young JC, Zhou T, Yu H, Zhu H and Gong J 2007. Degradation of trichothecene mycotoxins by chicken intestinal microbes. Food and Chemical Toxicology 45, 136–143.

Yu H, Zhou T, Gong J, Young C, Su X, Li X-Z, Zhu H, Tsao R and Yang R 2010. Isolation of deoxynivalenol-transforming bacteria from the chicken intestines using the approach of PCR-DGGE guided microbial selection. BMC Microbiology 10, 182.