# The influence of phylogeny and life history on telomere lengths and telomere rate of change among bird species: A meta-analysis

François Criscuolo, F. Stephen Dobson, Quentin Schull

► To cite this version:

François Criscuolo, F. Stephen Dobson, Quentin Schull. The influence of phylogeny and life history on telomere lengths and telomere rate of change among bird species: A meta-analysis. Ecology and Evolution, 2021, 10.1002/ece3.7931. hal-03364072

HAL Id: hal-03364072

https://hal.science/hal-03364072

Submitted on 17 Nov 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# The influence of phylogeny and life history on telomere lengths and telomere rate of change among bird species: A meta-analysis

François Criscuolo<fnref n="1" /> ID | F. Stephen Dobson<fnref n="1" />,<fnref n="2" /> ID | Quentin Schull<fnref n="3" /> ID

## Correspondence

François Criscuolo and F. Stephen Dobson, CNRS, Institut Pluridisciplinaire Hubert Curien, UMR 7178, University of Strasbourg, 67000 Strasbourg, France.

Emails: francois.criscuolo@iphc.cnrs.fr (FC); fsdobson@msn.com (FSD)

## Funding information

The present study was supported by a USIAS fellowship and by a Gutenberg Chair awarded to F.S. Dobson.

## Abstract

Longevity is highly variable among animal species and has coevolved with other life-history traits, such as body size and rates of reproduction. Telomeres, through their erosion over time, are one of the cell mechanisms that produce senescence at the cell level and might even have an influence on the rate of aging in whole organisms. However, uneroded telomeres are also risk factors of cell immortalization. The associations of telomere lengths, their rate of change, and life-history traits independent of body size are largely underexplored for birds. To test associations of life-history traits and telomere dynamics, we conducted a phylogenetic meta-analysis using studies of 53 species of birds. We restricted analyses to studies that applied the telomere restriction fragment length (TRF) method, and examined relationships between mean telomere length at the chick (Chick TL) and adult (Adult TL) stages, the mean rate of change in telomere length during life (TROC), and life-history traits. We examined 3 principal components of 12 life-history variables that represented: body size (PC1), the slow-fast continuum of pace of life (PC2), and postfledging parental care (PC3). Phylogeny had at best a small-to-medium influence on Adult and Chick TL (r<fnref n="2" /> = .190 and .138, respectively), but a substantial influence on TROC (r<fnref n="2" /> = .688). Phylogeny strongly influenced life histories: PC1 (r<fnref n="2" /> = .828), PC2 (.838), and PC3 (.613). Adult TL and Chick TL were poorly associated with the life-history variables. TROC, however, was negatively and moderate-to-strongly associated with PC2 (unadjusted r = -.340; with phylogenetic correction, r = -.490). Independent of body size, long-lived species with smaller clutches, and slower embryonic rate of growth may exhibit less change in telomere length over their lifetimes. We suggest that telomere lengths may have diverged, even among closely avian-related species, yet telomere dynamics are strongly linked to the pace of life.

## KEYWORDS

evolution, life history, MCMCglmm, multivariate analysis, phylogenetic correlation, telomere

<figure><img src="image_2.png" /></figure>

<footnote n="1">¹CNRS, Institut Pluridisciplinaire Hubert Curien, UMR 7178, University of Strasbourg, Strasbourg, France</footnote>

<footnote n="2">²Department of Biological Sciences, Auburn University, Auburn, AL, USA</footnote>

<footnote n="3">³MARBEC, IFREMER, IRD, CNRS, University of Montpellier, Sète, France</footnote>
# 1 | INTRODUCTION

Senescence affects most living organisms (Gaillard et al., 2020). With advancing age, individual performance in reproduction, physical activities, or cognition diminish, while concomitantly the rate of actuarial senescence (i.e., mortality rate with age) increases (Nussey et al., 2013). The age at which senescence occurs varies among species and individuals (Gaillard et al., 2020), while energy gathered from the environment may be differently allocated to promote self-maintenance and thus delay aging (Bouwhuis et al., 2010; Jones et al., 2008; Ricklefs, 2008). The view of senescence as an outcome of energy trade-offs that underlie life-history trajectories has shaped our understanding of the evolution of animal longevities. Fast-living species, with rapid growth, early sexual maturity, and high reproduction rates, experience senescence earlier (Jones et al., 2008). However, whether trade-offs are the only determinants of the rate of senescence has been debated, and environmental and biological constraints may also be influential, particularly in comparisons among species (Cohen et al., 2019). In addition, senescence might differently affect phenotypic traits in relation to their impact on fitness (Hayward et al., 2015), with senescence in reproduction beginning at different points in the lifecycles of different species, compared to senescence in other traits that are directly linked to survival (Gaillard & Lemaître, 2017). Studying aging among species will reveal how longevities covary with other organismal traits, and under which trade-offs among patterns of growth, reproduction, and life span (Boonekamp et al., 2020; Cohen et al., 2019). Further, studies among species may also reveal underlying cellular and physiological mechanisms of interest that sustained the diversity of animal longevities.

When comparing species, there is a strongly supported pattern of larger body size associated with longer life expectancy (Lindstedt & Calder, 1976, 1981). This association should reflect physiological traits of large species, like lower mass-specific metabolic rates, which have not been fully characterized (Speakman, 2005a). Metabolism produces damage, as from radical oxygen species (ROS), that are by-products of aerobic respiration. Radical oxygen species were proposed as the main underlying mechanism of aging, thus linking together the coevolution of body size and life span via the higher ROS production of fast-living, and generally small-sized species, with higher rates of oxygen consumption (Beckman & Ames, 1998; Buttemer et al., 2010; Delhaye et al., 2016). However, there are caveats to the so-called free-radical damage theory of aging (e.g., lower ROS production at high metabolism due to uncoupling of mitochondrial respiration; (Speakman et al., 2015). Cohen et al. (2019) called for additional study of age-related mechanistic processes that may be integral to the important trade-offs in play. Thereby, such research might yield a more complete explanation of interspecific variation in life span (Hulbert et al., 2007).

Shortening of telomere length is well known for its role in cell senescence in vitro (Harley et al., 1990). Within species, an increasing number of studies have examined the link between telomere length or its rate of erosion (i.e., telomere rate of change, TROC, expressed

as bp lost per year, a mostly negative variable) with aging and life span (e.g., a seminal paper by Haussmann et al. (2003)). For instance, in recent comparative analyses, short telomeres have been associated with higher risks of mortality in nonmodel and wild animals (Wilbourn et al., 2018) and TROC was found inversely related to maximum life span in birds and mammals (Sudyka et al., 2015; Tricola et al., 2018; Whittemore et al., 2019). However, the sample size in the latter study was small, and a negative relationship between TROC and life span was lacking in intraspecific studies (Beaulieu et al., 2011; Sudyka et al., 2014). In addition, experimental evidence is lacking for a causal link between telomere length or TROC and life span (Simons, 2015; Young, 2018; but see Varela et al., 2016). Telomere length and TROC need not be highly correlated (Tricola et al., 2018), since the telomere length an individual starts with in life with may have a strong heritable basis (Vedder et al., 2021; but see Becker et al., 2015; Viblanc et al., 2020), and both variables may be differently regulated by up-stream cell mechanisms (i.e., oxidative stress or telomerase activity; Grasman et al., 2011). As a consequence, they may not display similar associations with life span when comparing species.

Among species, TROC seems more strongly and negatively associated with mean longevity compared to a lack of significant association of mean telomere length and longevity (Dantzer & Fletcher, 2015; Haussmann et al., 2003; Sudyka et al., 2015; Tricola et al., 2018). Those studies also showed that TROC is phylogenetically conserved within bird families, thereby suggesting that mechanisms modulating telomere loss coevolved with life span similarly among closely related species (Tricola et al., 2018). However, Dantzer and Fletcher (2015) extended analyses to associations of TROC and life-history traits, finding that among 14 bird species, those with a more rapid pace of life also had greater TROC (Dantzer & Fletcher, 2015). The pace of life was calculated using a principal component analysis (PCA) that merged 9 variables on a single PCA axis, including adult body mass. However, body size effects must be removed from pace of life analysis, because the slow-fast continuum in life histories is statistically independent of body size (Gaillard et al., 1989; Read & Harvey, 1989; reviewed by Harvey & Purvis, 1999; Dobson & Oli, 2007). The results of Dantzer and Fletcher (2015), therefore, could not differentiate associations of body size and the slow-fast pace of life with TROC. Body mass was further used as covariate in comparative analysis of interspecific link between TROC and life span among 19 species of birds, and this analysis likely adjusted for variation among species in body size (via body mass versus life span and TROC; (Tricola et al., 2018). Body size may also modulate the interspecific variation of telomere length through its interaction with telomere maintenance enzyme, that is, telomerase, which appears inoperative in larger rodent species (Gorbunova & Seluanov, 2009), perhaps with an associated risk of cell immortalization (Risques & Promislov, 2018; Seluanov et al., 2007; Tian et al., 2018). In addition, body mass may be a key determinant of survival in animals (Miller et al., 2002), including bird species (e.g., Briga, 2016; Stier et al., 2014; but see Briga et al., 2019), and has a significant association with senescence rate among species (Jones et al., 2008). Thus, one needs to
statistically account for differences in body mass among species, to avoid spurious correlative conclusions about relationships between TROC or telomere length and longevity (Speakman, 2005b).

Telomere length at the end of growth seems to predict individual life span (Asghar et al., 2015; Heidinger et al., 2012; Lieshout et al., 2019). Early growth is the period when the greatest amount of telomere loss occurs, due to either rapid cell proliferation or high metabolism during early growth (Monaghan & Ozanne, 2018). Alternatively, longer telomeres early in life may inevitably erode faster (Grasman et al., 2011). Early growth is also the period during which physiological maturation of tissues occurs (Cornell et al., 2017; Ricklefs et al., 1994; Starck & Ricklefs, 1998), a phenomenon influenced by both genetic and environmental factors. From such cell and physiological mechanisms, ontogeny produces the adult phenotype. Rapid growth or stressful early-life conditions may increase telomere loss and shorten individual life span (Boonekamp et al., 2014; Herborn et al., 2014; Metcalfe & Monaghan, 2003; Tarry-Adkins et al., 2009). One hypothesis is that growth and longevity are connected via the processes that shape telomere maintenance during early life. Still, only a tenuous and indirect relationship presently relates short telomeres to reduced long-term survival (Quque et al., 2021) and more rapid growth to increased telomere shortening (Vedder et al., 2018), both in wild birds. However, in a related study, interindividual differences in telomere loss were only moderately sensitive to early-life conditions (i.e., hatching order and resource acquisition ability), termed “canalization of resource allocation” toward telomere maintenance (Vedder et al., 2017). If we hypothesized that this extends to the interspecific level, it suggests selection for reduced erosion of telomere length during growth in long-lived species. Since those species generally exhibit extended developmental periods, slow growth patterns might have been associated with reduced shortening of telomere length. This hypothesis needs to be tested at the phylogenetic level, as well as within long-lived avian groups of species like procellariforms, where body size effects can be statistically controlled due to a broad range in body mass (about 30 g to nearly 10 kg (Dobson & Jouventin, 2007). Further, the relatively small procellariform species *Oceanodroma leucorhroa* exhibits an intriguing apparent telomere lengthening over life and maintains telomerase activity in adult somatic cells (Haussmann et al., 2003, 2004, 2005; Tricola et al., 2018).

In the present study, we tested the association of telomere length and TROC with life span, reproduction, and growth patterns among 53 bird species. Our study described the correlative relationships between life-history traits and telomere variables at the interspecific level. We applied a principal component analysis to typify major elements of life histories, including body size, the pace of life syndrome, and parental care (after, e.g., Bennett & Owens, 2002; Dobson & Jouventin, 2007; Dobson & Oli, 2007). This approach has the advantage of separating body size effects from influences of the pace of life on telomeres and will nicely complement previous study on evolution of telomere length and dynamics in birds (Tricola et al., 2018). We then examined the association of these life-history indices and telomere variables, both without and with statistical

adjustments for phylogenetic influence. Because life span is more closely associated with TROC than with telomere length (i.e., TROC being estimated using between-individual data; Tricola et al., 2018), we predicted that after statistically controlling for body size, species with a fast pace of life should be characterized by higher TROC (i.e., greater telomere loss *per* time unit) and little pattern of change in adult telomere length (*Adult TL*). In addition, we tested the predicted negative association of body mass and telomere length, previously found in rodents (Seluanov et al., 2007). We expected that *Adult TL* would be negatively associated with body size, and TROC would be positively associated with body size (TROC is a negative value, and large-bodied species lose less TL). Predictions for *Chick TL* are more difficult to make. Since early growth is when telomeres are shortened the most and long-lived species have slower developmental rates, other things being equal, fledgling telomere lengths (*Chick TL*), once controlled for body size, might be relatively longer in species with a slow pace of life because telomere shortening is drawn out over a longer period (based on low growth rates and reduced telomere shortening). But *Chick TL* should also be longer in those with a rapid pace of life (based on greater TROC, so long *Chick TL* and short *Adult TL*). Following results of Tricola et al. (2018), we did not expect phylogeny to have a strong influence on *Adult TL* or *Chick TL*, even after controlling for the possible influence of body size (via body mass). On the other hand, phylogeny was found to explain a significant part of interspecific variation of TROC (Tricola et al., 2018). We hypothesized that extracting the influence of body size from the pace of life would lead to a comparable result, if closely related species with identical paces of life have evolved similar TROC independently of their body masses. This expectation seemed reasonable because some orders, like the Procellariiformes, have evolved long life spans along with considerable variation in body masses among species (Dobson & Jouventin, 2007), so they should present different telomere lengths (in relation to body size) but similar TROC (in relation to life span).

## 2 | MATERIAL AND METHODS

Our meta-analysis followed the recommendations of the Preferred Reporting Items for Systematic Reviews and Meta-Analysis (PRISMA) statement (Liberati et al., 2009).

### 2.1 | Data collection

Complete data and references list are presented in the Electronic Supplementary Material (ESM) 1.

We first identified published articles in peer-reviewed journals that reported telomere length or rate of telomere erosion in birds. We (a) screened our own reference lists, (b) backward searched references sections of selected articles from the list, and (c) collected publications using Google Scholar with the following terms: “telomere” OR “telomere shortening” OR “telomere erosion” AND “birds”
OR "avian." For each retrieved paper, methods were individually appraised by one author (FC) to restrict the selection to studies using the telomere restriction fragment (TRF) method to assess telomere length or erosion. Our meta-analysis therefore only relies on published telomere length data expressed in absolute length (kilobases, kb). Using absolute values of telomere length (TL) allowed us to do interspecific comparisons of telomere length relationships with other life-history traits, based on the calculation for each species of a mean adult telomere length (hereafter *Adult TL*), mean chick telomere length (hereafter *Chick TL*), and mean telomere length rate of change over the lifecycle (hereafter *TROC*), allowing interspecific comparisons. We primarily used studies that sampled erythrocytes, with one exception of a study that used tissue attached to feathers. Data extraction is detailed in ESM 2.

Data on life-history traits were compiled for species for which TRF measurements of *Adult TL* were available. ESM 1 references the principal sources used to collect TRF values and 12 life-history variables for each species. We included publications in peer-reviewed journals, online peer-reviewed databases, and books. The protocol for collecting life-history data is described in ESM 2. In particular, we examined (a) species body mass and size, (b) longevity (age at maturity and life span), (c) reproductive patterns (e.g., using egg or clutch sizes), (d) growth patterns (pre- and posthatching growth rates), and (e) postfledging parental care.

We adjusted statistical analyses for sample sizes using the standard procedure for meta-analyses (Liberati et al., 2009). We could not control for the mean age at which telomeres were measured in adult individuals of each species. However, age is likely to be related to life-history variables included in our analysis (e.g., body size, longevity) and likely to scale with body mass. As such, risks of multicollinearity of variables were high. Thus, we used principal component analyses as a variable reduction method (see details below). We conducted a separate analysis of Chick TL to check for specific relationships between life-history traits and telomere length early in the species' lifecycles. In addition, we analyzed associations of TROC and life-history traits. As in previous studies, TROC was estimated as the slope of the linear regression between TL and age over life (i.e., including chick and adults life stages; Haussmann et al., 2003; Tricola et al., 2018).

## 2.2 | Phylogenetic data

The 53 species included in our study comprised 13 orders and 29 families (Figure 1). Maximal life span ranged from 6 to 60 years, female adult body mass ranged from 11 to 8,060 g, and telomere length ranged from 1.7 to 49.3 kb. Using BirdTree, we downloaded 100 trees from http://www.bird.tree.org (De Magalhaes & Costa, 2009; Jetz et al., 2012), focusing on the 53 species of our dataset (producing a .nex tree file). We checked the validity of this first phylogenetic representation by building an additional tree based on DNA sequences of cytochrome b (ESM3). Both alternatives produced congruent trees, but we based our subsequent

analyses on the bird.tree.org based tree, because we lacked a reliable cytochrome b sequence for *Larus andouinii*, *Larus crassirostris*, *Acrocephalus sechellensis*, and *Erythrura gouldiae*. A consensus tree was produced using the Ape package (Paradis et al., 2004) and Phyltools (Revell, 2012). Branch lengths were estimated using Grafen's method (Grafen, 1989). Markov chain Monte Carlo generalized linear mixed model (MCMCglmm in R) (Hadfield, 2010; Hadfield & Nakagawa, 2010) packages were used to estimate the relative amount of variation in adult telomere length and life-history traits that could be explained by shared ancestry (viz., the phylogenetic tree). These last packages were also used to estimate the phylogeny content of associations between telomere and life-history variables.

## 2.3 | Statistical analyses

Full codes for statistical analyses are presented in ESM 4. Life-history variables exhibited a variety of correlations that suggested substantial collinearity, so variable reduction was necessary. We thus conducted a principal component analysis (PCA) using the R package FactoMineR (Lê et al., 2008) to extract independent orthogonal eigenvectors (PC axes) for later phylogenetic and generalized linear mixed modeling analyses (ESM 4). We investigated the interspecific relationships between TL variables and life-history axes as represented by the first three PC axes. We evaluated how phylogeny accounted for overall variance of TL variables and the species scores on PC1, PC2, and PC3, using linear mixed models in a Markov chain Monte Carlo environment (MCMCglmm) (Hadfield, 2010). MCMCglmm allows inclusion of phylogenetic tree files in glmm models, in a Bayesian framework. MCMCglmm also statistically "controls" for phylogenetic dependence among variables that is due to the pattern of shared evolutionary ancestry of the species, so that "phylogeny-adjusted" patterns among variables are revealed (see (Dunn & Møller, 2014). An additional random factor that adjusted for variance in sample sizes of measurements of variables among species was included, by inserting sampling variances into the mev argument of MCMCglmm.

We conducted univariate analyses of associations of telomere variables and each of the principal component axes (body size, pace of life, and parental care) and then multivariate analyses of all 3 life-history axes. The life-history axes were applied in MCMCglmm models that separately examined *Adult TL*, *Chick TL*, and *TROC* (see ESM 4). Phylogeny was entered as a random factor: TL variables = [PC1 or PC2 or PC3] + 1|phylogeny + 1|mESD [adjustment for standard errors of sample sizes, mESD = √(1/(N-3))] + error. These models produced estimates of phylogenetic associations of *Adult TL*, *Chick TL*, *TROC*, and of each of the life-history axes in the univariate analyses. Also, we examined the magnitude and significance of correlations between the 3 telomere variables and all 3 life-history axes, in multivariate analyses. Posterior distributions of covariances were transformed into correlations (using the *posterior.cor()* function in package MCMCglmm). We were then able to evaluate the part of
<figure><img src="image_3.png" /></figure>

FIGURE 1 Phylogenetic tree of the 53 bird species for which adult and chick telomere length (TL) and telomere length rate of change (TROC) were collected from published papers on avian telomeres using the TRF methodology. The consensus phylogenetic tree was obtained from BirdTree.org (100 trees) using ape, apTreeshape, and caper R packages (see ESM 4)

the correlation explained by the phylogenetic tree (i.e., phylogenetic correlation) or not (i.e., residual correlation) for each paired association of life-history PC axes and telomere variables. Finally, using Pearson's correlation test we examined associations of TL and the PC axes in 3 different bird orders (where the sample size >8, i.e., Passeriformes, Procellariiformes, Charadriiformes). These last analyses allowed a preliminary examination of correlative relationships of TL and life-history PC axes within specific categories of birds that shared both history and ecological features. Effect sizes of correlations followed Cohen's (1988) suggested criteria: *r* = .1, small; *r* = .3, medium; *r* = .5, large.

3 | RESULTS

Our literature search produced 1,450 publications. Based on an initial TRF selection (n = 155), we extracted 36 references with TRF values not corrected for an associated variable such as body mass or age. We added two unpublished TRF data measurements obtained from colleagues (also from blood, see ESM 1 and 2). We ended with 53 species identified with mean adult TRF values (among 38 studies, some reported more than one species (e.g., Haussmann et al., 2003; Tricola et al., 2018; Whittemore et al., 2019), 29 chick TRF values, and 35 TROC values (including 19 species from Tricola et al., 2018).
## 3.1 | Life-history traits and Principal Component Analysis

Among the 53 avian species, life-history traits exhibited several medium to strong intercorrelations (ESM 4). In particular, measures of body size and mass were associated with most variables. In order to deal with these strong multiple collinearities, we conducted a principal component analysis. The first principal component (PC1) accounted for nearly 60% of the variance among life-history traits and as expected was strongly associated with adult female mass and other variables that reflect the body size of the species (Table 1). The second principal component (PC2) accounted for about 11% of the variance, and positively reflected reproductive and prehatching development variables, and a negative association with life span. With reproduction and life span negatively associated and statistically independent of body size, PC2 likely reflected the "slow-fast continuum" or pace of life syndrome (reviewed by Dobson & Oli, 2007). Scores of the species on PC2 were ordered from relatively fast life histories (greater reproduction and faster prehatching growth and shorter lives) to relatively slow life histories (longer lives, lower reproduction, and lower prehatching growth). The third principal component (PC3) accounted for about 9% of the variance and primarily reflected the number of days postfledging chicks stay with their parents. We interpreted PC3 as postfledging parental investment in offspring.

## 3.2 | Analyses of TL and life-history traits

TROC was significantly but moderately negatively associated with Adult TL and Chick TL (respectively; *r* = -.349, *p* = .04, *n* = 35; *r* = -.488, *p* = .015, *n* = 24). However, the association of TROC

and Adult TL was strongly influenced by two extreme points, *Larus Andouini* and *Parus major*, which presented long telomere lengths and high rates of annual erosion (ESM 1). If we exclude these two extremes, the association changed in sign and became nonsignificant (*r* = .127, *p* = .48, *n* = 33). Adult TL and Chick TL showed a strong positive association (*r* = .882, *p* < .0001, *n* = 28).

Using the scores of the different avian species on the principal components, we investigated the associations of TL with the three statistically independent life-history axes. The three telomere variables (Adult TL, TROC, and Chick TL) showed small-to-medium and nonsignificant associations with phylogeny (Figure 2). Body size and the slow-fast continuum of pace of life (PC1, PC2) were strongly and significantly associated with the phylogenetic pattern. Parental care (PC3) exhibited a small-to-medium association with phylogeny, but a significant one.

We then investigated the associations of bird TL with the three principal component axes (that displayed phylogenetic patterns) using three different approaches. First, we explored direct Pearson's correlations (without considering any phylogenetic effect) between TL variables and the principal components. This suggested a small-to-medium and nonsignificant correlation of Adult TL with body size (PC1; Table 2a, Figure 3a). TROC exhibited a medium positive association with body size (PC1) and a large negative association with the slow-fast pace of life (PC2; species with greater reproduction and shorter life exhibited more strongly negative TROC; Table 2b, Figure 4a,b). Correlations were at best small and nonsignificant between Chick TL and the life-history axes (Table 2c, Figure 5).

Next, univariate analyses that statistically adjusted both telomere variables and life-history data for the phylogenetic pattern confirmed only some of the ahistorical correlational analyses. First, it revealed that the relatively weak pattern of association of adult TL and body size was likely due to the phylogenetic pattern (the residual

<table><caption>TABLE 1 Principal component (PC) analysis of life-history traits of 53 bird species</caption><thead><tr><th></th><th>PC 1</th><th>PC 2</th><th>PC 3</th></tr></thead><tbody><tr><td>Eigenvalue</td><td>7.070</td><td>1.266</td><td>1.087</td></tr><tr><td>% variance explained</td><td>58.913</td><td>10.552</td><td>9.061</td></tr><tr><td>Incubation (days)</td><td><strong>0.886</strong></td><td>-0.197</td><td>-0.035</td></tr><tr><td>Age at fledging (days)</td><td><strong>0.886</strong></td><td>-0.116</td><td>0.012</td></tr><tr><td>Mean age at sexual maturity (days)</td><td><strong>0.792</strong></td><td>-0.354</td><td>0.001</td></tr><tr><td>Adult structural size</td><td><strong>0.766</strong></td><td>0.287</td><td>0.215</td></tr><tr><td>Female adult body mass (g)</td><td><strong>0.813</strong></td><td>0.326</td><td>-0.010</td></tr><tr><td>Maximal life span (years)</td><td><strong>0.779</strong></td><td><strong>-0.365</strong></td><td>0.159</td></tr><tr><td>Posthatching growth rate (KG)</td><td>-0.688</td><td>-0.152</td><td>-0.091</td></tr><tr><td>Embryonic growth rate (g/day)</td><td>0.723</td><td><strong>0.427</strong></td><td>-0.086</td></tr><tr><td>Clutch size</td><td>-0.532</td><td><strong>0.702</strong></td><td>0.189</td></tr><tr><td>Egg mass (g)</td><td><strong>0.924</strong></td><td>0.119</td><td>-0.150</td></tr><tr><td>Hatching mass (g)</td><td><strong>0.920</strong></td><td>0.245</td><td>-0.141</td></tr><tr><td>Postfledging parental care (days)</td><td>0.114</td><td>-0.067</td><td><strong>0.966</strong></td></tr></tbody></table>

Note: Eigenvalues and the percentage of variance explained by each component are given for the three first PC axes, representing those with an eigenvalue larger than 1 (total variance explained is 79.09%). For each variable, correlation coefficients with the principal component axis are given, with the significant values in bold (contributions > average aleatory contribution, 100/n, *n* = number of variables). *Posthatching growth rate* had a large effect size on PC1, but did not load significantly on any of the three PC axis (indicated in italics).
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** Effect sizes (Rho values) and 95% confidence intervals for phylogenetic effects on bird TL variables (Adult telomere length (Adult TL), telomere length rate of change (TROC), and chick telomere length (Chick TL), and the three principal component axes that reflected bird species' life-history traits (PCA 1—body size, PCA 2—slow-fast continuum of the pace of life, and PCA 3—postfledging parental care). Effects are considered significant when not overlapping zero</figcaption></figure>

“phylogeny-adjusted” correlation was -0.153; Table 2a). The same was evident for the moderate positive association of *TROC* with body size (the residual correlation was -0.108; Table 2b). The negative association of *TROC* and the slow-fast pace of life was confirmed as a strong pattern. Association between *TROC* and life-history axes exhibited especially strong phylogenetic patterns.

We next examined these potential associations while considering the effect of phylogeny on both variance and covariance of *TL* and *PC* variables, using a multivariate phylogenetic framework. The correlation between PCs and *TL* variables due to phylogeny was marginal, since their posterior distributions mostly were broadly centered on or largely overlapped zero (Table 2, Figure 6). Only the negative association of *TROC* and slow-fast pace of life approached significance, and it again exhibited a strong phylogenetic pattern (Table 2b; Figure 6, center panels).

Finally, we investigated whether relationships among PC axis and TL variables might differ among avian orders. We found that the association of *Adult TL* and the slow-fast pace of life (PC2) was significantly positive in the Passeriiformes (*r* = .690, *t* = 3.321, *p* = .006, *n* = 14). However, in the Procellariiformes the association was negative (*r* = -.754, *t* = -2.812, *p* = .031, *n* = 8), that is, opposite to the general nonsignificant interspecific trend (Figure 3b). *TROC* and PC1 and PC2 axes were significantly and negatively correlated in the Procellariiformes (Figure 4; PC1: *r* = -.951, *t* = -6.872, *p* < .001; PC2: *r* = -.890, *t* = -4.368, *p* = .007; both *n* = 7). No significant correlations were found within orders for PC axes and *Chick TL* (Figure 5, all *p* > .25).

## 4 | DISCUSSION

Interest in telomeres has steadily increased in the last decade. Telomere length, as a result of the interaction between telomere loss and DNA maintenance processes like telomerase activity (an enzyme

which actively rebuilds the ever-shortening telomeres; Greider & Blackburn, 1985), has been linked to individual health, cellular aging, and organismal senescence (e.g., Bernardes de Jesus & Blasco, 2012; Verhulst et al., 2016). This has triggered additional interest in understanding how telomere length and its rate of change are related, and the extent to which both may predict longevity or fitness at different periods of life (Bichet et al., 2020; Salomons et al., 2009; Vera et al., 2012). Such links at the species level have been previously studied (Dantzer & Fletcher, 2015; Tricola et al., 2018; Whittemore et al., 2019). Our study aimed to disentangle the phylogenetic and nonphylogenetic correlations among *Adult TL*, *Chick TL*, *TROC*, and three general bird life-history traits.

### 4.1 | Phylogenetic patterns in life-history traits

For our 53 species of birds, life-history traits were largely represented on three orthogonal axes: *body size* (PC1), reproductive and longevity features along a *slow-fast continuum* of the *pace of life* (PC2), and *postfledging parental care* (PC3). Our interpretations of life-history axes followed well-accepted principal based on previous studies of vertebrate species (Gaillard et al., 1989; Harvey & Purvis, 1999; Stearns, 1992). The first axis of the PCA clearly reflected the magnitude of collinear life-history traits, usually interpreted as reflecting variation among species in body size. This conclusion is supported by positive loadings over PC1 for adult structural size and body mass, as well as egg and hatching mass. Some body size-related traits that reflect duration of life also had strong positive loadings on PC1, such as incubation duration, age at fledging, age at sexual maturity, and maximal life span (Speakman, 2005a). The second axis of the PCA appeared to reflect the *slow-fast continuum*. Eigenvectors built within a PC analysis are orthogonal to each other; thus, the species-specific trade-off between reproduction (clutch size, a positive factor loading) and survival (age at sexual maturity and maximal life span, negative factor loadings) on PC2 were statistically independent of the species' adult body sizes (the latter reflected by loadings on PC1). As investment in reproduction decreases, age at maturity increases, and life span increases, an interspecific pattern that is consistent with a fundamental life-history trade-off. Longevity thus has two components, that associated with body size and a residual component associated with the pace of life. Finally, the third axis of the PCA had a strong factor loading only for the duration of post-fledging parental care, that is, the number of days parents care for their offspring once they have left the nest. This axis was statistically independent of both body size and the slow-fast continuum.

Unsurprisingly, variance among species for the three PC axes of life-history traits was strongly associated with the species' phylogenetic histories. Associations of life history with phylogeny are well documented in vertebrate animals (Dobson & Jouventin, 2007; Dobson & Oli, 2007; de Magalhaes et al., 2007). As “body plans” are associated traits that change over evolutionary time, life-history traits are expected to coevolve (Bennett & Owens, 2002; Brown & Sibly, 2006; Roff, 2002; Stearns, 1992). In mammals, a relatively strong trade-off is found between growth and longevity,
<table><caption>**TABLE 2** For 53 bird species, analyses of associations of telomere traits (adult mean telomere length (a, Adult TL), telomere rate of changes (b, TROC), and chick mean telomere length (c, Chick TL)) and life-history traits (body size, PC1; the slow-fast continuum of pace of life, PC2; and parental care, PC3)</caption><thead><tr><th rowspan="2">Unadjusted Pearson's correlation</th><th colspan="2">Univariate MCMC model correcting for phylogeny</th><th colspan="4">Multivariate MCMC model correcting for phylogeny</th></tr><tr><th>Phylogenetic r</th><th>Residual r</th><th>Beta</th><th>p value</th><th>Phylogenetic r</th><th>Residual r</th></tr></thead><tbody><tr><td colspan="7">(a) Adult TL (53 species)</td></tr><tr><td>PC1</td><td>-0.228 (-0.469 to 0.045)</td><td>-0.153</td><td>-0.857</td><td>.14</td><td>-0.797 (-0.948 to 0.603)</td><td>-0.154 (-0.462 to 0.248)</td></tr><tr><td>PC2</td><td>0.147 (-0.128 to 0.401)</td><td>0.384</td><td>1.851</td><td>.19</td><td>0.075 (-0.824 to 0.806)</td><td>0.196 (-0.147 to 0.596)</td></tr><tr><td>PC3</td><td>-0.036 (-0.303 to 0.237)</td><td>-0.082</td><td>-0.148</td><td>.91</td><td>-0.225 (-0.805 to 0.866)</td><td>-0.035 (-0.418 to 0.218)</td></tr><tr><td colspan="7">(b) TROC (35 species)</td></tr><tr><td>PC1</td><td>0.344 (0.012 to 0.608)</td><td>0.905</td><td>-0.108</td><td>23.256</td><td>.067</td><td>0.746 (-0.468 to 0.967)</td><td>-0.186 (-0.530 to 0.336)</td></tr><tr><td>PC2</td><td>-0.541 (-0.740 to -0.253)</td><td>-0.952</td><td>-0.531</td><td>-154.465</td><td>&lt;.001</td><td>-0.886 (-0.979 to 0.358)</td><td>-0.465 (-0.733 to 0.008)</td></tr><tr><td>PC3</td><td>0.117 (-0.225 to 0.434)</td><td>-0.648</td><td>0.199</td><td>154.465</td><td>.351</td><td>-0.650 (-0.869 to 0.740)</td><td>0.250 (-0.203 to 0.538)</td></tr><tr><td colspan="7">(c) Chick TL (29 species)</td></tr><tr><td>PC1</td><td>-0.198 (-0.526 to 0.182)</td><td>-0.846</td><td>-0.057</td><td>-0.362</td><td>.416</td><td>-0.439 (-0.949 to 0.665)</td><td>-0.059 (-0.465 to 0.336)</td></tr><tr><td>PC2</td><td>0.053 (-0.319 to 0.412)</td><td>0.771</td><td>0.178</td><td>0.421</td><td>.732</td><td>-0.451 (-0.862 to 0.846)</td><td>0.179 (-0.365 to 0.438)</td></tr><tr><td>PC3</td><td>-0.004 (-0.370 to 0.363)</td><td>-0.030</td><td>0.028</td><td>0.125</td><td>.924</td><td>0.306 (-0.830 to 0.860)</td><td>0.004 (-0.353 to 0.365)</td></tr></tbody></table>

Note: Pearson's correlation coefficients and both univariate and multivariate MCMC models that measure effects of phylogeny and statistically "adjust" for shared ancestry are shown. Proportion of variance among species attributed to phylogeny (phylogenetic r), and unexplained associations (residual r, associations statistically "controlled" for phylogeny) are presented, with significant results in bold. Beta coefficients show changes in TL measures in kb for each standardized unit change in the principal component axis values.
<figure><img src="image_5.png" /></figure>

FIGURE 3 Adult telomere length (TL, in kb) in relation to life-history axes PC1 (a), PC2 (b), and PC3 (c). Note that none of the relationships turned out to be significant (Table 2). Dashed lines refer to Pearson's linear correlation (PC1, r = -.23; PC2, r = .15; PC3, r = .04, Table 2). Symbols indicate the distribution of the 13 bird orders. Note that two unadjusted correlations reached significance for adult TL/PC2: Procellariiformes, r = -.75, t = -2.18, p = .031, n = 8; Passeriformes, r = .69, t = 3.32, p = .006, n = 14

<figure><img src="image_6.png" /></figure>

FIGURE 4 Telomere length rate of change (TROC, in kb/year) in relation to life-history axes PC1 (a), PC2 (b), and PC3 (c). Only the TROC and PC2 relationship was statistically significant (Table 2). Dashed (nonsignificant) and plain (significant) lines refer to unadjusted for phylogeny Pearson's linear correlation (PC1, r = .34; PC2, r = -.54; PC3, r = .12, Table 2). Symbols indicate the distribution of the 13 bird orders. Note that all unadjusted correlations reached significance for the Procellariiformes, n = 7; PC1: r = -.95, t = -6.97, p < .001; PC2: r = -.92, t = -5.06, p = .004; PC3: r = -.79, t = -2.86, p = .035

long-lived species showing slower postnatal growth rates. In birds, there is strong selection for rapid growth after hatching, in order to quickly reach emancipation (thermal and nutritional) (de Magalhaes et al., 2007). While we found that embryonic growth rate showed the expected negatively related pattern of variation with life span along the slow-fast continuum, posthatching growth rate was

relatively low for the largest species and not associated with the slow-fast continuum. This suggests that embryonic growth rate may have a larger effect on life span evolution than posthatching growth rate. Short-lived birds were characterized by faster embryonic development, a relationship previously characterized for vertebrates in general (Ricklefs, 2006). Such observations need to be remembered
<figure><img src="image_7.png" /><figcaption>**FIGURE 5** Chick telomere length (TL, in kb) in relation to life-history axes PC1 (a), PC2 (b), and PC3 (c). Note that none of the relationships turned out to be significant (Table 2). Dashed lines refer to unadjusted for phylogeny Pearson's linear correlation (n = 29; PC1, r = -.20; PC2, r = .05; PC3, r = -.01, Table 2). Symbols indicates the distribution of the 13 bird orders</figcaption></figure>

when further focusing at how life-history axes relate to *Adult TL* (see below), since the mechanisms of the delayed cost of rapid growth in terms of impaired somatic maintenance and shortened life span are not well understood (Monaghan & Ozanne, 2018).

## 4.2 | Phylogenetic patterns in telomere lengths and TROC

In analyses of associations of TL values and life-history axes, strong phylogenetic values indicate that correlations are not independent over the history of the species. Species that are closely related may show similar patterns. A good example was the correlation between TROC and the body size axis (PC1), for which the “phylogeny-free” correlation was actually different in sign from the unadjusted correlations (note the high phylogenetic correlation; Table 2b). Alternatively, different clades may show different patterns, such as the *Adult TL* and pace of life (PC2) correlation in the Passeriiformes (positive) and in the Procellariiformes (negative) (note the high phylogenetic correlation; Table 2a). After adjustment for phylogeny, residual correlations reveal ahistorical associations, such as the fairly strong negative association of TROC and PC2 (see below and Table 2b).

*Adult* and *Chick TL* among species were poorly explained by phylogeny, while *TROC* showed the strongest phylogenetic signal (Figure 2). Life-history axes, however, were strongly associated with phylogenetic history. Our result from 53 avian species confirms and extends previous work done on 19 bird species (Tricola et al., 2018). These results are important to our understanding of how life histories and telomere length on the one hand, and the proximate factors that modulate telomere erosion on the other

hand, are inter-related. It suggests that similarity in telomere lengths between species is not primarily due to shared ancestry. Moreover, there may be independent evolutionary processes that have operated within different avian orders and families, leading to the present diversity in telomere lengths. If the evolution of life-history traits is strongly constrained by phylogeny, but *Adult* and *Chick TL* are not, it may not be surprising that PC1, 2, and 3 showed only weak associations with telomere lengths measured at both stages of life. However, even after controlling for phylogeny, we found that the remaining variance in PC axes and *Adult* or *Chick TL* was poorly related (Table 2). This supports the hypothesis that telomere length evolved short and long forms at both life stages, even in closely phylogenetically related species that share comparable life-history traits. Our results seem to accord with a low phylogenetic pattern of telomere lengths in mammals (Gomes et al., 2011; Seluanov et al., 2007). For *TROC*, the apparent high variability in telomerase activity in closely related species (Gomes et al., 2011; Tian et al., 2018) suggests that there should be little phylogenetic pattern, as we found.

We found no relationship between PC axes and *Chick TL*. We expected slow-living species to exhibit more preserved and then longer telomere lengths at the end of growth, given that rapid growth rate may enhance telomere shortening (Monaghan & Ozanne, 2018; Vedder et al., 2018). However, it is too early for conclusions about the coevolution of growth and telomere dynamics at the chick stage. Given that telomere lengths in early life are a complex mix of inherited traits and environmental conditions (Dugdale & Richardson, 2018; Eisenberg, 2019; Vedder et al., 2017), more data are needed on embryonic telomere lengths. Further comparative analyses should include ecological variables that may highlight how specific local conditions can influence *Chick TL*.
<figure><img src="image_8.png" /><figcaption>FIGURE 6 Frequency distributions of phylogenetic (a) and residual (b) correlation values produced by the iterations of the multivariate linear models using Monte Carlo methods and Bayesian Markov chain sampling (MCMC) between TROC and life-history axes of bird life-history traits, that is, PC1 (body size, left panels), PC2 (pace of life, middle panels), and PC3 (parental care, right panels). This multivariate model controlled for shared ancestry simultaneously for all variables. Black lines indicate the peak values reported in Table 2 (two last right columns)</figcaption></figure>

TROC (calculated over the chick to adult stages) appeared moderately (our data) or unrelated to mean adult telomere length, such that most species with long telomere lengths after growth were not showing high rates of *bp* loss per year (Tricola et al., 2018). This result did not support the hypothesis that long telomeres erode faster because they are more sensitive to damage, nor that species with short telomeres should exhibit lower TROC over adulthood (as suggested by intraspecific research (Bize et al., 2009; Grasman et al., 2011; Salomons et al., 2009)). Whether this could be the case for TROC

calculated from earlier development remains an open question, and the fact that we found a negative relationship between TROC and *Chick TL* suggests that, during growth, the mechanisms regulating telomere dynamics are different than those in play during adulthood. As for TROC within adults and chicks, these may reflect distinct and nonadditive processes defining species and trade-offs of individual growth, short-term chick survival, and adult life span (Boonekamp et al., 2014), through a tightly regulated balance between energy investment in telomere maintenance and somatic growth. In addition,
our results suggest a coevolution of low TROC with increasing body size along the bird phylogenetic tree. This conflict indicates that more detailed studies are needed on how body size, embryonic development, and reproduction investments may be related to specific telomere dynamics. However, the strong phylogenetic correlation of TROC and PC1 left little evidence of a relationship between these variables under current environmental conditions. In short, the strong phylogenetic pattern rendered the apparent positive phenotypic correlations of TROC and PC1 nonsignificant and slightly negative.

Adult TROC was closer to zero in species that produced smaller clutches and exhibited low embryonic growth rates on the slow-fast continuum (more negative values on the PC2 axis). Of all the potential associations of life histories and telomere dynamics, the association of slow growth and low reproductive rate with less TROC was the only consistent and significant pattern that we found, whether phylogenetically adjusted or not, and univariate or multivariate. Little TROC in slow-growing species with low clutch sizes might be maintained by genetic correlations among these traits. Such negative associations between life-history traits that divert energy from body maintenance and TROC may reflect two mechanisms of selection. First, the selective disappearance of individuals that were submitted to rapid telomere loss early in life, either due to unfavorable growth trajectories or unbalanced reproductive investments. This has been extensively discussed previously (Dantzer & Fletcher, 2015; Tricola et al., 2018): Within individuals, TROC values were significantly higher (i.e., a greater rate of loss) than between-individuals, suggesting the selective disappearance of individuals with short telomeres. However, comparisons of between- and within-individuals TROC relationships with life span were found to be consistent, suggesting that selective disappearance might have a weak effect at the interspecific level. The second possibility is that less TROC actually reflects the coselection of physiological or cellular mechanisms that protect telomere ends from a rapid erosion. In our study, the TROC-PC2 relationship was independent of body size effects, suggesting that it is not specifically in the largest birds that a coevolution of slow growth and low reproductive rates with less TROC has taken place. Because telomerase is expected to be inactive in somatic cells soon after birth in large-bodied species (Seluanov et al., 2007), the underlying mechanisms should be independent of telomerase.

## 4.3 | TROC in the Procellariiformes and the pace of life

In mammals, shorter telomeres prevailed in large-bodied and long-lived species (Gomes et al., 2011; Risques & Promislov, 2018), and short telomeres are associated with repressed telomerase activity in somatic cells at adulthood (Tian et al., 2018). In contrast, we found that bird Adult TL was not related to the body size PC1 axis. The variance in body sizes within birds is less extensive than for mammals, which might partially explain the difference observed between these two major clades of vertebrates; that is, too low a range in

body size of birds to reveal strong negative selection on telomere length. Adjustment for the phylogenetic tree reduced the association to a negligible value (Table 2). Still, our within-orders regression data suggested negative relationships between Adult TL and body size, and TROC and both body size and the slow-fast continuum in Procellariiformes, an order with the extensive range of body sizes of long-lived seabirds. The former relationship suggests that, as previously found in rodents, the small- and long-lived species in the *Hydrobatidae* may have evolved specific mechanisms promoting long life span while preserving telomerase activity and long telomeres over life (Tian et al., 2018). The latter link reveals that the relatively long-living and small-bodied Procellariiformes exhibit long telomeres in adulthood (*Oceanodroma leucorhoa* and *Hydrobates pelagicus* (both *Hydrobatidae*) and to a lesser extent *Fulmarus glacialis*)) along with particularly slow embryonic growth rates (respectively, 0.104, 0.166, and 1.181 vs. a mean procellariform value of 1.938 g/day; note that variation in clutch size along the slow-fast continuum is negligible for procellariiform species because they lay only a single egg each breeding season; ESM1). Thus, long telomeres in this order were not significantly related to body size (PC1), but were related to the slow end of the slow-fast continuum (PC2) and thus to species longevity. This relationship is reversed in the Passeriformes, where long telomeres were related to the fast (short life span) end of the slow-fast continuum of life history.

One of the pathways that could decelerate the rate of cell division, and hence the rate of telomere erosion (Monaghan & Ozanne, 2018), is a prolonged incubation duration (Vedder et al., 2018). Prolonged incubation might also be associated with the maintenance of high telomerase activity throughout life in procellariiform species (Haussmann et al., 2004). An alternative explanation could be that the *Hydrobatidae*, which are cavity nesters where chicks grow in a very contained and stable environment, even after hatching, may exhibit metabolism that favors the maintenance of the lengths of the telomeres (Stier et al., 2019). This pattern supports the antagonistic pleiotropy of telomere length processes, the main hypothesis proposed to explain why species in different phylogenetic clades have evolved large differences in telomere lengths (Hemann & Greider, 2000). This perspective deserves future exploration, notably on telomerase activity in the soma among the species of Procellariformes.

## ACKNOWLEDGMENTS

We are grateful to Dr C. Saraux and Dr B. Warren for helpful comments on the phylogenetic and statistical analyses. The authors wish to thank Dr A. Pauliny, Dr D. Lank, Dr M. Enstipp, Dr P. Bize, and Dr D. Gremillet for providing us with personal information on their biological models. FSD thanks the Institute of Advanced Studies of the University of Strasbourg for their financial support through an USIAS fellowship, and the Région Grand Est and the Eurométropole de Strasbourg for the award of a Gutenberg Excellence Chair.

## CONFLICT OF INTEREST

None declared.
## AUTHOR CONTRIBUTIONS

**François Criscuolo:** Conceptualization (equal); data curation (lead); formal analysis (equal); writing-original draft (equal); writing-review & editing (equal). **F. Stephen Dobson:** Conceptualization (equal); formal analysis (equal); writing-original draft (equal); writing-review & editing (equal). **Quentin Schull:** Conceptualization (equal); formal analysis (equal); writing-original draft (equal); writing-review & editing (equal).

## DATA AVAILABILITY STATEMENT

Data file is provided in the ESM1.

## ORCID

François Criscuolo

<figure><img src="image_9.png" /></figure>

https://orcid.org/0000-0001-8997-8184

F. Stephen Dobson

<figure><img src="image_10.png" /></figure>

https://orcid.org/0000-0001-5562-6316

Quentin Schull

<figure><img src="image_11.png" /></figure>

https://orcid.org/0000-0001-9297-3376

## REFERENCES

Asghar, M., Hasselquist, D., Hansson, B., Zehtindjiev, P., Westerdahl, H., & Bensch, S. (2015). Hidden costs of infection: Chronic malaria accelerates telomere degradation and senescence in wild birds. *Science*, 347(6220), 436–438.

Beaulieu, M., Reichert, S., Le Maho, Y., Ancel, A., & Criscuolo, F. (2011). Oxidative status and telomere length in a long-lived bird facing a costly reproductive event. *Functional Ecology*, 25(3), 577–585. https://doi.org/10.1111/j.1365-2435.2010.01825.x

Becker, P. J. J., Reichert, S., Zahn, S., Hegelbach, J., Massemin, S., Keller, L. F., Postma, E., & Criscuolo, F. (2015). Mother-offspring and nest-mate resemblance but no heritability in early-life telomere length in white-throated dippers. *Proceedings of the Royal Society B: Biological Sciences*, 282(1807), 20142924. https://doi.org/10.1098/rspb.2014.2924

Beckman, K. B., & Ames, B. N. (1998). The free radical theory of aging matures. *Physiological Reviews*, 78(2), 547–581. https://doi.org/10.1152/physrev.1998.78.2.547

Bennett, P. M., & Owens, I. P. F. (2002). *Evolutionary ecology of birds: Life histories, mating systems and extinction*. Oxford University Press.

Bernardes de Jesus, B., & Blasco, M. A. (2012). Potential of telomerase activation in extending health span and longevity. *Current Opinion in Cell Biology*, 24(6), 739–743. https://doi.org/10.1016/j.ceb.2012.09.004

Bichet, C., Bouwhuis, S., Bauch, C., Verhulst, S., Becker, P. H., & Vedder, O. (2020). Telomere length is repeatable, shortens with age and reproductive success, and predicts remaining lifespan in a long-lived seabird. *Molecular Ecology*, 29(2), 429–441. https://doi.org/10.1111/mec.15331

Bize, P., Criscuolo, F., Metcalfe, N. B., Nasir, L., & Monaghan, P. (2009). Telomere dynamics rather than age predict life expectancy in the wild. *Proceedings of the Royal Society B: Biological Sciences*, 276(1662), 1679–1683. https://doi.org/10.1098/rspb.2008.1817

Boonekamp, J. J., Bauch, C., & Verhulst, S. (2020). Experimentally increased brood size accelerates actuarial senescence and increases subsequent reproductive effort in a wild bird population. *Journal of Animal Ecology*, 89(6), 1395–1407. https://doi.org/10.1111/1365-2656.13186

Boonekamp, J. J., Mulder, G. A., Salomons, H. M., Dijkstra, C., & Verhulst, S. (2014). Nestling telomere shortening, but not telomere length, reflects developmental stress and predicts survival in wild birds. *Proceedings of the Royal Society B: Biological Sciences*, 281(1785), 20133287. https://doi.org/10.1098/rspb.2013.3287

Bouwhuis, S., Charmantier, A., Verhulst, S., & Sheldon, B. C. (2010). Individual variation in rates of senescence: Natal origin effects and disposable soma in a wild bird population. *Journal of Animal Ecology*, 79(6), 1251–1261. https://doi.org/10.1111/j.1365-2656.2010.01730.x

Briga, M. (2016). *Growing up and growing old: A longitudinal study of aging in zebra finches*. University of Groningen. https://doi.org/10.13140/rg.2.2.12893.00481

Briga, M., Jimeno, B., & Verhulst, S. (2019). Coupling lifespan and aging? The age at onset of body mass decline associates positively with sex-specific lifespan but negatively with environment-specific lifespan. *Experimental Gerontology*, 119, 111–119. https://doi.org/10.1016/j.exger.2019.01.030

Brown, J. H., & Sibly, R. M. (2006). Inaugural Article: Life-history evolution under a production constraint. *Proceedings of the National Academy of Sciences*, 103(47), 17595–17599. https://doi.org/10.1073/pnas.0608522103

Buttemer, W. A., Abele, D., & Costantini, D. (2010). From bivalves to birds: Oxidative stress and longevity. *Functional Ecology*, 24(5), 971–983. https://doi.org/10.1111/j.1365-2435.2010.01740.x

Cohen, A. A., Coste, C. F. D., Li, X.-Y., Bourg, S., Pavard, S., & Gaillard, J.-M. (2019). Are trade-offs really the key drivers of ageing and life span? *Functional Ecology*, 34(1), 153–166. https://doi.org/10.1111/1365-2435.13444

Cohen, J. (1988). *Statistical power analysis for the behavioral sciences* (2nd ed.). Erlbaum.

Cornell, A., Gibson, K. F., Williams, T. D., & Portugal, S. (2017). Physiological maturity at a critical life-history transition and flight ability at fledging. *Functional Ecology*, 31(3), 662–670. https://doi.org/10.1111/1365-2435.12777

Dantzer, B., & Fletcher, Q. E. (2015). Telomeres shorten more slowly in slow-aging wild animals than in fast-aging ones. *Experimental Gerontology*, 71, 38–47. https://doi.org/10.1016/j.exger.2015.08.012

De Magalhaes, J. P., & Costa, J. (2009). A database of vertebrate longevity records and their relation to other life-history traits. *Journal of Evolutionary Biology*, 22(8), 1770–1774. https://doi.org/10.1007/s11357-016-9940-z

de Magalhaes, J. P., Costa, J., & Church, G. M. (2007). An analysis of the relationship between metabolism, developmental schedules, and longevity using phylogenetic independent contrasts. *Journal of Gerontology*, 62A(2), 149–160. https://doi.org/10.1093/gerona/62.2.149

Delhaye, J., Salamin, N., Roulin, A., Criscuolo, F., Bize, P., & Christe, P. (2016). Interspecific correlation between red blood cell mitochondrial ROS production, cardiolipin content and longevity in birds. *Age*, 38(5-6), 433–443. https://doi.org/10.1007/s11357-016-9940-z

Dobson, F. S., & Jouventin, P. (2007). How slow breeding can be selected in seabirds: Testing Lack's hypothesis. *Proceedings of the Royal Society B: Biological Sciences*, 274(1607), 275–279. https://doi.org/10.1098/rspb.2006.3724

Dobson, F. S., & Oli, M. K. (2007). Fast and slow life histories of mammals. *Ecoscience*, 14(3), 292–299.

Dugdale, H. L., & Richardson, D. S. (2018). Heritability of telomere variation: It is all about the environment! *Philosophical Transactions of the Royal Society B: Biological Sciences*, 373(1741), 20160450. https://doi.org/10.1098/rstb.2016.0450

Dunn, P. O., & Møller, A. P. (2014). Changes in breeding phenology and population size of birds. *Journal of Animal Ecology*, 83(3), 729–739. https://doi.org/10.1111/1365-2656.12162

Eisenberg, D. T. A. (2019). Paternal age at conception effects on offspring telomere length across species—What explains the variability? *PLoS Genetics*, 15(2), e1007946. https://doi.org/10.1371/journal.pgen.1007946.t001

Gaillard, J.-M., & Lemaître, J.-F. (2017). The Williams' legacy: A critical reappraisal of his nine predictions about the evolution of senescence. *Evolution*, 71(12), 2768–2785. https://doi.org/10.1111/evo.13379

Gaillard, J.-M., Lemaître, J.-F., & Fox, C. (2020). An integrative view of senescence in nature. *Functional Ecology*, 34(1), 4–16. https://doi.org/10.1111/1365-2435.13506

Gaillard, J.-M., Pontier, D., Allainé, D., Lebreton, J. D., Trouvilliez, J., Clobert, J., & Allaine, D. (1989). An analysis of demographic tactics in birds and mammals. *Oikos*, 56, 59–76. https://doi.org/10.2307/3566088

Gomes, N. M. V., Ryder, O. A., Houck, M. L., Charter, S. J., Walker, W., Forsyth, N. R., Austad, S. N., Venditti, C., Pagel, M., Shay, J. W., &
Wright, W. E. (2011). Comparative biology of mammalian telomeres: Hypotheses on ancestral states and the roles of telomeres in longevity determination. *Aging Cell*, **10**, 761-768. https://doi.org/10.1111/j.1474-9726.2011.00718.x

Gorbunova, V., & Seluanov, A. (2009). Coevolution of telomerase activity and body mass in mammals: From mice to beavers. *Mechanisms of Ageing and Development*, **130**(1-2), 3-9. https://doi.org/10.1016/j.mad.2008.02.008

Grafen, A. (1989). The phylogenetic regression. *Philosophical Transactions of the Royal Society of London Series B: Biological Sciences*, **326**, 119-157.

Grasman, J., Salomons, H. M., & Verhulst, S. (2011). Stochastic modeling of length-dependent telomere shortening in *Corvus monedula*. *Journal of Theoretical Biology*, **282**(1), 1-6. https://doi.org/10.1016/j.jtbi.2011.04.026

Greider, C. W., & Blackburn, E. H. (1985). Identification of a specific telomere terminal transferase-activity in tetrahymena extracts. *Cell*, **43**(2), 405-413. https://doi.org/10.1016/0092-8674(85)90170-9

Hadfield, J. D. (2010). MCMC methods for multi-response generalized linear mixed models: The MCMCglmm R package. *Journal of Statistical Software*, **33**(2), 1-22.

Hadfield, J. D., & Nakagawa, S. (2010). General quantitative genetic methods for comparative biology: Phylogenies, taxonomies and multi-trait models for continuous and categorical characters. *Journal of Evolutionary Biology*, **23**(3), 494-508. https://doi.org/10.1111/j.1420-9101.2009.01915.x

Harley, C. B., Futcher, A. B., & Greider, C. W. (1990). Telomeres shorten during ageing of human fibroblasts. *Nature*, **345**, 458-460. https://doi.org/10.1038/345458a0

Harvey, P. H., & Purvis, A. (1999). Understanding the ecological and evolutionary reasons for life history variation: Mammals as a case study. In J. McGlade (Ed.), *Advanced ecological theory: Principles and applications* (pp. 232-348). Blackwell Science.

Haussmann, M. F., Winkler, D. W., Huntington, C. E., Nisbet, I. C. T., & Vleck, C. M. (2004). Telomerase expression is differentially regulated in birds of differing life span. *Annals of the New York Academy of Sciences*, **1019**(1), 186-190. https://doi.org/10.1196/annals.1297.029

Haussmann, M. F., Winkler, D. W., O'Reilly, K. M., Huntington, C. E., Nisbet, I. C. T., & Vleck, C. M. (2003). Telomeres shorten more slowly in long-lived birds and mammals than in short-lived ones. *Proceedings of the Royal Society B: Biological Sciences*, **270**(1522), 1387-1392. https://doi.org/10.1098/rspb.2003.2385

Haussmann, M. F., Winkler, D. W., & Vleck, C. M. (2005). Longer telomeres associated with higher survival in birds. *Biology Letters*, **1**(2), 212-214. https://doi.org/10.1098/rbl.2005.0301

Hayward, A. D., Moorad, J., Regan, C. E., Berenos, C., Pilkington, J. G., Pemberton, J. M., & Nussey, D. H. (2015). Asynchrony of senescence among phenotypic traits in a wild mammal population. *Experimental Gerontology*, **71**, 56-68. https://doi.org/10.1016/j.exger.2015.08.003

Heidinger, B. J., Blount, J. D., Boner, W., Griffiths, K., Metcalfe, N. B., & Monaghan, P. (2012). Telomere length in early life predicts lifespan. *Proceedings of the National Academy of Sciences*, **109**, 1742-1748. https://doi.org/10.1073/pnas.1113306109

Hemann, M. T., & Greider, C. W. (2000). Wild-derived inbred mouse strains have short telomeres. *Nucleic Acids Research*, **28**(22), 4474-4478. https://doi.org/10.1093/nar/28.22.4474

Herborn, K. A., Heidinger, B. J., Boner, W., Noguera, J. C., Adam, A., Daunt, F., & Monaghan, P. (2014). Stress exposure in early post-natal life reduces telomere length: An experimental demonstration in a long-lived seabird. *Proceedings of the Royal Society B: Biological Sciences*, **281**(1782), 2013-3151. https://doi.org/10.1098/rspb.2013.3151

Hulbert, A., Pamplona, R., Buffenstein, R., & Buttemer, W. (2007). Life and death: Metabolic rate, membrane composition, and life span of animals. *Physiological Reviews*, **87**, 1176-1213. https://doi.org/10.1152/physrev.00047.2006

Jetz, W., Thomas, G. H., Joy, J. B., Hartmann, K., & Mooers, A. O. (2012). The global diversity of birds in space and time. *Nature*, **491**(7424), 444-448. https://doi.org/10.1038/nature11631

Jones, O. R., Gaillard, J.-M., Tuljapurkar, S., Alho, J. S., Armitage, K. B., Becker, P. H., Bize, P., Brommer, J., Charmantier, A., Charpentier, M., Clutton-Brock, T., Dobson, F. S., Festa-Bianchet, M., Gustafsson, L., Jensen, H., Jones, C. G., Lillandt, B.-G., McCleery, R., Merilä, J., ... Coulson, T. (2008). Senescence rates are determined by ranking on the fast-slow life-history continuum. *Ecology Letters*, **11**(7), 664-673. https://doi.org/10.1111/j.1461-0248.2008.01187.x

Lë, S., Josse, J., & Husson, F. (2008). FactoMineR: An R package for multivariate analysis. *Journal of Statistical Software*, **25**(1), 2-18.

Liberati, A., Altman, D. G., Tetzlaff, J., Mulrow, C., Götzsche, P. C., Ioannidis, J. P. A., Clarke, M., Devereaux, P. J., Kleijnen, J., & Moher, D. (2009). The PRISMA statement for reporting systematic reviews and meta-analyses of studies that evaluate health care interventions: Explanation and elaboration. *PLoS Med*, **6**(7), e1000100. https://doi.org/10.1371/journal.pmed.1000100

Lieshout, S. H. J., Bretman, A., Newman, C., Buesching, C. D., Macdonald, D. W., & Dugdale, H. L. (2019). Individual variation in early-life telomere length and survival in a wild mammal. *Molecular Ecology*, **28**, 4152-4165. https://doi.org/10.1111/mec.15212

Lindstedt, S. L., & Calder, W. A. (1976). Body size and longevity in birds. *The Condor*, **78**, 91-145. https://doi.org/10.2307/1366920

Lindstedt, S. L., & Calder, W. A. (1981). Body size, physiological time, and longevity of homeothermic animals. *The Quarterly Review of Biology*, **56**(1), 1-16. https://doi.org/10.1086/412080

Metcalfe, N., & Monaghan, P. (2003). Growth versus lifespan: Perspectives from evolutionary ecology. *Experimental Gerontology*, **38**(9), 935-940. https://doi.org/10.1016/s0531-5565(03)00159-1

Miller, R. A., Harper, J. M., Galecki, A., & Burke, D. T. (2002). Big mice die young: Early life body weight predicts longevity in genetically heterogeneous mice. *Aging Cell*, **1**, 22-29. https://doi.org/10.1046/j.1474-9728.2002.00006.x

Monaghan, P., & Ozanne, S. E. (2018). Somatic growth and telomere dynamics in vertebrates: Relationships, mechanisms and consequences. *Philosophical Transactions of the Royal Society B: Biological Sciences*, **373**(1741), 2016-0446. https://doi.org/10.1098/rstb.2016.0446

Nussey, D. H., Froy, H., Lemaitre, J.-F., Gaillard, J.-M., & Austad, S. N. (2013). Senescence in natural populations of animals: Widespread evidence and its implications for bio-gerontology. *Ageing Research Reviews*, **12**(1), 214-225. https://doi.org/10.1016/j.arr.2012.07.004

Paradis, E., Claude, J., & Strimmer, K. (2004). APE: Analyses of phylogenetics and evolution in R language. *Bioinformatics*, **20**, 289-290. https://doi.org/10.1093/bioinformatics/btg412

Quque, M., Paquet, M., Zahn, S., Théron, F., Faivre, B., Sueur, C., & Covas, R. (2021). Contrasting associations between nestling telomere length and pre and postnatal helpers' presence in a cooperatively breeding bird. *Oecologia*, **196**(1), 37-51.

Read, A. F., & Harvey, P. H. (1989). Life history differences among the eutherian radiations. *Journal of Zoology*, **219**, 329-353. https://doi.org/10.1111/j.1469-7998.1989.tb02584.x

Revell, L. J. (2012). phytools: An R package for phylogenetic comparative biology (and other things). *Methods in Ecology and Evolution*, **3**(2), 217-223. https://doi.org/10.1111/j.2041-210X.2011.00169.x

Ricklefs, R. E. (2006). Embryo development and ageing in birds and mammals. *Proceedings of the Royal Society B: Biological Sciences*, **273**(1597), 2077-2082. https://doi.org/10.1098/rspb.2006.3544

Ricklefs, R. E. (2008). The evolution of senescence from a comparative perspective. *Functional Ecology*, **22**(3), 379-392. https://doi.org/10.1111/j.1365-2435.2008.01420.x

Ricklefs, R. E., Shea, R. E., & Choi, I.-H. (1994). Inverse relationship between functional maturity and exponential growth rate of avian skeletal muscle: A constraint on evolutionary response. *Evolution*, **48**(4), 1080-1088. https://doi.org/10.1111/j.1558-5646.1994.tb05295.x
Risques, R. A., & Promislov, D. E. L. (2018). All's well that ends well: Why large species have short telomeres. *Philosophical Transactions of the Royal Society B: Biological Sciences*, 373, 20160448. https://doi.org/10.1098/rstb.2016.0448

Roff, D. A. (2002). *Life history evolution*. Sinauer Associates.

Salomons, H. M., Mulder, G. A., van de Zande, L., Haussmann, M. F., Linskens, M. H. K., & Verhulst, S. (2009). Telomere shortening and survival in free-living corvids. *Proceedings of the Royal Society B: Biological Sciences*, 276(1670), 3157-3165. https://doi.org/10.1098/rspb.2009.0517

Seluanov, A., Chen, Z., Hine, C., Sasahara, T. H. C., Ribeiro, A. A. C. M., Catania, K. C., Presgraves, D. C., & Gorbunova, V. (2007). Telomerase activity coevolves with body mass, not lifespan. *Aging Cell*, 6, 45-52. https://doi.org/10.1111/j.1474-9726.2006.00262.x

Simons, M. J. P. (2015). Questioning causal involvement of telomeres in aging. *Ageing Research Reviews*, 24, 191-196. https://doi.org/10.1016/j.arr.2015.08.002

Speakman, J. R. (2005a). Body size, energy metabolism and lifespan. *Journal of Experimental Biology*, 208(9), 1717-1730. https://doi.org/10.1242/jeb.01556

Speakman, J. R. (2005b). Correlations between physiology and lifespan – Two widely ignored problems with comparative studies. *Aging Cell*, 4, 167-175. https://doi.org/10.1111/j.1474-9726.2005.00162.x

Speakman, J. R., Blount, J. D., Bronikowski, A. M., Buffenstein, R., Isaksson, C., Kirkwood, T. B. L., Monaghan, P., Ozanne, S. E., Beaulieu, M., Briga, M., Carr, S. K., Christensen, L. L., Cochemé, H. M., Cram, D. L., Dantzer, B., Harper, J. M., Jurk, D., King, A., Noguera, J. C., ... Selman, C. (2015). Oxidative stress and life histories: Unresolved issues and current needs. *Ecology and Evolution*, 5(24), 5745-5757. https://doi.org/10.1002/ece3.1790

Starck, J. M., & Ricklefs, R. E. (1998). Patterns of development: The altricial-precocial spectrum. *Oxford Ornithology Series*, 8, 3-30.

Stearns, S. C. (1992). *The evolution of life histories*. Oxford University Press.

Stier, A., Metcalfe, N. B., & Monaghan, P. (2019). Ageing before birth: Pace and stability of prenatal growth affect telomere dynamics. *bioRxiv* 809087, https://doi.org/10.1101/809087

Stier, A., Viblanc, V. A., Massemin-Challet, S., Handrich, Y., Zahn, S., Rojas, E. R., Saraux, C., Le Vaillant, M., Prud'homme, O., Grosbellet, E., Robin, J.-P., Bize, P., & Criscuolo, F. (2014). Starting with a handicap: Phenotypic differences between early- and late-born king penguin chicks and their survival correlates. *Functional Ecology*, 28(3), 601-611. https://doi.org/10.1111/1365-2435.12204

Sudyka, J., Arct, A., Drobniak, S., Dubiec, A., Gustafsson, L., & Cichon, M. (2014). Experimentally increased reproductive effort alters telomere length in the blue tit (*Cyanistes caeruleus*). *Journal of Evolutionary Biology*, 27(10), 2258-2264. https://doi.org/10.1111/jeb.12479

Sudyka, J., Arct, A., Drobniak, S., Gustafsson, L., & Cichoń, M. (2015). Longitudinal studies confirm faster telomere erosion in short-lived bird species. *Journal of Ornithology*, 157(1), 373-375. https://doi.org/10.1007/s10336-015-1304-4

Tarry-Adkins, J. L., Chen, J. H., Smith, N. S., Jones, R. H., Cherif, H., & Ozanne, S. E. (2009). Poor maternal nutrition followed by accelerated postnatal growth leads to telomere shortening and increased markers of cell senescence in rat islets. *The FASEB Journal*, 23(5), 1521-1528. https://doi.org/10.1096/fj.08-122796

Tian, X., Doerig, K., Park, R., Can Ran Qin, A., Hwang, C., Neary, A., Gilbert, M., Seluanov, A., & Gorbunova, V. (2018). Evolution of telomere maintenance and tumour suppressor mechanisms across mammals. *Philosophical Transactions of the Royal Society B: Biological Sciences*, 373, 20160443. https://doi.org/10.1098/rstb.2016.0443

Tricola, G. M., Simons, M. J. P., Atema, E., Boughton, R. K., Brown, J. L., Dearborn, D. C., Divoky, G., Eimes, J. A., Huntington, C. E., Kitaysky, A. S., Juola, F. A., Lank, D. B., Litwa, H. P., Mulder, E. G. A., Nisbet,

I. C. T., Okanoya, K., Safran, R. J., Schoech, S. J., Schreiber, E. A., ... Haussmann, M. F. (2018). The rate of telomere loss is related to maximum lifespan in birds. *Philosophical Transactions of the Royal Society B: Biological Sciences*, 373(1741), 20160445. https://doi.org/10.1098/rstb.2016.0445

Varela, E., Muñoz-Lorente, M. A., Tejera, A. M., Ortega, S., & Blasco, M. A. (2016). Generation of mice with longer and better preserved telomeres in the absence of genetic manipulations. *Nature Communications*, 7, 11739. https://doi.org/10.1038/ncomms11739

Vedder, O., Moiron, M., Bichet, C., Bauch, C., Verhulst, S., Becker, P. H., & Souwhuis, S. (2021). Telomere length is heritable and genetically correlated with lifespan in a wild bird. *Molecular Ecology*, 1-11. https://doi.org/10.1111/mec.15807

Vedder, O., Verhulst, S., Bauch, C., & Bouwhuis, S. (2017). Telomere attrition and growth: A life-history framework and case study in common terns. *Journal of Evolutionary Biology*, 30(7), 1409-1419. https://doi.org/10.1111/jeb.13119

Vedder, O., Verhulst, S., Zuidersma, E., & Bouwhuis, S. (2018). Embryonic growth rate affects telomere attrition: An experiment in a wild bird. *The Journal of Experimental Biology*, 221(15), jeb181586. https://doi.org/10.1242/jeb.181586

Vera, E., Bernardes de Jesus, B., Foronda, M., Flores, J. M., & Blasco, M. A. (2012). The rate of increase of short telomeres predicts longevity in mammals. *Cell Reports*, 2(4), 732-737. https://doi.org/10.1016/j.celrep.2012.08.023

Verhulst, S., Dalgård, C., Labat, C., Kark, J. D., Kimura, M., Christensen, K., Toupance, S., Aviv, A., Kyvik, K. O., & Benetos, A. (2016). A short leucocyte telomere length is associated with development of insulin resistance. *Diabetologia*, 59(6), 1258-1265. https://doi.org/10.1007/s00125-016-3915-6

Viblanc, V. A., Schull, Q., Stier, A., Durand, L., Lefol, E., Robin, J.-P., Zahn, S., Bize, P., & Criscuolo, F. (2020). Foster rather than biological parental telomere length predicts offspring survival and telomere length in king penguins. *Molecular Ecology*, 29(16), 3154-3166. https://doi.org/10.1111/mec.15485

Whittemore, K., Vera, E., Martínez-Nevado, E., Sanpera, C., & Blasco, M. A. (2019). Telomere shortening rate predicts species life span. *Proceedings of the National Academy of Sciences*, 116(30), 15122-15127. https://doi.org/10.1073/pnas.1902452116

Wilbourn, R. V., Moatt, J. P., Froy, H., Walling, C. A., Nussey, D. H., & Boonekamp, J. J. (2018). The relationship between telomere length and mortality risk in non-model vertebrate systems: A meta-analysis. *Philosophical Transactions of the Royal Society B: Biological Sciences*, 373, 20160447. https://doi.org/10.1098/rstb.2016.0447

Young, A. J. (2018). The role of telomeres in the mechanisms and evolution of life-history trade-offs and ageing. *Philosophical Transactions of the Royal Society B: Biological Sciences*, 373(1741), 20160452. https://doi.org/10.1098/rstb.2016.0452

## SUPPORTING INFORMATION

Additional supporting information may be found online in the Supporting Information section.

**How to cite this article:** Criscuolo, F., Dobson, F. S., & Schull, Q. (2021). The influence of phylogeny and life history on telomere lengths and telomere rate of change among bird species: A meta-analysis. *Ecology and Evolution*, 11, 12908-12922. https://doi.org/10.1002/ece3.7931