Animal ethics
All the programs performed in our animal experiment were endorsed and authorized by the Institutional Animal Care and Use Committee of China Agricultural University (No.AW10601202-1-2, Beijing, China).
MOA product and antibiotic
The product of MOA combination named PORCINATTM was supplied by Jefo (Jefagro, Canada), which is a selected formulation of essential oils primarily containing thymol, vanillin and eugenol and organic acid mainly containing fumaric, citric, butyric and sorbic acid microencapsulated in the triglyceride matrix of hydrogenated vegetable oils. The chlortetracycline was sourced from Tongli Xingke (Beijing Tonglixingke, China).
Experimental design and diets
A total of 120 piglets (Duroc × [Landrace × Yorkshire]; weighted 7.66 ± 1.79 kg, weaned at d 28) were selected and randomly allocated to 3 treatments with 4 replicates per group and 10 piglets per replicate according to the initial body weight and gender. The dietary treatment was as follows: 1) corn-soybean meal basal diets (Ctrl); 2) Ctrl + 75 mg/kg chlortetracycline (AGP); 3) Ctrl + 1,500 mg/kg MOA (MOA). The experiment period was lasted for 21 d. Table 1 lists the composition and nutritional levels of basal diets, which satisfied or excelled the NRC (2012) requirements.
Table 1
Composition and nutrient profile of the basal diets (as fed basis, %)
Ingredients | Content | Nutritional level | Content |
Corn, CP 7.6% | 60.24 | Calculated values | |
Whey powder, CP 3.8% | 10.00 | Digestive energy, kcal/kg | 3,487 |
Soybean meal, CP 43.6% | 13.90 | Crude protein | 18.50 |
Fermented soybean meal, CP 51% | 5.00 | Ether extract | 4.20 |
Angel yeast, CP 51.3% | 3.00 | Lactose | 8.00 |
Fish meal, CP 65.3% | 3.00 | Calcium | 0.60 |
Soybean oil | 1.00 | Phosphorus | 0.48 |
Salt | 0.40 | Sodium | 0.31 |
Dicalcium phosphate | 0.54 | Lysine | 1.54 |
Limestone | 0.28 | Methionine | 0.56 |
L-Lysine HCl, 78% | 0.71 | Threonine | 0.83 |
DL-Methionine, 98% | 0.33 | Tryptophan | 0.27 |
L-Threonine, 98% | 0.33 | |
L-Tryptophan, 98% | 0.11 | Analyzed values 2 | |
Valine | 0.18 | Gross energy, kcal/kg | 3,867 |
Zinc oxide | 0.20 | Crude protein | 18.42 |
10,000-IU Phytase | 0.03 | Ether extract | 4.37 |
Chromic oxide | 0.25 | Dry matter | 87.29 |
Non-antibiotic premix1 | 0.50 | Organic matter | 94.80 |
Total | 100.00 | | |
Control (Ctrl) = a corn soybean meal-based diet. AGP = Ctrl + 75 mg/kg chlortetracycline. MOA = Ctrl + 1,500 mg/kg MOA. |
1 Non-antibiotic premix for per kilogram diet included: vitamin A, 12,000 IU; vitamin D3, 2,000 IU; vitamin E, 24 IU; vitamin K3, 2.0 mg; vitamin B1, 2.0 mg; riboflavin, 6.0 mg; vitamin B6, 3 mg; vitamin B12, 24 µg; nicotinic acid, 30 mg; pantothenic acid, 20 mg; folic acid, 3.6 mg; biotin, 0.1 mg; choline chloride, 0.4 mg; iron, 96 mg; copper, 8.0 mg; zinc, 120 mg; manganese, 40 mg; iodine, 0.56 mg; selenium, 0.4 mg. |
2 The analyzed values were the average of the 3 feed nutrient levels measured in the Ctrl antibiotic growth promoter (AGP) and microencapsulated essential oils and organic acids (MOA) groups. |
The nursery house was thoroughly cleaned and disinfected in advance of the experiment, while the temperature (28ºC to 30 ºC for the first week, then lowered by 1 to 2ºC weekly up to 23ºC to 25ºC), humidity (60–70%) and CO2 concentration (below 0.15%) in the nursery house were controlled by an automatic monitoring system. The piglets were fed and drunk ad libitum overall. Also, the piglets were vaccinated and dewormed routinely following the standard procedures of pig farms, the pigsty were periodically cleaned, the feeding and health status of the weaner were observed and recorded anytime.
Sampling collection and detection method
On d 11 and d 21 of the experiment, 10 mL of blood was gathered from the anterior vena cava of piglets into a vacutainer, centrifuged at 3,000 × g for 15 min, the serum was separated and stored at -20℃ in 0.5 mL centrifuge tubes for analysis of serum immune function and antioxidant characteristic.
On the d 18 of the experiment, feces of piglets in the nursery were cleaned, and fecal samples in each replicate were gathered from the d 19 to 21, the sample collection twice per day without contamination, then stored at -20ºC for determination of the apparent digestibility of nutrients. Moreover, approximately 1 kg of representative feed samples were harvested weekly during the experiment.
On d 21 of the experiment, one piglet with mean weight was sampled for slaughter per replicate, the liver was gathered and approximately 1 to 2 cm of intestinal samples were collected from the duodenum, jejunum and ileum at the 1/3 of the posterior segment respectively, the intestinal contents were washed off with 0.9% saline gently, placed in 10 mL cryovials and pre-stored in a liquid nitrogen tank, then transferred to -80°C for determination of antioxidant enzyme properties, digestive enzyme activity and the expression of intestinal tight junction protein gene. Meanwhile, the tissues of approximately 2 cm small intestine were gathered, washed with sterile saline and fixed in 4% paraformaldehyde for determining intestinal morphology. Furthermore, the cecum and colon contents of piglets were collected for volatile fatty acid analysis and 16srRNA gene sequencing. The surgical trays, scalpel scissors and other instruments used in the sample collection process and the operating table were disinfected with 75% alcohol.
Growth performance measurement
Piglets were weighed on d 11 and d 21 of the experiment as well as recorded the feed consumption to calculate average daily gain (ADG), average daily feed intake (ADFI) and feed conversion ratio (FCR = ADFI/ADG). The piglets' anuses were checked one by one at 09:00 and 17:00 daily during the experiment to observe and recorded any fecal contamination and redness. The number of piglets with diarrhea per treatment was counted at the end of the experiment and the diarrhea rate was calculated with the following formulation:
Diarrhea rate (%) = 100% × total number of piglets with diarrhea/(total number of piglets × number of days)
Serum biochemical immunity and intestinal enzyme activity
The serum was defrosted at 4°C and mixed well before analysis. The enzyme activities of malondialdehyde (MDA), superoxide dismutase (SOD), glutathione peroxidase (GSH-Px) and total antioxidant capacity (T-AOC) in serum were analyzed by automatic biochemical analyzer (GF-D200, Gaomi Rainbow Analytical Instrument, China). Immunoglobulins (IgA, IgG, IgM), interleukin-1β (IL-1β), interleukin-10 (IL-10), gamma-interferon (IFN-γ) tumor necrosis factor-α (TNF-α), D-lactic acid (DLA) and diamine oxidase (DAO) levels in serum were tested by ELISA with Multiskan Ascent enzyme marker (Thermo Scientific, USA). The activities of small intestinal amylase, lipase, trypsin and chymotrypsin were measured by Immunoturbidimetry. The commercial kits were sourced from Nanjing Jiancheng Institute of Biological Engineering (Nanjing, China).
RNA extraction and real-time PCR
The total RNA was extracted from small intestine of piglets by means of Trizol method according to the manufacture instruction, the concentration and quality of RNA were detected by protein nucleic acid assay (ND-2000UV, Thermo Fisher, USA) and 1% agarose gel electrophoresis. The RNA was reverse-transcribed into cDNA using the TransScript All-in-One First-Strand cDNA Synthesis SuperMIX for qPCR kit (QIAGEN, Germany). The reverse transcription system: total RNA, 0.5 µg; 5×TransScript All-in-one SuperMix for qPCR, 5 µL; gDNA Remover, 0.5 µL; Nuclease-free H2O was added to 10 µL. The reaction procedure: 42 ℃ for 15 min, 85℃ for 5 s. Addition of 90 µL Nuclease-free H2O after reverse transcription, then held at -20℃.
Real-time PCR was conducted by a LightCycler® 480 Ⅱ Real-time PCR Instrument (Roche, Swiss) (PCR efficiency: 94–105%) with 10 µL of PCR reaction mixture, which included 1 µL of cDNA, 5 µL of 2× PerfectStartTM Green qPCR SuperMix, 0.2 µL of forward primer, 0.2 µL of reverse primer and 3.6 µL of nuclease-free water. The reactions were incubated in 384-well optical plates (Roche, Switzerland) for 30 s at 94°C, followed by 45 cycles of 5 s at 94°C, 30 s at 60°C. The melting curve analysis was performed at the end of the PCR cycle to verify the specific generation of the expected PCR product. Triplicate analyses were performed for each sample. The primer sequences shown in Table 2 were designed by Ouyi Biotech (Shanghai, China) and synthesized by TsingKe Biotech (Beijing, China) based on the mRNA sequences available from the NCBI database. The expression levels of mRNAs were normalized to the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and were calculated using the 2−ΔΔCt method.
Table 2
Primer sequences of housekeeping and target genes concerned with intestinal barrier function
Item1 | Sequences (5' to 3')2 | Length, bp | Tm, ℃ |
Occludin | F: GTGGGACAAGGAACGTATT R: TCTCTCCGCATAGTCCGAA | 115 | 60 |
Claudin-1 | F: ATACAGGAGGGAAGCCAT R: ATATATTTAAGGACCGCCCTCT | 89 | 60 |
ZO-1 | F: GCTCAGCCCTATCCATCT R: GGACGGGACCTGCTCATAA | 90 | 60 |
Mucin-1 | F: GTGCCGCTGCCCACAACCTG R: AGCCGGGTACCCCAGACCCA | 141 | 60 |
Mucin-2 | F: CAGACCTACTCAGAGTTCCT R: CTCGGGCTTGTTGATCTT | 84 | 60 |
GAPDH | F: CAGCAATGCCTCCTGTACCA | 72 | 60 |
R: ACGATGCCGAAGTTGTCATG |
1 ZO-1 = zonula occludens-1; GAPDH = glyceraldehyde-3-phosphate dehydrogenase. |
2 F, forward primer; R, Reverse primer. |
Apparent total tract digestibility (ATTD) of nutrients
The fecal samples were thawed at 4℃, then dried at 65℃ for 72 h, the samples of feed and feces were ground through a 40-mesh (425 micrometers) screen before analysis. Dry matter (DM), crude protein (CP), ether extract (EE) and ash were analyzed in accordance with the Association of Official Analytical Chemists [19]; Neutral detergent fiber (NDF) and Acid detergent fiber (ADF) were determined with reference to the method of Vansoest et al. [20] (A2000i fiber analyzer, Ankom, USA). The gross energy in feed and feces samples was analyzed using an oxygen bomb calorimeter (PARR 6400, Moline, USA). The chromium levels in feed and feces were analyzed by atomic absorption spectrophotometer (Z-5000; Hitachi, Japan) based on the methodology of Williams et al. [21] for calculating the ATTD with the following equation.
Apparent total tract digestibility (ATTD, %) = 1 - (Cr feed × Nutrient feces) / (Cr feces × Nutrient feed)
Intestinal morphology
After the small intestine samples were fixed in 4% paraformaldehyde solution for 48 h, the samples were rinsed, excised, and dehydrated with ethanol for 24 h, then paraffin-embedded, sliced in 4 cross-sections and stained with hematoxylin-eosin. Lastly, the samples were morphologically detected by light microscopy (Olympus CX31, Tokyo, Japan), and 10 intact, well-oriented villi-crypt units were measured in each section. The villi height (VH) was evaluated from the top of the villi to the junction of the villi and crypt, the crypt depth (CD) was determined as the depth of the villi invagination, then calculated the ratio of villi height to crypt depth (VH:CD).
Cecal volatile fatty acids
The cecal contents of piglets were thawed at 4ºC and mixed, approximately 0.5 g of the sample was weighed into a 10-mL centrifuge tube, adding 8 mL of ultrapure water, sonicated in an ice water bath for 30 min (mixing every 10 min), then centrifuged at 15,000 × g for 10 min. The supernatant was diluted 50 times with ultrapure water, filtered through a 0.22 mm membrane, transferred to the 2 mL injection vial and analyzed the volatile fatty acid using a high-performance ion chromatography analyzer (ICS-3000, Thermo Scientific, USA).
Pyrosequencing of 16S rRNA Amplicons
The cecum and colon contents were removed from the − 80℃ refrigerator and the total bacterial DNA was extracted in accordance with the manufacturer's instructions of FastDNA® SPIN for soil kit (MP Biomedicals, USA). The integrity and purity of the DNA was determined using 1% agarose gel electrophoresis and a NanoDrop 2000 spectrophotometer (Thermo Scientific, USA). The DNA concentration was quantified precisely by Qubit Fluorometer (Thermo Scientific, USA), and the samples were diluted to 1 ng/µL with sterile water. The specific primers with marker sequences (338F: 5'-ACTCCTACGGGAGGCAGCAG-3' and 806R :5'-GGACTACHVGGGTWTCTAAT-3') were employed for PCR amplification of the variable region of 16S rRNA gene V3-V4. The total volume of PCR reaction system is 30 µL, consisting of Phusion® High-Fidelity PCR Master Mix with GC Buffer 15 µL, Phusion® High-Fidelity DNA polymerase 0.5 µL (New England Biolabs, USA), 0.2 µmol/L upstream and downstream 1 µL each, 10 ng/µL genomic DNA 2 µL as well as 10.5 µL of sterilized ultrapure water. The amplification procedure was performed as follows: 98°C pre-denaturation for 1 min, 98°C denaturation for 10 s, 50°C annealing for 30 s, 72°C extension for 30 s totaling 30 cycles, 72°C stable extension for 5 min, and lastly stored at 4°C (PCR instrument: ABI GeneAmp® Model 9700, USA).
The PCR products from the samples were mixed and examined by 2% agarose gel electrophoresis; then quantified by Quantus™ Fluor-ST Fluorescence Quantification System (Promega, USA). The PCR products were mixed in equal amounts following the concentration of the obtained PCR products, re-electrophoresed on 2% agarose gels, the target product bands were recovered using QIAquick Gel Extraction Kit (Axygen, USA). The library construction was performed by TruSeq® DNA PCR-Free DNA library kit (Illumina, USA), and sequenced by HiSeq2500 PE250 after Qubit and qPCR quality control. The raw tags were obtained by splicing the sample reads with FLASH software (http://www.cbcb.umd.edu/software/flash, V1.2.10); the raw tags were filtered and processed to obtain the high-quality effective tags by following the tags quality control process of QIIME [22]. a) truncate the raw tags from the first low-quality site where the number of consecutive low-quality bases (≤ 19) reaches a set length of three; b) filter out the tags with consecutive high-quality bases less than 75% of the tags length from the intercepted tags dataset. The sequences were OTU (operational taxonomic unit) clustered and chimeras were removed based on 97% similarity [23] using UPARSE [24] software (http://drive5.com/uparse/, V7.1) to generate OTU. The RDP classifier [25] (http://sourceforge.net/projects/rdp-classifier/, V2.2) Bayesian algorithm was employed to annotate the OTU sequences against the Silva 16S rRNA database (http://www.arb-silva.de, V138) for species classification (comparison threshold of 70%) and the community composition of each sample at each taxonomic level was counted.
Statistical analysis
All data were initially organized by Excel software and then statistically analyzed using the GLM model data in SAS software (V9.2, 2008). The dietary treatment as a fixed effect and the block as a random effect. The replicate as a unit to analyze the data concerning growth performance and apparent digestibility of nutrients in piglets, the diarrhea rate of piglets was analyzed by the chi-square test. The data on piglet serum, digestive enzymes, intestinal morphology and so on were analyzed based on the individual as a unit. Multiple comparisons and significance of differences among groups were performed using the Turkey method, means and standard errors were calculated using LSMEANS method; P < 0.05 was regarded as statistically significant, 0.05 < P ≤ 0.10 was indicative of a differential trend.
For microbiota profiling, statistical analysis of α-diversity, including the Shannon and Simpson indices reflecting microbial diversity and the Chao and ACE indices indicating bacterial abundance [26], was performed by mothur (http://www.mothur.org/wiki/, V1.30.1). The Circos-0.67-7 software (http://circos.ca/) was employed to analyze the relationship between samples and species by visualizing circle plots. The Unifrac distance was computed based on the species differences between the sample sequences and calculated the β-diversity distance by QIIME (http://qiime.org/), the principal coordinates analysis (PCoA) and partial least squares discriminant analysis (PLS-DA) were conducted using R software (V3.4.4). The analysis of between-treatment variance (One-way ANOVA) of α-diversity and β-diversity was performed using R software. The Analysis of similarities (ANOSIM) was used to compare the significance of differences in community structure between treatments. The LEfSe tool (Linear discriminant analysis [LDA] and effect size) was applied to analyze the core flora (LDA score > 3.0) in the cecum and colon of piglets (http://huttenhower.sph.harvard.edu/galaxy/root?tool_id=lefse_upload). The intergroup species differences were analyzed at each taxonomic level using R software based on the obtained bacterial community abundance data and plotted with the " vegan" and "ggplot2" packages.