Predicting breeding habitat for amphibians: a spatiotemporal analysis across Yellowstone National Park.

The ability to predict amphibian breeding across landscapes is important for informing land management decisions and helping biologists better understand and remediate factors contributing to declines in amphibian populations. We built geospatial models of likely breeding habitats for each of four amphibian species that breed in Yellowstone National Park (YNP). We used field data collected in 2000-2002 from 497 sites among 16 basins and predictor variables from geospatial models produced from remotely sensed data (e.g., digital elevation model, complex topographic index, landform data, wetland probability, and vegetative cover). Except for 31 sites in one basin that were surveyed in both 2000 and 2002, all sites were surveyed once. We used polytomous regression to build statistical models for each species of amphibian from (1) field survey site data only, (2) field data combined with data from geospatial models, and (3) data from geospatial models only. Based on measures of receiver operating characteristic (ROC) scores, models of the second type best explained likely breeding habitat because they contained the most information (ROC values ranged from 0.70 to 0.88). However, models of the third type could be applied to the entire YNP landscape and produced maps that could be verified with reserve field data. Accuracy rates for models built for single years were highly variable, ranging from 0.30 to 0.78. Accuracy rates for models built with data combined from multiple years were higher and less variable, ranging from 0.60 to 0.80. Combining results from the geospatial multiyear models yielded maps of "core" breeding areas (areas with high probability values for all three years) surrounded by areas that scored high for only one or two years, providing an estimate of variability among years. Such information can highlight landscape options for amphibian conservation. For example, our models identify alternative areas that could be protected for each species, including 6828-10 764 ha for tiger salamanders, 971-3017 ha for western toads, 4732-16 696 ha for boreal chorus frogs, and 4940-19 690 ha for Columbia spotted frogs.

Distribution and pathogenicity of Batrachochytrium dendrobatidis in boreal toads from the Grand Teton area of western Wyoming.

The pathogen Batrachochytrium dendrobatidis (Bd), which causes the skin disease chytridiomycosis, has been linked to amphibian population declines and extinctions worldwide. Bd has been implicated in recent declines of boreal toads, Bufo boreas boreas, in Colorado but populations of boreal toads in western Wyoming have high prevalence of Bd without suffering catastrophic mortality. In a field and laboratory study, we investigated the prevalence of Bd in boreal toads from the Grand Teton ecosystem (GRTE) in Wyoming and tested the pathogenicity of Bd to these toads in several environments. The pathogen was present in breeding adults at all 10 sites sampled, with a mean prevalence of 67%. In an experiment with juvenile toadlets housed individually in wet environments, 10(6) zoospores of Bd isolated from GRTE caused lethal disease in all Wyoming and Colorado animals within 35 days. Survival time was longer in toadlets from Wyoming than Colorado and in toadlets spending more time in dry sites. In a second trial involving Colorado toadlets exposed to 35% fewer Bd zoospores, infection peaked and subsided over 68 days with no lethal chytridiomycosis in any treatment. However, compared with drier aquaria with dry refuges, Bd infection intensity was 41% higher in more humid aquaria and 81% higher without dry refuges available. Our findings suggest that although widely infected in nature, Wyoming toads may escape chytridiomycosis due to a slight advantage in innate resistance or because their native habitat hinders Bd growth or provides more opportunities to reduce pathogen loads behaviorally than in Colorado.

Saprolegniaceae identified on amphibian eggs throughout the Pacific Northwest, USA, by internal transcribed spacer sequences and phylogenetic analysis.

We assessed the diversity and phylogeny of Saprolegniaceae on amphibian eggs from the Pacific Northwest, with particular focus on Saprolegnia ferax, a species implicated in high egg mortality. We identified isolates from eggs of six amphibians with the internal transcribed spacer (ITS) and 5.8S gene regions and BLAST of the GenBank database. We identified 68 sequences as Saprolegniaceae and 43 sequences as true fungi from at least nine genera. Our phylogenetic analysis of the Saprolegniaceae included isolates within the genera Saprolegnia, Achlya and Leptolegnia. Our phylogeny grouped S. semihypogyna with Achlya rather than with the Saprolegnia reference sequences. We found only one isolate that grouped closely with S. ferax, and this came from a hatchery-raised salmon (Idaho) that we sampled opportunistically. We had representatives of 7-12 species and three genera of Saprolegniaceae on our amphibian eggs. Further work on the ecological roles of different species of Saprolegniaceae is needed to clarify their potential importance in amphibian egg mortality and potential links to population declines.

Landscape genetics of the blotched tiger salamander (Ambystoma tigrinum melanostictum).

The field of landscape genetics has great potential to identify habitat features that influence population genetic structure. To identify landscape correlates of genetic differentiation in a quantitative fashion, we developed a novel approach using geographical information systems analysis. We present data on blotched tiger salamanders (Ambystoma tigrinum melanostictum) from 10 sites across the northern range of Yellowstone National Park in Montana and Wyoming, USA. We used eight microsatellite loci to analyse population genetic structure. We tested whether landscape variables, including topographical distance, elevation, wetland likelihood, cover type and number of river and stream crossings, were correlated with genetic subdivision (F(ST)). We then compared five hypothetical dispersal routes with a straight-line distance model using two approaches: (i) partial Mantel tests using Akaike's information criterion scores to evaluate model robustness and (ii) the BIOENV procedure, which uses a Spearman rank correlation to determine the combination of environmental variables that best fits the genetic data. Overall, gene flow appears highly restricted among sites, with a global F(ST) of 0.24. While there is a significant isolation-by-distance pattern, incorporating landscape variables substantially improved the fit of the model (from an r2 of 0.3 to 0.8) explaining genetic differentiation. It appears that gene flow follows a straight-line topographic route, with river crossings and open shrub habitat correlated with lower F(ST) and thus, decreased differentiation, while distance and elevation difference appear to increase differentiation. This study demonstrates a general approach that can be used to determine the influence of landscape variables on population genetic structure.

A model for optimal reaction norms: the case of the pregnant garter snake and her temperature-sensitive embryos.

We present a model to test Osgood's ( 1978 ) proposition that viviparous snakes have optimal reaction norms for temperature-sensitive meristic traits, such as scale counts. Our model predicts that traits that are subject to temperature effects during development will evolve a flat or [Formula: see text]-shaped reaction norm (average scale count as a function of developmental temperature). We tested this prediction by maintaining 67 female garter snakes (Thamnophis elegans) at eight different constant temperatures (21 degrees -33 degrees C) during pregnancy and making a series of scale counts on their newborn offspring (n = 491). To insure that the experimental temperatures were ecologically relevant, we used automated radiotelemetry to record the body temperature of pregnant, free-ranging females. The resulting temperature data allowed us to test the prediction that the inflection points of reaction norms would correspond to the average temperature experienced by embryos in the field. In line with predictions of the Osgood model, reaction norms were flat or U-shaped. In the case of U-shaped reaction norms, the inflection point of the curves corresponded to the average temperature imposed on embryos by free-ranging females. In contrast to some past studies, none of the standard scale scores (ones commonly used in systematics) showed significant temperature effects in either sex. Reaction norms were flat. In contrast, incidences of various abnormalities showed U-shaped reaction norms. Temperature effects were more pronounced in males than in females. The results have implications for systematics and for the evolution of canalization and phenotypic plasticity.