Abstract
Background:
Clinical consensus differs as to when blood vitamin D (VD) levels should be measured in children. Obesity and metabolic syndrome are risk factors for low VD levels and are also associated with acanthosis nigricans (AN).
Objectives:
To test whether the clinical diagnosis of AN is a strong predictor for vitamin D deficiency (VDD) in children.
Methods:
Within the study period (2015–2020), we identified 677 consecutive individuals (age <18 years) with available calcidiol measurements and compared those with (n = 273) and without (n = 404) AN. Bivariate associations and the occurrence of AN were tested using the chi-squared test. Multivariate logistic regression was performed to control for confounding variables, and adjusted odds ratios with 95% confidence intervals (CI) were reported. Multiple regression analysis was performed, and unstandardized beta coefficients, standard errors, and standardized beta coefficients were reported.
Results:
Individuals with AN had 3.6 times higher odds of VDD than those without (95% CI: 1.38–9.51, P = 0.009). Males had 0.41 times lower odds of having AN than females (95% CI: 0.21–0.79, P = 0.008). Individuals with vitamin D sufficiency (VDS) were much less likely to be diagnosed with metabolic syndrome compared with those who were vitamin D deficient (P = 0.011), even after adjusting for body mass index z-scores.
Conclusion:
Children and adolescents with AN are at a higher risk of VDD and should likely be tested for low calcidiol levels.
Introduction
Acanthosis nigricans (AN) is a common pediatric diagnosis among African American and Latin American children and adolescents and a frequent referral to pediatric clinics from school nurses. Several types of AN have been described in the literature, including those associated with obesity (most common), acral AN in people with dark skin, syndromic AN (including metabolic syndrome and fibroblast growth factor receptor defects), nonsyndromic malignancy-associated AN, familial AN, and drug-induced AN, among others. 1 Although the causes of AN are numerous, it is most commonly associated with insulin resistance. 2 Excessive insulin production causes overstimulation of epidermal keratinocytes and dermal fibroblast growth receptors, leading to abnormalities in cell proliferation and AN. 1
Vitamin D deficiency (VDD) is also common in children, and it is associated with obesity, diabetes, and other components of metabolic syndrome. However, screening guidelines differ as to when to measure calcidiol (25-hydroxyvitamin D [25(OH)D]) levels in children. This knowledge gap is partly explained by the lack of a clear definition of “sufficiency” or acceptable 25(OH)D levels at different stages of puberty, during growth spurts, or among individuals with obesity, especially those with the metabolic syndrome phenotype. Due to the discrepancy between different clinical guidelines as to when to measure 25(OH)D levels, it is essential to determine if there are other strong predictors for VDD.
Interestingly, VDD and AN are both associated with obesity and metabolic syndrome. While studies have shown that AN and VDD occur in the same patient population, 3 very few studies have been conducted to determine if AN can predict the presence of VDD. One clinical study proposed that children and adolescents with obesity and AN had lower VD levels than those who had obesity without AN. 3 However, this study was conducted in a predominately non-Hispanic White population and only included individuals with obesity. Identifying if AN is predictive across minority groups, as well as if other clinical markers can predict VDD, would be an important tool for pediatricians as they decide when VD testing should be performed. Therefore, the present study aimed to test whether the clinical diagnosis of AN is a strong predictor for VDD in children.
Materials and Methods
Patients
A total of 12,959 children and adolescents (95% African American and Latin American) were cared for at the Kids ‘N Teens Clinics P.A., a community health center in Houston, Texas (latitude 29°N), over a 5-year period (April 2015–April 2020). Most high-risk individuals, especially those with AN and obesity [body mass index (BMI) ≥95th tile or ≥30 kg/m2], were routinely recommended to undergo testing for a fasting lipid profile, blood glucose level, insulin level, liver function, and calcidiol (25(OH)D; as measured by the QuestAssureD™ test) level. Individuals who presented with significant chronic idiopathic musculoskeletal complaints, recent long bone fractures, chronic idiopathic fatigue, mood swings, or depression were also routinely recommended for 25(OH)D testing. Written informed consent was not required considering that the study individual's identity was coded and deidentified by the principal investigator. This study was approved by the Baylor College of Medicine Institutional Review Board.
Definitions
VD status was defined by 25(OH)D levels as follows: VDD, <20 ng/mL; vitamin D insufficiency (VDI), 20–29 ng/mL; and vitamin D sufficiency (VDS), ≥30 ng/mL. 4 Although other guidelines have suggested that a 25(OH)D level of <12 ng/mL should be used to define VDD, 5 all agree that a level <20 ng/mL is insufficient. 5 Thus, we have used a broader definition endorsed by the Endocrine Society and frequently used to define VD status in nonbone health outcomes. 6 –8
The definition for the metabolic syndrome was modified from the original definition to include BMI z-scores (BMI adjusted for child age and sex) instead of waist circumference as a measure of obesity.
9
Metabolic syndrome was defined as meeting at least three of the following critical phenotypes: Systolic or diastolic Blood pressure (BP) >95th percentile, >130 mmHg (systolic), or >90 mmHg (diastolic). Obesity, BMI z-score ≥2.0, or BMI >97th percentile. Diabetes or fasting glycaemia ≥100 mg/dL. Oral glucose tolerance test 2-hr blood glucose of 140–200 mg/dL. Dyslipidemia, which is, low serum levels of high-density lipoprotein (HDL) cholesterol and high serum levels of low-density lipoprotein (LDL) cholesterol. Serum HDL cholesterol (mg/dL) <5th percentile, or <50 mg/dL (females) and <40 mg/dL (males).
10
High triglycerides (mg/dL): >95th percentile or ≥150 mg/dL.
10
Clinical measurement
We used a stadiometer to measure height in inches and a balanced standard mechanical scale for weight measurement of all subjects. BMI (kg/m2), its standard deviation z-scores, and percentile scores were calculated for everyone using the Center for Disease Control 2–20 years' weight and height for age data and clinical information as close to the 25(OH)D testing date as possible. The Children's BMI-percentile-for-age Calculator (
Laboratory methods
VD testing was ordered when the primary care provider judged it to be of clinical benefit (“targeted testing”) or when the subject's parents requested testing (“routine testing”). All laboratory assessments were performed using the Quest Diagnostics Laboratory in Houston, TX. We collected clinical information on the following: total blood levels of 25(OH)D (measured using the QuestAssureD vitamin D assay 11,12 ), lipids, insulin, glucose, alkaline phosphatase, calcium, and creatinine, as well as liver function tests. Assay methods and sensitivity in ng/mL have been previously published. 11,12 Homeostasis model assessment of insulin resistance (HOMA-IR) was calculated as follows: fasting plasma insulin level (μU/mL) x fasting plasma glucose (mM/L)/22.5. 13 –15 HOMA-IR status 14 was defined as follows: normal, <2.6; elevated, 2.6–3.8; and significantly elevated, >3.8.
Statistical methods
We used chi-squared tests to compare VD status groups (deficient, insufficient, and sufficient) with categorical parameters. As a categorical variable, we categorized age into three groups to account for possible pubertal differences: <10 years, 10–14 years, and >14 years. We also compared VD status among several continuous variables, including age. For continuous variables, we compared groups using one-way analysis of variance and then, if a significant main effect was detected, we performed a post hoc test (Tukey's test).
We also used chi-squared tests to analyze hypothetical bivariate associations between age, race, sex, and VD status with the occurrence of AN. Separate multivariate logistic regression analyses were conducted to control for confounding variables in outcomes related to VD status (sufficient vs. insufficient/deficient), AN, metabolic syndrome, and low HDL. We also performed a multiple regression analysis to test for multivariate associations between AN, total VD level, age, binary season, binary race, and sex on HOMA-IR values. All analyses were performed using SPSS Version 26 (IBM Corp., Armonk, NY), and statistical significance was considered at an alpha value of 0.05.
Results
Of the 12,959 individuals cared for during the 5-year study period (2015–2020), 745 received 25(OH)D testing results. We excluded 68 (9.1%) individuals from the study whose medical records noted that they had been on VD supplements before testing and those with no BP or height measurement during the intake examination. Furthermore, individuals with disorders that are known to affect the blood VD concentration, such as those with chronic kidney disease (high creatinine), kidney transplant, liver transplant after hepatic failure, syndromes causing intestinal malabsorption, hyperparathyroidism, sarcoidosis, granulomatous diseases, and cancer, as well as those prescribed certain medications (antiepileptics, glucocorticoids, ketoconazole, or cholestyramine), were also excluded. Finally, a total of 677 (90.9%) individuals were included in the study.
Among these, VD levels were sufficient (≥30 ng/mL) in 6.9%, insufficient (20–29 ng/mL) in 30.4%, moderately deficient (10–19 ng/mL) in 52.6%, and severely deficient (0–9 ng/mL) in 10.1% individuals. Primary diagnosis and 25(OH)D test results of the study sample at an initial clinical encounter, including those of patients with and without AN diagnosis, are presented in the supporting Supplementary Table S1 for reference.
Multiple parameters were associated with VD status in our study groups, potentially confounding the association of VD status with AN (Table 1). Among the VD status groups, we found statistically significant main effects of age, race, sex, season, AN, metabolic syndrome, low HDL, BMI, BMI z-score, systolic BP, diastolic BP, calcium, and alkaline phosphatase. No significant differences between the VD status groups for high LDL or high triglycerides were found. Participants in the VDD group were significantly older than those in the VDI (P = 0.01) and VDS (P = 0.002) groups. Patients with VDD had higher BMI z-score values than those in the VDI and VDS groups (P = 0.002) and significantly higher systolic and diastolic BP values than those in the VDI group (P < 0.001 and P = 0.019, respectively). Serum calcium levels were significantly lower in the VDD group than in the VDI (P = 0.002) and VDS (P = 0.001) groups. No differences were found among the VD status groups for LDL, HDL, triglycerides, creatinine, aspartate aminotransferase, alanine aminotransferase, and hemoglobin A1C levels (Table 1).
Comparison for Vitamin D Status Groups
Statistically significant.
ALT, alanine aminotransferase; AST, aspartate aminotransferase; BMI, body mass index; BP, blood pressure; HDL, high-density lipoprotein cholesterol; LDL, low-density lipoprotein cholesterol; SD, standard deviation; VD, vitamin D.
We also determined parameters associated with AN (Table 2). The group with AN had significant difference in terms of age, VD status, and BMI z-scores. The highest rates of AN were found in individuals 14–18 years of age and those with VDD. There were no differences among the groups when considering race or sex.
Comparison of Acanthosis Nigricans Groups
Statistically significant.
AN, acanthosis nigricans.
Given the multiple factors associated with AN, we determined if the presence of AN was a significant predictor of VDI or VDD. Indeed, after adjusting for age, race, season (winter/spring vs. summer/fall), sex, and BMI z-scores, we found that those with AN had a 3.62 times higher risk of having VDD/VDI (adjusted odds ratio (AOR): 3.62; 95% confidence intervals (CIs): 1.38–9.51, P = 0.009] (Table 3). Despite AN being associated with BMI z-scores, BMI and BMI z-scores were not independent predictors of VDS versus VDI/VDD.
Logistic Regression Model Parameters
Statistically significant.
AOR, adjusted odds ratio; CI, confidence interval; HDL, high-density lipoprotein.
For AN, VDS status, BMI z-score, and age were statistically significant predictors, whereas race, season, and sex were not (Table 3). Individuals in the VDS group had 0.20 times lower odds of having AN when compared with those in the VDD group (P = 0.010) after adjusting for potential confounders. VDI status (P = 0.528) was not a significant predictor of AN when compared with VDD (Table 3).
Our logistic regression model without adjusting for BMI z-scores suggested that VDI individuals were 62% less likely than VDD individuals to have lower HDL. However, after adjusting for BMI z-scores in the model, we found that VDI was a weak predictor (P = 0.06), and VDS was not a significant predictor. This suggests that obesity (high BMI z-scores) is the strongest risk factor (AOR: 2.3; 95% CI: 1.82 − 2.96, P < 0.001) for low HDL levels (Table 3).
In our study, metabolic syndrome was significantly associated with VD status, age, and BMI z-scores (Table 3). However, race, season, sex, and VDI status were not significant predictors. Importantly, individuals in the VDS group were less likely to have metabolic syndrome when compared with those in the VDD group (P = 0.011), even after adjusting for BMI z-scores.
Since AN has been linked to insulin resistance, we tested whether AN and VD levels were significant predictors for HOMA-IR and insulin levels. Only 187 individuals had appropriate laboratory results for HOMA-IR calculation. Thus, AN was a significant predictor of HOMA-IR as well as insulin levels (Supplementary Table S2). This association between AN and HOMA-IR was lost in a further sensitivity analysis when BMI-z scores were added (resulting in neither AN nor BMI z-scores being significant independent predictors of HOMA-IR). These data indicate that the effects of AN on HOMA-IR are likely influenced, partly, by BMI-z score status. Black race was a significant predictor for HOMA-IR but not for insulin levels. Lower VD levels did not predict higher HOMA-IR scores or insulin levels (Supplementary Table S2).
Discussion
VDD, as determined through the measurement of 25(OH)D levels, is highly prevalent in our society. In a study of 15,787 adults in the USA, 12.7% were found to have VDD (<20 ng/mL). 8 The 2006 National Health and Nutrition Examination Survey determined that 18% of children 1–11 years of age had VDD (<20 ng/mL). 16 In a recent study conducted in Dallas, Texas (USA-latitude 34°N), where adequate sun exposure is not a problem, ∼50% of the participants were found to have VDD. Additionally, this study noted that VDD was more clinically significant among African American and Latin American participants. 17 In our study, individuals without AN had a slightly higher prevalence of VDD than the Dallas cohort (57.7% vs. 50%), but those with AN had a much higher prevalence (71.4%). However, global testing for VD is not recommended. 18 Thus, determining possible clinical markers to guide testing would be useful.
There are a few clinical studies that consider AN as a phenotypic predictor of VDD. This is clinically important considering that there is a much higher prevalence of VDD among people with obesity; hence, determining the target population for VD testing would be useful. A previous study showed that AN was associated with lower VD levels but did not determine if AN truly predicted VDI or VDD. 3 Our results showed that the clinical diagnosis of AN significantly predicts VDD and thus suggests that the presence of AN may be a proxy for VDD. Only 6 (2.2%) individuals diagnosed with AN had a sufficient VD level (≥30 ng/mL). Even with using a 25(OH)D level of 20 ng/mL to define sufficiency, almost three-quarters of individuals with AN had VDI in our study. This is not secondary to the presence of obesity or race because after adjusting for potential confounders, including BMI-z scores, those with AN were 3.6 times more likely to have VDD or insufficiency than those without AN.
As expected, individuals with darker skin color were also more likely to be VDD or insufficient. A small study of 150 individuals with obesity showed that AN was present in 59.3%, and of those 89 individuals with obesity and AN, 39.3% also had VDD. 19 Our study has a much higher prevalence of VDD with AN likely due to the different racial composition. While our study consisted of mostly individuals of darker skin complexity, the former study consisted of mostly white adolescent individuals.
The association of VDD with obesity 17,19 and insulin resistance 3,19 has been reported in the literature. Williams et al 19 found that “hyperinsulinism” was a statistically significant covariate predictor of VDD (<20 ng/mL). However, high insulin levels were a much weaker predictor of VDD (OR = 1.04, P = 0.01) when compared with race (OR = 4.96, P = < 0.01), season of the year (OR = 3.69, P = 0.01), and location (OR = 5.91, P = 0.05). 18 AN has also been reported in the literature to be an important phenotypic expression of hyperinsulinism, 20 –22 although this relationship appears to be complex and likely mediated by the presence of obesity in children and adolescents. 21 Interestingly, a study found that HOMA-IR could only predict variations in 25(OH)D levels when AN was included in the model. 3 Similar to Soliman et al, 20 we found that the presence of AN was a strong predictor of higher HOMA-IR scores and higher insulin levels. While this significance partly depends on BMI z-scores, 21 it is likely that we were underpowered to detect independent associations of both AN and BMI z-scores in the same model, given the smaller sample size for HOMA-IR measurements.
Evidence exists to suggest a causal relationship between AN and VDD. VD receptors (VDR) are found in the skin as well as other multiple organs, including the liver, kidney, and pancreas. VD is known to have anti-inflammatory effects by modulating the function of T and B lymphocytes, 23 which can augment insulin action in the cells and regulate brain neuronal excitability through the insulin-sensitive phosphoinositide 3-kinase pathway. 24 Moreover, it has been suggested that 1,25-dihydroxyvitamin D (the most potent endogenous VDR ligand) and the VDR promote gene signals that inhibit keratinocyte proliferation and stimulate the differentiation (keratinization) of basal keratinocytes. 25 Taken together, the literature suggests that low VD levels may promote excessive keratinocyte proliferation, causing the appearance of AN, which could be worsened in the presence of excess insulin signaling.
The cross-sectional design of our study did not permit us to look at AN severity as a predictor of VDD and/or lower VD levels. A previous study suggests that the severity of AN predicts lower 25(OH)D levels, 3 but since they did not categorize by VD status, it is not clear which level of AN severity would help to predict VDD. A large prospective study utilizing validated severity scores is needed to further address this question.
Our clinical practice of targeting VD testing to high-risk subjects uncovered a high prevalence of VDD (62.6%). Interestingly, there were very few subjects with VDD in those tested for routine reasons, validating our clinical approach. Demographic variables, age (as a proxy for pubertal status), obesity, and time of the year when the VD D test was conducted, 3 are also important variables that were adjusted for in the final analysis. We found that the prevalence rate of VDD was much higher in winter than in summer. Additionally, children 10–14 years of age were also at higher risk of VDD/VDI. It is conceivable that pubertal changes in these individuals may be contributing to their pathology; however, since we did not have accurate sexual maturation ratings (Tanner staging) at each visit, we could not test this hypothesis. Furthermore, we detected high levels of VDD in our younger age group, which likely reflects their increased risk due to obesity and darker skin.
There is controversy regarding the clinical diagnosis of metabolic syndrome among children and adolescents. We modified existing criteria using a more practical clinical parameter (BMI z-scores ≥2) as a “proxy” for abdominal obesity, since waist circumference is not a routine vital sign in our clinic. This methodology, although not validated for the diagnosis of metabolic syndrome, is practical, simple, and may be useful for future pediatric obesity research. Our adjusted logistic regression model suggests that children and adolescents with VDD were more likely to have metabolic syndrome than those with VDS. Since this study was observational, we could not test whether VDD causes the development of metabolic syndrome (or vice versa). Obesity is an important phenotype of metabolic syndrome 26 and has also been inversely associated with blood levels of bioavailable VD. 27
Similarly, in our study, individuals in the VDD group had higher BMI z-scores than those in the VDS and VDI groups. Although there are many theories as to why VDD is more likely to occur as BMI values increase, there is no definitive explanation. Despite the association of low VD levels with increasing BMI, we did not observe associations of low VD with abnormal lipid, liver, or glucose parameters. We found an association (independent of BMI) between low VD levels and higher BP readings, but this analysis cannot explain if this relationship is causal or not. Several studies have suggested an important association between VDD and metabolic syndrome in adults, 28 –30 and an inverse association between VD levels and metabolic syndrome. 29,30 Thus, further study is required to understand the relationship between these disorders.
Overall, our results suggest that AN is a reliable proxy marker for VDD, especially in minority young individuals. However, there are many interesting unresolved clinical questions pertaining to our research. Further study is required to explore whether abnormalities in lipoprotein metabolism (low HDL phenotype) contribute to low VD levels, possibly through sequestration in the adipose tissue. Additionally, causal links and possible therapeutic uses of VD for other abnormal parameters, such as BP, obesity, AN, and lipids in children and adolescents with metabolic syndrome, still need to be elucidated.
Study limitations and strengths
Our study was observational and thus could not prove causality. Additionally, our population was skewed toward high-risk groups, especially those with obesity (more likely to have insulin resistance) and darker skin color. Thus, our study may not be generalizable to the entire US population. However, our study targeted the most at-risk populations for VDD. Due to the retrospective and cross-sectional nature of our research, we were unable to document sunlight exposure, skin phototypes, severity of AN, accurate diet history, or adherence to VD supplementation. Additionally, as noted above, we were unable to obtain other laboratory values, such as parathyroid hormone levels. Although we used standard measurement practices for vital signs, no interobserver variability was measured, and thus minor vital sign observer bias cannot be excluded. On the other hand, subject selection bias was reduced by including all subjects who underwent a 25(OH)D test during the 5-year study period.
Conclusions
Our results suggest that obesity-associated AN appears to be a reliable proxy for VDD. Future prospective studies to confirm these findings and determine whether VD treatment in patients with AN would improve the AN itself and/or its related comorbidities are required.
Footnotes
Acknowledgments
Kids ‘N Teens Clinics PA (
Authors' Contributions
F.A.I. conceived and designed the study, acquired funding, administered the project, and drafted the article. F.J.I.-I. participated in the conception and design, data acquisition, and managed resources. E.R.H. performed the analysis and interpretation of data. S.S. participated in the conception and design of the study, revised the article for key critically important intellectual content, and mentored coauthors. All authors participated in the article drafting, editing, and approval processes.
Author Disclosure Statement
S.S. receives speaking and consulting fees from Rhythm Pharmaceuticals. F.A.I., F.J.I.-I., and E.R.H. report no conflicts of interest.
Funding Information
S.S. receives salary support in part by agreement no. 58-3092-0-001 of the United States Department of Agriculture (USDA), Agricultural Research Service, and from the American Diabetes Association grant no. 1-17-JDF-037.
Supplementary Material
Supplementary Table S1
Supplementary Table S2
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
