This animal study was approved by the Institutional Animal Care and Use Committee of Qinghai University, China (Approval No.: QUA-2020-0710). Both the animal care and experimental protocols were reviewed and approved by the Committee. The study was conducted in compliance with local laws and regulations as well as institutional requirements.
Experimental Design
A feeding experiment was conducted at a Tibetan sheep breeding base in Qieji Township, Gonghe County, Qinghai Province (36°19′N, 99°41′E) from April 2024 to July 2024. Sixty 2-month-old plateau-type Tibetan sheep rams with similar body weights (17.72 ± 0.19 kg) and good body condition were selected and randomly divided into 2 groups: the control group (Ctrl) and the experimental group (Bet), with 30 sheep in each group (5 replicates, 6 sheep per replicate). The Ctrl group was fed a basal diet, while the Bet group was fed the basal diet supplemented with 0.08% betaine (on an air-dry basis). The basal diet consisted of a concentrate supplement and roughage at a concentrate-to-roughage ratio of 7:3. The roughage was a mixture of oat hay and oat silage in a 1:1 ratio on a dry matter basis. The betaine was provided by Jinan Tai Fei Animal Husbandry Technology Co., Ltd. with a purity of 98%. The composition and nutritional levels of the experimental diets are shown in Table 1. The total experimental period was 100 days, including a 10-day pre-test period and a 90-day formal test period. Feeding was carried out at 8:00 and 17:00 daily, and the sheep had free access to feed and water during the experiment. The nutritional components in the feed used in this experiment, such as dry matter (Thiex et al., 2002a), crude protein (Thiex et al., 2002b), acid detergent fiber (HIfiFLHQF et al.), neutral detergent fiber (Hiraoka et al., 2012), Ca and P (Babos et al., 2018), were all determined with reference to relevant standards.
Table 1
Composition of the basic diet.
Item | Ctrl | Bet |
|---|
Corn | 48.40 | 48.32 |
Wheat | 9.00 | 9.00 |
Palm meal | 16.00 | 16.00 |
Soybean meal | 4.00 | 4.00 |
Rapeseed meal | 15.00 | 15.00 |
Nacl | 1.00 | 1.00 |
Limestone | 1.00 | 1.00 |
Baking soda | 1.00 | 1.00 |
1% Premix1 | 0.60 | 0.60 |
4% Concentrate | 4.00 | 4.00 |
Betaine | 0.00 | 0.08 |
Total | 100.00 | 100.00 |
Nutrient levels2 | | |
Digestible energy (MJ·kg⁻¹) | 12.66 | 12.65 |
Crude protein | 13.69 | 13.69 |
Ether extract | 4.61 | 4.61 |
Crude fiber | 21.18 | 21.16 |
Neutral detergent fiber | 21.18 | 21.18 |
Acid detergent fiber | 14.72 | 14.72 |
Ca | 0.85 | 0.85 |
P | 0.41 | 0.41 |
Ctrl, basal diet; Bet, 0.08% betaine is added to the basal diet.
1The premix provides the following per kilogram of the feed: Cu, 18 mg; Fe, 66 mg; Zn, 30 mg; Mn, 48 mg; Se, 0.36 mg; I, 0.6 mg; Co, 0.24 mg; vitamin A, 24,000 IU; vitamin D, 4,800 IU; vitamin E, 48 IU.
2Digestible energy is a calculated value, while the others are measured values.
Sample Collection
At the end of the formal test period, 6 Tibetan sheep were randomly selected from each group and slaughtered at the slaughterhouse (n = 12). The Tibetan sheep were fasted for 12 hours. After slaughter, the intestines were quickly separated, and the ileocecal junction was clamped with hemostatic forceps. The collected ileal tissue samples were about 10 cm away from the ileocecal orifice, with a length of about 1 cm. Meanwhile, 4 mL of ileal contents were collected into a 5 mL cryopreservation tube, quickly frozen in liquid nitrogen, stored and transported in dry ice to the laboratory, and then preserved in a -80℃ refrigerator for subsequent untargeted metabolomics and 16S rDNA sequencing. Finally, the collected tissue samples were rinsed clean with physiological saline and placed in 4% paraformaldehyde for tissue observation.
H&E staining
The ileal tissues were fixed in 4% paraformaldehyde solution for no less than 48 hours, followed by gradient dehydration with ethanol and paraffin embedding. Finally, sections with a thickness of 3 µm were cut to prepare paraffin sections, which were then stained with hematoxylin-eosin (H&E) for analysis. The indicators measured included villus height (VH), villus width, crypt depth (CD), crypt number, muscular layer thickness, and the VH/CD ratio. A microscope (OLYMPUS, DP26, Tokyo, Japan) was used to observe the stained sections and collect images. Finally, Image-Pro Plus 6.0 software (Media Cybernetics, Inc., Rockville, MD, USA) was employed to measure all indicators of the ileal tissues in the collected images at magnification levels of 500× and 200×.
Determination of antioxidant capacity and lipopolysaccharide content
The contents of total antioxidant capacity (T-AOC), superoxide dismutase (SOD), glutathione peroxidase (GSH-Px), catalase (CAT), malondialdehyde (MDA) and the level of lipopolysaccharide (LPS) in ileal tissues were determined using ELISA kits (Shanghai Meilian Biotechnology Co., Ltd., Shanghai, China). Take 0.1 g of ileal tissue and add 900 µL of physiological saline. After crushing with a crusher (SM2000, Retsch, Haan, Germany), transfer to a low-temperature high-speed centrifuge (5430 R, Eppendorf, Hamburg, Germany) and centrifuge at 2500 × g for 15 minutes at 4°C to obtain the supernatant. The determination was carried out strictly in accordance with the steps provided in the kit, and finally the absorbance value was measured at a wavelength of 450 nm on an ELISA Analyzer (ELx808, BioTek, Winooski, VT, USA).
Determination of SCFAs
The ileal contents were preserved in dry ice and sent to Beijing Novogene Bioinformatics Technology Co., Ltd. (Beijing, China). The composition and content of SCFAs were quantitatively determined using an ultra-high performance liquid chromatography-tandem mass spectrometry system (Vanquish™ Flex UHPLC-TSQ Altis™, Thermo Scientific Corp., Germany). The sample extraction and determination procedures were carried out with reference to the method described by Yan et al. (Yan et al., 2022).
Real Time Quantitative (RT-qPCR)
Total RNA was extracted from ileal tissues using the Transzol Up kit (TRAN, Beijing, China), and its concentration was detected. Subsequently, the RNA was reverse-transcribed into cDNA using a kit. The subsequent fluorescent quantitative steps and qPCR reaction procedures were strictly performed according to the method described by Ji et al (Ji et al., 2024). The relative expression level of the target gene mRNA was calculated by the 2−ΔΔCt method, and the primer sequences are shown in Table 2.
Table 2
Items1 | Primer sequence (5’-3’)2 | Tm (℃) | Product length | NCBI Gene ID |
|---|
Claudin-1 | F-CCTGCTGTGCTGCTCCTGTC | 61.6 | 75bp | 780473 |
R-GAAGGTGCTGGCTTGGGATAGG | 61.4 |
Occludin | F-CGAGAAGCGACCGTATCCAGAG | 61.4 | 129bp | 101118304 |
R-TCCAAGTTACCACTGCTGCTGTAG | 59.6 |
ZO-1 | F-GGGCAAGTTAAAGATGGTGGTTCAG | 59.6 | 93bp | 443200 |
R-GAGGCGTCAGCAGAGTGGATG | 61.5 |
TNF-α | F-ACGGCGTGGAGCTGAAAGAC | 59.5 | 79bp | 443540 |
R-CTGAAGAGGACCTGCGAGTAGATG | 61.3 |
Muc-2 | F-ACGACTCCTACGCCCTCCTG | 61.6 | 130bp | 780488 |
R-ACGCTGCCATCCGACTTGAAG | 51.5 |
IL-6 | F-TCTAATAACCACTCCAGCCACACAC | 59.6 | 77bp | 443406 |
R-TTGCGTTCTTTACCCACTCGTTTG | 57.9 |
β-Actin | F-AGCAAGCGTGGCATCCTAACC | 59.5 | 87bp | 443340 |
R-ATCTTCTCCATGTCGTCCCAGTTG | 59.6 | 77bp |
1ZO-1, zonula occluden-1; MUC-2, mucoprotein-2; TNF-α, tumor necrosis factor-α; IL-6, interleukin-6.
2F, forward primer; R, reverse primer.
Analysis of microbial composition
Microbial DNA in ileal contents was extracted using the DNA extraction kit from Novogene (Beijing, China), and the V3-V4 region of 16S rDNA was amplified by PCR with specific primers 341F (5′-CCTAYGGGRBGCASCAG-3′) and 806R (5′-GGACTACNNGGGTATCTAAT-3′). Sample reads were spliced using FLASH (Version 1.2.11) to generate Raw Tags data. Then, Cutadapt (Version 3.3) software was used to match the reverse primer sequences and remove redundant parts, and fastp (Version 0.23.1) software was employed to obtain Clean Tags (Bokulich et al., 2013). Subsequently, the Tags sequences were compared with the species annotation database, and chimeric sequences were detected and removed to obtain Effective Tags (Edgar et al., 2011). The SVG function of Perl (Version 5.26.2) software was used to draw the relative abundance distribution histogram. QIIME2 (Version 2022.02) software was used to calculate the Alpha diversity index. Principal coordinate analysis was performed using the ade4 package and ggplot2 package of R (Version 4.0.3) software.
Untargeted metabolomics analysis
The extraction of metabolites referred to the study by Want et al. (Want et al., 2013). The chromatographic column used was Hypersil Gold C18. The off-machine data were processed using Compound Discoverer 3.3 software(Thermo Fisher Scientific Inc., Waltham, Massachusetts, USA), with data processing based on the Linux system (CentOS version 6.6) as well as R (Version 4.0.3) and Python (Version 3.12.1) software. The identified metabolites were annotated using the KEGG database (https://www.genome.jp/kegg/pathway.html). For multivariate statistical analysis, after data transformation using metaX software, principal component analysis and orthogonal partial least squares discriminant analysis were performed to obtain VIP values; univariate analysis calculated P values and FC values based on t-tests, and the criteria for screening differential metabolites were VIP > 1, P < 0.05, and FC > 1.5. The R package ggplot2 was used to draw match plots and bubble plots, the corrplot package was used for correlation analysis, and the KEGG database was used to study metabolite functions and pathways. The enrichment of metabolic pathways was determined according to x/n > y/n and P < 0.05.
Statistical analysis
Experimental data were preprocessed using Microsoft Excel 2024, and all data were subjected to independent samples t-test using SPSS 26.0 software (IBM Corp., Armonk, NY, USA). Results were expressed as mean ± standard error. Graphs were plotted using GraphPad Prism 8.0 software, with P < 0.05 set as the criterion for judging significant differences with statistical significance. Wilcoxon rank-sum test was used for data α-diversity index analysis, Unweighted Unifrac distance was used for PCoA analysis, and NMDS was used to reflect the degree of differences between samples.
Spearman rank correlation analysis was performed using the R language Psych package and R language Vegan package in R (Version 4.0.3) to calculate the correlations between ileal microbiota abundance, metabolite levels and short-chain fatty acids, antioxidant indices, tissue morphological parameters, ileal barrier function, and lipopolysaccharide levels.