2.1 Animal Selection Criteria
The trial was performed at the Experimental Animal Research and Application Centre of University of Milan and was authorized by the Italian Health Ministry (authorization n° 711/-PR) in accordance with EU regulations (Directive 2010/63/EU).
Animals enrolled in the experimental trial were selected from a conventional herd free from contagious diseases (Ex A-list International Office of Epizootic, porcine reproductive and respiratory syndrome, atrophic rhinitis, Aujeszky’s disease, transmissible gastroenteritis, salmonellosis) and without a history of PWD or oedema disease. Sows were assessed for genetic susceptibility to Escherichia coli carrying F18 adhesive fimbriae (F18 E. coli) by screening the fucosyltransferase 1 (FUT1) genotypes using polymerase chain reaction (PCR) reaction according to Luise et al. (2019a, b). Briefly, genomic DNA was extracted from hair samples of sows and genotyped to identify polymorphic variants. Sows carrying the GG genotypes at FUT1 gene were considered for piglet enrolment. A further selection criterion was the absence of hemolytic E. coli in piglets feces. Microbiological analyses of selective mediums (Agar MacConkey) (Hayer et al. 2020; Li et al. 2020; Remfry et al. 2020) were thus carried out before transport and upon arrival on fecal samples collected from enrolled piglets.
2.2 Animals and Experimental Design
Twenty-seven weaned piglets (28±2 days) balanced per weight (9.79±1.25 kg) and sex, after a 7-day adaptation period, were randomly allotted in four experimental groups in randomized complete block design and fed ad libitum for the entire experimental period: phytogenic additive 1 (PHY1, n=7), phytogenic additive 2 (PHY2, n=7), positive control (CTRL+, n=6), and negative control (CTRL -, n=7). After the adaptation period, the basal diet was added by 200g/100kg phytogenic additive (PHY1) and by phytogenic plus 2000 ppm of short and medium chain fatty acids premix (PHY2). Control groups received basal diet supplemented with a premix carrier (95% wheat meal and 5% of coconut oil) to achieve the same nutrient concentrations. The iso-energetic and iso-proteic diets (Table S1) were formulated (Plurimix; Fabermatica, CR, Italy) according to animal requirements for the post weaning phase defined by the US National Research Council (National Research Council 2012). The phytogenic feed additive (FRESTA®F, Delacon Biotechnik GmBH), approved by EU regulation (Reg. CE 1831/2003), as zootechnical additive, was composed of essential oil from caraway oil (d-carvone 3.5-6.0 mg/g) and lemon (limonene: 2.3 - 9.0 mg/g), dried herbs and spices (1.5% clove powder, 10% cinnamon powder, 1.5% nutmeg powder, 5% onion powder, 2% pimento powder, 5% orange peel powder, 12.5% peppermint powder and 12.5% chamomile powder). The SCFA and MCFA premix was composed by butyric (C4), caprylic (C8), capric (C10) and lauric acid (C12). The phytogenic products (with or without SCFA and MCFA) or premix carrier were mixed with the basal diet for 30 minutes.
Piglets were housed in two environmentally controlled rooms, in individual pens, with a plastic slatted floor and constant temperature (27° C) and humidity (60%) for the entire experimental period. The trial was divided into a pre- and post-challenge, considering the challenge as day 0 (Figure 1).
PHY1: treatment group fed basal diet supplemented with 200g/100kg of phytogenic additive; PHY2: treatment group fed basal diet supplemented with 200g/100kg of phytogenic additive supplemented with 2000 ppm of short and medium chain fatty acids premix; phytogenic additive composed by caraway oil, lemon oil, clove, cinnamon, nutmeg, onion, pimento, orange peel, peppermint and chamomile powder; short chain and medium chain fatty acids blend included butyric (C4), caprylic (C8), capric (C10) and lauric acid (C12) at dose levels of 2000 ppm. CTRL: group fed basal diet supplemented with premix carrier divided into negative control (CTRL-) and positive control (CTRL+) challenged at day 0. Diets were analyzed for proximate analysis, including moisture, crude protein (CP), crude fibre (CF), ether extract (EE), and ash. The moisture determination was performed by oven-drying at 65°C for 24 h (Regulament EC 152/2009). Crude protein content was measured according to the Kjeldahl method (Association of Official Analytical Chemists method 2001.11). Crude fiber was determined by the filter bags technique (American Oil Chemistry Society 2009). Ether extract content was determined in a Soxhlet system after hydrolysis (Association of Official Analytical Chemists method 2003.05). Ash was measured using a muffle furnace at 550°C (Association of Official Analytical Chemists method 942.05).
2.3 Experimental Challenge
E. coli challenger strain was genetically characterized by polymerase chain reaction (PCR) (Applied Biosystem 7500) in order to detect the presence of the two important virulence profile: subunit B of verocytotoxin type 2 and F18 adhesive fimbriae (Table 1).
Twenty piglets, except for piglets in CTRL- group (n = 7), on day 0 (challenge day) were orally infected with O138 Escherichia coli F18+ strain obtained from a permanent collection of the University of Milan and previously characterized (Rossi et al. 2014b; Dell'Anno et al. 2020; Rossi et al. 2021).
Sixty minutes before the challenge, the piglets were sedated with azaperone (Stresnil TM, Janssen Cilag Spa, 2 mL/pig), thereafter 30 mL of a 10% bicarbonate solution was orally administered to neutralize gastric acid and to increase the survival rate of the challenger strain in the stomach. After 10-15 min, the inoculum was given orally in a single dose of 5 mL of bacterial medium with 2x109 colony-forming unit (CFU) of challenger strain, using a 16G catheter (Rossi et al. 2021). Animals were fasted 3h before and 3h after the challenge. At the same time, piglets in CTRL- were orally inoculated with 5 mL of Luria Bertani (LB) medium to balance the level of stress associated with the oral challenge.
Table 1
PCR conditions and oligonucleotide sequences of F18 adhesive fimbriae and VTe2 (B-subunit) encoding genes.
Gene | Accession number (GenBank) | Size (pb) | Primer sequence (5’ to 3’) | PCR conditions |
F18 adhesive fimbriae | AJ308332.1 | 519 | 5’GATCCATGAAAAGACTAGTGTTTATTTCTTTTG 3’CGAATGCGCCAATGAATGTTCATTCTCGAG | Den.95°C 1’ ann.56°C 1’20’’ ext.72°C 1’30’’ 35 cycles |
VTe2 (B-subunit) | GU459254.1 | 270 | 5’GGATCCATGAAGAAGATGTTTATAGCGG 3’AACGGGTCCACTTCAAATGATTCTCGAG | Den.95°C 1’ ann.50°C 1’20’’ ext.72°C 1’30’’ 35 cycles |
2.4 Zootechnical Performance, Clinical and Fecal Score
Average daily feed intake (ADFI) was recorded daily from day -6 to day 7 by measuring the refusals. Body weight (BW) was recorded on day -6 (first day of experimental diets), day 0 (challenge day), day 4 and day 7 (sacrifice day). Average daily gain (ADG) and feed efficiency were also calculated.
Piglets were individually evaluated throughout the trial by clinical examination, including observation of behavioral disturbances. In particular, oedema, epiphora, respiratory and hair scores were evaluated through three-point scales (oedema score: 0=normal, 1=mild, 2=severe; epiphora score: 0=normal, 1=mild, 2=severe; respiratory score: 0=normal, 1=slightly quick, 2=quick; hair/bristles score: 0=smooth, 1=lightly brushy, 2=highly brushy) (Rossi et al. 2021). In addition, cyanosis, a blue or red discoloration of the skin, which may or may not be localized to small areas, was considered not as a specific skin condition but as a symptom of disease. From day -6 to day 7, all piglets were evaluated for the fecal score. Clinical signs of the disease were identified according to the point scale score described by Rossi et al. (Rossi et al. 2014b). A four-point scale was adopted to score fecal consistency: 0=normal, 1=soft consistency, 2=mild diarrhea, 3=severe diarrhea; considering >1 as an indicative of diarrhea. Fecal color was evaluated using a three-point scale: 1 = yellow, 2 = green; 3 = brown.
2.5 Microbiological Evaluation of Fecal Samples
Individual fecal samples were collected from rectal ampulla from each piglet, on days -1, 1, 2, 3 and 4 to perform microbiological analysis and evaluate the challenger strain shedding. For each sample, 1 g of feces was homogenized with 1 ml of saline solution and incubated overnight at 37°C on sheep blood agar plates 5% (Blood Agar Base No. 2-Oxoid) in order to examine the presence of hemolytic colonies. The total hemolytic bacteria count was performed by counting the number of colonies cultured from serial dilutions of each fecal sample in order to evaluate the presence of hemolytic E. coli in relation to the total bacteria population.
2.6 Necropsy, Intestinal Samples, and Histopathology
At day 7 post-challenge, sixteen animals (n=4/treatment) were sacrificed, tissue samples were collected for histopathological and molecular analyses of intestinal tissues. The intestine of each animal was weighed, and intestinal samples of ileum were harvested. For the histological evaluation, intestinal samples were fixed in 10% neutral formalin buffer and stored at 4°C. Tissues were rinsed with sterile saline solution and transferred into 2 mL cryotubes, snap-frozen in liquid nitrogen and stored at -80°C until further analysis.
Histological examinations of collected intestinal samples for each piglet were carried out. The fixed samples fixed were embedded in paraffin, and 5µm thick histological sections were performed with a microtome. Cross sections were stained with hematoxylin and eosin and were blind evaluated by light microscopy. A four-point scale was adopted for inflammatory infiltrates, epithelial regeneration, fusion of villi, oedema, hyperemia, necrosis of mucosa, T atrophy, stroma, and follicular hyperplasia; considering: 0=no evidence; 1=slight presence; 2=moderate; 3=severe. Samples of duodenum were collected and frozen in liquid nitrogen for gene expression analysis.
2.7 Duodenum Gene Expression
Total RNA was extracted from the duodenum using FastGene Scriptase Basic (Nippon genetics) according to the manufacturer’s instructions. The integrity of total RNA was assessed by gel electrophoresis to detect the 18S and 28S rRNA bands. A combination of oligo-dT and random primers was used to reverse transcribe 100 ng of total duodenal RNA to cDNA (cDNA synthesis kit, FastGene Scriptase Basic, Nippon Genetics). Primer pairs were first tested for their specificity in qualitative PCR, using the pooled cDNA as a template. The cycling profile for the assay consisted of initial denaturation of RNA (65°C x 5’), then the annealing of random primers (25°C x10’), followed by the annealing of oligo-dT and transcription (42°C x 60’). At the end of the cycle, the enzyme deactivation (90°C x 5’) was performed. The abundance of cytochrome c oxidase subunit I (COX1), cytochrome c oxidase subunit II (COX2), interleukin 10 (IL-10), interleukin 6 (IL-6), lysyl oxidase (LOX), glutathione peroxidase 2 (GPX2), NAD (P) H quinone dehydrogenase 1 (NQ01) claudin domain containing 1 (CLDND1) and occludin (OCLN) (Table 2) mRNA was determined using SYBR Green-based real-time quantitative PCR assays (7500 Fast Dx, Applied Biosystems). Only reaction efficiencies that were near to 100% were considered for further analysis. The mean values for the transcripts were normalized to the arithmetic mean of mRNA abundance of βactin as the reference gene within each sample. The comparative CT method was used to determine fold changes in gene expression, calculated as 2-∆∆CT. The final results were presented as the fold changes of target gene expression in a target sample relative to a reference sample, normalized to βactin rRNA (Livak and Schmittgen, 2001). The βactin rRNA was used to calculate the threshold cycles, since it previously showed constant values under all the conditions adopted.
Table 2
Primer sequences and relative amplicon dimensions.
Gene1 | Accession number (GenBank) | Size (pb) | Primer sequence (5’ to 3’) |
βactin F | DQ845171 | 76 bp | CTACGTCGCCCTGGACTTC |
βactin R | DQ845172 | | GCAGCTCGTAGCTCTTCTCC |
IL-6 F | JQ839263 | 112 bp | TGGGTTCAATCAGGAGACCT |
IL-6 R | JQ839264 | | CAGCCTCGACATTTCCCTTA |
IL-10 F | L2001 | 105 bp | TGAAGAGTGCCTTTAGCAAGCTC |
IL-10 R | L2002 | | CTCATCTTCATCGTCATGTAGGC |
COX1 F | EF568726 | 102 bp | GGAGCGGGTACTGGATGAAC |
COX1 R | EF568726 | | CACCTGCAAGGGTGTAGGGAGL |
COX2 F | AF304201 | 141 bp | AAGACGCCACTTCACCCATC |
COX2 R | AF304201 | | TCCATTGTGCTAGTGTGTGTCA |
GPx2 F | DQ898282 | 103 bp | GGAGATCCTGAACAGCCTCA |
GPx2 R | DQ898282 | | GCGAAGACAGGATGCTCATT |
LOX F | NM_001164001 | 112 bp | GTGGAGCACGAAAGCAAGACCC |
LOX R | NM_001164001 | | AAGGTGGGGTATGCATCGACAC |
NQ01 F | NM_001159613 | 118 bp | ATCACAGGTAAACTGAAGGACCC |
NQ01 R | NM_001159613 | | GCGGCTTCCACCTTCTTTTG |
CLAUDIN1 F | NM_001244539 | 90 bp | TCTTTCTTATTTCAGGTCTGGCT |
CLAUDIN1 R | NM_001244539 | | ACTGGGGTCATGGGGTCATA |
OCCLUDIN F | NM_001163647 | 106 bp | GTCCACCTCCTTATAGGCCTGATG |
OCCLUDIN R | NM_001163647 | | CGCTGGCTGAGAAAGCATTGG |
1CTB: actin beta; IL-6: interleukin-6; IL-10: interleukin 10; COX1: cytochrome c oxidase subunit I; COX2; cytochrome c oxidase subunit II; LOX: lysyl oxidase; GPX2: glutathione peroxidase 2; NQ01: NAD (P) H quinone dehydrogenase 1; CLDND1: claudin domain containing 1; OCLN: occludin. |
2.8 Blood Samples, Serum Metabolite Profile and Serum Acute Phase Proteins
Blood was collected from the jugular vein of each animal on day -1, day 3 and day 7 through vacutainer tubes containing ethylenediaminetetraacetic acid (EDTA) and heparin as anticoagulants. Plasma was collected after centrifugation (3000 rpm, 10 min, 4°C), aliquoted and stored at -20°C for further analysis. Hematocrit was evaluated on whole blood using the microhematocrit method. The concentration of total protein (g/L), albumin (g/L), globulin (g/L), albumin/globulin (A/G ratio), alanine aminotransferase (ALT-GPT; IU/L), aspartate aminotransferase (AST-GOT; IU/L), phosphatase alkaline (ALP; IU/L), glucose (mmol/L), urea (mmol/L), creatinine (µmol/L), total bilirubin (µmol/l), total cholesterol (mmol/L), triglycerides (mmol/L), high-density lipoprotein (HDL; mmol/L), low-density lipoproteins (LDL; mmol/L), calcium (mmol/L), phosphorus (mmol/L) and magnesium (mmol/L) were analyzed via standard enzymatic colorimetric analysis through a multiparametric autoanalyzer for clinical chemistry (ILab 650; Instrumentation Laboratory Company, Lexington, MA, USA) at 37°C by the Lombardy and Emilia Romagna Experimental Zootechnic Institute (IZSLER). Porcine C-reactive protein (CRP) concentration was determined with a commercial sandwich immunoassay Kit (Mybiosource, San Diego, CA, USA) following the manufacturer’s instructions. The results were read at 450 nm using a microplate reader (Model 680, Bio-Rad Laboratories, CA, USA). Haptoglobin (HP) concentrations were measured through a colorimetric kit (PhaseTM Range porcine Haptoglobin Assay; Tridelta Development Ltd) according to the manufacturer’s instructions. The results were read at 630 nm on a microplate reader (Model 680, Bio-Rad Laboratories, CA, USA).
2.9 Statistical Analysis
Zootechnical performance and fecal microbiological analysis were analyzed using a linear model after testing the normality of data through Shapiro-Wilk test using JMP Pro 15® (SAS Inst. Inc., Cary, NC, USA). The model included the fixed effect of treatments (Trt), the effect of time (Time), and the interaction between treatment and time (Trt x Time).
Serum metabolites were evaluated performing analysis of covariance (ANCOVA) to adjust the initial variability of the pre-challenge period after testing the normality of data through Shapiro-Wilk test using JMP Pro 15® (SAS Inst. Inc., Cary, NC, USA).
Clinical score data were converted into a dichotomous variable (normal/pathological), and observed frequencies were assessed using the Chi-squared Test. Histological scores, intestinal weight and relative gene expression were analyzed using Kruskal-Wallis test (PROC NPAR1WAY of SAS 9.4 software) for non-parametric data due to the small sample size of euthanized animals at day 7. Multiple comparisons for parametric statistics were evaluated with the Tukey’s Honestly Significant Difference test (Tukey’s HSD) or Tukey-Kramer test and Steel-Dwass test was used for non-parametric multiple comparisons. The results were presented as least square means (LSMEANS) ± standard error (SE) for parametric data and as medians and range (minimum-maximum) for non-parametric results. Means or medians were considered statistically different when p≤0.050 and statistical tendency was considered when p<0.100.