Analysis of gut microbiota in chinese donkey in different regions using metagenomic sequencing
BMC Genomics volume 24, Article number: 524 (2023)
Gut microbiota plays a significant role in host survival, health, and diseases; however, compared to other livestock, research on the gut microbiome of donkeys is limited.
In this study, a total of 30 donkey samples of rectal contents from six regions, including Shigatse, Changdu, Yunnan, Xinjiang, Qinghai, and Dezhou, were collected for metagenomic sequencing. The results of the species annotation revealed that the dominant phyla were Firmicutes and Bacteroidetes, and the dominant genera were Bacteroides, unclassified_o_Clostridiales (short for Clostridiales) and unclassified_f_Lachnospiraceae (short for Lachnospiraceae). The dominant phyla, genera and key discriminators were Bacteroidetes, Clostridiales and Bacteroidetes in Tibet donkeys (Shigatse); Firmicutes, Clostridiales and Clostridiales in Tibet donkeys (Changdu); Firmicutes, Fibrobacter and Tenericutes in Qinghai donkeys; Firmicutes, Clostridiales and Negativicutes in Yunnan donkeys; Firmicutes, Fibrobacter and Fibrobacteres in Xinjiang donkeys; Firmicutes, Clostridiales and Firmicutes in Dezhou donkeys. In the functional annotation, it was mainly enriched in the glycolysis and gluconeogenesis of carbohydrate metabolism, and the abundance was the highest in Dezhou donkeys. These results combined with altitude correlation analysis demonstrated that donkeys in the Dezhou region exhibited strong glucose-conversion ability, those in the Shigatse region exhibited strong glucose metabolism and utilization ability, those in the Changdu region exhibited a strong microbial metabolic function, and those in the Xinjiang region exhibited the strongest ability to decompose cellulose and hemicellulose.
According to published literature, this is the first study to construct a dataset with multi-regional donkey breeds. Our study revealed the differences in the composition and function of gut microbes in donkeys from different geographic regions and environmental settings and is valuable for donkey gut microbiome research.
Gut microflora is closely associated with host health and has achieved significant attention owing to its crucial role in the body’s metabolism, immunity, and physiological functions [1,2,3]. All equines are members of the family of herbivorous mammals and possess a specific hindgut (cecum and colon) microbiota, which can utilize forage for optimal nutrition. These microorganisms contribute to a significant portion of a horse’s daily energy needs by fermenting plant material into short-chain fatty acids such as acetate, propionate, and butyrate . The hindgut of a donkey shares some similarities with ruminants and monogastrics, like ruminants, the contents of the cecum and colon are combined, but unlike ruminants, donkeys do not regurgitate digestion for further chewing and do not have a complex forestomach to help digest plant carbohydrates and proteins . Donkeys are highly dependent on their gut microbiota for nutrition, as these provide glycoside hydrolases and polysaccharide lyases to digest complex polysaccharides that make up the bulk of their carbohydrate intake, and gut microbes convert “difficult” substrates such as cellulose, hemicellulose, and pectin into useful nutrients such as short-chain fatty acids .
Li et al.  used metagenomics to sequence horse intestinal samples and observed that a high abundance and diversity of antibiotic resistance genes (ARG) were identified in the myelin-associated glycoprotein, which demonstrated the horse gut as a reservoir of ARGs. Although donkeys and horses are closely related, their gut microbiomes markedly differ in composition and structure. On comparing the composition and function of the microbiomes of Przewalski’s proto antelope and Asiatic wild ass using a metagenomic shotgun sequencing method, Tang et al.  discovered that homologous wild horses have different microbial compositions while having stable microbial functional cores, which may help them survive in challenging habitats. Li et al.  observed a complex relationship between the gut microbiome and gene expression in their study on the hindgut microbial community of donkeys, especially in the immune system, and that the peroxisome proliferator-activated receptor pathway is mainly enriched in the cecum. The hindgut of donkeys shares some similarities with that of ruminants and monogastric animals, such as the combined contents of the cecum and colon. However, contrasting with ruminants, donkeys do not regurgitate the digesta for further chewing and have no complex foregut to aid the digestion of plant carbohydrates and proteins. Donkeys are highly dependent on their gut microbiota for nutrition since these microbiota secrete glycoside hydrolases (GH) and polysaccharide lyases to digest the complex polysaccharides, which account for the majority of their carbohydrate intake. “Difficult” substrates, such as hemicellulose and pectin, are converted into essential nutrients such as short-chain fatty acids . The gut microbiota constitute a complex ecosystem that is involved not only in digestive function  but also in immune function [12, 13], gut-brain connectivity and behavior , diabetes, and obesity [4, 15]. Although in some cases it can also become an enteric pathogen . China has various ecology types including the plateau continental climate, the temperate monsoon climate of the North China Plain, and the subtropical plateau monsoon climate, which comprises abundant donkey species resources. Particularly, Qinghai-Tibet Plateau is famous for its high-altitude, low-temperature, and hypoxic conditions. Moreover, the emergence of metagenomics has provided microbial ecology a technical boost, and human society has benefited from the progressive revelation of the roles of gut microbes.
In this study, metagenomic sequencing was used to conduct species annotation and gene annotation research on the donkey rectal microbiome in six regions and provide a data platform for the subsequent mining of functional genes. The importance of this paper lies in studying the differences in the composition and function of donkey gut microorganisms under different regional and environmental conditions, understanding the influence of intestinal flora on the metabolic homeostasis of donkeys, and providing important theoretical guidance for improving the microecological environment of donkey gastrointestinal tract and preventing and controlling the occurrence of diseases.
Materials and methods
Sample collection, DNA extraction and metagenomic sequencing
We collected five rectal contents of donkeys from each region, with a total of 20 samples from four areas, including Shigatse located in Tibet (RKZ, 5,000 m above sea level), Changdu located in Tibet (CD, 3,500 m above sea level), Yunnan (YN, 2,000 m above sea level), and Xinjiang (XJ, 1 200 m above sea level) (Fig. 1). All sampling subjects were adult healthy donkeys aged 5–10 years old with medium condition, feeding with local feed. All collected samples were guaranteed to be consistent with other environmental elements like temperature, humidity, etc. A rectal sampling method was employed. Long disposable gloves were worn to aseptically collect fresh rectal contents from the anus in 2-mL sterile Corning cryovials (Corning Incorporated, New York, USA) after cleaning, disinfection, and stabilization. The samples were labeled, quickly frozen in liquid nitrogen, kept in a -80 °C refrigerator.
DNA was extracted from the samples using the E.Z.N.A.® Soil DNA Kit following the manufacturer’s instructions and stored in a -80 ℃. NanoDrop2000 was used to detect DNA purity and concentration, and 1% agarose gel electrophoresis was performed to detect DNA integrity. The voltage was 5 V/cm, and the time was 20 min.
The extracted genomic DNA was submitted to a testing laboratory (Shanghai Majorbio Bio-pharm Technology Co., Ltd) for metagenomic sequencing. DNA fragments with an average size of approximately 400 bp were obtained using sonication. Paired-end sequencing was performed on an Illumina NovaSeq platform.
Additionally, we obtained metagenomic data from Qinghai (QH, 2 500 m above sea level) and Dezhou area (DZ, altitude ＜30 m) published in the previous article . In total, data from 30 samples from six regions were used for the subsequent analysis.
The software fastp  (https://github.com/OpenGene/fastp, version 0.20.0) was used to perform quality shearing on the adapter sequences at the 3′ and 5′ ends of the reads. The sheared length was < 50 bp; the average base quality value was < 20; the reads containing nitrogenous bases were removed to retain high-quality pair-end reads and single-end reads.
The optimized sequences were assembled using the splicing software MEGAHIT  (https://github.com/voutcn/megahit, version 1.1.2) based on the principle of succinct de Bruijn graphs. Contigs with a length of > 300 bp in the assembly results were screened as the final assembly results. Open Reading Frame (ORF) prediction of contigs in the splicing results was performed using MetaGene  (http://metagene.cb.k.u-tokyo.ac.jp/). Genes with nucleic acid lengths of ≥ 100 bp were selected and translated into amino acid (AA) sequences.
CD-HIT (http://www.bioinformatics.org/cd-hit/, version 4.6.1) was used to cluster the predicted gene sequences of all samples (parameters: 90% identity, 90% coverage), and the longest gene in each class was regarded as the representative sequence to construct a non-redundant (NR) gene set. Using the SOA Paligner  software (http://soap.genomics.org.cn/, version 2.21), the high-quality reads of each sample were separately aligned with the NR gene set (95% identity) and statistical gene abundance information in corresponding samples.
Species and functional annotations
The AA sequences of the NR gene sets were aligned with the NR database using Diamond  with an expected e-value of 1e-5, and species annotations were obtained through the taxonomic information database corresponding to the NR database. The sum of the gene abundances corresponding to the species was used to calculate the abundance of the species.
Diamond  was used to compare the eggNOG (version 4.5.1), Kyoto Encyclopedia of Genes and Genomes (KEGG) (version 94.2), and Carbohydrate-Active enZYme (CAZy) databases; the expected e-value was 1e-5 to obtain the corresponding function and calculate the abundance of functional categories.
Redundancy analysis (RDA) correlation analysis of environmental factors
RDA was a constrained Principal Component Analysis, which can reflect samples and environmental factors on the same two-dimensional ranking graph, from which we intuitively observed the relationship between sample distribution and environmental factors. It comprised the following three steps: (1) The Bioenv function was used to determine the maximum Pearson correlation coefficient between the environmental factors and the distribution difference of species/function of the sample, and the subset of environmental factors was obtained by the maximum correlation coefficient. (2) RDA analysis was performed between the species/function distribution table of the sample and environmental factors or their subsets. (3) The significance of RDA analysis was determined by Permutest analysis, which is similar to the Analysis of Variance.
After filtering, a total of 2,341,055,166 reads and 353,195,288,355 bp were obtained. After the assembly, a total of 28,097,975 contigs were received with an average length of 21,921,556,352 bp and N50 of 26,337 bp. Finally, 198,223,248,57 open reading frames were generated. Afterward, analysis was performed using the NR, eggNOG, KEGG, and CAZymes databases.
Compositional and differential analysis based on species abundance
At the phylum level, 227 common phyla were identified (Fig. 2A), and Firmicutes and Bacteroidetes were the dominant phyla in all intestinal sites (40.19% and 28.03% on average, respectively) (Fig. 2B). The dominant phyla were Bacteroidetes in RKZ and Firmicutes in other five donkey populations. Firmicutes accounted for the highest proportion in DZ; Bacteroidetes were the most abundant in RKZ. The abundance of Fibrobacteres rectal contents was significantly higher in XJ than that in DZ (P < 0.05). Spirochaetes had the highest abundance in XJ, which was significantly higher than that in RKZ (P < 0.05), CD (P < 0.01), and DZ (P < 0.01) (Fig. 3A).
At the genus level, 4,725 common genera were identified (Fig. 2D), and the dominant genera were Clostridiales and unclassified_o_ Bacteroidales (short for Bacteroidales) in all intestinal sites (8.39% and 7.43% on average, respectively) (Fig. 2E). The dominant genera were Clostridiales in RKZ, CD, YN and DZ; Fibrobacter in QH and XJ. The abundance of rectal contents of Clostridiales was the highest in CD, which was significantly higher than that in QH (P < 0.05) and XJ (P < 0.05). The abundance of Fibrobacter’s rectal contents was the highest in XJ, which was significantly higher than that in DZ (P < 0.05) (Fig. 3B).
It was worth noting that the proportions of Fibrobacteres in the XJ and QH regions and Spirochaetes in the QH, YN, and XJ regions were higher than those in other regions; Lachnospiraceae had the highest proportion in the QH region (Fig. 2C, F). As presented in the Latent Dirichlet Allocation (LDA) diagram in Fig. 3C, the key discriminators were Bacteroidetes in RKZ, Clostridiales in CD, Tenericutes in QH, Negativicutes in YN, Fibrobacteres in XJ, and Firmicutes in DZ. These results demonstrated significant differences in the microbiota composition among donkeys from six regions.
Compositional and Differential Analysis based on functional abundance
Among the clusters of orthologous genes (COG) functions, 24 functions were shared by the six groups (Fig. 4A). After removing the function unknown [S] (Function unknown, 30.12%), the top two with the highest abundance were replication, recombination, and repair [L] (9.93%) and carbohydrate transport and metabolism [G] (8.49%) (Fig. 4D). The functional abundance of [L] was highest in CD. The abundance of the [G] function was the highest in RKZ, and among this function, Lachnospiraceae had the highest contribution in the QH group. The functional abundance of cell wall/membrane/envelope biogenesis [M] was highest in DZ. Among [M] functions, Bacteroidales and Bacteria, Clostridiales, and Fibrobacter had high contributions in the DZ, CD, and XJ groups, respectively (Fig. 5A, D).
In the KEGG function, 424 functions were shared by the six groups (Fig. 4B), of which the top two with the highest abundance were metabolic pathways (17.76%) and biosynthesis of secondary metabolites (8.09%) (Fig. 4E). The functional abundance of glycolysis/gluconeogenesis was highest in DZ and significantly higher than that in XJ (P < 0.05); there were no significant differences in the remaining groups. The functional abundance of galactose metabolism was highest in RKZ, which was significantly higher than that of QH and DZ (P < 0.01) and XJ (P < 0.05). Pentose phosphate pathway abundance was highest in DZ, which was significantly higher than that in RKZ (P < 0.05) (Fig. 5B). Additionally, Clostridiales exhibited a high contribution to the function of microbial metabolism in diverse environments and CD groups. Bacteroidales presented a high contribution to the function of biosynthesis of secondary metabolites and DZ. Verrucomicrobia contributed more to the biosynthesis of secondary metabolites and RKZ. Bacteria had higher microbial metabolism in diverse environments and contributed more in DZ; Lachnospiraceae had a higher contribution in the biosynthesis of AA and QH; the contribution of Fibrobacter was significantly higher in the biosynthesis of AA and the XJ group than in other groups (Fig. 5E).
In the CAZy analysis, the six groups of gut microbiota covered six major enzymes and shared six functions (Fig. 4C), of which the first two with the highest abundance were glycoside hydrolase (GH) and glycosyl transferase (GT) (Fig. 4F).
At the class level, the functional abundance of GH was the highest in RKZ, which was significantly higher than that of QH (P < 0.01) and DZ (P < 0.001). Bacteroides contributed more to GH function and donkey’s gut microbes in the RKZ region. The functional abundance of GT was the highest in DZ, which was significantly higher than that in CD (P < 0.05), YN (P < 0.05), and RKZ (P < 0.001), and Clostridium had the highest contribution in GT function and donkey’s gut microbes in the QH region. The functional abundance of carbohydrate-binding modules (CBM) was significantly higher in XJ than that in DZ (P < 0.01) and QH than that in DZ (P < 0.05), respectively. The contribution of Fibrobacter in the CBM function and XJ group was significantly higher than that in other groups (Fig. 5 C, F).
Based on the KEGG annotation results, the differential detection and visual analysis of differentially abundant enzymes were performed for a certain pathway. At Level 2, metabolism had the highest functional enrichment, and the top three functional abundances were global and overview maps, carbohydrate metabolism, and AA metabolism (Fig. 6A). The most abundant carbohydrate pathway was glycolysis/gluconeogenesis (Fig. 6B), and the differences in their metabolic pathways are presented in Fig. 6B C. Among them, the abundance of pyruvate ferredoxin oxidoreductase (EC: 220.127.116.11, including four subunits, porA, porB, porG, and porD), in QH was the highest, which was higher than that in CD, YN, and DZ (P < 0.05); the abundance of pyruvate ferredoxin oxidoreductase in RKZ was higher than that in CD (P < 0.05) and DZ (P < 0.01). The abundance of 2-oxoglutarate/2-oxoacid ferredoxin oxidoreductase subunit alpha (EC: 18.104.22.168, including korA and korB subunits) in RKZ was the highest, which was higher than that in XJ (P < 0.01), QH, and YN (P < 0.05). The abundance of Aldose 1-epimerase (EC: 22.214.171.124, including galM subunits) in RKZ was the highest, which was higher than that in CD, QH, YN, and DZ (P < 0.05). The abundance of Acetyl-CoA synthetase (EC: 126.96.36.199, including two subunits of ACSS and AAE7) was highest in QH and the lowest in DZ, which was significantly lower than that in CD, YN, and XJ (P < 0.05), RKZ (P < 0.001), and QH (P < 0.01).
Altitude Association Analysis based on species and functional abundance
The top five species or functions in the samples from six different areas were examined for connection with altitude using RDA analysis. At the species level, Bacteroidetes was closest to RKZ (5 000 m) in the area with the highest altitude (Fig. 7A). In the COG function, the function of carbohydrate transport and metabolism [G] was the highest in the RKZ (5,000 m) group, and the maximum functional abundance of replication, recombination, and repair [L] was observed in the CD (3,500 m) (Fig. 7B). In the KEGG analysis, the function abundance of the metabolic pathway exhibited the maximum correlation with the RKZ (5,000 m) region (Fig. 7C). At the class level in the CAZy function, the RKZ (5,000 m) region had the highest correlation with the function of the GH (Fig. 7D).
Analysis of gut microbiota composition based on species abundance
Firmicutes, followed by Bacteroidetes, was the most prevalent phylum at the phylum level in the six regions. The maximum abundance of Firmicutes, Bacteroidetes, and Fibrobacteres were observed in DZ donkeys, RKZ, and XJ, respectively. Previous studies have shown that Firmicutes and Bacteroidetes are the dominant microbial communities in the gut of horses . Research by M Denise Dearing et al. showed that Firmicutes is the main microbial phylum that promotes the breakdown of fibers in the gastrointestinal tract of herbivores . Cheryl Spence et al. discovered that Bacteroidetes is the main microbial phylum for carbohydrate metabolism in herbivores . A study by Wu et al. demonstrated that Fibrobacteres microorganisms can facilitate the digestion of cellulose and hemicellulose in the diet . This has also been supported by previous studies on the diversity of mammalian gut microbia [29, 30].
The host’s diet and mucosal secretions were mostly sources of polysaccharides in the gut. Most Firmicutes produce butyric acid, which is not absorbed by the intestine, and one of the final metabolites of polysaccharides is also butyric acid. Butyric acid was crucial for physiological activities such as providing energy and promoting the development of intestinal epithelial cells. Although ATP-binding transporters can be employed for carbohydrate conversion and transport, Firmicutes encode fewer carbohydrate-degrading enzymes . Bacteroides can break down carbohydrates in food and provide about one-tenth of the energy to the host. It had been reported that Bacteroides can also resist the invasion of pathogens and enhance the body’s resistance by colonizing the intestinal mucosal epithelium . About 20% of the genes in Bacteroidetes are used for the decomposition of sugars, which can efficiently degrade polysaccharides and produce a significant amount of short-chain fatty acids, was one of the most studied bacteria on polysaccharides .
At the genus level, the most abundant genus was Clostridiales, followed by Bacteroidales. Clostridiales and Bacteroidales had the highest abundance in CD and RKZ, respectively. Clostridiales is a cellulose-degrading bacteria capable of breaking down cellulose and providing fermentation substrates for acid-producing bacteria [34, 35]. Clostridiales species act as fiber fermenters to produce acetate and propionate from succinic acid produced by other bacteria [36, 37]. By improving the gastrointestinal flora, Clostridiales can stabilize the environment and intestinal pH, provide a favorable environment for the growth of gastrointestinal anaerobic flora, enhance the decomposition, digestion, absorption, and utilization of the gastrointestinal tract, improve the growth performance of animals, and enhance the body’s immune function. 16 S ribosomal RNA gene sequencing revealed that cocci bacteria and Clostridium were dominant in the rumen of Tibetan antelope, Tibetan wild donkey, and Tibetan sheep . As ubiquitous hindgut symbionts in humans and other animals, many species of Bacteroidetes help digest various plant polysaccharides, including fibrous matter [39, 40]. In addition to their degradative activity, Bacteroides also play a crucial role in developing epithelial immunity and maintaining intestinal microecological balance . Compared with that in the low-altitude wild mice, the proportion of Bacteroides in the intestinal tract of the plateau wild mice was significantly increased . These results indicated that donkeys in the CD area had strong microbial metabolic function and those in the Shigatse region had strong glucose metabolism and utilization ability to adapt to the plateau environment.
Analysis of gut microbiota composition based on functional abundance
The results of COG functional annotation revealed that replication, recombination, and repair [L], carbohydrate transport and metabolism [G], and cell wall/membrane/envelope biogenesis [M] had the highest functional abundance in CD, RKZ, and DZ, respectively. The main pathways annotated in the KEGG database include metabolic pathways and biosynthesis of secondary metabolites, among which glycolysis/gluconeogenesis had the highest functional abundance.
GH and GT were the most abundant among the CAZy functional annotations. GH accounted for the highest proportion in RKZ, and according to the LDA discrimination results, the GH2 function had a significant effect on RKZ. The GH2 family contains a large number of glycosidases that cleave oligosaccharides and non-reducing carbohydrates in the side chains of hemicellulose and pectin. Furthermore, in the KEGG functional annotation, the top two including glycan degradation and galactose metabolism was located in RKZ, indicating that RKZ samples exhibited a strong sugar metabolism. Among various carbohydrate enzymes, GHs, including glycosidases and transglycosidases, which were the most abundant and diverse group of enzymes responsible for the hydrolysis or transglycosylation of glycosidic bonds, accounted for half of the enzymes in the CAZymes database , and these can break glycosidic bonds of polysaccharides in plants . Studies have confirmed that there are 72 different GHs families in buffalo and 78 different GHs in the rumen of cattle, which are responsible for the hydrolysis of various cell wall components . In the CAZy database annotations, GT accounted for the highest proportion in DZ. Previous studies have reported that GT was responsible for catalyzing the cleavage of glycosidic bonds, accounting for 19–24% of the total CAZyme, and was the second most abundant enzyme in the CAZy family . Carbohydrate-binding molecules exhibit no enzymatic activity but can bind GHs, CEs, and AAs to polysaccharides, thus, increasing their activity . CEs can remove the ester group modification present in polysaccharides, thus, promoting the effect of GHs on complex polysaccharides. This may be associated with the fact that donkeys in RKZ have a unique intestinal flora to adapt to the low temperature, hypoxia, and high altitude of the plateau. Meanwhile, the correlation analysis based on environmental factors (Fig. 7) revealed that the altitude of the RKZ area had the highest correlation with the donkey’s intestinal microbes. It exhibited that the composition and function of intestinal flora will be affected by the severe environment conditions such as low temperature and hypoxia in high altitude area, however, the specific mechanism of action still needs additional research and data support.
The glycolysis/gluconeogenesis pathways for carbohydrates were enriched in KEGG, with pyruvate ferredoxin oxidoreductase being the most abundant in QH. 2-oxoglutarate/2-oxoacid ferredoxin oxidoreductase subunit alpha and aldose 1-epimerase were most abundant in RKZ, and Acetyl-CoA synthetase is most abundant in QH. All tested archaea contained pyruvate ferredoxin oxidoreductase, which is involved in either catabolism or anabolism. Archaea are single-celled microorganisms that do not have a nucleus and any other membrane-bound organelles, and in the metabolic process, there are many special coenzymes, most of which are strictly anaerobic, facultative anaerobic, and obligate aerobic . Thus, pyruvate ferredoxin oxidoreductase represents the only mechanism for pyruvate acetyl-CoA conversion in the archaeal domain . Aldose 1-epimerase or mutarotase was a key enzyme in carbohydrate metabolism that catalyzes the interaction between α-D-glucose and β-D-glucose and their mutual transformation. This enzyme is necessary for normal lactose metabolism in E. coli. Strains, in which the gene that had been deleted grew more slowly on a medium containing any sugar. Mutarotases have been reported to convert various sugars, including D-glucose, D-fucose, D-quinose, L-arabinose, and D-xylose. A crucial metabolite in the metabolism of carbon and energy is acetyl-CoA. In mammalian cells, there are many ways for carbohydrates to be converted into various biomolecules. One of these included the breakdown of carbohydrates into acetyl-CoA, which was subsequently employed as a precursor for anabolism. This was a crucial technique for transforming sugars into other biomolecules. The three main pathways of mitochondrial acetyl-CoA synthesis include the oxidative decarboxylation of glycolysis product pyruvate, fatty acid β-oxidation, and branched-chain AA decomposition. Mitochondrial acetyl-CoA normally progresses to the tricarboxylic acid cycle for further metabolism . This indicated that donkeys in high-altitude regions had relatively strong carbohydrate metabolism and sugar conversion, which can fulfil the consumption of basal metabolism in the high-altitude environment and improve hypoxia endurance to help them adapt.
In this study, we analyzed the gut flora of six donkey populations by metagenomic high-throughput sequencing. We observed that Firmicutes and Bacteroidetes were the dominant phyla, and the dominant genera were Clostridiales and Bacteroidales. Additionally, the donkey flora structure and function had significant differences in different areas. We hope that this study will serve as a specific reference for understanding the gut microbiome structure of donkeys in China and close a research gap in donkey gut microbiome research.
The datasets generated in this study have been uploaded to National Center for Biotechnology Information repositories and will be made available after the article is received. Shigatse, Changdu, Yunnan, and Xinjiang areas of donkey gut microbe data repository/repository name and registration number can be found on the following site: https://www.ncbi.nlm.nih.gov/, BioProject ID: PRJNA894235. Qinghai and Dezhou areas of donkey gut microbe data repository/repository name and registration number can be found on the following site: https://www.ncbi.nlm.nih.gov/, BioProject ID: PRJNA843204.
Turnbaugh PJ, Ley RE, Mahowald MA, Magrini V, Mardis ER, Gordon JI. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature. 2006;444(7122):1027–31.
Kamada N, Kim Y-G, Sham HP, Vallance BA, Puente JL, Martens EC, Núñez G. Regulated virulence controls the ability of a pathogen to compete with the gut microbiota. Science. 2012;336(6086):1325–9.
Osbelt L, Thiemann S, Smit N, Lesker TR, Schröter M, Gálvez EJ, Schmidt-Hohagen K, Pils MC, Mühlen S, Dersch P. Variations in microbiota composition of laboratory mice influence Citrobacter rodentium infection via variable short-chain fatty acid production. PLoS Pathog. 2020;16(3):e1008448.
Biddle AS, Tomb JF, Fan Z. Microbiome and blood analyte differences point to Community and metabolic signatures in lean and obese horses. Front Vet Sci. 2018;5:225.
Dehority BA. Gastrointestinal tracts of herbivores, particularly the ruminant: anatomy, physiology and microbial digestion of plants. J Appl Anim Res. 2002;21(2):145–60.
Ley RE, Lozupone CA, Hamady M, Knight R, Gordon JI. Worlds within worlds: evolution of the vertebrate gut microbiota. Nat Rev Microbiol. 2008;6(10):776–88.
Li C, Li X, Guo R, Ni W, Liu K, Liu Z, Dai J, Xu Y, Abduriyim S, Wu Z, et al. Expanded catalogue of metagenome-assembled genomes reveals resistome characteristics and athletic performance-associated microbes in horse. Microbiome. 2023;11(1):7.
Tang L, Gao Y, Yan L, Jia H, Chu H, Ma X, He L, Wang X, Li K, Hu D et al. Comparative analysis of Microbiome Metagenomics in Reintroduced Wild Horses and Resident Asiatic Wild Asses in the Gobi Desert Steppe. Microorganisms 2022, 10(6).
Li Y, Ma Q, Shi X, Liu G, Wang C. Integrated multi-omics reveals novel microbe-host lipid metabolism and immune interactions in the donkey hindgut. Front Immunol. 2022;13:1003247.
Ley RE, Hamady M, Lozupone C, Turnbaugh PJ, Ramey RR, Bircher JS, Schlegel ML, Tucker TA, Schrenzel MD, Knight R, et al. Evolution of mammals and their gut microbes. Science. 2008;320(5883):1647–51.
Lerner A, Aminov R, Matthias T. Dysbiosis May trigger Autoimmune Diseases via Inappropriate Post-Translational modification of host proteins. Front Microbiol. 2016;7:84.
Manfredo VS, Hiltensperger M, Kumar V, Zegarra-Ruiz D, Dehner C, Khan N, Costa FRC, Tiniakou E, Greiling T, Ruff W, et al. Translocation of a gut pathobiont drives autoimmunity in mice and humans. Science. 2018;359(6380):1156–61.
Paun A, Danska JS. Immuno-ecology: how the microbiome regulates tolerance and autoimmunity. Curr Opin Immunol. 2015;37:34–9.
Cryan JF, Dinan TG. Mind-altering microorganisms: the impact of the gut microbiota on brain and behaviour. Nat Rev Neurosci. 2012;13(10):701–12.
Morrison PK, Newbold CJ, Jones E, Worgan HJ, Grove-White DH, Dugdale AH, Barfoot C, Harris PA, Argo CM. The equine gastrointestinal microbiome: impacts of age and obesity. Front Microbiol. 2018;9:3017.
Hornef M. Pathogens, commensal Symbionts, and Pathobionts: Discovery and Functional Effects on the host. ILAR J. 2015;56(2):159–62.
Guo R, Zhang S, Chen J, Shen W, Zhang G, Wang J, Zhang F, Pan Q, Xie T, Ai D. Comparison of gut microflora of donkeys in high and low altitude areas. Front Microbiol 2022, 13.
Chen S, Zhou Y, Chen Y, Gu J. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34(17):i884–90.
Li D, Liu CM, Luo R, Sadakane K, Lam TW. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de bruijn graph. Bioinformatics. 2015;31(10):1674–6.
Noguchi H, Park J, Takagi T. MetaGene: prokaryotic gene finding from environmental genome shotgun sequences. Nucleic Acids Res. 2006;34(19):5623–30.
Fu L, Niu B, Zhu Z, Wu S, Li W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 2012;28(23):3150–2.
Li R, Li Y, Kristiansen K, Wang J. SOAP: short oligonucleotide alignment program. Bioinformatics. 2008;24(5):713–4.
Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12(1):59–60.
Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51(D1):D587–d592.
Al Jassim RA, Andrews FM. The bacterial community of the horse gastrointestinal tract and its relation to fermentative acidosis, laminitis, colic, and stomach ulcers. Veterinary Clin North Am Equine Pract. 2009;25(2):199–215.
Dearing MD, Kohl KD. Beyond fermentation: other important services provided to endothermic herbivores by their gut microbiota. Integr Comp Biol. 2017;57(4):723–31.
Spence C, Wells WG, Smith CJ. Characterization of the primary starch utilization operon in the obligate anaerobe Bacteroides fragilis: regulation by carbon source and oxygen. J Bacteriol. 2006;188(13):4663–72.
Wu X, Zhang H, Chen J, Shang S, Wei Q, Yan J, Tu X. Comparison of the fecal microbiota of dholes high-throughput Illumina sequencing of the V3–V4 region of the 16S rRNA gene. Appl Microbiol Biotechnol. 2016;100(8):3577–86.
Mariat D, Firmesse O, Levenez F, Guimarăes V, Sokol H, Doré J, Corthier G, Furet J. The Firmicutes/Bacteroidetes ratio of the human microbiota changes with age. BMC Microbiol. 2009;9(1):1–6.
Middelbos IS, Vester Boler BM, Qu A, White BA, Swanson KS, Fahey GC Jr. Phylogenetic characterization of fecal microbial communities of dogs fed diets with or without supplemental dietary fiber using 454 pyrosequencing. PLoS ONE. 2010;5(3):e9768.
Mahowald MA, Rey FE, Seedorf H, Turnbaugh PJ, Fulton RS, Wollam A, Shah N, Wang C, Magrini V, Wilson RK. Characterizing a model human gut microbiota composed of members of its two dominant bacterial phyla. Proc Natl Acad Sci. 2009;106(14):5859–64.
Marcobal A, Barboza M, Sonnenburg ED, Pudlo N, Martens EC, Desai P, Lebrilla CB, Weimer BC, Mills DA, German JB. Bacteroides in the infant gut consume milk oligosaccharides via mucus-utilization pathways. Cell host microbe. 2011;10(5):507–14.
Singh RP. Glycan utilisation system in Bacteroides and Bifidobacteria and their roles in gut stability and health. Appl Microbiol Biotechnol. 2019;103(18):7287–315.
Perea K, Perz K, Olivo SK, Williams A, Lachman M, Ishaq SL, Thomson J, Yeoman CJ. Feed efficiency phenotypes in lambs involve changes in ruminal, colonic, and small-intestine-located microbiota. J Anim Sci. 2017;95(6):2585–92.
Han R, Zhu D, Xing J, Li Q, Li Y, Chen L. The effect of temperature fluctuation on the microbial diversity and community structure of rural household biogas digesters at Qinghai Plateau. Arch Microbiol. 2020;202(3):525–38.
Connors J, Dawe N, Van Limbergen J. The role of succinate in the regulation of intestinal inflammation. Nutrients. 2018;11(1):25.
Kauter A, Epping L, Semmler T, Antao E-M, Kannapin D, Stoeckle SD, Gehlen H, Lübke-Becker A, Günther S, Wieler LHJAM. The gut microbiome of horses: current research on equine enteral microbiota and future perspectives. 2019, 1(1):1–15.
Ma Y, Ma S, Chang L, Wang H, Ga Q, Ma L, Bai Z, Shen Y, Ge R-L. Gut microbiota adaptation to high altitude in indigenous animals. Biochem Biophys Res Commun. 2019;516(1):120–6.
Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, Nielsen T, Pons N, Levenez F, Yamada T, et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010;464(7285):59–65.
Zhang H, Shao M, Huang H, Wang S, Ma L, Wang H, Hu L, Wei K, Zhu R. The dynamic distribution of small-tail han sheep microbiota across different intestinal segments. Front Microbiol. 2018;9:32.
Zhang W, Li N, Tang X, Liu N, Zhao W. Changes in intestinal microbiota across an altitudinal gradient in the lizard Phrynocephalus vlangalii. Ecol Evol. 2018;8(9):4695–703.
Suzuki TA, Martins FM, Nachman MW. Altitudinal variation of the gut microbiota in wild house mice. Mol Ecol. 2019;28(9):2378–90.
Luo J, Huang W, Guo W, Ge R, Zhang Q, Fang F, Feng Q, Cao J, Wu Y. Novel strategy to stimulate the food wastes anaerobic fermentation performance by eggshell wastes conditioning and the underlying mechanisms. Chem Eng J. 2020;398:125560.
Cantarel BL, Coutinho PM, Rancurel C, Bernard T, Lombard V, Henrissat B. The carbohydrate-active EnZymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res. 2009;37(suppl1):D233–8.
Jose VL, Appoothy T, More RP, Arun AS. Metagenomic insights into the rumen microbial fibrolytic enzymes in indian crossbred cattle fed finger millet straw. AMB Express. 2017;7(1):1–11.
Mota TR, Oliveira D, Marchiosi R, Ferrarese-Filho O, Santos W. Plant cell wall composition and enzymatic deconstruction. AIMS Bioeng. 2018;5(1):63–77.
Bernardes A, Pellegrini V, Curtolo F, Camilo C, Mello B, Johns M, Scott J, Guimaraes F, Polikarpov I. Carbohydrate binding modules enhance cellulose enzymatic hydrolysis by increasing access of cellulases to the substrate. Carbohydr Polym. 2019;211:57–68.
Baker BJ, De Anda V, Seitz KW, Dombrowski N, Santoro AE, Lloyd KG. Diversity, ecology and evolution of Archaea. Nat Microbiol. 2020;5(7):887–900.
Ikeda T, Ochiai T, Morita S, Nishiyama A, Yamada E, Arai H, Ishii M, Igarashi Y. Anabolic five subunit-type pyruvate: ferredoxin oxidoreductase from Hydrogenobacter thermophilus TK-6. Biochem Biophys Res Commun. 2006;340(1):76–82.
Boroughs LK, DeBerardinis RJ. Metabolic pathways promoting cancer cell survival and growth. Nat Cell Biol. 2015;17(4):351–9.
All authors reviewed the article and approved the submitted version. We also thank all the individuals who contributed to this article and the colleges and schools for providing their platforms. The authors would like to thank all the reviewers who participated in the review, as well as MJEditor (www.mjeditor.com) for providing English editing services.
This research was funded by the Donkey Innovation Team of Shandong Modern Agricultural Industry Technology System, grant number SDAIT-27; Major Agricultural Application Technology Innovation Projects of Shandong Province, grant number SD2019 XM 008; Shandong Province Natural Science Foundation, grant number ZR2022QC091; Experimental Technology Research Program of Qingdao Agriculture University, grant number 20210021; Hebei Provincial Science and Technology Planning Project, grant number 225A6601D. “Systematic Evaluation and Screening of donkey germplasm Resources in the Yellow River Basin”, grant number 2022DZ01; Autonomous region science and technology branch Xinjiang project, grant number 2021E02035. The authors declare that the research was conducted in the absence of any commercial or personal relationships that could be construed as a potential conflict of interest.
Ethics approval and consent to participate
This study was approved by the Animal Care and Use Committee of Qingdao Agricultural University (Qingdao, Shandong Province, China).And This study was reported in accordance with ARRIVE guidelines.
Consent for publication
The authors declare that the research was conducted in the absence of any commercial or personal relationships that could be construed as a potential conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Guo, R., Zhang, W., Shen, W. et al. Analysis of gut microbiota in chinese donkey in different regions using metagenomic sequencing. BMC Genomics 24, 524 (2023). https://doi.org/10.1186/s12864-023-09575-z
- Chinese donkey
- Gut microbes