Molecular Vision 2023; 29:140-152
<http://www.molvis.org/molvis/v29/140>
Received 20 April 2022 |
Accepted 01 August 2023 |
Published 17 August 2023
Tuan Tran,1,2, 3 Paul Sanfilippo,2 Jonathon Goh,2,3 Lyndell L. Lim,2, 3 Sukhpal Sandhu,2,3 Sanjeewa Wickremasinghe2,3
1University of Sydney, Save Sight Institute, Sydney NSW, Australia; 2Centre for Eye Research Australia, Royal Victorian Eye and Ear Hospital, Victoria, Australia; 3Royal Victorian Eye and Ear Hospital, East Melbourne, Victoria, Australia
Correspondence to: Sanjeewa Wickremasinghe, Centre for Eye Research Australia, Level 1, St Andrews Place, East Melbourne. VIC 3002, Australia; Phone: +61 (03) 9929 8360 email: swi@unimelb.edu.au
Purpose: This study aimed to assess the association between cytokine expression and optical coherence tomography (OCT) features of diabetic macular edema (DME).
Methods: This is a post-hoc analysis of the baseline data of a prospective study designed to evaluate the effect of intravitreal ranibizumab injections on aqueous humor cytokine levels in patients with center-involved DME. OCT images were evaluated for baseline imaging features of DME, including central macular thickness (CMT), macular volume (MV), subretinal fluid (SRF) height, morphological pattern of DME, intraretinal cyst (IRC) size, status of the ellipsoid zone/external limiting membrane (EZ/ELM), presence of disorganization of retinal inner layers (DRIL), number of hyperreflective foci (HRF), and the presence of epiretinal membrane (ERM) or vitreomacular traction (VMT). Hierarchical cluster analysis (HCA) was used to identify two groups (clusters) of subjects that shared similar cytokine characteristics (low and high expression profiles), and differences in OCT imaging parameters across clusters were compared using multivariable logistic regression models.
Results: Aqueous cytokine concentrations and OCT images were obtained from 30 eyes of 25 patients recruited. The HCA demonstrated that eyes in the low-concentration cytokine profile group were associated with better BCVA (OR = 0.90, [95% CI 0.81–0.99], p = 0.05), thinner CMT (OR = 1.08, [95% CI 1.01–1.17], p = 0.03), and lower MV (OR = 2.08, [95% CI 1.10–3.90], p = 0.02). Reduced BCVA, greater CMT, and larger MV were individually correlated with multiple cytokines, including IL-7, IL-9, MIP-1α, and TNF-α. Increased SRF was significantly associated with IL-15 (r = 0.40, p = 0.03) and eotaxin (r = 0.47, p = 0.008). Increased levels of IL-6 (H = 39, p = 0.04), IL-13 (H = 44, p = 0.02), and sVEGFR-1 (H = 25.5, p = 0.05) were associated with the presence of DRIL, whereas GM-CSF was protective (H = 131, p = 0.01).
Conclusions: These findings associate higher cytokine expression with worse BCVA, CMT, and MV. Although it is likely that multiple cytokines are associated with the functional and anatomic features of DME, we identified specific cytokines that may drive the process, such as TNF-α, IL-6, IL-7, and IL-9.
Diabetic retinopathy is the leading cause of blindness among working-aged adults with diabetic macular edema (DME), the most common vision-impairing complication. DME is caused by changes in the retinal microvasculature encompassing endothelial cell dysfunction, capillary basement membrane thickening, and a reduced number of pericytes. The resulting vascular permeability and leakage cause fluid accumulation manifested as DME, culminating in retinal dysfunction and vision impairment [1].
There is growing emphasis on the role of inflammation and immune dysregulation in the pathophysiology of DME. Metabolic and oxidative stress from chronic hyperglycemia induces leukostasis and inflammatory breakdown of the blood–retina barrier driven by the coordination of cytokines and growth factors further promoting inflammation [2]. Elevated inflammatory cytokines, such as TNF-α, IL-1β, IL-6 and MCP-1, have been regularly associated with the severity of DME [3-6]. In the DISCERN study [7], we demonstrated changes in the concentrations of multiple cytokines after treatment with intravitreal ranibizumab (IL-1β, IL-7, IL-8, IL-10, IL-12, IL-17, MCP-1, and TNF-α). More recent studies have associated various aqueous cytokines with visual or anatomic responses to treatment with anti-VEGF agents [8-10].
Anti-VEGF-resistant cases of DME classified by worsening or lack of improvement in central macular thickness (CMT) or best-corrected visual acuity (BCVA) can be seen in up to 32%–66% of patients, despite intensive treatment regimes [11-13]. The underlying mechanisms by which these cases differ in their response to anti-VEGF agents is still unclear; however, alternative treatments, such as intravitreal corticosteroids, have proven effective in providing anatomic and functional improvement in such cases [14,15], thereby implicating a significant role of inflammation in DME. Furthermore, intravitreal dexamethasone and triamcinolone have been shown to improve DME by not only reducing VEGF levels but also by targeting the synthesis of proinflammatory mediators involved in DME, such as IL-6, IL-8, MCP-1, ICAM-1, TNF-α, VEGF, and Ang 2 [16-18].
A wide range of phenotypic manifestations of DME are considered imaging biomarkers that may predict treatment response and outcomes. Such OCT features include central macular thickness (CMT), macular volume (MV), subretinal fluid (SRF), morphological patterns, size of intraretinal cysts (IRC), hyper-reflective foci (HRF), disruption of the ellipsoid zone or external limiting membrane, disorganization of retinal inner layers (DRIL), and presence of epiretinal membrane (ERM) [19]. Recent studies have suggested that certain imaging biomarkers may be associated with raised aqueous or vitreous cytokine levels, and thus be contributed by inflammatory components involved in the pathophysiology of DME [20,21]. Analyzing inflammatory cytokines and potential links to imaging biomarkers may improve our understanding of the mechanisms of DME and support individualized therapeutic decision making. The objective of this study was to investigate the association of intraocular cytokines with the imaging features of eyes with DME.
This study was approved by the Human Research and Ethics Committee of the Royal Victorian Eye and Ear Hospital, Melbourne, Australia as part of the DISCERN (Diabetic macular edema: aqueous and Serum Cytokine profiling to determine the Efficacy of RaNibizumab) study and registered with the Australian New Zealand Clinical Trial Registry (ANZCTR identification number 12615000830594). The research adhered to the tenets of the Declaration of Helsinki, and written informed consent was obtained from all participants before enrollment in the study.
This post hoc analysis is based on the cohort of patients enrolled in the DISCERN study. A detailed description of the study design, patient treatment schedule, and aqueous sampling and analysis has been published [7], but in brief, patients with center-involved DME were recruited from the medical retina clinics at the Royal Victorian Eye and Ear Hospital, Melbourne, Australia, between September 2015 and May 2016. All recruited patients had eyes with DME affecting the foveal center and best-corrected visual acuity (BCVA) of 17–72 logarithms of the minimum angle of resolution (LogMAR) letters (Snellen equivalent, 20/400 to 20/40). Although the DISCERN study’s treatment schedule involved treating patients with intravitreal ranibizumab from baseline up to 48 weeks, as per the RESTORE protocol [22], this study only assessed baseline cytokine measurements with baseline OCT characteristics.
Aqueous humor samples (0.05–0.1 ml) were obtained from all patients immediately before intravitreal treatment with ranibizumab. Undiluted samples underwent multiplex immunoassay based on Luminex xMAP technology, and the levels of 32 different cytokines were measured using multiplex detection kits from Bio-Rad (Hercules, CA) and Millipore (Billerica, MA). These cytokines included IL-1b, IL-1ra, IL-2, IL-4, IL-5, IL-6, IL-7, IL-8, IL-9, IL-10, IL-12, IL-13, IL-15, IL-17, eotaxin, granulocyte colony-stimulating factor (G-CSF), granulocyte-macrophage CSF (GM-CSF), monocyte chemoattractant protein-1 (MCP-1), platelet-derived growth factor subunit B, TNF-α, epidermal growth factor (EGF), erythropoietin (EPO), IFN-c, basic FGF (FGF-β), IFN-c-induced protein 10 (IP-10), macrophage inflammatory protein (MIP)-1α, MIP-1β, angopoietin-2 (Ang-2), regulated on activation, normal T-cell expressed and secreted (RANTES), VEGF, soluble VEGF receptor-1 (sVEGFR-1), and soluble VEGF receptor 2 (sVEGFR-2). Standard curves were generated with the reference cytokine sample supplied in the kit using the Bio-Plex200 System (software version 5.1.1; Bio-Rad Laboratories) and were used to calculate the cytokine concentrations (pg/mL) of aqueous humor samples. Each cytokine had a different standard curve and a lower limit of detection on the standard curve of cytokine concentration. Cytokine levels that were below the lower limits of detection were recoded, assigned a numerical value of 0 pg/ml, and included in statistical analysis. We were unable to detect IL-2, RANTES, EPO, and Ang-2 in any sample; therefore, they were excluded from further analysis.
Spectral domain optical coherence tomography (OCT) scans were obtained at each patient visit using Heidelberg Spectralis (Heidelberg Engineering, Heidelberg, Germany). A trained masked grader ensured that all scans were adequately centered on the fovea and were re-centered if necessary. For quantitative and qualitative analysis of OCT images, 49 high-resolution horizontal raster B-scans (spanning 20° × 20° [5.7 mm × 5.7 mm]) centered on the fovea with 120 μm distance between each scan.
The foveal centered SD-OCT scans provided quantitative measurements of the central macular thickness (CMT) and macular volume (MV) within the Early Treatment Diabetic Retinopathy Study (ETDRS), defined as a central zone of 500 μm radius. The Spectralis® software automatically obtained measurements by the in-built “thickness map” function.
For assessing the qualitative features of DME, the horizontal area considered was outlined by the width of horizontal scans (11.4 mm), while the vertical area assessed included five scans above and below the central horizontal scan through the foveola (12 mm height). Each scan was graded by two independent masked retinal specialists (TT and JG), and discrepancies were adjudicated by a third grader (SW).
The morphological pattern of DME was graded as either (1) diffuse retinal thickening (DRT), (2) cystoid macular edema (CME), or (3) serous retinal detachment (SRD), as previously defined by other studies [23,24]. Intraretinal cyst (IRC) size was graded as absent, small, intermediate, or large with the use of reference figures [19], as defined elsewhere [24,25]. The status of the ellipsoid zone (EZ) and the external limiting membrane (ELM) were evaluated within a 1-mm area of the central fovea and graded as either intact or disrupted, as previously defined [19,23,26].
Disorganization of the inner retinal layers (DRIL) was graded as either absent or present [27]. Hyper-reflective foci (HRF) were defined as discrete well circumscribed dots of equal or higher reflectivity than the RPE band. The inclusion and exclusion criteria of HRF were analyzed as defined by Vujosevic et al. [28] HRF were manually counted within the central square area (11.4 × 12 mm) and were graded as follows: ‘0’ = none, ‘1’ = 1–20 and ‘2’ ≥21. Scans were graded as either having the presence or absence of either vitreomacular traction (VMT) or epiretinal membrane (ERM), as described elsewhere [29].
BCVA, CMT, MV, and SRF were measured as continuous variables. pattern of IRC, EZ/ELM integrity, HRF, DRIL, and the presence of ERM/ VMT were considered as categorical variables. For the continuous variables (BCVA, CMT, MV, SRF), Spearman correlation was estimated for each variable with the corresponding individual cytokine, and the results are graphed as scatter plots with a loess curve fitted. For categorical variables (pattern, IRC, EZ/ELM, HRF, ERM/VMT, and DRIL), Mann–Whitney (2 levels) or Kruskal–Wallis test (>2 levels) tests were performed, followed by pairwise Wilcoxon tests using the Benjami & Hochberg (false-discovery rate) adjustment for multiple comparisons.
Hierarchical cluster analysis (HCA) was used to identify two groups (clusters) of subjects that shared similar cytokine characteristics (expression profiles). This was a separate analysis in which the goal of HCA was to build a tree diagram of the cytokine expression profiles, with profiles most similar across subjects placed on branches that were close together, thus forming natural clusters. Cytokine values were first rescaled to z-scores by subtracting the mean cytokine value from each subject’s recorded measurement and then dividing it by the standard deviation of the particular cytokine. The Euclidean distance metric and Ward’s (minimum variance) method were used as the dissimilarity metric and clustering method, respectively. Dendrograms and heatmaps were generated to visualize the clusters and scaled cytokine values. Warmer (more intensely red) colors on the heatmap represent high cytokine values, white colors represent average values, while cooler colors (more intensely blue) represent low values.
Statistical analyses were performed to compare differences in BCVA and OCT imaging parameters (CMT, MV, and SRF) across high or low cytokine expression cluster profiles. Separate multivariable logistic regression models were created with BCVA and each OCT parameter specified (CMT, MV, and SRF) as predictors of each cluster of either high or low expression profiles (adjusted for age and sex). All statistical analyses were conducted in the R statistical environment with HCA performed using the heatmap package (R v4.0.5; R Foundation for Statistical Computing, Vienna, Austria). All tests were two-sided, with a p value <0.05 considered statistically significant.
A total of 30 eyes of 25 patients with center-involving DME were recruited from the DISCERN study [7]. The mean (SD) age at presentation was 63.0 (19.0) years, with a range of 41–82 years, while the mean HbA1c was 7.5% (±1.2). The DR severity according to the ETDRS DRSS ranged from 31 to 81, with 15 eyes (51.7%) being classified as mild to moderate non-proliferative diabetic retinopathy (DRSS score 31–43), 3 eyes (10.3%) were classified as having moderate non-proliferative diabetic retinopathy (DRSS 47), and the remainder had severe non-proliferative diabetic retinopathy or quiescent proliferative DR (12 eyes, 38.0%) and had a DRSS >47.
The mean BCVA was 59.3 letters (± 10.3), mean CMT was 490.6 μm (±134.3), mean MV was 11.2 mm3 (±2.0) and mean SRF was 45 μm (±71.1). Regarding the morphological pattern of DME, 4 eyes had DRT, 14 had CME, and 12 had SRD at baseline. In total, 4 eyes were considered as having small IRC, whereas 9 were graded as intermediate, and 17 were graded as large. Considering the integrity of the EZ/ELM at baseline, 6 patients had intact EZ/ELM, 19 had disrupted EZ/ELM, and 5 had the structure classified as absent. DRIL was present in 7 eyes, while ERM was present in 4 eyes. The presence of HRF was graded as ‘0’ in 2 eyes, ‘1’ in 11 eyes, ‘3’ in 4 eyes, ‘4’ in 3 eyes, and ‘5’ in 5 eyes. Samples with measurable concentrations of cytokines and chemokines were included in the analysis (Table 1).
HCA was performed to classify the patient cohort into two relatively similar groups demonstrating either ‘high’ or ‘low’ levels of cytokine expression. Figure 1 shows the heatmap representing raw data; individuals are segregated into two clusters with each cell colored in proportion to the value of standardized cytokine expression. There was no significant difference between the two cluster groups when adjusted for age and sex when each quantitative OCT parameter was assessed. At baseline, eyes in the high cytokine profile group were associated with worse BCVA, larger CMT, and greater MV. Eyes in the low cytokine profile group had better BCVA (odds ratio, OR = 0.90, [95% confidence interval (CI) 0.81–0.99], p = 0.05), indicating that with one unit worsening of BCVA, there was an average 10% reduction in odds being in the high profile group.
The low cytokine profile group also had thinner CMT (OR = 1.08, [95% CI 1.01–1.17], p = 0.03), indicating that with every 10 μm increase in CMT, there was an average 8% increase in the odds of being in the high profile group. The low cytokine profile group had better MV (OR = 2.08, [95% CI 1.10–3.90], p = 0.02); with every 1.0 mm3 increase in MV, there was an average 108% increase in odds being in the high profile group. There was no significant association with SRF (OR = 1.01, [95% CI 0.99–1.02], p = 0.13). Box plots depicting the significant differences in BCVA, CMT, and MV between low and high cytokine expression profiles are shown in Figure 2.
Worsened BCVA was significantly associated with increased aqueous levels of IL-7 (r = -0.45, p = 0.001), IL-9 (r = -0.41, p = 0.02), MIP-1α (r = -0.4, p = 0.03), and TNF-α (r = 0.33, p = 0.03).
With respect to quantitative OCT features, greater CMT was significantly associated with increased aqueous concentrations of IL-7 (r = 0.4, p = 0.03), IL-9 (r = 0.4, p = 0.03), and TNF-α (r = 0.39, p = 0.03). Similarly, a higher MV was also associated with increased levels of IL-7 (r = 0.37, p = 0.04), IL-9 (r = 0.39, p = 0.03), and TNF-α (r = 0.43, p = 0.02). MV was also associated with other cytokines, including IL-8 (r = 0.44, p = 0.01), IL-10 (r = 0.43, p = 0.02), IL-12 (r = 0.41, p = 0.02), IL-13 (r = 0.38, p = 0.04), IL-15 (r = 0.4, p = 0.03), MIP1α (r = 0.39, p = 0.03), and VEGF (r = 0.43, p = 0.04). The increased height of the SRF was significantly associated with raised aqueous levels of IL-15 (r = 0.40, p = 0.03) and eotaxin (r = 0.47, p = 0.008). Although Spearman’s correlation may be found to be weak, such as with an r = 0.2–0.4 or moderate r = 0.4–0.6, with sufficient sample size, this correlation may be found to be statistically significant (p <0.05). The lower the p-value, the higher the likelihood of the association, whether weak, moderate, or strong. These moderate correlations suggest that the individual cytokine had an impact on the OCT parameter, but other cytokines were also involved. These correlations are shown in Table 2 and depicted in Appendix 1.
The association of aqueous cytokine concentrations and the qualitative features is shown in Table 3, with respective median and interquartile ranges shown in Appendix 2. In keeping with the association of TNF-α with worse BCVA, CMT, and MV at baseline, smaller but more numerous IRC were also found to be associated with greater levels of TNF-α (H = 8.17, p = 0.02), whereas larger IRC were more likely to be observed with lower levels of TNF-α. A similar finding was observed with IFN-γ, with higher concentrations (H = 7.39, p = 0.03) being associated with the presence of smaller and numerous IRC. FGF-β (H = 6.74, p = 0.03) and GM-CSF (H = 6.99, p = 0.03) had more variable associations, with higher concentrations being associated with intermediate size cysts with lower concentrations visible in smaller and larger IRC.
Higher aqueous concentrations of IL-6 (W = 39, p = 0.04), IL-13 (W = 44, p = 0.02), and sVEGFR-1 (W = 25.5, p = 0.05) were more likely to be observed in eyes with DRIL, while higher levels of GM-CSF (W = 131, p = 0.01) were protective. The significant associations are illustrated in Figure 3.
This study demonstrated an association of several inflammatory cytokines with quantitative features of DME, such as BCVA, CMT, MV, and SRF, and qualitative features seen on OCT. From the hierarchical cluster analysis, we found that higher concentrations of inflammation-associated cytokines were associated with worse BCVA, greater CMT, and larger MV at baseline. Specifically, we found that higher concentrations of TNF-α, IL-7, and IL-9 were associated with worse BCVA, greater CMT, and a larger MV. Similarly, higher MIP-1α concentrations were associated with greater MV and worse BCVA but not CMT. Although CMT has evolved as the most commonly used OCT measurement of DME in both research and clinical practice, features greatly affecting CMT, such as large IRC, may not accurately reflect changes in MV, generalized macular thickening, disease severity, or BCVA [30,31]. Thus, MV has been suggested as the most appropriate parameter for assessing the anatomic severity of DME. Accordingly, multiple other cytokines, including IL-8, IL-10, IL-12, IL-13, IL-15, and VEGF, were also positively associated with greater MV. Interestingly, increased aqueous VEGF concentrations were only associated with increased MV, with no correlation with BCVA or other OCT features.
Our study also demonstrated the association of several cytokines with the qualitative OCT features of DME at baseline. Smaller but more numerous intraretinal cysts were associated with greater levels of TNF-α and IFN-γ, whereas increased levels of IL-6, IL-13, and sVEGFR-1 were associated with an increased likelihood of DRIL.
Previous reports have suggested a wide range of aqueous cytokines involved in DME, with a few cytokines such as TNF-α, IL-1β, IL-6 and MCP-1 frequently associated with DME or its features [3,4,20]. However, this has not been consistent across previous studies, which demonstrate conflicting and heterogeneous results with various other cytokines often identified [6,7,32,33].
Most groups performed their analyses by examining a certain range of cytokines from ELISA or multiplex-based assays, which were then individually correlated to clinical features or outcomes [3,8-10,17,18,21,34-43]. However, with generally small numbers of subjects in these studies, it may be difficult to ascertain the involvement or contribution of individual cytokines to a wide range of OCT manifestations or treatment responses. This is because cytokines are involved in complex and unclear mechanisms in DME, where both pro- and anti-inflammatory cytokines are involved in multiple networks and cascades, with regulatory feedback mechanisms functioning in autocrine or paracrine loops [44]. Inflammatory cytokines also exhibit multiple and diverse functions, and their activity may fluctuate according to the pathological state of the disease, whether acute, chronic, worsening, or improving [45]. Thus, it would be helpful to first explore the association of DME with an overview of cytokines by clustering those concordantly expressed in higher concentrations.
The HCA analysis enabled this study to demonstrate that high cytokine expression is associated with DME severity. Performing HCA analysis enables a novel method of clustering inflammatory cytokines, BCVA, CMT, and MVT into relatively similar groups called clusters. The correlation of the high cytokine expression group with the main outcome measures of DME demonstrates the importance of the expression of multiple inflammatory cytokines in the pathophysiological process of DME. The analysis of a large range of cytokines has enabled the identification of multiple inflammatory cytokines associated with disease severity, including IL-7 and IL-9, which have not been previously assessed in DME but have been shown to induce the secretion of TNF-α, IL-1β and IL-6 in bone marrow-derived dendritic cells [46].
TNF-α appears to be a critical element in the pathogenesis of DR. We found that increased aqueous concentrations of TNF-α were not only associated with BCVA, CMT, and MV but were also specifically associated with the presence of microcystic edema. This may relate to the presence of a more acute inflammatory phenotype of DME, as opposed to the larger IRC typically seen in the chronic structural morphology of cystoid degeneration in DME [47,48].
TNF-α is a potent inducer of VEGF, which plays a central role in promoting retinal neovascularization and vascular leakage. TNF-α mediated retinal leukostasis and apoptosis may also contribute to DME by the progressive breakdown in the blood–retina barrier and the increase in retinal endothelial cell permeability via disruption of tight junction complexes [49,50]. Treatment with intravitreal anti-VEGF agents has also been shown to reduce intraocular TNF-α levels, suggesting that the mechanisms of VEGF and inflammatory pathways involving TNF-α are intimately tied in DR [7,51]. The correlation of TNF-α to important quantitative outcome measures of DME, including BCVA, CMT, MV, and IRC, supports its prominent and central role in DME. In the clinical setting, measuring intraocular cytokines may reveal a more inflammatory process or phenotype, thus allowing clinicians to adjust their treatment accordingly.
Our study has limitations that should be acknowledged. First, aqueous sampling was used as a measure of posterior segment pathology. Analysis of vitreous fluid may be more accurate in reflecting DME severity; however, sampling is more invasive and is less preferred by our group and other researchers [3,8,20,52]. Cytokine concentrations found below the lower limit of detection were recoded as zero values and included in the statistical analysis. Although there are various methods of handling these values, with none that are particularly agreed upon [53], they may be of some importance when measuring cytokines or biomarkers in the eye, where concentrations may be relatively low. Recoding these values as zero (instead of as missing values) and including them for statistical analysis may underestimate type II statistical errors, or reduce false positives, which we conceived to be more evaluative for our exploratory analyses.
In correlating cytokines with the OCT features of DME, our study had a small sample size, limiting statistical power. Despite this, the results were statistically significant. We are reassured by the results of the HCA that high expression levels of cytokines are influenced by the main clinical parameters of DME, as assessed based on BCVA, CMT, and MV. Furthermore, multivariate regression analysis found overlap of similar cytokines significantly associated with these parameters. Lastly, analyzing a large number of cytokines raises the possibility of statistical errors related to multiple testing. Although we applied the Benjami & Hochberg method to address this, there are no agreed criteria to apply multiple testing correction methods [54], and multiple comparisons in this study need to be interpreted with caution.
In conclusion, we note that correlating individual cytokines with features of DME is difficult, given the complexity of cytokine functions within inflammatory cascades and the considerable disease overlap of imaging features. The HCA in our study demonstrated that in addition to VEGF, there is a significant contribution of inflammation-associated cytokines to the severity of DME. Specifically, TNF-α, IL-7 and IL-9, appear to be associated with OCT features of worse DME severity. Correlating imaging biomarkers of DME with inflammatory cytokines may allow clinicians to easily assess cases on OCT that portray a greater inflammatory phenotype and thus influence individualized therapeutic decision making.
This work was supported by a grant from Novartis Pharma AG. The Centre for Eye Research Australia receives Operational Infrastructure Support from the Victorian Government.
Medline reports the first author should be "Deák GG" not "Deak GG" in reference 30 "Deak, Schmidt-Erfurth, Jampol, 2018".
Medline reports the first author should be "Yalçın G" not "Yalcin G" in reference 47 "Yalcin, Ozdek, Baran Aksakal, 2019".