1
|
Báez-Yáñez MG, Siero JCW, Curcic V, van Osch MJP, Petridou N. On the influence of the vascular architecture on Gradient Echo and Spin Echo BOLD fMRI signals across cortical depth: a simulation approach based on realistic 3D vascular networks. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.30.596593. [PMID: 38853905 PMCID: PMC11160811 DOI: 10.1101/2024.05.30.596593] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2024]
Abstract
GE-BOLD contrast stands out as the predominant technique in functional MRI experiments for its high sensitivity and straightforward implementation. GE-BOLD exhibits rather similar sensitivity to vessels independent of their size at submillimeter resolution studies like those examining cortical columns and laminae. However, the presence of nonspecific macrovascular contributions poses a challenge to accurately isolate neuronal activity. SE-BOLD increases specificity towards small vessels, thereby enhancing its specificity to neuronal activity, due to the effective suppression of extravascular contributions caused by macrovessels with its refocusing pulse. However, even SE-BOLD measurements may not completely remove these macrovascular contributions. By simulating hemodynamic signals across cortical depth, we gain insights into vascular contributions to the laminar BOLD signal. In this study, we employed four realistic 3D vascular models to simulate oxygen saturation states in various vascular compartments, aiming to characterize both intravascular and extravascular contributions to GE and SE signals, and corresponding BOLD signal changes, across cortical depth at 7T. Simulations suggest that SE-BOLD cannot completely reduce the macrovascular contribution near the pial surface. Simulations also show that both the specificity and signal amplitude of BOLD signals at 7T depend on the spatial arrangement of large vessels throughout cortical depth and on the pial surface.
Collapse
Affiliation(s)
- Mario Gilberto Báez-Yáñez
- Translational Neuroimaging Group, Center for Image Sciences, University Medical Center Utrecht, Utrecht, the Netherlands
| | - Jeroen C W Siero
- Translational Neuroimaging Group, Center for Image Sciences, University Medical Center Utrecht, Utrecht, the Netherlands
- Spinoza Centre for Neuroimaging Amsterdam, Amsterdam, the Netherlands
| | - Vanja Curcic
- Translational Neuroimaging Group, Center for Image Sciences, University Medical Center Utrecht, Utrecht, the Netherlands
| | - Matthias J P van Osch
- C.J. Gorter MRI Center, Department of Radiology, Leiden University Medical Center, Leiden, the Netherlands
| | - Natalia Petridou
- Translational Neuroimaging Group, Center for Image Sciences, University Medical Center Utrecht, Utrecht, the Netherlands
| |
Collapse
|
2
|
Báez-Yáñez MG, Schellekens W, Bhogal AA, Roefs ECA, van Osch MJP, Siero JCW, Petridou N. A fully synthetic three-dimensional human cerebrovascular model based on histological characteristics to investigate the hemodynamic fingerprint of the layer BOLD fMRI signal formation. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.24.595716. [PMID: 38826311 PMCID: PMC11142244 DOI: 10.1101/2024.05.24.595716] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2024]
Abstract
Recent advances in functional magnetic resonance imaging (fMRI) at ultra-high field (≥7 tesla), novel hardware, and data analysis methods have enabled detailed research on neurovascular function, such as cortical layer-specific activity, in both human and nonhuman species. A widely used fMRI technique relies on the blood oxygen level-dependent (BOLD) signal. BOLD fMRI offers insights into brain function by measuring local changes in cerebral blood volume, cerebral blood flow, and oxygen metabolism induced by increased neuronal activity. Despite its potential, interpreting BOLD fMRI data is challenging as it is only an indirect measurement of neuronal activity. Computational modeling can help interpret BOLD data by simulating the BOLD signal formation. Current developments have focused on realistic 3D vascular models based on rodent data to understand the spatial and temporal BOLD characteristics. While such rodent-based vascular models highlight the impact of the angioarchitecture on the BOLD signal amplitude, anatomical differences between the rodent and human vasculature necessitate the development of human-specific models. Therefore, a computational framework integrating human cortical vasculature, hemodynamic changes, and biophysical properties is essential. Here, we present a novel computational approach: a three-dimensional VAscular MOdel based on Statistics (3D VAMOS), enabling the investigation of the hemodynamic fingerprint of the BOLD signal within a model encompassing a fully synthetic human 3D cortical vasculature and hemodynamics. Our algorithm generates microvascular and macrovascular architectures based on morphological and topological features from the literature on human cortical vasculature. By simulating specific oxygen saturation states and biophysical interactions, our framework characterizes the intravascular and extravascular signal contributions across cortical depth and voxel-wise levels for gradient-echo and spin-echo readouts. Thereby, the 3D VAMOS computational framework demonstrates that using human characteristics significantly affects the BOLD fingerprint, making it an essential step in understanding the fundamental underpinnings of layer-specific fMRI experiments.
Collapse
Affiliation(s)
- Mario Gilberto Báez-Yáñez
- Translational Neuroimaging Group, Center for Image Sciences, University Medical Center Utrecht, Utrecht, the Netherlands
| | - Wouter Schellekens
- Translational Neuroimaging Group, Center for Image Sciences, University Medical Center Utrecht, Utrecht, the Netherlands
- Donders Centre for Cognitive Neuroimaging, Radboud UMC, Nijmegen, Netherlands
| | - Alex A Bhogal
- Translational Neuroimaging Group, Center for Image Sciences, University Medical Center Utrecht, Utrecht, the Netherlands
| | - Emiel C A Roefs
- C.J. Gorter MRI Center, Department of Radiology, Leiden University Medical Center, Leiden, the Netherlands
| | - Matthias J P van Osch
- C.J. Gorter MRI Center, Department of Radiology, Leiden University Medical Center, Leiden, the Netherlands
| | - Jeroen C W Siero
- Translational Neuroimaging Group, Center for Image Sciences, University Medical Center Utrecht, Utrecht, the Netherlands
- Spinoza Centre for Neuroimaging Amsterdam, Amsterdam, the Netherlands
| | - Natalia Petridou
- Translational Neuroimaging Group, Center for Image Sciences, University Medical Center Utrecht, Utrecht, the Netherlands
| |
Collapse
|
3
|
Choi S, Hike D, Pohmann R, Avdievich N, Gomez-Cid L, Man W, Scheffler K, Yu X. Alpha-180 spin-echo based line-scanning method for high resolution laminar-specific fMRI. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2023.05.09.540065. [PMID: 37214920 PMCID: PMC10197646 DOI: 10.1101/2023.05.09.540065] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
Laminar-specific functional magnetic resonance imaging (fMRI) has been widely used to study circuit-specific neuronal activity by mapping spatiotemporal fMRI response patterns across cortical layers. Hemodynamic responses reflect indirect neuronal activity given limit of spatial and temporal resolution. Previous gradient-echo based line-scanning fMRI (GELINE) method was proposed with high temporal (50 ms) and spatial (50 µm) resolution to better characterize the fMRI onset time across cortical layers by employing 2 saturation RF pulses. However, the imperfect RF saturation performance led to poor boundary definition of the reduced region of interest (ROI) and aliasing problems outside of the ROI. Here, we propose α (alpha)-180 spin-echo-based line-scanning fMRI (SELINE) method to resolve this issue by employing a refocusing 180° RF pulse perpendicular to the excitation slice. In contrast to GELINE signals peaked at the superficial layer, we detected varied peaks of laminar-specific BOLD signals across deeper cortical layers with the SELINE method, indicating the well-defined exclusion of the large drain-vein effect with the spin-echo sequence. Furthermore, we applied the SELINE method with 200 ms TR to sample the fast hemodynamic changes across cortical layers with a less draining vein effect. In summary, this SELINE method provides a novel acquisition scheme to identify microvascular-sensitive laminar-specific BOLD responses across cortical depth.
Collapse
|
4
|
Chen HF, Lambers H, Nagelmann N, Sandbrink M, Segelcke D, Pogatzki-Zahn E, Faber C, Pradier B. Generation of a whole-brain hemodynamic response function and sex-specific differences in cerebral processing of mechano-sensation in mice detected by BOLD fMRI. Front Neurosci 2023; 17:1187328. [PMID: 37700753 PMCID: PMC10493293 DOI: 10.3389/fnins.2023.1187328] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Accepted: 07/05/2023] [Indexed: 09/14/2023] Open
Abstract
BOLD fMRI has become a prevalent method to study cerebral sensory processing in rodent disease models, including pain and mechanical hypersensitivity. fMRI data analysis is frequently combined with a general-linear-model (GLM) -based analysis, which uses the convolution of a hemodynamic response function (HRF) with the stimulus paradigm. However, several studies indicated that the HRF differs across species, sexes, brain structures, and experimental factors, including stimulation modalities or anesthesia, and hence might strongly affect the outcome of BOLD analyzes. While considerable work has been done in humans and rats to understand the HRF, much less is known in mice. As a prerequisite to investigate mechano-sensory processing and BOLD fMRI data in male and female mice, we (1) designed a rotating stimulator that allows application of two different mechanical modalities, including innocuous von Frey and noxious pinprick stimuli and (2) determined and statistically compared HRFs across 30 brain structures and experimental conditions, including sex and, stimulus modalities. We found that mechanical stimulation lead to brain-wide BOLD signal changes thereby allowing extraction of HRFs from multiple brain structures. However, we did not find differences in HRFs across all brain structures and experimental conditions. Hence, we computed a whole-brain mouse HRF, which is based on 88 functional scans from 30 mice. A comparison of this mouse-specific HRF with our previously reported rat-derived HRF showed significantly slower kinetics in mice. Finally, we detected pronounced differences in cerebral BOLD activation between male and female mice with mechanical stimulation, thereby exposing divergent processing of noxious and innocuous stimuli in both sexes.
Collapse
Affiliation(s)
- Hui-Fen Chen
- Clinic of Radiology, Translational Research Imaging Center (TRIC), University of Münster, Münster, Germany
| | - Henriette Lambers
- Clinic of Radiology, Translational Research Imaging Center (TRIC), University of Münster, Münster, Germany
| | - Nina Nagelmann
- Clinic of Radiology, Translational Research Imaging Center (TRIC), University of Münster, Münster, Germany
| | - Martin Sandbrink
- Clinic of Radiology, Translational Research Imaging Center (TRIC), University of Münster, Münster, Germany
| | - Daniel Segelcke
- Department of Anesthesiology, Intensive Care and Pain Medicine, University of Münster, Münster, Germany
| | - Esther Pogatzki-Zahn
- Department of Anesthesiology, Intensive Care and Pain Medicine, University of Münster, Münster, Germany
| | - Cornelius Faber
- Clinic of Radiology, Translational Research Imaging Center (TRIC), University of Münster, Münster, Germany
| | - Bruno Pradier
- Clinic of Radiology, Translational Research Imaging Center (TRIC), University of Münster, Münster, Germany
- Department of Anesthesiology, Intensive Care and Pain Medicine, University of Münster, Münster, Germany
| |
Collapse
|
5
|
Choi S, Chen Y, Zeng H, Biswal B, Yu X. Identifying the distinct spectral dynamics of laminar-specific interhemispheric connectivity with bilateral line-scanning fMRI. J Cereb Blood Flow Metab 2023; 43:1115-1129. [PMID: 36803280 PMCID: PMC10291453 DOI: 10.1177/0271678x231158434] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/02/2022] [Revised: 11/30/2022] [Accepted: 12/12/2022] [Indexed: 02/23/2023]
Abstract
Despite extensive efforts to identify interhemispheric functional connectivity (FC) with resting-state (rs-) fMRI, correlated low-frequency rs-fMRI signal fluctuation across homotopic cortices originates from multiple sources. It remains challenging to differentiate circuit-specific FC from global regulation. Here, we developed a bilateral line-scanning fMRI method to detect laminar-specific rs-fMRI signals from homologous forepaw somatosensory cortices with high spatial and temporal resolution in rat brains. Based on spectral coherence analysis, two distinct bilateral fluctuation spectral features were identified: ultra-slow fluctuation (<0.04 Hz) across all cortical laminae versus Layer (L) 2/3-specific evoked BOLD at 0.05 Hz based on 4 s on/16 s off block design and resting-state fluctuations at 0.08-0.1 Hz. Based on the measurements of evoked BOLD signal at corpus callosum (CC), this L2/3-specific 0.05 Hz signal is likely associated with neuronal circuit-specific activity driven by the callosal projection, which dampened ultra-slow oscillation less than 0.04 Hz. Also, the rs-fMRI power variability clustering analysis showed that the appearance of L2/3-specific 0.08-0.1 Hz signal fluctuation is independent of the ultra-slow oscillation across different trials. Thus, distinct laminar-specific bilateral FC patterns at different frequency ranges can be identified by the bilateral line-scanning fMRI method.
Collapse
Affiliation(s)
- Sangcheon Choi
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Harvard Medical School, Massachusetts General Hospital, Charlestown, MA, USA
- Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| | - Yi Chen
- Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| | - Hang Zeng
- Max Planck Institute for Biological Cybernetics, Tübingen, Germany
- Graduate Training Centre of Neuroscience, University of Tübingen, Tübingen, Germany
| | - Bharat Biswal
- Department of Biomedical Engineering, NJIT, Newark, NJ, USA
| | - Xin Yu
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Harvard Medical School, Massachusetts General Hospital, Charlestown, MA, USA
| |
Collapse
|