Skip to main content

Large mammal telomere length variation across ecoregions

Abstract

Background

Telomere length provides a physiological proxy for accumulated stress in animals. While there is a growing consensus over how telomere dynamics and their patterns are linked to life history variation and individual experience, knowledge on the impact of exposure to different stressors at a large spatial scale on telomere length is still lacking. How exposure to different stressors at a regional scale interacts with individual differences in life history is also poorly understood. To better understand large-scale regional influences, we investigated telomere length variation in moose (Alces alces) distributed across three ecoregions. We analyzed 153 samples of 106 moose representing moose of both sexes and range of ages to measure relative telomere lengths (RTL) in white blood cells.

Results

We found that average RTL was significantly shorter in a northern (montane) and southern (sarmatic) ecoregion where moose experience chronic stress related to severe summer and winter temperatures as well as high anthropogenic land-use compared to the boreal region. Our study suggests that animals in the northern boreal forests, with relatively homogenous land use, are less disturbed by environmental and anthropogenic stressors. In contrast, animals in areas experiencing a higher rate of anthropogenic and environmental change experience increased stress.

Conclusion

Although animals can often adapt to predictable stressors, our data suggest that some environmental conditions, even though predictable and ubiquitous, can generate population level differences of long-term stress. By measuring RTL in moose for the first time, we provide valuable insights towards our current understanding of telomere biology in free-ranging wildlife in human-modified ecosystems.

Background

Human-induced rapid environmental change is creating novel stressors for animals and their populations [1]. These external changes cascade via physiological mechanisms affecting long-term survival and fitness in wild animals. In particular, exposure to anthropogenic perturbations (resource extraction, infrastructural developments, hunting, and pollution) combined with environmental stressors (competition over resources, disease, or thermal stress) may activate the hypothalamic–pituitary–adrenal (HPA) axis of animals resulting in increased stress hormone levels [2,3,4]. Continued activation of the HPA axis beyond baseline levels can affect the metabolic system of the organism via increased oxidative damage from reactive oxygen species (ROS), and induce a state of chronic stress [5]. The (TTAGGG)n repeats that constitute vertebrate telomeres are particularly vulnerable to oxidative attack [6]. Telomeres, i.e., the non-coding ends of linear chromosomes, are considered to play a fundamental role in the protection of the structural integrity of chromosomal DNA and in the regulation of cellular senescence [7, 8]. Thus, they have the potential to serve as a molecular biomarker to determine individual physiological state and past environmental experiences [9, 10]. Shorter telomeres and elevated shortening rates are typically associated with stress and senescence [7, 9, 11, 12]. Angelier et al. [13] reviewed studies determining the relationships between how different stressors can influence telomere associations in wild vertebrates. Specifically, environmental factors such as water temperature [14], weather [15], habitat quality [16,17,18] as well as infectious diseases [19] were linked to altered telomere length in wildlife.

Ecoregions provide an ideal spatial scale to examine differences in metabolic expenditure and chronic stress expression as they offer a global categorization representing distinct units of biological diversity and its association with climatic conditions [20]. Distinct ecoregions also encompass differences in anthropogenic pressures, food availability and weather. Differences in the degree of exposure to different environmental conditions (including an array of stressors) can potentially cause chronic stress in organisms occupying ecoregions where they experience repeated triggering of the HPA axis beyond full recovery during the annual and seasonal cycles. Yet, comparative studies of chronic stress responses, or its indicators, across biogeographic regions are largely absent. This is especially true as data on multiple individuals and populations distributed across large spatial scales are not often compared.

The main objective of this study is to compare relative telomere length (RTL) across ecoregions and therefore identify how levels of anthropogenic and environmental stress may correlate with RTL of individuals across multiple populations. Our focal study species is the moose (Alces alces) across the main three ecoregions in Sweden. Shorter term stress in response to anthropogenic and environmental stressors have been demonstrated in moose previously [21] and their longevity (up to ~ 20 years) makes them an ideal model species to also evaluate accumulated stress across an individual’s life span and to compare these across ecoregions. In addition, moose are a cold-adapted species and are susceptible to heat-stress at ambient temperatures above 14–17 °C [22, 23] during summer and above − 5–0 °C during winter [22], meaning they may be particularly susceptible to temperature changes brought about by climate change. In combination with a known higher parasite burden [24, 25], higher hunting pressure and higher inter-species competition [26, 27] that moose are exposed to in their southern range, we expect the chronic stress burden of moose to decrease with increasing latitude. However, moose at high altitudes, i.e., montane tundra habitat, may experience stress due to other factors, such as high snow depth [28]. With this study we analyze RTL in moose for the first time and examine how it reflects chronic stress of individuals experiencing varying levels of environmental factors and anthropogenic impacts across large spatial scales.

Results

Geographic variation in climate and land use

In line with our general hypothesis, we observed a marked difference in GPS-collar recorded temperature [F(2,149) = 111.4, P-value < 0.001] and land-use intensity [F(2,150) = 404.4, P-value < 0.001] that moose experienced in each ecoregion based on their annual movements (Fig. 1). Mean annual temperature (based on GPS-collar temperature: Tc) was 5.44 ± 4.08 °C in the montane, 9.80 ± 2.71 °C in the boreal and 14.16 ± 2.10 °C in the sarmatic ecoregion (Fig. 1). Land use intensity follows a similar trend with low mean global Human Modification (gHM, [29]) values encountered by moose in the montane (0.06 ± 0.05) and boreal (0.05 ± 0.02) but high mean values (0.34 ± 0.08) encountered in the sarmatic ecoregion (Fig. 1).

Fig. 1
figure 1

A Capture locations of 106 moose in three ecoregions (dark grey = montane birch forest and grasslands, grey = boreal forest, light grey = sarmatic mixed forest). The map was created in Quantum GIS, version 6.10.6 (QGIS.org, 2020). Right: Mean annual GPS-collar temperature B and corresponding mean global human modification (gHM) values extracted based on the annual GPS track C of all individuals distinguished by ecoregion

Relative telomere length

The ecoregions variable explained significant differences in RTL as per the final model (Fig. 2; Table 1). Compared to the boreal region, RTLs were significantly shorter in the sarmatic study areas (1.42 [1.31, 1.53] 95% CI) in southern Sweden. Additionally, shorter RTLs were also observed in the northern montane area (1.35 [1.20, 1.44] 95% CI) compared to the boreal region (1.63 [1.49, 1.76] 95% CI). Sample storage time was negatively correlated with RTL. Based on linear mixed effect model selection, sex and age of animals did not influence RTL significantly and were subsequently removed as explanatory variables (Additional file 1: Table S1, Fig. S1–2). Pregnancy and the number of calves at heel did not affect RTL (Additional file 1: Fig. S3).

Fig. 2
figure 2

Average observed relative telomere length (RTL) of 153 samples in three ecoregions (nmontane = 58, nboreal = 29, nsarmatic = 66) from 106 animals. Animal age was included for each sample as black dots, despite having been excluded from the final model

Table 1 The best linear mixed effect model showing the relationship between relative telomere length of moose individuals (N = 106), the three considered ecoregions in Sweden, and storage time

Discussion

Our results show how RTL, an indicator of chronic stress, can vary across different ecoregions. Our study provides the first assessment of telomere measurement in moose, and after controlling for sample storage duration, we show that moose from the sarmatic and montane ecoregions had shorter RTL than moose from the boreal ecoregion. These findings align with our hypothesis that moose in ecoregions encompassing higher levels of anthropogenic and environmental stress would have significantly shorter RTLs.

The characteristics of the two ecoregions with shorter RTLs vary substantially, and therefore reflect potentially different stressors that drive variation in RTL. The montane ecoregion is characterized by relatively lower land use intensity (Fig. 1), deep and extensive snow cover limiting locomotion ability, and low forage availability during winter [28]. All of which are likely responsible for elevated metabolic expenditure resulting in shorter RTL of animals in the montane ecoregion [29]. Conversely, the sarmatic ecoregion is characterized by higher moose hunting pressure from humans, competition with sympatric ungulate species and a higher land use intensity through human population density, traffic infrastructure (barriers and direct stress) and forestry activities [26,27,28]. The combined effect of these factors likely contributes to the shorter RTLs in this ecoregion. In addition to these anthropogenic stressors, mean annual Tc in the sarmatic ecoregion is substantially higher than the suggested upper critical temperature of 0 °C during winter, where moose were observed to experience increased metabolic rates and behavioural adaptation, such as altered habitat use and activity patterns [22]. When Tc is corrected to reflect actual ambient temperature experienced by moose in the southern ecoregion (by a conservative mean of 7.2 °C [31]), animals are on average exposed to temperatures exceeding their thermoneutral zone by approximately 7 °C during winter. This finding emphasises the concerns that moose in the southern limit of their range are heat stressed during winter [32] (Singh N. J. personal communication). Ultimately, chronic thermal stress [14] and trade-offs influenced by selection of suboptimal habitats [33, 34] (Singh N. J. personal communication) may therefore contribute towards determining RTL. Pathogen prevalence is also higher at lower latitudes with warmer climate [35, 36] and Beirne et al. [19] have demonstrated that European badgers Meles meles exhibit higher telomere attrition rate post infection with bovine tuberculosis. In accordance with our results, Spong et al. [21] have demonstrated that, moose hair cortisol levels—a shorter-term stress proxy than RTL—were higher in southern Sweden than in north.

RTLs of animals in the boreal ecoregion were longer compared to the other two regions. This can be attributed to a number of factors. First, the boreal region is generally more homogenous in vegetation, dominated by conifers that are interspersed with deciduous species. Commercial forestry is the main form of land use in this region, characterized by large tracts of monoculture and clear cutting being the most common method of timber harvest. Moose are known to prefer clear cuts and young pine forest < 5 m in height [37, 38]. Secondly, the proportion of migratory moose is higher in this region [28, 39], which allows the population to evade stressful periods of low food availability and deep snow, and provides food access all year round. Thirdly, the year round- availability of food through conifers being green, reduces starvation related stress. Fohringer et al. [30] identified several metabolites linked to high metabolic expenditure (e.g., several amino acids and ketone bodies) in moose in the corresponding montane area, while animals in the boreal region did not show elevated concentrations of such biomarkers that indicated starvation responses due to limiting winter diets. Moose in the montane region were observed to have a lower propensity to migrate, move shorter distances and have smaller seasonal home ranges compared to those in the boreal region [28]. This reduced migratory propensity and relatively higher and prolonged exposure to environmental stressors and a lack of abundant winter forage likely causes a higher chronic stress. The fact that all our captures were carried out during peak winter suggests that animals do not evade the environmental stressors experienced in this region at least during this period of limited browse availability [30]. Shorter RTL was also determined for roe deer [18] experiencing poor environmental conditions compared to a population in less harsh environments. Similarly, Hoelzl et al. [40] detected shorter RTL in edible dormice Glis glis that were not provided food ad libitum compared to individuals that were, and suggested that forage availability could be a major factor in determining telomere length in a wild species subject to highly variable resource availability.

The lack of significant results in relation to moose RTL and age in this study could be related to the fact that only animals in good body condition with the vast majority past their developmental phase, i.e., adults, were captured (Additional file 1: Figs. S1 and S2). Adult vertebrates (beyond significant additional growth) were shown to exhibit less variation in RTL than during the developmental (growing) phase [9, 13, 41]. Changes of RTL with age might, therefore, be less pronounced in adult individuals, such as those included in this study (aligning with the results of Wilbourn et al. [18] and Fairlie et al. [41] who reasoned that a selective disappearance of individuals with short telomeres increases average RTL with age in wild mammals). The onset of cellular and reproductive senescence effects in moose has been observed after the age of 10 for males [42] and 12 for females [43], however the management strategy of maintaining a moose population in prime condition, to maximise the number of individuals that can be hunted, means that few(er) individuals achieve ages at which senescence occurs. Moose management strategies in Sweden may therefore also partly explain the absence of a relationship between age and RTL. To better understand the role of animal age in telomere dynamics, individuals of all age groups would have to be examined, ideally in a longitudinal experiment [41]. Despite having observed insignificant changes of RTL with age, variation of RTL withing age groups was high and could be driven by regional effects, that may be attributed to differing degrees of environmental stress exposure and/or genetic differences. Moreover, sex was shown to not be a significant predictor for RTL in our study, which is also in line with other studies performed on free-ranging mammals (reviewed by [18, 19, 44]), but see, for instance, Watson et al. having found sex-differences in wild Soay sheep Ovis aries [45].

Despite the known caveats in using mammalian blood as a source material for RTL quantification, most notably due to potential immune responses causing shifts of the leukocyte profile (see [46]), we were able to rely on this sample type by streamlining lab work and careful statistical examination of potential bias-inducing variables. We were therefore able to produce comparable results in line with several other studies that relied on leukocyte DNA [e.g., 11, 18, 19, 41, 45, 47]. The strong effect of storage time highlights that telomere studies should always control for this issue if varying storage periods cannot be avoided. Reichert et al. [48] found that the storage method of blood affected RTL, indicating that storage duration will also have an effect on RTL. The effect of storage duration did not impact our study as storage time was randomly distributed throughout ecoregions and the other variables. Our study was not able to investigate whether RTL is a suitable biomarker for age (in this species and in the developmental stage tested). Future studies may benefit from the inclusion of telomerase activity estimates as suggested by several authors (e.g., [49, 50]) in order to better understand the associations of telomere length with environmental variables in the examined study system and beyond.

Due to known genetic differences between moose in southern and northern Sweden [51,52,53], we cannot rule out potential population effects that might contribute to differing telomere length between northern and southern ecoregions. Our finding that differences in RTL were not consistent over a latitudinal gradient is in line with Kärkkäinen et al. [16], suggesting that regional variation of telomere length may mirror local environmental conditions and/or genetic differences. By measuring heritability and including more (known) populations in their analysis, future studies should account for the effects of population pedigree [15, 17] and between-population differences on RTL [54], thereby enabling the disentanglement of potential genetic differences from environmental conditions.

Conclusions

Animals that are highly adaptable to land use change likely face environmental constraints beyond high land use intensity that lead to an accumulation of stressors driving chronic stress and ultimately RTL. Increased encroachment via the accumulation and extension of different forms of land use and impacts of accelerated climate change at northern latitudes can limit the potential of animals to evade stressful environmental conditions via, for example, migration and will likely exacerbate metabolic demand and negative consequences on animal health. Our study emphasises that it is crucial to consider distinct biogeographic scales that encompass cumulative impacts affecting organisms holistically. Future analysis of chronic stress effects in free-ranging species should focus on the continuous resampling of cohorts of animals to understand inter and intra-individual telomere dynamics in wild animals at the life history scale.

Methods

Study area

The study area covers the three major ecoregions in Sweden, i.e. montane birch forest and grasslands (‘montane’), boreal forest (‘boreal’) and sarmatic mixed forest (‘sarmatic’) (Fig. 1; [55]). Moose were captured in all three ecoregions. Ecoregion assignment was based on the winter capture location. The ‘montane’ ecoregion is characterized by high-elevation tundra vegetation and mountain birch Betula pubescens belt. Duration of snow cover in the capture area within the montane ecoregion lasts approximately 210 days and mean snow depth is approximately 45 cm. Accordingly, the duration of the vegetation-growing season lasts less than 100 days in this capture area. The ‘boreal’ ecoregion occupies the largest portion of Sweden’s biomes and is dominated by coniferous trees, interspersed with patches of deciduous forest. Despite mean snow depths in capture area within the boreal being similar to the montane ecoregion, snow cover lasts less than 190 days and the growing season is extended to approximately 120 days. The ‘sarmatic’ ecoregion in southern Sweden consists of a mixed conifer-broadleaf plant association. The climate in the two capture areas within the sarmatic ecoregions is comparably mild, ranging between 90 and 200 days of snow cover, 10–15 cm snow depth and a vegetation growing period of 180–220 days. For detailed habitat characterization of moose capture areas see [28].

Forestry is the prevailing form of land use occurring throughout northern Sweden except for the montane ecoregion, where forestry is unfeasible. Generally, forestry is expected to be more intensive in the southern study area, where more commercial tree species occur and turn-over rate is higher [56]. While the landscape in the south is forest dominated, it is also highly fragmented with clear cuts, settlements and agriculture. In contrast, agriculture and settlements occur only sporadically in the boreal capture area, and are virtually absent in the montane region.

While moose and roe deer Capreolus capreolus occur throughout Sweden, the distribution of red deer Cervus elaphus, fallow deer Dama dama and wild boar Sus scrofa is limited to southern Sweden. Hunting pressure remained relatively stable for moose and roe deer in recent decades [57] but southern latitudes are experiencing higher hunting pressure due to the higher diversity of sympatric game species [26, 27]. Prevalence of disease and parasites affecting moose health was also shown to be higher in southern Sweden moose populations compared to those in the north [24, 25].

Data collection and sampling

From 2009 to 2018, 153 samples of free-ranging adult moose were collected during winter (Jan–April) within the framework of the national moose research. Animals were immobilized from a helicopter via dart injection [58] with a CO2-powered rifle (Dan-Inject, Børkop, Denmark) with the drug combination of 4.5 mg etorphine (Captivon® 98 Etorphine HCl 9.8 mg/ml, Wildlife Pharmaceuticals (PTY) Ltd., 38 Wilkens St., Rocky Drift, White River, South Africa) and 50 mg xylazine (Xylased® 500 mg, Bioveta, a.s., Komenského 212, 68,323 Ivanovice na Hané, Česká Republica) [59,60,61]. During immobilization, all animals were fitted with GPS-collars including a temperature receiver (Vectronic-Aerospace, Berlin, Germany). Pregnancy status was determined by a veterinarian via rectal palpation in sarmatic and montane areas [62]. Age was estimated based on tooth wear [42, 63]. The number of calves at heel was determined visually from the helicopter. Blood samples were collected into 9 ml S-Monovette® Z-Gel dry collection tubes (Sarstedt, Germany) by jugular venipuncture of the fully immobilized animals. Collection tubes were processed according to the manufacturer’s instructions and stored at − 20 °C until DNA extraction. Data on GPS positions, ancillary Tc, sex, pregnancy status, and number of calves at heel was stored and accessed via the Wireless Remote Animal Monitoring (WRAM) database [64].

Since RTL was compared across ecoregions to evaluate chronic stress, we estimated the ambient temperature and level of human impacts experienced by moose in each of our sample areas based on their GPS tracks. Anthropogenic impacts on the landscape were measured using the global Human Modification map (gHM), which provides a cumulative measure of human modification of terrestrial lands across the globe at a 1-km resolution [29]. The mean gHM value was estimated for each individual based on one year of movement post (re-)capture. The individual movement track was standardized to eight locations per day and used to estimate the mean gHM value from the underlying raster. Moose generally show fidelity to their winter and summer ranges [28, 65] and we therefore assume that movements post-capture also reflect environmental conditions pre-capture. Similarly, mean annual Tc (as a proxy for ambient temperature; [31]) was based on GPS-locations post capture. We used R packages amt [66], SDLfilter [67], trajr [68], adehabitatLT and adehabitatHR [69] for GPS- and Tc-data preparation as well as raster [70] and rgdal [71] for gHM value extraction.

DNA extraction

Prior to DNA extraction, blood samples were thawed simultaneously at 4 °C for 4 h and the serum fraction and the gel layer were discarded. Per sample, approximately 40 mg of the coagulated blood fraction was incubated at 56 °C with 30 µl proteinase K (20 mg/ml; Qiagen, Germany) for one hour with repeated inverting and shaking of samples. A liquid state of the sample was attained by subsequent addition of 190 ml PBS pH 7.4 (2.7 mM KCl, 140 mM NaCl, 10 mM Phosphate), pipetting up and down and vortexing for 30 s. DNA extraction and purification were carried-out on a QIASymphony SP platform using the DSP DNA minikit (Qiagen, Germany) according to the manufacturer’s instructions. DNA yield and quality were quantified using a NanoDrop 2000 spectrometer (Thermo Fisher Scientific, USA; Additional file 2). Purified DNA was stored at − 20 °C for up to one month until further processing via qPCR, wherefore DNA was refrigerated at 4 °C for up to two days.

Relative telomere length (RTL) assessment

For measuring RTL, we used the real-time PCR approach [72] adapted for moose for the first time. A 54 bp fraction of the beta-lactoglobulin (BLG) gene was used as non-variable copy number (non-VCN) gene (tested for non-variability as described by Cawthon [73], Smith et al. [74] and Turbill et al. [75]). Primer sequences for the non-VCN gene were 5′- GCA GCT GTC TTT CAG GGA GAA TG -3′ (rt_BLG F) and 5′- CCC GAC ACT TAC CAT CGA TCT TG -3′ (rt_BLG R). Telomeric primer sequences were 5′-CGG TTT GTT TGG GTT TGG GTT TGG GTT TGG GTT TGG GTT-3′ (tel 1b) and 5′-GGC TTG CCT TAC CCT TAC CCT TAC CCT TAC CCT TAC CCT-3′ (tel 2b). Telomere and non-VCN gene PCRs were carried out in 9 separate runs with 20 ng DNA per reaction, 400 nmol l−1 of each primer combination (Tel1b/Tel2b or rt_BLG F/ rt_BLG R) in a final volume of 20 μl containing 10 μl of GoTaq® qPCR Master Mix (Promega). Samples were randomized per run based on sex, capture area, and capture year (see Additional file 2). PCR conditions for the telomere runs were 2 min at 95 °C followed by 40 cycles of 15 s at 95 °C, 20 s at 58 °C and 20 s at 72 °C. For non-VCN runs, PCR conditions were 2 min at 95 °C followed by 45 cycles of 15 s at 95 °C, 20 s at 58 °C and 20 s at 72 °C. A final melting step was included in each run with the temperature ramping from 65 to 95 °C in 1 °C steps. Each run contained a negative (non-template) control and two DNA extracts from moose livers as standard samples (to assess inter-run variability). All samples and controls were run in triplicates. Reactions were prepared using the Qiagility PCR robot (Qiagen, Germany) to minimize pipetting errors, and cycling was performed on a Rotorgene Q quantitative thermocycler (Qiagen, Germany). We used the software LinRegPCR (2012.0) [76] for analysis of non-baseline-corrected raw qPCR data, exported from the instrument. RTL was calculated using the method described by Ruijter et al. [77], modified by Hoelzl et al. [78].

The mean qPCR efficiency was calculated via the amplification plot method [76] which gives lower but more accurate estimates of efficiency than standard curve based methods [79, 80]. The estimates were 76.9% and 86.7% for the non-VCN gene and telomere reactions, respectively.

The intraclass correlation coefficient (ICC) was calculated as a measure of reliability within and between the runs, as suggested by Koo and Li [81]. ICC estimates and their 95% confident intervals for sample triplicates were calculated in R Version 3.5.2 [82]. Intra-rater ICC was calculated on all included data points based on a single-rating, absolute-agreement, 2-way mixed-effects model (ICC in library ‘irr’, [83]). Intra-assay ICC for Ct values for telomere assay was 0.85 [p < 0.0001, 95% (CI 0.82–0.88)] and for BLG 0.96 [p < 0.0001, 95% (CI 0.94–0.97)] showing a good and an excellent degree of reliability respectively. The ICC for inter-assay reliability was calculated for the standard samples based on a mean rating (k = 3), agreement, 2-way mixed-effects model. Interrater ICC for Ct values for the telomere assay was 0.94 [p < 0.0001, 95% (CI 0.54–1.0)] and for BLG 0.99 [p < 0.0001, 95% (CI 0.97–1.0)] showing an excellent degree of reliability for both. As all samples per individual were run on the same plate, inter-assay variability should have minimal effect on our longitudinal results.

The intra-assay coefficient of variation among replicates (intra-assay variation), an estimate of system precision, was further used to assess reproducibility. Mean intra-assay CV for Ct values of the non-VCN gene and telomere assay were 0.35 and 0.86%, respectively. The mean coefficient of variation among replicates (intra-assay variation) for Ct values of the non-VCN gene and telomere assay were 0.35 and 0.86%, respectively. Among runs (inter-assay variation), the mean coefficient of variation for Ct values of the non-VCN gene was 0.94%, and this was 2.76% for the telomere reaction.

Statistics

All statistical analyses were carried out using R 3.5.2 [82]. To explain variation among individuals in RTL, linear mixed effects models and postHoc test with Tukey adjustment for multiple comparisons were used (library lme4; [84], library emmeans). The initial model contained the two-way interaction between animal age (continuous) and sex, ecoregion, as well as storage time (to control for potential effects of sample storage duration, since time of storage have been associated with change in RTL [79]) as explanatory variables. To account for potential pseudoreplication among samples from recaptured individuals, individual ID was included as a random effect. Capture location was not included in the models as they highly correlate with ecoregions that animals were captured in. Due to the limited number of samples from recaptured individuals (n = 39) and the absence of recaptures in the boreal ecoregion, intra-individual telomere dynamics were not considered in our analysis. Additionally, we ran a model on a subset of the data containing only females (n = 68), accounting for the explanatory variables mentioned above, and we also included pregnancy status and number of calves at heel as additional variables. Model selection was carried out using the R function dredge (library MuMIn; [85]) which evaluates all possible candidate models, from which the best-fit model was selected based on AICc. Coefficients, their standard errors (s.e.), degrees of freedom (df), t and corresponding P-values of the models are reported using the lmerTest package [86]. All means are given together with their standard error.

Availability of data and materials

The datasets generated and analysed during the current study are available in the Dryad repository, https://doi.org/10.5061/dryad.44j0zpcd0

Abbreviations

BLG:

Beta-lactoglobulin

DNA:

Deoxyribonucleic acid

gHM:

Global human modification

GPS:

Global positioning system

HPA:

Hypothalamo-pituitary-adrenal

ICC:

Intraclass correlation coefficient

PCR:

Polymerase chain reaction

qPCR:

Quantitative polymerase chain reaction

RTL:

Relative telomere length

VCN:

Variable copy number

References

  1. Sih A, Ferrari MCO, Harris DJ. Evolution and behavioural responses to human-induced rapid environmental change. Evol Appl. 2011;4(2):367–87.

    PubMed  PubMed Central  Article  Google Scholar 

  2. Sapolsky RM, Romero LM, Munck AU. How do glucocorticoids influence stress responses? Integrating permissive, suppressive, stimulatory, and preparative actions. Endocr Rev. 2000;21:55–89.

    CAS  PubMed  Google Scholar 

  3. Romero LM. Physiological stress in ecology: lessons from biomedical research. Trends Ecol Evol. 2004;19(5):249–55.

    PubMed  Article  Google Scholar 

  4. Silverman MN, Pearce BD, Biron CA, Miller AH. Immune modulation of the hypothalamic-pituitary-adrenal (HPA) axis during viral infection. Viral Immunol. 2005;18(1):41–78.

    CAS  PubMed  Article  Google Scholar 

  5. Constantini D, Marasco V, Møller PA. A meta-analysis of glucocorticoids as modulators of oxidative stress in vertebrates. J Comp Physiol B. 2011;181:447–56.

    Google Scholar 

  6. Wang Z, Rhee DB, Lu J, Bohr CT, Zhou F, Vallabhaneni H, de Souza-Pinto NC, Liu Y. Characterization of oxidative guanine damage and repair in mammalian telomeres. PLoS Genet. 2010;6(5): e1000951.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  7. Shalev I, Entringer S, Wadhwa PD, Wolkowitz OM, Puterman E, Lin J, Epel ES. Stress and telomere biology: a lifespan perspective. Psychoneuroendocrinology. 2013;38:1835–42.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  8. Nussey DH, Baird D, Barrett E, Boner W, Fairlie J, Gemmell N, Hartmann N, Horn T, Haussmann M, Olsson M, Turbill C, Verhulst S, Zahn S, Monaghan P. Measuring telomere length and telomere dynamics in evolutionary biology and ecology. Methods Ecol Evol. 2014;5:299–310.

    PubMed  PubMed Central  Article  Google Scholar 

  9. Monaghan P. Organismal stress, telomeres and life histories. J Exp Biol. 2014;217:57–66.

    PubMed  Article  Google Scholar 

  10. Dugdale HL, Richardson DS. Heritability of telomere variation: it is all about the environment! Philos Trans R Soc B. 2018;373:20160450.

    Article  Google Scholar 

  11. Kotrschal A, Ilmonen P, Penn DJ. Stress impacts telomere dynamics. Biol Lett. 2007;3:128–30.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  12. Haussmann MF, Marchetto NM. Telomeres: linking stress and survival, ecology and evolution. Curr Biol. 2010;56(6):714–27.

    CAS  Google Scholar 

  13. Angelier F, Costantini D, Blévin P, Chastel O. Do glucocorticoids mediate the link between environmental conditions and telomere dynamics in wild vertebrates? A review. Gen Comp Endocrinol. 2018;256:99–111.

    CAS  PubMed  Article  Google Scholar 

  14. Debes PV, Visse M, Panda B, Ilmonen P, Vasemägi A. Is telomere length a molecular marker of past thermal stress in wild fish? Mol Ecol. 2016;25:5412–24.

    CAS  PubMed  Article  Google Scholar 

  15. Foley NM, Petit EJ, Brazier T, Finarelli JA, Hughes GM, Touzalin F, Puechmaille SJ, Teeling EC. Drivers of longitudinal telomere dynamics in a long-lived bat species, Myotis myotis. Mol Ecol. 2020;29:2963–77.

    CAS  PubMed  Article  Google Scholar 

  16. Kärkkäinen T, Laaksonen T, Burgess M, Cantarero A, Martínez-Padilla J, Potti J, Moreno J, Thomson RL, Tilgar V, Stier A. Population differences in the length and early-life dynamics of telomeres among European pied flycatchers. Mol Ecol. 2021;00:1–13.

    Google Scholar 

  17. Spurgin LG, Bebbington K, Fairfield EA, Hammers M, Komedeur J, Burke T, Dugdale HL, Richardson DS. Spatio-temporal variation in lifelong telomere dynamics in a long-term ecological study. J Anim Ecol. 2018;87:187–98.

    PubMed  Article  Google Scholar 

  18. Wilbourn RV, Froy H, McManus M-C, Cheynel L, Gaillard J-M, Gillot-Fromont E, Regis C, Rey B, Pellerin M, Jean-François Lemaître J-M, Nussey DH. Age-dependent associations between telomere length and environmental conditions in roe deer. Biol Lett. 2017;13:20170434.

    PubMed  PubMed Central  Article  Google Scholar 

  19. Beirne C, Delahay R, Hares M, Young A. Age-related declines and disease-associated variation in immune cell telomere length in a wild mammal. PLoS ONE. 2014;9(9): e108964.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  20. Olson DM, Dinerstein E, Wikramanayake ED, Burgess ND, Powell GVN, Underwood EC, D’amico JA, Itoua I, Strand HE, Morrison JC, Loucks CJ, Allnutt TF, Ricketts TH, Kura Y, Lamoreux JF, Wettengel WW, Hedao P, Kassem KR. Terrestrial ecoregions of the world: a new map of life on Earth. Bioscience. 2001;51:933–8.

    Article  Google Scholar 

  21. Spong G, Gould NP, Sahlén E, Cromsigt JPGM, Kindberg J, DePerno CS. Large-scale spatial variation of chronic stress signals in moose. PLoS ONE. 2020;15(1): e0225990.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  22. Renecker LA, Hudson RJ. Seasonal energy expenditures and thermoregulatory responses of moose. Can J Zool. 1986;64:322–7.

    Article  Google Scholar 

  23. McCann N, Moen R, Harris T. Warm-season heat stress in moose (Alces alces). Can J Zool. 2013;91:893–8.

    Article  Google Scholar 

  24. Malmsten J, Jakubek E-B, Björkman C. Prevalence of antibodies against Toxoplasma gondii and Neospora caninum in moose (Alces alces) and roe deer (Capreolus capreolus) in Sweden. Vet Parasitol. 2011;177(3–4):275–80.

    PubMed  Article  Google Scholar 

  25. Välimäki P, Madslien K, Malmsten J, Harkonen L, Harkonen S, Kaitala A, Kortet R, Laaksonen S, Mehl R, Redford L, Ylönen H, Ytrehus B. Fennoscandian distribution of an important parasite of cervids, the deer ked (Lipoptena cervi), revisited. Parasitol Res. 2010;107:117–25.

    PubMed  Article  Google Scholar 

  26. Dressel S, Ericsson G, Sandström C. Mapping social-ecological systems to understand the challenges underlying wildlife management. Environ Sci Policy. 2018;84:105–12.

    Article  Google Scholar 

  27. Neumann W, Widemo F, Singh NJ, Seiler A, Cromsigt JPGM. Strength of correlation between wildlife collision data and hunting bags varies among ungulate species and with management scale. Eur J Wildl Res. 2020;66:86.

    Article  Google Scholar 

  28. Allen AM, Månsson J, Sand H, Malmsten J, Ericsson G, Singh NJ. Scaling up movements: from individual space use to population patterns. Ecosphere. 2016;7(10): e01524.

    Article  Google Scholar 

  29. Kennedy CM, Oakleaf JR, Theobald DM, Baruch-Mordo S, Kiesecker J. Managing the middle: a shift in conservation priorities based on the global human modification gradient. Glob Change Biol. 2019;25:811–26.

    Article  Google Scholar 

  30. Fohringer C, Dudka I, Spitzer R, Stenbacka R, Rzhepishevska O, Cromsigt JPGM, Gröbner G, Ericsson E, Singh NJ. Integrating omics to characterize eco-physiological adaptations: how moose diet and metabolism differ across biogeographic zones. Ecol Evol. 2021;00:1–25.

    Google Scholar 

  31. Ericsson G, Dettki H, Neumann W, Arnemo JM, Singh NJ. Offset between GPS collar-recorded temperature in moose and ambient weather station data. Eur J Wildl Res. 2015;61:919–22.

    Article  Google Scholar 

  32. Monteith KL, Bleich VC, Stephenson TR, Pierce BM, Conner MM, Klaver RW, Bowyer RT. Timing of seasonal migration in mule deer: effects of climate, plant phenology, and life-history characteristics. Ecosphere. 2011;2(4):art47.

    Article  Google Scholar 

  33. Schwab FE, Pitt MD. Moose selection of canopy cover types related to operative temperature, forage, and snow depth. Can J Zool. 1991;69(12):3071–7.

    Article  Google Scholar 

  34. van Beest FM, Van Moorter B, Milner JM. Temperature-mediated habitat use and selection by a heat-sensitive northern ungulate. Anim Behav. 2012;84:723–35.

    Article  Google Scholar 

  35. Hofmeister E, Moede-Rogall G, Wesenberg K, Abbott R, Work T, Schuler K, Sleeman J, Winton J. Climate Change and Wildlife Health: Direct and Indirect Effects. U.S. Geological Survey Fact Sheet Fact Sheet 2010–3017. Madison, WI: U.S. Geological Survey, National Wildlife Health Center; 2010.

  36. Malmsten J, Dalin A-M, Moutailler S, Devillers E, Gondard M, Felton A. Vector-Borne zoonotic pathogens in Eurasian Moose (Alces alces alces). Vector-Borne Zoonotic Dis. 2019;19(3):207–11.

    PubMed  Article  Google Scholar 

  37. Edenius L, Månsson J, Hjortstråle T, Roberge J-M, Ericsson G. Effects of high-stump pre-commercial thinning on moose browsing and damage in young Scots pine stands. Scand J For Res. 2015;30(5):382–7.

    Google Scholar 

  38. Allen AM, Dorey A, Malmsten J, Edenius L, Ericsson G, Singh NJ. Habitat–performance relationships of a large mammal on a predator-free island dominated by humans. Ecol Evol. 2017;7:305–19.

    PubMed  Article  Google Scholar 

  39. Singh NJ, Börger L, Dettki H, Bunnefeld N, Ericsson G. From migration to nomadism: movement variability in a northern ungulate across its latitudinal range. Ecol Appl. 2012;22:2007–20.

    PubMed  Article  Google Scholar 

  40. Hoelzl F, Cornils JS, Smith S, Moodley Y, Ruf T. Telomere dynamics in free-living edible dormice (Glis glis): the impact of hibernation and food supply. J Exp Biol. 2016;219:2469–74.

    PubMed  PubMed Central  Article  Google Scholar 

  41. 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:140–8.

    CAS  PubMed  Article  Google Scholar 

  42. Ericsson G, Wallin K. Age-specific moose (Alces alces) mortality in a predator-free environment: evidence for senescence in females. ÉcoScience. 2001;8(2):157–63.

    Article  Google Scholar 

  43. Ericsson G, Wallin K, Ball JP, Broberg M. Age-related reproductive effort and senescence in free-ranging moose. Alces alces Ecology. 2001;82:1613–20.

    Article  Google Scholar 

  44. Barrett ELB, Richardson DS. Sex differences in telomeres and lifespan. Aging Cell. 2011;10:913–21.

    CAS  PubMed  Article  Google Scholar 

  45. Watson RL, Bird EJ, Underwood S, Wilbourn RV, Fairlie J, Watt K, Salvo-Chirnside E, Pilkington JG, Pemberton JM, McNeilly TN, Froy H, Nussey DH. Sex differences in leucocyte telomere length in a free-living mammal. Mol Ecol. 2017;26(12):3230–40.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  46. Davis AK, Maney DL, Maerz JC. The use of leukocyte profiles to measure stress in vertebrates: a review for ecologists. Funct Ecol. 2008;22:760–72.

    Article  Google Scholar 

  47. Seeker LA, Underwood SL, Wilbourn RV, Dorrens J, Froy H, Holland R, Ilska JJ, Psifidi A, Bagnall A, Whitelaw B, Coffey M, Banos G, Nussey DH. Telomere attrition rates are associated with weather conditions and predict productive lifespan in dairy cattle. Sci Rep. 2021;11:5589.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  48. Reichert S, Froy H, Boner W, Burg TM, Daunt F, Gillespie R, Griffiths K, Lewis S, Phillips RA, Nussey DH, Monaghan P. Telomere length measurement by qPCR in birds is affected by storage method of blood samples. Oecologia. 2017;184:341–50.

    PubMed  PubMed Central  Article  Google Scholar 

  49. Epel ES, Lin J, Dhabhar FS, Wolkowitz OM, Puterman E, Karan L, Blackburn EH. Dynamics of telomerase activity in response to acute psychological stress. Brain Behav Immun. 2010;24(4):531–9.

    CAS  PubMed  Article  Google Scholar 

  50. Criscuolo F, Smith S, Zahn S, Heidinger BJ, Haussmann MF. Experimental manipulation of telomere length: does it reveal a corner-stone role for telomerase in the natural variability of individual fitness? Philos Trans R Soc B. 2018;373:20160440.

    Article  CAS  Google Scholar 

  51. Niedziałkowska M, Hundertmark KJ, Jedrzejewska B, Sidorovich VE, Zalewska H, Veeroja R, Solberg EJ, Laaksonen S, Sand H, Solovyev VA, Sagaydak A, Tiainen J, Juškaitis R, Done G, Borodulin VA, Tulandin EA, Niedziałkowski K. The contemporary genetic pattern of European moose is shaped by postglacial recolonization, bottlenecks, and the geographical barrier of the Baltic Sea. Biol J Linn Soc. 2016;117(4):879–94.

    Article  Google Scholar 

  52. Wennerstrom L, Ryman N, Tison J-L, Hasslow A, Dalen L, Laikre L. Genetic landscape with sharp discontinuities shaped by complex demographic history in moose (Alces alces). J Mammal. 2016;97(1):1–13.

    Article  Google Scholar 

  53. Blåhed I-M, Königsson H, Ericsson G, Spong G. Discovery of SNPs for individual identification by reduced representation sequencing of moose (Alces alces). PLoS ONE. 2018;13(5): e0197364.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  54. Burraco P, Lucas PM, Salmón P. Telomeres in a spatial context: a tool for understanding ageing pattern variation in wild populations. Ecography. 2022;202:e05565.

    Google Scholar 

  55. Dinerstein E, Olson D, Joshi A, Vynne C, Burgess ND, Wikramanayake E, Hahn N, Palminteri S, Hedao P, Noss R, Hansen M, Locke H, Ellis EC, Jones B, Barber CV, Hayes R, Kormos C, Martin V, Crist E, Sechrest W, Price L, Baillie JEM, Weeden D, Suckling K, Davis C, Sizer N, Moore R, Thau D, Birch T, Potapov P, Turubanova S, Tyukavina A, de Souza N, Pintea L, Brito JC, Llewellyn OA, Miller AG, Patzelt A, Ghazanfar SA, Timberlake J, Klöser H, Shennan-Farpón Y, Kindt R, Lillesø J-PB, van Breugel P, Graudal L, Voge M, Al-Shammari KF, Saleem M. An ecoregion-based approach to protecting half the terrestrial realm. J BioScience. 2017;67(6):534–45.

    Article  Google Scholar 

  56. Fridman J, Walheim M. Amount, structure, and dynamics of dead wood on managed forestland in Sweden. Forest Ecol Manag. 2000;131:23–36.

    Article  Google Scholar 

  57. Bergström R, Danell K. Mer vilt i dag än för 50 år sedan. Vilt och fisk fakta 4. Vilt & Fisk Adaptiv Förvalting. Swedish University of Agricultural Science [In Swedish]; 2009.

  58. Kreeger TJ, Arnemo JM. Handbook of Wildlife Chemical Immobilization (5th ed.). Published by authors; 2018.

  59. Evans AL, Fahlman A, Ericsson G, Haga HA, Arnemo JM. Physiological evaluation of free-ranging moose (Alces alces) immobilized with etorphine-xylazine-acepromazine in Northern Sweden. Acta Vet Scand. 2012;54(77):1.

    Google Scholar 

  60. Lian M, Evans AL, Bertelsen MF, Fahlman A, Haga HA, Ericsson G, Arnemo JM. Improvement of arterial oxygenation in free-ranging moose (Alces alces) immobilized with etorphine-acepromazine-xylazine. Acta Vet Scand. 2014;56(51):1.

    Google Scholar 

  61. Græsli A-R, Thiel A, Fuchs B, Singh NJ, Stenbacka F, Ericsson G, Neumann W, Arnemo JM, Evans AL. Seasonal hypometabolism in female moose. Front Ecol Evol. 2020;8:107.

    Article  Google Scholar 

  62. Solberg EJ, Heim M, Arnemo JM, Sæther B-E, Os Ø. Does rectal palpation of pregnant moose cows affect pre-and neo-natal mortality of their calves? Alces. 2003;39:65–77.

    Google Scholar 

  63. Rolandsen CM, Solberg EJ, Heim M, Holmstrøm F, Solem MI, Sæther B-E. Accuracy and repeatability of moose (Alces alces) age as estimated from dental cement layers. Eur J Wildl Res. 2008;54:6–14.

    Article  Google Scholar 

  64. Dettki H, Brode M, Giles T, Hallgren J. Wireless remote animal monitoring (WRAM)–A new international database e-infrastructure for management and sharing of telemetry sensor data from fish and wildlife. In: Proceedings of the International Congress on Environmental Modelling and Software, Colorado; 2014.

  65. Ofstad E. Seasonal Variation in Site Fidelity of Moose (Alces alces) (Master thesis). Norges teknisk-naturvitenskapelige universitet, Fakultet for naturvitenskap og teknologi, Institutt for biologi. 2013. http://hdl.handle.net/11250/245241. Accessed 29 Jun 2021.

  66. Signer J, Fieberg J, Avgar T. Animal movement tools (amt): R package for managing tracking data and conducting habitat selection analyses. Ecol Evol. 2019;9:880–90.

    PubMed  PubMed Central  Article  Google Scholar 

  67. Shimada T, Jones R, Limpus C, Hamann M. Improving data retention and home range estimates by data-driven screening. Mar Ecol Prog Ser. 2012;457:171–80.

    Article  Google Scholar 

  68. McLean DJ, Volponi MAS. trajr: an R package for characterisation of animal trajectories. Ethology. 2018;124(6):440.

    Article  Google Scholar 

  69. Calenge C. The package adehabitat for the R software: tool for the analysis of space and habitat use by animals. Ecol Model. 2006;197:1035.

    Article  Google Scholar 

  70. Hijmans RJ, van Etten J. raster: Geographic analysis and modeling with raster data. R package version 2.0-12. 2012. http://CRAN.R-project.org/package=raster. Accessed 29 Jun 2021.

  71. Bivand R, Keitt T, Rowlingson B, Pebesma E, Sumner M, Hijmans M, Baston D, Rouault E, Warmerdam F, Ooms J, Rundel C. rgdal: Bindings for the ’Geospatial’ Data Abstraction Library. 2020. https://r-forge.r-project.org/projects/rgdal/. Accessed 29 Jun 2021.

  72. Thomas P, O’Callaghan NJ, Fenech M. Telomere length in white blood cells, buccal cells and brain tissue and its variation with ageing and Alzheimer’s disease. Mech Ageing Dev. 2008;129:183–90.

    CAS  PubMed  Article  Google Scholar 

  73. Cawthon RM. Telomere measurement by quantitative PCR. Nucleic Acids Res. 2002;30:E47.

    PubMed  PubMed Central  Article  Google Scholar 

  74. Smith S, Turbill C, Penn DJ. Chasing telomeres, not red herrings, in evolutionary ecology. Heredity. 2011;107:372–3.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  75. Turbill C, Smith S, Deimel C, Ruf T. Daily torpor is associated with telomere length change over winter in Djungarian hamsters. Biol Lett. 2012;8:304–7.

    PubMed  Article  Google Scholar 

  76. Ramakers C, Ruijter JM, Deprez RHL, Moorman AFM. Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci Lett. 2003;339:62–6.

    CAS  PubMed  Article  Google Scholar 

  77. Ruijter JM, Ramakers C, Hoogaars WMH, Karlen Y, Bakker O, van den Hoff MJB, Moorman AFM. Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data. Nucleic Acids Res. 2009;37:E45.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  78. Hoelzl F, Smith S, Cornils JS, Aydinonat D, Bieber C, Ruf T. Telomeres are elongated in older individuals in a hibernating rodent, the edible dormouse (Glis glis). Sci Rep. 2016;6:36856.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  79. Morinha F, Magalhães P, Blanco G. Standard guidelines for the publication of telomere qPCR results in evolutionary ecology. Mol Ecol Resour. 2020;20:635–48.

    Article  Google Scholar 

  80. Spießberger M, Hoelzl F, Smith S, Vetter S, Ruf T, Nowack J. The tarnished silver spoon? Trade-off between prenatal growth and telomere length in wild boar. J Evol Biol. 2022;35:81–90.

    PubMed  Article  CAS  Google Scholar 

  81. Koo TK, Li MY. A guideline of selecting and reporting intraclass correlation coefficients for reliability research. J Chiropr Med. 2016;15(2):155–63.

    PubMed  PubMed Central  Article  Google Scholar 

  82. R Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2018.

  83. Gamer M, Lemon J, Singh IFP. Irr: Various coefficients of interrater reliability and agreement. 2019. https://CRAN.R-project.org/package=irr. Accessed 31 May 2022.

  84. Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67(1):1–48.

    Article  Google Scholar 

  85. Bartoń K. MuMIn: Multi-Model Inference. R package version 1.43.17. 2020. https://CRAN.R-project.org/package=MuMIn. Accessed 29 Jun 2021.

  86. Kuznetsova A, Brockhoff PB, Christensen RHB. lmerTest package: tests in linear mixed effects models. J Stat Softw. 2017;82(13):1–26.

    Article  Google Scholar 

Download references

Acknowledgements

We thank Fredrik Stenbacka for assistance during fieldwork and in accessing the biobank. Moreover, we thank Helena Königsson at the genetics lab at Department of Wildlife, Fish, and Environmental Studies at the Swedish University of Agricultural Sciences and staff at the genetics lab Evolution at University of Veterinary Medicine Vienna for facilitating DNA extraction and RTL quantification, respectively.

Funding

Open access funding provided by Swedish University of Agricultural Sciences. This study was financed by the project “Resource Extraction and Sustainable Arctic Communities (REXSAC), which is funded by Nordforsk—a “Nordic Centre of Excellence”—(project number 76938) and NJS was partially supported by another Nordforsk programme REIGN—Reindeer Husbandry in a globalizing north. The Seth M Kempe Foundation funded travel and accommodation costs to carry out lab work at the Konrad Lorenz Institute of Ethology, University of Veterinary Medicine, Vienna. Consumables and running costs for lab work was covered by a grant from the Helge Ax:son Johnsons stiftelse (Grant number F18-0363).

Author information

Authors and Affiliations

Authors

Contributions

CF, GS, and SS conceived the ideas and designed methodology; GE authorised sample acquisition from the in-house biobank; CF collected additional samples; CF extracted DNA; CF and FH ran the qPCR assay and assessed telomere length under supervision of SS; CF and FH compiled data and analysed data and together with CC, AMA and NJS; CF and NJS led the writing of the manuscript; AMA, NJS and GE contributed with information regarding the study species and movement ecology; SS and FH contributed with information on telomere dynamics. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Christian Fohringer.

Ethics declarations

Ethics approval and consent to participate

All moose captures and handling are in line with the ethical permits: A124-05, A77-06, A116-09, A50-12 and A14-15, granted by the Swedish Animal Ethics Committee.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1

is the model selection table for models included in this manuscript. Figures S1–3 are showing variables that did not pass model selection.

Additional file 2.

DNA quality and sample distribution per qPCR run. [https://doi.org/10.5061/dryad.44j0zpcd0].

Additional file 3.

RTL dataframe. [https://doi.org/10.5061/dryad.44j0zpcd0].

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Fohringer, C., Hoelzl, F., Allen, A.M. et al. Large mammal telomere length variation across ecoregions. BMC Ecol Evo 22, 105 (2022). https://doi.org/10.1186/s12862-022-02050-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12862-022-02050-5

Keywords

  • Alces alces
  • Biomarker
  • Chronic stress
  • Human modification
  • Life history
  • Telomere associations