|
|
||||||||
1From the Departments of Psychology and 2Ophthalmology, Columbia University, New York, New York; and the 3Department of Ophthalmology, University of Iowa and Veterans Administration, Iowa City, Iowa.
| Abstract |
|---|
|
|
|---|
METHODS. Fifteen patients with asymmetrical glaucoma, whose better eye was normal or near normal (mean deviations better than –3 dB) on SAP, were tested. SITA 24-2 standard and OCT RNFL thickness measures were made on three to five different occasions and the mean values were obtained. For each eye, the mean SAP loss was calculated for an upper and lower arcuate field region by averaging the loss in relative sensitivity on a linear scale. The average RNFL thickness for corresponding arcuate sectors of the lower and upper optic disc was obtained for each eye. A linear model was fitted to the plots of RNFL thickness versus SAP loss. According to the linear model, the RNFL thickness R = soT + b, where T is the SAP sensitivity loss relative to age-matched normal eyes (linear scale), (so + b) is the RNFL thickness in the healthy/normal state (T = 1), and b is the residual RNFL thickness measured when all sensitivity and all axons are lost.
RESULTS. The model provided a reasonable fit to the data with best fitting values of (so; b) of (upper field: 80.6 µm; 50.5 µm) and (lower field: 67.4 µm; 50.5 µm) and (upper field: 78.8 µm; 54.9 µm; r = 0.82) and (lower field: 59.2 µm; 61.5 µm; r = 0.70) for two different methods of best fit.
CONCLUSIONS. A linear model that relates RNFL thickness to losses in SAP sensitivity describes the results for arcuate regions of glaucomatous visual fields. The linear model provides a framework for assessing the relative efficacy of structural and functional tests throughout the course of the disease.
The advent of automated, noninvasive techniques for measuring retinal nerve fiber layer thickness (RNFL) has sparked renewed interest in the relationship between structural and functional glaucomatous losses. In the past 6 years, a large number of empiric studies have related functional losses measured with static automated perimetry (SAP) to RNFL losses measured with these automated techniques. The analyses of these data have been largely descriptive, and there is a debate about the form of the function relating functional and structural losses. (For references, see the review by Garway-Heath2 and a recent article by Bowd et al.3 ) However, the theoretical treatment of these data has been relatively limited.2 Recently, Hood4 proposed a simple model that predicts the relationship between RNFL thickness and loss of SAP sensitivity. This model, based on a model5 proposed to explain the relationship between multifocal visual evoked amplitudes and SAP losses, assumes that the local loss in RNFL thickness is linearly related to the loss in SAP sensitivity, when SAP sensitivity loss is expressed on a linear scale. However, a complete loss in sensitivity does not result in an RNFL thickness of zero; rather, it is associated with a finite RNFL thickness. In particular, there is some minimum value beyond which the thickness cannot be reduced because glial cells or other factors do not contribute to the number of axons. According to this model, the portion of the RNFL thickness due to retinal ganglion cell (RGC) axons is linearly related to the proportion of RGC axons ("response") remaining.
Hood4 reviewed the data from published studies6 7 8 9 10 that used optical coherence tomography (OCT) to measure RNFL thickness and concluded that these data were consistent with the proposed model. However, there are problems with evaluating this model with the published data. First, the results from both OCT and SAP may vary from day-to-day, due to measurement variability. Existing studies typically present the cross-sectional results for a single OCT and SAP test. Second, the best way to evaluate the model is to compare local sensitivity to local RNFL loss. Most of the studies in the literature present the SAP data averaged over a large region (e.g., hemifield or full-field of 24-2 SAP). Finally, in existing studies, the SAPs value for a region of the visual field are usually expressed as an arithmetic average of dB units for each location within the region considered, but, as detailed in the Methods section, according to the model, these values should be antilogged before averaging and then logged again after averaging.1 5 11
In the current study, we fit the linear model to data from patients with asymmetrical glaucoma tested on multiple occasions with RNFL thickness (OCT) and SAP. To minimize the effects of measurement variability, multiple measures of OCT retinal nerve fiber layer thickness recorded on different test days were averaged, as were the antilog values of the SAP measures for the same days. Regional changes were examined by using the mapping between SAP arcuate regions and corresponding retinal nerve fiber layer sectors proposed by Garway-Heath et al.12 The data have been presented in abstract form (Hood DC, et al. IOVS 2007;48:ARVO E-Abstract 490).
| Methods |
|---|
|
|
|---|
|
For each test date, the visual field threshold data and the OCT RNFL data were divided into six sectors, based on the schema proposed by Garway-Heath et al.12 (Fig. 1) . For visual field threshold data, the decibel levels in each location of the total deviation field were converted to a linear scale (e.g., 0 dB converted to 1.0 and –30 dB to 0.001) before averaging the data within each sector.1 5 11 To understand the rationale for averaging the values on a linear scale rather than on the log (decibel) scale, one must consider a region of the retina in which the left half is normal (0 dB total deviation and a normal complement of RGCs) and in which the right half has no RGCs remaining and a maximum loss in sensitivity (a total deviation of
–30 dB). If we take the average of the decibel levels, we get –15 dB, which is 1/30 on a linear scale. In contrast, if we take the average on the linear scale, we get 0.5 (1.001/2), which is –3 dB on the decibel scale. The model predicts that the OCT RNFL thickness due to RGC axons for the entire region should be one half the normal thickness.
|
| Results |
|---|
|
|
|---|
|
![]() | (1) |
–30 dB). To derive a predicted curve from equation 1 , we estimated the starting (so + b) and ending (b) points of the curve. In particular, a curve was fitted to the data by setting (so + b) equal to the median of the arcuate RNFL thicknesses for the better eye. (so + b) was calculated separately for the superior (131.1 µm) and inferior (117.9 µm) fields. b was set equal to the median (50.5 µm) of the RFNL thickness for field loss greater than –10 dB for both superior and inferior field regions. (For field losses worse than –10 dB, the RNFL thickness reaches an asymptotic value.)
Figures 2C and 2D show the same data on linear coordinates. RNFL thickness is plotted against the antilog of the total SAP deviations (Fig. 2D) . The straight lines are the fit of
![]() | (2) |
The structural–functional analysis is based on a mapping of specific visual field locations to a particular sector of the circular RNFL scan, namely the Garway-Heath et al.2 map in Figure 1 . To get a measure of whether this map was optimal for the eyes in this study, the following analyses were performed. First, the width of the Garway-Heath et al. arcuate disc sector was held constant (40°) and rotated in 5° steps around the optic disc. For each sector location, a correlation coefficient r between the SAP field loss and the RNFL for that location was determined with our linear model. Figure 3 shows the r-value on the y-axis as a function of the position of the center of the 40° sector for correlations with the superior arcuate SAP field values and inferior arcuate SAP field values. The best correlations are close to the center of the sectors in the map by Garway-Heath et al. For the superior arcuate field region, the best correlations are centered approximately 5° nasal from the center of the inferior arcuate disc sector in Garway-Heath et al. The agreement with the map of Garway-Heath et al. is excellent.
|
|
| Discussion |
|---|
|
|
|---|
First, the model does not predict a single curve, but rather it predicts that there will be a family of curves. Figure 5 shows the same data as in Figure 2 , with the dashed curves representing the boundaries of a family of predicted curves. To understand how these boundaries were obtained, one must recall that the model assumes that the measure of RNFL thickness has two components, one, soT, which represents the thickness of the retinal axons associated with a given relative sensitivity, T, and the other, b, which is the residual RNFL thickness measured when all the axons are lost. This residual portion includes glial cells and perhaps limitations imposed by the algorithm that determines the RNFL layer. In individuals with normal visual sensitivity (T = 1.0), the RNFL thickness is the sum of so + b, where so is the thickness of the axon portion in the normal healthy eye. There is a wide range of values for (so + b) as defined by the 95% confidence interval for normal RNFL thickness as shown by the green region in Figure 4C . Assuming for the moment that b is the same in different eyes, this confidence interval provides a range of normal (so + b) values. (For the purposes of this example, the effects of age on so, and thus on the confidence interval, were not taken into consideration.) The upper and lower boundaries of the confidence interval each provide the parameter (so + b) for a theoretical curve. Figure 5 shows the data from Figure 2 with the predicted curves associated with the upper and lower limits of (so + b) estimated from the green region in Figure 4E . In particular, the upper curve describes the predicted course of glaucomatous progression in a patient who started with a relatively large so, whereas the lower curve shows the predicted curve in a patient with an so that was relatively small when normal. Regardless of the initial RNFL thickness, all curves have the same common shape, meaning that there is a loss in SAP sensitivity that is proportional to RNFL thickness attributable to RGC axons. Note that the linear model, combined with the normal confidence interval, predicts that most of our data points should fall between these curves.
|
Although the model (equations 1 and 2) explicitly relates RNFL thickness to SAP loss, there is an implicit assumption that the local RNFL thickness is directly proportional to the number of RGC axons.4 That is, if local RNFL thickness is reduced by one half, then the number of RGC axons is reduced by one half. If we assume that the local number of axons is proportional to RNFL thickness, then the model predicts that the loss in the relative number of local RGCs is linearly related to local sensitivity loss (in linear terms). On the one hand, this assumption is in apparent contradiction with the conclusions in the monkey13 14 15 and human15 16 17 studies in which RGCs were counted in postmortem tissue. On the other hand, it is in agreement with theoretical treatments of the human data that argue for a linear relationship.1 2 5 18 19 For example, the best articulated model of the human data by Swanson et al.19 predicts that the relationship is linear for eccentricities beyond 5° to 10°. There are several differences among the methodologies in these studies,2 19 and further work is needed before we can account for the differences in conclusions.
A related question is the relationship between RNFL thickness and SAP sensitivity in individuals with normal vision. The linear model articulated in this study is specified only for losses in sensitivity (i.e., D
0 dB or T
1.0); it makes no prediction about individuals who may have greater than normal sensitivity. We know that both SAP sensitivity and OCT RNFL thickness decrease with age. What remains to be determined is the quantitative relationship between these changes in the arcuate regions, as well as the nature of the changes, if any, within groups of normal individuals of the same age.
Clinical Implications
The model provides a framework in which to determine the stage of a diseased eye. In particular, as glaucoma progresses, the data points associated with a particular eye should move along a single theoretical curve. We intend to test this prediction by observing the patients over time. In addition, at any given time point, the model combined with information about the variability of the SAP and RNFL thickness measures provides a forum for considering the optimal use of these tests in the clinic. For detecting early glaucomatous damage, it is often said that SAP does not show a statistically significant defect until structural changes have taken place. For example, Kerrigan-Baumrind et al.17 reported that a loss of 25% to 35% of the RGCs was associated with a statistically significant SAP field loss. According to our model, a loss of 25% to 35% of the RGCs would result in a loss in local sensitivity between –1.2 and –4.6 dB, within the 5% confidence interval (
–5 dB for local SAP points). However, this assumption does not necessarily mean that the RNFL measure is more sensitive than SAP for detecting early damage. In practice, which test is more sensitive to early damage will depend on the relative variation of each normal control measure and the initial RNFL thickness (so + b) when the eye is healthy.
Consider two patients, each with a 50% RGC loss, but with different so values. According to the linear model, each patient would show a proportional –3-dB field loss and normal SAP test results. In contrast, a 50% RGC loss and the consequent 50% decrease in so in the patient with the smaller so would result in a much smaller RNFL thickness than that in the patient with the larger so. Thus, the patient with the smaller initial RNFL thickness may show an abnormal RNFL thickness, whereas the other is still in the normal range. Thus, a significant change in the OCT RNFL thickness may occur before a significant change in the SAP in the patient with the smaller initial RNFL thickness. Conversely, the patient with the thicker initial RNFL may show abnormalities on the SAP before the RNFL thickness drops below the 5th percentile of the normal range. In general, to determine which of two tests is the more sensitive, one must know both the function that relates the two tests and the variability in test scores in healthy individuals.1 5 14
The model also has implications for the use of the OCT RNFL thickness to detect progression of glaucoma. The model indicates that the RNFL thickness approaches the asymptotic value, b, as the SAP loss approaches –10 to –15 dB. The value of b, approximately 50 µm, provides an estimate for the lower bound of the OCT RNFL thickness. On average, the thickest part of the RNFL profile is <200 µm (Fig. 1 , bottom). Thus, the largest range typically encountered spans a factor of less than 4. According to the model, a –15-dB loss would reduce a 200-µm thickness to 54.7 µm {[150 µm x 0.03 (antilog of –15 dB)] + 50 µm}. This value is within the measurement error of the lower bound of 50 µm. Further, this is close to a best-case scenario; in practice, local SAP losses of more than –10 dB will not yield an RNFL thickness detectably different from the base level, b. Thus, the RNFL thickness is of limited use for observing regions with extensive damage. In contrast, the SD of SAP under conditions of a –10-dB loss is very large.20 For both tests, progression is probably best observed by examining relatively healthier regions of the same eye.
| Summary |
|---|
|
|
|---|
| Acknowledgements |
|---|
| Footnotes |
|---|
Submitted for publication November 22, 2006; revised January 6 and February 9, 2007; accepted May 11, 2007.
Disclosure: D.C. Hood, None; S.C. Anderson, None; M. Wall, None; R.H. Kardon, None
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: Donald C. Hood, Department of Psychology, 406 Schermerhorn Hall, Columbia University, New York, NY, 10027; dch3{at}columbia.edu.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M A El-Dairi, S Holgado, S G Asrani, L B Enyedi, and S F Freedman Correlation between optical coherence tomography and glaucomatous optic nerve head damage in children Br J Ophthalmol, October 1, 2009; 93(10): 1325 - 1330. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. C. Hood, S. C. Anderson, M. Wall, A. S. Raza, and R. H. Kardon A Test of a Linear Model of Glaucomatous Structure-Function Loss Reveals Sources of Variability in Retinal Nerve Fiber and Visual Field Measurements Invest. Ophthalmol. Vis. Sci., September 1, 2009; 50(9): 4254 - 4266. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. A. Lagreze, M. Gaggl, M. Weigel, J. Schulte-Monting, A. Buhler, M. Bach, R. D. Munk, and T. A. Bley Retrobulbar Optic Nerve Diameter Measured by High-Speed Magnetic Resonance Imaging as a Biomarker for Axonal Loss in Glaucomatous Optic Atrophy Invest. Ophthalmol. Vis. Sci., September 1, 2009; 50(9): 4223 - 4228. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Turpin, G. P. Sampson, and A. M. McKendrick Combining Ganglion Cell Topology and Data of Patients with Glaucoma to Determine a Structure-Function Map Invest. Ophthalmol. Vis. Sci., July 1, 2009; 50(7): 3249 - 3256. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Wang, D. C. Hood, J.-S. Cho, Q. Ghadiali, G. V. De Moraes, X. Zhang, R. Ritch, and J. M. Liebmann Measurement of Local Retinal Ganglion Cell Layer Thickness in Patients With Glaucoma Using Frequency-Domain Optical Coherence Tomography Arch Ophthalmol, July 1, 2009; 127(7): 875 - 881. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. K. Horn, C. Y. Mardin, R. Laemmer, D. Baleanu, A. M. Juenemann, F. E. Kruse, and R. P. Tornow Correlation between Local Glaucomatous Visual Field Defects and Loss of Nerve Fiber Layer Thickness Measured with Polarimetry and Spectral Domain OCT Invest. Ophthalmol. Vis. Sci., May 1, 2009; 50(5): 1971 - 1977. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. A. Medeiros, L. M. Alencar, L. M. Zangwill, C. Bowd, G. Vizzeri, P. A. Sample, and R. N. Weinreb Detection of Progressive Retinal Nerve Fiber Layer Loss in Glaucoma Using Scanning Laser Polarimetry with Variable Corneal Compensation Invest. Ophthalmol. Vis. Sci., April 1, 2009; 50(4): 1675 - 1681. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. S. Harwerth, J. L. Wheat, and N. V. Rangaswamy Age-Related Losses of Retinal Ganglion Cells and Axons Invest. Ophthalmol. Vis. Sci., October 1, 2008; 49(10): 4437 - 4443. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Ferreras, L. E. Pablo, D. F. Garway-Heath, P. Fogagnolo, and J. Garcia-Feijoo Mapping Standard Automated Perimetry to the Peripapillary Retinal Nerve Fiber Layer in Glaucoma Invest. Ophthalmol. Vis. Sci., July 1, 2008; 49(7): 3018 - 3025. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Machida, Y. Gotoh, Y. Toba, A. Ohtaki, M. Kaneko, and D. Kurosaka Correlation between Photopic Negative Response and Retinal Nerve Fiber Layer Thickness and Optic Disc Topography in Glaucomatous Eyes Invest. Ophthalmol. Vis. Sci., May 1, 2008; 49(5): 2201 - 2207. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |