Perfusion-Based fmri Thomas T. Liu Center for Functional MRI University of California San Diego May 7, 2006 Goal To provide a basic understanding of the theory and application of arterial spin labeling in functional MRI. 1
Topics 1. What is Cerebral Blood Flow? 2. Arterial Spin Labeling 3. Data Processing 4. Applications Cerebral blood flow (CBF) is a fundamental physiological quantity. Closely related to brain function. From C. Iadecola 2004 2
Cerebral Blood Flow (CBF) CBF = Perfusion = Rate of delivery of arterial blood to a capillary bed in tissue. Units: (ml of Blood) (100 grams of tissue)(minute) Typical value is 60 ml/(100g-min) or 60 ml/(100 ml-min) = 0.01 s -1, assuming average density of brain equals 1 gm/ml High CBF Low CBF Time 3
Bereczki et al 1992 Methods for Measuring CBF 1. Quantitative Autoradiography 2. Laser Doppler 3. Diffusible and Radioactive Tracer Based Methods 4. PET and SPECT 5. Contrast-Based MRI 6. Arterial Spin Labeling 4
Arterial spin labeling (ASL) 1: Tag by Magnetic Inversion Acquire image 2: Control Acquire image Control - Tag = ΔM CBF Pulsed ASL Continuous ASL Velocityselective ASL Imaging Region Tag Spatial inversion Mz(blood) control ΔM Flow-driven inversion Mz(blood) control Velocity-dependent saturation ΔM TI tag t t tag TI tt 5
PASL/VSASL TR ASL Pulse Sequences Tag Acquire Control Acquire CASL TI = Inversion Time Tag Acquire Control Acquire Labeling Time Post Labeling Delay Conventional Pulsed ASL (PASL) imaging slice presaturation slice tag control off-resonance IR pulse EPISTAR PICORE FAIR Courtesy of Wen-Ming Luh 6
Continuous ASL (CASL) Imaging Plane Inversion Planes Blood Magnetization B 0 Tag Conventional Control Amplitude Modulated Control Velocity Selective ASL (VSASL) 1 Control M z Image 0 { 0.1 Physiological Motion Tag 1 10 Velocity (cm/s) Wong et al 2006 7
CASL PASL PASL VSASL Credit: E. Wong Mz(blood) control TI tag ASL Signal Equation ΔM t t - = ΔM M= CBF A eff A eff is the effective area of the arterial bolus. It depends on both physiology and pulse sequence parameters. 8
Quantitative ASL M= CBF A eff Goal: Make A eff a well-controlled parameter that is robust to assumptions about physiological parameters. Major Sources of Error for ASL Transit Delays Bolus Width in PASL Relaxation Effects - different relaxation rates for blood and tissue, time of exchange. Intravascular signal -- blood destined to perfuse more distal slices Offset bias to due imperfect subtraction of static tissue -- slice profiles, magnetization transfer effects 9
Transit Delays CASL PASL ~ 3 cm t < 1000ms ~ 1 cm t < 700ms Controlling for Transit Delays in CASL Tagging Plane A B A B Voxels A and B have the same CBF, but voxel B will appear to have lower CBF if the measurement is made too early. Proper selection of post-labeling delay is key time 10
Amplifying Transit Delays Effects in CASL Tagging Plane Baseline Signal Activation Signal Baseline Signal Activation Signal Acquiring the signal at an earlier TI amplifies the difference between the activated state and the baseline state. time Arterial Bolus Width CASL Baseline Global flow increase Temporal Width of bolus determined by the pulse sequence PASL Baseline Global flow increase time Temporal Width of bolus determined by arterial velocity and size of tagging slab. Underestimates global flow changes. 11
Defining Bolus Width in PASL (QUIPSS II) Tag the spins TI 1 Saturate spins still in the slab Baseline Global flow increase Bolus temporal width = TI 1 QUIPSS II Modification Saturation Pulses Tag Acquire Control Acquire TI1 12
Controlling for Transit Delays in PASL Tagging Slab A B A B TI 1 TI 2 > transit delay + TI 1 Topics 1. What is Cerebral Blood Flow? 2. Arterial Spin Labeling 3. Data Processing 4. Applications 13
ASL Time Series Tag Control Tag Control Wait Tag by Magnetic Inversion Tag by Magnetic Inversion Image 1 Image 2 Image 3 Image 4-0.5 1 +0.5-0.5-1 +0.5 Perfusion Images CBF and BOLD Time Series 14
Static Tissue M[ n] BOLD Weighting b[ n] Perfusion " + " 1" # exp ("TI p /T 1 ) q[ n] " Measurements e[ n] 1"# ( 1+ ("1) n )exp("ti /T 1B ) y[ n] + " Modulate ("1) n+1 BOLD Estimate g[ n] b ˆ [n] g[ n] ˆ q [ n] Perfusion Estimate Demodulate 15
Static Tissue M[ n] BOLD Weighting b[ n] Perfusion " + " 1" # exp ("TI p /T 1 ) q[ n] " Measurements + e[ n] 1"# ( 1+ ("1) n )exp("ti /T 1B ) y[ n] Modulate Hernandez-Garcia et al 2005; Roller et al 2006, Abstract 540 16
Removing Physiological Noise Modeling a fmri signal using a General Linear Model (GLM) = + + + + + Data Design Matrix Nuisance Matrix Physiological Noise Matrix Additive Gaussian Noise y = Xh + Sb + Pc + n Removing Physiological Noise Perfusion Time Series Perfusion time series: Before Correction (cc = 0.15) Signal Intensity Perfusion time series: After Correction (cc = 0.71) Image Number (4 Hz) Pulse Ox (finger) Cardiac component estimated in brain 17
Removing Physiological Noise Young Subject (c) Older Subject (c) 0.500 0.425 (d) Corrected (d) Corrected 0.350 0.275 (a) Uncorrected (a) Uncorrected 0.200 (e) (e) (b) Corrected (b) Corrected Principal Components Based Correction mretroicor CompCor Behzadi et al Abstract 235 18
Tag Background Suppression Inversion Pulses Acquire Pre-saturation Pulse 1 0.8 0.6 0.4 QUIPSS II Saturation Pulse Background Suppression with Double Inversion GM WM CSF M z 0.2 0-0.2-0.4-0.6-0.8-1 0 200 400 600 800 1000 1200 1400 time (ms) Nulling of static tissue element Background Suppression Visual Study Results St. Lawrence et al 2005 Behzadi et al Abstract 3295 19
Single Shot Perfusion Labeling (SSPL) From JA Duyn et al, MRM 2001 Single Shot Perfusion Labeling (SSPL) From JA Duyn et al, MRM 2001 20
Topics 1. What is Cerebral Blood Flow? 2. Arterial Spin Labeling 3. Data Processing 4. Applications Applications of ASL 1. Quantitative Measures of both baseline and functional CBF 2. Multimodal measures of CBF, BOLD, (and CBV) to estimate functional changes in oxygen metabolism. 3. Experiments with long task periods 4. Mapping of functional activity. 21
Baseline CBF Modulates the BOLD Response Effect of Baseline CBF on BOLD Response 22
Effect of Acetazolamide on BOLD and CBF Percent BOLD Changes CBF Changes Acetazolamide 32 24 Drug Fee 28 Acetazolamide! S (Local Units) 18 12! CBF (Local Units) 24 20 16 Acetazolamide Administration Drug Free 12 6 8 Acetazolamide Administration 0 0 1 2 3 4 5 6 7 8 TRIAL 4 0 0 1 2 3 4 5 6 7 8 TRIAL Brown et al 2003 Effect of Hypercapnia on BOLD and CBF Stefanovic et al 2006 23
Applications of ASL 1. Quantitative Measures of both baseline and functional CBF 2. Multimodal measures of CBF, BOLD, (and CBV) to estimate functional changes in oxygen metabolism. 3. Experiments with long task periods 4. Mapping of functional activity. BOLD Dynamics Stimulus ATP ADP CMRO 2 Neural Activity CO ADP 2 NO, K + CBF ADP Lac CMRGlc CBV BOLD Signal Credit: Rick Buxton 24
Effect of Age on CBF and BOLD Responses in the Hippocampus Restom et al, Abstract 377 Estimation of CMRO 2 Changes with Combined CBF and BOLD measures Hypercapnia Contralateral Stefanovic et al 2004 Ipsilateral 25
Simultaneous measure of CBF, BOLD, and VASO (CBV) Yang et al 2004 VASO CBF BOLD Bayesian Analysis of dual -echo ASL data (Woolrich et al. Abstract 538) TE=9 msecs TE=30 msecs data model fit delivered blood static 26
Bayesian Analysis of ASL data (Woolrich et al. Abstract 538) CBF static magnetisation, M * R 2 Applications of ASL 1. Quantitative Measures of both baseline and functional CBF 2. Multimodal measures of CBF, BOLD, (and CBV) to estimate functional changes in oxygen metabolism. 3. Experiments with long task periods 4. Mapping of functional activity. 27
ASL with very low task frequencies - Wang et al., MRM 2003 28
ASL Imaging of Stress- Wang et al., PNAS 2005 ASL Imaging of Natural Vision; Rao et al, Abstract 167 CBF and BOLD activations during viewing of Road Runner Gee W-h-I-z-z-z-z-z-z Cartoon 29
Applications of ASL 1. Quantitative Measures of both baseline and functional CBF 2. Multimodal measures of CBF, BOLD, (and CBV) to estimate functional changes in oxygen metabolism. 3. Experiments with long task periods 4. Mapping of functional activity. ASL Mapping of Cortical Columns in Cat Visual Cortex Duong et al, PNAS, 2001. FAIR sequence, TI = 1500 ms, TR 3000 ms 30
Z. Liu and X. Hu, Abstract 891 Summary 1. Arterial spin labeling can provide non-invasive quantitative measures of baseline and functional CBF. 2. Since CBF is a well-defined physiological quantity, measures of CBF can be used to improve the interpretation of fmri studies, especially for clinical populations. 3. The combination of CBF measures with BOLD (and CBV) measures can be used to estimate changes in oxygen metabolism. 31
Acknowledgements Yashar Behzadi Rick Buxton Wen-Ming Luh Joanna Perthen Khaled Restom Eric Roller Mark Woolrich Eric Wong 32