|
|
||||||||
1From the Center for Visual Science, the 2Department of Brain and Cognitive Sciences, and the 4Institute of Optics, University of Rochester, Rochester, New York; and the 5University of North Carolina School of Medicine, Chapel Hill, North Carolina.
| Abstract |
|---|
|
|
|---|
METHODS. Frequency-of-seeing curves were measured with 0.75' and 7.5' spots. On each trial, from 0 to 4 stimuli were randomly presented at any of eight equally spaced loci 0.5° from fixation. By correcting the aberrations of the eye, adaptive optics produced retinal images of the 0.75' spot that were 3.0 µm wide at half height, small enough to be almost entirely confined within the typical cone diameter at this eccentricity. Data were collected from a patient with deuteranopia (AOS1) whose retina, imaged with adaptive optics, suggested that
30% of his cones were missing or abnormal. Patients with protanomalous trichromacy (1 subject), deuteranopia (1 subject), and trichromacy (5 subjects) served as controls (all had normal cone density and complete cone mosaics). Psychophysical results were modeled by a Monte Carlo simulation incorporating measured properties of the cone mosaic.
RESULTS. Frequency-of-seeing curves for AOS1 obtained with 0.75' spots showed lower asymptote, slope, and sensitivity than for controls. The 7.5' results showed that these differences were the result of the small spot size, which on some trials was confined mostly to the locus of the putatively missing cones. A two-parameter model satisfactorily described the data and was highly sensitive to the proportion of missing cones simulated.
CONCLUSIONS. Adaptive-optics microperimetry is a powerful psychophysical test for assessing the loss of neural elements, even in retinas that appear otherwise normal in standard clinical tests. This technique may prove useful in estimating the proportion of missing cones in different patients and in detecting other visual losses such as those associated with glaucoma.
The person whose retina we tested (AOS1) represented a subclass of persons with dichromacy who carry genetic mutations that disrupt the function of one of the cone photopigments.12 Previous results with high-resolution retinal imaging13 have shown that the retina of AOS1 contained lacunae, equal in size to one or more cones, that the authors hypothesized represented the loci at which cones expressing the mutant gene have been severely damaged or lost. A micrograph of his retina is shown in the upper part of Figure 1 ; for comparison, the retina of a healthy trichromat (AOS2) is shown below it. Note the relative sparseness of visible cones in the upper micrograph. Although the hypothesized missing cones occupy almost a third of the cone mosaic, ophthalmic examination revealed no abnormalities aside from the dichromacy. For example, visual acuity was 20/16 (0.10 logMAR), standard visual fields were normal, and funduscopic examination results were unremarkable. Test stimuli used here were small enough to fall largely within the lacunae observed photographically and thereby made it possible to detect the associated microscotomas psychophysically. Use of large spots, approximating the size of those in clinical use, failed to produce the signs we attributed to the missing cones.
|
| Methods |
|---|
|
|
|---|
Stimulus Presentation
One drop each of phenylephrine hydrochloride (2.5%) and tropicamide (1%) dilated the subjects pupils and suspended accommodation in their right eyes. Their heads were stabilized by a bite bar containing their dental impression. Stimuli were 0.75' or 7.5' spots produced by a digital lightprocessing (DLP) video projector (MP1600; Compaq; Hewlett Packard, Mountain View, CA)15 at any of eight equally spaced locations 0.5° from a dim, 2' fixation target (Fig. 2) . The light passed through an interference filter with a 10-nm pass band centered at 550 nm, a wavelength to which M and L cones are equally sensitive. At maximum, the test spot illuminated the retina at 5.06 log troland. The projector also provided a steady background field of 3.08 log troland (its minimum output) covering a circular region slightly greater than 1°. Sensitivity to the 7.5' spot was so much greater than that to the 0.75' spot that a 1.11 log filter was required in the projector beam to present spots dim enough to measure the low end of the frequency-of-seeing curve. This brought the maximum log troland value down to 3.95. Because the filter also reduced the background field, we added a steady, 550-nm background field that brought the total background field for the 7.5' spots back up to 3.44 log troland, comparable to that for the 0.75' spots.
|
|
|
Because the time available for testing AOS1 was limited, we increased the rate of data acquisition by simultaneously presenting 0, 1, 2, 3, or 4 stimuli with equal frequency as is sometimes done in clinical perimetry.16 17 18 Although one can further increase efficiency by presenting four stimuli on a higher proportion of trials, we chose not to so as to avoid the response biases that might result. The subjects task was to report, by a key-press, the number of stimuli presented. This increase in the number of stimuli and alternative responses only slightly delays the response19 20 and has no noticeable effect on the reliability of detection19 21 22 23 (our results also show no relationship between detection and number of stimuli). The subjects response initiated the next trial.
The stimuli described were presented through an adaptive optical system that minimized optical aberrations.24 25 Measured point-spread functions with and without adaptive optics are shown in Figure 3 . Although cone sizes varied over the areas tested and among subjects, as a rough estimate, without adaptive optics approximately 9% of the incident light in a 0.75' test spot fell within the geometric boundaries of a single cone, whereas with adaptive optics approximately 55% did, so that adaptive optics increased the amount of light in this area approximately sixfold. Because of the Gaussian shape of the cone entrance aperture, this optical correction increased the relative amount of light that actually entered the cone even more.
Preliminary observations established the range spanned by each observers frequency-of-seeing curve, and from five to eight points on the curve, spaced at 0.1 intervals on a log troland scale, were tested to describe the curve. A block of trials consisted of 50 trials at each stimulus magnitude, in random order, and each observer had four such blocks. Measurement and correction of the aberrations of each subjects eye was performed between each block of trials, instead of during the blocks, to prevent the superluminescent diode used for wavefront sensing from affecting the subjects sensitivity. Although the measurement took less than 800 ms, more than 1 minute was required for the subjects eye to recover from exposure to the light from the diode.
Model
To determine whether the differences between the 0.75' data of AOS1 and of the control subjects could be explained on the basis of the evident differences in their retinal mosaic, we simulated the experiment with a model (implemented in MATLAB; MathWorks, Natick, MA) that incorporates the retinal mosaic and optical point-spread functions of AOS1 and that of one of the trichromatic subjects (AOS2).
The model consisted of a Monte Carlo simulation of 10,000 stimuli of varying energy, presented to a given subjects retinal cone mosaic. For each stimulus presentation, the model used a set of detection rules to determine whether the stimulus was seen. The output of the model was a frequency-of-seeing curve.
According to the model, a given stimulus was detected if the visual excitation exceeded the response criterion for that flash. Visual excitation was calculated by summing the products of each point of the retinal image of the spot with the corresponding point of the sensitivity profile of the photoreceptor mosaic (i.e., by finding the inner product). The retinal image of the spot was found by convolving the geometric image of the stimulus with the subjects own point spread function, measured at an eccentricity of 0.5°, after the aberrations of that subjects eye had been optically compensated.24 25
The sensitivity profile of the photoreceptor mosaic was obtained by convolving the optical aperture of the cones with a set of points representing locations of the observed centers of the cones in a patch of retina approximately 0.75° square that was imaged in vivo at 0.5° eccentricity by an adaptive optics ophthalmoscope.26 Individual cones were assigned a spectral sensitivity based on the relative numbers of short-wavelength-sensitive (S), M, and L cones in the retina. Because these vary across healthy subjects, we used previous densitometry measurements with an adaptive optics ophthalmoscope13 24 27 28 to determine the appropriate ratio of S-, M-, and L-cone assignments to use for each subject. In addition, the optical aperture of each cone was represented by a two-dimensional Gaussian profile with a width at half height of 0.615 times the intercone distance, on the basis of interferometric estimates.29 30 31 32 We truncated the Gaussian at the physical boundary of the cone (meaning that light falling on the cone mosaic outside the cones physical aperture would not be absorbed by photopigment within that particular cone), but doing so had a negligible effect on the results (data not shown).
Because of fixational eye movements, the absolute location of a stimulus on the retina varied from trial to trial. Therefore, we incorporated fixational instability in the model to determine the location of the retinal mosaic on which the stimulus fell on each trial. The position of the retinal mosaic with respect to the retinal image of the spot was sampled randomly from a probability density function representing the statistics of eye position measurements. The two-dimensional distribution of fixations of AOS2 and four other subjects were observed by using high-resolution imaging with adaptive optics to record the absolute position of the retina while fixating.24 The fixational patterns of AOS2 so obtained were used in the simulations of his psychophysical performance, whereas the mean of all five subjects was used in the simulation the performance of AOS1. Consistent with previous reports using a different methodology, in both cases the fixational patterns were elongated along the horizontal axis and had an SD of approximately 5'.
The model therefore yields on each trial an estimate of the relative amount of light absorbed by each cone in the vicinity of the test flash. Although previous models of the detection of small spots by the cone system33 34 35 36 37 38 39 40 41 42 have assumed that each cone has a threshold, in the sense that stimulation of a cone below a fixed value makes no contribution to detection by the subject, we found such a threshold unnecessary; excluding it simplifies the model while improving the fit slightly (though not reliably). We assumed instead that the total number of photopigment molecules activated (visual excitation) in all the cones sums linearly to determine whether the flash was seen: if the total number of activations exceeded a given value, the response criterion, we assumed that the subject detected the stimulus. (Nonlinear summation of cone excitation, such as root-mean-square summation,43 44 45 46 yields frequency-of-seeing curves with too great a slope.) To obtain frequency-of-seeing curves, we simply counted the proportion of trials on which the visual excitation exceeded the response criterion at each of the stimulus intensities.
The model thus has two free parameters that were separately adjusted to fit each subjects data: the observers response criterion, just mentioned, and sensitivity, the proportion of quanta in the stimulus that ultimately activates a photopigment molecule. Sensitivity simply relates such visual excitation to the amount of light absorbed, and changing its value translates the frequency-of-seeing curve along the log intensity axis. The response criterion is the visual excitation required for the subject to respond to the stimulus; changing its value translates the frequency-of-seeing curve along the absolute intensity axis and varies the slope and shape of the frequency-of-seeing curve when plotted on semilog axes. Thus, in the model, two parameters serve where three were required for the curves in Figure 4 .
In our simulations, variation in the locations of the spots relative to the cones produced by eye movements caused variation in visual excitation from trial to trial, depending on where the spots fell with respect to the cones and gaps between cones in the mosaic. Increasing the response criterion increased the number of trials in which too little light was absorbed to be detected. This affected both the height and the slope of the frequency-of-seeing curve. Hence, the topography of the individuals retinal mosaic determined the trial-by-trial variability, and the subjects response criterion determined the effect of that variability on the frequency-of-seeing curve.
The variability of light absorbed qualifies as noise, and so the model is based on the theory of signal detection.47 48 49 However, in the experiments using the 7.5' spot, the variability of absorbed light resulting from variation in spot location was negligible. The stochastic noise associated with quantal absorptions was also small, approximately 0.3% of the mean (SD). Hence, the trial-by-trial variability of the response has to be attributed to some other source, such as noise intrinsic to the visual system. For AOS2 this inferred noise was 22 times smaller than that for the smaller spot, and for AOS1 it was 115 times smaller. The contribution of such visual noise to detection of the small spot was therefore small relative to that associated with stimulus location. Nevertheless, for the sake of consistency, we included this noise in the model for detection of the small spot, assuming that the amount of intrinsic noise involved in detecting the small spot was the same as that for the large spot. This was incorporated into the model by adding a random sample from a Gaussian distribution to the total visual excitation on each trial.
The model is of course somewhat simplified. For example, although the intercone distances in our model faithfully reflected those of the individual retina tested, we assumed a fixed cone size. Similarly, the stimuli actually fell at eight locations, but our model represented the cone mosaic at only one location: there can be small differences in cone density on different meridians,50 51 and the effects of any such asymmetries could be exaggerated by the tendency of observers to fixate off the anatomic center of the fovea.24
The model also omitted absorption of light by the retinal blood vessels. Approximately 90% of the retina is covered with a capillary bed that absorbs some 40% of incident light.52 Placing the stimuli 0.5° from fixation, an area that in most persons is devoid of retinal vasculature, mitigates this problem. However the size of this avascular zone varies widely, and its radius in healthy persons can be as small as 0.3°.53 Indeed, retinal imaging reveals that subject AOS6 in this experiment has extensive invasion of the central fovea by retinal vessels (data not shown). Therefore, one would expect the slope of his frequency-of-seeing curve to be shallower than those of the other subjects and also shallower than predicted by the model. That this subject does have the shallowest frequency-of-seeing curve is evident in Table 1 and Figure 4 , and it is significantly shallower than that generated by the model for his retina (data not shown).
|
| Results |
|---|
|
|
|---|
![]() | (1) |
![]() | (2) |
|
|
The larger test spot produced a threshold difference between AOS1 and the control subjects that was 40% greater (P < 0.01) than for the smaller spot. Hence, AOS1 shows less spatial summation than the control subjects, possibly because of lack of a complete cone mosaic, but we cannot rule out other abnormalities of his visual system that are independent of this loss of cones.
The data from the other color-deficient control subjects (AOS7 and AOS8) do not differ reliably from those of the trichromatic subjects in any way. When AOS7 and AOS8 are compared with AOS1, only the asymptotes for the 0.75' data and the means for the 7.5' data differ reliably (t = 17.3 and 29.2, respectively; P < 0.05 in both cases); the decrease in the number of reliable differences may be due to the fact that the differences between AOS1 and two subjects (the two dichromats) must be greater to be statistically reliable than differences between AOS1 and seven subjects (all the control subjects).
These results show that some of the sensitivity differences between AOS1 and the control subjects show up only when small spots are used. Next we examined whether the difference could be accounted for on the basis of the differences in retinal topography.
Results of the Model
Figure 5 shows that the differences between the frequency-of-seeing curves of AOS1 and AOS2 for the small spot could be quantitatively accounted for by a model that incorporated the specific cone topographies of the respective retinas. The only differences between the models that produced these two curves were the respective subjects cone topographies, point-spread functions, and pattern of eye movements and the associated values of the two free parameters. The subjects point-spread functions and eye movements could be interchanged without eliminating the difference between the curves; hence, the difference in the curves was attributed to the difference in cone topographies.
The next question was whether the particular difference between retinal topographies that accounted for the difference in results lay in the number of microscotomas. Figure 6a shows the effects of varying the number of microscotomas: this family of curves resulted from running the model on the representation of the retina of AOS1 after removing varying numbers of cones (i.e., in the top two curves, missing cones have been replaced for the purposes of these tests). Free parameters are identical for all curves. Finer spacing of the curves near 72% remaining cones showed that an increase to 78% or a decrease to 67.7% produced a curve that differed reliably from the data (0.05 level, F test). This corresponded to a maximal difference in frequency-of-seeing of 0.17 log unit. According to the model, then, the proportion of the full complement of cones remaining in the retina of AOS1 was within this range; the value estimated anatomically was near the center of this range.13
For comparison, Figure 6b shows the corresponding family of curves for the retina of AOS2. These families of curves show that the differences between the frequency-of-seeing curves of AOS1 and AOS2 in Figure 4 were quantitatively attributable to the microscotomas in the retina of AOS1.
The two families of curves in Figures 6a and 6b differed from one another. In other words, even if AOS1 had a full complement of cones, his frequency-of-seeing curve would not be the same as that of AOS2. Such differences limit the accuracy with which one can use such curves to estimate the proportion of cones missing. To get an idea of how greatly these differences affect attempts to estimate the proportion of missing cones from the curves derived from psychophysical observations, we have plotted in Figure 6d the proportion of missing cones as a function of the maximum of the frequency-of-seeing curves for the two subjects. The greatest difference between corresponding curves is 0.05 log unit. A greater number of retinas will be needed to estimate the distribution of differences among subjects, but these preliminary data suggest that one can estimate the proportion of missing cones from the frequency-of-seeing curves to within approximately 5%.
| Discussion |
|---|
|
|
|---|
The next question is how small they have to be. To evaluate this, we ran our model simulation using the natural point-spread function of AOS1, without correction of the higher-order aberrations, to estimate AOS1s performance without adaptive optics. Results are shown as the dashed curves in Figure 6c , along with the simulations obtained when adaptive optics were used (the solid curves, taken from Fig. 6a ). The right curve of each pair shows the results using AOS1s actual retinal mosaic, and the other was obtained by filling the lacunae with cones, representing a normal retina with a 100% complement of cones. The differences between the right and left curves in each pair are attributed to the missing cones, and it is only on the basis of these differences that the gaps left by the missing cones can be detected psychophysically.
To quantify these results and to compare them with the data from the control subjects, we fitted a cumulative Gaussian curve to the curves in Figure 6c (nearly a perfect fit), as was done with the psychophysical data. These results show that the main benefit of adaptive optics shows up in the asymptote of the frequency-of-seeing curve. Without adaptive optics, the difference between the asymptotes of the curves for 72% and 100% mosaics was 0.027, less than the variability of asymptotes among the control subjects, as represented by the SD of the asymptotes in Figure 4 . However, when adaptive optics were used, the difference between asymptotes was 0.17, more than six times greater than without adaptive optics and more than four times the SD for the control subjects. However, the means and the spread parameters (
) of the curves were comparable regardless of whether adaptive optics were used. Thus, although it is possible that without adaptive optics the loss of 30% of the cones could have been detected on the basis of the means and spread parameters, adaptive optics does offer an enormous increase in sensitivity and precision (discussed in "Adaptive-Optics Microperimetry: Potential as a Research Tool") through its effect on the height of the frequency-of-seeing curves.
Moreover, the hope of detecting such losses without adaptive optics is damped by the results of Geller and Sieving.37 Using a test spot that subtended 1.125' before blurring by the eye, they observed decreased sensitivity and a decreased slope of the frequency-of-seeing curves with these patients, as we observed with AOS1, but no measurable effect on the asymptote. However, quantitatively, the effects on the slopes were much smaller than those reported. The slopes of curves from two of their patients, who were estimated to have lost 85% to 88% of their cones, corresponded to a 15% loss in our study (Figure 6) . In the other two patients, an estimated loss of 85% to 88% of the cones yielded slopes corresponding to 28%, and 55% losses in our study (Figure 6) . One cannot know whether the slope of the frequency-of-seeing curve of AOS1 would have been abnormal in their experiment, but if the curves of two patients who have lost 85% of their cones are barely abnormal (i.e., corresponding to a 15% loss), then the abnormality of someone with only a 30% loss is unlikely to be detected. The smaller effects observed by Geller and Sieving37 might have resulted from optical aberrations they were unable to correct.
Adaptive-Optics Microperimetry: Potential as a Research Tool
The technique outlined here serves not only to detect the loss of cone function but holds promise as a means of measuring the amount of loss. Given that the shape of the frequency-of-seeing curve depends on the proportion of cones lost (Fig. 6) , it can be used to estimate the magnitude of functional loss, perhaps with a precision within 5%. This might be valuable in trying to monitor progressive loss of cone function in patients with cone degeneration, especially in those in whom functional loss might precede visible signs of loss in retinal images. Note that by testing with a wavelength to which both M- and L-cones are equally sensitive, there is no difference between them as far as this test is concerned. Any cone loss, whether specific to a particular cone class or randomly distributed, should reveal itself to these test stimuli.
This technique may also hold promise as a means of detecting small scotomas caused by other anomalies, such as the loss of ganglion cells that occurs in glaucoma. We speculate here that most techniques for detecting missing cones are hampered by redundancy in the representation of the stimuli in the receptor layer. The detection of missing ganglion cells is hindered not only by the same source of redundancy but also by redundancy associated with overlapping receptive fields. If one ganglion cell is lost, others continue to respond. The total signal transmitted may be smaller, and the loss of a particular ganglion cell may change the properties of the signal transmitted, but the loss of a ganglion cell fails to produce a complete scotoma, as the loss of a single cone does. Nevertheless, the stimuli used here reduce the redundancy that limits the effectiveness of such tests and are likely to increase sensitivity to the loss of ganglion cell activity.
| Acknowledgements |
|---|
| Footnotes |
|---|
Presented, in part, at the 2005 OSA Fall Vision Meeting in Tuscon, Arizona, October 2123.
Supported by National Eye Institute Grants R01-EY04367, P30-EY01319, and F32-EY014749; and in part by the National Science Foundation Science and by the Technology Center for Adaptive Optics, managed by the University of California at Santa Cruz under Cooperative Agreement AST-9876783.
Submitted for publication September 7, 2005; revised February 23 and April 3, 2006; accepted July 3, 2006.
Disclosure: W. Makous, None; J. Carroll, None; J.I. Wolfing, None; J. Lin, None; N. Christie, None; D.R. Williams, P, C
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: Walter Makous, Center for Visual Science, University of Rochester, Box 270268, Rochester, NY 14627-0270; walt{at}cvs.rochester.edu.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M Nilsson, G von Wendt, P Wanger, and L Martin Early detection of macular changes in patients with diabetes using Rarebit Fovea Test and optical coherence tomography Br. J. Ophthalmol., December 1, 2007; 91(12): 1596 - 1598. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |