Reduced skin bacterial diversity correlates with increased pathogen infection intensity in an endangered amphibian host

Gene functions of the Ambystoma altamirani skin microbiome vary across space and time but potential antifungal genes are widespread and prevalent

Amphibian skin microbiomes can play a critical role in host survival against emerging diseases by protecting their host against pathogens. While a plethora of biotic and abiotic factors have been shown to influence the taxonomic diversity of amphibian skin microbiomes it remains unclear whether functional genomic diversity varies in response to temporal and environmental factors. Here we applied a metagenomic approach to evaluate whether seasonality, distinct elevations/sites, and pathogen presence influenced the functional genomic diversity of the A. altamirani skin microbiome. We obtained a gene catalogue of 92 107 nonredundant annotated genes and a set of 50 unique metagenome assembled genomes (MAGs). Our analysis showed that genes linked to general and potential antifungal traits significantly differed across seasons and sampling locations at different elevations. Moreover, we found that the functional genomic diversity of A. altamirani skin microbiome differed between B. dendrobatidis infected and not infected axolotls only during winter, suggesting an interaction between seasonality and pathogen infection. In addition, we identified the presence of genes and biosynthetic gene clusters (BGCs) linked to potential antifungal functions such as biofilm formation, quorum sensing, secretion systems, secondary metabolite biosynthesis, and chitin degradation. Interestingly genes linked to these potential antifungal traits were mainly identified in Burkholderiales and Chitinophagales MAGs. Overall, our results identified functional traits linked to potential antifungal functions in the A. altamirani skin microbiome regardless of variation in the functional diversity across seasons, elevations/sites, and pathogen presence. Our findings suggest that potential antifungal traits found in Burkholderiales and Chitinophagales taxa could be related to the capacity of A. altamirani to survive in the presence of Bd, although further experimental analyses are required to test this hypothesis.

Association between the skin microbiome and MHC class II diversity in an amphibian

Microbiomes play an important role in determining the ecology and behaviour of their hosts. However, questions remain pertaining to how host genetics shape microbiomes, and how microbiome composition influences host fitness. We explored the effects of geography, evolutionary history and host genetics on the skin microbiome diversity and structure in a widespread amphibian. More specifically, we examined the association between bacterial diversity and composition and the major histocompatibility complex class II exon 2 diversity in 12 moor frog (Rana arvalis) populations belonging to two geographical clusters that show signatures of past and ongoing differential selection. We found that while bacterial alpha diversity did not differ between the two clusters, MHC alleles/supertypes and genetic diversity varied considerably depending on geography and evolutionary history. Bacterial alpha diversity was positively correlated with expected MHC heterozygosity and negatively with MHC nucleotide diversity. Furthermore, bacterial community composition showed significant variation between the two geographical clusters and between specific MHC alleles/supertypes. Our findings emphasize the importance of historical demographic events on hologenomic variation and provide new insights into how immunogenetic host variability and microbial diversity may jointly influence host fitness with consequences for disease susceptibility and population persistence.

Low-level pathogen infection and geographic location correlate with the skin microbiomes of Columbia spotted frogs (Rana luteiventris) in a montane landscape

The skin microbiome of amphibians can influence host susceptibility towards the fungal pathogen Batrachochytrium dendrobatidis (Bd), while simultaneously having the potential to be altered by Bd. Severe Bd infections are known to alter the amphibian skin microbiome; however, little is known about microbiome interactions in amphibians with low infection intensity. In addition to disease dynamics, environmental factors may influence the microbiome. To test for patterns in bacterial diversity based on pathogen infection and environmental factors, 399 Columbia spotted frogs (Rana luteiventris) were sampled throughout northern Idaho and northeastern Washington across two years. Bd prevalence and intensity were measured in 376 frogs, revealing a prevalence of 69%, but generally low infection intensity (Mean = 127 Bd zoospore equivalents among infected frogs). Skin bacterial communities were characterized in 92 frogs using 16S rRNA gene amplicon sequencing. Our results indicated correlations of decreasing Shannon diversity and evenness as infection intensity increased. Latitude was correlated with bacterial richness and Faith's Phylogenetic Diversity measures, indicating increased diversity in northern locations. Beta diversity (UniFrac) analyses revealed that skin microbiomes were distinct between infected and uninfected frogs, and infection intensity had a significant effect on microbiome composition. Site explained the majority of microbiome variation (weighted UniFrac: 57.5%), suggesting a combination of local habitat conditions explain variation, as only small proportions of variation could be explained by year, month, temperature, elevation, and latitude individually. Bacterial genera with potential for Bd-inhibitory properties were found with differential relative abundance in infected and uninfected frogs, with higher Stenotrophomonas and lower Pseudomonas relative abundance observed in infected frogs. Further study may indicate if Bd inhibition by members of the skin microbiome is an influence behind the low infection intensities observed and whether low Bd infection intensities are capable of altering skin microbiome composition.

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.

Energy-related wastewater contamination alters microbial communities of sediment, water, and amphibian skin

To inform responsible energy development, it is important to understand the ecological effects of contamination events. Wastewaters, a common byproduct of oil and gas extraction, often contain high concentrations of sodium chloride (NaCl) and heavy metals (e.g., strontium and vanadium). These constituents can negatively affect aquatic organisms, but there is scarce information for how wastewaters influence potentially distinct microbiomes in wetland ecosystems. Additionally, few studies have concomitantly investigated effects of wastewaters on the habitat (water and sediment) and skin microbiomes of amphibians or relationships among these microbial communities. We sampled microbiomes of water, sediment, and skin of four larval amphibian species across a gradient of chloride contamination (0.04-17,500 mg/L Cl) in the Prairie Pothole Region of North America. We detected 3129 genetic phylotypes and 68 % of those phylotypes were shared among the three sample types. The most common shared phylotypes were Proteobacteria, Firmicutes, and Bacteroidetes. Salinity of wastewaters increased dissimilarity within all three microbial communities, but not the diversity or richness of water and skin microbial communities. Strontium was associated with lower diversity and richness of sediment microbial communities, but not those of water or amphibian skin, likely because metal deposition occurs in sediment when wetlands dry. Based on Bray Curtis distance matrices, sediment microbiomes were similar to those of water, but neither had substantial overlap with amphibian microbiomes. Species identity was the strongest predictor of amphibian microbiomes; frog microbiomes were similar but differed from that of the salamander, whose microbiome had the lowest richness and diversity. Understanding how effects of wastewaters on the dissimilarity, richness, and diversity of microbial communities also influence the ecosystem function of communities will be an important next step. However, our study provides novel insight into the characteristics of, and associations among, different wetland microbial communities and effects of wastewaters from energy production.

Fine cassava fibre utilization as a dietary fibre source for dogs: Effects on kibble characteristics, diet digestibility and palatability, faecal metabolites and microbiota

The objective was to evaluate through three experiments the effects of a fine cassava fibre (CA: 106 µm) on kibble characteristics, coefficients of total tract apparent digestibility (CTTAD) of macronutrients, diet palatability and faecal metabolites and microbiota of dogs. Dietary treatments consisted of a control diet (CO), without an added fibre source and with 4.3% total dietary fibre (TDF), and a diet with 9.6% CA (106 µm), with 8.4% TDF. Experiment I evaluated the physical characteristics of the kibbles. The palatability test was evaluated in experiment II, which compared the diets CO versus CA. In experiment III, 12 adult dogs were randomly assigned to one of the two dietary treatments for 15 days, totalling six replicates/treatment, to assess the CTTAD of macronutrients; faecal characteristics, faecal metabolites and microbiota. The expansion index, kibble size and friability of diets with CA were higher than the CO (p < 0.05). Additionally, the CA diet presented higher palatability than the CO (p < 0.05) but did not affect CTTAD except for those of fibre (p > 0.05). Moreover, a greater faecal concentration of acetate, butyrate and total short-chain fatty acids (SCFA) and a lower faecal concentration of phenol, indole and isobutyrate were observed in dogs fed the CA diet (p < 0.05). Dogs fed with the CA diet presented a greater bacterial diversity and richness and a greater abundance of genera considered to be beneficial for gut health, such as Blautia, Faecalibacterium and Fusobacterium when compared to the CO group (p < 0.05). The inclusion of 9.6% of a fine CA improves the expansion of kibbles and diet palatability without affecting most of the CTTAD of nutrients. Besides, it improves the production of some SCFA and modulates the faecal microbiota of dogs.

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).

Environment is associated with chytrid infection and skin microbiome richness on an amphibian rich island (Taiwan)

Growing evidence suggests that the origins of the panzootic amphibian pathogens Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal) are in Asia. In Taiwan, an island hotspot of high amphibian diversity, no amphibian mass mortality events linked to Bd or Bsal have been reported. We conducted a multi-year study across this subtropical island, sampling 2517 individuals from 30 species at 34 field sites, between 2010 and 2017, and including 171 museum samples collected between 1981 and 2009. We analyzed the skin microbiome of 153 samples (6 species) from 2017 in order to assess any association between the amphibian skin microbiome and the probability of infection amongst different host species. We did not detect Bsal in our samples, but found widespread infection by Bd across central and northern Taiwan, both taxonomically and spatially. Museum samples show that Bd has been present in Taiwan since at least 1990. Host species, geography (elevation), climatic conditions and microbial richness were all associated with the prevalence of infection. Host life-history traits, skin microbiome composition and phylogeny were associated with lower prevalence of infection for high altitude species. Overall, we observed low prevalence and burden of infection in host populations, suggesting that Bd is enzootic in Taiwan where it causes subclinical infections. While amphibian species in Taiwan are currently threatened by habitat loss, our study indicates that Bd is in an endemic equilibrium with the populations and species we investigated. However, ongoing surveillance of the infection is warranted, as changing environmental conditions may disturb the currently stable equilibrium.

Perturbation by Antimicrobial Bacteria of the Epidermal Bacterial Flora of Rainbow Trout in Flow-Through Aquaculture

The bacterial flora of the epidermal mucus of fish is closely associated with the host’s health and susceptibility to pathogenic infections. In this study, we analyzed the epidermal mucus bacteria of rainbow trout (Oncorhynchus mykiss) reared in flow-through aquaculture under environmental perturbations. Over ~2 years, the bacteria present in the skin mucus and water were analyzed based on the 16S rDNA sequences. The composition of the mucus bacterial community showed significant monthly fluctuations, with frequent changes in the dominant bacterial species. Analysis of the beta- and alpha-diversity of the mucus bacterial flora showed the fluctuations of the composition of the flora were caused by the genera Pseudomonas, Yersinia, and Flavobacterium, and some species of Pseudomonas and Yersinia in the mucus were identified as antimicrobial bacteria. Examination of the antimicrobial bacteria in the lab aquarium showed that the natural presence of antimicrobial bacteria in the mucus and water, or the purposeful addition of them to the rearing water, caused a transition in the mucus bacteria community composition. These results demonstrate that specific antimicrobial bacteria in the water or in epidermal mucus comprise one of the causes of changes in fish epidermal mucus microflora.

Genomics of Serrasalmidae teleosts through the lens of microbiome fingerprinting

Associations between host genotype and host-associated microbiomes have been shown in a variety of animal clades, but studies on teleosts mostly show weak associations. Our study aimed to explore these relationships in four sympatric Serrasalmidae (i.e. piranha) teleosts from an Amazonian lake, using datasets from the hosts genomes (SNPs from GBS), skin and gut microbiomes (16S rRNA gene metataxonomics), and diets (COI metabarcoding) from the same fish individuals. Firstly, we investigated whether there were significant covariations of microbiome and fish genotypes at the inter and intraspecific levels. We also assessed the extent of co-variation between Serrasalmidae diet and microbiome, to isolate genotypic from dietary effects on community structure. We observed a significant covariation of skin microbiomes and host genotypes at interspecific (R² =24.4%) and intraspecific (R² =6.2%) levels, whereas gut microbiomes correlated poorly with host genotypes. Serrasalmidae diet composition was significantly correlated to fish genotype only at the interspecific level (R² =5.4%), but did not covary with gut microbiome composition (mantel R=-0.04). Secondly, we investigated whether the study of interspecific differentiation could benefit from considering host associated microbial communities in addition to host genotypes. By using a NMDS ordination-based approach, we observed that ordinations from skin and gut species-specific bacterial biomarkers identified through a random forest algorithm, could significantly increase the average interspecific differentiation detected through host genotype data alone. Although future studies encompassing additional species and environments are needed, our results suggest Serrasalmidae microbiomes could constitute an insightful trait to be considered when studying the interspecific differences between members of this clade.

Factors Influencing Bacterial and Fungal Skin Communities of Montane Salamanders of Central Mexico

Host microbial communities are increasingly seen as an important component of host health. In amphibians, the first land vertebrates that are threatened by a fungal skin disease globally, our understanding of the factors influencing the microbiome of amphibian skin remains incomplete because recent studies have focused almost exclusively on bacteria, and little information exists on fungal communities associated with wild amphibian species. In this study, we describe the effects of host phylogeny, climate, geographic distance, and infection with a fungal pathogen on the composition and structure of bacterial and fungal communities in seven tropical salamander species that occur in the Trans-Mexican Volcanic Belt of Central Mexico. We find that host phylogenetic relatedness is correlated with bacterial community composition while a composite climatic variable of temperature seasonality and precipitation is significantly associated with fungal community composition. We also estimated co-occurrence networks for bacterial and fungal taxa and found differences in the degree of connectivity and the distribution of negative associations between the two networks. Our results suggest that different factors may be responsible for structuring the bacterial and fungal communities of amphibian skin and that the inclusion of fungi in future studies could shed light on important functional interactions within the microbiome.

Predicting fungal infection rate and severity with skin-associated microbial communities on amphibians

Pathogen success (risk and severity) is influenced by host-associated microbiota, but the degree to which variation in microbial community traits predict future infection presence/absence (risk) and load (severity) for the host is unknown. We conducted a time-series experiment by sampling the skin-associated bacterial communities of five amphibian species before and after exposure to the fungal pathogen, Batrachochytrium dendrobaditis (Bd). We ask whether microbial community traits are predictors of, or are affected by, Bd infection risk and intensity. Our results show that richness of putative Bd-inhibitory bacteria strongly predicts infection risk, while the proportion of putative Bd-inhibitory bacteria predicts future infection intensity. Variation in microbial community composition is high across time and individual, and bacterial prevalence is low. Our findings demonstrate how ecological community traits of host-associated microbiota may be used to predict infection risk by pathogenic microbes.

Effect of multiple stressors and population decline of frogs

The ongoing decline in anuran populations is linked primarily to the effects of stressor agents such as pathogens, pesticides, alterations of natural landscapes, and the introduction of exotic species. Most studies that have evaluated the effects of these stressors have focused on a single component, which is the opposite of the reality of most natural environments, where anuran populations tend to suffer the influence of multiple agents simultaneously. Studies of the effects of the interaction between these components are extremely important, given that one agent may potentialize (synergistic effect) or weaken another (antagonistic effect) or, in some cases, have a neutral effect. The present study is based on the scientometric analysis of three bibliographic databases (ISI Web of Science, Scopus, and PubMed), which identified 1376 papers that reported on the global decline of anuran populations, although only 172 of these studies focused on the interactive effects of environmental stressors. Synergistic effects were the most frequent type of interaction, followed by antagonistic effects, and a small number of studies that found no clear interaction between the stressors. Pathogens and pesticides were the classes of stressor studied most frequently, while climate-pathogen and pathogen-pesticide interactions were the combinations that featured in the largest number of studies. Overall, we would recommend a more systematic focus on the dynamics of the interactions among the stressors that impact anuran populations, in particular for the elaboration of conservation programs, given that these agents tend to have complex combined effects.

Microbiome associated with the tetrodotoxin-bearing anuran Brachycephalus pitanga

The genus Brachycephalus includes small species of aposematic anurans known as microendemic, occurring in the mountains of the Atlantic Forest. Brachycephalus ephippium, B. nodoterga and B. pernix have been reported to contain the neurotoxin tetrodotoxin in skin and viscera. The biological conservation of several Brachycephalus species is currently threatened by climate change, deforestation, and the pandemic caused by the fungus Batrachochytrium dendrobatidis (Bd). Despite the well-known importance of amphibians' associated bacteria in the defensive role against pathogens, there is still a poor understanding of amphibian microbiome composition. The present study investigated the composition of B. pitanga microbial community and the presence of TTX in the host and in cultures of bacterial isolates, using a combination of metagenomics, bacterial culture isolation, mass spectrometry and metabolomic analyses. Results of culture-dependent and -independent analyses characterized the microbial communities associated with the skin and viscera of B. pitanga. Mass spectrometry analysis indicated the presence of TTX in host tissues, while bacterial production of TTX was not observed under the experimental conditions used in this investigation. This is the first report confirming the occurrence of TTX in B. pitanga.

Longitudinal patterns in the skin microbiome of wild, individually marked frogs from the Sierra Nevada, California

The amphibian skin microbiome has been the focus of numerous studies because of the protective effects that some bacteria provide against the pathogenic fungus Batrachochytrium dendrobatidis, which has caused a global panzootic among amphibians. However, the mechanisms driving community structure and function in the amphibian skin microbiome are still poorly understood, and longitudinal analyses of the skin microbiome have not yet been conducted in wild populations. In this study, we investigate longitudinal patterns in the skin microbiome of 19 individually marked adult frogs from two wild populations of the endangered Sierra Nevada yellow-legged frog (Rana sierrae), sampled over the course of 2 years. We found that individuals with low bacterial diversity (dominated by order Burkhorderiales) had significantly more stable bacterial communities than those with higher diversity. Amplicon sequence variants (ASVs) with high relative abundance were significantly less transient than those with low relative abundance, and ASVs with intermediate-level relative abundances experienced the greatest volatility over time. Based on these results, we suggest that efforts to develop probiotic treatments to combat B. dendrobatidis should focus on bacteria that are found at high relative abundances in some members of a population, as these strains are more likely to persist and remain stable in the long term.

Assembly of the amphibian microbiome is influenced by the effects of land-use change on environmental reservoirs

A growing focus in microbial ecology is understanding how beneficial microbiome function is created and maintained through various assembly mechanisms. This study explores the role of both the environment and disease in regulating the composition of microbial species in the soil and on amphibian hosts. We compared the microbial communities of Plethodon cinereus salamanders along a land-use gradient in the New York metropolitan area and paired these with associated soil cores. Additionally, we characterized the diversity of bacterial and fungal symbionts that putatively inhibit the pathogenic fungus Batrachochytrium dendrobatidis. We predicted that variation in skin microbial community composition would correlate with changes seen in the soil which functions as the regional species pool. We found that salamanders and soil share many microbial taxa but that these two communities exhibit differences in the relative abundances of the bacterial phyla Acidobacteria, Actinobacteria, and Proteobacteria and the fungal phyla Ascomycota and genus Basidiobolus. Microbial community composition varies with changes in land-use associated factors creating site-specific compositions. By employing a quantitative, null-based assembly model, we identified that dispersal limitation, variable selection, and drift guide assembly of microbes onto their skin, creating high dissimilarity between individuals with likely consequences in disease preventative function.

Diet-Microbiota Interactions Alter Mosquito Development

Gut microbes and diet can both strongly affect the biology of multicellular animals, but it is often difficult to disentangle microbiota–diet interactions due to the complex microbial communities many animals harbor and the nutritionally variable diets they consume. While theoretical and empirical studies indicate that greater microbiota diversity is beneficial for many animal hosts, there have been few tests performed in aquatic invertebrates. Most mosquito species are aquatic detritivores during their juvenile stages that harbor variable microbiotas and consume diets that range from nutrient rich to nutrient poor. In this study, we produced a gnotobiotic model that allowed us to examine how interactions between specific gut microbes and diets affect the fitness of Aedes aegypti, the yellow fever mosquito. Using a simplified seven-member community of bacteria (ALL7) and various laboratory and natural mosquito diets, we allowed larval mosquitoes to develop under different microbial and dietary conditions and measured the resulting time to adulthood and adult size. Larvae inoculated with the ALL7 or a more complex community developed similarly when fed nutrient-rich rat chow or fish food laboratory diets, whereas larvae inoculated with individual bacterial members of the ALL7 community exhibited few differences in development when fed a rat chow diet but exhibited large differences in performance when fed a fish food diet. In contrast, the ALL7 community largely failed to support the growth of larvae fed field-collected detritus diets unless supplemented with additional protein or yeast. Collectively, our results indicate that mosquito development and fitness are strongly contingent on both diet and microbial community composition.

Co-infecting pathogen lineages have additive effects on host bacterial communities

Amphibian skin bacteria may confer protection against the fungus Batrachochytrium dendrobatidis (Bd), but responses of skin bacteria to different Bd lineages are poorly understood. The global panzootic lineage (Bd-GPL) has caused amphibian declines and extinctions globally. However, other lineages are enzootic (Bd-Asia-2/Brazil). Increased contact rates between Bd-GPL and enzootic lineages via globalization pose unknown consequences for host-microbiome-pathogen dynamics. We conducted a laboratory experiment and used 16S rRNA amplicon-sequencing to assess: (i) whether two lineages (Bd-Asia-2/Brazil and Bd-GPL) and their recombinant, in single and mixed infections, differentially affect amphibian skin bacteria; (ii) and the changes associated with the transition to laboratory conditions. We determined no clear differences in bacterial diversity among Bd treatments, despite differences in infection intensity. However, we observed an additive effect of mixed infections on bacterial alpha diversity and a potentially antagonistic interaction between Bd genotypes. Additionally, observed changes in community composition suggest a higher ability of Bd-GPL to alter skin bacteria. Lastly, we observed a drastic reduction in bacterial diversity and a change in community structure in laboratory conditions. We provide evidence for complex interactions between Bd genotypes and amphibian skin bacteria during coinfections, and expand on the implications of experimental conditions in ecological studies.

The skin microbiome of Xenopus laevis and the effects of husbandry conditions

BACKGROUND

Historically the main source of laboratory Xenopus laevis was the environment. The increase in genetically altered animals and evolving governmental constraints around using wild-caught animals for research has led to the establishment of resource centres that supply animals and reagents worldwide, such as the European Xenopus Resource Centre. In the last decade, centres were encouraged to keep animals in a "low microbial load" or "clean" state, where embryos are surface sterilized before entering the housing system; instead of the conventional, "standard" conditions where frogs and embryos are kept without prior surface treatment. Despite Xenopus laevis having been kept in captivity for almost a century, surprisingly little is known about the frogs as a holobiont and how changing the microbiome may affect resistance to disease. This study examines how the different treatment conditions, "clean" and "standard" husbandry in recirculating housing, affects the skin microbiome of tadpoles and female adults. This is particularly important when considering the potential for poor welfare caused by a change in husbandry method as animals move from resource centres to smaller research colonies.

RESULTS

We found strong evidence for developmental control of the surface microbiome on Xenopus laevis; adults had extremely similar microbial communities independent of their housing, while both tadpole and environmental microbiome communities were less resilient and showed greater diversity.

CONCLUSIONS

Our findings suggest that the adult Xenopus laevis microbiome is controlled and selected by the host. This indicates that the surface microbiome of adult Xenopus laevis is stable and defined independently of the environment in which it is housed, suggesting that the use of clean husbandry conditions poses little risk to the skin microbiome when transferring adult frogs to research laboratories. This will have important implications for frog health applicable to Xenopus laevis research centres throughout the world.

Post-epizootic microbiome associations across communities of neotropical amphibians

Microbiome-pathogen interactions are increasingly recognized as an important element of host immunity. While these host-level interactions will have consequences for community disease dynamics, the factors which influence host microbiomes at larger scales are poorly understood. We here describe landscape-scale pathogen-microbiome associations within the context of post-epizootic amphibian chytridiomycosis, a disease caused by the panzootic chytrid fungus Batrachochytrium dendrobatidis. We undertook a survey of Neotropical amphibians across altitudinal gradients in Ecuador ~30 years following the observed amphibian declines and collected skin swab-samples which were metabarcoded using both fungal (ITS-2) and bacterial (r16S) amplicons. The data revealed marked variation in patterns of both B. dendrobatidis infection and microbiome structure that are associated with host life history. Stream breeding amphibians were most likely to be infected with B. dendrobatidis. This increased probability of infection was further associated with increased abundance and diversity of non-Batrachochytrium chytrid fungi in the skin and environmental microbiome. We also show that increased alpha diversity and the relative abundance of fungi are lower in the skin microbiome of adult stream amphibians compared to adult pond-breeding amphibians, an association not seen for bacteria. Finally, stream tadpoles exhibit lower proportions of predicted protective microbial taxa than pond tadpoles, suggesting reduced biotic resistance. Our analyses show that host breeding ecology strongly shapes pathogen-microbiome associations at a landscape scale, a trait that may influence resilience in the face of emerging infectious diseases.

Composition of the North American Wood Frog (Rana sylvatica) Bacterial Skin Microbiome and Seasonal Variation in Community Structure

While a number of amphibian skin microbiomes have been characterized, it is unclear how these communities might vary in response to seasonal changes in the environment and the corresponding behaviors that many amphibians exhibit. Given recent studies demonstrating the importance of the skin microbiome in frog innate immune defense against pathogens, investigating how changes in the environment impact the microbial species present will provide a better understanding of conditions that may alter host susceptibility to pathogens in their environment. We sampled the bacterial skin microbiome of North American wood frogs (Rana sylvatica) from two breeding ponds in the spring, along with the bacterial community present in their vernal breeding pools, and frogs from the nearby forest floor in the summer and fall to determine whether community composition differs by sex, vernal pond site, or temporally across season (spring, summer, fall). Taxon relative abundance data reveals a profile of bacterial phyla similar to those previously described on anuran skin, with Proteobacteria, Bacteroidetes, and Actinobacteria dominating the wood frog skin microbiome. Our results indicate that sex had no significant effect on skin microbiota diversity; however, this may be due to our limited female frog sample size. Vernal pool site had a small but significant effect on skin microbiota, but skin-associated communities were more similar to each other than to the communities observed in the frogs' respective pond water. Across seasons, diversity analyses suggest that there are significant differences between the bacterial skin microbiome of frogs from spring and summer/fall groups while the average α-diversity per frog remained consistent. These results illustrate seasonal variation in wood frog skin microbiome structure and highlight the importance of considering temporal trends in an amphibian microbiome, particularly for species whose life history requires recurrent shifts in habitat and behavior.

Skin microbiome correlates with bioclimate and Batrachochytrium dendrobatidis infection intensity in Brazil's Atlantic Forest treefrogs

In Brazil’s Atlantic Forest (AF) biodiversity conservation is of key importance since the fungal pathogen Batrachochytrium dendrobatidis (Bd) has led to the rapid loss of amphibian populations here and worldwide. The impact of Bd on amphibians is determined by the host's immune system, of which the skin microbiome is a critical component. The richness and diversity of such cutaneous bacterial communities are known to be shaped by abiotic factors which thus may indirectly modulate host susceptibility to Bd. This study aimed to contribute to understanding the environment-host–pathogen interaction determining skin bacterial communities in 819 treefrogs (Anura: Hylidae and Phyllomedusidae) from 71 species sampled across the AF. We investigated whether abiotic factors influence the bacterial community richness and structure on the amphibian skin. We further tested for an association between skin bacterial community structure and Bd co-occurrence. Our data revealed that temperature, precipitation, and elevation consistently correlate with richness and diversity of the skin microbiome and also predict Bd infection status. Surprisingly, our data suggest a weak but significant positive correlation of Bd infection intensity and bacterial richness. We highlight the prospect of future experimental studies on the impact of changing environmental conditions associated with global change on environment-host–pathogen interactions in the AF.

Preparatory immunity: Seasonality of mucosal skin defences and Batrachochytrium infections in Southern leopard frogs

Accurately predicting the impacts of climate change on wildlife health requires a deeper understanding of seasonal rhythms in host-pathogen interactions. The amphibian pathogen, Batrachochytrium dendrobatidis (Bd), exhibits seasonality in incidence; however, the role that biological rhythms in host defences play in defining this pattern remains largely unknown. The aim of this study was to examine whether host immune and microbiome defences against Bd correspond with infection risk and seasonal fluctuations in temperature and humidity. Over the course of a year, five populations of Southern leopard frogs (Rana [Lithobates] sphenocephala) in Tennessee, United States, were surveyed for host immunity, microbiome and pathogen dynamics. Frogs were swabbed for pathogen load and skin bacterial diversity and stimulated to release stored antimicrobial peptides (AMPs). Secretions were analysed to estimate total hydrophobic peptide concentrations, presence of known AMPs and effectiveness of Bd growth inhibition in vitro. The diversity and proportion of bacterial reads with a 99% match to sequences of isolates known to inhibit Bd growth in vitro were used as an estimate of predicted anti-Bd function of the skin microbiome. Batrachochytrium dendrobatidis dynamics followed the expected seasonal fluctuations-peaks in cooler months-which coincided with when host mucosal defences were most potent against Bd. Specifically, the concentration and expression of stored AMPs cycled synchronously with Bd dynamics. Although microbiome changes followed more linear trends over time, the proportion of bacteria that can function to inhibit Bd growth was greatest when risk of Bd infection was highest. We interpret the increase in peptide storage in the fall and the shift to a more anti-Bd microbiome over winter as a preparatory response for subsequent infection risk during the colder periods when AMP synthesis and bacterial growth is slow and pathogen pressure from this cool-adapted fungus is high. Given that a decrease in stored AMP concentrations as temperatures warm in spring likely means greater secretion rates, the subsequent decrease in prevalence suggests seasonality of Bd in this host may be in part regulated by annual immune rhythms, and dominated by the effects of temperature.

Habitat disturbance influences the skin microbiome of a rediscovered neotropical-montane frog

BACKGROUND

The skin microbiome serves as a first line defense against pathogens in vertebrates. In amphibians, it has the potential to protect against the chytrid fungus Batrachochytrium dendrobatis (Bd), a likely agent of amphibian declines. Alteration of the microbiome associated with unfavorable environmental changes produced by anthropogenic activities may make the host more susceptible to pathogens. Some amphibian species that were thought to be "extinct" have been rediscovered years after population declines in the late 1980s probably due to evolved Bd-resistance and are now threatened by anthropogenic land-use changes. Understanding the effects of habitat disturbance on the host skin microbiome is relevant for understanding the health of these species, along with its susceptibility to pathogens such as Bd. Here, we investigate the influence of habitat alteration on the skin bacterial communities as well as specifically the putative Bd-inhibitory bacterial communities of the montane frog Lithobates vibicarius. This species, after years of not being observed, was rediscovered in small populations inhabiting undisturbed and disturbed landscapes, and with continuous presence of Bd.

RESULTS

We found that cutaneous bacterial communities of tadpoles and adults differed between undisturbed and disturbed habitats. The adults from disturbed habitats exhibited greater community dispersion than those from undisturbed habitats. We observed a higher richness of putative Bd-inhibitory bacterial strains in adults from disturbed habitats than in those from undisturbed habitats, as well as a greater number of these potential protective bacteria with a high relative abundance.

CONCLUSIONS

Our findings support the microbial "Anna Karenina principle", in which disturbance is hypothesized to cause greater microbial dispersion in communities, a so-called dysbiosis, which is a response of animal microbiomes to stress factors that decrease the ability of the host or its microbiome to regulate community composition. On the positive side, the high richness and relative abundance of putative Bd-inhibitory bacteria may indicate the development of a defense mechanism that enhances Bd-protection, attributed to a co-occurrence of more than 30-years of host and pathogen in these disturbed habitats. Our results provide important insight into the influence of human-modified landscapes on the skin microbiome and health implications of Bd-survivor species.

Environmental Factors and Host Microbiomes Shape Host-Pathogen Dynamics

Microorganisms are increasingly recognized as ecosystem-relevant components because they affect the population dynamics of hosts. Functioning at the interface of the host and pathogen, skin and gut microbiomes are vital components of immunity. Recent work reveals a strong influence of biotic and abiotic environmental factors (including the environmental microbiome) on disease dynamics, yet the importance of the host-host microbiome-pathogen-environment interaction has been poorly reflected in theory. We use amphibians and the disease chytridiomycosis caused by the fungal pathogen Batrachochytrium dendrobatidis to show how interactions between host, host microbiome, pathogen, and the environment all affect disease outcome. Our review provides new perspectives that improve our understanding of disease dynamics and ecology by incorporating environmental factors and microbiomes into disease theory.

Variation in the Slimy Salamander (Plethodon spp.) Skin and Gut-Microbial Assemblages Is Explained by Geographic Distance and Host Affinity

A multicellular host and its microbial communities are recognized as a metaorganism-a composite unit of evolution. Microbial communities have a variety of positive and negative effects on the host life history, ecology, and evolution. This study used high-throughput amplicon sequencing to characterize the complete skin and gut microbial communities, including both bacteria and fungi, of a terrestrial salamander, Plethodon glutinosus (Family Plethodontidae). We assessed salamander populations, representing nine mitochondrial haplotypes ('clades'), for differences in microbial assemblages across 13 geographic locations in the Southeastern United States. We hypothesized that microbial assemblages were structured by both host factors and geographic distance. We found a strong correlation between all microbial assemblages at close geographic distances, whereas, as spatial distance increases, the patterns became increasingly discriminate. Network analyses revealed that gut-bacterial communities have the highest degree of connectedness across geographic space. Host salamander clade was explanatory of skin-bacterial and gut-fungal assemblages but not gut-bacterial assemblages, unless the latter were analyzed within a phylogenetic context. We also inferred the function of gut-fungal assemblages to understand how an understudied component of the gut microbiome may influence salamander life history. We concluded that dispersal limitation may in part describe patterns in microbial assemblages across space and also that the salamander host may select for skin and gut communities that are maintained over time in closely related salamander populations.

Moving Beyond the Host: Unraveling the Skin Microbiome of Endangered Costa Rican Amphibians

Some neotropical amphibians, including a few species in Costa Rica, were presumed to be "extinct" after dramatic population declines in the late 1980s but have been rediscovered in isolated populations. Such populations seem to have evolved a resistance/tolerance to Batrachochytrium dendrobatidis (Bd), a fungal pathogen that causes a deadly skin disease and is considered one of the main drivers of worldwide amphibian declines. The skin microbiome is an important component of the host's innate immune system and is associated with Bd-resistance. However, the way that the bacterial diversity of the skin microbiome confers protection against Bd in surviving species remains unclear. We studied variation in the skin microbiome and the prevalence of putatively anti-Bd bacterial taxa in four co-habiting species in the highlands of the Juan Castro Blanco National Park in Costa Rica using 16S rRNA amplicon sequencing. Lithobates vibicarius, Craugastor escoces, and Isthmohyla rivularis have recently been rediscovered, whereas Isthmohyla pseudopuma has suffered population fluctuations but has never disappeared. To investigate the life stage at which the protective skin microbiome is shaped and when shifts occur in the diversity of putatively anti-Bd bacteria, we studied the skin microbiome of tadpoles, juveniles and adults of L. vibicarius. We show that the skin bacterial composition of sympatric species and hosts with distinct Bd-infection statuses differs at the phyla, family, and genus level. We detected 94 amplicon sequence variants (ASVs) with putative anti-Bd activity pertaining to distinct bacterial taxa, e.g., Pseudomonas spp., Acinetobacter johnsonii, and Stenotrophomonas maltophilia. Bd-uninfected L. vibicarius harbored 79% more putatively anti-Bd ASVs than Bd-infected individuals. Although microbiome composition and structure differed across life stages, the diversity of putative anti-Bd bacteria was similar between pre- and post-metamorphic stages of L. vibicarius. Despite low sample size, our results support the idea that the skin microbiome is dynamic and protects against ongoing Bd presence in endangered species persisting after their presumed extinction. Our study serves as a baseline to understand the microbial patterns in species of high conservation value. Identification of microbial signatures linked to variation in disease susceptibility might, therefore, inform mitigation strategies for combating the global decline of amphibians.