Relating Retinal Vascular Oxygen Saturation and Microvasculature Morphology at Progressive Stages of Diabetic Retinopathy.

Purpose: Diabetic retinopathy (DR) is a common cause of vision loss in working age adults and presents changes in retinal vessel oxygenation and morphology. The purpose of this study was to test the hypothesis that there is an association of retinal vessel oxygen saturation with vessel density (VD) and tortuosity in DR.
Methods: Ninety-five subjects were classified in the following groups: nondiabetic control (N = 25), no DR (N = 28), mild nonproliferative DR (NPDR; N = 21), moderate to severe NPDR (N = 14), or treated proliferative DR (PDR; N = 7). Retinal oximetry was performed to measure arterial and venous oxygen saturation (SO2A and SO2V) and calculate oxygen extraction fraction (OEF). Optical coherence tomography angiography (OCTA) was performed for measurements of VD and vessel tortuosity index (VTI).
Results: There were statistically significant differences in SO2A and SO2V among groups (P ≤ 0.004). SO2A and SO2V were higher in the PDR group compared to the control group and SO2V was also higher in the moderate to severe NPDR group. VD differed significantly among groups (P = 0.003), whereas VTI was not significantly different (P = 0.22). Compared to the control group, VD was lower in moderate to severe NPDR and PDR groups. VD was also lower in the PDR group than that in the no DR group (P = 0.03). There was a significant correlation of VTI with SO2V (r = 0.32, P = 0.002) and OEF (r = -0.35, P = 0.001). Conclusions: Retinal vessel morphology, oxygenation, and tissue oxygen extraction were associated with each other in a cohort of subjects with and without DR. Translational Relevance: The findings of this study have the potential to improve clinical management of DR by providing better understanding of human disease pathophysiology and propelling future studies to identify multiple image-based biomarkers for improved disease diagnosis and monitoring.

Introduction

Among diabetes complications, diabetic retinopathy (DR) has a significant impact on vision and can eventually lead to blindness. In fact, DR is one of the most common causes of visual loss in working age adults.– An underlying cause of vision-threatening DR pathologies is thought to be retinal hypoxia due to nonperfusion of the retinal capillaries during the progression of the disease.– Currently, there are no clinical methods available to measure retinal tissue oxygen content. However, previous studies have shown increased retinal arterial and venous oxygen saturation (SO2A and SO2V) at stages of DR.– Other studies have reported decreased difference between SO2A and SO2V or oxygen extraction fraction (OEF) in subjects with non-proliferative DR (NPDR), mild or no DR.–,– Some studies have attributed the decrease in oxygen extraction to changes in the retinal vessel endothelium and regulatory function of pericytes.–

Morphological changes in retinal vessels have been shown to be present in DR. Retinal vessel density (VD) was reported to be decreased with progression of DR, primarily due to capillary nonperfusion.– Additionally, tortuosity of retinal vessels was found to increase with progression of DR,, and eventually decrease in the severe proliferative stage. Although the mechanism of increased tortuosity is not entirely clear, studies suggested that it may be caused by damage to endothelial and autoregulatory cells, limiting the ability of the vessel walls to adjust to changes in blood flow,,, under normal physiology or hypoxic conditions.

Alterations in retinal VD and vessel wall elasticity may impair the delivery to and uptake of oxygen by the tissue, thereby affecting vascular oxygen content and tissue oxygen extraction. The purpose of the current study was to test the hypothesis that there is an association between retinal vessel oxygen saturation and tortuosity in DR.

Methods

Subjects

This study was approved by an institutional review board of the University of Southern California (USC). Subjects were recruited from the USC retina clinic. Before subjects were enrolled, the research study was explained, and informed consents were obtained according to the tenets of the Declaration of Helsinki. A total of 95 subjects participated in the study. Based on retina examination by experienced retina specialists, subjects were classified as diabetic with no DR (N = 28), mild (N = 21), or moderate to severe (N = 14) with or without macular edema, proliferative DR (PDR; N = 7), or healthy nondiabetic control (N = 25). Subjects with PDR had received no treatment (N = 1), antivascular endothelium growth factor treatment (N = 1), panretinal photocoagulation treatment (N = 3), or both treatments (N = 2). Nondiabetic subjects with a history of any eye disease and subjects with diabetes with a history of eye diseases, other than DR, were excluded from this study.

Retinal findings were documented using our previously established retinal oximetry system, and commercially available optical coherence tomography angiography (OCTA) instrument (Optovue Inc., Fremont, CA). Hematocrit (HCT) and hemoglobin A1C (HbA1c) were measured from blood samples obtained from a finger-prick and using a micro-hematocrit centrifuge (Unico, Dayton, NJ) and HbA1C analyzer (Siemens DCA Vantage Analyzer, Frimley, Camberley, UK). HCT and HbA1c data were not available in one subject in the no DR group and two subjects in the control group, respectively. Blood pressure (BP) was recorded using a wrist cuff to obtain the mean arterial pressure (MAP) = 1/3 (systolic BP) + 2/3 (diastolic BP). Three MAP measurements were averaged.

Oxygen Saturation and Extraction Fraction

Using our previously established retinal oximetry system,, 9 retinal reflectance images were acquired at oxygen-sensitive and oxygen-insensitive imaging wavelengths of 606 and 570 nm, within 3 seconds. Images acquired at each wavelength were registered using the stack-reg plug-in of ImageJ (National Institutes of Health, Bethesda, MD) to compensate for eye motion movements and then averaged. Retinal vessels within a circumpapillary region of interest (Fig. A) were automatically segmented and centerlines were generated by our previously published image analysis algorithm., From the mean image at each wavelength, optical densities (ODs) were calculated for each vessel as log (I/I), where I and I represent the average pixel intensity inside and outside the vessel, respectively. OD ratios (OD/OD), were adjusted and converted to SO2 values using a linear regression model, as previously described. Measurements in individual arteries and veins were averaged separately to calculate SO2A and SO2V per eye. OEF was calculated as: OEF = (SO2A – SO2V) / SO2A. OEF is the fraction of the total oxygen supplied by the retinal circulation that is extracted by the tissue. It is also defined by the ratio of oxygen metabolism to delivery, and thus independent of blood flow.

Figure.

Examples of retinal images obtained in a (A) control subject and (B) moderate to severe nonproliferative diabetic retinopathy subject. Retinal vascular oxygen saturation (SO2) measurements in major arteries and veins are depicted in pseudo-color within a circumpapillary region bounded by green circles. Examples of optical coherence tomography angiography images acquired in the same (C) control subject and (D) moderate to severe nonproliferative diabetic retinopathy subject. Vessel centerlines (red lines) are overlaid on the vessels and vessel end points are indicated by yellow circles. Vessel tortuosity index (VTI) was determined by averaging measurements from all vessels.

Vessel Density and Tortuosity

Using a commercially available OCTA system (Avanti Optovue, Inc., Fremont, CA), images were acquired at a wavelength of 840 nm with an axial scan rate and axial scan depth resolution of 70 kHz and 5 µm, respectively. The instrument generated OCTA images in a 6 mm × 6 mm perifoveal region based on detection of red blood cell motion from repeated B-scans acquired in identical retinal locations. The instrument's software defined the superficial vascular complex within the nerve fiber and ganglion cell layers.

OCTA images with quality scores better than 3 out of 10 were selected for analysis. VD in the superficial vascular complex was determined using our previously published fractal analysis method. VD values ranged from zero to one, corresponding to the fraction of the imaged area occupied by vessels. OCTA images were also analyzed to determine vessel tortuosity along each vessel segment using our previously published method and calculate a mean vessel tortuosity index (VTI). The minimum value of VTI is zero, representing a straight line. VTI increases with higher variation in angles, number of critical points, and amplitude of curvature along the vessel centerline.

Statistical Analysis

Statistical Analysis System (SAS, version 9.4; SAS Institute Inc., Cary, NC) was used to analyze the data, and significance was accepted at P < 0.05. One eye of each subject was selected for the study. For subjects with data available in both eyes, one eye was selected at random by alternating the selection between the right and left eyes. Assumption of normality of data distribution was confirmed through the Shapiro-Wilk tests and graphical visualization of quantile-quantile plots. Demographics among groups (control, no DR, mild NPDR, moderate to severe NPDR, and PDR) were compared by ANOVA (age) and χ2 tests (eyes and sex). Race was compared using Fisher's exact test because the expected counts in more than 20% of the cells was less than 5. Comparison of all other continuous variables among groups was performed using ANOVA (HbA1c, VD, SO2A, SO2V, and OEF), or Kruskal-Wallis tests (MAP and VTI) for those that were not normally distributed. Homogeneity of variance of ANOVA analyses was confirmed by Levene's test. Post hoc analysis was performed using Tukey tests. Spearman rank correlation analysis was used to determine associations of VTI with SO2A, SO2V, and OEF. Partial correlations of VTI with SO2A, SO2V, and OEF were conducted while adjusting for HCT and HbA1c. The study had 90% power to detect a correlation coefficient of 0.32 or higher at the alpha level of 0.05 with a sample size of 95.

Results

Demographic characteristics of subjects stratified by group are shown in Table 1. There were no significant differences in sex (P = 0.64), eyes (P = 0.09), race (P = 0.52), or age (P = 0.43) among the control, no DR, mild NPDR, moderate to severe NPDR, and PDR groups.

Table 1.

Comparison of Demographics Among Groups

	Control	No DR	Mild NPDR	Moderate to	PDR
Variable	(N = 25)	(N = 28)	(N = 21)	Severe NPDR (N = 14)	(N = 7)	P Value
Sex
Male	11	11	9	7	5	0.64
Female	14	17	12	7	2
Eye
Right	11	12	16	6	5	0.09
Left	14	16	5	8	2
Race
Asian	8	3	3	2	0	0.52
African American	2	3	1	0	0
White	4	5	6	6	2
Hispanic	11	17	11	6	5
Age (years)
	50 ± 14	55 ± 15	58 ± 15	56 ± 17	52 ± 12	0.43

Data presented in control, no diabetic retinopathy (no DR), mild nonproliferative diabetic retinopathy (mild NPDR), moderate to severe NPDR, and proliferative diabetic retinopathy (PDR) groups.

Table 2 shows measurements of MAP, HbA1c, and HCT stratified by group. There was no significant difference in MAP among groups (P = 0.80). HbA1c was significantly different among groups (P < 0.0001). HbA1c levels were higher in the no DR, mild NPDR, moderate to severe NPDR, and PDR groups compared to the control group (P < 0.0001) and higher in the PDR group compared to the no DR group (P = 0.03). There was a statistically significant difference in HCT among groups (P = 0.001). HCT was lower in the mild NPDR group compared to the control group (P = 0.0001).

Table 2.

Comparison of Mean Arterial Pressure (MAP), Hemoglobin A1c (HbA1c), and Hematocrit (HCT) Among Groups

	Control	No DR	Mild NPDR	Moderate to
Metric	(N = 25)	(N = 28)	(N = 21)	Severe NPDR (N = 14)	PDR (N = 7)	P Value
MAP (mm Hg)	103 ± 14	102 ± 16	100 ± 11	103 ± 17	108 ± 15	0.80
HbA1c (%)	5.4 ± 0.47a	7.0 ± 1.1	7.7 ± 1.4	8.0 ± 1.4	8.4 ± 1.6	<0.0001
HCT (%)	45 ± 4	43 ± 3b	40 ± 5	43 ± 4	41 ± 5	0.001

Data presented in control, no diabetic retinopathy (no DR), mild nonproliferative diabetic retinopathy (mild NPDR), moderate to severe NPDR, and proliferative diabetic retinopathy (PDR) groups.

Bolded values indicate P < 0.05. Superscripts represent groups that had missing HbA1c or HCT data and thus a lower N value of subjects in that group,

N = 23.

N = 27.

Measurements of SO2A, SO2V, and OEF stratified by group are displayed in Table 3. SO2A was significantly different among groups (P = 0.002). SO2A was higher in the PDR group than the control group (P = 0.001). SO2v was significantly different among groups (P = 0.004). SO2v was higher in the moderate to severe NPDR (P = 0.01) and PDR (P = 0.04) groups compared to the control group. There was a marginal difference in OEF among groups (P = 0.07).

Table 3.

Comparison of Arterial and Venous Oxygen Saturation (SO2A and SO2V), Oxygen Extraction Fraction (OEF), Vessel Density (VD), and Vessel Tortuosity Index (VTI) Among Groups

	Control	No DR	Mild NPDR	Moderate to Severe	PDR
Metric	(N = 25)	(N = 28)	(N = 21)	NPDR (N = 14)	(N = 7)	P Value
SO2A (%)	89 ± 7	94 ± 8	94± 10	95 ± 8	104 ± 12	0.002
SO2V (%)	53 ± 8	58 ± 9	56 ± 10	64 ± 13	64 ± 9	0.004
OEF	0.41 ± 0.08	0.39 ± 0.08	0.40 ± 0.10	0.33 ± 0.12	0.38 ± 0.08	0.07
VD	0.49 ± 0.04	0.48 ± 0.03	0.47 ± 0.05	0.45 ± 0.04	0.43 ± 0.05	0.003
VTI	0.50 ± 0.15	0.53 ± 0.14	0.51 ± 0.13	0.62 ± 0.18	0.55 ± 0.18	0.22

Data presented in control, no diabetic retinopathy (no DR), mild nonproliferative diabetic retinopathy (mild NPDR), moderate to severe NPDR, and proliferative diabetic retinopathy (PDR) groups.

Bolded values indicate P < 0.05.

Measurements of VD and VTI stratified by group are shown in Table 3. VD differed significantly among groups (P = 0.003). VD was lower in the moderate to severe NPDR (P = 0.04) and PDR (P = 0.01) groups compared to the control group. VD was also lower in the PDR group than in the no DR group (P = 0.03). There was not a statistically significant difference in VTI among groups (P = 0.22).

There was a positive correlation between VTI and SO2V (rs = 0.32, P = 0.002, N = 95). The correlation between VTI and SO2A was not significant (rs = 0.05, P = 0.62, N = 95). VTI and OEF were negatively correlated (rs = −0.35, P = 0.001, N = 95). The significance of these relationships was unaffected by adjustment for HCT and HbA1c.

Discussion

Previous studies have shown increased retinal vascular oxygen saturation and tortuosity with progression of DR.– Here, alterations in retinal vascular morphologic and oxygenation metrics at stages of DR were demonstrated and associations between these metrics were reported for the first time. The results confirmed the hypothesis that higher perifoveal vessel tortuosity was correlated with higher retinal venous oxygen saturation and lower fraction of extracted oxygen.

In the current study, an increase in retinal SO2A was shown in PDR compared to less advanced stages of DR. This finding is in agreement with previous studies that showed elevation of SO2A with increasing severity of DR.,,,– Previous studies in DR have reported increased SO2V,, and decreased retinal arteriovenous SO2 difference and OEF., In the current study, OEF was not found to be significantly decreased, although a trend in OEF reduction was observed. Maintained OEF suggests the relative rate of oxygen metabolism to delivery was not altered, despite potential changes in either parameter. Previous studies have shown reduced OEF in NPDR., Heterogeneity of DR, as well as separate grouping of stages of NPDR and limited sample size in the current study, likely contributed to an undetected statistically significant reduction in OEF. Several factors may account for the observed elevation retinal vascular SO2 in DR. First, increased glycosylated hemoglobin has been shown to exhibit higher affinity to oxygen.– Therefore, oxygen may be more tightly bound to hemoglobin and diffuses less readily to the tissue, which can account for higher SO2 as oxygen remains in the lumen and attached to hemoglobin. Second, endothelial dysfunction and loss of regulatory cells due to DR,, may adversely affect the maintenance of homeostasis and the permeability to oxygen. Third, cellular dysfunction or loss due to hyperglycemia or vascular supply may result in a reduction in oxygen demand by the retinal tissue.

VD was reduced in moderate to severe NPDR and PDR, consistent with previous studies.,, The reduction in vessel density is due to capillary nonperfusion and the progressive loss of capillary cells. Previous studies have shown increased retinal vessel tortuosity during progression of DR.,, An increase in advanced glycation end products causes pericyte loss in DR, and impairs the ability of cells to respond to certain biochemical signals, such as nitric oxide, which is important in the functioning of smooth muscle cell relaxation during vasodilation. Studies have shown reduced flow-mediated vasodilation in DR, suggesting endothelial dysfunction., Overall, increased vessel tortuosity in DR is likely due to impairment of the retinal vessel endothelial support and function, as well as due to tissue hypoxia.

The results showed correlations of increased VTI with increased SO2V and decreased OEF, although a causal relation between these metrics is not indicated, requiring a different study design. Increased vessel tortuosity can occur in DR due to impairment of the retinal vessel endothelial support and function,, coupled with increased blood flow., Moreover, retinal vessels have been shown to be tortuous and elongated under DR, hypobaric hypoxia, and central retinal vein occlusion.,, On the other hand, decreased OEF (ratio of oxygen metabolism to oxygen delivery) is due to a decrease in oxygen metabolism and / or an increase in oxygen delivery. A decrease in OEF in DR can result with increased SO2V, which is suggestive of reduced oxygen metabolism coupled with increased blood flow, causing an increase in oxygen delivery.

The current study presented certain limitations. The relatively small sample size in some DR groups may have limited the statistical power to detect some differences among groups. Nevertheless, the sample size of all subjects had adequate power to establish correlations that were not influenced by DR stage. Additionally, the findings were based on data obtained in treated and untreated DR, hence heterogeneity of subjects may have had a variable effect on other aspects of retinal pathophysiology. Under normal physiology, the retinal and choroidal circulations predominately supply oxygen to the inner retinal neurons and photoreceptors, respectively. However, oxygen supply by the dual circulations may be altered under pathologic conditions. Using enhanced depth imaging OCT, several studies have investigated changes in the choroid due to DR. Overall, changes in the choroidal thickness were shown to be variable, whereas choroidal vascularity metrics were abnormal and decreased with worsening of DR. One limitation of the current study was that vascular oxygen saturation was only evaluated in the retinal circulation, thus the potential contribution of choroidal circulation to the findings was not considered. The effect of history of uncontrolled hypertension on the results was not considered, although at the time of imaging, subjects had similar BP. Because oxygenation and tortuosity measurements were obtained from different retinal regions, future studies are needed to assess the effect of retinal region and vessel size on these relationships. In addition, VTI measurements consisted of manual selection of vessel segments on a binary image based on an observer's subjective identification of locations of vessel branching. Therefore, in future studies, development of automatic segmentation for VTI measurements will improve the methodology. As technology advances, the current measurement methods may be improved, however, to our knowledge, the methods used in this study are among the best-established techniques in the field today.

In conclusion, the study showed a positive correlation of VTI with SO2V and a negative correlation with OEF, suggesting changes in these parameters develop under similar hyperglycemic situations and the pathological processes common in DR. The findings provide a better understanding of DR and may eventually lead to improved management of people with diabetes.