Origin of Synchronized Low-Frequency Blood Oxygen Level–Dependent Fluctuations in the Primary Visual Cortex

Synchronized BOLD fluctuations. A, Time series data from 1 subject. Six minutes of data are shown. Traces represent average BOLD signal intensity from regions of interest in the bilateral primary visual cortex (V1), left LGN, and bilateral frontal eye fields (FEF). B, Cross-correlation analysis represents full 14-minute resting state data from regions listed previously. Bilateral V1 shows strong focal correlation, with similar but diminished correlation between V1 and LGN. No correlation is seen between V1 and FEF. Bilateral FEF show much broader temporal correlation. C, Summary resting state BOLD cross-correlation from 6 subjects. Green traces show cross-correlograms for individual subjects, and black traces show an average cross-correlogram across subjects (bilateral V1, n = 6; left V1-left LGN, n = 3; left V1-left FEF, n = 2, left FEF-right FEF: n = 2).

Vascular low-pass filtering of neurovascular coupling. A, Data in the first column represent BOLD responses to a 6-Hz flashing radial checkerboard pattern turned ON and OFF at frequencies ranging from 0.01 to 0.5 Hz. Data in the second column represent BOLD responses to single flashes of 1 screen refresh duration (13 ms) of the same radial checkerboard pattern presented at frequencies ranging from 0.1 to 4 Hz. B, Peak-to-peak response for 1 period of the stimulus. Individual points plotted are responses from each subject of the amplitude of the BOLD response averaged across all cycles of stimulus presentation. Left and right hemifield responses are plotted as separate data points on the left and are averaged for each data point on the right. Error bars represent 1 standard error of the mean. C, Mean BOLD responses compared with baseline (blank screen stimulus), averaged across the entire trial for each frequency condition. Data points and error bars are as in Fig 3B.

Frequencies present in spontaneous BOLD fluctuations in the primary visual cortex. A, Full-field (inset) and zoomed-in power spectra of resting-state BOLD fluctuations. Spectra are computed for each subject, normalized to the unit area under the curve, and averaged across subjects. Shaded regions represent ±1 standard error of the mean as a function of frequency. B, Superimposed vascular filtering data of Fig 3B and the resting-state power spectrum from Fig 4A. Frequencies between 0.05 and 0.2 Hz are under-represented if frequency distribution of spontaneous BOLD fluctuations arises from vascular filtering alone.

Power in spontaneous BOLD fluctuations in resting and visual stimulation conditions. Resting state (blue) condition shows the shift of frequency distribution to higher frequencies compared with natural visual stimulation condition (red). Spectra are averaged across 6 subjects for each condition. Shaded areas represent standard error of the mean as a function of frequency across the population. The pair-wise t test for each frequency range (0–0.02 Hz and 0.02–0.05 Hz) indicates statistically significant differences between the 2 conditions for each of the 2 ranges.

Amplitude modulation of MEG signal intensity. A, Procedure for calculating amplitude modulation for each frequency band of the MEG signal intensity. Raw data were first bandpass-filtered into δ, α, β, and γ bands. Data were then full-wave-rectified, and an envelope function was calculated to measure the variation in amplitude with time. Approximately 5 seconds of data are shown. B, Envelope function from Fig 6A for 1 subject. Approximately 7 minutes of data are shown. Note the similarity in appearance of low-frequency fluctuations to BOLD fluctuations shown in Fig 2.