Contrast sensitivity perimetry and clinical measures of glaucomatous damage

Predicting perimetric defects from en face maps of retinal nerve fibre layer reflectance

PURPOSE

To develop criteria to predict visual hemifields with deep perimetric defects based on retinal nerve fibre layer (RNFL) reflectance, in a transparent process whose components can be assessed by independent laboratories analysing data from their own small groups.

METHODS

The analysis was carried out in four stages, using three independent groups of patients-30, 33 and 62 participants-with glaucoma and age-similar controls. The first stage used Group 1 to develop a criterion for RNFL reflectance images at 24, 36 or 48 μm below the inner limiting membrane (ILM). The second stage evaluated the criterion using Group 2. The third stage developed a second criterion to improve performance for Groups 1 and 2 combined. The fourth stage evaluated the second criterion with Group 3. Confidence intervals for sensitivity and specificity were then computed by combining results from all three groups.

RESULTS

The first criterion identified all hemifields with deep defects and no hemifields from controls, using a within-eye reference for healthy RNFL. For Group 2, specificity remained high but sensitivity was reduced. The second criterion improved sensitivity by using location-specific reference values. For Group 3, sensitivity remained high but reduced specificity was found. Confidence intervals showed substantial overlap for the two criteria.

CONCLUSIONS

We developed two criteria to identify patients with deep perimetric defects with high specificity and sensitivity. Several improvements are warranted: automated identification of the fovea-disc angle and optic disc locations, evaluation of normal variation in patterns of RNFL thickness, improved segmentation of ILM and major vasculature, reduction of within-eye variability in RNFL reflectance of healthy eyes, assessment of effects of image quality, assessment of effects of comorbidity and effectiveness of other devices.

Online Circular Contrast Perimetry: A Comparison to Standard Automated Perimetry

PURPOSE

The aim was to validate and compare the diagnostic accuracy of a novel 24-degree, 52-loci online circular contrast perimetry (OCCP) application to standard automated perimetry (SAP).

DESIGN

Prospective cohort study.

METHODS

Two hundred and twenty participants (125 normal controls, 95 open angle glaucoma patients) were included. Agreement, correlation, sensitivity, specificity, and area under receiver operating curves (AUC) were compared for parameters of OCCP, SAP, and optical coherence tomography (OCT) for the retinal nerve fiber layer and macular ganglion cell complex inner plexiform layer.

RESULTS

Pointwise sensitivity for OCCP was greater than SAP by 1.02 log units (95% CI: 0.95-1.08); 95% limits of agreement 0.860 to 1.17. Correlation and agreement for global indices and regional zones between OCCP and SAP were strong. OCCP mean deviation (MD) AUC was 0.885±0.08, similar to other instruments' parameters with the highest AUC: SAP MD (0.851±0.08), OCT retinal nerve fiber layer inferior thickness (0.908±0.07), OCT ganglion cell complex inner plexiform layer inferior thickness (0.849±0.08), P>0.05. At best cutoff, OCCP MD sensitivity/specificity were comparable to SAP MD (90/74 vs 94/65%).

CONCLUSIONS

OCCP demonstrates similar perimetric sensitivities to SAP and similar AUC to SAP and OCT in distinguishing glaucoma patients from controls. OCCP holds promise as a glaucoma surveillance and screening tool, with the potential to be utilized for in-clinic and at-home perimetry and expand community testing.

Central contrast sensitivity perimetry discriminates between glaucomatous and non-glaucomatous eyes

Background

Glaucoma is an optic neuropathy which causes irreversible vision loss. Standard perimetry, which is essential for glaucoma diagnosis, can only detect glaucomatous visual filed loss when considerable structural damage has occurred. Contrast sensitivity is one of the visual function tests that is reduced in eyes with glaucoma. It is known to be affected in pre-perimetric stages of glaucoma.

Objective

The objective of this study was to investigate the discriminating ability of central contrast sensitivity perimetry in eyes with and without glaucoma.

Design

The study employed a cross-sectional study design.

Methods

The study participants were made of two groups; eyes diagnosed with glaucoma by an ophthalmologist based on visual field test and optical coherence tomography (OCT) and age- and sex-matched controls who were declared free from glaucoma. Static contrast sensitivity (CS) was measured in the central 10° of visual field using a custom psychophysical test.

Results

There were 45 eyes with glaucoma and 45 age- and sex-matched controls in this study. The static CS in the glaucoma group was significantly reduced in 9 out of the 13 tested locations in the central 10° of the visual field. The mean static CS at 5°, 10°, superior hemifield and inferior hemifield were all significantly reduced in the glaucoma patients compared to the controls.

Conclusion

Static CS measurement is a sensitive approach that can be utilized to aid in the detection of glaucoma. The use of static CS can be adopted in the development of a cost-effective yet sensitive screening tool for the detection of glaucoma.

The Open Perimetry Initiative: A framework for cross-platform development for the new generation of portable perimeters

The Open Perimetry Initiative was formed in 2010 with the aim of reducing barriers to clinical research with visual fields and perimetry. Our two principal tools are the Open Perimetry Interface (OPI) and the visualFields package with analytical tools. Both are fully open source. The OPI package contains a growing number of drivers for commercially available perimeters, head-mounted devices, and virtual reality headsets. The visualFields package contains tools for the analysis and visualization of visual field data, including methods to compute deviation values and probability maps. We introduce a new frontend, the opiApp, that provides tools for customization for visual field testing and can be used as a frontend to run the OPI. The app can be used on the Octopus 900 (Haag-Streit), the Compass (iCare), the AP 7000 (Kowa), and the IMO (CREWT) perimeters, with permission from the device manufacturers. The app can also be used on Android phones with virtual reality headsets via a new driver interface, the PhoneHMD, implemented on the OPI. The use of the tools provided by the OPI library is showcased with a custom static automated perimetry test for the full visual field (up to 50 degrees nasally and 80 degrees temporally) developed with the OPI driver for the Octopus 900 and using visualFields for statistical analysis. With more than 60 citations in clinical and translational science journals, this initiative has contributed significantly to expand research in perimetry. The continued support of researchers, clinicians, and industry are key in transforming perimetry research into an open science.

Parasympathetic innervation of emmetropization

Emmetropization is affected by the temporal parameters of visual stimulation and the spectral composition of light, as well as by autonomic innervation. The goal of the current experiments is to test the hypothesis that different types of visual stimulation interact with ocular innervation in the process of emmetropization. For that, selective lesions of the autonomic nervous system were performed in chickens: involving transection of parasympathetic input to the eye from either the ciliary ganglion, innervating accommodation and pupil responses (CGX; n = 32), or pterygopalatine ganglion, innervating choroidal blood vessels and cornea (PPGX; n = 26). After 1 week of recovery, chicks were exposed to sinusoidally modulated light (3 days, 2 Hz, 680 lux) that was either achromatic (black to white [RGB], or black to yellow [RG]), or chromatic (blue to yellow [B/Y] or red to green [R/G]). Exposure to light stimulation was followed by ocular biometry (Lenstar and a Hartinger refractometer). Surgical conditions revealed a small reduction in anterior chamber depth with CGX but no other significant changes in ocular biometry/refraction under standard light conditions. With RGB achromatic stimulation, CGX eyes produced an effect on ocular components, with a further reduction in anterior chamber depth and an increase in vitreous chamber depth, while RG stimulation showed no effect. No effect was detected in PPGX under both achromatic protocols. With chromatic stimulation, CGX with R/G modulation increased eye length, while PPGX with B/Y modulation decreased eye length. We conclude that the two different types of parasympathetic innervations have antagonistic effects on eye growth and the anterior eye when challenged with the appropriate stimulus, with possible implications for the role of choroidal blood flow in emmetropization.

Interpreting Retinal Nerve Fiber Layer Reflectance Defects Based on Presence of Retinal Nerve Fiber Bundles

SIGNIFICANCE

Adaptive-optics scanning-laser-ophthalmoscopy (AOSLO) retinal imaging of the retinal nerve fiber layer (RNFL) helps predict the severity of perimetric damage based on absence of fibers and projection of the defects in en face images of the RNFL from spectral-domain optical coherence tomography (SD-OCT).

PURPOSE

En face images of the RNFL reveal reflectance defects in patients with glaucoma and predict locations of perimetric defects. These defects could arise from either loss of retinal nerve fiber bundles or reduced bundle reflectance. This study used AOSLO to assess presence of bundles in areas with RNFL reflectance defects on SD-OCT.

METHODS

Adaptive-optics scanning laser ophthalmoscopy was used to image a vertical strip of RNFL measuring approximately 30 × 3° between the optic disc and the fovea. Fifteen patients with glaucoma who had SD-OCT reflectance defects that passed through this region were chosen. Four patients had reflectance defects in both superior and inferior hemifields, so presence of bundles on AOSLO was assessed for 19 hemifields. Where bundles were present, the hemifield was scored for whether bundles seemed unusual (low contrast and/or low density). Perimetric defects were considered deep when sensitivity was below 15 dB.

RESULTS

Ten hemifields had a region with no fibers present on AOSLO; all had a corresponding deep perimetric defect. The other nine hemifields had no region in the AOSLO image without fibers: four with normal fibers and five with unusual fibers. The only one of these nine hemifields with a deep perimetric defect was one with low-contrast fibers and overall thin RNFL.

CONCLUSIONS

Retinal nerve fiber layer reflectance defects, which were associated with deep perimetric defects, usually had a region with absence of fibers on AOSLO images of RNFL. Ability to predict severity of perimetric damage from en face SD-OCT RNFL reflectance images could benefit from quantification that differentiated between absence of fibers and unusual fibers.

A Novel Stimulus to Improve Perimetric Sampling within the Macula in Patients with Glaucoma

SIGNIFICANCE

Identifying glaucomatous damage to the macula has become important for diagnosing and managing patients with glaucoma. In this study, we presented an approach that provides better perimetric sampling for the macular region, by testing four locations, with a good structure-function agreement.

PURPOSE

We previously presented a basis for customizing perimetric locations within the macula. In this study, we aimed to improve perimetric sampling within the macula by presenting a stimulus at four locations, with maintaining a good structure-function agreement.

METHODS

We tested one eye each of 30 patients (aged 50 to 88 years). Patients were selected based on observed structural damage to the macula, whereas perimetric defect (using 24-2) did not reflect the locations and extent of this damage. We used en face images to visualize retinal nerve fiber bundle defects. To measure perimetric sensitivities, we used a blob stimulus (standard deviation of 0.25°) at the 10-2 locations. A perimetric defect for a location was defined as any value equal to or deeper than -4, -5, and -6 dB below the mean sensitivity for 37 age-similar controls (aged 47 to 78 years). We also presented an elongated sinusoidal stimulus for 20 patients at four locations within the macula, in which we defined a perimetric defect as any value below the 2.5th percentile from controls.

RESULTS

The -4, -5, and -6 dB criteria identified perimetric defects in 14, 13, and 11 patients, respectively. When testing with the elongated stimulus, 18 patients were identified with perimetric defect. The perimetric defects were consistent with the structural damage.

CONCLUSIONS

The elongated stimulus showed a good structure-function agreement with only four testing locations as compared with 68 locations used with the blob stimulus. This demonstrates a clinical potential for this new stimulus in the next generation of perimetry.

Mapping the Contrast Sensitivity of the Visual Field With Bayesian Adaptive qVFM

Current clinical evaluation, which focuses on central vision, could be improved through characterization of residual vision with peripheral testing of visual acuity, contrast sensitivity, color vision, crowding, and reading speed. Assessing visual functions in addition to light sensitivity, a comprehensive visual field map (VFM) would be valuable for detecting and managing eye diseases. In a previous study, we developed a Bayesian adaptive qVFM method that combines a global module for preliminary assessment of the VFM's shape and a local module for assessment at individual retinal locations. The method was validated in measuring the light sensitivity VFM. In this study, we extended the qVFM method to measure contrast sensitivity across the visual field. In both simulations and psychophysics, we sampled 64 visual field locations (48 x 48 deg) and compared the qVFM method with a procedure that tested each retinal location independently (qFC; Lesmes et al., 2015). In each trial, subjects were required to identify a single optotype (size: 2.5 x 2.5 deg), one of 10 filtered Sloan letters. To compare the accuracy and precision of the two methods, three simulated eyes were tested in 1,280 trials with each method. In addition, data were collected from 10 eyes (5 OS, 5 OD) of five normal observers. For simulations, the average RMSE of the estimated contrast sensitivity with the qVFM and qFC methods were 0.057 and 0.100 after 320 trials, and 0.037 and 0.041 after 1,280 trials [all in log10 units, represent as log(sensitivity)], respectively. The average SD of the qVFM and qFC estimates were 0.054 and 0.096 after 320 trials, and 0.032 and 0.041 after 1,280 trials, respectively. The within-run variability (68.2% HWCIs) were comparable to the cross-run variability (SD). In the psychophysics experiment, the average HWCI of the estimated contrast sensitivity from the qVFM and qFC methods across the visual field decreased from 0.33 on the first trial to 0.072 and 0.16 after 160, and to 0.060 and 0.10 after 320 trials. The RMSE between the qVFM and qFC estimates started at 0.26, decreased to 0.12 after 160 and to 0.11 after 320 qVFM trials. The qVFM provides an accurate, precise, and efficient mapping of contrast sensitivity across the entire visual field. The method might find potential clinical applications in monitoring vision loss, evaluating therapeutic interventions, and developing effective rehabilitation for visual diseases.

A novel Bayesian adaptive method for mapping the visual field

Measuring visual functions such as light and contrast sensitivity, visual acuity, reading speed, and crowding across retinal locations provides visual-field maps (VFMs) that are extremely valuable for detecting and managing eye diseases. Although mapping light sensitivity is a standard glaucoma test, the measurement is often noisy (Keltner et al., 2000). Mapping other visual functions is even more challenging. To improve the precision of light-sensitivity mapping and enable other VFM assessments, we developed a novel hybrid Bayesian adaptive testing framework, the qVFM method. The method combines a global module for preliminary assessment of the VFM's shape and a local module for assessing individual visual-field locations. This study validates the qVFM method in measuring light sensitivity across the visual field. In both simulation and psychophysics studies, we sampled 100 visual-field locations (60° × 60°) and compared the performance of qVFM with the qYN procedure (Lesmes et al., 2015) that measured light sensitivity at each location independently. In the simulations, a simulated observer was tested monocularly for 1,000 runs with 1,200 trials/run, to compare the accuracy and precision of the two methods. In the experiments, data were collected from 12 eyes (six left, six right) of six human subjects. Subjects were cued to report the presence or absence of a target stimulus, with the luminance and location of the target adaptively selected in each trial. Both simulations and a psychological experiment showed that the qVFM method can provide accurate, precise, and efficient mapping of light sensitivity. This method can be extended to map other visual functions, with potential clinical signals for monitoring vision loss, evaluating therapeutic interventions, and developing effective rehabilitation for low vision.

Insights for mfVEPs from perimetry using large spatial frequency-doubling and near frequency-doubling stimuli in glaucoma

PURPOSE

To compare two forms of perimetry that use large contrast-modulated grating stimuli in terms of: their relative diagnostic power, their independent diagnostic information about glaucoma and their utility for mfVEPs. We evaluated a contrast-threshold mfVEP in normal controls using the same stimuli as one of the tests.

METHODS

We measured psychophysical contrast thresholds in one eye of 16 control subjects and 19 patients aged 67.8 ± 5.65 and 71.9 ± 7.15, respectively, (mean ± SD). Patients ranged in disease severity from suspects to severe glaucoma. We used the 17-region FDT-perimeter C20-threshold program and a custom 9-region test (R9) with similar visual field coverage. The R9 stimuli scaled their spatial frequencies with eccentricity and were modulated at lower temporal frequencies than C20 and thus did not display a clear spatial frequency-doubling (FD) appearance. Based on the overlapping areas of the stimuli, we transformed the C20 results to 9 measures for direct comparison with R9. We also compared mfVEP-based and psychophysical contrast thresholds in 26 younger (26.6 ± 7.3 y, mean ± SD) and 20 older normal control subjects (66.5 ± 7.3 y) control subjects using the R9 stimuli.

RESULTS

The best intraclass correlations between R9/C20 thresholds were for the central and outer regions: 0.82 ± 0.05 (mean ± SD, p ≤ 0.0001). The areas under receiver operator characteristic plots for C20 and R9 were as high as 0.99 ± 0.012 (mean ± SE). Canonical correlation analysis (CCA) showed significant correlation (r = 0.638, p = 0.029) with 1 dimension of the C20 and R9 data, suggesting that the lower and higher temporal frequency tests probed the same neural mechanism(s). Low signal quality made the contrast-threshold mfVEPs non-viable. The resulting mfVEP thresholds were limited by noise to artificially high contrasts, which unlike the psychophysical versions, were not correlated with age.

CONCLUSION

The lower temporal frequency R9 stimuli had similar diagnostic power to the FDT-C20 stimuli. CCA indicated the both stimuli drove similar neural mechanisms, possibly suggesting no advantage of FD stimuli for mfVEPs. Given that the contrast-threshold mfVEPs were non-viable, we used the present and published results to make recommendations for future mfVEP tests.

Using Small Samples to Evaluate Normative Reference Ranges for Retinal Imaging Measures

SIGNIFICANCE

Retinal nerve fiber layer (RNFL) deviation maps often incorrectly score healthy eyes as having wedge defects. This study shows how to identify such problems early in the development of normative databases.

PURPOSE

After reference values are embedded in devices, clinicians and researchers often learn about issues that cause false-positive rates in healthy eyes. Here we show a way to detect and address such issues early on.

METHODS

The thickness of the RNFL was measured for both eyes of 60 healthy younger adults aged 20 to 31 years and one eye each of 30 healthy older adults aged 54 to 82 years. Deviation maps were developed from the left eyes of the first 30 younger adults, and between-subject variability in the shape of the RNFL was assessed. This was repeated in their right eyes, in the second group of younger adults and in the older adults.

RESULTS

For the first group of 30 healthy young adults, between-subject variability in the location of the region of greatest thickness meant that 58% of the pixels below the fifth percentile in the left eyes were from four people whose deviation maps had wedge-shaped patterns, as did the deviation maps for the nine right eyes with 87% of the pixels below the fifth percentile. Wedge patterns were also seen in deviation maps for 8 left eyes and 11 right eyes of the second group of young adults and for 9 eyes of the older adults.

CONCLUSIONS

Evaluation of RNFL thickness maps from 30 young adults was sufficient to determine that between-subject variability in the shape of the RNFL can cause wedge patterns in RNFL deviation maps in many healthy eyes.

Within-subject variability in human retinal nerve fiber bundle width

With the growing availability of high-resolution imaging there has been increased interest in developing new metrics for integrity of the retinal nerve fiber layer. In particular, it has been suggested that measurement of width of retinal nerve fiber bundles (RNFBs) may be useful in glaucoma, due to low between-subject variability in mean RNFB width. However, there have also been reports of substantial within-subject variability in the width of individual RNFBs. To assess within-subject variability as a potential source of selection bias in measurements of RNFB width, we used an adaptive optics scanning laser ophthalmoscope (AOSLO) to measure widths of individual RNFBs in one eye each of 11 young adults in good ocular health. In a pilot study we analyzed a large AOSLO image of RNFL in one participant then, based on those findings, in the main study we used AOSLO to image a smaller region in 10 additional healthy young adults. The pilot study of one eye found RNFB widths ranging from 10 μm to 44 μm. This suggested that biological variability was too high for measuring small changes arising from disease processes. This was confirmed in measurements of 10 eyes in the main study, RNFB widths ranged from 9 μm to 55 μm and every eye had large within-subject variability (exceeding 19 μm in all eyes) in RNFB width for nearby bundles. The within-subject variability in RNFB width, as well as variation in the width of single RNFBs over relatively short distances (<300 um) depending on the precise location of measurement, suggests that bundle width measurements would be highly susceptible to selection bias and therefore of limited clinical use.

Novel Technique for Quantifying Retinal Nerve Fiber Bundle Abnormality in the Temporal Raphe

SIGNIFICANCE

Glaucomatous nasal visual field abnormalities correspond to damage in the temporal raphe-where individual nerve bundles can be visualized. The ability to quantify structural abnormality in the raphe, with a clinically applicable protocol, sets the stage for investigating the raphe as a potential site for assessing early glaucoma.

PURPOSE

To develop a clinically applicable imaging and analysis technique for identifying retinal nerve fiber bundle abnormalities in the temporal raphe.

METHODS

Spectralis optical coherence tomography scans customized for the temporal raphe were gathered from 30 younger controls, 30 older controls, and 29 patients with glaucoma. An analysis technique was developed based on the reflectance of the nerve fiber bundles. The technique was first developed in the younger controls, and then applied to the older controls to generate normative data for quantifying nerve fiber bundle reflectance abnormalities in the patients with glaucoma. Matrix perimetric data were gathered in the patients with glaucoma to evaluate the reflectance technique's findings. Reflectance abnormality in the patients was defined when the fraction of enface area showing reflectance abnormality was greater than the 95th percentile estimated from controls. Spearman's rho was used to quantify the relation between the total deviation at the perimetric testing locations and the fraction of corresponding enface area showing reflectance abnormality.

RESULTS

Twenty-five of the 29 patients had reflectance abnormalities. Eight of these had mild to no perimetric mean deviation abnormality. Similar results were found when perimetric total deviations were compared to reflectance abnormalities in the corresponding enface locations. Spearman's rho comparing the total deviations to reflectance abnormalities found rs(174) = -0.72, P < .001.

CONCLUSIONS

The technique typically identified reflectance abnormality when perimetric abnormality was present. It also identified reflectance abnormalities even when perimetric abnormality was mild or absent. The findings support the potential of raphe imaging in detecting early glaucomatous damage.

Comparison of defect depths for sinusoidal and circular perimetric stimuli in patients with glaucoma

Purpose

Clinical use of perimetric testing in patients with glaucoma typically assumes that perimetric defects will be less deep for larger than smaller stimuli. However, studies have shown that very large sinusoidal stimuli can yield similar defects as small circular stimuli. In order to provide guidelines for new perimetric stimuli, we tested patients with glaucoma using five different stimuli and compared defects to their patterns of retinal nerve fibre layer (RNFL) damage.

Methods

Twenty subjects with glaucoma were imaged with optical coherence tomography (OCT) volume scans to allow for en face RNFL images and were also tested on a custom perimetry station with five stimuli: Goldmann sizes III and V, a two‐dimensional Gaussian blob (standard deviation 0.5°) and a 0.5 cycle degree⁻¹ sinusoidal grating presented two ways: flickered at 5 Hz, and pulsed for 200 ms instead of flickered. En face RNFL images were reviewed with the visual field locations overlaid, and each location was labelled for a patient as either no visible RNFL defect or as wedge, slit, edge, or diffuse defect. Nineteen age‐similar controls were tested with the same stimuli to define depth of defect as difference from mean normal. Bland‐Altman analysis was used to test three predictions of neural modelling by making five comparisons.

Results

Bland‐Altman analysis confirmed the three predictions. The flickered sinusoid gave deeper defects in damaged areas than the pulsed sinusoid (r = 0.25, p < 0.0001). When comparing data for sizes III and V there was increased spread of the data in deeper defects in the direction of size III having deeper defect (r = 0.35, p < 0.0001). The size V stimulus yielded shallower defects than a stimulus of similar size but with blurred edges (r = 0.20, p = 0.0004).

Conclusions

On average, all stimuli produced similar results comparing across type of RNFL damage. However, there were systematic patterns consistent with predictions of neural modelling: in damaged areas, depth of defect tended to be greater for the flickered sinusoid than the pulsed sinusoid, greater for the size III stimulus than the size V stimulus, and greater for the Gaussian blob than for the size V stimulus.

Identifying Glaucomatous Damage to the Macula

SIGNIFICANCE

Measurements of the macula have been increasingly used to diagnose and manage patients with glaucoma. Asymmetry analysis was clinically introduced to assess damage to the macular ganglion cells in patients with glaucoma, but its effectiveness is limited by high normal between-subject variability.

PURPOSE

We aimed to reduce the high normal between-subject variability and improve the potential of asymmetry analysis to identify glaucomatous damage to the macula.

METHODS

Twenty patients with glaucoma (aged 57 to 85 years) and 30 age-similar control subjects (aged 53 to 89 years) were recruited from a longitudinal glaucoma study. Participants were imaged with the Spectralis OCT using the posterior pole protocol; measurements of the averaged retinal thickness and ganglion cell layer (GCL) thickness were obtained. We established three zones per hemifield within the central ±9°, based on the lowest between-subject variability that we previously found and the course of retinal nerve fiber layer projections. The criteria for flagging abnormality were at least two contiguous zones when P < 5% or one zone when P < 1% with two-tailed tests.

RESULTS

Between-subject variability of the asymmetry analysis for both retinal and GCL thicknesses remained lower than that of the average thickness across each zone in control subjects (F > 2.52, P < .01). Asymmetry analysis of retinal and GCL thicknesses flagged 16 and 18 of 20 patients, respectively.

CONCLUSIONS

Between-subject variability was reduced in control subjects using the three zones; our criteria identified glaucomatous damage to the macula in most of the patients. We used high-density B-scans to confirm the patterns of the glaucomatous damage we found in this study.

Prediction Accuracy of the Dynamic Structure-Function Model for Glaucoma Progression Using Contrast Sensitivity Perimetry and Confocal Scanning Laser Ophthalmoscopy

PURPOSE

The purpose of this study was to determine whether combining a structural measure with contrast sensitivity perimetry (CSP), which has lower test-retest variability than static automated perimetry (SAP), reduces prediction error with 2 models of glaucoma progression.

METHODS

In this retrospective analysis, eyes with 5 visits with rim area (RA), SAP, and CSP measures were selected from 2 datasets. Twenty-six eyes with open-angle glaucoma were included in the analyses. For CSP and SAP, mean sensitivity (MS) was obtained by converting the sensitivity values at each location from decibel (SAP) or log units (CSP) to linear units, and then averaging all values. MS and RA values were expressed as percent of mean normal based on independent normative data. Data from the first 3 and 4 visits were used to calculate errors in prediction for the fourth and fifth visits, respectively. Prediction errors were obtained for simple linear regression and the dynamic structure-function (DSF) model.

RESULTS

With linear regression, the median prediction errors ranged from 6% to 17% when SAP MS and RA were used and from 9% to 17% when CSP MS and RA were used. With the DSF model, the median prediction errors ranged from 6% to 11% when SAP MS and RA were used and from 7% to 16% when CSP MS and RA were used.

CONCLUSIONS

The DSF model had consistently lower prediction errors than simple linear regression. The lower test-retest variability of CSP in glaucomatous defects did not, however, result in lower prediction error.

Associations between structure and function are different in healthy and glaucomatous eyes

Purpose

To assess if there are differences in the structure-function associations between healthy and glaucomatous eyes.

Methods

Structure-function associations were assessed in healthy and glaucomatous eyes in three datasets, globally and in the six sectors of the optic nerve head. Structural parameters included rim area (RA) and retinal nerve fiber layer thickness (RNFLT). Functional parameters included unweighted mean of sensitivity thresholds (MS) and unweighted mean of total deviation values (MD), assessed with standard automated perimetry, short-wavelength automated perimetry, frequency-doubling technology perimetry, or contrast sensitivity perimetry. All structural and functional parameters were expressed as percent of mean normal. SF associations were assessed with correlation analyses (Pearson, Spearman and Kendall). We also assessed the SF associations with linear regression analyses: the generalized estimating equation (GEE) was used to adjust for inter-eye correlations and ordinary least squares (OLS) linear models were used when these adjustments were not necessary. We applied Bonferroni corrections to adjust for the impact of multiple comparisons.

Results

Overall, none of the Pearson correlations tested in healthy eyes were significant (correlations ranged from -0.17 to 0.37), whereas 77% of the correlations tested in glaucomatous eyes were significant (correlations ranged from 0.01 to 0.79). Similarly, none of the slopes obtained with GEE and OLS were significant in healthy eyes (slopes ranged from -0.30 to 0.87), whereas 82% of the slopes obtained in glaucomatous eyes were significant (slopes ranged from 0.02 to 1.38).

Conclusions

Significant associations between structure and function were consistently observed in glaucomatous eyes, but not in healthy eyes. These differences in association should be considered in the design of structure-function models for progression.

Evaluating glaucomatous abnormality in peripapillary optical coherence tomography enface visualisation of the retinal nerve fibre layer reflectance

Purpose

Optical coherence tomography (OCT) enface visualisation of the retinal nerve fibre layer (RNFL) reflectance has been found to have some advantages over retinal thickness measures. However, it is not yet clear how abnormalities on enface images relate to findings of abnormalities from other clinical measures such as the circumpapillary retinal nerve fibre layer thickness (cRNFLT). We developed a technique to analyse the RNFL reflectance on the OCT enface images, and to investigate its relation with the cRNFLT.

Methods

Spectralis (http://www.heidelbergengineering.com) OCT scans of the central retinal ±24° were analysed in the study eye of 31 controls and 33 patients, ages 61 (±9) and 69 (±8) years respectively. Enface slab‐images were extracted at 16–24, 24–36, and 24–52 μm from the inner limiting membrane in the temporal raphe, perifoveal and disc regions respectively. Reflectance probability maps were generated for the patients based on the control data. Glaucomatous abnormality was defined on the slab‐images when the slab‐area with reflectance abnormality was greater than the 95th percentile, and on the cRNFLT when the thickness measure was less than the fifth percentile, of that found in controls. The fraction of slab‐image showing reflectance abnormality was compared to cRNFLT in the patient group, using Spearman's rho. Agreement between the findings of abnormality based on cRNFLT and slab‐image reflectance was assessed using Cohen's kappa.

Results

Slab‐image and cRNFLT findings were in agreement for 26/33 eyes; four subjects showed cRNFLT abnormality but not slab‐image abnormality, and three subjects showed slab‐image abnormality but not cRNFLT abnormality. Spearman's rho found r_s(31) = −0.82. The reflectance findings and cRNFLT findings were consistent in 27/33 for both the superior temporal (ST) and inferior temporal (IT) sectors, and Cohen's kappa found 0.53 and 0.61 respectively.

Conclusion

The surface area of enface slab‐images showing RNFL reflectance were strongly related to the cRNFLT measures, and the classification of a subject with glaucoma based on enface reflectance findings and cRNFLT findings had a generally good agreement. The larger retinal area assessed by the enface method preserves the spatial location of the RNFL abnormalities, and makes the technique a useful approach for identifying regions of potential RNFL abnormality for targeted perimetry.

Customizing Perimetric Locations Based on En Face Images of Retinal Nerve Fiber Bundles With Glaucomatous Damage

Purpose

Prior studies suggested the use of customized perimetric locations in glaucoma; these studies were limited by imaging only the superficial depths of the retinal nerve fiber layer (RNFL) and by prolonged perimetric testing. We aimed to develop a rapid perimetric test guided by high-resolution images of RNFL bundles.

Methods

We recruited 10 patients with glaucoma, ages 56 to 80 years, median 68 years, and 10 controls, ages 55 to 77 years, median 68 years. The patients were selected based on discrepancies between locations of glaucomatous damage for perimetric and structural measures. Montaging was used to produce optical coherence tomography en face images of the RNFL covering much of the 24-2 grid locations. In experiment 1, we presented the Goldmann size III stimulus at preselected retinal locations of glaucomatous damage, using just two contrasts. In experiment 2, we developed an elongated sinusoidal stimulus, aligned within the defect, to measure contrast sensitivities; abnormalities were defined based on lower 95% reference limits derived from the controls.

Results

The percentage of predicted locations where size III was not seen at 28 dB ranged from 16% to 80%, with a median of 48%. Contrast sensitivity for the sinusoidal stimulus was below the 95% reference range for 37 of 44 stimuli aligned within the defects.

Conclusions

We developed methods for rapid perimetric testing guided by en face images of the RNFL bundles in patients with glaucoma. Results indicated ganglion cell damage under all of the visible RNFL defects.

Translational Relevance

Customized perimetric locations have potential to improve clinical assessment of glaucoma.

A basis for customising perimetric locations within the macula in glaucoma

PURPOSE

It has been recognised that the 24-2 grid used for perimetry may poorly sample the macula, which has been recently identified as a critical region for diagnosing and managing patients with glaucoma. We compared data derived from patients and controls to investigate the efficacy of a basis for customising perimetric locations within the macula, guided by en face images of retinal nerve fibre layer (RNFL) bundles.

METHODS

We used SD-OCT en face montages (www.heidelbergengineering.com) of the RNFL in 10 patients with glaucoma (ages 56-80 years, median 67.5 years) and 30 age-similar controls (ages 47-77, median 58). These patients were selected because of either the absence of perimetric defect while glaucomatous damage to the RNFL bundles was observed, or because of perimetric defect that did not reflect the extent and locations of the glaucomatous damage that appeared in the RNFL images. We used a customised blob stimulus for perimetric testing (a Gaussian blob with 0.25° standard deviation) at 10-2 grid locations, to assess the correspondence between perimetric defects and damaged RNFL bundles observed on en face images and perimetric defects. Data from the age-similar controls were used to compute total deviation (TD) and pattern deviation (PD) values at each location; a perimetric defect for a location was defined as a TD or PD value of -0.5 log unit or deeper. A McNemar's test was used to compare the proportions of locations with perimetric defects that fell outside the damaged RNFL bundles, with and without accounting for displacement of ganglion cell bodies.

RESULTS

All patients but one had perimetric defects that were consistent with the patterns of damaged RNFL bundles observed on the en face images. We found six abnormal perimetric locations of 2040 tested in controls and 132 abnormal perimetric locations of 680 tested in patients. The proportions of abnormal locations that fell outside the damaged RNFL bundles, with and without accounting for displacement of the ganglion cell bodies were 0.08 and 0.07, respectively. The difference between the two proportions did not reach statistical significance (p = 0.5 for a one-tailed test).

CONCLUSIONS

We demonstrated that it is effective to customise perimetric locations within the macula, guided by en face images of the RNFL bundles. The perimetric losses found with a 10-2 grid demonstrated similar patterns as the damaged RNFL bundles observed on the en face images.

Retinal putative glial alterations: implication for glaucoma care

PURPOSE

Gliosis-like retinal alterations, presumed to be activated retinal astrocytes and Müller cells (ARAM), have been reported to occur frequently in patients with glaucoma but rarely in controls. We investigated the association between glaucomatous abnormality and the presence, the extent of retinal region, and the spatial distribution, of hyperreflective retinal alterations on optical coherence tomography (OCT) en-face images, presumed to be ARAM.

METHODS

Findings of hyperreflective structures, presumed to be ARAM, in the central retinal ±24 degrees of OCT en-face images (acquired with the SPECTRALIS^® OCT) were compared between 35 younger controls, 42 older controls and 38 patients with glaucoma. Presumed ARAM was defined as reflective structures on the en-face images other than retinal vasculature and retinal nerve fibre bundles. Chi-square tests were used to compare the proportion of younger controls vs older controls with presumed ARAM to investigate the effect of ageing, and the proportion of patients vs age-similar older controls with presumed ARAM to investigate the effect of disease. We also investigated the effect of glaucoma on the retinal area with presumed ARAM when it was present; we used an analysis of covariance (ancova) to compare the retinal area with hyperreflectivity in patients vs controls, adjusting for the effects of age and axial length.

RESULTS

The mean (S.D.) age of the younger controls, older controls, and patients with glaucoma was 26 (3), 62 (10) and 69 (8) years, respectively. The median (25th quartile, 75th quartile) of the retinal region with the hyperreflective structures, presumed to be ARAM, was zero (0,0), 1 (0,6), and 11 (0,43) degrees square in the younger controls, older controls and patients with glaucoma respectively. The chi-square test investigating the effect of ageing found χ² (1, N = 77) = 24.8, p < 0.001, and that investigating the effect of disease found χ² (1, N = 80) = 2.3, p = 0.1. The ancova found F(1, 46) = 10.32, P = 0.02.

CONCLUSIONS

There was an effect of ageing on the presence of the hyperreflective structures, presumed to be ARAM, on OCT images. Compared to the presence of hyperreflective structures, the extent of retinal region with the hyperreflective structures has a greater potential of being an indicator of glaucomatous degeneration. Further study is needed to investigate the nature of the relation between glaucomatous abnormality and the extent of the retina with the hyperreflective structures, presumed to be ARAM.

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.

Contrast sensitivity perimetry data from adults free of eye disease

This data article contains data referenced in “Individual Differences in the Shape of the Nasal Visual Field” [1]. The data were gathered from volunteers free of eye disease ages 21–85 who were tested with Contrast Sensitivity Perimetry (CSP), which uses a stimulus resistant to effects of defocus and reduced retinal illumination. Some subjects were tested only once or a few times, and others were part of a longitudinal cohort with as many as 10 tests. Parameters from maximum likelihood estimation of psychophysical threshold at each tested location are included in the data file, along with the participant׳s sex, age at time of test, the center of their physiological blind spot, the duration of test, the time of day that the test was begun, and the starting contrast used for the psychophysical staircases.

Individual differences in the shape of the nasal visual field

Between-subject differences in the shape of the nasal visual field were assessed for 103 volunteers 21-85years of age and free of visual disorder. Perimetry was conducted with a stimulus for which contrast sensitivity is minimally affected by peripheral defocus and decreased retinal illumination. One eye each was tested for 103 volunteers free of eye disease in a multi-center prospective longitudinal study. A peripheral deviation index was computed as the difference in log contrast sensitivity at outer (25-29° nasal) and inner (8° from fixation) locations. Values for this index ranged from 0.01 (outer sensitivity slightly greater than inner sensitivity) to -0.7 log unit (outer sensitivity much lower than inner sensitivity). Mean sensitivity for the inner locations was independent of the deviation index (R²<1%), while mean sensitivity for the outer locations was not (R²=38%, p<0.0005). Age was only modestly related to the index, with a decline by 0.017 log unit per decade (R²=10%). Test-retest data for 21 volunteers who completed 7-10 visits yielded standard deviations for the index from 0.04 to 0.17 log unit, with a mean of 0.09 log unit. Between-subject differences in peripheral deviation persisted over two years of longitudinal testing. Peripheral deviation indices were correlated with indices for three other perimetric stimuli used in a subset of 24 volunteers (R² from 20% to 49%). Between-subject variability in shape of the visual field raises concerns about current clinical visual field indices, and further studies are needed to develop improved indices.

Between-Subject Variability in Healthy Eyes as a Primary Source of Structural-Functional Discordance in Patients With Glaucoma

Purpose

To test with an independent data set the finding that between-subject variability in healthy eyes is the primary source of structural–functional discordance in patients with glaucoma.

Methods

Neuroretinal rim area, retinal nerve fiber layer thickness, and perimetric data were analyzed for one eye in each of 55 control subjects and for 245 right eyes of patients in the United Kingdom Glaucoma Treatment Study. Data were gathered with the Heidelberg Retina Tomograph (HRT), Stratus Optical Coherence Tomograph (OCT), and Humphrey Field Analyzer (HFA). Discordance was quantified as width of the limits of agreement from a Bland-Altman analysis of depth of defect. The ratio of variances (F test) for the patient and control groups was computed for comparisons of HFA-OCT, HFA-HRT, and OCT-HRT. Bonferroni adjustment required P less than 0.017 for statistical significance. The discordance in the patients was also quantified as the 95% prediction interval computed from the discordance in controls using the Hood-Kardon model for the HFA-OCT comparison.

Results

The F ratio comparing discordance in patients and controls was 0.77, 1.43, and 1.32 for the HFA-OCT, HFA-HRT, and OCT-HRT comparisons with P values 0.88, 0.06, and 0.11, respectively. For the Hood-Kardon model, 4.7% of the patients had discordance outside the 95% prediction interval computed from the discordance in controls. Similar results were obtained when all comparisons were repeated for left eyes of patients.

Conclusions

These results confirm previous findings that between-subject variability in healthy eyes is the primary source of structural–functional discordance in patients with glaucoma, and extends this finding to a structural–structural comparison.

Imaging Glaucomatous Damage Across the Temporal Raphe

PURPOSE

To image and analyze anatomical differences at the temporal raphe between normal and glaucomatous eyes using adaptive optics scanning laser ophthalmoscopy (AOSLO) and optical coherence tomography (OCT), and to relate these differences to visual field measurements.

METHODS

Nine glaucomatous eyes of 9 patients (age 54-78 years, mean deviation of visual field [MD] -5.03 to -0.20 dB) and 10 normal eyes of 10 controls (age 54-81, MD -1.13 to +1.39 dB) were enrolled. All the participants were imaged in a region that was centered approximately 9° temporal to the fovea. The size of imaging region was at least 10° vertically by 4° horizontally. The raphe gap, defined as the distance between the superior and inferior retinal nerve fiber layer (RNFL) bundles, was measured. A bundle index was computed to quantify the relative reflectivity and density of the nerve fiber bundles. We also measured thickness of the ganglion cell complex (GCC) and RNFL.

RESULTS

The raphe gap was larger in glaucomatous eyes than control eyes. Specifically, eight glaucomatous eyes with local averaged field loss no worse than -3.5 dB had larger raphe gaps than all control eyes. The bundle index, GCC thickness, and RNFL thickness were on average reduced in glaucomatous eyes, with the first two showing statistically significant differences between the two groups.

CONCLUSIONS

Structural changes in the temporal raphe were observed and quantified even when local functional loss was mild. These techniques open the possibility of using the raphe as a site for glaucoma research and clinical assessment.

Assessing assumptions of a combined structure-function index

Purpose

Medeiros et al. developed a combined structure-function index for glaucoma by combining two ganglion cell models developed by Harwerth et al. The current study assessed assumptions of the Medeiros combined structure-function index by evaluating whether the two Harwerth models gave similar distributions of ganglion cells in an independent dataset.

Methods

The Harwerth models were applied to our previously published data for retinal nerve fibre layer (RNFL) thickness (Stratus OCT 3.4) and visual field sensitivities (24-2 SITA Standard) from one eye each of 51 patients with glaucoma and 62 age-similar control subjects free of eye disease. RNFL thicknesses and perimetric sensitivities were converted to ganglion cell numbers using the Harwerth model for perimetry and the Harwerth model for RNFL. These two estimates of ganglion cell number were compared for the inferior temporal (IT) and superior temporal (ST) sectors of the optic disc and the corresponding visual field locations. Comparisons were made with 14 visual field locations per sector (including a point in the macula for the inferior temporal sector) and with 13 locations (no point in the macula). Data for controls and patients were analysed separately, comparing mean values for RNFL perimetry models. Bonferroni correction was applied to control for repeated tests of significance. The difference between mean values for the RNFL and perimetry models was quantified by equating the means for controls through reduction of the assumed axon diameter used by the RNFL model.

Results

For the control group, the Harwerth RNFL model yielded smaller mean number of retinal ganglion cells than the Harwerth perimetry model, 23–47% lower (t > 13, p < 0.0001). This corresponded to mean axon diameters from 0.48 to 0.69 μm, with the smallest axons when the 14th location in the macula was included. With these new axon diameters, estimates of ganglion cell numbers for patients were still lower than for the RNFL model, by 19–28% (t > 6.5, p < 0.0001).

Conclusions

The Harwerth RNFL model consistently gave lower ganglion cell numbers than the Harwerth perimetry model, and this discordance persisted in patients even after reducing assumed axon diameter for controls. This finding contradicts the assumptions of the Medeiros structure-function index.

Choice of Stimulus Range and Size Can Reduce Test-Retest Variability in Glaucomatous Visual Field Defects

PURPOSE

To develop guidelines for engineering perimetric stimuli to reduce test-retest variability in glaucomatous defects.

METHODS

Perimetric testing was performed on one eye for 62 patients with glaucoma and 41 age-similar controls on size III and frequency-doubling perimetry and three custom tests with Gaussian blob and Gabor sinusoid stimuli. Stimulus range was controlled by values for ceiling (maximum sensitivity) and floor (minimum sensitivity). Bland-Altman analysis was used to derive 95% limits of agreement on test and retest, and bootstrap analysis was used to test the hypotheses about peak variability.

RESULTS

Limits of agreement for the three custom stimuli were similar in width (0.72 to 0.79 log units) and peak variability (0.22 to 0.29 log units) for a stimulus range of 1.7 log units. The width of the limits of agreement for size III decreased from 1.78 to 1.37 to 0.99 log units for stimulus ranges of 3.9, 2.7, and 1.7 log units, respectively (F = 3.23, P < 0.001); peak variability was 0.99, 0.54, and 0.34 log units, respectively (P < 0.01). For a stimulus range of 1.3 log units, limits of agreement were narrowest with Gabor and widest with size III stimuli, and peak variability was lower (P < 0.01) with Gabor (0.18 log units) and frequency-doubling perimetry (0.24 log units) than with size III stimuli (0.38 log units).

CONCLUSIONS

Test-retest variability in glaucomatous visual field defects was substantially reduced by engineering the stimuli.

TRANSLATIONAL RELEVANCE

The guidelines should allow developers to choose from a wide range of stimuli.