Evolution of resistance to chytridiomycosis is associated with a robust early immune response

Amphibian (Xenopus laevis) Macrophage Subsets Vary in Their Responses to the Chytrid Fungus Batrachochytrium dendrobatidis

The chytrid fungus, Batrachochytrium dendrobatidis (Bd), infects amphibian skin, causing chytridiomycosis, which is a contributing cause of worldwide declines and extinctions of amphibians. Relatively little is known about the roles of amphibian skin-resident immune cells, such as macrophages, in these antifungal defenses. Across vertebrates, macrophage differentiation is controlled through the activation of colony-stimulating factor-1 (CSF1) receptor by CSF1 and interleukin-34 (IL34) cytokines. While the precise roles of these respective cytokines in macrophage development remain to be fully explored, our ongoing studies indicate that frog (Xenopus laevis) macrophages differentiated by recombinant forms of CSF1 and IL34 are functionally distinct. Accordingly, we explored the roles of X. laevis CSF1- and IL34-macrophages in anti-Bd defenses. Enriching cutaneous IL34-macrophages, but not CSF1-macrophages, resulted in significant anti-Bd protection. In vitro analysis of frog macrophage-Bd interactions indicated that both macrophage subsets phagocytosed Bd. However, IL34-macrophages cocultured with Bd exhibited greater pro-inflammatory gene expression, whereas CSF1-macrophages cocultured with Bd showed greater immunosuppressive gene expression profiles. Concurrently, Bd-cocultured with CSF1-macrophages, but not IL34-macrophages, possessed elevated expression of genes associated with immune evasion. This work marks a step forward in our understanding of the roles of frog macrophage subsets in antifungal defenses.

Amphibian cellular immune response to chytridiomycosis at metamorphic climax

The fungal disease chytridiomycosis (caused by Batrachochytrium dendrobatidis [Bd]) is a primary contributor to amphibian declines. The frog metamorphic stages, characterised by extensive physiological reorganisation and energy expenditure, have heightened susceptibility to Bd. However, little is known about how these metamorphic stages respond immunologically to Bd infection. In this study, we examined Bd infection and the cellular immune response of Mixophyes fleayi at Gosner stages 40, 42 and 45, using blood smears and skin and liver histology. Although proportional differences were observed, the impact of Bd exposure appeared negligible prior to Gosner stage 45 (onset of morbidity), with no significant differences observed in absolute leukocyte counts for blood or liver samples between control and Bd-exposed groups at Gosner stages 40 and 42. Animals exhibiting clinical signs at Gosner stage 45 demonstrated significant elevation in liver leukocyte counts, blood neutrophil and monocyte counts and neutrophil-to-lymphocyte ratios. These findings are reminiscent of the amplified inflammatory response characteristic of immunopathology in clinically infected amphibians. Interestingly, a subset of exposed animals that had apparently cleared infections at Gosner stage 45 had similar blood leukocyte counts but reduced liver leukocyte counts compared to naïve controls. This could be a consequence of prior cellular consumption during pathogen removal or effective immune regulation via anti-inflammatory protective feedback mechanisms. We recommend targeted gene expression analyses (e.g. immunomodulatory cytokines) to establish the mechanisms responsible for the varied immune expression and infection outcomes across metamorphosis.

Chytridiomycosis in a colony of hellbenders Cryptobranchus alleganiensis

OBJECTIVE

The hellbender Cryptobranchus alleganiensis is a fully aquatic, long-lived, and state endangered salamander endemic to the eastern and central United States. Chytridiomycosis is a fungal skin disease of salamanders that can be caused by either Batrachochytrium dendrobatidis (Bd) or by recently described B. salamandrivorans. B. salamandrivorans is responsible for massive mortality of salamanders in Europe but, to date, has not been reported in the United States, in contrast to Bd. Common signs of Bd infection in salamanders are excessive skin shedding; cutaneous discoloration; mainly in ventral parts of the abdomen and thighs; and eventually death as a result of osmotic imbalance.

METHODS

Eastern hellbenders C. a. alleganiensis in the Aquaculture Research Laboratory at Purdue University experienced 33% mortality. Three of the deceased animals were autopsied, and samples for histopathology, microbiology, and Taqman quantitative polymerase chain reaction were taken.

RESULT

Autopsy revealed an abundant, thick, yellow to green mucoid material on the surface of the skin, mainly on the tail and flanks. Histopathology of the skin samples revealed infection with chytrid fungus consistent with chytridiomycosis. Taqman quantitative polymerase chain reaction confirmed presence of Bd, and bacterial culture confirmed co-infection with Aeromonas hydrophila. Infected hellbenders were treated with ciprofloxacin, itraconazole, and ceftazidime.

CONCLUSION

The chytridiomycosis outbreak was most likely caused by sudden environmental stress due to filtration failure and secondary bacterial infection. The surviving animals have not shown any clinical signs of chytridiomycosis for more than 1 year after the treatment and were released the wild river system.

Reintroduction of resistant frogs facilitates landscape-scale recovery in the presence of a lethal fungal disease

Vast alteration of the biosphere by humans is causing a sixth mass extinction, driven in part by an increase in infectious diseases. The emergence of the lethal fungal pathogen Batrachochytrium dendrobatidis (Bd) has devastated global amphibian biodiversity. Given the lack of any broadly applicable methods to reverse these impacts, the future of many amphibians appears grim. The Sierra Nevada yellow-legged frog (Rana sierrae) is highly susceptible to Bd infection and most R. sierrae populations are extirpated following disease outbreaks. However, some populations persist and eventually recover, and frogs in these recovering populations have increased resistance against infection. Here, we conduct a 15-year reintroduction study and show that frogs collected from recovering populations and reintroduced to vacant habitats can reestablish populations despite the presence of Bd. In addition, the likelihood of establishment is influenced by site, cohort, and frog attributes. Results from viability modeling suggest that many reintroduced populations have a low probability of extinction over 50 years. These results provide a rare example of how reintroduction of resistant individuals can allow the landscape-scale recovery of disease-impacted species, and have broad implications for amphibians and other taxa that are threatened with extinction by novel pathogens.

What Can Frogs Teach Us about Resilience? Adaptive Renewal in Amphibian and Academic Ecosystems

Examples of resilience in nature give us hope amid a growing biodiversity crisis. While resilience has many definitions across disciplines, here I discuss resilience as the ability to continue to adapt and persist. Naturally, as biologists, we seek to uncover the underlying mechanisms that can help us explain the secrets of resilience across scales, from individuals to species to ecosystems and beyond. Perhaps we also ponder what the secrets to resilience are in our own lives, in our own research practices, and academic communities. In this paper, I highlight insights gained through studies of amphibian resilience following a global disease outbreak to uncover shared patterns and processes linked to resilience across amphibian communities. I also reflect on how classical resilience heuristics could be more broadly applied to these processes and to our own academic communities. Focusing on the amphibian systems that I have worked in-the Golden Frogs of Panama (Atelopus zeteki/varius) and the Mountain Yellow-Legged Frogs of California (Rana muscosa/sierrae)-I highlight shared and unique characteristics of resilience across scales and systems and discuss how these relate to adaptive renewal cycles. Reflecting on this work and previous resilience scholarship, I also offer my own thoughts about academia and consider what lessons we could take from mapping our own adaptive trajectories and addressing threats to our own community resilience.

Advances in Managing Chytridiomycosis for Australian Frogs: Gradarius Firmus Victoria

Extensive knowledge gains from research worldwide over the 25 years since the discovery of chytridiomycosis can be used for improved management. Strategies that have saved populations in the short term and/or enabled recovery include captive breeding, translocation into disease refugia, translocation from resistant populations, disease-free exclosures, and preservation of disease refuges with connectivity to previous habitat, while antifungal treatments have reduced mortality rates in the wild. Increasing host resistance is the goal of many strategies under development, including vaccination and targeted genetic interventions. Pathogen-directed strategies may be more challenging but would have broad applicability. While the search for the silver bullet solution continues, we should value targeted local interventions that stop extinction and buy time for evolution of resistance or development of novel solutions. As for most invasive species and infectious diseases, we need to accept that ongoing management is necessary. For species continuing to decline, proactive deployment and assessment of promising interventions are more valid than a hands-off, do-no-harm approach that will likely allow further extinctions.

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.

Genetic mechanisms and biological processes underlying host response to ophidiomycosis (snake fungal disease) inferred from tissue-specific transcriptome analyses

Emerging infectious diseases in wildlife species caused by pathogenic fungi are of growing concern, yet crucial knowledge gaps remain for diseases with potentially large impacts. For example, there is detailed knowledge about host pathology and mechanisms underlying response for chytridiomycosis in amphibians and white-nose syndrome in bats, but such information is lacking for other more recently described fungal infections. One such disease is ophidiomycosis, caused by the fungus Ophidiomyces ophidiicola, which has been identified in many species of snakes, yet the biological mechanisms and molecular changes occurring during infection are unknown. To gain this information, we performed a controlled experimental infection in captive Prairie rattlesnakes (Crotalus viridis) with O. ophidiicola at two different temperatures: 20 and 26°C. We then compared liver, kidney, and skin transcriptomes to assess tissue-specific genetic responses to O. ophidiicola infection. Given previous histopathological studies and the fact that snakes are ectotherms, we expected highest fungal activity on skin and a significant impact of temperature on host response. Although we found fungal activity to be localized on skin, most of the differential gene expression occurred in internal tissues. Infected snakes at the lower temperature had the highest host mortality whereas two-thirds of the infected snakes at the higher temperature survived. Our results suggest that ophidiomycosis is likely a systemic disease with long-term effects on host response. Our analysis also identified candidate protein coding genes that are potentially involved in host response, providing genetic tools for studies of host response to ophidiomycosis in natural populations.

Gene expression reveals immune response strategies of naïve Hawaiian honeycreepers experimentally infected with introduced avian malaria

The unprecedented rise in the number of new and emerging infectious diseases in the last quarter century poses direct threats to human and wildlife health. The introduction to the Hawaiian archipelago of Plasmodium relictum and the mosquito vector that transmits the parasite has led to dramatic losses in endemic Hawaiian forest bird species. Understanding how mechanisms of disease immunity to avian malaria may evolve is critical as climate change facilitates increased disease transmission to high elevation habitats where malaria transmission has historically been low and the majority of the remaining extant Hawaiian forest bird species now reside. Here, we compare the transcriptomic profiles of highly susceptible Hawai'i 'amakihi (Chlorodrepanis virens) experimentally infected with P. relictum to those of uninfected control birds from a naïve high elevation population. We examined changes in gene expression profiles at different stages of infection to provide an in-depth characterization of the molecular pathways contributing to survival or mortality in these birds. We show that the timing and magnitude of the innate and adaptive immune response differed substantially between individuals that survived and those that succumbed to infection, and likely contributed to the observed variation in survival. These results lay the foundation for developing gene-based conservation strategies for Hawaiian honeycreepers by identifying candidate genes and cellular pathways involved in the pathogen response that correlate with a bird's ability to recover from malaria infection.

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

Molecular diversity and functional implication of amphibian interferon complex: Remarking immune adaptation in vertebrate evolution

Cross-species comparison of vertebrate genomes has unraveled previously unknown complexities of interferon (IFN) systems in amphibian species. Recent genomic curation revealed that amphibian species have evolved expanded repertoires of four types of intron-containing IFN genes akin to those seen in jawed fish, intronless type I IFNs and intron-containing type III IFNs akin to those seen in amniotes, as well as uniquely intronless type III IFNs. This appears to be the case with at least ten analyzed amphibian species; with distinct species encoding diverse repertoires of these respective IFN gene subsets. Amphibians represent a key stage in vertebrate evolution, and in this context offer a unique perspective into the divergent and converged pathways leading to the emergence of distinct IFN families and groups. Recent studies have begun to unravel the roles of amphibian IFNs during these animals' immune responses in general and during their antiviral responses, in particular. However, the pleiotropic potentials of these highly expanded amphibian IFN repertoires warrant further studies. Based on recent reports and our omics analyses using Xenopus models, we posit that amphibian IFN complex may have evolved novel functions, as indicated by their extensive molecular diversity. Here, we provide an overview and an update of the present understanding of the amphibian IFN complex in the context of the evolution of vertebrate immune systems. A greater understanding of the amphibian IFN complex will grant new perspectives on the evolution of vertebrate immunity and may yield new measures by which to counteract the global amphibian declines.

Host immune responses to enzootic and invasive pathogen lineages vary in magnitude, timing, and efficacy

Infectious diseases of wildlife continue to pose a threat to biodiversity worldwide, yet pathogens are far from uniform in virulence or host disease outcome. Within the same pathogen species, virulence can vary considerably depending on strain or lineage, in turn eliciting variable host responses. One pathogen that has caused extensive biodiversity loss is the amphibian-killing fungus, Batrachochytrium dendrobatidis (Bd), which is comprised of a globally widespread hypervirulent lineage (Bd-GPL), and multiple geographically restricted, enzootic lineages. Whereas host immunogenomic responses to Bd-GPL have been characterized in a number of amphibian species, immunogenomic responses to geographically restricted, enzootic Bd lineages are less clear. To examine lineage-specific host immune responses to Bd, we exposed a species of pumpkin toadlet, Brachycephalus pitanga, which is endemic to Brazil's Southern Atlantic Forest, to either the Bd-GPL or the enzootic Bd-Asia-2/Brazil (hereafter Bd-Brazil) lineage. Using temporal samples from early, mid, and late infection stages, we quantified functional immunogenomic responses over the course of infection using differential gene expression tests and coexpression network analyses. Host immune responses varied significantly with Bd lineage. Relative to controls, toadlet responses to Bd-Brazil were weak at early infection (25 genes significantly differentially expressed), peaked by mid-stage infection (414 genes), and were nearly fully resolved by late-stage infection (nine genes). In contrast, responses to Bd-GPL were magnified and delayed; toadlets significantly differentially expressed 111 genes early, 87 genes at mid-stage infection, and 726 genes by late-stage infection relative to controls. Given that infection intensity did not vary between mid- and late-stage disease in either Bd-Brazil or Bd-GPL treatments, this suggests that pumpkin toadlets may be at least partially tolerant to the enzootic Bd-Brazil lineage. In contrast, late-stage immune activation against Bd-GPL was consistent with immune dysregulation previously observed in other species. Our results demonstrate that both the timing of immune response and the particular immune pathways activated are specific to Bd lineage. Within regions where multiple Bd lineages co-occur, and given continued global Bd movement, these differential host responses may influence not only individual disease outcome, but transmission dynamics at the population and community levels.

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.

Do immune system changes at metamorphosis predict vulnerability to chytridiomycosis? An update

Amphibians are among the vertebrate groups suffering great losses of biodiversity due to a variety of causes including diseases, such as chytridiomycosis (caused by the fungal pathogens Batrachochytrium dendrobatidis and B. salamandrivorans). The amphibian metamorphic period has been identified as being particularly vulnerable to chytridiomycosis, with dramatic physiological and immunological reorganisation likely contributing to this vulnerability. Here, we overview the processes behind these changes at metamorphosis and then perform a systematic literature review to capture the breadth of empirical research performed over the last two decades on the metamorphic immune response. We found that few studies focused specifically on the immune response during the peri-metamorphic stages of amphibian development and fewer still on the implications of their findings with respect to chytridiomycosis. We recommend future studies consider components of the immune system that are currently under-represented in the literature on amphibian metamorphosis, particularly pathogen recognition pathways. Although logistically challenging, we suggest varying the timing of exposure to Bd across metamorphosis to examine the relative importance of pathogen evasion, suppression or dysregulation of the immune system. We also suggest elucidating the underlying mechanisms of the increased susceptibility to chytridiomycosis at metamorphosis and the associated implications for population persistence. For species that overlap a distribution where Bd/Bsal are now endemic, we recommend a greater focus on management strategies that consider the important peri-metamorphic period.

Whole genome resequencing reveals signatures of rapid selection in a virus-affected commercial fishery

Infectious diseases are recognized as one of the greatest global threats to biodiversity and ecosystem functioning. Consequently, there is a growing urgency to understand the speed at which adaptive phenotypes can evolve and spread in natural populations to inform future management. Here we provide evidence of rapid genomic changes in wild Australian blacklip abalone (Haliotis rubra) following a major population crash associated with an infectious disease. Genome scans on H. rubra were performed using pooled whole genome resequencing data from commercial fishing stocks varying in historical exposure to haliotid herpesvirus-1 (HaHV-1). Approximately 25,000 single nucleotide polymorphism loci associated with virus exposure were identified, many of which mapped to genes known to contribute to HaHV-1 immunity in the New Zealand pāua (Haliotis iris) and herpesvirus response pathways in haliotids and other animal systems. These findings indicate genetic changes across a single generation in H. rubra fishing stocks decimated by HaHV-1, with stock recovery potentially determined by rapid evolutionary changes leading to virus resistance. This is a novel example of apparently rapid adaptation in natural populations of a nonmodel marine organism, highlighting the pace at which selection can potentially act to counter disease in wildlife communities.

Effectiveness of antifungal treatments during chytridiomycosis epizootics in populations of an endangered frog

The recently-emerged amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd) has had an unprecedented impact on global amphibian populations, and highlights the urgent need to develop effective mitigation strategies. We conducted in-situ antifungal treatment experiments in wild populations of the endangered mountain yellow-legged frog during or immediately after Bd-caused mass die-off events. The objective of treatments was to reduce Bd infection intensity ("load") and in doing so alter frog-Bd dynamics and increase the probability of frog population persistence despite ongoing Bd infection. Experiments included treatment of early life stages (tadpoles and subadults) with the antifungal drug itraconazole, treatment of adults with itraconazole, and augmentation of the skin microbiome of subadults with Janthinobacterium lividum, a commensal bacterium with antifungal properties. All itraconazole treatments caused immediate reductions in Bd load, and produced longer-term effects that differed between life stages. In experiments focused on early life stages, Bd load was reduced in the 2 months immediately following treatment and was associated with increased survival of subadults. However, Bd load and frog survival returned to pre-treatment levels in less than 1 year, and treatment had no effect on population persistence. In adults, treatment reduced Bd load and increased frog survival over the entire 3-year post-treatment period, consistent with frogs having developed an effective adaptive immune response against Bd. Despite this protracted period of reduced impacts of Bd on adults, recruitment into the adult population was limited and the population eventually declined to near-extirpation. In the microbiome augmentation experiment, exposure of subadults to a solution of J. lividum increased concentrations of this potentially protective bacterium on frogs. However, concentrations declined to baseline levels within 1 month and did not have a protective effect against Bd infection. Collectively, these results indicate that our mitigation efforts were ineffective in causing long-term changes in frog-Bd dynamics and increasing population persistence, due largely to the inability of early life stages to mount an effective immune response against Bd. This results in repeated recruitment failure and a low probability of population persistence in the face of ongoing Bd infection.

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.

Invasion history shapes host transcriptomic response to a body-snatching parasite

By shuffling biogeographical distributions, biological invasions can both disrupt long-standing associations between hosts and parasites and establish new ones. This creates natural experiments with which to study the ecology and evolution of host-parasite interactions. In estuaries of the Gulf of Mexico, the white-fingered mud crab (Rhithropanopeus harrisii) is infected by a native parasitic barnacle, Loxothylacus panopaei (Rhizocephala), which manipulates host physiology and behaviour. In the 1960s, L. panopaei was introduced to the Chesapeake Bay and has since expanded along the southeastern Atlantic coast, while host populations in the northeast have so far been spared. We use this system to test the host's transcriptomic response to parasitic infection and investigate how this response varies with the parasite's invasion history, comparing populations representing (i) long-term sympatry between host and parasite, (ii) new associations where the parasite has invaded during the last 60 years and (iii) naïve hosts without prior exposure. A comparison of parasitized and control crabs revealed a core response, with widespread downregulation of transcripts involved in immunity and moulting. The transcriptional response differed between hosts from the parasite's native range and where it is absent, consistent with previous observations of increased susceptibility in populations lacking exposure to the parasite. Crabs from the parasite's introduced range, where prevalence is highest, displayed the most dissimilar response, possibly reflecting immune priming. These results provide molecular evidence for parasitic manipulation of host phenotype and the role of gene regulation in mediating host-parasite interactions.

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.

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.

Ecoimmunology at spatial scales

IN FOCUS

Becker, D. J., Albery, G. F., Kessler, M. K., Lunn, T. J., Falvo, C. A., Czirják, G. Á., Martin, L. B., & Plowright, R. K. (2020). Macroimmunology: The drivers and consequences of spatial patterns in wildlife immune defence. Journal of Animal Ecology, 89, 972-995. Ecoimmunology seeks to identify and explain natural variation in immune function. Most research so far has focused on differences among individuals within populations, which are often driven by trade-offs in resource allocation between energetically costly immunity and competing processes such as reproduction. In their review article, Becker et al. (2020) have proposed a framework to explicitly address habitat- and population-level differences in wildlife immune phenotypes. Termed macroimmunology, this concept integrates principles from ecoimmunology and macroecology. Becker et al. (2020) have highlighted three non-mutually exclusive habitat features that are likely to vary at spatial scales and influence immune function: (a) parasite pressure, (b) abiotic and biotic factors and (c) anthropogenic changes. However, a large and robust body of literature suitable for synthesis to detect macroimmunology patterns and effect sizes is not yet available. Through their systematic review and critical assessment, Becker et al. (2020) identified common problems in existing research that hinders spatial inferences, such as a need for spatial replication in study design and statistical analyses that account for spatial dependence. Overall, macroimmunology has the potential to identify and even predict spatial patterns in immune phenotypes that form the mechanistic underpinnings of important wildlife disease processes, and this review represents an important step to realizing these goals.

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.

Mechanisms underlying host persistence following amphibian disease emergence determine appropriate management strategies

Emerging infectious diseases have caused many species declines, changes in communities and even extinctions. There are also many species that persist following devastating declines due to disease. The broad mechanisms that enable host persistence following declines include evolution of resistance or tolerance, changes in immunity and behaviour, compensatory recruitment, pathogen attenuation, environmental refugia, density-dependent transmission and changes in community composition. Here we examine the case of chytridiomycosis, the most important wildlife disease of the past century. We review the full breadth of mechanisms allowing host persistence, and synthesise research on host, pathogen, environmental and community factors driving persistence following chytridiomycosis-related declines and overview the current evidence and the information required to support each mechanism. We found that for most species the mechanisms facilitating persistence have not been identified. We illustrate how the mechanisms that drive long-term host population dynamics determine the most effective conservation management strategies. Therefore, understanding mechanisms of host persistence is important because many species continue to be threatened by disease, some of which will require intervention. The conceptual framework we describe is broadly applicable to other novel disease systems.

The amphibian complement system and chytridiomycosis

Understanding host immune function and ecoimmunology is increasingly important at a time when emerging infectious diseases (EIDs) threaten wildlife. One EID that has emerged and spread widely in recent years is chytridiomycosis, caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), which is implicated unprecedented amphibian declines around the world. The impacts of Bd have been severe for many amphibian species, but some populations have exhibited signs of persistence, and even recovery, in some regions. Many mechanisms may underpin this pattern and amphibian immune responses are likely one key component. Although we have made great strides in understanding amphibian immunity, the complement system remains poorly understood. The complement system is a nonspecific, innate immune defense that is known to enhance other immune responses. Complement activation can occur by three different biochemical pathways and result in protective mechanisms, such as inflammation, opsonization, and pathogen lysis, thereby providing protection to the host. We currently lack an understanding of complement pathway activation for chytridiomycosis, but several studies have suggested that it may be a key part of an early and robust immune response that confers host resistance. Here, we review the available research on the complement system in general as well as amphibian complement responses to Bd infection. Additionally, we propose future research directions that will increase our understanding of the amphibian complement system and other immune responses to Bd. Finally, we suggest how a deeper understanding of amphibian immunity could enhance the conservation and management of amphibian species that are threatened by chytridiomycosis.

Habitat disturbance influences the skin microbiome of a rediscovered neotropical-montane frog

BACKGROUND

The skin microbiome serves as a first line defense against pathogens in vertebrates. In amphibians, it has the potential to protect against the chytrid fungus Batrachochytrium dendrobatis (Bd), a likely agent of amphibian declines. Alteration of the microbiome associated with unfavorable environmental changes produced by anthropogenic activities may make the host more susceptible to pathogens. Some amphibian species that were thought to be "extinct" have been rediscovered years after population declines in the late 1980s probably due to evolved Bd-resistance and are now threatened by anthropogenic land-use changes. Understanding the effects of habitat disturbance on the host skin microbiome is relevant for understanding the health of these species, along with its susceptibility to pathogens such as Bd. Here, we investigate the influence of habitat alteration on the skin bacterial communities as well as specifically the putative Bd-inhibitory bacterial communities of the montane frog Lithobates vibicarius. This species, after years of not being observed, was rediscovered in small populations inhabiting undisturbed and disturbed landscapes, and with continuous presence of Bd.

RESULTS

We found that cutaneous bacterial communities of tadpoles and adults differed between undisturbed and disturbed habitats. The adults from disturbed habitats exhibited greater community dispersion than those from undisturbed habitats. We observed a higher richness of putative Bd-inhibitory bacterial strains in adults from disturbed habitats than in those from undisturbed habitats, as well as a greater number of these potential protective bacteria with a high relative abundance.

CONCLUSIONS

Our findings support the microbial "Anna Karenina principle", in which disturbance is hypothesized to cause greater microbial dispersion in communities, a so-called dysbiosis, which is a response of animal microbiomes to stress factors that decrease the ability of the host or its microbiome to regulate community composition. On the positive side, the high richness and relative abundance of putative Bd-inhibitory bacteria may indicate the development of a defense mechanism that enhances Bd-protection, attributed to a co-occurrence of more than 30-years of host and pathogen in these disturbed habitats. Our results provide important insight into the influence of human-modified landscapes on the skin microbiome and health implications of Bd-survivor species.

Sustained immune activation is associated with susceptibility to the amphibian chytrid fungus

The disease chytridiomycosis caused by the fungus Bd has devastated amphibian populations worldwide. Functional genomic contributions to host susceptibility remain enigmatic and vary between species and populations. We conducted experimental Bd infections in Rana yavapaiensis, a species with intraspecific variation in chytridiomycosis susceptibility, to assess the skin and spleen transcriptomic response to infection over time. We predicted that increased immune gene expression would be associated with a positive disease outcome, but we instead found that surviving frogs had significantly reduced immune gene expression compared to susceptible frogs and to uninfected controls. MHC class IIβ gene expression was also significantly higher in susceptible frogs compared to surviving frogs. Furthermore, susceptible frogs expressed a significantly larger number of distinct class IIβ alleles, demonstrating a negative correlation between class IIβ expression, functional diversity, and survival. Expression of the MHC class IIβ locus previously associated with Bd disease outcomes was a significant predictor of Bd infection intensity at early infection stages but not at late infection stages, suggesting initial MHC-linked immune processes are important for ultimate disease outcomes. We infer through disease association and phylogenetic analysis that certain MHC variants are linked to the immune expression that was negatively associated with survival, and we hypothesize that frogs that did not express these alleles could better survive infections. Our study finds that MHC expression at early and late infection stages predicts Bd infection intensity, and suggests that generating a sustained immune response against Bd may be counterproductive for surviving chytridiomycosis in this partially susceptible species.

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.

Xenopus Interferon Complex: Inscribing the Amphibiotic Adaption and Species-Specific Pathogenic Pressure in Vertebrate Evolution?

Several recent studies have revealed previously unknown complexity of the amphibian interferon (IFN) system. Being unique in vertebrate animals, amphibians not only conserve and multiply the fish-like intron-containing IFN genes, but also rapidly evolve amniote-like intronless IFN genes in each tested species. We postulate that the amphibian IFN system confers an essential model to study vertebrate immune evolution in molecular and functional diversity to cope with unprecedented pathophysiological requirement during terrestrial adaption. Studies so far have ascribed a potential role of these IFNs in immune regulation against intracellular pathogens, particularly viruses; however, many knowledge gaps remain elusive. Based on recent reports about IFN's multifunctional properties in regulation of animal physiological and defense responses, we interpret that amphibian IFNs may evolve novel function pertinent to their superior molecular diversity. Such new function revealed by the emerging studies about antifungal and developmental regulation of amphibian IFNs will certainly promote our understanding of immune evolution in vertebrates to address current pathogenic threats causing amphibian decline.

Moving Beyond the Host: Unraveling the Skin Microbiome of Endangered Costa Rican Amphibians

Some neotropical amphibians, including a few species in Costa Rica, were presumed to be "extinct" after dramatic population declines in the late 1980s but have been rediscovered in isolated populations. Such populations seem to have evolved a resistance/tolerance to Batrachochytrium dendrobatidis (Bd), a fungal pathogen that causes a deadly skin disease and is considered one of the main drivers of worldwide amphibian declines. The skin microbiome is an important component of the host's innate immune system and is associated with Bd-resistance. However, the way that the bacterial diversity of the skin microbiome confers protection against Bd in surviving species remains unclear. We studied variation in the skin microbiome and the prevalence of putatively anti-Bd bacterial taxa in four co-habiting species in the highlands of the Juan Castro Blanco National Park in Costa Rica using 16S rRNA amplicon sequencing. Lithobates vibicarius, Craugastor escoces, and Isthmohyla rivularis have recently been rediscovered, whereas Isthmohyla pseudopuma has suffered population fluctuations but has never disappeared. To investigate the life stage at which the protective skin microbiome is shaped and when shifts occur in the diversity of putatively anti-Bd bacteria, we studied the skin microbiome of tadpoles, juveniles and adults of L. vibicarius. We show that the skin bacterial composition of sympatric species and hosts with distinct Bd-infection statuses differs at the phyla, family, and genus level. We detected 94 amplicon sequence variants (ASVs) with putative anti-Bd activity pertaining to distinct bacterial taxa, e.g., Pseudomonas spp., Acinetobacter johnsonii, and Stenotrophomonas maltophilia. Bd-uninfected L. vibicarius harbored 79% more putatively anti-Bd ASVs than Bd-infected individuals. Although microbiome composition and structure differed across life stages, the diversity of putative anti-Bd bacteria was similar between pre- and post-metamorphic stages of L. vibicarius. Despite low sample size, our results support the idea that the skin microbiome is dynamic and protects against ongoing Bd presence in endangered species persisting after their presumed extinction. Our study serves as a baseline to understand the microbial patterns in species of high conservation value. Identification of microbial signatures linked to variation in disease susceptibility might, therefore, inform mitigation strategies for combating the global decline of amphibians.

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.

Metabolites Involved in Immune Evasion by Batrachochytrium dendrobatidis Include the Polyamine Spermidine

Amphibians have been declining around the world for more than four decades. One recognized driver of these declines is the chytrid fungus Batrachochytrium dendrobatidis, which causes the disease chytridiomycosis. Amphibians have complex and varied immune defenses against B. dendrobatidis, but the fungus also has a number of counterdefenses. Previously, we identified two small molecules produced by the fungus that inhibit frog lymphocyte proliferation, methylthioadenosine (MTA) and kynurenine (KYN). Here, we report on the isolation and identification of the polyamine spermidine (SPD) as another significant immunomodulatory molecule produced by B. dendrobatidis SPD and its precursor, putrescine (PUT), are the major polyamines detected, and SPD is required for growth. The major pathway of biosynthesis is from ornithine through putrescine to spermidine. An alternative pathway from arginine to agmatine to putrescine appears to be absent. SPD is inhibitory at concentrations of ≥10 μM and is found at concentrations between 1 and 10 μM in active fungal supernatants. Although PUT is detected in the fungal supernatants, it is not inhibitory to lymphocytes even at concentrations as high as 100 μM. Two other related polyamines, norspermidine (NSP) and spermine (SPM), also inhibit amphibian lymphocyte proliferation, but a third polyamine, cadaverine (CAD), does not. A suboptimal (noninhibitory) concentration of MTA (10 μM), a by-product of spermidine synthesis, enhances the inhibition of SPD at 1 and 10 μM. We interpret these results to suggest that B. dendrobatidis produces an "armamentarium" of small molecules that, alone or in concert, may help it to evade clearance by the amphibian immune system.

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.

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

Chytridiomycosis is among several recently emerged fungal diseases of wildlife that have caused decline or extinction of naïve populations. Despite recent advances in understanding pathogenesis, host response to infection remains poorly understood. Here we modelled a total of 162 metabolites across skin and liver tissues of 61 frogs from four populations (three long-exposed and one naïve to the fungus) of the Australian alpine tree frog (Litoria verreauxii alpina) throughout a longitudinal exposure experiment involving both infected and negative control individuals. We found that chytridiomycosis dramatically altered the organism-wide metabolism of clinically diseased frogs. Chytridiomycosis caused catastrophic failure of normal homeostatic mechanisms (interruption of biosynthetic and degradation metabolic pathways), and pronounced dysregulation of cellular energy metabolism. Key intermediates of the tricarboxylic acid cycle were markedly depleted, including in particular α-ketoglutarate and glutamate that together constitute a key nutrient pathway for immune processes. This study was the first to apply a non-targeted metabolomics approach to a fungal wildlife disease and specifically to dissect the host-pathogen interface of Bd-infected frogs. The patterns of metabolite accumulation we have identified reveal whole-body metabolic dysfunction induced by a fungal skin infection, and these findings have broad relevance for other fungal diseases.

Survival, gene and metabolite responses of Litoria verreauxii alpina frogs to fungal disease chytridiomycosis

The fungal skin disease chytridiomycosis has caused the devastating decline and extinction of hundreds of amphibian species globally, yet the potential for evolving resistance, and the underlying pathophysiological mechanisms remain poorly understood. We exposed 406 naïve, captive-raised alpine tree frogs (Litoria verreauxii alpina) from multiple populations (one evolutionarily naïve to chytridiomycosis) to the aetiological agent Batrachochytrium dendrobatidis in two concurrent and controlled infection experiments. We investigated (A) survival outcomes and clinical pathogen burdens between populations and clutches, and (B) individual host tissue responses to chytridiomycosis. Here we present multiple interrelated datasets associated with these exposure experiments, including animal signalment, survival and pathogen burden of 355 animals from Experiment A, and the following datasets related to 61 animals from Experiment B: animal signalment and pathogen burden; raw RNA-Seq reads from skin, liver and spleen tissues; de novo assembled transcriptomes for each tissue type; raw gene expression data; annotation data for each gene; and raw metabolite expression data from skin and liver tissues. These data provide an extensive baseline for future analyses.