Skip to main content

Wildlife gut microbiomes of sympatric generalist species respond differently to anthropogenic landscape disturbances

Abstract

Background

Human encroachment into nature and the accompanying environmental changes are a big concern for wildlife biodiversity and health. While changes on the macroecological scale, i.e. species community and abundance pattern, are well documented, impacts on the microecological scale, such as the host’s microbial community, remain understudied. Particularly, it is unclear if impacts of anthropogenic landscape modification on wildlife gut microbiomes are species-specific. Of special interest are sympatric, generalist species, assumed to be more resilient to environmental changes and which often are well-known pathogen reservoirs and drivers of spill-over events. Here, we analyzed the gut microbiome of three such sympatric, generalist species, one rodent (Proechimys semispinosus) and two marsupials (Didelphis marsupialis and Philander opossum), captured in 28 study sites in four different landscapes in Panama characterized by different degrees of anthropogenic disturbance.

Results

Our results show species-specific gut microbial responses to the same landscape disturbances. The gut microbiome of P. semispinosus was less diverse and more heterogeneous in landscapes with close contact with humans, where it contained bacterial taxa associated with humans, their domesticated animals, and potential pathogens. The gut microbiome of D. marsupialis showed similar patterns, but only in the most disturbed landscape. P. opossum, in contrast, showed little gut microbial changes, however, this species’ absence in the most fragmented landscapes indicates its sensitivity to long-term isolation.

Conclusion

These results demonstrate that wildlife gut microbiomes even in generalist species with a large ecological plasticity are impacted by human encroachment into nature, but differ in resilience which can have critical implications on conservation efforts and One Health strategies.

Introduction

In the Anthropocene, a plethora of challenges have forced wildlife to cope with environmental changes, such as habitat fragmentation and modification [1]. While species that are highly susceptible to changes in their environment are likely to disappear, species that are more resilient to environmental change increase in numbers and will probably become more dominant in local species assemblages [2]. The breaking apart of natural habitat into patches can have negative effects on the abundance of species and their health [3, 4]. These so-called generalist species often thrive in or at edges of altered landscapes and have contact with the matrix surrounding the habitat fragments, often leading to closer contact to humans, livestock and domesticated animals [5]. As a consequence of inhabiting a wide range of habitats, generalist species are often important pathogen reservoirs and act as vectors of zoonotic diseases, despite appearing phenotypically healthy [4]. Increasing contact rates with humans and domesticated animals, as well as the exploitation of wildlife for food or medicine, amplify the risk of spillover events in both directions [6]. However, assessing wildlife health of generalist species coping with anthropogenic disturbances can be challenging because of the problems posed in macroecological studies in assessing the internal (health) states in wild animals. Even common fitness indicators, such as body mass can point into wrong directions if land use change requires dietary shifts of the surviving, apparently adaptable resilient species [7].

The gut microbiome, the assemblage of microbes and their genomes inhabiting the gut of a host, is an integral part of an animal’s well-being and is considered as a good indicator of host health [8, 9]. Gut microbial communities carry a wide range of important functions ranging from food processing, access to nutrients, production of antibiotics, to the influence on behavior and the immune system [10], and thus can play a key role in host adaptation to rapidly changing environments [11]. Shifts in bacterial diversity beyond the ‘normal’ range can result in dysbiosis, i.e., the depletion of commensal microbes and increase in pathogenic ones [12]. Both can have negative consequences on host immunity and health. Recent studies in humans and wildlife highlighted that gut homeostasis as subject to disturbances by viral infections, facilitating co-infections and the spread of zoonotic diseases [13,14,15]. Dysbiosis has been linked to environmental change, with further influencing factors being stress, shifts in diet, and species assemblage on the macroecological scale [13, 14, 16,17,18].

Indeed, recent studies have emphasized the negative effects of landscape modifications and habitat fragmentation on host gut microbiome homeostasis. For example, the gut microbiome of howler monkeys living in fragmented areas or in captivity had a lowered bacterial diversity [19], with several follow-up studies confirming these findings in different settings and for different species [e.g. 20, 21]. However, these studies were not able to differentiate between true landscape effects and the effects of change in diet and other anthropogenic disturbances. Specifically, whether habitat fragmentation per se or the combination of habitat fragmentation with additional anthropogenic disturbances such as contact with humans, domesticated animals, and their pathogens dominate these changes and disrupts natural gut microbiota homeostasis was so far studied only in the Tome’s spiny rat (Proechimys semispinosus) [17]. In this study, spiny rats inhabiting isolated but otherwise undisturbed rainforest fragments showed a similar gut microbial composition as conspecifics living in undisturbed continuous forests, indicating that habitat fragmentation per se is not the driving factor behind changes in the gut microbiome. However, in general, gut microbiomes were less diverse and at the same time, individual microbiomes were more heterogeneous in spiny rats inhabiting fragmented and disturbed landscapes, with differences driven by microbes associated with humans and domesticated animals. However, whether these results are restricted to this species or a common feature of generalist species living in sympatry and independent of phylogeny is largely unknown, but of high importance [22].

In the present study, we analyzed the gut microbiomes of three sympatric, generalist, neotropical, mammalian species (two marsupials Didelphis marsupialis and Philander opossum, which might be renamed Philander melanurus in the future [23], and one rodent P. semispinosus) inhabiting four landscapes in central Panama with differing degrees of anthropogenic disturbances: protected continuous tropical forests, protected forested islands in the Panama Canal, nearby unprotected forested fragments embedded in an agricultural matrix and teak plantations. The three study species represent abundant generalists living in rainforests with a high tolerance for habitat modifications [24, 25]. P. semispinosus is a strictly terrestrial rodent that mainly feeds on fruits, seeds and mycorrhiza and occupies rather small home ranges, often representing one of the most abundant mammalian species within its geographic range [26, 27]. D. marsupialis is a common terrestrial, omnivorous marsupial with a large home range [28, 29] and P. opossum is a semi-arboreal marsupial with an omnivorous diet including fruits, vertebrates, and invertebrates with a home range of 0.34 ha in Panama [29,30,31]. All three species represent important pathogen reservoirs and hosts for zoonotic diseases of clinical relevance, caused by Trypanosoma [32], Hepacivirus [24], mammalian delta virus [33], Schaalia [34] and Venezuelan equine encephalitis virus [35]. P. semispinosus and D. marsupialis are known to be consumed by humans as a source of protein, making them particularly interesting in terms of their potential roles in spillover events [36, 37]. This highlights the importance of understanding any potential health impacts anthropogenic disturbances may have on these species. With this study design, we aimed to investigate (1) to what degree, if any, do host species identity (i.e., phylogeny) and landscape modifications shape the gut microbiomes of three sympatric, generalist species; (2) to which extent do the host species differ in their resilience, i.e. responses to different degrees of anthropogenic disturbances; and finally (3) which gut bacteria are involved in driving any changes, and are they similar for all three host species? Gaining a better understanding of the gut microbial response to anthropogenic disturbances is of high importance in maintaining wildlife health, especially in generalist species that can act as vectors transmitting diseases among wildlife, but also livestock, domestic animals and ultimately humans.

Results

Species community and abundance pattern

During three field seasons, 1523 animals of 16 different species were captured in 28 study sites distributed across four landscapes (protected continuous tropical forests (=C), protected forested islands in the Panama Canal (=I), nearby unprotected forested fragments embedded in an agricultural matrix (=A) and teak plantations (=P)) differing in degree of anthropogenic disturbance in Panama (Additional file 1: Figs. S1 and S2). The three most common animals captured were P. semispinosus (n = 1235), D. marsupialis (n = 137) and P. opossum (n = 72). P. semispinosus and D. marsupialis were captured in all four landscapes, while P. opossum was absent on landscape I. Fecal samples from 397 P. semispinosus, 104 D. marsupialis, and 68 P. opossum were available for microbiome investigation (Additional file 2: Table S1).

Gut bacterial composition of the three generalist species

After sequencing the 16S rRNA gene from fecal matter and rarefication a total of 6404 different bacterial ASVs were identified in the 569 samples kept after quality filtering (Additional files 1 and 2: Fig. S3 and Table S1). The two marsupial species, P. opossum and D. marsupials, had 1454 ASVs in common, whilst the ground-dwelling species D. marsupialis and P. semispinosus shared 445 ASVs, and the microbiome of P. opossum and P. semispinosus overlapped in only 252 ASVs (Additional file 1: Fig. S4). The two opossum species shared more dominant bacterial families with each other than with the rodent species, P. semispinosus (Fig. 1). In the marsupials, the gut microbiomes were dominated by the bacterial families Lachnospiraceae, Bacteroidaceae and Acidaminococcacae. Considering less frequent bacterial families in D. marsupialis, Diplorickettsiaceae were mainly found in landscape C, while Ruminococcaceae were mainly present in the landscapes I, A and P. Oscillospiraceae were detected in landscape I, while Clostridiaceae were mainly present in the disturbed landscapes A and P. Fusobacteriaceae and Peptostreptococcaceae were mainly detected in landscape P. In P. opossum, the gut microbiomes of individuals trapped in landscape C were characterized by the presence of Clostridia_UCG-014, Enterobacteriaceae and Morganellacea, while Peptostreptococcaceae and Clostridiaceae were mainly found in the gut microbiome of individuals inhabiting landscape A and the presence of Prevotellaceae was characteristic for landscape P. Muribaculaceae (formerly known as S24-7 [38]), Lachnospiraceae and Erysipelotrichaceae were the dominant bacterial families in the gut of P. semispinosus. The gut microbial composition of P. semispinosus differed between the landscapes C and I by the addition of Oscillospiraceae in landscape I. Rodents trapped in the landscapes A and P both lacked Clostridia_UCG-014 but harbored Gastranaerophilales instead (Fig. 1).

Fig. 1
figure 1

Gut bacterial composition (at the family level) of the two marsupials (D. marsupialis, P. opossum) and the spiny rat (P. semispinosus). Details on the landscapes C, I, A and P are provided in the methods and their locations are shown in Additional file 1: Fig. S1. Note that P. opossum was not captured in landscape I

Gut bacterial composition is shaped by both host phylogeny and landscape type

We first investigated if species identity and/or landscape type shaped gut microbial beta diversity (both composition and dispersion) across all three sympatric species. Using one model with host species identity and landscape type as fixed factors per distance matrix (weighted and unweighted UniFrac distances), both species identity (PERMANOVA: weighted Unifrac: df = 2, R2 = 0.23, p ≤ 0.001; unweighted UniFrac: df = 2, R2 = 0.21, p ≤ 0.001) and landscape type had significant effects on gut microbial composition (PERMANOVA: weighted Unifrac: df = 3, R2 = 0.02, p ≤ 0.001; unweighted UniFrac: df = 3, R2 = 0.02, p ≤ 0.001; Fig. 2a, b). Based on effect sizes (Cohen’s d), the effect of phylogeny was larger than the effect of landscape type (range based on weighted UniFrac distances: host phylogeny: 1.04–7.06, landscape: 0–0.76; range based on unweighted UniFrac distances: host phylogeny: 2.82–6.19, landscape: 0.33–1.34; Fig. 2c, d). Moreover, pairwise comparisons of each host species pair showed the smallest differences in gut microbial composition between both marsupials, and larger differences between either marsupial and the rodent (Fig. 2c, d). The effect size for pairwise comparisons of landscape type showed that each landscape was distinct from one another, except for the comparison between the disturbed landscapes A and P using weighted UniFrac distances (Fig. 2c). In addition to the bacterial composition, gut microbial dispersion was also significantly affected by both host species identity (weighted UniFrac: df = 2, F value = 174.93, p ≤ 0.001; unweighted UniFrac: df = 2, F value = 122.47, p ≤ 0.001, Fig. 2a, b) and landscape type (weighted UniFrac: df = 3, F value = 29.62, p ≤ 0.001; unweighted UniFrac: df = 2, F value = 33.38, p ≤ 0.001).

Fig. 2
figure 2

Differences in beta diversity between three sympatric generalist species, the marsupials D. marsupialis (Dm, green), P. opossum (Po, orange), and the rodent P. semispinosus (Ps, blue). The top row shows NMDS plots of a weighted and b unweighted UniFrac distances. The bottom row c, d displays corresponding forest plots with Cohen’s d effect sizes highlighting the phylogenetic signal and landscape effects, i.e. the different levels of anthropogenic disturbance. Details on the landscapes C, I, A and P are provided in the methods and their locations are shown in Additional file 1: Fig. S1

Species-specific effects of landscape type on gut microbial alpha and beta diversity of sympatric generalist species

To determine if landscape type had similar effects on gut microbial alpha and beta diversity across three different, yet sympatric species, we analyzed the effect of landscape type on gut microbial alpha diversity, composition, and dispersion for each species separately. In the ground-dwelling marsupial D. marsupialis, all alpha diversity metrics were significantly impacted by landscape type (Additional file 2: Table S2), with individuals inhabiting landscape P having a significantly lower gut microbial alpha diversity than those living in the other three landscape types (Additional file 2: Table S2, depicted for Faith’s PD, Fig. 3). Beta diversity based on weighted UniFrac distances did not show a significant effect of landscape type on gut microbial composition and dispersion (PERMANOVA: df = 3, R2 = 0.04, p = 0.1; PERMDISP: df = 3, F value = 0.43, p = 0.74, Fig. 4). However, using unweighted UniFrac distances, gut microbial composition (but not dispersion; PERMDISP: df = 3, F value = 1.22, p = 0.30) was significantly affected by landscape type (PERMANOVA: df = 3, R2 = 0.04, p ≤ 0.01, Fig. 4) with plantations (landscape P) clearly differing from the other three landscapes C, I and A based on pairwise comparisons (Additional file 2: Table S2).

Fig. 3
figure 3

Landscape effects on alpha diversity (measured as Faith´s PD) of three sympatric generalist species, the marsupials D. marsupialis and P. opossum, and the spiny rat P. semispinosus. Individuals were trapped in protected continuous tropical forests (landscape C, green); protected forested islands in the Panama Canal (landscape I, blue); in the nearby unprotected forested fragments embedded in an agricultural matrix (landscape A, yellow) and in teak plantations (landscape P, red)

Fig. 4
figure 4

NMDS plots visualizing landscape effects on beta diversity. Beta diversity is measured by a weighted and b unweighted UniFrac distances in the three sympatric generalist species, the marsupials D. marsupialis (top) and P. opossum (middle), and the spiny rat P. semispinosus (bottom). Distances to the group centroids are depicted in the inserted graphs in the top right corners and ellipses indicate 95% confident intervals. Individuals were trapped in protected continuous tropical forests (landscape C, green); protected forested islands in the Panama Canal (landscape I, blue); in the nearby unprotected forested fragments embedded in an agricultural matrix (landscape A, yellow) and in teak plantations (landscape P, red)

In the second, partially arboreal, marsupial species P. opossum that was never captured on landscape I, landscape type did not have a significant effect on Faith’s PD (Additional file 2: Table S3, Fig. 3) or on gut bacterial community dispersion (PERMDISP: weighted UniFrac: df = 2, F value = 0.16, p = 0.21; unweighted UniFrac: df = 2, F value = 2.72, p = 0.07). It did, however, significantly affect gut microbial composition with regards to unweighted UniFrac distances (PERMANOVA: df = 2, R2 = 0.04, p ≤ 0.01, Fig. 4), with all landscapes differing from one another based on post-hoc pairwise comparisons, but not with regards to weighted UniFrac distances (PERMANOVA: df = 2, R2 = 0.04, p = 0.21; Fig. 4, Additional file 2: Table S3).

In the ground-dwelling, rodent species P. semispinosus, landscape type had a significant effect on Faith’s PD (Additional file 2: Table S4, Fig. 3) and pairwise comparisons showed that each landscape type significantly differed from one another, except for the two protected landscapes C and I as well as the two unprotected landscapes A and P (Additional file 2: Table S4, Fig. 3). Similarly, gut microbial beta diversity (composition and dispersion) was significantly impacted by landscape type, regardless of distance metric (weighted UniFrac: PERMANOVA: df = 3, R2 = 0.07, p ≤ 0.001; PERMDISP: df = 3, F value = 13.21, p ≤ 0.001, pairwise post-hoc test Additional file 2: Table S4; unweighted UniFrac: PERMANOVA: df = 3, R2 = 0.05, p ≤ 0.001; PERMDISP: df = 3, F value = 7.17, p ≤ 0.001, Additional file 2: Table S4), with a clear separation between protected (landscapes C and I, without contact to humans, livestock and domestic animals) and unprotected (landscapes A and P, with contact to humans, livestock and domestic animals) landscape types (Fig. 4). Thus, the gut microbiomes of generalists are affected by anthropogenic disturbances, but sympatric species show species-specific differences in their responses.

Landscape-driven species-specific taxonomic differences in gut microbial beta diversity

To gain a better understanding of which bacterial taxa are driving the differences in beta diversity between landscape types within each host species, we ran three different differential abundance analyses, one for each species, using ANCOM [39]. Because this is a pairwise test, we grouped together individuals from landscape types that did not significantly differ in the aforementioned beta diversity analyses, i.e., showed similar responses to landscape effects on beta diversity.

First, for D. marsupialis, we compared gut microbiomes from individuals inhabiting landscapes C + I + A against P and identified 27 differentially abundant ASVs from the bacterial classes Alphaproteobacteria, Bacilli, Bacteroida, Clostridia, and Gammaproteobacteria (Fig. 5). Second, for P. opossum, we compared the gut microbial communities of individuals living in the control landscape C against those from the disturbed landscapes A and P (this species was not captured on landscape I). Because no clear separation of the microbial alpha- or beta-diversity based on landscape type was observed, we decided to use our hypothesis-driven assumption, comparing the two disturbed landscapes against the protected control (landscape C). In doing so, we detected 24 differentially abundant ASVs from the bacterial classes Alphaproteobacteria, Bacteroida, Clostridia, Coriobacteriia and Vampirivibrionia (Fig. 6). Lastly, for P. semispinosus, we compared the gut microbial compositions of rodents inhabiting the protected landscapes C and I to those living in the unprotected landscapes A and P and found 481 differentially abundant ASVs from the bacterial classes Actinobacteria, Alphaproteobacteria, Bacilli, Bacteroida, Clostridia, Coriobacteriia, Elusimicrobia, Gammaproteobacteria, Negativicutes, Spirochaetia, Vampirivibrionia, and Verrucomicrobiae (Figs. 7 and 8).

Fig. 5
figure 5

Differentially abundant ASVs (using ANCOM) in pairwise comparisons with landscape groupings based on similar responses in beta diversity to landscape modifications (Fig. 4) in D. marsupialis comparing landscape C + I + A against P. Colors represent the different bacterial classes. a Differential abundant ASVs according to the taxonomical assignment; b Volcano plot of differential abundant ASVs, depicting the F statistics and W-value

Fig. 6
figure 6

Differentially abundant ASVs (using ANCOM) in pairwise comparisons with landscape groupings based on similar responses in beta diversity to landscape modifications (Fig. 4) in P. opossum comparing landscape C against A + P. Colors represent the different bacterial classes. a Differential abundant ASVs according to the taxonomical assignment; b Volcano plot of differential abundant ASVs, depicting the F statistics and W-value

Fig. 7
figure 7

First half of differentially abundant ASVs (using ANCOM) in pairwise comparisons with landscape groupings based on similar responses in beta diversity to landscape modifications (Fig. 4) in P. semispinosus comparing landscape C + I against A + P. Results split into two graphs. Colors represent the different bacterial classes. a Differential abundant ASVs according to the taxonomical assignment; b Volcano plot of differential abundant ASVs, depicting the F statistics and W-value

Fig. 8
figure 8

Second half of differentially abundant ASVs (using ANCOM) in pairwise comparisons with landscape groupings based on similar responses in beta diversity to landscape modifications (Fig. 4) in P. semispinosus comparing landscape C + I against A + P. Results split into two graphs. Colors represent the different bacterial classes. a Differential abundant ASVs according to the taxonomical assignment; b Volcano plot of differential abundant ASVs, depicting the F statistics and W-value

Among the ASVs identified were representatives, such as Chrisensenellaceae R-7 group (D. marsupialis), Chlostridium senu stricto I (D. marsupialis), Rikenella (D. marsupialis), Vagococcus (D. marsupialis), Butyricoccus (P. opossum), Coprococcus (P. opossum), Peptococcus (P. opossum), Rikenella (P. opossum), Bacteroides (P. semispinosus), Allobaculum (P. semispinosus) and Dubsosiella (P. semispinosus).

Discussion

Impact of habitat change on species with high ecological plasticity

In the Anthropocene, the selection pressures faced by species worldwide have changed to favor species possessing ecological plasticity, with a higher tolerance to anthropogenically modified environments [40], the so-called generalist species [41]. Failure to adapt can especially impact animals with narrow ecological niches, leading to changes in the composition of animal communities [42]. Increases in the abundance of generalists in disturbed landscapes at the cost of the decline of less resilient species can be accompanied by negative consequences for the ecosystem and increase the risk of disease spread and spillover [43]. As anthropogenic selection pressures will affect these pathogen reservoirs and disease vectors best adapted to living in close proximity to human settlements, studies on generalists’ role in emerging infectious diseases from wildlife and their adaptability to global environmental changes are key in order to understand the threat they are under and the threat they may pose [4, 44]. In fact, deforestation is considered one of the main drivers of zoonotic diseases [45], yet, current reforestation programs often rely on a few tree species or even monocultures, such as eucalyptus or teak plantations, which might accelerate the risk [46].

Phylogeny and landscape influence the microbiome of sympatric species

In central Panama, we investigated four landscapes differing in their degree of anthropogenic disturbance which have been shown to have severe impacts on different levels of diversity: species assemblages and abundance patterns [24, 42, 47], pathogen diversity and host infections [18, 25, 48], neutral and adaptive genetic diversity [18, 25] along with gut microbial community patterns [17, 18] differed according to landscape type. As a crucial factor in host health, the gut microbiome has been put forward as a potentially vital component in helping wildlife to adapt to fast-paced global changes [11]. However, whether this holds true for a wider range of species facing the same degrees of anthropogenic disturbance was the focus of the present study.

Here, we showed that three sympatric, generalist species faced with the same landscape-level anthropogenic disturbances each responded differently with regards to gut microbial diversity and that the gut microbiome of sympatric species is determined by both phylogeny and landscape. We observed species-specific gut microbial responses to different landscape types, with P. semispinosus showing the highest sensitivity to proximity to humans. In D. marsupials, however, only the landscape with a combination of human contact, fragmentation, and complete changes in forest type (from lowland tropical rainforest to teak plantation) led to changes in the gut microbiome. In contrast, P. opossum’s gut microbiome showed the highest resilience, without a distinct pattern of change in the gut microbial alpha diversity or beta dispersion. These species-specific responses were mirrored in our differential abundance analyses that showed the highest number of differentially abundant ASVs in P. semispinosus versus the two other species. The gut microbiomes of our three analyzed species were affected by both species identity and landscape type, which is supported by other studies [49, 50].

Species identity can influence the gut microbiome not only because of species-specific genetic constitutions, but also because species identity determines suitable environment, food preferences and processing potential, pathogen resistance, as well as adaptive potential [51]. Similarly, landscape types can differ in food availability as well as exposure to pathogens, pesticides, and different plant and animal communities. The gut microbiomes of the two opossum species were more similar to one another than they were to the spiny rat. A similar trend was shown in large African mammals, where gut microbial composition was closely correlated with phylogeny [52], as did the gut microbiome in some primates [53]. While some studies show the effect of species identity and environment on several species’ gut microbiomes, [22, 49], they are often limited to extreme environmental situations [22, 54], thus only comparing extremely contrasting conditions for the animals. However, landscape plays a large role in shaping the microbial community in humans and apes [55], and habitat fragmentation can have an influence on the gut microbiome through dietary changes [56]. Similar results have been discovered in different howler monkey species, where host species was a key predictor of the gut microbiota, but forest type, habitat, and season explained species-specific variances [57].

Only extreme habitat perturbation changes D. marsupialis gut microbiomes

Analyzing the three species separately showed that changes in gut microbial composition in D. marsupialis occurred only in individuals inhabiting the most disturbed landscape, namely teak plantations. Teak plantations are defined by a different forest type in addition to being fragmented and embedded in a matrix with human contact. Because they are monocultures, teak plantations can alter the food availability of the inhabiting animals and provide a harsher environment, e.g. due to increased risk of forest fires [58]. Additionally, the conversion of natural habitats to teak plantations exposes animal communities to changes in microbial soil composition [59]. D. marsupialis individuals living in these plantations either need to adapt to these changes and the potential stress associated with them or use the plantations as a corridor for travel.

High resilience of P. opossum gut microbiomes to anthropogenic changes

The second marsupial, P. opossum, however, showed a different gut microbial response. Despite only occurring in three of the four sampled landscapes, P. opossum individuals showed little change in gut microbial composition with no clear patterns across the different landscape types. Changes in gut microbial composition between individuals inhabiting all three landscapes were only observed when using unweighted UniFrac distances, indicating that rare ASVs are likely responsible for these changes. As a semi-arboreal habitat generalist, P. opossum individuals could escape bacterial taxa introduced by humans and livestock to soils. This could explain the weakened gut microbial response compared to the more closely-related D. marsupialis. Further risk of coming into contact with microbes from humans or domesticated animals is reduced because part of P. opossum's diet consists of insects [60], meaning this host might be less exposed to bacteria in the ground or associated with plants compared to, say, herbivorous hosts. The microbiomes of animals and their environment are linked and shape each other in both directions [61, 62]. Still, the fact that prolonged isolation and fragmentation possibly lead to local extinction (as seen by their absence on islands) could indicate that P. opossum is not insensitive to landscape changes. Because islands provide the harshest matrix based on accessibility and usability as a corridor for movement, their extinction on islands could be due to a past bottleneck event. However, there are interesting dynamics of opossum species in the Panama Canal area with local disappearance as well as re-colonization of some small islands close to the mainland, thus opening the possibility for re-colonization [63].

Large impact of human vicinity on P. semispinosus gut microbiomes

The gut microbiome of the common P. semispinosus is very sensitive to human contact, as was shown by the reduced microbial diversity and shifts in composition, driven by bacterial taxa associated with humans and domesticated animals, as well as increased dispersion of the bacterial community in individuals from fragments surrounded by agriculture [17]. All these features were not detected in individuals sampled in continuous forests or protected forested islands [17]. While the impact of protected landscape types with no human contact (continuous forests and islands) versus fragmented landscapes with human contact (forest fragments embedded in an agricultural matrix) on P. semispinosus gut microbiomes have been discussed elsewhere [17], the apparent similar behavior of the gut microbiome from individuals living in teak plantation to those inhabiting forest fragments is noteworthy. Contrarily to forested fragments, teak plantations also differ in the natural assembly of trees, with the dominant species being the from Asia introduced teak tree (Tectona grandis). The gut microbiome of D. marsupialis was only impacted in this landscape type, which we consider to be the most extreme landscape type, and therefore, we expected an additional impact on the P. semispinosus gut microbiome on top of the changes caused by the forest fragments. However, the lack of a compounded impact could mean two things: (1) individuals from fragmented forests already display the most extreme gut microbial perturbation tolerable by this species, or (2) teak plantations provide a suitable forest habitat for P. semispinosus, despite being monocultures and lacking trees typical for rainforests. It also indicates that the presence of humans and domesticated animals in the matrix alone shapes the microbiome, and not fragmentation or forest tree composition per se.

Human-vicinity-associated ASVs drive gut microbial changes

Differential abundance analysis revealed a large difference in the number of ASVs determined to be differently abundant between the landscapes when considering each of the three species separately. While for the two opossum species only a handful of differentially abundant bacterial taxa were detected, many were discovered for P. semispinosus. In general, detected bacteria often showed similarity to ones characterized in more detail in human or domesticated animal microbiomes (e.g. Allobaculum [64]). Most of the detected ASVs have already been discussed [17], theorizing that a large portion of the identified bacteria could have been introduced by humans and their domesticated animals. In addition, in D. marsupialis, we detected ASVs assigned to the genus Vagococcus, which was first isolated from chicken feces [65], animals that D. marsupials are known to get into close contact with because they prey on their eggs [36]. Landscapes A and P are in close vicinity to human settlements, therefore uptake of these human-driven taxa into their gut microbiomes is not unexpected and could cause harm to the host by changing gut microbial composition.

Potential consequences of microbiome changes induced by anthropogenic landscape disturbance

D. marsupialis is known to be consumed by humans, both for nutritional and cultural reasons [36]. Thus, not only proximity, but also direct contact with humans can be an important pathway for pathogen transmission, and a perturbed microbiome might have negative effects on host immunity, given the microbiome’s interplay and crosstalk with the host immune system [66]. Altered microbial communities can facilitate pathogen infection [48] and vice versa pathogen infection can change the microbiome [13,14,15]. These perturbed microbiomes further increase the risk of horizontal gene transfer which could lead to pathogenic bacteria [67]. Close contact between humans and domesticated animals can be dangerous for wildlife and humans, as demonstrated with the spillover of Nipah virus, which made the jump from flying foxes to humans through pigs as an intermediate host [68, 69]. Because land-use change can cause pandemics and the emergence of new diseases [70] and because rodents and marsupials represent a significant zoonotic disease risk in the future [71], these findings become even more important regarding the potential for the origins and emergence of zoonotic diseases, and future studies will reveal how these microbiome changes impact animal’s fitness in detail. Moreover, the loss of microbial diversity has been recognized as a potential threat to the discovery of new drugs or therapeutics in the field of microbial biotechnology [72].

Conclusion

Overall, we could show that the gut microbiomes of sympatric species inhabiting landscapes with differing degrees of anthropogenic disturbance are mainly shaped by host species identity, with landscape type playing a smaller but significant role. Interestingly, there is a species-specific gut microbial response to landscape type, indicating that findings from one species cannot always be generalized to other species, even not to those living in the same habitat or to those that are closely related. This shows that these generalists' gut microbiomes are sensitive to landscape-level changes, a fact not detected by biodiversity monitoring of vertebrates, and that these changes are not uniform across host species. Thus, even for generalists, environmental changes can pose a big impact, which is an important finding as this might directly cause consequences and the risk of emerging zoonoses not only for wildlife health, but also for the health of domesticated animals intended for human consumption, and for humans themselves.

Material and methods

Study area and sampling

This study was carried out in the Panama Canal area, Panama. Animals were captured in 28 study sites grouped into four different landscapes based on their degree of anthropogenic disturbance (Additional file 1: Fig. S1, map created with ggmap [73]): (1) continuous rainforest (=C, five capture sites), i.e., undisturbed and protected lowland tropical rainforest within the Barro Colorado Nature Monument; (2) forested islands (=I, six capture sites) situated in the Gatún Lake and also protected by the Barro Colorado Nature Monument, i.e., fragmented but otherwise undisturbed landscape; (3) fragmented and disturbed (i.e. contact to humans and domesticated animals) tropical lowland rainforest embedded in an agricultural matrix (=A, nine capture sites); and (4) fragmented and disturbed teak plantations (=P, seven capture sites) planted by humans and mainly consisting of Tectona grandis.

Field work took place during three field seasons (October 2013 to May 2014, October 2014 to May 2015 and September 2016 to April 2017, alternating the order of capture sites between seasons). At each capture site a trapping grid consisting of 100 stations was set up. Each station was separated by 20 m and consisted of three live traps [one Tomahawk trap (size: 15.2 × 15.2 × 48.3 cm, www.livetrap.com) and two Sherman traps (size: 10.2 × 11.4 × 38.1 cm, www.shermantraps.com)], one of which was placed above ground, if possible, i.e., on trees, lianas or similar, if available, to include arboreal species. Traps were opened at dusk and baited with a mixture of peanut butter, oatmeal, bird seeds, banana, and dog food to attract species with various dietary preferences and controlled at dawn of the next day. Then, captured animals were measured, individually marked to recognize recaptures, and fecal samples were taken from animals during sampling. Afterwards, the animals were released at the capturing location (further details see [24]). Fecal samples were stored in collection tubes containing RNAlater at − 20 °C until DNA extraction.

DNA extraction, 16S rRNA gene amplification and sequencing

A detailed summary of sample processing is described elsewhere [17]. In brief, we extracted DNA from fecal samples from a total of 793 samples (including extraction blanks and PCR controls) using NucleoSpin®Soil-Extraction Kit from Macherey–Nagel (Germany). The final elution step was performed twice with 50 µl of elution buffer each time, resulting in a total volume of 100 µl. Following extraction, we amplified the 291 nucleotide-long V4 region of the 16S rRNA gene using the 515 F and 806 R primers [74, 75] applying a two-step polymerase chain reaction (PCR). The first step was an initial denaturation of 600 s at 95 °C followed by 30 cycles with 95 °C for 30 s, 60° for 30 s and 72 °C for 45 s, followed by a final elongation of 72 °C for 600 s. The second step consisted of ten cycles for the barcoding with the same conditions as described above. The samples were sequenced on six runs on an Illumina MiSeq at our Institute of Evolutionary Ecology and Conservation Genomics, Ulm University, Germany.

Data processing

Reads from all six Illumina runs were analyzed in QIIME 2 (Version 2020.6, August 2020) [76] with a total of 14,489,625 sequences. DADA2 [77] was used to process the sequences and assemble these into amplicon sequence variants (ASVs) and SILVA (version 1.38) [78] was used as a taxonomic reference database. We removed a total of 194,432 sequences (roughly 0.83%) annotated as Archaea, mitochondria and chloroplast. Data were transferred into a phyloseq object [79] within the R environment (version 4.0.2) [80] for further analyses. ASVs identified in the blanks and controls (a total of 55,463 ASVs, equivalent to 0.38%) were removed from samples to avoid false results. We applied an additional filter to remove rare ASVs with fewer than 50 reads across the entire dataset and which occurred in only 2% of all the samples. Finally, we rarefied the data in order to control for uneven sequencing depth. Rarefaction was performed using the rarefy_even_depth function from the phyloseq package and a sequencing depth of 10,000 reads was chosen based on rarefaction curves (Additional file 1: Fig. S3). After performing all the bioinformatic quality filtering steps at rarefaction, 569 samples (P. semispinosus: n = 397 individuals; D. marsupialis: n = 104 individuals; P. opossum: n = 68) remained for subsequent analysis (Additional file 2: Table S1).

Statistical analyses

Statistical analyses were separated into two parts based on the study questions. First, to determine if phylogeny and/or landscape affect the gut microbial beta diversity in the three species, we used the whole dataset, consisting of all three species and all four landscapes. Using weighted and unweighted UniFrac distances [81], we tested for differences in gut microbial composition and homogeneity according to landscape and host species using PERMANOVA (Permutational Analysis of Variance, with 9999 permutations) and PERMDISP2 (Permutational Analysis of Multivariate Dispersions) [82] using the vegan package [83]. These tests were followed by post-hoc pairwise comparisons. Effect sizes (Cohen’s d) [84] for pairwise comparisons were calculated by extracting the results from the PERMANOVA of the first two NMDS axes.

Second, to determine if the gut microbiome in each of the three host species is similarly impacted by landscape type, we subset our data for each species. Both alpha diversity (observed number of ASVs, Shannon Diversity and Faith’s phylogenetic diversity, PD) and beta diversity metrics (weighted and unweighted UniFrac distances) were calculated. For alpha diversity, we constructed generalized linear models (GLMs) with landscape, season, and sex as factors explaining the alpha diversity indices and afterwards used contrasts to analyze pairwise comparisons. For beta diversity, we applied PERMANOVAs and PERMDISPs as described above, followed by post-hoc pairwise comparisons calculated for each landscape.

Finally, to investigate which taxa were driving the differences in beta diversity, we performed differential abundance analyses using ANCOM (Analysis of Composition of Microbiomes) [39] which allows two-level factor comparisons [17]. We chose a w0 of 0.7 as originally described in [39]. Based on the results of the species-specific landscape effects on beta diversity (see Results), we compared the landscapes C + I + A versus P for D. marsupials, C versus A + P for P. opossum and C + I versus A + P for P. semispinosus for our two-level factor comparisons. All graphs were plotted in the R environment using the ggplot2 package [85].

Availability of data and materials

The datasets generated and/or analyzed during the current study will be made available in upon acceptance.

References

  1. Jones KE, Patel NG, Levy MA, Storeygard A, Balk D, Gittleman JL, et al. Global trends in emerging infectious diseases. Nature. 2008;451:990–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Ostfeld RS, Keesing F. Biodiversity series: the function of biodiversity in the ecology of vector-borne zoonotic diseases. Can J Zool. 2000;78:2061–78.

    Article  Google Scholar 

  3. Fahrig L. Effects of habitat fragmentation on biodiversity. Annu Rev Ecol Evol Syst. 2003;34:487–515.

    Article  Google Scholar 

  4. Gibb R, Redding DW, Chin KQ, Donnelly CA, Blackburn TM, Newbold T, et al. Zoonotic host diversity increases in human-dominated ecosystems. Nature. 2020;584:398–402.

    Article  CAS  PubMed  Google Scholar 

  5. Ostfeld RS, Keesing F. Effects of host diversity on infectious disease. Annu Rev Ecol Evol Syst. 2012;43:157–82.

    Article  Google Scholar 

  6. Goldberg TL, Gillespie TR, Rwego IB, Estoff EL, Chapman CA. Forest fragmentation as cause of bacterial transmission among nonhuman primates, humans, and livestock. Uganda Emerg Infect Dis. 2008;14:1375–82.

    Article  PubMed  Google Scholar 

  7. Alpízar P, Risely A, Tschapka M, Sommer S. Agricultural fast food: bats feeding in banana monocultures are heavier but have less diverse gut microbiota. Front Ecol Evol. 2021;9:746783.

    Article  Google Scholar 

  8. Marietta E, Rishi A, Taneja V. Immunogenetic control of the intestinal microbiota. Immunology. 2015;145:313–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Oliveira RA, Ng KM, Correia MB, Cabral V, Shi H, Sonnenburg JL, et al. Klebsiella michiganensis transmission enhances resistance to Enterobacteriaceae gut invasion by nutrition competition. Nat Microbiol. 2020;5:630–41.

    Article  CAS  PubMed  Google Scholar 

  10. Barko PC, McMichael MA, Swanson KS, Williams DA. The gastrointestinal microbiome: a review. J Vet Intern Med. 2018;32:9–25.

    Article  CAS  PubMed  Google Scholar 

  11. Alberdi A, Aizpurua O, Bohmann K, Zepeda-Mendoza ML, Gilbert MTP. Do vertebrate gut metagenomes confer rapid ecological adaptation? Trends Ecol Evol. 2016;31:689–99.

    Article  PubMed  Google Scholar 

  12. Schippa S, Conte MP. Dysbiotic events in gut microbiota: impact on human health. Nutrients. 2014;6:5786–805.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Wasimuddin, Brändel SD, Tschapka M, Page R, Rasche A, Corman VM, et al. Astrovirus infections induce age-dependent dysbiosis in gut microbiomes of bats. ISME J. 2018;12:2883–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Wasimuddin, Corman VM, Ganzhorn JU, Rakotondranary J, Ratovonamana YR, Drosten C, et al. Adenovirus infection is associated with altered gut microbial communities in a non-human primate. Sci Rep. 2019;9:13410.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Schmid DW, Fackelmann G, Wasimuddin, Rakotondranary J, Ratovonamana YR, Montero BK, et al. A framework for testing the impact of co-infections on host gut microbiomes. Anim Microbiome. 2022;4:48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Ingala MR, Becker DJ, Holm JB, Kristiansen K, Simmons NB. Habitat fragmentation is associated with dietary shifts and microbiota variability in common vampire bats. Ecol Evol. 2019;9(11):6508–23.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Fackelmann G, Gillingham MAF, Schmid J, Heni AC, Wilhelm K, Schwensow N, et al. Human encroachment into wildlife gut microbiomes. Commun Biol. 2021;4:1–11.

    Article  Google Scholar 

  18. Schwensow NI, Heni AC, Schmid J, Montero BK, Brändel SD, Halczok TK, et al. Disentangling direct from indirect effects of habitat disturbance on multiple components of biodiversity. J Anim Ecol. 2022;91:2200–34.

    Article  Google Scholar 

  19. Amato KR, Yeoman CJ, Kent A, Righini N, Carbonero F, Estrada A, et al. Habitat degradation impacts black howler monkey (Alouatta pigra) gastrointestinal microbiomes. ISME J. 2013;7:1344–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Frankel JS, Mallott EK, Hopper LM, Ross SR, Amato KR. The effect of captivity on the primate gut microbiome varies with host dietary niche. Am J Primatol. 2019;81: e23061.

    Article  PubMed  Google Scholar 

  21. Jiménez RR, Alvarado G, Sandoval J, Sommer S. Habitat disturbance influences the skin microbiome of a rediscovered neotropical-montane frog. BMC Microbiol. 2020;20:292.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Lavrinienko A, Hämäläinen A, Hindström R, Tukalenko E, Boratyński Z, Kivisaari K, et al. Comparable response of wild rodent gut microbiome to anthropogenic habitat contamination. Mol Ecol. 2021;30:3485–99.

    Article  CAS  PubMed  Google Scholar 

  23. Voss RS, Díaz-Nieto JF, Jansa SA. A revision of Philander (Marsupialia: Didelphidae), part 1: P. quica, P. canus, and a new species from Amazonia. Novi. 2018;2018:1–70.

    Article  Google Scholar 

  24. Schmid J, Rasche A, Eibner G, Jeworowski L, Page RA, Corman VM, et al. Ecological drivers of Hepacivirus infection in a neotropical rodent inhabiting landscapes with various degrees of human environmental change. Oecologia. 2018;188:289–302.

    Article  PubMed  Google Scholar 

  25. Heni AC, Schmid J, Rasche A, Corman VM, Drosten C, Sommer S. Pathogen-associated selection on innate immunity genes (TLR4, TLR7) in a neotropical rodent in landscapes differing in anthropogenic disturbance. Heredity. 2020;125:184–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Adler GH. Fruit and seed exploitation by Central American spiny rats, Proechimys semispinosus. Stud Neotrop Fauna Environ. 1995;30:237–44.

    Article  Google Scholar 

  27. Mangan SA, Adler GH. Consumption of arbuscular mycorrhizal fungi by spiny rats (Proechimys semispinosus) in eight isolated populations. J Trop Ecol. 1999;15:779–90.

    Article  Google Scholar 

  28. Sunquist ME, Austad SN, Sunquist F. Movement patterns and home range in the common opossum (Didelphis marsupialis). J Mammal. 1987;68:173–6.

    Article  Google Scholar 

  29. Julien-Laferriere D, Atramentowicz M. Feeding and reproduction of three didelphid marsupials in two neotropical forests (French Guiana). Biotropica. 1990;22:404–15.

    Article  Google Scholar 

  30. de la Sancha N, Pérez-Hernandez R, Costa LP, Brito D, Cáceres N. Philander opossum. The IUCN Red List of Threatened Species 2016:e.T40516A22176779. https://doi.org/10.2305/IUCN.UK.2016-2.RLTS.T40516A22176779.en. Accessed 6 Mar 2023.

  31. Hunsaker D, editor. The biology of marsupial. In: The biology of marsupials. Academic Press; 1977. pp. iii.

  32. Travi BL, Arteaga LT, León AP, Adler GH. Susceptibility of spiny rats (Proechimys semispinosus) to Leishmania (Viannia) panamensis and Leishmania (Leishmania) chagasi. Mem Inst Oswaldo Cruz. 2002;97:887–92.

    Article  CAS  PubMed  Google Scholar 

  33. Paraskevopoulou S, Pirzer F, Goldmann N, Schmid J, Corman VM, Gottula LT, et al. Mammalian deltavirus without hepadnavirus coinfection in the neotropical rodent Proechimys semispinosus. PNAS. 2020;117:17977–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Knoepfler S, Schauer A, Thomann A, Feyer S, Rüegg-van den Broek P, Glardon OJ, et al. Actinomycosis in a gray four-eyed opossum (Philander opossum) caused by a novel species of Schaalia. BMC Veterinary Res. 2021;17:243.

    Article  CAS  Google Scholar 

  35. Salas RA, Garcia CZ, Liria J, Barrera R, Navarro JC, Medina G, et al. Ecological studies of enzootic Venezuelan equine encephalitis in north-central Venezuela, 1997–1998. Am J Trop Med Hyg. 2001;64:84–92.

    Article  CAS  PubMed  Google Scholar 

  36. Barros FB, de Aguiar AP. Common opossum (Didelphis marsupialis Linnaeus, 1758): food and medicine for people in the Amazon. J Ethnobiol Ethnomed. 2014;10:65.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Asprilla-Perea J, Díaz-Puente JM, Martín-Fernández S. Estimating the potential of wild foods for nutrition and food security planning in tropical areas: experimentation with a method in Northwestern Colombia. Ambio. 2022;51:955–71.

    Article  PubMed  Google Scholar 

  38. Lagkouvardos I, Lesker TR, Hitch TCA, Gálvez EJC, Smit N, Neuhaus K, et al. Sequence and cultivation study of Muribaculaceae reveals novel species, host preference, and functional potential of this yet undescribed family. Microbiome. 2019;7:28.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Mandal S, Van Treuren W, White RA, Eggesbø M, Knight R, Peddada SD. Analysis of composition of microbiomes: a novel method for studying microbial composition. Microb Ecol Health Dis. 2015. https://doi.org/10.3402/mehd.v26.27663.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Glidden CK, Nova N, Kain MP, Lagerstrom KM, Skinner EB, Mandle L, et al. Human-mediated impacts on biodiversity and the consequences for zoonotic disease spillover. Curr Biol. 2021;31:R1342–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Otto SP. Adaptation, speciation and extinction in the Anthropocene. Proc Biol Sci. 2018;285:20182047.

    PubMed  PubMed Central  Google Scholar 

  42. Brändel SD, Hiller T, Halczok TK, Kerth G, Page RA, Tschapka M. Consequences of fragmentation for neotropical bats: the importance of the matrix. Biol Cons. 2020;252: 108792.

    Article  Google Scholar 

  43. Johnson CK, Hitchens PL, Pandit PS, Rushmore J, Evans TS, Young CCW, et al. Global shifts in mammalian population trends reveal key predictors of virus spillover risk. Proc Biol Sci. 2020;287:20192736.

    PubMed  PubMed Central  Google Scholar 

  44. Keesing F, Ostfeld RS. Impacts of biodiversity and biodiversity loss on zoonotic diseases. PNAS. 2021;118:e2023540118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Morand S, Lajaunie C. Outbreaks of vector-borne and zoonotic diseases are associated with changes in forest cover and oil palm expansion at global scale. Front Vet Sci. 2021;8: 661063.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Shah V, Shah A, Joshi V. Predicting the origins of next forest-based emerging infectious disease. Environ Monit Assess. 2018;190:337.

    Article  PubMed  Google Scholar 

  47. Marklewitz M, Dutari L, Paraskevopoulou S, Page R, Loaiza J, Junglen S. Diverse novel phleboviruses in sandflies from the Panama Canal area, Central Panama. J General Virol. 2019;100:938–49.

    Article  CAS  Google Scholar 

  48. Fleischer R, Schmid DW, Wasimuddin, Brändel SD, Rasche A, Corman VM, et al. Interaction between MHC diversity and constitution, gut microbiota and Astrovirus infections in a neotropical bat. Mol Ecol. 2022;31:3342–59.

    Article  PubMed  Google Scholar 

  49. Knowles SCL, Eccles RM, Baltrūnaitė L. Species identity dominates over environment in shaping the microbiota of small mammals. Ecol Lett. 2019;22:826–37.

    Article  CAS  PubMed  Google Scholar 

  50. Fu H, Zhang L, Fan C, Liu C, Li W, Cheng Q, et al. Environment and host species identity shape gut microbiota diversity in sympatric herbivorous mammals. Microb Biotechnol. 2021;14:1300–15.

    Article  CAS  PubMed  Google Scholar 

  51. Youngblut ND, Reischer GH, Walters W, Schuster N, Walzer C, Stalder G, et al. Host diet and evolutionary history explain different aspects of gut microbiome diversity among vertebrate clades. Nat Commun. 2019;10:2200.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Kartzinel TR, Hsing JC, Musili PM, Brown BRP, Pringle RM. Covariation of diet and gut microbiome in African megafauna. Proc Natl Acad Sci U S A. 2019;116:23588–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Amato KR, Sanders J, Song SJ, Nute M, Metcalf JL, Thompson LR, et al. Evolutionary trends in host physiology outweigh dietary niche in structuring primate gut microbiomes. ISME J. 2019;13:576–87.

    Article  CAS  PubMed  Google Scholar 

  54. Lavrinienko A, Mappes T, Tukalenko E, Mousseau TA, Møller AP, Knight R, et al. Environmental radiation alters the gut microbiome of the bank vole Myodes glareolus. ISME J. 2018;12:2801–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Amato KR, Mallott EK, McDonald D, Dominy NJ, Goldberg T, Lambert JE, et al. Convergence of human and old world monkey gut microbiomes demonstrates the importance of human ecology over phylogeny. Genome Biol. 2019;20:201.

    Article  PubMed  PubMed Central  Google Scholar 

  56. McManus N, Holmes SM, Louis EE, Johnson SE, Baden AL, Amato KR. The gut microbiome as an indicator of habitat disturbance in a critically endangered lemur. BMC Ecol Evol. 2021;21:222.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Amato KR, Martinez-Mota R, Righini N, Raguet-Schofield M, Corcione FP, Marini E, et al. Phylogenetic and ecological factors impact the gut microbiota of two Neotropical primate species. Oecologia. 2016;180:717–33.

    Article  PubMed  Google Scholar 

  58. Bell TIW. Erosion in the trinidad teak plantations. Commonw For Rev. 1973;52(153):223–33.

    Google Scholar 

  59. de Gannes V, Bekele I, Dipchansingh D, Wuddivira MN, De Cairies S, Boman M, et al. Microbial community structure and function of soil following ecosystem conversion from native forests to teak plantation forests. Front Microbiol. 2016;7:1976.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Atramentowicz M. Dynamique de population chez trois Marsupiaux Didelphides de Guyane. Biotropica. 1986;18:136–49.

    Article  Google Scholar 

  61. Dutton CL, Subalusky AL, Sanchez A, Estrela S, Lu N, Hamilton SK, et al. The meta-gut: community coalescence of animal gut and environmental microbiomes. Sci Rep. 2021;11:23117.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. McFrederick QS, Wcislo WT, Taylor DR, Ishak HD, Dowd SE, Mueller UG. Environment or kin: Whence do bees obtain acidophilic bacteria? Mol Ecol. 2012;21:1754–68.

    Article  PubMed  Google Scholar 

  63. Boyett WD, Endries MJ, Adler GH. Colonization-extinction dynamics of opossums on small islands in Panama. Can J Zool. 2000;78:1972–9.

    Article  Google Scholar 

  64. Greetham HL, Gibson GR, Giffard C, Hippe H, Merkhoffer B, Steiner U, et al. Allobaculum stercoricanis gen. nov., sp. nov., isolated from canine feces. Anaerobe. 2004;10:301–7.

    Article  CAS  PubMed  Google Scholar 

  65. Giannattasio-Ferraz S, Ene A, Maskeri L, Oliveira AP, Barbosa-Stancioli EF. Putonti C. Vagococcus fluvialis isolation and sequencing from urine of healthy cattle. G3 (Bethesda). 2021;11:034.

    Article  Google Scholar 

  66. Lambring CB, Siraj S, Patel K, Sankpal UT, Mathew S, Basha R. Impact of the microbiome on the immune system. Crit Rev Immunol. 2019;39:313–28.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Stecher B, Maier L, Hardt W-D. “Blooming” in the gut: how dysbiosis might contribute to pathogen evolution. Nat Rev Microbiol. 2013;11:277–84.

    Article  CAS  PubMed  Google Scholar 

  68. Field H, Young P, Yob JM, Mills J, Hall L, Mackenzie J. The natural history of Hendra and Nipah viruses. Microbes Infect. 2001;3:307–14.

    Article  CAS  PubMed  Google Scholar 

  69. Hooper P, Zaki S, Daniels P, Middleton D. Comparative pathology of the diseases caused by Hendra and Nipah viruses. Microbes Infect. 2001;3:315–22.

    Article  CAS  PubMed  Google Scholar 

  70. Tonissen P. IPBES #PandemicsReport: Escaping the “Era of Pandemics.” IPBES secretariat. 2020. http://ipbes.net/pandemics. Accessed 21 Nov 2021.

  71. Shivaprakash KN, Sen S, Paul S, Kiesecker JM, Bawa KS. Mammals, wildlife trade, and the next global pandemic. Curr Biol. 2021;31:3671–7.

    Article  CAS  PubMed  Google Scholar 

  72. Williams CL, Caraballo-Rodríguez AM, Allaband C, Zarrinpar A, Knight R, Gauglitz JM. Wildlife-microbiome interactions and disease: exploring opportunities for disease mitigation across ecological scales. Drug Discov Today Dis Model. 2018;28:105–15.

    Article  Google Scholar 

  73. Kahle D, Wickham H. ggmap: spatial visualization with ggplot2. The R Journal. 2013;5:144–61.

    Article  Google Scholar 

  74. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci U S A. 2011;108(Suppl 1):4516–22.

    Article  CAS  PubMed  Google Scholar 

  76. 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:852–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: High resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, et al. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007;35:7188–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. McMurdie PJ, Holmes S. phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8: e61217.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org/. 2015.

  81. Lozupone C, Knight R. UniFrac: a new phylogenetic method for Ccmparing microbial communities. Appl Environ Microbiol. 2005;71:8228–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Anderson MJ, Walsh DCI. PERMANOVA, ANOSIM, and the mantel test in the face of heterogeneous dispersions: What null hypothesis are you testing? Ecol Monogr. 2013;83:557–74.

    Article  Google Scholar 

  83. Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: community ecology package. R package version. 2019;2(5–6):2020.

    Google Scholar 

  84. Cohen J. Statistical power analysis for the behavioral sciences. 2nd ed. New York: Psychology Press; 2009.

    Google Scholar 

  85. Wickham H. ggplot2. New York, NY: Springer, New York; 2009.

    Book  Google Scholar 

Download references

Acknowledgements

We are grateful for access to the infrastructure and logistical support provided by the Smithsonian Tropical Research Institute (STRI, Panama), especially to Rachel Page, Oris Acevedo, Belkis Jimenez, and Melissa Cano. The whole study could not be performed without the immense help in the field, we would like to thank Tatjana Stooß, Robin Lietz, Manuel Stech Domene, Ruben Gunzenhäuser, Matthias Hammer, José Cascante Jiménez, José Ramírez Fernnández, José Ureña, Paula Ledezma Campos, Ricardo Sanchéz Calderón, and Stefan Dominik Brändel for dealing with permit paperwork. Ulrike Stehle and Gerd Maier were huge assets in the lab and we are thankful for the support of Luis Víquez and Randall Jiménez with bioinformatics and statistics. We are really grateful for the helpful comments of two anonymous reviewers.

Funding

Open Access funding enabled and organized by Projekt DEAL. This project is part of and funded by the DFG (German Science Foundation) Priority Program SPP 1596/2 Ecology and Species Barriers in Emerging Infectious Diseases (SO 428/9-1, 9-2 and DR 772/8-1).

Author information

Authors and Affiliations

Authors

Contributions

ACH performed study site selection and carried out fieldwork, designed present study, curated and analyzed the data, and wrote the manuscript; GF supported data analysis and edited the manuscript; GE contributed to study site selection and performed fieldwork; SW carried out lab work; JS performed study site selection, carried out fieldwork and curated data; NIS curated data; JW carried out lab work; KW carried out lab work, and supervised the molecular lab; SS initiated the overall project design, acquired funding, designed and supervised the study, supported data analysis, and edited the manuscript. All authors approved the final version.

Corresponding authors

Correspondence to Alexander Christoph Heni or Simone Sommer.

Ethics declarations

Ethics approval and consent to participate

Ethical approval was provided by the Smithsonian IACUC protocol 2013-0401-2016-A1-A7 and 2016-0627-2019-A1-A2 and the exportation of the samples was approved by the Panamanian Government (SE/A-21-14–SE/A-12-18, SEX/A-22-15, SEX/A-24-17, SEX/A-120-16, and SEX/A-52-17).

Consent for publication

All authors approved the final version.

Competing interests

The authors declare that there are no competing interests.

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: Fig. S1.

Location of the study area and 28 capture sites distributed across four landscapes differing in their anthropogenic impact in central Panama. Capture sites in the protected continuous tropical forests (landscape C) are marked in green; sites in the protected forested islands in the Panama Canal (landscape I) are marked in blue; sites in the nearby unprotected forested fragments embedded in an agricultural matrix (landscape A) are marked in yellow; and sites in teak plantations (landscape P) are marked in red. Map created with the R package ggmap (Kahle and Wickham, 2013) with the origin of the map material being Google Maps. Fig. S2. Distribution of the captured species across the four landscapes. Details on the landscapes C, I, A and P are provided in the methods and their locations are shown in Additional file 1: Fig. S1. Fig. S3. Rarefaction curves showing the number of detected ASVs in relation to 16S rRNA gene sequencing depth (i.e. total number of reads obtained per individual after quality filtering) for Didelphis marsupialis (turquoise), Philander opossum (orange) and Proechimys semispinosus (light-blue). The maximum diversity is reached at around 10,000 reads (vertical line). Fig. S4. Shared ASVs between Didelphis marsupialis (turquoise), Philander opossum (orange) and Proechimys semispinosus (light-blue).

Additional file 2: Table S1.

Number of samples per species and landscape in the final dataset. Details on the landscapes C, I, A and P are provided in the methods and their locations are shown in Additional file 1: Fig. S1. Table S2. Effects of landscape type on the gut bacterial diversity of D. marsupialis. Results from generalized linear models indicating the effects of landscape type, field season and sex on alpha diversity using a Faith’s PD; b Number of ASVs and c Shannon diversity. Results from pairwise comparisons (Contrasts) of landscapes using d Faith’s PD; e Number of ASVs and f Shannon diversity. g Results from PERMANOVA for pairwise comparisons of landscapes on beta diversity (weighted and unweighted UniFrac). SE Standard error; df Degrees of freedom. Table S3. Effects of landscape type on the gut bacterial diversity of P. opossum. Results from generalized linear models indicating the effects of landscape type, field season and sex on alpha diversity using a Faith’s PD; b Number of ASVs and c Shannon diversity. Results from pairwise comparisons (Contrasts) of landscapes using d Faith’s PD; e Number of ASVs and f Shannon diversity. g Results from PERMANOVA for pairwise comparisons of landscapes on beta diversity (weighted and unweighted UniFrac). SE Standard error; df Degrees of freedom. Table S4. Effects of landscape type on the gut bacterial diversity of P. semispinosus. Results from generalized linear models indicating the effects of landscape type, field season and sex on alpha diversity using a Faith’s PD; b Number of ASVs and c Shannon diversity. Results from pairwise comparisons (Contrasts) of landscapes using d Faith’s PD; e Number of ASVs and f Shannon diversity. g Results from PERMANOVA for pairwise comparisons of landscapes on beta diversity (weighted and unweighted UniFrac). SE Standard error; df Degrees of freedom.

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 http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Heni, A.C., Fackelmann, G., Eibner, G. et al. Wildlife gut microbiomes of sympatric generalist species respond differently to anthropogenic landscape disturbances. anim microbiome 5, 22 (2023). https://doi.org/10.1186/s42523-023-00237-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s42523-023-00237-9

Keywords