Temporal changes in cutaneous bacterial communities of terrestrial‐ and aquatic‐phase newts (Amphibia)

Rice straw ash and amphibian health: A deep dive into microbiota changes and potential ecological consequences

Rice straw is burned as a result of agricultural practices and technical limitations, generating significant volumes of ash that might have environmental and ecological consequences; however, the effects on organisms have not been researched. Amphibians depend on their gut and skin microbiomes. Ash exposure may cause inflammation and changes in microbial diversity and function in frogs' skin and gut microbiota due to its chemical composition and physical presence, but the implications remain unclear. Rana dybowskii were exposed to five aqueous extracts of ashes (AEA) concentrations for 30 days to study survival, metal concentrations, and microbial diversity, analyzing the microbiota of the cutaneous and gut microbiota using Illumina sequencing. Dominant elements in ash: K > Ca > Mg > Na > Al > Fe. In AEA, K > Na > Ca > Mg > As > Cu. Increased AEA concentrations significantly reduced frog survival. Skin microbiota alpha diversity varied significantly among all treatment groups, but not gut microbiota. Skin microbiota differed significantly across treatments via Bray-Curtis and weighted UniFrac; gut microbiota was only affected by Bray-Curtis. Skin microbiota varied significantly with AEA levels in Proteobacteria, Bacteroidetes, Actinobacteria, and Firmicutes, while the gut microbiota's dominant phyla, Firmicutes, Bacteroidetes, and Proteobacteria, remained consistent across all groups. Lastly, the functional prediction showed that the skin microbiota had big differences in how it worked and looked, which were linked to different health and environmental adaptation pathways. The gut microbiota, on the other hand, had smaller differences. In conclusion, AEA exposure affects R. dybowskii survival and skin microbiota diversity, indicating potential health and ecological impacts, with less effect on gut microbiota.

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.

Ontogeny drives shifts in skin bacterial communities in facultatively paedomorphic salamanders

Microbiomes are major determinants of host growth, development and survival. In amphibians, host-associated bacteria in the skin can inhibit pathogen infection, but many processes can influence the structure and composition of the community. Here we quantified the shifts in skin-associated bacteria across developmental stages in the striped newt (Notophthalmus perstriatus), a threatened salamander species with a complex life history and vulnerable to infection by the amphibian chytrid fungus Batrachochytrium dendrobatidis and ranavirus. Our analyses show that pre-metamorphic larval and paedomorphic stages share similar bacterial compositions, and that the changes in the microbiome coincided with physiological restructuring during metamorphosis. Newts undergoing metamorphosis exhibited microbiome compositions that were intermediate between paedomorphic and post-metamorphic stages, further supporting the idea that metamorphosis is a major driver of host-associated microbes in amphibians. We did not find support for infection-related disruption of the microbiome, though infection replicates were small for each respective life stage.

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

Life stage and proximity to roads shape the skin microbiota of eastern newts (Notophthalmus viridescens)

Host-associated microbiomes play an essential role in the health of organisms, including immune system activation, metabolism and energy uptake. It is well established that microbial communities differ depending on the life stage and natural history of the organism. However, the effects of life stage and natural history on microbial communities may also be influenced by human activities. We investigated the effects of amphibian life stage (terrestrial eft vs. aquatic adult) and proximity to roadways on newt skin bacterial communities. We found that the eft and adult life stages differed in bacterial community composition; however, the effects of roads on community composition were more evident in the terrestrial eft stage compared to the aquatic adult stage. Terrestrial efts sampled close to roads possessed richer communities than those living further away from the influence of roads. When accounting for amplicon sequence variants with predicted antifungal capabilities, in the adult life stage, we observed a decrease in anti-fungal bacteria with distance to roads. In contrast, in the eft stage, we found an increase in anti-fungal bacteria with distance to roads. Our results highlight the need to consider the effects of human activities when evaluating how host-associated microbiomes differ across life stages of wildlife.

Temporal variation in skin microbiota of cohabitating amphibians

Temporal changes and transmission patterns in host-associated microbial communities have important implications for host health. The diversity of amphibian skin microbial communities is associated with disease outcome in amphibians exposed to the fungal pathogen Batrachochytrium dendrobatidis (Bd). To successfully develop conservation strategies against Bd, we need a comprehensive understanding of how skin microbes are maintained and transmitted over time within populations. We used 16S rRNA sequence analysis to compare Epipedobates anthonyi frogs housed with one conspecific to frogs housed singly at four time points over the course of 1 year. We found that both α and β diversity of frog skin bacterial communities changed significantly over the course of the experiment. Specifically, we found that bacterial communities of cohabitating frogs became more similar over time. We also observed that some bacterial taxa were differentially abundant between frogs housed singly and frogs housed with a conspecific. These results suggest that conspecific contact may play a role in mediating amphibian skin microbial diversity and that turnover of skin microbial communities can occur across time. Our findings provide rationale for future studies exploring horizontal transmission as a potential mechanism of host-associated microbial maintenance in amphibians.

Seasonal changes and the unexpected impact of environmental disturbance on skin bacteria of individual amphibians in a natural habitat

Amphibians host diverse skin bacteria that have a role in pathogen defense, but these skin communities could change over time and impact this function. Here, we monitored individual Eastern red-spotted newts (Notophthalmus viridescens; N = 17) for 2 years in a field pond enclosure and assessed the effects of season and disturbance on skin bacterial community dynamics. We created disturbances by adding additional pond substrate to the enclosure at two timepoints. We planned to sample the skin bacterial community and metabolite profiles of each newt every 6 weeks; we ultimately sampled eight individuals at least six times. We used 16S rRNA gene amplicon sequencing to characterize the bacterial communities and HPLC-MS for metabolite profiling. We found that disturbance had a dramatic effect on skin bacterial communities and metabolite profiles, while season had an effect only using select metrics. There were seven core bacterial taxa (97% OTUs) that were found on all newts in all seasons, pre- and post-disturbance. Lastly, there was a correlation between bacterial and metabolite profiles post-disturbance, which was not observed pre-disturbance. This longitudinal study suggests that environmental disturbances can have lasting effects on skin bacterial communities that overwhelm seasonal changes, although the core bacteria remain relatively consistent over time.

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.

Comparison of the gut microbiota of Rana amurensis and Rana dybowskii under natural winter fasting conditions

Rana amurensis and R. dybowskii occupy similar habitats. As temperatures decrease with the onset of winter, both species migrate to ponds for hibernation. Our goal was to determine whether different species possess different intestinal microbiota under natural winter fasting conditions. We used high-throughput Illumina sequencing of 16S rRNA gene sequences to analyse the diversity of intestinal microbes in the two species. The dominant gut bacterial phyla in both species were Bacteroidetes, Proteobacteria and Firmicutes. Linear discriminant analysis (LDA) effect size revealed significant enrichment of Proteobacteria in R. amurensis and Firmicutes in R. dybowskii. There were significant differences in the gut microbiota composition between the species. The core operational taxonomic unit numbers in R. amurensis and R. dybowskii shared by the two species were 106, 100 and 36. This study indicates that the intestinal bacterial communities of the two frog species are clearly different. Phylum-level analysis showed that R. amurensis was more abundant in Proteobacteria and Verrucomicrobia than R. dybowskii was This is the first study of the composition and diversity of the gut microbiota of these two species, providing important insights for future research on the gut microbiota and the role of these bacterial communities in frogs.

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.

Effects of Seasonal Hibernation on the Similarities Between the Skin Microbiota and Gut Microbiota of an Amphibian (Rana dybowskii)

Both the gut and skin microbiotas have important functions for amphibians. The gut microbiota plays an important role in both the health and evolution of the host species, whereas the role of skin microbiota in disease resistance is particularly important for amphibians. Many studies have examined the effects of environmental factors on the skin and gut microbiotas, but no study has yet explored the similarities between the skin and gut microbiotas. In this study, the gut and skin microbiotas of Rana dybowskii in summer and winter were investigated via high-throughput Illumina sequencing. The results showed that the alpha diversity of gut and skin microbiotas decreased significantly from summer to winter. In both seasons, the microbial composition and structure differed significantly between the gut and skin, and the similarities between these microbiotas differed between seasons. The pairwise distances between the gut and skin microbiotas were greater in winter than in summer. The ratio of core OTUs and shared OTUs to the sum of the OTUs in the gut and skin microbiotas in summer was significantly higher than that in winter. The similarities between the gut and skin microbiotas are important for understanding amphibian ecology and life history.

The Skin Microbiota of Eleutherodactylus Frogs: Effects of Host Ecology, Phylogeny, and Local Environment

Amphibian skin microbiota has a potential protective role against diseases. However, the effects of environmental and host factors on symbiotic bacterial communities are not well understood. Caribbean frogs in the genus Eleutherodactylus represent a case of congeneric species that differ in ecological specialization by the process of adaptive radiation. For a small clade of Eleutherodactylus from Puerto Rico, we investigated the role of local environments, host species, and microhabitat in the composition of their skin microbiome. The potential congruence between microbial communities in hosts that are most closely related phylogenetically was also addressed. We hypothesized that the skin microbiota of Eleutherodactylus frogs would be mostly associated to microhabitat use, but also differ according to locality, and to a lesser extent to host species. To test this hypothesis, we swabbed the skin of a total of 98 adult individuals of seven Eleutherodactylus species distributed in two nearby localities in Puerto Rico, and sequenced the V4 region of the 16S rRNA gene. Results showed that locality had the greatest effect on determining skin bacterial communities of amphibian hosts, but this effect was stronger on the composition (based on presence/absence) than on its structure (based on sequence abundance). The most ecologically distinct host, E. cooki, and the generalist E. coqui presented, respectively, the most dissimilar and similar microbiota compared to other hosts. Host phylogeny showed a weak influence on skin microbiota. Results suggest that both local environment and ecological specialization are structuring the skin bacterial community in these Eleutherodactylus species, but that characteristics intrinsic to species may also render unique hosts the ability to maintain distinct microbiotas.

Skin bacterial communities of neotropical treefrogs vary with local environmental conditions at the time of sampling

The amphibian skin microbiome has been the focus of recent studies aiming to better understand the role of these microbial symbionts in host defense against disease. However, host-associated microbial communities are complex and dynamic, and changes in their composition and structure can influence their function. Understanding temporal variation of bacterial communities on amphibian skin is critical for establishing baselines from which to improve the development of mitigation techniques based on probiotic therapy and provides long-term host protection in a changing environment. Here, we investigated whether microbial communities on amphibian skin change over time at a single site. To examine this, we collected skin swabs from two pond-breeding species of treefrogs, Agalychnis callidryas and Dendropsophus ebraccatus, over 4 years at a single lowland tropical pond in Panamá. Relative abundance of operational taxonomic units (OTUs) based on 16S rRNA gene amplicon sequencing was used to determine bacterial community diversity on the skin of both treefrog species. We found significant variation in bacterial community structure across long and short-term time scales. Skin bacterial communities differed across years on both species and between seasons and sampling days only in D. ebraccatus. Importantly, bacterial community structures across days were as variable as year level comparisons. The differences in bacterial community were driven primarily by differences in relative abundance of key OTUs and explained by rainfall at the time of sampling. These findings suggest that skin-associated microbiomes are highly variable across time, and that for tropical lowland sites, rainfall is a good predictor of variability. However, more research is necessary to elucidate the significance of temporal variation in bacterial skin communities and their maintenance for amphibian conservation efforts.

Characterization of the Skin Microbiota of the Cane Toad Rhinella cf. marina in Puerto Rico and Costa Rica

Rhinella marina is a toad native to South America that has been introduced in the Antilles, likely carrying high loads of microorganisms, potentially impacting local community diversity. The amphibian skin is involved in pathogen defense and its microbiota has been relatively well studied, however, research focusing on the cane toad microbiota is lacking. We hypothesize that the skin microbial communities will differ between toads inhabiting different geographical regions in Central America and the Caribbean. To test our hypothesis, we compared the microbiota of three populations of R. cf. marina toads, two from Costa Rican (native) and one Puerto Rican (exotic) locations. In Costa Rica, we collected 11 toads, 7 in Sarapiquí and 4 from Turrialba while in Puerto Rico, 10 animals were collected in Santa Ana. Separate swab samples were collected from the dorsal and ventral sites resulting in 42 samples. We found significant differences in the structure of the microbial communities between Puerto Rico and Costa Rica. We detected as much as 35 different phyla; however, communities were dominated by Proteobacteria, Bacteroidetes, Firmicutes, and Actinobacteria. Alpha diversity and richness were significantly higher in toads from Puerto Rico and betadiversity revealed significant differences between the microbiota samples from the two countries. At the genus level, we found in Santa Ana, Puerto Rico, a high dominance of Kokuria, Niabella, and Rhodobacteraceae, while in Costa Rica we found Halomonas and Pseudomonas in Sarapiquí, and Acinetobacter and Citrobacter in Turrialba. This is the first report of Niabella associated with the amphibian skin. The core microbiome represented 128 Operational Taxonomic Units (OTUs) mainly from five genera shared among all samples, which may represent the symbiotic Rhinella's skin. These results provide insights into the habitat-induced microbial changes facing this amphibian species. The differences in the microbial diversity in Puerto Rican toads compared to those in Costa Rica provide additional evidence of the geographically induced patterns in the amphibian skin microbiome, and highlight the importance of discussing the microbial tradeoffs in the colonization of new ecosystems.

Temporal Variation of the Skin Bacterial Community and Batrachochytrium dendrobatidis Infection in the Terrestrial Cryptic Frog Philoria loveridgei

In animals and plants, symbiotic bacteria can play an important role in disease resistance of host and are the focus of much current research. Globally, amphibian population declines and extinctions have occurred due to chytridiomycosis, a skin disease caused by the pathogen Batrachochytrium dendrobatidis (Bd). Currently amphibian skin bacteria are increasingly recognized as important symbiont communities with a relevant role in the defense against pathogens, as some bacteria can inhibit the growth of B. dendrobatidis. This study aims to document the B. dendrobatidis infection status of wild populations of a terrestrial cryptic frog (Philoria loveridgei), and to determine whether infection status is correlated with changes in the skin microbial communities. Skin samples of P. loveridgei were collected along an altitudinal range within the species distribution in subtropical rainforests in southeast Australia. Sampling was conducted in two years during two breeding seasons with the first classified as a “La Niña” year. We used Taqman real-time PCR to determine B. dendrobatidis infection status and 16S amplicon sequencing techniques to describe the skin community structure. We found B. dendrobatidis-positive frogs only in the second sampling year with low infection intensities, and no correlation between B. dendrobatidis infection status and altitude, frog sex or size. Skin bacterial diversity was significantly higher in P. loveridgei frogs sampled in the 1st year than in the 2nd year. In addition, 7.4% of the total OTUs were significantly more abundant in the 1st year compared to the 2nd year. We identified 67 bacterial OTUs with a significant positive correlation between infection intensity and an OTU’s relative abundance. Forty-five percent of these OTUs belonged to the family Enterobacteriaceae. Overall, temporal variation was strongly associated with changes in B. dendrobatidis infection status and bacterial community structure of wild populations of P. loveridgei.