Race/ethnicity differences in vitamin D levels and impact on cardiovascular disease, bone health, and oral health

Cardiovascular disease: Abdominal aortic calcification

Vitamin D is linked to CVD risk. But the results of studies looking for a link between vitamin D and markers of CVD and heart disease have been nuanced and mixed. A nationwide, randomized control study of vitamin D3 supplementation of 2000 IU per day, the VITAL study, showed no improvement in cardiovascular risk, with myocardial infarctions (MCIs), strokes, and death from CVD staying statistically the same across control and experimental groups even after a five-year follow-up. The results, or the lack of them, did not vary by race/ethnicity (Manson et al., 2019). NHANES-III data, collected from 1988 to 1994, showed an inverse association between 25(OH)D levels compared between highest and lowest quartiles and death from coronary heart disease (CHD), with no effect when all quartiles were considered. The effect remained after controls were added for sex, diabetes, income, age, BMI, blood pressure (BP), dietary elements, and similar lifestyle, socioeconomic, and physiological parameters (Daraghmeh et al., 2016).

A correlation between vitamin D levels and CVD may well indicate that low vitamin D is a marker of CVD than its cause. Vitamin D levels are affected by VDR gene polymorphisms, sex hormones and hence sex, age, use of corticosteroids, obesity, with fat sequestering the lipophilic vitamin, and kidney disease (Marquina et al., 2018). Thus, variations in the CYP24A1 and VDR genes have been associated with coronary artery calcification (CAC) (Shen et al., 2010); low 25(OH)D was correlated to heart failure rates in white but not AA adults in a study that spanned 20 years (Lutsey et al., 2015). Complementing that, Mendelian randomization studies have frequently failed to find a correlation between various SNPs in CYP genes and coronary artery disease, indicating vitamin D may be a biomarker instead of a causal element (Marquina et al., 2018).

Balanced against a null result for a relation between 25(OH)D and heart failure in AAs, the Jackson Heart Study, of 5000 AA adults in Jackson, Mississippi, showed a direct correlation of 25(OH)D and HDL cholesterol, no correlation with LDL, and inverse correlations with CVD risk factors such as BMI, waist circumference, BP, and fasting plasma glucose, all mediated partly through C-reactive protein and aldosterone. This is a pathway through which 25(OH)D can lower BP by downregulating the renin-angiotensin-aldosterone (RAAS) system by inhibiting renin synthesis, the start of the RAAS pathway. Lowering aldosterone reduces how much water and sodium the body retains, lowering blood pressure. The results of the path analysis held after controlling for age, sex, education, smoking, BMI, season, and other behavioral, socioeconomic and physiological elements (Khan et al., 2016). Large parts of these correlations were replicated in the overall population, including with no race/ethnicity control, with three cycles of NHANES data from 2001 to 2006. That study did not perform path analysis and so could not detect mediating influences (Ganji et al., 2020). Still another study with NHANES data from 2003 to 2006 found an inverse correlation between treatment-resistant hypertension, a major CVD risk, and 25(OH)D levels in the overall population (Alagacone et al., 2020). A third study found an inverse correlation between 25(OH)D levels and CVD mortality for the 1988 to 1994 NHANES sample, with participant mortality followed till 2000 (Melamed et al., 2008). Recent reviews have concluded that a J-shaped or U-shaped inverse correlation exists between CVD elements and 25(OH)D, that is, in some range there is a protective effect (Durup et al., 2015), but this relationship may well indicate unmodeled elements than any true biological effect (Grant et al., 2016). Vitamin D also regulates parathyroid hormone (PTH) as part of the calcium homeostasis loop involving renal tubules, the thyroid gland, and the parathyroid. Low levels of 25(OH)D lead to high levels of PTH; this can cause myocardial hypertrophy and secretion of more aldosterone, all raising CVD risk (Marquina et al., 2018).

The effects of vitamin D are, however, complex. As noted before, randomized control trials have largely failed to show CVD benefits from vitamin D supplementation (Zitterman, 2018). A recent meta analysis of vitamin D supplementation and CVD risk across 21 randomized clinical trials found that supplementation had no benefits with regards to cardiovascular events and related deaths, strokes, MCI, or even overall mortality (Barbarawi et al., 2019). A case-cohort study of 79,705 postmenopausal females in the multi-ethnic Women’s Health Initiative Observational Study (WHI-OS) found no correlation between CVD risk and total 25(OH)D, free 25(OH)D, or bioavailable 25(OH)D in either white or AA females. The researchers found no PTH link to CVD risk in AAs (though they did find one in white females) (Zhang et al., 2019). The 6436-participant Multi Ethnic Study of Atherosclerosis (MESA), which weighted ethnic/racial minorities up in its sample, found an inverse link between 25(OH)D and CHD event applying to white adults but not AAs. The authors speculate that the cause may be an interruption of the RAAS sequence that leads to more aldosterone with lower 25(OH)D. AAs have higher circulating 1,25-dihydroxyvitamin D, likely because of increased activity by the hydroxylase coded for by the CYP27B1 gene, polymorphic partly by ethnicity. On the catabolism side, the hydroxylase responsible has reduced activity in AAs. VDR binding strength too may differ by race/ethnicity, with some genes here polymorphic by race/ethnicity. In addition, AAs have statistically higher rates of PTH, added to their lower levels of 25(OH)D (Robinson-Cohen et al., 2013). Nevertheless, they have lower levels of arterial calcification, particularly CAC, one CVD risk element. This is despite their having higher scores on CVD risk factors such as BP, diabetes, and smoking rates (see Table 1 of Orimoloye et al. (2018)).

For measuring CVD risk, CAC has proved a durable indicator. The standard measure for CVD risk, the Framingham Risk Score (FRS), which incorporates sex, age, BP, total cholesterol, LDL, HDL, smoking status, and diabetes status, predicts CVD moderately well for white and AA subjects (D’Agostino et al., 2001). The coronary calcium score, obtained with CT scans, provides even more predictive information for CVD events and CHD for both white Americans and AAs, both males and females. This despite the fact that the incidence and the score itself differ statistically between race/ethnicities, but within a race/ethnicity it still provides a good CVD-prediction marker when added to the FRS (Detrano et al., 2008). It also allows the researcher/clinician to steer clear of self-report data, with its biases and inaccuracies. Another study with the MESA cohort has concluded that what matters for CVD risk is the absolute CAC score within a group and not relative amounts for the individual compared to the group as a whole (Budoff et al., 2009). Calcified atherosclerotic plaque, in general, has been linked to higher 25(OH)D and 1,25(OH)D levels in AAs, with overall levels of atherosclerotic plaque much lower in them than in white population samples (Freedman and Register, 2012). A study of a MESA subsample showed that AAs with greater share of European ancestry had higher prevalence and amounts of CAC (though with lower scores on carotid intima-media thickness, another marker of CVD), indicating a genetic link (Wassel et al., 2009).

Arterial calcification in the abdominal aorta (AAC) works similar to CAC in predicting cardiovascular risk. AAC predicts both CVD risk and fracture risk, providing a mechanistic explanation for a noted link between CVD and bone fragility. The link is particularly strong between hip and vertebral fractures and MCI, stroke, and heart failure (Szulc, 2016). AAC has been linked to 25-OH(D) deficiency, albeit mediated by age, in a sample of 769 Asian females, 75% of whom were menopausal. This study confirmed an inverse association between AAC and bone loss (Kim et al., 2012). Another study in Rabat, Morocco failed to find a link between vitamin D status and AAC in 429 postmenopausal females. However, this study used threshold values to categorize participants as hypovitaminosis, insufficient, and sufficient categories, with 10 ng/mL and 20 ng/mL as thresholds, and then used chi-squared tests and ANOVA to check the means (El Maghraoui et al., 2018). In AA specifically, vitamin D supplementation was found to have no influence on AAC in females over 60, in the PODA trial (Brahmbhatt et al., 2020). However, in a study in Tartu, Estonia, of 152 males over the age of 60, 25(OH)D levels were found to be correlated to aortic calcification, with the entire aorta (thoracic, abdominal, and pelvic) checked up to the bifurcation into the iliac arteries (Zagura et al., 2011).

The literature and the data thus lead to a hypothesis filling a gap testable with NHANES 2011–14 data. Hypothesis 1: AAC is directly linked to 25(OH)D levels in white Americans but unlinked in AAs.

Bone health: Bone mineral density and fracture risks

Bone fracture risk has many components. A risk score, the FRAX score, based on cohort studies in multiple countries and ethnicities, was developed by the University of Sheffield, UK, in 2008, based on data collected in collaboration with the World Health Organization. This tool uses an algorithm to weight risk factors using t-scores and z-scores, that is, deviations from an average, to come up with a risk number. For the US, the risk elements for people not taking any medication include age, sex, BMI, individual and family history of fracture, smoking status, three or more alcoholic drinks a day, rheumatoid arthritis, and areal BMD (aBMD). The tool weights the risk differently for races/ethnicities (Center for Metabolic Bone Diseases, 2008). FRAX predicts hip fractures well but not necessarily other types of fractures, because some such as spine fractures do not get even reported or diagnosed. Biochemical markers such as 25OH(D) are not included; all the measures included are detected noninvasively (Unnanuntana et al., 2010).

The effect of vitamin D on bone is from its role in calcium balance. When serum CA²⁺ levels are low, the parathyroid gland detects this and increases production of PTH. PTH, in turn, stimulates renal CYP27B1, which hydroxylates 25(OH)D to 1,25(OH)D, which, in turn, inhibits both CYP27B1 and PTH in a negative feedback loop. In another feedback loop, 1,25(OH)D stimulates FGF23 in bone, which too inhibits CYP27B1. The catabolism of 1,25(OH)D, which usually exists for just a few hours, is triggered by the dehydroxylase CYP24A1. This dehydroxylase exists in all cells with VDRs (a nuclear receptor), and is stimulated by the presence of 1,25(OH)D. In the kidneys alone, CYP24A1 is inhibited by low CA²⁺ and PTH, that is, inversely with CYP27B1, the anabolic hydroxylase (Christakos et al., 2020).

Vitamin D facilitates CA²⁺ absorption in the intestine and, along with PTH, in the distal renal tubules and collecting tubules. Absorption is via the CA²⁺ channel TRPV5 and calbindin and excretion via CA²⁺ ATPases and Na⁺/CA²⁺ antiporters. VDRs are also expressed in osteoblasts and osteocytes but to a lesser extent in osteoclasts and chondrocytes. When serum CA²⁺ levels fall, both PTH and 1,25(OH)D levels rise, and this leads to bone resorption to replenish CA²⁺ levels. The resorption is mediated by expression of RANKL, which stimulates osteoclast differentiation and transfer of CA²⁺ from cortical and trabecular bone to serum. In parallel, bone remineralization is inhibited by the rise in pyrophosphate levels induced by 1,25(OH)D (Christakos et al., 2020). Most 25(OH)D in the blood is bound to VDBP; it is the free 25(OH)D that binds to VDRs and affects bone remodeling (LeBoff et al., 2020).

Bone mineral density, specifically aBMD, the area-wise density measured by DEXA, is strongly linked to fracture risk. Some studies indicate volumetric BMD too is linked to fracture risk, but this element is harder to measure, requiring quantitative computerized tomography (qCT). A recent meta-analysis of aBMD and fractures in studies looking at patients undergoing osteoporosis treatments showed a strong association of hip, femoral-neck, and lumbar-spine BMD with vertebral fractures and a less strong but significant association with hip fractures. For nonvertebral fractures, the heterogeneity of the definitions of such fractures and the occurrence of trauma fractures made the analysis not useful (Bouxsein et al., 2019). BMD itself depends, among other elements, on 25(OH)D levels: An analysis of NHANES III data found BMD to rise with 25(OH)D in all races/ethnicities (Bischoff-Ferrari et al., 2004). Even though BMD and Vitamin-D are linked, vitamin D supplementation, has had no effect on either aBMD measured with DEXA scans or vBMD and cortical and trabecular microstructure measured with peripheral quantitative CT scans. These null results, from the VITAL study among others, have held across age, race/ethnicity, and sex (LeBoff et al., 2020). Even higher doses of vitamin D, 4000 IU/day, led to reduced vBMD in the radius bone, and still higher doses in the tibia as well. This, however, must be considered more an instance of vitamin D poisoning from overdose (Burt et al., 2019).

Vitamin D interacts with race/ethnicity in its effect on BMD. AAs have higher BMD and lower risk of osteoporosis than white Americans though AAs have statistically lower levels of 25(OH)D. Though environmental and physical elements such as higher levels of obesity play a role in the lower 25(OH)D levels in AAs, as adipose tissue can sequester vitamin D or volumetrically dilute it, a large part is considered to be genetic. AAs also have statistically lower consumption of calcium, both from social elements and from the absence of the lactase enzyme making it harder to consume lactose-containing milk (Freedman and Register, 2012). Despite that, BMD as measured by DEXA, which measures largely aBMD, and vBMD measured by qCT, are higher in AAs than in others (Bachrach et al., 2019). QCT has shown higher cancellous, that is, trabecular, bone density in the spine in AA children than in white children; at the femoral shaft, AA children show greater bone area but no difference in cortical thickness. The active biochemical vitamin D marker measured in that study, 1,25(OH)D, showed a sex but no race/ethnicity difference and did not correlate with any bone measure (Gilsanz et al., 1998). Both larger bone mass and across-the-board volumetric density (vBMD) act as protective elements in AAs, lowering their osteoporosis and fracture risk.

The NHANES-III association between vitamin D and BMD was weaker than for the rest in AAs and found only among older adults. Its statistical significance vanished when physical activity was added to the list of confounders (Bischoff-Ferrari et al., 2004). This led the researchers to label BMD as the better endpoint marker for measuring bone health against 25(OH)D for non-AA populations, leaving the issue open for AA, especially older AA, populations (Bischoff-Ferrari et al., 2006). The higher BMD and lower 25(OH)D for AAs and the concomitant lower fracture risk does not bode well for a single 25(OH)D threshold for bone health across races/ethnicities, independent of whether the differences are environmental or inherited.

Gutiérrez et al. (2011) analyzed the relationship between vitamin D, PTH levels, and aBMD in NHANES data from 2003 to 2006. They found that BMD and 25(OH)D were directly related in non-AAs but not in AAs; PTH was inversely related to 25(OH)D in all groups but was maximally suppressed at a 25(OH)D level of 26 ng/mL in AAs. No measures of fracture risk were checked against either BMD or 25(OH)D.

These open issues and the availability of unanalyzed aBMD and FRAX data in NHANES 2011–12 and 2013–14 lead to the next hypothesis: Levels of 25(OH)D are directly, that is, positively, related to spine and femur BMD and spine and hip FRAX in white populations, but no relationship exists in AAs in general.

Oral health: Periodontitis vulnerability

Periodontitis (PD) arises from loss of gum and supporting tissues for teeth, eventually leading to tooth loss. It is preceded by gingivitis—inflammation of the gum—and tissue loss follows from immune response to the oral microbiome, that is, the body’s response to microbial infection of the gums (Machado et al., 2020). PD is detected by probing of the pockets in gums and then adjusting that against recession of the gingiva to obtain a measure of attachment loss (AL), loss of tooth-supporting tissue. Thresholds are applied for a diagnosis of PD. The etiology of PD is multifactorial, with strong environmental and some genetic components. One strong environmental component is smoking. Studies have estimated 52% of PD in the US as due to current or former smoking, the 52% figure being an attributable fraction, that is, what proportion will not have PD if no one smoked. Almost 3/4^ths of PD risk in current smokers is from smoking, again considered as an attributable fraction (Tomar and Asma, 2000).

Meta studies have shown higher vitamin D reduces periodontitis, though the issue is still open. Vitamin D supplementation has unclear effect, again across studies. Vitamin D acts probably by binding to VDRs in gingival fibroblasts and periodontal ligament cells, inducing cathelicidin, an antimicrobial protein. This immune response is protective against PD (Machado et al., 2020). The association of cathelicidin and another antimicrobial, human beta defensin 2, with vitamin D is believed to be mediated by the binding of 1,25(OH)D to the VDR that then binds to the nuclear DNA sequence called the vitamin D response element (Wang et al., 2004). Vitamin D, however, has complex interactions with both the humoral and the cell-based immune systems, and these anti-inflammatory actions too might influence its effects on periodontitis (Dragonas et al., 2020).

An analysis of NHANES data from 1999 to 2004 estimated that higher levels of Vitamin D decreased the odds of periodontitis by 39% and 24% in former and non-smoker groups, respectively. Several demographic and environmental elements mediated the influence in the negative direction, that is, increasing the risk: AA race/ethnicity, older age, smoking status, male sex, and higher blood levels of retinyl stearate, retinyl palmitate, lead, non-coplanar polychlorinated biphenyls (PCBs) PCB105, PCB146, PCB172, PCB177, PCB178, PCB183, and PCB206, and higher urine levels of antimony. Higher levels of serum carotenoids and vitamin D had a protective influence (Emecen-Huja et al., 2019). Another analysis of NHANES III showed vitamin D and periodontitis correlated in both males and females over the age of 50, even after controlling for BMD, indicating mediation via the immunobeneficent role of vitamin D. This version of NHANES had just a half-mouth exam (random upper and lower quadrant) and only two sites per tooth. The researchers used vitamin D quintiles as the independent variable (IV) and raw mean AL as the dependent variable (DV), with a host of control variables (Dietrich et al., 2004).

Opposing this, the Dragonas et al. (2020) meta study mentions several studies showing a lack of effect for both 25(OH)D and 1,25(OH)D on periodontitis markers, as opposed to raw AL. However, this metastudy found a lack of result in studies largely from non-US populations, specifically from Finland; Korea; Pomerania, Germany; and Beijing, China. The one US study with a null result was on postmenopausal women. In some cases, the results were based on a linear regression against mean AL taken as a continuous variable. However, mean AL is not normally distributed and neither would be the measurement error, violating a requirement for ordinary least-squares linear regression. The environmental influences on periodontitis, and on vitamin D levels as well, are so significant that the heterogeneity in the meta review makes conclusions difficult. Another recent review by Millen and Pavlesen (2020) also found some support for a vitamin D role in periodontitis. This survey focused on cross-sectional studies, noting the paucity of longitudinal studies on the subject. Cross-sectional studies cannot establish the temporal direction of the connection; there are some indications gingivitis induces gum tissue to produce more 25(OH)D, which confounds both the overall expected result and the directionality of the influence.

Race/ethnicity differences in the occurrence of periodontitis are known. Adult AAs develop periodontitis at higher rates than white Americans do, with a large part of the causation likely environmental. An analysis of the 2009–10 NHANES data, which conducted full-mouth 6-site periodontal exams, showed higher rates, though not the highest rates, of periodontitis in AAs. Age, sex, smoking status, and socioeconomic-status elements such as education and income to federal poverty level (FPL) ratio were all associated with AL: older, smoking, less educated, and poorer males had higher AL (Eke et al., 2012). This study did not look at the influence of vitamin D on periodontitis; the earlier analysis of NHANES III by Dietrich et al. (2004) did look at that and found no differences across races/ethnicities, that is, stratification by race/ethnicity yielded results similar to that of the total. However, one notes that the half-mouth 3-site periodontal exam of NHANES III significantly underestimated overall periodontitis rates, assuming rates did not dramatically jump from 1994 to 2009 (Eke et al., 2012).

This leads to our final hypothesis: Mean attachment loss is linked to vitamin D in both white and AA adult males and females, after controlling for age, sex, smoking status, education level, and family income.

DVs and models

For hypothesis 1, the DV is AAC24 in file DXXAAC, available for 2013–14. It has been found to be a reliable predictor for abdominal aortic calcification more accurately measured through other means, that is, via qCT scans. Note that AAC may contain both intimal and medial calcification, and medial calcification is not believed to be related to CVD risk. An advantage of AAC over qCT measures is that AAC is readily available as DEXA scans of the spine are recommended by several medical organizations to detect vertebral fractures in those older than 65 or 70 (DeSapri and Brook, 2020). AAC information is a byproduct of lateral scans of the spine. The variable itself is discrete with 24 values. It can be considered as count data, and its distribution does show a Poisson shape (see Fig. 1). The mean and variance are not equal, and hence a negative binomial distribution is a good model. The data have zero inflation, that is, bunching at the 0s beyond what the probability distribution would predict; hence the right GLM model to use is the zero-inflated negative binomial. Figure 2 shows the shape of the data after zero removal, a Poisson/negative binomial shape. For svyVGAM package, that would be the svy vglm function with GLM function as zinegbinomialff. The design object for this call comes from the survey package’s svydesign() call with id being SDMVPSU, and strata being SDMVSTRA. The setting of weights, required in this function, is explained later.

For hypothesis 2, there exist four DVs: spine BMD DXXOSBMD in file DXXSPN, femur BMD DXXOFBMD in file DXXFEM, and hip and osteoporotic FRAX calculated from files DXXFRX. All are available for 2013–14 though only for those above the age of 40. Both spine and femur BMDs are roughly normal shaped (Fig. 3, Fig. 4) and the QQ plots are roughly lines (not shown), and so direct linear regression (svyglm call from survey) is used. This analysis eschews K-S and Shapiro Wilk tests for normality as these are mostly ineffective with large sample sizes. For the BMD variables, only aBMD information is available as this is a DEXA scan; vBMD and cortical and trabecular bone geometry require qCT scans. Hip FRAX is neither normal nor log-normal, so it was used as is (Fig. 5). Osteoporotic FRAX is Poisson-shaped but with unequal mean and variance and with no zero inflation (Fig. 6); hence survey vglm with negbinomial as the GLM is used. The hip FRAX is picked up from either DXXFRAX1 or DXXFRAX3 depending on the value in DXXPRVFX (prior fracture present or not). Osteoporotic FRAX is calculated similarly except that DXXFRAX2 or DXXFRAX4 is used.

• Download figure
• Open in new tab

Figure 3.

Histogram of spine-BMD raw scores

• Download figure
• Open in new tab

Figure 4.

Histogram of femur-BMD raw scores

• Download figure
• Open in new tab

Figure 5.

Histogram of hip FRAX raw scores

• Download figure
• Open in new tab

Figure 6.

Histogram of osteoporotic FRAX raw scores

For hypothesis 3, clinical AL is the DV. This variable is available for 2011–14, and is calculated from raw data on probing depth and gingival recession by direct subtraction, with recession coded as negative when needed to get the math straight (Tomar and Asma, 2000). PD is diagnosed based on threshold criteria for both AL and pocket depth, but just AL may be a better measure of PD risk as opposed to prevalence. In NHANES 2009–10 periodontitis severity held steady with age, but AL severity rose. As risk is known to rise with age, AL with no thresholding seems the better measure for non-diagnostic purposes (Eke et al., 2012). Total AL is problematic because not all examined subjects have values for all sites. Hence, we must use the mean. Next, the variable is rounded to the nearest 0.1 to render it discrete. It now displays a Poisson shape (Fig. 7). There is no zero inflation (zero pockets for the whole mouth are rare), mean does not equal variance, and so the right regression call is svy vglm with negbinomial.

• Download figure
• Open in new tab

Figure 7.

Histogram of mean attachment loss measures

IVs

The IVs are all available for the years 2011–14. However, which years are used for a specific analysis depends on the DV. Clinical AL, associated with periodontitis, is the only DV spanning the whole range; the rest cover 2013–14.

Weighting, point estimates, standard errors

NHANES has a stratified, clustered, 4-stage design: 1. primary sampling unit (PSU) sampling counties, 2. census block segments, 3. dwelling units, and 4. individuals. Data collection is every year but release every two years. The sample is meant to represent both the noninstitutionalized, resident, civilian population of the 50 states of the US plus Washington, D.C., and its demographic breakdowns, for our purposes, the non-Hispanic white and AA populations and sex and age stratifications. NHANES used the Census Bureau’s 2011 and 2013 one-year American Community Survey (ACS) demographic numbers as the reference (Johnson et al., 2014). To ensure that representativeness holds in the analysis, each sample must be weighted by a number provided in the data. The variable that holds the weights varies. For interview data, the variable is WTINT2YR; for lab and physical exam data, it is WTMEC2YR and for fasting data, in our case plasma glucose and LDL levels, it is WTSAF2YR. In any analysis, the weights to be used are the most restrictive of all the variables in the analysis, with the most restrictive corresponding to the least frequent. Thus, WTSAF2YR is the most restrictive and WTMEC2YR next. We do not use WTINT2YR as no analysis has just interview data (Chen et al., 2018).

Point estimates, such as mean or variance of data or regression coefficients, must be adjusted to reflect the survey design. This is achieved by using the variables SDMVPSU and SDMVSTRA in the DEMO file along with the weights (Chen et al., 2018). SDMVPSU may not be unique across strata. Sampling is assumed to be with replacement. The R variable survey.lonely.psu is set to “adjust,” a conservative correction that centers a stratum with just one PSU at the grand mean instead of the stratum mean. If this were a certainty PSU the variance should be 0; the correction overestimates the variance here and is hence conservative. The variance estimator for the regression coefficients is the default provided by the R svyVGAM package, which is a Horvitz-Thompson-type sandwich estimator. This ensures the standard error of the coefficient is largely independent of the constant-variance assumption for error at different values, and that the stratified design is accounted for.

Hypothesis 1

Hypothesis 1 is confirmed by the data, though with several niceties. AAC24 does depend on vit. D, with more vit. D raising AAC scores in white Americans, with AAs showing no effect. However, the lack of an effect in AAs is due to higher standard errors for the coefficients, presumably from the smaller sample size. Absence of proof in this case does not point to proof of likely absence. In the overall populace, the race/ethnicity effect vanishes when vit. D is added to the regression. This means lower vit. D is protective when it comes to AAC; the race/ethnicity difference is due to lower vit. D levels in AAs. Incidentally, the inverse effects for BMI and LDL on AAC24 in the overall data are hard to explain. The data also make it appear that AAC24 scores are lower in summer; this is probably due to the geographic variation in where NHANES samples data in summer and winter. As the full model with all our IVs thrown in shows no race/ethnicity influence on AAC24, we must conclude that lower vit. D, in and of itself, is protective for AAC. Further studies are required to see whether this result applies to CAC as well. A mechanistic interpretation is that higher 25(OH)D facilitates higher calcium absorption in the intestines and deposition of that calcium elsewhere. However, this finding contrasts with others that found low 25(OH)D to produce more CAC in white Americans.

Hypothesis 2

The second hypothesis was only partly validated. In sum, higher levels of 25(OH)D reduced FBMD in both groups, increased hip FRAX in white Americans only, and decreased osteoporotic FRAX in white Americans. The effects of 25(OH)D on osteoporotic FRAX was similar in white Americans and AAs after stratification with no control variables: in both cases higher vit. D led to higher risk of osteoporotic fractures, possibly through reducing BMD. Addition of all the control variables led to an inversion of the relationship in white Americans, with higher vit. D leading to lower osteoporosis risk. In AAs, statistical significance was lost, presumably due to the smaller sample size. The results for osteoporotic FRAX held in white Americans after adding controls for systolic BP, education, BMI, smoking, fasting plasma glucose levels, income to FPL ratio, age, and sex, along with LDL and season, both nonsignificant. The increased hip FRAX with vit. D for white Americans is likely an artifact; adding all the control variables made the effect vanish (result not shown). To recapitulate, vitamin D increase reduces BMD in both groups but reduces FRAX just in white Americans. In AAs, with no controls and hence a larger sample, 25(OH)D rise increases FRAX; this effect vanishes when controls are added. This may be because the controls do explain the 25(OH)D link or because the sample size gets reduced.

Hypothesis 3

Hypothesis 3 was borne out unequivocally. Increased vit. D levels reduce periodontal risk in white Americans, AAs, and the total sample. The results are robust after controlling for systolic BP, LDL, education, season, BMI, cigarette smoking, plasma glucose level, income to FPL ratio, age, and sex. All those variables except for BP, LDL, and season show relationships with AL in white Americans; in AAs, the above plus fasting plasma glucose and income/FPL show no relationship, all else does. BMI is a surprise, with higher BMI leading to lower AL scores. Males seemingly are more at risk than females, p < .001 in both groups; this coefficient is in fact the largest in both stratified models and the combined model (as units differ, coefficients are not directly comparable, but as the unit distance is the largest, 0 to 1, for sex, this high coefficient is all the more remarkable). In the combined model, AAs are significantly at more risk after all the additional variables are controlled for. This may indicate residual risk from ancestral genetics but more so from unmodeled environmental elements such as stress and unequal receipt of treatment even after controlling for socio-economic status.

Differential impact of vit. D

For both aortic calcification, at least in the abdominal aorta, and for bone health, measured by FRAX scores and BMDs, the influence of vit. D is stratified by race/ethnicity. The alternative, that the better numbers in AAs for both elements is from protective elements overriding the vit. D deficiency effect finds no support in this analysis. The various physiological, biochemical, and environmental factors included in the controls do not exhaust all possibilities but do cover enough ground to conclude that the vit. D effect is a direct one. Many unmodeled environmental elements such as stress are not likely to be protective in AAs.

Periodontitis risk shows no such race/ethnicity effect; the risk is inversely proportional to vit. D levels across the board, once again with all the controls added. This does argue for the risk in AAs being directly from vit. D deficiency. One notes that the unmodeled environmental elements may weaken that conclusion. However, stress and such psychosocial elements are not easy to measure quantitatively and hence are likely to be out of the purview of survey analysis.

We thus see two diametrically opposite conclusions for CVD/bone health and oral health. For the first set, differing vit. D thresholds for different groups do look warranted. As race and ethnicity are both social constructs, how exactly to construct such differing thresholds is a challenge. Most likely, the assessment has to be individualized, based on ancestral genetic markers, context of the living environment, dietary elements, and the like. A single threshold for the entire populace has negative relevance for the CVD risk and bone health elements of AAs, as their risks rise with vit. D levels. The periodontitis data argue against this conclusion. However, PD is less of a burden on quality of life and longevity than CVD or fractures. PD is also more avoidable with improved and admittedly more inconvenient dental hygiene. Vitamin D level recommendations must hence probably focus more on CVD and fracture risks as the relatively greater risks.

This analysis, like most others in the literature, cannot extract causality from the data. NHANES provides cross-sectional data and so only correlations between 25(OH)D and other variables can be estimated. In some cases, such as gum inflammation, there is some evidence that the disease increases local production of vit. D. In CVD and bone fracture risks, there is a possibility that some other element contributes to both vit. D levels and the risk factor. Physiological markers such as blood pressure and BMI; biochemical markers such as fasting plasma glucose level and LDL level; climatic elements such as the season; the demographic factors of age and sex; and socioeconomic variables such as income to FPL ratio and educational level, are all included in this analysis to account for the influence of such hidden variables. The central conclusion of inverse relation across races/ethnicities between 25(OH)D and CVD and bone health risk holds even with all those variables controlled for, and after stratification of data by the groups in question.

Vitamin D levels can be mapped via multiple biomarkers. Free 25(OH)D is believed to be what gets converted to the bioactive 1,25(OH)D form. NHANES, in common with most other surveys, used total 25(OH)D as a measure of vitamin D levels. This marker has been standardized, initially by the US National Institute of Standards, and then by Ghent University, Belgium (Brown et al., 2018). No such reference exists for other markers. The NHANES data used in this study conforms to the reference standard. Free and total BMD are highly correlated Zhang et al. (2019), and hence perhaps this issue is not as relevant as it appears.

The vit. D paradox in AAs has been known for some time. Brown et al. (2018), a report on an expert panel meeting, notes that overall mortality was strongly and inversely linked to 25(OH)D for white Americans but not so for AAs in the 1988 to 2010 NHANES data. The report also mentions that the MESA trial showed lower fracture rates overall for AAs, and a direct relationship between CVD risk and 25(OH)D in AAs as opposed to an inverse relationship in white Americans and some other groups. This study confirms that finding in the 2011–14 NHANES data. Low vit. D levels can lead to osteomalacia in adults, but the level at which the protective effect levels off is 12 ng/mL, well below the recommendation of 30 ng/mL. The expert panel report explicitly notes that the “higher levels of vitamin D [present] hazards for Black Americans, resulting in increasing falls and fractures” (Brown et al., 2018, p. 11). This study confirms the fracture risk increase and general bone health issues, such as decreasing BMD, with increasing 25(OH)D levels in AAs.

For clinical recommendations, information on other effects of vitamin D levels must also be considered, specifically its effects on cancer incidence and growth, reproductive health, metabolic syndrome, COPD, acute respiratory distress, and the general immune system. Nevertheless, the specific harm implied by the statistical relation between higher serum 25(OH)D and worse bone and cardiovascular health in AAs merits a closer look. NHANES no longer collects vit. D levels. Directly assessing the relationship between vit. D levels and markers of the health and illness conditions mentioned earlier requires a nationally representative sample with the scale and scope of NHANES.