Effects of phase regression on high-resolution functional MRI of the primary visual cortex

High-resolution functional MRI studies have become a powerful tool to non-invasively probe the sub-millimeter functional organization of the human cortex. Advances in MR hardware, imaging techniques and sophisticated post-processing methods have allowed high resolution fMRI to be used in both the clinical and academic neurosciences. However, consensus within the community regarding the use of gradient echo (GE) or spin echo (SE) based acquisition remains largely divided. On one hand, GE provides a high temporal signal-to-noise ratio (tSNR) technique sensitive to both the macro- and micro-vascular signal while SE based methods are more specific to microvasculature but suffer from lower tSNR and specific absorption rate limitations, especially at high field and short repetition times. Fortunately, the phase of the GE-EPI signal is sensitive to vessel size and this provides a potential avenue to reduce the macrovascular weighting of the signal (phase regression, Menon 2002). In order to determine the efficacy of this technique at high-resolution, phase regression was applied to GE-EPI timeseries and compared to SE-EPI to determine if GE-EPI's specificity to the microvascular compartment improved. To do this, functional data was collected from seven subjects on a neuro-optimized 7 T system at 800 μm isotropic resolution with both GE-EPI and SE-EPI while observing an 8Hz contrast reversing checkerboard. Phase data from the GE-EPI was used to create a microvasculature-weighted time series (GE-EPI-PR). Anatomical imaging (MP2RAGE) was also collected to allow for surface segmentation so that the functional results could be projected onto a surface. A multi-echo gradient echo sequence was collected and used to identify venous vasculature. The GE-EPI-PR surface activation maps showed a high qualitative similarity with SE-EPI and also produced laminar activity profiles similar to SE-EPI. When the GE-EPI and GE-EPI-PR distributions were compared to SE-EPI it was shown that GE-EPI-PR had similar distribution characteristics to SE-EPI (p<0.05) across the top 60% of cortex. Furthermore, it was shown that GE-EPI-PR has a higher contrast-to-noise ratio (0.5 ± 0.2, mean ± std. dev. across layers) than SE-EPI (0.27 ± 0.07) demonstrating the technique has higher sensitivity than SE-EPI. Taken together this evidence suggests phase regression is a useful method in low SNR studies such as high-resolution fMRI.

[1]  Robin M Heidemann,et al.  Generalized autocalibrating partially parallel acquisitions (GRAPPA) , 2002, Magnetic resonance in medicine.

[2]  R. Gruetter,et al.  Mapping and characterization of positive and negative BOLD responses to visual stimulation in multiple brain regions at 7T , 2018, Human brain mapping.

[3]  Laurentius Huber,et al.  Techniques for blood volume fMRI with VASO: From low-resolution mapping towards sub-millimeter layer-dependent applications , 2018, NeuroImage.

[4]  M A Fernández-Seara,et al.  Postprocessing technique to correct for background gradients in image‐based R*2 measurements , 2000, Magnetic resonance in medicine.

[5]  D. O. Walsh,et al.  Adaptive reconstruction of phased array MR imagery , 2000, Magnetic resonance in medicine.

[6]  Satrajit S. Ghosh,et al.  Nipype: A Flexible, Lightweight and Extensible Neuroimaging Data Processing Framework in Python , 2011, Front. Neuroinform..

[7]  Klaus Scheffler,et al.  Functional MRI in human subjects with gradient‐echo and spin‐echo EPI at 9.4 T , 2014, Magnetic resonance in medicine.

[8]  Essa Yacoub,et al.  Strategies and prospects for cortical depth dependent T2 and T2* weighted BOLD fMRI studies , 2019, NeuroImage.

[9]  Gang Chen,et al.  Layer-specific BOLD activation in awake monkey V1 revealed by ultra-high spatial resolution functional magnetic resonance imaging , 2013, NeuroImage.

[10]  An T. Vu,et al.  Pushing the limits of ultra-high resolution human brain imaging with SMS-EPI demonstrated for columnar level fMRI , 2018, NeuroImage.

[11]  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.

[12]  Pieter R. Roelfsema,et al.  Benchmarking laminar fMRI: Neuronal spiking and synaptic activity during top-down and bottom-up processing in the different layers of cortex , 2017, NeuroImage.

[13]  Kamil Ugurbil,et al.  A critical assessment of data quality and venous effects in sub-millimeter fMRI , 2019, NeuroImage.

[14]  César Caballero-Gaudes,et al.  Methods for cleaning the BOLD fMRI signal , 2016, NeuroImage.

[15]  Essa Yacoub,et al.  Retinotopic mapping with spin echo BOLD at 7T. , 2010, Magnetic resonance imaging.

[16]  Jonathan R. Polimeni,et al.  Analysis strategies for high-resolution UHF-fMRI data , 2017, NeuroImage.

[17]  Essa Yacoub,et al.  Sub-millimeter T2 weighted fMRI at 7 T: comparison of 3D-GRASE and 2D SE-EPI , 2015, Front. Neurosci..

[18]  Anders M. Dale,et al.  Cortical Surface-Based Analysis I. Segmentation and Surface Reconstruction , 1999, NeuroImage.

[19]  Natalia Petridou,et al.  Laminar imaging of positive and negative BOLD in human visual cortex at 7T , 2018, NeuroImage.

[20]  Klaus Scheffler,et al.  The impact of vessel size, orientation and intravascular contribution on the neurovascular fingerprint of BOLD bSSFP fMRI , 2017, NeuroImage.

[21]  D. Feinberg,et al.  Halving MR imaging time by conjugation: demonstration at 3.5 kG. , 1986, Radiology.

[22]  K. Uğurbil,et al.  Microvascular BOLD contribution at 4 and 7 T in the human brain: Gradient‐echo and spin‐echo fMRI with suppression of blood effects , 2003, Magnetic resonance in medicine.

[23]  Kay Nehrke,et al.  On the steady‐state properties of actual flip angle imaging (AFI) , 2009, Magnetic resonance in medicine.

[24]  Vasily L Yarnykh,et al.  Actual flip‐angle imaging in the pulsed steady state: A method for rapid three‐dimensional mapping of the transmitted radiofrequency field , 2007, Magnetic resonance in medicine.

[25]  Steffen Bollmann,et al.  The challenge of bias‐free coil combination for quantitative susceptibility mapping at ultra‐high field , 2018, Magnetic resonance in medicine.

[26]  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.

[27]  Lawrence L. Wald,et al.  Laminar analysis of 7T BOLD using an imposed spatial activation pattern in human V1 , 2010, NeuroImage.

[28]  A. Dale,et al.  Cortical Surface-Based Analysis II: Inflation, Flattening, and a Surface-Based Coordinate System , 1999, NeuroImage.

[29]  G. Woodman,et al.  Microcircuitry of Agranular Frontal Cortex: Testing the Generality of the Canonical Cortical Microcircuit , 2014, The Journal of Neuroscience.

[30]  Markus Barth,et al.  A cortical vascular model for examining the specificity of the laminar BOLD signal , 2016, NeuroImage.

[31]  Claudine Joëlle Gauthier,et al.  Cortical lamina-dependent blood volume changes in human brain at 7T , 2015, NeuroImage.

[32]  Thomas T. Liu,et al.  A component based noise correction method (CompCor) for BOLD and perfusion based fMRI , 2007, NeuroImage.

[33]  Ravi S. Menon Postacquisition suppression of large‐vessel BOLD signals in high‐resolution fMRI , 2002, Magnetic resonance in medicine.

[34]  K. Uğurbil,et al.  Layer-Specific fMRI Reflects Different Neuronal Computations at Different Depths in Human V1 , 2012, PloS one.

[35]  J. Gore,et al.  Enhanced phase regression with savitzky‐golay filtering for high‐resolution BOLD fMRI , 2014, Human brain mapping.

[36]  R. Douglas,et al.  Neuronal circuits of the neocortex. , 2004, Annual review of neuroscience.

[37]  Keiji Tanaka,et al.  Human Ocular Dominance Columns as Revealed by High-Field Functional Magnetic Resonance Imaging , 2001, Neuron.

[38]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[39]  Amir Shmuel,et al.  A more accurate account of the effect of k-space sampling and signal decay on the effective spatial resolution in functional MRI , 2016, bioRxiv.

[40]  R. Goebel,et al.  Cortical Depth Dependent Functional Responses in Humans at 7T: Improved Specificity with 3D GRASE , 2013, PloS one.

[41]  Laurentius Huber,et al.  Comparison of BOLD and CBV using 3D EPI and 3D GRASE for cortical layer functional MRI at 7 T , 2020, Magnetic resonance in medicine.

[42]  H. Duvernoy,et al.  Cortical blood vessels of the human brain , 1981, Brain Research Bulletin.

[43]  Ravi S. Menon,et al.  Phase based venous suppression in resting-state BOLD GE-fMRI , 2014, NeuroImage.

[44]  Denis G. Pelli,et al.  ECVP '07 Abstracts , 2007, Perception.

[45]  Essa Yacoub,et al.  Sensitivity and specificity considerations for fMRI encoding, decoding, and mapping of auditory cortex at ultra-high field , 2018, NeuroImage.

[46]  Essa Yacoub,et al.  Signal and noise characteristics of Hahn SE and GE BOLD fMRI at 7 T in humans , 2005, NeuroImage.

[47]  David G Norris,et al.  Dissociable laminar profiles of concurrent bottom-up and top-down modulation in the human visual cortex , 2018, bioRxiv.

[48]  P. Batchelor,et al.  International Society for Magnetic Resonance in Medicine , 1997 .

[49]  P. Boesiger,et al.  SENSE: Sensitivity encoding for fast MRI , 1999, Magnetic resonance in medicine.

[50]  Yves Rosseel,et al.  On the Definition of Signal-To-Noise Ratio and Contrast-To-Noise Ratio for fMRI Data , 2013, PloS one.

[51]  Martin Havlicek,et al.  A dynamical model of the laminar BOLD response , 2019, NeuroImage.

[52]  Michael Schacht Hansen,et al.  Gadgetron: An open source framework for medical image reconstruction , 2013, Magnetic resonance in medicine.

[53]  D. Norris,et al.  Layer‐specific BOLD activation in human V1 , 2010, Human brain mapping.

[54]  R W Cox,et al.  AFNI: software for analysis and visualization of functional magnetic resonance neuroimages. , 1996, Computers and biomedical research, an international journal.

[55]  Federico D'Agata,et al.  Pushing functional MRI spatial and temporal resolution further: High‐density receive arrays combined with shot‐selective 2D CAIPIRINHA for 3D echo‐planar imaging at 7 T , 2020, NMR in biomedicine.

[56]  Laurentius Huber,et al.  High-Resolution CBV-fMRI Allows Mapping of Laminar Activity and Connectivity of Cortical Input and Output in Human M1 , 2017, Neuron.

[57]  Ravi S. Menon,et al.  Slice‐by‐slice B1+ shimming at 7 T , 2012, Magnetic resonance in medicine.

[58]  D H Brainard,et al.  The Psychophysics Toolbox. , 1997, Spatial vision.

[59]  Klaus Scheffler,et al.  Impact of prospective motion correction, distortion correction methods and large vein bias on the spatial accuracy of cortical laminar fMRI at 9.4 Tesla , 2019, NeuroImage.

[60]  N. Ramsey,et al.  BOLD Specificity and Dynamics Evaluated in Humans at 7 T: Comparing Gradient-Echo and Spin-Echo Hemodynamic Responses , 2013, PloS one.

[61]  Klaus Scheffler,et al.  Dependence of resting-state fMRI fluctuation amplitudes on cerebral cortical orientation relative to the direction of B0 and anatomical axes , 2019, NeuroImage.

[62]  Ravi S. Menon,et al.  Ocular dominance in human V1 demonstrated by functional magnetic resonance imaging. , 1997, Journal of neurophysiology.

[63]  Mark W. Woolrich,et al.  FSL , 2012, NeuroImage.

[64]  Jack L. Gallant,et al.  Using a novel source-localized phase regressor technique for evaluation of the vascular contribution to semantic category area localization in BOLD fMRI , 2015, Front. Neurosci..

[65]  Laurentius Huber,et al.  Techniques for blood volume fMRI with VASO , 2017 .

[66]  Wolfgang Bogner,et al.  Combining phase images from array coils using a short echo time reference scan (COMPOSER) , 2015, Magnetic resonance in medicine.

[67]  Ravi S. Menon,et al.  Submillimeter functional localization in human striate cortex using BOLD contrast at 4 Tesla: Implications for the vascular point‐spread function , 1999, Magnetic resonance in medicine.

[68]  Jun Hua,et al.  Noninvasive functional imaging of cerebral blood volume with vascular‐space‐occupancy (VASO) MRI , 2013, NMR in biomedicine.

[69]  Essa Yacoub,et al.  Evaluating the Columnar Stability of Acoustic Processing in the Human Auditory Cortex , 2018, The Journal of Neuroscience.

[70]  J. Lund,et al.  Anatomical organization of macaque monkey striate visual cortex. , 1988, Annual review of neuroscience.

[71]  Arno Klein,et al.  A reproducible evaluation of ANTs similarity metric performance in brain image registration , 2011, NeuroImage.

[72]  M. Tosetti,et al.  Characterization of high-resolution Gradient Echo and Spin Echo EPI for fMRI in the human visual cortex at 7T. , 2017, Magnetic resonance imaging.

[73]  Junjie Liu,et al.  Laminar profiles of functional activity in the human brain , 2007, NeuroImage.

[74]  Nick Todd,et al.  Phase reconstruction from multiple coil data using a virtual reference coil , 2014, Magnetic resonance in medicine.

[75]  Kendrick Kay,et al.  A temporal decomposition method for identifying venous effects in task-based fMRI , 2020, Nature Methods.

[76]  R S Menon,et al.  Highcor: A novel data-driven regressor identification method for BOLD fMRI , 2014, NeuroImage.

[77]  Ravi S. Menon The great brain versus vein debate , 2012, NeuroImage.

[78]  Essa Yacoub,et al.  Robust detection of ocular dominance columns in humans using Hahn Spin Echo BOLD functional MRI at 7 Tesla , 2007, NeuroImage.

[79]  Ninon Burgos,et al.  New advances in the Clinica software platform for clinical neuroimaging studies , 2019 .