Bacteria from the skin of amphibians promote growth of Arabidopsis thaliana and Solanum lycopersicum by modifying hormone-related transcriptome response.

Plants and microorganisms establish beneficial associations that can improve their development and growth. Recently, it has been demonstrated that bacteria isolated from the skin of amphibians can contribute to plant growth and defense. However, the molecular mechanisms involved in the beneficial effect for the host are still unclear. In this work, we explored whether bacteria isolated from three tropical frogs species can contribute to plant growth. After a wide screening, we identified three bacterial strains with high biostimulant potential, capable of modifying the root structure of Arabidopsis thaliana plants. In addition, applying individual bacterial cultures to Solanum lycopersicum plants induced an increase in their growth. To understand the effect that these microorganisms have over the host plant, we analysed the transcriptomic profile of A. thaliana during the interaction with the C32I bacterium, demonstrating that the presence of the bacteria elicits a transcriptional response associated to plant hormone biosynthesis. Our results show that amphibian skin bacteria can function as biostimulants to improve agricultural crops growth and development by modifying the plant transcriptomic responses.

Host Species and Environment Shape the Skin Microbiota of Mexican Axolotls.

Skin microbiomes in amphibians are complex systems that can be influenced by biotic and abiotic factors. In this study, we examined the effect of host species and environmental conditions on the skin bacterial and fungal microbiota of four obligate paedomorphic salamander species, commonly known as axolotls (Ambystoma andersoni, A. dumerilii, A. mexicanum, and A. taylori), all of them endemic to the Trans-Mexican Volcanic Belt. We found that despite their permanent aquatic lifestyle, these species present a host-specific skin microbiota that is distinct from aquatic communities. We identified skin-associated taxa that were unique to each host species and that differentiated axolotl species based on alpha and beta diversity metrics. Moreover, we identified a set of microbial taxa that were shared across hosts with high relative abundances across skin samples. Specifically, bacterial communities were dominated by Burkholderiales and Pseudomonadales bacterial orders and Capnodiales and Pleosporales fungal orders. Host species and environmental variables collectively explained more microbial composition variation in bacteria (R2 = 0.46) in comparison to fungi (R2 = 0.2). Our results contribute to a better understanding of the factors shaping the diversity and composition of skin microbial communities in Ambystoma. Additional studies are needed to disentangle the effects of specific host associated and environmental factors that could influence the skin microbiome of these endangered species.

Metagenomic analysis of carbohydrate-active enzymes and their contribution to marine sediment biodiversity.

Marine sediments constitute the world's most substantial long-term carbon repository. The microorganisms dwelling in these sediments mediate the transformation of fixed oceanic carbon, but their contribution to the carbon cycle is not fully understood. Previous culture-independent investigations into sedimentary microorganisms have underscored the significance of carbohydrates in the carbon cycle. In this study, we employ a metagenomic methodology to investigate the distribution and abundance of carbohydrate-active enzymes (CAZymes) in 37 marine sediments sites. These sediments exhibit varying oxygen availability and were isolated in diverse regions worldwide. Our comparative analysis is based on the metabolic potential for oxygen utilisation, derived from genes present in both oxic and anoxic environments. We found that extracellular CAZyme modules targeting the degradation of plant and algal detritus, necromass, and host glycans were abundant across all metagenomic samples. The analysis of these results indicates that the oxic/anoxic conditions not only influence the taxonomic composition of the microbial communities, but also affect the occurrence of CAZyme modules involved in the transformation of necromass, algae and plant detritus. To gain insight into the sediment microbial taxa, we reconstructed metagenome assembled genomes (MAG) and examined the presence of primary extracellular carbohydrate active enzyme (CAZyme) modules. Our findings reveal that the primary CAZyme modules and the CAZyme gene clusters discovered in our metagenomes were prevalent in the Bacteroidia, Gammaproteobacteria, and Alphaproteobacteria classes. We compared those MAGs to organisms from the same taxonomic classes found in soil, and we found that they were similar in its CAZyme repertoire, but the soil MAG contained a more abundant and diverse CAZyme content. Furthermore, the data indicate that abundant classes in our metagenomic samples, namely Alphaproteobacteria, Bacteroidia and Gammaproteobacteria, play a pivotal role in carbohydrate transformation within the initial few metres of the sediments.

Genetic variation of Batrachochytrium dendrobatidis is linked to skin bacterial diversity in the Pacific treefrog Hyliola regilla (hypochondriaca).

Symbiotic bacterial communities are crucial to combating infections and contribute to host health. The amphibian skin microbiome plays an important role in protecting their hosts against pathogens such as Batrachochytrium dendrobatidis (Bd), one of the causative agents of chytridiomycosis, which is responsible for dramatic amphibian population declines worldwide. Although symbiotic skin bacteria are known to inhibit Bd growth, an understanding of the relationship between Bd genetic variability, environmental conditions, and skin bacterial communities is limited. Therefore, we examined the associations between Bd infection load, Bd genetic diversity and skin bacterial communities in five populations of Hyliola regilla (hypochondriaca) from environmentally contrasting sites in Baja California, Mexico. We observed differences in Bd genetics and infection load among sites and environments. Genetic analysis of Bd isolates revealed patterns of spatial structure corresponding to the five sites sampled. Amphibian skin bacterial diversity and community structure differed among environments and sites. Bacterial community composition was correlated with Bd genetic differences and infection load, with specific bacterial taxa enriched on infected and un-infected frogs. Our results indicate that skin-associated bacteria and Bd strains likely interact on the host skin, with consequences for microbial community structure and Bd infection intensity.

Comparative Analysis of Skin Bacterial Diversity and Its Potential Antifungal Function Between Desert and Pine Forest Populations of Boreal Toads Anaxyrus boreas.

The skin microbiome in amphibians has gained a lot of attention as some of its members play a protective role against pathogens such as the fungus Batrachochytrium dendrobatidis (Bd). The composition of skin bacterial communities has been suggested as one of the factors explaining differences in susceptibility to Bd among amphibian species and populations. The boreal toad Anaxyrus boreas is known to be susceptible to Bd, and severe population declines in its southeastern range have been documented. However, throughout A. boreas distribution, populations present differences in susceptibility to Bd infections which may be associated with differences in skin microbial diversity. This study compared the skin bacterial diversity and Bd infection levels of A. boreas in one desert population and one pine forest population from Baja California, Mexico. We found that desert and pine forest toad populations exhibit differences in skin bacterial community structure but show similar Bd infection levels. Using a predictive method, we found that the abundance of bacteria with potential Bd-inhibitory properties differed between uninfected and infected individuals but not between populations. Our data suggest that several bacteria in the skin community may be offering protection from Bd infections in these A. boreas populations. This study provides foundational evidence for future studies seeking to understand the skin-microbial variation among boreal toads' populations and its relation with Bd susceptibility.

Molecular detection of Bifidobacterium spp. in faeces of black howler monkeys (Alouatta pigra).

BACKGROUND: Bifidobacterium genus are considered to be beneficial bacteria for their hosts; however, knowledge about the specific species that are part of the gut microbiome of howler monkeys is scarce. Polymerase chain reaction (PCR) is a useful technique for the identification of non-cultivable or difficult to grow bacterial species. With the goal of detecting species of the genus Bifidobacterium in black howler monkeys, we used PCR on DNA derived from faecal samples. METHODS: We collected and extracted DNA from 40 faecal samples. Using specific primers, we performed PCR and nested PCR to detect members of the Bifidobacterium genus and a subset of species: Bifidobacterium adolescentis, Bifidobacterium bifidum, Bifidobacterium infantis, Bifidobacterium longum and Bifidobacterium animalis subsp. animalis. RESULTS: 97.5% (39/40) of the samples were positive for Bifidobacterium spp. We found B longum in 100% of the analysed samples. CONCLUSIONS: This is the first report of B longum in black howler monkey faeces.

Prevalence and pathogen load estimates for the fungus Batrachochytrium dendrobatidis are impacted by ITS DNA copy number variation.

The ribosomal gene complex is a multi-copy region that is widely used for phylogenetic analyses of organisms from all 3 domains of life. In fungi, the copy number of the internal transcribed spacer (ITS) is used to detect abundance of pathogens causing diseases such as chytridiomycosis in amphibians and white nose syndrome in bats. Chytridiomycosis is caused by the fungi Batrachochytrium dendrobatidis (Bd) and B. salamandrivorans (Bsal), and is responsible for declines and extinctions of amphibians worldwide. Over a decade ago, a qPCR assay was developed to determine Bd prevalence and pathogen load. Here, we demonstrate the effect that ITS copy number variation in Bd strains can have on the estimation of prevalence and pathogen load. We used data sets from different amphibian species to simulate how ITS copy number affects prevalence and pathogen load. In addition, we tested 2 methods (gBlocks® synthetic standards and digital PCR) to determine ITS copy number in Bd strains. Our results show that assumptions about the ITS copy number can lead to under- or overestimation of Bd prevalence and pathogen load. The use of synthetic standards replicated previously published estimates of ITS copy number, whereas dPCR resulted in estimates that were consistently lower than previously published estimates. Standardizing methods will assist with comparison across studies and produce reliable estimates of prevalence and pathogen load in the wild, while using the same Bd strain for exposure experiments and zoospore standards in qPCR remains the best method for estimating parameters used in epidemiological studies.

Direct and Indirect Horizontal Transmission of the Antifungal Probiotic Bacterium Janthinobacterium lividum on Green Frog (Lithobates clamitans) Tadpoles.

Amphibian populations worldwide are being threatened by the disease chytridiomycosis, which is caused by Batrachochytrium dendrobatidis To mitigate the effects of B. dendrobatidis, bioaugmentation of antifungal bacteria has been shown to be a promising strategy. One way to implement bioaugmentation is through indirect horizontal transmission, defined as the transfer of bacteria from a host to the environment and to another host. In addition, direct horizontal transmission among individuals can facilitate the spread of a probiotic in a population. In this study, we tested whether the antifungal bacterium Janthinobacterium lividum could be horizontally transferred, directly or indirectly, in a laboratory experiment using Lithobates clamitans tadpoles. We evaluated the ability of J. lividumto colonize the tadpoles' skin and to persist through time using culture-dependent and culture-independent techniques. We also tested whether the addition of J. lividum affected the skin community in L. clamitans tadpoles. We found that transmission occurred rapidly by direct and indirect horizontal transmission, but indirect transmission that included a potential substrate was more effective. Even though J. lividum colonized the skin, its relative abundance on the tadpole skin decreased over time. The inoculation of J. lividum did not significantly alter the skin bacterial diversity of L. clamitans tadpoles, which was dominated by Pseudomonas Our results show that indirect horizontal transmission can be an effective bioaugmentation method. Future research is needed to determine the best conditions, including the presence of substrates, under which a probiotic can persist on the skin so that bioaugmentation becomes a successful strategy to mitigate chytridiomycosis.

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

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

Amphibian skin bacteria display antifungal activity and induce plant defense mechanisms against Botrytis cinerea.

Botrytis cinerea is the causal agent of gray mold, which affects a wide variety of plant species. Chemical agents have been used to prevent the disease caused by this pathogenic fungus. However, their toxicity and reduced efficacy have encouraged the development of new biological control alternatives. Recent studies have shown that bacteria isolated from amphibian skin display antifungal activity against plant pathogens. However, the mechanisms by which these bacteria act to reduce the effects of B. cinerea are still unclear. From a diverse collection of amphibian skin bacteria, three proved effective in inhibiting the development of B. cinerea under in vitro conditions. Additionally, the individual application of each bacterium on the model plant Arabidopsis thaliana, Solanum lycopersicum and post-harvest blueberries significantly reduced the disease caused by B. cinerea. To understand the effect of bacteria on the host plant, we analyzed the transcriptomic profile of A. thaliana in the presence of the bacterium C32I and the fungus B. cinerea, revealing transcriptional regulation of defense-related hormonal pathways. Our study shows that bacteria from the amphibian skin can counteract the activity of B. cinerea by regulating the plant transcriptional responses.

An insulator embedded in the chicken α-globin locus regulates chromatin domain configuration and differential gene expression.

Genome organization into transcriptionally active domains denotes one of the first levels of gene expression regulation. Although the chromatin domain concept is generally accepted, only little is known on how domain organization impacts the regulation of differential gene expression. Insulators might hold answers to address this issue as they delimit and organize chromatin domains. We have previously identified a CTCF-dependent insulator with enhancer-blocking activity embedded in the 5' non-coding region of the chicken α-globin domain. Here, we demonstrate that this element, called the αEHS-1.4 insulator, protects a transgene against chromosomal position effects in stably transfected cell lines and transgenic mice. We found that this insulator can create a regulated chromatin environment that coincides with the onset of adult α-globin gene expression. Furthermore, such activity is in part dependent on the in vivo regulated occupancy of CTCF at the αEHS-1.4 element. Insulator function is also regulated by CTCF poly(ADP-ribosyl)ation. Our results suggest that the αEHS-1.4 insulator contributes in organizing the chromatin structure of the α-globin gene domain and prevents activation of adult α-globin gene expression at the erythroblast stage via CTCF.

Water-sediment niche differentiation in ancient marine lineages of Exiguobacterium endemic to the Cuatro Cienegas Basin.

The evolutionary history and ecological differentiation of the genus Exiguobacterium was characterized within natural communities from the Cuatro Cienegas Basin, Mexico. Exiguobacterium comprises both halophilic and alkaliphilic bacteria that are abundant among the aquatic systems of the Cuatro Cienegas Basin. We obtained complete sequences of the 16srRNA gene and partial sequences of four housekeeping genes (citC, rpoB, recA and hsp70) in 183 Exiguobacterium isolates retrieved from distinct aquatic systems. We defined three main phylogroups that are closely related to marine and thermophilic species of the genus. These phylogroups were neither specific to a given aquatic system nor to a particular salinity. Phylogenetic reconstruction indicated the presence of several small clusters within the phylogroups. These clusters consisted of isolates predominantly retrieved from sediment or water. Unifrac and AdaptML analyses confirmed this observation, pointing towards a clear pattern of differentiation linked to either sediment or water habitats. Our results are in line with the concept that niche differentiation is one of the main factors shaping prokaryotic populations and leading to evolutionary divergence.

The skin microbiota of the axolotl Ambystoma altamirani is highly influenced by metamorphosis and seasonality but not by pathogen infection.

BACKGROUND: Microbiomes have been increasingly recognized as major contributors to host health and survival. In amphibians, bacterial members of the skin microbiota protect their hosts by inhibiting the growth of the fungal pathogen Batrachochytrium dendrobatidis (Bd). Even though several studies describe the influence of biotic and abiotic factors over the skin microbiota, it remains unclear how these symbiotic bacterial communities vary across time and development. This is particularly relevant for species that undergo metamorphosis as it has been shown that host physiology and ecology drastically influence diversity of the skin microbiome. RESULTS: We found that the skin bacterial communities of the axolotl A. altamirani are largely influenced by the metamorphic status of the host and by seasonal variation of abiotic factors such as temperature, pH, dissolved oxygen and conductivity. Despite high Bd prevalence in these samples, the bacterial diversity of the skin microbiota did not differ between infected and non-infected axolotls, although relative abundance of particular bacteria were correlated with Bd infection intensity. CONCLUSIONS: Our work shows that metamorphosis is a crucial process that shapes skin bacterial communities and that axolotls under different developmental stages respond differently to environmental seasonal variations. Moreover, this study greatly contributes to a better understanding of the factors that shape amphibian skin microbiota, especially in a largely underexplored group like axolotls (Mexican Ambystoma species).

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,andCraugastor 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.

The fungal pathogen Batrachochytrium salamandrivorans is not detected in wild and captive amphibians from Mexico.

The recent emergence of the pathogen Batrachochytrium salamandrivorans (Bsal) is associated with rapid population declines of salamanders in Europe and its arrival to new areas could cause dramatic negative effects on other amphibian populations and species. Amphibian species, present in areas with high amphibian diversity such as Mexico, could be highly threatened due to the arrival of Bsal, particularly salamander species which are more vulnerable to chytridiomycosis caused by this pathogen. Thus, immediate surveillance is needed as a strategy to efficiently contend with this emerging infectious disease. In this study, we analyzed 490 wild and captive amphibians from 48 species across 76 sites in the North, Central, and South of Mexico to evaluate the presence of Bsal. Amphibians were sampled in sites with variable degrees of amphibian richness and suitability for Bsal according to previous studies. From the 76 sampling sites, 10 of them were located in areas with high amphibian richness and potential moderate to high Bsal habitat suitability. We did not detect Bsal in any of the samples, and no signs of the disease were observed in any individual at the time of sampling. Our results suggest that Bsal has not yet arrived at the sampled sites or could be at low prevalence within populations with low occurrence probability. This is the first study that evaluates the presence of Bsal in different regions and amphibian species in Mexico, which is the second most diverse country in salamander species in the world. We highlight the risk and the importance of continuing surveillance of Bsal in Mexico and discuss control strategies to avoid the introduction and spread of Bsal in the country.

Transcriptional Basis for Haustorium Formation and Host Establishment in Hemiparasitic Psittacanthus schiedeanus Mistletoes.

The mistletoe Psittacanthus schiedeanus, a keystone species in interaction networks between plants, pollinators, and seed dispersers, infects a wide range of native and non-native tree species of commercial interest. Here, using RNA-seq methodology we assembled the whole circularized quadripartite structure of P. schiedeanus chloroplast genome and described changes in the gene expression of the nuclear genomes across time of experimentally inoculated seeds. Of the 140,467 assembled and annotated uniGenes, 2,000 were identified as differentially expressed (DEGs) and were classified in six distinct clusters according to their expression profiles. DEGs were also classified in enriched functional categories related to synthesis, signaling, homoeostasis, and response to auxin and jasmonic acid. Since many orthologs are involved in lateral or adventitious root formation in other plant species, we propose that in P. schiedeanus (and perhaps in other rootless mistletoe species), these genes participate in haustorium formation by complex regulatory networks here described. Lastly, and according to the structural similarities of P. schiedeanus enzymes with those that are involved in host cell wall degradation in fungi, we suggest that a similar enzymatic arsenal is secreted extracellularly and used by mistletoes species to easily parasitize and break through tissues of the host.

Genomic consequences of dietary diversification and parallel evolution due to nectarivory in leaf-nosed bats.

BACKGROUND: The New World leaf-nosed bats (Phyllostomids) exhibit a diverse spectrum of feeding habits and innovations in their nutrient acquisition and foraging mechanisms. However, the genomic signatures associated with their distinct diets are unknown. RESULTS: We conducted a genomic comparative analysis to study the evolutionary dynamics related to dietary diversification and specialization. We sequenced, assembled, and annotated the genomes of five Phyllostomid species: one insect feeder (Macrotus waterhousii), one fruit feeder (Artibeus jamaicensis), and three nectar feeders from the Glossophaginae subfamily (Leptonycteris yerbabuenae, Leptonycteris nivalis, and Musonycteris harrisoni), also including the previously sequenced vampire Desmodus rotundus. Our phylogenomic analysis based on 22,388 gene families displayed differences in expansion and contraction events across the Phyllostomid lineages. Independently of diet, genes relevant for feeding strategies and food intake experienced multiple expansions and signatures of positive selection. We also found adaptation signatures associated with specialized diets: the vampire exhibited traits associated with a blood diet (i.e., coagulation mechanisms), whereas the nectarivore clade shares a group of positively selected genes involved in sugar, lipid, and iron metabolism. Interestingly, in fruit-nectar-feeding Phyllostomid and Pteropodids bats, we detected positive selection in two genes: AACS and ALKBH7, which are crucial in sugar and fat metabolism. Moreover, in these two proteins we found parallel amino acid substitutions in conserved positions exclusive to the tribe Glossophagini and to Pteropodids. CONCLUSIONS: Our findings illuminate the genomic and molecular shifts associated with the evolution of nectarivory and shed light on how nectar-feeding bats can avoid the adverse effects of diets with high glucose content.

Potential risk of Batrachochytrium salamandrivorans in Mexico.

The recent decline in populations of European salamanders caused by the chytrid fungus Batrachochytrium salamandrivorans (Bsal) has generated worldwide concern, as it is a major threat to amphibians. Evaluation of the areas most suitable for the establishment of Bsal combined with analysis of the distribution of salamander species could be used to generate and implement biosecurity measures and protect biodiversity at sites with high salamander diversity. In this study, we identified the areas most suitable for the establishment of Bsal in Mexico. Mexico has the second-highest salamander species diversity in the world; thus, we identified areas moderately to highly suitable for the establishment of Bsal with high salamander diversity as potential hotspots for surveillance. Central and Southern Mexico were identified as high-risk zones, with 13 hotspots where 30% of Mexican salamander species occur, including range-restricted species and endangered species. We propose that these hotspots should be thoroughly monitored for the presence of Bsal to prevent the spread of the pathogen if it is introduced to the country.

Tadpole body size and behaviour alter the social acquisition of a defensive bacterial symbiont.

Individual differences in host phenotypes can generate heterogeneity in the acquisition and transmission of microbes. Although this has become a prominent factor of disease epidemiology, host phenotypic variation might similarly underlie the transmission of microbial symbionts that defend against pathogen infection. Here, we test whether host body size and behaviour influence the social acquisition of a skin bacterium, Janthinobacterium lividum, which in some hosts can confer protection against infection by Batrachochytrium dendrobatidis, the causative agent of the amphibian skin disease chytridiomycosis. We measured body size and boldness (time spent in an open field) of green frog tadpoles and haphazardly constructed groups of six individuals. In some groups, we exposed one individual in each group to J. lividum and, in other groups, we inoculated a patch of aquarium pebbles to J. lividum. After 24 h, we swabbed each individual to estimate the presence of J. lividum on their skin. On average, tadpoles acquired nearly four times more bacteria when housed with an exposed individual compared to those housed with a patch of inoculated substrate. When tadpoles were housed with an exposed group-mate, larger and 'bolder' individuals acquired more bacteria. These data suggest that phenotypically biased acquisition of defensive symbionts might generate biased patterns of mortality from the pathogens against which they protect.

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.