|
|
||||||||
1From the University of Pittsburgh Medical Center Eye Center, Eye and Ear Institute, Ophthalmology and Visual Science Research Center, Department of Ophthalmology, University of Pittsburgh School of Medicine, Pittsburgh, Pennsylvania; the 3Department of Electrical Engineering and Computer Science and Research Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, Massachusetts; the 4Institute of Physics, Nicolaus Copernicus University, Torun, Poland; and the 5New England Eye Center, Tufts-New England Medical Center, Tufts University School of Medicine, Boston, Massachusetts.
| Abstract |
|---|
|
|
|---|
METHODS. Three-dimensional, Fourier/spectral domain optical coherence tomography (OCT) data were obtained as raster scan data (512 x 180 axial scans in a 6 x 6-mm region centered on the optic nerve head [ONH]) with high-speed, ultrahigh-resolution OCT (hsUHR-OCT) from 12 healthy subjects. RNFL thickness was measured on this three-dimensional data set with an in-house software program. The disc margin was defined subjectively in each image and RNFL thickness profiles relative to distance from the disc center were computed for quadrants and clock hours. A mixed-effects model was used to characterize the slope of the profiles.
RESULTS. Thickness profiles in the superior, inferior, and temporal quadrants showed an initial increase in RNFL thickness, an area of peak thickness, and a linear decrease as radial distance from the disc center increased. The nasal quadrant showed a constant linear decay without the initial RNFL thickening. A mixed-effects model showed that the slopes of the inferior, superior, and nasal quadrants differed significantly from the temporal slope (P = 0.0012, P = 0.0003, and P = 0.0004, respectively).
CONCLUSIONS. RNFL thickness is generally inversely related to the distance from the ONH center in the peripapillary region of healthy subjects, as determined by hsUHR-OCT. However, several areas showed an initial increase in RNFL, followed by a peak and a gradual decrease.
Histologic studies of human and primate eyes have sampled the RNFL at various locations on the retina, both near the ONH and peripherally, and have shown that the convergence of ganglion cell axons from the retinal periphery toward the optic disc gives rise to an increasing RNFL thickness as the nerve head is approached.8 9 10 11 12 Skaf et al.13 recently determined four RNFL thickness profiles (inferonasal, superonasal, superotemporal, and inferotemporal) in healthy eyes using conventional optical coherence tomography (OCT) with a limited number of sampling points.13
Improvements in OCT technology have recently been introduced.14 15 16 17 Among these new iterations is high-speed, ultrahigh-resolution OCT (hsUHR-OCT), which uses Fourier/spectral domain detection to provide increased resolution and scanning speed compared with conventional time-domain OCT.14 18 Cross-sectional retinal images with an axial resolution up to five times higher than conventional OCT and imaging speeds 60 times faster than conventional OCT have been acquired in vivo.19 20 This increase in resolution and scanning speed permits high-density raster scanning of retinal tissue while minimizing eye motion artifacts. It is then possible to detect and segment the RNFL in each raster OCT image and use these data to construct a detailed RNFL thickness map.21
The purpose of this study was to use hsUHR-OCT RNFL thickness information to analyze the peripapillary RNFL architecture by creating a profile of RNFL thickness as a function of increasing distance from the ONH.
| Methods |
|---|
|
|
|---|
Inclusion criteria were best corrected visual acuity of 20/40 or better, refractive error between ± 6.0 D, no media opacity, normal clinical ocular examination with no evidence of peripapillary atrophy and reliable and normal 24-2 standard Swedish interactive thresholding algorithm (SITA) perimetry. A reliable SITA was defined by less than 30% fixation losses and false-positive and -negative responses. Normal visual field was defined as glaucoma hemifield test (GHT) result within normal limits. If both eyes were eligible, one eye was randomly selected for the study.
Instrument
All subjects had hsUHR-OCT raster scanning (512 x 180 axial scans in a 6 x 6-mm area, 2560 A-scans/mm2) of the ONH region without dilating the pupil. A schematic diagram and description of this prototype hsUHR-OCT instrument has been published.19 The device used in this study, however, had a superluminescent diode light source (Broadlighter; Superlum Diodes, Ltd., Moscow, Russia) with 100 nm bandwidth (full width at half maximum) and a center wavelength of 840 nm, corresponding to an axial resolution of
3.5 µm. The system had an A-scan acquisition rate of 24 kHz, resulting in an acquisition time of 3.84 seconds for each raster three-dimensional data set. An hsUHR-OCT en face fundus image was created for each raster scan by using custom software. The software created this image by summing all pixel intensity values along individual A-scans. The resultant sum along each A-scan was the intensity value used for the pixel corresponding to that A-scans location within the two-dimensional en face hsUHR-OCT fundus image (Fig. 1A) . hsUHR-OCT fundus images were used for evaluation of eye motion during the scans and defining disc margins. Multiple raster scans centered on the ONH were acquired for each subject and the best-quality image was subjectively chosen and used for subsequent analysis. Criteria for acceptable hsUHR-OCT fundus images included no large eye movements, defined as an abrupt shift in a large retinal vessel that completely disconnected the vessel in greater than three consecutive frames. In addition, we required that there be no black bands (caused by blinking during acquisition) and that there be consistent signal intensity level across the entire scan.
|
To ensure that a sufficient number of data points were present for the profile analysis, clock hour segments were excluded if the average intrasubject profile did not extend at least 200 pixels (
2.4 mm) from the center of the ONH, which could occur if a disc was off-center during scanning. Mean RNFL thickness profiles for the entire study population, at each quadrant and clock hour, were averaged, and the region in which all qualified subjects had measurements was plotted (average intersubject profile). This method ensured that all profile lines were equally represented by all available subjects.
To investigate the decaying portion of RNFL profile, thickness measurements for each subject for each quadrant were smoothed using a mean filter. The location of the distal edge of the most distal peak of RNFL thickness was subjectively identified along the profile. Thickness measurements distal to this location were used for this analysis. Thus, the decay portion of filtered RNFL data could start at a different location for each subject. The measurements were averaged for each quadrant and the region in which all qualified subjects had measurements was plotted. The slope of the RNFL thickness profile was computed for each subject and each sector from these data. A linear mixed-effects model was fitted to the profile slope as a function of radius from the disc center, taking into account clustering within subjects.22
To investigate the variation along the RNFL thickness profile among the entire study population, the variance at each radial point (±50 µm) was plotted for each quadrants mean-filtered data. An area of stable measurements with low variability among subjects has a variance approaching zero.
All the analyses were also conducted when the RNFL profiles were aligned, such that they all began at the disc margin without taking into account distance from the disc center. This alignment was made to ensure that variability in disc size within the study group did not affect the tissue profiles.
| Results |
|---|
|
|
|---|
RNFL thickness profile graphs of all quadrants and clock hours are shown in Figure 2 . Thickness profiles in the superior, inferior, and temporal quadrants showed an initial increase in RNFL thickness, an area of peak thickness, and finally a decrease as radial distance from the center of the disc increased up to 2.5, 2.8, and 2.5 mm from the disc center, respectively. The nasal quadrant showed a consistent decrease in RNFL thickness up to a distance of 2.3 mm from the disc center. Three clock hour segments (2, 3, and 4 oclock) did not exhibit an area of RNFL thickening and instead demonstrated a constant decrease in RNFL thickness only with increasing radius. The location of the scanning beam of the commercially available StratusOCT system is indicated on each profile (1.7 mm radius) and the thickness measurements are listed in Table 1 . The overall mean RNFL thickness using StratusOCT was 103 and 153 µm using hsUHR-OCT, a difference of 50 µm between the two devices.
|
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
The current clinical standard for obtaining RNFL thickness measurements with StratusOCT is a peripapillary circular scan centered on the ONH with a radius of 1.7 mm. Measurements at this location have been shown to be the most reproducible compared with two other scan circles of different diameters.23 24 RNFL thickness measurements at this diameter have also shown high intrasession and intersession reproducibility in healthy25 and glaucomatous eyes in the commercially available device.26 However, the OCT scanning position was originally arbitrarily chosen to avoid intersecting tissue within the ONH margin in large discs and to avoid areas with peripapillary atrophy, but simultaneously to be close enough to the disc to allow a dense sampling covering the entire distribution of RNFL. Our data showed that in this cohort of healthy subjects at a distance of 1.7 mm from the disc center, the RNFL was thickest in the inferior, superior, and temporal segments. This is coincidentally the region of low thickness variability in the superior and inferior quadrants (Fig. 5) . Thus, a decrease in the RNFL thickness in these thicker and more stable areas of the RNFL might be easier to detect. Measurements obtained from such locations would be more reproducible than in a constantly decaying region and would also be less prone to registration error. This is consistent with findings regarding reproducibility of OCT circumpapillary scans at various scan diameters.23
Mean RNFL thicknesses, shown in Table 1 , were noticeably higher than those obtained with commercial time-domain OCT, probably because a different segmentation algorithm is used and the axial calibration of the hsUHR-OCT is different from the commercial device.
In the slope analysis, the nasal quadrant showed a slope that was similar to both the superior and inferior quadrants. However, the temporal quadrant showed shallower decay than the other quadrants (Fig. 4 , Table 2 ), perhaps because the temporal RNFL bundle covers a narrower area in the retina and is a more dense structure than the other quadrants.
Findings similar to those reported herein were noted in each of the individual participants profiles and when the RNFL profiles were aligned to the disc margin instead of the disc center (Fig. 3) . Therefore, the pattern observed is not due to an averaging artifact but is instead due to tissue properties. This study was not designed to assess the effect of disc size on the RNFL profile and therefore further investigation is required to evaluate this effect.
Our findings agree with those of Skaf et al.,13 who created a profile of the RNFL using eight concentric circles around the ONH using StratusOCT. However, there are several limitations associated with the methods used in that study. Because the profile was created using only eight radial sampling distances, interpolation was necessary to fill in gaps between measurements. When considering the same retinal area used in our study (36 mm2), the number of thickness measurements in the study done by Skaf et al. was much lower than in our study (eight A-scans/mm2 vs the 2560 A-scans/mm2, respectively), as only nine points at each radial location (inferonasal, superonasal, superotemporal and inferotemporal) were used. In addition, data acquired using concentric circles with a separation of 0.2 mm are prone to eye motion artifacts during scanning with a slower time domain system. Processing of our three-dimensional data allowed virtual OCT fundus images to be created for our study. These images made it possible to exclude images in which large eye movements occurred and ensured there was a more accurate sampling registration.
The mean disc area in this study (4.3 ± 0.7 mm2) was larger than the commonly reported area, probably the result of forcing the disc margin to be defined by a perfect circle, thus overestimating the disc size in the temporal and nasal regions. The axial length range of subjects participating in this study was 22.89 to 25.57 mm. We did not use a magnification correction for the axial length because the range in this study had a negligible effect on the size of the scanning area (Piette S et al. IOVS 2002;43:ARVO E-Abstract 255) In addition, in a previous study in which OCT 2 was used, Bayraktar et al.27 found that in subjects with a wide range of axial lengths, the difference in RNFL thickness measurements before and after adjusting for magnification was less than the variability range inherent to the device. Therefore, for the range observed in our group, the effect was considered to be negligible.
The RNFL border-detection algorithm used in our study is influenced by neighboring measurements, and areas close to the ONH margin may have been affected due the rapid changes in RNFL at that location. Although data within the ONH region were excluded from analysis, the RNFL thickness profile measurement variability noted near the ONH margin and at the periphery in all quadrants and clock hours was most likely due to higher segmentation algorithm variability in these areas (Fig. 6) . This algorithm also tended to overestimate the RNFL thickness when vessel shadows were present. The ONH margin often contained large blood vessels, and both the inner and outer RNFL borders had rapid contour changes (Fig. 6) . This resulted in less reliable segmentation by the software in these specific locations. One other possible limitation of our study is the uneven sampling density of raster scanning since more axial scans were acquired in the horizontal direction (512 axial scans) than in the vertical (180 axial scans).
|
In conclusion, RNFL thickness is generally inversely related to the distance from the ONH center in the peripapillary region of healthy subjects, as determined with hsUHR-OCT. Aside from the nasal segment, all areas show an initial increase in RNFL, followed by a peak and a gradual linear decrease. In the nasal segment a linear decay appears without initial RNFL thickening.
| Footnotes |
|---|
Presented at the annual meeting of the Association for Research in Vision and Ophthalmology, Fort Lauderdale, Florida, May, 2006.
Supported in part by National Eye Institute Grants R01-EY13178-06, R01-EY11289-20, P30-EY08098, and P30-EY13078; National Science Foundation Grants ECS-0119452 and BES-0522845; Air Force Office of Scientific Research Grant FA9550040-1-0011; Medical Free Electron Laser Program Grant FA9550-040-1-0046; The Eye and Ear Foundation, Pittsburgh, PA, and an unrestricted grant from Research to Prevent Blindness, Inc.
Submitted for publication November 29, 2006; revised January 24, 2007; accepted April 13, 2007.
Disclosure: M.L. Gabriele, None; H. Ishikawa, None; G. Wollstein, None; R.A. Bilonick, None; L. Kagemann, None; M. Wojtkowski, None; V.J. Srinivasan, None; J.G. Fujimoto, None; J.S. Duker, None; J.S. Schuman, 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: Gadi Wollstein, UPMC Eye Center, Department of Ophthalmology, University of Pittsburgh School of Medicine, 203 Lothrop Street, Eye and Ear Institute, Suite 834, Pittsburgh, PA 15213; wollsteing{at}upmc.edu.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. N. Menke, P. Knecht, V. Sturm, S. Dabov, and J. Funk Reproducibility of Nerve Fiber Layer Thickness Measurements Using 3D Fourier-Domain OCT Invest. Ophthalmol. Vis. Sci., December 1, 2008; 49(12): 5386 - 5391. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. J. Srinivasan, D. C. Adler, Y. Chen, I. Gorczynska, R. Huber, J. S. Duker, J. S. Schuman, and J. G. Fujimoto Ultrahigh-Speed Optical Coherence Tomography for Three-Dimensional and En Face Imaging of the Retina and Optic Nerve Head Invest. Ophthalmol. Vis. Sci., November 1, 2008; 49(11): 5103 - 5110. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. L. Gabriele, H. Ishikawa, G. Wollstein, R. A. Bilonick, K. A. Townsend, L. Kagemann, M. Wojtkowski, V. J. Srinivasan, J. G. Fujimoto, J. S. Duker, et al. Optical Coherence Tomography Scan Circle Location and Mean Retinal Nerve Fiber Layer Measurement Variability Invest. Ophthalmol. Vis. Sci., June 1, 2008; 49(6): 2315 - 2321. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |