Data from: Rapid emergence of an amphibian pathogen coincided with historic amphibian declines in the Neotropics
Data files
Apr 09, 2026 version files 9.49 MB
-
Ecuador_historicBd.zip
9.48 MB
-
README.md
7.71 KB
Abstract
The emergence of the fungal pathogen Batrachochytrium dendrobatidis (Bd) in South America has been attributed to multiple introductions followed by a bidirectional latitudinal spread along the primary Andean cordilleras. In Ecuador, this hypothesis is supported by anecdotal reports of population declines in Bd-susceptible genera such as Telmatobius and Atelopus during the mid-1980s. To evaluate this pattern, we combined published Bd records with retrospective screening of museum specimens collected between 1950 and 2010, along with surveys of contemporary populations, to characterize the spatial, temporal, and host-specific dynamics of Bd emergence in Ecuador. We compiled a comprehensive database of Bd records from Ecuadorian anurans comprising over 7,000 samples and analyzed patterns of emergence using generalized linear mixed models and spatial clustering techniques. The dataset was also used to estimate Bd prevalence across four representative genera: Atelopus, Telmatobius, Gastrotheca, and Pristimantis. Additionally, Bayesian methods were applied to assess the probability that Bd was present in Ecuador prior to its first documented detection. Our results overturn the prevailing narrative of Bd emergence in the Neotropics and provide a rare, long-term view of an invasive pathogen’s trajectory from introduction to persistence. We found that: (1) Bd likely represents a novel pathogen to Ecuador; (2) Bd was first detected in 1974 and had become widespread by 1980; (3) there is no evidence supporting a bidirectional wave of spread; and (4) the genera Atelopus, Telmatobius, and Gastrotheca experienced epizootic disease dynamics, with declines peaking in the 1990s. We found that Bd is likely a novel pathogen that was introduced in the early 1970s and peaked in prevalence during the late 1980s and 1990s. The available evidence supports the conclusion that an outbreak of chytridiomycosis may lead to catastrophic declines in Ecuadorian amphibians. Disease dynamics have since subsided to enzootic coexistence in extant amphibian communities.
Dataset DOI: 10.5061/dryad.hqbzkh1xf
Description of the data and file structure
We assembled a database of Bd records in Ecuadorian anurans (>7000 samples) and assessed its spatiotemporal emergence using generalized linear mixed models and spatial clustering approaches. Data were also used to identify Bd prevalence within four representative genera of Ecuadorian anurans (Atelopus, Telmatobius, Gastrotheca, and Pristimantis). We also used Bayesian approaches to test the probability of the presence of Bd in Ecuador prior to the first positive record.
Files and variables
File: Ecuador_historicBd.zip
Description:
Historic_Ecuador_Bd.R contains the R script for reproducing all the figures and analyses within the manuscript
Declines.csv contains all the reported data on Ecuadorian amphibian declines. Data sourced from https://doi.org/10.1371/journal.pbio.0060072
- Species or Events- Name of reported decline event; in some cases, more than one decline is reported per species.
- DOD- date of decline
- LYO- Last year observed
- EDOD- Earliest decade of decline
- LDOD- latest decade of decline
- Status- Declining, data-deficient, or unknown
- Elevation- elevation above sea level in metres
- Longitude- decimal degrees
- Latitude- decimal degrees
- Source- literature sources, as cited from Lips et al 2008
Declines_Update_DeOD.csv contains all the data on amphibian declines. Same as above, but with an added decade column.
- Species or Events- Name of reported decline event; in some cases, more than one decline is reported per species.s
- DOD- date of decline
- DeOD- Decade of decline
- LYO- Last year observed
- EDOD- Earliest decade of decline
- LDOD- latest decade of decline
- Status- Declining, data-deficient, or unknown
- Elevation- elevation above sea level in metres
- Longitude- decimal degrees
- Latitude- decimal degrees
- Source- literature sources, as cited from Lips et al 2008
EcuadorChytrid.csv raw data on all amphibian Bd data (literature, contemporary swabs, and swabs from museum specimens), including Life stage, decade, species, genus, GPS coordinates, and diagnostic method
- Rownumber- numeric
- Identification- Specimen code or sample code from source material
- Year- year of sample collection (living animals) or date of specimen collection (museum specimens)
- Decade- as above, but rounded to a decade
- Life_stage- Adult (postmetamorphic) or tadpole (premetamorphic)
- Familia- amphibian family
- SPECIES- amphibian species
- Genero- amphibian genus
- Repr_mode- reproductive mode of amphibian species
- LOCALITY- site name
- Latitude- decimal degrees
- Longitude- decimal degrees
- Source- origin of data from authors or published material
- Diagnosis- Bd diagnosis method
- Binary- disease diagnostic result (positive or negative)
- Source Type- from living animal (contemporary), museum specimens (archived specimen), or Literature
- Field series- sample code from QCAZ collection
Satscan folder: Contains all inputs/outputs from Satscan analyses reported within the manuscript
All_SatScan_Clusters.csv- data on all clusters generated by satscan analysis
- ANALYSIS- space or space-time models in SatScan
- DECADE- samples grouped by decade, 1950s-2010s
- Cluster- Cluster ID number for each decade
- Cluster_type- Low or High rate disease cluster
- LOC_ID- ID number assigned to cluster
- LATITUDE- latitude of cluster in decimal degrees
- LONGITUDE- longitude of cluster in decimal degrees
- RADIUS- radius of cluster (in km)
- Span- maximum cluster size (in km)
- START_DATE- Start date of scanning window
- END_DATE- End date of scanning window
- NUMBER_LOC- number of locations within cluster
- LLR- Log-likelihood ratio
- P_VALUE- significance of cluster
- Observed- number of observed cases
- Expected- number of expected cases
- ODE- Observed/Expected ratio
- REL_RISK- relative risk of infection within cluster
- POPULATION- population within cluster
- GINI_CLUST- significant GINI coefficient? TRUE or FALSE
SigClusters.csv- data on all clusters generated by satscan analysis
- ANALYSIS- space or space-time models in SatScan
- DECADE- samples grouped by decade, 1950s-2010s
- Cluster- Cluster ID number for each decade
- Cluster_type- Low or High rate disease cluster
- LOC_ID- ID number assigned to cluster
- LATITUDE- latitude of cluster in decimal degrees
- LONGITUDE- longitude of cluster in decimal degrees
- RADIUS- radius of cluster (in km)
- Span- maximum cluster size (in km)
- START_DATE- Start date of scanning window
- END_DATE- End date of scanning window
- NUMBER_LOC- number of locations within cluster
- LLR- Log-likelihood ratio
- P_VALUE- significance of cluster
- Observed- number of observed cases
- Expected- number of expected cases
- ODE- Observed/Expected ratio
- REL_RISK- relative risk of infection within cluster
- POPULATION- population within cluster
- GINI_CLUST- significant GINI coefficient? TRUE or FALSE
SubDirectory: Space
Subfolder: Cases
CasesR1_1970.csv, CasesR1_1980.csv, CasesR1_1990.csv, CasesR1_2000.csv, CasesR1_2010.cs.v Files with all Bd cases used in Space models in Satscan
- Rownumber
- ID2- site name
- Decade- Decade of detections
- Number of cases- number of Bd cases
Subfolder: Controls
ControlsR1_1970.csv, ControlsR1_1980.csv, ControlsR1_1990.csv, ControlsR1_2000.csv, ControlsR1_2010.c.sv Files with all controls used in Space models in Satscan
- Rownumber
- ID2- site name
- Decade- Decade of detections
- Number of cases- number of controls
Subfolder: Coordinates
CoordinatesR1_1970.csv, CoordinatesR1_1980.csv, CoordinatesR1_1990.csv, CoordinatesR1_2000.csv, CoordinatesR1_2010.csv. Files with all Coordinates used in Space models in Satscan
- Rownumber
- ID2- site name
- Longitude- decimal degrees
- Latitude- decimal degrees
Subfolder: Outputs
Contains 1970, 1980, 1990, 2000, 2010 folders, all containing output files from Satscan
See details on https://www.satscan.org/techdoc.html for contents of output files.
SubDirectory: SpaceTime
Subfolder: Cases
CasesR1ST_1970.csv, CasesR1ST_1980.csv, CasesR1ST_1990.csv, CasesR1ST_2000.csv, CasesR1ST_2010.csv. Files with all Bd cases used in SpaceTime models in Satscan
- Rownumber
- ID2- site name
- Decade- Decade of detections
- Number of cases- number of Bd cases
Subfolder: Controls
ControlsR1ST_1970.csv, ControlsR1ST_1980.csv, ControlsR1ST_1990.csv, ControlsR1ST_2000.csv, ControlsR1ST_2010.csv. Files with all controls used in SpaceTime models in Satscan
- Rownumber
- ID2- site name
- Decade- Decade of detections
- Number of cases- number of controls
Subfolder: Coordinates
CoordinatesR1ST_1970.csv, CoordinatesR1ST_1980.csv, CoordinatesR1ST_1990.csv, CoordinatesR1ST_2000.csv, CoordinatesR1ST_2010.cs.v Files with all Coordinates used in SpaceTime models in Satscan
- Rownumber
- ID2- site name
- Longitude- decimal degrees
- Latitude- decimal degrees
Subfolder: Outputs
Contains 1970, 1980, 1990, 2000, 2010, all folders containing output files from Satscan
See details on https://www.satscan.org/techdoc.html for contents of output files.
Code/software
All data were processed in R, and spatial clustering was performed in SatScan.
