DNA methylome profiling identifies novel methylated genes in African American patients with colorectal neoplasia

The identification of genes that are differentially methylated in colorectal cancer (CRC) has potential value for both diagnostic and therapeutic interventions specifically in high-risk populations such as African Americans (AAs). However, DNA methylation patterns in CRC, especially in AAs, have not been systematically explored and remain poorly understood. Here, we performed DNA methylome profiling to identify the methylation status of CpG islands within candidate genes involved in critical pathways important in the initiation and development of CRC. We used reduced representation bisulfite sequencing (RRBS) in colorectal cancer and adenoma tissues that were compared with DNA methylome from a healthy AA subject’s colon tissue and peripheral blood DNA. The identified methylation markers were validated in fresh frozen CRC tissues and corresponding normal tissues from AA patients diagnosed with CRC at Howard University Hospital. We identified and validated the methylation status of 355 CpG sites located within 16 gene promoter regions associated with CpG islands. Fifty CpG sites located within CpG islands—in genes ATXN7L1 (2), BMP3 (7), EID3 (15), GAS7 (1), GPR75 (24), and TNFAIP2 (1)—were significantly hypermethylated in tumor vs. normal tissues (P < 0.05). The methylation status of BMP3, EID3, GAS7, and GPR75 was confirmed in an independent, validation cohort. Ingenuity pathway analysis mapped three of these markers (GAS7, BMP3 and GPR) in the insulin and TGF-β1 network—the two key pathways in CRC. In addition to hypermethylated genes, our analysis also revealed that LINE-1 repeat elements were progressively hypomethylated in the normal-adenoma-cancer sequence. We conclude that DNA methylome profiling based on RRBS is an effective method for screening aberrantly methylated genes in CRC. While previous studies focused on the limited identification of hypermethylated genes, ours is the first study to systematically and comprehensively identify novel hypermethylated genes, as well as hypomethylated LINE-1 sequences, which may serve as potential biomarkers for CRC in African Americans. Our discovered biomarkers were intimately linked to the insulin/TGF-B1 pathway, further strengthening the association of diabetic disorders with colon oncogenic transformation.     

Prenatal antidepressant exposure associated with CYP2E1 DNA methylation change in neonates

Some but not all neonates are affected by prenatal exposure to serotonin reuptake inhibitor antidepressants (SRI) and maternal mood disturbances. Distinguishing the impact of these 2 exposures is challenging and raises critical questions about whether pharmacological, genetic, or epigenetic factors can explain the spectrum of reported outcomes. Using unbiased DNA methylation array measurements followed by a detailed candidate gene approach, we examined whether prenatal SRI exposure was associated with neonatal DNA methylation changes and whether such changes were associated with differences in birth outcomes. Prenatal SRI exposure was first associated with increased DNA methylation status primarily at CYP2E1(β_Non-exposed = 0.06, β_SRI-exposed = 0.30, FDR = 0); however, this finding could not be distinguished from the potential impact of prenatal maternal depressed mood. Then, using pyrosequencing of CYP2E1 regulatory regions in an expanded cohort, higher DNA methylation status—both the mean across 16 CpG sites (P < 0.01) and at each specific CpG site (P < 0.05)—was associated with exposure to lower 3rd trimester maternal depressed mood symptoms only in the SRI-exposed neonates, indicating a maternal mood x SRI exposure interaction. In addition, higher DNA methylation levels at CpG2 (P = 0.04), CpG9 (P = 0.04) and CpG10 (P = 0.02), in the interrogated CYP2E1 region, were associated with increased birth weight independently of prenatal maternal mood, SRI drug exposure, or gestational age at birth. Prenatal SRI antidepressant exposure and maternal depressed mood were associated with altered neonatal CYP2E1 DNA methylation status, which, in turn, appeared to be associated with birth weight.     

DNA methylation provides insight into intergenerational risk for preterm birth in African Americans

African Americans are at increased risk for spontaneous preterm birth (PTB). Though PTB is heritable, genetic studies have not identified variants that account for its intergenerational risk, prompting the hypothesis that epigenetic factors may also contribute. The objective of this study was to evaluate DNA methylation from maternal leukocytes to identify patterns specific to PTB and its intergenerational risk. DNA from peripheral leukocytes from African American women that delivered preterm (24–34 weeks; N = 16) or at term (39–41 weeks; N = 24) was assessed for DNA methylation using the HumanMethylation450 BeadChip. In maternal samples, 17,829 CpG sites associated with PTB, but no CpG site remained associated after correction for multiple comparisons. Examination of paired maternal-fetal samples identified 5,171 CpG sites in which methylation of maternal samples correlated with methylation of her respective fetus (FDR < 0.05). These correlated sites were enriched for association with PTB in maternal leukocytes. The majority of correlated CpG sites could be attributed to one or more genetic variants. They were also significantly more likely to be in genes involved in metabolic, cardiovascular, and immune pathways, suggesting a role for genetic and environmental contributions to PTB risk and chronic disease. The results of this study may provide insight into the factors underlying intergenerational risk for PTB and its consequences.     

DNA methylation profiles in African American prostate cancer patients in relation to disease progression

This study examined whether differential DNA methylation is associated with clinical features of more aggressive disease at diagnosis and prostate cancer recurrence in African American men, who are more likely to die from prostate cancer than other populations. Tumor tissues from 76 African Americans diagnosed with prostate cancer who had radical prostatectomy as their primary treatment were profiled for epigenome-wide DNA methylation levels. Long-term follow-up identified 19 patients with prostate cancer recurrence. Twenty-three CpGs were differentially methylated (FDR q ≤ 0.25, mean methylation difference ≥ 0.10) in patients with vs. without recurrence, including CpGs in GCK, CDKL2, PRDM13, and ZFR2. Methylation differences were also observed between men with metastatic-lethal prostate cancer vs. no recurrence (five CpGs), regional vs. local pathological stage (two CpGs), and higher vs. lower tumor aggressiveness (one CpG). These results indicate that differentially methylated CpG sites identified in tumor tissues of African American men may contribute to prostate cancer aggressiveness.     

Race/Ethnicity-Specific Association of Vitamin D and Global DNA Methylation: Cross-Sectional and Interventional Findings

ObjectivesUnderstanding of the influence of vitamin D deficiency on epigenome will provide novel insights into the chronic disease risk. We tested our hypotheses that 1) vitamin D deficiency is associated with global hypomethylation and this association may be race/ethnicity dependent; and 2) vitamin D supplementation will increase global DNA methylation level.MethodsA two-stage design, cross-sectional observation followed by a 16 week randomized, double- blinded, placebo-controlled trial (RCT) of vitamin D3 supplementation, was undertaken. Global DNA methylation level (percentage of 5-methylcytosine, %5-mC) was quantified using leukocyte DNA with the MethylFlash^TM Methylated DNA Quantification kit (Epigentek). Global methylation data was obtained from 454 Caucasians and African Americans (42%) in the observation cohort and 58 African Americans with vitamin D deficiency in the dose responsive RCT.ResultsIn the cross-sectional study, African Americans had lower %5-mC than Caucasians (P = 0.04). A significant interaction was detected between plasma 25(OH)D and race on %5-mC (P = 0.05), as a positive association was observed between plasma 25(OH)D and %5-mC in African Americans (β = 0.20, p<0.01), but not in Caucasians (β = 0.03, p = 0.62). In the 16-week RCT, a dose-response benefit of vitamin D₃ supplementation was observed for %5-mC, as indicated by a significant linear upward trend (-0.01 ± 0.01%, placebo; 0.11 ± 0.01%, ~600 IU/day; 0.30 ± 0.01%, ~2,000 IU/day; and 0.65 ± 0.01%, ~4,000 IU/day group; P-trend = 0.04).ConclusionsVitamin D deficiency is associated with global hypomethylation in African Americans. Vitamin D3 supplementation increases global DNA methylation in a dose-response manner in African Americans with vitamin D deficiency.     

DNA methylation in neonates born to women receiving psychiatric care

Prenatal exposure both to maternal psychiatric illness and psychiatric medication has been linked with adverse child outcomes that affect physiological, emotional and psychiatric development. Studies suggest that epigenetic mechanisms, such as DNA methylation, may facilitate these effects. In this report, we explore the association between maternal psychiatric illness and treatment during pregnancy and neonatal DNA methylation patterns in a prospectively-characterized clinical cohort of 201 dyads. Associations between the percent of umbilical cord blood DNA methylated at 27,578 CpG sites and maternal psychiatric diagnosis, symptoms and antidepressant use were evaluated by fitting a separate linear mixed effects model for each CpG site. There were no significant changes in neonatal DNA methylation attributable to maternal psychiatric diagnosis or depressive symptoms during pregnancy. Exposure to an antidepressant medication was associated with differential methylation of CpG sites in TNFRSF21 and CHRNA2 (false discovery rate < 0.05), but the average difference in methylation for both CpG sites was less than 3% between each group. The results were not specific to type of antidepressant or duration of the exposure. This study suggests that there are no large effects of maternal psychiatric illness, depressive symptoms or prenatal exposure to antidepressants on neonatal DNA methylation. Delineation of the influence of maternal psychiatric illness and pharmacological exposures on the developing fetuses has critical implications for clinical care during pregnancy.     

Genetic Control of Individual Differences in Gene-Specific Methylation in Human Brain

We have observed extensive interindividual differences in DNA methylation of 8590 CpG sites of 6229 genes in 153 human adult cerebellum samples, enriched in CpG island “shores” and at further distances from CpG islands. To search for genetic factors that regulate this variation, we performed a genome-wide association study (GWAS) mapping of methylation quantitative trait loci (mQTLs) for the 8590 testable CpG sites. cis association refers to correlation of methylation with SNPs within 1 Mb of a CpG site. 736 CpG sites showed phenotype-wide significant cis association with 2878 SNPs (after permutation correction for all tested markers and methylation phenotypes). In trans analysis of methylation, which tests for distant regulation effects, associations of 12 CpG sites and 38 SNPs remained significant after phenotype-wide correction. To examine the functional effects of mQTLs, we analyzed 85 genes that were with genetically regulated methylation we observed and for which we had quality gene expression data. Ten genes showed SNP-methylation-expression three-way associations—the same SNP simultaneously showed significant association with both DNA methylation and gene expression, while DNA methylation was significantly correlated with gene expression. Thus, we demonstrated that DNA methylation is frequently a heritable continuous quantitatively variable trait in human brain. Unlike allele-specific methylation, genetic polymorphisms mark both cis- and trans-regulatory genetic sites at measurable distances from their CpG sites. Some of the genetically regulated DNA methylation is directly connected with genetically regulated gene expression variation.     

Genetic and non-genetic influences during pregnancy on infant global and site specific DNA methylation: role for folate gene variants and vitamin B12

Inter-individual variation in patterns of DNA methylation at birth can be explained by the influence of environmental, genetic and stochastic factors. This study investigates the genetic and non-genetic determinants of variation in DNA methylation in human infants. Given its central role in provision of methyl groups for DNA methylation, this study focuses on aspects of folate metabolism. Global (LUMA) and gene specific (IGF2, ZNT5, IGFBP3) DNA methylation were quantified in 430 infants by Pyrosequencing®. Seven polymorphisms in 6 genes (MTHFR, MTRR, FOLH1, CβS, RFC1, SHMT) involved in folate absorption and metabolism were analysed in DNA from both infants and mothers. Red blood cell folate and serum vitamin B₁₂ concentrations were measured as indices of vitamin status. Relationships between DNA methylation patterns and several covariates viz. sex, gestation length, maternal and infant red cell folate, maternal and infant serum vitamin B₁₂, maternal age, smoking and genotype were tested. Length of gestation correlated positively with IGF2 methylation (rho = 0.11, p = 0.032) and inversely with ZNT5 methylation (rho = −0.13, p = 0.017). Methylation of the IGFBP3 locus correlated inversely with infant vitamin B₁₂ concentration (rho = −0.16, p = 0.007), whilst global DNA methylation correlated inversely with maternal vitamin B₁₂ concentrations (rho = 0.18, p = 0.044). Analysis of common genetic variants in folate pathway genes highlighted several associations including infant MTRR 66G>A genotype with DNA methylation (χ² = 8.82, p = 0.003) and maternal MTHFR 677C>T genotype with IGF2 methylation (χ² = 2.77, p = 0.006). These data support the hypothesis that both environmental and genetic factors involved in one-carbon metabolism influence DNA methylation in infants. Specifically, the findings highlight the importance of vitamin B₁₂ status, infant MTRR genotype and maternal MTHFR genotype, all of which may influence the supply of methyl groups for DNA methylation. In addition, gestational length appears to be an important determinant of infant DNA methylation patterns.     

Reconstructing the Population Genetic History of the Caribbean

The Caribbean basin is home to some of the most complex interactions in recent history among previously diverged human populations. Here, we investigate the population genetic history of this region by characterizing patterns of genome-wide variation among 330 individuals from three of the Greater Antilles (Cuba, Puerto Rico, Hispaniola), two mainland (Honduras, Colombia), and three Native South American (Yukpa, Bari, and Warao) populations. We combine these data with a unique database of genomic variation in over 3,000 individuals from diverse European, African, and Native American populations. We use local ancestry inference and tract length distributions to test different demographic scenarios for the pre- and post-colonial history of the region. We develop a novel ancestry-specific PCA (ASPCA) method to reconstruct the sub-continental origin of Native American, European, and African haplotypes from admixed genomes. We find that the most likely source of the indigenous ancestry in Caribbean islanders is a Native South American component shared among inland Amazonian tribes, Central America, and the Yucatan peninsula, suggesting extensive gene flow across the Caribbean in pre-Columbian times. We find evidence of two pulses of African migration. The first pulse—which today is reflected by shorter, older ancestry tracts—consists of a genetic component more similar to coastal West African regions involved in early stages of the trans-Atlantic slave trade. The second pulse—reflected by longer, younger tracts—is more similar to present-day West-Central African populations, supporting historical records of later transatlantic deportation. Surprisingly, we also identify a Latino-specific European component that has significantly diverged from its parental Iberian source populations, presumably as a result of small European founder population size. We demonstrate that the ancestral components in admixed genomes can be traced back to distinct sub-continental source populations with far greater resolution than previously thought, even when limited pre-Columbian Caribbean haplotypes have survived.     

Relationship between tumor DNA methylation status and patient characteristics in African-American and European-American women with breast cancer

Aberrant DNA methylation is critical for development and progression of breast cancer. We investigated the association of CpG island methylation in candidate genes and clinicopathological features in 65 African-American (AA) and European-American (EA) breast cancer patients. Quantitative methylation analysis was carried out on bisulfite modified genomic DNA and sequencing (pyrosequencing) for promoter CpG islands of p16, ESR1, RASSF1A, RARβ2, CDH13, HIN1, SFRP1 genes and the LINE1 repetitive element using matched paired non-cancerous and breast tumor specimen (32 AA and 33 EA women). Five of the genes, all known tumor suppressor genes (RASSF1A, RARβ2, CDH13, HIN1 and SFRP1), were found to be frequently hypermethylated in breast tumor tissues but not in the adjacent non-cancerous tissues. Significant differences in the CDH13 methylation status were observed by comparing DNA methylation between AA and EA patients, with more obvious CDH13 methylation differences between the two patient groups in the ER- disease and among young patients (age<50). In addition, we observed associations between CDH13, SFRP1, and RASSF1A methylation and breast cancer subtypes and between SFRP1 methylation and patient's age. Furthermore, tumors that received neoadjuvant therapy tended to have reduced RASSF1A methylation when compared with chemotherapy na?ve tumors. Finally, Kaplan Meier survival analysis showed a significant association between methylation at 3 loci (RASSF1A, RARβ2 and CDH13) and reduced overall disease survival. In conclusion, the DNA methylation status of breast tumors was found to be significantly associated with clinicopathological features and race/ethnicity of the patients.     

The roles of DNA methylation of NR3C1 and 11β-HSD2 and exposure to maternal mood disorder in utero on newborn neurobehavior

Exposure to maternal mood disorder in utero may program infant neurobehavior via DNA methylation of the glucocorticoid receptor (NR3C1) and 11β-hydroxysteroid dehydrogenase type 2 (11β-HSD-2), two placental genes that have been implicated in perturbations of the hypothalamic pituitary adrenocortical (HPA) axis. We tested the relations among prenatal exposure to maternal depression or anxiety, methylation of exon 1F of NR3C1 and 11β-HSD-2, and newborn neurobehavior. Controlling for relevant covariates, infants whose mothers reported depression during pregnancy and showed greater methylation of placental NR3C1 CpG2 had poorer self-regulation, more hypotonia, and more lethargy than infants whose mothers did not report depression. On the other hand, infants whose mothers reported anxiety during pregnancy and showed greater methylation of placental 11β-HSD-2 CpG4 were more hypotonic compared with infants of mothers who did not report anxiety during pregnancy. Our results support the fetal programming hypothesis and suggest that fetal adjustments to cues from the intrauterine environment, in this case an environment that could be characterized by increased exposure to maternal cortisol, may lead to poor neurodevelopmental outcomes.     

Gene-Based Sequencing Identifies Lipid-Influencing Variants with Ethnicity-Specific Effects in African Americans

Although a considerable proportion of serum lipids loci identified in European ancestry individuals (EA) replicate in African Americans (AA), interethnic differences in the distribution of serum lipids suggest that some genetic determinants differ by ethnicity. We conducted a comprehensive evaluation of five lipid candidate genes to identify variants with ethnicity-specific effects. We sequenced ABCA1, LCAT, LPL, PON1, and SERPINE1 in 48 AA individuals with extreme serum lipid concentrations (high HDLC/low TG or low HDLC/high TG). Identified variants were genotyped in the full population-based sample of AA (n = 1694) and tested for an association with serum lipids. rs328 (LPL) and correlated variants were associated with higher HDLC and lower TG. Interestingly, a stronger effect was observed on a “European” vs. “African” genetic background at this locus. To investigate this effect, we evaluated the region among West Africans (WA). For TG, the effect size among WA was the same in AA with only African local ancestry (2–3% lower TG), while the larger association among AA with local European ancestry matched previous reports in EA (10%). For HDLC, there was no association with rs328 in AA with only African local ancestry or in WA, while the association among AA with European local ancestry was much greater than what has been observed for EA (15 vs. ∼5 mg/dl), suggesting an interaction with an environmental or genetic factor that differs by ethnicity. Beyond this ancestry effect, the importance of African ancestry-focused, sequence-based work was also highlighted by serum lipid associations of variants that were in higher frequency (or present only) among those of African ancestry. By beginning our study with the sequence variation present in AA individuals, investigating local ancestry effects, and seeking replication in WA, we were able to comprehensively evaluate the role of a set of candidate genes in serum lipids in AA.     

Maternal genome-wide DNA methylation patterns and congenital heart defects

The majority of congenital heart defects (CHDs) are thought to result from the interaction between multiple genetic, epigenetic, environmental, and lifestyle factors. Epigenetic mechanisms are attractive targets in the study of complex diseases because they may be altered by environmental factors and dietary interventions. We conducted a population based, case-control study of genome-wide maternal DNA methylation to determine if alterations in gene-specific methylation were associated with CHDs. Using the Illumina Infinium Human Methylation27 BeadChip, we assessed maternal gene-specific methylation in over 27,000 CpG sites from DNA isolated from peripheral blood lymphocytes. Our study sample included 180 mothers with non-syndromic CHD-affected pregnancies (cases) and 187 mothers with unaffected pregnancies (controls). Using a multi-factorial statistical model, we observed differential methylation between cases and controls at multiple CpG sites, although no CpG site reached the most stringent level of genome-wide statistical significance. The majority of differentially methylated CpG sites were hypermethylated in cases and located within CpG islands. Gene Set Enrichment Analysis (GSEA) revealed that the genes of interest were enriched in multiple biological processes involved in fetal development. Associations with canonical pathways previously shown to be involved in fetal organogenesis were also observed. We present preliminary evidence that alterations in maternal DNA methylation may be associated with CHDs. Our results suggest that further studies involving maternal epigenetic patterns and CHDs are warranted. Multiple candidate processes and pathways for future study have been identified.     

Folic Acid Supplementation during Pregnancy Induces Sex-Specific Changes in Methylation and Expression of Placental 11β-Hydroxysteroid Dehydrogenase 2 in Rats

In the placenta, 11β-hydroxysteroid dehydrogenase type 2 (11β-HSD2) limits fetal glucocorticoid exposure and its inhibition has been associated to low birth weight. Its expression, encoded by the HSD11B2 gene is regulated by DNA methylation. We hypothesized that maternal diets supplemented with folic acid (FA) during pregnancy modify the expression of placental HSD11B2 through gene methylation. Wistar rats were fed with high (8 mg/kg) or normal low (1mg/kg, control) levels of FA during pregnancy. Concentrations of mRNA and protein in placentas were determined by qRT-PCR and Western blot respectively. Methylation in five CpG sites of the placental HSD11B2 promoter (−378 to −275) was analyzed by bacterial cloning and subsequent sequencing. In the FA-supplemented group, mRNA and protein levels of 11β-HSD2 decreased by 58% and increased by 89%, respectively, only in placentas attached to males. In controls, most CpG sites were not methylated except for the CpG2 site which was 80% methylated. CpG2 methylation level increased under the FA treatment; however, only in placentas attached to females was this increase significant (113%). This change was not related to HSD11B2 expression. Fetal weight of females from FA- supplemented mothers was 6% higher than females from control mothers. In conclusion, this is the first study reporting that FA over supplementation during pregnancy modifies the placental HSD11B2 gene expression and methylation in a sex-dependent manner, suggesting that maternal diets with high content of FA can induce early sex-specific responses, which may lead to long-term consequences for the offspring.     

Postnatal stability,tissue, and time specific effects of AHRR methylation change in response to maternal smoking in pregnancy

The intrauterine environment has the potential to “program” the developing fetus in a way that can be potentially deleterious to later health. While in utero environmental/stochastic factors are known to influence DNA methylation profile at birth, it has been difficult to assign specific examples of epigenetic variation to specific environmental exposures. Recently, several studies have linked exposure to smoking with DNA methylation change in the aryl hydrocarbon receptor repressor (AHRR) gene in blood. This includes hypomethylation of AHRR in neonatal blood in response to maternal smoking in pregnancy. The role of AHRR as a negative regulator of pathways involved in pleiotropic responses to environmental contaminants raises the possibility that smoking-induced hypomethylation is an adaptive response to an adverse in utero environmental exposure. However, the tissue specificity of the response to maternal smoking, and the stability of the methylation changes early in life remain to be determined. In this study we analyzed AHRR methylation in three cell types—cord blood mononuclear cells (CBMCs), buccal epithelium, and placenta tissue—from newborn twins of mothers who smoked throughout pregnancy and matched controls. Further, we explored the postnatal stability of this change at 18 months. Our results confirm the previous association between maternal smoking and AHRR methylation in neonatal blood. In addition, this study expands the region of AHRR methylation altered in response to maternal smoking during pregnancy and reveals the tissue-specific nature of epigenetic responses to environmental exposures in utero. Further, the evidence for postnatal stability of smoking-induced epigenetic change supports a role for epigenetics as a mediator of long-term effects of specific in utero exposures in humans. Longitudinal analysis of further specific exposures in larger cohorts is required to examine the extent of this phenomenon in humans.     

Associations between Vascular Health Indices and Serum Total,Free and Bioavailable 25-Hydroxyvitamin D in Adolescents

Objective

The role of vitamin D in cardiovascular health remains debated as results have been inconsistent. Previous studies have not considered the bioavailability of 25-hydroxy vitamin D [25(OH)D]. Objectives of our study were to investigate the association between serum concentrations of total, free and bioavailable 25(OH)D and independent predictors of cardiovascular risk such as flow mediated dilatation (FMD) and augmentation index (AIx).

Design

This cross-sectional study included 47 post-menarchal, adolescent females [31 African American (AA) and 16 European American (EA)].

Methods

AIx was standardized to a heart rate of 75 beats/min (AIx75). Free and bioavailable 25(OH)D concentrations were calculated from standard formulas.

Results and Conclusions

Mean age of the participants was 15.8±1.4 years and mean body mass index was 23.1±4.0 kg/m². Serum total 25(OH)D was not associated with FMD, but was positively associated with AIx75 in the adjusted model (rho = 0.4, P = 0.03). AIx75 was positively associated with bioavailable 25(OH)D (rho = 0.4, P = 0.004) and free 25(OH)D (rho = 0.4, P = 0.009) and the associations persisted after adjusting for covariates. In race-specific analyses, total, free and bioavailable 25(OH)D were strongly positively associated with AIx75 in AA (rho = 0.5, 0.4, 0.4, respectively), which persisted even after adjusting for covariates. Whereas in EA there was an inverse association between total 25(OH)D and AIx75 in EA (rho = −0.6), which attenuated after adjusting for covariates.

Conclusion

Circulating total, free and bioavailable 25(OH)D were associated with arterial stiffness in adolescent girls, and these associations were race dependent. Notwithstanding, the implications of associations between vascular function indices and 25(OH)D remains unclear.