Research Paper Volume 13, Issue 6 pp 8369—8379

A dynamic nomogram for predicting diabetic macular edema in type 2 diabetes patients based on plasma cytokines

Ning Zhang1,2, , Jing Ke1,2, , Dawei Zhang3, , Yuanyuan Zhang1,2, , Ying Fu1,2, , Bin Cao1,2, , Dong Zhao1,2, ,

  • 1 Center for Endocrine Metabolism and Immune Diseases, Beijing Luhe Hospital, Capital Medical University, Beijing 101149, China
  • 2 Beijing Key Laboratory of Diabetes Research and Care, Beijing 101149, China
  • 3 Department of Ophthalmology, Beijing Luhe Hospital, Capital Medical University, Beijing 101149, China

Received: November 13, 2020       Accepted: January 22, 2021       Published: March 3, 2021      

https://doi.org/10.18632/aging.202647
How to Cite

Copyright: © 2021 Zhang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Objective: This study investigated changes of plasma cytokines and aimed to build a dynamic nomogram for diabetic macular edema (DME) in type 2 diabetes mellitus (T2DM).

Methods: In a pilot cohort, plasma samples were selected from 9 T2DM patients and 9 DME patients to screen for cytokine differences. The screening cytokines were then validated by enzyme-linked immunoassay in a cohort, which contained 100 DME (DME group) and 100 T2DM patients without DME (T2DM group). A dynamic nomogram for predicting DME was developed, based on the plasma cytokines.

Results: In the pilot cohort, 11 plasma cytokines were significantly increased in the DME group. In the validation cohort, platelet-derived growth factor (PDGF)-BB, tissue inhibitors of metalloproteinase (TIMP)-1, angiopoietin (ANG-1), and vascular endothelial cell growth factor receptor (VEGFR)-2 were confirmed to be significantly elevated in the DME group. The dynamic nomogram demonstrated good calibration and discrimination, with an area under the receiver operating characteristic curve (AUC) of 0.88. In the test set, sensitivity, specificity, and AUC were 73.3%, 80.0%, and 0.84, respectively.

Conclusion: Plasma cytokines were closely associated with DME. A novel dynamic monogram including ANG-1, PDGF-BB, TIMP-1, and VEGFR2 was a novel tool for predicting DME.

Introduction

Diabetic macular edema (DME) is one of the most common microvascular complications in patients with type 2 diabetes mellitus (T2DM) and is a major cause of vision loss worldwide. The prevalence of DME depends on the type and duration of diabetes, and the incidence of DME is higher in individuals with type 1 diabetes than in those with T2DM. Approximately 27% of patients develop DME within 9 years of diabetes onset [1]. Optical coherence tomography (OCT) is an effective and sensitive imaging tool for detecting DME as long as no other causes of cystoid macular edema are present [2]. The fact that patients with macular edema may be asymptomatic and easy to be neglected provides strong support for screening to detect DME [3]. Although the OCT imaging-based features of DME are well-known, knowledge about its protein phenotype is limited.

Angiogenesis and inflammation response are key mechanisms of DME development [4]. ANG-1 expression in preretinal cells activates transmembrane tyrosine kinase Tie2, which is highly expressed in the endothelium and promotes endothelial intercellular links through multiple pathways, thereby reducing endothelial cell permeability [58]. MMP-9/TIMP-1 are critical for maintaining the integrity and impermeability of the blood–retina barrier, and retinal damage occurs when the MMP-9/TIMP-1 ratio is out of balance [9]. PDGF-BB can upregulate VEGF expression and promote angiogenesis after binding to the receptor on the surface of human retinal vascular endothelial cells (hRVECs) [10]. However, data on the association between plasma cytokines and DME is limited.

To explore the potential role of plasma angiogenesis and inflammation cytokines in DME, we investigated the levels of vascular and inflammatory cytokines by protein-chip screening and enzyme-linked immunosorbent assays (ELISAs). A dynamic nomogram for predicting the risk of DME in patients with T2DM was developed and was based on the DME-associated plasma cytokines.

Results

Study subjects

The general characteristics of the study subjects are shown in Table 1. For plasma protein profiling, 9 DME patients and 9 T2DM patients were selected as a pilot cohort. For validation, 100 DME patients and 100 T2DM patients were enrolled. No significant differences in age, body mass index, duration of diabetes, fasting plasma glucose, hemoglobin (Hb)A1c, fasting C-peptide, 2-h post prandial C-peptide, and gender were found between the two groups in either the pilot or the validation cohort. There were also no between-group differences in diabetes-related complications.

Table 1. Comparisons of clinical characteristics of the study cohort.

Clinical characteristicsPilot cohortValidation cohort
T2DM (n=9)
(Mean ± SD)
DME (n=9)
(Mean± SD)
pT2DM (n=100)
(Mean ± SD)
DME (n=100)
(Mean ± SD)
p
Age
(years)
59.44±8.8261.00±10.480.73862.41±12.4359.76±11.860.141
BMI
(Kg/m2)
26.17±1.7426.36±3.890.89525.48±3.8225.80±3.600.573
Duration of diabetes (years)6.56±9.658.78±6.850.6827.40±8.637.48±7.310.951
Fasting plasma glucose (mmol/L)6.58±2.386.89±1.680.7588.73 ±3.268.26 ±4.180.437
HbA1c
(%)
10.2±2.349.66±1.790.5739.90 ±2.029.52 ±2.550.269
Fasting C peptide (mIU/L)1.47±0.661.63±1.060.6971.73 ±1.181.79 ±1.190.737
2-h post prandial C-peptide (mIU/L)5.80±4.303.68±2.010.2053.96 ± 3.284.60 ± 2.760.183
Triglyceride
(mmol/L)
2.09±1.442.05±1.910.9631.96 ± 2.321.65 ± 0.950.279
Total cholesterol (mmol/L)5.04±2.694.24±1.500.4814.37 ± 1.394.38 ± 1.030.933
Low-density lipoprotein
(mmol/L)
3.27±1.912.48±1.140.3342.77 ± 1.072.85 ± 0.840.606
Gender,
male (%)
5(55.6%)4 (44.4%)1.00053 (53.0%)46 (46.0%)0.396
Hypertension,
number (%)
4 (44.4%)6(66.7%)1.00050 (50.0%)44 (44.0%)0.479
*Diabetic nephropathy,
number (%)
3 (33.3%)5 (55.6%)0.63540 (42.6%)29 (31.5%)0.160
Diabetic peripheral neuropathy,
number (%)
0 (0%)0 (0%)1.0001 (1.0%)1 (1.0%)1.000
*, 14 missing data of diabetic nephropathy in validation cohort.

Multiple cytokine alterations in DME plasma

To profile plasma cytokines, blinded screening by protein microarray analysis was performed, and semi-quantitative results were obtained for 60 plasma proteins. The relative changes of plasma cytokines are shown in the heatmap in Figure 1A. Compared with the T2DM group, the levels of 15 cytokines were significantly lower in the DME group and 11 were significant higher, with more than a four-fold change (adjusted p < 0.05, Figure 1B). The increases in 5 plasma proteins, PDGF-BB, TIMP-1, ANG-1, CXCL16, and VEGFR2 in the DME group were greater than six-fold. Principal component analysis (PCA) found a clear distinction between the two groups, suggesting that these 5 plasma proteins might be helpful to distinguish T2DM patients with and without DME (Supplementary Figure 1).

Semi-quantitative results of cytokine analysis in the pilot cohort. Heatmaps of the relative changes of 60 plasma cytokines (A) and 26 cytokines with a fold change > 4 or B).

Figure 1. Semi-quantitative results of cytokine analysis in the pilot cohort. Heatmaps of the relative changes of 60 plasma cytokines (A) and 26 cytokines with a fold change > 4 or < 0.25 (B).

ELISA validation of five markedly increased cytokines

To further determine whether ANG-1, CXCL16, PDGF-BB, TIMP-1, and VEGFR2 levels were significantly increased in DME, 100 DME patients and 100 T2DM patients were recruited to validate the array data using ELISA. As shown in Figure 2, the mean levels of ANG-1 (334.62 ng/ml), PDGF-BB (37.86 pg/ml), TIMP-1 (7.39 ng/ml), and VEGFR2 (12.47 ng/ml) in the DME group were significantly higher than those in the T2DM group (268.94 ng/ml, 28.60 pg/ml, 6.27 ng/ml, and 10.71 ng/ml), respectively (p < 0.001). Between-group differences in CXCL16 (3776.43 pg/ml vs. 3794.39 pg/ml) were not significant (p = 0.791).

Comparison of the plasma concentrations of ANG-1, CXCL16, PDGF-BB, TIMP-1, and VEGFR2 in the validation cohort. ANG-1, PDGF-BB, TIMP-1 and VEGFR2 concentrations were significantly increased in DME samples. The between-group difference of CXCL16 concentration was not significant.

Figure 2. Comparison of the plasma concentrations of ANG-1, CXCL16, PDGF-BB, TIMP-1, and VEGFR2 in the validation cohort. ANG-1, PDGF-BB, TIMP-1 and VEGFR2 concentrations were significantly increased in DME samples. The between-group difference of CXCL16 concentration was not significant.

Correlation of cytokine levels and clinical characteristics in DME group

As shown in the heatmap (Figure 3A), Pearson correlation analysis of the relationships of cytokine levels and clinical features found no significant positive correlations between clinical characteristics, including age, body mass index, waist-to-hip ratio, diabetes duration, fasting blood glucose, HbA1c, fasting C-peptide, 2-h postprandial C-peptide, and the 5 plasma cytokines, including ANG-1, PDGF-BB, TIMP-1, VEGFR2 and CXCL16 (r < 0.3). Focusing on the inner relationship of cytokines, there was also no obvious correlation among these 5 cytokines (r < 0.3, Figure 3B).

Pearson correlation analysis of DME group. There were no significant positive correlations between plasma cytokines and clinical characteristics (A) among plasma cytokines (B).

Figure 3. Pearson correlation analysis of DME group. There were no significant positive correlations between plasma cytokines and clinical characteristics (A) among plasma cytokines (B).

PCA for cytokine selection

PCA was performed to determine the relative contribution of each cytokine to the ability to distinguish DME- and T2DM-group patients. The first and second principal components of the PCA plot (Dim1 and Dim2) accounted for 41.7%, and 20.3% of the variation, respectively, in the dataset. The projection of samples in the PCA revealed relatively little overlapping of areas. CXCL16 contributed more to the second than to the first principal component. ANG-1, PDGF-BB, TIMP-1, and VEGFR2 contributed more to the first principal component (Figure 4A). As shown in Figure 4B, the percentage contributions of cytokines to first principal component were ANG-1 (31.80%), PDGF-BB (27.99%), TIMP-1 (24.52%), VEGFR2 (14.43%), and CXCL16 (1.26%). According to contribution percent, ANG-1, PDGF-BB, TIMP-1, and VEGFR2 contribute more to distinguish T2DM patients with and without DME. To reduce overfitting, ANG-1, PDGF-BB, TIMP-1, and VEGFR2 were included in a logistic regression analysis.

Principal component analysis (PCA) for feature selection and computation of the relative contribution of each cytokine to distinguishing DME and T2DM (A) and the contribution of each cytokines to the first principal component (B).

Figure 4. Principal component analysis (PCA) for feature selection and computation of the relative contribution of each cytokine to distinguishing DME and T2DM (A) and the contribution of each cytokines to the first principal component (B).

Development and validation of logistic regression and nomogram development

Multivariate logistic regression was used to develop a model to predict DME. As shown in Figure 5A, ANG-1 (odds ratio [OR] = 1.01, 95% confidence interval [CI]: 1.00–1.02; p = 0.026], PDGF-BB (OR = 1.16, 95% CI: 1.08–1.24, p<0.001), and TIMP-1 (OR = 1.76, 95% CI: 1.16–2.66, p = 0.008) were identified as significant risk factors. VEGFR (p =0.061) did not reach statistical significance. However, as Akaike information criterion (AIC) increased from 129.65 to 131.32 when removing VEGFR2 from the model, it was retained in the analysis. Figure 5B shows that the logistic regression model performed well in the training dataset, with an area under the receiver operating characteristic curve (AUC) of 0.88. In the test set, the sensitivity was 73.3%, specificity was 80.0% and the AUC was 0.84.

Forest plot and receiver operating characteristic curve of logistic regression. ANG-1, PDGF-BB, and TIMP-1 were identified as risk factors. VEGFR2 did not reach statistical significance, (p-value slightly larger than 0.05). Because the Akaike Information Criterion increased from 129.65 to 131.32 when removing VEGFR2 from the model, it was retained (A). Receiver operating characteristic curve for the prediction model. The area under the receiver operating characteristic curve of the training set (red line) and the test set (olive drab line) were 0.88 and 0.84, respectively (B).

Figure 5. Forest plot and receiver operating characteristic curve of logistic regression. ANG-1, PDGF-BB, and TIMP-1 were identified as risk factors. VEGFR2 did not reach statistical significance, (p-value slightly larger than 0.05). Because the Akaike Information Criterion increased from 129.65 to 131.32 when removing VEGFR2 from the model, it was retained (A). Receiver operating characteristic curve for the prediction model. The area under the receiver operating characteristic curve of the training set (red line) and the test set (olive drab line) were 0.88 and 0.84, respectively (B).

Based on the results of multivariate logistic regression, a nomogram including ANG-1, PDGF-BB, TIMP-1, and VEGFR2 was constructed for DME risk prediction (Figure 6A). The point score of each variable was based on its contribution to the model and the point total corresponds to the risk and predicted likelihood of DME. Calibration curves demonstrated good consistency between the predicted risk and the actual probability (Figure 6B). The apparent curve confirmed the good prediction capability of the nomogram.

Nomogram and calibration curve for predicting DME. The DME-risk nomogram includes ANG-1, PDGF-BB, TIMP-1, and VEGFR2 as predictors (A). The calibration curve for predicting DME (B). The nomogram-predicted probability of DME is plotted on the X-axis and the actual probability is plotted on the Y-axis. Calibration curves demonstrated satisfactory consistency between the predicted risk and the actual probability.

Figure 6. Nomogram and calibration curve for predicting DME. The DME-risk nomogram includes ANG-1, PDGF-BB, TIMP-1, and VEGFR2 as predictors (A). The calibration curve for predicting DME (B). The nomogram-predicted probability of DME is plotted on the X-axis and the actual probability is plotted on the Y-axis. Calibration curves demonstrated satisfactory consistency between the predicted risk and the actual probability.

Clinical usefulness of the nomogram

The decision curve demonstrated that in both the training and test datasets, the net benefit was greater than with an all-or-none patient intervention scheme if the threshold probability was less than 60%, which supports use of the nomogram in clinical practice (Figure 7A). A clinical impact curve was plotted to predict the number of high-risk patients in a population of 1000. In both the training and test datasets, the predicted high-risk number was close to the actual number of event cases when the threshold probability was greater than 0.3, and the cost-benefit ratio was close to 2:5 (Figures 7B, 7C). The dynamic nomogram is accessible online as user-friendly digital interface (https://doctorcao.shinyapps.io/DynNomapp/).

Decision curve analysis (DCA) and clinical impact plot of the nomograms for DME in both the training and test sets. The net benefit is greater than that achieved with an all-or-none patient intervention scheme if the threshold probability is less than 60% (A). Predictive high-risk estimates were similar to the actual number of event cases when the threshold probability was > 0.3 in both the training (B) and test set (C).

Figure 7. Decision curve analysis (DCA) and clinical impact plot of the nomograms for DME in both the training and test sets. The net benefit is greater than that achieved with an all-or-none patient intervention scheme if the threshold probability is less than 60% (A). Predictive high-risk estimates were similar to the actual number of event cases when the threshold probability was > 0.3 in both the training (B) and test set (C).

Discussion

DME is a type of retinal vascular disease. Angiogenic and inflammatory processes mediate the genesis of DME, and it is important to define which processes are enhanced or decreased in clinical practice. In the pilot cohort, we analyzed the expression changes of 60 angiogenic and inflammatory factors in DME by protein-chip screening. Eleven cytokines, including PDGF-BB, TIMP-1, ANG-1, and VEGFR2 were significantly increased in DME compared with those T2DM group. We expanded the sample size in the validation cohort and obtained consistent results that were used to develop a dynamic online nomogram for DME prediction. The user-friendly digital interface makes it easier to predict risk and to make better clinical decisions.

In this study, we highlighted the significance of increased expression of inflammatory and angiogenic factors in DME. Studies found that cytokines were elevated or decreased in the DME group and were closely related to the pathogenesis of vascular permeability and inflammation during the development of DME [1113]. Our study confirmed that ANG-1, VEGFR2, TIMP1, and PDGF-BB, which were generally correlated with inflammatory process, were increased in DME. However, the pathological mechanisms of these four plasma cytokines during the progression of DME is unclear, and further studies are needed.

In this study, the characteristics of elevated vascular inflammatory factors in DME were used to establish a prediction model for DME. Previous clinical studies have reported elevated VEGF, MCP1, IL-6, and ICAM-1 were detected in vitreous fluids of DME patients, which is consistent with our results [1422]. There are few published studies on the joint utilization of cytokines to predict DME. In the investigation of the pathogenesis of DME, we screened 60 angiogenic and inflammatory cytokines and finally selected four plasma cytokines for joint modeling.

The use of OCT-imaging biomarkers for evaluating therapeutic response has been described. Ceravolo et al. [23] demonstrated that the number of hyper-reflective spots (HRSs) and serous detachment of the neuroepithelium were specific noninvasive biomarkers for monitoring treatment response. Vujosevic et al. [24] reported that HRS number, extent of disorganization of the inner retinal layers, central macular thickness, and the cyst area of the deep capillary plexus were retinal biomarkers that may help in evaluating treatment response. Few studies have focused on the use of plasma proteins in the diagnosis of DME. In this study, we confirmed that plasma cytokines not only help to evaluate the risk of DME with good sensitivity and specificity, but may also provide new insights into the pathological mechanism of DME. The results have clinical significance for risk assessment and diagnosis when the use of OCT imaging is limited or not universally available.

We did not observe positive correlations between cytokines and indexes of glucose metabolism such as HbA1c or C-peptide in the DME group. The possible cause was that the study subjects were all hospital inpatients in the department of endocrinology for improved glucose management and control, which may have eliminated the effects of glucose level on the results. The association of cytokines with the course of disease or diabetic complications needs further investigation.

There were some limitations. First, we included about 200 subjects for model development and verification. Although the results were clear, there was no external cohort to validate the performance of the model. A large prospective study is needed for further validation. Second, coronary heart disease (CHD) may also result in abnormal expression of plasma proteins. However, the diagnosis of CHD was based on a history of disease, and as CHD in diabetes patients may be asymptomatic, a few CHD patients may have been included in this study.

In conclusion, plasma ANG-1, TIMP-1, PDGF-BB, and VEGFR2 were increased and were comprehensive predicators of DME when included in a dynamic nomogram. The nomogram needs further confirmation in large populations.

Materials and Methods

Subjects

The study received ethical approval from the competent Institutional Review Boards of Beijing Luhe Hospital. All procedures complied with the ethical principles of Helsinki Declaration for studies of human subjects. The study was registered on clinicaltrials.gov (NCT03811470).

Patients were recruited at the Center for Endocrine Metabolism and Immune Diseases of Beijing Luhe Hospital, Capital Medical University, Beijing, China from January 2017 to December 2018. All patients participated in the program of the National Metabolic Management Center (MMC) [25]. 9 T2DM patients and 9 DME patients were recruited for a pilot study in which a protein antibody array was used to screen aberrantly expressed plasma protein. 100 DME patients (DME group) and 100 T2DM patients without DME or diabetic retinopathy (T2DM group) were recruited for a validation study using ELISA kits to verify the aberrantly expressed plasma proteins that were identified in the pilot study. T2DM patients with and without DME were eligible for inclusion. Patients with T1DM or other type of DM, with other types of retinopathy, with a history of any previous intravitreal injection or any other treatment for DME, or a history of cardiovascular diseases or stroke were excluded. The diagnosis of DME was determined by OCT [26], and confirmed by two senior ophthalmologists.

Plasma sample collection

For the antibody arrays and ELISA, blood samples were collected at our center with ethylenediaminetetraacetate as the anticoagulant. The plasma was isolated by standard blood processing and the aliquots were frozen and stored at −80° C, avoiding freeze-thaw cycles.

Cytokine antibody array screening

Cytokines were assayed in duplicate with RayBio R Human Cytokine Antibody Array (RayBiotech G-Series Human Angiogenesis Array 2 and G-Series Human Angiogenesis Array 3, RayBiotech). Spiking and recovery tests were performed to ensure a linear concentration response. As determined by densitometry, the interarray coefficient of variation of spot signal intensities was 20%. Cytokine results were presented as a heatmap using the “pheatmap” package.

ELISA validation

Plasma PDGF-BB, TIMP-1, ANG-1, CXCL16, VEGFR2 levels were determined with human ELISA kits (Mlbio) following the manufacturer’s instructions. The intra-assay coefficient of variation was 20%, and the intra-assay coefficient of variation was 12%. No significant cross reactivity or interference was observed.

Development of a dynamic nomogram

The patients in the validation cohort were randomly divided into training and test sets, which included 70% and 30% of the data, respectively. The cytokines most predictive of DME were selected by PCA using the data included in the training cohort. Multivariate logistic regression identified the cytokines significantly associated with DME development and they were used to construct a nomogram presenting a specific system for calculating the risk of DME. The performance of nomogram was evaluated by sensitivity, specificity, discrimination, and calibration. Discrimination, which is the predictive accuracy to distinguish patients with DME from those without DME, was measured by the receiver operating characteristic (AUC). A calibration curve, which reflected the consistency between the predicted probability and the actual probability, was plotted using 1000 bootstrap resamples. Decision curve analysis (DCA) was used to assess clinical usefulness of the nomogram. For access to risk estimation, a dynamic nomogram having a user-friendly digital interface was created online using the DynNom package [27].

Statistical analysis

All data analysis and visualization were performed with in R software, version 3.6.3 (The R Project for Statistical Computing). Significance was evaluated by t-tests for normally distributed data; otherwise, the nonparametric Mann-Whitney test was used to analyze the data. Chi-square tests were performed for categorical variables. Pearson’s correlation analysis was performed to assess the relationships of the plasma cytokines. In all cases, p-values < 0.05 were considered statistically significant.

Ethical approval

The study received ethical approval from the competent Institutional Review Boards of Beijing Luhe Hospital.

Supplementary Materials

Supplementary Figures

Author Contributions

The study was conceived and designed the experiments by Bin Cao and Dong Zhao. Ning Zhang, Ying Fu, Jing Ke, Dawei Zhang and Yuanyuan Zhang performed the experiments. Bin Cao analyzed the data and together with Ning Zhang interpreted the results and wrote the manuscript. Dong Zhao supervised the statistical computations, interpreted the results, commented on successive drafts of the manuscript. Dong Zhao approved the final version.

Acknowledgments

The authors would like to thank all the participants who took part in the study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Funding

This study was funded by Science and Technology Committee of Tongzhou District (Number: KJ2019CX014-25, KJ2020CX001-11) and the National Natural Science Foundation of China (81800723).

References

  • 1. Li JQ, Welchowski T, Schmid M, Letow J, Wolpers C, Pascual-Camps I, Holz FG, Finger RP. Prevalence, incidence and future projection of diabetic eye disease in Europe: a systematic review and meta-analysis. Eur J Epidemiol. 2020; 35:11–23. https://doi.org/10.1007/s10654-019-00560-z [PubMed]
  • 2. Flaxel CJ, Adelman RA, Bailey ST, Fawzi A, Lim JI, Vemulakonda GA, Ying GS. Diabetic retinopathy preferred practice pattern®. Ophthalmology. 2020; 127:P66–145. https://doi.org/10.1016/j.ophtha.2019.09.025 [PubMed]
  • 3. American Diabetes Association. 11. Microvascular complications and foot care: Standards of medical care in diabetes-2020. Diabetes Care. 2020 (Suppl 1); 43:S135–51. https://doi.org/10.2337/dc20-S011 [PubMed]
  • 4. Romero-Aroca P, Baget-Bernaldiz M, Pareja-Rios A, Lopez-Galvez M, Navarro-Gil R, Verges R. Diabetic Macular Edema Pathophysiology: Vasogenic versus Inflammatory. J Diabetes Res. 2016; 2016:2156273. https://doi.org/10.1155/2016/2156273 [PubMed]
  • 5. Thurston G, Rudge JS, Ioffe E, Zhou H, Ross L, Croll SD, Glazer N, Holash J, McDonald DM, Yancopoulos GD. Angiopoietin-1 protects the adult vasculature against plasma leakage. Nat Med. 2000; 6:460–63. https://doi.org/10.1038/74725 [PubMed]
  • 6. Gavard J, Patel V, Gutkind JS. Angiopoietin-1 prevents VEGF-induced endothelial permeability by sequestering Src through mDia. Dev Cell. 2008; 14:25–36. https://doi.org/10.1016/j.devcel.2007.10.019 [PubMed]
  • 7. Parikh SM. Angiopoietins and Tie2 in vascular inflammation. Curr Opin Hematol. 2017; 24:432–38. https://doi.org/10.1097/MOH.0000000000000361 [PubMed]
  • 8. Takagi H, Koyama S, Seike H, Oh H, Otani A, Matsumura M, Honda Y. Potential role of the angiopoietin/tie2 system in ischemia-induced retinal neovascularization. Invest Ophthalmol Vis Sci. 2003; 44:393–402. https://doi.org/10.1167/iovs.02-0276 [PubMed]
  • 9. Ziablitsev SV, Korobova AV, Petrenko OV, Serduk VN, Mogilevsky SU. Role of matrix metalloproteinase 9 and its tissue inhibitor 1 in development and prognosis of diabetic retinopathy. Fiziol Zh. 2016; 62:37–44. https://doi.org/10.15407/fz62.05.037 [PubMed]
  • 10. Kernt M, Neubauer AS, Liegl RG, Hirneiss C, Alge CS, Wolf A, Ulbig MW, Kampik A. Sorafenib prevents human retinal pigment epithelium cells from light-induced overexpression of VEGF, PDGF and PlGF. Br J Ophthalmol. 2010; 94:1533–39. https://doi.org/10.1136/bjo.2010.182162 [PubMed]
  • 11. Das A, McGuire PG, Rangasamy S. Diabetic macular edema: pathophysiology and novel therapeutic targets. Ophthalmology. 2015; 122:1375–94. https://doi.org/10.1016/j.ophtha.2015.03.024 [PubMed]
  • 12. Cohen SR, Gardner TW. Diabetic retinopathy and diabetic macular edema. Dev Ophthalmol. 2016; 55:137–46. https://doi.org/10.1159/000438970 [PubMed]
  • 13. Mesquida M, Drawnel F, Fauser S. The role of inflammation in diabetic eye disease. Semin Immunopathol. 2019; 41:427–45. https://doi.org/10.1007/s00281-019-00750-7 [PubMed]
  • 14. Lally DR, Shah CP, Heier JS. Vascular endothelial growth factor and diabetic macular edema. Surv Ophthalmol. 2016; 61:759–68. https://doi.org/10.1016/j.survophthal.2016.03.010 [PubMed]
  • 15. Dong N, Xu B, Chu L, Tang X. Study of 27 aqueous humor cytokines in type 2 diabetic patients with or without macular edema. PLoS One. 2015; 10:e0125329. https://doi.org/10.1371/journal.pone.0125329 [PubMed]
  • 16. Funatsu H, Noma H, Mimura T, Eguchi S, Hori S. Association of vitreous inflammatory factors with diabetic macular edema. Ophthalmology. 2009; 116:73–79. https://doi.org/10.1016/j.ophtha.2008.09.037 [PubMed]
  • 17. Chen H, Zhang X, Liao N, Wen F. Assessment of biomarkers using multiplex assays in aqueous humor of patients with diabetic retinopathy. BMC Ophthalmol. 2017; 17:176. https://doi.org/10.1186/s12886-017-0572-6 [PubMed]
  • 18. Noma H, Yasuda K, Shimura M. Cytokines and the Pathogenesis of Macular Edema in Branch Retinal Vein Occlusion. J Ophthalmol. 2019; 2019:5185128. https://doi.org/10.1155/2019/5185128 [PubMed]
  • 19. Chung YR, Kim YH, Lee SY, Byeon HE, Lee K. Insights into the pathogenesis of cystoid macular edema: leukostasis and related cytokines. Int J Ophthalmol. 2019; 12:1202–08. https://doi.org/10.18240/ijo.2019.07.23 [PubMed]
  • 20. Urias EA, Urias GA, Monickaraj F, McGuire P, Das A. Novel therapeutic targets in diabetic macular edema: beyond VEGF. Vision Res. 2017; 139:221–27. https://doi.org/10.1016/j.visres.2017.06.015 [PubMed]
  • 21. Klaassen I, Van Noorden CJ, Schlingemann RO. Molecular basis of the inner blood-retinal barrier and its breakdown in diabetic macular edema and other pathological conditions. Prog Retin Eye Res. 2013; 34:19–48. https://doi.org/10.1016/j.preteyeres.2013.02.001 [PubMed]
  • 22. Matsushima R, Noma H, Yasuda K, Goto H, Shimura M. Role of cytokines in ranibizumab therapy for macular edema in patients with central retinal vein occlusion. J Ocul Pharmacol Ther. 2019; 35:407–12. https://doi.org/10.1089/jop.2019.0011 [PubMed]
  • 23. Ceravolo I, Oliverio GW, Alibrandi A, Bhatti A, Trombetta L, Rejdak R, Toro MD, Trombetta CJ. The Application of Structural Retinal Biomarkers to Evaluate the Effect of Intravitreal Ranibizumab and Dexamethasone Intravitreal Implant on Treatment of Diabetic Macular Edema. Diagnostics (Basel). 2020; 10:413. https://doi.org/10.3390/diagnostics10060413 [PubMed]
  • 24. Vujosevic S, Toma C, Villani E, Muraca A, Torti E, Florimbi G, Leporati F, Brambilla M, Nucci P, De Cilla’ S. Diabetic macular edema with neuroretinal detachment: OCT and OCT-angiography biomarkers of treatment response to anti-VEGF and steroids. Acta Diabetol. 2020; 57:287–96. https://doi.org/10.1007/s00592-019-01424-4 [PubMed]
  • 25. Zhang Y, Wang W, Ning G. Metabolic management center: an innovation project for the management of metabolic diseases and complications in China. J Diabetes. 2019; 11:11–13. https://doi.org/10.1111/1753-0407.12847 [PubMed]
  • 26. Treatment techniques and clinical guidelines for photocoagulation of diabetic macular edema. Early treatment diabetic retinopathy study report number 2. Early treatment diabetic retinopathy study research group. Ophthalmology. 1987; 94:761–74. https://doi.org/10.1016/s0161-6420(87)33527-4 [PubMed]
  • 27. Jalali A, Alvarez-Iglesias A, Roshan D, Newell J. Visualising statistical models using dynamic nomograms. PLoS One. 2019; 14:e0225253. https://doi.org/10.1371/journal.pone.0225253 [PubMed]