Amphibian diversity in Serranía de Majé, an isolated mountain range in eastern Panamá.

Eastern Panamá is within the Mesoamerican biodiversity hotspot and supports an understudied amphibian fauna. Here we characterize the amphibian diversity across an elevational gradient in one of the least studied mountain ranges in eastern Panamá, Serranía de Majé. A total of 38 species were found, which represent 17% of all species reported for Panamá. Based on expected richness function and individual-based rarefaction curves, it is estimated that this is an underestimate and that at least 44 amphibian species occur in this area. Members of all three amphibian orders were encountered, represented by ten families and 22 genera, including five species endemic to Central America. Estimated species richness decreased with elevation, and the mid-elevation site supported both lowland and highland species. Our study provides a baseline for understanding the distribution pattern of amphibians in Panamá, for conservation efforts, and for determining disease-induced changes in amphibian communities.

Eco-evolutionary rescue promotes host-pathogen coexistence.

Emerging infectious pathogens are responsible for some of the most severe host mass mortality events in wild populations. Yet, effective pathogen control strategies are notoriously difficult to identify, in part because quantifying and forecasting pathogen spread and disease dynamics is challenging. Following an outbreak, hosts must cope with the presence of the pathogen, leading to host-pathogen coexistence or extirpation. Despite decades of research, little is known about host-pathogen coexistence post-outbreak when low host abundances and cryptic species make these interactions difficult to study. Using a novel disease-structured N-mixture model, we evaluate empirical support for three host-pathogen coexistence hypotheses (source-sink, eco-evolutionary rescue, and spatial variation in pathogen transmission) in a Neotropical amphibian community decimated by Batrachochytrium dendrobatidis (Bd) in 2004. During 2010-2014, we surveyed amphibians in Parque Nacional G. D. Omar Torríjos Herrera, Coclé Province, El Copé, Panama. We found that the primary driver of host-pathogen coexistence was eco-evolutionary rescue, as evidenced by similar amphibian survival and recruitment rates between infected and uninfected hosts. Average apparent monthly survival rates of uninfected and infected hosts were both close to 96%, and the expected number of uninfected and infected hosts recruited (via immigration/reproduction) was less than one host per disease state per 20-m site. The secondary driver of host-pathogen coexistence was spatial variation in pathogen transmission as we found that transmission was highest in areas of low abundance but there was no support for the source-sink hypothesis. Our results indicate that changes in the host community (i.e., through genetic or species composition) can reduce the impacts of emerging infectious disease post-outbreak. Our disease-structured N-mixture model represents a valuable advancement for conservation managers trying to understand underlying host-pathogen interactions and provides new opportunities to study disease dynamics in remnant host populations decimated by virulent pathogens.

External Reinfection of a Fungal Pathogen Does not Contribute to Pathogen Growth.

Chytridiomycosis is an emerging infectious disease of amphibians caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), which has led to devastating declines in amphibian populations worldwide. Current theory predicts that Bd infections are maintained through both reproduction on the host's skin and reinfection from sources outside of the host. To investigate the importance of external reinfection on pathogen burden, we infected captive-bred individuals of the highly susceptible Panamanian Golden Frog, Atelopus glyphus, and wild-caught glass frogs, Espadarana prosoblepon, with Bd. We housed the animals in one of three treatments: individually, in heterospecific pairs, and in conspecific pairs. For 8 weeks, we measured the Bd load and shedding rate of all frogs. We found that Atelopus had high rates of increase in both Bd load and shedding rate, but pathogen growth rates did not differ among treatments. The infection intensity of Espadarana co-housed with Atelopus was indistinguishable from those housed singly and those in conspecific pairs, despite being exposed to a large external source of Bd zoospores. Our results indicate that Bd load in both species is driven by pathogen replication within an individual, with reinfection from outside the host contributing little to the amplification of host fungal load.

Altitudinal distribution and advertisement call of Colostethus latinasus (Amphibia: Dendrobatidae), endemic species from eastern Panama and type species of Colostethus , with a molecular assessment of similar sympatric species.

We conducted a molecular assessment of Colostethus-like frogs along an elevational gradient in the Serranía de Pirre, above Santa Cruz de Cana, eastern Panama, aiming to establish their species identity and to determine the altitudinal distribution of C. latinasus. Our findings confirm the view of C. latinasus as an endemic species restricted to the highlands of this mountain range, i.e., 1350-1475 m.a.s.l., considered to be type locality of this species. We described the advertisement call of C. latinasus that consists of a series of 4-18 single, short and relatively loud "peep"-like notes given in rapid succession, and its spectral and temporal features were compared with calls of congeneric species. For the first time, DNA sequences from C. latinasus were obtained, since previously reported sequences were based on misidentified specimens. This is particularly important because C. latinasus is the type species of Colostethus, a genus considered paraphyletic according to recent phylogenetic analyses based on molecular data.

Disassembly of a tadpole community by a multi-host fungal pathogen with limited evidence of recovery.

Emerging infectious diseases can cause host community disassembly, but the mechanisms driving the order of species declines and extirpations following a disease outbreak are unclear. We documented the community disassembly of a Neotropical tadpole community during a chytridiomycosis outbreak, triggered by the generalist fungal pathogen, Batrachochytrium dendrobatidis (Bd). Within the first 11 months of Bd arrival, tadpole density and occupancy rapidly declined. Species rarity, in terms of tadpole occupancy and adult relative abundance, did not predict the odds of tadpole occupancy declines. But species losses were taxonomically selective, with glassfrogs (Family: Centrolenidae) disappearing the fastest and tree frogs (Family: Hylidae) and dart-poison frogs (Family: Dendrobatidae) remaining the longest. We detected biotic homogenization of tadpole communities, with post-decline communities resembling one another more strongly than pre-decline communities. The entire tadpole community was extirpated within 22 months following Bd arrival, and we found limited signs of recovery within 10 years post-outbreak. Because of imperfect species detection inherent to sampling species-rich tropical communities and the difficulty of devising a single study design protocol to sample physically complex tropical habitats, we used simulations to provide recommendations for future surveys to adequately sample diverse Neotropical communities. Our unique data set on tadpole community composition before and after Bd arrival is a valuable baseline for assessing amphibian recovery. Our results are of direct relevance to conservation managers and community ecologists interested in understanding the timing, magnitude, and consequences of disease outbreaks as emerging infectious diseases spread globally.

Evidence for the persistence of food web structure after amphibian extirpation in a Neotropical stream.

Species losses are predicted to simplify food web structure, and disease-driven amphibian declines in Central America offer an opportunity to test this prediction. Assessment of insect community composition, combined with gut content analyses, was used to generate periphyton-insect food webs for a Panamanian stream, both pre- and post-amphibian decline. We then used network analysis to assess the effects of amphibian declines on food web structure. Although 48% of consumer taxa, including many insect taxa, were lost between pre- and post-amphibian decline sampling dates, connectance declined by less than 3%. We then quantified the resilience of food web structure by calculating the number of expected cascading extirpations from the loss of tadpoles. This analysis showed the expected effects of species loss on connectance and linkage density to be more than 60% and 40%, respectively, than were actually observed. Instead, new trophic linkages in the post-decline food web reorganized the food web topology, changing the identity of "hub" taxa, and consequently reducing the effects of amphibian declines on many food web attributes. Resilience of food web attributes was driven by a combination of changes in consumer diets, particularly those of insect predators, as well as the appearance of generalist insect consumers, suggesting that food web structure is maintained by factors independent of the original trophic linkages.

Testing the role of ecology and life history in structuring genetic variation across a landscape: a trait-based phylogeographic approach.

Hypotheses to explain phylogeographic structure traditionally invoke geographic features, but often fail to provide a general explanation for spatial patterns of genetic variation. Organisms' intrinsic characteristics might play more important roles than landscape features in determining phylogeographic structure. We developed a novel comparative approach to explore the role of ecological and life-history variables in determining spatial genetic variation and tested it on frog communities in Panama. We quantified spatial genetic variation within 31 anuran species based on mitochondrial DNA sequences, for which hierarchical approximate Bayesian computation analyses rejected simultaneous divergence over a common landscape. Regressing ecological variables, on genetic divergence allowed us to test the importance of individual variables revealing that body size, current landscape resistance, geographic range, biogeographic origin and reproductive mode were significant predictors of spatial genetic variation. Our results support the idea that phylogeographic structure represents the outcome of an interaction between organisms and their environment, and suggest a conceptual integration we refer to as trait-based phylogeography.

DNA barcoding applied to ex situ tropical amphibian conservation programme reveals cryptic diversity in captive populations.

Amphibians constitute a diverse yet still incompletely characterized clade of vertebrates, in which new species are still being discovered and described at a high rate. Amphibians are also increasingly endangered, due in part to disease-driven threats of extinctions. As an emergency response, conservationists have begun ex situ assurance colonies for priority species. The abundance of cryptic amphibian diversity, however, may cause problems for ex situ conservation. In this study we used a DNA barcoding approach to survey mitochondrial DNA (mtDNA) variation in captive populations of 10 species of Neotropical amphibians maintained in an ex situ assurance programme at El Valle Amphibian Conservation Center (EVACC) in the Republic of Panama. We combined these mtDNA sequences with genetic data from presumably conspecific wild populations sampled from across Panama, and applied genetic distance-based and character-based analyses to identify cryptic lineages. We found that three of ten species harboured substantial cryptic genetic diversity within EVACC, and an additional three species harboured cryptic diversity among wild populations, but not in captivity. Ex situ conservation efforts focused on amphibians are therefore vulnerable to an incomplete taxonomy leading to misidentification among cryptic species. DNA barcoding may therefore provide a simple, standardized protocol to identify cryptic diversity readily applicable to any amphibian community.

Epidemic disease decimates amphibian abundance, species diversity, and evolutionary history in the highlands of central Panama.

Amphibian populations around the world are experiencing unprecedented declines attributed to a chytrid fungal pathogen, Batrachochytrium dendrobatidis. Despite the severity of the crisis, quantitative analyses of the effects of the epidemic on amphibian abundance and diversity have been unavailable as a result of the lack of equivalent data collected before and following disease outbreak. We present a community-level assessment combining long-term field surveys and DNA barcode data describing changes in abundance and evolutionary diversity within the amphibian community of El Copé, Panama, following a disease epidemic and mass-mortality event. The epidemic reduced taxonomic, lineage, and phylogenetic diversity similarly. We discovered that 30 species were lost, including five undescribed species, representing 41% of total amphibian lineage diversity in El Copé. These extirpations represented 33% of the evolutionary history of amphibians within the community, and variation in the degree of population loss and decline among species was random with respect to the community phylogeny. Our approach provides a fast, economical, and informative analysis of loss in a community whether measured by species or phylogenetic diversity.

Batrachochytrium dendrobatidis infection patterns among Panamanian amphibian species, habitats and elevations during epizootic and enzootic stages.

The pathogenic fungus Batrachochytrium dendrobatidis (Bd) has caused declines of many amphibian populations, yet the full course of the epizootic has rarely been observed in wild populations. We determined effects of elevation, habitat, and aquatic index (AI) on prevalence of infection among Panamanian amphibians sampled along 2 elevational transects. Amphibian populations on the Santa Fé transect (SFT) had declined in 2002, while those on the El Copé transect (ECT) were healthy until September 2004. In 2004 we sampled Bd along both transects, surveying the SFT 2 yr after decline, and surveying the ECT 4 mo prior to the arrival of Bd, during the epizootic, and 2 mo later. Overall prevalence of Bd along the ECT increased from 0.0 (95% CI 0.00-0.0003) to 0.51 (95% CI 0.48-0.55) over a 3 mo period, accompanied by significant decreases in amphibian abundance and species richness in all habitats. Prevalence of infection on the ECT was highest along riparian transects and at higher elevations, but not among levels of AI. Prevalence of infection on the SFT was highest in pool transects, and at higher elevations, but not among levels of AI. Riparian amphibian abundance and species richness also declined at SFT following detection of Bd in 2002. Variation among species, microenvironmental conditions, and the length of coexistence with Bd may contribute to observed differences in prevalence of Bd and in population response.

Fine scale gene flow and individual movements among subpopulations of Centrolene prosoblepon (Anura: Centrolenidae).

Dispersal capabilities determine and maintain local gene flow, and this has implications for population persistence and/or recolonization following environmental perturbations (natural or anthropogenic), disease outbreaks, or other demographic collapses. To predict recolonization and understand dispersal capacity in a stream-breeding frog, we examined individual movement patterns and gene flow among four subpopulations of the Neotropical glassfrog, Centrolene prosoblepon, at a mid-elevation cloud forest site at El Copé, Panama. We measured male movement directly during a two year mark-recapture study, and indirectly with gene flow estimates from mitochondrial DNA sequences (mtDNA). Individuals of this species showed strong site fidelity: over two years, male frogs in all four headwater streams moved very little (mean = 2.33 m; mode = 0 m). Nine individuals changed streams within one or two years, moving 675-1,108 m. For those males moving more than 10 m, movement was biased upstream (p < 0.001). Using mtDNA ND1 gene sequences, we quantified gene flow within and among headwater streams at two spatial scales: among headwater streams within two adjacent watersheds (2.5 km2) and among streams within a longitudinal gradient covering 5.0 km2. We found high gene flow among headwater streams (phi(ST) = 0.007, p = 0.325) but gene flow was more limited across greater distances (phi(CT) = 0.322, p = 0.065), even within the same drainage network. Lowland populations of C. prosoblepon potentially act as an important source of colonists for upland populations in this watershed.

A new species of Celestus from west-central Panama, with consideration of the status of the genera of the Anguidae: Diploglossinae (Squamata)

A recently discovered new species of diploglossine lizard is described from west-central Panama. The distinctiveness of the nominal genera Celestus and Diploglossus is confirmed; the new form represents the southernmost record for the genus Celestus. A summary of selected characteristics and general distribution is presented for all recent species of diploglossines, including members of the Antillean genera Saurisia and Wetmorea and the South American genus Ophiodes. A systematic key to mainland members of the genus Celestus is provided. Rev. Biol. Trop. 56 (2): 845-859. Epub 2008 June 30.

Se describe una especie nueva de lagarto diploglossino recientemente descubierta en el centro-oeste de Panamá. Confirmamos la peculiaridad de los géneros nominales Celestus y Diploglossus. La nueva especie representa el registro más al sur del género Celestus. Se presenta un resumen de características selectas y de la distribución general de todas las especies recientes de diploglossinos, incluyendo a los miembros de los géneros antillanos Saurisia y Wetmorea y del género suramericano Ophiodes. Se provee una clave sistemática para los miembros de tierra firme del género Celestus.

Riding the wave: reconciling the roles of disease and climate change in amphibian declines.

We review the evidence for the role of climate change in triggering disease outbreaks of chytridiomycosis, an emerging infectious disease of amphibians. Both climatic anomalies and disease-related extirpations are recent phenomena, and effects of both are especially noticeable at high elevations in tropical areas, making it difficult to determine whether they are operating separately or synergistically. We compiled reports of amphibian declines from Lower Central America and Andean South America to create maps and statistical models to test our hypothesis of spatiotemporal spread of the pathogen Batrachochytrium dendrobatidis (Bd), and to update the elevational patterns of decline in frogs belonging to the genus Atelopus. We evaluated claims of climate change influencing the spread of Bd by including error into estimates of the relationship between air temperature and last year observed. Available data support the hypothesis of multiple introductions of this invasive pathogen into South America and subsequent spread along the primary Andean cordilleras. Additional analyses found no evidence to support the hypothesis that climate change has been driving outbreaks of amphibian chytridiomycosis, as has been posited in the climate-linked epidemic hypothesis. Future studies should increase retrospective surveys of museum specimens from throughout the Andes and should study the landscape genetics of Bd to map fine-scale patterns of geographic spread to identify transmission routes and processes.

Fine scale gene flow and individual movements among subpopulations of Centrolene prosoblepon (Anura: Centrolenidae)

Dispersal capabilities determine and maintain local gene flow, and this has implications for population persistence and/or recolonization following environmental perturbations (natural or anthropogenic), disease outbreaks, or other demographic collapses. To predict recolonization and understand dispersal capacity in a stream-breeding frog, we examined individual movement patterns and gene flow among four subpopulations of the Neotropical glassfrog, Centrolene prosoblepon, at a mid-elevation cloud forest site at El Copé, Panama. We measured male movement directly during a two year mark-recapture study, and indirectly with gene flow estimates from mitochondrial DNA sequences (mtDNA). Individuals of this species showed strong site fidelity: over two years, male frogs in all four headwater streams moved very little (mean = 2.33 m; mode = 0 m). Nine individuals changed streams within one or two years, moving 675-1,108 m. For those males moving more than 10 m, movement was biased upstream (p < 0.001). Using mtDNA ND1 gene sequences, we quantified gene flow within and among headwater streams at two spatial scales: among headwater streams within two adjacent watersheds (2.5 km2) and among streams within a longitudinal gradient covering 5.0 km2. We found high gene flow among headwater streams (phi(ST) = 0.007, p = 0.325) but gene flow was more limited across greater distances (phi(CT) = 0.322, p = 0.065), even within the same drainage network. Lowland populations of C. prosoblepon potentially act as an important source of colonists for upland populations in this watershed.

La capacidad de dispersión determina y mantiene el flujo genético local, y esto tiene implicaciones para la persistencia poblacional y/o la recolonización que sigue a perturbaciones ambientales. Examinamos patrones individuales de movimiento y flujo genético entre subpoblaciones de Centrolene prosoblepon (Anura: Centrolenidae) en un sitio de elevación media en El Copé, Panamá. Medimos directamente el movimiento de los machos durante un estudio de marcado-recaptura, e indirectamente con estimaciones de flujo genético a partir de secuencias de ADN mitocondrial (mtDNA). Los individuos mostraron fuerte fidelidad a su lugar: por más de dos años, las ranas macho de los cuatro arroyos al inicio del río se movieron muy poco (promedio = 2.33 m; moda = 0 m). Nueve individuos cambiaron de corriente de agua en uno o dos años, moviéndose 675-1 108 m. Usando la secuencia genética ND1 del ADN mitocondrial, medimos el flujo genético en dos escalas espaciales: entre arroyos que originan el río (2.5 km2) y entre arroyos con un gradiente longitudinal en 5.0 km2. Encontramos un flujo genético alto entre los arroyos al inicio del río (f = 0.007, p = 0.325 y otro más limitado en distancias mayores (f = 0.322, p = 0.065).

A new species of Celestus from west-central Panama, with consideration of the status of the genera of the Anguidae: Diploglossinae (Squamata).

A recently discovered new species of diploglossine lizard is described from west-central Panama. The distinctiveness of the nominal genera Celestus and Diploglossus is confirmed; the new form represents the southernmost record for the genus Celestus. A summary of selected characteristics and general distribution is presented for all recent species of diploglossines, including members of the Antillean genera Saurisia and Wetmorea and the South American genus Ophiodes. A systematic key to mainland members of the genus Celestus is provided.

Emerging infectious disease and the loss of biodiversity in a Neotropical amphibian community.

Pathogens rarely cause extinctions of host species, and there are few examples of a pathogen changing species richness and diversity of an ecological community by causing local extinctions across a wide range of species. We report the link between the rapid appearance of a pathogenic chytrid fungus Batrachochytrium dendrobatidis in an amphibian community at El Copé, Panama, and subsequent mass mortality and loss of amphibian biodiversity across eight families of frogs and salamanders. We describe an outbreak of chytridiomycosis in Panama and argue that this infectious disease has played an important role in amphibian population declines. The high virulence and large number of potential hosts of this emerging infectious disease threaten global amphibian diversity.

Morphological variation in the tropical anole, Anolis casildae (Squamata: Polychrotidae)

We describe morphological variation (scalation and coloration) observed among eight individuals of the Panamanian lizard species Anolis casildae. This variation was not observed in the holotype and aids in identification of this recently described species (originally described on the basis of a single, male specimen). This species occurs only in the Reserva Forestal Fortuna (ChiriquÝ Province) and the adjacent Bosque Protector Palo Seco (Bocas del Toro Province) in western Panama. Anolis casildae can be distinguished from all other Panamanian anole species via six features: (1) two enlarged superciliary scales (the first larger than the second); (2) an anterior nasal scale in contact with the rostral scale or separated from the rostral by one scale; (3) 6-8 sublabial scales to the center of the eye; (4) 3-4 scales between the supraobital semicircles; (5) unique coloration (4-6 oblique brown bands interspersed by blue-outlined yellow patches; dewlap is a dirty cream color with broad yellow scale rows irregularly interspersed with smaller emerald green scales) and (6) A. casildae occurs from 1,050 to 1,400 m in the Cordillera Central. We also compare our natural history observations of A. casildae to a similar large anole, A. frenatus, a species which we believe A. casildae to be closely related.