Adaptive tolerance to a pathogenic fungus drives major histocompatibility complex evolution in natural amphibian populations

The Amphibian Major Histocompatibility Complex-A Review and Future Outlook

The major histocompatibility complex (MHC) is a cluster of functionally related genes encoding proteins which, among other functions, mediate immune system activation. While the MHC of many vertebrates has been extensively studied, less is known about the amphibian MHC. This represents an important knowledge gap because amphibians mark the evolutionary transition from an aquatic to a terrestrial lifestyle and often maintain a biphasic lifestyle. Hence, they tend to be exposed to both aquatic and terrestrial pathogen communities, providing opportunities to gain fundamental insights into how the immune system responds to different environmental challenges. Moreover, amphibians are globally threatened by invasive pathogens and the MHC may play a role in combating population decline. In this review, we summarize the current state of knowledge regarding the amphibian MHC and identify the major differences with other vertebrates. We also review how the number of MHC gene copies varies across amphibian groups and how MHC-based variation relates to amphibian ontogeny, behaviour, disease, and phylogeography. We conclude by identifying knowledge gaps and proposing priorities for future research.

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.

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.

Varying conjunctival immune response adaptations of house finch populations to a rapidly evolving bacterial pathogen

Pathogen adaptations during host-pathogen co-evolution can cause the host balance between immunity and immunopathology to rapidly shift. However, little is known in natural disease systems about the immunological pathways optimised through the trade-off between immunity and self-damage. The evolutionary interaction between the conjunctival bacterial infection Mycoplasma gallisepticum (MG) and its avian host, the house finch (Haemorhous mexicanus), can provide insights into such adaptations in immune regulation. Here we use experimental infections to reveal immune variation in conjunctival tissue for house finches captured from four distinct populations differing in the length of their co-evolutionary histories with MG and their disease tolerance (defined as disease severity per pathogen load) in controlled infection studies. To differentiate contributions of host versus pathogen evolution, we compared house finch responses to one of two MG isolates: the original VA1994 isolate and a more evolutionarily derived one, VA2013. To identify differential gene expression involved in initiation of the immune response to MG, we performed 3'-end transcriptomic sequencing (QuantSeq) of samples from the infection site, conjunctiva, collected 3-days post-infection. In response to MG, we observed an increase in general pro-inflammatory signalling, as well as T-cell activation and IL17 pathway differentiation, associated with a decrease in the IL12/IL23 pathway signalling. The immune response was stronger in response to the evolutionarily derived MG isolate compared to the original one, consistent with known increases in MG virulence over time. The host populations differed namely in pre-activation immune gene expression, suggesting population-specific adaptations. Compared to other populations, finches from Virginia, which have the longest co-evolutionary history with MG, showed significantly higher expression of anti-inflammatory genes and Th1 mediators. This may explain the evolution of disease tolerance to MG infection in VA birds. We also show a potential modulating role of BCL10, a positive B- and T-cell regulator activating the NFKB signalling. Our results illuminate potential mechanisms of house finch adaptation to MG-induced immunopathology, contributing to understanding of the host evolutionary responses to pathogen-driven shifts in immunity-immunopathology trade-offs.

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.

MassArray Genotyping as a Selection Tool for Extending the Shelf-Life of Fresh Gilthead Sea Bream and European Seabass

In modern aquaculture, genomics-driven breeding programs have emerged as powerful tools for optimizing fish quality. This study focused on two emblematic Mediterranean fish species, the European seabass (Dicentrarchus labrax) and the gilthead sea bream (Sparus aurata), with a primary aim of exploring the genetic basis of white muscle/fillet degradation in fresh fish following harvest. We identified 57 and 44 missense SNPs in gilthead sea bream and European seabass, respectively, located within genes encoding for endogenous proteases responsible for fillet quality. These SNPs were cherry-picked based on their strategic location within the catalytic/regulatory domains of endogenous proteases that are expressed in the white muscle. Using MassArray technology, we successfully associated differentiated enzymatic activity of those endogenous proteases post-harvest as a phenotypic trait with genetic polymorphism of six SNPs in gilthead sea bream and nine in European seabass. These findings can be valuable attributes in selective breeding programs toward the extension of freshness and shelf life of these species. The integration of MassArray technology into breeding programs offers a cost-effective strategy for harnessing the potential of these genetic variants to enhance the overall quality of the final product. Recognizing that fresh fish perishability is a challenge, extending shelf-life is pivotal in reducing losses and production costs.

Body size mediates latitudinal population differences in the response to chytrid fungus infection in two amphibians

Factors behind intraspecific variation in sensitivity to pathogens remain poorly understood. We investigated how geographical origin in two North European amphibians affects tolerance to infection by the chytrid fungus Batrachochytrium dendrobatidis (Bd), a generalist pathogen which has caused amphibian population declines worldwide. We exposed newly metamorphosed individuals of moor frog Rana arvalis and common toad Bufo bufo from two latitudinal regions to two different BdGPL strains. We measured survival and growth as infections may cause sub-lethal effects in fitness components even in the absence of mortality. Infection loads were higher in B. bufo than in R. arvalis, and smaller individuals had generally higher infection loads. B. bufo had high mortality in response to Bd infection, whereas there was little mortality in R. arvalis. Bd-mediated mortality was size-dependent and high-latitude individuals were smaller leading to high mortality in the northern B. bufo. Bd exposure led to sub-lethal effects in terms of reduced growth suggesting that individuals surviving the infection may have reduced fitness mediated by smaller body size. In both host species, the Swedish Bd strain caused stronger sublethal effects than the British strain. We suggest that high-latitude populations can be more vulnerable to chytrids than those from lower latitudes and discuss the possible mechanisms how body size and host geographical origin contribute to the present results.

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.

Widespread amphibian Perkinsea infections associated with Ranidae hosts, cooler months and Ranavirus co-infection

Amphibians suffer from large-scale population declines globally, and emerging infectious diseases contribute heavily to these declines. Amphibian Perkinsea (Pr) is a worldwide anuran pathogen associated with mass mortality events, yet little is known about its epidemiological patterns, especially in comparison to the body of literature on amphibian chytridiomycosis and ranavirosis. Here, we establish Pr infection patterns in natural anuran populations and identify important covariates including climate, host attributes and co-infection with Ranavirus (Rv). We used quantitative (q)PCR to determine the presence and intensity of Pr and Rv across 1234 individuals sampled throughout central Florida in 2017-2019. We then implemented random forest ensemble learning models to predict infection with both pathogens based on physiological and environmental characteristics. Perkinsea infected 32% of all sampled anurans, and Pr prevalence was significantly elevated in Ranidae frogs, cooler months, metamorphosed individuals and frogs co-infected with Rv, while Pr intensity was significantly higher in ranid frogs and individuals collected dead. Ranavirus prevalence was 17% overall and was significantly higher in Ranidae frogs, metamorphosed individuals, locations with higher average temperatures, and individuals co-infected with Pr. Perkinsea prevalence was significantly higher than Rv prevalence across months, regions, life stages and species. Among locations, Pr prevalence was negatively associated with crayfish prevalence and positively associated with relative abundance of microhylids, but Rv prevalence did not associate with any tested co-variates. Co-infections were significantly more common than single infections for both pathogens, and we propose that Pr infections may propel Rv infections because seasonal Rv infection peaks followed Pr infection peaks and random forest models found Pr intensity was a leading factor explaining Rv infections. Our study elucidates epidemiological patterns of Pr in Florida and suggests that Pr may be under-recognized as a cause of anuran declines, especially in the context of pathogen co-infection.

Metabolites of Xenorhabdus bacteria are potent candidates for mitigating amphibian chytridiomycosis

Chytridiomycosis, caused by the chytrid fungus Batrachochytrium dendrobatidis (Bd), has caused extreme losses in amphibian biodiversity. Finding bacteria that produce metabolites with antifungal properties may turn out to be invaluable in the fight against this devastating disease. The entomopathogenic bacteria, Xenorhabdus szentirmaii and X. budapestensis produce secondary metabolites that are effective against a wide range of fungal plant pathogens. To assess whether they may also be effective against Bd, we extracted cell-free culture media (CFCM) from liquid cultures of X. szentirmaii and X. budapestensis and tested their ability to inhibit Bd growth in vitro. As a second step, using juvenile common toads (Bufo bufo) experimentally infected with Bd we also tested the in vivo antifungal efficacy of X. szentirmaii CFCM diluted to 2 and 10% (v/v), while also assessing possible malign side effects on amphibians. Results of the in vitro experiment documented highly effective growth inhibition by CFCMs of both Xenorhabdus species. The in vivo experiment showed that treatment with CFCM of X. szentirmaii applied at a dilution of 10% resulted in infection intensities reduced by ca. 73% compared to controls and to juvenile toads treated with CFCM applied at a dilution of 2%. At the same time, we detected no negative side effects of treatment with CFCM on toad survival and development. Our results clearly support the idea that metabolites of X. szentirmaii, and perhaps of several other Xenorhabdus species as well, may prove highly useful for the treatment of Bd infected amphibians.

Pathogenic selection promotes adaptive immune variations against serious bottlenecks in early invasions of bullfrogs

Adaptive genetic variations are key for understanding evolutionary processes influencing invasions. However, we have limited knowledge on how adaptive genetic diversity in invasive species responds to new pathogenic environments. Here, we compared variations in immune major histocompatibility complex (MHC) class-II β gene and neutral loci in relation to pathogenic chytrid fungus (Batrachochytrium dendrobatidis, Bd) infection across invasive and native populations of American bullfrog between China and United States (US). Chinese invasive populations show a 60% reduction in neutral cytb variations relative to US native populations, and there were similar MHC variation and functional diversity between them. One MHC allele private to China was under recent positive selection and associated with decreased Bd infection, partly explaining the lower Bd prevalence for Chinese populations than for native US populations. These results suggest that pathogen-mediated selection favors adaptive MHC variations and functional diversity maintenance against serious bottlenecks during the early invasions (within 15 generations) of bullfrogs.

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.

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.

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.

Diversifying selection and climatic effects on major histocompatibility complex class II gene diversity in the greater horseshoe bat

Heterogeneous pathogenic stress can shape major histocompatibility complex (MHC) diversity by influencing the functional plasticity of the immune response. Therefore, MHC diversity could reflect environmental stress, demonstrating its importance in uncovering the mechanisms of adaptive genetic variation. In this study, we combined neutral microsatellite loci, an immune-related MHC II-DRB locus, and climatic factors to unravel the mechanisms affecting the diversity and genetic differentiation of MHC genes in the greater horseshoe bat (Rhinolophus ferrumequinum), a species with a wide geographical distribution that has three distinct genetic lineages in China. First, increased genetic differentiation at the MHC locus among populations compared using microsatellites indicated diversifying selection. Second, the genetic differentiation of MHC and microsatellites were significantly correlated, suggesting that demographic processes exist. However, MHC genetic differentiation was significantly correlated with geographical distance among populations, even after controlling for the neutral markers, suggesting a major effect of selection. Third, although the MHC genetic differentiation was larger than that for microsatellites, there was no significant difference in the genetic differentiation between the two markers among genetic lineages, indicating the effect of balancing selection. Fourth, combined with climatic factors, MHC diversity and supertypes showed significant correlations with temperature and precipitation, but not with the phylogeographic structure of R. ferrumequinum, suggesting an effect of local adaptation driven by climate on MHC diversity. Moreover, the number of MHC supertypes varied between populations and lineages, suggesting regional characteristics and support for local adaptation. Taken together, the results of our study provide insights into the adaptive evolutionary driving forces at different geographic scales in R. ferrumequinum. In addition, climate factors may have played a vital role in driving adaptive evolution in this species.

Assessing the causes and consequences of gut mycobiome variation in a wild population of the Seychelles warbler

BACKGROUND

Considerable research has focussed on the importance of bacterial communities within the vertebrate gut microbiome (GM). However, studies investigating the significance of other microbial kingdoms, such as fungi, are notably lacking, despite their potential to influence host processes. Here, we characterise the fungal GM of individuals living in a natural population of Seychelles warblers (Acrocephalus sechellensis). We evaluate the extent to which fungal GM structure is shaped by environment and host factors, including genome-wide heterozygosity and variation at key immune genes (major histocompatibility complex (MHC) and Toll-like receptor (TLR)). Importantly, we also explore the relationship between fungal GM differences and subsequent host survival. To our knowledge, this is the first time that the genetic drivers and fitness consequences of fungal GM variation have been characterised for a wild vertebrate population.

RESULTS

Environmental factors, including season and territory quality, explain the largest proportion of variance in the fungal GM. In contrast, neither host age, sex, genome-wide heterozygosity, nor TLR3 genotype was associated with fungal GM differences in Seychelles warblers. However, the presence of four MHC-I alleles and one MHC-II allele was associated with changes in fungal GM alpha diversity. Changes in fungal richness ranged from between 1 and 10 sequencing variants lost or gained; in some cases, this accounted for 20% of the fungal variants carried by an individual. In addition to this, overall MHC-I allelic diversity was associated with small, but potentially important, changes in fungal GM composition. This is evidenced by the fact that fungal GM composition differed between individuals that survived or died within 7 months of being sampled.

CONCLUSIONS

Our results suggest that environmental factors play a primary role in shaping the fungal GM, but that components of the host immune system-specifically the MHC-may also contribute to the variation in fungal communities across individuals within wild populations. Furthermore, variation in the fungal GM can be associated with differential survival in the wild. Further work is needed to establish the causality of such relationships and, thus, the extent to which components of the GM may impact host evolution. Video Abstract.

Over 25 Years of Partnering to Conserve Chiricahua Leopard Frogs (Rana chiricahuensis) in Arizona, Combining Ex Situ and In Situ Strategies

The Phoenix Zoo has partnered with US Fish and Wildlife Service, Arizona Game and Fish Department, US Forest Service, and other organizations for more than 25 years to help recover Chiricahua leopard frogs (Rana [=Lithobates] chiricahuensis) in Arizona, USA. This federally threatened species faces declines due to habitat loss and degradation, long-term drought, disease, and invasive species. Over 26,000 larvae, froglets, and adults, as well as 26 egg masses produced by adults held at the Phoenix Zoo have been released to the wild, augmenting and/or re-establishing wild populations. Chiricahua leopard frog-occupied sites in Arizona have increased from 38 in 2007, when the species’ recovery plan was published, to a high of 155 in the last five years, as a result of ex situ and in situ conservation efforts. As one of the longest-running programs of its kind in the United States, communication among partners has been key to sustaining it. Recovery strategies and complex decisions are made as a team and we have worked through numerous management challenges together. Though Chiricahua leopard frogs still face significant threats and a long road to recovery, this program serves as a strong example of the positive effects of conservation partnerships for native wildlife.

Habitat fragmentation in the Brazilian Atlantic Forest is associated with erosion of frog immunogenetic diversity and increased fungal infections

Habitat fragmentation and infectious diseases threaten wildlife globally, but the interactions of these threats are poorly understood. For instance, while habitat fragmentation can impact genetic diversity at neutral loci, the impacts on disease-relevant loci are less well-studied. We examined the effects of habitat fragmentation in Brazil's Atlantic Forest on amphibian genetic diversity at an immune locus related to antigen presentation and detection (MHC IIB Exon 2). We used a custom high-throughput assay to sequence a fragment of MHC IIB and quantified Batrachochytrium dendrobatidis (Bd) infections in six frog species in two Atlantic Forest regions. Habitat fragmentation was associated with genetic erosion at MHC IIB Exon 2. This erosion was most severe in forest specialists. Significant Bd infections were detected only in one Atlantic Forest region, potentially due to relatively higher elevation. In this region, forest specialists showed an increase in both Bd prevalence and infection loads in fragmented habitats. Reduced population-level MHC IIB diversity was associated with increased Bd infection risk. On the individual level, MHC IIB heterozygotes exhibited a trend toward reduced Bd infection risk, although this was marginally non-significant. Our results suggest that habitat fragmentation increases Bd infection susceptibility in amphibians, mediated at least in part through erosion of immunogenetic diversity. Our findings have implications for management of fragmented populations in the face of emerging infectious diseases.

Challenging a host-pathogen paradigm: Susceptibility to chytridiomycosis is decoupled from genetic erosion

The putatively positive association between host genetic diversity and the ability to defend against pathogens has long attracted the attention of evolutionary biologists. Chytridiomycosis, a disease caused by the chytrid fungus Batrachochytrium dendrobatidis (Bd), has emerged in recent decades as a cause of dramatic declines and extinctions across the amphibian clade. Bd susceptibility can vary widely across populations of the same species, but the relationship between standing genetic diversity and susceptibility has remained notably underexplored so far. Here, we focus on a putatively Bd-naive system of two mainland and two island populations of the common toad (Bufo bufo) at the edge of the species' range and use controlled infection experiments and dd-RAD sequencing of >10 000 SNPs across 95 individuals to characterize the role of host population identity, genetic variation and individual body mass in mediating host response to the pathogen. We found strong genetic differentiation between populations and marked variation in their susceptibility to Bd. This variation was not, however, governed by isolation-mediated genetic erosion, and individual heterozygosity was even found to be negatively correlated with survival. Individual survival during infection experiments was strongly positively related to body mass, which itself was unrelated to population of origin or heterozygosity. Our findings underscore the general importance of context-dependency when assessing the role of host genetic variation for the ability of defence against pathogens.

Adaptive evolution of major histocompatibility complex class I immune genes and disease associations in coastal juvenile sea turtles

Characterizing polymorphism at the major histocompatibility complex (MHC) genes is key to understanding the vertebrate immune response to disease. Despite being globally afflicted by the infectious tumour disease fibropapillomatosis (FP), immunogenetic variation in sea turtles is minimally explored. We sequenced the α 1 peptide-binding region of MHC class I genes (162 bp) from 268 juvenile green (Chelonia mydas) and 88 loggerhead (Caretta caretta) sea turtles in Florida, USA. We recovered extensive variation (116 alleles) and trans-species polymorphism. Supertyping analysis uncovered three functional MHC supertypes corresponding to the three well-supported clades in the phylogeny. We found significant evidence of positive selection at seven amino acid sites in the class I exon. Random forest modelling and risk ratio analysis of Ch. mydas alleles uncovered one allele weakly associated with smooth FP tumour texture, which may be associated with disease outcome. Our study represents the first characterization of MHC class I diversity in Ch. mydas and the largest sample of sea turtles used to date in any study of adaptive genetic variation, revealing tremendous genetic variation and high adaptive potential to viral pathogen threats. The novel associations we identified between MHC diversity and FP outcomes in sea turtles further highlight the importance of evaluating genetic predictors of disease, including MHC and other functional markers.

Predicting fungal infection rate and severity with skin-associated microbial communities on amphibians

Pathogen success (risk and severity) is influenced by host-associated microbiota, but the degree to which variation in microbial community traits predict future infection presence/absence (risk) and load (severity) for the host is unknown. We conducted a time-series experiment by sampling the skin-associated bacterial communities of five amphibian species before and after exposure to the fungal pathogen, Batrachochytrium dendrobaditis (Bd). We ask whether microbial community traits are predictors of, or are affected by, Bd infection risk and intensity. Our results show that richness of putative Bd-inhibitory bacteria strongly predicts infection risk, while the proportion of putative Bd-inhibitory bacteria predicts future infection intensity. Variation in microbial community composition is high across time and individual, and bacterial prevalence is low. Our findings demonstrate how ecological community traits of host-associated microbiota may be used to predict infection risk by pathogenic microbes.

Bacterial biofilm thickness and fungal-inhibitory bacterial richness both prevent establishment of the amphibian fungal pathogen, Batrachochytrium dendrobatidis

Host-associated microbial biofilms can provide protection against pathogen establishment. In many host-microbe symbioses (including, but not limited to: humans, plants, insects, and amphibians), there is a correlation between host-associated microbial diversity and pathogen infection risk. Diversity may prevent infection by pathogens through sampling effects and niche complementarity- but an alternative hypothesis may be that microbial biomass is confounded with diversity, and that host-associated biofilms are deterring pathogen establishment through space pre-emption. In this study, we use the amphibian system as a model for host-microbe-pathogen interactions to ask two questions: (1) is bacterial richness confounded with biofilm thickness or cell density, and (2) to what extent does biofilm thickness, cell density, and bacterial richness each deter the establishment of the amphibian fungal pathogen, Batrachochytrium dendrobatidis (Bd)? To answer these questions, we built a custom biofilm microcosm that mimics the host-environment interface by allowing nutrients to diffuse out of a fine-pore biofilm scaffolding. This created a competitive environment in which bacteria and the fungal pathogen compete for colonization space. We then challenged bacterial biofilms ranging in community richness, biofilm thickness, bacterial cell density, and Bd-inhibitory metabolite production with live Bd zoospores to determine how Bd establishment success on membranes vary. We found that biofilm thickness and Bd-inhibitory isolate richness work in complement to reduce Bd establishment success. This work underscores that physical aspects of biofilm communities can play a large role in pathogen inhibition and in many studies, these traits are not studied. IMPORTANCE Our finding highlights the fact that diversity, as measured through 16S rDNA sequencing, may obscure the true mechanisms behind microbe-mediated pathogen defence, and that physical space occupation by biofilm-forming symbionts may significantly contribute to pathogen protection. These findings have implications across a wide range of host-microbe systems, since 16S rDNA sequencing is a standard tool used across many microbial systems. Further, our results are potentially relevant to many host-pathogen systems, since host-associated bacterial biofilms are ubiquitous.

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.

Virulence and Pathogenicity of Chytrid Fungi Causing Amphibian Extinctions

Ancient enzootic associations between wildlife and their infections allow evolution to innovate mechanisms of pathogenicity that are counterbalanced by host responses. However, erosion of barriers to pathogen dispersal by globalization leads to the infection of hosts that have not evolved effective resistance and the emergence of highly virulent infections. Global amphibian declines driven by the rise of chytrid fungi and chytridiomycosis are emblematic of emerging infections. Here, we review how modern biological methods have been used to understand the adaptations and counteradaptations that these fungi and their amphibian hosts have evolved. We explore the interplay of biotic and abiotic factors that modify the virulence of these infections and dissect the complexity of this disease system. We highlight progress that has led to insights into how we might in the future lessen the impact of these emerging infections. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

Host defense or parasite cue: Skin secretions mediate interactions between amphibians and their parasites

Amphibian skin secretions (substances produced by the amphibian plus microbiota) plausibly act as a first line of defense against parasite/pathogen attack, but may also provide chemical cues for pathogens. To clarify the role of skin secretions in host-parasite interactions, we conducted experiments using cane toads (Rhinella marina) and their lungworms (Rhabdias pseudosphaerocephala) from the range-core and invasion-front of the introduced anurans' range in Australia. Depending on the geographical area, toad skin secretions can reduce the longevity and infection success of parasite larvae, or attract lungworm larvae and enhance their infection success. These striking differences between the two regions were due both to differential responses of the larvae, and differential effects of the skin secretions. Our data suggest that skin secretions play an important role in host-parasite interactions in anurans, and that the arms race between a host and parasite can rapidly generate spatial variation in critical features of that interaction.

Host-multiparasite interactions in amphibians: a review

Parasites, including viruses, bacteria, fungi, protists, helminths, and arthropods, are ubiquitous in the animal kingdom. Consequently, hosts are frequently infected with more than one parasite species simultaneously. The assessment of such co-infections is of fundamental importance for disease ecology, but relevant studies involving non-domesticated animals have remained scarce. Many amphibians are in decline, and they generally have a highly diverse parasitic fauna. Here we review the literature reporting on field surveys, veterinary case studies, and laboratory experiments on co-infections in amphibians, and we summarize what is known about within-host interactions among parasites, which environmental and intrinsic factors influence the outcomes of these interactions, and what effects co-infections have on hosts. The available literature is piecemeal, and patterns are highly diverse, so that identifying general trends that would fit most host–multiparasite systems in amphibians is difficult. Several examples of additive, antagonistic, neutral, and synergistic effects among different parasites are known, but whether members of some higher taxa usually outcompete and override the effects of others remains unclear. The arrival order of different parasites and the time lag between exposures appear in many cases to fundamentally shape competition and disease progression. The first parasite to arrive can gain a marked reproductive advantage or induce cross-reaction immunity, but by disrupting the skin and associated defences (i.e., skin secretions, skin microbiome) and by immunosuppression, it can also pave the way for subsequent infections. Although there are exceptions, detrimental effects to the host are generally aggravated with increasing numbers of co-infecting parasite species. Finally, because amphibians are ectothermic animals, temperature appears to be the most critical environmental factor that affects co-infections, partly via its influence on amphibian immune function, partly due to its direct effect on the survival and growth of parasites. Besides their importance for our understanding of ecological patterns and processes, detailed knowledge about co-infections is also crucial for the design and implementation of effective wildlife disease management, so that studies concentrating on the identified gaps in our understanding represent rewarding research avenues.

Chytridiomycosis in Asian Amphibians, a Global Resource for Batrachochytrium dendrobatidis (Bd) Research

Chytridiomycosis is an emerging infectious disease affecting amphibians globally and it is caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd). Chytridiomycosis has caused dramatic declines and even extinctions in wild amphibian populations in Europe, Australia, Central and North America. Spanning over two and a half decades, extensive research has led to discovery of epizootic and enzootic lineages of this pathogen. However, the Bd–amphibian system had garnered less attention in Asia until recently when an ancestral Bd lineage was identified in the Korean peninsula. Amphibians co-exist with the pathogen in Asia, only sub-lethal effects have been documented on hosts. Such regions are ‘coldspots’ of infection and are an important resource to understand the dynamics between the enzootic pathogen—Bd and its obligate host—amphibians. Insights into the biology of infection have provided new knowledge on the multi-faceted interaction of Bd in a hyperdiverse Asian amphibian community. We present the findings and highlight the knowledge gap that exists, and propose the ways to bridge them. We emphasize that chytridiomycosis in Asia is an important wildlife disease and it needs focussed research, as it is a dynamic front of pathogen diversity and virulence.

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.

The inflated significance of neutral genetic diversity in conservation genetics

The current rate of species extinction is rapidly approaching unprecedented highs, and life on Earth presently faces a sixth mass extinction event driven by anthropogenic activity, climate change, and ecological collapse. The field of conservation genetics aims at preserving species by using their levels of genetic diversity, usually measured as neutral genome-wide diversity, as a barometer for evaluating population health and extinction risk. A fundamental assumption is that higher levels of genetic diversity lead to an increase in fitness and long-term survival of a species. Here, we argue against the perceived importance of neutral genetic diversity for the conservation of wild populations and species. We demonstrate that no simple general relationship exists between neutral genetic diversity and the risk of species extinction. Instead, a better understanding of the properties of functional genetic diversity, demographic history, and ecological relationships is necessary for developing and implementing effective conservation genetic strategies.

Genotyping and differential bacterial inhibition of Batrachochytrium dendrobatidis in threatened amphibians in Costa Rica

Amphibians have declined around the world in recent years, in parallel with the emergence of an epidermal disease called chytridiomycosis, caused by the chytrid fungus Batrachochytrium dendrobatidis (Bd). This disease has been associated with mass mortality in amphibians worldwide, including in Costa Rica, and Bd is considered an important contributor to the disappearance of this group of vertebrates. While many species are susceptible to the disease, others show tolerance and manage to survive infection with the pathogen. We evaluated the pathogen Bd circulating in Costa Rica and the capacity of amphibian skin bacteria to inhibit the growth of the pathogen in vitro. We isolated and characterized - genetically and morphologically - several Bd isolates from areas with declining populations of amphibians. We determined that the circulating chytrid fungus in Costa Rica belongs to the virulent strain Bd-GPL-2, which has been related to massive amphibian deaths worldwide; however, the isolates obtained showed genetic and morphological variation. Furthermore, we isolated epidermal bacteria from 12 amphibian species of surviving populations, some in danger of extinction, and evaluated their inhibitory activity against the collection of chytrid isolates. Through bioassays we confirmed the presence of chytrid-inhibitory bacterial genera in Costa Rican amphibians. However, we observed that the inhibition varied between different isolates of the same bacterial genus, and each bacterial isolation inhibited fungal isolation differently. In total, 14 bacterial isolates belonging to the genera Stenotrophomonas, Streptomyces, Enterobacter, Pseudomonas and Klebsiella showed inhibitory activity against all Bd isolates. Given the observed variation both in the pathogen and in the bacterial inhibition capacity, it is highly relevant to include local isolates and to consider the origin of the microorganisms when performing in vivo infection tests aimed at developing and implementing mitigation strategies for chytridiomycosis.

Sarcoptic mange severity is associated with reduced genomic variation and evidence of selection in Yellowstone National Park wolves (Canis lupus)

Population genetic theory posits that molecular variation buffers against disease risk. Although this "monoculture effect" is well supported in agricultural settings, its applicability to wildlife populations remains in question. In the present study, we examined the genomics underlying individual-level disease severity and population-level consequences of sarcoptic mange infection in a wild population of canids. Using gray wolves (Canis lupus) reintroduced to Yellowstone National Park (YNP) as our focal system, we leveraged 25 years of observational data and biobanked blood and tissue to genotype 76,859 loci in over 400 wolves. At the individual level, we reported an inverse relationship between host genomic variation and infection severity. We additionally identified 410 loci significantly associated with mange severity, with annotations related to inflammation, immunity, and skin barrier integrity and disorders. We contextualized results within environmental, demographic, and behavioral variables, and confirmed that genetic variation was predictive of infection severity. At the population level, we reported decreased genome-wide variation since the initial gray wolf reintroduction event and identified evidence of selection acting against alleles associated with mange infection severity. We concluded that genomic variation plays an important role in disease severity in YNP wolves. This role scales from individual to population levels, and includes patterns of genome-wide variation in support of the monoculture effect and specific loci associated with the complex mange phenotype. Results yielded system-specific insights, while also highlighting the relevance of genomic analyses to wildlife disease ecology, evolution, and conservation.

Why disease ecology needs life-history theory: a host perspective

When facing an emerging infectious disease of conservation concern, we often have little information on the nature of the host-parasite interaction to inform management decisions. However, it is becoming increasingly clear that the life-history strategies of host species can be predictive of individual- and population-level responses to infectious disease, even without detailed knowledge on the specifics of the host-parasite interaction. Here, we argue that a deeper integration of life-history theory into disease ecology is timely and necessary to improve our capacity to understand, predict and mitigate the impact of endemic and emerging infectious diseases in wild populations. Using wild vertebrates as an example, we show that host life-history characteristics influence host responses to parasitism at different levels of organisation, from individuals to communities. We also highlight knowledge gaps and future directions for the study of life-history and host responses to parasitism. We conclude by illustrating how this theoretical insight can inform the monitoring and control of infectious diseases in wildlife.

Chytridiomycosis-induced mortality in a threatened anuran

Effectively planning conservation introductions involves assessing the suitability of both donor and recipient populations, including the landscape of disease risk. Chytridiomycosis, caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), has caused extensive amphibian declines globally and may hamper reintroduction attempts. To determine Bd dynamics in potential source populations for conservation translocations of the threatened California red-legged frog (Rana draytonii) to Yosemite National Park, we conducted Bd sampling in two populations in the foothills of the Sierra Nevada Mountains, California, U.S.A. At one of two sites, we observed lethally high Bd loads in early post-metamorphic life stages and confirmed one chytridiomycosis-induced mortality, the first such report for this species. These results informed source population site selection for subsequent R. draytonii conservation translocations. Conservation efforts aimed at establishing new populations of R. draytonii in a landscape where Bd is ubiquitous can benefit from an improved understanding of risk through disease monitoring and ex situ infection studies.

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.

Population demographic history and adaptability of the vulnerable Lolokou Sucker Frog

Amphibians are experiencing worldwide declines due to increasing anthropogenetic disturbances. However, the genetic variability and hence adaptability are still unknown for most frogs. We integrated the mitochondrial (ND2 gene), nuclear (TYR gene) and major histocompatibility complex (MHC) loci, to clarify the demographic patterns and immune-gene diversity of the Lolokou Sucker Frog (Amolops loloensis). Demographic analysis of the ND2 and TYR genes suggested that the Lolokou Sucker Frog experienced a population expansion within the last 10,000 years. High MHC diversity was detected, which has likely resulted from positive selection, indicating the current diversity bodes well for the species' adaptive potential to pathogenic challenges. These findings broaden our knowledge on the population history and evolution adaptation of the reclusive torrent frog, and conservation implications are provided.

Major histocompatibility complex variation is similar in little brown bats before and after white-nose syndrome outbreak

White-nose syndrome (WNS), caused by the fungal pathogen Pseudogymnoascus destructans (Pd), has driven alarming declines in North American hibernating bats, such as little brown bat (Myotis lucifugus). During hibernation, infected little brown bats are able to initiate anti-Pd immune responses, indicating pathogen-mediated selection on the major histocompatibility complex (MHC) genes. However, such immune responses may not be protective as they interrupt torpor, elevate energy costs, and potentially lead to higher mortality rates. To assess whether WNS drives selection on MHC genes, we compared the MHC DRB gene in little brown bats pre- (Wisconsin) and post- (Michigan, New York, Vermont, and Pennsylvania) WNS (detection spanning 2014-2015). We genotyped 131 individuals and found 45 nucleotide alleles (27 amino acid alleles) indicating a maximum of 3 loci (1-5 alleles per individual). We observed high allelic admixture and a lack of genetic differentiation both among sampling sites and between pre- and post-WNS populations, indicating no signal of selection on MHC genes. However, post-WNS populations exhibited decreased allelic richness, reflecting effects from bottleneck and drift following rapid population declines. We propose that mechanisms other than adaptive immunity are more likely driving current persistence of little brown bats in affected regions.

Tracking Batrachochytrium dendrobatidis Infection Across the Globe

Infection records of Batrachochytrium dendrobatidis (Bd), a pathogen that has devastated amphibian populations worldwide, have rapidly increased since the pathogen's discovery. Dealing with so many records makes it difficult to (a) know where, when and in which species infections have been detected, (b) understand how widespread and pervasive Bd is and (c) prioritize study and management areas. We conducted a systematic review of papers and compiled a database with Bd infection records. Our dataset covers 71 amphibian families and 119 countries. The data revealed how widespread and adaptable Bd is, being able to infect over 50% of all tested amphibian species, with over 1000 confirmed host species and being present in 86 countries. The distribution of infected species is uneven among and within countries. Areas where the distributions of many infected species overlap are readily visible; these are regions where Bd likely develops well. Conversely, areas where the distributions of species that tested negative overlap, such as the Atlantic Coast in the USA, suggest the presence of Bd refuges. Finally, we report how the number of tested and infected species has changed through time, and provide a list of oldest detection records per country.

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.

Evolutionary principles guiding amphibian conservation

The Anthropocene has witnessed catastrophic amphibian declines across the globe. A multitude of new, primarily human-induced drivers of decline may lead to extinction, but can also push species onto novel evolutionary trajectories. If these are recognized by amphibian biologists, they can be engaged in conservation actions. Here, we summarize how principles stemming from evolutionary concepts have been applied for conservation purposes, and address emerging ideas at the vanguard of amphibian conservation science. In particular, we examine the consequences of increased drift and inbreeding in small populations and their implications for practical conservation. We then review studies of connectivity between populations at the landscape level, which have emphasized the limiting influence of anthropogenic structures and degraded habitat on genetic cohesion. The rapid pace of environmental changes leads to the central question of whether amphibian populations can cope either by adapting to new conditions or by shifting their ranges. We gloomily conclude that extinction seems far more likely than adaptation or range shifts for most species. That said, conservation strategies employing evolutionary principles, such as selective breeding, introduction of adaptive variants through translocations, ecosystem interventions aimed at decreasing phenotype-environment mismatch, or genetic engineering, may effectively counter amphibian decline in some areas or for some species. The spread of invasive species and infectious diseases has often had disastrous consequences, but has also provided some premier examples of rapid evolution with conservation implications. Much can be done in terms of setting aside valuable amphibian habitat that should encompass both natural and agricultural areas, as well as designing protected areas to maximize the phylogenetic and functional diversity of the amphibian community. We conclude that an explicit consideration and application of evolutionary principles, although certainly not a silver bullet, should increase effectiveness of amphibian conservation in both the short and long term.

Investigating the potential use of an ionic liquid (1-Butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide) as an anti-fungal treatment against the amphibian chytrid fungus, Batrachochytrium dendrobatidis

The disease chytridiomycosis, caused by the pathogenic chytrid fungus, Batrachochytrium dendrobatidis (Bd), has contributed to global amphibian declines. Bd infects the keratinized epidermal tissue in amphibians and causes hyperkeratosis and excessive skin shedding. In individuals of susceptible species, the regulatory function of the amphibian’s skin is disrupted resulting in an electrolyte depletion, osmotic imbalance, and eventually death. Safe and effective treatments for chytridiomycosis are urgently needed to control chytrid fungal infections and stabilize populations of endangered amphibian species in captivity and in the wild. Currently, the most widely used anti-Bd treatment is itraconazole. Preparations of itraconazole formulated for amphibian use has proved effective, but treatment involves short baths over seven to ten days, a process which is logistically challenging, stressful, and causes long-term health effects. Here, we explore a novel anti-fungal therapeutic using a single application of the ionic liquid, 1-Butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide (BMP-NTf2), for the treatment of chytridiomycosis. BMP-NTf2 was found be effective at killing Bd in vitro at low concentrations (1:1000 dilution). We tested BMP-NTf2 in vivo on two amphibian species, one that is relatively tolerant of chytridiomycosis (Pseudacris regilla) and one that is highly susceptible (Dendrobates tinctorius). A toxicity trial revealed a surprising interaction between Bd infection status and the impact of BMP-NTf2 on D. tinctorius survival. Uninfected D. tinctorius tolerated BMP-NTf2 (mean ± SE; 96.01 ± 9.00 μl/g), such that only 1 out of 30 frogs died following treatment (at a dose of 156.95 μL/g), whereas, a lower dose (mean ± SE; 97.45 ± 3.52 μL/g) was not tolerated by Bd-infected D. tinctorius, where 15 of 23 frogs died shortly upon BMP-NTf2 application. Those that tolerated the BMP-NTf2 application did not exhibit Bd clearance. Thus, BMP-NTf2 application, under the conditions tested here, is not a suitable option for clearing Bd infection in D. tinctorius. However, different results were obtained for P. regilla. Two topical applications of BMP-NTf2 on Bd-infected P. regilla (using a lower BMP-NTf2 dose than on D. tinctorius, mean ± SE; 9.42 ± 1.43 μL/g) reduced Bd growth, although the effect was lower than that obtained by daily doses of itracanozole (50% frogs exhibited complete clearance on day 16 vs. 100% for itracanozole). Our findings suggest that BMP-NTf2 has the potential to treat Bd infection, however the effect depends on several parameters. Further optimization of dose and schedule are needed before BMP-NTf2 can be considered as a safe and effective alternative to more conventional antifungal agents, such as itraconazole.

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.