Major histocompatibility complex variation and the evolution of resistance to amphibian chytridiomycosis

Evolutionary analysis of major histocompatibility complex variants in chytrid-resistant and susceptible amphibians

An amphibian emerging infectious disease (EID), chytridiomycosis, caused by Batrachochytrium dendrobatidis (Bd), originated in Asia but primarily led to declines and extinctions in amphibian populations outside of Asia. Host major histocompatibility complex (MHC) molecules exhibit high polymorphism, and the evolution of MHC can be influenced by recombination and pathogens. Previous studies have indicated that host MHC class II is associated with Bd resistance. In this study, I conducted recombination and selection tests on functional MHC IIß1 alleles from an Asian Bd-resistant anuran species (Bufo gargarizans) and an Australasian Bd-susceptible species (Litoria caerulea). Recombination at the same site was identified in both species, supporting the hypothesis that recombination contributes to MHC IIß1 diversity in amphibians. Positive selection was observed in MHC IIß1 alleles in both species. In L. caerulea, at least four amino acid sites were identified under significant positive selection in the MHC IIß1, whereas these sites were either negatively selected or conserved in B. gargarizans. This suggests these sites might be selected for Bd resistance. Hydrophobicity was detected in certain amino acid sites relating to Bd resistance, suggesting this physicochemical property may be a factor selected to counteract Bd infection. These findings of this study provide an evolutionary basis for understanding how amphibian MHC IIß1 may undergo selection in response to chytrid infection.

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.

The wood frog (Rana sylvatica): An emerging comparative model for anuran immunity and host-ranavirus interactions

The wood frog (Rana sylvatica) is widely distributed across North America and is the only amphibian found north of the Arctic Circle due to its remarkable ability to tolerate whole-body freezing. Recent mass mortalities attributable to Ranavirus spp. (family Iridoviridae) in wild juvenile wood frogs, coupled with the apparent high susceptibility of wood frogs to experimental infection with frog virus 3 (FV3), the type species of the Ranavirus genus, or FV3-like isolates underscore the serious threat ranaviruses poses to wood frog populations. Despite the ecological relevance and unique life history of wood frogs, our understanding of the wood frog immune system and antiviral response to ranaviral infections is in its infancy. Here we aim to (1) synthesize the limited knowledge of wood frog immune defences, (2) review recent progress in establishing the wood frog as a study system for ranavirus infection, and (3) highlight the future use of wood frogs as a model anuran to provide insight into the evolution of anuran immune systems and antiviral responses.

Identification of major histocompatibility complex genotypes associated with resistance to an amphibian emerging infectious disease

Amphibian chytridiomycosis, caused by Batrachochytrium dendrobatidis (Bd), emerged from Asia and spread globally. By comparing functional MHC IIß1 alleles from an Asian Bd-resistant anuran species (Bufo gargarizans) with those of an Australasian Bd-susceptible species (Litoria caerulea), we identified MHC genotypes associated with Bd resistance. These alleles encode a glycine deletion (G90β1) and adjacent motifs in the deepest pathogen-derived peptide-binding groove. Every Bd-resistant individual, but no susceptible individuals, possessed at least one allele encoding the variant. We detected trans-species polymorphism at the end of the MHC IIβ1 sequences. The G90β1 deletion was encoded by different alleles in the two species, suggesting it may have evolved independently in each species rather than having been derived from a common ancestor. These results are consistent with a scenario by which MHC adaptations that confer resistance to the pathogen have evolved by convergent evolution. Immunogenetic studies such as this are critical to ongoing conservation efforts.

The amphibian immune system

Amphibians are at the forefront of bridging the evolutionary gap between mammals and more ancient, jawed vertebrates. Currently, several diseases have targeted amphibians and understanding their immune system has importance beyond their use as a research model. The immune system of the African clawed frog, Xenopus laevis, and that of mammals is well conserved. We know that several features of the adaptive and innate immune system are very similar for both, including the existence of B cells, T cells and innate-like T cells. In particular, the study of the immune system at early stages of development is benefitted by studying X. laevis tadpoles. The tadpoles mainly rely on innate immune mechanisms including pre-set or innate-like T cells until after metamorphosis. In this review we lay out what is known about the innate and adaptive immune system of X. laevis including the lymphoid organs as well as how other amphibian immune systems are similar or different. Furthermore, we will describe how the amphibian immune system responds to some viral, bacterial and fungal insults. This article is part of the theme issue 'Amphibian immunity: stress, disease and ecoimmunology'.

Divergent allele advantage in the MHC and amphibian emerging infectious disease

Genetic variation in the major histocompatibility complex (MHC) may be associated with resistance to the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). The pathogen originated in Asia, then spread worldwide, causing amphibian population declines and species extinctions. We compared the expressed MHC IIβ1 alleles of a Bd-resistant toad species, Bufo gargarizans, from South Korea with those of a Bd-susceptible Australasian frog species, Litoria caerulea. We found at least six expressed MHC IIβ1 loci in each of the two species. Amino acid diversity encoded by these MHC alleles was similar between species, but the genetic divergence of those alleles known for broader pathogen-derived peptide binding was greater in the Bd-resistant species. In addition, we found a potentially rare allele in one resistant individual from the Bd-susceptible species. Deep next-generation sequencing recovered approximately triple the genetic resolution accessible from traditional cloning-based genotyping. Targeting more than one MHC IIβ1 expressed locus enables us to better understand how host MHC may adapt to emerging infectious diseases.

Disease outbreaks select for mate choice and coat color in wolves

We know much about pathogen evolution and the emergence of new disease strains, but less about host resistance and how it is signaled to other individuals and subsequently maintained. The cline in frequency of black-coated wolves (Canis lupus) across North America is hypothesized to result from a relationship with canine distemper virus (CDV) outbreaks. We tested this hypothesis using cross-sectional data from wolf populations across North America that vary in the prevalence of CDV and the allele that makes coats black, longitudinal data from Yellowstone National Park, and modeling. We found that the frequency of CDV outbreaks generates fluctuating selection that results in heterozygote advantage that in turn affects the frequency of the black allele, optimal mating behavior, and black wolf cline across the continent.

Novel chytrid pathogen variants and the global amphibian pet trade

Global wildlife trade spreads emerging infectious diseases that threaten biodiversity. The amphibian chytrid pathogen Batrachochytrium dendrobatidis (Bd) has caused population declines and species extinctions worldwide except in Asia. Fire-bellied toads (Bombina orientalis), exported in large numbers from Asia, are tolerant of Bd and carry hypervirulent ancestral chytrid BdAsia-1 variants. We assayed the virulence of a new isolate of BdAsia-1 on the model Australasian frog host Litoria caerulea. Infected individuals (n = 15) all showed rapid disease progression culminating in death, whereas sham-inoculated individuals (n = 10) presented no clinical signs of disease and all survived (log rank test, χ² = 15.6, df = 1, p < 0.0001). The virulence of the new isolate of BdAsia-1 is comparable to the one we assayed previously (χ² = 0.0, df = 1, p = 0.91). Internationally traded wildlife, even when they appear healthy, can carry hypervirulent variants of pathogens. Once new pathogen variants escape into the environment, native species that have had no opportunity to evolve resistance to them may perish. Our study suggests that hypervirulent pathogens are being spread by the international pet trade. Notifiable wildlife diseases attributable to locally endemic pathogens often fail to generate conservation concern so are rarely subject to border surveillance or import controls. Because of the danger novel variants pose, national border control agencies need to implement disease screening and quarantine protocols to ensure the safety of their endemic fauna.

Low neutral and immunogenetic diversity in northern fringe populations of the green toad Bufotes viridis: implications for conservation

Genetic variation is often lower at high latitudes, which may compromise the adaptability and hence survival of organisms. Here we show that genetic variability is negatively correlated with northern latitude in European green toads (Bufotes viridis). The result holds true for both putatively neutral microsatellite variation and supposedly adaptive MHC Class IIB variation. In particular, our findings have bearing on the conservation status of this species in Sweden, on the northern limit of its distribution where local populations are small and fragmented. These genetically impoverished populations are closely related to other populations found around the Baltic Sea basin. The low neutral and adaptive variation in these fringe populations compared to population at central ranges confirms a pattern shared across all other amphibians so far studied. In Sweden, the situation of green toads is of concern as the remaining populations may not have the evolutionary potential to cope with present and future environmental challenges.

Spatio-temporal analyses reveal infectious disease-driven selection in a free-ranging ungulate

Infectious diseases play an important role in wildlife population dynamics by altering individual fitness, but detecting disease-driven natural selection in free-ranging populations is difficult due to complex disease-host relationships. Chronic wasting disease (CWD) is a fatal infectious prion disease in cervids for which mutations in a single gene have been mechanistically linked to disease outcomes, providing a rare opportunity to study disease-driven selection in wildlife. In Wyoming, USA, CWD has gradually spread across mule deer (Odocoileus hemionus) populations, producing natural variation in disease history to evaluate selection pressure. We used spatial variation and a novel temporal comparison to investigate the relationship between CWD and a mutation at codon 225 of the mule deer prion protein gene that slows disease progression. We found that individuals with the 'slow' 225F allele were less likely to test positive for CWD, and the 225F allele was more common in herds exposed to CWD longer. We also found that in the past 2 decades, the 225F allele frequency increased more in herds with higher CWD prevalence. This study expanded on previous research by analysing spatio-temporal patterns of individual and herd-based disease data to present multiple lines of evidence for disease-driven selection in free-ranging wildlife.

Towards a more healthy conservation paradigm: integrating disease and molecular ecology to aid biological conservation†

Parasites, and the diseases they cause, are important from an ecological and evolutionary perspective because they can negatively affect host fitness and can regulate host populations. Consequently, conservation biology has long recognized the vital role that parasites can play in the process of species endangerment and recovery. However, we are only beginning to understand how deeply parasites are embedded in ecological systems, and there is a growing recognition of the important ways in which parasites affect ecosystem structure and function. Thus, there is an urgent need to revisit how parasites are viewed from a conservation perspective and broaden the role that disease ecology plays in conservation-related research and outcomes. This review broadly focusses on the role that disease ecology can play in biological conservation. Our review specifically emphasizes on how the integration of tools and analytical approaches associated with both disease and molecular ecology can be leveraged to aid conservation biology. Our review first concentrates on disease-mediated extinctions and wildlife epidemics. We then focus on elucidating how host–parasite interactions has improved our understanding of the eco-evolutionary dynamics affecting hosts at the individual, population, community and ecosystem scales. We believe that the role of parasites as drivers and indicators of ecosystem health is especially an exciting area of research that has the potential to fundamentally alter our view of parasites and their role in biological conservation. The review concludes with a broad overview of the current and potential applications of modern genomic tools in disease ecology to aid biological conservation.

Infection with Batrachochytrium dendrobatidis reduces salamander capacity to mount a cell-mediated immune response

The vertebrate immune system is a costly defense system that is responsible for preventing and eliminating parasites and pathogens. Theory predicts that hosts experience tradeoffs associated with immune deployment and other physiological functions. Although empirical evidence for immune-physiology tradeoffs are well documented in the literature, fewer studies have examined tradeoffs within the immune system in wild vertebrates. We explored the topic of concomitant immune challenges in amphibians by exposing salamanders (Plethodon cinereus) to a fungal pathogen Batrachochytrium dendrobatidis (hereafter "Bd") and then to phytohemagglutinin (hereafter "PHA"). We measured Bd infection using quantitative PCR and used measurements of the tail thickness at the PHA injection site as an estimate of skin swelling. We tested whether Bd reduced the salamander's capacity to mount an immune response towards PHA or whether Bd would stimulate immune activity and thereby increase the response towards PHA. Salamanders that were infected with Bd had a reduced skin-swelling when injected with PHA compared to noninfected salamanders, a result that is consistent with the hypothesis that Bd-infected salamanders have lower immunocompetence than noninfected salamanders. We also found that PHA-induced swelling response was negatively associated with Bd infection abundance (i.e., the infection burden of all exposed salamanders, including those that were exposed but not infected), indicating that salamanders with a higher infection abundance had the lowest swelling response to PHA. Our results suggest that individuals of P. cinereus might experience an energetic tradeoff between successfully fighting off Bd and mounting an immune response towards PHA.

Spatiotemporal adaptive evolution of an MHC immune gene in a frog-fungus disease system

Genetic diversity of major histocompatibility complex (MHC) genes is linked to reduced pathogen susceptibility in amphibians, but few studies also examine broad spatial and temporal patterns of MHC and neutral genetic diversity. Here, we characterized range-wide MHC diversity in the Northern leopard frog, Rana pipiens, a species found throughout North America that is experiencing disease-related declines. We used previously sequenced neutral markers (mitochondrial DNA and microsatellites), sequenced an expressed MHC class IIß gene fragment, and measured infection prevalence and intensity of the global fungal pathogen Batrachochytrium dendrobatidis (Bd) across 14 populations. Four populations were sampled across two decades, enabling temporal comparisons of selection and demography. We recovered 37 unique MHC alleles, including 17 that were shared across populations. Phylogenetic and population genetic patterns between MHC and neutral markers were incongruent, and five MHC codon positions associated with peptide binding were under positive selection. MHC heterozygosity, but not neutral marker heterozygosity, was a significant factor explaining spatial patterns of Bd prevalence, whereas only environmental variables predicted Bd intensity. MHC allelic richness (AR) decreased significantly over time but microsatellite-based AR did not, highlighting a loss of functional immunogenetic diversity that may be associated with Bd selective pressures. MHC supertype 4 was significantly associated with an elevated risk of Bd infection, whereas one supertype 2 allele was associated with a nearly significant reduced risk of Bd. Taken together, these results provide evidence that positive selection contributes to MHC class IIß evolution in R. pipiens and suggest that functional MHC differences across populations may contribute to disease adaptation.

Locality, time and heterozygosity affect chytrid infection in yellow-bellied toads

The chytrid fungus Batrachochytrium dendrobatidis (Bd) infects numerous amphibian species worldwide and is suggested to drive population declines and extinction events. We report a study of Bd infection at the northernmost distribution of the European yellow-bellied toad Bombina variegata. A total of 577 individuals from ponds in 16 study sites were sampled for DNA and Bd throughout the breeding season. Microsatellite genotyping revealed 3 genetic clusters for the host B. variegata with an overall low genetic diversity. One of the clusters displayed a low microsatellite heterozygosity, a high inbreeding coefficient as well as high Bd infection prevalence and intensities. Multi-model estimates identified site, time of sampling, and heterozygosity to be important predictors of an individual's Bd infection status, and identified a strong effect of site on individual Bd infection intensity. The study site effects are suggestive of localized infection peaks, and the increase of individual Bd infection probabilities towards the end of the sampling period suggests cumulative infection during the breeding season. This study highlights the need for regular monitoring of Bd infection variables at multiple localities and times to gain insights into Bd dynamics. Due to the detected relationship between individual Bd infection status and heterozygosity, conservation measures should focus on the maintenance of high genetic diversity and connectivity within and among amphibian populations.

Connectivity over a disease risk gradient enables recovery of rainforest frogs

Chytridiomycosis has been a key driver of global frog declines and extinctions, particularly for high-altitude populations across Australia and the Americas. While recent evidence shows some species are recovering, the extent of such recoveries and the mechanisms underpinning them remain poorly resolved. We surveyed the historical latitudinal and elevational range of four Australian rainforest frogs that disappeared from upland sites between 1989 and 1994 to establish their contemporary distribution and elevational limits, and investigate factors affecting population recovery. Five rainforest streams were surveyed from mountain-base to summit (30 sites in total), with swabs collected from the target species (Litoria dayi, L. nannotis, L. rheocola, and L. serrata) to determine their infection status, and data loggers deployed to measure microclimatic variation across the elevational gradient. Infection probability increased with elevation and canopy cover as it was tightly and inversely correlated with stream-side air temperature. Occupancy patterns suggest varying responses to this disease threat gradient. Two species, L. rheocola and L. serrata, were found over their full historical elevational range (≥1,000 m above sea level [asl]), while L. dayi was not detected above 400 m (formerly known up to 900 m asl) and L. nannotis was not detected above 800 m (formerly known up to 1,200 m asl). Site occupancy probability was negatively related to predicted infection prevalence for L. dayi, L. nannotis, and L. rheocola, but not L. serrata, which appears to now tolerate high fungal burdens. This study highlights the importance of environmental refuges and connectivity across disease risk gradients for the persistence and natural recovery of frogs susceptible to chytridiomycosis. Likewise, in documenting both interspecific variation in recovery rates and intraspecific differences between sites, this study suggests interactions between disease threats and host selection exist that could be manipulated. For example, translocations may be warranted where connectivity is poor or the increase in disease risk is too steep to allow recolonization, combined with assisted selection or use of founders from populations that have already undergone natural selection.

Selected Emerging Infectious Diseases of Amphibians

This article updates the understanding of two extirpation-driving infectious diseases, Batrachochytrium dendrobatidis and Batrachochytrium salamandrivorans, and Ranavirus. Experimental studies and dynamic, multifactorial population modeling have outlined the epidemiology and future population impacts of B dendrobatidis, B salamandrivorans, and Ranavirus. New genomic findings on divergent fungal and viral pathogens can help optimize control and disease management strategies. Although there have been major advances in knowledge of amphibian pathogens, controlled studies are needed to guide population recovery to elucidate and evaluate transmission routes for several pathogens, examine environmental control, and validate new diagnostic tools to confirm the presence of disease.

Depauperate major histocompatibility complex variation in the endangered reticulated flatwoods salamander (Ambystoma bishopi)

Reticulated flatwoods salamander (Ambystoma bishopi) populations began decreasing dramatically in the 1900s. Contemporary populations are small, isolated, and may be susceptible to inbreeding and reduced adaptive potential because of low genetic variation. Genetic variation at immune genes is especially important as it influences disease susceptibility and adaptation to emerging infectious pathogens, a central conservation concern for declining amphibians. We collected samples from across the extant range of this salamander to examine genetic variation at major histocompatibility complex (MHC) class Iα and IIβ exons as well as the mitochondrial control region. We screened tail or toe tissue for ranavirus, a pathogen associated with amphibian declines worldwide. Overall, we found low MHC variation when compared to other amphibian species and did not detect ranavirus at any site. MHC class Iα sequencing revealed only three alleles with a nucleotide diversity of 0.001, while MHC class IIβ had five alleles with a with nucleotide diversity of 0.004. However, unique variation still exists across this species' range with private alleles at three sites. Unlike MHC diversity, mitochondrial variation was comparable to levels estimated for other amphibians with nine haplotypes observed, including one haplotype shared across all sites. We hypothesize that a combination of a historic disease outbreak and a population bottleneck may have contributed to low MHC diversity while maintaining higher levels of mitochondrial DNA variation. Ultimately, MHC data indicated that the reticulated flatwoods salamander may be at an elevated risk from infectious diseases due to low levels of immunogenetic variation necessary to combat novel pathogens.

Serratia marcescens shapes cutaneous bacterial communities and influences survival of an amphibian host

Ongoing investigations into the interactions between microbial communities and their associated hosts are changing how emerging diseases are perceived and ameliorated. Of the numerous host-microbiome-disease systems of study, the emergence of chytridiomycosis (caused by Batrachochytrium dendrobatidis, hereafter Bd) has been implicated in ongoing declines and extinction events of amphibians worldwide. Interestingly, there has been differential survival among amphibians in resisting Bd infection and subsequent disease. One factor thought to contribute to this resistance is the host-associated cutaneous microbiota. This has raised the possibility of using genetically modified probiotics to restructure the host-associated microbiota for desired anti-fungal outcomes. Here, we use a previously described strain of Serratia marcescens (Sm) for the manipulation of amphibian cutaneous microbiota. Sm was genetically altered to have a dysfunctional pathway for the production of the extracellular metabolite prodigiosin. This genetically altered strain (Δpig) and the functional prodigiosin producing strain (wild-type, WT) were compared for their microbial community and anti-Bd effects both in vitro and in vivo. In vitro, Bd growth was significantly repressed in the presence of prodigiosin. In vivo, the inoculation of both Sm strains was shown to significantly influence amphibian microbiota diversity with the Δpig-Sm treatment showing increasing alpha diversity, and the WT-Sm having no temporal effect on diversity. Differences were also seen in host mortality with Δpig-Sm treatments exhibiting significantly decreased survival probability when compared with WT-Sm in the presence of Bd. These results are an important proof-of-concept for linking the use of genetically modified probiotic bacteria to host microbial community structure and disease outcomes, which in the future may provide a way to ameliorate disease and address critical frontiers in disease and microbial ecology.

Ancestral chytrid pathogen remains hypervirulent following its long coevolution with amphibian hosts

Many amphibian species around the world, except in Asia, suffer morbidity and mortality when infected by the emerging infectious pathogen Batrachochytrium dendrobatidis (Bd). A lineage of the amphibian chytrid fungus isolated from South Korean amphibians (BdAsia-1) is evolutionarily basal to recombinant global pandemic lineages (BdGPL) associated with worldwide amphibian population declines. In Asia, the Bd pathogen and its amphibian hosts have coevolved over 100 years or more. Thus, resilience of Asian amphibian populations to infection might result from attenuated virulence of endemic Bd lineages, evolved immunity to the pathogen or both. We compared susceptibilities of an Australasian amphibian, Litoria caerulea, known to lack resistance to BdGPL, with those of three Korean species, Bufo gargarizans, Bombina orientalis and Hyla japonica, after inoculation with BdAsia-1, BdGPL or a blank solution. Subjects became infected in all experimental treatments but Korean species rapidly cleared themselves of infection, regardless of Bd lineage. They survived with no apparent secondary effects. By contrast, L. caerulea, after infection by either BdAsia-1 or BdGPL, suffered deteriorating body condition and carried progressively higher Bd loads over time. Subsequently, most subjects died. Comparing their effects on L. caerulea, BdAsia-1 induced more rapid disease progression than BdGPL. The results suggest that genomic recombination with other lineages was not necessary for the ancestral Bd lineage to evolve hypervirulence over its long period of coevolution with amphibian hosts. The pathogen's virulence may have driven strong selection for immune responses in endemic Asian amphibian host species.

Functional variation at an expressed MHC class IIβ locus associates with Ranavirus infection intensity in larval anuran populations

Infectious diseases are causing catastrophic losses to global biodiversity. Iridoviruses in the genus Ranavirus are among the leading causes of amphibian disease-related mortality. Polymorphisms in major histocompatibility complex (MHC) genes are significantly associated with variation in amphibian pathogen susceptibility. MHC genes encode two classes of polymorphic cell-surface molecules that can recognize and bind to diverse pathogen peptides. While MHC class I genes are the classic mediators of viral-acquired immunity, larval amphibians do not express them. Consequently, MHC class II gene diversity may be an important predictor of Ranavirus susceptibility in larval amphibians, the life stage most susceptible to Ranavirus. We surveyed natural populations of larval wood frogs (Rana sylvatica), which are highly susceptible to Ranavirus, across 17 ponds and 2 years in Maryland, USA. We sequenced the peptide-binding region of an expressed MHC class IIβ locus and assessed allelic and genetic diversity. We converted alleles to functional supertypes and determined if supertypes or alleles influenced host responses to Ranavirus. Among 381 sampled individuals, 26% were infected with Ranavirus. We recovered 20 unique MHC class IIβ alleles that fell into two deeply diverged clades and seven supertypes. MHC genotypes were associated with Ranavirus infection intensity, but not prevalence. Specifically, MHC heterozygotes and supertype ST1/ST7 had significantly lower Ranavirus infection intensity compared to homozygotes and other supertypes. We conclude that MHC class IIβ functional genetic variation is an important component of Ranavirus susceptibility. Identifying immunogenetic signatures linked to variation in disease susceptibility can inform mitigation strategies for combatting global amphibian declines.

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.

Effects of host species and environmental factors on the prevalence of Batrachochytrium dendrobatidis in northern Europe

The fungal pathogen Batrachochytrium dendrobatidis (Bd) poses a major threat to amphibian populations. To assist efforts to address such threats, we examined differences in Bd host infection prevalence among amphibian species and its relations to both local environmental factors in breeding habitats and landscape variables measured at three scales (500, 2000 and 5000 m radii) around breeding sites in southernmost Sweden. We sampled 947 anurans of six species in 31 ponds and assessed their infection status. We then examined correlations of infection prevalence with canopy cover, pond perimeter and pH (treated as local-scale pond characteristics), and the number of ponds, area of arable land, area of mature forest, number of resident people and presence of sea within the three radii (treated as landscape variables). The Bd infection prevalence was very low, 0.5–1.0%, in two of the six anuran species (Bufo bufo and Rana temporaria), and substantially higher (13–64%) in the other four (Bombina bombina, Bufotes variabilis, Epidalea calamita, Rana arvalis). In the latter four species Bd infection prevalence was positively associated with ponds’ pH (site range: 5.3–8.1), and negatively associated with areas of mature forest and/or wetlands in the surroundings. Our results show that the infection dynamics of Bd are complex and associated with host species, local pond characteristics and several landscape variables at larger spatial scales. Knowledge of environmental factors associated with Bd infections and differences in species’ susceptibility may help to counter further spread of the disease and guide conservation action plans, especially for the most threatened species.

Toward an integrative molecular approach to wildlife disease

Pathogens pose serious threats to human health, agricultural investment, and biodiversity conservation through the emergence of zoonoses, spillover to domestic livestock, and epizootic outbreaks. As such, wildlife managers are often tasked with mitigating the negative effects of disease. Yet, parasites form a major component of biodiversity that often persist. This is due to logistical challenges of implementing management strategies and to insufficient understanding of host-parasite dynamics. We advocate for an inclusive understanding of molecular diversity in driving parasite infection and variable host disease states in wildlife systems. More specifically, we examine the roles of genetic, epigenetic, and commensal microbial variation in disease pathogenesis. These include mechanisms underlying parasite virulence and host resistance and tolerance, and the development, regulation, and parasite subversion of immune pathways, among other processes. Case studies of devil facial tumor disease in Tasmanian devils (Sarcophilus harrisii) and chytridiomycosis in globally distributed amphibians exemplify the broad range of questions that can be addressed by examining different facets of molecular diversity. For particularly complex systems, integrative molecular analyses present a promising frontier that can provide critical insights necessary to elucidate disease dynamics operating across scales. These insights enable more accurate risk assessment, reconstruction of transmission pathways, discernment of optimal intervention strategies, and development of more effective and ecologically sound treatments that minimize damage to the host population and environment. Such measures are crucial when mitigating threats posed by wildlife disease to humans, domestic animals, and species of conservation concern.

Gene expression differs in susceptible and resistant amphibians exposed to Batrachochytrium dendrobatidis

Chytridiomycosis, the disease caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), has devastated global amphibian biodiversity. Nevertheless, some hosts avoid disease after Bd exposure even as others experience near-complete extirpation. It remains unclear whether the amphibian adaptive immune system plays a role in Bd defence. Here, we describe gene expression in two host species-one susceptible to chytridiomycosis and one resistant-following exposure to two Bd isolates that differ in virulence. Susceptible wood frogs (Rana sylvatica) had high infection loads and mortality when exposed to the more virulent Bd isolate but lower infection loads and no fatal disease when exposed to the less virulent isolate. Resistant American bullfrogs (R. catesbeiana) had high survival across treatments and rapidly cleared Bd infection or avoided infection entirely. We found widespread upregulation of adaptive immune genes and downregulation of important metabolic and cellular maintenance components in wood frogs after Bd exposure, whereas American bullfrogs showed little gene expression change and no evidence of an adaptive immune response. Wood frog responses suggest that adaptive immune defences may be ineffective against virulent Bd isolates that can cause rapid physiological dysfunction. By contrast, American bullfrogs exhibited robust resistance to Bd that is likely attributable, at least in part, to their continued upkeep of metabolic and skin integrity pathways as well as greater antimicrobial peptide expression compared to wood frogs, regardless of exposure. Greater understanding of these defences will ultimately help conservationists manage chytridiomycosis.