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.

Non-Lethal Detection of Frog Virus 3-Like (RUK13) and Common Midwife Toad Virus-Like (PDE18) Ranaviruses in Two UK-Native Amphibian Species.

Ranaviruses have been involved in amphibian mass mortality events worldwide. Effective screening to control this pathogen is essential; however, current sampling methods are unsuitable for the detection of subclinical infections. Non-lethal screening is needed to prevent both further spread of ranavirus and losses of at-risk species. To assess non-lethal sampling methods, we conducted two experiments: bath exposing common frogs to RUK13 ranavirus at three concentrations, and exposing common toads to RUK13 or PDE18. Non-lethal sampling included buccal, digit, body and tank swabs, along with toe clips and stool taken across three time-points post-exposure. The presence/load of ranavirus was examined using quantitative PCR in 11 different tissues obtained from the same euthanised animals (incl. liver, gastro-intestinal tract and kidney). Buccal swab screening had the highest virus detection rate in both species (62% frogs; 71% toads) and produced consistently high virus levels compared to other non-lethal assays. The buccal swab was effective across multiple stages of infection and differing infection intensities, though low levels of infection were more difficult to detect. Buccal swab assays competed with, and even outperformed, lethal sampling in frogs and toads, respectively. Successful virus detection in the absence of clinical signs was observed (33% frogs; 50% toads); we found no difference in detectability for RUK13 and PDE18. Our results suggest that buccal swabbing could replace lethal sampling for screening and be introduced as standard practice for ranavirus surveillance.

Climate anomalies and competition reduce establishment success during island colonization.

Understanding the factors that facilitate or constrain establishment of populations in novel environments is crucial for conservation biology and the study of adaptive radiation. Important questions include: (1) Does the timing of colonization relative to stochastic events, such as climatic perturbations, impact the probability of successful establishment? (2) To what extent does community context (e.g., the presence of competitors) change the probability of establishment? (3) How do sources of intrapopulation variance, such as sex differences, affect success at an individual level during the process of establishment? Answers to these questions are rarely pursued in a field-experimental context or on the same time scales (months to years) as the processes of colonization and establishment. We introduced slender anole lizards (Anolis apletophallus) to eight islands in the Panama Canal and tracked them over multiple generations to investigate the factors that mediate establishment success. All islands were warmer than the mainland (ancestral) environment, and some islands had a native competitor. We transplanted half of these populations only 4 months before the onset of a severe regional drought and the other half 2 years (two generations) before the drought. We found that successful establishment depended on both the intensity of interspecific competition and the timing of colonization relative to the drought. The islands that were colonized shortly before the drought went functionally extinct by the second generation, and regardless of time before the drought, the populations on islands with interspecific competition declined continuously over the study period. Furthermore, the effect of the competitor interacted with sex, with males suffering, and females benefitting, from the presence of a native competitor. Our results reveal that community context and the timing of colonization relative to climactic events can combine to determine establishment success and that these factors can generate opposite effects on males and females.

Host Identity Matters-Up to a Point: The Community Context of Batrachochytrium dendrobatidis Transmission.

AbstractThe level of detail on host communities needed to understand multihost parasite invasions is an unresolved issue in disease ecology. Coarse community metrics that ignore functional differences between hosts, such as host species richness, can be good predictors of invasion outcomes. Yet if host species vary in the extent to which they maintain and transmit infections, then explicitly accounting for those differences may be important. Through controlled mesocosm experiments and modeling, we show that interspecific differences between host species are important for community-wide infection dynamics of the multihost fungal parasite of amphibians (Batrachochytrium dendrobatidis [Bd]), but only up to a point. The most abundant host species in our system, fire salamander larvae (Salamandra salamandra), did not maintain or transmit infections. Rather, two less abundant "auxiliary" host species, Iberian tree frog (Hyla molleri) and spiny toad (Bufo spinosus) larvae, maintained and transmitted Bd. Frogs had the highest mean rates of Bd shedding, giving them the highest contributions to the basic reproduction number, R0. Toad contributions to R0 were substantial, however, and when examining community-level patterns of infection and transmission, the effects of frogs and toads were similar. Specifying more than just host species richness to distinguish salamanders from auxiliary host species was critical for predicting community-level Bd prevalence and transmission. Distinguishing frogs from toads, however, did not improve predictions. These findings demonstrate limitations to the importance of host species identities in multihost infection dynamics. Host species that exhibit different functional traits, such as susceptibility and infectiousness, may play similar epidemiological roles in the broader community.

Microbiome function predicts amphibian chytridiomycosis disease dynamics.

BACKGROUND: The fungal pathogen Batrachochytrium dendrobatidis (Bd) threatens amphibian biodiversity and ecosystem stability worldwide. Amphibian skin microbial community structure has been linked to the clinical outcome of Bd infections, yet its overall functional importance is poorly understood. METHODS: Microbiome taxonomic and functional profiles were assessed using high-throughput bacterial 16S rRNA and fungal ITS2 gene sequencing, bacterial shotgun metagenomics and skin mucosal metabolomics. We sampled 56 wild midwife toads (Alytes obstetricans) from montane populations exhibiting Bd epizootic or enzootic disease dynamics. In addition, to assess whether disease-specific microbiome profiles were linked to microbe-mediated protection or Bd-induced perturbation, we performed a laboratory Bd challenge experiment whereby 40 young adult A. obstetricans were exposed to Bd or a control sham infection. We measured temporal changes in the microbiome as well as functional profiles of Bd-exposed and control animals at peak infection. RESULTS: Microbiome community structure and function differed in wild populations based on infection history and in experimental control versus Bd-exposed animals. Bd exposure in the laboratory resulted in dynamic changes in microbiome community structure and functional differences, with infection clearance in all but one infected animal. Sphingobacterium, Stenotrophomonas and an unclassified Commamonadaceae were associated with wild epizootic dynamics and also had reduced abundance in laboratory Bd-exposed animals that cleared infection, indicating a negative association with Bd resistance. This was further supported by microbe-metabolite integration which identified functionally relevant taxa driving disease outcome, of which Sphingobacterium and Bd were most influential in wild epizootic dynamics. The strong correlation between microbial taxonomic community composition and skin metabolome in the laboratory and field is inconsistent with microbial functional redundancy, indicating that differences in microbial taxonomy drive functional variation. Shotgun metagenomic analyses support these findings, with similar disease-associated patterns in beta diversity. Analysis of differentially abundant bacterial genes and pathways indicated that bacterial environmental sensing and Bd resource competition are likely to be important in driving infection outcomes. CONCLUSIONS: Bd infection drives altered microbiome taxonomic and functional profiles across laboratory and field environments. Our application of multi-omics analyses in experimental and field settings robustly predicts Bd disease dynamics and identifies novel candidate biomarkers of infection. Video Abstract.

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

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

Exposure to Batrachochytrium dendrobatidis affects chemical defences in two anuran amphibians, Rana dalmatina and Bufo bufo.

BACKGROUND: Batrachochytrium dendrobatidis (Bd) is the causative agent of chytridiomycosis, one of the major causes of worldwide amphibian biodiversity loss. Many amphibians exhibit skin-based chemical defences, which may play an important role against invading pathogens, but whether the synthesis of these chemical compounds is enhanced or suppressed in the presence of pathogens is largely unknown. Here we investigated direct and indirect effects of larval exposure to the globally distributed and highly virulent Bd-GPL strain on skin secreted chemical defences and life history traits during early ontogeny of agile frogs (Rana dalmatina) and common toads (Bufo bufo). RESULTS: Exposure to Bd during the larval stage did not result in enhanced synthesis of the antimicrobial peptide Brevinin-1 Da in R. dalmatina tadpoles or in increased production of bufadienolides in B. bufo tadpoles. However, exposure to Bd during the larval stage had a carry-over effect reaching beyond metamorphosis: both R. dalmatina and B. bufo froglets contained smaller quantities of defensive chemicals than their Bd-naïve conspecifics in the control treatment. Prevalence of Bd and infection intensities were very low in both larvae and metamorphs of R. dalmatina, while in B. bufo we observed high Bd prevalence and infection intensities, especially in metamorphs. At the same time, we did not find a significant effect of Bd-exposure on body mass or development rate in larvae or metamorphs in either species. CONCLUSIONS: The lack of detrimental effect of Bd-exposure on life history traits, even parallel with high infection intensities in the case of B. bufo individuals, is surprising and suggests high tolerance of local populations of these two species against Bd. However, the lowered quantity of defensive chemicals may compromise antimicrobial and antipredatory defences of froglets, which may ultimately contribute to population declines also in the absence of conspicuous mass-mortality events.

Discriminating lineages of Batrachochytrium dendrobatidis using quantitative PCR.

The ability to detect and monitor infectious disease in a phylogenetically informative manner is critical for their management. Phylogenetically informative diagnostic tests enable patterns of pathogen introduction or changes in the distribution of genotypes to be measured, enabling research into the ecology of the pathogen. Batrachochytrium dendrobatidis (Bd), a causative agent of chytridiomycosis in amphibian populations, emerged worldwide in the 21st century and is composed of six lineages which are display varying levels of virulence in their hosts. Research into the distribution, ecology and pathogenicity of these lineages has been hampered by an inability to type lineage efficiently. Here, we describe a lineage-specific TaqMan qPCR assay that differentiates the two lineages of Bd most commonly associated with chytridiomycosis: BdGPL and BdCAPE. We demonstrate how this assay can be used for the surveillance of wild populations of amphibians in Southern Africa using skin swabs, tissue samples and cultured isolates.

Significant reductions of host abundance weakly impact infection intensity of Batrachochytrium dendrobatidis.

Infectious diseases are considered major threats to biodiversity, however strategies to mitigate their impacts in the natural world are scarce and largely unsuccessful. Chytridiomycosis is responsible for the decline of hundreds of amphibian species worldwide, but an effective disease management strategy that could be applied across natural habitats is still lacking. In general amphibian larvae can be easily captured, offering opportunities to ascertain the impact of altering the abundance of hosts, considered to be a key parameter affecting the severity of the disease. Here, we report the results of two experiments to investigate how altering host abundance affects infection intensity in amphibian populations of a montane area of Central Spain suffering from lethal amphibian chytridiomycosis. Our laboratory-based experiment supported the conclusion that varying density had a significant effect on infection intensity when salamander larvae were housed at low densities. Our field experiment showed that reducing the abundance of salamander larvae in the field also had a significant, but weak, impact on infection the following year, but only when removals were extreme. While this suggests adjusting host abundance as a mitigation strategy to reduce infection intensity could be useful, our evidence suggests only heavy culling efforts will succeed, which may run contrary to objectives for conservation.

Occurrence of Batrachochytrium dendrobatidis in Sweden: higher infection prevalence in southern species.

The chytrid fungus Batrachochytrium dendrobatidis (Bd) has caused worldwide declines in amphibian populations. While Bd is widespread in southern and central Europe, its occurrence and distribution in northernmost Europe is mostly unknown. We surveyed for Bd in breeding anurans in Sweden by sampling 1917 amphibians from 101 localities and 3 regions in Sweden (southern, northern and central). We found that Bd was widespread in southern and central Sweden, occurring in all 9 investigated species and in 45.5% of the 101 localities with an overall prevalence of 13.8%. No infected individuals were found in the 4 northern sites sampled. The records from central Sweden represent the northernmost records of Bd in Europe. While the proportion of sites positive for Bd was similar between the southern and central regions, prevalence was much higher in the southern region. This was because southern species with a distribution mainly restricted to southernmost Sweden had a higher prevalence than widespread generalist species. The nationally red-listed green toad Bufotes variabilis and the fire-bellied toad Bombina bombina had the highest prevalence (61.4 and 48.9%, respectively). Across species, Bd prevalence was strongly positively, correlated with water temperature at the start of egg laying. However, no individuals showing visual signs of chytridiomycosis were found in the field. These results indicate that Bd is widespread and common in southern and central Sweden with southern species, breeding in higher temperatures and with longer breeding periods, having higher prevalence. However, the impact of Bd on amphibian populations in northernmost Europe remains unknown.

Genetic and demographic vulnerability of adder populations: Results of a genetic study in mainland Britain.

Genetic factors are often overlooked in conservation planning, despite their importance in small isolated populations. We used mitochondrial and microsatellite markers to investigate population genetics of the adder (Vipera berus) in southern Britain, where numbers are declining. We found no evidence for loss of heterozygosity in any of the populations studied. Genetic diversity was comparable across sites, in line with published levels for mainland Europe. However, further analysis revealed a striking level of relatedness. Genetic networks constructed from inferred first degree relationships suggested a high proportion of individuals to be related at a level equivalent to that of half-siblings, with rare inferred full-sib dyads. These patterns of relatedness can be attributed to the high philopatry and low vagility of adders, which creates high local relatedness, in combination with the polyandrous breeding system in the adder, which may offset the risk of inbreeding in closed populations. We suggest that reliance on standard genetic indicators of inbreeding and diversity may underestimate demographic and genetic factors that make adder populations vulnerable to extirpation. We stress the importance of an integrated genetic and demographic approach in the conservation of adders, and other taxa of similar ecology.

Response to Comment on "Amphibian fungal panzootic causes catastrophic and ongoing loss of biodiversity".

Lambert et al question our retrospective and holistic epidemiological assessment of the role of chytridiomycosis in amphibian declines. Their alternative assessment is narrow and provides an incomplete evaluation of evidence. Adopting this approach limits understanding of infectious disease impacts and hampers conservation efforts. We reaffirm that our study provides unambiguous evidence that chytridiomycosis has affected at least 501 amphibian species.

Chytrid fungi and global amphibian declines.

Discovering that chytrid fungi cause chytridiomycosis in amphibians represented a paradigm shift in our understanding of how emerging infectious diseases contribute to global patterns of biodiversity loss. In this Review we describe how the use of multidisciplinary biological approaches has been essential to pinpointing the origins of amphibian-parasitizing chytrid fungi, including Batrachochytrium dendrobatidis and Batrachochytrium salamandrivorans, as well as to timing their emergence, tracking their cycles of expansion and identifying the core mechanisms that underpin their pathogenicity. We discuss the development of the experimental methods and bioinformatics toolkits that have provided a fuller understanding of batrachochytrid biology and informed policy and control measures.

Spatiotemporal heterogeneity decouples infection parameters of amphibian chytridiomycosis.

Emerging infectious diseases are responsible for declines in wildlife populations around the globe. Mass mortality events associated with emerging infectious diseases are often associated with high number of infected individuals (prevalence) and high pathogen loads within individuals (intensity). At the landscape scale, spatial and temporal variation in environmental conditions can alter the relationship between these infection parameters and blur the overall picture of disease dynamics. Quantitative estimates of how infection parameters covary with environmental heterogeneity at the landscape scale are scarce. If we are to identify wild populations at risk of disease epidemics, we must elucidate the factors that shape, and potentially decouple, the link between pathogen prevalence and intensity of infection over complex ecological landscapes. Using a network of 41 populations of the amphibian host Rana pipiens in Ontario, Canada, we present the spatial and temporal heterogeneity in pathogen prevalence and intensity of infection of the chytrid fungus Batrachochytrium dendrobatidis (Bd), across a 3-year period. We then quantify how covariation between both infection parameters measured during late summer is modified by previously experienced spatiotemporal environmental heterogeneity across 14 repeat sampled populations. Late summer Bd infection parameters are governed, at least in part, by different environmental factors operating during separate host life-history events. Our results provide evidence for a relationship between Bd prevalence and thermal regimes prior to host breeding at the site level, and a relationship between intensity of infection and aquatic conditions (precipitation, hydroshed size and river density) throughout host breeding period at the site level. This demonstrates that microclimatic variation within temporal windows can drive divergent patterns of pathogen dynamics within and across years, by effecting changes in host behaviour which interfere with the pathogen's ability to infect and re-infect hosts. A clearer understanding of the role that spatiotemporal heterogeneity has upon infection parameters will provide valuable insights into host-pathogen epidemiology, as well as more fundamental aspects of the ecology and evolution of interspecific interactions.

Outbreaks of an Emerging Viral Disease Covary With Differences in the Composition of the Skin Microbiome of a Wild United Kingdom Amphibian.

There is growing appreciation of the important role of commensal microbes in ensuring the normal function and health of their hosts, including determining how hosts respond to pathogens. A range of infectious diseases are threatening amphibians worldwide, and evidence is accumulating that the host-associated bacteria that comprise the microbiome may be key in mediating interactions between amphibian hosts and infectious pathogens. We used 16S rRNA amplicon sequencing to quantify the skin microbial community structure of over 200 individual wild adult European common frogs (Rana temporaria), from ten populations with contrasting history of the lethal disease ranavirosis, caused by emerging viral pathogens belonging to the genus Ranavirus. All populations had similar species richness irrespective of disease history, but populations that have experienced historical outbreaks of ranavirosis have a distinct skin microbiome structure (beta diversity) when compared to sites where no outbreaks of the disease have occurred. At the individual level, neither age, body length, nor sex of the frog could predict the structure of the skin microbiota. Our data potentially support the hypothesis that variation among individuals in skin microbiome structure drive differences in susceptibility to infection and lethal outbreaks of disease. More generally, our results suggest that population-level processes are more important for driving differences in microbiome structure than variation among individuals within populations in key life history traits such as age and body size.

Modelling Ranavirus Transmission in Populations of Common Frogs (Rana temporaria) in the United Kingdom.

Ranaviruses began emerging in common frogs (Rana temporaria) in the United Kingdom in the late 1980s and early 1990s, causing severe disease and declines in the populations of these animals. Herein, we explored the transmission dynamics of the ranavirus(es) present in common frog populations, in the context of a simple susceptible-infected (SI) model, using parameters derived from the literature. We explored the effects of disease-induced population decline on the dynamics of the ranavirus. We then extended the model to consider the infection dynamics in populations exposed to both ulcerative and hemorrhagic forms of the ranaviral disease. The preliminary investigation indicated the important interactions between the forms. When the ulcerative form was present in a population and the hemorrhagic form was later introduced, the hemorrhagic form of the disease needed to be highly contagious, to persist. We highlighted the areas where further research and experimental evidence is needed and hope that these models would act as a guide for further research into the amphibian disease dynamics.

Effects of historic and projected climate change on the range and impacts of an emerging wildlife disease.

The global trend of increasing environmental temperatures is often predicted to result in more severe disease epidemics. However, unambiguous evidence that temperature is a driver of epidemics is largely lacking, because it is demanding to demonstrate its role among the complex interactions between hosts, pathogens, and their shared environment. Here, we apply a three-pronged approach to understand the effects of temperature on ranavirus epidemics in UK common frogs, combining in vitro, in vivo, and field studies. Each approach suggests that higher temperatures drive increasing severity of epidemics. In wild populations, ranavirosis incidents were more frequent and more severe at higher temperatures, and their frequency increased through a period of historic warming in the 1990s. Laboratory experiments using cell culture and whole animal models showed that higher temperature increased ranavirus propagation, disease incidence, and mortality rate. These results, combined with climate projections, predict severe ranavirosis outbreaks will occur over wider areas and an extended season, possibly affecting larval recruitment. Since ranaviruses affect a variety of ectothermic hosts (amphibians, reptiles, and fish), wider ecological damage could occur. Our three complementary lines of evidence present a clear case for direct environmental modulation of these epidemics and suggest management options to protect species from disease.

Common midwife toad ranaviruses replicate first in the oral cavity of smooth newts (Lissotriton vulgaris) and show distinct strain-associated pathogenicity.

Ranavirus is the second most common infectious cause of amphibian mortality. These viruses affect caudates, an order in which information regarding Ranavirus pathogenesis is scarce. In the Netherlands, two strains (CMTV-NL I and III) were suspected to possess distinct pathogenicity based on field data. To investigate susceptibility and disease progression in urodeles and determine differences in pathogenicity between strains, 45 adult smooth newts (Lissotriton vulgaris) were challenged via bath exposure with these ranaviruses and their detection in organs and feces followed over time by PCR, immunohistochemistry and in situ hybridization. Ranavirus was first detected at 3 days post infection (p.i.) in the oral cavity and upper respiratory mucosa. At 6 days p.i, virus was found in connective tissues and vasculature of the gastrointestinal tract. Finally, from 9 days p.i onwards there was widespread Ranavirus disease in various organs including skin, kidneys and gonads. Higher pathogenicity of the CMTV-NL I strain was confirmed by higher correlation coefficient of experimental group and mortality of challenged animals. Ranavirus-exposed smooth newts shed virus in feces intermittently and infection was seen in the absence of lesions or clinical signs, indicating that this species can harbor subclinical infections and potentially serve as disease reservoirs.

Amphibian fungal panzootic causes catastrophic and ongoing loss of biodiversity.

Anthropogenic trade and development have broken down dispersal barriers, facilitating the spread of diseases that threaten Earth's biodiversity. We present a global, quantitative assessment of the amphibian chytridiomycosis panzootic, one of the most impactful examples of disease spread, and demonstrate its role in the decline of at least 501 amphibian species over the past half-century, including 90 presumed extinctions. The effects of chytridiomycosis have been greatest in large-bodied, range-restricted anurans in wet climates in the Americas and Australia. Declines peaked in the 1980s, and only 12% of declined species show signs of recovery, whereas 39% are experiencing ongoing decline. There is risk of further chytridiomycosis outbreaks in new areas. The chytridiomycosis panzootic represents the greatest recorded loss of biodiversity attributable to a disease.

Diversity-Stability Dynamics of the Amphibian Skin Microbiome and Susceptibility to a Lethal Viral Pathogen.

Variation among animals in their host-associated microbial communities is increasingly recognized as a key determinant of important life history traits including growth, metabolism, and resistance to disease. Quantitative estimates of the factors shaping the stability of host microbiomes over time at the individual level in non-model organisms are scarce. Addressing this gap in our knowledge is important, as variation among individuals in microbiome stability may represent temporal gain or loss of key microbial species and functions linked to host health and/or fitness. Here we use controlled experiments to investigate how both heterogeneity in microbial species richness of the environment and exposure to the emerging pathogen Ranavirus influence the structure and temporal dynamics of the skin microbiome in a vertebrate host, the European common frog (Rana temporaria). Our evidence suggests that altering the bacterial species richness of the environment drives divergent temporal microbiome dynamics of the amphibian skin. Exposure to ranavirus effects changes in skin microbiome structure irrespective of total microbial diversity, but individuals with higher pre-exposure skin microbiome diversity appeared to exhibit higher survival. Higher diversity skin microbiomes also appear less stable over time compared to lower diversity microbiomes, but stability of the 100 most abundant ("core") community members was similar irrespective of microbiome richness. Our study highlights the importance of extrinsic factors in determining the stability of host microbiomes over time, which may in turn have important consequences for the stability of host-microbe interactions and microbiome-fitness correlations.