Using perimetric data to estimate ganglion cell loss for detecting progression of glaucoma: a comparison of models

Examining the concordance of retinal ganglion cell counts generated using measures of structure and function

PURPOSE

There are several indirect methods used to estimate retinal ganglion cell (RGC) count in an individual eye, but there is limited information as to the agreement between these methods. In this work, RGC receptive field (RGC-RF) count underlying a spot stimulus (0.43°, Goldmann III) was calculated and compared using three different methods.

METHODS

RGC-RF count was calculated at a retinal eccentricity of 2.32 mm for 44 healthy adult participants (aged 18-58 years, refractive error -9.75 DS to +1.75 DS) using: (i) functional measures of achromatic peripheral grating resolution acuity (PGRA), (ii) structural measures of RGC-layer thickness (OCT-model, based on the method outlined by Raza and Hood) and (iii) scaling published histology density data to simulate a global expansion in myopia (Histology-Balloon).

RESULTS

Whilst average RGC-RF counts from the OCT-model (median 105.3, IQR 99.6-111.0) and the Histology-Balloon model (median 107.5, IQR 97.7-114.6) were similar, PGRA estimates were approximately 65% lower (median 37.7, IQR 33.8-46.0). However, there was poor agreement between all three methods (Bland-Altman 95% limits of agreement; PGRA/OCT: 55.4; PGRA/Histology-Balloon 59.3; OCT/Histology-Balloon: 52.4). High intersubject variability in RGC-RF count was evident using all three methods.

CONCLUSIONS

The lower PGRA RGC-RF counts may be the result of targeting only a specific subset of functional RGCs, as opposed to the coarser approach of the OCT-model and Histology-Balloon, which include all RGCs, and also likely displaced amacrine cells. In the absence of a 'ground truth', direct measure of RGC-RF count, it is not possible to determine which method is most accurate, and each has limitations. However, what is clear is the poor agreement found between the methods prevents direct comparison of RGC-RF counts between studies utilising different methodologies and highlights the need to utilise the same method in longitudinal work.

Effects of Criterion Bias on Perimetric Sensitivity and Response Variability in Glaucoma

Purpose

The purpose of this study was to isolate and quantify the effects of observer response criterion on perimetric sensitivity, response variability, and maximum response probability.

Methods

Twelve people with glaucoma were tested at three locations in the visual field (age = 47-77 years, mean deviation = -0.61 to -14.54 dB, test location Humphrey field analyzer [HFA] sensitivities = 1 to 30 dB). Frequency of seeing (FoS) curves were measured using a method of constant stimuli with two response paradigms: a "yes-no" paradigm similar to static automated perimetry and a criterion-free two interval forced choice (2IFC) paradigm. Comparison measures of sensitivity, maximum response probability, and response variability were derived from the fitted FoS curves.

Results

Sensitivity differences between the tasks varied widely (range = -11.3 dB to 21.6 dB) and did not correlate with visual field sensitivity nor whether the visual field location was in an area of steep sensitivity gradient within the visual field. Due to the wide variation in differences between the methods, there was no significant difference in mean sensitivity between the 2IFC task relative to the yes-no task, but a trend for higher sensitivity (mean = 1.9 dB, SD = 6.0 dB, P = 0.11). Response variability and maximum response probability did not differ between the tasks (P > 0.99 and 0.95, respectively).

Conclusions

Perimetric sensitivity estimates are demonstrably altered by observer response criterion but the effect varies widely and unpredictably, even within a single test. Response bias should be considered a factor in perimetric test variability and when comparing sensitivities to nonperimetric data.

Translational Relevance

The effect of response criterion on perimetric response variability varies widely and unpredictably, even within a single test.

Increased Equivalent Input Noise in Glaucomatous Central Vision: Is it Due to Undersampling of Retinal Ganglion Cells?

Purpose

Recent evidence shows that macular damage is common even in early stages of glaucoma. Here we investigated whether contrast sensitivity loss in the central vision of glaucoma patients is due to an increase in equivalent input noise (Neq), a decrease in calculation efficiency, or both. We also examined how retinal undersampling resulting from loss of retinal ganglion cells (RGCs) may affect Neq and calculation efficiency.

Methods

This study included 21 glaucoma patients and 23 age-matched normally sighted individuals. Threshold contrast for orientation discrimination was measured with a sinewave grating embedded in varying levels of external noise. Data were fitted to the linear amplifier model (LAM) to factor contrast sensitivity into Neq and calculation efficiency. We also correlated macular RGC counts estimated from structural (spectral-domain optical coherence tomography) and functional (standard automated perimetry Swedish interactive thresholding algorithm 10-2) data with either Neq or efficiency. Furthermore, using analytical and computer simulation approach, the relative effect of retinal undersampling on Neq and efficiency was evaluated by adding the RGC sampling module into the LAM.

Results

Compared with normal controls, glaucoma patients exhibited a significantly larger Neq without significant difference in efficiency. Neq was significantly correlated with Pelli-Robson contrast sensitivity and macular RGC counts. The results from analytical derivation and model simulation further demonstrated that Neq can be expressed as a function of internal noise and retinal sampling.

Conclusions

Our results showed that equivalent input noise is significantly elevated in glaucomatous vision, thereby impairing foveal contrast sensitivity. Our findings further elucidated how undersampling at the retinal level may increase equivalent input noise.

Topographic correlation and asymmetry analysis of ganglion cell layer thinning and the retinal nerve fiber layer with localized visual field defects

Purpose

To evaluate the accuracy of the measurement of the ganglion cell layer (GCL) of the posterior pole analysis (PPA) software of the Spectralis spectral-domain (SD) optical coherence tomography (OCT) device (Heidelberg Engineering, Inc., Heidelberg, Germany), the asymmetry of paired GCL sectors, the total retinal thickness asymmetry (RTA), and the peripapillary retinal nerve fiber layer (pRNFL) test to discriminate between healthy, early and advanced glaucoma eyes.

Methods

Three hundred eighteen eyes of 161 individuals with reliable visual fields (VF) were enrolled in this study. All participants were examined using the standard posterior pole and the pRNFL protocols of the Spectralis OCT device. VF impairment was graded in hemifields, and the GCL sectors were correlated with this damage. Thicknesses of each GCL, the GCL map deviation asymmetry and the pRNFL were compared between control and glaucomatous eyes. The area under the receiver operating characteristic curve (AUC) of these analyses was assessed.

Results

Fourteen of the 16 sectors of the GCL and pRNFL were significantly thinner in eyes with glaucoma than in control eyes (p<0.006). Similarly, the GCL map deviation showed a significant difference between these eyes and both the control eyes as well as the eyes with early glaucoma (p = 0.001 and p = 0.039, respectively). The highest values of AUC to diagnose both early and advanced glaucoma corresponded to the average pRNFL analysis and the GCL map deviation (AUC>0.823, p<0.040 and AUC>0.708, p<0.188, respectively).

Conclusions

Although 16 central sectors of the GCL observed with PPA showed good correlation with VF damage, the pRNFL and the GCL map deviation were more effective for discrimination of glaucomatous damage.

The Topographical Relationship between Visual Field Loss and Peripapillary Retinal Nerve Fibre Layer Thinning Arising from Long-Term Exposure to Vigabatrin

BACKGROUND

The antiepileptic drug vigabatrin is associated with characteristic visual field loss (VAVFL) and thinning of the peripapillary retinal nerve fibre layer (PPRNFL); however, the relationship is equivocal.

OBJECTIVE

The aim of this study was to determine the function-structure relationship associated with long-term exposure to vigabatrin, thereby improving the risk/benefit analysis of the drug.

METHODS

A cross-sectional observational design identified 40 adults who had received long-term vigabatrin for refractory seizures, who had no evidence of co-existing retino-geniculo-cortical visual pathway abnormality, and who had undergone a standardized protocol of perimetry and of optical coherence tomography (OCT) of the PPRNFL. Vigabatrin toxicity was defined as the presence of VAVFL. The function-structure relationship for the superior and inferior retinal quadrants was evaluated by two established models applicable to other optic neuropathies.

RESULTS

The function-structure relationship for each model was consistent with an optic neuropathy. PPRNFL thinning, expressed in micrometres, asymptoted at an equivalent visual field loss of worse than approximately - 10.0 dB, thereby preventing assessment of more substantial thinning. Transformation of the outcomes to retinal ganglion cell soma and axon estimates, respectively, resulted in a linear relationship.

CONCLUSIONS

Functional and structural abnormality is strongly related in individuals with vigabatrin toxicity and no evidence of visual pathway comorbidity, thereby implicating retinal ganglion cell dysfunction. OCT affords a limited measurement range compared with perimetry: severity cannot be directly assessed when the PPRNFL quadrant thickness is less than approximately 65 µm, depending on the tomographer. This limitation can be overcome by transformation of thickness to remaining axons, an outcome requiring input from perimetry.

Imaging and quantifying ganglion cells and other transparent neurons in the living human retina

Ganglion cells (GCs) are fundamental to retinal neural circuitry, processing photoreceptor signals for transmission to the brain via their axons. However, much remains unknown about their role in vision and their vulnerability to disease leading to blindness. A major bottleneck has been our inability to observe GCs and their degeneration in the living human eye. Despite two decades of development of optical technologies to image cells in the living human retina, GCs remain elusive due to their high optical translucency. Failure of conventional imaging-using predominately singly scattered light-to reveal GCs has led to a focus on multiply-scattered, fluorescence, two-photon, and phase imaging techniques to enhance GC contrast. Here, we show that singly scattered light actually carries substantial information that reveals GC somas, axons, and other retinal neurons and permits their quantitative analysis. We perform morphometry on GC layer somas, including projection of GCs onto photoreceptors and identification of the primary GC subtypes, even beneath nerve fibers. We obtained singly scattered images by: (i) marrying adaptive optics to optical coherence tomography to avoid optical blurring of the eye; (ii) performing 3D subcellular image registration to avoid motion blur; and (iii) using organelle motility inside somas as an intrinsic contrast agent. Moreover, through-focus imaging offers the potential to spatially map individual GCs to underlying amacrine, bipolar, horizontal, photoreceptor, and retinal pigment epithelium cells, thus exposing the anatomical substrate for neural processing of visual information. This imaging modality is also a tool for improving clinical diagnosis and assessing treatment of retinal disease.

Using perimetric data to estimate ganglion cell loss for detecting progression of glaucoma: a comparison of models

PURPOSE

Models relating perimetric sensitivities to ganglion cell numbers have been proposed for combining structural and functional measures from patients with glaucoma. Here we compared seven models for ability to differentiate progressing and stable patients, testing the hypothesis that the model incorporating local spatial scale would have the best performance.

METHODS

The models were compared for the United Kingdom Glaucoma Treatment Study (UKGTS) data for the right eyes of 489 patients recently diagnosed with glaucoma. The SITA 24-2 program was utilised for perimetry and Stratus OCT fast scanning protocol for thickness of circumpapillary retinal nerve fibre layer (RNFL). The first analysis defined progression in terms of decline in RNFL thickness. The highest and lowest quintiles (22 subjects per group) were identified for change in thickness of inferior temporal (IT), superior temporal (ST), and global RNFL (μm year^-1 ); a two-way anova was used to look for differences between the models in ability to discriminate the two quintiles. The second analysis defined a 'progression group' as those who were flagged by the UKGTS criteria as having progressive loss in perimetric sensitivity, and a 'no progression' group as those with rate of change in Mean Deviation (MD) closest to 0 dB year^-1 (87 subjects per group). The third analysis characterised variability of retinal ganglion cell (RGC) models for the two groups in the second analysis, using the standard deviation of residuals from linear regression of ganglion cell number over time to compute Coefficient of Variation (CoV).

RESULTS

The first analysis produced a negative result because the three anovas found no effect of model or interaction of model and group (F_6,294 < 3.1, p > 0.08). There was an effect of group only for the anova with the ST sector (F_6,294 = 12.2, p < 0.001). The second analysis also produced a negative result, because ROC areas were in the range 0.69-0.72 for all models. The third analysis found that even when variability in MD was low, the CoV was so large that test-retest variation could include 100% loss of ganglion cells.

CONCLUSIONS

Two very different approaches for testing the hypothesis both gave a negative result. For all seven ganglion cell models, rates of ganglion cell loss were highly affected by fluctuations in height of the hill of vision. Methods for reducing effects of between-visit variability are needed in order to assess progression by relating perimetric sensitivities and ganglion cell numbers.