Plethodontid salamanders show variable disease dynamics in response to Batrachochytrium salamandrivorans chytridiomycosis

Geographic risk assessment of Batrachochytrium salamandrivorans invasion in Costa Rica as a means of informing emergence management and mitigation

Remotely-sensed risk assessments of emerging, invasive pathogens are key to targeted surveillance and outbreak responses. The recent emergence and spread of the fungal pathogen, Batrachochytrium salamandrivorans (Bsal), in Europe has negatively impacted multiple salamander species. Scholars and practitioners are increasingly concerned about the potential consequences of this lethal pathogen in the Americas, where salamander biodiversity is higher than anywhere else in the world. Although Bsal has not yet been detected in the Americas, certain countries have already proactively implemented monitoring and detection plans in order to identify areas of greatest concern and enable efficient contingency planning in the event of pathogen detection. To predict areas in Costa Rica with a high Bsal transmission risk, we employed ecological niche modeling combined with biodiversity and tourist visitation data to ascertain the specific risk to a country with world renowned biodiversity. Our findings indicate that approximately 23% of Costa Rica's landmass provides suitable conditions for Bsal, posing a threat to 37 salamander species. The Central and Talamanca mountain ranges, in particular, have habitats predicted to be highly suitable for the pathogen. To facilitate monitoring and mitigation efforts, we identified eight specific protected areas that we believe are at the greatest risk due to a combination of high biodiversity, tourist visitation, and suitable habitat for Bsal. We advise regular monitoring utilizing remotely-sensed data and ecological niche modeling to effectively target in-situ surveillance and as places begin implementing educational efforts.

Invasibility of a North American soil ecosystem to amphibian-killing fungal pathogens

North American salamanders are threatened by intercontinental spread of chytridiomycosis, a deadly disease caused by the fungal pathogen Batrachochytrium salamandrivorans (Bsal). To predict potential dispersal of Bsal spores to salamander habitats, we evaluated the capacity of soil microbial communities to resist invasion. We determined the degree of habitat invasibility using soils from five locations throughout the Great Smoky Mountains National Park, a region with a high abundance of susceptible hosts. Our experimental design consisted of replicate soil microcosms exposed to different propagule pressures of the non-native pathogen, Bsal, and an introduced but endemic pathogen, B. dendrobatidis (Bd). To compare growth and competitive interactions, we used quantitative PCR, live/dead cell viability assays, and full-length 16S rRNA sequencing. We found that soil microcosms with intact bacterial communities inhibited both Bsal and Bd growth, but inhibitory capacity diminished with increased propagule pressure. Bsal showed greater persistence than Bd. Linear discriminant analysis (LDA) identified the family Burkolderiaceae as increasing in relative abundance with the decline of both pathogens. Although our findings provide evidence of environmental filtering in soils, such barriers weakened in response to pathogen type and propagule pressure, showing that habitats vary their invasibility based on properties of their local microbial communities.

Ecological Barriers for an Amphibian Pathogen: A Narrow Ecological Niche for Batrachochytrium salamandrivorans in an Asian Chytrid Hotspot

The chytrid fungal pathogens Batrachochytrium salamandrivorans (Bsal) and B. dendrobatidis (Bd) are driving amphibian extinctions and population declines worldwide. As their origins are believed to be in East/Southeast Asia, this region is crucial for understanding their ecology. However, Bsal screening is relatively limited in this region, particularly in hotspots where Bd lineage diversity is high. To address this gap, we conducted an extensive Bsal screening involving 1101 individuals from 36 amphibian species, spanning 17 natural locations and four captive facilities in the biodiversity-rich Guangxi Zhuang Autonomous Region (GAR). Our PCR assays yielded unexpected results, revealing the complete absence of Bsal in all tested samples including 51 individuals with Bd presence. To understand the potential distribution of Bsal, we created niche models, utilizing existing occurrence records from both Asia and Europe. These models estimated potential suitable habitats for Bsal largely in the northern and southwestern parts of the GAR. Although Bsal was absent in our samples, the niche models identified 10 study sites as being potentially suitable for this pathogen. Interestingly, out of these 10 sites, Bd was detected at 8. This suggests that Bsal and Bd could possibly co-exist in these habitats, if Bsal were present. Several factors seem to influence the distribution of Bsal in Asia, including variations in temperature, local caudate species diversity, elevation, and human population density. However, it is climate-related factors that hold the greatest significance, accounting for a notable 60% contribution. The models propose that the specific climatic conditions of arid regions, primarily seen in the GAR, play a major role in the distribution of Bsal. Considering the increased pathogenicity of Bsal at stable and cooler temperatures (10-15 °C), species-dependent variations, and the potential for seasonal Bd-Bsal interactions, we emphasize the importance of periodic monitoring for Bsal within its projected range in the GAR. Our study provides deeper insights into Bsal's ecological niche and the knowledge generated will facilitate conservation efforts in amphibian populations devastated by chytrid pathogens across other regions of the world.

Bacterial Skin Assemblages of Sympatric Salamanders Are Primarily Shaped by Host Genus

Bacterial assemblages on the skins of amphibians are known to influence pathogen resistance and other important physiological functions in the host. Host-specific factors and the environment play significant roles in structuring skin assemblages. This study used high-throughput 16S rRNA sequencing and multivariate analyses to examine differences in skin-bacterial assemblages from 246 salamanders belonging to three genera in the lungless family Plethodontidae along multiple spatial gradients. Composition and α- and β-diversity of bacterial assemblages were defined, indicator species were identified for each host group, and the relative influences of host- versus environment-specific ecological factors were evaluated. At the broadest spatial scale, host genus, host species, and sampling site were predictive of skin assemblage structure, but host genus and species were more influential after controlling for the marginal effects of site, as well as nestedness of site. Furthermore, assemblage similarity within each host genus did not change with increasing geographic distance. At the smallest spatial scale, site-specific climate analyses revealed different relationships to climatic variables for each of the three genera, and these relationships were determined by host ecomode. Variation in bacterial assemblages of terrestrial hosts correlated with landscape-level climatic variability, and this pattern decayed with increasing water dependence of the host. Results from this study highlight host-specific considerations for researchers studying wildlife diseases in co-occurring, yet ecologically divergent, species.

Inhibitory Bacterial Diversity and Mucosome Function Differentiate Susceptibility of Appalachian Salamanders to Chytrid Fungal Infection

Mucosal defenses are crucial in animals for protection against pathogens and predators. Host defense peptides (antimicrobial peptides, AMPs) as well as skin-associated microbes are key components of mucosal immunity, particularly in amphibians. We integrate microbiology, molecular biology, network-thinking, and proteomics to understand how host and microbially derived products on amphibian skin (referred to as the mucosome) serve as pathogen defenses. We studied defense mechanisms against chytrid pathogens, Batrachochytrium dendrobatidis (Bd) and B. salamandrivorans (Bsal), in four salamander species with different Batrachochytrium susceptibilities. Bd infection was quantified using qPCR, mucosome function (i.e., ability to kill Bd or Bsal zoospores in vitro), skin bacterial communities using 16S rRNA gene amplicon sequencing, and the role of Bd-inhibitory bacteria in microbial networks across all species. We explored the presence of candidate-AMPs in eastern newts and red-backed salamanders. Eastern newts had the highest Bd prevalence and mucosome function, while red-back salamanders had the lowest Bd prevalence and mucosome function, and two-lined salamanders and seal salamanders were intermediates. Salamanders with highest Bd infection intensity showed greater mucosome function. Bd infection prevalence significantly decreased as putative Bd-inhibitory bacterial richness and relative abundance increased on hosts. In co-occurrence networks, some putative Bd-inhibitory bacteria were found as hub-taxa, with red-backs having the highest proportion of protective hubs and positive associations related to putative Bd-inhibitory hub bacteria. We found more AMP candidates on salamanders with lower Bd susceptibility. These findings suggest that salamanders possess distinct innate mechanisms that affect chytrid fungi. IMPORTANCE How host mucosal defenses interact, and influence disease outcome is critical in understanding host defenses against pathogens. A more detailed understanding is needed of the interactions between the host and the functioning of its mucosal defenses in pathogen defense. This study investigates the variability of chytrid susceptibility in salamanders and the innate defenses each species possesses to mediate pathogens, thus advancing the knowledge toward a deeper understanding of the microbial ecology of skin-associated bacteria and contributing to the development of bioaugmentation strategies to mediate pathogen infection and disease. This study improves the understanding of complex immune defense mechanisms in salamanders and highlights the potential role of the mucosome to reduce the probability of Bd disease development and that putative protective bacteria may reduce likelihood of Bd infecting skin.

Widespread occurrence of Batrachochytrium dendrobatidis in Ontario, Canada, and predicted habitat suitability for the emerging Batrachochytrium salamandrivorans

Chytridiomycosis, caused by the fungi Batrachochytrium dendrobatidis and Batrachochytrium salamandrivorans, is associated with massive amphibian mortality events worldwide and with some species’ extinctions. Previous ecological niche models suggest that B. dendrobatidis is not well‐suited to northern, temperate climates, but these predictions have often relied on datasets in which northern latitudes are underrepresented. Recent northern detections of B. dendrobatidis suggest that these models may have underestimated the suitability of higher latitudes for this fungus. We used qPCR to test for B. dendrobatidis in 1,041 non‐invasive epithelial swab samples from 18 species of amphibians collected across 735,345 km² in Ontario and Akimiski Island (Nunavut), Canada. We detected the pathogen in 113 samples (10.9%) from 11 species. Only one specimen exhibited potential clinical signs of disease. We used these data to produce six Species Distribution Models of B. dendrobatidis, which classified half of the study area as potential habitat for the fungus. We also tested each sample for B. salamandrivorans, an emerging pathogen that is causing alarming declines in European salamanders, but is not yet detected in North America. We did not detect B. salamandrivorans in any of the samples, providing a baseline for future surveillance. We assessed the potential risk of future introduction by comparing salamander richness to temperature‐dependent mortality, predicted by a previous exposure study. Areas with the highest species diversity and predicted mortality risk extended 60,530 km² across southern Ontario, highlighting the potential threat B. salamandrivorans poses to northern Nearctic amphibians. Preventing initial introduction will require coordinated, transboundary regulation of trade in amphibians (including frogs that can carry and disperse B. salamandrivorans), and surveillance of the pathways of introduction (e.g., water and wildlife). Our results can inform surveillance for both pathogens and efforts to mitigate the spread of chytridiomycosis through wild populations.