Low Blood Levels of Selenium, Selenoprotein P and GPx3 are Associated with Accelerated Biological Aging: Results from the Berlin Aging Study II (BASE-II) ========================================================================================================================================================= * Valentin Max Vetter * Kamil Demircan * Jan Homann * Thilo Samson Chillon * Michael Mülleder * Orr Shomroni * Elisabeth Steinhagen-Thiessen * Markus Ralser * Christina M. Lill * Lars Bertram * Lutz Schomburg * Ilja Demuth ## Abstract **Introduction** Biological age reflects inter-individual differences in biological function and capacity beyond chronological age. Biological age can be estimated by DNA methylation age (DNAmA) and its deviation from chronological age, DNAmA acceleration (DNAmAA). Low levels of serum selenium, selenoprotein P (SELENOP), and the selenocysteine-containing glutathione peroxidase 3 (GPx3) are associated with adverse health outcomes and selenium supplementation is discussed as an anti-aging intervention. **Methods** In this study we analyzed 1,568 older participants from the Berlin Aging Study II (mean age +/− SD: 68.8 +/− 3.7 years, 51% women). DNAmA was estimated from genome-wide DNA methylation data using the Horvath, GrimAge, and DunedinPACE algorithms. Serum selenium levels were measured by total reflection X-ray fluorescence (TXRF) spectroscopy. SELENOP was measured by ELISA and GPx3 was derived from a larger set of mass spectrometry proteomics data. **Results** Participants with deficient serum selenium levels (<90μg/L) had a higher rate of biological aging (DunedinPACE, p=0.01, n=865). This association remained statistically significant after adjustment for age, sex, BMI, smoking, and genetic ancestry (β=-0.02, SE=0.01, 95%CI: - 0.034 to −0.004, n=757). Compared to the highest quartile, participants in the lowest quartile of SELENOP levels showed an accelerated biological aging rate (DunedinPACE, β=-0.03, SE=0.01, 95%CI: −0.051 to −0.008, n=740, fully adjusted model). Similarly, after adjustment for covariates, accelerated biological age was found in participants within the lowest GPx3 quartile compared to participants in the fourth quartile (DunedinPACE, p=<0.001 and GrimAge, p<0.001). **Conclusion** Our study suggests that low levels of selenium biomarkers are associated with accelerated biological aging measured as DNAmA. This effect was not substantially changed after adjustment for known covariates. ## Introduction Over the past few decades the average lifespan increased faster than the healthspan, which resulted in a growing burden of late-life diseases (1, 2). However, there are substantial inter-individual differences in the pace of aging and some people remain disease-free and in good health until old age. Although some lifestyle habits were shown to impact some areas of aging (3) they do not fully explain the observed between-person differences. Hence, the identification of additional factors promoting a “healthy aging” trajectory would potentially benefit many. One recent, and often considered very useful (4, 5), biomarker to quantify difference in biological aging is DNA methylation age (DNAmA, epigenetic age) estimated from DNA methylation data by so called epigenetic clocks (6). The difference between epigenetic and chronological age, DNAmA acceleration (DNAmAA), was shown in previous studies to be associated with mortality as well as with numerous age-related diseases and phenotypes (6). Due to the ability of these clocks to assess the biological age of an individual they are increasingly used as endpoints for the evaluation of anti-aging interventions (reviewed in (7–11)). In related research, one factor discussed in the context of healthy aging is selenium, an essential trace element that regulates thyroid hormone metabolism, antioxidative, and redox processes through incorporation into selenoproteins (12), which needs to be present in sufficient levels to exert its physiological function. Selenium status and its cellular availability (determined by dietary intake and its absorption) directly regulate expression and activity of selenoproteins (13). Selenium status can be determined by assessment of serum selenium biomarkers such as total selenium, the selenium transporter selenoprotein P (SELENOP), or glutathione peroxidase 3 (GPx3), a selenocysteine-containing extracellular antioxidant enzyme (14). Emerging clinical evidence indicates an association between a deficient status of these selenium biomarkers and outcomes of mortality and chronic diseases of older age, i.e., measures of an increased pace of aging (15–19). The direct association between selenium biomarkers and epigenetic aging, however, remains understudied. The relationship between blood selenium levels and DNAmA was analyzed in two previous studies, which revealed promising links between selenium levels and biological aging (20, 21). However, both studies were of small sample size (n=93 (20) and n=276 (21)), and analyzed only total selenium, but none of the other selenium-related variables. In this study, we analyzed the association between biological aging as determined with three DNAmA estimators (i.e. Horvath (22), GrimAge (23), and DunedinPACE (24)) total serum selenium, and two complementary protein biomarkers, SELENOP, GPx3, in a large cohort of 1,568 older participants from the Berlin Aging Study II (BASE-II). ## Methods ### Study population The Berlin Aging Study II (BASE-II) is a multi-disciplinary study that aims to identify factors promoting healthy aging trajectories. BASE-II comprises 1,671 participants (60 – 85 years), who were recruited in the metropolitan area of Berlin, Germany from 2009 to 2014. Older participants who reported the intake of selenium supplements were excluded (n=103) resulting in a final dataset of up to 1,568 participants. A detailed description of the study’s cohort profile was published before and can be found in ref. (25). ### DNA methylation age (DNAmA), DNAmA acceleration (DNAmAA) and pace of aging The participants’ epigenetic age was estimated from DNA methylation data utilizing three “first-generation” clocks (Horvath (22), Hannum (26), 7-CpG (27)), two “second-generation” clocks (PhenoAge (28), GrimAge (23)) and DunedinPACE (24), and a third-generation clock. DNAm data for the 7-CpG clock was measured by Single Nucleotide Prime Extension (SNuPE) in 1471 samples. The laboratory protocol is described in detail elsewhere (27). All other clocks were derived from DNAm data measured by the “Infinium MethylationEPIC” array (Illumina, Inc., USA), data of which was available for n=1,030 participants. Additional information on laboratory procedures and data processing can be found in reference (29). The DunedinPACE clock was calculated based using the methods described in the original publication (24). DNA methylation age acceleration (DNAmAA) was calculated for all clocks (except for DunedinPACE) as residuals of a linear regression of DNAmA on chronological age accounting for leukocyte cell distribution (neutrophils, monocytes, lymphocytes, and eosinophils in G/l). ### Assessment of selenium, SELENOP and GPx3 The selenium biomarkers *SELENOP* and *GPx3* were measured at the Institute of Experimental Endocrinology, Charité University Hospital, Berlin, from serum samples that were taken during the same blood draw that produced the sample for the DNA isolation. Selenium status was measured by total reflection X-ray fluorescence (TXRF) spectroscopy as described in detail before (30, 31). SELENOP was measured from serum samples with a sandwich ELISA method using monoclonal human antibodies. Additional information on how selenium and SELENOP were measured in BASE-II were previously described in detail (32). GPx3 intensities were derived from proteomics data at baseline was measured by Liquid chromatography-mass spectrometry (LC-MS), using data independent acquisition (Supplementary Methods). ### Covariates Chronological age, sex, and smoking behavior (in packyears) were assessed during a one-on-one interview with trained study personnel. Body weight and height were measured with the measuring station (763, SECA, Germany). Genetic ancestry was controlled for by using the first four principal components from a principal component analysis on genome-wide SNP genotyping data. A detailed description of the procedures applied to calculate this co-variate were published previously (33). Leukocyte cell distribution was measured from samples obtained at the same blood-draw that was used for DNA isolation in an accredited clinical biochemistry laboratory (MVZ Labor 28 GmbH, Berlin, Germany) by automated standard methods (flow cytometry). ### Statistical analysis Statistical analyses were conducted in R (version 4.3.1) (34). Descriptive statistics were calculated with R’s tableone package (35). Linear regression analyses were calculated with R’s *lm* and *lm.beta* function (lm.beta package (36)) to examine the relationship between DNAmA estimators (dependent variable) and selenium biomarkers (independent variables). Figures were drawn with the ggplot2 package. Selenium biomarker outliers were defined as <25th percentile −1.5*SD and >the 75th percentile +1.5*SD. This resulted in the exclusion of n=46 (serum selenium), n=111 (SELENOP), and n=151 (GPx3) observations. Participants with missing values in phenotypic data were excluded from the respective analyses (available case analyses). Statistical significance was defined at α=5%. Due to the exploratory nature of this study, no adjustment was made for multiple testing. ## Results ### Study population and characteristics The sample analyzed in this study comprised 1,568 BASE-II participants with a mean age of 69 years (SD=3.4, range 60-85 years). Sex was almost equally distributed (51% women). 656 (48.3%) participants had a selenium level of 90μg/L or lower, i.e. displayingselenium-deficiency. Study population characteristics are shown in Table 1. View this table: [Table 1:](http://medrxiv.org/content/early/2024/04/05/2024.04.04.24305314/T1) Table 1: Baseline characteristics of analyzed sample. ### Association of selenium with biological age acceleration Visual inspections of the distribution of serum selenium levels vs. DNAmAA estimated by the GrimaAge and DunedinPACE clocks suggested possible non-linear associations (Figure 1). Therefore, in the next step we examined epigenetic age stratified by selenium deficiency (<90μg/L) and found a statistically significantly higher pace of aging measured by DunedinPACE in selenium-deficient participants (SMD=0.2, p=0.01, t-test, Figure 1 F). DunedinPACE remained higher in selenium-deficient participants after the adjustment for known covariates (chronological age, sex, BMI, smoking, genetic ancestry) in linear regression analyses (β=-0.02, SE=0.007, p=0.012, n=757, Table 2). Sex-stratified analyses showed that these results were mainly driven by the subgroup of men (Supplementary Table 2). ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/04/05/2024.04.04.24305314/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2024/04/05/2024.04.04.24305314/F1) Figure 1: A-C: Scatterplots of serum selenium levels at baseline and biological age estimators. The x-axis is log-scaled. All available participants of the older age group are included. D-F: Boxplots of biological age estimators stratified by selenium status (deficient vs. sufficient. cut-off 90 μg/L) at baseline. Statistical significance of difference between group means was assessed by t-test. Note: DNAmAA = DNA methylation age acceleration; biolog. = biological years; chron. = chronological years. Boxplots show the median, hinges (corresponding to 25th and 75th percentile) and whiskers (1.5*inter-quartile range (IQR)). View this table: [Table 2:](http://medrxiv.org/content/early/2024/04/05/2024.04.04.24305314/T2) Table 2: Linear regression analysis of epigenetic age estimators on selenium status (deficient vs. sufficient) at baseline. Model 1 is unadjusted. Model 2 is adjusted for chronological age, sex, BMI, smoking (packyears), and the first four genetic principal components (PC1 to PC4). ### Association of SELENOP with biological age Low SELENOP levels were similarly associated with accelerated GrimAge DNAmAA and DunedinPACE (Figure 2). Assuming a non-linear relationship based on the scatterplot (Figure 2 A-C), SELENOP values were categorized into quartiles. A statistically significantly higher pace of biological aging (DunedinPACE) was measured in participants of the lowest quartile compared to participants in the highest quartile (SMD = 0.21, p=0.032, Figure 2 F). This association remained statistically significant after covariate adjustment (β=-0.03, SE=0.011, p=0.007, Table 3). After covariate adjustment the findings were only statistically significant in the subgroup of women (Supplementary Table 4 and 5). ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/04/05/2024.04.04.24305314/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2024/04/05/2024.04.04.24305314/F2) Figure 2: A-C: Scatterplots of SELENOP levels and biological age estimators. The x-axis is log-scaled. D-F: Boxplots of biological age estimators stratified by quartiles of SELENOP. Statistical significance of difference between group means was assessed by t-test. Note: DNAmAA = DNA methylation age acceleration; biolog. = biological years; chron. = chronological years. Boxplots show the median, hinges (corresponding to 25th and 75th percentile) and whiskers (1.5*inter-quartile range (IQR)). View this table: [Table 3:](http://medrxiv.org/content/early/2024/04/05/2024.04.04.24305314/T3) Table 3: Linear regression analysis of epigenetic age estimators on quartiles of SELENOP at baseline. Model 1 is unadjusted. Model 2 is adjusted for chronological age, sex, BMI, smoking (packyears), and the first four genetic principal components (PC1 to PC4). The first quartile is used as reference. ### Association of GPx3 with biological age Finally, compared to the lowest quartile, participants in the fourth quartile had a higher age acceleration estimated from the GrimAge clock (SMD = 0.41, p<0.001) and DunedinPACE clock (SMD = 0.54, p<0.001, Figure 3). These associations persisted to be statistically significant after covariate adjustment (p<=0.002, Table 4). Sex-specific analyses are presented in Supplementary Table 7 and 8. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/04/05/2024.04.04.24305314/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2024/04/05/2024.04.04.24305314/F3) Figure 3: A-C: Scatterplots of GPx3 intensity (proteomics data) and biological age estimators. The x-axis is log-scaled. D-F: Boxplots of biological age estimators stratified by quartiles of GPx3 intensity (proteomics data). Statistical significance of difference between group means was assessed by t-test. Note: DNAmAA = DNA methylation age acceleration; biolog. = biological years; chron. = chronological years. Boxplots show the median, hinges (corresponding to 25th and 75th percentile) and whiskers (1.5*inter-quartile range (IQR)). View this table: [Table 4:](http://medrxiv.org/content/early/2024/04/05/2024.04.04.24305314/T4) Table 4: Linear regression analysis of epigenetic age estimators on quartiles of GPx3 intensity (proteomics data) at baseline. Model 1 is unadjusted. Model 2 is adjusted for chronological age, sex, BMI, smoking (packyears), and the first four genetic principal components (PC1 to PC4). The first quartile is used as reference group. ## Discussion In this study, we analyzed the association between serum biomarkers, namely total serum selenium, SELENOP, GPx3 and biological age measured by epigenetic clocks in 865 participants of BASE-II. Lower values in all three selenium biomarkers were associated with an increased pace of aging measured with the DunedinPACE clock. This association was independent from covariates included in the linear regression analysis. The stepwise increase in the strength of association from the first- to the third-generation clock is likely due to the methodological approach used for their development with the later clocks potentially including CpG sites that better reflect the underlying aging process (37). Furthermore, newer DNAmA clocks are more sensitive towards exogenous effects (20). This trend was also observed with clinical phenotypes (7, 24, 38). Our findings are in line with results published by Cheng and colleagues (20) who found a statistically significant association between serum selenium levels and DunedinPACE in a sample of 93 participants. In contrast, we could not validate the association between selenium levels and PhenoAge DNAmAA reported by the same authors (Supplementary Table 2). Furthermore, regarding our non-significant associations of selenium levels with earlier DNAmA clocks, our results are consistent with findings from a cohort of 276 older Chinese men and women that showed no statistically significant association between selenium and the first- and second-generation clocks (21). Boyer and colleagues did not find a consistent association between urine selenium and any epigenetic age estimation in linear regression models in an American Indian population of 2,301 participants (Strong Heart Study, SHS) (37). However, when analyzing the effect of the total metal mixture using Bayesian Kernel Machine Regression, a non-linear association with PhenoAge, GrimAge and DunedinPACE was found (37). The difference in these results compared to our study may arise from the choice of matrix used to determine selenium status. Specifically, the selenium status in the SHS was replete (37), whereas BASE-II participants showed a borderline selenium-deficiency. The lack of effects observed in SHS may also be attributable to the saturated expression of selenoproteins, and thus likely reflecting a threshold effect as commonly observed for selenium and health outcomes. Our analyses do not allow to draw any conclusions on cause-effect relationships between selenium levels and accelerated biological aging. However, our results corroborate recent findings on aging phenotypes assessed by other clinical and phenotypical outcomes that show an association between selenium biomarkers and mortality (15, 19, 39), cardiovascular outcomes (16, 17, 39), and several cancer entities (14, 40–42). Moreover, an effect of selenium levels on epigenetic age appears biologically plausible since it is well established that selenium levels in rodents, cell-lines (human, mouse) and human tissue impact the methylome (reviewed in (43)). Strengths of this study include the large sample size with an equal representation of men and women and a high degree of data completeness for the study parameters and potential confounders. Furthermore, this is the first study that examines additional selenium biomarkers, i.e., SELENOP and GPx3, in addition to selenium levels. The differing half-lives of these biomarkers ensure low risk of misclassification due to e.g., prior supplementation or enriched dietary selenium intake. However, our study also has a few limitations. Firstly, as outlined above, the cross-sectional analyses presented in this study do not allow to draw any conclusions about causality or direction of effect. Secondly, the analyzed sample is of above-average health. Therefore, it is possible that the effects reported here are underestimated since BASE-II participants are generally healthier. In addition, this limits the generalizability of our results to the background population of mainly European elderly subjects with borderline sufficient selenium status residing in a metropolitan area. Thirdly, due to the explorative nature of this analysis, no adjustment for multiple testing was made. Fourthly, no detailed information about dosage and frequency of potential additional selenium intake by nutritional supplements was available. Therefore, all participants who reported selenium supplementation were excluded and evaluating effects of selenium supplementation on biological age was precluded. ## Conclusion Our study shows that selenium biomarkers are moderately associated with accelerated biological aging in a European-based dataset. These results support previous findings that suggest health benefits of sufficient selenium levels. Randomized controlled trials are needed to evaluate the potential therapeutic benefit of selenium supplementation and its effect on biological age acceleration. ## Supporting information Supplementary Material [[supplements/305314_file02.pdf]](pending:yes) ## Funding This work was supported by grants of the Deutsche Forschungsgemeinschaft (grant number 460683900 to ID and LB, and CRC/TR 296 “LocoTact” to LS), the ERC (as part of the “Lifebrain” project to LB), and the Cure Alzheimer’s Fund (as part of the “CIR-CUITS” consortium to LB), the German Federal Ministry of Education and Research as part of the National Research Initiative ‘Mass Spectrometry in Systems Medicine’ (MSCoreSys), under grant agreement number 01EP2201 (to MR) and 16LW0239K (to MM), the Berlin University Alliance (BUA Link Lab, 501_Massenspektrometrie, 501_Linklab), as well as the German Cancer Consortium (DKTK) under grant BE01 1020000483 (to MR). C.M.L. was supported by the Heisenberg program of the German Research Foundation (DFG; LI 2654/4-1). This article uses data from the Berlin Aging Study II (BASE-II). BASE-II was supported by the German Federal Ministry of Education and Research under grant numbers #01UW0808; #16SV5536K, #16SV5537, #16SV5538, #16SV5837, #01GL1716A, and #01GL1716B. ## Conflict of interest LS holds shares of selenOmed GmbH, a company involved in Se status assessment; no other relationships or activities that could appear to have influenced the submitted work are indicated. Other authors: none. ## Data Availability Due to concerns for participant privacy, data are available only upon reasonable request. Please contact Ludmila MuCller, scientific coordinator, at lmueller@mpib-berlin.mpg.de, for additional information. ## Ethics All participants gave written informed consent. The Ethics Committee of the Charité – Universitätsmedizin Berlin approved the study (approval numbers EA2/029/09 and EA2/144/16). The study was conducted in accordance with the Declaration of Helsinki and was registered in the German Clinical Trials Registry as DRKS00009277. ## Author contributions Conceptualization: V.M.V., K.D., T.S.C.; L.S., I.D.; Data curation: V.M.V., J.H., M.M., O.S., L.B., I.D.; Formal analysis: V.M.V.; Funding acquisition: E.S.-T.; M.R., L.B., L.S., I.D.; Resources: E.S.-T.; M.R., L.B., L.S., I.D.; Supervision: C.L., M.R., L.B., L.S., I.D.; Visualization: V.M.V.; Writing: V.M.V, K.D., I.D.; Reviewed and approved the manuscript: all authors. ## Acknowledgments ## Footnotes * a Valentin Max Vetter and Kamil Demircan are joint first authors. * b Lutz Schomburg and Ilja Demuth are joint last authors * Received April 4, 2024. * Revision received April 4, 2024. * Accepted April 5, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Partridge L, Deelen J, Slagboom PE. Facing up to the global challenges of aging. Nature. 2018;561(7721):45–56. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-018-0457-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30185958&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F05%2F2024.04.04.24305314.atom) 2. 2.Beard JR, Officer A, De Carvalho IA, Sadana R, Pot AM, Michel J-P, et al. The World report on aging and health: a policy framework for healthy aging. The lancet. 2016;387(10033):2145–54. 3. 3.Li Y, Schoufour J, Wang DD, Dhana K, Pan A, Liu X, et al. Healthy lifestyle and life expectancy free of cancer, cardiovascular disease, and type 2 diabetes: prospective cohort study. BMJ (Clinical research ed). 2020;368:l6669. 4. 4.Jylhava J, Pedersen NL, Hagg S. Biological Age Predictors. EBioMedicine. 2017;21:29–36. 5. 5.Li X, Ploner A, Wang Y, Magnusson PK, Reynolds C, Finkel D, et al. Longitudinal trajectories, correlations and mortality associations of nine biological ages across 20-years follow-up. Elife. 2020;9:e51507. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.51507&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F05%2F2024.04.04.24305314.atom) 6. 6.Horvath S, Raj K. DNA methylation-based biomarkers and the epigenetic clock theory of aging. Nature Reviews Genetics. 2018:1. 7. 7.Waziry R, Ryan CP, Corcoran DL, Huffman KM, Kobor MS, Kothari M, et al. Effect of long-term caloric restriction on DNA methylation measures of biological aging in healthy adults from the CALERIE trial. Nature Aging. 2023;3(3):248–57. 8. 8.Vetter VM, Sommerer Y, Kalies CH, Spira D, Bertram L, Demuth I. Vitamin D supplementation is associated with slower epigenetic aging. GeroScience. 2022. 9. 9.Fahy GM, Brooke RT, Watson JP, Good Z, Vasanawala SS, Maecker H, et al. Reversal of epigenetic aging and immunosenescent trends in humans. Aging cell. 2019;18(6):e13028. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/acel.13028&link_type=DOI) 10. 10.Fiorito G, Caini S, Palli D, Bendinelli B, Saieva C, Ermini I, et al. DNA methylation[based biomarkers of aging were slowed down in a two[year diet and physical activity intervention trial: the DAMA study. Aging cell. 2021;20(10):e13439. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/acel.13439&link_type=DOI) 11. 11.Moqri M, Herzog C, Poganik JR, Justice J, Belsky DW, Higgins-Chen A, et al. Biomarkers of aging for the identification and evaluation of longevity interventions. Cell. 2023;186(18):3758–75. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2023.08.003&link_type=DOI) 12. 12.Labunskyy VM, Hatfield DL, Gladyshev VN. Selenoproteins: molecular pathways and physiological roles. Physiol Rev. 2014;94(3):739–77. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/physrev.00039.2013&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24987004&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F05%2F2024.04.04.24305314.atom) 13. 13.Xia Y, Hill KE, Li P, Xu J, Zhou D, Motley AK, et al. Optimization of selenoprotein P and other plasma selenium biomarkers for the assessment of the selenium nutritional requirement: a placebo-controlled, double-blind study of selenomethionine supplementation in selenium-deficient Chinese subjects. The American journal of clinical nutrition. 2010;92(3):525–31. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czo4OiI5Mi8zLzUyNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA0LzA1LzIwMjQuMDQuMDQuMjQzMDUzMTQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 14. 14.Demircan K, Bengtsson Y, Sun Q, Brange A, Vallon-Christersson J, Rijntjes E, et al. Serum selenium, selenoprotein P and glutathione peroxidase 3 as predictors of mortality and recurrence following breast cancer diagnosis: A multicentre cohort study. Redox Biol. 2021;47:102145. 15. 15.Schöttker B, Holleczek B, Hybsier S, Köhrle J, Schomburg L, Brenner H. Strong associations of serum selenoprotein P with all-cause mortality and mortality due to cancer, cardiovascular, respiratory and gastrointestinal diseases in older German adults. Eur J Epidemiol. 2024. 16. 16.Jujic A, Molvin J, Schomburg L, Hartmann O, Bergmann A, Melander O, et al. Selenoprotein P deficiency is associated with higher risk of incident heart failure. Free Radic Biol Med. 2023;207:11–6. 17. 17.Schomburg L, Orho-Melander M, Struck J, Bergmann A, Melander O. Selenoprotein-P Deficiency Predicts Cardiovascular Disease and Death. Nutrients. 2019;11(8). 18. 18.Qiu Z, Geng T, Wan Z, Lu Q, Guo J, Liu L, et al. Serum selenium concentrations and risk of all-cause and heart disease mortality among individuals with type 2 diabetes. The American journal of clinical nutrition. 2022;115(1):53–60. 19. 19.Bleys J, Navas-Acien A, Guallar E. Serum selenium levels and all-cause, cancer, and cardiovascular mortality among US adults. Arch Intern Med. 2008;168(4):404–10. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archinternmed.2007.74&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18299496&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F05%2F2024.04.04.24305314.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000253374600011&link_type=ISI) 20. 20.Cheng X, Wei Y, Wang R, Jia C, Zhang Z, An J, et al. Associations of essential trace elements with epigenetic aging indicators and the potential mediating role of inflammation. Redox Biology. 2023;67:102910. 21. 21.Xiao L, Zan G, Feng X, Bao Y, Huang S, Luo X, et al. The associations of multiple metals mixture with accelerated DNA methylation aging. Environmental Pollution. 2021;269:116230. 22. 22.Horvath S. DNA methylation age of human tissues and cell types. Genome biology. 2013;14(10):R115. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/gb-2013-14-10-r115&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24138928&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F05%2F2024.04.04.24305314.atom) 23. 23.Lu AT, Quach A, Wilson JG, Reiner AP, Aviv A, Raj K, et al. DNA methylation GrimAge strongly predicts lifespan and healthspan. Aging. 2019;11(2):303. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18632/aging.101684&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F05%2F2024.04.04.24305314.atom) 24. 24.Belsky DW, Caspi A, Corcoran DL, Sugden K, Poulton R, Arseneault L, et al. DunedinPACE, a DNA methylation biomarker of the pace of aging. eLife. 2022;11:e73420. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.73420&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35029144&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F05%2F2024.04.04.24305314.atom) 25. 25.Bertram L, Bockenhoff A, Demuth I, Duzel S, Eckardt R, Li SC, et al. Cohort profile: The Berlin Aging Study II (BASE-II). International journal of epidemiology. 2014;43(3):703–12. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyt018&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23505255&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F05%2F2024.04.04.24305314.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000338127000010&link_type=ISI) 26. 26.Hannum G, Guinney J, Zhao L, Zhang L, Hughes G, Sadda S, et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Molecular cell. 2013;49(2):359–67. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.molcel.2012.10.016&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23177740&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F05%2F2024.04.04.24305314.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000314379400018&link_type=ISI) 27. 27.Vetter VM, Meyer A, Karbasiyan M, Steinhagen-Thiessen E, Hopfenmuller W, Demuth I. Epigenetic clock and relative telomere length represent largely different aspects of aging in the Berlin Aging Study II (BASE-II). The journals of gerontology Series A, Biological sciences and medical sciences. 2018. 28. 28.Levine ME, Lu AT, Quach A, Chen BH, Assimes TL, Bandinelli S, et al. An epigenetic biomarker of aging for lifespan and healthspan. Aging. 2018;10(4):573. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18632/AGING.101414&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F05%2F2024.04.04.24305314.atom) 29. 29.Vetter VM, Spieker J, Sommerer Y, Buchmann N, Kalies CH, Regitz-Zagrosek V, et al. DNA methylation age acceleration is associated with risk of diabetes complications. Communications Medicine. 2023;3(1):21. 30. 30.Hughes DJ, Duarte-Salles T, Hybsier S, Trichopoulou A, Stepien M, Aleksandrova K, et al. Prediagnostic selenium status and hepatobiliary cancer risk in the European Prospective Investigation into Cancer and Nutrition cohort. The American journal of clinical nutrition. 2016;104(2):406–14. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czo5OiIxMDQvMi80MDYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNC8wNS8yMDI0LjA0LjA0LjI0MzA1MzE0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 31. 31.Stoppe C, McDonald B, Meybohm P, Christopher KB, Fremes S, Whitlock R, et al. Effect of high-dose selenium on postoperative organ dysfunction and mortality in cardiac surgery patients: the SUSTAIN CSX randomized clinical trial. JAMA surgery. 2023;158(3):235–44. 32. 32.Hybsier S, Schulz T, Wu Z, Demuth I, Minich WB, Renko K, et al. Sex-specific and inter-individual differences in biomarkers of selenium status identified by a calibrated ELISA for selenoprotein P. Redox Biology. 2017;11:403–14. 33. 33.Vetter VM, Drewelies J, Sommerer Y, Kalies CH, Regitz-Zagrosek V, Bertram L, et al. Epigenetic aging and perceived psychological stress in old age. Translational Psychiatry. 2022;12(1):410. 34. 34.Team RC. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria URL https://wwwR-projectorg. 2022. 35. 35.Yoshida K. Bartel, A. tableone: Create “Table 1” to Describe Baseline Characteristics with or without Propensity Score Weights, version 0.13. 2, R Studio package; RStudio. Inc: Vienna, Austria. 2022. 36. 36.Behrendt S. lm. beta Add standardized regression coefficients to linear-model-objects. R package version; 2023. 37. 37.Boyer K, Domingo-Relloso A, Jiang E, Haack K, Goessler W, Zhang Y, et al. Metal mixtures and DNA methylation measures of biological aging in American Indian populations. Environment International. 2023;178:108064. 38. 38.Sugden K, Caspi A, Elliott ML, Bourassa KJ, Chamarti K, Corcoran DL, et al. Association of pace of aging measured by blood-based DNA methylation with age-related cognitive impairment and dementia. Neurology. 2022;99(13):e1402–e13. 39. 39.Al-Mubarak AA, Grote Beverborg N, Suthahar N, Gansevoort RT, Bakker SJL, Touw DJ, et al. High selenium levels associate with reduced risk of mortality and new-onset heart failure: data from PREVEND. Eur J Heart Fail. 2022;24(2):299–307. 40. 40.Baker JR, Umesh S, Jenab M, Schomburg L, Tjønneland A, Olsen A, et al. Prediagnostic Blood Selenium Status and Mortality among Patients with Colorectal Cancer in Western European Populations. Biomedicines. 2021;9(11). 41. 41.Sandsveden M, Nilsson E, Borgquist S, Rosendahl AH, Manjer J. Prediagnostic serum selenium levels in relation to breast cancer survival and tumor characteristics. Int J Cancer. 2020;147(9):2424–36. 42. 42.Pietrzak S, Wójcik J, Scott RJ, Kashyap A, Grodzki T, Baszuk P, et al. Influence of the selenium level on overall survival in lung cancer. J Trace Elem Med Biol. 2019;56:46–51. 43. 43.Speckmann B, Grune T. Epigenetic effects of selenium and their implications for health. Epigenetics. 2015;10(3):179–90. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/15592294.2015.1013792&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25647085&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F05%2F2024.04.04.24305314.atom)