Population genetics of the frog-killing fungus Batrachochytrium dendrobatidis

Sequence capture identifies fastidious chytrid fungi directly from host tissue

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

Predicted Batrachochytrium dendrobatidis infection sites in Guyana, Suriname, and French Guiana using the species distribution model maxent

The fungal pathogen Batrachochytrium dendrobatidis (Bd) which causes that amphibian disease chytridiomycosis is expanding its worldwide range from an Asian origin, infecting amphibians in a growing number of countries. Modelling the potential range of this amphibian pathogen using environmental variables and presence data could advance our understanding of at-risk areas and species in locations with limited surveillance to date. We used a species distribution model to assess Bd habitat suitability in the three Guiana's (Guyana, Suriname, and French Guiana) in South America. The model output showed that all three countries have substantial areas where Bd could grow and proliferate, and maximum temperature of the warmest month was the top predictor of suitable Bd habitat, inversely correlated with modeled Bd occurrence. Predicted Bd infection areas in Guyana and French Guiana were large and localized whereas possible sites in Suriname were more scattered throughout the country. The areas projected as potential suitable in Suriname were mostly high elevation regions. These results could help inform efficiencies for development of a proactive monitoring program that could alert managers of novel Bd outbreaks for focused mitigation actions to forestall the spread of this amphibian disease.

Panzootic chytrid fungus exploits diverse amphibian host environments through plastic infection strategies

While some pathogens are limited to single species, others can colonize many hosts, likely contributing to the emergence of novel disease outbreaks. Despite this biodiversity threat, traits associated with host niche expansions are not well understood in multi-host pathogens. Here, we aimed to uncover functional machinery driving multi-host invasion by focusing on Batrachochytrium dendrobatidis (Bd), a pathogen that infects the skin of hundreds of amphibians worldwide. We performed a meta-analysis of Bd gene expression using data from published infection experiments and newly generated profiles. We analyzed Bd transcriptomic landscapes across the skin of 14 host species, reconstructed Bd isolates phylogenetic relationships, and inferred the origin and evolutionary history of differentially expressed genes under a phylogenetic framework comprising other 12 zoosporic fungi. Bd displayed plastic infection strategies when challenged by hosts with different disease susceptibility. Our analyses identified sets of differentially expressed genes under host environments with similar infection outcome. We stressed nutritional immunity and gene silencing as important processes required to overcome challenging skin environments in less susceptible hosts. Overall, Bd genes expressed during amphibian skin exploitation have arisen mainly via gene duplications with great family expansions, increasing the gene copy events previously described for this fungal species. Finally, we provide a comprehensive gene dataset that can be used to further examine eco-evolutionary hypotheses for this host-pathogen system. Our study supports the idea that host environments exert contrasting selective pressures, such that gene expression plasticity could be one of the evolutionary keys leading to the success of multi-host pathogens.

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.

Social Behavior, Community Composition, Pathogen Strain, and Host Symbionts Influence Fungal Disease Dynamics in Salamanders

The emerging fungal pathogen, Batrachochytrium dendrobatidis (Bd), which can cause a fatal disease called chytridiomycosis, is implicated in the collapse of hundreds of host amphibian species. We describe chytridiomycosis dynamics in two co-occurring terrestrial salamander species, the Santa Lucia Mountains slender salamander, Batrachoseps luciae, and the arboreal salamander, Aneides lugubris. We (1) conduct a retrospective Bd-infection survey of specimens collected over the last century, (2) estimate present-day Bd infections in wild populations, (3) use generalized linear models (GLM) to identify biotic and abiotic correlates of infection risk, (4) investigate susceptibility of hosts exposed to Bd in laboratory trials, and (5) examine the ability of host skin bacteria to inhibit Bd in culture. Our historical survey of 2,866 specimens revealed that for most of the early 20th century (~1920–1969), Bd was not detected in either species. By the 1990s the proportion of infected specimens was 29 and 17% (B. luciae and A. lugubris, respectively), and in the 2010s it was 10 and 17%. This was similar to the number of infected samples from contemporary populations (2014–2015) at 10 and 18%. We found that both hosts experience signs of chytridiomycosis and suffered high Bd-caused mortality (88 and 71% for B. luciae and A. lugubris, respectively). Our GLM revealed that Bd-infection probability was positively correlated with intraspecific group size and proximity to heterospecifics but not to abiotic factors such as precipitation, minimum temperature, maximum temperature, mean temperature, and elevation, or to the size of the hosts. Finally, we found that both host species contain symbiotic skin-bacteria that inhibit growth of Bd in laboratory trials. Our results provide new evidence consistent with other studies showing a relatively recent Bd invasion of amphibian host populations in western North America and suggest that the spread of the pathogen may be enabled both through conspecific and heterospecific host interactions. Our results suggest that wildlife disease studies should assess host-pathogen dynamics that consider the interactions and effects of multiple hosts, as well as the historical context of pathogen invasion, establishment, and epizootic to enzootic transitions to better understand and predict disease dynamics.

Tracking, Synthesizing, and Sharing Global Batrachochytrium Data at AmphibianDisease.org

Emerging infectious diseases have been especially devastating to amphibians, the most endangered class of vertebrates. For amphibians, the greatest disease threat is chytridiomycosis, caused by one of two chytridiomycete fungal pathogens Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal). Research over the last two decades has shown that susceptibility to this disease varies greatly with respect to a suite of host and pathogen factors such as phylogeny, geography (including abiotic factors), host community composition, and historical exposure to pathogens; yet, despite a growing body of research, a comprehensive understanding of global chytridiomycosis incidence remains elusive. In a large collaborative effort, Bd-Maps was launched in 2007 to increase multidisciplinary investigations and understanding using compiled global Bd occurrence data (Bsal was not discovered until 2013). As its database functions aged and became unsustainable, we sought to address critical needs utilizing new technologies to meet the challenges of aggregating data to facilitate research on both Bd and Bsal. Here, we introduce an advanced central online repository to archive, aggregate, and share Bd and Bsal data collected from around the world. The Amphibian Disease Portal (https://amphibiandisease.org) addresses several critical community needs while also helping to build basic biological knowledge of chytridiomycosis. This portal could be useful for other amphibian diseases and could also be replicated for uses with other wildlife diseases. We show how the Amphibian Disease Portal provides: (1) a new repository for the legacy Bd-Maps data; (2) a repository for sample-level data to archive datasets and host published data with permanent DOIs; (3) a flexible framework to adapt to advances in field, laboratory, and informatics technologies; and (4) a global aggregation of Bd and Bsal infection data to enable and accelerate research and conservation. The new framework for this project is built using biodiversity informatics best practices and metadata standards to ensure scientific reproducibility and linkages across other biological and biodiversity repositories.

Distribution and Genetic Diversity of the Amphibian Chytrid in Japan

While research on frog chytrid fungus Batrachochytrium dendrobatidis (Bd), an infectious disease that threatens amphibian diversity, continues to advance worldwide, little progress has been made in Japan since around 2010. The reason for this is, which we pointed out in 2009, that the origin of frog chytrid fungus may be in the East Asian region, including Japan based on the Bd ITS-DNA variation, and as few cases of mass mortality caused by this fungus have been observed in wild amphibian populations in Japan, the interest of the Japanese government and the general public in Bd has waned. However, we believe that organizing the data obtained so far in Japan and distributing the status of frog chytrid fungus in Japan to the world will provide useful insight for future risk management of this pathogen. We collected more than 5500 swab samples from wild amphibians throughout Japan from 2009 to 2010. Then, we investigated the infection status using the Nested-PCR method. We sequenced the obtained DNA samples and constructed a maximum-parsimony (MP) tree to clarify the phylogenetic diversity of Bd. We detected Bd infection in 11 (nine native and two alien) amphibian species in Japan and obtained 44 haplotypes of Bd ITS-DNA. The MP tree showed a high diversity of Bd strains in Japan, suggesting that some strains belong to Bd-GPL and Bd-Brazil. Except for local populations of the Japanese giant salamanders Andrias japonicus in Honshu Island and the sword tail newts Cynops ensicauda in Okinawa Island, the Bd infection prevalence in native amphibian species was very low. The alien bullfrog Aquarana catesbeiana had high Bd infection rates in all areas where they were sampled. No Bd infection was detected in other native amphibians in the areas where giant salamanders, sword tail newts, and bullfrogs were collected, suggesting that many native amphibians are resistant to Bd infection. The sword tail newt of Okinawa Island had both the highest infectious incidence and greatest number of haplotypes. The giant salamanders also showed relatively high infection prevalence, but the infected strains were limited to those specific to this species. These two Caudata species are endemic to a limited area of Japan, and it was thought that they may have been refugia for Bd, which had been distributed in Japan Islands for a long time.

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.

Low MSP-1 haplotype diversity in the West Palearctic population of the avian malaria parasite Plasmodium relictum

Background

Although avian Plasmodium species are widespread and common across the globe, limited data exist on how genetically variable their populations are. Here, the hypothesis that the avian blood parasite Plasmodium relictum exhibits very low genetic diversity in its Western Palearctic transmission area (from Morocco to Sweden in the north and Transcaucasia in the east) was tested.

Methods

The genetic diversity of Plasmodium relictum was investigated by sequencing a portion (block 14) of the fast-evolving merozoite surface protein 1 (MSP1) gene in 75 different P. relictum infections from 36 host species. Furthermore, the full-length MSP1 sequences representing the common block 14 allele was sequenced in order to investigate if additional variation could be found outside block 14.

Results

The majority (72 of 75) of the sequenced infections shared the same MSP1 allele. This common allele has previously been found to be the dominant allele transmitted in Europe.

Conclusion

The results corroborate earlier findings derived from a limited dataset that the globally transmitted malaria parasite P. relictum exhibits very low genetic diversity in its Western Palearctic transmission area. This is likely the result of a recent introduction event or a selective sweep.

Chytridiomycosis in Asian Amphibians, a Global Resource for Batrachochytrium dendrobatidis (Bd) Research

Chytridiomycosis is an emerging infectious disease affecting amphibians globally and it is caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd). Chytridiomycosis has caused dramatic declines and even extinctions in wild amphibian populations in Europe, Australia, Central and North America. Spanning over two and a half decades, extensive research has led to discovery of epizootic and enzootic lineages of this pathogen. However, the Bd–amphibian system had garnered less attention in Asia until recently when an ancestral Bd lineage was identified in the Korean peninsula. Amphibians co-exist with the pathogen in Asia, only sub-lethal effects have been documented on hosts. Such regions are ‘coldspots’ of infection and are an important resource to understand the dynamics between the enzootic pathogen—Bd and its obligate host—amphibians. Insights into the biology of infection have provided new knowledge on the multi-faceted interaction of Bd in a hyperdiverse Asian amphibian community. We present the findings and highlight the knowledge gap that exists, and propose the ways to bridge them. We emphasize that chytridiomycosis in Asia is an important wildlife disease and it needs focussed research, as it is a dynamic front of pathogen diversity and virulence.

HIGHER INFECTION PREVALENCE IN AMPHIBIANS INHABITING HUMAN-MADE COMPARED TO NATURAL WETLANDS

It is unclear how suitable human-made wetlands are for supporting wildlife and how they impact wildlife disease risk. Natural wetlands (those that were created without human actions) can support more diverse and resilient communities that are at lower risk of disease outbreaks. We compared frog community composition and infection with the pathogenic fungus Batrachochytrium dendrobatidis (Bd) between human-made and natural wetlands in Tippecanoe County, Indiana, US. We conducted visual encounter surveys of frog communities and quantified Bd infection prevalence at four natural and five human-made wetlands. Water parameters associated with human practices (e.g., pH, salinity) and surrounding land use were also compared across sites. We found higher Bd infection prevalence at human-made sites than at natural sites, with monthly differences showing highest infection in spring and fall, and decreasing infection with increasing water temperature. However, we found no differences between human-made and natural sites regarding amphibian community composition, water quality, or surrounding land use. Further, we found frog density increased with distance to nearest roads among both human-made and natural sites. These findings might suggest that human-made wetlands can support frog communities similar to natural wetlands, but pose a greater risk of Bd infection.

Genotyping and differential bacterial inhibition of Batrachochytrium dendrobatidis in threatened amphibians in Costa Rica

Amphibians have declined around the world in recent years, in parallel with the emergence of an epidermal disease called chytridiomycosis, caused by the chytrid fungus Batrachochytrium dendrobatidis (Bd). This disease has been associated with mass mortality in amphibians worldwide, including in Costa Rica, and Bd is considered an important contributor to the disappearance of this group of vertebrates. While many species are susceptible to the disease, others show tolerance and manage to survive infection with the pathogen. We evaluated the pathogen Bd circulating in Costa Rica and the capacity of amphibian skin bacteria to inhibit the growth of the pathogen in vitro. We isolated and characterized - genetically and morphologically - several Bd isolates from areas with declining populations of amphibians. We determined that the circulating chytrid fungus in Costa Rica belongs to the virulent strain Bd-GPL-2, which has been related to massive amphibian deaths worldwide; however, the isolates obtained showed genetic and morphological variation. Furthermore, we isolated epidermal bacteria from 12 amphibian species of surviving populations, some in danger of extinction, and evaluated their inhibitory activity against the collection of chytrid isolates. Through bioassays we confirmed the presence of chytrid-inhibitory bacterial genera in Costa Rican amphibians. However, we observed that the inhibition varied between different isolates of the same bacterial genus, and each bacterial isolation inhibited fungal isolation differently. In total, 14 bacterial isolates belonging to the genera Stenotrophomonas, Streptomyces, Enterobacter, Pseudomonas and Klebsiella showed inhibitory activity against all Bd isolates. Given the observed variation both in the pathogen and in the bacterial inhibition capacity, it is highly relevant to include local isolates and to consider the origin of the microorganisms when performing in vivo infection tests aimed at developing and implementing mitigation strategies for chytridiomycosis.

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

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

Isolation and maintenance of Batrachochytrium salamandrivorans cultures

Discovered in 2013, the chytrid fungus Batrachochytrium salamandrivorans (Bsal) is an emerging amphibian pathogen that causes ulcerative skin lesions and multifocal erosion. A closely related pathogen, B. dendrobatidis (Bd), has devastated amphibian populations worldwide, suggesting that Bsal poses a significant threat to global salamander biodiversity. To expedite research into this emerging threat, we seek to standardize protocols across the field so that results of laboratory studies are reproducible and comparable. We have collated data and experience from multiple labs to standardize culturing practices of Bsal. Here we outline common culture practices including a medium for standardized Bsal growth, standard culturing protocols, and a method for isolating Bsal from infected tissue.

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.

Fungal clones win the battle, but recombination wins the war

Clonal reproduction is common in fungi and fungal-like organisms during epidemics and invasion events. The success of clonal fungi shaped systems for their classification and some pathogens are tacitly treated as asexual. We argue that genetic recombination driven by sexual reproduction must be a starting hypothesis when dealing with fungi for two reasons: (1) Clones eventually crash because they lack adaptability; and (2) fungi find a way to exchange genetic material through recombination, whether sexual, parasexual, or hybridisation. Successful clones may prevail over space and time, but they are the product of recombination and the next successful clone will inevitably appear. Fungal pathogen populations are dynamic rather than static, and they need genetic recombination to adapt to a changing environment.

Widespread chytrid infection across frogs in the Peruvian Amazon suggests critical role for low elevation in pathogen spread and persistence

Outbreaks of emerging infectious diseases are becoming more frequent as climate changes wildlife communities at unprecedented rates, driving population declines and raising concerns for species conservation. One critical disease is the global pandemic of chytridiomycosis in frogs, which can be caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd). Although there is clear evidence for Bd-induced mortality across high-elevation frog communities, little attention is given to the role of lowlands in Bd’s persistence and spread because low elevations are assumed to be too warm to harbor significant levels of Bd. Here, we report widespread Bd infection across 80 frog species from three sites in the lowland Peruvian Amazon, an area with no documented Bd-related amphibian declines. Despite observing no clinical signs of infection in the field, we found that 24–46% of individuals were infected per site (up to ≈105,000 zoospore equivalents per frog) by three Bd strains from the global pandemic lineage (Bd-GPL). We also found collection site and seasonal effects to be only weak predictors of Bd prevalence and load, with lower elevation and drier habitats marginally decreasing both prevalence and load. We found no further effect of host phylogeny, ecotype, or body size. Our results showing high and widespread prevalence across a lowland tropical ecosystem contradict the expectations based on the global pattern of pathogenicity of Bd that is largely restricted to higher elevations and colder temperatures. These findings imply that the lowlands may play a critical role in the spread and persistence of Bd over time and space.

Evolutionary rescue in a host-pathogen system results in coexistence not clearance

The evolutionary rescue of host populations may prevent extinction from novel pathogens. However, the conditions that facilitate rapid evolution of hosts, in particular the population variation in host susceptibility, and the effects of host evolution in response to pathogens on population outcomes remain largely unknown. We constructed an individual-based model to determine the relationships between genetic variation in host susceptibility and population persistence in an amphibian-fungal pathogen (Batrachochytrium dendrobatidis) system. We found that host populations can rapidly evolve reduced susceptibility to a novel pathogen and that this rapid evolution led to a 71-fold increase in the likelihood of host-pathogen coexistence. However, the increased rates of coexistence came at a cost to host populations; fewer populations cleared infection, population sizes were depressed, and neutral genetic diversity was lost. Larger adult host population sizes and greater adaptive genetic variation prior to the onset of pathogen introduction led to substantially reduced rates of extinction, suggesting that populations with these characteristics should be prioritized for conservation when species are threatened by novel infectious diseases.

Signatures of historical selection on MHC reveal different selection patterns in the moor frog (Rana arvalis)

MHC genes are key components in disease resistance and an excellent system for studying selection acting on genetic variation in natural populations. Current patterns of variation in MHC genes are likely to be influenced by past and ongoing selection as well as demographic fluctuations in population size such as those imposed by post-glacial recolonization processes. Here, we investigated signatures of historical selection and demography on an MHC class II gene in 12 moor frog populations along a 1700-km latitudinal gradient. Sequences were obtained from 207 individuals and consecutively assigned into two different clusters (northern and southern clusters, respectively) in concordance with a previously described dual post-glacial colonization route. Selection analyses comparing the relative rates of non-synonymous to synonymous substitutions (dN/dS) suggested evidence of different selection patterns in the northern and the southern clusters, with divergent selection prevailing in the south but uniform positive selection predominating in the north. Also, models of codon evolution revealed considerable differences in the strength of selection: The southern cluster appeared to be under strong selection while the northern cluster showed moderate signs of selection. Our results indicate that the MHC alleles in the north diverged from southern MHC alleles as a result of differential selection patterns.

Epidemic and endemic pathogen dynamics correspond to distinct host population microbiomes at a landscape scale

Infectious diseases have serious impacts on human and wildlife populations, but the effects of a disease can vary, even among individuals or populations of the same host species. Identifying the reasons for this variation is key to understanding disease dynamics and mitigating infectious disease impacts, but disentangling cause and correlation during natural outbreaks is extremely challenging. This study aims to understand associations between symbiotic bacterial communities and an infectious disease, and examines multiple host populations before or after pathogen invasion to infer likely causal links. The results show that symbiotic bacteria are linked to fundamentally different outcomes of pathogen infection: host-pathogen coexistence (endemic infection) or host population extirpation (epidemic infection). Diversity and composition of skin-associated bacteria differed between populations of the frog, Rana sierrae, that coexist with or were extirpated by the fungal pathogen, Batrachochytrium dendrobatidis (Bd). Data from multiple populations sampled before or after pathogen invasion were used to infer cause and effect in the relationship between the fungal pathogen and symbiotic bacteria. Among host populations, variation in the composition of the skin microbiome was most strongly predicted by pathogen infection severity, even in analyses where the outcome of infection did not vary. This result suggests that pathogen infection shapes variation in the skin microbiome across host populations that coexist with or are driven to extirpation by the pathogen. By contrast, microbiome richness was largely unaffected by pathogen infection intensity, but was strongly predicted by geographical region of the host population, indicating the importance of environmental or host genetic factors in shaping microbiome richness. Thus, while both richness and composition of the microbiome differed between endemic and epidemic host populations, the underlying causes are most likely different: pathogen infection appears to shape microbiome composition, while microbiome richness was less sensitive to pathogen-induced disturbance. Because higher richness was correlated with host persistence in the presence of Bd, and richness appeared relatively stable to Bd infection, microbiome richness may contribute to disease resistance, although the latter remains to be directly tested.

Probiotic consortia are not uniformly effective against different amphibian chytrid pathogen isolates

Symbiotic bacterial communities can protect their hosts from infection by pathogens. Treatment of wild individuals with protective bacteria (probiotics) isolated from hosts can combat the spread of emerging infectious diseases. However, it is unclear whether candidate probiotic bacteria can offer consistent protection across multiple isolates of globally distributed pathogens. Here, we use the lethal amphibian fungal pathogen Batrachochytrium dendrobatidis to investigate whether probiotic richness (number of bacteria) or genetic distance among consortia members influences broad-scale in vitro inhibitory capabilities of probiotics across multiple isolates of the pathogen. We show that inhibition of multiple pathogen isolates by individual bacteria is rare, with no systematic pattern among bacterial genera in ability to inhibit multiple B. dendrobatidis isolates. Bacterial consortia can offer stronger protection against B. dendrobatidis compared to single strains, and this tended to be more pronounced for consortia containing multiple genera compared with those consisting of bacteria from a single genus (i.e., with lower genetic distance), but critically, this effect was not uniform across all B. dendrobatidis isolates. These novel insights have important implications for the effective design of bacterial probiotics to mitigate emerging infectious diseases.

Dr Jekyll and Mrs Hyde: Risky hybrid sex by amphibian-parasitizing chytrids in the Brazilian Atlantic Forests

In their article in this issue of Molecular Ecology, Jenkinson et al. () and colleagues address a worrying question-how could arguably the most dangerous pathogen known to science, Batrachochytrium dendrobatidis (Bd), become even more virulent? The answer: start having sex. Jenkinson et al. present a case for how the introduction into Brazil of the globally invasive lineage of Bd, BdGPL, has disrupted the relationship between native amphibians and an endemic Bd lineage, BdBrazil. BdBrazil is hypothesized to be native to the Atlantic Forest and so have a long co-evolutionary history with biodiverse Atlantic Forest amphibian community. The authors suggest that this has resulted in a zone of hybrid Bd genotypes which are potentially more likely to cause fatal chytridiomycosis than either parent lineage. The endemic-nonendemic Bd hybrid genotypes described in this study, and the evidence for pathogen translocation via the global amphibian trade presented, highlights the danger of anthropogenic pathogen dispersal. This research emphasizes that biosecurity regulations may have to refocus on lineages within species if we are to mitigate against the danger of new, possibly hypervirulent genotypes of pathogens emerging as phylogeographic barriers are breached.

Host size influences the effects of four isolates of an amphibian chytrid fungus

Understanding factors that influence host–pathogen interactions is key to predicting outbreaks in natural systems experiencing environmental change. Many amphibian population declines have been attributed to an amphibian chytrid fungus, Batrachochytrium dendrobatidis (Bd). While this fungus is widespread, not all Bd‐positive populations have been associated with declines, which could be attributed to differences in pathogen virulence or host susceptibility. In a laboratory experiment, we examined the effects of Bd isolate origin, two from areas with Bd‐associated amphibian population declines (El Copé, Panama, and California, USA) and two from areas without Bd‐related population declines (Ohio and Maine, USA), on the terrestrial growth and survival of American toad (Anaxyrus americanus) metamorphs reared in larval environments with low or high intraspecific density. We predicted that (1) Bd isolates from areas experiencing declines would have greater negative effects than Bd isolates from areas without declines, and (2) across all isolates, growth and survival of smaller toads from high‐density larval conditions would be reduced by Bd exposure compared to larger toads from low‐density larval conditions. Our results showed that terrestrial survival was reduced for smaller toads exposed to Bd with variation in the response to different isolates, suggesting that smaller size increased susceptibility to Bd. Toads exposed to Bd gained less mass, which varied by isolate. Bd isolates from areas with population declines, however, did not have more negative effects than isolates from areas without recorded declines. Most strikingly, our study supports that host condition, measured by size, can be indicative of the negative effects of Bd exposure. Further, Bd isolates’ impact may vary in ways not predictable from place of origin or occurrence of disease‐related population declines. This research suggests that amphibian populations outside of areas experiencing Bd‐associated declines could be impacted by this pathogen and that the size of individuals could influence the magnitude of Bd's impact.

Contribution of Multiple Inter-Kingdom Horizontal Gene Transfers to Evolution and Adaptation of Amphibian-Killing Chytrid, Batrachochytrium dendrobatidis

Amphibian populations are experiencing catastrophic declines driven by the fungal pathogen Batrachochytrium dendrobatidis (Bd). Although horizontal gene transfer (HGT) facilitates the evolution and adaptation in many fungi by conferring novel function genes to the recipient fungi, inter-kingdom HGT in Bd remains largely unexplored. In this study, our investigation detects 19 bacterial genes transferred to Bd, including metallo-beta-lactamase and arsenate reductase that play important roles in the resistance to antibiotics and arsenates. Moreover, three probable HGT gene families in Bd are from plants and one gene family coding the ankyrin repeat-containing protein appears to come from oomycetes. The observed multi-copy gene families associated with HGT are probably due to the independent transfer events or gene duplications. Five HGT genes with extracellular locations may relate to infection, and some other genes may participate in a variety of metabolic pathways, and in doing so add important metabolic traits to the recipient. The evolutionary analysis indicates that all the transferred genes evolved under purifying selection, suggesting that their functions in Bd are similar to those of the donors. Collectively, our results indicate that HGT from diverse donors may be an important evolutionary driver of Bd, and improve its adaptations for infecting and colonizing host amphibians.

Heritability of Batrachochytrium dendrobatidis burden and its genetic correlation with development time in a population of Common toad (Bufo spinosus)

Despite the important threat that emerging pathogens pose for the conservation of biodiversity as well as human health, very little is known about the adaptive potential of host species to withstand infections. We studied the quantitative genetic architecture responsible for the burden of the fungal pathogen Batrachochytrium dendrobatidis in a population of common toads in conjunction with other life-history traits (i.e., body size and development rate) that may be affected by common selective pressures. We found a significant heritable component that is associated with fungal burden, which may allow for local adaptation to this pathogen to proceed. In addition, the high genetic correlation found between fungal burden and development time suggests that both traits have to be taken into account in order to assess the adaptive response of host populations to this emerging pathogen.

Insights From Genomics Into Spatial and Temporal Variation in Batrachochytrium dendrobatidis

Advances in genetics and genomics have provided new tools for the study of emerging infectious diseases. Researchers can now move quickly from simple hypotheses to complex explanations for pathogen origin, spread, and mechanisms of virulence. Here we focus on the application of genomics to understanding the biology of the fungal pathogen Batrachochytrium dendrobatidis (Bd), a novel and deadly pathogen of amphibians. We provide a brief history of the system, then focus on key insights into Bd variation garnered from genomics approaches, and finally, highlight new frontiers for future discoveries. Genomic tools have revealed unexpected complexity and variation in the Bd system suggesting that the history and biology of emerging pathogens may not be as simple as they initially seem.

Comparative transcriptome analyses of seven anurans reveal functions and adaptations of amphibian skin

Animal skin, which is the tissue that directly contacts the external surroundings, has evolved diverse functions to adapt to various environments. Amphibians represent the transitional taxon from aquatic to terrestrial life. Exploring the molecular basis of their skin function and adaptation is important to understand the survival and evolutionary mechanisms of vertebrates. However, comprehensive studies on the molecular mechanisms of skin functions in amphibians are scarce. In this study, we sequenced the skin transcriptomes of seven anurans belonging to three families and compared the similarities and differences in expressed genes and proteins. Unigenes and pathways related to basic biological processes and special functions, such as defense, immunity, and respiration, were enriched in functional annotations. A total of 108 antimicrobial peptides were identified. The highly expressed genes were similar in species of the same family but were different among families. Additionally, the positively selected orthologous groups were involved in biosynthesis, metabolism, immunity, and defense processes. This study is the first to generate extensive transcriptome data for the skin of seven anurans and provides unigenes and pathway candidates for further studies on amphibian skin function and adaptation.

pophelper: an R package and web app to analyse and visualize population structure

The pophelper r package and web app are software tools to aid in population structure analyses. They can be used for the analyses and visualization of output generated from population assignment programs such as admixture, structure and tess. Some of the functions include parsing output run files to tabulate data, estimating K using the Evanno method, generating files for clumpp and functionality to create barplots. These functions can be streamlined into standard r analysis workflows. The latest version of the package is available on github (https://github.com/royfrancis/pophelper). An interactive web version of the pophelper package is available which covers the same functionalities as the r package version with features such as interactive plots, cluster alignment during plotting, sorting individuals and ordering of population groups. The interactive version is available at http://pophelper.com/.

The roles of sexual and asexual reproduction in the origin and dissemination of strains causing fungal infectious disease outbreaks

Sexual reproduction commonly refers to the reproductive process in which genomes from two sources are combined into a single cell through mating and then the zygote genomes are partitioned to progeny cells through meiosis. Reproduction in the absence of mating and meiosis is referred to as asexual or clonal reproduction. One major advantage of sexual reproduction is that it generates genetic variation among progeny which may allow for faster adaptation of the population to novel and/or stressful environments. However, adaptation to stressful or new environments can still occur through mutation, in the absence of sex. In this review, we analyzed the relative contributions of sexual and asexual reproduction in the origin and spread of strains causing fungal infectious diseases outbreaks. The necessity of sex and the ability of asexual fungi to initiate outbreaks are discussed. We propose a framework that relates the modes of reproduction to the origin and propagation of fungal disease outbreaks. Our analyses suggest that both sexual and asexual reproduction can play critical roles in the origin of outbreak strains and that the rapid spread of outbreak strains is often accomplished through asexual expansion.

Disentangling host, pathogen, and environmental determinants of a recently emerged wildlife disease: lessons from the first 15 years of amphibian chytridiomycosis research

The amphibian fungal disease chytridiomycosis, which affects species across all continents, recently emerged as one of the greatest threats to biodiversity. Yet, many aspects of the basic biology and epidemiology of the pathogen, Batrachochytrium dendrobatidis (Bd), are still unknown, such as when and from where did Bd emerge and what is its true ecological niche? Here, we review the ecology and evolution of Bd in the Americas and highlight controversies that make this disease so enigmatic. We explore factors associated with variance in severity of epizootics focusing on the disease triangle of host susceptibility, pathogen virulence, and environment. Reevaluating the causes of the panzootic is timely given the wealth of data on Bd prevalence across hosts and communities and the recent discoveries suggesting co-evolutionary potential of hosts and Bd. We generate a new species distribution model for Bd in the Americas based on over 30,000 records and suggest a novel future research agenda. Instead of focusing on pathogen "hot spots," we need to identify pathogen "cold spots" so that we can better understand what limits the pathogen's distribution. Finally, we introduce the concept of "the Ghost of Epizootics Past" to discuss expected patterns in postepizootic host communities.

DNA Extraction Method Affects the Detection of a Fungal Pathogen in Formalin-Fixed Specimens Using qPCR

Museum collections provide indispensable repositories for obtaining information about the historical presence of disease in wildlife populations. The pathogenic amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd) has played a significant role in global amphibian declines, and examining preserved specimens for Bd can improve our understanding of its emergence and spread. Quantitative PCR (qPCR) enables Bd detection with minimal disturbance to amphibian skin and is significantly more sensitive to detecting Bd than histology; therefore, developing effective qPCR methodologies for detecting Bd DNA in formalin-fixed specimens can provide an efficient and effective approach to examining historical Bd emergence and prevalence. Techniques for detecting Bd in museum specimens have not been evaluated for their effectiveness in control specimens that mimic the conditions of animals most likely to be encountered in museums, including those with low pathogen loads. We used American bullfrogs (Lithobates catesbeianus) of known infection status to evaluate the success of qPCR to detect Bd in formalin-fixed specimens after three years of ethanol storage. Our objectives were to compare the most commonly used DNA extraction method for Bd (PrepMan, PM) to Macherey-Nagel DNA FFPE (MN), test optimizations for Bd detection with PM, and provide recommendations for maximizing Bd detection. We found that successful detection is relatively high (80-90%) when Bd loads before formalin fixation are high, regardless of the extraction method used; however, at lower infection levels, detection probabilities were significantly reduced. The MN DNA extraction method increased Bd detection by as much as 50% at moderate infection levels. Our results indicate that, for animals characterized by lower pathogen loads (i.e., those most commonly encountered in museum collections), current methods may underestimate the proportion of Bd-infected amphibians. Those extracting DNA from archived museum specimens should ensure that the techniques they are using are known to provide high-quality throughput DNA for later analysis.

Clonal reproduction shapes evolution in the lizard malaria parasite Plasmodium floridense

The preponderant clonal evolution hypothesis (PCE) predicts that frequent clonal reproduction (sex between two clones) in many pathogens capable of sexual recombination results in strong linkage disequilibrium and the presence of discrete genetic subdivisions characterized by occasional gene flow. We expand on the PCE and predict that higher rates of clonal reproduction will result in: (1) morphologically cryptic species that exhibit (2) low within-species variation and (3) recent between-species divergence. We tested these predictions in the Caribbean lizard malaria parasite Plasmodium floridense using 63 single-infection samples in lizards collected from across the parasite's range, and sequenced them at two mitochondrial, one apicoplast, and five nuclear genes. We identified 11 provisionally cryptic species within P. floridense, each of which exhibits low intraspecific variation and recent divergence times between species (some diverged approximately 110,000 years ago). Our results are consistent with the hypothesis that clonal reproduction can profoundly affect diversification of species capable of sexual recombination, and suggest that clonal reproduction may have led to a large number of unrecognized pathogen species. The factors that may influence the rates of clonal reproduction among pathogens are unclear, and we discuss how prevalence and virulence may relate to clonal reproduction.

Chytrid fungus acts as a generalist pathogen infecting species-rich amphibian families in Brazilian rainforests

The fungus Batrachochytrium dendrobatidis (Bd) is among the main causes of declines in amphibian populations. This fungus is considered a generalist pathogen because it infects several species and spreads rapidly in the wild. To date, Bd has been detected in more than 100 anuran species in Brazil, mostly in the southern portion of the Atlantic forest. Here, we report survey data from some poorly explored regions; these data considerably extend current information on the distribution of Bd in the northern Atlantic forest region. In addition, we tested the hypothesis that Bd is a generalist pathogen in this biome. We also report the first positive record for Bd in an anuran caught in the wild in Amazonia. In total, we screened 90 individuals (from 27 species), of which 39 individuals (from 22 species) were Bd-positive. All samples collected in Bahia (2 individuals), Pernambuco (3 individuals), Pará (1 individual), and Minas Gerais (1 individual) showed positive results for Bd. We found a positive correlation between anuran richness per family and the number of infected species in the Atlantic forest, supporting previous observations that Bd lacks strong host specificity; of 38% of the anuran species in the Atlantic forest that were tested for Bd infection, 25% showed positive results. The results of our study exemplify the pandemic and widespread nature of Bd infection in amphibians.

Chytrid fungus infections in laboratory and introduced Xenopus laevis populations: assessing the risks for U.K. native amphibians

The chytrid fungus Batrachochytrium dendrobatidis (Bd) is notorious amongst current conservation biology challenges, responsible for mass mortality and extinction of amphibian species. World trade in amphibians is implicated in global dissemination. Exports of South African Xenopus laevis have led to establishment of this invasive species on four continents. Bd naturally infects this host in Africa and now occurs in several introduced populations. However, no previous studies have investigated transfer of infection into co-occurring native amphibian faunas. A survey of 27 U.K. institutions maintaining X. laevis for research showed that most laboratories have low-level infection, a risk for native species if animals are released into the wild. RT-PCR assays showed Bd in two introduced U.K. populations of X. laevis, in Wales and Lincolnshire. Laboratory and field studies demonstrated that infection levels increase with stress, especially low temperature. In the U.K., native amphibians may be exposed to intense transmission in spring when they enter ponds to spawn alongside X. laevis that have cold-elevated Bd infections. Exposure to cross-infection has probably been recurrent since the introduction of X. laevis, >20 years in Lincolnshire and 50 years in Wales. These sites provide an important test for assessing the impact of X. laevis on Bd spread. However, RT-PCR assays on 174 native amphibians (Bufo, Rana, Lissotriton and Triturus spp.), sympatric with the Bd-infected introduced populations, showed no foci of self-sustaining Bd transmission associated with X. laevis. The abundance of these native amphibians suggested no significant negative population-level effect after the decades of co-occurrence.

Unisexual reproduction

Sexual reproduction is ubiquitous throughout the eukaryotic kingdom, but the capacity of pathogenic fungi to undergo sexual reproduction has been a matter of intense debate. Pathogenic fungi maintained a complement of conserved meiotic genes but the populations appeared to be clonally derived. This debate was resolved first with the discovery of an extant sexual cycle and then unisexual reproduction. Unisexual reproduction is a distinct form of homothallism that dispenses with the requirement for an opposite mating type. Pathogenic and nonpathogenic fungi previously thought to be asexual are able to undergo robust unisexual reproduction. We review here recent advances in our understanding of the genetic and molecular basis of unisexual reproduction throughout fungi and the impact of unisex on the ecology and genomic evolution of fungal species.

The pathogen Batrachochytrium dendrobatidis disturbs the frog skin microbiome during a natural epidemic and experimental infection

Symbiotic microbial communities may interact with infectious pathogens sharing a common host. The microbiome may limit pathogen infection or, conversely, an invading pathogen can disturb the microbiome. Documentation of such relationships during naturally occurring disease outbreaks is rare, and identifying causal links from field observations is difficult. This study documented the effects of an amphibian skin pathogen of global conservation concern [the chytrid fungus Batrachochytrium dendrobatidis (Bd)] on the skin-associated bacterial microbiome of the endangered frog, Rana sierrae, using a combination of population surveys and laboratory experiments. We examined covariation of pathogen infection and bacterial microbiome composition in wild frogs, demonstrating a strong and consistent correlation between Bd infection load and bacterial community composition in multiple R. sierrae populations. Despite the correlation between Bd infection load and bacterial community composition, we observed 100% mortality of postmetamorphic frogs during a Bd epizootic, suggesting that the relationship between Bd and bacterial communities was not linked to variation in resistance to mortal disease and that Bd infection altered bacterial communities. In a controlled experiment, Bd infection significantly altered the R. sierrae microbiome, demonstrating a causal relationship. The response of microbial communities to Bd infection was remarkably consistent: Several bacterial taxa showed the same response to Bd infection across multiple field populations and the laboratory experiment, indicating a somewhat predictable interaction between Bd and the microbiome. The laboratory experiment demonstrates that Bd infection causes changes to amphibian skin bacterial communities, whereas the laboratory and field results together strongly support Bd disturbance as a driver of bacterial community change during natural disease dynamics.

Unisexual reproduction reverses Muller's ratchet

Cryptococcus neoformans is a pathogenic basidiomycetous fungus that engages in outcrossing, inbreeding, and selfing forms of unisexual reproduction as well as canonical sexual reproduction between opposite mating types. Long thought to be clonal, >99% of sampled environmental and clinical isolates of C. neoformans are MATα, limiting the frequency of opposite mating-type sexual reproduction. Sexual reproduction allows eukaryotic organisms to exchange genetic information and shuffle their genomes to avoid the irreversible accumulation of deleterious changes that occur in asexual populations, known as Muller's ratchet. We tested whether unisexual reproduction, which dispenses with the requirement for an opposite mating-type partner, is able to purge the genome of deleterious mutations. We report that the unisexual cycle can restore mutant strains of C. neoformans to wild-type genotype and phenotype, including prototrophy and growth rate. Furthermore, the unisexual cycle allows attenuated strains to purge deleterious mutations and produce progeny that are returned to wild-type virulence. Our results show that unisexual populations of C. neoformans are able to avoid Muller's ratchet and loss of fitness through a unisexual reproduction cycle involving α-α cell fusion, nuclear fusion, and meiosis. Similar types of unisexual reproduction may operate in other pathogenic and saprobic eukaryotic taxa.

Batrachochytrium dendrobatidis prevalence and haplotypes in domestic and imported pet amphibians in Japan

The international trade in amphibians is believed to have increased the spread of Batrachochytrium dendrobatidis (Bd), the fungal pathogen responsible for chytridiomycosis, which has caused a rapid decline in amphibian populations worldwide. We surveyed amphibians imported into Japan and those held in captivity for a long period or bred in Japan to clarify the Bd infection status. Samples were taken from 820 individuals of 109 amphibian species between 2008 and 2011 and were analyzed by a nested-PCR assay. Bd prevalence in imported amphibians was 10.3% (58/561), while it was 6.9% (18/259) in those in private collections and commercially bred amphibians in Japan. We identified the genotypes of this fungus using partial DNA sequences of the internal transcribed spacer (ITS) region. Sequencing of PCR products of all 76 Bd-positive samples revealed 11 haplotypes of the Bd ITS region. Haplotype A (DNA Data Bank of Japan accession number AB435211) was found in 90% (52/58) of imported amphibians. The results show that Bd is currently entering Japan via the international trade in exotic amphibians as pets, suggesting that the trade has indeed played a major role in the spread of Bd.