https://doi.org/10.5812/hepatmon.92992

The Association Between Serum Vitamin D Level and Nonalcoholic Fatty Liver Disease

authors:

avatar Elham Ehrampoush ORCID 1 , 2 , 3 , avatar Alireza Ghaemi ORCID 4 , avatar Saeed Osati 5 , avatar Maryam Maleki 6 , avatar Hadi Raeisi Shahraki ORCID 7 , avatar Behnam Honarvar ORCID 2 , avatar Maryam Bahmanyar ORCID 1 , avatar Reza Homayounfar ORCID 1 , 2 , 3 , *

Noncommunicable Diseases Research Center, Fasa University of Medical Sciences, Fasa, Iran
Health Policy Research Center, Institute of Health, Shiraz University of Medical Sciences, Shiraz, Iran
Department of Nutrition, Fasa University of Medical Sciences, Fasa, Iran
Department of Basic Sciences and Nutrition, Health Sciences Research Center, Addiction institute, Faculty of Public Health, Mazandaran University of Medical Sciences, Sari, Iran
National Nutrition and Food Technology Research Institute, Faculty of Nutrition Sciences and Food Technology, Shahid Beheshti University of Medical Sciences, Tehran, Iran
Student Research Committee, Fasa University of Medical Sciences, Fasa, Iran
Department of Epidemiology and Biostatistics, Faculty of Health, Shahrekord University of Medical Sciences, Shahrekord, Iran
Hepatitis Monthly: Vol.19, issue 9; e92992
published online: October 2, 2019
article type: Research Article
received: May 5, 2019
revised: August 17, 2019
accepted: September 19, 2019

How To Cite Ehrampoush E, Ghaemi A, Osati S, Maleki M , Raeisi Shahraki H, et al. The Association Between Serum Vitamin D Level and Nonalcoholic Fatty Liver Disease. Hepat Mon. 2019;19(9):e92992. https://doi.org/10.5812/hepatmon.92992.

Abstract

Background:

Nonalcoholic fatty liver disease (NAFLD) is a condition, in which triglyceride accumulates in hepatic cells without a history of alcohol consumption and is strongly associated with insulin resistance, obesity, diabetes mellitus, hypertension and dyslipidemia. The potential role of vitamin D in the pathogenesis of NAFLD has been reported.

Objectives:

The aim of this study was to determine the optimal vitamin D levels for prevention of NAFLD.

Methods:

In a cross-sectional study, 2,160 cases who referred to a university-affiliated health center were randomly selected and their demographic information, anthropometric and metabolic indices and also vitamin D levels were collected. Fatty liver was assessed by fatty liver index (FLI) and confirmed by FibroScan using controlled attenuation parameter (CAP). Based on the NAFLD score, the subjects were divided into two groups and the vitamin D cutoff point was calculated by ROC curve.

Results:

Based on the results, 745 patients (34.5%) had different degrees of fatty liver. Significant differences in the stiffness of liver tissue were observed between vitamin D categories (285.10 ± 30.56 for severe, 251.82 ± 42.63 for moderate and 201.02 ± 36.08for mild deficiency). According to the multivariate analysis, age, fasting insulin and vitamin D levels were found as the most significant factors in NAFLD pathogenesis. Vitamin D cutoff point was obtained 18 nmol/L in women and 21 nmol/L in men.

Conclusions:

The results indicated a significant association between vitamin D level and NAFLD score. Accordingly, increasing the public awareness to maintain a proper level of vitamin D may be a preventative strategy against NAFLD.

1. Background

Nonalcoholic fatty liver disease (NAFLD) is characterized by hepatic lipid accumulation in hepatocytes without a history of alcohol consumption (1). It is the most common liver disease worldwide with the prevalence of 25-35% in different populations (2). The total prevalence of NAFLD and also the prevalence of mild, moderate and severe fatty liver disease were estimated 33.9% (95% CI: 26.4% - 41.5%), 26.7% (95% CI: 21.7% - 31.7%), 7.6% (95% CI: 5.7% - 9.4%), and 0.5% (95% CI: 0.1% - 0.9%), respectively (3, 4). People with the age range of 40 - 60 years or in some cases younger, are more likely to develop NAFLD.

Epidemiological studies have demonstrated that regardless of the interferer factors, such as obesity and insulin resistance (IR), there is an association between hypovitaminosis D and the presence of NAFLD and steatohepatitis (NASH) (5-7). Furthermore, anti-fibrotic, anti-inflammatory and insulin-sensitizing effects of vitamin D on hepatocytes have been confirmed by several studies (8). Accordingly, it seems that vitamin D possibly plays a critical important role in the pathogenesis of NAFLD (1). It has been shown that the reduced level of 25-hydroxy vitamin D (25(OH)D) is associated with the progression of NAFLD (2, 9-11). In some studies, however, no association was found between vitamin D levels in patients with non-cirrhotic livers and those with liver cirrhosis or also between vitamin D levels in patients with NAFLD and healthy individuals (12). Nevertheless, previous studies have faced some limitations, including a small sample size and no assessment for major interferer factors, such as obesity, insulin resistance, physical activity and energy intake. A recently published study in Asia shows that the association between vitamin D level and NAFLD is varied among the Asian and non-Asian populations (13), which can be due to the differences in the vitamin D receptor alleles among the Asian and Western societies. The results of a recently published systematic review on 5,000 patients with NAFLD and 8,000 healthy individuals (control) demonstrated that there is a vitamin D deficiency in NAFLD patients, indicating that this vitamin could be a part of the pathogenesis of NAFLD. NAFLD patients are 26% more likely to have vitamin D deficiency. This relationship is more obvious by comparing the Western and Eastern societies (14). The increased number of new cases of NAFLD in Iran indicates the importance of raising public awareness and also developing appropriate preventive and therapeutic strategies for management of NAFLD.

2. Objectives

Due to the high prevalence of NAFLD, inconsistent results of different relevant studies and also the effects of environmental and nutritional factors on the pathogenesis of the disease, this study aimed at evaluating the prevalence of vitamin D deficiency in patients with NAFLD and healthy individuals.

3. Methods

This study is a part of a larger project to investigate the factors affecting the development of nonalcoholic fatty liver in an Iranian sample population (Tehran). For this purpose, the information of 2,160 people who referred to a university affiliated nutrition consolidation center following a recall from the early April 2016 to the late September 2017 (a period of little over a year) was randomly selected and investigated. The inclusion criteria included alcohol consumption less than 30 g/d for men and 20 g/d for women and the subjects’ willingness to participate in the study (15). Exclusion criteria also included those with diabetes, cancer, hypothyroidism or other disorders or taking medications that influence the metabolic status of the body. In addition, those who changed their diet over the last six months and consumed the supplements containing vitamin D, were also excluded.

The study protocol was in accordance with the Helsinki Declaration and confirmed by the Ethics Committee of Fasa University of Medical Sciences (approval code: IR.FUMS.REC.1396.267). The participants were informed about the research objectives and the written informed consent was obtained from the subjects before starting the survey. The study was supported by the Fasa University of Medical Sciences (grant No.: 96050).

The demographic characteristics, including age, gender and woman’s menopausal status were asked from the participants. The weight and height of the subjects were measured and recorded by a trained researcher using a stadiometer accurate to 0.1 cm and the Seca 767 digital scale (± 0.1 kg; Seca; Japan), respectively. Body mass index (BMI) was also calculated as weight in kilograms divided by height in meters squared.

Blood pressure was measured 15 minutes after rest in the sitting position two times within a 15-minute interval from the right arm of the subjects. The mean value of these two measurements was considered as the final blood pressure.

The subjects’ food intake data over the past year were collected by the semi-quantitative food frequency questionnaire (FFQ). FFQ includes a list of 168 food items and also portion size information as standardized portions. The subjects were asked if they had eaten a particular food item (standardized portions) within the last year. The collected values for each food were converted into daily intake (g/day) using a household scales guidebook, based on which the energy intake of the subjects was calculated.

To measure the levels of serum lipids, blood sugar and insulin, 10 cc venous blood sample was drawn after 12 hours of fasting. The blood sample was then transferred to the laboratory under appropriate conditions to perform biochemical tests. High-density lipoprotein (HDL), triglyceride (TG) and total cholesterol levels were measured using enzymatic method by commercially available kits (Pars Azmoon Laboratories, Iran). In addition, the Friedewald formula (FF) was employed to calculate low-density lipoprotein (LDL) level (16). Eventually, the total values of serum lipids were reported as mg/dl. Using 75 g glucose, we measured 2-hour plasma glucose and insulin.

Homeostatic model assessment (HOMA) is used for assessing IR and β-cell function (7).

25(OH)D was measured through the enzyme-linked immunosorbent assay (ELISA) kits (item No.: 501050; Cayman Chemical Co., USA). Its levels were then classified based on the results of Lips study (17). The values below 12.5 nmol/L, between 12.5 and 25 nmol/L, between 25 and 50 nmol/L and higher than 50 nmol/L were considered as severe vitamin D deficiency, moderate vitamin D deficiency, mild vitamin D deficiency and vitamin D replete, respectively.

Fatty liver index (FLI) was used to measure liver steatosis (18). The participants were assigned into three groups based on the FLI values: FLI < 30 was defined as not having NAFLD, FLI ≥ 60 was defined as having NAFLD and patients with FLI of 30 - 59 were defined as having intermediate FLI. The results were then confirmed (Pearson correlation = 0.698, P value < 0.001) by an Gastroenterologist with the extent of stiffness of liver mass using controlled attenuation parameter (CAP) measured with FibroScan XL probe, which is a more reliable diagnostic tool for the obese people and those who have high fat levels (19-21).

FLI = [e0.953 × loge (TG) + 0.139 × BMI + 0.718 × loge (GGT) + 0.053 × waist circumference - 15.745)]/[1 + e0.953 × loge (TG) + 0.139 × BMI + 0.718 × loge (GGT) + 0.053 × waist circumference - 15.745] × 100, with triglycerides measured in mmol/l, GGT in U/L, and waist circumference in cm.

3.1. Statistical Methods

The data were considered to be normally distributed. Based on the normal distribution results, for descriptive statistics, the mean ± SD and No. (%) were used for quantitative and qualitative variables, respectively. Univariate analysis was performed using the chi-Square test or independent samples t-test. To assess the simultaneous effect of variables on NAFLD, logistic regression was implemented and the linear regression was also applied to find a model based on effective factors on NAFLD. For logistic regression, we proceeded in two steps; initially, only the vitamin D variable was entered. In the next step, all related variables were analyzed by Forward: conditional method and the amount of variance explained in each step was investigated. Two model was selected. The first model with three variables and the second model with four variables explained 82.5% and 82.5% of the variance, respectively. Accordingly, the second model with three variables was selected as the optimal model. To assess the model accuracy, the probability of NAFLD for each participant was estimated using the proposed model and the optimal cut-off point was obtained using the receiver operating characteristic (ROC) curve analysis. All the statistical analyses were performed by MedCalc 14.0 (Medcalc software, Ostend, Belgium) or SPSS 20.0 (IBM Co., Armonk, NY, USA).

4. Results

2,160 subjects participated in the study, of whom 745 subjects (34.5%) had NAFLD, whereas 1,415 subjects (65.5%) were healthy. The demographic characteristics of the subjects based on the groups are shown in Table 1. Almost all studied variables, except for the gender and smoking showed a significant difference between two groups. The odds ratio values in the univariate model were calculated for all factors associated with fatty liver. It can be stated that all variables were found to increase the risk of NAFLD.

Table 1.

Demographic and Basic Information of the Participants in Healthy and Patient Groupsa

Healthy People (N = 1415)NAFLD Patient (N = 745)P ValueOR95% CI
LowerUpper
Age, y35 ± 949 ± 8< 0.0011.191.1731.207
Weight, kg80.88 ± 12.2796.16 ± 11.50< 0.0011.1031.0931.113
BMI, kg/m²25.43 ± 4.1030.95 ± 3.49< 0.0011.4131.3711.457
Waist circumference, cm96.49 ± 16.04118.10 ± 12.77< 0.0011.0971.0881.106
Physical activity, MET/24h24.9 ± 5.018.3 ± 4.5< 0.0010.7620.7440.781
Energy, kcal2405.49 ± 413.732664.98 ± 460.11< 0.0011.0011.0011.002
FBS, mg/dL99.39 ± 13.53118.29 ± 13.15< 0.0011.1021.0921.111
Fast Insulin, µU/mL8.16 ± 3.7814.17 ± 3.39< 0.0011.4961.4461.547
Two hour glucose, mg/dL115.22 ± 19.05130.51 ± 20.18< 0.0011.0391.0341.044
Two hour insulin, µU/mL40.28 ± 12.4454.67 ± 12.40< 0.0011.0931.0841.103
HOMA-IR0.92 ± 0.241.17 ± 0.23< 0.00170.60545.294110.060
HOMA-B68.18 ± 15.9079.28 ± 18.02< 0.0011.0401.0341.046
LDL, mg/dL93.72 ± 15.56113.71 ± 16.61< 0.0011.0771.0691.085
HDL, mg/dL48.18 ± 7.9538.88 ± 7.08< 0.0010.8580.8460.871
TG, mg/dL186.57 ± 26.15214.28 ± 25.03< 0.0011.0411.0371.045
Total cholesterol, mg/dL176.01 ± 17.11195.40 ± 16.18< 0.0011.0721.0641.079
Systolic blood pressure, mmHg12.47 ± 1.6313.09 ± 1.66< 0.0011.2561.1881.327
Diastolic blood pressure, mmHg8.04 ± 1.028.84 ± 0.96< 0.0012.2051.9972.435
ALT, IU/L39.77 ± 13.6559.38 ± 13.50< 0.0011.1041.0941.113
AST, IU/L32.40 ± 13.8749.51 ± 12.62< 0.0011.0941.0851.103
GGT, IU/L23.64 ± 10.1436.86 ± 12.98< 0.0011.1021.0921.113
Vitamin D, nmol/L26.77 ± 8.2615.84 ± 5.50< 0.0010.8040.7890.820
Gender0.0010.7410.6210.886
Female656 (62.1)401 (37.9)
Smoking0.9141.0130.7991.284
Yes240 (17)125 (16.8)

Abbreviations: ALT, alanine aminotransferase; AST, aspartate aminotransferase; BMI, body mass index; FBS, fasting blood sugar; GGT, gamma-glutamyl transferase; HDL, high-density lipoprotein; HOMA-IR, homeostasis model of insulin resistance; HOMA-β, Homeostatic model assessment of β cell function; LDL, low-density lipoprotein; TG, Triglyceride.

aValues are expressed as mean ± SD or No. (%).

In the studied subjects, no participant had adequate levels of vitamin D and 845, 1,055, and 260 subjects had mild, moderate, and severe vitamin D deficiency, respectively.

Using regression model, we tried to determine the role of vitamin D in the development of NAFLD. For this purpose, different variables (age, sex, woman’s menopausal status, BMI, waist circumference, waist-to-hip ratio (WHR), waist-to-height ratio (WHtR), physical activity, fasting blood sugar (FBS), 2 hours blood glucose, fasting insulin, 2 hours insulin, HOMA-IR, LDL, HDL, total cholesterol, TG, systolic blood pressure, diastolic blood pressure, alanine aminotransferase (ALT), aspartate aminotransferase (AST), gamma-glutamyl transferase (GGT), total energy intake and vitamin D) were introduced into a stepwise model and the level of inferential variance was examined by the models. Three regression models were chosen. In the first model, only vitamin D was introduced. In this model, by considering the mild vitamin D deficiency as the reference, the risk of NAFLD for the groups with moderate and severe vitamin D deficiencies obtained about 5 and 13, respectively. In the second model, in addition to vitamin D, age and fasting insulin variables were also introduced and the model explained about 82.5% of the variance. Those with moderate and severe vitamin D deficiencies were 4.2 and 2.4 times more at risk for developing NAFLD, respectively. In the third model, ALT was added and the model explained 82.5% of the variance. It was found that those with moderate and severe vitamin D deficiencies were 3.6 and 2.1 times more at risk for developing NAFLD compared to the mild deficiency (Table 2).

Table 2.

Relationship Between Vitamin D Levels and the State of Disease

Model 1aOR (95% CI)Model 2bOR (95% CI)Model 3cOR (95% CI)
Age, y-1.144 (1.114 - 1.175)1.079 (1.058 - 1.099)
Fasting insulin, µU/mL-1.284 (1.189 - 1.386)1.183 (1.131 - 1.238)
ALT1.029 (1.018 - 1.040)
Vitamin D mild deficiencyRefRefRef
Vitamin D moderate deficiency5.858 (4.144 - 8.281)4.257 (2.498 - 7.255)3.616 (2.104 - 6.214)
Vitamin D severe deficiency13.914 (10.149 - 19.075)2.416 (1.669 - 3.497)2.148 (1.476 - 3.126)

aLogistic regression with only one predictor variable.

bMultivariate logistic regression with 3 variables that explains 82.5% of model variance.

cMultivariate logistic regression with stepwise method that excluded Sex, Menopause, Waist cicumference, Waist to hip ratio, Waist to height ratio, physical activity, FBS, 2 hours Glucose, HOMA-IR, LDL, HDL, total cholestrol, systolic blood pressure, diastolic blood pressure, ALT, AST, GGT and total energy intake that explains 82.5% of models variance.

Table 3 shows the cut-off point for vitamin D in men and women, vitamin D levels lower than 18 (nmol/L) for women and 21 (nmol/L) for men significantly increase the risk of NAFLD (Figure 1).

Table 3.

Receiver Operating Characteristic Curve Analysis for Defining the Ideal Vitamin D Cut-Off Point Based on Nonalcoholic Fatty Liver Diseasea

Cut PointSensitivitySpecificityAUCPPVNPVPLRNLR
Female18.2576.0680.6486.566.9386.733.900.30
Male21.581.4077.7385.865.3189.023.660.24

Abbreviations: AUC, Area under the curve; PPV, Positive predictive value; NPV, Negative predictive value; PLR, Positive likelihood ratio; NLR, Negative likelihood ratio.

aOptimum cut-off was defined as the point on the ROC curve with the maximum sensitivity + specificity (Youden method).

ROC curves for NAFLD using vitamin D. (Left) ROC for males (right) ROC for females.
ROC curves for NAFLD using vitamin D. (Left) ROC for males (right) ROC for females.
Figure 1.

ROC curves for NAFLD using vitamin D. (Left) ROC for males (right) ROC for females.

5. Discussion

The results showed an association between 25(OH)D and NAFLD independent variables, including anthropometric indicators, energy intake, physical activity and biochemical indices. Furthermore, the developed models for predicting NAFLD based on the various variables suggested a significant correlation between NAFLD and age, fasting insulin and the level of 25(OH)D. In addition, our study indicated that 20 - 25 nmol/L of 25(OH)D is an acceptable level to prevent NAFLD in men and women.

A recently published systematic review including 45 cross-sectional studies, revealed that 29 studies (64.4%) reported an inverse relationship between vitamin D levels and NAFLD, whereas 16 studies (35.6%) were unable to demonstrate a correlation between these two (22). The optimal threshold for vitamin D in chronic liver disease has reported 10 to 32 ng/mL in (22), which is equal to 24 - 80 nmol/L.

NAFLD is a multifaceted condition, affecting various aspects of the patient’s life, in which several factors, including genetics (23), lifestyle (24), diet (25) and anthropometric factors (26) are involved. Vitamin D has recently been attracted attention as NAFLD risk factor as numerous studies have investigated the association between the levels of serum vitamin D and the risk of NAFLD. Consistent with our results, almost all studies have found that vitamin D deficiency is correlated with the increased risk of NAFLD (2, 6). Bril et al. study on the relationship between NAFLD and serum vitamin D level showed no correlation between these variables, which is not consistent with our results (27). However, since the majority of the relevant studies have not been able to study all affecting variables and modify them, no reliable conclusions have yet been provided. On the other hand, several studies evaluated only a limited number of factors, including ALT to diagnose NAFLD, which reduces the accuracy of their studies.

Various mechanisms have been proposed for the relationship between vitamin D and the incidence of NAFLD. Vitamin D plays its significant intracellular role through bonding with vitamin D receptor (VDR). It has been shown that VDR influences the expression of over 200 genes involved in metabolism, of which inflammation (28) or cellular differentiation (29) should be considered as inflammation plays a significant role in developing NAFLD (30). VDR is also expressed by macrophages. In animal studies, it has been shown that the active form of vitamin D is able to enhance stability of mRNA of the IκB-α (inhibitor of NF-κB) and reduce its phosphorylation, indicating the anti-inflammatory effect of vitamin D in macrophages (31).

The results of genome-wide association studies (GWAS) on patients with NAFLD have shown that the mutation in the vitamin D-binding protein gene (as the main vector of vitamin D) is one of the four major polymorphisms associated with the risk of NAFLD (32). It has been shown that VDR gene expression has an inverse relationship with the intensity of NAFLD (33). On the other hand, it has been shown that vitamin D levels have a positive and significant correlation with the level of adiponectin, independent of BMI (34). Adiponectin as an anti-inflammatory factor has an inverse relationship with NAFLD severity (35).

Our proposed model for the variables affecting the development of NAFLD suggests that age, fasting insulin level and vitamin D status are the strongest factors to predict incidence of NAFLD. In this regard, using the minimum number of variables, we could predict the risk of NAFLD. The area under the ROC curve for this model obtained 0.96 (95% CI: 0.95 - 0.97). Accordingly, we could predict the risk of NAFLD using fewer variables used in relevant studies. Lin (36) announced that BMI, hemoglobin, fasting glucose and triglyceride are the major variables to predict the risk of NAFLD. However, in another study, BMI, WHR, triglyceride, glucose, systolic blood pressure and ALT were reported as the major influential factors in the used model (37). It should be noted that compared to other relevant studies, in our study, no anthropometric variables were entered into the main model. This suggests that metabolic variables, including fasting insulin level possibly play a significant role in enhancing the susceptibility to NAFLD. Anthropometric variables, such as the waist size or BMI may exert their effect through altering the level of metabolic factors, including fasting insulin level (38, 39). Analysis of mediator variables are needed to confirm this finding (26).

According to our results, the cutoff point of vitamin D level to increase the risk of NAFLD was 18.25 and 21.5 nmol/L for women and men, respectively. Using the Youden’s index, which is calculated as the maximum sum of sensitivity and characteristics, the area under the curve for the cutoff point obtained 0.865 and 0.858 for women and men, respectively.

In many studies, the level of vitamin D has been compared in patient and healthy individuals. In a meta-analysis, Eliades et al. (14) investigated the level of vitamin D in patients with NAFLD and concluded that those with NAFLD had 0.36 ng/mL (95% CI: 0.32 - 0.40) lower serum vitamin D level compared to the healthy subjects and also they were 1.26 times more likely to suffer from vitamin D deficiency (OR: 1.26, 95% CI: 1.17 - 1.35). By converting to nmol/L, the difference in vitamin D level between two groups obtained 1 nmol/L. However, the mean difference between the two groups in our study was about 12 nmol/L.

This study faced some limitations. This research was a cross-sectional study, which limited us to draw a causative conclusion. Nevertheless, the studied population was larger than most of other studies. One of the strengths of our study is using volunteer healthy people who had not a history of NAFLD, so as a population-based study, its results can be generalized to the community. Our results showed that an appropriate level of vitamin D can prevent the development of NAFLD. Moreover, the samples were selected by random sampling from those who referred to the medical center following a recall through a year. Although it minimized the possibility of selection bias, however sampling in different seasons was another limitation of this study. Using FibroScan device as a verified method to study NAFLD (18, 40), confirming the results by a gastroenterologist using FibroScan device as well as the large number of factors affecting NAFLD can enhance the reliability of our research. However, conducting a longitudinal study at multicenter level is suggested to obtain more accurate results.

5.1. Conclusions

The results of various studies and also our findings indicate a significant relationship between vitamin D levels and the risk of NAFLD. Particularly, our study suggests that the values less than 18 nmol/L for women and 21 nmol/L for men can increase the risk of NAFLD. Accordingly, the level of 25 nmol/L can be introduced as the threshold for preventing NAFLD in the community. In conclusion, it is important to maintain the appropriate serum level of vitamin D in preventing NAFLD that can be provided by consuming a vitamin D-rich food basket or by more exposure to sunlight.

Acknowledgements

References

  • 1.

    Dey P, Sutradhar S, Barman T, Khan N, Hasan I, Haque M, et al. Risk factors of non-alcoholic fatty liver disease. Mymensingh Med J. 2013;22(4):649-54.

  • 2.

    Barchetta I, Angelico F, Del Ben M, Baroni MG, Pozzilli P, Morini S, et al. Strong association between non alcoholic fatty liver disease (NAFLD) and low 25(OH) vitamin D levels in an adult population with normal serum liver enzymes. BMC Med. 2011;9:85. [PubMed ID: 21749681]. [PubMed Central ID: PMC3148980]. https://doi.org/10.1186/1741-7015-9-85.

  • 3.

    Moghaddasifar I, Lankarani KB, Moosazadeh M, Afshari M, Ghaemi A, Aliramezany M, et al. Prevalence of Non-alcoholic fatty liver disease and its related factors in Iran. Int J Organ Transplant Med. 2016;7(3):149-60. [PubMed ID: 27721961]. [PubMed Central ID: PMC5054138].

  • 4.

    Fattahi MR, Niknam R, Safarpour A, Sepehrimanesh M, Lotfi M. The prevalence of metabolic syndrome in non-alcoholic fatty liver disease; A population-based study. Middle East J Dig Dis. 2016;8(2):131-7. [PubMed ID: 27252820]. [PubMed Central ID: PMC4885612]. https://doi.org/10.15171/mejdd.2016.18.

  • 5.

    Jaruvongvanich V, Ahuja W, Sanguankeo A, Wijarnpreecha K, Upala S. Vitamin D and histologic severity of nonalcoholic fatty liver disease: A systematic review and meta-analysis. Dig Liver Dis. 2017;49(6):618-22. [PubMed ID: 28274829]. https://doi.org/10.1016/j.dld.2017.02.003.

  • 6.

    Kwok RM, Torres DM, Harrison SA. Vitamin D and nonalcoholic fatty liver disease (NAFLD): is it more than just an association? Hepatology. 2013;58(3):1166-74. [PubMed ID: 23504808]. https://doi.org/10.1002/hep.26390.

  • 7.

    Osati S, Homayounfar R, Hajifaraji M. Metabolic effects of vitamin D supplementation in vitamin D deficient patients (a double-blind clinical trial). Diabetes Metab Syndr. 2016;10(2 Suppl 1):S7-S10. [PubMed ID: 27094871]. https://doi.org/10.1016/j.dsx.2016.01.007.

  • 8.

    Barchetta I, Cimini FA, Cavallo MG. Vitamin D supplementation and non-alcoholic fatty liver disease: Present and future. Nutrients. 2017;9(9):1015. [PubMed ID: 28906453]. [PubMed Central ID: PMC5622775]. https://doi.org/10.3390/nu9091015.

  • 9.

    Arteh J, Narra S, Nair S. Prevalence of vitamin D deficiency in chronic liver disease. Dig Dis Sci. 2010;55(9):2624-8. [PubMed ID: 19960254]. https://doi.org/10.1007/s10620-009-1069-9.

  • 10.

    Manco M, Ciampalini P, Nobili V. Low levels of 25-hydroxyvitamin D(3) in children with biopsy-proven nonalcoholic fatty liver disease. Hepatology. 2010;51(6):2229. author reply 2230. [PubMed ID: 20513013]. https://doi.org/10.1002/hep.23724.

  • 11.

    Targher G, Bertolini L, Scala L, Cigolini M, Zenari L, Falezza G, et al. Associations between serum 25-hydroxyvitamin D3 concentrations and liver histology in patients with non-alcoholic fatty liver disease. Nutr Metab Cardiovasc Dis. 2007;17(7):517-24. [PubMed ID: 16928437]. https://doi.org/10.1016/j.numecd.2006.04.002.

  • 12.

    Katz K, Brar PC, Parekh N, Liu YH, Weitzman M. Suspected nonalcoholic Fatty liver disease is not associated with vitamin d status in adolescents after adjustment for obesity. J Obes. 2010;2010:496829. [PubMed ID: 21331294]. [PubMed Central ID: PMC3038678]. https://doi.org/10.1155/2010/496829.

  • 13.

    Wang X, Li W, Zhang Y, Yang Y, Qin G. Association between vitamin D and non-alcoholic fatty liver disease/non-alcoholic steatohepatitis: Results from a meta-analysis. Int J Clin Exp Med. 2015;8(10):17221-34. [PubMed ID: 26770315]. [PubMed Central ID: PMC4694215].

  • 14.

    Eliades M, Spyrou E, Agrawal N, Lazo M, Brancati FL, Potter JJ, et al. Meta-analysis: Vitamin D and non-alcoholic fatty liver disease. Aliment Pharmacol Ther. 2013;38(3):246-54. [PubMed ID: 23786213]. https://doi.org/10.1111/apt.12377.

  • 15.

    European Association for the Study of the L, European Association for the Study of D, European Association for the Study of O. EASL-EASD-EASO clinical practice guidelines for the management of non-alcoholic fatty liver disease. Obes Facts. 2016;9(2):65-90. [PubMed ID: 27055256]. [PubMed Central ID: PMC5644799]. https://doi.org/10.1159/000443344.

  • 16.

    Friedewald WT, Levy RI, Fredrickson DS. Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem. 1972;18(6):499-502. [PubMed ID: 4337382].

  • 17.

    Lips P. Which circulating level of 25-hydroxyvitamin D is appropriate? J Steroid Biochem Mol Biol. 2004;89-90(1-5):611-4. [PubMed ID: 15225848]. https://doi.org/10.1016/j.jsbmb.2004.03.040.

  • 18.

    Bedogni G, Bellentani S, Miglioli L, Masutti F, Passalacqua M, Castiglione A, et al. The Fatty Liver Index: A simple and accurate predictor of hepatic steatosis in the general population. BMC Gastroenterol. 2006;6:33. [PubMed ID: 17081293]. [PubMed Central ID: PMC1636651]. https://doi.org/10.1186/1471-230X-6-33.

  • 19.

    Mikolasevic I, Orlic L, Franjic N, Hauser G, Stimac D, Milic S. Transient elastography (FibroScan((R))) with controlled attenuation parameter in the assessment of liver steatosis and fibrosis in patients with nonalcoholic fatty liver disease - Where do we stand? World J Gastroenterol. 2016;22(32):7236-51. [PubMed ID: 27621571]. [PubMed Central ID: PMC4997649]. https://doi.org/10.3748/wjg.v22.i32.7236.

  • 20.

    Bonder A, Afdhal N. Utilization of FibroScan in clinical practice. Curr Gastroenterol Rep. 2014;16(2):372. [PubMed ID: 24452634]. https://doi.org/10.1007/s11894-014-0372-6.

  • 21.

    Cassinotto C, Boursier J, de Ledinghen V, Lebigot J, Lapuyade B, Cales P, et al. Liver stiffness in nonalcoholic fatty liver disease: A comparison of supersonic shear imaging, FibroScan, and ARFI with liver biopsy. Hepatol. 2016;63(6):1817-27. [PubMed ID: 26659452]. https://doi.org/10.1002/hep.28394.

  • 22.

    Pacifico L, Osborn JF, Bonci E, Pierimarchi P, Chiesa C. Association between Vitamin D Levels and nonalcoholic fatty liver disease: Potential confounding variables. Mini Rev Med Chem. 2019;19(4):310-32. [PubMed ID: 30360708]. https://doi.org/10.2174/1389557518666181025153712.

  • 23.

    Dongiovanni P, Anstee QM, Valenti L. Genetic predisposition in NAFLD and NASH: Impact on severity of liver disease and response to treatment. Curr Pharm Des. 2013;19(29):5219-38. [PubMed ID: 23394097]. [PubMed Central ID: PMC3850262]. https://doi.org/10.2174/13816128113199990381.

  • 24.

    Nobili V, Marcellini M, Devito R, Ciampalini P, Piemonte F, Comparcola D, et al. NAFLD in children: A prospective clinical-pathological study and effect of lifestyle advice. Hepatol. 2006;44(2):458-65. [PubMed ID: 16871574]. https://doi.org/10.1002/hep.21262.

  • 25.

    de Wit NJ, Afman LA, Mensink M, Muller M. Phenotyping the effect of diet on non-alcoholic fatty liver disease. J Hepatol. 2012;57(6):1370-3. [PubMed ID: 22796155]. https://doi.org/10.1016/j.jhep.2012.07.003.

  • 26.

    Ghaemi A, Hosseini N, Osati S, Naghizadeh MM, Dehghan A, Ehrampoush E, et al. Waist circumference is a mediator of dietary pattern in non-alcoholic fatty liver disease. Sci Rep. 2018;8(1):4788. [PubMed ID: 29555959]. [PubMed Central ID: PMC5859081]. https://doi.org/10.1038/s41598-018-23192-x.

  • 27.

    Bril F, Maximos M, Portillo-Sanchez P, Biernacki D, Lomonaco R, Subbarayan S, et al. Relationship of vitamin D with insulin resistance and disease severity in non-alcoholic steatohepatitis. J Hepatol. 2015;62(2):405-11. [PubMed ID: 25195551]. https://doi.org/10.1016/j.jhep.2014.08.040.

  • 28.

    Chun RF, Liu PT, Modlin RL, Adams JS, Hewison M. Impact of vitamin D on immune function: Lessons learned from genome-wide analysis. Front Physiol. 2014;5:151. [PubMed ID: 24795646]. [PubMed Central ID: PMC4000998]. https://doi.org/10.3389/fphys.2014.00151.

  • 29.

    Jiang YJ, Teichert AE, Fong F, Oda Y, Bikle DD. 1alpha,25(OH)2-dihydroxyvitamin D3/VDR protects the skin from UVB-induced tumor formation by interacting with the beta-catenin pathway. J Steroid Biochem Mol Biol. 2013;136:229-32. [PubMed ID: 23026511]. [PubMed Central ID: PMC4102608]. https://doi.org/10.1016/j.jsbmb.2012.09.024.

  • 30.

    Crespo J, Cayon A, Fernandez-Gil P, Hernandez-Guerra M, Mayorga M, Dominguez-Diez A, et al. Gene expression of tumor necrosis factor alpha and TNF-receptors, p55 and p75, in nonalcoholic steatohepatitis patients. Hepatology. 2001;34(6):1158-63. [PubMed ID: 11732005]. https://doi.org/10.1053/jhep.2001.29628.

  • 31.

    Cohen-Lahav M, Shany S, Tobvin D, Chaimovitz C, Douvdevani A. Vitamin D decreases NFkappaB activity by increasing IkappaBalpha levels. Nephrol Dial Transplant. 2006;21(4):889-97. [PubMed ID: 16455676]. https://doi.org/10.1093/ndt/gfi254.

  • 32.

    Adams LA, White SW, Marsh JA, Lye SJ, Connor KL, Maganga R, et al. Association between liver-specific gene polymorphisms and their expression levels with nonalcoholic fatty liver disease. Hepatol. 2013;57(2):590-600. [PubMed ID: 23213074]. https://doi.org/10.1002/hep.26184.

  • 33.

    Barchetta I, Carotti S, Labbadia G, Gentilucci UV, Muda AO, Angelico F, et al. Liver vitamin D receptor, CYP2R1, and CYP27A1 expression: Relationship with liver histology and vitamin D3 levels in patients with nonalcoholic steatohepatitis or hepatitis C virus. Hepatol. 2012;56(6):2180-7. [PubMed ID: 22753133]. https://doi.org/10.1002/hep.25930.

  • 34.

    Vaidya A, Williams JS, Forman JP. The independent association between 25-hydroxyvitamin D and adiponectin and its relation with BMI in two large cohorts: The NHS and the HPFS. Obesity (Silver Spring). 2012;20(1):186-91. [PubMed ID: 21760630]. [PubMed Central ID: PMC3461263]. https://doi.org/10.1038/oby.2011.210.

  • 35.

    Hui JM, Hodge A, Farrell GC, Kench JG, Kriketos A, George J. Beyond insulin resistance in NASH: TNF-alpha or adiponectin? Hepatol. 2004;40(1):46-54. [PubMed ID: 15239085]. https://doi.org/10.1002/hep.20280.

  • 36.

    Lin YJ, Gao XM, Pan WW, Gao S, Yu ZZ, Xu P, et al. A model to predict the onset of non-alcoholic fatty liver disease within 2 years in elderly adults. J Gastroenterol Hepatol. 2017;32(10):1739-45. [PubMed ID: 28183156]. https://doi.org/10.1111/jgh.13760.

  • 37.

    Birjandi M, Ayatollahi SM, Pourahmad S, Safarpour AR. Prediction and diagnosis of non-alcoholic fatty liver disease (NAFLD) and identification of Its associated factors using the classification tree method. Iran Red Crescent Med J. 2016;18(11). e32858. [PubMed ID: 28191344]. [PubMed Central ID: PMC5292777]. https://doi.org/10.5812/ircmj.32858.

  • 38.

    Babai MA, Arasteh P, Hadibarhaghtalab M, Naghizadeh MM, Salehi A, Askari A, et al. Defining a BMI cut-off point for the Iranian population: The Shiraz heart study. PLoS One. 2016;11(8). e0160639. [PubMed ID: 27509026]. [PubMed Central ID: PMC4980035]. https://doi.org/10.1371/journal.pone.0160639.

  • 39.

    Ehrampoush E, Arasteh P, Homayounfar R, Cheraghpour M, Alipour M, Naghizadeh MM, et al. New anthropometric indices or old ones: Which is the better predictor of body fat? Diabetes Metab Syndr. 2017;11(4):257-63. [PubMed ID: 27578617]. https://doi.org/10.1016/j.dsx.2016.08.027.

  • 40.

    Alizadeh A, Mansour-Ghanaei F, Roozdar A, Joukar F, Sepehrimanesh M, Hojati SA, et al. Laboratory tests, liver vessels color doppler sonography, and FibroScan findings in patients with nonalcoholic fatty liver disease: An observation study. J Clin Imaging Sci. 2018;8:12. [PubMed ID: 29692949]. [PubMed Central ID: PMC5894278]. https://doi.org/10.4103/jcis.JCIS_93_17.