Decades of field data reveal that turtles senesce in the wild
Edited by Joan E. Strassmann, Washington University in St. Louis, St. Louis, MO, and approved March 30, 2016 (received for review January 2, 2016)
Commentary
May 25, 2016
Significance
Turtles are icons of longevity, popularly characterized as lacking aging and remaining robust as they get older. Indeterminate growth and a positive relationship between body size and fecundity suggest that a greater proportion of reproductive output could come from older, rather than younger, individuals. However, studies of turtle populations are typically too short to empirically test these assertions. We tracked >1,000 painted turtles for >20 y in a population in northwest Illinois, United States. Contrary to traditional thought, successful reproduction and survival declined as the turtles aged. Consistent with life-history theory, the observed senescence in reproduction and survival in this population may be attributable to relatively high extrinsic mortality as a result of human disruption.
Abstract
Lifespan and aging rates vary considerably across taxa; thus, understanding the factors that lead to this variation is a primary goal in biology and has ramifications for understanding constraints and flexibility in human aging. Theory predicts that senescence—declining reproduction and increasing mortality with advancing age—evolves when selection against harmful mutations is weaker at old ages relative to young ages or when selection favors pleiotropic alleles with beneficial effects early in life despite late-life costs. However, in many long-lived ectotherms, selection is expected to remain strong at old ages because reproductive output typically increases with age, which may lead to the evolution of slow or even negligible senescence. We show that, contrary to current thinking, both reproduction and survival decline with adult age in the painted turtle, Chrysemys picta, based on data spanning >20 y from a wild population. Older females, despite relatively high reproductive output, produced eggs with reduced hatching success. Additionally, age-specific mark–recapture analyses revealed increasing mortality with advancing adult age. These findings of reproductive and mortality senescence challenge the contention that chelonians do not age and more generally provide evidence of reduced fitness at old ages in nonmammalian species that exhibit long chronological lifespans.
Sign up for PNAS alerts.
Get alerts for new articles, or get an alert when an article is cited.
Why do some organisms show little to no signs of aging as they get older, whereas others exhibit substantial physiological deterioration and reproductive senescence with advancing age (1, 2)? This question has motivated many studies of aging in the wild, ranging from reviews of demographic aging (3, 4) to compilations of mechanistic studies (e.g., ref. 5 and references therein). Senescence should evolve when selection is weaker on deleterious traits expressed at an old age relative to those expressed at a young age (6). Indeed, mutations with senescent effects can persist because of a tradeoff between beneficial effects early in life and pleiotropic detrimental effects late in life (7, 8). Importantly, the persistence of age-specific deleterious mutations will therefore be influenced by levels and sources of extrinsic mortality (e.g., predation, resource scarcity, infectious disease) and the particular ages at which extrinsic sources cause death (9, 10). For example, the onset of physiological or reproductive deterioration is expected to occur relatively early in populations exposed to high levels of extrinsic mortality of adults. Accordingly, such organisms are expected to evolve rapid development/growth, high reproductive effort at young ages, and a shortened lifespan. For species or populations that experience low extrinsic mortality rates, physiological and reproductive function are expected to decline more slowly with advancing age or with a delay in the onset of senescent decline (7). Moreover, such age-specific selection dynamics can result in negligible (11) and even negative (12) senescence.
To clarify the role of natural selection in shaping both reproductive and general senescence and lifespans, an understanding of not only mortality but also reproductive parameters and their fitness consequences across adult ages is required. For instance, although declining fertility in old age indicates reproductive senescence (13), how the resultant age-related changes in reproductive output translate to lifetime fitness is largely unknown. Indeed, even constant reproductive output with age may not necessarily translate into increasing Darwinian fitness with advancing age; reproduction would actually need to increase if selection against deleterious mutations is to remain strong across the lifespan (14). Intuitively, one would expect that relatively low reproductive output (attributable to declining fertility) would result in reduced fitness, but such an effect may not hold if the quality of offspring (and hence the likelihood that offspring contribute to future generations) differs as individuals age. That is, the relative number of progeny produced by old versus young females may not necessarily reflect age differences in reproductive fitness if there is age-related variation in offspring quality and survival. Finally, whether a postreproductive life stage exists in any given species and how such a postreproductive stage affects an individual’s lifetime reproductive fitness are pressing questions in evolutionary biology (15, 16).
Here, we analyze 24 y of individual-based birth, death, and reproductive data gathered for more than 1,000 marked individuals since 1988 from a population of painted turtles (Chrysemys picta) inhabiting the backwaters of the Mississippi River in the state of Illinois, United States (17). We quantify age-related changes in mortality and reproduction from capture–mark–recapture and reproductive/nesting behavior data of known-age individuals. We assessed hatching success of eggs each September (18, 19), which we use as a measure of fitness. With these unique long-term data we ask the following. (i) How does reproductive output and reproductive behavior change with advancing age? (ii) How does fitness change with advancing age? (iii) How does mortality change with advancing age? We then interpret our findings in the context of evolutionary theories for rates of senescence.
Results
Age-Specific Reproduction.
Survival and physical data were collected on 2,234 nests constructed by 600 painted turtles with known identities and reproductive ages (Table S1). Although reproductive and nest microhabitat variables generally increased with age, only egg mass and nest vegetation cover exhibited a statistically significant increase (Fig. 1). Importantly, however, our analyses detected a significant interactive effect of age × plastron length on egg mass (F14,891 = 6.3; P < 0.001) such that the generally positive relationship between female body size and egg size waned with older ages: the slopes declined and the y intercepts increased with advancing age (Fig. 2). These findings imply that investment into egg size is similar between relatively small and large late-age females, but large late-age females invest more than smaller females of all other age classes [sensu (20)]. Thus, the relatively large eggs produced by old females resulted in larger hatchlings than those produced by young females, irrespective of maternal body size (Fig. 1B).
Fig. 1.
Fig. 2.
Fig. S1.
Table S1.
Reproductive age (y) | Excluding depredated nests | Including depredated nests | Nests predated (%) | ||
---|---|---|---|---|---|
Females | Nests | Females | Nests | ||
1 | 226 | 276 | 590 | 718 | 61.6 |
2 | 105 | 151 | 231 | 319 | 52.7 |
3 | 59 | 87 | 168 | 245 | 64.5 |
4 | 55 | 75 | 166 | 233 | 67.8 |
5 | 35 | 51 | 133 | 184 | 72.3 |
6 | 42 | 68 | 100 | 145 | 53.1 |
7 | 40 | 51 | 85 | 119 | 57.1 |
8 | 17 | 35 | 58 | 89 | 60.7 |
9 | 6 | 11 | 39 | 59 | 81.4 |
10 | 5 | 10 | 21 | 31 | 67.7 |
11 | 8 | 15 | 20 | 32 | 53.1 |
12 | 7 | 7 | 13 | 17 | 58.8 |
13 | 1 | 1 | 9 | 15 | 93.3 |
14 | 1 | 2 | 13 | 17 | 88.2 |
15 | 4 | 5 | 7 | 11 | 54.6 |
Total unique females or nests | 355 | 845 | 600 | 2,234 | 62.2 |
Data reported for the no. of females represent the total no. of unique nesting females at each age, and the total in the bottom row represents the no. of unique females used in this study, rather than the sum total down each column. Reproductive ages 1–14 reflect actual ages of 5–7 to 18–20 y old, and reproductive age 15 includes individuals that are ≥20 y old.
Age-Specific Fitness.
Females that produced large eggs and frequent clutches yielded more successfully hatched offspring than individuals that produced small eggs and infrequent clutches (P values, <0.05; Table 1). Fitness (number of eggs hatched per year) declined in relatively old turtles despite a moderate age-related increase in total egg output (number of eggs produced each year; Fig. 3). Moreover, fitness was influenced by the interactions of total egg output and the number of clutches produced at each reproductive age (Fig. 3). Specifically, fitness declined with age, particularly for individuals with relatively high seasonal egg output (i.e., defined as either total egg output or clutch frequency per season); in animals with overall low egg output, an age-related decline in fitness was not evident. Because depredated nests were excluded from these analyses, extrinsic mortality of eggs attributable to nest predation did not drive this age-related variation in hatching success (patterns were generally similar when depredated nests were included in the analyses; Figs. S2–S4).
Table 1.
Effect | Parameter estimate, SE | ndf, ddf | F | P |
---|---|---|---|---|
Excluding depredated nests | ||||
Reproductive age | −0.0520 (0.01) | 1, 166 | 18.5 | <0.001 |
Plastron length | 0.0751 (0.04) | 1, 166 | 3.6 | 0.059 |
No. of clutches produced | −0.2828 (0.13) | 1, 166 | 5.1 | 0.025 |
Egg mass | 0.0913 (0.04) | 1, 166 | 5.3 | 0.023 |
Clutch size | −0.0183 (0.06) | 1, 166 | 0.1 | 0.763 |
Total egg output per season | 0.5483 (0.13) | 1, 166 | 16.7 | <0.001 |
Nest vegetation cover | 0.0611 (0.03) | 1, 166 | 3.5 | 0.063 |
Nest distance to water | −0.0312 (0.03) | 1, 166 | 0.9 | 0.439 |
Age × number of clutches produced | 0.0453 (0.02) | 1, 166 | 6.9 | 0.009 |
Age × total egg output per season | −0.0400 (0.02) | 1, 166 | 6.8 | 0.010 |
Including depredated nests | ||||
Reproductive age | −0.0432 (0.02) | 1, 763 | 6.0 | 0.014 |
Plastron length | 0.0694 (0.07) | 1, 763 | 1.1 | 0.293 |
No. of clutches produced | 0.0500 (0.13) | 1, 763 | 0.2 | 0.691 |
Egg mass | 0.0323 (0.06) | 1, 763 | 0.3 | 0.616 |
Clutch size | −0.1237 (0.08) | 1, 763 | 2.2 | 0.140 |
Total egg output per season | 0.4019 (0.14) | 1, 763 | 8.1 | 0.005 |
Nest vegetation cover | 0.1580 (0.05) | 1, 763 | 10.0 | 0.002 |
Nest distance to water | −0.0097 (0.05) | 1, 763 | 0.0 | 0.852 |
The first model excluded depredated nests, and thus fitness (dependent variable) represented variation in relative egg-hatching success solely attributable to factors other than nest predation. The second model included depredated nests (relative fitness, 0 for depredated nests). Both models included maternal identity as a random effect. Relationships are illustrated in Fig. 3 and Figs. S2–S4. ddf, denominator degrees of freedom; ndf, numerator degrees of freedom. P values in bold face are statistically significant.
Fig. 3.
Fig. S2.
Fig. S3.
Fig. S4.
Age-Specific Mortality.
Mortality analysis that included encounter histories (capture–mark–recapture) for 1,031 females from 1993 to 2012 revealed measurable mortality senescence [i.e., accelerating mortality probability with advancing age (Fig. 4)]. The Gompertz function [ux = Aebx, where ux is the age-specific mortality hazard, A is the initial mortality rate at first nesting (i.e., the first year of the adult stage), and b is the rate of increasing mortality] had the lowest Deviance Information Criterion (DIC) value among those models considered (Table 2). Estimated parameters were A (initial mortality rate/year) = 0.102 (95% credible interval: 0.090–0.116) and b = 0.050 (95% credible interval: 0.034–0.067). This result corresponds to a doubling of mortality every 13.8 y (10.3–20.2 y). Two additional mortality functions are included in Fig. 4 to provide comparison with mammalian-type curves for baboons [ux = 0.028e0.123x (21)] and domestic dogs [ux = 0.02e0.23x (22)].
Fig. 4.
Table 2.
Model | DIC | ΔDIC | pD | k | ||
---|---|---|---|---|---|---|
Function | Makeham? | Temporal? | ||||
Exponential | No | no | 8,198.8 | 833.7 | 292.8 | 2 |
Exponential | No | yes | 8,004.9 | 639.8 | 328.1 | 21 |
Gompertz | No | no | 7,466.0 | 100.9 | 242.2 | 3 |
Gompertz | No | yes | 7,365.1 | 0.0 | 300.3 | 22 |
Gompertz | Yes | no | 7,440.9 | 75.8 | 220.4 | 4 |
Gompertz | Yes | yes | 7,435.6 | 70.5 | 381.6 | 23 |
Logistic | No | no | 7,473.6 | 108.5 | 246.6 | 4 |
Logistic | No | yes | 7,427.5 | 62.4 | 368.8 | 23 |
Logistic | Yes | no | * | * | * | 5 |
Logistic | Yes | yes | * | * | * | 24 |
Weibull | No | no | 8,170.8 | 805.6 | 395.6 | 3 |
Weibull | No | yes | 7,880.2 | 515.1 | 294.7 | 22 |
Weibull | Yes | no | * | * | * | 4 |
Weibull | Yes | yes | * | * | * | 23 |
The models were fitted in a Bayesian framework and used to select among alternative survival models to describe mortality patterns for turtles. Functions included those where senescence occurs (Gompertz, logistic, and Weibull) and the exponential where no senescence occurs. In addition, support for a Makeham parameter and for temporal variation in detection among years was examined. The model with the smallest DIC value had the greatest support. pD is the effective number of parameters; k is the number of fixed parameters in the model. An asterisk indicates that the model did not converge.
Discussion
Survival senescence has generally been thought to be negligible or absent in various reptiles (refs. 23–27) but see ref. 3) because traits that either protect against predators (such as external ribcages and venom; reviewed in ref. 28) or that increase reproductive success with advancing age (such as indeterminate growth and fecundity) may overcome the declining power of natural selection with advancing age that characterizes natural populations (10, 14). In the specific case of C. picta, which has both adult-protective morphology and indeterminate growth and fecundity, lifetime fitness would increase if individuals maintained successful and increasing reproduction through adulthood. Such steady or even increasing reproduction with advancing age is particularly likely for species that exhibit low rates of extrinsic mortality at old ages, indeterminate growth, or a positive relationship between body size and fecundity (e.g., ref. 26) [i.e., many reptiles, including turtles (27)]. Although our aim was to understand whether and how indeterminate growth and fecundity may offset declining late-life selection, our 24-y dataset on C. picta revealed instead a significant decline in both reproductive success and adult survival for females. Moreover, when we tested among competing models for mortality acceleration, some of which would allow deceleration later in life, the Gompertz model had the most support (Table 2), suggesting increasing mortality acceleration over the entire adult lifespan (Fig. 4).
Our results differ from those for another C. picta population, which exhibits both overall low mortality and no differences in survivorship of old- versus middle-aged females (23). Although the survivorship curves of ref. 23 are consistent with gently increasing mortality with advancing age in that population, for direct tests of senescence, it is important to estimate and analyze age-specific mortality hazards as we have done here (29). Thus, estimating the rate and form of mortality senescence by computing instantaneous mortality rates revealed that the present study population of C. picta is characterized by mortality senescence (i.e., accelerating mortality with advancing age). Our study population differs from that of ref. 23 in that ours has experienced a relatively long history of human-driven mortality. The extent to which this human-induced mortality is in any way age-dependent is unknown. Specific human-induced mortality derives primarily from injuries to nesting females (e.g., when crossing roads) and injuries to adults of both sexes in the water column (e.g., by boating activity). Were these injuries to manifest in an age-dependent manner across the lifespan, the age-specific selection gradients could be altered. In contrast, if the increased mortality affected all adult ages equally, we would expect acceleration of the entire life history, because the population age structure would become skewed toward young ages. A comparison of the high mortality population (this study) and a lower mortality population (23) supports this contention. Specifically, here we detected rapid young-adult growth rates (Fig. S5), increased reproductive effort at young adult ages, and reduced lifespan of turtles in our population. Moreover, extensive demographic analyses of prematuration growth in our population indicated faster individual growth and earlier maturation than in less human-impacted populations (30). Thus, these two painted turtle populations fit strikingly well with predictions from life-history theory based on the force and sensitive ages of mortality (14, 31) (Fig. S4). Overall then, these studies support the notion that the life-history strategy is a set of correlated traits that respond predictably to natural selection.
Fig. S5.
Interestingly, in contrast to reproduction and survival, measures related to body size and behavior did not exhibit senescence (Fig. 1 and Fig. S5). This result accords with evolutionary senescence theory in that such traits (i.e., growth and learning, respectively) could buffer selective forces that accelerate death and infertility rates (9, 32). At the same time, hatching success and adult mortality—traits that are substantially influenced by physiological processes (27)—exhibited senescence. Long-term research that not only assesses shifts in survival and reproductive parameters with age but also quantifies age-related changes in physiology and offspring fitness may reveal that senescent patterns in long-lived reptiles are common (see also refs. 33, 34). In addition, long-term studies that explore tradeoffs between early- and late-life performance will provide crucial information for better understanding these declines in reproduction or survival as individuals age. For example, greater resource allocation to early-life growth or reproduction is expected to come at a cost to late-life performance that can lead to age-related declines in reproduction and survival [i.e., disposable soma theory of aging (35)]. These patterns have been observed in many vertebrates (36), including reptiles (37, 38), and could contribute to the reproductive and survival senescence observed here in C. picta.
Determining the universality of reproductive aging has received great recent interest both in terms of documenting its existence and calculating whether reproductive aging occurs at similar (vs. faster) rates than general senescence (13). Analyses have taken the form of either measuring the rate of declining fertility and reproductive success with advancing age (39, 40) or measuring the distribution of ages of last successful reproduction and comparing that to the distribution of deaths (16). Our analysis is similar to the former of these approaches, and we indeed document reproductive decline with advancing age, although not in measures traditionally used (such as rates of conception) but rather by calculating egg-hatching success. Similarly, Congdon et al. (24) reported greater embryo mortality attributable to arrested development in eggs produced by the oldest female age class in Blanding’s turtles. Other long-lived species tend to exhibit reproductive senescence at rates that are faster than somatic aging and that ends at an earlier age than death, which yields a postreproductive lifespan [primates including humans (16), killer whales (41), but see elephants (42, 43)]. These mammalian species may exhibit reproductive senescence attributable to a shelf-life of their primary oocytes, but this possibility is unlikely the case in C. picta because reptiles are not oocyte-limited (44). One firm conclusion is that reproductive senescence in turtles and other ectothermic amniotes cannot be attributable to adaptive explanations invoked in the mammalian literature [e.g., the “grandmother hypothesis,” which largely rests on the role of parental care to increase one’s own fitness or the fitness of relatives (45)].
A critical question is whether senescent changes in mortality (as in primates) or in reproduction drive declines in fitness. In reptiles, which are often characterized as indeterminate growers, few studies document either mortality senescence or the decline in some physiological trait with advancing adult age (3). Because many species with indeterminate growth increase fecundity with increasing age (attributable in part to increasing body size), the question becomes, at what level of increased fecundity will selection cease declining with advancing age? Although we demonstrate mortality senescence in our population of C. picta, the age-specific increase in mortality is considerably less rapid than that reported for mammalian taxa (e.g., Fig. 4), which is consistent with the view that rates of senescence in these long-lived reptiles are relatively slow. These divergent patterns may be explained by very different rates of extrinsic mortality in turtles (or their ancestors) compared with mammals. Indeed, because of their protective morphology (shell) and lack of senescence in immune function (46), turtles likely experience considerably less late-age mortality attributable to extrinsic factors than do mammals. This comparison between turtles and mammals fits with theoretical expectations for patterns of survival senescence in these taxa, but our results challenge the traditional view that long-lived ectotherms, such as turtles, do not exhibit declining reproductive function with advancing age. Indeed, declines in fitness with age will be exacerbated by a combination of increased mortality rate and decreased reproductive success, as we demonstrate in these long-lived ectotherms.
Methods
Since 1988, nesting patterns of painted turtles (C. picta) have been continually monitored at the Thomson Causeway Recreation Area in northwestern Illinois (Carroll County; 41°57′N, 90°7′W) (10–12). Unmarked females exhibiting 5–7 annuli on the pectoral scutes were classified as primiparous (i.e., a reproductive age of one, chronological age 5–7 y old). Thus, all future recaptures of these animals were at a known age (20). Individuals that could not be assigned ages confidently were excluded from analyses of reproduction unless the turtles were recaptured across a 15+ y time span, in which case, data collected from captures ≥15 y after initial capture were included in the oldest age class (i.e., chronological age ≥20 y old). During each nesting season (mid-May to early July), the study area was monitored intensively for nesting turtles. All females at the nesting grounds were individually marked, and their eggs were weighed and counted (typically within 4 h of oviposition). All eggs were returned to their nests and allowed to incubate naturally, and nests were checked for predation almost daily until the end of June of each year. Microhabitat data (shade cover and distance from the water) were also collected for each nest to assess age-related changes in nesting behavior. Nesting behaviors are important because maternal choice of nest microhabitat differs among age classes (47) and can influence egg survival (18). These data are available in the Dryad Digital Repository (dx.doi.org/10.5061/dryad.g2r87).
This research has been approved by the Iowa State University Institutional Animal Care and Use Committee and the Illinois Department of Natural Resources.
Age-Specific Reproduction.
Analyses of age-specific fitness included 600 individuals with known ages captured from 1997 to 2010. Of these individuals, 65% were observed nesting more than once, and most were first-time nesters (Table S1). Thus, our comparisons across ages were largely based on longitudinally sampled individuals. Linear mixed models were used to evaluate age-specific variation in reproductive variables (i.e., egg mass, clutch size, clutch frequency, total egg output per female per season, hatchling mass) and nest microhabitat variables (i.e., nest shade cover, distance from nest site to nearest water). In cases where females nested two or three times per season, within-season mean female values were used for egg mass, clutch size, nest shade cover, and nest distance to water. Each dependent variable was standardized to a mean of 0 and unit variance within each year before analysis. Reproductive age (1–15 y) was defined as a fixed main effect, and plastron length was a covariate. The interaction term of plastron length × reproductive age was removed from models when it was not significant, which occurred in all cases except in the analysis of egg mass. Maternal identity was included as a random effect to account for repeated measurements on females across ages. Linear and second-order polynomial regressions were used to detect age-specific trends in reproductive and nest microhabitat variables.
Age-specific Fitness.
We used generalized linear mixed models to test for linear selection on each reproductive variable and interactions with reproductive age using egg-hatching success per female as a measure of fitness, standardized within season [i.e., ωrelative = ω/ωmean (14)]. This measure of fitness (number of successfully hatched eggs per clutch) strongly correlates with the likelihood of offspring recruiting into future adult age classes (17) and is therefore a meaningful proxy for maternal reproductive success (i.e., fitness). Independent variables (standardized to a mean of 0 with unit variance) included plastron length, egg mass, clutch size, clutch frequency, nest shade cover, distance from the nest site to nearest water, and total egg output per female per season, as well as all two-way interactions with reproductive age. Three-way and higher-order interactions were not included in the model because they were not relevant to our specific questions, and models with higher-level interactions often would not converge, because degrees of freedom were reduced by including these factors. Nonsignificant interaction terms were sequentially removed from the initial model, and final models were selected based on the lowest Akaike Information Criterion (AIC) score. Again, maternal identity was included as a random effect to account for repeated measurements on females across ages. Parameter estimates for fixed effects were calculated using restricted maximum likelihood and statistical significance was determined with F tests.
The selection analyses described above were performed twice. First, to assess non–predation-related variation in fitness, all depredated nests were excluded from our analyses. Second, to assess how nest predation impacts the relationship between the reproductive parameters and fitness (i.e., selection), all depredated nests (depredated nests had a fitness of 0) were included in the second set of analyses. Standardizations of all traits (including fitness) were recalculated accordingly for the subsets of nests used in each analysis. Results from analyses that included versus excluded depredated nests were compared qualitatively. Because the potential fate of eggs in depredated nests (had they been left intact) is unknown, this qualitative comparison assumes that nest predation is random with respect to the variables measured. A generalized linear mixed model using reproductive traits (egg mass, clutch size, clutch number) and nest microhabitat (vegetation cover, nest distance to water) as independent variables and maternal identity as a random effect showed that this assumption was met (Table S2).
Table S2.
Trait | ndf | ddf | F | P |
---|---|---|---|---|
Egg mass | 1 | 1,876 | 0.01 | 0.925 |
Clutch size | 1 | 1,876 | 0.14 | 0.707 |
Clutch no. | 1 | 1,876 | 2.54 | 0.112 |
Nest vegetation cover | 1 | 1,876 | 0.37 | 0.545 |
Nest distance to water | 1 | 1,876 | 3.02 | 0.083 |
Analyses were performed with a generalized linear mixed model including maternal identity as a random variable. ddf, denominator degrees of freedom; ndf, numerator degrees of freedom.
Additional analyses were performed to confirm and visualize significant effects. Regression analyses (weighted by sample size because of repeated measures per female) were used to assess the relationship between fitness and reproductive and nest variables. Second-order polynomial regressions were used to assess age-specific trends in fitness. Lastly, to visualize two-way interactions between age and reproductive output on fitness, selection surfaces were calculated with a first-degree polynomial locally weighted scatterplot smoothing technique (for total egg output) and a negative exponential smoothing technique (for number of clutches produced) in Sigma Plot.
Age-Specific Mortality.
Analyses of age-specific mortality included 1,031 individuals. These individuals included the 600 used in the analyses above but also included turtles with unknown ages that were captured from 1993 to 2012. To test for a signature of demographic senescence in mortality, accelerating mortality models were fit to mark–recapture data for both known age and unknown age females using a Bayesian hierarchical analysis (48, 49). DIC was used to choose among alternative models for the survival function fit using the BaSTA (Bayesian Survival Trajectory Analysis) package (48) (version 1.3) in R (version 2.14; R development Core Team 2011) (Table 2). Four survival functions were considered: exponential, Gompertz, logistic, and Weibull functions (49, 50). In the exponential model, survival is constant across ages, whereas the other three models all allow for decreasing survival as individuals become older. Thus, the alternative models permitted assessment of whether there was evidence for demographic senescence and, if so, what form the senescence took. In addition, models that included an age-independent constant to the survival function (Makeham coefficient) were considered for the Gompertz, logistic, and Weibull models. The approach does not allow for temporal or cohort effects to be incorporated. However, previous work detected no evidence for an overall temporal trend in survival across the study period (51) or for cohort effects on survival related to developmental conditions in the population (38). Finally, alternatives that differed in whether detection was constant or variable among all years were examined.
Data Availability
Data deposition: The reproductive data and encounter histories reported in this paper have been deposited in the Dryad Digital Repository (https://doi.org/10.5061/dryad.g2r87).
Acknowledgments
We thank the many volunteers, students, and postdoctoral researchers who have been involved with fieldwork at Turtle Camp over the past decades. We thank members of F.J.J.’s laboratory and T. Schwartz for helpful comments on earlier drafts of this paper and J. Sherwood for assistance with ArcGIS software. We acknowledge ongoing support from the US Fish and Wildlife Service and the US Army Corp of Engineers. Funding for this research was provided by National Science Foundation Grants DEB-9629529, DEB-0089680, DEB-0640932, and DEB-1242510 (to F.J.J.). Additional support was provided by a National Institutes of Health Grant R01AG049416.
Supporting Information
Supporting Information (PDF)
Supporting Information
- Download
- 856.66 KB
References
1
SN Austad, Methusaleh’s Zoo: How nature provides us with clues for extending human health span. J Comp Pathol 142, S10–S21 (2010).
2
DJ Holmes, MA Ottinger, RE Ricklefs, CE Finch, Sosa-2: Introduction to the proceedings of the second symposium on organisms with slow aging. Exp Gerontol 38, 721–722 (2003).
3
DH Nussey, H Froy, J-F Lemaitre, J-M Gaillard, SN Austad, Senescence in natural populations of animals: Widespread evidence and its implications for bio-gerontology. Ageing Res Rev 12, 214–225 (2013).
4
DA Roach, JR Carey, Population biology of aging in the wild. Annu Rev Ecol Evol Syst 45, 421–443 (2014).
5
QE Fletcher, C Selman, Aging in the wild: Insights from free-living and non-model organisms. Exp Gerontol 71, 1–3 (2015).
6
B Charlesworth, Fisher, Medawar, Hamilton and the evolution of aging. Genetics 156, 927–931 (2000).
7
WD Hamilton, The moulding of senescence by natural selection. J Theor Biol 12, 12–45 (1966).
8
GC Williams, Pleiotropy, natural selection, and the evolution of senescence. Evolution 11, 398–411 (1957).
9
PD Williams, T Day, Antagonistic pleiotropy, mortality source interactions, and the evolutionary theory of senescence. Evolution 57, 1478–1488 (2003).
10
A Baudisch, Hamilton’s indicators of the force of selection. Proc Natl Acad Sci USA 102, 8263–8268 (2005).
11
CE Finch, Update on slow aging and neglibible senescence - a mini-review. Gerontology 55, 307–313 (2009).
12
JW Vaupel, A Baudisch, M Dölling, DA Roach, J Gampe, The case for negative senescence. Theor Popul Biol 65, 339–351 (2004).
13
TBL Kirkwood, DP Shanley, The connections between general and reproductive senescence and the evolutionary basis of menopause. Ann N Y Acad Sci 1204, 21–29 (2010).
14
BC Charlesworth Evolution in Age-Structured Populations (Cambridge Univ Press, 2nd Ed, Cambridge, UK, 1994).
15
AA Cohen, Female post-reproductive lifespan: A general mammalian trait. Biol Rev Camb Philos Soc 79, 733–750 (2004).
16
SC Alberts, et al., Reproductive aging patterns in primates reveal that humans are distinct. Proc Natl Acad Sci USA 110, 13440–13445 (2013).
17
LE Schwanz, RJ Spencer, RM Bowden, FJ Janzen, Climate and predation dominate juvenile and adult recruitment in a turtle with temperature-dependent sex determination. Ecology 91, 3016–3026 (2010).
18
DA Warner, CF Jorgensen, FJ Janzen, Maternal and abiotic effects on egg mortality and hatchling size of turtles: Temporal variation in selection over seven years. Funct Ecol 24, 857–866 (2010).
19
FJ Janzen, DA Warner, Parent-offspring conflict and selection on egg size in turtles. J Evol Biol 22, 2222–2230 (2009).
20
RM Bowden, HK Harms, RT Paitz, FJ Janzen, Does optimal egg size vary with demographic stage because of a physiological constraint? Funct Ecol 18, 522–529 (2004).
21
AM Bronikowski, et al., Aging in the natural world: Comparative data reveal similar mortality patterns across primates. Science 331, 1325–1328 (2011).
22
CE Finch, MC Pike, M Witten, Slow mortality rate accelerations during aging in some animals approximate that of humans. Science 249, 902–905 (1990).
23
JD Congdon, et al., Testing hypotheses of aging in long-lived painted turtles (Chrysemys picta). Exp Gerontol 38, 765–772 (2003).
24
JD Congdon, RD Nagle, OM Kinney, RC van Loben Sels, Hypotheses of aging in a long-lived vertebrate, Blanding’s turtle (Emydoidea blandingii). Exp Gerontol 36, 813–827 (2001).
25
OR Jones, et al., Diversity of ageing across the tree of life. Nature 505, 169–173 (2014).
26
M Olsson, R Shine, Does reproductive success increase with age or with size in species with indeterminate growth? A case study using sand lizards (Lacerta agilis). Oecologia 105, 175–178 (1996).
27
TS Schwartz, AM Bronikowski, Molecular stress pathways and the evolution of life histories in reptiles. Molecular Mechanisms of Life History Evolution: The Genetics and Physiology of Life History Traits and Trade-Offs, eds T Flatt, A Heyland (Oxford Univ Press, Oxford), pp. 193–209 (2011).
28
MA Blanco, PW Sherman, Maximum longevities of chemically protected and non-protected fishes, reptiles, and amphibians support evolutionary hypotheses of aging. Mech Ageing Dev 126, 794–803 (2005).
29
AM Bronikowski, T Flatt, Aging and its demographic measurement. Nat Educ Knowl 3, 3 (2010).
30
RJ Spencer, FJ Janzen, Demographic consequences of adaptive growth and the ramifications for conservation of long-lived organisms. Biol Conserv 143, 1951–1959 (2010).
31
SC Stearns, M Ackermann, M Doebeli, M Kaiser, Experimental evolution of aging, growth, and reproduction in fruitflies. Proc Natl Acad Sci USA 97, 3309–3313 (2000).
32
AM Bronikowski, DEL Promislow, Testing evolutionary theories of aging in wild populations. Trends Ecol Evol 20, 271–273 (2005).
33
M Massot, et al., An integrative study of ageing in a wild population of common lizards. Funct Ecol 25, 848–858 (2011).
34
DH Nussey, T Coulson, M Festa-Bianchet, JM Gaillard, Measuring senescence in wild animal populations: Towards a longitudinal approach. Funct Ecol 22, 393–406 (2008).
35
TBL Kirkwood, MR Rose, Evolution of senescence: Late survival sacrificed for reproduction. Philos Trans R Soc Lond B Biol Sci 332, 15–24 (1991).
36
JF Lemaitre, et al., Early-late life trade-offs and the evolution of ageing in the wild. Proc R Soc Lond B 282, 20150209 (2015).
37
M Massot, P Aragón, Phenotypic resonance from a single meal in an insectivorous lizard. Curr Biol 23, 1320–1323 (2013).
38
DAW Miller, FJ Janzen, GM Fellers, PM Kleeman, A Bronikowski, Biodemography of ectothermic tetrapods provides insights into the evolution and plasticity of mortality trajectories. Sociality, Hierarchy, Health: Comparative Demography Advances in Biodemography: Cross-Species Comparisons of Social Environments and Social Behaviors, and Their Effects on Health and Longevity, eds M Weinstein, MA Lane (Natl Acad Press, Washington, DC), pp. 295–314 (2014).
39
C Packer, M Tatar, A Collins, Reproductive cessation in female mammals. Nature 392, 807–811 (1998).
40
AM Sparkman, SJ Arnold, AM Bronikowski, An empirical test of evolutionary theories for reproductive senescence and reproductive effort in the garter snake Thamnophis elegans. Proc Biol Sci 274, 943–950 (2007).
41
EJ Ward, K Parsons, EE Holmes, 3rd KC Balcomb, JKB Ford, The role of menopause and reproductive senescence in a long-lived social mammal. Front Zool 6, 4 (2009).
42
CJ Moss, The demography of an African elephant (Loxodonta africana) population in Amboseli, Kenya. J Zool (Lond) 255, 145–156 (2001).
43
MR Robinson, KU Mar, V Lummaa, Senescence and age-specific trade-offs between reproduction and survival in female Asian elephants. Ecol Lett 15, 260–266 (2012).
44
SM Jones, Hormonal regulation of ovarian function in reptiles. Hormones and Reproduction of Vertebrates, eds DO Norris, KH Lopez (Academic, Amsterdam) Vol 3, 89–115 (2011).
45
K Hawkes, JF O’Connell, NGB Jones, H Alvarez, EL Charnov, Grandmothering, menopause, and the evolution of human life histories. Proc Natl Acad Sci USA 95, 1336–1339 (1998).
46
L Schwanz, DA Warner, S McGaugh, R Di Terlizzi, A Bronikowski, State-dependent physiological maintenance in a long-lived ectotherm, the painted turtle (Chrysemys picta). J Exp Biol 214, 88–97 (2011).
47
HK Harms, RT Paitz, RM Bowden, FJ Janzen, Age and season impact resource allocation to eggs and nesting behavior in the painted turtle. Physiol Biochem Zool 78, 996–1004 (2005).
48
F Colchero, OR Jones, M Rebke, BaSTA: An R package for Bayesian estimation of age-specific survival from incomplete mark-recapture/recovery data with covariates. Methods Ecol Evol 3, 466–470 (2012).
49
F Colchero, JS Clark, Bayesian inference on age-specific survival for censored and truncated data. J Anim Ecol 81, 139–149 (2012).
50
SD Pletcher, Model fitting and hypothesis testing for age-specific mortality data. J Evol Biol 12, 430–439 (1999).
51
A Jergenson, DAW Miller, LA Neuman-Lee, DA Warner, FJ Janzen, Swimming against the tide: Resilience of a riverine turtle to extreme environment events. Biol Lett 10, 20130782 (2014).
52
AM Bronikowski, et al., Data from: Aging in the natural world: Comparative data reveal similar mortality patterns across primates. Dryad Digital Repository. Available at dx.doi.org/10.5061/dryad.8682. (2011).
53
A Comfort, The Biology of Senescence (Elsevier, New York), 3rd Ed. (1979).
54
JW Gibbons, Population structure and survivorship in painted turtle Chrysemys picta. Copeia 1968, 260–268 (1968).
55
JW Gibbons, Variation in growth rates in three populations of the painted turtle, Chrysemys picta. Herpetologica 23, 296–303 (1967).
56
HM Wilbur, A growth model for the painted turtle Chrysemys picta. Copeia 1975, 337–343 (1975).
Information & Authors
Information
Published in
Classifications
Data Availability
Data deposition: The reproductive data and encounter histories reported in this paper have been deposited in the Dryad Digital Repository (https://doi.org/10.5061/dryad.g2r87).
Submission history
Published online: May 2, 2016
Published in issue: June 7, 2016
Keywords
Acknowledgments
We thank the many volunteers, students, and postdoctoral researchers who have been involved with fieldwork at Turtle Camp over the past decades. We thank members of F.J.J.’s laboratory and T. Schwartz for helpful comments on earlier drafts of this paper and J. Sherwood for assistance with ArcGIS software. We acknowledge ongoing support from the US Fish and Wildlife Service and the US Army Corp of Engineers. Funding for this research was provided by National Science Foundation Grants DEB-9629529, DEB-0089680, DEB-0640932, and DEB-1242510 (to F.J.J.). Additional support was provided by a National Institutes of Health Grant R01AG049416.
Notes
This article is a PNAS Direct Submission.
See Commentary on page 6328.
Authors
Competing Interests
The authors declare no conflict of interest.
Metrics & Citations
Metrics
Citation statements
Altmetrics
Citations
Cite this article
113 (23) 6502-6507,
Export the article citation data by selecting a format from the list below and clicking Export.
Cited by
Loading...
View Options
View options
PDF format
Download this article as a PDF file
DOWNLOAD PDFLogin options
Check if you have access through your login credentials or your institution to get full access on this article.
Personal login Institutional LoginRecommend to a librarian
Recommend PNAS to a LibrarianPurchase options
Purchase this article to access the full text.