Search

Global biodiversity loss and mass extinction of species are two of the most critical environmental issues the world is currently facing, resulting in the disruption of various ecosystems central to environmental functions and human health. Microbiome-targeted interventions, such as probiotics and microbiome transplants, are emerging as potential options to reverse deterioration of biodiversity and increase the resilience of wildlife and ecosystems. However, the implementation of these interventions is urgently needed. We summarize the current concepts, bottlenecks and ethical aspects encompassing the careful and responsible management of ecosystem resources using the microbiome (termed microbiome stewardship) to rehabilitate organisms and ecosystem functions. We propose a real-world application framework to guide environmental and wildlife probiotic applications. This framework details steps that must be taken in the upscaling process while weighing risks against the high toll of inaction. In doing so, we draw parallels with other aspects of contemporary science moving swiftly in the face of urgent global challenges. published

Lactobacillaceae are an important family of lactic acid bacteria that play key roles in the gut microbiome of many animal species. In the honey bee (Apis mellifera) gut microbiome, many species of Lactobacillaceae are found, and there is functionally important strain-level variation in the bacteria. In this study, we completed whole-genome sequencing of 3 unique Lactobacillaceae isolates collected from hives in Virginia, USA. Using 107 genomes of known bee-associated Lactobacillaceae and Limosilactobacillus reuteri as an outgroup, the phylogenetics of the 3 isolates was assessed, and these isolates were identified as novel strains of Apilactobacillus kunkeei, Lactobacillus kullabergensis, and Bombilactobacillus mellis. Genome rearrangements, conserved orthologous genes (COG) categories and potential prophage regions were identified across the 3 novel strains. The new A. kunkeei strain was enriched in genes related to replication, recombination and repair, the L. kullabergensis strain was enriched for carbohydrate transport, and the B. mellis strain was enriched in transcription or transcriptional regulation and in some genes with unknown functions. Prophage regions were identified in the A. kunkeei and L. kullabergensis isolates. These new bee-associated strains add to our growing knowledge of the honey bee gut microbiome, and to Lactobacillaceae genomics more broadly.

Laboratory experiments have advanced our understanding of honey bee (Apis mellifera) responses to environmental factors, but removal from the hive environment may also impact physiology. To examine whether the laboratory environment alters the honey bee gut bacterial community and immune responses, we compared bacterial community structure (based on amplicon sequence variant relative abundance), total bacterial abundance, and immune enzyme (phenoloxidase and glucose oxidase) activity of cohort honey bee workers kept under laboratory and hive conditions. Workers housed in the laboratory showed differences in the relative abundance of their core gut taxa, an increase in total gut bacterial abundance, and reduced phenoloxidase activity, compared to bees housed in hives. Virginia Agricultural Foundation [677] Published version We would like to thank Philene Vu for performing the immunological assays, and Hannah Miko for assisting with lab and field work. This work was supported by the Virginia Agricultural Foundation (Project number 677). We thank Virginia Tech's Open Access Subvention Fund for assistance with open access fees.

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.

Animal-associated microbiomes are integral to host health, yet key biotic and abiotic factors that shape host-associated microbial communities at the global scale remain poorly understood. We investigated global patterns in amphibian skin bacterial communities, incorporating samples from 2,349 individuals representing 205 amphibian species across a broad biogeographic range. We analysed how biotic and abiotic factors correlate with skin microbial communities using multiple statistical approaches. Global amphibian skin bacterial richness was consistently correlated with temperature-associated factors. We found more diverse skin microbiomes in environments with colder winters and less stable thermal conditions compared with environments with warm winters and less annual temperature variation. We used bioinformatically predicted bacterial growth rates, dormancy genes and antibiotic synthesis genes, as well as inferred bacterial thermal growth optima to propose mechanistic hypotheses that may explain the observed patterns. We conclude that temporal and spatial characteristics of the host’s macro-environment mediate microbial diversity. National Science Foundation/[DEB-1146284]/NSF/Estados Unidos National Science Foundation/[IOS-1121758]/NSF/Estados Unidos National Science Foundation/[DEB-1310036]/NSF/Estados Unidos John Templeton Foundation/[]/JTF/Estados Unidos Deutsche Forschungsgemeinschaft/[]/DFG/Alemania Deutsche Forschungsgemeinschaft/[VE247/9-1]/DFG/Alemania Coordenação de Aperfeiçoamento de Pessoal de Nível Superior/[]/CAPES/Brasil Fundação de Amparo à Pesquisa do Estado de São Paulo/[#2013/50741-7]/FAPESP/Brasil Conselho Nacional de Desenvolvimento Científico e Tecnológico/[]/CNPq/Brasil Simons Foundation/[429440, WTW]//Estados Unidos Deutscher Akademischer Austauschdienst/[]/DAAD/Alemania Universidad de Costa Rica/[801-B2-029]/UCR/Costa Rica Ministerio de Ciencia, Tecnología y Telecomunicaciones/[849-PINN-2015]/MICITT/Costa Rica National Research Foundation of Korea/[2015R1D1A1A01057282]/NRF/Corea del Sur UCR::Vicerrectoría de Investigación::Unidades de Investigación::Ciencias Básicas::Centro de Investigación en Biología Celular y Molecular (CIBCM) UCR::Vicerrectoría de Investigación::Unidades de Investigación::Ciencias Básicas::Centro de Investigación en Estructuras Microscópicas (CIEMIC)

Interactions between pathogens and the symbiotic microbial communities that reside in hosts can impact disease processes. Symbiotic microbial communities can act to prevent pathogen infection in some instances, while in other instances, pathogens can disrupt these symbiotic communities. We sought to address these interactions between the fungal skin pathogen Batrachochytrium dendrobatidis (Bd) and the skin bacterial communities of adult Spring Peepers (Pseudacris crucifer). In the laboratory, frogs were exposed to Bd zoospores for a 12-hour period and then monitored for five days. Bacterial community composition on the skin was assessed upon initial collection, 24 hours post-exposure, and at the conclusion of the experiment using 16S rRNA gene amplicon sequencing. Little impact on the bacterial community was observed as a result of Bd exposure, suggesting that the timeline was either too short to observe the effects of Bd, or that the skin communities are resilient in the face of such short exposure periods. In contrast, laboratory housing was associated with changes in skin bacterial community composition, in terms of both relative abundances and bacterial taxa that were present. These findings suggest that even a short time away from their natural environment can have strong effects on the composition of amphibian skin communities and that Bd exposure may not always disrupt these skin communities.

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.

Host-associated microbial communities can influence the overall health of their animal hosts, and many factors, including behavior and physiology, can impact the formation of these complex communities. Bacteria within these communities can be transmitted socially between individuals via indirect (e.g., shared environments) or direct (e.g., physical contact) pathways. Limited research has been done to investigate how social interactions that occur in the context of mating shape host-associated microbial communities. To gain a better understanding of these interactions and, more specifically, to assess how mating behavior shapes an animal's microbiome, we studied the cloacal bacterial communities of a socially monogamous yet genetically polygynous songbird, the North American tree swallow (Tachycineta bicolor). We address two questions: (1) do the cloacal bacterial communities differ between female and male tree swallows within a population? and (2) do pair-bonded social partners exhibit more similar cloacal bacterial communities than expected by chance? To answer these questions, we sampled the cloacal microbiome of adults during the breeding season and then used culture-independent, 16S rRNA gene amplicon sequencing to assess bacterial communities. Overall, we found that the cloacal bacterial communities of females and males were similar, and that the communities of pair-bonded social partners were not more similar than expected by chance. Our results suggest that social monogamy does not correlate with an increased similarity in cloacal bacterial community diversity or structure. As social partners were not assessed at the same time, it is possible that breeding stage differences masked social effects on bacterial community diversity and structure. Further, given that tree swallows exhibit high variation in rates of extra-pair activity, considering extra-pair activity when assessing cloacal microbial communities may be important for understanding how these bacterial communities are shaped. Further insight into how bacterial communities are shaped will ultimately shed light on potential tradeoffs associated with alternative behavioral strategies and socially-transmitted microbes.

Summary Testosterone mediates several key aspects of male reproduction, but maintaining high testosterone levels can reduce long-term survival. One of the proposed pathways by which testosterone can influence survival is through increased risk of parasite infection. The hormone has the potential to affect the transmission of sexually transmitted infections by promoting behaviours that increase sexual contact rates and/or by decreasing immune function. We hypothesized that males with high concentrations of testosterone would increase their chances of being infected with sexually transmitted bacteria, which would manifest as increased diversity of cloacal bacteria. To test this hypothesis, we measured circulating testosterone concentrations in breeding male rufous-collared sparrows (Zonotrichia capensis) and collected cloacal swabs to quantify bacterial diversity using 16S rRNA gene amplicon sequencing. There was a positive correlation between testosterone concentrations and the phylogenetic diversity of cloacal bacteria. In addition, individuals with high and medium testosterone concentrations had cloacal bacterial communities that were more similar to each other than to those of low testosterone individuals. Finally, when considering bacterial taxa that are potential avian pathogens, we found that the relative abundance of Chlamydiae, a class of obligate intracellular parasites, was positively correlated with testosterone concentrations. Two non-exclusive explanations for these results are that testosterone affects behaviours that lead to increased sexual contacts and thus the exposure and acquisition of additional phylogenetically diverse bacteria, and/or that testosterone is altering the immune system or the cloacal environment, thus making it easier for bacteria to colonize. Either way, these data suggest that increased exposure to sexually transmitted pathogens in the form of cloacal bacteria could be a cost of maintaining high testosterone levels. A Lay Summary is available for this article.

Amphibian population declines caused by the fungus Batrachochytrium dendrobatidis (Bd) have prompted studies on the bacterial community that resides on amphibian skin. However, studies addressing the fungal portion of these symbiont communities have lagged behind. Using ITS1 amplicon sequencing, we examined the fungal portion of the skin microbiome of temperate and tropical amphibian species currently coexisting with Bd in nature. We assessed cooccurrence patterns between bacterial and fungal OTUs using a subset of samples for which bacterial 16S rRNA gene amplicon data were also available. We determined that fungal communities were dominated by members of the phyla Ascomycota and Basidiomycota, and also by Chytridiomycota in the most aquatic amphibian species. Alpha diversity of the fungal communities differed across host species, and fungal community structure differed across species and regions. However, we did not find a correlation between fungal diversity/community structure and Bd infection, though we did identify significant correlations between Bd and specific OTUs. Moreover, positive bacterial-fungal cooccurrences suggest that positive interactions between these organisms occur in the skin microbiome. Understanding the ecology of amphibian skin fungi, and their interactions with bacteria will complement our knowledge of the factors influencing community assembly and the overall function of these symbiont communities.

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.

Currently, it is estimated that only 0.001% to 15% of bacteria in any given system can be cultured by use of commonly used techniques and media, yet culturing is critically important for investigations of bacterial function. Despite this situation, few studies have attempted to link culture-dependent and culture-independent data for a single system to better understand which members of the microbial community are readily cultured. In amphibians, some cutaneous bacterial symbionts can inhibit establishment and growth of the fungal pathogen Batrachochytrium dendrobatidis , and thus there is great interest in using these symbionts as probiotics for the conservation of amphibians threatened by B. dendrobatidis . The present study examined the portion of the culture-independent bacterial community (based on Illumina amplicon sequencing of the 16S rRNA gene) that was cultured with R2A low-nutrient agar and whether the cultured bacteria represented rare or dominant members of the community in the following four amphibian species: bullfrogs ( Lithobates catesbeianus ), eastern newts ( Notophthalmus viridescens ), spring peepers ( Pseudacris crucifer ), and American toads ( Anaxyrus americanus ). To determine which percentage of the community was cultured, we clustered Illumina sequences at 97% similarity, using the culture sequences as a reference database. For each amphibian species, we cultured, on average, 0.59% to 1.12% of each individual's bacterial community. However, the average percentage of bacteria that were culturable for each amphibian species was higher, with averages ranging from 2.81% to 7.47%. Furthermore, most of the dominant operational taxonomic units (OTUs), families, and phyla were represented in our cultures. These results open up new research avenues for understanding the functional roles of these dominant bacteria in host health.

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.

One current challenge in microbial ecology is elucidating the functional roles of the large diversity of free-living and host-associated bacteria identified by culture-independent molecular methods. Importantly, the characterization of this immense bacterial diversity will likely require merging data from culture-independent approaches with work on bacterial isolates in culture. Amphibian skin bacterial communities have become a recent focus of work in host-associated microbial systems due to the potential role of these skin bacteria in host defense against the pathogenic fungus Batrachochytrium dendrobatidis (Bd), which is associated with global amphibian population declines and extinctions. As there is evidence that some skin bacteria may inhibit growth of Bd and prevent infection in some cases, there is interest in using these bacteria as probiotic therapy for conservation of at-risk amphibians. In this study, we used skin swabs from American toads (Anaxyrus americanus) to: (1) assess the diversity and community structure of culturable amphibian skin bacteria grown on high and low nutrient culture media, (2) determine which culture media recover the highest proportion of the total skin bacterial community of individual toads relative to culture-independent data, and (3) assess whether the plated communities from the distinct media types vary in their ability to inhibit Bd growth in in-vitro assays. Overall, we found that culture media with low nutrient concentrations facilitated the growth of more diverse bacterial taxa and grew distinct communities relative to media with higher nutrient concentrations. Use of low nutrient media also resulted in culturing proportionally more of the bacterial diversity on individual toads relative to the overall community defined using culture-independent methods. However, while there were differences in diversity among media types, the variation among individual hosts was greater than variation among media types, suggesting that swabbing more individuals in a population is the best way to maximize culture collections, regardless of media type. Lastly, the function of the plated communities against Bd did not vary across culture media type or between high and low nutrient media. These results inform current efforts for developing a probiotic-based approach for amphibian conservation and help to ensure that culture collections are capturing the majority of the important diversity in these systems.

Summary Some amphibian skin bacteria inhibit growth of a fungal amphibian pathogen, Batrachochytrium dendrobatidis (Bd), but it is unclear how dominant these anti-Bd bacteria are in skin communities. Using in vitro co-culture challenge assays, we quantified Bd inhibition by bacterial isolates collected from the skin of four amphibian species: bullfrogs, Eastern newts, spring peepers, and American toads. The 16S rRNA sequences for each isolate were matched to culture-independent amplicon sequences from the same individuals to assess inhibitory function versus relative abundance. Dominant bacteria had higher Bd inhibition than rare bacteria in bullfrog and newt populations, in which Bd was prevalent (>25%). Dominant and rare bacteria did not differ in Bd inhibition in spring peeper and toad populations, in which Bd was absent or at low prevalence (

Host-microbe symbioses rely on the successful transmission or acquisition of symbionts in each new generation. Amphibians host a diverse cutaneous microbiota, and many of these symbionts appear to be mutualistic and may limit infection by the chytrid fungus, Batrachochytrium dendrobatidis, which has caused global amphibian population declines and extinctions in recent decades. Using bar-coded 454 pyrosequencing of the 16S rRNA gene, we addressed the question of symbiont transmission by examining variation in amphibian skin microbiota across species and sites and in direct relation to environmental microbes. Although acquisition of environmental microbes occurs in some host-symbiont systems, this has not been extensively examined in free-living vertebrate-microbe symbioses. Juvenile bullfrogs (Rana catesbeiana), adult red-spotted newts (Notophthalmus viridescens), pond water and pond substrate were sampled at a single pond to examine host-specificity and potential environmental transmission of microbiota. To assess population level variation in skin microbiota, adult newts from two additional sites were also sampled. Cohabiting bullfrogs and newts had distinct microbial communities, as did newts across the three sites. The microbial communities of amphibians and the environment were distinct; there was very little overlap in the amphibians' core microbes and the most abundant environmental microbes, and the relative abundances of OTUs that were shared by amphibians and the environment were inversely related. These results suggest that, in a host species-specific manner, amphibian skin may select for microbes that are generally in low abundance in the environment.

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 Panama. 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. National Science Foundation [DEB-1136602, DEB-1136640] This project was supported by the National Science Foundation Grants: DEB-1136602 (to Reid N. Harris) and DEB-1136640 (Lisa K. Belden). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

All multicellular organisms are host to microbial symbionts that constitute the microbiome and can have significant impacts on the host, including altering development, behavior, and health [1]. In turn, aspects of the host and their environment can influence the microbiome [2]. Here, we briefly summarize current knowledge of the amphibian skin microbiome and its role in heath and disease. Given the increase in fungal diseases that now threaten amphibians and other wildlife—including bees, bats, snakes, and corals, as well as a variety of economically important crops [3]—we hope that lessons learned from amphibian host–microbe interactions can also ultimately be applied in other systems (Fig 1). Funding was provided by the National Science Foundation (DEB-1136640). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Both the structure and function of host-associated microbial communities are potentially impacted by environmental conditions, just as the outcomes of many free-living species interactions are context-dependent. Many amphibian populations have declined around the globe due to the fungal skin pathogen, Batrachochytrium dendrobatidis (Bd), but enivronmental conditions may influence disease dynamics. For instance, in Panama, the most severe Bd outbreaks have occurred at high elevation sites. Some amphibian species harbor bacterial skin communities that can inhibit the growth of Bd, and therefore, there is interest in understanding whether environmental context could also alter these host-associated microbial communities in a way that might ultimately impact Bd dynamics. In a field survey in Panama, we assessed skin bacterial communities (16S rRNA amplicon sequencing) and metabolite profiles (HPLC-UV/Vis) of Silverstoneia flotator from three high- and three low-elevation populations representing a range of environmental conditions. Across elevations, frogs had similar skin bacterial communities, although one lowland site appeared to differ. Interestingly, we found that bacterial richness decreased from west to east, coinciding with the direction of Bd spread through Panama. Moreover, metabolite profiles suggested potential functional variation among frog populations and between elevations. While the frogs have similar bacterial community structure, the local environment might shape the metabolite profiles. Ultimately, host-associated community structure and function could be dependent on environmental conditions, which could ultimately influence host disease susceptibility across sites.

Disruptions to the microbiome can impact host health as can exposure to environmental contaminants. However, few studies have addressed how environmental contaminants impact the microbiome. We explored this question for frogs that breed in wetlands contaminated with fly ash, a by-product of coal combustion that is enriched in trace elements. We found differences in the bacterial communities among a fly ash-contaminated site and several reference wetlands. We then experimentally assessed the impacts of fly ash on the skin microbiome of adult spring peepers ( Pseudacris crucifer ). Frogs were exposed to fly ash in the laboratory for 12 h, the duration of a typical breeding event, and the skin microbiome was assessed after 5 days (experiment 1) or after 5 and 15 days (experiment 2). We examined bacterial community structure using 16S rRNA gene amplicon sequencing and metabolite profiles using high-pressure liquid chromatography-mass spectrometry (HPLC-MS). We found little impact as the result of acute exposure to fly ash on the bacterial communities or metabolite profiles in either experiment, suggesting that the bacterial symbiont communities of adults may be relatively resistant to brief contaminant exposure. However, housing frogs in the laboratory altered bacterial community structure in the two experiments, which supports prior research suggesting that environmental source pools are important for maintaining the amphibian skin microbiome. Therefore, for contaminants like fly ash that may alter the potential source pool of symbionts, we think it may be important to explore how contaminants affect the initial assembly of the amphibian skin microbiome in larval amphibians that develop within contaminated sites. IMPORTANCE Animals are hosts to many symbiotic microorganisms, collectively called the microbiome, that play critical roles in host health. Therefore, environmental contaminants that alter the microbiome may impact hosts. Some of the most widespread contaminants, produced worldwide, are derived from the mining, storage, and combustion of coal for energy. Fly ash, for example, is a by-product of coal combustion. It contains compounds such as arsenic, selenium, cadmium, and strontium and is a recognized source of ground and surface water contamination. Here, we experimentally investigated the impacts of short-term fly ash exposure on the skin microbiome of spring peepers, one of many species of amphibian that sometimes breed in open fly ash disposal ponds. This research provides a look into the potential impacts of fly ash on an animal's microbiome and suggests important future directions for research on the effects of environmental contaminants on the microbiome.

Emerging infectious diseases in wildlife are responsible for massive population declines. In amphibians, chytridiomycosis caused by Batrachochytrium dendrobatidis, Bd, has severely affected many amphibian populations and species around the world. One promising management strategy is probiotic bioaugmentation of antifungal bacteria on amphibian skin. In vivo experimental trials using bioaugmentation strategies have had mixed results, and therefore a more informed strategy is needed to select successful probiotic candidates. Metagenomic, transcriptomic, and metabolomic methods, colloquially called "omics," are approaches that can better inform probiotic selection and optimize selection protocols. The integration of multiple omic data using bioinformatic and statistical tools and in silico models that link bacterial community structure with bacterial defensive function can allow the identification of species involved in pathogen inhibition. We recommend using 16S rRNA gene amplicon sequencing and methods such as indicator species analysis, the Kolmogorov-Smirnov Measure, and co-occurrence networks to identify bacteria that are associated with pathogen resistance in field surveys and experimental trials. In addition to 16S amplicon sequencing, we recommend approaches that give insight into symbiont function such as shotgun metagenomics, metatranscriptomics, or metabolomics to maximize the probability of finding effective probiotic candidates, which can then be isolated in culture and tested in persistence and clinical trials. An effective mitigation strategy to ameliorate chytridiomycosis and other emerging infectious diseases is necessary; the advancement of omic methods and the integration of multiple omic data provide a promising avenue toward conservation of imperiled species.

Amphibian embryos are at risk of microbial infection. Here we find evidence that innate immune defenses, both antimicrobial skin peptides and mutualistic microbiota, of adult glass frogs, Hyalinobatrachium colymbiphyllum, can be transmitted to embryos deposited on leaves above rain forest streams in Panama and can inhibit the fungal pathogen Batrachochytrium dendrobatidis. Abstract in Spanish is available at http://www.blackwell-synergy.com/loi/btp.

Vertebrates, including amphibians, host diverse symbiotic microbes that contribute to host disease resistance. Globally, and especially in montane tropical systems, many amphibian species are threatened by a chytrid fungus, Batrachochytrium dendrobatidis (Bd), that causes a lethal skin disease. Bd therefore may be a strong selective agent on the diversity and function of the microbial communities inhabiting amphibian skin. In Panama, amphibian population declines and the spread of Bd have been tracked. In 2012, we completed a field survey in Panama to examine frog skin microbiota in the context of Bd infection. We focused on three frog species and collected two skin swabs per frog from a total of 136 frogs across four sites that varied from west to east in the time since Bd arrival. One swab was used to assess bacterial community structure using 16S rRNA amplicon sequencing and to determine Bd infection status, and one was used to assess metabolite diversity, as the bacterial production of anti fungal metabolites is an important disease resistance function. The skin microbiota of the three Panamanian frog species differed in OTU (operational taxonomic unit, bacterial species) community composition and metabolite profiles, although the pattern was less strong for the metabolites. Comparisons between frog skin bacterial communities from Panama and the US suggest broad similarities at the phylum level, but key differences at lower taxonomic levels. In our field survey in Panama, across all four sites, only 35 individuals (similar to 26%) were Bd infected. There was no clustering of OTUs or metabolite profiles based on Bd infection status and no clear pattern of west east changes in OTUs or metabolite profiles across the four sites. Overall, our field survey data suggest that different bacterial communities might be producing broadly similar sets of metabolites across frog hosts and sites. Community structure and function may not be as tightly coupled in these skin symbiont microbial systems as it is in many macro systems.

Symbiotic microbes can dramatically impact host health and fitness, and recent research in a diversity of systems suggests that different symbiont community structures may result in distinct outcomes for the host. In amphibians, some symbiotic skin bacteria produce metabolites that inhibit the growth of Batrachochytrium dendrobatidis (Bd), a cutaneous fungal pathogen that has caused many amphibian population declines and extinctions. Treatment with beneficial bacteria (probiotics) prevents Bd infection in some amphibian species and creates optimism for conservation of species that are highly susceptible to chytridiomycosis, the disease caused by Bd. In a laboratory experiment, we used Bd-inhibitory bacteria from Bd-tolerant Panamanian amphibians in a probiotic development trial with Panamanian golden frogs, Atelopus zeteki , a species currently surviving only in captive assurance colonies. Approximately 30% of infected golden frogs survived Bd exposure by either clearing infection or maintaining low Bd loads, but this was not associated with probiotic treatment. Survival was instead related to initial composition of the skin bacterial community and metabolites present on the skin. These results suggest a strong link between the structure of these symbiotic microbial communities and amphibian host health in the face of Bd exposure and also suggest a new approach for developing amphibian probiotics.

The introduction of next-generation sequencing has allowed for greater understanding of community composition of symbiotic microbial communities. However, determining the function of individual members of these microbial communities still largely relies on culture-based methods. Here, we present results on the phylogenetic distribution of a defensive functional trait of cultured symbiotic bacteria associated with amphibians. Amphibians are host to a diverse community of cutaneous bacteria and some of these bacteria protect their host from the lethal fungal pathogen Batrachochytrium dendrobatidis (Bd) by secreting antifungal metabolites. We cultured over 450 bacterial isolates from the skins of Panamanian amphibian species and tested their interactions with Bd using an in vitro challenge assay. For a subset of isolates, we also completed coculture experiments and found that culturing isolates with Bd had no effect on inhibitory properties of the bacteria, but it significantly decreased metabolite secretion. In challenge assays, approximately 75% of the bacterial isolates inhibited Bd to some extent and these inhibitory isolates were widely distributed among all bacterial phyla. Although there was no clear phylogenetic signal of inhibition, three genera, Stenotrophomonas, Aeromonas and Pseudomonas, had a high proportion of inhibitory isolates (100%, 77% and 73%, respectively). Overall, our results demonstrate that antifungal properties are phylogenetically widespread in symbiotic microbial communities of Panamanian amphibians and that some functional redundancy for fungal inhibition occurs in these communities. We hope that these findings contribute to the discovery and development of probiotics for amphibians that can mitigate the threat of chytridiomycosis.

The emergence of the infectious disease chytridiomycosis is a major factor responsible for amphibian extinctions in pristine habitats. However, some populations coexist with Batrachochytrium dendrobatidis (Bd), the agent of chytridiomycosis, and others go extinct when Bd arrives. Variation of pathogen severity among populations may be explained by differences of antimicrobial skin peptides and anti-Bd skin bacteria. Previous work showed that a population of the frog Rana sierrae had a high proportion of individuals with at least one species of cultured anti-Bd skin bacteria and was coexisting with Bd for more than 6 years. A population of the closely related sister species Rana muscosa had a significantly lower proportion of individuals with anti-Bd bacteria, and the following year, it went extinct due to Bd. We extended previous work to include another sampling of the R. sierrae population coexisting with Bd and found that, although the anti-Bd bacterial community somewhat differed, both populations had a high proportion of individuals with antifungal bacteria. We also included a population of R. muscosa that was naive of Bd and predicted to go extinct once Bd emerged since that was the fate of neighboring populations. However, the naive R. muscosa population had a high proportion of antifungal bacteria, which suggested that the population might survive if Bd emerged. Two years later, the population was surviving with Bd endemic. Variation in the proportion of individuals with anti-Bd skin bacteria appears to explain why some populations persist when Bd emerges in a population.

Emerging infectious diseases threaten human and wildlife populations. Altered ecological interactions between mutualistic microbes and hosts can result in disease, but an understanding of interactions between host, microbes and disease-causing organisms may lead to management strategies to affect disease outcomes. Many amphibian species in relatively pristine habitats are experiencing dramatic population declines and extinctions due to the skin disease chytridiomycosis, which is caused by the chytrid fungus Batrachochytrium dendrobatidis. Using a randomized, replicated experiment, we show that adding an antifungal bacterial species, Janthinobacterium lividum, found on several species of amphibians to the skins of the frog Rana muscosa prevented morbidity and mortality caused by the pathogen. The bacterial species produces the anti-chytrid metabolite violacein, which was found in much higher concentrations on frog skins in the treatments where J. lividum was added. Our results show that cutaneous microbes are a part of amphibians' innate immune system, the microbial community structure on frog skins is a determinant of disease outcome and altering microbial interactions on frog skins can prevent a lethal disease outcome. A bioaugmentation strategy may be an effective management tool to control chytridiomycosis in amphibian survival assurance colonies and in nature.

The vertebrate microbiome contributes to disease resistance, but few experiments have examined the link between microbiome community structure and disease resistance functions. Chytridiomycosis, a major cause of amphibian population declines, is a skin disease caused by the fungus, Batrachochytrium dendrobatidis (Bd). In a factorial experiment, bullfrog skin microbiota was reduced with antibiotics, augmented with an anti-Bd bacterial isolate (Janthinobacterium lividum), or unmanipulated, and individuals were then either exposed or not exposed to Bd. We found that the microbial community structure of individual frogs prior to Bd exposure influenced Bd infection intensity one week following exposure, which, in turn, was negatively correlated with proportional growth during the experiment. Microbial community structure and function differed among unmanipulated, antibiotic-treated, and augmented frogs only when frogs were exposed to Bd. Bd is a selective force on microbial community structure and function, and beneficial states of microbial community structure may serve to limit the impacts of infection. This work was funded by the Virginia Tech Fralin Life Science Institute and the National Science Foundation (DEB-1136640 and CHE-0958973). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Published version

Microbial symbionts of vertebrate skin have an important function in defense of the host against pathogens. In particular, the emerging chytrid fungus Batrachochytrium dendrobatidis, causes widespread disease in amphibians but can be inhibited via secondary metabolites produced by many different skin-associated bacteria. Similarly, the fungal pathogens of terrestrial salamander eggs Mariannaea elegans and Rhizomucor variabilis are also inhibited by a variety of skin-associated bacteria. Indeed, probiotic therapy against fungal diseases is a recent approach in conservation medicine with growing experimental support. We present a comprehensive Antifungal Isolates Database of amphibian skin-associated bacteria that have been cultured, isolated, and tested for antifungal properties. At the start, this database includes nearly 2000 cultured bacterial isolates from 37 amphibian host species across 18 studies on five continents: Africa, Oceania, Europe, and North and South America. As the research community gathers information on additional isolates, the database will be updated periodically. The resulting database can serve as a conservation tool for amphibians and other organisms, and provides empirical data for comparative and bioinformatic studies. The database consists of a FASTA file containing 16S rRNA gene sequences of the bacterial isolates, and a metadata file containing information on the host species, life-stage, geographic region, and antifungal capacity and taxonomic identity of the isolate.