Longitudinal Macular Structure-Function Relationships in Glaucoma

Detecting Fast Progressors: Comparing a Bayesian Longitudinal Model to Linear Regression for Detecting Structural Changes in Glaucoma

PURPOSE

Demonstrate that a novel Bayesian hierarchical spatial longitudinal (HSL) model identifies macular superpixels with rapidly deteriorating ganglion cell complex (GCC) thickness more efficiently than simple linear regression (SLR).

DESIGN

Prospective cohort study.

SETTING

Tertiary Glaucoma Center.

SUBJECTS

One hundred eleven eyes (111 patients) with moderate to severe glaucoma at baseline and ≥4 macular optical coherence tomography scans and ≥2 years of follow-up.

OBSERVATION PROCEDURE

Superpixel-patient-specific GCC slopes and their posterior variances in 49 superpixels were derived from our latest Bayesian HSL model and Bayesian SLR. A simulation cohort was created with known intercepts, slopes, and residual variances in individual superpixels.

MAIN OUTCOME MEASURES

We compared HSL and SLR in the fastest progressing deciles on (1) proportion of superpixels identified as significantly progressing in the simulation study and compared to SLR slopes in cohort data; (2) root mean square error (RMSE), and SLR/HSL RMSE ratios.

RESULTS

Cohort- In the fastest decile of slopes per SLR, 77% and 80% of superpixels progressed significantly according to SLR and HSL, respectively. The SLR/HSL posterior SD ratio had a median of 1.83, with 90% of ratios favoring HSL. Simulation- HSL identified 89% significant negative slopes in the fastest progressing decile vs 64% for SLR. SLR/HSL RMSE ratio was 1.36 for the fastest decile of slopes, with 83% of RMSE ratios favoring HSL.

CONCLUSION

The Bayesian HSL model improves the estimation efficiency of local GCC rates of change regardless of underlying true rates of change, particularly in fast progressors.

Study of the ganglion cell complex of the macula by optical coherence tomography in the diagnosis of glaucoma progression

INTRODUCTION

The aim of this work is to evaluate the usefulness of the study of the ganglion cell complex of the macula using the OCT technique to estimate the progression of glaucoma according to its severity.

MATERIAL AND METHODS

This is a retrospective cross-sectional study. It includes 205 eyes of 131 patients with glaucoma or ocular hypertension followed for a mean of 5.7 years. The parameters and rates of three tests have been analyzed using the progression software of each instrument: visual field, optical coherence tomography (OCT) in the ganglion cell complex of the macula and in the nerve fiber layer of the optic nerve. The results of each test, the concordance between them and how they differ according to severity stage have been evaluated.

RESULTS

Visual field classifies more cases of progression in moderate-advanced glaucoma, while in mild glaucoma its capacity is limited. Optic nerve fiber layer OCT classifies more cases of progression in mild glaucoma than in moderate-advanced glaucoma, as it is artifacted by the floor effect. OCT of the macular ganglion cell complex is the test that classifies more cases of progression and has the highest agreement with visual field, regardless of severity.

CONCLUSION

In both mild and moderate-advanced glaucoma, OCT of the macula ganglion cell complex may be a better biomarker of progression than OCT of the macula ganglion cell complex.

Assessment of the relationship between structural and functional tests in patients with idiopathic intracranial hypertension

PURPOSE

To assess the relationship between structural and functional tests in mild and moderate idiopathic intracranial hypertension (IIH).

METHODS

Patients with mild and moderate IIH and a control group were enrolled. Best-corrected visual acuity (BCVA), macular ganglion cell layer (MGCL) thickness, peripapillary retinal nerve fiber layer (pp RNFL) thickness, perimetric mean deviation (MD), and photopic negative responses (PhNR) of the electroretinogram were recorded. The associations between structural (pp RNFL and MGCL thickness) and functional (PhNR amplitude, MD and BCVA) parameters were assessed.

RESULTS

154 eyes from 78 subjects (74 eyes from IIH patients and 80 eyes from healthy subjects) were included in this comparative observational study. The MGCL thickness, VA, pp RNFL, and PhNR base-to-trough (BT) amplitude were significantly worse in moderate IIH. The BCVA and MD were associated with MGCL thickness only in moderate IIH. The relationship between MD and MGCL thickness started when MD fell below -5.7 dB.

CONCLUSIONS

The association between functional and structural parameters varies between mild and moderate IIH. The MD and MGCL thickness outperformed in assessing disease severity in mild and moderate IIH, respectively. The association between MD and MGCL thickness could be considered in IIH severity categorization.

Impacts of Chronic Kidney Disease on Retinal Neurodegeneration: A Cross-Cohort Analysis

PURPOSE

To assess the cross-sectional and longitudinal associations between chronic kidney disease (CKD) and ganglion cell-inner plexiform layer (GCIPL) thickness in a UK Biobank population and a Chinese cohort.

DESIGN

Prospective observational cohort study and cross-sectional study.

METHOD

This study included 23,014 individuals without neurodegenerative diseases from the UK Biobank, and 3 years of annual follow-up data of 2197 individuals from a Chinese cohort. Three groups were defined by estimated glomerular filtration rate (eGFR) based on serum creatinine classifying CKD severity as no CKD, mild CKD, and moderate to severe CKD (MS-CKD). GCIPL thickness, measured using optical coherence tomography, was analyzed through linear regression over time to determine its decline rate in micrometers per year. Linear regression models were used to assess the correlation between renal function and both the baseline GCIPL thickness and the GCIPL decline rate.

RESULTS

The cross-sectional analysis in a largely white population showed that poorer renal function negatively correlated with GCIPL thickness with a mean of 0.15 µm thinner (95% confidence interval [CI] -0.30 to -0.01; P = .042) in mild CKD and 0.83 µm thinner (95% CI -1.34 to -0.32; P = .001) in MS-CKD compared with that of control subjects without CKD. Longitudinal analysis in the Chinese cohort showed that the GCIPL decreased more rapidly in persons with poorer renal function. After correcting for all confounding factors, the rate of GCIPL thinning was 0.30 µm/year (95% CI -0.41 to -0.19; P < .001) more in the mild CKD group and 0.52 µm/year (95% CI -0.79 to -0.26; P < .001) more in the MS-CKD group compared with control subjects without CKD. This relationship also occurred in individuals with diabetes or hypertension.

CONCLUSIONS

Poor renal function was associated with a lower baseline GCIPL thickness in the UK population and a faster decline rate in Chinese participants. However, the detailed underlying mechanisms still need further exploration.

A Bayesian Hierarchical Spatial Longitudinal Model Improves Estimation of Local Macular Rates of Change in Glaucomatous Eyes

Purpose

Demonstrate that a novel Bayesian hierarchical spatial longitudinal (HSL) model improves estimation of local macular ganglion cell complex (GCC) rates of change compared to simple linear regression (SLR) and a conditional autoregressive (CAR) model.

Methods

We analyzed GCC thickness measurements within 49 macular superpixels in 111 eyes (111 patients) with four or more macular optical coherence tomography scans and two or more years of follow-up. We compared superpixel-patient-specific estimates and their posterior variances derived from the latest version of a recently developed Bayesian HSL model, CAR, and SLR. We performed a simulation study to compare the accuracy of intercept and slope estimates in individual superpixels.

Results

HSL identified a significantly higher proportion of significant negative slopes in 13/49 superpixels and a significantly lower proportion of significant positive slopes in 21/49 superpixels than SLR. In the simulation study, the median (tenth, ninetieth percentile) ratio of mean squared error of SLR [CAR] over HSL for intercepts and slopes were 1.91 (1.23, 2.75) [1.51 (1.05, 2.20)] and 3.25 (1.40, 10.14) [2.36 (1.17, 5.56)], respectively.

Conclusions

A novel Bayesian HSL model improves estimation accuracy of patient-specific local GCC rates of change. The proposed model is more than twice as efficient as SLR for estimating superpixel-patient slopes and identifies a higher proportion of deteriorating superpixels than SLR while minimizing false-positive detection rates.

Translational Relevance

The proposed HSL model can be used to model macular structural measurements to detect individual glaucoma progression earlier and more efficiently in clinical and research settings.

Comparison of the Circumpapillary Structure-Function and Vasculature-Function Relationships at Different Glaucoma Stages Using Longitudinal Data

Purpose

This study investigated the global and regional correlations between longitudinal structure-function (S-F) and vasculature-function (V-F) data using circumpapillary retinal nerve fiber layer thickness (cpRNFLT) measurements from optical coherence tomography (OCT), circumpapillary vessel density (cpVD) from OCT angiography (OCTA), and the corresponding visual field mean sensitivities at different glaucoma stages.

Methods

A total of 107 eyes from 107 glaucoma patients with progressive visual field (VF) changes followed up for an average of 3.33 ± 1.39 years were enrolled, including early-to-moderate (51 eyes) and advanced (56 eyes) stages. The rates of longitudinal change in the VF mean deviation (MD), cpRNFLT, and cpVD were evaluated using linear mixed-effects models and compared between different glaucoma stages. Longitudinal global and regional S-F and V-F relationships were assessed by repeated measures correlation analysis by glaucoma stage.

Results

No significant differences were found in the rates of VF MD and cpVD changes (P > 0.05) between the two glaucoma stage groups. CpRNFLT decreased more rapidly in the early-to-moderate stage group (P < 0.001) in which significant longitudinal global and regional correlations were found in both S-F and V-F relationships (all P < 0.05), except for the nasal sector. Significant global and regional correlations were only found in V-F relationship in advanced stage cases (all P < 0.05).

Conclusions

Significant longitudinal V-F relationships exist globally and regionally regardless of glaucoma stage but no longitudinal S-F relationship is present in advanced glaucoma. Longitudinal follow-up of cpVD parameters may be useful for monitoring glaucomatous VF progression at all disease stages.

Deep Learning Estimation of 10-2 Visual Field Map Based on Macular Optical Coherence Tomography Angiography Measurements

PURPOSE

To develop deep learning (DL) models estimating the central visual field (VF) from optical coherence tomography angiography (OCTA) vessel density (VD) measurements.

DESIGN

Development and validation of a deep learning model.

METHODS

A total of 1051 10-2 VF OCTA pairs from healthy, glaucoma suspects, and glaucoma eyes were included. DL models were trained on en face macula VD images from OCTA to estimate 10-2 mean deviation (MD), pattern standard deviation (PSD), 68 total deviation (TD) and pattern deviation (PD) values and compared with a linear regression (LR) model with the same input. Accuracy of the models was evaluated by calculating the average mean absolute error (MAE) and the R2 (squared Pearson correlation coefficients) of the estimated and actual VF values.

RESULTS

DL models predicting 10-2 MD achieved R2 of 0.85 (95% confidence interval [CI], 74-0.92) for 10-2 MD and MAEs of 1.76 dB (95% CI, 1.39-2.17 dB) for MD. This was significantly better than mean linear estimates for 10-2 MD. The DL model outperformed the LR model for the estimation of pointwise TD values with an average MAE of 2.48 dB (95% CI, 1.99-3.02) and R2 of 0.69 (95% CI, 0.57-0.76) over all test points. The DL model outperformed the LR model for the estimation of all sectors.

CONCLUSIONS

DL models enable the estimation of VF loss from OCTA images with high accuracy. Applying DL to the OCTA images may enhance clinical decision making. It also may improve individualized patient care and risk stratification of patients who are at risk for central VF damage.

Cost-efficient and glaucoma-specifical model by exploiting normal OCT images with knowledge transfer learning

Monitoring the progression of glaucoma is crucial for preventing further vision loss. However, deep learning-based models emphasize early glaucoma detection, resulting in a significant performance gap to glaucoma-confirmed subjects. Moreover, developing a fully-supervised model is suffering from insufficient annotated glaucoma datasets. Currently, sufficient and low-cost normal OCT images with pixel-level annotations can serve as valuable resources, but effectively transferring shared knowledge from normal datasets is a challenge. To alleviate the issue, we propose a knowledge transfer learning model for exploiting shared knowledge from low-cost and sufficient annotated normal OCT images by explicitly establishing the relationship between the normal domain and the glaucoma domain. Specifically, we directly introduce glaucoma domain information to the training stage through a three-step adversarial-based strategy. Additionally, our proposed model exploits different level shared features in both output space and encoding space with a suitable output size by a multi-level strategy. We have collected and collated a dataset called the TongRen OCT glaucoma dataset, including pixel-level annotated glaucoma OCT images and diagnostic information. The results on the dataset demonstrate our proposed model outperforms the un-supervised model and the mixed training strategy, achieving an increase of 5.28% and 5.77% on mIoU, respectively. Moreover, our proposed model narrows performance gap to the fully-supervised model decreased by only 1.01% on mIoU. Therefore, our proposed model can serve as a valuable tool for extracting glaucoma-related features, facilitating the tracking progression of glaucoma.

Prediction of Central Visual Field Measures From Macular OCT Volume Scans With Deep Learning

Purpose

Predict central 10° global and local visual field (VF) measurements from macular optical coherence tomography (OCT) volume scans with deep learning (DL).

Methods

This study included 1121 OCT volume scans and 10-2 VFs from 289 eyes (257 patients). Macular scans were used to estimate 10-2 VF mean deviation (MD), threshold sensitivities (TS), and total deviation (TD) values at 68 locations. A three-dimensional (3D) convolutional neural network based on the 3D DenseNet121 architecture was used for prediction. We compared DL predictions to those from baseline linear models. We carried out 10-fold stratified cross-validation to optimize generalizability. The performance of the DL and baseline models was compared based on correlations between ground truth and predicted VF measures and mean absolute error (MAE; ground truth - predicted values).

Results

Average (SD) MD was -9.3 (7.7) dB. Average (SD) correlations between predicted and ground truth MD and MD MAE were 0.74 (0.09) and 3.5 (0.4) dB, respectively. Estimation accuracy deteriorated with worsening MD. Average (SD) Pearson correlations between predicted and ground truth TS and MAEs for DL and baseline model were 0.71 (0.05) and 0.52 (0.05) (P < 0.001) and 6.5 (0.6) and 7.5 (0.5) dB (P < 0.001), respectively. For TD, correlation (SD) and MAE (SD) for DL and baseline models were 0.69 (0.02) and 0.48 (0.05) (P < 0.001) and 6.1 (0.5) and 7.8 (0.5) dB (P < 0.001), respectively.

Conclusions

Macular OCT volume scans can be used to predict global central VF parameters with clinically relevant accuracy.

Translational Relevance

Macular OCT imaging may be used to confirm and supplement central VF findings using deep learning.

Comparing Rates of Change in Moderate to Advanced Glaucoma: Retinal Nerve Fiber Layer Versus Bruch Membrane Opening-Minimum Rim Width

PURPOSE

To compare rates of change (RoC) of peripapillary retinal nerve fiber layer (RNFL) and Bruch membrane opening-based minimum rim width (BMO-MRW) thickness in moderate-to-advanced glaucoma.

DESIGN

Prospective cohort study.

METHODS

Longitudinal optical coherence tomography (OCT) optic nerve head volume scans of 113 eyes of 113 glaucoma patients with moderate-to-advanced or central damage were exported. This study estimated and compared global and sectoral RoC with linear mixed effects models and simple linear regression (SLR) of RNFL and BMO-MRW thickness. Permutation analyses were used to test significance of RoC in the SLR model. It also compared longitudinal signal-to-noise ratios (LSNR) defined as RoC divided by residual standard deviation (SD) between the two groups.

RESULTS

Mean (SD) follow-up and median (IQR) OCT scan sessions were 5.2 (1.3) years and 10 (8-11), respectively. Baseline average (SD) visual field mean deviation was -9.2 (5.8) dB. Based on SLR, a higher proportion of significant negative RNFL RoC was observed compared to BMO-MRW in the inferotemporal (35% vs 20%; P = .015) and inferonasal (42% vs 17%; P < .001) sectors. Permutation analyses also demonstrated a higher proportion of worsening RNFL RoC than BMO-MRW in the inferotemporal (P = .026) and inferonasal (P < .001) sectors along with overall lower positive RoC. Longitudinal signal-to-noise ratios for RNFL were significantly more negative than for BMO-MRW globally, and in the inferotemporal, inferonasal, and superonasal sectors (P ≤ .01).

CONCLUSIONS

Longitudinal RNFL OCT measurements are more likely to detect structural change and demonstrate better LSNR compared with BMO-MRW in eyes with central or moderate-to-advanced glaucoma damage at baseline.

Relationships of Macular Functional Impairment With Structural and Vascular Changes According to Glaucoma Severity

Purpose

To determine the pointwise relationships of central visual field (VF) defects with macular ganglion cell loss and macular vessel density (VD) loss during various stages of glaucoma.

Methods

Eyes with primary open-angle glaucoma (POAG) were subjected to optical coherence tomography (OCT) and OCT angiography (OCTA) to evaluate macular ganglion cell layer (GCL) thickness and macular VD in the superficial and deep vascular complexes (SVC and DVC). OCT, OCTA, and VF locations were matched after correcting for retinal ganglion cell (RGC) displacement. Pointwise correlations of GCL thickness and VDs of the SVC and DVC with central VF sensitivity (VFS) were evaluated by Pearson's correlation analysis and compared in eyes with early and advanced POAG by Meng's test.

Results

Of the 100 eyes, 52 and 48 were classified as early and advanced POAG. Macular VD showed overall better correlation with central VFS than GCL thickness in both the early and advanced groups. SVC density showed the strongest correlation with central VFS in all groups (R = 0.327 in early group, R = 0.325 in advanced group, all P < 0.001). Although DVC density showed better correlation with VFS (R = 0.311) than GCL thickness (R = 0.212) in the early group (P < 0.001), the correlation was comparable in the advanced group (R = 0.199 and 0.176, respectively, P = 0.254).

Conclusions

After adjustment for RGC displacement, macular SVC density was better correlated with central VFS than macular GCL thickness in both early and advanced POAG. Macular DVC density showed better correlation with VFS than GCL thickness in early but not in advanced POAG, indicating that DVC loss may be involved in early central VF loss.

Importance of Pattern Standard Deviation of Humphrey 10-2 Visual Field to Evaluate Central Visual Function in Patients with Early-Stage Glaucoma

To explore various parameters that can evaluate the central visual impairment in patients with early-stage glaucoma, we included patients into a study with central visual impairments with an MD value greater than -6.0 dB on the 24-2 VF test. A possible association between structural parameters acquired by OCT and functional parameters of VF and PERG was determined. A total of 70 eyes of patients with suspected glaucoma or NTG underwent VF, OCT, and PERG examinations. The patients were classified into two groups according to the MD of the 24-2 VF test. We used Pearson correlation analysis to evaluate the relationships between GCIPL thickness/RNFL thickness and visual functional parameters, such as PERG and perimetry. Linear regression analyses were conducted to evaluate the significant factors affecting the PSD of VF 10-2. In the low MD group, the P50 amplitude presented significant correlations (r = 0.346, p = 0.048) with GCIPL thickness. In the correlation analysis of the high MD group, it was found that only the PSD of 10-2 uniquely presented borderline significant correlations with GCIPL thickness (r = -0.327, p = 0.055), and no other functional parameter showed significant correlation. Univariate and multivariate analyses revealed that GCIPL thickness was significantly associated with a PSD of 10-2 VF (p < 0.001 and 0.013, respectively). Among various parameters, the P50 amplitude and 10-2 PSD demonstrated statistically borderline significant structure-function relationships with GCIPL thickness in early-stage glaucoma.

GABA decrease is associated with degraded neural specificity in the visual cortex of glaucoma patients

Glaucoma is an age-related neurodegenerative disease of the visual system, affecting both the eye and the brain. Yet its underlying metabolic mechanisms and neurobehavioral relevance remain largely unclear. Here, using proton magnetic resonance spectroscopy and functional magnetic resonance imaging, we investigated the GABAergic and glutamatergic systems in the visual cortex of glaucoma patients, as well as neural specificity, which is shaped by GABA and glutamate signals and underlies efficient sensory and cognitive functions. Our study shows that among the older adults, both GABA and glutamate levels decrease with increasing glaucoma severity regardless of age. Further, our study shows that the reduction of GABA but not glutamate predicts the neural specificity. This association is independent of the impairments on the retina structure, age, and the gray matter volume of the visual cortex. Our results suggest that glaucoma-specific decline of GABA undermines neural specificity in the visual cortex and that targeting GABA could improve the neural specificity in glaucoma.

Vision-Related Quality of Life Among Healthy, Preperimetric Glaucoma, and Perimetric Glaucoma Patients

PURPOSE

To investigate the association of vision-related quality of life (VRQOL) with the central visual field and macular ganglion cell complex (GCC) thickness in healthy control participants, patients with preperimetric glaucoma, and patients with perimetric glaucoma.

DESIGN

Retrospective cross-sectional study.

METHODS

A total of 39 healthy, 34 preperimetric glaucoma, and 145 perimetric glaucoma patients completed the 25-item National Eye Institute Visual Function Questionnaire (NEI-VFQ). A linear mixed-effect models was used to investigate the association between the glaucoma stage as measured by binocular 10-2 visual field mean sensitivity (VFMS) and GCC thickness with the Rasch-calibrated NEI-VFQ score.

RESULTS

A total of 436 eyes from 218 participants (mean age = 67.2 [95% CI = 65.1 to 69.2] years) were enrolled. VRQOL calculated by the NEI-VFQ Rasch-calibrated score was worst for patients with perimetric glaucoma (50.7 [95% CI = 47.2 to 54.2]), followed by patients with preperimetric glaucoma (41.2 [95% CI = 34.5 to 47.9]) and healthy controls (29.3 [95% CI = 24.0 to 34.7]. Worse VRQOL had a moderate association with a worse global binocular 10-2 VFMS (-3.4 [95% CI = -5.0 to -1.9] dB per 1 score; P < .001; adjusted R² = 0.27), but not with a thinner global GCC in the better eye (-0.1 [95% CI = -0.2 to 0.1] µm per 1 score; P =.0485; adjusted R² = 0.17).

CONCLUSIONS

These findings suggest that patients with perimetric and preperimetric glaucoma have worse VRQOL than patients with healthy eyes. As compared to macular thickness measurements, the central visual field is more strongly associated with VRQOL and may better help to identify patients in need of intervention.

Retinal Neurodegeneration and Visual Acuity Decline in Patients with Chronic Kidney Disease

INTRODUCTION

Chronic kidney disease (CKD) has been associated with accelerated retinal neurodegeneration. The purpose of this study is to evaluate the association between retinal neurodegeneration and the best-corrected visual acuity (BCVA) decline in patients with CKD.

METHODS

Post hoc analysis of two prospective studies. Patients with CKD stage ≥ 3 were enrolled. Macular thickness, peripapillary retinal nerve fiber layer (pRNFL) thickness, and macular ganglion cell complex (GCC) thickness were measured by optical coherence tomography. Eyes were classified into three groups: Group 1, no GCC defect; Group 2, GCC defect confined to parafoveal area; and Group 3, GCC defects extending beyond the parafoveal area. Each group was matched for age, sex, axial length, lens status, and cataract grading.

RESULTS

A total of 120 eyes (40 eyes in each group) from 120 patients (age 63.0 ± 10.3 years) were included. The logMAR BCVA was 0.076 ± 0.101, 0.100 ± 0.127, and 0.196 ± 0.191 in Group 1, 2, and 3, respectively. Group 3, but not Group 2, had a significantly worse BCVA than Group 1. In simple linear regression, parafoveal inner retinal thickness, pRNFL thickness, presence of pRNFL defect, GCC thickness, GCC global loss volume, GCC focal loss volume, and GCC defect extending beyond parafoveal area were associated with BCVA. Central subfield retinal thickness (CRT), parafoveal full retinal thickness, and parafoveal outer retinal thickness were not associated with BCVA. In backward stepwise linear regression, age and GCC defects extending beyond the parafoveal area were factors associated with BCVA. Moreover, GCC defect extending beyond parafoveal area was connected with worse BCVA in both phakic and pseudophakic subgroups.

CONCLUSIONS

GCC defect extending beyond parafoveal area could be an independent biomarker associated with decreased BCVA in patients with CKD. However, macular thinning measured by CRT or parafoveal full retinal thickness might have low discriminative power in determining BCVA.

Effect of Corneal Hysteresis on the Rates of Microvasculature Loss in Glaucoma

PURPOSE

To investigate the association between corneal hysteresis (CH) and rates of optic nerve head whole image capillary density (wiCD) loss over time in open-angle glaucoma (OAG).

DESIGN

Observational cohort.

PARTICIPANTS

One hundred seventy-four eyes (122 OAG and 52 glaucoma suspect eyes) from 112 patients over more than 2 years and 4 visits or more.

METHODS

Baseline CH measurements were acquired with the Ocular Response Analyzer. Linear mixed-effect models were designed to investigate the effect of CH, average intraocular pressure (IOP) during follow-up, and baseline visual field (VF) mean deviation (MD) on the rates of wiCD loss and circumpapillary retinal nerve fiber layer (cpRNFL) thinning over time, while adjusting for confounders. Interaction between CH or baseline MD and average IOP during follow-up were included in final models to evaluate the effect of baseline MD or average IOP during follow-up on structural changes for different values of CH.

MAIN OUTCOME MEASURE

Effect of CH, IOP, and baseline MD on the rates of wiCD loss and cpRNFL thinning over time.

RESULTS

The average follow-up time was 3.9 years. In the multivariable model, non-Black race, higher average IOP during follow-up, lower baseline CH, lower baseline VF MD, and higher numbers of IOP-lowering medications were associated with faster rates of wiCD loss over time. For CH values 6 mmHg and 12 mmHg, every 1-mmHg increase in average IOP during follow-up was associated with 0.23% per year faster and 0.07% per year slower rates of wiCD loss over time, respectively. While every 1-mmHg decrease in CH was associated with 1.89% per year faster rate of wiCD loss for MD of -12 dB, it was associated with 0.81% per year faster rate of wiCD loss for MD of -3 dB.

CONCLUSION

Lower CH values were significantly associated with faster rates of wiCD loss over time. In eyes with lower CH, both higher average IOP during follow-up and more severe glaucoma damage at baseline were associated with faster rates of wiCD loss and cpRNFL thinning. These results support CH as a useful parameter for risk assessment of glaucoma progression.

FINANCIAL DISCLOSURE(S)

Proprietary or commercial disclosure may be found after the references.

Glaucoma progression. Clinical practice guide

OBJECTIVE

To provide general recommendations that serve as a guide for the evaluation and management of glaucomatous progression in daily clinical practice based on the existing quality of clinical evidence.

METHODS

After defining the objectives and scope of the guide, the working group was formed and structured clinical questions were formulated following the PICO (Patient, Intervention, Comparison, Outcomes) format. Once all the existing clinical evidence had been independently evaluated with the AMSTAR 2 (Assessment of Multiple Systematic Reviews) and Cochrane "Risk of bias" tools by at least two reviewers, recommendations were formulated following the Scottish Intercollegiate Guideline network (SIGN) methodology.

RESULTS

Recommendations with their corresponding levels of evidence that may be useful in the interpretation and decision-making related to the different methods for the detection of glaucomatous progression are presented.

CONCLUSIONS

Despite the fact that for many of the questions the level of scientific evidence available is not very high, this clinical practice guideline offers an updated review of the different existing aspects related to the evaluation and management of glaucomatous progression.

Measurement of the Inner Macular Layers for Monitoring of Glaucoma: Confounding Effects of Age-Related Macular Degeneration

OBJECTIVE

To investigate the confounding effect of nonexudative age-related macular degeneration (AMD), specifically drusen and outer retinal atrophy, on the architecture and automated segmentation of the inner retinal layers as measured with OCT.

DESIGN

Observational cross-sectional study.

SUBJECTS

Two hundred sixty-three consecutive eyes with nonexudative AMD were identified through a retrospective chart review. Exclusion criteria were a diagnosis of glaucoma or glaucoma suspect, other retinal pathology affecting the macula, axial length > 26.5 mm or spherical equivalent less than -6 diopters, any other optic nerve or neurologic disorders, or poor image quality.

METHODS

Drusen were automatically segmented on macular OCT B-scans with a publicly available and validated deep learning approach. Automated segmentation of the inner plexiform layer (IPL)/inner nuclear layer (INL) boundary was carried out with the device's proprietary software.

MAIN OUTCOME MEASURES

Quality of segmentation of the IPL/INL boundary as a function of drusen size and presence of inner retinal layer displacement in the area of macular pathology (drusen or atrophy).

RESULTS

One hundred twenty-five eyes (65 patients) met the inclusion criteria. Drusen size varied between 16 and 272 μm (mean, 118 μm). Automated segmentation had a 22% chance of failure if the drusen height was between 145 and 185 μm and was most likely to fail with drusen heights above 185 μm. When drusen height was normalized by total retinal thickness, segmentation failed 36% of the time when the drusen to total retinal thickness ratio was 0.45 or above. Images were likely to show displacement of inner retinal layers with drusen heights above 176 μm and a normalized drusen height ratio of 0.5 or higher. Eighty-seven percent of images with outer retinal atrophy displayed incorrect segmentation.

CONCLUSIONS

Outer retinal diseases can alter the retinal topography and affect the segmentation accuracy of the inner retinal layers. Large drusen may cause segmentation error and compression of the inner macular layers. Geographic atrophy confounds automated segmentation in a high proportion of eyes. Clinicians should be cognizant of the effects of outer retinal disease on the inner retinal layer measurements when interpreting the results of macular OCT imaging in patients with glaucoma.

Advanced glaucoma. Clinical practice guideline

OBJECTIVE

To present an update clinical practice guideline that serve as a guide for the detection, evaluation and treatment of adults patients with advanced glaucoma.

METHODS

After defining the objectives and scope of the guide, the working group was formed and structured clinical questions were formulated following the PICO (Patient, Intervention, Comparison, Outcomes) format. Once all the existing clinical evidence had been independently evaluated with the AMSTAR 2 (Assessment of Multiple systematic Rewiews) and Cochrane "Risk of bias" tools by at least two reviewers, recommendations were formulated following the Scottish Intercollegiate methodology. Guideline Network (SIGN).

RESULTS

Recommendations with their corresponding levels of evidence that may be useful in the diagnosis, monitoring and treatment of adults patients with advanced glaucoma.

CONCLUSIONS

Despite the fact that for many of the questions the level of scientific evidence available is not very high, this clinical practice guideline offers an updated review of the different existing aspects related to the evaluation and management of advanced glaucoma.

Association Between Rate of Ganglion Cell Complex Thinning and Rate of Central Visual Field Loss

Importance

Whether rapid ganglion cell complex (GCC) thinning during an initial follow-up period is associated with rates of central visual field loss over time is unclear but important to understand because risk of glaucoma progression can help guide treatment intensity.

Objective

To investigate the association between the rate of GCC thinning during initial follow-up and the rate of central visual field loss.

Design, Setting, and Participants

This retrospective cohort study assessed patients older than 18 years with glaucoma at a tertiary glaucoma center who were followed up from June 18, 2014, to January 11, 2019. Data analysis for the current study was undertaken in March 2022.

Main Outcomes and Measures

Initial rates of GCC thinning were obtained from global GCC thickness values of the first 3 optical coherence tomography (OCT) scans. Rates of central visual field loss were assessed as the change in central (10-2) visual field mean deviation during the 4.7-year follow-up period by univariable and multivariable linear mixed-effects models. Eyes were categorized as slow (>-1 μm/y) or fast (≤-1 μm/y) progressors based on rates of GCC thinning.

Results

The cohort consisted of 202 eyes of 139 patients (mean [SD] age, 68.7 [10.0] years; 72 male [51.8%]); 44 African American patients (31.7%), 13 Asian patients (9.4%), 80 White patients (57.6%), and 2 patients who identified as other race and ethnicity (1.4%) were analyzed. The rate of GCC change was -0.56 μm/y (95% CI, -0.66 to -0.46 μm/y) during a mean initial follow-up of 1.8 years (95% CI, 1.7-2.0 years). A total of 163 eyes (80.7%) were slow OCT progressors, and 39 (19.3%) were fast OCT progressors, with rates of GCC thinning of -0.3 μm/y (95% CI, -0.4 to -0.2 μm/y) and -1.6 μm/y (-1.8 to -1.3 μm/y), respectively. The rates of 10-2 visual field mean deviation worsening among slow and fast OCT progressors were -0.10 dB/y (95% CI, -0.16 to 0.00 dB/y) and -0.34 dB/y (95% CI, -0.51 to -0.16 dB/y), respectively (difference, -0.26 dB/y; 95% CI, -0.45 to -0.07 dB/y; P = .008).

Conclusions and Relevance

In this cohort study, rapid GCC thinning during an initial follow-up period was associated with faster rates of central visual field decline. These findings support use of longitudinal macular OCT scans assisting clinical decision-making for glaucoma and also may guide possible intensification of therapy in high-risk patients.

Comparison of Ganglion Cell Layer and Ganglion Cell/Inner Plexiform Layer Measures for Detection of Early Glaucoma

PURPOSE

To test the hypothesis that macular ganglion cell layer (GCL) measurements detect early glaucoma with higher accuracy than ganglion cell/inner plexiform layer (GCIPL) thickness measurements.

DESIGN

Cross-sectional study.

PARTICIPANTS

The first cohort included 58 glaucomatous eyes with visual field mean deviation (MD) ≥ -6 dB and 125 normal eyes. The second cohort included 72 glaucomatous and 73 normal/glaucoma suspect (GS) eyes with scans able to create GCL/GCIPL deviation maps.

METHODS

In the first cohort, 8 × 8 GCL and GCIPL grids were exported and 5 superior and inferior sectors were defined. Global and sectoral GCL and GCIPL measures were used to predict glaucoma. In the second cohort, proportions of scan areas with abnormal (< 5% and < 1% cutoffs) and supernormal (> 95% and > 99% cutoffs) thicknesses on deviation maps were calculated. The extents of GCL and GCIPL abnormal areas were used to predict glaucoma.

MAIN OUTCOME MEASURES

Extents of abnormal GCL/GCIPL regions and areas under receiver operating characteristic curves (AUROC) for prediction of glaucoma were compared between GCL or GCIPL measures.

RESULTS

The average ± standard deviation MDs were -3.7 ± 1.6 dB and -2.7 ± 1.8 dB in glaucomatous eyes in the first and second cohorts, respectively. Global GCIPL thickness measures (central 18° × 18° macular region) performed better than GCL for early detection of glaucoma (AUROC, 0.928 vs. 0.884, respectively; P = 0.004). Superior and inferior sector 3 thickness measures provided the best discrimination with both GCL and GCIPL (inferior GCL AUROC, 0.860 vs. GCIPL AUROC, 0.916 [P = 0.001]; superior GCL AUROC, 0.916 vs. GCIPL AUROC, 0.900 [P = 0.24]). The extents of abnormal GCL regions at a 1% cutoff in the central elliptical area were 17.5 ± 22.2% and 6.4 ± 10.8% in glaucomatous and normal/GS eyes, respectively, versus 17.0 ± 22.2% and 5.7 ± 10.5%, respectively, for GCIPL (P = 0.06 for GCL and 0.002 for GCIPL). The extents of GCL and GCIPL supernormal regions were mostly similar in glaucomatous and normal eyes. The best performance for prediction of glaucoma in the second cohort was detected at a P value of < 1% within the entire scan for both GCL and GCIPL (AUC, 0.681 vs. 0.668, respectively; P = 0.29).

CONCLUSIONS

Macular GCL and GCIPL thicknesses are equivalent for identifying early glaucoma with current OCT technology. This is likely explained by limitations of inner macular layer segmentation and concurrent changes within the inner plexiform layer in early glaucoma.

Concordance of Objectively Detected Retinal Nerve Fiber Bundle Defects in En Face OCT Images with Conventional Structural and Functional Changes in Glaucoma

PURPOSE

To assess how objectively detected defects in retinal nerve fiber bundle (RNFB) reflectance on en face OCT images relate to circumpapillary retinal nerve fiber layer thickness (cpRNFLT) and visual field defects.

DESIGN

Cross-sectional study.

PARTICIPANTS

Sixteen participants with early glaucoma and 29 age-matched healthy controls, of whom 22 had usable en face images for the establishment of normative levels of RNFB reflectance.

METHODS

All the participants underwent cpRNFLT scans, visual field examination, and wide-field OCT. En face reflectivity was assessed objectively using the Summary of Multiple Anatomically Adjusted Slabs method. En face defects were deemed concordant with cpRNFLT when they had at least 1 cpRNFLT point with P < 0.01, within ± 15° of the predicted insertion on the optic disc. Visual fields were examined using custom suprathreshold perimetry and SITA Standard 24-2. For each visual field location, the corresponding reflectance was deemed abnormal if any en face superpixel within ± 1° was abnormal. The overall, positive, and negative agreements were measured in each participant.

MAIN OUTCOME MEASURES

Proportion of concordant defects between en face reflectance analysis and cpRNFLT (%) as well as overall, positive, and negative agreements between en face reflectance analysis and visual field results.

RESULTS

Most en face abnormalities had concordant cpRNFLT defects in the mapped sector (median proportion concordant, 0.85; interquartile range, 0.74-0.95). In eyes with glaucoma, a median of 8.1% (range, 2.4%-23.7%) and 14.9% (range, 3.5%-29.1%) locations showed corresponding en face and visual field defects using 24-2 and custom perimetry, respectively. Both the perimetric strategies had moderate-to-good raw agreement with en face analysis (0.66-0.68), with stronger agreement on normal findings than on defects (0.77-0.78 and 0.4-0.44).

CONCLUSIONS

Objectively extracted reflectance defects showed strong concordance with conventional cpRNFLT damage and good agreement with perimetry, which could be enhanced by further minimization of image artifacts.

Superpixel-Based Optic Nerve Head Segmentation Method of Fundus Images for Glaucoma Assessment

Glaucoma disease is the second leading cause of blindness in the world. This progressive ocular neuropathy is mainly caused by uncontrolled high intraocular pressure. Although there is still no cure, early detection and appropriate treatment can stop the disease progression to low vision and blindness. In the clinical practice, the gold standard used by ophthalmologists for glaucoma diagnosis is fundus retinal imaging, in particular optic nerve head (ONH) subjective/manual examination. In this work, we propose an unsupervised superpixel-based method for the optic nerve head (ONH) segmentation. An automatic algorithm based on linear iterative clustering is used to compute an ellipse fitting for the automatic detection of the ONH contour. The tool has been tested using a public retinal fundus images dataset with medical expert ground truths of the ONH contour and validated with a classified (control vs. glaucoma eyes) database. Results showed that the automatic segmentation method provides similar results in ellipse fitting of the ONH that those obtained from the ground truth experts within the statistical range of inter-observation variability. Our method is a user-friendly available program that provides fast and reliable results for clinicians working on glaucoma screening using retinal fundus images.

Comparison of Ganglion Cell Layer and Inner Plexiform Layer Rates of Change in Suspected and Established Glaucoma

PURPOSE

We compared ganglion cell layer (GCL) and inner plexiform layer (IPL) rates of change (RoC) in patients with glaucoma suspect (GS) and established glaucoma (EG) to test the hypothesis that IPL thickness changes would occur earlier than GCL changes in eyes with early damage.

DESIGN

Prospective, cohort study.

METHODS

A total of 64 GS eyes (46 patients) and 112 EG eyes (112 patients) with ≥2 years of follow-up and ≥3 macular optical coherence tomography scans were included. GCL and IPL superpixel thickness measurements were exported. A Bayesian hierarchical model with random intercepts/slopes and random residual variances was fitted to estimate RoC in individual superpixels. Normalized RoC and proportions of superpixels with significantly negative and positive GCL and IPL RoC were compared within the groups.

RESULTS

The average (SD) follow-up time and number of scans were 3.5 (0.7) years and 4.2 (1.0), respectively, in the GS group and 3.6 (0.4) years and 7.3 (1.1) in the EG group. Mean (SD) normalized RoC was faster for GCL than IPL (-0.69 [0.05] vs -0.33 [0.04]) in the GS group, whereas it was faster for IPL (-0.47 [0.03] vs -0.28 [0.02]) in EG eyes. GCL RoC were significantly negative in 24 of 36 superpixels compared with 8 of 36 for IPL (P < .001) in GS eyes. In the EG group, 23 of 36 superpixels had significant negative IPL RoC compared with 13 of 36 superpixels for GCL (P = .006).

CONCLUSIONS

GCL thickness is more likely to demonstrate change over time compared with IPL in glaucoma suspects. There is no evidence of preferential IPL thinning in eyes with suspected early glaucoma damage.

Longitudinal Structure-Function Relationship between Macular Vessel Density and Thickness and Central Visual Field in Early Glaucoma

PURPOSE

To investigate the relationship of longitudinal changes in macular vessel density (VD) from OCT angiography and in ganglion cell complex (GCC) from OCT with central visual field (VF) in eyes with early glaucoma.

DESIGN

Observational cohort.

PARTICIPANTS

A total of 95 eyes, 37 preperimetric and 58 with early glaucoma (24-2 VF mean deviation [MD] ≥ -6 decibels), with an average follow-up of 3.8 years and 5.3 visits, were included.

METHODS

Whole-image VD (wiVD) and whole-image GCC (wiGCC) and parafoveal scans, as well as localized regions of interest (LROIs), hemiretinae of whole images, and superior, inferior, temporal, and nasal sectors of parafoveal maps, were matched with central VF locations. Age-adjusted rates of change of VD, GCC, mean sensitivity of VF locations, and 10-2 VF MD were calculated using linear mixed-effect models. Normalized rates of change were calculated for comparison of change rates in wiVD and wiGCC.

MAIN OUTCOME MEASURES

Structure-function (SF) correlations of VD and GCC with central VF measurement change rates and comparison of different correlations of SF relationships after bootstrapping the difference of the correlation coefficients.

RESULTS

Vessel density loss and GCC thinning demonstrated significant correlations with central VF damage, globally and with most LROIs. The SF correlation (r, 95% confidence interval [CI]) between wiVD and 10-2 VF MD change rates was 0.42 [0.24, 0.58], whereas it was 0.27 [0.08, 0.45] between wiGCC and 10-2 VF MD changes rates (all P < 0.05). In contrast to GCC thinning, VD loss in the parafoveal sectors demonstrated significant correlations with central VF damage in inferior and temporal sectors. Differences in the relationship of SF with central VF damage were not significant between VD loss and GCC thinning. The mean (95% CI) normalized change rates of wiVD (-7.40 [-7.71 to 7.09] %/year) was faster than that of wiGCC (-2.39 [-2.94 to 1.84] %/year) (P < 0.05).

CONCLUSIONS

Rates of VD loss and GCC thinning are associated with central VF loss over time. Assessment of both macular VD and GCC thickness should be considered for evaluation of glaucoma progression.

Comparison of the Structure-Function Relationship Between Advanced Primary Open Angle Glaucoma and Normal Tension Glaucoma

PURPOSE

The aim was to investigate and compare the characteristics of visual field (VF) defects in primary open angle glaucoma (POAG) and normal-tension glaucoma (NTG) with advanced glaucomatous damage and to determine whether the structure-function relationships found in advanced glaucoma differ based on their glaucoma classification.

PATIENTS AND METHODS

Ninety-seven eyes of 97 patients (59 eyes with POAG and 38 eyes with NTG) with advanced glaucoma were included in this cross-sectional study. Scores at each test point of the 30-2 VF total deviation map were recorded, and average values at each test point were point-wise compared between the groups. Peripapillary retinal nerve fiber layer (RNFL) and macular thickness (total, RNFL, ganglion cell layer, and inner plexiform layer thickness) were measured. The structure-function relationship based on the map of Garway-Heath was determined and compared between the 2 groups.

RESULTS

At advanced stage of glaucoma, POAG eyes demonstrated more diffusely distributed VF defects, whereas NTG eyes had more severe VF defects at the superior nasal quadrant, showing increased asymmetry. Overall, peripapillary RNFL, macular ganglion cell layer, and macular inner plexiform layer thickness showed good relationships with 30-2 VF parameters in both groups. However, in total macula and macular RNFL thickness, the structure-function relationships tended to show different characteristics depending on the glaucoma classification; NTG eyes showed overall better relationships.

CONCLUSIONS

In advanced glaucoma, differences in patterns of VF damage were found between POAG and NTG eyes. Conventional peripapillary RNFL and macular measurements showed generally good performance for estimating functional status, particularly in NTG eyes.

Association Between Ganglion Cell Complex Thinning and Vision-Related Quality of Life in Glaucoma

Importance

Faster structural changes may be associated with worse vision-related quality of life in patients with glaucoma.

Objectives

To evaluate the association between the rate of ganglion cell complex thinning and the Vision Function Questionnaire in glaucoma.

Design, Setting, and Participants

This retrospective analysis of a longitudinal cohort was designed in October 2021. Patients were enrolled from the Diagnostic Innovations in Glaucoma Study and the African Descent and Glaucoma Evaluation Study. Two hundred thirty-six eyes of 118 patients with diagnosed or suspected glaucoma were followed up with imaging for a mean of 4.1 years from September 2014 to March 2020.

Main Outcomes and Measures

The Vision Function Questionnaire was evaluated using the 25-item National Eye Institute Visual Function at the last follow-up visit. Ganglion cell complex thickness was derived from macular optical coherence tomography scans and averaged within 3 circular areas (3.4°, 5.6°, and 6.8° from the fovea) and superior and inferior hemiregions. Linear mixed-effects models were used to investigate the association between the rate of ganglion cell complex thinning and Rasch-calibrated Vision Function Questionnaire score.

Results

The mean (SD) age was 73.2 (8.7) years, 65 participants (55.1%) were female, and 53 participants (44.9%) were African American. Race was self-reported by the participants. Mean composite Rasch-calibrated National Eye Institute Visual Function Questionnaire score was 50.3 (95% CI, 45.9-54.6). A faster annual rate of global ganglion cell complex thinning in the better eye was associated with a higher disability reflected by the composite National Eye Institute Visual Function Questionnaire score (-15.0 [95% CI, -28.4 to -1.7] per 1 μm faster; P = .03). When stratified by degrees from the fovea, the 5.6° and 6.8° areas were associated with the composite National Eye Institute Visual Function Questionnaire Rasch-calibrated score (-14.5 [95% CI, -27.0 to -2.0] per 1 μm faster; R2 = 0.201; P = .03; and -23.7 [95% CI, -45.5 to -1.9] per 1 μm faster; R2 = 0.196; P = .02, respectively), and -8.0 (95% CI, -16.8 to 0.8) per 1 μm faster for the 3.4° area (R2 = 0.184; P = .07) after adjusting for confounding factors.

Conclusions and Relevance

These findings suggest that faster and sectoral central location of ganglion cell complex thinning provides useful information in determining the risk of vision-related quality of life in glaucoma. Monitoring macular structure may be useful for determining the risk of functional impairment in glaucoma.

Elucidating macular structure-function correlations in glaucoma

Correlation between structural data from optical coherence tomography and functional data from the visual field may be suboptimal because of poor mapping of OCT measurement locations to VF stimuli. We tested the hypothesis that stronger structure–function correlations in the macula can be achieved with fundus-tracking perimetery, by precisely mapping OCT measurements to VF sensitivity at the same location. The conventional 64 superpixel (3° × 3°) OCT grid was mapped to VF sensitivities averaged in 40 corresponding VF units with standard automated perimetry (conventional mapped approach, CMA) in 38 glaucoma patients and 10 healthy subjects. Similarly, a 144 superpixel (2° × 2°) OCT grid was mapped to each of the 68 locations with fundus-tracking perimetry (localized mapped approach, LMA). For each approach, the correlation between sensitivity at each VF unit and OCT superpixel was computed. Vector maps showing the maximum correlation between each VF unit and OCT pixel was generated. CMA yielded significantly higher structure–function correlations compared to LMA. Only 20% of the vectors with CMA and < 5% with LMA were within corresponding mapped OCT superpixels, while most were directed towards loci with structural damage. Measurement variability and patterns of structural damage more likely impact correlations compared to precise mapping of VF stimuli.

Multivariate Longitudinal Modeling of Macular Ganglion Cell Complex: Spatiotemporal Correlations and Patterns of Longitudinal Change

Purpose

To investigate spatiotemporal correlations among ganglion cell complex (GCC) superpixel thickness measurements and explore underlying patterns of longitudinal change across the macular region.

Design

Longitudinal cohort study.

Subjects

One hundred eleven eyes from 111 subjects from the Advanced Glaucoma Progression Study with ≥ 4 visits and ≥ 2 years of follow-up.

Methods

We further developed our proposed Bayesian hierarchical model for studying longitudinal GCC thickness changes across macular superpixels in a cohort of glaucoma patients. Global priors were introduced for macular superpixel parameters to combine data across superpixels and better estimate population slopes and intercepts.

Main Outcome Measures

Bayesian residual analysis to inspect cross-superpixel correlations for subject random effects and residuals. Principal component analysis (PCA) to explore underlying patterns of longitudinal macular change.

Results

Average (standard deviation [SD]) follow-up and baseline 10-2 visual field mean deviation were 3.6 (0.4) years and -8.9 (5.9) dB, respectively. Superpixel-level random effects and residuals had the greatest correlations with nearest neighbors; correlations were higher in the superior than in the inferior region and strongest among random intercepts, followed by random slopes, residuals, and residual SDs. PCA of random intercepts showed a first large principal component (PC) across superpixels that approximated a global intercept, a second PC that contrasted the superior and inferior macula, and a third PC, contrasting inner and nasal superpixels with temporal and peripheral superpixels. PCs for slopes, residual SDs, and residuals were remarkably similar to those of random intercepts.

Conclusions

Introduction of cross-superpixel random intercepts and slopes is expected to improve estimation of population and subject parameters. Further model enhancement may be possible by including cross-superpixel random effects and correlations to address spatiotemporal relationships in longitudinal data sets.

Correlation Between Enlargement of Retinal Nerve Fiber Defect Angle in En Face Imaging and Visual Field Progression

Purpose

Retinal nerve fiber layer defects (RNFLDs) become enlarged with glaucoma progression. We measured the RNFLD angle and investigated whether it was correlated with deterioration of the visual field in patients with glaucoma.

Methods

This study included 84 eyes of 84 patients with open-angle glaucoma (mean deviation [MD] = −6.51 ± 5.91 dB, follow-up period = 2.82 ± 0.74 years) with the RNFLDs, who underwent en face swept-source optical coherence tomography (SS-OCT) wide scans (12 × 9 mm) at least 6 times. The RNFLD angle was measured as the intersection between the RNFLD and a circle centered on the disc with a radius half the distance between the disc and the fovea. Slopes for the RNFLD angle, macular ganglion cell layer thickness (GCCT), and circumpapillary RNFL thickness (cpRNFLT) were compared with the MD slope, as measured with the Humphrey field analyzer 24-2 program.

Results

The correlation coefficients with MD slope were −0.67 for the RNFLD angle slope (P < 0.001), 0.15 for the macular GCCT slope (P = 0.163), and 0.04 for the cpRNFLT slope (P = 0.719). The RNFLD angle tended to increase as the number of disc hemorrhage occurrences increased (rs = 0.31, P = 0.004). The RNFLD angle slope also had good predictive power for glaucoma progression (area under the receiver operating characteristic curve = 0.88, 95% confidence interval = 0.81–0.95).

Conclusions

We found that the RNFLD angle slope was more closely associated with the MD slope than were other OCT parameters. This suggests that measurement of the RNFLD angle with en face OCT images could be effective in evaluating glaucoma progression.

Translational Relevance

Our study provides a method for monitoring glaucoma progression with SS-OCT.

Prediction of Retinal Ganglion Cell Counts Considering Various Displacement Methods From OCT-Derived Ganglion Cell-Inner Plexiform Layer Thickness

Purpose

To compare various displacement models using midget retinal ganglion cell to cone (mRGC:C) ratios and to determine viability of estimating RGC counts from optical coherence tomography (OCT)–derived ganglion cell–inner plexiform layer (GCIPL) measurements.

Methods

Four Drasdo model variations were applied to macular visual field (VF) stimulus locations: (1) using meridian-specific Henle fiber length along the stimulus circumference; (2) using meridian-specific differences in RGC receptive field and counts along the stimulus circumference; (3) per method (2), averaged across principal meridians; and (4) per method (3), with the stimulus center displaced only. The Sjöstrand model was applied (5) along the stimulus circumference and (6) to the stimulus center only. Eccentricity-dependent mRGC:C ratios were computed over displaced areas, with comparisons to previous models using sum of squares of the residuals (SSR) and root mean square error (RMSE). RGC counts estimated from OCT-derived ganglion cell layer (GCL) and GCIPL measurements, from 143 healthy participants, were compared using Bland–Altman analyses.

Results

Methods 1, 2, and 5 produced mRGC:C ratios most consistent with previous models (SSR 3.82, 4.07, and 3.02; RMSE 0.22, 0.23, and 0.20), while central mRGC:C ratios were overestimated by method 3 and underestimated by methods 4 and 6. RGC counts predicted from GCIPL measurements were within 16% of GCL-based counts, with no notable bias with increasing RGC counts.

Conclusions

Sjöstrand displacement and meridian-specific Drasdo displacement applied to VF stimulus circumferences produce mRGC:C ratios consistent with previous models. RGC counts can be estimated from OCT-derived GCIPL measurements.

Translational Relevance

Implementing appropriate displacement methods and deriving RGC estimates from relevant OCT parameters enables calculation of the number of RGCs responding to VF stimuli from commercial instrumentation.

Ganglion Cell Complex: The Optimal Measure for Detection of Structural Progression in the Macula

PURPOSE

To test the hypothesis that macular ganglion cell complex (GCC) thickness from optical coherence tomography (OCT) provides a stronger change signal regardless of glaucoma severity compared with other macular measures.

DESIGN

Prospective cohort study.

METHODS

Eyes were from 112 patients with moderate to severe glaucoma at baseline from a tertiary glaucoma center. In each 3° × 3° macular superpixel, a hierarchical Bayesian random intercept and slope model with random residual variance was fit to longitudinal full macular thickness (FMT), outer retina layers, GCC, ganglion cell-inner plexiform layer (GCIPL), and ganglion cell layer (GCL) measurements. We estimated population- and individual-level slopes and intercepts. Proportions of substantial worsening and improving superpixel slopes were compared between layers and in superpixels with mild to moderate vs severe damage (total deviation of corresponding visual field location ≥ -8 vs < -8 dB).

RESULTS

Mean (SD) follow-up time and baseline 10-2 visual field mean deviation were 3.6 (0.4) years and -8.9 (5.9) dB, respectively. FMT displayed the highest proportion of significant negative slopes (1932/3519 [54.9%]), followed by GCC (1286/3519 [36.5%]), outer retina layers (1254/3519 [35.6%]), (GCIPL) (1075/3518 [30.6%]), and (GCL) (698/3518 [19.8%]). Inner macular measures detected less worsening in the severe glaucoma group; yet GCC (223/985 [22.6%]) identified the highest proportion (GCIPL: 183/985 [18.6%]; GCL: 106/985 [10.8%]). Proportions of positive rates were small and comparable among all measures.

CONCLUSIONS

GCC is the optimal macular measure for detection of structural change in eyes with moderate to severe glaucoma. Although a higher proportion of worsening superpixels was observed for FMT, a large portion of FMT change could be attributed to changes in outer retina layers.

Detection of Longitudinal Ganglion Cell/Inner Plexiform Layer Change: Comparison of Two Spectral-Domain Optical Coherence Tomography Devices

PURPOSE

We compared rates of change of macular ganglion cell/inner plexiform (GCIPL) thickness and proportion of worsening and improving rates from 2 optical coherence tomography (OCT) devices in a cohort of eyes with glaucoma.

DESIGN

Longitudinal cohort study.

METHODS

In a tertiary glaucoma clinic we evaluated 68 glaucoma eyes with ≥2 years of follow-up and ≥4 OCT images. Macular volume scans from 2 OCT devices were exported, coregistered, and segmented. Global and sectoral GCIPL data from the central 4.8 × 4.0-mm region were extracted. GCIPL rates of change were estimated with linear regression. Permutation analyses were used to control specificity with the 2.5 percentile cutoff point used to define "true" worsening. Main outcome measures included differences in global/sectoral GCIPL rates of change between 2 OCT devices and the proportion of negative vs positive rates of change (P < .05).

RESULTS

Average (standard deviation) 24-2 visual field mean deviation, median (interquartile range) follow-up time, and number of OCT images were -9.4 (6.1) dB, 3.8 (3.3-4.2) years, and 6 (5-8), respectively. GCIPL rates of thinning from Spectralis OCT were faster (more negative) compared with Cirrus OCT; differences were significant in superonasal (P = .03) and superotemporal (P = .04) sectors. A higher proportion of significant negative rates was observed with Spectralis OCT both globally and in inferotemporal/superotemporal sectors (P < .04). Permutation analyses confirmed the higher proportion of global and sectoral negative rates of change with Spectralis OCT (P < .001).

CONCLUSIONS

Changes in macular GCIPL were detected more frequently on Spectralis' longitudinal volume scans than those of Cirrus OCT. OCT devices are not interchangeable with regard to detection of macular structural progression.

Estimating Ganglion Cell Complex Rates of Change With Bayesian Hierarchical Models

Purpose

Develop a hierarchical longitudinal regression model for estimating local rates of change of macular ganglion cell complex (GCC) measurements with optical coherence tomography (OCT).

Methods

We enrolled 112 eyes with four or more macular OCT images and ≥2 years of follow-up. GCC thickness measurements within central 6 × 6 superpixels were extracted from macular volume scans. We fit data from each superpixel separately with several hierarchical Bayesian random-effects models. Models were compared with the Watanabe–Akaike information criterion. For our preferred model, we estimated population and individual slopes and intercepts (baseline thickness) and their correlation.

Results

Mean (SD) follow-up time and median (interquartile range) baseline 24-2 visual field mean deviation were 3.6 (0.4) years and −6.8 (−12.2 to −4.3) dB, respectively. The random intercepts and slopes model with random residual variance was the preferred model. While more individual and population negative slopes were observed in the paracentral and papillomacular superpixels, superpixels in the superotemporal and inferior regions displayed the highest correlation between baseline thickness and rates of change (r = –0.43 to –0.50 for the top five correlations).

Conclusions

A Bayesian linear hierarchical model with random intercepts/slopes and random variances is an optimal initial model for estimating GCC slopes at population and individual levels. This novel model is an efficient method for estimating macular rates of change and probability of glaucoma progression locally.

Translational Relevance

The proposed Bayesian hierarchical model can be applied to various macular outcomes from different OCT devices and to superpixels of variable sizes to estimate local rates of change and progression probability.

Prediction of Visual Field Progression from OCT Structural Measures in Moderate to Advanced Glaucoma

PURPOSE

To test the hypothesis that visual field (VF) progression can be predicted from baseline and longitudinal optical coherence tomography (OCT) structural measurements.

DESIGN

Prospective cohort study.

METHODS

A total of 104 eyes (104 patients) with ≥3 years of follow-up and ≥5 VF examinations were enrolled. We defined VF progression based on pointwise linear regression on 24-2 VF (≥3 locations with slope less than or equal to -1.0 dB/year and P < .01). We used elastic net logistic regression (ENR) and machine learning to predict VF progression with demographics, baseline circumpapillary retinal nerve fiber layer (RNFL), macular ganglion cell/inner plexiform layer (GCIPL) thickness, and RNFL and GCIPL change rates at central 24 superpixels and 3 eccentricities, 3.4°, 5.5°, and 6.8°, from fovea and hemimaculas. Areas-under-ROC curves (AUC) were used to compare models.

RESULTS

Average ± SD follow-up and VF examinations were 4.5 ± 0.9 years and 8.7 ± 1.6, respectively. VF progression was detected in 23 eyes (22%). ENR selected rates of change of superotemporal RNFL sector and GCIPL change rates in 5 central superpixels and at 3.4° and 5.6° eccentricities as the best predictor subset (AUC = 0.79 ± 0.12). Best machine learning predictors consisted of baseline superior hemimacular GCIPL thickness and GCIPL change rates at 3.4° eccentricity and 3 central superpixels (AUC = 0.81 ± 0.10). Models using GCIPL-only structural variables performed better than RNFL-only models.

CONCLUSIONS

VF progression can be predicted with clinically relevant accuracy from baseline and longitudinal structural data. Further refinement of proposed models would assist clinicians with timely prediction of functional glaucoma progression and clinical decision making.

Characterizing and quantifying the temporal relationship between structural and functional change in glaucoma

PURPOSE

To characterize and quantify the temporal relationship between structural and functional change in glaucoma.

METHODS

120 eyes of 120 patients with ocular hypertension or primary open-angle glaucoma were selected from the Diagnostic Innovations in Glaucoma Study or the African Descent and Glaucoma Evaluation Study. Patients had 11 visits, separated by at least 3 months over 5 to 10 years. Each visit had rim area (RA) and mean sensitivity (MS) measurements taken within a 30-day period. The structure-function (SF) relationship was summarized using conventional and modified cross-correlation functions (CCFs), which identified the strongest absolute and positive correlation, respectively. Patients were categorized in one of the following three groups: RA and MS evolved simultaneously (lag = 0), RA preceded MS (lag<0), and MS preceded RA (lag>0). Lagging regression analysis was used to examine the variations of the SF relationship within groups.

RESULTS

The number of participants, mean visit lag, and mean correlation (standard deviation) were, for the conventional and modified CCFs, respectively: lag = 0 [16, 0, 0.53 (0.10) and 16, 0, 0.46 (0.11)]; lag<0 [50, -2.94, 0.51 (0.11) and 55, -3.45, 0.44 (0.12)], and lag>0 [54, 3.35, 0.53 (0.13) and 49, 3.78, 0.45 (0.12)]. A significant difference of the visit lag relation within groups was identified using lagging regression analysis (p<0.0001).

CONCLUSIONS

The strongest relationship between structure and function was obtained at different visit lags in different patients. This finding also suggests that the SF relationship should be addressed at the subject level when using both measurements jointly to model glaucoma progression.

Detection of Glaucoma Deterioration in the Macular Region with Optical Coherence Tomography: Challenges and Solutions

PURPOSE

Macular imaging with optical coherence tomography (OCT) measures the most critical retinal ganglion cells (RGCs) in the human eye. The goal of this perspective is to review the challenges to detection of glaucoma progression with macular OCT imaging and propose ways to enhance its performance.

DESIGN

Perspective with review of relevant literature.

METHODS

Review of challenges and issues related to detection of change on macular OCT images in glaucoma eyes. The primary outcome measures were confounding factors affecting the detection of change on macular OCT images.

RESULTS

The main challenges to detection of structural progression in the macula consist of the magnitude of and the variable amount of test-retest variability among patients, the confounding effect of aging, lack of a reliable and easy-to-measure functional outcome or external standard, the confounding effects of concurrent macular conditions including myopia, and the measurement floor of macular structural outcomes. Potential solutions to these challenges include controlling head tilt or torsion during imaging, estimating within-eye variability for individual patients, improved data visualization, the use of artificial intelligence methods, and the implementation of statistical approaches suitable for multidimensional longitudinal data.

CONCLUSIONS

Macular OCT imaging is a crucial structural imaging modality for assessing central RGCs. Addressing the current shortcomings in acquisition and analysis of macular volume scans can enhance its utility for measuring the health of central RGCs and therefore assist clinicians with timely institution of appropriate treatment.

Local Macular Thickness Relationships between 2 OCT Devices

PURPOSE

To compare local ganglion cell-inner plexiform layer (GCIPL) thickness measurements between 2 OCT devices and to explore factors that may influence the difference in measurements.

DESIGN

Cross-sectional study.

PARTICIPANTS

Sixty-nine glaucoma eyes (63 patients) with evidence of central damage or mean deviation (MD) of -6.0 dB or worse on a 24-2 visual field (VF).

METHODS

Cirrus and Spectralis OCT macular volume scans were exported, data from the central 20° of both OCT devices were centered and aligned, and 50 × 50 arrays of 0.4° × 0.4° superpixels were created. We estimated nonparametric (Spearman's) correlations and used Bland-Altman plots to compare GCIPL thickness measurements between the two OCTs at the superpixel level. Factors that may have influenced the differences between thickness measurements between the two devices were explored with linear mixed models.

MAIN OUTCOME MEASURES

Pooled and individual-eye Spearman's correlation and agreement between thickness measurements from the two devices.

RESULTS

The median 24-2 VF MD was -6.8 dB (interquartile range [IQR], -4.9 to -12.3 dB). The overall pooled Spearman's correlation between the two devices for all superpixels and eyes was 0.97 (P < 0.001). The median within-eye correlation coefficient was 0.72 (IQR, 0.59-0.79). Bland-Altman plots demonstrated a systematic bias in most individual eyes, with Spectralis GCIPL measurements becoming larger than Cirrus measurements with increasing superpixel thickness. The average superpixel thickness and distance to the fovea influenced the thickness difference between the two devices in multivariate models (P < 0.001).

CONCLUSIONS

Local macular thickness measurements from the Spectralis and Cirrus devices are highly correlated, but not interchangeable. Differences in thickness measurements between the two devices are influenced by the location of superpixels and their thickness.

Macular imaging with optical coherence tomography in glaucoma

With the advent of spectral-domain optical coherence tomography, imaging of the posterior segment of the eye can be carried out rapidly at multiple anatomical locations, including the optic nerve head, circumpapillary retinal nerve fiber layer, and macula. There is now ample evidence to support the role of spectral-domain optical coherence tomography imaging of the macula for detection of early glaucoma. Macular spectral-domain optical coherence tomography measurements demonstrate high reproducibility, and evidence on its utility for detection of glaucoma progression is accumulating. We present a comprehensive review of macular spectral-domain optical coherence tomography imaging emerging as an essential diagnostic tool in glaucoma.