Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Metabolic impact of the VDR rs1544410 in diabetic retinopathy

Abstract

Aims

To investigate the association between BsmI and DM2 in patients with and without DR and to correlate with clinical parameters in a population in northeastern Brazil.

Methods

Cross-sectional case-control study in which data were collected from 285 individuals, including 128 patients with DM2 and 157 with DR. Clinical, biochemical and anthropometric parameters were analyzed, in addition to the single nucleotide polymorphism (SNP) BsmI of the VDR gene (rs1544410), genotyped by PCR-RFLP.

Results

In the DR group we found a greater number of patients using insulin therapy (p = 0.000) and with longer duration of DM2 (p = 0.000), in addition to higher serum creatinine values (p = 0.001). Higher fasting glucose levels and higher frequency of insulinoterapy were independently observed in patients with DR and b allele carriers, when compared to BB.

Conclusion

The association of the bb/Bb genotypes (rs1544410) of the VDR gene with increased blood glucose levels and insulinoterapy may represent worse glicemic control in rs1544410 b allele carriers in DR Latin American individuals.

1. Introduction

Diabetes mellitus is global epidemic health issue that affects around 415 million adults [1] Diabetic retinopathy (DR) is a frequent complication in type 2 Diabetes Mellitus (DM2) [2, 3], and clinical and metabolic factors are associated with the development and progression of DR [46].

There is evidence that deficiency of Vitamin D (VD) is related to DM2 [7, 8]. Vitamin D may have a direct effect on the function of pancreatic beta cells, mediated by the vitamin D receptor (VDR) [9], in addition, altered VDR gene transcription may influence fasting glucose levels by two potentially additive effects of vitamin D on adipocytes and pancreas cells [10, 11]. There is evidence that Vitamin D deficiency affects the pathogenesis and progression of DR and that patients with proliferative diabetic retinopathy (PDR) have lower levels of 25(OH)D than those without diabetes [8, 12]. A meta-analysis encompassing 17,000 patients from several continents concluded that vitamin D deficiency increased the risk of DR in T2DM [13].

The relationship between vitamin D levels and DM2 and DR raises interest in investigating the effect of genetic aspects linked to this vitamin regarding its involvement in the etiology or modulation of DR, including single nucleotide polymorphisms (SNPs) in the Vitamin D receptor (VDR). The B allele of the BsmI polymorphism (rs1544410) of the VDR gene was associated with a lower risk of DR in Korean patients with DM2 [14], and the bb genotype was associated with a decrease in (25 [OH] D) in micro and macrovascular complications of DM2 in an Indian population [15]. However, there is no consensus in the literature on these relationships [13, 1618], and the results seem to be influenced by the to geographic origin of patients.

Given the facts presented in the absence of data in the Brazilian population, we propose here to investigate the association between BsmI and DM2 in patients with and without DR and to correlate with clinical parameters using a case-control approach, in a population of northeastern Brazil.

2. Methods

2.1 Study design and ethical aspects

This cross-sectional case-control study was conducted in 285 patients recruited at the University Hospital Lauro Wanderley and the Basic Health Units of João Pessoa/Paraíba/Brazil. Data collection took place from 2014 to 2021 and patients invited to participate in the study signed an informed consent form. The project was approved by the Ethics Committee for Research with Human Beings of the Federal University of Paraíba (UFPB; Opinions No.: 796.459 (26/08/2014) and 3053068 (03/12/2018) ards of the institution and Resolution 466/2012 of the National Health Council.

Inclusion criteria were diagnosis of DM2 for at least 5 years, being in outpatient care. Exclusion criteria: diagnosis of DM1, insufficient DNA sample or with an inconclusive result in the genotypic analysis.

2.1.1 Clinical characterization.

The diagnosis of DR was based on ophthalmoscopy after pupil dilation with 0.5% tropicamide. Images of the retina (macula and central disc) were captured at a 45° angle by a background camera. The images were analyzed according to the standards and recommendations of ACCORD (Action to Control Cardiovascular Risk in Diabetes) and the Early Treatment Diabetic Retinopathy Study (ETDRS). DM2 patients without DR (n = 128) constituted the control group and the group with some degree of retinopathy constituted the DR (n = 157). In the RD group, 49 patients had been diagnosed with DM2 for up to ten years; the other 108 patients had been diagnosed with DM2 for more than ten years. Regarding the classification according to DR stage, 103 participants had non-proliferative diabetic retinopathy (NPDR) and 54 of them had proliferative diabetic retinopathy (PDR). Patients with a medical diagnosis of hypertension or who reported presenting this condition and were taking long-term medication for blood pressure control were considered hypertensive. The classification was performed according to the American Heart Association Guideline, American College of Cardiology, and the American Society of Hypertension (2015) [18].

2.2. Blood sampling

The blood samples were collected by venipuncture after night fasting. For biochemical analysis in general, blood in the presence of clot activator, for the determination of HbA1C in the presence of anticoagulant K3EDTA and for the determination of glucose in the presence of anticoagulant sodium fluoride, the samples were centrifuged at a speed of 3000 rpm for 10 minutes at room temperature for the separation of serum or plasma and subjected to analysis within 2 hours after collection, except for the sample of HbA1C that was analyzed in hemolyzed whole blood.

For DNA extraction, blood collection was performed by a venous puncture in sterile tubes containing 7.2 mg of K3 EDTA. Blood samples were stored for up to 20 days at -20°C until DNA extraction was performed.

2.3. Biochemical analysis and anthropometric measurements

Enzymatic methods were used for total cholesterol, high-density lipoprotein (HDL) and triglyceride analyses. All the analyses were performed using an automated analyzer (Lab-Max 240, Labtest, Lagoa Santa, MG, Brazil) and standardized kits according to the manufacturer’s instructions (Labtest, Lagoa Santa, MG, Brazil).

Low-density lipoprotein (LDL) concentration was calculated using the Friedewald formula: [LDL] = [total cholesterol] − [HDL] − [triglycerides ÷ 5] [19]. Glycated hemoglobin (HbA1c), C-reactive protein (CRP) levels, were determined by an immunoturbidimetry technique (LabMax 240, Labtest, Lagoa Santa, MG, Brazil) using commercial kits (Labtest, Lagoa Santa, MG, Brazil).

For the classification of dyslipidemia, the criteria were established according to the Guidelines of the Brazilian Diabetes Society (2019): Total Cholesterol (<190mg / dL), LDL-C (<130mg / dL), HDL-C (≤40mg / dL) and Triglycerides (≥150mg / dL) [20]. For the values of Hba1c (≤6,5%), total cholesterol (≤190mg / dL), HDL (≥40mg / dL) and LDL (<130mg / dL), the cutoff points adopted followed the recommendations of the Guidelines of the Brazilian Diabetes Society (2019) [20]. For Triglycerides (<150mg / dL), based on the Brazilian Archives of Cardiology (2017) [21]. Serum creatinine according to sex, (women ≤0.995mg / dL; men ≤1.20mg / dL), according to the Brazilian Society of Nephrology (2011) [22]. C Reactive Protein (≤3mg / dL) according to the Brazilian Guidelines on Dyslipidemias and Atherosclerosis Prevention [23], and for Malondialdehyde, the reference value adopted was (≤3.31μM) [24].

The anthropometric variables were body weight (kg), which was measured with the use of a scale, and height (cm). The body mass index (BMI) was calculated by dividing body weight by height squared (in meters) and the individuals were classified according to the presence of overweight and obesity, according to Body Mass Index (BMI) in kg / m2. The values for the adult age were: overweight when BMI = 25–29.9 kg/m2 and obesity when BMI ≥ 30kg /m2. For the elderly patients, overweight was defined as BMI = ≥27 kg / m2 [25]. Clinical and anthropometric information such as gender, age, T2DM time, BMI and blood pressure were obtained in the clinical evaluation by a nutritionist from the research team and by the team in the endocrinology service.

2.4. Isolation of leukocyte DNA

To obtain leukocyte DNA, appropriate protocols were used. The samples were diluted in an initial lysis solution containing 10mM Tris-HCl pH, 8.5mM EDTA, 0.3M sucrose and 1% Triton-X-100. Centrifugation was performed at 3.200 rpm, and the supernatant was discarded. This process was repeated 3 times to obtain a leukocyte precipitate free from hemoglobin remnants. The precipitate was resuspended in a lysis solution containing 10mM Tris-HCl pH8.0, 0.5% sodium dodecyl sulfate (SDS), 5mM EDTA and 0.2μg proteinase K (Invitrogen, Carlsbad, CA, USA) and incubated at 55°C in a water bath for 7h. Then, 500μl of an aqueous solution of 1mM EDTA and 7.5M ammonium acetate was added and mixed for 30 seconds. The mixture was centrifuged for 10min at 14.000g at 4°C, and 700μl of the supernatant was transferred to a new tube where DNA precipitation was performed with 540μl of iced isopropanol. The DNA precipitate was washed with 70% ethanol, centrifuged (12.000g for 5min), dried and resuspended in Tris-EDTA buffer pH 8.0 [26]. The samples were kept at -20° C until genetic analysis.

2.5. rs1544410 genotyping

Genotypes were determined by PCR-RFLP. Appropriate primers were used to amplify the region of the gene containing the polymorphism [27], and amplification occurred under the following conditions: denaturation 94° C for 5 minutes, 30 cycles of denaturation (1 minute at 94° C), annealing (1 minute at 58° C) and extension (3 minutes at 72° C) with an extra 10 minute extension step minutes. The 825 bp product was digested with BsmI which recognizes and cleaves the polymorphic allele (b) generating two fragments (650 bp and another 175 bp) while the wild allele (B) remains at 825 bp. The genotypes were analyzed by 15% polyacrylamide gel electrophoresis and 0.5% silver nitrate staining.

2.6. Statistical analysis

SPSS 26.0 software (SPSS Inc., Chicago, IL, USA) was used for statistical analysis. Normality in continuous variables was assessed using the Kolmogorov-Smirnov test, and distributions with p> 0.05 were accepted as variables with normal distribution and expressed as mean ± standard deviation values and evaluated by the t test of independent samples. The ’non normally distributed’ were expressed as median values and 95% confidence intervals and compared using the Mann-Whitney U test between groups. Nominal categorical variables were expressed as total number and percentage, analyzed by chi-square. Hardy-Weinberg equilibria were calculated to assess expected and observed genotypic and allelic frequencies were tested by chi-square. Data from the RD group were analyzed according to rs1544410 genotypes by independent T test and its non-parametric counterpart, and then a logistic regression model was used to establish which variables (fasting glucose, insulin therapy, sex, age and time of DM2) could be influenced by rs1544410 genotypes. A p value < 0.05 was considered significant for all analyses.

3. Results

The control (DM) and RD groups were similar regarding the frequency of hypertension, dyslipidemia, family history of DM, smoking, sedentary lifestyle, age, HbA1c, fasting glucose, LDL-c, triglycerides and BMI, however, in the group with DR, male gender (p = 0.001), insulin therapy (p = 0.000), and serum creatinine (p = 0.000) were higher in addition to a longer duration of DM2 (p = 0.000). In the DM group, patients had higher values of total cholesterol (p = 0.017) and HDL-c fraction (p = 0.008) (Table 1).

thumbnail
Table 1. Clinical, biochemical and metabolic parameters of the studied sample.

https://doi.org/10.1371/journal.pone.0263346.t001

Table 2 shows the distributions of genotypes (BB + Bb / bb; Bb + bb / BB; BB + bb / Bb) and alleles (B / b) as a function of the DM control group compared to different DR categorizations; DR with ≤10 years of DM2 diagnosis; DR with >10 years of DM2 diagnosis; groups with RDNP and RDP, respectively. In all groups there was no difference in the genotypic or allelic distribution of the groups. The minor allelle (B) frequency (MAF) was 0.46 in control DM group and 0.43 in the DM RD group.

thumbnail
Table 2. Genotypic and allelic distribution of the polymorphism rs1544410 among type 2 diabetic patients without DR and RD, <10 years DM diagnostic, ≥10 years DM diagnostic, RDNP and RDP.

https://doi.org/10.1371/journal.pone.0263346.t002

Clinical parameters of the RD group (n = 157) were analyzed according to BsmI genotypes (BB x Bb+ bb). In the group with the Bb+bb genotype, insulin therapy was significantly more frequent (p = 0.014) and fasting glucose had a higher mean value in the same group (p = 0.047) (Table 3). A logistic regression model, adopting the BsmI genotypes as the dependent variable and with correction/adjustment variables, revealed that fasting glucose (p = 0.018) and insulin therapy (p = 0.036) were statistically significant predictors of increased odds of risk in the group of patients with the Bb+bb genotype (Table 4). The glycemia difference between the genotypic groups presents a Cohen’s effect size d = 0.5, both comparing only the BB x bb group and Bb+bb x BB. According to the literature, this value means an effect of medium magnitude [28].

thumbnail
Table 3. Clinical, biochemical and metabolic parameters of the RD group according to rs1544410 genotypes.

https://doi.org/10.1371/journal.pone.0263346.t003

thumbnail
Table 4. Logistic regression model analyzing the influence of rs1544410 genotypes and metabolic variables in patients with DR.

https://doi.org/10.1371/journal.pone.0263346.t004

4. Discussion

Vitamin D in its active form acts through its specific receptor (VDR), expressed in human tissues including the retina [29]. In addition to the association of lower levels of Vitamin D and increased risk of DM2 and RD [8, 13], several studies have also related VDR gene polymorphisms as risk factors for DM2 [30], pathogenesis [14, 15, 31], and DR progression [32].

The relationship of polymorphisms with the occurrence of chronic complications in diabetics has also been explored [3335] and has a pronounced importance in the early identification of individuals at higher risk of developing specific complications. The SNP polymorphism rs1544410 is located in intron 8 in the 3’ regulatory region, which is involved in gene expression, especially through the regulation of messenger RNA stability [36, 37], which can result in reduced expression of the VDR gene [38]. Another possible mechanism is that the alternation of BsmI in the intronic sequence may influence protein expression [39].

The MAF found in this study was 0.47 in the diabetic group and 0.43 in DR patients while in another Brazilian sample from Minas Gerais State (Southeast region) of the country, the MAF found was 0.40 in the group of type 2 diabetics. In the samples of the 1000 genome project, obtained from healthy population, the MAF (A) for the SNP is 0.29. According to the SNP database (ncbi.nlm.gov) in two different Latin American populations the MAF found was 0.36 (n = 798) and 0.23 (n = 3896), in Europeans it was 0.39 and in African Americans 0.26. Interestingly, in Asian populations the MAF is significantly lower, showing that it is a SNP that is highly influenced by ethnic aspects of the population under analysis. These informations evidenciate similarity between the MAF found in this study and another sample of Brazilian diabetics and a Latin American sample.

The genotypic distribution found deviated from the Hardy Weinberg Equilibrium in both groups, DM2 and DM2 with DR. This finding has not been uncommon in studies involving the polymorphism in question and diabetes. The HWE imbalance was present in a study conducted in the USA [40] and another in Chile [41]. In a recent meta-analysis that evaluated the effect of rs1544410 on the risk of DM2, of the 37 studies analyzed, 10 showed a deviation from the HWE in samples from different locations around the world [42].

The “b” allele was associated with lower HDL levels and increased risk for T2DM and obesity in Arabs [43], and the occurrence of DR in Croatian population with DM1 [32] the bb genotype conferred greater insulin resistance in Caucasians [40].

Despite the evidence, there are few studies exploring the association of the SNP rs1544410 and RD; one of them in DM1 [32] and four in DM2 [14, 15, 44, 45]. In these studies, one found an association of the B allele with lower risk of DR in DM2 [14] and two studies found association between RD and the bb genotype in DM1 [32] and DM2 [15], whereas in other studies no genotypic or allelic association was found [44, 45]. The geographic regions involved in these results cover Europe [32, 44] and Asia [14, 15, 45], however there are no studies analyzing the relationship in the Latin American population, which characterizes this study as pioneer. The results of the present study do not show an association of the rs1544410 polymorphism with DR in any of the studied scenarios, either by classifying patients by duration of diabetes, by DR stage or by simple diagnosis of the condition.

The time of diabetes diagnosis was also evaluated by Búcan (2009), Hong (2015), Zhong (2015) and Cyganek (2006), and the results found reinforce the data obtained in the study presented here. A single study [14] used the RD stage classification variable controlling the time since DM2 diagnosis (mean of 13 years) between the groups. In this study, the BB, Bb and the B allele genotypes were related to a lower risk of DR.

Meta-analyses [17, 29] point to the absence of association between rs1544410 and DR in both DM1 and DM2. However, the meta-analysis conducted by Song (2019) indicates an association restricted to the Indian population, but that is the result of a single study [15]. Therefore, we consider it important to explore this possibility in the Brazilian population. As the association did not result in positive data, we decided to investigate the effect of different rs1544410 genotypes in patients with DR on clinical and metabolic parameters.

Among the existing studies with DM2 and DR, only one controlled the duration of DM2 between the groups, with no statistically significant differences [14], while Zhong (2015), Cyganek (2006) and Ezhilarasi (2018) present difference or do not report the comparison. Aware of the great influence of this variable on the studied aspect, in this work, we performed different comparison scenarios, including separating patients with DR up to 10 years of diagnosis and more than 10 years of DM2 diagnosis, comparing genotypic and allelic distributions with the control group. In none of these scenarios were the genotypes or alleles statistically different from the control.

As the association analysis found conflitant results in comparison to some previous studies [14, 15], we decided to investigate the effect of different rs1544410 genotypes on metabolic parameters, in order to assess whether despite not impacting the risk of DR, having knowledge of the genotype would have an impact on the knowledge of the individual’s metabolic profile. We compared the genotypic groups within the profile of DR carriers and identified that carriers of the b allele (Bb and bb) had higher blood glucose levels and a greater chance of insulin therapy than carriers of the BB genotype.

Persistent hyperglycemia triggers most of the processes that promote the development of DR and high blood glucose concentrations and high levels of HbA1C characterize poor glycemic control and associate them with an independent risk factor for DR [5, 6, 4649].

The association between VDR variants (BsmI) and insulin resistance-related diseases, demonstrate controversial conclusions. While the BB+Bb genotypes and the B allele were related to insulin resistance-related diseases in darkly pigmented Caucasians [50], in a prospective cohort, the B allele was associated with increased insulin secretion in women with a history of gestational DM [51]. The BsmI polymorphism of the VDR gene was also associated with glycemic outcomes in other populations. BsmI was tested in a homogeneous population of healthy young men, at an age before the typical onset of DM2, in which physical activity was also measured. Individuals with the BB genotype had significantly higher fasting glucose levels than individuals with the Bb or bb genotype. This finding was associated with low physical activity (≤3h per week), but the effect was absent in men with a high degree of physical activity [52].

In addition to associations related to the vitamin D receptor gene and BsmI genotypes with glycemic markers, low serum concentrations of 25OHD3 have been associated with reduced insulin sensitivity, impaired glucose metabolism, and metabolic syndrome [53, 54]. The mechanism by which vitamin D affects insulin sensitivity is still unknown. Vitamin D can stimulate the expression of insulin receptors in peripheral tissues and, thus, increase glucose transport [55]. Insulin-mediated glucose uptake is calcium dependent and therefore vitamin D status may indirectly influence glucose uptake [56].

The b allele was considered a risk factor for low serum concentrations of 25-hydroxyvitamin D, in addition the same study reports that BsmI influenced glucose and HOMA-IR concentrations even after adjustments in Brazilian (Amazon State) children [57].

These evidences demonstrate the importance of going beyond the simple genotype-clinical condition relationship. It is possible that the lack of association found in some populations, including this study, is a reflection of the applied experimental design, and the influencing effect of other aspects not measured or analyzed, including other SNPs. The association between higher glucose levels and the b allele is relevant in the context studied, as this same genetic characteristic seems to influence vitamin D levels. The results found here encourage further analysis, taking into account clinical, metabolic and vitamin D levels of patients when evaluating the association of VDR genotypes and DR.

This is the first study to analyze a Latin American population for the influence of rs1544410 in the context of DR. The study has limitations related to its cross-sectional nature. The results found here encourage the performance of studies with other populations and with different experimental models to unravel the effects of other important components for DR.

5. Conclusion

We did not identify the BsmI SNP (rs1544410) of the VDR gene as a risk factor for DR, but the b allele may be indicative of increased glycemic levels and insulinoterapy in the population with DR. The relationship between VDR BsmI and susceptibility to DR in a Latin American population should be further studied.

Acknowledgments

We thank at the Lauro Wanderley University Hospital, Federal University of Paraíba and all participants in the study for participation and the financial support institutions, National Council for Scientific and Technological Development (CNPq, Brasilia, Brazil), the Coordination for the Improvement of Higher Education Personnel (CAPES, Brasília, Brazil) and the State of Paraíba Research Support Foundation (FAPESQ, Paraíba, Brazil). The authors declare that they have no competing interests.

References

  1. 1. Ogurtsova K, da Rocha Fernandes JD, Huang Y, et al. IDF Diabetes Atlas: Global estimates for the prevalence of diabetes for 2015 and 2040. Diabetes Res Clin Pract. 2017 Jun; 128:40–50. pmid:28437734
  2. 2. Wong TY, Cheung CM, Larsen M, et al. Diabetic retinopathy. Nat Rev Dis Primers. 2016 Mar 17;2:16012. pmid:27159554
  3. 3. Duh EJ, Sun JK, Stitt AW. Diabetic retinopathy: current understanding, mechanisms, and treatment strategies. JCI Insight. 2017 Jul 20;2(14):e93751. pmid:28724805
  4. 4. Chen YH, Chen HS, Tarng DC. More impact of microalbuminuria on retinopathy than moderately reduced GFR among type 2 diabetic patients. Diabetes Care. 2012 Apr;35(4):803–8. pmid:22338100
  5. 5. Kawasaki R, Tanaka S, Tanaka S, et al. Japan Diabetes Complications Study Group. Incidence and progression of diabetic retinopathy in Japanese adults with type 2 diabetes: 8 year follow-up study of the Japan Diabetes Complications Study (JDCS). Diabetologia. 2011 Sep;54(9):2288–94. pmid:21630126
  6. 6. Stratton IM, Kohner EM, Aldington SJ, et al. UKPDS 50: risk factors for incidence and progression of retinopathy in Type II diabetes over 6 years from diagnosis. Diabetologia. 2001 Feb;44(2):156–63. pmid:11270671
  7. 7. Maddaloni E, Cavallari I, Napoli N, et al. Vitamin D and Diabetes Mellitus. Front Horm Res. 2018;50:161–176. pmid:29597238
  8. 8. Zhang J, Upala S, Sanguankeo A. Relationship between vitamin D deficiency and diabetic retinopathy: a meta-analysis. Can J Ophthalmol. 2017 Nov;52 Suppl 1:S39–S44. pmid:29074012
  9. 9. Johnson JA, Grande JP, Roche PC, et al. Immunohistochemical localization of the 1,25(OH)2D3 receptor and calbindin D28k in human and rat pancreas. Am J Physiol. 1994 Sep;267(3 Pt 1):E356–60. pmid:7943215
  10. 10. Zemel MB, Hang S, Greer B, et al. Regulation of adiposity by dietary calcium. FASEB J 2000 Jun; 14: 1132–1138. pmid:10834935
  11. 11. Norman AW, Frankel JB, Heldt AM, et al. Vitamin D deficiency inhibits pancreatic secretion of insulin. Science. 1980 Aug 15;209(4458):823–5. pmid:6250216
  12. 12. Payne JF, Ray R, Watson DG, et al. Vitamin D insufficiency in diabetic retinopathy. Endocr Pract. 2012 Mar-Apr;18(2):185–93. pmid:21940279
  13. 13. Luo BA, Gao F, Qin LL. The Association between Vitamin D Deficiency and Diabetic Retinopathy in Type 2 Diabetes: A Meta-Analysis of Observational Studies. Nutrients. 2017 Mar 20;9(3):307. pmid:28335514
  14. 14. Hong YJ, Kang ES, Ji MJ, et al. Association between Bsm1 Polymorphism in Vitamin D Receptor Gene and Diabetic Retinopathy of Type 2 Diabetes in Korean Population. Endocrinol Metab (Seoul). 2015 Dec;30(4):469–74. pmid:26790383
  15. 15. Ezhilarasi K, Dhamodharan U, Vijay V. BSMI single nucleotide polymorphism in vitamin D receptor gene is associated with decreased circulatory levels of serum 25-hydroxyvitamin D among micro and macrovascular complications of type 2 diabetes mellitus. Int J Biol Macromol. 2018 Sep;116:346–353. pmid:29738868
  16. 16. Song N, Yang S, Wang YY, et al. The Impact of Vitamin D Receptor Gene Polymorphisms on the Susceptibility of Diabetic Vascular Complications: A Meta-Analysis. Genet Test Mol Biomarkers. 2019; 23(8):533–556. pmid:31373850
  17. 17. Jiao J, Li Y, Xu S, et al. Association of FokI, TaqI, BsmI and ApaI polymorphisms with diabetic retinopathy: a pooled analysis of case-control studies. Afr Health Sci. 2018 Dec;18(4):891–899. pmid:30766552
  18. 18. Rosendorff C, Lackland DT, Allison M, et al. American Heart Association, American College of Cardiology, and American Society of Hypertension. Treatment of hypertension in patients with coronary artery disease: a scientific statement from the American Heart Association, American College of Cardiology, and American Society of Hypertension. Hypertension. 2015 Jun;65(6):1372–407. pmid:25828847
  19. 19. 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 Jun; 18(6): 499–502. pmid:4337382
  20. 20. Brazilian Diabetes Society. GUIDELINES OF THE BRAZILIAN DIABETES SOCIETY CLANNAD 2020. [http://www.saude.ba.gov.br/wp-content/uploads/2020/02/Diretrizes-Sociedade-Brasileira-de-Diabetes-2019-2020.pdf]
  21. 21. Précoma DB, Oliveira GMM, Simão AF, et al. Updated Cardiovascular Prevention Guideline of the Brazilian Society of Cardiology—2019. Arq Bras Cardiol. 2019 Nov 4;113(4):787–891. Erratum in: Arq Bras Cardiol. 2021 Apr;116(4):855. pmid:31691761
  22. 22. Brazilian Society of Nephrology e-book Biomarkers in Nephrology 2011. [https://arquivos.sbn.org.br/pdf/biomarcadores.pdf]. Accessed on 05/04/2021
  23. 23. Faludi AA, Izar MCO, Saraiva JFK, et al. Update of the Brazilian Dyslipidemia and Atherosclerosis Prevention Directive- 2017. Arq Bras Cardiol 2017; 109(2) (Suppl. 1): 1–76.
  24. 24. Antunes MV, Lazzaretti C, Gamaro GD, et al. Estudo pré-analítico e de validação para determinação de malondialdeído em plasma humano por cromatografia líquida de alta eficiência, após derivatização com 2,4-dinitrofenilhidrazina. Rev Bras Cien Farm 2008 Ago; 44: 279–87
  25. 25. WHO Consultation on Obesity (1999: Geneva, Switzerland) & World Health Organization. (2000). Obesity: preventing and managing the global epidemic: report of a WHO consultation. World Health Organization. [https://apps.who.int/iris/handle/10665/42330]
  26. 26. Miller SA, Dykes DD, Polesky HF. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988 Feb 11;16(3):1215. pmid:3344216
  27. 27. Pani MA, Knapp M, Donner H, et al. Vitamin D receptor allele combinations influence genetic susceptibility to type 1 diabetes in Germans. Diabetes. 2000 Mar;49(3):504–7. pmid:10868975
  28. 28. Cohen J. Statistical power analysis for the behavioral sciences 2d ed. New York: Academic Press. 1988.
  29. 29. Zhang Y, Xia W, Lu P, et al. The Association between VDR Gene Polymorphisms and Diabetic Retinopathy Susceptibility: A Systematic Review and Meta-Analysis. Biomed Res Int. 2016 Nov; 2016:5305282. pmid:27891515
  30. 30. Safar HA, Chehadeh SEH, Abdel-Wareth L, et al. Vitamin D receptor gene polymorphisms among Emirati patients with type 2 diabetes mellitus. J Steroid Biochem Mol Biol. 2018 Jan;175:119–124. Epub 2017 Mar 18. pmid:28323045
  31. 31. Motohashi Y, Yamada S, Yanagawa T, et al. Vitamin D receptor gene polymorphism affects onset pattern of type 1 diabetes. J Clin Endocrinol Metab. 2003 Jul;88(7):3137–40. pmid:12843155
  32. 32. Bućan K, Ivanisević M, Zemunik T, et al. Retinopathy and nephropathy in type 1 diabetic patients—association with polymorphysms of vitamin D-receptor, TNF, Neuro-D and IL-1 receptor 1 genes. Collegium antropologicum [Internet]. 2009 Dec [cited 2021 Sep 2];33 Suppl 2:99–105. Available from: https://hrcak.srce.hr/index.php?show=clanak&id_clanak_jezik=79096 pmid:20120526
  33. 33. Vedralová M, Kotrbova-Kozak A, Zelezníková V, et al. Polymorphisms in the vitamin D receptor gene and parathyroid hormone gene in the development and progression of diabetes mellitus and its chronic complications, diabetic nephropathy and non-diabetic renal disease. Kidney & Blood Pressure Research. 2012;36(1):1–9. pmid:22777106
  34. 34. Ferrarezi DAF, Bellili- Muñoz N, Dubois-Laforgue D, et al. Allelic variations of the vitamin D receptor (VDR) gene are associated with increased risk of coronary artery disease in type 2 diabetics: the DIABHYCAR prospective study. Diabetes Metab 2013 May;39:263–270. pmid:23352876
  35. 35. Zhang H, Wang J, Yi B, et al. BsmI polymorphisms in vitamin D receptor gene are associated with diabetic nephropathy in type 2 diabetes in the Han Chinese population. Gene. 2012 Mar 10;495(2):183–8. pmid:22245613
  36. 36. Durrin LK, Haile RW, Ingles SA, et al. Vitamin D receptor 3’-untranslated region polymorphisms: lack of effect on mRNA stability. Biochim Biophys Acta. 1999 Mar 30;1453(3):311–20. pmid:10101249
  37. 37. Uitterlinden AG, Fang Y, Van Meurs JB, et al. Genetics and biology of vitamin D receptor polymorphisms. Gene. 2004 Sep 1;338(2):143–56. pmid:15315818
  38. 38. Carvalho C, Marinho A, Leal B, et al. Association between vitamin D receptor (VDR) gene polymorphisms and systemic lupus erythematosus in Portuguese patients. Lupus. 2015 Jul;24(8):846–53. pmid:25661837
  39. 39. Györffy B, Vásárhelyi B, Krikovszky D, et al. Gender-specific association of vitamin D receptor polymorphism combinations with type 1 diabetes mellitus. Eur J Endocrinol. 2002 Dec;147(6):803–8. pmid:12457456
  40. 40. Oh J. L; Barrett-Connor E. Association BetweenVitamin D Receptor Polymorphism and Type 2 Diabetes or Metabolic Syndrome in Community-Dwelling Older Adults: The Rancho Bernardo Study. Metabolism. 2002 May; 51 (3):356–359. pmid:11887173
  41. 41. Angel B, Lera L, Márquez C, Albala C. The association of VDR polymorphisms and type 2 diabetes in older people living in community in Santiago de Chile. Nutr Diabetes. 2018 May 25;8(1):31. pmid:29795525
  42. 42. Liu Y, Guo X, Huang SY, Gong L, Cui JH, Shen HW, et al. Evaluation of association studies and a systematic review and meta-analysis of VDR polymorphisms in type 2 diabetes mellitus risk. Medicine (Baltimore). 2021 Jul 16;100(28):e25934. pmid:34260520
  43. 43. Al-Daghri NM, Al-Attas OS, Alkharfy KM, et al. Association of VDR-gene variants with factors related to the metabolic syndrome, type 2 diabetes and vitamin D deficiency. Gene. 2014 Jun;542(2):129–133. pmid:24680778
  44. 44. Cyganek K, Mirkiewicz-Sieradzka B, Malecki MT, et al. Clinical risk factors and the role of VDR gene polymorphisms in diabetic retinopathy in Polish type 2 diabetes patients. Acta Diabetol. 2006 Dec;43(4):114–9. pmid:17211561
  45. 45. Zhong X, Du Y, Lei Y, et al. Effects of vitamin D receptor gene polymorphism and clinical characteristics on risk of diabetic retinopathy in Han Chinese type 2 diabetes patients. Gene. 2015 Jul 25;566(2):212–6. Epub 2015 Apr 18. pmid:25899017
  46. 46. Klein R, Knudtson MD, Lee KE, et al. The Wisconsin Epidemiologic Study of Diabetic Retinopathy: XXII the twenty-five-year progression of retinopathy in persons with type 1 diabetes. Ophthalmology. 2008 Nov;115(11):1859–68. pmid:19068374
  47. 47. Kim JH, Kwon HS, Park YM, et al. Prevalence and associated factors of diabetic retinopathy in rural Korea: the Chungju metabolic disease cohort study. J Korean Med Sci. 2011 Aug;26(8):1068–73. pmid:21860558
  48. 48. Salinero-Fort MÁ, San Andrés-Rebollo FJ, de Burgos-Lunar C, et al. Four-year incidence of diabetic retinopathy in a Spanish cohort: the MADIABETES study. PLoS One. 2013 Oct 17;8(10):e76417. pmid:24146865
  49. 49. Xu J, Xu L, Wang YX, et al. Ten-year cumulative incidence of diabetic retinopathy. The Beijing Eye Study 2001/2011. PLoS One. 2014 Oct 27;9(10):e111320. pmid:25347072
  50. 50. Han FF, Lv YL, Gong LL, et al. VDR Gene variation and insulin resistance related diseases. Lipids Health Dis. 2017 Aug 19;16(1):157. pmid:28822353
  51. 51. Shaat N, Katsarou A, Shahida B, et al. Association between the rs1544410 polymorphism in the vitamin D receptor (VDR) gene and insulin secretion after gestational diabetes mellitus. PLoS One. 2020 May 14;15(5):e0232297. pmid:32407388
  52. 52. Ortlepp JR, Metrikat J, Albrecht M, et al. The vitamin D receptor gene variant and physical activity predicts fasting glucose levels in healthy young men. Diabet Med. 2003 Jun;20(6):451–4. pmid:12786678
  53. 53. Alvarez JA, Ashraf A. Role of vitamin d in insulin secretion and insulin sensitivity for glucose homeostasis. Int J Endocrinol. 2010;2010:351385. pmid:20011094
  54. 54. Pittas AG, Lau J, Hu FB, et al. The role of vitamin D and calcium in type 2 diabetes. A systematic review and meta-analysis. J Clin Endocrinol Metab. 2007 Jun;92(6):2017–29. Epub 2007 Mar 27. pmid:17389701
  55. 55. Maestro B, Campión J, Dávila N, et al. Stimulation by 1,25-dihydroxyvitamin D3 of insulin receptor expression and insulin responsiveness for glucose transport in U-937 human promonocytic cells. Endocr J. 2000 Aug;47(4):383–91. pmid:11075718
  56. 56. Ojuka EO. Role of calcium and AMP kinase in the regulation of mitochondrial biogenesis and GLUT4 levels in muscle. Proc Nutr Soc. 2004 May;63(2):275–8. pmid:15294043
  57. 57. Cobayashi F, Lourenço BH, Cardoso MA. 25-Hydroxyvitamin D3 Levels, BsmI Polymorphism and Insulin Resistance in Brazilian Amazonian Children. Int J Mol Sci. 2015 Jun 3;16(6):12531–46. pmid:26047339