|Year : | Volume
| Issue : | Page :
|Serum 25-hydroxyvitamin D levels and the risk of non-alcoholic fatty liver: A two-sample Mendelian randomization study
Qi Sheng1, Huanchen Shi2, Shousheng Liu3, Likun Zhuang3, Zhenzhen Zhao3, Yongning Xin4
1 Department of Infectious Disease, Qingdao Municipal Hospital, Cheeloo College of Medicine, Shandong University, Jinan; Department of Nutrition, The Affiliated Hospital of Qingdao University, Qingdao, Shangdong, China
2 School of Basic Medicine, College of Medicine, Qingdao University, Qingdao, Shangdong, China
3 Clinical Research Center, Qingdao, China
4 Department of Infectious Disease, Qingdao Municipal Hospital, Cheeloo College of Medicine, Shandong University, Jinan, Shangdong; Department of Infectious Disease, Qingdao Municipal Hospital, Qingdao University, Qingdao, China
Click here for correspondence address and email
|Date of Submission||04-Jul-2022|
|Date of Acceptance||08-Sep-2022|
|Date of Web Publication||10-Oct-2022|
| Abstract|| |
Background: Accumulated studies have shown that low expression of 25-hydroxyvitamin D [25(OH)D] was significantly associated with the risk of non-alcoholic fatty liver disease (NAFLD). However, the exact causality is still unknown. The aim of this study was to investigate whether levels of 25(OH)D are associated with risk of NAFLD, using a two-sample Mendelian randomization (MR).
Methods: Data from a recent large vitamin D genome-wide association study (GWAS) on 417,580 Europeans were utilized, and the largest published histology-based NAFLD GWAS study (1,483 cases and 17,781 healthy controls) for genetic variants predicted to cause NAFLD were searched. All genetic datasets for the MR analyses were obtained using publicly available summary statistics based on individuals of European ancestry from the MR-Base and NHGRI-EBI GWAS Catalog database. Inverse-variance weighted (IVW) MR approach was used to estimate causal effects in the main analysis, complemented by 4 additional methods to control for pleiotropy. Sensitivity analyses were conducted to verify whether heterogeneity and pleiotropy can bias the MR results.
Results: The MR analysis did not provide strong evidence for the causal association of circulating 25(OH)D with NAFLD by IVW method (OR = 0.746, 95%CI 0.517–1.078; P = 0.119). The results were consistent using four other MR methods. Sensitivity analysis using all different analytical approaches yielded similar results. There was no evidence for pleiotropy (MR-Egger intercept: −0.0003758, P = 0.970). The replication process also showed consistent results using IVW method (P = 0.710).
Conclusion: This study indicates that serum 25(OH)D levels did not possess an obvious effect on the risk of NAFLD. The associations in previous studies may be due to residual confounding or reverse causation.
Keywords: 25-hydroxyvitamin D, GWAS, Mendelian randomization, NAFLD, SNPs
|How to cite this URL:|
Sheng Q, Shi H, Liu S, Zhuang L, Zhao Z, Xin Y. Serum 25-hydroxyvitamin D levels and the risk of non-alcoholic fatty liver: A two-sample Mendelian randomization study. Saudi J Gastroenterol [Epub ahead of print] [cited 2022 Dec 6]. Available from: https://www.saudijgastro.com/preprintarticle.asp?id=358009
Authors Qi Sheng and Huanchen Shi contributed equally to this work.
| Introduction|| |
Non-alcoholic fatty liver disease (NAFLD) represents the most common cause of chronic liver disease and affects almost a quarter of the world's population., It has become an epidemic, which causes a tremendous health and economic burden. Moreover, NAFLD can evolve towards steatohepatitis (NASH), eventually cirrhosis and hepatocellular carcinoma. There is still no Food and Drug Administration (FDA) approved medications available for the treatment of NAFLD. Therefore, it is particularly urgent to identify-risk factors of this disease.
Vitamin D is a sterol derivative that affects a few common chronic diseases such as bone metabolic disorders, tumors, cardiovascular diseases, diabetes and so on. Serum 25-hydroxyvitamin D [25(OH)D] is a primary circulating form of vitamin D, which can best reflect vitamin D stores., In the last few years, several studies explored the association between vitamin D and NAFLD but arrived at controversial conclusions. An early meta-analysis linked the risk for NAFLD to lower circulating concentrations of 25(OH)D. Two recent meta-analyses concluded that 25(OH)D level was not related to either NAFLD activity score or fibrosis., Whether vitamin D deficiency is a contributing factor to NAFLD remains unclear. However, observational studies and RCTs are limited to confounding factors, such as obesity, which affects both serum 25(OH)D concentrations and incidence of NAFLD.,, Further, it is difficult to judge the sequence of exposure and outcome which may lead to the reverse causality.
Mendelian randomization (MR) has been increasingly used to assess causal inference between exposure and outcome. The principle is based on the surmise that alleles are randomly allocated at meiosis and can be used as instrumental variables (IVs) to test causation. This method can strengthen the exposure-outcome association inference by robustly reducing the risk of confounding and eliminating reverse causality. Two-sample MR is one particularly valuable MR method using two independent GWAS datasets. Through dividing the dataset into two (or more) samples for estimation, it can mitigate the phenomenon of the winner's curse leading to bias in MR.,
Mendelian randomization has been used elsewhere to investigate the effect of circulating 25(OH)D levels on NAFLD risk. However, the study was limited by the fact that all participants were of Asian origin and the diagnosis was based on ultrasound. In this study, we examined whether 25(OH)D is causally related to NAFLD risk using a two-sample MR approach, with recent GWAS studies of the Europeans population. Moreover, the diagnosis of NAFLD was histology-based and more IVs were extracted from recent GWAS datasets. Several sensitivity analyses were performed to critically examine the underlying assumptions of the MR approach. We also performed MR with another two GWAS datasets for replication purposes to make our results more robust.
| Materials and Methods|| |
Software selection and data source
This two-sample MR was conducted using the “Two Sample MR” R package following the guidelines provided by the developers (https://mrcieu.github.io/TwoSampleMR), and in-house developed R scripts (available upon request). The model based on MR is shown in [Figure 1] and needs to satisfy three assumptions [Figure 2]. In order to avoid the potential bias caused by population stratification, only the subjects with genetic background of European origin were selected. GWAS data of 25(OH)D from the Integrative Epidemiology Unit (IEU) GWAS database based on the large UK Biobank (UKB) sample (N = 417,580) were obtained. The datasets of Revez et al.'s study was obtained through observing data from UK Biobank. The database was developed by the MRC IEU at the University of Bristol. It is a collection of complete GWAS summary datasets, which can be downloaded as open-source files. The NAFLD dataset was obtained from the largest histology-based NAFLD GWAS which involved 1,483 European NAFLD cases and 17,781 genetically matched controls. This dataset was obtained from NHGRI-EBI GWAS Catalog database. Another large GWAS study associated with 25(OH)D (including 79,366 European-ancestry individuals) and GWAS open summary level datasets associated with NAFLD from FinnGen study were included for replication purposes. All datasets were publicly available summarized results from published genome-wide association studies.
|Figure 1: Forest plot of the MR sensitivity analyses excluding SNPs with possible pleiotropic effects|
Click here to view
|Figure 2: Graph illustrating Mendelian randomization assumptions. (1) SNPs are strongly correlated with exposure. (2) SNPs should only affect the outcomes through affecting 25(OH)D. (3) SNPs should not be associated with potential confounder|
Click here to view
Selection of the genetic instruments
Instrumental variables associated with 25(OH)D concentration were extracted from IEU GWAS database as previously mentioned, including 417,580 individuals of European ancestry. One hundred and forty three loci associated with 25(OH)D concentration were identified. Instrument SNPs associated with 25(OH)D levels were chosen at a level of genome-wide significance (P < 5 × 10−8), so as to satisfy the first MR assumption that the instrument SNPs robustly associates with the exposure. Independent SNPs were selected with the clump function in the “Two Sample MR” R package, which used the 1,000 genomes project as a reference panel. Clumping process (R2 < 0.001, window size = 10,000 kb) was conducted with the samples to ensure the instruments were truly in linkage equilibrium. PhenoScanner database was explored to exclude SNPs with traits inducing potential horizontal pleiotropic associations. Genetic variants of palindromic and incompatible alleles should be removed when harmonizing exposure and outcome. 25(OH)D associated variants which possess the potential effects on NAFLD were extracted and the Proxy SNPs were excluded. F-statistics was calculated to examine the strength of the instruments using an equation used in previous MR studies. If the F-statistic is greater than 10, the weak instrument bias appears unlikely., The variance (R2) of 25(OH)D, which is explained by the instruments, was calculated using the approximation described by Bowden et al., and the calculation method was showed in the Supplementary Material Text S1.
Mendelian randomization analysis
A two-sample MR was performed for testing the causal relationship between 25(OH)D and NAFLD. Analyses were performed using Two Sample MR R package of MR-Base (https://github.com/MRCIEU/TwoSampleMR). Causal effect was estimated in the primary analyses with the inverse-variance weighted (IVW) method, which uses a meta-analysis approach to combine the Wald ratios of the causal effects of each SNP and can provide the most precise estimates., When the directional pleiotropy (P for MR-Egger intercept >0.05) was not found among the selected IVs, the IVW method was regarded as reliable. MR-Egger, weighted median, weighted mode methods and simple mode analysis were also performed to estimate the causal effect of exposure on outcome to make our results more robust. Furthermore, the PhenoScanner tool was used to search for other phenotypes (including BMI, alcohol, triglyceride, cholesterol, liver enzyme, waist circumference, LDL) which are related to the selected instrumental variable SNPs at risk of affecting NAFLD phenotypes. Then these SNPs were manually removed from the MR analysis to rule out possible pleiotropic effects.
Heterogeneity indicates potential violation of the IV assumptions and displayed the results using forest plots. Cochran's Q statistics was used to estimate the heterogeneity of IVW. Cochran's Q test follows χ2 distribution, and P < 0.01 is defined as significant heterogeneity.
To evaluate the sensitivity of IVs in MR analysis, a leave-one-out analysis was performed by discarding each SNP in turn and repeating the IVW analysis to examine whether a single SNP could drive the association. In the main MR analysis, sensitivity analysis was estimated by excluding variants associated with potential confounders. Therefore, the PhenoScanner database was used to search GWAS traits of each 25(OH)D-related SNP as instrument, in order to identify the potential confounders. When the GWAS P value of the variant for a trait was below 5 × 10−8, the associations were considered as positive. Certain traits were classified into 10 larger trait categories and an unclassified group, according to the number of traits-associated SNPs [Table 1]. Although no confounders were found which may cause the horizontal pleiotropy between the vitamin D and NAFLD, we evaluated parameters that may be associated with NAFLD in sensitivity analysis. Through searching by keywords in PhenoScanner database, we identified all the SNPs of possible traits that may associate with NAFLD. The traits included lifestyle (alcohol consumption), obesity-related (waist circumference, body mass index, etc.), lipid-related (triglyceride, cholesterol, etc.), and so on. Then traits with SNPs intersecting with 25(OH)D were obtained to exclude the interference [Table 1]. Certain parameters were classified into 10 larger trait categories and an unclassified group. The top three trait categories which were ranged according to the relevancy of phenotype and SNPs were selected for sensitivity analysis. Blood cell counts, body composition and anthropometrics were included [Table 1]. In addition, a sensitivity analysis was performed by excluding all SNPs that were associated with any trait with P value below 5 × 10−8 in the PhenoScanner database. The aim of performing special sensitivity analyses by excluding traits clustered SNPs is to try to explore how these SNPs impact on 25(OH)D levels and distinguish vertical pleiotropy to adjust causal estimation.
|Table 1: Results of the MR study testing causal association between low 25OHD and NAFLD|
Click here to view
| Results|| |
Selection of genetic instrumental variables
Following the method described above, 113 significant SNPs (P < 5 × 10−8) were selected accounting for approximately 3.5% of vitamin D variation, while 10 SNPs were absent in the NAFLD GWAS datasets and 26 SNPs with NAFLD-associated traits were excluded based on PhenoScanner database [Supplementary Materials Table S1[Additional file 1]]. We also conducted another harmonizing process and excluded 10 SNPs as a result of palindromic sequences. Finally, 63 SNPs were included after removing linkage disequilibrium and palindrome sequence. Detailed information of the 63 IVs is available in the supplementary material. In replication practice, the data showed six SNPs associated with genome-wide significance (P < 5 × 10–8) with 25(OH)D. One SNP (rs12785878) did not exist in NAFLD GWAS datasets. A palindromic SNP (rs8018720) was not included. Finally, four intersected SNPs (rs10741657, rs10745742, rs17216707, rs3755967) were included in the analysis [Supplementary Material Table S2[Additional file 2]].
The influence of 25(OH)D on the risk of non-alcoholic fatty liver disease
The results of the MR analyses are shown in [Table 1] and [Figure 1]. No significant association of low vitamin D levels with increased risk of NAFLD was identified (OR = 0.746, 95%CI 0.517–1.078; P = 0.119). A potential heterogeneity was identified using the Cochran's Q test (Q = 84.9; P = 0.0283). IVW with multiplicative random effects method was utilized in case of potential heterogeneity. The total F-statistic was 198.82 and each single F statistic was larger than 10, indicating strong instrumental variables. There was no evidence for pleiotropy (MR-Egger intercept: −0.0003758, se = 0.010079; P = 0.97). Similar results were obtained using the other 4 MR methods [Figure 1]. A scatter plot was also generated [Supplementary Material Figure S1[Additional file 4]].
In replication practice, the result also showed a negative causal link between vitamin D and NAFLD using the IVW method (OR = 0.609, 95%CI: 0.045–8.307; P = 0.710), simple mode method (OR = 0.905, 95%CI: 0.0393–20.844; P = 0.954), weighted median method (OR = 1.148, 95%CI: 0.121–10.925; P = 0.954) and MR-Egger method (OR = 7.498, 95%CI: 0.201–280.004; P = 0.389). No significantly causal relationship between vitamin D and NAFLD was found under the different MR methods [Supplementary Material Figures S2[Additional file 5]] and [Supplementary Material Figures S3[Additional file 6]].
Although the intercept obtained from the MR-Egger regression did not show significant unbalanced horizontal pleiotropy (MR-Egger intercept: −0.0003758, se = 0.010079; P = 0.97), evidence for pleiotropy based on MR-PRESSO analysis (P value global test 0.001) was found. The result of MR estimates did not alter the inference after removing outlier SNPs [Supplementary Materials Table S3[Additional file 3]]. Leave-one-out sensitivity analysis by discarding each SNP one by one was conducted and the results remained unchanged (P > 0.05) [Supplementary Materials Figure S4[Additional file 7]]. Weighted Median, Simple Mode, Weighted Mode and MR-Egger method confirmed the lack of association [Supplementary Material Figure S5[Additional file 8]]. All the results indicated that none of the genetic variants were sensitive.
Additionally, sensitivity analyses were conducted after removing several traits-clustered SNPs, as the above methods described, referring to previous studies. SNPs associated with blood cell counts, body composition and anthropometrics were included in this sensitivity analysis. Moreover, SNPs with any distinctive traits based on PhenoScanner databases were removed (P < 5 × 10−8). Sensitivity analyses demonstrated a consistent association, suggesting that our conclusions were robust [Figure 1].
In replication practice, leave-one-out sensitivity analysis by discarding each SNP in turn was conducted and the results were also consistent, [Supplementary Material Figure S6[Additional file 9]].
| Discussion|| |
In this study, the Two-sample MR was used to evaluate the association between 25(OH)D and NAFLD. The results showed that there was no genetic evidence for a causal effect of 25(OH)D levels on the risk of NAFLD, and this conclusion was consistent after the sensitivity analyses and replication. To our knowledge, this is the first two-sample MR study to examine the effect of 25(OH)D levels on histology-based NAFLD risk on a European population.
There has been a significant scientific interest in the relationship between vitamin D status and NAFLD. However, whether the serum level of vitamin D is a risk factor for NAFLD is controversial, where previous findings have been contradictory., Many previous studies have investigated the association between vitamin D deficiency and NAFLD/NASH. Observational studies found that vitamin D levels decreased in NAFLD.,, Whether vitamin D deficiency is a contributing factor to NAFLD remains unclear. Some studies showed that low vitamin D levels were associated with the presence of NAFLD, suggesting that vitamin D may be a target for the treatment of NAFLD and NASH.,,,37],, In a meta-analysis, patients in NAFLD were 26% more likely to have lower 25(OH)D level, reflecting deficiency in vitamin D, when compared with controls. However, the results of other studies were inconsistent.,,,,
Differences between conclusions may be attributed to a number of reasons. Initially, the failure in controlling the confounders included different sample sizes, races, BMI, lack of a control group without NAFLD on liver biopsy and so on. Some confounders themselves such as BMI can apparently influence the 25(OH)D serum levels, making it difficult to distinguish whether the confounders or the 25(OH)D serum levels caused the NAFLD. On the other hand, observational studies can be hampered by confounding or reverse causation. NAFLD contributes to lower concentrations of serum 225(OH)D. The limitation of current works prevent us from establishing any temporality.
MR analyses reduce confounding and reverse causality due to the random allocation of genotypes from parents to offspring. To the best of our knowledge, the present MR study is the first two-sample MR to evaluate the causal role of 25(OH)D for NAFLD in the European population. We have used the largest GWAS of serum 25(OH)D traits and diagnosed NAFLD disease to investigate causal relationships. Horizontal pleiotropy means that genetic variation has many independent ways to affect the results and reject the IV assumptions of the MR model. In this study, MR-Egger test showed the low intercept and high P value (−0.00038; P = 0.970). It proved that the results were un-affected by horizontal multidirectionality.
There has been a MR analysis limited to less genetic variation tools sample incorporating four vitamin D-related SNPs (rs12785878, rs10741657, rs2282679, rs6013897), which were chosen on the basis of a genome-wide association study on 25(OH)D. All participants were of Asian origin and the diagnosis of NAFLD was made by ultrasound. These factors may bias the results. This Mendelian randomization study used the largest histology-based NAFLD GWAS datasets as outcome data. GWAS data of 25(OH)D from the IEU GWAS database based on the large UKB sample (N = 417,580) were obtained. The research chose a conservative threshold of 25 nmol/L for vitamin D deficiency and applied a rank-based inverse-normal transformation (RINT) to the phenotype (25(OH)D levels) and fit age at time of assessment, sex, assessment month, assessment centre, supplement-intake information, genotyping batch and the first 40 ancestry PCs as covariates in the model. The above aspects made the datasets more reliable and authoritative. The results of our Mendelian analysis were consistant and reassuring by using five MR methods (IVW, MR-Egger, weighted median, weighted mode MR and simple mode analysis) and sensitivity analyses by excluding SNPs with pleiotropic effects. In addition, another two GWAS samples were used for the replication to validate the conclusion.
This study has several limitations. Primarily, the exact function of some SNPs remains unknown, making residual bias possible while examining pleiotropy. However, we obtained consistent results using 5 MR methods making sure of pleiotropy, robust methodology, and sensitivity analyses excluding SNPs with pleiotropic effects is reassuring. Secondarily, the MR results cannot be generalized to other race origins. Although both exposure and outcome GWAS were restricted to European ancestry participants, residual confounding from population stratification cannot be completely excluded, due to ethnic differences in the exposure and outcome GWAS populations.
In summary, this study provides new evidence that serum 25(OH)D levels may not be causally related to the risk of NAFLD. Our findings imply that the associations based on some previous observational studies between 25(OH)D and risk of NAFLD might be due to reverse causality and confounders, such as environment correlated with exposure to sun and skin pigmentation. This assumption cannot be eliminated assuming that low concentration of vitamin D may just play an intermediary role as a real risk factor for NAFLD. This conclusion needs to be confirmed by large RCTs or future MR evidence based on updated GWAS samples.
Financial support and sponsorship
This project was supported by the National Natural Science Foundation of China (Grant No. 32171277).
Conflicts of interest
There are no conflicts of interest.
| Supplementary Material|| |
Text S1: Equation For F-statistic
The total F statistic of the IVs sample in the exposure was explained by the genetic variants (R2), sample size (N) and number of instruments (K), by the formula:
We characterized instrument strength in the univariable MR setting by generating mean F-statistic. For a single variant, we used the approximation described by Bowden et al :
Where γj is the SNP-exposure association and σxj is the standard deviation for SNP-exposure association for variant j.
- Bowden, J., Del Greco, M.F., Minelli, C., Davey Smith, G., Sheehan, N.A., and Thompson, J.R. (2016). Assessing the suitability of summary data for two-sample Mendelian randomization analyses using MR-Egger regression: the role of the I2 statistic. Int J Epidemiol 45, 1961-1974.
| References|| |
Younossi Z, Anstee QM, Marietti M, Hardy T, Henry L, Eslam M, et al.
Global burden of NAFLD and NASH: trends, predictions, risk factors and prevention. Nat Rev Gastroenterol Hepatol 2018;15:11-20.
Younossi ZM, Koenig AB, Abdelatif D, Fazel Y, Henry L, Wymer M, Global epidemiology of nonalcoholic fatty liver disease-Meta-analytic assessment of prevalence, incidence, and outcomes. Hepatology 2016;64:73-84.
Cholankeril G, Kanwal F. Nonalcoholic fatty liver disease and hepatocellular carcinoma: Time to bridge the gap. Hepatology 2021;74:2336-8.
Wang H, Chen W, Li D, Yin X, Zhang X, Olsen N, et al.
Vitamin D and chronic diseases. Aging Dis 2017;8:346-53.
Holick MF, Chen TC. Vitamin D deficiency: A worldwide problem with health consequences. Am J Clin Nutr 2008;87:1080s-6s.
Jiang X, O'Reilly PF, Aschard H, Hsu YH, Richards JB, Dupuis J, et al.
Genome-wide association study in 79,366 European-ancestry individuals informs the genetic architecture of 25-hydroxyvitamin D levels. Nat Commun 2018;9:260.
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:246-54.
Saberi B, Dadabhai AS, Nanavati J, Wang L, Shinohara RT, Mullin GE. Vitamin D levels do not predict the stage of hepatic fibrosis in patients with non-alcoholic fatty liver disease: A PRISMA compliant systematic review and meta-analysis of pooled data. World J Hepatol 2018;10:142-54.
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:618-22.
Pereira-Santos M, Costa PR, Assis AM, Santos CA, Santos DB. Obesity and vitamin D deficiency: A systematic review and meta-analysis. Obes Rev 2015;16:341-9.
Polyzos SA, Kountouras J, Mantzoros CS. Obesity and nonalcoholic fatty liver disease: From pathophysiology to therapeutics. Metabolism 2019;92:82-97.
Fan JG, Kim SU, Wong VW. New trends on obesity and NAFLD in Asia. J Hepatol 2017;67:862-73.
Lawlor DA, Harbord RM, Sterne JA, Timpson N, Davey Smith G. Mendelian randomization: Using genes as instruments for making causal inferences in epidemiology. Stat Med 2008;27:1133-63.
Emdin CA, Khera AV, Kathiresan S. Mendelian randomization. JAMA 2017;318:1925-6.
Richmond RC, Davey Smith G, Ness AR, den Hoed M, McMahon G, Timpson NJ. Assessing causality in the association between child adiposity and physical activity levels: A Mendelian randomization analysis. PLoS Med 2014;11:e1001618.
Davey Smith G, Hemani G. Mendelian randomization: Genetic anchors for causal inference in epidemiological studies. Hum Mol Genet 2014;23:R89-98.
Wang N, Chen C, Zhao L, Chen Y, Han B, Xia F, et al.
Vitamin D and nonalcoholic fatty liver disease: Bi-directional mendelian randomization analysis. EBioMedicine 2018;28:187-93.
Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al.
The MR-Base platform supports systematic causal inference across the human phenome. Elife 2018;7:e34408. doi: 10.7554/eLife. 34408.
Revez JA, Lin T, Qiao Z, Xue A, Holtz Y, Zhu Z, et al.
Genome-wide association study identifies 143 loci associated with 25 hydroxyvitamin D concentration. Nat Commun 2020;11:1647.
Anstee QM, Darlay R, Cockell S, Meroni M, Govaere O, Tiniakos D, et al.
Genome-wide association study of non-alcoholic fatty liver and steatohepatitis in a histologically characterised cohort. J Hepatol 2020;73:505-15.
Collins A, Allelic association: Linkage disequilibrium structure and gene mapping. Mol Biotechnol 2009;41:83-9.
Yarmolinsky J, Bonilla C, Haycock PC, Langdon RJ, Lotta LA, Langenberg C, et al.
Circulating selenium and prostate cancer risk: A mendelian randomization analysis. J Natl Cancer Inst 2018;110:1035-8.
Davies NM, Holmes MV, Davey Smith G. Reading Mendelian randomisation studies: A guide, glossary, and checklist for clinicians. BMJ 2018;362:k601. doi: 10.1136/bmj.k601.
Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan NA, Thompson JR. Assessing the suitability of summary data for two-sample Mendelian randomization analyses using MR-Egger regression: The role of the I2 statistic. Int J Epidemiol 2016;45:1961-74.
Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol 2013;37:658-65.
Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan N, Thompson J. A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat Med 2017;36:1783-802.
Qiu S, Cao P, Guo Y, Lu H, Hu Y. Exploring the causality between hypothyroidism and non-alcoholic fatty liver: A mendelian randomization study. Front Cell Dev Biol 2021;9:643582. doi: 10.3389/fcell. 2021.643582.
Manousaki D, Harroud A, Mitchell RE, Ross S, Forgetta V, Timpson NJ, et al.
Vitamin D levels and risk of type 1 diabetes: A Mendelian randomization study. PLoS Med 2021;18:e1003536.
Kwok RM, Torres DM, Harrison SA. Vitamin D and nonalcoholic fatty liver disease (NAFLD): Is it more than just an association? Hepatology 2013;58:1166-74.
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:310-32.
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:517-24.
Nseir W, Taha H, Khateeb J, Grosovski M, Assy N. Fatty liver is associated with recurrent bacterial infections independent of metabolic syndrome. Dig Dis Sci 2011;56:3328-34.
Mohamed Ahmed A, Abdel Ghany M, Abdel Hakeem GL, Kamal A, Khattab R, Abdalla A, et al.
Assessment of Vitamin D status in a group of Egyptian children with non alcoholic fatty liver disease (multicenter study). Nutr Metab (Lond) 2016;13:53.
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.
Black LJ, Jacoby P, She Ping-Delfos WC, Mori TA, Beilin LJ, Olynyk JK, et al.
Low serum 25-hydroxyvitamin D concentrations associate with non-alcoholic fatty liver disease in adolescents independent of adiposity. J Gastroenterol Hepatol 2014;29:1215-22.
Hao YP, Ma XJ, Luo YQ, Ni J, Dou JX, Hu YQ, et al.
Serum vitamin D is associated with non-alcoholic fatty liver disease in Chinese males with normal weight and liver enzymes. Acta Pharmacol Sin 2014;35:1150-6.
Lu Z, Pan X, Hu Y, Hao Y, Luo Y, Hu X, et al.
Serum vitamin D levels are inversely related with non-alcoholic fatty liver disease independent of visceral obesity in Chinese postmenopausal women. Clin Exp Pharmacol Physiol 2015;42:139-45.
Zhai HL, Wang NJ, Han B, Li Q, Chen Y, Zhu CF, et al.
Low vitamin D levels and non-alcoholic fatty liver disease, evidence for their independent association in men in East China: A cross-sectional study (Survey on Prevalence in East China for Metabolic Diseases and Risk Factors (SPECT-China)). Br J Nutr 2016;115:1352-9.
Chung GE, Kim D, Kwak MS, Yang JI, Yim JY, Lim SH, et al.
The serum vitamin D level is inversely correlated with nonalcoholic fatty liver disease. Clin Mol Hepatol 2016;22:146-51.
Jeong DW, Lee HW, Cho YH, Yi DW, Lee SY, Son SM, et al.
Comparison of serum ferritin and vitamin d in association with the severity of nonalcoholic Fatty liver disease in korean adults. Endocrinol Metab (Seoul) 2014;29:479-88.
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.
Hourigan SK, Abrams S, Yates K, Pfeifer K, Torbenson M, Murray K, et al.
Relation between vitamin D status and nonalcoholic fatty liver disease in children. J Pediatr Gastroenterol Nutr 2015;60:396-404.
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:405-11.
Patel YA, Henao R, Moylan CA, Guy CD, Piercy DL, Diehl AM, et al.
Vitamin D is not associated with severity in NAFLD: Results of a paired clinical and gene expression profile analysis. Am J Gastroenterol 2016;111:1591-8.
El Khoudary SR, Samargandy S, Zeb I, Foster T, de Boer IH, Li D, et al.
Serum 25-hydroxyvitamin-D and nonalcoholic fatty liver disease: Does race/ethnicity matter? Findings from the MESA cohort. Nutr Metab Cardiovasc Dis 2020;30:114-22.
Andrews SJ, Goate A, Anstey KJ. Association between alcohol consumption and Alzheimer's disease: A Mendelian randomization study. Alzheimers Dement 2020;16:345-53.
Hemani G, Bowden J, Davey Smith G. Evaluating the potential role of pleiotropy in Mendelian randomization studies. Hum Mol Genet 2018;27:R195-208.
Department of Infectious Disease, Qingdao Municipal Hospital, 1 Jiaozhou Road, Qingdao 266011, Shandong Province
Source of Support: None, Conflict of Interest: None
[Figure 1], [Figure 2]
| Article Access Statistics|
| Viewed||1086 |
| PDF Downloaded||10 |