Investigative Ophthalmology & Visual Science Cover Image for Volume 66, Issue 4
April 2025
Volume 66, Issue 4
Open Access
Retina  |   April 2025
Structure-Function Associations Between Quantitative Contrast Sensitivity Function And Peripapillary Optical Coherence Tomography Angiography in Diabetic Retinopathy
Author Affiliations & Notes
  • Sierra K. Ha
    Harvard Retinal Imaging Lab, Massachusetts Eye and Ear, Harvard Medical School, Boston, Massachusetts, United States
  • Xinyi Ding
    Harvard Retinal Imaging Lab, Massachusetts Eye and Ear, Harvard Medical School, Boston, Massachusetts, United States
    Retina Service, Department of Ophthalmology, Massachusetts Eye and Ear, Boston, Massachusetts, United States
  • Francesco Romano
    Harvard Retinal Imaging Lab, Massachusetts Eye and Ear, Harvard Medical School, Boston, Massachusetts, United States
    Retina Service, Department of Ophthalmology, Massachusetts Eye and Ear, Boston, Massachusetts, United States
  • Katherine M. Overbey
    Harvard Retinal Imaging Lab, Massachusetts Eye and Ear, Harvard Medical School, Boston, Massachusetts, United States
  • Filippos Vingopoulos
    Harvard Retinal Imaging Lab, Massachusetts Eye and Ear, Harvard Medical School, Boston, Massachusetts, United States
  • Itika Garg
    Harvard Retinal Imaging Lab, Massachusetts Eye and Ear, Harvard Medical School, Boston, Massachusetts, United States
  • Cade F. Bennett
    Harvard Retinal Imaging Lab, Massachusetts Eye and Ear, Harvard Medical School, Boston, Massachusetts, United States
  • Isabella Stettler
    Harvard Retinal Imaging Lab, Massachusetts Eye and Ear, Harvard Medical School, Boston, Massachusetts, United States
  • Ioanna Ploumi
    Harvard Retinal Imaging Lab, Massachusetts Eye and Ear, Harvard Medical School, Boston, Massachusetts, United States
  • Grace Baldwin
    Harvard Retinal Imaging Lab, Massachusetts Eye and Ear, Harvard Medical School, Boston, Massachusetts, United States
  • Matthew J. Finn
    Harvard Retinal Imaging Lab, Massachusetts Eye and Ear, Harvard Medical School, Boston, Massachusetts, United States
  • Peyman Razavi
    Harvard Retinal Imaging Lab, Massachusetts Eye and Ear, Harvard Medical School, Boston, Massachusetts, United States
  • Demetrios G. Vavvas
    Retina Service, Department of Ophthalmology, Massachusetts Eye and Ear, Boston, Massachusetts, United States
  • Deeba Husain
    Retina Service, Department of Ophthalmology, Massachusetts Eye and Ear, Boston, Massachusetts, United States
  • Nimesh A. Patel
    Retina Service, Department of Ophthalmology, Massachusetts Eye and Ear, Boston, Massachusetts, United States
  • Leo A. Kim
    Retina Service, Department of Ophthalmology, Massachusetts Eye and Ear, Boston, Massachusetts, United States
  • John B. Miller
    Harvard Retinal Imaging Lab, Massachusetts Eye and Ear, Harvard Medical School, Boston, Massachusetts, United States
    Retina Service, Department of Ophthalmology, Massachusetts Eye and Ear, Boston, Massachusetts, United States
  • Correspondence: John B. Miller, Massachusetts Eye and Ear Infirmary Harvard Medical School, 243 Charles Street, Boston MA 02114, USA; [email protected]
  • Footnotes
     SKH and XD contributed equally to this work.
Investigative Ophthalmology & Visual Science April 2025, Vol.66, 69. doi:https://doi.org/10.1167/iovs.66.4.69
  • Views
  • PDF
  • Share
  • Tools
    • Alerts
      ×
      This feature is available to authenticated users only.
      Sign In or Create an Account ×
    • Get Citation

      Sierra K. Ha, Xinyi Ding, Francesco Romano, Katherine M. Overbey, Filippos Vingopoulos, Itika Garg, Cade F. Bennett, Isabella Stettler, Ioanna Ploumi, Grace Baldwin, Matthew J. Finn, Peyman Razavi, Demetrios G. Vavvas, Deeba Husain, Nimesh A. Patel, Leo A. Kim, John B. Miller; Structure-Function Associations Between Quantitative Contrast Sensitivity Function And Peripapillary Optical Coherence Tomography Angiography in Diabetic Retinopathy. Invest. Ophthalmol. Vis. Sci. 2025;66(4):69. https://doi.org/10.1167/iovs.66.4.69.

      Download citation file:


      © ARVO (1962-2015); The Authors (2016-present)

      ×
  • Supplements
Abstract

Purpose: To assess changes in radial peripapillary capillary (RPC) microvasculature and their impact on visual function, measured by visual acuity (VA) and contrast sensitivity, in diabetic retinopathy (DR).

Methods: This was a cross-sectional study in 96 eyes of 67 patients, including controls, diabetes without DR (DMnoDR), nonproliferative DR (NPDR), and proliferative DR (PDR) groups. Participants underwent same-day quantitative contrast sensitivity function (qCSF) and 6 × 6 mm OCT angiography (OCTA) centered on the optic disc. The Peripapillary Nerve Fiber Layer Microvasculature Density algorithm (ARI Network) was used to calculate capillary perfusion density (total area of perfused microvasculature per unit area), capillary flux index (CFI, total weighted area of perfused microvasculature per unit area), and retinal nerve fiber layer (RNFL) thickness surrounding the optic disc. Mixed-effects multivariable regression models, controlling for age, hypertension, and lens status, evaluated associations between RPC OCTA metrics, DR severity, VA, and qCSF.

Results: Significant RPC microvascular changes were observed across DR stages. Capillary perfusion density decreased with DR severity and even before retinopathy onset in DMnoDR versus controls (βavg = −0.42, P = 0.021). PDR compared to NPDR showed a significant decrease in CFI (β = −1.02 to −0.92, P < 0.01) and in RNFL (βavg = −0.71, P = 0.033). CFI had significant associations with qCSF at various spatial frequencies (β = 0.20 to 0.34, P = 0.002 to 0.042), but not with VA.

Conclusions: Radial peripapillary capillary perfusion density worsens with onset of diabetes and increasing severity of DR while capillary flux index is more significantly affected later in disease. Structure-function associations suggest that DR-induced peripapillary microvascular changes are more strongly associated with contrast sensitivity changes than with visual acuity.

Diabetic retinopathy (DR) is a leading cause of preventable blindness among working-aged adults worldwide, with approximately one third of people with diabetes having concurrent DR.1 The pathophysiology of DR involves both retinal microvascular and neural alterations, with many studies reporting on these alterations within the macula.24 However, emerging evidence suggests that the peripapillary area may also be affected early in the course of DR, with microvascular alterations potentially serving as preclinical indicators of diabetic microvascular disease.5 
The radial peripapillary capillary (RPC) plexus, a network of long, straight vessels with minimal anastomoses, supplies the retinal nerve fiber layer (RNFL) and plays a crucial role in maintaining its structural integrity and function.6,7 Optical coherence tomography angiography (OCTA) provides a noninvasive method of measuring RPC microvasculature through two key metrics: peripapillary capillary perfusion density and capillary flux index (CFI). Capillary perfusion density is defined as the area of perfused vasculature per unit area and is analogous to the commonly used vessel density (VD), offering a structural measure of perfusion. In contrast, CFI is defined as the area of perfused vasculature per unit area weighted by normalized pixel flow intensity, providing a functional assessment of perfusion and flow dynamics.8 
Previous studies suggest that damage to the RPC microvasculature is associated with increasing DR severity and correlates with changes in thickness of the RNFL.5,810 However, findings on the relationship between RNFL thickness and DR severity remain inconsistent, with some studies identifying significant associations5,8 and others reporting no clear relationship.9,10 These discrepancies highlight the need for further investigation into the structural and functional changes in DR. 
The resulting impact of microvascular and RNFL alterations on visual function remains largely unexplored in DR. Traditionally, visual acuity (VA) has been the primary measure of visual function in clinical trials.11 However, growing evidence supports the use of alternative tests for visual function, such as contrast sensitivity (CS), given it is affected earlier in the course of neurodegenerative disorders and may be more reflective of vision-related quality of life and subtle changes in visual function than VA.1215 Although limitations in testing modalities have previously prevented the integration of CS into routine clinical practice, the advent of the quantitative contrast sensitivity function (qCSF) computational approach has facilitated the assessment of CS in various retinal conditions.1618 Our group has previously demonstrated that DR-induced microvascular changes in the superficial capillary plexus and deep capillary plexus measured by expanded field swept-source OCTA centered on fovea are more strongly associated with changes in CS than VA. Additionally, CS is affected earlier than VA in the course of DR and performs better in discriminating between controls, diabetes mellitus without DR, and across DR stages.19 
No studies to date have explored how changes in RPC vasculature correlate with various visual function metrics, including CS and VA, in DR. Clarifying the relationship between RPC alterations and qCSF may help elucidate how visual function in DR is affected in regions beyond the macula. Understanding these associations is crucial for the early diagnosis and effective monitoring of DR progression. Therefore the current study aims to investigate the association between alterations in RPC microvasculature and RNFL thickness in the peripapillary area with functional changes in visual acuity and qCSF in DR. 
Methods
Study Design
This is a cross-sectional observational, single-center study including patients seen in the outpatient retina clinics at Massachusetts Eye and Ear between June 2019 and January 2024. The study was approved by the Institutional Review Board of Massachusetts General Brigham (2019P001863), and informed consent was obtained from all subjects. All procedures were conducted in accordance with the tenets set forth in the Declaration of Helsinki and Accountability Act regulations. 
Participants
Participants were included based on the following criteria: (1) age of 18 years or older; (2) diagnosis of Type 1 or Type 2 diabetes mellitus (DM) using the American Diabetes Association criteria20; (3) current diagnosis of DM without DR (DMnoDR) or DR without diabetic macular edema (DME) by a senior retina specialist (LAK, DH, DGV, NAP, or JBM) on the day of study visit; (4) visual acuity of at least 20/200 as measured by a standard Snellen chart; (5) same-day quantitative contrast sensitivity function (qCSF) testing and imaging with 6 × 6 mm OCTA centered on the optic disc. Exclusion criteria were as follows: (1) other ocular comorbidities (e.g., glaucoma or age-related macular degeneration); (2) significant media opacities including corneal scars and lens opacity greater than grade > 3+ following the Lens Opacities Classification System (LOCS) III21; (3) significant artifacts or signal strength less than 7 on OCTA imaging; (4) neovascularization of the disc, due to potential confounding of OCTA metrics. A control group of patients with same-day qCSF and peripapillary OCTA was derived from our previously published normative database, which includes healthy eyes with no ocular pathology other than non-visually significant cataracts.22 
All subjects underwent a comprehensive ophthalmologic evaluation, including slit lamp examination, measurement of VA using a standard Snellen chart, qCSF testing (Manifold Platform; Adaptive Sensory Technology, San Diego, CA), measurement of intraocular pressure, and dilated fundus examination and photography (Optos California; Optos plc., Dunfermline, UK), and 6 × 6 mm OCTA scans of the optic disc (PLEX Elite 9000; Carl Zeiss Meditec, Dublin, CA, USA). Lens status grading followed the LOCS III and was simplified for the purposes of the multivariable regression analysis as follows: clear lens was graded as “clear”; NO1NC1 graded as 1+ NS (nuclear sclerosis); NO2NC1, NO1NC2 or NO2NC2 was graded as 2+ NS; and NO3NC1-3 and NO1-2NC3 graded as 3+ NS. Staging of DR into nonproliferative diabetic retinopathy (NPDR) or proliferative diabetic retinopathy (PDR) was established by a retina specialist based on the Early Treatment Diabetic Retinopathy Study grading system.23 
Contrast Sensitivity Testing
Contrast sensitivity was evaluated using the qCSF method on the AST platform (Adaptive Sensory Technology, San Diego, CA, USA) as previously described by our group.24 In brief, the qCSF method evaluates the ability to distinguish different contrast ratios across a range of spatial frequencies. Patients are initially presented with three letters at a single spatial frequency and decreasing levels of contrast. A technician records the patient's responses (correct, incorrect, or no answer) for each letter. Using a Bayesian active learning algorithm, the next combination of contrast and frequency is chosen to maximize information gain and reduce testing time, with a total of 25 triplets of optotypes tested. The time required for testing is 2-5 minutes per eye. Various qCSF values derived include: (1) area under the logarithm of contrast sensitivity function (AULCSF) represents CS estimates integrated across spatial frequencies ranging from 1.5 to 18 cycles per degree (cpd); (2) contrast acuity (CA) corresponds with the intersection of CSF curve with the x-axis and is the smallest optotype seen at full (100%) contrast; (3) contrast sensitivity thresholds at 1, 1.5, 3, 6, 12, and 18 cpd provide the lowest contrast visible at each spatial frequency in logCS units. The qCSF method has demonstrated good test-retest reliability as compared to other contrast sensitivity tests that evaluate thresholds at various spatial frequencies.16,25 
OCTA Imaging and Analysis
Trained research technicians acquired 6 × 6 mm OCTA images centered on the optic disc using a 100-kHz SS-OCTA instrument (Plex Elite 9000, Carl Zeiss Meditec Inc., Dublin, CA). Images were manually assessed for quality, and those with significant artifacts or signal strength less than 7 were excluded. The Peripapillary Nerve Fiber Layer Microvasculature Density algorithm v0.9, a prototype research algorithm available on the ARI Network (Zeiss Portal v5.4–1206) was utilized to analyze imaging. This algorithm performs segmentation of the inner limiting membrane and the RNFL to create an en face image of the RPC vasculature (Fig. 1). Because this algorithm is a prototype, the automated segmentation and all output images and values were reviewed manually for accuracy of the analysis, and manual segmentation was performed in cases with segmentation errors. The algorithm outputs include capillary perfusion density, capillary flux index (CFI), and retina nerve fiber layer (RNFL) thickness in four quadrants (temporal, superior, nasal, inferior) of an annulus with an outer diameter of 6 mm and inner diameter of 2 mm centered on the optic disc, as well as the average of the four quadrants. Capillary perfusion density represents the total area of perfused microvasculature per unit area in a region of measurement, whereas CFI represents the total weighted area of perfused microvasculature per unit area in a region of measurement, with the weight being the normalized flow intensity corresponding to each pixel. RNFL thickness is provided in micrometers for the same analyzed regions of interest. 
Figure 1.
 
Representative images of RPC vasculature enface, structural enface, RPC perfusion heat map, and RNFL thickness map by DR stage.
Figure 1.
 
Representative images of RPC vasculature enface, structural enface, RPC perfusion heat map, and RNFL thickness map by DR stage.
Statistical Analysis
Statistical analysis was performed using RStudio Version 2023.12.1+402 and Stata/BE 18.0. Demographics and clinical characteristics were described using traditional descriptive methods. Non-normally distributed data were reported as medians with interquartile ranges. Snellen visual acuity was converted to the logarithm of the minimum angle of resolution (logMAR) equivalent. Initial variables for multivariate modeling were selected based on their potential to confound outcomes related to DR microvascular progression and visual function.2628 These included age, gender, lens status, duration of diabetes, coexisting hypertension, and type of diabetes. Variables were initially chosen based on their unadjusted statistical significance in univariate analyses and further refined using backward stepwise elimination to develop the final multivariable regression models. Participants were divided into groups including controls, DMnoDR, NPDR (including mild, moderate, and severe NPDR), and PDR. Mixed-effects multivariable regression models controlling for age and hypertension were used to compare RPC microvascular metrics between groups. Mixed-effects multivariable regression models controlling for age and lens status were used to evaluate associations between OCTA metrics and both VA and qCSF metrics. In all models, a random effect for each participant was included to account for the potential correlation between two eyes from the same individual. Associations were reported as standardized beta coefficients by refitting standardized data. Statistical significance was defined as a P value <0.05. 
Results
Demographics and Clinical Characteristics
The study included 96 eyes from 67 participants, comprising 23 control eyes, 27 DMnoDR eyes, 27 NPDR eyes, and 19 PDR eyes. The demographics and clinical characteristics are summarized in Table 1. The median age was 59 (47–67) years, and 40.3% of participants were female. Among the eyes with DM, 61.6% had type 2 diabetes. 
Table 1.
 
Clinical and Demographic Features of the Studied Cohort
Table 1.
 
Clinical and Demographic Features of the Studied Cohort
Visual Acuity, Contrast Sensitivity, and OCTA Metrics
VA and qCSF outcomes are summarized in Table 1. The median LogMAR VA was 0.04 (0.00–0.16). The median AULCSF was 1.10 (0.91–1.32). Figure 1 illustrates the output of the Peripapillary Nerve Fiber Layer Microvasculature Density algorithm v0.9 across different stages of DR. The RPC OCTA microvascular metrics and RNFL thickness results are summarized across different stages of DR in Table 2. Median CFI averaged across the four quadrants was 0.40 (0.39–0.41) in controls, 0.41 (0.39–0.42) in DMnoDR, 0.40 (0.37–0.41) in NPDR, and 0.38 (0.36–0.39) in the PDR group (Fig. 2). Median capillary perfusion density averaged across the four quadrants was 0.57 (0.56–0.58) in controls, 0.57 (0.56–0.57) in DMnoDR, 0.56 (0.55–0.57) in NPDR, and 0.55 (0.54–0.56) in the PDR group (Fig. 2). Median RNFL thickness averaged across the four quadrants was 85.25 (79.83–93.19) in controls, 87.19 (82.37–91.25) in DMnoDR, 87.20 (78.88–93.13) in NPDR, and 74.32 (71.39–81.50) in PDR. 
Table 2.
 
Radial Peripapillary Capillary Optical Coherence Tomography Angiography Metrics Among All Eyes
Table 2.
 
Radial Peripapillary Capillary Optical Coherence Tomography Angiography Metrics Among All Eyes
Figure 2.
 
Peripapillary microvascular metrics across stages of DR. (A) Capillary perfusion density, mean and by quadrant, across DR stages. (B) CFI, mean and by quadrant, across DR stages. *P < 0.05.
Figure 2.
 
Peripapillary microvascular metrics across stages of DR. (A) Capillary perfusion density, mean and by quadrant, across DR stages. (B) CFI, mean and by quadrant, across DR stages. *P < 0.05.
Peripapillary Microvasculature Between Group Comparisons
Mixed-effects multivariable regression models controlling for age and hypertension were used to compare RPC microvascular metrics across various stages of DR (Table 3). Eyes with PDR showed a significant decrease in CFI in all quadrants compared to NPDR (β = −1.02 to −0.92, P = 0.003 to <0.001) (Table 3). There were no significant differences in CFI between DMnoDR, NPDR, and control eyes. Eyes with increasing severity of DR were associated with decreased capillary perfusion density in multiple quadrants. Specifically, eyes with DMnoDR showed decreased average, temporal, and inferior capillary perfusion density compared to controls (βavg = −0.42, P = 0.021; βtemp = −0.65, P = 0.001; βinf =-0.63, P = 0.007) (Table 3). Differences in capillary perfusion density were also found between control eyes and NPDR, DMnoDR and NPDR, and NPDR and PDR eyes. The largest effect size was detected for the average capillary perfusion density between eyes with NPDR and PDR (β = −0.89, P = 0.009) (Table 3). Eyes with DMnoDR were associated with an increased RNFL in the superior quadrant compared to control eyes (β = 0.61, P = 0.015). Eyes with NPDR did not show a significant difference in RNFL compared to DMnoDR. However, eyes with PDR compared to NPDR showed a decrease in average, nasal, and inferior RNFL (βavg = −0.71, P = 0.033; βnas = −0.58, P = 0.047; βinf = −0.87, P = 0.006) (Table 3). 
Table 3.
 
Mixed Effects Multivariable Regression Analysis, Controlling for Age and Hypertension, Comparing Radial Peripapillary Capillary OCTA Metrics Across Different Stages of Diabetic Retinopathy
Table 3.
 
Mixed Effects Multivariable Regression Analysis, Controlling for Age and Hypertension, Comparing Radial Peripapillary Capillary OCTA Metrics Across Different Stages of Diabetic Retinopathy
Structure-Function Analysis
Mixed-effects multivariable regression analyses controlling for age and lens status were used to evaluate the associations between OCTA metrics and visual function (Table 4). CFI in only the temporal quadrant showed a significant association with decreased logMAR VA (β = −0.22, P = 0.049), suggesting that higher CFI in this quadrant is correlated with better visual acuity. Conversely, CFI in all quadrants showed significant associations with multiple CS metrics, including CA and CS at thresholds of 6–18 cpd. Additionally, CFI on average and in all quadrants with the exception of the nasal quadrant was associated with AULCSF. For temporal CFI, standardized β coefficients for AULCSF, CA, and CS at 6, 12, and 18 cpd were all greater than or equal to effect size for VA (0.22, 0.25, 0.22, 0.25, 0.33 vs. −0.22, respectively) (Table 4). 
Table 4.
 
Mixed Effects Multivariable Regression Analysis Controlling for Age and Lens Status Evaluating Associations Between Radial Peripapillary Capillary OCTA Metrics and Visual Acuity (VA) and Contrast Sensitivity (CS) Outcome Measures
Table 4.
 
Mixed Effects Multivariable Regression Analysis Controlling for Age and Lens Status Evaluating Associations Between Radial Peripapillary Capillary OCTA Metrics and Visual Acuity (VA) and Contrast Sensitivity (CS) Outcome Measures
Capillary perfusion density on average or in each quadrant individually was not significantly associated with VA and CS metrics (all p-values > 0.05) (Table 4). Additionally, average and nasal RNFL thickness were not associated with VA and CS metrics. Temporal RNFL thickness was significantly associated with decreased CS at 6 cpd (β = −0.22, P = 0.022). Similar negative associations were observed between superior RNFL thickness and CS at 1 cpd (β = −0.21, P = 0.028) and inferior RNFL thickness and CS at 3 cpd (β = −0.24, P = 0.001) (Table 4). 
Discussion
Here, we present the first study investigating the structure-function associations between peripapillary microvascular changes on OCTA and contrast sensitivity in diabetic patients without DR and across various stages of DR. In our cohort, which underwent same-day qCSF-measured CS and OCTA centered on the optic disc, we found significant reductions in peripapillary microvascular metrics, as measured by CFI and capillary perfusion density, with increasing DR severity. Notably, peripapillary microvascular alterations due to DR correlated more with CS than with VA, underscoring the importance of peripapillary microvascular health in regulating visual function, particularly CS, in patients with DR. 
The RPC plexus in the peripapillary region is crucial for maintaining the health of the RNFL by supporting the neural metabolic demands around the optic nerve head. Histologic studies reveal that the RPC plexus consist of long, straight vessels with few anastomoses.7 Non-myelinated axons of retinal ganglion cells in the RNFL have high energy demands, making them highly vulnerable to ischemic insults.29 Consequently, damage to the RPC plexus can lead to significant visual impairment. Early changes in peripapillary microvasculature may precede clinically detectable DR, highlighting the importance of advanced imaging and functional biomarkers.5 
Mixed-effects multivariable regression analysis demonstrated progressive reductions in RPC capillary perfusion density after onset of DM and across DR stages. Capillary perfusion density in the temporal quadrant and average of all quadrants was the most affected. Meanwhile, eyes with PDR versus NPDR showed significant reduction in CFI in all quadrants and on average, suggesting that advanced DR stages are associated with substantial RPC functional microvascular compromise. These findings align with other studies that have also demonstrated early changes occurring in the peripapillary microvasculature of patients with DM without DR compared to controls and progressive reduction with increasing DR severity.5,3032 Prior data from Mitamura et al.8 described changes in RPC CFI occurring at the onset of diabetes and changes in RPC capillary perfusion density occurring at the onset of DR. Although these findings appear to differ from our study, where capillary perfusion density is affected before CFI, several factors may account for this discrepancy. Differences in the OCTA devices (PLEX Elite 9000 in our study vs. Cirrus HD-OCT 5000 with AngioPlex OCT Angiography in Mitamura et al.8) and differences in the size of the examined retinal area (6 × 6 mm in our study versus 4.5 × 4.5 mm in Mitamura et al.8) may have contributed to the contrasting findings. Aditionally, controlling for concurrent hypertension in our multivariable regression models, which was not done in the study by Mitamura et al.,8 may also help explain the differences between the studies.33,34 Capillary perfusion density, which measures total area of perfused vasculature per unit area, serves as an indicator of microvascular structure. In contrast, CFI, which measures capillary perfusion density weighted by the brightness of flow signal, provides a more qualitative assessment of microvascular function and perfusion effectiveness.8,35 The initial changes in capillary perfusion density may be reflective of structural damage to vascular cells, including endothelial cells and pericytes. This ultimately leads to reduced perfusion effectiveness in the microvasculature as measured by CFI reduction. 
Few changes in peripapillary RNFL were present across stages of DR. Our findings show significant thinning of the RNFL in eyes with PDR compared to NPDR, particularly in the nasal and inferior quadrants, as well as the average of all quadrants. Potential mechanisms of RNFL thinning in eyes with PDR may include primary neural degeneration, secondary effects after panretinal photocoagulation,36 or thinning secondary to RPC microvascular changes, as indicated by decreases in capillary perfusion density and CFI in the corresponding quadrants. Changes in RPC hemodynamics may result in decreased oxygen and nutrient delivery to retinal ganglion cells, leading to axonal loss and RNFL thinning.37 Interestingly, the superior quadrant RNFL was thicker in DMnoDR eyes compared to controls, contrary to previous studies that have demonstrated thinning of the RNFL in patients with early stages of diabetes without retinopathy.38,39 The current study did not account for other variables such as ethnicity, which has been associated with changes in RNFL thickness and may explain differences in our findings.40 Moreover, studies on RNFL thickness in the progression of DR have demonstrated mixed results, with findings ranging from thinning, to no change, to increased RNFL thickness with increased severity of DR.5,9,31,41 Our findings suggest dynamic changes in RNFL thickness in the context of both DM and DR. Although not designed as a longitudinal study, these results indicate a more complex progression of RNFL changes in DR, with an initial increase in the early phases followed by subsequent thinning. 
Our study also assessed the relationships between DM-induced microvascular changes using peripapillary OCTA and functional changes in CS and VA. After controlling for age and lens status, mixed-effects regression models revealed a stronger effect of peripapillary OCTA metrics on CS, as evidenced by the larger standardized beta coefficients observed for CFI to predict various qCSF metrics compared to VA. Although CFI was associated with VA only in the temporal quadrant, CFI across all quadrants demonstrated robust associations with most qCSF metrics. Significant metrics included AULCSF, CA, and CS values at 6–18 cpd, with greater standardized beta coefficients observed for CS at higher spatial frequencies (6–18 cpd). Interestingly, capillary perfusion density was not associated with any measurements of visual function. This finding may be due to the fact that capillary perfusion density reflects the overall area of perfused microvasculature but does not capture specific functional impairments as sensitively as other metrics. We also found that RNFL thickness was associated with fewer qCSF metrics and showed smaller effect sizes. Specifically, RNFL thickness in the superior, inferior, and temporal quadrants were inversely associated with CS at thresholds of 1, 3, and 6 cpd, respectively. The clinical implications of these findings are significant. Early detection of microvascular changes is crucial for managing DR and improving outcomes. The sensitivity of capillary perfusion density and CFI to microvascular and functional changes across stages of DR highlight their potential as biomarkers for DR. These OCTA-derived metrics could be combined with functional testing such as qCSF to identify subtle abnormalities at earlier disease stages. This could guide earlier intervention strategies, such as tighter metabolic control or closer monitoring, to delay DR progression and visual decline. In addition, these findings further support the incoporation of qCSF as an endpoint in clinical trials, as it provides a more significant correlation with microvascular markers than visual acuity. 
The clinical implications of these findings are significant. Early detection of microvascular changes remains a cornerstone in managing DR, and our results suggest that integrating peripapillary OCTA metrics, such as CFI and capillary perfusion density, with functional testing like qCSF may enhance the ability to identify early structural and functional abnormalities in DR. 
To our knowledge, no other studies have evaluated associations between peripapillary OCTA and CS in DR. However, in glaucoma, RPC density has been associated with multiple qCSF metrics, including AULCSF, CA, and CS at 1–18 cpd.42 The lack of association of CS with capillary perfusion density, despite robust associations with CFI, could be attributed to the fact that although capillary perfusion density is an indicator of structural microvascular change, it is not as sensitive as CFI in measuring functional changes in perfusion effectiveness. Reduced capillary perfusion density may indicate early structural damage to the microvasculature, whereas reductions in CFI, which accounts for the brightness of the flow signal, reflect more advanced and functionally significant impairment in blood flow and perfusion.8 Although our findings of some associations between RNFL thinning and qCSF metrics are intriguing, some studies on other types of neurodegeneration have shown that AULCSF remains unchanged, despite progressive loss of RNFL, until neuronal loss surpasses a certain threshold.43 If such a threshold effect exists in the context of DR, our regression analyses may not capture the full impact. Additionally, increased RNFL thickness may reflect inflammation, glial cell activation, axonal swelling, or other secondary neuronal changes in response to underlying pathology that already compromises visual function.44,45 The findings of our study differ from those of Dong et al.,46 who reported associations between CS and greater RNFL and GCC thickness but found no correlations with angiographic features in a population of older adults. Several factors may explain these discrepancies, including differences in the methods used to measure CS (MARS chart in Dong et al.46 vs. qCSF in our study) and variations in patient demographics (the study by Dong et al.46 had a patient population with 46% of the participants identifying as Black). Additionally, the study described by Dong et al.46 only assessed angiographic features surrounding the macula. Whereas Dong et al.46 focused on older adults without significant microvascular compromise, our cohort included diabetic patients across various DR stages, where there may be disruptions in the normal structural and functional neurovascular coupling within the retina. Further understanding of the influence of RNFL on CS in DR may require studies examining these outcomes and their correlation with microvascular metrics, given that disruptions in the neurovascular coupling, or the interactions between the neurosensory retina and its blood vessels, are known to play a key role in the development of DR.37,41,47 
Limitations of the current study include its relatively small sample size and single-center setting, which may have masked subtle changes between different stages of DR. The cross-sectional design precludes causal inferences; therefore longitudinal studies are needed to elucidate the temporal relationships between microvascular changes, RNFL thinning, and functional decline. Additionally, although the study controlled for age and hypertension, other potential confounders, such as macular status and cerebrovascular disease, which may affect peripapillary microvascular status, were not considered because of limitations of our study database.4850 Further studies incorporating other imaging and molecular biology evidence may also provide a more comprehensive understanding of the mechanisms underlying microvascular structural and functional changes. 
The absence of axial length measurements is a critical limitation of this study, because individual differences in axial length can alter the scaling of OCTA images, resulting in discrepancies in the actual physical dimensions of the scanned retinal area. This scaling variability may introduce systematic inaccuracies in the quantitative OCTA measurements, because the nominal scan size may not accurately represent the actual retinal area being analyzed.5153 Future studies should incorporate direct axial length measurements or surrogate metrics, such as spherical equivalent or corneal curvature, to correct OCTA and RNFL data, ensuring the accuracy and comparability of these measurements.54 Alternatively, axial length could be estimated from macular optical coherence tomography images using a deep learning-based approach,55 or a Monte Carlo simulation approach could be used to adjust for image scaling given the known distribution of axial length.56 
In summary, the novel structure-function relationships examined in this study by combining peripapillary OCTA with a comprehensive measurement of CS provide insights into the underlying mechanisms of how peripapillary retinal dysfunction affects visual function. This study reveals that, in addition to macular alterations, the microvascular structure of the peripapillary region may also significantly impact CS in DR. Additionally, qCSF shows greater sensitivity to early peripapillary structural changes caused by DR compared to VA, suggesting that incorporating qCSF testing into routine clinical practice could facilitate a more comprehensive assessment of visual function in the clinical management of DR. 
Acknowledgments
Disclosure: S.K. Ha, None; X. Ding, None; F. Romano, None; K.M. Overbey, None; F. Vingopoulos, None; I. Garg, None; C.F. Bennett, None; I. Stettler, None; I. Ploumi, None; G. Baldwin, None; M.J. Finn, None; P. Razavi, None; D.G. Vavvas, Olix Pharma (C), Valitor (C), TwentyTwenty (C), Sumitomo/Sunovion (C), Cambridge Polymer Group (C), National Eye Institute (F), Research to Prevent Blindness (F), Loefflers Family Foundation (F), Yeatts Family Foundation (F), Alcon Research Institute (F); D. Husain, Allergan (C), Genentech (C), Omeicos Therapeutics (C), National Eye Institute (F), Lions Vision Gift (F), Commonwealth Grant (F), Lions International (F), Syneos LLC (F), the Macular Society (F); N.A. Patel, Alimera Sciences (C), Alcon (C), Allergan (C), Genentech (C), Eyepoint (C), Lifesciences Guidepoint (C), GLG (C); L.A. Kim, National Eye Institute (F), CureVac AG (F), Pykus Therapeutics (F); J.B. Miller, Alcon (C), Allergan (C), Carl Zeiss (C), Sunovion (C), Topcon (C), Genentech (C) 
References
Cheung N, Mitchell P, Wong TY. Diabetic retinopathy. Lancet. 2010; 376(9735): 124–136. [CrossRef] [PubMed]
Curtis TM, Gardiner TA, Stitt AW. Microvascular lesions of diabetic retinopathy: clues towards understanding pathogenesis? Eye. 2009; 23: 1496–1508. [CrossRef] [PubMed]
Baldwin G, Vingopoulos F, Garg I, et al. Structure–function associations between contrast sensitivity and widefield swept–source optical coherence tomography angiography in diabetic macular edema. Graefes Arch Clin Exp Ophthalmol. 2023; 261: 3113–3124. [CrossRef] [PubMed]
Kim K, Kim ES, Yu SY. Optical coherence tomography angiography analysis of foveal microvascular changes and inner retinal layer thinning in patients with diabetes. Br J Ophthalmol. 2018; 102: 1226–1231. [CrossRef] [PubMed]
Vujosevic S, Muraca A, Gatti V, et al. Peripapillary microvascular and neural changes in diabetes mellitus: an OCT-angiography study. Invest Ophthalmol Vis Sci. 2018; 59: 5074–5081. [CrossRef] [PubMed]
Liu G, Wang Y, Gao P. Distributions of radial peripapillary capillary density and correlations with retinal nerve fiber layer thickness in normal subjects. Med Sci Monit. 2021; 27: e933601-1–e933601-10.
Campbell JP, Zhang M, Hwang TS, et al. Detailed vascular anatomy of the human retina by projection-resolved optical coherence tomography angiography. Sci Rep. 2017; 7(1): 42201. [CrossRef] [PubMed]
Mitamura M, Endo H, Kase S, et al. Peripapillary circulatory dysfunction precedes structural loss in treatment-naive diabetic retinopathy. Graefes Arch Clin Exp Ophthalmol. 2023; 261: 85–95. [CrossRef] [PubMed]
Liu L, Wang Y, Liu HX, Gao J. Peripapillary region perfusion and retinal nerve fiber layer thickness abnormalities in diabetic retinopathy assessed by OCT angiography. Transl Vis Sci Technol. 2019; 8(4): 14. [CrossRef] [PubMed]
Huang J, Zheng B, Lu Y, Gu X, Dai H, Chen T. Quantification of microvascular density of the optic nerve head in diabetic retinopathy using optical coherence tomographic angiography. J Ophthalmol. 2020; 2020: e5014035.
Levine SR, Myers MG, Barunas R, et al. Report From the 2022 Mary Tyler Moore Vision Initiative Diabetic Retinal Disease Clinical Endpoints Workshop. Transl Vis Sci Technol. 2023; 12(11): 33. [CrossRef] [PubMed]
Pondorfer SG, Terheyden JH, Heinemann M, Wintergerst MWM, Holz FG, Finger RP. Association of vision-related quality of life with visual function in age-related macular degeneration. Sci Rep. 2019; 9(1): 15326. [CrossRef] [PubMed]
Sokol S, Moskowitz A, Skarf B, Evans R, Molitch M, Senior B. Contrast sensitivity in diabetics with and without background retinopathy. Arch Ophthalmol. 1985; 103: 51–54. [CrossRef] [PubMed]
Wai KM, Vingopoulos F, Garg I, et al. Contrast sensitivity function in patients with macular disease and good visual acuity. Br J Ophthalmol. 2022; 106: 839–844. [CrossRef] [PubMed]
Vingopoulos F, Bannerman A, Zhou P, et al. Towards the validation of quantitative contrast sensitivity as a clinical endpoint: correlations with vision-related quality of life in bilateral AMD. Br J Ophthalmol. 2024; 108: 846–851. [CrossRef] [PubMed]
Finn M, Vingopoulos F, Zhao Y, et al. Test-retest repeatability and agreement of the quantitative contrast sensitivity function test: towards the validation of a new clinical endpoint. Graefes Arch Clin Exp Ophthalmol. 2024; 262: 813–822. [CrossRef] [PubMed]
Thomas M, Silverman RF, Vingopoulos F, et al. Active learning of contrast sensitivity to assess visual function in macula-off retinal detachment. J Vitreoretin Dis. 2020; 5: 313–320. [CrossRef] [PubMed]
Joltikov KA, Sesi CA, de Castro VM, et al. Disorganization of retinal inner layers (DRIL) and Neuroretinal Dysfunction in Early Diabetic Retinopathy. Invest Ophthalmol Vis Sci. 2018; 59: 5481–5486. [CrossRef] [PubMed]
Vingopoulos F, Baldwin G, Katz R, et al. Associations of quantitative contrast sensitivity with vascular metrics on widefield swept-source OCT angiography across stages of diabetic retinopathy. Br J Ophthalmol. 2024; 108: 1708–1715. [CrossRef] [PubMed]
American Diabetes Association. 2. Classification and diagnosis of diabetes: standards of medical care in diabetes–2021. Diabetes Care. 2021; 44(Suppl 1): S15–S33. [PubMed]
Chylack LT, Jr, Wolfe JK, Singer DM, et al. The Lens Opacities Classification System III. Arch Ophthalmol. 1993; 111: 831–836. [CrossRef] [PubMed]
Choi H, Vingopoulos F, Razavi P, et al. Quantitative contrast sensitivity function and the effect of aging in healthy adult eyes: a normative database. Ophthalmic Surg Lasers Imaging Retina. 2024; 55: 212–219. [CrossRef] [PubMed]
Early Treatment Diabetic Retinopathy Study Research Group. Grading diabetic retinopathy from stereoscopic color fundus photographs–an extension of the modified Airlie House classification. ETDRS report number 10. Ophthalmology. 1991; 98(5 Suppl): 786–806. [PubMed]
Vingopoulos F, Wai KM, Katz R, Vavvas DG, Kim LA, Miller JB. Measuring the contrast sensitivity function in non-neovascular and neovascular age-related macular degeneration: the Quantitative Contrast Sensitivity Function Test. J Clin Med. 2021; 10: 2768. [CrossRef] [PubMed]
Richman J, Spaeth GL, Wirostko B. Contrast sensitivity basics and a critique of currently available tests. J Cataract Refract Surg. 2013; 39: 1100–1106. [CrossRef] [PubMed]
Wei Y, Jiang H, Shi Y, et al. Age-related alterations in the retinal microvasculature, microcirculation, and microstructure. Invest Ophthalmol Vis Sci. 2017; 58: 3804–3817. [CrossRef] [PubMed]
Voigt M, Schmidt S, Lehmann T, et al. Prevalence and progression rate of diabetic retinopathy in type 2 diabetes patients in correlation with the duration of diabetes. Exp Clin Endocrinol Diabetes. 2018; 126: 570–576. [PubMed]
Estacio RO, Jeffers BW, Gifford N, Schrier RW. Effect of blood pressure control on diabetic microvascular complications in patients with hypertension and type 2 diabetes. Diabetes Care. 2000; 23(Suppl 2): B54–B64. [PubMed]
Yu DY, Cringle SJ, Balaratnasingam C, Morgan WH, Yu PK, Su EN. Retinal ganglion cells: Energetics, compartmentation, axonal transport, cytoskeletons and vulnerability. Prog Retin Eye Res. 2013; 36: 217–246. [CrossRef] [PubMed]
Cao D, Yang D, Yu H, et al. Optic nerve head perfusion changes preceding peripapillary retinal nerve fibre layer thinning in preclinical diabetic retinopathy. Clin Exp Ophthalmol. 2019; 47: 219–225. [CrossRef] [PubMed]
Zhang M, Jia F, Li N, et al. Quantitative analysis of the RPC vessel density and the RNFL thickness in patients with type 2 diabetes mellitus by using OCT angiography. Ophthalmic Res. 2021; 64: 951–959. [CrossRef] [PubMed]
Shin YI, Nam KY, Lee SE, et al. Peripapillary microvasculature in patients with diabetes mellitus: An optical coherence tomography angiography study. Sci Rep. 2019; 9: 15814. [CrossRef] [PubMed]
Lee WH, Park JH, Won Y, et al. Retinal microvascular change in hypertension as measured by optical coherence tomography angiography. Sci Rep. 2019; 9: 156. [CrossRef] [PubMed]
Peng Q, Hu Y, Huang M, et al. Retinal neurovascular impairment in patients with essential hypertension: an optical coherence tomography angiography study. Invest Ophthalmol Vis Sci. 2020; 61(8): 42. [CrossRef] [PubMed]
Lee JC, Grisafe DJ, Burkemper B, et al. Intrasession repeatability and intersession reproducibility of peripapillary OCTA vessel parameters in non-glaucomatous and glaucomatous eyes. Br J Ophthalmol. 2021; 105: 1534–1541. [CrossRef] [PubMed]
Kim J, Woo SJ, Ahn J, Park KH, Chung H, Park KH. Long-term temporal changes of peripapillary retinal nerve fiber layer thickness before and after panretinal photocoagulation in severe diabetic retinopathy. Retina. 2012; 32: 2052–2060. [CrossRef] [PubMed]
Simó R, Stitt AW, Gardner TW. Neurodegeneration in diabetic retinopathy: does it really matter? Diabetologia. 2018; 61: 1902–1912. [CrossRef] [PubMed]
Oshitari T, Hanawa K, Adachi-Usami E. Changes of macular and RNFL thicknesses measured by Stratus OCT in patients with early stage diabetes. Eye (Lond). 2009; 23: 884–889. [CrossRef] [PubMed]
Sugimoto M, Sasoh M, Ido M, Wakitani Y, Takahashi C, Uji Y. Detection of early diabetic change with optical coherence tomography in type 2 diabetes mellitus patients without retinopathy. Ophthalmologica. 2005; 219: 379–385. [CrossRef] [PubMed]
Nousome D, Mckean-Cowdin R, Richter GM, et al. Retinal nerve fiber layer thickness in healthy eyes of Black, Chinese, and Latino Americans: a population-based multiethnic study. Ophthalmology. 2021; 128: 1005–1015. [CrossRef] [PubMed]
Li X, Xie J, Zhang L, et al. Identifying microvascular and neural parameters related to the severity of diabetic retinopathy using optical coherence tomography angiography. Invest Ophthalmol Vis Sci. 2020; 61(5): 39. [CrossRef] [PubMed]
Pang R, Peng J, Cao K, et al. Association between contrast sensitivity function and structural damage in primary open-angle glaucoma. Br J Ophthalmol. 2024; 108: 801–806. [CrossRef] [PubMed]
Rosenkranz SC, Gutmann L, Silemek ACH, et al. Visual function resists early neurodegeneration in the visual system in primary progressive multiple sclerosis. J Neurol Neurosurg Psychiatry. 2023; 94: 924–933. [CrossRef] [PubMed]
Lecleire-Collet A, Tessier LH, Massin P, et al. Advanced glycation end products can induce glial reaction and neuronal degeneration in retinal explants. Br J Ophthalmol. 2005; 89: 1631–1633. [CrossRef] [PubMed]
Tso MO, Fine BS. Electron microscopic study of human papilledema. Am J Ophthalmol. 1976; 82: 424–434. [CrossRef] [PubMed]
Dong Y, Guo X, Arsiwala-Scheppach LT, et al. Association of optical coherence tomography and optical coherence tomography angiography retinal features with visual function in older adults. JAMA Ophthalmol. 2022; 140: 809–817. [PubMed]
van Dijk HW, Verbraak FD, Kok PHB, et al. Early neurodegeneration in the retina of type 2 diabetic patients. Invest Ophthalmol Vis Sci. 2012; 53: 2715–2719. [CrossRef] [PubMed]
Song Y, Lai Z, Ding K, Sun Y, Zeng L. Peripapillary vessels density is closely related to cerebral white matter hyperintensities: An OCTA study. PLoS One. 2024; 19(10): e0312534. [CrossRef] [PubMed]
Langner SM, Terheyden JH, Geerling CF, et al. Structural retinal changes in cerebral small vessel disease. Sci Rep. 2022; 12: 9315. [CrossRef] [PubMed]
Cooper LS, Wong TY, Klein R, et al. Retinal microvascular abnormalities and MRI-defined subclinical cerebral infarction: the Atherosclerosis Risk in Communities Study. Stroke. 2006; 37: 82–86. [CrossRef] [PubMed]
Llanas S, Linderman RE, Chen FK, Carroll J. Assessing the use of incorrectly scaled optical coherence tomography angiography images in peer-reviewed studies: a systematic review. JAMA Ophthalmol. 2020; 138: 86–94. [CrossRef] [PubMed]
Linderman R, Salmon AE, Strampe M, Russillo M, Khan J, Carroll J. Assessing the accuracy of foveal avascular zone measurements using optical coherence tomography angiography: segmentation and scaling. Transl Vis Sci Technol. 2017; 6(3): 16. [CrossRef] [PubMed]
Sampson DM, Gong P, An D, et al. Axial length variation impacts on superficial retinal vessel density and foveal avascular zone area measurements using optical coherence tomography angiography. Invest Ophthalmol Vis Sci. 2017; 58: 3065–3072. [CrossRef] [PubMed]
Dutt DDCS, Yazar S, Charng J, Mackey DA, Chen FK, Sampson DM. Correcting magnification error in foveal avascular zone area measurements of optical coherence tomography angiography images with estimated axial length. Eye Vis. 2022; 9: 29. [CrossRef]
Liu J, Li H, Zhou Y, et al. Deep learning-based estimation of axial length using macular optical coherence tomography images. Front Med. 2023; 10: 1308923. [CrossRef]
Fotedar R, Wang JJ, Burlutsky G, et al. Distribution of axial length and ocular biometry measured using partial coherence laser interferometry (IOL Master) in an older white population. Ophthalmology. 2010; 117: 417–423. [CrossRef] [PubMed]
Figure 1.
 
Representative images of RPC vasculature enface, structural enface, RPC perfusion heat map, and RNFL thickness map by DR stage.
Figure 1.
 
Representative images of RPC vasculature enface, structural enface, RPC perfusion heat map, and RNFL thickness map by DR stage.
Figure 2.
 
Peripapillary microvascular metrics across stages of DR. (A) Capillary perfusion density, mean and by quadrant, across DR stages. (B) CFI, mean and by quadrant, across DR stages. *P < 0.05.
Figure 2.
 
Peripapillary microvascular metrics across stages of DR. (A) Capillary perfusion density, mean and by quadrant, across DR stages. (B) CFI, mean and by quadrant, across DR stages. *P < 0.05.
Table 1.
 
Clinical and Demographic Features of the Studied Cohort
Table 1.
 
Clinical and Demographic Features of the Studied Cohort
Table 2.
 
Radial Peripapillary Capillary Optical Coherence Tomography Angiography Metrics Among All Eyes
Table 2.
 
Radial Peripapillary Capillary Optical Coherence Tomography Angiography Metrics Among All Eyes
Table 3.
 
Mixed Effects Multivariable Regression Analysis, Controlling for Age and Hypertension, Comparing Radial Peripapillary Capillary OCTA Metrics Across Different Stages of Diabetic Retinopathy
Table 3.
 
Mixed Effects Multivariable Regression Analysis, Controlling for Age and Hypertension, Comparing Radial Peripapillary Capillary OCTA Metrics Across Different Stages of Diabetic Retinopathy
Table 4.
 
Mixed Effects Multivariable Regression Analysis Controlling for Age and Lens Status Evaluating Associations Between Radial Peripapillary Capillary OCTA Metrics and Visual Acuity (VA) and Contrast Sensitivity (CS) Outcome Measures
Table 4.
 
Mixed Effects Multivariable Regression Analysis Controlling for Age and Lens Status Evaluating Associations Between Radial Peripapillary Capillary OCTA Metrics and Visual Acuity (VA) and Contrast Sensitivity (CS) Outcome Measures
×
×

This PDF is available to Subscribers Only

Sign in or purchase a subscription to access this content. ×

You must be signed into an individual account to use this feature.

×