Substrate-specific gene expression in Batrachochytrium dendrobatidis, the chytrid pathogen of amphibians

Small Interfering RNA Mediated Messenger RNA Knockdown in the Amphibian Pathogen Batrachochytrium dendrobatidis

RNA interference (RNAi) has not been tested in the pandemic amphibian pathogen, Batrachochytrium dendrobatidis, but developing this technology could be useful to elucidate virulence mechanisms, identify therapeutic targets, and may present a novel antifungal treatment option for chytridiomycosis. To manipulate and decipher gene function, rationally designed small interfering RNA (siRNA) can initiate the destruction of homologous messenger RNA (mRNA), resulting in the "knockdown" of target gene expression. Here, we investigate whether siRNA can be used to manipulate gene expression in B. dendrobatidis via RNAi using differing siRNA strategies to target genes involved in glutathione and ornithine synthesis. To determine the extent and duration of mRNA knockdown, target mRNA levels were monitored for 24-48 h after delivery of siRNA targeting glutamate-cysteine ligase, with a maximum of ~56% reduction in target transcripts occurring at 36 h. A second siRNA design targeting glutamate-cysteine ligase also resulted in ~53% knockdown at this time point. siRNA directed toward a different gene target, ornithine decarboxylase, achieved 17% reduction in target transcripts. Although no phenotypic effects were observed, these results suggest that RNAi is possible in B. dendrobatidis, and that gene expression can be manipulated in this pathogen. We outline ideas for further optimization steps to increase knockdown efficiency to better harness RNAi techniques for control of B. 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.

Conserved Evolution of MHC Supertypes among Japanese Frogs Suggests Selection for Bd Resistance

The chytrid fungus Batrachochytrium dendrobatidis (Bd) is a major threat to amphibians, yet there are no reports of major disease impacts in East Asian frogs. Genetic variation of the major histocompatibility complex (MHC) has been associated with resistance to Bd in frogs from East Asia and worldwide. Using transcriptomic data collated from 11 Japanese frog species (one individual per species), we isolated MHC class I and IIb sequences and validated using molecular cloning. We then compared MHC from Japanese frogs and other species worldwide, with varying Bd susceptibility. Supertyping analysis, which groups MHC alleles based on physicochemical properties of peptide binding sites, identified that all examined East Asian frogs contained at least one MHC-IIb allele belonging to supertype ST-1. This indicates that, despite the large divergence times between some Japanese frogs (up to 145 million years), particular functional properties in the peptide binding sites of MHC-II are conserved among East Asian frogs. Furthermore, preliminary analysis using NetMHCIIpan-4.0, which predicts potential Bd-peptide binding ability, suggests that MHC-IIb ST-1 and ST-2 have higher overall peptide binding ability than other supertypes, irrespective of whether the peptides are derived from Bd, other fungi, or bacteria. Our findings suggest that MHC-IIb among East Asian frogs may have co-evolved under the same selective pressure. Given that Bd originated in this region, it may be a major driver of MHC evolution in East Asian frogs.

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

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

Immune priming prior to pathogen exposure sheds light on the relationship between host, microbiome and pathogen in disease

Dynamic interactions between host, pathogen and host-associated microbiome dictate infection outcomes. Pathogens including Batrachochytrium dendrobatidis (Bd) threaten global biodiversity, but conservation efforts are hindered by limited understanding of amphibian host, Bd and microbiome interactions. We conducted a vaccination and infection experiment using Eastern hellbender salamanders (Cryptobranchus alleganiensis alleganiensis) challenged with Bd to observe infection, skin microbial communities and gene expression of host skin, pathogen and microbiome throughout the experiment. Most animals survived high Bd loads regardless of their vaccination status and vaccination did not affect pathogen load, but host gene expression differed based on vaccination. Oral vaccination (exposure to killed Bd) stimulated immune gene upregulation while topically and sham-vaccinated animals did not significantly upregulate immune genes. In early infection, topically vaccinated animals upregulated immune genes but orally and sham-vaccinated animals downregulated immune genes. Bd increased pathogenicity-associated gene expression in late infection when Bd loads were highest. The microbiome was altered by Bd, but there was no correlation between anti-Bd microbe abundance or richness and pathogen burden. Our observations suggest that hellbenders initially generate a vigorous immune response to Bd, which is ineffective at controlling disease and is subsequently modulated. Interactions with antifungal skin microbiota did not influence disease progression.

Panzootic chytrid fungus exploits diverse amphibian host environments through plastic infection strategies

While some pathogens are limited to single species, others can colonize many hosts, likely contributing to the emergence of novel disease outbreaks. Despite this biodiversity threat, traits associated with host niche expansions are not well understood in multi-host pathogens. Here, we aimed to uncover functional machinery driving multi-host invasion by focusing on Batrachochytrium dendrobatidis (Bd), a pathogen that infects the skin of hundreds of amphibians worldwide. We performed a meta-analysis of Bd gene expression using data from published infection experiments and newly generated profiles. We analyzed Bd transcriptomic landscapes across the skin of 14 host species, reconstructed Bd isolates phylogenetic relationships, and inferred the origin and evolutionary history of differentially expressed genes under a phylogenetic framework comprising other 12 zoosporic fungi. Bd displayed plastic infection strategies when challenged by hosts with different disease susceptibility. Our analyses identified sets of differentially expressed genes under host environments with similar infection outcome. We stressed nutritional immunity and gene silencing as important processes required to overcome challenging skin environments in less susceptible hosts. Overall, Bd genes expressed during amphibian skin exploitation have arisen mainly via gene duplications with great family expansions, increasing the gene copy events previously described for this fungal species. Finally, we provide a comprehensive gene dataset that can be used to further examine eco-evolutionary hypotheses for this host-pathogen system. Our study supports the idea that host environments exert contrasting selective pressures, such that gene expression plasticity could be one of the evolutionary keys leading to the success of multi-host pathogens.

In Silico Predictions of Ecological Plasticity Mediated by Protein Family Expansions in Early-Diverging Fungi

Early-diverging fungi (EDF) are ubiquitous and versatile. Their diversity is reflected in their genome sizes and complexity. For instance, multiple protein families have been reported to expand or disappear either in particular genomes or even whole lineages. The most commonly mentioned are CAZymes (carbohydrate-active enzymes), peptidases and transporters that serve multiple biological roles connected to, e.g., metabolism and nutrients intake. In order to study the link between ecology and its genomic underpinnings in a more comprehensive manner, we carried out a systematic in silico survey of protein family expansions and losses among EDF with diverse lifestyles. We found that 86 protein families are represented differently according to EDF ecological features (assessed by median count differences). Among these there are 19 families of proteases, 43 CAZymes and 24 transporters. Some of these protein families have been recognized before as serine and metallopeptidases, cellulases and other nutrition-related enzymes. Other clearly pronounced differences refer to cell wall remodelling and glycosylation. We hypothesize that these protein families altogether define the preliminary fungal adaptasome. However, our findings need experimental validation. Many of the protein families have never been characterized in fungi and are discussed in the light of fungal ecology for the first time.

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.

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.

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.

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.

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.

The Synchytrium endobioticum AvrSen1 Triggers a Hypersensitive Response in Sen1 Potatoes While Natural Variants Evade Detection

Synchytrium endobioticum is an obligate biotrophic fungus of division Chytridiomycota. It causes potato wart disease, has a worldwide quarantine status and is included on the Health and Human Services and United States Department of Agriculture Select Agent list. S. endobioticum isolates are grouped in pathotypes based on their ability to evade host resistance in a set of differential potato varieties. Thus far, 39 pathotypes are reported. A single dominant gene (Sen1) governs pathotype 1 (D1) resistance and we anticipated that the underlying molecular model would involve a pathogen effector (AvrSen1) that is recognized by the host. The S. endobioticum-specific secretome of 14 isolates representing six different pathotypes was screened for effectors specifically present in pathotype 1 (D1) isolates but absent in others. We identified a single AvrSen1 candidate. Expression of this candidate in potato Sen1 plants showed a specific hypersensitive response (HR), which cosegregated with the Sen1 resistance in potato populations. No HR was obtained with truncated genes found in pathotypes that evaded recognition by Sen1. These findings established that our candidate gene was indeed Avrsen1. The S. endobioticum AvrSen1 is a single-copy gene and encodes a 376-amino-acid protein without predicted function or functional domains, and is the first effector gene identified in Chytridiomycota, an extremely diverse yet underrepresented basal lineage of fungi.

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.

Infection mechanisms and putative effector repertoire of the mosquito pathogenic oomycete Pythium guiyangense uncovered by genomic analysis

Pythium guiyangense, an oomycete from a genus of mostly plant pathogens, is an effective biological control agent that has wide potential to manage diverse mosquitoes. However, its mosquito-killing mechanisms are almost unknown. In this study, we observed that P. guiyangense could utilize cuticle penetration and ingestion of mycelia into the digestive system to infect mosquito larvae. To explore pathogenic mechanisms, a high-quality genome sequence with 239 contigs and an N50 contig length of 1,009 kb was generated. The genome assembly is approximately 110 Mb, which is almost twice the size of other sequenced Pythium genomes. Further genome analysis suggests that P. guiyangense may arise from a hybridization of two related but distinct parental species. Phylogenetic analysis demonstrated that P. guiyangense likely evolved from common ancestors shared with plant pathogens. Comparative genome analysis coupled with transcriptome sequencing data suggested that P. guiyangense may employ multiple virulence mechanisms to infect mosquitoes, including secreted proteases and kazal-type protease inhibitors. It also shares intracellular Crinkler (CRN) effectors used by plant pathogenic oomycetes to facilitate the colonization of plant hosts. Our experimental evidence demonstrates that CRN effectors of P. guiyangense can be toxic to insect cells. The infection mechanisms and putative virulence effectors of P. guiyangense uncovered by this study provide the basis to develop improved mosquito control strategies. These data also provide useful knowledge on host adaptation and evolution of the entomopathogenic lifestyle within the oomycete lineage. A deeper understanding of the biology of P. guiyangense effectors might also be useful for management of other important agricultural pests.

Utilization of biocontrol agents has emerged as a promising mosquito control strategy, and Pythium guiyangense has wide potential to manage diverse mosquitoes with high efficiency. However, the molecular mechanisms underlying pathological processes remain almost unknown. We observed that P. guiyangense invades mosquito larvae through cuticle penetration and through ingestion of mycelia via the digestive system, jointly accelerating mosquito larvae mortality. We also present a high-quality genome assembly of P. guiyangense that contains two distinct genome complements, which likely resulted from a hybridization of two parental species. Our analyses revealed expansions of kinases, proteases, kazal-type protease inhibitors, and elicitins that may be important for adaptation of P. guiyangense to a mosquito-pathogenic lifestyle. Moreover, our experimental evidence demonstrated that some Crinkler effectors of P. guiyangense can be toxic to insect cells. Our findings suggest new insights into oomycete evolution and host adaptation by animal pathogenic oomycetes. Our new genome resource will enable better understanding of infection mechanisms, with the potential to improve the biological control of mosquitoes and other agriculturally important pests.

Phenotypic Divergence along Geographic Gradients Reveals Potential for Rapid Adaptation of the White-Nose Syndrome Pathogen, Pseudogymnoascus destructans, in North America

White-nose syndrome (WNS) is an ongoing epizootic affecting multiple species of North American bats, caused by epidermal infections of the psychrophilic filamentous fungus Pseudogymnoascus destructans Since its introduction from Europe, WNS has spread rapidly across eastern North America and resulted in high mortality rates in bats. At present, the mechanisms behind its spread and the extent of its adaptation to different geographic and ecological niches remain unknown. The objective of this study was to examine the geographic patterns of phenotypic variation and the potential evidence for adaptation among strains representing broad geographic locations in eastern North America. The morphological features of these strains were evaluated on artificial medium, and the viability of asexual arthroconidia of representative strains was investigated after storage at high (23°C), moderate (14°C), and low (4°C) temperatures at different lengths of time. Our analyses identified evidence for a geographic pattern of colony morphology changes among the clonal descendants of the fungus, with trait values correlated with increased distance from the epicenter of WNS. Our genomic comparisons of three representative isolates revealed novel genetic polymorphisms and suggested potential candidate mutations that might be related to some of the phenotypic changes. These results show that even though this pathogen arrived in North America only recently and reproduces asexually, there has been substantial evolution and phenotypic diversification during its rapid clonal expansion.IMPORTANCE The causal agent of white-nose syndrome in bats is Pseudogymnoascus destructans, a filamentous fungus recently introduced from its native range in Europe. Infections caused by P. destructans have progressed across the eastern parts of Canada and the United States over the last 10 years. It is not clear how the disease is spread, as the pathogen is unable to grow above 23°C and ambient temperature can act as a barrier when hosts disperse. Here, we explore the patterns of phenotypic diversity and the germination of the fungal asexual spores, arthroconidia, from strains across a sizeable area of the epizootic range. Our analyses revealed evidence of adaptation along geographic gradients during its expansion. The results have implications for understanding the diversification of P. destructans and the limits of WNS spread in North America. Given the rapidly expanding distribution of WNS, a detailed understanding of the genetic bases for phenotypic variations in growth, reproduction, and dispersal of P. destructans is urgently needed to help control this disease.

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

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

Development and worldwide use of non-lethal, and minimal population-level impact, protocols for the isolation of amphibian chytrid fungi

Parasitic chytrid fungi have emerged as a significant threat to amphibian species worldwide, necessitating the development of techniques to isolate these pathogens into culture for research purposes. However, early methods of isolating chytrids from their hosts relied on killing amphibians. We modified a pre-existing protocol for isolating chytrids from infected animals to use toe clips and biopsies from toe webbing rather than euthanizing hosts, and distributed the protocol to researchers as part of the BiodivERsA project RACE; here called the RML protocol. In tandem, we developed a lethal procedure for isolating chytrids from tadpole mouthparts. Reviewing a database of use a decade after their inception, we find that these methods have been applied across 5 continents, 23 countries and in 62 amphibian species. Isolation of chytrids by the non-lethal RML protocol occured in 18% of attempts with 207 fungal isolates and three species of chytrid being recovered. Isolation of chytrids from tadpoles occured in 43% of attempts with 334 fungal isolates of one species (Batrachochytrium dendrobatidis) being recovered. Together, these methods have resulted in a significant reduction and refinement of our use of threatened amphibian species and have improved our ability to work with this group of emerging pathogens.

What We Do Not Know about Fungal Cell Adhesion Molecules

There has been extensive research on structure and function of fungal cell adhesion molecules, but the most of the work has been about adhesins in Candida albicans and Saccharomyces cerevisiae. These yeasts are members of a single ascomycete order, and adhesion molecules from the six other fungal phyla are only sparsely described in the literature. In these other phyla, most of the research is at the cellular level, rather than at the molecular level, so there has been little characterization of the adhesion molecules themselves. A catalog of known adhesins shows some common features: high Ser/Thr content, tandem repeats, N- and O-glycosylations, GPI anchors, dibasic sequence motifs, and potential amyloid-forming sequences. However, none of these features is universal. Known ligands include proteins and glycans on homologous cells and host cells. Existing and novel tools can exploit the availability of genome sequences to identify and characterize new fungal adhesins. These include bioinformatics tools and well-established yeast surface display models, which could be coupled with an adhesion substrate array. Thus, new knowledge could be exploited to answer key questions in fungal ecology, animal and plant pathogenesis, and roles of biofilms in infection and biomass turnover.

Lung epithelial cells have virus-specific and shared gene expression responses to infection by diverse respiratory viruses

The severity of respiratory viral infections is partially determined by the cellular response mounted by infected lung epithelial cells. Disease prevention and treatment is dependent on our understanding of the shared and unique responses elicited by diverse viruses, yet few studies compare host responses to viruses from different families while controlling other experimental parameters. Murine models are commonly used to study the pathogenesis of respiratory viral infections, and in vitro studies using murine cells provide mechanistic insight into the pathogenesis observed in vivo. We used microarray analysis to compare changes in gene expression of murine lung epithelial cells infected individually by three respiratory viruses causing mild (rhinovirus, RV1B), moderate (coronavirus, MHV-1), and severe (influenza A virus, PR8) disease in mice. RV1B infection caused numerous gene expression changes, but the differential effect peaked at 12 hours post-infection. PR8 altered an intermediate number of genes whose expression continued to change through 24 hours. MHV-1 had comparatively few effects on host gene expression. The viruses elicited highly overlapping responses in antiviral genes, though MHV-1 induced a lower type I interferon response than the other two viruses. Signature genes were identified for each virus and included host defense genes for PR8, tissue remodeling genes for RV1B, and transcription factors for MHV-1. Our comparative approach identified universal and specific transcriptional signatures of virus infection that can be used to distinguish shared and virus-specific mechanisms of pathogenesis in the respiratory tract.

Characterization of Batrachochytrium dendrobatidis Inhibiting Bacteria from Amphibian Populations in Costa Rica

Global amphibian declines and extinction events are occurring at an unprecedented rate. While several factors are responsible for declines and extinction, the fungal pathogen Batrachochytrium dendrobatidis (Bd) has been cited as a major constituent in these events. While the effects of this chytrid fungus have been shown to cause broad scale population declines and extinctions, certain individuals and relict populations have shown resistance. This resistance has been attributed in part to the cutaneous bacterial microbiome. Here, we present the first study characterizing anti-Bd bacterial isolates from amphibian populations in Costa Rica, including the characterization of two strains of Serratia marcescens presenting strong anti-Bd activity. Transcriptome sequencing was utilized for delineation of shifts in gene expression of the two previously uncharacterized strains of S. marcescens grown in three different treatments comprising Bd, heat-killed Bd, and a no Bd control. These results revealed up- and down-regulation of key genes associated with different metabolic and regulatory pathways. This information will be valuable in continued efforts to develop a bacterial-based approach for amphibian protection as well as providing direction for continued mechanistic inquiries of the bacterial anti-Bd response.

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.

Glycoside hydrolases family 20 (GH20) represent putative virulence factors that are shared by animal pathogenic oomycetes, but are absent in phytopathogens

BACKGROUND

Although interest in animal pathogenic oomycetes is increasing, the molecular basis mediating oomycete-animal relationships remains virtually unknown. Crinkler (CRN) genes, which have been traditionally associated with the cytotoxic activity displayed by plant pathogenic oomycetes, were recently detected in transcriptome sequences from the entomopathogenic oomycete Lagenidium giganteum, suggesting that these genes may represent virulence factors conserved in both animal and plant pathogenic oomycetes. In order to further characterize the L. giganteum pathogenome, an on-going genomic survey was mined to reveal novel putative virulence factors, including canonical oomycete effectors Crinkler 13 (CRN13) orthologs. These novel sequences provided a basis to initiate gene expression analyses and determine if the proposed L. giganteum virulence factors are differentially expressed in the presence of mosquito larvae (Aedes aegypti).

RESULTS

Sequence analyses revealed that L. giganteum express CRN13 transcripts. The predicted proteins, like other L. giganteum CRNs, contained a conserved LYLA motif at the N terminal, but did not display signal peptides. In contrast, other potential virulence factors, such as Glycoside Hydrolases family 20 (hexosaminidase) and 37 (trehalase) proteins (GH20 and GH37), contained identifiable signal peptides. Genome mining demonstrated that GH20 genes are absent from phytopathogenic oomycete genomes, and that the L. giganteum GH20 sequence is the only reported peronosporalean GH20 gene. All other oomycete GH20 homologs were retrieved from animal pathogenic, saprolegnialean genomes. Furthermore, phylogenetic analyses demonstrated that saprolegnialean and peronosporalean GH20 protein sequences clustered in unrelated clades. The saprolegnialean GH20 sequences appeared as a strongly supported, monophyletic group nested within an arthropod-specific clade, suggesting that this gene was acquired via a lateral gene transfer event from an insect or crustacean genome. In contrast, the L. giganteum GH20 protein sequence appeared as a sister taxon to a plant-specific clade that included exochitinases with demonstrated insecticidal activities. Finally, gene expression analyses demonstrated that the L. giganteum GH20 gene expression level is significantly modulated in the presence of mosquito larvae. In agreement with the protein secretion predictions, CRN transcripts did not show any differential expression.

CONCLUSIONS

These results identified GH20 enzymes, and not CRNs, as potential pathogenicity factors shared by all animal pathogenic oomycetes.

Contribution of Multiple Inter-Kingdom Horizontal Gene Transfers to Evolution and Adaptation of Amphibian-Killing Chytrid, Batrachochytrium dendrobatidis

Amphibian populations are experiencing catastrophic declines driven by the fungal pathogen Batrachochytrium dendrobatidis (Bd). Although horizontal gene transfer (HGT) facilitates the evolution and adaptation in many fungi by conferring novel function genes to the recipient fungi, inter-kingdom HGT in Bd remains largely unexplored. In this study, our investigation detects 19 bacterial genes transferred to Bd, including metallo-beta-lactamase and arsenate reductase that play important roles in the resistance to antibiotics and arsenates. Moreover, three probable HGT gene families in Bd are from plants and one gene family coding the ankyrin repeat-containing protein appears to come from oomycetes. The observed multi-copy gene families associated with HGT are probably due to the independent transfer events or gene duplications. Five HGT genes with extracellular locations may relate to infection, and some other genes may participate in a variety of metabolic pathways, and in doing so add important metabolic traits to the recipient. The evolutionary analysis indicates that all the transferred genes evolved under purifying selection, suggesting that their functions in Bd are similar to those of the donors. Collectively, our results indicate that HGT from diverse donors may be an important evolutionary driver of Bd, and improve its adaptations for infecting and colonizing host amphibians.

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.

Transposons to toxins: the provenance, architecture and diversification of a widespread class of eukaryotic effectors

Enzymatic effectors targeting nucleic acids, proteins and other cellular components are the mainstay of conflicts across life forms. Using comparative genomics we identify a large class of eukaryotic proteins, which include effectors from oomycetes, fungi and other parasites. The majority of these proteins have a characteristic domain architecture with one of several N-terminal 'Header' domains, which are predicted to play a role in trafficking of these effectors, including a novel version of the Ubiquitin fold. The Headers are followed by one or more diverse C-terminal domains, such as restriction endonuclease (REase), protein kinase, HNH endonuclease, LK-nuclease (a RNase) and multiple distinct peptidase domains, which are predicted to carry their toxicity determinants. The most common types of these proteins appear to have originated from prokaryotic transposases (e.g. TN7 and Mu) and combine a CDC6/ORC1-STAND clade NTPase domain with a C-terminal REase domain. Other than the so-called Crinkler effectors of oomycetes and fungi, these effectors are encoded by other eukaryotic parasites such as trypanosomatids (the RHS proteins) and the rhizarian Plasmodiophora, and symbionts like Capsaspora Remarkably, we also find these proteins in free-living eukaryotes, including several viridiplantae, fungi, amoebozoans and animals. These versions might either still be transposons or function in other poorly understood eukaryote-specific inter-organismal and inter-genomic conflicts. These include the Medea1 selfish element of Tribolium that spreads via post-zygotic killing. We present a unified mechanism for the recombination-dependent diversification and action of this widespread class of molecular weaponry deployed across diverse conflicts ranging from parasitic to free-living forms.

CRN13 candidate effectors from plant and animal eukaryotic pathogens are DNA-binding proteins which trigger host DNA damage response

To successfully colonize their host, pathogens produce effectors that can interfere with host cellular processes. Here we investigated the function of CRN13 candidate effectors produced by plant pathogenic oomycetes and detected in the genome of the amphibian pathogenic chytrid fungus Batrachochytrium dendrobatidis (BdCRN13). When expressed in Nicotiana, AeCRN13, from the legume root pathogen Aphanomyces euteiches, increases the susceptibility of the leaves to the oomycete Phytophthora capsici. When transiently expressed in amphibians or plant cells, AeCRN13 and BdCRN13 localize to the cell nuclei, triggering aberrant cell development and eventually causing cell death. Using Förster resonance energy transfer experiments in plant cells, we showed that both CRN13s interact with nuclear DNA and trigger plant DNA damage response (DDR). Mutating key amino acid residues in a predicted HNH-like endonuclease motif abolished the interaction of AeCRN13 with DNA, the induction of DDR and the enhancement of Nicotiana susceptibility to P. capsici. Finally, H2AX phosphorylation, a marker of DNA damage, and enhanced expression of genes involved in the DDR were observed in A. euteiches-infected Medicago truncatula roots. These results show that CRN13 from plant and animal eukaryotic pathogens promotes host susceptibility by targeting nuclear DNA and inducing DDR.

An Oomycete CRN Effector Reprograms Expression of Plant HSP Genes by Targeting their Promoters

Oomycete pathogens produce a large number of CRN effectors to manipulate plant immune responses and promote infection. However, their functional mechanisms are largely unknown. Here, we identified a Phytophthora sojae CRN effector PsCRN108 which contains a putative DNA-binding helix-hairpin-helix (HhH) motif and acts in the plant cell nucleus. Silencing of the PsCRN108 gene reduced P. sojae virulence to soybean, while expression of the gene in Nicotiana benthamiana and Arabidopsis thaliana enhanced plant susceptibility to P. capsici. Moreover, PsCRN108 could inhibit expression of HSP genes in A. thaliana, N. benthamiana and soybean. Both the HhH motif and nuclear localization signal of this effector were required for its contribution to virulence and its suppression of HSP gene expression. Furthermore, we found that PsCRN108 targeted HSP promoters in an HSE- and HhH motif-dependent manner. PsCRN108 could inhibit the association of the HSE with the plant heat shock transcription factor AtHsfA1a, which initializes HSP gene expression in response to stress. Therefore, our data support a role for PsCRN108 as a nucleomodulin in down-regulating the expression of plant defense-related genes by directly targeting specific plant promoters.

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

Emerging infectious diseasespose a significant threat to global health, but predicting disease outcomes for particular species can be complicated when pathogen virulence varies across space, time, or hosts. The pathogenic chytrid fungus Batrachochytrium dendrobatidis (Bd) has caused worldwide declines in frog populations. Not only do Bd isolates from wild populations vary in virulence, but virulence shifts can occur over short timescales when Bd is maintained in the laboratory. We leveraged changes in Bd virulence over multiple generations of passage to better understand mechanisms of pathogen virulence. We conducted whole-genome resequencing of two samples of the same Bd isolate, differing only in passage history, to identify genomic processes associated with virulence attenuation. The isolate with shorter passage history (and greater virulence) had greater chromosome copy numbers than the isolate maintained in culture for longer, suggesting that virulence attenuation may be associated with loss of chromosome copies. Our results suggest that genomic processes proposed as mechanisms for rapid evolution in Bd are correlated with virulence attenuation in laboratory culture within a single lineage of Bd. Moreover, these genomic processes can occur over extremely short timescales. On a practical level, our results underscore the importance of immediately cryo-archiving new Bd isolates and using fresh isolates, rather than samples cultured in the laboratory for long periods, for laboratory infection experiments. Finally, when attempting to predict disease outcomes for this ecologically important pathogen, it is critical to consider existing variation in virulence among isolates and the potential for shifts in virulence over short timescales.

Quantitative Proteomics of an Amphibian Pathogen, Batrachochytrium dendrobatidis, following Exposure to Thyroid Hormone

Batrachochytrium dendrobatidis (Bd), a chytrid fungus, has increasingly been implicated as a major factor in the worldwide decline of amphibian populations. The fungus causes chytridiomycosis in susceptible species leading to massive die-offs of adult amphibians. Although Bd infects the keratinized mouthparts of tadpoles and negatively affects foraging behavior, these infections are non-lethal. An important morphogen controlling amphibian metamorphosis is thyroid hormone (T3). Tadpoles may be infected with Bd and the fungus may be exposed to T3 during metamorphosis. We hypothesize that exposure of Bd to T3 may induce the expression of factors associated with host colonization and pathogenicity. We utilized a proteomics approach to better understand the dynamics of the Bd-T3 interaction. Using liquid chromatography-mass spectrometry (LC-MS), we generated a data set of a large number of cytoplasmic and membrane proteins following exposure of Bd to T3. From these data, we identified a total of 263 proteins whose expression was significantly changed following T3 exposure. We provide evidence for expression of an array of proteins that may play key roles in both genomic and non-genomic actions of T3 in Bd. Additionally, our proteomics study shows an increase in several proteins including proteases and a class of uncommon crinkler and crinkler-like effector proteins suggesting their importance in Bd pathogenicity as well as those involved in metabolism and energy transfer, protein fate, transport and stress responses. This approach provides insights into the mechanistic basis of the Bd-amphibian interaction following T3 exposure.

Correlates of virulence in a frog-killing fungal pathogen: evidence from a California amphibian decline

The fungal pathogen Batrachochytrium dendrobatidis (Bd) has caused declines and extinctions in amphibians worldwide, and there is increasing evidence that some strains of this pathogen are more virulent than others. While a number of putative virulence factors have been identified, few studies link these factors to specific epizootic events. We documented a dramatic decline in juvenile frogs in a Bd-infected population of Cascades frogs (Rana cascadae) in the mountains of northern California and used a laboratory experiment to show that Bd isolated in the midst of this decline induced higher mortality than Bd isolated from a more stable population of the same species of frog. This highly virulent Bd isolate was more toxic to immune cells and attained higher density in liquid culture than comparable isolates. Genomic analyses revealed that this isolate is nested within the global panzootic lineage and exhibited unusual genomic patterns, including increased copy numbers of many chromosomal segments. This study integrates data from multiple sources to suggest specific phenotypic and genomic characteristics of the pathogen that may be linked to disease-related declines.

Chromosomal copy number variation, selection and uneven rates of recombination reveal cryptic genome diversity linked to pathogenicity

Pathogenic fungi constitute a growing threat to both plant and animal species on a global scale. Despite a clonal mode of reproduction dominating the population genetic structure of many fungi, putatively asexual species are known to adapt rapidly when confronted by efforts to control their growth and transmission. However, the mechanisms by which adaptive diversity is generated across a clonal background are often poorly understood. We sequenced a global panel of the emergent amphibian pathogen, Batrachochytrium dendrobatidis (Bd), to high depth and characterized rapidly changing features of its genome that we believe hold the key to the worldwide success of this organism. Our analyses show three processes that contribute to the generation of de novo diversity. Firstly, we show that the majority of wild isolates manifest chromosomal copy number variation that changes over short timescales. Secondly, we show that cryptic recombination occurs within all lineages of Bd, leading to large regions of the genome being in linkage equilibrium, and is preferentially associated with classes of genes of known importance for virulence in other pathosystems. Finally, we show that these classes of genes are under directional selection, and that this has predominantly targeted the Global Panzootic Lineage (BdGPL). Our analyses show that Bd manifests an unusually dynamic genome that may have been shaped by its association with the amphibian host. The rates of variation that we document likely explain the high levels of phenotypic variability that have been reported for Bd, and suggests that the dynamic genome of this pathogen has contributed to its success across multiple biomes and host-species.

Pathogenic fungi constitute a growing threat to both plant and animal species on a global scale. However, many features of the fungal genome that enable them to successfully adapt to infect diverse hosts and ecological niches remain cryptic, especially for newly evolved emerging lineages. In this paper, we report three novel features of genome diversity linked to pathogenicity in the emerging amphibian pathogen, Batrachochytrium dendrobatidis (Bd). Firstly, we identified widespread chromosome copy number variation (CCNV) across our lineages, with individual isolates harboring between 2 to 5 copies of each chromosome and rapid rates of CCNV occurring in culture. In addition, by using in vitro divergence of replicate lines of Bd, we showed that changes in ploidy can occur within as few as 40 generations. Secondly, we identified uneven rates of recombination across the genomes and lineages, revealing hot spots in known classes of virulence factors. Finally we identified significant evidence of diversifying selection across the secretome of Bd, and showed that selection also targets putative virulence factors. These findings add to our knowledge of genome-dynamicity and modes of evolution manifested by eukaryote microbial pathogens, and may explain the varied phenotypic responses observed in Bd.

Experimental evidence for American bullfrog (Lithobates catesbeianus) susceptibility to chytrid fungus (Batrachochytrium dendrobatidis)

The emerging fungal pathogen, Batrachochytrium dendrobatidis (Bd), has been associated with global amphibian population declines and extinctions. American bullfrogs (Lithobates catesbeianus) are widely reported to be a tolerant host and a carrier of Bd that spreads the pathogen to less tolerant hosts. Here, we examined whether bullfrogs raised from eggs to metamorphosis in outdoor mesocosms were susceptible to Bd. We experimentally exposed metamorphic juveniles to Bd in the laboratory and compared mortality rates of pathogen-exposed animals to controls (non-exposed) in two separate experiments; one using a Bd strain isolated from a Western toad and another using a strain isolated from an American bullfrog. We wanted to examine whether metamorphic bullfrogs were susceptible to either of these strains. We show that bullfrogs were susceptible to one strain of Bd and not the other. In both experiments, infection load detected in the skin decreased over time, suggesting that metamorphic bullfrogs from some populations may be inefficient long-term carriers of Bd.

Complex history of the amphibian-killing chytrid fungus revealed with genome resequencing data

Understanding the evolutionary history of microbial pathogens is critical for mitigating the impacts of emerging infectious diseases on economically and ecologically important host species. We used a genome resequencing approach to resolve the evolutionary history of an important microbial pathogen, the chytrid Batrachochytrium dendrobatidis (Bd), which has been implicated in amphibian declines worldwide. We sequenced the genomes of 29 isolates of Bd from around the world, with an emphasis on North, Central, and South America because of the devastating effect that Bd has had on amphibian populations in the New World. We found a substantial amount of evolutionary complexity in Bd with deep phylogenetic diversity that predates observed global amphibian declines. By investigating the entire genome, we found that even the most recently evolved Bd clade (termed the global panzootic lineage) contained more genetic variation than previously reported. We also found dramatic differences among isolates and among genomic regions in chromosomal copy number and patterns of heterozygosity, suggesting complex and heterogeneous genome dynamics. Finally, we report evidence for selection acting on the Bd genome, supporting the hypothesis that protease genes are important in evolutionary transitions in this group. Bd is considered an emerging pathogen because of its recent effects on amphibians, but our data indicate that it has a complex evolutionary history that predates recent disease outbreaks. Therefore, it is important to consider the contemporary effects of Bd in a broader evolutionary context and identify specific mechanisms that may have led to shifts in virulence in this system.