Aquatic Soundscape Recordings Reveal Diverse Vocalizations and Nocturnal Activity of an Endangered Frog

AbstractAutonomous sensors provide opportunities to observe organisms across spatial and temporal scales that humans cannot directly observe. By processing large data streams from autonomous sensors with deep learning methods, researchers can make novel and important natural history discoveries. In this study, we combine automated acoustic monitoring with deep learning models to observe breeding-associated activity in the endangered Sierra Nevada yellow-legged frog (Rana sierrae), a behavior that current surveys do not measure. By deploying inexpensive hydrophones and developing a deep learning model to recognize breeding-associated vocalizations, we discover three undocumented R. sierrae vocalization types and find an unexpected temporal pattern of nocturnal breeding-associated vocal activity. This study exemplifies how the combination of autonomous sensor data and deep learning can shed new light on species' natural history, especially during times or in locations where human observation is limited or impossible.

Host density has limited effects on pathogen invasion, disease-induced declines and within-host infection dynamics across a landscape of disease

1. Host density is hypothesized to be a major driver of variability in the responses and outcomes of wildlife populations following pathogen invasion. While the effects of host density on pathogen transmission have been extensively studied, these studies are dominated by theoretical analyses and small-scale experiments. This focus leads to an incomplete picture regarding how host density drives observed variability in disease outcomes in the field. 2. Here, we leveraged a dataset of hundreds of replicate amphibian populations that varied by orders of magnitude in host density. We used these data to test the effects of host density on three outcomes following the arrival of the amphibian-killing fungal pathogen Batrachochytrium dendrobatidis (Bd): the probability that Bd successfully invaded a host population and led to a pathogen outbreak, the magnitude of the host population-level decline following an outbreak and within-host infection dynamics that drive population-level outcomes in amphibian-pathogen systems. 3. Based on previous small-scale transmission experiments, we expected that populations with higher densities would be more likely to experience Bd outbreaks and would suffer larger proportional declines following outbreaks. To test these predictions, we developed and fitted a Hidden Markov Model that accounted for imperfectly observed disease outbreak states in the amphibian populations we surveyed. 4. Contrary to our predictions, we found minimal effects of host density on the probability of successful Bd invasion, the magnitude of population decline following Bd invasion and the dynamics of within-host infection intensity. Environmental conditions, such as summer temperature, winter severity and the presence of pathogen reservoirs, were more predictive of variability in disease outcomes. 5. Our results highlight the limitations of extrapolating findings from small-scale transmission experiments to observed disease trajectories in the field and provide strong evidence that variability in host density does not necessarily drive variability in host population responses following pathogen arrival. In an applied context, we show that feedbacks between host density and disease will not necessarily affect the success of reintroduction efforts in amphibian-Bd systems of conservation concern.

Effectiveness of antifungal treatments during chytridiomycosis epizootics in populations of an endangered frog

The recently-emerged amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd) has had an unprecedented impact on global amphibian populations, and highlights the urgent need to develop effective mitigation strategies. We conducted in-situ antifungal treatment experiments in wild populations of the endangered mountain yellow-legged frog during or immediately after Bd-caused mass die-off events. The objective of treatments was to reduce Bd infection intensity ("load") and in doing so alter frog-Bd dynamics and increase the probability of frog population persistence despite ongoing Bd infection. Experiments included treatment of early life stages (tadpoles and subadults) with the antifungal drug itraconazole, treatment of adults with itraconazole, and augmentation of the skin microbiome of subadults with Janthinobacterium lividum, a commensal bacterium with antifungal properties. All itraconazole treatments caused immediate reductions in Bd load, and produced longer-term effects that differed between life stages. In experiments focused on early life stages, Bd load was reduced in the 2 months immediately following treatment and was associated with increased survival of subadults. However, Bd load and frog survival returned to pre-treatment levels in less than 1 year, and treatment had no effect on population persistence. In adults, treatment reduced Bd load and increased frog survival over the entire 3-year post-treatment period, consistent with frogs having developed an effective adaptive immune response against Bd. Despite this protracted period of reduced impacts of Bd on adults, recruitment into the adult population was limited and the population eventually declined to near-extirpation. In the microbiome augmentation experiment, exposure of subadults to a solution of J. lividum increased concentrations of this potentially protective bacterium on frogs. However, concentrations declined to baseline levels within 1 month and did not have a protective effect against Bd infection. Collectively, these results indicate that our mitigation efforts were ineffective in causing long-term changes in frog-Bd dynamics and increasing population persistence, due largely to the inability of early life stages to mount an effective immune response against Bd. This results in repeated recruitment failure and a low probability of population persistence in the face of ongoing Bd infection.

Divergent regional evolutionary histories of a devastating global amphibian pathogen

Emerging infectious diseases are a pressing threat to global biological diversity. Increased incidence and severity of novel pathogens underscores the need for methodological advances to understand pathogen emergence and spread. Here, we use genetic epidemiology to test, and challenge, key hypotheses about a devastating zoonotic disease impacting amphibians globally. Using an amplicon-based sequencing method and non-invasive samples we retrospectively explore the history of the fungal pathogen Batrachochytrium dendrobatidis (Bd) in two emblematic amphibian systems: the Sierra Nevada of California and Central Panama. The hypothesis in both regions is the hypervirulent Global Panzootic Lineage of Bd (BdGPL) was recently introduced and spread rapidly in a wave-like pattern. Our data challenge this hypothesis by demonstrating similar epizootic signatures can have radically different underlying evolutionary histories. In Central Panama, our genetic data confirm a recent and rapid pathogen spread. However, BdGPL in the Sierra Nevada has remarkable spatial structuring, high genetic diversity and a relatively older history inferred from time-dated phylogenies. Thus, this deadly pathogen lineage may have a longer history in some regions than assumed, providing insights into its origin and spread. Overall, our results highlight the importance of integrating observed wildlife die-offs with genetic data to more accurately reconstruct pathogen outbreaks.

Drosophila melanogaster and worker honeybees (Apis mellifera) do not require olfaction to be susceptible to honeybee queen mandibular pheromone

Eusociality is characterised by the reproductive division of labour; a dominant female (queen) or females are responsible for the majority of reproduction, and subordinate females are reproductively constrained. Reproductive constraint can be due to behavioural aggression and/or chemical cues, so-called queen pheromones, produced by the dominant females. In the honeybee, Apis mellifera, this repressive queen pheromone is queen mandibular pheromone (QMP). The mechanism by which honeybee workers are susceptible to QMP is not yet completely understood, however it is thought to be through olfaction via the antennae and/or gustation via trophallaxis. We have investigated whether olfaction is key to sensing of QMP, using both Drosophila melanogaster- a tractable non-eusocial insect which is also reproductively repressed by QMP- and the target species, A. mellifera worker honeybees. D. melanogaster are still capable of sensing and responding to QMP without their antenna and maxillary palps, and therefore without olfactory receptors. When worker honeybees were exposed to QMP but unable to physically interact with it, therefore required to use olfaction, they were similarly not reproductively repressed. Combined, these findings support either a non-olfactory based mechanism for the repression of reproduction via QMP, or redundancy via non-olfactory mechanisms in both D. melanogaster and A. mellifera. This study furthers our understanding of how species are susceptible to QMP, and provides insight into the mechanisms governing QMP responsiveness in these diverse species.

Stepping into the past to conserve the future: Archived skin swabs from extant and extirpated populations inform genetic management of an endangered amphibian

Moving animals on a landscape through translocations and reintroductions is an important management tool used in the recovery of endangered species, particularly for the maintenance of population genetic diversity and structure. Management of imperiled amphibian species rely heavily on translocations and reintroductions, especially for species that have been brought to the brink of extinction by habitat loss, introduced species, and disease. One striking example of amphibian declines and associated management efforts is in California's Sequoia and Kings Canyon National Parks with the mountain yellow-legged frog species complex (Rana sierrae/muscosa). Mountain yellow-legged frogs have been extirpated from more than 93% of their historic range, and limited knowledge of their population genetics has made long-term conservation planning difficult. To address this, we used 598 archived skin swabs from both extant and extirpated populations across 48 lake basins to generate a robust Illumina-based nuclear amplicon data set. We found that samples grouped into three main genetic clusters, concordant with watershed boundaries. We also found evidence for historical gene flow across watershed boundaries with a north-to-south axis of migration. Finally, our results indicate that genetic diversity is not significantly different between populations with different disease histories. Our study offers specific management recommendations for imperiled mountain yellow-legged frogs and, more broadly, provides a population genetic framework for leveraging minimally invasive samples for the conservation of threatened species.

Cryptic diversity of a widespread global pathogen reveals expanded threats to amphibian conservation

Biodiversity loss is one major outcome of human-mediated ecosystem disturbance. One way that humans have triggered wildlife declines is by transporting disease-causing agents to remote areas of the world. Amphibians have been hit particularly hard by disease due in part to a globally distributed pathogenic chytrid fungus (Batrachochytrium dendrobatidis [Bd]). Prior research has revealed important insights into the biology and distribution of Bd; however, there are still many outstanding questions in this system. Although we know that there are multiple divergent lineages of Bd that differ in pathogenicity, we know little about how these lineages are distributed around the world and where lineages may be coming into contact. Here, we implement a custom genotyping method for a global set of Bd samples. This method is optimized to amplify and sequence degraded DNA from noninvasive skin swab samples. We describe a divergent lineage of Bd, which we call BdASIA3, that appears to be widespread in Southeast Asia. This lineage co-occurs with the global panzootic lineage (BdGPL) in multiple localities. Additionally, we shed light on the global distribution of BdGPL and highlight the expanded range of another lineage, BdCAPE. Finally, we argue that more monitoring needs to take place where Bd lineages are coming into contact and where we know little about Bd lineage diversity. Monitoring need not use expensive or difficult field techniques but can use archived swab samples to further explore the history-and predict the future impacts-of this devastating pathogen.

Reduced skin bacterial diversity correlates with increased pathogen infection intensity in an endangered amphibian host

The fungal pathogen Batrachochytrium dendrobatidis (Bd) infects the skin of amphibians and has caused severe declines and extinctions of amphibians globally. In this study, we investigate the interaction between Bd and the bacterial skin microbiome of the endangered Sierra Nevada yellow-legged frog, Rana sierrae, using both culture-dependent and culture-independent methods. Samples were collected from two populations of R. sierrae that likely underwent Bd epizootics in the past, but that continue to persist with Bd in an enzootic disease state, and we address the hypothesis that such "persistent" populations are aided by mutualistic skin microbes. Our 16S rRNA metabarcoding data reveal that the skin microbiome of highly infected juvenile frogs is characterized by significantly reduced species richness and evenness, and by strikingly lower variation between individuals, compared to juveniles and adults with lower infection levels. Over 90% of DNA sequences from the skin microbiome of highly infected frogs were derived from bacteria in a single order, Burkholderiales, compared to just 54% in frogs with lower infection levels. In a culture-dependent Bd inhibition assay, the bacterial metabolites we evaluated all inhibited the growth of Bd. Together, these results illustrate the disruptive effects of Bd infection on host skin microbial community structure and dynamics, and suggest possible avenues for the development of anti-Bd probiotic treatments.

Epidemic and endemic pathogen dynamics correspond to distinct host population microbiomes at a landscape scale

Infectious diseases have serious impacts on human and wildlife populations, but the effects of a disease can vary, even among individuals or populations of the same host species. Identifying the reasons for this variation is key to understanding disease dynamics and mitigating infectious disease impacts, but disentangling cause and correlation during natural outbreaks is extremely challenging. This study aims to understand associations between symbiotic bacterial communities and an infectious disease, and examines multiple host populations before or after pathogen invasion to infer likely causal links. The results show that symbiotic bacteria are linked to fundamentally different outcomes of pathogen infection: host-pathogen coexistence (endemic infection) or host population extirpation (epidemic infection). Diversity and composition of skin-associated bacteria differed between populations of the frog, Rana sierrae, that coexist with or were extirpated by the fungal pathogen, Batrachochytrium dendrobatidis (Bd). Data from multiple populations sampled before or after pathogen invasion were used to infer cause and effect in the relationship between the fungal pathogen and symbiotic bacteria. Among host populations, variation in the composition of the skin microbiome was most strongly predicted by pathogen infection severity, even in analyses where the outcome of infection did not vary. This result suggests that pathogen infection shapes variation in the skin microbiome across host populations that coexist with or are driven to extirpation by the pathogen. By contrast, microbiome richness was largely unaffected by pathogen infection intensity, but was strongly predicted by geographical region of the host population, indicating the importance of environmental or host genetic factors in shaping microbiome richness. Thus, while both richness and composition of the microbiome differed between endemic and epidemic host populations, the underlying causes are most likely different: pathogen infection appears to shape microbiome composition, while microbiome richness was less sensitive to pathogen-induced disturbance. Because higher richness was correlated with host persistence in the presence of Bd, and richness appeared relatively stable to Bd infection, microbiome richness may contribute to disease resistance, although the latter remains to be directly tested.

Resistance, tolerance and environmental transmission dynamics determine host extinction risk in a load-dependent amphibian disease

While disease-induced extinction is generally considered rare, a number of recently emerging infectious diseases with load-dependent pathology have led to extinction in wildlife populations. Transmission is a critical factor affecting disease-induced extinction, but the relative importance of transmission compared to load-dependent host resistance and tolerance is currently unknown. Using a combination of models and experiments on an amphibian species suffering extirpations from the fungal pathogen Batrachochytrium dendrobatidis (Bd), we show that while transmission from an environmental Bd reservoir increased the ability of Bd to invade an amphibian population and the extinction risk of that population, Bd-induced extinction dynamics were far more sensitive to host resistance and tolerance than to Bd transmission. We demonstrate that this is a general result for load-dependent pathogens, where non-linear resistance and tolerance functions can interact such that small changes in these functions lead to drastic changes in extinction dynamics.

Increased accuracy of species lists developed for alpine lakes using morphology and cytochrome oxidase I for identification of specimens

The first step in many community ecology studies is to produce a species list from a sample of individuals. Community ecologists now have two viable ways of producing a species list: morphological and barcode identification. In this study, we compared the taxonomic resolution gained by a combined use of both methods and tested whether a change in taxonomic resolution significantly impacted richness estimates for benthic macroinvertebrates sampled from ten lakes in Sequoia National Park, USA. Across all lakes, 77 unique taxa were identified and 42% (32) were reliably identified to species using both barcode and morphological identification. Of the 32 identified to species, 63% (20) were identified solely by comparing the barcode sequence from cytochrome oxidase I to the Barcode of Life reference library. The increased resolution using a combined identification approach compared to identifications based solely on morphology resulted in a significant increase in estimated richness within a lake at the order, family, genus and species levels of taxonomy (P < 0.05). Additionally, young or damaged individuals that could not be identified using morphology were identified using their COI sequences to the genus or species level on average 75% of the time. Our results demonstrate that a combined identification approach improves accuracy of benthic macroinvertebrate species lists in alpine lakes and subsequent estimates of richness. We encourage the use of barcodes for identification purposes and specifically when morphology is insufficient, as in the case of damaged and early life stage specimens of benthic macroinvertebrates.

Temporal and spatial variation of atmospherically deposited organic contaminants at high elevation in Yosemite National Park, California, USA

Contaminants used at low elevation, such as pesticides on crops, can be transported tens of kilometers and deposited in adjacent mountains in many parts of the world. Atmospherically deposited organic contaminants in the Sierra Nevada Mountains of California, USA, have exceeded some thresholds of concern, but the spatial and temporal distributions of contaminants in the mountains are not well known. The authors sampled shallow-water sediment and tadpoles (Pseudacris sierra) for pesticides, polycyclic aromatic hydrocarbons (PAHs), and polychlorinated biphenyls in four high-elevation sites in Yosemite National Park in the central Sierra Nevada twice during the summers of 2006, 2007, and 2008. Both historic- and current-use pesticides showed a striking pattern of lower concentrations in both sediment and tadpoles in Yosemite than was observed previously in Sequoia-Kings Canyon National Parks in the southern Sierra Nevada. By contrast, PAH concentrations in sediment were generally greater in Yosemite than in Sequoia-Kings Canyon. The authors suggest that pesticide concentrations tend to be greater in Sequoia-Kings Canyon because of a longer air flow path over agricultural lands for this park along with greater pesticide use near this park. Concentrations for DDT-related compounds in some sediment samples exceeded guidelines or critical thresholds in both parks. A general pattern of difference between Yosemite and Sequoia-Kings Canyon was not evident for total tadpole cholinesterase activity, an indicator of harmful exposure to organophosphorus and carbamate pesticides. Variability of chemical concentrations among sites, between sampling periods within each year, and among years, contributed significantly to total variation, although the relative contributions differed between sediment and tadpoles.

Pesticide distributions and population declines of California, USA, alpine frogs, Rana muscosa and Rana sierrae

Atmospherically deposited pesticides from the intensively cultivated Central Valley of California, USA, have been implicated as a cause for population declines of several amphibian species, with the strongest evidence for the frogs Rana muscosa and Rana sierrae at high elevation in the Sierra Nevada mountains. Previous studies on these species have relied on correlations between frog population status and either a metric for amount of upwind pesticide use or limited measurements of pesticide concentrations in the field. The present study tested the hypothesis that pesticide concentrations are negatively correlated with frog population status (i.e., fraction of suitable water bodies occupied within 2 km of a site) by measuring pesticide concentrations in multiple media twice at 28 sites at high elevation in the southern Sierra Nevada. Media represented were air, sediment, and Pseudacris sierra tadpoles. Total cholinesterase (ChE), which has been used as an indicator for organophosphorus and carbamate pesticide exposure, was also measured in P. sierra tadpoles. Results do not support the pesticide-site occupancy hypothesis. Among 46 pesticide compounds analyzed, nine were detected with ≥ 30% frequency, representing both historically and currently used pesticides. In stepwise regressions with a chemical metric and linear distance from the Central Valley as predictor variables, no negative association was found between frog population status and the concentration of any pesticide or tadpole ChE activity level. By contrast, frog population status showed a strong positive relationship with linear distance from the Valley, a pattern that is consistent with a general west-to-east spread across central California of the amphibian disease chytridiomycosis observed by other researchers.

Nonnative trout impact an alpine-nesting bird by altering aquatic-insect subsidies

Adjacent food webs may be linked by cross-boundary subsidies: more-productive donor systems can subsidize consumers in less-productive neighboring recipient systems. Introduced species are known to have direct effects on organisms within invaded communities. However, few studies have addressed the indirect effects of nonnative species in donor systems on organisms in recipient systems. We studied the direct role of introduced trout in altering a lake-derived resource subsidy and their indirect effects in altering a passerine bird's response to that subsidy. We compared the abundance of aquatic insects and foraging Gray-crowned Rosy-Finches (Leucosticte tephrocotis dawsoni, "Rosy-Finch") at fish-containing vs. fishless lakes in the Sierra Nevada Mountains of California (USA). Introduced trout outcompeted Rosy-Finches for emerging aquatic insects (i.e., mayflies). Fish-containing lakes had 98% fewer mayflies than did fishless lakes. In lakes without fish, Rosy-Finches showed an aggregative response to emerging aquatic insects with 5.9 times more Rosy-Finches at fishless lakes than at fish-containing lakes. Therefore, the introduction of nonnative fish into the donor system reduced both the magnitude of the resource subsidy and the strength of cross-boundary trophic interactions. Importantly, the timing of the subsidy occurs when Rosy-Finches feed their young. If Rosy-Finches rely on aquatic-insect subsidies to fledge their young, reductions in the subsidy by introduced trout may have decreased Rosy-Finch abundances from historic levels. We recommend that terrestrial recipients of aquatic subsidies be included in conservation and restoration plans for ecosystems with alpine lakes.

Allee effect limits colonization success of sexually reproducing zooplankton

Understanding the dynamics of populations at low density and the role of Allee effects is a priority due to concern about the decline of rare species and interest in colonization/invasion dynamics. Despite well-developed theory and observational support, experimental examinations of the Allee effect in natural systems are rare, partly because of logistical difficulties associated with experiments at low population density. We took advantage of fish introduction and removal in alpine lakes to experimentally test for the Allee effect at the whole-ecosystem scale. The large copepod Hesperodiaptomus shoshone is often extirpated from the water column by fish and sometimes fails to recover following fish disappearance, despite the presence of a long-lived egg bank. Population growth rate of this dioecious species may be limited by mate encounter rate, such that below some critical density a colonizing population will fail to establish. We conducted a multi-lake experiment in which H. shoshone was stocked at densities that bracketed our hypothesized critical density of 0.5-5 copoepods/m3. Successful recovery by the copepod was observed only in the lake with the highest initial density (3 copepods/m3). Copepods stocked into small cages at 3000 copepods/m3 survived and reproduced at rates comparable to natural populations, confirming that the lakes were suitable habitat for this species. In support of mate limitation as the mechanism underlying recovery failure, we found a significant positive relationship between mating success and density across experimental and natural H. shoshone populations. Furthermore, a mesocosm experiment provided evidence of increased per capita population growth rate with increasing population density in another diaptomid species, Skistodiaptomus pallidus. Together, these lines of evidence support the importance of the Allee effect to population recovery of H. shoshone in the Sierra Nevada, and to diaptomid copepods in general.

Recovery after local extinction: factors affecting re-establishment of alpine lake zooplankton

The introduction of fishes into naturally fishless mountain lakes often results in the extirpation of large-bodied zooplankton species. The ability to predict whether or not particular species will recover following fish removal is critically important for the design and implementation of lake restoration efforts but is currently not possible because of a lack of information on what factors affect recovery. The objective of this study was to identify the factors influencing recovery probability in two large-bodied zooplankton species following fish removal. We predicted that (1) Daphnia melanica would have a higher probability of recovery than Hesperodiaptomus shoshone due to differences in reproductive mode (D. melanica is parthenogenetic, H. shoshone is obligately sexual), (2) recovery probability would be a decreasing function of fish residence time due to the negative relationship between fish residence time and size of the egg bank, and (3) recovery probability would be an increasing function of lake depth as a consequence of a positive relationship between lake depth and egg bank size. To test these predictions, we sampled contemporary zooplankton populations and collected paleolimnological data from 44 naturally fishless lakes that were stocked with trout for varying lengths of time before reverting to a fishless condition. D. melanica had a significantly higher probability of recovery than did H. shoshone (0.82 vs. 0.54, respectively). The probability of recovery for H. shoshone was also significantly influenced by lake depth, fish residence time, and elevation, but only elevation influenced the probability of recovery in D. melanica. These results are consistent with between-species differences in reproductive mode combined with the much greater longevity of diapausing eggs in D. melanica than in H. shoshone. Our data also suggest that H. shoshone will often fail to recover in lakes with fish residence times exceeding 50 years.

Population genetics of the frog-killing fungus Batrachochytrium dendrobatidis

Global amphibian decline by chytridiomycosis is a major environmental disaster that has been attributed to either recent fungal spread or environmental change that promotes disease. Here, we present a population genetic comparison of Batrachochytrium dendrobatidis isolates from an intensively studied region of frog decline, the Sierra Nevada of California. In support of a novel pathogen, we find low diversity, no amphibian-host specificity, little correlation between fungal genotype and geography, local frog extirpation by a single fungal genotype, and evidence of human-assisted fungus migration. In support of endemism, at a local scale, we find some diverse, recombining populations. Therefore neither epidemic spread nor endemism alone explains this particular amphibian decline. Recombination raises the possibility of resistant sporangia and a mechanism for rapid spread as well as persistence that could greatly complicate global control of the pathogen.

Multiple stressors and amphibian declines: dual impacts of pesticides and fish on yellow-legged frogs

More than 40% of Earth's 5700+ amphibian species have undergone recent declines. Despite the likely involvement of multiple factors in driving these declines, most studies continue to focus on single stressors. In California (USA), separate studies have implicated either introduced fish or pesticides as causal agents. To date, however, no study has simultaneously evaluated the respective roles of these two potential stressors nor attempted to assess their relative importance, information critical for the development of effective conservation efforts and environmental policies. We examined the role and relative effect of fish and pesticides on the mountain yellow-legged frog (Rana muscosa) using unusually detailed data sets for a large portion of R. muscosa's historic range in California's Sierra Nevada. Habitat characteristics and presence/absence of R. muscosa and fish were quantified at each of 6831 sites during field surveys. Pesticide use upwind of each site was calculated from pesticide application records and predominant wind directions. Using generalized additive models, we found that, after accounting for habitat effects, the probability of R. muscosa presence was significantly reduced by both fish and pesticides, with the landscape-scale effect of pesticides much stronger than that of fish. The degree to which a site was sheltered from the predominant wind (and associated pesticides) was also a significant predictor of R. muscosa presence. Taken together, these results represent the strongest evidence to date that windborne pesticides are contributing to amphibian declines in pristine locations. Our results suggest that amphibian declines may have complex multi-factorial causes, and caution that single-factor studies that demonstrate the importance of one factor should not be used as evidence against the importance of other factors.

Emerging infectious disease as a proximate cause of amphibian mass mortality

A newly discovered infectious disease of amphibians, chytridiomycosis, caused by the fungal pathogen Batrachochytrium dendrobatidis, is implicated in population declines and possible extinctions throughout the world. The purpose of our study was to examine the effects of B. dendrobatidis on the mountain yellow-legged frog (Rana muscosa) in the Sierra Nevada of California (USA). We (1) quantified the prevalence and incidence of B. dendrobatidis through repeat surveys of several hundred R. muscosa populations in the southern Sierra Nevada; (2) described the population-level effects of B. dendrobatidis on R. muscosa population abundance; and (3) compared the mortality rates of infected and uninfected R. muscosa individuals from pre- through post-metamorphosis using both laboratory and field experiments. Mouthpart inspections conducted in 144 and 132 R. muscosa populations in 2003 and 2004, respectively, indicated that 19% of R. muscosa populations in both years showed indications of chytridiomycosis. Sixteen percent of populations that were uninfected in 2003 became infected by 2004. Rana muscosa population sizes were reduced by an average of 88% following B. dendrobatidis outbreaks at six sites, but at seven B. dendrobatidis-negative sites, R. muscosa population sizes increased by an average of 45% over the same time period. In the laboratory, all infected R. muscosa developed fatal chytridiomycosis after metamorphosis, while all uninfected individuals remained healthy. In the field experiment in which R. muscosa tadpoles were caged at infected and uninfected sites, 96% of the individuals that metamorphosed at infected sites died vs. 5% at the uninfected sites. These studies indicate that chytridiomycosis causes high mortality in post-metamorphic R. muscosa, that this emerging disease is the proximate cause of numerous observed R. muscosa population declines, and that the disease threatens this species with extirpation at numerous sites in California's Sierra Nevada.