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

Characterising the Tasmanian devil (Sarcophilus harrisii) pouch microbiome in lactating and non-lactating females

Wildlife harbour a diverse range of microorganisms that affect their health and development. Marsupials are born immunologically naïve and physiologically underdeveloped, with primary development occurring inside a pouch. Secretion of immunological compounds and antimicrobial peptides in the epithelial lining of the female's pouch, pouch young skin, and through the milk, are thought to boost the neonate's immune system and potentially alter the pouch skin microbiome. Here, using 16S rRNA amplicon sequencing, we characterised the Tasmanian devil pouch skin microbiome from 25 lactating and 30 non-lactating wild females to describe and compare across these reproductive stages. We found that the lactating pouch skin microbiome had significantly lower amplicon sequence variant richness and diversity than non-lactating pouches, however there was no overall dissimilarity in community structure between lactating and non-lactating pouches. The top five phyla were found to be consistent between both reproductive stages, with over 85% of the microbiome being comprised of Firmicutes, Proteobacteria, Fusobacteriota, Actinobacteriota, and Bacteroidota. The most abundant taxa remained consistent across all taxonomic ranks between lactating and non-lactating pouch types. This suggests that any potential immunological compounds or antimicrobial peptide secretions did not significantly influence the main community members. Of the more than 16,000 total identified amplicon sequence variants, 25 were recognised as differentially abundant between lactating and non-lactating pouches. It is proposed that the secretion of antimicrobial peptides in the pouch act to modulate these microbial communities. This study identifies candidate bacterial clades on which to test the activity of Tasmanian devil antimicrobial peptides and their role in pouch young protection, which in turn may lead to future therapeutic development for human diseases.

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.

Sex, health status and habitat alter the community composition and assembly processes of symbiotic bacteria in captive frogs

BACKGROUND

Frogs are critical economic animals essential to agricultural ecosystem equilibrium. However, Meningitis-like Infectious Disease (MID) often affects them in agricultural settings. While frog-associated microbiota contribute to elemental cycling and immunity, the effects of frog sex and health on gut bacteria remain understudied, and the relationship between frog habitat and soil microbes is unclear. We aimed to determine how frog sex, health status and habitat influence symbiotic bacteria and community assembly mechanism to provide guidance for sustainable frog farming and conservation.

RESULTS

We employed 16S rRNA sequencing to investigate gut microbiota differences in relation to frog sex and health status. We also compared symbiotic communities in frog-aggregation, native and soybean soil on the farm. Results showed that gut bacterial β-diversity and taxonomy were markedly influenced by frog sex and health. Healthy frogs had more robust gut bacterial metabolism than frogs infected with MID. Cooccurrence network analysis revealed that healthy female frogs had more complex microbial network structure than males; however, diseased males showed the greatest network complexity. The assembly mechanism of gut bacteria in male frogs was dominated by deterministic processes, whereas in female frogs it was dominated by stochastic processes. Among symbiotic bacteria in frog habitat soils, deterministic processes predominantly shaped the community assembly of soybean soil. In particular, soybean soil was enriched in pathogens and nitrogen functions, whereas frog-aggregation soil was markedly increased in sulphur respiration and hydrocarbon degradation.

CONCLUSION

Our study reveals that sex mainly alters the interaction network and assembly mechanism of frog intestinal bacteria; MID infection significantly inhibits the metabolic functions of intestinal bacteria. Furthermore, diverse frog habitat soils could shape more symbiotic bacteria to benefit frog farming. Our findings provide new horizons for symbiotic bacteria among frogs, which could contribute to sustainable agriculture and ecological balance.

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.

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

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

Immune priming prior to pathogen exposure sheds light on the relationship between host, microbiome and pathogen in disease

Dynamic interactions between host, pathogen and host-associated microbiome dictate infection outcomes. Pathogens including Batrachochytrium dendrobatidis (Bd) threaten global biodiversity, but conservation efforts are hindered by limited understanding of amphibian host, Bd and microbiome interactions. We conducted a vaccination and infection experiment using Eastern hellbender salamanders (Cryptobranchus alleganiensis alleganiensis) challenged with Bd to observe infection, skin microbial communities and gene expression of host skin, pathogen and microbiome throughout the experiment. Most animals survived high Bd loads regardless of their vaccination status and vaccination did not affect pathogen load, but host gene expression differed based on vaccination. Oral vaccination (exposure to killed Bd) stimulated immune gene upregulation while topically and sham-vaccinated animals did not significantly upregulate immune genes. In early infection, topically vaccinated animals upregulated immune genes but orally and sham-vaccinated animals downregulated immune genes. Bd increased pathogenicity-associated gene expression in late infection when Bd loads were highest. The microbiome was altered by Bd, but there was no correlation between anti-Bd microbe abundance or richness and pathogen burden. Our observations suggest that hellbenders initially generate a vigorous immune response to Bd, which is ineffective at controlling disease and is subsequently modulated. Interactions with antifungal skin microbiota did not influence disease progression.

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.

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

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

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

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

Batrachochytrium fungi: stealth invaders in amphibian skin

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

Host-Specificity and Core Taxa of Seagrass Leaf Microbiome Identified Across Tissue Age and Geographical Regions

The seagrass Zostera marina is a widespread foundational species in temperate coastal ecosystems that supports diverse communities of epiphytes and grazers. Bacteria link the production of seagrass to higher trophic levels and are thought to influence seagrass biology and health. Yet, we lack a clear understanding of the factors that structure the seagrass microbiome, or whether there is a consistent microbial community associated with seagrass that underpins functional roles. We sampled surface microbiome (epibiota) from new and old growth seagrass leaves and the surrounding seawater in eight meadows among four regions along the Central Coast of British Columbia, Canada to assess microbiome variability across space and as leaves age. We found that the seagrass leaf microbiome differs strongly from seawater. Microbial communities in new and old growth leaves are different from each other and from artificial seagrass leaves we deployed in one meadow. The microbiome on new leaves is less diverse and there is a small suite of core OTUs (operational taxonomic units) consistently present across regions. The overall microbial community for new leaves is more dispersed but with little regional differentiation, while the epiphytes on old leaves are regionally distinct. Many core OTUs on old leaves are commonly associated with marine biofilms. Together these observations suggest a stronger role for host filtering in new compared to old leaves, and a stronger influence of the environment and environmental colonization in old leaves. We found 11 core microbial taxa consistently present on old and new leaves and at very low relative abundance on artificial leaves and in the water column. These 11 taxa appear to be strongly associated with Z. marina. These core taxa may perform key functions important for the host such as detoxifying seagrass waste products, enhancing plant growth, and controlling epiphyte cover.

Advances in Microbiome Research for Animal Health

Host-associated microbiomes contribute in many ways to the homeostasis of the metaorganism. The microbiome's contributions range from helping to provide nutrition and aiding growth, development, and behavior to protecting against pathogens and toxic compounds. Here we summarize the current knowledge of the diversity and importance of the microbiome to animals, using representative examples of wild and domesticated species. We demonstrate how the beneficial ecological roles of animal-associated microbiomes can be generally grouped into well-defined main categories and how microbe-based alternative treatments can be applied to mitigate problems for both economic and conservation purposes and to provide crucial knowledge about host-microbiota symbiotic interactions. We suggest a Customized Combination of Microbial-Based Therapies to promote animal health and contribute to the practice of sustainable husbandry. We also discuss the ecological connections and threats associated with animal biodiversity loss, microorganism extinction, and emerging diseases, such as the COVID-19 pandemic.