An emerging viral pathogen truncates population age structure in a European amphibian and may reduce population viability

Changes in capture availability due to infection can lead to detectable biases in population-level infectious disease parameters

Correctly identifying the strength of selection that parasites impose on hosts is key to predicting epidemiological and evolutionary outcomes of host-parasite interactions. However, behavioral changes due to infection can alter the capture probability of infected hosts and thereby make selection difficult to estimate by standard sampling techniques. Mark-recapture approaches, which allow researchers to determine if some groups in a population are less likely to be captured than others, can be used to identify infection-driven capture biases. If a metric of interest directly compares infected and uninfected populations, calculated detection probabilities for both groups may be useful in identifying bias. Here, we use an individual-based simulation to test whether changes in capture rate due to infection can alter estimates of three key metrics: 1) reduction in the reproductive success of infected parents relative to uninfected parents, 2) the relative risk of infection for susceptible genotypes compared to resistant genotypes, and 3) changes in allele frequencies between generations. We explore the direction and underlying causes of the biases that emerge from these simulations. Finally, we argue that short series of mark-recapture sampling bouts, potentially implemented in under a week, can yield key data on detection bias due to infection while not adding a significantly higher burden to disease ecology studies.

A novel tamanavirus (Flaviviridae) of the European common frog (Rana temporaria) from the UK

Flavivirids are small, enveloped, positive-sense RNA viruses from the family Flaviviridae with genomes of ~9-13 kb. Metatranscriptomic analyses of metazoan organisms have revealed a diversity of flavivirus-like or flavivirid viral sequences in fish and marine invertebrate groups. However, no flavivirus-like virus has been identified in amphibians. To remedy this, we investigated the virome of the European common frog (Rana temporaria) in the UK, utilizing high-throughput sequencing at six catch locations. De novo assembly revealed a coding-complete virus contig of a novel flavivirid ~11.2 kb in length. The virus encodes a single ORF of 3456 aa and 5' and 3' untranslated regions (UTRs) of 227 and 666 nt, respectively. We named this virus Rana tamanavirus (RaTV), as BLASTp analysis of the polyprotein showed the closest relationships to Tamana bat virus (TABV) and Cyclopterus lumpus virus from Pteronotus parnellii and Cyclopterus lumpus, respectively. Phylogenetic analysis of the RaTV polyprotein compared to Flavivirus and Flavivirus-like members indicated that RaTV was sufficiently divergent and basal to the vertebrate Tamanavirus clade. In addition to the Mitcham strain, partial but divergent RaTV, sharing 95.64-97.39 % pairwise nucleotide identity, were also obtained from the Poole and Deal samples, indicating that RaTV is widespread in UK frog samples. Bioinformatic analyses of predicted secondary structures in the 3'UTR of RaTV showed the presence of an exoribonuclease-resistant RNA (xrRNA) structure standard in flaviviruses and TABV. To examine this biochemically, we conducted an in vitro Xrn1 digestion assay showing that RaTV probably forms a functional Xrn1-resistant xrRNA.

Widespread amphibian Perkinsea infections associated with Ranidae hosts, cooler months and Ranavirus co-infection

Amphibians suffer from large-scale population declines globally, and emerging infectious diseases contribute heavily to these declines. Amphibian Perkinsea (Pr) is a worldwide anuran pathogen associated with mass mortality events, yet little is known about its epidemiological patterns, especially in comparison to the body of literature on amphibian chytridiomycosis and ranavirosis. Here, we establish Pr infection patterns in natural anuran populations and identify important covariates including climate, host attributes and co-infection with Ranavirus (Rv). We used quantitative (q)PCR to determine the presence and intensity of Pr and Rv across 1234 individuals sampled throughout central Florida in 2017-2019. We then implemented random forest ensemble learning models to predict infection with both pathogens based on physiological and environmental characteristics. Perkinsea infected 32% of all sampled anurans, and Pr prevalence was significantly elevated in Ranidae frogs, cooler months, metamorphosed individuals and frogs co-infected with Rv, while Pr intensity was significantly higher in ranid frogs and individuals collected dead. Ranavirus prevalence was 17% overall and was significantly higher in Ranidae frogs, metamorphosed individuals, locations with higher average temperatures, and individuals co-infected with Pr. Perkinsea prevalence was significantly higher than Rv prevalence across months, regions, life stages and species. Among locations, Pr prevalence was negatively associated with crayfish prevalence and positively associated with relative abundance of microhylids, but Rv prevalence did not associate with any tested co-variates. Co-infections were significantly more common than single infections for both pathogens, and we propose that Pr infections may propel Rv infections because seasonal Rv infection peaks followed Pr infection peaks and random forest models found Pr intensity was a leading factor explaining Rv infections. Our study elucidates epidemiological patterns of Pr in Florida and suggests that Pr may be under-recognized as a cause of anuran declines, especially in the context of pathogen co-infection.

Age structure of amphibian populations with endemic chytridiomycosis, across climatic regions with markedly different infection risk

Threatening processes, such as disease, can drive major changes in population demographics of the host. Chytridiomycosis, caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), has led to the decline of at least 500 amphibian species across the globe and has been shown to truncate host age structure by lowering adult survival rates. This results in heavy reliance on annual recruitment and the inability to recover in the event of periodic recruitment failure. We used skeletochronology to determine the age structure, growth, and survival rates of populations of an endangered amphibian, Litoria raniformis, with endemic chytridiomycosis, across two climatically disparate regions in south‐eastern Australia: semi‐arid and temperate. Contrary to predictions, populations in the semi‐arid region (in which chytrid prevalence is substantially lower due to high temperatures) displayed a more truncated age structure than populations in the temperate study regions. Maximum recorded age was only two years in the semi‐arid region compared with up to four years in the temperate region. Wetland hydroperiod and average seasonal air temperature were correlated with age, and males had a slightly higher survival rate than females (0.31 for males and 0.27 for females). Despite the previously documented differences in chytrid prevalence between the two climatic regions, water availability and wetland hydroperiods appear the over‐riding determinants of the age structure and survival rates of L. raniformis. Targeted management which ensures water availability and improves survival of 1‐year‐old frogs into their second and third breeding season would reduce the impact of stochastic events on L. raniformis, and this may be true for numerous frog species susceptible to chytridiomycosis.

A New Disease Caused by an Unidentified Etiological Agent Affects European Salamanders

New pathologies are causing dramatic declines and extinctions of multiple amphibian species. In 2013, in one fire salamander population of Northern Italy, we found individuals with undescribed cysts at the throat level, a malady whose existence has not previously been reported in amphibians. With the aim of describing this novel disease, we performed repeated field surveys to assess the frequency of affected salamanders from 2014 to 2020, and integrated morphological, histological, and molecular analyses to identify the pathogen. The novel disease affected up to 22% of salamanders of the study population and started spreading to nearby populations. Cysts are formed by mucus surrounding protist-like cells about 30 µm long, characterized by numerous cilia/undulipodia. Morphological and genetic analyses did not yield a clear match with described organisms. The existence of this pathogen calls for the implementation of biosecurity protocols and more studies on the dynamics of transmission and the impact on wild populations.

Age, sex, and inflammatory markers predict chronic conditions, cardiac disease, and mortality among captive western lowland gorillas (Gorilla gorilla gorilla)

In humans, inflammatory markers predict health risks. As great apes experience many similar conditions, measuring inflammation may provide valuable health information. We examined four serum inflammatory markers in zoo-housed gorillas (n = 48): albumin, CRP, IL-6, and TNF-α. We first analyzed age- and sex-associated patterns, then used multimodel inference to evaluate models with age, sex, and inflammatory markers as predictors of all-cause morbidity, cardiac disease, and mortality. Older gorillas had lower albumin and higher IL-6, and males had higher albumin, lower CRP, and lower TNF-α. All-cause morbidity was best predicted by age, sex, and TNF-α, but the second model containing only age and sex was equivalent. Cardiac disease was best predicted by TNF-α alongside age and sex, with lower levels associated with increased risk. When outliers were removed, the model with TNF-α was second to the model containing only age and sex. Finally, mortality risk was best predicted by the model with only age and sex. Other models containing individual inflammatory markers were within top model sets for each health outcome. Our results indicate that age and sex are robust for predicting all-cause morbidity and mortality risk in gorillas; while models which include individual inflammatory markers also predict risk, they may not improve predictions over age and sex alone. However, given the prevalence of cardiac disease in great apes, these results suggest that TNF-α warrants further investigation. With their potential to provide valuable health information, data on inflammatory markers may contribute to the care and management of gorillas in human care.

Age, size and body condition do not equally reflect population response to habitat change in the common spadefoot toad Pelobates fuscus

Urbanization impacts biodiversity both directly through physical expansion over land, and indirectly due to land use conversion and human behaviors associated with urban areas. We assessed the response of a common spadefoot toad population (Pelobates fuscus) to habitat loss and fragmentation resulting from urban development by studying changes in size, body condition and age parameters. We compared samples collected in the early 2000s (sample A) and later on during 2012-2014 (sample B). The terrestrial habitats in the study area were severely reduced and fragmented due to the expansion of the human settlement. We found no significant differences in the age parameters between the two sampling periods; the median lifespan shortened from 3.5 (sample A) to 3.0 years (sample B), while the other age parameters were similar in both samples. In contrast, snout-vent length, body mass and body condition experienced a significant decrease over time. Our results suggest that changes in body size and body condition, rather than age parameters, better reflect the response of the common spadefoot toad population to declining habitat quality. Therefore, body measurements can provide reliable estimates of the impact of habitat degradation in amphibian populations.

Amphibian ranaviruses in Europe: important directions for future research

Ranaviruses are an emerging group of pathogens capable of infecting all cold-blooded vertebrates. In Europe, ranaviruses pose a particularly potent threat to wild amphibian populations. Since the 1980s research on amphibian-infecting ranaviruses in Europe has been growing. The wide distribution of amphibian populations in Europe, the ease with which many are monitored, and the tractable nature of counterpart ex situ experimental systems have provided researchers with a unique opportunity to study many aspects of host–ranavirus interactions in the wild. These characteristics of European amphibian populations will also enable researchers to lead the way as the field of host–ranavirus interactions progresses. In this review, we provide a summary of the current key knowledge regarding amphibian infecting ranaviruses throughout Europe. We then outline important areas of further research and suggest practical ways each could be pursued. We address the study of potential interactions between the amphibian microbiome and ranaviruses, how pollution may exacerbate ranaviral disease either as direct stressors of amphibians or indirect modification of the amphibian microbiome. Finally, we discuss the need for continued surveillance of ranaviral emergence in the face of climate change.

Diversity-Stability Dynamics of the Amphibian Skin Microbiome and Susceptibility to a Lethal Viral Pathogen

Variation among animals in their host-associated microbial communities is increasingly recognized as a key determinant of important life history traits including growth, metabolism, and resistance to disease. Quantitative estimates of the factors shaping the stability of host microbiomes over time at the individual level in non-model organisms are scarce. Addressing this gap in our knowledge is important, as variation among individuals in microbiome stability may represent temporal gain or loss of key microbial species and functions linked to host health and/or fitness. Here we use controlled experiments to investigate how both heterogeneity in microbial species richness of the environment and exposure to the emerging pathogen Ranavirus influence the structure and temporal dynamics of the skin microbiome in a vertebrate host, the European common frog (Rana temporaria). Our evidence suggests that altering the bacterial species richness of the environment drives divergent temporal microbiome dynamics of the amphibian skin. Exposure to ranavirus effects changes in skin microbiome structure irrespective of total microbial diversity, but individuals with higher pre-exposure skin microbiome diversity appeared to exhibit higher survival. Higher diversity skin microbiomes also appear less stable over time compared to lower diversity microbiomes, but stability of the 100 most abundant ("core") community members was similar irrespective of microbiome richness. Our study highlights the importance of extrinsic factors in determining the stability of host microbiomes over time, which may in turn have important consequences for the stability of host-microbe interactions and microbiome-fitness correlations.

The fast-slow continuum of longevity among yellow-bellied toad populations (Bombina variegata): intrinsic and extrinsic drivers of variation

Yellow-bellied toad populations (Bombina variegata) show a wide fast-slow continuum of the life-history trait longevity ranging from 5 to 23 years. We investigated populations in Germany (n = 8) and Austria (n = 1) to determine their position within the continuum of longevity and the potential drivers of adult survival at the local and the continental scale. Intrinsic and extrinsic factors considered were local weather, nutritional state, allocation of ingested energy to somatic growth, pathogen prevalence, and geographical clines (latitude, altitude, and longitude). Capture-mark-recapture (CMR) monitoring and direct age assessment by skeletochronology allowed for reliable estimates of longevity and adult survival. Raw and corrected recapture rates as well as a probabilistic estimate of the lifespan of the eldest 1% adults of a cohort (CMR data) were used as surrogates for adult survival and thus longevity in a population. Additionally, survival rates were calculated from static life tables based on the age structure (skeletochronological data) of eight populations. Populations in Germany were short-lived with a maximum lifespan of annual cohorts varying from 5 to 8 years, whereas the population in Austria was long-lived with a cohort longevity of 13 to 23 years. We provide evidence that annual survival rates and longevity differ among years and between short- and long-lived populations, but there was no decrease of survival in older toads (i.e. absence of senescence). Variation of weather among years accounted for 90.7% of variance in annual survival rates of short-lived populations, whereas the sources of variation in the long-lived population remained unidentified. At the continental scale, longevity variation among B. variegata populations studied so far did not correspond to geographical clines or climate variation. Therefore, we propose that a population's position within the fast-slow continuum integrates the response to local environmental stochasticity (extrinsic source of variation) and the efficiency of chemical antipredator protection determining the magnitude of longevity (intrinsic source of variation).

LETHAL AND SUBLETHAL AMPHIBIAN HOST RESPONSES TO BATRACHOCHYTRIUM DENDROBATIDIS EXPOSURE ARE DETERMINED BY THE ADDITIVE INFLUENCE OF HOST RESOURCE AVAILABILITY

Host species may differ in their responses to pathogen exposures based on host energy reserves, which could be important for long-term trends in host population growth. Batrachochytrium dendrobatidis (BD) is a pathogen associated with amphibian population declines but also occurs without causing mass mortalities. The impact of BD in populations without associated declines is not well understood, and food abundance could play a role in determining the magnitude of its effects. We exposed American toad (Anaxyrus americanus), northern leopard frog (Lithobates pipiens), and cricket frog (Acris blanchardi) metamorphs to BD under low or high food treatments. Overall, anuran species responded differently to BD exposure and the combined effect of BD exposure and food abundance was additive. American toad survival was lowered by BD exposure and low food availability. Based on these results, we developed a population model for American toads to estimate how reductions in survival could influence population growth. We found that BD could reduce population growth by 14% with high food availability and 21% with low food availability. In contrast, survival of northern leopard frogs was high across all treatments, but their growth was negatively impacted by the additive effects of BD exposure and low food availability. Cricket frog growth and survival were unaffected by BD exposure, suggesting that this species is not sensitive to the effects of this pathogen in terms of growth and survival across environments of different quality in the time period examined. Our results showed that low food availability additively increased the species-specific lethal and sublethal impacts of BD on hosts, which could have implications for long-term host population dynamics.

Outbreaks of an Emerging Viral Disease Covary With Differences in the Composition of the Skin Microbiome of a Wild United Kingdom Amphibian

There is growing appreciation of the important role of commensal microbes in ensuring the normal function and health of their hosts, including determining how hosts respond to pathogens. A range of infectious diseases are threatening amphibians worldwide, and evidence is accumulating that the host-associated bacteria that comprise the microbiome may be key in mediating interactions between amphibian hosts and infectious pathogens. We used 16S rRNA amplicon sequencing to quantify the skin microbial community structure of over 200 individual wild adult European common frogs (Rana temporaria), from ten populations with contrasting history of the lethal disease ranavirosis, caused by emerging viral pathogens belonging to the genus Ranavirus. All populations had similar species richness irrespective of disease history, but populations that have experienced historical outbreaks of ranavirosis have a distinct skin microbiome structure (beta diversity) when compared to sites where no outbreaks of the disease have occurred. At the individual level, neither age, body length, nor sex of the frog could predict the structure of the skin microbiota. Our data potentially support the hypothesis that variation among individuals in skin microbiome structure drive differences in susceptibility to infection and lethal outbreaks of disease. More generally, our results suggest that population-level processes are more important for driving differences in microbiome structure than variation among individuals within populations in key life history traits such as age and body size.

Modelling Ranavirus Transmission in Populations of Common Frogs (Rana temporaria) in the United Kingdom

Ranaviruses began emerging in common frogs (Rana temporaria) in the United Kingdom in the late 1980s and early 1990s, causing severe disease and declines in the populations of these animals. Herein, we explored the transmission dynamics of the ranavirus(es) present in common frog populations, in the context of a simple susceptible-infected (SI) model, using parameters derived from the literature. We explored the effects of disease-induced population decline on the dynamics of the ranavirus. We then extended the model to consider the infection dynamics in populations exposed to both ulcerative and hemorrhagic forms of the ranaviral disease. The preliminary investigation indicated the important interactions between the forms. When the ulcerative form was present in a population and the hemorrhagic form was later introduced, the hemorrhagic form of the disease needed to be highly contagious, to persist. We highlighted the areas where further research and experimental evidence is needed and hope that these models would act as a guide for further research into the amphibian disease dynamics.