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

Batrachochytrium dendrobatidis strain affects transcriptomic response in liver but not skin in latitudinal populations of the common toad (Bufo bufo)

Batrachochytrium dendrobatidis (Bd) is a fungal pathogen that has decimated amphibian populations worldwide for several decades. We examined the changes in gene expression in response to Bd infection in two populations of the common toad, Bufo bufo, in a laboratory experiment. We collected B. bufo eggs in southern and northern Sweden, and infected the laboratory-raised metamorphs with two strains of the global panzoonotic lineage Bd-GPL. Differential expression analysis showed significant differences between infected and control individuals in both liver and skin. The skin samples showed no discernible differences in gene expression between the two strains used, while liver samples were differentiated by strain, with one of the strains eliciting no immune response from infected toads. Immune system genes were overexpressed in skin samples from surviving infected individuals, while in liver samples the pattern was more diffuse. Splitting samples by population revealed a stronger immune response in northern individuals. Differences in transcriptional regulation between populations are particularly relevant to study in Swedish amphibians, which may have experienced varying exposure to Bd. Earlier exposure to this pathogen and subsequent adaptation or selection pressure may contribute to the survival of some populations over others, while standing genetic diversity in different populations may also affect the infection outcome.

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.

Body size mediates latitudinal population differences in the response to chytrid fungus infection in two amphibians

Factors behind intraspecific variation in sensitivity to pathogens remain poorly understood. We investigated how geographical origin in two North European amphibians affects tolerance to infection by the chytrid fungus Batrachochytrium dendrobatidis (Bd), a generalist pathogen which has caused amphibian population declines worldwide. We exposed newly metamorphosed individuals of moor frog Rana arvalis and common toad Bufo bufo from two latitudinal regions to two different BdGPL strains. We measured survival and growth as infections may cause sub-lethal effects in fitness components even in the absence of mortality. Infection loads were higher in B. bufo than in R. arvalis, and smaller individuals had generally higher infection loads. B. bufo had high mortality in response to Bd infection, whereas there was little mortality in R. arvalis. Bd-mediated mortality was size-dependent and high-latitude individuals were smaller leading to high mortality in the northern B. bufo. Bd exposure led to sub-lethal effects in terms of reduced growth suggesting that individuals surviving the infection may have reduced fitness mediated by smaller body size. In both host species, the Swedish Bd strain caused stronger sublethal effects than the British strain. We suggest that high-latitude populations can be more vulnerable to chytrids than those from lower latitudes and discuss the possible mechanisms how body size and host geographical origin contribute to the present results.

Association between the skin microbiome and MHC class II diversity in an amphibian

Microbiomes play an important role in determining the ecology and behaviour of their hosts. However, questions remain pertaining to how host genetics shape microbiomes, and how microbiome composition influences host fitness. We explored the effects of geography, evolutionary history and host genetics on the skin microbiome diversity and structure in a widespread amphibian. More specifically, we examined the association between bacterial diversity and composition and the major histocompatibility complex class II exon 2 diversity in 12 moor frog (Rana arvalis) populations belonging to two geographical clusters that show signatures of past and ongoing differential selection. We found that while bacterial alpha diversity did not differ between the two clusters, MHC alleles/supertypes and genetic diversity varied considerably depending on geography and evolutionary history. Bacterial alpha diversity was positively correlated with expected MHC heterozygosity and negatively with MHC nucleotide diversity. Furthermore, bacterial community composition showed significant variation between the two geographical clusters and between specific MHC alleles/supertypes. Our findings emphasize the importance of historical demographic events on hologenomic variation and provide new insights into how immunogenetic host variability and microbial diversity may jointly influence host fitness with consequences for disease susceptibility and population persistence.

Amphibians and reptiles in North Sweden: distribution, habitat affinities, and abundance (Classes: Amphibia and Reptilia)

North Sweden ('Norrland' in Swedish) covers 243 000 km2 and lies mainly in the boreal biome. The herpetofauna comprises five amphibian and four reptile species: Lissotriton vulgaris, Triturus cristatus, Bufo bufo, Rana temporaria, Rana arvalis, Zootoca vivipara, Anguis fragilis, Natrix natrix, and Vipera berus. Successful conservation and management of amphibians and reptiles depend on accurate information about distribution, habitat affinities, and abundance. Such knowledge is also essential as a benchmark to assess changes in distribution and abundance that may come about as a result of climate change and human habitat alteration. This paper aims to present accurate distribution maps, describe habitat affinities, and provide abundance estimates for the herpetofauna of North Sweden for the period 1970-2022. Distribution data are presented by traditional faunistic provinces, as well as by biotic regions and alpine life zones. Separate sections address post-glacial colonization and a herpetological perspective on anthropogenic changes in relation to species´ present status. Bufo bufo, Rana temporaria, Rana arvalis, Zootoca vivipara, and Vipera berus are widely distributed throughout the boreal sub-regions. Rana temporaria, Zootoca vivipara, and Vipera berus also extend into the alpine region. Triturus cristatus, Anguis fragilis, and Natrix natrix occur mainly in the coastal parts of the Southern Boreal region. There are no signs of recent changes in distribution range, but Lissotriton vulgaris, Triturus cristatus, and Rana arvalis have been largely overlooked in the past and have a much wider occurrence than previously recognized. Most species are found in habitats usually not described in all-European field guides. Nearly all anurans hibernate in water. Abundance estimates suggest that some species are more common in the boreal than thought, supporting the notion that a large share of their total European population occurs there. Although local extinctions and declines are known, there are no signs of widespread population decline for any species during the study period.

Small-scale population divergence is driven by local larval environment in a temperate amphibian

Genomic variation within and among populations is shaped by the interplay between natural selection and the effects of genetic drift and gene flow. Adaptive divergence can be found in small-scale natural systems even when population sizes are small, and the potential for gene flow is high, suggesting that local environments exert selection pressures strong enough to counteract the opposing effects of drift and gene flow. Here, we investigated genomic differentiation in nine moor frog (Rana arvalis) populations in a small-scale network of local wetlands using 16,707 ddRAD-seq SNPs, relating levels of differentiation with local environments, as well as with properties of the surrounding landscape. We characterized population structure and differentiation, and partitioned the effects of geographic distance, local larval environment, and landscape features on total genomic variation. We also conducted gene-environment association studies using univariate and multivariate approaches. We found small-scale population structure corresponding to 6-8 clusters. Local larval environment was the most influential component explaining 2.3% of the total genetic variation followed by landscape features (1.8%) and geographic distance (0.8%), indicative of isolation-by-environment, -by-landscape, and -by-distance, respectively. We identified 1000 potential candidate SNPs putatively under divergent selection mediated by the local larval environment. The candidate SNPs were involved in, among other biological functions, immune system function and development. Our results suggest that small-scale environmental differences can exert selection pressures strong enough to counteract homogenizing effects of gene flow and drift in this small-scale system, leading to observable population differentiation.