Correlation analysis of intestinal ATPase activity and antioxidant indexes with growth traits and feed conversion efficiency of Hu sheep
The correlations between growth traits, feed intake, and feed efficiency in Hu sheep with the ATPase activity and antioxidant indicators in the jejunal tissue are presented in Fig. 1A. The activities of Na+ K+-ATPase, Ca2+ Mg2+-ATPase, and SOD showed significant negative correlations with birth weight (P < 0.05). Na+ K+-ATPase and Ca2+ Mg2+-ATPase activities correlated significantly and negatively with the daily weight gain of Hu sheep at 120–140 days old (P < 0.05). SOD activity exhibited a significant positive correlation with the daily weight gain at 0–80 days old (P < 0.05). The MDA content demonstrated an extremely significant positive correlation with the daily weight gain of Hu sheep at 0–80 days old. Furthermore, from 80 days old until the end of the 180-day trial, the MDA content correlated significantly and positively with the body weights at various stages (P < 0.05). In terms of feed efficiency (Fig. 1B), the jejunal tissue total antioxidant capacity (T-AOC) showed a significant negative correlation with the residual feed intake at 160–180 days old (P < 0.05). The Na+ K+‑ATPase activity exhibited a significant positive correlation with feed conversion efficiency at 120–140 days old (P < 0.05). The MDA content showed a significant positive correlation with the average daily feed intake at 80–100 days old (P < 0.05).
Effects of intestinal ATPase activity and MDA content on growth traits and feed conversion efficiency of Hu sheep
Based on the results of the correlation analysis, the relationship between the MDA content and antioxidant indicators and the growth and feeding traits of Hu sheep was found to be significant. Na+ K+-ATPase is considered the most crucial ATPase on the cell membrane. Therefore, in this experiment, the top 10 individuals with the highest and lowest ATPase activity and MDA content were selected, forming the H-ATPase and L-ATPase groups, H-MDA and L-MDA groups. The differences in growth traits and feed conversion efficiencies between these extreme groups were compared (Table 1). The results showed that the birth weight of Hu sheep in the L-ATP group was significantly higher than that in the H-ATP group (P < 0.05). At 80 d, 100 d, and 120 d, the body weight of the H-MDA group was significantly higher than that of the L-MDA group (P < 0.05). From 0–80 d, the ADG of the H-MDA group was significantly higher than that of the L-MDA group (P < 0.05). At 120–140 d, the FCR of the H-MDA group was significantly higher than that of the L-MDA group (P < 0.05).
Table 1 Effects of high and low ATPase activity and MDA content on growth performance and feed efficiency of Hu sheep.
Items
|
Groups
|
SEM
|
P-value
|
Groups
|
SEM
|
P-value
|
H-ATP
|
L-ATP
|
H-MDA
|
L-MDA
|
BW, Kg
|
Birth
|
3.59
|
4.17
|
0.22
|
0.017
|
4.30
|
4.00
|
0.28
|
0.313
|
80 d
|
17.99
|
17.80
|
2.07
|
0.926
|
20.88
|
16.84
|
1.62
|
0.023
|
100 d
|
22.76
|
23.10
|
2.21
|
0.881
|
26.24
|
22.29
|
1.83
|
0.044
|
120 d
|
28.58
|
29.32
|
2.44
|
0.766
|
32.04
|
27.63
|
1.92
|
0.033
|
140 d
|
33.99
|
35.32
|
2.53
|
0.607
|
37.36
|
33.36
|
2.13
|
0.077
|
160 d
|
40.03
|
41.12
|
2.52
|
0.669
|
43.37
|
39.38
|
2.19
|
0.086
|
180 d
|
45.83
|
46.42
|
2.52
|
0.817
|
48.55
|
44.67
|
2.43
|
0.128
|
ADG, Kg/d
|
0–80 d
|
0.18
|
0.17
|
0.03
|
0.721
|
0.21
|
0.16
|
0.02
|
0.022
|
80–100 d
|
0.24
|
0.27
|
0.03
|
0.321
|
0.27
|
0.27
|
0.03
|
0.880
|
120–140 d
|
0.27
|
0.30
|
0.02
|
0.148
|
0.27
|
0.29
|
0.03
|
0.446
|
160–180 d
|
0.29
|
0.27
|
0.03
|
0.419
|
0.26
|
0.27
|
0.03
|
0.900
|
80–180 d
|
0.28
|
0.29
|
0.02
|
0.594
|
0.28
|
0.28
|
0.01
|
0.945
|
ADFI, Kg/d
|
80–100 d
|
1.01
|
1.09
|
0.12
|
0.486
|
1.17
|
1.06
|
0.10
|
0.284
|
120–140 d
|
1.72
|
1.68
|
0.13
|
0.789
|
1.84
|
1.69
|
0.13
|
0.242
|
160–180 d
|
2.02
|
1.92
|
0.13
|
0.446
|
2.04
|
2.00
|
0.13
|
0.752
|
80–180 d
|
1.60
|
1.63
|
0.10
|
0.765
|
1.70
|
1.60
|
0.10
|
0.366
|
FCR
|
80–100 d
|
4.21
|
4.24
|
0.44
|
0.957
|
4.41
|
3.94
|
0.37
|
0.213
|
120–140 d
|
6.46
|
5.68
|
0.55
|
0.171
|
7.40
|
5.91
|
0.68
|
0.042
|
160–180 d
|
7.10
|
7.73
|
0.85
|
0.468
|
8.03
|
8.08
|
0.83
|
0.954
|
80–180 d
|
5.73
|
5.73
|
0.33
|
0.995
|
6.17
|
5.75
|
0.28
|
0.150
|
RFI, Kg/d
|
80–100 d
|
-0.04
|
0.01
|
0.05
|
0.379
|
-0.06
|
0.00
|
0.05
|
0.236
|
120–140 d
|
0.04
|
-0.08
|
0.06
|
0.061
|
0.02
|
0.02
|
0.06
|
0.935
|
160–180 d
|
0.01
|
-0.10
|
0.12
|
0.366
|
-0.04
|
0.05
|
0.07
|
0.219
|
80–180 d
|
-0.02
|
-0.02
|
0.05
|
1.000
|
-0.02
|
0.02
|
0.04
|
0.315
|
Abbreviations: BW = body weight, ADG = average daily gain, ADFI = average daily feed intake, FCR = feed conversion ratio, RFI = residual feed intake, H-ATP = high ATP level group, L-ATP = low ATP level group, H-MDA = high MDA level group, L- MDA = low MDA level group, SEM = standard error of mean.
Correlation analysis of ATPase activity and antioxidant indexes in intestinal tissues with the development of internal organs of Hu sheep
Correlation analysis was conducted between the ATPase activity and antioxidant indicators in intestinal tissues and the visceral organ weights of Hu sheep (Fig. 2). The results indicated a significant positive correlation between the Na+ K+-ATPase, Ca2+ Mg2+-ATPase, and SOD activities in intestinal tissues and the weight of the Hu sheep heart (P < 0.05). Additionally, there was a significant positive correlation between MDA contents and the Hu sheep heart weight, lung weight, and cecum weight (P < 0.05).
Effects of intestinal ATPase activity and MDA content on the intestinal tract development of Hu sheep
We conducted an analysis to investigate the impact of intestinal ATPase activity and MDA content on the development of the intestinal tract in Hu sheep (Table 2). The results indicate that in the H-ATP group, the length of the ileum in Hu sheep is significantly higher than that in the L-ATP group (P < 0.05), with no significant impact observed on other indicators due to ATP activity (P > 0.05). However, in the H-MDA group, the relative duodenum weight (%body), jejunum weight, relative jejunum weight (%body), relative jejunum weight (%intestinal tract), ileum weight, and relative colon weight (%intestinal tract) are all significantly higher compared to the L-MDA group (P < 0.05).
Table 2 Effects of high and low ATPase activity and MDA content on development of intestinal tract of Hu sheep
Traits
|
Groups
|
SEM
|
P-value
|
Groups
|
SEM
|
P-value
|
H-ATP
|
L-ATP
|
H-MDA
|
L-MDA
|
Duodenum
|
Weight (g)
|
35.35
|
36.20
|
1.779
|
0.739
|
40.22
|
36.95
|
2.598
|
0.386
|
Length (cm)
|
65.20
|
64.50
|
3.136
|
0.876
|
64.05
|
64.80
|
2.458
|
0.832
|
Relative weight (%body)
|
0.08
|
0.07
|
0.000
|
0.895
|
0.09
|
0.07
|
0.000
|
0.035
|
Relative length (%IT)
|
1.66
|
1.79
|
0.001
|
0.273
|
1.78
|
1.79
|
0.001
|
0.911
|
Jejunum
|
Weight (g)
|
899.80
|
877.10
|
30.360
|
0.603
|
931.36
|
853.85
|
21.743
|
0.021
|
Length (m)
|
29.36
|
27.94
|
0.779
|
0.213
|
29.51
|
28.03
|
0.780
|
0.195
|
Relative weight (%body)
|
1.86
|
1.79
|
0.001
|
0.538
|
1.89
|
1.67
|
0.001
|
0.050
|
Relative length (%IT)
|
76.25
|
77.49
|
0.006
|
0.192
|
78.13
|
76.73
|
0.007
|
0.174
|
Ileum
|
weight (g)
|
29.80
|
24.70
|
2.187
|
0.116
|
28.75
|
23.20
|
1.664
|
0.030
|
length (cm)
|
51.70
|
40.60
|
3.684
|
0.047
|
47.40
|
44.20
|
4.630
|
0.631
|
Relative weight (%body)
|
0.06
|
0.05
|
0.000
|
0.220
|
0.05
|
0.05
|
0.000
|
0.281
|
Relative weight (%IT)
|
2.02
|
1.77
|
0.001
|
0.162
|
2.03
|
1.72
|
0.001
|
0.053
|
Colon
|
weight (g)
|
385.05
|
409.35
|
13.136
|
0.207
|
383.80
|
360.95
|
19.666
|
0.422
|
length (cm)
|
742.90
|
670.10
|
33.237
|
0.139
|
691.30
|
679.00
|
39.965
|
0.830
|
Relative weight (%body)
|
0.80
|
0.84
|
0.000
|
0.432
|
0.74
|
0.78
|
0.000
|
0.545
|
Relative length (%IT)
|
20.52
|
18.60
|
0.007
|
0.081
|
19.66
|
18.00
|
0.008
|
0.144
|
Cecum
|
weight (g)
|
61.05
|
56.85
|
3.146
|
0.358
|
61.70
|
54.10
|
4.549
|
0.253
|
length (cm)
|
39.90
|
35.40
|
2.389
|
0.199
|
37.70
|
34.50
|
2.852
|
0.438
|
Relative weight (%body)
|
0.13
|
0.12
|
0.000
|
0.305
|
0.12
|
0.12
|
0.000
|
0.735
|
Relative length (%IT)
|
1.04
|
0.98
|
0.001
|
0.506
|
1.06
|
0.91
|
0.001
|
0.115
|
The total weight of IT (kg)
|
1.46
|
1.40
|
0.045
|
0.406
|
1.42
|
1.31
|
0.038
|
0.052
|
The total length of IT (m)
|
38.76
|
36.05
|
0.735
|
0.017
|
37.16
|
36.00
|
1.079
|
0.457
|
Abbreviations: IT = intestinal tract, H-ATP = high ATP level group, L-ATP = low ATP level group, H-MDA = high MDA level group, L- MDA = low MDA level group, SEM = standard error of mean.
Effects of ATPase activity and MDA content on the intestinal morphology of Hu sheep
We conducted a further analysis of the differences in intestinal tissue morphological indicators between the H-ATP group and L-ATP group, as well as between the H-MDA group and L-MDA group (Table 3). The results indicated that the crypt depth in the H-ATP group was significantly higher than that in the L-ATP group (P < 0.05), and the circular muscle thickness was significantly higher than that in the L-ATP group (P < 0.05). Additionally, in the H-MDA group, the villus height, crypt depth, and longitudinal muscle thickness were significantly higher than those in the L-MDA group. There were no significant differences observed for the other indicators.
Table 3 Effects of high and low ATP and MDA on jejunum tissue morphology of Hu sheep
Items
|
Groups
|
SEM
|
P-value
|
Groups
|
SEM
|
P-value
|
H-ATP
|
L-ATP
|
H-MDA
|
L-MDA
|
VH, µm
|
813.58
|
807.23
|
27.46
|
0.818
|
910.78
|
733.52
|
36.09
|
0.000
|
VW, µm
|
288.01
|
260.45
|
18.85
|
0.149
|
272.42
|
269.42
|
22.71
|
0.896
|
CD, µm
|
577.13
|
530.66
|
18.01
|
0.012
|
589.31
|
472.99
|
18.80
|
0.000
|
AMT, µm
|
285.57
|
223.00
|
13.23
|
0.000
|
276.64
|
267.10
|
20.09
|
0.636
|
LMT, µm
|
86.08
|
76.91
|
5.51
|
0.101
|
102.00
|
86.46
|
4.89
|
0.002
|
VH:CD, %
|
142.94
|
153.33
|
5.70
|
0.073
|
155.67
|
158.22
|
7.44
|
0.732
|
VW:CD, %
|
50.20
|
50.02
|
3.65
|
0.961
|
47.30
|
58.37
|
5.10
|
0.034
|
Abbreviations: VH = Villus height, VW = Villus width, CD = Crypt depth, AMT = Annular muscle thickness, LMT = Longitudinal muscle thickness, VH:CD = Villus height: Crypt depth, VW:CD = Villus width: Crypt depth, H-ATP = high ATP level group, L-ATP = low ATP level group, SEM = standard error of mean, H-MDA = high MDA level group, L- MDA = low MDA level group.
Effects of ATPase activity and MDA content on intestinal microbial diversity of Hu sheep
To investigate the impact of intestinal tissue ATPase activity and antioxidant function on the diversity of the microbial community in Hu sheep, this experiment employed 16S rRNA amplicon sequencing technology to compare the differences in the gut microbiota between the high H-ATP and low L-ATP groups, as well as between the high H-MDA and low L-MDA groups. After filtering and quality control of the raw data, 74,568 to 141,026 effective sequences were obtained, with Q20 and Q30 high quality data exceeding 97.41% and 91.75%, respectively. The statistical results from the data processing are presented in Supplementary Tables S2 and S3. Dilution curves indicated ample sequencing depth, meeting the requirements for subsequent data analysis (Supplementary Fig. S1). As shown in Fig. 3A, 8711 ASVs were detected collectively in the H-ATP and L-ATP groups, with only 1808 ASVs shared between the two groups, and the unique ASVs in each group were 2681 and 4222, respectively. Similarly, Fig. 3B shows that 8191 ASVs were detected collectively in the H-MDA and L-MDA groups, with only 1532 ASVs shared between the two groups. The H-MDA group had 2458 unique ASVs, while the L-MDA group had 4201 unique ASVs.
The results of alpha diversity analysis (Table 4) indicated a significant increase in both observed features and Chao1 index in the L-ATP group compared with those in the H-ATP group (P < 0.05). Furthermore, a noticeable upward trend was observed in the L-MDA group concerning observed features and Chao1 index compared with that in the H-MDA group, although this difference fell just outside the conventional significance threshold (0.05 < P < 0.1). However, PCoA analysis revealed no distinct clustering based on weighted unifrac measurements (Fig. 4A-B). Additionally, ANOSIM analysis using the Bray–Curtis metric demonstrate an non-significant difference between the groups (P > 0.05).
Table 4 Abundance and diversity index of jejunal fecal microflora of Hu sheep
Diversity indices
|
Groups
|
SEM
|
P-Value
|
Groups
|
SEM
|
P-Value
|
H-ATP
|
L-ATP
|
H-MDA
|
L-MDA
|
Chao1 index
|
680.06
|
928.76
|
113.36
|
0.042
|
650.85
|
902.65
|
128.86
|
0.066
|
Observed index
|
657.50
|
898.20
|
112.34
|
0.046
|
621.70
|
863.00
|
124.71
|
0.069
|
Shannon index
|
5.05
|
5.65
|
0.77
|
0.442
|
4.83
|
5.70
|
0.65
|
0.193
|
Simpson index
|
0.82
|
0.87
|
0.08
|
0.572
|
0.81
|
0.90
|
0.07
|
0.254
|
Abbreviations: H-ATP = high ATP level group, L-ATP = low ATP level group, SEM = standard error of mean, H-MDA = high MDA level group, L- MDA = low MDA level group.
Effects of ATPase activity and MDA content on the intestinal microbial composition of Hu sheep
At the phylum level, the predominant microbial taxa in the jejunum (relative abundance > 5%) were Firmicutes, Proteobacteria, Actinobacteriota, and Euryarchaeota in each group (Fig. 5A-B). Among the top 10 most abundant phyla, the relative abundance of Euryarchaeota in the H-MDA group was significantly lower than that in the L-MDA group (P < 0.05), while there were no significant differences in the relative abundance of the other phyla (P > 0.05). At the genus level (Fig. 5C-D), the relative abundance of Olsenella in the H-ATP group was significantly higher than that in the L-ATP group (P < 0.05), while Eubacterium hallii group and Blautia showed significantly lower relative abundance in the H-ATP group (P < 0.05). In the H-MDA group, the relative abundances of Methanobrevibacter and Clostridia UCG-014 were significantly lower than those in the L-MDA group (P < 0.05). At the species level (Fig. 5E), the relative abundance of Olsenella umbonata in the H-MDA group was significantly higher than that in the L-MDA group (P < 0.05), while the abundance of Methanobrevibacter ruminantium was significantly lower than that in the L-MDA group (P < 0.05).
In vitro validation of the interaction between Olsenella umbonata and MDA
To directly validate the interactions between Olsenella umbonata and MDA, the anaerobic growth of Olsenella umbonata and two other common intestinal bacteria was characterized at different concentrations of MDA. The results indicated that with increasing MDA concentrations, the OD600 of Olsenella umbonata was higher after 4 hours of cultivation, demonstrating that MDA indeed promoted the proliferation of Olsenella umbonata (Fig. 6A). By contrast, the anaerobic cultivation results for Selenomonas bovis and Acidaminococcus intestini (Fig. 6B-C) showed inconsistent effects of different MDA concentrations on the growth curves of these two species, with specific MDA concentrations inhibiting their proliferation (P < 0.05).