|
|
||||||||
1From the Department of Ophthalmology and Visual Sciences, University of Iowa Hospitals and Clinics, Iowa City, Iowa; the 2Department of Veterans Affairs, Iowa City VA Medical Center, Iowa City, Iowa; the 3Department of Neurosurgery, SUNY Health Sciences Center, Syracuse, New York; 4Visionquest Biomedical Inc., Albuquerque, New Mexico; the 5Department of Biostatistics, College of Public Health, University of Iowa, Iowa City, Iowa; and the 6Department of Ophthalmology and Visual Science, University of Chicago, Chicago, Illinois.
| Abstract |
|---|
|
|
|---|
METHODS. Five normal human subjects were studied. A modified fundus camera was used to image changes in retinal reflectance of 780-nm near-infrared light imaged onto a 12-bit charge-coupled device (CCD) camera in response to a green (540 nm) visual stimulus. During 60 seconds of recording (frame rate, 3 Hz) 10 cycles were recorded, during each of which 3 seconds of blank and then 3 seconds of either vertical bar or blank stimulus was projected. The change in the average near-infrared reflectance of the stimulated retinal region relative to an equal-sized nonstimulated region (r is the ratio of reflectance between the two retinal areas) was analyzed with a mixed model for repeated measures.
RESULTS. The mixed model showed a significant average decrease in r of 0.14% (95% CI, 0.25 to 0.03) over all subjects induced by bar stimulus cycles, with a gradual return to baseline after stimulus offset, compared with only a 0.04% (95% CI, 0.11+0.20) decrease in r induced by blank, nonstimulated cycles. The mixed model for individuals showed a decreasing linear trend in r over time during bar stimulation, but no decrease for blank cycles in three of five subjects.
CONCLUSIONS. There was a localized decrease in reflectance in response to 780-nm near-infrared light in the retinal region exposed to a visual stimulus, which was significant in three of five subjects. It is presumed that the reflectance change represents the functional activity of the retina in response to a visual stimulus.
Several publications have documented changes in retinal reflectance with light exposure. Riva et al.,6 provided an extensive review of visually evoked hemodynamic retinal changes, including changes in vessel diameter, laser Doppler flowmetry, oximetry, and functional magnetic resonance imaging signals. Crittin and Riva,7 demonstrated the existence of visually evoked light reflectance changes at 570 and 600 nm around the optic nerve head in response to a whole-field flicker stimulus in humans. Tsunoda et al.8 reported a visually evoked reflectance change in the macula of macaque monkeys elicited by a whole-field flash stimulus. De Lint et al.9 showed retinal reflectance changes in humans on a time scale of minutes, with the fundus reflectance during light adaptation and decreasing during dark adaptation. From cortical optical imaging, visually evoked reflectance changes of light in the 400- to 600-nm range are dominated by hemodynamic components, whereas the reflectance changes of near-infrared light over 630 nm may be significantly influenced by changes in the light-scattering properties of cells.4 10 Our group (Tso D et al. IOVS 2004;45:ARVO zE-Abstract 3495) has reported, both in an anesthetized feline preparation and in humans, that stationary localized visual stimuli can evoke focal reflectance changes on the retina by 780-nm near-infrared light (Abràmoff MD et al. IOVS 2004;45:ARVO E-Abstract 5472).11 12
Optical imaging of retinal intrinsic signals in awake humans could offer a novel means by which the temporal and spatial activation properties of the retinal layers can be revealed in normal and pathologic states. The measurements in humans are technically more challenging than those obtained in stabilized animal preparations in which the retina is immobilized and heart rate and breathing can be regulated. The purpose of this study was to characterize the changes in amount of near-infrared 780-nm light reflected from the retina in awake humans in response to a focal visual stimulus.
| Materials and Methods |
|---|
|
|
|---|
|
StimulatorJ software, written by one of the authors (MDA) was used to control the visual stimulus on the miniature LCD display. It is based on the ImageJ program from the National Institutes of Health (Bethesda, MD; developed by Wayne Rasband and available at rsb.info.nih.gov/ij),13 and can display a variety of solid, checkerboard, flickering and static stimuli with configurable timing (described later). For the present study, homogeneous blank-field and vertical bar stimuli were used.
For recording, the tested eye was dilated with a 2.5%phenylephrine0.5% cyclopentolate combination mydriatic at least 30 minutes before the experiment. After corneal anesthesia with topical 0.5% proparacaine hydrochloride ophthalmic solution and insertion of the Burian-Allen lens, the subject dark adapted for at least 5 minutes. The CCD camera was focused on the retinal vessels under visual inspection of the reflectance image and the near-infrared illumination was adjusted to a level to attain an optimally uniform reflectance image. The recording area covered approximately 45° in diameter and was centered on the fovea to include the vascular arcades and optic disc.
Stimuli were presented in a randomized sequence paradigm, as depicted in Figure 2 . At the start of the run, the CCD camera was triggered to begin image collection, while a prestimulus blank period was presented for 3 seconds. After the prestimulus period, a 3-second stimulus period occurred, consisting of either a stationary bright vertical bar or a blank (no stimulus given), as determined by the randomized block design. Throughout the run, a small central fixation mark was flickered at 8 Hz so that the subject could maintain a steady eye position. The prestimulus and stimulus sequence (referred to as a cycle) was repeated 10 times without interruption, typically resulting in 60 seconds (10 cycles) of continuous blank and vertical bar stimulus conditions. The randomized block sequence design ensured that five vertical bars and five blanks were randomly projected in each 60-second period. The CCD image sequence that was recorded as a 60-second run contained 180 frames, each 333 ms in length, and was stored to disc for off-line analysis (described later). The rationale for this protocol was to have an extended recording time without interruption, minimizing the total testing time and reducing potential sources of error associated with the subject exiting and then repositioning in front of the instrument. Eight runs a day could typically be obtained without undue discomfort and fatigue. Runs were limited to 60 seconds, to minimize head and eye movement. In preliminary experiments, longer continuous recording led to unacceptable head movement.
|
Volunteer Inclusion and Exclusion Criteria
All experiments were performed under a protocol that adhered to the tenets of the Declaration of Helsinki and was approved by the Institutional Review Board of the University of Iowa Hospitals and Clinics, including informed consent. In this study, five normal subjects were enrolled. The subjects included two women and three men of ages 51, 50, 41, 41, and 30 years. Subjects 1, 2, 3 and 5 were white; subject 4 was Asian. To satisfy inclusion criteria, subjects underwent a complete ophthalmic history and slit lamp and retinal examination. All subjects had normal eyes with refractive errors less than ±6 D. One eye of each subject was chosen at random for the experiment.
Data Analysis
For off-line analysis, each image was aligned to the first image in each sequence using an affine mutual information-registration algorithm designed by one of the authors (MDA) in ImageJ (available at http://bij.isi.uu.nl).14 This algorithm segments the vessel ridges in each retinal image using the eigenvalues of the Hessian matrix at each pixel and then finds the affine transformation containing the minimum amount of mutual information with the first image in the sequence.12 The amount of x- and y-shift and rotation that had to be applied to each frame for optimal registration was recorded.
As explained earlier, each run in a subject consisted of a 60-second recording period containing 10 cycles (6 sec/cycle), of 18 frames each. Each cycle consisted of a 3-second (nine frames) baseline recording in which no stimulus was presented followed by another 3 seconds (nine frames) in which either a vertical bar stimulus (bar cycle) or no stimulus was presented (blank cycle). The pixel intensities in each recorded frame were analyzed to determine the reflectance ratio (r) of two equal-sized rectangular regions of the retina; one (S) corresponding to the projection of the vertical bar stimulus and the other (N) corresponding to a nonstimulated area oriented perpendicular to the first region (Fig. 3) . The relative reflectance or reflectance ratio r of each image at time t was calculated by dividing the average pixel intensity of S by the average pixel intensity of N, for each frame, for time t from the start of the cycle. Both the S and N regions of interest spared the optic disc and the parafoveal region (Fig. 3) . r(t) was normalized with respect to the frame immediately preceding the stimulus onset by dividing each rj(t) by r(3000 ms). For subsequent analysis, windows for r,
, were defined as the average rof three consecutive frames each, one window before and three after stimulus onset,
and
, for frames 6, 7, 8; 10, 11, 12; 13, 14, 15; and 16, 17, 18 respectively, for both blank and bar stimulus conditions. Corresponding windows were also computed for x- and y-position. A mixed-effect model for repeat measures was used to determine the effect of the presence or absence of a vertical bar visual stimulus (bar versus blank cycles) on the reflectance ratio
.15 The advantage of the mixed model over less robust methods, such as multiple linear regression is that it allows an analysis of all sources of variation separately and independently. This is in contrast to the multiple linear regression, in which all variation is lumped into one error term and which can only account for multiple sources of variation, with a considerable loss of generality. The sources of variation in this study were stimulus, time (window), cycle, cycles in runs, cycles in runs across subjects, and interactions between these sources. This overall mixed model was designed to determine the existence of a change in
over four time windows (from prestimulus window 0 to poststimulus window 3) in all experiments and in all subjects in both bar and blank cycles. The mixed-model analysis was performed using the PROC MIXED feature (SAS ver. 9.1; SAS Institute Inc., Cary, NC). Fixed effects were stimulus state (vertical bar or blank), window, and stimuluswindow interaction. The random effects were runs within subject, cycle within run within subject, and subject. Within each stimulus cycle (bar or blank), a statistical test for linear mean contrast was performed to determine whether mean changes in retinal reflectance, as measured by the mean ratio, decreases linearly across the four sequential time windows (see preceding paragraph). The model described is designed to measure the existence of a change in
in all subjects in all cycles, but it cannot determine the relative size of these changes in individual subjects. Therefore, a second, "individual" mixed model was designed to test for the presence of a change in individual subjects. The presence of a linear trend was defined as the individual mixed model showing a change in reflectance from before to after the stimulus
to
. The effect of changes in eye fixation position on the changes in r was tested by Pearsons product moment calculation. The correlation and its 95% confidence interval (CI) of the reflectance ratio (r) windows with the corresponding small x- and y-position shifts of the fundus that are inevitable during image recording were calculated, to evaluate whether small changes in eye position could account for any of the significant changes in the reflectance ratio (r) that were induced by the bar stimulus.
|
| Results |
|---|
|
|
|---|
To study the significance of this decrease in individual subjects, the individual mixed model was used (Table 1) . In three of five subjects a decrease in reflectance ratio was found after the onset of the stimulus; however, the 95% CI bracketed zero, except in one subject (subject 3). In four of five subjects, no decrease in reflectance ratio was found in blank stimulus cycles, whereas one subject showed a small decrease.
|
Figure 4 shows the time course of the reflectance ratio averaged over all stimulus cycles in a run, for subject 3, the subject with lowest signal-to-noise ratio of the subjects studied. It shows the gradual decrease of r during the first 3 to 4 seconds of a bar cycle, but no change during a blank cycle. After the vertical bar stimulus is turned off, the reflectance ratio gradually returns to baseline before the onset of the next stimulus.
|
Figure 5 shows the spatial distribution of the decrease in reflectance of near-infrared light in a bar and a blank cycle in one of the subjects. In this case, there was a clear spatial correlation between the intrinsic signal of the retina and the retinal location where the stimulus bar was placed, with appropriate control. However, the occurrence of a distinct spatial signal correlating with the shape and location of the stimulus was not frequently observed in most runs, because of a relatively low signal-to-noise ratio, in contrast to our observations in the cat. The location of the stimulus bar in Figure 5 is slightly offset from the spatial distribution of the optical response, which is similar to the response obtained in the anesthetized feline preparation by our group (Tso D et al. IOVS 2004;45:ARVO E-Abstract 3495; Tso D et al. IOVS 2003;44:ARVO E-Abstract 43).
|
| Discussion |
|---|
|
|
|---|
First, slight eye movements could conceivably change the retinal near-infrared reflectance pattern. For example, if the subject is not fixating properly before the onset of the stimulus and then consistently makes a slight, saccadic fixation movement after the stimulus onset of the bar but not the blank stimulus, there could be a change in retinal reflectance time locked to the stimulus condition. We analyzed the eye movement patterns and found no significant correlation of the x- and y-retinal positions with the reflectance changes r, indicating that eye movements did not contribute to the change in reflectance observed. The time course of the reflectance ratio of the intrinsic signal showed a slow decrease after the stimulus was turned on followed by a slow return to baseline after the stimulus was turned off (Fig. 4) . These results are also inconsistent with reflectance changes due to microsaccadic eye movements.
Second, inadequate filtering allowing the visible stimulus light reflected from the retina to "bleed" through the dichroic cold mirror and subsequent barrier filter to reach the CCD sensor could cause artifactual stimulus-related changes in retinal reflectance time, locked to the stimulus. We found no evidence for this artifact in the artificial eye experiments, and the intrinsic optical signal that we measured showed a gradual decrease in reflectance and not the abrupt increase that would be expected from a visible light leak onto the sensor.
Third, nonretinal physiologic factors, such as variation in heart rate, respiration, or blood pressure, could be influenced by an alerting visual stimulus such as the onset of the bar stimulus and could alter the retinal reflectance of near-infrared light. These physiologic factors were not routinely monitored for the data reported herein. However, we performed a separate series of experiments in which the cardiac cycle was synchronized to stimulus presentation/data collection and found the results to be similar to those that were not synchronized (data not shown). It is also difficult to conceive of how such variables would give rise to the localized changes in retinal reflectance that we observed, as opposed to the global changes uniformly spread over the retina.
Because there are important similarities between the retina and brain cortex, the origin of the signal described in this study may be similar to the intrinsic optical signal from visual cortex. In the primate and human neocortex, the intrinsic optical signal is thought to be derived from stimulus-related changes in local hemoglobin volume, the oxygenated state of hemoglobin, and/or the light-scattering properties of activated tissue. Light-scattering changes are thought to be related to stimulus-induced changes in extracellular and intracellular fluids of neurons and glial cells and to changes in the optical properties of membranes due to electrical activity.2 4 In the retina, evidence has also been found recently for the existence of blood flow and light-scattering changes in response to visual stimuli. Investigators in several studies have described in vivo blood flow changes and visible light reflectance changes at the optic nerve head in response to visual stimuli projected onto the retina,7 16 17 Duarte et al.18 have described in vitro increases in 543- and 632-nm laser light-scattering (i.e., decreased reflectance) in the inner plexiform layer in response to electrical activation of the chicken retina, and Tsunoda et al.8 reported the existence of a visually evoked reflectance change in macaque monkeys. The novelty of the present study lies in the documentation of a localized reflectance change in the macula in response to a localized visual stimulus in awake humans. At this time, the exact substrate(s) for the intrinsic optical signals recorded in this study has yet to be defined.
The reflectance decrease in response to the visual stimulus is small and difficult to detect in this first series of experiments, even when using sophisticated image and signal processing techniques, and may not reveal statistically significant stimulus-related changes in some subjects. We are currently improving the acquisition device and image- and signal-analysis techniques to increase the signal-to-noise ratio so that the intrinsic optic signal in response to visual stimulation can be more consistently detected and the origin and location of the signal source within the retinal layers can be defined.
In summary, we recorded from the retina of normal, awake humans a small, localized decrease in reflectance in response to a localized visual stimulus. The signal showed a linearly decreasing (but variable) trend in three out of five subjects. The spatial and temporal properties of these signals suggest that they reflect the functional activity of the retina in a manner parallel to previous optical imaging studies of the sensory neocortex. If the signal-to-noise ratio can be improved, imaging these stimulus-induced changes in near-infrared fundus reflectance in awake humans may, in the future, offer novel methods for the study of the physiology of the human visual system and for the early detection of disease states.
| Acknowledgements |
|---|
| Footnotes |
|---|
Supported by a K-12 Career Development Award, University of Iowa; the Netherlands Organization for Health Related Research (MDA); Research to Prevent Blindness (RPB) (MDA, YK, RK); a National Eye Institute grant for Phase II Small Business Innovative Research (SBIR) (YK); a grant from the National Institute of Biomedical Imaging and BioEngineering (DT); grants from the National Eye Institute (YK, DT, PS, JP, RK); and a grant from Veterans Administration Merit Review and Rehabilitation (RK). JP is an RPB Senior Scientific Investigator.
Submitted for publication January 4, 2005; revised May 26, September 28, and October 28, 2005; accepted December 14, 2005.
Disclosure: M.D. Abràmoff, None; Y.H. Kwon, (P); D. Tso, (P); P. Soliz, Visionquest Biomedical Inc. (I); B. Zimmerman, None; J. Pokorny, None; R. Kardon, (P)
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: Michael D. Abràmoff, Department of Ophthalmology and Visual Sciences, University of Iowa Hospitals and Clinics, 200 Hawkins Drive, Iowa City, IA 52242; michael-abramoff{at}uiowa.edu.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
G. Hanazono, K. Tsunoda, Y. Kazato, K. Tsubota, and M. Tanifuji Evaluating Neural Activity of Retinal Ganglion Cells by Flash-Evoked Intrinsic Signal Imaging in Macaque Retina Invest. Ophthalmol. Vis. Sci., October 1, 2008; 49(10): 4655 - 4663. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Inomata, K. Tsunoda, G. Hanazono, Y. Kazato, K. Shinoda, M. Yuzawa, M. Tanifuji, and Y. Miyake Distribution of Retinal Responses Evoked by Transscleral Electrical Stimulation Detected by Intrinsic Signal Imaging in Macaque Monkeys Invest. Ophthalmol. Vis. Sci., May 1, 2008; 49(5): 2193 - 2200. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Grieve and A. Roorda Intrinsic Signals from Human Cone Photoreceptors Invest. Ophthalmol. Vis. Sci., February 1, 2008; 49(2): 713 - 719. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Okawa, T. Fujikado, T. Miyoshi, H. Sawai, S. Kusaka, T. Mihashi, Y. Hirohara, and Y. Tano Optical Imaging to Evaluate Retinal Activation by Electrical Currents Using Suprachoroidal-Transretinal Stimulation Invest. Ophthalmol. Vis. Sci., October 1, 2007; 48(10): 4777 - 4784. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Hanazono, K. Tsunoda, K. Shinoda, K. Tsubota, Y. Miyake, and M. Tanifuji Intrinsic Signal Imaging in Macaque Retina Reveals Different Types of Flash-Induced Light Reflectance Changes of Different Origins Invest. Ophthalmol. Vis. Sci., June 1, 2007; 48(6): 2903 - 2912. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |