1. Introduction
Vitamin D is known for its role in supporting bone growth, a few studies have also shown beneficial links between vitamin D and other non-skeletal health outcomes, including immune function, cardiovascular health and cancer[
1]. Vitamin D functions in the body through both an endocrine mechanism (regulation of calcium absorption) and an autocrine mechanism (facilitation of gene expression)[
2]. However, both vitamin D deficiency and excess can lead to serious health problems. Vitamin D deficiency is the most common nutritional deficiency worldwide. Insufficient vitamin D has been implicated in a range of disorders such as metabolic, autoimmune, psychiatric, and cardiovascular diseases, cancers (especially colon, prostate, and breast cancer), and chronic pain diseases [
3,
4,
5,
6]. Vitamin D as a fat-soluble steroid hormone raised concerns about the toxicity of over-supplementation. Early symptoms of vitamin D toxicity include drowsiness, continuous headaches, irregular heartbeat, and loss of appetite[
2,
7].
The synthesis, absorption and metabolic processes of vitamin D can be influenced by a lot of factors, including environmental factors such as sun exposure, latitude, and season and individual factors such as adiposity, genetics, skin complexion, and age [
8,
9]. Environmental pollutants are also suggested to be able to interfere with human vitamin D[
10]. Epidemiologic studies suggest an association between air pollution and vitamin D deficiency. In cities with high levels of air pollutants (ozone, particulate matter, and sulfur dioxide), UV light is efficiently absorbed by these pollutants, thus reducing the skin synthesis of pro-vitamin D3[
11]. In addition, due to toxic metals are widely distributed in the environment and their potential role in endocrine disruption, it is necessary to investigate the influences of exposure to toxic metals on human vitamin D status[
10,
12]. A few population-based studies have been demonstrated that cadmium (Cd) and lead (Pb) are capable of interfering human vitamin D status. Several studies conducted in cadmium-contaminated areas of Japan found that exposure to environmental cadmium resulted in vitamin D deficiency[
13,
14]. As an endocrine modulator in human populations, lead is also considered as one of the risk factors of vitamin D deficiency. One cohort study of pregnant women shown that vitamin D intake is negatively correlated with maternal blood lead and cord blood lead[
15]. Another cross-sectional study demonstrated a significant negative correlation between blood lead level and serum 1,25(OH)2D in the children exposure to lead[
16].
Arsenic, a well-recognized environmental contaminant, is listed by the Agency for Toxic Substances and Disease Registry (ATSDR) as the first on the list of hazardous substances[
17]. Arsenic is widely distributed in environmental media such as air, water, and soil[
18]. It is worth noting that water-borne arsenic contamination has led to the greatest concern for human health on a global scale[
19]. Long term drinking water with excessive arsenic can induce liver and kidney damage, cardiovascular disease, diabetes and cancer[
18]. Similar to heavy metals such as cadmium and lead, arsenic may affect vitamin D status. It was found that As
2O
3 increased the transcriptional activity of paricalcitol, a noncalcemic vitamin D analogue, probably by decreasing the function of the mitochondrial enzyme 24-hydroxylase, which functions to metabolize the active vitamin D in acute myeloid leukemia (AML) cells[
20]. A cross-sectional study conducted near the lead smelter in city of Torreón in Mexico found positive associations between urinary arsenic and serum 1,25(OH)
2D, but no association between urinary arsenic and serum 25(OH)D[
21]. Another hospital-based prospective cohort study focused on the associations of urine metals during pregnancy with cord serum vitamin D levels and found that urinary arsenic was positively associated with serum 25(OH)D[
22]. However, no studies on the correlation between urinary arsenic and serum vitamin D have been conducted in areas of water-borne arsenicosis.
The key factor in an individual's susceptibility to arsenic is the ability to metabolize and efficiently excrete arsenic[
23]. Inorganic arsenic (iAs), namely the sum of [As (III)+As (V)], accepts a methyl group from S-adenosylmethionine and methylated into monomethylarsonic acid (MMA) or dimethylarsonic acid (DMA) after ingested by humans [
19,
24]. In humans, efficient methylation from iAs to DMA is associated with reduced responsiveness and increased urinary arsenic excretion[
25]. The extent of DMA methylation is different within and between individuals and populations, suggesting a genetic component to arsenic metabolism[
26]. In general, high fractions of MMA and residual non-methylated inorganic arsenic in urine are associated with a higher risk of arsenic-related adverse health effects, suggesting that more efficient arsenic metabolism capacity (a high proportion of DMA) protects the body against arsenic toxicity [
27,
28]. However, it is still not known whether the arsenic metabolism capacity can affect the association of arsenic and vitamin D levels.
As we know, skin is the only organ that has the capacity to produce vitamin D [
29]. Most vitamin D is obtained from ultraviolet B (UVB) exposure of the skin, which stimulates the conversion of 7-dehydrocholesterol into vitamin D
3 (cholecalciferol)[
30,
31]. The skin is also a toxic target organ for arsenic, and skin hyperpigmentation and hyperkeratosis have long been known to be the hallmark signs of chronic arsenic exposure. Therefore, the association of arsenic with vitamin d in a population with skin lesions deserves further study.
Serum 25(OH)D is considered as a reliable marker of vitamin D sufficiency because it has the highest abundance in all vitamin D metabolites in serum and its level remains stable for almost two weeks[
32,
33]. In the present study, we analyzed the association of urinary arsenic species and arsenic metabolism capacity to serum 25(OH)D. We further assessed the association between urinary arsenic species and serum 25(OH)D according to the stratification of skin keratinization due to arsenic exposure.
2. Materials and Methods
2.1. Study site and population
Data for this study was derived from a cross-sectional study of the association between arsenic exposure and diabetes conducted in 2019 in Wenshui County, Shanxi Province, People's Republic of China, which was identified as an area of water-borne arsenicosis according to long-term monitoring data provided by the Shanxi Institute of Endemic Disease Prevention and Control. All the participants were given face-to-face questionnaires and a general physical examination by well-trained staff. General demographic data including general demographic characteristics (age, gender, height, and weight), socioeconomic status (education and occupation), and milk consumption (≤1 glass/week and >1glass/week), because milk is fortified with vitamin D. The body mass index (BMI) of the participants was further calculated with height and weight. Skin hyperkeratosis (a hyperkeratotic thickening of the skin typically on the palms and soles) was ascertained by diagnosis of endemic arsenism (WS/T 211-2015) by study physicians. 762 participants were included for the final analysis after excluding participants with incomplete basic information and abnormal liver and kidney function. This study was approved by the Ethics Committee of Center for Endemic Disease Control, Chinese Center for Disease Control and Prevention, Harbin Medical University, the number of ethical review is HRBMUECDC20210410.
2.2. Sample collection
Blood and urine samples of the participants were collected by professional nurses from the Affiliated Hospital of Shanxi Institute of Endemic Disease Prevention and Control. The venous blood of the cubital fossa was collected through blood collection vessels and the serum and blood cells were centrifuged on the site. The instant urine was collected and divided into 1.5 mL EP tubes after centrifugation. The samples were frozen and stored in the -80°C refrigerator until analysis.
2.3. Determination of urinary arsenic concentrations
Urinary arsenic species were determined using Liquid Chromatography Atomic Fluorescence Spectrometer (LC-AFS6500, Beijing Haiguang Instrument Co., Ltd., China). Chromatographic separation of the arsenic compounds was carried out in accordance with a previously validated method[
34]. The limit of detection (LOD) of As
3+, As
5+, DMA, and MMA were 0.36μg/L, 0.42μg/L, 0.45μg/L and 0.81μg/L, respectively. When the arsenic content in the urine form of the sample is lower than LOD, the value is LOD
.
2.4. Determination of serum 25(OH)D concentrations
Serum 25(OH)D was detected by biochemical analyzer 3100 (Hitachi High-Technologies Corporation, Japan). The vitamin D status of the participants was categorized according to the definition of the kit specification, which defines 25(OH)D <20 ng/mL as deficient, 21-29 ng/mL as inadequate, 30–100 ng/mL as adequate and >100 ng/mL as excessive. More information about reagents is available at
https://www. nbmedicalsystem.com/.
2.5. Determination of blood glucose concentrations
Participants were asked to fast overnight. Fasting blood glucose was measured on the day of blood collection using the glucose tester (Acon Biotechnology co., Hangzhou, China).
2.6. Quality assurance and quality control
Training manual and work manual were developed to train the investigators and standardize the survey methods according to the relevant contents of the questionnaire. Before the test, the instruments and equipment used in the laboratory were corrected in time to reduce the measurement error. Frozen human serum samples were equilibrated at room temperature for 30 minutes prior to use, followed by centrifugation at 12000g for 5 minutes to filter out the precipitate. 1 mL urine stored at −80 °C was thawed at room temperature and centrifuged for 5 min at 13000 rpm. Filtered urine with MCE membrane was been analyzed after centrifugation.
2.7. Statistical analysis
Continuous variables were expressed as mean (± standard deviation) or median (P25-P75). Continuous variables with normal distribution were compared by student t-test, while continuous variables with non-normal distribution were compared by Mann-Whitney U test. Categorical variables were compared by the chi-square test. The concentration of different forms of arsenic (As3+, As5+, MMA, and DMA) was calculated by standard curve. iAs and tAs are calculated by the following formula: iAs = As3++As5+; tAs= iAs + MMA + DMA. All concentrations of urinary arsenic species were transformed using an extension of the Box-Cox transformation to approach normality (tAs: λ = 0.298, iAs: λ = 0.145, MMA: λ = 0.021, DMA: λ = 0.330). Multivariable linear regressions were used to assess changes in serum 25(OH)D for each one-unit increment in the Box-Cox transformed urinary arsenic species. The baseline information and levels of these covariates were adjusted in these equations.
Model 1: Y= β (tAs or iAs or MMA or DMA) + ε
Model 2: Y= β (tAs or iAs or MMA or DMA) + β1 (age) + β2 (gender) + β3 (BMI) + ε
Model 3: Y= β (tAs or iAs or MMA or DMA) + β1 (age) + β2 (gender) + β3 (BMI) + β4 (occupation) + β5 (education) + β6 (skin hyperkeratosis) + β7 (milk consumption) + β8 (blood glucose) + ε
Binary logistic regression models were used to further evaluate the relationship between urinary arsenic species and vitamin D status, in which vitamin D status was determined with serum 25(OH)D and divided into two categories: vitamin D excess (>100 ng/mL) and normal (30-100 ng/mL) according to the instructions. 6 participants with serum 25(OH)D levels below 30ng/mL were excluded because the sample size was too small for follow-up analysis. The effect estimates were presented as β or odds ratios (ORs) with their 95% confidence intervals (95% CIs). The baseline information and levels of these covariates were adjusted in these equations.
Model 1: Y= β (tAs or iAs or MMA or DMA) + ε
Model 2: Y= β (tAs or iAs or MMA or DMA) + β1 (gender) + β2 (skin hyperkeratosis) + β3 (milk consumption) + ε
Stratification of arsenic methylation metabolic capacity by median iAs%, MMA%, and DMA%. All statistical analyses were performed in SPSS 26.0 software and p value < 0.05 (two sides) was considered significant.
3. Results
3.1. Basic characteristics of the study participants
There was a total of 762 participants and 86.1% of participants were farmers. Of these, 256 participants (33.7%) were male and 506 participants (66.3%) were female. The mean (± SD) age of enrolled participants was 57.92 ± 10.80 years and the mean (± SD) body mass index was 25.77 ± 4.01 kg/m
2. The median concentrations of urinary tAs, iAs, MMA, and DMA were 69.81, 3.46, 4.78, and 51.15 μg/L, respectively. The median concentration of blood glucose was 5.70 mmol/L. The mean serum 25(OH)D concentration was 74.03 ± 22.67 ng/mL. These results are presented in
Table 1.
3.2. Multivariate linear regression analysis between urinary arsenic species and serum 25(OH)D
Table 2 reflects the association between urinary arsenic species and serum 25(OH)D. For the urinary tAs, multivariate linear regression analysis found that the Box-Cox transformed urinary tAs was positively associated with serum 25(OH)D in the non-adjusted model (β=0.046, 95%CI: 0.020 to 0.071, p < 0.01). And this positive relationship still existed after adjusted other potential confounders (Model 2: β=0.044, 95%CI: 0.021 to 0.067, p < 0.01; Model 3: β=0.044, 95%CI: 0.020 to 0.069, p < 0.01).
For the urinary iAs, we found a positively associated between Box-Cox transformed urinary iAs and serum 25(OH)D in the non-adjusted model (β=0.330, 95%CI: 0.035 to 0.624, p < 0.01). However, when we adjusted other potential confounders, the association between urinary iAs and serum 25(OH)D became null in model 2 (β=0.155, 95%CI: -0.115 to 0.425, p = 0.260) and model 3 (β=0.100, 95%CI: -0.178 to 0.377, p = 0.482).
For the urinary MMA, our results demonstrated that the Box-Cox transformed urinary MMA was positively associated with serum 25(OH)D in the non-adjusted model (β=0.447, 95%CI: 0.145 to 0.748, p < 0.01). And this positive relationship still existed after we adjusted for age, gender, and BMI (β=0.276, 95%CI: 0.001 to 0.551, p = 0.049). However, when we added potential confounders including occupation, education, skin lesions, and milk consumption, this positive relationship was not existed (β=0.272, 95%CI: -0.013 to 0.556).
For the urinary DMA, our multivariate linear regression analysis found that the Box-Cox transformed urinary DMA was positively associated with serum 25(OH)D in the non-adjusted model (β=0.057, 95%CI: 0.024 to 0.089, p < 0.01). After we adjusted other potential confounders, this positive relationship still existed (Model 2: β=0.060, 95%CI: 0.031 to 0.091, p < 0.01; Model 3: β=0.062, 95%CI: 0.030 to 0.094, p < 0.01).
3.3. Associations between urinary arsenic species and vitamin D status
We excluded 6 participants with vitamin deficiency (serum 25(OH)D < 30ng/mL), and defined participants with serum 25(OH)D levels of 30-100ng/mL as a normal group and those with higher than 100ng/mL as a vitamin D excess group.
Table S1 reveals the basic characteristics of the population between two groups. No significant difference was noted in age, BMI, education, occupations, urinary iAs, urinary MMA, and blood glucose between two groups (all p > 0.05). But we found the gender, skin hyperkeratosis, and milk consumption were varied between two groups (all p < 0.05). We noticed that the populations in the vitamin D excess group had higher urinary tAs and urinary DMA levels compared to the normal group (all p < 0.05).
Table 3 shows the relationship between urinary arsenic species and vitamin D status. For urinary tAs, our binary logical regression analysis showed that the significant increases of 0.4% in the risk of vitamin D excess for every 1-unit increment in the Box-Cox transformed urinary tAs in the non-adjusted model (OR=1.004, 95%CI: 1.001 to1.007, p < 0.01). After adjusting for gender, skin health and milk consumption, this positive relationship was still existed (OR=1.004, 95%CI: 1.001 to1.007, p = 0.045). For urinary DMA, we found that the significant increases of 0.6% in the risk of vitamin D excess for every 1-unit increment in the Box-Cox transformed urinary DMA in the non-adjusted model (OR=1.006, 95%CI: 1.001 to1.010, p < 0.01). After adjusting for above mentioned covariates, this positively relationship was still existed (OR=1.006, 95%CI: 1.001 to1.011, p = 0.017).
3.4. Association between arsenic metabolism efficiency and serum vitamin D
We performed the linear regression analyses to further identify the effect of arsenic metabolism efficiency on the association between urinary tAs and serum 25(OH)D in different populations (
Table 4). Interestingly, no significant correlation was observed between urinary tAs and serum 25(OH)D (β=0.015, 95%CI: -0.027 to 0.057, p = 0.484) in the populations with iAs% below the median. However, in the populations with iAs% above the median, each one-unit increase in the Box-Cox transformed tAs level was associated with increases of 0.064 (95%CI: 0.032 to 0.096, p < 0.01) in serum 25(OH)D. This result suggests arsenic methylation capacity might affect the positive relationship between urinary tAs and serum 25(OH)D levels.
3.5. Association between urinary arsenic species and serum 25 (OH)D in subgroups stratified by skin hyperkeratosis
The skin is both an organ for vitamin D synthesis and a target organ for arsenic toxicity. Therefore, stratified analyses were performed to determine the relationship between urinary arsenic species and serum 25(OH)D in populations with normal and arsenic-induced hyperkeratotic skin.
Table S2 shows the basic characteristics of the population in two groups. We found that the serum 25(OH)D in the skin hyperkeratosis group was higher than that in the normal group.
Table 5 shows the associations between urinary arsenic species and serum 25(OH)D in two groups. In the normal group, the positive associations were found between tAs (β=0.041, 95%CI: 0.013 to 0.069, p < 0.01), MMA (β=0.425, 95%CI: 0.075 to 0.776, p < 0.01), and DMA (β=0.054, 95%CI: 0.017 to 0.090, p < 0.01) with serum 25(OH)D. However, we only found the positive association between urinary iAs (β=0.592, 95%CI: 0.041 to 1.143, p = 0.129) and serum 25(OH)D in the arsenic-induced hyperkeratotic skin group.
4. Discussion
Vitamin D is a group of biologically active steroid compounds. It is found in nature as vitamin D2 (ergocalciferol) and vitamin D3 (cholecalciferol)[
35,
36]. Vitamin D2 is considered a plant form of vitamin D, because it produced by fungi and yeast by UVB-exposure of ergosterol (provitamin D2)[
37]. For most vertebrate animals, 7-dehydrocholesterol is converted under ultraviolet B (UVB) rays to vitamin D3[
38]. When vitamin D2 and vitamin D3 are ingested, they are both metabolized in the liver to form 25(OH)D. As the main storage form of vitamin D, total serum or plasma 25(OH)D concentrations reflect exposure to skin synthesis, dietary intake, and mobilization of tissue stores, therefore 25(OH)D is often considered as reliable biomarker of the nutritional status of the vitamin D endocrine system[
39].
Studies on the association between urinary arsenic and serum vitamin D are scarce. A prospective cohort study focused on association of maternal metals exposure during pregnancy with newborns' vitamin D status in China found a positive correlation between urinary arsenic and serum 25(OH)D.[
22]. In Torreón, a town in Mexico, no association of urinary arsenic with serum 25(OH)D was reported in a case-control study focused on adolescents[
21]. In this study, we found that urinary iAs and MMA was positively associated with serum 25(OH)D in the incompletely adjusted models. However, the positive associations between urinary tAs and DMA with serum 25(OH)D were exist in all linear regression models. The results of binary logical regression analysis also demonstrated urinary tAs and DMA was positively associated with the risk of vitamin D excess. These findings suggest that arsenic metabolism efficiency may affect serum 25(OH)D levels.
Once inorganic arsenic (iAs) is ingested into the body, it will be metabolized in the liver. iAs then accepted a methyl group from S-adenosylmethionine and metabolized to monomethylarsonic acid (MMA) and dimethylarsonic acid (DMA) through a series of oxidation and methylation processes[
19,
40]. These arsenic species (iAs, MMA, and DMA) are excreted from the urine and can be measured and expressed as a percentage of total urinary arsenic (iAs%, MMA%, and DMA%)[
41]. Generally, low iAs% indicates enhanced methylation capacity, while high iAs% indicates decreased methylation capacity. It is hypothesized that increased methylation capacity may reduce susceptibility to arsenic-related toxicity[
41,
42,
43]. The results in present study showed the significantly positively association between urinary tAs and serum 25(OH)D in the populations with iAs% above the median. However, no significant correlation was observed between urinary tAs and serum 25(OH)D in the populations with iAs% below the median. Therefore, our study support that the positive association between urinary arsenic and serum 25(OH)D exist in the incompletely methylated population and this correlation might be altered by efficiency of arsenic metabolism. However, further work is needed to expand sample size for relevant studies.
The maintenance of the superficial structure of the skin is dependent on the balance between the differentiation and proliferation of keratinocytes. Chronic exposure to inorganic arsenic dramatically disrupts this balance, ultimately leading to hyperkeratosis of the skin[
44,
45,
46]. The epidermis is also the major source of vitamin D for the body, because the keratinocytes within the epidermis are capable of producing 25 (OH)D3[
29,
47,
48]. The results of our stratified analyses showed the positive association between urinary iAs and serum 25(OH)D in the arsenic-induced hyperkeratotic skin group. Future work is required to explore the mechanisms.
One potential limitation of our study is its cross-sectional nature, which restricts inferences about causality, and more longitudinal studies are needed. Further, we didn’t collect the information of time spent outside because the sun exposure is a potential confounder, but our study population was homogenous, consisting entirely of adult, 86.1% of whom were farmers and the lifestyle of them is similar. Finally, further in vitro experiments are recommended to verify whether arsenic leads to elevated vitamin D levels by promoting hyperproliferation of keratinocytes.
Author Contributions
Kunyu Zhang: Conceptualization, Methodology, Visualization, Formal analysis, Writing–original draft. Yunyi Yin: Methodology, Software, Formal analysis. Man LV: Formal analysis, Resources, Validation. Xin Zhang: Visualization, Investigation. Meichen Zhang: Methodology, Software, Formal analysis. Jia Cui: Validation. Ziqiao Guan: Validation. Xiaona Liu: Software, Validation, Investigation. Yang Liu: Software, Validation, Investigation. Yanmei Yang: Writing-review & editing, supervision, project administration. Yanhui Gao: Writing-review & editing, supervision, project administration, Funding acquisition.