Genomic Correlates of Virulence Attenuation in the Deadly Amphibian Chytrid Fungus, Batrachochytrium dendrobatidis

Advances in Managing Chytridiomycosis for Australian Frogs: Gradarius Firmus Victoria

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

Coinfection with chytrid genotypes drives divergent infection dynamics reflecting regional distribution patterns

By altering the abundance, diversity, and distribution of species-and their pathogens-globalization may inadvertently select for more virulent pathogens. In Brazil's Atlantic Forest, a hotspot of amphibian biodiversity, the global amphibian trade has facilitated the co-occurrence of previously isolated enzootic and panzootic lineages of the pathogenic amphibian-chytrid (Batrachochytrium dendrobatidis, 'Bd') and generated new virulent recombinant genotypes ('hybrids'). Epidemiological data indicate that amphibian declines are most severe in hybrid zones, suggesting that coinfections are causing more severe infections or selecting for higher virulence. We investigated how coinfections involving these genotypes shapes virulence and transmission. Overall, coinfection favored the more virulent and competitively superior panzootic genotype, despite dampening its transmission potential and overall virulence. However, for the least virulent and least competitive genotype, coinfection increased both overall virulence and transmission. Thus, by integrating experimental and epidemiological data, our results provide mechanistic insight into how globalization can select for, and propel, the emergence of introduced hypervirulent lineages, such as the globally distributed panzootic lineage of Bd.

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.

Major Histocompatibility Complex Variation and Haplotype Associated Survival in Response to Experimental Infection of Two Bd-GPL Strains Along a Latitudinal Gradient

While both innate and adaptive immune system mechanisms have been implicated in resistance against the chytrid fungus Batrachochytrium dendrobatidis (Bd), studies on the role of specific MHC haplotypes on Bd infection are rare. Here, we studied variation in MHC Class IIB loci in the common toad Bufo bufo along a latitudinal gradient across Sweden. In general, Swedish toad populations had few MHC Class IIB haplotypes and MHC diversity declined from south (13 haplotypes) to the north (four haplotypes). The low diversity may compromise the ability of northern populations to fight emerging disease, such as Bd. In a laboratory experiment, we infected newly metamorphosed toads with two strains of the Global Pandemic Lineage of the fungus (Bd-GPL) and compared survival with sham controls. Bd-infected toads had lower survival compared to controls. Moreover, survival was dependent on the Bd-strain and northern toads had lower Bd-mediated survival than southern individuals. MHC diversity was lower in northern toads. All northern experimental animals were monomorphic for a single MHC haplotype, whereas we found seven different haplotypes in southern experimental animals. In southern toads, survival was dependent on both Bd-strain and MHC haplotype suggesting differential infection dynamics depending on both Bd-strain and host immune system characteristics.

Host species is linked to pathogen genotype for the amphibian chytrid fungus (Batrachochytrium dendrobatidis)

Host-pathogen specificity can arise from certain selective environments mediated by both the host and pathogen. Therefore, understanding the degree to which host species identity is correlated with pathogen genotype can help reveal historical host-pathogen dynamics. One animal disease of particular concern is chytridiomycosis, typically caused by the global panzootic lineage of the amphibian chytrid fungus (Batrachochytrium dendrobatidis, Bd), termed Bd-GPL. This pathogen lineage has caused devastating declines in amphibian communities around the world. However, the site of origin for the common ancestor of modern Bd-GPL and the fine-scale transmission dynamics of this lineage have remained a mystery. This is especially the case in North America where Bd-GPL is widespread, but disease outbreaks occur sporadically. Herein, we use Bd genetic data collected throughout the United States from amphibian skin swabs and cultured isolate samples to investigate Bd genetic patterns. We highlight two case studies in Pennsylvania and Nevada where Bd-GPL genotypes are strongly correlated with host species identity. Specifically, in some localities bullfrogs (Rana catesbeiana) are infected with Bd-GPL lineages that are distinct from those infecting other sympatric amphibian species. Overall, we reveal a previously unknown association of Bd genotype with host species and identify the eastern United States as a Bd diversity hotspot and potential site of origin for Bd-GPL.

Social Behavior, Community Composition, Pathogen Strain, and Host Symbionts Influence Fungal Disease Dynamics in Salamanders

The emerging fungal pathogen, Batrachochytrium dendrobatidis (Bd), which can cause a fatal disease called chytridiomycosis, is implicated in the collapse of hundreds of host amphibian species. We describe chytridiomycosis dynamics in two co-occurring terrestrial salamander species, the Santa Lucia Mountains slender salamander, Batrachoseps luciae, and the arboreal salamander, Aneides lugubris. We (1) conduct a retrospective Bd-infection survey of specimens collected over the last century, (2) estimate present-day Bd infections in wild populations, (3) use generalized linear models (GLM) to identify biotic and abiotic correlates of infection risk, (4) investigate susceptibility of hosts exposed to Bd in laboratory trials, and (5) examine the ability of host skin bacteria to inhibit Bd in culture. Our historical survey of 2,866 specimens revealed that for most of the early 20th century (~1920–1969), Bd was not detected in either species. By the 1990s the proportion of infected specimens was 29 and 17% (B. luciae and A. lugubris, respectively), and in the 2010s it was 10 and 17%. This was similar to the number of infected samples from contemporary populations (2014–2015) at 10 and 18%. We found that both hosts experience signs of chytridiomycosis and suffered high Bd-caused mortality (88 and 71% for B. luciae and A. lugubris, respectively). Our GLM revealed that Bd-infection probability was positively correlated with intraspecific group size and proximity to heterospecifics but not to abiotic factors such as precipitation, minimum temperature, maximum temperature, mean temperature, and elevation, or to the size of the hosts. Finally, we found that both host species contain symbiotic skin-bacteria that inhibit growth of Bd in laboratory trials. Our results provide new evidence consistent with other studies showing a relatively recent Bd invasion of amphibian host populations in western North America and suggest that the spread of the pathogen may be enabled both through conspecific and heterospecific host interactions. Our results suggest that wildlife disease studies should assess host-pathogen dynamics that consider the interactions and effects of multiple hosts, as well as the historical context of pathogen invasion, establishment, and epizootic to enzootic transitions to better understand and predict disease dynamics.

Diversity, multifaceted evolution, and facultative saprotrophism in the European Batrachochytrium salamandrivorans epidemic

While emerging fungi threaten global biodiversity, the paucity of fungal genome assemblies impedes thoroughly characterizing epidemics and developing effective mitigation strategies. Here, we generate de novo genomic assemblies for six outbreaks of the emerging pathogen Batrachochytrium salamandrivorans (Bsal). We reveal the European epidemic currently damaging amphibian populations to comprise multiple, highly divergent lineages demonstrating isolate-specific adaptations and metabolic capacities. In particular, we show extensive gene family expansions and acquisitions, through a variety of evolutionary mechanisms, and an isolate-specific saprotrophic lifecycle. This finding both explains the chytrid's ability to divorce transmission from host density, producing Bsal's enigmatic host population declines, and is a key consideration in developing successful mitigation measures.

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.

Presence of low virulence chytrid fungi could protect European amphibians from more deadly strains

Wildlife diseases are contributing to the current Earth’s sixth mass extinction; one disease, chytridiomycosis, has caused mass amphibian die-offs. While global spread of a hypervirulent lineage of the fungus Batrachochytrium dendrobatidis (BdGPL) causes unprecedented loss of vertebrate diversity by decimating amphibian populations, its impact on amphibian communities is highly variable across regions. Here, we combine field data with in vitro and in vivo trials that demonstrate the presence of a markedly diverse variety of low virulence isolates of BdGPL in northern European amphibian communities. Pre-exposure to some of these low virulence isolates protects against disease following subsequent exposure to highly virulent BdGPL in midwife toads (Alytes obstetricans) and alters infection dynamics of its sister species B. salamandrivorans in newts (Triturus marmoratus), but not in salamanders (Salamandra salamandra). The key role of pathogen virulence in the complex host-pathogen-environment interaction supports efforts to limit pathogen pollution in a globalized world.

The pathogen Batrachochytrium dendrobatidis (BD) associated with widespread amphibian declines is present in Europe but has not consistently caused disease-induced declines in that region. Here, the authors suggest that an endemic strain of BD with low virulence may protect the hosts upon co-infection with more virulent strains.

Experimental methodologies can affect pathogenicity of Batrachochytrium salamandrivorans infections

Controlled experiments are one approach to understanding the pathogenicity of etiologic agents to susceptible hosts. The recently discovered fungal pathogen, Batrachochytrium salamandrivorans (Bsal), has resulted in a surge of experimental investigations because of its potential to impact global salamander biodiversity. However, variation in experimental methodologies could thwart knowledge advancement by introducing confounding factors that make comparisons difficult among studies. Thus, our objective was to evaluate if variation in experimental methods changed inferences made on the pathogenicity of Bsal. We tested whether passage duration of Bsal culture, exposure method of the host to Bsal (water bath vs. skin inoculation), Bsal culturing method (liquid vs. plated), host husbandry conditions (aquatic vs. terrestrial), and skin swabbing frequency influenced diseased-induced mortality in a susceptible host species, the eastern newt (Notophthalmus viridescens). We found that disease-induced mortality was faster for eastern newts when exposed to a low passage isolate, when newts were housed in terrestrial environments, and if exposure to zoospores occurred via water bath. We did not detect differences in disease-induced mortality between culturing methods or swabbing frequencies. Our results illustrate the need to standardize methods among Bsal experiments. We provide suggestions for future Bsal experiments in the context of hypothesis testing and discuss the ecological implications of our results.

Comparative Analysis of Putative Virulence-Associated Factors of Microsporum canis Isolates from Human and Animal Patients

Microsporum canis is a zoophilic dermatophyte and the most common fungus isolated from dogs and cats worldwide. To invade skin, this pathogen uses different enzymes, which may be associated with virulence, that contribute to the fungal pathogenicity. The aim of this study is to compare the expression of enzymes that may be associated with virulence, and thermotolerance of M. canis strains isolated from dogs, cats, and humans. The in vitro expression of the enzymes keratinase, catalase, urease, hemolysin, and aspartic protease was evaluated in 52 M. canis strains recently isolated from 14 human patients, 12 dogs, 15 symptomatic, and 11 asymptomatic cats. In addition, thermotolerance was assessed by comparative analysis of fungal growth at 25 °C and 35 °C. Keratinase activity was low in 34 and moderate in 18 strains. Aspartic-protease activity was low in 7, moderate in 33, and high in 12 strains. Hemolysin activity was low in 44 and moderate in 8 strains. All strains were classified as low producers of catalase. All but three strains produced urease in vitro, with a broad range of activity. The strains presented in vitro growth at the two studied temperatures were classified as presenting low (36.5%), medium (44.3%), or high (19.2%) thermotolerance. There was no statistically significant difference in the new putative virulence-associated factors studied among the different hosts, which suggests that they may have a similar role on human, cat, and dog infection. Also, no difference was observed between strains isolated from symptomatic and asymptomatic cats. This suggests that these factors have a limited impact on the fate of feline dermatophytosis caused by M. canis.

Isolation and maintenance of Batrachochytrium salamandrivorans cultures

Discovered in 2013, the chytrid fungus Batrachochytrium salamandrivorans (Bsal) is an emerging amphibian pathogen that causes ulcerative skin lesions and multifocal erosion. A closely related pathogen, B. dendrobatidis (Bd), has devastated amphibian populations worldwide, suggesting that Bsal poses a significant threat to global salamander biodiversity. To expedite research into this emerging threat, we seek to standardize protocols across the field so that results of laboratory studies are reproducible and comparable. We have collated data and experience from multiple labs to standardize culturing practices of Bsal. Here we outline common culture practices including a medium for standardized Bsal growth, standard culturing protocols, and a method for isolating Bsal from infected tissue.

Chytrid fungi and global amphibian declines

Discovering that chytrid fungi cause chytridiomycosis in amphibians represented a paradigm shift in our understanding of how emerging infectious diseases contribute to global patterns of biodiversity loss. In this Review we describe how the use of multidisciplinary biological approaches has been essential to pinpointing the origins of amphibian-parasitizing chytrid fungi, including Batrachochytrium dendrobatidis and Batrachochytrium salamandrivorans, as well as to timing their emergence, tracking their cycles of expansion and identifying the core mechanisms that underpin their pathogenicity. We discuss the development of the experimental methods and bioinformatics toolkits that have provided a fuller understanding of batrachochytrid biology and informed policy and control measures.

Hybridization Facilitates Adaptive Evolution in Two Major Fungal Pathogens

Hybridization is increasingly recognized as an important force impacting adaptation and evolution in many lineages of fungi. During hybridization, divergent genomes and alleles are brought together into the same cell, potentiating adaptation by increasing genomic plasticity. Here, we review hybridization in fungi by focusing on two fungal pathogens of animals. Hybridization is common between the basidiomycete yeast species Cryptococcus neoformans × Cryptococcus deneoformans, and hybrid genotypes are frequently found in both environmental and clinical settings. The two species show 10-15% nucleotide divergence at the genome level, and their hybrids are highly heterozygous. Though largely sterile and unable to mate, these hybrids can propagate asexually and generate diverse genotypes by nondisjunction, aberrant meiosis, mitotic recombination, and gene conversion. Under stress conditions, the rate of such genetic changes can increase, leading to rapid adaptation. Conversely, in hybrids formed between lineages of the chytridiomycete frog pathogen Batrachochytrium dendrobatidis (Bd), the parental genotypes are considerably less diverged (0.2% divergent). Bd hybrids are formed from crosses between lineages that rarely undergo sex. A common theme in both species is that hybrids show genome plasticity via aneuploidy or loss of heterozygosity and leverage these mechanisms as a rapid way to generate genotypic/phenotypic diversity. Some hybrids show greater fitness and survival in both virulence and virulence-associated phenotypes than parental lineages under certain conditions. These studies showcase how experimentation in model species such as Cryptococcus can be a powerful tool in elucidating the genotypic and phenotypic consequences of hybridization.

Gene expression varies within and between enzootic and epizootic lineages of Batrachochytrium dendrobatidis (Bd) in the Americas

While much research focus is paid to hypervirulent fungal lineages during emerging infectious disease outbreaks, examining enzootic pathogen isolates can be equally fruitful in delineating infection dynamics and determining pathogenesis. The fungal pathogen of amphibians, Batrachochytrium dendrobatidis (Bd), exhibits markedly different patterns of disease in natural populations, where it has caused massive amphibian declines in some regions, yet persists enzootically in others. Here we compare in vitro gene expression profiles of a panel of Bd isolates representing both the enzootic Bd-Brazil lineage, and the more recently diverged, panzootic lineage, Bd-GPL. We document significantly different lineage-specific and intralineage gene expression patterns, with Bd-Brazil upregulating genes with aspartic-type peptidase activity, and Bd-GPL upregulating CBM18 chitin-binding genes, among others. We also find pronounced intralineage variation in membrane integrity and transmembrane transport ability within our Bd-GPL isolates. Finally, we highlight unexpectedly divergent expression profiles in sympatric panzootic isolates, underscoring microgeographic functional variation in a largely clonal lineage. This variation in gene expression likely plays an important role in the relative pathogenesis and host range of Bd-Brazil and Bd-GPL isolates. Together, our results demonstrate that functional genomics approaches can provide information relevant to studies of virulence evolution within the Bd clade.

The return of the frogs: The importance of habitat refugia in maintaining diversity during a disease outbreak

Recent decades have seen the emergence and spread of numerous infectious diseases, often with severe negative consequences for wildlife populations. Nevertheless, many populations survive the initial outbreaks, and even undergo recoveries. Unfortunately, the long-term effects of these outbreaks on host population genetics are poorly understood; to increase this understanding, we examined the population genetics of two species of rainforest frogs (Litoria nannotis and Litoria serrata) that have largely recovered from a chytridiomycosis outbreak at two national parks in the Wet Tropics of northern Australia. At the wetter, northern park there was little evidence of decreased genetic diversity in either species, and all of the sampled sites had high minor allele frequencies (mean MAF = 0.230-0.235), high heterozygosity (0.318-0.325), and few monomorphic markers (1.4%-4.0%); however, some recovered L. nannotis populations had low N_e values (59.3-683.8) compared to populations that did not decline during the outbreak (1,537.4-1,756.5). At the drier, southern park, both species exhibited lower diversity (mean MAF = 0.084-0.180; heterozygosity = 0.126-0.257; monomorphic markers = 3.7%-43.5%; N_e  = 18.4-676.1). The diversity patterns in this park matched habitat patterns, with both species having higher diversity levels and fewer closely related individuals at sites with higher quality habitat. These patterns were more pronounced for L. nannotis, which has lower dispersal rates than L. serrata. These results suggest that refugia with high quality habitat are important for retaining genetic diversity during disease outbreaks, and that gene flow following disease outbreaks is important for re-establishing diversity in populations where it was reduced.

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.

Pooling skin swabs does not inhibit qPCR detection of amphibian chytrid infection

Immediate and reliable pathogen detection in large numbers of samples is essential in wildlife disease monitoring and is often realized by DNA-based techniques. Pooling samples increases processing efficiency and reduces processing costs, and has been suggested as a viable technique for quantitative PCR detection of fungal amphibian pathogens of the genus Batrachochytrium. For these fungi, this diagnostic method has been validated by in vitro set ups that provided controlled test conditions but did not take into account potential effects from amphibian skin compounds (e.g. skin secretions and Microbiota) on the approach. Some of these skin compounds are known to cause PCR inhibition in single sample applications and could lead to false negative reactions and thereby hamper pathogen detection. In this study we examined the effect of skin compounds on the pooled extraction method by swabbing individuals of seven amphibian species (one Anura and six Caudata) prior to the inoculation of the swabs with chytrid zoospores. For each species, swabs were extracted in pools of different sizes (from one to four swabs) with only one swab per pool being inoculated with zoospores. There were no significant differences regarding the ability to detect zoospores when comparing pool sizes for any species, with a tendency for more false negatives when the inoculated swab had been inoculated with a single zoospore. This study provides further in vivo evidence for the viability of the pooled extraction method for DNA-based detection of pathogens.

Diverse genotypes of the amphibian-killing fungus produce distinct phenotypes through plastic responses to temperature

Phenotypes are the target of selection and affect the ability of organisms to persist in variable environments. Phenotypes can be influenced directly by genes and/or by phenotypic plasticity. The amphibian-killing fungus Batrachochytrium dendrobatidis (Bd) has a global distribution, unusually broad host range, and high genetic diversity. Phenotypic plasticity may be an important process that allows this pathogen to infect hundreds of species in diverse environments. We quantified phenotypic variation of nine Bd genotypes from two Bd lineages (Global Pandemic Lineage [GPL] and Brazil) and a hybrid (GPL-Brazil) grown at three temperatures (12, 18 and 24°C). We measured five functional traits including two morphological traits (zoospore and zoosporangium sizes) and three life history traits (carrying capacity, time to fastest growth and exponential growth rate) in a phylogenetic framework. Temperature caused highly plastic responses within each genotype, with all Bd genotypes showing phenotypic plasticity in at least three traits. Among genotypes, Bd generally showed the same direction of plastic response to temperature: larger zoosporangia, higher carrying capacity, longer time to fastest growth and slower exponential growth at lower temperatures. The exception was zoospore size, which was highly variable. Our findings indicate that Bd genotypes have evolved novel phenotypes through plastic responses to temperature over very short timescales. High phenotypic variability likely extends to other traits and may facilitate the large host range and rapid spread of Bd.

Globally invasive genotypes of the amphibian chytrid outcompete an enzootic lineage in coinfections

Competition between genotypes is likely to be a key driver of pathogen evolution, particularly following a geographical invasion by distant strains. Theory predicts that competition between disease strains will result in the most virulent strain persisting. Despite its evolutionary implications, the role of strain competition in shaping populations remains untested for most pathogens. We experimentally investigated the in vivo competitive differences between two divergent lineages of the amphibian-killing chytrid fungus ( Batrachochytrium dendrobatidis, Bd). These Bd lineages are hypothesized to have diverged in allopatry but been recently brought back into secondary contact by human introduction. Prior studies indicate that a panzootically-distributed, global lineage of Bd was recently introduced into southern Brazil, and is competitively excluding enzootic lineages in the southern Atlantic Forest. To test for differences in competitive ability between invasive and enzootic Brazilian Bd isolates, we coinfected a model host frog system which we developed for this study ( Hymenochirus curtipes). We tracked isolate-specific zoospore production over the course of the coinfection experiment with chip-based digital PCR (dPCR). The globally invasive panzootic lineage had a competitive advantage in spore production especially during the first one to four weeks of infection, and on frogs that eventually succumbed to Bd infection. Our study provides new evidence that competitive pressure resulting from the human movement of pathogen strains can rapidly alter the genetics, community dynamics and spatial epidemiology of pathogens in the wild.

Evolutionary rescue in a host-pathogen system results in coexistence not clearance

The evolutionary rescue of host populations may prevent extinction from novel pathogens. However, the conditions that facilitate rapid evolution of hosts, in particular the population variation in host susceptibility, and the effects of host evolution in response to pathogens on population outcomes remain largely unknown. We constructed an individual-based model to determine the relationships between genetic variation in host susceptibility and population persistence in an amphibian-fungal pathogen (Batrachochytrium dendrobatidis) system. We found that host populations can rapidly evolve reduced susceptibility to a novel pathogen and that this rapid evolution led to a 71-fold increase in the likelihood of host-pathogen coexistence. However, the increased rates of coexistence came at a cost to host populations; fewer populations cleared infection, population sizes were depressed, and neutral genetic diversity was lost. Larger adult host population sizes and greater adaptive genetic variation prior to the onset of pathogen introduction led to substantially reduced rates of extinction, suggesting that populations with these characteristics should be prioritized for conservation when species are threatened by novel infectious diseases.

Opening the file drawer: Unexpected insights from a chytrid infection experiment

Infection experiments are critical for understanding wildlife disease dynamics. Although infection experiments are typically designed to reduce complexity, disease outcomes still result from complex interactions between host, pathogen, and environmental factors. Cryptic variation across factors can lead to decreased repeatability of infection experiments within and between research groups and hinder research progress. Furthermore, studies with unexpected results are often relegated to the “file drawer” and potential insights gained from these experimental outcomes are lost. Here, we report unexpected results from an infection experiment studying the response of two differentially-susceptible but related frogs (American Bullfrog Rana catesbeiana and the Mountain yellow-legged frog Rana muscosa) to the amphibian-killing chytrid fungus (Batrachochytrium dendrobatidis, Bd). Despite well-documented differences in susceptibility between species, we found no evidence for antibody-mediated immune response and no Bd-related mortality in either species. Additionally, during the study, the sham-inoculated R. catesbeiana control group became unexpectedly Bd-positive. We used a custom genotyping assay to demonstrate that the aberrantly-infected R. catesbeiana carried a Bd genotype distinct from the inoculation genotype. Thus R. catesbeiana individuals were acquired with low-intensity infections that could not be detected with qPCR. In the Bd-inoculated R. catesbeiana treatment group, the inoculated genotype appeared to out-compete the cryptic infection. Thus, our results provide insight into Bd coinfection dynamics, a phenomenon that is increasingly relevant as different pathogen strains are moved around the globe. Our experiment highlights how unexpected experimental outcomes can serve as both cautionary tales and opportunities to explore unanswered research questions. We use our results as a case study to highlight common sources of anomalous results for infection experiments. We argue that understanding these factors will aid researchers in the design, execution, and interpretation of experiments to understand wildlife disease processes.

Shifts in disease dynamics in a tropical amphibian assemblage are not due to pathogen attenuation

Infectious diseases rarely end in extinction. Yet the mechanisms that explain how epidemics subside are difficult to pinpoint. We investigated host-pathogen interactions after the emergence of a lethal fungal pathogen in a tropical amphibian assemblage. Some amphibian host species are recovering, but the pathogen is still present and is as pathogenic today as it was almost a decade ago. In addition, some species have defenses that are more effective now than they were before the epidemic. These results suggest that host recoveries are not caused by pathogen attenuation and may be due to shifts in host responses. Our findings provide insights into the mechanisms underlying disease transitions, which are increasingly important to understand in an era of emerging infectious diseases and unprecedented global pandemics.

Sources of Fungal Genetic Variation and Associating It with Phenotypic Diversity

The first eukaryotic genome to be sequenced was fungal, and there continue to be more sequenced genomes in the kingdom Fungi than in any other eukaryotic kingdom. Comparison of these genomes reveals many sources of genetic variation, from single nucleotide polymorphisms to horizontal gene transfer and on to changes in the arrangement and number of chromosomes, not to mention endofungal bacteria and viruses. Population genomics shows that all sources generate variation all the time and implicate natural selection as the force maintaining genome stability. Variation in wild populations is a rich resource for associating genetic variation with phenotypic variation, whether through quantitative trait locus mapping, genome-wide association studies, or reverse ecology. Subjects of studies associating genetic and phenotypic variation include model fungi, e.g., Saccharomyces and Neurospora, but pioneering studies have also been made with fungi pathogenic to plants, e.g., Pyricularia (= Magnaporthe), Zymoseptoria, and Fusarium, and to humans, e.g., Coccidioides, Cryptococcus, and Candida.

Amphibian species traits, evolutionary history and environment predict Batrachochytrium dendrobatidis infection patterns, but not extinction risk

The fungal pathogen Batrachochytrium dendrobatidis (B. dendrobatidis) has emerged as a major agent of amphibian extinction, requiring conservation intervention for many susceptible species. Identifying susceptible species is challenging, but many aspects of species biology are predicted to influence the evolution of host resistance, tolerance, or avoidance strategies towards disease. In turn, we may expect species exhibiting these distinct strategies to differ in their ability to survive epizootic disease outbreaks. Here, we test for phylogenetic and trait-based patterns of B. dendrobatidis infection risk and infection intensity among 302 amphibian species by compiling a global data set of B. dendrobatidis infection surveys across 95 sites. We then use best-fit models that associate traits, taxonomy and environment with B. dendrobatidis infection risk and intensity to predict host disease mitigation strategies (tolerance, resistance, avoidance) for 122 Neotropical amphibian species that experienced epizootic B. dendrobatidis outbreaks, and noted species persistence or extinction from these events. Aspects of amphibian species life history, habitat use and climatic niche were consistently linked to variation in B. dendrobatidis infection patterns across sites around the world. However, predicted B. dendrobatidis infection risk and intensity based on site environment and species traits did not reveal a consistent pattern between the predicted host disease mitigation strategy and extinction outcome. This suggests that either tolerant or resistant species may have no advantage in ameliorating disease during epizootic events, or that other factors drive the persistence of amphibian populations during chytridiomycosis outbreaks. These results suggest that using a trait-based approach may allow us to identify species with resistance or tolerance to endemic B. dendrobatidis infections, but that this approach may be insufficient to ultimately identify species at risk of extinction from epizootics.

Ploidy dynamics and evolvability in fungi

Rapid responses to acute stresses are essential for stress survival and are critical to the ability of fungal pathogens to adapt to new environments or hosts. The rapid emergence of drug resistance is used as a model for how fungi adapt and survive stress conditions that inhibit the growth of progenitor cells. Aneuploidy and loss of heterozygosity (LOH), which are large-scale genome shifts involving whole chromosomes or chromosome arms, occur at higher frequency than point mutations and have the potential to mediate stress survival. Furthermore, the stress of exposure to an antifungal drug can induce elevated levels of LOH and can promote the formation of aneuploids. This occurs via mitotic defects that first produce tetraploid progeny with extra spindles, followed by chromosome mis-segregation. Thus, drug exposure induces elevated levels of aneuploidy, which can alter the copy number of genes that improve survival in a given stress or drug. Selection then acts to increase the proportion of adaptive aneuploids in the population. Because aneuploidy is a common property of many pathogenic fungi, including those posing emerging threats to plants, animals and humans, we propose that aneuploid formation and LOH often accompanying it contribute to the rapid generation of diversity that can facilitate the emergence of fungal pathogens to new environmental niches and/or new hosts, as well as promote antifungal drug resistance that makes emerging fungal infections ever more difficult to contain.This article is part of the themed issue 'Tackling emerging fungal threats to animal health, food security and ecosystem resilience'.

Resistance, tolerance and environmental transmission dynamics determine host extinction risk in a load-dependent amphibian disease

While disease-induced extinction is generally considered rare, a number of recently emerging infectious diseases with load-dependent pathology have led to extinction in wildlife populations. Transmission is a critical factor affecting disease-induced extinction, but the relative importance of transmission compared to load-dependent host resistance and tolerance is currently unknown. Using a combination of models and experiments on an amphibian species suffering extirpations from the fungal pathogen Batrachochytrium dendrobatidis (Bd), we show that while transmission from an environmental Bd reservoir increased the ability of Bd to invade an amphibian population and the extinction risk of that population, Bd-induced extinction dynamics were far more sensitive to host resistance and tolerance than to Bd transmission. We demonstrate that this is a general result for load-dependent pathogens, where non-linear resistance and tolerance functions can interact such that small changes in these functions lead to drastic changes in extinction dynamics.

Batrachochytrium dendrobatidis and the Decline and Survival of the Relict Leopard Frog

Epizootic disease caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd) is a major driver of amphibian declines, yet many amphibians declined before the pathogen was described. The Relict Leopard Frog, Rana onca (=Lithobates onca), was nearly extinct, with the exception of populations within a few geothermal springs. Growth of Bd, however, is limited by high water temperature, and geothermal springs may have provided refuge during outbreaks of chytridiomycosis. We conducted field surveys and laboratory experiments to assess the susceptibility of R. onca to Bd. In the field, we found Bd at one of the two areas where remnant populations of R. onca still occur, but not in the other. In the laboratory, we infected juvenile frogs from these two areas with two hypervirulent Bd isolates associated with declines in other ranid species. In our experiments, these Bd isolates did not affect survivorship of R. onca and most infections (64%) were cleared by the end of the experiments. We propose that R. onca either has inherent resistance to Bd or has recently evolved such resistance. These results may be important for conservation efforts aimed at establishing new populations of R. onca across a landscape where Bd exists. Resistance, however, varies among life stages, and we also did not assess Bd from the local environment. We caution that the resistance we observed for young frogs under laboratory conditions may not translate to the situation for R. onca in the wild.

Modulation of Aneuploidy in Leishmania donovani during Adaptation to Different In Vitro and In Vivo Environments and Its Impact on Gene Expression

Aneuploidy is usually deleterious in multicellular organisms but appears to be tolerated and potentially beneficial in unicellular organisms, including pathogens. Leishmania, a major protozoan parasite, is emerging as a new model for aneuploidy, since in vitro-cultivated strains are highly aneuploid, with interstrain diversity and intrastrain mosaicism. The alternation of two life stages in different environments (extracellular promastigotes and intracellular amastigotes) offers a unique opportunity to study the impact of environment on aneuploidy and gene expression. We sequenced the whole genomes and transcriptomes of Leishmania donovani strains throughout their adaptation to in vivo conditions mimicking natural vertebrate and invertebrate host environments. The nucleotide sequences were almost unchanged within a strain, in contrast to highly variable aneuploidy. Although high in promastigotes in vitro, aneuploidy dropped significantly in hamster amastigotes, in a progressive and strain-specific manner, accompanied by the emergence of new polysomies. After a passage through a sand fly, smaller yet consistent karyotype changes were detected. Changes in chromosome copy numbers were correlated with the corresponding transcript levels, but additional aneuploidy-independent regulation of gene expression was observed. This affected stage-specific gene expression, downregulation of the entire chromosome 31, and upregulation of gene arrays on chromosomes 5 and 8. Aneuploidy changes in Leishmania are probably adaptive and exploited to modulate the dosage and expression of specific genes; they are well tolerated, but additional mechanisms may exist to regulate the transcript levels of other genes located on aneuploid chromosomes. Our model should allow studies of the impact of aneuploidy on molecular adaptations and cellular fitness.IMPORTANCE Aneuploidy is usually detrimental in multicellular organisms, but in several microorganisms, it can be tolerated and even beneficial. Leishmania-a protozoan parasite that kills more than 30,000 people each year-is emerging as a new model for aneuploidy studies, as unexpectedly high levels of aneuploidy are found in clinical isolates. Leishmania lacks classical regulation of transcription at initiation through promoters, so aneuploidy could represent a major adaptive strategy of this parasite to modulate gene dosage in response to stressful environments. For the first time, we document the dynamics of aneuploidy throughout the life cycle of the parasite, in vitro and in vivo We show its adaptive impact on transcription and its interaction with regulation. Besides offering a new model for aneuploidy studies, we show that further genomic studies should be done directly in clinical samples without parasite isolation and that adequate methods should be developed for this.

First in Vivo Batrachochytrium dendrobatidis Transcriptomes Reveal Mechanisms of Host Exploitation, Host-Specific Gene Expression, and Expressed Genotype Shifts

For generalist pathogens, host species represent distinct selective environments, providing unique challenges for resource acquisition and defense from host immunity, potentially resulting in host-dependent differences in pathogen fitness. Gene expression modulation should be advantageous, responding optimally to a given host and mitigating the costs of generalism. Batrachochytrium dendrobatidis (Bd), a fungal pathogen of amphibians, shows variability in pathogenicity among isolates, and within-strain virulence changes rapidly during serial passages through artificial culture. For the first time, we characterize the transcriptomic profile of Bd in vivo, using laser-capture microdissection. Comparison of Bd transcriptomes (strain JEL423) in culture and in two hosts (Atelopus zeteki and Hylomantis lemur), reveals >2000 differentially expressed genes that likely include key Bd defense and host exploitation mechanisms. Variation in Bd transcriptomes from different amphibian hosts demonstrates shifts in pathogen resource allocation. Furthermore, expressed genotype variant frequencies of Bd populations differ between culture and amphibian skin, and among host species, revealing potential mechanisms underlying rapid changes in virulence and the possibility that amphibian community composition shapes Bd evolutionary trajectories. Our results provide new insights into how changes in gene expression and infecting population genotypes can be key to the success of a generalist fungal pathogen.

Copy number variation contributes to cryptic genetic variation in outbreak lineages of Cryptococcus gattii from the North American Pacific Northwest

Background

Copy number variants (CNVs) are a class of structural variants (SVs) and are defined as fragments of DNA that are present at variable copy number in comparison with a reference genome. Recent advances in bioinformatics methodologies and sequencing technologies have enabled the high-resolution quantification of genome-wide CNVs. In pathogenic fungi SVs have been shown to alter gene expression, influence host specificity, and drive fungicide resistance, but little attention has focused specifically on CNVs. Using publicly available sequencing data, we identified 90 isolates across 212 Cryptococcus gattii genomes that belong to the VGII subgroups responsible for the recent deadly outbreaks in the North American Pacific Northwest. We generated CNV profiles for each sample to investigate the prevalence and function of CNV in C. gattii.

Results

We identified eight genetic clusters among publicly available Illumina whole genome sequence data from 212 C. gattii isolates through population structure analysis. Three clusters represent the VGIIa, VGIIb, and VGIIc subgroups from the North American Pacific Northwest. CNV was bioinformatically predicted and affected ~300–400 Kilobases (Kb) of the C. gattii VGII subgroup genomes. Sixty-seven loci, encompassing 58 genes, showed highly divergent patterns of copy number variation between VGII subgroups. Analysis of PFam domains within divergent CN variable genes revealed enrichment of protein domains associated with transport, cell wall organization and external encapsulating structure.

Conclusions

CNVs may contribute to pathological and phenotypic differences observed between the C. gattii VGIIa, VGIIb, and VGIIc subpopulations. Genes overlapping with population differentiated CNVs were enriched for several virulence related functional terms. These results uncover novel candidate genes to examine the genetic and functional underpinnings of C. gattii pathogenicity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3044-0) contains supplementary material, which is available to authorized users.

Projecting the Global Distribution of the Emerging Amphibian Fungal Pathogen, Batrachochytrium dendrobatidis, Based on IPCC Climate Futures

Projected changes in climate conditions are emerging as significant risk factors to numerous species, affecting habitat conditions and community interactions. Projections suggest species range shifts in response to climate change modifying environmental suitability and is supported by observational evidence. Both pathogens and their hosts can shift ranges with climate change. We consider how climate change may influence the distribution of the emerging infectious amphibian chytrid fungus, Batrachochytrium dendrobatidis (Bd), a pathogen associated with worldwide amphibian population losses. Using an expanded global Bd database and a novel modeling approach, we examined a broad set of climate metrics to model the Bd-climate niche globally and regionally, then project how climate change may influence Bd distributions. Previous research showed that Bd distribution is dependent on climatic variables, in particular temperature. We trained a machine-learning model (random forest) with the most comprehensive global compilation of Bd sampling records (~5,000 site-level records, mid-2014 summary), including 13 climatic variables. We projected future Bd environmental suitability under IPCC scenarios. The learning model was trained with combined worldwide data (non-region specific) and also separately per region (region-specific). One goal of our study was to estimate of how Bd spatial risks may change under climate change based on the best available data. Our models supported differences in Bd-climate relationships among geographic regions. We projected that Bd ranges will shift into higher latitudes and altitudes due to increased environmental suitability in those regions under predicted climate change. Specifically, our model showed a broad expansion of areas environmentally suitable for establishment of Bd on amphibian hosts in the temperate zones of the Northern Hemisphere. Our projections are useful for the development of monitoring designs in these areas, especially for sensitive species and those vulnerable to multiple threats.

Insights From Genomics Into Spatial and Temporal Variation in Batrachochytrium dendrobatidis

Advances in genetics and genomics have provided new tools for the study of emerging infectious diseases. Researchers can now move quickly from simple hypotheses to complex explanations for pathogen origin, spread, and mechanisms of virulence. Here we focus on the application of genomics to understanding the biology of the fungal pathogen Batrachochytrium dendrobatidis (Bd), a novel and deadly pathogen of amphibians. We provide a brief history of the system, then focus on key insights into Bd variation garnered from genomics approaches, and finally, highlight new frontiers for future discoveries. Genomic tools have revealed unexpected complexity and variation in the Bd system suggesting that the history and biology of emerging pathogens may not be as simple as they initially seem.

Ploidy plasticity: a rapid and reversible strategy for adaptation to stress

Organisms must be able to grow in a broad range of conditions found in their normal growth environment and for a species to survive, at least some cells in a population must adapt rapidly to extreme stress conditions that kill the majority of cells.Candida albicans, the most prevalent fungal pathogen of humans resides as a commensal in a broad range of niches within the human host. Growth conditions in these niches are highly variable and stresses such exposure to antifungal drugs can inhibit population growth abruptly. One of the mechanisms C. albicans uses to adapt rapidly to severe stresses is aneuploidy-a change in the total number of chromosomes such that one or more chromosomes are present in excess or are missing. Aneuploidy is quite common in wild isolates of fungi and other eukaryotic microbes. Aneuploidy can be achieved by chromosome nondisjunction during a simple mitosis, and in stress conditions it begins to appear after two mitotic divisions via a tetraploid intermediate. Aneuploidy usually resolves to euploidy (a balanced number of chromosomes), but not necessarily to diploidy. Aneuploidy of a specific chromosome can confer new phenotypes by virtue of the copy number of specific genes on that chromosome relative to the copies of other genes. Thus, it is not aneuploidy per se, but the relative copy number of specific genes that confers many tested aneuploidy-associated phenotypes. Aneuploidy almost always carries a fitness cost, as cells express most proteins encoded by genes on the aneuploid chromosome in proportion to the number of DNA copies of the gene. This is thought to be due to imbalances in the stoichiometry of different components of large complexes. Despite this, fitness is a relative function-and if stress is severe and population growth has slowed considerably, then even small growth advantages of some aneuploidies can provide a selective advantage. Thus, aneuploidy appears to provide a transient solution to severe and sudden stress conditions, and may promote the appearance of more stable solutions as well. Importantly, in many clinical and environmental isolates of different fungal species aneuploidy does not appear to have a high fitness cost, and is well-tolerated. Thus, rapid changes in ploidy may provide the opportunity for rapid adaptation to stress conditions in the environment, host niches or in response to antifungal drugs.