Skip to main content
Advertisement

Main menu

  • Home
  • Content
    • Current Issue
    • Accepted Manuscripts
    • Article Preview
    • Past Issue Archive
    • Video Articles
    • AJNR Case Collection
    • Case of the Week Archive
    • Case of the Month Archive
    • Classic Case Archive
  • Special Collections
    • AJNR Awards
    • Low-Field MRI
    • Alzheimer Disease
    • ASNR Foundation Special Collection
    • Photon-Counting CT
    • View All
  • Multimedia
    • AJNR Podcasts
    • AJNR SCANtastic
    • Trainee Corner
    • MRI Safety Corner
    • Imaging Protocols
  • For Authors
    • Submit a Manuscript
    • Submit a Video Article
    • Submit an eLetter to the Editor/Response
    • Manuscript Submission Guidelines
    • Statistical Tips
    • Fast Publishing of Accepted Manuscripts
    • Graphical Abstract Preparation
    • Imaging Protocol Submission
    • Author Policies
  • About Us
    • About AJNR
    • Editorial Board
    • Editorial Board Alumni
  • More
    • Become a Reviewer/Academy of Reviewers
    • Subscribers
    • Permissions
    • Alerts
    • Feedback
    • Advertisers
    • ASNR Home

User menu

  • Alerts
  • Log in

Search

  • Advanced search
American Journal of Neuroradiology
American Journal of Neuroradiology

American Journal of Neuroradiology

ASHNR American Society of Functional Neuroradiology ASHNR American Society of Pediatric Neuroradiology ASSR
  • Alerts
  • Log in

Advanced Search

  • Home
  • Content
    • Current Issue
    • Accepted Manuscripts
    • Article Preview
    • Past Issue Archive
    • Video Articles
    • AJNR Case Collection
    • Case of the Week Archive
    • Case of the Month Archive
    • Classic Case Archive
  • Special Collections
    • AJNR Awards
    • Low-Field MRI
    • Alzheimer Disease
    • ASNR Foundation Special Collection
    • Photon-Counting CT
    • View All
  • Multimedia
    • AJNR Podcasts
    • AJNR SCANtastic
    • Trainee Corner
    • MRI Safety Corner
    • Imaging Protocols
  • For Authors
    • Submit a Manuscript
    • Submit a Video Article
    • Submit an eLetter to the Editor/Response
    • Manuscript Submission Guidelines
    • Statistical Tips
    • Fast Publishing of Accepted Manuscripts
    • Graphical Abstract Preparation
    • Imaging Protocol Submission
    • Author Policies
  • About Us
    • About AJNR
    • Editorial Board
    • Editorial Board Alumni
  • More
    • Become a Reviewer/Academy of Reviewers
    • Subscribers
    • Permissions
    • Alerts
    • Feedback
    • Advertisers
    • ASNR Home
  • Follow AJNR on Twitter
  • Visit AJNR on Facebook
  • Follow AJNR on Instagram
  • Join AJNR on LinkedIn
  • RSS Feeds

AJNR Awards, New Junior Editors, and more. Read the latest AJNR updates

Research ArticleBRAIN

Mapping Functionally Related Regions of Brain with Functional Connectivity MR Imaging

Dietmar Cordes, Victor M. Haughton, Konstantinos Arfanakis, Gary J. Wendt, Patrick A. Turski, Chad H. Moritz, Michelle A. Quigley and M. Elizabeth Meyerand
American Journal of Neuroradiology October 2000, 21 (9) 1636-1644;
Dietmar Cordes
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Victor M. Haughton
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Konstantinos Arfanakis
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gary J. Wendt
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Patrick A. Turski
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Chad H. Moritz
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michelle A. Quigley
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
M. Elizabeth Meyerand
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • Responses
  • References
  • PDF
Loading

Abstract

BACKGROUND AND PURPOSE: In subjects who are performing no prescribed cognitive task, functional connectivity mapped with MR imaging (fcMRI) shows regions with synchronous fluctuations of cerebral blood flow. When specific tasks are performed, functional MR imaging (fMRI) can map locations in which regional cerebral blood flow increases synchronously with the performance of the task. We tested the hypothesis that fcMRI maps, based on the synchrony of low-frequency blood flow fluctuations, identify brain regions that show activation on fMRI maps of sensorimotor, visual, language, and auditory tasks.

METHODS: In four volunteers, task-activation fMRI and functional connectivity (resting-state) fcMRI data were acquired. A small region of interest (in an area that showed maximal task activation) was chosen, and the correlation coefficient of the corresponding resting-state signal with the signal of all other voxels in the resting data set was calculated. The correlation coefficient was decomposed into frequency components and its distribution determined for each fcMRI map. The fcMRI maps were compared with the fMRI maps.

RESULTS: For each task, fcMRI maps based on one to four seed voxel(s) produced clusters of voxels in regions of eloquent cortex. For each fMRI map a closely corresponding fcMRI map was obtained. The frequencies that predominated in the cross-correlation coefficients for the functionally related regions were below 0.1 Hz.

CONCLUSION: Functionally related brain regions can be identified by means of their synchronous slow fluctuations in signal intensity. Such blood flow synchrony can be detected in sensorimotor areas, expressive and receptive language regions, and the visual cortex by fcMRI. Regions identified by the slow synchronous fluctuations are similar to those activated by motor, language, or visual tasks.

Slow fluctuations of regional cerebral blood flow characterize the normal mammalian brain when it is performing no prescribed task (1, 2). In functionally related regions of the brain, even in those regions separated by substantial distances, these fluctuations are synchronous. The synchrony of the blood flow fluctuations in functionally related brain regions implies the existence of neuronal connections that facilitate coordinated activity. To describe this interdependence of functionally related brain regions, the term functional connectivity has been used.

While functional connectivity has been studied by means of electrodes placed in the brain (3–10), by radioisotopes administered intravenously (11–13), and by EEG (14), functional connectivity can now be studied with functional MR imaging (fMRI). Functionally related regions of the brain have been identified by means of their synchronous fluctuations on fMRI studies (15–17). The highly correlated fluctuations in signal intensity in fMRI reflect synchronous changes in the blood oxygen–level dependent (BOLD) effect in functionally related brain regions. Studies with functional connectivity MR imaging (fcMRI) provide a new means of measuring the interdependency of brain regions.

fcMRI uses rapidly acquired images to map the location of changing signal intensity in the brain, as conventional fMRI uses the same technique to map the location of increased blood flow when a stimulus is applied or a cognitive task is performed. In fMRI, signal intensity changes temporally correlated with the performance of a task are identified to localize hemodynamic changes resulting from the task or the stimulus. In fcMRI, regions of synchronous blood flow change are detected to map functionally connected regions. Synchronous fluctuations in signal intensity in the brain are used to examine the strength of neural connections between different brain regions while they are not activated by a specific task. For the fcMRI study, the subject refrains from any cognitive activity as much as possible and receives no experimental stimuli from the investigator. Studies with fcMRI have been used to show connectivity in motor, auditory, and visual cortices in the brain (13, 15–18). Data acquisition using fcMRI provides a means to obtain functional connectivity noninvasively in the human brain.

While the feasibility of fcMRI has been demonstrated, methods are under investigation and their reliability is unknown. Functional connectivity measurements may be subject to artifacts resulting from bulk motion and fluctuations (or aliases) in signal intensity due to cardiac or respiratory cycles. Also, the resting state cannot be specifically defined or controlled during data acquisition and might be quite different during different scan sessions in the same subject.

The purpose of this study was to obtain fcMRI maps based on slow periodic fluctuation in signal intensity and compare them with fMRI maps of activation. For this purpose we derived an expression for the spectral distribution of the cross-correlation coefficient. The spectral decomposition of the correlation coefficient showed which frequencies contribute to interregional connectivity and that only frequencies less than 0.1 Hz contribute significantly to interregional functional connectivity. Although these results may have been assumed, they have not to our knowledge been reported for fcMRI. We compared functional connectivity MR imaging (fcMRI) maps with functional MR imaging (fMRI) activation maps for sensorimotor, language, auditory, and visual tasks in a series of volunteers.

Methods

Four volunteers between the ages of 20 and 40 years, claiming to be in good health with no history of neurologic illness, were enrolled after consent was obtained according to institutional regulations. Using a 1.5-T commercial MR scanner equipped with high-speed gradients, we obtained high-resolution whole-brain anatomic images with 2D multislice, spin-echo sequences in the coronal plane with the following parameters: 500/8 (TR/TE), 18 slices, 256 × 128 matrix, 24-cm field of view (FOV), 7-mm section thickness with a 2-mm gap. In addition, a 3D spoiled gradient-recalled whole-brain volume scan was acquired in the axial plane with the following parameters: 21/7, 124 slices, 256 × 128 matrix, 40° flip angle, 24-cm FOV, 1.2-mm section thickness with no gap. The subject's head was positioned within a prototype quadrature bird-cage coil specifically developed for functional imaging of the brain. Foam padding was provided for comfort and to minimize head movements. For fMRI, a series of 148 time frames was acquired with single-shot gradient-recalled echo-planar imaging in the coronal plane. For fcMRI, the number of time frames was increased to 260. In each series, the images belonging to the first four time frames were discarded because the MR signal was not stationary. Technical parameters for these acquisitions were as follows: 2000/50, 18 slices, 64 × 64 matrix, 90° flip angle, 24-cm FOV, 125-kHz receiver bandwidth, 7-mm section thickness with a 2-mm gap. Voxel dimensions were 3.75 × 3.75 × 7 mm (x,y,z). The cardiac and respiratory rates were monitored using a pulse oximeter and a flexible respiratory belt.

For task activation, the subject performed a specific task for four 32-second periods, each separated by 32 seconds of rest. Each volunteer performed each of four tasks: finger tapping, silent word generation, listening to text, and viewing of a screen with an alternating checkerboard pattern. These tasks have been described in detail elsewhere (19–24). In brief, the finger-tapping task was a self-paced apposition of the thumb and other fingers on cue from the investigator. For the word-generation task, the subject was instructed to think of as many words as possible that begin with a specific letter cued by the investigator. For the text-listening task, the subject listened to a recorded text. For visual cortex activation, the subject viewed a screen on which an alternating checkerboard pattern was presented with a frequency of 8 Hz. During fcMRI, the subject was instructed to rest, be as motionless as possible, keep the eyes closed, and perform no specific cognitive exercise during the entire acquisition (resting acquisition). fcMRI was performed after fMRI.

The echo-planar raw data were filtered in the spatial frequency domain with a low-pass Hanning filter to increase signal-to-noise ratio (25). Signal intensities for each acquired slice were time-corrected by using a three-point Hanning filter in the temporal domain (26). Images were checked for motion using a 3D registration program included in automated functional neuroimaging (AFNI; Robert Cox, Medical College of Wisconsin). The study was discarded and repeated if the detected motion gave a displacement larger than 0.4 mm in x, y, or z. Two studies were repeated because of motion.

For fMRI, the time course was analyzed with a locally developed program (27, 28). A least-squares fitting algorithm was used to compare the observed data on a voxel-by-voxel basis to a constant (baseline signal intensity), a ramp (signal drift), and a boxcar function (idealized expected response to the task or stimulus). The boxcar function had unit amplitude and a period matching the on/off cycles of the task or stimulus that was smoothed by convolution with a Poisson function with a 6-second mean to account for hemodynamic delay. The amplitude and the uncertainty of the fit to the boxcar reference function were calculated by least-squares minimization. A t-score was calculated and converted to a z-score. An overall adjusted P value was determined to address levels of significance. All voxels with P less than .05 were merged with coregistered anatomic images.

For fcMRI, “seed” voxels were selected in the rolandic cortex, primary auditory cortex, Broca's area, and striate cortex by referring to the maximum activation seen on the fMRI study. For comparison, several other seed voxels were selected in other regions where the cardiac and respiratory effects were expected to be different. These regions were in the CSF (left lateral ventricle), in a branch of the middle cerebral artery (left sylvian fissure), and in the left jugular vein. The correlation coefficients of the seed voxels and every other voxel in the brain were calculated.

To compute the statistical significance of fcMRI maps, a nonparametric multiple comparisons procedure (29, 30) was designed. This procedure maintained strong control over type 1 errors. To arrive at a null statistic for resting-state data, the time course of the seed voxel or cluster was altered by randomization of the phase of its complex Fourier components, an approach that does not change the power spectrum. The null distribution of correlation coefficients for a particular seed voxel or cluster was then obtained by calculating all possible cross-correlation coefficients between the seed voxel after phase randomization and all other voxels in the brain. This approach gave an assignment of overall adjusted P values to each cross-correlation value when carried out for each seed voxel or cluster. All fcMRI maps were then thresholded at the P = .05 level. The corresponding cross-correlation coefficients (using 256 time frames) were in the range of .35 to .55, depending on the seed voxel or cluster. The degrees of freedom of the seed voxel or cluster as computed by Friston's method (31) were in the range of 10 to 120. Voxels correlating significantly were merged with anatomic images. The locations of the voxels on the fcMRI maps were compared with the location of task activation on the fMRI maps.

The average frequency spectrum of the correlation coefficients was computed for each fcMRI map according to Equation 17 in the Appendix. Briefly, Equation 17 was derived by expressing the voxel time courses in terms of their inverse Fourier transforms and calculating the cross-correlation coefficient in frequency (ie, Fourier) space. The resulting expression defines the average magnitude of each frequency component in an fcMRI map. The so-defined spectral distribution of the correlation coefficient has the property that the sum over all frequencies will yield the average correlation coefficient in an fcMRI map.

Results

Technically adequate task-activation maps were obtained for all volunteers. Activation was identified in the sensorimotor cortices, the supplementary motor areas, the dentate nucleus (finger tapping), the superior temporal lobe (text listening), the left inferior frontal gyrus (word generation), and the striate cortex (screen viewing). Other regions of activation were identified in some subjects.

In each subject, the fcMRI maps based on a seed voxel from each of four eloquent brain locations resembled the task-activation maps for the region. On average, the number of voxels showing functional connectivity in the resting acquisition was more than the number of voxels showing task activation at the P = .05 significance threshold. In general, more slices had voxels that passed the threshold on the fcMRI maps than on the fMRI maps. Some voxels passing threshold on the fcMRI maps were imperfectly correlated with task-activation fMRI maps, but similar regions were denoted by both methods.

For the seed voxels selected in the sensorimotor cortex, significantly correlated voxels were found in the contralateral and ipsilateral sensorimotor cortex regions. For these regions, the fcMRI map in the resting acquisition resembled the fMRI map of activation produced by the finger-tapping task. Figure 1 shows a typical example of task activation and functional connectivity in the sensorimotor cortex. Functional connectivity with the dentate nucleus in the cerebellum was not observed on the fcMRI map. For a seed voxel selected in the superior temporal gyrus in the region of the auditory cortex, significantly correlated voxels were present in the resting acquisition in the ipsilateral and contralateral auditory cortex regions (Fig 2B), and the locations were similar to the regions activated by text listening (Fig 2A). Figure 3 shows frontal lobe activation in or near Broca's area in a subject performing the text-listening task. An fcMRI map based on a pixel from the activated regions produced a map of functional connectivity that resembled the task-activation map. An fMRI map showing activation in or near Broca's area in a subject performing the word-generation task is shown in Figure 4A. A seed voxel selected in the resting data set produced a functional connectivity map that resembled the task-activation map (Fig 4B). Task-activation and fcMRI maps of the visual cortices also compared well (Fig 5). The checkerboard viewing task produced activation in the striate cortex and related regions (Fig 5A). The fcMRI map, based on a pixel in the striate cortex, produced a similar map (Fig 5B).

fig 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
fig 1.

A and B, Selected coronal and axial fMRI maps of voxels activated by the bilateral finger-tapping task (A) and of the voxels functionally connected to a seed voxel (crosshairs) in the sensorimotor cortex in a resting acquisition (B). The sensorimotor cortex is identified on both the fMRI and the fcMRI maps. Voxels within the dentate nucleus of the cerebellum were identified with fMRI but not with the functional connectivity study

fig 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
fig 2.

A and B, Consecutive coronal fMRI maps of voxels activated by a text-listening task (A) and of the voxels functionally connected with a 2 × 2 seed voxel ROI (crosshairs) within the auditory cortex (B). Task activation was evident in the region of the primary and association auditory cortices in the superior temporal lobe. Many of the voxels functionally connected with the seed voxel in the left auditory cortex (crosshairs) have a similar distribution. A few other voxels with lower correlation coefficients (red) were identified outside the auditory cortices. For the connectivity study, data were collected without the patient performing a specific cognitive task

fig 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
fig 2.

fig 3. A and B, A coronal image showing task activation for the text-listening task (A) and connectivity for a seed voxel (crosshairs) in the left frontal lobe in a resting data set (B). The task produces activation in the superior temporal gyral regions bilaterally and in the left frontal lobe in or near Broca's region. Activation is also identified in the right frontal lobe. The voxels with connectivity to the seed voxel (crosshairs) appear to lie in similar locations. fig 4. A and B, Coronal images showing task activation for the word-generation task (A) and connectivity for the left frontal cortex in the region that showed task activation (B). The task produced activation in the left frontal lobe in or near Broca's area and in the superior temporal lobes bilaterally. The voxels with connectivity to the seed voxel (crosshairs) appear to lie in similar locations to the task activation produced by word generation and to the locations with connectivity in figure 3

fig 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
fig 5.

A and B, Consecutive coronal images showing task activation for the visual task (A) and connectivity for the visual cortex (B). The task produces activation in the striate cortex region and in the posterior parietal lobes bilaterally. The voxels with connectivity to the seed voxel (crosshairs) appear to lie in similar locations in the striate cortex region and in the posterior parietal lobe

Not every seed voxel selected in a region produced the same map. Some seed voxels, selected from voxels showing considerable activation, produced fcMRI maps with less resemblance to the task-activation map, whereas adjacent voxels produced maps closely resembling the activation map. In some cases, averaging several voxels to produce a time course for cross correlation slightly improved consistency. An fcMRI map closely resembling a task-activation map was obtained when the investigator selected among voxels in regions with maximum task activation. Voxels selected randomly outside regions that were activated by a task did not, in general, show significant correlation with voxels in the sensorimotor, auditory, language, and visual regions.

The spectral decomposition analysis for the functional connectivity showed predominantly low frequency contributing to the correlation coefficient in functionally related brain regions (Figs 6–10). The recorded respiratory rate for all volunteers (below 15 breaths per minute for the resting state study with an average of 12 breaths per minute) did not alias into the lower frequency range at TR = 2000. For seed voxels taken from regions containing a major blood vessel or CSF, the spectral distributions of the correlation maps did not display the same predominance of low frequencies. In general, tissue near arteries showed more power in frequencies related to the cardiac cycle whereas tissue near veins showed more evidence of frequencies related to the respiratory cycle (32). For example, a seed voxel in the jugular vein showed a large contribution at a frequency of 0.2 Hz; otherwise, the frequency distribution was about equal across the full range (Fig 11). For a seed voxel from the middle cerebral artery branch, several dominant peaks were present across the full range. In the low-frequency domain (<0.1 Hz), the major contribution came from a peak at 0.002 Hz, whereas in the respiratory frequency range (0.15–0.25 Hz), some contributions from bands centered at 0.19 and 0.24 Hz appeared (Fig 12). A seed voxel representing CSF in the lateral ventricle showed major contributions in the low-frequency range (<0.1 Hz) but also some contributions in the interval between 0.21 and 0.25 Hz (Fig 13).

fig 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
fig 6.

Spectral decomposition of the average correlation coefficient for interregional connectivity in the motor cortex (same subject as in fig 1). Only frequency components between 0 and 0.05 Hz contribute significantly to the correlation coefficient

fig 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
fig 7.

Spectral decomposition of the average correlation coefficient for interregional connectivity in the auditory cortex (same subject as in fig 2). Only frequency components less than 0.05 Hz contribute significantly to the correlation coefficient

fig 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
fig 8.

Spectral decomposition of the average correlation coefficient for interregional connectivity in the prefrontal cortex (same subject as in fig 3). Only the same low frequency components contribute significantly to the correlation coefficient

fig 9.
  • Download figure
  • Open in new tab
  • Download powerpoint
fig 9.

Spectral decomposition of the average correlation coefficient for interregional connectivity in the prefrontal cortex (same subject as in fig 4). Only the same low-frequency components contribute significantly to the correlation coefficient

fig 10.
  • Download figure
  • Open in new tab
  • Download powerpoint
fig 10.

Spectral decomposition of the average correlation coefficient for interregional connectivity in the visual cortex (same subject as in fig 5). Only low-frequency components contribute significantly to the correlation coefficient

fig 11.
  • Download figure
  • Open in new tab
  • Download powerpoint
fig 11.

Spectral decomposition of the average correlation coefficient from a seed voxel in the left jugular vein. Note the relative paucity of frequencies in the 0 to 0.05 Hz range, the spread over many frequencies, and the peak at 0.2 Hz

fig 12.
  • Download figure
  • Open in new tab
  • Download powerpoint
fig 12.

Spectral decomposition of the average correlation coefficient using a seed voxel in the right middle cerebral artery. Low frequencies do not predominate. Peaks are present at multiple frequencies above 0.05 Hz

fig 13.
  • Download figure
  • Open in new tab
  • Download powerpoint
fig 13.

Spectral decomposition of the average correlation coefficient using a seed voxel in the left ventricle. Frequencies from 0.025 to 0.075 predominate. Since no BOLD effect can be invoked in CSF, these frequencies represent likely aliased cardiac cycle effects

Discussion

In previous studies, functional connectivity was detected in the sensorimotor cortex (15, 16, 18–20, 22), the auditory cortex/language regions (21), and the visual cortex (23). This study shows that by using only one voxel or a small region of interest (ROI) (four voxels or less) in a region of eloquent cortex, functionally related brain regions can be identified by means of their synchronous slow fluctuations in blood flow in the absence of a prescribed cognitive task. Such blood flow synchrony can be detected in sensorimotor areas, expressive language regions, receptive language regions, and visual cortex by fcMRI.

fcMRI is distinguished from fMRI, which depicts blood-flow changes (activation) temporally correlated with the performance of a task. Tasks such as finger tapping, text listening, word generation, or picture viewing typically produce activation in the sensorimotor cortices, auditory cortices, language cortex, and visual cortex, respectively. The regions identified with fMRI as activated correlate well with regions identified as having functional significance with intraoperative mapping or positron emission tomography. However, the fcMRI and fMRI maps do not correlate perfectly. The regions of brain that have functional connectivity are not necessarily the same as the regions activated by the performance of a task. fcMRI is not a replacement for fMRI. That fMRI and fcMRI identify similar regions of brain has both investigative and practical significance.

We applied a novel, nonparametric statistical method to the fcMRI data in order not to assume a specific probability distribution for the voxel values in the image. In particular, a null hypothesis is not as easily defined for fcMRI maps as for task-activation data. Parametric approaches based on gaussian random fields, which are valid for the analysis of task-activation data, cannot be applied to fcMRI maps, because each seed voxel chosen shows a large variability in autocorrelations of its signal time course. For example, the seed voxels in our study had 10 to 120 degrees of freedom, as calculated by the method of Friston (31), for a series of 256 time frames. To define a proper null hypothesis for fcMRI as a measure of signal randomness or signal nonsynchronicity, our approach was to construct first a function from each seed voxel time course by randomizing the phase of its complex Fourier components. The null distribution of correlation coefficients for a particular seed voxel or cluster was then obtained by calculating all possible cross-correlation coefficients between the seed voxel (or cluster) after phase randomization with all other voxels in the brain. This method retained the frequency content of the seed voxel, since varying only the phase of its complex Fourier components did not change the power spectrum. The obtained P values were very reasonable. Equally thresholded fcMRI and fMRI maps had similar intensity distributions. Our nonparametric method maintained strong control over type 1 errors, meaning that voxels that lacked functional connectivity were detected with a probability, at most, of P = .05 in any region of the brain, regardless of connectivity elsewhere.

The decomposition of the cross-correlation maps into frequency components reveals several interesting facts. Only frequencies less than 0.1 Hz contributed significantly to the correlation coefficient in functionally connected brain regions. This finding is consistent with other reports (15, 16) in which low-frequency fluctuations (<0.1 Hz) in functionally related brain regions showed a high degree of temporal correlation. Because our TR was 2000, it is uncertain how much the aliased cardiac rate contributed to the connectivity. A spectrum for CSF shows several dominant peaks in the low-frequency range, which are very likely due to the cardiac rate. Similarly, in the spectrum of the middle cerebral artery, there is some contribution from frequencies less than 0.1 Hz, possibly due to cardiac aliasing. However, to filter out the cardiac contribution is not trivial because of its bandlike appearance and unknown region of aliasing. The spread of the cardiac rate occurs as a result of normal changes in cardiac frequency during the scan and depends on the resting state of the volunteer. In our studies, the cardiac rate, monitored with a pulse oximeter, varied by approximately three beats per minute, leading to a frequency band rather than a single peak. To estimate the contribution of the cardiac rate to the correlation map, it is necessary to use a short TR so that the cardiac rate does not alias and can easily be identified. A short TR, however, limits the study to a single slice (or a few slices).

On the other hand, at TR = 2000, it is possible to rule out respiratory effects on interregional connectivity. The observed respiratory rate, limited to a narrow region around 0.2 Hz, did not alias in our study. Based on the calculated cross-correlation spectrum, only CSF and blood vessels showed contributions in the respiratory range. We did not find any significant amount of respiratory frequencies contributing to interregional connectivity of motor, language, and visual cortices. The application of finite impulse response filters to exclude respiratory contaminations in the frequency range above 0.1 Hz did not quantify the effect of respiratory aliasing on the correlation map. On the basis of our finding, the correlation map for interregional connectivity appears to be relatively insensitive to the respiratory cycle.

This preliminary result cannot be generalized without further studies. A weakness of the study is the selection bias in choosing a seed voxel. If the seed voxels partially overlap with CSF or a blood vessel, results can vary substantially between voxels selected from within a small region in eloquent cortex. In our study, the cross-correlation analysis was not tabulated for every voxel in the regions that showed significant activation on task-acquired fMRI maps. Whether functional connectivity can be measured consistently and routinely in patients with neurologic disorders is under investigation. Also, individual variation and the physiology of the functional connectivity in the human brain have not yet been determined. The significance of regions that are identified in one and not the other mapping technique needs more study.

Nevertheless, our method may be useful for evaluating the strength of particular neuronal connections in healthy individuals and in patients with neurologic disorders. The fluctuations in signal intensity, on which fcMRI maps are based, are affected by blood oxygen and carbon dioxide concentrations (33). Their application to neurologic disorders that are thought to stem from disconnection syndromes may be especially important (34).

Conclusion

We found a high degree of temporal correlation of physiological fluctuations between functionally related neuronal sites in subjects who were performing no specified cognitive task. Technically adequate correlation maps showing interregional connectivity in the motor cortex, the expressive and receptive language areas, and the visual cortex were found by selecting only a small ROI (one to four voxels) in eloquent cortex as a probe for functional connectivity. Our observations indicate that low-frequency fluctuations are a general phenomenon in eloquent cortex; therefore, fcMRI may be a useful method for evaluating the strength of neuronal connections in healthy individuals and in patients with neurologic disorders.

Appendix: Derivation of the Spectral Decomposition of the Correlation Coefficient

The correlation coefficient can be decomposed into frequency components by Fourier transformation of the voxel time courses. Let f(t) and g(t) represent the signal intensity of two voxels as a function of time, defined over the interval 0, N. Furthermore, the mean signal of each function has been subtracted. Using the inverse Fourier transform, we can write Embedded Image where N−1 refers to the number of frequency components and ωn, λn to the complex-valued frequencies of the Fourier transform of f(t) and g(t), respectively. Note, that ω0 = λ0 = 0 in our case. The correlation coefficient, cc, involving f(t) and g(t) is defined as Embedded Image where the denominator, D, is given by Embedded Image

Inserting Equations 1 and 2 into Equation 3 yields Embedded Image

Rearranging the order of summation and combining terms gives Embedded Image The sum over t simplifies owing to the orthogonality condition Embedded Image where δi,j is the Kronecker delta symbol and m can be any integer. Since in our case Embedded Image m = 1 is the only possibility in Equation 7 for a nonvanishing contribution. Hence, we can write Equation 6 as Embedded Image which simplifies to Embedded Image Since the original voxel-time courses are real-valued functions, the complex Fourier components ωn, λn satisfy the symmetry condition Embedded Image The expression for the correlation coefficient in Equation 10 becomes Embedded Image and hence Embedded Image We are now able to define the spectral decomposition of the correlation coefficient by Embedded Image

In the following we include superscripts to label the two voxels that are correlated. For an ROI, the average spectral decomposition of the correlation coefficient between a seed voxel, s, and all other voxels in the brain, j, can be defined by Embedded Image where the upper indexes s and j refer to the seed voxel and to all other voxels on the fcMRI map, respectively. The index j labels all voxels in the brain that have significant correlation; that is, Embedded Image for all j voxels on the fcMRI map. The term thr represents a threshold at the appropriate significance level.

Footnotes

  • ↵1 Address reprint requests to Victor M. Haughton, MD, Department of Radiology, The University of Wisconsin Hospitals and Clinics, 600 Highland Ave, Madison, WI 53792.

References

  1. ↵
    Haughton VM, Biswal B. Clinical application of basal regional cerebral blood flow fluctuation measurements by FMRI. In: Hudetz A, Bruley D, eds: Oxygen Transport to Tissue. New York: Plenum 1998;:583-590
  2. Hudetz AG. Blood flow in the cerebral capillary network: a review emphasizing observations with intravital microscopy. Microcirculation 1997;4:233-252
    PubMed
  3. ↵
    Dora E, Kovach AGB. Metabolic and vascular volume oscillation in the cat brain cortex. Acta Physiol Hung 1980;57:261-275
  4. Fasano VA, Urciuoli R, Bolognese P, Mostert M. Introperative use of laser Doppler in the study of cerebral microvascular circulation. Acta Neurochir (Wien) 1988;95:40-48
    CrossRefPubMed
  5. Gray CM, Konig P, Engel AK, Singer W. Oscillatory responses in cat visual cortex exhibit inter-columnar synchronization which reflects global stimulus properties. Nature 1989;338:334-337
    CrossRefPubMed
  6. Hudetz AG, Roman RJ, Harder DR. Spontaneous flow oscillations in the cerebral cortex during acute changes in mean arterial pressure. J Cereb Blood Flow Metab 1992;12:491-499
    PubMed
  7. Mayevsky A, Ziv I. Oscillations of cortical oxidative metabolism and microcirculation in the ischaemic brain. Neurol Res 1991;13:39-47
    PubMed
  8. Mitra PP, Ogawa S, Hu X, Ugurbil K. The nature of spatiotemporal changes in cerebral hemodynamics as manifested in functional magnetic resonance imaging. Magn Reson Med 1997;37:511-518
    PubMed
  9. Tomita M, Gotoh F, Sato T, Tanahashi N, Tanaka K. 4–6 cycle per minute fluctuation in cerebral blood volume of feline cortical tissue in situ. J Cereb Blood Flow Metab 1981;1::S443-S444
  10. Vern BA, Schuette WH, Leheta B, Juel VC, Radulovacki M. Low frequency oscillations of cortical oxidative metabolism in waking and sleep. J Cereb Blood Flow Metab 1998;8:215-226
  11. ↵
    Friston KJ, Frith CD, Liddle PF, Frackowiak RSJ. Functional connectivity: the principal-component analysis of large (PET) data sets. J Cereb Blood Flow Metab 1993;13:5-14
    PubMed
  12. Gerstein GL, Perkel DH, Subramanian KN. Identification of functionally related neural assemblies. Brain Res 1978;140:43-62
    CrossRefPubMed
  13. ↵
    Horwitz B. Functional interactions in the brain: use of correlation between regional metabolic rates. J Cereb Blood Flow Metab 1991;11:A114-A120
    PubMed
  14. ↵
    Wada Y, Nanbu Y, Jiang ZY, Koshino Y, Hashimoto T. Interhemispheric EEG coherence in never-medicated patients with paranoid schizophrenia: analysis at rest and during photic stimulation. Clin Electroencephalogr 1998;29:170-176
    PubMed
  15. ↵
    Biswal B, Yetkin FZ, Haughton VM, Hyde JS. Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn Reson Med 1995;34:537-541
    PubMed
  16. Lowe MJ, Mock BJ, Sorenson JA. Functional connectivity in single and multislice echoplanar imaging using resting state fluctuations. Neuroimage 1998;7:119-132
    CrossRefPubMed
  17. Lowe MJ, Mock BJ, Sorenson JA. Resting state fMRI signal correlation in multi-slice. Neuroimage 1996;3:S257
  18. Xiong J, Parsons LM, Gao JH, Fox PT. Interregional connectivity to primary motor cortex revealed using MRI resting state images. Hum Brain Mapping 1999;8:151-156
    CrossRefPubMed
  19. ↵
    Hammeke TA, Yetkin FZ, Mueller WM, Morris GL III, Haughton VM. Functional magnetic resonance imaging of somatosensory stimulation. Neurosurgery 1994;35:677-681
    PubMed
  20. Mueller WM, Yetkin FZ, Hammeke TA, et al. Functional magnetic resonance mapping of the motor cortex in patients with cerebral tumors. Neurosurgery 1996;39:515-521
    CrossRefPubMed
  21. ↵
    Yetkin FZ, Hammeke TA, Swanson SJ, et al. A comparison of functional MR activation patterns during silent and audible language tasks. AJNR Am J Neuroradiol 1995;16:1087-1092
    Abstract
  22. Yetkin FZ, Mueller WM, Hammeke TA, Morris GL III, Haughton VM. Functional magnetic resonance image mapping of the sensorimotor cortex with tactile stimulation. Neurosurgery 1995;36:921-925
    PubMed
  23. ↵
    DeYoe EA, Bandettini P, Neitz J, Miller D, Winans P. Functional magnetic resonance imaging (FMRI) of the human brain. J Neurosci Methods 1994;54:171-187
    CrossRefPubMed
  24. Turski PA, Cordes D, Mock B, et al. Basic concepts of functional magnetic resonance imaging and arteriovenous malformations. Neuroimaging Clin N Am 1998;8:371-381
    PubMed
  25. ↵
    Lowe MJ, Sorenson JA. Spatially filtering functional magnetic resonance imaging data. Magn Reson Med 1997;37:723-729
    PubMed
  26. ↵
    McKeown MJ, Makeig S, Brown GG, et al. Analysis of fMRI data by blind separation into independent spatial components. Hum Brain Mapping 1988;6:160-188
  27. ↵
    Russel DP. A generalized approach to time-course analysis of functional MRI of the human brain. In: Proceeding of the 2nd Annual Meeting of the Society of Magnetic Resonance, San Francisco, 1994. Society of Magnetic Resonance 1994;:636
  28. Lowe MJ, Russel DP. Treatment of baseline drifts in fMRI time series analysis. J Comput Assist Tomogr 1999;23:463-473
    CrossRefPubMed
  29. ↵
    Sprent P. Applied Nonparametric Statistical Methods. London: Chapman and Hall 1993;
  30. Holmes AP, Blair RC, Watson DG, Ford I. Nonparametric analysis of statistic images from functional mapping experiments. J Cereb Blood Flow Metab 1996;16:7-22
    CrossRefPubMed
  31. ↵
    Friston KJ, Holmes AP, Poline JB, et al. Analysis of fMRI time-series revisited. Neuroimage 1995;2:45-53
    CrossRefPubMed
  32. ↵
    Lowe MJ. Gram-Schmidt orthogonalization to reduce aliased physiologic noise in low sampling rate fMRI data. In: Proceedings of the 7th Annual Meeting of the International Society of Magnetic Resonance in Medicine, Philadelphia, 1999. International Society of Magnetic Resonance in Medicine 1999;:1711
  33. ↵
    Biswal B, Hudetz AG, Yetkin FZ, Haughton VM, Hyde JS. Hypercapnia reversibly suppresses low-frequency fluctuations in the human motor cortex during rest using echo-planar MRI. J Cereb Blood Flow Metab 1997;17:301-308
    CrossRefPubMed
  34. ↵
    Geschwind N. Disconnection syndromes in animals and man. Brain 1965;88:231-294
  • Received July 6, 1999.
  • Accepted after revision April 11, 2000.
  • Copyright © American Society of Neuroradiology
View Abstract
PreviousNext
Back to top

In this issue

American Journal of Neuroradiology
Vol. 21, Issue 9
1 Oct 2000
  • Table of Contents
  • Index by author
Advertisement
Print
Download PDF
Email Article

Thank you for your interest in spreading the word on American Journal of Neuroradiology.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Mapping Functionally Related Regions of Brain with Functional Connectivity MR Imaging
(Your Name) has sent you a message from American Journal of Neuroradiology
(Your Name) thought you would like to see the American Journal of Neuroradiology web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Cite this article
Dietmar Cordes, Victor M. Haughton, Konstantinos Arfanakis, Gary J. Wendt, Patrick A. Turski, Chad H. Moritz, Michelle A. Quigley, M. Elizabeth Meyerand
Mapping Functionally Related Regions of Brain with Functional Connectivity MR Imaging
American Journal of Neuroradiology Oct 2000, 21 (9) 1636-1644;

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
0 Responses
Respond to this article
Share
Bookmark this article
Mapping Functionally Related Regions of Brain with Functional Connectivity MR Imaging
Dietmar Cordes, Victor M. Haughton, Konstantinos Arfanakis, Gary J. Wendt, Patrick A. Turski, Chad H. Moritz, Michelle A. Quigley, M. Elizabeth Meyerand
American Journal of Neuroradiology Oct 2000, 21 (9) 1636-1644;
del.icio.us logo Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Methods
    • Results
    • Discussion
    • Conclusion
    • Appendix: Derivation of the Spectral Decomposition of the Correlation Coefficient
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • Responses
  • References
  • PDF

Related Articles

  • No related articles found.
  • PubMed
  • Google Scholar

Cited By...

  • Dynamic fusion of structural and functional connectivity via joint connectivity matrix ICA
  • Aging and the Spectral Properties of Brain Hemodynamics
  • Prediction of dynamic balance state and recovery following stroke using fMRI graph analysis
  • Unraveling the Neural Landscape of Mental Disorders using Double Functional Independent Primitives (dFIPs)
  • Functional brain network dynamics of brooding in depression: insights from real-time fMRI neurofeedback
  • Altered dynamic functional network connectivity in rheumatoid arthritis associated with peripheral inflammation and neuropsychiatric disorders
  • Functional connectivity across the human subcortical auditory system using an autoregressive matrix-Gaussian copula graphical model approach with partial correlations
  • Resting-state functional connectivity in children cooled for neonatal encephalopathy
  • How can graph theory inform the dual-stream model of speech processing? a resting-state fMRI study of post-stroke aphasia
  • Efficient estimation of time-dependent functional connectivity using Structural Connectivity constraints
  • Single Ventricle Reconstruction III: Brain Connectome and Neurodevelopmental Outcomes: Design, Recruitment, and Technical Challenges of a Multicenter, Observational Neuroimaging Study
  • Resting state networks of awake adolescent and adult squirrel monkeys using ultra-high field (9.4T) functional magnetic resonance imaging
  • Efficient estimation of time-dependent functional connectivity using Structural Connectivity constraints
  • Joint cmICA: auto-linking of structural and functional connectivity
  • Acquisition and processing methods of whole-brain layer-fMRI VASO and BOLD: The Kenshu dataset
  • A structure-function substrate of memory for spatial configurations in medial and lateral temporal cortices
  • Resting-State Functional Networks of Different Topographic Representations in the Somatosensory Cortex of Macaque Monkeys and Humans
  • A Paradigm for Measuring Resting State Functional Connectivity in Young Children Using fNIRS and Freeplay
  • MOU-EC: model-based whole-brain effective connectivity to extract biomarkers for brain dynamics from fMRI data and study distributed cognition
  • Trait-like variants in human functional brain networks
  • Studying the precuneus reveals structure-function-affect correlation in long-term meditators
  • The Relationship between BOLD and Neural Activity Arises from Temporally Sparse Events
  • Large-scale intrinsic connectivity is consistent across varying task demands
  • Extracting orthogonal subject- and behavior-specific signatures from fMRI data using whole-brain effective connectivity
  • Analysis of fMRI data using noise-diffusion network models: a new covariance-coding perspective
  • Stereotypical modulations in dynamic functional connectivity explained by changes in BOLD variance
  • Effective connectivity inferred from fMRI transition dynamics during movie viewing points to a balanced reconfiguration of cortical interactions
  • Resting State Networks in empirical and simulated dynamic functional connectivity
  • Impact of gray matter signal regression in resting state and language task functional networks
  • State-dependent modulation of functional connectivity in early blind individuals
  • High spatial correspondence at a columnar level between activation and resting state fMRI signals and local field potentials
  • Brain Age: A State-Of-Mind? On the Stability of Functional Connectivity across Behavioral States
  • Functional connectivity arises from a slow rhythmic mechanism
  • Imperceptible Somatosensory Stimulation Alters Sensorimotor Background Rhythm and Connectivity
  • Noninvasive Functional and Anatomical Imaging of the Human Medial Temporal Lobe
  • Intrinsically organized resting state networks in the human spinal cord
  • Low-frequency calcium oscillations accompany deoxyhemoglobin oscillations in rat somatosensory cortex
  • Time-varying functional network information extracted from brief instances of spontaneous brain activity
  • Reorganization of Functional Connectivity of the Language Network in Patients with Brain Gliomas
  • Unravelling the Intrinsic Functional Organization of the Human Lateral Frontal Cortex: A Parcellation Scheme Based on Resting State fMRI
  • Functional Connectivity during Resting-State Functional MR Imaging: Study of the Correspondence between Independent Component Analysis and Region-of-Interest-Based Methods
  • Prefrontal Transcranial Direct Current Stimulation Changes Connectivity of Resting-State Networks during fMRI
  • Global Functional Connectivity Deficits in Schizophrenia Depend on Behavioral State
  • Comparison of Hypothesis- and a Novel Hybrid Data/Hypothesis-Driven Method of Functional MR Imaging Analysis in Patients with Brain Gliomas
  • Functional Connectivity MR Imaging of the Language Network in Patients with Drug-Resistant Epilepsy
  • Growing Together and Growing Apart: Regional and Sex Differences in the Lifespan Developmental Trajectories of Functional Homotopy
  • Neural basis of global resting-state fMRI activity
  • Investigating the Functional Heterogeneity of the Default Mode Network Using Coordinate-Based Meta-Analytic Modeling
  • Correspondence of the brain's functional architecture during activation and rest
  • Decoupling of the brain's default mode network during deep sleep
  • Distinct Cerebellar Contributions to Intrinsic Connectivity Networks
  • The Blood Oxygen Level-Dependent Functional MR Imaging Signal Can Be Used to Identify Brain Tumors and Distinguish Them from Normal Tissue
  • Regional Variation in Interhemispheric Coordination of Intrinsic Hemodynamic Fluctuations
  • Resting-state networks in the infant brain
  • Dissociable Intrinsic Connectivity Networks for Salience Processing and Executive Control
  • Consistent resting-state networks across healthy subjects
  • Spontaneous neuronal activity distinguishes human dorsal and ventral attention systems
  • From The Cover: The human brain is intrinsically organized into dynamic, anticorrelated functional networks
  • Frequencies Contributing to Functional Connectivity in the Cerebral Cortex in ""Resting-state"" Data
  • Crossref
  • Google Scholar

This article has not yet been cited by articles in journals that are participating in Crossref Cited-by Linking.

More in this TOC Section

  • Statin Therapy Does Not Affect the Radiographic and Clinical Profile of Patients with TIA and Minor Stroke
  • Usefulness of Quantitative Susceptibility Mapping for the Diagnosis of Parkinson Disease
  • White Matter Alterations in the Brains of Patients with Active, Remitted, and Cured Cushing Syndrome: A DTI Study
Show more Brain

Similar Articles

Advertisement

Indexed Content

  • Current Issue
  • Accepted Manuscripts
  • Article Preview
  • Past Issues
  • Editorials
  • Editor's Choice
  • Fellows' Journal Club
  • Letters to the Editor
  • Video Articles

Cases

  • Case Collection
  • Archive - Case of the Week
  • Archive - Case of the Month
  • Archive - Classic Case

More from AJNR

  • Trainee Corner
  • Imaging Protocols
  • MRI Safety Corner
  • Book Reviews

Multimedia

  • AJNR Podcasts
  • AJNR Scantastics

Resources

  • Turnaround Time
  • Submit a Manuscript
  • Submit a Video Article
  • Submit an eLetter to the Editor/Response
  • Manuscript Submission Guidelines
  • Statistical Tips
  • Fast Publishing of Accepted Manuscripts
  • Graphical Abstract Preparation
  • Imaging Protocol Submission
  • Evidence-Based Medicine Level Guide
  • Publishing Checklists
  • Author Policies
  • Become a Reviewer/Academy of Reviewers
  • News and Updates

About Us

  • About AJNR
  • Editorial Board
  • Editorial Board Alumni
  • Alerts
  • Permissions
  • Not an AJNR Subscriber? Join Now
  • Advertise with Us
  • Librarian Resources
  • Feedback
  • Terms and Conditions
  • AJNR Editorial Board Alumni

American Society of Neuroradiology

  • Not an ASNR Member? Join Now

© 2025 by the American Society of Neuroradiology All rights, including for text and data mining, AI training, and similar technologies, are reserved.
Print ISSN: 0195-6108 Online ISSN: 1936-959X

Powered by HighWire