- Research article
- Open Access
Stress-related changes in leukocyte profiles and telomere shortening in the shortest-lived tetrapod, Furcifer labordi
BMC Evolutionary Biology volume 20, Article number: 160 (2020)
Life history theory predicts that during the lifespan of an organism, resources are allocated to either growth, somatic maintenance or reproduction. Resource allocation trade-offs determine the evolution and ecology of different life history strategies and define an organisms’ position along a fast–slow continuum in interspecific comparisons. Labord’s chameleon (Furcifer labordi) from the seasonal dry forests of Madagascar is the tetrapod species with the shortest reported lifespan (4–9 months). Previous investigations revealed that their lifespan is to some degree dependent on environmental factors, such as the amount of rainfall and the length of the vegetation period. However, the intrinsic mechanisms shaping such a fast life history remain unknown. Environmental stressors are known to increase the secretion of glucocorticoids in other vertebrates, which, in turn, can shorten telomeres via oxidative stress. To investigate to what extent age-related changes in these molecular and cellular mechanisms contribute to the relatively short lifetime of F. labordi, we assessed the effects of stressors indirectly via leukocyte profiles (H/L ratio) and quantified relative telomere length from blood samples in a wild population in Kirindy Forest. We compared our findings with the sympatric, but longer-lived sister species F. cf. nicosiai, which exhibit the same annual timing of reproductive events, and with wild-caught F. labordi that were singly housed under ambient conditions.
We found that H/L ratios were consistently higher in wild F. labordi compared to F. cf. nicosiai. Moreover, F. labordi already exhibited relatively short telomeres during the mating season when they were 3–4 months old, and telomeres further shortened during their post-reproductive lives. At the beginning of their active season, telomere length was relatively longer in F. cf. nicosiai, but undergoing rapid shortening towards the southern winter, when both species gradually die off. Captive F. labordi showed comparatively longer lifespans and lower H/L ratios than their wild counterparts.
We suggest that environmental stress and the corresponding accelerated telomere attrition have profound effects on the lifespan of F. labordi in the wild, and identify physiological mechanisms potentially driving their relatively early senescence and mortality.
Life history theory is based on the premise that during the lifetime of an individual, energy and resources are allocated to either growth, somatic maintenance, or reproduction [1,2,3]. For example, resources, which are invested into fast growth and early reproduction, cannot be used for somatic maintenance, which may lead to shorter lifespans. The trade-offs between traits shape life history strategies and the distribution of species along a fast–slow continuum of life history speeds [4,5,6]. In spite of the supposed significance of extrinsic factors in shaping life histories, aging research is still largely biased towards captive animals living under standardized, optimal conditions (e.g., ). In the wild, studies of senescence have largely focused on long-lived animals that face relatively low levels of extrinsic mortality (e.g. sea turtles , birds , Soay sheep ). However, studies focusing on age-related changes in short-lived species in the wild are rare. Hence, studies of wild populations with high extrinsic mortality are essential for testing hypotheses on the evolution of lifespan and senescence.
Oxidative stress and its damage to macromolecules is one of the most cited causes of aging . The oxidative damage is a byproduct of aerobic respiration  and intensified by chronic stress conditions characterized by a persistent release of glucocorticoids (GCs) in vertebrates . Physiological stress is an important mediator in the trade-off between survival and reproduction [14, 15]. GCs are released in response to a wide range of stressful stimuli (e.g., ), and several of their effects parallel those observed during aging, suggesting that chronic stress has a potential to accelerate the aging process [17, 18].
The immunosuppressive effects of chronic GC elevation and their consequences for morbidity and mortality have been studied intensively [19, 20]. Alterations in key immunological parameters during chronic stress parallel those during normal immunosenescence to a large degree . These hormones are important regulators of carbohydrate, lipid, and protein metabolism , and several earlier studies linked poor body condition to elevated GC concentrations (e.g., ). The direct measurement of baseline GC levels in wildlife via blood plasma can be challenging as stress hormones can rise immediately following capture . However, leukocyte profiles are a suitable tool to indirectly assess stress levels as these hormones increase the number of heterophils and decrease the number of lymphocytes. Leukocyte responses to stress take about 12 h to several days in ectotherms (reviewed in ). Heterophils are the primary phagocytic leukocyte, which proliferate in circulation in response to infections, inflammation and stress [26,27,28,29,30]. Lymphocytes are involved in a variety of immunological functions such as the production of immunoglobulin and modulation of immune defense .
At the cellular level, telomere length (TL) and shortening are thought to be significant proximate contributors to the aging process. Telomeres are short, tandem-repeated sequences of DNA found at the ends of linear eukaryotic chromosomes, whose sequence (TTAGGG) is highly conserved among vertebrates . Telomeres function in stabilizing chromosomal end integrity , inhibiting aberrant fusions and rearrangements that occur on broken chromosomes , and aiding in the completion of duplication . During each cell cycle, telomeric repeats are lost because DNA polymerase is unable to completely replicate the 3′end of linear DNA .
There is great variation among species in age-specific TL . Sexual differences in TL and attrition have been suggested to contribute to sex-specific disease and mortality patterns in humans [37, 38], where women typically have longer telomeres and are longer-lived (e.g., ). Telomerase, the enzyme that countervails telomere shortening was found to be active in stem cells, gametes and most cancer cells, but normally absent from or at very low levels in most somatic cells . However, some studies in reptiles suggested that telomerase may not be turned off in adult somatic cells . Besides cell division dependent telomere shortening, elevated levels of corticosterone can further affect TL via increased oxidative damage by reactive oxygen species (ROS) [42, 43]. Elevated GCs, particularly during long-term physiological or psychological stress, have been linked to increased oxidative stress and concomitant telomere shortening and reduced telomerase activity [43, 44]. As the nucleobase guanine is a major oxidation target for ROS, the (TTAGGG) repeats are particularly exposed to oxidative damage .
Telomeres may also act as sentinels of the general level of DNA damage in a given cell. High levels of telomere damage would be indicative of high levels of damage to the coding sequences. Thus, telomeres could offer a mechanism to ensure that cells with high levels of DNA damage soon terminate division . Overall, demanding life history stages and harsh environmental conditions seem to be linked to a rapid rate of telomere degradation, and there is also a clear connection between physiological stress and telomere attrition in humans, laboratory rodents and wild vertebrates [44, 47,48,49,50]. This evidence suggests that telomere dynamics could be closely related to stress in wild vertebrates (reviewed in ), and Houben et al.  emphasized that telomeres are a promising biomarker for chronic oxidative stress.
Labord’s chameleon (Furcifer labordi) from the seasonal deciduous dry forests in western and southwestern Madagascar has a lifespan of only 4–9 months [53, 54]. This extreme life history makes this species an interesting model for studying potential mechanisms of accelerated senescence, especially because longer-lived sympatric congeners are available for comparative studies. During their short lives, these chameleons hatch at the beginning of the wet season in November, passing through subsequent fast juvenile growth, maturation and courtship, followed by the death of both sexes during the early dry season in May [53, 54]. Wild females tend to live slightly longer, whereas no sex difference in lifespan was found in caged individuals kept under ambient conditions . Fast growth rates, high reproductive rates and intense mating competition might proximately contribute to increased stress levels and telomere shortening, which in turn may facilitate the decrease of physiological functioning, ultimately leading to death (e.g., [55, 56]).
To investigate whether the ratio of heterophils and lymphocytes (H/L ratio) and telomere shortening are associated with the early die-off in F. labordi in the wild, we determined their telomere dynamics as well as their leukocyte profiles as an indict measure of physiological stress. Our study included two comparisons; one between wild F. labordi and their sympatric and longer-lived congener F. cf. nicosiai, and one with F. labordi kept in single cages under ambient conditions, shielding them substantially from environmental stressors, like hunger or predation risk. We predicted an increase in H/L ratios as well as rapidly shortening telomeres in post-reproductive wild F. labordi as well as lower H/L ratios and decelerated telomere attrition in F. cf. nicosiai. Furthermore, as age-related changes should be delayed in the longer-lived females of both species, we predicted females to exhibit comparatively slower rate of senescence than males. Finally, caged F. labordi, which were shielded from extrinsic mortality and from a substantial part of the costs of reproduction and starvation, were expected to exhibit slower correlates of aging compared to their wild conspecifics.
In both species, heterophils were the most abundant leukocyte type, followed by lymphocytes azurophils and basophils. Heterophiles exhibited a spherical shape with an eccentric mostly lobed nucleus containing clumpy basophilic purplish chromatin. Most lymphocytes contained a large nucleus with coarse chromatin, leaving only a small visible band of cytoplasma around it. Basophils were only found sporadically. On average, the H/L ratio of F. labordi (2.45 ± 0.97 SD, n = 319) was significantly higher compared to that of F. cf. nicosiai (1.51 ± 0.47 SD, n = 103, t = − 9.921, p < 0.001). Moreover, we detected an increase of the H/L ratio in both species between February and May (Fig. 1, Table 1), reflecting the cessation of mating activities. In captive specimen, we found an average H/L profile of (1.42 ± 0.14 SD, n = 40) and no significant sex differences (Table 2). As in their wild conspecifics, the H/L ratio of captive chameleons increased significantly from February until June (Table 2).
During our sampling period, we did not detect any significant sex and age-related changes in TL in F. labordi. Average TL was significantly longer in F. cf. nicosiai (t = 6.438, p < 0.001). Furthermore, TL of F. cf. nicosiai was comparatively long in March (1.87 ± 0.77 SD, n = 14) and decreased dramatically until May (1.14 ± 0.33 SD, n = 10, t = − 2.686, p < 0.01). Moreover, TL of F. cf. nicosiai males was significantly shorter compared to females (t = − 2.67, p < 0.01, df = 38). For statistical analyses (Table 3), the months June and July were excluded due to small sample sizes (but June is included in Fig. 2), and we found a negative correlation between the H/L ratio and TL in F. labordi (\(r\)= − 0.556, df = 65, p < 0.01) and in F. cf. nicosiai (\(r\) = − 0.687, df = 38, p \(<0.01\); see Fig. 3).
Our study revealed that H/L ratios were consistently higher in wild F. labordi compared to F. cf. nicosiai, hinting at higher stress levels in the shorter-lived species. Furcifer labordi already exhibited relatively short telomeres when they were 3–4 months old. TL was initially comparatively longer in F. cf. nicosiai, but undergoing rapid shortening after the mating season. In this species, we also detected intersexual differences in H/L ratio and TL, with shorter living males exhibiting higher H/L ratios and shorter telomeres. Interestingly, heterophils were the most common leucocyte type in both wild and captive chameleons. Captive F. labordi exhibited comparatively longer lifespans and lower H/L profiles than their wild conspecifics. In planning this study, we assumed that the captive chameleons would be buffered from some environmental stressors, like starvation, desiccation and predation risk. Our data therefore indicate that relatively long-lived wild F. labordi individuals were, on average, more stressed and lived shorter lives than their captive conspecifics, indicating a link between stress and longevity.
Baseline stress levels and leukocyte profiles
Investigations in other reptile species indicated large differences between hematology values of different species as well as intraspecific variation as a function of season and sex [57, 58]. In their study of blood chemistry and hematology in captive panther chameleons (Furcifer pardalis), Laube et al.  found that lymphocytes were the predominant leukocyte type in both summer and winter. In contrast, Cuadrado et al.  reported that heterophils were the most frequently found leucocyte type in dystoic and healthy post-reproductive females of the common chameleon (Chamaeleo chamaeleon). The H/L ratio from that study (2.24) resembled the values reported here for F. labordi (2.45). More recently, Eshar et al.  found that heterophils were the most abundant leukocytes type in wild common chameleons. As part of their study of leukocyte profiles of an iguanid species, Davis et al.  reviewed several studies of white blood cell profiles of iguanids and other lizard species. They extracted data on the relative numbers of all cell types (mean percentages) and categorized the studies based on whether lizards were from captivity or the wild. They showed that all wild animals had higher H/L ratios than the captive conspecifics. In fact, the relative abundance of lymphocytes and heterophils was completely opposite in both groups, with lymphocytes being the most abundant leukocyte type in captive lizards and heterophils being most common one in wild specimens. Thus, either wild lizards naturally have higher baseline stress levels (and thus higher H/L ratios) than captive ones, or trapping of wild animals induced stress-related alterations in the animals’ leukocyte profiles, a notion also supported by the elevated H/L ratios of the captive F. labordi in our study.
During a stress response, GC secretion increases partly to mobilize more metabolic energy to deal with the stressor. While this stress response provides obvious short-term benefits, chronic elevation of GCs is harmful [19, 63,64,65]. In the present study, we observed stress-related changes in leukocyte profiles in both chameleon species, which may indirectly contribute to their rapid senescence after the reproductive season. Captive F. labordi showed comparatively lower, but in relation to other captive lizards, elevated H/L ratios , indicating that they perceived these captive conditions as mildly stressful, but that they were also buffered from major environmental stressors. It is possible that the brief biweekly handling to obtain blood samples might have contributed to the perceived stress level of caged individuals, but this manipulation did most likely not impact the measurements of H/L ratios because such effects were found only after 12 h in other species .
Any interpretation of the potential physiological effects of variable H/L ratios should take into account that a review published after our field work found inconsistent relationships between GC profiles and leukocyte profiles across studies ). In gopher tortoises, Gopherus polyphemus, both GC levels and leucocyte profiles changed across seasons, but the changes were not correlated . Moreover, in two studies of garter snakes, Thamnophis sirtalis, conducted by the same research group, but on different populations and in different years, one study revealed a positive correlation between GC levels and H/L ratio , whereas the other did not . Furthermore, the interpretation of leukocyte dynamics relies on baseline data for the taxon of interest . Reports about leukocyte profiles in chameleons in the wild [60, 61] and captivity  are rare and based on relatively small sample sizes. Our study therefore contributes valuable comparative data based on large samples of two wild chameleon species, but future studies may want to assess stress levels more directly, e.g. by measuring GC levels from fecal samples.
Telomere dynamics differed between the two chameleon species. Telomeres were relatively longer in F. cf. nicosiai, but shortened rapidly with the disappearance of the adult cohort. In contrast, the telomeres of F. labordi were relatively short, but a deterioration over time was not detectable. The first 3 months in the life of F. labordi are characterized by fast growth rates, whereas juvenile F. cf. nicosiai show much slower growth and reach maturity at an age of 11–12 months . The lifespan of F. cf. nicosiai is longer, but both species mate at the same time and die off afterwards. A study of wild jackdaws (Corvus monedula) revealed that long telomeres shorten more rapidly than short ones, regardless of the individual’s age . Additionally, telomere degradation was highest in humans with long telomeres . These studies suggest that mechanisms for telomere maintenance exist in vivo, which potentially protect the shortest telomeres from further attrition and might explain why we could not detect any significant TL reduction in F. labordi. It would therefore seem interesting to also examine telomerase activity in these species. In ectothermic vertebrates, the expression of telomerase is frequently found in somatic tissues and is thought to be due to the indeterminate growth . Thus, regulation by this enzyme might enable F. labordi to maintain its TL up to a certain level.
Whether TL is a universal predictor of longevity is still up for debate. Whittemore et al.  found that the telomere shortening rate, but not the initial telomere length alone, is a powerful predictor of life span in several bird and mammal species. These results support the notion that critical telomere shortening and the consequent onset of telomeric DNA damage and cellular senescence are a general determinant of species life span. In humans, telomere attrition is also more rapid in the first decade of life, stabilizes in adulthood and is followed by a gradual loss at old age . We could not study telomere dynamics because of low recapture rates and a lack of data on juveniles, but a relatively large male juvenile F. cf. nicosiai was sampled at approx. 4 months of age and showed a TL of 3.44, which was the highest measured in this species. In contrast, TL of hatchling pythons (Lisais fuscus) was significantly shorter than that of older snakes, increasing during their first year of life and subsequently decreasing with age . Similar curvilinear telomere dynamics were found in frilled-necked lizards (Chlamydosaurus kingii) .
In F. cf. nicosiai, we also observed sexual dimorphism in telomere length across the sampling period, with females having longer average telomeres. The associated longer female survival may be adaptive as the maturation of eggs after insemination takes several weeks, and female chameleons are capable of producing additional clutches from stored sperm (, FE pers. observation). In several other species, including sand lizards (Lacerta agilis) , Medaka fish ) and humans , females also live longer and have longer telomeres. The actual mechanisms contributing to sex-specific telomere patterns are unknown, however. Previous work on humans suggested that the difference in TL stems from larger body mass in men compared to women , leading to the assumption that larger bodies require more tissue growth and cell division. However, female sand lizards are larger than males  and have longer telomeres. Gopalakrishnan et al.  postulated that estrogen is a key factor contributing to the decelerated telomere shortening in female Medaka fish, but corresponding data from other species are lacking. Thus, telomere attrition probably depends on multiple factors that remain to be identified.
Nowadays, telomere attrition is widely recognized as one of the hallmarks of aging (e.g. ), and telomeric assessments are widely used in evolutionary biology as biomarkers of somatic integrity. However, limited attention has been paid to addressing the fundamental question raised by these relationships: Which role do telomeres play in shaping the evolution of life history trade-offs and senescence ? While it is broadly accepted that telomere degradation can have causal effects on cell fates, the extent to which it contributes to age-related declines on organismal level is less clear. A proximate causal role for telomeres would more possibly reflect an adaptive strategy, born out of telomere maintenance costs and/or a function for telomere attrition (e.g. in counteracting cancer), the relative importance of which is currently unclear. Nevertheless, it is frequently mentioned that telomere length as a predictor of overall health could instead reflect it acting as a non-causal biomarker of accumulated damage to other biological structures that themselves have causal deleterious effects on the organismal performance (e.g. ). While it is mechanistically conceivable that telomere dynamics are one proximate cause of current–future trade-offs and senescence, whether telomeres play a significant proximate causal role relative to alternative mechanisms, such as oxidative damage to other biological structures, is currently uncertain . Finally, advances in understanding of the selection pressures that might have shaped a proximate causal role for telomeres according to life history trade-offs have the potential to shed light on the nature of the evolutionary restrictions at play in life history evolution and help explain the form of the current–future trade-offs and ageing trajectories .
Stress-related leukocyte profiles and telomere shortening
In both species, we found a negative correlation between average H/L ratio and TL. Chronical stress has potentially negative consequence through an increase in oxidative damage [42, 43] and ultimately telomere shortening . Oxidative stress also dramatically decreases telomerase activity [86, 87]. Therefore, oxidative stress not only accelerates telomere shortening by direct damage to telomeres, but also by inhibiting telomere restoration as well. Even though we are well aware of the correlational nature of our study, we suggest that physiological stress negatively affected TL in our two study species. Although our findings and additional studies suggest a strong association between stress and telomere shortening [88, 89], we cannot discard other mechanisms that could affect TL, like alterations of early growth rates (e.g. ). More direct future studies should acknowledge that the link between stress response and telomere degradation is probably not straightforward and depends on the benefits and costs of activating an emergency life history state that is species- and context-dependent.
At an ultimate level, rates of extrinsic mortality are thought to determine where a species falls on the slow-fast continuum, with high rates of extrinsic mortality selecting for fast life histories . The results of our previous capture-mark-recapture study  also suggest that extrinsic mortality rates in both chameleon species are presumably high in adults. Williams also postulated that juvenile mortality has no influence on the evolution of senescence; predicting that senescence should be associated with extrinsic mortality rates . However, formal, mathematical theory [92,93,94] showed that this particular prediction is wrong. Accordingly, selection leading to senescence does not directly depend on survival to old age, but rather on the shape of the stable age distribution. The aim of evolutionary theories of aging is to clarify why organismal fitness mechanisms decline with age. Moorad et al.  therefore proposed to investigate the actual phenomenon of aging, not its proxies. More theory and careful physiological measurements from many species under many different environmental conditions are therefore required to further illuminate factors that shape life histories. Remarkably, Reznick and colleagues  even found that guppies (Poecilia reticulata) derived from natural populations with high levels of predation live the longest in the laboratory. This study demonstrates that our understanding of the evolution of senescence will profit from modeling numerous aging parameters, traits other than age at death as well as the causes of mortality.
Although there are many examples of negative correlations between lifespan and the apparent extrinsic risk of death faced by organisms, this risk is more often deduced than measured. In our study species, besides extrinsic mortality at old age, several factors might impact the short lifespan of this species. High juvenile mortality in F. labordi might lead to the extremely high investment in reproduction that in turn facilitates the pronounced stress response and relatively short telomeres. As physiological stress also has a strong influence on immune responses , the increasing gastrointestinal—and blood parasite burden observed in both species in the wild towards the dry season  might reflect an unavoidable consequence of this adaptation. This notion about the physiological processes contributing to such a short life span in F. labordi is also supported by a maximum lifespan in caged individuals of 16 months, indicating that their lifetime is indeed bounded by molecular and cellular mechanisms of aging.
The results of our study provide rare information about leukocyte profiles and telomere dynamics in relation to senescence and mortality patterns of two chameleon species in the wild. The results of this study suggest that the presumably energetically demanding reproductive season in the short-lived species contributes to environmental stress ensued by increased oxidative damage and subsequent accelerated telomere shortening. To fully understand telomere dynamics and their relation to stress-related measures (H/L profiles) in these species, repeated samples from wild specimens and samples from younger life stages are necessary, however.
Study site and study species
This study was conducted in Kirindy Forest, which is located in the region of Menabe Central, Western Madagascar, ca. 60 km northeast of Morondava (44°39′E, 20°03′S, 30–60 m asl). It is one of the largest remaining Malagasy dry deciduous forest fragments. The local climate is characterized by a hot wet season between November and April, followed by a cool dry season from May to October . Kirindy Forest is located near the northern end of the range of Furcifer labordi, a medium-sized and sexually highly dimorphic chameleon from the western and southwestern regions of Madagascar . Furcifer cf. nicosiai is a relatively larger species, also sexually dimorphic , and appears to be limited to intact dry forests [100, 101].
Wild chameleons were located at night using LED flashlights. The capture location was marked and GPS data were taken. We sampled alternating along two transects of 3 km length each. Animals were transported to the nearby research station in cloth bags and handled the following morning. They were sexed, age categorized (hatchling, juvenile, adult), and their snout vent length (SVL) was measured. Animals were individually marked by visual implant elastomers (VIE; Northwest Marine Technology Inc., Shaw Island, WA) . Hatchlings and small juveniles were individually marked with nail polish on the toes. All chameleons were released at their point of capture within 12 h. Sampling took place over three field seasons: November 19, 2013–July 8, 2014, and January 11, 2015–July 15, 2015, and October 12, 2015–December 17, 2015.
We collected a total of 20 male and 20 female juveniles of F. labordi in early January, at approximately two months of age, and kept them individually without visual contact in cylindrical nylon cages (90 cm height, 60 cm diameter) inside a large outdoor forest enclosure. Chameleons received a standardized amount of insects, adjusted to their age and size to match growth and final size of the wild population. Specifically, animals were fed five times per week with two grasshoppers, crickets or butterflies. Water was offered daily with a spray flask. Between February and June, the caged animals were handled biweekly to obtain blood samples. Only captive animals were sampled repeatedly.
To measure leukocyte profiles, a drop of blood was taken by lateral puncture of the caudal vein. No blood was taken from females that were obviously gravid. The drop of blood was placed on a microscope slide and distributed as a blood smear. After air-drying, blood smears were processed with a rapid differential haematology staining utilizing the Diff-Quik staining solution system (Medion Diagnostics AG, Düdingen, Switzerland). We determined the ratio of heterophils to lymphocytes (H/L ratio) in at least 200 cells per slide, using the 100 × oil objective. For identification, the general description of reptilian blood cells and terminology was used [103, 104]. Counting of the leukocytes started at the most distal edge of the feather end of the smear and proceeded one field of view at a time, across the entire smear in an ‘S’ fashion. Only fields of view with > 15 erythrocytes in a monolayer were considered . All cell counting was conducted by FE. Blood samples were taken from February onwards until mid-July in 2014 and 2015. In total, 319 samples from wild F. labordi, 103 samples of F. cf. nicosiai and 278 samples from 40 captive animals were analysed.
Telomere length (TL)
Blood samples for telomere analysis were taken between March and mid-July in 2015. In total, 66 blood samples of F. labordi and 39 of F. cf. nicosiai were obtained. Captive animals were excluded from this analysis because not enough blood samples were available.
For the determination of TL, we took approx. 5–10 µl blood from the caudal vein after lateral puncture with a capillary and transferred it into a 1.5 ml tube containing 0.5 ml SET buffer. Samples were directly frozen at − 20 °C. To avoid melting during transportation, samples were stored in a compressor cooling box. Telomere length was measured using real‐time quantitative PCR (qPCR) using SensiMix SYBR No‐ROX Kit (Bioline, Sydney, NSW, Australia) and a Rotor‐gene 6000 thermocycler (Qiagen, Chadstone, VIC, Australia) according to published protocols [106, 107] using techniques developed by  with the 18S ribosomal RNA (18S) gene as the non-variable copy number reference gene. The telomere primers used were Telb1 (5′‐CGGTTTGTTTGGGTTTGGGTTTGGGTTTGGGTTTGGGTT‐3′) and Telb2 (5′‐GGCTTGCCTTACCCTTACCCTTACCCTTACCCTTACCCT‐3′, (109)). The 18S gene (92 bp amplicon in Anolis) was selected as the reference gene as it had previously been validated in reptiles [8, 107, 109]. The primer sequences used were 18S‐F (5′‐GAGGTGAAATTCTTGGACCGG‐3′) and 18S‐R (5′‐CGAACCTCCGACTTTCGTTCT‐3′). Reactions were run in triplicate for each sample, with each run containing either Telb or 18S primers. Amplifications were carried out in a Rotor‐Gene 6000 thermocycler (Qiagen, Australia) using an initial Taq activation step at 95 °C for 10 min and a total of 40 cycles of 95 °C for 15 s, 60 °C for 15 s, and 72 °C for 15 s. Each reaction had a final volume of 20 μl with 10 ng of DNA. A melt curve was generated after each run over the temperature range of 60 to 95 °C to ensure that there was no nonspecific product amplification.
All of the DNA samples for a given individual were included in the same run. No‐template control reactions were run in triplicate for each primer set during every qPCR run to ensure that there was no contamination.
Linear models (LM) were used to examine the influence of leukocyte profiles on TL in wild F. labordi and F. cf. nicosiai. As fixed factors, we added month (age), sex and species. For captive F. labordi, we used linear mixed models (LMM). As fixed factors, we added month (age), and sex, while ID was included as a random factor for recaptured samples. For all models, we compared the respective full model with the null model by using a likelihood ratio test. In addition, we visually inspected normality and homoscedasticity with residual plots. For model analysis, we used the package lme4 . All data analysis was conducted in R (R-Code Team 2017) . To check for correlation between H/L profile and TL, we calculated the Pearson correlation coefficient.
Availability of data and materials
All data generated or analyzed during this study are included in this published article [and its supplementary information files].
- H/L ratio:
Reactive oxygen species
Stearns SC. Life-history tactics: a review of the ideas. Q Rev Biol. 1976;51:3–47.
Stearns SC. The evolution of life histories. Oxford Univ. 1992.
Stearns SC. Life history evolution: successes, limitations and prospects. Naturwissenschaften. 2000;87:476–86.
Sæther BE, Bakke Ø. Avian life history variation and contribution of demographic traits to the population growth rate. Ecology. 2000;81(3):642–53.
Sæther BE, Engen S, Pape Møller A, Weimerskirch H, Visser ME, Fiedler W, et al. Life-history variation predicts the effects of demographic stochasticity on avian population dynamics. Am Nat. 2004;164(6):793–802.
Salguero-Gómez R, Jones OR, Jongejans E, Blomberg SP, Hodgson DJ, Mbeau-Ache C, et al. Fast–slow continuum and reproductive strategies structure plant life-history variation worldwide. PNAS. 2016;113(1):230–5.
Languille S, Blanc S, Blin O, Canale CI, Dal-Pan A, Devau G, et al. The grey mouse lemur: a non-human primate model for ageing studies. Aging Res Rev. 2012;11(1):150–62.
Plot V, Criscuolo F, Zahn S, Georges JY. Telomeres, age and reproduction in a long-lived reptile. PloS One. 2012;7(7):e40855.
Barrett EL, Burke TA, Hammers M, Komdeur J, Richardson DS. Telomere length and dynamics predict mortality in a wild longitudinal study. Mol Ecol. 2013;22(1):249–59.
Fairlie J, Holland R, Pilkington JG, Pemberton JM, Harrington L, Nussey DH. Lifelong leukocyte telomere dynamics and survival in a free-living mammal. Aging Cell. 2016;15(1):140–8.
Sohal RS, Weindruch R. Oxidative stress, caloric restriction, and aging. Science. 1996;273(5271):59–63.
Sena LA, Chandel NS. Physiological roles of mitochondrial reactive oxygen species. Mol Cell. 2012;48(2):158–67.
Salmon AB, Richardson A, Pérez VI. Update on the oxidative stress theory of aging: does oxidative stress play a role in aging or healthy aging? Free Radic Bio Med. 2010;48(5):642–55.
Boonstra A, Barrat FJ, Crain C, Heath VL, Savelkoul HF, O’Garra A. 1α, 25-Dihydroxyvitamin D3 has a direct effect on naive CD4+ T cells to enhance the development of Th2 cells. J Immunol. 2001;167(9):4974–80.
Ricklefs RE, Wikelski M. The physiology/life-history nexus. Trends Ecol Evol. 2002;17(10):462–8.
Romero LM, Reed JM, Wingfield JC. Effects of weather on corticosterone responses in wild free-living passerine birds. Gen Comp Endocrinol. 2000;118(1):113–22.
Sapolsky R, Rivier C, Yamamoto G, Plotsky P, Vale W. Interleukin-1 stimulates the secretion of hypothalamic corticotropin-releasing factor. Science. 1987;238(4826):522–4.
Veldhuis JD, Sharma A, Roelfsema F. Age-dependent and gender-dependent regulation of hypothalamic-adrenocorticotropic-adrenal axis. Endocrin Metab Clin. 2013;42(2):201–25.
Sapolsky RM, Romero LM, Munck AU. How do glucocorticoids influence stress responses? Integrating permissive, suppressive, stimulatory, and preparative actions. Endocr Rev. 2000;21(1):55–89.
Dhabhar FS. A hassle a day may keep the doctor away: stress and the augmentation of immune function. Integr Comp Biol. 2002;42:556–64.
Bauer ME. Chronic stress and immunosenescence: a review. Neuroimmunomodulat. 2008;15:241–50.
Dallman MF. Stress update: adaptation of the hypothalamic-pituitary-adrenal axis to chronic stress. Trend Endocrin Met. 1993;4(2):62–9.
Wingfield JC, Romero LM. The endocrine system. In: McEwen BS, Goodman HM, editors. Handbook of physiology. 2001. p. 211–234.
Romero LM, Reed JM. Collecting baseline corticosterone samples in the field: is under 3 min good enough? Comp Biochem Phys A. 2005;140(1):73–9.
Davis AK, Maney DL, Maerz JC. The use of leukocyte profiles to measure stress in vertebrates: a review for ecologists. Funct Ecol. 2008;22(5):760–72.
Jain NC. Essentials of veterinary hematology. Philadelphia: Blackwell Publishing; 1993.
Campbell TW. Avian Hematology and Cytology. Ames, Iowa: Iowa State University Press; 1995.
Rupley AE. Manual of avian practice. Bangkok: Univ Bangkok Press; 1997.
Harmon BG. Avian heterophils in inflammation and disease resistance. Poultry Sci. 1998;77:972–7.
Thrall MA. Hematology of amphibians, veterinary hematology and clinical chemistry: text and clinical case presentations. Philadelphia: Lippincott Williams & Wilkins; 2004.
Campbell TW. Clinical pathology. In: Mader DR, editor. Reptile medicine and surgery. W.B. Saunders Company: Philadelphia; 1996. p. 248–57.
Meyne J, Ratliff RL, MoYzIs RK. Conservation of the human telomere sequence (TTAGGG) n among vertebrates. PNAS. 1989;86(18):7049–53.
Prowse KR, Greider CW. Developmental and tissue-specific regulation of mouse telomerase and telomere length. PNAS. 1995;92(11):4818–22.
McClintock B. The stability of broken ends of chromosomes in Zea mays. Genetics. 1941;26(2):234.
Watson JD. Origin of concatemeric T7DNA. Nat New Biol. 1972;239(94):197–201.
Gomes NM, Shay JW, Wright WE. Telomere biology in metazoa. FEBS Lett. 2010a;584(17):3741–51.
Stindl R. Is telomere erosion a mechanism of species extinction? J Exp Zool Part B. 2004;302(2):111–20.
Eskes T, Haanen C. Why do women live longer than men? Eur J Obstet Gynecol Reprod Biol. 2007;133(2):126–33.
Benetos A, Okuda K, Lajemi M, Kimura M, Thomas F, Skurnick J, et al. Telomere length as an indicator of biological aging: the gender effect and relation with pulse pressure and pulse wave velocity. Hypertension. 2001;37(2):381–5.
Tanaka H, Horikawa I, Barrett JC, Oshimura M. Evidence for inactivation of distinct telomerase repressor genes in different types of human cancers. Int J Cancer. 2005;115(4):653–7.
Dantzer B, Fletcher QE. Telomeres shorten more slowly in slow-aging wild animals than in fast-aging ones. Exp Gerontol. 2015;71:38–47.
Agostinho P, Cunha RA, Oliveira C. Neuroinflammation, oxidative stress and the patho-genesis of Alzheimer’s disease. Curr Pharm Des. 2010;16:2766–78.
Constantini D, Marasco V, Møller AP. A meta-analysis of glucocorticoids as modulators of oxidative stress in vertebrates. J Comp Physiol B. 2011;181:447–56.
Epel ES, Blackburn EH, Lin J, Dhabhar FS, Adler NE, Morrow JD, Cawthon RM. Accelerated telomere shortening in response to life stress. PNAS. 2004;101(49):17312–5.
Wang Z, Rhee DB, Lu J, Bohr CT, Zhou F, Vallabhaneni H, et al. Characterization of oxidative guanine damage and repair in mammalian telomeres. PLoS Genet. 2010;6(5):e1000951.
Von Zglinicki T. Replicative senescence and the art of counting. Exp Gerontol. 2003;38(11–12):1259–64.
Kotrschal A, Ilmonen P, Penn DJ. Stress impacts telomere dynamics. Biol Lett. 2007;3(2):128–30.
Haussmann MF, Marchetto NM. Telomeres: linking stress and survival, ecology and evolution. Curr Zool. 2010;56(6):714–27.
Shalev I, Moffitt TE, Sugden K, Williams B, Houts RM, Danese A, et al. Exposure to violence during childhood is associated with telomere erosion from 5 to 10 years of age: a longitudinal study. Mol Psychiatr. 2013;18(5):576–81.
Bateson M. Cumulative stress in research animals: telomere attrition as a biomarker in a welfare context? BioEssays. 2016;38(2):201–12.
Angelier F, Costantini D, Blevin P, Chastel O. Do glucocorticoids mediate the link between environmental conditions and telomere dynamics in wild vertebrates? A review. Gen Comp Endocr. 2018;256:99–111.
Houben JM, Moonen HJ, van Schooten FJ, Hageman GJ. Telomere length assessment: biomarker of chronic oxidative stress? Free Radic Biol Med. 2008;44(3):235–46.
Karsten KB, Andriamandimbiarisoa LN, Fox SF, Raxworthy CJ. A unique life history among tetrapods: an annual chameleon living mostly as an egg. PNAS. 2008;105(26):8980–4.
Eckhardt F, Kappeler PM, Kraus C. Highly variable lifespan in an annual reptile, Labord’s chameleon (Furcifer labordi). Sci Rep. 2017;7(1):11397.
Braithwaite RW, Lee AK. A mammalian example of semelparity. Am Nat. 1979;113(1):151–5.
Oakwood M, Bradley AJ, Cockburn A. Semelparity in a large marsupial. Proc R Soc Lond B Biol. 2001;268(1465):407–11.
Campbell W. Ellis Avian and exotic animal hematology and cytology. 3rd ed. Ames: Blackwell Publishing; 2007. p. 287.
Strik NI, Alleman R, Harr KE. Circulating inflammatory cells. In: Jacobson E, editor. Infectious diseases and pathology of reptiles, color atlas and text. Boca Raton: CRC Press; 2007. p. 167–218.
Laube A, Pendl H, Clauss M, Altherr B, Hatt JM. Plasma biochemistry and hematology reference values of captive panther chameleons (Furcifer pardalis) with special emphasis on seasonality and gender differences. J Zoo Wildl Med. 2016;47(3):743–53.
Cuadrado M, Díaz-Paniagua C, Quevedo MA, Aguilar JM, Prescott IM. Hematology and clinical chemistry in dystocic and healthy post-reproductive female chameleons. J Wildl Dis. 2002;38(2):395–401.
Eshar D, Ammersbach M, Shacham B, Katzir G, Beaufrère H. Venous blood gases, plasma biochemistry, and hematology of wild-caught common chameleons (Chamaeleo chamaeleon). Can J Vet Res. 2018;82(2):106–14.
Davis AK, Ruyle LE, Maerz JC. Effect of trapping method on leukocyte profiles of black-chested spiny-tailed iguanas (Ctenosaura melanosterna): implications for zoologists in the field. ISRN Zool. 2011. https://doi.org/10.5402/2011/384825.
McEwen BS. Protection and damage from acute and chronic stress: allostasis and allostatic overload and relevance to the pathophysiology of psychiatric disorders. Ann NY Acad Sci. 2004;1032(1):1–7.
Romero LM, Wikelski M. Corticosterone levels predict survival probabilities of Galapagos marine iguanas during El Nino events. PNAS. 2001;98(13):7366–70.
Romero LM, Wikelski M. Severe effects of low-level oil contamination on wildlife predicted by the corticosterone-stress response: preliminary data and a research agenda. Spill Sci Technol B. 2002;7(5–6):309–13.
Davis AK, Maney DL. The use of glucocorticoid hormones or leucocyte profiles to measure stress in vertebrates: what’s the difference? Meth Ecol Evol. 2018;9(6):1556–68.
Goessling JM, Guyer C, Mendona MT. Seasonal acclimation of constitutive immunity in gopher tortoises Gopherus polyphemus. Physiol Biochem Zool. 2016;89:487–97.
Sparkman AM, Bronikowski AM, Williams S, Parsai S, Manhart W, Palacios MG. Physiological indices of stress in wild and captive garter snakes: correlations, repeatability, and ecological variation. Comp Biochem Phys A. 2014;174:11–7.
Gangloff EJ, Sparkman AM, Holden KG, Corwin CJ, Topf M, Bronikowski AM. Geographic variation and within-individual correlations of physiological stress markers in a widespread reptile, the common garter snake (Thamnophis sirtalis). Comp Biochem Phys A. 2017;205:68–76.
Eckhardt F, Kraus C, Kappeler PM. Life histories, demographies and population dynamics of three sympatric chameleon species (Furcifer spp from western Madagascar. Amphibia-Reptilia. 2019;40(1):41–54.
Salomons HM, Mulder GV, van de Zande L, Haussmann MF, Linskens MH, Verhulst S. Telomere shortening and survival in free-living corvids. Proc R Soc B: Biol Sci. 2009;276(1670):3157–65.
Nordfjäll K, Svenson U, Norrback KF, Adolfsson R, Lenner P, Roos G. The individual blood cell telomere attrition rate is telomere length dependent. PLoS Genet. 2009;5(2):e1000375.
Gomes NM, Shay JW, Wright WE. Telomeres and telomerase. In: Wolf NS, editor. The comparative biology of aging. New York: Springer; 2010b. p. 227–58.
Whittemore K, Vera E, Martínez-Nevado E, Sanpera C, Blasco MA. Telomere shortening rate predicts species life span. PNAS. 2019;116(30):25122–15127.
Frenck RW Jr, Blackburn EH, Shannon KM. The rate of telomere sequence loss in human leukocytes varies with age. PNAS. 1998;95(10):5607–10.
Ujvari B, Madsen T. Short telomeres in hatchling snakes: erythrocyte telomere dynamics and longevity in tropical pythons. PloS One. 2009;4(10):e7493.
Ujvari B, Biro PA, Charters JE, Brown G, Heasman K, Beckmann C, Madsen T. Curvilinear telomere length dynamics in a squamate reptile. Funct Ecol. 2017;31(3):753–9.
Tolley KA, Chauke LF, Jackson JC, Feldheim KA. Multiple paternity and sperm storage in the Cape dwarf chameleon (Bradypodion pumilum). Afr J Herpetol. 2014;63(1):47–56.
Olsson M, Pauliny A, Wapstra E, Uller T, Schwartz T, Miller E, Blomqvist D. Sexual differences in telomere selection in the wild. Mol Ecol. 2011;20:2085–99.
Gopalakrishnan S, Cheung NK, Yip BW, Au DW. Medaka fish exhibits longevity gender gap, a natural drop in estrogen and telomere shortening during aging: a unique model for studying sex-dependent longevity. Front Zool. 2013;10(1):78.
Seluanov A, Chen Z, Hine C, Sasahara TH, Ribeiro AA, Catania KC, et al. Telomerase activity coevolves with body mass not lifespan. Aging Cell. 2007;6(1):45–52.
Gullberg A, Olsson M, Tegelström H. Male mating success, reproductive success and multiple paternity in a natural population of sand lizards: behavioural and molecular genetics data. Mol Ecol. 1997;6(2):105–12.
Lopez-Otin C, Blasco MA, Partridge L, Serrano M, Kroemer G. The hallmarks of aging. Cell. 2013;153:1194–217.
Young AJ. The role of telomeres in the mechanisms and evolution of life-history trade-offs and ageing. Philos Trans R Soc B. 2018;373(1741):20160452.
Simons MJP. Questioning causal involvement of telomeres in aging. Ageing Res Rev. 2015;24:191–6.
Borrás C, Esteve JM, Viña JR, Sastre J, Viña J, Pallardó FV. Glutathione regulates telomerase activity in 3T3 fibroblasts. J Biol Chem. 2004;279(33):34332–5.
Kurz DJ, Decary S, Hong Y, Trivier E, Akhmedov A, Erusalimsky JD. Chronic oxidative stress compromises telomere integrity and accelerates the onset of senescence in human endothelial cells. J Cell Sci. 2004;117(11):2417–26.
Haussmann MF, Longenecker AS, Marchetto NM, Juliano SA, Bowden RM. Embryonic exposure to corticosterone modifies the juvenile stress response, oxidative stress and telomere length. Proc R Soc B: Biol Sci. 2011;279(1732):1447–56.
Quirici V, Guerrero CJ, Krause JS, Wingfield JC, Vásquez RA. The relationship of telomere length to baseline corticosterone levels in nestlings of an altricial passerine bird in natural populations. Front Zool. 2016;13(1):1.
Schultner J, Moe B, Chastel O, Bech C, Kitaysky AS. Migration and stress during reproduction govern telomere dynamics in a seabird. Biol Lett. 2014;10:20130889.
Williams GC. Pleiotropy, natural selection, and the evolution of senescence. Evolution. 1957;11:398–411.
Abrams PA. Does increased mortality favor the evolution of more rapid senescence? Evolution. 1993;47(3):877–87.
Caswell H. Extrinsic mortality and the evolution of senescence. Trend Ecol Evol. 2007;4(22):173–4.
Wensink MJ, Caswell H, Baudisch A. The rarity of survival to old age does not drive the evolution of senescence. Evol Biol. 2017;44(1):5–10.
Moorad J, Promislow D, Silvertown J. Evolutionary ecology of senescence and a reassessment of Williams’ ‘Extrinsic Mortality’ Hypothesis. Trend Ecol Evol. 2019;34(6):519–30.
Reznick DN, Bryant MJ, Roff D, Ghalambor CK, Ghalambor DE. Effect of extrinsic mortality on the evolution of senescence in guppies. Nature. 2004;431:1095–9.
Padgett DA, Glaser R. How stress influences the immune response. Trend Immunol. 2003;24(8):444–8.
Eckhardt F, Strube C, Mathes KA, Mutschmann F, Thiesler H, Kraus C, Kappeler PM. Parasite burden in a short-lived chameleon, Furcifer labordi. Int J Parasitol PAR. 2019;10:231–40.
Kappeler PM, Fichtel C. A 15-year perspective on the social organization and life history of sifaka in Kirindy Forest. In: Long-term field studies of primates. Berlin: Springer. 2012. p. 101–121
Glaw F, Vences M. A field guide to the amphibians and reptiles of Madagascar. Vences & Glaw. 2007.
Jesu R, Mattioli FSG. On the discovery of a new large chameleon inhabiting the limestone outcrops of western Madagascar: Furcifer nicosiai sp. nov. (Reptilia, Chamaeleonidae). Doriana. 1999;7(311):1–14.
MacNeil JE, Dharmarajan GUHA, Williams RN. Salamarker: a code generator and standardized marking system for use with visible implant elastomers. Herpetol Conserv Bio. 2011;6(2):260–5.
Origgi FC. Reptile immunology. In: Jacobson E, editor. Infectious diseases and pathology of reptiles, color atlas and text. Boca Raton: CRC Press; 2007. p. 131–66.
Tamukai K, Takami Y, Akabane Y, Kanazawa Y, Une Y. Plasma biochemical reference values in clinically healthy captive bearded dragons (Pogona vitticeps) and the effects of sex and season. Vet Clin Path. 2011;40:368–73.
Davis AK, Maerz JC. Comparison of hematological stress indicators in recently captured and captive paedomorphic mole salamanders, Ambystoma talpoideum. Copeia. 2008;2008(3):613–7.
Rollings N, Friesen CR, Sudyka J, Whittington C, Giraudeau M, Wilson M, Olsson M. Telomere dynamics in a lizard with morph-specific reproductive investment and self-maintenance. Ecol Evol. 2017;7:5163–9.
Rollings N, Uhrig EJ, Krohmer RW, Waye HL, Mason RT, Olsson M, et al. Age-related sex differences in body condition and telomere dynamics of red-sided garter snakes. Proc R Soc Lond Ser B Biol Sci. 2017;284(9):20162146.
Cawthon RM. Telomere measurement by quantitative PCR. Nucleic Acids Res. 2002;30:e47–e47.
Rollings N, Friesen CR, Whittington CM, Johansson R, Shine R, Olsson M. Sex- and tissue-specific differences in telomere length in a reptile. Ecol Evol. 2019;9:6211–19.
Team RC. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org. 2017.
We thank the CAFF/CORE of the Direction Général des Eaux et Forêts, and CNFEREF Morondava for authorization and support for this study. We are thankful for support in the field to Camilla Cenni, Hauke Thiesler, Wanda von Bremen and the team of field assistants in Kirindy Forest. We thank the DGHT Wilhelm-Peters Fonds for funding parts of the telomere analysis at the University of Gothenburg in Sweden.
Open Access funding enabled and organized by Projekt DEAL. The DGHT Wilhelm-Peters Fonds funded parts of the telomere analyses at the University of Gothenburg in Sweden. This funding body played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
All work conducted in Madagascar was done with the written approval by the Commission Ad hoc Faune et Flore/Commité d’Orientation sur la Recherche Environnementale (CAFF/CORE) of the Direction Général des Eaux et Forêts of the Malagasy Ministry of the Environment and Sustainable Development (research permission No. 053/13 issued on 21 Fenruary 2013) and the Centre National de Formation, d’Etudes et de Recherche en Environnement et Forestier (CNFEREF) Morondava, Madagascar.
Consent for publication
The authors declare that they have no competing interests. Author Mats Olsson is an Associate Editor of this journal.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Eckhardt, F., Pauliny, A., Rollings, N. et al. Stress-related changes in leukocyte profiles and telomere shortening in the shortest-lived tetrapod, Furcifer labordi. BMC Evol Biol 20, 160 (2020). https://doi.org/10.1186/s12862-020-01724-2
- Furcifer labordi
- Life history
- H/L ratio
- Body condition