Email updates

Keep up to date with the latest news and content from BMC Ecology and BioMed Central.

Open Access Research article

Developmental dynamics of Ambystoma tigrinum in a changing landscape

Sarah K McMenamin12* and Elizabeth A Hadly1

Author Affiliations

1 Department of Biology, Stanford University, Stanford CA 94305-5020, USA

2 Present address: Department of Biology, University of Washington, Box 351800, Seattle WA 98195-1800, USA

For all author emails, please log on.

BMC Ecology 2010, 10:10  doi:10.1186/1472-6785-10-10


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1472-6785/10/10


Received:15 April 2009
Accepted:3 April 2010
Published:3 April 2010

© 2010 McMenamin and Hadly; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Loss of pond habitat is catastrophic to aquatic larval amphibians, but even reduction in the amount of time a breeding site holds water (hydroperiod) can influence amphibian development and limit reproductive success. Using the landscape variation of a glacial valley in the Greater Yellowstone Ecosystem as the context for a natural experiment, we examined variation in growth pattern and life history of the salamander Ambystoma tigrinum melanostictum and determined how these developmental characteristics varied with hydroperiod over several summers.

Results

In ponds that dried early in the season, maximum larval size was reduced relative to the sizes achieved in permanent ponds. Ephemeral ponds were associated with early metamorphosis at small body sizes, while permanent ponds facilitated longer larval periods and later metamorphosis. Paedomorphosis resulted from indefinite metamorphic postponement, and was identified only in the most permanent environments. Patterns of growth and allometry were similar between ponds with different hydroperiods, but considerable life history variation was derived from modulating the timing of and size at metamorphosis. Considering maximum rates of growth and inferring the minimum size at metamorphosis across 25 ponds over the course of three years, we calculated that hydroperiods longer than three months are necessary to support these populations through metamorphosis and/or reproductive maturity.

Conclusions

Landscape heterogeneity fosters life history variation in this natural population. Modulation of the complex ambystomatid life cycle allows this species to survive in unpredictable environments, but current trends towards rapid pond drying will promote metamorphosis at smaller sizes and could eliminate the paedomorphic phenotype from this region. Metamorphosis at small size is has been linked to altered fitness traits, including reduced survival and fecundity. Thus, widespread environmental truncation of larval periods may lead to decreased population persistence. We found that the hydroperiods of many ponds in this region are now shorter than the developmental period required for larvae to reach the minimum size for metamorphosis; these locations serve as reproductive sinks that may be detrimental for persistence of the species in the region.

Background

Phenotype is determined by the interaction of genotype and developmental environment, but the influence of environment on variation within species is sometimes under-emphasized. While canalization produces consistent phenotypic outcomes under a variety of conditions [1], phenotypic plasticity produces developmental outcomes that vary with environment [2]. Vertebrate development, and particularly development of larval amphibians, can be modified by many factors, including temperature [3,4], density [5-7] and resource availability [4,8,9]. These extrinsic signals are transduced into endocrinological signals that ultimately control rates of growth and timing of developmental events [8,10].

The optimal phenology of developmental events depends on ecological context and conditions [11], and plasticity in developmental timing can allow organisms to maximally exploit a range of unpredictable environments [12-15]. When a complex life cycle involves larval and adult stages that inhabit different ecological niches, species can particularly benefit from plasticity in the timing of the metamorphic transition. Most amphibians require an aquatic environment during larval development, and developmental plasticity allows populations to maximally exploit both permanent and ephemeral larval habitats, with rapidly drying aquatic habitats associated with rapid metamorphosis [16-19]. However, the timing of metamorphosis affects phenotype, survival, age at first reproduction and reproductive success, and shorter larval periods may lead to decreased lifetime fecundity [5,18,20]. Moreover, there is thought to be a minimum size at which metamorphosis can occur (SVLm) [11,12,21,22], and if larvae cannot grow fast enough to reach SVLm before an environment becomes uninhabitable, larvae perish before recruitment. There is thus a apparent trade-off between metamorphosis at large sizes (maximizing terrestrial survival and fecundity) and metamorphosis at small sizes (minimizing larval mortality) [22,23].

If the resources necessary to support growth to SVLm are not available, individuals of facultatively paedomorphic species are hypothesized to become reproductively mature in a larval form ("best of a bad lot" progenic paedomorphosis) [21,24]. At the opposite end of the plasticity spectrum, if the larval environment is particularly resource-rich, individuals may permanently remain in the larval niche, becoming paedomorphic even after SVLm is reached. These "paedomorph advantage" individuals are expected to show higher fecundity than both progenic paedomorphs and metamorphosed individuals [21,24].

With alternative developmental strategies expressed under different environmental conditions, habitat heterogeneity across geographic space promotes diversity in phenotype and life history [16,25,26]. Developmental responses to environmental conditions are observed to vary through time as well, and millennial-scale climatic events catalyzed changes in the growth and developmental timing of Ambystoma tigrinum melanostictum populations in Yellowstone National Park (YNP) [27].

In this study, we characterized life histories, metamorphic timing and patterns of growth in A. t. melanostictum populations across a glacially modified landscape in the Greater Yellowstone Ecosystem (GYE). Recent evidence that ponds in this landscape are changing to the detriment of local amphibian populations [28] makes investigating the relationships between population status and landscape conditions particularly urgent. Focusing on populations in 25 heterogeneous pond environments, we characterized (1) the rates of larval growth, (2) the maximum sizes attained by larvae and the minimum sizes at maturity, and (3) the patterns of allometric growth as each varied with the hydroperiods of the different ponds. Considering the fitness associations of different developmental strategies and the limits of developmental flexibility of this species, we discuss the long-term implications of the climatic and environmental changes taking place in northern YNP.

Results

During the summer months of 2006-2008, we measured 609 A. t. melanostictum individuals captured from 25 pond environments in the GYE, including 22 within YNP (Additional file 1). Of the individuals sampled, 519 (85%) were aquatic larvae, 71 (12%) were terrestrial metamorphs and 19 (3%) were mature aquatic paedomorphs.

Additional file 1. Hydroperiods and collection dates. Total number of A. t. melanostictum individuals (larvae, metamorphs and paedomorphs) collected each day shown for every pond-year. Approximate hydroperiods (± 7 days) are shown in light blue; permanent ponds are blue throughout. Individuals collected outside of the hydroperiod represent metamorphs collected from the location after drying.

Format: PDF Size: 75KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Growth rate

We modelled the rate of larval growth in ephemeral (n = 405 larvae) and permanent ponds (n = 81 larvae) within YNP, and fit von Bertalanffy models of growth [29] to each group (Figure 1). The calculated tangent of the model was higher for larvae from permanent ponds (L= 75.1) than larvae from ephemeral ponds (L= 58.7), reflecting the fact that the larvae from permanent reached larger sizes. The rate constants and x-intercepts were similar between both models.

thumbnailFigure 1. Growth of salamanders over time during summers 2006-2008. A. t. melanostictum larvae collected in YNP during the summers of 2006-2008. Blue rings denote larvae collected from permanent environments (n = 81); gray rings denote larvae collected from ephemeral environments (n = 405). Curves show Von Bertalanffy models of growth: growth of larvae from permanent environments shown with blue line; growth of larvae from ephemeral ponds shown with gray. On the y axis, a shows date when the first ephemeral pond dried (30 June); b shows first date that an emerging metamorph was found in an ephemeral pond (10 July); c shows earliest date that an emerging metamorph was found in a permanent pond (19 July); d shows earliest date a recently mature paedomorph was found in a permanent pond (30 July). On the y axis, e shows the size of the smallest metamorph captured (45 mm ≅ SVLm).

Modelled as linear functions, larval growth was statistically indistinguishable between the two environments (Permanent: y = 21.6 + 0.87x, r2 = 0.72; Ephemeral: y = 18.2 + 0.87x, r2 = 0.57). We used linear models to describe early growth rates during 2008 prior to 30 June before any of the ponds dried completely. Considering the nine ponds with sufficient sampling to create realistic linear models of early larval growth, we found slopes ranging from 0.81 (Pond 2, dry in July) to 1.48 mm/day (Pond 37, dry in August). Based on these nine locations, we found no statistical differences in early growth rate between Hydroperiod categories (Figure 2).

thumbnailFigure 2. Early growth rates are similar in ponds of different hydroperiods. Slope of growth rates (mm/day) of YNP larval salamanders in 2008 prior to 30 June.

Hatching dates and larval period

We used the nine independent models of early growth to infer the dates at which pond populations hatched in 2008. The projected hatching dates ranged from 1 June (Pond 2, dry in July) to 14 June (Pond 38, permanent; Table 1). The smallest metamorph identified during the course of this study was 45 mm SVL, and we interpret this size as close to the minimum size at which metamorphosis can occur (SVLm). Growing at the most rapid recorded linear growth rate (1.48 mm/day, Pond 37), it takes 33 days from hatching at 5 mm to reach SVLm (~45 mm).

Table 1. Characteristics of different ponds in 2006 through 2008

Snowmelt and spawning occur in early April and hatching occurs approximately 55 days later in early June. Thus, in this region, a hydroperiod of more than 88 days (from April snowmelt to early July) is necessary to support this species through metamorphosis. Our model predicts that few if any individuals will reach metamorphosis before July. Indeed, the earliest young-of-the-year metamorph was captured on 11 July.

Size of Larvae and Size at Metamorphosis

The average size of larvae from permanent ponds was more than 15 mm larger than the average size of larvae from ephemeral ponds. However, this comparison was biased by the fact that collection occurred earlier in ephemeral ponds. We employed a generalized linear mixed effects model [30] considering SVL of YNP larvae as a function of Year of Collection, Date of Collection and Hydroperiod, with random slopes for Date within Pond and Year within Pond (AIC = 6161, log likelihood = -3069). This model performed somewhat better than the same model without the Hydroperiod variable (AIC = 6170, log likelihood = -3076), but an ANOVA showed that the two models were not statistically different from one another. Thus, when Date of Collection and differences in growth between Ponds and Years were factored out, Hydroperiod had no measurable effect on larval size. Results were unchanged when Hydroperiod was assessed as a binary variable (permanent/ephemeral).

We compared the sizes of the largest larvae found in each pond during each year, including largest larvae from all permanent pond-years and largest larvae from ephemeral pond-years within three weeks of pond drying (Figure 3, black crosses). ANOVA revealed significant differences between Hydroperiod categories (df = 3, f = 3.0, p-value < 0.05), although Tukey honest significant differences showed no significance between the four categories. However, when the largest larvae from permanent ponds were compared to the largest larvae from ephemeral ponds (Hydroperiod categories June, July and August), larvae from permanent ponds were ~14 mm larger than larvae from ephemeral ponds (t = -2.6, df = 14.1, p-value = 0.02). Paedomorphs were found only in permanent ponds (Figure 3, purple inverted triangles).

thumbnailFigure 3. Maximum larval size and minimum size at maturity increase with hydroperiod. Each point represents size category for a single pond-year 2006-2008. Black crosses show SVL of largest larva found in each pond-year; red triangles show SVL of smallest metamorph found in each pond-year; inverted blue triangles show SVL of smallest paedomorph in each pond-year. Pond-years categorized by hydroperiod into locations that dried in June, July, August, and locations that did not dry (Permanent).

Allometry and shape change

We performed a principal components analysis using the natural logarithms of eight different physical parameters (Table 2). All measurements loaded strongly onto PC1, with SVL contributing the strongest loading. PC2, the shape component, was dominated by tail height. Although PC2 explained relatively little of the total variance in shape (4.9%), PC2 along with PC1 strongly separated aquatic from terrestrial individuals (similar to Figure 4). Aquatic individuals (both larvae and paedomorphs) showed a negatively allometric relationship between SVL and tail height. Larvae and paedomorphs showed approximately a 3:1 relationship between SVL and tail height; this relationship changed to 8:1 after metamorphosis (Figure 4).

thumbnailFigure 4. Allometric relationships in aquatic and terrestrial salamanders. Open circles represent larvae (n = 373); filled circles show paedomorphs (n = 19); open triangles show terrestrial individuals (n = 39). Blue symbols denote individuals collected from permanent locations; gray symbols denote those collected from ephemeral ponds. Solid line and adjacent formula indicate best fit linear model for all aquatic individuals (larvae and paedomorphs); dashed line and adjacent formula indicate the best fit line for all terrestrial individuals.

Table 2. Principal components analysis of salamander physiology

Discussion

The A. t. melanostictum pond populations in this region of the GYE are closely related to one another, exchange gene flow, and exhibit low overall genetic diversity [31-33]. Therefore, majority of the phenotypic variation encountered in this region is due to phenotypic plasticity rather than genetic differentiation and adaptation. Measuring 609 individuals collected from 25 ponds within the GYE, we documented substantial variation in the duration of larval periods and the sizes at metamorphosis, differences which correlated with differing pond hydroperiods.

Larvae in more permanent environments reached greater sizes than larvae in ponds with brief hydroperiods (Figures 1 and 3). While most larvae from ephemeral ponds had either metamorphosed or perished before reaching 60 mm SVL, larvae in permanent environments were able to remain under aquatic growth conditions up to nearly 80 mm SVL, after which they either became paedomorphic or metamorphic. The smallest metamorph collected from a permanent environment was 62 mm SVL, 17 mm larger than the smallest metamorph found in an ephemeral environment (45 mm). Indeed nearly half of metamorphs collected from ephemeral environments were smaller than the smallest metamorph collected from any permanent environment. Furthermore, while the first emerging metamorph was captured from an ephemeral location on 10 July (Figure 1: b), the first metamorph emerged from a permanent pond more than a week later (and this individual was more than 80 mm SVL, suggesting that it may have not been young-of-the-year; Figure 1: c).

Many permanent environments contain over-wintering larvae, which postpone metamorphosis until their second or third years. Individuals under this developmental regime grow as aquatic larvae for several years, and some eventually forego metamorphosis to become reproductively mature neotenic paedomorphs ("paedomorph advantage") [21]. Only the two most stable environments examined, Ice Lake and Rainbow Lake, contained paedomorphs. These were the only two YNP locations that never dried over the course of the three year study period. We did not identify any small progenic paedomorphs (i.e. "best of a bad lot" paedomorphs [21]), possibly because the resource-poor locations that might have promoted precocious paedomorphosis were ephemeral and would have destroyed lingering aquatic individuals.

Larval growth rate determines the amount of time required to reach SVLm, and therefore defines the minimum larval period required prior to metamorphosis. Environmental conditions greatly influence larval growth rates in controlled settings [12,34]; however, we found average larval growth to be extremely similar between ponds and ponds of different hydroperiods. The earlier onset of metamorphosis observed in environments with shorter hydroperiods was not accompanied by accelerated larval growth rate. Between 16 June and 30 June, 2008 (during which all ponds were hydrated), individuals grew a linear average of more than 1 mm per day (Figure 2).

We used the early larval growth models to extrapolate dates of hatching in each of these ponds, all of which fell between 1 June and 14 June (Table 1). Assuming an incubation period of 55 days (the length of time between spawning and hatching in the permanent Ice Lake location [35]), we further extrapolated that spawning likely occurred in each of these ponds during the first two weeks of April, soon after snow melted from the region [36]. These inferred hatching dates are a month earlier than the hatching dates recorded at Ice Lake in 1993 [35]. Northern YNP has experienced warmer spring temperatures and earlier snowmelt dates in recent decades [28,36], suggesting that long-term warming has contributed to phenological changes and earlier breeding of this amphibian population [see also [37-39]].

At metamorphosis, amphibians undergo major restructuring of internal and external organ systems, remodelling somatic tissues to adopt a terrestrial lifestyle. In addition to the loss of external gills, the lateral tail swim fin is lost, resulting in a decrease in tail height. Tail height dominated the second component (shape) of the principal components analysis (PCA; Table 2), and tail height compared to size reliably separated aquatic from terrestrial individuals (Figure 4). After loss of the swim fin, terrestrial individuals showed far narrower tails relative to their length. Metamorphosis furthermore modified the relationship between size and tail shape. For every centimetre SVL growth gained by an aquatic salamander, it gained approximately 30 mm of tail height; while every centimetre of snout-vent growth gained by a terrestrial salamander was accompanied by only 12 mm of tail height. Paedomorphs maintained the aquatic SVL/tail height relationship throughout their lives. We have developed these observations into a conceptual framework of allometric relationships though different developmental pathways (Figure 5). As shown, the metamorphosis of an individual at a small size (Figure 5: A), at a large size (B) or not at all (C) is functionally determined by the hydroperiod of the pond environment.

thumbnailFigure 5. Allometric relationships through ontogeny. Aquatic and terrestrial individuals show different relationships between tail morphology (y axis) and size (x axis), represented by the two large grey arrows. A: Larvae in ephemeral environments undergo early metamorphosis and transform into terrestrial individuals at a small size. SVLm shows the minimum size for metamorphosis. B: Larvae in more permanent environments delay metamorphosis and transform into terrestrial adults at larger body sizes. C: In permanent aquatic environments, some larvae forgo metamorphosis to become reproductively mature neotenic paedomorphs.

Individuals that are able to delay metamorphosis (Figure 5: B) may be more likely to survive and reproduce than individuals that metamorphose at or close to SVLm (Figure 5: A) [5,18-20,40]. Although early metamorphosis or progenic paedomorphosis allow immediate reproduction [21,41], studies show that salamanders that delay metamorphosis to larger body sizes enjoy overall benefits in survival and fecundity [18,40]. Larger individuals are at less risk of desiccation and predation in terrestrial environments [12], are able to more efficiently capture terrestrial food, have higher mating success, larger clutch sizes and generally higher fecundity than smaller individuals [24,40]. As ponds in this region become increasingly ephemeral [28], resident amphibians may be forced to adopt ever more rapid and potentially less fecund developmental profiles.

Although differences in adult survival and fecundity have not yet been tested in this region, we determined that many ponds in the area are now entirely unsuitable for Ambystoma development, leaving the population vulnerable to continued warming. Larvae must reach a certain minimum size in order to undergo successful metamorphosis [19,21,24,42], and based on field observations, the SVLm for this population is close to 45 mm. If the hydroperiod is shorter than the minimum larval period required to reach SVLm, larvae have no opportunity to escape a drying pond. Growing at the most rapid documented growth rate, the minimum larval period is approximately 33 days from hatching, and the incubation period of eggs is approximately 55 days. Thus, these salamanders need a minimum of 88 days of pond hydration achieve SVLm. Snowmelt and egg deposition occur in early April, and ponds must therefore contain water through approximately the first week in July to support populations through metamorphosis. As predicted by this model, the first young-of-the-year metamorph was captured on 11 July. Earlier snowmelt [36] could extend the functional hydroperiod, potentially buffering the effects of earlier drying, however, such changes in seasonality may negatively impact the population in other ways [37-39].

Over the course of three years, we observed six locations to dry during the month of June, and hundreds of bodies of young salamanders were found trapped in the drying wetlands [28]. These ponds with hydroperiods shorter than the required larval period serve as reproductive sinks, where few if any larvae survive to recruitment. Larvae in ephemeral environments do not appear to accelerate larval growth rate, and this species may thus be poorly equipped to thrive in newly ephemeral, rapidly drying environments.

A. t. melanostictum is a long-lived, opportunistically breeding organism, and populations are able to recover from dry years [43]. Indeed, although past droughts likely contributed to population bottlenecks and loss of genetic diversity [31], the GYE population has recovered from both ancient climatic fluctuations [27] and more recent dry periods during the late 1980s and the early 2000s [28]. Nonetheless, long term drought conditions can prevent breeding over the course of decades, exceeding both the lifespan of the organism and the capacity of populations to recover.

Conclusions

The adaptable, complex lifecycle of A. t. melanostictum allows populations to exploit a variety of hydrologically diverse habitats. Nonetheless, there are fitness costs associated with the rapid metamorphosis necessary for development in ephemeral environments. Even though the timing of metamorphosis is highly plastic in this species, we found that the rate of larval growth is consistent between environments, and that there is an absolute minimum amount of time required for larvae to reach the minimum size for metamorphosis. Several ponds in northern YNP now have hydroperiods too short to support the minimum larval period, and serve as reproductive sinks. If this landscape continues to dry and more locations become functionally uninhabitable, the population may be pushed beyond the limits of its adaptive and reproductive capacity.

Methods

Study area and collection

We focused this study on 22 bodies of water in northern YNP, as well as 3 lakes northwest of the northern park boundary but within the GYE (Axolotl 1, 2, and 3). These 25 locations represented a variety of developmental habitats for larval amphibians, ranging from small, ephemeral kettle ponds fed by local aquifers that dried early in the summer, to permanent ponds larger than seven hectares. Elevation and coordinates of each location was recorded using a handheld GPS unit; maximum area was estimated from satellite photos (characteristics and locations in Table 1). During every summer 2006-2008, each pond was characterized as having dried to 3 cm or less of standing water during June, July, August or never (permanent; see Additional file 1). Not all ponds were sampled every year of the study.

Each summer we captured and released between 1 and 52 (mean = 16.5) immature larval, paedomorph and terrestrial A. t. melanostictum individuals from in and around each sampling location. Most individuals were captured by active dip-netting, some were captured using aquatic funnel or terrestrial pitfall trapping. SVL of each individual was measured to the nearest mm. Individuals captured during 2007 and 2008 were also measured for limb length, tail length, tail height, tail width, head width and body width.

Metamorphic status was determined based on gill morphology: individuals with fully brachiated gills were considered larvae if they were small and immature and were considered paedomorphs if their external cloacae appeared mature. Individuals lacking gills or with gills in the process of resorption were considered to be terrestrial adults.

Size and growth assessments

All analyses were performed in R 2.8.1 [44]. We assessed average growth of all larvae from the 22 ponds within YNP. Ponds outside of the park were excluded from growth rate analyses because they were at a higher elevation and spawned later in the year. Average growth was fit to the von Bertalanffy model of growth [11,29]:

where L(t) is size at time t, Lis the mean maximum size, t0 is the x-axis intercept and k is the rate constant. We used linear models to describe early growth (before 30 June) of nine individual ponds during 2008. We used these linear models to approximate the date at which the populations hatched in 2008. This species is generally 5 mm SVL at hatching, so we extrapolated the model to the date at which SVL = 5.

We fit the sizes of YNP larvae to a generalized linear mixed model using the lmer function in the R package lme4 [30]. This mixed model assessed larval SVL as it varied with the categorical variable Year, the categorical variable Hydroperiod (June, July, August, Never), and with Date of Collection as a continuous variable. Since observations of salamanders were clustered by Ponds and by Years within Ponds, we included a random intercept for Pond and a random interaction between Year and Pond. In order to capture the possibility of different slopes for date across ponds, we also included random slopes for Date within Ponds.

The largest larva from each pond-year was recorded for ephemeral ponds from which larvae had been collected within 3 weeks (21 days) of drying. Thus, the largest larva from Pond 34 in 2008 was excluded because the last larva was collected from that location more than 5 weeks before the pond dried. Largest larvae from permanent ponds were considered each year regardless of when individuals were collected. Smallest metamorphs and paedomorphs were determined for every pond-year.

To characterize differences in growth patterns through ontogeny, we plotted tail height against SVL. We performed a PCA using as variables the natural log-transformed measurements of SVL, average front leg length, average rear leg length, tail length, tail height, tail width, head width and body width.

Abbreviations

(GYE): Greater Yellowstone Ecosystem; (YNP): Yellowstone National Park; (SVL): Snout vent length; (PCA): Principal components analysis.

Authors' contributions

EAH and SKM designed the study. SKM collected specimens, analyzed data and prepared the manuscript. Both authors have read and approved the final manuscript.

Acknowledgements

We would like to thank all current and former members of the Hadly lab: Jessica Blois, Lily Cheng Li, Rebecca Terry and Andrew Romiger for feedback and comments. We thank Bill Anderegg, Lauren Palumbi, Laura Ha, Yu-Jin Lee and Alice Eamsherangkoon for assistance with fieldwork. We also acknowledge and thank the staff and rangers at YNP, particularly Christie Hendrix, for access and assistance. Thanks to Veronika Skrivankova and John Boik statistical advice, to S. Randal Voss, Steve Palumbi and Terry Root for valuable discussion, and to two anonymous reviewers for thoughtful comments. Funding was provided through grants from the Sigma Xi Grant-In-Aid program, the Stanford Center for Evolutionary Biology, the Gifford Fund and the Stanford Office of the Vice Provost.

References

  1. Waddington CH: Canalization of Development and the Inheritance of Acquired Characters.

    Nature 1942, 150(3811):563. Publisher Full Text OpenURL

  2. Gould SJ: Ontogeny and phylogeny. Belknap Press of Harvard University Press Cambridge, Mass; 1977.

  3. Bizer JR: Growth rates and size at metamorphosis of high elevation populations of Ambystoma tigrinum.

    Oecologia 1978, 34(2):175-184. Publisher Full Text OpenURL

  4. Sprules WG: Environmental factors and the incidence of neoteny in Ambystoma gracile (Baird) (Amphibia: Caudata).

    Canadian Journal of Zoology 1974, 52:1545-1552. PubMed Abstract | Publisher Full Text OpenURL

  5. Semlitsch RD: Density-Dependent Growth and Fecundity in the Paedomorphic Salamander Ambystoma talpoideum.

    Ecology 1987, 68(4):1003-1008. Publisher Full Text OpenURL

  6. Harris RN: Density-dependent paedomorphosis in the salamander Notophthalmus viridescens dorsalis.

    Ecology 1987, 68:705-712. Publisher Full Text OpenURL

  7. Provenzano S, Boone M: Effects of Density on Metamorphosis of Bullfrogs in a Single Season.

    J Herpetol 2009, 43(1):49-54. Publisher Full Text OpenURL

  8. Ryan TJ, Semlitsch RD: Growth and the expression of alternative life cycles in the salamander Ambystoma talpoideum (Caudata: Ambystomatidae).

    Biol J Linn Soc 2003, 80(4):639-646. Publisher Full Text OpenURL

  9. Voss SR: Genetic Basis of Paedomorphosis in the Axolotl, Ambystoma mexicanum: A Test of the Single-Gene Hypothesis.

    J Hered 1995, 86(6):441. OpenURL

  10. Denver RJ: Environmental stress as a developmental cue: corticotropin-releasing hormone is a proximate mediator of adaptive phenotypic plasticity in amphibian metamorphosis.

    Horm Behav 1997, 31:169-179. PubMed Abstract | Publisher Full Text OpenURL

  11. Day T, Rowe L: Developmental thresholds and the evolution of reaction norms for age and size at life-history transitions.

    Am Nat 2002, 159(4):338-350. PubMed Abstract | Publisher Full Text OpenURL

  12. Wilbur HM, Collins JP: Ecological Aspects of Amphibian Metamorphosis: Nonnormal distributions of competitive ability reflect selection for facultative metamorphosis.

    Science 1973, 182(4119):1305-1314. PubMed Abstract | Publisher Full Text OpenURL

  13. Griffiths RA: Temporary ponds as amphibian habitats.

    Aquat Conserv: Mar Freshwat Ecosyst 1997, 7(2):119-126. Publisher Full Text OpenURL

  14. Sultan SE, Spencer HG: Metapopulation Structure Favors Plasticity over Local Adaptation.

    Am Nat 2002, 160(2):271-283. PubMed Abstract | Publisher Full Text OpenURL

  15. Newman RA: Adaptive Plasticity in Amphibian Metamorphosis.

    Bioscience 1992, 42(9):671-678. Publisher Full Text OpenURL

  16. Márquez-García M, Correa-Solis M, Sallaberry M, Méndez M: Effects of pond drying on morphological and life-history traits in the anuran Rhinella spinulosa (Anura: Bufonidae).

    Evol Ecol Res 2009, 11:803-815. OpenURL

  17. Richter-Boix A, Llorente G, Montori A: Effects of phenotypic plasticity on post-metamorphic traits during pre-metamorphic stages in the anuran Pelodytes punctatus.

    Evol Ecol Res 2006, 8(2):309-320. OpenURL

  18. Semlitsch RD, Scott DE, Pechmann JHK: Time and Size at Metamorphosis Related to Adult Fitness in Ambystoma Talpoideum.

    Ecology 1988, 69(1):184-192. Publisher Full Text OpenURL

  19. Semlitsch RD, Wilbur HM: Effects of pond drying time on metamorphosis and survival in the salamander Ambystoma talpoideum.

    Copeia 1988, 1988(4):978-983. Publisher Full Text OpenURL

  20. Scott DE: The Effect of Larval Density on Adult Demographic Traits in Ambystoma Opacum.

    Ecology 1994, 75(5):1383-1396. Publisher Full Text OpenURL

  21. Whiteman HH: Evolution of Facultative Paedomorphosis in Salamanders.

    Q Rev Biol 1994, 69(2):205-221. Publisher Full Text OpenURL

  22. Rudolf VHW, Rödel MO: Phenotypic plasticity and optimal timing of metamorphosis under uncertain time constraints.

    Evol Ecol 2007, 21(1):121-142. Publisher Full Text OpenURL

  23. Werner EE: Amphibian Metamorphosis: Growth Rate, Predation Risk, and the Optimal Size at Transformation.

    Am Nat 1986, 128(3):319-341. Publisher Full Text OpenURL

  24. Denoël M, Joly P, Whiteman HH: Evolutionary ecology of facultative paedomorphosis in newts and salamanders.

    Biol Rev 2005, 80(04):663-671. PubMed Abstract | Publisher Full Text OpenURL

  25. Morrison C, Hero JM: Geographic variation in life-history characteristics of amphibians: a review.

    J Anim Ecol 2003, 72(2):270-279. Publisher Full Text OpenURL

  26. Relyea RA, Werner EE: Morphological Plasticity in Four Larval Anurans Distributed along an Environmental Gradient.

    Copeia 2000, 2000(1):178-190. Publisher Full Text OpenURL

  27. Bruzgul JE, Long W, Hadly EA: Temporal response of the tiger salamander (Ambystoma tigrinum) to 3,000 years of climatic variation.

    BMC Ecol 2005, 5(1):7. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  28. McMenamin SK, Hadly EA, Wright C: Climatic change and wetland desiccation cause amphibian decline in Yellowstone National Park.

    Proc Natl Acad Sci USA 2008, 105(44):16988-16993. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. von Bertalanffy L: Quantitative laws in metabolism and growth.

    Q Rev Biol 1957, (32):217-231. PubMed Abstract OpenURL

  30. Bates D, Maechler M, Dai B: lme4: Linear mixed-effects models using S4 classes. [http://lme4.r-forge.r-project.org/] webcite

    2008.

  31. Spear SF, Peterson CR, Matocq MD, Storfer A: Molecular evidence for historical and recent population size reductions of tiger salamanders (Ambystoma tigrinum) in Yellowstone National Park.

    Conserv Genet 2006, 7(4):605-611. Publisher Full Text OpenURL

  32. Spear SF, Peterson CR, Matocq MD, Storfer A: Landscape genetics of the blotched tiger salamander (Ambystoma tigrinum melanostictum).

    Mol Ecol 2005, 14(8):2553-2564. PubMed Abstract | Publisher Full Text OpenURL

  33. McMenamin SK: Landscape & Yellowstone amphibians: Insights from genetic, developmental and community variation through space and time. PhD thesis. Stanford University, Department of Biology; 2009. OpenURL

  34. Semlitsch RD, Harris RN, Wilbur HM: Paedomorphosis in Ambystoma talpoideum: Maintenance of Population Variation and Alternative Life-History Pathways.

    Evolution 1990, 44(6):1604-1613. Publisher Full Text OpenURL

  35. Hill SR: Migratory chronology of adult tiger salamanders (Ambystoma tigrinum) and survey of larvae of the tiger salamander in the northern range of Yellowstone National Park. In MS thesis. Montana State University; 1995. OpenURL

  36. Wilmers CC, Getz WM: Gray wolves as climate change buffers in Yellowstone.

    PLoS Biol 2005, 3(4):e92. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Corn PS: Amphibian Breeding and Climate Change: Importance of Snow in the Mountains.

    Conserv Biol 2003, 17(2):622-625. Publisher Full Text OpenURL

  38. Beebee TJC: Amphibian breeding and climate.

    Nature 1995, 374(6519):219-220. Publisher Full Text OpenURL

  39. Blaustein AR, Belden LK, Olson DH, Green DM, Root TL, Kiesecker JM: Amphibian Breeding and Climate Change.

    Conserv Biol 2001, 15(6):1804-1809. Publisher Full Text OpenURL

  40. Semlitsch RD, Gibbons JW: Phenotypic Variation in Metamorphosis and Paedomorphosis in the Salamander Ambystoma talpoideum.

    Ecology 1985, 66(4):1123-1130. Publisher Full Text OpenURL

  41. Denoël M, Joly P: Neoteny and progenesis as two heterochronic processes involved in paedomorphosis in Triturus alpestris (Amphibia: Caudata).

    Proc R Soc Lond, Ser B: Biol Sci 2000, 267:1481-1485. Publisher Full Text OpenURL

  42. Petranka JW: Evolution of complex life cycles of amphibians: bridging the gap between metapopulation dynamics and life history evolution.

    Evol Ecol 2007, 21(6):751-764. Publisher Full Text OpenURL

  43. Taylor B, Scott D, Gibbons J: Catastrophic Reproductive Failure, Terrestrial Survival, and Persistence of the Marbled Salamander.

    Conserv Biol 2006, 20(3):792. PubMed Abstract | Publisher Full Text OpenURL

  44. R Development Core Team: R: A language and environment for statistical computing. [http://www.R-project.org] webcite

    Vienna, Austria: R Foundation for Statistical Computing; 2008.