Correlation between optic nerve head structural parameters and glaucomatous visual field indices

Abstract

We examined associations between optic nerve head structural parameters and glaucomatous visual field indices. The study population included patients with glaucomatous optic neuropathy who were evaluated at Aichi Medical University, Nagakute, Aichi, Japan, from October 2010 to January 2011. A total of 57 eyes from 33 patients were assessed. We measured visual field using a Humphrey field analyzer, peripapillary retinal nerve fiber layer thickness (RNFL-T), and Bruch’s membrane opening-minimum rim width (BMO-MRW) using spectral domain optical coherence tomography, and rim area with referring three-dimensional photography. Spearman’s rank correlation coefficients were calculated between the threshold of visual sensitivity or total deviation of visual field and the following five optic nerve head structural parameters: RNFL-T length, BMO-MRW length, rim area, and calculated RNFL-T and BMO-MRW volumes (each length multiplied by rim area). The Akaike information criterion was calculated to determine which structural parameter was the best predictor of each visual field index. Threshold of visual sensitivity had correlation coefficients of 0.23 with global sector of RNFL-T, 0.32 with BMO-MRW, 0.14 with rim area, 0.21 with RNFL-T volume, and 0.26 with BMO-MRW volume. The correlation coefficients for each parameter with total deviation of visual field were 0.22, 0.33, 0.28, 0.36, and 0.37, respectively. The Akaike information criterion of BMO-MRW showed the smallest values in analyses of both threshold of visual sensitivity and total deviation. The present results show that RNFL-T volume and BMO-MRW volume were more strongly correlated with total deviation than BMO-MRW, but BMO-MRW appeared to be the best predictor of the two glaucomatous visual field indices.

Keywords:

• Bruch’s membrane opening-minimum rim width
• glaucoma
• retinal nerve fiber layer thickness
• rim area
• spectral domain optical coherence tomography
• visual field

Introduction

Markers for glaucoma detection and risk profiling include several structural parameters within the optic nerve head (ONH). Recent use of spectral domain optical coherence tomography (SDOCT) to study the structure–function relationship in glaucoma reveals an association between loss of visual function in glaucoma and peripapillary retinal nerve fiber layer thickness (RNFL-T).^1–3 Furthermore, a case control study in early glaucoma shows high sensitivity of the newly defined Bruch’s membrane opening-minimum rim width (BMO-MRW): ie, the shortest distance from the edge of the retinal pigment epithelium–Bruch’s membrane complex to the inner limiting membrane in the ONH, determined using SDOCT.^4,5

Rim area in the ONH is another certain landmark for glaucoma detection and profiling. The accuracy of rim area measurement with SDOCT has been investigated,⁶ and data confirm the usefulness of optic disc rim area measured by SDOCT for early glaucoma detection.⁷ Morgan et al⁸ reported the accuracy of diagnosing glaucomatous optic disc damage based on neuroretinal rim area. Overall, these findings indicate that SDOCT is a promising technique, and that rim area is a structural parameter of absolute importance for glaucoma detection. The available data may also imply another type of important structural parameter in the ONH that seems to be closely associated with glaucomatous functional disorders – namely, RNFL-T volume and BMO-MRW volume.

In the present study we defined RNFL-T volume and BMO-MRW volume, which were calculated by multiplying rim area by RNFL-T and BMO-MRW, respectively. We next calculated Spearman’s rank correlations between the threshold of visual sensitivity or total deviation of visual field and the following five structural parameters in the ONH: RNFL-T and BMO-MRW lengths, rim area, and the calculated RNFL-T and BMO-MRW volumes. Finally, we calculated the Akaike information criterion (AIC)⁹ to determine which of the five ONH structural parameters was the best predictor of the threshold of visual sensitivity or total deviation of visual field among patients with glaucomatous optic neuropathy.

Methods

Patients

This study included patients with primary open-angle glaucoma who were seen at Aichi Medical University, Nagakute, Aichi, Japan, from October 2010 to January 2011. Exclusion criteria were history of any other ocular disease, previous ocular surgery, tilted disc, micropapillae, myopia lower than −10.0 diopter, hyperopia greater than +4.0 diopter, best-corrected visual acuity of less than 20/40, visual field test with more than 20% fixation losses, and false positive or false negative errors greater than 15%. The study protocol were approved by the internal review board of our institution, and all patients provided informed consent prior to study entry.

RNFL-T and BMO-MRW measurement by SDOCT

The SDOCT RS-3000 (Nidek Co., Ltd, Aichi, Japan) was used, along with the Nidek Advanced Vision Information System (NAVIS) EX version 1.3.0.3. (Nidek Co., Ltd). We used the optic disc radial 6 program (5.7 mm) and the optic disc map program (5.1×5.1 mm) as a scanning pattern. This version of the software automatically draws the contour line of the disc margin. RNFL-T was calculated automatically by the optic disc map program (5.1×5.1 mm) in NAVIS. We divided the ONH into six sectors as described by Garway-Heath et al¹⁰ and we calculated the mean, which was considered the RNFL-T for each sector.

The BMO-MRW was defined as the shortest distance from the edge of the retinal pigment epithelium–Bruch’s membrane complex to the inner limiting membrane at the rim.^4,5 We measured this distance from the captured image using the optic disc radial 6 program, 12 regions of the optic disc at 30° intervals, on the manual setting using a measurement tool in the NAVIS program. We divided the ONH into six sectors as reported by Garway-Heath et al¹⁰ and calculated the mean, which was considered the BMO-MRW for each sector.

Rim area

An ophthalmologist with over 20 years of experience evaluated the optic disc stereographs, which were obtained using a Kowa nonmyd WX 3D Retinal Camera^® (Kowa Company Ltd, Aichi, Japan) in stereo mode (angle of view: 34°, width: 20°, length: 27°). Based on the stereograph analysis, the ophthalmologist defined the optic disc and the cup borders on a scanning laser ophthalmoscope image from the optic disc map program (5.1×5.1 mm). The centroid of the optic disc (Figure 1A, green point) was determined on the scanning laser ophthalmoscope image using ImageJ 1.46r (Wayne Rasband National Institutes of Health, http://rsb.info.nih.gov/ij/). We added a layer to represent the boundary lines for six sectors using Adobe^© Photoshop^© CS5 Extended (Adobe Systems Incorporated, San Jose, CA, USA) (Figure 1A, red lines). We also added layers for the optic disc area (Figure 1B, blue oval) and the cup area (Figure 1C, yellow oval). These three layers were combined, and the resultant blue ring was defined as the rim area (Figure 1D). We defined 715 pixels as 5.1 mm because the image size of 715×715 pixels corresponded to the scan size of 5.1×5.1 mm on the optic disc map. Finally, each rim area, divided into six sectors, was calculated using Adobe Photoshop CS5 Extended (Adobe Systems Incorporated).

Figure 1 Measurement of the rim area of the right eye on images from a scanning laser ophthalmoscope.

Notes: (A) Centroid of the optic disc (green point) and boundary lines dividing the six sectors (red lines). (B) The optic disc area (blue oval). (C) The cup area (yellow oval). (D) A combination of the three layers (A–C), with the resultant blue ring defined as the rim area. Numbers indicate the six sectors: 1, temporal; 2, superotemporal; 3, superonasal; 4, nasal; 5, inferonasal; 6, inferotemporal.

Visual field analysis

White-on-white standard automated perimetry was performed using a Humphrey Field Analyzer^® Model 750 (Carl Zeiss Meditec, Jena, Germany) and the C-30-2 SITA-standard strategy. We used threshold decibel (dB) values and total deviation as glaucomatous indices for corresponding topographic locations of visual field. We analyzed the relationship between 52 data points of the right eye obtained from the threshold static perimeter and the six corresponding sectors of the optic disc. These points were separated from each other by 6° and from the horizontal and vertical meridians by 3°. A mirror image version was prepared using data obtained from examinations of the left eye. A geometric mean of threshold and a mean of total deviation were calculated for each individual sector.

Statistical analysis

All patients enrolled in the study were included in the statistical analysis. Spearman’s rank correlation coefficients were calculated to analyze the associations between the threshold of visual sensitivity or total deviation and the following five structural parameters: RNFL-T, BMO-MRW, rim area, RNFL-T volume, and BMO-MRW volume. Fisher’s transformation was applied to calculate P-values to test the hypotheses of zero correlation and differences between correlations. The threshold of visual sensitivity and total deviation for the global sector were analyzed using a linear mixed effect model, including each structural parameter as a fixed effect, and eye and sector as random effects, considering the repeated measurements on each patient. The AIC⁹ was calculated to determine which of the five structural parameters was the best predictor of the threshold of visual sensitivity or total deviation of visual field based on the aforementioned mixed effect model. Smaller AIC values indicated better models. A P-value of <0.05 was considered statistically significant. All statistical analyses were performed using SAS 9.3 (SAS Institute, Inc., Cary, NC, USA).

Results

The study included a total of 57 eyes from 33 patients (20 female and 13 male; mean age, 53.7±16.3 years; range, 20–78 years). Mean intraocular pressure measured with Goldmann applanation tonometry was 14.8±3.47 mmHg (range, 7–22 mmHg).

The threshold of visual sensitivity was significantly correlated with the global sectors of RNFL-T (r=0.23, P<0.001), BMO-MRW (r=0.32, P<0.001), rim area (r=0.14, P=0.007), RNFL-T volume (r=0.21, P<0.001), and BMO-MRW volume (r=0.26, P<0.001; Table 1). Pairwise comparison analysis showed that threshold of visual sensitivity was more significantly correlated with BMO-MRW than with RNFL-T, rim area, and RNFL-T volume. Threshold of visual sensitivity was also more significantly correlated with BMO-MRW volume than with rim area and RNFL-T volume.

Table 1 Associations between threshold of visual sensitivity and five structural parameters

Sector	Statistics	RNFL-T	BMO-MRW	Rim area	RNFL-T × rim area	BMO-MRW × rim area
1	ra	0.31	0.24	0.19	0.37	0.32
	P-valueb	0.019	0.076	0.163	0.005	0.014
2	ra	0.54	0.54	0.52	0.57	0.56
	P-valueb	<0.001	<0.001	<0.001	<0.001	<0.001
3	ra	0.33	0.38	0.40	0.41	0.42
	P-valueb	0.012	0.003	0.002	0.001	<0.001
4	ra	−0.01	0.15	0.08	0.07	0.12
	P-valueb	0.95	0.283	0.549	0.592	0.387
5	ra	0.25	0.37	−0.01	0.18	0.25
	P-valueb	0.058	0.004	0.914	0.189	0.056
6	ra	0.39	0.47	0.25	0.41	0.44
	P-valueb	0.003	<0.001	0.057	0.002	<0.001
Global (n=342)	ra	0.23	0.32	0.14	0.21	0.26
	P-valueb	<0.001	<0.001	0.007	<0.001	<0.001
P-value for pairwise comparison of correlation coefficientsc
	RNFL-T		0.030	0.302	0.694	0.587
	BMO-MRW			0.006	0.015	0.148
	Rim area				0.101	<0.001
	RNFL-T × rim area					0.044

Notes:

a Spearman’s rank correlation coefficient.

b Correlation coefficient from testing the null hypothesis of zero correlation.

c Pairwise comparison of the correlation coefficients between the structural parameters for the global sector. Numbers indicate the six sectors: 1, temporal; 2, superotemporal; 3, superonasal; 4, nasal; 5, inferonasal; 6, inferotemporal.

Abbreviations: BMO-MRW, Bruch’s membrane opening-minimum rim width; RNFL-T, retinal nerve fiber layer thickness.

The global sectors of RNFL-T (r=0.22, P<0.001), BMO-MRW (r=0.33, P<0.001), rim area (r=0.28, P<0.001), RNFL-T volume (r=0.36, P<0.001), and BMO-MRW volume (r=0.37, P<0.001) also demonstrated significant correlation coefficients with total deviation (Table 2). Pairwise comparison analysis showed that total deviation was more significantly correlated with BMO-MRW than with RNFL-T. The analysis also revealed that total deviation was more significantly correlated with both BMO-MRW volume and RNFL-T volume than with RNFL-T and rim area.

Table 2 Associations between total deviation and five structural parameters

Sector	Statistics	RNFL-T	BMO-MRW	Rim area	RNFL-T × rim area	BMO-MRW × rim area
1	ra	0.22	0.08	0.30	0.38	0.27
	P-valueb	0.095	0.543	0.021	0.003	0.039
2	ra	0.51	0.53	0.44	0.53	0.53
	P-valueb	<0.001	<0.001	<0.001	<0.001	<0.001
3	ra	0.36	0.39	0.45	0.46	0.46
	P-valueb	0.006	0.003	<0.001	<0.001	<0.001
4	ra	−0.01	0.19	0.07	0.07	0.14
	P-valueb	0.922	0.159	0.596	0.609	0.295
5	ra	0.23	0.29	−0.03	0.17	0.20
	P-valueb	0.079	0.027	0.817	0.213	0.133
6	ra	0.49	0.46	0.35	0.52	0.50
	P-valueb	<0.001	<0.001	0.007	<0.001	<0.001
Global (n=342)	ra	0.22	0.33	0.28	0.36	0.37
	P-valueb	<0.001	<0.001	<0.001	<0.001	<0.001
P-value for pairwise comparison of correlation coefficientsc
	RNFL-T		0.013	0.470	0.006	0.019
	BMO-MRW			0.391	0.494	0.333
	Rim area				0.014	0.001
	RNFL-T × rim area					0.771

Notes:

a Spearman’s rank correlation coefficient.

b Correlation coefficient from testing the null hypothesis of zero correlation.

c Pairwise comparison of the correlation coefficients between the structural parameters for the global sector. Numbers indicate the six sectors: 1, temporal; 2, superotemporal; 3, superonasal; 4, nasal; 5, inferonasal; 6, inferotemporal.

Abbreviations: BMO-MRW, Bruch’s membrane opening-minimum rim width; RNFL-T, retinal nerve fiber layer thickness.

Among the five structural parameters, the standardized regression coefficient divided by its standard error of BMO-MRW for the threshold of visual sensitivity was the largest (Table 3). AIC under the model based on BMO-MRW as a predictor of the threshold of visual sensitivity was smaller than that under the other structural parameters. Similar results were found in analyses of the total deviation of visual field.

Table 3 Relationships between threshold of visual sensitivity or total deviation and five structural parameters for the global sector

	RNFL-T	BMO-MRW	Rim area	RNFL-T × rim area	BMO-MRW × rim area
Threshold of visual sensitivity
Estimatea	4.66	6.34	1.17	2.73	1.46
P-valueb	<0.001	<0.001	0.245	0.007	0.146
AICc	1,940	1,924d	1,947	1,950	1,958
Total deviation
Estimatea	4.62	6.77	4.38	6.14	4.30
P-valueb	<0.001	<0.001	<0.001	<0.001	<0.001
AICc	2,076	2,054d	2,065	2,057	2,077

Notes:

a Standardized regression coefficient (estimated regression coefficient divided by its standard error) based on mixed effect model.

b Testing for the regression coefficient (null hypothesis for “regression coefficient=0”).

c Smaller values of AIC indicate better models.

d The best model.

Abbreviations: AIC, Akaike information criterion; BMO-MRW, Bruch’s membrane opening-minimum rim width; RNFL-T, retinal nerve fiber layer thickness.

Discussion

Our present data indicated that, compared with BMO-MRW length, the RNFL-T and BMO-MRW volumes showed higher correlation coefficients with total deviations. However, BMO-MRW appeared to be the best predictor of both the threshold of visual sensitivity and the total deviation of visual field in patients with glaucomatous optic neuropathy.

For each sector, Spearman’s rank correlation with threshold of visual sensitivity and with total deviation of visual field showed a different trend in the relationship between each sector and the corresponding visual field function. Leite et al³ previously reported that the strongest relationship was found between superotemporal RNFL-T and inferonasal retinal sensitivity (r=0.314, P<0.001); however, the nasal thickness sector was not significantly related to the corresponding visual field function. Our present data also showed that the superotemporal sector had the strongest relationship, with threshold of visual sensitivity and with total deviation of visual field in all five parameters examined except rim area in total deviation. Furthermore, the nasal sector tended to show the weakest relationships with the examined visual field indices.

Ocular structures in the ONH may be influenced by optic disc size. Although RNFL-T is not affected by optic disc size,^11–13 myopia reportedly produces variable RNFL-T.¹⁴ Using SDOCT, Rao et al¹⁵ found that rim area sensitivity increased and specificity decreased in large optic discs in cases of glaucoma. Similarly, the diagnostic performances of the GDx VCC, HRT II, and Stratus OCT were significantly influenced by optic disc size.¹⁶ Further investigations are necessary to examine whether optic disc size also affects BMO-MRW.

We speculated that volume (length multiplied by area) might be more closely correlated with glaucomatous indices than with length or area separately. Our present results partly supported this hypothesis. Association between the total deviation of visual field and five structural parameters revealed that total deviation was more strongly correlated with RNFL-T volume than with RNFL-T and rim area, and also more strongly correlated with BMO-MRW volume than with rim area (Table 2). Similarly, threshold of visual sensitivity was more strongly correlated with BMO-MRW volume than with rim area (Table 1). However, threshold of visual sensitivity was more significantly correlated with BMO-MRW than with RNFL-T volume.

SDOCT is a useful tool to detect and evaluate glaucoma progression based on retinal structures.^2,3,17 In diagnostic performance for glaucoma detection, the measurement of RNFL-T using SDOCT shows greater sensitivity than Heidelberg retina tomograph (HRT) optic disc measurement.¹⁸ Girkin et al¹³ reported no variations in the glaucoma diagnostic performance of SDOCT between racial groups, and did not find that race was associated with differences in disc area regarding structural parameters of the optic nerve, RNFL-T, and macula.

BMO-MRW may most accurately reflect the amount of retinal nerve fiber from retinal ganglion cells. This para meter is neither perpendicular nor parallel to the peripapillary retinal surface; thus, the BMO-MRW cannot be measured in a plane parallel to the optic disc. However, most programs equipped with SDOCT and HRT, which analyze rim at the ONH, measure the area based on the reference plane parallel to the peripapillary retinal surface. We believe that this measurement method has been a limitation in prior analyses, as has been previously suggested.^4,5

In addition to RNFL-T, macular ganglion cell complex (GCC) thickness has recently been used as a supplement in the diagnosis and evaluation of glaucoma. Although GCC thickness and RNFL-T show similar results in glaucoma detection,^19,20 here we used only RNFL-T as a control for comparison. Further investigations of values derived from BMO-MRW with GCC may further explain the structure–function relationship in respect of the macular GCC.

Conclusion

Our study data confirm that total deviation was more strongly correlated with RNFL-T volume and BMO-MRW volume than with BMO-MRW. However, BMO-MRW appeared to be the best predictor of both the threshold of visual sensitivity and total deviation of visual field, and may be more useful for detecting and evaluating the progression of glaucoma.

Disclosure

The authors report no conflicts of interest in this work.