Next Article in Journal
Immediate Effects of Hurricanes on a Diverse Coral/Mangrove Ecosystem in the U.S. Virgin Islands and the Potential for Recovery
Next Article in Special Issue
Biogeography, Systematics, and Ecomorphology of Pacific Island Anoles
Previous Article in Journal
Parageobacillus thermantarcticus, an Antarctic Cell Factory: From Crop Residue Valorization by Green Chemistry to Astrobiology Studies
Previous Article in Special Issue
Phylogenomic Reconstruction of the Neotropical Poison Frogs (Dendrobatidae) and Their Conservation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Diversity
Volume 11
Issue 8
10.3390/d11080129
diversity-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Endemic Infection of Batrachochytrium dendrobatidis in Costa Rica: Implications for Amphibian Conservation at Regional and Species Level

by
Héctor Zumbado-Ulate
1,*,
Kiersten N. Nelson
1,
Adrián García-Rodríguez
2,3,
Gerardo Chaves
2,
Erick Arias
2,3,
Federico Bolaños
2,
Steven M. Whitfield
4 and
Catherine L. Searle
1
1
Department of Biological Sciences, Purdue University, West Lafayette, IN 47907, USA
2
Escuela de Biología, Universidad de Costa Rica, San Pedro, 11501–2060 San José, Costa Rica
3
Departamento de Zoología, Instituto de Biología, Universidad Nacional Autónoma de México, 04510 Ciudad de México, México
4
Conservation and Research Department, Zoo Miami, Miami, FL 33177, USA
*
Author to whom correspondence should be addressed.
Diversity 2019, 11(8), 129; https://doi.org/10.3390/d11080129
Submission received: 28 June 2019 / Revised: 6 August 2019 / Accepted: 7 August 2019 / Published: 9 August 2019
(This article belongs to the Special Issue Systematics and Conservation of Neotropical Amphibians and Reptiles)
Browse Figures
Versions Notes

Abstract

:
Batrachochytrium dendrobatidis (Bd) has been associated with the severe declines and extinctions of amphibians in Costa Rica that primarily occurred during the 1980s and 1990s. However, the current impact of Bd infection on amphibian species in Costa Rica is unknown. We aimed to update the list of amphibian species in Costa Rica and evaluate the prevalence and infection intensity of Bd infection across the country to aid in the development of effective conservation strategies for amphibians. We reviewed taxonomic lists and included new species descriptions and records for a total of 215 amphibian species in Costa Rica. We also sampled for Bd at nine localities from 2015–2018 and combined these data with additional Bd occurrence data from multiple studies conducted in amphibian communities across Costa Rica from 2005–2018. With this combined dataset, we found that Bd was common (overall infection rate of 23%) across regions and elevations, but infection intensity was below theoretical thresholds associated with mortality. Bd was also more prevalent in Caribbean lowlands and in terrestrial amphibians with an aquatic larval stage; meanwhile, infection load was the highest in direct-developing species (forest and stream-dwellers). Our findings can be used to prioritize regions and taxonomic groups for conservation strategies.

1. Introduction

Anthropogenic threats including habitat destruction, pollution, climate change, introduction of invasive species, and pathogens are causing a rapid and severe decline in global biodiversity [1]. Scientific consensus states that we are in the midst of a sixth mass extinction event [2,3]. Within vertebrates, amphibians are the most endangered taxonomic class with approximately 41% of described species classified as “globally threatened” [4,5]. The majority of the amphibian declines have occurred in the tropics of Australia, Central America, and South America [6,7], and have been observed even in seemingly pristine and protected environments [8,9]. However, information is still lacking regarding which species are suffering the greatest declines and which abiotic and biotic factors are contributing the most [10]. Identifying threatened species and factors contributing to global amphibian declines is vital for effective conservation and management efforts [11,12].
Costa Rica, with an area of only 51,100 km2, is home to a great diversity of amphibians [13]. More than 200 of the approximately 8000 described amphibian species are present in Costa Rica [14], and new species continue to be described. The vast species richness confined to a relatively small area is due to complex biogeographic events and climatic conditions throughout the country, and a long history of work has been done by in-country taxonomic specialists [13,15]. Costa Rica is also an example of a country where numerous amphibian population declines have occurred in response to multiple environmental threats [16], highlighted by the enigmatic disappearance of the golden toad (Incilius periglenes) [17]. However, several species that catastrophically declined in the last thirty years, such as the harlequin frog, the Golfito Robber frog, and the Holdridge’s toad, have been recently rediscovered in viable populations [18,19,20]. These findings suggest that highly susceptible species can recover from or at least persist when faced with deadly threats [21]. Thus, Costa Rica is an excellent location to study not only how amphibian communities are affected by environmental threats but also their resistance and resilience from declines [22].
One widespread cause of amphibian declines is the introduction of the pathogen Batrachochytrium dendrobatidis (Bd) [23]. This fungus causes chytridiomycosis [24], a potentially deadly skin disease that has contributed to the decline of at least 500 amphibian species globally [10]. In Central America, amphibian declines peaked during the 1980s and 1990s and have been linked to the introduction of Bd, which caused deadly outbreaks of chytridiomycosis (i.e., epizootics) [16,25,26]. It has been suggested that Bd-driven epizootic declines mostly affected species in highland lotic environments because moisture and temperature in these sites matches the optimal conditions for Bd growth in the lab [27,28]. However, it is also well known that some amphibian species suffered unexpected and unexplained declines in lowland environments (<700 m above sea level) during the 1980s and 1990s, likely due to chytridiomycosis [16,19,29,30,31]. After the declines, the evolution of resistance and tolerance mechanisms in amphibian communities [32], and/or the evolution of less-pathogenic strains of Bd [33], might have allowed susceptible amphibians to persist with endemic Bd infection (i.e., enzootics) [22,34,35,36]. However, susceptible species are still at a high risk of extinction under endemic infection if conditions shift in favor of the pathogen. For example, the introduction of an invasive species that is also a competent reservoir might amplify infection in the environment to epizootic levels [37,38,39]. Thus, examining the life history traits and conditions that may favor outbreaks of Bd is the key to understanding the underlying mechanisms behind why some infected species declined more severely than others and which species are most vulnerable to future outbreaks [40,41].
In this study, we present an updated list of all the amphibian species of Costa Rica, quantifying species diversity in each herpetological province and describing their conservation status. We also identified the effect of geography (herpetological province and altitudinal belt) and life-history traits associated with foraging and reproduction on current infection with Bd. For this, we sampled for Bd at nine tropical localities across Costa Rica from 2015–2018. In addition, we built a robust dataset by adding records from studies that detected Bd across Costa Rica from 2005–2018 in multi-species amphibian assemblages. We hypothesized that Bd is widespread across herpetological provinces and altitudinal belts in Costa Rica and would exhibit an infection intensity below theoretical thresholds associated with mass mortalities [42]. To compare across life-history traits, we developed an index that combines foraging habitat, reproductive habitat, and type of development. We hypothesized that Bd infection would vary across habitats, with the highest prevalence and infection intensity found in species with the greatest use of cool and humid environments [43,44]. The knowledge from this work will aid policy-makers in identifying the most threatened regions and taxonomic clades to develop better conservation strategies in Costa Rica [22,45,46].

2. Materials and Methods

2.1. Species Assessment

We updated the last official list of amphibian species in Costa Rica published in 2011 [47] by consulting the Herpetological Database (“Herp Database”) of the Museo de Zoología at Universidad de Costa Rica (http://museo.biologia.ucr.ac.cr/) and taxonomists’ lists [48,49]. In addition, we georeferenced the distribution of all amphibian species within the five Costa Rican herpetological provinces (see Section 2.3.1). For this, we extracted all collection points for each species from the “Herp Database” (Datum WGS1984) and mapped them using a shapefile of the Costa Rican herpetological provinces and QGIS software 3.8.1 (QGIS Development Team, http://qgis.osgeo.org).
For every species, we showed their status in Costa Rica [50] according to the International Union of Conservation of Nature (IUCN) [51] as follows: NA = “not applicable,” DD = “data deficient,” LC = “least concerned,” NT = “near threatened,” VU = “vulnerable,” EN = “endangered,” CR = “critically endangered,” and EX = “extinct in the wild” (for additional details see http://www.iucnredlist.org/). We also included environmental vulnerability scores (EVS) [52], a regional vulnerability index that classifies amphibians and reptiles into four levels of risk: “no immediate risk” (EVS < 3), “low vulnerability” (EVS of 3–9), “medium vulnerability” (EVS of 10–13), and “high vulnerability” (EVS of 14–17). A high EVS indicates species that are restricted in distribution, occur in a single life zone, and have a highly derived reproductive mode. The EVS for Costa Rican amphibians reported here were extracted from Sasa et al. [49]. Finally, we compiled a list of all the species that have been screened for Bd and the methods used for detection: histology or polymerase chain reaction (PCR).

2.2. Field Dataset

To add to existing datasets of amphibian distribution and Bd infection, we surveyed nine amphibian assemblages across Costa Rica in both versants (Caribbean and Pacific) and at elevations ranging from sea level to 1385 m (Figure S1). All surveys were conducted during the months of June and July between 2016–2018, except in the locality of Alto Lari, which was sampled in March 2015. At each site, we conducted visual and acoustic encounter surveys searching for amphibians in streams, ponds/puddles, and forest (leaflitter and canopy), and then caught individuals to screen them for Bd (see below). In total, we screened for Bd from 267 amphibians from 33 species (see Tables S1 and S2, Figure S2). Four of those species were classified in threatened categories: Oophaga granulifera (VU), Ptychohyla legleri (EN), Craugastor ranoides (CR), and C. taurus (CR).
All observed amphibians were collected using nitrile gloves and temporally placed individually in clean, unused plastic bags. Each individual was inspected for visible signs of chytridiomycosis, such as hyperplasia, hyperkeratosis, abnormal shedding, depigmentation, and lethargic behavior [24,53]. We swabbed (using MW-113 swabs) each individual’s skin to detect Bd as follows: five strokes on one hand, five strokes on the ventral patch, five strokes on one foot, and five strokes along the inner thigh. The swabs were placed into 1.5 mL screw-cap tubes and stored dry at −20 °C until fungal DNA extraction. Once swabbed, all animals were released back to the site of their collection. During this study we followed field protocols [54,55] approved by the National System of Conservation Areas of Costa Rica (SINAC), the Comisión Nacional para la Gestión de la Biodiversidad (CONAGEBIO), and animal care protocols from the Purdue Animal Care and Use Committee (PACUC 1604001392), ensuring that all animals were being cared for in accordance with standard protocols and treated in an ethical manner (research permits 001-2012-SINAC, R-019-2016-OT-CONAGEBIO, R-023-2016-OT-CONAGEBIO, R-057-2016-OT-CONAGEBIO, R-060-2016-OT-CONAGEBIO).
We conducted diagnostic quantitative polymerase chain reaction (qPCR) on each swab to quantify Bd infection load following standard protocols [56], with the following modifications: (1) the fungal DNA was extracted using 60 µL of PrepMan Ultra, and (2) an internal positive control (IPC) was used to detect inhibitors [57]. Fungal DNA was diluted 1:10 in 0.25× TE buffer and run in singlicate [58] using a Step One Plus (Applied Biosystems, Woburn, MA, USA). Negative controls (DNase/RNase-free water) were run in triplicate on every 96-well qPCR plate. We classified samples as positive when both dyes (Bd probe and IPC) amplified in each well. Samples absent of IPC amplification were considered inhibited. In order to eliminate inhibitors, we diluted 5 µL of a new dilution in 0.25×TE buffer in a proportion of 1:100. Ten samples were classified as inhibited and then determined to be negative after dilution. Quantification curves for genomic equivalents were constructed using 1000, 100, 10, and 1 zoospore quantification standards derived from a gBlock® Gene fragment (Integrated DNA Technologies, Coralville, IA, USA). In order to calculate the zoospore genomic equivalents in the original sample, we multiplied the qPCR score by the dilution factor of 120 (dilution factor = 60 × 20 × 1/10). We estimated prevalence with 95% binomial confidence intervals (CIs) by locality.

2.3. Combined Dataset

We generated a dataset from multiple studies that screened for Bd in multiple amphibian assemblages in Costa Rica after the year 2000 using conventional PCR and qPCR methods [59,60,61,62,63,64,65] (Figure 1, Table 1) (including the 267 individuals from the 33 species we tested in the “field dataset” (see Section 2.2 and supporting data). In total, this “combined dataset” consisted of 1750 individual records from 79 species and 20 localities at elevations ranging from sea level to 2000 m. We identified the year 2000 as the starting of post-decline because most epizootic outbreaks of Bd occurred during the 1980s and early 1990s [26,30,66]. We also assumed that Bd expanded its range across Costa Rica by 2000 due to the rapid rate of spreading that this pathogen exhibits in tropical locations [67,68]. Although Bd was detected in 405 swabs in this dataset, quantification through qPCR was conducted only in 351 Bd-positive swabs (from the “field dataset” and three of the seven reviewed studies [61,62,63]). We did not consider studies that used histology as a method of detection because most of these studies evaluated samples that were taken before 2000. We also excluded records of individuals that were identified only at the genus level, e.g., Craugastor spp. [64] and Agalychnis spp. [61], and cases where only one species was screened for Bd (e.g., Atelopus varius in the locality of Uvita [65]). Finally, we classified all sampled amphibians according to herpetological province, altitudinal belt, and life history traits (foraging habitat, reproductive habitat, and type of development).

2.3.1. Herpetological Provinces

We classified all surveyed assemblages within the five herpetological provinces proposed by Savage [13] and modified by Sasa and colleagues [49].
Caribbean Lowlands: This faunal area represents 30% of Costa Rica and includes the lowlands of the Caribbean versant and the northernmost region of the country, predominantly consisting of lowland wet forest. Sampling for Bd through PCR has been conducted in the localities of La Virgen de Sarapiquí, Puerto Viejo de Sarapiquí, Tres Equis de Turrialba, Guayacán de Siquirres, Kekoldi, and the remote Alto Lari. Alto Lari was surveyed as part of a recent expedition following an enigmatic path that connects the Caribbean Lowlands with the highlands of Cordillera de Talamanca and is known as “the Gabb’s route” [69].
Pacific Northwest: This herpetological province includes the lowlands of the Pacific Northwest and extends into the western side of the Central Valley, in the Meseta Central Occidental (Central Valley) up to the base of Cerros de Ochomogo. The Pacific Northwest consists of predominantly Lowland Dry Forest vegetation and constitutes 24% of Costa Rica’s area. This province contains a distinctive dry season that lasts five to six months. Within the Pacific Northwest, sampling for Bd has been conducted in the tropical dry forest at Guanacaste National Park (Santa Rosa and Santa Elena Peninsula stations).
Pacific Southwest: Encompassing the lowlands of the Pacific central and south, the herpetological provinces consist primarily of lowland wet forest and lowland moist forest and accounts for 15% of the country’s area. This herpetological province is biogeographically related to the Caribbean Lowlands and species have more recently been differentiated between these herpetological provinces due to isolation caused by the uplifting of the Cordillera de Talamanca. Within the Pacific Southwest, sampling for Bd through PCR has been conducted in the localities of Punta Banco–Burica, Golfito, Rincón de Osa, Piro, Corcovado, and Uvita.
Montane Slopes and Cordillera Central: This area represents 23% of Costa Rica and occurs along all of Costa Rica’s mountain ranges from 500–2100 m elevation in Cordillera de Guanacaste, 500–3400 in Cordillera Central, and 500–1600 in Cordillera de Talamanca (Lower Talamanca). The Montane Slopes and Cordillera Central province includes regions that receive the highest annual precipitation in the country. Sampling for Bd through PCR has been conducted in the localities of Monteverde, San Vito de Coto Brus, Las Tablas, Tinamastes de Pérez Zeledón, San Rafael de Heredia, and Santo Domingo de Heredia.
Cordillera de Talamanca: Found at elevations above 1600 m (Upper Talamanca). This is the smallest faunal area (8% of Costa Rica) and consists primarily of montane rainforest and subalpine pluvial paramo. This faunal area is the least explored herpetological province of Costa Rica and there is no published data for Bd detection through PCR or qPCR in this faunal province.

2.3.2. Altitudinal Belt

We classified species according to the Holdridge’s life zone system [70], which divides Costa Rica into five altitudinal belts: lowland, premontane, lower montane, montane, and subalpine. Due to the elevational limits for altitudinal belts being slightly different among regions in Costa Rica, we established the limits of each belt as follows: lowland (0–700 m), premontane (700–1500 m), lower montane (1500–2700 m), montane (2700–3500 m), and subalpine (>3500 m).

2.3.3. Foraging-Reproduction Habitat Index

To compare Bd infection dynamics across taxonomic groups, we developed a foraging–reproduction habitat index (FRHI). The FRHI was created to classify species with a system of three letters that represented life history traits associated with foraging and reproduction (Table 2). First, we classified species according to their development into indirect- (I) or direct-developing amphibians (D). Second, we classified species according to their foraging habitat into terrestrial (T), arboreal (A), pond/puddle-dwellers (P), stream-dwellers (R), and phytotelma (F). Finally, species were classified according to their reproductive habitat into terrestrial (T), arboreal (A), pond/puddle-breeders (P), stream-breeders (R), and phytotelma (F).

2.4. Statistical Analysis

We reduced our “combined dataset” to 1741 samples from 74 species and 20 localities because there was insufficient information to accurately classify nine records of the species Diasporus tigrillo, D. vocator, Hyloscirtus palmeri, Triprion spinosus, and Cruziohyla calcarifer in the FRHI. For our analyses, we pooled all species together instead of using species as a predictor because the samples sizes per species were highly variable (from 1–177), which could produce significant models that may be an artifact of opportunistic sampling instead of a real pattern. Instead, we used the FRHI, which is highly correlated with taxonomic group. We were unable to include time of sampling as a predictor in our analyses because these data were missing in several of the amphibian assemblages sampled. All our analyses were performed with the R package “stats” [71].
We analyzed Bd prevalence with fixed-effects generalized linear models (GLMs) using infected status as a binomial response variable (uninfected or infected) and herpetological province, altitudinal belt, and the FRHI as predictors. Ranking of the candidate GLMs followed the Akaike’s information criterion (AIC) where the model with the lowest AIC was considered the most robust [72]. To analyze Bd infection intensity (estimated as the number of genomic equivalents), we analyzed the 351 Bd-positive swabs where Bd infection intensity was quantified through qPCR (see Section 2.3). We used linear models (LMs) to compare Bd infection intensity (response variable) across herpetological provinces, altitudinal belts, and FRHI (predictor variables). We log-transformed the Bd infection intensity to reduce skewness. Statistical significance of models was tested with ANOVA. For both, GLMs and LMs, we conducted post hoc, pairwise comparisons (Tukey’s honestly significant difference; HSD-test) to confirm where the differences occurred between significant predictors. We were unable to run mix-effects models or fixed-effects interaction models because some combinations of predictors presented missing or low values, causing convergence difficulties.

3. Results

3.1. Species Assessment

The previous list of Costa Rican amphibians from 2011 included 196 species; however, we excluded the species Pristimantis educatoris [73] due to taxonomic uncertainty [74] for a total of 195 species. Our new list of amphibians in Costa Rica included a total of 215 species grouped in three orders, 16 families, and 48 genera (Table S3). This represented and addition of 20 species (10 anurans, 9 salamanders, and 1 caecilian; Table 3). The order Anura (frogs and toads) is the most diverse in Costa Rica, being 72% of the total species (13 families and 41 genera). Salamanders (order Caudata) are represented by only one family (Plethodontidae), three genera, and 53 species. Finally, caecilians (order Gymnophiona) are represent by two families, four genera, and eight species. A total of 63 species are endemic to Costa Rica (36 salamanders, 24 anurans, and 3 caecilians). We also included five anuran species that are not native to Costa Rica (Eleutherodactylus coqui, E. johnstoni, E. planirostris, Lithobates catesbeianus, and Osteopilus septentrionalis).
Regionally (Figure 2a, Table S3), the Cordillera de Talamanca is the most diverse province in terms of species per area (88 species, 2.2 species/100 km2). It contains 23 species of amphibians that only occurred within this herpetological province (e.g., Diasporus ventrimaculatus, Nototriton costaricense). The Montane Slopes and Cordillera Central is the most diverse herpetological province (188 species, 1.3 species/100 km2), with 27 species that are exclusively found within this herpetological province (e.g., Cochranella euknemos, Nototriton picadoi). The Caribbean Lowlands (101 species, 0.7 species/100 km2) includes 20 species that are only found within this herpetological province (e.g., Cruziohyla calcarifer, Caecilia volcani). The Pacific Southwest (71 species, 0.9 species/100 km2), has five species that only occur within this province (e.g., Craugastor taurus, Oophaga granulifera). Finally, the Pacific Northwest (66 species, 0.5 species/100 km2) includes only two species that are found exclusively within this province (Rhinophrynus dorsalis and Eleutherodactylus johnstonei). A total of 20 species occur in all five herpetological provinces (e.g., Craugastor fitzingeri, Diasporus diastema, Dendropsophus ebraccatus, Hyalinobatrachium fleischmanni, Lithobates warszewitschii, Rhinella horribilis, Smilisca sordida).
Overall, 200 species have been classified into IUCN categories at the country level and 15 species need future assessment. A total of 155 species do not fulfill the criteria to be considered in the threatened categories, including 136 species listed as LC, 18 as DD, and one as NA (a category for taxa that occur in the region but have been excluded from the regional Red List for a specific reason). Within threatened categories, two species are listed as EX, 24 as CR, ten as EN, seven as VU, and two as NT. Regionally, lowlands exhibited the lowest percentage (0–2%) of DD species (Figure 2b). Similarly, about 75–80% of species occurring in lowlands are listed as LC (Figure 2c). In highlands, 6–10% of species are categorized as DD (Figure 2b) and 26% of species in Cordillera de Talamanca are classified within threatened categories (Figure 2d). According to EVS, a total of 81 species were classified as “no immediate risk,” four species at “low vulnerability,” 50 species at “medium vulnerability,” and 48 species at “high vulnerability” (Table S3).
In our review, we found a total of 105 amphibian species (49%) that have been screened for Bd in Costa Rica (103 anurans and only 2 species of salamanders) (Table S4). In the field, the most common method used to detect Bd was qPCR, especially after 2005. Conventional PCR was used only in one study in the Caribbean Lowlands [64]. Histology and qPCR have also been used in retrospective studies on preserved specimens from declined and extinct species.

3.2. Endemic Dynamics

Overall, Bd prevalence in Costa Rica was estimated to be 0.23 (60% of species tested positive for Bd) (Table S5). The most robust GLM found both herpetological province and the FRHI as significant predictors of Bd prevalence (AIC = 1700, p < 0.001, Table S6). Among herpetological provinces (Figure 3), the highest percentage of infected individuals was found in the Caribbean Lowlands (34%) and the lowest in the Pacific Northwest (4%). The Mountain Slopes, Cordillera Central, and Pacific Southwest had a similar percentage of infected individuals (≈23%). Furthermore, Bd was proportionally more prevalent in amphibians with terrestrial foraging and larval stage in phytotelma (ITF), pond-breeding treefrogs (IAP), and direct-developing species that breed in the forest (leaf-litter frogs DTT, rain frogs, DAT) (Figure 4a).
The species Craugastor taurus (the Golfito robber frog) was the species that had the highest average infection load (average Bd load of 11632.4 versus 571.6 genomic equivalents or 2.51 versus 1.18 after log transformation) (Table 4). We found an effect of the FRHI on infection load (F8342 = 7.91, p < 0.01, Table S6). Direct-developing frogs with terrestrial reproduction (robber frogs and leaf-litter frogs; DTR and DTT respectively) had the highest Bd loads (Figure 4b, Table 4).

4. Discussion

4.1. Species Assessment

We presented the first updated list of Costa Rican amphibians since 2011 [47]. Compared to the last list, we added 10 anurans, 9 salamanders, and 1 caecilian (Table 3) for a total of 215 species (Table S3). As is common throughout the world, anurans exhibited the highest amphibian species richness in Costa Rica, with 72% of listed species. However, the richness of salamander species is also high (25%). In Costa Rica, the diversity and endemism of amphibians (especially salamanders) increase with elevation and complex mountain topography [13]. Proportionally, in terms of number of species per unit area (km2), the richest herpetological province was the Cordillera de Talamanca. In this herpetological province, the number of species continued to increase and most of the newly described species in our report came from this remote and almost inaccessible province [69]. The Montane Slopes and Cordillera Central presented the highest number of species (158 species). Within this province, numerous mountain ranges provide multiple microhabitats for niche differentiation and further speciation [13,49]. In lowlands, the highest number of amphibian species occurred in the Caribbean Lowlands with 101 species. However, the Pacific Southwest presented more species per unit area. The Pacific Northwest only had 66 species, which was also the lowest number of species per unit area. This pattern may be attributed to the warm and dry conditions that occur in most part of this herpetological province [13,49].
According IUCN, the species Craugastor escoces and Incilius periglenes are classified as EX in Costa Rica. However, C. escoces was recently rediscovered [89]. Similarly, several species that remained undetected after the 1980s and 1990s such as Incilius holdridgei [20], Craugastor taurus [19], and Atelopus varius [18] have been rediscovered in peripheral populations during the last few years. However, the number of extinct species could be higher because multiple threatened species still remain undetected in the field (e.g., Craugastor andi, Incilius fastidiosus, Atelopus senex). We recommend expedition surveys to find populations of declined and data-deficient species [90] and captive-breeding for species where ex situ reproduction has been successful (e.g., harlequin frogs.) [91]. Although we acknowledge that there are limited funds available for these types of conservation efforts, knowledge from these sites is essential to be able to identify conditions that favor the persistence of threatened species and identify species that should be targeted for future conservation efforts [92].
Lowlands of Costa Rica exhibited the lowest proportion of DD species (0–2%; Figure 2b) and the highest proportion of LC species (75–80%; Figure 2c). On the other hand, highlands exhibited the highest percentage of DD species (6–10%; Figure 2b). Similarly, Cordillera de Talamanca had the highest percentage of threatened species (26%; Figure 2d). Based on these findings, we strongly recommend increasing the sampling effort in the montane and subalpine altitudinal belts (>2800 m) that exclusively occur in Cordillera de Talamanca and Montane Slopes and Cordillera Central. These herpetological provinces present the highest rate of endemism (especially for salamanders) and contain several of the recently described species [85,86]. Conducting expeditions and long-term studies in highlands will aid in monitoring threatened species and reducing information gaps, allowing for more accurate assessments of amphibian species.
To better evaluate the vulnerability of amphibian species, we utilized EVS (Table S3). This index relies on ecological information for categorizing threat levels, which makes application easy for most species in a specific region [52]. Unlike the IUCN Red List of Threatened Species, previous evaluations of threats are not considered by this index [49]. For that reason, species that are classified as LC by IUCN can be classified as highly vulnerable in this index (e.g., Duellmanohyla rufioculis). We categorized 48 species (24 salamanders, 21 anurans, and 3 caecilians) in “high vulnerability” (e.g., Atelopus senex, Craugastor andi, Bolitoglossa pesrubra, Nototriton guanacaste, and Oscaecilia osae). These species exhibited the highest EVS values because their habitats are restricted and because they exhibit complex reproductive modes. Quantifying environmental threats and combining information from both indexes will help policy-makers to prioritize conservation actions for threatened species.
In Central America, habitat destruction is the most important threat impacting amphibian populations [49,52]. Although approximately 30% of Costa Rica remains forested and protected, rapid urbanization, extensive agriculture, excessive pesticide use, illegal traffic, and inappropriate waste management negatively affect numerous amphibian populations. However, even in seemingly pristine locations, amphibian declines have occurred [93]. Additionally, climate change has been associated with the decline of several amphibian species in Costa Rica by affecting their reproduction and likely increasing susceptibility to pathogens [28]. Although it has not been found in Central America, we recommend screening for the recently emerged fungus Batrachochytrium salamandrivorans [94], which causes chytridiomycosis in salamanders. Conveniently, swabbing methods and qPCR allow for accurate detection of both fungi species in the same assay, which may facilitate rapid population assessments. For a fully detailed review of the environmental threats for amphibian communities in Costa Rica, we recommend the work of Sasa et al. [49].

4.2. Post-Decline Dynamics

In this study, we found strong evidence that Bd is widespread in Costa Rica. Our results also suggest that post-decline Bd exhibited enzootic dynamics, characterized by high prevalence of infection across regions and pathogenic loads below thresholds associated with mass mortalities [39,42]. We found Bd in all the herpetological provinces and altitudinal belts surveyed (Figure 3 and Figure 4, Tables S5 and S6), for a total infection rate of 23%. We also found that Bd was more prevalent in terrestrial amphibians with an aquatic larval stage and direct-developing frogs exhibited the highest pathogenic loads.
We found the lowest infection rate (<5%) in the Pacific Northwest and the highest (33%) in the Caribbean Lowlands (Figure 3). However, these values may be an effect of sampling periods. A study conducted at La Selva Biological Station in the Caribbean Lowlands [95] reported infection rates varying from <5% during the warmest months (May to early November) to 35% in the coolest months (mid-November to January) in three common amphibian species. In addition, the gradual population declines observed at La Selva over several decades [31] and opportunistic observations of small-scale mortality events during cold periods [96] suggest that Bd may be causing mortality in amphibians long after its initial invasion. Similar mortality events in response to seasonality could be occurring in amphibian communities in “refuges from decline” in the Pacific Northwest [62,97] and Pacific Southwest [63] of Costa Rica; however, they have not yet been documented. There, seasonal changes in precipitation and temperature caused Bd prevalence to vary from >5% in the peak of the dry season (March and April) to 80% in the coldest months (November–December). Therefore, we recommend follow-up studies at these sites to identify whether seasonal disease dynamics are causing mortality events in regions that have been considered unsuitable for Bd [98].
Bd was found across all altitudinal belts across Costa Rica. Similar results of high Bd prevalence across all elevations has been found in Panama [99,100,101]. These findings suggest that current environmental conditions are suitable for Bd at most elevations in Central America [63]. It is also plausible that Bd-driven declines during the 1980s and 1990s were not exclusively restricted to highlands [30] but were relatively undetected at lower elevations. Another hypothesis is that species with high susceptibility historically occupied high elevations sites, but severely declined or went extinct after Bd was introduced, leaving only species with mid-to-low susceptibility across elevations [102]. On the other hand, the absence of samples from montane and subalpine belts and uncontrolled variables (e.g., changes in species composition, climatic disturbances) could have reduced the statistical power to determine changes in Bd prevalence across altitudinal belts. We suggest that future studies increase sampling at high elevations (>2700 m) to better understand the local spatial dynamics of Bd across elevations.
Our results showed that infection with Bd was common in amphibians across all life-history traits evaluated in the FRHI (Table 2). However, Bd was significantly more prevalent in terrestrial amphibians with a larval stage (Figure 4a), especially those that complete metamorphosis in phytotelma (ITF). All the species within the ITF category belonged to the family Dendrobatidae, which have previously been shown to easily acquire Bd infection [62,95]. The high susceptibly of the Dendrobatidae family is likely due to their preferred habitat (e.g., water-filled bromeliads for many species), as it offers suitable conditions for Bd infection [103]. In addition, dendrobatid adults forage in the tropical forest floor and stream-associated low vegetation, which are environments that can sustain Bd [104]. Regarding infection intensity (Figure 4b), the FRHI showed similar results to studies that have used the aquatic index [43]. We found that direct-developing species with terrestrial reproduction had significantly higher infection load than other species with different life-history traits (Figure 4b, Table 4). This life-history trait is exhibited by leaf-litter frogs and all the species within the Craugastor punctariolus clade (robber frogs), which is one of the clades most affected by chytridiomycosis in Central America [29,105]. Robber frogs spend a majority of their life cycle along fast-flowing streams [106], an aquatic environment that seems highly suitable for Bd in Central America [107]. In addition, these frogs appear to be highly susceptible to Bd-driven mass mortalities outside warm and dry ecosystems [19,30,108]. Our results suggest the FRHI is particularly useful for identifying taxonomic units that are more susceptible to Bd [92].

5. Conclusions

Our results demonstrated that the number of identified amphibian species in Costa Rica is still growing, and there may be potential future additions (e.g., Bolitoglossa anthracina [109] and B. indio [110]). A continuous assessment of species and regions is needed to identify continuing threats to amphibian biodiversity. We found that Bd was widespread across species, herpetological provinces, and altitudinal belts in samples collected since 2000. Conducting more studies in remote regions, such as Cordillera de Talamanca, may help to better describe spatial dynamics of both amphibian hosts and Bd. In addition, future studies should test whether seasonal disease dynamics are causing mortality events in regions that are considered unsuitable for Bd. Under potential scenarios of climate change, environmental conditions may shift to ideal ranges for Bd infection [28] and seasonal regions that sustain critically endangered species (e.g., tropical dry forest) may experience future outbreaks of chytridiomycosis [111]. We also recommend continuous surveillance of invasive species, which might amplify Bd in the environment, causing future epizootics [37]. This vital information will aid in the development of more effective conservation strategies for amphibians across a broader range of habitats [46,112,113,114].

Supplementary Materials

The following are available online at https://www.mdpi.com/1424-2818/11/8/129/s1. Figure S1: Map of Costa Rica showing elevational gradient and nine study sites surveyed for Batrachochytrium dendrobatidis, Figure S2: Prevalence and intensity of infection of Batrachochytrium dendrobatidis in amphibian assemblages from the nine surveyed sites in Costa Rica, Table S1: List of species where Batrachochytrium dendrobatidis was surveyed in Costa Rica for our “field dataset”, Table S2: Prevalence and infection intensity of Batrachochytrium dendrobatidis at nine sites in Costa Rica, Table S3: Distribution of amphibians in Costa Rica according to herpetological province and elevation, Table S4: List of species that have been screened for Batrachochytrium dendrobatidis in Costa Rica, Table S5: List of species where Batrachochytrium dendrobatidis was surveyed in Costa Rica in our combined dataset, Table S6: Candidacy generalized linear models (GLMs) and linear models (LMs) used to determine the best predictors of prevalence and infection intensity of Batrachochytrium dendrobatidis in amphibian assemblages from Costa Rica.

Author Contributions

Conceptualization, H.Z.; Data curation, H.Z., K.N., A.G., G.C., E.A., F.B., S.W., and C.S.; Formal analysis, H.Z., K.N., A.G., and G.C.; Methodology, H.Z., K.N., A.G., G.C., E.A., F.B., S.W., and C.S.; Project administration, C.S.; Supervision, H.Z. and C.S.; Writing—Original draft, H.Z. and K.N.; Writing—Review and editing, H.Z., K.N., A.G., G.C., E.A., F.B., S.W., and C.S.

Funding

H.Z. was funded by (1) the Department of Biological Sciences at Purdue University through the Lindsay Fellowship and the Andrews Environmental Travel Grant, (2) the Osa Conservation through the Álvaro Ugalde Scholarship, and (3) IdeaWild. K.N. was funded by the Department of Biological Sciences at Purdue University through the McAtee Scholarship and a Cable-Silkman award. A.G. is currently supported by a postdoctoral fellowship from Dirección General de Asuntos del Personal Académico (DGAPA) at Instituto de Biología, Universidad Nacional Autónoma de México. The field trip following the Gabb’s route was funded by the National Geographic Society (grant number W-346-14).

Acknowledgments

We thank B. Araya, F. Ureña, and, P. Piza for field assistance, and D. Vasquez for statistical counseling. We also acknowledge the Costa Rican Ministry of Environment and Energy (MINAE)) for providing the corresponding scientific collecting permits for this work (research permits 001-2012-SINAC, R-019-2016-OT-CONAGEBIO, R-023-2016-OT-CONAGEBIO, R-057-2016-OT-CONAGEBIO, R-060-2016-OT-CONAGEBIO).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Novacek, M.J.; Cleland, E.E. The current biodiversity extinction event: Scenarios for mitigation and recovery. Proc. Natl. Acad. Sci. USA 2001, 98, 5466–5470. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Barnosky, A.D.; Matzke, N.; Tomiya, S.; Wogan, G.O.U.; Swartz, B.; Quental, T.B.; Marshall, C.; McGuire, J.L.; Lindsey, E.L.; Maguire, K.C.; et al. Has the Earth’s sixth mass extinction already arrived? Nature 2011, 471, 51–57. [Google Scholar] [CrossRef] [PubMed]
  3. Wake, D.B.; Vredenburg, V.T. Are we in the midst of the sixth mass extinction? A view from the world of amphibians. Proc. Natl. Acad. Sci. USA 2008, 105, 11466–11473. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Monastersky, R. Life—A status report. Nature 2014, 516, 158–161. [Google Scholar] [CrossRef] [PubMed]
  5. Stuart, S.N.; Chanson, J.; Cox, N.A.; Young, B.E.; Rodrigues, A.S.L.; Fischman, D.; Waller, R. Status and trends of amphibian declines and extinctions worldwide. Science 2004, 306, 1783–1786. [Google Scholar] [CrossRef]
  6. Catenazzi, A. State of the World’s Amphibians. Annu. Rev. Environ. Resour. 2015, 40, 91–119. [Google Scholar] [CrossRef]
  7. Daszak, P.; Berger, L.; Cunningham, A.A.; Hyatt, A.D.; Green, D.E.; Speare, R. Emerging infectious diseases and amphibian population declines. Emerg. Infect. Dis. 1999, 5, 735–748. [Google Scholar] [CrossRef] [PubMed]
  8. Collins, J.P. Amphibian decline and extinction: What we know and what we need to learn. Dis. Aquat. Organ. 2010, 92, 93–99. [Google Scholar] [CrossRef] [PubMed]
  9. La Marca, E.; Lips, K.R.; Lotters, S.; Puschendorf, R.; Ibanez, R.; Rueda-Almonacid, J.V.; Schulte, R.; Marty, C.; Castro, F.; Manzanilla-Puppo, J.; et al. Catastrophic population declines and extinctions in Neotropical harlequin frogs (Bufonidae: Atelopus). Biotropica 2005, 37, 190–201. [Google Scholar] [CrossRef]
  10. Scheele, B.C.; Pasmans, F.; Skerratt, L.F.; Berger, L.; Martel, A.; Beukema, W.; Acevedo, A.A.; Burrowes, P.A.; Carvalho, T.; Catenazzi, A.; et al. Amphibian fungal panzootic causes catastrophic and ongoing loss of biodiversity. Science 2019, 363, 1459–1463. [Google Scholar] [CrossRef]
  11. Gerber, B.D.; Converse, S.J.; Muths, E.; Crockett, H.J.; Mosher, B.A.; Bailey, L.L. Identifying species conservation strategies to reduce disease-associated declines. Conserv. Lett. 2018, 11, e12393. [Google Scholar] [CrossRef]
  12. Meredith, H.; Van Buren, C.; Antwis, R.E. Making amphibian conservation more effective. Conserv. Evid. 2016, 13, 1–6. [Google Scholar]
  13. Savage, J.M. The Amphibians and Reptiles of Costa Rica: A Herpetofauna between Two Continents, between Two Deas; University of Chicago Press: Chicago, IL, USA, 2002. [Google Scholar]
  14. Frost, D.R. Amphibian Species of the World: An Online Reference. Version 6.0. Available online: http://research.amnh.org/vz/herpetology/amphibia/ (accessed on 27 June 2019).
  15. Bagley, J.C.; Johnson, J.B. Phylogeography and biogeography of the lower Central American Neotropics: Diversification between two continents and between two seas. Biol. Rev. 2014, 89, 767–790. [Google Scholar] [CrossRef] [PubMed]
  16. Bolaños, F. Situación de los anfibios de Costa Rica. Biocenosis 2009, 22, 95–108. [Google Scholar]
  17. Pounds, J.A.; Crump, M.L. Amphibian declines and climate disturbance: The case of the golden toad and the harlequin frog. Conserv. Biol. 1994, 8, 72–85. [Google Scholar] [CrossRef]
  18. González-Maya, J.F.; Belant, J.L.; Wyatt, S.A.; Schipper, J.; Cardenal, J.; Corrales, D.; Cruz-Lizano, I.; Hoepker, A.; Escobedo-Galván, A.H.; Castañeda, F.; et al. Renewing hope: The rediscovery of Atelopus varius in Costa Rica. Amphib. Reptil. 2013, 34, 573–578. [Google Scholar] [CrossRef]
  19. Chaves, G.; Zumbado-Ulate, H.; García-Rodríguez, A.; Gómez, E.; Vredenburg, V.T.; Ryan, M.J. Rediscovery of the critically endangered streamside frog, Craugastor taurus (Craugastoridae), in Costa Rica. Trop. Conserv. Sci. 2014, 7, 628–638. [Google Scholar] [CrossRef]
  20. Abarca, J.; Chaves, G.; García-Rodríguez, A.; Vargas, R. Reconsidering extinction: Rediscovery of Incilius holdridgei (Anura: Bufonidae) in Costa Rica after 25 years. Herpetol. Rev. 2010, 41, 150. [Google Scholar]
  21. Hero, J.-M.; Williams, S.E.; Magnusson, W.E. Ecological traits of declining amphibians in upland areas of eastern Australia. J. Zool. 2005, 267, 221–232. [Google Scholar] [CrossRef]
  22. Mendelson, J.R.; Whitfield, S.M.; Sredl, M.J. A recovery engine strategy for amphibian conservation in the context of disease. Biol. Conserv. 2019, 236, 188–191. [Google Scholar] [CrossRef]
  23. Longcore, J.E.; Pessier, A.P.; Nichols, D.K. Batrachochytrium dendrobatidis gen. et sp. nov., a chytrid pathogenic to amphibians. Mycologia 1999, 91, 219–227. [Google Scholar] [CrossRef]
  24. Berger, L.; Speare, R.; Daszak, P.; Green, D.E.; Cunningham, A.A.; Goggin, C.L.; Slocombe, R.; Ragan, M.A.; Hyatt, A.D.; McDonald, K.R.; et al. Chytridiomycosis causes amphibian mortality associated with population declines in the rain forests of Australia and Central America. Proc. Natl. Acad. Sci. USA 1998, 95, 9031–9036. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Lips, K.R.; Green, D.E.; Papendick, R. Chytridiomycosis in wild frogs from southern Costa Rica. J. Herpetol. 2003, 37, 215–218. [Google Scholar] [CrossRef]
  26. Puschendorf, R.; Bolaños, F.; Chaves, G. The amphibian chytrid fungus along an altitudinal transect before the first reported declines in Costa Rica. Biol. Conserv. 2006, 132, 136–142. [Google Scholar] [CrossRef]
  27. Piotrowski, J.S.; Annis, S.L.; Longcore, J.E. Physiology of Batrachochytrium dendrobatidis, a chytrid pathogen of amphibians. Mycologia 2004, 96, 9–15. [Google Scholar] [CrossRef] [PubMed]
  28. Pounds, A.J.; Bustamante, M.R.; Coloma, L.A.; Consuegra, J.A.; Fogden, M.P.L.; Foster, P.N.; La Marca, E.; Masters, K.L.; Merino-Viteri, A.; Puschendorf, R.; et al. Widespread amphibian extinctions from epidemic disease driven by global warming. Nature 2006, 439, 161–167. [Google Scholar] [CrossRef] [PubMed]
  29. Zumbado-Ulate, H.; Bolaños, F.; Willink, B.; Soley-Guardia, F. Population status and natural history notes on the critically endangered stream dwelling frog Craugastor ranoides (Craugastoridae) in a Costa Rican tropical dry forest. Herpetol. Conserv. Biol. 2011, 6, 455–464. [Google Scholar]
  30. Puschendorf, R.; Carnaval, A.C.; VanDerWal, J.; Zumbado-Ulate, H.; Chaves, G.; Bolaños, F.; Alford, R.A. Distribution models for the amphibian chytrid Batrachochytrium dendrobatidis in Costa Rica: Proposing climatic refuges as a conservation tool. Divers. Distrib. 2009, 15, 401–408. [Google Scholar] [CrossRef]
  31. Whitfield, S.M.; Bell, K.E.; Philippi, T.; Sasa, M.; Bolaños, F.; Chaves, G.; Savage, J.M.; Donnelly, M.A. Amphibian and reptile declines over 35 years at La Selva, Costa Rica. Proc. Natl. Acad. Sci. USA 2007, 104, 8352–8356. [Google Scholar] [CrossRef] [Green Version]
  32. Christie, M.R.; Searle, C.L. Evolutionary rescue in a host–pathogen system results in coexistence not clearance. Evol. Appl. 2018, 11, 681–693. [Google Scholar] [CrossRef]
  33. Retallick, R.W.R.; Miera, V. Strain differences in the amphibian chytrid Batrachochytrium dendrobatidis and non-permanent, sub-lethal effects of infection. Dis. Aquat. Organ. 2007, 75, 201–207. [Google Scholar] [CrossRef] [PubMed]
  34. Briggs, C.J.; Vredenburg, V.T.; Knapp, R.A.; Rachowicz, L.J. Investigating the population-level effects of chytridiomycosis: An emerging infectious disease of amphibians. Ecology 2005, 86, 3149–3159. [Google Scholar] [CrossRef]
  35. Rachowicz, L.J.; Knapp, R.A.; Morgan, J.A.T.; Stice, M.J.; Vredenburg, V.T.; Parker, J.M.; Briggs, C.J. Emerging infectious disease as a proximate cause of amphibian mass mortality. Ecology 2006, 87, 1671–1683. [Google Scholar] [CrossRef]
  36. Retallick, R.W.R.; McCallum, H.; Speare, R. Endemic infection of the amphibian chytrid fungus in a frog community post-decline. PLoS Biol. 2004, 2, e351. [Google Scholar] [CrossRef] [PubMed]
  37. Searle, C.L.; Biga, L.M.; Spatafora, J.W.; Blaustein, A.R. A dilution effect in the emerging amphibian pathogen Batrachochytrium dendrobatidis. Proc. Natl. Acad. Sci. USA 2011, 108, 16322–16326. [Google Scholar] [CrossRef] [PubMed]
  38. O’Brien, V.A.; Moore, A.T.; Young, G.R.; Komar, N.; Reisen, W.K.; Brown, C.R. An enzootic vector-borne virus is amplified at epizootic levels by an invasive avian host. Proc. R. Soc. B Biol. Sci. 2011, 278, 239–246. [Google Scholar] [CrossRef] [PubMed]
  39. Briggs, C.J.; Knapp, R.A.; Vredenburg, V.T. Enzootic and epizootic dynamics of the chytrid fungal pathogen of amphibians. Proc. Natl. Acad. Sci. USA 2010, 107, 9695–9700. [Google Scholar] [CrossRef] [Green Version]
  40. Miller, C.A.; Tasse Taboue, G.C.; Ekane, M.M.P.; Robak, M.; Sesink Clee, P.R.; Richards-Zawacki, C.; Fokam, E.B.; Fuashi, N.A.; Anthony, N.M. Distribution modeling and lineage diversity of the chytrid fungus Batrachochytrium dendrobatidis (Bd) in a central African amphibian hotspot. PLoS ONE 2018, 13, e0199288. [Google Scholar] [CrossRef]
  41. Hitchman, S.M.; Mather, M.E.; Smith, J.M.; Fencl, J.S. Identifying keystone habitats with a mosaic approach can improve biodiversity conservation in disturbed ecosystems. Glob. Chang. Biol. 2018, 24, 308–321. [Google Scholar] [CrossRef]
  42. Vredenburg, V.T.; Knapp, R.A.; Tunstall, T.S.; Briggs, C.J. Dynamics of an emerging disease drive large-scale amphibian population extinctions. Proc. Natl. Acad. Sci. USA 2010, 107, 9689–9694. [Google Scholar] [CrossRef] [Green Version]
  43. Brem, F.; Lips, K. Batrachochytrium dendrobatidis infection patterns among Panamanian amphibian species, habitats and elevations during epizootic and enzootic stages. Dis. Aquat. Organ. 2008, 81, 189–202. [Google Scholar] [CrossRef] [PubMed]
  44. Kriger, K.M.; Hero, J.-M. The chytrid fungus Batrachochytrium dendrobatidis is non-randomly distributed across amphibian breeding habitats. Divers. Distrib. 2007, 13, 781–788. [Google Scholar] [CrossRef]
  45. Heard, G.W.; Scroggie, M.P.; Ramsey, D.S.L.; Clemann, N.; Hodgson, J.A.; Thomas, C.D. Can habitat management mitigate disease impacts on threatened amphibians? Conserv. Lett. 2018, 11, e12375. [Google Scholar] [CrossRef]
  46. Scheele, B.C.; Hunter, D.A.; Grogan, L.F.; Berger, L.; Kolby, J.E.; Mcfadden, M.S.; Marantelli, G.; Skerratt, L.F.; Driscoll, D.A. Interventions for reducing extinction risk in chytridiomycosis-threatened amphibians. Conserv. Biol. 2014, 28, 1195–1205. [Google Scholar] [CrossRef] [PubMed]
  47. Bolanos, F.; Savage, J.M.; Chaves, G. Anfibios y Reptiles de Costa Rica. Listas Zoológicas Actualizadas UCR. Available online: http://museo.biologia.ucr.ac.cr/Listas/Anteriores/HerpCREsp.htm (accessed on 27 June 2019).
  48. Savage, J.M.; Bolaños, F. A checklist of the amphibians and reptiles of Costa Rica: Additions and nomenclatural revisions. Zootaxa 2009, 2005, 1–23. [Google Scholar] [CrossRef]
  49. Sasa, M.; Chaves, G.; Porras, L.W. The Costa Rican herpetofauna: Conservation status and future perspectives. In Conservation of Mesoamerican Amphibians and Reptiles; Townsend, J.H., Johnson, J.D., Eds.; Eagle Mountain Press: Salt Lake City, UT, USA, 2010; pp. 510–603. [Google Scholar]
  50. Chaves, G.; Bolaños, F.; Rodríguez, J.E.; Matamoros, Y. Actualización de las Listas Rojas Nacionales de Costa Rica. Anfibios y Reptiles; Conservation Breeding Specialist Group (SSC/IUCN)/CBSG Mesoamerica): San José, Costa Rica, 2014. [Google Scholar]
  51. IUCN The IUCN Red List of Threatened Species. Version 2019-1. Available online: https://www.iucnredlist.org/en (accessed on 27 June 2019).
  52. Wilson, L.D.; McCranie, J.R. The conservation status of the herpetofauna of Honduras. Amphib. Reptile Conserv. 2004, 3, 6–33. [Google Scholar] [PubMed]
  53. Voyles, J.; Young, S.; Berger, L.; Campbell, C.; Voyles, W.F.; Dinudom, A.; Cook, D.; Webb, R.; Alford, R.A.; Skerratt, L.F.; et al. Pathogenesis of chytridiomycosis, a cause of catastrophic amphibian declines. Science 2009, 326, 582–585. [Google Scholar] [CrossRef]
  54. Kriger, K.M.; Hines, H.B.; Hyatt, A.D.; Boyle, D.G.; Hero, J.-M. Techniques for detecting chytridiomycosis in wild frogs: Comparing histology with real-time Taqman PCR. Dis. Aquat. Organ. 2006, 71, 141–148. [Google Scholar] [CrossRef]
  55. Skerratt, L.; Berger, L.; Hines, H.; McDonald, K.; Mendez, D.; Speare, R. Survey protocol for detecting chytridiomycosis in all Australian frog populations. Dis. Aquat. Organ. 2008, 80, 85–94. [Google Scholar] [CrossRef] [Green Version]
  56. Boyle, D.G.; Boyle, D.B.; Olsen, V.; Morgan, J.A.T.; Hyatt, A.D. Rapid quantitative detection of chytridiomycosis (Batrachochytrium dendrobatidis) in amphibian samples using real-time Taqman PCR assay. Dis. Aquat. Organ. 2004, 60, 141–148. [Google Scholar] [CrossRef]
  57. Hyatt, A.D.; Boyle, D.G.; Olsen, V.; Boyle, D.B.; Berger, L.; Obendorf, D.; Dalton, A.; Kriger, K.; Hero, M.; Hines, H.; et al. Diagnostic assays and sampling protocols for the detection of Batrachochytrium dendrobatidis. Dis. Aquat. Organ. 2007, 73, 175–192. [Google Scholar] [CrossRef] [PubMed]
  58. Kriger, K.M.; Hero, J.-M.; Ashton, K.J. Cost efficiency in the detection of chytridiomycosis using PCR assay. Dis. Aquat. Organ. 2006, 71, 149–154. [Google Scholar] [CrossRef] [PubMed]
  59. Picco, A.M.; Collins, J.P. Fungal and viral pathogen occurrence in Costa Rican amphibians. J. Herpetol. 2007, 41, 746–749. [Google Scholar] [CrossRef]
  60. Goldberg, C.S.; Hawley, T.J.; Waits, L.P. Local and regional patterns of amphibian chytrid prevalence on the Osa Peninsula, Costa Rica. Herpetol. Rev. 2009, 40, 309–311. [Google Scholar]
  61. Whitfield, S.M.; Geerdes, E.; Chacon, I.; Ballestero Rodriguez, E.; Jimenez, R.; Donnelly, M.; Kerby, J. Infection and co-infection by the amphibian chytrid fungus and ranavirus in wild Costa Rican frogs. Dis. Aquat. Organ. 2013, 104, 173–178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Zumbado-Ulate, H.; Bolaños, F.; Gutiérrez-Espeleta, G.; Puschendorf, R. Extremely low prevalence of Batrachochytrium dendrobatidis in frog populations from Neotropical dry forest of Costa Rica supports the existence of a climatic refuge from disease. EcoHealth 2014, 11, 593–602. [Google Scholar] [CrossRef] [PubMed]
  63. Zumbado-Ulate, H.; García-Rodríguez, A.; Vredenburg, V.T.; Searle, C.L. Infection with Batrachochytrium dendrobatidis is common in tropical lowland habitats: Implications for amphibian conservation. Ecol. Evol. 2019, 9, 4917–4930. [Google Scholar] [CrossRef] [PubMed]
  64. Saenz, D.; Adams, C.K.; Pierce, J.B.; Laurencio, D. Occurrence of Batrachochytrium dendrobatidis in an anuran community in the southeastern Talamanca region of Costa Rica. Herpetol. Rev. 2009, 40, 311–313. [Google Scholar]
  65. Abarca, J.G. Quitridiomicosis en Costa Rica: Aislamiento y Descripción de Cepas Circulantes del Patógeno y Análisis de la Microbiota del Hospedero como Posible Factor en la Incidencia de la Enfermedad. Master’s Thesis, Universidad de Costa Rica, Heredia, Costa Rica, 2018. [Google Scholar]
  66. Lips, K.R.; Diffendorfer, J.; Mendelson, J.R.; Sears, M.W. Riding the wave: Reconciling the roles of disease and climate change in amphibian declines. PLoS Biol. 2008, 6, e72. [Google Scholar] [CrossRef]
  67. Lips, K.R.; Burrowes, P.A.; Mendelson, J.R.; Parra-Olea, G. Amphibian declines in Latin America: Widespread population declines, extinctions, and impacts. Biotropica 2005, 37, 163–165. [Google Scholar] [CrossRef]
  68. Whitfield, S.M.; Lips, K.R.; Donnelly, M.A. Amphibian decline and conservation in Central America. Copeia 2016, 104, 351–379. [Google Scholar] [CrossRef]
  69. Arias, E.; Chaves, G. 140 years after William M. Gabb’s climb to Cerro Pico Blanco. Mesoam. Herpetol. 2014, 1, 176–180. [Google Scholar]
  70. Holdridge, L.R. Life Zone Ecology; Tropical Science Center: San Jose, Costa Rica, 1967. [Google Scholar]
  71. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019; Available online: https://www.r-project.org/ (accessed on 22 February 2019).
  72. Burnham, K.P.; Anderson, D.R. Multimodel inference: Understanding AIC and BIC in model selection. Sociol. Methods Res. 2004, 33, 261–304. [Google Scholar] [CrossRef]
  73. Ryan, M.; Lips, K.; Giermakowski, J.T. New species of Pristimantis (Anura: Terrarana: Strabomantinae) from lower Central America. J. Herpetol. 2010, 44, 193–200. [Google Scholar] [CrossRef]
  74. Batista, A.; Hertz, A.; Köhler, G.; Mebert, K.; Vesely, M. Morphological variation and phylogeography of frogs related to Pristimantis caryophyllaceus (Anura: Terrarana: Craugastoridae) in Panama. Salamandra 2014, 50, 155–171. [Google Scholar]
  75. Kubicki, B.; Salazar, S.; Puschendorf, R. A new species of glassfrog, genus Hyalinobatrachium (Anura: Centrolenidae), from the Caribbean foothills of Costa Rica. Zootaxa 2015, 3920, 069–084. [Google Scholar] [CrossRef]
  76. Arias, E.; Chaves, G.; Parra-Olea, G. A new species of Craugastor (Anura: Craugastoridae) from the montane rainforest of the Cordillera de Talamanca, Costa Rica. Phyllomedusa J. Herpetol. 2018, 17, 211–232. [Google Scholar] [CrossRef]
  77. Arias, E.; Chaves, G.; Crawford, A.J.; Parra-Olea, G. A new species of the Craugastor podiciferus species group (Anura: Craugastoridae) from the premontane forest of southwestern Costa Rica. Zootaxa 2016, 4132, 347–363. [Google Scholar] [CrossRef]
  78. Arias, E.; Hertz, A.; Parra-Olea, G. Taxonomic assessment of Craugastor podiciferus (Anura: Craugastoridae) in lower Central America with the description of two new species. Amphib. Reptile Conserv. 2019, 13, 173–197. [Google Scholar]
  79. Arias, E.; Chaves, G.; Salazar, S.; Salazar-Zúñiga, J.A.; García-Rodríguez, A. A new species of dink frog, genus Diasporus (Anura: Eleutherodactylidae), from the Caribbean foothills of the Cordillera de Talamanca, Costa Rica. Zootaxa 2019, 4609, 269–288. [Google Scholar] [CrossRef]
  80. Barquero, M.D.; Araya, M.F. First record of the Greenhouse frog, Eleutherodactylus planirostris (Anura: Eleutherodactylidae), in Costa Rica. Herpetol. Notes 2016, 9, 145–147. [Google Scholar]
  81. Kubicki, B.; Salazar, S. Discovery of the golden-eyed fringe-limbed treefrog, Ecnomiohyla bailarina (Anura: Hylidae), in the Caribbean foothills of southeastern Costa Rica. Mesoam. Herpetol. 2015, 2, 76–86. [Google Scholar]
  82. McCranie, J.R. Morphological and systematic comments on the Caribbean lowland population of Smilisca baudinii (Anura: Hylidae: Hylinae) in northeastern Honduras, with the resurrection of Hyla manisorum. Mesoam. Herpetol. 2017, 4, 513–526. [Google Scholar]
  83. Gray, A.R. Review of the genus Cruziohyla (Anura: Phyllomedusidae), with description of a new species. Zootaxa 2018, 4450, 401–426. [Google Scholar] [CrossRef] [PubMed]
  84. Kubicki, B.; Arias, E. A beautiful new yellow salamander, genus Bolitoglossa (Caudata: Plethodontidae), from the northeastern slopes of the Cordillera de Talamanca, Costa Rica. Zootaxa 2016, 4184, 329–346. [Google Scholar] [CrossRef] [PubMed]
  85. Boza-Oviedo, E.; Rovito, S.M.; Chaves, G.; Garcia-Rodriguez, A.; Artavia, L.G.; Bolaños, F.; Wake, D.B. Salamanders from the eastern Cordillera de Talamanca, Costa Rica, with descriptions of five new species (Plethodontidae: Bolitoglossa, Nototriton, and Oedipina) and natural history notes from recent expeditions. Zootaxa 2012, 3309, 36–61. [Google Scholar] [CrossRef]
  86. Arias, E.; Kubicki, B. A new moss salamander, genus Nototriton (Caudata: Plethodontidae), from the Cordillera de Talamanca, in the Costa Rica-Panama border region. Zootaxa 2018, 4369, 487–500. [Google Scholar] [CrossRef] [PubMed]
  87. Kubicki, B. A new species of salamander (Caudata: Plethodontidae: Oedipina) from the central Caribbean foothills of Costa Rica. Mesoam. Herpetol. 2016, 3, 819–840. [Google Scholar]
  88. Kubicki, B.; Arias, E. Vulcan’s Slender Caecilian, Caecilia volcani, in Costa Rica. Mesoam. Herpetol. 2017, 4, 488–492. [Google Scholar]
  89. Jiménez, R.; Alvarado, G. Craugastor escoces (Anura: Craugastoridae) reappears after 30 years: Rediscovery of an “extinct” Neotropical frog. Amphib. Reptil. 2017, 38, 257–259. [Google Scholar] [CrossRef]
  90. García-Rodríguez, A.; Chaves, G.; Benavides-Varela, C.; Puschendorf, R. Where are the survivors? Tracking relictual populations of endangered frogs in Costa Rica. Divers. Distrib. 2012, 18, 204–212. [Google Scholar] [CrossRef]
  91. Lewis, C.H.R.; Richards-Zawacki, C.L.; Ibáñez, R.; Luedtke, J.; Voyles, J.; Houser, P.; Gratwicke, B. Conserving Panamanian harlequin frogs by integrating captive-breeding and research programs. Biol. Conserv. 2019, 236, 180–187. [Google Scholar] [CrossRef]
  92. Searle, C.L.; Gervasi, S.S.; Hua, J.; Hammond, J.I.; Relyea, R.A.; Olson, D.H.; Blaustein, A.R. Differential host susceptibility to Batrachochytrium dendrobatidis, an emerging amphibian pathogen. Conserv. Biol. 2011, 25, 965–974. [Google Scholar] [CrossRef] [PubMed]
  93. Young, B.E.; Lips, K.R.; Reaser, J.K.; Ibáñez, R.; Salas, A.W.; Cedeño, J.R.; Coloma, L.A.; Ron, S.; La Marca, E.; Meyer, J.R.; et al. Population declines and priorities for amphibian conservation in Latin America. Conserv. Biol. 2001, 15, 1213–1223. [Google Scholar] [CrossRef]
  94. Martel, A.; Spitzen-van der Sluijs, A.; Blooi, M.; Bert, W.; Ducatelle, R.; Fisher, M.C.; Woeltjes, A.; Bosman, W.; Chiers, K.; Bossuyt, F.; et al. Batrachochytrium salamandrivorans sp. nov. causes lethal chytridiomycosis in amphibians. Proc. Natl. Acad. Sci. USA 2013, 110, 15325–15329. [Google Scholar] [CrossRef] [PubMed]
  95. Whitfield, S.M.; Kerby, J.; Gentry, L.R.; Donnelly, M.A. Temporal variation in infection prevalence by the amphibian chytrid fungus in three species of frogs at La Selva, Costa Rica. Biotropica 2012, 44, 779–784. [Google Scholar] [CrossRef]
  96. Whitfield, S.M.; Zoo Miami, FL, USA. Personal Communication, 2019.
  97. Whitfield, S.M.; Alvarado, G.; Abarca, J.; Zumbado-Ulate, H.; Zuñiga, I.; Wainwright, M.; Kerby, J. Differential patterns of Batrachochytrium dendrobatidis infection in relict amphibian populations following severe disease-associated declines. Dis. Aquat. Organ. 2017, 126, 33–41. [Google Scholar] [CrossRef]
  98. Puschendorf, R.; Hoskin, C.J.; Cashins, S.D.; McDonald, K.; Skerratt, L.F.; Vanderwal, J.; Alford, R.A. Environmental refuge from disease-driven amphibian extinction. Conserv. Biol. 2011, 25, 956–964. [Google Scholar] [CrossRef]
  99. Perez, R.; Richards-Zawacki, C.L.; Krohn, A.R.; Robak, M.; Griffith, E.J.; Ross, H.; Gratwicke, B.; Ibanez, R.; Voyles, J. Field surveys in Western Panama indicate populations of Atelopus varius frogs are persisting in regions where Batrachochytrium dendrobatidis is now enzootic. Amphib. Reptile Conserv. 2014, 8, 30–35. [Google Scholar]
  100. Woodhams, D.C.; Kilburn, V.L.; Reinert, L.K.; Voyles, J.; Medina, D.; Ibáñez, R.; Hyatt, A.D.; Boyle, D.G.; Pask, J.D.; Green, D.M.; et al. Chytridiomycosis and amphibian population declines continue to spread eastward in Panama. EcoHealth 2008, 5, 268–274. [Google Scholar] [CrossRef]
  101. Kilburn, V.L.; Ibáñez, R.; Sanjur, O.; Bermingham, E.; Suraci, J.P.; Green, D.M. Ubiquity of the pathogenic chytrid fungus, Batrachochytrium dendrobatidis, in anuran communities in Panamá. EcoHealth 2010, 7, 537–548. [Google Scholar] [CrossRef] [PubMed]
  102. Acosta-Chaves, V.J.; Madrigal-Elizondo, V.; Chaves, G.; Morera-Chacón, B.; Garcia-Rodriguez, A.; Bolaños, F. View of shifts in the diversity of an amphibian community from a premontane forest of San Ramón, Costa Rica. Rev. Biol. Trop. 2019, 67, 259–273. [Google Scholar]
  103. Garner, T.W.J.; Walker, S.; Bosch, J.; Leech, S.; Marcus Rowcliffe, J.; Cunningham, A.A.; Fisher, M.C. Life history tradeoffs influence mortality associated with the amphibian pathogen Batrachochytrium dendrobatidis. Oikos 2009, 118, 783–791. [Google Scholar] [CrossRef]
  104. Kolby, J.E.; Ramirez, S.D.; Berger, L.; Richards-Hrdlicka, K.L.; Jocque, M.; Skerratt, L.F. Terrestrial dispersal and potential environmental transmission of the amphibian chytrid fungus (Batrachochytrium dendrobatidis). PLoS ONE 2015, 10, e0125386. [Google Scholar] [CrossRef] [PubMed]
  105. Ryan, M.J.; Lips, K.R.; Eichholz, M.W. Decline and extirpation of an endangered Panamanian stream frog population (Craugastor punctariolus) due to an outbreak of chytridiomycosis. Biol. Conserv. 2008, 141, 1636–1647. [Google Scholar] [CrossRef]
  106. Campbell, J.A.; Savage, J.M. Taxonomic reconsideration of Middle American frogs of the Eleutherodactylus rugulosus group (Anura: Leptodactylidae): A reconnaissance of subtle nuances among frogs. Herpetol. Monogr. 2000, 14, 186–292. [Google Scholar] [CrossRef]
  107. Lips, K.R.; Reeve, J.D.; Witters, L.R. Ecological traits predicting amphibian population declines in Central America. Conserv. Biol. 2003, 17, 1078–1088. [Google Scholar] [CrossRef]
  108. Köhler, G.; Batista, A.; Carrizo, A.; Hertz, A. Field notes on Craugastor azueroensis (Savage, 1975) (Amphibia: Anura: Craugastoridae). Herpetol Notes 2012, 5, 157. [Google Scholar]
  109. Hanken, J.; Wake, D.B.; Savage, J.M. A solution to the large black salamander problem (genus Bolitoglossa) in Costa Rica and Panamá. Copeia 2005, 2005, 227–245. [Google Scholar] [CrossRef]
  110. Sunyer, J.; Wake, D.B.; Obando, L. Distributional data for Bolitoglossa (Amphibia, Caudata, Plethodontidae) from Nicaragua and Costa Rica. Herpetol. Rev. 2012, 43, 564–568. [Google Scholar]
  111. AlMutairi, B.S.; Grossmann, I.; Small, M.J. Climate model projections for future seasonal rainfall cycle statistics in Northwest Costa Rica. Int. J. Climatol. 2019, 39, 2933–2946. [Google Scholar] [CrossRef]
  112. Grenyer, R.; Orme, C.D.L.; Jackson, S.F.; Thomas, G.H.; Davies, R.G.; Davies, T.J.; Jones, K.E.; Olson, V.A.; Ridgely, R.S.; Rasmussen, P.C.; et al. Global distribution and conservation of rare and threatened vertebrates. Nature 2006, 444, 93–96. [Google Scholar] [CrossRef] [PubMed]
  113. Woodhams, D.C.; Bosch, J.; Briggs, C.J.; Cashins, S.; Davis, L.R.; Lauer, A.; Muths, E.; Puschendorf, R.; Schmidt, B.R.; Sheafor, B.; et al. Mitigating amphibian disease: Strategies to maintain wild populations and control chytridiomycosis. Front. Zool. 2011, 8, 8. [Google Scholar] [CrossRef] [PubMed]
  114. Garner, T.W.J.; Schmidt, B.R.; Martel, A.; Pasmans, F.; Muths, E.; Cunningham, A.A.; Weldon, C.; Fisher, M.C.; Bosch, J. Mitigating amphibian chytridiomycosis in nature. Philos. Trans. R. Soc. B Biol. Sci. 2016, 371, 20160207. [Google Scholar] [CrossRef] [PubMed]
Diversity 11 00129 g001 550
Figure 1. Map of 20 survey sites across Costa Rica. Sites are color-coded by herpetological province.
Figure 1. Map of 20 survey sites across Costa Rica. Sites are color-coded by herpetological province.
Diversity 11 00129 g001
Diversity 11 00129 g002 550
Figure 2. Map of Costa Rica showing (a) amphibian species richness for each herpetological province and percentage of amphibian species classified as (b) data deficient, (c) least concerned, and (d) threatened categories (near threatened, vulnerable, endangered, critically endangered, and extinct in the wild) for each herpetological province according to the Red List of Threatened Species from the International Union of Conservation of Nature (IUCN). Symbology: CL—Caribbean Lowlands, MSCC—Mountain Slopes and Cordillera Central, PN—Pacific Northwest, and PS—Pacific Southwest.
Figure 2. Map of Costa Rica showing (a) amphibian species richness for each herpetological province and percentage of amphibian species classified as (b) data deficient, (c) least concerned, and (d) threatened categories (near threatened, vulnerable, endangered, critically endangered, and extinct in the wild) for each herpetological province according to the Red List of Threatened Species from the International Union of Conservation of Nature (IUCN). Symbology: CL—Caribbean Lowlands, MSCC—Mountain Slopes and Cordillera Central, PN—Pacific Northwest, and PS—Pacific Southwest.
Diversity 11 00129 g002
Diversity 11 00129 g003 550
Figure 3. Mean prevalence of infection with Batrachochytrium dendrobatidis (Bd) in amphibian assemblages at four herpetological provinces in Costa Rica (with 95% binomial CI). Means followed by a common letter are not significantly different according to the Tukey’s honestly significant difference (HSD) test at the 5% level of significance. The plot does not display results for Cordillera de Talamanca because no sampling has been conducted for Bd in that province. Symbology: CL—Caribbean Lowlands, MSCC—Mountain Slopes and Cordillera Central, PN—Pacific Northwest, and PS—Pacific Southwest.
Figure 3. Mean prevalence of infection with Batrachochytrium dendrobatidis (Bd) in amphibian assemblages at four herpetological provinces in Costa Rica (with 95% binomial CI). Means followed by a common letter are not significantly different according to the Tukey’s honestly significant difference (HSD) test at the 5% level of significance. The plot does not display results for Cordillera de Talamanca because no sampling has been conducted for Bd in that province. Symbology: CL—Caribbean Lowlands, MSCC—Mountain Slopes and Cordillera Central, PN—Pacific Northwest, and PS—Pacific Southwest.
Diversity 11 00129 g003
Diversity 11 00129 g004 550
Figure 4. (a) Mean prevalence of infection with Batrachochytrium dendrobatidis (Bd) in amphibian assemblages (with 95% binomial CI) according to the foraging–reproduction habitat index (FRHI); (b) box plots with whiskers and notches that describe infection intensity for the 351 individuals where Bd was quantified using qPCR in Costa Rica between 2000–2018 according to the foraging–reproduction habitat index (FRHI). The box displays the inter-quantile range (25th–75th percentiles) with a center line representing the median (50th percentile). Notches show the median confidence region, and whiskers display the highest and lowest points. Means followed by a common letter are not significantly different according to the Tukey’s honestly significant difference (HSD) test at the 5% level of significance. Symbology: First letter represents development: (I) indirect or (D) direct-developing amphibians. Second letter represents foraging habitat: terrestrial (T), arboreal (A), pond/puddle-dwellers (P), stream-breeders (R), and phytotelma (F). Third letter represents reproductive habitat: terrestrial (T), arboreal (A), pond/puddle-dwellers (P), stream-breeders (R), and phytotelma (F).
Figure 4. (a) Mean prevalence of infection with Batrachochytrium dendrobatidis (Bd) in amphibian assemblages (with 95% binomial CI) according to the foraging–reproduction habitat index (FRHI); (b) box plots with whiskers and notches that describe infection intensity for the 351 individuals where Bd was quantified using qPCR in Costa Rica between 2000–2018 according to the foraging–reproduction habitat index (FRHI). The box displays the inter-quantile range (25th–75th percentiles) with a center line representing the median (50th percentile). Notches show the median confidence region, and whiskers display the highest and lowest points. Means followed by a common letter are not significantly different according to the Tukey’s honestly significant difference (HSD) test at the 5% level of significance. Symbology: First letter represents development: (I) indirect or (D) direct-developing amphibians. Second letter represents foraging habitat: terrestrial (T), arboreal (A), pond/puddle-dwellers (P), stream-breeders (R), and phytotelma (F). Third letter represents reproductive habitat: terrestrial (T), arboreal (A), pond/puddle-dwellers (P), stream-breeders (R), and phytotelma (F).
Diversity 11 00129 g004
Table
Table 1. Summary of studies where Batrachochytrium dendrobatidis (Bd) was detected in multi-species amphibian assemblages using conventional PCR and quantitative PCR (qPCR) in Costa Rica between 2005–2018. The table shows surveyed localities, herpetological province, sampling period, percentage of infection, and Holdridge’s altitudinal belt. Symbology: CL—Caribbean Lowlands, MSCC—Montane Slopes and Cordillera Central, PN—Pacific Northwest, and PS—Pacific Southwest.
Table 1. Summary of studies where Batrachochytrium dendrobatidis (Bd) was detected in multi-species amphibian assemblages using conventional PCR and quantitative PCR (qPCR) in Costa Rica between 2005–2018. The table shows surveyed localities, herpetological province, sampling period, percentage of infection, and Holdridge’s altitudinal belt. Symbology: CL—Caribbean Lowlands, MSCC—Montane Slopes and Cordillera Central, PN—Pacific Northwest, and PS—Pacific Southwest.
Study Site (Elevation m) and Herpetological ProvincesSampling Period% of Infection (n Sampled)Altitudinal BeltReference
Monteverde (1400–2000), MSCCJuly 200512.2 (41)Lower montane[59]
San Vito de Coto Brus (1120–1385), MSCC9.3 (43)Premontane
Rincón de Osa (0–100), PSMay–June 20060.1 (91)Lowland[60]
Piro (0–100), PS0.0 (62)Lowland
Corcovado (0–100), PS0.1 (25)Lowland
Kekoldi (0–100), CLJanuary 20087.9 (126)Lowland[64]
La Virgen de Sarapiquí (0–200), CLJanuary–March 201121.3 (253)Lowland[61]
Santa Elena Peninsula (0–200), PNJanuary–March 2007–20080.0 (310)Lowland[62]
Santa Rosa (0–200), PN9.0 (100)Lowland
Punta Banco–Burica (0–100), PSNovember–December 201168.6 (35)Lowland[63]
Rincón de Osa (0–100), PS0.0 (25)Lowland
Puerto Viejo de Sarapiquí (0–200), CL67.4 (144)Lowland
Guayacán de Siquirres (400–600), CL47.9 (144)Lowland
San Vito de Coto Brus (1120–1385), MSCCUnknown/not indicated10.5 (19)Lowland[65]
Punta Banco–Burica (0–100), PS0.0 (20)Lowland
Guayacán de Siquirres (400–600), CL5.3 (19)Lowland
San Rafael de Heredia (1800), MSCC66.7 (15)Lower montane
Santo Domingo de Heredia (1000–1200), MSCC45.5 (11)Premontane
Las Tablas (1350), MSCC28.6 (14)Lower montane
Table
Table 2. Categories and taxonomic examples of the foraging–reproduction habitat index (FRHI) that we developed for this study to analyze current prevalence of Batrachochytrium dendrobatidis across taxonomic groups. Symbology: First letter represents development: (I) indirect- or (D) direct-developing amphibians. Second letter represents foraging habitat: terrestrial (T), arboreal (A) pond/puddle-dwellers (P), stream-breeders (R), and phytotelma (F). Third letter represents reproductive habitat: terrestrial (T), arboreal (A), pond/puddle-dwellers (P), stream-breeders (R), and phytotelma (F).
Table 2. Categories and taxonomic examples of the foraging–reproduction habitat index (FRHI) that we developed for this study to analyze current prevalence of Batrachochytrium dendrobatidis across taxonomic groups. Symbology: First letter represents development: (I) indirect- or (D) direct-developing amphibians. Second letter represents foraging habitat: terrestrial (T), arboreal (A) pond/puddle-dwellers (P), stream-breeders (R), and phytotelma (F). Third letter represents reproductive habitat: terrestrial (T), arboreal (A), pond/puddle-dwellers (P), stream-breeders (R), and phytotelma (F).
FRHISpeciesTaxonomic Group (Examples)
DAA2Diasporus spp. (dink frogs, e.g., Diasporus diastema)
DAT5Pristimantis spp. (rain frogs, e.g., Pristimantis cerasinus)
DRT3Craugastor punctariolus clade (robber frogs, e.g., Craugastor taurus)
C. fitzingeri (Pacific side)
DTT13Craugastor spp. (leaf-litter frogs. e.g., Craugastor bransfordi)
C. fitzingeri (Caribbean side)
Plethodontidae (e.g., Oedipina gracilis)
IAP17Hylidae (pond-breeding treefrogs, e.g., Boana rufitela)
IAR15Centrolenidae (glass frogs, e.g., Teratohyla pulverata)
Hylidae (stream-breeding treefrogs, e.g., Duellmanohyla rufioculis)
ITF4Dendrobatidae (Poison-dart frogs, e.g., Oophaga pumilio)
ITP12Leptodactylidae (Leptodactylid frogs, e.g., Leptodactylus melanonotus)
Microhylidae (sheep frogs, e.g., Hypopachus variolosus)
Ranidae (Ranid frogs, e.g., Lithobates forreri)
Bufonidae (toads, e.g., Incilius coccifer)
ITR7Bufonidae (river toads, e.g., Rhaebo haematiticus)
Table
Table 3. List of new additions to the updated list of amphibians of Costa Rica.
Table 3. List of new additions to the updated list of amphibians of Costa Rica.
Order FamilySpeciesSource
AnuraCentrolenidaeHyalinobatrachium dianae[75]
CraugastoridaeCraugastor aenigmaticus[76]
Craugastor gabbi[77]
Craugastor zunigai[78]
EleutherodactylidaeDiasporus amirae[79]
Eleutherodactylus planirostris *[80]
HylidaeEcnomiohyla bailarina[81]
Ecnomiohyla veraguensisUnpublished
Smilisca manisorum[82]
PhyllomedusidaeCruziohyla sylviae[83]
CaudataPlethodontidaeBolitoglossa aurae[84]
Bolitoglossa aureogularis[85]
Bolitoglossa kamuk[85]
Bolitoglossa pygmaeaUnpublished
Bolitoglossa splendida[85]
Nototriton costaricense[86]
Nototriton matama[85]
Oedipina berlini[87]
Oedipina nimaso[85]
GymnophionaCaeciliidaeCaecilia volcani[88]
* Eleutherodactylus planirostris is an invasive species that have been found in the Caribbean Lowlands of Costa Rica.
Table
Table 4. Infection intensity in the 351 individuals where Batrachochytrium dendrobatidis (Bd) was quantified using qPCR in Costa Rica between 2000–2018. For every species, the table shows the foraging–reproduction habitat index (FRHI), the number of Bd positive swabs, the average (SE), and log10(SE) of genomic equivalents of Bd zoospores quantified per species estimated from the Bd + swabs. Symbology: First letter represents development: (I) indirect or (D) direct-developing amphibians. Second letter represents foraging habitat: terrestrial (T), arboreal (A), pond/puddle-dwellers (P), stream-breeders (R), and phytotelma (F). Third letter represents reproductive habitat: terrestrial (T), arboreal (A), pond/puddle-dwellers (P), stream-breeders (R), and phytotelma (F).
Table 4. Infection intensity in the 351 individuals where Batrachochytrium dendrobatidis (Bd) was quantified using qPCR in Costa Rica between 2000–2018. For every species, the table shows the foraging–reproduction habitat index (FRHI), the number of Bd positive swabs, the average (SE), and log10(SE) of genomic equivalents of Bd zoospores quantified per species estimated from the Bd + swabs. Symbology: First letter represents development: (I) indirect or (D) direct-developing amphibians. Second letter represents foraging habitat: terrestrial (T), arboreal (A), pond/puddle-dwellers (P), stream-breeders (R), and phytotelma (F). Third letter represents reproductive habitat: terrestrial (T), arboreal (A), pond/puddle-dwellers (P), stream-breeders (R), and phytotelma (F).
Species (FRHI)Bd + SwabsBd Load Average (SE)Log10 Bd Load Average (SE)
Agalychnis callidryas (IAP)48.19 (3.81)0.77 (0.21)
Agalychnis spurrelli (IAP)539.83 (32.47)1.10 (0.30)
Boana rufitela (IAP)88.41 (4.12)0.53 (0.23)
Bolitoglossa colonnea (DTT)11.83 (0.00)0.26 (0.00)
Cochranella granulosa (IAR)13.95 (0.00)0.60 (0.00)
Craugastor bransfordi (DTT)231007.06 (483.50)1.78 (0.25)
Craugastor crassidigitus (DTT)51636.74 (1583.43)1.69 (0.64)
Craugastor fitzingeri (DTT, DRT)44951.48 (310.27)1.97 (0.17)
Craugastor megacephalus (DTT)10.62 (0.00)−0.21 (0.00)
Craugastor mimus (DTT)9125.48 (74.69)1.01 (0.46)
Craugastor ranoides (DRT)3187.40 (174.13)1.65 (0.55)
Craugastor stejnegerianus (DTT)22.18 (0.95)0.29 (0.20)
Craugastor taurus (DRT)1211,632.50 (6564.67)2.51 (0.41)
Dendropsophus ebraccatus (IAP)34315.85 (194.09)1.00 (0.18)
Diasporus diastema (DAA)214.44 (4.12)1.14 (0.13)
Duellmanohyla rufioculis (IAR)13.65 (0.00)0.56 (0.00)
Engystomops pustulosus (ITP)1134.83 (13.26)1.11 (0.3)
Espadarana prosoblepon (IAR)33691.59 (3684.73)1.06 (1.75)
Hyalinobatrachium colymbiphyllum (IAR)10.01 (0.00)−2.00 (0.00)
Hyalinobatrachium valerioi (IAR)28.38 (2.04)0.91 (0.11)
Incilius melanochlorus (ITR)223.27 (19.93)1.08 (0.56)
Leptodactylus melanonotus (ITP)411.86 (0.66)1.07 (0.02)
Leptodactylus poecilochilus (ITP)11073.45 (0.00)3.03 (0.00)
Leptodactylus savagei (ITP)133.49 (0.00)1.52 (0.00)
Lithobates forreri (ITP)2569.24 (241.10)2.71 (0.20)
Lithobates warszewitschii (ITR)14978.92 (801.60)1.47 (0.31)
Oophaga granulifera (ITF)923.92 (11.31)1.20 (0.11)
Oophaga pumilio (ITF)341765.81 (778.67)1.71 (0.25)
Pristimantis cerasinus (DAT)914.82 (10.97)0.47 (0.32)
Pristimantis ridens (DAT)748.37 (32.34)0.69 (0.50)
Rhaebo haematiticus (ITR)22239.20 (178.56)0.70 (0.26)
Scinax boulengeri (IAP)1195.20 (0.00)2.29 (0.00)
Scinax elaeochroa (IAP)51384.15 (1350.51)1.78 (0.58)
Smilisca phaeota (IAP)437.25 (19.15)1.44 (0.18)
Smilisca sordida (IAP)4614.96 (9.27)0.24 (0.16)
Teratohyla pulverata (IAR)234.53 (22.92)1.41 (0.35)
Teratohyla spinosa (IAR)5937.99 (825.54)1.90 (0.57)
Tlalocohyla loquax (IAP)11144.66 (107.28)1.22 (0.30)

Share and Cite

MDPI and ACS Style

Zumbado-Ulate, H.; Nelson, K.N.; García-Rodríguez, A.; Chaves, G.; Arias, E.; Bolaños, F.; Whitfield, S.M.; Searle, C.L. Endemic Infection of Batrachochytrium dendrobatidis in Costa Rica: Implications for Amphibian Conservation at Regional and Species Level. Diversity 2019, 11, 129. https://doi.org/10.3390/d11080129

AMA Style

Zumbado-Ulate H, Nelson KN, García-Rodríguez A, Chaves G, Arias E, Bolaños F, Whitfield SM, Searle CL. Endemic Infection of Batrachochytrium dendrobatidis in Costa Rica: Implications for Amphibian Conservation at Regional and Species Level. Diversity. 2019; 11(8):129. https://doi.org/10.3390/d11080129

Chicago/Turabian Style

Zumbado-Ulate, Héctor, Kiersten N. Nelson, Adrián García-Rodríguez, Gerardo Chaves, Erick Arias, Federico Bolaños, Steven M. Whitfield, and Catherine L. Searle. 2019. "Endemic Infection of Batrachochytrium dendrobatidis in Costa Rica: Implications for Amphibian Conservation at Regional and Species Level" Diversity 11, no. 8: 129. https://doi.org/10.3390/d11080129

APA Style

Zumbado-Ulate, H., Nelson, K. N., García-Rodríguez, A., Chaves, G., Arias, E., Bolaños, F., Whitfield, S. M., & Searle, C. L. (2019). Endemic Infection of Batrachochytrium dendrobatidis in Costa Rica: Implications for Amphibian Conservation at Regional and Species Level. Diversity, 11(8), 129. https://doi.org/10.3390/d11080129

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop