Next Article in Journal
Traditional Meat Products—A Mycotoxicological Review
Previous Article in Journal
Conservation Priorities and Demographic History of Saussurea involucrata in the Tianshan Mountains and Altai Mountains
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Different Feeds on the Fungi Microbiome of Suffolk Crossed with Tibetan Sheep

1
Institute of Livestock Research, Tibet Academy of Agricultural and Animal Husbandry Sciences, Lhasa 850000, China
2
Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lhasa 850000, China
3
Institute of Traditional Chinese Veterinary Medicine, College of Veterinary Medicine, Nanjing Agricultural University, Nanjing 210095, China
4
MOE Joint International Research Laboratory of Animal Health and Food Safety, College of Veterinary Medicine, Nanjing Agricultural University, Nanjing 210095, China
5
Faculty of Veterinary and Animal Sciences, The Islamia University of Bahawalpur, Bahawalpur 6300, Pakistan
6
Department of Biochemistry, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
*
Author to whom correspondence should be addressed.
Life 2023, 13(11), 2210; https://doi.org/10.3390/life13112210
Submission received: 12 October 2023 / Revised: 7 November 2023 / Accepted: 9 November 2023 / Published: 14 November 2023
(This article belongs to the Section Microbiology)

Abstract

:
The gut microbiome plays an important role in the metabolism, nutrient absorption and immunocompetency of animals. The dynamics of the microbiota can be influenced by modulatory factors that involve nutrition, environment, health, diseases, etc. Few reports have been documented regarding the effects of different feeds on the fungi microbiome of Suffolk crossed with Tibetan sheep. A total of 30 Suffolk crossed with Tibetan sheep (ST sheep) were selected for the study and randomly divided into five equal groups (n = 6): AZ, BZ, CZ, DZ and EZ. Group AZ was fed with alfalfa and oat grass, whereas group BZ was fed with mixture of concentrated feed, alfalfa and oat grass. Groups CZ, DZ and EZ were fed with concentrated feed #1, #2 and #3, respectively. All experimental animals were fed twice a day for four months, and rectum samples were collected for microbiota analysis. Results revealed that 2,781,461 raw reads and 2,333,239 clean reads were achieved in the ST sheep. When compared with the sheep of groups AZ and BZ (164), the shared amplicon sequence variants (ASVs) between AZ and CZ (109), AZ (113) and DZ (118) as well as AZ along with EZ were fewer. Conspicuous different phyla (8) and genera (56) were examined and compared with free-range sheep in AZ. Genera including Xeromyces, Kazachstania, Cordyceps, Rhodotorula, Pichia, Spor, etc. were found higher in animals in the CZ, DZ and EZ groups. The results of this study provide new insights regarding the effects of different feeds on the fungi microbiome of sheep farmed on the plateau. We concluded that the differences in feed in Suffolk crossed with Tibetan sheep altered their gut microbiota.

1. Introduction

Sheep are an important ruminant that provide mutton, milk and cashmere to human beings [1,2]. With the improvement in living standards, there is a growing demand for mutton and its products. The Qinghai–Tibetan Plateau in China is home to many wild and domestic animals [1], of which sheep and yaks are considered important economically, as they provide indispensable elementary resources to local pastoralists [2].
Evolutionary changes in Tibetan sheep have made them adaptable to a harsh environment, as the plateau is has limited forage [3]. Thes plateau sheep not only produce economic and agricultural products but also have significant cultural importance [4,5]. The Suffolk sheep is a worldwide distributed mutton breed, and Suffolk crossed with Tibetan sheep (ST sheep) are popular animals on the plateau because of their fast weight gain and adaptability to the plateau environment.
It has been documented that multiple species of fungi, protozoans and bacteria make up the microbiome [6], which contributes greatly to a host’s metabolism, nutrient absorption and immune status [7]. The abundance of fungi in the gut biome is lower than that of bacteria but still is of significant importance to the health of animal [8]. Previous studies have revealed that colitis and inflammatory bowel and liver diseases are highly co-related to the dysbiosis of fungal microbiota [9,10,11,12]. Microbiota can be influenced by modulatory factors that include nutrition, environment, health, diseases, medications, etc. [13]. Diet plays a crucial role in altering the microbiota by changing the abundance of bacterial species [14]. It has been reported that animals fed with different forages had different gut microbiota [15,16,17]. However, few reports have been documented regarding the effects of different feeds on the fungal microbiome of Suffolk crossed with Tibetan sheep on the plateau. Therefore, the current study aimed to compare the fungal microbiome of Suffolk crossed with Tibetan sheep fed with different feeds.

2. Materials and Methods

2.1. Experiment Design

A total of thirty (30) three-month-old male Suffolk crossed with Tibetan sheep (ST sheep) with weights of 17.21 ± 1.34 Kg were chosen and divided into five groups, namely AZ, BZ, CZ, DZ and EZ. All the selected sheep were healthy, non-castrated meat animals, and the average weights among all sheep groups were almost same. The ST sheep in group AZ were fed with alfalfa and oat grass; those in group BZ were fed with a mixture of concentrated feed, alfalfa and oat grass; those in groups CZ, DZ and EZ were fed with concentrated feed #1, #2 and #3 (Table S1), respectively. The mixture of concentrated feed in group B was mixed with concentrated feed 1#, 2# and 3# in an equal amount. All the ST sheep were offered feed two times a day for four months, while animals had ad libitum access to clean, fresh drinking water. The animals remained healthy during the experimental period, and no specific treatment strategy was adopted. The type of breeding system of the animals was semi-wild.
The body weight of the ST sheep was measured, and fecal samples from the rectum of all the sheep were collected with cotton swabs and stored in sterile tubes till examination. All samples were transferred to the Institute of Livestock Research of the Tibet Academy of Agricultural and Animal Husbandry Sciences on Drikold and kept at −80 °C in a freezer until further analysis.

2.2. Suffolk Cross with Tibetan Sheep Microbiome Analysis

The genomic DNA extraction from each ST sheep group (n = 6) were performed using a DNA Stool Mini Kit (QIAamp, Hildan, Germany) according to the guidelines. The quality inspection and quantity detection of the DNA products were carried out through 1.2% agarose gel electrophoresis and a NanoDrop one (Thermo Scientific, Waltham, Massachusetts, USA). The ITS1 region was amplified by using specific primers pairs (ITS1 F:5′-CTTGGTCATTTAGAGGAAGTAA3′; ITS1 R:5′-GCTGCGTTCTTCATCGATGC3′) as per a previous study [18] (Table A2), and the PCR products’ quality and quantity were also examined. Then, library construction was performed by piloting a TIANSeq DirectFast Library Kit (Tiangen Biotech Co., Ltd., Beijing, China). After that, the present library was sent for sequencing via the Illlumina MiSeq platform (Bioyi Biotechnology Co., Ltd., Beijing, China).
The generated raw data from the Illlumina platform were initially passed through quality control DADA2 to achieve the amplicon sequence variants [19]. The ASVs were aligned to the UNITE database (ITS) to produce the taxonomy table [20]. Alpha diversity indexes observed operational taxonomic units (OTUs), and the Chao1, Shannon and Faith’s phylogenetic were calculated using QIIME2. Moreover, the beta diversity indexes, which included principal coordinate analysis, non-metric multidimensional scaling and the unweighted pair group method were measured via arithmetic means. Different fungi among the ST sheep groups were examined through analysis of variance (ANOVA), LEfSe and DEseq2 [21,22]. Function predication of the KEGG Ortholog of the microbiome of the ST sheep was performed through PICRUSt [23].

2.3. Statistical Analysis

All the data from the ST sheep were evaluated via ANOVA and Dunn’s test by utilizing IBM SPSS (27.0). Data are presented as the means ± SD and were considered statistically significant at p < 0.05.

3. Results

3.1. Generated Sequencing Data in the ST Sheep in Different Feeding Groups

At the end of this experiment, the body weights of the ST animals in the CZ (p < 0.05), DZ (p < 0.0001) and EZ (p < 0.05) groups were significantly higher than those of the ST sheep in the AZ group (Figure S1). In the present study, a total of 2,781,461 raw data and 2,333,239 filtered data were generated (Table S3). These data produced 3, 577 ASVs, with 44 shared ASVs among the current sheep groups. Group AZ shared 164, 109, 113 and 118 ASVs, respectively, compared with groups BZ, CZ, DZ and EZ (Figure 1a). Then, the ASVs were aligned to different taxonomies, and 15, 15, 16, 14 and 14 phyla were found in the AZ, BZ, CZ, DZ and EZ groups, respectively. A total of 138, 184, 181, 201 and 221 genera were found in the five ST sheep groups, respectively (Figure 1).

3.2. Comparing Analyses of the Fungal Microbiota in ST Sheep in Different Feeding Groups

Rarefaction curves of faith_pd and features of shannon_entropy were observed in all of the five ST sheep groups with saturated curves, which demonstrated that all the sequencing data were adequate and that the richness and evenness of the microbiota in the ST sheep were sufficient (Figure 2). The Alpha diversity index analysis found that there was no significant difference in the chao1, faith_pd, observed_features, shannon_entropy and Simpson, respectively, among all groups (Table S4, Figure S2). The beta diversity index analysis found that there were no significant differences among the groups through the analyses of PCA and PCoA, whereas significant differences were found via NMDS and Qiime2 (Figure 3). The pairwise ANOSIM analysis found that the distances between the AZ, BZ (p < 0.05), CZ (p < 0.05), DZ (p < 0.01) and EZ (p < 0.05) groups were statistically significant. Furthermore, significant differences were found between the BZ, CZ (p < 0.05), DZ (p < 0.01) and EZ (p < 0.01) groups (Figure S3). The primary species in different ST sheep groups for different taxa are shown in Table 1 and Figure 4.

3.3. Revealing Marker Fungi Species in ST Sheep in Different Feeding Groups

A grouping and clustering heat map analysis found a higher abundance of phyla of Neocallimastigomycota and Chytridiomycota in AZ; Neocallimastigomycota and Bryophyta in BZ; Chlorophyta and Bryophyta in CZ; Ascomycota, Rozellomycota, Anthophyta, Cnidaria and Olpidiomycota in DZ; and Ascomycota, Cercozoa and Blastocladiomycota in EZ (Figure 5a). At the class level, the abundance of Neocallimastigomycetes was higher in AZ and BZ. The abundances of Sordariomycetes, Dothideomycetes, Tremellomycetes, Cystobasidiomycetes and Rhizophydiomycetes were higher in CZ; Saccharomycetes, Sordariomycetes, Microbotryomycetes, Ustilaginomycetes and Rhizophydiomycetes were higher in DZ; and Eurotiomycetes, Sordariomycetes, Wallemiomycetes, Leotiomycetes and Ustilaginomycetes were higher in EZ (Figure 5b). At the order level, Neocallimastigales was higher in AZ and BZ. Higher abundances of Pleosporales, Hypocreales, Filobasidiales, Sordariales and Microascales in CZ; Saccharomycetales, Sordariales, Capnodiales and Cantharellales in DZ; and Eurotiales, Wallemiales, Capnodiales and Helotiales in EZ were observed (Figure 5c). At the family level, Neocallimastigaceae was higher in AZ and BZ. Higher abundances of Filobasidiaceae, Nectriaceae, Trichocomaceae, Didymellaceae and Microascaceae in CZ; Debaryomycetaceae, Trichocomaceae, Cladosporiaceae, Microascaceae and Hypocreales_fam_Incertae_sedis in DZ; and Aspergillaceae, Wallemiaceae and Cladosporiaceae in EZ were observed (Figure 5d). At the genus level, higher abundances of Orpinomyces and Pecoramyces in AZ; Piromyces in BZ; Aspergillus, Naganishia, Talaromyces and Fusarium in CZ; Scheffersomyces, Cladosporium, Xeromyces and Talaromyces in DZ; and Aspergillus, Wallemia and Cladosporium in EZ were detected (Figure 5e). An evolutionary tree of species with a heat map revealed higher abundances of the classes of Agaricostilbomycetes in AZ; Neocallimastigomycetes in BZ; Bryopsida, Tritirachiomycetes and Trebouxiophyceae in CZ; Anthozoa, Archaeosporomycetes, Exobasidiomycetes, Olpidiomycetes, GS13 and Rhizophydiomycetes in DZ; and Ustilaginomycetes, Laboulbeniomycetes, Orbiliomycetes, Cercozoa_cls_Incertae_sedis and Wallemiomycetes in EZ (Figure 6a). At the genus level, higher abundances of Piromyces in BZ; Myrothecium, Aphanoascus, Ectophoma and Coniothyrium in CZ; Lophotrichus, Acremonium, Neurospora, Meyerozyma, Scheffersomyces and Xeromyces in DZ; and Cladosporium, Aspergillus and Wallemia in EZ were found (Figure 6b).
A LefSe analysis revealed that five and 34 observably higher abundances of phyla and genus, respectively, were observed among the sheep groups (Table 2, Figure 7). Furthermore, an analysis of the DESeq2 volcano map was performed, and it was found that compared with the AZ group, eight phyla and 56 genera were obviously different in the other ST sheep groups (Table 3, Figure 8). A network analysis showed that the phyla of Basidiomycota, Glomeromycota, Nematoda, Bryophyta and Marchantiophyta were positively related to the microbiota of the ST sheep, while Ascomycota, Mortierellomycota and Mucoromycota were negative phyla (Figure 9a). At the genus level, Piromyces, Caecomyces, Penicillium, Aspergillus, Wallemia, Naganishia, Sebacina, Mortierella, Inocybe, Cladosporium, Xeromyces, Talaromyces, Alternaria, Fusarium, Acremonium, Archaeorhizomyces, Tausonia, Trichoderma and Thermomyces were positively related to the fungi microbiota of the ST sheep, while Scheffersomyces, Orpinomyces, Neocallimastix and Pecoramyces were negative genera (Figure 9b).

3.4. Functional Analysis of Fungal Microbiota in the ST Sheep in Different Feeding Groups

Function comparison through a metaCys pathways analysis showed that there were 63 significantly different pathways among the five ST sheep groups (Figure 10). An enzyme abundance analysis found that there were 770 markedly different enzymes among the current ST sheep groups (Table S5).

4. Discussion

In the plateau region of China, food animals are considered to be important agricultural resources. Sheep are among the popular food animals, with more than 326 million heads that provided 5.25 million tons of mutton in China in 2022 [24]. Therefore, due to an increase in the demand for quality meat in China, efficient farming is of great importance. Gut microbiome analysis is considered to be an indicator of animal health [25]. Similar to bacteria in gut microbiota, fungi also play a vital role in homoeostasis, health status and assembling co-residing bacteria [9].
In this study, we found that the average weights of the ST sheep fed with concentrated feeds were higher, especially in the DZ group (Figure S1). After performing a fungal microbiome analysis of Suffolk crossed with Tibetan sheep with different feeds, 2,781,461 raw reads and 2,333,239 clean reads were obtained (Table 1). Compared with the sheep in groups AZ and BZ (164), the shared ASVs between AZ and CZ (109), AZ (113), DZ (118), and AZ and EZ were found to be fewer (Figure 1), which indicated a difference in the microbiota of the ST sheep fed with concentrated feed. No remarkable differences in the Alpha diversity index were found among the ST sheep groups, which is in accordance with Hu lambs fed with different diets [26] and dogs fed with different diets [27]. However, these results are not in line with those reported in Tibetan sheep fed with different forage types [28]. The current results and previous findings indicate that diet is one of the factors that may affect the richness and diversity of the microbiota. The ANOSIM analysis demonstrated prominent differences among the sheep groups (Figure S2) that were in accordance with the results in dogs fed with different diets [27] and yaks supplemented with dietaries [29]. These findings reveal that diet affects the intestinal microbiota. Structural differences in the microbiota among the five ST sheep groups were revealed through an analysis of a grouped percentage chart (Figure 4), a grouping and clustering heat map (Figure 5) and an evolutionary tree of species with a heat map (Figure 6). Moreover, different fungi species were examined through LEfSe (Figure 7) and a DESeq2 volcano map (Figure 8). The shifted fungi of the microbiota also affected its function (Figure 10 and Table S2).
In the ST sheep, the principal phyla in AZ (48.51%, 22.81%) and BZ (57.76%, 14.09%) were Neocallimastigomycota and Ascomycota, whereas Ascomycota and Neocallimastigomycota were the main phyla in CZ (33.32%, 22.24%). Ascomycota and Basidiomycota were the primary phyla in DZ (73.28%, 5.90%) and EZ (51.40%, 13.26%), respectively (Figure 4a). However, Ascomycota was the dominant phyla in all ST sheep, but its higher abundance was found in CZ, DZ and EZ. These results demonstrated that differences in feed altered the dominant fungi phyla in the sheep, as the fungal mycobiome was less stable and affected by the environment [30,31]. The current results are not in agreement with the results found in yaks using different feed models [32]. The phylum Ascomycota is commonly known due to its function in degrading lignin and keratin, which explains the higher weight gain in the ST sheep in the concentrated feed groups. Among the ST sheep groups eight phyla and 56 genera were examined and compared with free-range sheep in AZ (Figure 8). Among those genera, a higher abundance of Xeromyces was reported in healthy children compared with depressed people [33]. The higher abundance of this genus in the CZ, DZ and EZ groups indicates that the concentrated feed promoted ruminant health by increasing the abundance of those fungi. Kazachstania is related to epithelial glycolysis [34], and the higher abundance of this genera in CZ means that concentrated feed #1 promoted glycometabolism in the ST sheep. Species from Cordyceps are widely used in medicines that increase the abundance of useful gut bacteria [35,36]. The abundance of this genera in sheep in the CZ and DZ groups shows that concentrated feeds #1 and #2 promoted sheep health by regulating the intestinal microbiota. Rhodotorula produces nutrients like proteins and vitamins, which benefit host health [37]; the higher abundance of this genera in the DZ and EZ groups revealed that concentrated feeds #2 and #3 promoted the colonization of this useful microflora. Pichia and Sporisorium are two major genera in healthy humans [31,38], while higher amounts of Lecanicillium were found in healthy Tibetan pigs compared with diarrheal animals [39]; the higher abundances of these genera in the DZ groups means that concentrated feed #2 increased the abundance of beneficial fungi in the sheep. A higher abundance of Mortierella was detected in healthy yaks compared with ruminants infected with Cryptosporidium parvum [40]. A lower abundance of Plectosphaerella was examined in fluoride-induced mice [41,42], Tomentella was positively related to the generation of short-chain fatty acids [43,44,45] and Filobasidium reduced ulcerative colitis [46]; the higher abundance of this genus in the EZ group revealed that concentrated feed #3 promoted sheep health by mediating fungi abundance. The genera of Aphanoascus, Gibberella, Thelonectria and Sebacina are commonly reported in the environment [47,48,49,50]. However, these genera are not pathogenic, but the abundant changes could affect the microbiota function, which eventually affect host status and weight gains. So, these genera are commonly recognized as beneficial or dominant genera in healthy animals. Therefore, it is recommended that ST sheep may be fed with concentrated feed to increase weight gain through the regulation of fungi.

5. Conclusions

Based on the above, we concluded that differences in feed in Suffolk crossed with Tibetan sheep altered their gut microbiota. It was also found that different concentrations of phyla (eight) and genera (56) were observed among the different groups. Higher abundances of the phyla Ascomycota and the genera Kazachstania and Rhodotorula were found in the ST sheep fed with concentrated feed, which revealed that the concentrated feeds promoted weight gain in the sheep through the regulation of fungi. The results of this study may provide new insights regarding sheep farming on the plateau in China.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/life13112210/s1, Table S1: Nutritional ingredient information for the feeds used in the current study; Table S2: Detailed PCR information for this study; Table S3: Sequencing data generated in the current study; Table S4: Statistical analysis of the Alpha diversity index in the ST sheep groups; Table S5: Comparison of fungi microbiota function of enzyme abundances in ST sheep in different feeding groups; Figure S1. Body weights of the ST sheep in different groups. Significance is presented as * p < 0.05 and *** p < 0.001; data are presented as the mean ± SEM (n = 6); Figure S2. Comparing analyses of the Alpha diversity index in the ST sheep groups: (a) chao1; (b) faith_pd; (c) observed_features; (d) shannon_entropy; (e) Simpson. ns: non-significant; Figure S3. Pairwise ANOSIM analysis of Qiime2 in the five ST sheep groups.

Author Contributions

Y.R. and K.L.: research idea and methodology; Y.R., R.C., J.W. and K.L.: reagents, materials and analysis tools; K.L.: writing—original draft and preparation; K.M., F.S.A., Y.R. and K.L.: writing—review and editing; Y.R. and K.L.: visualization and supervision. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Major Projects for the Tibet Autonomous Region (XZ202101ZD0001N).

Institutional Review Board Statement

All the experiment procedures were conducted under the guidelines of the ethics committee of Nanjing Agricultural University (NJAU.No20220520108) approved on 20 May 2022.

Data Availability Statement

All raw sequence data were deposited in the NCBI Sequence Read Archive database under accession number PRJNA1008147.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cai, R.; Cui, X.; Zhang, S.; Xu, C. Effects of Regular Water Replenishment on Enzyme Activities and Fungal Metabolic Function of Sheep Manure Composting on the Qinghai–Tibet Plateau. Int. J. Environ. Res. Public Health 2022, 19, 12143. [Google Scholar] [CrossRef] [PubMed]
  2. Ge, Q.; Guo, Y.; Zheng, W.; Cai, Y.; Qi, X.; Zhao, S. A comparative analysis of differentially expressed mRNAs, miRNAs and circRNAs provides insights into the key genes involved in the high-altitude adaptation of yaks. BMC Genom. 2021, 22, 744. [Google Scholar] [CrossRef] [PubMed]
  3. Liu, J.; Lu, Z.; Yuan, C.; Wang, F.; Yang, B. Phylogeography and Phylogenetic Evolution in Tibetan Sheep Based on MT-CYB Sequences. Animals 2020, 10, 1177. [Google Scholar] [CrossRef] [PubMed]
  4. Jing, X.; Wang, W.; Degen, A.; Guo, Y.; Kang, J.; Liu, P.; Ding, L.; Shang, Z.; Fievez, V.; Zhou, J.; et al. Tibetan sheep have a high capacity to absorb and to regulate metabolism of SCFA in the rumen epithelium to adapt to low energy intake. Brit. J. Nutr. 2020, 123, 721–736. [Google Scholar] [CrossRef] [PubMed]
  5. Zhao, Y.; Zhao, E.; Zhang, N.; Duan, C. Mitochondrial DNA diversity, origin, and phylogenic relationships of three Chinese large-fat-tailed sheep breeds. Trop. Anim. Health Prod. 2011, 43, 1405–1410. [Google Scholar] [CrossRef]
  6. El-Sayed, A.; Aleya, L.; Kamel, M. Microbiota’s role in health and diseases. Environ. Sci. Pollut. Res. 2021, 28, 36967–36983. [Google Scholar] [CrossRef]
  7. Ottman, N.; Smidt, H.; de Vos, W.M.; Belzer, C. The function of our microbiota: Who is out there and what do they do? Front. Cell. Infect. Microbiol. 2012, 2, 104. [Google Scholar] [CrossRef]
  8. Lapiere, A.; Richard, M.L. Bacterial-fungal metabolic interactions within the microbiota and their potential relevance in human health and disease: A short review. Gut Microbes 2022, 14, 2105610. [Google Scholar] [CrossRef]
  9. Zhang, F.; Aschenbrenner, D.; Yoo, J.Y.; Zuo, T. The gut mycobiome in health, disease, and clinical applications in association with the gut bacterial microbiome assembly. Lancet Microbe 2022, 3, e969–e983. [Google Scholar] [CrossRef]
  10. Das, A.; Herlihy, E.O.; Shanahan, F.; Toole, P.W.O.; Jeffery, I.B. The fecal mycobiome in patients with irritable bowel syndrome. Sci. Rep. 2021, 11, 124. [Google Scholar] [CrossRef]
  11. Li, X.V.; Leonardi, I.; Putzel, G.G.; Semon, A.; Fiers, W.D.; Kusakabe, T.; Lin, W.; Gao, I.H.; Doron, I.; Gutierrez-Guerrero, A.; et al. Immune regulation by fungal strain diversity in inflammatory bowel disease. Nature 2022, 603, 672–678. [Google Scholar] [CrossRef] [PubMed]
  12. Zeng, S.; Schnabl, B. Roles for the mycobiome in liver disease. Liver Int. 2022, 42, 729–741. [Google Scholar] [CrossRef] [PubMed]
  13. Scarpellini, E.; Ianiro, G.; Attili, F.; Bassanelli, C.; De Santis, A.; Gasbarrini, A. The human gut microbiota and virome: Potential therapeutic implications. Dig. Liver Dis. 2015, 47, 1007–1012. [Google Scholar] [CrossRef]
  14. Anwar, H.; Iftikhar, A.; Muzaffar, H.; Almatroudi, A.; Allemailem, K.S.; Navaid, S.; Saleem, S.; Khurshid, M. Biodiversity of Gut Microbiota: Impact of Various Host and Environmental Factors. Biomed. Res. Int. 2021, 2021, 5575245. [Google Scholar] [CrossRef]
  15. Zhang, Y.K.; Zhang, X.X.; Li, F.D.; Li, C.; Li, G.Z.; Zhang, D.Y.; Song, Q.Z.; Li, X.L.; Zhao, Y.; Wang, W.M. Characterization of the rumen microbiota and its relationship with residual feed intake in sheep. Animal 2021, 15, 100161. [Google Scholar] [CrossRef] [PubMed]
  16. Kiarie, E.; Romero, L.F.; Nyachoti, C.M. The role of added feed enzymes in promoting gut health in swine and poultry. Nutr. Res. Rev. 2013, 26, 71–88. [Google Scholar] [CrossRef]
  17. Wang, C.; Wei, S.; Chen, N.; Xiang, Y.; Wang, Y.; Jin, M. Characteristics of gut microbiota in pigs with different breeds, growth periods and genders. Microb. Biotechnol. 2022, 15, 793–804. [Google Scholar] [CrossRef]
  18. Li, Z.; Zu, C.; Wang, C.; Yang, J.; Yu, H.; Wu, H. Different responses of rhizosphere and non-rhizosphere soil microbial communities to consecutive Piper nigrum L. monoculture. Sci. Rep. 2016, 6, 35825. [Google Scholar] [CrossRef]
  19. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef]
  20. Bokulich, N.A.; Kaehler, B.D.; Rideout, J.R.; Dillon, M.; Bolyen, E.; Knight, R.; Huttley, G.A.; Gregory Caporaso, J. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome 2018, 6, 90. [Google Scholar] [CrossRef]
  21. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef]
  22. Segata, N.; Izard, J.; Waldron, L.; Gevers, D.; Miropolsky, L.; Garrett, W.S.; Huttenhower, C. Metagenomic biomarker discovery and explanation. Genome Biol. 2011, 12, R60. [Google Scholar] [CrossRef]
  23. Langille, M.G.I.; Zaneveld, J.; Caporaso, J.G.; McDonald, D.; Knights, D.; Reyes, J.A.; Clemente, J.C.; Burkepile, D.E.; Vega Thurber, R.L.; Knight, R.; et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 2013, 31, 814–821. [Google Scholar] [CrossRef]
  24. Available online: https://data.stats.gov.cn/easyquery.htm?cn=C01 (accessed on 2 November 2023).
  25. Gao, B.; Chi, L.; Zhu, Y.; Shi, X.; Tu, P.; Li, B.; Yin, J.; Gao, N.; Shen, W.; Schnabl, B. An Introduction to Next Generation Sequencing Bioinformatic Analysis in Gut Microbiome Studies. Biomolecules 2021, 11, 530. [Google Scholar] [CrossRef]
  26. Xiong, Y.; Wang, X.; Li, X.; Guo, L.; Yang, F.; Ni, K. Exploring the rumen microbiota of Hu lambs in response to diet with paper mulberry. Appl. Microbiol. Biotechnol. 2023, 107, 4961–4971. [Google Scholar] [CrossRef] [PubMed]
  27. Schmidt, M.; Unterer, S.; Suchodolski, J.S.; Honneffer, J.B.; Guard, B.C.; Lidbury, J.A.; Steiner, J.M.; Fritz, J.; Kölle, P. The fecal microbiome and metabolome differs between dogs fed Bones and Raw Food (BARF) diets and dogs fed commercial diets. PLoS ONE 2018, 13, e201279. [Google Scholar] [CrossRef] [PubMed]
  28. Cui, X.; Wang, Z.; Yan, T.; Chang, S.; Wang, H.; Hou, F. Rumen bacterial diversity of Tibetan sheep (Ovis aries) associated with different forage types on the Qinghai-Tibetan Plateau. Can. J. Microbiol. 2019, 65, 859–869. [Google Scholar] [CrossRef] [PubMed]
  29. Ahmad, A.A.; Zhang, J.B.; Liang, Z.; Yang, C.; Kalwar, Q.; Shah, T.; Du, M.; Muhammad, I.; Zheng, J.; Yan, P.; et al. Dynamics of rumen bacterial composition of yak (Bos grunniens) in response to dietary supplements during the cold season. PeerJ 2021, 9, e11520. [Google Scholar] [CrossRef] [PubMed]
  30. Hallen-Adams, H.E.; Suhr, M.J. Fungi in the healthy human gastrointestinal tract. Virulence 2016, 8, 352–358. [Google Scholar] [CrossRef]
  31. Szóstak, N.; Handschuh, L.; Samelak-Czajka, A.; Tomela, K.; Schmidt, M.; Pruss, Ł.; Milanowska-Zabel, K.; Kozlowski, P.; Philips, A.; Sangwan, N. Host Factors Associated with Gut Mycobiome Structure. mSystems 2023, 8, e98622. [Google Scholar] [CrossRef]
  32. Zhu, Y.; Cidan, Y.; Sun, G.; Li, X.; Shahid, M.A.; Luosang, Z.; Suolang, Z.; Suo, L.; Basang, W. Comparative analysis of gut fungal composition and structure of the yaks under different feeding models. Front. Vet. Sci. 2023, 10, 1193558. [Google Scholar] [CrossRef] [PubMed]
  33. Hao, S.; Zhang, Z.; Zhou, Y.; Zhang, X.; Sun, W.; Yang, Z.; Zhao, J.; Jiang, H. Altered gut bacterial–fungal interkingdom networks in children and adolescents with depression. J. Affect. Disord. 2023, 332, 64–71. [Google Scholar] [CrossRef] [PubMed]
  34. Hu, J.; Chen, J.; Hou, Q.; Xu, X.; Ren, J.; Ma, L.; Yan, X. Core-predominant gut fungus Kazachstania slooffiae promotes intestinal epithelial glycolysis via lysine desuccinylation in pigs. Microbiome 2023, 11, 31. [Google Scholar] [CrossRef] [PubMed]
  35. Tian, J.; Wang, X.; Zhang, X.; Chen, X.; Rui, X.; Zhang, Q.; Dong, M.; Li, W. Simulated digestion and fecal fermentation behaviors of exopolysaccharides from Paecilomyces cicadae TJJ1213 and its effects on human gut microbiota. Int. J. Biol. Macromol. 2021, 188, 833–843. [Google Scholar] [CrossRef] [PubMed]
  36. Tian, J.; Zhao, X.; Tang, C.; Wang, X.; Zhang, X.; Xiao, L.; Li, W. Protective effect of Paecilomyces cicadae TJJ11213 exopolysaccharide on intestinal mucosa and regulation of gut microbiota in immunosuppressed mice. Food Res. Int. 2023, 165, 112477. [Google Scholar] [CrossRef]
  37. Hof, H. Rhodotorula spp. in the gut—Foe or friend? GMS Infect. Dis. 2019, 7, 6734584. [Google Scholar] [CrossRef]
  38. Nash, A.K.; Auchtung, T.A.; Wong, M.C.; Smith, D.P.; Gesell, J.R.; Ross, M.C.; Stewart, C.J.; Metcalf, G.A.; Muzny, D.M.; Gibbs, R.A.; et al. The gut mycobiome of the Human Microbiome Project healthy cohort. Microbiome 2017, 5, 153. [Google Scholar] [CrossRef]
  39. Kong, Q.; Liu, S.; Li, A.; Wang, Y.; Zhang, L.; Iqbal, M.; Jamil, T.; Shang, Z.; Suo, L.; Li, J. Characterization of fungal microbial diversity in healthy and diarrheal Tibetan piglets. BMC Microbiol. 2021, 21, 204. [Google Scholar] [CrossRef]
  40. Lu, S.; Zou, W.; Chen, X.; Sun, G.; Cidan, Y.; Almutairi, M.H.; Dunzhu, L.; Nazar, M.; Mehmood, K.; Zhu, Y.; et al. Effects of Cryptosporidium parvum infection on intestinal fungal microbiota in yaks (Bos grunniens). Microb. Pathog. 2023, 183, 106322. [Google Scholar] [CrossRef]
  41. Cao, Q.; Li, R.; Fu, R.; Zhang, X.; Yue, B.; Wang, J.; Sun, Z.; Niu, R. Intestinal fungal dysbiosis in mice induced by fluoride. Chemosphere 2020, 245, 125617. [Google Scholar] [CrossRef]
  42. Li, J.; Chen, D.; Yu, B.; He, J.; Huang, Z.; Mao, X.; Zheng, P.; Yu, J.; Luo, J.; Tian, G.; et al. The fungal community and its interaction with the concentration of short-chain fatty acids in the faeces of Chenghua, Yorkshire and Tibetan pigs. Microb. Biotechnol. 2019, 13, 509–521. [Google Scholar] [CrossRef] [PubMed]
  43. Zaki, M.G.; Barka, T.A.; Tayeb, F.A.E.F. Effect of ginger powder (Zingiber officinale) on acid-base balance, rumen and blood constituents in healthy Egyptian sheep. Int. J. Vet. Sci. 2021, 10, 55–58. [Google Scholar] [CrossRef]
  44. Oraby, M.I.; Baraka, T.A.; Rakha, G.H. Impact of cadmium intoxication on health status, rumen and blood constituents in Egyptian Ossimi sheep. Int. J. Vet. Sci. 2021, 10, 102–106. [Google Scholar] [CrossRef]
  45. Naibaho, F.G.; Hartanto, A.; Bintang, M.; Jamilah, I.; Priyani, N.; Putra, E.D. GC-MS analysis antimicrobial activity of the aqueous extract from the bulbs of Allium chinense G. Don. cultivated in North Sumatra, Indonesia. Asian J. Agric. Biol. 2021, 2, 201912562. [Google Scholar] [CrossRef]
  46. Van Thiel, I.A.M.; Rahman, S.; Hakvoort, T.B.M.; Davids, M.; Verseijden, C.; van Hamersveld, P.H.P.; Bénard, M.V.; Lodders, M.H.; Boekhout, T.; van den Wijngaard, R.M.; et al. Fecal Filobasidium Is Associated with Clinical Remission and Endoscopic Response following Fecal Microbiota Transplantation in Mild-to-Moderate Ulcerative Colitis. Microorganisms 2022, 10, 737. [Google Scholar] [CrossRef] [PubMed]
  47. Novak Babič, M.; Gunde-Cimerman, N.; Breskvar, M.; Džeroski, S.; Brandão, J. Occurrence, Diversity and Anti-Fungal Resistance of Fungi in Sand of an Urban Beach in Slovenia—Environmental Monitoring with Possible Health Risk Implications. J. Fungi 2022, 8, 860. [Google Scholar] [CrossRef]
  48. Undugoda, L.; Kannangara, S. Nature and activities of microfungi associated with the decomposition of rice straw in Sri Lanka. Asian J. Agric. Biol. 2022, 2022, 202103115. [Google Scholar] [CrossRef]
  49. Chafai, W.; El Gabardi, S.; Douira, A.; Khalid, A. Diversity and mycorrhizal potential of arbuscular mycorrhizal fungi in two natural soils in the eastern region of Morocco. Asian J. Agric. Biol. 2022, 2022, 202102101. [Google Scholar] [CrossRef]
  50. He, M.; Sun, W.; Cui, S.; Mu, G.; Liu, L.; Guo, W. Analysis of Microbial Diversity and Community Structure of Peanut Pod and Its Surrounding Soil in Peanut Rot Epidemic Area. Curr. Microbiol. 2021, 78, 2173–2182. [Google Scholar] [CrossRef]
Figure 1. Sequence data analysis: (a) Venn map; (b) annotation statistics. Groups: AZ, BZ, CZ, DZ and EZ. In (b), the different colors represent different taxa.
Figure 1. Sequence data analysis: (a) Venn map; (b) annotation statistics. Groups: AZ, BZ, CZ, DZ and EZ. In (b), the different colors represent different taxa.
Life 13 02210 g001
Figure 2. Alpha rarefaction curves for the ST sheep groups: (a) faith_pd; (b) observed_features; (c) shannon_entropy. The different colors represent the different groups of AZ, BZ, CZ, DZ and EZ.
Figure 2. Alpha rarefaction curves for the ST sheep groups: (a) faith_pd; (b) observed_features; (c) shannon_entropy. The different colors represent the different groups of AZ, BZ, CZ, DZ and EZ.
Life 13 02210 g002
Figure 3. Comparing analyses of the beta diversity index in the ST sheep groups: (a) PCA; (b) PCoA; (c) NMDS; (d) Qiime2. The different colors represent the different groups of AZ, BZ, CZ, DZ and EZ.
Figure 3. Comparing analyses of the beta diversity index in the ST sheep groups: (a) PCA; (b) PCoA; (c) NMDS; (d) Qiime2. The different colors represent the different groups of AZ, BZ, CZ, DZ and EZ.
Life 13 02210 g003
Figure 4. The microbiota structure of the ST sheep in different taxa: (a) phylum; (b) class; (c) order; (d) family; (e) genus. Groups: AZ, BZ, CZ, DZ and EZ. The different colors represent different species.
Figure 4. The microbiota structure of the ST sheep in different taxa: (a) phylum; (b) class; (c) order; (d) family; (e) genus. Groups: AZ, BZ, CZ, DZ and EZ. The different colors represent different species.
Life 13 02210 g004
Figure 5. Comparison of the microbiota structures of the ST sheep in different taxa via grouping and a clustering heat map: (a) phylum; (b) class; (c) order; (d) family; (e) genus.
Figure 5. Comparison of the microbiota structures of the ST sheep in different taxa via grouping and a clustering heat map: (a) phylum; (b) class; (c) order; (d) family; (e) genus.
Life 13 02210 g005
Figure 6. Comparison of the microbiota structure of the ST sheep in different taxa via an evolutionary tree of species with a heat map: (a) class; (b) genus.
Figure 6. Comparison of the microbiota structure of the ST sheep in different taxa via an evolutionary tree of species with a heat map: (a) class; (b) genus.
Life 13 02210 g006
Figure 7. Revealing significant different fungi species among the five ST sheep groups in different taxa via LEfSe. (a) Phylum, (b) Genera.
Figure 7. Revealing significant different fungi species among the five ST sheep groups in different taxa via LEfSe. (a) Phylum, (b) Genera.
Life 13 02210 g007
Figure 8. Significant different fungi species among the five ST sheep groups for different taxa according to the DESeq2 volcano map: (a) phylum; (b) genus.
Figure 8. Significant different fungi species among the five ST sheep groups for different taxa according to the DESeq2 volcano map: (a) phylum; (b) genus.
Life 13 02210 g008
Figure 9. Network analysis of fungi microbiome of the ST sheep for different taxa: (a) phylum; (b) genus.
Figure 9. Network analysis of fungi microbiome of the ST sheep for different taxa: (a) phylum; (b) genus.
Life 13 02210 g009
Figure 10. Comparison of the fungi microbiota function of metaCys pathways in the ST sheep for different taxa. Different lowercase letters indicate significant difference.
Figure 10. Comparison of the fungi microbiota function of metaCys pathways in the ST sheep for different taxa. Different lowercase letters indicate significant difference.
Life 13 02210 g010
Table 1. The main species in different ST sheep groups for different taxa.
Table 1. The main species in different ST sheep groups for different taxa.
TaxonomyGroupSpeciesPercent (%)
PhylumAZAscomycota22.81%
AZNeocallimastigomycota48.51%
AZChytridiomycota7.10%
BZAscomycota14.09%
BZNeocallimastigomycota57.76%
BZChytridiomycota4.55%
CZAscomycota33.32%
CZNeocallimastigomycota22.24%
CZBasidiomycota14.57%
DZAscomycota73.28%
DZNeocallimastigomycota4.29%
DZBasidiomycota5.90%
EZAscomycota51.40%
EZNeocallimastigomycota4.29%
EZBasidiomycota13.26%
ClassAZNeocallimastigomycetes48.51%
AZSaccharomycetes16.30%
AZSordariomycetes1.36%
BZNeocallimastigomycetes57.76%
BZSaccharomycetes8.50%
BZSordariomycetes1.23%
CZNeocallimastigomycetes22.24%
CZEurotiomycetes 9.73%
CZSordariomycetes 9.50%
DZSaccharomycetes44.07%
DZEurotiomycetes9.43%
DZSordariomycetes7.43%
EZSaccharomycetes22.15%
EZEurotiomycetes11.17%
EZSordariomycetes8.11%
OrderAZNeocallimastigales48.51%
AZSaccharomycetales16.30%
AZEurotiales 9.60%
BZNeocallimastigales57.76%
BZSaccharomycetales8.50%
BZPleosporales 0.46%
CZNeocallimastigales22.24%
CZFilobasidiales9.15%
CZEurotiales8.79%
DZSaccharomycetales 44.07%
DZEurotiales9.08%
DZHypocreales3.40%
EZSaccharomycetales22.15%
EZEurotiales10.84%
EZWallemiales8.59%
FamilyAZNeocallimastigaceae48.51%
AZDebaryomycetaceae16.28%
AZAspergillaceae0.59%
BZNeocallimastigaceae 57.76%
BZDebaryomycetaceae8.02%
BZPhaffomycetaceae0.37%
CZNeocallimastigaceae22.24%
CZFilobasidiaceae9.15%
CZAspergillaceae7.25%
DZDebaryomycetaceae43.14%
DZAspergillaceae7.57%
DZWallemiaceae1.63%
EZDebaryomycetaceae21.76%
EZAspergillaceae10.26%
EZWallemiaceae8.59%
GeneraAZPiromyces18.77%
AZScheffersomyces16.12%
AZCaecomyces9.65%
BZScheffersomyces7.97%
BZOrpinomyces4.14%
BZPiromyces4.12%
CZCaecomyces10.57%
CZNaganishia9.04%
CZAspergillus5.29%
DZScheffersomyces42.80%
DZAspergillus3.26%
DZXeromyces2.85%
EZScheffersomyces21.55%
EZWallemia8.59%
EZAspergillus6.12%
Table 2. Significant different species among the sheep groups according to LefSe.
Table 2. Significant different species among the sheep groups according to LefSe.
TaxonomyGroupSpeciesp-ValueSignificance
PhylumAZChytridiomycota<0.01
BZNeocallimastigomycota<0.01
DZAscomycota<0.01
EZBasidiomycota<0.05
EZMortierellomycota<0.05
GeneraAZOrpinomyces<0.01
AZNeocallimastix<0.01
BZNeocallimastigaceae<0.01
BZNeocallimastigomycetes<0.01
BZNeocallimastigales<0.01
BZNeocallimastigomycota<0.01
BZPiromyces<0.001
CZSebacina<0.05
CZArchaeorhizomycetales<0.05
CZBasidiomycota<0.05
CZArchaeorhizomyces<0.05
CZArchaeorhizomycetes<0.05
CZSebacinales<0.05
CZAphanoascus<0.05
CZOnygenaceae<0.05
CZArchaeorhizomycetaceae<0.05
CZNaganishia<0.05
CZFilobasidiales<0.05
CZSebacinaceae<0.05
CZFilobasidiaceae<0.05
CZOnygenales<0.05
DZAscomycota<0.01
DZAgaricales<0.05
DZLophotrichus<0.05
DZXeromyces<0.05
EZEurotiales<0.05
EZEurotiomycetes<0.05
EZWallemiales<0.05
EZWallemia<0.05
EZPenicillium<0.05
EZAspergillaceae<0.05
EZWallemiomycetes<0.05
EZWallemiaceae<0.05
EZAgaricomycetes<0.01
Table 3. Significant different species among the sheep groups according to DESeq2.
Table 3. Significant different species among the sheep groups according to DESeq2.
TaxonomyGroupSpeciesp-ValueSignificance
PhylumCZBasidiomycota<0.01
CZMucoromycota <0.01
CZChytridiomycota<0.0001
DZMortierellomycota<0.0001
DZMucoromycota<0.05
DZAscomycota<0.05
DZRozellomycota<0.05
DZChytridiomycota<0.0001
DZNeocallimastigomycota<0.0001
DZStreptophyta<0.05
EZBasidiomycota <0.001
EZMortierellomycota<0.001
EZChytridiomycota<0.001
EZChytridiomycota<0.01
GeneraCZPhaeosphaeria<0.0001
CZAphanoascus<0.0001
CZGibberella<0.0001
CZPlenodomus<0.0001
CZThelonectria<0.0001
CZMucor<0.001
CZSebacina<0.001
CZXeromyces<0.001
CZKazachstania<0.001
CZCordyceps<0.001
CZChrysosporium<0.001
CZBjerkandera<0.001
CZLophotrichus<0.001
CZNaganishia<0.001
CZArchaeorhizomyces<0.001
CZCephaliophora<0.01
CZScopulariopsis<0.01
CZClavispora<0.01
CZAmphinema<0.01
CZFusarium<0.01
CZLophiostoma<0.01
CZPenicillium<0.05
CZPhoma<0.05
CZAspergillus<0.05
CZPecoramyces<0.001
CZOccultifur<0.001
CZTuber<0.01
CZPapiliotrema<0.05
CZSaitozyma<0.05
DZXeromyces<0.0001
DZHyphopichia<0.0001
DZLophotrichus<0.0001
DZPichia<0.0001
DZSebacina<0.001
DZChrysosporium<0.001
DZPenicillium<0.001
DZClonostachys<0.01
DZRhodotorula<0.01
DZCortinarius<0.01
DZSporisorium<0.01
DZFusarium<0.01
DZAmphinema<0.01
DZClavispora<0.01
DZVerticillium<0.01
DZClavulina<0.05
DZIlyonectria<0.05
DZCordyceps<0.05
DZSetophoma<0.05
DZAcremonium<0.05
DZLecanicillium<0.05
DZPecoramyces<0.0001
DZCaecomyces<0.0001
DZNeocallimastix<0.0001
DZPiromyces<0.0001
DZOccultifur<0.0001
DZStephanonectria<0.001
DZPapiliotrema<0.01
DZTuber<0.05
DZColletotrichum<0.05
DZSaitozyma<0.05
EZXeromyces<0.0001
EZClavispora<0.0001
EZPenicillium<0.0001
EZMortierella<0.0001
EZPlectosphaerella<0.0001
EZGibberella<0.0001
EZRussula<0.0001
EZWallemia<0.0001
EZSebacina<0.001
EZArchaeorhizomyces<0.001
EZHirsutella<0.001
EZTomentella<0.001
EZAphanoascus<0.001
EZPseudallescheria<0.001
EZRhodotorula<0.01
EZVerticillium<0.01
EZLophiostoma<0.01
EZFilobasidium<0.05
EZAspergillus<0.05
EZThelonectria<0.05
EZIlyonectria<0.05
EZInocybe<0.05
EZNeocallimastix<0.0001
EZOccultifur<0.001
EZStephanonectria<0.001
EZPhlebiopsis<0.01
EZSetophoma<0.05
EZXenopolyscytalum<0.05
EZTuber<0.05
EZZopfiella<0.05
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ren, Y.; Ciwang, R.; Wang, J.; Mehmood, K.; Ataya, F.S.; Li, K. Effect of Different Feeds on the Fungi Microbiome of Suffolk Crossed with Tibetan Sheep. Life 2023, 13, 2210. https://doi.org/10.3390/life13112210

AMA Style

Ren Y, Ciwang R, Wang J, Mehmood K, Ataya FS, Li K. Effect of Different Feeds on the Fungi Microbiome of Suffolk Crossed with Tibetan Sheep. Life. 2023; 13(11):2210. https://doi.org/10.3390/life13112210

Chicago/Turabian Style

Ren, Yue, Renzeng Ciwang, Jia Wang, Khalid Mehmood, Farid Shokry Ataya, and Kun Li. 2023. "Effect of Different Feeds on the Fungi Microbiome of Suffolk Crossed with Tibetan Sheep" Life 13, no. 11: 2210. https://doi.org/10.3390/life13112210

APA Style

Ren, Y., Ciwang, R., Wang, J., Mehmood, K., Ataya, F. S., & Li, K. (2023). Effect of Different Feeds on the Fungi Microbiome of Suffolk Crossed with Tibetan Sheep. Life, 13(11), 2210. https://doi.org/10.3390/life13112210

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop