Advanced MRI Techniques: Diffusion and Functional Imaging

Diffusion and Functional MRI Techniques

Introduction

Tom Barrack's lecture comprehensively covers advanced MRI techniques, specifically focusing on diffusion and functional MRI (fMRI) for neuroimaging. The objectives include understanding:

• The fundamental origin of the diffusion signal in MRI, elaborating on how molecular motion influences detectable signals.

• Diffusion Tensor Imaging (DTI) technique, detailing the mathematical and physical principles.

• Clinical applications of DTI, with specific examples such as stroke diagnosis, white matter disease characterization, and surgical planning.

• Diffusion tractography techniques, including algorithms and limitations in visualizing neural pathways.

• Functional Magnetic Resonance Imaging (fMRI) and its diverse applications in cognitive neuroscience and clinical diagnostics.

Diffusion Weighted Imaging

Diffusion

Diffusion in MRI refers to the random translational motion of water molecules, which is crucial for understanding tissue microstructure.

• Einstein Equation: This equation quantitatively links mean squared displacement to the diffusion coefficient (D) and diffusion time (T), providing a basis for measuring diffusion parameters. $Mean Squared Displacement = D \cdot T$

• Diffusion's time-dependent nature means that as observation time increases, molecules diffuse further, affecting MRI signal attenuation.

• The root mean square displacement of a molecule increases proportionally with time, illustrating continuous molecular dispersion.

Diffusion Coefficient

The diffusion coefficient is influenced by temperature and the specific molecule, both of which affect molecular mobility and signal characteristics.

• For water at 37 degrees Celsius (body temperature), the diffusion coefficient is approximately $3 \times 10^{-3} m^2/s$. This applies to free water without obstructions.

Types of Diffusion

• Isotropic Diffusion:

  • In free water, diffusion is uniform in all directions, such as in cerebrospinal fluid (CSF), where water molecules move without directional preference.

  • Molecules diffuse similarly in all directions over a given diffusion time if unobstructed, leading to predictable signal behavior.

  • Root mean square displacement is spherical, indicating equal diffusion probabilities in all orientations.

• Anisotropic Diffusion:

  • In environments like white matter, diffusion is greater along the cylinder than across it due to structures like myelin sheaths.

  • Root mean square displacement is ellipsoidal, reflecting directional diffusion biases.

Diffusion in Tissue

• Cerebrospinal Fluid (CSF): Contains free water, resulting in isotropic diffusion due to the unconstrained movement of water molecules.

• Grey Matter: Shows hindered diffusion due to complex tissue microstructure, affecting diffusion coefficients. Water diffusion occurs within cellular structures, with slight permeability variations between intracellular and extracellular spaces. The observed diffusion coefficient is about $0.8 \times 10^{-3} m^2/s$, lower than free water because of the hindering tissue.

• White Matter: Exhibits hindered and restricted diffusion because of myelin sheaths, causing preferential diffusion along axonal structures. This anisotropic diffusion results in an apparent diffusion coefficient around $0.8 \times 10^{-3} m^2/s$.

Diffusion MRI

Diffusion MRI quantifies water diffusion in tissue microstructure, measured in mm²/s, providing detailed insights into tissue integrity.

• Interaction of water with cellular structures helps in understanding tissue state and pathology.

• It infers the health of white matter microstructure, detecting damage and disease non-invasively without radiation.

• Spatial resolution is approximately 2 mm³ voxels, allowing detailed microstructural analysis.

Diffusion Weighted Imaging (DWI)

• In DWI, signal loss magnitude in a voxel indicates diffusion extent, providing a quantitative measure of water molecule displacement.

• The diffusion coefficient is calculated from this signal loss, enabling precise quantification of tissue properties.

• Diffusion is quantified in different spatial directions to model isotropy and anisotropy, revealing tissue organization.

• This leads to Diffusion Tensor Imaging (DTI) for detailed analysis.

Diffusion Weighted Sequences

DWI sequences use T2 weighted acquisition with added diffusion weighting using pulse gradients to enhance sensitivity to water diffusion.

The process:

1. T2 Weighted Acquisition: Standard pulse sequence to generate baseline anatomical images.

2. Diffusion Weighting: Pulse gradients manipulate the MR signal to be sensitive to diffusion.

3. Magnetization diffusion alters the refocusing of water molecules, which affects signal intensity.

4. The refocusing effect is diminished by water diffusion, leading to smaller echo signals where diffusion occurs.

Increasing the B value enhances diffusion weighting, resulting in greater signal attenuation in areas of high diffusion. Typical B values range from 0 to 1000 s/mm².

Diffusion Coefficient Calculation

$Signal(B) = Signal_0 \cdot e^{-B \cdot D}$

Where:

• Signal(B) is the signal intensity at a specific B value.

• $Signal_0$ is the initial signal intensity without diffusion weighting.

• B is the diffusion weighting factor, determining gradient strength and duration.

• D is the diffusion coefficient, quantifying molecular mobility.

This equation can be rearranged to solve for D, allowing for quantitative assessment of diffusion characteristics.

Water diffusion in tissue is hindered by macromolecules and restricted by cell membranes, affecting the diffusion signal.

A large diffusion coefficient causes high signal loss, evident in CSF, whereas tissue shows reduced diffusion.

The apparent diffusion coefficient (ADC) is calculated voxel-by-voxel, mapping spatial diffusion variations.

Clinical Diffusion Weighting

Diffusion weighting is applied in three orthogonal directions (x, y, and z) to capture comprehensive diffusion information.

The average diffusion coefficient is computed to provide an overall measure of tissue diffusion.

• Acute stroke example: High signal on DWI indicates restricted diffusion, coupled with low ADC values.

Stroke and Diffusion Changes

• Normal Tissue: Shows typical cellular structure with normal extracellular space.

• Acute Stroke: Cell swelling reduces extracellular space and increases intracellular space, leading to significant diffusion restriction and a reduced ADC. T2 signal remains relatively constant initially.

• Chronic Stroke: Cell damage and destruction expand extracellular space. Membrane damage and cell loss increase ADC and T2 values.

Anisotropic Diffusion

• In white matter, myelinated axons enhance diffusion along neurons but restrict it across, leading to anisotropic diffusion.

• Gradient direction across the neuronal structure results in a low ADC.

• Gradient along the neuronal structure results in a higher ADC.

• This anisotropy is used to build 3D diffusion models of the brain.

Diffusion Tensor Imaging (DTI)

DTI constructs 3D diffusion models by applying diffusion-weighted images in at least six non-collinear directions, along with an image without diffusion sensitization (B=0).

• This creates a 3D model for each voxel, depicting either anisotropy (ellipsoid) or isotropy (sphere), crucial for detailed tissue analysis.

• Mean Diffusivity (MD): Represents average diffusivity in a voxel, uniform across white matter, quantifying isotropic diffusion (about $0.8 \times 10^{-3} m^2/s$). Structural damage increases MD.

• Fractional Anisotropy (FA): Measures anisotropy at a brain point, indicating diffusion differences along and across neural structures (ranging from 0 to 1). FA indicates structural integrity, structural complexity and tissue damage.

Sensitivity to Tissue Damage

FA detects white matter changes with age and disease.

• FA decreases with aging, indicating structural decline.

• Decline is more severe in cerebrovascular small vessel disease and multiple sclerosis, indicating pronounced damage.

• Mean diffusivity increases with aging because of changes in the tissue microenvironment. During white matter development, FA increases, and overall mean diffusivity decreases.

Tissue Properties

• Aging Effects: FA increases from 10 to 18 years, then decreases, mirroring developmental and degenerative processes.

• Small Vessel Disease: Leads to perivascular structure breakdown, increasing mean diffusivity and decreasing FA because of axonal degeneration.

• Multiple Sclerosis: Detects remyelination and demyelination stages, with Gadolinium contrast enhancing lesion visualization. FA drops and mean diffusivity increases during demyelination; during remyelination, mean diffusivity decreases, and FA increases.

Diffusion Tractography

Tractography traces axonal structures using principal diffusion directions to visualize white matter connectivity.

The process:

1. Determine the principal diffusion directions across the long axis of diffusion.

2. Follow these directions step by step from a given voxel.

3. Join lines to infer white matter structural connectivity.

Applications of Tractography

• Extract white matter structures and bundles to assess tissue microstructure changes in diseases.

• Example: Extraction of frontal and occipital lobes or the superior longitudinal fasciculus (SLF) to study connectivity between Broca's area and Wernicke's area.

• Presurgical planning: Determine motor pathway distortion by tumors with low fractional anisotropy to guide surgical strategies and minimize damage.

Limitations of DTI and Tractography

• DTI provides only a single principal diffusion direction, limiting detection of multiple fiber orientations within a voxel.

• Fiber crossings are challenging to resolve without advanced techniques like constrained spherical deconvolution (CSD).

• Structural tissue changes are inferred from diffusion characteristics, not directly observed.

• It's a macroscopic measurement of microscopic diffusion effects.

• It's difficult to observe tissue pathological processes directly.

Clinical Applications

Primary clinical applications include:

• Three-directional DWI for detecting acute stroke.

• DTI for surgical planning and vascular disease observation.

Functional MRI (fMRI)

Introduction to fMRI

fMRI measures brain activity by detecting blood oxygenation and flow changes tied to neural activity.

• Active brain areas increase oxygen consumption, boosting local blood flow.

• fMRI generates activation maps showing brain regions involved in mental processes.

Advantages

• Non-invasive, ensuring patient safety.

• Good spatial resolution (3 mm³ voxels) for detailed mapping.

• Adequate temporal resolution (about 2 seconds per whole-brain image) for tracking dynamic brain activity.

Applications

• Imaging normal brain function, mapping cognitive and sensory processes.

• Providing insights into memory, language, pain, learning, and emotion.

• Assessing functional connectivity between regions.

Hemoglobin and BOLD

Measures oxygen delivery to neurons through hemoglobin in red blood cells.

• Neuronal activity increases oxygen demand, boosting local blood flow.

• Hemoglobin is diamagnetic when oxygenated and paramagnetic when deoxygenated, creating magnetic signal variations.

• These variations detect brain activity via blood oxygen level dependence (BOLD).

Hemodynamic Response

Following increased neural activity:

• Initial dip reflects oxygen consumption before blood flow increases.

• Blood flow peaks at about six seconds post-stimulation to meet demand.

• Blood oxygen levels then dictate temporal signal patterns.

fMRI Image Acquisition

Uses gradient echo images sensitive to T2* relaxation.

• T2* decreases with magnetic field uniformity loss, caused by deoxygenated blood's paramagnetic properties.

• Increased deoxyhemoglobin decreases the T2* signal, darkening the image.

• The contrast in T2* highlights differences between oxygenated and deoxygenated blood.

Physiological Process

1. Oxygenated hemoglobin in arterial blood has a long T2.

2. Oxygen exchange during metabolism produces deoxygenated hemoglobin, which is paramagnetic.

3. The brain vasodilates to increase blood flow and oxygen supply.

4. Deoxyhemoglobin briefly increases after arterial blood increases.

5. A larger arterial blood increase dilutes deoxyhemoglobin.

6. Gradient echo intensity increases by approximately 2% at 1.5 Tesla.

Visual Stimulus Example

• A subject alternates between a visual stimulus ('on' condition) and a dark screen ('off' condition) every 30 seconds to activate visual processing areas.

• Signal in the active brain area increases ('on') and decreases ('off').

• Voxels with a significant BOLD response are considered activated.

• The stimulus activates the visual cortex, with blue voxels indicating deactivation.

Applications of fMRI

• Task-based activation maps chart motor and sensory-motor cortices.

• In vivo determination of functional activation before surgery.

• Assessment of similarity between sensory-evoked potential and fMRI data.

• Pre- and post-surgical function assessment around a tumor.

Limitations of fMRI

• The hemodynamic response indirectly measures neural activity.

• Comparing responses between individuals cannot differentiate neural versus physiological differences.

• fMRI should map functional networks, not localize functions alone.

• Single brain regions participate in multiple tasks, requiring network analysis.

Hemoglobin - Revisited

With activation, the key is that oxygenated hemoglobin is diamagnetic, hence we observe a signal increase. The amount of deoxyhemoglobin decreases during activation.

Additional Insights

Simultaneous functional and diffusion image acquisition explores functional effects and structural connectivity (via diffusion tractography). This helps