1. Introduction
Pork has long been considered the staple meat source for Chinese daily life [
1]. Traditional Chinese black pigs generally exhibited higher backfat, intramuscular fat (IMF), and protein content, which provided a flavorful taste compared with the lean-type European and American pigs and gained the raising popularity in modern lifestyle [
2,
3,
4,
5]. Yushan black pig is a quintessential representation of traditional Chinese pigs, characterized by a big belly, higher protein (20.1–24.3%), intramuscular fat (2.01–3.61%), and unsaturated fatty acids (47.0–57.5%) contents [
6]. However, the thickness of backfat and abdominal fat restricted feed efficiency and significantly increased feed consumption during the rearing process. Therefore, proper strategies that effectively reduce backfat and abdominal fat content significantly promote the utilization and popularity of Yushan pigs.
Dietary energy provision is the key factor that regulates the fat deposition process, and it varies efficiently due to the ingredient composition of feedstuffs. Previous studies showed a higher percentage of unsaturated fatty acids, such as C18:1 n-9, C18:2 n-6, C18:3 n-3, and total polyunsaturated fatty acids (PUFA), and a lower percentage of saturated fatty acids, including C16:0, C18:0, and C20:0, in the lower energy provision dietary treatment [
7,
8]. In addition, changes of feed nutritional components, particularly the dietary energy and fiber content, significantly impacted the nutrient digestibility and metabolic transformation and further impacted the fat deposition through modulating the physiological energy expenditure process [
9,
10].
Generally, brans and husks are representative sources of crude fiber, which effectively decline dietary energy content, remain resistant to digestive enzyme catalyzation, and reduce the absorption of fat and carbohydrates in animal GIT [
11]. Feedstuff consumed by pigs that contains higher dietary fiber has been reported to effectively reduce the risk of diabetes, obesity, and gastrointestinal disorders [
12]. Higher content of dietary fiber also functionally shaped the intestinal bacterial community and regulated the digestibility of other nutritional components. Alterations of bacterial communities in the pig gastro-intestinal tract (GIT) substantially influenced host physiological and immunological processes, which causatively regulated the growth and shaped the meat quality and flavor [
13,
14,
15]. Additionally, an increase in dietary fiber content has been reported with a notable impact on hepatic lipo-metabolic-related gene expression and the fat deposition and composition of pigs [
16]. Despite these findings, the underlying mechanism of increasing dietary fiber content on body composition was still of limited information. Further research is needed to fully understand the intricate relationship between dietary fiber, body composition, and overall health outcomes.
Therefore, this study aims to investigate the regulatory effects of lowering the energy provision induced by the higher content wheat bran, rice bran, and corn bran on the productive performance, meat quality, and fat acid composition. We hypothesized that higher dietary fiber content effectively decreased the backfat and abdominal fat content of Yushan pigs through enhancing intestinal fiber-degradable bacteria, regulating metabolic circulation, and up-regulating the hepatic lipolysis-related gene expression.
2. Materials and Methods
Animal care and procedures followed The Chinese Guidelines for Animal Welfare, which were approved by the Animal Care and Use Committee of the Jiangxi Academy of Agricultural Sciences; approval number 2024-JXAAS-XM-16.
2.1. Experiment Animals and Management
A total of 18 120-day-old Yushan pigs with similar initial body weight were selected and randomly allotted into low dietary fiber content (high energy, HE) and high dietary fiber content (low energy, LE) treatments for a 130-day-long feeding process, which contained a 10-day adaptive stage and a 120-day feeding stage. Each treatment contained 9 pigs, with each pig considered a replicate and reared in an individual pen that was 220 cm long × 90 cm wide × 130 cm heigh. Daily feed intake and daily weight gain were recorded individually. All pigs were reared in the same piggery, with all pens in a bidirectional arrangement. Feed was provided for each treatment twice per day at 07:30 and 18:30, and ingredients and chemical nutrients are displayed in
Table 1. All pigs were allowed to get feed and water ad libitum. Productive procedures, which included immunization and cleaning, were conducted based on the regular rearing program. The relative humidity and temperature of the feeding piggery were maintained at 60–65% and 20–27 °C humidity and temperature, respectively.
2.2. Productive Performances and Carcass Performances Measurement
The feed intake (FI) of each pig was recorded and calculated by the deviation between the feeding scale and the residue before the next morning feed, and the average daily feed intake (ADFI) was calculated as the rate between total feed intake and the feeding days. The body weight of each pig was weighed at the end of the trial, and the average daily weight gain (ADG) was calculated through the following equation:
Next, the feed conversion ratio (FCR) between feed consumption and egg production was calculated based on the following equation:
At the end of the trial, all 18 pigs were slaughtered using CO
2 stunning and bleeding in accordance with Chinese regulations for determining carcass performance of breeding pigs (NY/T 822-2004) [
17] before a 12 h fast but with free access to water.
Carcass-related indexes, including carcass and length and backfat thickness, were measured. Carcass length was measured from the anterior edge of the first cervical vertebra to the posterior edge of the last lumbar vertebra. Backfat thickness was measured between the last cervical and first thoracic vertebrae and between the last thoracic and first lumbar vertebrae. The average of each measurement was calculated as backfat thickness. Organs including the heart, liver, spleen, and kidney were acquired and weighed; the organ indexes were calculated as the rate between organ weight and carcass weight. Furthermore, abdominal fat, subcutaneous fat, lean meat, skin, and bones were separated and weighed. The fat, lean, skin, and bone percentages were calculated using the following equation:
2.3. Meat Quality and Lipid Composition Measurement
Meat characteristics, which include the shear force, dripping loss, meat protein content (%), and intramuscular fat (IMF) (%), were measured based on the methods introduced in Chinese technical regulation for the determination of pork quality (NY/T 821-2019) [
18]. Specifically, the meat was first cooked at a temperature of 72 °C and cut into rectangular, cooked meat sections (1 × 1 × 3 cm) when cooling to room temperature. Shear force was further measured perpendicular to the direction of fibers using a texture analyzer TA HD Plus (Stable Micro Systems Ltd., Surrey, UK) equipped with a Warner–Bratzler V-shaped shear blade (1.2 mm thick). Meat protein content was determined by the Dumas combustion method (method 992.15) [
19]. IMF was determined through the ether extraction method, which was introduced by J. Folch [
20].
The lipid composition was further analyzed according to the acetyl chloride methanol–methyl esterification method presented by S. Jaturasitha [
21]. Both saturated fatty acid (SFA) and unsaturated fatty acid (UFA) were measured using the gas chromatograph Agilent 8860 GC, Santa Clara, CA, USA. Parameters of chromatographic measurement were set as follows: A dicyanopropyl polysiloxane column (100 m × 0.25 mm, 0.20 mm) was applied with a temperature of the column oven set at 140 °C for 5 min, and then gradually increased to 240 °C at 4 °C/min. The injector and detector temperatures were set as 260 °C and 280 °C, respectively. The chromatographic peak area was obtained by the integral method and further used for quantification as the content of lipid composition.
2.4. Blood Metabolomic Measurements
Five milliliters of each blood sample were collected during exsanguinations (using 10 mL EDTA tubes (Yuli medical instrument Co., Ltd., Nanjing, China)). Plasma was obtained by centrifugation at 2500× g, 4 °C for 10 min, and followingly applied for the blood metabolic measurement through the LC/MS analysis method.
To be simply stated, an internal standard was primarily made through the dissolution of 100 µL plasma by 400 µL 80% methanol and 0.02 mg/mL L-2-chlorophenyl alanine. All samples were following cleaned using ultrasound at 40 kHz for 30 min at 5 °C and carefully transferred to sample vials for LC-MS/MS analysis after the 13,000× g centrifugation at 4 °C for 15 min. Chromatographic identification of the metabolites was performed using a Thermo UHPLC system equipped with an ACQUITY UPLC HSS T3 (100 mm × 2.1 mm i.d., 1.8 µm; Waters, Milford, CT, USA). Mass spectrometric data were further collected using a Thermo UHPLC-Q Exactive HF-X mass spectrometer equipped with an electrospray ionization (ESI) source in either positive or negative ion mode. All raw data were imported into the Progenesis QI 2.3 (Nonlinear Dynamics, Waters, CT, USA) for peak detection and alignment and displayed into a data matrix that consisted of retention time (RT), mass-to-charge ratio (m/z) values, and peak intensity.
2.5. Gut Microbiota Analysis
Contents in the cecal section of each pig were sampled for 16S rRNA microbial diversity analysis. Briefly, cecal content DNA from each sample was extracted using the cetyltrimethylammonium bromide and sodium dodecyl sulfate (CTAB/SDS) method described by F. Xue [
22], followed by the PCR amplification process using the primers designed through the V4 region of 520F (F: GTGCCAGCMGCCGCGGTAA) and 802R (R: GGACTACHVGGGTWTCTAAT). PCR amplification output of each sample was sequenced under the Illumina HiSeq 4000 platform (Illumina Inc., San Diego, CA, USA), and to control the raw tag quality a Quantitative Insights Into Microbial Ecology 2 (QIIME, V 2.0) package was used. Sequences within similarity > 97% were assigned to the same operational taxonomic unit (OTU). Alpha diversity, beta diversity, and functional prediction analyses were further conducted.
2.6. Liver Energy-Related Metabolic Gene Expression and Validation
Total RNA was first extracted from the liver samples of each pig using RNAiso Plus reagent (code No. 9109, Takara, Dalian, China). Purity and contamination of RNA were detected using 1% agarose gels and a NanoPhotometer® spectrophotometer (IMPLEN, Westlake Village, CA, USA). RNA concentration was measured using a Qubit® RNA Assay Kit with a Qubit® 2.0 Fluorometer (Life Technologies, Carlsbad, CA, USA), followed by selecting 3 μg of RNA for transcriptome sequencing.
Then, the transcriptomic sequencing processes, which contained the amplification, library construction, sequencing, filtration, and gene identification, were conducted. PCR products were purified using the AMPure XP system, followed by the construction and quality assessment of the library using the Agilent 2100 Bioanalyzer system. The library was sequenced using the Illumina NovaSeq 6000 platform (Illumina Inc., San Diego, CA, USA). The quality of the RNA sequences was checked using FastQC (v0.11.9, Babraham Institute,
www.bioinformatics.babraham.ac.uk (accessed on 21 May 2024)) and sequence adapters. Low-quality reads (read quality < 20) were removed, and the filtered sequenced reads were further mapped to the pig genome, followed by the quantification of the expression using the FeatureCounts (version: 1.5.0-p3) software.
The expression level of genes in each sample was first normalized and further conducted the differential analysis by the DESeq2 analysis package of R software (version 4.1.3, R Core Team, Vienna, Austria). Differential expressed genes (DEG) were filtered based on the criteria of |log2FoldChange|≥1 and
p ≤ 0.05. KOBAS 3.0 (kobas.cbi.pku.edu.cn) was used to perform gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG,
http://www.genome.jp/kegg/, accessed on 22 May 2024) pathway enrichment analyses of the DEGs and association genes. A
p value < 0.05 was considered significant in GO term enrichment and pathway analysis.
Total RNA was extracted using the RNA Simple Total RNA Kit (Tiangen, Tiangen Biochemical Technology Co., Ltd., Beijing, China) according to the manufacturer’s protocol. One microgram of total RNA was used to carry out reverse transcription using the Fast Quant RT Kit (Tiangen). Gene expression was quantitatively analyzed by qRT-PCR in a QuantStudio™ 5 Real-Time PCR System using AceQ qPCR SYBR Green Master Mix (Low ROX Premixed; Vazyme, Nanjing, China). The cycling conditions were as follows: 95 °C for 5 min, followed by 40 cycles of 95 °C for 10 s and 60 °C for 30 s. The relative mRNA expression of the target genes was calculated using the 2−ΔΔCt method after normalization by the levels of GAPDH (a constitutively expressed gene that was used as the internal control).
2.7. Statistical Analysis
Productive performance, meat quality-related parameters, and intramuscular lipid content were firstly assessed for normal distribution using the SAS procedure “proc univariate data = test normal” and subsequently carried out the student’s t-test using SAS 9.4 (SAS Institute, Inc., Cary, NC, USA). Significance would be considered when p < 0.05, while a tendency towards significance was considered when 0.05 ≤ p < 0.10.
The normal distribution of relative OTU abundances of each sample was first assessed using the SAS procedure “proc univariate data = test normal”, followed by the differential analysis using the student’s t-test (SAS version 9.4, SAS Institute Inc., Cary, NC, USA). Alpha diversity and beta diversity of our samples were calculated with QIIME (version 1.7.0) and displayed with R software (version 4.1.3, R Core Team, Vienna, Austria). Principal coordinate analysis (PCoA) was displayed using the ggplot2 package in R software. Spearman correlations between bacteria communities and production performance were assessed using the PROC CORR procedure of SAS 9.4. Finally, the correlation matrix was created and visualized in a heatmap format using R software (version 4.1.3, R Core Team, Vienna, Austria).
Multivariate analyses on plasma metabolomic results, including principal component analysis (PCA) and orthogonal correction partial least squares discriminant analysis (OPLS-DA), were conducted using SIMCA-P software (V 14.0, Umetrics, Umea, Sweden). Differentially expressed metabolites between two treatments were identified based on variable importance in projection (VIP) from OPLS-DA analysis and statistical analysis (VIP > 1 and
p < 0.05). Kyoto Encyclopedia of Genes and Genomes (KEGG,
http://www.genome.jp/kegg/, accessed on 22 May 2024) was used to view the enriched pathways of different metabolites.
3. Results
3.1. Differential Analysis of Dietary Fiber Content on Growth Performance, Carcass Performance, and Meat Quality
The productive performance, including the ADFI, ADG, and the FCR, were calculated and are displayed in
Table 2. A significant increment was discovered in ADFI (
p < 0.05) in LE treatment compared with HE, while a decreased tendency of FCR was observed in the high-energy feeding treatment compared with the low-energy level feeding treatment (0.05 <
p < 0.10). No other significant alterations were discovered for initial body weight, final body weight, or ADG.
All experimental pigs were chosen for the carcass performance measurement, and the results are shown in
Table 3. The carcass weight and fat percentage were significantly higher, while the skin percentage was significantly lower in the HE treatment compared with LE (
p < 0.05). No significant alterations were observed in leg weight, dressing percentage, lean percentage, and bone weight percentage (
p > 0.05). Further, meat qualities such as shear force, dropping loss, meat protein content, and intramuscular fat (IMF) content were measured and results are shown in
Table 3. The IMF content in HE treatment was significantly higher than that in LE treatment (
p < 0.05). No significant changes were found for other parameters (
p > 0.05).
3.2. Differential Analysis of High Fiber Content Dietary on Sarcous Fatty Acids Composition
The sarcous fatty acid composition, which includes both saturated and unsaturated fatty acids, was investigated, and the results are shown in
Table 4.
C16:0 content was the only saturated fatty acid that significantly increased in HE treatment compared with LE (p < 0.05). In addition, the higher fiber content treatment significantly increased the content of C18:1(cis) while significantly decreasing the C18:1(trans) content compared with the higher concentrate diet treatment. No other significant differences are observed for the residual fatty acids between LE and HE treatments (p > 0.05).
3.3. Plasma Metabolic Responses to High Fiber Content Dietary Treatment of Yushan Pigs
Metabolic determination revealed a total of 623 across all samples through the filtering method, all these metabolites are listed in
Table S1. All identified metabolites were chosen for the differential analysis of integrative metabolic alteration between the HE and LE treatment through PCA and OPLS-DA analysis. As shown in
Figure 1, PC1 and PC2 accounted for 31.5% and 14.7% of the total variation, respectively. Metabolites identified in the HE and LE treatment samples clustered into two clearly separated sections in both PCA and OPLS-DA analysis results, which means a significant difference between HE and LE treatments regarding plasma metabolites.
Furthermore, specifically differentially expressed metabolites between HE and LE treatments were identified based on the statistical standard of fold change > 2, VIP > 1, and
p < 0.05. As shown in
Table 5, a total of 10 down-regulated and 17 up-regulated metabolites were identified in HE treatment compared with LE. The up-regulated metabolites mainly contained cyclic ADP-ribose, 7-methylguanine, S-adenosylhomocysteine, adenosine, triiodothyronine, and hippuric acid, while the down-regulated metabolites mainly consist of DL-Arginine, ciprostene, propionylcarnitine, and desloratadine. No other significantly altered metabolites were found between the HE and LE treatments.
Functional enrichment and pathway analysis were applied based on the differentially identified compounds. All results are shown in
Figure 2 and
Figure 3. As
Figure 2 shows, differentially expressed metabolites are functionally enriched in the carbohydrate metabolism, energy metabolism, amino acid metabolism, cofactors, and nucleotide metabolism. Specifically, energy metabolism showed a significant decline in higher fiber dietary treatment compared with HE treatment (
p < 0.05). KEGG results showed that glycerophospholipid metabolism, nicotinate and nicotinamide metabolism, and bile secretion are the most enriched three pathways based on the differential metabolic content. Specifically, the pathway of regulation of the actin cytoskeleton is more enriched between LE and HE differential metabolites; however, the enriched number was low.
3.4. Effects of High Fiber Content Dietary Treatment on Gastrointestinal Bacteria Communities
Cecal bacterial communities were detected for the causative investigation on the productive and metabolic differences between HE and LE treatments. A total of 17 phyla and more than 250 genera were identified, and all information is shown in
Table S2. All identified microbial communities were selected for further analysis.
3.4.1. α-diversity
Alpha diversity was first applied to detect the modulative effects of energy on bacterial abundances, and results are shown in
Table 6. Bacterial diversities showed no significant alterations between HE and LE treatments through all parameters, including the Sobs, Shannon, Simpson, Chao, ACE, Pielou, and PD indexes (
p > 0.05).
3.4.2. β-diversity
Differential analyses on cecal bacteria between HE and LE treatments were followingly applied, and the result is displayed through PCoA. As shown in
Figure 4, PCoA axes 1 and 2 accounted for 71.19% and 14.69% of the total variation, respectively. Bacteria communities between HE and LE treatments could be clearly separated by PCo1 and PCo2, except HE-5.
Differential analysis of the relative abundances of cecal bacteria at phyla and genera levels was calculated, and the results are shown in
Table 7 and
Table 8. At the phyla level,
Firmicutes,
Bacteroidota, and
Spirochaetota accounted for the top three abundant phyla. Relative abundances of
Bacteroidota and
Spirochaetota significantly decreased, while
Firmicutes significantly increased in HE treatment compared with LE treatment (
p < 0.05). No other significantly altered phylum was observed between HE and LE.
At the level of genera, Bacteroides, Lactobacillus, Rummeliibacillus, Prevotellaceae, and Faecalibacterium contributed to the top 5 most abundant genera. Abundances of Bacteroides, Prevotellaceae_UCG, and Clostridium significantly increased in HE treatment (p < 0.05). Meanwhile, abundances of Lactobacillus, Ruminococcus, Romboutsia, Succinivibrio, Phascolarctobacterium, and Bifidobacterium are significantly proliferated in LE feeding treatment (p < 0.05). No significant changes were detected among other genera (p > 0.05).
3.5. Functional Prediction on the Differential Gut Microbiota
Predictive functions of the above-mentioned differential microbiota between HE and LE treatments were conducted through the Tax4Fun process [
23], and the result is shown in
Figure 5.
The predictive functional results include the metabolism, genetic information processing, environmental information processing, cellular processes, and organismal systems. Metabolism processes are the most impacted pathways, which include the carbohydrate metabolism, energy metabolism, amino acid, cofactors, and lipid metabolism. In addition, translation in genetic information processing and membrane transport in environmental information processing are also the most impacted functional processes enriched by the differential bacterial communities between high-fiber dietary treatment and high-concentrate diet treatments.
Moreover, the interactive effects between significantly altered gut bacteria and the significantly altered metabolites were selected for a correlation analysis, and all results are shown in
Figure 6.
Microbial communities gathered into two main clusters. One was mainly composed of Bacteroides, Rummeliibacillus, Enterococcus, Escherichia-Shigella, and Clostridium, which showed a significantly positive correlation with ciprostene and 7-aminoflunitrazepam-d7 and a significantly negative correlation with 7-methylguanine, S-adenosylhomocysteine, adenosine, and hippuric acid. The other cluster mainly consisted of Lactobacillus, Faecalibacterium, Ruminococcus, Bifidobacterium, and Fibrobacter, which showed a completely inverse correlation compared with the former, as evidenced by the significantly positive correlation with hippuric acid, carbamazepine-d10, diflucortolone pivalate, and diacerein, as well as the significantly negative correlation with DL-arginine, desloratadine, and propionylcarnitine. Specifically, Bifidobacterium, which was considered the probiotic bacteria, showed a significantly positive correlation with all up-regulated metabolites in LE treatment compared with HE. No other significant correlations were observed.
3.6. Liver Energy-Related Metabolic Gene Expression Responses to High Fiber Content Dietary Treatment
Differential analyses of the liver energy-related metabolic genes were carried out. A total of 293 significantly differentially expressed genes based on the filter standard of fold change > 2, and
p < 0.05 were identified, including 165 significant down-regulated and 129 significant up-regulated genes between HE and LE treatments. All these genes are list in
Table S3.
In addition, four energy metabolism-related genes, including two up-regulated and two down-regulated genes, were selected to verify the validity of the transcriptomic results (
Figure 7). As the results show, the up-regulated genes of the protein phosphatase, Mg
2+/Mn
2+ dependent 1 K (
PPM1K), and the protein phosphatase 1 inhibitor 3C (
PPR3C) in the transcriptomic results of the HE treatment were significantly higher compared with the LE treatment. Similarly, the down-regulated genes of GTPase activating Rap/RanGAP domain like 3 (
GARNL3) and growth arrest and DNA damage-inducible protein GADD45 beta (
GADD45B) also showed the same expressing alteration (
p < 0.05).
Finally, all differentially expressed genes were selected for the functional analysis, and the results are shown in
Figure 8. As results show, the up-regulated genes in high-fiber dietary-treated livers are mainly enriched in metabolism pathways, which include the fatty acid metabolism, PPAR signaling pathway, glutathione metabolism, and fatty acid degradation. Additionally, the down-regulated genes of low-energy-treated livers are mainly enriched into the circadian rhythm, glycolysis, and apelin signaling pathways. Specifically speaking, the highest number of genes (140) enriched into the metabolic pathway; however, not significant. No other significant pathways are enriched by the significantly differentially expressed genes.