Transmission of the gut microbiome in cohousing goats and pigs

Social interaction facilitates the horizontal transmission of the microbiota between different individuals. However, little is known about the level of microbiota transmission in different livestock animals and different digestive tracts. The Hainan black goat and Wuzhishan pig are typical tropical local breeds on Hainan Island in China. Thus, we sampled and analyzed the gut microbiome in Hainan black goats (cecum and rumen) and Wuzhishan pigs (cecum) to study horizontal transmission by rearing them in the same pen (six goats and six pigs) or separate pens (nine goats and nine pigs). De novo assembly and binning recovered 3,262 strain-level and 2,488 species-level metagenome-assembled genomes (MAGs) using ∼1.3 Tb sequencing data. Of these MAGs, 1,856 MAGs were identified as novel strain. Compared with goats living in separate pens, social interaction in the same pen promotes community homogeneity in the rumen microbiome (P < 0.05) and the cecum microbiome (P < 0.05), respectively. Notably, approximately 7.08% (231/3262) of the gut microbial population could transmit during cohousing, 12 strains only in inter-species transmission, versus 190 strains only in intra-species transmission, and 10 strains only in foregut and hindgut transmission. In addition, the social contact group has high transmitted strain abundance, which is correlated with community composition. This study provided a new insight into the influence of social interaction on the animal gut microbiota.

Social interaction facilitates the horizontal transmission of the microbiota between different individuals. However, little is known about the level of microbiota transmission in different livestock animals and different digestive tracts. The Hainan black goat and Wuzhishan pig are typical tropical local breeds on Hainan Island in China. Thus, we sampled and analyzed the gut microbiome in Hainan black goats (cecum and rumen) and Wuzhishan pigs (cecum) to study horizontal transmission by rearing them in the same pen (six goats and six pigs) or separate pens (nine goats and nine pigs). De novo assembly and binning recovered 3,262 strain-level and 2,488 species-level metagenome-assembled genomes (MAGs) using ∼1.3 Tb sequencing data. Of these MAGs, 1,856 MAGs were identified as novel strain. Compared with goats living in separate pens, social interaction in the same pen promotes community homogeneity in the rumen microbiome (P < 0.05) and the cecum microbiome (P < 0.05), respectively. Notably, approximately 7.08% (231/3262) of the gut microbial population could transmit during cohousing, 12 strains only in inter-species transmission, versus 190 strains only in intra-species transmission, and 10 strains only in foregut and hindgut transmission. In

Introduction
The digestive system harbors a diverse and complex microbiota which collectively modulates host health, development, and physiology (Wang et al., 2018;Guo et al., 2020). Colonization of the gut microbiota is a complex process that is affected by both host genetics (Zhou et al., 2018) and environmental factors (Eisenstein, 2020). In certain circumstances, the influence of environmental factors appears to outweigh host genetics in shaping microbiome composition (Henderson et al., 2015;Ji et al., 2018). Several recent studies in both humans and non-human primates provide strong evidence for the contribution of social interaction to microbiome assembly (Song et al., 2013;Tung et al., 2015;Li et al., 2016;Grieneisen et al., 2017;Sharma et al., 2019). Song et al. (2013) found that household members shared more microbiota than individuals from different households, thereby decreasing interindividual variation. Tung et al. (2015) showed that social interaction could predict microbiome structure in wild baboons. In addition, a study examined gut microbiome composition in nine sympatric wild non-human primate species and found that microbiota varied with host species but, importantly, also by social groups within species (Gogarten et al., 2018). Despite the social interaction playing an important role in gut microbiome composition, few studies focused on non-primate animals.
Use of the 16S rRNA gene analyses does not have sufficient taxonomic and functional resolution, while metagenomic shotgun sequencing has the potential to reveal microbial heterogeneity at the strain-level, which may be crucial to identifying and assessing the transmission of strains with distinct genetic repertoires (Nayfach et al., 2016;Yassour et al., 2018). Compared with bumblebees, honeybees often possess a higher level of microbial strain variation, which also exhibited more complex gene repertoires linked to polysaccharide metabolism (Su et al., 2021). In addition, the newly emerged worker bees share more strains with their queens, suggesting a vertical transmission of strains from queens to the newborn workers (Su et al., 2021). Yassour et al. (2018) found two patterns of mother-to-child microbial transmission at strainlevel. A Bacteroides uniformis strain lacking utilization gene cluster in the mother's strain was inherited to the child, implying that the monther's strain has a selective advantage in the infant gut (Yassour et al., 2018). Therefore, strain-level metagenomic profiling is needed to infer the transmission at higher taxonomic resolutions.
Hainan Island, in the South China Sea, has a tropical climate that is hot and humid, prompting the local animals to have differences from intensive commercial species in physical and life modes (Zhu, 2016). Hainan black goats and Wuzhishan pigs are native breeds that have long been bred in the unique natural ecological environment of Hainan Island (Jiang et al., 2020;Shao et al., 2021). In addition, the Wuzhishan pig is a famous miniature pig breed and is regarded as an ideal experimental animal model (Shao et al., 2021). Hainan black goats as ruminants have tender and delicious meat, making them a popular choice among consumers (do Rosario et al., 2016). Thus, to detect the microbiome transmission during cohousing, we chose two kinds of animals with different feeding habits and selected Hainan black goat to analyze the transmitted strain between the foregut and hindgut.
Using metagenome-assembled genomes (MAGs), we examined changes in the gut microbiome community and microbial transmission from cohousing goats and pigs. Our results provide insights into the gut microbiota during social interaction and expand our understanding of microbiome transmission. Among them, 1,856 new strainlevel MAGs broaden the useful microbial resources for the tropical climate system.

Materials and methods
Study design, sample collection, and sequencing In this study, 15 newborn Hainan black goats and 15 newborn Wuzhishan pigs were used, and they weaned at 3 and 1.5 months, respectively. After weaning, six goats and six pigs were cohoused in a new barn, while left goats and pigs still lived in the separate barns. Rumen and cecum samples were collected at weaning (0 day), 3 and 12 months after co-housing, respectively (Supplementary Figure 1). Three goats and pigs from each group were slaughtered 3 h later after the morning feeding at sample collecting day. All fresh samples were frozen in liquid nitrogen before being transferred to the laboratory in a dryice pack and promptly kept at −80 • C before total DNA extraction. DNA extractions followed the instructions using the DNeasy PowerLyzer PowerSoil Kit (Qiagen, United Kingdom). Shotgun sequencing was performed on an Illumina NovaSeq platform.

Metagenome assembly and annotation
Illumina sequencing data adaptors were trimmed using Trimmomatic (Bolger et al., 2014) and the subsequent trimmed reads were aligned to the reference genome to remove host-derived DNA contamination. The clean reads per sample were performed using MegaHit , which is included in MetaWRAP (Uritskiy et al., 2018). The contigs were retained for the subsequent analysis with the length ≥ 1 kb. Metagenomic binning was applied to recover individual genomes using MateBAT2, Maxbin2 and CONCOCT in the binning module. Three sets of draft bins were further analyzed using the bin_refinement module of the MetaWRAP pipeline with options -c 50 -x 10. The bins ≥ 50% complete were re-assembled using the SPAdes. To determine whether the same bins have been reconstructed via different samples, all bins were dereplicated using dRep at a threshold of 99% average nucleotide identity (ANI). The CheckM estimate results of MAGs meeting the following thresholds were retained: Genome Quality ≥ 50; contamination ≤ 10%; and completeness ≥ 50% [genome Quality = Completeness−5 × Contamination + log(N50)].
ORFs were predicted by Prodigal (v2.6) based on the default parameters and the function assignment was performed using the Prokka software (Seemann, 2014). The 16S RNA genes in MAGs sequences were predicted using Barrnap (v 0.9). Then, the MAGs-derived proteins were searched against the carbohydrate-active (CAZy 1 ) database using dbCAN2  and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (Kanehisa et al., 2017) using KofamKOALA (Aramaki et al., 2020). Reads from each sample were aligned against the contigs of MAGs by the Salmon software (Patro et al., 2017). The abundance of each contig was calculated in contigs per million (CPM), and then the median CPM of the contig was used as the abundance of MAGs.
Taxonomic classification, phylogenetic trees, and species-level clustering of metagenome-assembled genomes The Genome Taxonomy Database Toolkit (GTDB-Tk 2 ) pipeline was used for the taxonomic classification of the 3,262 MAGs (Parks et al., 2018). In short, the GTDB-Tk classified the bins based on ANI to the reference genome, placement in the bacterial or archaeal reference genome tree, and relative evolutionary distance. MAGs were clustered into species-level genome bins at the threshold of 95% ANI. Phylogenetic relationships among the 2,488 species-level genome bins were computed by the total branch length based on proteins. The phylogenetic tree was produced from concatenated protein sequences using PhyloPhlAn and re-rooted manually using iTOL at the branch between archaea and bacteria. The tree was used to calculate the phylogenetic gain at different taxonomic levels using the pd_clade routine in GenomeTreeTk (v0.1.6). 3 We downloaded 48,622 microbial genomes from GTDB (31,910 genomes) (Parks et al., 2018), ruminant gastrointestinal tract (10,737 genomes) (Xie et al., 2021), and pig gut (6,339 genomes) . All these genomes were used as reference genomic datasets for the identification of novel microbial genomes from 3,262 MAGs using FastANI (Jain et al., 2018). Genomes were determined as novel strains based on < 99% ANI and defined as novel species based on < 95% ANI. All the sequenced reads from our study were mapped to two datasets using Kraken2 to estimate the read classification. Two combined databases: a common database consisting of the bacterial, archaeal, fungal, and protozoan genomes in GenBank, and the GenBank database plus the MAGs.

Strains filtering, identification of microbial single nucleotide variants and detection of microbial transmission
To determine the presence or absence of strains in samples, the strains were filtered to satisfy the following criteria: horizontal coverage (breadth) of ≥ 10%, average vertical coverage (depth) ≥ 0.15 x (Schmidt et al., 2019). Microbial Single Nucleotide Variants (SNVs) were called using metaSNV software (Costea et al., 2017). In short, the transmission score (S T ) quantifies how much the similarity between two sample SNV profiles from social contact group, and between foregut and hindgut SNV profiles within an individual. Each potential SNV requires support by at least two reads at a base call quality of Phred ≥ 15. To identify the microbial transmission event, we applied a text mining approach (Schmidt et al., 2019) to calculate the S T . Next, we calculated the incidence of each allele across screening samples. For these, we calculated the probabilities in paired samples observations for each of the four possible cases: any given allele i could either be present in both samples (p 1 , 1 ), absent in both samples (p 0 , 0 ), or present in one but absent in the other sample (p 1 , 0 and p 0 , 1 ): For each pair of aligned samples, we calculated the raw and logarithmic probability of SNV profile overlap (L obs ) observed, which quantifies how likely the observed average allele profile agreement between two samples.
Similarly, we computed the log-likelihood of the least likely agreement case (L min ) per allele: Next, a raw probability score (P raw ) for the observed allele agreement between a given paired samples is calculated: P raw = L obs / L min . We defined the S T (t, s) for starin t in subject s as a standard Z score of the P raw : S T = (P raw (s) − µ raw ) / σ raw .
Microbial transmission events that were identified in pair of samples from social contact group and the foregut-hindgut pair of samples from the same individual have been divided into three types: intra-species transmission events, in which goats or pigs spread microbiota within their own species; inter-species transmission events, in which microbiota transmit between goats and pigs; and foregut-hindgut transmission events, in which microbiota spread between the foregut and hindgut within an individual.

Statistical analysis
Permutational multivariate analysis of variance (PERMANOVA) quantifies the variation among the samples that partitions social contact group and control group distances to permit the assessment of the effect of an exposure or intervention (grouping factor) on the sampled microbiome (Kelly et al., 2015). PERMANOVA was performed using the adonis2 function from the Vegan package with default settings with 999 unrestricted permutations and the Monte Carlo P value was calculated. P < 0.05 was considered statistically significant. Spearman correlation analysis was used to assess relationships between stain richness and community composition and transmission strain abundance in the control group or social contact group (goat rumen, goat cecum, and pig cecum).
Using species-level thresholds (≥ 95% ANI and ≥ 65% alignment), 1,856 of the MAGs did not match any available databases (Supplementary Table 1). Of these, 663 were highquality genomes, 835 were medium-quality genomes, and 358 were low-quality genomes. The 1,856 MAGs were assigned to 27 phyla, 70 orders, 131 families, and 336 genera, with 17.6% of these MAGs unable to be assigned to a recognized genus, implying that a significant fraction of the MAGs are likely novel genera. Prevalent among these MAGs classified at the order level were Oscillospirales (28.0%), Bacteroidales (27.9%), and Lachnospirales (10.5%), while the top genera all belonged to the order Bacteroidales, including Prevotella (5.6%) and RC9 (4.8%). In addition, to understand the phylogenetic position of the uncultured strains, we placed the 3,262 MAGs in a maximum-likelihood tree and found that the unknown MAGs improved phylogenetic diversity by an average of 49.8 and 11.3% for bacterial and archaeal lineages, respectively ( Figure 1B).

Functional characteristics of metagenome-assembled genomes in different gastrointestinal tracts
We explored the proteomic contents of our datasets and their putative functions. A total of proteins that could be annotated according to one or both approaches were 2,753,800, of which 2,589,938 and 295,216 proteins were annotated by using the KEGG database and CAZy, respectively (Supplementary Tables 2, 3).
To investigate the ability of the identified microbial to produce short-chain fatty acids (SCFAs), we explored the genes that are required to produce some components such as acetate, propionate, and butyrate (Supplementary Table 4). As a result, genes required for the complete acetate-producing pathway were found in 1,248 MAGs, which were generally Taxonomic annotation and phylogenetic tree of 3,262 metagenome-assembled genomes (MAGs). (A) The maximum-likelihood tree of the 3,262 MAGs identified in this study was produced from concatenated protein sequences using PhyloPhlAn. (B) The phylogenetic gain contributed by the microbial tree of the gut provided by the unknown species-level genome bins is shown as proportional increases in branch length per phylum (left) and absolute branch lengths (right). (C) Comparison of the read classification rates of all the digesta samples using the following datasets: a common database consisting of all complete microbial genomes in RefSeq and the GenBank database plus MAGs. The Wilcoxon rank-sum test was used to assess the differences. assigned in the Firmicutes_A and Bacteroidota phyla. Complete enzymes for the propionate production (via propionate pathway by the conversion of L-Lactate, succinate, and acetyl-CoA to propanoyl-CoA) were found in 239 MAGs of the Bacteroidota and Firmicutes_A. Enzymes encoded by genes for butyryl-CoA dehydrogenase, methylmalonyl-CoA mutase, methylmalonyl-CoA epimerase, enoyl-CoA hydratase, formate C-acetyltransferase, pyruvate ferredoxin oxidoreductase alpha subunit, acetate kinase, and phosphate acetyltransferase were common in these genomes. Sixty MAGs were identified to form butyrate (from acetyl-CoA through the typical pathway via butyrate kinase or acetate CoA-transferase).
Carbohydrate active enzymes (CAZymes) are enzymes involved in the metabolism and binding of carbohydrates. The majority of classes belonged to glycoside hydrolases (GHs), with 141,404 genes identified. The next glycosyltransferases (GTs; 74,345) and carbohydrate esterases (CEs; 50,099) were the most annotated groups (Supplementary Table 5). And then, the remaining three categories including carbohydrate-binding modules (CBMs; 47,009), polysaccharide lyases (PLs; 15,498), and auxiliary activities (AAs; 13,240), contained fewer proteins. Bacteroidota, Firmicutes_A, and Verrucomicrobiota showed high proportions of GHs, while Bacteroidota and Firmicutes_A showed high proportions of GTs. Notably, the most commonly identified GH was the GH109 family, which combines many enzymes responsible for the degradation of glycoproteins. Furthermore, is followed by GH2 and GH3 families, which are mostly related to cellulose and hemicellulose degradation.
Using permutational multivariate analysis of variance, we investigated the effects of host and environmental variables on the microbiota. The environmental factor exerted the main effects in both the single-and multi-factor analyses, followed by the developmental stage ( Figure 2B). Since the sampling strategy applied in this study was optimized to compare microbial compositions across the control and social contact groups, we focused on the comparative analysis between the control group and the social contact group and left the comparisons of other factors for future studies. Moreover, we analyzed the relative abundances of the MAGs in the control group and social contact group, and observed no significant differences in alpha diversities between the two groups in goat rumen, goat cecum, and pig cecum, as evaluated by the Shannon and richness indexes (Supplementary Figure 4). The social contact microbial communities exhibit significantly similar features in goat rumen and goat cecum (P = 0.014, P = 5.4e-05, Figure 2C). However, the social contact microbial communities of the pig cecum were similar but not significant ( Figure 2C). Overall, these results indicate that microbiota may be transmitted during the microbiota in co-housed animals.

Transmission of microbes during the social interaction
We further investigated the transmission of strain-level MAGs following cohousing by comparing different species of animals that cohabited. In addition, we also determined the transmission strains in the foregut and hindgut in the same sample. The 7.08% of the microbiota (231 out of a total of 3,262) were detected as transmission microbiota (Figure 3A). The majority of these strains were only assigned to withinspecies microbial transmission (190,82.25% out of 231), after which the between-species group (12, 5.19%) was the most assigned group, and the foregut-hindgut transmission events (10, 4.33%) contained fewer strains ( Figure 3B). The transmitted strains were assigned to 16 phyla, and the dominant phyla were Bacteroidota (49.78%) and Firmicutes_A (22.51%) ( Figure 3C). Transmissible strains accounted for a high fraction of classifiable microbial abundance in the social contact group (Figure 3D). The abundance of transmission microbiome in the social contact group was found to negatively correlate with community composition in goat rumen (ρ Spearman = −0.65). Goat cecum transmitted strains abundance correlated with strain richness (social contact group: ρ Spearman = −0.77; control group: ρ Spearman = 0.77), and with community composition (social contact group: ρ Spearman = 0.54; control group: ρ Spearman = 0.65). Pig cecum transmitted strains abundance (control group) correlated with strain richness (ρ Spearman = 0.82) and community composition (ρ Spearman = −0.71) (Figure 3E).

Discussion
Social interaction, in an enriched social environment, plays an important role in an individual's health (Spor et al., 2011;Iida et al., 2013). One mechanism for this effect is the interaction of animal-associated microbiota (Rothschild et al., 2018). This study on the effects of cohabitation with different species on the gut microbiome in non-primates. From our dataset, we constructed 3,262 MAGs and 1,856 novel MAGs, considerably extending prior MAGs (Parks et al., 2018;Chen et al., 2021;Xie et al., 2021). The dataset increased the reads classification rate by over 30% and provides the possibility for analysis of microorganism transmission. The gut microbiome similarity within the social contact group increased during cohabitation. Moreover, we found that within-species transmission is the main type in the gut microbiome. Change in the gut microbiota following environmental conversion. (A) Relative abundance of microbial phyla for the goat rumen, goat cecum, and pig cecum samples. (B) A bar plot summarizing the results of single and multiple factors permutational multivariate analysis of variance (PERMANOVA) results in the goat rumen, goat cecum, and pig cecum samples. The variations were derived from the between-sample Bray-Curtis distance. (C) Stability of strain-level composition profiles between control and social contact group as measured by Bray-Curtis dissimilarity. The Wilcoxon rank-sum test was used to assess the differences.
Our results highlight the importance of cohousing mediated transmission in shaping gut microbiomes. Previous studies have reported that social interaction can lead to a similarity in the structures of the gut microbiome (Tung et al., 2015;Zhu et al., 2020Zhu et al., , 2021. Although these studies have great ecological validity, it is difficult to understand the effects of social interactions on microbiome similarity from the combined effects of shared dietary intake or exposure to local environmental microbes. Animals in our study were fed identical diets in each species. Our findings indicated that cohabitation with diverse species had an effect on the gut microbiome community by increasing gut microbiota similarity. Despite the small number of transmitted trains in the gut microbiome (goat rumen, goat cecum, and pig cecum), we discovered that transmitted strains were frequent throughout microbial phyla. Previous research reported that humans have a greater gut microbiome among households and between couples than non-related people, and there was no consistent indication of transmission across any specific phyla (Brito et al., 2019). All these results showed that transmission may be largely driven by chance events and indirect transfer. Although we could not determine whether or not direct transmission occurred and the direction of transmission, our results showed that the number of intra-species transmission strains was higher than that of inter-species transmission strains, which may be correlated with diet (Zhu et al., 2021). Based on the abundance analysis of the putatively transmitted gut microbiome, we observed that the social contact group has a high proportion of transmitted strains, and correlates with the microbial community. One possible reason for this pattern is that when microbiota colonizes a new host, they become more abundant to occupy the available gut microbial niches (Tung et al., 2015).
Using culture-independent approaches to recover the genomes, which have become commonly used for species discovery and characterization (Almeida et al., 2019). In this study, we constructed 1,856 novel strains from our dataset. Among all the MAGs, the novel bacterial taxa derived from the Clostridia class were more abundant, which were mainly distributed in the ten families such as Lachnospiraceae, Acutalibacteraceae, Oscillospiraceae, CAG-272, Ruminococcaceae, Borkfalkiaceae, CAG-74, Saccharofermentanaceae, and CAG-138. Many commensal Clostridia species belong to a large group of obligate anaerobic and highly diverse bacteria, which are thought to be responsible for the maintenance of gut homeostasis. The high diversity of unclassified genomes affiliated with Clostridia has also been reported in the cecum microbiome (Glendinning et al., 2020) and the bovine rumen microbiome (Stewart et al., 2019). Phylogenetic analysis revealed that unknown genomes expand the diversity of bacterial and archaeal lineages. In addition, we The proportion of transmission microbiota in total microbiome abundance in the control and social contact groups from goat rumen, goat cecum, and pig cecum. The Wilcoxon rank-sum test was used to assess the differences. (E) Tests for the relationship between the abundance of transmission strains in the control or social contact groups (goat rumen, goat cecum, and pig cecum) and stain richness (ρ Spearman ) and community composition (ρ Spearman ).
provided a dataset of gut microbial genomes, which not only increases the reads classification but also enables us to analyze the transmission of microorganisms.
The Short-chain fatty acids (SCFAs), mainly acetate, propionate, and butyrate, are the principal products of the microbial fermentative activity in the gastrointestinal tract and have significant effects on intestinal health and immunity (Koh et al., 2016;Baaske et al., 2020). Propionate that is mainly transported to the liver is an important energy source for the host (Koh et al., 2016). The majority of bacteria belong to the Lachnospiraceae family, carrying complete genes for propionate production. Moreover, butyrate is an important energy source for the ruminal epithelium enterocytes and is involved in the maintenance of gut homeostasis (Koh et al., 2016;McNabney and Henagan, 2017;Baaske et al., 2020). In this study, the majority of taxa encoding enzymes required for butyrate production were Bacteroidota and Firmicutes_A. Notably, the Bacteroidaceae and Lachnospiraceae families carrying genes from the butyrate pathway have the metabolic capability to degrade and utilize plant-derived fibers as nutrients (Vital et al., 2014;Onrust et al., 2015;Segura-Wang et al., 2021). Furthermore, all genes required for reductive acetate in the Lachnospiraceae family, have also been reported in other studies (Bhattacharya et al., 2015). Therefore, the identification of strains carrying the enzymes necessary for acetate, propionate, and butyrate production, is an important step toward future modulation of the gut microbiota to improve the growth performance of animals.
The gastrointestinal tract microbiome carries a large repertoire of genes encoding enzymes that can breakdown polysaccharides (Lillington et al., 2020;Wardman et al., 2022). In the current study, GHs were discovered in the gut microbiota of a variety of species that can catalyze the hydrolysis of glycosidic bonds found in complex carbohydrates (Bohra et al., 2019). The high abundance observed for the GH109 family underlines its importance, which could cleave the terminal N-acetylgalactosamine from the A antigen (Wardman et al., 2022). Furthermore, the GH2 and GH3 families contain widely distributed enzymes with a range of degradative activity on plant cell wall polysaccharides (Dai et al., 2015;Chang et al., 2018). This provides a basis for more effort to explore the gastrointestinal tract microbiome, and we will be able to identify a large number of enzymes of industrial value.

Conclusion
This report offers novel insight into the relevant microbial transmission patterns within cohousing goats and pigs. It suggests that within-species microbial transmission is the main transmission pattern. Transmission of potential microbiome communities is relevant to the overall microbiome composition of individuals, and thus social interaction-derived changes in microbiomes may potentially impact the health status of individuals. In addition, a total of 3,262 strain-level MAGs, 2,488 species-level MAGs, and 1,856 unknown strain-level MAGs from our dataset were reconstructed and detected from the gut microbiome.

Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI SRA BioProject, accession no. PRJNA792486.

Ethics statement
The animal study was reviewed and approved by Northwest A&F University. Written informed consent was obtained from the owners for the participation of their animals in this study.

Author contributions
TZ, ZL, XW, and YJ designed the research project. TZ, TS, and YY analyzed the data. TZ and ML collected animal samples. TZ wrote the manuscript with significant input from all authors. TZ, ZL, ZN, and YJ edited the language of the manuscript. All authors contributed to the article and approved the submitted version.

Funding
This study was supported by the National Natural Science Foundation of China (31902126 and 31822052), the Key Research and Development Program of Hainan Province (ZDYF2022XDNY153), the Central Public-interest Scientific Institution Basal Research Fund for Chinese Academy of Tropical Agricultural Sciences (No.1630032022006), and the Project of Agricultural Germplasm Resources Protection in Hainan Province. We thank the High-Performance Computing Platform of Northwest A&F University for providing the computing resources. SUPPLEMENTARY FIGURE 4 (A-C): Boxplots are used to summarize the abundance distributions of the abundance of taxa of interest in the gut sections (goat rumen, goat cecum, and pig cecum). The y-axis in these panels (except C) indicates the relative abundances, while that in C indicates the abundance ratios between Firmicutes_all and Bacteroidota. Firmicutes_all consists of the following phyla: Firmicutes, Firmicutes_A, Firmicutes_B, and Firmicutes_C. (D,E) Microbiota richness and diversity (Shannon index) analyses for each host organ.