Birds and conditions
A total of 240-pair parent White King pigeons, 60 wk of age, were obtained from a commercial pigeon farm (Wenzhou, China). An artificial aviary with a perch and a nest was provided for each pair. Parent pigeons were randomly allocated to four treatments which consisted of 6 replications of 10-pair pigeons each. The four treatments are: control group, 1% LA supplementation group (LA1%), 2% LA supplementation group (LA2%), 4% LA supplementation group (LA4%). The LA (༞99.0% purity) was obtained from a commercial supplier (Hebei bawei biotechnology co., LTD, Handan, China). All experimental diets were isonitrogenous and isocaloric. The LA source was substituted for colleseed oil by the same weight to equalize the total fat level among diets [17]. The ingredients, nutrient levels and analyzed LA content of experimental diets for parent pigeons were shown in Table 1. The fatty acid composition of each diet analyzed according to Xie et al. [18] were shown in Additional file 1. The pigeons were given to water ad libitum and were fed twice daily (7:00 a.m. and 3:00 p.m.) throughout the experiment.
Table 1
Ingredient compositions and nutrient levels of experimental diets1 for maternal pigeons.2 (on as-fed basis)
Items | Control | LA1% | LA2% | LA4% |
Ingredients of a whole-grain form feed (%) | | | | |
Corn | 56.27 | 56.27 | 56.27 | 56.27 |
Pea | 28.12 | 28.12 | 28.12 | 28.12 |
Wheat | 5.63 | 5.63 | 5.63 | 5.63 |
Sorghum | 5.63 | 5.63 | 5.63 | 5.63 |
Colleseed oil | 4.35 | 3.26 | 2.17 | 0 |
Linoleic acid | 0 | 1.09 | 2.18 | 4.35 |
Total | 100.00 | 100.00 | 100.00 | 100.00 |
Calculated nutrients3 (%) | | | | |
Metabolizable energy (MJ/kg)4 | 13.84 | 13.84 | 13.84 | 13.84 |
Crude protein | 12.63 | 12.63 | 12.63 | 12.63 |
Analyzed linoleic acid | 1.82 | 2.75 | 3.68 | 5.53 |
Ingredients of grit meal (%) | | | | |
Limestone | 52.93 | 52.93 | 52.93 | 52.93 |
Shell meal | 28.10 | 28.10 | 28.10 | 28.10 |
Yellow mud | 14.05 | 14.05 | 14.05 | 14.05 |
Salt | 1.41 | 1.41 | 1.41 | 1.41 |
Ferrous sulfate (monohydrate) | 0.23 | 0.23 | 0.23 | 0.23 |
Premix5 | 3.28 | 3.28 | 3.28 | 3.28 |
Total | 100.00 | 100.00 | 100.00 | 100.00 |
Calculated nutrients2 (%) | | | | |
Calcium | 27.79 | 27.79 | 27.79 | 27.79 |
Phosphorus | 0.01 | 0.01 | 0.01 | 0.01 |
Sodium | 0.03 | 0.03 | 0.03 | 0.03 |
Chlorine | 0.59 | 0.59 | 0.59 | 0.59 |
Ferrum | 0.30 | 0.30 | 0.30 | 0.30 |
1Control = control group; LA1% = 1% linoleic acid addition group; LA2% = 2% linoleic acid addition group; LA4% = 4% linoleic acid addition group. |
2All feeds were fed in a whole-grain form at 07:00 and 15:00 a day, and grit meal was offered to the birds on a continuous basis. |
3Nutrient values were calculated from tables of feed composition and nutritive values in China (twenty-eighth edition, 2017). |
4Metabolizable energy values determined in pigeons were calculated from those reported for chickens in accordance with a previous study (Hullar et al., 1999), which observed that the metabolizable energy values of feed in pigeons did not differ significantly from those in chickens. |
5The premix provided the following per kg of diet: Vitamin A 5000 IU, Vitamin E 50 IU, Vitamin D3 2000 IU, copper sulfate 15 mg, manganese sulfate 45 mg, zinc sulfate 90 mg. |
Each pair of parent pigeons laid two eggs in a nest. Eggs were picked out by hand and meanwhile fake eggs were put into nests to meet parents’ broodiness. The picked eggs were transferred to an incubator for an 18-day artificially incubation. On the day of hatching, 480 artificially hatched pigeon squabs with similar BW were selected from this pigeon farm. These squabs were randomly pair-matched and assigned into the nests of the selected parent pigeons to replace the fake eggs. Each parent pair adopted a pair of squabs. Parent pigeons were fed with experimental diets continually in the next 21 days. Squabs were fed with crop milk which was secreted by parent pigeons. The ambient temperature was 18 to 26 °C. The relative humidity was 60 to 70%. And the photoperiod was 12 L:12 D throughout the total experiment period.
Sample collection
On the day 21 posthatch, six squabs per treatment (one for each replication) were selected for sampling. The squabs were fasted for 12 h before weighing and slaughter. And they were sedated before cervical dislocation. Ileal contents were sampled into sterile plastic tubes with a sterile spatula, immediately frozen in liquid nitrogen, then stored at -80 °C for subsequent analyses. Three segments of the small intestine (duodenum, jejunum and ileum) were sampled, flushed with 0.9% saline to remove all the content, and processed for morphological examination. The collected segments were the loop of the duodenum (duodenum), the tract before Meckel’s diverticulum (jejunum), the final segment of small intestine (ileum). The mucosa of each intestinal segment was scraped off with a glass slide carefully, frozen in liquid nitrogen rapidly, and then stored at -80 °C for subsequent analyses.
Secretory IgA and cytokine analysis
The homogenates of intestinal mucosa were prepared with phosphate buffered saline for secretory IgA (sIgA) and cytokine analyses. The sIgA concentration and the concentrations of tumor necrosis factor-α (TNF-α), interleukin-1β (IL-1β), IL-6, IL-10, and IL-4 were determined by absorbance changes at a wavelength of 450 nm with commercial ELISA kits (Shanghai Mlbio Institute, Shanghai, China) according to the manufacture’s protocol. The final sIgA and cytokine concentrations were expressed as units per milligram of protein.
RNA extraction and quantitative PCR analysis
Total RNA of intestinal mucosa was extracted using TRIzol procedure (Invitrogen, Carlsbad, CA) according to the instructions of the manufacturer. The extracted RNA was quantified by NanoDrop ND-2000 spectrophotometer (NanoDrop Technologies Inc., Wilmington, USA). The RNA integrity was verified by native RNA electrophoresis on 1.0% agarose gel. The complementary DNA was synthesized from 2-µg total RNA by M-MLV reverse transcriptase (Takara, Dalian, China) with oligo dT-adaptor primer at 42 °C for 60 min following the protocol of the manufacturer. The abundance of mRNA was assayed on a StepOne Plus Real-Time PCR system (ABI 7500, Applied Biosystems, Foster City, CA). The specific primers uesd were indicated in Table 2. The SYBR Green Realtime PCR Master Mix (Toyobo Co., Ltd., Osaka, Japan) was used for PCR reaction, consisting of an initial DNA denaturation of 95 °C for 60 s, followed by 40 cycles of 95 °C for 15 s, and 60 °C for 60 s. The GAPDH was considered as an appropriate endogenous reference. The average gene expression relative to the endogenous reference for each sample was calculated according to the 2−ΔΔCt method [19]. The calibrator for each gene in study was the average ΔCt value of control group.
Table 2
Primers used for quantitative real-time PCR analysis of gene expression in domestic pigeons.
Gene1 | Primer | Sequence (5’→3’) | Amplicon size (bp) | Accession number |
CLDN2 | CLDN2-F | GTGCAGATGGGAACAAGGT | 119 | XM_021283269.1 |
CLDN2-R | GAGCCAAGGAAGCTACGG |
CLDN3 | CLDN3-F | ACCTCATCCCCGTCTCCT | 109 | XM_005515008.2 |
CLDN3-R | CAGCCCACGTAGAGCGA |
OCLN | OCLN-F | CAGGACGTGGCAGAGGA | 105 | XM_005509325.2 |
OCLN-R | GTGGAAGAGCTTGTTGCGT |
GAPDH | GAPDH-F | CCCAATGTCTCTGTTGTGG | 116 | NM_001282835.1 |
GAPDH-R | TATGCCAGGATGCCCTT |
1CLDN2 = Claudin2; CLDN3 = Claudin 3; OCLN = Occludin; GAPDH = glyceraldehyde-3-phosphate dehydrogenase. |
Intestinal morphological examination
Approximately 1-cm intestinal samples of each segment were collected and fixed in 10% neutral-buffered formalin solution. Each sample was dehydrated, cleared, and embedded in paraffin. Serial sections (5 µm) were placed on glass slides and submitted to two different staining. The hematoxylin and eosin staining was for the identification of villus height, villus area, crypt depth, and enterocyte. The periodic-acid Schiff staining was for the identification of goblet cell (GC). The examination was performed with light microscopy (Nikon Corp., Tokyo, Japan) using the Image-pro plus 6.0 (Media Cybernetics, Inc., Rockville, MD, USA). The density of enterocyte was expressed as number per millimeter of epithelial layer length. The density of GC was expressed as number per 100 enterocytes according to Tamura et al. [20].
16S rDNA high-throughput sequencing
The total microbial genomic DNA was extracted from the ileal contents of squabs from each treatment using the E.Z.N.A.®Stool DNA Kit (Omega Bio-Tek, USA). The concentration and quality of the extracted DNA were assessed with 2% agarose gel electrophoresis and a NanoDrop ND-2000 spectrophotometer (NanoDrop Technologies Inc., Wilmington, USA). The V3-V4 hyper-variable region of the 16S rRNA genes was amplified by specific degenerate primers (338F: ACTCCTACGGGAGGCAGCAG, 806R: GGACTACHVGGGTWTCTAAT). The PCR program was conducted following the conditions: initial denaturation step, 98 °C, 30 s; denaturation, 35 cycles, 98 °C, 10 s; annealing, 54 °C, 30 s; elongation, 72 °C, 45 s; and final extension, 72 °C, 10 min. The PCR amplification products were detected by 2% agarose gel electrophoresis. The target fragments were recovered using AMPure XT Beads (Beckman Coulter Genomics, Danvers, MA, USA). Qubit was used to quantify the purified PCR products, and the qualified library concentration should be more than 2 nM. Sequencing was performed with V3 chemistry and generated 2 × 300 bp double-end reads using MiSeq Reagent Kit V3 with the MiSeq Illumina platform (Illumina, Santa Clara, CA, USA) at Hangzhou Personal Biotechnology Co., Ltd. (Hangzhou, China).
Statistical analyses
The data obtained from sIgA and cytokine analysis, quantitative PCR analyses, and intestinal morphological examination were subjected to one-way analysis of variance in SPSS 24.0 (SPSS Inc., Chicago, USA) for Windows. The significance level was set at P < 0.05. When significant differences were found, Tukey post-hoc tests were performed. Plotting was placed with GraphPad Prism 7.0 (GraphPad Software Inc., San Diego, CA, USA).
The raw sequence data obtained from the Illumina MiSeq platform were quality-filtered and demultiplexed using Quantitative Insights into Microbial Ecology (QIIME), version 1.8.0-dev (http://qiime.org/index.html). After trimming and assembling, the final clean data were obtained and assigned into bacterial operational taxonomic units (OTUs) using the UCLUST function in QIIME (http://qiime.org/scripts/pick_otus) according to a 97% similarity threshold. The data were categorized with the Ribosomal Database Project (RDP), version 2.2 (http://sourceforge.net/projects/rdp-classifier/). The alpha diversity measures, including the observed speacies, Chao 1, and Shannon indices were calculated using MicrobiomeAnalyst (https://www.microbiomeanalyst.ca/MicrobiomeAnalyst/faces/home.xhtml). Venn diagram was performed in Venny 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/index.html). Principal coordinate analysis (PCoA) was conducted using Galaxy (http://huttenhower.sph.harvard.edu/galaxy/) according to Unweighted Unifra representing Beta diversity values. Kolmogorov-Smirnov program is used to verify the normality of data. Data in line with normal distribution were counted by one-way ANOVA, and multiple comparisons were conducted by Tukey post-hoc tests. Nonnormal distribution data were analyzed by Kruskal-Wallis test and Duncan’s multiple range tests. P < 0.05 means significant difference.