key: cord-0019041-tgmh5nju
authors: Sawaswong, Vorthon; Praianantathavorn, Kesmanee; Chanchaem, Prangwalai; Khamwut, Ariya; Kemthong, Taratorn; Hamada, Yuzuru; Malaivijitnond, Suchinda; Payungporn, Sunchai
title: Comparative analysis of oral-gut microbiota between captive and wild long-tailed macaque in Thailand
date: 2021-07-12
journal: Sci Rep
DOI: 10.1038/s41598-021-93779-4
sha: ef7aba5759ae28a4e521ee21c823d63a12c3e21a
doc_id: 19041
cord_uid: tgmh5nju

Long-tailed macaques (Macaca fascicularis), distributed in Southeast Asia, are generally used in biomedical research. At present, the expansion of human communities overlapping of macaques’ natural habitat causes human-macaque conflicts. To mitigate this problem in Thailand, the National Primate Research Center of Thailand, Chulalongkorn University (NPRCT-CU), was granted the permit to catch the surplus wild-born macaques and transfer them to the center. Based on the fact that the diets provided and the captive environments were different, their oral-gut microbiota should be altered. Thus, we investigated and compared the oral and fecal microbiome between wild-born macaques that lived in the natural habitats and those transferred to and reared in the NPRCT-CU for 1 year. The results from 16S rRNA high-throughput sequencing showed that the captive macaques had distinct oral-gut microbiota profiles and lower bacterial richness compared to those in wild macaques. The gut of wild macaques was dominated by Firmicutes which is probably associated with lipid absorption and storage. These results implicated the effects of captivity conditions on the microbiome that might contribute to crucial metabolic functions. Our study should be applied to the animal health care program, with respect to microbial functions, for non-human primates.

www.nature.com/scientificreports/ for humans due to their similar genetic and physiological characteristics 15 . In addition, the NHPs have typically been used in biomedical research, particularly in pharmaceutical development 16, 17 . Recent studies have suggested that the microbiota also contributes in the modulation of vaccine and drug responses 18, 19 . Accordingly, to utilize the NHPs either in microbiome study or drug and vaccine testing, the baseline microbial profiles of these animals should first be explored. The long-tailed macaque (Macaca fascicularis), also known as cynomolgus macaque, is one of the NHPs commonly used in biomedical research 18, 20 , notably in the recent COVID-19 vaccine research 21 . It is widely distributed in Southeast Asian countries 22 and can be found in all regions of Thailand. Regarding their overpopulation status which has caused the human-macaque conflict 23 , the National Primate Research Center of Thailand, Chulalongkorn University (NPRCT-CU) was granted the permit from the Department of the National Parks, Wildlife and Plant Conservation to catch the wild-borne macaques and use them as breeding founders of the center. These macaques were reared under a well-controlled and hygienic condition and the NPRCT-CU has been awarded the AAALAC International Accreditation since February 2020. Nevertheless, the captivity restricted the roaming and foraging behaviors of captive macaques. In addition, rearing them in the cage may enhance the close contact among macaques and allow the transmission of microbes via oral-fecal route which could affect their microbiome. The shift of diets and the captive living environment could also cause the alteration in their microbiota, as previously observed 24 . A recent study suggested that the microbiota in wild animals could enhance the host response in inflammatory stimuli and improve disease resistance 25 . Concordantly, several primate centers have an incidence of idiopathic diarrhea which was possibly caused by dysbiosis in microbiota 26 . Thus, the studies of microbiome in captive (wild-born) macaques living in NPRCT-CU might elucidate the captivity effects on microbiota composition and their roles in influencing the health of animals. This study aimed to explore the oral and fecal microbiome of the captive macaques compared to those in the wild. The knowledge of this study should be beneficial for translation to animal healthcare and management in the primate centers.

Characteristics of macaque and sequencing summary. The characteristics of long-tailed macaques were shown in Table 1 . The sex ratio between the two populations had no difference, but the difference was detected for the age which was due to the uncontrollable limitation of subject recruitment. The wild group had significantly older age (P = 4.0e−3) and heavier bodyweight (P = 0.01) than the captive macaques (tested by t-test, P < 0.05). The age, sex and body mass sightly influent the beta diversity (Jaccard dissimilarity) of microbiome as showed in Supplementary Fig. S1 , but these effects were less than the effect of living condition (wild versus captive). Thus, the comparison of microbiome between wild and captive macaques was the main focus in this study.

The number of samples and the sequencing output summary per group including raw reads, filtered reads and identified reads were described in Table 2 . The samples which were unable to amplify for 16S library preparation were excluded from the study. Over 99% of cleaned reads were successfully classified by the QIIME2 pipeline 27 . The rarefaction analysis was also examined to check the sampling coverage by sequencing depth. The results indicated that the bacterial classification was rarified in all study groups as shown in Fig. 1a. Alpha diversity of oral and fecal bacteria. Bacterial diversity was evaluated based on Shannon's diversity index (Fig. 1b) while the richness of bacterial operational taxonomic units (OTUs) was determined by www.nature.com/scientificreports/ Chao1 index (Fig. 1c) . The statistical comparisons of indices between groups were carried out by Mann-Whitney U test with P < 0.05. The Shannon's diversity of oral microbiome in captive (CO) macaques (3.15 ± 0.60) exhibited significantly higher (P = 0.03) compared to those in the wild (WO) group (2.88 ± 0.75). However, the Chao1 index indicated the greater OTUs richness (P = 8.63e−4) in WO (189 ± 31.81) than the CO (160. 8 ± 30.49) . The Shannon's diversity of fecal microbiome in wild (WF) macaques (3.68 ± 0.46) was significantly higher (P = 2.41e−8) than the fecal microbiome of the captive (CF) group (2.96 ± 0.36) whereas the Chao1 richness was not significantly different (P = 0.92) between WF (137.2 ± 34.1) and CF (137.5 ± 34.7) group. Between oral and fecal microbiota, the Shannon diversity of CO was slightly higher (P = 5.64e−3) than CF samples but it was significantly lower (P = 2.61e−5) in WO compared to WF. We found significantly higher Chao1 richness in the oral microbiome than the fecal microbiome in both wild group (P = 1.15e−6) and captive groups (P = 1.03e−4).

To investigate the shared OTUs between groups, the shared OTUs were counted and illustrated in the Venn Euler diagram shown in Fig. 1d . The results showed that about 76 OTUs were shared across all groups suggesting that they constitute the common bacterial populations in the oral-gut environment. Interestingly, 192 sharing OTUs between WO and CO were found, whereas fewer overlapping taxa of WF and CF were noted (comprising only 144 common OTUs). The WO microbiome contained the highest number of unique taxa compared to other groups.

Beta diversity of microbiota in long-tailed macaque. The beta diversity was investigated based on Jaccard dissimilarity index and presented by principal coordinate analysis (PCoA) plot in Fig. 1e . The distance of centroids and dispersion of the groups were statistically evaluated by permutational multivariate analysis of variance (PERMANOVA) with P < 0.05. The findings indicated that the bacterial profiles between wild and captive macaques were significantly different (P = 9.99e−4). It also showed a very high dissimilarity (P = 2.00e−3) between the oral and fecal microbiome. The Fig. 1f presented the comparisons of average distance between groups and within-group which were statistically tested by Mann-Whitney U test, P < 0.05. The result presented significantly higher community distances (P = 8.76e−37) in WO group compared to those in CO group while they were also higher (P = 5.48e−87) in WF than those in CF microbiome. Rarefaction curve indicated the observed OTUs (mean ± SD) against sequencing depth in each group: CO captive-oral, WO wild-oral, CF captive-fecal, WF wild-fecal. The alpha diversity was determined by two indices: (b) Shannon's diversity index and (c) Chao1 index which were plotted by mean ± SD and statistically tested by Mann-Whitney U test (*P < 0.05, **P < 0.01, ***P < 0.001, **** P < 0.0001, ns not significant). (d) The Venn diagram showed the numbers of shared taxa among groups. (e) The beta diversity was plotted in principal coordinate analysis (PCoA) plots based on Jaccard dissimilarity (tested by PERMANOVA analysis). (f) The box plot represented the Jaccard distances between communities which were statistically tested by Mann-Whitney U test. (Fig. 2a) . The second most enriched phylum was Proteobacteria which accounted for 38.7 ± 14.0% in CO and 24.6 ± 13.0% in the WO group. The Bacteroidetes was the third most abundant phylum in CO (14.0 ± 8.4%), while Cyanobacteria was the third most enriched bacteria in WO (18.9 ± 25.5%). It was observed that some of WO have a distinctively high abundance of Cyanobacteria. The remaining community included less abundant phyla such as Fusobacteria, Actinobacteria, Spirochaetes, and Tenericutes. The comparisons of phylum tested by Mann-Whitney U test (P < 0.05) showed that relative abundances of Bacteroidetes and Proteobacteria were significantly greater in the CO than in the WO (P = 1.51e−11 and P = 1.08e−3, respectively) while Cyanobacteria was significantly higher (P = 1.80e−5) in WO. The most dominant phyla in the fecal microbiome were Bacteroidetes and Firmicutes. This finding shows that the captive (CF) macaques have a significantly higher abundance (P = 4.20e−14) of Bacteroidetes (64.1 ± 6.8%) than wild (WF) macaques (34.4 ± 11.9%). In contrast, the Firmicutes contribute a significantly greater abundance (P = 6.16e−6) in WF (52.8 ± 11.5%) than in the CF (27.6 ± 7.8%). The significant difference in phyla abundance also indicated that the CF has greater proportions of Proteobacteria (P = 1.31e−12) but lower levels of Tenericutes (P = 3.53e−11) than the WF group.

To observe the core structure of the bacterial community in both groups of macaques, the top 18 most dominant taxa were identified and presented in Fig. 2b .

These results indicate the distinct major population of microbial profiles between groups. Streptococcus was the most abundant oral bacteria, contributing the proportion of 31.3 ± 17.2% and 36.1 ± 19.6% in CO and WO, respectively. The other predominant bacteria in the oral microbiome were Actinobacillus, Aggregatibacter, Porphyromonas, Fusobacterium, Gamella, and Haemophilus. The proportion of these bacteria were varied among samples and different between wild and captive groups. The major population of the fecal microbiome in captive (CF) macaques was Prevotella from phylum Bacteroidetes, which made up to 60.5 ± 7.7% of bacterial abundance while it accounted for only 27.4 ± 14.4% in the WF group. The remaining enriched fecal bacteria included the family Ruminococcaceae, Faecalibacterium, family Lachnospiraceae, family S24-7, Bifidobacterium, and Blautia. Relative abundance(%) www.nature.com/scientificreports/ Differential abundance analysis of wild and captive macaque. The differential abundance analysis was performed using the LDA Effective Size (LEfSe) tool 28 . The significant taxa with LDA score > 3 and P < 0.01 were illustrated in Fig. 3 . The oral microbiome of wild (WO) macaques had a greater number of significantly enriched taxa than the captive (CO) macaques, as presented in Fig. 3a . The phylum Bacteroidetes with a high proportion of Prevotella and phylum Spirochaetes were more enriched in the CO microbiome. In addition, the CO group also included high levels of bacteria in the Lactobacillales order, Clostridiales, Campylobaterales, and Pasteurellales. In contrast, the most prominent phyla in WO were Verrucomicrobia, Nitrospirae, and Fusobacteria. The various bacteria in the Alphaproteobacteria class were also more abundant than in the CO. As shown in Fig. 3b , comparing the fecal microbial community, revealed that the fecal microbial community has a lower number of significant taxa than the oral community. We identified the enriched bacteria in CF which were mainly Prevotella, Ruminobacter, Succinivibrio, and order Lactobacillales. In the WF samples, Spirochaetes and Tenericutes were significantly abundant phyla. A significantly greater abundance of Catenibacterium, family Ruminococcaceae and family Clostridiaceae were observed in the WF microbiome.

Functional inference of microbiome in wild and captive macaques. The functional analysis was carried out using PICRUSt2 29 following the standard pipeline implemented in QIIME2 27 . Differential pathway abundance comparison was statistically tested by multiple t-test (P < 0.05) with 1% false discovery rate (FDR) correction. The significant differences of KEGG ortholog (KO) pathway abundance in the oral microbiome were shown in Fig. 4a . The pathway enriched in the WO group is mostly involved with amino acid metabolism while that of the CO group was significantly high in the glycan metabolism. The WO microbiome was also enriched with pathways in lipid metabolism, particularly such as the fatty acid biosynthesis pathway and synthesis and degradation of ketone bodies. The differential functional pathways of the fecal microbiome ( Fig. 4b) showed that most pathways related to energy metabolism were enriched in CF macaques. In contrast, the WF group had high abundance of carbohydrate metabolism including pentose phosphate pathway, fructose and mannose metabolism, galactose metabolism, propanoate metabolism, and glycolysis/gluconeogenesis. Moreover, the WF group also had higher lipid metabolism such as synthesis/degradation of ketone bodies and glycerolipid metabolism than in the CF macaques.

The relationship of Firmicutes/Bacteroidetes ratio and body fat. The Firmicute/Bacteroidetes ratio in fecal microbiome of wild and captive macaques were also compared by t-test with P < 0.05. The wild macaques exhibit the significant higher Firmicute/Bacteroidetes ratio (P = 4.64e−12) as shown in Fig. 5a . The average Firmicute/Bacteroidetes ratios (mean ± SD) were 1.82 ± 0.94 for wild macaques and 0.45 ± 0.19 for captive macaques. The proportion of Firmicute and Bacteroidetes is usually associated with lipid metabolism. Therefore, the skinfolds measured from 3 body sites (belly, suprailiac and subscapular) were analyzed to investigate the relationship between body fat and Firmicute/Bacteroidetes ratio. The skinfold belly, suprailiac and subscapular of wild macaques were significantly higher than those in the captive group (t-test, P < 0.05) as shown in Fig. 5b . The skinfolds of belly and suprailiac (mean ± SD) of wild macaques were marginally higher than the captive macaques 

The comparison of microbiome between wild and captive long-tailed macaques presents the perspective knowledge for exploring the bacterial diversity as well as the effects of habitat and captive environments on their oral and fecal microbiota. In this study, the wild macaques showed a greater richness of oral bacteria than the captive macaques, suggesting the higher microbial exposure in their natural habitat. However, the oral microbiota of both macaques was similarly dominated by Streptococcus which is the most common bacteria in oral healthy macaques 11, 30 . The Streptococci group is composed of various species that contribute to different roles such as carbohydrate metabolizer, acid-base regulator, and inhibitor of pathogenic species 31 . These complex species typically colonize in different sites in the oral cavity depending on their specific niches 30 . Other common oral bacteria identified in these macaques were Actinobacillus, Aggregatibacter, Porphyromonas and Fusobacterium, some of which are pathogenic to humans 30, 32 . However, they are frequently observed in oral microbiota of animals such as chimpanzee and dog 33, 34 . Thus, they may not cause diseases to animals, but further investigation is needed to confirm their pathogenesis abilities.

Our study revealed the different abundance of oral microbes between wild and captive macaques. The Bacteroidetes particularly Prevotella spp. and order Lactobacillales were more abundant in the captive macaque. Another large group of the order Pasteurellales including Actinobacillus, Aggregatibacter, Haemophilus, and Pasteurella were also more enriched in the captive group. The order Pasteurellales are normally commensal microbes on mucous membranes of the upper respiratory tract 35 . They also had been recently identified in oral microbiota of captive baboons 36 . In contrast, the oral microbiota of wild macaques was enriched by the group of environmental bacteria such as phylum Verrucomicrobia, phylum Nitrospirae, Fusobacteria and various members of class Alphaproteobacteria. These were possibly influenced by the environmental exposure allowing the wild macaque to obtain more natural microbes 37 The graph on the left shows the pathway abundances while the right shows 95% confidence intervals of differential abundance. www.nature.com/scientificreports/ were typical bacteria found in soil or freshwater which were mostly autotrophs or symbionts in the rhizosphere or bacterioplanktons. There was no evidence supported that these bacteria can be colonized on the epithelia of oral cavity. A recent publication suggested that the permanent colonizers require the specific environment (such as suitable pH and O 2 level), binding to the host receptor, and crosstalk with other colonizers 30 . Hence, it was hypothesized that several mentioned bacteria were passing microbes which could not colonize and may not be survived in the oral environment.

Since the oral cavity is the entry point of microbial acquisition from the environment, the richness of the oral microbiome was higher than those in fecal microbiota. The limited taxa in fecal microbes might be explained by the selection of the gut environment 44 . Our findings highlighted the lower Shannon diversity of fecal microbiome in captive macaques as a result of the Prevotella enrichment. Nevertheless, the fecal microbiota of both groups shared some common characteristics with the gut microbiome of humans and other NHPs 11, 45, 46 . Therefore, the long-tailed macaque might be a suitable model for investigating the complex interaction and important physiological functions of these gut microbes.

Following the described results, the fecal microbiome of captive macaques was enriched with Prevotella, one of the most predominant genera in the human gut which is also a key feature for the classification of three human enterotypes 47 . A previous study also reported that the captive environment which shifted the diversity and type of diet to polysaccharides can cause the transition of primate microbiome to be dominated with Prevotella 45 . This had also been observed in other captive mammals 24 . The Prevotella spp. was associated with the consumption of high-fiber diets and plant-based carbohydrates especially those found in fruits 48 . Even though all captive macaques were healthy, enrichment of Prevotella spp. could increase host susceptibility to mucosal inflammation as reported previously 49 . This raised an interesting point of concern about the effect of dietary change in captive macaques on animal health. Another predominant genus in captive macaques was Succinivibrio which was also related to the high fiber and complex carbohydrates intake which resembled that of the Prevotella and Treponema 50 . Interestingly, the enrichment of order Lactobacillales and Lactobacillus spp. in the www.nature.com/scientificreports/ captive group may be related to the promoting of bile salt hydrolase activity as well as lactic acid and simple sugar fermentation 51 . In addition, they might contribute to the role in the recycling of urea nitrogen that could facilitate the maintenance of the bodyweight of animals fed with high-fiber and low protein diets 52 .

In contrast, the gut microbiota of wild macaques was mainly dominated by several phyla such as Tenericutes, Spirochaetes and Firmicutes, including especially family Ruminococcaceae and Catenebacterium. The phylum Tenericutes identified in wild macaques mainly in order RF39 associated with fiber digestion as reported in swine 53 . Similarly, Spirochaetes such as Treponema also have cellulolytic capabilities commonly identified in termites 54 and the rumen of grass-fed animals 55 . The higher abundance of family Ruminococcaceae in wild macaques was related to fibrolytic activity in the gut 56 . The members of the family Ruminococcaceae such as Fecalibacterium prausnitzii are known to be a butyrate-producing microorganism that enhances anti-inflammatory effects in the mucosal environment 57 . Depletion of F. prausnitzii was known to be associated with inflammatory bowel disease 58 . The Catenibacterium, another abundant genus in wild macaques in this study, was known to produce the short-chain fatty acids (SCFAs) including acetate, butyrate, and lactate from carbohydrate fermentation that may extract a high yield of energy from carbohydrate-rich diets for the host 59 . The increase in Catenibacteria could induce the energy storage and obesity in humans and mice 60 .

These findings imply that both groups of macaques mainly consumed high-fiber diets supporting the observation that both groups of animals took various kinds of tropical fruits. Basically, the wild macaques foraged for some seasonal fruits in their habitats such as mangoes, figs and bananas while the captive macaques were fed with monkey pellet and fresh fruits. The distinct sources and types of these diets may differently affect the gut microbiota causing the alteration of energy usage and storage.

Generally, the abundance of phylum Firmicutes in the fecal microbiome is known to be associated with fatty acid absorption and lipid metabolism 61 . Interestingly, wild long-tailed macaques were significantly enriched by Firmicutes, so the Firmicute/Bacteroidetes ratio was higher than in the captive macaques. The high abundance of Firmicutes was observed in other wildlife such as sika deer, and rhinoceros 62, 63 . This might probably be the co-adaptation of microbiota and wild animals to enhance the efficacy of energy storage because wild animals can face the situation of starvation or survival in low food abundance. Supporting this hypothesis, the significant association between a greater Firmicutes/Bacteroidetes ratio and higher body fat in wild macaques was detected. Interestingly, the higher Firmicutes/Bacteroidetes ratio in oral microbiome of wild macaque was also observed. The ratio of Firmicutes/Bacteroidetes in oral microbiota was positively correlated with skinfold data similar to that in fecal microbiome. However, there was no previous study related to the association between Firmicutes/ Bacteroidetes ratio in oral microbiota and body fat. Therefore, these results should be carefully interpreted. Further analysis of the relationship between bacteria in oral cavity and host metabolism may be helpful for better elucidating the functions of oral microbiota.

The inference of functional analysis showed the different activities of biochemical pathways in the microbiota of wild and captive macaques. The oral microbiota of captive macaques had a high preference for glycan metabolism such as glycosaminoglycans (GAGs) degradation pathway. A previous study reported that the GAGs degradation properties of bacteria could facilitate their binding to the host epithelial cells 64 . Therefore, most oral bacteria in captive macaques might have more potential to colonize and persist in the oral environment. In wild macaques, the oral microbiota was enriched with amino acid metabolism and lipid biosynthesis, but the sequential outcomes have not been elucidated.

Interestingly, the fecal microbiome of captive macaques increased the citrate cycle (tricarboxylic acid; TCA cycle) and glyoxylate and dicarboxylate metabolism which were suggested to reduce energy consumption and promote energy production 65 . Additionally, the increment of the TCA cycle could oxidize the SCFAs to generate ATP 66 . Therefore, this might cause a lower SCFA uptake into the host and limit energy storage. In contrast, the microbiota in the wild group were involved in ketone body metabolism. It has been reported that gut microbes could enhance ketone metabolism and use it as an energy source during starvation (low blood glucose) in humans 67 . This mechanism could also be the consequence of SCFA production for energy storage as suggested in a previous report 68 . Hence, the microbiota might play important roles in energy reserves and fat accumulation in wild macaques.

The findings of this study provided the preliminary overviews of microbiota in wild long-tailed macaques living in the natural habitat and wild macaques that were transferred from the natural habitat and reared in hygienic housing conditions in the captivity. This study can be conducted only in the countries where the distribution range of long-tailed macaques is and where the primate center in compliance with the national standard (such as AAALAC International) is established. As mentioned above, several captive macaque colonies bred in other primate centers, especially in the US, have usually been facing with the diarrhea problems which might be caused by microbiota dysbiosis 26, 69 . From our field observation, we had not seen the diarrhea in wild animals. Besides, there was no comparison data of the oral-gut microbiota between wild animals living in natural habitat and that wild animals after transferred to hygienic captive condition and depicted how the oral-gut microbiota altered which can help to guide for the good practice of the animal care and management. The finding of this study would be useful for selection the beneficial bacteria for development of probiotics to mitigate the diarrhea in captive long-tailed macaques. For example, the bacteria of family Ruminococcaceae were enriched in wild macaque's gut comparing to that of the captive animals (see Figs. 2, 3) . It was suggested that these bacteria were butyrate-producer that may enhance anti-inflammation in the mucosal environment 57 . Thus, the future study is that the probiotics containing the bacteria complex of family Ruminococcaceae will be produced, fed to diarrhea captive long-tailed macaques, and observe if the symptom is resolved.

In summary, this study presented an overview of the effects of the captive environments on microbiome diversity and the abundance of wild-borne Thai long-tailed macaques. The findings implied that the changes of the environmental exposure and types of diet from the natural habitat to the captivity for only 1 year could alter oral-gut microbiota diversity and abundance in long-tailed macaques. The microbial profiles also suggested their www.nature.com/scientificreports/ putative metabolic functions. However, the interpretation should be cautious due to the lack of metagenomics and metabolomics data in this study. Moreover, the differences of origin and genetic background between the populations could also be a limitation of our study. Nevertheless, the information gained in this study represents a basis for functional microbiome research in long-tailed macaques, implementing the improvement of animal health care. The proper diet feeding program, or a novel probiotics development may be useful to modulate the microbiota protecting against dysbiosis and pathogen infections in these animals in the captive environment.

Animal experiment statement. All methods were performed in accordance with the relevant guidelines and regulations. The animal use protocol was examined and approved by the committee of animal care and use at NPRCT-CU (Protocol Review No. 1775005). Study design and animal specimen collection were carried out and reported following the ARRIVE guidelines (https:// arriv eguid elines. org/ arrive-guide lines).

Animal cohort. Two groups of 80 male and female long-tailed macaques (Macaca fascicularis); wild (n = 37) and captive (n = 43) at the age between 1.5 and 18 years were recruited for this study. The age and body weights of each group of animals were presented in Table 1 . The ages of animals were estimated based on the dental eruption pattern established by Smith et al. 70 Sample collections. The fecal and oral swabs collection of wild macaques was performed in the field site, while that of the captive (wild-born) macaques was performed during the annual health check. The specimens were collected when the macaques were anesthetized with a mixture of dexmedetomidine hydrochloride (Zoetis, USA) (0.03-0.05 mg/kg) and Zoletil (Virbac, New Zealand) (3-5 mg/kg). The polyester tipped swabs (Puritan, USA) was used for specimen collection performed by experienced veterinarians. The oral swab of each individual was obtained by swabbing on both sides of buccal surface as well as lower and upper gum. The fecal specimen was collected by rectal swab directly from animals. The collected swabs were transferred into 3 mL of viral transport media which was fleshly prepared mixtures comprised of Hank's balanced salt solution (Gibco, ThermoFisher Scientific, USA), 1% (w/v) of bovine serum albumin (Gibco, ThermoFisher Scientific, USA), 100 U/mL of penicillin G (Calbiochem, Merck, Germany), 50 μg/mL of streptomycin (Calbiochem, Merck, Germany), and 15 μg/mL of amphotericin B (Supelco, Merck, Germany). The collected samples stored at 4 °C, transferred to the laboratory within 12 h. Then, they were aliquot and kept at − 80 °C until used. In addition, the skinfold measurement was also conducted to record the body fat accumulation in these macaques following the procedure described in previous study 71 .

DNA preparation and 16S microbiome sequencing. Total DNA extraction was performed as previously described in Sawaswong et al. 37 The V4 region of the bacterial 16S small subunit ribosomal RNA (16S rRNA) gene was amplified based on specific primers 16S-515F (5′-GTG CCA GCMGCC GCG GTAA-3′) and 16S-806R (5′-GGA CTA CHVGGG TWT CTAAT-3′) containing phasing sequences and TruSeq adaptor modified from previous study 72 sequences were analyzed using the QIIME2 pipeline (version 2018.8). Briefly, the mate paired reads were joined and quality trimmed (> Q30). The deduplication and de-novo clustering (97% similarity) were performed by the VSEARCH algorithm. The chimeric sequences were eliminated by the UCHIME algorithm. Sequences of OTUs were then classified by comparing against Greengenes bacterial rRNA database (version 13.8) using VSEARCH. Differential abundance comparison of bacteria between captive and wild macaques was investigated by Linear discriminant analysis effect size (LEfSe). The significantly different taxa (P < 0.01) with Linear discriminant analysis (LDA) score > 3 were selected to plot the cladogram. The functional inference of bacterial microbiome was analyzed using the Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt2) software. The KEGG ortholog (KO) 73 abundances of wild and captive macaques were statistically tested by multiple t-tests with 1% false discovery rate (FDR) correction. www.nature.com/scientificreports/

Human gut microbiome: Hopes, threats and promises

The functionality of the gastrointestinal microbiome in non-human animals

Role of the gut microbiota in nutrition and health

Role of the normal gut microbiota

Role of gut microbiota-generated short-chain fatty acids in metabolic and cardiovascular health

An obesity-associated gut microbiome with increased capacity for energy harvest

Activities of compound G2 isolated from alfalfa roots against dermatophytes

Role of the gut microbiota in inflammatory bowel disease pathogenesis: What have we learnt in the past 10 years?

Gut microbiota-mediated inflammation in obesity: A link with gastrointestinal cancer

Environment dominates over host genetics in shaping human gut microbiota

Diversity of macaque microbiota compared to the human counterparts

Microbial ecology along the gastrointestinal tract

Oral microbiome: Unveiling the fundamentals

Experimental approaches for defining functional roles of microbes in the human gut

Comparative microbiome signatures and short-chain fatty acids in mouse, rat, non-human primate, and human feces

The contribution of non-human primate models to the development of human vaccines

The role of the microbiome in drug response

Utility of non-human primates in drug development: Comparison of non-human primate and human drug-metabolizing cytochrome P450 enzymes

Role of the microbiota in the modulation of vaccine immune responses

Survey of the long-tailed macaques introduced onto Tinjil Island, Indonesia

NVX-CoV2373 vaccine protects cynomolgus macaque upper and lower airways against SARS-CoV-2 challenge

Tail length variation in Macaca fascicularis and M. mulatta

Current situation and status of long-tailed macaques (Macaca fascicularis) in Thailand

The effects of captivity on the mammalian gut microbiome

Wild mouse gut microbiota promotes host fitness and improves disease resistance

Idiopathic chronic diarrhea associated with dysbiosis in a captive cynomolgus macaque (Macaca fascicularis)

Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2

Metagenomic biomarker discovery and explanation

PICRUSt2 for prediction of metagenome functions

Role of oral microbiome signatures in diagnosis and prognosis of oral cancer

A practical guide to the oral microbiome and its relation to health and disease

Oral microbiome diversity in chimpanzees from Gombe National Park

Cultivable oral microbiota of domestic dogs

Ferrets, rabbits, and rodents

The microbiome of captive hamadryas baboons

Oral-fecal mycobiome in wild and captive cynomolgus macaques (Macaca fascicularis)

Gemmata obscuriglobus: A connecting link between prokaryotic and eukaryotic cell

Metagenomic achievements in microbial diversity determination in croplands

The Prokaryotes: Alphaproteobacteria and Betaproteobacteria

The diversity of the Limnohabitans genus, an important group of freshwater bacterioplankton, by characterization of 35 isolated strains

The intestinal microbiota: Antibiotics, colonization resistance, and enteric pathogens

Captivity humanizes the primate microbiome

A great-ape view of the gut microbiome

Enterotypes of the human gut microbiome

Review article: Dietary fibre-microbiota interactions

Perturbation of the gut microbiome by Prevotella spp. enhances host susceptibility to mucosal inflammation

Diet, environments, and gut microbiota. A preliminary investigation in children living in rural and urban

Conservation implications of shifting gut microbiomes in captive-reared endangered voles intended for reintroduction into the wild. Microorganisms 6, 94

Experimental evolution on a wild mammal species results in modifications of gut microbial communities

Dynamic distribution of the gut microbiota and the relationship with apparent crude fiber digestibility and growth stages in pigs

Fiber-associated spirochetes are major agents of hemicellulose degradation in the hindgut of wood-feeding higher termites

Persistence of cellulolytic bacteria fibrobacter and treponema after short-term corn stover-based dietary intervention reveals the potential to improve rumen fibrolytic function

Untangling the genetic basis of fibrolytic specialization by lachnospiraceae and ruminococcaceae in diverse gut communities

Obesity and the gastrointestinal microbiota: A review of associations and mechanisms

Association between Faecalibacterium prausnitzii reduction and Inflammatory Bowel Disease: A meta-analysis and systematic review of the literature

In vitro characterization of the impact of different substrates on metabolite production, energy extraction and composition of gut microbiota from lean and obese subjects

The effect of diet on the human gut microbiome: A metagenomic analysis in humanized gnotobiotic mice

The core gut microbiome, energy balance and obesity

Comparison of the gut microbiota composition between wild and captive sika deer (Cervus nippon hortulorum) from feces by high-throughput sequencing

Gut microbiome differences between wild and captive black rhinoceros-Implications for rhino health

Probiotics in human gut microbiota can degrade host glycosaminoglycans

Caloric restriction remodels energy metabolic pathways of gut microbiota and promotes host autophagy

Metabolic tinkering by the gut microbiome: Implications for brain development and function

Regulation of myocardial ketone body metabolism by the gut microbiota during nutrient deprivation

Gut microbiota modulate the metabolism of brown adipose tissue in mice

Longitudinal characterization of Escherichia coli in healthy captive non-human primates

Ages of eruption of primate teeth: A compendium for aging individuals and comparing life histories

Morphological and body color variation in Thai Macaca fascicularis fascicularis North and South of the Isthmus of Kra

Phasing amplicon sequencing on Illumina Miseq for robust environmental microbial community analysis

KEGG as a reference resource for gene and protein annotation

Authors would like to thank the staff of the National Primate Research Center of Thailand, Chulalongkorn University (NPRCT-CU) for macaques' sample collection.

V.S. conceived the study design, performed the analysis, and wrote the manuscript. K.P., P.C. and A.K. performed the laboratory works in 16S sequencing. T.K. captured the macaques and collected the samples for analysis. Y.H. contributed body measurement data of macaques. S.M. supervised the project, edited the manuscript and provided the macaques' specimen. S.P. conceived the study design, supervised the study, and edited the manuscript. All authors reviewed and approved the final manuscript. 

The authors declare no competing interests.

Supplementary Information The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-021-93779-4.Correspondence and requests for materials should be addressed to S.M. or S.P.

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

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:// creat iveco mmons. org/ licen ses/ by/4. 0/.