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.

The amphibian microbiome exhibits poor resilience following pathogen-induced disturbance.

Infectious pathogens can disrupt the microbiome in addition to directly affecting the host. Impacts of disease may be dependent on the ability of the microbiome to recover from such disturbance, yet remarkably little is known about microbiome recovery after disease, particularly in nonhuman animals. We assessed the resilience of the amphibian skin microbial community after disturbance by the pathogen, Batrachochytrium dendrobatidis (Bd). Skin microbial communities of laboratory-reared mountain yellow-legged frogs were tracked through three experimental phases: prior to Bd infection, after Bd infection (disturbance), and after clearing Bd infection (recovery period). Bd infection disturbed microbiome composition and altered the relative abundances of several dominant bacterial taxa. After Bd infection, frogs were treated with an antifungal drug that cleared Bd infection, but this did not lead to recovery of microbiome composition (measured as Unifrac distance) or relative abundances of dominant bacterial groups. These results indicate that Bd infection can lead to an alternate stable state in the microbiome of sensitive amphibians, or that microbiome recovery is extremely slow-in either case resilience is low. Furthermore, antifungal treatment and clearance of Bd infection had the additional effect of reducing microbial community variability, which we hypothesize results from similarity across frogs in the taxa that colonize community vacancies resulting from the removal of Bd. Our results indicate that the skin microbiota of mountain yellow-legged frogs has low resilience following Bd-induced disturbance and is further altered by the process of clearing Bd infection, which may have implications for the conservation of this endangered amphibian.

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.

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.

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.

Large-scale recovery of an endangered amphibian despite ongoing exposure to multiple stressors.

Amphibians are one of the most threatened animal groups, with 32% of species at risk for extinction. Given this imperiled status, is the disappearance of a large fraction of the Earth's amphibians inevitable, or are some declining species more resilient than is generally assumed? We address this question in a species that is emblematic of many declining amphibians, the endangered Sierra Nevada yellow-legged frog (Rana sierrae). Based on >7,000 frog surveys conducted across Yosemite National Park over a 20-y period, we show that, after decades of decline and despite ongoing exposure to multiple stressors, including introduced fish, the recently emerged disease chytridiomycosis, and pesticides, R. sierrae abundance increased sevenfold during the study and at a rate of 11% per year. These increases occurred in hundreds of populations throughout Yosemite, providing a rare example of amphibian recovery at an ecologically relevant spatial scale. Results from a laboratory experiment indicate that these increases may be in part because of reduced frog susceptibility to chytridiomycosis. The disappearance of nonnative fish from numerous water bodies after cessation of stocking also contributed to the recovery. The large-scale increases in R. sierrae abundance that we document suggest that, when habitats are relatively intact and stressors are reduced in their importance by active management or species' adaptive responses, declines of some amphibians may be partially reversible, at least at a regional scale. Other studies conducted over similarly large temporal and spatial scales are critically needed to provide insight and generality about the reversibility of amphibian declines at a global scale.

Detecting the influence of rare stressors on rare species in Yosemite National Park using a novel stratified permutation test.

Statistical models often use observational data to predict phenomena; however, interpreting model terms to understand their influence can be problematic. This issue poses a challenge in species conservation where setting priorities requires estimating influences of potential stressors using observational data. We present a novel approach for inferring influence of a rare stressor on a rare species by blending predictive models with nonparametric permutation tests. We illustrate the approach with two case studies involving rare amphibians in Yosemite National Park, USA. The endangered frog, Rana sierrae, is known to be negatively impacted by non-native fish, while the threatened toad, Anaxyrus canorus, is potentially affected by packstock. Both stressors and amphibians are rare, occurring in ~10% of potential habitat patches. We first predict amphibian occupancy with a statistical model that includes all predictors but the stressor to stratify potential habitat by predicted suitability. A stratified permutation test then evaluates the association between stressor and amphibian, all else equal. Our approach confirms the known negative relationship between fish and R. sierrae, but finds no evidence of a negative relationship between current packstock use and A. canorus breeding. Our statistical approach has potential broad application for deriving understanding (not just prediction) from observational data.

A network extension of species occupancy models in a patchy environment applied to the Yosemite Toad (Anaxyrus canorus).

A central challenge of conservation biology is using limited data to predict rare species occurrence and identify conservation areas that play a disproportionate role in regional persistence. Where species occupy discrete patches in a landscape, such predictions require data about environmental quality of individual patches and the connectivity among high quality patches. We present a novel extension to species occupancy modeling that blends traditional predictions of individual patch environmental quality with network analysis to estimate connectivity characteristics using limited survey data. We demonstrate this approach using environmental and geospatial attributes to predict observed occupancy patterns of the Yosemite toad (Anaxyrus (= Bufo) canorus) across >2,500 meadows in Yosemite National Park (USA). A. canorus, a Federal Proposed Species, breeds in shallow water associated with meadows. Our generalized linear model (GLM) accurately predicted ~84% of true presence-absence data on a subset of data withheld for testing. The predicted environmental quality of each meadow was iteratively 'boosted' by the quality of neighbors within dispersal distance. We used this park-wide meadow connectivity network to estimate the relative influence of an individual Meadow's 'environmental quality' versus its 'network quality' to predict: a) clusters of high quality breeding meadows potentially linked by dispersal, b) breeding meadows with high environmental quality that are isolated from other such meadows, c) breeding meadows with lower environmental quality where long-term persistence may critically depend on the network neighborhood, and d) breeding meadows with the biggest impact on park-wide breeding patterns. Combined with targeted data on dispersal, genetics, disease, and other potential stressors, these results can guide designation of core conservation areas for A. canorus in Yosemite National Park.

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.

Dynamics of an emerging disease drive large-scale amphibian population extinctions.

Epidemiological theory generally suggests that pathogens will not cause host extinctions because the pathogen should fade out when the host population is driven below some threshold density. An emerging infectious disease, chytridiomycosis, caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd) is directly linked to the recent extinction or serious decline of hundreds of amphibian species. Despite continued spread of this pathogen into uninfected areas, the dynamics of the host-pathogen interaction remain unknown. We use fine-scale spatiotemporal data to describe (i) the invasion and spread of Bd through three lake basins, each containing multiple populations of the mountain yellow-legged frog, and (ii) the accompanying host-pathogen dynamics. Despite intensive sampling, Bd was not detected on frogs in study basins until just before epidemics began. Following Bd arrival in a basin, the disease spread to neighboring populations at approximately 700 m/yr in a wave-like pattern until all populations were infected. Within a population, infection prevalence rapidly reached 100% and infection intensity on individual frogs increased in parallel. Frog mass mortality began only when infection intensity reached a critical threshold and repeatedly led to extinction of populations. Our results indicate that the high growth rate and virulence of Bd allow the near-simultaneous infection and buildup of high infection intensities in all host individuals; subsequent host population crashes therefore occur before Bd is limited by density-dependent factors. Preventing infection intensities in host populations from reaching this threshold could provide an effective strategy to avoid the extinction of susceptible amphibian species in the wild.

Enzootic and epizootic dynamics of the chytrid fungal pathogen of amphibians.

Chytridiomycosis, the disease caused by the chytrid fungus, Batrachochytrium dendrobatidis (Bd), has contributed to amphibian population declines and extinctions worldwide. The impact of this pathogen, however, varies markedly among amphibian species and populations. Following invasion into some areas of California's Sierra Nevada, Bd leads to rapid declines and local extinctions of frog populations (Rana muscosa, R. sierrae). In other areas, infected populations of the same frog species have declined but persisted at low host densities for many years. We present results of a 5-year study showing that infected adult frogs in persistent populations have low fungal loads, are surviving between years, and frequently lose and regain the infection. Here we put forward the hypothesis that fungal load dynamics can explain the different population-level outcomes of Bd observed in different areas of the Sierra Nevada and possibly throughout the world. We develop a model that incorporates the biological details of the Bd-host interaction. Importantly, model results suggest that host persistence versus extinction does not require differences in host susceptibility, pathogen virulence, or environmental conditions, and may be just epidemic and endemic population dynamics of the same host-pathogen system. The different disease outcomes seen in natural populations may result solely from density-dependent host-pathogen dynamics. The model also shows that persistence of Bd is enhanced by the long-lived tadpole stage that characterize these two frog species, and by nonhost Bd reservoirs.

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.