|
|
||||||||
urbil3
1 Advanced Imaging Laboratories, The John P. Robarts Research Institute, London, Ontario N6A 5K8, Canada; 2 Lucent Technologies Bell Laboratories, Murray Hill, New Jersey 07974; and 3 Centre for Magnetic Resonance Research, University of Minnesota, Minneapolis, Minnesota 55455
| |
ABSTRACT |
|---|
|
|
|---|
Menon, Ravi S., Seiji Ogawa, John P. Strupp, and Kâmil U
urbil. Ocular dominance in human V1 demonstrated by functional magnetic resonance imaging. J. Neurophysiol. 77: 2780-2787, 1997. Very high resolution functional magnetic resonance imaging (fMRI) at a 4 Tesla (T) magnetic field was used to map ocular dominance regions in the human visual cortical layers using the blood oxygen level dependent (BOLD) contrast mechanism. The fMRI response from primary visual cortex (V1) exhibited a distribution of ocular dominance reminiscent of the single-cell recordings of Hubel and Wiesel. Pixels could be grouped into seven categories varying from left-only response to binocular-only response to right-only responses. Nonspecific responses were found in the MRI-visible draining veins as well as in the parenchyma. Although large vessel BOLD signals are easily detectable, regardless of field strength, they demonstrate a fMRI response to photic input that could not be used to distinguish ocular dominance. The difference in BOLD response between a region activated by one eye and that activated by the other is only 2.9% on average. This necessitates the use of a difference paradigm to visualize the regions of ocular dominance accurately. The data show that BOLD-based fMRI is sensitive to neuronal activity in cortical columns when using differential techniques, opening up the possibility of mapping specialized populations of neurons in humans that are not accessible to electrophysiological or other methods of invasive mapping.
In primates, axons from the left and right eyes terminate in monocular laminae of the lateral geniculate body. From this nucleus, geniculostriate projections to primary visual cortex (V1) continue to reflect either left or right eye input and terminate in layer IVC of V1 where they are arranged in a system of roughly parallel alternating stripes known as ocular dominance columns (ODCs). In cortical layers above and below IVC, cortical neurons vary in the strength of their response to inputs from the two eyes. This response can vary from exclusive domination by one eye to equal influence of both as one makes a tangential penetration with an electrode through the striate cortex (Hubel and Wiesel 1962 Subjects
Five normal subjects [4 right eye and right hand dominant (3 male, 1 female), 1 left eye and left hand dominant (1 male), ages 25-32] gave informed consent before participating in this study. All had previous fMRI experience. Approval for this protocol was given by the University of Western Ontario Review Board for Health Sciences Research involving human subjects, and radio frequency power deposition guidelines established for clinical scanners by the Food and Drug Administration were adhered to.
Activation tasks
A single round red (635 nm, 300 mcd, Radio Shack "Jumbo") light-emitting diode (LED), placed above the subject's head and driven by a stimulator (GRASS Instruments, Quincy, MA) at 8 Hz was used for photic stimulation. The LED subtended ~10° of visual field. The subjects were directed to fixate on the position of the LED, which was gated on or off by the scanner, but to generate monocular input, subjects were instructed via a special headset (Resonance Technology, Van Nuys, CA) when to close either the left or the right eyelid gently as the scanning progressed. Binocular stimulation ("B" state), sandwiched between periods of darkness ("D" state), was used to identify primary visual cortex. Monocular stimulation using the open left eye ("L" state) or open right eye ("R" state) was used to delineate ocular dominance regions. Each of the states, B, D, L, and R were 15 s in duration.
MRI studies
All MRI experiments were performed on a 4 T whole-body human imager [varian (Palo Alto, CA)/Siemens (Erlangen, Germany)] with a 7.6-cm diam double-balanced, distributed capacitance transmit/receive RF coil built into the bottom of a rigid Plexiglas head-holder. Subject head motion was restrained with an integrated foam padded vice that pressed against the sides of the head, around the headphones. In some cases, the surface coil was removed and a head coil inserted for anatomic imaging, without disturbing the subject.
FMRI analysis
All image analysis was done using Stimulate v5.0 (Strupp 1996
An oblique T1-weighted anatomic image parallel to the calcarine fissure, acquired as described above, is shown in Fig. 2A. Overlayed on this image, in Fig. 2B, is the functional map determined for binocular stimulation. This high-resolution functional map of binocular activation was derived from the first four periods of the visual paradigm shown in Fig. 1A using the cross-correlation function in Fig. 1B. The activation was observed to be almost exclusively confined to the occipital pole and to coregister extremely well with the cortical gray matter in the primary visual areas. Areas that appeared yellow in the color scale (>10% change on binocular stimulation) were in regions of cortex with visible venous vasculature.
To assess the significance of the changes measured, we first examine their magnitude relative to the noise. The mean fractional image intensity change of the pixels we characterize as responding during binocular stimulation is 2.7 ± 0.5% (mean ± SE; n = 5). The mean fractional change in the same regions of cortex during monocular stimulation was similar to the binocular case at 2.9 ± 0.5% (n = 5). We have excluded visible veins from this measurement. In the local draining veins, changes >20% can be observed. We observe the maximum gray matter changes between binocular and dark stimulus conditions to be ~5% with our slice thickness of 4 mm, and, taking this as the maximum expected parenchymal change, we require the fMRI mapping procedure to be sensitive to changes of <5% in a single pixel to detect ocular dominance regions. Given the stability of the head and our instrument, the signal-to-noise ratio in the images and our desired confidence levels, differences of <1.5% between states, are not significant in a single pixel, using FLASH, at the P = 0.05 level. Therefore we sort our L-R differences into bins of 0-1.5% (binocular or no significant difference), 1.6-2.5%, 2.6-3.5%, and 3.5-4.5% (monocular). The difference could be positive or negative depending on which eye is dominant. A positive difference indicated right eye dominance.
![]()
INTRODUCTION
Abstract
Introduction
Methods
Results
Discussion
References
). Only in layer IVC do the cells receive innervation from exclusively one eye, and hence the classical periodic columnar pattern is only observed in layer IVC. In nonhuman primates, the organization of ODCs has been studied by histological stains and autoradiography (Hubel and Wiesel 1977
; LeVay et al. 1985
; Kennedy et al. 1976
), by microelectrode recordings (Hubel and Wiesel 1976), real-time optical imaging using voltage-sensitive dyes (Salzberg et al. 1973
), and optical imaging of intrinsic signals (Grinvald et al. 1986
, 1991
; Tso et al. 1990
). In humans, the ODCs have been demonstrated as interdigitated stripes of ~1 mm in width by postmortem histochemical staining for cytochrome oxidase in striate cortex by Horton et al. (1984, 1990), but a noninvasive technique for examining human striate cortex organization on the scale of cortical functional subunits has not been available.
, 1991
; Malonek and Grinvald 1996
; Tso et al. 1990
). The blood oxygen level dependent (BOLD) technique (Ogawa et al. 1990
a,b; Turner et al. 1991
) on which the vast majority of cortical mapping using functional magnetic resonance imaging (fMRI) is based (Bandettini et al. 1992
; Kwong et al. 1992
; Ogawa et al. 1992
), is also sensitive to the hemodynamic changes in the local vasculature, which suggests that, in principle, cortical column organization could be mapped noninvasively with the use of fMRI as well. Although the optical data demonstrate that the capillary bed hemodynamic response is sufficiently confined to layer IVC of the cortical column (the ODC), localization of the mapping signal in fMRI with respect to the active cortical area and vascular tree is still quite controversial (Cohen and Bookheimer 1994
; Duyn et al. 1994
; Frahm et al. 1994
; Kim et al. 1994
; Kwong 1995
; Lai et al. 1993
; Menon et al. 1993
). At issue here is the belief that BOLD signals coming from large vessels may dominate those coming from the microvasculature, particularly on conventional MRI scanners. This is problematic, because it raises concern that macrovascular changes distal to the actual site of neuronal activity can occur. This would place a fundamental limit on correlation of fMRI activation maps and neuronal activity. Because cells above and below the ODC defined in layer IVC can also respond in varying degrees as mentioned above, the vascular response may not be confined to layer IVC, but might traverse the entire depth of the cortical ribbon. Thus it may be expected, particularly with the resolution used in fMRI, that each pixel would contain variable degrees of left or right eye dominance, and that confinement of the fMRI mapping signal to layer IVC may not be feasible.
; Malonek and Grinvald 1996
) with considerably lower resolution fMRI experiments performed at a magnetic field strength of 4 T (Menon et al. 1995
), the temporal response and sign of the optical signal from the capillary beds in V1 of awake monkeys parallels that of the fMRI time course in certain regions of striate cortex in humans. The remarkable similarity of the results obtained from these two techniques, along with other theoretical and experimental evidence for increasing capillary bed BOLD contributions at higher magnetic field strengths (Bandettini et al. 1994
; Gati et al. 1997
; Menon et al. 1993
, 1994
; Ogawa et al. 1993
; Song et al. 1994
) suggest that a component of the fMRI mapping signal arises from microvasculature, at least at the highest magnetic fields available for human research (4 T).
). Using a high-resolution fMRI pulse-sequence with imaging hardware and parameters optimized at three different field strengths, we have found that the SNR at 4 T is at least three times higher than at the much more commonly available 1.5 T field strength (Gati et al. 1997
). This increase is sufficiently large to attempt imaging of ocular dominance in human V1, where ODCs are ~0.8-1 mm on a side for a column and 5-10 mm long in humans (Horton and Hedley-White 1984
; Horton et al. 1990
). Using a simple visual paradigm in combination with an optimized radio-frequency (RF) coil, head restraint, and the enhanced SNR provided by 4 T, we have been able to demonstrate adjacent image pixels in human V1 that respond predominantly to left or right eye photic input, as well as many that respond only to binocular input.
![]()
METHODS
Abstract
Introduction
Methods
Results
Discussion
References
; Menon et al. 1993
; Ogawa et al. 1992
). Briefly, 256 centrically ordered phase-encoding steps segmented in 4 blocks of 64 interleaved steps with a magnetization preparation inversion time (TI) of 1.2 s were used for this high-resolution FLASH sequence (Haase et al. 1986
). For these anatomic images, the echo time (TE) was 5 ms, the repetition time (TR) was 10 ms, the slice thickness was 4 mm, and the segment interleave time was 4 s. From these multislice sagittal images the calcarine fissure was easily identified (Menon et al. 1993
). Five equally spaced and abutting oblique slices of 4-mm thickness were chosen parallel to the calcarine fissure at the posterior occipital pole. This location and orientation was chosen on the basis of previous MRI myeloarchitectonic analysis of layer IVC in striate cortex (Clark et al. 1992
). The chosen orientation allows the columns to run perpendicular to the slice plane. Because the columns are expected to be 5-10 mm in length in humans and 0.8-1.2 mm on a side (Horton et al. 1990
), it is expected that they will be perpendicular to the 4-mm-thick slice in at least a few local regions, but not necessarily in the whole slice, given that the calcarine sulcus is rarely straight in human subjects. Anatomic imaging in the same manner as described above was also performed on these five slices.
; Tootell et al. 1995
), this simple paradigm allows rapid identification of activated primary visual cortex, which is essential to set up a single slice study.

View larger version (18K):
[in a new window]
FIG. 1.
Visual paradigm and binocular reference vector. A: the visual paradigm consisted of 24 serial states of darkness (D), binocular stimulation (B), and monocular stimulation of either the left (L) or right (R) eye as shown at top. Each state lasted 15 s during which time 1 image was acquired. B: to identify pixels that responded to photic stimulation in either or both eyes, a cross-correlation of images 1 to 9 was performed with the reference vector shown, shifted by 1 image to account for the hemodynamic lag time.
) provided by the University of Minnesota, running on a Sun UltraSparc 140. To initially identify primary visual cortex in the pilot multislice scanning protocol, a Student's t-test was done, on a pixel by pixel basis, to compare the dark and binocularly driven states for each slice. Only pixels with a statistically significant difference between B and D states (P < 0.05) were included. The functional maps were overlaid on the corresponding anatomic slices. The ocular dominance mapping experiment was then performed on the single slice determined most suitable from this analysis.
) shown in Fig. 1B with slightly Gaussian blurred image data. An extremely high correlation value of 0.75 was used. Because a centric phase-encoding scheme was used, the image intensity is dominated by the first few lines of the acquisition, which occur well within the hemodynamic lag time of 5-8 s. Therefore the reference waveform was shifted by one image to account for the hemodynamic lag in response to visual stimulation as we have demonstrated previously (Menon et al. 1995
). The net result is a clean map ("mask") of the cortical ribbon and vasculature as shown in Fig. 2B in which the pixels that we detect are activated by either or both eyes. The blurring ensures that all layers of the striate cortex are included in the mask. This mask was then used as a template on the original (unfiltered) image data to generate time courses from each 0.55 mm by 0.55 mm by 4 mm voxel that was deemed activated by this procedure. No other image or background thresholding, pixel clustering, or region-of-interest (ROI) limitation was used.

View larger version (150K):
[in a new window]
FIG. 2.
Functional magnetic resonance imaging (fMRI) maps of striate cortex to binocular and monocular visual input. A: T1-weighted MRI image of subject 1. B: map of pixels that satisfy the cross-correlation criteria shown in Fig. 1B at a correlation value of 0.75. Activation is coded by color [red (1%)
yellow (>20%)] overlaid on a high-resolution T1 weighted MRI image of the same subject. No significant negative changes were observed. C: map demonstrating pixels responding predominantly to left eye monocular stimulation in reds [red (1%)
yellow (>10%)], whereas those responding primarily to right eye monocular stimulation are shown in blues [blue (1%)
violet (>10%)]. Pixels were chosen on the basis of the ocular dominance histograms shown in Fig. 4. D: expansion of part of C. In A-C the field of view is 14 by 14 cm.
). Binning was accomplished by doing a Student's t-test for significant differences between right and left eye stimulation periods at the P = 0.05 level. Those pixels that had significant differences were then binned according to the percentage difference between the two monocular conditions. Differences in right minus left response or left minus right response that were <1.5% were not significant at this P value and were considered binocularly responding, those between 1.5 and 4.5% were deemed partially dominated by the appropriate eye, and those exceeding 4.5% were considered monocular. Maps of ocular dominance were made based on these histograms for each subject.
![]()
RESULTS
Abstract
Introduction
Methods
Results
Discussion
References

View larger version (40K):
[in a new window]
FIG. 3.
Time course of image intensity corresponding to pixels activated during the binocular states. The ROI from which they are extracted is shown in Fig. 2B. Contributions from visible vessels have been excluded. The fractional signal change in this region was 2.48%.

View larger version (48K):
[in a new window]
FIG. 4.
Histogram of ocular dominance for all 5 subjects and the average across subjects. Positive differences indicate right eye dominance, and negative changes indicate left eye dominance. The bins are 0-1.5% (binocular or no significant difference), 1.6-2.5%, 2.6-3.5%, and 3.5-4.5% (monocular). At a P value of 0.05, differences under 1.5% could not be detected in single pixels. The few pixels whose changes were >4.5% were placed in the 4.5% categories.

View larger version (41K):
[in a new window]
FIG. 5.
Time course of image intensity corresponding to pixels activated during the monocular states. The average time course of the 4 red pixels (left eye) and the 4 blue pixels (right eye) marked in Fig. 2D are shown. These pixels are chosen on the basis of the ocular dominance histogram for subject 1 in Fig. 4 and have changes >2.5%.

View larger version (17K):
[in a new window]
FIG. 6.
Effects of filtering on ocular dominance response. For pixels such as those shown in Fig. 2D, whose time courses are seen in Fig. 5, the difference between left and right responses was observed as a function of blurring. Half the contrast between eyes is lost at a Gaussian filter full-width at half-maximum (FWHM) of 1.6 pixels (875 µm).
![]()
DISCUSSION
Abstract
Introduction
Methods
Results
Discussion
References
; Horton et al. 1990
). Our blurring procedure shown in Fig. 6 demonstrates that 50% of left-right differentiation is wiped out by a Gaussian filter of full-width at half-maximum of 875 µm (1.6 pixels). Thus, on average, the ocular dominance regions are 875 µm on a side. Third, the ocular dominance areas appear confined to the primary visual cortex as expected. Fourth, we can demonstrate that the temporal behavior of adjacent pixels of assigned ocular dominance are consistent with that of the stimulus paradigm (Fig. 5). These four observations lend strong support to our identification of the activated regions as regions of ocular dominance and not random noise.
| |
ACKNOWLEDGEMENTS |
|---|
We thank J. Gati and P. Andersen for technical assistance.
This work was supported by Grant MT13350 and a salary support award from the Medical Research Council of Canada to R. S. Menon and National Institutes of Health Grants 1 RO1 EY-11551-01 to R. S. Menon and RR-08079 to K. U
urbil.
| |
FOOTNOTES |
|---|
Address for reprint requests: R. S. Menon, Advanced Imaging Laboratories, The John P. Robarts Research Institute, PO Box 5015, 100 Perth Drive, London, Ontario N6A 5K8, Canada.
Received 6 February 1996; accepted in final form 24 January 1997.
| |
REFERENCES |
|---|
|
|
|---|
R2* and functional contrast to noise ratio.
Proc. Soc. Mag. Reson.
2: 434, 1994.This article has been cited by other articles:
![]() |
E. Yacoub, N. Harel, and K. Ugurbil High-field fMRI unveils orientation columns in humans PNAS, July 29, 2008; 105(30): 10607 - 10612. [Abstract] [Full Text] [PDF] |
||||
![]() |
S.-G. Kim and M. Fukuda Lessons from fMRI about Mapping Cortical Columns Neuroscientist, June 1, 2008; 14(3): 287 - 299. [Abstract] [PDF] |
||||
![]() |
D. L. Adams, L. C. Sincich, and J. C. Horton Complete Pattern of Ocular Dominance Columns in Human Primary Visual Cortex J. Neurosci., September 26, 2007; 27(39): 10391 - 10403. [Abstract] [Full Text] [PDF] |
||||
![]() |
C.-H. Moon, M. Fukuda, S.-H. Park, and S.-G. Kim Neural Interpretation of Blood Oxygenation Level-Dependent fMRI Maps at Submillimeter Columnar Resolution J. Neurosci., June 27, 2007; 27(26): 6892 - 6902. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. O. Duncan, P. A. Sample, R. N. Weinreb, C. Bowd, and L. M. Zangwill Retinotopic Organization of Primary Visual Cortex in Glaucoma: A Method for Comparing Cortical Function with Damage to the Optic Disk Invest. Ophthalmol. Vis. Sci., February 1, 2007; 48(2): 733 - 744. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Fukuda, C.-H. Moon, P. Wang, and S.-G. Kim Mapping Iso-Orientation Columns by Contrast Agent-Enhanced Functional Magnetic Resonance Imaging: Reproducibility, Specificity, and Evaluation by Optical Imaging of Intrinsic Signal. J. Neurosci., November 15, 2006; 26(46): 11821 - 11832. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Vanzetta, R. Hildesheim, and A. Grinvald Compartment-Resolved Imaging of Activity-Dependent Dynamics of Cortical Blood Volume and Oximetry J. Neurosci., March 2, 2005; 25(9): 2233 - 2244. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. A. Sheth, M. Nemoto, M. Guiou, M. Walker, N. Pouratian, N. Hageman, and A. W. Toga Columnar Specificity of Microvascular Oxygenation and Volume Responses: Implications for Functional Brain Mapping J. Neurosci., January 21, 2004; 24(3): 634 - 641. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Turner and T. Jones Techniques for imaging neuroscience Br. Med. Bull., March 1, 2003; 65(1): 3 - 20. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. C. Silva and A. P. Koretsky Laminar specificity of functional MRI onset times during somatosensory stimulation in rat PNAS, November 12, 2002; 99(23): 15182 - 15187. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Shoham and A. Grinvald The Cortical Representation of the Hand in Macaque and Human Area S-I: High Resolution Optical Imaging J. Neurosci., September 1, 2001; 21(17): 6820 - 6835. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Q. Duong, D.-S. Kim, K. Ugurbil, and S.-G. Kim Localized cerebral blood flow response at submillimeter columnar resolution PNAS, August 23, 2001; (2001) 191101098. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. M. Ramsden, C. P. Hung, and A. W. Roe Real and Illusory Contour Processing in Area V1 of the Primate: a Cortical Balancing Act Cereb Cortex, July 1, 2001; 11(7): 648 - 665. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Hess, D. Stiller, T. Kaulisch, P. Heil, and H. Scheich New Insights into the Hemodynamic Blood Oxygenation Level-Dependent Response through Combination of Functional Magnetic Resonance Imaging and Optical Recording in Gerbil Barrel Cortex J. Neurosci., May 1, 2000; 20(9): 3328 - 3338. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Vanzetta and A. Grinvald Increased Cortical Oxidative Metabolism Due to Sensory Stimulation: Implications for Functional Brain Imaging Science, November 19, 1999; 286(5444): 1555 - 1558. [Abstract] [Full Text] |
||||
![]() |
W. Chen, X.-H. Zhu, K. R. Thulborn, and K. Ugurbil Retinotopic mapping of lateral geniculate nucleus in humans using functional magnetic resonance imaging PNAS, March 2, 1999; 96(5): 2430 - 2434. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Q. Duong, D.-S. Kim, K. Ugurbil, and S.-G. Kim Localized cerebral blood flow response at submillimeter columnar resolution PNAS, September 11, 2001; 98(19): 10904 - 10909. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME |