Sources of systematic error in calibrated BOLD based mapping of baseline oxygen extraction fraction

AbstractRecently a new class of calibrated blood oxygen level dependent (BOLD) functional magnetic resonance imaging (MRI) methods were introduced to quantitatively measure the baseline oxygen extraction fraction (OEF). These methods rely on two respiratory challenges and a mathematical model of the resultant changes in the BOLD functional MRI signal to estimate the OEF. However, this mathematical model does not include all of the effects that contribute to the BOLD signal, it relies on several physiological assumptions and it may be affected by intersubject physiological variability. The aim of this study was to investigate these sources of systematic error and their effect on estimating the OEF. This was achieved through simulation using a detailed model of the BOLD signal. Large ranges for intersubject variability in baseline physiological parameters such as haematocrit and cerebral blood volume were considered. Despite this the uncertainty in the relationship between the measured BOLD signals and the OEF was relatively low. Investigations of the physiological assumptions that underlie the mathematical model revealed that OEF measurements are likely to be overestimated if oxygen metabolism changes during hypercapnia or cerebral blood flow changes under hyperoxia. Hypoxic hypoxia was predicted to result in an underestimation of the OEF, whilst anaemic hypoxia was found to have only a minimal effect

A novel method of combining blood oxygenation and blood flow sensitive magnetic resonance imaging techniques to measure the cerebral blood flow and oxygen metabolism responses to an unknown neural stimulus.

Simultaneous implementation of magnetic resonance imaging methods for Arterial Spin Labeling (ASL) and Blood Oxygenation Level Dependent (BOLD) imaging makes it possible to quantitatively measure the changes in cerebral blood flow (CBF) and cerebral oxygen metabolism (CMRO(2)) that occur in response to neural stimuli. To date, however, the range of neural stimuli amenable to quantitative analysis is limited to those that may be presented in a simple block or event related design such that measurements may be repeated and averaged to improve precision. Here we examined the feasibility of using the relationship between cerebral blood flow and the BOLD signal to improve dynamic estimates of blood flow fluctuations as well as to estimate metabolic-hemodynamic coupling under conditions where a stimulus pattern is unknown. We found that by combining the information contained in simultaneously acquired BOLD and ASL signals through a method we term BOLD Constrained Perfusion (BCP) estimation, we could significantly improve the precision of our estimates of the hemodynamic response to a visual stimulus and, under the conditions of a calibrated BOLD experiment, accurately determine the ratio of the oxygen metabolic response to the hemodynamic response. Importantly we were able to accomplish this without utilizing a priori knowledge of the temporal nature of the neural stimulus, suggesting that BOLD Constrained Perfusion estimation may make it feasible to quantitatively study the cerebral metabolic and hemodynamic responses to more natural stimuli that cannot be easily repeated or averaged

Multiparametric measurement of cerebral physiology using calibrated fMRI

The ultimate goal of calibrated fMRI is the quantitative imaging of oxygen metabolism (CMRO2), and this has been the focus of numerous methods and approaches. However, one underappreciated aspect of this quest is that in the drive to measure CMRO2, many other physiological parameters of interest are often acquired along the way. This can significantly increase the value of the dataset, providing greater information that is clinically relevant, or detail that can disambiguate the cause of signal variations. This can also be somewhat of a double-edged sword: calibrated fMRI experiments combine multiple parameters into a physiological model that requires multiple steps, thereby providing more opportunity for error propagation and increasing the noise and error of the final derived values. As with all measurements, there is a trade-off between imaging time, spatial resolution, coverage, and accuracy. In this review, we provide a brief overview of the benefits and pitfalls of extracting multiparametric measurements of cerebral physiology through calibrated fMRI experiments

Magnetic resonance imaging of resting cerebral oxygen metabolism : applications in Alzheimer’s disease

The BOLD contrast employed in functional MRI studies is an ambiguous signal composed of changes in blood flow, blood volume and oxidative metabolism. In situations where the vasculature and metabolism may have been affected, such as in aging and in certain diseases, the dissociation of the more physiologically-specific components from the BOLD signal becomes crucial. The latest generation of calibrated functional MRI methods allows the estimation of both resting blood flow and absolute oxygen metabolism. The work presented here is based on one such proof-of-concept approach, dubbed QUO2, whereby taking into account, within a generalized model, both arbitrary changes in blood flow and blood O2 content during a combination of hypercapnia and hyperoxia breathing manipulations, yields voxel-wise estimates of resting oxygen extraction fraction and oxidative metabolism. In the first part of this thesis, the QUO2 acquisition protocol and data analysis were revisited in order to enhance the temporal stability of individual blood flow and BOLD responses, consequently improving reliability of the model-derived estimates. Thereafter, an assessment of the within and between-subject variability of the optimized QUO2 measurements was performed on a group of healthy volunteers. In parallel, an analysis was performed of the sensitivity of the model to different sources of random and systematic errors, respectively due to errors in measurements and choice of assumed parameters values. Moreover, the various impacts of the oxygen concentration administered during the hyperoxia manipulation were evaluated through a simulation and experimentally, indicating that a mild hyperoxia was beneficial. Finally, the influence of Alzheimer’s disease in vascular and metabolic changes was explored for the first time by applying the QUO2 approach in a cohort of probable Alzheimer’s disease patients and age-matched control group. Voxel-wise and region-wise differences in resting blood flow, oxygen extraction fraction, oxidative metabolism, transverse relaxation rate constant R2* and R2* changes during hypercapnia were identified. A series of limitations along with recommended solutions was given with regards to the delayed transit time, the susceptibility artifacts and the challenge of performing a hypercapnia manipulation in cohorts of elderly and Alzheimer’s patients.Le contraste BOLD employé dans les études d’imagerie par résonance magnétique fonctionnelle (IRMf) provient d’une combinaison ambigüe de changements du flux sanguin cérébral, du volume sanguin ainsi que du métabolisme oxydatif. Dans un contexte où les fonctions vasculaires ou métaboliques du cerveau ont pu être affectées, tel qu’avec l’âge ou certaines maladies, il est crucial d’effectuer une décomposition du signal BOLD en composantes physiologiquement plus spécifiques. La dernière génération de méthodes d’IRMf calibrée permet d’estimer à la fois le flux sanguin cérébral et le métabolisme oxydatif au repos. Le présent travail est basé sur une telle technique, appelée QUantitative O2 (QUO2), qui, via un model généralisé, prend en considération les changements du flux sanguin ainsi que ceux en concentrations sanguine d’O2 durant des périodes d’hypercapnie et d’hyperoxie, afin d’estimer, à chaque voxel, la fraction d’extraction d’oxygène et le métabolisme oxydatif au repos. Dans la première partie de cette thèse, le protocole d’acquisition ainsi que la stratégie d’analyse de l’approche QUO2 ont été revus afin d’améliorer la stabilité temporelle des réponses BOLD et du flux sanguin, conséquemment, afin d’accroître la fiabilité des paramètres estimés. Par la suite, une évaluation de la variabilité intra- et inter-sujet des différentes mesures QUO2 a été effectuée auprès d’un groupe de participants sains. En parallèle, une analyse de la sensibilité du model à différentes sources d’erreurs aléatoires (issues des mesures acquises) et systématiques (dues aux assomptions du model) a été réalisée. De plus, les impacts du niveau d’oxygène administré durant les périodes d’hyperoxie ont été évalués via une simulation puis expérimentalement, indiquant qu’une hyperoxie moyenne était bénéfique. Finalement, l’influence de la maladie d’Alzheimer sur les changements vasculaires et métaboliques a été explorée pour la première fois en appliquant le protocole QUO2 à une cohorte de patients Alzheimer et à un groupe témoin du même âge. Des différences en terme de flux sanguin, fraction d’oxygène extraite, métabolisme oxydatif, et taux de relaxation transverse R2* au repos comme en réponse à l’hypercapnie, ont été identifiées au niveau du voxel, ainsi qu’au niveau de régions cérébrales vulnérables à la maladie d’Alzheimer. Une liste de limitations accompagnées de recommandations a été dressée en ce qui a trait au temps de transit différé, aux artéfacts de susceptibilité magnétique, de même qu’au défi que représente l’hypercapnie chez les personnes âgées ou atteintes de la maladie d’Alzheimer

Field strength dependence of grey matter R2* on venous oxygenation

The relationship between venous blood oxygenation and change in transverse relaxation rate (ΔR2 *) plays a key role in calibrated BOLD fMRI. This relationship, defined by the parameter β, has previously been determined using theoretical simulations and experimental measures. However, these earlier studies have been confounded by the change in venous cerebral blood volume (CBV) in response to functional tasks. This study used a double-echo gradient echo EPI scheme in conjunction with a graded isocapnic hyperoxic sequence to assess quantitatively the relationship between the fractional venous blood oxygenation (1-Yv) and transverse relaxation rate of grey matter (ΔR2 * GM), without inducing a change in vCBV. The results demonstrate that the relationship between ΔR2 * and fractional venous oxygenation at all magnet field strengths studied was adequately described by a linear relationship. The gradient of this relationship did not increase monotonically with field strength, which may be attributed to the relative contributions of intravascular and extravascular signals which will vary with both field strength and blood oxygenation

Dynamic Assessment of Cerebral Metabolic Rate of Oxygen (cmro2) With Magnetic Resonance Imaging

The brain is almost entirely dependent on oxidative metabolism to meet its energy requirements. As such, the cerebral metabolic rate of oxygen (CMRO2) is a direct measure of brain energy use. CMRO2 provides insight into brain functional architecture and has demonstrated potential as a clinical tool for assessing many common neurological disorders. Recent developments in magnetic resonance imaging (MRI)-based CMRO2 quantification have shown promise in spatially resolving CMRO2 in clinically feasible scan times. However, brain energy requirements are both spatially heterogeneous and temporally dynamic, responding to rapid changes in oxygen supply and demand in response to physiologic stimuli and neuronal activation. Methods for dynamic quantification of CMRO2 are lacking, and this dissertation aims to address this gap. Given the fundamental tradeoff between spatial and temporal resolution in MRI, we focus initially on the latter. Central to each proposed method is a model-based approach for deriving venous oxygen saturation (Yv) – the critical parameter for CMRO2 quantification – from MRI signal phase using susceptometry-based oximetry (SBO). First, a three-second-temporal-resolution technique for whole-brain quantification of Yv and CMRO2 is presented. This OxFlow method is applied to measure a small but highly significant increase in CMRO2 in response to volitional apnea. Next, OxFlow is combined with a competing approach for Yv quantification based on blood T2 relaxometry (TRUST). The resulting interleaved-TRUST (iTRUST) pulse sequence greatly improves T2-based CMRO2 quantification, while allowing direct, simultaneous comparison of SBO- and T2-based Yv. iTRUST is applied to assess the CMRO2 response to hypercapnia – a topic of great interest in functional neuroimaging – demonstrating significant biases between SBO- and T2-derived Yv and CMRO2. To address the need for dynamic and spatially resolved CMRO2 quantification, we explore blood-oxygen-level-dependent (BOLD) calibration, introducing a new calibration model and hybrid pulse sequence combining OxFlow with standard BOLD/CBF measurement. Preliminary results suggest Ox-BOLD provides improved calibration “M-maps” for converting BOLD signal to CMRO2. Finally, OxFlow is applied clinically to patients with obstructive sleep apnea (OSA). A small clinical pilot study demonstrates OSA-associated reductions in CMRO2 at baseline and in response to apnea, highlighting the potential utility of dynamic CMRO2 quantification in assessing neuropathology

Direct Estimation of Evoked Hemoglobin Changes by Multimodality Fusion Imaging

In the last two decades, both diffuse optical tomography (DOT) and blood oxygen level dependent (BOLD)-based functional magnetic resonance imaging (fMRI) methods have been developed as noninvasive tools for imaging evoked cerebral hemodynamic changes in studies of brain activity. Although these two technologies measure functional contrast from similar physiological sources, i.e., changes in hemoglobin levels, these two modalities are based on distinct physical and biophysical principles leading to both limitations and strengths to each method. In this work, we describe a unified linear model to combine the complimentary spatial, temporal, and spectroscopic resolutions of concurrently measured optical tomography and fMRI signals. Using numerical simulations, we demonstrate that concurrent optical and BOLD measurements can be used to create cross-calibrated estimates of absolute micromolar deoxyhemoglobin changes. We apply this new analysis tool to experimental data acquired simultaneously with both DOT and BOLD imaging during a motor task, demonstrate the ability to more robustly estimate hemoglobin changes in comparison to DOT alone, and show how this approach can provide cross-calibrated estimates of hemoglobin changes. Using this multimodal method, we estimate the calibration of the 3tesla BOLD signal to be −0.55%±0.40% signal change per micromolar change of deoxyhemoglobin

Applicability of Quantitative Functional MRI Techniques for Studies of Brain Function at Ultra-High Magnetic Field

This thesis describes the development, implementation and application of various quantitative functional magnetic resonance imaging (fMRI) approaches at ultra-high magnetic field including the assessment with regards to applicability and reproducibility. Functional MRI (fMRI) commonly uses the blood oxygenation level dependent (BOLD) contrast to detect functionally induced changes in the oxy-deoxyhaemoglobin composition of blood which reflect cerebral neural activity. As these blood oxygenation changes do not only occur at the activation site but also downstream in the draining veins, the spatial specificity of the BOLD signal is limited. Therefore, the focus has moved towards more quantitative fMRI approaches such as arterial spin labelling (ASL), vascular space occupancy (VASO) or calibrated fMRI which measure quantifiable physiologically and physically relevant parameters such as cerebral blood flow (CBF), cerebral blood volume (CBV) or cerebral metabolic rate of oxygen (CMRO2), respectively. In this thesis a novel MRI technique was introduced which allowed the simultaneous acquisition of multiple physiological parameters in order to beneficially utilise their spatial and temporal characteristics. The advantages of ultra-high magnetic field were utilised to achieve higher signal-to-noise and contrast-to-noise ratios compared to lower field strengths. This technique was successfully used to study the spatial and temporal characteristics of CBV, CBF and BOLD in the visual cortex. This technique is the first one that allows simultaneous acquisition of CBV, CBF and BOLD weighted fMRI signals in the human brain at 7 Tesla. Additionally, this thesis presented a calibrated fMRI technique which allowed the quantitative estimation of changes in cerebral oxygen metabolism at ultra-high field. CMRO2 reflects the amount of thermodynamic work due to neural activity and is therefore a significant physical measure in neuroscience. The calibrated fMRI approach presented in this thesis was optimised for the use at ultra-high field by adjusting the MRI parameters as well as implementing a specifically designed radio-frequency (RF) pulse. A biophysical model was used to calibrate the fMRI data based on the simultaneous acquisition of BOLD and CBF weighted MRI signals during a gas-breathing challenge. The reproducibility was assessed across multiple brain regions and compared to that of various physiologically relevant parameters. The results indicate that the degree of intra-subject variation for calibrated fMRI is lower than for the classic BOLD contrast or ASL. Consequently, calibrated fMRI is a viable alternative to classic fMRI contrasts with regards to spatial specificity as well as functional reproducibility. This calibrated fMRI approach was also compared to a novel direct calibration technique which relies on complete venous oxygenation saturation during the calibration scan via a gas-breathing challenge. This thesis introduced several reliable quantitative fMRI approaches at 7 Tesla and the results presented are a step forward to the wider application of quantitative fMRI.:1 Introduction 3 2 Background to Functional Magnetic Resonance Imaging 7 2.1 Magnetic Resonance 7 2.1.1 Quantum Mechanics 7 2.1.2 The Classical Point of View 10 2.1.3 Radio Frequency Pulses 12 2.1.4 Relaxation Effects 13 2.1.5 The Bloch Equations 15 2.2 Magnetic Resonance Imaging 16 2.2.1 Data Acquisition 16 2.2.2 Image Formation 17 2.2.2.1 Slice Selection 17 2.2.2.2 Frequency Encoding 18 2.2.2.3 Phase Encoding 19 2.2.2.4 Mathematics of Image Formation 20 2.2.2.5 Signal Formation 22 2.3 Advanced Imaging Methods 24 2.3.1 Echo-Planar Imaging (EPI) 24 2.3.2 Partial Fourier Acquisition 25 2.3.3 Generalised Autocalibrating Partially Parallel Acquisition (GRAPPA) 25 2.3.4 Inversion Recovery (IR) 26 2.3.5 Adiabatic Inversion 26 2.3.5.1 Hyperbolic Secant (HS) RF pulses 28 2.3.5.2 Time Resampled Frequency Offset Corrected Inversion (tr-FOCI) RF Pulses 28 2.4 Physiological Background 29 2.4.1 Neuronal Activity 30 2.4.2 Energy Metabolism 31 2.4.3 Physiological Changes During Brain Activation 32 2.4.4 The BOLD Contrast 34 2.4.5 Disadvantages of the BOLD Contrast 35 2.5 Arterial Spin Labelling (ASL) 35 2.5.1 Pulsed Arterial Spin Labelling 37 2.5.2 Arterial Spin Labelling at Ultra-High Field 41 2.6 Vascular Space Occupancy (VASO) 42 2.6.1 VASO at Ultra-High Field 44 2.6.2 Slice-Saturation Slab-Inversion (SS-SI) VASO 45 2.7 Calibrated Functional Magnetic Resonance Imaging 47 2.7.1 The Davis Model 47 2.7.2 The Chiarelli Model 50 2.7.3 The Generalised Calibration Model (GCM) 52 3 Materials and Methods 53 3.1 Scanner Setup 53 3.2 Gas Delivery and Physiological Monitoring System 53 3.3 MRI Sequence Developments 55 3.3.1 Tr-FOCI Adiabatic Inversion 55 3.3.2 Optimisation of the PASL FAIR QUIPSSII Sequence Parameters 60 3.3.3 Multi-TE Multi-TI EPI 64 4 Experiment I: Comparison of Direct and Modelled fMRI Calibration 68 4.1 Background Information 68 4.2 Methods 69 4.2.1 Experimental Design 69 4.2.2 Visuo-Motor Task 70 4.2.3 Gas Manipulations 71 4.2.4 Scanning Parameters 71 4.2.5 Data Analysis 72 4.2.6 M-value Modelling 72 4.2.7 Direct M-Value Estimation 73 4.3 Results 74 4.4 Discussion 79 4.4.1 M-value Estimation 79 4.4.2 BOLD Time Courses 82 4.4.3 M-Maps and Single Subject Analysis 82 4.4.4 Effects on CMRO2 Estimation 83 4.4.5 Technical Limitations and Implications for Calibrated fMRI 84 4.5 Conclusion 89 5 Experiment II: Reproducibility of BOLD, ASL and Calibrated fMRI 90 5.1 Background Information 90 5.2 Methods 91 5.2.1 Experimental Design 91 5.2.2 Data Analysis 91 5.2.3 Reproducibility 93 5.2.4 Learning and Habituation Effects 95 5.3 Results 95 5.4 Discussion 101 5.4.1 Breathing Manipulations 102 5.4.2 Functional Reproducibility 107 5.4.3 Habituation Effects on Reproducibility 109 5.4.4 Technical Considerations for Calibrated fMRI 110 5.5 Conclusion 112 6 Experiment III: Simultaneous Acquisition of BOLD, ASL and VASO Signals 113 6.1 Background Information 113 6.2 Methods 114 6.2.1 SS-SI VASO Signal Acquisition 114 6.2.2 ASL and BOLD Signal Acquisition 114 6.2.3 Experimental Design 114 6.2.4 Data Analysis 115 6.3 Results 115 6.4 Discussion 116 6.5 Conclusion 120 7 Conclusion and Outlook 12

Elimination of visually evoked BOLD responses during carbogen inhalation: Implications for calibrated MRI

Breathing a mixture of 10% CO2 with 90% O2 (referred to here as carbogen-10) increases blood flow due to the vasodilatory effect of CO2, and raises blood O2 saturation due to the enriched oxygen level. These effects both tend to reduce the level of deoxygenated hemoglobin in brain tissues, thereby reducing the potential for further increases in BOLD contrast. In the present study, blocks of intense visual stimulation (60 s) were presented amid longer blocks (180 s) during which subjects breathed various fractional concentrations (0–100%) of carbogen-10 diluted with medical air. When breathing undiluted carbogen-10, the BOLD response to visual stimulation was reduced below the level of noise against the background of the carbogen-10 response. At these concentrations, the total (visual+carbogen) BOLD response amplitude (7.5±1.0%, n=6) converged toward that seen with carbogen alone (7.5 ± 1.0%, n = 6). In spite of the almost complete elimination of the visual BOLD response, pseudo-continuous arterial spin-labeling on a separate cohort indicated a largely preserved perfusion response (89±34%, n=5) to the visual stimulus during inhalation of carbogen-10. The previously discussed observations suggest that venous saturation can be driven to very high levels during carbogen inhalation, a finding which has significant implications for calibrated MRI techniques. The latter methods involve estimation of the relative change in venous O2 saturation by expressing activation-induced BOLD signal increases as a fraction of the maximal BOLD signal M that would be observed as venous saturation approaches 100%. While the value of M has generally been extrapolated from much smaller BOLD responses induced using hypercapnia or hyperoxia, our results suggest that these effects could be combined through carbogen inhalation to obtain estimates of M based on larger BOLD increases. Using a hybrid BOLD calibration model taking into account changes in both blood flow and arterial oxygenation, we estimated that inhalation of carbogen-10 led to an average venous saturation of 91%, allowing us to compute an estimated M value of 9.5%

Global intravascular and local hyperoxia contrast phase-based blood oxygenation measurements

AbstractThe measurement of venous cerebral blood oxygenation (Yv) has potential applications in the study of patient groups where oxygen extraction and/or metabolism are compromised. It is also useful for fMRI studies to assess the stimulus-induced changes in Yv, particularly since basal Yv partially accounts for inter-subject variation in the haemodynamic response to a stimulus. A range of MRI-based methods of measuring Yv have been developed recently. Here, we use a method based on the change in phase in the MR image arising from the field perturbation caused by deoxygenated haemoglobin in veins. We build on the existing phase based approach (Method I), where Yv is measured in a large vein (such as the superior sagittal sinus) based on the field shift inside the vein with assumptions as to the vein's shape and orientation. We demonstrate two novel modifications which address limitations of this method. The first modification (Method II), maps the actual form of the vein, rather than assume a given shape and orientation. The second modification (Method III) uses the intra and perivascular phase change in response to a known change in Yv on hyperoxia to measure normoxic Yv in smaller veins. Method III can be applied to veins whose shape, size and orientation are not accurately known, thus allowing more localised measures of venous oxygenation. Results demonstrate that the use of an overly fine spatial filter caused an overestimation in Yv for Method I, whilst the measurement of Yv using Method II was less sensitive to this bias, giving Yv=0.62±0.03. Method III was applied to mapping of Yv in local veins across the brain, yielding a distribution of values with a mode of Yv=0.661±0.008