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

Phylogenomics as baseline for taxonomy description: Amphibiibacter pelophylacis gen. nov., sp. nov., a novel taxon of the family Sphaerotilaceae, class Betaproteobacteria, isolated from the skin microbiota of Pelophylax perezi from different populations

Bacterial strain SL12-8T was characterized and isolated from the skin microbiota of Pelophylax perezi, the Perez's frog. Strain SL12-8T stained Gram-negative and formed rod-shaped cells that grew optimally at 25 °C and pH 7.0-7.5. The G+C content of the DNA was 66.2 mol%. Ubiquinone 8 was the respiratory quinone identified in the studied strain and the most closely related taxon. The major fatty acids were summed in feature 3 and summed in feature 8 and C16:0, representing 84% of the total fatty acids. Phylogenetic analyses based on the 16S rRNA gene sequence placed strain SL12-8T within the order Burkholderiales in a distinct lineage. The 16S rRNA gene sequence similarities of strain SL12-8T to that of Rubrivivax albus, Scleromatobacter humisilvae, Piscinibacter aquaticus, Azohydromonas caseinilytica and Aquincola agrisoli were 94.41, 94.08, 93.72, 93.72 and 93.64%, respectively. The draft genome sequence of strain SL12-8T comprises 3,115,197 bases with a 313-fold mapped genome coverage. The assembled genome consists of 53 large contigs with more than 500 bp, and the genome encodes 2,814 putative coding sequences. The analysis of the available genomes from the closest genera showed 124 core genes that reveal a novel genus-level clade including the strain SL12-8T. Analysis of the SL12-8T genome revealed the presence of the beta-lactone and terpene biosynthetic gene clusters. The phylogenomic, phylogenetic, phenotypic and chemotaxonomic data showed that strain SL12-8T (=UCCCB 131T=CECT 30762T) represents the type of a novel species and genus, for which we propose the name Amphibiibacter pelophylacis gen. nov., sp. nov.

Amphibian skin bacteria contain a wide repertoire of genes linked to their antifungal capacities

Emerging diseases caused by fungi are a serious threat to wildlife biodiversity. The widespread fungal pathogen Batrachochytrium dendrobatidis (Bd) has caused dramatic amphibian population declines and species extinctions worldwide. While many amphibians have been negatively affected by Bd, some populations/species have persisted despite its presence. One factor contributing to amphibian protection against this fungus is the host-associated skin microbiome. In this study, we aimed to identify gene clusters associated with the antifungal activity of amphibian skin bacteria. Specifically, we explored skin bacteria isolated from species that have persisted in the wild despite the presence of Bd: the frogs Agalychnis callidryas, Craugastor fitzingeri, Dendropsophus ebraccatus, and the axolotl Ambystoma altamirani. Bacterial isolates were tested in vitro for their capacity to inhibit the growth of two fungal pathogens: Bd and the phytopathogen Botrytis cinerea (Bc). Genome mining of these bacterial isolates revealed a diverse repertoire of Biosynthetic Gene Clusters (BGCs) and chitin-degrading gene families (ChDGFs) whose composition and abundance differed among bacterial families. We found specific BGCs and ChDGFs that were associated with the capacity of bacteria to inhibit the growth of either Bd or Bc, suggesting that distinct fungi could be inhibited by different molecular mechanisms. By using similarity networks and machine learning, we identified BGCs encoding known antifungal compounds such as viscosin, fengycin, zwittermicin, as well as siderophores and a novel family of beta-lactones. Finally, we propose that the diversity of BGCs found in amphibian skin bacteria comprise a substantial genetic reservoir that could collectively explain the antifungal properties of the amphibian skin microbiome.

Priority effects alter microbiome composition and increase abundance of probiotic taxa in treefrog tadpoles

Host-associated microbial communities, like other ecological communities, may be impacted by the colonization order of taxa through priority effects. Developing embryos and their associated microbiomes are subject to stochasticity during colonization by bacteria. For amphibian embryos, often developing externally in bacteria-rich environments, this stochasticity may be particularly impactful. For example, the amphibian microbiome can mitigate lethal outcomes from disease for their hosts; however, this may depend on microbiome composition. Here, we examined the assembly of the bacterial community in spring peeper (Pseudacris crucifer) embryos and tadpoles. First, we reared embryos from identified mating pairs in either lab or field environments to examine the relative impact of environment and parentage on embryo and tadpole bacterial communities. Second, we experimentally inoculated embryos to determine if priority effects (i) could be used to increase the relative abundance of Janthinobacterium lividum, an amphibian-associated bacteria capable of preventing fungal infection, and (ii) would lead to observed differences in the relative abundances of two closely related bacteria from the genus Pseudomonas. Using 16S rRNA gene amplicon sequencing, we observed differences in community composition based on rearing location and parentage in embryos and tadpoles. In the inoculation experiment, we found that priority inoculation could increase the relative abundance of J. lividum, but did not find that either Pseudomonas isolate was able to prevent colonization by the other when given priority. These results highlight the importance of environmental source pools and parentage in determining microbiome composition, while also providing novel methods for the administration of a known amphibian probiotic.

IMPORTANCE

Harnessing the functions of host-associated bacteria is a promising mechanism for managing disease outcomes across different host species. In the case of amphibians, certain frog-associated bacteria can mitigate lethal outcomes of infection by the fungal pathogen Batrachochytrium dendrobatidis. Successful probiotic applications require knowledge of community assembly and an understanding of the ecological mechanisms that structure these symbiotic bacterial communities. In our study, we show the importance of environment and parentage in determining bacterial community composition and that community composition can be influenced by priority effects. Further, we provide support for the use of bacterial priority effects as a mechanism to increase the relative abundance of target probiotic taxa in a developing host. While our results show that priority effects are not universally effective across all host-associated bacteria, our ability to increase the relative abundance of specific probiotic taxa may enhance conservation strategies that rely on captive rearing of endangered vertebrates.

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.