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.

Response to "Specimen collection is essential for modern science".

Natural history museums are vital repositories of specimens, samples and data that inform about the natural world; this Reply responds to a Formal Comment that queried whether it will ever be possible to completely do away with whole animal specimen collection, inviting open conversations about the ethical implications of specimen collection.

Reimagining the future of natural history museums with compassionate collection.

Compassionate collection involves minimizing harm while collecting museum data in the field. By adopting this practice, natural history museums could better maintain existing collections, accommodate more nonlethal specimens and data, and foster an inclusive community.

The fungal pathogen Batrachochytrium salamandrivorans is not detected in wild and captive amphibians from Mexico.

The recent emergence of the pathogen Batrachochytrium salamandrivorans (Bsal) is associated with rapid population declines of salamanders in Europe and its arrival to new areas could cause dramatic negative effects on other amphibian populations and species. Amphibian species, present in areas with high amphibian diversity such as Mexico, could be highly threatened due to the arrival of Bsal, particularly salamander species which are more vulnerable to chytridiomycosis caused by this pathogen. Thus, immediate surveillance is needed as a strategy to efficiently contend with this emerging infectious disease. In this study, we analyzed 490 wild and captive amphibians from 48 species across 76 sites in the North, Central, and South of Mexico to evaluate the presence of Bsal. Amphibians were sampled in sites with variable degrees of amphibian richness and suitability for Bsal according to previous studies. From the 76 sampling sites, 10 of them were located in areas with high amphibian richness and potential moderate to high Bsal habitat suitability. We did not detect Bsal in any of the samples, and no signs of the disease were observed in any individual at the time of sampling. Our results suggest that Bsal has not yet arrived at the sampled sites or could be at low prevalence within populations with low occurrence probability. This is the first study that evaluates the presence of Bsal in different regions and amphibian species in Mexico, which is the second most diverse country in salamander species in the world. We highlight the risk and the importance of continuing surveillance of Bsal in Mexico and discuss control strategies to avoid the introduction and spread of Bsal in the country.

Host species is linked to pathogen genotype for the amphibian chytrid fungus (Batrachochytrium dendrobatidis).

Host-pathogen specificity can arise from certain selective environments mediated by both the host and pathogen. Therefore, understanding the degree to which host species identity is correlated with pathogen genotype can help reveal historical host-pathogen dynamics. One animal disease of particular concern is chytridiomycosis, typically caused by the global panzootic lineage of the amphibian chytrid fungus (Batrachochytrium dendrobatidis, Bd), termed Bd-GPL. This pathogen lineage has caused devastating declines in amphibian communities around the world. However, the site of origin for the common ancestor of modern Bd-GPL and the fine-scale transmission dynamics of this lineage have remained a mystery. This is especially the case in North America where Bd-GPL is widespread, but disease outbreaks occur sporadically. Herein, we use Bd genetic data collected throughout the United States from amphibian skin swabs and cultured isolate samples to investigate Bd genetic patterns. We highlight two case studies in Pennsylvania and Nevada where Bd-GPL genotypes are strongly correlated with host species identity. Specifically, in some localities bullfrogs (Rana catesbeiana) are infected with Bd-GPL lineages that are distinct from those infecting other sympatric amphibian species. Overall, we reveal a previously unknown association of Bd genotype with host species and identify the eastern United States as a Bd diversity hotspot and potential site of origin for Bd-GPL.

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.

Thermal Performance Curves of Multiple Isolates of Batrachochytrium dendrobatidis, a Lethal Pathogen of Amphibians.

Emerging infectious disease is a key factor in the loss of amphibian diversity. In particular, the disease chytridiomycosis has caused severe declines around the world. The lethal fungal pathogen that causes chytridiomycosis, Batrachochytrium dendrobatidis (Bd), has affected amphibians in many different environments. One primary question for researchers grappling with disease-induced losses of amphibian biodiversity is what abiotic factors drive Bd pathogenicity in different environments. To study environmental influences on Bd pathogenicity, we quantified responses of Bd phenotypic traits (e.g., viability, zoospore densities, growth rates, and carrying capacities) over a range of environmental temperatures to generate thermal performance curves. We selected multiple Bd isolates that belong to a single genetic lineage but that were collected across a latitudinal gradient. For the population viability, we found that the isolates had similar thermal optima at 21°C, but there was considerable variation among the isolates in maximum viability at that temperature. Additionally, we found the densities of infectious zoospores varied among isolates across all temperatures. Our results suggest that temperatures across geographic point of origin (latitude) may explain some of the variation in Bd viability through vertical shifts in maximal performance. However, the same pattern was not evident for other reproductive parameters (zoospore densities, growth rates, fecundity), underscoring the importance of measuring multiple traits to understand variation in pathogen responses to environmental conditions. We suggest that variation among Bd genetic variants due to environmental factors may be an important determinant of disease dynamics for amphibians across a range of diverse environments.

Divergent regional evolutionary histories of a devastating global amphibian pathogen.

Emerging infectious diseases are a pressing threat to global biological diversity. Increased incidence and severity of novel pathogens underscores the need for methodological advances to understand pathogen emergence and spread. Here, we use genetic epidemiology to test, and challenge, key hypotheses about a devastating zoonotic disease impacting amphibians globally. Using an amplicon-based sequencing method and non-invasive samples we retrospectively explore the history of the fungal pathogen Batrachochytrium dendrobatidis (Bd) in two emblematic amphibian systems: the Sierra Nevada of California and Central Panama. The hypothesis in both regions is the hypervirulent Global Panzootic Lineage of Bd (BdGPL) was recently introduced and spread rapidly in a wave-like pattern. Our data challenge this hypothesis by demonstrating similar epizootic signatures can have radically different underlying evolutionary histories. In Central Panama, our genetic data confirm a recent and rapid pathogen spread. However, BdGPL in the Sierra Nevada has remarkable spatial structuring, high genetic diversity and a relatively older history inferred from time-dated phylogenies. Thus, this deadly pathogen lineage may have a longer history in some regions than assumed, providing insights into its origin and spread. Overall, our results highlight the importance of integrating observed wildlife die-offs with genetic data to more accurately reconstruct pathogen outbreaks.

Early presence of Batrachochytrium dendrobatidis in Mexico with a contemporary dominance of the global panzootic lineage.

Chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), is a devastating infectious disease of amphibians. Retrospective studies using museum vouchers and genetic samples supported the hypothesis that Bd colonized Mexico from North America and then continued to spread into Central and South America, where it led to dramatic losses in tropical amphibian biodiversity (the epizootic wave hypothesis). While these studies suggest that Bd has been in Mexico since the 1970s, information regarding the historical and contemporary occurrence of different pathogen genetic lineages across the country is limited. In the current study, we investigated the historical and contemporary patterns of Bd in Mexico. We combined the swabbing of historical museum vouchers and sampling of wild amphibians with a custom Bd genotyping assay to assess the presence, prevalence, and genetic diversity of Bd over time in Mexico. We found Bd-positive museum specimens from the late 1800s, far earlier than previous records and well before recent amphibian declines. With Bd genotypes from samples collected between 1975-2019, we observed a contemporary dominance of the global panzootic lineage in Mexico and report four genetic subpopulations and potential for admixture among these populations. The observed genetic variation did not have a clear geographic signature or provide clear support for the epizootic wave hypothesis. These results provide a framework for testing new questions regarding Bd invasions and their temporal relationship to observed amphibian declines in the Americas.

Whole exome sequencing identifies the potential for genetic rescue in iconic and critically endangered Panamanian harlequin frogs.

Avoiding extinction in a rapidly changing environment often relies on a species' ability to quickly adapt in the face of extreme selective pressures. In Panamá, two closely related harlequin frog species (Atelopus varius and Atelopus zeteki) are threatened with extinction due to the fungal pathogen Batrachochytrium dendrobatidis (Bd). Once thought to be nearly extirpated from Panamá, A. varius have recently been rediscovered in multiple localities across their historical range; however, A. zeteki are possibly extinct in the wild. By leveraging a unique collection of 186 Atelopus tissue samples collected before and after the Bd outbreak in Panama, we describe the genetics of persistence for these species on the brink of extinction. We sequenced the transcriptome and developed an exome-capture assay to sequence the coding regions of the Atelopus genome. Using these genetic data, we evaluate the population genetic structure of historical A. varius and A. zeteki populations, describe changes in genetic diversity over time, assess the relationship between contemporary and historical individuals, and test the hypothesis that some A. varius populations have rapidly evolved to resist or tolerate Bd infection. We found a significant decrease in genetic diversity in contemporary (compared to historical) A. varius populations. We did not find strong evidence of directional allele frequency change or selection for Bd resistance genes, but we uncovered a set of candidate genes that warrant further study. Additionally, we found preliminary evidence of recent migration and gene flow in one of the largest persisting A. varius populations in Panamá, suggesting the potential for genetic rescue in this system. Finally, we propose that previous conservation units should be modified, as clear genetic breaks do not exist beyond the local population level. Our data lay the groundwork for genetically informed conservation and advance our understanding of how imperiled species might be rescued from extinction.

Fungal infection, decline and persistence in the only obligate troglodytic Neotropical salamander.

The fungal pathogen Batrachochytrium dendrobatidis (Bd) is implicated in global mass die-offs and declines in amphibians. In Mesoamerica, the Bd epidemic wave hypothesis is supported by detection of Bd in historic museum specimens collected over the last century, yet the timing and impact of the early stages of the wave remain poorly understood. Chiropterotriton magnipes, the only obligate troglodytic Neotropical salamander, was abundant in its small range in the decade following its description in 1965, but subsequently disappeared from known localities and was not seen for 34 years. Its decline is roughly coincident with that of other populations of Neotropical salamanders associated with the invasion and spread of Bd. To determine the presence and infection intensity of Bd on C. magnipes and sympatric amphibian species (which are also Bd hosts), we used a noninvasive sampling technique and qPCR assay to detect Bd on museum specimens of C. magnipes collected from 1952 to 2012, and from extant populations of C. magnipes and sympatric species of amphibians. We also tested for the presence of the recently discovered Batrachochytrium salamandivorans (Bsal), another fungal chytridiomycete pathogen of salamanders, using a similar technique specific for Bsal. We did not detect Bd in populations of C. magnipes before 1969, while Bd was detected at low to moderate prevalence just prior to and during declines. This pattern is consistent with Bd-caused epizootics followed by host declines and extirpations described in other hosts. We did not detect Bsal in any extant population of C. magnipes. We obtained one of the earliest positive records of the fungus to date in Latin America, providing additional historical evidence consistent with the Bd epidemic wave hypothesis. Genotyping results show that at least one population is currently infected with the Global Panzootic Lineage of Bd, but our genotyping of the historical positive samples was unsuccessful. The lack of large samples from some years and the difficulty in genotyping historical Bd samples illustrate some of the difficulties inherent in assigning causality to historical amphibian declines. These data also provide an important historical baseline for actions to preserve the few known remaining populations of C. magnipes.

Invasive vegetation affects amphibian skin microbiota and body condition.

Invasive plants are major drivers of habitat modification and the scale of their impact is increasing globally as anthropogenic activities facilitate their spread. In California, an invasive plant genus of great concern is Eucalyptus. Eucalyptus leaves can alter soil chemistry and negatively affect underground macro- and microbial communities. Amphibians serve as excellent models to evaluate the effect of Eucalyptus invasion on ground-dwelling species as they predate on soil arthropods and incorporate soil microbes into their microbiotas. The skin microbiota is particularly important to amphibian health, suggesting that invasive plant species could ultimately affect amphibian populations. To investigate the potential for invasive vegetation to induce changes in microbial communities, we sampled microbial communities in the soil and on the skin of local amphibians. Specifically, we compared Batrachoseps attenuatus skin microbiomes in both Eucalyptus globulus (Myrtaceae) and native Quercus agriflolia (Fagaceae) dominated forests in the San Francisco Bay Area. We determined whether changes in microbial diversity and composition in both soil and Batrachoseps attenuatus skin were associated with dominant vegetation type. To evaluate animal health across vegetation types, we compared Batrachoseps attenuatus body condition and the presence/absence of the amphibian skin pathogen Batrachochytrium dendrobatidis. We found that Eucalyptus invasion had no measurable effect on soil microbial community diversity and a relatively small effect (compared to the effect of site identity) on community structure in the microhabitats sampled. In contrast, our results show that Batrachoseps attenuatus skin microbiota diversity was greater in Quercus dominated habitats. One amplicon sequence variant identified in the family Chlamydiaceae was observed in higher relative abundance among salamanders sampled in Eucalyptus dominated habitats. We also observed that Batrachoseps attenuatus body condition was higher in Quercus dominated habitats. Incidence of Batrachochytrium dendrobatidis across all individuals was very low (only one Batrachochytrium dendrobatidis positive individual). The effect on body condition demonstrates that although Eucalyptus may not always decrease amphibian abundance or diversity, it can potentially have cryptic negative effects. Our findings prompt further work to determine the mechanisms that lead to changes in the health and microbiome of native species post-plant invasion.

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

Scheele et al (Reports, 29 March 2019, p. 1459) bring needed attention to the effects of amphibian infectious disease. However, the data and methods implicating the disease chytridiomycosis in 501 amphibian species declines are deficient. Which species are affected, and how many, remains a critical unanswered question. Amphibians are imperiled; protective actions require public support and robust science.

Cryptic diversity of a widespread global pathogen reveals expanded threats to amphibian conservation.

Biodiversity loss is one major outcome of human-mediated ecosystem disturbance. One way that humans have triggered wildlife declines is by transporting disease-causing agents to remote areas of the world. Amphibians have been hit particularly hard by disease due in part to a globally distributed pathogenic chytrid fungus (Batrachochytrium dendrobatidis [Bd]). Prior research has revealed important insights into the biology and distribution of Bd; however, there are still many outstanding questions in this system. Although we know that there are multiple divergent lineages of Bd that differ in pathogenicity, we know little about how these lineages are distributed around the world and where lineages may be coming into contact. Here, we implement a custom genotyping method for a global set of Bd samples. This method is optimized to amplify and sequence degraded DNA from noninvasive skin swab samples. We describe a divergent lineage of Bd, which we call BdASIA3, that appears to be widespread in Southeast Asia. This lineage co-occurs with the global panzootic lineage (BdGPL) in multiple localities. Additionally, we shed light on the global distribution of BdGPL and highlight the expanded range of another lineage, BdCAPE. Finally, we argue that more monitoring needs to take place where Bd lineages are coming into contact and where we know little about Bd lineage diversity. Monitoring need not use expensive or difficult field techniques but can use archived swab samples to further explore the history-and predict the future impacts-of this devastating pathogen.

Opening the file drawer: Unexpected insights from a chytrid infection experiment.

Infection experiments are critical for understanding wildlife disease dynamics. Although infection experiments are typically designed to reduce complexity, disease outcomes still result from complex interactions between host, pathogen, and environmental factors. Cryptic variation across factors can lead to decreased repeatability of infection experiments within and between research groups and hinder research progress. Furthermore, studies with unexpected results are often relegated to the "file drawer" and potential insights gained from these experimental outcomes are lost. Here, we report unexpected results from an infection experiment studying the response of two differentially-susceptible but related frogs (American Bullfrog Rana catesbeiana and the Mountain yellow-legged frog Rana muscosa) to the amphibian-killing chytrid fungus (Batrachochytrium dendrobatidis, Bd). Despite well-documented differences in susceptibility between species, we found no evidence for antibody-mediated immune response and no Bd-related mortality in either species. Additionally, during the study, the sham-inoculated R. catesbeiana control group became unexpectedly Bd-positive. We used a custom genotyping assay to demonstrate that the aberrantly-infected R. catesbeiana carried a Bd genotype distinct from the inoculation genotype. Thus R. catesbeiana individuals were acquired with low-intensity infections that could not be detected with qPCR. In the Bd-inoculated R. catesbeiana treatment group, the inoculated genotype appeared to out-compete the cryptic infection. Thus, our results provide insight into Bd coinfection dynamics, a phenomenon that is increasingly relevant as different pathogen strains are moved around the globe. Our experiment highlights how unexpected experimental outcomes can serve as both cautionary tales and opportunities to explore unanswered research questions. We use our results as a case study to highlight common sources of anomalous results for infection experiments. We argue that understanding these factors will aid researchers in the design, execution, and interpretation of experiments to understand wildlife disease processes.

Shifts in disease dynamics in a tropical amphibian assemblage are not due to pathogen attenuation.

Infectious diseases rarely end in extinction. Yet the mechanisms that explain how epidemics subside are difficult to pinpoint. We investigated host-pathogen interactions after the emergence of a lethal fungal pathogen in a tropical amphibian assemblage. Some amphibian host species are recovering, but the pathogen is still present and is as pathogenic today as it was almost a decade ago. In addition, some species have defenses that are more effective now than they were before the epidemic. These results suggest that host recoveries are not caused by pathogen attenuation and may be due to shifts in host responses. Our findings provide insights into the mechanisms underlying disease transitions, which are increasingly important to understand in an era of emerging infectious diseases and unprecedented global pandemics.

Unlocking the story in the swab: A new genotyping assay for the amphibian chytrid fungus Batrachochytrium dendrobatidis.

One of the most devastating emerging pathogens of wildlife is the chytrid fungus, Batrachochytrium dendrobatidis (Bd), which affects hundreds of amphibian species around the world. Genomic data from pure Bd cultures have advanced our understanding of Bd phylogenetics, genomic architecture and mechanisms of virulence. However, pure cultures are laborious to obtain and whole-genome sequencing is comparatively expensive, so relatively few isolates have been genetically characterized. Thus, we still know little about the genetic diversity of Bd in natural systems. The most common noninvasive method of sampling Bd from natural populations is to swab amphibian skin. Hundreds of thousands of swabs have been collected from amphibians around the world, but Bd DNA collected via swabs is often low in quality and/or quantity. In this study, we developed a custom Bd genotyping assay using the Fluidigm Access Array platform to amplify 192 carefully selected regions of the Bd genome. We obtained robust sequence data for pure Bd cultures and field-collected skin swabs. This new assay has the power to accurately discriminate among the major Bd clades, recovering the basic tree topology previously revealed using whole-genome data. Additionally, we established a critical value for initial Bd load for swab samples (150 Bd genomic equivalents) above which our assay performs well. By leveraging advances in microfluidic multiplex PCR technology and the globally distributed resource of amphibian swab samples, noninvasive skin swabs can now be used to address critical spatial and temporal questions about Bd and its effects on declining amphibian populations.