Comparison of the gut microbiota and metabolites between Diannan small ear pigs and Diqing Tibetan pigs

Objective Host genetics and environment participate in the shaping of gut microbiota. Diannan small ear pigs and Diqing Tibetan pigs are excellent native pig breeds in China and live in different environments. However, the gut microbiota of Diannan small ear pigs and Diqing Tibetan pigs were still rarely understood. Therefore, this study aimed to analyze the composition characteristics of gut microbiota and metabolites in Diannan small ear pigs and Diqing Tibetan pigs. Methods Fresh feces of 6 pigs were randomly collected from 20 4-month-old Diannan small ear pigs (DA group) and 20 4-month-old Diqing Tibetan pigs (TA group) for high-throughput 16S rRNA sequencing and liquid chromatography-mass spectrometry (LC-MS) non-targeted metabolome analysis. Results The results revealed that Firmicutes and Bacteroidetes were the dominant phyla in the two groups. Chao1 and ACE indices differed substantially between DA and TA groups. Compared with the DA group, the relative abundance of Prevotellaceae, and Ruminococcus was significantly enriched in the TA group, while the relative abundance of Lachnospiraceae, Actinomyces, and Butyricicoccus was significantly reduced. Cholecalciferol, 5-dehydroepisterol, stigmasterol, adrenic acid, and docosahexaenoic acid were significantly enriched in DA group, which was involved in the steroid biosynthesis and biosynthesis of unsaturated fatty acids. 3-phenylpropanoic acid, L-tyrosine, phedrine, rhizoctin B, and rhizoctin D were significantly enriched in TA group, which was involved in the phenylalanine metabolism and phosphonate and phosphinate metabolism. Conclusion We found that significant differences in gut microbiota composition and metabolite between Diannan small ear pigs and Diqing Tibetan pigs, which provide a theoretical basis for exploring the relationship between gut microbiota and pig breeds.


Introduction
Pigs are an important livestock species, and pork is essential in human daily meat protein intake (Xiao et al., 2018). Nowadays, pork represents more than 30% of the worldwide meat market, second only to chicken (Oh and See, 2012). The National Commission for Livestock and Poultry Genetic Resources of China has recognized 48 local pig breeds (Fang et al., 2005;Peng et al., 2022). Yunnan province is situated in China's southwestern region. Its unusual geographical position, complicated topography, unique climate, and various national farming traditions have produced abundant native pig resources, of which the Diannan small ear pigs and the Diqing Tibetan pigs are noteworthy examples (Shang et al., 2020;Wu et al., 2020).
Diannan small ear pigs live at the low altitude area in the south of Yunnan province, with an altitude of 800-1,300 m, belonging to the subtropical climate (Yonggang, 2010). Compared to European pig breeds such as Duroc, Landrace, and Yorkshire pigs, the Diannan small ear pigs has exceptional resistance to coarse food and excellent meat quality (Wang et al., 2016). Diqing Tibetan pigs live at an average elevation of approximately 3,500 m on the Yunnan Plateau (Li et al., 2013). Diqing Tibetan pigs have developed complex phenotypic and physiological adaptations to high-altitude hypoxia compared with lowland pigs Zhang B. et al., 2015). In summary, Diannan small ear pigs and Diqing Tibetan pigs have abundant phenotypic traits and superior commercial traits, adaptability to different environments (Figure 1), so it is regarded as valuable genetic resource that deserves to be efficiently utilized for scientific conservation and commercial exploitation.
A dense, dynamic, and highly complex microbial community inhabits the animal gastrointestinal tract (Yang et al., 2016). Intestinal microorganisms are essential and mobile animal metabolic organs that maintain body homeostasis (LeBlanc et al., 2013). Intestinal microorganisms may be influenced by the host heredity and environmental (Ley et al., 2006;Sutherland et al., 2006;Miao et al., 2009;Guo et al., 2011;Wu et al., 2013). Presently, comparative studies on the intestinal microorganisms of Diannan small ear pigs and Banna small pig inbred lines, Tibetan Pigs and Duroc pigs, Jinhua Pigs, and Landrace pigs have been conducted by the predecessors (Xiao et al., 2018;Wei et al., 2020;Yang et al., 2022). However, there are few reports on the comparative study of intestinal microorganisms between the Diannan small ear pigs and the Diqing Tibetan Pigs. Exploring the effect of host genetics on gut microbiota and metabolites can help reveal the metabolic interactions between microorganisms and hosts in vivo, which is important for understanding the gut health of hosts and improving the productivity of pig industry.
This study aims to analyze the composition characteristics of gut microbiota and metabolites in Diannan small ear pigs and Diqing Tibetan pigs and investigate the relationship between pig germplasm features and gut microbiota and metabolites. To provide a theoretical basis for exploring the relationship between gut microbiota and pig breeds.

Materials and methods
Animals and sample collection All animal procedures were performed according to the Guide for Animal Care and Use of Laboratory Animals in the Institutional Animal Care and Use Committee of Yunnan Agricultural University. The Department Animal Ethics Committee approved the experimental protocol of the Yunnan Agricultural University. Twenty Diannan small ear pigs (4 months of age, 39 ± 7.35 kg) and twenty Diqing Tibetan pigs (4 months of age, 38 ± 6.96 kg) raised by the pig-breeding institute of the Yunnan Academy of Animal Husbandry and Veterinary Sciences were selected. Each pig is housed in a pig house with the same environment and raised for 30 days using a standardized feeding scheme (Supplementary Table S1). Then, randomly select 6 Diannan small ear sows (DA group) and 6 Diqing Tibetan sows (TA group), collect fresh feces through rectal swabs, quickly freeze them in liquid nitrogen, and store them at −80°C.
Genomic DNA extraction and 16S rRNA gene sequencing The fecal DNA extraction kit was utilized to extract the samples' total DNA (Beijing tiangen biochemical technology co., LTD.). Using agarose gel electrophoresis, the purity and concentration of DNA were assessed, and an adequate amount of sample DNA was deposited in a centrifuge tube. The sample was diluted with sterile water to 1 ng/ μL. The diluted genomic DNA was utilized as the template for PCR expansion, with the V3-V4 region of the 16S rRNA gene being Geographical distribution of Diannan small ear pigs and Diqing Tibetan pigs.
Frontiers in Microbiology 03 frontiersin.org amplified. The extended primers were 338F (5′-ACTCCTACG GGAGGCAGCA-3′) and 806R (5′-GGACTACHVGGGTWTC TAAT-3′). The Barcode-specific primers were used. PCR products were identified using 2% agarose gel electrophoresis. According to the concentration of the PCR product, the same amount of samples was combined with the PCR product. After carefully combining the samples, the PCR product was purified by 2% 1 × TAE gel electrophoresis, and the target band was cut and recovered (Thermo Scientific company GeneJET adhesive recovery kit). Thermofisher's Ion Plus Fragment Library Kit 48 RXNS was used to construct the Library. After the library was validated by Qubit quantitative and library testing, Thermofisher's Ion S5™XL was implemented for computer sequencing.

Sequence processing and data analysis
Cutadapt (Version 1.9.1 1 ) was utilized to remove low-quality portions of reads, which were then matched to the species annotation database to identify chimeric sequences (Aßhauer et al., 2015). Chimera sequences were eliminated in order to acquire the final reliable data (Clean Reads) (Uparse version 7.0.1001 2 ) (Haas et al., 2011), and by default, the sequence was clustered as Operational Taxonomical Units (OTUs) with 97% identity. Mothur technique and SSUrRNA database of SILVA132 (threshold set to 0.8) were utilized to perform species annotation on OTUs sequences (Edgar, 2013). Rapid multiple sequence alignment utilizing MUSCLE software (Version 3.8.31) to determine system occurrence relationships for all OTUs sequences (Quast and Pruesse, 2013). QIIME software was used to examine the samples for Alpha diversity and Beta diversity (Version 1.9.1). By calculating the unifrac distance between OTUs in samples (Lozupone and Knight, 2005;Lozupone et al., 2011), unweighted group average clustering analysis (UPGMA) was created. The screening value of the linear discrimination criterion (LDA Score) was 2 using the LEfSe program for inter-group differential species analysis (Segata et al., 2011). Tax4Fun may be used to forecast the function of all samples (Qin et al., 2012).

LC-MS analysis of fecal metabolites
Collect the liquid sample for analysis 100 μL (0.1 mg fecal tissue liquid nitrogen grinding) put in centrifuge tube, add in 400 including 80% of methanol aqueous solution (4 times the volume of methanol), vortex oscillation, in −20°C stand for 60 min, 14,000 g, 4°C centrifugal 20 min, take a certain amount of supernatant put in 1.5 mL centrifuge tube, vacuum freeze-drying, the residue with 100 μL complex solvents to dissolve, vortex oscillation, 14,000 g, 4°C centrifugal 15 min, take that into the supernatant fluid sample LC-MS analysis.
Load the output data (raw) file into compound discoverer (CD) search software for straightforward filtering of retention duration, mass charge ratio, and additional characteristics. Several samples were then aligned according to a retention time deviation of 0.2 min and a mass deviation of 5 parts per million to improve the accuracy of the identification. The peak was recovered based on the specified mass deviation of 5 ppm, signal strength deviation of 30%, signal-to-noise ratio of 3, minimum signal strength of 100,000, and other data. The peak area was quantified, the target ions were integrated, the molecular formula was predicted and compared with the mzCloud database, the background ions were removed with blank samples, the quantitative results were normalized with quality control (QC) samples, and the identification and quantitative results of the data were obtained.

Correlation analysis
On the basis of the correlation coefficient (r) of the different microbiota and metabolites, SPSS 22.0 was used to conduct a Pearson correlation analysis, and cytoscape (version 3.5.1) was utilized to create a correlation heat map.

Gut microbiota: analysis and function prediction
Comparing the OTUs of the two groups determined that the number of identical OTUs in the DA and TA groups was 1,132, while the number of special OTUs in the DA and TA groups was 147 and 244, respectively ( Figure 2).
At the phylum level, the top 10 microorganisms in the relative abundance of the DA group and TA group are Firmicutes, Bacteroidetes, Actinobacteria, Spirochaetes, Proteobacteria, Tenericutes, Euryarchaeota, Kiritimatiellaeota, unidentified_Bacteria, Verrucomicrobia, among which Firmicutes and Bacteroidetes were the major phylum of bacteria in both the DA and TA groups, accounting for more than 90% of total microbiota in both groups ( Figure 3A). The relative abundance of Firmicutes and Actinobacteria in the DA group was higher than that in the TA Group. The relative abundance of Bacteroidetes, Spirochaetes, Proteobacteria, and Tenericutes in the TA Group was higher than that in the DA group.
At the genus level, the top ten microorganisms in the relative abundance of the DA and TA groups are Terrisporobacter, Bacteroides, Lactobacillus, unidentified_Clostridiales, Bifidobacterium, Faecalibacterium, unidentified_Ruminococcaceae, Streptococcus, unidentified_Prevotellaceae, Parabacteroides ( Figure 3B). The relative abundance of Terrisporobacter, Bacteroides, Bifidobacterium, Faecalibacterium, and Streptococcus in the DA group was higher than that in the TA group. The relative abundance of Lactobacillus, Frontiers in Microbiology 04 frontiersin.org unidentified_Clostridiales, unidentified_Ruminococcaceae, unidentified_Prevotellaceae, and Parabacteroides in the TA group was higher than that in the DA group.
After 97% identity clustering, the DA group obtained an average of 741 OTUs, while the TA group obtained an average of 847 OTUs. OTUs are then further annotated to calculate and compare the Alpha diversity between samples (Table 1). The Chao1 and ACE indices of the TA group were significantly higher than those of the DA group (p < 0.05), indicating that the TA group had a greater diversity of intestinal microbial than the DA group. UPGMA with weighted mean revealed that DA and TA samples are clustered in distinct branches, indicating that the microbiota structure of these two groups is distinct ( Figure 4).
LEfSe multistage discriminant analysis showed that there were 40 biomarkers (LDA score ≥ 2) at different classification levels in the intestinal microbial of pigs in the DA and TA groups ( Figure 5). At the phylum level, the relative abundance of Tenericutes, Fibrobacteres, Venn diagram of OTU. Atribacteria, Kiritimatiellaeota, Gemmatimonadetes, and Euryarchaeota in the TA group was significantly higher than that in the DA group (p < 0.05); At the class level, the abundance of Mollicutes, Methanobacteria, Deltaproteobacteria, Kiritimatiellae, unidentified_ Elusimicrobia, unidentified_Gemmatimonadetes in the TA group was significantly higher than that in DA group (p < 0.05). While the abundance of Ignavibacteria in the DA group was significantly higher than that in the TA group (p < 0.05); At the order level, the abundance of Methanobacteriales, Fibrobacterales, Elusimicrobiales, Methanomassiliicoccales, Bradymonadales, Gemmatimonadales in the TA group was significantly higher than that in DA group (p < 0.05); At the family level, the abundance of Prevotellaceae, Methanobacteriaceae, Fibrobacteraceae, Elusimicrobiaceae, Gemmatimonadaceae in the TA group was significantly higher than that in DA group. While the abundance of Lachnospiraceae in the DA group was significantly higher than that in the TA group (p < 0.05). At the genus level, the abundance of unidentified_Prevotellaceae, Candidatus_Soleaferrea, Methanobrevibacter, Rickettsiella, Fibrobacter, Elusimicrobium in the TA group was significantly higher than that in the DA group (p < 0.05). While the abundance of Rothia and Butyricicoccus in the DA group was significantly higher than that in the TA group (p < 0.05); At the species level, the abundance of Ruminococcus_sp_YE281, bacterium_enrichment_culture_clone_ BBMC_9, Methanobrevibacter_ruminantium in the TA group was significantly higher than that in DA group. In contrast, the abundance of Actinomyces_graevenitzii in the DA group was significantly higher than that in the TA group (p < 0.05).
Tax4Fun software was used for functional prediction analysis ( Figure 6). The result showed that at the KEGG level 3, DA and TA groups were mainly involved in carbon fixation pathways in prokaryotes, glycine serine and threonine metabolism, citrate cycle (TCA cycle), methane metabolism, one carbon pooly folate, sulfur metabolism, protein phosphatase and associated proteins, D-Glutamine and D-glutamate metabolism. Among them, carbon fixation pathways in prokaryotes, glycine serine and threonine metabolism, TCA cycle, methane metabolism, one carbon pool by folate, D-Glutamine and D-glutamate metabolism was significantly enriched in TA group, while sulfur metabolism, protein phosphatase and associated proteins was significantly enriched in DA Group.

Metabolome data analysis
In the positive and negative ion modes, the correlation of QC samples and R 2 values are close to 1 (Figure 7). The contributions of PC1 and PC2 were 16.59 and 23.86% in the positive ion mode and 22.00 and 16.28% in the negative ion mode, respectively. Both groups were effectively distinguished in the positive and negative ion modes ( Figure 8). 2,505 and 1,848 metabolites were detected in fecal samples of DA group and TA group, respectively. Among the top 10 metabolites with relative expression, oleamide is the metabolite with the highest relative expression in the positive and negative ion mode of the two groups (Supplementary Tables S2-S5).
The threshold was set as VIP > 1.0, Fold Change (FC) > 2.0 and FC < 0.5, and p value < 0.05, and the differentially expressed metabolites were screened out ( Figure 9). In the positive ion mode, there were 124 metabolites with significant differences between the In the same row, values with different letter superscripts mean significant difference (p < 0.05), while with no letter superscripts mean no significant difference (p > 0.05).

FIGURE 4
The beta diversity results of different grouped samples-UPGMA clustering tree.
KEGG analysis (Figure 11) showed that in the positive ion mode, the metabolic pathways with significant differences were mainly steroid biosynthesis involved by cholecalciferol, 5-dehydroepisterol and stigmasterol; phenylalanine metabolism involved by 3-phenylpropanoic acid, L-tyrosine, phedrine; methane metabolism involved by L-tyrosine, 7-oxoheptanoic acid; Biosynthesis of alkaloids derived from shikimate pathway involved by L-tyrosine, tetrahydroamine and ephdrine. In the negative ion mode, the metabolic pathways with significant differences were mainly phosphonate and phosphinate metabolism involved by rhizoctin B, rhizoctin D; Porphyrin and chlorophyll metabolism involved by stercobilin, pyrophoride a; Biosynthesis of unsaturated fatty acids involved by adrenic acid, docosahexaenoic acid.

Discussion
The composition of intestinal microorganisms in different breeds of pigs is different, affecting the body's physiological characteristics and digestive function . Our findings revealed that Firmicutes and Bacteroidetes are the most abundant phylum in the intestine of Diannan small ear pigs and Diqing Tibetan pigs, accounting for more than 90% of the total sequence, which is consistent with previous study results (Kim et al., 2011;Gao et al., 2019). The relative abundance of Firmicutes in Diannan small ear pigs is higher than that in Diqing Tibetan pigs. The relative abundance of Bacteroidetes is lower than that in Diqing Tibetan pigs. Firmicutes degrade fiber and cellulose (Armougom and Raoult, 2008), and Bacteroidetes help digest complex carbohydrates and ferment organic matter (Xu et al., 2019). Biddle et al. (2013) have shown that individuals with a high ratio of Firmicutes/Bacteroidetes have a higher efficiency of absorbing energy from food, thus effectively promoting the host's energy absorption and increasing the host's body weight. Therefore, the results of this study are consistent with the trait of more fat deposition in Diannan small ear pigs .
At the genus level, Bacteroides is related to the metabolism of proteins, various amino acids, and saturated fats (Zhang et al., 2014). Lactobacillus maintains intestinal and body health by inhibiting the invasion and colonization of the intestinal tract by pathogenic bacteria. At the same time, can enhance the digestive function of the  The volcano map of differential metabolites. intestinal tract (Valeriano et al., 2017). Bifidobacterium has critical physiological functions such as improving immunity, anti-aging, antitumor, and maintaining gastrointestinal health (Luo et al., 2018). Through studies in mice, Sokol et al. (2008) confirmed that Clostridium flexis has anti-inflammatory properties and plays a central role in maintaining the balance of intestinal flora. Therefore, we assume that the intestinal microbial composition of Diannan small ear pigs and Diqing Tibetan pigs may be related to their potential resistance to stress and disease. The diversity of gut microbes can reflect the health of the host. The diversity of gut microbes are usually positively correlated with the microbiota's stability, disease resistance and stress resistance. High diversity benefit overall health, and productivity . The environmental differences in different regions of Yunnan have formed the local pig breeding resources of Diannan small ear pigs and Diqing Tibetan pigs with excellent germplasm and unique traits Zhang B. et al., 2015). The results of this study show that the diversity of intestinal microorganisms in Diqing Tibetan pigs is higher than that in Diannan small ear pigs, which may be related to the fact that Diqing Tibetan pigs can adapt to the cold environment on the plateau. The interaction between genotype and environment resulted in the strong disease resistance and cold resistance of Diqing Tibetan pigs (Ge et al., 2018).
Host genetics have essential effects on specific strains of bacteria Zhang B. et al., 2015). This study showed significant differences in the intestinal microbial structure between Diannan small ear pig and Diqing Tibetan Pig. The relative abundance of Lachnospiraceae, Actinomyces, and Butyricicoccus in the feces of the Diannan small ear pigs was significantly higher than that of the Diqing Tibetan Pigs. Lachnospiraceae is a crucial beneficial bacterium which can ferment and degrade carbohydrates into SCFAs (Duncan et al., 2007;McKnite et al., 2012). Actinomyces are gram-positive bacteria with branching hyphae whose metabolites can block particular pathogens and modulate the body's immunity (Forster et al., 1997). Ding et al. (2016) found that Actinomyces metabolites isolated from wild animal feces can decompose cellulose, enzyme activity, antibacterial activity, and anti-tumor activity. Relevant studies isolated Butyricicoccus from the rumen of cattle, sheep, and deer. It was found that Butyricicoccus can also regulate the secretion of hemicellulose degrading enzymes, degrade polysaccharides (Louis and Flint, 2009), and produce butyrate to proliferate gastric epithelial cells (Dunne et al., 2012). These microorganisms may be associated with the Diannan small ear pigs' resilience to roughage and stress (Wang et al., 2016).
Genetic and environmental factors interact in the formation of intestinal microbial composition (Dowd et al., 2007;Kurilshikov et al., 2021). Diqing Tibetan pigs have lived in pastoral plateau areas for  (Hildebrand et al., 2013). Cellulose-related bacteria colonized in the intestine can catabolize dietary fiber in plant cell walls and produce a large amount of SCFAs for energy. Help pigs better adapt to fiber-containing diets and develop roughage tolerance. Prevotelaceae and Ruminocus are associated with cellulose decomposition . Our study on Diqing Tibetan pigs also detected that Prevotellaceae, Ruminococcus, were significantly higher than those of Diannan small ear pigs, which may be related to the herbivorous habits of the Tibetan pigs.
The metabolomics results of this study found that oleamide is the metabolite with the highest relative expression in the positive and negative ion mode of Diannan small ear pigs and Diqing Tibetan pigs. Oleamide is an endogenous fatty acid amide recognized as an endogenous sleep-inducing drug after its discovery in the cerebrospinal fluid of sleep-deprived cats (Boger et al., 1998). Oleamide promotes multiple biological functions, including sleep induction, anxiolytic effect, anti-inflammatory activity, and skeletal muscle atrophy inhibition (Cravatt et al., 1995;Ano et al., 2015;Moon et al., 2018;Kobayashi et al., 2021Kobayashi et al., , 2022. In addition, hydroalcoholic Frontiers in Microbiology 11 frontiersin.org lard extract containing oleamide exhibits anti-inflammatory effects on lipopolysaccharide-stimulated neutrophils and peripheral mononuclear cells, lowering the ability to create TNF-α and the activity of antioxidants and pro-inflammatory enzymes (Capó et al., 2021). As one of the primary agents for eliminating metabolites in animals, feces can potentially communicate a large quantity of information regarding the internal physiology of animals and serve as a source of chemical signals (Guan et al., 2022). Therefore, this study suggested that oleamide may be a potential marker of Diannan small ear pigs and Diqing Tibetan pigs' feces. However, the function and mechanism of oleamide as a metabolite in Diannan small ear pigs and Diqing Tibetan pigs still need to be further studied.
The metabolomics results of this study found that the steroid biosynthesis and biosynthesis of unsaturated fatty acids of Diannan small pigs were significantly higher than those of Diqing Tibetan pigs. Steroid hormones are a class of cyclopentane polyhydrophenanthrene derivatives, mainly divided into corticosteroids and sex hormones (Bovo et al., 2016). Corticosteroids can regulate physiological processes such as glucose metabolism, water, salt metabolism, and stress response (Kapugi and Cunningham, 2019). Sex hormones are mainly related to biological sex, the development of secondary sexual characteristics, and reproductive regulation (Wonyoung and Chungsoo, 2010). The regular operation of physiological functions requires hormone metabolism in the body to maintain a stable balance and be appropriately regulated. When the body is under stress, glucocorticoids in the blood increase to regulate the metabolism of sugar, water and salt, thereby enhancing immunity (Kuo et al., 2015). Glucocorticoids can induce adipocytes to differentiate and produce free or non-esterified fatty acids (Macfarlane et al., 2008); After persistent glucocorticoid release, animals develop metabolic problems such as hyperglycemia and obesity (Slavin et al., 1994). Glucocorticoids contribute to lipid metabolism and influence the expression of hormone-sensitive lipase and fat triglyceride lipase, according to the research of Vllena (Villena et al., 2004). Unsaturated fatty acids affect lipid metabolism in animal tissues. Bessa et al. discovered that during the hydrogenation of dietary polyunsaturated fatty acids by rumen microbes of lambs, some of the intermediate products and unsaturated fatty acids produced might be deposited directly into the muscle influencing the nutritional value of the muscle (Bessa et al., 2007). The pork quality of Diannan small ear pigs is outstanding, with an unsaturated fatty acid content of about 53.22%, an intramuscular fat content of about 8.5%, and a low lean meat rate of about 35%, which is consistent with the high expression of related metabolic pathways.
In addition, the phosphonate and phosphinate metabolism and phenylalanine metabolism pathways of Diqing Tibetan pigs were significantly higher than those of Diannan small ear pigs. Aminophosphonate compounds have good antibacterial activity, and studies have found that aminophosphonate compounds can inhibit  harmful Escherichia coli to a certain extent (Elsherbiny et al., 2020). Phenylalanine is an aromatic non-polar alpha amino acid, which can be directly converted into tyrosine to participate in protein synthesis, and it also plays a vital role in maintaining body health (MacDonald et al., 2020). These two metabolic pathways are consistent with the germplasm characteristics of Diqing Tibetan pigs with strong disease and stress resistance. In conclusion, according to the study results, it is speculated that the differential metabolic pathways enriched by differential metabolites are related to the specific germplasm characteristics of Diannan small ear pigs and Diqing Tibetan pigs. The link between intestinal microbiota and metabolites shows that intestinal microbiota is crucial for bodily health (Chow et al., 2010). Interestingly, our study found that Clostridium_disporicum was significantly positively correlated with hydrocortisone succinate.
Hydrocortisone is a steroid hormone that regulates multiple metabolic and physiological processes and is involved in response to many stresses (Pruett, 2003;Timmermans et al., 2019). Prior to excretion, the hydrocortisone metabolism is often complicated (Santamaria et al., 2021). During irreversible enzymatic processes, hydrocortisone is conjugated with sulfate or glucuronides to improve its solubility and permit its excretion in urine via the kidney or feces via the colon (Palme et al., 2005;Keegan, 2019). In addition, microorganisms in the gut can undertake side-chain cleavage of cortisol metabolites, which ultimately leads to the formation of fecal cortisol metabolites (Morris and Ridlon, 2017). Clostridium_disporicum is a saccharolytic species within Firmicutes and is associated with the decomposition of complex organic matter, but its actual actions in the intestine remain unknown (Li et al., 2021). Therefore, our study predicted that Clostridium_disporicum could degrade hydrocortisone to form hydrocortisone succinate in the intestine of Diqing Tibetan pigs.

Data availability statement
The data presented in the study are deposited in the GenBank repository, accession numbers OR162609-OR164377.

Ethics statement
The animal study was reviewed and approved by Institutional Animal Care and Use Committee of Yunnan Agricultural University. Written informed consent was obtained from the owners for the participation of their animals in this study.

Author contributions
SZ and XG were responsible for acquisition of the data and participated in drafting the manuscript. JZ, LY, HS, MY, YH, HP, HW, HZ, and YZ participated in conception and design of the experiment and analyses of the results, performed statistical analyses, and agreed to be accountable for accuracy and integrity of the data. All authors contributed to the article and approved the submitted version.

Conflict of interest
The authors declared that they have no conflicts of interest to this work.

Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.