Functional magnetic resonance imaging
Functional magnetic resonance imaging or functional MRI (fMRI) is a type of specialized MRI scan used to measure the hemodynamic response (change in blood flow) related to neural activity in the brain or spinal cord of humans or other animals. It is one of the most recently developed forms of neuroimaging. Since the early 1990s, fMRI has come to dominate the brain mapping field due to its relatively low invasiveness, absence of radiation exposure, and relatively wide availability.
- 1 Background
- 2 Neural correlates of BOLD
- 3 Technique
- 4 fMRI studies draw from many disciplines
- 5 Advantages and Disadvantages of fMRI
- 6 Scanning in practice
- 7 Related techniques
- 8 Approaches to fMRI data analysis
- 9 Cost of fMRI
- 10 Commercial use
- 11 See also
- 12 References
- 13 External links
Since the 1890s it has been known that changes in blood flow and blood oxygenation in the brain (collectively known as hemodynamics) are closely linked to neural activity. When neural cells are active they increase their consumption of energy from glucose and switch to less energetically effective, but more rapid anaerobic glycolysis. The local response to this energy utilization is to increase blood flow to regions of increased neural activity, which occurs after a delay of approximately 1–5 seconds. This hemodynamic response rises to a peak over 4–5 seconds, before falling back to baseline (and typically undershooting slightly). This leads to local changes in the relative concentration of oxyhemoglobin and deoxyhemoglobin and changes in local cerebral blood volume and in local cerebral blood flow.
Blood-oxygen-level dependence (BOLD) is the MRI contrast of blood deoxyhemoglobin, first discovered in 1990 by Seiji Ogawa at AT&T Bell labs. Ogawa and colleagues had recognized the potential importance of BOLD for functional brain imaging with MRI, but the first successful fMRI study was reported by John W. Belliveau and colleagues in 1991 using an intraveneously administered paramagnetic contrast agent (Gadolinium). Using a visual stimulus paradigm, localized increases in blood volume (32 +/- 10 percent, n = 7 subjects) were detected in the primary visual cortex. In 1992, three papers were published using endogenous BOLD contrast MRI. One was submitted by Peter Bandettini at the Medical College of Wisconsin on February 5, revised March 31, accepted March 31 and published in the June 1992 issue of Magnetic Resonance in Medicine (MRM). The second by Kenneth Kwong and colleagues also applied BOLD to image human brain activities with MRI and was submitted on March 26 and published in the June issue of PNAS in 1992. In the same year, Dr. Ogawa submitted their result on March 31 and published in July issue of PNAS. In the following year, Dr. Ogawa published the biophysics model of BOLD contrast in Biophysical Journal. Dr. Bandettini also published a further paper in 1993 demonstrating quantitative determination of functional activation maps.
As neurons do not have internal reserves for glucose and oxygen, more neuronal activity requires more glucose and oxygen to be delivered rapidly through the blood stream. Through a process called the hemodynamic response, blood releases glucose to neurons and astrocytes at a greater rate than in the area of inactive neurons. It results in a surplus of oxyhemoglobin in the veins of the area and distinguishable change of the local ratio of oxyhemoglobin to deoxyhemoglobin, the "marker" of BOLD for MRI.
Hemoglobin is diamagnetic when oxygenated (oxyhemoglobin) but paramagnetic when deoxygenated (deoxyhemoglobin). The magnetic resonance (MR) signal of blood is therefore slightly different depending on the level of oxygenation. Higher BOLD signal intensities arise from increases in the concentration of oxygenated hemoglobin since the blood magnetic susceptibility now more closely matches the tissue magnetic susceptibility. By collecting data in an MRI scanner with sequence parameters sensitive to changes in magnetic susceptibility one can assess changes in BOLD contrast. These changes can be either positive or negative depending upon the relative changes in both cerebral blood flow (CBF) and oxygen consumption. Increases in CBF that outstrip changes in oxygen consumption will lead to increased BOLD signal, conversely decreases in CBF that outstrip changes in oxygen consumption will cause decreased BOLD signal intensity. The signal difference is very small, but given many repetitions of a thought, action or experience, statistical methods can be used to determine the areas of the brain which reliably show more of this difference as a result, and therefore which areas of the brain are active during that thought, action or experience.
Almost all current fMRI research uses BOLD as the method for determining where activity occurs in the brain as the result of various experiences, but because the signals are relative and not individually quantitative, some question its rigor. Other methods which propose to measure neural activity more directly have been attempted (for example measurement of the Oxygen Extraction Fraction (OEF) in regions of the brain, which measures how much of the oxyhemoglobin in the blood has been converted to deoxyhemoglobin or direct detection of magnetic fields generated by neuronal currents), but because the electromagnetic fields created by an active or firing neuron are so weak, the signal-to-noise ratio is extremely low and statistical methods used to extract quantitative data have been largely unsuccessful as of yet.
Neural correlates of BOLD
The precise relationship between neural signals and BOLD is under active research. In general, changes in BOLD signal are well correlated with changes in blood flow. Numerous studies during the past several decades have identified a coupling between blood flow and metabolic rate;  that is, the blood supply is tightly regulated in space and time to provide the nutrients for brain metabolism. However, neuroscientists have been seeking a more direct relationship between the blood supply and the neural inputs/outputs that can be related to observable electrical activity and circuit models of brain function.
While current data indicate that local field potentials, an index of integrated electrical activity, form a marginally better correlation with blood flow than the spiking action potentials that are most directly associated with neural communication , no simple measure of electrical activity to date has provided an adequate correlation with metabolism and the blood supply across a wide dynamic range. Presumably, this reflects the complex nature of metabolic processes, which form a superset with regards to electrical activity. Some recent results have suggested that the increase in cerebral blood flow (CBF) following neural activity is not causally related to the metabolic demands of the brain region, but rather is driven by the presence of neurotransmitters, like glutamate, serotonin, nitric oxide, acetylcholine, dopamine and noradrenaline.
Some other recent results suggest that an initial small, negative dip before the main positive BOLD signal is more highly localized and also correlates with measured local decreases in tissue oxygen concentration (perhaps reflecting increased local metabolism during neuron activation). Use of this more localized negative BOLD signal has enabled imaging of human ocular dominance columns in primary visual cortex, with resolution of about 0.5 mm. One problem with this technique is that the early negative BOLD signal is small and can only be seen using larger scanners with magnetic fields of at least 3 Tesla. Further, the signal is much smaller than the normal BOLD signal, making extraction of the signal from noise more difficult. Also, this initial dip occurs within 1–2 seconds of stimulus initiation, which may not be captured when signals are recorded at long repetition (TR). If the TR is sufficiently low, increased speed of the cerebral blood flow response due to consumption of vasoactive drugs (such as caffeine) or natural differences in vascular responsiveness may further obscure observation of the initial dip.
The BOLD signal is composed of CBF contributions from larger arteries and veins, smaller arterioles and venules, and capillaries. Experimental results indicate that the BOLD signal can be weighted to the smaller vessels, and hence closer to the active neurons, by using larger magnetic fields. For example, whereas about 70% of the BOLD signal arises from larger vessels in a 1.5 tesla scanner, about 70% arises from smaller vessels in a 7 tesla scanner. Furthermore, the size of the BOLD signal increases roughly as the square of the magnetic field strength. Hence there has been a push for larger field scanners to both improve localization and increase the signal. A few 7 tesla commercial scanners have become operational, and experimental 8 and 9 tesla scanners are under development.
BOLD effects are measured using rapid volumetric acquisition of images with contrast weighed by T1 or T2*. Such images can be acquired with moderately good spatial and temporal resolution; images are usually taken every 1–4 seconds, and the voxels in the resulting image typically represent cubes of tissue about 2–4 millimeters on each side in humans. Recent technical advancements, such as the use of high magnetic fields and multichannel RF reception, have advanced spatial resolution to the millimeter scale. Although responses to stimuli presented as close together as one or two seconds can be distinguished from one another, using a method known as event-related fMRI, the full time course of a BOLD response to a briefly presented stimulus lasts about 15 seconds for the robust positive response.
fMRI studies draw from many disciplines
fMRI is a highly interdisciplinary research area and many studies draw on knowledge in several fields:
- Physics: Physical principles underlie fMRI signals and many studies require an understanding of these underlying principles.
- Psychology: Almost all fMRI studies are essentially cognitive psychological, cognitive psychophysiological, and/or psychophysical experiments in which the MRI scanner is used to obtain an extra set of measurements in addition to behavioral or electroencephalographic measurements.
- Neuroanatomy: The fMRI signals can be put into the context of previous knowledge only with an understanding of the neuroanatomy.
- Statistics: Correct application of statistics is essential to "tease out" observations and avoid false-positive results.
- Electrophysiology: Familiarity with neuronal behavior at the electrophysiological level can help investigators design a useful fMRI study.
Advantages and Disadvantages of fMRI
Like any technique, fMRI has advantages and disadvantages, and in order to be useful, the experiments that employ it must be carefully designed and conducted to maximize its strengths and minimize its weaknesses.
Advantages of fMRI
- It can noninvasively record brain signals without risks of radiation inherent in other scanning methods, such as CT or PET scans.
- It has high spatial resolution. 2–3 mm is typical but resolution can be as good as 1mm.
- It can record signal from all regions of the brain, unlike EEG/MEG which are biased towards the cortical surface.
- fMRI is widely used and standard data-analysis approaches have been developed which allow researchers to compare results across labs.
- fMRI produces compelling images of brain "activation".
Disadvantages of fMRI
- The images produced must be interpreted carefully, since correlation does not imply causality, and brain processes are complex and often non-localized.
- Statistical methods must be used carefully because they can produce false positives. One team of researchers studying reactions to pictures of human emotional expressions reported a few activated voxels in the brain of a dead salmon when no correction for multiple comparisons was applied, illustrating the need for rigorous statistical analyses.
- The BOLD signal is only an indirect measure of neural activity, and is therefore susceptible to influence by non-neural changes in the body. This also means that it is difficult to interpret positive and negative BOLD responses
- BOLD signals are most strongly associated with the input to a given area rather than with the output. It is therefore possible (although unlikely) that a BOLD signal could be present in a given area even if there is no single unit activity.
- fMRI has poor temporal resolution. The BOLD response peaks approximately 5 seconds after neuronal firing begins in an area. This means that it is hard to distinguish BOLD responses to different events which occur within a short time window. Careful experimental design can reduce this problem. Also, some research groups are attempting to combine fMRI signals that have relatively high spatial resolution with signals recorded with other techniques, electroencephalography (EEG) or magnetoencephalography (MEG), which have higher temporal resolution but worse spatial resolution.
- fMRI has often been used to show activation localized to specific regions, thus minimizing the distributed nature of processing in neural networks. Several recent multivariate statistical techniques work around this issue by characterizing interactions between "active" regions found via traditional univariate techniques.
- The BOLD response can be affected by a variety of factors, including: drugs/substances; age, brain pathology; local differences in neurovascular coupling; attention; amount of carbon dioxide in the blood; etc.
Scanning in practice
Subjects participating in a fMRI experiment are asked to lie still and are usually restrained with soft pads to prevent movement from disturbing measurements. Some labs also employ bite bars to reduce motion, although these are unpopular as they can be uncomfortable. Small head movements can be corrected for in post-processing of the data, but large transient motion cannot be corrected. Motion in excess of around 3 millimeters results in unusable data. Motion is an issue for all populations, but most especially problematic for subjects with certain medical conditions (e.g. Alzheimer's Disease or schizophrenia) or with young children. Participants can be habituated to the scanning environment and trained to remain still in an MRI simulator.
An fMRI experiment usually lasts between 15 minutes and an hour. Depending on the purpose of study, subjects may view movies, hear sounds, smell odors, perform cognitive tasks such as n-back, memorization or imagination, press a few buttons, or perform other tasks. Researchers are required to give detailed instructions and descriptions of the experiment plan to each subject, who must sign a consent form before the experiment.
Safety is an important issue in all experiments involving MRI. Potential subjects must ensure that they are able to enter the MRI environment. The MRI scanner is built around an extremely strong magnet (1.5 teslas or more), so potential subjects must be thoroughly examined for any ferromagnetic objects (e.g. watches, glasses, hair pins, pacemakers, bone plates and screws, etc.) before entering the scanning environment.
Aside from BOLD fMRI, there are other related ways to probe brain activity using magnetic resonance properties:
Diffusion based functional MRI
Neuronal activity produces some immediate physical changes in cell shape that can be detected because they affect the compartment shape and size for water diffusion. A much improved spatial and temporal resolution for fMRI data collection has now been achieved by using diffusion MRI methodology that can detect these changes in neurons. The abrupt onset of increased neuron cell size occurs before the metabolic response commences, is shorter in duration and does not extend significantly beyond the area of the actual cell population involved. This technique is a diffusion weighted technique (DWI). There is some evidence that similar changes in axonal volume in white matter may accompany activity and this has been observed using a DTI (diffusion tensor imaging) technique. The future importance of diffusion-based functional techniques relative to BOLD techniques is not yet clear.
An injected contrast agent such as an iron oxide that has been coated by a sugar or starch (to hide from the body's defense system), causes a local disturbance in the magnetic field that is measurable by the MRI scanner. The signals associated with these kinds of contrast agents are proportional to the cerebral blood volume. While this semi-invasive method presents a considerable disadvantage in terms of studying brain function in normal subjects, it enables far greater detection sensitivity than BOLD signal, which may increase the viability of fMRI in clinical populations. Other methods of investigating blood volume that do not require an injection are a subject of current research, although no alternative technique in theory can match the high sensitivity provided by injection of contrast agent.
Arterial spin labeling
Arterial Spin Labelling (ASL), also known as arterial spin tagging, is an MRI technique capable of measuring cerebral blood flow (CBF) in vivo. ASL is capable of providing cerebral perfusion maps, without requiring the administration of a contrast agent or the use of ionising radiation, as it uses magnetically-labelled endogenous blood water as a freely-diffusible tracer. It was first proposed in 1992  and has since benefited from a number of modifications aimed at improving its robustness. ASL can monitor changes in CBF with activation and fMRI studies can therefore be conducted using ASL instead of relying on the BOLD effect. ASL fMRI is less popular than BOLD, as it suffers from a lower signal to noise ratio, can be less sensitive to weak stimuli and its temporal resolution is poorer than in BOLD studies. On the plus side, it can provide quantitative measures of a single well-defined parameter, CBF, whose baseline value can also be determined in the same experiment. It has also been found to outperform BOLD in terms of stability to slow signal drifts and localization of the activation area. The ASL activation signal is believed to be dominated by changes in the capillary bed of the activated area of the cortex, wheareas the BOLD signal is likely to be dominated by changes in the oxygenation of nearby veins.
Magnetic resonance spectroscopic imaging
Magnetic resonance spectroscopic imaging (MRS) is another, NMR-based process for assessing function within the living brain. MRS takes advantage of the fact that protons (hydrogen atoms) residing in differing chemical environments depending upon the molecule they inhabit (H2O vs. protein, for example) possess slightly different resonant properties (chemical shift). For a given volume of brain (typically > 1 cubic cm), the distribution of these H resonances can be displayed as a spectrum.
The area under the peak for each resonance provides a quantitative measure of the relative abundance of that compound. The largest peak is composed of H2O. However, there are also discernible peaks for choline, creatine, N-acetylaspartate (NAA) and lactate. Fortuitously, NAA is mostly inactive within the neuron, serving as a precursor to glutamate and as storage for acetyl groups (to be used in fatty acid synthesis) — but its relative levels are a reasonable approximation of neuronal integrity and functional status. Brain diseases (schizophrenia, stroke, certain tumors, multiple sclerosis) can be characterized by the regional alteration in NAA levels when compared to healthy subjects. Creatine is used as a relative control value since its levels remain fairly constant, while choline and lactate levels have been used to evaluate brain tumors.
Diffusion tensor imaging
Diffusion tensor imaging (DTI) is a related use of MR to measure anatomical connectivity between areas. Although it is not strictly a functional imaging technique because it does not measure dynamic changes in brain function, the measures of inter-area connectivity it provides are complementary to images of cortical function provided by BOLD fMRI. White matter bundles carry functional information between brain regions. The diffusion of water molecules is hindered across the axes of these bundles, such that measurements of water diffusion can reveal information about the location of large white matter pathways. Illnesses that disrupt the normal organization or integrity of cerebral white matter (such as multiple sclerosis) have a quantitative impact on DTI measures.
fMRI and EEG
Functional MRI has high spatial resolution but relatively poor temporal resolution (of the order of several seconds). Electroencephalography (EEG) directly measures the brain's electrical activity, giving high temporal resolution (~milliseconds) but low spatial resolution. The two techniques are therefore complementary and may be used simultaneously to record brain activity.
Recording an EEG signal inside an MRI system is technically challenging. The MRI system introduces artifacts into the EEG recording by inducing currents in the EEG leads via Faraday induction. This can happen through several different mechanisms. An imaging sequence applies a series of short radiofrequency pulses which induce a signal in the EEG system. The pulses are short and relatively infrequent, so interference may be avoided by blanking (switching off) the EEG system during their transmission. Magnetic field gradients used during imaging also induce a signal, which is harder to remove as it is in a similar frequency range to the EEG signal. Current is also induced when EEG leads move inside the magnet bore (i.e. when the patient moves during the exam). Finally, pulsed blood flow in the patient in the static magnetic field also induces a signal (called a ballistocardiographic artifact), which is also within the frequency range of interest. The EEG system also affects the MRI scan. Metal in the EEG leads and electrodes can introduce susceptibility artifacts into MR images. Care must also be taken to limit currents induced in the EEG leads via the MRI RF system, which could heat the leads sufficiently to burn the subject.
Having simultaneously recorded EEG and fMRI data, the final hurdle is to co-register the two datasets, as each is reconstructed using a different algorithm, subject to different distortions.
Before the advent of fMRI functional neuroimaging was typically performed with positron emission tomography (PET) scanners or more rarely with SPECT scanners. Niels A. Lassen and his coworkers lead the earliest efforts of functional neuroimaging, using radioactive gases to construct images of the working brain.
These nuclear imaging techniques do not use the nuclear magnetic resonance property and employ entirely different scanners.
Approaches to fMRI data analysis
The ultimate goal of fMRI data analysis is to detect correlations between brain activation and the task the subject performs during the scan. The BOLD signature of activation is relatively weak, however, so other sources of noise in the acquired data must be carefully controlled. This means that a series of processing steps must be performed on the acquired images before the actual statistical search for task-related activation can begin.
For a typical fMRI scan, the 3D volume of the subject's head is imaged every one or two seconds, producing a few hundred to a few thousand complete images per scanning session. The nature of MRI is such that these images are acquired in Fourier transform space, so they must be transformed back to image space to be useful. Because of practical limitations of the scanner the Fourier samples are not acquired on a grid, and scanner imperfections like thermal drift and spike noise introduce additional distortions. Small motions on the part of the subject and the subject's pulse and respiration will also affect the images.
The most common situation is that the researcher uses a pulse sequence supplied by the scanner vendor, such as an echo-planar imaging (EPI) sequence that allows for relatively rapid acquisition of many images. Software in the scanner platform itself then performs the reconstruction of images from Fourier transform space. During this stage some information is lost (specifically the complex phase of the reconstructed signal). Some types of artifacts, for example spike noise, become more difficult to remove after reconstruction, but if the scanner is working well these artifacts are thought to be relatively unimportant. For pulse sequences not provided by the vendor, for example spiral EPI, reconstruction may have to be done by software running on a separate platform.
After reconstruction the output of the scanning session consists of a series of 3D images of the brain. The most common corrections performed on these images are motion correction and correction for physiological effects. Outlier correction and spatial and/or temporal filtering may also be performed. If the task performed by the subject is thought to produce bursts of activation which are short compared to the BOLD response time (on the order of 6 seconds), temporal filtering may be performed at this stage to attempt to deconvolve out the BOLD response and recover the temporal pattern of activation.
At this point the data provides a time series of samples for each voxel in the scanned volume. A variety of methods are used to correlate these voxel time series with the task in order to produce maps of task-dependent activation.
There are many software packages available for analysing fMRI data.
Reconstruction of MRI data needs to be tested, calibrated and confirmed. MRI can suffer from numerous artifacts that include, geometric distortions, Nyquist ghosting, and signal dropout. Medical Imaging Phantoms are used to provide a consistent geometrical source for calibration and testing purposes. Minute tumor changes can require recalibration by use of a phantom to quantify the change.
Cost of fMRI
The major cost of an fMRI experiment is the MR scanner itself. New 1.5 tesla scanners often cost between $1,000,000 USD and $1,500,000 USD. New 3.0 tesla scanners often cost between $2,000,000 and $2,300,000 USD. Construction of MRI suites can cost $500,000 USD.
MRI procedures themselves can vary considerably in cost but generally fall somewhere between $400 and $3,500, depending on the facility and which region of the body is being scanned. Extremity scans (feet, hands, etc.) tend to be lower in price while body scans (including the brain) tend to be higher.
Most fMRI scans are for research or clinical use. Commercial use is limited. However, a few companies have been set up that attempt to sell fMRI specific hardware or services for research or clinical use.
In using fMRI techniques for use in lie detection, activated areas of the brain are observed while the subject is making a statement. Depending on what regions are the most active, the technician might determine whether a subject is telling the truth or not. Since a specific combination of brain functions are needed in order to tell a lie, the simultaneous activation of these regions often indicates deception. This technology is in its early stages of development, and many of its proponents hope to replace older lie detection techniques.
In clinical trials, the usage of fMRI as a method of lie detection has appeared reliable, with studies from 2005 by Kozel et al. indicating a 90% to 93% success rate.
However, there is still a fair amount of controversy over whether these techniques are reliable enough to be used in a legal setting. Some studies indicate that while there is an overall positive correlation, there is a great deal of variation between findings and in some cases considerable difficulty in replicating the findings.
- Brain Mapping
- Issues in fMRI
- Brain function
- Event related fMRI
- SEEP fMRI
- List of neuroscience databases
- Functional neuroimaging
- Roy CS, Sherrington CS (January 1890). "On the Regulation of the Blood-supply of the Brain". Journal of Physiology 11 (1–2): 85–158.17. PMC 1514242. PMID 16991945.
- Raichle, M. E., Mintun, M. A. (2006). "Brain Work and Brain Imaging". The Annual Review of Neuroscience.
- Laureys, S., Boly, M., Tononi, G.: "Functional Neuroimaging" in "The Neurology of Consciousness: Cognitive Neuroscience and Neuropathology" (Laureys, S., Tononi, G. eds.), Academic Press-Elsevier, 2009, ISBN 978-0-12-374168-4, pp.31-42
- Ogawa, S., Lee, T.M., Nayak, A.S., and Glynn, P. (1990). "Oxygenation-sensitive contrast in magnetic resonance image of rodent brain at high magnetic fields". Magnetic Resonance in Medicine 14 (1): 68–78. doi:10.1002/mrm.1910140108. PMID 2161986.
- Belliveau JW, Kennedy DN, McKinstry RC, Buchbinder BR, Weisskoff RM, Cohen MS, Vevea JM, Brady TJ, and Rosen BR (1991). "Functional mapping of the human visual cortex by magnetic resonance imaging". Science 254 (5032): 716–719. doi:10.1126/science.1948051. PMID 1948051.
- KK Kwong, JW Belliveau, DA Chesler, IE Goldberg, RM Weisskoff, BP Poncelet, DN Kennedy, BE Hoppel, MS Cohen, R Turner, H Cheng, TJ Brady, and BR Rosen (1992). "Dynamic Magnetic Resonance Imaging of Human Brain Activity During Primary Sensory Stimulation". PNAS 89 (12): 5675–79. doi:10.1073/pnas.89.12.5675.
- OGAWA S, TANK DW, MENON R, ELLERMANN JM, KIM SG, MERKLE H, UGURBIL K (1992). "Intrinsic signal changes accompanying sensory stimulation: Functional brain mapping with magnetic resonance imaging". PNAS 89 (13): 5951–55. doi:10.1073/pnas.89.13.5951. PMC 402116. PMID 1631079.
- Ogawa S, Menon RS, Tank DW, Kim SG, Merkle H, Ellermann JM, Ugurbil K. (1993). "Functional brain mapping by blood oxygenation level-dependent contrast magnetic resonance imaging. A comparison of signal characteristics with a biophysical model". Biophysical Journal 64 (3): 803–12. Bibcode 1993BpJ....64..803O. doi:10.1016/S0006-3495(93)81441-3.
- Bandettini, P.A.; Jesmanowicz, A.; Wong, E.C.; Hyde, J.S. (1993). "Processing strategies for time-course data sets in functional MRI of the human brain". Magnetic Resonance in Medicine 30 (2): 161–173. doi:10.1002/mrm.1910300204. PMID 8366797.
- L Pauling and CD Coryell (1936). "The Magnetic Properties and Structure of Hemoglobin, Oxyhemoglobin and Carbonmonoxyhemoglobin". PNAS 22 (4): 210–6. doi:10.1073/pnas.22.4.210. PMC 1076743. PMID 16577697.
- Gusnard DA, Raichle ME (2001). "Searching for a baseline: Functional imaging and the resting human brain". Nature Reviews Neuroscience 2 (10): 685–694. doi:10.1038/35094500. PMID 11584306.
- Yablonskiy DA, Haacke EM (1994). "Theory of NMR signal behavior in magnetically inhomogeneous tissues: the static dephasing regime". Magnetic Resonance in Medicine 32 (6): 749–63. doi:10.1002/mrm.1910320610. PMID 7869897.
- Konn D, Gowland P, Bowtell R (2003). "MRI detection of weak magnetic fields due to an extended current dipole in a conducting sphere: a model for direct detection of neuronal currents in the brain". Magnetic Resonance in Medicine 50 (1): 40–49. doi:10.1002/mrm.10494. PMID 12815677.
- Magistretti PJ, Pellerin L. (July 1999). "Cellular mechanisms of brain energy metabolism and their relevance to functional brain imaging". Philosophical Transactions of the Royal Cociety of London 29 (354(1387)): 1155–63. PMC 1692634. PMID 10466143.
- Logothetis NK, Pauls J, Augath M, Trinath T, Oeltermann A. (July 2001). "Neurophysiological investigation of the basis of the fMRI signal". Nature 412 (6843): 150–7. doi:10.1038/35084005. PMID 11449264.
- Yang, G. and Iadecola, C. (1996). "Glutamate microinjections in cerebellar cortex reproduce cerebrovascular effects of parallel fiber stimulation". Am J Physiol Regul Integr Comp Physiol 271 (6): R1568–1575. PMID 8997354. http://ajpregu.physiology.org/cgi/content/abstract/271/6/R1568.
- Bonvento, Giles; Cholet, Nathalie and Seylaz, Jacques (2000). "Sustained attenuation of the cerebrovascular response to a 10 min whisker stimulation following neuronal nitric oxide synthase inhibition". Neurosci Res 37 (2): 163–166. PMID 10867178. http://www.sciencedirect.com/science/article/pii/S0168010200001097.
- Behzadi, Y. et al. (2006). "Caffeine reduces the initial dip in the visual bold response at 3 t". Neuroimage 32 (1): 9–15. doi:10.1016/j.neuroimage.2006.03.005. PMID 16635577.
- van der Zwaag W, Francis S, Head K, Peters A, Gowland P, Morris P, Bowtell R. (May 2009). "fMRI at 1.5, 3 and 7 T: Characterising BOLD signal changes". Neuroimage 47 (4): 1425–34. doi:10.1016/j.neuroimage.2009.05.015. PMID 19446641.
- Roemer PB, Edelstein WA, Hayes CE, Souza SP, Mueller OM. (November 1990). "The NMR phased array". Magn Reson Med. 16 (2): 192–225. doi:10.1002/mrm.1910160203. PMID 2266841.
- Pruessmann KP, Weiger M, Scheidegger MB, Boesiger P. (Nov 1999). "SENSE: Sensitivity encoding for fast MRI". Magn Reson Med. 42 (5): 952–62. doi:10.1002/(SICI)1522-2594(199911)42:5<952::AID-MRM16>3.0.CO;2-S. PMID 10542355.
- Griswold MA, Jakob PM, Heidemann RM, Nittka M, Jellus V, Wang J, Kiefer B, Haase A. (June 2002). "Generalized autocalibrating partially parallel acquisitions (GRAPPA)". Magn Reson Med. 47 (6): 1202–10. doi:10.1002/mrm.10171. PMID 12111967.
- Shmuel et al. Negative fMRI response correlates with decreases in neuronal activity in monkey visual area V1. Nature Neurosci. 9(4):569-577 (2006).
- Logothetis, N.K.; Pauls, J; Augath, M; Trinath, T; Oeltermann, A (2001). "Neurophysiological investigation of the basis of the fMRI signal". Nature 412 (6843): 150. doi:10.1038/35084005. PMID 11449264. http://www.ssc.uwo.ca/psychology/culhamlab/fmri/pdfs/Logothetis.pdf.
- Magalhaes, A. (2005). Functional magnetic resonance and spectroscopy in drug and substance abuse. Top Magnetic Resonance Imaging. 3, 247-251.
- Chen, C., Hou, B., Holodny, A. (2008). Effect of age and tumor grade on BOLD functional MR imaging in pre-operative assessment of patients with glioma. Radiology. 3, 971-978.
- Aguirre, G., Zarahn, E., and D’esposito, M. (1998). The variability of human BOLD hemodynamic responses. Neuroimage. 8 (4), 360-369.
- Corbetta, M., Miezin, F., Dobmeyer, S., Shulman, G., Petersen, S. (1990). Attentional modulation of neural processing of shape, color, and velocity in humans. Science. 4962, 1556-1559.
- Haller, S., Bartsch, A. (2009). Pitfalls in fMRI. European Radiology. 19, 2689-2706.
- Le Bihan D, et al (2006). "Direct and fast detection of neuronal activation in the human brain with diffusion MRI". PNAS 103 (21): 8263–8268. doi:10.1073/pnas.0600644103. PMC 1472461. PMID 16702549.
- Kohno S et al. (2009). "Water diffusion slowdown in the human visual cortex on visual stimulation precedes vascular responses". Journal of Cerebral Blood Flow & Metabolism 29 (6): 1197–1207. doi:10.1038/jcbfm.2009.45. PMID 19384332.
- Mandl RCW et al. (2008). Warrant, Eric. ed. "Functional diffusion tensor imaging: Measuring task-related fractional anisotropy changes in the human brain along white matter tracts". PLoS ONE 3 (11): 1–10. doi:10.1371/journal.pone.0003631.
- Detre, J.A., Leigh, J.S., Williams, D.S., and Koretsky, A.P. (1992). Perfusion imaging. Magn Reson Med 23, 37-45.
- Williams, D.S., Detre, J.A., Leigh, J.S., and Koretsky, A.P. (1992). Magnetic resonance imaging of perfusion using spin inversion of arterial water. Proc Natl Acad Sci U S A 89, 212-216.
- Buxton, R.B. Introduction to Functional Magnetic Resonance Imaging, Principles and Techniques, 2009, Cambridge University Press, Cambridge, UK
- No Lie MRI Inc
- Cephos Corporation
- Kozel FA et al. "A Replication Study of the Neural Correlates of Deception." Behavioral Neuroscience, Vol. 118, No. 4, 2004 http://teaching.ust.hk/~econ695/Neural%20Correlates%20of%20Different%20Types%20of%20Deception.pdf
- Kozel FA et al. "Detecting Deception Using Functional Magnetic Resonance Imaging." Biol Psychiatry. 2005 Sep 23 http://www.musc.edu/pr/fmri.pdf
- Spence, Sean A. "Playing Devil's advocate: The case against fMRI lie detection." Legal and Criminological Psychology, Volume 13, Number 1, February 2008 , pp. 11-25(15) http://www.ingentaconnect.com/content/bpsoc/lcp/2008/00000013/00000001/art00002
- Scott A. Huettel, Allen W. Song, Gregory McCarthy, Functional Magnetic Resonance Imaging, Sinauer Associates, 2004, ISBN 0-87893-288-7
- Richard B. Buxton, An Introduction to Functional Magnetic Resonance Imaging: Principles and Techniques, Cambridge Univ Press, 2002, ISBN 0-521-58113-3
- Roberto Cabeza and Alan Kingstone, Editors, Handbook of Functional Neuroimaging of Cognition, Second Edition, MIT Press, 2006, ISBN 0-262-03344-5
- www.mri-tutorial.com- MRI-TUTORIAL.COM | A free learning repository about neuroimaging
- RadiologyInfo- The radiology information resource for patients: Functional Magnetic Resonance Imaging of the Brain
- Introduction to FMRI from the Oxford Centre for Functional Magnetic Resonance Imaging of the Brain, Oxford University
- BrainMapping.ORG project Community web site for information Brain Mapping and methods
- fMRI Videos RadiologyTube.com - a collection of fMRI videos
- fMRI Methods Wiki - with advice on how to report on fMRI studies
- fMRI entry at Scholarpedia