Skip to main content

The microbiome of captive hamadryas baboons



The hamadryas baboon (Papio hamadryas) is a highly social primate that lives in complex multilevel societies exhibiting a wide range of group behaviors akin to humans. In contrast to the widely studied human microbiome, there is a paucity of information on the host-associated microbiomes of nonhuman primates (NHPs). Here, our goal was to understand the microbial composition throughout different body sites of cohabiting baboons.


We analyzed 170 oral, oropharyngeal, cervical, uterine, vaginal, nasal and rectal samples from 16 hamadryas baboons via 16S rRNA gene sequencing. Additionally, raw Miseq sequencing data from 1041 comparable publicly available samples from the human oral cavity, gut and vagina were reanalyzed using the same pipeline. We compared the baboon and human microbiome of the oral cavity, gut and vagina, showing that the baboon microbiome is distinct from the human. Baboon cohabitants share similar microbial profiles in their cervix, uterus, vagina, and gut. The oral cavity, gut and vagina shared more bacterial amplicon sequence variants (ASVs) in group living baboons than in humans. The shared ASVs had significantly positive correlations between most body sites, suggesting a potential bacterial exchange throughout the body. No significant differences in gut microbiome composition were detected within the maternity line and between maternity lines, suggesting that the offspring gut microbiota is shaped primarily through bacterial exchange among cohabitants. Finally, Lactobacillus was not so predominant in baboon vagina as in the human vagina but was the most abundant genus in the baboon gut.


This study is the first to provide comprehensive analyses of the baboon microbiota across different body sites. We contrast this to human body sites and find substantially different microbiomes. This group of cohabitating baboons generally showed higher microbial diversity and remarkable similarities between body sites than were observed in humans. These data and findings from one group of baboons can form the basis of future microbiome studies in baboons and be used as a reference in research where the microbiome is expected to impact human modeling with baboons.


Humans and other primates are home to trillions of symbiotic microorganisms. Interactions between a host and its microbes affect host physiology, behavior, reproduction, immunity and evolution [1,2,3]. The Human Microbiome Project, through monitoring or manipulations of the human microbiome, helps us better understand the associations between microbes and human health [4]. In contrast to the widely studied human microbiome, there is a paucity of information on the host-associated microbiomes of nonhuman primates (NHPs). Information about NHP microbiota is essential for understanding the factors underlying microbial coevolution with their hosts [5, 6]. Broad primate microbiome surveys could also allow for the development of predictive biomarkers to improve nonhuman primate health and management.

Baboons (genus Papio) are one of the most biologically relevant research animal models due to their genetic and physiological similarities to humans [7]. Baboons are large-bodied, omnivorous, highly social, terrestrial Old World African monkeys that occupy a wide array of habitats similar to those of early hominins [8, 9]. Of the six recognized species [10], the social system of hamadryas baboons shares more similarities with humans than that of other baboons [9]. Like modern humans, the hierarchical social networks of hamadryas baboons connect individuals at multiple levels [8]. Frequent social interactions (mostly grooming) are necessary for baboons to maintain affiliative bonds [11].

To our knowledge, only the rectal and vaginal microbiota of baboons have been examined, likely because baboons can be used as a model in gastrointestinal and female reproductive studies due to the features shared with humans [12,13,14]. In this study, we investigated 170 samples of rectal, oral, oropharyngeal, cervical, uterine, vaginal and nasal microbiota from 16 captive hamadryas baboons by culture-independent sequencing of the 16S rRNA gene hypervariable V3-V4 region. Our study provides detailed insights into the baboon microbiome structure and ecology.


Microbial distribution in different body sites of captive baboons

Six main phyla were detected in the seven body sites of 16 captive hamadryas baboons (Papio hamadryas, 13 females, 3 males). The phylum Spirochaetae was found to be abundant in the baboon gut (Fig. 1). Firmicutes dominated in all body sites. The dominance of other phyla varied among body sites; for instance, Fusobacteria was dominant in the oral cavity and oropharynx, while Actinobacteria predominated in the nose. On the order level, we found Lactobacillales to be abundant in all seven body sites (Mean 27% ± SD 10%), especially in the oral cavity, oropharynx and nose, constituting around 40% of the total bacteria. In addition to Lactobacillales, Clostridiales were predominant in the vagina, cervix, uterus, and gut, accounting for 16–27% of the microbiota. The oral cavity and oropharynx shared a microbial profile mainly composed of Lactobacillales, Bacteroidales, and Pasteurellales, representing over 50% of the microbiota. Lactobacillales and Corynebacteriales were the two major bacteria in the nose, together comprising almost 70% of the microbiota.

Fig. 1
figure 1

Microbial composition on phylum and order levels of each anatomic site. Pie charts show the average microbial composition of seven body habitats on order (donut chart) and phylum level (pie chart), respectively. Lactobacillales were found in all the body sites (Mean 27% ± SD 10%)

Microbial characterization varied among the body sites

To investigate microbial features of different body sites, we analyzed microbial diversity between and within the different body sites, quantified the microbiome divergence within each body site, and compared gut microbial diversity between and within maternal lines (Fig. 2). Permutational Multivariate Analysis of Variance (PERMANOVA) using Bray-Curtis distances, visualized by Non-metric multidimensional scaling (NMDS) (Fig. 2a) showed that the oral cavity, oropharynx, and nose had unique microbial profiles but the microbial profiles in the cervix, uterus, vagina, and rectum did not have significant differences (pairwise comparison, adjusted p values are listed in Table S1). The same NMDS plot but showing the host source is shown in Supplementary Fig. 1. Microbial alpha-diversity in the baboon nose was significantly lower than in the other body sites (adjusted p values are listed in Table S2, TukeyHSD) (Fig. 2b). The divergence of microbiomes within each body site was quantified and extrapolated (Fig. 2c). The microbiomes in baboon nose and oropharynx had a smaller dispersion but were more heterogeneous in the reproductive tract and oral cavity (adjusted p values are listed in Table S3, Pairwise Wilcoxon Rank Sum Tests). The vertical inheritance of gut microbial communities was analyzed by comparing the weighted Unifrac distances between microbiomes of the same maternal line and microbiomes of different maternal lines (Fig. 2d). Statistical analysis showed no significant difference in weighted Unifrac distance within the maternal lines compared to between the maternal lines (t-test, p = 0.194), indicating the vital role of horizontal exchange in shaping the gut microbiota.

Fig. 2
figure 2

Microbial characterization in different body sites. a Non-metric multidimensional scaling (NMDS) based on Bray-Curtis distance of microbial communities from the baboon oral cavity, oropharynx, nose, cervix, uterus, vagina and rectum. The colored lines surrounding each sample type are covariance ellipsoids. b Alpha diversity in different body habitats, grouped by area, as measured using the Shannon index of ASV-level bacteria. Pharynx and nose had a significantly lower diverse microbiota than reproductive tract (TukeyHSD). c Divergence of microbes in a specific body site was quantified as the average dissimilarity of each sample from the group mean. d Weighted Unifrac distance within and between maternal lines for gut microbiota. ns means no significant difference detected between the same maternal line and different maternal lines by t-test (p = 0.19)

Baboons had significantly different vaginal, gut and oral microbiomes from humans

The sequencing data from 1041 human oral cavity, gut and vaginal samples (EBI databases available through study ID PRJEB14941) [15], processed using similar DNA extraction, DNA kit and sequencing platform, were analyzed using the same pipeline as used for the baboon. To ensure that the sample types were comparable between baboons and humans, we compared the homogeneity of dispersions for baboons and humans for each sample type. The average dispersion within sample types from baboons and human did not show significant differences (baboon gut: 0.4, human gut: 0.5, baboon oral: 0.3, human oral: 0.37, baboon vagina: 0.5, human vagina: 0.43, adjusted p value > 0.5, adjusted p values for all pairwise comparisons are shown in Supplementary Table S4). We found that the baboon microbiomes in the three body sites were significantly different from human microbiomes (pairwise comparison, adjusted p = 0.015 for oral vs oral microbiome, gut vs gut microbiome and vaginal vs vaginal microbiome, PERMANOVA, Fig. 3a).

Fig. 3
figure 3

Comparisons between the human and baboon microbiomes in the oral cavity, gut and vagina. a Non-metric multidimensional scaling (NMDS) based on a Bray-Curtis distance matrix of microbial communities from baboon and human oral, gut and vagina. The colored lines surrounding each sample type are covariance ellipsoids. b Microbial alpha diversity in the human and baboon gut, oral cavity and vagina, as measured using the Shannon index of ASV-level bacteria. c Heatmap showing the 30 most abundant bacterial genera in human and baboon oral, nose and vagina

The baboon vaginal microbiome had a higher alpha diversity than that of the human (adjusted p value <1e-7, TukeyHSD, Fig. 3b). The oral microbiome had a slightly higher alpha diversity in humans than in the baboon, while the gut microbiome showed similar alpha diversity in the two hosts. Lactobacillus (13% ± 12%) was not dominant in the baboon vagina, as it is in the human vagina (54% ± 37%), but it was the most abundant genus in the baboon gut (22% ± 17%) (Fig. 3c). In contrast, the human gut only contained a mean relative abundance of 0.35% Lactobacillus (Fig. 3c). Despite an overall significant difference in oral microbiome profiles between baboons and humans (Fig. 3a), they clustered together based on the abundances of the 30 most abundant genera overall (Fig. 3c). Streptococcus was the most abundant genus both in the human and baboon oral cavity.

The microbiomes of group living baboons shared more similarities across body sites than in humans

Bacterial exchange during group living is inevitable. We define shared ASVs as those that were present in at least two or more of the body sites sampled. The different body sites of group living baboons shared more ASVs than those of humans (Fig. 4a). The abundances of these shared ASVs were not significantly correlated between the body sites in humans but had a significant positive correlation between gut and vagina in the baboons (Fig. 4a), indicating potential bacterial exchange. Of the 15 ASVs shared between the human oral cavity, gut and vagina, 5 ASVs belonged to unclassified taxa (Fig. 4a). Most of the 94 ASVs shared by the baboon’s oral cavity, gut and vagina belonged to Firmicutes (Fig. 4a). In addition, the human oral cavity, gut and vagina had unique microbial compositions as shown in Non-metric multidimensional scaling (NMDS) plot (pairwise comparison, adjusted p = 0.003 for oral vs gut microbiome, oral vs vaginal microbiome and gut vs vaginal microbiome, PERMANOVA, Fig. 3a). However, the baboon gut microbiome was similar to the vaginal microbiome (Fig. 2a) despite a large variance within a single individual (Fig. 4b). From Fig. 2a, it is clear that multiple body sites of baboons had similar microbial profiles. Considering the potential bacterial exchange across body sites among co-habitants, we analyzed the ASVs (Fig. 4c) shared by the 7 body sites and their correlations (Fig. 4d). We identified 35 shared ASVs representing 33% (17–70% in each body habitat) of the mean relative abundance (Fig. 4c) and the majority belonged to the phylum Firmicutes (Fig. 4c). Among these 35 shared ASVs, 11 ASVs belong to the genus Lactobacillus and 2 ASVs belong to the genus Faecalibacterium. The relative abundance of these 35 shared ASVs was positively correlated between cervix, uterus, vagina and rectum, between nose, cervix, uterus, and vagina, as well as between oral and oropharynx.

Fig. 4
figure 4

Shared microbial ASVs and their correlation analysis among human/baboon oral, gut and vagina, and among all seven body sites of baboons. a Venn diagrams showing the number of shared ASVs and spearman correlation matrix for shared ASVs among oral, gut and vagina of humans and baboons. X in the correlation matrix means no significant correlations. All the significant correlations shown in the matrix have an adjusted p-value< 0.01. Color intensity and the size of the circle are proportional to the correlation coefficients. The shared ASVs by human and baboon samples were plotted in the circular barplots. b Bar plots showing the 40 most abundant bacteria ASVs of vagina and rectum in eight baboons. c Heatmap for 35 ASVs shared by all body sites. The ASVs belonging to genus Lactobacillus and Faecalibacterium are marked in red. The relative abundances of Lactobacillus and the 35 ASVs are listed. d Spearman correlation analysis for 35 ASVs shared by seven body sites. Only the Spearman correlation coefficients (PCC) with an adjusted p < 0.01 were plotted. X in the correlation matrix means no significant correlations. Color intensity and the size of the circle are proportional to the correlation coefficients


In our study, we found that the baboon microbiome has unique characteristics compared to the human microbiome. The microbiomes of group living baboons shared more similarities between body sites than observed in humans. A significant positive correlation of the relative abundance of shared ASVs among the different body sites suggested potential bacterial exchange throughout the body. We note that the baboons investigated in our study lived in the same environment and shared the same diet, and acknowledge that both diet [16] and to a lesser extent genetics [17, 18] are known to be important factors in shaping the human gut microbiome, and may have contributed towards the high similarity of baboon microbiomes found here.

Group living generally entails frequent social interactions, especially for highly social baboons [19], who use them to maintain affiliative bonds [11]. Bacterial exchange in a shared environment has previously been reported. Members of a shared environment have more similar gut and skin bacterial communities than individuals living separately, indicating that a shared lifestyle or environment affects the microbiome composition [5, 6, 18, 20]. In a study of the gut microbiota in wild chimpanzees, group living chimpanzees shared more of their gut microbiome than individuals from different groups, and dietary convergence did not explain the convergence among chimpanzee-gut microbiome community memberships during periods of high sociability [6]. In the study by Tung et al., the yellow baboon (Papio cynocephalus) gut microbiome was detected to have a high group specificity as well [5]. They found that even if two social groups of wild baboons shared almost the same diet, the gut microbiota of these two groups was still significantly different and shaped by the social interactions [5]. Therefore, we speculate that bacterial exchange was an important reason for microbial similarities between multiple body sites of the group living baboons investigated here. In addition, the gut microbiota within the three maternal lines did not show a significant difference in comparison to those of unrelated individuals, which indicated that the gut microbiota is more affected by the bacterial exchange between group members in general than transfer from mother to offspring. Our findings are consistent with the results previously reported [6]: inheritance of microbial communities across generations were primarily driven by horizontal transfer among interacting hosts. In our study, approximately one-third of all the shared ASVs were from the genera Lactobacillus and Faecalibacterium, indicating a higher tendency for these genera to be shared between body sites. Many species of Lactobacillus and Faecalibacterium are widely considered to be probiotics [21, 22], suggesting sociality may also foster commensal and mutualistic microbial diversity, which could potentially help explain the driving force for bacterial transmission.

In this study, we found that the baboon vaginal microbiome was similar to the gut microbiome, which is different from what has been seen in humans. For baboons, some grooming bouts, especially those directed from adult males towards estrous females, concentrate heavily on the anogenital region, increasing the probability of fecal-vagina transfer [5]. Besides, we found that the baboon vagina had a distinct microbiota profile from that of humans (Fig. 3a). Compared to humans, the polygynous mating system and promiscuity of the baboons probably boosted genital bacterial transfer among group members. The human vaginal samples analyzed in our study were from pregnant women and some were given a yogurt supplement [15]. There are differences in the composition and stability of the vaginal microbiota between pregnant and non-pregnant women [23], which can influence the comparison between human gut and vaginal microbial compositions. Nevertheless, normal healthy humans are known to have completely different gastrointestinal and urogenital microbiomes [24], which is different from what we find in the baboons. The baboon vagina had a relatively high microbial diversity (Fig. 3b) and Lactobacillus was not as predominant in the baboon vagina as in humans (Fig. 3c). Similar findings have been reported in other studies [13]. The human vagina is primarily colonized by Lactobacillus [25], which maintains an acidic environment and prevents the invasion of nonindigenous strains and potential pathogens and can account for 65.9 to 98.1% of the vaginal microbiota [26,27,28,29]. However, to some extent, the baboon vaginal microbial profiles were characterized by low Lactobacillus abundance, low lactic acid concentration and a higher, near-neutral vaginal pH [30], which are typically associated with bacterial vaginosis in the human vagina [31].

The baboon oral cavity had a unique microbial distribution and exhibited only minor overlap with the other body sites (Fig. 2a); it had somewhat lower microbial diversity than humans (Fig. 3b), despite the fact that humans exhibit more oral hygiene practices. Cleaner teeth could be argued to lead to higher diversity, as no bacteria will be allowed to grow to dominate the oral cavity, enabling a higher diversity of bacteria to be present at any given moment. The baboon gut microbiome is unique compared to the human gut microbiome. The baboon gut had high relative abundances of Lactobacillales and Clostridiales, which was consistent with previous findings [32]. In this previous study, the gut microbiome of captive olive baboons (Papio Anubis) was also reported to be significantly different from human gut microbiome. In our study, Spirochaetes, which was extremely rare in the modern human but enriched in ancient humans [33], was also enriched in the baboon gut (Fig. 1). Lactobacillus was the most abundant genus in the baboon gut but presented a low abundance in the human gut (Fig. 3c). Notably, the important human gut bacterium, Akkermansia [34], was not detected in the baboon gut. Therefore, Akkermansia was more human-specific and thus absent in the baboons in our study, no Akkermansia has been reported in NHPs [35, 36].


To our knowledge, this study is the first to investigate microbial compositions in cohabitating baboons across different body sites. Our results showed that baboons have a unique microbiome compared to humans. The microbial diversity of the baboon vagina was much higher than that of humans. Lactobacillus was not so predominant in baboon vagina as in the human vagina but was the most abundant genus in the baboon gut. The microbial compositions in the baboon reproductive tract and gut were similar. Oral cavity, vagina and gut in group living baboons shared more bacterial ASVs than humans. The significantly positive correlations of those shared ASVs between multiple body sites in this group of baboons combined with highly social characteristics of baboons indicated a potential bacterial exchange throughout the body. We reported that the probably transmitted bacteria across body sites tend to be bacteria known to be beneficial in humans, which may suggest that some modern human populations, due to changed social behaviors, may have lost an important source of beneficial microbiota with consequences for human health.

Materials and methods

Controls manage the risk of contamination during wet-lab processing and sterile surgical procedure manage the risk of contamination at sampling

To avoid contamination risks, we strictly controlled the sampling process, DNA extraction, PCR and sequencing. All DNA extractions strictly followed the aseptic operation process under a clean bench. We also have DNA extraction negative control (from the DNA extraction to the sequencing process), sequencing blank control (clean water for sequencing), and sequencing positive control (mock community, E.coli). The detailed information is included in the supplementary material.

Sample collection

Samples were collected from 16 captive baboons (Papio hamadryas, 13 females, 3 males) housed at the Copenhagen Zoo, Denmark. Seven different sites were sampled (Fig. 1) and all sample information is listed in Table S5. Animals were anesthetized for a full medical evaluation and physical examination. Non-invasive samples (vagina, nose, oral, oropharynx, and gut) were collected in a sheltered housing facility using sterile polyester swabs (cat no. 300263, Deltalab, Spain). Following thorough medical evaluations, eight of the animals were euthanized by a licensed veterinarian and a thorough postmortem examination was conducted in a separate necropsy room. Carcasses were opened ventrally to expose the organs, and all invasive sampling was performed sequentially from cranial to caudal. A new sterile scalpel was used for each organ; new gloves were donned and new surgical utensils were used for each of the carcasses. Animals were euthanized for reasons unrelated to this project.

DNA extraction and 16S rRNA gene amplicon sequencing

Genomic DNA was extracted from swab samples with the PowerLyzer® PowerSoil® DNA Isolation Kit (MO-BIO Laboratories, Inc., Carlsberg, CA, USA), and 50 μL of elution buffer was used for each sample. All operations were performed under aseptic conditions. Extracted DNA was stored at − 20 °C. Sterilized PBS solution and Molecular grade water (Sigma-Aldrich, United States) were used as DNA extraction and DNA amplification negative control, whereas mock community and only Escherichia.coli strain were included in all the following steps as a positive control. The 16S rRNA gene hypervariable V3-V4 region was amplified with 2 μL template DNA, using 0.25 μL Phusion high-fidelity (HF) DNA Polymerase (Thermo Fisher Scientific, Waltham, MA, USA), 5 μL 5 × Phusion buffer HF, 0.5 μL 10 mM dNTPs, 1 μL 10 μM of each primer (the modified broad primers 341F (5′-CCTAYGGGRBGCASCAG-3′) and Uni806R (5′-GGACTACNNGGGTATCTAAT-3′) [37] in a 25 μL PCR reaction volume. The first PCR program included 30 s at 98 °C, 30 cycles of 5 s at 98 °C, 15 s at 56 °C, and 72 °C for 10 s, and then 5 min at 72 °C. In the second PCR, sequencing primers and adaptors were attached to the amplicon library following the first PCR conditions with only 15 cycles. The size of the PCR product (≈466 bp) was evaluated using gel electrophoresis. The amplicon products were purified by use of Agencourt AMPure XP beads (Beckman Coulter Genomics, MA, USA) with the 96-well magnet stand, normalized by the SequalPrepTM Normalization Plate (96) kit (Invitrogen Ltd., Paisley, UK), pooled in equimolar concentrations and concentrated using the DNA Clean & Concentrator™-5 Kit (Zymo Research, Irvine, CA, USA). Sequencing of the amplicon library was performed on the Illumina MiSeq System with MiSeq reagent kit v2 (Illumina Inc., CA, USA), including 5.0% PhiX as an internal control.

Bioinformatics analysis and statistical analysis

The raw fastq files were demultiplexed using the Miseq Controller Software. Primers and diversity spacers were removed from fastq files using “Cutadapt” [38]. The data trimming and feature classification were done using QIIME 2 Core 2017.12 distribution microbiota analysis platform [39]. Paired-end sequences were merged by vsearch plugin [40] and then followed by filtering with the quality-filter plugin [41], both with default settings. Deblur plugin was then used to denoise the sequences with a trim length of 400 bp based on quality score plots [42]. Sequence alignments were generated using MAFFT and the aligned sequences were masked by MASK plugin [43]. FastTree and midpoint-root built-in phylogeny plugin were used to create a rooted phylogenetic tree [44]. Pre-fitted sklearn-based taxonomy classifier ( was used to blast representative sequences against silva 132 database for taxonomic classification of features [45]. Rarefaction curves were plotted by workflow in QIIME 2.

A histogram of average unweighted UniFrac distances between each sample and all the rest of the samples were plotted to confirm that PCR and Sequencing controls differed from our samples (supplementary material). The DNA amplification negative control with only 63 reads was filtered. The histogram showed that the three controls (DNA extraction blank control, E.coli strain and Mock control) were far away from the real samples, which meant controls were different from the real samples. The rarefaction curves (supplementary material) demonstrated the observed richness in a given count of sequences. Observed richness curves reached asymptotes after 4000 reads for most samples. With an average of 15,247 clean sequences per sample, sufficient sequences for all 170 samples were generated to characterize the microbial community in the seven body habitats.

The open-source statistical program “R” was used for data treatment and statistical analysis [46], predominantly the R-package “phyloseq” [47]. Alpha diversity between the groups was tested by analysis of variance using the function “ano”. If significant differences between the groups were present, multiple comparisons with the function “TukeyHSD” were performed in a pairwise manner between all groups (all two functions from R-package “stats”). Bray Curtis distance was used to explain differences among microbial communities and the dissimilarity was examined by permutational multivariate analysis of variance (PERMANOVA, “vegan” function “adonis”) [48]. R Function “pairwise.adonis” [49] was used for multiple comparisons and Bonferroni correction was used to account for multiple comparisons. Group Divergence was quantified as the average dissimilarity of each sample from the group mean by using the function “divergence” from R-package “microbiome” [50]. Venn diagram was plotted by R function “VennCounts” and “VennDiagram” in R package “VennDiagram” [51]. R function “rcorr” was used to compute the Spearman correlation analysis and the significance levels [52]. The correlation matrix and the significance test were visualized by R function “corrplot” [53]. Pie chart, donut chart, violin plot, box plot, heatmap, bar chart and circle bar chart were plotted using ggplot2 [54].

Availability of data and materials

The dataset analyzed during the current study will be available, upon publication, in the Sequence Read Archive (SRA) repository, All the scripts used in this study can also be found at,



Nonhuman primates


Amplicon sequence variant


Permutational Multivariate Analysis of Variance


Non-metric multidimensional scaling


Multiple sequence alignment program


Polymerase chain reaction


  1. Cho I, Blaser MJ. The human microbiome: at the interface of health and disease. Nat Rev Genet. 2012;13:260–70.

    Article  CAS  Google Scholar 

  2. Pflughoeft KJ, Versalovic J. Human microbiome in health and disease. Annu Rev Pathol Mech Dis. 2012;7:99–122. Available from:

    Article  CAS  Google Scholar 

  3. Ezenwa VO, Gerardo NM, Inouye DW, Medina M, Xavier JB. Animal behavior and the microbiome. Science. 2012;338(6104):198–9.

  4. Nih T, Working HMP. The NIH human microbiome project. Genome Res. 2009;19:2317–23 Available from:

    Article  Google Scholar 

  5. Tung J, Barreiro LB, Burns MB, Grenier JC, Lynch J, Grieneisen LE, et al. Social networks predict gut microbiome composition in wild baboons. Elife. 2015;4:e05224.

  6. Moeller AH, Foerster S, Wilson ML, Pusey AE, Hahn BH, Ochman H. Social behavior shapes the chimpanzee pan-microbiome. Sci Adv. 2016;2:e1500997.

    Article  Google Scholar 

  7. Rogers J, Hixson JE. Baboons as an animal model for genetic studies of common human disease. Am J Hum Genet. 1997;61:489–93.

    Article  CAS  Google Scholar 

  8. Barrett L, Henzi SP. Baboons. Curr Biol. 2008;18:404–6.

    Article  Google Scholar 

  9. Swedell L, Plummer T. A Papionin multilevel society as a model for hominin social evolution. Int J Primatol. 2012;33:1165–93.

    Article  Google Scholar 

  10. Newman TK, Jolly CJ, Rogers J. Mitochondrial phylogeny and systematics of baboons (Papio). Am J Phys Anthropol. 2004;124:17–27.

    Article  Google Scholar 

  11. Lehmann J, Korstjens AH, Dunbar RIM. Group size, grooming and social cohesion in primates. Anim Behav. 2007;74:1617–29.

    Article  Google Scholar 

  12. Mckenney EA, Ashwell M, Lambert JE, Fellner V. Fecal microbial diversity and putative function in captive western lowland gorillas (Gorilla gorilla gorilla), common chimpanzees (pan troglodytes), Hamadryas baboons (Papio hamadryas) and binturongs (Arctictis binturong). Integr Zool. 2014;9:557–69.

    Article  Google Scholar 

  13. Yildirim S, Yeoman CJ, Janga SC, Thomas SM, Ho M, Leigh SR, et al. Primate vaginal microbiomes exhibit species specificity without universal lactobacillus dominance. ISME J. 2014;8:2431–44.

    Article  CAS  Google Scholar 

  14. Miller EA, Livermore JA, Alberts SC, Tung J, Archie EA. Ovarian cycling and reproductive state shape the vaginal microbiota in wild baboons. Microbiome. Microbiome. 2017;5:1–14Available from:.

    Article  Google Scholar 

  15. Bisanz JE, Enos MK, PrayGod G, Seney S, Macklaim JM, Chilton S, et al. Microbiota at multiple body sites during pregnancy in a rural tanzanian population and effects of Moringa-supplemented probiotic yogurt. Appl Environ Microbiol. 2015;81:4965–75.

    Article  CAS  Google Scholar 

  16. David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, et al. Diet rapidly and reproducibly alters the human gut microbiome. NatureNature Publishing Group. 2014;505:559–63.

    Article  CAS  PubMed  Google Scholar 

  17. Goodrich JK, Waters JL, Poole AC, Sutter JL, Koren O, Blekhman R, et al. Human genetics shape the gut microbiome. CellElsevier Inc. 2014;159:789–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Rothschild D, Weissbrod O, Barkan E, Kurilshikov A, Korem T, Zeevi D, et al. Environment dominates over host genetics in shaping human gut microbiota. Nat Publ GrNature Publishing Group. 2018;555:210–5.

    Article  CAS  Google Scholar 

  19. Schreier AL, Swedell L. The fourth level of social structure in a multi-level society: ecological and social functions of clans in Hamadryas baboons. Am J Primatol. 2009;71(11):948–55.

  20. Song SJ, Lauber C, Costello EK, Lozupone CA, Humphrey G, Berg-Lyons D, et al. Cohabiting family members share microbiota with one another and with their dogs. Elife. 2013;2:e00458.

  21. Liévin-Le Moal V, Servin AL. Anti-infective activities of lactobacillus strains in the human intestinal microbiota: from probiotics to gastrointestinal anti- infectious biotherapeutic agents. Clin Microbiol Rev. 2014;27:167–99.

    Article  Google Scholar 

  22. Miquel S, Martín R, Rossi O, Bermúdez-Humarán LG, Chatel JM, Sokol H, et al. Faecalibacterium prausnitzii and human intestinal health. Curr Opin Microbiol. 2013;16:255–61.

    Article  CAS  Google Scholar 

  23. Romero R, Hassan SS, Gajer P, Tarca AL, Fadrosh DW, Nikita L, et al. The composition and stability of the vaginal microbiota of normal pregnant women is different from that of non-pregnant women. Microbiome. 2014;2(1):4.

  24. Huttenhower C, Gevers D, Knight R, Abubucker S, Badger JH, Chinwalla AT, et al. Structure, function and diversity of the healthy human microbiome. NatureNature Publishing Group. 2012;486:207–14.

    Article  CAS  Google Scholar 

  25. Srinivasan S, Fredricks DN. The human vaginal bacterial biota and bacterial vaginosis. Interdiscip Perspect Infect Dis. 2008;2008:1–22 Available from:

    Article  Google Scholar 

  26. Kim TK, Thomas SM, Ho M, Sharma S, Reich CI, Frank JA, et al. Heterogeneity of vaginal microbial communities within individuals. J Clin Microbiol. 2009;47:1181–9.

    Article  CAS  Google Scholar 

  27. Hummelen R, Fernandes AD, Macklaim JM, Dickson RJ, Changalucha J, Gloor GB, et al. Deep sequencing of the vaginal microbiota of women with HIV. PLoS One. 2010;5:e12078.

  28. Srinivasan S, Liu C, Mitchell CM, Fiedler TL, Thomas KK, Agnew KJ, et al. Temporal variability of human vaginal bacteria and relationship with bacterial vaginosis. PLoS One. 2010;5:e10197.

  29. Ravel J, Brotman RM, Gajer P, Ma B, Nandy M, Fadrosh DW, et al. Daily temporal dynamics of vaginal microbiota before, during and after episodes of bacterial vaginosis. Microbiome. 2013;1:29 Available from:

    Article  Google Scholar 

  30. Miller EA, Beasley DAE, Dunn RR, Archie EA. Lactobacilli dominance and vaginal pH: why is the human vaginal microbiome unique? Front Microbiol. 2016;7:1–13.

    Google Scholar 

  31. Aldunate M, Srbinovski D, Hearps AC, Latham CF, Ramsland PA, Gugasyan R, et al. Antimicrobial and immune modulatory effects of lactic acid and short chain fatty acids produced by vaginal microbiota associated with eubiosis and bacterial vaginosis. Front Physiol. 2015;6:164.

  32. Tsukayama P, Boolchandani M, Patel S, Pehrsson EC, Gibson MK, Chiou KL, et al. Characterization of wild and captive baboon gut microbiota and their antibiotic Resistomes. mSystems. 2018;3:e00016-18.

  33. Turroni S, Rampelli S, Centanni M, Schnorr SL, Consolandi C, Severgnini M, et al. Enterocyte-associated microbiome of the Hadza hunter-gatherers. Front MicrobiolNature Publishing Group. 2016;7:1–12.

    Article  CAS  Google Scholar 

  34. van Passel MWJ, Kant R, Zoetendal EG, Plugge CM, Derrien M, Malfatti SA, et al. The genome of Akkermansia muciniphila, a dedicated intestinal mucin degrader, and its use in exploring intestinal metagenomes. PLoS One. 2011;6:e16876.

  35. Li X, Liang S, Xia Z, Qu J, Liu H, Liu C, et al. Establishment of a Macaca fascicularis gut microbiome gene catalog and comparison with the human, pig, and mouse gut microbiomes. GigascienceOxford University Press. 2018;7:1–10.

    PubMed  Google Scholar 

  36. Firrman J, Liu LS, Tanes C, Bittinger K, Mahalak K, Rinaldi W. Metagenomic assessment of the Cebus apella gut microbiota. Am J Primatol. 2019;81(10–11):e23023.

  37. Takai K, Horikoshi K, Takai KEN. Rapid detection and quantification of members of the archaeal community by quantitative PCR using Fluorogenic probes rapid detection and quantification of members of the archaeal community by quantitative PCR using Fluorogenic probes. Appl Environ Microbiol. 2000;66:5066–72.

    Article  CAS  Google Scholar 

  38. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal. 2011;17:10 Available from:

    Article  Google Scholar 

  39. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852-7.

  40. Rognes T, Flouri T, Nichols B, Quince C, Mahé F. VSEARCH: a versatile open source tool for metagenomics. PeerJ. 2016;4:e2584 Available from:

    Article  Google Scholar 

  41. Bokulich NA, Subramanian S, Faith JJ, Gevers D, Gordon JI, Knight R, et al. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat Methods. 2013;10:57–9.

    Article  CAS  Google Scholar 

  42. Amir A, McDonald D, Navas-Molina JA, Kopylova E, Morton JT, Zech Xu Z, et al. Deblur rapidly resolves single-nucleotide community sequence patterns. mSystems. 2017;2:e00191–16 Available from:

    PubMed  PubMed Central  Google Scholar 

  43. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.

    Article  CAS  Google Scholar 

  44. Price MN, Dehal PS, Arkin AP. FastTree 2 - approximately maximum-likelihood trees for large alignments. PLoS One. 2010;5:e9490.

  45. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590-D596.

  46. core Team R. R: a language and environment for statistical computing. Vienna: R Found Stat Comput; 2018.

    Google Scholar 

  47. McMurdie PJ, Holmes S. Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8:e61217.

  48. Oksanen J, Kindt R, Legendre P, O’Hara B, Simpson GL, Solymos PM, et al. The vegan package. Community Ecol Packag. 2008;190 Available from:

  49. Arbizu M. Pairwiseadonis: pairwise multilevel comparison using Adonis. R Packag version 00; 2017. p. 1.

    Google Scholar 

  50. Lahti L, Shetty S, Blake T. Tools for microbiome analysis in R. microbiome Packag version 099; 2017.

    Google Scholar 

  51. Chen H, Boutros PC. VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinformatics. 2011;12:1–7.

  52. Miscellaneous TH, Yes L. Package ‘Hmisc’; 2018.

    Google Scholar 

  53. Wei T, Simko V. The corrplot package. R Core team; 2016. Available from:

    Google Scholar 

  54. Wilkinson L. ggplot2: elegant graphics for data analysis by WICKHAM, H. Biometrics. 2011;67:678–9. Available from:.

    Article  Google Scholar 

Download references


We express our deepest gratitude to Copenhagen Zoo for collaboration with the baboon samples. We acknowledge Luma Odish for the help and support with the construction of the 16S rRNA gene amplicon libraries and sequencing.


This work was supported by the China Scholarship Council and research grants by the Danish Council for Independent Research (DFF-1323-00235, SIMICOM), and Novo Nordisk Foundation (NNF17OC0025014).

Author information

Authors and Affiliations



The sampling was carried out by UT and MFB. XJL performed the DNA extraction and sequencing. XJL performed the microbiota analysis under the supervision of ADB and GV. SJS, GV and MSM helped interpret the data. XJL wrote most of the paper; the sampling process is written by UT and MFB. This project was conceived and designed by SJS, UT and XJL. The authors read, revised, and approved the final manuscript.

Corresponding author

Correspondence to Søren Johannes Sørensen.

Ethics declarations

Ethics approval and consent to participate

No baboons were euthanized for the purpose of this study. Rather, we were allowed to obtain samples opportunistically in conjunction with routine health checks. As for the invasive samples, they were obtained from individuals euthanized for management reasons (social restructuring).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no conflict of interest.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1. Supplementary table

: sample overview; statistical test tables.

Additional file 2. Supplementary material

: controls manage risk of contamination during wet-lab processing and sterile surgical procedure manage risk of contamination at sampling.

Additional file 3.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, X., Trivedi, U., Brejnrod, A.D. et al. The microbiome of captive hamadryas baboons. anim microbiome 2, 25 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: