Reframing wildlife disease management problems with decision analysis.

Contemporary wildlife disease management is complex because managers need to respond to a wide range of stakeholders, multiple uncertainties, and difficult trade-offs that characterize the interconnected challenges of today. Despite general acknowledgment of these complexities, managing wildlife disease tends to be framed as a scientific problem, in which the major challenge is lack of knowledge. The complex and multifactorial process of decision-making is collapsed into a scientific endeavor to reduce uncertainty. As a result, contemporary decision-making may be oversimplified, rely on simple heuristics, and fail to account for the broader legal, social, and economic context in which the decisions are made. Concurrently, scientific research on wildlife disease may be distant from this decision context, resulting in information that may not be directly relevant to the pertinent management questions. We propose reframing wildlife disease management challenges as decision problems and addressing them with decision analytical tools to divide the complex problems into more cognitively manageable elements. In particular, structured decision-making has the potential to improve the quality, rigor, and transparency of decisions about wildlife disease in a variety of systems. Examples of management of severe acute respiratory syndrome coronavirus 2, white-nose syndrome, avian influenza, and chytridiomycosis illustrate the most common impediments to decision-making, including competing objectives, risks, prediction uncertainty, and limited resources.

Replanteamiento del manejo de problemas por enfermedades de fauna mediante el análisis de decisiones Resumen El manejo actual de las enfermedades de la fauna es complejo debido a que los gestores necesitan responder a una amplia gama de actores, varias incertidumbres y compensaciones difíciles que caracterizan los retos interconectados del día de hoy. A pesar de que en general se reconocen estas complejidades, el manejo de las enfermedades tiende a plantearse como un problema científico en el que el principal obstáculo es la falta de conocimiento. El proceso complejo y multifactorial de la toma decisiones está colapsado dentro de un esfuerzo científico para reducir la incertidumbre. Como resultado de esto, las decisiones contemporáneas pueden estar simplificadas en exceso, depender de métodos heurísticos simples y no considerar el contexto legal, social y económico más amplio en el que se toman las decisiones. De manera paralela, las investigaciones científicas sobre las enfermedades de la fauna pueden estar lejos de este contexto de decisiones, lo que deriva en información que puede no ser directamente relevante para las preguntas pertinentes de manejo. Proponemos replantear los obstáculos para el manejo de enfermedades de fauna como problemas de decisión y abordarlos con herramientas analíticas de decisión para dividir los problemas complejos en elementos más manejables de manera cognitiva. En particular, las decisiones estructuradas tienen el potencial de mejorar la calidad, el rigor y la transparencia de las decisiones sobre las enfermedades de la fauna en una variedad de sistemas. Ejemplos como el manejo del coronavirus del síndrome de respiración agudo tipo 2, el síndrome de nariz blanca, la influenza aviar y la quitridiomicosis ilustran los impedimentos más comunes para la toma de decisiones, incluyendo los objetivos en competencia, riesgos, incertidumbre en las predicciones y recursos limitados.

Accommodating the role of site memory in dynamic species distribution models.

First-order dynamic occupancy models (FODOMs) are a class of state-space model in which the true state (occurrence) is observed imperfectly. An important assumption of FODOMs is that site dynamics only depend on the current state and that variations in dynamic processes are adequately captured with covariates or random effects. However, it is often difficult to understand and/or measure the covariates that generate ecological data, which are typically spatiotemporally correlated. Consequently, the non-independent error structure of correlated data causes underestimation of parameter uncertainty and poor ecological inference. Here, we extend the FODOM framework with a second-order Markov process to accommodate site memory when covariates are not available. Our modeling framework can be used to make reliable inference about site occupancy, colonization, extinction, turnover, and detection probabilities. We present a series of simulations to illustrate the data requirements and model performance. We then applied our modeling framework to 13 yr of data from an amphibian community in southern Arizona, USA. In this analysis, we found residual temporal autocorrelation of population processes for most species, even after accounting for long-term drought dynamics. Our approach represents a valuable advance in obtaining inference on population dynamics, especially as they relate to metapopulations.

Assessing amphibian disease risk across tropical streams while accounting for imperfect pathogen detection.

Ecologists studying emerging wildlife diseases need to confront the realism of imperfect pathogen detection across heterogeneous habitats to aid in conservation decisions. For example, spatial risk assessments of amphibian disease caused by Batrachochytrium dendrobatidis (Bd) has largely ignored imperfect pathogen detection across sampling sites. Because changes in pathogenicity and host susceptibility could trigger recurrent population declines, it is imperative to understand how pathogen prevalence and occupancy vary across environmental gradients. Here, we assessed how Bd occurrence, prevalence, and infection intensity in a diverse Neotropical landscape vary across streams in relation to abiotic and biotic predictors using a hierarchical Bayesian model that accounts for imperfect Bd detection caused by qPCR error. Our model indicated that the number of streams harboring Bd-infected frogs is higher than observed, with Bd likely being present at ~ 43% more streams than it was detected. We found that terrestrial-breeders captured along streams had higher Bd prevalence, but lower infection intensity, than aquatic-breeding species. We found a positive relationship between Bd occupancy probability and stream density, and a negative relationship between Bd occupancy probability and amphibian local richness. Forest cover was a weak predictor of Bd occurrence and infection intensity. Finally, we provide estimates for the minimum number of amphibian captures needed to determine the presence of Bd at a given site where Bd occurs, thus, providing guidence for cost-effective disease risk monitoring programs.

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

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

Tropical snake diversity collapses after widespread amphibian loss.

Biodiversity is declining at unprecedented rates worldwide. Yet cascading effects of biodiversity loss on other taxa are largely unknown because baseline data are often unavailable. We document the collapse of a Neotropical snake community after the invasive fungal pathogen Batrachochytrium dendrobatidis caused a chytridiomycosis epizootic leading to the catastrophic loss of amphibians, a food source for snakes. After mass mortality of amphibians, the snake community contained fewer species and was more homogeneous across the study site, with several species in poorer body condition, despite no other systematic changes in the environment. The demise of the snake community after amphibian loss demonstrates the repercussive and often unnoticed consequences of the biodiversity crisis and calls attention to the invisible declines of rare and data-deficient species.

Disease-structured N-mixture models: A practical guide to model disease dynamics using count data.

Obtaining inferences on disease dynamics (e.g., host population size, pathogen prevalence, transmission rate, host survival probability) typically requires marking and tracking individuals over time. While multistate mark-recapture models can produce high-quality inference, these techniques are difficult to employ at large spatial and long temporal scales or in small remnant host populations decimated by virulent pathogens, where low recapture rates may preclude the use of mark-recapture techniques. Recently developed N-mixture models offer a statistical framework for estimating wildlife disease dynamics from count data. N-mixture models are a type of state-space model in which observation error is attributed to failing to detect some individuals when they are present (i.e., false negatives). The analysis approach uses repeated surveys of sites over a period of population closure to estimate detection probability. We review the challenges of modeling disease dynamics and describe how N-mixture models can be used to estimate common metrics, including pathogen prevalence, transmission, and recovery rates while accounting for imperfect host and pathogen detection. We also offer a perspective on future research directions at the intersection of quantitative and disease ecology, including the estimation of false positives in pathogen presence, spatially explicit disease-structured N-mixture models, and the integration of other data types with count data to inform disease dynamics. Managers rely on accurate and precise estimates of disease dynamics to develop strategies to mitigate pathogen impacts on host populations. At a time when pathogens pose one of the greatest threats to biodiversity, statistical methods that lead to robust inferences on host populations are critically needed for rapid, rather than incremental, assessments of the impacts of emerging infectious diseases.

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.

Eco-evolutionary rescue promotes host-pathogen coexistence.

Emerging infectious pathogens are responsible for some of the most severe host mass mortality events in wild populations. Yet, effective pathogen control strategies are notoriously difficult to identify, in part because quantifying and forecasting pathogen spread and disease dynamics is challenging. Following an outbreak, hosts must cope with the presence of the pathogen, leading to host-pathogen coexistence or extirpation. Despite decades of research, little is known about host-pathogen coexistence post-outbreak when low host abundances and cryptic species make these interactions difficult to study. Using a novel disease-structured N-mixture model, we evaluate empirical support for three host-pathogen coexistence hypotheses (source-sink, eco-evolutionary rescue, and spatial variation in pathogen transmission) in a Neotropical amphibian community decimated by Batrachochytrium dendrobatidis (Bd) in 2004. During 2010-2014, we surveyed amphibians in Parque Nacional G. D. Omar Torríjos Herrera, Coclé Province, El Copé, Panama. We found that the primary driver of host-pathogen coexistence was eco-evolutionary rescue, as evidenced by similar amphibian survival and recruitment rates between infected and uninfected hosts. Average apparent monthly survival rates of uninfected and infected hosts were both close to 96%, and the expected number of uninfected and infected hosts recruited (via immigration/reproduction) was less than one host per disease state per 20-m site. The secondary driver of host-pathogen coexistence was spatial variation in pathogen transmission as we found that transmission was highest in areas of low abundance but there was no support for the source-sink hypothesis. Our results indicate that changes in the host community (i.e., through genetic or species composition) can reduce the impacts of emerging infectious disease post-outbreak. Our disease-structured N-mixture model represents a valuable advancement for conservation managers trying to understand underlying host-pathogen interactions and provides new opportunities to study disease dynamics in remnant host populations decimated by virulent pathogens.

External Reinfection of a Fungal Pathogen Does not Contribute to Pathogen Growth.

Chytridiomycosis is an emerging infectious disease of amphibians caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), which has led to devastating declines in amphibian populations worldwide. Current theory predicts that Bd infections are maintained through both reproduction on the host's skin and reinfection from sources outside of the host. To investigate the importance of external reinfection on pathogen burden, we infected captive-bred individuals of the highly susceptible Panamanian Golden Frog, Atelopus glyphus, and wild-caught glass frogs, Espadarana prosoblepon, with Bd. We housed the animals in one of three treatments: individually, in heterospecific pairs, and in conspecific pairs. For 8 weeks, we measured the Bd load and shedding rate of all frogs. We found that Atelopus had high rates of increase in both Bd load and shedding rate, but pathogen growth rates did not differ among treatments. The infection intensity of Espadarana co-housed with Atelopus was indistinguishable from those housed singly and those in conspecific pairs, despite being exposed to a large external source of Bd zoospores. Our results indicate that Bd load in both species is driven by pathogen replication within an individual, with reinfection from outside the host contributing little to the amplification of host fungal load.

Inhibition of Fungal Pathogens across Genotypes and Temperatures by Amphibian Skin Bacteria.

Symbiotic bacteria may dampen the impacts of infectious diseases on hosts by inhibiting pathogen growth. However, our understanding of the generality of pathogen inhibition by different bacterial taxa across pathogen genotypes and environmental conditions is limited. Bacterial inhibitory properties are of particular interest for the amphibian-killing fungal pathogens (Batrachochytrium dendrobatidis and Batrachochytrium salamandrivorans), for which probiotic applications as conservation strategies have been proposed. We quantified the inhibition strength of five putatively B. dendrobatidis-inhibitory bacteria isolated from woodland salamander skin against six Batrachochytrium genotypes at two temperatures (12 and 18°C). We selected six genotypes from across the Batrachochytrium phylogeny: B. salamandrivorans, B. dendrobatidis-Brazil and four genotypes of the B. dendrobatidis Global Panzootic Lineage (GPL1: JEL647, JEL404; GPL2: SRS810, JEL423). We performed 96-well plate challenge assays in a full factorial design. We detected a Batrachochytrium genotype by temperature interaction on bacterial inhibition score for all bacteria, indicating that bacteria vary in ability to inhibit Batrachochytrium depending on pathogen genotype and temperature. Acinetobacter rhizosphaerae moderately inhibited B. salamandrivorans at both temperatures (µ = 46-53%), but not any B. dendrobatidis genotypes. Chryseobacterium sp. inhibited three Batrachochytrium genotypes at both temperatures (µ = 5-71%). Pseudomonas sp. strain 1 inhibited all Batrachochytrium genotypes at 12°C and four Batrachochytrium genotypes at 18°C (µ = 5-100%). Pseudomonas sp. strain 2 and Stenotrophomonas sp. moderately to strongly inhibited all six Batrachochytrium genotypes at both temperatures (µ = 57-100%). All bacteria consistently inhibited B. salamandrivorans. Using cluster analysis of inhibition scores, we found that more closely related Batrachochytrium genotypes grouped together, suggesting that bacterial inhibition strength may be predictable based on Batrachochytrium relatedness. We conclude that bacterial inhibition capabilities change among bacterial strains, Batrachochytrium genotypes and temperatures. A comprehensive understanding of bacterial inhibitory function, across pathogen genotypes and temperatures, is needed to better predict the role of bacterial symbionts in amphibian disease ecology. For targeted conservation applications, we recommend using bacterial strains identified as strongly inhibitory as they are most likely to produce broad-spectrum antimicrobial agents at a range of temperatures.

Antifungal Bacteria on Woodland Salamander Skin Exhibit High Taxonomic Diversity and Geographic Variability.

Diverse bacteria inhabit amphibian skin; some of those bacteria inhibit growth of the fungal pathogen Batrachochytrium dendrobatidis Yet there has been no systematic survey of anti-B. dendrobatidis bacteria across localities, species, and elevations. This is important given geographic and taxonomic variations in amphibian susceptibility to B. dendrobatidis Our collection sites were at locations within the Appalachian Mountains where previous sampling had indicated low B. dendrobatidis prevalence. We determined the numbers and identities of anti-B. dendrobatidis bacteria on 61 Plethodon salamanders (37 P. cinereus, 15 P. glutinosus, 9 P. cylindraceus) via culturing methods and 16S rRNA gene sequencing. We sampled co-occurring species at three localities and sampled P. cinereus along an elevational gradient (700 to 1,000 meters above sea level [masl]) at one locality. We identified 50 anti-B. dendrobatidis bacterial operational taxonomic units (OTUs) and found that the degree of B. dendrobatidis inhibition was not correlated with relatedness. Five anti-B. dendrobatidis bacterial strains occurred on multiple amphibian species at multiple localities, but none were shared among all species and localities. The prevalence of anti-B. dendrobatidis bacteria was higher at Shenandoah National Park (NP), VA, with 96% (25/26) of salamanders hosting at least one anti-B. dendrobatidis bacterial species compared to 50% (7/14) at Catoctin Mountain Park (MP), MD, and 38% (8/21) at Mt. Rogers National Recreation Area (NRA), VA. At the individual level, salamanders at Shenandoah NP had more anti-B. dendrobatidis bacteria per individual (µ = 3.3) than those at Catoctin MP (µ = 0.8) and at Mt. Rogers NRA (µ = 0.4). All salamanders tested negative for B. dendrobatidis Anti-B. dendrobatidis bacterial species are diverse in central Appalachian Plethodon salamanders, and their distribution varied geographically. The antifungal bacterial species that we identified may play a protective role for these salamanders.IMPORTANCE Amphibians harbor skin bacteria that can kill an amphibian fungal pathogen, Batrachochytrium dendrobatidis Some amphibians die from B. dendrobatidis infection, whereas others do not. The bacteria that can kill B. dendrobatidis, called anti-B. dendrobatidis bacteria, are thought to influence the B. dendrobatidis infection outcome for the amphibian. Yet how anti-B. dendrobatidis bacterial species vary among amphibian species and populations is unknown. We determined the distribution of anti-B. dendrobatidis bacterial species among three salamander species (n = 61) sampled at three localities. We identified 50 unique anti-B. dendrobatidis bacterial species and found that all of the tested salamanders were negative for B. dendrobatidis Five anti-B. dendrobatidis bacterial species were commonly detected, suggesting a stable, functional association with these salamanders. The number of anti-B. dendrobatidis bacteria per individual varied among localities but not among co-occurring salamander species, demonstrating that environment is more influential than host factors in structuring the anti-B. dendrobatidis bacterial community. These anti-B. dendrobatidis bacteria may serve a protective function for their salamander hosts.

Disassembly of a tadpole community by a multi-host fungal pathogen with limited evidence of recovery.

Emerging infectious diseases can cause host community disassembly, but the mechanisms driving the order of species declines and extirpations following a disease outbreak are unclear. We documented the community disassembly of a Neotropical tadpole community during a chytridiomycosis outbreak, triggered by the generalist fungal pathogen, Batrachochytrium dendrobatidis (Bd). Within the first 11 months of Bd arrival, tadpole density and occupancy rapidly declined. Species rarity, in terms of tadpole occupancy and adult relative abundance, did not predict the odds of tadpole occupancy declines. But species losses were taxonomically selective, with glassfrogs (Family: Centrolenidae) disappearing the fastest and tree frogs (Family: Hylidae) and dart-poison frogs (Family: Dendrobatidae) remaining the longest. We detected biotic homogenization of tadpole communities, with post-decline communities resembling one another more strongly than pre-decline communities. The entire tadpole community was extirpated within 22 months following Bd arrival, and we found limited signs of recovery within 10 years post-outbreak. Because of imperfect species detection inherent to sampling species-rich tropical communities and the difficulty of devising a single study design protocol to sample physically complex tropical habitats, we used simulations to provide recommendations for future surveys to adequately sample diverse Neotropical communities. Our unique data set on tadpole community composition before and after Bd arrival is a valuable baseline for assessing amphibian recovery. Our results are of direct relevance to conservation managers and community ecologists interested in understanding the timing, magnitude, and consequences of disease outbreaks as emerging infectious diseases spread globally.

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.

Dead or alive? Viability of chytrid zoospores shed from live amphibian hosts.

Pathogens vary in virulence and rates of transmission because of many differences in the host, the pathogen, and their environment. The amphibian chytrid fungus, Batrachochytrium dendrobatidis (Bd), affects amphibian hosts differently, causing extinction and population declines in some species but having limited effects on others. Phenotypic differences in zoospore production rates among Bd lineages likely contribute to some of the variation observed among host responses, although no studies have quantified the viability of zoospores shed from live animals. We compared host survivorship, infection intensity, shedding rates, and zoospore viability between 2 species of endangered tropical frogs, Hylomantis lemur and Atelopus zeteki, when exposed to a highly virulent lineage of Bd (JEL 423). We applied a dye to zoospores 30 to 60 min following animal soaks, to estimate shedding rate and proportion of live zoospores shed by different species. The average infection intensity for A. zeteki was nearly 17 times higher (31,455 ± 10,103 zoospore genomic equivalents [ZGEs]) than that of H. lemur (1832 ± 1086 ZGEs), and A. zeteki died earlier than H. lemur. The proportion of viable zoospores was ~80% in both species throughout the experiment, although A. zeteki produced many more zoospores, suggesting it may play a disproportionate role in spreading disease in communities where it occurs, because the large number of viable zoospores they produce might increase infection in other species where they are reintroduced.

More than skin deep: functional genomic basis for resistance to amphibian chytridiomycosis.

The amphibian-killing chytrid fungus Batrachochytrium dendrobatidis (Bd) is one of the most generalist pathogens known, capable of infecting hundreds of species globally and causing widespread population declines and extinctions. However, some host species are seemingly unaffected by Bd, tolerating or clearing infections without clinical signs of disease. Variation in host immune responses is commonly evoked for these resistant or tolerant species, yet to date, we have no direct comparison of amphibian species responses to infection at the level of gene expression. In this study, we challenged four Central American frog species that vary in Bd susceptibility, with a sympatric virulent strain of the pathogen. We compared skin and spleen orthologous gene expression using differential expression tests and coexpression gene network analyses. We found that resistant species have reduced skin inflammatory responses and increased expression of genes involved in skin integrity. In contrast, only highly susceptible species exhibited suppression of splenic T-cell genes. We conclude that resistance to chytridiomycosis may be related to a species' ability to escape the immunosuppressive activity of the fungus. Moreover, our results indicate that within-species differences in splenic proteolytic enzyme gene expression may contribute to intraspecific variation in survival. This first comparison of amphibian functional immunogenomic architecture in response to Bd provides insights into key genetic mechanisms underlying variation in disease outcomes among amphibian species.

Fungal infection intensity and zoospore output of Atelopus zeteki, a potential acute chytrid supershedder.

Amphibians vary in their response to infection by the amphibian-killing chytrid fungus, Batrachochytrium dendrobatidis (Bd). Highly susceptible species are the first to decline and/or disappear once Bd arrives at a site. These competent hosts likely facilitate Bd proliferation because of ineffective innate and/or acquired immune defenses. We show that Atelopus zeteki, a highly susceptible species that has undergone substantial population declines throughout its range, rapidly and exponentially increases skin Bd infection intensity, achieving intensities that are several orders of magnitude greater than most other species reported. We experimentally infected individuals that were never exposed to Bd (n = 5) or previously exposed to an attenuated Bd strain (JEL427-P39; n = 3). Within seven days post-inoculation, the average Bd infection intensity was 18,213 zoospores (SE: 9,010; range: 0 to 66,928). Both average Bd infection intensity and zoospore output (i.e., the number of zoospores released per minute by an infected individual) increased exponentially until time of death (t50 = 7.018, p<0.001, t46 = 3.164, p = 0.001, respectively). Mean Bd infection intensity and zoospore output at death were 4,334,422 zoospores (SE: 1,236,431) and 23.55 zoospores per minute (SE: 22.78), respectively, with as many as 9,584,158 zoospores on a single individual. The daily percent increases in Bd infection intensity and zoospore output were 35.4% (SE: 0.05) and 13.1% (SE: 0.04), respectively. We also found that Bd infection intensity and zoospore output were positively correlated (t43 = 3.926, p<0.001). All animals died between 22 and 33 days post-inoculation (mean: 28.88; SE: 1.58). Prior Bd infection had no effect on survival, Bd infection intensity, or zoospore output. We conclude that A. zeteki, a highly susceptible amphibian species, may be an acute supershedder. Our results can inform epidemiological models to estimate Bd outbreak probability, especially as they relate to reintroduction programs.