Community Structure and Function of Amphibian Skin Microbes: An Experiment with Bullfrogs Exposed to a Chytrid Fungus

Dual stressors of infection and warming can destabilize host microbiomes

Climate change is causing extreme heating events and intensifying infectious disease outbreaks. Animals harbour microbial communities, which are vital for their survival and fitness under stressful conditions. Understanding how microbiome structures change in response to infection and warming may be important for forecasting host performance under global change. Here, we evaluated alterations in the microbiomes of several wild Caenorhabditis elegans isolates spanning a range of latitudes, upon warming temperatures and infection by the parasite Leucobacter musarum. Using 16S rRNA sequencing, we found that microbiome diversity decreased, and dispersion increased over time, with the former being more prominent in uninfected adults and the latter aggravated by infection. Infection reduced dominance of specific microbial taxa, and increased microbiome dispersion, indicating destabilizing effects on host microbial communities. Exposing infected hosts to warming did not have an additive destabilizing effect on their microbiomes. Moreover, warming during pre-adult development alleviated the destabilizing effects of infection on host microbiomes. These results revealed an opposing interaction between biotic and abiotic factors on microbiome structure. Lastly, we showed that increased microbiome dispersion might be associated with decreased variability in microbial species interaction strength. Overall, these findings improve our understanding of animal microbiome dynamics amidst concurrent climate change and epidemics. This article is part of the theme issue 'Sculpting the microbiome: how host factors determine and respond to microbial colonization'.

Baseline Skin Microbiota of the Leatherback Sea Turtle

The integumentary system of the leatherback sea turtle (Dermochelys coriacea) is the most visible and defining difference of the species, with its smooth and waxy carapace and finely scaled skin, distinguishing it from the other six sea turtle species. The skin is the body's largest organ and serves as a primary defense against the outside world and is thus essential to health. To date, we have begun to understand that the microorganisms located on the skin aid in these functions. However, many host-microbial interactions are not yet fully defined or understood. Prior to uncovering these crucial host-microbial interactions, we must first understand the communities of microorganisms present and how they differ through life-stage classes and across the body. Here, we present a comprehensive bacterial microbial profile on the skin of leatherbacks. Using next-generation sequencing (NGS), we identified the major groups of bacteria on the skin of neonates at emergence, neonates at 3-4 weeks of age (i.e., post-hatchlings), and nesting females. These data show that the predominant bacteria on the skin of the leatherback are different at each life-stage class sampled. This suggests that there is a shift in the microbial communities of the skin associated with life-stage class or even possibly age. We also found that different sample locations on the nesting female (i.e., carapace and front appendages = flipper) have significantly different communities of bacteria present. This is likely due to differences in the microhabitats of these anatomic locations and future studies should explore if this variation also holds true for neonates. These data define baseline skin microbiota on the leatherback and can serve as a foundation for additional work to broaden our understanding of the leatherbacks' host-microbial interactions, the impacts of environmental changes or stressors over time, and even the pathogenicity of disease processes.

Skin microbiome disturbance linked to drought-associated amphibian disease

The onset of global climate change has led to abnormal rainfall patterns, disrupting associations between wildlife and their symbiotic microorganisms. We monitored a population of pumpkin toadlets and their skin bacteria in the Brazilian Atlantic Forest during a drought. Given the recognized ability of some amphibian skin bacteria to inhibit the widespread fungal pathogen Batrachochytrium dendrobatidis (Bd), we investigated links between skin microbiome health, susceptibility to Bd and host mortality during a die-off event. We found that rainfall deficit was an indirect predictor of Bd loads through microbiome disruption, while its direct effect on Bd was weak. The microbiome was characterized by fewer putative Bd-inhibitory bacteria following the drought, which points to a one-month lagged effect of drought on the microbiome that may have increased toadlet susceptibility to Bd. Our study underscores the capacity of rainfall variability to disturb complex host-microbiome interactions and alter wildlife disease dynamics.

Body size mediates latitudinal population differences in the response to chytrid fungus infection in two amphibians

Factors behind intraspecific variation in sensitivity to pathogens remain poorly understood. We investigated how geographical origin in two North European amphibians affects tolerance to infection by the chytrid fungus Batrachochytrium dendrobatidis (Bd), a generalist pathogen which has caused amphibian population declines worldwide. We exposed newly metamorphosed individuals of moor frog Rana arvalis and common toad Bufo bufo from two latitudinal regions to two different BdGPL strains. We measured survival and growth as infections may cause sub-lethal effects in fitness components even in the absence of mortality. Infection loads were higher in B. bufo than in R. arvalis, and smaller individuals had generally higher infection loads. B. bufo had high mortality in response to Bd infection, whereas there was little mortality in R. arvalis. Bd-mediated mortality was size-dependent and high-latitude individuals were smaller leading to high mortality in the northern B. bufo. Bd exposure led to sub-lethal effects in terms of reduced growth suggesting that individuals surviving the infection may have reduced fitness mediated by smaller body size. In both host species, the Swedish Bd strain caused stronger sublethal effects than the British strain. We suggest that high-latitude populations can be more vulnerable to chytrids than those from lower latitudes and discuss the possible mechanisms how body size and host geographical origin contribute to the present results.

Climate match is key to predict range expansion of the world's worst invasive terrestrial vertebrates

Understanding the determinants of the range expansion of invasive alien species is crucial for developing effective prevention and control strategies. Nevertheless, we still lack a global picture of the potential factors influencing the invaded range expansion across taxonomic groups, especially for the world's worst invaders with high ecological and economic impacts. Here, by extensively collecting data on 363 distributional ranges of 19 of world's worst invasive terrestrial vertebrates across 135 invaded administrative jurisdictions, we observed remarkable variations in the range expansion across species and taxonomic groups. After controlling for taxonomic and geographic pseudoreplicates, model averaging analyses based on generalized additive mixed-effect models showed that species in invaded regions having climates more similar to those of their native ranges tended to undergo a larger range expansion. In addition, as proxies of propagule pressure and human-assisted transportation, the number of introduction events and the road network density were also important predictors facilitating the range expansion. Further variance partitioning analyses validated the predominant role of climate match in explaining the range expansion. Our study demonstrated that regions with similar climates to their native ranges could still be prioritized to prevent the spread of invasive species under the sustained global change.

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.

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.

Comparative genomics of bacteria from amphibian skin associated with inhibition of an amphibian fungal pathogen, Batrachochytrium dendrobatidis

Chytridiomycosis, caused by Batrachochytrium dendrobatidis (Bd), is a skin disease associated with worldwide amphibian declines. Symbiotic microbes living on amphibian skin interact with Bd and may alter infection outcomes. We completed whole genome sequencing of 40 bacterial isolates cultured from the skin of four amphibian species in the Eastern US. Each isolate was tested in vitro for the ability to inhibit Bd growth. The aim of this study was to identify genomic differences among the isolates and generate hypotheses about the genomic underpinnings of Bd growth inhibition. We identified sixty-five gene families that were present in all 40 isolates. Screening for common biosynthetic gene clusters revealed that this set of isolates contained a wide variety of clusters; the two most abundant clusters with potential antifungal activity were siderophores (N=17 isolates) and Type III polyketide synthases (N=22 isolates). We then examined various subsets of the 22 isolates in the phylum Proteobacteria for genes encoding specific compounds that may inhibit fungal growth, including chitinase and violacein. We identified differences in Agrobacterium and Sphingomonas isolates in the chitinase genes that showed some association with anti-Bd activity, as well as variation in the violacein genes in the Janthinobacterium isolates. Using a comparative genomics approach, we generated several testable hypotheses about differences among bacterial isolates from amphibian skin communities that could contribute to variation in the ability to inhibit Bd growth. Further work is necessary to explore and uncover the various mechanisms utilized by amphibian skin bacterial isolates to inhibit Bd.

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

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

The adaptive microbiome hypothesis and immune interactions in amphibian mucus

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

Metabolite compositions on skins of eastern hellbenders Cryptobranchus alleganiensis alleganiensis differ with location and captivity

Eastern hellbenders Cryptobranchus alleganiensis alleganiensis, large aquatic salamanders, are declining over most of their range. The amphibian-killing fungus Batrachochytrium dendrobatidis (Bd) has contributed to global amphibian declines and has been detected on eastern hellbenders, but infection intensities were lower than those of species that are more susceptible to Bd. The factors limiting Bd on hellbenders may include antifungal metabolites produced by their skin microbiota. We used a metabolite fingerprinting technique to noninvasively identify the presence, but not identity, of metabolites associated with eastern hellbenders. We surveyed the skin of wild eastern hellbenders to test whether the composition and richness (i.e. number of metabolites) of their metabolites are explained by Bd status or location. Furthermore, we surveyed for metabolites on captive eastern hellbenders to test whether metabolite compositions were different between captive and wild eastern hellbenders. Bd detection was not associated with either metabolite richness or composition. Both metabolite composition and richness differed significantly on hellbenders from different locations (i.e. states). For metabolite composition, there was a statistical interaction between location and Bd status. Metabolite richness was greater on captive eastern hellbenders compared to wild hellbenders, and metabolite compositions differed between wild and captive eastern hellbenders. The methods we employed to detect metabolite profiles effectively grouped individuals by location even though metabolite composition and richness have high levels of intraspecific variation. Understanding the drivers and functional consequences of assemblages of skin metabolites on amphibian health will be an important step toward understanding the mechanisms that result in disease vulnerability.

Strong restructuration of skin microbiota during captivity challenges ex-situ conservation of amphibians

In response to the current worldwide amphibian extinction crisis, conservation instances have encouraged the establishment of ex-situ collections for endangered species. The resulting assurance populations are managed under strict biosecure protocols, often involving artificial cycles of temperature and humidity to induce active and overwintering phases, which likely affect the bacterial symbionts living on the amphibian skin. However, the skin microbiota is an important first line of defense against pathogens that can cause amphibian declines, such as the chytrid Batrachochytrium dendrobatidis (Bd). Determining whether current husbandry practices for assurance populations might deplete amphibians from their symbionts is therefore essential to conservation success. Here, we characterize the effect of the transitions from the wild to captivity, and between aquatic and overwintering phases, on the skin microbiota of two newt species. While our results confirm differential selectivity of skin microbiota between species, they underscore that captivity and phase-shifts similarly affect their community structure. More specifically, the translocation ex-situ is associated with rapid impoverishment, decrease in alpha diversity and strong species turnover of bacterial communities. Shifts between active and overwintering phases also cause changes in the diversity and composition of the microbiota, and on the prevalence of Bd-inhibitory phylotypes. Altogether, our results suggest that current husbandry practices strongly restructure the amphibian skin microbiota. Although it remains to be determined whether these changes are reversible or have deleterious effects on their hosts, we discuss methods to limit microbial diversity loss ex-situ and emphasize the importance of integrating bacterial communities to applied amphibian conservation.

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.

What lives on and in the sea turtle? A literature review of sea turtle bacterial microbiota

Within the United States, all populations of sea turtles are listed as threatened or endangered under the Endangered Species Act. Identifying methods of supporting health, preventing disease, and treating disease is essential for conservation and management strategies for all species. Over the last few decades, advances in technology and high throughput sequencing have allowed a proliferation of research into core microbiota and microbiomes in humans and animals. Such investigations have proven that microbiota on and within a host can influence physiology, immunity, and development. Accordingly, a comprehensive understanding of microbiota is essential for unearthing the complex relationships within a microbiome and how those interactions and relationships can be used to promote health and prevent or treat disease. The goal of this review is to summarize the current microbiota research available across all species of sea turtles and identify any emerging trends. Methodological differences made it challenging to draw conclusions across studies, but it is apparent that each anatomical location investigated has a unique core microbiota with some potential overlap. In the future, unifying methodology across microbiota studies will allow broader conclusions to be drawn across all anatomic locations and species of sea turtles. These conclusions will then allow clinicians and conservationists to apply the research results in the field. Additionally, future efforts should include a wider range of organisms including fungi, viruses, parasites, epibiota, and archaea to unveil essential relationships among and between the organisms and host for maintenance of a healthy microbiome.

An experimental test of disease resistance function in the skin-associated bacterial communities of three tropical amphibian species

Variation in the structure of host-associated microbial communities has been correlated with the occurrence and severity of disease in diverse host taxa, suggesting a key role of the microbiome in pathogen defense. However, whether these correlations are typically a cause or consequence of pathogen exposure remains an open question, and requires experimental approaches to disentangle. In amphibians, infection by the fungal pathogen Batrachochytrium dendrobatidis (Bd) alters the skin microbial community in some host species, whereas in other species, the skin microbial community appears to mediate infection dynamics. In this study, we completed experimental Bd exposures in three species of tropical frogs (Agalychnis callidryas, Dendropsophus ebraccatus, Craugastor fitzingeri) that were sympatric with Bd at the time of the study. For all three species, we identified key taxa within the skin bacterial communities that were linked to Bd infection dynamics. We also measured higher Bd infection intensities in D. ebraccatus and C. fitzingeri that were associated with higher mortality in C. fitzingeri. Our findings indicate that microbially-mediated pathogen resistance is a complex trait that can vary within and across host species, and suggest that symbiont communities that have experienced prior selection for defensive microbes may be less likely to be disturbed by pathogen exposure.

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.

Skin bacterial metacommunities of San Francisco Bay Area salamanders are structured by host genus and habitat quality

Host-associated microbial communities can influence physiological processes of macroorganisms, including contributing to infectious disease resistance. For instance, some bacteria that live on amphibian skin produce antifungal compounds that inhibit two lethal fungal pathogens, Batrachochytrium dendrobatidis (Bd) and B. salamandrivorans (Bsal). Therefore, differences in microbiome composition among host species or populations within a species can contribute to variation in susceptibility to Bd/Bsal. This study applies 16S rRNA sequencing to characterize the skin bacterial microbiomes of three widespread terrestrial salamander genera native to the western United States. Using a metacommunity structure analysis, we identified dispersal barriers for these influential bacteria between salamander families and localities. We also analyzed the effects of habitat characteristics such as percent natural cover and temperature seasonality on the microbiome. We found that certain environmental variables may influence the skin microbial communities of some salamander genera more strongly than others. Each salamander family had a somewhat distinct community of putative anti-Bd skin bacteria, suggesting that salamanders may select for a functional assembly of cutaneous symbionts that could differ in its ability to protect these amphibians from disease. Our observations raise the need to consider host identity and environmental heterogeneity during the selection of probiotics to treat wildlife diseases.

Experimental test of microbiome protection across pathogen doses reveals importance of resident microbiome composition

The commensal microbes inhabiting a host tissue can interact with invading pathogens and host physiology in ways that alter pathogen growth and disease manifestation. Prior work in house finches (Haemorhous mexicanus) found that resident ocular microbiomes were protective against conjunctival infection and disease caused by a relatively high dose of Mycoplasma gallisepticum. Here, we used wild-caught house finches to experimentally examine whether protective effects of the resident ocular microbiome vary with the dose of invading pathogen. We hypothesized that commensal protection would be strongest at low M. gallisepticum inoculation doses because the resident microbiome would be less disrupted by invading pathogen. Our five M. gallisepticum dose treatments were fully factorial with an antibiotic treatment to perturb resident microbes just prior to M. gallisepticum inoculation. Unexpectedly, we found no indication of protective effects of the resident microbiome at any pathogen inoculation dose, which was inconsistent with the prior work. The ocular bacterial communities at the beginning of our experiment differed significantly from those previously reported in local wild-caught house finches, likely causing this discrepancy. These variable results underscore that microbiome-based protection in natural systems can be context dependent, and natural variation in community composition may alter the function of resident microbiomes in free-living animals.

Genetically modifying skin microbe to produce violacein and augmenting microbiome did not defend Panamanian golden frogs from disease

We designed two probiotic treatments to control chytridiomycosis caused by Batrachochytrium dendrobatidis (Bd) on infected Panamanian golden frogs (Atelopus zeteki), a species that is thought to be extinct in the wild due to Bd. The first approach disrupted the existing skin microbe community with antibiotics then exposed the frogs to a core golden frog skin microbe (Diaphorobacter sp.) that we genetically modified to produce high titers of violacein, a known antifungal compound. One day following probiotic treatment, the engineered Diaphorobacter and the violacein-producing pathway could be detected on the frogs but the treatment failed to improve frog survival when exposed to Bd. The second approach exposed frogs to the genetically modified bacterium mixed into a consortium with six other known anti-Bd bacteria isolated from captive A. zeteki, with no preliminary antibiotic treatment. The consortium treatment increased the frequency and abundance of three probiotic isolates (Janthinobacterium, Chryseobacterium, and Stenotrophomonas) and these persisted on the skin 4 weeks after probiotic treatment. There was a temporary increase in the frequency and abundance of three other probiotics isolates (Masillia, Serratia, and Pseudomonas) and the engineered Diaphorobacter isolate, but they subsequently disappeared from the skin. This treatment also failed to reduce frog mortality upon exposure.

Carryover effects minimized the positive effects of treated wastewater on anuran development

Constructed wetlands (CWs) are a potential solution for wastewater treatment due to their capacity to support native species and provide tertiary wastewater treatment. However, CWs can expose wildlife communities to excess nutrients and harmful contaminants, affecting their development, morphology, and behavior. To examine how wastewater CWs may affect wildlife, we raised Southern leopard frogs, Lithobates sphenocephalus, in wastewater from conventional secondary lagoon and tertiary CW treatments for comparison with pondwater along with the presence and absence of a common plant invader to these systems - common duckweed (Lemna minor) - and monitored their juvenile development for potential carryover effects into the terrestrial environment. The tertiary CW treatment did not change demographic or morphological outcomes relative to conventional wastewater treatment in our study. Individuals emerging from both wastewater treatments demonstrated lower terrestrial survival rates than those emerging from pondwater throughout the experiment though experiment-wide survival rates were equivalent among treatments. Individuals from wastewater treatments transformed at larger sizes relative to those in pondwater, but this advantage was minimized in the terrestrial environment. Individuals that developed with duckweed had consistent but marginally better performance in both environments. Our results suggest a potential trade-off between short-term benefits of development in treated effluent and long-term consequences on overall fitness. Overall, we demonstrate that CWs for the purpose of wastewater treatment may not be suitable replicates for wildlife habitat and could have consequences for local population dynamics.

Immune Gene Expression Covaries with Gut Microbiome Composition in Stickleback

Commensal microbial communities have immense effects on their vertebrate hosts, contributing to a number of physiological functions, as well as host fitness. In particular, host immunity is strongly linked to microbiota composition through poorly understood bi-directional links. Gene expression may be a potential mediator of these links between microbial communities and host function. However, few studies have investigated connections between microbiota composition and expression of host immune genes in complex systems. Here, we leverage a large study of laboratory-raised fish from the species Gasterosteus aculeatus (three-spined stickleback) to document correlations between gene expression and microbiome composition. First, we examined correlations between microbiome alpha diversity and gene expression. Our results demonstrate robust positive associations between microbial alpha diversity and expression of host immune genes. Next, we examined correlations between host gene expression and abundance of microbial taxa. We identified 15 microbial families that were highly correlated with host gene expression. These families were all tightly correlated with host expression of immune genes and processes, falling into one of three categories—those positively correlated, negatively correlated, and neutrally related to immune processes. Furthermore, we highlight several important immune processes that are commonly associated with the abundance of these taxa, including both macrophage and B cell functions. Further functional characterization of microbial taxa will help disentangle the mechanisms of the correlations described here. In sum, our study supports prevailing hypotheses of intimate links between host immunity and gut microbiome composition.

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

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

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.

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.

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

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

Cross-Disciplinary Genomics Approaches to Studying Emerging Fungal Infections

Emerging fungal pathogens pose a serious, global and growing threat to food supply systems, wild ecosystems, and human health. However, historic chronic underinvestment in their research has resulted in a limited understanding of their epidemiology relative to bacterial and viral pathogens. Therefore, the untargeted nature of genomics and, more widely, -omics approaches is particularly attractive in addressing the threats posed by and illuminating the biology of these pathogens. Typically, research into plant, human and wildlife mycoses have been largely separated, with limited dialogue between disciplines. However, many serious mycoses facing the world today have common traits irrespective of host species, such as plastic genomes; wide host ranges; large population sizes and an ability to persist outside the host. These commonalities mean that -omics approaches that have been productively applied in one sphere and may also provide important insights in others, where these approaches may have historically been underutilised. In this review, we consider the advances made with genomics approaches in the fields of plant pathology, human medicine and wildlife health and the progress made in linking genomes to other -omics datatypes and sets; we identify the current barriers to linking -omics approaches and how these are being underutilised in each field; and we consider how and which -omics methodologies it is most crucial to build capacity for in the near future.

Chytrid fungi shape bacterial communities on model particulate organic matter

Microbial colonization and degradation of particulate organic matter (POM) are important processes that influence the structure and function of aquatic ecosystems. Although POM is readily used by aquatic fungi and bacteria, there is a limited understanding of POM-associated interactions between these taxa, particularly for early-diverging fungal lineages. Using a model ecological system with the chitin-degrading freshwater chytrid fungus Rhizoclosmatium globosum and chitin microbeads, we assessed the impacts of chytrid fungi on POM-associated bacteria. We show that the presence of chytrids on POM alters concomitant bacterial community diversity and structure, including differing responses between chytrid life stages. We propose that chytrids can act as ecosystem facilitators through saprotrophic feeding by producing ‘public goods’ from POM degradation that modify bacterial POM communities. This study suggests that chytrid fungi have complex ecological roles in aquatic POM degradation not previously considered, including the regulation of bacterial colonization, community succession and subsequent biogeochemical potential.

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.

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.

Columbia spotted frogs (Rana luteiventris) have characteristic skin microbiota that may be shaped by cutaneous skin peptides and the environment

Global amphibian declines due to the fungal pathogen Batrachochytrium dendrobatidis (Bd) have led to questions about how amphibians defend themselves against skin diseases. A total of two amphibian defense mechanisms are antimicrobial peptides (AMPs), a component of amphibian innate immune defense and symbiotic skin bacteria, which can act in synergy. We characterized components of these factors in four populations of Columbia spotted frogs (Rana luteiventris) to investigate their role in disease defense. We surveyed the ability of their AMPs to inhibit Bd, skin bacterial community composition, skin metabolite profiles and presence and intensity of Bd infection. We found that AMPs from R. luteiventris inhibited Bd in bioassays, but inhibition did not correlate with Bd intensity on frogs. R. luteiventris had two prevalent and abundant core bacteria: Rhizobacter and Chryseobacterium. Rhizobacter relative abundance was negatively correlated with AMP's ability to inhibit Bd, but was not associated with Bd status itself. There was no relationship between metabolites and Bd. Bacterial communities and Bd differ by location, which suggests a strong environmental influence. R. luteiventris are dominated by consistent core bacteria, but also house transient bacteria that are site specific. Our emergent hypothesis is that host control and environmental factors shape the microbiota on R. luteiventris.

Antifungal Potential of the Skin Microbiota of Hibernating Big Brown Bats (Eptesicus fuscus) Infected With the Causal Agent of White-Nose Syndrome

Little is known about skin microbiota in the context of the disease white-nose syndrome (WNS), caused by the fungus Pseudogymnoascus destructans (Pd), that has caused enormous declines of hibernating North American bats over the past decade. Interestingly, some hibernating species, such as the big brown bat (Eptesicus fuscus), appear resistant to the disease and their skin microbiota could play a role. However, a comprehensive analysis of the skin microbiota of E. fuscus in the context of Pd has not been done. In January 2017, we captured hibernating E. fuscus, sampled their skin microbiota, and inoculated them with Pd or sham inoculum. We allowed the bats to hibernate in the lab under controlled conditions for 11 weeks and then sampled their skin microbiota to test the following hypotheses: (1) Pd infection would not disrupt the skin microbiota of Pd-resistant E. fuscus; and (2) microbial taxa with antifungal properties would be abundant both before and after inoculation with Pd. Using high-throughput 16S rRNA gene sequencing, we discovered that beta diversity of Pd-inoculated bats changed more over time than that of sham-inoculated bats. Still, the most abundant taxa in the community were stable throughout the experiment. Among the most abundant taxa, Pseudomonas and Rhodococcus are known for antifungal potential against Pd and other fungi. Thus, in contrast to hypothesis 1, Pd infection destabilized the skin microbiota but consistent with hypothesis 2, bacteria with known antifungal properties remained abundant and stable on the skin. This study is the first to provide a comprehensive survey of skin microbiota of E. fuscus, suggesting potential associations between the bat skin microbiota and resistance to the Pd infection and WNS. These results set the stage for future studies to characterize microbiota gene expression, better understand mechanisms of resistance to WNS, and help develop conservation strategies.

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.

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.

Predicting Microbiome Function Across Space Is Confounded by Strain-Level Differences and Functional Redundancy Across Taxa

Variation in the microbiome among individual organisms may play a critical role in the relative susceptibility of those organisms to infection, disease, and death. However, predicting microbiome function is difficult because of spatial and temporal variation in microbial diversity, and taxonomic diversity is not predictive of microbiome functional diversity. Addressing this issue may be particularly important when addressing pandemic diseases, such as the global amphibian die-off associated with Bd. Some of the most important factors in probiotic development for disease treatment are whether bacteria with desired function can be found on native amphibians in the local environment. To address this issue, we isolated, sequenced, and assayed the cutaneous bacterial communities of Plethodon cinereus along a gradient of land use change. Our results suggest that cutaneous community composition, but not overall diversity, change with changes in land use, but this does not correspond to significant change in Bd-inhibitory function. We found that Bd-inhibition is a functionally redundant trait, but that level of inhibition varies over phylogenetic, spatial, and temporal scales. This research provides further evidence for the importance of continued examination of amphibian microbial communities across environmental gradients, including biotic and abiotic interactions, when considering disease dynamics.

Functional Redundancy of Batrachochytrium dendrobatidis Inhibition in Bacterial Communities Isolated from Lithobates clamitans Skin

The cutaneous microbial community can influence the health of amphibians exposed to Batrachochytrium dendrobatidis (Bd), a fungal pathogen that has contributed to recent amphibian declines. Resistance to Bd in amphibian populations is correlated with the presence of anti-Bd cutaneous microbes, which confer disease resistance by inhibiting Bd growth. I aimed to determine if green frogs (Lithobates clamitans), an abundant and widely distributed species in New Jersey, harbored bacteria that inhibit Bd and whether the presence and identity of these microbes varied among sites. I used in vitro challenge assays to determine if bacteria isolated from green frog skin could inhibit or enhance the growth of Bd. I found that green frogs at all sites harbored anti-Bd bacteria. However, there were differences in Bd inhibition capabilities among bacterial isolates identified as the same operational taxonomic unit (OTU), lending support to the idea that phylogenetic relatedness does not always predict Bd inhibition status. Additionally, anti-Bd bacterial richness did not vary by site, but the composition of anti-Bd bacterial taxa was distinct at each site. This suggests that there is functional redundancy of Bd inhibition across unique communities of anti-Bd symbionts found on frogs at different sites. These findings highlight the need to better elucidate the structure-function relationship of microbiomes and their role in disease resistance.

Amphibian Host and Skin Microbiota Response to a Common Agricultural Antimicrobial and Internal Parasite

Holistic approaches that simultaneously characterize responses of both microbial symbionts and their hosts to environmental shifts are imperative to understanding the role of microbiotas on host health. Using the northern leopard frog (Lithobates pipiens) as our model, we investigated the effects of a common trematode (family Echinostomatidae), a common agricultural antimicrobial (Sulfadimethoxine; SDM), and their interaction on amphibian skin microbiota and amphibian health (growth metrics and susceptibility to parasites). In the trematode-exposed individuals, we noted an increase in alpha diversity and a shift in microbial communities. In the SDM-treated individuals, we found a change in the composition of the skin microbiota similar to those induced by the trematode treatment. Groups treated with SDM, echinostomes, or a combination of SDM and echinostomes, had higher relative abundances of OTUs assigned to Flavobacterium and Acinetobacter. Both of these genera have been associated with infectious disease in amphibians and the production of anti-pathogen metabolites. Similar changes in microbial community composition between SDM and trematode exposed individuals may have resulted from stress-related disruption of host immunity. Despite changes in the microbiota, we found no effect of echinostomes and SDM on host health. Given the current disease- and pollution-related threats facing amphibians, our study highlights the need to continue to evaluate the influence of natural and anthropogenic stressors on host-associated microbial communities.

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.

Ecological Correlates of Large-Scale Turnover in the Dominant Members of Pseudacris crucifer Skin Bacterial Communities

Animals host a wide diversity of symbiotic microorganisms that contribute important functions to host health, and our knowledge of what drives variation in the composition of these complex communities continues to grow. Microbiome studies at larger spatial scales present opportunities to evaluate the contribution of large-scale factors to variation in the microbiome. We conducted a large-scale field study to assess variation in the bacterial symbiont communities on adult frog skin (Pseudacris crucifer), characterized using 16S rRNA gene amplicon sequencing. We found that skin bacterial communities on frogs were less diverse than, and structurally distinct from, the surrounding habitat. Frog skin was typically dominated by one of two bacterial OTUs: at western sites, a Proteobacteria dominated the community, whereas eastern sites were dominated by an Actinobacteria. Using a metacommunity framework, we then sought to identify factors explaining small- and large-scale variation in community structure-that is, among hosts within a pond, and among ponds spanning the study transect. We focused on the presence of a fungal skin pathogen, Batrachochytrium dendrobatidis (Bd) as one potential driver of variation. We found no direct link between skin bacterial community structure and Bd infection status of individual frog hosts. Differences in pond-level community structure, however, were explained by Bd infection prevalence. Importantly, Bd infection prevalence itself was correlated with numerous other environmental factors; thus, skin bacterial diversity may be influenced by a complex suite of extrinsic factors. Our findings indicate that large-scale factors and processes merit consideration when seeking to understand microbiome diversity.

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

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

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.

Outbreaks of an Emerging Viral Disease Covary With Differences in the Composition of the Skin Microbiome of a Wild United Kingdom Amphibian

There is growing appreciation of the important role of commensal microbes in ensuring the normal function and health of their hosts, including determining how hosts respond to pathogens. A range of infectious diseases are threatening amphibians worldwide, and evidence is accumulating that the host-associated bacteria that comprise the microbiome may be key in mediating interactions between amphibian hosts and infectious pathogens. We used 16S rRNA amplicon sequencing to quantify the skin microbial community structure of over 200 individual wild adult European common frogs (Rana temporaria), from ten populations with contrasting history of the lethal disease ranavirosis, caused by emerging viral pathogens belonging to the genus Ranavirus. All populations had similar species richness irrespective of disease history, but populations that have experienced historical outbreaks of ranavirosis have a distinct skin microbiome structure (beta diversity) when compared to sites where no outbreaks of the disease have occurred. At the individual level, neither age, body length, nor sex of the frog could predict the structure of the skin microbiota. Our data potentially support the hypothesis that variation among individuals in skin microbiome structure drive differences in susceptibility to infection and lethal outbreaks of disease. More generally, our results suggest that population-level processes are more important for driving differences in microbiome structure than variation among individuals within populations in key life history traits such as age and body size.

Integrating the role of antifungal bacteria into skin symbiotic communities of three Neotropical frog species

Chytridiomycosis, caused by the pathogen Batrachochytrium dendrobatidis (Bd), has led to population declines and extinctions of frog species around the world. While it is known that symbiotic skin bacteria can play a protective role against pathogens, it is not known how these defensive bacteria are integrated into the bacterial community on amphibian skin. In this study, we used 16S rRNA gene amplicon sequencing, culturing and Bd inhibition bioassays to characterize the communities of skin bacteria on three Neotropical frog species that persist in a Bd-infected area in Panama and determined the abundance and integration of anti-Bd bacteria into the community. We found that the two treefrog species had a similar bacterial community structure, which differed from the more diverse community found on the terrestrial frog. Co-occurrence networks also revealed differences between frog species such that the treefrogs had a significantly higher number of culturable Bd-inhibitory OTUs with high centrality scores compared with the terrestrial frog. We found that culture-dependent OTUs captured between 21 and 39% of the total relative abundance revealed in culture-independent communities. Our results suggest different ecological strategies occurring within skin antifungal communities on host species that have not succumbed to Bd infections in the wild.

Genetic variability and ontogeny predict microbiome structure in a disease-challenged montane amphibian

Amphibian populations worldwide are at risk of extinction from infectious diseases, including chytridiomycosis caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd). Amphibian cutaneous microbiomes interact with Bd and can confer protective benefits to the host. The composition of the microbiome itself is influenced by many environment- and host-related factors. However, little is known about the interacting effects of host population structure, genetic variation and developmental stage on microbiome composition and Bd prevalence across multiple sites. Here we explore these questions in Amietia hymenopus, a disease-affected frog in southern Africa. We use microsatellite genotyping and 16S amplicon sequencing to show that the microbiome associated with tadpole mouthparts is structured spatially, and is influenced by host genotype and developmental stage. We observed strong genetic structure in host populations based on rivers and geographic distances, but this did not correspond to spatial patterns in microbiome composition. These results indicate that demographic and host genetic factors affect microbiome composition within sites, but different factors are responsible for host population structure and microbiome structure at the between-site level. Our results help to elucidate complex within- and among- population drivers of microbiome structure in amphibian populations. That there is a genetic basis to microbiome composition in amphibians could help to inform amphibian conservation efforts against infectious diseases.

Effects of Emerging Infectious Diseases on Amphibians: A Review of Experimental Studies

Numerous factors are contributing to the loss of biodiversity. These include complex effects of multiple abiotic and biotic stressors that may drive population losses. These losses are especially illustrated by amphibians, whose populations are declining worldwide. The causes of amphibian population declines are multifaceted and context-dependent. One major factor affecting amphibian populations is emerging infectious disease. Several pathogens and their associated diseases are especially significant contributors to amphibian population declines. These include the fungi Batrachochytrium dendrobatidis and B. salamandrivorans, and ranaviruses. In this review, we assess the effects of these three pathogens on amphibian hosts as found through experimental studies. Such studies offer valuable insights to the causal factors underpinning broad patterns reported through observational studies. We summarize key findings from experimental studies in the laboratory, in mesocosms, and from the field. We also summarize experiments that explore the interactive effects of these pathogens with other contributors of amphibian population declines. Though well-designed experimental studies are critical for understanding the impacts of disease, inconsistencies in experimental methodologies limit our ability to form comparisons and conclusions. Studies of the three pathogens we focus on show that host susceptibility varies with such factors as species, host age, life history stage, population and biotic (e.g., presence of competitors, predators) and abiotic conditions (e.g., temperature, presence of contaminants), as well as the strain and dose of the pathogen, to which hosts are exposed. Our findings suggest the importance of implementing standard protocols and reporting for experimental studies of amphibian disease.

Increased Numbers of Culturable Inhibitory Bacterial Taxa May Mitigate the Effects of Batrachochytrium dendrobatidis in Australian Wet Tropics Frogs

Symbiotic bacterial communities resident on amphibian skin can benefit their hosts. For example, antibiotic production by community members can control the pathogen Batrachochytrium dendrobatidis (Bd) and it is possible for these community members to be used as probiotics to reduce infection levels. In the early 1990s, the emergence of Bd caused declines and disappearances of frogs in the Australian Wet Tropics; the severity of its effects varied among species and sites. Some species have since recolonized despite enzootic Bd within their populations. This variation in history among species and sites provided an opportunity to investigate the role of anti-fungal cutaneous bacteria in protecting frogs against Bd infection. We collected cutaneous swab samples from three species of frogs at two upland and two lowland sites in the Wet Tropics, and used in vitro challenge assays to identify culturable Bd-inhibitory bacterial isolates for further analysis. We sequenced DNA from cultured inhibitory isolates to identify taxa, resulting in the classification of 16 Bd-inhibitory OTUs, and determined whether inhibitory taxa were associated with frog species, site, or intensity of infection. We present preliminary results showing that the upper limit of Bd infection intensity was negatively correlated with number of inhibitory OTUs present per frog indicating that increased numbers of Bd-inhibiting taxa may play a role in reducing the intensity of Bd infections, facilitating frog coexistence with enzootic Bd. One upland site had a significantly lower prevalence of Bd infection, a significantly higher proportion of frogs with one or more culturable Bd-inhibitory OTUs, a greater number of inhibitory bacterial genera present per frog, and statistically significant clustering of individual frogs with similar Bd-inhibitory signatures when compared to all other sites. This suggests that Bd-inhibitory taxa are likely to be particularly important to frogs at this site and may have played a role in their ability to recolonize following population declines. Our findings suggest that the use of multi-taxon Bd-inhibitory probiotics to support at-risk amphibian populations may be more effective than single-taxon alternatives.

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

Symbiotic microbial communities play key roles in the health and development of their multicellular hosts. Understanding why microbial communities vary among different host species or individuals is an important step toward understanding the diversity and function of the microbiome. The amphibian skin microbiome may affect resistance to the fungal pathogen Batrachochytrium dendrobatidis (Bd). Still, the factors that determine the diversity and composition of the amphibian skin microbiome, and therefore may ultimately contribute to disease resistance, are not well understood. We conducted a two-phase experiment to first test how host and environment shape the amphibian skin microbiome, and then test if the microbiome affects or is affected by Bd infection. Most lab experiments testing assembly of the amphibian skin microbiome so far have compared sterile to non-sterile environments or heavily augmented to non-augmented frogs. A goal of this study was to evaluate, in an experimental setting, realistic potential drivers of microbiome assembly that would be relevant to patterns observed in nature. We tested effects of frog genetic background (2 source populations) and 6 natural lake water sources in shaping the microbiome of the frog Rana sierrae. Water in which frogs were housed affected the microbiome in a manner that partially mimicked patterns observed in natural populations. In particular, frogs housed in water from disease-resistant populations had greater bacterial richness than frogs housed in water from populations that died out due to Bd. However, in the experiment this difference in microbiomes did not lead to differences in host mortality or rates of pathogen load increase. Frog source population also affected the microbiome and, although none of the frogs in this study showed true resistance to infection, host source population had a small effect on the rate of pathogen load increase. This difference in infection trajectories could be due to the observed differences in the microbiome, but could also be due to other traits that differ between frogs from the two populations. In addition to examining effects of the microbiome on Bd, we tested the effect of Bd infection severity on the microbiome. Specifically, we studied a time series of the microbiome over the course of infection to test if the effects of Bd on the microbiome are dependent on Bd infection severity. Although limited to a small subset of frogs, time series analysis suggested that relative abundances of several bacterial phylotypes changed as Bd loads increased through time, indicating that Bd-induced disturbance of the R. sierrae microbiome is not a binary effect but instead is dependent on infection severity. We conclude that both host and aquatic environment help shape the R. sierrae skin microbiome, with links to small changes in disease resistance in some cases, but in this study the effect of Bd on the microbiome was greater than the effect of the microbiome on Bd. Assessment of the microbiome differences between more distantly related populations than those studied here is needed to fully understand the role of the microbiome in resistance to Bd.

The Skin Microbiome of the Neotropical Frog Craugastor fitzingeri: Inferring Potential Bacterial-Host-Pathogen Interactions From Metagenomic Data

Skin symbiotic bacteria on amphibians can play a role in protecting their host against pathogens. Chytridiomycosis, the disease caused by Batrachochytrium dendrobatidis, Bd, has caused dramatic population declines and extinctions of amphibians worldwide. Anti-Bd bacteria from amphibian skin have been cultured, and skin bacterial communities have been described through 16S rRNA gene amplicon sequencing. Here, we present a shotgun metagenomic analysis of skin bacterial communities from a Neotropical frog, Craugastor fitzingeri. We sequenced the metagenome of six frogs from two different sites in Panamá: three frogs from Soberanía (Sob), a Bd-endemic site, and three frogs from Serranía del Sapo (Sapo), a Bd-naïve site. We described the taxonomic composition of skin microbiomes and found that Pseudomonas was a major component of these communities. We also identified that Sob communities were enriched in Actinobacteria while Sapo communities were enriched in Gammaproteobacteria. We described gene abundances within the main functional classes and found genes enriched either in Sapo or Sob. We then focused our study on five functional classes of genes: biosynthesis of secondary metabolites, metabolism of terpenoids and polyketides, membrane transport, cellular communication and antimicrobial drug resistance. These gene classes are potentially involved in bacterial communication, bacterial-host and bacterial-pathogen interactions among other functions. We found that C. fitzingeri metagenomes have a wide array of genes that code for secondary metabolites, including antibiotics and bacterial toxins, which may be involved in bacterial communication, but could also have a defensive role against pathogens. Several genes involved in bacterial communication and bacterial-host interactions, such as biofilm formation and bacterial secretion systems were found. We identified specific genes and pathways enriched at the different sites and determined that gene co-occurrence networks differed between sites. Our results suggest that skin microbiomes are composed of distinct bacterial taxa with a wide range of metabolic capabilities involved in bacterial defense and communication. Differences in taxonomic composition and pathway enrichments suggest that skin microbiomes from different sites have unique functional properties. This study strongly supports the need for shotgun metagenomic analyses to describe the functional capacities of skin microbiomes and to tease apart their role in host defense against pathogens.

Skin Microbiomes of California Terrestrial Salamanders Are Influenced by Habitat More Than Host Phylogeny

A multitude of microorganisms live on and within plant and animal hosts, yet the ecology and evolution of these microbial communities remains poorly understood in many taxa. This study examined the extent to which environmental factors and host taxonomic identity explain microbiome variation within two salamander genera, Ensatina and Batrachoseps, in the family Plethodontidae. In particular, we assessed whether microbiome differentiation paralleled host genetic distance at three levels of taxonomy: genus and high and low clade levels within Ensatina eschscholtzii. We predicted that more genetically related host populations would have more similar microbiomes than more distantly related host populations. We found that salamander microbiomes possess bacterial species that are most likely acquired from their surrounding soil environment, but the relative representation of those bacterial species is significantly different on the skin of salamanders compared to soil. We found differences in skin microbiome alpha diversity among Ensatina higher and lower clade groups, as well as differences between Ensatina and Batrachoseps. We also found that relative microbiome composition (beta diversity) did vary between Ensatina lower clades, but differences were driven by only a few clades and not correlated to clade genetic distances. We conclude this difference was likely a result of Ensatina lower clades being associated with geographic location and habitat type, as salamander identity at higher taxonomic levels (genus and Ensatina higher clades) was a weak predictor of microbiome composition. These results lead us to conclude that environmental factors are likely playing a more significant role in salamander cutaneous microbiome assemblages than host-specific traits.

Influence of immunogenetics, sex and body condition on the cutaneous microbial communities of two giant salamanders

The complex association between hosts and microbial symbionts requires the implementation of multiple approaches to evaluate variation in host physiology. Within amphibians, heterogeneity in immunogenetic traits and cutaneous microbiota is associated with variation in disease resistance. Ozark (Cryptobranchus alleganiensis bishopi) and eastern hellbenders (C. a. alleganiensis) provide a model system to assess variation in host traits and microbial communities. Ozark hellbenders have experienced declines throughout their range, are federally endangered and experience wound retardation that is absent in the eastern subspecies. Previous microbial investigations indicate differentiation in the composition of the skin microbiota of both hellbender subspecies, but it is not clear whether these patterns are concurrent with diversity in the major histocompatibility complex (MHC) genes. We characterized the MHC IIB and the skin microbiota of hellbenders in Missouri, where both subspecies co-occur though not sympatric. We compared the microbiota composition and MHC diversity between both subspecies and investigated whether individual-level MHC diversity, sex and body condition were associated with microbiota composition. Overall, MHC IIB diversity was lower in Ozark hellbenders compared to the eastern subspecies. Multivariate statistical comparisons identified microbiota differentiation between Ozark and eastern hellbenders. MHC IIB allele presence/absence, allele divergence, body composition and sex defined grouping of hellbender microbiotas within populations. Differentiation of the cutaneous microbiotas and MHC IIB genes between eastern and Ozark hellbenders suggests that differences exist in immunity between the two subspecies. This study demonstrates how simultaneous assessments of host genetic traits and microbiotas can inform patterns of microbial community structure in natural systems.