Predator mediated selection and the impact of developmental stage on viability in wood frog tadpoles (Rana sylvatica)
© Calsbeek and Kuchta; licensee BioMed Central Ltd. 2011
Received: 10 August 2011
Accepted: 7 December 2011
Published: 7 December 2011
Complex life histories require adaptation of a single organism for multiple ecological niches. Transitions between life stages, however, may expose individuals to an increased risk of mortality, as the process of metamorphosis typically includes developmental stages that function relatively poorly in both the pre- and post-metamorphic habitat. We studied predator-mediated selection on tadpoles of the wood frog, Rana sylvatica, to identify this hypothesized period of differential predation risk and estimate its ontogenetic onset. We reared tadpoles in replicated mesocosms in the presence of the larval odonate Anax junius, a known tadpole predator.
The probability of tadpole survival increased with increasing age and size, but declined steeply at the point in development where hind limbs began to erupt from the body wall. Selection gradient analyses indicate that natural selection favored tadpoles with short, deep tail fins. Tadpoles resorb their tails as they progress toward metamorphosis, which may have led to the observed decrease in survivorship. Path models revealed that selection acted directly on tail morphology, rather than through its indirect influence on swimming performance.
This is consistent with the hypothesis that tail morphology influences predation rates by reducing the probability a predator strikes the head or body.
Many organisms exploit different environments over the course of their life cycle. Perhaps the most extreme example of this shift in resource use is that which accompanies metamorphosis in animals with complex life cycles . Complex life cycles - hereafter referring to organisms with at least two discrete post-embryonic life-stages [2, 3] - are ubiquitous in animals, being expressed in at least 80% of all species [4, 5]. They may evolve for several reasons, such as trophic switching or specialized dispersal/breeding forms . The tradeoffs that accompany shifts in niche occupancy will typically be accompanied by divergent selective regimes and alternative adaptations. In part, this accounts for the large differences in morphology, physiology, behavior, and other aspects of the phenotype observed among life stages. Although dramatically divergent morphologies among different life stages allow individuals to exploit multiple kinds of resources throughout ontogeny, complex life cycles also involve functional trade-offs and thereby create a new problem: how to optimize the transition between life stages [7, 8].
The challenge of adapting to multiple adaptive peaks can be partially resolved by genetic and developmental decoupling among life stages . Nonetheless, it is often the case that genetic, developmental, and functional correlations persist across life stages (e.g., [9–14]). Moreover, even if there is complete adaptive decoupling of divergent life stages, the transitional period between life stages is still likely to be a performance trough that exposes individuals to increased risks. Indeed, the more differentiated the life stages, the more intense the risks are likely to be. Metaphorically, the transition from juvenile to adult may be viewed as movement between alternative peaks on an individual's adaptive landscape [15, 16], where peaks represent correspondence between an individual's phenotypic traits and the local maximum probability of survival.
Many amphibians exhibit a complex life cycle in which larval development (intervals of which are referred to as Gosner stages in frog tadpoles; ) is followed by metamorphosis into an adult form [5, 18]. Tadpoles are highly specialized for feeding, and the tadpole body plan consists mostly of a globose body and a sheet-like, laterally compressed tail . During metamorphosis, the tail is resorbed as hind and forelimbs emerge, thereby facilitating the transition from an aquatic swimming form (undulatory, axial locomotion) to a terrestrial hopping form (saltatory, appendicular locomotion). It is during the intermediate stages of metamorphosis that individuals are thought to experience increased predation risk . The hypothesized period of increased predation risk separating larval and adult forms derives from the observation that metamorphs are optimized for neither larval nor adult niches . For example, emergent hind limbs may impose drag and reduce swimming performance [19–21] and residual tail tissue may negatively impact saltatory locomotion . For instance, metamorphosing chorus frogs, Pseudacris triseriata, are more likely to be captured by predatory garter snakes, Thamnophus sirtalis, than are either tadpoles or adult frogs . Laboratory selection experiments on tadpoles likewise suggest the presence of a performance decline at metamorphosis . Field experiments designed to measure both natural selection and variation in viability during ontogenetic stages near the developmental switch between life-stages, however, are still lacking.
Here, we test a set of related hypotheses about variation in survival probability in the wood frog, Rana sylvatica. We begin broadly, by first testing whether fitness (i.e., survival) correlates with morphology across tadpole development [22–25]. This first question is designed to test the hypothesis of increased predation risk during tadpole metamorphosis. We next use path analytic models to compare alternative hypotheses regarding the causal structure underlying selection on tadpole morphology. This includes a test of the hypothesis that tail morphology is subject to selection via its effect on swimming performance [24, 26], which may be important for predator escape. We also consider the alternative hypothesis that tail shape may, as has been demonstrated previously, enhance survival by serving as a "lure" to attract predatory attacks towards the tail, thereby reducing the probability of mortal wounds to the head/body region [25, 27–30].
We collected tadpoles of the wood frog, Rana sylvatica, from a single pond near Randolph, Vermont, USA (43°54' N, 72°38' W) on June 11, 2010. The pool naturally contained larval odonates and other predatory invertebrates (Calsbeek and Kuchta, pers. obs.). Tadpoles were held overnight in 5 gallon food-grade plastic buckets with filtered pond water, and were fed an ad libitum diet of boiled lettuce. The morning after capture, we individually marked each tadpole with a unique color-coded combination of elastomer dyes (visible elastomer implants available from Northwest Marine Technologies, Shaw Island WA, U.S.A.) that we injected into the dorsal half of the tail fin, posterior to the body wall. Tadpoles were immobilized (but not anesthetized) during the marking procedure by holding them in a plastic multi-channel pipette well. We scored each tadpole's developmental stage  with the aid of a dissecting microscope just prior to the initiation of the selection experiment (mean Gosner stage = 34 ± 4.62 SD).
Tadpoles were then individually transferred to a V-shaped glass tank (which imposed a consistent orientation on the tadpoles) with a size standard, and were digitally photographed. We used digital images of each tadpole to make the following linear measurements: Head length: the distance from the anterior tip of the snout to the junction of the body with the tail; head height: the depth of head at its tallest point; tail length: from the junction of the tail with the body wall to the distal tip of the tail; tail muscle height: muscle height at the tallest point of the tail muscle; and tail fin height: fin height at the tallest part of the tail.
We measured swimming performance for half of the individuals in our selection study (N = 200 tadpoles) using a small (36L × 26W × 5H cm) tank containing filtered pond water and a size standard. Rapid development among the tadpoles held in buckets prevented us from measuring swimming performance for the remaining 200 individuals. Swimming trials were videotaped at 250 frames/sec using a high-definition digital camcorder (JVC Evario GZ-HM550-bu). Each tadpole was introduced to the swimming chamber and then motivated to initiate a "C-start" by touching the junction point between the tail fin and the body wall using a small metal pointer. We recorded three C-starts for each tadpole and used the fastest of these trials to estimate swimming performance, recording average speed over the 50 fastest frames. We chose to use this measure in our selection analyses because fifty frames was the average time required to swim one body length, and we assume that this is a good metric for predator avoidance. Burst speed was measured along the path of the tadpole movement using MotionAnalysis software (available from M. Chappell, University of California, Riverside, CA, U.S.A. http://warthog.ucr.edu/). We used the tadpole eye as a landmark for tracking individuals. All capturing, marking, photography, and swimming performance trials were conducted within 36 hours and the tadpoles were immediately transferred to cattle tanks for the selection experiment.
We conducted our selection experiment using eight 1136 L (300 gallon) cattle tanks that were randomly selected from an array of 49 tanks housed in an open field near the Dartmouth College campus. One month before introducing tadpoles, cattle tanks were cleaned and filled with ground water, 0.550 kg of dried Oak leaf litter, 15.4 g of rabbit chow (a nutrient source), and a three-liter aliquot of mixed zooplankton and phytoplankton collected from a pond near Norwich, VT (43.73° N, 72.31°W). We added five larval dragonflies (Anax junius) to each tank to serve as predators on tadpoles. To provide developing frogs with retreat sites, we placed three White Water Lilly (Nymphea ordorata) fronds on the water surface of each tank. To provide dragonfly larvae with perches, we used stones to anchor three to five tree branches (~100 cm) to the bottom of each tank. Finally, we randomly assigned 50 tadpoles to each of the eight tanks. Set up this way, the cattle tanks functioned as self-sustaining mesocosms that mimicked conditions experienced by tadpoles in nature . We covered each cattle tank with 0.5 × 0.5 cm hardware cloth pulled taut and secured with elastic cords. This functioned to shade the tanks, prevent predation, ensure that no metamorphosing individuals escaped, and preclude large, predatory insects from laying eggs.
We recorded the identity of individual surviving tadpoles in our selection experiment five and fourteen days following introduction to the artificial ponds. Survival was scored by removing all the leaf litter and filtering each tank to recover tadpoles with hand-held dip nets. We also verified the presence of all five dragonfly-larvae in each mesocosm (one dragonfly-larva in each of three tanks was replaced to account for single dead individuals). After the first census, we replaced the leaf litter, dragonfly-larvae, and tadpoles, and re-covered the tanks with the shade cloth. Following the second census, all tadpoles were brought back to the laboratory, sacrificed with an overdose of MS-222, and stored in 70% ethanol.
Comparison of alternative selection models showing the number of parameters (k) in each model.
Traits in the model
Tail Muscle Height
Tail Muscle Height2
Tail Fin Height2
Linear selection gradients were calculated from models that included only linear terms, whereas quadratic gradients (e.g., non-linear selection) and cross-product terms (i.e., correlational selection) were calculated from models that included both the linear and quadratic terms. Because the GLM underestimates quadratic terms by half [36, 37], quadratic gradients and their standard errors were doubled. Though parametric statistics provide robust estimates of selection gradients and other parameters [32, 38], these tests may be violated by survival data (live/die), which tend to have non-normally distributed errors [39, 40]. We computed significance values for selection gradients using generalized linear models including a logit link function . Prior to pooling data from individual tanks (i.e., replicates), we tested for any interaction between relevant terms and the factor for tank. None of these were significant, indicating that selection operated in the same way in all replicates. We dropped the interaction terms but retained a factor for "tank" in our models. The factor for tank explained a significant portion of the variance in all full models (0.02 > P < 0.03), but not in reduced models (0.06 > P < 0.08). All variables used in selection analyses were standardized to a mean of zero with unit standard deviation, except our fitness variable (survival), which was scaled by the mean [32, 42]. The degree of multi-colinearity among traits was assessed by estimating variance inflation factors (VIF; ), all of which were less than five. We visualized fitness surfaces using cubic splines .
In a second analysis, our Model 3 (Figure 1) was iteratively reduced to its significant components by sequentially setting causal paths with the lowest partial regression coefficients and the highest P values to zero . Alternative models were compared using Akaike's information criteria (AIC, ), including the difference between the preferred model and each subsequent model (Δi), normalized relative likelihoods (w i ), and evidence ratios [35, 51].
In addition to path analysis using maximum likelihood, we also conducted Bayesian analyses of the data. We did this to account for the binomial distribution of our fitness variable (survival), which likely violates the assumption of normal errors and multivariate normality in least-square calculations [52, 53]. Bayesian analysis in Amos 18  employs a Markov Chain Monte Carlo (MCMC) algorithm for estimating posterior distributions, and properly accounts for the binomial status of our fitness variable. Parameter estimates were estimated from 150,000 generations following a burn-in of 500 generations. Convergence of the MCMC algorithm was assessed using the convergence statistic developed by Gelman et al.  and implemented in AMOS . The significance of parameter estimates was assessed using 95% Bayesian credibility intervals. The results of the Bayesian analysis of the path coefficients were very similar to the maximum likelihood estimates (data not shown) and will not be presented.
Linear (β) and quadratic (γ) selection on all morphological traits measured in this study in a data set that included all individuals (Full, N = 172), as well as in a second dataset from which we excluded all individuals whose Gosner stage was > 39 at the time of release were excluded (Reduced, N = 151).
β ± SE
γ ± SE
0.08 ± 0.16
- 0.06 ± 0.16
0.10 ± 0.12
0.14 ± 0.14
- 0.16 ± 0.13
0.004 ± 0.16
Tail muscle Height
- 0.14 ± 0.15
- 0.10 ± 0.14
Tail Fin Height
0.22 ± 0.09*
- 0.06 ± 0.11
- 0.01 ± 0.07
- 0.06 ± 0.10
0.05 ± 0.19
- 0.10 ± 0.17
0.09 ± 0.13
0.14 ± 0.15
- 0.17 ± 0.14
0.04 ± 0.19
Tail muscle Height
- 0.08 ± 0.16
- 0.03 ± 0.17
Tail Fin Height
0.21 ± 0.10*
- 0.05 ± 0.12
- 0.05 ± 0.08
- 0.10 ± 0.11
Linear (β) and quadratic (γ) selection on all morphological traits except burst swimming speed.
β ± SE
γ ± SE
- 0.01 ± 0.11
- 0.18 ± 0.10
0.10 ± 0.09
0.04 ± 0.09
- 0.19 ± 0.07*
- 0.16 ± 0.06*
Tail muscle Height
0.02 ± 0.09
- 0.03 ± 0.09
Tail Fin Height
0.24 ± 0.06***
- 0.07 ± 0.06
- 0.08 ± 0.12
- 0.24 ± 0.11*
0.07 ± 0.10
0.07 ± 0.10
- 0.14 ± 0.08
- 0.11 ± 0.08
Tail muscle Height
0.07 ± 0.10
-0.03 ± 0.10
Tail Fin Height
0.16 ± 0.07*
-0.08 ± 0.08
Linear (β) and quadratic (γ) selection on morphological traits used in a model chosen based on AIC scores.
β ± SE
γ ± SE
0.05 ± 0.09
- 0.20 ± 0.07*
- 0.17 ± 0.07*
- 0.16 ± 0.06*
Tail Fin Height
0.27 ± 0.05***
- 0.007 ± 0.10
- 0.20 ± 0.08*
- 0.11 ± 0.08
- 0.10 ± 0.08
Tail Fin Height
0.20 ± 0.07***
In the full data set, the mean percentage survival (± SE) in each tank to the first census period (5 days) was 0.77 ± 0.02 (range 0.66-0.90). By the second census (14 days), mean survival had decreased to 0.57 ± 0.02 (range 0.48-0.66). Qualitatively, selection results during the two time periods were nearly identical (data not shown), but to maximize our power to detect selection, and to simplify the presentation of results, we use viability estimates from the second census as our measure of fitness. Frequent bite marks on the tails of surviving tadpoles suggest that dragonfly larvae were a key source of mortality in our study populations. We also recovered two complete tails during our census, with elastomer tags still intact, from tadpoles that did not survive. We conclude that mortality in the selection replicates was largely due to predation by dragonfly larvae.
Swimming speed was positively correlated with developmental stage (r2 = 0.14, df = 186, P < 0.0001) and tail length (r2 = 0.12; df = 186, P < 0.0001), and there was a weak quadratic relationship between swimming speed and tail fin height (individuals of intermediate tail fin height swam fastest: r2 = 0.07, df = 186, P = 0.051). However, we did not detect any selection on swimming speed in our experiment. In a model that included linear terms for tail length and tail fin height, the selection gradient for swimming speed was weakly negative and non-significant (β = -0.03 ± 0.07, P = 0.65). Even when we considered selection on swimming speed alone (i.e., the selection differential for swimming speed) we detected no variation in survival that was related to swimming performance (s = -0.07 ± 0.06, P = 0.23).
Comparison of path models.
Maximum Likelihood Analyses
Morphology -- > Fitness
3.16 × 10-4
Morphology -- > Performance -- > Fitness
1.48 × 10-12
6.74 × 1011
Comparison of path models.
Maximum Likelihood Analyses
Path constrained to zero
Head length → Fitness
Max. swim speed → Fitness
Head ht. → Fitness
Tail muscle ht. → Max. swim speed
Tail fin ht. → Max. swim speed
Tail length → Max. swim speed
Tail muscle ht. → Fitness
Head ht. → Max. swim speed
One of the most common hypotheses regarding the evolution of complex life cycles is that alternative morphological strategies are employed to exploit different resources throughout ontogeny. The transition between life stages, however, can be a vulnerable period in which individuals suffer higher rates of mortality. We have presented empirical evidence that tadpoles of the wood frog, Rana sylvatica, when facing predation by dragonfly larvae, experience a higher probability of mortality as they approach metamorphic climax. That mortality probabilities increase during metamorphosis is not unexpected, as a tadpole with emergent hind and forelimbs is well adapted for neither swimming nor jumping [7, 20, 21]. For instance, Arnold and Wassersug  showed across a large geographic range (Mexico to Washington state) that garter snakes, Thamnophis spp., were more likely to have consumed anuran metamorphs (tree frogs and toads) than either tadpoles or adults. They concluded that transforming anurans were highly susceptible to snake predation as a consequence of "locomotor ineptitude."
Our data further suggest that selection acts strongly on morphological traits, favoring tadpoles with short tails and deep tail fins, but that this selection acts largely independently of swimming performance. This latter result is surprising given that tail shape influences swimming performance . Indeed, in our data swimming performance was correlated with both tail length and tail fin depth, and larger values of both tail elements produced greater swimming speeds, consistent with patterns demonstrated elsewhere [56, 57]. Our analyses may have suffered from reduced power given that we could only measure swimming speed for half of our study animals. However, even when we removed all other terms from the model and measured selection differentials on swimming speed alone, the results were not significant. Moreover, path analyses revealed that the effects of morphology (tail length and tail fin height) were largely direct, acting to enhance survival probability per se, rather than serving as a functional link to swimming performance. We interpret this result as consistent with the hypothesis that short tails and deep tail fins are adaptive because they attract predatory strikes and increase the probability that a predator will strike tail tissue rather than sites on the head or body (i.e., "the caudal lure hypothesis"; [25, 27, 30]).
Tadpoles of many frog species exhibit developmental plasticity in response to chemical cues from potential predators, whereby they develop a relatively deep tail fin and a small body (e.g., [22, 58–61]). In particular, enlarged tail fins lead to enhanced survival in the presence of larval odonates (summarized in ). There is reason to believe, however, that differences in tail shape do not influence swimming performance effectively enough to have a large impact on survival in the presence of odonate larvae. This result is unexpected at first blush, given the high prevalence of causal relationships between morphology and performance in other animal systems [62–64]. Van Buskirk and McCollum  used experimental manipulation of tail fin morphology, trimming tissue to reduce both the total length and depth of the tail fin, to investigate the direct effects of changes in tail morphology on swimming performance. Their study revealed that changes in swimming performance were not apparent until one third of the tail was surgically removed, leading them to conclude that reduced susceptibility to predation must have been due to something other than enhanced swimming performance. Similarly, Wilbur and Semlitch  showed that damaged tails of R. utricularia incurred little survivorship cost in the presence of predatory newts (Notophthalmus viridescens). On the other hand, Van Buskirk et al  showed that tadpoles with predator-induced morphologies suffered fewer lethal strikes to the body, suggesting that enlarged tail fins may enhance survival via a "caudal lure" effect.
The approach adopted in this study was to quantify relative survival and selection across ontogeny. One challenge faced by such an approach is that changes in size and shape are confounded throughout the development of the tadpole. This is the phenomenon summarized by Gosner stages. In addition, we were only able to quantify swimming performance and morphometric variables at the start of the study. Depredated tadpoles, unfortunately, cannot be measured. Our analyses thus assume that fundamental elements of size and shape were captured in our initial measures, and that the quantitative signal is maintained to some degree throughout ontogeny. If this were not the case, it is unlikely that we would have obtained sensible results.
Though the number of studies of natural selection has grown rapidly in recent decades [66, 67], there are still fundamental gaps in our understanding of the selective process. This is, in part, owing to the fact that selection studies are rarely replicated either temporally or spatially  and when studies are replicated, selection estimates tend to be highly variable among replicates . Our study provides a rare example of repeatable selection, as replicate estimates of selection were highly congruent among mesocosms, suggesting that the changes that characterize metamorphosis are subject to strong and consistent patterns of selection among individuals.
Our study demonstrates an increase in mortality risk as tadpoles began to metamorphose. Owing to the nature of our experimental design, which focused on tadpole mortality, our data did not examine the effects of the transition from tadpole to froglet on survivorship (see ). As metamorphosis proceeds and the tail fin is resorbed, we expect that froglets would become better at hopping and thus less susceptible to predation. We suggest, as have others , that selection should thus favor individuals that minimize the transition time during metamorphic climax. This does not necessarily mean that selection should favor the most rapid possible development. Indeed, faster overall development often results in small adult body sizes, a condition that can have serious fitness consequences for adult anurans [69, 70]. Rather, the optimal strategy should be to metamorphose at a rate that maximizes the balance between the probability of surviving metamorphosis and later fitness costs. Future studies should aim to measure selection on the separate components of developmental timing to improve our understanding of the targets of selection, including the costs and benefits of pursuing alternative metamorphic strategies.
We thank David Herron, Arielle Filiberti, and Kelsey Pinson for help in the field. Thanks to Robert Cox and two anonymous referees for helpful suggestions that improved the manuscript. All research was conducted under IACUC protocol # 10-02-01 and funding was provided by Dartmouth College. RC was also supported by the National Evolutionary Synthesis Center.
- Ebenman B: Evolution in organisms that change their niches during the life-cycle. American Naturalist. 1992, 139 (5): 990-1021. 10.1086/285370.View ArticleGoogle Scholar
- Wilbur HM: Complex life-cycles. Annual Review of Ecology and Systematics. 1980, 11: 67-93. 10.1146/annurev.es.11.110180.000435.View ArticleGoogle Scholar
- Moran NA: Adaptation and constraint in the complex life-cycles of animals. Annual Review of Ecology and Systematics. 1994, 25: 573-600. 10.1146/annurev.es.25.110194.003041.View ArticleGoogle Scholar
- Werner EE, Hall DJ: Ontogenetic habitat shifts in bluegill: The foraging rate-predation risk trade-off. Ecology. 1988, 69: 1352-1366. 10.2307/1941633.View ArticleGoogle Scholar
- Hanken J, Hall BK, Wake MH: Larvae in amphibian development and evolution. The origin and evolution of larval forms. 1999, 61-108.View ArticleGoogle Scholar
- Istock CA: Evolution of complex life cycle phenomena- An ecological perspective. Evolution. 1967, 21 (3): 592-605. 10.2307/2406619.View ArticleGoogle Scholar
- Wassersug RJ, Sperry DG: Relationship of locomotion to differntial predatin on Pseudacric-tridseriata (Anura Hylidae). Ecology. 1977, 58 (4): 830-839. 10.2307/1936218.View ArticleGoogle Scholar
- Arnold SJ, Wassersug RJ: Differential predation on metamorphic anurans by gartner sankes (Thamnophis)- Social behavior as a possible defense. Ecology. 1978, 59 (5): 1014-1022. 10.2307/1938553.View ArticleGoogle Scholar
- Shaffer HB, Austin CC, Huey RB: The consequences of metamorphosis on salamander (Ambystoma) locomotor performance. Physiological Zoology. 1991, 64 (1): 212-231.View ArticleGoogle Scholar
- Blouin MS: Genetic correlations among morphometric traits and rates of growth and differentiation in the green tree frog, Hyla cinerea. Evolution. 1992, 46 (3): 735-744. 10.2307/2409641.View ArticleGoogle Scholar
- Van Buskirk J, Saxer G: Delayed costs of an induced defense in tadpoles? Morphology, hopping, and development rate at metamorphosis. Evolution. 2001, 55 (4): 821-829. 10.1554/0014-3820(2001)055[0821:DCOAID]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Watkins TB: A quantitative genetic test of adaptive decoupling across metamorphosis for locomotor and life-history traits in the Pacific tree frog, Hyla regilla. Evolution. 2001, 55 (8): 1668-1677.View ArticlePubMedGoogle Scholar
- Azizi E, Landberg T: Effects of metamorphosis on the aquatic escape response of the two-lined salamander (Eurycea bislineata). Journal of Experimental Biology. 2002, 205 (6): 841-849.PubMedGoogle Scholar
- Wilson RS, Kraft PG, Van Damme R: Predator-specific changes in the morphology and swimming performance of larval Rana lessonae. Functional Ecology. 2005, 19 (2): 238-244. 10.1111/j.1365-2435.2005.00958.x.View ArticleGoogle Scholar
- Dejongh HJ: Functional morphology of the jaw aparatus of larval and metamorphosing Rana temporaria. Netherlands Journal of Zoology. 1968, 18 (1): 1-103.View ArticleGoogle Scholar
- Svensson EI, Calsbeek R, (eds.): Evolution on Adaptive Landscapes. 2012, Oxford, UK: Oxford University Press
- Gosner KL: A simplified table for staging anuran embryos and larvae with notes on identification. Herpetologiea. 1960, 16: 183-190.Google Scholar
- McDiarmid RW, Altig R: Tadpoles: The biology of anuran larvae. Tadpoles: The biology of anuran larvae. 1999, i-xiv. 1-444Google Scholar
- Wassersug RJ: Locomotion in amphibian larvae (or why aren't tadpoles built like fishes). American Zoologist. 1989, 29 (1): 65-84.View ArticleGoogle Scholar
- Huey RB: Sprint velocity of tadpoles (Bufo boreas) through metamorphosis. Copeia. 1980, 537-540. 3
- Dudley R, King VA, Wassersug RJ: The implications of shape and metamorphosis for drag forces on a generalized pond tadpole (Rana catesbeiana). Copeia. 1991, 252-257. 1
- Van Buskirk J, Relyea RA: Selection for phenotypic plasticity in Rana sylvatica tadpoles. Biol J Linn Soc. 1998, 65 (3): 301-328.View ArticleGoogle Scholar
- Van Buskirk J, McCollum SA, Werner EE: Natural selection for environmentally induced phenotypes in tadpoles. Evolution. 1997, 51: 1983-1992. 10.2307/2411018.View ArticleGoogle Scholar
- Van Buskirk J, McCollum SA: Influence of tail shape on tadpole swimming performance. Journal of Experimental Biology. 2000, 203 (14): 2149-2158.PubMedGoogle Scholar
- Johnson JB, Burt DB, DeWitt TJ: Form, function, and fitness: Pathways to survival. Evolution. 2008, 62 (5): 1243-1251. 10.1111/j.1558-5646.2008.00343.x.View ArticlePubMedGoogle Scholar
- Teplitsky c, Plenet S, Lena JP, Mermet N, Malet E, Joly P: Escape behaviour and ultimate causes of specific induced defences in an anuran tadpole. Journal of Evolutionary Biology. 2005, 18: 180-190. 10.1111/j.1420-9101.2004.00790.x.View ArticlePubMedGoogle Scholar
- Caldwell JP: Disruptive selection- a tail color polymorphoism in Acris tadpoles in response to differential predation. Canadian Journal of Zoology-Revue Canadienne De Zoologie. 1982, 60 (11): 2818-2827. 10.1139/z82-361.View ArticleGoogle Scholar
- Van Buskirk J, Anderwald P, Lupold S, Reinhardt L, Schuler H: The lure effect, tadpole tail shape, and the target of dragonfly strikes. Journal of Herpetology. 2003, 37 (2): 420-424. 10.1670/0022-1511(2003)037[0420:TLETTS]2.0.CO;2.View ArticleGoogle Scholar
- Van Buskirk J, Aschwanden J, Buckelmuller I, Reolon S, Ruttiman S: Bold tail coloration protects tadpoles from dragonfly strikes. Copeia. 2004, 599-602. 3
- Johnson KE, Eidietis L: Tadpole body zones differ with regard to strike frequencies and kill rates by dragonfly naiads. Copeia. 2005, 909-913. 4
- Semlitsch RD, Boone MD, Dodd CK: Aquatic mesocosms. Amphibian ecology and conservation: a handbook of techniques [Techniques in Ecology & Conservation Series]. 2010, 87-104.Google Scholar
- Lande R, Arnold SJ: The measurement of selection on correlated characters. Evolution. 1983, 37: 1210-1226. 10.2307/2408842.View ArticleGoogle Scholar
- Brodie EDI, Moore AJ, Janzen FJ: Visualizing and quantifying natural selection. Trends Ecol Evol. 1995, 10 (8): 313-318. 10.1016/S0169-5347(00)89117-X.View ArticlePubMedGoogle Scholar
- Akaike H: Factor-Analysis and Aic. Psychometrika. 1987, 52 (3): 317-332. 10.1007/BF02294359.View ArticleGoogle Scholar
- Burnham KP, Anderson DR: Model selection and multimodel inference: a practical information-theoretical approach. 2002Google Scholar
- Phillips PC, Arnold SJ: Visualizing Multivariate Selection. Evolution. 1989, 43 (6): 1209-1222. 10.2307/2409357.View ArticleGoogle Scholar
- Stinchcombe JR, Agrawal AF, Hohenlohe PA, Arnold SJ, Blows MW: Estimating nonlinear selection gradients using quadratic regression coefficients: Double or nothing?. Evolution. 2008, 62 (9): 2435-2440. 10.1111/j.1558-5646.2008.00449.x.View ArticlePubMedGoogle Scholar
- Mitchell-Olds T, Shaw RG: Regression analysis of natural selection: statistical and biological interpretation. Evolution. 1987, 41: 1149-1161. 10.2307/2409084.View ArticleGoogle Scholar
- Brodie ED, Janzen FJ: On the assignment of fitness values in statistical analyses of selection. Evolution. 1996, 50: 437-442. 10.2307/2410813.View ArticleGoogle Scholar
- Janzen FJ, Stern HS: Logistic regression for empirical studies of multivariate selection. Evolution. 1998, 52 (6): 1564-1571. 10.2307/2411330.View ArticleGoogle Scholar
- Svensson EI, Eroukhmanoff F, Friberg M: Effects of natural and sexual selection on adaptive population divergence and premating isolation in a damselfly. Evolution. 2006, 60 (6): 1242-1253.View ArticlePubMedGoogle Scholar
- Brodie ED: Correlational Selection for Color Pattern and Antipredator Behavior in the Garter Snake Thamnophis ordinoides. Evolution. 1992, 46 (5): 1284-1298. 10.2307/2409937.View ArticleGoogle Scholar
- Petraitis PS, Dunham AE, Niewiarowski PH: Inferring multiple causality: The limitations of path analysis. Functional Ecology. 1996, 10 (421-431):
- Schluter D, Nychka D: Exploring fitness surfaces. Amer Nat. 1994, 143 (4): 597-616. 10.1086/285622.View ArticleGoogle Scholar
- Wright S: The method of path coefficients. Ann Math Stats. 1934, 5: 161-215. 10.1214/aoms/1177732676.View ArticleGoogle Scholar
- Li CC: Path Analysis - a primer. 1975, Pacific Grove, California: Boxwood PressGoogle Scholar
- Garland T, Losos JB: Ecological morphology of locomotor performance in squamate reptiles. Ecological Morphology: Integrative Organismal Biology. Edited by: Wainwright PC, Reilly SM. 1994, Chicago: University of Chicago PressGoogle Scholar
- Gomez JM, Zamora R: Spatial variation in the selective scenarios of Hormatophylla spinosa (Cruciferae). American Naturalist. 2000, 155: 630-648.View ArticleGoogle Scholar
- Arbuckle JL: Amos 18 User's Guide. 2009, Chicago, IL: SPSS Inc.Google Scholar
- Gomez JM, Perfectti F, Camacho JPM: Natural selection on Erysimum mediohispanicum flower shape: insights into the evolution of zygomorphy. American Naturalist. 2006, 168: 531-545. 10.1086/507048.View ArticlePubMedGoogle Scholar
- Carstens BC, Stoute HN, Reid NM: An information-theoretic approach to phylogeography. Molecular Ecology. 2009, 18: 4270-4282. 10.1111/j.1365-294X.2009.04327.x.View ArticlePubMedGoogle Scholar
- Tabachnick BG, Fidell LS: Using multivariate statistics. 1996, New York, NY: Harper Collins College, 3Google Scholar
- Arbuckle JL: AMOS - Analysis of Moment Structures. Psychometrika. 1994, 59 (1): 135-137. 10.1007/BF02294272.View ArticleGoogle Scholar
- Gelman A, Carlin JB, Stern HS, Rubin DB: Bayesian Data Analysis. 2004, Boca Raton: Chapman and Hall/CRC, 2Google Scholar
- Dayton GH, Saenz D, Baum KA, Langerhans RB, DeWitt TJ: Body shape, burst speed and escape behavior of larval anurans. Oikos. 2005, 111 (3): 582-591. 10.1111/j.1600-0706.2005.14340.x.View ArticleGoogle Scholar
- Wassersug RJ, Hoff K: The kinematics of swimming in anuran larvae. Journal of Experimental Biology. 1985, 119: 1-30.Google Scholar
- Liu H, Wassersug RJ, Kawachi K: A computational fluid dynamics study of tadpole swimming. Journal of Experimental Biology. 1996, 199 (6): 1245-1260.PubMedGoogle Scholar
- Van Buskirk J: Natural variation in morphology of larval amphibians: Phenotypic plasticity in nature?. Ecological Monographs. 2009, 79 (4): 681-705. 10.1890/08-1692.1.View ArticleGoogle Scholar
- Van Buskirk J, McCollum SA: Functional mechanisms of an inducible defence in tadpoles: morphology and behaviour influence mortality risk from predation. Journal of Evolutionary Biology. 2000, 13 (2): 336-347. 10.1046/j.1420-9101.2000.00173.x.View ArticleGoogle Scholar
- Relyea RA: The many faces of predation: How induction, selection, and thinning combine to alter prey phenotypes. Ecology. 2002, 83 (7): 1953-1964. 10.1890/0012-9658(2002)083[1953:TMFOPH]2.0.CO;2.View ArticleGoogle Scholar
- Relyea RA: Fine-tuned phenotypes: Tadpole plasticity under 16 combinations of predators and competitors. Ecology. 2004, 85 (1): 172-179. 10.1890/03-0169.View ArticleGoogle Scholar
- Arnold SJ: Morphology, performance and fitness. Am Zool. 1983, 23: 347-361.View ArticleGoogle Scholar
- Wainwright PC, Reilly SM: Ecological Morphology: Integrative Organismal Biology. 1999, Chicago, ILL: University of Chicago PressGoogle Scholar
- Schwenk K: Feeding: Form, Function, and Evolution in Tetrapod Vertebrates. 2000, San Diego, CA: Academic PressGoogle Scholar
- Wilbur HM, Semlitsch RD: Ecological consequences of tail injury in Rana tadpoles. Copeia. 1990, 18-24. 1
- Kingsolver JGH, Hoekstra HE, Hoekstra JM, Berrigan D, Vignieri SN, Hill CE, Hoang A, Gilbert P, Beerli P: The strength of phenotypic selection in natural populations. Am Nat. 2001, 157: 245-261. 10.1086/319193.View ArticlePubMedGoogle Scholar
- Hereford J, Hansen TF, Houle D: Comparing strengths of directional selection: how strong is strong?. Evolution. 2004, 58: 2133-2143.View ArticlePubMedGoogle Scholar
- Siepielski AM, DiBattista JD, Carlson SM: It's about time: the temporal dynamics of phenotypic selection in the wild. Ecology Letters. 2009, 12 (11): 1261-1276. 10.1111/j.1461-0248.2009.01381.x.View ArticlePubMedGoogle Scholar
- Merila J, Laurila A, Lindgren B: Variation in the degree and costs of adaptive phenotypic plasticity among Rana temporaria populations. Journal of Evolutionary Biology. 2004, 17 (5): 1132-1140. 10.1111/j.1420-9101.2004.00744.x.View ArticlePubMedGoogle Scholar
- Altwegg R, Reyer HU: Patterns of natural selection on size at metamorphosis in water frogs. Evolution. 2003, 57 (4): 872-882.View ArticlePubMedGoogle Scholar
This article is published under license to 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.