Skip to main content

Low microbiome diversity in threatened amphibians from two biodiversity hotspots


Microbial diversity positively influences community resilience of the host microbiome. However, extinction risk factors such as habitat specialization, narrow environmental tolerances, and exposure to anthropogenic disturbance may homogenize host-associated microbial communities critical for stress responses including disease defense. In a dataset containing 43 threatened and 90 non-threatened amphibian species across two biodiversity hotspots (Brazil’s Atlantic Forest and Madagascar), we found that threatened host species carried lower skin bacterial diversity, after accounting for key environmental and host factors. The consistency of our findings across continents suggests the broad scale at which low bacteriome diversity may compromise pathogen defenses in species already burdened with the threat of extinction.


Biodiversity loss is a major engine of ecosystem change in the Anthropocene [1]. The holobiont concept, in which the host and its microbiota function as a single organism, broadens the scales at which organisms and extinction drivers interact. The host microbiome is essential to immunity, development, metabolism, and stress responses, deeming it vitally interlinked with threatened species conservation. As a result, the targets of conservation biology have been expanding beyond preserving taxonomic and genetic diversity of species to include host-associated microbial diversity [2,3,4].

Species richness of host-associated bacterial assemblages (bacteriome) can exceed optimal levels, relative to host taxon and ecological context, if the microbiome is overwhelmed with transient or opportunistically pathogenic taxa [5]. However, in general, high bacterial diversity facilitates bacteriome stability, resilience, and function through mechanisms such as functional redundancy and synergy [6]. For instance, dramatic increases in metabolic, immune, and cognitive diseases since the twentieth century have been attributed to a consistent decline in human-associated microbiota diversity linked with industrialization [7].

Species more vulnerable to extinction may share certain biological traits, including habitat specialization, narrow environmental tolerances and geographic distributions, poor dispersal ability, and low genetic diversity [8]. These host traits, along with the habitat disturbances fueling species endangerment, may radically influence host bacteriomes, potentially skewing the bacteriome toward less diverse communities that may be less likely to support the symbiotic relationships and functions critical to host fitness and survival under the stresses of global change. However, the links between host endangerment and microbial diversity are unresolved across the host tree of life [9].

Amphibians stand apart from other animal groups as the most threatened taxon, with an estimated 41% of species listed in IUCN Red List threat categories, and are disproportionately impacted by emerging fungal, viral and protozoan diseases [10]. As in other taxa, amphibian bacteriome diversity often contributes to defenses against invading pathogens [11]. We tested the association between threat status and the taxonomic and phylogenetic diversity of skin-associated bacteria in 133 amphibian species (1454 samples) in two biodiversity hotspots nearly 10,000 km apart (Brazil’s Atlantic Forest and Madagascar). We compiled and analyzed data from 43 species listed under IUCN Red List threat categories (Brazil: 49 samples from nine species; Madagascar: 301 samples from 34 species) and 90 co-occurring non-threatened species (Brazil: 354 samples from 24 species; Madagascar: 750 samples from 66 species) across 63 localities (Brazil: 22 localities; Madagascar: 41 localities), defined as sampling areas less than 50 m in radius with consistent vegetation cover and microclimatic conditions. We targeted samples collected in natural forest habitat to control for direct effects of landcover. To standardize spatial variables, we averaged bacterial diversity metrics for each species within each sampling locality (Brazil: n = 56 species-localities; Madagascar: n = 278 species-localities; mean sample size per species-locality = 5.4 [Brazil/threatened], 7.5 [Brazil/non-threatened], 4.6 [Madagascar/threatened], 3.5 [Madagascar/non-threatened]). Threatened and non-threatened species spanned primary habitat use categories (aquatic, arboreal, terrestrial; see Methods) and were evenly distributed across the latitudinal sampling extent (Additional file 1: Fig. S1). We predicted decreasing skin bacteriome diversity with increasing threat status, based on risk factors for extinction such as specialization and narrow thermal breadth, and population characteristics of threatened species, including low genetic variation [12].

Our study reveals a cross-continental pattern of lower bacteriome diversity in the most threatened species within the most vulnerable animal taxon. Using a dual analysis approach including piecewise structural equation modeling and phylogenetic path analysis, we found that taxonomic diversity of host skin-associated bacteria (richness of sub-operational taxonomic units [sOTUs]) was negatively correlated with host threat status among amphibian communities in Brazil’s Atlantic Forest (Figs. 1A, C; Additional file 1: Table S1A) and Madagascar (Figs. 1B, D; Additional file 1: Table S1B). This correlation was robust to variation in climate, sample DNA extraction method (Fig. 1; Additional file 1: Table S1), host body length, vegetation density, and host geographic range area (Additional file 1: Table S2), and remained unaltered after accounting for host phylogeny (Additional file 1: Table S3) and host primary habitat type (Additional file 1: Fig. S2; see interaction term in Additional file 1: Table S4A). We found an identical pattern of lower Faith’s phylogenetic skin bacterial diversity in threatened compared to non-threatened species (Additional file 1: Fig. S3; Table S4B). In both focal geographic regions, threatened and non-threatened species carried differentially abundant skin bacterial taxa (Additional file 1: Figs. S4 and S5, Table S5), including higher abundances of the fungal pathogen-fighting skin bacterium Janthinobacterium lividum carried by non-threatened amphibian species in Brazil. Determining the drivers of these patterns will require disentangling host physiology, ecology and biogeography alongside the biotic and abiotic environmental stressors differentially affecting threatened and non-threatened species.

Fig. 1
figure 1

Association between skin bacterial diversity and threat status of amphibian species in two geographically distinct biodiversity hotspots. A, B Piecewise structural equation models accounting for environmental, host, and methodological factors. Amphibian skin bacterial diversity was estimated as number of detected sub-operational taxonomic units (sOTUs). Numbers are standardized coefficients (*p < 0.05). Unsupported paths shown as gray arrows were removed to improve model fit. Red arrows indicate the correlation of interest between host skin bacterial diversity and host threat status. Black arrows indicate other paths that were retained in the final pruned models. C, D Average skin sOTU richness between threatened (red) and non-threatened (gray) species. Threatened species carried lower skin sOTU richness in (C) Brazil’s Atlantic Forest (t = − 2.407, df = 15.436, p = 0.029) and D Madagascar (t = − 2.894, df = 131.981, p = 0.004). Error bars represent standard error

Ecological specialization may predispose species to endangerment. Many of the threatened species included in our study have relatively small distributions and specialized behaviors and physiologies. For example, the endangered treefrog Xenohyla truncata, endemic to coastal scrub forests in the Brazilian state of Rio de Janeiro, feeds on fruit (a unique trait among amphibians), in addition to small arthropods, and shelters in bromeliads [13]. Critically endangered toadlets Melanophryniscus admirabilis are confined to mossy rock pools along one river in southern Brazil [14]. The critically endangered Brazilian treefrog Nyctimantis pomba is found exclusively in bamboo groves within a single forest fragment [15]. In Madagascar, the microendemic ecological specialists Gephyromantis corvus and G. kintana are restricted to the Isalo sandstone massif, and the threatened montane specialists Boophis laurenti, Mantidactylus pauliani, M. madecassus, and Anodonthyla montana are each restricted to extremely narrow geographic ranges on the island's highest mountain peaks [16]. Together, these highly specialized host natural histories suggest high fidelity to specific microhabitats, dietary constraints, low dispersal capabilities and, in mild climates, narrow physiological tolerances that may limit the diversity of bacteria encountered in the environment and subsequently recruited to the skin microbiome.

Geographic range area was a poor predictor of skin bacterial diversity, potentially a reflection of our individual-level sampling scale. Compared to wide-ranging non-threatened species, threatened species with inherently small or shrinking geographic distributions may encounter smaller and potentially less diverse environmental source pools from which host microbiomes are seeded across the entire species range. However, our study controlled for direct spatial effects by sampling at a narrow, local scale (locality). Limiting our sampling to habitats with natural vegetation likewise accounted for direct effects of disturbance on host microbiomes and environmental bacterial source pools, such as microbiome depletion with increasing deforestation [17]. Indeed, threatened amphibian species are often only found in the most intact habitats [10] where environmental bacterial pools should be highly diverse [17], drawing an even stronger contrast between observed skin bacterial diversity of threatened species and expected diversity based on meso- or macro-scale habitat characteristics.

In the case of some endemic species with restricted geographic distributions, however, encroaching anthropogenic activity could impose physiological stress that may account for lower skin bacterial diversity, such as microclimatic edge effects or environmental contamination. For instance, the critically endangered Physalaemus soaresi is restricted to a single natural forest fragment less than 5 km2 in area within the suburbs of Rio de Janeiro [18]. Thus, factors associated with environmental degradation such as pollution, altered food web dynamics, and climatic shifts from nearby development could generate physiological stress through impacts on energetics, food quality and availability, behavior, and environmental microbial pools, with potential downstream shifts in skin microbiome composition [19]. Host population bottlenecks commonly experienced by threatened species in stress-generating environments may also constrain bacteriome diversity through sampling effects from generation to generation [20].

We found that threatened amphibian species from two continents carried lower skin bacterial diversity than non-threatened species. This pattern was consistent across all three primary host habitat use categories (arboreal, aquatic, terrestrial) and was robust to climate, sample DNA extraction method, host body length, vegetation density, host geographic range area, and host phylogeny. Regardless of the cause-and-effect associations between low bacteriome diversity and endangerment, this pattern has implications for health and fitness of threatened vertebrates. Diversity of the bacteriome may enhance resilience to invasion, environmental change, and other stressors through functional redundancy, partitioning use of limiting resources, or production of antimicrobial metabolites [21]. Another finding from the Brazil dataset that deserves further attention is that non-threatened amphibian species not only carried higher skin bacterial diversity than threatened species but also carried higher abundances of Janthinobacterium lividum. This amphibian skin bacterium has known inhibitory function against the causative agent of the amphibian fungal disease chytridiomycosis and has been successfully employed as a probiotic in amphibians susceptible to this disease [22]. Given that species are at increased risk of disease and other pressures as they move toward extinction [23], our findings suggest that the combined threats of low microbial diversity and vulnerability to pathogens may compound extinction risk, especially considering the slow pace of bacteriome recovery following pathogenic infections or other environmental disturbances [24, 25]. Our results raise a red flag for amphibians and threatened taxa from other animal groups and highlight the need for mechanistic insights to explain the geographically consistent and statistically robust inverse association between host threat status and bacteriome diversity.

There is some evidence that amphibian species driven to endangerment by disease can develop resistant microbiomes in response to pathogen pressure [26]. Alternatively, microbiome manipulation can mimic these processes interventionally [27]. However, whether naturally low bacteriome diversity hampers adaptive immune processes or efficacy of probiotic strategies is unresolved. Also, critical to threatened species conservation is the intersection of the host bacteriome and increasingly extreme temperature fluctuations associated with anthropogenic climate change, one of the most serious threats to listed species [28]. Our findings bring to light the urgency in characterizing microbial baselines for threatened species, not only as reference points for probiotics, but also to maximize success of wildlife ‘ark’ programs. Within these programs, bacteriome surveillance can guide captive husbandry protocols, antibiotic use, early detection of disease outbreaks and other environmental disturbances, translocation and reintroduction programs, and can be used as a metric for assessing efficacy of habitat restoration.


Study regions and sampling design

We conducted this study with data from two tropical biodiversity hotspots: Brazil’s Atlantic Forest and Madagascar (Additional file 2). The Brazil database comprised sequence data from nine previously published [29,30,31,32,33] and unpublished studies (SRA Bioproject accession PRJNA767814). This database consisted of samples from 403 individual post-metamorphic anurans representing 33 species (9 threatened; 24 non-threatened) and eight families from 22 localities (Additional file 1: Fig. S1). The Madagascar database comprised sequence data from one previously published study [1634] (SRA Bioproject accession PRJNA394790). This database consisted of samples from 1,051 individual post-metamorphic anurans representing 100 species (34 threatened; 66 non-threatened) and four families from 41 localities (Additional file 1: Fig. S1).

Both countries were sampled across a broad geographic area to capture a wide breadth of amphibian species diversity, with sampling localities haphazardly distributed (Additional file 1: Fig. S1). Threatened and non-threatened species were evenly distributed across the latitudinal sampling extent and frequently co-occurred within localities (Additional file 1: Fig. S1). We further ruled out spatial bias in our sampling by estimating potential differences in spatial autocorrelation of sOTU richness across distance bands for threatened versus non-threatened amphibian species using Moran’s I correlograms (Additional file 1: Fig. S6). Threatened and non-threatened species both spanned primary habitat use categories (species per category = 2 terrestrial, 3 arboreal, 4 aquatic [Brazil/threatened]; 8 terrestrial, 7 arboreal, 9 aquatic [Brazil/non-threatened]; 16 terrestrial, 9 arboreal, 9 aquatic [Madagascar/threatened]; 22 terrestrial, 32 arboreal, 12 aquatic [Madagascar/non-threatened]).

Sampling protocol

Anurans were captured in the wild in localities with primarily natural vegetation, handled using disposable gloves, rinsed with sterile water to remove transient microbes and debris, and the skin surface was swabbed using sterile swabs following a standard protocol [35]. Swabs were kept on ice in the field and stored frozen at − 20 °C until further processing. Samples were primarily collected during the Austral breeding season from September to March, with a small proportion of samples collected in northern Brazil during the boreal breeding season in June.

Bacterial sequencing and bioinformatics

Bacterial DNA was extracted from swabs from Brazil using the Qiagen DNeasy Blood and Tissue Kit or Prepman Ultra. DNA was extracted from swabs from Madagascar using the Qiagen DNeasy PowerSoil kit. DNA was amplified, purified, and sequenced on Illumina MiSeq sequencing platforms (2 × 250 or 2 × 150) following the Earth Microbiome Project 16S Illumina Amplicon Protocol. This protocol targets the V4 region of the 16S rRNA gene from the bacterial genome using barcoded primers 515F and 806R and a dual index approach.

The Brazil dataset was processed similarly to the previously published Madagascar dataset [1634]. Sequences were initially processed with Quantitative Insights into Microbial Ecology (Brazil: QIIME 2 version 2019.1 [36]; Madagascar: QIIME [34]). To maximize read quality [37], only forward reads trimmed to 150 base pairs were used for both datasets. Sequences were filtered based on quality score, using the q-score plugin for Brazil sequences and analogous criteria for Madagascar sequences [34]. Quality-filtered sequences were clustered into sub-operational taxonomic units (sOTUs) using the Deblur workflow [38]. Phylogenetic trees were constructed using FastTree and taxonomy was assigned with the classify‐sklearn naïve Bayes taxonomy classifier and the Greengenes 13.8 reference sequence database (Brazil) or the Ribosomal Database Project Classifier 5 with a custom script (Madagascar) [34]. Rare sOTUs were filtered by discarding sOTUs with less than 0.001% (53 reads) of total sequence reads across the dataset (Brazil) or less than 10 sequence reads across the dataset (Madagascar). We used conservative thresholds for removal of rare sOTUs because our datasets spanned large geographic areas and many host species, and thus a large proportion of sOTUS were relatively rare. For the Brazil samples, we discarded one sOTU that was considered a contaminant because it was abundant across control and template DNA samples. For the Madagascar samples, potential contaminants were identified with the script in QIIME and discarded [34].

After filtering and decontamination, the Brazil dataset contained 4,823,452 reads and an average of 11,969 reads per sample and the Madagascar dataset contained 18,084,933 reads and an average of 17,207 reads per sample. To normalize read counts across samples, samples were rarefied to 1500 (Brazil) or 2500 (Madagascar) and averaged 154 (Brazil) and 110 (Madagascar) sOTUs per sample [39]. We examined rarefaction curves to ensure sufficient sequence depth (Additional file 1: Figs. S7 and S8).

To ensure that our data were not biased by using only forward sequence reads, we repeated the Brazil bioinformatics using both forward and reverse sequence reads (QIIME 2 version 2021.8). Forward and reverse reads were joined for compatibility with the Deblur workflow [38]. We used the same filtering (0.001% of total sequence reads across dataset = 40) and decontamination methods as described for the Brazil dataset above. After filtering and decontamination, this dataset contained 3,632,608 reads and an average of 7762 reads per sample. To normalize read counts across samples, samples were rarefied to 500 and averaged 77 sOTUs per sample [39]. The lower rarefaction threshold (500) compared to the threshold used with only forward reads (1500) reflects consistently lower numbers of sequence reads detected per sample when using joined forward and reverse reads, likely due to lower quality of reverse reads. However, we examined rarefaction curves to ensure that this rarefaction level captured sufficient sequence depth (Additional file 1: Fig. S9). After rarefaction, the dataset contained 401 samples (49 threatened; 352 non-threatened), a similar sample size to the 403 samples (49 threatened; 354 non-threatened) analyzed when using only forward sequence reads.

Statistical analysis

We analyzed the Brazil and Madagascar datasets separately. We used piecewise structural equation models (psem function in the piecewiseSEM package in Program R) to test for correlations between IUCN threat status and sOTU richness while accounting for direct and indirect associations with key environmental and host factors [40, 41]. We averaged sOTU richness for each species within each sampling locality (Brazil: n = 56 species-localities; Madagascar: n = 278 species-localities; mean sample size per species-locality = 5.4 [Brazil/threatened], 7.5 [Brazil/non-threatened], 4.6 [Madagascar/threatened], 3.5 [Madagascar/non-threatened]). We classified species into threat categories using a Brazil-specific Red List [42] and the International Union for Conservation of Nature’s (IUCN) Red List of Threatened Species for Madagascar species [10]. We classified threat status as a binary variable with Least Concern species assigned a value of 0 and threatened (Vulnerable, Endangered, and Critically Endangered) and Near-Threatened (NT) species assigned a value of 1 [10]. We grouped threatened and NT species because almost all NT species are reported to be impacted by land use change [10, 43, 44].

Environmental covariates included annual mean temperature (BIO1) and annual precipitation (BIO12), extracted for each sampling locality from WorldClim at 5-min spatial resolution, and Normalized Difference Vegetation Index (NDVI), averaged over the sampling period for each locality. Host species covariates included body length (maximum male snout-vent length [SVL])[45] and geographic range area (log-transformed). For Brazil, we included DNA extraction method as an additional predictor of sOTU richness to account for variation attributable to differences in DNA extraction protocol among samples. Preliminary generalized linear models (Poisson error distribution, log link) revealed that NDVI and geographic range area were consistently poor predictors of sOTU richness for both the Brazil and Madagascar datasets and were excluded from SEMs (Additional file 1: Table S2). Results of preliminary models remained unaltered as standard least squares general linear models, so we report the more conservative Poisson models. For the SEMs, we first ran saturated models with all ecologically relevant paths and then ran simplified models with unsupported variables excluded to improve model fit based on Akaike’s Information Criterion (AIC; Brazil: saturated model AIC = 21.591, simplified model AIC = 15.108; Madagascar: saturated model AIC = 33.550, simplified model AIC = 14.612). To ensure data were not biased by using only forward sequence reads, we repeated SEMs for Brazil using joined forward and reverse sequence data and our results remained unaltered (Additional file 1: Table S6).

To determine if SEM results were robust to effects of possible non-independence due to shared trends of common ancestry among frog species, we also tested the association between skin bacterial richness and threat status with phylogenetic path analysis (PPA) under the d-separation method [46] as implemented in the R package phylopath [47]. Each causal model was formulated in terms of separate directed acyclic graphs where conditional independencies (i.e. d-separation statements) were translated into phylogenetic linear regressions for analysis using phylolm [48]. This R package allows for measuring the degree of phylogenetic signal embedded in the data when the causal parents of the models are either continuous (λ) [49, 50] or discrete (α) [51]. In contrast to SEM, which accommodates a bidirectional link between skin bacterial richness and threat status, here we modeled averaged path coefficients shared by the alternative directed acyclic graphs, weighing them under the C-statistic information criterion corrected for small sample sizes (CICc) which provides a measure of the strength of evidence for each whole path model [52].

The evolutionary correlations for PPA were accounted for by using the phylogenies from Hedges et al. [53] for the Brazil dataset, and Bletz et al. [16] for the Madagascar dataset. The former [53] corresponds to a time-calibrated amphibian tree synthesized from studies in molecular evolution and phylogenetics [54, 55], and pruned to include only our sampled Brazilian species. The latter [16] corresponds to the ultrametric time-tree computed with optimized branch lengths under the GTR model from 16S rRNA sequences. For the focal taxa, we compiled sequences of the 16S rRNA gene for all included taxa of Madagascar frogs, and aligned them with MAFFT v. 7 [56]. We then computed an initial phylogenetic tree under Maximum Likelihood and the general time reversible (GTR) substitution model in MEGA7 [57]. We then manually adjusted the topology of this tree to fit the most recent multi-gene phylogenetic trees available for subsets of these taxa, and for their deep relationships: Scherz et al. [58] for microhylids, Wollenberg et al. [59] for Mantidactylus, Kaffenberger et al. [60] for Gephyromantis, and Hutter et al. [61] for Boophis. We then optimized branch lengths of this user tree in MEGA7 under the GTR model. We also used the same program to compute an ultrametric time-tree from the same data set, using the RELTIME approach.

The negative correlation between threat status and skin bacterial diversity persisted after accounting for host phylogeny (Additional file 1: Table S3). The phylogenetic approach favored threat status as a predictor of bacterial diversity over the reverse. The phylogenetic analysis was consistent with the SEMs, except for one discrepancy for Madagascar. While SEM detected that abiotic effects on the correlation between threat status and skin bacterial diversity were primarily mediated through effects of temperature on threat status, PPA detected stronger effects of temperature on skin bacterial diversity.

We used general linear models as an additional tool to verify that known host factors did not influence the correlation between host threat status and skin bacterial diversity. For these models, we included host skin bacterial diversity (sOTU richness or Faith’s phylogenetic diversity) as the response and the following predictors: host primary habitat type (arboreal, aquatic, or terrestrial), threat status (threatened or non-threatened), country (Brazil or Madagasar), and the interaction between host primary habitat type and threat status.

We used linear discriminant analysis effect size (LEfSe) on the Galaxy platform to detect differentially abundant sOTUs between threatened and non-threatened species in each diversity hotspot [62, 63]. We used default parameters except increasing the threshold on the logarithmic LDA score (LogLDA) from 2.0 to 3.0 in order to remove weaker correlations. To visualize the results of this analysis we constructed heat maps using the heatmap.2 function in R (Heatplus package, R version 4.2.0) for Brazil and Madagascar data separately [64].

Availability of data and materials

The datasets generated and analyzed during the current study are available in the SRA repository (Brazil: Bioproject accession PRJNA767814; Madagascar: Bioproject accession PRJNA394790).



Akaike’s information criterion

CICc :

C-statistic information criterion corrected for small sample sizes


General time reversible


International union for conservation of Nature


Linear discriminant analysis effect size


Normalized difference vegetation index




Phylogenetic path analysis


Structural equation model


Sub-operational taxonomic unit


Snout-vent length


  1. Hooper DU, Adair EC, Cardinale BJ, Byrnes JEK, Hungate BA, Matulich KL, et al. A global synthesis reveals biodiversity loss as a major driver of ecosystem change. Nature. 2012;486:105–8.

    Article  CAS  Google Scholar 

  2. Trevelline BK, Fontaine SS, Hartup BK, Kohl KD. Conservation biology needs a microbial renaissance: a call for the consideration of host-associated microbiota in wildlife management practices. Proc Royal Soc B Biol Sci. 2019;286.

  3. Banerjee A, Cornejo J, Bandopadhyay R. Emergent climate change impact throughout the world: call for “microbiome conservation” before it’s too late. Biodivers Conserv. 2020;29:345–8.

    Article  Google Scholar 

  4. Redford KH, Segre JA, Salafsky N, Del Rio CM, Mcaloose D. Conservation and the microbiome. Conserv Biol. 2012;26:195–7.

    Article  Google Scholar 

  5. Greenspan SE, Lyra ML, Migliorini GH, Kersch-Becker MF, Bletz MC, Lisboa CS, et al. Arthropod-bacteria interactions influence assembly of aquatic host microbiome and pathogen defense. Proc Royal Soc B Biol Sci. 2019;286.

  6. Fassarella M, Blaak EE, Penders J, Nauta A, Smidt H, Zoetendal EG. Gut microbiome stability and resilience: elucidating the response to perturbations in order to modulate gut health. Gut. 2021;70:595–605.

    Article  CAS  Google Scholar 

  7. Dominguez Bello MG, Knight R, Gilbert JA, Blaser MJ. Preserving microbial diversity. Science. 2018;362:33–5.

    Article  CAS  Google Scholar 

  8. Foden WB, Butchart SHM, Stuart SN, Vié JC, Akçakaya HR, Angulo A, et al. Identifying the world’s most climate change vulnerable species: a systematic trait-based assessment of all birds, amphibians and corals. PLoS ONE. 2013;8: e65427.

    Article  CAS  Google Scholar 

  9. Peixoto RS, Voolstra CR, Sweet M, Duarte CM, Carvalho S, Villela H, et al. Harnessing the microbiome to prevent global biodiversity loss. Nat Microbiol. 2022.

    Article  Google Scholar 

  10. IUCN. IUCN Red List of Threatened Species. Version 2021–1. 2021. Accessed 11 May 2021.

  11. Bernardo-Cravo AP, Schmeller DS, Chatzinotas A, Vredenburg VT, Loyau A. Environmental factors and host microbiomes shape host–pathogen dynamics. Trends Parasitol. 2020;36:616–33.

    Article  CAS  Google Scholar 

  12. Colles A, Liow LH, Prinzing A. Are specialists at risk under environmental change? Neoecological, paleoecological and phylogenetic approaches. Ecol Lett. 2009;12:849–63.

    Article  Google Scholar 

  13. Da Silva HR, De Britto-Pereira MC. How much fruit do fruit-eating frogs eat? An investigation on the diet of Xenohyla truncata (Lissamphibia: Anura: Hylidae). J Zool. 2006;270:692–8.

    Article  Google Scholar 

  14. Di-Bernardo M, Maneyro R, Grillo H. New species of Melanophryniscus (Anura: Bufonidae) from Rio Grande do Sul, southern Brazil. J Herpetol. 2006;40:261–6.

    Article  Google Scholar 

  15. De Assis CL, Santana DJ, Da Silva FA, Quintela FM, Feio RN. A new and possibly critically endangered species of casque-headed tree frog Aparasphenodon Miranda-Ribeiro, 1920 (Anura, Hylidae) from southeastern Brazil. Zootaxa. 2013;3716:583–91.

    Article  Google Scholar 

  16. Bletz MC, Archer H, Harris RN, McKenzie VJ, Rabemananjara FCE, Rakotoarison A, et al. Host ecology rather than host phylogeny drives amphibian skin microbial community structure in the biodiversity hotspot of Madagascar. Front Microbiol. 2017;8:1530.

    Article  Google Scholar 

  17. Becker CG, Longo AV, Haddad CFB, Zamudio KR. Land cover and forest connectivity alter the interactions among host, pathogen and skin microbiome. Proc Royal Soc B Biol Sci. 2017;284:20170582.

    Article  Google Scholar 

  18. Silvano DL, Segalla MV. Conservation of Brazilian amphibians. Conserv Biol. 2005;19:653–8.

    Article  Google Scholar 

  19. Stothart MR, Bobbie CB, Schulte-Hostedde AI, Boonstra R, Palme R, Mykytczuk NCS, et al. Stress and the microbiome: linking glucocorticoids to bacterial community dynamics in wild red squirrels. Biol Let. 2016;12:20150875.

    Article  Google Scholar 

  20. Ørsted M, Yashiro E, Hoffmann AA, Kristensen TN. Population bottlenecks constrain host microbiome diversity and genetic variation impeding fitness. PLoS Genet. 2022;18: e1010206.

    Article  Google Scholar 

  21. Lozupone CA, Stombaugh JI, Gordon JI, Jansson JK, Knight R. Diversity, stability and resilience of the human gut microbiota. Nature. 2012;489:220–30.

    Article  CAS  Google Scholar 

  22. Harris RN, Brucker RM, Walke JB, Becker MH, Schwantes CR, Flaherty DC, et al. Skin microbes on frogs prevent morbidity and mortality caused by a lethal skin fungus. ISME J. 2009;3:818–24.

    Article  CAS  Google Scholar 

  23. Heard MJ, Smith KF, Ripp K, Berger M, Chen J, Dittmeier J, et al. Increased threat of disease as species move towards extinction. Conserv Biol. 2013;27:1378–88.

    Article  Google Scholar 

  24. Jani AJ, Bushell J, Arisdakessian CG, Belcaid M, Boiano DM, Brown C, et al. The amphibian microbiome exhibits poor resilience following pathogen-induced disturbance. ISME J. 2021;15:1628–40.

    Article  CAS  Google Scholar 

  25. Walke JB, Becker MH, Krinos A, Chang EAB, Santiago C, Umile TP, et al. Seasonal changes and the unexpected impact of environmental disturbance on skin bacteria of individual amphibians in a natural habitat. FEMS Microbiol Ecol. 2021;97:fiaa248.

  26. Jiménez RR, Alvarado G, Estrella J, Sommer S. Moving beyond the host: Unraveling the skin microbiome of endangered Costa Rican amphibians. Front Microbiol. 2019;10:2060.

    Article  Google Scholar 

  27. Jin Song S, Woodhams DC, Martino C, Allaband C, Mu A, Javorschi-Miller-Montgomery S, et al. Engineering the microbiome for animal health and conservation. Exp Biol Med. 2019;244:494–504.

    Article  CAS  Google Scholar 

  28. Delach A, Caldas A, Edson KM, Krehbiel R, Murray S, Theoharides KA, et al. Agency plans are inadequate to conserve US endangered species under climate change. Nat Clim Chang. 2019;9:999–1004.

    Article  Google Scholar 

  29. Mesquita AFC, Lambertini C, Lyra M, Malagoli LR, James TY, Toledo LF, et al. Low resistance to chytridiomycosis in direct-developing amphibians. Sci Rep. 2017;7:16605.

    Article  Google Scholar 

  30. Moura-Campos D, Greenspan SE, DiRenzo GV, Neely WJ, Toledo LF, Becker CG. Fungal disease cluster in tropical terrestrial frogs predicted by low rainfall. Biol Cons. 2021;261: 109246.

    Article  Google Scholar 

  31. Preuss JF, Greenspan SE, Rossi EM, Lucas Gonsales EM, Neely WJ, Valiati VH, et al. Widespread pig farming practice linked to shifts in skin microbiomes and disease in pond-breeding amphibians. Environ Sci Technol. 2020;54:11301–12.

    Article  CAS  Google Scholar 

  32. Neely WJ, Greenspan SE, Ribeiro LP, Carvalho T, Martins RA, Rodriguez D, et al. Synergistic effects of warming and disease linked to high mortality in cool-adapted terrestrial frogs. Biol Cons. 2019;2020(245): 108521.

    Google Scholar 

  33. Martins RA, Greenspan SE, Medina D, Buttimer S, Marshall VM, Neely WJ, et al. Signatures of functional bacteriome structure in a tropical direct-developing amphibian species. Animal Microbiome. 2022;4.

  34. Kueneman JG, Bletz MC, McKenzie VJ, Becker CG, Joseph MB, Abarca JG, et al. Community richness of amphibian skin bacteria correlates with bioclimate at the global scale. Nat Ecol Evolut. 2019;3:381–9.

    Article  Google Scholar 

  35. Hyatt A, Boyle D, Olsen V, Boyle D, Berger L, Obendorf D, et al. Diagnostic assays and sampling protocols for the detection of Batrachochytrium dendrobatidis. Dis Aquat Org. 2007;73:175–92.

    Article  CAS  Google Scholar 

  36. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37:848–57.

    Google Scholar 

  37. Kwon S, Park S, Lee B, Yoon S. In-depth analysis of interrelation between quality scores and real errors in illumina reads. In: Proceedings of the annual international conference of the IEEE engineering in medicine and biology society, EMBS. 2013;635–8.

  38. Amir A, McDonald D, Navas-Molina JA, Kopylova E, Morton JT, Xu ZZ, et al. Deblur rapidly resolves single-nucleotide community sequence patterns. mSystems. 2017;2:e00191–16.

  39. Hong J, Karaoz U, De Valpine P, Fithian W. To rarefy or not to rarefy: Robustness and efficiency trade-offs of rarefying microbiome data. Bioinformatics. 2022;38:2389–96.

    Article  CAS  Google Scholar 

  40. Lefcheck JS. piecewiseSEM: piecewise structural equation modelling in R for ecology, evolution, and systematics. Methods Ecol Evol. 2016;7:573–9.

    Article  Google Scholar 

  41. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria. URL 2018.

  42. Brasil. Portaria no. 444, de 17 de dezembro de 2014. Diário Oficial da União, 245(1). 2014.

  43. Wilcove DS, Rothstein D, Dubow J, Phillips A, Losos E. Quantifying threats to imperiled species in the United States: assessing the relative importance of habitat destruction, alien species, pollution, overexploitation, and disease. Bioscience. 1998;48:607–15.

    Article  Google Scholar 

  44. Vitousek PM, Mooney HA, Lubchenco J, Melillo JM. Human domination of Earth’s ecosystems. Science. 1997;277:494–9.

    Article  CAS  Google Scholar 

  45. Ruthsatz K, Lyra ML, Lambertini C, Belasen AM, Jenkinson TS, da Silva LD, et al. Skin microbiome correlates with bioclimate and Batrachochytrium dendrobatidis infection intensity in Brazil’s Atlantic Forest treefrogs. Sci Rep. 2020;10:22311.

    Article  CAS  Google Scholar 

  46. von Hardenberg A, Gonzalez-Voyer A. Disentangling evolutionary cause-effect relationships with phylogenetic confirmatory path analysis. Evolution. 2013;67:378–87.

    Article  Google Scholar 

  47. van der Bijl W. phylopath: easy phylogenetic path analysis in R. PeerJ. 2018;2018: e4718.

    Article  Google Scholar 

  48. Tung Ho LS, Ané C. A linear-time algorithm for gaussian and non-gaussian trait evolution models. Syst Biol. 2014;63:397–408.

    Article  Google Scholar 

  49. Freckleton RP, Harvey PH, Pagel M. Phylogenetic analysis and comparative data: a test and review of evidence. Am Nat. 2002;160:712–26.

    Article  CAS  Google Scholar 

  50. Pagel M. Inferring the historical patterns of biological evolution. Nature. 1999;401:877–84.

    Article  CAS  Google Scholar 

  51. Ives AR, Garland T. Phylogenetic logistic regression for binary dependent variables. Syst Biol. 2010;59:9–26.

    Article  Google Scholar 

  52. Gonzalez-Voyer A, von Hardenberg A. An introduction to phylogenetic path analysis. In: Garamszegi LS, editor. Modern phylogenetic comparative methods and their application in evolutionary biology. Berlin: Springer; 2014. p. 201–29.

    Chapter  Google Scholar 

  53. Hedges SB, Marin J, Suleski M, Paymer M, Kumar S. Tree of life reveals clock-like speciation and diversification. Mol Biol Evol. 2015;32:835–45.

    Article  CAS  Google Scholar 

  54. Hedges SB, Dudley J, Kumar S. TimeTree: a public knowledge-base of divergence times among organisms. Bioinformatics. 2006;22:2971–2.

    Article  CAS  Google Scholar 

  55. Hedges SB, Kumar S. The timetree of life. New York: Oxford University Press; 2009.

    Google Scholar 

  56. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.

    Article  CAS  Google Scholar 

  57. Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.

    Article  CAS  Google Scholar 

  58. Scherz MD, Vences M, Rakotoarison A, Andreone F, Köhler J, Glaw F, et al. Reconciling molecular phylogeny, morphological divergence and classification of Madagascan narrow-mouthed frogs (Amphibia: Microhylidae). Mol Phylogenet Evol. 2016;100:372–81.

    Article  Google Scholar 

  59. Wollenberg KC, Vieites DR, Glaw F, Vences M. Speciation in little: the role of range and body size in the diversification of Malagasy mantellid frogs. BMC Evol Biol. 2011;11:217.

    Article  Google Scholar 

  60. Kaffenberger N, Wollenberg KC, Köhler J, Glaw F, Vieites DR, Vences M. Molecular phylogeny and biogeography of Malagasy frogs of the genus Gephyromantis. Mol Phylogenet Evol. 2012;62:555–60.

    Article  Google Scholar 

  61. Hutter CR, Lambert SM, Andriampenomanana ZF, Glaw F, Vences M. Molecular phylogeny and diversification of Malagasy bright-eyed tree frogs (Mantellidae: Boophis). Mole Phylogenet Evol. 2018;127:568–78.

    Article  Google Scholar 

  62. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biology. 2011;12.

  63. Afgan E, Baker D, Batut B, van den Beek M, Bouvier D, Ech M. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucleic Acids. 2018;46:537–44.

    Article  Google Scholar 

  64. Ploner A. Heatplus: Heatmaps with row and/or column covariates and colored clusters. 2022. R package version 3.6.0,

Download references


Atlantic Forest Fieldwork was conducted through the “Documenting Threatened Species” project. We thank everyone who helped with fieldwork and lab work on the projects that contributed to this paper.


National Science Foundation (DEB-1947681) (CGB), (DEB-2225683) (JP). Fresno Chaffee Zoo Wildlife Conservation Fund (PP). Columbus Zoo Fund for Conservation (PP). Brazilian National Council for Scientific and Technological Development (CNPq #300984/2022-7) (PP) (CNPq #306623/2018-8) (CFBH), (CNPq #302834/2020-6) (LFT). São Paulo Research Foundation (FAPESP #2013/50741-7) (CFBH), (FAPESP #2014/16320-7) (CN) (FAPESP #2016/25358-3; #2019/18335-5) (LFT).

Author information

Authors and Affiliations



Conceptualization: SEG, CGB. Formal analysis: SEG, JAFG, WJN, CGB. Investigation: PP, MB, MLL, IFM, RAM, DM, DMC, WJN, JP, MJS, RIV, CAN, AMT, MV, DW, CFBH, JP, CGB. Funding acquisition: PP, LFT, CGB. Writing—original draft: SEG. Writing—review & editing: all co-authors. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Sasha E. Greenspan or C. Guilherme Becker.

Ethics declarations

Ethics approval and consent to participate

Animal swabbing procedures were approved by the Institutional Animal Care and Use Committees of the University of Alabama, James Madison University, Universidade Estadual de Campinas, Universidade Estadual Paulista, and Universidade Federal do Pará. Collecting permits were approved by authorities in both Brazil (SISBIO #17242-5, #22511-3, #29964, #57061-8, #67469-1, #67633, #70883-1) and Madagascar (N.182/13/MEF/SF/DGF/DCB.SAP/SCB and N.288/14/MEEF/SG/DGEF/DCB.SAP/SCB).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

 Supplemental figures and tables.

Additional file 2:

 List of samples and associated metadata.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Greenspan, S.E., Peloso, P., Fuentes-González, J.A. et al. Low microbiome diversity in threatened amphibians from two biodiversity hotspots. anim microbiome 4, 69 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Anuran
  • Endangered species
  • Host-associated microbial diversity
  • Brazil
  • Madagascar