Gut microbiota composition and metabolomic profiles of wild and captive Chinese monals (Lophophorus lhuysii)
Frontiers in Zoology volume 17, Article number: 36 (2020)
The Chinese monal (Lophophorus lhuysii) is an endangered bird species, with a wild population restricted to the mountains in southwest China, and only one known captive population in the world. We investigated the fecal microbiota and metabolome of wild and captive Chinese monals to explore differences and similarities in nutritional status and digestive characteristics. An integrated approach combining 16S ribosomal RNA (16S rRNA) gene sequencing and ultra-high performance liquid chromatography (UHPLC) based metabolomics were used to examine the fecal microbiota composition and the metabolomic profile of Chinese monals.
The results showed that the alpha diversity of gut microbes in the wild group were significantly higher than that in the captive group and the core bacterial taxa in the two groups showed remarkable differences at phylum, class, order, and family levels. Metabolomic profiling also revealed differences, mainly related to galactose, starch and sucrose metabolism, fatty acid, bile acid biosynthesis and bile secretion. Furthermore, strong correlations between metabolite types and bacterial genus were detected.
There were remarkable differences in the gut microbiota composition and metabolomic profile between wild and captive Chinese monals. This study has established a baseline for a normal gut microbiota and metabolomic profile for wild Chinese monals, thus allowing us to evaluate if differences seen in captive organisms have an impact on their overall health and reproduction.
The Chinese monal (Lophophorus lhuysii) belongs to the order Galliformes, family Phasianidae, distributed in the mountains of southwest China at an elevation of 3000 to 4900 m . It is an endemic bird species of China, which has been listed as endangered by the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES) , classified as a vulnerable species on the International Union for Conservation of Nature (IUCN) red list . There has been limited success with establishing captive breeding groups of Chinese monals throughout the world, the reason for which is not entirely clear . At present, the world’s only captive collection exists in the Fengtongzhai National Nature Reserve in Sichuan, China, with only 23 individuals as of early 2020. The IUCN states the wild population of Chinese monal continues to decrease with research needed to understand their ecology and threats to their livelihood (Available at https://www.iucnredlist.org/species/22679192/30181918). However, given the limited access to the species and their biological samples, few studies have been conducted in captive and wild Chinese monals . With a declining population, it is imperative to learn more about the Chinese monal to ensure a healthy existence of the species.
The gut microbiota is recognized as a coevolutionary partner with an ecosystem that is metabolically adaptable, rapidly renewable, and metabolically flexible [5, 6]. It performs numerous beneficial functions for the host, such as nutrient acquisition, immunomodulation and physiogenesis in response to profound lifestyle changes [7, 8]. The composition and diversity of gut microbiota is influenced by numerous factors, such as diet composition , social interactions  and health status . For example, wildlife that is maintained in a captive setting have significant differences in diet, social structure, and stress, all of which potentially affect their gut microbiota. Previous studies have demonstrated that captivity affects the composition of gut microbes in birds [11,12,13,14] and mammals [9, 10, 15, 16], with considerable differences were observed between wild and captive individuals. Thus, study on the gut microbiota in both wild and captive animals may provide a way to understand more about how to successfully maintain captive wildlife and in optimal health condition.
The fecal metabolomic profile is the product of functional activity of both host cells and gut bacteria, which provides information on their combined past activity . Any differences in the microbial community may have a significant effect on the metabolite profiles of the host , which can be explored through untargeted metabolomics. The combined study of fecal 16S ribosomal RNA (16S rRNA) gene sequencing and untargeted metabolomics is expected to uncover the inherent associations between the microbiota and the host metabolic phenotypes , and have been used to elucidate the differences in gut microbes and host metabolic phenotypes in the context of different food composition [19, 20].
Given the critical role of gut microbiota to the host and its close relationship with the host’s metabolic phenotype, the study of gut microbial composition and metabolomic profile could be helpful to understand the health and nutritional status of the captive Chinese monal (CCM) in relation to the wild Chinese monal (WCM). Therefore, the objectives of the present study were to investigate the fecal microbiota and metabolome of WCM and CCM, and discuss the causes for the differences. The study may provide a theoretical basis for the future breeding of CCM, a species whose wild population is declining and for which there are limited scientific publications regarding species health.
Monal species identification
There were six and 35 samples of the feces in Yuancaopo and Hongshanding belonging to the Chinese monal, respectively, according to the species identification. The remaining samples were identified as 24 blood pheasants (Ithaginis cruentus), six temminck’s tragopan (Tragopan temminckii), and 20 samples without amplification results.
Microbial community composition
A total of 1,466,455 high-quality reads with an average sequencing depth of 35,767 ± 1265 reads per WCM fecal sample was conducted and classified into 2132 operational taxonomic units (OTUs). Whereas, 580,736 high-quality reads with an average sequencing depth of 36,296 ± 1062 reads per CCM fecal sample was conducted and classified into 4899 OTUs. The number of OTUs present in both the wild and captive groups was 1234, with 898 unique OTUs in the wild group, and 3665 unique OTUs in the captive group (Fig. 1a). The rarefaction curves (Fig. S1A) and rank abundance curves (Fig. S1B) of the WCM and CCM fecal samples convey the richness and evenness of the microbial species in the samples, and the rarefaction curves had reached a plateau. These findings demonstrated that the sequencing data was reasonable and reflected the number of samples in this study was sufficient.
The gut microbiota alpha diversity indexes (Shannon and Chao1) in WCM were significantly higher than that of the CCM at the microbial species level (Kruskal-Wallis test: Effect sizes = 0.137, 0.358; P = 0.006, 7.63E-06. Wilcoxon test: W = 484, 580; P = 0.006, 1.299E-06). The discrepancy between the WCM and CCM was further analyzed via the partial least squares discrimination analysis (PLS-DA) plot as illustrated in Fig. 1b. The distance between the orange and blue color markers demonstrated a unique bacterial community structure in WCM vs. CCM, suggesting that there was a significant difference in gut microbial composition between the two groups.
The predominant microbial phyla identified included Proteobacteria (49.6%), Firmicutes (23.8%), Actinobacteria (8.7%) in the wild group, and Firmicutes (53.5%), Proteobacteria (32.0%), Cyanobacteria (6.0%) in the captive group. At the genus level, Ochrobactrum (12.1%) dominated the WCM gut microbiota, followed by Faecalitalea (9.0%), and Acinetobacter (6.7%), whereas, for captive individuals, Escherichia-Shigella (16.7%), Enterococcus (16.0%) and Streptococcus (6.3%) were primarily identified. The top 20 species for each group at phylum and genus levels were shown in Fig. 2a and b, respectively.
Metabolites in Chinese monal fecal samples were acquired in both positive (POS) and negative (NEG) ion modes. A total of 41,943 mass spectral features were detected, of these, 146 and 177 metabolites were positively identified in POS and NEG modes, respectively, from a library of 2500 known biochemical compounds. PLS-DA score plots showed a distinctive difference in the metabolic profile between the two groups (Fig. 3a & b). In the POS and NEG modes, 58 significant metabolites were identified whose VIP-value>1 and P-value ≤0.05. The heatmap of Hierarchical Clustering was shown in Fig. 4a, b. These metabolites were mainly grouped into fatty acids, bile acids derivatives, sugars, and indole derivatives. These metabolites were then submitted to the Kyoto Encyclopedia of Genes and Genomes (KEGG) website for relevant pathway analysis, the results of which showed that these metabolites were associated with 20 significant enrichment KEGG pathways. The top 10 most significant enrichment results are shown in Fig. 3c. These pathways were primarily related to galactose metabolism, starch and sucrose metabolism, metabolic pathways, biosynthesis of fatty acid, biosynthesis of bile acid, and bile secretion.
Strong associations between microbiota composition (at genus level) and metabolite perturbations were revealed by a Mantel test (r = 0.825, P = 0.001). Significant correlations could be detected between the fecal microbiota and metabolite based on Spearman correlation coefficients (r > 0.5 or r < − 0.5, P < 0.05) . Among the metabolites, fatty acids, bile acid derivatives, and carbohydrate were involved in more than half of the metabolic pathways. Metabolites that were remarkably different were highly correlated with microorganisms. The correlations between specific fecal metabolites and the top 10 most common fecal bacteria genera were shown in Fig. 5a.
Ochrobactrum, Faecalitalea, Acinetobacter, Pelomonas, Ralstonia were the dominant genera of the wild group, represented by F1, Escherichia-Shigella, Enterococcus, Peptoclostridium, Streptococcus and Clostridium_sensu_stricto_1 were dominant in the captive individuals, and Escherichia-Shigella, Enterococcus and Peptoclostridium were represented by F2. The results showed fatty acids (linoleic acid, oleic acid, palmitic acid, stearic acid, 2E-eicosenoic acid), and bile salt (chenodeoxycholate) were higher in the feces of the CCM and negatively correlated with F1, but mainly positively correlated with F2. However, carbohydrate (stachyose, galactinol, sucrose, and raffinose), bile acid (taurodeoxycholic acid, tauroursodeoxycholic acid) and bile salts (taurochenodeoxycholate, taurocholate) were higher in the feces of WCM and mainly negatively correlated with F2, but positively correlated with F1. Remarkably, 5-hydroxyindoleacetic acid and indole-2-carboxylic acid, which increased 113.3 and 295.7-fold, respectively, in CCM, negatively correlated with F1 but positively correlated with F2. Moreover, Sinigrin was decreased 3189.2-fold in the CCM, and positively correlated with F1 but negatively correlated with F2. There was no significant correlation between cellobiose and F1, F2 (Fig. 5a), while cellobiose was positively correlated with Enterorhabdus and negatively correlated with Arthrobacter, Pseudomonas and Paenibacillus (Fig. 5b).
Through fecal analysis, we studied the gut microbiota structure and the metabolomic profiles in both WCM and CCM. The results of this research indicated distinct differences of the fecal bacterial communities between the two groups, furthermore, we identified significant correlations between the fecal microbiota and metabolites.
The levels of alpha and beta-diversity of microbiota from WCM were significantly higher than CCM, indicating a more diverse microbial community in the wild population. These differences may be the result of differences in diet composition, geographical ranges, energy utilization, climate conditions, and stress exposure in the WCM vs. the CCM [10, 22].
Proteobacteria was the dominant phylum in the feces of WCM followed by Firmicutes, while Firmicutes was relatively more abundant than Proteobacteria in the feces of the captive group. The characterization of gut bacteria in several studies reveals a dominance of Firmicutes and Bacteroidetes, which appears to be a common feature of birds. The dominance of these phyla has been described in the ceca of turkeys (Meleagris gallopavo) , the cloacal samples of Chinstrap penguins (Pygoscelis antarctica) , the feces of Japanese quail (Coturnix japonica) , as well as in the crop of hoatzin (Opisthocomus hoazin) . However, previous studies also reported that Proteobacteria and Firmicutes are the dominant phyla in the feces of macaroni penguins (Eudyptes chrysolophus), little blue penguins (Eudyptula minor) , Procellariiform seabirds , Lady Amherst’s pheasant (Chrysolophus amherstiae), Reeves’s pheasant (Syrmaticus reevesii), and Cabot’s tragopan (Tragopan caboti) , which is more consistent with our observations of the microbiota composition in the Chinese monal. Moreover, Bacteroidetes just accounted for 3.8 and 1.6%, respectively, in the WCM and CCM unlike what was observed in the turkey, Chinstrap penguin, quail, and hoatzin studies. Based on the results of this study, it would appear that the gut microbial composition of the Chinese monal is most reflective of that which has been detected in other species of the family Phasianidae; however, in the aforementioned cited research only two samples per group were collected and study details are lacking which highlights the need for further research on this subject.
A large diversity of gut microbes has been described as adaptive and beneficial, with diet being considered one of the most critical factors shaping gut microbial structures. The Chinese monal is an omnivorous bird. For wild populations, food is scarce in the winter months, due to their natural environment being at a high altitude and low temperatures. Whereas, the small captive population of monals receives a relatively steady diet with no change in environment except for exposure to outdoor seasonal climate conditions. Our results indicated a high abundance of Proteobacteria in WCM. Proteobacteria have greatly variable morphology and versatile functions , previous studies have demonstrated that an increased richness of Proteobacteria in the gut microbial flora is mainly related to energy accumulation [31,32,33], and is more abundant when the host animal has prolonged exposure to cold climates . We suspect our findings in WCM may be in response to their comparatively complex dietary composition, as well as an increase in energy storage in preparation for cold weather and a reduction in food availability.
In contrast to WCM, Firmicutes was identified as the dominant bacteria phylum in CCM. Compared with the WCM, we suspect captive individuals are fed a diet with a comparatively higher lipid content. Diet composition, as indicated in Table 1 showed a 6.2% crude fat in the commercial pelleted feed, which provides the base of the diet, in addition to corn which contains 3.6 to 5.3% crude fat (The data provided at the China Feed-database Information Network Centre: http://www.chinafeeddata.org.cn/picture/pdf/CFIC2019_1.pdf). In a 2014 study, mice were fed a high-fat diet, and the Firmicutes: Bacteroidetes ratio in mice gut microbiota increased after 3 weeks, with the abundance of Firmicutes increased and Bacteroidetes depleted . This finding is in accordance with many former studies, suggesting that a high-fat content in diet is one of the primary factors responsible for changes observed in gut microbe composition [7, 35]. Thus, the increased abundances of Firmicutes and the Firmicutes: Bacteroidetes ratio in CCW group (CCM 33.45 versus WCM 6.25) is likely the result of a higher lipid content in the diet.
Concerning the metabolome, we detected 58 significant metabolites that associated with 20 metabolic pathways. We found that a significantly higher fatty-acid content in feces of the CCM was related to unsaturated fatty acid and fatty acid biosynthesis. A previous study found that different dietary triacylglyceride composition has an influence on fatty acid content in the feces of human neonates . Thus, our finding may be related to the nutrient composition of the diet, which contains a higher level of fat and a different lipid composition than what wild specimens are likely consuming. Otherwise, our results showed that the primary microbial genera in Chinese monal feces are significantly related to the metabolites we identified; therefore, microbial metabolism is an essential aspect that needs to be evaluated. Microbial biosynthesis of fatty acid may be due to cellular structure or storage, de novo synthesis from glucose, or incorporation of exogenous fatty acids directly into lipid structures [37, 38].
The content of metabolites in the feces relating to bile acid biosynthesis and bile secretion had a statistically significant difference between the WCM and CCM, including chenodeoxycholate, taurochenodeoxycholate, taurocholate, taurodeoxycholic acid, and tauroursodeoxycholic acid. Our results showed differences in metabolism of bile acid, an important component of bile, which plays a role in digestion and metabolism of dietary lipids and cholesterol [39, 40]. Chenodeoxycholate, taurochenodeoxycholate, and taurocholate are bile salts that are hydrolyzed into free bile acids by bile salt hydrolase . Chenodeoxycholate amount was higher in the feces of CCM, whereas taurochenodeoxycholate and taurocholate were more abundant in the WCM. Taurodeoxycholic acid and tauroursodeoxycholic acid are taurine conjugated bile acids that were expectedly more abundant in the feces of WCM. The type of bile salt that is most abundant may reflect how the corresponding bile acids play a major role in digestion which is dependent on the dietary composition of each individual bird . In the wild, Chinese monals primarily consume seeds and roots of shrubs in the winter when food is scarce, but they also consume a small amount of moss, earthworms and insect pupa [43, 44]. High concentrations of total non-structural carbohydrate (TNC) reserves are usually found in root tissues of plants [45,46,47], and may serves as one of the main energy sources for WCM. Because of the seed consumption, the WCM also need to secrete bile acids for fat digestion thus providing another valuable energy source. We also found a high correlation between the bile acid metabolites and gut microbial flora. There exists a mutual regulation between intestinal microbes and bile acids, for bile acid composition seems to be a critical regulator of microbiota structure, which in turn plays a key role in regulation of bile acid pool size, composition and metabolism [48, 49]. Kim G et al.  found that several bacterial species of Firmicutes and Actinobacteria can hydrolyze bile salt, among which Enterococcus was significantly related to the above metabolites in this study. A significantly higher abundance of carbohydrates was found in the feces of WCM as opposed to the CCM, including galactinol, raffinose, stachyose, sucrose, and cellobiose. In plants, raffinose and stachyose are oligosaccharides synthesized from sucrose, subsequently galactitol and galactose are also involved . Animals don’t secrete enzymes that utilize raffinose-series oligosaccharides, which are likely to be digested by microbial enzymes at the end of the gastrointestinal tract . Coon et al.  reported that ileal digestibility of raffinose and stachyose was less than 1.0% in roosters; however, the digestibility, determined with excreta collection, reached 90.5 and 83.8%, respectively. Previous studies have also indicated several Passerine bird species have a low digestibility of sucrose [53, 54] as fecal sugar content increased after consuming sucrose solutions . Cellobiose is an important hydrolytic product of cellulose degradation . Plants contain significant quantities of cellulose, which is difficult to naturally digest by most birds species, but can be utilized by bacteria inhabiting the intestinal tract, especially within the ceca . Yuzhang Wang  reported that the cecum of the Chinese monal is well developed, accounting for 24.09% of the total length of the intestine, which may help the species utilize cellulose and other indigestible components of the diet. Microorganisms in genus level that have significant positive correlation with these metabolites may play a big role in these processes.
Indole-containing metabolites were significantly different in fecal samples of CCM, such as 5-hydroxyindoleacetic acid (5-HIAA) and indole-2-carboxylic acid, which are nearly 100 and 300 times higher than that in the wild, respectively. Indole and its derivatives widely exist in animals, microorganisms, and plants , and are the main metabolites of serotonin . Degg et al.  found that excretion of indole containing metabolites, such as 5-HIAA, increased after humans ingested high serotonin-containing foods including tomatoes, bananas, etc. Our results which indicate higher levels of 5-HIAA and indole-2-carboxylic acid in CCM feces is likely associated with the incorporation of tomatoes in their routine diet. These two altered indole-containing metabolites were also highly correlated with fecal microbial flora. Previous studies have suggested that intestinal bacteria can convert tryptophan to indole through enzymatic processing, and then form indole-containing metabolites [40, 61]. These studies showed that metabolic processes in intestinal bacteria are essential for the synthesis of indole-containing metabolites, therefore, it may be a specific indicator of intestinal bacterial differences.
Finally, the last metabolite with a significant difference observed between the two groups is sinigrin (2-propenyl glucosinolate). Sinigrin was several thousand times higher in the feces of WCM than that observed in CCM. Sinigrin, which is enriched in plants such as cruciferous vegetables, leaf mustard, and horseradish, has been associated with carbohydrate regulation and lipid metabolism and has also been shown to have anti-neoplasia and anti-microbial properties . In vitro experiments confirmed that sinigrin was degraded by rat intestinal microbiota . Enterococcus and Bacteroides in the human intestine can degrade glucosinolates , and have a significant positive correlation with sinigrin in this study. Therefore, based on our findings of high levels of sinigrin in WCM feces, we can speculate that the diet of WCM also contains a large amount of sinigrin. Although the intestinal flora plays a role in the digestion process, sinigrin is more excreted through feces.
This is the first study to evaluate the fecal microbiota and the fecal metabolic profiles in the Chinese monal. The results clearly demonstrate the differences of the gut bacterial composition and metabolism of the host and microbiota between wild and captive Chinese monals. The results of this study, will without a doubt, enhance our knowledge of wild and captive Chinese monals which may help guide recommendations that will affect the health and fecundity of birds in the captive breeding program. Further studies are required to determine if there is a difference in the fecal microbiota and metabolic profiles in spring and summer months when the dietary composition is likely to change.
Birds selected, study site and sample collection
Samples were collected from October to December of 2018 at Baoxing Fengtongzhai National Nature Reserve in Sichuan Province, China. Fecal samples from the wild were collected at two research sites, located on two mountains: Yuancaopo (30°34′38′′-30°35′3′′N, 102°48′15′′- 102°49′57′′E) and Hongshanding (30°37′39′′- 30°37′46′′N, 102°54′57′′- 102°55′3′′E). Samples that appeared fresh (e.g. still wet) were collected at an altitude of 3223.7 to 3568.7 m, and a temperature of − 17 °C – 5 °C. In Yuancaopo and Hongshanding, 49 and 42 samples were collected, respectively.
In November 2018, fecal samples of 11 captive birds (3 males and 8 females) were collected at the Chinese Monal Conservation and Research Center, where the altitude is 1582.0 m and the temperature range from 9 °C – 13 °C. These monals are maintained in four adjacent enclosures, with three enclosures containing one male each paired with one or two females, while one enclosure contains only three sub-adult females. Feces was collected by placing clean cardboard beneath the nightly roosting perches, overnight. The following morning at 8 AM, four fecal samples per enclosure were collected from the cardboard with an attempt to collect feces with varying colors and shapes.
Approximately 1 g of fecal material was collected from the central section of each fecal dropping, from both wild and captive animals, and placed into 1.5 mL sterile collection tubes (Corning, New York, USA) and immediately frozen in liquid nitrogen. Thereafter, the samples were shipped to the Chengdu Research Base of Giant panda Breeding within 24 h., where they were stored at − 80 °C for subsequent analysis.
Measurement of the main components of the commercial pellet feed
The commercial pellet (New hope group, Chengdu, China) being fed to the CCM was analyzed to measure the primary nutritional components according to Chinese national standards. Moisture was determined by direct drying method . Crude fat was determined by petroleum ether extraction method . Crude protein was determined by Kjeldahl method . Crude fiber was determined by Soxhlet extraction method .
Monal species identification
Genomic DNA was isolated from the fecal samples using the TIANamp Micro DNA Kit (Tiangen, Beijing, China) following the manufacturer’s instructions. The DNA was amplified using cytochrome b (Cyt b) gene primers; forward: 5′-ACATTGGACGCGGCCTCTAC-3′ and reverse: 5′-GTGGGCGAAATGTTATGGTT-3′. PCR products were separated by electrophoresis in a 1% agarose gel, then recovered using the AxyPrep DNA gel recovery kit (Axygen, California, USA). The purified PCR products were sequenced using an ABI3730-XL genetic sequencer (Applied Biosystems, Foster City, CA, USA). Genetic sequencing data were blasted against the nucleic acid database of the National Center for Biotechnology Information (NCBI) using the Basic Local Alignment Search Tool (BLAST) program. Species identification was selected based on results that returned a BLANK % sequence similarity result.
16S rRNA microbial community analysis
The samples were homogenized separately, and then the total DNA was extracted from fecal samples (200 mg per sample) of all Chinese monal using the DNeasy PowerSoil Kit (QIAGEN, Inc. Netherlands), using their designated protocol. The V4-V5 region of the 16S rRNA gene was PCR amplified using the forward primer 515 F: 5′-GTGCCAGCMGCCGCGGTAA-3′ and the reverse primer 907 R: 5′-CCGTCAATTCMTTTRAGTTT-3′. Sample-specific 7-bp barcodes were incorporated into the primers for multiplex sequencing. Before being pooled, PCR amplicons were purified with Agencourt AMPure Beads (Beckman Coulter, Indianapolis, IN, USA), quantified using the PicoGreen dsDNA Assay Kit (Invitrogen, Carlsbad, CA, USA), and finally, paired-end sequencing with 2 × 300 bp read lengths was performed using the Illumina MiSeq platform with MiSeq Reagent Kit v3 at Shanghai Personal Biotechnology Co., Ltd. (Shanghai, China).
The raw sequencing data were processed using the Quantitative Insights Into Microbial Ecology (QIIME, v1.8.0) software package. Briefly, raw sequencing reads with exact matches to the barcodes were assigned to respective samples and identified as valid sequences. The low-quality sequences were filtered through criteria as previously described: sequences that had a length of < 150 bp, sequences that had average Phred scores of < 20, sequences that contained ambiguous bases, and sequences that contained mononucleotide repeats of > 8 bp . Potential chimreas were removed using the Flash usearch61 tool with the minsnp parameter set to 2, and the remaining high-quality sequences were clustered into OTUs with a threshold of 97% sequence similarity by UCLUST. OTUs containing less than 0.001% of total sequences across all samples were discarded, and OTU taxonomic classification was conducted via BLAST searching the representative sequences set against the Greengenes Database  using the best hit .
Sequence data analyses were primarily calculated using QIIME (v1.8.0) and R packages (v3.2.0). We generated OTU-level ranked abundance curves with R package “ggplot 2” to compare the richness and evenness of OTUs among samples. Venn diagram was obtained to visualize the mutual and unique OTUs between groups using R package “VennDiagram”. Alpha diversity was calculated with QIIME software and discrepancies analyzed between the two groups using the Kruskal-Wallis test and Wilcoxon test with R software. Beta diversity analysis was performed to investigate the structural variation of microbial communities across samples using UniFrac distance metrics  and visualized via PLS-DA using the mixOmics package in R software. The microbiota compositional profiles could be achieved at different taxonomic levels, including phylum, class, order, family and genus, wherewith the taxa abundances at these levels were statistically compared among samples or groups by Metastats . The sequences data available for our study were submitted to NCBI Sequence Read Archive under the BioProject ID PRJNA662165.
Nine samples were selected from WCM and CCM, respectively, for untargeted metabolomics analyses, which was performed using ultra-high performance liquid chromatography (UHPLC) (1290 Infinity LC, Agilent Technologies, Palo Alto, USA) coupled to a Triple TOF 6600 mass spectrometer (AB Sciex, Foster City, USA) equipped with an electrospray ionization (ESI) source in positive and negative ion modes. Samples (500 mg per sample) were pretreated and then separated by UHPLC system, followed by mass spectrometry analysis. Samples were mixed in equal amounts to prepare quality control (QC) samples, which were spaced evenly among the injections. All experimental samples were randomly distributed throughout the run in order to monitor the precision and stability of the method during its operation.
The raw data were converted to mzXML files using ProteoWizard MSConvert and processed using XCMS for peak detection, alignment, and data filtering [73, 74]. Mass accuracy tolerance within 25 ppm was used as the mass window for the database search and secondary spectral map matching, and the standards database of the laboratory was retrieved. For the data extracted by XCMS, ion peaks of group sum > 50% were deleted. Pattern recognition was carried out by applying the software SIMCA-P 14.1 (Umetrics, Umea, Sweden).
PLS-DA was performed for the supervised multivariate statistical analysis using MetaboAnalyst (www.metaboanalyst.ca) web-based system. The significantly different metabolites were determined based on the combination of a statistically significant threshold of variable influence on projection (VIP) values obtained from PLS-DA model and two-tailed Student’s t test (p value) on the raw data, and the metabolites with VIP values larger than 1.0 and p values less than 0.05 were considered as significant. Heat maps were generated using a hierarchical clustering algorithm to visualize the metabolite difference within the data set. The data were normalized using z-scores of the intensity areas of differential metabolites and were ploted by Pheatmap package in R language. The differential metabolites were submitted to the KEGG website (http://www.genome.jp/kegg/pathway.html) for related pathway analysis. The metabolic pathway enrichment of these metabolites was performed by Fisher’s exact test in Omicsbean, when P-value of metabolic pathway < 0.05, metabolic pathway was considered as statistically significant enrichment. The metabolomics data were submitted to MetaboLights (https://www.ebi.ac.uk/metabolights/) under accession number MTBLS2089.
The correlation matrix between the gut microflora–related metabolites and gut bacterial species was generated using Pearson’s correlation coefficient.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Convention on International Trade in Endangered Species of Wild Fauna and Flora
International Union for Conservation of Nature
Wild Chinese monal
Captive Chinese monal
- 16S rRNA:
16S ribosomal RNA
Ultra-high performance liquid chromatography
operational taxonomic unit
Partial least squares discrimination analysis
National Center for Biotechnology Information
Basic Local Alignment Search Tool
Kyoto Encyclopedia of Genes and Genomes
Wang B, Xu Y, Ran J. Predicting suitable habitat of the Chinese monal (Lophophorus lhuysii) using ecological niche modeling in the Qionglai Mountains, China. PeerJ. 2017;5:e3477.
CITES. Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES) Appendices I, II and III. Available at http://www.cites.org. Accessed 17 Mar 2020.
BirdLife International. Lophophorus lhuysii. The IUCN Red List of Threatened Species 2016. Available at https://doi.org/10.2305/IUCN.UK.2016-3.RLTS.T22679192A92806697.en. Accessed 17 Mar 2020.
Zheng G. Pheasants in China: the Chinese monal (Lophophorus lhuysii). In: Zheng G, Ding P, Lu X, Zhang Z, Zhang Y, editors. Higher education press, Beijing; 2015. p.338–351.
Martin F-PJ, Collino S, Rezzi S, Kochhar S. Metabolomic applications to decipher gut microbial metabolic influence in health and disease. Front Physio. 2012;3:113.
He F, Zhai J, Zhang L, Liu D, Ma Y, Rong K, et al. Variations in gut microbiota and fecal metabolic phenotype associated with Fenbendazole and Ivermectin tablets by 16S rRNA gene sequencing and LC/MS-based metabolomics in Amur tiger. Biochem Biophys Res Commun. 2018;499:447–53.
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:1027–31.
Wang J, Tang H, Zhang C, Zhao Y, Derrien M, Rocher E, et al. Modulation of gut microbiota during probiotic-mediated attenuation of metabolic syndrome in high fat diet-fed mice. ISME J. 2015;9:1–15.
Guan Y, Yang H, Han S, Feng L, Wang T, Ge J. Comparison of the gut microbiota composition between wild and captive sika deer (Cervus nippon hortulorum) from feces by high-throughput sequencing. AMB Expr. 2017;7:212.
Schmidt E, Mykytczuk N, Schulte-Hostedde AI. Effects of the captive and wild environment on diversity of the gut microbiome of deer mice (Peromyscus maniculatus). ISME J. 2019;13:1293–305.
Ushida K, Segawa T, Tsuchida S, Murata K. Cecal bacterial communities in wild Japanese rock ptarmigans and captive Svalbard rock ptarmigans. J Vet Med Sci. 2016;78:251–7.
Scupham AJ, Patton TG, Bent E, Bayles DO. Comparison of the Cecal microbiota of domestic and wild turkeys. Microb Ecol. 2008;56:322–31.
Xenoulis PG, Gray PL, Brightsmith D, Palculict B, Hoppes S, Steiner JM, et al. Molecular characterization of the cloacal microbiota of wild and captive parrots. Vet Microbiol. 2010;146:320–5.
Salgado-Flores A, Tveit AT, Wright A-D, Pope PB, Sundset MA. Characterization of the cecum microbiome from wild and captive rock ptarmigans indigenous to Arctic Norway. Chaves AV, editor. PLoS ONE. 2019;14:e0213503.
Uenishi G, Fujita S, Ohashi G, Kato A, Yamauchi S, Matsuzawa T, et al. Molecular analyses of the intestinal microbiota of chimpanzees in the wild and in captivity. Am J Primatol. 2007;69:367–76.
Gao H, Chi X, Qin W, Wang L, Song P, Cai Z, et al. Comparison of the gut microbiota composition between the wild and captive Tibetan wild ass (Equus kiang). J Appl Microbiol. 2019;126:1869–78.
Marcobal A, Kashyap PC, Nelson TA, Aronov PA, Donia MS, Spormann A, et al. A metabolomic view of how the human gut microbiota impacts the host metabolome using humanized and gnotobiotic mice. ISME J. 2013;7:1933–43.
Xie G, Zhang S, Zheng X, Jia W. Metabolomics approaches for characterizing metabolic interactions between host and its commensal microbes. Electrophoresis. 2013;19:2787–98.
Kieffer DA, Piccolo BD, Vaziri ND, Liu S, Lau WL, Khazaeli M, et al. Resistant starch alters gut microbiome and metabolomic profiles concurrent with amelioration of chronic kidney disease in rats. Am J Physiol Renal Physiol. 2016;310:F857–71.
Garcia-Mazcorro JF, Pedreschi R, Yuan J, Kawas JR, Chew B, Dowd SE, et al. Apple consumption is associated with a distinctive microbiota, proteomics and metabolomics profile in the gut of Dawley Sprague rats fed a high-fat diet. Blachier F, editor. PLoS ONE. 2019;14:e0212586.
Gomez A, Petrzelkova K, Yeoman CJ, Vlckova K, Mrázek J, Koppova I, et al. Gut microbiome composition and metabolomic profiles of wild western lowland gorillas (Gorilla gorilla gorilla) reflect host ecology. Mol Ecol. 2015;24:2551–65.
Nelson TM, Rogers TL, Carlini AR, Brown MV. Diet and phylogeny shape the gut microbiota of Antarctic seals: a comparison of wild and captive animals: the gut microbiota of wild and captive Antarctic seals. Environ Microbiol. 2013;15:1132–45.
D’Andreano S, Sànchez Bonastre A, Francino O, Cuscó Martí A, Lecchi C, Grilli G, et al. Gastrointestinal microbial population of Turkey (Meleagris gallopavo) affected by hemorrhagic enteritis virus. Poult Sci. 2017;96:3550–8.
Barbosa A, Balagué V, Valera F, Martínez A, Benzal J, Motas M, et al. Age-Related Differences in the Gastrointestinal Microbiota of Chinstrap Penguins (Pygoscelis antarctica). Peter H-U, editor. PLoS ONE. 2016;11:e0153215.
Wilkinson N, Hughes RJ, Aspden WJ, Chapman J, Moore RJ, Stanley D. The gastrointestinal tract microbiota of the Japanese quail, Coturnix japonica. Appl Microbiol Biotechnol. 2016;100:4201–9.
Godoy-Vitorino F, Ley RE, Gao Z, Pei Z, Ortiz-Zuazaga H, Pericchi LR, et al. Bacterial Community in the Crop of the hoatzin, a Neotropical Folivorous flying bird. Appl Environ Microbiol. 2008;74:5905–12.
Dewar ML, Arnould JPY, Dann P, Trathan P, Groscolas R, Smith S. Interspecific variations in the gastrointestinal microbiota in penguins. MicrobiologyOpen. 2013;2:195–204.
Dewar ML, Arnould JPY, Krause L, Dann P, Smith SC. Interspecific variations in the faecal microbiota of Procellariiform seabirds. FEMS Microbiol Ecol. 2014;89:47–55.
Shen J, Zhang X, Wu W, Hu W, Zhang C, Zhu L. The Gut Microbes’Diversity and Function of Phasianidae Provide Insights on the Adaptation to Their Diet. J Nanjing Normal Univ (Natural Science Edition). 2016;39:90–5.
Shin N-R, Whon TW, Bae J-W. Proteobacteria: microbial signature of dysbiosis in gut microbiota. Trends Biotechnol. 2015;33:496–503.
Koren O, Goodrich JK, Cullender TC, Spor A, Laitinen K, Kling Bäckhed H, et al. Host remodeling of the gut microbiome and metabolic changes during pregnancy. Cell. 2012;150:470–80.
Chevalier C, Stojanović O, Colin DJ, Suarez-Zamorano N, Tarallo V, Veyrat-Durebex C, et al. Gut microbiota orchestrates energy homeostasis during cold. Cell. 2015;163:1360–74.
Sun B, Wang X, Bernstein S, Huffman MA, Xia D-P, Gu Z, et al. Marked variation between winter and spring gut microbiota in free-ranging Tibetan macaques (Macaca thibetana). Sci Rep. 2016;6:26035.
Walker A, Pfitzner B, Neschen S, Kahle M, Harir M, Lucio M, et al. Distinct signatures of host–microbial meta-metabolome and gut microbiome in two C57BL/6 strains under high-fat diet. ISME J. 2014;8:2380–96.
Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457:480–4.
López-López A, Castellote-Bargalló AI, Campoy-Folgoso C, Rivero-Urgel M, Tormo-Carnicé R, Infante-Pina D, et al. The influence of dietary palmitic acid triacylglyceride position on the fatty acid, calcium and magnesium contents of at term newborn faeces. Early Hum Dev. 2001;65:S83–94.
Certik M, Shimizu S. Biosynthesis and regulation of microbial polyunsaturated fatty acid production. J Biosci Bioeng. 1999;87:1–14.
Steen EJ, Kang Y, Bokinsky G, Hu Z, Schirmer A, McClure A, et al. Microbial production of fatty-acid-derived fuels and chemicals from plant biomass. Nature. 2010;463:559–62.
Nguyen A, Bouscarel B. Bile acids and signal transduction: role in glucose homeostasis. Cell Signal. 2008;20:2180–97.
Lu K, Abo RP, Schlieper KA, Graffam ME, Levine S, Wishnok JS, et al. Arsenic exposure perturbs the gut microbiome and its metabolic profile in mice: an integrated Metagenomics and metabolomics analysis. Environ Health Perspect. 2014;122:284–91.
Kim G-B, Lee BH. Biochemical and molecular insights into bile salt hydrolase in the gastrointestinal microflora - a review. Asian Australas J Anim Sci. 2005;18:1505–12.
Yeh Y-H, Hwang D-F. High-performance liquid chromatographic determination for bile components in fish, chicken and duck. J Chromatogr B Biomed Sci Appl. 2001;751:1–8.
He F, Lu T. Ecology of the Chinese monal in winter. Zool Res. 1985;6:345–52.
Long Y, Shao K, Guo G, Cheng C, Zhou X, Hans L, et al. A follow-up study of the ecology of the Chinese monal in winter. Sichuan J Zool. 1998;17:104–5.
Loescher WH, McCamant T, Keller JD. Carbohydrate reserves, translocation, and storage in Woody Plant roots. HortSci. 1990;25:274–81.
Cruz A, Moreno JM. Seasonal course of Total non-structural carbohydrates in the Lignotuberous Mediterranean-type shrub Erica australis. Oecologia. 2001;128:343–50.
Knox KJE, Clarke PJ. Nutrient availability induces contrasting allocation and starch formation in resprouting and obligate seeding shrubs. Funct Ecol. 2005;19:690–8.
Ridlon JM, Alves JM, Hylemon PB, Bajaj JS. Cirrhosis, bile acids and gut microbiota: unraveling a complex relationship. Gut Microbes. 2013;4:382–7.
Gérard P. Metabolism of cholesterol and bile acids by the gut microbiota. Pathogens. 2013;3:14–24.
Sengupta S, Mukherjee S, Basak P, Majumder AL. Significance of galactinol and raffinose family oligosaccharide synthesis in plants. Front Plant Sci [Internet]. 2015 [cited 2020 May 12];6. Available from: http://journal.frontiersin.org/Article/10.3389/fpls.2015.00656/abstract.
Iji PA, Saki AA, Tivey DR. Intestinal structure and function of broiler chickens on diets supplemented with a mannan oligosaccharide. J Sci Food Agric. 2001;81:1186–92.
Coon CN, Leske KL, Akavanichan O, Cheng TK. Effect of oligosaccharide-free soybean meal on true Metabolizable energy and Fiber digestion in adult roosters. Poult Sci. 1990;69:787–93.
Martínez del Rio C, Baker HG, Baker I. Ecological and evolutionary implications of digestive processes: Bird preferences and the sugar constituents of floral nectar and fruit pulp. Experientia. 1992;48:544–51.
Apanius V. Sucrose intolerance in birds: simple nonlethal diagnostic methods and consequences for assimilation of complex carbohydrates. Auk. 1994;111:170–7.
Lila ZA, Mohammed N, Takahashi T, Tabata M, Yasui T, Kurihara M, et al. Increase of ruminal fiber digestion by cellobiose and a twin strain of Saccharomyces cerevisiae live cells in vitro. Animal Sci J. 2006;77:407–13.
Inman DL. Cellulose digestion in ruffed grouse, Chukar partridge, and bobwhite quail. J Wildl Manag. 1973;37:114–21.
Wang Y, Yue X, Hu J, Guo Y, Wu X. Morphology of digestive system of Lophophorus lhuysii. J China West Normal University (Natural Science). 2007;28:7–10.
Lu M, Hashimoto K, Uda Y. Rat intestinal microbiota digest desulfosinigrin to form allyl cyanide and 1-cyano-2,3-epithiopropane. Food Res Int. 2011;44:1023–8.
Mokhtari M, Rezaei A, Ghasemi A. Determination of urinary 5-Hydroxyindoleacetic acid as a metabolomics in gastric cancer. J Gastrointest Canc. 2015;46:138–42.
Degg TJ, Allen KR, Barth JH. Measurement of plasma 5-hydroxyindoleacetic acid in carcinoid disease: an alternative to 24-h urine collections? Ann Clin Biochem. 2000;37:724–6.
Sheflin A, Borresen E, Wdowik M, Rao S, Brown R, Heuberger A, et al. Pilot dietary intervention with heat-stabilized Rice bran modulates stool microbiota and metabolites in healthy adults. Nutrients. 2015;7:1282–300.
Patel D, Patel K, Gadewar M, Tahilyani V. A concise report on pharmacological and bioanalytical aspect of sinigrin. Asian Pac J Trop Biomed. 2012;2:S446–8.
Krul C, Humblot C, Philippe C, Vermeulen M, van Nuenen M, Havenaar R, et al. Metabolism of sinigrin (2-propenyl glucosinolate) by the human colonic microflora in a dynamic in vitro large-intestinal model. Carcinogenesis. 2002;23:1009–16.
GB/T 6435–2014. Determination of moisture in feedstuffs. Standards Press of China, Beijing; 2014.
GB/T 6433–2006. Determination of crude fat in Feeds. Standards Press of China, Beijing, China; 2006.
GB/T 6432–1994. Method for the determination of crude protein in Feedstuffs. Standards Press of China, Beijing; 1994.
GB/T 6434–2006. Determination of crude fiber in Feeds. Beijing: Standards Press of China; 2006.
Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27:2957–63.
DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. AEM. 2006;72:5069–72.
Altschul S. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.
Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. AEM. 2005;71:8228–35.
White JR, Nagarajan N, Pop M. Statistical Methods for Detecting Differentially Abundant Features in Clinical Metagenomic Samples. Ouzounis CA, editor. PLoS Comput Biol. 2009;5:e1000352.
Ivanisevic J, Epstein AA, Kurczy ME, Benton PH, Uritboonthai W, Fox HS, et al. Brain region mapping using global metabolomics. Chem Biol. 2014;21:1575–84.
Benton HP, Ivanisevic J, Mahieu NG, Kurczy ME, Johnson CH, Franco L, et al. Autonomous metabolomics for rapid metabolite identification in global profiling. Anal Chem. 2015;87:884–91.
We thank the Shanghai Personal Biotechnology Co., Ltd. for the bioinformatics analysis support provided.
This work was supported by the Chengdu Giant Panda Breeding Research Foundation (CPF2017–10); and the independent project of Chengdu Research Base of Giant Panda Breeding (2020CPB-C10).
Ethics approval and consent to participate
The experimental protocol was approved by the Chengdu Research Base of Giant Panda Breeding and in accordance with the current laws on animal welfare and research in China.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Jiang, D., He, X., Valitutto, M. et al. Gut microbiota composition and metabolomic profiles of wild and captive Chinese monals (Lophophorus lhuysii). Front Zool 17, 36 (2020). https://doi.org/10.1186/s12983-020-00381-x