{"paper_id":"c7bc94bb-234c-4fdb-9a23-d9b1e23074a0","body_text":"Diabetic retinopathy (DR) is a microvascular complication of type 2 diabetes mellitus (T2DM) and a leading cause of preventable blindness in adults [ 1 – 3 ]. DR is a spectrum that starts with non-proliferative DR (NPDR) that eventually progresses to proliferative DR (PDR) [ 1 – 3 ]. The worldwide prevalence of DR is 24.6% in patients with diabetes [ 4 ], the highest incidence rates are observed in African Americans, Hispanics, and South Asians [ 3 ]. In China, among patients with diabetes, the prevalence rates of DR and vision-threatening DR are 27.9 and 12.6%, respectively [ 5 ]. The management of DR includes glycemic control, lipidemia control, blood pressure control, laser therapy, intravitreal anti-vascular endothelial growth factor (VEGF) therapy, vitrectomy, and intravitreal steroids [ 2 , 6 ].\nVEGF is a key molecule in the development of DR [ 7 , 8 ]. Intravitreal VEGF levels are high in diabetic patients and correlate with neovascularization [ 9 ], especially in patients with PDR [ 3 , 10 ]. In addition to VEGF levels, the expression and distribution of the VEGF receptor (VEGFR)-2 are also crucial in developing DR and progressing from NPDR to PDR [ 11 – 13 ]. Indeed, VEGFR-2 controls endothelial cell proliferation, migration, survival, and osmosis via its activation by VEGF [ 8 ].\nSingle nucleotide polymorphisms (SNPs) in the  VEGFR-2  gene can affect the development and progression of multiple diseases by affecting the binding efficiency between VEGFR-2 with VEGF or by reducing the transcriptional activity of VEGFR-2 [ 14 – 16 ]. The involvement of some  VEGF  SNPs in the development of DR and response to treatment are well-known [ 17 – 20 ], but the role of  VEGFR-2  SNPs in the development of DR is poorly known. The rs2071559 SNP is located in the promoter of  VEGFR-2  and decreases the transcriptional activity of VEGFR-2 [ 14 , 21 , 22 ]. A previous study revealed that four SNPs in the  VEGF  gene and one (rs2071559 or -604T/C) in the  VEGFR-2  gene were associated with the development of DR [ 23 ], but it did not examine the association of the rs2071559 SNP with PDR. There are several studies on the polymorphisms of  VEGF  in DR [ 17 , 23 – 25 ], but only one study of the rs2071559 SNP in  VEGFR-2  was found [ 23 ]. Since DR is a complex trait affected by multiple genes and the environment, exploring the role of  VEGFR-2  SNPs in various populations is relevant.\nTherefore, the present study aimed to examine the association between the -604T/C locus polymorphism of the  VEGFR-2  gene and the risk of DR among T2DM patients and then of PDR among DR patients in a Han population in North China.\n\nThis retrospective case–control study included patients with T2DM who visited the Department of Ophthalmology and the Department of Endocrinology of the Fourth Hospital of Hebei Medical University between July 2014 and July 2016. T2DM was diagnosed according to the International World Health Organization (WHO) T2DM diagnostic criteria (1999) [ 26 ]. The healthy controls were unrelated individuals who received routine physical examinations between July 2014 and July 2016. All subjects were of Han ethnicity and from families living in Shijiazhuang, Hebei Province, China, for at least three generations.\nThe subjects were divided into four groups (one control group and three patient groups) according to the 2002 International Clinical Classification Criteria for DR [ 27 ]. The DM group included the T2DM patients without fundus examination abnormalities. The PDR group included the patients with one or more of the following: neovascularization, vitreous hemorrhage, and preretinal hemorrhage. The NPDR group included the patients with microangioma, intraretinal hemorrhage, and microvascular abnormalities, but without the PDR-specific findings.\nInformed consent:  Informed consent has been obtained from all individuals included in this study.\nEthical approval:  The research related to human use has been complied with all the relevant national regulations and institutional policies and in accordance with the tenets of the Helsinki Declaration, and has been approved by the ethics committee of The Fourth Hospital of Hebei Medical University.\nThe data were obtained from the medical records and included sex, age, fasting blood glucose levels, body height, body weight, systolic blood pressure, diastolic blood pressure, and T2DM duration. The body mass index (BMI) was calculated as BMI = body weight (kg)/height (m 2 ).\nDNA was routinely extracted from venous blood collected in sodium citrate tubes and stored at 4°C. Peripheral blood leukocyte DNA was extracted using the proteinase K digestion-saturated sodium chloride salting-out method within 1 week of blood collection [ 28 ].\nThe polymerase chain reaction-ligation detection reaction (PCR-LDR) method was used to determine the  VEGFR-2  -604T/C polymorphism. The forward primer for the amplification of rs2071559 was 5′-AAA TAT TTT GGG AAA TAG CGG GAA TG-3′, and the reverse primer was 5′-TGG CGA ACT GGG CAA GTG CGT TTT C-3′. The forward and reverse primers for -604T/C were 5′-GAG CAC GAT GGA CAA AAG CCT-3′ and 5′-GGG AGA AGC GGA TAC TCA GCC-3′, respectively. The PCR product length was 255 bp, and the GC content was 47.5%. The LDR probe design included two probes. The sequences of the probes for rs2071559 were TC: 5′-TTT TTT TTA TAT TTT GGG AAA TAG CGG GAA TGC-3′, TT: 5′-TTT TTT TTT TTA TAT TTT GGG AAA TAG CGG GAA TGT-3′, and TR: 5′-P-TGG CGA ACT GGG CAA GTG CGT TTT CTT TTT T-FAM-3′.\nThe size of the ligated products was 64/C and 67/T. The LDR reaction consisted of cycling through the following conditions for 25 cycles: 94°C for 30 s and 56°C for 1 min.\nSterile distilled water was used as a negative control. In addition, 10% of the DNA samples were randomly selected for duplicate experiments. All repeated results were consistent with the original ones.\nStatistical analysis was performed using SPSS 13.0 (SPSS, Chicago, IL, USA). The Hardy–Weinberg equilibrium was tested in all groups using the chi-square test. Normally distributed continuous variables were expressed as means ± standard deviations (and analyzed using one-way analysis of variance and the Bonferroni  post hoc  test. Non-normally distributed continuous variables were presented as medians (ranges) and compared using the Kruskal–Wallis test. Non-conditional logistic regression analysis was performed to calculate the odds ratios (ORs) and 95% confidence intervals (95% CIs).  P -values < 0.05 were considered statistically significant.\n\nFinally, 438 subjects were included: 114 controls, 123 patients in the DM group, 105 in the NPDR group, and 96 in the PDR group. There were no significant differences in age and sex among the four groups (all  P  > 0.05) ( Table A1 ).\nThe three genotypes (CC, CT, and TT) were present in all four groups and met the Hardy–Weinberg equilibrium (all  P  > 0.05) ( Table A2 ). The frequencies of CC, CT, and TT were 17.6, 45.6, and 36.8% in the control group, 12.2, 50.4, and 37.4% in the DM group, 7.6, 53.3, and 39.1% in the NPDR group, and 7.3, 46.9, and 45.8% in the PDR group. There were no significant differences among the four groups ( P  = 0.191) ( Table 1 ).\nDistribution of genotype and allele frequencies of  VEGFR -2 -604T/C\nDM, diabetes mellitus without DR group; NPDR, non-proliferative DR group; PDR, PDR group.\nThe C and T allele frequencies were 40.4 and 59.6% in the control group, 37.4 and 62.6% in the DM group, 34.3 and 65.7% in the NPDR group, and 30.7 and 69.3% in the PDR group ( P  = 0.198) ( Table 1 ).\nThe rs2071559 genotypes were not associated with DR in the univariable analysis ( P  > 0.05) ( Table 2 ). The duration of T2DM (OR = 1.064, 95% CI: 1.024–1.105,  P  = 0.002), blood glucose (OR = 1.801, 95% CI: 1.057–3.067,  P  = 0.030), systolic blood pressure (OR = 1.023, 95% CI: 1.009–1.038,  P  = 0.001), and diastolic blood pressure (OR = 1.047, 95% CI: 1.021–1.075,  P  < 0.001) were associated with DR. The rs2071559 SNP of the  VEGFR-2  gene was not associated with DR after adjustment for age, sex, duration of DM, blood glucose, systolic blood pressure, diastolic blood pressure, and BMI ( Table 3 ).\nUnivariable logistic regression analysis for factors associated with DR among type 2 diabetic patients\nDM, diabetes mellitus; SBP, systolic blood pressure; DBP, diastolic blood pressure; BMI, body weight index; OR, odds ratio; CI, confidence interval.\nMultivariable logistic regression analysis of the factors associated with DR among type 2 diabetic patients using different genetic models\nAdjusted for age, sex, duration of diabetes mellitus, blood glucose, SBP, DBP, and BMI.\nOR, odds ratio; CI, confidence interval.\nThe rs2071559 genotypes were not associated with PDR in the univariable analysis ( P  > 0.05) ( Table 4 ). Female sex (OR = 1.776, 95% CI: 1.010–3.125,  P  = 0.046), duration of T2DM (OR = 1.082, 95% CI: 1.032–1.134,  P  = 0.001), systolic blood pressure (OR = 1.018, 95% CI: 1.002–1.034,  P  = 0.032), and BMI > 24.0 kg/m 2  (OR = 0.358, 95% CI: 0.198–0.645,  P  = 0.001) were associated with the development of PDR. The rs2071559 SNP of the  VEGFR-2  gene was not associated with PDR in any genetic model after adjustment for age, sex, duration of DM, blood glucose, systolic blood pressure, diastolic blood pressure, and BMI ( Table 5 ).\nUnivariable logistic regression analysis for factors associated with PDR among DR patients\nPDR, proliferative DR; DM, diabetes mellitus; SBP, systolic blood pressure; DBP, diastolic blood pressure; BMI, body weight index; OR, odds ratio; CI, confidence interval.\nMultivariable logistic regression analysis of the factors associated with PDR among DR patients in different genetic models\nAdjusted for age, sex, duration of diabetes mellitus, blood glucose, SBP, DBP, and BMI.\nPDR, PDR; OR, odds ratio; CI, confidence interval.\n\nWhile the role of  VEGF  SNPs in the development of DR and response to treatment has been described [ 18 , 20 ], the role of  VEGFR-2  SNPs in the development of DR is poorly known. Therefore, this study aimed to examine the association between the rs2071559 SNP of the  VEGFR-2  gene and the risk of DR in a Han population in North China. The results suggest that the  VEGFR-2  -604T/C polymorphism (rs2071559) is not associated with DR among T2DM patients or PDR among patients with DR in the Han population of Shijiazhuang (China).\nDR is a complex trait affected by numerous genes and environmental factors [ 29 , 30 ], and the pathogenesis is still unclear. Of note, DR is a polygenic trait affected by numerous genes, and any polymorphisms within each of these genes will influence the development of DR [ 29 , 30 ]. VEGF is a key molecule in the development of DR [ 7 , 8 ]. Intravitreal VEGF levels are high in diabetic patients and correlate with neovascularization [ 9 ], especially in patients with PDR [ 3 , 10 ]. In addition to VEGF levels, the expression and distribution of the VEGF receptor (VEGFR)-2 are also important in developing DR and progressing from NPDR to PDR [ 11 – 13 ]. Indeed, VEGFR-2 controls endothelial cell proliferation, migration, survival, and osmosis via its activation by VEGF [ 8 ]. Still, other receptors are involved in DR. VEGFR belongs to the receptor tyrosine kinase superfamily and is a membrane mosaic protein. It contains three subtypes: VEGFR-1 (Flt-1), VEGFR-2 (KDR/Flk-1), and VEGFR-3 (Flt-4) [ 31 , 32 ]. In addition, there are some related receptors, such as neuropilin 1 and 2 [ 33 ]. Among them, VEGFR-1 is mainly expressed in monocytes, trophoblast cells, renal mesangial cells, and vascular endothelial cells [ 31 , 32 ]. It is known that VEGF-A, VEGF-B, and PIGF can be combined with it, which is mainly related to the growth regulation of hematopoietic stem cells. VEGFR-2 is mainly distributed on the surface of endothelial cells [ 31 , 32 ]. Other cells include hematopoietic stem cells, retinal precursor cells, and megakaryocytes, which bind to VEGF-A, VEGF-C, VEGF-D, and VEGF-E [ 31 , 32 ]. Combined with VEGF, it stimulates endothelial cell proliferation, increases vascular permeability, and promotes angiogenesis [ 31 , 32 ]. VEGFR-3 is mainly expressed in lymphatic endothelial cells [ 31 , 32 ]. It is known that VEGF-C and VEGF-D combine with VEGFR-3 to regulate the survival and function of lymphatic endothelial cells [ 34 ]. NP-1 and NP-2 mainly exist in nerve axons and vascular endothelial cells. They can regulate the signal transduction of VEGFR-2 [ 35 ]. They are specific receptors of VEGF165.\nAt present, it has been recognized that the effect of VEGF on the proliferation and differentiation of endothelial cells is mediated by VEGFR-2. VEGFR-2 regulates the survival of endothelial cells. VEGF-A depends on VEGFR-2 and its subsequently activated PI3K (phosphatidylinositol 3-kinase) to regulate the survival of human umbilical vein endothelial cells, which is conducive to angiogenesis [ 31 , 32 ]. Second, VEGFR-2 regulates the migration of endothelial cells. The migration of endothelial cells is of great significance to angiogenesis. Endothelial cells can migrate along the concentration gradient of VEGF and other growth factors through the basement membrane hydrolyzed by protease, which plays a key role in metastasis [ 31 , 32 ]. Third, VEGFR-2 regulates the proliferation of endothelial cells. VEGF-A is the VEGF level of many endothelial cells. VEGFR-2 on the cell membrane combines with it to transmit extracellular signals to cells and activate a series of downstream signal pathways to regulate the proliferation of endothelial cells [ 31 , 32 ]. Fourth, VEGFR-2 regulates the permeability of endothelial cells. After VEGFR-2 is activated by VEGF binding, it stimulates vascular endothelial cells to release nitric oxide and causes changes in vascular permeability. The above effects are helpful to promote the survival, proliferation, and migration of endothelial cells and induce neovascularization [ 31 , 32 ].\nBesides  VEGFR-2 , a previous genome-wide study revealed that the rs918519 SNP in  LOC285626 , rs1158314 in  NRXN3 , rs8004963 in  NRXN3 , rs11159428 in  NRXN3 , and rs918520 in  LOC285626  are associated with PDR [ 36 ]. Another study revealed that the -634G>C SNP in  VEGF  is associated with PDR [ 19 ]. Churchill et al. [ 17 ] showed that the -160C, rs13207351, and rs1570360 SNPs in  VEGF  are associated with PDR. So far, no SNPs in  VEGFR-2  have been associated with PDR. A previous study showed that the  VEGFR-2  rs2071559 SNP was associated with DR [ 23 ], but that previous study did not examine the association of the rs2071559 SNP with PDR. The present study examined the rs2071559 SNP with PDR, and it was not associated with DR in all subjects or with PDR among patients with DR. DR is a complex trait affected by the clinical characteristics of the T2DM patient, genetic background, and living environment [ 30 ]. The impact of the rs2071559 on DR can be affected by various other SNPs in the VEGF pathway, as well as in other pathways interacting with the VEGF pathway since DR is a polygenic disorder [ 29 , 30 ]. The relationship between the rs2071559 SNP and the development of PDR in Chinese may be weakened by the influence of other polymorphisms found in the Chinese population, but it will have to be examined. A Spanish study reported no association between the rs2071559 SNP and age-related macular degeneration [ 37 ]. Nevertheless, this SNP has been associated with various conditions, including age-related macular degeneration, stroke, endometriosis, hepatocellular carcinoma, and recurrent spontaneous abortion [ 15 , 16 , 20 , 38 – 40 ], but not with lymphoma [ 41 ]. In patients with T2DM, the rs2071559 SNP is associated with the extent of carotid atherosclerosis and with myocardial infarction [ 42 , 43 ].\nAnti-VEGF drugs have become the first-line treatment for many retinal neovascularization diseases, including age-related macular degeneration and PDR [ 44 , 45 ]. Anti-VEGFR-2 drugs, including the monoclonal antibody ramucirumab and a variety of small-molecule tyrosine kinase inhibitors, remain solely used for the systemic treatment of malignant tumors (e.g., gastric and gastroesophageal adenocarcinoma in the second-line setting) [ 46 ], and are only in early trials for eventual ophthalmologic applications [ 47 , 48 ]. Therefore, in the context of future applications of VEGFR-2-targeting drugs in ophthalmic diseases, understanding the factors related to VEGFR-2 that can influence the course of the diseases is of importance. Future studies should examine all SNPs in  VEGFR-2  and the various proteins involved in its pathway.\nIndeed, reliable SNPs could serve as predictive markers to make informed decisions about the management of specific diseases, the effectiveness of various drugs, and adverse reactions to specific drugs. The plasma sVEGFR-2 levels are elevated in severe DR, while plasma levels of sVEGFR-1 and sVEGFR-3 remain unaltered between patients and controls [ 49 ]. In this study, there were no associations between the  VEGFR-2  rs2071559 SNP and DR of PDR. VEGFR-2 controls endothelial cell proliferation, migration, survival, and osmosis via its activation by VEGF [ 8 ]. VEGFR-2 expression is higher in diabetic eyes than in controls [ 50 ]. The activation of VEGFR-2 by VEGF leads to nitric oxide production, blood vessel permeability, and endothelial proliferation [ 50 ]. An animal study showed that the specific knockdown of VEGFR-2 in the retina led to decreased neovascularization and improved the morphology of the endothelial cells, highlighting the central role of VEGFR-2 in retinal diseases [ 51 ]. Hence, in the context of future anti-VEGFR-2 therapies, a better understanding of the potential effects of VEGFR-2 in DR is essential.\nThis study has limitations. The sample size was small and from a single center. In addition, only one SNP was analyzed in this study because it was not studied before. Other  VEGFR-2  SNPs might yield different results.\nThis study showed no associations between the rs2071559 SNP and DR (among T2DM patients) or PDR (among DR patients). In the Han population in the Shijiazhuang area of China, the rs2071559 SNP of  VEGFR-2  is not associated with the risk of developing DR or PDR.","source_license":"CC-BY-4.0","license_restricted":false}