Abstract
Purpose
To evaluate the relationship between the structural damage as assessed by time-domain optical coherence tomography (OCT) and functional changes in glaucoma.
Methods
In total, 190 patients with normal tension glaucoma or primary open angle glaucoma were included in this study. The thickness of retinal nerve fiber layer (RNFL) around the optic disc and the area of RNFL defect were determined using OCT scans. The relationships between the RNFL thickness or area of the defect and visual field (VF) indices were assessed using the Lowess function, regression analysis and partial Spearman correlation. The differences between these associations depending on the stage of VF damage were further analyzed. Age, optic disc size, refraction, central corneal thickness and the presence of systemic disease were corrected for in order to exclude confounding factors.
Results
A logarithmic scale of RNFL thickness showed a negative linear relationship with VF indices. The area of the RNFL defect showed a weak correlation with the pattern of standard deviation, whereas the remnant RNFL thickness was moderately correlated with the pattern of standard deviation (partial Spearman correlation coefficient, 0.39, -0.47, respectively; p < 0.0001). Many outliers were detected in the Lowess-plotted graphs. Multiplication of the area and the inverted RNFL thickness showed a moderately correlated logarithmic relationship with the VF indices (partial Spearman correlation coefficient, 0.46; 95% confidence interval, 0.34 to 0.57; p < 0.0001). In the severe stage of VF damage, correlation between the area of the RNFL defect and mean deviation was significantly greater than in other stages (partial Spearman correlation coefficient, -0.66; p = 0.02).
Glaucoma is a progressive optic neuropathy characterized by structural changes of the optic nerve and retina that are associated with visual functional defects [1]. However, as Hood and Kardon [1] mentioned, there is no consensus regarding the relationship between the functional loss of sensitivity to light and the structural loss of retinal ganglion cell axons. Many studies have reported various results for this topic, depending on the imaging techniques and statistical methods used, as well as the subject groups investigated and other factors. Leung et al. [2] reported that the coefficient of correlation between visual field (VF) mean deviation (MD) and average retinal nerve fiber layer (RNFL) thickness was 0.79 and that this structure-function relationship was primarily investigated using linear regression analysis. However, a nonlinear function for describing the structure-function relationship was suggested to be better than linear functions, using a confocal laser ophthalmoscope and time-domain (TD) optical coherence tomography (OCT) to quantify the structures [3]. Garway-Heath et al. [4] reported a curvilinear relationship between the decibel differential light sensitivity and both the pattern electroretinogram amplitude and neuroretinal rim area, as well as a linear relationship between 1/Lambert differential light sensitivity and the pattern electroretinogram amplitude and neuroretinal rim area. Another previous study reported that the structure-function relationship is linear when the decibel scale from the VF is converted to a linear (1/Lambert) scale [1,4].
We sometimes encounter patients with similar retinal nerve fiber defect lesions who show different functional deficits, or those who show a clinical discrepancy between the structural damage and the VF defect. In order to account for these inconsistencies, initial quantitative RNFL defect-VF damage analysis is needed. An RNFL defect lesion has a three-dimensional structure, but most studies focused on the RNFL thickness only. Therefore, we divided the concept of the 'RNFL defect lesion' into two components, the depth and the area of the RNFL defect, and then analyzed the structure-function relationship via the two components.
This study was conducted in compliance with the tenets of the Declaration of Helsinki for the use of human subjects in biomedical research, and institutional review board approval was obtained. The medical records of patients with primary open angle glaucoma or normal tension glaucoma were reviewed. All patients underwent full ophthalmic examinations including assessment of the visual acuity, refraction and slit-lamp biomicroscopy, intraocular pressure with Goldmann applanation tonometry, gonioscopy, binocular optic disc evaluation, fundus photography, ultrasonic pachymetry (Humphrey ultrasonic pachometer model 850; Humphrey Instruments, San Leandro, CA, USA), Heidelberg retina tomography (HRT II; Heidelberg Engineering, Heidelberg, Germany), TD-OCT and a VF test (Humphrey Visual Field Analyzer [Carl Zeiss, Dublin, CA, USA] using the Swedish interactive thresholding algorithm standard 30-2 perimetry). Accompanying systemic diseases were evaluated in all patients using a questionnaire. The diagnostic criteria for glaucoma were as follows: open angle, glaucomatous optic disc damage and an abnormal VF test result consisting of a pattern of standard deviation (PSD) of <5%, glaucoma hemifield test results that were outside the normal limits, or both of these findings in at least two consecutive, reliable examinations. Patients under 18 years of age or with best corrected visual acuity of <20 / 40, inflammatory eye disease, exfoliation or pigment dispersion syndrome, ocular trauma, disease involving the retina and choroid or non-glaucomatous optic neuropathy or VF loss were excluded.
Optic disc size was measured using the HRT II. The depth of the RNFL defect was substituted with the remnant RNFL thickness measured with OCT, taking into consideration that the generally remnant thickness is suspected to be reciprocally related to the depth. The real depth of the RNFL defect is calculated by subtraction from the normal RNFL thickness and getting the remnant RNFL thickness. However, the normal RNFL thickness varies individually, and the use of the average RNFL value in previous studies could be one of the confounding factors in this study. The area of the RNFL defect was recorded by the clock hour unit from the OCT results. The perimetry and OCT examinations were performed within a maximum period of 1 month. Statistical analysis was executed using SAS ver. 9.1.3 (SAS Institute, Cary, NC, USA) and R 2.13.2 (Lucent Technologies, Vienna, Austria; http://www.R-project.org). The relationship between the VF indices and the area of the RNFL/RNFL thickness/multiplication of the area and inverted RNFL thickness were plotted using the Lowess function, respectively. Multiplication of the area and inverted thickness was hypothesized to reflect the cube of the RNFL defect. The point of the RNFL thickness or area that showed a rapid change in the slope was calculated. We further analyzed these correlations depending on the stage of the VF defect according to the Hodapp-Parrish-Anderson classifications [5]. Partial Spearman correlation analysis with 95% confidence intervals (CIs) calculated using Fisher's z transformation and the z-test were used to assess whether the extent of correlation showed a difference according to the stage of the VF defect. Multiple linear regression tests were conducted, correcting for optic disc size, age, sex, refractive error, central corneal thickness and systemic disease. When the normality test was not satisfied, the residuals plot was reconfirmed, and we calculated the partial Spearman correlation. A p-value was corrected by Bonferroni's method due to multiple testing and a p-value of less than 0.05 was considered statistically significant.
A total of one hundred and ninety patients with primary open angle glaucoma (n = 81) and normal tension glaucoma (109) were enrolled in the study (Table 1). The clinical findings of the subjects were as follows (average ± standard deviation): age, 56.16 ± 13.44 years; refractive error, -1.88 ± 2.74 diopters; central corneal thickness, 526.36 ± 37.29 µm; disc size, 2.50 ± 0.52; and the rate of accompanying systemic disease, 25.79%. The average RNFL thickness was 62.86 ± 13.99 and the average area of the RNFL defect was 3.65 ± 2.48. The MD index of the VF ranged from 1.6 to -30.86, and the PSD from 1.51 to 18.17. In the classification of the stage, seventy-eight eyes had early defects, forty-nine eyes had moderate defects and sixty-three eyes had severe defects.
In the correlation analysis with MD, the slope of the RNFL thickness rapidly changed around the numerical value of 81 (Fig. 1A). A logarithmic scale of the RNFL thickness showed a negative linear relationship with the MD index of VF (Fig. 1B). When the residual RNFL thickness was plotted, normality was not satisfied and many outliers were found (Fig. 1C and 1D). RNFL thickness and MD were moderately correlated (partial Spearman correlation coefficient, 0.52015; 95% CI, 0.4048 to 0.6175; p < 0.0001). The correlation between RNFL thickness and MD was not significantly different depending on the stage of the VF defect (adjusted p-value of early vs. moderate, early vs. severe and moderate vs. severe stages = 1, 1 and 0.82, respectively).
The slope of the line changed rapidly at the point when the RNFL thickness was 84.25 (Fig. 1E) and was moderately correlated (partial Spearman correlation coefficient, -0.47312; 95% CI, -0.577342 to -0.351454; p < 0.0001). The Lowess plot of the RNFL thickness and PSD suggested a logarithmic relationship (Fig. 1F). The correlation between RNFL thickness and PSD was also not significantly different depending on the stage of the VF defect (adjusted p-value of early vs. moderate, early vs. severe and moderate vs. severe stages = 0.99, 1 and 0.90, respectively).
The area of the RNFL defect showed a negative linear relationship with the MD and this correlation was moderate (partial Spearman correlation coefficient, -0.50897; 95% CI, -0.60793 to -0.392089; p < 0.0001) (Fig. 2A). However, in the relationship with the PSD, the area of the RNFL defect showed a weakly logarithmic scale correlation, as shown by the Lowess plot (partial Spearman correlation coefficient, 0.39766; 95% CI, 0.267438 to 0.512024; p < 0.0001) (Fig. 2B and 2C). The slope of the line became stable after the area of the RNFL defect lesion passed 6 o'clock. In the glaucomatous stages, the RNFL defect area in the severe stage was more negatively correlated with the MD index than in the moderate stage (partial Spearman correlation coefficient, -0.66; adjusted p-value, 0.02). There was no significant difference in the correlation between the early and moderate stages (adjusted p-value, 1). The correlation between the area of the RNFL defect and PSD was not significantly different compared to the glaucomatous stage (adjusted p-value of early vs. moderate, early vs. severe and moderate vs. severe stages = 0.99, 1 and 0.99, respectively).
The multiplication values showed a logarithmic relationship with MD and were moderately correlated (partial Spearman correlation coefficient, 0.56314; 95% CI, -0.653602 to -0.4543712; p < 0.0001). When calculated with PSD, the multiplication value also showed a logarithmic relationship and a moderate correlation (partial Spearman correlation coefficient, 0.46729; 95% CI, 0.344899 to 0.572344; p < 0.0001).
In the glaucoma field, there is no absolute agreement concerning which models, linear or nonlinear, yield the best correlation between structural and functional damage. The structure-function relationship has been evaluated in various ways, including confocal scanning laser ophthalmoscopy, scanning laser polarimetry, OCT and histologic studies [6,7,8]. The onset timing of RNFL defect and VF defects is still controversial. Several studies determined that the structural alterations appeared before the development of VF defects [9,10]. However, as the disease progresses, functional measures may prove more sensitive for detecting progressive damage because alterations in the structure of the optic disc approach a state of absolute cupping, which, when reached, indicates the end stage of the disease [11]. The concept of 'structure and functional correlation' is vague since it may simply reflect the correlation itself, the relationship between the preceding and the forwarding or a mathematical correlation.
In our study, the authors divided the concept of the 'RNFL defect' into two aspects, RNFL thickness and the area of the RNFL defect. We detected a moderately correlated logarithmic relationship between the RNFL thickness and the VF indices, when corrected for age, central corneal thickness, refractive error, optic disc size and systemic disease. Additionally, the PSD index was more strongly correlated with the RNFL thickness than with the area of the RNFL defect, even when considering the scale difference between the thickness and area. The multiplication value also demonstrated a moderately correlated logarithmic relationship with the VF indices. When analyzing these correlations according to glaucomatous stages, the correlation between the area of the RNFL defect and MD was significantly stronger in the severe stage than in the mild and moderate stages, whereas the correlation between RNFL thickness and MD did not differ according to glaucomatous stage. Kanamori et al. [12] reported a significant relationship between the MD and RNFL thickness in all parameters excluding the 3-o'clock area. Recently, Leite et al. [13] reported that the relationship was strongest for superotemporal RNFL thickness and inferonasal sensitivity (R2 = 0.314, p < 0.001), using the cirrus spectral domain (SD) OCT. In addition, Lee et al. [14] reported no correlation between retinal mean sensitivity (MS) and RNFL thickness in normal, glaucoma suspects. However, in glaucomatous eyes, the association with TD-OCT and the correlation of superior RNFL thickness with inferior MS was significantly better with SD-OCT than with TD-OCT in both linear and logarithmic regression models. In this study, the slope of MD more rapidly dropped at an RNFL thickness of 81, and PSD slope became stable after the RNFL thickness decreased below 84.25, implying that an RNFL thickness of around 80 was a meaningful value.
There are some limitations to this study: the cross-sectional study design and the small number of the populations involved. In addition, the RNFL thickness in the OCT results did not exactly reflect the changes in the pure RNFL because the OCT thickness was the sum of the thickness of the RNFL and the thickness of other tissues, such as blood vessels and glial tissue [1]. Moreover, attempting to correct this can result in other confounding factors, and the RNFL thickness measured using OCT has been known to be useful in accurately and objectively detecting glaucoma, even in its early stages [15]. The multiplication value also did not reflect the real value of the cube of the RNFL defect. In this study, however, we focused on the concept of the 'three-dimensional RNFL defect' and compared the influence of the depth and area of the RNFL defect alone. Both primary open angle glaucoma and normal tension glaucoma patients were enrolled in this study. It will be difficult to determine any difference between these two groups in correlation analysis because the position or pattern of retinal nerve fiber damage has been suggested to be different in the two groups. Patients with normal tension glaucoma appear to have deeper, more localized scotomas, even though there have been conflicting reports [16].
Interestingly, we detected many outliers on a Lowess plot. There are some possible reasons for this: 1) the number of study subjects, 2) OCT measurement error, 3) individual variation of RNFL thickness, 4) the characteristic that VF damage may present later in subjects with a thick RNFL than in those with a thin RNFL, and 5) selection bias. Based on these results, caution should be exercised because more patients with an RNFL-VF discrepancy may exist than expected in the clinical setting; further research regarding the RNFL-VF discrepancy is necessary. Most studies that have evaluated the structure-functional relationship in glaucoma were cross-sectional comparative studies, but longitudinal studies are needed to define the relationship, considering the disease characteristics.
In conclusion, RNFL thickness had a negative logarithmic correlation with the VF indices, and was more strongly correlated with the VF indices than the area of the RNFL defect, as measured by OCT. An RNFL thickness of around 80 microns was the point at which the slope of the MD and PSD changed. In addition, there existed many outliers in the structure-function relationship graphs.
Figures and Tables
References
1. Hood DC, Kardon RH. A framework for comparing structural and functional measures of glaucomatous damage. Prog Retin Eye Res. 2007; 26:688–710.
2. Leung CK, Chan WM, Hui YL, et al. Analysis of retinal nerve fiber layer and optic nerve head in glaucoma with different reference plane offsets, using optical coherence tomography. Invest Ophthalmol Vis Sci. 2005; 46:891–899.
3. Leung CK, Chong KK, Chan WM, et al. Comparative study of retinal nerve fiber layer measurement by StratusOCT and GDx VCC, II: structure/function regression analysis in glaucoma. Invest Ophthalmol Vis Sci. 2005; 46:3702–3711.
4. Garway-Heath DF, Holder GE, Fitzke FW, Hitchings RA. Relationship between electrophysiological, psychophysical, and anatomical measurements in glaucoma. Invest Ophthalmol Vis Sci. 2002; 43:2213–2220.
5. Hodapp E, Parrish RK, Anderson DR. Clinical decisions in glaucoma. St Louis: Mosby;1993. p. 84–126.
6. Chauhan BC, Blanchard JW, Hamilton DC, LeBlanc RP. Technique for detecting serial topographic changes in the optic disc and peripapillary retina using scanning laser tomography. Invest Ophthalmol Vis Sci. 2000; 41:775–782.
7. Harwerth RS, Carter-Dawson L, Shen F, et al. Ganglion cell losses underlying visual field defects from experimental glaucoma. Invest Ophthalmol Vis Sci. 1999; 40:2242–2250.
8. Bellezza AJ, Rintalan CJ, Thompson HW, et al. Deformation of the lamina cribrosa and anterior scleral canal wall in early experimental glaucoma. Invest Ophthalmol Vis Sci. 2003; 44:623–637.
9. Sommer A, Katz J, Quigley HA, et al. Clinically detectable nerve fiber atrophy precedes the onset of glaucomatous field loss. Arch Ophthalmol. 1991; 109:77–83.
10. Quigley HA, Katz J, Derick RJ, et al. An evaluation of optic disc and nerve fiber layer examinations in monitoring progression of early glaucoma damage. Ophthalmology. 1992; 99:19–28.
11. Cioffi GA, Liebmann JM, Johnson CA, Weinreb RN. Structural-functional relationships of the optic nerve in glaucoma. J Glaucoma. 2000; 9:3–4.
12. Kanamori A, Nakamura M, Escano MF, et al. Evaluation of the glaucomatous damage on retinal nerve fiber layer thickness measured by optical coherence tomography. Am J Ophthalmol. 2003; 135:513–520.
13. Leite MT, Zangwill LM, Weinreb RN, et al. Structure-function relationships using the Cirrus spectral domain optical coherence tomograph and standard automated perimetry. J Glaucoma. 2012; 21:49–54.
14. Lee JR, Jeoung JW, Choi J, et al. Structure-function relationships in normal and glaucomatous eyes determined by time- and spectral-domain optical coherence tomography. Invest Ophthalmol Vis Sci. 2010; 51:6424–6430.
15. Sihota R, Sony P, Gupta V, et al. Diagnostic capability of optical coherence tomography in evaluating the degree of glaucomatous retinal nerve fiber damage. Invest Ophthalmol Vis Sci. 2006; 47:2006–2010.
16. Gramer E, Althaus G, Leydhecker W. Site and depth of glaucomatous visual field defects in relation to the size of the neuroretinal edge zone of the optic disk in glaucoma without hypertension, simple glaucoma, pigmentary glaucoma. A clinical study with the Octopus perimeter 201 and the optic nerve head analyzer. Klin Monbl Augenheilkd. 1986; 189:190–198.