4.1 Animals
The Institutional Animal Care and Use Committee of Huazhong Agricultural University (HZAUCH-2018-008), Wuhan, China) approved all the animal procedures, and all methods were performed in accordance with the relevant guidelines and regulations.
Newly hatched chickens (Turpan cockfighting × White Leghorn) were reared under similar husbandry conditions on the chicken farm of Huazhong Agricultural University. At the age of 1, 4, and 12 months, 120 chickens were randomly selected for each time point. Based on the abdominal fat index, the chickens at each time point were categorized into two groups, namely the high abdominal fat deposition group (H) and the low abdominal fat deposition group (L) (n = 10, 5 males and 5 females). For the fecal microbiota transplantation (FMT) experiment, the chickens with high body weight and low abdominal fat deposition were selected as FMT donors. 60 one-day-old white feather broilers were selected as recipients.
4.2 Selection of FMT donors
Two adult female white Leghorn chicken × Turpan fighting chicken possibly having high or low abdominal fat deposition were scanned with computed tomography (CT) instrument (Aquilion PRIME Tsx-303A, Canon Medical, Japan). Pari software was used to mark the abdominal fat in different frame images of each chicken (Fig. S1), and then Python language was used to write programs to analyze the images and calculate the volume of the body and abdominal fat of each chicken. The volume of the body was 2.22 dm3 and 2.50 dm3, and the volume of abdominal fat was 0.06 dm3 and 0.15 dm3, respectively. Similarly, the volume percentage of abdominal fat was 2.66% and 5.92%. The chicken with less abdominal fat volume percentage was selected as FMT donor. After FMT experiment, the two chickens were dissected to get the abdominal fat weight and index. The abdominal fat weight was 74.3 g and 161.2 g, and the abdominal fat index was 3.12% and 5.78%, respectively, which are consistent with the CT results and indicated that the FMT donor has low abdominal fat deposition.
4.3 Preparation of fecal suspension
Every morning, once the donor chickens defecated, the white part of the fecal materials was removed as it contains uric acid. Then 10 g of feces were collected in the sterile tube (50 mL) and gently mixed with 60 mL of 0.75% normal saline. The mixture was kept on the ice for settling down the precipitates. The supernatant was obtained and filtered with the sterile gauze to get fecal suspension.
4.4 Animal treatment
60 one-day-old white feather broilers were selected as recipients and randomly divided into FMT group and control group (n = 30). Broilers in the FMT group were orally administrated with 1 mL fecal microbiota suspension, while 1 mL 0.75% saline was used as a substitute in the control group for 28 days. At the age of 42 days, they were sacrificed and samples were collected.
4.5 Sample collection
After fasting for 12 hours, the chickens were weighed and sacrificed, then blood, liver, abdominal adipose tissue, and left cecum were collected. The abdominal adipose tissue was measured as well. For gut microbiota analysis, the cecal content (1 to 1.5 g per bird) was collected into two sterilized centrifuge tubes (1.5 mL) and snap-frozen in liquid nitrogen, then stored at -80°C for sequencing. For analysis of lipometabolic parameter, blood samples (3 mL per bird) were centrifuged at 3,000× g at 4°C for 15 min to get the serum, and then it was stored at − 80 ℃ for subsequent blood biochemical index analysis. For histo-morphological analysis, freshly harvested liver and abdominal adipose tissues were fixed in 4% paraformaldehyde solution. For gene expression analysis, parts of freshly harvested liver and abdominal adipose tissues were snap-frozen in liquid nitrogen and then stored at -80°C.
4.6 Muscle or abdominal fat index calculation
The muscle or abdominal fat index was calculated using the following formula: muscle index = muscle weight (g)/ body weight (g) × 100%, abdominal fat index = abdominal fat weight (g)/ body weight (g) × 100%.
4.7 16S rRNA and Metagenomic genes sequencing
Microbial genomic DNA was extracted from the chicken’s cecal content using Fast DNA SPIN extraction kit (MP Biomedicals, Santa Ana, CA, USA), according to manufacturer’s instructions. The hypervariable region V3-V4 of the bacterial 16S rRNA gene was amplified with primer pairs 338F (5'-ACTCCTACGGGAGGCAGCAG-3') and 806R (5'-GGACTACHVGGGTWTCTAAT-3'). The PCR amplification of the 16S rRNA gene was performed as follows: an initial denaturation (3 min) at 95 ℃ following 27 cycles of denaturing (30 s) at 95 ℃, annealing (30 s) at 55 ℃, extension (45 s) at 72 ℃, and single extension (10 min) at 72 ℃, and ended at 4 ℃. The PCR product was extracted from 2% agarose gel and purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA) according to manufacturer’s instructions and quantified using Quantus™ Fluorometer (Promega, USA). Illumina MiSeq PE300 platform (Illumina, San Diego, USA) was used for 16S rRNA gene sequencing. For 20 chickens at the age of 4 months with metagenomic sequencing, the same DNA extract was fragmented to an average size of about 400 bp using Covaris M220 (Gene Company Limited, China) for paired-end library construction, which was constructed using NEXTFLEX Rapid DNA-Seq (Bioo Scientific, Austin, TX, USA). Illumina NovaSeq platform (Illumina, San Diego, CA, USA) was used for metagenomic sequencing.
4.8 16S rRNA gene sequencing data processing
The raw 16S rRNA gene sequencing reads were demultiplexed, quality-filtered by fastp version 0.20.0, and merged by FLASH version 1.2.7. Operational taxonomic units (OTUs) with 97% similarity cutoff were clustered using UPARSE version 7.1, and chimeric sequences were identified and removed. The taxonomy of each OTU representative sequence was analyzed by RDP Classifier version 2.2 against the 16S rRNA database (Silva 132) using a confidence threshold of 0.7. For α and β diversity measurements, the sequencing depth was minimized by subsampling the readings of each sample. The lowest valid reads of cecal microbiota of high and low abdominal fat deposition chickens at the age of 1 month were 25,339, the lowest effective reading of cecal microbiota of high and low abdominal fat deposition chickens at the age of 4 months was 30,671, and the lowest effective reading of cecal microbiota of high and low abdominal fat deposition chickens at the age of 12 months was 45,053. Similarly, the lowest valid reads of cecal microbiota in the control and FMT chickens were 14,960. The α-diversity was described using the Shannon index and Chao index. Principal coordinates analysis (PCoA) based on Bray-Curtis was used to estimate the dissimilarity in the community structure. The community composition at the phylum level and the change of abundance at the genus level were visualized by bar chart and histogram. Linear discriminant analysis effect size (LEfSe) was performed to detect differentially abundant taxa across groups using the default parameters linear discriminant analysis (LDA > 2).
4.9 Metagenomic sequencing data processing
The low-quality reads (length < 50 bp or with a quality value < 20 or having N bases) were removed by fastp (https://github.com/OpenGene/fastp, version 0.20.0). Reads were aligned to the chicken genome by burrows-wheeler alignment (BWA) tool (http://bio-bwa.sourceforge.net, version 0.7.9a), and any hit associated with the reads and their mated reads were removed. The optimized sequence was spliced and assembled, and contigs ≥ 300 bp were selected as the final assembly result, and then the contigs were used for further gene prediction and annotation. Open reading frames (ORFs) from each assembled contig were predicted using MetaGene (http://metagene.cb.k.u-tokyo.ac.jp/). The predicted ORFs with length ≥ 100 bp were retrieved and translated into amino acid sequences. A non-redundant gene catalog was constructed using CD-HIT (http://www.bioinformatics.org/cd-hit/, version 4.6.1) with 90% sequence identity and 90% coverage. Reads after quality control were mapped to the non-redundant gene catalog with 95% identity using SOAPaligner (http://soap.genomics.org.cn/, version 2.21), and gene abundance in each sample was evaluated. Public data used for taxonomic analysis and gene functional classification included the integrated NCBI-NR database, KEGG database, and CAZy database. Amino acid sequence of non-redundant gene was aligned to NR database and KEGG database respectively with an e-value cutoff of 1e− 5 using Diamond (http://www.diamondsearch.org/index.php, version 0.8.35), and obtained the species annotation and KEGG function corresponding to the gene. Carbohydrate-active enzymes annotation was conducted using hmmscan (http://hmmer.janelia.org/search/hmmscan) against CAZy database (http://www.cazy.org/) with an e-value cutoff of 1e− 5.
4.10 Blood parameters analysis
For the analysis of different blood parameters, the serum concentrations of triglycerides (TG), total cholesterol (TC), high density lipoprotein cholesterol (HDL-C), and low density lipoprotein cholesterol (LDL-C) were determined using a Rayto Chemistry Analyzer (Chemray 800, China) according to the manufacturer’s instructions with the commercial diagnostic kits (Shenzhen Rayto Life Science Co., Ltd). Briefly, the serum samples were thoroughly mixed with the reaction solution in the recommended proportion and maintained at 37 ° C for 10 minutes. Finally, the absorbance for each sample was measured, and the total concentrations were calculated according to the following formula. Total concentrations = Absorbance of sample / Absorbance of calibration solution × Calibration concentrations (mmol per liter).
4.11 Hematoxylin and eosin (HE) staining
For morphological observation, liver and abdominal fat tissue samples were embedded in paraffin, and the sections were prepared. Liver tissues were cut into 3 µm thick sections, and abdominal fat tissues were cut into 7 µm thick sections with a rotary slicer (LEICA 819, Leica, Germany). HE staining was performed according to the routine protocol, and the stained tissue sections were examined with the light microscope (BH-2, Olympus, Japan) using a digital camera (DP72, Olympus, Japan).
4.12 Real-time quantitative polymerase chain reaction (qPCR)
In order to detect the expression of fat metabolism related genes on mRNA level, total RNA was extracted from abdominal adipose and liver tissues using Trizol reagent (Takara, Japan) following the instructions of the manufacturer. RNA (1µg) from each sample was reverse transcribed into cDNA using the PrimeScript™ RT reagent Kit with gDNA Eraser (Takara, Japan). The qPCR reaction mixture (10 µL) consisted of 5 µL of SYBR (Takara, Japan), 0.4 µL of forward and reverse primer, 3.2 µL of ddH2O, and 1 µL of template cDNA. The qPCR reaction is carried out on Bio-Rad CFX Connect real-time qPCR detection system (Bio-Rad, Hercules, CA, USA). The steps are as follows: 5 min pre-denaturation at 95°C, following 30 s denaturation at 95°C (40 cycles), 30 s annealing at 60°C, and 15 s elongation at 72°C. The sequences of primers were listed in Table 1 with reference gene (β-actin). Gene expression levels were quantified using the 2 −ΔΔCT method.
Table 1
Primers used for real-time qPCR
Gene
|
Forward Sequence (5’-3’)
|
Reverse Sequence (5’-3’)
|
Gene Bank No.
|
β-actin
|
TTGTTGACAATGGCTCCGGT
|
TCTGGGCTTCATCACCAACG
|
NM_205518.2
|
ACC
|
TCCAGCAGAACCGCATTGACAC
|
GTATGAGCAGGCAGGACTTGGC
|
NM_205505.2
|
FAS
|
GCTCTGCGTCTGCTTCAGTCTAC
|
GGTACAGGACTCTGCCATCAATGC
|
NM_205155.4
|
LPL
|
TGGACATTGGTGACCTGCTTATGC
|
TCGCCTGACTTCACTCTGACTCTC
|
NM_205282.2
|
ACSL1
|
GACTAATGGTCACAGGAGCAGCAC
|
CCAGGCATTGACAGTGAGCATCC
|
NM_001012578.2
|
FADS1
|
CCGTGCCACTGTGGAGAAGATG
|
GCCTAGAAGCAACGCAGAGAAGAG
|
XM_040673219.1
|
CYP2C45
|
AACAAGCACCACCACACGATACG
|
GGTCAGCCACGCAAGGTCTTC
|
NM_001001752.3
|
APOAI
|
GTGACCCTCGCTGTGCTCTT
|
CACTCAGCGTGTCCAGGTTGT
|
NM_205525.5
|
PPARα
|
TGCTGTGGAGATCGTCCTGGTC
|
CTGTGACAAGTTGCCGGAGGTC
|
XM_040699549.1
|
CPT-1
|
ACAGCGAATGAAAGCAGGGT
|
GCCATGGCTAAGGTTTTCGT
|
NM_001012898.1
|
LEPR
|
CACTCGCTGGGAACACTTGA
|
TTCAGCAGCCCATCGTTTCT
|
NM_204323.2
|
JAK2
|
GAGCGTGAGAATGCCACTGAC
|
TGGAGGACAGCACTTGATGAAC
|
NM_001030538.3
|
STAT3
|
GCCGAATCACAACTACAGACTC
|
CTGACTTTGGTGGTGAACTGC
|
NM_001030931.3
|
HSL
|
GAGGCACAGCGTCTTCTTTAGG
|
GGCACGAACTGGAACCCGAG
|
XM_040695201.1
|
4.13 Immunohistochemistry (IHC)
Following the steps described in earlier studies, immunohistochemical staining was performed to observe the protein distribution and expression in the liver. Briefly, the sections were dewaxed twice in xylene and rehydrated in graded series ethanol. The antigen was retrieved in sodium citrate buffer (pH 6.0) using a microwave oven (MYA-2270M, Haier, Qingdao, China) for 18 min, i.e., three min at 700 W and fifteen min at 116 W, and then cooled for 2 to 3 h at room temperature. Endogenous peroxidase was inactivated with 3% hydrogen peroxide (H2O2), and tissue sections were incubated with 5% bovine serum albumin (BSA) (boster, China) at 37 ℃ for 30 minutes to block nonspecific binding sites. Then, the sections were incubated with primary antibodies of rabbit anti-JAK2 (1:100) (A11497, ABclonal Technology, Wuhan, China), rabbit anti-p-JAK2 (1:100) (AP0531, ABclonal Technology, Wuhan, China), rabbit anti-STAT-3 (1:100) (A1192, ABclonal Technology, Wuhan, China) and rabbit anti-p-STAT3 (1:100) (AP0474, ABclonal Technology, Wuhan, China). Subsequently, the horseradish peroxidase (HRP)-conjugated secondary antibody (Proteintech, China) was used to incubate the tissue sections for 30 min at 37°C. After diaminobenzidine (DAB) (Proteintech, China) staining, the sections were counterstained with hematoxylin, cleaned and dehydrated until they became transparent, and finally sealed with neutral gum and coverslips. Finally, we used a light microscope (BH-2, Olympus, Japan) with a digital camera (DP72, Olympus, Japan) to examine the sections.
4.14 Statistical analysis
Under 10 x 20 microscope, 10 abdominal fat HE stained sections were selected from each group, and 5 visual fields were randomly selected for the image acquisition. The average diameter of abdominal fat adipocytes was measured with image pro plus 6.0 (Media Cybernetics, USA). Under a 10 x 40 microscope, 10 liver immunohistochemical sections were selected from each group and five positive visual fields were randomly selected from each section for the image acquisition. Image Pro Plus 6.0 was used to calculate the integral optical density of positive signals. GraphPad Prism 6.0 (Media Cybernetics, USA) was used to analyze the test data. The measurement data were expressed as mean ± standard error of the mean (mean ± SEM). The statistical significance of the mean values in the comparisons of two groups was determined using Student’s t-test. Value of p < 0.05 was considered statistically significant.