Antibiotic Resistance Profiles of Gut Microbiota Across Various Primate Species in Guangxi

doi:10.21203/rs.3.rs-3334710/v1

Download PDF

Research Article

Antibiotic Resistance Profiles of Gut Microbiota Across Various Primate Species in Guangxi

https://doi.org/10.21203/rs.3.rs-3334710/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Understanding the gut microbiota and antibiotic resistance gene (ARG) profiles in non-human primates (NHPs) is crucial for evaluating their potential impact on human health and the environment. In this study, we performed metagenomic analysis of 203 primate fecal samples, including nine NHP species and humans, to comprehensively characterize their gut microbiota and ARGs. Our study reveals the prevailing phyla in primates as Firmicutes, Bacteroidetes, Euryarchaeota, and Proteobacteria. The captive NHPs exhibited higher ARG abundance compared to their wild counterparts, with tetracycline and beta-lactam resistance genes prevailing. Notably, ARG subtypes in Trachypithecus leucocephalus (T. leucocephalus) residing in karst limestone habitats displayed a more dispersed distribution compared to other species. Interestingly, ARG profiles of NHPs clustered based on geographic location and captivity status. Co-occurrence network analysis revealed intricate correlations between ARG subtypes and bacterial taxa. Procrustes analysis unveiled a significant correlation between ARGs and microbial phylogenetic community structure. Taxonomic composition analysis further highlighted differences in microbial abundance among NHPs and humans. Our study underscores the impact of lifestyle and geographical location on NHP gut microbiota and ARGs, providing essential insights into the potential risks posed by NHPs to antibiotic resistance dissemination. This comprehensive analysis enhances our understanding of the interplay between NHPs and the gut resistome, offering a critical reference for future research on antibiotic resistance and host-microbe interactions.

Metagenomic sequencing

Antibiotic resistance genes

Gut microbiota

Non-human primates

Network analysis

Antibiotic resistance is a significant global threat[1–5]. More than 100,000 tons of antibiotics are produced worldwide each year and are extensively used in aquaculture, healthcare, and clinical laboratories, leading to widespread drug resistance in many bacteria and the emergence of ARGs[6–11]. The World Bank has warned that if we do not take strong actions, antimicrobial resistance is expected to kill 10 million people by 2050, with accumulative economic losses reaching 100 trillion dollars[12, 13]. What’s worse, ARGs are widely distributed in sewage, soil, surface water, medical waste, and animal feces, and they can be transmitted to humans through agriculture products, water, or air in both captive and wild environments, which will increase the difficulty of dealing with antibiotic resistance[14, 15]. From the perspective of “One Health”, an effective way to accelerate harmless and resourceful treatment of ARGs is crucial[16].

NHPs are genetically and evolutional relatives of humans, providing a good model for studying the spread of ARGs[17]. At present, there are many kinds of NHPs lives in Guangxi, China, including T. leucocephalus, Trachypithecus francoisi (T. francoisi), Macaca multta (M. multta), etc. Studies on them have focused on genome adaptation, conservation, and behavior[18–22], but studies on the diversity and composition of ARGs in them have been lacking. Identification of ARGs in NHPs can be used as a proactive method for monitoring antibiotic resistance and as a control measure against infection[23–25]. A study showed that seasonal variations in the gut microbiota of free-range white-faced capuchins, but their ARGs composition remained unknown[26]. The composition and abundance of ARGs carried by the gut microbiota of cynomolgus monkeys were significantly affected by dietary changes[27]. Some antibiotics, such as beta-lactam, aminoglycoside, macrolides, and bacitracin in captive pandas were significantly higher than those in wild pandas[28]. Nevertheless, few studies have been conducted on the gut microbiota of multiple NHPs, and their gut microbiota and composition of the ARGs they carried are unknown. Therefore, there is an urgent need to fully understand the gut microbiota of NHPs and the composition of ARGs carried by NHPs, especially in the endangered primate T. leucocephalus.

Karst limestone is one of the typical geographical features of Guangxi, and it is the habitat of T. leucocephalus. Whether karst limestone has an impact on the gut microbiota and ARGs of primates, including T. leucocephalus, remains to be further investigated. Our previous findings showed that the gut microbiota of T. leucocephalus is mainly related to heat tolerance and cellulose degradation. Xiang et al.[29] found that the diversity, and composition of microbiota in karst water in southwest China were closely related to ARGs and antibiotic usage. C. Stange et al.[30] detected beta-lactams, aminoglycosides, tetracyclines, macrolides, and beta-lactams in karst lava springs in Germany, the frequency of ermB (42.1%), tet (C) (40.8%), sul2 (39.5%), and sul1 (36.8%) was higher. With increasing detection levels of ARGs detected in the environment, we speculate that the karst limestone in Guangxi has an impact on the gut microbiota and ARGs of local primates, and they play an important role in the development and human-to-human transmission of ARGs[31].

In this study, 150 fecal samples of 9 NHPs were collected in Guangxi. At the same time, metagenomic data of Chinese population was downloaded from the database. In this study, high-throughput sequencing and metagenomic analysis of ARGs in NHPs in Guangxi was conducted to explore the driving factors of lifestyle, geographical location and host species on ARGs, and to study the symbiotic network between ARG subtypes and related microbiota. This study will provide a comprehensive understanding of the distribution of ARGs in the NHPs gut microbiota in Guangxi and reveal their potential risks to human health.

2.1. Sample collection and public datasets acquisition

A total of 150 fecal samples were collected, including 27 T. leucocephalus (1 from Terrestrial Wildlife Course Research and Epidemic Diseases Monitor Center and 26 from Chongzuo White-headed Langur National Nature Reserve), 20 T. francoisi (6 from Wuzhou Langur Breeding and Research Center and 14 from Nanning Zoon), 5 Papio, 5 Hylobatidae, 9 Trachypithecus cristatus (T. cristatus), 23 Lemur catta (from Nanning Zoon), 8 Nycticebus pygmaeus (N. pygmaeus), 19 Nycticebus coucang (N. coucang) (from Terrestrial Wildlife Medical-aid Monitoring Epidemic Diseases Research Center), and 34 M. mulatta samples (from Hezhou Gupo Mountain) (Table S1, Figure S1). All the samples were collected during the spring of 2018 to mitigate potential seasonal variations in microbial content.

To obtain wild fecal samples, we strategically placed food on the ground and concealed ourselves in vegetation. Upon defecation, we used latex gloves to select fresh samples based on visual cues related to sample integrity and color. For captive specimens, healthy NHPs were randomly chosen, and their fecal samples were promptly collected in an aseptic manner post-defecation. Using a spoon, feces from the center portion were collected and deposited into sterile containers, which were then transported to the laboratory in a low-temperature box accompanied by dry ice. Subsequently, samples were immediately frozen and stored at − 80°C for subsequent processing. DNA was extracted within 48 hour following samples delivery. The research adhered to protocols established by the China Wildlife Conservation Association and the legal requirements of China.

For comparative analysis, metagenomic data (PRJEB21612) of 53 health humans in China were downloaded from the European Nucleotide Archive Database (https://www.ebi.ac.uk/ena/browser/view/). These datasets were incorporated into our study to enhance the breadth of our analysis.

2.2. DNA extraction and sequencing

DNA extraction was conducted using the TIANamp Stool DNA Kit (Germany) according to the manufacturer’s instructions. DNA concentration and purity were assessed using the Qubit® 2.0 Flurometer (Life Technologies, CA, USA) and NanoPhotometer (IMPLEN, CA, USA), respectively. DNA integrity was evaluated through 1% agarose gel electrophoresis. Metagenomic DNA was fragmented using ultrasonic Processor, yielding an average size of approximately 350 bp. A paired-end library was generated using the NEBNext® Ultra™ DNA Library Prep Kit for Illumina (NEB, USA), with adapters ligation performed on blunt-end fragments containing full sequencing primer hybridization sites. High-throughput paired-end sequencing was performed on the Illumina HiSeq X-ten platform utilizing the index PE150 sequencing strategy. Subsequently, sequences with low-quality bases (Q ≤ 19) at the terminal regions were trimmed to ensure high-quality data. Reads containing substantial "N" (≥ 5% of the read), adapter contamination, or those falling below a specified length (> 100-bp) were filtered out. Following these steps, a set of clean data possessing high quality were obtained for further analysis.

2.3. Profiling of ARGs and microorganism composition

We quantified the profiles of ARGs in the 203 metagenomic datasets using the ARGs-OAP pipeline (version 2.0) (http://smile.hku.hk/SARGs/) [13, 32–34]. The datasets were aligned to a structured ARG (SARG) database (http://smile.hku.hk/SARGs) following established protocols outlined in prior publications. Specifically, we combined the metadata file of samples with the clean reads files in a single directory within our local Linux system. The initial Perl script was executed to identify ARG-like sequences and 16S ribosomal RNA (rRNA) gene sequences, with parameters set at an alignment length of 25 amino acids, 80% sequence similarity, and an e-value of 1 × 10^− 7[32, 34]. Subsequently, a candidate ARG sequences file and an online metadata file were matched against the integrated structure to align with SARG reference databases, encompassing 24 types and 1244 subtypes of ARGs, using the second Perl script within our local Linux system. ARG abundances are normalized with respect to total 16S reads (unit: copies of ARG per copy of 16S rRNA), the number of cells (unit: ARG copies per cell), and ppm (reads carrying ARGs per million reads), a unit normalized by dividing the total reads carrying ARGs by a million read counts[35]. This process yielded the ARG profile abundance at different levels, type, subtype level, and gene (reference sequence).

For taxonomic composition analysis, metagenomes data underwent taxonomic classification and relative abundance quantification (in terms of the number of cells, not just fraction of reads) through HUMAnN2 pipeline (v0.11.1) integrated with the MetaPhlAn database of marker genes mpa_v20_m200[36, 37]. Differential abundances of taxa were identified using LEfSe[38].

2.4. Procrustes analysis

We conducted correlation analyses between ARGs-ARGs and ARGs-genus taxa using ‘vegan’ package in R (v.4.3.1) [39]. For comparisons of M2 values with 999 permutations, Spearman's correlation coefficients were calculated using the protest function[13]. The M2 value, representing the sum of the squared distances between matched sample pairs, was computed. To obtain the P-value, the original M2 value was compared to the simulated distribution of M2 values. When the P-value was below 0.05, a significant correlation between antibiotic resistance and bacterial species was established[32].

2.5. Statistical and network analysis

Non-metric multidimensional scaling (NMDS) and principal component analysis (PCA) were conducted using the abundance correlation matrix of ARG subtypes. A heat map illustrating the abundance of ARG types and subtypes was generated using Hemi1.0. Column diagram depicting the abundance of ARG types were drawn using GraphPad prism 9.0. Petal figure showing specific genes was drawn using Venn (http://www.ehbio.com/test/venn/#/)[40]. Simpson index, Shannon index and Chao1 index (measuring within-sample diversity) were calculated based on ARG relative abundance profiles using the PAST, and visualization was accomplished using origin 2022.

To reduce false-positive results, the False Discovery Rate (FDR) method was applied to adjust P-value using JMP 16[41]. Network analysis was performed using Gephi (0.10.1) based on significant Spearman's correlations (r > 0.6, p < 0.01)[42–45].

3.1. Microbial community annotation of different NHPs.

In all primate species, bacterial overwhelmingly dominated, contributing a greater degree of phylogenetic diversity compared to archaea, eukaryotes, and viruses (Figure S2). At the phylum level, all samples were allocated to 17 distinct phyla, encompassing 14 for NHPs and 12 for humans. Striking similarities emerged between NHPs and humans in terms of their bacterial composition patterns. Notably, within Homo sapiens, T. leucocephalus, Hylobatidae, Papio, M. multta, T. francoisi, N. pygmaeus, T. cristatus, and N. coucang, the predominant phylum was Firmicutes (51%, 47%, 95%, 69%, 71%, 53%, 77%, 81% and 65%, respectively), while in Lemur catta, the most abundant phylum was Spirochaetes (58.0%). Broadly, Firmicutes, Bacteroidetes, Euryarchaeota, and Proteobacteria constituted the four prevailing phyla across all species, cumulatively accounting for approximately 85.83% of the total bacterial community (Fig. 1A).

At the genus level, 228 genera were identified across all samples, showcasing diverse representative genera in each species. A substantial proportion of the Subdoligranulum genus were identified in Presbytis (T. leucocephalus, T. francoisi, T. cristatus) (56.53%, 31.03% and 42.36 respectively), while the predominant genus Megamonas was observed in Nycticebus (N. pygmaeus, N. coucang) (47.93% and 27.14% respectively). Notably, humans exhibited a uniquely enriched genus, Bacteroides (22%), whereas Papio harbored the most abundant genus Megasphaera (17%). In Hylobatidae, the leading genus was Faecalibacterium (33%), whereas in Lemur catta, the genus Treponema was particularly enriched (38%) (Fig. 1B). Furthermore, a total of 85 species-specific genera were identified, each characterized by relatively lower abundances. For instance, Papio had 20 such genera, T. leucocephalus had 19, Homo sapiens had 14, N. coucang had 12, Lemur catta had 9, N. pygmaeus had 4, T. cristatus had 3, M. multta had 2, and both Hylobatidae and T. francoisi had 1 each. Additionally, 18 genera were detected as shared across all species, including Subdoligranulum, Bacteroides, Phascolarctobacterium and Methanobrevibacter, which were included in the top 50 most abundant genera (Figure S3).

3.2. Distinct gut bacterial communities across NHPs.

The phylogenetic relationships among the ten primate species are depicted in Fig. 2A. Utilizing a genus-level relative abundance profile, a principal coordination analysis (PCoA) plot elucidated that axis 1 (PCoA1) accounted for 26.77% of the variance, while axis 2 (PCoA2) accounted for 16.27%. Homo sapiens, N. coucang and N. pygmaeus samples exhibited clustering tendencies, distinct from the remaining samples (Fig. 2B). Intriguingly, species sharing genetic host backgrounds demonstrated resemblant gut microbial structures, diverging from prior findings[41]. Employing LEfSe, we compared the relative abundance of gut bacteria between humans and seven other NHPs species (Figure S4). Specifically, 2 phyla (Actinobacteria and Bacteroidetes), 2 classes (Actinobacteria and Bacteroidia), 2 orders (Bacteroidales and Bifidobacteriales), 4 families (Bacteroidaceae, Bifidobacteriaceae, Eubacteriaceae and Rikenellaceae), 7 genera (Alistipes, Bacteroides, Bifidobacterium, Blautia, Eubacterium, Faecalibacterium and Ruminococcus), and 10 species (Ruminococcus_bromii, Roseburia_inulinivorans, Roseburia_intestinalis, Faecalibacterium_prausnitzii, Eubacterium_rectale, Eubacterium_eligens, Bifidobacterium_longum, Bacteroides_vulgatus, Bacteroides_stercoris) exhibited significantly higher abundance in homo sapiens (P < 0.01 and LDA > 4) (Fig. 2D).

In total, 107 opportunistic pathogens were detected across all species through comparison with a previously published pathogen list[46]. Among these, 10 pathogens were shared across all ten primates (Fig. 2C). Escherichia coli and Bacteroides eggerthii prevailed most frequently, with Homo sapiens (59), N. coucang (55) and N. pygmaeus (45) hosting the broadest spectrum of pathogens. Clinically pertinent pathogens were identified in NHPs, including Clostridium perfringens, Klebsiella pneumoniae, and Clostridium difficile, which were also found in the migratory birds.

3.3. The ARGs and resistant phenotype among NHPs.

A total of 18 distinct ARG types were identified across various primate populations. Tetracycline resistance genes comprised the largest proportion at 40% of all detected ARGs, followed by MLS resistance genes (21%), beta-lactam resistance genes (11%), and multidrug resistance genes (8%) ((Figure S5 and Figure S6). The ARGs for Tetracycline, bacitracin, MLS, multidrug, and vancomycin were consistently detected in all samples (Fig. 3A). The total abundance of ARGs in the different species exhibited the following ranges: Homo sapiens (0.384–1.597), T. leucocephalus (9.587×10^− 6–4.480×10^− 1), M. multta (1.023×10^− 5–1.040×10^− 1), N. coucang (1.438×10^− 5–3.24×10^− 1), N. pygmaeus (1.007×10^− 5–2.890×10^− 1), Lemur catta (6.671×10^− 6–1.250×10^− 1), T. francoisi (1.004×10^− 5–2.670×10^− 1), T. cristatus (1.057×10^− 5–1.880×10^− 1), Hylobatidae (2.018×10^− 5–2.410×10^− 1), and Papio (4.513×10^− 5–2.20×10^− 1) (Fig. 3C).

More specifically, within M. multta, the dominant resistance genes were tetracycline, MLS, vancomycin, and aminoglycoside, accounting for 65.71%, 9.00%, 5.23%, and 4.61% respectively (Figure S7 A). For T. leucocephalus, the foremost resistance genes were multidrug, tetracycline, MLS, and bacitracin, constituting 46.47%, 16.26%, 9.96%, and 8.05%, respectively (Figure S6 B). In captive N. coucang and N. pygmaeus, 16 ARGs were observed, with prominent types being tetracycline (49.13%, 47.39%), MLS (18.07%, 13.96%), and beta-lactam (9.07%, 7.02%) (Fig. 7C, D). Tetracycline, aminoglycoside, MLS, and vancomycin accounted for 42.20%, 12.01%, 6.69%, and 10.82% respectively in T. francoisi, and 62.90%, 12.15%, 6.95%, and 5.03% respectively in T. cristatus (Figure S6 G, H). Among Hylobatidae and Papio, which exhibited 18 ARGs, tetracycline (46.19%), beta-lactam (32.02%), MLS (5.98%), and bacitracin (5.03%) were predominant in Hylobatidae (Figure S6 F). In Papio, tetracycline (39.77%), multidrug (20.41%), MLS (11.95%), and beta-lactam (5.92%) dominated (Figure S6 I, J). For Lemur catta, 18 ARGs were evident, with tetracycline (36.51%), beta-lactam (24.3%), MLS (18.55%), and chloramphenicol (6.64%) as primary resistance genes (Figure S7 E). Notably, all ARGs annotated in human were also present in the NHPs (Figure S7 K).

A total of 580 ARGs were identified across all samples, with relative abundances spanning from 3.01x10^− 6 (TEM-201) to 10.07 (tetQ). The heatmap showcased the top 100 most abundant ARGs, comprising 34 classified as multidrug resistance, 16 as tetracycline, 11 as MLS, and 8 as aminoglycoside (Fig. 2B). vanR and macB were universally shared across all samples, indicating their ubiquitous presence among primates. Tetracycline-related ARGs, including tet35, tet36, tet37, tet39, tet40, tet44, tetA, tetL, tetM, tetO, tetP, tetQ, tetracycline_resistance_protein, tetW, tetX1, and tetX2, exhibited high abundances in captive primates. Furthermore, a total 93 ARGs were shared among the ten primate species, with Papio harboring the most diverse ARGs (355) (Figure S8).

3.4 Comparative analysis of resistome between NPHs and humans

Previous studies have underscored the close resemblance between the gut microbiomes of NHPs and humans, suggesting the potential utility of NHPs in biomedical research as model[38, 47]. Building on this premise, we conducted a comparison of antibiotic resistome composition and diversity indices between nine NHPs and humans adhering to a non-westernized diet. Our analysis revealed the existence of 16 ARGs shared between both humans and NHPs. Moreover, rifamycin resistance genes were exclusively identified in NHPs, albeit at a limited frequency. Among these genes, tetracycline stood out as the most prevalent and dominant ARGs in metagenomic resistome profiles of both humans and NHPs, followed by MLS, beta-lactam, and multidrug (Fig. 4A).

Diversity index analysis demonstrated that human exhibited relatively higher Simpson, Shannon, and Chao1 indices, compared to NHPs. Notably, captive Papio displayed the highest diversity indexes. PCA plots based on ARGs relative abundance indicated that nearly all NHPs samples (particularly those from T. francoisi, T. cristatus, N. pygmaeus, and M. mulatta) were closely situated to human samples (Fig. 4D, Figure S9 B-D, F). However, some T. leucocephalus and Lemur catta samples diverged from human samples in the PCA plot (Figure S10 A and E). A considerable abundance of ARGs was identified in samples from the Wuzhou Langur Breeding and Research Center and Nanning Zoon’s T. francoisi (Figure S9). consequently, it can be inferred that lifestyles and geography exert an influence on similarities and abundance of ARGs among the NHPs.

Thus, NHPs resistomes represent a potentially crucial reservoir for studying antibiotic resistance within human gut microbiota. This warrants attention as captive NHPs frequently interact with humans, raising concerns for potential public health implications. Furthermore, it is noteworthy that both wild T. leucocephalus and M. mulatta harbor numerous ARGs that align with those found in human beings.

3.5. Co-occurrence network analysis of ARGs in NHPs.

Procrustes analysis was employed to elucidate correlations between ARG subtypes and bacteria taxa at the genus level. Deviation square M2 yielded a value of 0.9374, with a corresponding P value 0.001 for ARGs and ARGs (Fig. 5A). Additionally, the deviation square M2 for ARGs and microbiota amounted to 0.9708, with P value of 0.001 (Fig. 5B).

To unravel the co-occurrence patterns of ARGs within the ten NHPs, a co-occurrence network was constructed based on robust and statistically significant correlations among different NPHs found in individuals (|ρ|>0.6, P < 0.01). The co-occurrence network diagrams for ARGs contained 12 antibiotic types encompassing 79 resistance genes (Fig. 5C). Positive correlations emerged among similar ARGs indicative of the shared selective pressure imposed by specific antibiotics. For example, tetracycline ARGs tetO, tetM, tet32, tetX2, tetX and tetX3 exhibited such correlations. Similar patterns were observed for tetracycline-related ARGs vanR and vanS. Only one negative correlation was observed between different ARGs, such as bacA (bacitracin ARG) and vatB (MLS ARG). This could be attributed to competition among host bacteria or inhibition-related dynamics[48].

Illustrated in Fig. 5D, the network comprised 61 nodes and 84 edges (with a modularity index of 0.602), including Bifidobacterium, Clostridium, Subdoligranulum, and Ruminococcaceae_noname. Furthermore, 38 species levels were linked to 10 resistance gene types, encompassing tetracycline, vancomycin, aminoglycoside, bacitracin, multidrug, MLS, beta-lactamase, beta-lactam, trimethoprim and chloramphenicol. Notably, Treponema exhibited a positive correlation with vatB (MLS ARG), indicating that its potential role as a host for ARGs. Concurrently, certain species displayed negative correlation with ARGs, exemplified by the Spirochaetes phylum, Treponema, and bacA (bacitracin ARG). Thus, intricate interaction between various ARGs and microbes were evident, underscoring the intimate connection between NHP gut microbes and ARGs. Besides, based on these results we can tract the hosts for specific ARG subtypes[49].

NHPs, situated proximately to humans within the animal kingdom, share a substantial degree of genetic and physiological resemblance[27]. Exploring and contrasting the gut microbiota and antibiotic resistome within captive and wild primate populations can yield valuable insights into human interaction and captivity on NHPs[47]. Consequently, a comprehensive understanding of the gut microbiota and ARGs across diverse NHPs and their distinctions from humans is of notable import. Nevertheless, the existing research in this domain remains limited. In this study, we observed a greater similarity in gut microbial composition particularly between captive NHPs (N. coucang, N. pygmaeus) and humans. Our investigation revealed 107 opportunistic pathogens and 580 distinct ARGs, unveiling active correlation between the ARGs to ARGs, as well as ARGs and microbiota at genus level, across ten primate species.

We established that Firmicutes, Bacteroidetes, Euryarchaeota and Proteobacteria emerged as the dominant phyla in NHPs. Notably, the gut microbiome clustering exhibited species affiliations within specific geographic region. Human gut microbial composition exhibited a closer affinity to that of N. coucang and N. pygmaeus compared to other NHPs. As such, we infer that NHPs captivity may foster similarities in gut microbiomes akin to humans. However, discernible differences persist in the gut microbiota of captive NHPs and modern humans following diverse diets. This observation aligns with our identification of 10 species (Ruminococcus_bromii, Roseburia_inulinivorans, Roseburia_intestinalis, Faecalibacterium_prausnitzii, Eubacterium_rectale, Eubacterium_eligens, Bifidobacterium_longum, Bacteroides_vulgatus, Bacteroides_stercoris) that exhibited significantly greater abundance within the human gut, as determined through LEfSe analysis.

Furthermore, our investigation unveiled a notable trend of tetracycline and beta-lactam resistance among captive primates. Remarkably, tetracycline resistance was also discerned in T. leucocephalus, albert at a proportion lower than that in M. mulatta. Tetracycline and beta-lactam, recognized as broad-spectrum antibiotics, find primary use as first-line treatments for non-communicable diseases resulting from bacterial infections[50, 51]. Their antibacterial spectrum encompasses Gram-positive and Gram-negative bacteria, chlamydia, mycoplasma, rickettsia, spirochetes, mycobacteria, and some protozoa. Consequently, tetracycline and beta-lactam are extensively employed in veterinary clinics for dogs, cats, cattle, sheep, pigs, turkey, and chickens. This widespread application potentially underpins the elevated prevalence of tetracycline and beta-lactam resistance among captive NHPs. Among the 17 identified tetracycline-related ARGs, tetX stands out as an enzyme modification gene, capable of inactivating antibiotic enzymes through enzymatic hydrolysis. Notable, tetQ, tetW, tetO, tet32, and tet44 genes, which are also prevalent resistance genes in pigs and cattle, manifest substantial abundance[52]. For instance, tet32, initially identified in Clostridium difficile-associated human colonic anaerobe K10, finds wide distribution in bovine rumen and pig feces[53]. The identification of numerous tetracycline-related ARGs underscores the necessity for heightened resistance surveillance and judicious antibiotic usage. Strategies such as alternative antibiotic selection or modifications in the living and dietary conditions of NHPs warrant consideration as control measures.

Co-occurrence network analysis offers a high-resolution depiction of correlations between ARGs and microbiota; however, it does not elucidate definitive relationships. Nonrandom symbiosis models can be used to elucidate potential hosts of ARGs within animal feces and environment samples. Essentially, this nonrandom symbiosis pattern suggests that specific microbiota may harbor distinct ARGs. In our study, we postulated that when the abundance trends of ARGs align significantly with coexisting microbiota across different species, a nonrandom co-occurrence pattern potentially exists between ARGs and microbiota. This led us to infer that Clostridium (with 11 ARG subtypes), Bifidobacterium (with 10 ARG subtypes), and T4likevirus (with 7 ARG subtypes) might emerge as primary potential hosts for tetracycline and MLS. Given that these ARG-carrying bacteria may pose substantial risks to animal health, and potentially transmit to humans, it becomes imperative to implement stringent control measures to mitigate these concerns and curtail the risk of therapeutic failure. Several potential ARGs hosts were confirmed through experimental validation, such as mdtO in Escherichia, underscoring the value of co-occurrence network analysis as a robust analytical tool capable of unearthing fresh insights into ARGs and their potential hosts within intricate environmental systems. Nevertheless, these potential symbiosis relationships necessitate further exploration through alternate methods, such as metagenomic assembly and microbiological cultivation.

Compared to human, we observed a resemblance in human gut ARGs to those present in captive NHPs, with numerous overlapping ARG detected in both populations. An investigation involving chimpanzees, gorillas, and humans with different lifestyles reported a segregation of ARGs in wild chimpanzees and gorillas, while captive individuals exhibited a clustering closer to humans[47]. Our result further revealed an elevated ARGs abundance in captive primates compared to their wild counterparts. Concurrently, primate species within the same geographic region, such as N. coucang and N. pygmaeus, exhibited similarities in their ARG profiles, whereas significant differences in ARGs abundance were evident in T. francoisi across distinct regions. Remarkably, the detected ARG in humans were also observed in NHPs, underscoring the notable impact of lifestyle and geographic location on gut microbiota-associated ARGs. Moreover, the dispersion of ARGs in T. leucocephalus, with certain samples approximating the human pattern, potentially signifies interactions with humans, possibly through contact with agricultural products. In totality, our findings suggest that lifestyle and geographical location exert substantial influence on NHPs ARGs dynamics.

Recognizing areas warranting refinement, we acknowledge the necessity for multivariate analysis to unveil the distinct contributions of lifestyle and geographic factors to the variability of ARG profiles. While our study presents a pioneering exploration of the ARG composition in T. leucocephalus and its comparison to humans, the inclusion of only one captive T. leucocephalus precludes a comprehensive analysis of the specific impact of the karst environment on ARGs. In our analyses, spearman correlation elucidated associations between ARG and genus-level microbiota. However, it is crucial to note that numerous ARGs can transfer among bacteria through self-proliferation and horizontal gene transfer in real-world settings. Therefore, certain genes, particularly those widely distributed across bacterial species, might not be solely discerned through correlation analysis[54]. Consequently, further investigations are warranted to deepen our insights in subsequent studies. Regardless, our study unveils the gut microbiota and resistome of different primates by metagenomic sequencing, providing essential evidence for NHPs and humans to combat the threat of antibiotic resistance.

Through comprehensive metagenomic analysis of ten primates sequencing datasets, we have established an extensive and diverse repository of primate gut microbiota and ARGs. Our study underscores the potential influence of factors like captivity, geographic location, and species on the composition of gut microbiota and ARGs in NHPs. Notably, captive NHPs exhibited elevated relative abundance of ARGs compared to their wild counterparts. Moreover, our co-occurrence network analysis revealed discernible correlation between ARGs and distinct microbiota components. This research enriches our comprehension of the parallels and distinctions in gut microbiota and ARGs profiles between NHPs and humans, offering a pivotal reference for future investigation in this domain.

Acknowledgements

We are grateful to the Guangxi Zhuang Autonomous Region Terrestrial Wildlife Course Research and Epidemic Diseases Monitor Center, the Chongzuo White-headed Langur National Nature Reserve, the Nanning Zoo, the Wuzhou Langur Breeding and Research Center, and the Hezhou Gupo Mountain for providing a valuable source of samples and for assisting in sample collection. We acknowledge the providers who submitted the metagenomics-seq data to the public databases. We also acknowledge Lei Yu for technical support during the conduct of the experimental analysis. We sincerely thank the editors and anonymous reviewers for their insightful and constructive feedback that helped to improve this manuscript.

Authors’ contributions

Yanling Hu designed and supervised the project. Hongli Huang and Xianwu Pang performed the bioinformatics analyses and generated the figures and tables; Tengcheng Que, Panyu Chen, Shousheng Li, Aiqiong Wu, and Meihong He collected and prepared samples; Tengcheng Que, Hong Qiu and Panyu Chen supported and performed laboratory experiments and sequencing. Hongli Huang and Xianwu Pang wrote and revised the manuscript. All authors read and approved the final manuscript.

Funding

This study was funded by Guangxi Key Research and Development Program (No: Guike AB20059002).

Availability of data and materials

The data of 27 T. leucocephalus and 53 human individuals are available at the Genome Sequence Archive (Genomics, Proteomics, and Bioinformatics 2017) in BIG Data Center (Nucleic Acids Res 2019), Beijing Institute of Genomics (BIG), Chinese Academy of Sciences (https://bigd.big.ac.cn/gsa) and European Nucleotide Archive Database (https://www.ebi.ac.uk/ena/browser/view/) with accession number CRA002281 and PRJEB21612. All data generated or analyzed during this study are included in this published article [and its supplementary information files].

Ethics approval and consent to participate

The animal study was reviewed and approved by the Ethics Committee of the Guangxi Chongzuo White-Headed Langur National Nature Reserve Management Center.

Competing interests

The authors declare that they have no competing of interests.

Bolan GA, Sparling PF, Wasserheit JN. The Emerging Threat of Untreatable Gonococcal Infection. N Engl J Med. 2012;366(6):485–7.
Brinkac L, et al. The Threat of Antimicrobial Resistance on the Human Microbiome. Microb Ecol. 2017;74(4):1001–8.
Mustafa GR et al. Metagenomic analysis revealed a wide distribution of antibiotic resistance genes and biosynthesis of antibiotics in the gut of giant pandas. BMC Microbiol, 2021. 21(1).
Gil-Gil T, et al. Antibiotic resistance: Time of synthesis in a post-genomic age. Comput Struct Biotechnol J. 2021;19:3110–24.
Qiao M, et al. Review of antibiotic resistance in China and its environment. Environ Int. 2018;110:160–72.
Van Boeckel TP, et al. Global trends in antimicrobial use in food animals. Proc Natl Acad Sci USA. 2015;112(18):5649–54.
Zhang QQ, et al. Comprehensive Evaluation of Antibiotics Emission and Fate in the River Basins of China: Source Analysis, Multimedia Modeling, and Linkage to Bacterial Resistance. Volume 49. Environmental Science & Technology; 2015. pp. 6772–82. 11.
Xu Y, et al. Effect of the selective pressure of sub-lethal level of heavy metals on the fate and distribution of ARGs in the catchment scale. Environ Pollut. 2017;220:900–8.
Sun YM, et al. Inorganic and organic fertilizers application enhanced antibiotic resistome in greenhouse soils growing vegetables. Ecotoxicol Environ Saf. 2019;179:24–30.
Wang CL et al. Characterization of the Pig Gut Microbiome and Antibiotic Resistome in Industrialized Feedlots in China. Msystems, 2019. 4(6).
Qiu QW et al. Metagenomic Analysis Reveals the Distribution of Antibiotic Resistance Genes in a Large-Scale Population of Healthy Individuals and Patients With Varied Diseases. Front Mol Biosci, 2020. 7.
He LY, et al. Microbial diversity and antibiotic resistome in swine farm environments. Sci Total Environ. 2019;685:197–207.
Qian X et al. Metagenomic analysis reveals the shared and distinct features of the soil resistome across tundra, temperate prairie, and tropical ecosystems. Microbiome, 2021. 9(1).
Chi T et al. Characteristics of the antibiotic resistance genes in the soil of medical waste disposal sites. Sci Total Environ, 2020. 730.
Yang JT et al. Metagenomic Insights into Chicken Gut Antibiotic Resistomes and Microbiomes. Microbiol Spectr, 2022. 10(2).
Su ZG et al. Industrial effluents boosted antibiotic resistome risk in coastal environments. Environ Int, 2023. 171.
Yildirim S et al. Characterization of the Fecal Microbiome from Non-Human Wild Primates Reveals Species Specific Microbial Communities. PLoS ONE, 2010. 5(11).
Huang CM, Li YB. How does the white-headed langur (Trachypithecus leucocephalus) adapt locomotor behavior to its unique limestone hill habitat? Primates. 2005;46(4):261–7.
Li YB, Huang XH, Huang ZH. Behavioral adjustments and support use of Francois' langur in limestone habitat in Fusui, China: Implications for behavioral thermoregulation. Ecol Evol. 2020;10(11):4956–67.
Liu ZJ, et al. Genomic Mechanisms of Physiological and Morphological Adaptations of Limestone Langurs to Karst Habitats. Mol Biol Evol. 2020;37(4):952–68.
Fan Y, Pedersen O. Gut microbiota in human metabolic health and disease. Nat Rev Microbiol. 2021;19(1):55–71.
Que TC et al. The reference genome and transcriptome of the limestone langur, Trachypithecus leucocephalus, reveal expansion of genes related to alkali tolerance. BMC Biol, 2021. 19(1).
Roca I et al. The global threat of antimicrobial resistance: science for intervention. New Microbes and New Infections, 2015. 6: p. 22–29.
Schuchat A, et al. Active bacterial core surveillance of the emerging infections program network. Emerg Infect Dis. 2001;7(1):92–9.
Liu YY, et al. Emergence of plasmid-mediated colistin resistance mechanism MCR-1 in animals and human beings in China: a microbiological and molecular biological study. Lancet Infect Dis. 2016;16(2):161–8.
Orkin JD, et al. Seasonality of the gut microbiota of free-ranging white-faced capuchins in a tropical dry forest. ISME J. 2019;13(1):183–96.
Yan YY et al. Metagenomic and network analysis revealed wide distribution of antibiotic resistance genes in monkey gut microbiota. Microbiol Res, 2022. 254.
Guo W et al. Comparative Study of Gut Microbiota in Wild and Captive Giant Pandas (Ailuropoda melanoleuca). Genes, 2019. 10(10).
Xiang SZ et al. Response of microbial communities of karst river water to antibiotics and microbial source tracking for antibiotics. Sci Total Environ, 2020. 706.
Stange C, Tiehm A. Occurrence of antibiotic resistance genes and microbial source tracking markers in the water of a karst spring in Germany. Sci Total Environ, 2020. 742.
Allen HK, et al. Call of the wild: antibiotic resistance genes in natural environments. Nat Rev Microbiol. 2010;8(4):251–9.
Feng J, et al. Antibiotic resistome in a large-scale healthy human gut microbiota deciphered by metagenomic and network analyses. Environ Microbiol. 2018;20(1):355–68.
Yang YY, et al. Metagenomic insights into the abundance and composition of resistance genes in aquatic environments: Influence of stratification and geography. Environ Int. 2019;127:371–80.
Yin XL, et al. ARGs-OAP v2.0 with an expanded SARG database and Hidden Markov Models for enhancement characterization and quantification of antibiotic resistance genes in environmental metagenomes. Bioinformatics. 2018;34(13):2263–70.
Zheng D, et al. Global biogeography and projection of soil antibiotic resistance genes. Sci Adv. 2022;8(46):eabq8015.
Gacesa R, et al. Environmental factors shaping the gut microbiome in a Dutch population. Nature. 2022;604(7907):732–9.
Li B, et al. Metagenomic and network analysis reveal wide distribution and co-occurrence of environmental antibiotic resistance genes. ISME J. 2015;9(11):2490–502.
Yang S, et al. The gut microbiome and antibiotic resistome of chronic diarrhea rhesus macaques (Macaca mulatta) and its similarity to the human gut microbiome. Microbiome. 2022;10(1):29.
Munk P et al. Abundance and diversity of the faecal resistome in slaughter pigs and broilers in nine European countries (vol 3, pg 898, 2018). Nature Microbiology, 2018. 3(10): p. 1186–1186.
Chen T, et al. EVenn: Easy to create repeatable and editable Venn diagrams and Venn networks online. J Genet Genomics. 2021;48(9):863–6.
Cao J, et al. Metagenomic analysis reveals the microbiome and resistome in migratory birds. Microbiome. 2020;8(1):26.
Li B, et al. Metagenomic and network analysis reveal wide distribution and co-occurrence of environmental antibiotic resistance genes. ISME J. 2015;9(11):2490–502.
Liu JW, et al. Spatiotemporal dynamics of the archaeal community in coastal sediments: assembly process and co-occurrence relationship. ISME J. 2020;14(6):1463–78.
Liu YX, et al. A practical guide to amplicon and metagenomic analysis of microbiome data. Protein Cell. 2021;12(5):315–30.
Yang FX et al. Cow manure simultaneously reshaped antibiotic and metal resistome in the earthworm gut tract by metagenomic analysis. Sci Total Environ, 2023. 856.
Li B, et al. Profile and Fate of Bacterial Pathogens in Sewage Treatment Plants Revealed by High-Throughput Metagenomic Approach. Environ Sci Technol. 2015;49(17):10492–502.
Campbell TP, et al. The microbiome and resistome of chimpanzees, gorillas, and humans across host lifestyle and geography. ISME J. 2020;14(6):1584–99.
Yan Y, et al. Metagenomic and network analysis revealed wide distribution of antibiotic resistance genes in monkey gut microbiota. Microbiol Res. 2022;254:126895.
Li ZH, et al. Metagenome sequencing to unveil the occurrence and distribution of antibiotic resistome and in a wastewater treatment plant. Environmental Technology; 2023.
Hu T, et al. Geographic pattern of antibiotic resistance genes in the metagenomes of the giant panda. Microb Biotechnol. 2021;14(1):186–97.
Fair RJ, Tor Y. Antibiotics and Bacterial Resistance in the 21st Century. Perspect Med Chem, 2014. 6.
Lim SK et al. Antibiotic resistomes discovered in the gut microbiomes of Korean swine and cattle. Gigascience, 2020. 9(5).
Cao J et al. Metagenomic analysis reveals the microbiome and resistome in migratory birds. Microbiome, 2020. 8(1).
Liu JX et al. The fecal resistome of dairy cattle is associated with diet during nursing. Nat Commun, 2019. 10.

No competing interests reported.

SupplementaryInformation.docx

Download PDF

Version 1

posted

You are reading this latest preprint version

Antibiotic Resistance Profiles of Gut Microbiota Across Various Primate Species in Guangxi

Status:

Version 1

Abstract

Figures

1. Introduction

2. Methods

2.1. Sample collection and public datasets acquisition

2.2. DNA extraction and sequencing

2.3. Profiling of ARGs and microorganism composition

2.4. Procrustes analysis

2.5. Statistical and network analysis

3. Results

3.1. Microbial community annotation of different NHPs.

3.2. Distinct gut bacterial communities across NHPs.

3.3. The ARGs and resistant phenotype among NHPs.

3.4 Comparative analysis of resistome between NPHs and humans

3.5. Co-occurrence network analysis of ARGs in NHPs.

4. Discussion

5. Conclusion

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1