Geography, Host Genetics, and Cross-Domain Microbial Networks Structure the Skin Microbiota of Fragmented Brazilian Atlantic Forest Frog Populations

Identification of keystone taxa in rhizosphere microbial communities using different methods and their effects on compounds of the host Cinnamomum migao

Exploring keystone taxa affecting microbial community stability and host function is crucial for understanding ecosystem functions. However, identifying keystone taxa from humongous microbial communities remains challenging. We collected 344 rhizosphere and bulk soil samples from the endangered plant C. migao for 2 years consecutively. Used high-throughput sequencing 16S rDNA and ITS to obtain the composition of bacterial and fungal communities. We explored keystone taxa and the applicability and limitations of five methods (SPEC-OCCU, Zi-Pi, Subnetwork, Betweenness, and Module), as well as the impact of microbial community domain, time series, and rhizosphere boundary on the identification of keystone taxa in the communities. Our results showed that the five methods, identified abundant keystone taxa in rhizosphere and bulk soil microbial communities. However, the keystone taxa shared by the rhizosphere and bulk soil microbial communities over time decreased rapidly decrease in the five methods. Among five methods on the identification of keystone taxa in the rhizosphere community, Module identified 113 taxa, SPEC-OCCU identified 17 taxa, Betweenness identified 3 taxa, Subnetwork identified 3 taxa, and Zi-Pi identified 4 taxa. The keystone taxa are mainly conditionally rare taxa, and their ecological functions include chemoheterotrophy, aerobic chemoheterotrophy, nitrate reduction, and anaerobic photoautotrophy. The results of the random forest model and structural equation model predict that keystone taxa Mortierella and Ellin6513 may have an effects on the accumulation of 1, 4, 7, - Cycloundecatriene, 1, 5, 9, 9-tetramethyl-, Z, Z, Z-, beta-copaene, bicyclogermacrene, 1,8-Cineole in C. migao fruits, but their effects still need further evidence. Our study evidence an unstable microbial community in the bulk soil, and the definition of microbial boundary and ecologically functional affected the identification of keystone taxa in the community. Subnetwork and Module are more in line with the definition of keystone taxa in microbial ecosystems in terms of maintaining community stability and hosting function.

Association between the skin microbiome and MHC class II diversity in an amphibian

Microbiomes play an important role in determining the ecology and behaviour of their hosts. However, questions remain pertaining to how host genetics shape microbiomes, and how microbiome composition influences host fitness. We explored the effects of geography, evolutionary history and host genetics on the skin microbiome diversity and structure in a widespread amphibian. More specifically, we examined the association between bacterial diversity and composition and the major histocompatibility complex class II exon 2 diversity in 12 moor frog (Rana arvalis) populations belonging to two geographical clusters that show signatures of past and ongoing differential selection. We found that while bacterial alpha diversity did not differ between the two clusters, MHC alleles/supertypes and genetic diversity varied considerably depending on geography and evolutionary history. Bacterial alpha diversity was positively correlated with expected MHC heterozygosity and negatively with MHC nucleotide diversity. Furthermore, bacterial community composition showed significant variation between the two geographical clusters and between specific MHC alleles/supertypes. Our findings emphasize the importance of historical demographic events on hologenomic variation and provide new insights into how immunogenetic host variability and microbial diversity may jointly influence host fitness with consequences for disease susceptibility and population persistence.

Host phylogeny and environment shape the diversity of salamander skin bacterial communities

The composition and diversity of animal-associated microbial communities are shaped by multiple ecological and evolutionary processes acting at different spatial and temporal scales. Skin microbiomes are thought to be strongly influenced by the environment due to the direct interaction of the host's skin with the external media. As expected, the diversity of amphibian skin microbiomes is shaped by climate and host sampling habitats, whereas phylogenetic effects appear to be weak. However, the relative strength of phylogenetic and environmental effects on salamander skin microbiomes remains poorly understood. Here, we analysed sequence data from 1164 adult salamanders of 44 species to characterise and compare the diversity and composition of skin bacteria. We assessed the relative contribution of climate, host sampling habitat, and host phylogeny to the observed patterns of bacterial diversity. We found that bacterial alpha diversity was mainly associated with host sampling habitat and climate, but that bacterial beta diversity was more strongly associated with host taxonomy and phylogeny. This phylogenetic effect predominantly occurred at intermediate levels of host divergence (0-50 Mya). Our results support the importance of environmental factors shaping the diversity of salamander skin microbiota, but also support host phylogenetic history as a major factor shaping these bacterial communities.

Geographic variation in bacterial assemblages on cane toad skin is influenced more by local environments than by evolved changes in host traits

Bacterial assemblages on amphibian skin may play an important role in protecting hosts against infection. In hosts that occur over a range of environments, geographic variation in composition of bacterial assemblages might be due to direct effects of local factors and/or to evolved characteristics of the host. Invasive cane toads (Rhinella marina) are an ideal candidate to evaluate environmental and genetic mechanisms, because toads have evolved major shifts in physiology, morphology, and behavior during their brief history in Australia. We used samples from free-ranging toads to quantify site-level differences in bacterial assemblages and a common-garden experiment to see if those differences disappeared when toads were raised under standardised conditions at one site. The large differences in bacterial communities on toads from different regions were not seen in offspring raised in a common environment. Relaxing bacterial clustering to operational taxonomic units in place of amplicon sequence variants likewise revealed high similarity among bacterial assemblages on toads in the common-garden study, and with free-ranging toads captured nearby. Thus, the marked geographic divergence in bacterial assemblages on wild-caught cane toads across their Australian invasion appears to result primarily from local environmental effects rather than evolved shifts in the host.

Assessing the causes and consequences of gut mycobiome variation in a wild population of the Seychelles warbler

BACKGROUND

Considerable research has focussed on the importance of bacterial communities within the vertebrate gut microbiome (GM). However, studies investigating the significance of other microbial kingdoms, such as fungi, are notably lacking, despite their potential to influence host processes. Here, we characterise the fungal GM of individuals living in a natural population of Seychelles warblers (Acrocephalus sechellensis). We evaluate the extent to which fungal GM structure is shaped by environment and host factors, including genome-wide heterozygosity and variation at key immune genes (major histocompatibility complex (MHC) and Toll-like receptor (TLR)). Importantly, we also explore the relationship between fungal GM differences and subsequent host survival. To our knowledge, this is the first time that the genetic drivers and fitness consequences of fungal GM variation have been characterised for a wild vertebrate population.

RESULTS

Environmental factors, including season and territory quality, explain the largest proportion of variance in the fungal GM. In contrast, neither host age, sex, genome-wide heterozygosity, nor TLR3 genotype was associated with fungal GM differences in Seychelles warblers. However, the presence of four MHC-I alleles and one MHC-II allele was associated with changes in fungal GM alpha diversity. Changes in fungal richness ranged from between 1 and 10 sequencing variants lost or gained; in some cases, this accounted for 20% of the fungal variants carried by an individual. In addition to this, overall MHC-I allelic diversity was associated with small, but potentially important, changes in fungal GM composition. This is evidenced by the fact that fungal GM composition differed between individuals that survived or died within 7 months of being sampled.

CONCLUSIONS

Our results suggest that environmental factors play a primary role in shaping the fungal GM, but that components of the host immune system-specifically the MHC-may also contribute to the variation in fungal communities across individuals within wild populations. Furthermore, variation in the fungal GM can be associated with differential survival in the wild. Further work is needed to establish the causality of such relationships and, thus, the extent to which components of the GM may impact host evolution. Video Abstract.

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.

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).

Signatures of functional bacteriome structure in a tropical direct-developing amphibian species

Background

Host microbiomes may differ under the same environmental conditions and these differences may influence susceptibility to infection. Amphibians are ideal for comparing microbiomes in the context of disease defense because hundreds of species face infection with the skin-invading microbe Batrachochytrium dendrobatidis (Bd), and species richness of host communities, including their skin bacteria (bacteriome), may be exceptionally high. We conducted a landscape-scale Bd survey of six co-occurring amphibian species in Brazil’s Atlantic Forest. To test the bacteriome as a driver of differential Bd prevalence, we compared bacteriome composition and co-occurrence network structure among the six focal host species.

Results

Intensive sampling yielded divergent Bd prevalence in two ecologically similar terrestrial-breeding species, a group with historically low Bd resistance. Specifically, we detected the highest Bd prevalence in Ischnocnema henselii but no Bd detections in Haddadus binotatus. Haddadus binotatus carried the highest bacteriome alpha and common core diversity, and a modular network partitioned by negative co-occurrences, characteristics associated with community stability and competitive interactions that could inhibit Bd colonization.

Conclusions

Our findings suggest that community structure of the bacteriome might drive Bd resistance in H. binotatus, which could guide microbiome manipulation as a conservation strategy to protect diverse radiations of direct-developing species from Bd-induced population collapses.

Metabarcoding approaches in amphibian disease ecology: Disentangling the functional contributions of skin bacteria on disease outcome

Molecular technologies have revolutionized the field of wildlife disease ecology, allowing the detection of outbreaks, novel pathogens, and invasive strains. In particular, metabarcoding approaches, defined here as tools used to amplify and sequence universal barcodes from a single sample (e.g., 16S rRNA for bacteria, ITS for fungi, 18S rRNA for eukaryotes), are expanding our traditional view of host-pathogen dynamics by integrating microbial interactions that modulate disease outcome. Here, I provide an analysis from the perspective of the field of amphibian disease ecology, where the emergence of multi-host pathogens has caused global declines and species extinctions. I re-analyzed an experimental mesocosm dataset to infer the functional profiles of the skin microbiomes of coqui frogs (Eleutherodactylus coqui), an amphibian species that is consistently found infected with the fungal pathogen Batrachochytrium dendrobatidis and has high turnover of skin bacteria driven by seasonal shifts. I found that the metabolic activities of microbiomes operate at different capacities depending on the season. Global enrichment of predicted functions was more prominent during the warm-wet season, indicating that microbiomes during the cool-dry season were either depauperate, resistant to new bacterial colonization, or that their functional space was more saturated. These findings suggest important avenues to investigate how microbes regulate population growth and contribute to host physiological processes. Overall, this study highlights the current challenges and future opportunities in the application of metabarcoding to investigate the causes and consequences of disease in wild systems.