Linking genetic and environmental factors in amphibian disease risk

Batrachochytrium dendrobatidis strain affects transcriptomic response in liver but not skin in latitudinal populations of the common toad (Bufo bufo)

Batrachochytrium dendrobatidis (Bd) is a fungal pathogen that has decimated amphibian populations worldwide for several decades. We examined the changes in gene expression in response to Bd infection in two populations of the common toad, Bufo bufo, in a laboratory experiment. We collected B. bufo eggs in southern and northern Sweden, and infected the laboratory-raised metamorphs with two strains of the global panzoonotic lineage Bd-GPL. Differential expression analysis showed significant differences between infected and control individuals in both liver and skin. The skin samples showed no discernible differences in gene expression between the two strains used, while liver samples were differentiated by strain, with one of the strains eliciting no immune response from infected toads. Immune system genes were overexpressed in skin samples from surviving infected individuals, while in liver samples the pattern was more diffuse. Splitting samples by population revealed a stronger immune response in northern individuals. Differences in transcriptional regulation between populations are particularly relevant to study in Swedish amphibians, which may have experienced varying exposure to Bd. Earlier exposure to this pathogen and subsequent adaptation or selection pressure may contribute to the survival of some populations over others, while standing genetic diversity in different populations may also affect the infection outcome.

Chytrid infections exhibit historical spread and contemporary seasonality in a declining stream-breeding frog

Species with extensive geographical ranges pose special challenges to assessing drivers of wildlife disease, necessitating collaborative and large-scale analyses. The imperilled foothill yellow-legged frog (Rana boylii) inhabits a wide geographical range and variable conditions in rivers of California and Oregon (USA), and is considered threatened by the pathogen Batrachochytrium dendrobatidis (Bd). To assess drivers of Bd infections over time and space, we compiled over 2000 datapoints from R. boylii museum specimens (collected 1897-2005) and field samples (2005-2021) spanning 9° of latitude. We observed a south-to-north spread of Bd detections beginning in the 1940s and increase in prevalence from the 1940s to 1970s, coinciding with extirpation from southern latitudes. We detected eight high-prevalence geographical clusters through time that span the species' geographical range. Field-sampled male R. boylii exhibited the highest prevalence, and juveniles sampled in autumn exhibited the highest loads. Bd infection risk was highest in lower elevation rain-dominated watersheds, and with cool temperatures and low stream-flow conditions at the end of the dry season. Through a holistic assessment of relationships between infection risk, geographical context and time, we identify the locations and time periods where Bd mitigation and monitoring will be critical for conservation of this imperilled species.

Conservation genomics of an endangered montane amphibian reveals low population structure, low genomic diversity and selection pressure from disease

Wildlife diseases are a major global threat to biodiversity. Boreal toads (Anaxyrus [Bufo] boreas) are a state-endangered species in the southern Rocky Mountains of Colorado and New Mexico, and a species of concern in Wyoming, largely due to lethal skin infections caused by the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). We performed conservation and landscape genomic analyses using single nucleotide polymorphisms from double-digest, restriction site-associated DNA sequencing in combination with the development of the first boreal toad (and first North American toad) reference genome to investigate population structure, genomic diversity, landscape connectivity and adaptive divergence. Genomic diversity (π = 0.00034-0.00040) and effective population sizes (Ne = 8.9-38.4) were low, likely due to post-Pleistocene founder effects and Bd-related population crashes over the last three decades. Population structure was also low, likely due to formerly high connectivity among a higher density of geographically proximate populations. Boreal toad gene flow was facilitated by low precipitation, cold minimum temperatures, less tree canopy, low heat load and less urbanization. We found >8X more putatively adaptive loci related to Bd intensity than to all other environmental factors combined, and evidence for genes under selection related to immune response, heart development and regulation and skin function. These data suggest boreal toads in habitats with Bd have experienced stronger selection pressure from disease than from other, broad-scale environmental variations. These findings can be used by managers to conserve and recover the species through actions including reintroduction and supplementation of populations that have declined due to Bd.

Importance of Genetic-Fitness Correlations for the Conservation of Amphibians

Endangered animals suffer from isolation of their habitats. Isolation leads to a reduction in population size as well as a decrease in genetic diversity and a concomitant increase in the risk of extinction. Amphibians are the most endangered vertebrate class. Besides habitat loss, fragmentation and isolation, amphibians are threatened by emerging diseases e.g., chytrid fungus or Ranavirus. By employing experiments, researchers investigate whether changes in genetic diversity within or among isolated populations affect amphibian fitness. While genetic diversity estimates are based on molecular markers, typically microsatellites, fitness is mostly measured as tadpole performance in rearing experiments often under varying environmental conditions. Tadpole performances (e.g., body mass, growth rate and survival) have been found to be negatively affected by low genetic diversity, as several studies have found a positive association between genetic diversity and these fitness traits. Moreover, infection with pathogens also seems to be more likely in individuals or populations with lower genetic diversity. Overall, these genetic-fitness correlations seem to be more pronounced or detectable in smaller, declining populations but not in larger populations. Genomic studies, which sample a larger fraction of the genome, are still scarce in the conservation genetic literature on amphibians. These are likely to increase in upcoming years and may reveal adaptive variants that protect against dangerous pathogens or environmental changes. Altogether, genetic-fitness correlation studies should be a priority in order to develop effective management plans for the genetic rescue of isolated, imperilled amphibian populations.

Recovered frog populations coexist with endemic Batrachochytrium dendrobatidis despite load-dependent mortality

Novel infectious diseases, particularly those caused by fungal pathogens, pose considerable risks to global biodiversity. The amphibian chytrid fungus (Batrachochytrium dendrobatidis, Bd) has demonstrated the scale of the threat, having caused the greatest recorded loss of vertebrate biodiversity attributable to a pathogen. Despite catastrophic declines on several continents, many affected species have experienced population recoveries after epidemics. However, the potential ongoing threat of endemic Bd in these recovered or recovering populations is still poorly understood. We investigated the threat of endemic Bd to frog populations that recovered after initial precipitous declines, focusing on the endangered rainforest frog Mixophyes fleayi. We conducted extensive field surveys over 4 years at three independent sites in eastern Australia. First, we compared Bd infection prevalence and infection intensities within frog communities to reveal species-specific infection patterns. Then, we analyzed mark-recapture data of M. fleayi to estimate the impact of Bd infection intensity on apparent mortality rates and Bd infection dynamics. We found that M. fleayi had lower infection intensities than sympatric frogs across the three sites, and cleared infections at higher rates than they gained infections throughout the study period. By incorporating time-varying individual infection intensities, we show that healthy M. fleayi populations persist despite increased apparent mortality associated with infrequent high Bd loads. Infection dynamics were influenced by environmental conditions, with Bd prevalence, infection intensity, and rates of gaining infection associated with lower temperatures and increased rainfall. However, mortality remained constant year-round despite these fluctuations in Bd infections, suggesting major mortality events did not occur over the study period. Together, our results demonstrate that while Bd is still a potential threat to recovered populations of M. fleayi, high rates of clearing infections and generally low average infection loads likely minimize mortality caused by Bd. Our results are consistent with pathogen resistance contributing to the coexistence of M. fleayi with endemic Bd. We emphasize the importance of incorporating infection intensity into disease models rather than infection status alone. Similar population and infection dynamics likely exist within other recovered amphibian-Bd systems around the globe, promising longer-term persistence in the face of endemic chytridiomycosis.

Recovered frog populations coexist with endemic Batrachochytrium dendrobatidis despite load-dependent mortality

Novel infectious diseases, particularly those caused by fungal pathogens, pose considerable risks to global biodiversity. The amphibian chytrid fungus (Batrachochytrium dendrobatidis, Bd) has demonstrated the scale of the threat, having caused the greatest recorded loss of vertebrate biodiversity attributable to a pathogen. Despite catastrophic declines on several continents, many affected species have experienced population recoveries after epidemics. However, the potential ongoing threat of endemic Bd in these recovered or recovering populations is still poorly understood. We investigated the threat of endemic Bd to frog populations that recovered after initial precipitous declines, focusing on the endangered rainforest frog Mixophyes fleayi. We conducted extensive field surveys over 4 years at three independent sites in eastern Australia. First, we compared Bd infection prevalence and infection intensities within frog communities to reveal species-specific infection patterns. Then, we analyzed mark-recapture data of M. fleayi to estimate the impact of Bd infection intensity on apparent mortality rates and Bd infection dynamics. We found that M. fleayi had lower infection intensities than sympatric frogs across the three sites, and cleared infections at higher rates than they gained infections throughout the study period. By incorporating time-varying individual infection intensities, we show that healthy M. fleayi populations persist despite increased apparent mortality associated with infrequent high Bd loads. Infection dynamics were influenced by environmental conditions, with Bd prevalence, infection intensity, and rates of gaining infection associated with lower temperatures and increased rainfall. However, mortality remained constant year-round despite these fluctuations in Bd infections, suggesting major mortality events did not occur over the study period. Together, our results demonstrate that while Bd is still a potential threat to recovered populations of M. fleayi, high rates of clearing infections and generally low average infection loads likely minimize mortality caused by Bd. Our results are consistent with pathogen resistance contributing to the coexistence of M. fleayi with endemic Bd. We emphasize the importance of incorporating infection intensity into disease models rather than infection status alone. Similar population and infection dynamics likely exist within other recovered amphibian-Bd systems around the globe, promising longer-term persistence in the face of endemic chytridiomycosis.

Recovered frog populations coexist with endemic Batrachochytrium dendrobatidis despite load-dependent mortality

Novel infectious diseases, particularly those caused by fungal pathogens, pose considerable risks to global biodiversity. The amphibian chytrid fungus (Batrachochytrium dendrobatidis, Bd) has demonstrated the scale of the threat, having caused the greatest recorded loss of vertebrate biodiversity attributable to a pathogen. Despite catastrophic declines on several continents, many affected species have experienced population recoveries after epidemics. However, the potential ongoing threat of endemic Bd in these recovered or recovering populations is still poorly understood. We investigated the threat of endemic Bd to frog populations that recovered after initial precipitous declines, focusing on the endangered rainforest frog Mixophyes fleayi. We conducted extensive field surveys over 4 years at three independent sites in eastern Australia. First, we compared Bd infection prevalence and infection intensities within frog communities to reveal species-specific infection patterns. Then, we analyzed mark-recapture data of M. fleayi to estimate the impact of Bd infection intensity on apparent mortality rates and Bd infection dynamics. We found that M. fleayi had lower infection intensities than sympatric frogs across the three sites, and cleared infections at higher rates than they gained infections throughout the study period. By incorporating time-varying individual infection intensities, we show that healthy M. fleayi populations persist despite increased apparent mortality associated with infrequent high Bd loads. Infection dynamics were influenced by environmental conditions, with Bd prevalence, infection intensity, and rates of gaining infection associated with lower temperatures and increased rainfall. However, mortality remained constant year-round despite these fluctuations in Bd infections, suggesting major mortality events did not occur over the study period. Together, our results demonstrate that while Bd is still a potential threat to recovered populations of M. fleayi, high rates of clearing infections and generally low average infection loads likely minimize mortality caused by Bd. Our results are consistent with pathogen resistance contributing to the coexistence of M. fleayi with endemic Bd. We emphasize the importance of incorporating infection intensity into disease models rather than infection status alone. Similar population and infection dynamics likely exist within other recovered amphibian-Bd systems around the globe, promising longer-term persistence in the face of endemic chytridiomycosis.

Linking pathogen-microbiome-host interactions to explain amphibian population dynamics

Symbiotic interactions can determine the evolutionary trajectories of host species, influencing genetic variation through selection and changes in demography. In the context of strong selective pressures such as those imposed by infectious diseases, symbionts providing defences could contribute to increase host fitness upon pathogen emergence. Here, we generated genome-wide data of an amphibian species to find evidence of evolutionary pressures driven by two skin symbionts: a batrachochytrid fungal pathogen and an antifungal bacterium. Using demographic modelling, we found evidence of decreased effective population size, probably due to pathogen infections. Additionally, we investigated host genetic associations with infection status, antifungal bacterium abundance and overall microbiome diversity using structural equation models. We uncovered relatively lower nucleotide diversity in infected frogs and potential heterozygote advantage to recruit the candidate beneficial symbiont and fight infections. Our models indicate that environmental conditions have indirect effects on symbiont abundance through both host body traits and microbiome diversity. Likewise, we uncovered a potential offsetting effect among host heterozygosity-fitness correlations, plausibly pointing to different ecological and evolutionary processes among the three species due to dynamic interactions. Our findings revealed that evolutionary pressures not only arise from the pathogen but also from the candidate beneficial symbiont, and both interactions shape the genetics of the host. Our results advance knowledge about multipartite symbiotic relationships and provide a framework to model ecological and evolutionary dynamics in wild populations. Finally, our study approach can be applied to inform conservation actions such as bioaugmentation strategies for other imperilled amphibians affected by infectious diseases.

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

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

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

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

Top @ Mohd. Tah M, Shohaimi S, Ab Ghani NI, Fletcher C. Morphometric study of Kalophrynus palmatissimus at two forest reserves in Malaysia

A research study on morphometrics of Kalophrynus palmatissimus (commonly known as Lowland Grainy Frog) at Ayer Hitam Forest Reserve (AHFR), Selangor and Pasoh Forest Reserve (PFR), Negeri Sembilan was carried out from 12 November 2016 to 13 September 2017. The study was to examine data on the morphometric traits of K. palmatissimus at the two forest reserves. 15 morphometric traits of K. palmatissimus that were taken by using vernier calipers. Frog surveys were done by using 15 and 18 nocturnal 400 m transect lines with an interval distance of 20 m at AHFR and PFR, respectively. The GPS coordinates for all frog samples were recorded to ensure the precise geographic location. In addition, five climatic data were recorded. The results showed that most morphometric traits in AHFR (n = 34) and PFR (n = 31) were positively correlated with each other. On the other hand, climatic factor, which was soil pH, had a significant positive influence on most of the morphometric traits (p < .01), except for tympanum diameter and upper eyelid width (p ≥ .05). Meanwhile, the temperature had a significantly negative influence on all morphometric traits (p < .01). General linear model (GLM) analysis showed that snout-vent length (SVL) influenced most morphometric traits (F ≤ 80.86, p < .01), except for hand length (HAL: F = 0.299, p > .05). Later, it was found that the snout-vent length of K. palmatissimus at AHFR was slightly larger than at PFR (AHFR: μ = 37.00 mm, SE = 1.16 c.f. PFR: μ = 30.29 mm, SE = 1.07). It showed that there were variations in morphometric traits of K. palmatissimus at AHFR and PFR. From PCA analysis, morphometric traits are grouped into two components for AHFR and PFR, respectively. In AHFR, head length, eye diameter, head width, internarial distance, interorbital distance, forearm length, tibia length, foot length, and thigh length were strongly correlated, while snout length and eye-nostril distance were strongly correlated. In PFR, eye diameter, head width, internarial distance, interorbital distance, foot length, and thigh length were strongly correlated, though snout length and eye-nostril distance were strongly correlated, hence, suggested that all morphometric traits grow simultaneously in K. palmatissimus with eye-nostril distance (EN), and snout length (SL) growing almost simultaneously at AHFR (r = .91) and PFR (r = .97). There is still a lack of available information regarding the distribution and morphometric studies of K. palmatissimus in Malaysia, especially at AHFR and PFR. This study showed 15 different morphometric traits of K. palmatisssimus between AHFR and PFR, with K. palmatissimus at AHFR were found to be slightly larger than at PFR.

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

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

Metapopulation Structure Predicts Population Dynamics in the Cakile maritima-Alternaria brassicicola Host-Pathogen Interaction

AbstractIn symbiotic interactions, spatiotemporal variation in the distribution or population dynamics of one species represents spatial and temporal heterogeneity of the landscape for the other. Such interdependent demographic dynamics result in situations where the relative importance of biotic and abiotic factors in determining ecological processes is complicated to decipher. Using a detailed survey of three metapopulations of the succulent plant Cakile maritima and the necrotrophic fungus Alternaria brassicicola located along the southeastern Australian coast, we developed a series of statistical analyses-namely, synchrony analysis, patch occupancy dynamics, and a spatially explicit metapopulation model-to understand how habitat quality, weather conditions, dispersal, and spatial structure determine metapopulation dynamics. Climatic conditions are important drivers, likely explaining the high synchrony among populations. Host availability, landscape features facilitating dispersal, and habitat conditions also impact the occurrence and spread of disease. Overall, we show that the collection of extensive data on host and pathogen population dynamics, in combination with spatially explicit epidemiological modeling, makes it possible to accurately predict disease dynamics-even when there is extreme variability in host population dynamics. Finally, we discuss the importance of genetic information for predicting demographic dynamics in this pathosystem.

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

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

Does genetic diversity protect host populations from parasites? A meta-analysis across natural and agricultural systems

If parasites transmit more readily between closely related hosts, then parasite burdens should decrease with increased genetic diversity of host populations. This important hypothesis is often accepted at face value—notorious epidemics of crop monocultures testify to the vulnerability of host populations that have been purged of diversity. Yet the relationship between genetic diversity and parasitism likely varies across contexts, differing between crop and noncrop hosts and between experimental and natural host populations. Here, we used a meta‐analytic approach to ask if host diversity confers protection against parasites over the range of contexts in which it has been tested. We synthesized the results of 102 studies, comprising 2004 effect sizes representing a diversity of approaches and host‐parasite systems. Our results validate a protective effect of genetic diversity, while revealing significant variation in its strength across biological and empirical contexts. In experimental host populations, genetic diversity reduces parasitism by ∼20% for noncrop hosts and by ∼50% for crop hosts. In contrast, observational studies of natural host populations show no consistent relationship between genetic diversity and parasitism, with both strong negative and positive correlations reported. This result supports the idea that, if parasites preferentially attack close relatives, the correlation of genetic diversity with parasitism could be positive or negative depending upon the potential for host populations to evolve in response to parasite selection. Taken together, these results reinforce genetic diversity as a priority for both conservation and agriculture and emphasize the challenges inherent to drawing comparisons between controlled experimental populations and dynamic natural populations.

Preparing for invasion: Assessing risk of infection by chytrid fungi in southeastern plethodontid salamanders

Understanding the responses of naïve communities to the invasion of multihost pathogens requires accurate estimates of susceptibility across taxa. In the Americas, the likely emergence of a second amphibian pathogenic fungus (Batrachochytrium salamandrivorans, Bsal) calls for new ways of prioritizing disease mitigation among species due to the high diversity of naïve hosts with prior B. dendrobatidis (Bd) infections. Here, we applied the concept of pathogenic potential to quantify the virulence of chytrid fungi on naïve amphibians and evaluate species for conservation efforts in the event of an outbreak. The benefit of this measure is that it combines and summarizes the variation in disease effects into a single numerical index, allowing for comparisons across species, populations or groups of individuals that may inherently exhibit differences in susceptibility. As a proof of concept, we obtained standardized responses of disease severity by performing experimental infections with Bsal on five plethodontid salamanders from southeastern United States. Four out of five species carried natural infections of Bd at the start of the experiments. We showed that Bsal exhibited its highest value of pathogenic potential in a species that is already declining (Desmognathus auriculatus). We find that this index provides additional information beyond the standard measures of disease prevalence, intensity, and mortality, because it leveraged these disease parameters within each categorical group. Scientists and practitioners could use this measure to justify research, funding, trade, or conservation measures.

Chytridiomycosis-induced mortality in a threatened anuran

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

Mechanisms underlying host persistence following amphibian disease emergence determine appropriate management strategies

Emerging infectious diseases have caused many species declines, changes in communities and even extinctions. There are also many species that persist following devastating declines due to disease. The broad mechanisms that enable host persistence following declines include evolution of resistance or tolerance, changes in immunity and behaviour, compensatory recruitment, pathogen attenuation, environmental refugia, density-dependent transmission and changes in community composition. Here we examine the case of chytridiomycosis, the most important wildlife disease of the past century. We review the full breadth of mechanisms allowing host persistence, and synthesise research on host, pathogen, environmental and community factors driving persistence following chytridiomycosis-related declines and overview the current evidence and the information required to support each mechanism. We found that for most species the mechanisms facilitating persistence have not been identified. We illustrate how the mechanisms that drive long-term host population dynamics determine the most effective conservation management strategies. Therefore, understanding mechanisms of host persistence is important because many species continue to be threatened by disease, some of which will require intervention. The conceptual framework we describe is broadly applicable to other novel disease systems.

Connectivity over a disease risk gradient enables recovery of rainforest frogs

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

Sustained immune activation is associated with susceptibility to the amphibian chytrid fungus

The disease chytridiomycosis caused by the fungus Bd has devastated amphibian populations worldwide. Functional genomic contributions to host susceptibility remain enigmatic and vary between species and populations. We conducted experimental Bd infections in Rana yavapaiensis, a species with intraspecific variation in chytridiomycosis susceptibility, to assess the skin and spleen transcriptomic response to infection over time. We predicted that increased immune gene expression would be associated with a positive disease outcome, but we instead found that surviving frogs had significantly reduced immune gene expression compared to susceptible frogs and to uninfected controls. MHC class IIβ gene expression was also significantly higher in susceptible frogs compared to surviving frogs. Furthermore, susceptible frogs expressed a significantly larger number of distinct class IIβ alleles, demonstrating a negative correlation between class IIβ expression, functional diversity, and survival. Expression of the MHC class IIβ locus previously associated with Bd disease outcomes was a significant predictor of Bd infection intensity at early infection stages but not at late infection stages, suggesting initial MHC-linked immune processes are important for ultimate disease outcomes. We infer through disease association and phylogenetic analysis that certain MHC variants are linked to the immune expression that was negatively associated with survival, and we hypothesize that frogs that did not express these alleles could better survive infections. Our study finds that MHC expression at early and late infection stages predicts Bd infection intensity, and suggests that generating a sustained immune response against Bd may be counterproductive for surviving chytridiomycosis in this partially susceptible species.

Evolutionary principles guiding amphibian conservation

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

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

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

A review on the conservation genetic studies of Indian amphibians and their implications on developing strategies for conservation

Amphibians show a very high level of diversity and endemism and are facing global declines from the past few decades. Studies have shown that the molecular tools can be helpful in their conservation efforts. In India, more than 80% of amphibians are endemic and most show a narrow range of distribution. Most of the Indian amphibians lack information on their genetic diversity. In this study, were view the overall trend on amphibian studies in India with the specific focus on conservation genetics. Overall, of the 173 studies, only 14 dealt with the conservation of amphibians through genetic tools and five studies estimated the genetic diversity or gene structure. Here, we discuss the gaps and provide future directions on how genetic studies can be helpful in Indian amphibian conservation.

The Influence of Habitat and Phylogeny on the Skin Microbiome of Amphibians in Guatemala and Mexico

Microbial symbionts are increasingly recognized as playing a critical role in organismal health across a wide range of hosts. Amphibians are unique hosts in that their skin helps to regulate the exchange of water, ions, and gases, and it plays an active role in defense against pathogens through the synthesis of anti-microbial peptides. The microbiome of amphibian skin includes a diverse community of bacteria known to defend against pathogens, including the global pandemic lineage of Batrachochytrium dendrobatidis associated with mass amphibian die-offs. The relative influence of host phylogeny and environment in determining the composition of the amphibian skin microbiome remains poorly understood. We collected skin swabs from montane amphibians in Mexico and Guatemala, focusing on two genera of plethodontid salamanders and one genus of frogs. We used high throughput sequencing to characterize the skin bacterial microbiome and tested the impact of phylogeny and habitat on bacterial diversity. Our results show that phylogenetic history strongly influences the diversity and community structure of the total bacterial microbiome at higher taxonomic levels (between orders), but on lower scales (within genera and species), the effect of habitat predominates. These results add to a growing consensus that habitat exerts a strong effect on microbiome structure and composition, particularly at shallow phylogenetic scales.

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

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

Pathogen Dynamics in an Invasive Frog Compared to Native Species

Emerging infectious diseases threaten the survival of wildlife populations and species around the world. In particular, amphibians are experiencing population declines and species extinctions primarily in response to two pathogens, the fungus Batrachochytrium dendrobatidis (Bd) and the iridovirus Ranavirus (Rv). Here, we use field surveys and quantitative (q)PCR to compare infection intensity and prevalence of Bd and Rv across species and seasons on Jekyll Island, a barrier island off the coast of Georgia, USA. We collected oral and skin swabs for 1 year from four anuran species and three families, including two native hylids (Hyla cinerea and Hyla squirella), a native ranid (Rana sphenocephala), and the invasive rain frog Eleutherodactylus planirostris. Bd infection dynamics did not vary significantly over sampling months, but Rv prevalence and intensity were significantly higher in fall 2014 compared to spring 2015. Additionally, Rv prevalence and intensity were significantly higher in E. planirostris than in the other three species. Our study highlights the potential role of invasive amphibians as drivers of disease dynamics and demonstrates the importance of pathogen surveillance across multiple time periods and species to accurately capture the infectious disease landscape.

Longitudinal study of Amphibiocystidium sp. infection in a natural population of the Italian stream frog (Rana italica)

Mesomycetozoean-induced infections (order Dermocystida, genus Amphibiocystidium) in European and North American amphibians are causing alarm. To date, the pathogenicity of these parasites in field conditions has been poorly studied, and demographic consequences on amphibian populations have not been explored. In this study, an Amphibiocystidium sp. infection is reported in a natural population of the Italian stream frog (Rana italica) of Central Italy, over a 7-year period from 2008 to 2014. Light and electron microscope examinations, as well as partial 18S rDNA sequence analysis were used to characterize the parasite. Moreover, a capture-mark-recapture study was conducted to assess the frog demographics in response to infection. Negative effects of amphibiocystidiosis on individual survival and population fitness were absent throughout the sampling period, despite the high estimates of disease prevalence. This might have been due to resistance and/or tolerance strategies developed by the frogs in response to the persistence of Amphibiocystidium infection in this system. We hypothesized that in the examined R. italica population, amphibiocystidiosis is an ongoing endemic/epidemic infection. However, ecological and host-specific factors, interacting in a synergistic fashion, might be responsible for variations in the susceptibility to Amphibiocystidium infection of both conspecific populations and heterospecific individuals of R. italica.

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

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

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.

Hybrids of amphibian chytrid show high virulence in native hosts

Hybridization of parasites can generate new genotypes with high virulence. The fungal amphibian parasite Batrachochytrium dendrobatidis (Bd) hybridizes in Brazil’s Atlantic Forest, a biodiversity hotspot where amphibian declines have been linked to Bd, but the virulence of hybrid genotypes in native hosts has never been tested. We compared the virulence (measured as host mortality and infection burden) of hybrid Bd genotypes to the parental lineages, the putatively hypovirulent lineage Bd-Brazil and the hypervirulent Global Pandemic Lineage (Bd-GPL), in a panel of native Brazilian hosts. In Brachycephalus ephippium, the hybrid exceeded the virulence (host mortality) of both parents, suggesting that novelty arising from hybridization of Bd is a conservation concern. In Ischnocnema parva, host mortality in the hybrid treatment was intermediate between the parent treatments, suggesting that this species is more vulnerable to the aggressive phenotypes associated with Bd-GPL. Dendropsophus minutus showed low overall mortality, but infection burdens were higher in frogs treated with hybrid and Bd-GPL genotypes than with Bd-Brazil genotypes. Our experiment suggests that Bd hybrids have the potential to increase disease risk in native hosts. Continued surveillance is needed to track potential spread of hybrid genotypes and detect future genomic shifts in this dynamic disease system.

Genetic variation and selection of MHC class I loci differ in two congeneric frogs

Major histocompatibility complex (MHC) genes encode proteins in the acquired immune response pathway that often show distinctive selection-driven patterns in wild vertebrate populations. We examined genetic variation and signatures of selection in the MHC class I alpha 1 (A1)- and alpha 2 (A2)-domain encoding exons of two frog congeners [Agalychnis callidryas (n = 20) and A. lemur (n = 20)] from a single locality in Panama. We also investigated how historical demographic processes may have impacted MHC genetic diversity by analyzing a neutral mitochondrial marker. We found that both MHC domains were highly variable in both species, with both species likely expressing three loci. Our analyses revealed different signatures of selection between the two species, most notably that the A. callidryas A2 domain had experienced positive selection while the A2 domain of A. lemur had not. Diversifying selection acted on the same number of A1 and A2 allelic lineages, but on a higher percentage of A1 sites compared to A2 sites. Neutrality tests of mitochondrial haplotypes predominately indicated that the two species were at genetic equilibrium when the samples were collected. In addition, two historical tests of demography indicated both species have had relatively stable population sizes over the past 100,000 years; thus large population size changes are unlikely to have greatly influenced MHC diversity in either species during this time period. In conclusion, our results suggest that the impact of selection on MHC diversity varied between these two closely related species, likely due to a combination of distinct ecological conditions and past pathogenic pressures.

Joint effects of habitat, zooplankton, host stage structure and diversity on amphibian chytrid

Why does the severity of parasite infection differ dramatically across habitats? This question remains challenging to answer because multiple correlated pathways drive disease. Here, we examined habitat-disease links through direct effects on parasites and indirect effects on parasite predators (zooplankton), host diversity and key life stages of hosts. We used a case study of amphibian hosts and the chytrid fungus, Batrachochytrium dendrobatidis, in a set of permanent and ephemeral alpine ponds. A field experiment showed that ultraviolet radiation (UVR) killed the free-living infectious stage of the parasite. Yet, permanent ponds with more UVR exposure had higher infection prevalence. Two habitat-related indirect effects worked together to counteract parasite losses from UVR: (i) UVR reduced the density of parasite predators and (ii) permanent sites fostered multi-season host larvae that fuelled parasite production. Host diversity was unlinked to hydroperiod or UVR but counteracted parasite gains; sites with higher diversity of host species had lower prevalence of infection. Thus, while habitat structure explained considerable variation in infection prevalence through two indirect pathways, it could not account for everything. This study demonstrates the importance of creating mechanistic, food web-based links between multiple habitat dimensions and disease.

Using a Bayesian network to clarify areas requiring research in a host-pathogen system

Bayesian network analyses can be used to interactively change the strength of effect of variables in a model to explore complex relationships in new ways. In doing so, they allow one to identify influential nodes that are not well studied empirically so that future research can be prioritized. We identified relationships in host and pathogen biology to examine disease-driven declines of amphibians associated with amphibian chytrid fungus (Batrachochytrium dendrobatidis). We constructed a Bayesian network consisting of behavioral, genetic, physiological, and environmental variables that influence disease and used them to predict host population trends. We varied the impacts of specific variables in the model to reveal factors with the most influence on host population trend. The behavior of the nodes (the way in which the variables probabilistically responded to changes in states of the parents, which are the nodes or variables that directly influenced them in the graphical model) was consistent with published results. The frog population had a 49% probability of decline when all states were set at their original values, and this probability increased when body temperatures were cold, the immune system was not suppressing infection, and the ambient environment was conducive to growth of B. dendrobatidis. These findings suggest the construction of our model reflected the complex relationships characteristic of host-pathogen interactions. Changes to climatic variables alone did not strongly influence the probability of population decline, which suggests that climate interacts with other factors such as the capacity of the frog immune system to suppress disease. Changes to the adaptive immune system and disease reservoirs had a large effect on the population trend, but there was little empirical information available for model construction. Our model inputs can be used as a base to examine other systems, and our results show that such analyses are useful tools for reviewing existing literature, identifying links poorly supported by evidence, and understanding complexities in emerging infectious-disease systems.

Phenotypes in phylogeography: Species' traits, environmental variation, and vertebrate diversification

Almost 30 y ago, the field of intraspecific phylogeography laid the foundation for spatially explicit and genealogically informed studies of population divergence. With new methods and markers, the focus in phylogeography shifted to previously unrecognized geographic genetic variation, thus reducing the attention paid to phenotypic variation in those same diverging lineages. Although phenotypic differences among lineages once provided the main data for studies of evolutionary change, the mechanisms shaping phenotypic differentiation and their integration with intraspecific genetic structure have been underexplored in phylogeographic studies. However, phenotypes are targets of selection and play important roles in species performance, recognition, and diversification. Here, we focus on three questions. First, how can phenotypes elucidate mechanisms underlying concordant or idiosyncratic responses of vertebrate species evolving in shared landscapes? Second, what mechanisms underlie the concordance or discordance of phenotypic and phylogeographic differentiation? Third, how can phylogeography contribute to our understanding of functional phenotypic evolution? We demonstrate that the integration of phenotypic data extends the reach of phylogeography to explain the origin and maintenance of biodiversity. Finally, we stress the importance of natural history collections as sources of high-quality phenotypic data that span temporal and spatial axes.

Low resistance to chytridiomycosis in direct-developing amphibians

Host-generalist pathogens sporadically infect naive hosts, potentially triggering epizootics. The waterborne fungus Batrachochytrium dendrobatidis (Bd) is linked to declines of hundreds of amphibian species with aquatic larvae. Although several population declines and extinctions attributed to Bd have been reported among cryptic species undergoing direct development away from water, epidemiological studies focused on these terrestrial frogs are lacking. Our field data support that terrestrial direct-developing hosts are less exposed to Bd during their ontogeny than species with aquatic larvae, and thus they might lack adaptive responses against waterborne chytrids. Using controlled laboratory experiments, we exposed wild-caught amphibian species with terrestrial and aquatic life histories to Bd and found that direct developers showed more rapid increases in infection loads and experienced higher mortality rates than species with aquatic larvae. Our findings provide novel information about host responses to generalist pathogens and specifically show that our focal direct developing species have low resistance to Bd infections. Finally, our results underscore that we should not ignore Bd as a potential threat to direct developing species simply because they are less exposed to Bd in nature; instead future amphibian conservation plans should include efforts to safeguard hundreds of direct-developing amphibian species globally.

Differential patterns of Batrachochytrium dendrobatidis infection in relict amphibian populations following severe disease-associated declines

Global amphibian biodiversity has declined dramatically in the past 4 decades, and many amphibian species have declined to near extinction as a result of emergence of the amphibian chytrid fungus, Batrachochytrium dendrobatidis (Bd). However, persistent or recovering populations of several amphibian species have recently been rediscovered, and such populations may illustrate how amphibian species that are highly susceptible to chytridiomycosis may survive in the presence of Bd. We conducted field surveys for Bd infection in 7 species of Costa Rican amphibians (all species that have declined to near extinction but for which isolated populations persist) to characterize infection profiles in highly Bd-susceptible amphibians post-decline. We found highly variable patterns in infection, with some species showing low prevalence (~10%) and low infection intensity and others showing high infection prevalence (>80%) and either low or high infection intensity. Across sites, infection rates were negatively associated with mean annual precipitation, and infection intensity across sites was negatively associated with mean average temperatures. Our results illustrate that even the most Bd-susceptible amphibians can persist in Bd-enzootic ecosystems, and that multiple ecological or evolutionary mechanisms likely exist for host-pathogen co-existence between Bd and the most Bd-susceptible amphibian species. Continued monitoring of these populations is necessary to evaluate population trends (continuing decline, stability, or population growth). These results should inform efforts to mitigate impacts of Bd on amphibians in the field.

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

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

Cryptic chytridiomycosis linked to climate and genetic variation in amphibian populations of the southeastern United States

North American amphibians have recently been impacted by two major emerging pathogens, the fungus Batrachochytrium dendrobatidis (Bd) and iridoviruses in the genus Ranavirus (Rv). Environmental factors and host genetics may play important roles in disease dynamics, but few studies incorporate both of these components into their analyses.

Here, we investigated the role of environmental and genetic factors in driving Bd and Rv infection prevalence and severity in a biodiversity hot spot, the southeastern United States.

We used quantitative PCR to characterize Bd and Rv dynamics in natural populations of three amphibian species: Notophthalmus perstriatus, Hyla squirella and Pseudacris ornata. We combined pathogen data, genetic diversity metrics generated from neutral markers, and environmental variables into general linear models to evaluate how these factors impact infectious disease dynamics. Occurrence, prevalence and intensity of Bd and Rv varied across species and populations, but only one species, Pseudacris ornata, harbored high Bd intensities in the majority of sampled populations. Genetic diversity and climate variables both predicted Bd prevalence, whereas climatic variables alone predicted infection intensity.

We conclude that Bd is more abundant in the southeastern United States than previously thought and that genetic and environmental factors are both important for predicting amphibian pathogen dynamics. Incorporating both genetic and environmental information into conservation plans for amphibians is necessary for the development of more effective management strategies to mitigate the impact of emerging infectious diseases.

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

Amphibians have been affected globally by the disease chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), and we are just now beginning to understand how immunogenetic variability contributes to disease susceptibility. Lineages of an expressed major histocompatibility complex (MHC) class II locus involved in acquired immunity are associated with chytridiomycosis susceptibility in controlled laboratory challenge assays. Here, we extend these findings to natural populations that vary both in exposure and response to Bd. We find that MHC alleles and supertypes associated with Bd survival in the field show a molecular signal of positive selection, while those associated with susceptibility do not, supporting the hypothesis that heritable Bd tolerance is rapidly evolving. We compare MHC supertypes to neutral loci to demonstrate where selection versus demography is shaping MHC variability. One population with Bd tolerance in nature shows a significant signal of directional selection for the same allele (allele Q) that was significantly associated with survival in an earlier laboratory study. Our findings indicate that selective pressure for Bd survival drives rapid immunogenetic adaptation in some natural populations, despite differences in environment and demography. Our field-based analysis of immunogenetic variation confirms that natural amphibian populations have the evolutionary potential to adapt to chytridiomycosis.

Prediction and Prevention of Parasitic Diseases Using a Landscape Genomics Framework

Substantial heterogeneity exists in the dispersal, distribution and transmission of parasitic species. Understanding and predicting how such features are governed by the ecological variation of landscape they inhabit is the central goal of spatial epidemiology. Genetic data can further inform functional connectivity among parasite, host and vector populations in a landscape. Gene flow correlates with the spread of epidemiologically relevant phenotypes among parasite and vector populations (e.g., virulence, drug and pesticide resistance), as well as invasion and re-invasion risk where parasite transmission is absent due to current or past intervention measures. However, the formal integration of spatial and genetic data ('landscape genetics') is scarcely ever applied to parasites. Here, we discuss the specific challenges and practical prospects for the use of landscape genetics and genomics to understand the biology and control of parasitic disease and present a practical framework for doing so.

Reduced immune function predicts disease susceptibility in frogs infected with a deadly fungal pathogen

The relationship between amphibian immune function and disease susceptibility is of primary concern given current worldwide declines linked to the pathogenic fungus Batrachochytrium dendrobatidis (Bd). We experimentally infected lowland leopard frogs (Lithobates yavapaiensis) with Bd to test the hypothesis that infection causes physiological stress and stimulates humoral and cell-mediated immune function in the blood. We measured body mass, the ratio of circulating neutrophils to lymphocytes (a known indicator of physiological stress) and plasma bacterial killing ability (BKA; a measure of innate immune function). In early exposure (1-15 days post-infection), stress was elevated in Bd-positive vs. Bd-negative frogs, whereas other metrics were similar between the groups. At later stages (29-55 days post-infection), stress was increased in Bd-positive frogs with signs of chytridiomycosis compared with both Bd-positive frogs without disease signs and uninfected control frogs, which were similar to each other. Infection decreased growth during the same period, demonstrating that sustained resistance to Bd is energetically costly. Importantly, BKA was lower in Bd-positive frogs with disease than in those without signs of chytridiomycosis. However, neither group differed from Bd-negative control frogs. The low BKA values in dying frogs compared with infected individuals without disease signs suggests that complement activity might signify different immunogenetic backgrounds or gene-by-environment interactions between the host, Bd and abiotic factors. We conclude that protein complement activity might be a useful predictor of Bd susceptibility and might help to explain differential disease outcomes in natural amphibian populations.

Community Structure and Function of Amphibian Skin Microbes: An Experiment with Bullfrogs Exposed to a Chytrid Fungus

The vertebrate microbiome contributes to disease resistance, but few experiments have examined the link between microbiome community structure and disease resistance functions. Chytridiomycosis, a major cause of amphibian population declines, is a skin disease caused by the fungus, Batrachochytrium dendrobatidis (Bd). In a factorial experiment, bullfrog skin microbiota was reduced with antibiotics, augmented with an anti-Bd bacterial isolate (Janthinobacterium lividum), or unmanipulated, and individuals were then either exposed or not exposed to Bd. We found that the microbial community structure of individual frogs prior to Bd exposure influenced Bd infection intensity one week following exposure, which, in turn, was negatively correlated with proportional growth during the experiment. Microbial community structure and function differed among unmanipulated, antibiotic-treated, and augmented frogs only when frogs were exposed to Bd. Bd is a selective force on microbial community structure and function, and beneficial states of microbial community structure may serve to limit the impacts of infection.

Seasonal and ontogenetic variation of skin microbial communities and relationships to natural disease dynamics in declining amphibians

Recently, microbiologists have focused on characterizing the probiotic role of skin bacteria for amphibians threatened by the fungal disease chytridiomycosis. However, the specific characteristics of microbial diversity required to maintain health or trigger disease are still not well understood in natural populations. We hypothesized that seasonal and developmental transitions affecting susceptibility to chytridiomycosis could also alter the stability of microbial assemblages. To test our hypothesis, we examined patterns of skin bacterial diversity in two species of declining amphibians (Lithobates yavapaiensis and Eleutherodactylus coqui) affected by the pathogenic fungus Batrachochytrium dendrobatidis (Bd). We focused on two important transitions that affect Bd susceptibility: ontogenetic (from juvenile to adult) shifts in E. coqui and seasonal (from summer to winter) shifts in L. yavapaiensis. We used a combination of community-fingerprinting analyses and 16S rRNA amplicon sequencing to quantify changes in bacterial diversity and assemblage composition between seasons and developmental stages, and to investigate the relationship between bacterial diversity and pathogen load. We found that winter-sampled frogs and juveniles, two states associated with increased Bd susceptibility, exhibited higher diversity compared with summer-sampled frogs and adult individuals. Our findings also revealed that hosts harbouring higher bacterial diversity carried lower Bd infections, providing support for the protective role of bacterial communities. Ongoing work to understand skin microbiome resilience after pathogen disturbance has the potential to identify key taxa involved in disease resistance.