Association of FokI, TaqI, BsmI and ApaI polymorphisms with diabetic retinopathy: a pooled analysis of case-control studies

Background To assess synthetically the association between polymorphisms in the vitamin D receptor (VDR) gene (FokI, BsmI, ApaI, and TaqI) and diabetic retinopathy (DR). Methods Pubmed, Embase, ISI Web of Science, Google-scholar and CBMDisc, CNKI and Chongqing VIP databases were searched. A meta-analysis was performed. Results Six studies with 636 cases and 1,035 controls were included in this meta-analysis. The outcomes showed that the FokI polymorphism (F allele) of VDR gene had no statistical protective relationship with DR in overall studies. Interestingly, stratification analysis showed that the FokI polymorphism (Fallele) was significantly associated with decreased DR risk in the Chinese population, among included studies without publication bias, during a comparison analysis between normal subjects and DR patients, and among articles published after 2010. However, the TaqI, BsmI and ApaI polymorphisms of VDR gene had no significant association with the risk of DR. Conclusion This meta-analysis of case-control studies revealed that the VDR-FokI polymorphism (F allele) decreased the risk of DR in Chinese people, among included studies without publication bias, during a comparison analysis between normal subjects and DR patients, and among articles published after 2010. Further rigorous and prospective studies with large sample size are needed to confirm our findings.


Introduction
To date, diabetic retinopathy (DR) is one of most common causes of visual impairment (VI) in adults starting from the age of 20 years until 74 years 1 . According to the reports by the World Health Organization (WHO), there were 4.8% of the 37 million annual cases of blindness by Although many strategies are performed to prevent DR 4 , DR is more prevalent with the increasing prevalence of diabetes mellitus (DM), especially in China 5 . The mechanisms of DR are sophisticated and there are many studies to identify a number of potential susceptibility genes for DR including those coding for vascular endothelial growth factor (VEGF) 6 , angiotensin converting enzyme (ACE) 7 , aldose reductase (ALR) 8 and receptor for advanced glycation end products (RAGE) 9 .
Furthermore, many studies have shown that Vitamin D plays an important role in calcium metabolism, anti-angiogenic activities and inhibition of proliferation [10][11][12] . The vitamin D receptor (VDR) is a ligand-dependent transcription factor that in complex with hormonally active Vitamin D, and extensively express in vascular endothelial cells of humans retina 13 . The VDR gene is located on chromosome 12q13-12q14, and it has four frequently single nucleotide polymorphisms (SNPs). These are FokI G/A (rs2228570) change in exon 2, TaqI T/C (rs731236) change in exon 9, BsmI A/G (rs1544410), and ApaI G/T (rs7975232) changes in intron 8. Generally, the FokI polymorphism is located in the start codon, but the BsmI, ApaI and TaqI polymorphisms are located in the 3` UTR regions. Therefore, there is an alternative, a transcription initiation site, resulting in a VDR protein with three additional amino acids. A Previous study revealed a novel association between the functional VDR gene (FokI) polymorphism and a risk of DR in patients with Type 1 diabetes (T1DM) 14 . However, another study suggested that FokI polymorphism of VDR gene did not increase risk of DR in Type 2 diabetes (T2DM) in Poland 15 . Hence, previous studies concerning the association between the VDR-FokI gene polymorphism and the risk of DR are limited and rather conflicting. In order to assess the powerful-strength of the association between the VDR gene polymorphism and susceptibility to DR, we conducted this meta-analysis, which would obtain a greater possibility of acquiring or finding reasonably conclusions.

Search strategy
Two independent reviewers (WJY and LL) screened all studies on the relationship between VDR gene polymorphism and DR. Disagreements were dealt with by discussion. If the results of the discussion still had no agreement, another researcher (SY) would make the fi-nal decision. Both English-language databases including Pubmed, Embase, ISI Web of Science, Google (scholar), and Chinese-language databases including China National Knowledge Infrastructure (CNKI), Chinese Biological Medicine Disc (CBMDisc), Wan Fang (Chinese), Chongqing VIP database were searched using the terms "Vitamin D receptor or VD receptor or VDR", "diabetic retinopathy or DR" and "polymorphism, variant or mutation". The last search was updated on June 10, 2017. Researchers checked reference lists manually to identify additional literatures. Searching was performed without any limits of language, but only English or Chinese articles included.
Inclusion and exclusion criteria of the study Studies were included based on the following criteria: (1) case-control studies; (2) complete information for VDR gene mutations with the risk of DR for calculating the odds ratio (OR) with 95% CI directly and indirectly. The exclusion criteria were as follows: (1) editorials, abstracts, case reports, reviews and studies that used an animal or cell line, (2) any data that was insufficient or overlapping, (3) genotypes not in Hardy Weinberg equilibrium (HWE).

Data extraction and quality assessment
After selection process, data was extracted within a standard form by those two investigators. The following information was extracted: name of the First author, location, publication year, ethnicity, sample sizes and genotype distributions of cases and controls. Newcastle-Ottawa Quality Assessment Scale (NOS) 16 (Table 1) was used to evaluate the quality of included studies. The score of study with more than 6 stars was considered of relatively high quality. This meta-analysis of case-control studies was guided and performed by the PRISMA statement 17 and MOOSE guidelines 18 .

Statistical analysis
The HWE was conducted for testing the distribution of genotypes in controls within each included study, and a P-value > 0.05 suggested that the study sample was representative of the population. The ORs with 95% CI for allele and genotype models (homozygous model, heterozygous model, recessive model, and dominant model) were used to evaluate the association between the VDR gene polymorphisms and the risk of DR. I2 value was used to examine the Heterogeneity among eligible studies. When the heterogeneity was significant (I2 > 50%), the random-effects model (based on the DerSimonian and Laird method) would be used to evaluate the pooled ORs, if not, pooled ORs could then be examined by the fixed-effects model (based on the Mantel-Haenszel method). The potential publication bias was assessed with a funnel plot and Egger regression asymmetry test. An asymmetric plot suggested a possible publication bias and the P value of Egger's test less than 0.05 was considered representative of statistically significant publication bias. If there was publication bias, stratification analysis was carried in order to assess the reason for bias and strength the association between VDRgene polymorphisms and DR. All statistical tests were performed with Stata 11.0 software (StataCorp LP, College Station, TX, USA). All the P-values were two-sided and less than 0.05 were considered to be statistically significant.

Study characteristics
After detailed screening the databases, six case-control articles (11 study-samples) 14,15,[19][20][21][22] with 636 cases of any DR and 1,035 controls (including DM and normal subjects) were included in this meta-analysis. There were five study-samples researched on the FokI polymorphisms including 535 cases and 936 controls, and two study-samples studied on the BsmI polymorphisms including 179 cases and 408 controls, and two study-samples studied on the ApaI polymorphisms including 179 cases and 408 controls, and two study-samples studied on the TaqI polymorphisms including 186 cases and 281 controls. The selection process of eligible studies was shown in Figure  1. The qualities of the included studies were evaluated by NOS; the quality scores of six included studies were presented in Table 1. In addition, the characteristics of six articles were shown in Table 2. All of including studies were consistent with HWE. All studies indicated that the controls were DM without DR or non-DM subjects except two articles by Taverna et al. involved non-severe DR subjects as controls 14,22 .   Cyganek [15]

Meta-analysis
Overall, there was no link between VDR-FokI polymorphism and DR [F vs. f: OR = 0.68 (95% CI: 0.42-1.11); FF vs. ff: OR = 0.55 (95% CI: 0.14-2.11); Ff vs. ff:  Table 3. The results in Table 4, Table 5 and Table 6 indicated that there was no statistical significant association between polymorphisms in the VDR gene including BsmI, TaqI, ApaI and the risk of DR.

Publication bias
The potential publication bias was examined by using a funnel plot for the allele F vs. f against the reciprocal of its standard error ( Figure 2). Egger`s test provided con-sistent results in the overall (t=−0.22, P<0.05); When the study in which control group including non-severe DR patients was excluded, there was no significant publication bias according to funnel plot ( Figure 3).

Discussion
We performed a meta-analysis of case-control studies to evaluate the association between VDR gene polymorphisms and the risk of DR. To the best of our knowledge, Vitamin D supplements will lower the risk of DR in DM patients 23 . However, metabolic activation of vitamin D is mediated through binding to a VDR. Therefore, the expression of the VDR gene might be involved in the pathogenesis of DR. There were some studies conducted on VDR gene polymorphisms that had been investigated and shown to be associated with presence of any DR 24,25 . However, the published results were not consistent. Meta-analysis is needed to strengthen the link between VDR gene polymorphisms and DR. VDR gene is located in chromosomes 12q13-12q14 26 , and the allele of FokI polymorphism is located in the 5` end of the VDR gene. Overall, there was no significant association between VDR-FokI polymorphism and risk of DR in this meta-analysis. Notably, we detected a significantly (F allele) decreased susceptibility for DR among Chinese population, cases with DR compared with normal controls, among studies without publication bias and articles published after 2010 year for all models. In contrast, there was no significant association that was found in Caucasians, cases with DR compared with DM controls, and even cases with DR compared with T2DM controls. Another meta-analysis revealed that FokI polymorphism was possibly a risk factor for T2DM, mostly in Chinese subjects 27 . However, it was inexplicable that F allele within FokI polymorphism lowered risk for DR in normal subjects rather than DM or T2DM.
BsmI, ApaI and TaqI are the other three polymorphisms for VDR gene. In our meta-analysis, there was no significant association between these three polymorphisms and DR, which may be due to different locations of VDR gene on its chromosome. These three polymorphisms are located in the 3` UTR regions of VDR gene, but another single-nucleotide polymorphism (FokI) located in the start codon of VDR gene. According to the results, there was significant heterogeneity detection in some genetic models. Although we have carried out stratification analysis, significant heterogeneity was also detected. After reviewing all included studies, different ethnicity and control subjects may contribute to the heterogeneity. Therefore, the conclusions should be interpreted with caution. Another previous meta-analysis 28 revealed that DR was significantly associated with VDR-FokI polymorphism both in overall and sub-group analysis which was not consistent with our results. Zhang et al. study has some limitations. There were seven studies included in that meta-analysis but lacking of the sub-groups analysis stratified by control population (diabetic and normal subjects) or publication date. Moreover, there is lacking of detailed description on study design in the research by Bućan K et al. 25 which was also included in the Zhang`s meta-anal-ysis. So we excluded this study and only six case-control studies were included in our current meta-analysis. This meta-analysis has detected the association between VDR gene polymorphism and DR in stratification analysis. However, there were still several limitations in this study. First, relatively limited number of case-control studies was included in meta-analysis. Thus, further studies involving large sample sizes are needed to archive more reliable conclusions on association between VDR gene polymorphisms and DR. Secondly, no analysis of combined genotype effects was done due to the small sample size. Third, the control groups within two eligible studies were including non-severe DR patients, which could lead to potential publication bias. Last but not least,, we did not include articles published in other languages only in English or Chinese. Although we have tried our best to review related articles in all databases, there were many unpublished, insufficient data studies missing.

Conclusion
This meta-analysis revealed a significant association between DR and VDR-FokI polymorphism in Chinese subjects, among included studies without publication bias, during a comparison analysis between normal subjects and DR patients, and among articles published after 2010. However, owing to the limitations of this current meta-analysis, further prospective studies with larger sample sizes are needed to confirm our findings.