|
|
||||||||
1From the Department of Ophthalmology, Caritas Medical Centre, Hong Kong, Peoples Republic of China; the 2Departments of Ophthalmology and Visual Sciences and 4Physiology, Faculty of Medicine, The Chinese University of Hong Kong, Hong Kong, Peoples Republic of China; and the 3Department of Industrial and Manufacturing Systems Engineering, The University of Hong Kong, Hong Kong, Peoples Republic of China.
| Abstract |
|---|
|
|
|---|
METHODS. Eighty-nine subjects (27 who had healthy eyes, 21 who were glaucoma suspect, 41 who had glaucoma) were enrolled in this cross-sectional study. RNFL thickness was measured using the StratusOCT and the GDx VCC, and visual field (VF) was examined using the Humphrey VF analyzer. The relationship between RNFL thickness and VF sensitivityexpressed in terms of mean deviation (MD) in decibel (dB) scale, unlogged 1/lambert (L), and Advanced Glaucoma Intervention Study (AGIS) and Collaborative Initial Glaucoma Treatment Study (CIGTS) VF scoreswere evaluated with linear and nonlinear regression models. Coefficient of determination (R2) was calculated, and regression models were compared using the Akaike information criterion and the F test.
RESULTS. In plotting MD against RNFL thickness, curvilinear regression models demonstrated the best fit, whereas linear regression attained the best associations when VF sensitivity was expressed in 1/L. However, when healthy subjects were excluded from the analyses, the second-order polynomial was better than linear regression in describing the relation between 1/L and GDx VCCmeasured RNFL thickness. Regression profiles between AGIS/CIGTS VF scores and RNFL thickness were best described in the linear and the first-order inverse models for GDx VCC and StratusOCT RNFL measurements, respectively. In general, StratusOCT RNFL measurements achieved higher associations with visual function in all the respective regression analyses than did GDx VCC.
CONCLUSIONS. Description of structure/function relationships in glaucoma depends on the choice of perimetry scale, the type of RNFL measuring device, and the characteristics of the studied groups. The higher association with visual function in StratusOCT RNFL measurements compared with that in GDx VCC suggested optical coherence tomography might be a better approach for evaluating structure/function relationships. Curvilinear regression profiles found between StratusOCT RNFL thickness and MD/VF scores provide an explanation for those longitudinal observations, showing that VFs with higher AGIS/CIGTS VF scores or worse MD at baseline are at higher risk for deterioration. Regression analysis of the structure/function profile could provide important information in the assessment of the trend and pattern of glaucoma progression.
Information regarding the nature of progression in glaucoma has been essentially derived from the longitudinal VF analyses.6 7 8 9 10 11 12 To avoid the confounding effect of variability, a number of analytical procedures were suggested to determine the true progression. The Humphrey visual field analyzer (Humphrey Field Analyzer II; Humphrey Instruments, Dublin, CA) has a built-in "change analysis" option that performs linear regression analysis of the VF MD. Major clinical trials such as the Advanced Glaucoma Intervention Study (AGIS) and the Collaborative Initial Glaucoma Treatment Study (CIGTS) use defect classification scoring systems to provide a discrete score, between 0 and 20, for each VF result.13 14 The AGIS and the CIGTS define progression as an increase in VF score of 4 and 3, respectively. However, the scoring systems can be arbitrary, and the scale may not be linear. A change in score from 1 to 5 may not be equal to a change from 16 to 20. Analyzing the structure/function profiles in terms of VF scores may reveal the nature of the scaling of these scoring systems and may provide better understanding regarding the trend of disease progression observed in these trials. Using different nonlinear and linear regression models, we investigated and compared the relationship between VF sensitivityexpressed in terms of decibel unit (MD), unlogged (1/Lambert [1/L]), and AGIS/CIGTS VF scoresand RNFL thickness, measured with StratusOCT and GDx VCC, in subjects with glaucoma.
| Materials and Methods |
|---|
|
|
|---|
All subjects underwent full ophthalmic examination, including visual acuity, refraction, intraocular pressure measurement with Goldmann tonometry, dilated fundus examination with stereoscopic biomicroscopy of optic nerve head under slit-lamp, and indirect ophthalmoscopy. Inclusion criteria were best-corrected visual acuity of not worse than 20/40 and spherical refractive error within the range of 6.00 diopters (D) to +3.00 D. Subjects were excluded if they had history of retinal disease, surgery, laser procedures, diabetes mellitus, or neurologic disease. Healthy subjects visited the clinic during the same recruitment period and had no ocular or intraocular diseases. In particular, they had no VF defects based on the Humphrey VF results, no structural optic disc abnormalities, and no history of intraocular pressure exceeding 21 mm Hg. Glaucoma-suspect subjects were those with ocular hypertension, pre-perimetric glaucoma, or both. Ocular hypertension was defined as intraocular pressure greater than 21 mm Hg measured in at least three separate visits. Pre-perimetric glaucoma was diagnosed when subjects had been found to have asymmetric cup disc ratios of >0.2 and early glaucomatous optic disc changes, including thinning of the neuroretinal rim and notching. However, all subjects in the glaucoma-suspect group had normal VF results, as did the control group with healthy eyes. Glaucoma had been diagnosed in subjects based on the presence of VF defects. Among the 41 patients with glaucoma, 21 had primary open-angle glaucoma, 13 had normal-tension glaucoma, 6 had primary angle-closure glaucoma, and 1 had uveitic glaucoma.
Visual Sensitivity Measurements
Visual fields were evaluated with the Humphrey visual field analyzer (central 24-2 SITA fast program; Humphrey Instruments). Humphrey VF test results were considered reliable when the fixation loss was <20% and the false-positive and false-negative error rates were <25%. A VF defect was defined as three or more significant (P < 0.05) nonedge-contiguous points, with at least one at the P < 0.01 level on the same side of the horizontal meridian in the pattern deviation plot and classified outside normal limits in the glaucoma hemifield test. Any detected field defect had to be confirmed in at least one other attempt to be classified as abnormal. Average visual sensitivity was expressed in 3 different forms: MD, unlogged 1/L scale, and AGIS and CIGTS VF scores. MD represents an age-matched average visual sensitivity index expressed in decibels. It is calculated based on the total deviation plot in the Humphrey VF analysis. Differential light sensitivity at each tested location is measured in decibels, where the differential light sensitivity (dB) = 10 x log10 (Lmax/(Lt Lb), with Lmax is the maximal stimulus luminance, Lt is the stimulus luminance at threshold, and Lb is the background luminance. For simplicity, this relationship can be written as 10 x log10 (1/L). The unlogged 1/L at each tested location was calculated by dividing the decibel unit by 10 and then unlogging it. The average value was then evaluated. Each VF in the glaucoma group was also analyzed to determine the AGIS/CIGTS VF scores. Both scoring systems use 20-interval scales, with 0 representing no defect and 20 representing severe damage. The scorings of the VFs were computed according to the original descriptions in the respective trials.13 14
Optical Coherence Tomography Measurements
StratusOCT (Carl Zeiss Meditec) with software version 3.0.1 was used in the study. Detailed descriptions of the optical principles and applications of optical coherence tomography (OCT) have been described.15 The RNFL (3.4) scan (with 512 scan points) was selected for RNFL measurement. RNFL thickness was measured by averaging the results of 3 sequential circular scans with 3.4-mm diameters centered at the optic nerve head. The parameter, average RNFL thickness, is obtained in the analysis printout. A good-quality scan was defined as one with a signal-to-noise ratio of >35, 100% accepted A-scans, and well-delineated anatomic boundaries. A subject would be excluded from the study if the OCT image was not of good quality after 3 attempts. All OCT scans were of good quality, and no subject was excluded.
Scanning Laser Polarimetry Measurements
Scanning laser polarimetry (SLP) measurements were performed with the use of GDx VCC (Laser Diagnostic Technologies). Images were analyzed with software version 5.5.0. The principles and applications of SLP have been described elsewhere.16 The parameter TSINT average, equivalent to the average RNFL thickness, in the analysis printout was used in the study. Nine eyes were excluded in this study because of suboptimal scanning quality secondary to poor fixation, motion artifacts, or overilluminated images.
All measurements (VF sensitivity and OCT/SLPmeasured RNFL thickness) were obtained within the same study period (AugustNovember 2004).
Statistical Analysis
Differences in age, refraction, and VF MD among the diagnostic groups were evaluated with one-way ANOVA. Bonferroni correction was used for multiple comparisons. In all analyses, P < 0.05 was considered statistically significant.
Relationships between average RNFL thickness and visual sensitivity were evaluated with linear and nonlinear regression analyses. The linear model (y = ax + b) was compared with four common nonlinear models, including the second-order polynomial (y = ax2 + bx + c), the third-order polynomial (y = ax3 + bx2 + cx + d), the first-order inverse (y = a ln(x) + b), and the logarithmic regressions (y = a log(x) + b). In regression analysis, the goodness-of-fit of any particular regression model is expressed as the coefficient of determination, R2, which indicates how much of the total variation in the dependent variable can be accounted for by the regression function. However, it is not possible to determine whether model A is more correct than model B in describing the relationship profile based on the value of R2 because the model with more parameters often has a higher R2 than the model with fewer parameters. The extra-sum-of-square F test and the Akaike information criterion (AIC) are two mathematical approaches that take model complexity (the number of data points and the number of parameters) into account in calculating the F ratio and the AIC difference, respectively, to determine which regression model to accept.
The F ratio is calculated by comparing the relative difference in the sum-of-square with the relative difference in the degrees of freedom of two nested (related) models. When comparing model A (the simpler model) and model B (the more complicated model),
![]() |
The AIC adopts another mathematical approach that can be used to compare nested and nonnested models. The concepts and applications of AIC have been described in detail by Burnham and Anderson.17 Because the logic behind AIC is not hypothesis testing, there is no P value in the calculation. The model with the lowest AIC is more likely to be correct, and the difference in corrected AIC (AICc) between model A and model B is defined as:
![]() |
![]() |
The evidence ratio is defined as the probability that model A is correct divided by the probability that model B is correct. By convention, it is arranged to give a value greater than 1. Therefore, the evidence ratio should be interpreted with reference to the values of the probability of correct fitting. An evidence ratio of 2 describes one model as twice as likely to be correct as another. The difference in AICc is linearly correlated with the evidence ratio. In the present study, the linear model was compared with each of the nonlinear regression models. The probability of correct fitting, the difference in AICc, the evidence ratio, and the P value of the F test (in nested models only) were computed. All these parameters, though expressed in different forms, are in agreement with each other to determine which one is the best-fit regression model.
| Results |
|---|
|
|
|---|
|
|
|
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
The relationship between RNFL thickness and VF MD was found to be better fit with all the nonlinear models examined (best fit with the third- and the second-order polynomials for GDx VCC and StratusOCT RNFL measurements, respectively) compared with the linear model (Fig. 1) . However, because healthy subjects and subjects with glaucoma were included in the analysis, one could argue that the flatter portion of the curve may represent the wide variations in RNFL thickness among healthy persons and may give the false impression of a functional reserve in the early stage of the disease (using OCT as an example, average RNFL thicknesses have been reported from 80 to 150 µm in various studies18 19 20 21 ). To eliminate this factor, therefore, we excluded healthy persons in the second analysis (Fig. 3) . It was confirmed that the nonlinear models were better than the linear model (best fit with the logarithmic and the second-order polynomials for the GDx VCC and the StratusOCT RNFL measurements, respectively) in analyzing the relationship between RNFL thickness and VF MD. These regression models describe a curvilinear structure/function relationship, suggesting that the progression of visual field loss, when it is expressed in MD (dB), increases during the course of the disease. This observation is also consistent with postmortem histologic measurement in patients with glaucoma, indicating that at least 25% to 35% of retinal ganglion cells were lost before abnormalities were statistically detected through automated visual field testing.22
When the visual field scale was expressed in unlogged 1/L, the relationship was best fit with linear models in GDx VCC and StratusOCT RNFL measurements (Fig. 2) . However, when healthy subjects were excluded, a quadratic fit was better than a linear fit when the RNFL thickness was measured by GDx VCC, whereas StratusOCT-measured RNFL thickness remained best fit with linear regression (Fig. 4) . The linear relationship has been described between visual field sensitivity in the unlogged 1/L scale and other parameters, including pattern-evoked electroretinogram (PERG) amplitude,23 Heidelberg retina tomography (HRT)measured neuroretinal rim area,23 and GDx VCCmeasured RNFL thickness.5 It was proposed that the structure/function change is linear and that it was the decibel scaling in the VF sensitivity that gave the impression of a functional reserve in the curvilinear relationship. In these clinical studies, however, only subjects with glaucoma with mild to moderate defect were included (the mean visual field MD ± SD in the glaucoma group was 4.15 ± 3.05 dB in the study by Garway-Heath et al.23 and 6.9 ± 5.6 dB in that by Schlottmann et al.5 ). Examining a full spectrum of glaucomatous optic neuropathy, from glaucoma suspect to end-stage glaucoma, would be essential to reveal the complete profile of the structure/function relationship. (In the present study, the mean visual field MD ± SD in the glaucoma group was 11.1 ± 7.74 dB; the range was 2.36 to 30.00 dB). The better fitting with a quadratic model in GDx VCC RNFL measurement among the glaucoma and glaucoma-suspect groups (Fig. 4A) suggested this relationship may vary with the composition of the study group. On the other hand, though first-order inverse regression best described the relationship between StratusOCT-measured RNFL and AGIS/CIGTS VF scores, it was the linear model that best fit in the relationships when RNFL thickness was measured by GDx VCC (Figs. 5 6) . Taking all these results into consideration, we concluded that the pattern of the structure/function relationship in glaucoma not only varied with the scale/method of expression in visual sensitivity, it depended on the choice of the RNFL imaging devices and the composition of the study groups.
The discrepancy in the regression profiles between OCT and SLP measurements in describing the relationship between RNFL thickness and VF sensitivity may stem from the intrinsic differences in the imaging principles of OCT and SLP, thus leading to the different estimations of RNFL thickness. The measurement of RNFL thickness in OCT is based on the differential reflectivity signals of the RNFL obtained using a scanning interferometer. It is calculated by measuring the difference between the vitreoretinal interface and the posterior border of the nerve fiber layer (which is identified with a predefined signal threshold algorithm). In contrast to OCT, SLP estimated the RNFL thickness by measuring the retardation of polarized scanning beams generated as a result of the form birefringence of the nerve fiber layer. It has been validated that good correlation was found between retardation and the histopathologic measurement of RNFL thickness, and a constant proportion was adopted for calculation of RNFL thickness in micrometer in relation to the measured retardation in nanometer.24 Yet it was recently reported that the birefringence retardation detected by SLP actually varied with the position around the optic nerve head.25 As a result, using a constant proportion for calculation of RNFL thickness may not reflect the true thickness in some sectors around the optic nerve head. Furthermore, changes in the SLP measurements could result from the changes in either RNFL thickness or RNFL birefringence. Morphologic alteration like gliosis, partial loss of organelles, or shrinkage of ganglion cells might lead to changes in the birefringence pattern before the irreversible loss of axons.26 27 Therefore, variations in the regression profiles of the structure/function relationship may signify the difference in the inherent measuring capacity between OCT and SLP. In the present study, the OCT RNFL measurements attained higher associations with visual function in all the respective regression analyses compared with the GDx VCC RNFL measurements. In addition, plotting of the OCT RNFL measurements against MD (dB) and the unlogged 1/L scale demonstrated consistent curvilinear and linear relationships, respectively (Figs. 1B 2B 3B 4B) , which is in greater agreement with the results from the histologic study in experimental glaucoma investigating the association between ganglion cell loss and decrease in VF sensitivity.28 29 It was reported that the relationship between sensitivity loss (in decibels) and ganglion cell loss (in percentages) was curvilinear and that it changed to a linear function when both VF sensitivity and ganglion cell loss were expressed in decibels.28 29 Collectively, we believe the use of StratusOCT may provide a better approach to understanding the structure/function relationship in glaucoma than the use of GDx VCC.
One of the most controversial risk factors for predicting progression is the severity of the baseline VF. Despite a number of prospective studies performed to investigate whether the initial VF had any predictive power for subsequent progression, different conclusions were reached. In the CIGTS, it was found that patients with higher VF scores at baseline had higher VF scores during follow-up.30 This was repeated by the results in the Early Manifest Glaucoma Trial (EMGT), which reported an increased risk for progression in patients with MD worse than 4 dB.9 Conversely, the AGIS found that lesser VF defect increased the risk for additional VF loss.31 Interestingly, using the same AGIS VF scoring system, Chen and Park10 showed that an increased initial AGIS score was associated with progression. Although most studies supported the finding that increased severity of initial visual field was associated with further visual field worsening,9 10 11 12 30 32 others showed no association,6 7 8 and a few arrived at the same conclusions reported in the AGIS.33 34 In the present study, regression analysis revealed that first-order inverse regression best described the relationship between the AGIS/CIGTS VF scores and StratusOCT-measured RNFL thickness (Figs. 5B 6B) . Assuming the rate of loss of the retinal ganglion cell and its nerve fiber in glaucoma is constant during progression, our results provide an explanation for the findings in those longitudinal studies showing that increased severity of VF defect (documented by the AGIS/CIGTS VF scores or MD) at baseline indicates higher risk for progression. For the same degree of structural damage (reduction in RNFL thickness), the decrease in VF score or decibel is more dramatic in the advanced stages of the disease than in the early stages. For example, based on the regression functions in the study (Figs. 5B 6B) , a reduction in RNFL thickness from 80 to 70 µm (in the early stage) would lead to a corresponding increase in the AGIS score of 2.0 and in the CIGTS score of 2.1. When the RNFL thickness is reduced from 50 to 40 µm (in the advanced stage), the corresponding increase in the AGIS and the CIGTS scores would be 5.7 and 5.9, respectively. Therefore, it is easier to detect the change in VF progression in VF scores during the moderate/advanced stages of disease, leaving the impression of increased risk for progression when the baseline VF scores are high. Therefore, the current AGIS and CIGTS scoring systems are considered less sensitive for detecting progression in the early stages of glaucoma because the steps for progression (an increase in score of 4 in the AGIS and 3 in the CIGTS) is defined independently of the disease stage. Consistent with earlier investigations, we demonstrated the CIGTS VF scores were systematically higher than the AGIS VF scores, and both scorings were highly correlated with each other and with the VF MD.35 36 Our results are also consistent with the findings of longitudinal studies showing that the CIGTS scoring method leads to higher rates of detection of disease progression than the AGIS scoring system.35 37
StratusOCT and GDx VCC are the latest versions of commercially available nerve fiber analyzers, but longitudinal data are still lacking regarding the change in RNFL during glaucoma progression. It is unknown whether the rate of loss of RNFL is constant at different stages of glaucoma; in addition, the temporal relationship in structure/function change is not completely understood. Knowing these factors would have great impact on the interpretation on the structure/function regression profile. A prospective follow-up study is essential to cross-validate the structure/function regression models.
| Conclusions |
|---|
|
|
|---|
| Footnotes |
|---|
Disclosure: C.K. Leung, None; K.K. Chong, None; W. Chan, None; C.K. Yiu, None; M. Tso, None; J. Woo, None; M.-K. Tsang, None; K. Tse, None; W. Yung, 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: Christopher Kai-shun Leung, Department of Ophthalmology, Caritas Medical Centre, 111 Wing Hong Street, Sham Shui Po, Hong Kong, Peoples Republic of China; tlims00{at}hotmail.com.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. Kanamori, M. Naka, A. Nagai-Kusuhara, Y. Yamada, M. Nakamura, and A. Negi Regional Relationship Between Retinal Nerve Fiber Layer Thickness and Corresponding Visual Field Sensitivity in Glaucomatous Eyes Arch Ophthalmol, November 1, 2008; 126(11): 1500 - 1506. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Jehle, K. Wingert, C. Dimitriu, W. Meschede, J. Lasseck, M. Bach, and W. A. Lagreze Quantification of Ischemic Damage in the Rat Retina: A Comparative Study Using Evoked Potentials, Electroretinography, and Histology Invest. Ophthalmol. Vis. Sci., March 1, 2008; 49(3): 1056 - 1064. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Cheng, M. Laron, J. S. Schiffman, R. A. Tang, and L. J. Frishman The Relationship between Visual Field and Retinal Nerve Fiber Layer Measurements in Patients with Multiple Sclerosis Invest. Ophthalmol. Vis. Sci., December 1, 2007; 48(12): 5798 - 5805. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. K.-s. Leung, F. A. Medeiros, L. M. Zangwill, P. A. Sample, C. Bowd, D. Ng, C. Y. L. Cheung, D. S. C. Lam, and R. N. Weinreb American Chinese Glaucoma Imaging Study: A Comparison of the Optic Disc and Retinal Nerve Fiber Layer in Detecting Glaucomatous Damage Invest. Ophthalmol. Vis. Sci., June 1, 2007; 48(6): 2644 - 2652. [Abstract] [Full Text] [PDF] |
||||
![]() |
J.-C. Tsai, P.-W. Lin, M.-C. Teng, and I.-C. Lai Longitudinal Changes in Retinal Nerve Fiber Layer Thickness after Acute Primary Angle Closure Measured with Optical Coherence Tomography Invest. Ophthalmol. Vis. Sci., April 1, 2007; 48(4): 1659 - 1664. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Ajtony, Z. Balla, S. Somoskeoy, and B. Kovacs Relationship between Visual Field Sensitivity and Retinal Nerve Fiber Layer Thickness as Measured by Optical Coherence Tomography Invest. Ophthalmol. Vis. Sci., January 1, 2007; 48(1): 258 - 263. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Bowd, L. M. Zangwill, F. A. Medeiros, I. M. Tavares, E. M. Hoffmann, R. R. Bourne, P. A. Sample, and R. N. Weinreb Structure-function relationships using confocal scanning laser ophthalmoscopy, optical coherence tomography, and scanning laser polarimetry. Invest. Ophthalmol. Vis. Sci., July 1, 2006; 47(7): 2889 - 2895. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK |