Abstract
Are shorter telomeres causal risk factors for facial aging on a large population level? To examine if longer, genetically predicted telomeres were causally associated with less facial aging using Mendelian randomization analysis. Two-sample Mendelian randomization methods were applied to the summary statistics of a genome-wide association study (GWAS) for self-reported facial aging from 417, 772 participants of the UK Biobank data. Twenty single-nucleotide polymorphisms (SNPs) that were of genome-wide significance were selected as instrumental variables for leukocyte telomere length. The main analyses were performed primarily using the random-effects inverse-variance weighted method and were complemented with the MR-Egger regression, weighted median, and weighted mode approaches. The intercept of MR-Egger regression was used to assess horizontal pleiotropy. Longer genetically predicted telomeres were associated with a lower likelihood of facial aging (β = − 0.02, 95% confidence interval: − 0.04, − 0.002). Comparable results were obtained using MR-Egger regression, weighted median, and weighted mode approaches. The intercept of MR-Egger regression was close to zero (0.002) that was not suggestive of horizontal pleiotropy. Our findings provided evidence to support a potential causal relationship between longer genetically predicted telomeres and less facial aging.
Avoid common mistakes on your manuscript.
Introduction
Telomeres are comprised of hexanucleotide DNA repeats and a complex structure of surrounding proteins at the end of chromosomes, protecting genetic information by maintaining the stability of chromosomes during cellular divisions [1]. Each time a cell divides, a small amount of telomeric DNA is lost because of the inability of the polymerase to fully elongate the ends of DNA. Consequently, telomeres are shorten with each cell division and are therefore recognized as a potential biological marker for cellular aging [2]. Telomere attrition over time results in critically short telomere lengths and leads to cellular senescence and apoptosis in normal cells. Skin cells, whose turnover cycle could be driven in part by telomeres replenishment, have a proliferative capacity that varies from a few days in newborns to a few weeks in adults [3, 4]. And facial skin cells are of particular interest. Maintaining a fast renewal of skin cells could slow down facial aging and make people look younger.
Different molecular mechanisms have been suggested to explain facial aging, where telomeres as a marker for biological aging have recently been gaining momentum [5]. For example, an in vitro experiment showed that skin cells with telomeres lengthened by a procedure that delivered a modified mRNA encoding TERT to cells were able to divide more times than untreated cells [6]. However, studies like this one need to be further investigated in large population-based analyses.
Therefore, in this study, we aimed to examine the potential causal association between leukocyte telomere length (TL), which was correlated with TL in skin cells (correlation coefficient r = 0.83) [7], and facial aging using a Mendelian randomization (MR) approach [8] in data collected from 417,772 participants of the UK Biobank [9].
Materials and methods
Instrumental variable selection
The European Network for Genetic and Genomic Epidemiology (ENGAGE) conducted a genome-wide association study (GWAS) for leukocyte TL in 78,592 individuals of European ancestry [10]. Mean leukocyte TL was measured in a mixed population of leukocytes, and measurements were conducted using an established quantitative polymerase chain reaction technique which expressed TL as a ratio of the telomere repeat number (T) to a single-copy gene (S) [11]. Leukocyte TL measurements were standardized either by using a calibrator sample or by quantifying against a standard curve. In total, 20 single-nucleotide polymorphisms (SNPs) at 17 genomic loci were independently associated with leukocyte TL at a level of genome-wide statistical significance (P < 5 × 10−8). In our study, we use these 20 SNPs as instrumental variables, and included proxy SNPs through LDlink if SNPs found to be palindromic (Table 1).
Facial aging GWAS
In the UK Biobank, the facial aging domain was measured by the question (field code 1757)—Do people say that you look younger than you are, older than you are, about your age, do not know, and prefer not to answer? Analyses of this variable were performed in 417,772 British participants using a mixed linear model-based tool (fastGWA) that has treated facial aging as an ordinal categorical variable while controls for population stratification by principal components and for relatedness by a sparse genetic relationship matrix as well as age and sex [9]. These data were analyzed and made publicly available by the Complex Trait Genomics lab (http://fastgwa.info).
Statistical analysis
The data on the instrumental variables and the facial aging GWAS were harmonized by the respective chromosomes and positions (human genome build 37). The random-effects inverse-variance weighted (IVW) method was used as the primary estimator for the MR analysis [12]. This method has a higher statistical power with the assumption that all SNPs are valid instrumental variables. The weighted median approach [13], MR-Egger regression [14], and weighted mode method [15] were used as complementary analysis. The weighted median approach yields consistent estimates when at least 50% of the weights in the analysis are from valid instrumental variables. The MR-Egger regression can adjust for directional pleiotropy but is of low power. The intercept of MR-Egger regression is used as a test for horizontal pleiotropy. The weighted mode approach is consistent when the largest number of similar individual-instrument causal effect estimates comes from valid instruments, even if most instruments are invalid. To examine if there was a reverse causation between TL and facial aging, we conducted an additional MR analysis as bi-directional MR analysis using 98 SNPs of genome-wide significance for facial aging as instrumental variables in UKB. Then, we harmonized their respective effect sizes and standard errors for these SNPs in ENGAGE TL GWAS. All statistical analyses were conducted in R 3.6 and TwoSampleMR package [16].
The present study only used publicly available summary-level statistics. No individual-level data was analyzed. Ethical approval is therefore not required.
Results
Table 1 describes 20 genetic variants as instrumental variables and their associations with leukocyte TL. Six palindromic SNPs were replaced by their corresponding proxy SNPs (r2 = 1). The scatter plot for the effects of these SNPs on leukocyte TL and facial aging is shown in Fig. 1.
The MR analysis shows that longer genetically predicted leukocyte TL was associated with a lower likelihood of facial aging (β = − 0.02, 95% confidence interval [CI]: − 0.04, − 0.002) using the IVW method (Table 2). Similar results were obtained through MR-Egger regression (β = − 0.06, 95% CI: − 0.12, − 0.002), weighted median (β = − 0.03, 95% CI: − 0.05, − 0.007), and weighted mode approaches (β = − 0.03, 95% CI: − 0.05, − 0.007). These estimates were also plotted in Fig. 1.
We additionally performed a leave-one-out analysis, which yielded comparable results and did not find noticeable effects of any single SNP that could dominate the results (Fig. 2). The funnel plot did not imply that there were heterogeneous SNPs (Fig. 3). We did not find evidence for directional pleiotropy from the MR-Egger regression (intercept: 0.002, P = 0.16).
Additional bi-directional MR analysis was performed by using 98 SNPs that were of genome-wide significance for facial aging as instrumental variables. We did not find a significant effect of genetically predicted facial aging on leukocyte TL (β = − 0.08, 95% CI: − 0.22, 0.05).
Discussion
In the present study, for the first time, we examined the association of genetically predicted leukocyte TL and facial aging in a large population-based cohort, the UK Biobank, using the MR method. By leveraging several MR estimation approaches, we found that longer genetically predicted leukocyte TL was associated with a lower likelihood of facial aging. Our study, corroborating previous experimental studies of skin cells, provides further evidence to support a causal role of leukocyte TL in facial aging.
To the best of our knowledge, no studies have been published to investigate the role of TL in facial aging using population-based cohorts until now. A previous study examined TL in cells of the sun-protected and sun-exposed skins [17]. Further studies employed ultraviolet-exposed skin cells to study the roles of TL in photoaging [18,19,20,21]. These studies, taken together, were suggestive of the functions of TL in skin aging at the cellular level [5, 22]. Further evidence at the population level, however, is lacking. The limited number of studies on this topic at the population level is partly due to the difficulties in facial aging measurements. For example, although the three-dimensional human facial morphology assessment instruments [23], among others [24], can offer objective and comprehensive features for facial aging, they have not been widely available for aging researchers because of infeasibility and high cost.
To overcome the challenges in facial aging measurement and to make the best use of questionnaire-based instruments, a subjective assessment of facial aging could be valuable in this regard measurement [25,26,27]. In this study, we examined the role of TL by taking advantage of questionnaire-based measurement for facial aging in the UK Biobank. An advantage of this type of data is that the large sample size in the UK Biobank could outweigh the concerns in measurement error for questionnaire-based facial aging assessment. Indeed, we obtained consistent and significant estimates by using various MR approaches in this study. Further studies are warranted to explore the biological mechanisms of TL in facial aging.
Despite the advantages of the large sample size, our study is prone to several limitations. First, the MR assumptions, particularly the no-pleiotropy assumption, must be satisfied in order to yield a valid estimate. The intercept of MR-Egger regression analysis approaches zero suggesting no strong evidence for directional pleiotropy. Second, the MR methods used in this study also assume a linear relationship between TL and facial aging. We cannot examine further if there was a non-linear relationship. Third, as facial aging GWAS was analyzed using mixed linear models, the effect sizes cannot be interpreted using a more intuitive way. The results could rather be used for testing purposes. Fourth, as alluded above, the facial aging was measured using a single question in the UK Biobank, which could lead to measurement error. Future studies are encouraged in order to develop a more detailed questionnaire-based instrument to measure diverse dimensions of facial aging. Fifth, in this study, facial aging was assessed only using qualitative data and not objectively assessed using skin turgor [28] or facial wrinkles per area of skin [29]. Lastly, TL was measured in leukocyte rather than skin cells. However, a previous study demonstrated that the correlation of TL in leukocyte and skin cells was high (correlation coefficient r = 0.83) [7], and a recent study found that TL varies by tissue type but is generally correlated among tissue types [30].
In summary, our study provided novel evidence to support a causal role of genetically predicted leukocyte TL in facial aging. Further studies are warranted to explore the biological mechanisms of TL in facial aging.
Data availability
The data are publicly available.
References
Blackburn EH, Epel ES, Lin J. Human telomere biology: a contributory and interactive factor in aging, disease risks, and protection. Science. 2015;350(6265):1193–8. https://doi.org/10.1126/science.aab3389.
Jylhava J, Pedersen NL, Hagg S. Biological age predictors. EBioMedicine. 2017;21:29–36. https://doi.org/10.1016/j.ebiom.2017.03.046.
Iizuka H. Epidermal turnover time. J Dermatol Sci. 1994;8(3):215–7. https://doi.org/10.1016/0923-1811(94)90057-4.
Weinstein GD, McCullough JL, Ross P. Cell proliferation in normal epidermis. J Invest Dermatol. 1984;82(6):623–8. https://doi.org/10.1111/1523-1747.ep12261462.
Zhang S, Duan E. Fighting against skin aging: the way from bench to bedside. Cell Transplant. 2018;27(5):729–38. https://doi.org/10.1177/0963689717725755.
Ramunas J, Yakubov E, Brady JJ, Corbel SY, Holbrook C, Brandt M, et al. Transient delivery of modified mRNA encoding TERT rapidly extends telomeres in human cells. FASEB J. 2015;29(5):1930–9. https://doi.org/10.1096/fj.14-259531.
Daniali L, Benetos A, Susser E, Kark JD, Labat C, Kimura M, et al. Telomeres shorten at equivalent rates in somatic tissues of adults. Nat Commun. 2013;4:1597. https://doi.org/10.1038/ncomms2602.
Smith GD, Ebrahim S. ‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32(1):1–22. https://doi.org/10.1093/ije/dyg070.
Jiang L, Zheng Z, Qi T, Kemper KE, Wray NR, Visscher PM, et al. A resource-efficient tool for mixed model association analysis of large-scale data. Nat Genet. 2019;51(12):1749–55. https://doi.org/10.1038/s41588-019-0530-8.
Li C, Stoma S, Lotta LA, Warner S, Albrecht E, Allione A, et al. Genome-wide association analysis in humans links nucleotide metabolism to leukocyte telomere length. Am J Hum Genet. 2020;106(3):389–404. https://doi.org/10.1016/j.ajhg.2020.02.006.
Cawthon RM. Telomere length measurement by a novel monochrome multiplex quantitative PCR method. Nucleic Acids Res. 2009;37(3):e21. https://doi.org/10.1093/nar/gkn1027.
Burgess S, Scott RA, Timpson NJ, Davey Smith G, Thompson SG, Consortium E-I. Using published data in Mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur J Epidemiol. 2015;30(7):543–52. https://doi.org/10.1007/s10654-015-0011-z.
Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40(4):304–14. https://doi.org/10.1002/gepi.21965.
Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44(2):512–25. https://doi.org/10.1093/ije/dyv080.
Hartwig FP, Davey Smith G, Bowden J. Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol. 2017;46(6):1985–98. https://doi.org/10.1093/ije/dyx102.
Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7. https://doi.org/10.7554/eLife.34408.
Sugimoto M, Yamashita R, Ueda M. Telomere length of the skin in association with chronological aging and photoaging. J Dermatol Sci. 2006;43(1):43–7. https://doi.org/10.1016/j.jdermsci.2006.02.004.
Ma HM, Liu W, Zhang P, Yuan XY. Human skin fibroblast telomeres are shortened after ultraviolet irradiation. J Int Med Res. 2012;40(5):1871–7. https://doi.org/10.1177/030006051204000526.
Kosmadaki MG, Gilchrest BA. The role of telomeres in skin aging/photoaging. Micron. 2004;35(3):155–9. https://doi.org/10.1016/j.micron.2003.11.002.
Rochette PJ, Brash DE. Human telomeres are hypersensitive to UV-induced DNA damage and refractory to repair. PLoS Genet. 2010;6(4):e1000926. https://doi.org/10.1371/journal.pgen.1000926.
Holck DE, Ng JD. Facial skin rejuvenation. Curr Opin Ophthalmol. 2003;14(5):246–52. https://doi.org/10.1097/00055735-200310000-00004.
Buckingham EM, Klingelhutz AJ. The role of telomeres in the ageing of human skin. Exp Dermatol. 2011;20(4):297–302. https://doi.org/10.1111/j.1600-0625.2010.01242.x.
Chen W, Qian W, Wu G, Chen W, Xian B, Chen X, et al. Three-dimensional human facial morphologies as robust aging markers. Cell Res. 2015;25(5):574–87. https://doi.org/10.1038/cr.2015.36.
Windhager S, Mitteroecker P, Rupic I, Lauc T, Polasek O, Schaefer K. Facial aging trajectories: a common shape pattern in male and female faces is disrupted after menopause. Am J Phys Anthropol. 2019;169(4):678–88. https://doi.org/10.1002/ajpa.23878.
Dykiert D, Bates TC, Gow AJ, Penke L, Starr JM, Deary IJ. Predicting mortality from human faces. Psychosom Med. 2012;74(6):560–6. https://doi.org/10.1097/PSY.0b013e318259c33f.
Christensen K, Thinggaard M, McGue M, Rexbye H, Hjelmborg JV, Aviv A, et al. Perceived age as clinically useful biomarker of ageing: cohort study. BMJ. 2009;339:b5262. https://doi.org/10.1136/bmj.b5262.
Rexbye H, Petersen I, Johansens M, Klitkou L, Jeune B, Christensen K. Influence of environmental factors on facial ageing. Age Ageing. 2006;35(2):110–5. https://doi.org/10.1093/ageing/afj031.
Fortes MB, Owen JA, Raymond-Barker P, Bishop C, Elghenzai S, Oliver SJ, et al. Is this elderly patient dehydrated? Diagnostic accuracy of hydration assessment using physical signs, urine, and saliva markers. J Am Med Dir Assoc. 2015;16(3):221–8. https://doi.org/10.1016/j.jamda.2014.09.012.
Hamer MA, Pardo LM, Jacobs LC, Ikram MA, Laven JS, Kayser M, et al. Lifestyle and physiological factors associated with facial wrinkling in men and women. J Invest Dermatol. 2017;137(8):1692–9. https://doi.org/10.1016/j.jid.2017.04.002.
Demanelis K, Jasmine F, Chen LS, Chernoff M, Tong L, Shinkle J et al. Determinants of telomere length across human tissues. BioArxiv. 2019. doi: https://doi.org/10.1101/793406.
Acknowledgments
The authors thank ENGAGE telomere consortium for sharing GWAS summary statistics.
Funding
Open Access funding enabled and organized by Projekt DEAL. This is study was supported by a Junior Group Leader Startup Grants at DZNE. SH is supported by grants from the Swedish Research Council (2019-01272, 2015-03255), the Strategic Research Program in Epidemiology at Karolinska Institutet, and King Gustaf V:s and Queen Victoria’s Freemason Foundation.
Author information
Authors and Affiliations
Contributions
YZ and SH contributed to study design, data analysis, manuscript writing, and results interpretation.
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflicts of interest.
Ethical approval
The present study only used summary statistics rather than individual-level data. IRB is exempted according to the rules at Karolinska Institutet, Sweden.
Consent to participate
Not applicable.
Consent for publication
Not applicable.
Code availability
The codes are from TwoSampleMR package.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
This article is published under an open access license. Please check the 'Copyright Information' section either on this page or in the PDF for details of this license and what re-use is permitted. If your intended use exceeds what is permitted by the license or if you are unable to locate the licence and re-use information, please contact the Rights and Permissions team.
About this article
Cite this article
Zhan, Y., Hägg, S. Association between genetically predicted telomere length and facial skin aging in the UK Biobank: a Mendelian randomization study. GeroScience 43, 1519–1525 (2021). https://doi.org/10.1007/s11357-020-00283-0
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11357-020-00283-0