Chytridiomycosis causes catastrophic organism-wide metabolic dysregulation including profound failure of cellular energy pathways

Chytridiomycosis disrupts metabolic responses in amphibians at metamorphic climax

The fungal disease chytridiomycosis (causative agent Batrachochytrium dendrobatidis [Bd]) is a primary contributor to amphibian species declines. The morphological and physiological reorganization that occurs during amphibian metamorphosis likely increases the vulnerability of metamorphs to Bd. To address this, we exposed pro-metamorphic tadpoles of Fleay's barred frog (Mixophyes fleayi) to Bd and sampled skin and liver sections from control and exposed animals throughout metamorphosis (Gosner stages 40, 42 and 45). We used an untargeted metabolomics approach to assess the metabolic impacts of Bd infection during the critical metamorphic stages, extracting metabolites from sampled tissues and analysing them via Nuclear Magnetic Resonance spectrometry. Most exposed animals became moribund at Gosner stage 45, while a subset seemingly cleared their infections. Metabolite abundance varied throughout development, with Gosner stage 45 samples distinct from previous stages. Clinically infected animals at Gosner stage 45 exhibited profound metabolic dysregulation (e.g., upregulation of amino acid biosynthesis and degradation) in comparison to uninfected groups (negative controls and 'cleared' animals). Despite showing parallels with previous metabolomic analyses of Bd-infected adult frogs, we identified variations in our results that could be attributed to the dramatic changes that characterise metamorphosis and may be driving the heightened vulnerability observed in metamorphic amphibians.

Batrachochytrium dendrobatidis strain affects transcriptomic response in liver but not skin in latitudinal populations of the common toad (Bufo bufo)

Batrachochytrium dendrobatidis (Bd) is a fungal pathogen that has decimated amphibian populations worldwide for several decades. We examined the changes in gene expression in response to Bd infection in two populations of the common toad, Bufo bufo, in a laboratory experiment. We collected B. bufo eggs in southern and northern Sweden, and infected the laboratory-raised metamorphs with two strains of the global panzoonotic lineage Bd-GPL. Differential expression analysis showed significant differences between infected and control individuals in both liver and skin. The skin samples showed no discernible differences in gene expression between the two strains used, while liver samples were differentiated by strain, with one of the strains eliciting no immune response from infected toads. Immune system genes were overexpressed in skin samples from surviving infected individuals, while in liver samples the pattern was more diffuse. Splitting samples by population revealed a stronger immune response in northern individuals. Differences in transcriptional regulation between populations are particularly relevant to study in Swedish amphibians, which may have experienced varying exposure to Bd. Earlier exposure to this pathogen and subsequent adaptation or selection pressure may contribute to the survival of some populations over others, while standing genetic diversity in different populations may also affect the infection outcome.

Toxicities of Legacy and Current-Use PFAS in an Anuran: Do Larval Exposures Influence Responses to a Terrestrial Pathogen Challenge?

Legacy polyfluoroalkyl substances (PFAS) [perfluorooctanesulfonate (PFOS) and perfluorooctanoic acid (PFOA)] are being replaced by various other fluorinated compounds, such as hexafluoropropylene oxide dimer acid (GenX). These alternatives are thought to be less bioaccumulative and, therefore, less toxic than legacy PFAS. Contaminant exposures occur concurrently with exposure to natural stressors, including the fungal pathogen Batrachocytrium dendrobatidis (Bd). Despite evidence that other pollutants can increase the adverse effects of Bd on anurans, no studies have examined the interactive effects of Bd and PFAS. This study tested the growth and developmental effects of PFOS, PFOA, and GenX on gray treefrog (Hyla versicolor) tadpoles, followed by a Bd challenge after metamorphosis. Despite PFAS exposure only occurring during the larval stage, carry-over effects on growth were observed post metamorphosis. Further, PFAS interacted with Bd exposure to influence growth; Bd-exposed animals had significantly shorter SVL [snout-vent length (mm)] with significantly increased body condition, among other time-dependent effects. Our data suggest that larval exposure to PFAS can continue to impact growth in the juvenile stage after exposure has ended. Contrary to predictions, GenX affected terrestrial performance more consistently than its legacy congener, PFOA. Given the role of Bd in amphibian declines, further investigation of interactions of PFAS with Bd and other environmentally relevant pathogens is warranted.

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.

Microbiome function predicts amphibian chytridiomycosis disease dynamics

BACKGROUND

The fungal pathogen Batrachochytrium dendrobatidis (Bd) threatens amphibian biodiversity and ecosystem stability worldwide. Amphibian skin microbial community structure has been linked to the clinical outcome of Bd infections, yet its overall functional importance is poorly understood.

METHODS

Microbiome taxonomic and functional profiles were assessed using high-throughput bacterial 16S rRNA and fungal ITS2 gene sequencing, bacterial shotgun metagenomics and skin mucosal metabolomics. We sampled 56 wild midwife toads (Alytes obstetricans) from montane populations exhibiting Bd epizootic or enzootic disease dynamics. In addition, to assess whether disease-specific microbiome profiles were linked to microbe-mediated protection or Bd-induced perturbation, we performed a laboratory Bd challenge experiment whereby 40 young adult A. obstetricans were exposed to Bd or a control sham infection. We measured temporal changes in the microbiome as well as functional profiles of Bd-exposed and control animals at peak infection.

RESULTS

Microbiome community structure and function differed in wild populations based on infection history and in experimental control versus Bd-exposed animals. Bd exposure in the laboratory resulted in dynamic changes in microbiome community structure and functional differences, with infection clearance in all but one infected animal. Sphingobacterium, Stenotrophomonas and an unclassified Commamonadaceae were associated with wild epizootic dynamics and also had reduced abundance in laboratory Bd-exposed animals that cleared infection, indicating a negative association with Bd resistance. This was further supported by microbe-metabolite integration which identified functionally relevant taxa driving disease outcome, of which Sphingobacterium and Bd were most influential in wild epizootic dynamics. The strong correlation between microbial taxonomic community composition and skin metabolome in the laboratory and field is inconsistent with microbial functional redundancy, indicating that differences in microbial taxonomy drive functional variation. Shotgun metagenomic analyses support these findings, with similar disease-associated patterns in beta diversity. Analysis of differentially abundant bacterial genes and pathways indicated that bacterial environmental sensing and Bd resource competition are likely to be important in driving infection outcomes.

CONCLUSIONS

Bd infection drives altered microbiome taxonomic and functional profiles across laboratory and field environments. Our application of multi-omics analyses in experimental and field settings robustly predicts Bd disease dynamics and identifies novel candidate biomarkers of infection. Video Abstract.

Ranavirus Amplification in Low-Diversity Amphibian Communities

In an era where emerging infectious diseases are a serious threat to biodiversity, epidemiological patterns need to be identified, particularly the complex mechanisms driving the dynamics of multi-host pathogens in natural communities. Many amphibian species have faced unprecedented population declines associated with diseases. Yet, specific processes shaping host-pathogen relationships within and among communities for amphibian pathogens such as ranaviruses (RV) remain poorly understood. To address this gap, we conducted a comprehensive study of RV in low-diversity amphibian communities in north-western Canada to assess the effects of biotic factors (species identity, species richness, abundance) and abiotic factors (conductivity, pH) on the pathogen prevalence and viral loads. Across 2 years and 18 sites, with communities of up to three hosts (wood frog, Rana sylvatica; boreal chorus frog, Pseudacris maculata; Canadian toad, Anaxyrus hemiophrys), we observed that RV prevalence nearly doubled with each additional species in a community, suggesting an amplification effect in aquatic, as well as terrestrial life-history stages. Infection intensity among infected wood frogs and boreal chorus frogs also significantly increased with an increase in species richness. Interestingly, we did not observe any effects of host abundance or abiotic factors, highlighting the importance of including host identity and species richness when investigating multi-host pathogens. Ultimately, only such a comprehensive approach can improve our understanding of complex and often highly context-dependent host-pathogen interactions.

Host gene expression in wildlife disease: making sense of species-level responses

Emerging infectious diseases are significant threats to wildlife conservation, yet the impacts of pathogen exposure and infection can vary widely among host species. As such, conservation biologists and disease ecologists have increasingly aimed to understand species-specific host susceptibility using molecular methods. In particular, comparative gene expression assays have been used to contrast the transcriptomic responses of disease-resistant and disease-susceptible hosts to pathogen exposure. This work usually assumes that the gene expression responses of disease-resistant species will reveal the activation of molecular pathways contributing to host defence. However, results often show that disease-resistant hosts undergo little gene expression change following pathogen challenge. Here, we discuss the mechanistic implications of these "null" findings and offer methodological suggestions for future molecular studies of wildlife disease. First, we highlight that muted transcriptomic responses with minimal immune system recruitment may indeed be protective for nonsusceptible hosts if they limit immunopathology and promote pathogen tolerance in systems where susceptible hosts suffer from genetic dysregulation. Second, we argue that overly narrow investigation of responses to pathogen exposure may overlook important, constitutively active molecular pathways that underlie species-specific defences. Finally, we outline alternative study designs and approaches that complement interspecific transcriptomic comparisons, including intraspecific gene expression studies and genomic methods to detect signatures of selection. Collectively, these insights will help ecologists extract maximal information from conservation-relevant transcriptomic data sets, leading to a deeper understanding of host defences and, ultimately, the implementation of successful conservation interventions.

Declining amphibians might be evolving increased reproductive effort in the face of devastating disease

The devastating infectious disease chytridiomycosis has caused declines of amphibians across the globe, yet some populations are persisting and even recovering. One understudied effect of wildlife disease is changes in reproductive effort. Here, we aimed to understand if the disease has plastic effects on reproduction and if reproductive effort could evolve with disease endemism. We compared the effects of experimental pathogen exposure (trait plasticity) and population-level disease history (evolution in trait baseline) on reproductive effort using gametogenesis as a proxy in the declining and endangered frog Litoria verreauxii alpina. We found that unexposed males from disease-endemic populations had higher reproductive effort, which is consistent with an evolutionary response to chytridiomycosis. We also found evidence of trait plasticity, where males and females were affected differently by infection: pathogen exposed males had higher reproductive effort (larger testes), whereas females had reduced reproductive effort (smaller and fewer developed eggs) regardless of the population of origin. Infectious diseases can cause plastic changes in the reproductive effort at an individual level, and population-level disease exposure can result in changes to baseline reproductive effort; therefore, individual- and population-level effects of disease should be considered when designing management and conservation programs for threatened and declining species.

Exposure to a fungal pathogen increases the critical thermal minimum of two frog species

The ability of an organism to tolerate seasonal temperature changes, such as extremely cold temperatures during the winter, can be influenced by their pathogens. We tested how exposure to a virulent fungal pathogen, Batrachochytrium dendrobatidis (Bd), affected the critical thermal minimum (CTmin) of two frog species, Hyla versicolor (gray treefrog) and Lithobates palustris (pickerel frog). The CTmin is the minimum thermal performance point of an organism, which we estimated via righting response trials. For both frog species, we compared the righting response of Bd-exposed and Bd-unexposed individuals in either a constant (15ºC) environment or with decreasing temperatures (-1°C/2.5 min) starting from 15°C. The CTmin for both species was higher for Bd-exposed frogs than unexposed frogs, and the CTmin of H. versicolor was higher than L. palustris. We also found that Bd-exposed frogs of both species righted themselves significantly fewer times in both decreasing and constant temperature trials. Our findings show that pathogen exposure can reduce cold tolerance and limit the thermal performance range of hosts, which may lead to increased overwintering mortality.

Differential liver histopathological responses to amphibian chytrid infection

Amphibians have been facing a pandemic caused by the deadly fungus Batrachochytrium dendrobatidis (Bd). Although studies have elucidated cutaneous and homeostatic disturbances, it is still unknown if the hepatic function can be affected or if hepatic effects differ among host species. Thus, we evaluated the effects of an experimental Bd infection on the liver (histopathology and the hepatosomatic index) of 2 anuran species (Xenopus laevis and Physalaemus albonotatus) with different susceptibilities to Bd infection and compared them to uninfected controls. Bd infection increased the melanomacrophage cell area and induced leukocyte infiltration in both species. The effects were more pronounced in the sensitive species, P. albonotatus, which showed severe reduction in glycogen stores and liver atrophy, due to energetic imbalance. Hepatocytes of P. albonotatus also showed ballooning degeneration (vacuolization), which could lead to cell death and liver failure. Our results provide evidence that although the sensitive species showed more severe effects, the tolerant species also had hepatic responses to the infection. These findings indicate that hepatic function can play an important role in detoxification and in immune responses to chytridiomycosis, and that it may be used as a new biomarker of health status in chytrid infections.

Immunological Aspects of Chytridiomycosis

Amphibians are currently the most threatened vertebrate class, with the disease chytridiomycosis being a major contributor to their global declines. Chytridiomycosis is a frequently fatal skin disease caused by the fungal pathogens Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal). The severity and extent of the impact of the infection caused by these pathogens across modern Amphibia are unprecedented in the history of vertebrate infectious diseases. The immune system of amphibians is thought to be largely similar to that of other jawed vertebrates, such as mammals. However, amphibian hosts are both ectothermic and water-dependent, which are characteristics favouring fungal proliferation. Although amphibians possess robust constitutive host defences, Bd/Bsal replicate within host cells once these defences have been breached. Intracellular fungal localisation may contribute to evasion of the induced innate immune response. Increasing evidence suggests that once the innate defences are surpassed, fungal virulence factors suppress the targeted adaptive immune responses whilst promoting an ineffectual inflammatory cascade, resulting in immunopathology and systemic metabolic disruption. Thus, although infections are contained within the integument, crucial homeostatic processes become compromised, leading to mortality. In this paper, we present an integrated synthesis of amphibian post-metamorphic immunological responses and the corresponding outcomes of infection with Bd, focusing on recent developments within the field and highlighting future directions.

A meta-analysis reveals temperature, dose, life stage, and taxonomy influence host susceptibility to a fungal parasite

Complex ecological relationships, such as host-parasite interactions, are often modeled with laboratory experiments. However, some experimental laboratory conditions, such as temperature or infection dose, are regularly chosen based on convenience or convention, and it is unclear how these decisions systematically affect experimental outcomes. Here, we conducted a meta-analysis of 58 laboratory studies that exposed amphibians to the pathogenic fungus Batrachochytrium dendrobatidis (Bd) to understand better how laboratory temperature, host life stage, infection dose, and host species affect host mortality. We found that host mortality was driven by thermal mismatches: hosts native to cooler environments experienced greater Bd-induced mortality at relatively warm experimental temperatures and vice versa. We also found that Bd dose positively predicted Bd-induced host mortality and that the superfamilies Bufonoidea and Hyloidea were especially susceptible to Bd. Finally, the effect of Bd on host mortality varied across host life stages, with larval amphibians experiencing lower risk of Bd-induced mortality than adults or metamorphs. Metamorphs were especially susceptible and experienced mortality when inoculated with much smaller Bd doses than the average dose used by researchers. Our results suggest that when designing experiments on species interactions, researchers should carefully consider the experimental temperature, inoculum dose, and life stage, and taxonomy of the host species.

In vitro modeling of Batrachochytrium dendrobatidis infection of the amphibian skin

The largest current disease-induced loss of vertebrate biodiversity is due to chytridiomycosis and despite the increasing understanding of the pathogenesis, knowledge unravelling the early host-pathogen interactions remains limited. Batrachochytrium dendrobatidis (Bd) zoospores attach to and invade the amphibian epidermis, with subsequent invasive growth in the host skin. Availability of an in vitro assay would facilitate in depth study of this interaction while reducing the number of experimental animals needed. We describe a fluorescent cell-based in vitro infection model that reproduces host-Bd interactions. Using primary keratinocytes from Litoria caerulea and the epithelial cell line A6 from Xenopus laevis, we reproduced different stages of host cell infection and intracellular growth of Bd, resulting in host cell death, a key event in chytridiomycosis. The presented in vitro models may facilitate future mechanistic studies of host susceptibility and pathogen virulence.

Susceptibility of frogs to chytridiomycosis correlates with increased levels of immunomodulatory serotonin in the skin

Chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), is a skin disease responsible for the global decline of amphibians. Frog species and populations can vary in susceptibility, but this phenomenon remains poorly understood. Here, we investigated serotonin in the skin of infected and uninfected frogs. In more susceptible frog populations, skin serotonin rose with increasing infection intensity, but decreased in later stages of the disease. The more resistant population maintained a basal level of skin serotonin. Serotonin inhibited both Bd sporangial growth and Jurkat lymphocyte proliferation in vitro. However, serotonin accumulates in skin granular glands, and this compartmentalisation may prevent inhibition of Bd growth in vivo. We suggest that skin serotonin increases in susceptible frogs due to pathogen excretion of precursor tryptophan, but that resistant frogs are able to control the levels of serotonin. Overall, the immunosuppressive effects of serotonin may contribute to the susceptibility of frogs to chytridiomycosis.

Review of the Amphibian Immune Response to Chytridiomycosis, and Future Directions

The fungal skin disease, chytridiomycosis (caused by Batrachochytrium dendrobatidis and B. salamandrivorans), has caused amphibian declines and extinctions globally since its emergence. Characterizing the host immune response to chytridiomycosis has been a focus of study with the aim of disease mitigation. However, many aspects of the innate and adaptive arms of this response are still poorly understood, likely due to the wide range of species' responses to infection. In this paper we provide an overview of expected immunological responses (with inference based on amphibian and mammalian immunology), together with a synthesis of current knowledge about these responses for the amphibian-chytridiomycosis system. We structure our review around four key immune stages: (1) the naïve immunocompetent state, (2) immune defenses that are always present (constitutive defenses), (3) mechanisms for recognition of a pathogen threat and innate immune defenses, and (4) adaptive immune responses. We also evaluate the current hot topics of immunosuppression and immunopathology in chytridiomycosis, and discuss their respective roles in pathogenesis. Our synthesis reveals that susceptibility to chytridiomycosis is likely to be multifactorial. Susceptible amphibians appear to have ineffective constitutive and innate defenses, and a late-stage response characterized by immunopathology and Bd-induced suppression of lymphocyte responses. Overall, we identify substantial gaps in current knowledge, particularly concerning the entire innate immune response (mechanisms of initial pathogen detection and possible immunoevasion by Bd, degree of activation and efficacy of the innate immune response, the unexpected absence of innate leukocyte infiltration, and the cause and role of late-stage immunopathology in pathogenesis). There are also gaps concerning most of the adaptive immune system (the relative importance of B and T cell responses for pathogen clearance, the capacity and extent of immunological memory, and specific mechanisms of pathogen-induced immunosuppression). Improving our capacity for amphibian immunological research will require selection of an appropriate Bd-susceptible model species, the development of taxon-specific affinity reagents and cell lines for functional assays, and the application of a suite of conventional and emerging immunological methods. Despite current knowledge gaps, immunological research remains a promising avenue for amphibian conservation management.