Prevalence and pathogen load estimates for the fungus Batrachochytrium dendrobatidis are impacted by ITS DNA copy number variation

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

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.

Extraction-free detection of amphibian pathogens from water baths

Detecting and quantifying pathogens with quick, cost-efficient and sensitive methods is needed across disease systems for addressing pertinent epidemiological questions. Typical methods rely on extracting DNA from collected samples. Here we develop and test an extraction-free method from water bath samples that is both sensitive and efficient for 2 major amphibian pathogens-Batrachochytrium dendrobatidis and B. salamandrivorans. We tested mock samples with known pathogen quantities as well as comparatively assessed detection from skin swabs and water baths from field sampled amphibians. Quantitative PCR (qPCR) directly on lyophilized water baths was able to reliably detect low loads of 10 and 1 zoospores for both pathogens, and detection rates were greater than those of swabs from field samples. Further concentration of samples did not improve detection, and collection container type did not influence pathogen load estimates. This method of lyophilization (i.e. freeze-drying) followed by direct qPCR offers an effective and efficient tool from detecting amphibian pathogens, which is crucial for surveillance efforts and estimating shedding rates for robust epidemiological understanding of transmission dynamics. Furthermore, water bath samples have multiple functions and can be used to evaluate mucosal function against pathogens and characterize mucosal components. The multifunctionality of water bath samples and reduced monetary costs and time expenditures make this method an optimal tool for amphibian disease research and may also prove to be useful in other wildlife disease systems.

Co-infecting pathogen lineages have additive effects on host bacterial communities

Amphibian skin bacteria may confer protection against the fungus Batrachochytrium dendrobatidis (Bd), but responses of skin bacteria to different Bd lineages are poorly understood. The global panzootic lineage (Bd-GPL) has caused amphibian declines and extinctions globally. However, other lineages are enzootic (Bd-Asia-2/Brazil). Increased contact rates between Bd-GPL and enzootic lineages via globalization pose unknown consequences for host-microbiome-pathogen dynamics. We conducted a laboratory experiment and used 16S rRNA amplicon-sequencing to assess: (i) whether two lineages (Bd-Asia-2/Brazil and Bd-GPL) and their recombinant, in single and mixed infections, differentially affect amphibian skin bacteria; (ii) and the changes associated with the transition to laboratory conditions. We determined no clear differences in bacterial diversity among Bd treatments, despite differences in infection intensity. However, we observed an additive effect of mixed infections on bacterial alpha diversity and a potentially antagonistic interaction between Bd genotypes. Additionally, observed changes in community composition suggest a higher ability of Bd-GPL to alter skin bacteria. Lastly, we observed a drastic reduction in bacterial diversity and a change in community structure in laboratory conditions. We provide evidence for complex interactions between Bd genotypes and amphibian skin bacteria during coinfections, and expand on the implications of experimental conditions in ecological studies.

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.

see also the Perspective by Claudio Azat

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.

Batrachochytrium salamandrivorans ELICITS ACUTE STRESS RESPONSE IN SPOTTED SALAMANDERS BUT NOT INFECTION OR MORTALITY

The emerging fungal pathogen Batrachochytrium salamandrivorans (Bsal) is a major threat to amphibian species worldwide with potential to infect many species if it invades salamander biodiversity hotspots in the Americas. Bsal can cause the disease chytridiomycosis, and it is important to assess the risk of Bsal-induced chytridiomycosis to species in North America. We evaluated the susceptibility to Bsal of the common and widespread spotted salamander, Ambystoma maculatum, across life history stages and monitored the effect of Bsal exposure on growth rate and response of the stress hormone, corticosterone. We conclude that spotted salamanders appear resistant to Bsal because they showed no indication of disease or infection, and experienced minor effects on growth upon exposure. While we focused on a single population for this study, results were consistent across conditions of exposure including high or repeated doses of Bsal, life-stage at exposure, environmental conditions including two temperatures and two substrates, and promoting pathogen infectivity by conditioning Bsal cultures with thyroid hormone. Exposure to high levels of Bsal elicited an acute but not chronic increase in corticosterone in spotted salamanders, and reduced growth. We hypothesize that the early acute increase in corticosterone facilitated mounting an immune response to the pathogen, perhaps through immunoredistribution to the skin, but further study is needed to determine immune responses to Bsal. These results will contribute to development of appropriate Bsal management plans to conserve species at risk of emerging disease.

Population-Level Resistance to Chytridiomycosis is Life-Stage Dependent in an Imperiled Anuran

Amphibian declines caused by chytridiomycosis have been severe, but some susceptible populations have persisted or even recovered. Resistance to the causal agent Batrachochytrium dendrobatidis (Bd) could result from alleles of the adaptive immune system. During metamorphosis, however, immune systems may not be fully functional, implying that an effective immune response to Bd may be life-stage dependent. We evaluated the susceptibility of the relict leopard frog (Rana onca) sourced from two areas where Bd was present or absent, and where the populations appeared to show differences in pathogen resistance. We evaluated whether population-level resistance manifested across life stages using challenge experiments with late-stage tadpoles (Gosner stage 31-38), metamorphs (stage 45-46), and juvenile frogs. We used three different Bd isolates including one from wild R. onca to challenge juvenile frogs and focused on the isolate from R. onca to challenge tadpoles and resulting metamorphs. We found that juveniles from the Bd exposed population were 5.5 times more likely to survive Bd infection and 10 times more likely to clear infections than those from the area without Bd. In contrast, and regardless of the source area, we observed 98% survivorship of tadpoles, but only 19% survivorship of resulting metamorphs following re-exposure. Given the low survivorship of exposed metamorphs in the laboratory, we speculate on how resistance characteristics, whether adaptive or innate, that do not manifest at each life stage could develop in the wild. We suggest that seasonal high temperatures during times when metamorphosis appears common may modulate the effects of the pathogen during this most susceptible life stage.

Patterns of PCR Amplification Artifacts of the Fungal Barcode Marker in a Hybrid Mushroom

The polymerase chain reaction (PCR) is widely used in modern biology and medicine. However, PCR artifacts can complicate the interpretation of PCR-based results. The internal transcribed spacer (ITS) region of the ribosomal RNA gene cluster is the consensus fungal barcode marker and suspected PCR artifacts have been reported in many studies, especially for the analyses of environmental fungal samples. At present, the patterns of PCR artifacts in the whole fungal ITS region (ITS1+5.8S+ITS2) are not known. In this study, we analyzed the error rates of PCR at three template complexity levels using the divergent copies of ITS from the mushroom Agaricus subrufescens. Our results showed that PCR using the Phusion^® High-Fidelity DNA Polymerase has a per nucleotide error rate of about 4 × 10^–6 per replication. Among the detected mutations, transitions were much more frequent than transversions, insertions, and deletions. When divergent alleles were mixed as templates in the same reaction, a significant proportion (∼30%) of recombinant molecules were detected. The in vitro mixed-template results were comparable to those obtained from using the genomic DNA of the original mushroom specimen as template. Our results indicate that caution should be in place when interpreting ITS sequences from individual fungal specimens, especially those containing divergent ITS copies. Similar results could also happen to PCR-based analyses of other multicopy DNA fragments as well as single-copy DNA sequences with divergent alleles in diploid organisms.

Serratia marcescens shapes cutaneous bacterial communities and influences survival of an amphibian host

Ongoing investigations into the interactions between microbial communities and their associated hosts are changing how emerging diseases are perceived and ameliorated. Of the numerous host-microbiome-disease systems of study, the emergence of chytridiomycosis (caused by Batrachochytrium dendrobatidis, hereafter Bd) has been implicated in ongoing declines and extinction events of amphibians worldwide. Interestingly, there has been differential survival among amphibians in resisting Bd infection and subsequent disease. One factor thought to contribute to this resistance is the host-associated cutaneous microbiota. This has raised the possibility of using genetically modified probiotics to restructure the host-associated microbiota for desired anti-fungal outcomes. Here, we use a previously described strain of Serratia marcescens (Sm) for the manipulation of amphibian cutaneous microbiota. Sm was genetically altered to have a dysfunctional pathway for the production of the extracellular metabolite prodigiosin. This genetically altered strain (Δpig) and the functional prodigiosin producing strain (wild-type, WT) were compared for their microbial community and anti-Bd effects both in vitro and in vivo. In vitro, Bd growth was significantly repressed in the presence of prodigiosin. In vivo, the inoculation of both Sm strains was shown to significantly influence amphibian microbiota diversity with the Δpig-Sm treatment showing increasing alpha diversity, and the WT-Sm having no temporal effect on diversity. Differences were also seen in host mortality with Δpig-Sm treatments exhibiting significantly decreased survival probability when compared with WT-Sm in the presence of Bd. These results are an important proof-of-concept for linking the use of genetically modified probiotic bacteria to host microbial community structure and disease outcomes, which in the future may provide a way to ameliorate disease and address critical frontiers in disease and microbial ecology.

Ancestral chytrid pathogen remains hypervirulent following its long coevolution with amphibian hosts

Many amphibian species around the world, except in Asia, suffer morbidity and mortality when infected by the emerging infectious pathogen Batrachochytrium dendrobatidis (Bd). A lineage of the amphibian chytrid fungus isolated from South Korean amphibians (BdAsia-1) is evolutionarily basal to recombinant global pandemic lineages (BdGPL) associated with worldwide amphibian population declines. In Asia, the Bd pathogen and its amphibian hosts have coevolved over 100 years or more. Thus, resilience of Asian amphibian populations to infection might result from attenuated virulence of endemic Bd lineages, evolved immunity to the pathogen or both. We compared susceptibilities of an Australasian amphibian, Litoria caerulea, known to lack resistance to BdGPL, with those of three Korean species, Bufo gargarizans, Bombina orientalis and Hyla japonica, after inoculation with BdAsia-1, BdGPL or a blank solution. Subjects became infected in all experimental treatments but Korean species rapidly cleared themselves of infection, regardless of Bd lineage. They survived with no apparent secondary effects. By contrast, L. caerulea, after infection by either BdAsia-1 or BdGPL, suffered deteriorating body condition and carried progressively higher Bd loads over time. Subsequently, most subjects died. Comparing their effects on L. caerulea, BdAsia-1 induced more rapid disease progression than BdGPL. The results suggest that genomic recombination with other lineages was not necessary for the ancestral Bd lineage to evolve hypervirulence over its long period of coevolution with amphibian hosts. The pathogen's virulence may have driven strong selection for immune responses in endemic Asian amphibian host species.

Asymptomatic infection of the fungal pathogen Batrachochytrium salamandrivorans in captivity

One of the most important factors driving amphibian declines worldwide is the infectious disease, chytridiomycosis. Two fungi have been associated with this disease, Batrachochytrium dendrobatidis and B. salamandrivorans (Bsal). The latter has recently driven Salamandra salamandra populations to extirpation in parts of the Netherlands, and Belgium, and potentially also in Germany. Bsal has been detected in the pet trade, which has been hypothesized to be the pathway by which it reached Europe, and which may continuously contribute to its spread. In the present study, 918 amphibians belonging to 20 captive collections in Germany and Sweden were sampled to explore the extent of Bsal presence in captivity. The fungus was detected by quantitative Polymerase Chain Reaction (qPCR) in ten collections, nine of which lacked clinical symptoms. 23 positives were confirmed by independent processing of duplicate swabs, which were analysed in a separate laboratory, and/or by sequencing ITS and 28 S gene segments. These asymptomatic positives highlight the possibility of Bsal being widespread in captive collections, and is of high conservation concern. This finding may increase the likelihood of the pathogen being introduced from captivity into the wild, and calls for according biosecurity measures. The detection of Bsal-positive alive specimens of the hyper-susceptible fire salamander could indicate the existence of a less aggressive Bsal variant or the importance of environmental conditions for infection progression.