Host and Aquatic Environment Shape the Amphibian Skin Microbiome but Effects on Downstream Resistance to the Pathogen Batrachochytrium dendrobatidis Are Variable

The commensal skin microbiome of amphibian mountain populations and its association with the pathogen Batrachochytrium dendrobatidis

Microbial assemblages naturally living on the skin are an integral part of immunity. In amphibians, this skin microbiota may hold a mitigation solution against the fungal pathogen Batrachochytrium dendrobatidis (Bd), which causes the panzootic disease chytridiomycosis. We used 16S rRNA gene metabarcoding to test the adaptive microbiome hypothesis. We compared the community composition, richness, and putative Bd-inhibitory function of the skin microbiome of three amphibian host species in the Pyrenees, as well as three species in Taiwan, in both Bd-positive and negative mountain populations. In both geographical regions, the amphibian host species played a decisive role in shaping the microbial assemblage and putative anti-Bd properties. In the Pyrenees, the species most susceptible to chytridiomycosis, Alytes obstetricans, had the lowest relative abundances of putative protective bacteria. In Bd-positive and negative sites, individuals had different skin microbiomes, with all anuran species showing increased relative abundances of potential anti-Bd bacteria, while the Taiwanese caudata Hynobius sonani showed the opposite pattern. Our results suggest that, in response to exposure to the pathogen, the skin microbiota shifted to a defensive state with increased anti-Bd function, which may contribute to promoting disease resistance, as proposed by the adaptive microbiome hypothesis.

From Skin to Gut: Understanding Microbial Diversity in Rana amurensis and R. dybowskii

Amphibians face the threat of decline and extinction, and their health is crucially affected by the microbiota. Their health and ecological adaptability essentially depend on the diverse microbial communities that are shaped by unique host traits and environmental factors. However, there is still limited research on this topic. In this study, cutaneous (C) and gut (G) microbiota in Rana amurensis (A) and R. dybowskii (D) was analyzed through 16S amplicon sequencing. Groups AC and DC significantly differed in alpha diversity, while the gut groups (AG and DG) showed no such differences. Analyses of Bray-Curtis dissimilarity matrix and unweighted UniFrac distances showed significant differences in cutaneous microbiota between groups AC and DC, but not between groups AG and DG. Stochastic processes significantly influenced the assembly of cutaneous and gut microbiota in amphibians, with a notably higher species dispersal rate in the gut. The predominant phyla in the skin of R. amurensis and R. dybowskii were Bacteroidetes and Proteobacteria, respectively, with significant variations in Bacteroidota. Contrarily, the gut microbiota of both species was dominated by Firmicutes, Proteobacteria, and Bacteroidetes, without significant phylum-level differences. Linear discriminant analysis effect size (LEfSe) analysis identified distinct microbial enrichment in each group. Predictive analysis using phylogenetic investigation of communities by reconstruction of unobserved states 2 (PICRUSt2) revealed the significant functional pathways associated with the microbiota, which indicates their potential roles in immune system function, development, regeneration, and response to infectious diseases. This research underscores the critical impact of both host and environmental factors in shaping amphibian microbial ecosystems and emphasizes the need for further studies to explore these complex interactions for conservation efforts.

Host-associated helminth diversity and microbiome composition contribute to anti-pathogen defences in tropical frogs impacted by forest fragmentation

Habitat fragmentation can negatively impact wildlife populations by simplification of ecological interactions, but little is known about how these impacts extend to host-associated symbiotic communities. The symbiotic communities of amphibians play important roles in anti-pathogen defences, particularly against the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). In this study, we analyse the role of macroparasitic helminth communities in concert with microbial communities in defending the host against Bd infection within the context of forest fragmentation. We found that skin microbial and helminth communities are disrupted at fragmented habitats, while gut microbiomes appear more resilient to environmental change. We also detected potential protective roles of helminth diversity and anti-pathogen microbial function in limiting Bd infection. Microbial network analysis revealed strong patterns of structure in both skin and gut communities, with helminths playing central roles in these networks. We reveal consistent roles of microbial and helminth diversity in driving host-pathogen interactions and the potential implications of fragmentation on host fitness.

Nor climate nor human impact factors: Chytrid infection shapes the skin bacterial communities of an endemic amphibian in a biodiversity hotspot

The bacterial communities of the amphibian skin (i.e., the bacteriome) are critical to the host's innate immune system. However, it is unclear how different drivers can alter this function by modulating the bacteriome's structure. Our aim was to assess the extent to which different host attributes and extrinsic factors influence the structure of the bacterial communities of the skin. Skin bacterial diversity was examined in 148 individuals of the four-eyed frog (Pleurodema thaul) from 16 localities spanning almost 1800 km in latitude. The richness and beta diversity of bacterial families and the richness and abundance of Bd-inhibitory bacterial genera were used to describe their structure. Predictors associated with the host (developmental stage, genetic lineage, individual Batrachochytrium dendrobatidis [Bd] infection status) and the landscape (current climate, degree of anthropogenic disturbance) were used in the statistical modeling in an information theoretical approach. Bd infection and host developmental stage were the only predictors affecting bacteriome richness, with Bd+ individuals and postmetamorphic stages (adults and juveniles) having higher richness than Bd- ones and tadpoles. High diversity in Bd+ individuals is not driven by bacterial genera with known anti-Bd properties. Beta diversity was not affected by Bd infection and was mainly a consequence of bacterial family turnover rather than nestedness. Finally, for those bacterial genera known to have inhibitory effects on chytrid, Bd+ individuals had a slightly higher diversity than Bd- ones. Our study confirms an association between Bd infection and the host developmental stage with the bacterial communities of the skin of P. thaul. Unexpectedly, macroclimate and human impact factors do not seem to play a role in shaping the amphibian skin microbiome. Our study exemplifies that focusing on a single host-parasite system over a large geographic scale can provide essential insights into the factors driving host-parasite-bacteriome interactions.

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.

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.

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.

Wildfires disturb the natural skin microbiota of terrestrial salamanders

Environmental change can disturb natural associations between wildlife and microbial symbionts, in many cases to the detriment of host health. We used a North American terrestrial salamander system to assess how the skin microbiota of amphibians responds to wildfires. In northern California's redwood/oak forests, we assessed how recent wildfires affected the skin microbiota of three different salamander species (Taricha sp., Batrachoseps attenuatus, and Ensatina eschscholtzii) over two different sampling seasons in 2018 and 2021. We found species-specific responses to wildfire disturbance on the alpha diversity of the skin microbiota of terrestrial salamanders, although burning in general altered the composition of the skin microbiota. The effect of burning on alpha diversities and body condition indices varied by sampling season, suggesting an additional effect of annual climatic conditions on body condition and skin microbiota response. We tested all salamanders for Batrachochytrium dendrobatidis and found four infected individuals in 2018 and none in 2021. Our study documents correlations in the skin microbiota response to an increasing source of disturbance in western North American ecosystems. In addition, our results highlight the need to consider the effects of increased wildfire regimes/intensities and longitudinal effects on wildlife-associated microbiota and animal health.

Skin and gut microbiomes of tadpoles vary differently with host and water environment: a short-term experiment using 16S metabarcoding

The host-microbiome community is influenced by several host and environmental factors. In order to disentangle the individual effects of host and environment, we performed a laboratory experiment to assess the effects of the exposure to different water sources on the skin and gut microbiome of two amphibian species (Pelophylax perezi and Bufo spinosus). We observed that the bacterial communities greatly varied with water environment and host identity. Tadpoles of B. spinosus collected from a waterbody with poorer bacterial diversity exhibited a more diverse skin and gut microbiome after exposed to a richer water source. Tadpoles of P. perezi, originally collected from a richer water environment, exhibited less marked alterations in diversity patterns independently of the water source but showed alterations in gut composition. These results highlight that environment alterations, such as the water source, combined with the host effect, impact the microbiome of amphibian species in different ways; the population history (e.g., previous water environment and habitat) of the host species may also influence future alterations on tadpole microbiome.

Host genotype and microbiome associations in co-occurring clonal and non-clonal kelp, Ecklonia radiata

A fundamental question in holobiont biology is the extent to which microbiomes are determined by host characteristics regulated by their genotype. Studies on the interactions of host genotype and microbiomes are emerging but disentangling the role that host genotype has in shaping microbiomes remains challenging in natural settings. Host genotypes tend to be segregated in space and affected by different environments. Here we overcome this challenge by studying an unusual situation where host asexual (5 clonal lineages) and sexual genotypes (15 non-clonal lineages) of the same species co-occur under the same environment. This allowed us to partition the influence of morphological traits and genotype in shaping host-associated bacterial communities. Lamina-associated bacteria of co-occurring kelp sexual non-clonal (Ecklonia radiata) and asexual clonal (E. brevipes) morphs were compared to test whether host genotype influences microbiomes beyond morphology. Similarity of bacterial composition and predicted functions were evaluated among individuals within a single clonal genotype or among non-clonal genotypes of each morph. Higher similarity in bacterial composition and inferred functions were found among identical clones of E. brevipes compared to other clonal genotypes or unique non-clonal E. radiata genotypes. Additionally, bacterial diversity and composition differed significantly between the two morphs and were related with one morphological trait in E. brevipes (haptera). Thus, factors regulated by the host genotype (e.g. secondary metabolite production) likely drive differences in microbial communities between morphs. The strong association of genotype and microbiome found here highlights the importance of genetic relatedness of hosts in determining variability in their bacterial symbionts.

Skin bacterial community differences among three species of co-occurring Ranid frogs

Skin microbial communities are an essential part of host health and can play a role in mitigating disease. Host and environmental factors can shape and alter these microbial communities and, therefore, we need to understand to what extent these factors influence microbial communities and how this can impact disease dynamics. Microbial communities have been studied in amphibian systems due to skin microbial communities providing some resistance to the amphibian chytrid fungus, Batrachochytrium dendrobatidis. However, we are only starting to understand how host and environmental factors shape these communities for amphibians. In this study, we examined whether amphibian skin bacterial communities differ among host species, host infection status, host developmental stage, and host habitat. We collected skin swabs from tadpoles and adults of three Ranid frog species (Lithobates spp.) at the Mianus River Gorge Preserve in Bedford, New York, USA, and used 16S rRNA gene amplicon sequencing to determine bacterial community composition. Our analysis suggests amphibian skin bacterial communities change across host developmental stages, as has been documented previously. Additionally, we found that skin bacterial communities differed among Ranid species, with skin communities on the host species captured in streams or bogs differing from the communities of the species captured on land. Thus, habitat use of different species may drive differences in host-associated microbial communities for closely-related host species.

Colonization order of bacterial isolates on treefrog embryos impacts microbiome structure in tadpoles

Priority effects, or impacts of colonization order, may have lasting influence on ecological community composition. The embryonic microbiome is subject to stochasticity in colonization order of bacteria. Stochasticity may be especially impactful for embryos developing in bacteria-rich environments, such as the embryos of many amphibians. To determine if priority effects experienced as embryos impacted bacterial community composition in newly hatched tadpoles, we selectively inoculated the embryos of laboratory-raised hourglass treefrogs, Dendropsophus ebraccatus, with bacteria initially isolated from the skin of wild D. ebraccatus adults over 2 days. First, embryos were inoculated with two bacteria in alternating sequences. Next, we evaluated the outcomes of priority effects in an in vitro co-culture assay absent of host factors. We then performed a second embryo experiment, inoculating embryos with one of three bacteria on the first day and a community of five target bacteria on the second. Through 16S rRNA gene amplicon sequencing, we observed relative abundance shifts in tadpole bacteria communities due to priority effects. Our results suggest that the initial bacterial source pools of embryos shape bacterial communities at later life stages; however, the magnitude of those changes is dependent on the host environment and the identity of bacterial colonists.

Geographic variation in bacterial assemblages on cane toad skin is influenced more by local environments than by evolved changes in host traits

Bacterial assemblages on amphibian skin may play an important role in protecting hosts against infection. In hosts that occur over a range of environments, geographic variation in composition of bacterial assemblages might be due to direct effects of local factors and/or to evolved characteristics of the host. Invasive cane toads (Rhinella marina) are an ideal candidate to evaluate environmental and genetic mechanisms, because toads have evolved major shifts in physiology, morphology, and behavior during their brief history in Australia. We used samples from free-ranging toads to quantify site-level differences in bacterial assemblages and a common-garden experiment to see if those differences disappeared when toads were raised under standardised conditions at one site. The large differences in bacterial communities on toads from different regions were not seen in offspring raised in a common environment. Relaxing bacterial clustering to operational taxonomic units in place of amplicon sequence variants likewise revealed high similarity among bacterial assemblages on toads in the common-garden study, and with free-ranging toads captured nearby. Thus, the marked geographic divergence in bacterial assemblages on wild-caught cane toads across their Australian invasion appears to result primarily from local environmental effects rather than evolved shifts in the host.

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.

Neglected skin-associated microbial communities: a unique immune defense strategy of Bufo raddei under environmental heavy metal pollution

Amphibians defend against pathogens using skin microbial communities, in addition to innate and adaptive immunity. Despite skin microbial communities play a key role in the immune function of amphibians, few studies have focused on the changes in its composition and function. In the present study, we identified the variation in adaptive immunity, as well as the corresponding changes in skin microbiome of Bufo raddei living in a heavy metal polluted area. The adaptive immunity of B. raddei in heavy metal polluted area was significantly lower than that in relatively unpolluted area. Further, different skin bacterial communities were found in the two areas. In the heavy metal polluted area, Actinobacteria and Microbacterium were the dominant bacteria in the skin microbiome of B. raddei, which showed broad-spectrum antibacterial activity. Besides, the antibiotic synthesis was also increased in metabolic pathways. The present study suggested that the adaptive immunity of B. raddei was weakened under long-term heavy metal stress. However, the toads increased the abundance of bacteriostatic bacteria by regulating the composition of skin microbiome, which released a large number of bacteriostatic metabolites and enhanced the host resistance to external pathogens in turn.

Phylogenetically under-dispersed gut microbiomes are not correlated with host genomic heterozygosity in a genetically diverse reptile community

While key elements of fitness in vertebrate animals are impacted by their microbiomes, the host genetic characteristics that factor into microbiome composition are not fully understood. Here, we correlate host genomic heterozygosity and gut microbiome phylogenetic diversity across a community of reptiles in southwestern New Mexico to test hypotheses about the behaviour of host genes that drive microbiome assembly. We find that microbiome communities are phylogenetically under-dispersed relative to random expectations, and that host heterozygosity is not correlated with microbiome diversity. Our analyses reinforce results from functional genomic work that identify conserved host immune and nonimmune genes as key players in microbiome assembly, rather than gene families that rely on heterozygosity for their function.

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.

Environmental and Anthropogenic Factors Shape the Skin Bacterial Communities of a Semi-Arid Amphibian Species

The amphibian skin microbiome is important in maintaining host health, but is vulnerable to perturbation from changes in biotic and abiotic conditions. Anthropogenic habitat disturbance and emerging infectious diseases are both potential disrupters of the skin microbiome, in addition to being major drivers of amphibian decline globally. We investigated how host environment (hydrology, habitat disturbance), pathogen presence, and host biology (life stage) impact the skin microbiome of wild Dhofar toads (Duttaphrynus dhufarensis) in Oman. We detected ranavirus (but not Batrachochytrium dendrobatidis) across all sampling sites, constituting the first report of this pathogen in Oman, with reduced prevalence in disturbed sites. We show that skin microbiome beta diversity is driven by host life stage, water source, and habitat disturbance, but not ranavirus infection. Finally, although trends in bacterial diversity and differential abundance were evident in disturbed versus undisturbed sites, bacterial co-occurrence patterns determined through network analyses revealed high site specificity. Our results therefore provide support for amphibian skin microbiome diversity and taxa abundance being associated with habitat disturbance, with bacterial co-occurrence (and likely broader aspects of microbial community ecology) being largely site specific.

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.

Short-term continuous and pulse Pb exposure causes negative effects on skin histomorphological structure and bacterial composition of adult Pelophylax nigromaculatus

Lead (Pb) is a toxic heavy metal often present in the environment as a pulse in water. Traditional toxicity tests are usually carried out under conditions of continuous concentration, without considering the impact of pulse exposure on aquatic organisms. This study aimed to evaluate the effects of short-term continuous and pulse Pb exposures on the skin bacteria and histomorphological structure of Pelophylax nigromaculatus. Results showed that compared to the control (CON) and Pb continuous exposure group (CEPb), the Pb pulse exposure group (PEPb) showed the smallest size of granular glands, which would interfere with the permeability and secretory function of skin, making the individual more sensitive to external pollution. Lead exposure significantly changed the composition and diversity of skin bacteria. Compared to the CON and CEPb groups, the PEPb group showed a significant increase in the abundance of harmful bacteria (e.g., Bacteroidetes and Chryseobacterium) and a decrease in the abundance of beneficial bacteria (e.g., Pseudomonas). PICRUSt software showed that there were differences in the metabolic pathway of skin bacteria among the three groups (CON, CEPb, and PEPb). Overall, this study indicates that Pb pulse exposure can aggravate the toxicity of Pb for frog skin, providing a new framework for simulating short-term heavy metal exposure in the context of frog health.

First Report of Culturable Skin Bacteria in Melanophryniscus admirabilis (Admirable Redbelly Toad)

Melanophryniscus admirabilis is a small toad, critically endangered with a microendemic distribution in the Atlantic Forest in southern Brazil. The amphibian skin microbiome is considered one of the first lines of defense against pathogenic infections, such as Batrachochytrium dendrobatidis (Bd). The knowledge of skin amphibian microbiomes is important to numerous fields, including species conservation, detection, and quantification of environmental changes and stressors. In the present study, we investigated, for the first time, cultivable bacteria in the skin of wild M. admirabilis, and detected Bd fungus by nested polymerase chain reaction (PCR) technique. Skin swab samples were collected from 15 wild M. admirabilis, and the isolation of bacteria was performed by means of different culture strategies. A total of 62 bacterial isolates being Bacillus (n = 22; 34.48%), Citrobacter (n = 10; 16.13%), and Serratia (n = 12; 19.35%) were more frequently isolated genera. Interestingly, all skin samples tested were Bd negative. Some bacterial genera identified in our study might be acting in a synergic relationship and protecting them against the Bd fungus. In addition, these bacteria may play an essential role in maintaining this species in an environment modulated by anthropic actions. This first report of skin cultivable bacteria from M. admirabilis natural population improves our knowledge of skin amphibian microbiomes, contributing to a better understanding of their ecology and how this species has survived in an environment modulated by anthropic action.

Integrated Analysis of lncRNA and circRNA Mediated ceRNA Regulatory Networks in Skin Reveals Innate Immunity Differences Between Wild-Type and Yellow Mutant Rainbow Trout (Oncorhynchus mykiss)

Fish skin is a vital immune organ that forms the first protective barrier preventing entry of external pathogens. Rainbow trout is an important aquaculture fish species that is farmed worldwide. However, our knowledge of innate immunity differences between wild-type (WR_S) and yellow mutant rainbow trout (YR_S) remains limited. In this study, we performed whole transcriptome analysis of skin from WR_S and YR_S cultured in a natural flowing water pond. A total of 2448 mRNAs, 1630 lncRNAs, 22 circRNAs and 50 miRNAs were found to be differentially expressed (DE). Among these DEmRNAs, numerous key immune-related genes, including ifih1, dhx58, trim25, atp6v1e1, tap1, tap2, cd209, hsp90a.1, nlrp3, nlrc3, and several other genes associated with metabolism (gstp1, nampt, naprt and cd38) were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEmRNAs revealed that many were significantly enriched in innate immune-related GO terms and pathways, including NAD+ADP-ribosyltransferase activity, complement binding, immune response and response to bacterium GO terms, and RIG-I-like receptor signaling, NOD-like receptor signaling and phagosome KEGG pathways. Furthermore, the immune-related competing endogenous RNA networks were constructed, from which we found that lncRNAs MSTRG.11484.2, MSTRG.32014.1 and MSTRG.29012.1 regulated at least three immune-related genes (ifih1, dhx58 and irf3) through PC-5p-43254_34, PC-3p-28352_70 and bta-miR-11987_L-1R-1_1ss8TA, and tap2 was regulated by two circRNAs (circRNA5279 and circRNA5277) by oni-mir-124a-2-p5_1ss13GA. The findings expand our understanding of the innate immune system of rainbow trout, and lay the foundation for further study of immune mechanisms and disease resistance breeding.

Skin Microbiota Variation Among Bat Species in China and Their Potential Defense Against Pathogens

Host-associated skin bacteria are essential for resisting pathogen infections and maintaining health. However, we have little understanding of how chiropteran skin microbiota are distributed among bat species and their habitats, or of their putative roles in defending against Pseudogymnoascus destructans in China. In this study, we characterized the skin microbiomes of four bat species at five localities using 16S rRNA gene amplicon sequencing to understand their skin microbial composition, structure, and putative relationship with disease. The alpha- and beta-diversities of skin microbiota differed significantly among the bat species, and the differences were affected by environmental temperature, sampling sites, and host body condition. The chiropteran skin microbial communities were enriched in bacterial taxa that had low relative abundances in the environment. Most of the potential functions of skin microbiota in bat species were associated with metabolism. Focusing on their functions of defense against pathogens, we found that skin microbiota could metabolize a variety of active substances that could be potentially used to fight P. destructans. The skin microbial communities of bats in China are related to the environment and the bat host, and may be involved in the host's defense against pathogens.

Molecular basis of a bacterial-amphibian symbiosis revealed by comparative genomics, modeling, and functional testing

The molecular bases for the symbiosis of the amphibian skin microbiome with its host are poorly understood. Here, we used the odor-producer Pseudomonas sp. MPFS and the treefrog Boana prasina as a model to explore bacterial genome determinants and the resulting mechanisms facilitating symbiosis. Pseudomonas sp. MPFS and its closest relatives, within a new clade of the P. fluoresens Group, have large genomes and were isolated from fishes and plants, suggesting environmental plasticity. We annotated 16 biosynthetic gene clusters from the complete genome sequence of this strain, including those encoding the synthesis of compounds with known antifungal activity and of odorous methoxypyrazines that likely mediate sexual interactions in Boana prasina. Comparative genomics of Pseudomonas also revealed that Pseudomonas sp. MPFS and its closest relatives have acquired specific resistance mechanisms against host antimicrobial peptides (AMPs), specifically two extra copies of a multidrug efflux pump and the same two-component regulatory systems known to trigger adaptive resistance to AMPs in P. aeruginosa. Subsequent molecular modeling indicated that these regulatory systems interact with an AMP identified in Boana prasina through the highly acidic surfaces of the proteins comprising their sensory domains. In agreement with a symbiotic relationship and a highly selective antibacterial function, this AMP did not inhibit the growth of Pseudomonas sp. MPFS but inhibited the growth of another Pseudomonas species and Escherichia coli in laboratory tests. This study provides deeper insights into the molecular interaction of the bacteria-amphibian symbiosis and highlights the role of specific adaptive resistance toward AMPs of the hosts.

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.

Variability in environmental persistence but not per capita transmission rates of the amphibian chytrid fungus leads to differences in host infection prevalence

Heterogeneities in infections among host populations may arise through differences in environmental conditions through two mechanisms. First, environmental conditions may alter host exposure to pathogens via effects on survival. Second, environmental conditions may alter host susceptibility, making infection more or less likely if contact between a host and pathogen occurs. Further, host susceptibility might be altered through acquired resistance, which hosts can develop, in some systems, through exposure to dead or decaying pathogens and their metabolites. Environmental conditions may alter the rates of pathogen decomposition, influencing the likelihood of hosts developing acquired resistance. The present study primarily tests how environmental context influences the relative contributions of pathogen survival and per capita transmission on host infection prevalence using the amphibian chytrid fungus (Batrachochytrium dendrobatidis; Bd) as a model system. Secondarily, we evaluate how environmental context influences the decomposition of Bd because previous studies have shown that dead Bd and its metabolites can illicit acquired resistance in hosts. We conducted Bd survival and infection experiments and then fit models to discern how Bd mortality, decomposition and per capita transmission rates vary among water sources [e.g. artificial spring water (ASW) or water from three ponds]. We found that infection prevalence differed among water sources, which was driven by differences in mortality rates of Bd, rather than differences in per capita transmission rates. Bd mortality rates varied among pond water treatments and were lower in ASW compared to pond water. These results suggest that variation in Bd infection dynamics could be a function of environmental factors in waterbodies that result in differences in exposure of hosts to live Bd. In contrast to the persistence of live Bd, we found that the rates of decomposition of dead Bd did not vary among water sources, which may suggest that exposure of hosts to dead Bd or its metabolites might not commonly vary among nearby sites. Ultimately, a mechanistic understanding of the environmental dependence of free-living pathogens could lead to a deeper understanding of the patterns of outbreak heterogeneity, which could inform surveillance and management strategies.

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.

Characterization of the microbiome of the invasive Asian toad in Madagascar across the expansion range and comparison with a native co-occurring species

Biological invasions are on the rise, with each invader carrying a plethora of associated microbes. These microbes play important, yet poorly understood, ecological roles that can include assisting the hosts in colonization and adaptation processes or as possible pathogens. Understanding how these communities differ in an invasion scenario may help to understand the host's resilience and adaptability. The Asian common toad, Duttaphrynus melanostictus is an invasive amphibian, which has recently established in Madagascar and is expected to pose numerous threats to the native ecosystems. We characterized the skin and gut bacterial communities of D. melanostictus in Toamasina (Eastern Madagascar), and compared them to those of a co-occurring native frog species, Ptychadena mascareniensis, at three sites where the toad arrived in different years. Microbial composition did not vary among sites, showing that D. melanostictus keeps a stable community across its expansion but significant differences were observed between these two amphibians. Moreover, D. melanostictus had richer and more diverse communities and also harboured a high percentage of total unique taxa (skin: 80%; gut: 52%). These differences may reflect the combination of multiple host-associated factors including microhabitat selection, skin features and dietary preferences.

Batrachochytrium fungi: stealth invaders in amphibian skin

Amphibian populations around the world have been affected by two pathogenic fungi within the phylum Chytridiomycota. Batrachochytrium dendrobatidis (Bd) has infected hundreds of species and led to widespread declines and some species extinctions. Batrachochytrium salamandrivorans (Bsal) has devastated some native European salamanders, especially the iconic fire salamanders (Salamandra salamandra). Comparative genomic studies show that Bd is more diverse and widespread than previously thought, and global lineages occur together allowing for the development of hybrid lineages. New studies raise the concern of greater pathogenesis if both Bd and Bsal infect the same host. Although amphibians possess robust immune defenses, co-infected and many single-infected hosts seem unable to mount effective immune responses. A strong defense may actually be harmful. Analysis of Bd and Bsal secretions documents small metabolites that signal high density to limit their growth and to suppress adaptive immune defenses, thus enabling a stealth presence in the skin compartment.

Quantity and Quality of Aquaculture Enrichments Influence Disease Epidemics and Provide Ecological Alternatives to Antibiotics

Environmental heterogeneity is a central component influencing the virulence and epidemiology of infectious diseases. The number and distribution of susceptible hosts determines disease transmission opportunities, shifting the epidemiological threshold between the spread and fadeout of a disease. Similarly, the presence and diversity of other hosts, pathogens and environmental microbes, may inhibit or accelerate an epidemic. This has important applied implications in farming environments, where high numbers of susceptible hosts are maintained in conditions of minimal environmental heterogeneity. We investigated how the quantity and quality of aquaculture enrichments (few vs. many stones; clean stones vs. stones conditioned in lake water) influenced the severity of infection of a pathogenic bacterium, Flavobacterium columnare, in salmonid fishes. We found that the conditioning of the stones significantly increased host survival in rearing tanks with few stones. A similar effect of increased host survival was also observed with a higher number of unconditioned stones. These results suggest that a simple increase in the heterogeneity of aquaculture environment can significantly reduce the impact of diseases, most likely operating through a reduction in pathogen transmission (stone quantity) and the formation of beneficial microbial communities (stone quality). This supports enriched rearing as an ecological and economic way to prevent bacterial infections with the minimal use of antimicrobials.

Increased tropospheric ozone levels enhance pathogen infection levels of amphibians

As a result of anthropogenic activities, changes to the chemistry of Earth's atmosphere pose a threat to ecosystem health and biodiversity. One such change is the increase in tropospheric ozone (O₃), which is particularly severe in the Mediterranean basin area, where the levels of this pollutant are chronically high during spring and summer time. Within this region, Mediterranean mountain ecosystems are hot spots for biodiversity which may be especially vulnerable to changes in O₃ levels. Declines in montane amphibian populations have been recorded worldwide, including the Mediterranean basin. A significant driver of these declines is the emerging infection disease, chytridiomycosis, caused by the aquatic fungus Batrachochytrium dendrobatidis (Bd). Chytridiomycosis has negatively affected populations of several amphibian species in the Spanish Central Range, including in the Sierra Guadarrama, and interactions with other biotic and abiotic factors are an important part of these declines. However, there is little evidence or knowledge of whether tropospheric O₃ levels may be another factor in the outbreaks of this disease. To test the hypothesis that O₃ levels are another interactive driver of Bd infection dynamics, two different approaches were followed: 1) an experimental study in open top chambers was used to quantify the aspects of how Bd infection progressed throughout the metamorphic process under four different O₃ levels; and 2) a field epidemiological study was used to analyse the relationship between the Bd infection load in the Sierra de Guadarrama and tropospheric O₃ levels during a 9 year period. Our results suggest that high O₃ levels significantly delayed the rate of development of tadpoles and increased Bd infection, providing empirical evidence of two new separate ways that may explain population declines of montane amphibians.

The amphibian microbiome exhibits poor resilience following pathogen-induced disturbance

Infectious pathogens can disrupt the microbiome in addition to directly affecting the host. Impacts of disease may be dependent on the ability of the microbiome to recover from such disturbance, yet remarkably little is known about microbiome recovery after disease, particularly in nonhuman animals. We assessed the resilience of the amphibian skin microbial community after disturbance by the pathogen, Batrachochytrium dendrobatidis (Bd). Skin microbial communities of laboratory-reared mountain yellow-legged frogs were tracked through three experimental phases: prior to Bd infection, after Bd infection (disturbance), and after clearing Bd infection (recovery period). Bd infection disturbed microbiome composition and altered the relative abundances of several dominant bacterial taxa. After Bd infection, frogs were treated with an antifungal drug that cleared Bd infection, but this did not lead to recovery of microbiome composition (measured as Unifrac distance) or relative abundances of dominant bacterial groups. These results indicate that Bd infection can lead to an alternate stable state in the microbiome of sensitive amphibians, or that microbiome recovery is extremely slow—in either case resilience is low. Furthermore, antifungal treatment and clearance of Bd infection had the additional effect of reducing microbial community variability, which we hypothesize results from similarity across frogs in the taxa that colonize community vacancies resulting from the removal of Bd. Our results indicate that the skin microbiota of mountain yellow-legged frogs has low resilience following Bd-induced disturbance and is further altered by the process of clearing Bd infection, which may have implications for the conservation of this endangered amphibian.

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.

Recovery and resiliency of skin microbial communities on the southern leopard frog (Lithobates sphenocephalus) following two biotic disturbances

BACKGROUND

Microorganisms have intimate functional relationships with invertebrate and vertebrate taxa, with the potential to drastically impact health outcomes. Perturbations that affect microbial communities residing on animals can lead to dysbiosis, a change in the functional relationship, often associated with disease. Batrachochytrium dendrobatidis (Bd), a fungal pathogen of amphibians, has been responsible for catastrophic amphibian population declines around the globe. Amphibians harbor a diverse cutaneous microbiome, including some members which are known to be antagonistic to Bd (anti-Bd). Anti-Bd microorganisms facilitate the ability of some frog populations to persist in the presence of Bd, where other populations that lack anti-Bd microorganisms have declined. Research suggests disease-antagonistic properties of the microbiome may be a function of microbial community interactions, rather than individual bacterial species. Conservation efforts have identified amphibian-associated bacteria that exhibit anti-fungal properties for use as 'probiotics' on susceptible amphibian populations. Probiotic application, usually with a single bacterial species, may benefit from a greater understanding of amphibian species-specific microbiome responses to disturbances (e.g. dysbiosis vs. recovery). We assessed microbiome responses to two microbial disturbance events over multiple time points.

RESULTS

Exposing Lithobates sphenocephalus (southern leopard frog) adults to the biopesticidal bacteria Bacillus thuringiensis, followed by exposure to the fungal pathogen Bd, did not have long term impacts on the microbiome. After initial shifts, microbial communities recovered and returned to a state that resembled pre-disturbance.

CONCLUSIONS

Our results indicate microbial communities on L. sphenocephalus are robust and resistant to permanent shifts from some disturbances. This resiliency of microbial communities may explain why L. sphenocephalus is not experiencing the population declines from Bd that impacts many other species. Conservation efforts may benefit from studies outlining amphibian species-specific microbiome responses to disturbances (e.g. dysbiosis vs. recovery). If microbial communities on a threatened amphibian species are unlikely to recover following a disturbance, additional measures may be implemented to ameliorate the impacts of physical and chemical stressors on host-associated microbial communities.

Widespread Pig Farming Practice Linked to Shifts in Skin Microbiomes and Disease in Pond-Breeding Amphibians

Farming practices may reshape the structure of watersheds, water quality, and the health of aquatic organisms. Nutrient enrichment from agricultural pollution increases disease pressure in many host-pathogen systems, but the mechanisms underlying this pattern are not always resolved. For example, nutrient enrichment should strongly influence pools of aquatic environmental bacteria, which has the potential to alter microbiome composition of aquatic animals and their vulnerability to disease. However, shifts in the host microbiome have received little attention as a link between nutrient enrichment and diseases of aquatic organisms. We examined nutrient enrichment through the widespread practice of integrated pig-fish farming and its effects on microbiome composition of Brazilian amphibians and prevalence of the globally distributed amphibian skin pathogen Batrachochytrium dendrobatidis (Bd). This farming system drove surges in fecal coliform bacteria, disturbing amphibian skin bacterial communities such that hosts recruited higher proportions of Bd-facilitative bacteria and carried higher Bd prevalence. Our results highlight previously overlooked connections between global trends in land use change, microbiome dysbiosis, and wildlife disease. These interactions may be particularly important for disease management in the tropics, a region with both high biodiversity and continually intensifying anthropogenic pressures on aquatic wildlife habitats.

Frog Skin Microbiota Vary With Host Species and Environment but Not Chytrid Infection

Describing the structure and function of the amphibian cutaneous microbiome has gained importance with the spread of Batrachochytrium dendrobatidis (Bd), the fungal pathogen that can cause the skin disease chytridiomycosis. Sampling amphibian skin microbiota is needed to characterize current infection status and to help predict future susceptibility to Bd based on microbial composition since some skin microbes have antifungal capabilities that may confer disease resistance. Here, I use 16S rRNA sequencing to describe the composition and structure of the cutaneous microbiota of six species of amphibians. Frog skin samples were also tested for Bd, and I found 11.8% Bd prevalence among all individuals sampled (n = 76). Frog skin microbiota varied by host species and sampling site, but did not differ among Bd-positive and Bd-negative individuals. These results suggest that bacterial composition reflects host species and the environment, but does not reflect Bd infection among the species sampled here. Of the bacterial OTUs identified using an indicator species analysis as strongly associated with amphibians, significantly more indicator OTUs were putative anti-Bd taxa than would be expected based on the proportion of anti-Bd OTUs among all frog OTUs, suggesting strong associations between host species and anti-Bd OTUs. This relationship may partially explain why some of these frogs are asymptomatic carriers of Bd, but more work is needed to determine the other factors that contribute to interspecific variation in Bd susceptibility. This work provides important insights on inter- and intra-specific variation in microbial community composition, putative function, and disease dynamics in populations of amphibians that appear to be coexisting with Bd.

Changes in physiology and microbial diversity in larval ornate chorus frogs are associated with habitat quality

Environmental change associated with anthropogenic disturbance can lower habitat quality, especially for sensitive species such as many amphibians. Variation in environmental quality may affect an organism's physiological health and, ultimately, survival and fitness. Using multiple health measures can aid in identifying populations at increased risk of declines. Our objective was to measure environmental variables at multiple spatial scales and their effect on three indicators of health in ornate chorus frog (Pseudacris ornata) tadpoles to identify potential correlates of population declines. To accomplish this, we measured a glucocorticoid hormone (corticosterone; CORT) profile associated with the stress response, as well as the skin mucosal immune function (combined function of skin secretions and skin bacterial community) and bacterial communities of tadpoles from multiple ponds. We found that water quality characteristics associated with environmental variation, including higher water temperature, conductivity and total dissolved solids, as well as percent developed land nearby, were associated with elevated CORT release rates. However, mucosal immune function, although highly variable, was not significantly associated with water quality or environmental factors. Finally, we examined skin bacterial diversity as it aids in immunity and is affected by environmental variation. We found that skin bacterial diversity differed between ponds and was affected by land cover type, canopy cover and pond proximity. Our results indicate that both local water quality and land cover characteristics are important determinants of population health for ornate chorus frogs. Moreover, using these proactive measures of health over time may aid in early identification of at-risk populations that could prevent further declines and aid in management decisions.

Assessing similarities and disparities in the skin microbiota between wild and laboratory populations of house mice

The house mouse is a key model organism in skin research including host–microbiota interactions, yet little is known about the skin microbiota of free-living mice. It is similarly unclear how closely laboratory mice, which typically live under exceptionally hygienic conditions, resemble the ancestral state of microbial variation in the wild. In this study, we sampled an area spanning 270 km² in south-west France and collected 203 wild Mus musculus domesticus. We profiled the ear skin microbiota on standing and active communities (DNA-based and RNA-based 16 rRNA gene sequencing, respectively), and compared multiple community aspects between wild-caught and laboratory-reared mice kept in distinct facilities. Compared to lab mice, we reveal the skin microbiota of wild mice on the one hand to be unique in their composition within the Staphylococcus genus, with a majority of sequences most closely matching known novobiocin-resistant species, and display evidence of a rare biosphere. On the other hand, despite drastic disparities between natural and laboratory environments, we find that shared taxa nonetheless make up the majority of the core skin microbiota of both wild- and laboratory skin communities, suggesting that mammalian skin is a highly specialized habitat capable of strong selection from available species pools. Finally, the influence of environmental factors suggests RNA-based profiling as a preferred method to reduce environmental noise.

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.

Composition of the cutaneous bacterial community of a cave amphibian, Proteus anguinus

The European cave salamander Proteus anguinus is a charismatic amphibian endemic to the concealed and inaccessible subterranean waters of the Dinaric Karst. Despite its exceptional conservation importance not much is known about its ecology and interactions with the groundwater microbiome. The cutaneous microbiota of amphibians is an important driver of metabolic capabilities and immunity, and thus a key factor in their wellbeing and survival. We used high-throughput 16S rRNA gene sequencing based on seven variable regions to examine the bacteriome of the skin of five distinct evolutionary lineages of P. anguinus and in their groundwater environment. The skin bacteriomes turned out to be strongly filtered subsamples of the environmental microbial community. The resident microbiota of the analyzed individuals was dominated by five bacterial taxa. Despite an indicated functional redundancy, the cutaneous bacteriome of P. anguinus presumably provides protection against invading microbes by occupying the niche, and thus could serve as an indicator of health status. Besides conservation implications for P. anguinus, our results provide a baseline for future studies on other endangered neotenic salamanders.

Microbial associations and spatial proximity predict North American moose (Alces alces) gastrointestinal community composition

Microbial communities are increasingly recognized as crucial for animal health. However, our understanding of how microbial communities are structured across wildlife populations is poor. Mechanisms such as interspecific associations are important in structuring free-living communities, but we still lack an understanding of how important interspecific associations are in structuring gut microbial communities in comparison with other factors such as host characteristics or spatial proximity of hosts. Here, we ask how gut microbial communities are structured in a population of North American moose Alces alces. We identify key microbial interspecific associations within the moose gut and quantify how important they are relative to key host characteristics, such as body condition, for predicting microbial community composition. We sampled gut microbial communities from 55 moose in a population experiencing decline due to a myriad of factors, including pathogens and malnutrition. We examined microbial community dynamics in this population utilizing novel graphical network models that can explicitly incorporate spatial information. We found that interspecific associations were the most important mechanism structuring gut microbial communities in moose and detected both positive and negative associations. Models only accounting for associations between microbes had higher predictive value compared to models including moose sex, evidence of previous pathogen exposure or body condition. Adding spatial information on moose location further strengthened our model and allowed us to predict microbe occurrences with ~90% accuracy. Collectively, our results suggest that microbial interspecific associations coupled with host spatial proximity are vital in shaping gut microbial communities in a large herbivore. In this case, previous pathogen exposure and moose body condition were not as important in predicting gut microbial community composition. The approach applied here can be used to quantify interspecific associations and gain a more nuanced understanding of the spatial and host factors shaping microbial communities in non-model hosts.

Diversity-Stability Dynamics of the Amphibian Skin Microbiome and Susceptibility to a Lethal Viral Pathogen

Variation among animals in their host-associated microbial communities is increasingly recognized as a key determinant of important life history traits including growth, metabolism, and resistance to disease. Quantitative estimates of the factors shaping the stability of host microbiomes over time at the individual level in non-model organisms are scarce. Addressing this gap in our knowledge is important, as variation among individuals in microbiome stability may represent temporal gain or loss of key microbial species and functions linked to host health and/or fitness. Here we use controlled experiments to investigate how both heterogeneity in microbial species richness of the environment and exposure to the emerging pathogen Ranavirus influence the structure and temporal dynamics of the skin microbiome in a vertebrate host, the European common frog (Rana temporaria). Our evidence suggests that altering the bacterial species richness of the environment drives divergent temporal microbiome dynamics of the amphibian skin. Exposure to ranavirus effects changes in skin microbiome structure irrespective of total microbial diversity, but individuals with higher pre-exposure skin microbiome diversity appeared to exhibit higher survival. Higher diversity skin microbiomes also appear less stable over time compared to lower diversity microbiomes, but stability of the 100 most abundant ("core") community members was similar irrespective of microbiome richness. Our study highlights the importance of extrinsic factors in determining the stability of host microbiomes over time, which may in turn have important consequences for the stability of host-microbe interactions and microbiome-fitness correlations.

Fungal infection alters the selection, dispersal and drift processes structuring the amphibian skin microbiome

Symbiotic microbial communities are important for host health, but the processes shaping these communities are poorly understood. Understanding how community assembly processes jointly affect microbial community composition is limited because inflexible community models rely on rejecting dispersal and drift before considering selection. We developed a flexible community assembly model based on neutral theory to ask: How do dispersal, drift and selection concurrently affect the microbiome across environmental gradients? We applied this approach to examine how a fungal pathogen affected the assembly processes structuring the amphibian skin microbiome. We found that the rejection of neutrality for the amphibian microbiome across a fungal gradient was not strictly due to selection processes, but was also a result of species-specific changes in dispersal and drift. Our modelling framework brings the qualitative recognition that niche and neutral processes jointly structure microbiomes into quantitative focus, allowing for improved predictions of microbial community turnover across environmental gradients.

Amphibian skin-associated Pigmentiphaga: Genome sequence and occurrence across geography and hosts

The bacterial communities colonizing amphibian skin have been intensively studied due to their interactions with pathogenic chytrid fungi that are causing drastic amphibian population declines. Bacteria of the family Alcaligenaceae, and more specifically of the genus Pigmentiphaga, have been found to be associated specifically to arboreal frogs. Here we analyze their occurrence in a previously assembled global skin microbiome dataset from 205 amphibian species. Pigmentiphaga made up about 5% of the total number of reads in this global dataset. They were mostly found in unrelated arboreal frogs from Madagascar (Mantellidae and Hyperoliidae), but also occurred at low abundances on Neotropical frogs. Based on their 16S sequences, most of the sequences belong to a clade within Pigmentiphaga not assignable to any type strains of the five described species of the genus. One isolate from Madagascar clustered with Pigmentiphaga aceris (>99% sequence similarity on 16S rRNA gene level). Here, we report the full genome sequence of this bacterium which, based on 16S sequences of >97% similarity, has previously been found on human skin, floral nectar, tree sap, stream sediment and soil. Its genome consists of a single circular chromosome with 6,165,255 bp, 5,300 predicted coding sequences, 57 tRNA genes, and three rRNA operons. In comparison with other known Pigmentiphaga genomes it encodes a higher number of genes associated with environmental information processing and cellular processes. Furthermore, it has a biosynthetic gene cluster for a nonribosomal peptide syntethase, and bacteriocin biosynthetic genes can be found, but clusters for β-lactones present in other comparative Pigmentiphaga genomes are lacking.