Growth performance
Dietary DMG-Na supplementation improved the growth performance of IUGR piglets during the suckling period (Fig. 1 and Table 1). Throughout the trial, the body weights of all piglets increased as the piglets grew. The interaction effect (group × time × DMG-Na) was significant (P < 0.001). Throughout the trial, group I had significantly lower body weight (50.33%, 35.92%, 27.82%, 30.68%, 31.06%, 28.67%, 23.22%) than did group N. In addition, groups ND and ID had significantly higher body weight [(26.07%, 29.68%, 18.07%, 15.83%, 15.43%) and (12.50%, 15.80%, 16.78%, 23.00%, 18.75%)] than did groups N and I, respectively. From day 7 to day 21 during the suckling period, the interaction effect (group × DMG-Na) was significant for ADG (P = 0.009), ADFI (P = 0.017) and diarrhea (P < 0.001). From day 7 to day 21 during the suckling period, group I had significantly lower ADG (12.77%) and ADFI (2.50%), but higher diarrhea (201.12%) than did group N. In addition, groups ND and ID had significantly higher ADG (28.01%, 30.08%), but lower ADFI (8.70%, 7.86%) and diarrhea (10.11%,
Table 1
Dietary DMG-Na supplementation improved the growth performance of IUGR piglets from day 7 to day 21 during the suckling period.1
| Groups2 | P value |
Item3 | N | ND | I | ID | PG | PD | PG×D |
ADG (g/d) | 282 ± 1.28c | 361 ± 1.26a | 246 ± 1.21d | 320 ± 1.30b | < 0.001 | < 0.001 | 0.009 |
ADFI (g/d) | 494.73 ± 1.52a | 451.67 ± 1.55c | 482.35 ± 1.60b | 444.44 ± 1.28d | < 0.001 | < 0.001 | 0.017 |
G:F | 0.57 ± 0.02 | 0.80 ± 0.03 | 0.51 ± 0.01 | 0.72 ± 0.01 | < 0.001 | < 0.001 | 0.397 |
Diarrhea (%) | 0.89 ± 0.02b | 0.80 ± 0.01b | 2.68 ± 0.07a | 0.89 ± 0.01b | < 0.001 | < 0.001 | < 0.001 |
1 Data are expressed as the mean ± SD, n = 10. Different superscripts a, b, c, d (N, ND, I, ID group) represent significant differences (P < 0.05). |
2 NBW represents normal birth weight newborn piglets; IUGR represents intrauterine growth restriction newborn piglets; N represents NBW piglets fed a basic milk diet; ND represents NBW piglets fed a basic milk diet supplemented with 0.1% DMG-Na; I represents IUGR piglets fed a basic milk diet; ID represents IUGR piglets fed a basic milk diet supplemented with 0.1% DMG-Na. |
3 ADG, average daily gain; ADFI, average daily feed intake; G:F, gain: feed intake. |
Histological morphology assay
Dietary DMG-Na supplementation improved the histological morphology of the small intestine of IUGR piglets (Figs. 2, 3). The interaction effect (groups × intestines) was significant for villus width (P = 0.001) and villus area (P < 0.001) in trial A. However, the interaction effects in trials B (groups × time × intestines) and C (groups × DMG-Na × intestines) were not significant (P > 0.05) for any of the tested histological morphology traits (villus length, villus width, crypt depth, and villus area). In trial A, small intestinal villus length, crypt depth, villus width [(25.93%, P < 0.001), (32.73%, P < 0.001), (28.00%, P < 0.001)], and villus area [(35.90%, P < 0.001), (42.55%, P = 0.032), 34.17%] in IUGR newborn piglets were decreased relative to NBW newborn piglets. In trial B, small intestinal villus length, crypt depth, and villus area were increased in the N and I groups relative to the NBW and IUGR groups, respectively. In addition, small intestinal villus length, crypt depth, villus width, and villus area were decreased in the I group relative to the N group. In trial C, small intestinal villus length, crypt depth, villus width, and villus area were increased in the ND and ID groups relative to the N and I groups, respectively.
Redox status assay
Dietary DMG-Na supplementation improved the redox status of the small intestine of IUGR piglets (Fig. 4). The interaction effect (groups × intestines) was significant for GSH-Px (P = 0.046) and GSH (P = 0.016) in trial A. However, the interaction effects in trials B (groups × time × intestines) and C (groups × DMG-Na × intestines) were not significant (P > 0.05) for any of the tested redox status values (SOD, GSH-Px, GSH, CAT, and MDA). In trial A, decreases in small intestinal SOD, GSH-Px [(16.48%, P = 0.046), (0.43%, P = 0.008), (16.67%, P = 0.002)], GSH [(32.60%), (36.89%, P = 0.001), (31.62%, P = 0.024)], CAT, and an increase in small intestinal MDA were observed in IUGR newborn piglets relative to NBW newborn piglets. In trial B, increases in small intestinal SOD, GSH-Px, GSH, CAT, and MDA were observed in the N and I groups relative to the NBW and IUGR groups, respectively. In addition, a worse small intestinal redox status was found in the I group than in the N group. In trial C, an improvement in the small intestinal redox status was observed in the ND and ID groups relative to the N and I groups.
Mitochondrial redox status assay
Dietary DMG-Na supplementation improved the mitochondrial redox status of the small intestine of IUGR piglets (Fig. 5). The interaction effects were not significant (P > 0.05) in any trial for any of the tested mitochondrial redox status values (MnSOD, GSPx, GSH, and γ-GCL). In trial A, the mitochondrial redox status of the small intestine was worse in IUGR newborn piglets than in NBW newborn piglets. In trial B, the mitochondrial redox status of the small intestine was increased in the N and I groups relative to the NBW and IUGR groups. In addition, a worse small intestine mitochondrial redox status was found in the I group than in the N group. In trial C, the mitochondrial redox status of the small intestine was increased in the ND and ID groups relative to the N and I groups.
Oxidative damage assay
Dietary DMG-Na supplementation improved the oxidative damage values in the small intestines of IUGR piglets (Fig. 6). The interaction effects were not significant (P > 0.05) in any trial for any of the tested oxidative damage values (ROS, PC, 8-OHdG, MMP, apoptosis cells, necrotic cells, ATP, and mtDNA). In trial A, oxidative damage of the small intestine was enhanced in IUGR newborn piglets relative to NBW newborn piglets. In trial B, oxidative damage of the small intestine was increased in the N and I groups relative to the NBW and IUGR groups. In addition, the oxidative damage level in the small intestine was higher in the I group than in the N group. In trial C, small intestinal oxidative damage was improved in the ND and ID groups relative to the N and I groups.
RNA-seq and WGCNA assay
From the above results, we found that IUGR most affected the jejunum. Therefore, we identified the specific mechanism by which IUGR affected jejunum structure and function through RNA-seq and WGCNA methods.
To identify gene differentiation in the jejunum between the NBW and IUGR groups were constructed with total RNA and subjected to Illumina deep sequencing. Overviews of the sequencing and assembly results are presented in Table S2. After discarding the low-quality raw reads, 276,210,270 clean reads remained. All 19.935 assembled genes were referenced against the Swiss-Prot (12.630 genes; 63.36%), Nr (18.172 genes; 91.16%), Pfam (16.602 genes; 83.28%), KEGG (6.314 genes; 31.67%), KOG (10.976 genes; 55.06%), and GO (7.487 genes; 37.56%) databases (Table S3). To explore the molecular mechanisms of the jejunum in response to IUGR, RPKM analysis was performed to identify differentially expressed genes. We found that 14 mRNAs were upregulated (P < 0.05), while 25 genes were downregulated (P < 0.05) in response to IUGR (Table S4). We then performed hierarchical clustering of the differentially expressed genes based on the six samples’ log10 (RPKM + 1). The results indicated that the samples could be sorted into two distinct groups (Fig. 7A). Overall, IUGR treatment had a significant impact on the global gene expression profile of the jejunum in newborn piglets.
According to the GO classification system, genes involved in “cellular processes” (25 genes) and “metabolic processes” (21 genes) were notably represented in the biological process category. Among the cellular components, “cell” (21 genes) was the most commonly represented, followed by “cell part” (21 genes), and “organelle” (15 genes). In the category of molecular function, a significant proportion of clusters were assigned to “binding” (21 genes) and “catalytic activity” (15 genes) (Fig. S1). To classify orthologous gene products, the differentially expressed genes were subdivided into 25 KOG classifications. Among them, the cluster “General function prediction only” (10 genes) represented the largest group, followed by “Signal transduction mechanisms” (9 genes) and “Transcription” (4 genes) (Fig. S2). The Kyoto Encyclopedia of Genes and Genomes (KEGG) classification was found for the differentially expressed genes that were further classified into six biochemical pathways, including cellular processes (5 genes), environmental information processing (17 genes), genetic information processing (2 genes), human diseases (31 genes), metabolism (10 genes), and organismal systems (50 genes) (Fig. S3).
By performing the KEGG pathway analyses, a total of 15 pathways were identified that differed significantly (P < 0.05) between the IUGR newborn piglets and the NBW newborn piglets (Fig. 7B). Among these pathways, “Bile secretion,” “Pancreatic secretion,” and “Salivary secretion” were included in the “Digestive system” sub-class. “Regulation of lipolysis in adipocytes” and “PPAR signaling pathway” were included in the “Endocrine system” sub-class. In addition, some important subclasses have been significantly enriched, including “Signal transduction,” “Nervous system,” “Metabolism of cofactors and vitamins,” “Substance dependence,” “Membrane transport,” “Nervous system,” “Immune system,” “Folding, sorting and degradation,” “Environmental adaptation,” and “Endocrine system.” These results imply that the genes involved in these pathways may play crucial roles in the jejunum of newborn piglets in response to IUGR.
The WGCNA analysis found nine modules (power = 12) (Fig. S4). The 1,159 gene modules identified by WGCNA are shown in a cluster dendrogram (Fig. S5), in which the branches correspond to modules and each leaf in the branch represents one probe. The WGCNA R package allowed us to quantify the correlations between genes in each module and the relevant phenotypes, thus quantifying the module-trait associations. Figure 7C shows module-trait associations using a heatmap plot, which graphically represents Pearson’s correlation coefficients measured between every single module and trait (Table S5). The yellow module (P = 0.05, r = 0.81) was most significantly correlated with the trait. According to the GO classification system, the 254 genes in the yellow module were involved in “cellular process” (65 genes), and “single-organism process” (57 genes) were notably represented in the biological process category. Among the cellular components, “cell” (71 genes) and “cell part” (71 genes) were the most commonly represented, followed by “organelle” (50 genes). In the category of molecular function, a significant proportion of clusters were assigned to “binding” (73 genes) and “catalytic activity” (39 genes) (Fig. S6). According to the KEGG classification, the 254 genes of the yellow module were further classified into four biochemical pathways (Fig S7); cellular processes (25 genes), environmental information processing (59 genes), human diseases (58 genes), and organismal systems (118 genes). A total of 13 hub genes were obtained from the co-analysis of RNA-seq and WGCNA (Table S6).
Gene expression
Dietary DMG-Na supplementation improved the jejunum redox status-related, cell adhesion-related, and mitochondrial function-related gene expression in IUGR piglets (Fig. 8). There was a significant difference (P < 0.05) between the NBW and IUGR groups in trial A in redox status-related gene expression (Nrf2, HO1, SOD1, GSH-Px, Sirt1, PGC1α, SOD2, γ-GCL, Trx2, Trx-R2, Prx3), cell adhesion-related gene expression (OCLN, CLDN2, CLDN3, and ZO1), and mitochondrial function-related gene expression (MCD, MCAD, SDH, UCP2, COX2, CS, COX1, Cyt C, ATP8, MHC1, mtTFA, Ndufa2, NRF1, UCP1, POLG1, POLG2, SSBP1, Drp1, Fis1, and Mfn2). The interaction effect (groups × time) was significant for Nrf2 (P = 0.035), SOD1 (P = 0.032), PGC1α (P = 0.028), SIRT1 (P = 0.009), SOD2 (P = 0.026), Trx-R2 (P = 0.008), Prx3 (P = 0.010), OCLN (P = 0.005), MCD (P = 0.002), MCAD (P = 0.014), SDH (P = 0.010), UCP2 (P = 0.001), COX2 (P = 0.010), CS (P = 0.011), COX1 (P = 0.018), Cyt C (P = 0.038), ATP8 (P = 0.016), MHC1 (P = 0.046), mtTFA (P = 0.039), Ndufa2 (P = 0.005), NRF1 (P = 0.014), UCP1 (P = 0.023), POLG1 (P = 0.009), POLG2 (P = 0.018), SSBP1 (P = 0.007), Drp1 (P = 0.010), Fis1 (P = 0.009), and Mfn2 (P = 0.012) in trial B. In addition, the interaction effect (groups × time × intestines) was significant for Prx3 (P = 0.028), SDH (P = 0.031), UCP2 (P = 0.008), COX2 (P = 0.035), mtTFA (P = 0.036), and Mfn2 (P = 0.049) in trial C. In trial B, better values of redox status-related gene expression [among which, Nrf2 (39.00%, 43.26%), SOD1 (43.00%, 36.36%), PGC1α (43.00%, 33.93%), SIRT1 (43.00%, 31.58%), SOD2 (45.00%, 40.00%), Trx-R2 (46.00%, 30.36%), and Prx3 (46.00%, 32.76%) were significant], cell adhesion-related gene expression [among which, OCLN (43.00%, 22.41%) was significant], and mitochondrial function-related gene expression [among which, MCD (49.00%, 32.73%), MCAD (44.00%, 26.92%), SDH (47.00%, 29.41%), UCP2 (46.00%, 24.53%), COX2 (48.00%, 40.35%), CS (44.00%, 27.12%), COX1 (48.00%, 37.10%), Cyt C (48.00%, 39.22%), ATP8 (45.00%, 34.55%), MHC1 (42.00%, 33.93%), mtTFA (49.00%, 45.28%), Ndufa2 (48.00%, 21.05%), NRF1 (48.00%, 28.85%), UCP1 (46.00%, 26.79%), POLG1 (46.00%, 27.12%), POLG2 (45.00%, 31.37%), SSBP1 (48.00%, 29.31%), Drp1 (45.00%, 25.00%), Fis1 (44.00%, 28.57%), and Mfn2 (49.00%, 31.03%) were significant] were observed in N and I groups than in NBW and IUGR groups. In trial C, an improvement in redox status-related gene expression [among which, Prx3 (29.45%, 23.38%) was significant], cell adhesion-related gene expression, and mitochondrial function-related gene expression [among which, SDH (23.81%, 16.67%), UCP2 (27.40%, 19.70%), COX2 (27.70%, 18.75%), mtTFA (26.85%, 19.48%) were significant] was observed in the ND and ID groups compared to the N and I groups.
Protein expression
Dietary DMG-Na supplementation improved the jejunum redox status-related, cell adhesion-related, and mitochondrial function-related protein expression in IUGR piglets (Fig. 9). There was a significant difference (P < 0.05) between the NBW and IUGR groups in trial A in redox status-related protein expression (Nrf2, HO1, SOD, GSH-Px, Sirt1, and PGC1α), cell adhesion-related protein expression (OCLN and ZO1), and mitochondrial function-related protein expression (Cyt C, mtTFA, Mfn2, Drp1, and Fis1). The interaction effect (groups × time) was significant for Nrf2 (P = 0.002), SOD (P < 0.001), GSH-Px (P = 0.014), Sirt1 (P = 0.004), PGC1α (P = 0.004), OCLN (P < 0.001), ZO1 (P < 0.001), Cyt C (P = 0.011), mtTFA (P = 0.041), Mfn2 (P = 0.037), Drp1 (P = 0.001), and Fis1 (P = 0.004) in trial B. In addition, the interaction effect (groups × time × intestines) was significant for mtTFA (P = 0.007) in trial C. In trial B, better values of redox status-related protein expression [among which, Nrf2 (42.00%, 69.92%), SOD (59.00%, 20.97%), GSH-Px (51.00%, 48.15%), Sirt1 (62.00%, 52.94%), and PGC1α (63.00%, 46.15%) were significant], cell adhesion-related protein expression [among which, OCLN (59.00%, 31.48%) and ZO1 (58.00%, 39.22%) was significant], and mitochondrial function-related protein expression [among which, Cyt C (48.00%, 39.22%), mtTFA (35.00%, 33.96%), Mfn2 (52.00%, 46.15%), Drp1 (54.00%, 36.36%), and Fis1 (53.00%, 41.18%) were significant] were observed in N and I groups than in NBW and IUGR groups. In trial C, an improvement in redox status-related protein expression, cell adhesion-related protein expression, and mitochondrial function-related protein expression [among which, mtTFA (36.30%, 29.58%) was significant] was observed in the ND and ID groups compared to the N and I groups.