- Open Access
Comparative analysis of gut microbiota in healthy and diarrheic yaks
Microbial Cell Factories volume 21, Article number: 111 (2022)
Yak (Bos grunniens) mainly inhabiting Tibet Plateau, displayed a high incidence of diarrhea due to harsh living environment and nutritional deficit. Gut microbial community has been reported to be closely related to many diseases including diabetes, obesity and inflammatory bowel disease, but information regarding diarrheic influence on gut microbiota in yaks remains scarce. Here, this study was performed to investigate the gut bacterial and fungal alternations of diarrheic yaks.
Results revealed that the gut bacterial and fungal communities of diarrheic yaks showed a distinct decline in alpha diversity, accompanied by significant shifts in taxonomic compositions. Specifically, diarrhea caused a distinct increase in the relative abundance of 1 phylum and 8 genera as well as a distinct decrease in 3 phyla and 30 genera. Fungal taxonomic analysis indicated that the relative richness of 1 phylum and 2 genera dramatically increased, whereas the relative richness of 2 phylum and 43 genera significantly decreased during diarrhea. Surprisingly, 2 bacterial genera and 5 fungal genera even cannot be detected in the gut microbiota of diarrheic yaks.
In summary, this study indicated that the gut bacterial and fungal compositions and diversities of yaks altered significantly during diarrhea. Moreover, these findings also contribute to understanding the gut microbial composition and diversity of yaks and developing strategies to alleviate and prevent diarrhea from gut microbial perspective.
Ruminant intestines harbor trillions of microbes, that serve key roles in metabolism, digestive absorption, intestinal homeostasis and host health [1, 2]. Additionally, gut microbiota has also been demonstrated to function in immune system maturation, permeability and epithelial differentiation of intestine [3, 4]. Gut-residing beneficial bacteria can limit the colonization of foreign pathogens in the intestine via secreting antimicrobial peptide, regulating intestinal environment and competing nutrition, which was deemed as a natural barrier against pathogenic bacteria invasion [5, 6]. Statistically, the normal intestine contains more than 1014 microbes, approximately 10 times the total quantity of host cells . Among them, intestinal bacteria account for approximately 98% of the total microbial community, whereas the rest contains fungi (0.1%), viruses and protozoa [7,8,9]. The stabilized gut microbiota is required for various complex physiological and metabolic processes, whereas gut microbial dysbiosis may result in multiple gastrointestinal diseases including diarrhea, enteritis and irritable bowel syndrome [10, 11]. Although these microbes colonize the intestine, they may result in systemic effects . Increasing evidence suggests that gut microbial dysbiosis has an impact on intestinal functions with negative influence beyond the gastrointestinal system, impairing the functioning of other organs including liver and brain [13, 14].
Diarrhea is a common disease that is often accompanied by gastrointestinal dysfunction and may cause decreased productivity, weakened immunity and even mortality [15, 16]. Numerous studies demonstrated that diarrhea occurred in almost all species, especially in newborn piglets, goats and cattle with immature gut microbiota, which was regarded as one of the important constraint to livestock sector [17, 18]. Certain gut-residing bacteria and fungi change between dominant and weak populations accompanied by diarrhea, indicating some inevitable relationships may exist between gut microbial dysbiosis and diarrhea, but specific links and laws are still unclear [19, 20].
Currently, high-throughput sequencing technology has been successfully applied to investigate gut microbial alterations after the onset of various diseases, making it possible to deeply characterize the potential relationship between gut microbiota and some diseases [6, 21]. Furthermore, in-depth analysis of the complicated gut microbiota is beneficial to further understand the mechanisms contributing to ill health and timely formulate strategies to minimize the collateral damage . Yak is an aboriginal breed mainly inhabits in the Tibet Plateau (average elevation above 4000 m), characterized by cold resistance and strong adaptability [23, 24]. This indigenous breed has resided in the high-altitude hypoxic environment for thousands of years and has evolved unique digestive characteristics and gut microbiota that contribute to adapting to high-fiber diet, but also make them susceptible to various gastrointestinal diseases, especially diarrhea [25, 26]. Consequently, the gut microbiota of yaks plays a more noticeable role in various physiological functions compared with poultry and other mammal. However, the potential relationship between the gut microbiota in yaks and diarrhea remain to be determined. In this study, we investigated the gut bacterial and fungal compositions and variabilities of healthy and diarrheic yaks.
Materials and methods
In the present study, a total of 12 12-month-old free-range yaks with similar blood profile and weight (6 healthy and 6 diarrheic) were selected for sample collection from Qinghai Province, China, including 3 male and female yaks in each group. The diarrheic yaks were diagnosed by professional veterinarian and didn’t received any treatment before sample acquisition. On the day of sample acquisition, each yak was placed in a separate enclosure to prevent potential contamination between samples. The rectal feces (approximately 200 g) were collected from each selected yak using sampler. The obtained fresh fecal samples were sub-sampled from the central proportion to decrease contamination via flooring and bedding. The final samples were snap-frozen using liquid nitrogen and stored at − 80 ℃ for further analysis.
DNA extraction and high-throughput sequencing
Total bacterial and fungal genomic DNA were extracted from 12 frozen feces (approximately 200 mg) of healthy and diarrheic yaks utilizing QIAamp DNA Mini Kit (QIAGEN, Hilden, Germany) based on manufacturer’s protocol recommendations. To further ensure the quality of extracted total genomic DNA, 0.8% agarose gel electrophoresis and UV–Vis spectrophotometer (NanoDrop 2000, United States) were used for evaluating the integrity and concentration of the extract, respectively. The universal primers including bacterial 16S rDNA (338F: ACTCCTACGGGAGGCAGCA and 806R: GGACTACHVGGGTWTCTAAT) and fungal ITS (ITS5F: GGAAG TAAAAGTCGTAACAAGG and ITS2R: GCTGCGTTCTTCATCGA TGC) gene primers were used for amplifying V3/V4 hypervariable and ITS2 regions, respectively. Subsequently, the amplified products were evaluated by 2.0% agarose gel electrophoresis and then subjected for purification and recycle. The recycled product were conducted fluorescent quantitation and the samples were mixed on a pro-rata basis following the fluorescence quantitative result and sequencing amount requirement. Based on the manufacturer’s protocol, the purified products were applied to construct sequencing library utilizing Illumina TruSeq (Illumina, United States). Prior to the sequencing, the prepared libraries required to be further processed including purification, quality control and fluorescence quantification. The libraries that passed the quality inspection and showed only one peak were considered qualified. Finally, the qualified libraries were diluted and denatured to single-stranded and then performed 2 × 300 bp paired-end sequencing.
Bioinformatics and statistical analysis
The raw data need to be preprocessed to obtain reliable results for subsequent analysis. The raw reads were performed quality filtering using Trimmomatic software (v0.33). Afterwards, the Cutadapt software (1.9.1) was employed to recognize and remove primer sequences to obtain clean reads. The clean reads were spliced using Usearch software (v10) and the sequences less than 200 bp were abandoned. Additionally, the Uchime software (v4.2) was used for identifying and eliminating chimera to obtain the final effective reads. The effective reads were obtained through clustering and OTU discrimination based on 97% nucleotide-sequence similarity. The Ribosome Database Program classifier and MUSCLE software were used for recognizing the representative sequence and performing phylogenetic analysis of each OTU, respectively. Prior to conducting the bacterial and fungal diversities analysis, the rank abundance and rarefaction curve were structured to dissect sequencing depth. The gut bacterial and fungal diversities indices including Good's coverage, ACE, Simpson and Shannon were generated based on the relative abundance of OTU in each sample. Additionally, beta diversity analysis was also conducted to dissect the differences of the main components in different samples. LEfSe was used for identifying the taxa related to diarrhea in gut bacterial and fungal communities. Statistical analysis was conducted to assess the differences of both groups using GraphPad Prism (version 8.0c) and SPSS statistical program (v19.0). Probability values (means ± SD) < 0.05 were considered statistically significant.
To investigate the gut microbial alterations of yaks during diarrhea, 12 (6 control samples and 6 diarrheic samples) fecal samples were subjected to amplicon sequencing. We totally acquired 907,882 (CY = 451,996, DY = 455,886) and 959,714 (CY = 479,401, DY = 480,313) raw sequences from the V3/4 and ITS2 regions of collected samples, respectively (Table 1). After quality evaluation, a total of 1813,360 (CY = 863,743, DY = 949,617) high-quality sequences were achieved, with a median read count of 71,978 (varying from 49,427 to 77,974) and 79,134 (ranging from 78,367 to 79,540) reads from bacterial and fungal communities per sample, respectively (Table 2). The rarefaction curves and shannon curves of each sample was invariably extended to the right end of the x-axis and showed a saturated tendency, indicating that the quantity and depth of sequencing met the demands of further analysis (Fig. 1A–C, G–I). Following taxonomic assignment, the sequences generated from the V3/4 and ITS2 regions were clustered into 946 bacterial OTUs and 716 fungal OTUs on the basis of 97% sequence similarity (Fig. 1E, F, K, L). Furthermore, 880 OTUs were shared in the bacterial community as well as 458 core OTUs were recognized in the fungal community, accounting for approximately 93.02% and 63.97% of the overall OTUs, respectively (Fig. 1D, J).
Bacterial and fungal diversities in fecal microbiota associated with diarrhea
To further characterize the shifts of gut microbial community in yaks during diarrhea, we calculated the multiple alpha and beta diversity indices of the gut microbiota. Results demonstrated that the Good's coverage estimates in the bacterial and fungal populations ranged from 99.79% to 99.93% and 99.93% to 99.95%, respectively, implying that the most of bacterial and fungal phenotypes in the samples could be detected (Fig. 2A, E). Moreover, statistical analysis revealed that the bacterial ACE index (854.09 17.22 versus 799.86 10.72, p = 0.032) differed significantly, but the Simpson (0.96 ± 0.0041 versus 0.95 ± 0.0059, p = 0.209) and Shannon (6.9466 ± 0.12 versus 6.59 ± 0.17, p = 0.13) indices were not significantly different between CY and DY groups (Fig. 2B–D). These results demonstrated that the gut bacterial diversity yaks was strongly influenced by the diarrhea. Additionally, we also observed a obvious decrease in the gut fungal diversity index during diarrhea including ACE, Simpson and Shannon, suggesting that diarrhea observably reduce the gut fungal diversity and abundance of yaks (Fig. 2F–H). The principal coordinate analysis (PCoA) was used for evaluating the similarity and variability between intergroup and intra-group samples. PCoA scatterplot of gut bacterial and fungal communities revealed a separation of individuals in the CY and DY groups, which was in line with the UPGMA analytical results, suggesting that the gut microbial principal compositions of yaks were strongly influenced by the diarrhea (Fig. 2I–N).
Significant alterations of bacterial taxonomic composition in diarrheic yaks
The relative abundances of preponderant taxa at the phylum and genus level were determined using microbial taxon assignment and distinct shifts were observed in the taxonomic composition between both groups during diarrhea. At the phylum level, 16 phyla were recognized in the gut bacterial community, ranging from 13 to 16 phyla per sample. The phyla Firmicutes (72.95%), Bacteroidetes (14.80%) and Verrucomicrobia (4.58%) were the three most dominant phyla in CY group, which together consisted of total 92.33% of the bacterial composition (Fig. 3A). Moreover, Firmicutes (68.95%) was the most preponderant bacterial phylum in the DY groups, followed by the Bacteroidetes (13.67%) and Verrucomicrobia (14.09%), accounting for approximately 96.71% of all bacterial taxa. Other phyla such as Cyanobacteria (CY = 0.39%, DY = 0.56%), Tenericutes (CY = 0.24%, DY = 0.26%), Spirochaetes (CY = 0.18%, DY = 0.035%) and Kiritimatiellaeota (CY = 0.089%, DY = 0.10%) in CY and DY groups were detected in lower abundances. To further evaluate the shifts of gut bacterial compositions during diarrhea, a total of 191 genera were recognized. Among them, Ruminococcaceae_UCG-005 (16.11%) was the most dominant genus in the CY group, followed by Clostridium_sensu_stricto_1 (14.08%) and Christensenellaceae_R-7_group (5.46%) (Fig. 3B). Meanwhile, Akkermansia (14.09%), Solibacillus (10.92%) and Ruminococcaceae_UCG-005 (10.23%) were abundantly present in the DY group, which accounted for approximately 35.24% of the total taxonomic groups identified. The clustering heatmap also revealed the distribution of the identified bacterial genus as well as a variability of the gut bacterial community in diarrheic yaks (Additional file 1: Fig. S1).
To further characterize the changes of taxonomic composition in yaks during diarrhea, Metastats analysis was conducted for different classification levels (Fig. 4 and Additional file 1:Fig. S3). At the phylum level, Proteobacteria, Chloroflexi and Fibrobacteres were dramatically more preponderant in the CY group than in the DY group, whereas the Verrucomicrobia was lower (P < 0.05 or P < 0.01). Moreover, a comparison of the DY and CY groups indicated a distinct decline in the richness of 30 bacterial genus as well as a significant increase in the richness of 8 bacterial genus. Given that this discriminant analysis may not detect the whole taxa, LEfSe coupled with LDA scores was applied to recognize the specific bacteria associated with diarrhea (Fig. 6A, B). Besides the above-mentioned differential taxa, we also found that several bacterial genera including uncultured_bacterium_f_Planococcaceae was the most preponderant in the DY group, whereas Escherichia_Shigella was dramatically overrepresented in the CY group.
Significant shifts of gut fungal compositions in yaks during diarrhea
In this study, a total of 8 phyla and 248 genera were identified in gut fungal community and the main preponderant phyla and genera were shown in Fig. 3C, D. Specifically, the phyla Ascomycota (CY = 69.20%, DY = 81.14), Basidiomycota (CY = 16.57%, DY = 15.15%) and Neocallimastigomycota (CY = 5.34%, DY = 0.60%) were abundantly present in both groups, regardless of health status, accounting for over 90% of total fungi taxa. Other fungal phyla including Rozellomycota (CY = 0.050%, DY = 0.065%), Mortierellomycota (CY = 0.075%, DY = 0.035%), Mucoromycota (CY = 0.083%, DY = 0.014%) and Olpidiomycota (CY = 0.014%, DY = 0.0011%) in CY and DY groups were identified in a low richness. Among identified genera, Thelebolus (13.04%) was the most dominant fungal genera in the CY group, followed by Preussia (10.75%) and Podospora (6.26%). However, the preponderant fungal genera recognized in the DY group were Thelebolus (30.62%), Preussia (12.74%) and Ascobolus (7.43%), which were different from the CY group. The distribution and relative richness of identified fungal genera were further investigated by clustering analysis indicated by the heatmap (Additional file 1: Fig. S2).
Using Metastats analysis to compare the differences in the gut fungal community of both groups (Fig. 5 and Additional file 1: Fig. S4). At the phylum level, Chytridiomycota and Neocallimastigomycota in DY group were significantly reduced, whereas Ascomycota was significantly increased as compared to CY group (p < 0.05 or p < 0.01). Additionally, 45 fungal genera were detected to be significantly different between CY and DY groups. Of these differential taxa, the relative richness of 43 fungal genera distinctly decreased, while 2 fungal genera obviously increased in diarrheic yaks. Among decreased fungus, 4 genera including Humicola, Mucor, Ramularia and Zoellneria even could not be detected in the gut fungal community. LEfSe combined with LDA scores were used for further dissecting the alternations of gut fungal composition (Fig. 6C, D).
Correlation network analysis
Network analysis was performed to dissect the correlations between different bacteria or fungi of gut microbial community (Fig. 7). Results demonstrated that Dorea was positively associated with Coprococcus_3 (0.8531) and Erysipelatoclostridium (0.8951). Caecomyces was positively correlated with Neoascochyta (0.8112), Piromyces (0.8042), Paraphaeosphaeria (0.8671), Itersonilia (0.9231), Erythrobasidium (0.7972), Paraphaeosphaeria (0.8671), Phaeosphaeria (0.8112) and Trichosporon (0.8601). Ruminococcaceae_UCG-004 was positively associated with Lachnospiraceae_UCG-010 (0.8671) and Erysipelatoclostridium (0.8881). Ruminococcaceae_UCG-005 was positively associated with Tyzzerella_4 (0.8951), Coprococcus_3 (0.8881), Erysipelatoclostridium (0.8881), Dorea (0.9301) and Ruminococcaceae_UCG-004 (0.8322). Parasola was positively correlated with Itersonilia (0.8042) and Plectosphaerella (0.8112). Neoascochyta was positively related to Piromyces (0.7273), Itersonilia (0.8741), Truncatella (0.7902) and Buwchfawromyces (0.5873). Truncatella was positively related to Erythrobasidium (0.8322), Neocallimastix (0.8049) and Trichosporon (0.8462). Periconia was positively related to Didymella (0.6434), Caecomyces (0.8741), Paraphaeosphaeria (0.9441), Neosetophoma (0.8371), Trichosporon (0.8182), Ustilago (0.7273), Truncatella (0.6783), Erythrobasidium (0.8462) and Itersonilia (0.9441). Itersonilia was positively associated with Neosetophoma (0.8336), Neocallimastix (0.7904), Plectosphaerella (0.8182), Erythrobasidium (0.8252) and Ustilago (0.8881). Piromyces was positively correlated with Didymella (0.8601), Cortinarius (0.8811), Itersonilia (0.8182), Buwchfawromyces (0.8194), Trichosporon (0.9021), Neocallimastix (0.8121), Kondoa (0.8112) and Cercospora (0.7902). Caecomyces was positively related to Phaeosphaeria (0.8112), Neoascochyta (0.8112), Itersonilia (0.9231), Erythrobasidium (0.7972), Ustilago (0.9091), Pilidium (0.8392), Trichosporon (0.8601), Neocallimastix (0.8774) and Paraphaeosphaeria (0.8671). Pilidium was positively associated with Phaeosphaeria (0.8951), Anaeromyces (0.8531) and Kondoa (0.8601). Paraphaeosphaeria was positively related to Cortinarius (0.8671), Itersonilia (0.9021), Erythrobasidium (0.9371), Truncatella (0.8601) and Trichosporon (0.8322).
Diarrhea is one of the most prevalent diseases of farm animals regardless of species and age, which was regarded as a vital factor causing the decreased global animal productivity . The yaks mainly inhabiting the Tibet plateau possess a high rate of diarrhea, causing substantial economic loss to the yak industry . However, multiple factors including hostile environment, nutritional imbalance and stress response cause the control of diarrhea in yaks particularly difficult . Gut microbial significance has been extensively acknowledged due to its positive roles in immunity, metabolism and intestinal barrier [3, 4]. Moreover, recent studies on gut microbiota have revealed its important role in the control of diarrhea [5, 29]. Therefore, the high diarrhea rate of the yaks may not only be related to their habitat environment but also involve gut microbiota. Systematically understanding the gut microbial information is beneficial to further investigate the disease etiology and develop the potential treatment and prevention options to minimize collateral damage . Currently, research into the gut microbial community of diarrheic ruminant has covered many aspects including goat, cattle and giraffe, but knowledge regarding the gut microbiota in diarrheic yaks remains scarce [16, 18]. In this study, we deeply analyzed and compared the differences of gut microbiota between healthy and diarrheic yaks and revealed a high variability of gut microbial composition and diversity of diarrheic yaks.
Anecdotal evidence indicated that the gut microbiota vary dynamically within limits and affected by species, age and feed, but these physiological fluctuations cannot impair the normal intestinal functions [31, 32]. However, the ecological balance of of gut microbiota can be disrupted by multiple factors including diabetes, antibiotics and diarrhea [33, 34]. Wang et al. reported that the gut bacterial diversity of diarrheic Boer goats decreased dramatically as compared with healthy populations . Moreover, Li et al. also indicated a decreased gut bacterial diversity of giraffes during diarrhea . Consistent with previous investigations, the present study revealed a significantly decreased gut bacterial diversity of diarrheic yaks, suggesting the gut bacterial dysbiosis. The intestine is the leading absorption site, which depends on the normal gut microbial composition and diversity . Early investigations have indicated that the normal gut microbial composition and diversity were the prerequisite for performing its complex physiological functions, while gut microbial dysbiosis may be the central or driving factor of multiple diseases [37, 38]. Previous research demonstrated that the gut microbial dysbiosis might be one of the reasons for the high mortality of diarrheic goats . Additionally, diarrheic yaks are characterized by high mortality and weight loss. Several previous studies indicated that gut microbial dysbiosis can affect host immunity and intestinal permeability, thereby increasing the susceptibility to pathogens [36, 39]. Moreover, some conditioned pathogen may also display pathogenicity under circumstances of gut microbial dysbiosis . Therefore, we speculated that the gut microbial dysbiosis may be one of the reasons for the high mortality and weight loss of diarrheic yaks. As a part of gut microbiota, the intestinal fungus is considered as a vital participant in intestinal functions and host health . Similarly, we also observed a obvious reduction in the gut fungal diversity in diarrheic yaks, indicating the gut fungal dysbiosis. Furthermore, beta diversity analysis showed that despite sharing the same habitat and diet, the major components of healthy and diarrheic yaks' gut bacterial and fungal communities were substantially different, demonstrating that the diarrhea may be a fundamental driving factor for shifts in gut bacterial and fungal communities.
This study indicated that Firmicutes, Bacteroidetes, Ascomycota and Basidiomycota were the most preponderant microbial phyla in yaks, regardless of health status, which was in line with previous investigations on yaks . Additionally, these dominant phyla were also demonstrated to be extensively existed in the goats, giraffes and cattle, indicating their importance in intestinal ecology and functions in ruminants [31, 41]. Some specific bacterial and fungal alternations may reflect the potential relationship between diarrhea and gut microbiota, thus we further characterized the intestinal bacteria and fungi associated with diarrhea. Results indicated a significant increase in 1 bacterial phylum (Verrucomicrobia) and 1 fungal phylum (Ascomycota) as well as a decrease in 3 bacterial phyla (Proteobacteria, Chloroflexi and Fibrobacteres) and 2 fungal phyla (Chytridiomycota and Neocallimastigomycota) in diarrheic yaks. The phylum Proteobacteria has been demonstrated to be abundantly present in the gastrointestinal tract of yaks and contribute to meeting their high nutrient and energy requirements, due to the highly diverse metabolic functions [24, 42]. Moreover, anecdotal evidence indicated that the most members of phylum Fibrobacteres can degrade cellulose . Consistent with the present observations, Wang et al. also revealed that the relative abundance of Verrucomicrobia in diarrheic goats was dramatically increased . Most members of phylum Chytridiomycota can decompose cellulose and chitin. Similarly, the members of Neocallimastigomycota mainly inhabit the digestive tract of mammals, which have the ability to decompose cellulose [44, 45]. Diarrhea is a common gastrointestinal diseases in yaks, which is inevitably accompanied by weight loss. As a strict herbivorous ruminant, the yaks need to consume a large amount of forage to maintain their energy consumption and growth on the high-altitude hypoxic environment of Tibetan plateau . However, we observed that some bacteria associated with cellulose degradation decreased significantly during diarrhea, indicating a decreased ability to digest and degrade food. We presumed that this may be one of the important reasons for the weight loss of yaks during diarrhea.
Importantly, we also found considerable variability in some bacterial and fungal taxonomic taxa of diarrheic yaks, which may play vital roles in gut microbial balance and intestinal functions. Interestingly, most of the altered bacterial phyla and genera in diarrhea yaks showed a downward trend and even 2 bacterial genera could not be detected, implying that these bacteria cannot adapt to the current intestinal environment. We speculated that the intestinal environment of the diarrheic yaks may undergo significant alternations, which in turn selectively inhibited the colonization of some bacteria. Additionally, some decreased bacterial genera (Prevotellaceae_UCG-004, Alloprevotella, Ruminococcaceae_UCG-011, Ruminococcaceae_UCG-005, Ruminococcaceae_UCG-004, Lachnospiraceae_UCG-010, Ruminiclostridium_5, Ruminococcus_2, Butyricicoccus, Coprococcus_3, Butyrivibrio_2, Fibrobacter and Lachnospiraceae_XPB1014_group) were considered as intestinal beneficial bacteria, which play crucial roles in improving the digestion, metabolism, immunity and gut microbiota. Ruminococcaceae, a potential beneficial bacteria commonly found in colon and caecum, has been demonstrated to be involved in degrading cellulose and actively regulating intestinal environment and immunity . Moreover, increased intestinal permeability, non-alcoholic fatty liver, and liver cirrhosis have all been linked to decreased Ruminococcaceae [48, 49]. Previous studies indicated that Ruminiclostridium residing in the intestine was conducive to increase the growth performance and decrease the gastrointestinal diseases . Ruminococcus that mainly inhabit the rumen and hindgut showed the characteristics of degrading cellulose and starch . The higher abundances of Butyricicoccus and Lachnospiraceae in the gut bacterial community is beneficial to alleviate intestinal inflammation [52, 53]. Prevotellaceae displayed the ability to degrade hemicellulose and high carbohydrate . Alloprevotella was previously reported to reduce the risk of cardiovascular disease [55, 56]. As a butyrate-producing bacteria, Butyrivibrio not only decompose polysaccharides, cellulose and starch but also reduce the obesity-induced diabetes and cardiovascular dysfunctions through brain-gut axis [57, 58]. Fibrobacter can degrade cellulose . Notably, the above-mentioned bacterial genera such as Ruminococcaceae_UCG-011, Coprococcus_3, Ruminococcus_2, Ruminococcaceae_UCG-005, Ruminococcaceae_UCG-004, Butyrivibrio_2 and Ruminiclostridium_5 were potential producers of short-chain fatty acids (SCFAs) in the intestine. SCFAs not only improve host immunity but also inhibit the colonization of pathogenic bacteria by regulating the pH of the intestines . Moreover, SCFAs also play important roles in maintaining intestinal function and gut microbial balance . Recent studies on the SCFAs also revealed their important role in anti-inflammatory, anti-cancer aspects and regulating energy intake. Eubacterium ventriosum can cause bacteremia and endocarditis.
Similar to the intestinal bacteria, intestinal fungi are also an important part of the gut microbiota, which play vital roles in intestinal ecosystem and host health . Intestinal fungi have been demonstrated to induce intestinal inflammation and increase intestinal mucosal permeability . Moreover, Li et al. indicated that gut fungal community was the promoter and participant of diarrhea in giraffe . However, the gut fungal role and importance in yaks were chronically neglected due to their lower abundance. In this study, we dissected the shifts of the gut fungal community in yaks during diarrhea, which contribute to providing an insight into gut fungal community in yaks. Similar to the altered gut bacterial community, we observed that the taxonomic composition of gut fungal community changed significantly during diarrhea, characterized by an increase in the abundance of 3 fungal genera and a decline in 46 fungal genera. Among decreased fungi, some of them (Neocallimastix, Periconia, Tilletiopsis and Mortierella) were closely related to host digestion and health. The Neocallimastix can degrade polysaccharide and cellulose . The metabolite of Periconia have antibacterial activity . Tilletiopsis can produce antifungal compounds . Mortierella can produce arachidonic acid . Arachidonic acid plays important roles in the prevention of cardiovascular diseases, diabetes and tumors . Moreover, some fungal genera (Humicola, Mucor, Ramularia and Zoellneria) in the gut fungal community of diarrheic yaks even cannot be detected, suggesting that their growth was significantly restricted.
Gut microbial dysbiosis has been deemed as the pathological mediators of many diseases. Under normal physiological conditions, these microbes inhabiting the gut can engage in a commensal, synergetic or antagonistic relationship, maintaining intestinal homeostasis . Therefore, some gut bacterial and fungal alternations can affect the functions of other bacteria and fungi through interactions between microorganisms. Correlation network analysis of this study revealed a significant correlation between some dramatically altered bacteria and fungi affected by diarrhea, which may further affect overall intestinal functions. Moreover, these changed bacteria or fungi can also affect some bacteria or fungi that do not significantly change during diarrhea through the interaction between microorganisms, thereby expanding the influence of diarrhea on the gut microbiota and intestinal functions. These results indicated that diarrhea not only directly altered the gut bacterial and fungal compositions and diversities but also indirectly affected some bacteria and fungi through the interaction of microorganisms, which may further destroy the intestinal homeostasis and induce intestinal dysfunction.
Taken together, this study dissected the shifts of gut bacterial and fungal communities in diarrheic yaks. Results revealed that the gut bacterial and fungal diversities of diarrheic yak were significantly decreased and accompanied by significant changes in taxonomic composition. This study also extended the understanding of gut bacterial and fungal characteristics in yaks with different health states and released a crucial message that the gut bacterial and fungal dysbiosis may be one of the causes of diarrhea in yaks. Furthermore, the present study may provide a theoretical basis for establishing diarrhea control system in yaks from the gut microbial perspective.
Availability of data and materials
Han H, Jiang Y, Wang M, Melaku M, Liu L, Zhao Y, et al. Intestinal dysbiosis in nonalcoholic fatty liver disease (NAFLD): focusing on the gut-liver axis. Crit Rev Food Sci Nutr. 2021. https://doi.org/10.1080/10408398.2021.1966738.
Li A, Wang Y, He Y, Liu B, Iqbal M, Mehmood K, et al. Environmental fluoride exposure disrupts the intestinal structure and gut microbial composition in ducks. Chemosphere. 2021;277: 130222.
Kayama H, Okumura R, Takeda K. Interaction between the microbiota, epithelia, and immune cells in the intestine. Annu Rev Immunol. 2020;38:23–48.
Huang P, Jiang A, Wang X, Zhou Y, Tang W, Ren C, et al. NMN maintains intestinal homeostasis by regulating the gut microbiota. Front Nutr. 2021;8: 714604.
Dong H, Liu B, Li A, Iqbal M, Mehmood K, Jamil T, et al. Microbiome analysis reveals the attenuation effect of lactobacillus from yaks on diarrhea via modulation of gut microbiota. Front Cell Infect Microbiol. 2020;10: 610781.
Liu Z, Li A, Wang Y, Iqbal M, Zheng A, Zhao M, et al. Comparative analysis of microbial community structure between healthy and Aeromonas veronii-infected Yangtze finless porpoise. Microb Cell Fact. 2020;19:123.
Manichanh C, Borruel N, Casellas F, Guarner F. The gut microbiota in IBD. Nat Rev Gastroenterol Hepatol. 2012;9:599–608.
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:10416.
Cao Q, Li R, Fu R, Zhang X, Yue B, Wang J, et al. Intestinal fungal dysbiosis in mice induced by fluoride. Chemosphere. 2020;245: 125617.
Qi M, Cao Z, Shang P, Zhang H, Hussain R, Mehmood K, et al. Comparative analysis of fecal microbiota composition diversity in Tibetan piglets suffering from diarrheagenic Escherichia coli (DEC). Microb Pathog. 2021;158: 105106.
Yan R, Andrew L, Marlow E, Kunaratnam K, Devine A, Dunican IC, et al. Dietary fibre intervention for gut microbiota, sleep, and mental health in adults with irritable bowel syndrome: a scoping review. Nutrients. 2021;13:2159.
Kanmani P, Suganya K, Kim H. The gut microbiota: how does it influence the development and progression of liver diseases. Biomedicines. 2020;8:501.
Zhong G, Wan F, Lan J, Jiang X, Wu S, Pan J, et al. Arsenic exposure induces intestinal barrier damage and consequent activation of gut-liver axis leading to inflammation and pyroptosis of liver in ducks. Sci Total Environ. 2021;788: 147780.
Zhou Y, Qiu W, Wang Y, Wang R, Takano T, Li X, et al. Beta-elemene suppresses obesity-induced imbalance in the microbiota-gut-brain axis. Biomedicines. 2021;9:704.
Karasova D, Crhanova M, Babak V, Jerabek M, Brzobohaty L, Matesova Z, et al. Development of piglet gut microbiota at the time of weaning influences development of postweaning diarrhea—a field study. Res Vet Sci. 2021;135:59–65.
Wang Y, Zhang H, Zhu L, Xu Y, Liu N, Sun X, et al. Dynamic distribution of gut microbiota in goats at different ages and health states. Front Microbiol. 2018;9:2509.
Gresse R, Chaucheyras-Durand F, Fleury MA, Van de Wiele T, Forano E, Blanquet-Diot S. Gut microbiota dysbiosis in postweaning piglets: understanding the keys to health. Trends Microbiol. 2017;25:851–73.
Kwon MS, Jo HE, Lee J, Choi KS, Yu D, Oh YS, et al. Alteration of the gut microbiota in post-weaned calves following recovery from bovine coronavirus-mediated diarrhea. J Anim Sci Technol. 2021;63:125–36.
Yang Q, Huang X, Zhao S, Sun W, Yan Z, Wang P, et al. Structure and function of the fecal microbiota in diarrheic neonatal piglets. Front Microbiol. 2017;8:502.
He K, Yan W, Sun C, Liu J, Bai R, Wang T, et al. Alterations in the diversity and composition of gut microbiota in weaned piglets infected with Balantioides coli. Vet Parasitol. 2020;288: 109298.
Zhu L, Xu F, Wan W, Yu B, Tang L, Yang Y, et al. Gut microbial characteristics of adult patients with allergy rhinitis. Microb Cell Fact. 2020;19:171.
Wang J, Zhu G, Sun C, Xiong K, Yao T, Su Y, et al. TAK-242 ameliorates DSS-induced colitis by regulating the gut microbiota and the JAK2/STAT3 signaling pathway. Microb Cell Fact. 2020;19:158.
Ding XZ, Liang CN, Guo X, Wu XY, Wang HB, Johnson KA, et al. Physiological insight into the high-altitude adaptations in domesticated yaks (Bos grunniens) along the Qinghai-Tibetan Plateau altitudinal gradient. Livest Sci. 2014;162:233–9.
Zhang L, Jiang X, Li A, Waqas M, Gao X, Li K, et al. Characterization of the microbial community structure in intestinal segments of yak (Bos grunniens). Anaerobe. 2020;61: 102115.
Qiu Q, Zhang G, Ma T, Qian W, Wang J, Ye Z, et al. The yak genome and adaptation to life at high altitude. Nat Genet. 2012;44:946–9.
Zhang Z, Xu D, Wang L, Hao J, Wang J, Zhou X, et al. Convergent evolution of rumen microbiomes in high-altitude mammals. Curr Biol. 2016;26:1873–9.
Rehman MU, Zhang H, Iqbal MK, Mehmood K, Huang S, Nabi F, et al. Antibiotic resistance, serogroups, virulence genes, and phylogenetic groups of Escherichia coli isolated from yaks with diarrhea in Qinghai Plateau, China. Gut Pathog. 2017;9:24.
Li K, Jiang J, Shahzad M, Hui Z, Mehmood K, Jiang X, et al. Revealing the parasitic infection in diarrheic yaks by piloting high-throughput sequencing. Microb Pathog. 2018;117:153–6.
Wang Y, Li A, Liu J, Mehmood K, Wangdui B, Shi H, et al. L. pseudomesenteroides and L. johnsonii isolated from yaks in Tibet modulate gut microbiota in mice to ameliorate enteroinvasive Escherichia coli-induced diarrhea. Microb Pathog. 2019;132:1–9.
Yang Y, Misra BB, Liang L, Bi D, Weng W, Wu W, et al. Integrated microbiome and metabolome analysis reveals a novel interplay between commensal bacteria and metabolites in colorectal cancer. Theranostics. 2019;9:4101–14.
Li A, Yang Y, Qin S, Lv S, Jin T, Li K, et al. Microbiome analysis reveals gut microbiota alteration of early-weaned Yimeng black goats with the effect of milk replacer and age. Microb Cell Fact. 2021;20:78.
Hu J, Nie Y, Chen J, Zhang Y, Wang Z, Fan Q, et al. Gradual changes of gut microbiota in weaned miniature piglets. Front Microbiol. 2016;7:1727.
Iapichino G, Callegari ML, Marzorati S, Cigada M, Corbella D, Ferrari S, et al. Impact of antibiotics on the gut microbiota of critically ill patients. J Med Microbiol. 2008;57:1007–14.
Larsen N, Vogensen FK, van den Berg FW, Nielsen DS, Andreasen AS, Pedersen BK, et al. Gut microbiota in human adults with type 2 diabetes differs from non-diabetic adults. PLoS ONE. 2010;5: e9085.
Li A, Liu B, Li F, He Y, Wang L, Fakhar-E-Alam KM, et al. Integrated bacterial and fungal diversity analysis reveals the gut microbial alterations in diarrheic giraffes. Front Microbiol. 2021;12: 712092.
Liu J, Wang HW, Lin L, Miao CY, Zhang Y, Zhou BH. Intestinal barrier damage involved in intestinal microflora changes in fluoride-induced mice. Chemosphere. 2019;234:409–18.
Baothman OA, Zamzami MA, Taher I, Abubaker J, Abu-Farha M. The role of gut microbiota in the development of obesity and diabetes. Lipids Health Dis. 2016;15:108.
Manco M, Putignani L, Bottazzo GF. Gut microbiota, lipopolysaccharides, and innate immunity in the pathogenesis of obesity and cardiovascular risk. Endocr Rev. 2010;31:817–44.
Li H, Chen X, Liu J, Chen M, Huang M, Huang G, et al. Ethanol extract of Centella asiatica alleviated dextran sulfate sodium-induced colitis: restoration on mucosa barrier and gut microbiota homeostasis. J Ethnopharmacol. 2021;267: 113445.
Suhr MJ, Hallen-Adams HE. The human gut mycobiome: pitfalls and potentials–a mycologist’s perspective. Mycologia. 2015;107:1057–73.
Aricha H, Simujide H, Wang C, Zhang J, Lv W, Jimisi X, et al. Comparative analysis of fecal microbiota of grazing Mongolian cattle from different regions in Inner Mongolia, China. Animals. 2021;11:1938.
Sun B, Wang X, Bernstein S, Huffman MA, Xia DP, Gu Z, et al. Marked variation between winter and spring gut microbiota in free-ranging Tibetan Macaques (Macaca thibetana). Sci Rep. 2016;6:26035.
Ransom-Jones E, Jones DL, Edwards A, McDonald JE. Distribution and diversity of members of the bacterial phylum Fibrobacteres in environments where cellulose degradation occurs. Syst Appl Microbiol. 2014;37:502–9.
Saye L, Navaratna TA, Chong J, O’Malley MA, Theodorou MK, Reilly M. The anaerobic fungi: challenges and opportunities for industrial lignocellulosic biofuel production. Microorganisms. 2021;9:694.
Wilken SE, Monk JM, Leggieri PA, Lawson CE, Lankiewicz TS, Seppala S, et al. Experimentally validated reconstruction and analysis of a genome-scale metabolic model of an anaerobic Neocallimastigomycota fungus. Msystems. 2021;6:e00002.
Xin J, Chai Z, Zhang C, Zhang Q, Zhu Y, Cao H, et al. Comparing the microbial community in four stomach of dairy cattle, yellow cattle and three yak herds in Qinghai-Tibetan Plateau. Front Microbiol. 2019;10:1547.
Shang Q, Shan X, Cai C, Hao J, Li G, Yu G. Dietary fucoidan modulates the gut microbiota in mice by increasing the abundance of Lactobacillus and Ruminococcaceae. Food Funct. 2016;7:3224–32.
Kakiyama G, Pandak WM, Gillevet PM, Hylemon PB, Heuman DM, Daita K, et al. Modulation of the fecal bile acid profile by gut microbiota in cirrhosis. J Hepatol. 2013;58:949–55.
Alvares-da-Silva MR, Oliveira CP, Fagan A, Longo L, Thoen RU, Yoshimura ZP, et al. Interaction of microbiome, diet, and hospitalizations between Brazilian and American patients with cirrhosis. Clin Gastroenterol Hepatol. 2021;20(4):930–40.
Ren Z, You W, Wu S, Poetsch A, Xu C. Secretomic analyses of Ruminiclostridium papyrosolvens reveal its enzymatic basis for lignocellulose degradation. Biotechnol Biofuels. 2019;12:183.
Miller TL, Currenti E, Wolin MJ. Anaerobic bioconversion of cellulose by Ruminococcus albus, Methanobrevibacter smithii, and Methanosarcina barkeri. Appl Microbiol Biotechnol. 2000;54:494–8.
Ducatelle R, Maes S, Haesebrouck F, Van Immerseel F, Vanderheyden S, Eeckhaut V. Anti-inflammatory effects of Butyricicoccus pullicaecorum: a normal inhabitant of the chicken intestinal tract. J Comp Pathol. 2010;143:319.
Konturek PC, Konturek K, Brzozowski T, Wojcik D, Magierowski M, Targosz A, et al. Participation of the intestinal microbiota in the mechanism of beneficial effect of treatment with synbiotic Syngut on experimental colitis under stress conditions. J Physiol Pharmacol. 2020. https://doi.org/10.26402/jpp.2020.3.03.
Cann I, Dodd D, Kiyonari S, Moon YH, Schroeder C, Nair S, et al. The genus prevotella, a resource of enzymes for hemicellulose degradation. Biophys J. 2010;981:210A.
Yue S, Zhao D, Peng C, Tan C, Wang Q, Gong J. Effects of theabrownin on serum metabolites and gut microbiome in rats with a high-sugar diet. Food Funct. 2019;10:7063–80.
Wang G, Zhang Y, Zhang R, Pan J, Qi D, Wang J, et al. The protective effects of walnut green husk polysaccharide on liver injury, vascular endothelial dysfunction and disorder of gut microbiota in high fructose-induced mice. Int J Biol Macromol. 2020;162:92–106.
Noureldein MH, Bitar S, Youssef N, Azar S, Eid AA. Butyrate modulates diabetes-linked gut dysbiosis: epigenetic and mechanistic modifications. J Mol Endocrinol. 2020;64:29–42.
Hamer HM, Jonkers D, Venema K, Vanhoutvin S, Troost FJ, Brummer RJ. Review article: the role of butyrate on colonic function. Aliment Pharmacol Ther. 2008;27:104–19.
Koike S, Kobayashi Y. Development and use of competitive PCR assays for the rumen cellulolytic bacteria: Fibrobacter succinogenes, Ruminococcus albus and Ruminococcus flavefaciens. Fems Microbiol Lett. 2001;204:361–6.
Pratt VC, Tappenden KA, McBurney MI, Field CJ. Short-chain fatty acid-supplemented total parenteral nutrition improves nonspecific immunity after intestinal resection in rats. JPEN J Parenter Enteral Nutr. 1996;20:264–71.
Li Q, Chen H, Zhang M, Wu T, Liu R. Altered short chain fatty acid profiles induced by dietary fiber intervention regulate AMPK levels and intestinal homeostasis. Food Funct. 2019;10:7174–87.
Gosiewski T, Salamon D, Szopa M, Sroka A, Malecki MT, Bulanda M. Quantitative evaluation of fungi of the genus Candida in the feces of adult patients with type 1 and 2 diabetes—a pilot study. Gut Pathog. 2014;6:43.
Dalrymple BP, Cybinski DH, Layton I, McSweeney CS, Xue GP, Swadling YJ, et al. Three Neocallimastix patriciarum esterases associated with the degradation of complex polysaccharides are members of a new family of hydrolases. Microbiology. 1997;143(Pt 8):2605–14.
Kim S, Shin DS, Lee T, Oh KB. Periconicins, two new fusicoccane diterpenes produced by an endophytic fungus Periconia sp. with antibacterial activity. J Nat Prod. 2004;67:448–50.
Urquhart EJ, Punja ZK. Hydrolytic enzymes and antifungal compounds produced by Tilletiopsis species, phyllosphere yeasts that are antagonists of powdery mildew fungi. Can J Microbiol. 2002;48:219–29.
Zhang H, Wang Z, Feng Y, Cui Q, Song X. Phytohormones as stimulators to improve arachidonic acid biosynthesis in Mortierella alpina. Enzyme Microb Technol. 2019;131: 109381.
Suresh Y, Das UN. Protective action of arachidonic acid against alloxan-induced cytotoxicity and diabetes mellitus. Prostaglandins Leukot Essent Fatty Acids. 2001;64:37–52.
The current research was supported by Establishment and application of a new mode for prevention and control of mixed infection of main animal diseases, Shandong Province agricultural major application technology innovation project (SD2019XM007); Evaluation test of animal disease detection reagents of Hangzhou Bori Technology Co., LTD (LYDX-BIOER-202011); Development and application of animal disease monitoring and early warning system (LYDX-SPRING-202107); Start-up fund of Nanjing Agricultural University (804131) and Start-up Fund for Distinguished Scholars of Nanjing Agricultural University (80900219).
Ethics approval and consent to participate
The animal study was reviewed and approved by the Ethics Committee of the Nanjing Agricultural University. Accession numbers: The original sequence data was submitted to the Sequence Read Archive (SRA) (NCBI, USA) with the accession no. PRJNA761293.
Consent for publication
There is no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1:
Fig. S1. Clustered heatmap of yaks in different health status at the bacterial genus level. The color values of the heatmap indicate the normalized relative richness of each species. Fig. S2. Clustered heatmap of yaks in different health status at the fungal genus level. The color values of the heatmap indicate the normalized relative richness of each species. Fig. S3. The gut bacterial comparisons between healthy and diarrheic yaks in phylum and genus levels. Metastats analysis was applied to identify the significantly differentially abundant bacterial genera between both groups and all of the data represent means ± SD. * p < 0.05, ** p < 0.01. Fig. S4. The gut fungal comparisons between healthy and diarrheic yaks in phylum and genus levels. Metastats analysis was applied to identify the significantly differentially abundant fungal genera between both groups and all of the data represent means ± SD. * p < 0.05, ** p < 0.01.
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.
About this article
Cite this article
Liu, J., Wang, X., Zhang, W. et al. Comparative analysis of gut microbiota in healthy and diarrheic yaks. Microb Cell Fact 21, 111 (2022). https://doi.org/10.1186/s12934-022-01836-y
- Tibet Plateau
- Gut microbiota