Volumetric BOLD fMRI simulation: from neurovascular coupling to multivoxel imaging

BackgroundThe blood oxygenation-level dependent (BOLD) functional magnetic resonance imaging (fMRI) modality has been numerically simulated by calculating single voxel signals. However, the observation on single voxel signals cannot provide information regarding the spatial distribution of the signals. Specifically, a single BOLD voxel signal simulation cannot answer the fundamental question: is the magnetic resonance (MR) image a replica of its underling magnetic susceptibility source? In this paper, we address this problem by proposing a multivoxel volumetric BOLD fMRI simulation model and a susceptibility expression formula for linear neurovascular coupling process, that allow us to examine the BOLD fMRI procedure from neurovascular coupling to MR image formation.MethodsSince MRI technology only senses the magnetism property, we represent a linear neurovascular-coupled BOLD state by a magnetic susceptibility expression formula, which accounts for the parameters of cortical vasculature, intravascular blood oxygenation level, and local neuroactivity. Upon the susceptibility expression of a BOLD state, we carry out volumetric BOLD fMRI simulation by calculating the fieldmap (established by susceptibility magnetization) and the complex multivoxel MR image (by intravoxel dephasing). Given the predefined susceptibility source and the calculated complex MR image, we compare the MR magnitude (phase, respectively) image with the predefined susceptibility source (the calculated fieldmap) by spatial correlation.ResultsThe spatial correlation between the MR magnitude image and the magnetic susceptibility source is about 0.90 for the settings of TE = 30 ms, B0 = 3 T, voxel size = 100 micron, vessel radius = 3 micron, and blood volume fraction = 2%. Using these parameters value, the spatial correlation between the MR phase image and the susceptibility-induced fieldmap is close to 1.00.ConclusionOur simulation results show that the MR magnitude image is not an exact replica of the magnetic susceptibility source (spatial correlation ≈ 0.90), and that the MR phase image conforms closely with the susceptibility-induced fieldmap (spatial correlation ≈ 1.00).

[1]  E. Haacke,et al.  Theory of NMR signal behavior in magnetically inhomogeneous tissues: The static dephasing regime , 1994, Magnetic resonance in medicine.

[2]  V. Calhoun,et al.  Computed Inverse Resonance Imaging for Magnetic Susceptibility Map Reconstruction , 2012, Journal of computer assisted tomography.

[3]  R W Bowtell,et al.  Functional magnetic resonance imaging: imaging techniques and contrast mechanisms. , 1999, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[4]  N. Logothetis.  The neural basis of the blood-oxygen-level-dependent functional magnetic resonance imaging signal. , 2002, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[5]  R. Buxton,et al.  A Model for the Coupling between Cerebral Blood Flow and Oxygen Metabolism during Neural Stimulation , 1997, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[6]  R. Buxton.  Neuroenergetics Review Article , 2022 .

[7]  B. Rosen,et al.  Microscopic susceptibility variation and transverse relaxation: Theory and experiment , 1994, Magnetic resonance in medicine.

[8]  O. Arthurs,et al.  How well do we understand the neural origins of the fMRI BOLD signal? , 2002, Trends in Neurosciences.

[9]  Essa Yacoub,et al.  Linearity of blood-oxygenation-level dependent signal at microvasculature , 2009, NeuroImage.

[10]  C. Moonen,et al.  A fast calculation method for magnetic field inhomogeneity due to an arbitrary distribution of bulk susceptibility , 2003 .

[11]  Bin He,et al.  Noninvasive Study of Neurovascular Coupling during Graded Neuronal Suppression , 2008, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[12]  Vince D Calhoun,et al.  Magnitude and phase behavior of multiresolution BOLD signal. , 2010, Concepts in magnetic resonance. Part B, Magnetic resonance engineering.

[13]  Kamil Ugurbil,et al.  An integrative model for neuronal activity-induced signal changes for gradient and spin echo functional imaging , 2009, NeuroImage.

[14]  Ravi S. Menon,et al.  Functional brain mapping by blood oxygenation level-dependent contrast magnetic resonance imaging. A comparison of signal characteristics with a biophysical model. , 1993, Biophysical journal.

[15]  R. Bowtell,et al.  Application of a Fourier‐based method for rapid calculation of field inhomogeneity due to spatial variation of magnetic susceptibility , 2005 .

[16]  B R Rosen,et al.  Mr contrast due to intravascular magnetic susceptibility perturbations , 1995, Magnetic resonance in medicine.

[17]  Arvind Caprihan,et al.  Effect of surrounding vasculature on intravoxel BOLD signal. , 2010, Medical physics.

[18]  Bejoy Thomas,et al.  Principles, techniques, and applications of T2*-based MR imaging and its special applications. , 2009, Radiographics : a review publication of the Radiological Society of North America, Inc.

[19]  Vince D. Calhoun,et al.  Two pitfalls of BOLD fMRI magnitude-based neuroimage analysis: Non-negativity and edge effect , 2011, Journal of Neuroscience Methods.

[20]  L. Toth,et al.  How accurate is magnetic resonance imaging of brain function? , 2003, Trends in Neurosciences.

[21]  G. Crelier,et al.  Linear coupling between cerebral blood flow and oxygen consumption in activated human cortex. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[22]  Vince D. Calhoun,et al.  Voxel magnetic field disturbance from remote vasculature in BOLD fMRI , 2011, Medical Imaging.

[23]  José P Marques,et al.  Using forward calculations of the magnetic field perturbation due to a realistic vascular model to explore the BOLD effect , 2008, NMR in biomedicine.

[24]  Ruola Ning,et al.  Breast volume denoising and noise characterization by 3D wavelet transform. , 2004, Computerized medical imaging and graphics : the official journal of the Computerized Medical Imaging Society.

[25]  Vince D. Calhoun,et al.  Multiresolution voxel decomposition of complex-valued BOLD signals reveals phasor turbulence , 2011, Medical Imaging.

[26]  T. Carpenter,et al.  Linear coupling between functional magnetic resonance imaging and evoked potential amplitude in human somatosensory cortex , 2000, Neuroscience.

[27]  Seong-Gi Kim,et al.  Neural Interpretation of Blood Oxygenation Level-Dependent fMRI Maps at Submillimeter Columnar Resolution , 2007, The Journal of Neuroscience.

[28]  D. Norris.  Principles of magnetic resonance assessment of brain function , 2006, Journal of magnetic resonance imaging : JMRI.

[29]  Vince D. Calhoun,et al.  A Computational Multiresolution BOLD fMRI Model , 2011, IEEE Transactions on Biomedical Engineering.

[30]  J. R. Baker,et al.  The intravascular contribution to fmri signal change: monte carlo modeling and diffusion‐weighted studies in vivo , 1995, Magnetic resonance in medicine.

[31]  D. Attwell,et al.  The neural basis of functional brain imaging signals , 2002, Trends in Neurosciences.

[32]  B. Douglas Ward,et al.  A novel technique for modeling susceptibility-based contrast mechanisms for arbitrary microvascular geometries: The finite perturber method , 2008, NeuroImage.

[33]  Xenophon Papademetris,et al.  Rapid calculations of susceptibility-induced magnetostatic field perturbations for in vivo magnetic resonance , 2006, Physics in medicine and biology.

[34]  Ying Zheng,et al.  Theory and generalization of monte carlo models of the BOLD signal source , 2008, Magnetic resonance in medicine.

[35]  Wei Chen,et al.  Investigating the source of BOLD nonlinearity in human visual cortex in response to paired visual stimuli , 2008, NeuroImage.