FIRST RECORD OF BATRACHOCHYTRIUM DENDROBATIDIS IN THE NORTHERN NETHERLANDS.

Batrachochytrium dendrobatidis (Bd) infects amphibians and has been linked to the decline of hundreds of anuran amphibians all over the world. In the province of Groningen in the Netherlands, this fungal pathogen was not detected before this study. To determine whether Groningen was Bd-free, we surveyed 12 locations in this province in 2020 and 2021. Samples were then used to quantify the presence of Bd with a qPCR assay. In total, 2 out of 110 (∼0.02%) collected in 2020 and 11 out of 249 samples collected in 2021 tested positive for Bd. Infected amphibians were found in 4 out of the 12 sites, and the prevalence of Bd was estimated at 4% for both years combined. Our study provides the first record of Bd in Groningen, and we hypothesize that Bd is present throughout the Netherlands in regions currently considered "Bd-free." Furthermore, we warn scientists and policymakers to be apprehensive when calling a site free from Bd when sampling is limited or not recent.

Divergent population responses following salamander mass mortalities and declines driven by the emerging pathogen Batrachochytrium salamandrivorans.

Understanding wildlife responses to novel threats is vital in counteracting biodiversity loss. The emerging pathogen Batrachochytrium salamandrivorans (Bsal) causes dramatic declines in European salamander populations, and is considered an imminent threat to global amphibian biodiversity. However, real-life disease outcomes remain largely uncharacterized. We performed a multidisciplinary assessment of the longer-term impacts of Bsal on highly susceptible fire salamander (Salamandra salamandra) populations, by comparing four of the earliest known outbreak sites to uninfected sites. Based on large-scale monitoring efforts, we found population persistence in strongly reduced abundances to over a decade after Bsal invasion, but also the extinction of an initially small-sized population. In turn, we found that host responses varied, and Bsal detection remained low, within surviving populations. Demographic analyses indicated an ongoing scarcity of large reproductive adults with potential for recruitment failure, while spatial comparisons indicated a population remnant persisting within aberrant habitat. Additionally, we detected no early signs of severe genetic deterioration, yet nor of increased host resistance. Beyond offering additional context to Bsal-driven salamander declines, results highlight how the impacts of emerging hypervirulent pathogens can be unpredictable and vary across different levels of biological complexity, and how limited pathogen detectability after population declines may complicate surveillance efforts.

Compensatory recruitment allows amphibian population persistence in anthropogenic habitats.

Habitat anthropization is a major driver of global biodiversity decline. Although most species are negatively affected, some benefit from anthropogenic habitat modifications by showing intriguing life-history responses. For instance, increased recruitment through higher allocation to reproduction or improved performance during early-life stages could compensate for reduced adult survival, corresponding to "compensatory recruitment". To date, evidence of compensatory recruitment in response to habitat modification is restricted to plants, limiting understanding of its importance as a response to global change. We used the yellow-bellied toad (Bombina variegata), an amphibian occupying a broad range of natural and anthropogenic habitats, as a model species to test for and to quantify compensatory recruitment. Using an exceptional capture-recapture dataset composed of 21,714 individuals from 67 populations across Europe, we showed that adult survival was lower, lifespan was shorter, and actuarial senescence was higher in anthropogenic habitats, especially those affected by intense human activities. Increased recruitment in anthropogenic habitats fully offset reductions in adult survival, with the consequence that population growth rate in both habitat types was similar. Our findings indicate that compensatory recruitment allows toad populations to remain viable in human-dominated habitats and might facilitate the persistence of other animal populations in such environments.

Landscape epidemiology of Batrachochytrium salamandrivorans: reconciling data limitations and conservation urgency.

Starting in 2010, rapid fire salamander (Salamandra salamandra) population declines in northwestern Europe heralded the emergence of Batrachochytrium salamandrivorans (Bsal), a salamander-pathogenic chytrid fungus. Bsal poses an imminent threat to global salamander diversity owing to its wide host range, high pathogenicity, and long-term persistence in ecosystems. While there is a pressing need to develop further research and conservation actions, data limitations inherent to recent pathogen emergence obscure necessary insights into Bsal disease ecology. Here, we use a hierarchical modeling framework to describe Bsal landscape epidemiology of outbreak sites in light of these methodological challenges. Using model selection and machine learning, we find that Bsal presence is associated with humid and relatively cool, stable climates. Outbreaks are generally located in areas characterized by low landscape heterogeneity and low steepness of slope. We further find an association between Bsal presence and high trail density, suggesting that human-mediated spread may increase risk for spillover between populations. We then use distribution modeling to show that favorable conditions occur in lowlands influenced by the North Sea, where increased survey effort is needed to determine how Bsal impacts local newt populations, but also in hill- and mountain ranges in northeastern France and the lower half of Germany. Finally, connectivity analyses suggest that these hill- and mountain ranges may act as stepping stones for further spread southward. Our results provide initial insight into regional environmental conditions underlying Bsal epizootics, present updated invasibility predictions for northwestern Europe, and lead us to discuss a wide variety of potential survey and research actions needed to advance future conservation and mitigation efforts.

Ranavirus genotypes in the Netherlands and their potential association with virulence in water frogs (Pelophylax spp.).

Ranaviruses are pathogenic viruses for poikilothermic vertebrates worldwide. The identification of a common midwife toad virus (CMTV) associated with massive die-offs in water frogs (Pelophylax spp.) in the Netherlands has increased awareness for emerging viruses in amphibians in the country. Complete genome sequencing of 13 ranavirus isolates collected from ten different sites in the period 2011-2016 revealed three CMTV groups present in distinct geographical areas in the Netherlands. Phylogenetic analysis showed that emerging viruses from the northern part of the Netherlands belonged to CMTV-NL group I. Group II and III viruses were derived from the animals located in the center-east and south of the country, and shared a more recent common ancestor to CMTV-amphibian associated ranaviruses reported in China, Italy, Denmark, and Switzerland. Field monitoring revealed differences in water frog host abundance at sites where distinct ranavirus groups occur; with ranavirus-associated deaths, host counts decreasing progressively, and few juveniles found in the north where CMTV-NL group I occurs but not in the south with CMTV-NL group III. Investigation of tandem repeats of coding genes gave no conclusive information about phylo-geographical clustering, while genetic analysis of the genomes revealed truncations in 17 genes across CMTV-NL groups II and III compared to group I. Further studies are needed to elucidate the contribution of these genes as well as environmental variables to explain the observed differences in host abundance.

Post-epizootic salamander persistence in a disease-free refugium suggests poor dispersal ability of Batrachochytrium salamandrivorans.

Lack of disease spill-over between adjacent populations has been associated with habitat fragmentation and the absence of population connectivity. We here present a case which describes the absence of the spill-over of the chytrid fungus Batrachochytrium salamandrivorans (Bsal) between two connected subpopulations of fire salamanders (Salamandra salamandra). Based on neutrally evolving microsatellite loci, both subpopulations were shown to form a single genetic cluster, suggesting a shared origin and/or recent gene flow. Alpine newts (Ichthyosaura alpestris) and fire salamanders were found in the landscape matrix between the two sites, which are also connected by a stream and separated by no obvious physical barriers. Performing a laboratory trial using alpine newts, we confirmed that Bsal is unable to disperse autonomously. Vector-mediated dispersal may have been impeded by a combination of sub-optimal connectivity, limited dispersal ability of infected hosts and a lack of suitable dispersers following the rapid, Bsal-driven collapse of susceptible hosts at the source site. Although the exact cause remains unclear, the aggregate evidence suggests that Bsal may be a poorer disperser than previously hypothesized. The lack of Bsal dispersal between neighbouring salamander populations opens perspectives for disease management and stresses the necessity of implementing biosecurity measures preventing human-mediated spread.

Fragile coexistence of a global chytrid pathogen with amphibian populations is mediated by environment and demography.

Unravelling the multiple interacting drivers of host-pathogen coexistence is crucial in understanding how an apparently stable state of endemism may shift towards an epidemic and lead to biodiversity loss. Here, we investigate the apparent coexistence of the global amphibian pathogen Batrachochytrium dendrobatidis (Bd) with Bombina variegata populations in The Netherlands over a 7-year period. We used a multi-season mark-recapture dataset and assessed potential drivers of coexistence (individual condition, environmental mediation and demographic compensation) at the individual and population levels. We show that even in a situation with a clear cost incurred by endemic Bd, population sizes remain largely stable. Current environmental conditions and an over-dispersed pathogen load probably stabilize disease dynamics, but as higher temperatures increase infection probability, changing environmental conditions, for example a climate-change-driven rise in temperature, could unbalance the current fragile host-pathogen equilibrium. Understanding the proximate mechanisms of such environmental mediation and of site-specific differences in infection dynamics can provide vital information for mitigation actions.

Comparison of photo-matching algorithms commonly used for photographic capture-recapture studies.

Photographic capture-recapture is a valuable tool for obtaining demographic information on wildlife populations due to its noninvasive nature and cost-effectiveness. Recently, several computer-aided photo-matching algorithms have been developed to more efficiently match images of unique individuals in databases with thousands of images. However, the identification accuracy of these algorithms can severely bias estimates of vital rates and population size. Therefore, it is important to understand the performance and limitations of state-of-the-art photo-matching algorithms prior to implementation in capture-recapture studies involving possibly thousands of images. Here, we compared the performance of four photo-matching algorithms; Wild-ID, I3S Pattern+, APHIS, and AmphIdent using multiple amphibian databases of varying image quality. We measured the performance of each algorithm and evaluated the performance in relation to database size and the number of matching images in the database. We found that algorithm performance differed greatly by algorithm and image database, with recognition rates ranging from 100% to 22.6% when limiting the review to the 10 highest ranking images. We found that recognition rate degraded marginally with increased database size and could be improved considerably with a higher number of matching images in the database. In our study, the pixel-based algorithm of AmphIdent exhibited superior recognition rates compared to the other approaches. We recommend carefully evaluating algorithm performance prior to using it to match a complete database. By choosing a suitable matching algorithm, databases of sizes that are unfeasible to match "by eye" can be easily translated to accurate individual capture histories necessary for robust demographic estimates.

Investigation of Amphibian Mortality Events in Wildlife Reveals an On-Going Ranavirus Epidemic in the North of the Netherlands.

In the four years following the first detection of ranavirus (genus Ranavirus, family Iridoviridae) infection in Dutch wildlife in 2010, amphibian mortality events were investigated nationwide to detect, characterize and map ranaviruses in amphibians over time, and to establish the affected host species and the clinico-pathological presentation of the disease in these hosts. The ultimate goal was to obtain more insight into ranavirus disease emergence and ecological risk. In total 155 dead amphibians from 52 sites were submitted between 2011 and 2014, and examined using histopathology, immunohistochemistry, virus isolation and molecular genetic characterization. Ranavirus-associated amphibian mortality events occurred at 18 sites (35%), initially only in proximity of the 2010 index site. Specimens belonging to approximately half of the native amphibian species were infected, including the threatened Pelobates fuscus (spadefoot toad). Clustered massive outbreaks involving dead adult specimens and ranavirus genomic identity indicated that one common midwife toad virus (CMTV)-like ranavirus strain is emerging in provinces in the north of the Netherlands. Modelling based on the spatiotemporal pattern of spread showed a high probability that this emerging virus will continue to be detected at new sites (the discrete reproductive power of this outbreak is 0.35). Phylogenetically distinct CMTV-like ranaviruses were found in the south of the Netherlands more recently. In addition to showing that CMTV-like ranaviruses threaten wild amphibian populations not only in Spain but also in the Netherlands, the current spread and risk of establishment reiterate that understanding the underlying causes of CMTV-like ranavirus emergence requires international attention.

Expanding Distribution of Lethal Amphibian Fungus Batrachochytrium salamandrivorans in Europe.

Emerging fungal diseases can drive amphibian species to local extinction. During 2010-2016, we examined 1,921 urodeles in 3 European countries. Presence of the chytrid fungus Batrachochytrium salamandrivorans at new locations and in urodeles of different species expands the known geographic and host range of the fungus and underpins its imminent threat to biodiversity.

Environmental determinants of recent endemism of Batrachochytrium dendrobatidis infections in amphibian assemblages in the absence of disease outbreaks.

The inconsistent distribution of large-scale infection mediated die-offs and the subsequent population declines of several animal species, urges us to understand how, when, and why species are affected by disease. It is often unclear when or under what conditions a pathogen constitutes a threat to a host. Often, variation of environmental conditions plays a role. Globally Batrachochytrium dendrobatidis (Bd) causes amphibian declines; however, host responses are inconsistent and this fungus appears equally capable of reaching a state of endemism and subsequent co-existence with native amphibian assemblages. We sought to identify environmental and temporal factors that facilitate host-pathogen coexistence in northern Europe. To do this, we used molecular diagnostics to examine archived and wild amphibians for infection and general linear mixed models to explore relationships between environmental variables and prevalence of infection in 5 well-sampled amphibian species. We first detected infection in archived animals collected in 1999, and infection was ubiquitous, but rare, throughout the study period (2008-2010). Prevalence of infection exhibited significant annual fluctuations. Despite extremely rare cases of lethal chytridiomycosis in A. obstetricans, Bd prevalence was uncorrelated with this species' population growth. Our results suggest context dependent and species-specific host susceptibility. Thus, we believe recent endemism of Bd coincides with environmentally driven Bd prevalence fluctuations that preclude the build-up of Bd infection beyond the critical threshold for large-scale mortality and host population crashes.

The absence of zoonotic agents in invasive bullfrogs (Lithobates catesbeianus) in Belgium and The Netherlands.

Exotic invasive bullfrogs (Lithobates catesbeianus) are considered to exert a considerable negative impact on native amphibian communities. This can be due to competition and predation, but they are also a notorious source of the infectious diseases chytridiomycosis and ranavirosis, affecting amphibian populations globally. Little is known regarding their carriage of other microbial agents that might be transferred to humans or other animals. In this study we determined the occurrence of the amphibian pathogens Ranavirus and Batrachochytrium dendrobatidis and of the zoonotic agents Coxiella burnetii, Neospora caninum, Leptospira sp., Toxoplasma gondii, Mycoplasma sp., Campylobacter sp., Salmonella sp. and extended-spectrum beta-lactamase producing Escherichia coli in 164 bullfrogs from three populations in Belgium and The Netherlands. Although B. dendrobatidis was present at a high prevalence of 63%, mean infection loads were low with an average of 10.9 genomic equivalents (SD 35.5), confirming the role of bullfrogs as B. dendrobatidis carriers, but questioning their role as primary reservoirs for B. dendrobatidis transmission to native amphibian communities. All tested samples were negative for the other infectious agents examined. These results suggest a limited role of bullfrogs as carrier of these pathogens.

Batrachochytrium salamandrivorans sp. nov. causes lethal chytridiomycosis in amphibians.

The current biodiversity crisis encompasses a sixth mass extinction event affecting the entire class of amphibians. The infectious disease chytridiomycosis is considered one of the major drivers of global amphibian population decline and extinction and is thought to be caused by a single species of aquatic fungus, Batrachochytrium dendrobatidis. However, several amphibian population declines remain unexplained, among them a steep decrease in fire salamander populations (Salamandra salamandra) that has brought this species to the edge of local extinction. Here we isolated and characterized a unique chytrid fungus, Batrachochytrium salamandrivorans sp. nov., from this salamander population. This chytrid causes erosive skin disease and rapid mortality in experimentally infected fire salamanders and was present in skin lesions of salamanders found dead during the decline event. Together with the closely related B. dendrobatidis, this taxon forms a well-supported chytridiomycete clade, adapted to vertebrate hosts and highly pathogenic to amphibians. However, the lower thermal growth preference of B. salamandrivorans, compared with B. dendrobatidis, and resistance of midwife toads (Alytes obstetricans) to experimental infection with B. salamandrivorans suggest differential niche occupation of the two chytrid fungi.

The novel 'Candidatus Amphibiichlamydia ranarum' is highly prevalent in invasive exotic bullfrogs (Lithobates catesbeianus).

Knowledge concerning microbial infectious diseases in the current amphibian crisis is rudimentary and largely limited to ranavirosis and chytridiomycosis. The family Chlamydiaceae is gaining attention as a common cause of disease in amphibians and may harbour new and emerging amphibian pathogens. We identified a novel species of Chlamydiales (Candidatus Amphibiichlamydia ranarum) with a prevalence of 71% in exotic invasive bullfrog tadpoles (Lithobates catesbeianus) from an introduced population in the Netherlands. The sequence of a 1474 bp 16S rRNA gene fragment showed that the novel taxon forms a well-defined clade with 'Candidatus Amphibiichlamydia salamandrae' within the Chlamydiaceae family. Although none of the tadpoles examined showed signs of clinical disease, urgent evaluation of its pathogenic potential for native amphibian species is required.

Role of sand lizards in the ecology of Lyme and other tick-borne diseases in the Netherlands.

BACKGROUND: Lizards are considered zooprophylactic for almost all Borrelia burgdorferi species, and act as dilution hosts in parts of North America. Whether European lizards significantly reduce the ability of B. burgdorferi to maintain itself in enzootic cycles, and consequently decrease the infection rate of Ixodes ricinus ticks for B. burgdorferi and other tick-borne pathogens in Western Europe is not clear. RESULTS: Ticks were collected from sand lizards, their habitat (heath) and from the adjacent forest. DNA of tick-borne pathogens was detected by PCR followed by reverse line blotting. Tick densities were measured at all four locations by blanket dragging. Nymphs and adult ticks collected from lizards had a significantly lower (1.4%) prevalence of B. burgdorferi sensu lato, compared to questing ticks in heath (24%) or forest (19%). The prevalence of Rickettsia helvetica was significantly higher in ticks from lizards (19%) than those from woodland (10%) whereas neither was significantly different from the prevalence in ticks from heather (15%). The prevalence of Anaplasma and Ehrlichia spp in heather (12%) and forest (14%) were comparable, but significantly lower in ticks from sand lizards (5.4%). The prevalence of Babesia spp in ticks varied between 0 and 5.3%. Tick load of lizards ranged from 1 - 16. Tick densities were ~ 5-fold lower in the heather areas than in woodlands at all four sites. CONCLUSIONS: Despite their apparent low reservoir competence, the presence of sand lizards had insignificant impact on the B. burgdorferi s.l. infection rate of questing ticks. In contrast, sand lizards might act as reservoir hosts for R. helvetica. Remarkably, the public health risk from tick-borne diseases is approximately five times lower in heather than in woodland, due to the low tick densities in heather.