High-resolution mapping of iso-orientation columns by fMRI

Blood-oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) is an important tool for localizing brain functions in vivo. However, the ability of BOLD fMRI to map cortical columnar structures is highly controversial, as the ultimate functional specificity of BOLD remains unknown. Here we report a biphasic BOLD response to visual stimulation in the primary visual cortex of cats. In functional imaging, the initial BOLD signal decrease accurately labeled individual iso-orientation columns. In contrast, the delayed positive BOLD changes indicated the pattern of overall activation in the visual cortex, but were less suited to discriminate active from inactive columns.     

Neurometabolic coupling in cerebral cortex reflects synaptic more than spiking activity

In noninvasive neuroimaging, neural activity is inferred from local fluctuations in deoxyhemoglobin. A fundamental question of functional magnetic resonance imaging (fMRI) is whether the inferred neural activity is driven primarily by synaptic or spiking activity. The answer is critical for the interpretation of the blood oxygen level-dependent (BOLD) signal in fMRI. Here, we have used well-established visual-system circuitry to create a stimulus that elicits synaptic activity without associated spike discharge. In colocalized recordings of neural and metabolic activity in cat primary visual cortex, we observed strong coupling between local field potentials (LFPs) and changes in tissue oxygen concentration in the absence of spikes. These results imply that the BOLD signal is more closely coupled to synaptic activity.     

Quantitative study of changes in oxidative metabolism during visual stimulation using absolute relaxation rates

In the context of quantitative functional MRI (fMRI), deoxyhemoglobin (dHb) content is the essential physiological parameter for calibrating the blood oxygenation level-dependent (BOLD) signal. In studies on humans, the baseline dHb content or its equivalent has been evaluated indirectly by means of carbon dioxide breathing as a physiological reference condition. In this study with normal volunteers, quantitative mapping of baseline dHb content was performed in a direct manner by measuring the reversible contribution of the effective transverse relaxation rate. The BOLD signal change in the visual cortex during 8 Hz flicker visual stimulation was calibrated based on the quantitative map of baseline dHb content. The calibrated relaxation rate change that represents the stimulation-induced fractional change of dHb content decreased by 14% within the activated visual cortex. Simultaneous measurement of cerebral blood flow (CBF) with BOLD showed an increase of 59%. From the calibrated relaxation rate and CBF changes, the cerebral metabolic rate of oxygen (CMRO2) was calculated to increase by 19-28% within the activated visual cortex. The ratio of the CBF increase to the CMRO2 increase was 2-3:1, which agreed well with results of similar quantitative fMRI studies for humans. The method proposed here for quantitative evaluation of the BOLD signal may be applicable not only to fMRI for normal human subjects, but also to physiologically altered or diseased states, because it requires no physiological perturbation.     

Maximizing sensitivity for fast GABA edited spectroscopy in the visual cortex at 7 T

The combination of functional MRI (fMRI) and MRS is a promising approach to relate BOLD imaging to neuronal metabolism, especially at high field strength. However, typical scan times for GABA edited spectroscopy are of the order of 6‐30 min, which is long compared with functional changes observed with fMRI. The aim of this study is to reduce scan time and increase GABA sensitivity for edited spectroscopy in the human visual cortex, by enlarging the volume of activated tissue in the primary visual cortex. A dedicated setup at 7 T for combined fMRI and GABA MRS is developed. This setup consists of a half volume multi‐transmit coil with a large screen for visual cortex activation, two high density receive arrays and an optimized single‐voxel MEGA‐sLASER sequence with macromolecular suppression for signal acquisition. The coil setup performance as well as the GABA measurement speed, SNR, and stability were evaluated. A 2.2‐fold gain of the average SNR for GABA detection was obtained, as compared with a conventional 7 T setup. This was achieved by increasing the viewing angle of the participant with respect to the visual stimulus, thereby activating almost the entire primary visual cortex, allowing larger spectroscopy measurement volumes and resulting in an improved GABA SNR. Fewer than 16 signal averages, lasting 1 min 23 s in total, were needed for the GABA fit method to become stable, as demonstrated in three participants. The stability of the measurement setup was sufficient to detect GABA with an accuracy of 5%, as determined with a GABA phantom. In vivo, larger variations in GABA concentration are found: 14‐25%. Overall, the results bring functional GABA detections at a temporal resolution closer to the physiological time scale of BOLD cortex activation.     

Blood oxygenation level dependent contrast resting state networks are relevant to functional activity in the neocortical sensorimotor system

The relevance of correlations between blood oxygenation level dependent (BOLD) signal changes across the brain acquired at rest (resting state networks, or RSN) to functional networks was tested using two quantitative criteria: (1) the localisation of major RSN correlation clusters and the task-related maxima defined in BOLD fMRI signal changes from the same subjects; and (2) the relative hemispheric lateralisation (LI) of BOLD fMRI signal changes in sensorimotor cortex. RSN were defined on the basis of signal changes correlated with that of a “seed” voxel in the primary sensorimotor cortex. We found a generally close spatial correspondence between clusters of correlated BOLD signal change in RSN and activation maxima associated with hand movement. Conventional BOLD fMRI during active hand movement showed the expected wide variation in relative hemispheric lateralisation of LI for sensorimotor cortex across the subjects. There was a good correlation between LIs for the active hand movement task and the RSN (r=0.74, p<0.001). The RSN thus define anatomically relevant regions of motor cortex and change with functionally relevant variations in hemispheric lateralisation of sensorimotor cortical interactions with hand movement.     

Prolonged fasting impairs neural reactivity to visual stimulation

Previous literature has shown that hypoglycemia influences the intensity of the BOLD signal. A similar but smaller effect may also be elicited by low normal blood glucose levels in healthy individuals. This may not only confound the BOLD signal measured in fMRI, but also more generally interact with cognitive processing, and thus indirectly influence fMRI results. Here we show in a placebo-controlled, crossover, double-blind study on 40 healthy subjects, that overnight fasting and low normal levels of glucose contrasted to an activated, elevated glucose condition have an impact on brain activation during basal visual stimulation. Additionally, functional connectivity of the visual cortex shows a strengthened association with higher-order attention-related brain areas in an elevated blood glucose condition compared to the fasting condition. In a fasting state visual brain areas show stronger coupling to the inferior temporal gyrus. Results demonstrate that prolonged overnight fasting leads to a diminished BOLD signal in higher-order occipital processing areas when compared to an elevated blood glucose condition. Additionally, functional connectivity patterns underscore the modulatory influence of fasting on visual brain networks. Patterns of brain activation and functional connectivity associated with a broad range of attentional processes are affected by maturation and aging and associated with psychiatric disease and intoxication. Thus, we conclude that prolonged fasting may decrease fMRI design sensitivity in any task involving attentional processes when fasting status or blood glucose is not controlled.     

A direct quantitative relationship between the functional properties of human and macaque V5

The nature of the quantitative relationship between single-neuron recordings in monkeys and functional magnetic resonance imaging (fMRI) measurements in humans is crucial to understanding how experiments in these different species are related, yet it remains undetermined. We measured brain activity in humans attending to moving visual stimuli, using blood oxygenation level-dependent (BOLD) fMRI. Responses in V5 showed a strong and highly linear dependence on increasing strength of motion signal (coherence). These population responses in human V5 had a remarkably simple mathematical relationship to previously observed single-cell responses in macaque V5. We provided an explicit quantitative estimate for the interspecies comparison of single-neuron activity and BOLD population responses. Our data show previously unknown dissociations between the functional properties of human V5 and other human motion-sensitive areas, thus predicting similar dissociations for the properties of single neurons in homologous areas of macaque cortex.     

Effects of isoflurane and alpha-chloralose anesthesia on BOLD fMRI responses to ingested l-glutamate in rats

It is important to investigate the effect of anesthesia on blood oxygenation level-dependent (BOLD) signals in an animal model. Many researchers have investigated the BOLD response to visual, sensory, and chemical stimuli in anesthetized rats. There are no reports, however, comparing the differences in the BOLD signal change between anesthetized and conscious rats when a visceral nutrient signal arises. Here, using functional magnetic resonance imaging (fMRI), we investigated the differences in the BOLD signal changes after intragastric administration of l-glutamate (Glu) under three anesthesia conditions: conscious, alpha-chloralose-anesthetized, and isoflurane-anesthetized condition. Under the conscious and alpha-chloralose condition, we observed the significant BOLD signal increase in the medial prefrontal cortex (mPFC), insular cortex (IC), hippocampus, and several hypothalamic regions including the lateral and ventromedial nucleus. In chloralose group, however, gut Glu stimulation induced BOLD signal increase in the prelimbic cortex and orbital cortex, which did not activate in conscious condition. Meanwhile, under isoflurane-anesthetized condition, we did not observe the BOLD signal increase in these areas. BOLD signal intensity in the nucleus of the solitary tract (NTS), to which vagus nerve transmits the visceral information from the gastrointestinal tract, increased in all conditions. Importantly, under conscious condition, we observed increased BOLD signal intensity in several regions related to the metabolic state (i.e. hunger or satiety), such as the mPFC, ventromedial and lateral hypothalamus (LH). Our results suggest that alpha-chloralose and isoflurane anesthesia caused distinct effects on BOLD response to the gut l-Glu stimulation in several brain regions.     

The role of the GABAergic and dopaminergic systems in the brain response to an intragastric load of alcohol in conscious rats

The brain's response to ethanol intake has been extensively investigated using electrophysiological recordings, brain lesion techniques, and c-Fos immunoreactivity. However, few studies have investigated this phenomenon using functional magnetic resonance imaging (fMRI). In the present study, we used fMRI to investigate the blood oxygenation level-dependent (BOLD) signal response to an intragastric (IG) load of ethanol in conscious, ethanol-naive rats. An intragastrically infused 10% ethanol solution induced a significant decrease in the intensity of the BOLD signal in several regions of the brain, including the bilateral amygdala (AMG), nucleus accumbens (NAc), hippocampus, ventral pallidum, insular cortex, and cingulate cortex, and an increase in the BOLD signal in the ventral tegmental area (VTA) and hypothalamic regions. Treatment with bicuculline, which is an antagonist of the gamma-aminobutyric acid A (GABA_A) receptor, increased the BOLD signal intensity in the regions that had shown decreases in the BOLD signal after the IG infusion of 10% ethanol solution, but it did not affect the BOLD signal increase in the hypothalamus. Treatment with SCH39166, which is an antagonist of D1-like receptors, eliminated the increase in the BOLD signal intensity in the hypothalamic areas but did not affect the BOLD signal decrease following the 10% ethanol infusion. These results indicate that an IG load of ethanol caused both a GABA_A receptor–mediated BOLD decrease in the limbic system and the cortex and a D1-like receptor-mediated BOLD increase in the hypothalamic regions in ethanol-naive rats.     

Analysis of fMRI data by blind separation of data in a tiny spatial domain into independent temporal component

Independent Component Analysis (ICA) is a promising tool for the analysis of functional magnetic resonance imaging (fMRI) time series. In these studies, mostly assumed is a spatially independent component map of fMRI data (spatial ICA). In this paper, we assume that the temporal courses of the signal and noises are independent within a Tiny spatial domain (temporal ICA). Then with fast-ICA algorithm, spatially neighboring fMRI data were blindly separated into several temporal courses and were preassumed to be formed by a signal time course and several noise time courses where the signal has the largest correlation coefficient with the reference signal. The final functional imaging was completed for the signals obtained from each voxel. Simulations showed that compared with the spatial ICA method, the new temporal ICA method is more effective than the spatial ICA in detecting weak signal in a fMRI dataset. As background noise, the simulations include simulated Gaussian noise and fMRI data without stimulation. Finally, vivo fMRI tests showed that the excited areas evoked by a visual stimuli are mainly in the region of the primary visual cortex and that evoked by auditory stimuli are mainly in the region of the primary temporal cortex.     

Linking visual gamma to task‐related brain networks—a simultaneous EEG‐fMRI study

There is a growing interest in human gamma‐band oscillatory activity due to its direct link to neuronal populations, its associations with many cognitive processes, and its positive relationship with fMRI BOLD signal. Visual gamma has been successfully detected using concurrent EEG‐fMRI recordings and linked to activity in the visual cortex using voxel‐wise regression analysis. As gamma‐band oscillations reflect predominantly feedforward projections between brain regions, its inclusion in functional connectivity analysis is highly recommended; however, very few studies have investigated this line of research. In the current study, we aimed to explore this gap by asking which fMRI brain network is related to gamma activity induced by the color discrimination task. Advanced denoising strategies and multitaper spectral decomposition were applied to EEG data to detect gamma oscillations, and group independent component analysis was performed on fMRI data to identify task‐related neural networks. Despite using only trials without motor response (50% of the trials), the two neural measures were successfully coupled. One of the six task‐related networks, the occipito‐parietal network, exhibited significant trial‐by‐trial covariations with gamma oscillations. In addition to the expected extrastriate visual cortex, the network encompasses extensive brain activations in the precuneus, bilateral intraparietal, and anterior insular cortices. We argue that the visual cortex is the source of gamma, whereas the remaining brain regions exhibit feedforward and feedback connections related to this oscillatory activity. Our findings provide evidence for the electrophysiological basis of the connectivity revealed by BOLD signal and impart novel insights into the neural mechanism of color discrimination.     

Independent component-cross correlation-sequential epoch (ICS) analysis of high field fMRI time series: direct visualization of dual representation of the primary motor cortex in human

A new technique for functional magnetic resonance imaging (fMRI) time series analysis is presented. The technique referred to here as independent component-cross correlation-sequential epoch (ICS) analysis is a hybrid technique of two standard methodologies of biological signal analysis, namely, data driven methods, represented by independent component analysis, and hypothesis driven methods, represented by a general linear model. The technique successfully identified four functionally discrete areas within the primary sensorimotor cortex (SMI) in normal human subjects based on blood oxygenation level dependent (BOLD) contrast functional magnetic resonance imaging (fMRI) time series performed on a high field (3.0 T) system. Each of the four areas identified corresponded to the four physiological subdivisions of SMI, recognized in primates to be essential for voluntary hand motion, namely, 4 anterior (MI-4a) and 4 posterior (MI-4p) of the primary motor cortex, and 3a and the 'classical' (Brodmann areas 1, 2, and 3b) primary sensory cortex, respectively. ICS analysis appears to be a highly reliable and versatile technique for fMRI time series analysis.     

Quantitative functional imaging of the brain: towards mapping neuronal activity by BOLD fMRI.

Quantitative magnetic resonance imaging (MRI) and spectroscopy (MRS) measurements of energy metabolism (i.e. cerebral metabolic rate of oxygen consumption, CMR(O2)), blood circulation (i.e. cerebral blood flow, CBF, and volume, CBV), and functional MRI (fMRI) signal over a wide range of neuronal activity and pharmacological treatments are used to interpret the neurophysiologic basis of blood oxygenation level dependent (BOLD) image-contrast at 7 T in glutamatergic neurons of rat cerebral cortex. Multi-modal MRI and MRS measurements of CMR(O2), CBF, CBV and BOLD signal (both gradient-echo and spin-echo) are used to interpret the neuroenergetic basis of BOLD image-contrast. Since each parameter that can influence the BOLD image-contrast is measured quantitatively and separately, multi-modal measurements of changes in CMR(O2), CBF, CBV, BOLD fMRI signal allow calibration and validation of the BOLD image-contrast. Good agreement between changes in CMR(O2) calculated from BOLD theory and measured by (13)C MRS, reveals that BOLD fMRI signal-changes at 7 T are closely linked with alterations in neuronal glucose oxidation, both for activation and deactivation paradigms. To determine the neurochemical basis of BOLD, pharmacological treatment with lamotrigine, which is a neuronal voltage-dependent Na(+) channel blocker and neurotransmitter glutamate release inhibitor, is used in a rat forepaw stimulation model. Attenuation of the functional changes in CBF and BOLD with lamotrigine reveals that the fMRI signal is associated with release of glutamate from neurons, which is consistent with a link between neurotransmitter cycling and energy metabolism. Comparisons of CMR(O2) and CBF over a wide dynamic range of neuronal activity provide insight into the regulation of energy metabolism and oxygen delivery in the cerebral cortex. The current results reveal the energetic and physiologic components of the BOLD fMRI signal and indicate the required steps towards mapping neuronal activity quantitatively by fMRI at steady-state. Consequences of these results from rat brain for similar calibrated BOLD fMRI studies in the human brain are discussed.     

Temporal feature of BOLD responses varies with temporal patterns of movement

Which brain sites represent the final form of motor commands that encode temporal patterns of muscle activities? Here, we show the possible brain sites which have activity equivalent to the motor commands with functional magnetic resonance imaging (fMRI). We hypothesized that short-temporal patterns of movements or stimuli are reflected in blood-oxygenation-level-dependent (BOLD) responses and we searched for regions representing the response. Participants performed two temporal patterns of tapping and/or listened to the same patterns of auditory stimuli in a 3T fMRI. The patterns were designed to have the same number (11) of events and the same duration, but different temporal distribution of events. The 11 events were divided into two parts (10 repetitive taps and one stand-alone tap) and the interval of the two parts was 3s. The two patterns had reverse order of the two parts. The results revealed that different temporal patterns of auditory stimuli were represented in different temporal features of BOLD responses in the bilateral auditory cortex, whereas different temporal patterns of tapping were reflected in contralateral primary motor cortex and the ipsilateral anterior cerebellum. In bilateral premotor cortex, supplementary motor area, visual cortex, and posterior cerebellum, task-related BOLD responses were exhibited, but their responses did not reflect the temporal patterns of the movement and/or stimuli. One possible explanation is that the neuronal activities were similar for the two patterns in these regions. The sensitivity of the BOLD response to the temporal patterns reflects local differences in functional contributions to the tasks. The present experimental design and analysis may be useful to reveal particular brain regions that participate in multiple functions.     

Evoked local field potentials can explain temporal variation in blood oxygenation level-dependent responses in rat somatosensory cortex

The aim of this study was to explain the temporal variations between subjects in the blood oxygenation level-dependent (BOLD) response. Somatosensory responses were elicited with the electrical forepaw stimulus at a frequency of 10 Hz in urethane-anesthetized rats, and functional magnetic resonance imaging (fMRI) with BOLD contrast and local field potential (LFP) measurements were performed simultaneously. BOLD fMRI activation was evaluated by two different models, one based on the stimulus paradigm (the block model) and the other on the simultaneously measured evoked LFP responses. In the initial analysis, the LFP model captured the BOLD activation in the primary somatosensory cortex in all cases, and the block model in 10 of 12 rats. A statistical comparison of the two models revealed that the LFP-derived model was able to explain additional BOLD variation over the block model in the somatosensory cortex in nine of 12 rats. These results suggest that there is more information regarding neuronal activity in the BOLD signal than can be exploited using the block model alone. Furthermore, the hemodynamic coupling remains unchanged in the case of temporally variable BOLD signals.     

Alcohol induced region-dependent alterations of hemodynamic response: implications for the statistical interpretation of pharmacological fMRI studies

Worldwide, ethanol abuse causes thousands of fatal accidents annually as well as innumerable social dysfunctions and severe medical disorders. Yet, few studies have used the blood oxygenation level dependent functional magnetic resonance imaging method (BOLD fMRI) to map how alcohol alters brain functions, as fMRI relies on neurovascular coupling, which may change due to the vasoactive properties of alcohol. We monitored the hemodynamic response function (HRF) with a high temporal resolution. In both motor cortices and the visual cortex, alcohol prolonged the time course of the HRF, indicating an overall slow-down of neurovascular coupling rather than an isolated reduction in neuronal activity. However, in the supplementary motor area, alcohol-induced changes to the HRF suggest a reduced neuronal activation. This may explain why initiating and coordinating complex movements, including speech production, are often impaired earlier than executing basic motor patterns. Furthermore, the present study revealed a potential pitfall associated with the statistical interpretation of pharmacological fMRI studies based on the general linear model: if the functional form of the HRF is changed between the conditions data may be erroneously interpreted as increased or decreased neuronal activation. Thus, our study not only presents an additional key to how alcohol affects the network of brain functions but also implies that potential changes to neurovascular coupling have to be taken into account when interpreting BOLD fMRI. Therefore, measuring individual drug-induced HRF changes is recommended for pharmacological fMRI.     

Gain control in the response of human visual cortex to plaids

A recent intrinsic signal optical imaging study in tree shrew showed, surprisingly, that the population response of V1 to plaid patterns comprising grating components of equal contrast is predicted by the average of the responses to the individual components (MacEvoy SP, Tucker TR, Fitzpatrick D. Nat Neurosci 12: 637-645, 2009). This prompted us to compare responses to plaids and gratings in human visual cortex as a function of contrast and orientation. We found that the functional MRI (fMRI) blood oxygenation level-dependent (BOLD) responses of areas V1-V3 to a plaid comprising superposed grating components of equal contrast are significantly higher than the responses to a single component. Furthermore, the orientation response profile of a plaid is poorly predicted from a linear combination of the responses to its components. Together, these results indicate that the model of MacEvoy et al. (2009) cannot, without modification, account for the fMRI BOLD response to plaids in human visual cortex.     

Cortical and subcortical correlates of electroencephalographic alpha rhythm modulation

Neural correlates of electroencephalographic (EEG) alpha rhythm are poorly understood. Here, we related EEG alpha rhythm in awake humans to blood-oxygen-level-dependent (BOLD) signal change determined by functional magnetic resonance imaging (fMRI). Topographical EEG was recorded simultaneously with fMRI during an open versus closed eyes and an auditory stimulation versus silence condition. EEG was separated into spatial components of maximal temporal independence using independent component analysis. Alpha component amplitudes and stimulus conditions served as general linear model regressors of the fMRI signal time course. In both paradigms, EEG alpha component amplitudes were associated with BOLD signal decreases in occipital areas, but not in thalamus, when a standard BOLD response curve (maximum effect at approximately 6 s) was assumed. The part of the alpha regressor independent of the protocol condition, however, revealed significant positive thalamic and mesencephalic correlations with a mean time delay of approximately 2.5 s between EEG and BOLD signals. The inverse relationship between EEG alpha amplitude and BOLD signals in primary and secondary visual areas suggests that widespread thalamocortical synchronization is associated with decreased brain metabolism. While the temporal relationship of this association is consistent with metabolic changes occurring simultaneously with changes in the alpha rhythm, sites in the medial thalamus and in the anterior midbrain were found to correlate with short time lag. Assuming a canonical hemodynamic response function, this finding is indicative of activity preceding the actual EEG change by some seconds.     

The human brain representation of odor identification

Odor identification (OI) tests are increasingly used clinically as biomarkers for Alzheimer's disease and schizophrenia. The aim of this study was to directly compare the neuronal correlates to identified odors vs. nonidentified odors. Seventeen females with normal olfactory function underwent a functional magnetic resonance imaging (fMRI) experiment with postscanning assessment of spontaneous uncued OI. An event-related analysis was performed to compare within-subject activity to spontaneously identified vs. nonidentified odors at the whole brain level, and in anatomic and functional regions of interest (ROIs) in the medial temporal lobe (MTL). Parameter estimate values and blood oxygenated level-dependent (BOLD) signal curves for correctly identified and nonidentified odors were derived from functional ROIs in hippocampus, entorhinal, piriform, and orbitofrontal cortices. Number of activated voxels and max parameter estimate values were obtained from anatomic ROIs in the hippocampus and the entorhinal cortex. At the whole brain level the correct OI gave rise to increased activity in the left entorhinal cortex and secondary olfactory structures, including the orbitofrontal cortex. Increased activation was also observed in fusiform, primary visual, and auditory cortices, inferior frontal plus inferior temporal gyri. The anatomic MTL ROI analysis showed increased activation in the left entorhinal cortex, right hippocampus, and posterior parahippocampal gyri in correct OI. In the entorhinal cortex and hippocampus the BOLD signal increased specifically in response to identified odors and decreased for nonidentified odors. In orbitofrontal and piriform cortices both identified and nonidentified odors gave rise to an increased BOLD signal, but the response to identified odors was significantly greater than that for nonidentified odors. These results support a specific role for entorhinal cortex and hippocampus in OI, whereas piriform and orbitofrontal cortices are active in both smelling and OI. Moreover, episodic as well as semantic memory systems appeared to support OI.     

Light stimulus frequency dependence of activity in the rat visual system as studied with high-resolution BOLD fMRI

The neurophysiology of the rodent visual system has mainly been investigated by invasive and ex-vivo techniques providing fragmented data. This area of research has been deprived of functional MRI studies based on blood oxygenation level dependent (BOLD) contrast, which allows a whole brain approach with a high spatial and temporal resolution. In the present study, we looked at the neurovascular response properties of the visual system of the pigmented rat, focusing on the visual cortex (VC), the superior colliculus (SC) and the flocculus-paraflocculus of the cerebellum (FL-PFL), using BOLD fMRI under domitor anesthesia. Visual stimulation was performed monocularly or binocularly while flashing light from a strobe unit was presented. For each structure, we assessed the flashing frequency that evoked the optimal BOLD response: Neither the VC nor the FL-PFL displayed frequency dependence during monocular visual stimulation, but were most sensitive to low frequencies (1-5 Hz) when flashing light was provided binocularly. The SC responded optimally to high flashing rates (8-12 Hz) during both monocular and binocular stimulation. The signal intensity changes in the VC and FL-PFL were locked to the stimulation period, whereas the BOLD response in the SC showed a similar onset but a very slow recovery at offset. The VC and FL-PFL, but not the SC, showed signs of binocular competition. The observed correlation between frequency-dependent responses of different visual areas during binocular visual presentation suggests a functional relationship between the VC and FL-PFL rather than between the SC and FL-PFL.