fMRI Flashcards
What is another name for fMRI?
Blood Oxygenation Level Dependent Magnetic Resonance Imaging
== BOLD MRI
What is the sole purpose of fMRI?
We want to use MRI to produce activation maps in response to some kind of task
–> Taking the brain and finding out which regions are activated under a specific task condition
What is the mechanism of fMRI?
The fMRI BOLD signal is altered due to the increase of blood flow in response to brain activity
Stimulus = Neuronal activity = NVC = Haemodynamic response = Detection by MRI scanner = fMRI BOLD response
What did Angelo Mosso study?
Subject laid on a delicately balanced table which could tip downward either at the head or at the foot if the weight of either end were increased
The moment emotional or intellectual activity began in the subject, down went the balance at the head end, in consequence of the redistribution of blood in his system
How does blood flow change the fMRI signal?
The BOLD sequence isn’t measuring the blood flow, but is sensitive to the different magnetic properties of oxygenated and deoxygenated red blood cells
The increased magnetic field in the deoxygenated blood causes a magnetic field gradient resulting in adjacent water molecules precessing at different frequencies (larmour equation –> increase in magnetic field = protons speed up)
As water molecules are precessing at different frequencies in the gradient field, they dephase causing a loss of signal around deoxygenated blood compared to oxygenated blood
What is magnetic susceptibility?
Indicates whether a material is attracted into or repelled out of a magnetic field
Paramagnetic materials align with the applied field and are attracted to regions of greater magnetic field
Diamagnetic materials are anti-aligned and are pushed away, toward regions of lower magnetic fields
What is the magnetic susceptibility of oxygenated and deoxygenated blood?
Oxygenated blood = Diamagnetic = slightly reduces the magnetic field (it is oriented in the opposite direction to the main magnetic field)
Deoxygenated blood = Paramagnetic = adds to the magnetic field (Oriented in line with the main magnetic field)
Explain activation contrast
At rest, oxygenated blood travels through the arteries and exchanges in the capillary bed, deoxygenated blood leaves via the veins
The deoxygenated blood increases dephasing, reducing T2* signal
During activation, an excess of oxygenated blood flows into the activated region, swamping the deoxygenated blood
The excess of oxygenated blood reduces dephasing, increasing the T2* signal
What is the reason for the oversupply of oxygenated blood?
The reasons for the oversupply of oxygenated blood are still unclear, but may be related to the supply of other nutrients such as glucose
It could also be to provide a strong enough concentration gradient across the capillary wall so that the oxygenated blood can diffuse efficiently
An alternative hypothesis is that it is related to blood flow being controlled at the arteriole level - could be due to the control of the blood flow being higher up the vascular tree so to be safe of getting sufficient blood flow it increases it higher up the tree
What are we actually imaging in fMRI?
Not imaging individual red blood cells or even vessels
Imaging voxels are an average signal over mm3, typically 4 x 4 x 4 = 64 mm3 in fMRI
An fMRI voxel will contain many different vessels, but these still only make up around 3% of the voxel volume
What can be found in an individual voxel in fMRI?
Approx. =
6 million neurons
250km of axons
5 x 10 (to power 10) synapses
25km of dendrites
How do we image through k space?
We move through k-space by the negative phase encode and readout gradients
- Initially we have slice select at 90 degree pulse, selects a slice in the z direction
- The first line of k-space is acquired with a positive readout gradient during the first TR after applying the largest negative phase encode gradient
- During the second TR, the phase encode gradient is slightly less negative to read the second line of k-space
- We continue the process, incrementing the phase encoding gradient strength to read a different line of k-space during each TR
- Finally we apply the strongest positive phase encoding gradient to read the last line of k-space
What is a problem with acquiring an image in this way?
This process is far too slow for functional brain imaging
What method is used to image the BOLD signal which is faster?
Echo planar imaging
Explain echo planar imaging
- Negative y gradient, negative x gradient (moves southwest)
- Positive x gradient (moves along k space to the right)
- Apply a gradient blip in y axis (moves up k space tiny amount)
- Negative x gradient (moves back along k space to the left)
- Another gradient blip to move up k space
- Positive x gradient (move to right of k space)
- Continue until we have filled k space
What are the advantages of EPI?
Can be done in one TR
Means we can image the BOLD effect with a spin echo sequence
Single slice in under 100ms, so whole brain volumes in the order of a second or two
The whole of k space is acquired after a single 90 degree pulse - single shot imaging
What is EPI acquisition?
Echo planar imaging acquisition
What is a disadvantage of BOLD EPI acquisition?
Very low quality image
EPI is notorious for artefacts
Why does EPI create such a low quality image with artefacts?
Arises due to the imperfections in the rephasing-dephasing cycle of the rapidly switching bipolar frequency-encode gradient
We do it all after one excitation, so we don’t get a lot of signal at the two extremes, and those spaces at the edge of k space give a sharper image with more detail
What is the difference between sequential and interleaved slice acquisition?
Sequential slice acquisition acquires each adjacent slice consecutively, either bottom-to-top or top-to-bottom. Interleaved slice acquisition acquires every other slice, and then fills in the gaps on the second pass
Why do we normally interleave the slices?
Interleaving allows the use of a lower sampling bandwidth with a significant increase in signal-to-noise. The method also has the advantages of relative ease of implementation, no need for postprocessing to remove image distortion, and no need for shimming on a case-by-case basis
What causes signal drop out in BOLD MRI?
Huge susceptibility gradient between air and tissue which can cause signal drop out
How can we solve the problem of signal drop out in BOLD MRI?
Can reduce TE (echo time) and flip angle to manage signal drop out
What is a problem with reducing TE and flip angle?
It also reduced T2* weighting of the sequence and therefore reduces the BOLD signal that we measure
Why do we acquire dummy scans during BOLD imaging?
The first sequences have different brightness and contrast so we acquire dummy scans allowing the signal to reach steady state
What is the next step after acquiring the BOLD signal images?
We now have a BOLD EPI temporal imaging series and now we want to investigate the signal intensity time course voxel by voxel
How do we model the BOLD signal?
Model based
1. Linear regression
2. Nonlinear regression
Model free/data driven
1. PCA (principal component analysis)
2. ICA ( independent component analysis)
What is the purpose of using a model in BOLD analysis?
We want to fit a model to our data so that we can analyse which voxels show significant activation
Use a model to fit voxel signal time course
What is a problem with the model?
Model is a poor fit to the data - the data is clearly shifted in time with the response delayed compared to the stimulus
What is HRF?
Haemodynamic response function = the change in MR signal triggered by instantaneous neuronal activity
How do we account for the shift in the data?
Try and model the physiology - this is done using the haemodynamic response function which characterises the response to an impulse stimulus
It was derived experimentally and shows that the vascular response is delayed by a few seconds, peaking at around 5s before decaying
We can convolve the haemodynamic response function (HRF) with our experimental model to produce a more realistic estimate of the expected signal change
We can then fit the HRF convolved model to our data - it provides a much better fit to the voxel signal time course, but there appears to be some residual signal drift
How do we account for the residual drift in signal over time?
Add a linear ramp to the model
Next we need a value that tells us how well our model fits the data - how is this achieved?
Our data is a series of discrete signal measurements taken every few seconds
First thing is to estimate the amplitude - this is generated by the GLM fitting process
We then need to estimate the residual noise - calculated from the distance between the model function and acquired data at each point