Vitamin D Supplementation is Associated with Slower Epigenetic Aging
====================================================================

* Valentin Max Vetter
* Yasmine Sommerer
* Christian Humberto Kalies
* Dominik Spira
* Lars Bertram
* Ilja Demuth

## Abstract

Adverse effects of low vitamin D level on mortality and morbidity are controversially discussed. Especially older people are at risk for vitamin D deficiency and therefore exposed to its potentially harmful influence on the aging process. A way of measuring differences in the biological age is through DNA methylation age (DNAm age) and its deviation from chronological age, DNAm age acceleration (DNAmAA). We previously reported on an association between vitamin D deficiency and higher 7-CpG DNAmAA in participants of the Berlin Aging Study II (BASE-II).

In this study, we employ a quasi-interventional study design to assess the relationship between DNAmAA of five epigenetic clocks and vitamin D supplementation. Longitudinal data were available for 1,036 participants of BASE-II that were reexamined on average 7.4 years later in the GendAge study (mean age at follow-up: 75.6 years, SD = 3.8 years, age range: 64.9 – 94.1 years, 51.9 % female). DNAmAA was estimated with the 7-CpG clock, Horvath’s clock, Hannum’s clock, PhenoAge and GrimAge. Methylation data were obtained through methylation-sensitive single nucleotide primer extension (MS-SNuPE) or Illumina’s Infinium “MethylationEPIC” array.

Vitamin D deficient participants who chose to start vitamin D supplementation after baseline examination showed a 2.6 year lower 7-CpG DNAmAA (p=0.011) and 1.3-year lower Horvath DNAmAA (p=0.042) compared to untreated and vitamin D deficient participants. DNAmAA did not statistically differ between participants with successfully treated vitamin D deficiency and healthy controls (p>0.16).

Therefore, we conclude that intake of vitamin D supplement is associated with lower DNAmAA in participants with vitamin D deficiency. Additionally, our findings suggest that sufficient vitamin D supplementation can compensate and potentially reverse the increase in 7-CpG DNAmAA that we found in in vitamin D deficient participants.

## Introduction

The health impact of vitamin D deficiency and benefits of its supplementation are subject to an ongoing controversial discussion [1]. Estimations suggest that 1 billion people worldwide have a serum vitamin D level below 50 nmol/L, the widely accepted cut-off for vitamin D deficiency [2]. Risk of vitamin D deficiency is particularly high for older adults because of their lower capability for cutaneous synthetization and less exposure to sunlight (homebound, institutionalized) [3]. Schöttker and colleagues found a constant decrease of on average 3 nmol/L serum vitamin D for every 10 years of age in a large population-based cohort (n=9,940, age range: 50-74 years) and therefore proposed vitamin D as a marker of aging [4].

The consequences of vitamin D deficiency are not fully understood. Prolonged and severe vitamin D deficiency as cause of rickets in children [5] and osteomalacia in adults [6] is uncontested. However, its role in prevention of falls [7-9] and fractures [2, 10-18] remains unclear. Besides its potential effect on musculoskeletal health, associations with numerous diseases such as respiratory tract infections [19, 20], depressive symptoms [21], rheumatoid arthritis [22, 23], multiple sclerosis [24], colorectal cancer [25, 26], breast cancer [27] and mortality [28] were reported, as well. These findings, however, are controversially discussed [29] and whether vitamin D deficiency is cause or result of these diseases is still unclear [30].

Considering the wide range of diseases associated with or promoted by vitamin D deficiency, a better understanding of its role in the process of healthy (biological) aging has far-reaching individual and socio-economic implications. Several biomarkers of aging are available that allow objectivation of biological age [31, 32]. The latest and often considered most promising [32] marker is the DNA methylation age (DNAm age) and its deviation from chronological age, DNAm age acceleration (DNAmAA). Different versions of epigenetic clocks that are used for its estimation are available. They differ mostly in number and location of analyzed cytosine-phosphate-guanin (CpG) sites as well as the domain of aging they represent best [33]. For most of these clocks the association between DNAm age and mortality (reviewed in [34]) and morbidity (reviewed in [32, 33]) is well documented. However, studies that examine the relationship between epigenetic age and vitamin D status are limited and inconclusive [4, 35, 36].

In this study we make use of two-wave longitudinal data on 1,036 older adults of the Berlin Aging Study II (BASE-II) that were reexamined on average 7.4 years later as part of the GendAge study. We previously reported on a cross-sectional association between vitamin D deficiency and high 7-CpG DNAmAA in this cohort at baseline [36]. We now repeated this analysis in the same participants at follow-up examination and compared it with DNAmAA obtained from four other epigenetic clocks, including GrimAge and PhenoAge [36]. Additionally, we employed a quasi-interventional (non-randomized) analysis scheme [37] to compare vitamin D deficient participants that chose to start vitamin D supplementation during the follow-up period with participants with “untreated” vitamin D deficiency at both time-points. Subsequently, we selected participants with sufficient vitamin D level to analyze DNAmAA in “treated” participants in comparison with “healthy” controls. We hypothesized, that the treatment of vitamin D deficiency would be associated with slower epigenetic aging and treated participants would not differ from participants without any signs for vitamin D deficiency.

## Methods

### Berlin Aging Study II (BASE-II) and GendAge Study

The Berlin Aging Study II (BASE-II) cohort consists of a convenience sample of 1,671 residents (age range: 60-85 years) of the greater metropolitan area of Berlin, Germany [38, 39] that was first examined between 2009 and 2014 (baseline, T0). The younger age group of BASE-II (n=500, age range: 20-37 years) has not been considered in this study.

On average 7.4 years later, 1,083 participants of the older group were re-examined as part of the GendAge study (follow-up, T1) [40]. Seventeen additional participants that were not assessed in the baseline examination were included at follow-up, resulting in a total dataset of 1,100 cases at follow-up.

All participants gave written informed consent. The medical assessments at baseline and follow-up were conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee of the Charité – Universitätsmedizin Berlin (approval numbers EA2/029/09 and EA2/144/16) and were registered in the German Clinical Trials Registry as DRKS00009277 and DRKS00016157.

### DNA Methylation Age (DNAm age)

We employed five different epigenetic clocks in this study. The 7-CpG clock [41] incorporates seven cytosine-phosphate-guanine (CpG) sites that were chosen to be most informative on chronological age and feasible to be measured by single nucleotide primer extension (SNuPE) [41-44]. Briefly, genomic DNA was isolated from whole blood samples and subsequently bisulfite converted. The segments of interest were amplified by a multiplex polymerase chain reaction (mPCR) and a single nucleotide primer extension (SNuPE) was performed.

Horvath’s clock [45], Hannum’s clock [46], PhenoAge [47] and GrimAge [48] were estimated from DNA methylation data that was obtained from the same DNA samples via Illumina’s “Infinium MethylationEPIC” array. Loading and quality control (QC) of the DNAm data was performed with the R package bigmelon [49] using default settings as described previously [50]. Probes that showed ≥1% of samples with a detection p-value of 0.05 or a bead count below three in >5% of samples were excluded from all analyses. Samples with a bisulfite conversion efficiency below 80% (estimated by the *bscon* function) and outliers (identified by the *outlyx* and *pcout* function [51]) were removed. After elimination of the outliers, the samples were reloaded and normalized with the *dasen* function. A root-mean-square deviation (RMSD) of ≥ 0.1 in DNAm-values after normalization, as determined by the *qual* function, led to exclusion of the sample from the dataset and loading and normalization were repeated with the new set of samples. Methylation age was calculated according to the manual of Steve Horvath’s website ([https://horvath.genetics.ucla.edu/html/dnamage/](https://horvath.genetics.ucla.edu/html/dnamage/)). Of the 513 CpGs of Levine’s clock and 71 CpGs of Hannum’s clock, 512 CpGs and 64 CpGs were available in this dataset. For a more detailed description of methods, please refer to ref. [52].

### DNA Methylation Age Acceleration (DNAmAA)

The deviation of DNAm age from chronological age, DNAm age acceleration (DNAmAA), was calculated as residuals of a linear regression analysis of DNAm age on chronological age and cell counts (neutrophils, monocytes, lymphocytes, eosinophils) [53]. Cell counts were measured by flow-cytometry in an accredited standard laboratory (LaborBerlin, Berlin, Germany).

### Vitamin D

Vitamin D levels were measured as 25-hydroxycholecalciferol (25-OHD) in blood serum samples that were taken during the same blood draw as the whole blood samples used for DNA isolation. A standard laboratory (LaborBerlin, Berlin, Germany) determined the serum vitamin D level through chemiluminescence immunoassays (DS-iSYS 25-hydroxvitamin D, Immunodiagnostic Systems, UK). Participants with serum 25-OHD levels <50nmol/L were classified as “deficient” ([54-56], reviewed in [2, 57]). Please note that there is no consensus on the optimal vitamin D level and numerous definitions of vitamin D “deficiency”, “insufficiency” and “inadequacy” exist [2].

Vitamin D supplementation was determined by searching the participants individual medication lists for the ATC-code for vitamin D and its analogs (“A11CC”). Additionally, participants were specifically asked for vitamin D supplementation at a different point during the examination. Only participants who gave consistent answers are kept in the supplemented (“treated”) group and cases with inconsistent answers were ignored during the group matching process.

### Quasi-interventional study design

A quasi-interventional study design was employed to analyze the effect of vitamin D supplementation on DNAmAA. Inclusion criteria for the “treated” group were: (i) vitamin D deficient at baseline, (ii) no vitamin D supplementation at baseline, (iii) start of vitamin D supplementation during follow-up period, (iv) sufficient vitamin D level at follow-up. To identify valid groups for comparison [58], optimal pair matching was performed based on chronological age, sex and a modified version of Charlson’s morbidity index [59, 60]. It was executed with R 3.6.2 [61] and the “MatchIt” package [62] which relies on functions from the “optmatch” package [63]. The “optimal pair matching” method selects matched samples with the smallest mean of the absolute pair distance across all pairs. It is very similar to the often used “nearest neighbor” method but reduces the distance within each pair better [64]. Absolute standardized mean differences before and after the matching process are displayed in Supplementary Figure 3 and 5. Both comparison groups (“untreated” and “healthy”) were determined according to the same procedure.

### Covariates

We included sex and chronological age (years since birth) as covariates in the regression analyses and the group matching process. To be able to adjust for seasonal changes in vitamin D level, we defined the period between October and April as “winter”. This seasonal stratification is in consistency with the analyses at baseline examination [36].

### Statistical Analyses

All statistical analyses and figures were produced with R 3.6.2 [61] and the packages “ggplot” [65] and “ggalluvial” [66]. Linear regression analyses were performed with the “lm” function. We defined participants whose DNAm age parameters differed more than 3 SD from the mean as outliers. The regarding cut-off values were calculated for each individual epigenetic clock parameter considered in this study. A p-value below 0.05 was considered statistically significant.

## Results

### Study population and descriptive statistics

The analyzed longitudinal cohort consists of 1,036 participants of BASE-II that provided information on their vitamin D status at baseline examination (mean age: 68.3 years, SD = 3.5 years) and at reexamination on average 7.4 years later (as part of the GendAge study, mean age: 75.6 years, SD = 3.8 years). Cohort characteristics are displayed in Table 1 and sex-stratified information is available in Supplementary Table 1 and 2.

View this table:
[Table 1:](http://medrxiv.org/content/early/2021/12/11/2021.12.10.21267597/T1)

Table 1: Cohort characteristics of the longitudinally analyzed sample of participants of the BASE-II and GendAge study.

DNAm age estimation differed between 0.01 years (GrimAge) and 14.4 years (Levine’s clock) from actual chronological age (Table 2). Correlation between DNAm age and chronological age was between Pearson’s r = 0.3 and Pearson’s r = 0.6. The highest correlation between clocks was found for Horvath’s and Hannum’s clock (Pearson’s r = 0.59 [67]). The test-retest correlation of the longitudinally available 7-CpG clock was in this cohort previously shown to be high (Pearson’s r=0.81) [67]. Participants aged epigenetically slower compared to their chronological age (mean slope of change = 0.75, SD = 0.64, [67]). A detailed cross-sectional and longitudinal description of the included epigenetic clocks can be found in ref. [67].

View this table:
[Table 2:](http://medrxiv.org/content/early/2021/12/11/2021.12.10.21267597/T2)

Table 2: Descriptive Statistics of the Five Available Epigenetic Clocks at follow-up.
DNAm age acceleration (DNAmAA) was calculated in years as residuals of a linear regression of DNAm age on chronological age and leukocyte cell count (neutrophils, monocytes, lymphocytes, eosinophils).

Of all longitudinally available participants, 46.4% had deficient vitamin D levels and 6.9% reported intake of vitamin D supplements at baseline. At follow-up we found 22.2% of all participants to be vitamin D deficient and the proportion of participants that used supplements increased by 13.5 percentage points to 20.4% (Table 1, Figure 1). This might be explained by a high fraction of participants that was examined during summer. Only for 55 participants the transition from deficient (at T0) to sufficient (at T1) vitamin D levels cannot be explained by supplementation or season of blood draw (Supplementary Figure 1).

![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/12/11/2021.12.10.21267597/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2021/12/11/2021.12.10.21267597/F1)

Figure 1: Alluvial plot of vitamin D status and information on supplementation and season of blood draw.
The individual trajectories of the 1,036 participants with longitudinal vitamin D data are displayed. Compared to baseline examination, 140 more participants reported intake of vitamin D supplements at follow-up.

Note: VitD supplement. = vitamin D supplementation, season = season of blood draw, T0 = baseline examination, T1 = follow-up examination

Of the 71 participants that reported vitamin D supplementation at baseline, 32 stopped within the follow up period. Therefore, 82% of the participants that reported vitamin D supplementation at follow-up started it after baseline examination. The deviation of detected vitamin D deficiency over calendar months during the follow-up examination period showed the expected higher frequency in winter (Supplementary Figure 2).

### Quasi-interventional study: Vitamin D supplementation is associated with lower DNAmAA

We previously reported on a 0.9 year higher 7-CpG DNAmAA in vitamin D deficient BASE-II participants compared to participants with a sufficient vitamin D level at baseline examination [36]. To assess whether this difference persisted after successful treatment of vitamin D deficiency, we employed a quasi-interventional analysis scheme in which the intervention (vitamin D supplementation) was non-randomly assigned by means of self-selection [37]. We compared vitamin D deficient participants that started vitamin D supplementation during follow-up time (and successfully treated vitamin D deficiency as indicated by sufficient serum vitamin level at follow-up) with those who did not and remained with deficient vitamin D level at follow-up. Sixty-three participants met the criteria to be included in the “treated” group (Figure 2 A and B). To minimize selection bias, optimal pair matching was used to select a control group of vitamin D deficient participants without vitamin D supplementation based on chronological age, sex, and morbidity index (Figure 2 A, Supplementary Figure 3).

![Figure 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/12/11/2021.12.10.21267597/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2021/12/11/2021.12.10.21267597/F2)

Figure 2: Overview on the longitudinal cohort (A) and of inclusion criteria for the “treated” group (B) and the “untreated” and “healthy” control group (B).
VitD suppl. = Vitamin D Supplementation; Supp. = Supplementation., suff. = sufficient, def. = deficient.

We found participants in the treated group to have 2.61 years lower 7-CpG DNAmAA (*t-test*, p=0.011) and 1.27 years lower Horvath DNAmAA (*t-test*, p=0.042) than the participants of the control group of untreated vitamin D deficient participants (Figure 3 B, Supplementary Figure 4). The association between vitamin D supplementation and 7-CpG DNAmAA remained significant after adjustment for covariates (season of blood draw, chronological age, sex) in a multiple linear regression analysis (β = −2.55, SE = 0.99, p = 0.011, n = 126, Supplementary Table 3). Sex-stratified subgroup analyses showed significant associations between 7-CpG DNAmAA and vitamin D supplementation in women (n=86) and between Horvath’s DNAmAA and vitamin D supplementation in men (n=40, Supplementary Table 3).

### Quasi-interventional study: participants with treated vitamin D deficiency do not differ in DNAmAA from vitamin D sufficient participants

We compared participants with successfully treated vitamin D deficiency with participants that reached vitamin D sufficiency without any supplementation at both baseline and follow-up (Figure 2 A and C). These participants were used as “healthy” control group to assess whether the vitamin D deficiency associated accelerated DNAmAA at baseline would still be present after treatment. The same matching process (Supplementary Figure 5) that we described above was employed to select the “healthy” control group (Figure 2 A). No significant difference in DNAmAA was found between participants with treated vitamin D deficiency and the healthy control group (Table 3).

View this table:
[Table 3:](http://medrxiv.org/content/early/2021/12/11/2021.12.10.21267597/T3)

Table 3: Comparison between participants with successfully treated vitamin D deficiency (“treated”), untreated vitamin D deficient participants and healthy controls.

### Cross-sectional association between vitamin D level and DNAmAA at follow-up

All participants who provided data on their vitamin D level at follow-up were included in the cross-sectional analyses. Too high serum vitamin D levels are harmful [68] and u-shaped association between vitamin D level and health parameters [69, 70] were reported before. Therefore, we excluded 34 participants with potentially harmful vitamin D levels of 125 nmol/L or above [68] from the cross-sectional regression analyses. Individual linear regression analyses of the DNAmAA of five different epigenetic clocks on serum vitamin D level and covariates (season of blood draw, chronological age, sex) were conducted.

We found a significant association between higher levels of vitamin D and 7-CpG DNAmAA (β = −0.02, SE = 0.01, p<0.05, n=1,008) and GrimAge DNAmAA (β = −0.01, SE = 0.004, p=0.03, Table 4). This associations seems to be mainly driven by women, as the all-women subgroup revealed these associations in sex-stratified analyses as well (Table 4).

View this table:
[Table 4:](http://medrxiv.org/content/early/2021/12/11/2021.12.10.21267597/T4)

Table 4: Multiple Linear Regression Analysis of DNAmAA of five different epigenetic clocks on Vitamin D level (nmol/L).
The model is adjusted for season of blood draw (winter/summer), chronological age (years), and sex (if applicable). Only participants with serum vitamin D levels below 125 nmol/L are included in these analyses.

## Discussion

In this study we report on a quasi-interventional analyses scheme that revealed significant associations between lower 7-CpG and Horvath DNAmAA and vitamin D supplementation in vitamin D deficient participants. Participants with successfully treated vitamin D deficiency had on average 2.6 years (7-CpG clock, p=0.011) and 1.3 years (Horvath’s clock, p=0.042) lower DNAmAA compared to the group of participants with untreated vitamin D deficiency. Further, we were able to show, that successfully treated participants do not differ in biological age from participants that have sufficient vitamin D levels.

Although vitamin D supplementation in the context of disease and mortality is subject to an intensive discussion, the body of literature analyzing its relationship to DNAmAA is limited. To our knowledge, only one study is available that employs an interventional design. Chen and colleagues conducted a small randomized controlled trial with 51 overweight and vitamin D deficient African Americans (26.1 ± 9.3 years old) that were assigned to four groups (600 IU/d, 2000 IU/d, 4000 IU/d vitamin D3 and placebo). Participants that took 4000 IU/d for 16-weeks showed 1.85 years decrease in Horvath aging compared to the placebo group (p=0.046). A 1.9-year decrease in Hannum aging was found in the group that took 2000 IU per day (p=0.044) [35]. Two cross-sectional studies analyzing the relationship between DNAmAA and vitamin D are available. Schöttker and colleagues found no cross-sectional association between Horvath’s or Hannum’s clock and vitamin D levels in a large cohort of 9,940 participants between 50 and 74 years [4]. In contrast to these findings, we reported on a cross-sectional association between lower 7-CpG DNAmAA and vitamin D sufficiency before. These cross-sectional findings [36] were repeatable in the same participants at follow-up on average 7.4 years later. The additional analyses with DNAmAA estimated by four epigenetic clocks, that rely on epigenomewide methylation data and are available only for follow-up examination, revealed a cross-sectional association between low GrimAge DNAmAA and high vitamin D level (β = −0.01, SE = 0.004, p=0.03). These associations were found in the all-female but not in the all-male subgroups. Men and women seem to epigenetically age at a different pace which was shown in this study (Table 2 and ref. [67, 71]) as well as others [72-74]. Therefore, sex differences in the association between DNAmAA and vitamin D level seem plausible.

Furthermore, it is known that the epigenetic clocks analyzed here differ in the aspects of aging they reflect [33]. This presumably results from the difference in CpG sites included and whether the epigenetic clocks were trained to predict chronological age (first-generation clocks) or phenotypic measures (“second-generation clocks”). Therefore, differences in association with vitamin D between different clocks are expected and help to further differentiate the areas of aging they represent best.

This study has several limitations we want to point out. First, the participants of this study are above average health and the proportion of vitamin D deficient participants (22.2%) is comparatively low. In a large study of the Robert-Koch Institute (DEGS, n=6,995, age range: 18-79 years), 69.9% of men and 62.6% of women between 65 and 79 years had a vitamin D level below 50 nmol/L [75]. The high proportion of participants with sufficient vitamin D levels at follow-up of our study might in part be explained by the high percentage of vitamin D supplemented participants. The frequency of participants with vitamin D supplementation increased by 13.5 percentage points at baseline to 20.37% at follow-up. This rise might be partly explained by the increasing popularity of vitamin D supplements that are easily accessible and comparatively cheap. At baseline examination each participant was given a medical report that included information on the vitamin D level, which might have added to this effect. Therefore, it would be of great interest to repeat these analyses in sample sets with a higher proportion of vitamin D deficient participants. Second, although the sample size for the cross-sectional analyses is comparatively large (n=1,008), the number of participants that meets the criteria to be included in the “treated” group” is rather small. We therefore cannot rule out that missing statistical significance in our study reflects a lack of power. Third, a limitation of this study is the missing information on the exact doses of vitamin D supplementation and duration of intake.

However, the dose for optimal treatment depends on several factors to be considered such as the initial vitamin D level, time spent in the sun, diet and chronological age and needs to be determined individually and to be adjusted if conditions change. Therefore, we used successful treatment (controlled by serum vitamin D level) rather than a specific dose as definition for the “treated” group. A cautious approach to vitamin D supplementation is important because uncontrolled intake of high doses of vitamin D supplements can lead to severe adverse effects [76]. Vitamin D serum levels above 374 nmol/L are associated with hyperphosphatemia and hypercalcemia [2] which can lead to renal failure and cardiac arrest [57]. The Institute of Medicine reported that a vitamin D serum level above 125 nmol/L over a long time period “should raise concerns among clinicians about potential adverse effects” [68]. Therefore, the supplementation of vitamin D, although at least to some extent undoubtfully useful, should be done attentively and only after medical consultation to avoid excessive (as well as insufficient) treatment. Fourth, no causal conclusion can be drawn from the analyses conducted here because we cannot confidently assume that all possible (known and unknown) covariates were considered. It is likely that some of the participants who started vitamin D supplementation during the follow-up period also tried to improve their lifestyle in other ways. As these possible lifestyle changes are not known and therefore not controlled for, they could confound the results reported in this study. To minimize potential confounding effects and selection bias we applied optimal pair matching with carefully selected variables to form valid comparison groups [58]. Nevertheless, only randomized controlled trials are able to control for all confounding factors (known or unknown). Fifth, due to the quasi-interventional study design we do not have a placebo control group. Again, this would be needed to fully understand the potential causal effect of vitamin D on DNAmAA but can only be obtained from randomized controlled trials. Sixth, due to the very limited knowledge of the pathways leading to age related epigenetic changes, we do not know whether vitamin D has an actual influence on biological aging or rather just modifies epigenetic signatures and thus interferes with the measurement accuracy of DNAmAA as an aging marker. The current body of literature makes the former most likely [77, 78], but only sufficiently sized randomized controlled trials can clarify this assumption. Seventh, in line with most other studies that evaluate more than one epigenetic clock, we did not adjust our analyses for multiple testing. Therefore, an increase in the rate of false-positive findings cannot be ruled out and future work needs to independently validate our findings before any further reaching conclusions can be drawn.

Strengths of this study include its longitudinal design and the large, well characterized cohort. It allows to test for associations between vitamin D supplementation and DNAmAA in well balanced intervention and control groups. We compare results across five different epigenetic clocks that are known to represent different domains of aging. Therefore, this study helps to distinguish the domains that are represented best by the different versions of DNAmAA and adds valuable information on the consequences of vitamin D deficiency and how it possibly could be reversed.

## Conclusion

We found DNAmAA to be on average 2.6 years (7-CpG clock, p = 0.011) and 1.3 years (Horvath clock, p = 0.042) lower in vitamin D deficient participants who were successfully supplemented compared to vitamin D-deficient participants without supplementation. Additionally, no difference in DNAmAA of all five clocks was found between treated participants and the control group of participants who are vitamin D sufficient without supplementation. Therefore, a sufficient supplementation of vitamin D-deficient patients seems to be beneficial with regards to epigenetic aging, at least when measured with the 7-CpG or Horvath’s clock. Sufficiently sized randomized controlled trials need to further investigate this potentially causal effect.

## Supporting information

Supplement [[supplements/267597_file03.pdf]](pending:yes)

## Data Availability

Data are available upon reasonable request. Interested investigators are invited to contact the study coordinating PI Ilja Demuth at ilja.demuth{at}charite.de to obtain additional information about the GendAge study and the data-sharing application form.

## Funding statement

This work was supported by grants of the Deutsche Forschungsgemeinschaft (grant number DE 842/7-1 to ID), the ERC (as part of the “Lifebrain” project to LB), and the Cure Alzheimer’s Fund (as part of the “CIRCUITS” consortium to LB). This article uses data from the Berlin Aging Study II (BASE-II) and the GendAge study which were supported by the German Federal Ministry of Education and Research under grant numbers #01UW0808; #16SV5536K, #16SV5537, #16SV5538, #16SV5837, #01GL1716A and #01GL1716B. We thank all probands of the BASE-II/GendAge study for their participation in this research.

## Author contributions

Conceived and designed the study: VMV and ID. Contributed study specific data: all authors. Analyzed the data: VMV and YS. Wrote the manuscript: VMV and ID. All authors revised and approved the manuscript.

## Conflict of interest

None declared.

## Data availability statement

Data are available upon reasonable request. Interested investigators are invited to contact the study coordinating PI Ilja Demuth at ilja.demuth{at}charite.de to obtain additional information about the GendAge study and the data-sharing application form.

## Acknowledgments

-

*   Received December 10, 2021.
*   Revision received December 10, 2021.
*   Accepted December 11, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  1.Spiro, A. and  J.L. Buttriss, Vitamin D: An overview of vitamin D status and intake in Europe. Nutr Bull, 2014. 39(4): p. 322–350.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/nbu.12108&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25635171&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 

2.  2.Holick, M.F., Vitamin D deficiency. N Engl J Med, 2007. 357(3): p. 266–81.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMra070553&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17634462&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248115100008&link_type=ISI) 

3.  3.Holick, M.F. High prevalence of vitamin D inadequacy and implications for health. In Mayo Clinic Proceedings. 2006. Elsevier.
    
    

4.  4.Schöttker, B., et al., Serum 25-hydroxyvitamin D levels as an aging marker: strong associations with age and all-cause mortality independent from telomere length, epigenetic age acceleration, and 8-isoprostane levels. The Journals of Gerontology: Series A, 2019. 74(1): p. 121–128.
    
    

5.  5.Pettifor, J.M.,  K. Thandrayen, and  T.D. Thacher, Vitamin D deficiency and nutritional rickets in children, in Vitamin D. 2018, Elsevier. p. 179–201.
    
    

6.  6.Aaron, J., et al., Frequency of osteomalacia and osteoporosis in fractures of the proximal femur. The Lancet, 1974. 303(7851): p. 229–233.
    
    

7.  7.Bischoff-Ferrari, H., et al., Estimation of optimal serum concentrations of 25-hydroxyvitamin D for multiple health outcomes (vol 84, pg 18, 2006). American Journal of Clinical Nutrition, 2006. 84(5): p. 1253–1253.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czo5OiI4NC81LzEyNTMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8xMi8xMS8yMDIxLjEyLjEwLjIxMjY3NTk3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

8.  8.Bischoff-Ferrari, H.A., et al., Estimation of optimal serum concentrations of 25-hydroxyvitamin D for multiple health outcomes. The American journal of clinical nutrition, 2006. 84(1): p. 18–28.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czo3OiI4NC8xLzE4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMTEvMjAyMS4xMi4xMC4yMTI2NzU5Ny5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

9.  9.Bolland, M.J., et al., Vitamin D supplementation and falls: a trial sequential meta-analysis. The lancet Diabetes & endocrinology, 2014. 2(7): p. 573–580.
    
    

10. 10.Reid, I.R.,  M.J. Bolland, and  A. Grey, Effects of vitamin D supplements on bone mineral density: a systematic review and meta-analysis. The Lancet, 2014. 383(9912): p. 146–155.
    
    

11. 11.Bischoff-Ferrari, H.A., et al., A pooled analysis of vitamin D dose requirements for fracture prevention. N Engl J Med, 2012. 367: p. 40–49.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1109617&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22762317&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000305979400004&link_type=ISI) 

12. 12.Bolland, M.J., et al., The effect of vitamin D supplementation on skeletal, vascular, or cancer outcomes: a trial sequential meta-analysis. The lancet Diabetes & endocrinology, 2014. 2(4): p. 307–320.
    
    

13. 13.Group, D., Patient level pooled analysis of 68 500 patients from seven major vitamin D fracture trials in US and Europe. BMJ, 2011. 343: p. d5245.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzNDMvYXVnMzBfMS9kNTI0NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEyLzExLzIwMjEuMTIuMTAuMjEyNjc1OTcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

14. 14.Group, D., Patient level pooled analysis of 68 500 patients from seven major vitamin D fracture trials in US and Europe. The BMJ, 2010. 340.
    
    

15. 15.Grossman, D.C., et al., Vitamin D, Calcium, or Combined Supplementation for the Primary Prevention of Fractures in Community-Dwelling Adults: US Preventive Services Task Force Recommendation Statement. Jama, 2018. 319(15): p. 1592–1599.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2018.3185&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 

16. 16.Weaver, C.M., et al., Calcium plus vitamin D supplementation and risk of fractures: an updated meta-analysis from the National Osteoporosis Foundation. Osteoporos Int, 2016. 27(1): p. 367–76.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 

17. 17.Zhao, J.G., et al., Association Between Calcium or Vitamin D Supplementation and Fracture Incidence in Community-Dwelling Older Adults: A Systematic Review and Meta-analysis. Jama, 2017. 318(24): p. 2466–2482.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2017.19344&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29279934&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 

18. 18.Bischoff-Ferrari, H.A., et al., Vitamin D supplementation and musculoskeletal health. Lancet Diabetes Endocrinol, 2019. 7(2): p. 85.
    
    

19. 19.Ginde, A.A.,  J.M. Mansbach, and  C.A. Camargo, Association between serum 25-hydroxyvitamin D level and upper respiratory tract infection in the Third National Health and Nutrition Examination Survey. Archives of internal medicine, 2009. 169(4): p. 384–390.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archinternmed.2008.560&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19237723&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000263599600009&link_type=ISI) 

20. 20.Bergman, P., et al., Vitamin D and respiratory tract infections: a systematic review and meta-analysis of randomized controlled trials. PloS one, 2013. 8(6): p. e65835.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0065835&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23840373&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 

21. 21.Anglin, R.E., et al., Vitamin D deficiency and depression in adults: systematic review and meta-analysis. The British journal of psychiatry, 2013. 202(2): p. 100–107.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImJqcHJjcHN5Y2giO3M6NToicmVzaWQiO3M6OToiMjAyLzIvMTAwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMTEvMjAyMS4xMi4xMC4yMTI2NzU5Ny5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

22. 22.Merlino, L.A., et al., Vitamin D intake is inversely associated with rheumatoid arthritis: results from the Iowa Women’s Health Study. Arthritis & Rheumatism: Official Journal of the American College of Rheumatology, 2004. 50(1): p. 72–77.
    
    

23. 23.Salesi, M. and  Z. Farajzadegan, Efficacy of vitamin D in patients with active rheumatoid arthritis receiving methotrexate therapy. Rheumatology international, 2012. 32(7): p. 2129–2133.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21523344&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 

24. 24.Soilu-Hänninen, M., et al., A longitudinal study of serum 25-hydroxyvitamin D and intact parathyroid hormone levels indicate the importance of vitamin D and calcium homeostasis regulation in multiple sclerosis. Journal of Neurology, Neurosurgery & Psychiatry, 2008. 79(2): p. 152–157.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiam5ucCI7czo1OiJyZXNpZCI7czo4OiI3OS8yLzE1MiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEyLzExLzIwMjEuMTIuMTAuMjEyNjc1OTcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

25. 25.Feskanich, D., et al., Plasma vitamin D metabolites and risk of colorectal cancer in women. Cancer Epidemiology and Prevention Biomarkers, 2004. 13(9): p. 1502–1508.
    
    

26. 26.Wei, M.Y., et al., Vitamin D and prevention of colorectal adenoma: a meta-analysis. Cancer Epidemiology and Prevention Biomarkers, 2008. 17(11): p. 2958–2969.
    
    

27. 27.Garland, C.F., et al., The role of vitamin D in cancer prevention. American journal of public health, 2006. 96(2): p. 252–261.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2105/AJPH.2004.045260&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16380576&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000234991400014&link_type=ISI) 

28. 28.Bjelakovic, G., et al., Vitamin D supplementation for prevention of mortality in adults. Cochrane Database of Systematic Reviews, 2014(1).
    
    

29. 29.Allan, G.M., et al., Vitamin D: A Narrative Review Examining the Evidence for Ten Beliefs. Journal of General Internal Medicine, 2016. 31(7): p. 780–791.
    
    

30. 30.Autier, P., et al., Vitamin D status and ill health: a systematic review. The lancet Diabetes & endocrinology, 2014. 2(1): p. 76–89.
    
    

31. 31.Lopez-Otin, C., et al., The hallmarks of aging. Cell, 2013. 153(6): p. 1194–217.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2013.05.039&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23746838&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000319979200007&link_type=ISI) 

32. 32.Jylhava, J.,  N.L. Pedersen, and  S. Hagg, Biological Age Predictors. EBioMedicine, 2017. 21: p. 29–36.
    
    

33. 33.Horvath, S. and  K. Raj, DNA methylation-based biomarkers and the epigenetic clock theory of ageing. Nature Reviews Genetics, 2018: p. 1.
    
    

34. 34.Fransquet, P.D., et al., The epigenetic clock as a predictor of disease and mortality risk: a systematic review and meta-analysis. Clinical epigenetics, 2019. 11(1): p. 62.
    
    

35. 35.Chen, L., et al., Effects of Vitamin D3 Supplementation on Epigenetic Aging in Overweight and Obese African Americans With Suboptimal Vitamin D Status: A Randomized Clinical Trial. J Gerontol A Biol Sci Med Sci, 2019. 74(1): p. 91–98.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/gerona/gly223&link_type=DOI) 

36. 36.Vetter, V.M., et al., Epigenetic clock and leukocyte telomere length are associated with vitamin D status, but not with functional assessments and frailty in the Berlin Aging Study II. The Journals of Gerontology: Series A, 2020.
    
    

37. 37.Cook, T.D.,  D.T. Campbell, and  W. Shadish, Experimental and quasi-experimental designs for generalized causal inference. 2002: Houghton Mifflin Boston, MA.
    
    

38. 38.Bertram, L., et al., Cohort profile: The Berlin Aging Study II (BASE-II). Int J Epidemiol, 2014. 43(3): p. 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%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000338127000010&link_type=ISI) 

39. 39.Gerstorf, D., et al., Editorial. Gerontology, 2016. 62(3): p. 311–5.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000441495&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 

40. 40.Demuth, I., et al., Cohort profile: follow-up of a Berlin Aging Study II (BASE-II) subsample as part of the GendAge study. BMJ Open, 2021. 11(6): p. e045576.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYm1qb3BlbiI7czo1OiJyZXNpZCI7czoxMjoiMTEvNi9lMDQ1NTc2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMTEvMjAyMS4xMi4xMC4yMTI2NzU5Ny5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

41. 41.Vetter, V.M., et al., Epigenetic clock and relative telomere length represent largely different aspects of aging in the Berlin Aging Study II (BASE-II). J Gerontol A Biol Sci Med Sci, 2018.
    
    

42. 42.Vidal-Bralo, L.,  Y. Lopez-Golan, and  A. Gonzalez, Simplified Assay for Epigenetic Age Estimation in Whole Blood of Adults. Front Genet, 2016. 7: p. 126.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fgene.2016.00126&link_type=DOI) 

43. 43.Vidal-Bralo, L.,  Y. Lopez-Golan, and  A. Gonzalez, Corrigendum: simplified assay for epigenetic age estimation in whole blood of adults. Frontiers in genetics, 2017. 8: p. 51.
    
    

44. 44.Kaminsky, Z.A., et al., Single nucleotide extension technology for quantitative site-specific evaluation of metC/C in GC-rich regions. Nucleic Acids Res, 2005. 33(10): p. e95.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gni094&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15958788&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 

45. 45.Horvath, S., DNA methylation age of human tissues and cell types. Genome Biol, 2013. 14(10): p. 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%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 

46. 46.Hannum, G., et al., Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell, 2013. 49(2): p. 359–367.
    
    [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%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000314379400018&link_type=ISI) 

47. 47.Levine, M.E., et al., An epigenetic biomarker of aging for lifespan and healthspan. Aging (Albany NY), 2018. 10(4): p. 573.
    
    

48. 48.Lu, A.T., et al., DNA methylation GrimAge strongly predicts lifespan and healthspan. Aging (Albany NY), 2019. 11(2): p. 303.
    
    

49. 49.Gorrie-Stone, T.J., et al., Bigmelon: tools for analysing large DNA methylation datasets. Bioinformatics, 2019. 35(6): p. 981–986.
    
    

50. 50.Sommerer, Y., et al., Entorhinal cortex EWAS meta-analysis highlights four novel loci showing differential methylation in Alzheimer’s disease. bioRxiv, 2021: p. 2021.07.02.450878.
    
    

51. 51.Filzmoser, P.,  R. Maronna, and  M. Werner, Outlier identification in high dimensions. Computational statistics & data analysis, 2008. 52(3): p. 1694–1711.
    
    

52. 52.Vetter, V.M., et al., Seven-CpG DNA Methylation Age determined by Single Nucleotide Primer Extension and Illumina’s Infinium MethylationEPIC array provide highly comparable results. bioRxiv, 2021: p. 2021.08.13.456213.
    
    

53. 53.Quach, A., et al., Epigenetic clock analysis of diet, exercise, education, and lifestyle factors. Aging (Albany NY), 2017. 9(2): p. 419–446.
    
    

54. 54.Organization, W.H., Prevention and management of osteoporosis. World Health Organ Tech Rep Ser, 2003. 921: p. 1–164.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15293701&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 

55. 55.Society, G.N., New reference values for vitamin D. 2012.
    
    

56. 56.Ross, A.C., et al., Institute of Medicine (US) Committee to Review Dietary Reference Intakes for Vitamin D and Calcium. Dietary reference intakes for calcium and vitamin D. 2011, Washington, DC: National Academies Press.
    
    

57. 57.Spiro, A. and  J. Buttriss, Vitamin D: an overview of vitamin D status and intake in E urope. Nutrition bulletin, 2014. 39(4): p. 322–350.
    
    

58. 58.White, H. and  S. Sabarwal, Quasi-experimental design and methods. Methodological briefs: impact evaluation, 2014. 8: p. 1–16.
    
    

59. 59.Charlson, M.E., et al., A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. J Chronic Dis, 1987. 40(5): p. 373–83.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0021-9681(87)90171-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3558716&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1987G855900002&link_type=ISI) 

60. 60.Meyer, A., et al., Leukocyte telomere length is related to appendicular lean mass: cross-sectional data from the Berlin Aging Study II (BASE-II). Am J Clin Nutr, 2016. 103(1): p. 178–83.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czo5OiIxMDMvMS8xNzgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8xMi8xMS8yMDIxLjEyLjEwLjIxMjY3NTk3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

61. 61.Team, R.C., R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL [https://www.R-project.org](https://www.R-project.org), 2017.
    
    

62. 62.Ho, D.E., et al., MatchIt: nonparametric preprocessing for parametric causal inference. J Stat Softw, 2011. 42(8): p. 1–28.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/JSS.V042.I08&link_type=DOI) 

63. 63.Hansen, B.B. and  S.O. Klopfer, Optimal full matching and related designs via network flows. Journal of computational and Graphical Statistics, 2006. 15(3): p. 609–627.
    
    

64. 64.Gu, X.S. and  P.R. Rosenbaum, Comparison of multivariate matching methods: Structures, distances, and algorithms. Journal of Computational and Graphical Statistics, 1993. 2(4): p. 405–420.
    
    

65. 65.Wickham, H., Elegant graphics for data analysis (ggplot2). 2009, New York, NY: Springer-Verlag.
    
    

66. 66.Brunson, J., ggalluvial: Alluvial Plots in ‘ggplot2’. R package version 0.11. 3. 2020.
    
    

67. 67.Vetter, V.M., et al., Relationship between five Epigenetic Clocks, Telomere Length and Functional Capacity assessed in Older Adults: Cross-sectional and Longitudinal Analyses. medRxiv, 2021: p. 2021.10.05.21264547.
    
    

68. 68.Ross, A.C., et al., The 2011 report on dietary reference intakes for calcium and vitamin D from the Institute of Medicine: what clinicians need to know. The Journal of Clinical Endocrinology & Metabolism, 2011. 96(1): p. 53–58.
    
    

69. 69.Anagnostis, P., et al., Vitamin D supplementation and fracture risk: Evidence for a U-shaped effect. Maturitas, 2020. 141: p. 63–70.
    
    

70. 70.Mayne, S.T.,  L.M. Ferrucci, and  B. Cartmel, Lessons Learned from Randomized Clinical Trials of Micronutrient Supplementation for Cancer Prevention. Annual Review of Nutrition, 2012. 32(1): p. 369–390.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev-nutr-071811-150659&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22524186&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000307963100018&link_type=ISI) 

71. 71.Lemke, E., et al., Cardiovascular health is associated with the epigenetic clock in the Berlin Aging Study II (BASE-II). Mechanisms of Ageing and Development, 2021: p. 111616.
    
    

72. 72.Horvath, S., et al., An epigenetic clock analysis of race/ethnicity, sex, and coronary heart disease. Genome biology, 2016. 17(1): p. 1–23.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-015-0866-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26753840&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 

73. 73.Simpkin, A.J., et al., Prenatal and early life influences on epigenetic age in children: a study of mother–offspring pairs from two cohort studies. Human molecular genetics, 2016. 25(1): p. 191–201.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddv456&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26546615&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 

74. 74.Banszerus, V.L., et al., Exploring the relationship of relative telomere length and the epigenetic clock in the LipidCardio cohort. International journal of molecular sciences, 2019. 20(12): p. 3032.
    
    

75. 75.Rabenberg, M. and  G. Mensink, Vitamin-D-Status in Deutschland. 2016, Robert Koch-Institut, Epidemiologie und Gesundheitsberichterstattung.
    
    

76. 76.Vieth, R., Critique of the considerations for establishing the tolerable upper intake level for vitamin D: critical need for revision upwards. The Journal of nutrition, 2006. 136(4): p. 1117–1122.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToibnV0cml0aW9uIjtzOjU6InJlc2lkIjtzOjEwOiIxMzYvNC8xMTE3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMTEvMjAyMS4xMi4xMC4yMTI2NzU5Ny5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

77. 77.Zhu, H., et al., A genome-wide methylation study of severe vitamin D deficiency in African American adolescents. The Journal of pediatrics, 2013. 162(5): p. 1004–1009. e1.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jpeds.2012.10.059&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23219444&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F11%2F2021.12.10.21267597.atom) 

78. 78.Beckett, E.L., et al., Relationship between methylation status of vitamin D-related genes, vitamin D levels, and methyl-donor biochemistry. Journal of Nutrition & Intermediary Metabolism, 2016. 6: p. 8–15.