Nor climate nor human impact factors: Chytrid infection shapes the skin bacterial communities of an endemic amphibian in a biodiversity hotspot.

The bacterial communities of the amphibian skin (i.e., the bacteriome) are critical to the host's innate immune system. However, it is unclear how different drivers can alter this function by modulating the bacteriome's structure. Our aim was to assess the extent to which different host attributes and extrinsic factors influence the structure of the bacterial communities of the skin. Skin bacterial diversity was examined in 148 individuals of the four-eyed frog (Pleurodema thaul) from 16 localities spanning almost 1800 km in latitude. The richness and beta diversity of bacterial families and the richness and abundance of Bd-inhibitory bacterial genera were used to describe their structure. Predictors associated with the host (developmental stage, genetic lineage, individual Batrachochytrium dendrobatidis [Bd] infection status) and the landscape (current climate, degree of anthropogenic disturbance) were used in the statistical modeling in an information theoretical approach. Bd infection and host developmental stage were the only predictors affecting bacteriome richness, with Bd+ individuals and postmetamorphic stages (adults and juveniles) having higher richness than Bd- ones and tadpoles. High diversity in Bd+ individuals is not driven by bacterial genera with known anti-Bd properties. Beta diversity was not affected by Bd infection and was mainly a consequence of bacterial family turnover rather than nestedness. Finally, for those bacterial genera known to have inhibitory effects on chytrid, Bd+ individuals had a slightly higher diversity than Bd- ones. Our study confirms an association between Bd infection and the host developmental stage with the bacterial communities of the skin of P. thaul. Unexpectedly, macroclimate and human impact factors do not seem to play a role in shaping the amphibian skin microbiome. Our study exemplifies that focusing on a single host-parasite system over a large geographic scale can provide essential insights into the factors driving host-parasite-bacteriome interactions.

Interpopulation differences in male reproductive effort drive the population dynamics of a host exposed to an emerging fungal pathogen.

Compensatory recruitment is a key demographic mechanism that has allowed the coexistence of populations of susceptible amphibians with Batrachochytrium dendrobatidis (Bd), a fungus causing one of the most devastating emerging infectious disease ever recorded among vertebrates. However, the underlying processes (e.g. density-dependent increase in survival at early life stages, change in reproductive traits) as well as the level of interpopulation variation in this response are poorly known. We explore potential mechanisms of compensatory recruitment in response to Bd infection by taking advantage of an amphibian system where male reproductive traits are easy to quantify in free-living populations. The Southern Darwin's frog Rhinoderma darwinii is a vocal sac-brooding species that exhibits a high susceptibility to lethal Bd infection. Using a 7-year capture-recapture study at four populations with contrasting Bd infection status (one high prevalence, one low prevalence and two Bd-free populations), we evaluated whether Bd-positive populations exhibited a higher adult recruitment and a higher male reproductive effort than Bd-negative populations. We also estimated population growth rates to explore whether recruitment compensated for the negative impacts of Bd on the survival of adults. In addition, we evaluated a potential demographic signal of compensatory recruitment (i.e. positive relationship between the proportion of juveniles and Bd prevalence) in response to Bd infection using raw count data from 13 R. darwinii populations. The high Bd prevalence population exhibited the highest male reproductive effort and the highest recruitment among the four monitored populations. This led to a growing population during the study period despite high mortality of adult hosts. In contrast, males from the population with low Bd prevalence had a low reproductive effort and this population, which had the lowest adult recruitment, was declining during the study period despite adults having a higher survival in comparison to the high Bd prevalence population. We also found a demographic signal of compensatory recruitment in response to Bd infection in our broader analysis of 13 R. darwinii populations. Our study underlines the importance of interpopulation variation in life-history strategies on the fate of host populations after infectious disease emergence. Our results also suggest that an increase in reproductive effort can be one of the processes underlying compensatory recruitment in populations of Bd-susceptible amphibians.

Sex-related differences in aging rate are associated with sex chromosome system in amphibians.

Sex-related differences in mortality are widespread in the animal kingdom. Although studies have shown that sex determination systems might drive lifespan evolution, sex chromosome influence on aging rates have not been investigated so far, likely due to an apparent lack of demographic data from clades including both XY (with heterogametic males) and ZW (heterogametic females) systems. Taking advantage of a unique collection of capture-recapture datasets in amphibians, a vertebrate group where XY and ZW systems have repeatedly evolved over the past 200 million years, we examined whether sex heterogamy can predict sex differences in aging rates and lifespans. We showed that the strength and direction of sex differences in aging rates (and not lifespan) differ between XY and ZW systems. Sex-specific variation in aging rates was moderate within each system, but aging rates tended to be consistently higher in the heterogametic sex. This led to small but detectable effects of sex chromosome system on sex differences in aging rates in our models. Although preliminary, our results suggest that exposed recessive deleterious mutations on the X/Z chromosome (the "unguarded X/Z effect") or repeat-rich Y/W chromosome (the "toxic Y/W effect") could accelerate aging in the heterogametic sex in some vertebrate clades.

Chytridiomycosis Outbreak in a Chilean Giant Frog (Calyptocephalella gayi) Captive Breeding Program: Genomic Characterization and Pathological Findings.

Emerging infectious diseases in wildlife are increasingly associated with animal mortality and species declines, but their source and genetic characterization often remains elusive. Amphibian chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), has been associated with catastrophic and well-documented amphibian population declines and extinctions at the global scale. We used histology and whole-genome sequencing to describe the lesions caused by, and the genetic variability of, two Bd isolates obtained from a mass mortality event in a captive population of the threatened Chilean giant frog (Calyptocephalella gayi). This was the first time an association between Bd and high mortality had been detected in this charismatic and declining frog species. Pathological examinations revealed that 30 dead metamorphosed frogs presented agnathia or brachygnathia, a condition that is reported for the first time in association with chytridiomycosis. Phylogenomic analyses revealed that Bd isolates (PA1 and PA2) from captive C. gayi group with other Bd isolates (AVS2, AVS4, and AVS7) forming a single highly supported Chilean Bd clade within the global panzootic lineage of Bd (BdGPL). These findings are important to inform the strengthening of biosecurity measures to prevent the impacts of chytridiomycosis in captive breeding programs elsewhere.

Bioclimatic and anthropogenic variables shape the occurrence of Batrachochytrium dendrobatidis over a large latitudinal gradient.

Amphibian chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), has caused the greatest known loss of biodiversity due to an infectious disease. We used Bd infection data from quantitative real-time PCR (qPCR) assays of amphibian skin swabs collected across Chile during 2008-2018 to model Bd occurrence with the aim to determine bioclimatic and anthropogenic variables associated with Bd infection. Also, we used Bd presence/absence records to identify geographical Bd high-risk areas and compare Bd prevalence and infection loads between amphibian families, ecoregions, and host ecology. Data comprised 4155 Bd-specific qPCR assays from 162 locations across a latitudinal gradient of 3700 km (18º to 51ºS). Results showed a significant clustering of Bd associated with urban centres and anthropogenically highly disturbed ecosystems in central-south Chile. Both Bd prevalence and Bd infection loads were higher in aquatic than terrestrial amphibian species. Our model indicated positive associations of Bd prevalence with altitude, temperature, precipitation and human-modified landscapes. Also, we found that macroscale drivers, such as land use change and climate, shape the occurrence of Bd at the landscape level. Our study provides with new evidence that can improve the effectiveness of strategies to mitigate biodiversity loss due to amphibian chytridiomycosis.

Towards understanding the effects of oceanic plastic pollution on population growth for a South American fur seal (Arctocephalus australis australis) colony in Chile.

Entanglement of pinnipeds with plastic debris is an emerging conservation and animal welfare issue worldwide. However, the origins and long-term population level consequences of these entanglements are usually unknown. Plastic entanglement could produce a combination of wounds, asphyxiation, or inability to feed that results in the death of a certain percentage of individuals from the total population. In this research, we report on the consequent effect of plastic entanglement on population growth demographics in a South American fur seal (Arctocephalus australis australis) colony on Guafo Island, southern Chile. Using a stochastic matrix population model structured according to age and sex, and assuming an otherwise stable population, we explored population growth rates under five scenarios with differing rates of entanglement: A) a zero rate of plastic entanglement, B) entanglement rates (number of entangled individuals as a proportion of the total number of individuals) as observed in our study population (overall entanglement ratio of 1.2 × 10-3); and for the other scenarios, entanglement ratios as reported in the literature for other pinniped colonies around the world: C) 3.04 × 10-3, D) 4.42 × 10-2, and E) 8.39 × 10-2. Over the 30 years forecasting period and starting with a population size of ∼2950 individuals, the population growth rate was lower under all scenarios with rates of entanglement greater than zero (scenarios B-E). In comparison with scenario A, at the end of the 30-year period forecasted, we calculated a projected decrease in population size of between 20.34% (scenario B) and 91.38% (scenario E). These results suggest that even the lowest levels of entanglement in pinnipeds as reported in the literature might have significant effects over time on population-level dynamics. Our research offers potential insight when devising policy for the management and limitation of plastic pollution in the oceans, and indeed for the conservation and management policy of affected marine species. Furthermore, whilst there are some limitations to our methodology, it offers a straightforward and potentially useful approach for the standardized prediction of impacts at a population level of different rates of plastic pollution and entanglement and could be applied in distinct populations of the same species around the world.

Why disease ecology needs life-history theory: a host perspective.

When facing an emerging infectious disease of conservation concern, we often have little information on the nature of the host-parasite interaction to inform management decisions. However, it is becoming increasingly clear that the life-history strategies of host species can be predictive of individual- and population-level responses to infectious disease, even without detailed knowledge on the specifics of the host-parasite interaction. Here, we argue that a deeper integration of life-history theory into disease ecology is timely and necessary to improve our capacity to understand, predict and mitigate the impact of endemic and emerging infectious diseases in wild populations. Using wild vertebrates as an example, we show that host life-history characteristics influence host responses to parasitism at different levels of organisation, from individuals to communities. We also highlight knowledge gaps and future directions for the study of life-history and host responses to parasitism. We conclude by illustrating how this theoretical insight can inform the monitoring and control of infectious diseases in wildlife.

Recent Asian origin of chytrid fungi causing global amphibian declines.

Globalized infectious diseases are causing species declines worldwide, but their source often remains elusive. We used whole-genome sequencing to solve the spatiotemporal origins of the most devastating panzootic to date, caused by the fungus Batrachochytrium dendrobatidis, a proximate driver of global amphibian declines. We traced the source of B. dendrobatidis to the Korean peninsula, where one lineage, BdASIA-1, exhibits the genetic hallmarks of an ancestral population that seeded the panzootic. We date the emergence of this pathogen to the early 20th century, coinciding with the global expansion of commercial trade in amphibians, and we show that intercontinental transmission is ongoing. Our findings point to East Asia as a geographic hotspot for B. dendrobatidis biodiversity and the original source of these lineages that now parasitize amphibians worldwide.

Development and worldwide use of non-lethal, and minimal population-level impact, protocols for the isolation of amphibian chytrid fungi.

Parasitic chytrid fungi have emerged as a significant threat to amphibian species worldwide, necessitating the development of techniques to isolate these pathogens into culture for research purposes. However, early methods of isolating chytrids from their hosts relied on killing amphibians. We modified a pre-existing protocol for isolating chytrids from infected animals to use toe clips and biopsies from toe webbing rather than euthanizing hosts, and distributed the protocol to researchers as part of the BiodivERsA project RACE; here called the RML protocol. In tandem, we developed a lethal procedure for isolating chytrids from tadpole mouthparts. Reviewing a database of use a decade after their inception, we find that these methods have been applied across 5 continents, 23 countries and in 62 amphibian species. Isolation of chytrids by the non-lethal RML protocol occured in 18% of attempts with 207 fungal isolates and three species of chytrid being recovered. Isolation of chytrids from tadpoles occured in 43% of attempts with 334 fungal isolates of one species (Batrachochytrium dendrobatidis) being recovered. Together, these methods have resulted in a significant reduction and refinement of our use of threatened amphibian species and have improved our ability to work with this group of emerging pathogens.

High prevalence of chigger mite infection in a forest-specialist frog with evidence of parasite-related granulomatous myositis.

Amphibians are hosts for a wide variety of micro- and macro-parasites. Chigger mites from the Hannemania genus are known to infect a wide variety of amphibian species across the Americas. In Chile, three species (H. pattoni, H. gonzaleacunae and H. ortizi) have been described infecting native anurans; however, neither impacts nor the microscopic lesions associated with these parasites have been described. Here, we document 70% prevalence of chigger mite infection in Eupsophus roseus and absence of infection in Rhinoderma darwinii in the Nahuelbuta Range, Chile. Additionally, we describe the macroscopic and microscopic lesions produced by H. ortizi in one of these species, documenting previously undescribed lesions (granulomatous myositis) within the host's musculature. These findings highlight that further research to better understand the impacts of chigger mite infection on amphibians is urgently required in Chile and elsewhere.

Cryptic disease-induced mortality may cause host extinction in an apparently stable host-parasite system.

The decline of wildlife populations due to emerging infectious disease often shows a common pattern: the parasite invades a naive host population, producing epidemic disease and a population decline, sometimes with extirpation. Some susceptible host populations can survive the epidemic phase and persist with endemic parasitic infection. Understanding host-parasite dynamics leading to persistence of the system is imperative to adequately inform conservation practice. Here we combine field data, statistical and mathematical modelling to explore the dynamics of the apparently stable Rhinoderma darwinii-Batrachochytrium dendrobatidis (Bd) system. Our results indicate that Bd-induced population extirpation may occur even in the absence of epidemics and where parasite prevalence is relatively low. These empirical findings are consistent with previous theoretical predictions showing that highly pathogenic parasites are able to regulate host populations even at extremely low prevalence, highlighting that disease threats should be investigated as a cause of population declines even in the absence of an overt increase in mortality.

Dispersal and extrapolation on the accuracy of temporal predictions from distribution models for the Darwin's frog.

Climate change is a major threat to biodiversity; the development of models that reliably predict its effects on species distributions is a priority for conservation biogeography. Two of the main issues for accurate temporal predictions from Species Distribution Models (SDM) are model extrapolation and unrealistic dispersal scenarios. We assessed the consequences of these issues on the accuracy of climate-driven SDM predictions for the dispersal-limited Darwin's frog Rhinoderma darwinii in South America. We calibrated models using historical data (1950-1975) and projected them across 40 yr to predict distribution under current climatic conditions, assessing predictive accuracy through the area under the ROC curve (AUC) and True Skill Statistics (TSS), contrasting binary model predictions against temporal-independent validation data set (i.e., current presences/absences). To assess the effects of incorporating dispersal processes we compared the predictive accuracy of dispersal constrained models with no dispersal limited SDMs; and to assess the effects of model extrapolation on the predictive accuracy of SDMs, we compared this between extrapolated and no extrapolated areas. The incorporation of dispersal processes enhanced predictive accuracy, mainly due to a decrease in the false presence rate of model predictions, which is consistent with discrimination of suitable but inaccessible habitat. This also had consequences on range size changes over time, which is the most used proxy for extinction risk from climate change. The area of current climatic conditions that was absent in the baseline conditions (i.e., extrapolated areas) represents 39% of the study area, leading to a significant decrease in predictive accuracy of model predictions for those areas. Our results highlight (1) incorporating dispersal processes can improve predictive accuracy of temporal transference of SDMs and reduce uncertainties of extinction risk assessments from global change; (2) as geographical areas subjected to novel climates are expected to arise, they must be reported as they show less accurate predictions under future climate scenarios. Consequently, environmental extrapolation and dispersal processes should be explicitly incorporated to report and reduce uncertainties in temporal predictions of SDMs, respectively. Doing so, we expect to improve the reliability of the information we provide for conservation decision makers under future climate change scenarios.

Geographic body size variation in ectotherms: effects of seasonality on an anuran from the southern temperate forest.

BACKGROUND: Body size variation has played a central role in biogeographical research, however, most studies have aimed to describe trends rather than search for underlying mechanisms. In order to provide a more comprehensive understanding of the causes of intra-specific body size variation in ectotherms, we evaluated eight hypotheses proposed in the literature to account for geographical body size variation using the Darwin's frog (Rhinoderma darwinii), an anuran species widely distributed in the temperate forests of South America. Each of the evaluated hypotheses predicted a specific relationship between body size and environmental variables. The level of support for each of these hypotheses was assessed using an information-theoretic approach and based on data from 1015 adult frogs obtained from 14 sites across the entire distributional range of the species. RESULTS: There was strong evidence favouring a single model comprising temperature seasonality as the predictor variable. Larger body sizes were found in areas of greater seasonality, giving support to the "starvation resistance" hypothesis. Considering the known role of temperature on ectothermic metabolism, however, we formulated a new, non-exclusive hypothesis, termed "hibernation hypothesis": greater seasonality is expected to drive larger body size, since metabolic rate is reduced further and longer during colder, longer winters, leading to decreased energy depletion during hibernation, improved survival and increased longevity (and hence growth). Supporting this, a higher post-hibernation body condition in animals from areas of greater seasonality was found. CONCLUSIONS: Despite largely recognized effects of temperature on metabolic rate in ectotherms, its importance in determining body size in a gradient of seasonality has been largely overlooked so far. Based on our results, we present and discuss an alternative mechanism, the "hibernation hypothesis", underlying geographical body size variation, which can be helpful to improve our understanding of biogeographical patterns in ectotherms.

Is chytridiomycosis driving Darwin's frogs to extinction?

Darwin's frogs (Rhinoderma darwinii and R. rufum) are two species of mouth brooding frogs from Chile and Argentina that have experienced marked population declines. Rhinoderma rufum has not been found in the wild since 1980. We investigated historical and current evidence of Batrachochytrium dendrobatidis (Bd) infection in Rhinoderma spp. to determine whether chytridiomycosis is implicated in the population declines of these species. Archived and live specimens of Rhinoderma spp., sympatric amphibians and amphibians at sites where Rhinoderma sp. had recently gone extinct were examined for Bd infection using quantitative real-time PCR. Six (0.9%) of 662 archived anurans tested positive for Bd (4/289 R. darwinii; 1/266 R. rufum and 1/107 other anurans), all of which had been collected between 1970 and 1978. An overall Bd-infection prevalence of 12.5% was obtained from 797 swabs taken from 369 extant individuals of R. darwinii and 428 individuals representing 18 other species of anurans found at sites with current and recent presence of the two Rhinoderma species. In extant R. darwinii, Bd-infection prevalence (1.9%) was significantly lower than that found in other anurans (7.3%). The prevalence of infection (30%) in other amphibian species was significantly higher in sites where either Rhinoderma spp. had become extinct or was experiencing severe population declines than in sites where there had been no apparent decline (3.0%; x(2) = 106.407, P<0.001). This is the first report of widespread Bd presence in Chile and our results are consistent with Rhinoderma spp. declines being due to Bd infection, although additional field and laboratory investigations are required to investigate this further.

The population decline and extinction of Darwin's frogs.

Darwin's frogs (Rhinoderma darwinii and R. rufum) are two species of mouth-brooding frogs from Chile and Argentina. Here, we present evidence on the extent of declines, current distribution and conservation status of Rhinoderma spp.; including information on abundance, habitat and threats to extant Darwin's frog populations. All known archived Rhinoderma specimens were examined in museums in North America, Europe and South America. Extensive surveys were carried out throughout the historical ranges of R. rufum and R. darwinii from 2008 to 2012. Literature review and location data of 2,244 archived specimens were used to develop historical distribution maps for Rhinoderma spp. Based on records of sightings, optimal linear estimation was used to estimate whether R. rufum can be considered extinct. No extant R. rufum was found and our modelling inferred that this species became extinct in 1982 (95% CI, 1980-2000). Rhinoderma darwinii was found in 36 sites. All populations were within native forest and abundance was highest in Chiloé Island, when compared with Coast, Andes and South populations. Estimated population size and density (five populations) averaged 33.2 frogs/population (range, 10.2-56.3) and 14.9 frogs/100 m(2) (range, 5.3-74.1), respectively. Our results provide further evidence that R. rufum is extinct and indicate that R. darwinii has declined to a much greater degree than previously recognised. Although this species can still be found across a large part of its historical range, remaining populations are small and severely fragmented. Conservation efforts for R. darwinii should be stepped up and the species re-classified as Endangered.