Nor climate nor human impact factors: Chytrid infection shapes the skin bacterial communities of an endemic amphibian in a biodiversity hotspot.

The bacterial communities of the amphibian skin (i.e., the bacteriome) are critical to the host's innate immune system. However, it is unclear how different drivers can alter this function by modulating the bacteriome's structure. Our aim was to assess the extent to which different host attributes and extrinsic factors influence the structure of the bacterial communities of the skin. Skin bacterial diversity was examined in 148 individuals of the four-eyed frog (Pleurodema thaul) from 16 localities spanning almost 1800 km in latitude. The richness and beta diversity of bacterial families and the richness and abundance of Bd-inhibitory bacterial genera were used to describe their structure. Predictors associated with the host (developmental stage, genetic lineage, individual Batrachochytrium dendrobatidis [Bd] infection status) and the landscape (current climate, degree of anthropogenic disturbance) were used in the statistical modeling in an information theoretical approach. Bd infection and host developmental stage were the only predictors affecting bacteriome richness, with Bd+ individuals and postmetamorphic stages (adults and juveniles) having higher richness than Bd- ones and tadpoles. High diversity in Bd+ individuals is not driven by bacterial genera with known anti-Bd properties. Beta diversity was not affected by Bd infection and was mainly a consequence of bacterial family turnover rather than nestedness. Finally, for those bacterial genera known to have inhibitory effects on chytrid, Bd+ individuals had a slightly higher diversity than Bd- ones. Our study confirms an association between Bd infection and the host developmental stage with the bacterial communities of the skin of P. thaul. Unexpectedly, macroclimate and human impact factors do not seem to play a role in shaping the amphibian skin microbiome. Our study exemplifies that focusing on a single host-parasite system over a large geographic scale can provide essential insights into the factors driving host-parasite-bacteriome interactions.

Genomic characterization of two ranavirus isolates identified from a gopher frog (Lithobates capito) and a striped newt (Notophthalmus perstriatus) during a mass mortality event in Florida.

Two ranavirus isolates were recovered from anuran and salamander samples collected during an amphibian mass mortality event in North-Central Florida in 2021. Phylogenetic analyses of the full genomes confirmed that the two isolates were nearly identical and strains of the species Frog virus 3.

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.

Sequence capture identifies fastidious chytrid fungi directly from host tissue.

The chytrid fungus Batrachochytrium dendrobatidis (Bd) was discovered in 1998 as the cause of chytridiomycosis, an emerging infectious disease causing mass declines in amphibian populations worldwide. The rapid population declines of the 1970s-1990s were likely caused by the spread of a highly virulent lineage belonging to the Bd-GPL clade that was introduced to naïve susceptible populations. Multiple genetically distinct and regional lineages of Bd have since been isolated and sequenced, greatly expanding the known biological diversity within this fungal pathogen. To date, most Bd research has been restricted to the limited number of samples that could be isolated using culturing techniques, potentially causing a selection bias for strains that can grow on media and missing other unculturable or fastidious strains that are also present on amphibians. We thus attempted to characterize potentially non-culturable genetic lineages of Bd from distinct amphibian taxa using sequence capture technology on DNA extracted from host tissue and swabs. We focused our efforts on host taxa from two different regions that likely harbored distinct Bd clades: (1) wild-caught leopard frogs (Rana) from North America, and (2) a Japanese Giant Salamander (Andrias japonicus) at the Smithsonian Institution's National Zoological Park that exhibited signs of disease and tested positive for Bd using qPCR, but multiple attempts failed to isolate and culture the strain for physiological and genetic characterization. We successfully enriched for and sequenced thousands of fungal genes from both host clades, and Bd load was positively associated with number of recovered Bd sequences. Phylogenetic reconstruction placed all the Rana-derived strains in the Bd-GPL clade. In contrast, the A. japonicus strain fell within the Bd-Asia3 clade, expanding the range of this clade and generating additional genomic data to confirm its placement. The retrieved ITS locus matched public barcoding data from wild A. japonicus and Bd infections found on other amphibians in India and China, suggesting that this uncultured clade is widespread across Asia. Our study underscores the importance of recognizing and characterizing the hidden diversity of fastidious strains in order to reconstruct the spatiotemporal and evolutionary history of Bd. The success of the sequence capture approach highlights the utility of directly sequencing pathogen DNA from host tissue to characterize cryptic diversity that is missed by culture-reliant approaches.

Ontogeny drives shifts in skin bacterial communities in facultatively paedomorphic salamanders.

Microbiomes are major determinants of host growth, development and survival. In amphibians, host-associated bacteria in the skin can inhibit pathogen infection, but many processes can influence the structure and composition of the community. Here we quantified the shifts in skin-associated bacteria across developmental stages in the striped newt (Notophthalmus perstriatus), a threatened salamander species with a complex life history and vulnerable to infection by the amphibian chytrid fungus Batrachochytrium dendrobatidis and ranavirus. Our analyses show that pre-metamorphic larval and paedomorphic stages share similar bacterial compositions, and that the changes in the microbiome coincided with physiological restructuring during metamorphosis. Newts undergoing metamorphosis exhibited microbiome compositions that were intermediate between paedomorphic and post-metamorphic stages, further supporting the idea that metamorphosis is a major driver of host-associated microbes in amphibians. We did not find support for infection-related disruption of the microbiome, though infection replicates were small for each respective life stage.

Evolutionary ecology of host competence after a chytrid outbreak in a naive amphibian community.

Naive multi-host communities include species that may differentially maintain, transmit and amplify novel pathogens; therefore, we expect species to fill distinct roles during infectious disease emergence. Characterizing these roles in wildlife communities is challenging because most disease emergence events are unpredictable. Here, we used field-collected data to investigate how species-specific attributes influenced the degree of exposure, probability of infection, and pathogen intensity, during the emergence of the fungal pathogen Batrachochytrium dendrobatidis (Bd) in a highly diverse tropical amphibian community. Our findings confirmed that ecological traits commonly evaluated as correlates of decline were positively associated with infection prevalence and intensity at the species level during the outbreak. We identified key hosts that disproportionally contributed to transmission dynamics in this community and found a signature of phylogenetic history in disease responses associated with increased pathogen exposure via shared life-history traits. Our findings establish a framework that could be applied in conservation efforts to identify key species driving disease dynamics under enzootics before reintroducing amphibians back into their original communities. Reintroductions of supersensitive hosts that are unable to overcome infections will limit the success of conservation programmes by amplifying the disease at the community level. This article is part of the theme issue 'Amphibian immunity: stress, disease and ecoimmunology'.

Linking pathogen-microbiome-host interactions to explain amphibian population dynamics.

Symbiotic interactions can determine the evolutionary trajectories of host species, influencing genetic variation through selection and changes in demography. In the context of strong selective pressures such as those imposed by infectious diseases, symbionts providing defences could contribute to increase host fitness upon pathogen emergence. Here, we generated genome-wide data of an amphibian species to find evidence of evolutionary pressures driven by two skin symbionts: a batrachochytrid fungal pathogen and an antifungal bacterium. Using demographic modelling, we found evidence of decreased effective population size, probably due to pathogen infections. Additionally, we investigated host genetic associations with infection status, antifungal bacterium abundance and overall microbiome diversity using structural equation models. We uncovered relatively lower nucleotide diversity in infected frogs and potential heterozygote advantage to recruit the candidate beneficial symbiont and fight infections. Our models indicate that environmental conditions have indirect effects on symbiont abundance through both host body traits and microbiome diversity. Likewise, we uncovered a potential offsetting effect among host heterozygosity-fitness correlations, plausibly pointing to different ecological and evolutionary processes among the three species due to dynamic interactions. Our findings revealed that evolutionary pressures not only arise from the pathogen but also from the candidate beneficial symbiont, and both interactions shape the genetics of the host. Our results advance knowledge about multipartite symbiotic relationships and provide a framework to model ecological and evolutionary dynamics in wild populations. Finally, our study approach can be applied to inform conservation actions such as bioaugmentation strategies for other imperilled amphibians affected by infectious diseases.

Panzootic chytrid fungus exploits diverse amphibian host environments through plastic infection strategies.

While some pathogens are limited to single species, others can colonize many hosts, likely contributing to the emergence of novel disease outbreaks. Despite this biodiversity threat, traits associated with host niche expansions are not well understood in multihost pathogens. Here, we aimed to uncover functional machinery driving multihost invasion by focusing on Batrachochytrium dendrobatidis (Bd), a pathogen that infects the skin of hundreds of amphibians worldwide. We performed a meta-analysis of Bd gene expression using data from published infection experiments and newly generated profiles. We analysed Bd transcriptomic landscapes across the skin of 14 host species, reconstructed Bd isolates phylogenetic relationships, and inferred the origin and evolutionary history of differentially expressed genes under a phylogenetic framework comprising other 12 zoosporic fungi. Bd displayed plastic infection strategies when challenged by hosts with different disease susceptibility. Our analyses identified sets of differentially expressed genes under host environments with similar infection outcome. We stressed nutritional immunity and gene silencing as important processes required to overcome challenging skin environments in less susceptible hosts. Overall, Bd genes expressed during amphibian skin exploitation have arisen mainly via gene duplications with great family expansions, increasing the gene copy events previously described for this fungal species. Finally, we provide a comprehensive gene data set that can be used to further examine eco-evolutionary hypotheses for this host-pathogen system. Our study supports the idea that host environments exert contrasting selective pressures, such that gene expression plasticity could be one of the evolutionary keys leading to the success of multihost pathogens.

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

Sustained Ranavirus Outbreak Causes Mass Mortality and Morbidity of Imperiled Amphibians in Florida.

A persistent 2-month long outbreak of Ranavirus in a natural community of amphibians contributed to a mass die-off of gopher frog tadpoles (Lithobates capito) and severe disease in striped newts (Notophthalmus perstriatus) in Florida. Ongoing mortality in L. capito and disease signs in N. perstriatus continued for 5 weeks after the first observation. Hemorrhagic disease and necrosis were diagnosed from pathological examination of L. capito tadpoles. We confirmed detection of a frog virus 3 (FV3)-like Ranavirus via quantitative PCR in all species. Our findings highlight the susceptibility of these species to Rv and the need for long-term disease surveillance during epizootics.

Preparing for invasion: Assessing risk of infection by chytrid fungi in southeastern plethodontid salamanders.

Understanding the responses of naïve communities to the invasion of multihost pathogens requires accurate estimates of susceptibility across taxa. In the Americas, the likely emergence of a second amphibian pathogenic fungus (Batrachochytrium salamandrivorans, Bsal) calls for new ways of prioritizing disease mitigation among species due to the high diversity of naïve hosts with prior B. dendrobatidis (Bd) infections. Here, we applied the concept of pathogenic potential to quantify the virulence of chytrid fungi on naïve amphibians and evaluate species for conservation efforts in the event of an outbreak. The benefit of this measure is that it combines and summarizes the variation in disease effects into a single numerical index, allowing for comparisons across species, populations or groups of individuals that may inherently exhibit differences in susceptibility. As a proof of concept, we obtained standardized responses of disease severity by performing experimental infections with Bsal on five plethodontid salamanders from southeastern United States. Four out of five species carried natural infections of Bd at the start of the experiments. We showed that Bsal exhibited its highest value of pathogenic potential in a species that is already declining (Desmognathus auriculatus). We find that this index provides additional information beyond the standard measures of disease prevalence, intensity, and mortality, because it leveraged these disease parameters within each categorical group. Scientists and practitioners could use this measure to justify research, funding, trade, or conservation measures.

Incapacitating effects of fungal coinfection in a novel pathogen system.

As globalization lowers geographic barriers to movement, coinfection with novel and enzootic pathogens is increasingly likely. Novel and enzootic pathogens can interact synergistically or antagonistically, leading to increased or decreased disease severity. Here we examine host immune responses to coinfection with two closely related fungal pathogens: Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal). Both pathogens have had detrimental effects on amphibian populations, with Bd now largely enzootic, while Bsal is currently spreading and causing epizootics. Recent experimental work revealed that newts coinfected with Bd and Bsal had significantly higher mortality than those infected with either pathogen alone. Here we characterize host immunogenomic responses to chytrid coinfection relative to single infection. Across several classes of immune genes including pattern recognition receptors, cytokines, and MHC, coinfected host gene expression was weakly upregulated or comparable to that seen in single Bd infection, but significantly decreased when compared to Bsal infection. Combined with strong complement pathway downregulation and keratin upregulation, these results indicate that coinfection with Bd and Bsal compromises immune responses active against Bsal alone. As Bsal continues to invade naïve habitats where Bd is enzootic, coinfection will be increasingly common. If other Bd-susceptible species in the region have similar responses, interactions between the two pathogens could cause severe population and community-level declines.

Response to Comment on "Amphibian fungal panzootic causes catastrophic and ongoing loss of biodiversity".

Lambert et al question our retrospective and holistic epidemiological assessment of the role of chytridiomycosis in amphibian declines. Their alternative assessment is narrow and provides an incomplete evaluation of evidence. Adopting this approach limits understanding of infectious disease impacts and hampers conservation efforts. We reaffirm that our study provides unambiguous evidence that chytridiomycosis has affected at least 501 amphibian species.

Differentiating Batrachochytrium dendrobatidis and B. salamandrivorans in Amphibian Chytridiomycosis Using RNAScope® in situ Hybridization.

Batrachochytrium dendrobatidis and B. salamandrivorans are important amphibian pathogens responsible for morbidity and mortality in free-ranging and captive frogs, salamanders, and caecilians. While B. dendrobatidis has a widespread global distribution, B. salamandrivorans has only been detected in amphibians in Asia and Europe. Although molecular detection methods for these fungi are well-characterized, differentiation of the morphologically similar organisms in the tissues of affected amphibians is incredibly difficult. Moreover, an accurate tool to identify and differentiate Batrachochytrium in affected amphibian tissues is essential for a specific diagnosis of the causative agent in chytridiomycosis cases. To address this need, an automated dual-plex chromogenic RNAScope® in situ hybridization (ISH) assay was developed and characterized for simultaneous detection and differentiation of B. dendrobatidis and B. salamandrivorans. The assay, utilizing double Z target probe pairs designed to hybridize to 28S rRNA sequences, was specific for the identification of both organisms in culture and in formalin-fixed paraffin-embedded amphibian tissues. The assay successfully identified organisms in tissue samples from five salamander and one frog species preserved in formalin for up to 364 days and was sensitive for the detection of Batrachochytrium in animals with qPCR loads as low as 1.1 × 102 zoospores/microliter. ISH staining of B. salamandrivorans also highlighted the infection of dermal cutaneous glands, a feature not observed in amphibian B. dendrobatidis cases and which may play an important role in B. salamandrivorans pathogenesis in salamanders. The developed ISH assay will benefit both amphibian chytridiomycosis surveillance projects and pathogenesis studies by providing a reliable tool for Batrachochytrium differentiation in tissues.

Amphibian fungal panzootic causes catastrophic and ongoing loss of biodiversity.

Anthropogenic trade and development have broken down dispersal barriers, facilitating the spread of diseases that threaten Earth's biodiversity. We present a global, quantitative assessment of the amphibian chytridiomycosis panzootic, one of the most impactful examples of disease spread, and demonstrate its role in the decline of at least 501 amphibian species over the past half-century, including 90 presumed extinctions. The effects of chytridiomycosis have been greatest in large-bodied, range-restricted anurans in wet climates in the Americas and Australia. Declines peaked in the 1980s, and only 12% of declined species show signs of recovery, whereas 39% are experiencing ongoing decline. There is risk of further chytridiomycosis outbreaks in new areas. The chytridiomycosis panzootic represents the greatest recorded loss of biodiversity attributable to a disease.

Community richness of amphibian skin bacteria correlates with bioclimate at the global scale.

Animal-associated microbiomes are integral to host health, yet key biotic and abiotic factors that shape host-associated microbial communities at the global scale remain poorly understood. We investigated global patterns in amphibian skin bacterial communities, incorporating samples from 2,349 individuals representing 205 amphibian species across a broad biogeographic range. We analysed how biotic and abiotic factors correlate with skin microbial communities using multiple statistical approaches. Global amphibian skin bacterial richness was consistently correlated with temperature-associated factors. We found more diverse skin microbiomes in environments with colder winters and less stable thermal conditions compared with environments with warm winters and less annual temperature variation. We used bioinformatically predicted bacterial growth rates, dormancy genes and antibiotic synthesis genes, as well as inferred bacterial thermal growth optima to propose mechanistic hypotheses that may explain the observed patterns. We conclude that temporal and spatial characteristics of the host's macro-environment mediate microbial diversity.

Eco-evolutionary rescue promotes host-pathogen coexistence.

Emerging infectious pathogens are responsible for some of the most severe host mass mortality events in wild populations. Yet, effective pathogen control strategies are notoriously difficult to identify, in part because quantifying and forecasting pathogen spread and disease dynamics is challenging. Following an outbreak, hosts must cope with the presence of the pathogen, leading to host-pathogen coexistence or extirpation. Despite decades of research, little is known about host-pathogen coexistence post-outbreak when low host abundances and cryptic species make these interactions difficult to study. Using a novel disease-structured N-mixture model, we evaluate empirical support for three host-pathogen coexistence hypotheses (source-sink, eco-evolutionary rescue, and spatial variation in pathogen transmission) in a Neotropical amphibian community decimated by Batrachochytrium dendrobatidis (Bd) in 2004. During 2010-2014, we surveyed amphibians in Parque Nacional G. D. Omar Torríjos Herrera, Coclé Province, El Copé, Panama. We found that the primary driver of host-pathogen coexistence was eco-evolutionary rescue, as evidenced by similar amphibian survival and recruitment rates between infected and uninfected hosts. Average apparent monthly survival rates of uninfected and infected hosts were both close to 96%, and the expected number of uninfected and infected hosts recruited (via immigration/reproduction) was less than one host per disease state per 20-m site. The secondary driver of host-pathogen coexistence was spatial variation in pathogen transmission as we found that transmission was highest in areas of low abundance but there was no support for the source-sink hypothesis. Our results indicate that changes in the host community (i.e., through genetic or species composition) can reduce the impacts of emerging infectious disease post-outbreak. Our disease-structured N-mixture model represents a valuable advancement for conservation managers trying to understand underlying host-pathogen interactions and provides new opportunities to study disease dynamics in remnant host populations decimated by virulent pathogens.

External Reinfection of a Fungal Pathogen Does not Contribute to Pathogen Growth.

Chytridiomycosis is an emerging infectious disease of amphibians caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), which has led to devastating declines in amphibian populations worldwide. Current theory predicts that Bd infections are maintained through both reproduction on the host's skin and reinfection from sources outside of the host. To investigate the importance of external reinfection on pathogen burden, we infected captive-bred individuals of the highly susceptible Panamanian Golden Frog, Atelopus glyphus, and wild-caught glass frogs, Espadarana prosoblepon, with Bd. We housed the animals in one of three treatments: individually, in heterospecific pairs, and in conspecific pairs. For 8 weeks, we measured the Bd load and shedding rate of all frogs. We found that Atelopus had high rates of increase in both Bd load and shedding rate, but pathogen growth rates did not differ among treatments. The infection intensity of Espadarana co-housed with Atelopus was indistinguishable from those housed singly and those in conspecific pairs, despite being exposed to a large external source of Bd zoospores. Our results indicate that Bd load in both species is driven by pathogen replication within an individual, with reinfection from outside the host contributing little to the amplification of host fungal load.

Prevalence and genetic diversity of Batrachochytrium dendrobatidis in Central African island and continental amphibian communities.

The fungal pathogen Batrachochytrium dendrobatidis (Bd) infects hundreds of amphibian species and is implicated in global amphibian declines. Bd is comprised of several lineages that differ in pathogenicity, thus, identifying which Bd strains are present in a given amphibian community is essential for understanding host-pathogen dynamics. The presence of Bd has been confirmed in Central Africa, yet vast expanses of this region have not yet been surveyed for Bd prevalence, and the genetic diversity of Bd is largely unknown in this part of the world. Using retrospective surveys of museum specimens and contemporary field surveys, we estimated the prevalence of Bd in Central African island and continental amphibian assemblages, and genotyped strains of Bd present in each community. Our sampling of museum specimens included just a few individuals collected in the Gulf of Guinea archipelago prior to 1998, yet one of these individuals was Bd-positive indicating that the pathogen has been on Bioko Island since 1966. We detected Bd across all subsequent sample years in our study and found modest support for a relationship between host life history and Bd prevalence, a positive relationship between prevalence and host community species richness, and no significant relationship between elevation and prevalence. The Global Panzootic Lineage (Bd GPL) was present in all the island and continental amphibian communities we surveyed. Our results are consistent with a long-term and widespread distribution of Bd in amphibian communities of Gabon and the Gulf of Guinea archipelago.

Temperature variation, bacterial diversity and fungal infection dynamics in the amphibian skin.

Host-associated bacterial communities on the skin act as the first line of defence against invading pathogens. Yet, for most natural systems, we lack a clear understanding of how temperature variability affects structure and composition of skin bacterial communities and, in turn, promotes or limits the colonization of opportunistic pathogens. Here, we examine how natural temperature fluctuations might be related to changes in skin bacterial diversity over time in three amphibian populations infected by the pathogenic fungus Batrachochytrium dendrobatidis (Bd). Our focal host species (Eleutherodactylus coqui) is a direct-developing frog that has suffered declines at some populations in the last 20 years, while others have not experienced any changes. We quantified skin bacterial alpha- and beta-diversity at four sampling time points, a period encompassing two seasons and ample variation in natural infections and environmental conditions. Despite the different patterns of infection across populations, we detected an overall increase in bacterial diversity through time, characterized by the replacement of bacterial operational taxonomic units (OTUs). Increased frog body temperatures possibly allowed the colonization of bacteria as well as the recruitment of a subset of indicator OTUs, which could have promoted the observed changes in diversity patterns. Our results suggest that natural environmental fluctuations might be involved in creating opportunities for bacterial replacement, potentially attenuating pathogen transmission and thus contributing to host persistence in E. coqui populations.