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.

Priority effects alter microbiome composition and increase abundance of probiotic taxa in treefrog tadpoles.

Host-associated microbial communities, like other ecological communities, may be impacted by the colonization order of taxa through priority effects. Developing embryos and their associated microbiomes are subject to stochasticity during colonization by bacteria. For amphibian embryos, often developing externally in bacteria-rich environments, this stochasticity may be particularly impactful. For example, the amphibian microbiome can mitigate lethal outcomes from disease for their hosts; however, this may depend on microbiome composition. Here, we examined the assembly of the bacterial community in spring peeper (Pseudacris crucifer) embryos and tadpoles. First, we reared embryos from identified mating pairs in either lab or field environments to examine the relative impact of environment and parentage on embryo and tadpole bacterial communities. Second, we experimentally inoculated embryos to determine if priority effects (i) could be used to increase the relative abundance of Janthinobacterium lividum, an amphibian-associated bacteria capable of preventing fungal infection, and (ii) would lead to observed differences in the relative abundances of two closely related bacteria from the genus Pseudomonas. Using 16S rRNA gene amplicon sequencing, we observed differences in community composition based on rearing location and parentage in embryos and tadpoles. In the inoculation experiment, we found that priority inoculation could increase the relative abundance of J. lividum, but did not find that either Pseudomonas isolate was able to prevent colonization by the other when given priority. These results highlight the importance of environmental source pools and parentage in determining microbiome composition, while also providing novel methods for the administration of a known amphibian probiotic. IMPORTANCE: Harnessing the functions of host-associated bacteria is a promising mechanism for managing disease outcomes across different host species. In the case of amphibians, certain frog-associated bacteria can mitigate lethal outcomes of infection by the fungal pathogen Batrachochytrium dendrobatidis. Successful probiotic applications require knowledge of community assembly and an understanding of the ecological mechanisms that structure these symbiotic bacterial communities. In our study, we show the importance of environment and parentage in determining bacterial community composition and that community composition can be influenced by priority effects. Further, we provide support for the use of bacterial priority effects as a mechanism to increase the relative abundance of target probiotic taxa in a developing host. While our results show that priority effects are not universally effective across all host-associated bacteria, our ability to increase the relative abundance of specific probiotic taxa may enhance conservation strategies that rely on captive rearing of endangered vertebrates.

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.

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.

High-mortality epizootic Mycobacterium ulcerans ecovar Liflandii in a colony of Zaire Dwarf Clawed Frogs (Hymenochirus boettgeri).

Mycobacterium ulcerans ecovar Liflandii (MuLiflandii) was identified as the causative agent of mycobacteriosis in a research colony of Zaire dwarf clawed frogs (Hymenochirus boettgeri) at the University of Michigan. Clinical presentation included lethargy, generalized septicemia, cutaneous granulomas, coelomic effusion, and acute mortality. Identification of the mycobacterial species was based on molecular, microbiological, and histopathologic characteristics. These findings indicate that MuLiflandii is a primary cause of morbidity and mortality in Zaire dwarf clawed frogs and should be considered in the differential diagnosis of sepsis and coelomic effusion in amphibians. Mycobacterial speciation is important given the variability in pathogenesis within the family Mycobacteriaceae and the implications for both animal and human health as potential zoonoses. The Zaire dwarf clawed frog is a species common in the pet trade, and these findings provide consideration for this pathogen as a potentially important public health concern. This is the first report of MuLiflandii infection in the genus Hymenochirus and illustrates the diagnostic challenges of differentiating among both mycolactone-producing mycobacteria and Mycobacterium marinum. Furthermore, we demonstrate the utility of environmental sampling for this pathogen within the tank system, suggesting this mode of sampling could replace the need for direct frog surveillance.

Widespread occurrence of the amphibian chytrid panzootic lineage in Uruguay is constrained by climate.

The amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd) causes chytridiomycosis, a disease among the main causes of amphibian declines worldwide. However, Bd studies on Neotropical amphibians from temperate areas are scarce. We present a comprehensive survey of Bd in Uruguay, in temperate central eastern South America, carried out between 2006 and 2014. Skin swabs of 535 specimens of 21 native and exotic frogs were tested by PCR. We used individual-level data to examine the relationship between infection, climatic variables, and their effects on body condition and the number of prey items found in stomach contents. Infection was widespread in free-ranging anurans with an overall prevalence of 41.9%, detected in 15 native species, wild American bullfrogs Aquarana catesbeiana, and captive specimens of Ceratophrys ornata and Xenopus laevis. Three haplotypes of the Bd ITS region were identified in native amphibians, all belonging to the global panzootic lineage (BdGPL), of which only one was present in exotic hosts. Despite high infection frequencies in different anurans, we found no evidence of morbidity or mortality attributable to chytridiomycosis, and we observed no discernible impact on body condition or consumed prey. Climatic conditions at the time of our surveys suggested that the chance of infection is associated with monthly mean temperature, mean humidity, and total precipitation. Temperatures below 21°C combined with moderate humidity and pronounced rainfall may increase the likelihood of infection. Multiple haplotypes of BdGPL combined with high frequencies of infection suggest an enzootic pattern in native species, underscoring the need for continued monitoring.

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

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

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.

Limited impact of chytridiomycosis on juvenile frogs in a recovered species.

The amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd) has caused catastrophic frog declines on several continents, but disease outcome is mediated by a number of factors. Host life stage is an important consideration and many studies have highlighted the vulnerability of recently metamorphosed or juvenile frogs compared to adults. The majority of these studies have taken place in a laboratory setting, and there is a general paucity of longitudinal field studies investigating the influence of life stage on disease outcome. In this study, we assessed the effect of endemic Bd on juvenile Mixophyes fleayi (Fleay's barred frog) in subtropical eastern Australian rainforest. Using photographic mark-recapture, we made 386 captures of 116 individuals and investigated the effect of Bd infection intensity on the apparent mortality rates of frogs using a multievent model correcting for infection state misclassification. We found that neither Bd infection status nor infection intensity predicted mortality in juvenile frogs, counter to the expectation that early life stages are more vulnerable to disease, despite average high infection prevalence (0.35, 95% HDPI [0.14, 0.52]). Additionally, we found that observed infection prevalence and intensity were somewhat lower for juveniles than adults. Our results indicate that in this Bd-recovered species, the realized impacts of chytridiomycosis on juveniles were apparently low, likely resulting in high recruitment contributing to population stability. We highlight the importance of investigating factors relating to disease outcome in a field setting and make recommendations for future studies.

Drivers of salamander extirpation mediated by Batrachochytrium salamandrivorans.

The recent arrival of Batrachochytrium salamandrivorans in Europe was followed by rapid expansion of its geographical distribution and host range, confirming the unprecedented threat that this chytrid fungus poses to western Palaearctic amphibians. Mitigating this hazard requires a thorough understanding of the pathogen's disease ecology that is driving the extinction process. Here, we monitored infection, disease and host population dynamics in a Belgian fire salamander (Salamandra salamandra) population for two years immediately after the first signs of infection. We show that arrival of this chytrid is associated with rapid population collapse without any sign of recovery, largely due to lack of increased resistance in the surviving salamanders and a demographic shift that prevents compensation for mortality. The pathogen adopts a dual transmission strategy, with environmentally resistant non-motile spores in addition to the motile spores identified in its sister species B. dendrobatidis. The fungus retains its virulence not only in water and soil, but also in anurans and less susceptible urodelan species that function as infection reservoirs. The combined characteristics of the disease ecology suggest that further expansion of this fungus will behave as a 'perfect storm' that is able to rapidly extirpate highly susceptible salamander populations across Europe.

First evidence of ranavirus in native and invasive amphibians in Colombia.

Ranaviruses can cause mass mortality events in amphibians, thereby becoming a threat to populations that are already facing dramatic declines. Ranaviruses affect all life stages and persist in multiple amphibian hosts. The detrimental effects of ranavirus infections to amphibian populations have already been observed in the UK and in North America. In Central and South America, the virus has been reported in several countries, but the presence of the genus Ranavirus (Rv) in Colombia is unknown. To help fill this knowledge gap, we surveyed for Rv in 60 species of frogs (including one invasive species) in Colombia. We also tested for co-infection with Batrachochytrium dendrobatidis (Bd) in a subset of individuals. For Rv, we sampled 274 vouchered liver tissue samples collected between 2014 and 2019 from 41 localities covering lowlands to mountaintop páramo habitat across the country. Using quantitative polymerase chain reaction (qPCR) and end-point PCR, we detected Rv in 14 individuals from 8 localities, representing 6 species, including 5 native frogs of the genera Osornophryne, Pristimantis and Leptodactylus, and the invasive American bullfrog Rana catesbeiana. Bd was detected in 7 of 140 individuals, with 1 co-infection of Rv and Bd in an R. catesbeiana specimen collected in 2018. This constitutes the first report of ranavirus in Colombia and should set off alarms about this new emerging threat to amphibian populations in the country. Our findings provide some preliminary clues about how and when Rv may have spread and contribute to understanding how the pathogen is distributed globally.

Constant-temperature predictions underestimate growth of a fungal amphibian pathogen under individual host thermal profiles.

Ectotherm body temperatures fluctuate with environmental variability and host behavior, which may influence host-pathogen interactions. Fungal pathogens are a major threat to ectotherms and may be highly responsive to the fluctuating thermal profiles of individual hosts, especially cool-loving fungi exposed to high host temperatures. However, most studies estimate pathogen thermal performance based on averages of host or surrogate environmental temperatures, potentially missing effects of short-term host temperature shifts such as daily or hourly heat spikes. We recorded individual thermal profiles of Australian rainforest frogs using temperature-sensitive radio-transmitters. We then reproduced a subset of individual thermal profiles in growth chambers containing cultures of the near-global amphibian pathogen Batrachochytrium dendrobatidis (Bd) to investigate how realistic host temperature profiles affect Bd growth. We focused on thermal profiles that exceed the thermal optimum of Bd because the effects of realistic heat spikes on Bd growth are unresolved. Our laboratory incubation experiment revealed that Bd growth varied in response to relatively small differences in heat spike characteristics of individual frog thermal profiles, such as a single degree or a few hours, highlighting the importance of individual host behaviors in predicting population-level disease dynamics. The fungus also grew better than predicted under the most extreme and unpredictable frog temperature profile, recovering from two days of extreme (nearly 32 °C) heat spikes without negative effects on overall growth, suggesting we are underestimating the growth potential of the pathogen in nature. Combined with the previous finding that Bd reduces host heat tolerance, our study suggests that this pathogen may carry a competitive edge over hosts in the face of anthropogenic climate change.

Genetic variation of Batrachochytrium dendrobatidis is linked to skin bacterial diversity in the Pacific treefrog Hyliola regilla (hypochondriaca).

Symbiotic bacterial communities are crucial to combating infections and contribute to host health. The amphibian skin microbiome plays an important role in protecting their hosts against pathogens such as Batrachochytrium dendrobatidis (Bd), one of the causative agents of chytridiomycosis, which is responsible for dramatic amphibian population declines worldwide. Although symbiotic skin bacteria are known to inhibit Bd growth, an understanding of the relationship between Bd genetic variability, environmental conditions, and skin bacterial communities is limited. Therefore, we examined the associations between Bd infection load, Bd genetic diversity and skin bacterial communities in five populations of Hyliola regilla (hypochondriaca) from environmentally contrasting sites in Baja California, Mexico. We observed differences in Bd genetics and infection load among sites and environments. Genetic analysis of Bd isolates revealed patterns of spatial structure corresponding to the five sites sampled. Amphibian skin bacterial diversity and community structure differed among environments and sites. Bacterial community composition was correlated with Bd genetic differences and infection load, with specific bacterial taxa enriched on infected and un-infected frogs. Our results indicate that skin-associated bacteria and Bd strains likely interact on the host skin, with consequences for microbial community structure and Bd infection intensity.

Interpopulation differences in male reproductive effort drive the population dynamics of a host exposed to an emerging fungal pathogen.

Compensatory recruitment is a key demographic mechanism that has allowed the coexistence of populations of susceptible amphibians with Batrachochytrium dendrobatidis (Bd), a fungus causing one of the most devastating emerging infectious disease ever recorded among vertebrates. However, the underlying processes (e.g. density-dependent increase in survival at early life stages, change in reproductive traits) as well as the level of interpopulation variation in this response are poorly known. We explore potential mechanisms of compensatory recruitment in response to Bd infection by taking advantage of an amphibian system where male reproductive traits are easy to quantify in free-living populations. The Southern Darwin's frog Rhinoderma darwinii is a vocal sac-brooding species that exhibits a high susceptibility to lethal Bd infection. Using a 7-year capture-recapture study at four populations with contrasting Bd infection status (one high prevalence, one low prevalence and two Bd-free populations), we evaluated whether Bd-positive populations exhibited a higher adult recruitment and a higher male reproductive effort than Bd-negative populations. We also estimated population growth rates to explore whether recruitment compensated for the negative impacts of Bd on the survival of adults. In addition, we evaluated a potential demographic signal of compensatory recruitment (i.e. positive relationship between the proportion of juveniles and Bd prevalence) in response to Bd infection using raw count data from 13 R. darwinii populations. The high Bd prevalence population exhibited the highest male reproductive effort and the highest recruitment among the four monitored populations. This led to a growing population during the study period despite high mortality of adult hosts. In contrast, males from the population with low Bd prevalence had a low reproductive effort and this population, which had the lowest adult recruitment, was declining during the study period despite adults having a higher survival in comparison to the high Bd prevalence population. We also found a demographic signal of compensatory recruitment in response to Bd infection in our broader analysis of 13 R. darwinii populations. Our study underlines the importance of interpopulation variation in life-history strategies on the fate of host populations after infectious disease emergence. Our results also suggest that an increase in reproductive effort can be one of the processes underlying compensatory recruitment in populations of Bd-susceptible amphibians.

Rapid detection of amphibian chytrid fungus Batrachochytrium dendrobatidis using in situ DNA extraction and a handheld mobile thermocycler.

The amphibian chytrid fungus (Bd) has caused declines and some extinctions of amphibian populations worldwide. Early and accurate Bd detection is essential for management of susceptible anurans. We analyzed the effectiveness of in situ DNA extraction with a handheld mobile quantitative PCR (qPCR) thermocycler to detect Bd on frog skin swabs and in water samples using environmental DNA (eDNA). We collected duplicate eDNA samples and skin swabs from 3 Bd-positive Rana sierrae populations. We processed one set of samples using a field protocol (a handheld thermocycler) and the other half using a standard lab protocol. We detected Bd DNA in all R. sierrae swabbed using both the field and lab protocols. We also detected Bd DNA in eDNA samples at all sites, although the field and lab protocols failed to detect Bd eDNA at separate singular sites; results from the field and lab eDNA protocol did not match. The probability of detecting Bd DNA in the technical replicates was lower for the field protocol than the lab protocol, suggesting the field protocol has lower sensitivity and may not detect low quantities of DNA. Our results suggest that the field extraction protocol using a handheld qPCR platform is a promising tool for rapid detection of Bd in susceptible amphibian populations, yielding accurate results in less than 60 min. However, the applied field protocol may be prone to false negatives when analyzing low-quantity DNA samples such as eDNA water samples or frog swabs with low pathogen loads.

Chytrid fungus in amphibians from the lowland Brazilian Amazon.

Infectious diseases are one of the main threats to biodiversity. The fungus Batrachochytrium dendrobatidis (Bd) is associated with several amphibian losses around the globe, and environmental conditions may dictate the success of pathogen spread. The Brazilian Amazon has been considered climatically unsuitable for chytrid fungus, but additional information on Bd dynamics in this ecoregion is still lacking. We sampled 462 amphibians (449 anurans, 4 caudatans and 9 caecilians), representing 57 species from the Brazilian Amazon, and quantified Bd infections using qPCR. We tested whether abiotic variables predicted the risk of Bd infections, and tested for relationships between biotic variables and Bd. Finally, we experimentally tested the effects of Bd strains CLFT 156 and CLFT 102 (from the southern and northern Atlantic Forest, respectively) on Atelopus manauensis. We detected higher Bd prevalence than those previously reported for the Brazilian Amazon, and positive individuals in all 3 orders of amphibians sampled. Both biotic and abiotic predictors were related to prevalence, and no variable explained infection load. Moreover, we detected higher Bd prevalence in forested than open areas, while the host's reproductive biology was not a factor. We detected higher mortality in the experimental group infected with CLFT 156, probably because this strain was isolated from a region characterized by discrepant climatic conditions (latitudinally more distant) when compared with the host's sampling site in Amazon. The lowland Brazilian Amazon is still underexplored and future studies targeting all amphibian orders are essential to better understand Bd infection dynamics in this region.

Presence of invasive American bullfrogs may reduce infectious disease in a native frog species.

Amphibians breeding in aquatic environments may encounter a myriad of threats during their life cycle. One species known to prey on native amphibians in aquatic habitats is the invasive North American bullfrog Lithobates catesbeianus, which, besides being a voracious predator and competitor, often acts as a pathogen carrier and disease superspreader because it tolerates high infection loads of the frog-killing fungus Batrachochytrium dendrobatidis (Bd). Here, we hypothesized that the presence of the bullfrogs in microcosms should either (1) decrease Bd disease severity in native frog species by discouraging them from using the aquatic environment, or (2) increase the mortality of the native species. We tested these 2 mutually exclusive hypotheses by co-housing the snouted treefrog Scinax x-signatus (native to our study area) with L. catesbeianus in the laboratory, exposing them to Bd, and using qPCR analysis to quantify the resulting Bd infection loads in the native frogs. Our experiment had the following replicated treatments: (1) native-only treatment (3 individuals of S. x-signatus), (2) native-predominant treatment (2 S. x-signatus + 1 L. catesbeianus), and (3) exotic-predominant treatment (1 S. x-signatus + 2 L. catesbeianus). We found that Bd infection loads in the native S. x-signatus were highest in the native-only treatment, and lowest in the exotic-predominant treatment, indicating that bullfrogs may discourage native frogs from occupying the aquatic habitat, thus reducing encounter rates between native frogs and the waterborne pathogen. This effect could be driven by the bullfrogs' predatory behavior and their high philopatry to aquatic habitats. Our results highlight that predation risk adds to the complexity of host-species interactions in Bd epidemiology.

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.

Susceptibility of frogs to chytridiomycosis correlates with increased levels of immunomodulatory serotonin in the skin.

Chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), is a skin disease responsible for the global decline of amphibians. Frog species and populations can vary in susceptibility, but this phenomenon remains poorly understood. Here, we investigated serotonin in the skin of infected and uninfected frogs. In more susceptible frog populations, skin serotonin rose with increasing infection intensity, but decreased in later stages of the disease. The more resistant population maintained a basal level of skin serotonin. Serotonin inhibited both Bd sporangial growth and Jurkat lymphocyte proliferation in vitro. However, serotonin accumulates in skin granular glands, and this compartmentalisation may prevent inhibition of Bd growth in vivo. We suggest that skin serotonin increases in susceptible frogs due to pathogen excretion of precursor tryptophan, but that resistant frogs are able to control the levels of serotonin. Overall, the immunosuppressive effects of serotonin may contribute to the susceptibility of frogs to chytridiomycosis.

Infection dynamics, dispersal, and adaptation: understanding the lack of recovery in a remnant frog population following a disease outbreak.

Emerging infectious diseases can cause dramatic declines in wildlife populations. Sometimes, these declines are followed by recovery, but many populations do not recover. Studying differential recovery patterns may yield important information for managing disease-afflicted populations and facilitating population recoveries. In the late 1980s, a chytridiomycosis outbreak caused multiple frog species in Australia's Wet Tropics to decline. Populations of some species (e.g., Litoria nannotis) subsequently recovered, while others (e.g., Litoria dayi) did not. We examined the population genetics and current infection status of L. dayi, to test several hypotheses regarding the failure of its populations to recover: (1) a lack of individual dispersal abilities has prevented recolonization of previously occupied locations, (2) a loss of genetic variation has resulted in limited adaptive potential, and (3) L. dayi is currently adapting to chytridiomycosis. We found moderate-to-high levels of gene flow and diversity (Fst range: <0.01-0.15; minor allele frequency (MAF): 0.192-0.245), which were similar to previously published levels for recovered L. nannotis populations. This suggests that dispersal ability and genetic diversity do not limit the ability of L. dayi to recolonize upland sites. Further, infection intensity and prevalence increased with elevation, suggesting that chytridiomycosis is still limiting the elevational range of L. dayi. Outlier tests comparing infected and uninfected individuals consistently identified 18 markers as putatively under selection, and several of those markers matched genes that were previously implicated in infection. This suggests that L. dayi has genetic variation for genes that affect infection dynamics and may be undergoing adaptation.