Amphibian fungal panzootic causes catastrophic and ongoing loss of biodiversity

Gut pathobiome mediates behavioral and developmental disorders in biotoxin-exposed amphibians

Emerging evidence suggests a link between alterations in the gut microbiome and adverse health outcomes in the hosts exposed to environmental pollutants. Yet, the causal relationships and underlying mechanisms remain largely undefined. Here we show that exposure to biotoxins can affect gut pathobiome assembly in amphibians, which in turn triggers the toxicity of exogenous pollutants. We used Xenopus laevis as a model in this study. Tadpoles exposed to tropolone demonstrated notable developmental impairments and increased locomotor activity, with a reduction in total length by 4.37%-22.48% and an increase in swimming speed by 49.96%-84.83%. Fusobacterium and Cetobacterium are predominant taxa in the gut pathobiome of tropolone-exposed tadpoles. The tropolone-induced developmental and behavioral disorders in the host were mediated by assembly of the gut pathobiome, leading to transcriptome reprogramming. This study not only advances our understanding of the intricate interactions between environmental pollutants, the gut pathobiome, and host health but also emphasizes the potential of the gut pathobiome in mediating the toxicological effects of environmental contaminants.

Broad-scale pesticide screening finds anticoagulant rodenticide and legacy pesticides in Australian frogs

Pesticide contamination poses a significant threat to non-target wildlife, including amphibians, many of which are already highly threatened. This study assessed the extent of pesticide exposure in dead frogs collected during a mass mortality event across eastern New South Wales, Australia between July 2021 and March 2022. Liver tissue from 77 individual frogs of six species were analysed for >600 legacy and contemporary pesticides, including rodenticides. More than a third (36 %) of the liver samples contained at least one of the following pesticides: brodifacoum, dieldrin, DDE, heptachlor/heptachlor epoxide, fipronil sulfone, and 2-methyl-4-chlorophenoxyacetic acid (MCPA). Brodifacoum, a second-generation anticoagulant rodenticide, was found in four of the six frog species analysed: the eastern banjo frog (Limnodynastes dumerilii), cane toad (Rhinella marina), green tree frog (Litoria caerulea) and Peron's tree frog (Litoria peronii). This is the first report of anticoagulant rodenticide detected in wild amphibians, raising concerns about potential impacts on frogs and extending the list of taxa shown to accumulate rodenticides. Dieldrin, a banned legacy pesticide, was also detected in two species: striped marsh frog (Limnodynastes peronii) and green tree frog (Litoria caerulea). The toxicological effects of these pesticides on frogs are difficult to infer due to limited comparable studies; however, due to the low frequency of detection the presence of these pesticides was not considered a major contributing factor to the mass mortality event. Additional research is needed to investigate the effects of pesticide exposure on amphibians, particularly regarding the impacts of second-generation anticoagulant rodenticides. There is also need for continued monitoring and improved conservation management strategies for the mitigation of the potential threat of pesticide exposure and accumulation in amphibian populations.

The impacts of water quality on the amphibian chytrid fungal pathogen: A systematic review

The pathogenic fungus Batrachochytrium dendrobatidis has caused declines of amphibians worldwide. Yet our understanding of how water quality influences fungal pathogenicity is limited. Here, we reviewed experimental studies on the effect of water quality on this pathogen to determine which parameters impacted disease dynamics consistently. The strongest evidence for protective effects is salinity which shows strong antifungal properties in hosts at natural levels. Although many fungicides had detrimental effects on the fungal pathogen in vitro, their impact on the host is variable and they can worsen infection outcomes. However, one fungicide, epoxiconazole, reduced disease effects experimentally and likely in the field. While heavy metals are frequently studied, there is weak evidence that they influence infection outcomes. Nitrogen and phosphorous do not appear to impact pathogen growth or infection in the amphibian host. The effects of other chemicals, like pesticides and disinfectants on infection were mostly unclear with mixed results or lacking an in vivo component. Our study shows that water chemistry does impact disease dynamics, but the effects of specific parameters require more investigation. Improving our understanding of how water chemistry influences disease dynamics will help predict the impact of chytridiomycosis, especially in amphibian populations affected by land use changes.

Quantifying intraspecific variation in host resistance and tolerance to a lethal pathogen

Testing for intraspecific variation for host tolerance or resistance in wild populations is important for informing conservation decisions about captive breeding, translocation, and disease treatment. Here, we test the importance of tolerance and resistance in multiple populations of boreal toads (Anaxyrus boreas boreas) against Batrachochytrium dendrobatidis (Bd), the amphibian fungal pathogen responsible for the greatest host biodiversity loss due to disease. Boreal toads have severely declined in Colorado (CO) due to Bd, but toad populations challenged with Bd in western Wyoming (WY) appear to be less affected. We used a common garden infection experiment to expose post-metamorphic toads sourced from four populations (2 in CO and 2 in WY) to Bd and monitored changes in mass, pathogen burden and survival for 8 weeks. We used a multi-state modelling approach to estimate weekly survival and transition probabilities between infected and cleared states, reflecting a dynamic infection process that traditional approaches fail to capture. We found that WY boreal toads are more tolerant to Bd infection with higher survival probabilities than those in CO when infected with identical pathogen burdens. WY toads also appeared more resistant to Bd with a higher probability of infection clearance and an average of 5 days longer to reach peak infection burdens. Our results demonstrate strong intraspecific differences in tolerance and resistance that likely contribute to why population declines vary regionally across this species. Our multi-state framework allowed us to gain inference on typically hidden disease processes when testing for host tolerance or resistance. Our findings demonstrate that describing an entire host species as 'tolerant' or 'resistant' (or lack thereof) is unwise without testing for intraspecific variation.

Alpine salamanders at risk? The current status of an emerging fungal pathogen

Amphibians globally suffer from emerging infectious diseases like chytridiomycosis caused by the continuously spreading chytrid fungi. One is Batrachochytrium salamandrivorans (Bsal) and its disease ‒ the 'salamander plague' ‒ which is lethal to several caudate taxa. Recently introduced into Western Europe, long distance dispersal of Bsal, likely through human mediation, has been reported. Herein we study if Alpine salamanders (Salamandra atra and S. lanzai) are yet affected by the salamander plague in the wild. Members of the genus Salamandra are highly susceptible to Bsal leading to the lethal disease. Moreover, ecological modelling has shown that the Alps and Dinarides, where Alpine salamanders occur, are generally suitable for Bsal. We analysed skin swabs of 818 individuals of Alpine salamanders and syntopic amphibians at 40 sites between 2017 to 2022. Further, we compiled those with published data from 319 individuals from 13 sites concluding that Bsal infections were not detected. Our results suggest that the salamander plague so far is absent from the geographic ranges of Alpine salamanders. That means that there is still a chance to timely implement surveillance strategies. Among others, we recommend prevention measures, citizen science approaches, and ex situ conservation breeding of endemic salamandrid lineages.

Are novel or locally adapted pathogens more devastating and why? Resolving opposing hypotheses

There is a rich literature highlighting that pathogens are generally better adapted to infect local than novel hosts, and a separate seemingly contradictory literature indicating that novel pathogens pose the greatest threat to biodiversity and public health. Here, using Batrachochytrium dendrobatidis, the fungus associated with worldwide amphibian declines, we test the hypothesis that there is enough variance in "novel" (quantified by geographic and phylogenetic distance) host-pathogen outcomes to pose substantial risk of pathogen introductions despite local adaptation being common. Our continental-scale common garden experiment and global-scale meta-analysis demonstrate that local amphibian-fungal interactions result in higher pathogen prevalence, pathogen growth, and host mortality, but novel interactions led to variable consequences with especially virulent host-pathogen combinations still occurring. Thus, while most pathogen introductions are benign, enough variance exists in novel host-pathogen outcomes that moving organisms around the planet greatly increases the chance of pathogen introductions causing profound harm.

Rapid Evolution of Resistance and Tolerance Leads to Variable Host Recoveries following Disease-Induced Declines

AbstractRecoveries of populations that have suffered severe disease-induced declines are being observed across disparate taxa. Yet we lack theoretical understanding of the drivers and dynamics of recovery in host populations and communities impacted by infectious disease. Motivated by disease-induced declines and nascent recoveries in amphibians, we developed a model to ask the following question: How does the rapid evolution of different host defense strategies affect the transient recovery trajectories of hosts following pathogen invasion and disease-induced declines? We found that while host life history is predictably a major driver of variability in population recovery trajectories (including declines and recoveries), populations that use different host defense strategies (i.e., tolerance, avoidance resistance, and intensity-reduction resistance) experience notably different recoveries. In single-species host populations, populations evolving tolerance recovered on average four times slower than populations evolving resistance. Moreover, while populations using avoidance resistance strategies had the fastest potential recovery rates, these populations could get trapped in long transient states at low abundance prior to recovery. In contrast, the recovery of populations evolving intensity-reduction resistance strategies were more consistent across ecological contexts. Overall, host defense strategies strongly affect the transient dynamics of population recovery and may affect the ultimate fate of real populations recovering from disease-induced declines.

Immune and endocrine alterations at the early stage of inflammatory assemblage in toads after stimulation with heat-killed bacteria (Aeromonas hydrophila)

The red-leg syndrome in amphibians is a condition commonly associated with the bacteria Aeromonas hydrophila and has led to population declines. However, there is little information concerning the inflammatory assemblage in infected anurans. We evaluated immune and endocrine alterations induced by stimulation with heat-killed A. hydrophila injected in Rhinella diptycha toads. Control animals were not manipulated, while the others were separated into groups that received intraperitoneal injection of 300 μl of saline or heat-killed bacteria: groups A1 (3 × 107 cells), A2 (3 × 108 cells), and A3 (3 × 109 cells). Animals were bled and euthanized six hours post-injection. We evaluated neutrophil: lymphocyte ratio (NLR), plasma bacterial killing ability (BKA), testosterone (T), melatonin (MEL), and corticosterone (CORT) plasma levels. Heat-killed A. hydrophila increased CORT and NLR, and decreased MEL, especially at higher concentrations. There was no effect of treatment on T and BKA. We then selected the saline and A3 groups to conduct mRNA expression of several genes including glucocorticoid receptor (GR), toll-like receptor-4 (TLR-4), interferon-γ (IFN-γ), interleukin (IL)-1β, IL-6, and IL-10. We found higher expression of IL-6, IL-1β, IL-10, and IFN-γ in group A3 compared to the saline group. These results indicate the beginning of an inflammatory assemblage, notably at the two highest concentrations of bacteria, and give a better understanding of how anurans respond to an infection within an integrated perspective, evaluating different physiological aspects. Future studies should investigate later phases of the immune response to elucidate more about the inflammation in amphibians challenged with A. hydrophila.

Enhancing amphibian biomonitoring through eDNA metabarcoding

Surveying biodiversity has taken a quantum leap with environmental DNA (eDNA) metabarcoding, an immensely powerful approach lauded for its efficiency, sensitivity, and non-invasiveness. This approach emerges as a game-changer for the elusive realm of endangered and rare species-think nocturnal, environmentally elusive amphibians. Here, we have established a framework for constructing a reliable metabarcoding pipeline for amphibians, covering primer design, performance evaluation, laboratory validation, and field validation processes. The Am250 primer, located on the mitochondrial 16S gene, was optimal for the eDNA monitoring of amphibians, which demonstrated higher taxonomic resolution, smaller species amplification bias, and more extraordinary detection ability compared to the other primers tested. Am250 primer exhibit an 83.8% species amplification rate and 75.4% accurate species identification rate for Chinese amphibians in the in silico PCR and successfully amplified all tested species of the standard samples in the in vitro assay. Furthermore, the field-based mesocosm experiment showed that DNA can still be detected by metabarcoding even days to weeks after organisms have been removed from the mesocosm. Moreover, field mesocosm findings indicate that eDNA metabarcoding primers exhibit different read abundances, which can affect the relative biomass of species. Thus, appropriate primers should be screened and evaluated by three experimental approaches: in silico PCR simulation, target DNA amplification, and mesocosm eDNA validation. The selection of a single primer set or multiple primers' combination should be based on the monitoring groups to improve the species detection rate and the credibility of results.

Invasibility of a North American soil ecosystem to amphibian-killing fungal pathogens

North American salamanders are threatened by intercontinental spread of chytridiomycosis, a deadly disease caused by the fungal pathogen Batrachochytrium salamandrivorans (Bsal). To predict potential dispersal of Bsal spores to salamander habitats, we evaluated the capacity of soil microbial communities to resist invasion. We determined the degree of habitat invasibility using soils from five locations throughout the Great Smoky Mountains National Park, a region with a high abundance of susceptible hosts. Our experimental design consisted of replicate soil microcosms exposed to different propagule pressures of the non-native pathogen, Bsal, and an introduced but endemic pathogen, B. dendrobatidis (Bd). To compare growth and competitive interactions, we used quantitative PCR, live/dead cell viability assays, and full-length 16S rRNA sequencing. We found that soil microcosms with intact bacterial communities inhibited both Bsal and Bd growth, but inhibitory capacity diminished with increased propagule pressure. Bsal showed greater persistence than Bd. Linear discriminant analysis (LDA) identified the family Burkolderiaceae as increasing in relative abundance with the decline of both pathogens. Although our findings provide evidence of environmental filtering in soils, such barriers weakened in response to pathogen type and propagule pressure, showing that habitats vary their invasibility based on properties of their local microbial communities.

Medicinal Chemistry Gone Wild

Earth is currently experiencing a mass extinction event. The flora and fauna of our planet are experiencing mass die-offs from a multitude of factors, with wildlife disease emerging as one parameter where medicinal chemists are equipped to intervene. While contemporary medicinal chemistry focuses on human health, many traditional pharmaceutical companies have historic roots in human health, animal health, and plant health. This trifecta of health sciences perfectly maps to the current field of One Health, which recognizes that optimal health outcomes can only be achieved through the health of humans, plants, animals, and their shared environments. This Perspective imagines a world where state-of-the-art medicinal chemistry tactics are used to prevent the extinction of endangered species and points to preliminary work in the emerging area of conservation medicine.

Chytridiomycosis causes high amphibian mortality prior to the completion of metamorphosis

Amphibian populations are undergoing extensive declines globally. The fungal disease chytridiomycosis, caused by the pathogenic fungus Batrachochytrium dendrobatidis (Bd), is a primary contributor to these declines. The amphibian metamorphic stages (Gosner stages 42-46) are particularly vulnerable to a range of stressors, including Bd. Despite this, studies that explicitly examine host response to chytridiomycosis throughout the metamorphic stages are lacking. We aimed to determine how Bd exposure during the larval stages impacts metamorphic development and infection progression in the endangered Fleay's barred frog (Mixophyes fleayi). We exposed M. fleayi to Bd during pro-metamorphosis (Gosner stages 35-38) and monitored infection dynamics throughout metamorphosis. We took weekly morphological measurements (weight, total body length, snout-vent-length and Gosner stage) and quantified Bd load using qPCR. While we observed minimal impact of Bd infection on animal growth and development, Bd load varied throughout ontogeny, with an infection load plateau during the tadpole stages (Gosner stages 35-41) and temporary infection clearance at Gosner stage 42. Bd load increased exponentially between Gosner stages 42 and 45, with most exposed animals becoming moribund at Gosner stage 45, prior to the completion of metamorphosis. There was variability in infection outcome of exposed individuals, with a subgroup of animals (n = 5/29) apparently clearing their infection while the majority (n = 21/29) became moribund with high infection burdens. This study demonstrates the role that metamorphic restructuring plays in shaping Bd infection dynamics and raises the concern that substantial Bd-associated mortality could be overlooked in the field due to the often cryptic nature of these latter metamorphic stages. We recommend future studies that directly examine the host immune response to Bd infection throughout metamorphosis, incorporating histological and molecular methods to elucidate the mechanisms responsible for the observed trends.

Microplastics as an emerging threat to amphibians: Current status and future perspectives

Given their pervasiveness in the environment, particularly in aquatic ecosystems, plastics are posing a growing concern worldwide. Many vertebrates and invertebrates in marine, freshwater, and terrestrial ecosystems exhibit microplastic (MP) uptake and accumulation. Some studies have indicated the fatal impacts of MPs on animals and their possible transfer through food chains. Thus, it is crucial to study MP pollution and its impacts on environment-sensitive and globally threatened animal groups, such as amphibians, which also play an important role in the energy transfer between ecosystems. Unfortunately, research in this field is lacking and sources of organized information are also scarce. Hence, we systematically reviewed published literature on MPs in amphibians to fill the existing knowledge gap. Our review revealed that most of the previous studies have focused on MP bioaccumulation in amphibians, whereas, only a few research highlighted its impacts. We found that more than 80% of the studied species exhibited MP accumulation. MPs were reported to persist in different organs for a long time and get transferred to other trophic levels. They can also exhibit cytotoxic and mutagenic effects and may have fatal impacts. Moreover, they can increase the disease susceptibility of amphibians. Our study concludes the MPs as a potential threat to amphibians and urges increasing the scope and frequency of research on MP pollution and its impacts on this vulnerable animal group. We also provide a generalized method for studying MPs in amphibians with future perspectives and research directions. Our study is significant for extending the knowledge of MPs and their impacts on amphibians and guiding prospective research.

Bacteria from the skin of amphibians promote growth of Arabidopsis thaliana and Solanum lycopersicum by modifying hormone-related transcriptome response

Plants and microorganisms establish beneficial associations that can improve their development and growth. Recently, it has been demonstrated that bacteria isolated from the skin of amphibians can contribute to plant growth and defense. However, the molecular mechanisms involved in the beneficial effect for the host are still unclear. In this work, we explored whether bacteria isolated from three tropical frogs species can contribute to plant growth. After a wide screening, we identified three bacterial strains with high biostimulant potential, capable of modifying the root structure of Arabidopsis thaliana plants. In addition, applying individual bacterial cultures to Solanum lycopersicum plants induced an increase in their growth. To understand the effect that these microorganisms have over the host plant, we analysed the transcriptomic profile of A. thaliana during the interaction with the C32I bacterium, demonstrating that the presence of the bacteria elicits a transcriptional response associated to plant hormone biosynthesis. Our results show that amphibian skin bacteria can function as biostimulants to improve agricultural crops growth and development by modifying the plant transcriptomic responses.

An endogenous DNA virus in an amphibian-killing fungus associated with pathogen genotype and virulence

The global panzootic lineage (GPL) of the pathogenic fungus Batrachochytrium dendrobatidis (Bd) has caused severe amphibian population declines, yet the drivers underlying the high frequency of GPL in regions of amphibian decline are unclear. Using publicly available Bd genome sequences, we identified multiple non-GPL Bd isolates that contain a circular Rep-encoding single-stranded (CRESS)-like DNA virus, which we named Bd DNA virus 1 (BdDV-1). We further sequenced and constructed genome assemblies with long read sequences to find that the virus is integrated into the nuclear genome in some strains. Attempts to cure virus-positive isolates were unsuccessful; however, phenotypic differences between naturally virus-positive and virus-negative Bd isolates suggested that BdDV-1 decreases the growth of its host in vitro but increases the virulence of its host in vivo. BdDV-1 is the first-described CRESS DNA mycovirus of zoosporic true fungi, with a distribution inversely associated with the emergence of the panzootic lineage.

Morbidity in California giant salamander (Dicamptodon ensatus Eschscholtz, 1833) caused by Euryhelmis sp. Poche, 1926 (Trematoda: Heterophyiidae)

In the fall of 2021, California Department of Fish and Wildlife reported larval and adult California giant salamanders (Dicamptodon ensatus Eschscholtz, 1833) with skin lesions at multiple creeks in Santa Clara and Santa Cruz Counties, California, USA. Field signs in both stages included rough, lumpy textured skin, and larvae with tails that were disproportionately long, flat, wavy, and flaccid. Presence of large-bodied larvae suggested delayed metamorphosis, with some larvae having cloudy eyes and suspected blindness. To determine the cause of the disease, three first-of-the-year salamanders from one location were collected, euthanized with 20% benzocaine, and submitted for necropsy to the U.S. Geological Survey, National Wildlife Health Center. Upon gross examination, all salamanders were emaciated with no internal fat stores, and had multiple pinpoint to 1.5-mm diameter raised nodules in the skin over the body, including the head, gills, dorsum, ventrum, all four limbs, and the tail; one also had nodules in the oral cavity and tongue. Histologically all salamanders had multiple encysted metacercariae in the dermis, subcutis, and skeletal muscles of the head, body, and tail that were often associated with granulomatous and granulocytic inflammation and edema. A small number of encysted metacercariae or empty cysts were present in the gills with minimal inflammation, and rarely in the kidney with no associated inflammation. Morphology of live metacercariae (Trematoda: Heterophyiidae), and sequencing of the 28S rRNA gene identified a species of Euryhelmis (Poche, 1926). Artificial digestion of a 1.65 g, decapitated, eviscerated carcass yielded 773 metacercariae, all of similar size and morphology as the live specimens. Based on these findings, the poor body condition of these salamanders was concluded to be due to heavy parasite burden. Environmental factors such as drought, increased temperature, and overcrowded conditions may be exacerbating parasite infections in these populations of salamander.

Coevolution of a generalist pathogen with many hosts: the case of the amphibian chytrid Batrachochytrium dendrobatidis

Generalist pathogens maintain infectivity in numerous hosts; how this broad ecological niche impacts host-pathogen coevolution remains to be widely explored. Batrachochytrium dendrobatidis (Bd) is a highly generalist pathogenic fungus that has caused devastating declines in hundreds of amphibian species worldwide. This review examines amphibian chytridiomycosis host-pathogen interactions and available evidence for coevolution between Bd and its numerous hosts. We summarize recent evidence showing that Bd genotypes vary in geographic distribution and virulence, and that amphibian species also vary in Bd susceptibility according to their geographic distribution. How much variation can be explained by phenotypic plasticity or genetic differences remains uncertain. Recent research suggests that Bd genotypes display preferences for specific hosts and that some hosts are undergoing evolution as populations rebound from Bd outbreaks. Taken together, these findings suggest the potential for coevolution to occur and illuminate a path for addressing open questions through integrating historical and contemporary genetic data.

Field-based molecular detection of Batrachochytrium dendrobatidis in critically endangered Atelopus toads and aquatic habitats in Ecuador

Batrachochytrium dendrobatidis (Bd) is a lethal fungal species that parasitizes vertebrates and is associated with the worldwide decline of amphibian populations. The development of sensitive, rapid detection methods, particularly DNA-based techniques, is critical for effective management strategies. This study evaluates the efficacy of DNA extraction and a portable PCR device in a mountable field laboratory setup for detecting Bd near the habitats of three critically endangered Atelopus toad species in Ecuador. We collected skin swabs from Atelopus balios, A. nanay, and A. bomolochos, and environmental DNA (eDNA) samples from streams in Andean and coastal regions of Ecuador. For eDNA, a comparison was made with duplicates of the samples that were processed in the field and in a standard university laboratory. Our findings revealed Bd detection in eDNA and swabs from 6 of 12 water samples and 10 of 12 amphibian swab samples. The eDNA results obtained in the field laboratory were concordant with those obtained under campus laboratory conditions. These findings highlight the potential of field DNA-based monitoring techniques for detecting Bd in amphibian populations and their aquatic habitats, particularly in remote areas. Furthermore, this research aligns with the National Action Plan for the Conservation of Ecuadorian Amphibians and contributes to the global effort to control this invasive and deadly fungus.

Chemical disinfection as a simple and reliable method to control the amphibian chytrid fungus at breeding points of endangered amphibians

Chytridiomycosis caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd) is pushing amphibians towards extinction. Whilst mitigation methods were suggested a decade ago, we lack field trials testing their efficacy. We used the agrochemical fungicide, tebuconazole, to treat Bd infected breeding waterbodies of an endangered species that is highly susceptible to the fungus. Just two applications of tebuconazole led to a significant reduction in infection loads in the vast majority of sites, and at six sites the disinfection remained one/two-years post-application. Tebuconazole values drastically decreased in the waterbodies within a week after application, with no significant effects on their hydrochemical and hydrobiological characteristics. Although the use of chemicals in natural populations is undesirable, the growing existential threat to amphibians all over the world indicates that effective interventions in selected populations of endangered species are urgently needed.

Evolutionary analysis of major histocompatibility complex variants in chytrid-resistant and susceptible amphibians

An amphibian emerging infectious disease (EID), chytridiomycosis, caused by Batrachochytrium dendrobatidis (Bd), originated in Asia but primarily led to declines and extinctions in amphibian populations outside of Asia. Host major histocompatibility complex (MHC) molecules exhibit high polymorphism, and the evolution of MHC can be influenced by recombination and pathogens. Previous studies have indicated that host MHC class II is associated with Bd resistance. In this study, I conducted recombination and selection tests on functional MHC IIß1 alleles from an Asian Bd-resistant anuran species (Bufo gargarizans) and an Australasian Bd-susceptible species (Litoria caerulea). Recombination at the same site was identified in both species, supporting the hypothesis that recombination contributes to MHC IIß1 diversity in amphibians. Positive selection was observed in MHC IIß1 alleles in both species. In L. caerulea, at least four amino acid sites were identified under significant positive selection in the MHC IIß1, whereas these sites were either negatively selected or conserved in B. gargarizans. This suggests these sites might be selected for Bd resistance. Hydrophobicity was detected in certain amino acid sites relating to Bd resistance, suggesting this physicochemical property may be a factor selected to counteract Bd infection. These findings of this study provide an evolutionary basis for understanding how amphibian MHC IIß1 may undergo selection in response to chytrid infection.

From genomic threat assessment to conservation action

Genomic approaches have the potential to play a pivotal role in conservation, both to detect threats to species and populations and to restore biodiversity through actions. We here separate these approaches into two subdisciplines, vulnerability and restoration genomics, and discuss current applications, outstanding questions, and future potential.

Impact of Climate Change on the Distribution of Three Rare Salamanders (Liua shihi, Pseudohynobius jinfo, and Tylototriton wenxianensis) in Chongqing, China, and Their Conservation Implications

The Wushan Salamander (Liua shihi), Jinfo Salamander (Pseudohynobius jinfo), and Wenxian Knobby Salamander (Tylototriton wenxianensis) are rare national Class II protected wild animals in China. We performed MaxEnt modeling to predict and analyze the potential distribution and trends of these species in Chongqing under current and future climate conditions. Species distribution data were primarily obtained from field surveys, supplemented by museum collections and the existing literature. These efforts yielded 636 records, including 43 for P. jinfo, 23 for T. wenxianensis, and 570 for L. shihi. Duplicate records within the same 100 m × 100 m grid cell were removed using ENMTools, resulting in 10, 12, and 58 valid distribution points for P. jinfo, T. wenxianensis, and L. shihi, respectively. The optimization of feature class parameters (FC) and the regularization multiplier (RM) were applied using R package "ENMeval 2.0" to establish the optimal model with MaxEnt. The refined models were applied to simulate the suitable distribution areas for the three species. The results indicate that the current suitable habitat area for L. shihi accounted for 9.72% of the whole region of the Chongqing municipality. It is projected that, by 2050, the proportion of suitable habitat will increase to 12.54% but will decrease to 11.98% by 2070 and further decline to 8.80% by 2090. The current suitable habitat area for P. jinfo accounted for 1.08% of the whole region of the Chongqing municipality, which is expected to decrease to 0.31%% by 2050, 0.20% by 2070, and 0.07% by 2090. The current suitable habitat area for T. wenxianensis accounted for 0.81% of the whole region of the Chongqing municipality, which is anticipated to decrease to 0.37% by 2050, 0.21% by 2070, and 0.06% by 2090. Human disturbance, climate variables, and habitat characteristics are the primary factors influencing the distribution of three salamander species in Chongqing. The proximity to roads significantly impacts L. shihi, while climate conditions mainly affect P. jinfo, and the distance to water sources is crucial for T. wenxianensis. The following suggestions were made based on key variables identified for each species: (1) For L. shihi, it is imperative to minimize human disturbances and preserve areas without roads and the existing vegetation within nature reserves to ensure their continued existence. (2) For P. jinfo, the conservation of high-altitude habitats is of utmost importance, along with the reduction in disturbances caused by roads to maintain the species' ecological niche. (3) For T. wenxianensis, the protection of aquatic habitats is crucial. Additionally, efforts to mitigate the impacts of road construction and enhance public awareness are essential for the preservation of this species and the connectivity of its habitats.

Advances in Managing Chytridiomycosis for Australian Frogs: Gradarius Firmus Victoria

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

Using eDNA to survey amphibians: Methods, applications, and challenges

In recent years, environmental DNA (eDNA) has received attention from biologists due to its sensitivity, convenience, labor and material efficiency, and lack of damage to organisms. The extensive application of eDNA has opened avenues for the monitoring and biodiversity assessment of amphibians, which are frequently small and difficult to observe in the field, in areas such as biodiversity survey assessment and detection of specific, rare and threatened, or alien invasive species. However, the accuracy of eDNA can be influenced by factors such as ambient temperature, pH, and false positives or false negatives, which makes eDNA an adjunctive tool rather than a replacement for traditional surveys. This review provides a concise overview of the eDNA method and its workflow, summarizes the differences between applying eDNA for detecting amphibians and other organisms, reviews the research progress in eDNA technology for amphibian monitoring, identifies factors influencing detection efficiency, and discusses the challenges and prospects of eDNA. It aims to serve as a reference for future research on the application of eDNA in amphibian detection.

Sequence capture identifies fastidious chytrid fungi directly from host tissue

The chytrid fungus Batrachochytrium dendrobatidis (Bd) was discovered in 1998 as the cause of chytridiomycosis, an emerging infectious disease causing mass declines in amphibian populations worldwide. The rapid population declines of the 1970s-1990s were likely caused by the spread of a highly virulent lineage belonging to the Bd-GPL clade that was introduced to naïve susceptible populations. Multiple genetically distinct and regional lineages of Bd have since been isolated and sequenced, greatly expanding the known biological diversity within this fungal pathogen. To date, most Bd research has been restricted to the limited number of samples that could be isolated using culturing techniques, potentially causing a selection bias for strains that can grow on media and missing other unculturable or fastidious strains that are also present on amphibians. We thus attempted to characterize potentially non-culturable genetic lineages of Bd from distinct amphibian taxa using sequence capture technology on DNA extracted from host tissue and swabs. We focused our efforts on host taxa from two different regions that likely harbored distinct Bd clades: (1) wild-caught leopard frogs (Rana) from North America, and (2) a Japanese Giant Salamander (Andrias japonicus) at the Smithsonian Institution's National Zoological Park that exhibited signs of disease and tested positive for Bd using qPCR, but multiple attempts failed to isolate and culture the strain for physiological and genetic characterization. We successfully enriched for and sequenced thousands of fungal genes from both host clades, and Bd load was positively associated with number of recovered Bd sequences. Phylogenetic reconstruction placed all the Rana-derived strains in the Bd-GPL clade. In contrast, the A. japonicus strain fell within the Bd-Asia3 clade, expanding the range of this clade and generating additional genomic data to confirm its placement. The retrieved ITS locus matched public barcoding data from wild A. japonicus and Bd infections found on other amphibians in India and China, suggesting that this uncultured clade is widespread across Asia. Our study underscores the importance of recognizing and characterizing the hidden diversity of fastidious strains in order to reconstruct the spatiotemporal and evolutionary history of Bd. The success of the sequence capture approach highlights the utility of directly sequencing pathogen DNA from host tissue to characterize cryptic diversity that is missed by culture-reliant approaches.

Environment predicts Batrachochytrium dendrobatidis lineage distribution and zones of recombination in South Africa

The amphibian-infecting chytrid fungus, Batrachochytrium dendrobatidis (Bd), is widespread throughout Africa and is linked to declines of populations and species across the continent. While it is well established that the lineage of Bd encodes traits which determine disease severity, knowledge around how lineages are distributed according to environmental envelope is unclear. We here studied the distribution of Bd in South Africa based on the two lineages found, BdGPL and BdCAPE, in terms of their genome and environmental envelope statistically associated with their distribution. We used Bd surveillance data from published studies, as well as data collected during fieldwork from across South Africa, Lesotho, and eSwatini with samples collected along a transect spanning most of South Africa from Lesotho to the west coast. We utilized lineage-typing qPCR to resolve the spatial distribution of BdGPL and BdCAPE across South Africa and used the resulting surveillance data to create a predictive ecological niche model for Bd lineages in South Africa. Phylogenomic analyses were performed on isolates sourced from across the transect. We show that BdGPL demonstrates a strong isolation by distance suggestive of stepping-stone dispersal, while BdCAPE showed two distinct clusters within their genomic structure that appear geographically and temporally clustered, indicating two separate invasions. Our predictive niche model revealed that the two lineages tended to occur in different ecotypes; BdGPL was associated with lower altitude, arid regions while BdCAPE occurred across cooler, higher altitude environs. Niche predictions identified a zone of lineage contact, where genomics identified inter-lineage recombinants. We argue that this zone of recombination should be prioritized for disease surveillance as it is a potential hotspot for the evolution of variants of amphibian chytrid with novel traits that may be epidemiologically relevant.

Genomic footprints of (pre) colonialism: Population declines in urban and forest túngara frogs coincident with historical human activity

Urbanisation is rapidly altering ecosystems, leading to profound biodiversity loss. To mitigate these effects, we need a better understanding of how urbanisation impacts dispersal and reproduction. Two contrasting population demographic models have been proposed that predict that urbanisation either promotes (facilitation model) or constrains (fragmentation model) gene flow and genetic diversity. Which of these models prevails likely depends on the strength of selection on specific phenotypic traits that influence dispersal, survival, or reproduction. Here, we a priori examined the genomic impact of urbanisation on the Neotropical túngara frog (Engystomops pustulosus), a species known to adapt its reproductive traits to urban selective pressures. Using whole-genome resequencing for multiple urban and forest populations we examined genomic diversity, population connectivity and demographic history. Contrary to both the fragmentation and facilitation models, urban populations did not exhibit substantial changes in genomic diversity or differentiation compared with forest populations, and genomic variation was best explained by geographic distance rather than environmental factors. Adopting an a posteriori approach, we additionally found both urban and forest populations to have undergone population declines. The timing of these declines appears to coincide with extensive human activity around the Panama Canal during the last few centuries rather than recent urbanisation. Our study highlights the long-lasting legacy of past anthropogenic disturbances in the genome and the importance of considering the historical context in urban evolution studies as anthropogenic effects may be extensive and impact nonurban areas on both recent and older timescales.

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.

Interplay between acetylation and ubiquitination of imitation switch chromatin remodeler Isw1 confers multidrug resistance in Cryptococcus neoformans

Cryptococcus neoformans poses a threat to human health, but anticryptococcal therapy is hampered by the emergence of drug resistance, whose underlying mechanisms remain poorly understood. Herein, we discovered that Isw1, an imitation switch chromatin remodeling ATPase, functions as a master modulator of genes responsible for in vivo and in vitro multidrug resistance in C. neoformans. Cells with the disrupted ISW1 gene exhibited profound resistance to multiple antifungal drugs. Mass spectrometry analysis revealed that Isw1 is both acetylated and ubiquitinated, suggesting that an interplay between these two modification events exists to govern Isw1 function. Mutagenesis studies of acetylation and ubiquitination sites revealed that the acetylation status of Isw1K97 coordinates with its ubiquitination processes at Isw1K113 and Isw1K441 through modulating the interaction between Isw1 and Cdc4, an E3 ligase. Additionally, clinical isolates of C. neoformans overexpressing the degradation-resistant ISW1K97Q allele showed impaired drug-resistant phenotypes. Collectively, our studies revealed a sophisticated acetylation-Isw1-ubiquitination regulation axis that controls multidrug resistance in C. neoformans.

The Binding, Infection, and Promoted Growth of Batrachochytrium dendrobatidis by the Ranavirus FV3

Increasing reports suggest the occurrence of co-infection between Ranaviruses such as Frog Virus 3 (FV3) and the chytrid fungus Batrachochytrium dendrobatidis (Bd) in various amphibian species. However, the potential direct interaction of these two pathogens has not been examined to date. In this study, we investigated whether FV3 can interact with Bd in vitro using qPCR, conventional microscopy, and immunofluorescent microscopy. Our results reveal the unexpected ability of FV3 to bind, promote aggregation, productively infect, and significantly increase Bd growth in vitro. To extend these results in vivo, we assessed the impact of FV3 on Xenopus tropicalis frogs previously infected with Bd. Consistent with in vitro results, FV3 exposure to previously Bd-infected X. tropicalis significantly increased Bd loads and decreased the host's survival.

Effects of grassland controlled burning on symbiotic skin microbes in Neotropical amphibians

Climate change has led to an alarming increase in the frequency and severity of wildfires worldwide. While it is known that amphibians have physiological characteristics that make them highly susceptible to fire, the specific impacts of wildfires on their symbiotic skin bacterial communities (i.e., bacteriomes) and infection by the deadly chytrid fungus, Batrachochytrium dendrobatidis, remain poorly understood. Here, we address this research gap by evaluating the effects of fire on the amphibian skin bacteriome and the subsequent risk of chytridiomycosis. We sampled the skin bacteriome of the Neotropical species Scinax squalirostris and Boana leptolineata in fire and control plots before and after experimental burnings. Fire was linked with a marked increase in bacteriome beta dispersion, a proxy for skin microbial dysbiosis, alongside a trend of increased pathogen loads. By shedding light on the effects of fire on amphibian skin bacteriomes, this study contributes to our broader understanding of the impacts of wildfires on vulnerable vertebrate species.

eDNA-based monitoring of Batrachochytrium dendrobatidis and Batrachochytrium salamandrivorans with ddPCR in Luxembourg ponds: taking signals below the Limit of Detection (LOD) into account

BACKGROUND

Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal) are two pathogenic fungi that are a significant threat to amphibian communities worldwide. European populations are strongly impacted and the monitoring of the presence and spread of these pathogens is crucial for efficient decision-making in conservation management.

RESULTS

Here we proposed an environmental DNA (eDNA) monitoring of these two pathogenic agents through droplet digital PCR (ddPCR) based on water samples from 24 ponds in Luxembourg. In addition, amphibians were swabbed in eight of the targeted ponds in order to compare the two approaches at site-level detection. This study allowed the development of a new method taking below-Limit of Detection (LOD) results into account thanks to the statistical comparison of the frequencies of false positives in no template controls (NTC) and below-LOD results in technical replicates. In the eDNA-based approach, the use of this method led to an increase in Bd and Bsal detection of 28 and 50% respectively. In swabbing, this resulted in 8% more positive results for Bd. In some samples, the use of technical replicates allowed to recover above-LOD signals and increase Bd detection by 35 and 33% respectively for eDNA and swabbing, and Bsal detection by 25% for eDNA.

CONCLUSIONS

These results confirmed the usefulness of technical replicates to overcome high levels of stochasticity in very low concentration samples even for a highly sensitive technique such as ddPCR. In addition, it showed that below-LOD signals could be consistently recovered and the corresponding amplification events assigned either to positive or negative detection via the method developed here. This methodology might be particularly worth pursuing in pathogenic agents' detection as false negatives could have important adverse consequences. In total, 15 ponds were found positive for Bd and four for Bsal. This study reports the first record of Bsal in Luxembourg.

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.

What factors influence the rediscovery of lost tetrapod species?

We created a database of lost and rediscovered tetrapod species, identified patterns in their distribution and factors influencing rediscovery. Tetrapod species are being lost at a faster rate than they are being rediscovered, due to slowing rates of rediscovery for amphibians, birds and mammals, and rapid rates of loss for reptiles. Finding lost species and preventing future losses should therefore be a conservation priority. By comparing the taxonomic and spatial distribution of lost and rediscovered tetrapod species, we have identified regions and taxa with many lost species in comparison to those that have been rediscovered-our results may help to prioritise search effort to find them. By identifying factors that influence rediscovery, we have improved our ability to broadly distinguish the types of species that are likely to be found from those that are not (because they are likely to be extinct). Some lost species, particularly those that are small and perceived to be uncharismatic, may have been neglected in terms of conservation effort, and other lost species may be hard to find due to their intrinsic characteristics and the characteristics of the environments they occupy (e.g. nocturnal species, fossorial species and species occupying habitats that are more difficult to survey such as wetlands). These lost species may genuinely await rediscovery. However, other lost species that possess characteristics associated with rediscovery (e.g. large species) and that are also associated with factors that negatively influence rediscovery (e.g. those occupying small islands) are more likely to be extinct. Our results may foster pragmatic search protocols that prioritise lost species likely to still exist.

Gene functions of the Ambystoma altamirani skin microbiome vary across space and time but potential antifungal genes are widespread and prevalent

Amphibian skin microbiomes can play a critical role in host survival against emerging diseases by protecting their host against pathogens. While a plethora of biotic and abiotic factors have been shown to influence the taxonomic diversity of amphibian skin microbiomes it remains unclear whether functional genomic diversity varies in response to temporal and environmental factors. Here we applied a metagenomic approach to evaluate whether seasonality, distinct elevations/sites, and pathogen presence influenced the functional genomic diversity of the A. altamirani skin microbiome. We obtained a gene catalogue of 92 107 nonredundant annotated genes and a set of 50 unique metagenome assembled genomes (MAGs). Our analysis showed that genes linked to general and potential antifungal traits significantly differed across seasons and sampling locations at different elevations. Moreover, we found that the functional genomic diversity of A. altamirani skin microbiome differed between B. dendrobatidis infected and not infected axolotls only during winter, suggesting an interaction between seasonality and pathogen infection. In addition, we identified the presence of genes and biosynthetic gene clusters (BGCs) linked to potential antifungal functions such as biofilm formation, quorum sensing, secretion systems, secondary metabolite biosynthesis, and chitin degradation. Interestingly genes linked to these potential antifungal traits were mainly identified in Burkholderiales and Chitinophagales MAGs. Overall, our results identified functional traits linked to potential antifungal functions in the A. altamirani skin microbiome regardless of variation in the functional diversity across seasons, elevations/sites, and pathogen presence. Our findings suggest that potential antifungal traits found in Burkholderiales and Chitinophagales taxa could be related to the capacity of A. altamirani to survive in the presence of Bd, although further experimental analyses are required to test this hypothesis.

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

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

Association between the skin microbiome and MHC class II diversity in an amphibian

Microbiomes play an important role in determining the ecology and behaviour of their hosts. However, questions remain pertaining to how host genetics shape microbiomes, and how microbiome composition influences host fitness. We explored the effects of geography, evolutionary history and host genetics on the skin microbiome diversity and structure in a widespread amphibian. More specifically, we examined the association between bacterial diversity and composition and the major histocompatibility complex class II exon 2 diversity in 12 moor frog (Rana arvalis) populations belonging to two geographical clusters that show signatures of past and ongoing differential selection. We found that while bacterial alpha diversity did not differ between the two clusters, MHC alleles/supertypes and genetic diversity varied considerably depending on geography and evolutionary history. Bacterial alpha diversity was positively correlated with expected MHC heterozygosity and negatively with MHC nucleotide diversity. Furthermore, bacterial community composition showed significant variation between the two geographical clusters and between specific MHC alleles/supertypes. Our findings emphasize the importance of historical demographic events on hologenomic variation and provide new insights into how immunogenetic host variability and microbial diversity may jointly influence host fitness with consequences for disease susceptibility and population persistence.

A new cryptic species of terrestrial breeding frog of the Pristimantisdanae Group (Anura, Strabomantidae) from montane forests in Ayacucho, Peru

Based on morphological and molecular characters, we describe a new species of terrestrial breeding frog of the Pristimantisdanae Group from montane forests of La Mar Province, Ayacucho Department in southern Peru, at elevations from 1200 to 2000 m a.s.l. The phylogenetic analysis, based on concatenated sequences of gene fragments of 16S rRNA, RAG1, COI and TYR suggests that the new species is a sister taxon of a clade that includes one undescribed species of Pristimantis from Cusco, Pristimantispharangobates and Pristimantisrhabdolaemus. The new species is most similar to P.rhabdolaemus, which differs by lacking scapular tubercules and by its smaller size (17.0-18.6 mm in males [n = 5], 20.8-25.2 mm in females [n = 5] in the new species vs. 22.8-26.3 mm in males [n = 19], 26.0-31.9 mm in females [n = 30] of P.rhabdolaemus). Additionally, we report the prevalence of Batrachochytriumdendrobatidis (Bd) in this species.

The global importance and interplay of colour-based protective and thermoregulatory functions in frogs

Small-scale studies have shown that colour lightness variation can have important physiological implications in ectotherms, with darker species having greater heating rates, as well as protection against pathogens and photooxidative damage. Using data for 41% (3059) of all known frog and toad species (Anura) from across the world, we reveal ubiquitous and strong clines of decreasing colour lightness towards colder regions and regions with higher pathogen pressure and UVB radiation. The relative importance of pathogen resistance is higher in the tropics and that of thermoregulation is higher in temperate regions. The results suggest that these functions influence colour lightness evolution in anurans and filtered for more similarly coloured species under climatic extremes, while their concurrent importance resulted in high within-assemblage variation in productive regions. Our findings indicate three important functions of colour lightness in anurans - thermoregulation, pathogen and UVB protection - and broaden support for colour lightness-environment relationships in ectotherms.

Glutathione is required for growth and cadmium tolerance in the amphibian chytrid fungus, Batrachochytrium dendrobatidis

Batrachochytrium dendrobatidis (Bd) is a lethal amphibian pathogen, partly due to its ability to evade the immune system of susceptible frog species. In many pathogenic fungi, the antioxidant glutathione is a virulence factor that helps neutralise oxidative stressors generated from host immune cells, as well as other environmental stressors such as heavy metals. The role of glutathione in stress tolerance in Bd has not been investigated. Here, we examine the changes in the glutathione pool after stress exposure and quantify the effect of glutathione depletion on cell growth and stress tolerance. Depletion of glutathione repressed growth and release of zoospores, suggesting that glutathione is essential for life cycle completion in Bd. Supplementation with <2 mM exogenous glutathione accelerated zoospore development, but concentrations >2 mM were strongly inhibitory to Bd cells. While hydrogen peroxide exposure lowered the total cellular glutathione levels by 42 %, glutathione depletion did not increase the sensitivity to hydrogen peroxide. Exposure to cadmium increased total cellular glutathione levels by 93 %. Glutathione-depleted cells were more sensitive to cadmium, and this effect was attenuated by glutathione supplementation, suggesting that glutathione plays an important role in cadmium tolerance. The effects of heat and salt were exacerbated by the addition of exogenous glutathione. The impact of glutathione levels on Bd stress sensitivity may help explain differences in host susceptibility to chytridiomycosis and may provide opportunities for synergistic therapeutics.

Linking microbiome and stress hormone responses in wild tropical treefrogs across continuous and fragmented forests

The amphibian skin microbiome is an important component of anti-pathogen defense, but the impact of environmental change on the link between microbiome composition and host stress remains unclear. In this study, we used radiotelemetry and host translocation to track microbiome composition and function, pathogen infection, and host stress over time across natural movement paths for the forest-associated treefrog, Boana faber. We found a negative correlation between cortisol levels and putative microbiome function for frogs translocated to forest fragments, indicating strong integration of host stress response and anti-pathogen potential of the microbiome. Additionally, we observed a capacity for resilience (resistance to structural change and functional loss) in the amphibian skin microbiome, with maintenance of putative pathogen-inhibitory function despite major temporal shifts in microbiome composition. Although microbiome community composition did not return to baseline during the study period, the rate of microbiome change indicated that forest fragmentation had more pronounced effects on microbiome composition than translocation alone. Our findings reveal associations between stress hormones and host microbiome defenses, with implications for resilience of amphibians and their associated microbes facing accelerated tropical deforestation.

Critical review of the phytohemagglutinin assay for assessing amphibian immunity

Infectious diseases are a major driver of the global amphibian decline. In addition, many factors, including genetics, stress, pollution, and climate change can influence the response to pathogens. Therefore, it is important to be able to evaluate amphibian immunity in the laboratory and in the field. The phytohemagglutinin (PHA) assay is an inexpensive and relatively non-invasive tool that has been used extensively to assess immunocompetence, especially in birds, and more recently in amphibians. However, there is substantial variation in experimental methodology among amphibian PHA studies in terms of species and life stages, PHA doses and injection sites, and use of experimental controls. Here, we compile and compare all known PHA studies in amphibians to identify knowledge gaps and develop best practices for future work. We found that research has only been conducted on a limited number of species, which may not reflect the diversity of amphibians. There is also a lack of validation studies in most species, so that doses and timing of PHA injection and subsequent swelling measurements may not effectively evaluate immunocompetence. Based on these and other findings, we put forward a set of recommendations to make future PHA studies more consistent and improve the ability to utilize this assay in wild populations, where immune surveillance is greatly needed.

Low-level pathogen infection and geographic location correlate with the skin microbiomes of Columbia spotted frogs (Rana luteiventris) in a montane landscape

The skin microbiome of amphibians can influence host susceptibility towards the fungal pathogen Batrachochytrium dendrobatidis (Bd), while simultaneously having the potential to be altered by Bd. Severe Bd infections are known to alter the amphibian skin microbiome; however, little is known about microbiome interactions in amphibians with low infection intensity. In addition to disease dynamics, environmental factors may influence the microbiome. To test for patterns in bacterial diversity based on pathogen infection and environmental factors, 399 Columbia spotted frogs (Rana luteiventris) were sampled throughout northern Idaho and northeastern Washington across two years. Bd prevalence and intensity were measured in 376 frogs, revealing a prevalence of 69%, but generally low infection intensity (Mean = 127 Bd zoospore equivalents among infected frogs). Skin bacterial communities were characterized in 92 frogs using 16S rRNA gene amplicon sequencing. Our results indicated correlations of decreasing Shannon diversity and evenness as infection intensity increased. Latitude was correlated with bacterial richness and Faith's Phylogenetic Diversity measures, indicating increased diversity in northern locations. Beta diversity (UniFrac) analyses revealed that skin microbiomes were distinct between infected and uninfected frogs, and infection intensity had a significant effect on microbiome composition. Site explained the majority of microbiome variation (weighted UniFrac: 57.5%), suggesting a combination of local habitat conditions explain variation, as only small proportions of variation could be explained by year, month, temperature, elevation, and latitude individually. Bacterial genera with potential for Bd-inhibitory properties were found with differential relative abundance in infected and uninfected frogs, with higher Stenotrophomonas and lower Pseudomonas relative abundance observed in infected frogs. Further study may indicate if Bd inhibition by members of the skin microbiome is an influence behind the low infection intensities observed and whether low Bd infection intensities are capable of altering skin microbiome composition.

Toxicities of Legacy and Current-Use PFAS in an Anuran: Do Larval Exposures Influence Responses to a Terrestrial Pathogen Challenge?

Legacy polyfluoroalkyl substances (PFAS) [perfluorooctanesulfonate (PFOS) and perfluorooctanoic acid (PFOA)] are being replaced by various other fluorinated compounds, such as hexafluoropropylene oxide dimer acid (GenX). These alternatives are thought to be less bioaccumulative and, therefore, less toxic than legacy PFAS. Contaminant exposures occur concurrently with exposure to natural stressors, including the fungal pathogen Batrachocytrium dendrobatidis (Bd). Despite evidence that other pollutants can increase the adverse effects of Bd on anurans, no studies have examined the interactive effects of Bd and PFAS. This study tested the growth and developmental effects of PFOS, PFOA, and GenX on gray treefrog (Hyla versicolor) tadpoles, followed by a Bd challenge after metamorphosis. Despite PFAS exposure only occurring during the larval stage, carry-over effects on growth were observed post metamorphosis. Further, PFAS interacted with Bd exposure to influence growth; Bd-exposed animals had significantly shorter SVL [snout-vent length (mm)] with significantly increased body condition, among other time-dependent effects. Our data suggest that larval exposure to PFAS can continue to impact growth in the juvenile stage after exposure has ended. Contrary to predictions, GenX affected terrestrial performance more consistently than its legacy congener, PFOA. Given the role of Bd in amphibian declines, further investigation of interactions of PFAS with Bd and other environmentally relevant pathogens is warranted.

Population-level call properties of endangered Dryophytes suweonensissensu lato (Anura: Amphibia) in South Korea

Calling is one of the unique amphibian characteristics that facilitates social communication and shows individuality; however, it also makes them vulnerable to predators. Researchers use amphibian call properties to study their population status, ecology, and behavior. This research scope has recently broadened to species identification and taxonomy. Dryophytes flaviventris has been separated from the endangered anuran species, D. suweonensis, based on small variations in genetic, morphometric, and temporal call properties observed in South Korea. The Chilgap Mountain (CM) was considered as the potential geographic barrier for the speciation. However, it initiated taxonomic debates as CM has been hardly used and is considered a potential barrier for other species. The calls of populations from both sides are also apparently similar. Thus, to verify the differences in call properties among populations of D. suweonensis sensu lato (s.l.; both of the species), we sampled and analyzed call data from five localities covering its distribution range, including the southern (S) and northern (N) parts of CM. We found significant differences in many call properties among populations; however, no specific pattern was observed. Some geographically close populations, such as Iksan (S), Wanju (S), and Gunsan (S), had significant differences, whereas many distant populations, such as Pyeongtaek (N) and Wanju (S), had no significant differences. Considering the goal of this study was only to observe the call properties, we cautiously conclude that the differences are at the population level rather than the species level. Our study indicates the necessity of further investigation into the specific status of D. flaviventris using robust integrated taxonomic approaches, including genetic and morphological parameters from a broader array of localities.

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.

Mr. Toad's Wild Fungi: Fungal Isolate Diversity on Colorado Boreal Toads and their Capacity for Pathogen Inhibition

The amphibian skin pathogen Batrachochytrium dendrobatidis (Bd) has caused an ongoing biodiversity crisis, including in the locally endangered Colorado boreal toad (Anaxyrus boreas boreas). Although researchers have investigated the bacteria living on amphibian skin and how they interact with Bd, there is less information about fungal community members. This study describes (1) the diversity of culturable fungi from boreal toad skin, (2) which subset of these isolates is Bd-inhibitory, and (3) how Bd affects these isolates' growth and morphology. Most isolates were from the orders Capnodiales, Helotiales, and Pleosporales. Of 16 isolates tested for Bd-inhibition, two from the genus Neobulgaria and three from Pseudeurotium inhibited Bd. Fungal growth in co-culture with Bd varied with weak statistical support for Neobulgaria sp. (isolate BTF_36) and cf Psychrophila (isolate BTF_60) (p-values = 0.076 and 0.092, respectively). Fungal morphology remained unchanged in co-culture with Bd, however, these results could be attributed to low replication per isolate. Nonetheless, two fungal isolates' growth may have been affected by Bd, implying that fungal growth changes in Bd co-culture could be a variable worth measuring in the future (with higher replication). These findings add to the sparse but growing literature on amphibian-associated fungi and suggest further study may uncover the relevance of fungi to amphibian health and Bd infection.

Amphibian Hymenochirus boettgeri as an experimental model for infection studies with the chytrid fungus Batrachochytrium dendrobatidis

Model organisms are crucial in research as they can provide key insights applicable to other species. This study proposes the use of the amphibian species Hymenochirus boettgeri, widely available through the aquarium trade, as a model organism for the study of chytridiomycosis, a disease caused by the fungus Batrachochytrium dendrobatidis (Bd) and linked to amphibian decline and extinction globally. Currently, no model organisms are used in the study of chytridiomycosis, particularly because of the lack of availability and nonstandardized methods. Thus, laboratories around the world use wild local species to conduct Bd infection experiments, which prevents comparisons between studies and reduces reproducibility. Here, we performed a series of Bd infection assays that showed that H. boettgeri has a dose- and genotype-dependent response, can generalize previous findings on virulence estimates in other species, and can generate reproducible results in replicated experimental conditions. We also provided valuable information regarding H. boettgeri husbandry, including care, housing, reproduction, and heat treatment to eliminate previous Bd infections. Together, our results indicate that H. boettgeri is a powerful and low-ecological-impact system for studying Bd pathogenicity and virulence.

Batrachochytrium dendrobatidis in the Arid and Thermally Extreme Sonoran Desert

Batrachochytrium dendrobatidis (Bd), the causative agent of the devastating global amphibian disease chytridiomycosis, was not projected to threaten amphibians in hot and arid regions due to its sensitivity to heat and desiccation. However, Bd is being detected more frequently than ever in hot and arid regions of Australia and the USA, challenging our current understanding of the environmental tolerances of the pathogen under natural conditions. We surveyed for Bd in an extremely hot and arid portion of the Sonoran Desert, where the pathogen is not projected to occur, and related presence and prevalence of the pathogen to local environmental conditions. We collected eDNA samples from isolated desert water sites including six tinajas and 13 catchments in June and August of 2020 and swabbed a total of 281 anurans of three species (red-spotted toad Anaxyrus punctatus, Couch's spadefoot Scaphiopus couchii, and the Sonoran Desert toad Incillius alvarius) across five catchments and six tinajas from June to September of 2020. Overall, Bd occurred at 68.4% of sites, despite extreme heat and aridity routinely exceeding tolerances established in laboratory studies. Average summer maximum air and water temperatures were 40.7°C and 30.7°C, respectively, and sites received an average of just 16.9 mm of precipitation throughout the summer monsoon season. Prevalence was low (5.7%) across species and life stage. Our results demonstrate that Bd is capable of persisting and infecting amphibians beyond its projected range, indicating a need to account for higher thermal tolerances when quantifying risk of Bd presence and infection.