Skip to main content

Diversity and functional prediction of fungal communities in different segments of mongolian horse gastrointestinal tracts

Abstract

Background

Anaerobic fungi are effective fibre-degrading microorganisms in the digestive tract of horses. However, our understanding of their diversity and community structure is limited, especially in different parts of the gastrointestinal tract.

Results

For the first time, high-throughput sequencing technology was used to analyse and predict fungal microbial diversity in different parts of the gastrointestinal tract of Mongolian horses. The results revealed that the richness and diversity of fungi in the hindgut of Mongolian horses were much higher than those in the foregut. The foregut was dominated by Basidiomycota and Ascomycota, whereas the hindgut was dominated by Neocallimastigomycota and Basidiomycota. At the genus level, the relative abundance of many pathogenic fungi (Cryptococcus, Cladosporium, Alternaria, and Sarocladium) in the foregut was significantly higher than that in the posterior gut, indicating that Mongolian horses have strong disease resistance. The prediction of fungal function also showed significant differences in the fungal flora between the foregut and the hindgut. The fungi in Mongolian horses’ foreguts were mainly pathologically nutritive and contained many animal and plant pathogens, particularly in the small intestine (jejunum and ileum). This indicates that the foregut may be the most important immune site in the digestive system of Mongolian horses, which explains the high disease resistance of Mongolian horses. The number of unassigned functional groups in the posterior gut was significantly higher than that in the anterior gut, indicating that the functions of fungal groups in the posterior gut have not been fully explored, and further studies are required in the future.

Conclusions

Analysis of high-throughput sequencing results revealed that the fungal composition varied greatly among different gastrointestinal tract segments in Mongolian horses, whose hindgut contains many anaerobic fungi involved in plant cellulose degradation. This provides important basic data for studying fungal diversity in the digestive system of healthy horses, which can be used for the health assessment of horses and provides clues for further research on the disease resistance and digestive capacity of horses, as well as a reference for the early diagnosis of intestinal diseases and innovative treatment methods.

Peer Review reports

Background

The gastrointestinal tract (GIT) microbiota not only provides the host with nutrients such as amino acids, vitamins, and short-chain fatty acids that are necessary for growth but also plays important roles in digestion, absorption, metabolism, and the immune system [1, 2]. However, the stability of the GIT microbiota is easily influenced by intrinsic (host) and extrinsic (environmental) factors [3], such as species [4], age [5], geographic location [6], transportation [7, 8], exercise stress [8], diet [9], and season [10]. Impairment or dysregulation of the equine GIT and its microbiome may lead to pathological responses [10,11,12]. Studies have revealed that diarrhoea in horses is closely related to changes in the composition and structure of intestinal microbiota [13, 14]. In addition, alterations in the gut microbiome may lead to various diseases, such as obesity [15], laminitis [16], equine glandular gastric disease (EGGD) [17], colic [18], and colitis [19]. Consequently, it is essential to determine the gut microbiota composition in healthy horses to determine the effects of metabolic disorders and diseases.

Horses are classified as monogastric herbivores, and they have a specific GIT tissue, and a strong posterior intestinal system. Studies have demonstrated that the equine GIT can be divided into two main regions: the upper GIT (UG) and lower GIT (LG), and each segment of the GIT has a different ecosystem. The bacterial microbiota composition of adjacent compartments is more similar than that of far-apart compartments [20]. The microbiota structure of equine UG (stomach, jejunum, and ileum) is more complex and diverse because of the influence of many environmental bacteria in forage. In contrast, the microbiota composition of the lower equine gut (caecum and colon) appears to be relatively stable, despite variables such as individual history, breed, or age [20,21,22]. The equine GIT contains various microbial communities, including fungi, parasites, protozoa, archaea, viruses, and bacteria [23]. However, knowledge about the role of fungi in the GIT and their contribution to the microbiome of healthy horses is limited, as reflected in only a few studies on hindgut anaerobic fungi [24]. Anaerobic fungi (Neocallimastigomycota) are widely found in the hindgut of horses and play an important role in degrading plant cellulose, given that they possess a complete and highly effective group of plant cell wall-degrading enzymes [25]. To date, anaerobic fungi isolation and morphological identification mainly rely on purification culture and microscopy. Nine genera have been identified, and more than 20 species have been described [26]. Studies have indicated that the digestive tract of horses is dominated by novel, as-yet-uncultured anaerobic fungi that are different from previously described foregut herbivores [24].

In horses, the gut microbiota has been used as a molecular marker for bacteria and archaea, mainly by DNA sequencing, targeting variable regions of the 16 S rRNA gene. In contrast, eukaryotic communities (e.g., anaerobic fungi) have been described in the equine hindgut and have been identified to play an important role in digestion [11, 13]. Next-generation sequencing (NGS) technology has been applied to analyse fungal diversity copies in equine gut contents and faecal samples, but relevant research remains far from adequate [14, 24]. Indeed, stool reflects microbial changes only in the distal portion of the hindgut and not the entire GI microbiota. Therefore, evaluating the composition of anaerobic fungal communities in different parts of the gastrointestinal tract is important. This study aimed to characterize and compare the fungal composition in different parts of the gastrointestinal tract of Mongolian horses using NGS technology.

Materials and methods

Horse and sample collection

GIT content samples were obtained from five healthy Mongolian horses (three males and two females, mean age 4.4 years, 3–6 years, weight 292.8 ± 11.9 kg) that were euthanized by overdose of pentobarbital sodium (Table 1). As described in a previous study [20], these horses were fed on the same pasture in the grasslands of the Xilin Gol League, Inner Mongolia, China, maintaining the same grazing conditions without any history of intestinal diseases. The gastrointestinal sample collection sites included the stomach (S), jejunum (J), and ileum (I) for the upper GIT and the caecum (C), ventral colon (VC), and dorsal colon (DC) for the lower GIT. The intestinal segments were tied to prevent mixing between adjacent segments, and the contents were collected in the middle part of each segment after horizontal placement. Samples were collected in 50 mL sterile enzyme-free centrifuge tubes, snap-frozen in liquid nitrogen, and stored at − 80 °C for further study.

Table 1 Specific details of the experimental Mongolian horses

DNA extraction, PCR amplification, and library construction

Genomic DNA was extracted from the GIT samples by the cetyltrimethylammonium bromide (CTAB) method [27], and the concentration and purity of DNA were detected by 1% agarose gel electrophoresis. DNA diluted at 1 ng/µL was used as a template for PCR amplification using specific ITS1 region primers: 1737 F (5’-GGAAGTAAAAGTCGTAACAAGG-3’) and 2043 R (5’-GCTGCGTTCTTCATCGATGC-3’) [28, 29]. Fungal genomic library construction and on-machine sequencing were completed by Beijing Nuohe Zhiyuan Technology Co., Ltd. An Ion Plus Fragment Library Kit (No: 4471252, Thermo, USA) was used to construct the library. After Qubit quantification and library qualification, the constructed library was sequenced by the Ion S5TMXL platform.

Bioinformatics analysis

Sequencing data processing

Cutadapt (V1.9.1, http://cutadapt.readthedocs.io/en/stable/) was used for preliminary quality control to obtain raw reads, which were compared with the Unite database (https://unite.ut.ee/) to remove chimeric sequences and obtain clean reads [30,31,32].

OTU clustering and species annotation

Clean reads were clustered using UPARSE software (v7.0.1001, http://drive5.com/uparse/), and sequences were clustered into operational taxonomic units (OTUs) based on a 97% identity threshold to obtain the number of OTUs [33]. The OTU representative sequences were analysed by species annotation using the BLAST method in QIIME software (Version 1.9.1, http://qiime.org/scripts/assign_taxonomy.html) and the UNITE database [34, 35]. Finally, the data of each sample were homogenized for subsequent alpha and beta diversity analysis.

Alpha diversity analysis

QIIME software (Version 1.9.1) was used to calculate the ACE, Chao1, Shannon, observed-species, and Goods-coverage indices. R software (Version 2.15.3) was used to draw the dilution curve, rank abundance curve, and species accumulation curve and to analyse the differences in the alpha diversity index between groups. Nonmetric multidimensional scaling (NMDS) analysis was performed with the vegan package using R software (Version 2.15.3).

Beta diversity analysis

QIIME software (Version 1.9.1) was used to calculate the UniFrac distance and construct the UPGMA clustering tree. Linear discriminant analysis effect size (LEfSe) analysis was conducted using LEfSe software, and the default screening value of the LDA score was set to 4. LEfSe analysis was used to analyse different species between groups. FunGuild software was used to predict the functional composition of the fungal communities [36].

Data analysis

All parameters of different segments of the equine GIT are expressed as the mean ± standard deviation. Statistical significance was analysed using ANOVA in SPSS 22.0 software, and multiple comparisons were performed by the LSD test. A P value < 0.05 indicated a significant difference, and a P value < 0.01 indicated a highly significant difference.

Results

Sequencing results were clustered with OTUs

In this study, the fungal diversity in the GIT of Mongolian horses was completely consistent with that previously reported, which provided a reference for further analysis [20]. After the sequencing data of 30 content samples from different compartments of equine GIT were spliced, filtered, and screened, the average number of valid sequences was 74,477, and the average number of reads of constructed OTUs with annotation information was 58,739. The average number of OTUs was 525 (Additional File 1). Dilution curve analysis showed that the sequencing depth of all samples in different sections almost reached a plateau at the OTU level, indicating that the current sequencing depth covered most microorganisms, and the amount of undetected data had little impact on subsequent analysis. This curve also indirectly reflected the species richness in different segments of the GI tract; the cecum (C) exhibited greater species richness than other parts (Fig. 1A). The rank abundance curve can objectively reflect the richness and evenness of species in the sample [37], and the results showed that the species richness and evenness were similar in each segment of the GIT at the OTU level (Fig. 1B).

Fig. 1
figure 1

Microbial community diversity in the gastrointestinal tract of Mongolian horses. A: Dilution curve; B: Rank abundance curve; C: OTU Venn diagram of gastrointestinal flora. In each figure, S represents the stomach segment, J represents the jejunum segment, I represents the ileum segment, C represents the caecum segment, VC represents the ventral colon segment, and DC represents the dorsal colon segment

The common and unique OTUs among different GIT segments were analysed based on the results of OTU cluster analysis. According to the petal diagram between different gastrointestinal segments, 359 OTUs were shared by the six groups of samples, which could be divided into six phyla (Fig. 1C, Additional File 2 and 3). A Venn plot of the lower GIT showed that the proportion of unique OTUs in the caecum, ventral colon, and dorsal colon in the total OTUs was 30.11%, 14.06%, and 10.20%, respectively, and the number of OTUs accounted for 27.11% of the total number of OTUs (Additional File 4).

Alpha diversity analysis

At a sequencing depth of 31,744 reads per sample, the alpha diversity indices (ACE, Chao1, Shannon, observed-species, Good’s coverage) were calculated to analyse the diversity and richness of microbial communities within the samples (Table 2). The richness indices were significantly higher (P < 0.05) in the caecum and ventral colon groups than in the small intestine group (jejunum and ileum) (observed species: 396.60 ± 56.35 (I) ~ 530.00 ± 116.83 (C), Chao1: 456.06 ± 77.84 (I) ~ 592.96 ± 130.23 (C), ACE: 458.39 ± 68.45 (I) ~ 600.33 ± 122.78 (C)). There was no significant difference in the species diversity index (Shannon: 5.17 ± 0.50 (I) ~ 5.28 ± 0.47 (C)). It is widely acknowledged that the Good’s coverage index can represent the coverage of the test results. The larger the index, the higher the coverage, and the more it represents the actual situation in the sample. The coverage index of fungal flora in all samples was above 99.6%, and the dilution curve tended to be flat, indicating that the amount of data was reasonable and could cover most species in the samples.

Table 2 Comparison of microbial diversity in different digestive parts of Mongolian horses

Analysis of fungal composition structure based on OTUs

According to the OTU species annotation results, Basidiomycota and Ascomycota had the highest proportion in the foregut. The fungal composition of the hindgut contents changed at the boundary between the small and large intestines, with Neocallimastigomycota and Basidiomycota as the dominant phyla. There were significant differences in the abundance of Neocallimastigomycota and Basidiomycota in the foregut and hindgut (P < 0.01). The relative abundance of Ascomycota significantly differed in the foregut and hindgut (P < 0.05) (Fig. 2A; Additional File 2). The relative abundance of microbial species in different gastrointestinal segments differed at the phylum level. Neocallimastigomycota was more abundant in dorsal colon than in the stomach, jejunum, ileum, and caecum (P < 0.01, < 0.01, < 0.01, < 0.05, respectively). Basidiomycota was more abundant in the stomach than in the caecum, ventral colon, and dorsal colon (P = 0.01, < 0.01 and < 0.01, respectively). Ascomycota was more abundant in the ileum than in the stomach, ventral colon, and dorsal colon (P = 0.01, < 0.01, < 0.01, respectively). Chytridiomycota was more abundant in the jejunum than in other parts (all P < 0.01) (Fig. 2A; Additional File 3).

Fig. 2
figure 2

Relative abundance in luminal Mongolian horses GIT fungal composition. A: On the left is the UPGMA clustering tree structure, and on the right is the relative species abundance distribution of each sample at the phylum level; B: Subject level; C: Genus level. In each figure, S represents the stomach segment, J represents the jejunum segment, I represents the ileum segment, C represents the caecum segment, VC represents the ventral colon segment, and DC represents the dorsal colon segment

At the family level, the relative abundance of Neocallimastigaceae in the hindgut was significantly higher than that in the foregut (P < 0.01) (Fig. 2B). At the genus level, the fungal composition of luminal contents changed significantly during the transition from the small intestine to the large intestine. Although the abundance of some fungal genera was high in the foregut, no single genus exceeded 15% in any foregut samples (Fig. 2C; Additional File 5). The genera with the most abundant genera detected in the fungal community composition in the hindgut (relative abundance > 1%) included Cryptococcus (2.34%), Cladosporium (2.36%), and Alternaria (1.39%). The large proportion of “other genera” in the hindgut indicated that many genera had not been identified (Additional File 5). Analysis of significant differences among different groups showed that Cladosporium was more abundant in the ileum (I) than in the stomach, caecum, ventral colon, and dorsal colon (P < 0.01 or P < 0.05) (Fig. 2C, Additional File 6).

To better assess the structural differences between samples, all OTUs were subjected to NMDS analysis (Fig. 3). The samples were formed into two distinct clusters, UG and LG. The LG samples were more dispersed, indicating a greater difference in the fungal community across the segments. In contrast, the UG samples were relatively clustered, indicating higher compositional similarity. LDA and LEfSe were used to further search for biomarkers with statisticant differences between groups [38]. Twenty funfal genera had significant biological differences between different parts of the intestine, and 36 taxa at different taxonomic levels had LDA values greater than 4 (Additional File 7). At the phylum level, the abundance of Basidiomycota was the highest in the stomach, that of Ascomycota was highest in the ileum, and that of Neocallimastigomycota was highest in the dorsal colon compared with those of different parts of the GIT, and the ventral colon did not have a taxon with LDA > 4.

Fig. 3
figure 3

Nonmetric multidimensional scaling (NMDS) with clustering represents the fungal structure dissimilarity found among samples from Mongolian horse GIT compartments. In each figure, S represents the stomach segment, J represents the jejunum segment, I represents the ileum segment, C represents the caecum segment, VC represents the ventral colon segment, and DC represents the dorsal colon segment

Functional prediction of the FunGuild gene

FUNGuild was used to predict the function of fungal communities in different parts of the equine GIT [39], and based on the trophic mode, fungi were divided into symbiotrophs, saprotrophs, and pathotrophs. A total of 1,587 OTUs in this study were divided into nine nutrient types, and 43 functional groups were detected. The nine main functional groups with high abundance were analysed (Fig. 4A). The saprophyte-symbiotic trophic type in the ileum (I) group, pathological saprophyte-symbiotic trophic type in the jejunum (J) and symbiotic trophic type in the dorsal colon (DC) group were significantly higher than in other segments (Fig. 4B). The main functional groups of pathotrophic fungi in the GIT of Mongolian horses were animal and plant pathogens. The relative abundance of animal and plant pathogens in the foregut was significantly higher than that in the hindgut (P < 0.01), and the highest abundance was observed in the ileum (I) group (Fig. 4C).

Fig. 4
figure 4

Functional analysis of gastrointestinal fungi in Mongolian equine. A: Trophic mode histogram of relative abundance; B: Abundance clustering heatmap; C: Histogram of relative abundance based on the guilds. In each figure, S represents the stomach segment, J represents the jejunum segment, I represents the ileum segment, C represents the caecum segment, VC represents the ventral colon segment, and DC represents the dorsal colon segment

Discussion

Compared with traditional culture methods, high-throughput sequencing technology is an efficient, rapid, and accurate approach to analysing microbial community composition, especially for studying microbial species with harsh environmental requirements [40]. At present, studies on equine intestinal microbiota have mostly focused on the bacterial diversity of stool samples [4, 7, 12, 41, 42], but these studies have limitations. Indeed, it should be borne in mind that stool samples can only represent the microbial structure of the distal part of the hindgut (colon and rectum) and cannot truly reflect the composition of the entire gastrointestinal microbial community. In addition, significant emphasis has been placed on the role of fungal diversity in regulating gut microbiota function.

Fungi are an important component of the intestinal flora and affect the balance of intestinal microecology. Simultaneously, they also rely on and harmoniously coexist with gut bacteria to maintain normal intestinal function. In this study, we compared the structure and composition of the fungal community in different parts of the gastrointestinal tract of Mongolian horses by sequencing an internal transcribed spacer region. Currently, studies on the structure and composition of fungi in different parts of the Mongolian horse gut are limited to the analysis of clone libraries [24]. Studies using high-throughput sequencing technology to analyse the whole-gut fungal community and functional prediction in horses have been reported [14], but no systematic studies have been conducted on the whole-gut fungal community diversity in Mongolian horses.

Fungal community composition and hindgut fermentation characteristics in the gastrointestinal tract of Mongolian horses

The results of the present study showed that the dominant fungal phyla in the GIT of Mongolian horses were Basidiomycota, Ascomycota, and Neocallimastigomycota, which accounted for more than 97% of the total intestinal fungi, consistent with the findings of past literature on the microbial life of Mongolian horse feeces [43]. Basidiomycota and Ascomycota are non-strictly anaerobic fungi. Most fungi belonging to the Ascomycota are saprophytes, which decompose plant residues and degrade soil organic matter [44]. However, most fungi belonging to Basidiomycota have a strong ability to decompose lignocellulose [45]. Neocallimastigomycota plays an important role in the degradation of structural carbohydrates. These three fungal groups are abundant in the GIT of Mongolian horses, which is consistent with the results of a previously conducted study on equine faecal fungal groups [14]. This indicates that Mongolian horses are more adapted to consuming plant food resources, and the metabolic effect accounts for the largest proportion during the prediction of microbial function.

The results of this study revealed that Ascomycota is the dominant phylum in the GIT of Mongolian horses, with a relative abundance of 28.33% (Additional File 8). Saccharomycetes is an important member of Ascomycetes. Reportedly, feeding live yeast to ruminants can significantly increase the abundance of rumen microorganisms and improve production performance [46]. Interestingly, supplementing live yeast to horse diets can reportedly significantly improve feed digestibility and performance, increase feed intake, reduce the number of pathogenic intestinal bacteria, and reduce the negative impact of horse production on the environment [47]. In our study, we found that Neocallimastigomycota represented the dominant group of anaerobic fungi in the hindgut, especially in the dorsal colon, and its relative abundance was significantly higher than that in the foregut (P < 0.01). Current evidence suggests that Neocallimastigomycota plays an important role in host cellulose digestion, suggesting that horses are hindgut fermenters. Neocallimastigomycota is mainly found in the digestive tract of animals fed high-fibre vegetative diets. It destroys fibrous diets through highly active cellulolytic enzymes at the tip of the hyphae and mechanical penetration during growth. In turn, it increases the access and possibility of other types of cellulose-degrading microorganisms to obtain carbohydrates. A rich literature substantiates that, although the proportion of anaerobic fungi in intestinal microorganisms is not high, they can release more than 50% of fermentable polysaccharides from plant feed ingested by moving objects [46, 48]. Moreover, in vitro studies have shown that anaerobic fungi can release up to 95% fermentable polysaccharides from plant leaves during the culture period [49]. This study substantiated comparable fungal diversity in different parts of the GIT. The richness index of the hindgut fungal community was significantly higher than that of the foregut. Futhermore, fungi associated with higher abundance were found in the hindgut, which played an important role in cellulose degradation.

At the family level, the relative abundances of Neocallimastigaceae and Tuberaceae in the hindgut were significantly higher than those in the foregut(P < 0.01). Neocallimastigaceae is the only obligate anaerobic group of fungi and the only group currently known to contain six genera. At the genus level, the relative abundance of the six known genera of Neocallimastigaceae in the hindgut is extremely low, indicating that many genera of Neocallimastigaceae are found in the Mongolian horse hindgut, and the role of these genera in the gut is currently unclear, warranting further research. These fungi are not necessarily active in the GIT tract. Like humans, these may be present in their food [50].

Diversity and disease resistance of fungal communities in the gastrointestinal tract

The proportions of Basidiomycota and Ascomycota were the highest in the foregut contents of Mongolian horses. In humans, the relative abundance of Basidiomycota and Ascomycota is positively associated with the occurrence of inflammatory bowel disease (IBD) and colorectal cancer (CRC) [51, 52]. Therefore, the stability and balance of the fungal community structure in the foregut (especially the stomach and ileum) is an important index directly related to the intestinal health of horses.

At the family level, the relative abundances of Tremellales incertae sedis, Davidiellaceae, and Pleosporaceae in the foregut were significantly higher than those in the hindgut (P < 0.01). However, little is known about the specific role of these family-level fungi in equine physiology and metabolism. At the genus level, the relative abundances of Cryptococcus, Cladosporium, Alternaria, and Sarocladium in the foregut were significantly higher than those in the hindgut (P < 0.01 or P < 0.05). These fungi are opportunistic pathogens in the human body. For example, Cladosporium is not only significantly enriched in the gut of subjects with nonalcoholic fatty liver disease (NAFLD) but also correlated with the severity of the disease and can be used as a disease predictor of nonalcoholic steatohepatitis and significant fibrosis. Cladosporium is one of the most common fungal genera in the GIT [53] and has been reported to be a human pathogen in recent years [54, 55]. Meanwhile, studies have found that the abundance of Cladosporium in the intestinal tract is significantly increased in overweight and obese people. Moreover, it can directly or indirectly participate in the occurrence and development of inflammatory bowel disease and HIV infection through interaction with bacteria [56, 57]. Fusarium is also an opportunistic pathogen whose abundance is significantly increased in alcoholic liver disease and has been implicated in the progression of inflammatory bowel disease (IBD) and acquired immunodeficiency syndrome (AIDS) [58]. In short, Mongolian horses can still maintain good health and survive in harsh environments all year long despite the abundance of saprophytic and pathogenic fungi in their foregut, substantiating that Mongolian horses have strong disease resistance.

Fungal community diversity in different gastrointestinal tract segments

We found that Neocallimastigomycota and Basidiomycota were predominantly found in the contents of the hindgut of Mongolian horses, which is inconsistent with the results of studies on intestinal microbial communities of horses in faecal samples [43]. This finding suggests that stool may not fully represent the composition of the hindgut microbiome, or it may be related to the sequencing method (DNA extraction technology and sequencing platform), breed, geographical location, or diet. Indeed, further studies are needed to verify the relevance of these factors.

The composition of the fungal community in the contents of different parts of the Mongolian horses’ GIT was significantly different. The ecological environment in the whole intestinal tract was not unchanged but showed significant regional changes. Moreover, the fungal community in the GIT of Mongolian horses could be divided into two distinct regions, the foregut and the hindgut, consistent with our previous study on the bacterial diversity in the GIT of Mongolian horses. While the different parts of the foregut had similar fungal communities, those of the hindgut fungal community were highly variable among the specific parts, as well as in different horses. The NMDS plot shows that individual horses differed most in the caecum fungal community. This also indicates that there may be an interaction between different microbial or fungal groups, which jointly maintain the intestinal microbial balance and body health of Mongolian horses.

Functional analysis of fungal composition in the gastrointestinal tract

There are significant differences in fungal function between the foregut and hindgut in horses, and they were confirmed by FUNGuild prediction of fungal community function in the GIT of horses. The functional analysis results were consistent with the results of fungal diversity analysis in different parts. We found that fungi in the foregut of Mongolian horses were mainly pathologically nutritive and contained a large number of animal and plant pathogens, especially in the small intestine (jejunum and ileum), indicating that the foregut may be the most important immune site in the digestive system of Mongolian horses, which accounted for the high disease resistance of Mongolian horses. The unassigned functional group in the hindgut was significantly higher than that in the foregut, indicating that the functions of fungal groups in the hindgut have not been fully explored, and further research is warranted. The FUNGuild function prediction and analysis methods still have some limitations. In the future, metagenomic sequencing, related functional gene analysis, metabolome sequencing, and other methods should be combined to more accurately analyse and grasp the function of the Mongolian horse intestinal fungal community. Mongolian horses have strong adaptability and disease resistance, which is related to their excellent breed characteristics and long-term grazing and foraging, which may promote the formation of unique intestinal microbiota.

This study has several limitations, which decrease the robustness of our findings. First, compared with whole-genome shotgun sequencing, ITS amplicon sequencing has limited resolution and can only be annotated to the family or genus level. In addition, given the lack of fungal information in omics databases, less than half of the genes can be annotated after fungal sequencing. Moreover, many bacteria exist in the digestive tract microbiota, which increases the difficulty of obtaining sufficient fungal genomes in the digestive tract microbiota. Finally, the limited sample size in this study may have reduced the efficiency of detecting differences between gastrointestinal regions to some extent. Although this study provides a detailed composition of Mongolian horses’ entire gastrointestinal fungal flora, more studies are needed to determine the influence of other factors, such as age, geographical location, and diet.

Conclusion

Based on ITS amplicon high-throughput sequencing and FUNGuild function prediction, this study provides novel insights into fungal communities and functions in different parts of the digestive tract in healthy Mongolian horses. The following conclusions were obtained: (1) the Mongolian horse is a hindgut fermentative animal, which is related to its unique fungal community structure; (2) there are a large number of pathogenic fungi in the foregut of Mongolian horses, indicating that Mongolian horses have strong disease resistance; and (3) based on the composition and function of fungal groups, the digestive tract of Mongolian horses can be divided into the foregut and hindgut.

Data Availability

The datasets generated and analyzed during the current study are available in the NCBI Sequence Read Archive at https://www.ncbi.nlm.nih.gov/bioproject/PRJNA898862.

Abbreviations

GIT:

Gastrointestinal tracts

UG:

Upper GIT

LG:

Lower GIT

References

  1. Lynch SV, Pedersen O. The human intestinal microbiome in health and disease. N Engl J Med. 2016;375(24):2369–79. https://doi.org/10.1056/NEJMra1600266.

    Article  CAS  PubMed  Google Scholar 

  2. Garrett WS, Gordon JI, Glimcher LH. Homeostasis and inflammation in the intestine. Cell. 2010;140(6):859–70. https://doi.org/10.1016/j.cell.2010.01.023.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Arnold CE, Pilla R, Chaffin MK, Leatherwood JL, Wickersham TA, Callaway TR, et al. The effects of signalment, diet, geographic location, season, and colitis associated with antimicrobial use or Salmonella infection on the fecal microbiome of horses. J Vet Intern Med. 2021;35(5):2437–48. https://doi.org/10.1111/jvim.16206.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Massacci FR, Clark A, Ruet A, Lansade L, Costa M, Mach N. Inter-breed diversity and temporal dynamics of the faecal microbiota in healthy horses. J Anim Breed Genet. 2020;137(1):103–20. https://doi.org/10.1111/jbg.12441.

    Article  CAS  PubMed  Google Scholar 

  5. Garber A, Hastie P, Murray J-A. Factors influencing equine gut microbiota: current knowledge. J Equine Veterinary Sci. 2020;88:102943. https://doi.org/10.1016/j.jevs.2020.102943.

    Article  Google Scholar 

  6. Zhao Y, Li B, Bai D, Huang J, Shiraigo W, Yang L, et al. Comparison of fecal microbiota of Mongolian and Thoroughbred Horses by high-throughput sequencing of the V4 region of the 16S rRNA gene. Asian-Australasian J Anim Sci. 2016;29(9):1345. https://doi.org/10.5713/ajas.15.0587.

    Article  CAS  Google Scholar 

  7. Schoster A, Mosing M, Jalali M, Staempfli H, Weese J. Effects of transport, fasting and anaesthesia on the faecal microbiota of healthy adult horses. Equine Vet J. 2016;48(5):595–602. https://doi.org/10.1111/evj.12479.

    Article  CAS  PubMed  Google Scholar 

  8. Szemplinski KL, Thompson A, Cherry N, Guay K, Smith WB, Brady J, et al. Transporting and exercising unconditioned horses: effects on microflora populations. J Equine Veterinary Sci. 2020;90:102988. https://doi.org/10.1016/j.jevs.2020.102988.

    Article  Google Scholar 

  9. Colombino E, Raspa F, Perotti M, Bergero D, Vervuert I, Valle E, et al. Gut health of horses: effects of high fibre vs high starch diet on histological and morphometrical parameters. BMC Vet Res. 2022;18(1):338. https://doi.org/10.1186/s12917-022-03433-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Salem SE, Maddox TW, Berg A, Antczak P, Ketley JM, Williams NJ, et al. Variation in faecal microbiota in a group of horses managed at pasture over a 12-month period. Sci Rep. 2018;8(1):8510. https://doi.org/10.1038/s41598-018-26930-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Kauter A, Epping L, Semmler T, Antao E-M, Kannapin D, Stoeckle SD, et al. The gut microbiome of horses: current research on equine enteral microbiota and future perspectives. Anim Microbiome. 2019;1(1):1–15. https://doi.org/10.1186/s42523-019-0013-3.

    Article  Google Scholar 

  12. Perricone V, Sandrini S, Irshad N, Comi M, Lecchi C, Savoini G, et al. The role of yeast Saccharomyces cerevisiae in supporting Gut Health in Horses: an updated review on its Effects on Digestibility and Intestinal and Fecal Microbiota. Animals. 2022;12(24):3475. https://doi.org/10.3390/ani12243475.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Chaucheyras-Durand F, Sacy A, Karges K, Apper E. Gastro-intestinal microbiota in equines and its role in health and disease: the black box opens. Microorganisms. 2022;10(12):2517. https://doi.org/10.3390/microorganisms10122517.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Lan Y, Li Y, Yu G, Zhang Z, Irshad I. Dynamic changes of gut fungal community in horse at different health states. Front Veterinary Sci. 2022;9:1047412. https://doi.org/10.3389/fvets.2022.1047412.

    Article  Google Scholar 

  15. Morrison PK, Newbold CJ, Jones E, Worgan HJ, Grove-White DH, Dugdale AH, et al. The equine gastrointestinal microbiome: impacts of weight-loss. BMC Vet Res. 2020;16:1–18. https://doi.org/10.1186/s12917-020-02295-6.

    Article  Google Scholar 

  16. Moreau MM, Eades SC, Reinemeyer CR, Fugaro MN, Onishi JC. Illumina sequencing of the V4 hypervariable region 16S rRNA gene reveals extensive changes in bacterial communities in the cecum following carbohydrate oral infusion and development of early-stage acute laminitis in the horse. Vet Microbiol. 2014;168(2–4):436–41. https://doi.org/10.1016/j.vetmic.2013.11.017.

    Article  CAS  PubMed  Google Scholar 

  17. Paul LJ, Ericsson AC, Andrews FM, Keowen ML, Morales Yniguez F, Garza F Jr, et al. Gastric microbiome in horses with and without equine glandular gastric disease. J Vet Intern Med. 2021;35(5):2458–64. https://doi.org/10.1111/jvim.16241.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Lara F, Castro R, Thomson P. Changes in the gut microbiome and colic in horses: are they causes or consequences? Open Veterinary Journal. 2022;12(2):242–9. https://doi.org/10.5455/OVJ.2022.v12.i2.12.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Uzal FA, Arroyo LG, Navarro MA, Gomez DE, Asín J, Henderson E. Bacterial and viral enterocolitis in horses: a review. J Vet Diagn Invest. 2022;34(3):354–75. https://doi.org/10.1177/10406387211057.

    Article  PubMed  Google Scholar 

  20. Su S, Zhao Y, Liu Z, Liu G, Du M, Wu J, et al. Characterization and comparison of the bacterial microbiota in different gastrointestinal tract compartments of mongolian horses. Microbiologyopen. 2020;9(6):1085–101. https://doi.org/10.1002/mbo3.1020.

    Article  CAS  PubMed  Google Scholar 

  21. Ericsson AC, Johnson PJ, Lopes MA, Perry SC, Lanter HR. A microbiological map of the healthy equine gastrointestinal tract. PLoS ONE. 2016;11(11):e0166523. https://doi.org/10.1371/journal.pone.0166523.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Costa M, Silva G, Ramos R, Staempfli H, Arroyo L, Kim P, et al. Characterization and comparison of the bacterial microbiota in different gastrointestinal tract compartments in horses. Vet J. 2015;205(1):74–80. https://doi.org/10.1016/j.tvjl.2015.03.018.

    Article  CAS  PubMed  Google Scholar 

  23. Costa MC, Weese JS. Understanding the intestinal microbiome in health and disease. Veterinary Clinics: Equine Practice. 2018;34(1):1–12. https://doi.org/10.1016/j.cveq.2017.11.005.

    Article  CAS  PubMed  Google Scholar 

  24. Mura E, Edwards J, Kittelmann S, Kaerger K, Voigt K, Mrázek J, et al. Anaerobic fungal communities differ along the horse digestive tract. Fungal Biology. 2019;123(3):240–6. https://doi.org/10.1016/j.funbio.2018.12.004.

    Article  PubMed  Google Scholar 

  25. Edwards JE, Forster RJ, Callaghan TM, Dollhofer V, Dagar SS, Cheng Y, et al. PCR and omics based techniques to study the diversity, ecology and biology of anaerobic fungi: insights, challenges and opportunities. Front Microbiol. 2017;8:1657. https://doi.org/10.3389/fmicb.2017.01657.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Ariyawansa HA, Hyde KD, Jayasiri SC, Buyck B, Chethana KT, Dai DQ, et al. Fungal diversity notes 111–252—taxonomic and phylogenetic contributions to fungal taxa. Fungal Divers. 2015;75:27–74. https://doi.org/10.1007/s13225-015-0346-5.

    Article  Google Scholar 

  27. Zhang B-W, Li M, Ma L-C, Wei F-W. A widely applicable protocol for DNA isolation from fecal samples. Biochem Genet. 2006;44:494–503. https://doi.org/10.1007/s10528-006-9050-1.

    Article  CAS  Google Scholar 

  28. Schmidt V, Plenz B, Pfaff M, Pees M. Disseminated systemic mycosis in veiled chameleons (Chamaeleo calyptratus) caused by Chamaeleomyces granulomatis. Vet Microbiol. 2012;161(1–2). https://doi.org/10.1016/j.vetmic.2012.07.017. 145 – 52; doi.

  29. Zhao S, Liu D, Ling N, Chen F, Fang W, Shen Q. Bio-organic fertilizer application significantly reduces the Fusarium oxysporum population and alters the composition of fungi communities of watermelon Fusarium wilt rhizosphere soil. Biol Fertil Soils. 2014;50:765–74. https://doi.org/10.1007/s00374-014-0898-7.

    Article  Google Scholar 

  30. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17(1):10–2. https://doi.org/10.14806/ej.17.1.200.

    Article  Google Scholar 

  31. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200. https://doi.org/10.1093/bioinformatics/btr381.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Haas BJ, Gevers D, Earl AM, Feldgarden M, Ward DV, Giannoukos G, et al. Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 2011;21(3):494–504. https://doi.org/10.1101/gr.112730.110. http://www.genome.org/cgi/doi/.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10(10):996–8. https://doi.org/10.1038/nmeth.2604.

    Article  CAS  PubMed  Google Scholar 

  34. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–. https://doi.org/10.1016/S0022-2836(05)80360-2. 10; doi.

    Article  CAS  PubMed  Google Scholar 

  35. Kõljalg U, Nilsson RH, Abarenkov K, Tedersoo L, Taylor AF, Bahram M, et al. Towards a unified paradigm for sequence-based identification of fungi. In.: Wiley Online Library; 2013.

  36. White JR, Nagarajan N, Pop M. Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput Biol. 2009;5(4):e1000352. https://doi.org/10.1371/journal.pcbi.1000352.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Lundberg DS, Yourstone S, Mieczkowski P, Jones CD, Dangl JL. Practical innovations for high-throughput amplicon sequencing. Nat Methods. 2013;10(10):999–1002. https://doi.org/10.1038/nmeth.2634.

    Article  CAS  PubMed  Google Scholar 

  38. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12:1–18. https://doi.org/10.1186/gb-2011-12-6-r60.

    Article  Google Scholar 

  39. Nguyen NH, Song Z, Bates ST, Branco S, Tedersoo L, Menke J, et al. FUNGuild: an open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 2016;20:241–8. https://doi.org/10.1016/j.funeco.2015.06.006.

    Article  Google Scholar 

  40. Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31(9):814–21. https://doi.org/10.1038/nbt.2676.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Shepherd ML, Swecker WS Jr, Jensen RV, Ponder MA. Characterization of the fecal bacteria communities of forage-fed horses by pyrosequencing of 16S rRNA V4 gene amplicons. FEMS Microbiol Lett. 2012;326(1):62–. https://doi.org/10.1111/j.1574-6968.2011.02434.x. 8; doi.

    Article  CAS  PubMed  Google Scholar 

  42. Weese JS, Holcombe S, Embertson R, Kurtz K, Roessner H, Jalali M, et al. Changes in the faecal microbiota of mares precede the development of post partum colic. Equine Vet J. 2015;47(6):641–9. https://doi.org/10.1111/evj.12361.

    Article  CAS  PubMed  Google Scholar 

  43. XING ZC. Mongolian horse and thoroughbred intesti-nal fungal diversity preliminary research and analysis. In., vol. Master: Inner Mongolia Agricultural University; 2016.

    Google Scholar 

  44. Frey SD, Knorr M, Parrent JL, Simpson RT. Chronic nitrogen enrichment affects the structure and function of the soil microbial community in temperate hardwood and pine forests. For Ecol Manag. 2004;196(1):159–71. https://doi.org/10.1016/j.foreco.2004.03.018.

    Article  Google Scholar 

  45. Liggenstoffer AS, Youssef NH, Couger M, Elshahed MS. Phylogenetic diversity and community structure of anaerobic gut fungi (phylum neocallimastigomycota) in ruminant and non-ruminant herbivores. ISME J. 2010;4(10):1225–35. https://doi.org/10.1038/ismej.2010.49.

    Article  PubMed  Google Scholar 

  46. Solomon KV, Haitjema CH, Henske JK, Gilmore SP, Borges-Rivera D, Lipzen A, et al. Early-branching gut fungi possess a large, comprehensive array of biomass-degrading enzymes. Science. 2016;351(6278):1192–5. https://doi.org/10.1126/science.aad1431.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Elghandour M, Tan Z, Abu Hafsa S, Adegbeye M, Greiner R, Ugbogu E, et al. Saccharomyces cerevisiae as a probiotic feed additive to non and pseudo-ruminant feeding: a review. J Appl Microbiol. 2020;128(3):658–74. https://doi.org/10.1111/jam.14416.

    Article  CAS  PubMed  Google Scholar 

  48. Sijtsma L, Tan B. Degradation of perennial ryegrass leaf and stem cell walls by the anaerobic fungus Neocallimastix sp. strain CS3b. Appl Environ Microbiol. 1996;62(4):1437–40. https://doi.org/10.1128/aem.62.4.1437-1440.1996.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Wallace RJ. Rumen microbiology, biotechnology and ruminant nutrition: the application of research findings to a complex microbial ecosystem. FEMS Microbiol Lett. 1992;100(1–3). https://doi.org/10.1111/j.1574-6968.1992.tb14088.x. 529 – 34; doi.

  50. Auchtung TA, Fofanova TY, Stewart CJ, Nash AK, Wong MC, Gesell JR, et al. Investigating colonization of the healthy adult gastrointestinal tract by fungi. MSphere. 2018;3(2):e00092–18. https://doi.org/10.1128/mSphere.00092-18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Sokol H, Leducq V, Aschard H, Pham H-P, Jegou S, Landman C, et al. Fungal microbiota dysbiosis in IBD. Gut. 2017;66(6):1039–48. https://doi.org/10.1136/gutjnl-2015-310746.

    Article  CAS  PubMed  Google Scholar 

  52. Coker OO, Nakatsu G, Dai RZ, Wu WKK, Wong SH, Ng SC, et al. Enteric fungal microbiota dysbiosis and ecological alterations in colorectal cancer. Gut. 2019;68(4):654–62. https://doi.org/10.1136/gutjnl-2018-317178.

    Article  CAS  PubMed  Google Scholar 

  53. Hallen-Adams HE, Suhr MJ. Fungi in the healthy human gastrointestinal tract. Virulence. 2017;8(3):352–8. https://doi.org/10.1080/21505594.2016.1247140.

    Article  CAS  PubMed  Google Scholar 

  54. Ökmen B, Collemare J, Griffiths S, van der Burgt A, Cox R, de Wit PJ. Functional analysis of the conserved transcriptional regulator CfWor1 in C ladosporium fulvum reveals diverse roles in the virulence of plant pathogenic fungi. Mol Microbiol. 2014;92(1):10–27. https://doi.org/10.1111/mmi.12535.

    Article  CAS  PubMed  Google Scholar 

  55. Mercier E, Peters I, Billen F, Battaille G, Clercx C, Day M, et al. Potential role of Alternaria and Cladosporium species in canine lymphoplasmacytic rhinitis. J Small Anim Pract. 2013;54(4):179–83. https://doi.org/10.1111/jsap.12049.

    Article  CAS  PubMed  Google Scholar 

  56. Qiu X, Zhang F, Yang X, Wu N, Jiang W, Li X, et al. Changes in the composition of intestinal fungi and their role in mice with dextran sulfate sodium-induced colitis. Sci Rep. 2015;5(1):1–12. https://doi.org/10.1038/srep10416.

    Article  Google Scholar 

  57. Mukherjee PK, Chandra J, Retuerto M, Sikaroodi M, Brown RE, Jurevic R, et al. Oral mycobiome analysis of HIV-infected patients: identification of Pichia as an antagonist of opportunistic fungi. PLoS Pathog. 2014;10(3):e1003996. https://doi.org/10.1371/journal.ppat.1003996.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Ott SJ, Kühbacher T, Musfeldt M, Rosenstiel P, Hellmig S, Rehman A, et al. Fungi and inflammatory bowel diseases: alterations of composition and diversity. Scand J Gastroenterol. 2008;43(7):831–41. https://doi.org/10.1080/00365520801935434.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank the Academy of Agricultural and Animal Husbandry Sciences of Inner Mongolia Autonomous Region for their support in sample handling and preservation.

Funding

National Natural Science Foundation of China (31860636, 31960657, 31961143025); Natural Science Foundation of Inner Mongolia Autonomous Region of China (2020BS03035, 2021LHMS03012); China Postdoctoral Foundation “Postdoctoral Talents in Western China” (2020M673547XB); Innovation Fund of Inner Mongolia Academy of Agricultural and Animal Husbandry Science (2022CXJJM04).

Author information

Authors and Affiliations

Authors

Contributions

SS, XR: conceptualization, data analysis, and manuscript drafting. HW, CC, and LW: animal sampling and determination. HH, XZ, LY, and MD: data analyst. YH, JT, YL, and YL: supervision, review, and editing. MD, YZ: conception, design, and financial support. All authors contributed to the article and approved the final manuscript.

Corresponding authors

Correspondence to Shaofeng Su or Manglai Dugarjaviin.

Ethics declarations

Ethics approval and consent to participate

Animal experiments were approved by the Inner Mongolia Agriculture University Animal Welfare Ethics Committee, and all procedures were conducted in accordance with the guidelines of the Chinese Animal Protection Association. The study was carried out in compliance with the ARRIVE guidelines.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Yiping Zhao and Xiujuan Ren have contributed equally to this work and share the first authorship.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhao, Y., Ren, X., Wu, H. et al. Diversity and functional prediction of fungal communities in different segments of mongolian horse gastrointestinal tracts. BMC Microbiol 23, 253 (2023). https://doi.org/10.1186/s12866-023-03001-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12866-023-03001-w

Keywords