Functional magnetic resonance imaging

Functional MRI or functional Magnetic Resonance Imaging (fMRI) is a type of specialized MRI scan. It measures the haemodynamic response 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 low invasiveness, lack of radiation exposure, and relatively wide availability.

Background

Since the 1890s (Roy and Sherrington, 1890) [Cite journal
author = Roy CS, Sherrington CS
title = On the Regulation of the Blood-supply of the Brain
journal = Journal of Physiology
year = 1890
month = January
volume = 11
issue = 1-2
pages = 85-158.17] 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 nerve cells are active they consume oxygen carried by hemoglobin in red blood cells from local capillaries. The local response to this oxygen utilization is an increase in blood flow to regions of increased neural activity, occurring 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 in addition to this change in local cerebral blood flow.

Blood-oxygen-level dependent or BOLD fMRI is a method of observing which areas of the brain are active at any given time. It was first found by Dr. Seiji Ogawa [Cite journal
author = Ogawa, S., Lee, T.M., Nayak, A.S., and Glynn, P.
year = 1990
title = Oxygenation-sensitive contrast in magnetic resonance image of rodent brain at high magnetic fields
journal = Magn Reson Med
volume = 14
pages = 68-78] in 1990 and following by Dr. Kenneth Kwong in 1992. [Cite journal
author = 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
title = Dynamic Magnetic Resonance Imaging of Human Brain Activity During Primary Sensory Stimulation
journal = PNAS
volume = 89
pages = 5675-79
year = 1992]
Neurons do not have internal reserves of energy in the form of glucose and oxygen, so their firing causes a need for more energy to be brought in quickly. Through a process called the hemodynamic response, blood releases oxygen to them at a greater rate than to inactive neurons, and the difference in magnetic susceptibility between oxyhemoglobin and deoxyhemoglobin, and thus oxygenated or deoxygenated blood, leads to magnetic signal variation which can be detected using an MRI scanner. Given many repetitions of a thought, action or experience, statistical methods can be used to determine the areas of the brain which reliably have more of this difference as a result, and therefore which areas of the brain are active during that thought, action or experience.

Almost all 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 directly have been attempted (for example, measurement of the Oxygen Extraction Fraction, or OEF, in regions of the brain, which measures how much of the oxyhemoglobin in the blood has been converted to deoxyhemoglobin [ [http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=pubmed&list_uids=7869897&dopt=Citation Theory of NMR signal behavior in magnetically inho... [Magn Reson Med. 1994 - PubMed Result ] ] ), 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.

Hemoglobin is diamagnetic when oxygenated but paramagnetic when deoxygenated. The magnetic resonance (MR) signal of blood is therefore slightly different depending on the level of oxygenation. These differential signals can be detected using an appropriate MR pulse sequence as blood-oxygen-level dependent (BOLD) contrast. 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 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.

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, especially glutamate.

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 [cite journal
author=Behzadi, Y. "et al"
title=Caffeine reduces the initial dip in the visual bold response at 3 t.
journal=Neuroimage
year=2006
volume=32
pages=9–15
doi=10.1016/j.neuroimage.2006.03.005] )or natural differences in vascular responsivnesses 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 4 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.

Technique

BOLD effects are measured using rapid volumetric acquisition of images with contrast weighed by T2 or T2* (see MRI). 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 advanced "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

To use fMRI effectively, an investigator must have a firm grasp of the relevant principles from all of these fields:

* Physics: Researchers should have a reasonable understanding of the physical principles underlying fMRI.
* 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 and electroencephalographic measurements. This allows for more detailed theory testing and inference on perceptual and cognitive processes, and allows to relate these to specific brain structures.
* 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 (of humans and other animals) without risks of radiation inherent in other scanning methods, such as CT scans.
* It can record on a spatial resolution in the region of 3-6 millimeters, but with relatively poor temporal resolution (on the order of seconds) compared with techniques such as EEG. However, this is mainly because of the phenomena being measured, not because of the technique. EEG measures electrical/neural activity while fMRI measures blood activity, which has a longer response. The MRI equipment used for fMRI can be used for high temporal resolution, if one measures different phenomena.

Disadvantages of fMRI

* The BOLD signal is only an indirect measure of neural activity, and is therefore susceptible to influence by non-neural changes in the body.

* 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. [cite journal
author=Logothetis, N.K.
year=2001
title=Neurophysiological investigation of the basis of the fMRI signal.
journal=Nature
volume=412
url=http://www.ssc.uwo.ca/psychology/culhamlab/fmri/pdfs/Logothetis.pdf
pages=150
doi=10.1038/35084005]

* Different brain areas have different hemodynamic responses, which would not be accurately reflected by the general linear model often used to filter fMRI time signals.

* For a non-invasive scan, fMRI has good spatial resolution. However, the temporal response of the blood supply, which is the basis of fMRI, is poor relative to the electrical signals that define neuronal communication. Therefore, some research groups are working around this issue by combining fMRI with data collection techniques such as electroencephalography (EEG) or magnetoencephalography (MEG), which have much higher temporal resolution but rather poor 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. Such techniques might prove useful in the future.

* fMRI is usually used to try to determine "where" task-related activity occurs in the brain. Some scientists prefer models which explain "how" psychological mechanisms function. The counter-argument to this criticism is that knowing "where" a cognitive function is located is vitally important. Neuropsychology, neurophysiology, and functional imaging each give us different windows of understanding into what each brain region does and how.

Scanning in practice

Subjects participating in a fMRI experiment are asked to lie still and are usually restrained with soft pads to prevent small motions from disturbing measurements. Some labs also employ bite bars to reduce motion, although these are unpopular as they can cause some discomfort to subjects. It is possible to correct for some amount of head movement with post-processing of the data, but large transient motion can render these attempts futile. Generally motion in excess of 3 millimeters will result in unusable data. The issue of motion is present for all populations, but most notably within populations that are not physically or emotionally equipped for even short MRI sessions (e.g., those with Alzheimer's Disease or schizophrenia, or young children). In these populations, various and negative reinforcement strategies can be employed in an attempt to attenuate motion artifacts, but in general the solution lies in designing a compatible paradigm with these populations.

An fMRI experiment usually lasts between 15 minutes and 2 hours. 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 a very important issue in all experiments involving MRI. Potential subjects must ensure that they are able to enter the MRI environment. Due to the nature of the MRI scanner, there is an extremely strong magnetic field surrounding the MRI scanner (at least 1.5 teslas, possibly stronger). 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.

Related techniques

Aside from BOLD fMRI, there are other related ways to probe brain activity using magnetic resonance properties:

Contrast MR

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

By magnetic labeling the proximal blood supply using "arterial spin labeling" ASL, the associated signal is proportional to the cerebral blood flow, or perfusion. This method provides more quantitative physiological information than BOLD signal, and has the same sensitivity for detecting task-induced changes in local brain function

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 (H₂O vs. protein, for example) possess slightly different resonant properties. 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 H₂O. 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 time resolution (of the order of several seconds). Electroecephalography (EEG) directly measures the brain's electrical activity, giving high temporal resolution (~milliseconds) but low spatial resolution. The two techniques are therefore complimentary and may be used simultaneously to record brain activity.

Recording an EEG signal inside an MRI system is technically challenging. The MRI system introduces artefacts 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 artefacts 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.

Nuclear neuroimaging

Before the rise 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 gasses 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 must 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.

Some fMRI neuroimaging software:
* AFNI [http://afni.nimh.nih.gov]
* BrainVoyager [http://www.brainvoyager.com]
* CamBA [http://sourceforge.net/projects/camba]
* Fiasco/FIAT [http://www.stat.cmu.edu/~fiasco]
* FreeSurfer [http://surfer.nmr.mgh.harvard.edu]
* mrVista [http://white.stanford.edu/newlm/index.php/Software]
* FSL [http://www.fmrib.ox.ac.uk/fsl]
* 3D Slicer [http://www.slicer.org]
* SPM [http://www.fil.ion.ucl.ac.uk/spm]
* [http://www.imagilys.com/autospm.html AutoSPM: Automated SPM for Surgical Planning]
* [http://www.bioimagesuite.org BioImage Suite]
* [http://www.nordicimaginglab.com/ nordicICE]
* [http://sccn.ucsd.edu/fmrlab/index.html FMRLAB]

Realtime fMRI

Realtime fMRI attempts to acquire and process brain activation data while a scan is in progress. A biofeedback loop can then be created by presenting a subject with their own brain activation patterns while they are being scanned. Using this technique, it has been investigated whether patients can use awareness of localized brain activation patterns to decrease symptoms of social anxiety disorder and chronic pain, with some reported success [Cite conference
author = R. A. Adcock el al.
title = Real time fMRI during the psychotherapy session: toward a methodology to augment therapeutic benefit
conference = 11th Annual Meeting of the Organization For Human Brain Mapping. June 12-16 2005, Toronto, Ontario, Canada.
year = 2005
booktitle = NeuroImage] . Other research groups have used the method in proof-of-concept applications to train subjects to control a game of pong using only their brain [Cite news
author = Mark Peplow
title = Mental ping-pong could aid paraplegics
publisher = Nature news
year = 2004
month = August
date = 2004-08-27
doi = 10.1038/news040823-18
url = http://www.nature.com/news/2004/040823/full/news040823-18.html] .

To date, only BOLD fMRI has been used in real-time applications, which delays the signal by approximately 2-5 seconds due to the physiological delay of the hemodynamic response. In the future, other methods of fMRI which do not rely on a secondary messenger like blood flow may reduce the delay and allow more immediate signal generation.

Cost of fMRI

Cost of fMRI scanners are about the same as MRI scanners because fMRI is only a specialized type of MRI scan ("functional" MRI). MRI equipment is expensive, though the cost is expected to exponentially decline soon in spurt of recent breakthroughs. 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.

Currently, in the US, there is increasing interest in reducing the costs associated with fMRI services and simultaneously improving the ability to effectively and efficiently provide fMRI examination services to larger numbers of researchers/other people with the same equipment.

Commercial use

Most fMRI scans are for research, and clinical use is limited. Research scans are usually done by the researchers themselves. However, a few companies have been set up that attempts to sell fMRI specific hardware or services for research or clinical use, e.g.,
* "Omneuron" [ [http://www.omneuron.com/] ] is a US-based company that is researching potential practical and clinical applications of real-time fMRI.
* "Applied fMRI Institute" [ [http://www.appliedfmri.org] ] is a San Diego, CA based company offering commercial use of their Siemens 3T TIM Trio.
* "Neurognostics" [http://www.neurognostics.com/] is a US-based company that offers a standardized fMRI system
* "Imagilys" [ [http://www.imagilys.com/] ] is a European company specialized in clinical and research fMRI.

At least two companies have been set up to use fMRI in lie detection ("No Lie MRI, Inc" [ [http://www.noliemri.com/] ] and "Cephos Corporation" [ [http://www.cephoscorp.com/] ] ).The utility of lie detection with fMRI is questioned.In episode 109 of the popular science show Mythbusters, the three members of the build team attempted to fool an fMRI test.Although two of them were unsuccessful, the third was able to successfully fool the machine.The signals are extrapolated from the fMRI machine onto a screen, displaying the active regions of the brain. Depending on what regions are the most active, the technician might determine whether a subject is telling the truth or not.This technology is in its early stages of development, and many of its proponents hope to replace older lie detection techniques.

ee also

* Brain Mapping
* Brain function
* Event related fMRI
* Spinal fMRI
* SEEP fMRI
* EEG-fMRI
* Real-time fMRI
* Functional neuroimaging
* The fMRI Data Centre
* Linear transform model

References

Textbooks

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-52158-113-3

Journal articles

cite journal
author=Weiller C "et al"
title=Role of functional imaging in neurological disorders
journal=Journal of Magnetic Resonance Imaging
year=2006
volume=23
issue=6
pages= 840–850
doi=10.1002/jmri.20591

cite journal
author=Lin, Lyons, and Berkowitz
title=Somatotopic Identification of Language-SMA in Language Processing via fMRI
journal=Journal of Scientific and Practical Computing
year=2007
volume=1
issue=2
pages= 3–8 [http://www.spclab.com/publisher/journals/Vol1No2/L1.pdf]

External links

* [http://radiologyinfo.org/en/info.cfm?pg=fmribrain RadiologyInfo] - The radiology information resource for patients: Functional Magnetic Resonance Imaging of the Brain
* [http://www.fmrib.ox.ac.uk/fmri_intro/ Introduction to FMRI] from the [http://www.fmrib.ox.ac.uk/ Oxford Centre for Functional Magnetic Resonance Imaging of the Brain] , Oxford University
* [http://www.brainmapping.org BrainMapping.ORG project] Community web site for information Brain Mapping and methods
* [http://www.mri-tutorial.com/tutorial_fmri.html A list of the best introductions to fMRI on the web]
* [http://www.fmrimethods.org/ fMRI Methods Wiki] - with advice on how to report on fMRI studies
* [http://www.scholarpedia.org/article/Functional_magnetic_resonance_imaging fMRI entry at Scholarpedia]