- 1Section of Preventive Medicine and Epidemiology, Department of Medicine, Boston University School of Medicine, Boston, MA, United States
- 2Department of Epidemiology and Center for Social Epidemiology and Population Health, School of Public Health, University of Michigan, Ann Arbor, MI, United States
- 3Survey Research Center, Institute for Social Research, University of Michigan, Ann Arbor, MI, United States
- 4Department of Biostatistics, University of Washington, Seattle, WA, United States
- 5Department of Biostatistics, School of Public Health, University of Michigan, Ann Arbor, MI, United States
- 6Department of Medicine, Division of Geriatrics, UCLA, Los Angeles, CA, United States
Background: Short telomere length (TL) and telomere attrition (TA) have been associated with age-related diseases.
Objective: We assessed whether a genetic risk score for short TL (GRS-TL) combining seven TL-associated genetic variants identified in a European-ancestry genome-wide association study (GWAS) was associated with TL and TA over 10 years.
Methods: Relative TL (T/S ratio) was measured by the quantitative polymerase chain reaction method for a sample of white, African American, and Hispanic participants, who attended Exam 1 and/or 5 of the Multi-Ethnic Study of Atherosclerosis (MESA). Our final sample included 1,227 participants for the TL analysis and 1,138 for the TA analysis. Participants were 45–84 years at Exam 1. We used a linear mixed effects model and adjusted for age, sex, and population structure. Models were stratified by race/ethnicity.
Results: In the TL analysis, higher GRS-TL significantly predicted shorter TL (estimates = -0.18 [S.E. = 0.08], p = 0.02 for white; -0.18 [0.07], p < 0.01 for African American; and -0.13 [0.05], p = 0.02 for Hispanic) in fully adjusted models. In the TA analysis, no association between GRS-TL and TA over 10 years was found.
Conclusion: Although GRS-TL was developed in European-ancestry populations, it was significantly associated with TL (but not TA) in all three race/ethnic groups examined.
Introduction
Telomeric DNA shortens with each division of the cell (Fuster and Andres, 2006; Martinez and Blasco, 2011), and telomere shortening occurs throughout life. (Frenck et al., 1998). Telomere length (TL) has been proposed as a biomarker of aging-related diseases. (De Meyer et al., 2018). In vitro, telomere shortening correlates with cellular senescence and apoptosis. (Rudolph et al., 1999). In epidemiological studies, the attrition of telomeres has been linked to cardiovascular disease (CVD). Studies have found that shorter telomere length is associated with increased odds of mortality, (Cawthon et al., 2003; Epel et al., 2008; Kimura et al., 2008; Weischer et al., 2012), psychological stress, (Epel et al., 2004), and adverse socioeconomic position. (Needham et al., 2013; Geronimus et al., 2015).
TL shows important variability among individuals. (Frenck et al., 1998; Okuda et al., 2002). At birth, TL varies from 5,000 to 15,000 base pairs, (Sanders and Newman, 2013), with a fast rate of telomere attrition (TA) during growth and development, but by the age of 20, the rate of TA declines gradually. (Aubert and Lansdorp, 2008; Hjelmborg et al., 2015). Men and older people tend to have shorter TL than women and younger people. (Muezzinler et al., 2013). These individual differences appear to result from both environmental and genetic factors. (Okuda et al., 2002). A meta-analysis in over 19,000 subjects reported an estimate of TL heritability of 0.70, (Broer et al., 2013), and a genome-wide association study (GWAS) meta-analysis of more than 37,000 individuals (from 15 cohorts in Europe and Australia, all of European descent) identified seven loci associated with mean TL. (Codd, Nelson et al., 2013).
Although other TL GRSs have been estimated using TL measured by DNA sequencing, our GRS, is the only one that examined TL measured by quantitative polymerase chain reaction method (qPCR) which is the same phenotype as in our study. To our knowledge, only one previous study has examined a genetic risk score for short TL (GRS-TL) as a predictor of TL in a multi-ethnic sample that included whites, African Americans, and Hispanics. (Hamad et al., 2016). However, there is still scant literature about longitudinal measurements of TL. The purpose of the current study was to assess the association of the GRS-TL with both TL and TA over 10 years in the Multi-Ethnic Study of Atherosclerosis (MESA).
Materials and methods
Study population
Our sample was derived from the Multi-Ethnic Study of Atherosclerosis (MESA), (Bild et al., 2002), a prospective cohort study of the natural history of CVD. Between 2000 and 2002, MESA recruited 6,814 white, African American, Chinese American, and Hispanic participants aged 45–84 in the United States. Our analytical sample for the TL analysis included 1,162 participants with TL measurement at Exam 1, 1,203 participants with TL measurement at Exam 5, and 1,227 participants with TL measurement at Exams 1 and/or 5 (pooled sample). The analytic sample for the TA analysis included 1,138 participants who had TL measurement at both Exams 1 and 5.
Telomere length
Stored blood samples from MESA Exams 1 and 5 were used to assess TL. Extensive details of the protocol for blood collection and DNA extraction are described elsewhere. (Carroll et al., 2019). DNA was extracted and purified using sodium dodecyl sulfate cell lysis and salt precipitation methods. DNA was of high quality and high molecular weight. TL assays were performed at the Blackburn Laboratory at the University of California, San Francisco, using the quantitative polymerase chain reaction (qPCR) method to measure TL relative to standard reference DNA (T/S ratio). (Cawthon, 2002; Lin et al., 2010). Each sample was assayed 3 times on 3 different days. (Carroll et al., 2019; Needham et al., 2019). Briefly, sample plates were assayed in groups of three 96 well plates. Any assay run with 8 or more invalid control wells was considered a failed run. The mean of the T/S values was calculated, and values that deviated most from the mean were marked as a potential outlier. The mean of the T/S value was calculated without the potential outlier (99.8% of all samples was free of outliers). DNA samples were coded, and the lab personnel was blind to all other measurements of study participants. The interassay coefficient of variation was 2.9% ± 2.1%. (Carroll et al., 2019).
Genetic risk score for short telomere length
To create our GRS-TL, we used the effect estimates (betas) reported in the GWAS meta-analysis (Codd et al., 2013) that identified seven loci associated with mean TL including ACYP2 (rs11125529), TERC (rs10936599), NAF1 (rs7675998), TERT (rs2736100), OBFC1 (rs9420907), ZNF208 (rs8105767) and RTEL1 (rs755017). These 7 single nucleotide polymorphisms (SNPs) showed a consistent association with short TL at the genome-wide significance (P < 5 × 10−8) among individuals of European ancestry. Each SNP explained a very small proportion of the total variance in TL (0.08–0.36%), but the variance in TL explained by the seven SNPs jointly was greater. (Codd et al., 2013). MESA measured genotype data was available from genetic data collected between 2000 and 2013. For SNPs not directly genotyped in MESA (rs10936599, rs7675998, and rs755017), we used imputed genotype by obtaining SNP dosages from 1000 Genomes Project imputed data. The imputation quality was excellent (INFO score >0.92) for all imputed SNPs.
We calculated the GRS-TL by summing the number of risk alleles (from 0 to 2) in the seven loci associated with mean TL and weighted by their effect size in the GWAS meta-analysis. (Codd et al., 2013). A higher GRS-TL indicates a greater genetic risk for short TL.
We created histograms and violin plots to show the distribution of the GRS-TL in the whole sample and in different race/ethnic groups. For violin plots, TL was dichotomized at the 25th percentile, and participants in the bottom quartile (shortest TL) were compared to those in the top three quartiles (longest TL). We plotted the distribution of the GRS-TL in the bottom TL quartile vs. the other quartiles within each race/ethnic group.
Covariates
Our models were adjusted for study exam (Exam 1 vs. 5) in the TL analysis and for time since Exam 1 in the TA analysis. We also accounted for age in years (at Exam 1 or Exam 5 for the TL analysis and age at Exam 1 for the TA analysis), sex, and genetic principal components (PCs) calculated by MESA that were extracted from genome-wide data. Adding PCs as covariates in regression models accounts for population stratification that can cause type I error inflation and/or loss of power in genetic association studies. (Zhang et al., 2003). We used the PCs calculated within each ethnic group separately for primary analysis and those in the total MESA sample (all race/ethnicities combined) for secondary analysis. The first 3 PCs explained about 78% of the observed variation in the whole sample.
Statistical analyses
We conducted descriptive analyses to examine the average TL at Exam 1, Exam 5 and change over the 10-year follow-up from Exam 1 to Exam 5 for the full sample and by demographics.
To examine the associations of GRS-TL with TL and 10-year TA in our multi-ethnic sample, our primary analyses were race/ethnicity stratified, and secondarily we performed analyses with all race/ethnicities combined for both the TL and TA analyses.
TL analysis
We used a linear mixed-effects model (LMEM) that included all valid observations from Exams 1 and 5 to examine the association between the GRS-TL and TL comparing whites, African Americans, and Hispanics. The model included TL as the outcome and GRS-TL as the main exposure. GRS-TL and covariates were centered at their population means.
Models were adjusted for an exam indicator (reference: Exam 1), age at Exam 1 or 5, sex (reference: female), and PCs specific for each race/ethnic group. An individual-level random intercept was included in the model to account for the within-person correlation of repeated telomere measures (for those who have telomere measures at both exams). Robust standard errors were reported.
In a secondary analysis, we combined all race/ethnic groups and additionally included in the model race/ethnicity (dummy variables for African American and Hispanic with white as the reference category), the interactions between GRS-TL and race/ethnicity, and PCs combined for all race/ethnic groups.
TA analysis
We used LMEM to examine the association between the GRS-TL and 10-year TA. Follow-up time since Exam 1 variable was centered on the individual’s average follow-up time to obtain the estimate of the intra-individual change in TL over time. All other covariates were centered on the population mean. The model included follow-up time since Exam 1 and the interaction between follow-up time and the GRS-TL. The coefficient for the interaction term of follow-up time and GRS-TL represents the effect of the GRS-TL on the annual change in TL. The model also included interactions of time with age at Exam 1, sex, and PCs. An individual-level random intercept was included to account for the intra-individual variability of the baseline TL and the change in TL over time. Robust standard errors were reported. For the TA analysis, our primary analysis was also race/ethnicity stratified. In a secondary analysis, we combined all race/ethnic groups. We included three-way interaction terms of the follow-up time, GRS-TL, and race/ethnicity dummy variables. To account for population stratification, we included the first three race/ethnicity-specific PCs for the primary race-stratified analysis and the first three PCs for the combined sample for the secondary combined-race analysis.
Results
The sample included white (n = 323, 28%), African American (n = 330, 28%), and Hispanic participants (n = 509, 44%) at Exam 1 (Table 1) with an average age of 61 years at Exam 1 and 70 years at Exam 5.
We calculated the effect allele frequencies of the SNPs in each race/ethnic group and in all race/ethnicity groups combined and compared them to the effect allele frequencies from the GWAS (Codd et al., 2013), from which we obtained the effect size of the SNPs used to create our GRS-TL (Table 2). The calculated GRS-TL was slightly skewed (Shapiro Wilcoxon p = 0.07 and Kolmogorov-Smirnov p < 0.01) (Supplementary Figure S1) with a mean (SD) = 0.51 (0.13). Distribution of the GRS-TL for the TA sample and by race/ethnicity is shown in Supplementary Table S1. Hispanics and whites had a greater GRS-TL mean (0.55 and 0.54, respectively) than African American participants (0.42). Plots of GRS-TL and TL (T/S ratio) at Exams 1 and 5 by race/ethnicity and for the whole sample combined are presented in Figure 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021051/fragi-03-1021051-HTML/image_m/fragi-03-1021051-g001.jpg)
FIGURE 1. Plots of GRS-TL and TL (T/S ratio) at Exams 1 and 5 by race/ethnicity (white, African American, and Hispanic) and for the whole sample combined. Telomere length as a function of the GRS-TL in Exam 1 (red circles, n = 1162) and Exam 5 (blue circles, n = 1203). Telomere length is plotted as log T/S ratio. Panel (A): GRS-TL and log-transformed TL for whites. Panel (B): GRS-TL and log-transformed TL for Blacks. Panel (C): GRS-TL and log-transformed TL for Hispanics. Panel (D): GRS-TL and log-transformed TL for the whole sample.
As expected, we observed a decreasing TL over time. At Exam 1, the average TL (T/S ratio) was 0.92 (SD = 0.21), and at Exam 5 TL was 0.71 (SD = 0.14). Ten-year TA was –0.22 (SD = 0.19). The 10-year TA by demographics is shown in Table 3. Older participants and men had significantly shorter TL at both Exams 1 and 5, but the 10-year TA was not significantly different by age or sex. African Americans had significantly shorter TL than whites at Exam 1, but TL was not significantly different by race/ethnicity at Exam 5; however, 10-year TA was greater for whites and Hispanics than for African Americans. In the whole sample, 83.3% of participants experienced at least a 5% decrease in TL between Exams 1 and 5, whereas only 6.4% experienced at least a 5% increase in TL.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021051/fragi-03-1021051-HTML/image_m/fragi-03-1021051-t003.jpg)
TABLE 3. Mean Telomere Length (TL, T/S Ratio) and 10-year Telomere Attrition (10-year TA) for the Full Sample and by demographics.
Figure 2 shows violin plots of the GRS-TL distributions at Exam 1 (Panel A) and Exam 5 (Panel B) stratified by race/ethnicity, comparing those in the lower 25th percentile of TL (shortest) to those in the upper 75th percentile (longest). Although median GRS-TL was higher in those with the lower 25th percentile of TL for most race/ethnic groups at both exams, the GRS-TL distributions between the lower 25th percentile and upper 75th percentile did have substantial overlap.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021051/fragi-03-1021051-HTML/image_m/fragi-03-1021051-g002.jpg)
FIGURE 2. Violin Plots of the GRS-TL distributions stratified by race/ethnicity. Violin Plots of the GRS-TL distributions including a box plot with median, and quartiles stratified by race/ethnicity for those at or below the 25th percentile of TL (shorter TL) and those above the 25th percentile of TL (longer TL). Panel (A) GRS by TL at 25 Percentile by Race-Ethnicity at Exam 1; Panel (B) GRS by TL at 25 Percentile by Race-Ethnicity at Exam 5.
TL analysis
In stratified analysis by race/ethnicity, higher GRS-TL was significantly associated with shorter TL among whites, African Americans, and Hispanics. A one-unit increase in the GRS-TL from its mean was associated with a 0.18 unit decrease in T/S ratio (b = –0.18; SE = 0.08; p ≤ 0.02 and b = –0.18; SE = 0.07; p ≤ 0.01) in whites and African Americans, respectively, and a 0.13 unit decrease in T/S ratio (b = –0.13; SE = 0.05; p = 0.02) in Hispanics, after adjustment for exam year, age at exam, sex, and the first three PCs specific for each race/ethnicity (Table 4). In a secondary analysis, we combined all races/ethnicities. The combined sample also showed a significant inverse association between the GRS-TL and TL (estimate = –0.16 [SE = 0.04], p < 0.0001) after adjustment for exam year, age, sex, race/ethnicity, and the three PCs of the combined sample. Neither the interaction term between the GRS-TL and African American or GRS-TL and Hispanic were significant (Supplementary Table S2), suggesting that the association between GRS-TL and TL did not differ by race/ethnicity.
TA analysis
In stratified analysis by race/ethnicity, GRS-TL was not significantly associated with 10-year TA for any race/ethnic group (Table 5). Results were similar in the combined sample (Supplementary Table S3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021051/fragi-03-1021051-HTML/image_m/fragi-03-1021051-t005.jpg)
TABLE 5. Analysis of the association of GRS-TL with 10-year Telomere Attrition (10-year TA) stratified by race/ethnicity.
Discussion
The goals of the current study were 1) to examine the association between genetic predisposition to shorter TL, measured by a weighted GRS of seven SNPs, and TL and 2) to examine the association between the GRS-TL and within-person TA over 10 years. In our study, we observed an association between the GRS-TL and TL in each of the race/ethnic groups (white, African American, and Hispanic). When we analyzed the whole sample (all three race/ethnicity groups combined), the association between the GRS-TL and TL was also significant in the fully adjusted model. However, the GRS-TL was not associated with the risk of 10-year TA.
The heritability (h2) estimate of TL in a meta-analysis of six cohorts, all of European descent, was 0.70 (95% confidence interval 0.64–0.76). (Broer et al., 2013). Other cohorts have showed greater variation of h2 of TL. (Dugdale and Richardson, 2018). The h2 of TA has been less investigated, and it seems lower than the h2 of TL. In a study using same-sex twin data, h2 of TA was between 24% and 32% with low or no shared environmental effect. (Hjelmborg et al., 2015). The authors of this study suggested that heritability and early life environment may determine the rate of TA early in life, before adulthood, and this could be related to why TA is much faster in children than in adults. (Hjelmborg et al., 2015).
GRS for TL
In the current study, we observed a significant inverse association between a weighted GRS of seven SNPs for shorter TL, which was derived from a GWAS in samples of European descent, with TL in African Americans, Hispanics, and whites. In the combined sample, there was also a significant inverse association between the GRS-TL and TL in the fully adjusted model, and the association between GRS-TL and TL did not differ by race/ethnicity.
Replication studies in different race/ethnic backgrounds have showed mixed results. A GWAS meta-analysis in mostly white and a small group of African Americans participants identified associations of TL with two loci: OBFC1(rs4387287, rs9419958, and rs2487999) and CXCR4 (rs4452212 and 4954585). Further, evidence of replication was found for rs4387287 and rs9419958 in whites and a combined sample of whites and African Americans, but not in African Americans alone. The other SNPs were not replicated in any group (whites, African Americans, or both races combined). (Levy et al., 2010). Another study sought replication in only two out of their seven SNPs associated with short TL in a pooled cohort of participants of European descent. (Codd et al., 2013).
In a sample of 11,934 older adults from the US Health and Retirement Study (76% white, and 24% African Americans, Hispanics, or other race/ethnicity), a GRS-TL derived from the same GWAS we used in the current study was evaluated. In that study, compared to whites, African Americans were significantly more likely to have a lower GRS-TL, but Hispanics were more likely to have a higher GRS-TL. However, after adjustment for population structure, these associations were no longer significant. (Hamad et al., 2016). Hamad et al. also found differences in the association of individual TL-associated SNPs and TL by race/ethnicity for some SNPs. (Hamad et al., 2016).
Strengths and limitations
The current study has several strengths. First, this is a large community-based study that includes three race/ethnic groups (whites, African Americans, and Hispanics), with an overrepresentation of the last group as opposed to what is observed in other studies, whose percent of Hispanics were small. (Hamad et al., 2016). Second, the average inter-assay coefficient of variation for the telomere length assay was only 2.9%. Third, our GRS was estimated on qPCR, which is a better choice for our current study, but future work will examine sequencing-estimated scores. (Cawthon, 2002; Codd et al., 2013). Fourth, our analytical approach included LMEM, which allowed us to use our data more efficiently. Finally, the period between the first and second telomere length measurements was 10 years, which is about the time needed between measurements to observe true change in TL in adults, as reported in previous research. (Steenstrup et al., 2013). Whereas most participants in our sample experienced a decrease in TL between Exams 1 and 5, 6.4% experienced an increase in TL of at least 5%. Although we cannot rule out the possibility that measurement error accounts for this finding, a simulation study by Bateson and Nettle (Bateson and Nettle, 2017) indicates that scenarios in which true lengthening occurs provide a better fit for the empirical data than scenarios in which no lengthening occurs. In addition, prior experimental studies have shown that non-technical factors, such as weight loss (García-Calzón et al., 2014; Carulli et al., 2016) may contribute to telomere lengthening.
Our study also has several limitations. First, to create our GRS-TL, we used 7 SNPs derived from a meta-analysis GWAS in samples of European descent, which almost certainly has greater predictive utility for whites than for African Americans or Hispanics. Second, the genetic factors we used to create our GRS-TL may influence only telomere length, but the rate of TA may be influenced by different genetic factors, as mentioned by others. (Levy et al., 2010). Sanders et al. suggest that TL indicates somatic growth before the end of puberty and cellular senescence and oxidative stress after puberty. (Sanders and Newman, 2013). TL may be affected by DNA lost during replication and by oxidative stress, (Proctor and Kirkwood, 2002), but also by several environmental factors such as smoking and low socioeconomic status. (Needham et al., 2013). The genetic and environment correlates of TL may differ from the genetic and environment correlates of TA. Third, TA is difficult to study, particularly when using PCR, because it is measured with error, which could explain why the GRS-TL was not significant in the TA analysis. (Nettle et al., 2021). Fourth, the period between the first and second measurement may not be long enough to observe a meaningful change in TL given that the average annual rate of TL shortening in adults is slower than it is before the 20 years of age. Last, the range of age of our sample did not include very young participants where we could have observed greater TA.
Conclusion
In conclusion, our study focused on the ability to evaluate whether SNPs identified in the European GWAS were associated with TL and 10-year TA in a multi-ethnic sample. In the TL analysis, we observed an association between GRS-TL and TL in each of the race/ethnic groups studied, with a similar degree of association in whites and African Americans and a little lower among Hispanics. The association was also present and significant when we combined the race/ethnic groups, suggesting that the GRS-TL is associated with TL well in the three race/ethnic groups. However, we did not find an association between our GRS-TL and 10-year TA and this could be explained because the SNPs used in our GRS-TL reflect variants of telomere length but not of telomere shortening over some time. Larger longitudinal studies that include a wide range of age from birth to old age may be needed to better understand the dynamics of telomere length.
Data availability statement
The data analyzed in this study is subject to the following licenses/restrictions: “The MESA data used to support the findings of this study are available at the Multi-ethnic Study of Atherosclerosis website https://www.mesa-nhlbi.org at request after following procedures”. Requests to access these datasets should be directed to Belinda Needham, needhamb@umich.edu.
Author contributions
CC-D was the lead author, responsible for data analysis and writing the manuscript. JAS assisted in conceptualization and analysis. WZ and XW assisted in data analysis. BM and TS gave detailed comments on the study. BLN was the senior author, conceptualized the study, and obtained funding for this study. All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
MESA and the MESA SHARe project are conducted and supported by the National Heart, Lung, and Blood Institute (NHLBI) in collaboration with MESA investigators. Support for MESA is provided by contracts 75N92020D00001, HHSN268201500003I, N01-HC-95159, 75N92020D00005, N01-HC-95160, 75N92020D00002, N01-HC-95161, 75N92020D00003, N01-HC-95162, 75N92020D00006, N01-HC-95163, 75N92020D00004, N01-HC-95164, 75N92020D00007, N01-HC-95165, N01-HC-95166, N01-HC-95167, N01-HC-95168, N01-HC-95169, UL1-TR-000040, UL1-TR-001079, UL1-TR-001420, UL1-TR-001881, and DK063491. The authors thank the MESA participants and the MESA investigators and staff for their valuable contributions. A full list of participating MESA investigators and institutions can be found at http://www.mesa-nhlbi.org.
Funding for SHARe genotyping was provided by NHLBI Contract N02-HL-64278. Genotyping was performed at Affymetrix (Santa Clara, California, United State) and the Broad Institute of Harvard and MIT (Boston, Massachusetts, United States) using the Affymetrix Genome-Wide Human SNP Array 6.0. The MESA Stress Study was supported by R01HL076831 and R01HL101161 (PI: Ana Diez Roux). Work was also partially supported by funding from the UCLA Older Americans Independence Center, NIH/NIA Grant P30AG028748, and the USC/UCLA Biodemography Center through a P30 grant from the NIA (P30AG017265). CC-D was supported by the Multidisciplinary Training Program (T32) in Cardiovascular Epidemiology (5T32HL125232).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fragi.2022.1021051/full#supplementary-material
References
Aubert, G., and Lansdorp, P. M. (2008). Telomeres and aging. Physiol. Rev. 88 (2), 557–579. doi:10.1152/physrev.00026.2007
Bateson, M., and Nettle, D. (2017). The telomere lengthening conundrum – It could be biology. Aging Cell 16 (2), 312–319. doi:10.1111/acel.12555
Bild, D. E., Bluemke, D. A., Burke, G. L., Detrano, R., Diez Roux, A. V., Folsom, A. R., et al. (2002). Multi-ethnic study of Atherosclerosis: objectives and design. Am. J. Epidemiol. 156 (9), 871–881. doi:10.1093/aje/kwf113
Broer, L., Codd, V., Nyholt, D. R., Deelen, J., Mangino, M., Willemsen, G., et al. (2013). Meta-analysis of telomere length in 19, 713 subjects reveals high heritability, stronger maternal inheritance and a paternal age effect. Eur. J. Hum. Genet. 21 (10), 1163–1168. doi:10.1038/ejhg.2012.303
Carroll, J. E., Irwin, M. R., Seeman, T. E., Diez-Roux, A. V., Prather, A. A., Olmstead, R., et al. (2019). Obstructive sleep apnea, nighttime arousals, and leukocyte telomere length: the multi-ethnic study of Atherosclerosis. Sleep 42 (7), zsz089. doi:10.1093/sleep/zsz089
Carulli, L., Anzivino, C., Baldelli, E., Zenobii, M. F., Rocchi, M. B. L., and Bertolotti, M. (2016). Telomere length elongation after weight loss intervention in obese adults. Mol. Genet. Metab. 118 (2), 138–142. doi:10.1016/j.ymgme.2016.04.003
Cawthon, R. M., Smith, K. R., O'Brien, E., Sivatchenko, A., and Kerber, R. A. (2003). Association between telomere length in blood and mortality in people aged 60 years or older. Lancet 361 (9355), 393–395. doi:10.1016/S0140-6736(03)12384-7
Cawthon, R. M. (2002). Telomere measurement by quantitative PCR. Nucleic Acids Res. 30 (10), e47. doi:10.1093/nar/30.10.e47
Codd, V., Nelson, C. P., Albrecht, E., Mangino, M., Deelen, J., Buxton, J. L., et al. (2013). Identification of seven loci affecting mean telomere length and their association with disease. Nat. Genet. 45 (4), 422–427. doi:10.1038/ng.2528
De Meyer, T., Nawrot, T., Bekaert, S., De Buyzere, M. L., Rietzschel, E. R., and Andres, V. (2018). Telomere length as cardiovascular aging biomarker: JACC review topic of the week. J. Am. Coll. Cardiol. 72 (7), 805–813. doi:10.1016/j.jacc.2018.06.014
Dugdale, H. L., and Richardson, D. S. (2018). Heritability of telomere variation: it is all about the environment. Philos. Trans. R. Soc. Lond. B Biol. Sci. 373, 20160450. doi:10.1098/rstb.2016.0450
Epel, E. S., Blackburn, E. H., Lin, J., Dhabhar, F. S., Adler, N. E., Morrow, J. D., et al. (2004). Accelerated telomere shortening in response to life stress. Proc. Natl. Acad. Sci. U. S. A. 101 (49), 17312–17315. doi:10.1073/pnas.0407162101
Epel, E. S., Merkin, S. S., Cawthon, R., Blackburn, E. H., Adler, N. E., Pletcher, M. J., et al. (2008). The rate of leukocyte telomere shortening predicts mortality from cardiovascular disease in elderly men. Aging (Albany NY) 1 (1), 81–88. doi:10.18632/aging.100007
Frenck, R. W., Blackburn, E. H., and Shannon, K. M. (1998). The rate of telomere sequence loss in human leukocytes varies with age. Proc. Natl. Acad. Sci. U. S. A. 95 (10), 5607–5610. doi:10.1073/pnas.95.10.5607
Fuster, J. J., and Andres, V. (2006). Telomere biology and cardiovascular disease. Circ. Res. 99 (11), 1167–1180. doi:10.1161/01.RES.0000251281.00845.18
García-Calzón, S., Moleres, A., Marcos, A., Campoy, C., Moreno, L. A., Azcona-Sanjulián, M. C., et al. (2014). Telomere length as a biomarker for adiposity changes after a multidisciplinary intervention in overweight/obese adolescents: The EVASYON study. PLOS ONE 9 (2), e89828. doi:10.1371/journal.pone.0089828
Geronimus, A. T., Pearson, J. A., Linnenbringer, E., Schulz, A. J., Reyes, A. G., Epel, E. S., et al. (2015). Race-ethnicity, poverty, urban stressors, and telomere length in a detroit community-based sample. J. Health Soc. Behav. 56 (2), 199–224. doi:10.1177/0022146515582100
Hamad, R., Tuljapurkar, S., and Rehkopf, D. H. (2016). Racial and socioeconomic variation in genetic markers of telomere length: A cross-sectional study of U.S. Older adults. EBioMedicine 11, 296–301. doi:10.1016/j.ebiom.2016.08.015
Hjelmborg, J. B., Dalgard, C., Moller, S., Steenstrup, T., Kimura, M., Christensen, K., et al. (2015). The heritability of leucocyte telomere length dynamics. J. Med. Genet. 52 (5), 297–302. doi:10.1136/jmedgenet-2014-102736
Kimura, M., Hjelmborg, J. V., Gardner, J. P., Bathum, L., Brimacombe, M., Lu, X., et al. (2008). Telomere length and mortality: a study of leukocytes in elderly Danish twins. Am. J. Epidemiol. 167 (7), 799–806. doi:10.1093/aje/kwm380
Levy, D., Neuhausen, S. L., Hunt, S. C., Kimura, M., Hwang, S. J., Chen, W., et al. (2010). Genome-wide association identifies OBFC1 as a locus involved in human leukocyte telomere biology. Proc. Natl. Acad. Sci. U. S. A. 107 (20), 9293–9298. doi:10.1073/pnas.0911494107
Lin, J., Epel, E., Cheon, J., Kroenke, C., Sinclair, E., Bigos, M., et al. (2010). Analyses and comparisons of telomerase activity and telomere length in human T and B cells: insights for epidemiology of telomere maintenance. J. Immunol. Methods 352 (1-2), 71–80. doi:10.1016/j.jim.2009.09.012
Martinez, P., and Blasco, M. A. (2011). Telomeric and extra-telomeric roles for telomerase and the telomere-binding proteins. Nat. Rev. Cancer 11 (3), 161–176. doi:10.1038/nrc3025
Muezzinler, A., Zaineddin, A. K., and Brenner, H. (2013). A systematic review of leukocyte telomere length and age in adults. Ageing Res. Rev. 12 (2), 509–519. doi:10.1016/j.arr.2013.01.003
Needham, B. L., Adler, N., Gregorich, S., Rehkopf, D., Lin, J., Blackburn, E. H., et al. (2013). Socioeconomic status, health behavior, and leukocyte telomere length in the National Health and Nutrition Examination Survey, 1999-2002. Soc. Sci. Med. 85, 1–8. doi:10.1016/j.socscimed.2013.02.023
Needham, B. L., Wang, X., Carroll, J. E., Barber, S., Sanchez, B. N., Seeman, T. E., et al. (2019). Sociodemographic correlates of change in leukocyte telomere length during mid- to late-life: The Multi-Ethnic Study of Atherosclerosis. Psychoneuroendocrinology 102, 182–188. doi:10.1016/j.psyneuen.2018.12.007
Nettle, D., Gadalla, S. M., Lai, T.-P., Susser, E., Bateson, M., and Aviv, A. (2021). Telomere length measurement for longitudinal analysis: implications of assay precision. Am. J. Epidemiol. 190, 1406–1413. doi:10.1093/aje/kwab025
Okuda, K., Bardeguez, A., Gardner, J. P., Rodriguez, P., Ganesh, V., Kimura, M., et al. (2002). Telomere length in the newborn. Pediatr. Res. 52 (3), 377–381. doi:10.1203/00006450-200209000-00012
Proctor, C. J., and Kirkwood, T. B. (2002). Modelling telomere shortening and the role of oxidative stress. Mech. Ageing Dev. 123 (4), 351–363. doi:10.1016/s0047-6374(01)00380-3
Rudolph, K. L., Chang, S., Lee, H. W., Blasco, M., Gottlieb, G. J., Greider, C., et al. (1999). Longevity, stress response, and cancer in aging telomerase-deficient mice. Cell 96 (5), 701–712. doi:10.1016/s0092-8674(00)80580-2
Sanders, J. L., and Newman, A. B. (2013). Telomere length in epidemiology: a biomarker of aging, age-related disease, both, or neither? Epidemiol. Rev. 35, 112–131. doi:10.1093/epirev/mxs008
Steenstrup, T., Hjelmborg, J. V., Kark, J. D., Christensen, K., and Aviv, A. (2013). The telomere lengthening conundrum-artifact or biology? Nucleic Acids Res. 41 (13), e131. doi:10.1093/nar/gkt370
Weischer, M., Bojesen, S. E., Cawthon, R. M., Freiberg, J. J., Tybjaerg-Hansen, A., and Nordestgaard, B. G. (2012). Short telomere length, myocardial infarction, ischemic heart disease, and early death. Arterioscler. Thromb. Vasc. Biol. 32 (3), 822–829. doi:10.1161/ATVBAHA.111.237271
Keywords: genetic risk score, telomere length, telomere attrition, race/ethnicity, PCR
Citation: Castro-Diehl C, Smith JA, Zhao W, Wang X, Mukherjee B, Seeman T and Needham BL (2022) Prediction of telomere length and telomere attrition using a genetic risk score: The multi-ethnic study of atherosclerosis (MESA). Front. Aging 3:1021051. doi: 10.3389/fragi.2022.1021051
Received: 16 August 2022; Accepted: 22 September 2022;
Published: 11 October 2022.
Edited by:
Victoria Belancio, Tulane University, United StatesReviewed by:
Florencia Maria Barbé-Tuana, Pontifical Catholic University of Rio Grande do Sul, BrazilDaniel Trott, University of Texas at Arlington, United States
Copyright © 2022 Castro-Diehl, Smith, Zhao, Wang, Mukherjee, Seeman and Needham. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Cecilia Castro-Diehl, ocastro-diehl@bwh.harvard.edu