Inhibitory Bacterial Diversity and Mucosome Function Differentiate Susceptibility of Appalachian Salamanders to Chytrid Fungal Infection

Invasibility of a North American soil ecosystem to amphibian-killing fungal pathogens

North American salamanders are threatened by intercontinental spread of chytridiomycosis, a deadly disease caused by the fungal pathogen Batrachochytrium salamandrivorans (Bsal). To predict potential dispersal of Bsal spores to salamander habitats, we evaluated the capacity of soil microbial communities to resist invasion. We determined the degree of habitat invasibility using soils from five locations throughout the Great Smoky Mountains National Park, a region with a high abundance of susceptible hosts. Our experimental design consisted of replicate soil microcosms exposed to different propagule pressures of the non-native pathogen, Bsal, and an introduced but endemic pathogen, B. dendrobatidis (Bd). To compare growth and competitive interactions, we used quantitative PCR, live/dead cell viability assays, and full-length 16S rRNA sequencing. We found that soil microcosms with intact bacterial communities inhibited both Bsal and Bd growth, but inhibitory capacity diminished with increased propagule pressure. Bsal showed greater persistence than Bd. Linear discriminant analysis (LDA) identified the family Burkolderiaceae as increasing in relative abundance with the decline of both pathogens. Although our findings provide evidence of environmental filtering in soils, such barriers weakened in response to pathogen type and propagule pressure, showing that habitats vary their invasibility based on properties of their local microbial communities.

Integrating microbiome science and evolutionary medicine into animal health and conservation

Microbiome science has provided groundbreaking insights into human and animal health. Similarly, evolutionary medicine - the incorporation of eco-evolutionary concepts into primarily human medical theory and practice - is increasingly recognised for its novel perspectives on modern diseases. Studies of host-microbe relationships have been expanded beyond humans to include a wide range of animal taxa, adding new facets to our understanding of animal ecology, evolution, behaviour, and health. In this review, we propose that a broader application of evolutionary medicine, combined with microbiome science, can provide valuable and innovative perspectives on animal care and conservation. First, we draw on classic ecological principles, such as alternative stable states, to propose an eco-evolutionary framework for understanding variation in animal microbiomes and their role in animal health and wellbeing. With a focus on mammalian gut microbiomes, we apply this framework to populations of animals under human care, with particular relevance to the many animal species that suffer diseases linked to gut microbial dysfunction (e.g. gut distress and infection, autoimmune disorders, obesity). We discuss diet and microbial landscapes (i.e. the microbes in the animal's external environment), as two factors that are (i) proposed to represent evolutionary mismatches for captive animals, (ii) linked to gut microbiome structure and function, and (iii) potentially best understood from an evolutionary medicine perspective. Keeping within our evolutionary framework, we highlight the potential benefits - and pitfalls - of modern microbial therapies, such as pre- and probiotics, faecal microbiota transplants, and microbial rewilding. We discuss the limited, yet growing, empirical evidence for the use of microbial therapies to modulate animal gut microbiomes beneficially. Interspersed throughout, we propose 12 actionable steps, grounded in evolutionary medicine, that can be applied to practical animal care and management. We encourage that these actionable steps be paired with integration of eco-evolutionary perspectives into our definitions of appropriate animal care standards. The evolutionary perspectives proposed herein may be best appreciated when applied to the broad diversity of species under human care, rather than when solely focused on humans. We urge animal care professionals, veterinarians, nutritionists, scientists, and others to collaborate on these efforts, allowing for simultaneous care of animal patients and the generation of valuable empirical data.

Phylosymbiosis shapes skin bacterial communities and pathogen-protective function in Appalachian salamanders

Phylosymbiosis is an association between host-associated microbiome composition and host phylogeny. This pattern can arise via the evolution of host traits, habitat preferences, diets, and the co-diversification of hosts and microbes. Understanding the drivers of phylosymbiosis is vital for modelling disease-microbiome interactions and manipulating microbiomes in multi-host systems. This study quantifies phylosymbiosis in Appalachian salamander skin in the context of infection by the fungal pathogen Batrachochytrium dendrobatidis (Bd), while accounting for environmental microbiome exposure. We sampled ten salamander species representing >150M years of divergence, assessed their Bd infection status, and analysed their skin and environmental microbiomes. Our results reveal a significant signal of phylosymbiosis, whereas the local environmental pool of microbes, climate, geography, and Bd infection load had a smaller impact. Host-microbe co-speciation was not evident, indicating that the effect stems from the evolution of host traits influencing microbiome assembly. Bd infection is correlated with host phylogeny and the abundance of Bd-inhibitory bacterial strains, suggesting that the long-term evolutionary dynamics between salamander hosts and their skin microbiomes affect the present-day distribution of the pathogen, along with habitat-linked exposure risk. Five Bd-inhibitory bacterial strains showed unusual generalism: occurring in most host species and habitats. These generalist strains may enhance the likelihood of probiotic manipulations colonising and persisting on hosts. Our results underscore the substantial influence of host-microbiome eco-evolutionary dynamics on environmental health and disease outcomes.

Chytrid infections exhibit historical spread and contemporary seasonality in a declining stream-breeding frog

Species with extensive geographical ranges pose special challenges to assessing drivers of wildlife disease, necessitating collaborative and large-scale analyses. The imperilled foothill yellow-legged frog (Rana boylii) inhabits a wide geographical range and variable conditions in rivers of California and Oregon (USA), and is considered threatened by the pathogen Batrachochytrium dendrobatidis (Bd). To assess drivers of Bd infections over time and space, we compiled over 2000 datapoints from R. boylii museum specimens (collected 1897-2005) and field samples (2005-2021) spanning 9° of latitude. We observed a south-to-north spread of Bd detections beginning in the 1940s and increase in prevalence from the 1940s to 1970s, coinciding with extirpation from southern latitudes. We detected eight high-prevalence geographical clusters through time that span the species' geographical range. Field-sampled male R. boylii exhibited the highest prevalence, and juveniles sampled in autumn exhibited the highest loads. Bd infection risk was highest in lower elevation rain-dominated watersheds, and with cool temperatures and low stream-flow conditions at the end of the dry season. Through a holistic assessment of relationships between infection risk, geographical context and time, we identify the locations and time periods where Bd mitigation and monitoring will be critical for conservation of this imperilled species.

Glutathione is required for growth and cadmium tolerance in the amphibian chytrid fungus, Batrachochytrium dendrobatidis

Batrachochytrium dendrobatidis (Bd) is a lethal amphibian pathogen, partly due to its ability to evade the immune system of susceptible frog species. In many pathogenic fungi, the antioxidant glutathione is a virulence factor that helps neutralise oxidative stressors generated from host immune cells, as well as other environmental stressors such as heavy metals. The role of glutathione in stress tolerance in Bd has not been investigated. Here, we examine the changes in the glutathione pool after stress exposure and quantify the effect of glutathione depletion on cell growth and stress tolerance. Depletion of glutathione repressed growth and release of zoospores, suggesting that glutathione is essential for life cycle completion in Bd. Supplementation with <2 mM exogenous glutathione accelerated zoospore development, but concentrations >2 mM were strongly inhibitory to Bd cells. While hydrogen peroxide exposure lowered the total cellular glutathione levels by 42 %, glutathione depletion did not increase the sensitivity to hydrogen peroxide. Exposure to cadmium increased total cellular glutathione levels by 93 %. Glutathione-depleted cells were more sensitive to cadmium, and this effect was attenuated by glutathione supplementation, suggesting that glutathione plays an important role in cadmium tolerance. The effects of heat and salt were exacerbated by the addition of exogenous glutathione. The impact of glutathione levels on Bd stress sensitivity may help explain differences in host susceptibility to chytridiomycosis and may provide opportunities for synergistic therapeutics.

Linking microbiome and stress hormone responses in wild tropical treefrogs across continuous and fragmented forests

The amphibian skin microbiome is an important component of anti-pathogen defense, but the impact of environmental change on the link between microbiome composition and host stress remains unclear. In this study, we used radiotelemetry and host translocation to track microbiome composition and function, pathogen infection, and host stress over time across natural movement paths for the forest-associated treefrog, Boana faber. We found a negative correlation between cortisol levels and putative microbiome function for frogs translocated to forest fragments, indicating strong integration of host stress response and anti-pathogen potential of the microbiome. Additionally, we observed a capacity for resilience (resistance to structural change and functional loss) in the amphibian skin microbiome, with maintenance of putative pathogen-inhibitory function despite major temporal shifts in microbiome composition. Although microbiome community composition did not return to baseline during the study period, the rate of microbiome change indicated that forest fragmentation had more pronounced effects on microbiome composition than translocation alone. Our findings reveal associations between stress hormones and host microbiome defenses, with implications for resilience of amphibians and their associated microbes facing accelerated tropical deforestation.

Associations of Batrachochytrium dendrobatidis with skin bacteria and fungi on Asian amphibian hosts

Amphibian skin harbors microorganisms that are associated with the fungal pathogen Batrachochytrium dendrobatidis (Bd), which causes chytridiomycosis, one of the most significant wildlife diseases known. This pathogen originated in Asia, where diverse Bd lineages exist; hence, native amphibian hosts have co-existed with Bd over long time periods. Determining the nuances of this co-existence is crucial for understanding the prevalence and spread of Bd from a microbial context. However, associations of Bd with the natural skin microbiome remain poorly understood for Asian hosts, especially in relation to skin-associated fungi. We used 16 S rRNA and fungal internal transcribed spacer (ITS) gene sequencing to characterize the skin microbiome of four native Asian amphibian species and examined the relationships between Bd infection and their skin bacterial and fungal communities; we also analyzed the correlates of the putative anti-Bd bacteria. We show that both skin bacterial and fungal community structure and composition had significant associations with infection status (Bd presence/absence) and infection intensity (frequency of Bd sequence reads). We also found that the putative anti-Bd bacterial richness was correlated with Bd infection status and infection intensity, and observed that the relative abundance of anti-Bd bacteria roughly correspond with changes in both Bd prevalence and mean infection intensity in populations. Additionally, the microbial co-occurrence network of infected frogs was significantly different from that of uninfected frogs that were characterized by more keystone nodes (connectors) and larger proportions in correlations between bacteria, suggesting stronger inter-module bacterial interactions. These results indicate that the mutual effects between Bd and skin-associated microbiome, including the interplay between bacteria and fungi, might vary with Bd infection in susceptible amphibian species. This knowledge will help in understanding the dynamics of Bd from a microbial perspective, potentially contributing to mitigate chytridiomycosis in other regions of the world.

Bacterial communities of the threatened Western Pond Turtle may be impacted by land use

As semi-aquatic species that use both terrestrial and aquatic habitats, freshwater turtles and their microbial communities are especially sensitive to the impacts of habitat disturbance. In this study, we use 16S rRNA amplicon sequencing to characterize the shell and cloacal bacterial communities of turtles in the San Francisco Bay Area. We captured western pond turtles (Actinemys/Emys marmorata) across eight sites located in urban and rural environments, along with invasive red-eared sliders (Trachemys scripta elegans). We assessed differences in western pond turtle bacterial communities diversity/composition between shell and cloacal samples and evaluated how alpha/beta diversity metrics were influenced by habitat quality. We found phylum-level bacterial taxonomic turnover in the bacterial communities of western pond turtles relative to the host tissue substrate samples. Our findings indicate that location identity elicits a high degree of lower-level (i.e. species/genus) bacterial taxonomic turnover. Further, we found that samples originating from good quality habitat had poorer shell bacterial communities but more diverse cloacal ones. The shell bacterial communities of red-eared sliders overlapped with those western pond turtles suggesting the existence of microbial dispersal between these two species. Our results add to our current understanding of turtle symbiont microbial ecology by establishing patterns of bacterial symbiont variation in an urban to rural gradient.

Host phylogeny and environment shape the diversity of salamander skin bacterial communities

The composition and diversity of animal-associated microbial communities are shaped by multiple ecological and evolutionary processes acting at different spatial and temporal scales. Skin microbiomes are thought to be strongly influenced by the environment due to the direct interaction of the host's skin with the external media. As expected, the diversity of amphibian skin microbiomes is shaped by climate and host sampling habitats, whereas phylogenetic effects appear to be weak. However, the relative strength of phylogenetic and environmental effects on salamander skin microbiomes remains poorly understood. Here, we analysed sequence data from 1164 adult salamanders of 44 species to characterise and compare the diversity and composition of skin bacteria. We assessed the relative contribution of climate, host sampling habitat, and host phylogeny to the observed patterns of bacterial diversity. We found that bacterial alpha diversity was mainly associated with host sampling habitat and climate, but that bacterial beta diversity was more strongly associated with host taxonomy and phylogeny. This phylogenetic effect predominantly occurred at intermediate levels of host divergence (0-50 Mya). Our results support the importance of environmental factors shaping the diversity of salamander skin microbiota, but also support host phylogenetic history as a major factor shaping these bacterial communities.

The adaptive microbiome hypothesis and immune interactions in amphibian mucus

The microbiome is known to provide benefits to hosts, including extension of immune function. Amphibians are a powerful immunological model for examining mucosal defenses because of an accessible epithelial mucosome throughout their developmental trajectory, their responsiveness to experimental treatments, and direct interactions with emerging infectious pathogens. We review amphibian skin mucus components and describe the adaptive microbiome as a novel process of disease resilience where competitive microbial interactions couple with host immune responses to select for functions beneficial to the host. We demonstrate microbiome diversity, specificity of function, and mechanisms for memory characteristic of an adaptive immune response. At a time when industrialization has been linked to losses in microbiota important for host health, applications of microbial therapies such as probiotics may contribute to immunotherapeutics and to conservation efforts for species currently threatened by emerging diseases.

Major Emerging Fungal Diseases of Reptiles and Amphibians

Emerging infectious diseases (EIDs) are caused by pathogens that have undergone recent changes in terms of geographic spread, increasing incidence, or expanding host range. In this narrative review, we describe three important fungal EIDs with keratin trophism that are relevant to reptile and amphibian conservation and veterinary practice. Nannizziopsis spp. have been mainly described in saurians; infection results in thickened, discolored skin crusting, with eventual progression to deep tissues. Previously only reported in captive populations, it was first described in wild animals in Australia in 2020. Ophidiomyces ophidiicola (formely O. ophiodiicola) is only known to infect snakes; clinical signs include ulcerating lesions in the cranial, ventral, and pericloacal regions. It has been associated with mortality events in wild populations in North America. Batrachochytrium spp. cause ulceration, hyperkeratosis, and erythema in amphibians. They are a major cause of catastrophic amphibian declines worldwide. In general, infection and clinical course are determined by host-related characteristics (e.g., nutritional, metabolic, and immune status), pathogens (e.g., virulence and environmental survival), and environment (e.g., temperature, hygrometry, and water quality). The animal trade is thought to be an important cause of worldwide spread, with global modifications in temperature, hygrometry, and water quality further affecting fungal pathogenicity and host immune response.

Immune priming prior to pathogen exposure sheds light on the relationship between host, microbiome and pathogen in disease

Dynamic interactions between host, pathogen and host-associated microbiome dictate infection outcomes. Pathogens including Batrachochytrium dendrobatidis (Bd) threaten global biodiversity, but conservation efforts are hindered by limited understanding of amphibian host, Bd and microbiome interactions. We conducted a vaccination and infection experiment using Eastern hellbender salamanders (Cryptobranchus alleganiensis alleganiensis) challenged with Bd to observe infection, skin microbial communities and gene expression of host skin, pathogen and microbiome throughout the experiment. Most animals survived high Bd loads regardless of their vaccination status and vaccination did not affect pathogen load, but host gene expression differed based on vaccination. Oral vaccination (exposure to killed Bd) stimulated immune gene upregulation while topically and sham-vaccinated animals did not significantly upregulate immune genes. In early infection, topically vaccinated animals upregulated immune genes but orally and sham-vaccinated animals downregulated immune genes. Bd increased pathogenicity-associated gene expression in late infection when Bd loads were highest. The microbiome was altered by Bd, but there was no correlation between anti-Bd microbe abundance or richness and pathogen burden. Our observations suggest that hellbenders initially generate a vigorous immune response to Bd, which is ineffective at controlling disease and is subsequently modulated. Interactions with antifungal skin microbiota did not influence disease progression.

Spatiotemporal and ontogenetic variation, microbial selection, and predicted Bd-inhibitory function in the skin-associated microbiome of a Rocky Mountain amphibian

Host-associated microbiomes play important roles in host health and pathogen defense. In amphibians, the skin-associated microbiota can contribute to innate immunity with potential implications for disease management. Few studies have examined season-long temporal variation in the amphibian skin-associated microbiome, and the interactions between bacteria and fungi on amphibian skin remain poorly understood. We characterize season-long temporal variation in the skin-associated microbiome of the western tiger salamander (Ambystoma mavortium) for both bacteria and fungi between sites and across salamander life stages. Two hundred seven skin-associated microbiome samples were collected from salamanders at two Rocky Mountain lakes throughout the summer and fall of 2018, and 127 additional microbiome samples were collected from lake water and lake substrate. We used 16S rRNA and ITS amplicon sequencing with Bayesian Dirichlet-multinomial regression to estimate the relative abundances of bacterial and fungal taxa, test for differential abundance, examine microbial selection, and derive alpha diversity. We predicted the ability of bacterial communities to inhibit the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd), a cutaneous fungal pathogen, using stochastic character mapping and a database of Bd-inhibitory bacterial isolates. For both bacteria and fungi, we observed variation in community composition through time, between sites, and with salamander age and life stage. We further found that temporal trends in community composition were specific to each combination of salamander age, life stage, and lake. We found salamander skin to be selective for microbes, with many taxa disproportionately represented relative to the environment. Salamander skin appeared to select for predicted Bd-inhibitory bacteria, and we found a negative relationship between the relative abundances of predicted Bd-inhibitory bacteria and Bd. We hope these findings will assist in the conservation of amphibian species threatened by chytridiomycosis and other emerging diseases.

The skin microbiota of the axolotl Ambystoma altamirani is highly influenced by metamorphosis and seasonality but not by pathogen infection

BACKGROUND

Microbiomes have been increasingly recognized as major contributors to host health and survival. In amphibians, bacterial members of the skin microbiota protect their hosts by inhibiting the growth of the fungal pathogen Batrachochytrium dendrobatidis (Bd). Even though several studies describe the influence of biotic and abiotic factors over the skin microbiota, it remains unclear how these symbiotic bacterial communities vary across time and development. This is particularly relevant for species that undergo metamorphosis as it has been shown that host physiology and ecology drastically influence diversity of the skin microbiome.

RESULTS

We found that the skin bacterial communities of the axolotl A. altamirani are largely influenced by the metamorphic status of the host and by seasonal variation of abiotic factors such as temperature, pH, dissolved oxygen and conductivity. Despite high Bd prevalence in these samples, the bacterial diversity of the skin microbiota did not differ between infected and non-infected axolotls, although relative abundance of particular bacteria were correlated with Bd infection intensity.

CONCLUSIONS

Our work shows that metamorphosis is a crucial process that shapes skin bacterial communities and that axolotls under different developmental stages respond differently to environmental seasonal variations. Moreover, this study greatly contributes to a better understanding of the factors that shape amphibian skin microbiota, especially in a largely underexplored group like axolotls (Mexican Ambystoma species).

Bacterial Skin Assemblages of Sympatric Salamanders Are Primarily Shaped by Host Genus

Bacterial assemblages on the skins of amphibians are known to influence pathogen resistance and other important physiological functions in the host. Host-specific factors and the environment play significant roles in structuring skin assemblages. This study used high-throughput 16S rRNA sequencing and multivariate analyses to examine differences in skin-bacterial assemblages from 246 salamanders belonging to three genera in the lungless family Plethodontidae along multiple spatial gradients. Composition and α- and β-diversity of bacterial assemblages were defined, indicator species were identified for each host group, and the relative influences of host- versus environment-specific ecological factors were evaluated. At the broadest spatial scale, host genus, host species, and sampling site were predictive of skin assemblage structure, but host genus and species were more influential after controlling for the marginal effects of site, as well as nestedness of site. Furthermore, assemblage similarity within each host genus did not change with increasing geographic distance. At the smallest spatial scale, site-specific climate analyses revealed different relationships to climatic variables for each of the three genera, and these relationships were determined by host ecomode. Variation in bacterial assemblages of terrestrial hosts correlated with landscape-level climatic variability, and this pattern decayed with increasing water dependence of the host. Results from this study highlight host-specific considerations for researchers studying wildlife diseases in co-occurring, yet ecologically divergent, species.