Effects of amphibian phylogeny, climate and human impact on the occurrence of the amphibian-killing chytrid fungus

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.

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.

Spatial Risk Analysis of Batrachochytrium dendrobatidis, A Global Emerging Fungal Pathogen

Chytridiomycosis, a leading cause for the global decline in the number of amphibians, is caused by the fungal pathogen Batrachochytrium dendrobatidis. In this study, global distribution data of B. dendrobatidis were collected from January 2009 to May 2019. Space-time scan statistics and the maximum entropy (MaxEnt) model were used to analyze the epidemic trends and aggregation of the pathogen, and predict B. dendrobatidis distribution through its relationships with climate factors, wind speed, and solar radiation. The results of space-time scan statistics show seven clusters of data for the distribution of B. dendrobatidis. The time was mainly concentrated in 2009, 2013, 2015, and 2016, and the regions were primarily concentrated in southeastern Canada, southwestern France, Nigeria, Cameroon, eastern Brazil, southeastern Brazil, central Madagascar, and central and eastern Australia. MaxEnt showed that annual precipitation had the largest contribution percentage in the model, and annual mean temperature highly influenced the distribution of B. dendrobatidis. The global high-risk areas of B. dendrobatidis distribution were mainly observed in western Canada, southern Brazil, Chile, the United Kingdom, Japan, the Republic of Korea, eastern South Africa, eastern Madagascar, southeastern Australia, and southern China.

Chiggers (Acariformes: Trombiculoidea) do not increase rates of infection by Batrachochytrium dendrobatidis fungus in the endemic Dwarf Mexican Treefrog Tlalocohyla smithii (Anura: Hylidae)

Amphibian populations are globally declining at an alarming rate, and infectious diseases are among the main causes of their decline. Two micro-parasites, the fungus Batrachochytrium dendrobatidis (Bd) and the virus Ranavirus (RV) have caused mass mortality of amphibians and population declines. Other, less understood epizootics are caused by macro-parasites, such as Trombiculoidea chiggers. Infection with chiggers can affect frog behavior and survival. Furthermore, synergistic effects of co-infection with both macro and micro-parasites may lead to higher morbidity. To better understand these potential synergies, we investigated the presence and co-infection by chiggers, Bd and RV in the endemic frog Tlalocohyla smithii (T. smithii). Co-infection of Bd, RV, and/or chiggers is expected in habitats that are suitable for their co-occurrence; and if infection with one parasite facilitates infection with the others. On the other hand, co-infection could decrease if these parasites were to differ in their micro-environmental requirements (i.e. niche apportionment). A total of 116 frogs of T. smithii were studied during 2014 and 2016 in three streams within the Chamela-Cuixmala Biosphere Reserve in Jalisco, Mexico. Our results show that 31% of the frogs were infected with Trombiculoidea chiggers (Hannemania sp. and Eutrombicula alfreddugesi); Hannemania prevalence increased with air temperature and decreased in sites with high canopies and with water pH values above 8.5 and below 6.7. Bd prevalence was 2.6%, RV prevalence was 0%, and none of the frogs infected with chiggers were co-infected with Bd. Together, this study suggests that chiggers do not facilitate infection with Bd, as these are apportioned in different micro-habitats. Nevertheless, the statistical power to assure this is low. We recommend further epidemiological monitoring of multiple parasites in different geographical locations in order to provide insight on the true hazards, risks and conservation options for amphibian populations.

Amphibian fungal panzootic causes catastrophic and ongoing loss of biodiversity

Anthropogenic trade and development have broken down dispersal barriers, facilitating the spread of diseases that threaten Earth's biodiversity. We present a global, quantitative assessment of the amphibian chytridiomycosis panzootic, one of the most impactful examples of disease spread, and demonstrate its role in the decline of at least 501 amphibian species over the past half-century, including 90 presumed extinctions. The effects of chytridiomycosis have been greatest in large-bodied, range-restricted anurans in wet climates in the Americas and Australia. Declines peaked in the 1980s, and only 12% of declined species show signs of recovery, whereas 39% are experiencing ongoing decline. There is risk of further chytridiomycosis outbreaks in new areas. The chytridiomycosis panzootic represents the greatest recorded loss of biodiversity attributable to a disease.

How sensitive are temperate tadpoles to climate change? The use of thermal physiology and niche model tools to assess vulnerability

Ectotherms are vulnerable to climate change, given their dependence on temperature, and amphibians are particularly interesting because of their complex life cycle. Tadpoles may regulate their body temperature by using suitable thermal microhabitats. Thus, their physiological responses are the result of adjustment to the local thermal limits experienced in their ponds. We studied three anuran tadpole species present in Argentina and Chile: Pleurodema thaul and Pleurodema bufoninum that are seasonal and have broad geographic ranges, and Batrachyla taeniata, a geographically restricted species with overwintering tadpoles. Species with restricted distribution are more susceptible to climate change than species with broader distribution that may cope with potential climatic changes in the environments in which they occur. We aim to test whether these species can buffer the potential effects of climate warming. We used ecological niche models and the outcomes of their thermal attributes (critical thermal limits, optimal temperature, and locomotor performance breadth) as empirical evidence of their capacity. We found that Pleurodema species show broader performance curves, related to their occurrence, while the geographically restricted B. taeniata shows a narrower thermal breadth, but is faster in warmer conditions. The modeled distributions and empirical physiological results suggest no severe threats for these three anurans. However, the risk level is increasing and a retraction of their distribution range might be possible for Pleurodema species, and some local population extinctions may happen, particularly for the narrowly distributed B. taeniata.

Characterization of the Skin Microbiota of the Cane Toad Rhinella cf. marina in Puerto Rico and Costa Rica

Rhinella marina is a toad native to South America that has been introduced in the Antilles, likely carrying high loads of microorganisms, potentially impacting local community diversity. The amphibian skin is involved in pathogen defense and its microbiota has been relatively well studied, however, research focusing on the cane toad microbiota is lacking. We hypothesize that the skin microbial communities will differ between toads inhabiting different geographical regions in Central America and the Caribbean. To test our hypothesis, we compared the microbiota of three populations of R. cf. marina toads, two from Costa Rican (native) and one Puerto Rican (exotic) locations. In Costa Rica, we collected 11 toads, 7 in Sarapiquí and 4 from Turrialba while in Puerto Rico, 10 animals were collected in Santa Ana. Separate swab samples were collected from the dorsal and ventral sites resulting in 42 samples. We found significant differences in the structure of the microbial communities between Puerto Rico and Costa Rica. We detected as much as 35 different phyla; however, communities were dominated by Proteobacteria, Bacteroidetes, Firmicutes, and Actinobacteria. Alpha diversity and richness were significantly higher in toads from Puerto Rico and betadiversity revealed significant differences between the microbiota samples from the two countries. At the genus level, we found in Santa Ana, Puerto Rico, a high dominance of Kokuria, Niabella, and Rhodobacteraceae, while in Costa Rica we found Halomonas and Pseudomonas in Sarapiquí, and Acinetobacter and Citrobacter in Turrialba. This is the first report of Niabella associated with the amphibian skin. The core microbiome represented 128 Operational Taxonomic Units (OTUs) mainly from five genera shared among all samples, which may represent the symbiotic Rhinella's skin. These results provide insights into the habitat-induced microbial changes facing this amphibian species. The differences in the microbial diversity in Puerto Rican toads compared to those in Costa Rica provide additional evidence of the geographically induced patterns in the amphibian skin microbiome, and highlight the importance of discussing the microbial tradeoffs in the colonization of new ecosystems.

Host niche may determine disease-driven extinction risk

The fungal pathogen Batrachochytrium dendrobatidis (Bd) drives declines and extinctions in amphibian communities. However, not all regions and species are equally affected. Here, we show that association with amphibian aquatic habitat types (bromeliad phytotelmata versus stream) across Central America results in the odds of being threatened by Bd being five times higher in stream microhabitats. This differential threat of Bd was supported in our study by a significantly lower prevalence of Bd in bromeliad-associated amphibian species compared to riparian species in Honduran cloud forests. Evidence that the bromeliad environment is less favorable for Bd transmission is exemplified by significantly less suitable physicochemical conditions and higher abundance of Bd-ingesting micro-eukaryotes present in bromeliad water. These factors may inhibit aquatic Bd zoospore survival and the development of an environmental reservoir of the pathogen. Bromeliad phytotelmata thus may act as environmental refuges from Bd, which contribute to protecting associated amphibian communities against chytridiomycosis-driven amphibian declines that threaten the nearby riparian communities.