1
|
Thomas ER, Haarsma J, Nicholson J, Yon D, Kok P, Press C. Predictions and errors are distinctly represented across V1 layers. Curr Biol 2024; 34:2265-2271.e4. [PMID: 38697110 DOI: 10.1016/j.cub.2024.04.036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2024] [Revised: 04/09/2024] [Accepted: 04/13/2024] [Indexed: 05/04/2024]
Abstract
Popular accounts of mind and brain propose that the brain continuously forms predictions about future sensory inputs and combines predictions with inputs to determine what we perceive.1,2,3,4,5,6 Under "predictive processing" schemes, such integration is supported by the hierarchical organization of the cortex, whereby feedback connections communicate predictions from higher-level deep layers to agranular (superficial and deep) lower-level layers.7,8,9,10 Predictions are compared with input to compute the "prediction error," which is transmitted up the hierarchy from superficial layers of lower cortical regions to the middle layers of higher areas, to update higher-level predictions until errors are reconciled.11,12,13,14,15 In the primary visual cortex (V1), predictions have thereby been proposed to influence representations in deep layers while error signals may be computed in superficial layers. Despite the framework's popularity, there is little evidence for these functional distinctions because, to our knowledge, unexpected sensory events have not previously been presented in human laminar paradigms to contrast against expected events. To this end, this 7T fMRI study contrasted V1 responses to expected (75% likely) and unexpected (25%) Gabor orientations. Multivariate decoding analyses revealed an interaction between expectation and layer, such that expected events could be decoded with comparable accuracy across layers, while unexpected events could only be decoded in superficial laminae. Although these results are in line with these accounts that have been popular for decades, such distinctions have not previously been demonstrated in humans. We discuss how both prediction and error processes may operate together to shape our unitary perceptual experiences.
Collapse
Affiliation(s)
- Emily R Thomas
- Neuroscience Institute, New York University Medical Center, 435 East 30(th) Street, New York 10016, USA; Department of Psychological Sciences, Birkbeck, University of London, Malet Street, London WC1E 7HX, UK.
| | - Joost Haarsma
- Wellcome Centre for Human Neuroimaging, Institute of Neurology, University College London, 12 Queen Square, London WC1N 3AR, UK
| | - Jessica Nicholson
- Department of Psychological Sciences, Birkbeck, University of London, Malet Street, London WC1E 7HX, UK
| | - Daniel Yon
- Department of Psychological Sciences, Birkbeck, University of London, Malet Street, London WC1E 7HX, UK
| | - Peter Kok
- Wellcome Centre for Human Neuroimaging, Institute of Neurology, University College London, 12 Queen Square, London WC1N 3AR, UK
| | - Clare Press
- Department of Psychological Sciences, Birkbeck, University of London, Malet Street, London WC1E 7HX, UK; Wellcome Centre for Human Neuroimaging, Institute of Neurology, University College London, 12 Queen Square, London WC1N 3AR, UK; Department of Experimental Psychology, University College London, 26 Bedford Way, London WC1H 0AP, UK.
| |
Collapse
|
2
|
Kennedy B, Malladi SN, Tootell RBH, Nasr S. A previously undescribed scene-selective site is the key to encoding ego-motion in naturalistic environments. eLife 2024; 13:RP91601. [PMID: 38506719 PMCID: PMC10954307 DOI: 10.7554/elife.91601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/21/2024] Open
Abstract
Current models of scene processing in the human brain include three scene-selective areas: the parahippocampal place area (or the temporal place areas), the restrosplenial cortex (or the medial place area), and the transverse occipital sulcus (or the occipital place area). Here, we challenged this model by showing that at least one other scene-selective site can also be detected within the human posterior intraparietal gyrus. Despite the smaller size of this site compared to the other scene-selective areas, the posterior intraparietal gyrus scene-selective (PIGS) site was detected consistently in a large pool of subjects (n = 59; 33 females). The reproducibility of this finding was tested based on multiple criteria, including comparing the results across sessions, utilizing different scanners (3T and 7T) and stimulus sets. Furthermore, we found that this site (but not the other three scene-selective areas) is significantly sensitive to ego-motion in scenes, thus distinguishing the role of PIGS in scene perception relative to other scene-selective areas. These results highlight the importance of including finer scale scene-selective sites in models of scene processing - a crucial step toward a more comprehensive understanding of how scenes are encoded under dynamic conditions.
Collapse
Affiliation(s)
- Bryan Kennedy
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General HospitalCharlestownUnited States
| | - Sarala N Malladi
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General HospitalCharlestownUnited States
| | - Roger BH Tootell
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General HospitalCharlestownUnited States
- Department of Radiology, Harvard Medical SchoolBostonUnited States
| | - Shahin Nasr
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General HospitalCharlestownUnited States
- Department of Radiology, Harvard Medical SchoolBostonUnited States
| |
Collapse
|
3
|
Yun SD, Küppers F, Shah NJ. Submillimeter fMRI Acquisition Techniques for Detection of Laminar and Columnar Level Brain Activation. J Magn Reson Imaging 2024; 59:747-766. [PMID: 37589385 DOI: 10.1002/jmri.28911] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2023] [Revised: 07/07/2023] [Accepted: 07/07/2023] [Indexed: 08/18/2023] Open
Abstract
Since the first demonstration in the early 1990s, functional MRI (fMRI) has emerged as one of the most powerful, noninvasive neuroimaging tools to probe brain functions. Subsequently, fMRI techniques have advanced remarkably, enabling the acquisition of functional signals with a submillimeter voxel size. This innovation has opened the possibility of investigating subcortical neural activities with respect to the cortical depths or cortical columns. For this purpose, numerous previous works have endeavored to design suitable functional contrast mechanisms and dedicated imaging techniques. Depending on the choice of the functional contrast, functional signals can be detected with high sensitivity or with improved spatial specificity to the actual activation site, and the pertaining issues have been discussed in a number of earlier works. This review paper primarily aims to provide an overview of the subcortical fMRI techniques that allow the acquisition of functional signals with a submillimeter resolution. Here, the advantages and disadvantages of the imaging techniques will be described and compared. We also summarize supplementary imaging techniques that assist in the analysis of the subcortical brain activation for more accurate mapping with reduced geometric deformation. This review suggests that there is no single universally accepted method as the gold standard for subcortical fMRI. Instead, the functional contrast and the corresponding readout imaging technique should be carefully determined depending on the purpose of the study. Due to the technical limitations of current fMRI techniques, most subcortical fMRI studies have only targeted partial brain regions. As a future prospect, the spatiotemporal resolution of fMRI will be pushed to satisfy the community's need for a deeper understanding of whole-brain functions and the underlying connectivity in order to achieve the ultimate goal of a time-resolved and layer-specific spatial scale. EVIDENCE LEVEL: 1 TECHNICAL EFFICACY: Stage 1.
Collapse
Affiliation(s)
- Seong Dae Yun
- Institute of Neuroscience and Medicine 4, INM-4, Forschungszentrum Jülich, Jülich, Germany
| | - Fabian Küppers
- Institute of Neuroscience and Medicine 4, INM-4, Forschungszentrum Jülich, Jülich, Germany
- RWTH Aachen University, Aachen, Germany
- Institute of Neuroscience and Medicine 11, INM-11, JARA, Forschungszentrum Jülich, Jülich, Germany
| | - N Jon Shah
- Institute of Neuroscience and Medicine 4, INM-4, Forschungszentrum Jülich, Jülich, Germany
- Institute of Neuroscience and Medicine 11, INM-11, JARA, Forschungszentrum Jülich, Jülich, Germany
- JARA - BRAIN - Translational Medicine, Aachen, Germany
- Department of Neurology, RWTH Aachen University, Aachen, Germany
| |
Collapse
|
4
|
Nasr S, Skerswetat J, Gaier ED, Malladi SN, Kennedy B, Tootell RB, Bex P, Hunter DG. Using high-resolution functional MRI to differentiate impacts of strabismic and anisometropic amblyopia on evoked ocular dominance activity in humans. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.02.11.579855. [PMID: 38405701 PMCID: PMC10888796 DOI: 10.1101/2024.02.11.579855] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/27/2024]
Abstract
We employed high-resolution functional MRI (fMRI) to distinguish the impacts of anisometropia and strabismus (the two most frequent causes of amblyopia) on the evoked ocular dominance (OD) response. Sixteen amblyopic participants (8 females), comprising 8 individuals with strabismus, 7 with anisometropia, 1 with deprivational amblyopia, along with 8 individuals with normal visual acuity (1 female), participated in this study for whom, we measured the difference between the response to stimulation of the two eyes, across early visual areas (V1-V4). In controls, as expected from the organization of OD columns, the evoked OD response formed a striped pattern that was mostly confined to V1. Compared to controls, the OD response in amblyopic participants formed larger fused patches that extended into downstream visual areas. Moreover, both anisometropic and strabismic participants showed stronger OD responses in V1, as well as in downstream visual areas V2-V4. Although this increase was most pronounced in V1, the correlation between the OD response level and the interocular visual acuity difference (measured behaviorally) was stronger in higher-level visual areas (V2-V4). Beyond these common effects, and despite similar densities of amblyopia between the anisometropic and strabismic participants, we found a greater increase in the size of V1 portion that responded preferentially to fellow eye stimulation in anisometropic compared to strabismic individuals. We also found a greater difference between the amplitudes of the response to binocular stimulation, in those regions that responded preferentially to the fellow vs. amblyopic eye, in anisometropic compared to strabismic subjects. In contrast, strabismic subjects demonstrated increased correlation between the OD responses evoked within V1 superficial and deep cortical depths, whereas anisometropic subjects did not. These results provide some of the first direct functional evidence for distinct impacts of strabismus and anisometropia on the mesoscale functional organization of the human visual system, thus extending what was inferred previously about amblyopia from animal models.
Collapse
Affiliation(s)
- Shahin Nasr
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States
- Department of Radiology, Harvard Medical School, Boston, MA, United States
| | - Jan Skerswetat
- Department of Psychology, Northeastern University, Boston, MA, United States
| | - Eric D. Gaier
- Department of Ophthalmology, Harvard Medical School, Boston, MA, United States
- Department of Ophthalmology, Boston’s Children Hospital, Boston, MA, United States
- Picower Institute for Learning and Memory, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - Sarala N. Malladi
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States
| | - Bryan Kennedy
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States
| | - Roger B.H. Tootell
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States
- Department of Radiology, Harvard Medical School, Boston, MA, United States
| | - Peter Bex
- Department of Psychology, Northeastern University, Boston, MA, United States
| | - David G. Hunter
- Department of Ophthalmology, Harvard Medical School, Boston, MA, United States
- Department of Ophthalmology, Boston’s Children Hospital, Boston, MA, United States
| |
Collapse
|
5
|
Gomez DEP, Polimeni JR, Lewis LD. The temporal specificity of BOLD fMRI is systematically related to anatomical and vascular features of the human brain. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.02.01.578428. [PMID: 38352610 PMCID: PMC10862860 DOI: 10.1101/2024.02.01.578428] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 02/22/2024]
Abstract
The ability to detect fast responses with functional MRI depends on the speed of hemodynamic responses to neural activity, because hemodynamic responses act as a temporal low-pass filter smoothing out rapid changes. However, hemodynamic responses (their shape and timing) are highly variable across the brain and across stimuli. This heterogeneity of responses implies that the temporal specificity of fMRI signals, or the ability of fMRI to preserve fast information, should also vary substantially across the cortex. In this work we investigated how local differences in hemodynamic response timing impact the temporal specificity of fMRI. We conducted our research using ultra-high field (7T) fMRI at high spatiotemporal resolution, using the primary visual cortex (V1) as a model area for investigation. We used visual stimuli oscillating at slow and fast frequencies to probe the temporal specificity of individual voxels. As expected, we identified substantial variability in temporal specificity, with some voxels preserving their responses to fast neural activity more effectively than others. We investigated which voxels had the highest temporal specificity and related those to anatomical and vascular features of V1. We found that low temporal specificity is only weakly explained by the presence of large veins or cerebral cortical depth. Notably, however, temporal specificity depended strongly on a voxel's position along the anterior-posterior anatomical axis of V1, with voxels within the calcarine sulcus being capable of preserving close to 25% of their amplitude as the frequency of stimulation increased from 0.05-Hz to 0.20-Hz, and voxels nearest to the occipital pole preserving less than 18%. These results indicate that detection biases in high-resolution fMRI will depend on the anatomical and vascular features of the area being imaged, and that these biases will differ depending on the timing of the underlying neuronal activity. Importantly, this spatial heterogeneity of temporal specificity suggests that it could be exploited to achieve higher specificity in some locations, and that tailored data analysis strategies may help improve the detection and interpretation of fast fMRI responses.
Collapse
Affiliation(s)
- Daniel E. P. Gomez
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Boston, MA, United States
- Department of Radiology, Harvard Medical School, Boston, MA, United States
- Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - Jonathan R. Polimeni
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Boston, MA, United States
- Department of Radiology, Harvard Medical School, Boston, MA, United States
- Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - Laura D. Lewis
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Boston, MA, United States
- Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, MA, United States
- Institute for Medical Engineering and Science, Massachusetts Institute of Technology, Cambridge, MA, United States
| |
Collapse
|
6
|
Kennedy B, Malladi SN, Tootell RBH, Nasr S. A previously undescribed scene-selective site is the key to encoding ego-motion in naturalistic environments. RESEARCH SQUARE 2024:rs.3.rs-3378081. [PMID: 38260553 PMCID: PMC10802707 DOI: 10.21203/rs.3.rs-3378081/v2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/24/2024]
Abstract
Current models of scene processing in the human brain include three scene-selective areas: the Parahippocampal Place Area (or the temporal place areas; PPA/TPA), the restrosplenial cortex (or the medial place area; RSC/MPA) and the transverse occipital sulcus (or the occipital place area; TOS/OPA). Here, we challenged this model by showing that at least one other scene-selective site can also be detected within the human posterior intraparietal gyrus. Despite the smaller size of this site compared to the other scene-selective areas, the posterior intraparietal gyrus scene-selective (PIGS) site was detected consistently in a large pool of subjects (n=59; 33 females). The reproducibility of this finding was tested based on multiple criteria, including comparing the results across sessions, utilizing different scanners (3T and 7T) and stimulus sets. Furthermore, we found that this site (but not the other three scene-selective areas) is significantly sensitive to ego-motion in scenes, thus distinguishing the role of PIGS in scene perception relative to other scene-selective areas. These results highlight the importance of including finer scale scene-selective sites in models of scene processing - a crucial step toward a more comprehensive understanding of how scenes are encoded under dynamic conditions.
Collapse
Affiliation(s)
- Bryan Kennedy
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States
| | - Sarala N. Malladi
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States
| | - Roger B. H. Tootell
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States
- Department of Radiology, Harvard Medical School, Boston, MA, United States
| | - Shahin Nasr
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States
- Department of Radiology, Harvard Medical School, Boston, MA, United States
| |
Collapse
|
7
|
Kennedy B, Malladi SN, Tootell RBH, Nasr S. A previously undescribed scene-selective site is the key to encoding ego-motion in naturalistic environments. RESEARCH SQUARE 2024:rs.3.rs-3378081. [PMID: 38260553 PMCID: PMC10802707 DOI: 10.21203/rs.3.rs-3378081/v3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
Current models of scene processing in the human brain include three scene-selective areas: the Parahippocampal Place Area (or the temporal place areas; PPA/TPA), the restrosplenial cortex (or the medial place area; RSC/MPA) and the transverse occipital sulcus (or the occipital place area; TOS/OPA). Here, we challenged this model by showing that at least one other scene-selective site can also be detected within the human posterior intraparietal gyrus. Despite the smaller size of this site compared to the other scene-selective areas, the posterior intraparietal gyrus scene-selective (PIGS) site was detected consistently in a large pool of subjects (n=59; 33 females). The reproducibility of this finding was tested based on multiple criteria, including comparing the results across sessions, utilizing different scanners (3T and 7T) and stimulus sets. Furthermore, we found that this site (but not the other three scene-selective areas) is significantly sensitive to ego-motion in scenes, thus distinguishing the role of PIGS in scene perception relative to other scene-selective areas. These results highlight the importance of including finer scale scene-selective sites in models of scene processing - a crucial step toward a more comprehensive understanding of how scenes are encoded under dynamic conditions.
Collapse
Affiliation(s)
- Bryan Kennedy
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States
| | - Sarala N. Malladi
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States
| | - Roger B. H. Tootell
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States
- Department of Radiology, Harvard Medical School, Boston, MA, United States
| | - Shahin Nasr
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States
- Department of Radiology, Harvard Medical School, Boston, MA, United States
| |
Collapse
|
8
|
Feinberg DA, Beckett AJS, Vu AT, Stockmann J, Huber L, Ma S, Ahn S, Setsompop K, Cao X, Park S, Liu C, Wald LL, Polimeni JR, Mareyam A, Gruber B, Stirnberg R, Liao C, Yacoub E, Davids M, Bell P, Rummert E, Koehler M, Potthast A, Gonzalez-Insua I, Stocker S, Gunamony S, Dietz P. Next-generation MRI scanner designed for ultra-high-resolution human brain imaging at 7 Tesla. Nat Methods 2023; 20:2048-2057. [PMID: 38012321 PMCID: PMC10703687 DOI: 10.1038/s41592-023-02068-7] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2022] [Accepted: 10/09/2023] [Indexed: 11/29/2023]
Abstract
To increase granularity in human neuroimaging science, we designed and built a next-generation 7 Tesla magnetic resonance imaging scanner to reach ultra-high resolution by implementing several advances in hardware. To improve spatial encoding and increase the image signal-to-noise ratio, we developed a head-only asymmetric gradient coil (200 mT m-1, 900 T m-1s-1) with an additional third layer of windings. We integrated a 128-channel receiver system with 64- and 96-channel receiver coil arrays to boost signal in the cerebral cortex while reducing g-factor noise to enable higher accelerations. A 16-channel transmit system reduced power deposition and improved image uniformity. The scanner routinely performs functional imaging studies at 0.35-0.45 mm isotropic spatial resolution to reveal cortical layer functional activity, achieves high angular resolution in diffusion imaging and reduces acquisition time for both functional and structural imaging.
Collapse
Affiliation(s)
- David A Feinberg
- Erwin Hahn 7T MRI Laboratory, Henry H. Wheeler Brain Imaging Center, Helen Wills Neuroscience Institute, University of California, Berkeley, Berkeley, CA, USA.
- Advanced MRI Technologies, Sebastopol, CA, USA.
| | - Alexander J S Beckett
- Erwin Hahn 7T MRI Laboratory, Henry H. Wheeler Brain Imaging Center, Helen Wills Neuroscience Institute, University of California, Berkeley, Berkeley, CA, USA
- Advanced MRI Technologies, Sebastopol, CA, USA
| | - An T Vu
- Radiology and Biomedical Imaging, University of California, San Francisco, San Francisco, CA, USA
- San Francisco Veteran Affairs Health Care System, San Francisco, CA, USA
| | - Jason Stockmann
- A. A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA
- Harvard-MIT Health Sciences and Technology, MIT, Cambridge, MA, USA
| | - Laurentius Huber
- Faculty of Psychology and Neuroscience, Maastricht University, Maastricht, the Netherlands
| | | | | | - Kawin Setsompop
- Radiological Sciences Laboratory, Stanford University, Stanford, CA, USA
| | - Xiaozhi Cao
- Radiological Sciences Laboratory, Stanford University, Stanford, CA, USA
| | - Suhyung Park
- Erwin Hahn 7T MRI Laboratory, Henry H. Wheeler Brain Imaging Center, Helen Wills Neuroscience Institute, University of California, Berkeley, Berkeley, CA, USA
- Department of Computer Engineering, Chonnam National University, Gwangju, Republic of Korea
- Department of ICT Convergence System Engineering, Chonnam National University, Gwangju, Republic of Korea
| | - Chunlei Liu
- Erwin Hahn 7T MRI Laboratory, Henry H. Wheeler Brain Imaging Center, Helen Wills Neuroscience Institute, University of California, Berkeley, Berkeley, CA, USA
| | - Lawrence L Wald
- A. A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA
- Harvard-MIT Health Sciences and Technology, MIT, Cambridge, MA, USA
| | - Jonathan R Polimeni
- A. A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA
- Harvard-MIT Health Sciences and Technology, MIT, Cambridge, MA, USA
| | - Azma Mareyam
- A. A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA
- Harvard-MIT Health Sciences and Technology, MIT, Cambridge, MA, USA
| | - Bernhard Gruber
- A. A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA
- Harvard-MIT Health Sciences and Technology, MIT, Cambridge, MA, USA
- BARNLabs, Muenzkirchen, Austria
| | | | - Congyu Liao
- Radiological Sciences Laboratory, Stanford University, Stanford, CA, USA
| | - Essa Yacoub
- Center for Magnetic Resonance Research, University of Minnesota, Minneapolis, MN, USA
| | - Mathias Davids
- A. A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA
- Harvard-MIT Health Sciences and Technology, MIT, Cambridge, MA, USA
- Computer Assisted Clinical Medicine, Medical Faculty Mannheim, Heidelberg University, Heidelberg, Germany
| | - Paul Bell
- Siemens Medical Solutions, Malvern, PA, USA
| | | | | | | | | | | | - Shajan Gunamony
- Imaging Centre of Excellence, University of Glasgow, Glasgow, UK
- MR CoilTech Limited, Glasgow, UK
| | | |
Collapse
|
9
|
Chang WT, Lin W, Giovanello KS. Enabling brain-wide mapping of directed functional connectivity at 3T via layer-dependent fMRI with draining-vein suppression. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.10.24.563835. [PMID: 37961360 PMCID: PMC10634801 DOI: 10.1101/2023.10.24.563835] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/15/2023]
Abstract
Layer-dependent functional magnetic resonance imaging (fMRI) offers a compelling avenue for investigating directed functional connectivity (FC). To construct a comprehensive map of brain-wide directed FC, several technical criteria must be met, including sub-mm spatial resolution, adequate temporal resolution, functional sensitivity, global brain coverage, and high spatial specificity. Although gradient echo (GE)-based echo planar imaging (EPI) is commonly used for rapid fMRI acquisition, it faces significant challenges due to the draining-vein effect, particularly when utilizing blood oxygen level-dependent (BOLD) contrast. In this study, we mitigated this effect by incorporating velocity-nulling (VN) gradients into a GE-BOLD fMRI sequence, opting for a 3T magnetic field strength over 7T. We also integrated several advanced techniques, such as simultaneous multi-slice (SMS) acceleration and NORDIC denoising, to enhance temporal resolution, spatial coverage, and signal sensitivity. Collectively, the VN fMRI method exhibited notable spatial specificity, as evidenced by the identification of double-peak activation patterns within the primary motor cortex (M1) during a finger-tapping task. Additionally, the technique demonstrated BOLD sensitivity in the lateral geniculate nucleus (LGN). Furthermore, our VN fMRI technique displayed superior robustness when compared to conventional fMRI approaches across participants. Our findings of directed FC elucidate several layer-specific functional relationships between different brain regions and align closely with existing literature. Given the widespread availability of 3T scanners, this technical advancement has the potential for significant impact across multiple domains of neuroscience research.
Collapse
Affiliation(s)
- Wei-Tang Chang
- Biomedical Research Imaging Center, University of North Carolina at Chapel Hill, NC, USA
- Department of Radiology, University of North Carolina at Chapel Hill, NC, USA
- Department of Biomedical Engineering, University of North Carolina at Chapel Hill, NC, USA
| | - Weili Lin
- Biomedical Research Imaging Center, University of North Carolina at Chapel Hill, NC, USA
- Department of Radiology, University of North Carolina at Chapel Hill, NC, USA
| | - Kelly S. Giovanello
- Biomedical Research Imaging Center, University of North Carolina at Chapel Hill, NC, USA
- Department of Psychology & Neuroscience, University of North Carolina at Chapel Hill, NC, USA
| |
Collapse
|
10
|
Olman CA. What multiplexing means for the interpretation of functional MRI data. Front Hum Neurosci 2023; 17:1134811. [PMID: 37091812 PMCID: PMC10117671 DOI: 10.3389/fnhum.2023.1134811] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2022] [Accepted: 03/20/2023] [Indexed: 04/08/2023] Open
Abstract
Despite technology advances that have enabled routine acquisition of functional MRI data with sub-millimeter resolution, the inferences that cognitive neuroscientists must make to link fMRI data to behavior are complicated. Thus, a single dataset subjected to different analyses can be interpreted in different ways. This article presents two optical analogies that can be useful for framing fMRI analyses in a way that allows for multiple interpretations of fMRI data to be valid simultaneously without undermining each other. The first is reflection: when an object is reflected in a mirrored surface, it appears as if the reflected object is sharing space with the mirrored object, but of course it is not. This analogy can be a good guide for interpreting the fMRI signal, since even at sub-millimeter resolutions the signal is determined by a mixture of local and long-range neural computations. The second is refraction. If we view an object through a multi-faceted prism or gemstone, our view will change-sometimes dramatically-depending on our viewing angle. In the same way, interpretation of fMRI data (inference of underlying neuronal activity) can and should be different depending on the analysis approach. Rather than representing a weakness of the methodology, or the superiority of one approach over the other (for example, simple regression analysis versus multi-voxel pattern analysis), this is an expected consequence of how information is multiplexed in the neural networks of the brain: multiple streams of information are simultaneously present in each location. The fact that any one analysis typically shows only one view of the data also puts some parentheses around fMRI practitioners' constant search for ground truth against which to compare their data. By holding our interpretations lightly and understanding that many interpretations of the data can all be true at the same time, we do a better job of preparing ourselves to appreciate, and eventually understand, the complexity of the brain and the behavior it produces.
Collapse
Affiliation(s)
- Cheryl A. Olman
- Department of Psychology, University of Minnesota, Minneapolis, MN, United States
| |
Collapse
|
11
|
Knudsen L, Bailey CJ, Blicher JU, Yang Y, Zhang P, Lund TE. Improved sensitivity and microvascular weighting of 3T laminar fMRI with GE-BOLD using NORDIC and phase regression. Neuroimage 2023; 271:120011. [PMID: 36914107 DOI: 10.1016/j.neuroimage.2023.120011] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2022] [Revised: 03/03/2023] [Accepted: 03/09/2023] [Indexed: 03/15/2023] Open
Abstract
INTRODUCTION Functional MRI with spatial resolution in the submillimeter domain enables measurements of activation across cortical layers in humans. This is valuable as different types of cortical computations, e.g., feedforward versus feedback related activity, take place in different cortical layers. Laminar fMRI studies have almost exclusively employed 7T scanners to overcome the reduced signal stability associated with small voxels. However, such systems are relatively rare and only a subset of those are clinically approved. In the present study, we examined if the feasibility of laminar fMRI at 3T could be improved by use of NORDIC denoising and phase regression. METHODS 5 healthy subjects were scanned on a Siemens MAGNETOM Prisma 3T scanner. To assess across-session reliability, each subject was scanned in 3-8 sessions on 3-4 consecutive days. A 3D gradient echo EPI (GE-EPI) sequence was used for BOLD acquisitions (voxel size 0.82 mm isotopic, TR = 2.2 s) using a block design finger tapping paradigm. NORDIC denoising was applied to the magnitude and phase time series to overcome limitations in temporal signal-to-noise ratio (tSNR) and the denoised phase time series were subsequently used to correct for large vein contamination through phase regression. RESULTS AND CONCLUSION NORDIC denoising resulted in tSNR values comparable to or higher than commonly observed at 7T. Layer-dependent activation profiles could thus be extracted robustly, within and across sessions, from regions of interest located in the hand knob of the primary motor cortex (M1). Phase regression led to substantially reduced superficial bias in obtained layer profiles, although residual macrovascular contribution remained. We believe the present results support an improved feasibility of laminar fMRI at 3T.
Collapse
Affiliation(s)
- Lasse Knudsen
- Department of Clinical Medicine, Center of Functionally Integrative Neuroscience (CFIN), Aarhus University, Universitetsbyen 3, Aarhus C 8000, Denmark; Sino-Danish Center for Education and Research (SDC), University of Chinese Academy of Sciences, 380 Huaibeizhuang, Huairou District, Beijing 101400, PR China.
| | - Christopher J Bailey
- Department of Clinical Medicine, Center of Functionally Integrative Neuroscience (CFIN), Aarhus University, Universitetsbyen 3, Aarhus C 8000, Denmark; Sino-Danish Center for Education and Research (SDC), University of Chinese Academy of Sciences, 380 Huaibeizhuang, Huairou District, Beijing 101400, PR China
| | - Jakob U Blicher
- Department of Clinical Medicine, Center of Functionally Integrative Neuroscience (CFIN), Aarhus University, Universitetsbyen 3, Aarhus C 8000, Denmark; Department of Neurology, Aalborg University Hospital, Aalborg, Denmark
| | - Yan Yang
- Sino-Danish Center for Education and Research (SDC), University of Chinese Academy of Sciences, 380 Huaibeizhuang, Huairou District, Beijing 101400, PR China; Institute of Biophysics, Chinese Academy of Sciences, Beijing, PR China
| | - Peng Zhang
- Sino-Danish Center for Education and Research (SDC), University of Chinese Academy of Sciences, 380 Huaibeizhuang, Huairou District, Beijing 101400, PR China; Institute of Biophysics, Chinese Academy of Sciences, Beijing, PR China
| | - Torben E Lund
- Department of Clinical Medicine, Center of Functionally Integrative Neuroscience (CFIN), Aarhus University, Universitetsbyen 3, Aarhus C 8000, Denmark
| |
Collapse
|
12
|
Scheeringa R, Bonnefond M, van Mourik T, Jensen O, Norris DG, Koopmans PJ. Relating neural oscillations to laminar fMRI connectivity in visual cortex. Cereb Cortex 2023; 33:1537-1549. [PMID: 35512361 PMCID: PMC9977363 DOI: 10.1093/cercor/bhac154] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2021] [Revised: 03/30/2022] [Accepted: 03/31/2022] [Indexed: 11/13/2022] Open
Abstract
Laminar functional magnetic resonance imaging (fMRI) holds the potential to study connectivity at the laminar level in humans. Here we analyze simultaneously recorded electroencephalography (EEG) and high-resolution fMRI data to investigate how EEG power modulations, induced by a task with an attentional component, relate to changes in fMRI laminar connectivity between and within brain regions in visual cortex. Our results indicate that our task-induced decrease in beta power relates to an increase in deep-to-deep layer coupling between regions and to an increase in deep/middle-to-superficial layer connectivity within brain regions. The attention-related alpha power decrease predominantly relates to reduced connectivity between deep and superficial layers within brain regions, since, unlike beta power, alpha power was found to be positively correlated to connectivity. We observed no strong relation between laminar connectivity and gamma band oscillations. These results indicate that especially beta band, and to a lesser extent, alpha band oscillations relate to laminar-specific fMRI connectivity. The differential effects for alpha and beta bands indicate that they relate to different feedback-related neural processes that are differentially expressed in intra-region laminar fMRI-based connectivity.
Collapse
Affiliation(s)
- René Scheeringa
- Erwin L. Hahn Institute for Magnetic Resonance Imaging, UNESCO-Weltkulturerbe Zollverein, University of Duisburg-Essen, Kokereiallee 7, 45141 Essen, Germany.,High-Field and Hybrid MR Imaging, University Hospital Essen, University of Duisburg-Essen, Hufelandstr. 55, 45147 Essen, Germany.,Lyon Neuroscience Research Center; CRNL, INSERM U1028, CNRS UMR5292, University of Lyon 1, Université de Lyon, Bâtiment 462 - Neurocampus, 95 Bd Pinel, 69500 Bron, France.,Donders Centre for Cognitive Neuroimaging, Donders Institute for Brain, Cognition and Behaviour, Radboud University Nijmegen, Trigon 204, P.O. Box 9101, 6500 HB Nijmegen, The Netherlands
| | - Mathilde Bonnefond
- Lyon Neuroscience Research Center; CRNL, INSERM U1028, CNRS UMR5292, University of Lyon 1, Université de Lyon, Bâtiment 462 - Neurocampus, 95 Bd Pinel, 69500 Bron, France
| | - Tim van Mourik
- Donders Centre for Cognitive Neuroimaging, Donders Institute for Brain, Cognition and Behaviour, Radboud University Nijmegen, Trigon 204, P.O. Box 9101, 6500 HB Nijmegen, The Netherlands
| | - Ole Jensen
- School of Psychology, Centre for Human Brain Health, University of Birmingham, Hills Building, Birmingham B15 2TT, United Kingdom
| | - David G Norris
- Erwin L. Hahn Institute for Magnetic Resonance Imaging, UNESCO-Weltkulturerbe Zollverein, University of Duisburg-Essen, Kokereiallee 7, 45141 Essen, Germany.,Donders Centre for Cognitive Neuroimaging, Donders Institute for Brain, Cognition and Behaviour, Radboud University Nijmegen, Trigon 204, P.O. Box 9101, 6500 HB Nijmegen, The Netherlands
| | - Peter J Koopmans
- Erwin L. Hahn Institute for Magnetic Resonance Imaging, UNESCO-Weltkulturerbe Zollverein, University of Duisburg-Essen, Kokereiallee 7, 45141 Essen, Germany.,High-Field and Hybrid MR Imaging, University Hospital Essen, University of Duisburg-Essen, Hufelandstr. 55, 45147 Essen, Germany.,Donders Centre for Cognitive Neuroimaging, Donders Institute for Brain, Cognition and Behaviour, Radboud University Nijmegen, Trigon 204, P.O. Box 9101, 6500 HB Nijmegen, The Netherlands.,Department of Radiation Oncology, Radboud University Medical Center, Geert Grooteplein Zuid 10, 6525 GA, Nijmegen, The Netherlands
| |
Collapse
|
13
|
Kennedy B, Bex P, Hunter DG, Nasr S. Two fine-scale channels for encoding motion and stereopsis within the human magnocellular stream. Prog Neurobiol 2023; 220:102374. [PMID: 36403864 PMCID: PMC9832588 DOI: 10.1016/j.pneurobio.2022.102374] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2022] [Revised: 10/16/2022] [Accepted: 11/15/2022] [Indexed: 11/18/2022]
Abstract
In humans and non-human primates (NHPs), motion and stereopsis are processed within fine-scale cortical sites, including V2 thick stripes and their extensions into areas V3 and V3A that are believed to be under the influence of magnocellular stream. However, in both species, the relative functional organization (overlapping vs. none overlapping) of these sites remains unclear. Using high-resolution functional MRI (fMRI), we found evidence for two minimally-overlapping channels within human extrastriate areas that contribute to processing motion and stereopsis. Across multiple experiments that included different stimuli (random dots, gratings, and natural scenes), the functional selectivity of these channels for motion vs. stereopsis remained consistent. Furthermore, an analysis of resting-state functional connectivity revealed stronger functional connectivity within the two channels rather than between them. This finding provides a new perspective toward the mesoscale organization of the magnocellular stream within the human extrastriate visual cortex, beyond our previous understanding based on animal models.
Collapse
Affiliation(s)
- B Kennedy
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States
| | - P Bex
- Department of Psychology, Northeastern University, Boston, MA, United States
| | - D G Hunter
- Department of Ophthalmology, Harvard Medical School, Boston, MA, United States; Department of Ophthalmology, Boston's Children Hospital, Boston, MA, United States
| | - S Nasr
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States; Department of Radiology, Harvard Medical School, Boston, MA, United States.
| |
Collapse
|
14
|
Pfaffenrot V, Koopmans PJ. Magnetization transfer weighted laminar fMRI with multi-echo FLASH. Neuroimage 2022; 264:119725. [PMID: 36328273 DOI: 10.1016/j.neuroimage.2022.119725] [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: 08/10/2022] [Revised: 10/13/2022] [Accepted: 10/30/2022] [Indexed: 11/06/2022] Open
Abstract
Laminar functional magnetic resonance imaging (fMRI) using the gradient echo (GRE) blood oxygenation level dependent (BOLD) contrast is prone to signal changes arising from large unspecific venous vessels. Alternatives based on changes of cerebral blood volume (CBV) become more popular since it is expected that this hemodynamic response is dominant in microvasculature. One approach to sensitize the signal toward changes in CBV, and to simultaneously reduce unwanted extravascular (EV) BOLD blurring, is to selectively reduce gray matter (GM) signal via magnetization transfer (MT). In this work, we use off-resonant MT-pulses with a 3D FLASH readout to perform MT-prepared (MT-prep) laminar fMRI of the primary visual cortex (V1) at multiple echo times at 7 T. With a GRE-BOLD contrast without additional MT-weighting as reference, we investigated the influence of the MT-preparation on the shape and the echo time dependency of laminar profiles. Through numerical simulations, we optimized the sequence parameters to increase the sensitivity toward signal changes induced by changes in arterial CBV and to delineate the contributions of different compartments to the signal. We show that at 7 T, GM signals can be reduced by 30 %. Our laminar fMRI responses exhibit an increased signal change in the parenchyma at very short TE compared to a BOLD-only reference as a result of reduced EV signal intensity. By varying echo times, we could show that MT-prep results in less sensitivity toward unwanted signal changes based on changes in T2*. We conclude that when accounting for nuclear overhauser enhancement effects in blood, off-resonant MT-prep combined with efficient short TE readouts can become a promising method to reduce unwanted EV venous contributions in GRE-BOLD and/or to allow scanning at much shorter echo times without incurring a sensitivity penalty in laminar fMRI.
Collapse
Affiliation(s)
- Viktor Pfaffenrot
- Erwin L. Hahn Institute for Magnetic Resonance Imaging, University of Duisburg-Essen, 45141 Essen, Germany; High Field and Hybrid MR Imaging, University Hospital Essen, University of Duisburg-Essen, 45147 Essen, Germany.
| | - Peter J Koopmans
- Erwin L. Hahn Institute for Magnetic Resonance Imaging, University of Duisburg-Essen, 45141 Essen, Germany; High Field and Hybrid MR Imaging, University Hospital Essen, University of Duisburg-Essen, 45147 Essen, Germany; Department of Radiation Oncology, Radboud University Medical Center, Nijmegen, the Netherlands
| |
Collapse
|
15
|
Akbari A, Bollmann S, Ali TS, Barth M. Modelling the depth-dependent VASO and BOLD responses in human primary visual cortex. Hum Brain Mapp 2022; 44:710-726. [PMID: 36189837 PMCID: PMC9842911 DOI: 10.1002/hbm.26094] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2022] [Revised: 08/05/2022] [Accepted: 08/07/2022] [Indexed: 01/25/2023] Open
Abstract
Functional magnetic resonance imaging (fMRI) using a blood-oxygenation-level-dependent (BOLD) contrast is a common method for studying human brain function noninvasively. Gradient-echo (GRE) BOLD is highly sensitive to the blood oxygenation change in blood vessels; however, the spatial signal specificity can be degraded due to signal leakage from activated lower layers to superficial layers in depth-dependent (also called laminar or layer-specific) fMRI. Alternatively, physiological variables such as cerebral blood volume using the VAscular-Space-Occupancy (VASO) contrast have shown higher spatial specificity compared to BOLD. To better understand the physiological mechanisms such as blood volume and oxygenation changes and to interpret the measured depth-dependent responses, models are needed which reflect vascular properties at this scale. For this purpose, we extended and modified the "cortical vascular model" previously developed to predict layer-specific BOLD signal changes in human primary visual cortex to also predict a layer-specific VASO response. To evaluate the model, we compared the predictions with experimental results of simultaneous VASO and BOLD measurements in a group of healthy participants. Fitting the model to our experimental data provided an estimate of CBV change in different vascular compartments upon neural activity. We found that stimulus-evoked CBV change mainly occurs in small arterioles, capillaries, and intracortical arteries and that the contribution from venules and ICVs is smaller. Our results confirm that VASO is less susceptible to large vessel effects compared to BOLD, as blood volume changes in intracortical arteries did not substantially affect the resulting depth-dependent VASO profiles, whereas depth-dependent BOLD profiles showed a bias towards signal contributions from intracortical veins.
Collapse
Affiliation(s)
- Atena Akbari
- Centre for Advanced ImagingUniversity of QueenslandBrisbaneAustralia
| | - Saskia Bollmann
- Centre for Advanced ImagingUniversity of QueenslandBrisbaneAustralia
| | - Tonima S. Ali
- Centre for Advanced ImagingUniversity of QueenslandBrisbaneAustralia
| | - Markus Barth
- Centre for Advanced ImagingUniversity of QueenslandBrisbaneAustralia,ARC Training Centre for Innovation in Biomedical Imaging TechnologyThe University of QueenslandBrisbaneAustralia,School of Information Technology and Electrical EngineeringThe University of QueenslandBrisbaneQueenslandAustralia
| |
Collapse
|
16
|
Haarsma J, Kok P, Browning M. The promise of layer-specific neuroimaging for testing predictive coding theories of psychosis. Schizophr Res 2022; 245:68-76. [PMID: 33199171 PMCID: PMC9241988 DOI: 10.1016/j.schres.2020.10.009] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/30/2020] [Revised: 10/03/2020] [Accepted: 10/28/2020] [Indexed: 12/24/2022]
Abstract
Predictive coding potentially provides an explanatory model for understanding the neurocognitive mechanisms of psychosis. It proposes that cognitive processes, such as perception and inference, are implemented by a hierarchical system, with the influence of each level being a function of the estimated precision of beliefs at that level. However, predictive coding models of psychosis are insufficiently constrained-any phenomenon can be explained in multiple ways by postulating different changes to precision at different levels of processing. One reason for the lack of constraint in these models is that the core processes are thought to be implemented by the function of specific cortical layers, and the technology to measure layer specific neural activity in humans has until recently been lacking. As a result, our ability to constrain the models with empirical data has been limited. In this review we provide a brief overview of predictive processing models of psychosis and then describe the potential for newly developed, layer specific neuroimaging techniques to test and thus constrain these models. We conclude by discussing the most promising avenues for this research as well as the technical and conceptual challenges which may limit its application.
Collapse
Affiliation(s)
- J. Haarsma
- Wellcome Centre for Human Neuroimaging, University College London, London, United Kingdom,Department of Psychiatry, University of Oxford, Oxford, United Kingdom,Corresponding author at: Wellcome Centre for Human Neuroimaging, University College London, London, United Kingdom.
| | - P. Kok
- Wellcome Centre for Human Neuroimaging, University College London, London, United Kingdom
| | - M. Browning
- Department of Psychiatry, University of Oxford, Oxford, United Kingdom,Oxford Health NHS Trust, Oxford, United Kingdom
| |
Collapse
|
17
|
Pais-Roldán P, Yun SD, Shah NJ. Pre-processing of Sub-millimeter GE-BOLD fMRI Data for Laminar Applications. FRONTIERS IN NEUROIMAGING 2022; 1:869454. [PMID: 37555171 PMCID: PMC10406219 DOI: 10.3389/fnimg.2022.869454] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/04/2022] [Accepted: 03/31/2022] [Indexed: 08/10/2023]
Abstract
Over the past 30 years, brain function has primarily been evaluated non-invasively using functional magnetic resonance imaging (fMRI) with gradient-echo (GE) sequences to measure blood-oxygen-level-dependent (BOLD) signals. Despite the multiple advantages of GE sequences, e.g., higher signal-to-noise ratio, faster acquisitions, etc., their relatively inferior spatial localization compromises the routine use of GE-BOLD in laminar applications. Here, in an attempt to rescue the benefits of GE sequences, we evaluated the effect of existing pre-processing methods on the spatial localization of signals obtained with EPIK, a GE sequence that affords voxel volumes of 0.25 mm3 with near whole-brain coverage. The methods assessed here apply to both task and resting-state fMRI data assuming the availability of reconstructed magnitude and phase images.
Collapse
Affiliation(s)
- Patricia Pais-Roldán
- Institute of Neuroscience and Medicine 4, Medical Imaging Physics, Forschungszentrum Jülich, Jülich, Germany
| | - Seong Dae Yun
- Institute of Neuroscience and Medicine 4, Medical Imaging Physics, Forschungszentrum Jülich, Jülich, Germany
| | - N. Jon Shah
- Institute of Neuroscience and Medicine 4, Medical Imaging Physics, Forschungszentrum Jülich, Jülich, Germany
- Institute of Neuroscience and Medicine 11, Molecular Neuroscience and Neuroimaging, Jülich Aachen Research Alliance, Forschungszentrum Jülich, Jülich, Germany
- Jlich Aachen Research Alliance, Brain - Translational Medicine, Aachen, Germany
- Department of Neurology, Rheinisch-Westfälische Technische Hochschule Aachen University, Aachen, Germany
| |
Collapse
|
18
|
Choi S, Zeng H, Chen Y, Sobczak F, Qian C, Yu X. Laminar-specific functional connectivity mapping with multi-slice line-scanning fMRI. Cereb Cortex 2022; 32:4492-4501. [PMID: 35107125 PMCID: PMC9574235 DOI: 10.1093/cercor/bhab497] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2021] [Revised: 11/30/2021] [Accepted: 12/01/2021] [Indexed: 11/13/2022] Open
Abstract
Despite extensive studies detecting laminar functional magnetic resonance imaging (fMRI) signals to illustrate the canonical microcircuit, the spatiotemporal characteristics of laminar-specific information flow across cortical regions remain to be fully investigated in both evoked and resting conditions at different brain states. Here, we developed a multislice line-scanning fMRI (MS-LS) method to detect laminar fMRI signals in adjacent cortical regions with high spatial (50 μm) and temporal resolution (100 ms) in anesthetized rats. Across different trials, we detected either laminar-specific positive or negative blood-oxygen-level-dependent (BOLD) responses in the surrounding cortical region adjacent to the most activated cortex under the evoked condition. Specifically, in contrast to typical Layer (L) 4 correlation across different regions due to the thalamocortical projections for trials with positive BOLD, a strong correlation pattern specific in L2/3 was detected for trials with negative BOLD in adjacent regions, which indicated brain state-dependent laminar-fMRI responses based on corticocortical interaction. Also, in resting-state (rs-) fMRI study, robust lag time differences in L2/3, 4, and 5 across multiple cortices represented the low-frequency rs-fMRI signal propagation from caudal to rostral slices. In summary, our study provided a unique laminar fMRI mapping scheme to better characterize trial-specific intra- and inter-laminar functional connectivity in evoked and resting-state MS-LS.
Collapse
Affiliation(s)
- Sangcheon Choi
- Department of High-field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen 72076, Germany,Graduate Training Centre of Neuroscience, University of Tübingen, Tübingen 72074, Germany
| | - Hang Zeng
- Department of High-field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen 72076, Germany,Graduate Training Centre of Neuroscience, University of Tübingen, Tübingen 72074, Germany
| | - Yi Chen
- Department of High-field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen 72076, Germany
| | - Filip Sobczak
- Department of High-field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen 72076, Germany,Graduate Training Centre of Neuroscience, University of Tübingen, Tübingen 72074, Germany
| | - Chunqi Qian
- Address correspondence to Dr Xin Yu, 149, 13th Street, Charlestown, MA 02129, USA. ; Dr Chunqi Qian, 846 Service Rd, East Lansing, MI 48824, USA.
| | - Xin Yu
- Address correspondence to Dr Xin Yu, 149, 13th Street, Charlestown, MA 02129, USA. ; Dr Chunqi Qian, 846 Service Rd, East Lansing, MI 48824, USA.
| |
Collapse
|
19
|
Deshpande G, Wang Y, Robinson J. Resting state fMRI connectivity is sensitive to laminar connectional architecture in the human brain. Brain Inform 2022; 9:2. [PMID: 35038072 PMCID: PMC8764001 DOI: 10.1186/s40708-021-00150-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2021] [Accepted: 12/28/2021] [Indexed: 11/10/2022] Open
Abstract
Previous invasive studies indicate that human neocortical graymatter contains cytoarchitectonically distinct layers, with notable differences in their structural connectivity with the rest of the brain. Given recent improvements in the spatial resolution of anatomical and functional magnetic resonance imaging (fMRI), we hypothesize that resting state functional connectivity (FC) derived from fMRI is sensitive to layer-specific thalamo-cortical and cortico-cortical microcircuits. Using sub-millimeter resting state fMRI data obtained at 7 T, we found that: (1) FC between the entire thalamus and cortical layers I and VI was significantly stronger than between the thalamus and other layers. Furthermore, FC between somatosensory thalamus (ventral posterolateral nucleus, VPL) and layers IV, VI of the primary somatosensory cortex were stronger than with other layers; (2) Inter-hemispheric cortico-cortical FC between homologous regions in superficial layers (layers I-III) was stronger compared to deep layers (layers V-VI). These findings are in agreement with structural connections inferred from previous invasive studies that showed that: (i) M-type neurons in the entire thalamus project to layer-I; (ii) Pyramidal neurons in layer-VI target all thalamic nuclei, (iii) C-type neurons in the VPL project to layer-IV and receive inputs from layer-VI of the primary somatosensory cortex, and (iv) 80% of collosal projecting neurons between homologous cortical regions connect superficial layers. Our results demonstrate for the first time that resting state fMRI is sensitive to structural connections between cortical layers (previously inferred through invasive studies), specifically in thalamo-cortical and cortico-cortical networks.
Collapse
Affiliation(s)
- Gopikrishna Deshpande
- AU MRI Research Center, Department of Electrical & Computer Engineering, Auburn University, 560 Devall Dr, Suite 266D, Auburn, AL, 36849, USA. .,Department of Psychological Sciences, Auburn University, Auburn, AL, USA. .,Alabama Advanced Imaging Consortium, Birmingham, AL, USA. .,Center for Neuroscience, Auburn University, Auburn, AL, USA. .,Key Laboratory for Learning and Cognition, School of Psychology, Capital Normal University, Beijing, China. .,Department of Psychiatry, National Institute of Mental Health and Neurosciences, Bangalore, India. .,Centre for Brain Research, Indian Institute of Science, Bangalore, India.
| | - Yun Wang
- AU MRI Research Center, Department of Electrical & Computer Engineering, Auburn University, 560 Devall Dr, Suite 266D, Auburn, AL, 36849, USA.,Department of Psychiatry, Columbia University, New York, NY, USA
| | - Jennifer Robinson
- AU MRI Research Center, Department of Electrical & Computer Engineering, Auburn University, 560 Devall Dr, Suite 266D, Auburn, AL, 36849, USA.,Department of Psychological Sciences, Auburn University, Auburn, AL, USA.,Alabama Advanced Imaging Consortium, Birmingham, AL, USA.,Center for Neuroscience, Auburn University, Auburn, AL, USA
| |
Collapse
|
20
|
Lacy TC, Robinson PA, Aquino KM, Pang JC. Cortical depth-dependent modeling of visual hemodynamic responses. J Theor Biol 2021; 535:110978. [PMID: 34952032 DOI: 10.1016/j.jtbi.2021.110978] [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: 07/27/2021] [Revised: 10/18/2021] [Accepted: 12/01/2021] [Indexed: 11/15/2022]
Abstract
A physiologically based three-dimensional (3D) hemodynamic model is developed to predict the experimentally observed blood oxygen level dependent (BOLD) responses versus the cortical depth induced by visual stimuli. Prior 2D approximations are relaxed in order to analyze 3D blood flow dynamics as a function of cortical depth. Comparison of the predictions with experimental data for evoked stimuli demonstrates that the full 3D model performs at least as well as previous approaches while remaining parsimonious. In particular, the 3D model requires significantly fewer assumptions and model parameters than previous models such that there is no longer need to define depth-specific parameter values for spatial spreading, peak amplitude, and hemodynamic velocity.
Collapse
Affiliation(s)
- Thomas C Lacy
- School of Physics, University of Sydney, New South Wales, Australia; Center for Integrative Brain Function, University of Sydney, New South Wales, Australia
| | - Peter A Robinson
- School of Physics, University of Sydney, New South Wales, Australia; Center for Integrative Brain Function, University of Sydney, New South Wales, Australia
| | - Kevin M Aquino
- School of Physics, University of Sydney, New South Wales, Australia; The Turner Institute for Brain and Mental Health, School of Psychological Sciences, and Monash Biomedical Imaging, Monash University, Victoria, Australia
| | - James C Pang
- The Turner Institute for Brain and Mental Health, School of Psychological Sciences, and Monash Biomedical Imaging, Monash University, Victoria, Australia; QIMR Berghofer Medical Research Institute, Queensland, Australia.
| |
Collapse
|
21
|
The Global Configuration of Visual Stimuli Alters Co-Fluctuations of Cross-Hemispheric Human Brain Activity. J Neurosci 2021; 41:9756-9766. [PMID: 34663628 DOI: 10.1523/jneurosci.3214-20.2021] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/26/2020] [Revised: 09/11/2021] [Accepted: 10/07/2021] [Indexed: 11/21/2022] Open
Abstract
We tested how a stimulus gestalt, defined by the neuronal interaction between local and global features of a stimulus, is represented within human primary visual cortex (V1). We used high-resolution fMRI, which serves as a surrogate of neuronal activation, to measure co-fluctuations within subregions of V1 as (male and female) subjects were presented with peripheral stimuli, each with different global configurations. We found stronger cross-hemisphere correlations when fine-scale V1 cortical subregions represented parts of the same object compared with different objects. This result was consistent with the vertical bias in global processing and, critically, was independent of the task and local discontinuities within objects. Thus, despite the relatively small receptive fields of neurons within V1, global stimulus configuration affects neuronal processing via correlated fluctuations between regions that represent different sectors of the visual field.SIGNIFICANCE STATEMENT We provide the first evidence for the impact of global stimulus configuration on cross-hemispheric fMRI fluctuations, measured in human primary visual cortex. Our results are consistent with changes in the level of γ-band synchrony, which has been shown to be affected by global stimulus configuration, being reflected in the level fMRI co-fluctuations. These data help narrow the gap between knowledge of global stimulus configuration encoding at the single-neuron level versus at the behavioral level.
Collapse
|
22
|
Fracasso A, Dumoulin SO, Petridou N. Point-spread function of the BOLD response across columns and cortical depth in human extra-striate cortex. Prog Neurobiol 2021; 207:102187. [PMID: 34798198 DOI: 10.1016/j.pneurobio.2021.102187] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Columns and layers are fundamental organizational units of the brain. Well known examples of cortical columns are the ocular dominance columns (ODCs) in primary visual cortex and the column-like stripe-based arrangement in the second visual area V2. The spatial scale of columns and layers is beyond the reach of conventional neuroimaging, but the advent of high field magnetic resonance imaging (MRI) scanners (UHF, 7 Tesla and above) has opened the possibility to acquire data at this spatial scale, in-vivo and non-invasively in humans. The most prominent non-invasive technique to measure brain function is blood oxygen level dependent (BOLD) fMRI, measuring brain activity indirectly, via changes in hemodynamics. A key determinant of the ability of high-resolution BOLD fMRI to accurately resolve columns and layers is the point-spread function (PSF) of the BOLD response in relation to the spatial extent of neuronal activity. In this study we take advantage of the stripe-based arrangement present in visual area V2, coupled with sub-millimetre anatomical and gradient-echo BOLD (GE BOLD) acquisition at 7 T to obtain PSF estimates and along cortical depth in human participants. Results show that the BOLD PSF is maximal in the superficial part of the cortex (1.78 mm), and it decreases with increasing cortical depth (0.83 mm close to white matter).
Collapse
Affiliation(s)
- Alessio Fracasso
- University of Glasgow, Institute of Neuroscience and Psychology, Glasgow, Scotland, United Kingdom.
| | - Serge O Dumoulin
- Experimental Psychology, Helmholtz Institute, Utrecht University, Utrecht, the Netherlands; Spinoza Center for Neuroimaging, Amsterdam, the Netherlands; Experimental and Applied Psychology, VU University Amsterdam, the Netherlands
| | - Natalia Petridou
- Department of Radiology, Center for Image Sciences, University Medical Center Utrecht, 3584 CX, Utrecht, the Netherlands.
| |
Collapse
|
23
|
Advances in spiral fMRI: A high-resolution study with single-shot acquisition. Neuroimage 2021; 246:118738. [PMID: 34800666 DOI: 10.1016/j.neuroimage.2021.118738] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2020] [Revised: 10/23/2021] [Accepted: 11/15/2021] [Indexed: 01/15/2023] Open
Abstract
Spiral fMRI has been put forward as a viable alternative to rectilinear echo-planar imaging, in particular due to its enhanced average k-space speed and thus high acquisition efficiency. This renders spirals attractive for contemporary fMRI applications that require high spatiotemporal resolution, such as laminar or columnar fMRI. However, in practice, spiral fMRI is typically hampered by its reduced robustness and ensuing blurring artifacts, which arise from imperfections in both static and dynamic magnetic fields. Recently, these limitations have been overcome by the concerted application of an expanded signal model that accounts for such field imperfections, and its inversion by iterative image reconstruction. In the challenging ultra-high field environment of 7 Tesla, where field inhomogeneity effects are aggravated, both multi-shot and single-shot 2D spiral imaging at sub-millimeter resolution was demonstrated with high depiction quality and anatomical congruency. In this work, we further these advances towards a time series application of spiral readouts, namely, single-shot spiral BOLD fMRI at 0.8 mm in-plane resolution. We demonstrate that high-resolution spiral fMRI at 7 T is not only feasible, but delivers both excellent image quality, BOLD sensitivity, and spatial specificity of the activation maps, with little artifactual blurring. Furthermore, we show the versatility of the approach with a combined in/out spiral readout at a more typical resolution (1.5 mm), where the high acquisition efficiency allows to acquire two images per shot for improved sensitivity by echo combination.
Collapse
|
24
|
Piantoni G, Hermes D, Ramsey N, Petridou N. Size of the spatial correlation between ECoG and fMRI activity. Neuroimage 2021; 242:118459. [PMID: 34371189 PMCID: PMC10627020 DOI: 10.1016/j.neuroimage.2021.118459] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2021] [Revised: 07/13/2021] [Accepted: 08/04/2021] [Indexed: 10/20/2022] Open
Abstract
Electrocorticography (ECoG) is typically employed to accurately identify the seizure focus as well as the location of brain functions to be spared during surgical resection in participants with drug-resistant epilepsy. Increasingly, this technique has become a powerful tool to map cognitive functions onto brain regions. Cortical mapping is more commonly investigated with functional MRI (fMRI), which measures blood-oxygen level dependent (BOLD) changes induced by neuronal activity. The multimodal integration between typical 3T fMRI activity maps and ECoG measurements can provide unique insight into the spatiotemporal aspects of cognition. However, the optimal integration of fMRI and ECoG requires fundamental insight into the spatial smoothness of the BOLD signal under each electrode. Here we use ECoG as ground truth for the extent of activity, as each electrode is thought to record from the cortical tissue directly underneath the contact, to estimate the spatial smoothness of the associated BOLD response at 3T fMRI. We compared the high-frequency broadband (HFB) activity recorded with ECoG while participants performed a motor task. Activity maps were obtained with fMRI at 3T for the same task in the same participant prior to surgery. We then correlated HFB power with the fMRI BOLD signal change in the area around each electrode. This latter measure was quantified by applying a 3D Gaussian kernel of varying width (sigma between 1 mm and 20 mm) to the fMRI maps including only gray-matter. We found that the correlation between HFB and BOLD activity increased sharply up to the point when the kernel width was set to 4 mm, which we defined as the kernel width of maximal spatial specificity. After this point, as the kernel width increased, the highest level of explained variance was reached at a kernel width of 9 mm for most participants. Intriguingly, maximal specificity was also limited to 4 mm for low-frequency bands, such as alpha and beta, but the kernel width with the highest explained variance was less spatially limited than the HFB. In summary, spatial specificity is limited to a kernel width of 4 mm but explained variance keeps on increasing as you average over more and more voxels containing the relatively noisy BOLD signal. Future multimodal studies should choose the kernel width based on their research goal. For maximal spatial specificity, ECoG electrodes are best compared to 3T fMRI with a kernel width of 4 mm. When optimizing the correlation between modalities, highest explained variance can be obtained at larger kernel widths of 9 mm, at the expense of spatial specificity. Finally, we release the complete pipeline so that researchers can estimate the most appropriate kernel width from their multimodal datasets.
Collapse
Affiliation(s)
- Giovanni Piantoni
- Dept Neurology & Neurosurgery, UMC Utrecht, Heidelberglaan 100, Utrecht 3584 CX, the Netherlands.
| | - Dora Hermes
- Dept Physiology & Biomedical Engineering, Mayo Clinic, Rochester, MN, United States; Dept Neurology, Mayo Clinic, Rochester, MN, United States; Dept Radiology, Mayo Clinic, Rochester, MN, United States.
| | - Nick Ramsey
- Dept Neurology & Neurosurgery, UMC Utrecht, Heidelberglaan 100, Utrecht 3584 CX, the Netherlands.
| | - Natalia Petridou
- Dept Radiology, UMC Utrecht, Heidelberglaan 100, Utrecht, the Netherlands.
| |
Collapse
|
25
|
Raimondo L, Knapen T, Oliveira ĹAF, Yu X, Dumoulin SO, van der Zwaag W, Siero JCW. A line through the brain: implementation of human line-scanning at 7T for ultra-high spatiotemporal resolution fMRI. J Cereb Blood Flow Metab 2021; 41:2831-2843. [PMID: 34415208 PMCID: PMC8756483 DOI: 10.1177/0271678x211037266] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Functional magnetic resonance imaging (fMRI) is a widely used tool in neuroscience to detect neurally evoked responses, e.g. the blood oxygenation level-dependent (BOLD) signal. Typically, BOLD fMRI has millimeter spatial resolution and temporal resolution of one to few seconds. To study the sub-millimeter structures and activity of the cortical gray matter, the field needs an fMRI method with high spatial and temporal resolution. Line-scanning fMRI achieves very high spatial resolution and high sampling rate, at the cost of a sacrifice in volume coverage. Here, we present a human line-scanning implementation on a 7T MRI system. First, we investigate the quality of the saturation pulses that suppress MR signal outside the line. Second, we established the best coil combination for reconstruction. Finally, we applied the line-scanning method in the occipital lobe during a visual stimulation task, showing BOLD responses along cortical depth, every 250 µm with a 200 ms repetition time (TR). We found a good correspondence of t-statistics values with 2D gradient-echo echo planar imaging (GE-EPI) BOLD fMRI data with the same temporal resolution and voxel volume (R = 0.6 ± 0.2). In summary, we demonstrate the feasibility of line-scanning in humans and this opens line-scanning fMRI for applications in cognitive and clinical neuroscience.
Collapse
Affiliation(s)
- Luisa Raimondo
- Spinoza Centre for Neuroimaging, Amsterdam, Netherlands.,Experimental and Applied Psychology, VU University, Amsterdam, Netherlands
| | - Tomas Knapen
- Spinoza Centre for Neuroimaging, Amsterdam, Netherlands.,Experimental and Applied Psychology, VU University, Amsterdam, Netherlands
| | - Ĺcaro A F Oliveira
- Spinoza Centre for Neuroimaging, Amsterdam, Netherlands.,Experimental and Applied Psychology, VU University, Amsterdam, Netherlands
| | - Xin Yu
- MGH/MIT/HMS Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Harvard Medical School, Massachusetts General Hospital, Charlestown, MA, USA
| | - Serge O Dumoulin
- Spinoza Centre for Neuroimaging, Amsterdam, Netherlands.,Experimental and Applied Psychology, VU University, Amsterdam, Netherlands.,Experimental Psychology, 8125Utrecht University, Utrecht University, Utrecht, Netherlands
| | | | - Jeroen C W Siero
- Spinoza Centre for Neuroimaging, Amsterdam, Netherlands.,Radiology, Centre for Image Sciences, University Medical Centre Utrecht, Utrecht, Netherlands
| |
Collapse
|
26
|
Platt T, Ladd ME, Paech D. 7 Tesla and Beyond: Advanced Methods and Clinical Applications in Magnetic Resonance Imaging. Invest Radiol 2021; 56:705-725. [PMID: 34510098 PMCID: PMC8505159 DOI: 10.1097/rli.0000000000000820] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2021] [Revised: 08/07/2021] [Accepted: 08/07/2021] [Indexed: 12/15/2022]
Abstract
ABSTRACT Ultrahigh magnetic fields offer significantly higher signal-to-noise ratio, and several magnetic resonance applications additionally benefit from a higher contrast-to-noise ratio, with static magnetic field strengths of B0 ≥ 7 T currently being referred to as ultrahigh fields (UHFs). The advantages of UHF can be used to resolve structures more precisely or to visualize physiological/pathophysiological effects that would be difficult or even impossible to detect at lower field strengths. However, with these advantages also come challenges, such as inhomogeneities applying standard radiofrequency excitation techniques, higher energy deposition in the human body, and enhanced B0 field inhomogeneities. The advantages but also the challenges of UHF as well as promising advanced methodological developments and clinical applications that particularly benefit from UHF are discussed in this review article.
Collapse
Affiliation(s)
- Tanja Platt
- From the Medical Physics in Radiology, German Cancer Research Center (DKFZ)
| | - Mark E. Ladd
- From the Medical Physics in Radiology, German Cancer Research Center (DKFZ)
- Faculty of Physics and Astronomy
- Faculty of Medicine, University of Heidelberg, Heidelberg
- Erwin L. Hahn Institute for MRI, University of Duisburg-Essen, Essen
| | - Daniel Paech
- Division of Radiology, German Cancer Research Center (DKFZ), Heidelberg
- Clinic for Neuroradiology, University of Bonn, Bonn, Germany
| |
Collapse
|
27
|
Dowdle LT, Ghose G, Chen CCC, Ugurbil K, Yacoub E, Vizioli L. Statistical power or more precise insights into neuro-temporal dynamics? Assessing the benefits of rapid temporal sampling in fMRI. Prog Neurobiol 2021; 207:102171. [PMID: 34492308 DOI: 10.1016/j.pneurobio.2021.102171] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2021] [Revised: 08/09/2021] [Accepted: 09/02/2021] [Indexed: 01/25/2023]
Abstract
Functional magnetic resonance imaging (fMRI), a non-invasive and widely used human neuroimaging method, is most known for its spatial precision. However, there is a growing interest in its temporal sensitivity. This is despite the temporal blurring of neuronal events by the blood oxygen level dependent (BOLD) signal, the peak of which lags neuronal firing by 4-6 seconds. Given this, the goal of this review is to answer a seemingly simple question - "What are the benefits of increased temporal sampling for fMRI?". To answer this, we have combined fMRI data collected at multiple temporal scales, from 323 to 1000 milliseconds, with a review of both historical and contemporary temporal literature. After a brief discussion of technological developments that have rekindled interest in temporal research, we next consider the potential statistical and methodological benefits. Most importantly, we explore how fast fMRI can uncover previously unobserved neuro-temporal dynamics - effects that are entirely missed when sampling at conventional 1 to 2 second rates. With the intrinsic link between space and time in fMRI, this temporal renaissance also delivers improvements in spatial precision. Far from producing only statistical gains, the array of benefits suggest that the continued temporal work is worth the effort.
Collapse
Affiliation(s)
- Logan T Dowdle
- Center for Magnetic Resonance Research, University of Minnesota, 2021 6th St SE, Minneapolis, MN, 55455, United States; Department of Neurosurgery, University of Minnesota, 500 SE Harvard St, Minneapolis, MN, 55455, United States; Department of Neuroscience, University of Minnesota, 321 Church St SE, Minneapolis, MN, 55455, United States.
| | - Geoffrey Ghose
- Center for Magnetic Resonance Research, University of Minnesota, 2021 6th St SE, Minneapolis, MN, 55455, United States; Department of Neuroscience, University of Minnesota, 321 Church St SE, Minneapolis, MN, 55455, United States
| | - Clark C C Chen
- Department of Neurosurgery, University of Minnesota, 500 SE Harvard St, Minneapolis, MN, 55455, United States
| | - Kamil Ugurbil
- Center for Magnetic Resonance Research, University of Minnesota, 2021 6th St SE, Minneapolis, MN, 55455, United States
| | - Essa Yacoub
- Center for Magnetic Resonance Research, University of Minnesota, 2021 6th St SE, Minneapolis, MN, 55455, United States
| | - Luca Vizioli
- Center for Magnetic Resonance Research, University of Minnesota, 2021 6th St SE, Minneapolis, MN, 55455, United States; Department of Neurosurgery, University of Minnesota, 500 SE Harvard St, Minneapolis, MN, 55455, United States.
| |
Collapse
|
28
|
Pfaffenrot V, Voelker MN, Kashyap S, Koopmans PJ. Laminar fMRI using T 2-prepared multi-echo FLASH. Neuroimage 2021; 236:118163. [PMID: 34023449 DOI: 10.1016/j.neuroimage.2021.118163] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2021] [Revised: 05/03/2021] [Accepted: 05/11/2021] [Indexed: 02/06/2023] Open
Abstract
Functional magnetic resonance imaging (fMRI) using blood oxygenation level dependent (BOLD) contrast at a sub-millimeter scale is a promising technique to probe neural activity at the level of cortical layers. While gradient echo (GRE) BOLD sequences exhibit the highest sensitivity, their signal is confounded by unspecific extravascular (EV) and intravascular (IV) effects of large intracortical ascending veins and pial veins leading to a downstream blurring effect of local signal changes. In contrast, spin echo (SE) fMRI promises higher specificity towards signal changes near the microvascular compartment. However, the T2-weighted signal is typically sampled with a gradient echo readout imposing additional T2'-weighting. In this work, we used a T2-prepared (T2-prep) sequence with short GRE readouts to investigate its capability to acquire laminar fMRI data during a visual task in humans at 7 T. By varying the T2-prep echo time (TEprep) and acquiring multiple gradient echoes (TEGRE) per excitation, we studied the specificity of the sequence and the influence of possible confounding contributions to the shape of laminar fMRI profiles. By fitting and extrapolating the multi-echo GRE data to a TEGRE = 0 ms condition, we show for the first time laminar profiles free of T2'-pollution, confined to gray matter. This finding is independent of TEprep, except for the shortest one (31 ms) where hints of a remaining intravascular component can be seen. For TEGRE > 0 ms a prominent peak at the pial surface is observed that increases with longer TEGRE and dominates the shape of the profiles independent of the amount of T2-weighting. Simulations show that the peak at the pial surface is a result of static EV dephasing around pial vessels in CSF visible in GM due to partial voluming. Additionally, another, weaker, static dephasing effect is observed throughout all layers of the cortex, which is particularly obvious in the data with shortest T2-prep echo time. Our simulations show that this cannot be explained by intravascular dephasing but that it is likely caused by extravascular effects of the intracortical and pial veins. We conclude that even for TEGRE as short as 2.3 ms, the T2'-weighting added to the T2-weighting is enough to dramatically affect the laminar specificity of the BOLD signal change. However, the bulk of this corruption stems from CSF partial volume effects which can in principle be addressed by increasing the spatial resolution of the acquisition.
Collapse
Affiliation(s)
- Viktor Pfaffenrot
- Erwin L. Hahn Institute for Magnetic Resonance Imaging, University of Duisburg-Essen, 45141 Essen, Germany; High-Field and Hybrid MR Imaging, University Hospital Essen, University of Duisburg-Essen, 45147 Essen, Germany.
| | - Maximilian N Voelker
- Erwin L. Hahn Institute for Magnetic Resonance Imaging, University of Duisburg-Essen, 45141 Essen, Germany; High-Field and Hybrid MR Imaging, University Hospital Essen, University of Duisburg-Essen, 45147 Essen, Germany
| | - Sriranga Kashyap
- Department of Cognitive Neuroscience, Faculty of Psychology and Neuroscience, Maastricht University, 6229 Maastricht, Netherlands
| | - Peter J Koopmans
- Erwin L. Hahn Institute for Magnetic Resonance Imaging, University of Duisburg-Essen, 45141 Essen, Germany; High-Field and Hybrid MR Imaging, University Hospital Essen, University of Duisburg-Essen, 45147 Essen, Germany
| |
Collapse
|
29
|
Uludag K, Havlicek M. Determining laminar neuronal activity from BOLD fMRI using a generative model. Prog Neurobiol 2021; 207:102055. [PMID: 33930519 DOI: 10.1016/j.pneurobio.2021.102055] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2020] [Revised: 03/12/2021] [Accepted: 04/20/2021] [Indexed: 11/17/2022]
Abstract
Laminar fMRI using the BOLD contrast enables the non-invasive investigation of mesoscopic functional circuits in the human brain. However, the laminar neuronal activity is spatiotemporally biased in the observed cortical depth profiles of the BOLD signal. In this study, we propose a generative fMRI signal model, comprehensively covering the relationship between cortical depth-dependent changes in excitatory and inhibitory neuronal activity with the sampling of the BOLD signal with finite voxels. The generative model allowed us to investigate pertinent questions regarding the accuracy of the laminar BOLD signal relative to the neuronal activity, and we found that: a) condition differences in laminar BOLD signals may be more reflective of neuronal activity than single condition BOLD signal depth profiles; b) angular dependence of the BOLD signal induces significant signal variability, which can mask underlying activity profiles; c) even if only three neuronal depths are of interest, more BOLD signal depths should be considered in the analysis. In addition, we recommend that the laminar BOLD data should be displayed using the centroid method to appreciate its spatial distribution in the original resolution. Finally, we showed that Bayesian model inversion of the generative model can improve sensitivity and specificity of assessing depth-dependent neuronal changes both for steady-state and dynamically.
Collapse
Affiliation(s)
- Kamil Uludag
- Techna Institute & Koerner Scientist in MR Imaging, University Health Network, Toronto, Canada; Center for Neuroscience Imaging Research, Institute for Basic Science & Department of Biomedical Engineering, Sungkyunkwan University, Suwon, Republic of Korea.
| | - Martin Havlicek
- Department of Cognitive Neuroscience, Faculty of Psychology and Neuroscience, Maastricht University, Maastricht, Netherlands
| |
Collapse
|
30
|
Bennett MR, Farnell L, Gibson WG. Quantitative relations between BOLD responses, cortical energetics and impulse firing across cortical depth. Eur J Neurosci 2021; 54:4230-4245. [PMID: 33901325 DOI: 10.1111/ejn.15247] [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: 05/29/2020] [Accepted: 04/08/2021] [Indexed: 11/28/2022]
Abstract
The blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal arises as a consequence of changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen ( CMR O 2 ) that in turn are modulated by changes in neural activity. Recent advances in imaging have achieved sub-millimetre resolution and allowed investigation of the BOLD response as a function of cortical depth. Here, we adapt our previous theory relating the BOLD signal to neural activity to produce a quantitative model that incorporates venous blood draining between cortical layers. The adjustable inputs to the model are the neural activity and a parameter governing this blood draining. A three-layer version for transient neural inputs and a multi-layer version for constant or tonic neural inputs are able to account for a variety of experimental results, including negative BOLD signals.
Collapse
Affiliation(s)
- Maxwell R Bennett
- Brain and Mind Research Centre, University of Sydney, Camperdown, NSW, Australia
- Center for Mathematical Biology, University of Sydney, Sydney, NSW, Australia
| | - Leslie Farnell
- Center for Mathematical Biology, University of Sydney, Sydney, NSW, Australia
- The School of Mathematics and Statistics, University of Sydney, Sydney, NSW, Australia
| | - William G Gibson
- Center for Mathematical Biology, University of Sydney, Sydney, NSW, Australia
- The School of Mathematics and Statistics, University of Sydney, Sydney, NSW, Australia
| |
Collapse
|
31
|
Fracasso A, Dumoulin SO, Petridou N. Point-spread function of the BOLD response across columns and cortical depth in human extra-striate cortex. Prog Neurobiol 2021; 202:102034. [PMID: 33741401 DOI: 10.1016/j.pneurobio.2021.102034] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2020] [Revised: 03/08/2021] [Accepted: 03/14/2021] [Indexed: 10/21/2022]
Abstract
Columns and layers are fundamental organizational units of the brain. Well known examples of cortical columns are the ocular dominance columns (ODCs) in primary visual cortex and the column-like stripe-based arrangement in the second visual area V2. The spatial scale of columns and layers is beyond the reach of conventional neuroimaging, but the advent of high field magnetic resonance imaging (MRI) scanners (UHF, 7 T and above) has opened the possibility to acquire data at this spatial scale, in-vivo and non-invasively in humans. The most prominent non-invasive technique to measure brain function is blood oxygen level dependent (BOLD) fMRI, measuring brain activity indirectly, via changes in hemodynamics. A key determinant of the ability of high-resolution BOLD fMRI to accurately resolve columns and layers is the point-spread function (PSF) of the BOLD response in relation to the spatial extent of neuronal activity. In this study we take advantage of the stripe-based arrangement present in visual area V2, coupled with sub-millimetre anatomical and gradient-echo BOLD (GE BOLD) acquisition at 7 T to obtain PSF estimates and along cortical depth in human participants. Results show that the BOLD PSF is maximal in the superficial part of the cortex (1.78 mm), and it decreases with increasing cortical depth (0.83 mm close to white matter).
Collapse
Affiliation(s)
- Alessio Fracasso
- University of Glasgow, Institute of Neuroscience and Psychology, Glasgow, Scotland, United Kingdom.
| | - Serge O Dumoulin
- Experimental Psychology, Helmholtz Institute, Utrecht University, Utrecht, the Netherlands; Spinoza Center for Neuroimaging, Amsterdam, the Netherlands; Experimental and Applied Psychology, VU University Amsterdam, the Netherlands
| | - Natalia Petridou
- Department of Radiology, Center for Image Sciences, University Medical Center Utrecht, 3584 CX, Utrecht, the Netherlands.
| |
Collapse
|
32
|
Zooming-in on higher-level vision: High-resolution fMRI for understanding visual perception and awareness. Prog Neurobiol 2021; 207:101998. [PMID: 33497652 DOI: 10.1016/j.pneurobio.2021.101998] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2020] [Revised: 11/11/2020] [Accepted: 01/16/2021] [Indexed: 12/24/2022]
Abstract
One of the central questions in visual neuroscience is how the sparse retinal signals leaving our eyes are transformed into a rich subjective visual experience of the world. Invasive physiology studies, which offers the highest spatial resolution, have revealed many facts about the processing of simple visual features like contrast, color, and orientation, focusing on the early visual areas. At the same time, standard human fMRI studies with comparably coarser spatial resolution have revealed more complex, functionally specialized, and category-selective responses in higher visual areas. Although the visual system is the best understood among the sensory modalities, these two areas of research remain largely segregated. High-resolution fMRI opens up a possibility for linking them. On the one hand, it allows studying how the higher-level visual functions affect the fine-scale activity in early visual areas. On the other hand, it allows discovering the fine-scale functional organization of higher visual areas and exploring their functional connectivity with visual areas lower in the hierarchy. In this review, I will discuss examples of successful work undertaken in these directions using high-resolution fMRI and discuss where this method could be applied in the future to advance our understanding of the complexity of higher-level visual processing.
Collapse
|
33
|
Markuerkiaga I, Marques JP, Bains LJ, Norris DG. An in-vivo study of BOLD laminar responses as a function of echo time and static magnetic field strength. Sci Rep 2021; 11:1862. [PMID: 33479362 PMCID: PMC7820587 DOI: 10.1038/s41598-021-81249-w] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2019] [Accepted: 12/22/2020] [Indexed: 11/18/2022] Open
Abstract
Layer specific functional MRI requires high spatial resolution data. To compensate the associated poor signal to noise ratio it is common to integrate the signal from voxels at a given cortical depth. If the region is sufficiently large then physiological noise will be the dominant noise source. In this work, activation profiles in response to the same visual stimulus are compared at 1.5 T, 3 T and 7 T using a multi-echo, gradient echo (GE) FLASH sequence, with a 0.75 mm isotropic voxel size and the cortical integration approach. The results show that after integrating over a cortical volume of 40, 60 and 100 mm3 (at 7 T, 3 T, and 1.5 T, respectively), the signal is in the physiological noise dominated regime. The activation profiles obtained are similar for equivalent echo times. BOLD-like noise is found to be the dominant source of physiological noise. Consequently, the functional contrast to noise ratio is not strongly echo-time or field-strength dependent. We conclude that laminar GE-BOLD fMRI at lower field strengths is feasible but that larger patches of cortex will need to be examined, and that the acquisition efficiency is reduced.
Collapse
Affiliation(s)
- Irati Markuerkiaga
- Donders Centre for Cognitive Neuroimaging, Radboud University, Nijmegen, The Netherlands
| | - José P Marques
- Donders Centre for Cognitive Neuroimaging, Radboud University, Nijmegen, The Netherlands
| | - Lauren J Bains
- Donders Centre for Cognitive Neuroimaging, Radboud University, Nijmegen, The Netherlands
| | - David G Norris
- Donders Centre for Cognitive Neuroimaging, Radboud University, Nijmegen, The Netherlands. .,Erwin L. Hahn Institute for Magnetic Resonance Imaging, 45141, Essen, Germany.
| |
Collapse
|
34
|
Tootell RBH, Nasr S. Scotopic Vision Is Selectively Processed in Thick-Type Columns in Human Extrastriate Cortex. Cereb Cortex 2021; 31:1163-1181. [PMID: 33073288 PMCID: PMC7786355 DOI: 10.1093/cercor/bhaa284] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2020] [Revised: 07/25/2020] [Accepted: 08/17/2020] [Indexed: 11/26/2022] Open
Abstract
In humans, visual stimuli can be perceived across an enormous range of light levels. Evidence suggests that different neural mechanisms process different subdivisions of this range. For instance, in the retina, stimuli presented at very low (scotopic) light levels activate rod photoreceptors, whereas cone photoreceptors are activated relatively more at higher (photopic) light levels. Similarly, different retinal ganglion cells are activated by scotopic versus photopic stimuli. However, in the brain, it remains unknown whether scotopic versus photopic information is: 1) processed in distinct channels, or 2) neurally merged. Using high-resolution functional magnetic resonance imaging at 7 T, we confirmed the first hypothesis. We first localized thick versus thin-type columns within areas V2, V3, and V4, based on photopic selectivity to motion versus color, respectively. Next, we found that scotopic stimuli selectively activated thick- (compared to thin-) type columns in V2 and V3 (in measurements of both overlap and amplitude) and V4 (based on overlap). Finally, we found stronger resting-state functional connections between scotopically dominated area MT with thick- (compared to thin-) type columns in areas V2, V3, and V4. We conclude that scotopic stimuli are processed in partially segregated parallel streams, emphasizing magnocellular influence, from retina through middle stages of visual cortex.
Collapse
Affiliation(s)
- Roger B H Tootell
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Boston, MA 02114, USA.,Department of Radiology, Harvard Medical School, Boston, MA 02115, USA
| | - Shahin Nasr
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Boston, MA 02114, USA.,Department of Radiology, Harvard Medical School, Boston, MA 02115, USA
| |
Collapse
|
35
|
Barron HC, Mars RB, Dupret D, Lerch JP, Sampaio-Baptista C. Cross-species neuroscience: closing the explanatory gap. Philos Trans R Soc Lond B Biol Sci 2021; 376:20190633. [PMID: 33190601 PMCID: PMC7116399 DOI: 10.1098/rstb.2019.0633] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/20/2020] [Indexed: 12/17/2022] Open
Abstract
Neuroscience has seen substantial development in non-invasive methods available for investigating the living human brain. However, these tools are limited to coarse macroscopic measures of neural activity that aggregate the diverse responses of thousands of cells. To access neural activity at the cellular and circuit level, researchers instead rely on invasive recordings in animals. Recent advances in invasive methods now permit large-scale recording and circuit-level manipulations with exquisite spatio-temporal precision. Yet, there has been limited progress in relating these microcircuit measures to complex cognition and behaviour observed in humans. Contemporary neuroscience thus faces an explanatory gap between macroscopic descriptions of the human brain and microscopic descriptions in animal models. To close the explanatory gap, we propose adopting a cross-species approach. Despite dramatic differences in the size of mammalian brains, this approach is broadly justified by preserved homology. Here, we outline a three-armed approach for effective cross-species investigation that highlights the need to translate different measures of neural activity into a common space. We discuss how a cross-species approach has the potential to transform basic neuroscience while also benefiting neuropsychiatric drug development where clinical translation has, to date, seen minimal success. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.
Collapse
Affiliation(s)
- Helen C. Barron
- Medical Research Council Brain Network Dynamics Unit, Nuffield Department of Clinical Neurosciences, University of Oxford, Mansfield Road, Oxford OX1 3TH, UK
- Wellcome Centre for Integrative Neuroimaging, University of Oxford, FMRIB, John Radcliffe Hospital, Oxford OX3 9DU, UK
| | - Rogier B. Mars
- Wellcome Centre for Integrative Neuroimaging, University of Oxford, FMRIB, John Radcliffe Hospital, Oxford OX3 9DU, UK
- Donders Institute for Brain, Cognition and Behavior, Radboud University, 6525 AJ Nijmegen, The Netherlands
| | - David Dupret
- Medical Research Council Brain Network Dynamics Unit, Nuffield Department of Clinical Neurosciences, University of Oxford, Mansfield Road, Oxford OX1 3TH, UK
| | - Jason P. Lerch
- Wellcome Centre for Integrative Neuroimaging, University of Oxford, FMRIB, John Radcliffe Hospital, Oxford OX3 9DU, UK
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, CanadaM5G 1L7
| | - Cassandra Sampaio-Baptista
- Wellcome Centre for Integrative Neuroimaging, University of Oxford, FMRIB, John Radcliffe Hospital, Oxford OX3 9DU, UK
- Institute of Neuroscience and Psychology, University of Glasgow, Glasgow G12 8QB, UK
| |
Collapse
|
36
|
van Dijk JA, Fracasso A, Petridou N, Dumoulin SO. Validating Linear Systems Analysis for Laminar fMRI: Temporal Additivity for Stimulus Duration Manipulations. Brain Topogr 2021; 34:88-101. [PMID: 33210193 PMCID: PMC7803719 DOI: 10.1007/s10548-020-00808-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2020] [Accepted: 11/05/2020] [Indexed: 11/26/2022]
Abstract
Advancements in ultra-high field (7 T and higher) magnetic resonance imaging (MRI) scanners have made it possible to investigate both the structure and function of the human brain at a sub-millimeter scale. As neuronal feedforward and feedback information arrives in different layers, sub-millimeter functional MRI has the potential to uncover information processing between cortical micro-circuits across cortical depth, i.e. laminar fMRI. For nearly all conventional fMRI analyses, the main assumption is that the relationship between local neuronal activity and the blood oxygenation level dependent (BOLD) signal adheres to the principles of linear systems theory. For laminar fMRI, however, directional blood pooling across cortical depth stemming from the anatomy of the cortical vasculature, potentially violates these linear system assumptions, thereby complicating analysis and interpretation. Here we assess whether the temporal additivity requirement of linear systems theory holds for laminar fMRI. We measured responses elicited by viewing stimuli presented for different durations and evaluated how well the responses to shorter durations predicted those elicited by longer durations. We find that BOLD response predictions are consistently good predictors for observed responses, across all cortical depths, and in all measured visual field maps (V1, V2, and V3). Our results suggest that the temporal additivity assumption for linear systems theory holds for laminar fMRI. We thus show that the temporal additivity assumption holds across cortical depth for sub-millimeter gradient-echo BOLD fMRI in early visual cortex.
Collapse
Affiliation(s)
- Jelle A van Dijk
- Spinoza Centre for Neuroimaging, Amsterdam, The Netherlands.
- Experimental Psychology, Utrecht University, Utrecht, The Netherlands.
| | - Alessio Fracasso
- Spinoza Centre for Neuroimaging, Amsterdam, The Netherlands
- Institute of Neuroscience and Psychology, University of Glasgow, Glasgow, G12 8QB, UK
- Radiology Department, Imaging Division, Center for Image Sciences, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Natalia Petridou
- Radiology Department, Imaging Division, Center for Image Sciences, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Serge O Dumoulin
- Spinoza Centre for Neuroimaging, Amsterdam, The Netherlands
- Experimental Psychology, Utrecht University, Utrecht, The Netherlands
- Experimental and Applied Psychology, VU University, Amsterdam, The Netherlands
| |
Collapse
|
37
|
de Hollander G, van der Zwaag W, Qian C, Zhang P, Knapen T. Ultra-high field fMRI reveals origins of feedforward and feedback activity within laminae of human ocular dominance columns. Neuroimage 2020; 228:117683. [PMID: 33385565 DOI: 10.1016/j.neuroimage.2020.117683] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2020] [Revised: 11/02/2020] [Accepted: 12/14/2020] [Indexed: 11/25/2022] Open
Abstract
Ultra-high field MRI can functionally image the cerebral cortex of human subjects at the submillimeter scale of cortical columns and laminae. Here, we investigate both in concert, by imaging ocular dominance columns (ODCs) in primary visual cortex (V1) across different cortical depths. We ensured that putative ODC patterns in V1 (a) are stable across runs, sessions, and scanners located in different continents, (b) have a width (~1.3 mm) expected from post-mortem and animal work and (c) are absent at the retinotopic location of the blind spot. We then dissociated the effects of bottom-up thalamo-cortical input and attentional feedback processes on activity in V1 across cortical depth. Importantly, the separation of bottom-up information flows into ODCs allowed us to validly compare attentional conditions while keeping the stimulus identical throughout the experiment. We find that, when correcting for draining vein effects and using both model-based and model-free approaches, the effect of monocular stimulation is largest at deep and middle cortical depths. Conversely, spatial attention influences BOLD activity exclusively near the pial surface. Our findings show that simultaneous interrogation of columnar and laminar dimensions of the cortical fold can dissociate thalamocortical inputs from top-down processing, and allow the investigation of their interactions without any stimulus manipulation.
Collapse
Affiliation(s)
- Gilles de Hollander
- Department of Psychology, Vrije Universiteit Amsterdam, the Netherlands; Zurich Center for Neuroeconomics (ZNE), Department of Economics, University of Zurich, Zurich, Switzerland; Spinoza Centre for Neuroimaging, Royal Academy of Sciences, Amsterdam, the Netherlands
| | - Wietske van der Zwaag
- Spinoza Centre for Neuroimaging, Royal Academy of Sciences, Amsterdam, the Netherlands
| | - Chencan Qian
- Institute of Biophysics, Chinese Academy of Sciences, Beijing, China
| | - Peng Zhang
- Institute of Biophysics, Chinese Academy of Sciences, Beijing, China
| | - Tomas Knapen
- Department of Psychology, Vrije Universiteit Amsterdam, the Netherlands; Spinoza Centre for Neuroimaging, Royal Academy of Sciences, Amsterdam, the Netherlands
| |
Collapse
|
38
|
Stanley OW, Kuurstra AB, Klassen LM, Menon RS, Gati JS. Effects of phase regression on high-resolution functional MRI of the primary visual cortex. Neuroimage 2020; 227:117631. [PMID: 33316391 DOI: 10.1016/j.neuroimage.2020.117631] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2020] [Accepted: 12/04/2020] [Indexed: 12/14/2022] Open
Abstract
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 with 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 8 Hz 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.
Collapse
Affiliation(s)
- Olivia W Stanley
- Centre for Functional and Metabolic Mapping, Robarts Research Institute, The University of Western Ontario, London, Ontario, Canada; Medical Biophysics, The University of Western Ontario, London, Ontario, Canada.
| | - Alan B Kuurstra
- Centre for Functional and Metabolic Mapping, Robarts Research Institute, The University of Western Ontario, London, Ontario, Canada
| | - L Martyn Klassen
- Centre for Functional and Metabolic Mapping, Robarts Research Institute, The University of Western Ontario, London, Ontario, Canada; Medical Biophysics, The University of Western Ontario, London, Ontario, Canada
| | - Ravi S Menon
- Centre for Functional and Metabolic Mapping, Robarts Research Institute, The University of Western Ontario, London, Ontario, Canada; Medical Biophysics, The University of Western Ontario, London, Ontario, Canada
| | - Joseph S Gati
- Centre for Functional and Metabolic Mapping, Robarts Research Institute, The University of Western Ontario, London, Ontario, Canada; Medical Biophysics, The University of Western Ontario, London, Ontario, Canada
| |
Collapse
|
39
|
Weldon KB, Olman CA. Forging a path to mesoscopic imaging success with ultra-high field functional magnetic resonance imaging. Philos Trans R Soc Lond B Biol Sci 2020; 376:20200040. [PMID: 33190599 PMCID: PMC7741029 DOI: 10.1098/rstb.2020.0040] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Functional magnetic resonance imaging (fMRI) studies with ultra-high field (UHF, 7+ Tesla) technology enable the acquisition of high-resolution images. In this work, we discuss recent achievements in UHF fMRI at the mesoscopic scale, on the order of cortical columns and layers, and examine approaches to addressing common challenges. As researchers push to smaller and smaller voxel sizes, acquisition and analysis decisions have greater potential to degrade spatial accuracy, and UHF fMRI data must be carefully interpreted. We consider the impact of acquisition decisions on the spatial specificity of the MR signal with a representative dataset with 0.8 mm isotropic resolution. We illustrate the trade-offs in contrast with noise ratio and spatial specificity of different acquisition techniques and show that acquisition blurring can increase the effective voxel size by as much as 50% in some dimensions. We further describe how different sources of degradations to spatial resolution in functional data may be characterized. Finally, we emphasize that progress in UHF fMRI depends not only on scientific discovery and technical advancement, but also on informal discussions and documentation of challenges researchers face and overcome in pursuit of their goals. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.
Collapse
Affiliation(s)
- Kimberly B Weldon
- Department of Psychiatry and Behavioral Sciences, University of Minnesota, Minneapolis, MN 55455, USA.,Center for Magnetic Resonance Imaging, University of Minnesota, Minneapolis, MN 55455, USA
| | - Cheryl A Olman
- Center for Magnetic Resonance Imaging, University of Minnesota, Minneapolis, MN 55455, USA.,Department of Psychology, University of Minnesota, Minneapolis, MN 55455, USA
| |
Collapse
|
40
|
Bollmann S, Barth M. New acquisition techniques and their prospects for the achievable resolution of fMRI. Prog Neurobiol 2020; 207:101936. [PMID: 33130229 DOI: 10.1016/j.pneurobio.2020.101936] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2020] [Revised: 09/10/2020] [Accepted: 10/18/2020] [Indexed: 01/17/2023]
Abstract
This work reviews recent advances in technologies for functional magnetic resonance imaging (fMRI) of the human brain and highlights the push for higher functional specificity based on increased spatial resolution and specific MR contrasts to reveal previously undetectable functional properties of small-scale cortical structures. We discuss how the combination of MR hardware, advanced acquisition techniques and various MR contrast mechanisms have enabled recent progress in functional neuroimaging. However, these advanced fMRI practices have only been applied to a handful of neuroscience questions to date, with the majority of the neuroscience community still using conventional imaging techniques. We thus discuss upcoming challenges and possibilities for fMRI technology development in human neuroscience. We hope that readers interested in functional brain imaging acquire an understanding of current and novel developments and potential future applications, even if they don't have a background in MR physics or engineering. We summarize the capabilities of standard fMRI acquisition schemes with pointers to relevant literature and comprehensive reviews and introduce more recent developments.
Collapse
Affiliation(s)
- Saskia Bollmann
- Centre for Advanced Imaging, The University of Queensland, Brisbane, Australia
| | - Markus Barth
- Centre for Advanced Imaging, The University of Queensland, Brisbane, Australia; School of Information Technology and Electrical Engineering, The University of Queensland, Brisbane, Australia; ARC Training Centre for Innovation in Biomedical Imaging Technology, The University of Queensland, Brisbane, QLD, Australia.
| |
Collapse
|
41
|
Eye-selective fMRI activity in human primary visual cortex: Comparison between 3 T and 9.4 T, and effects across cortical depth. Neuroimage 2020; 220:117078. [DOI: 10.1016/j.neuroimage.2020.117078] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2019] [Revised: 06/08/2020] [Accepted: 06/18/2020] [Indexed: 12/14/2022] Open
|
42
|
Berman AJL, Grissom WA, Witzel T, Nasr S, Park DJ, Setsompop K, Polimeni JR. Ultra-high spatial resolution BOLD fMRI in humans using combined segmented-accelerated VFA-FLEET with a recursive RF pulse design. Magn Reson Med 2020; 85:120-139. [PMID: 32705723 DOI: 10.1002/mrm.28415] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2020] [Revised: 06/15/2020] [Accepted: 06/16/2020] [Indexed: 11/08/2022]
Abstract
PURPOSE To alleviate the spatial encoding limitations of single-shot echo-planar imaging (EPI) by developing multi-shot segmented EPI for ultra-high-resolution functional MRI (fMRI) with reduced ghosting artifacts from subject motion and respiration. THEORY AND METHODS Segmented EPI can reduce readout duration and reduce acceleration factors, however, the time elapsed between segment acquisitions (on the order of seconds) can result in intermittent ghosting, limiting its use for fMRI. Here, "FLEET" segment ordering, where segments are looped over before slices, was combined with a variable flip angle progression (VFA-FLEET) to improve inter-segment fidelity and maximize signal for fMRI. Scaling a sinc pulse's flip angle for each segment (VFA-FLEET-Sinc) produced inconsistent slice profiles and ghosting, therefore, a recursive Shinnar-Le Roux (SLR) radiofrequency (RF) pulse design was developed (VFA-FLEET-SLR) to generate unique pulses for every segment that together produce consistent slice profiles and signals. RESULTS The temporal stability of VFA-FLEET-SLR was compared against conventional-segmented EPI and VFA-FLEET-Sinc at 3T and 7T. VFA-FLEET-SLR showed reductions in both intermittent and stable ghosting compared to conventional-segmented and VFA-FLEET-Sinc, resulting in improved image quality with a minor trade-off in temporal SNR. Combining VFA-FLEET-SLR with acceleration, we achieved a 0.6-mm isotropic acquisition at 7T, without zoomed imaging or partial Fourier, demonstrating reliable detection of blood oxygenation level-dependent (BOLD) responses to a visual stimulus. To counteract the increased repetition time from segmentation, simultaneous multi-slice VFA-FLEET-SLR was demonstrated using RF-encoded controlled aliasing. CONCLUSIONS VFA-FLEET with a recursive RF pulse design supports acquisitions with low levels of artifact and spatial blur, enabling fMRI at previously inaccessible spatial resolutions with a "full-brain" field of view.
Collapse
Affiliation(s)
- Avery J L Berman
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, Massachusetts, USA.,Department of Radiology, Harvard Medical School, Boston, Massachusetts, USA
| | - William A Grissom
- Vanderbilt University Institute of Imaging Science, Vanderbilt University, Nashville, Tennessee, USA.,Department of Biomedical Engineering, Vanderbilt University, Nashville, Tennessee, USA
| | - Thomas Witzel
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, Massachusetts, USA.,Department of Radiology, Harvard Medical School, Boston, Massachusetts, USA
| | - Shahin Nasr
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, Massachusetts, USA.,Department of Radiology, Harvard Medical School, Boston, Massachusetts, USA
| | - Daniel J Park
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, Massachusetts, USA
| | - Kawin Setsompop
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, Massachusetts, USA.,Department of Radiology, Harvard Medical School, Boston, Massachusetts, USA.,Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
| | - Jonathan R Polimeni
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, Massachusetts, USA.,Department of Radiology, Harvard Medical School, Boston, Massachusetts, USA.,Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
| |
Collapse
|
43
|
Beckett AJS, Dadakova T, Townsend J, Huber L, Park S, Feinberg DA. Comparison of BOLD and CBV using 3D EPI and 3D GRASE for cortical layer functional MRI at 7 T. Magn Reson Med 2020; 84:3128-3145. [PMID: 32557752 DOI: 10.1002/mrm.28347] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2019] [Revised: 05/12/2020] [Accepted: 05/14/2020] [Indexed: 01/27/2023]
Abstract
PURPOSE Functional MRI (fMRI) at the mesoscale of cortical layers and columns requires both sensitivity and specificity, the latter of which can be compromised if the imaging method is affected by vascular artifacts, particularly cortical draining veins at the pial surface. Recent studies have shown that cerebral blood volume (CBV) imaging is more specific to the actual laminar locus of neural activity than BOLD imaging using standard gradient-echo EPI sequences. Gradient and spin-echo (GRASE) BOLD imaging has also shown greater specificity when compared with standard gradient-echo EPI BOLD. Here we directly compare CBV and BOLD contrasts in high-resolution imaging of the primary motor cortex for laminar functional MRI in four combinations of signal labeling, CBV using slice-selective slab-inversion vascular space occupancy (VASO) and BOLD, each with 3D gradient-echo EPI and zoomed 3D-GRASE image readouts. METHODS Activations were measured using each sequence and contrast combination during a motor task. Activation profiles across cortical depth were measured to assess the sensitivity and specificity (pial bias) of each method. RESULTS Both CBV imaging using gradient-echo 3D-EPI and BOLD imaging using 3D-GRASE show similar specificity and sensitivity and are therefore useful tools for mesoscopic functional MRI in the human cortex. The combination of GRASE and VASO did not demonstrate high levels of sensitivity, nor show increased specificity. CONCLUSION Three-dimensional EPI with VASO contrast and 3D-GRASE with BOLD contrast both demonstrate sufficient sensitivity and specificity for laminar functional MRI to be used by neuroscientists in a wide range of investigations of depth-dependent neural circuitry in the human brain.
Collapse
Affiliation(s)
- Alexander J S Beckett
- Helen Wills Neuroscience Institute, University of California, Berkeley, California, USA.,Advanced MRI Technologies, Sebastopol, California, USA
| | - Tetiana Dadakova
- Helen Wills Neuroscience Institute, University of California, Berkeley, California, USA
| | - Jennifer Townsend
- Helen Wills Neuroscience Institute, University of California, Berkeley, California, USA.,Advanced MRI Technologies, Sebastopol, California, USA
| | - Laurentius Huber
- Department of Cognitive Neuroscience, Faculty of Psychology and Neuroscience, Maastricht University, Maastricht, the Netherlands
| | - Suhyung Park
- Helen Wills Neuroscience Institute, University of California, Berkeley, California, USA.,Advanced MRI Technologies, Sebastopol, California, USA
| | - David A Feinberg
- Helen Wills Neuroscience Institute, University of California, Berkeley, California, USA.,Advanced MRI Technologies, Sebastopol, California, USA
| |
Collapse
|
44
|
Marquardt I, De Weerd P, Schneider M, Gulban OF, Ivanov D, Wang Y, Uludağ K. Feedback contribution to surface motion perception in the human early visual cortex. eLife 2020; 9:e50933. [PMID: 32496189 PMCID: PMC7314553 DOI: 10.7554/elife.50933] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2019] [Accepted: 06/03/2020] [Indexed: 01/03/2023] Open
Abstract
Human visual surface perception has neural correlates in early visual cortex, but the role of feedback during surface segmentation in human early visual cortex remains unknown. Feedback projections preferentially enter superficial and deep anatomical layers, which provides a hypothesis for the cortical depth distribution of fMRI activity related to feedback. Using ultra-high field fMRI, we report a depth distribution of activation in line with feedback during the (illusory) perception of surface motion. Our results fit with a signal re-entering in superficial depths of V1, followed by a feedforward sweep of the re-entered information through V2 and V3. The magnitude and sign of the BOLD response strongly depended on the presence of texture in the background, and was additionally modulated by the presence of illusory motion perception compatible with feedback. In summary, the present study demonstrates the potential of depth-resolved fMRI in tackling biomechanical questions on perception.
Collapse
Affiliation(s)
- Ingo Marquardt
- Department of Cognitive Neuroscience, Maastricht Brain Imaging Centre (MBIC) Faculty of Psychology and Neuroscience, Maastricht UniversityMaastrichtNetherlands
| | - Peter De Weerd
- Department of Cognitive Neuroscience, Maastricht Brain Imaging Centre (MBIC) Faculty of Psychology and Neuroscience, Maastricht UniversityMaastrichtNetherlands
- Maastricht Center of Systems Biology (MACSBIO), Faculty of Science & Engineering, Maastricht UniversityMaastrichtNetherlands
| | - Marian Schneider
- Department of Cognitive Neuroscience, Maastricht Brain Imaging Centre (MBIC) Faculty of Psychology and Neuroscience, Maastricht UniversityMaastrichtNetherlands
| | - Omer Faruk Gulban
- Department of Cognitive Neuroscience, Maastricht Brain Imaging Centre (MBIC) Faculty of Psychology and Neuroscience, Maastricht UniversityMaastrichtNetherlands
| | - Dimo Ivanov
- Department of Cognitive Neuroscience, Maastricht Brain Imaging Centre (MBIC) Faculty of Psychology and Neuroscience, Maastricht UniversityMaastrichtNetherlands
| | - Yawen Wang
- Department of Cognitive Neuroscience, Maastricht Brain Imaging Centre (MBIC) Faculty of Psychology and Neuroscience, Maastricht UniversityMaastrichtNetherlands
| | - Kâmil Uludağ
- Center for Neuroscience Imaging Research, Institute for Basic Science and Department of Biomedical Engineering, N Center, Sungkyunkwan UniversityJangan-guRepublic of Korea
- Techna Institute and Koerner Scientist in MR Imaging, University Health NetworkTorontoCanada
| |
Collapse
|
45
|
Layer-dependent functional connectivity methods. Prog Neurobiol 2020; 207:101835. [DOI: 10.1016/j.pneurobio.2020.101835] [Citation(s) in RCA: 39] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2020] [Revised: 03/30/2020] [Accepted: 05/21/2020] [Indexed: 12/16/2022]
|
46
|
Vizioli L, De Martino F, Petro LS, Kersten D, Ugurbil K, Yacoub E, Muckli L. Multivoxel Pattern of Blood Oxygen Level Dependent Activity can be sensitive to stimulus specific fine scale responses. Sci Rep 2020; 10:7565. [PMID: 32371891 PMCID: PMC7200825 DOI: 10.1038/s41598-020-64044-x] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2019] [Accepted: 04/08/2020] [Indexed: 12/25/2022] Open
Abstract
At ultra-high field, fMRI voxels can span the sub-millimeter range, allowing the recording of blood oxygenation level dependent (BOLD) responses at the level of fundamental units of neural computation, such as cortical columns and layers. This sub-millimeter resolution, however, is only nominal in nature as a number of factors limit the spatial acuity of functional voxels. Multivoxel Pattern Analysis (MVPA) may provide a means to detect information at finer spatial scales that may otherwise not be visible at the single voxel level due to limitations in sensitivity and specificity. Here, we evaluate the spatial scale of stimuli specific BOLD responses in multivoxel patterns exploited by linear Support Vector Machine, Linear Discriminant Analysis and Naïve Bayesian classifiers across cortical depths in V1. To this end, we artificially misaligned the testing relative to the training portion of the data in increasing spatial steps, then investigated the breakdown of the classifiers’ performances. A one voxel shift led to a significant decrease in decoding accuracy (p < 0.05) across all cortical depths, indicating that stimulus specific responses in a multivoxel pattern of BOLD activity exploited by multivariate decoders can be as precise as the nominal resolution of single voxels (here 0.8 mm isotropic). Our results further indicate that large draining vessels, prominently residing in proximity of the pial surface, do not, in this case, hinder the ability of MVPA to exploit fine scale patterns of BOLD signals. We argue that tailored analytical approaches can help overcoming limitations in high-resolution fMRI and permit studying the mesoscale organization of the human brain with higher sensitivities.
Collapse
Affiliation(s)
- Luca Vizioli
- CMRR, University of Minnesota, Minneapolis, MN, United States.
| | - Federico De Martino
- CMRR, University of Minnesota, Minneapolis, MN, United States.,Maastricht University, Maastricht, Netherlands
| | | | - Daniel Kersten
- Department of Psychology, University of Minnesota, Minneapolis, MN, United States
| | - Kamil Ugurbil
- CMRR, University of Minnesota, Minneapolis, MN, United States
| | - Essa Yacoub
- CMRR, University of Minnesota, Minneapolis, MN, United States
| | - Lars Muckli
- University of Glasgow, Glasgow, United Kingdom
| |
Collapse
|
47
|
van Dijk JA, Fracasso A, Petridou N, Dumoulin SO. Linear systems analysis for laminar fMRI: Evaluating BOLD amplitude scaling for luminance contrast manipulations. Sci Rep 2020; 10:5462. [PMID: 32214136 PMCID: PMC7096513 DOI: 10.1038/s41598-020-62165-x] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2019] [Accepted: 03/09/2020] [Indexed: 01/18/2023] Open
Abstract
A fundamental assumption of nearly all functional magnetic resonance imaging (fMRI) analyses is that the relationship between local neuronal activity and the blood oxygenation level dependent (BOLD) signal can be described as following linear systems theory. With the advent of ultra-high field (7T and higher) MRI scanners, it has become possible to perform sub-millimeter resolution fMRI in humans. A novel and promising application of sub-millimeter fMRI is measuring responses across cortical depth, i.e. laminar imaging. However, the cortical vasculature and associated directional blood pooling towards the pial surface strongly influence the cortical depth-dependent BOLD signal, particularly for gradient-echo BOLD. This directional pooling may potentially affect BOLD linearity across cortical depth. Here we assess whether the amplitude scaling assumption for linear systems theory holds across cortical depth. For this, we use stimuli with different luminance contrasts to elicit different BOLD response amplitudes. We find that BOLD amplitude across cortical depth scales with luminance contrast, and that this scaling is identical across cortical depth. Although nonlinearities may be present for different stimulus configurations and acquisition protocols, our results suggest that the amplitude scaling assumption for linear systems theory across cortical depth holds for luminance contrast manipulations in sub-millimeter laminar BOLD fMRI.
Collapse
Affiliation(s)
- Jelle A van Dijk
- Experimental Psychology, Utrecht University, Utrecht, NL, Netherlands.
- Spinoza Centre for Neuroimaging, Amsterdam, NL, Netherlands.
| | - Alessio Fracasso
- Spinoza Centre for Neuroimaging, Amsterdam, NL, Netherlands
- Institute of Neuroscience and Psychology, University of Glasgow, Glasgow, G12 8QB, UK
- Radiology Department, Imaging Division, Center for Image Sciences, University Medical Center Utrecht, Utrecht, NL, Netherlands
| | - Natalia Petridou
- Radiology Department, Imaging Division, Center for Image Sciences, University Medical Center Utrecht, Utrecht, NL, Netherlands
| | - Serge O Dumoulin
- Experimental Psychology, Utrecht University, Utrecht, NL, Netherlands
- Spinoza Centre for Neuroimaging, Amsterdam, NL, Netherlands
- Experimental and Applied Psychology, VU University, Amsterdam, NL, Netherlands
| |
Collapse
|
48
|
Bause J, Polimeni JR, Stelzer J, In MH, Ehses P, Kraemer-Fernandez P, Aghaeifar A, Lacosse E, Pohmann R, Scheffler K. Impact of prospective motion correction, distortion correction methods and large vein bias on the spatial accuracy of cortical laminar fMRI at 9.4 Tesla. Neuroimage 2020; 208:116434. [DOI: 10.1016/j.neuroimage.2019.116434] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2019] [Revised: 11/08/2019] [Accepted: 12/02/2019] [Indexed: 01/24/2023] Open
|
49
|
Ultra-high-resolution fMRI of Human Ventral Temporal Cortex Reveals Differential Representation of Categories and Domains. J Neurosci 2020; 40:3008-3024. [PMID: 32094202 DOI: 10.1523/jneurosci.2106-19.2020] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2019] [Revised: 02/10/2020] [Accepted: 02/12/2020] [Indexed: 11/21/2022] Open
Abstract
Human ventral temporal cortex (VTC) is critical for visual recognition. It is thought that this ability is supported by large-scale patterns of activity across VTC that contain information about visual categories. However, it is unknown how category representations in VTC are organized at the submillimeter scale and across cortical depths. To fill this gap in knowledge, we measured BOLD responses in medial and lateral VTC to images spanning 10 categories from five domains (written characters, bodies, faces, places, and objects) at an ultra-high spatial resolution of 0.8 mm using 7 Tesla fMRI in both male and female participants. Representations in lateral VTC were organized most strongly at the general level of domains (e.g., places), whereas medial VTC was also organized at the level of specific categories (e.g., corridors and houses within the domain of places). In both lateral and medial VTC, domain-level and category-level structure decreased with cortical depth, and downsampling our data to standard resolution (2.4 mm) did not reverse differences in representations between lateral and medial VTC. The functional diversity of representations across VTC partitions may allow downstream regions to read out information in a flexible manner according to task demands. These results bridge an important gap between electrophysiological recordings in single neurons at the micron scale in nonhuman primates and standard-resolution fMRI in humans by elucidating distributed responses at the submillimeter scale with ultra-high-resolution fMRI in humans.SIGNIFICANCE STATEMENT Visual recognition is a fundamental ability supported by human ventral temporal cortex (VTC). However, the nature of fine-scale, submillimeter distributed representations in VTC is unknown. Using ultra-high-resolution fMRI of human VTC, we found differential distributed visual representations across lateral and medial VTC. Domain representations (e.g., faces, bodies, places, characters) were most salient in lateral VTC, whereas category representations (e.g., corridors/houses within the domain of places) were equally salient in medial VTC. These results bridge an important gap between electrophysiological recordings in single neurons at a micron scale and fMRI measurements at a millimeter scale.
Collapse
|
50
|
Jorge J, Gretsch F, Najdenovska E, Tuleasca C, Levivier M, Maeder P, Gallichan D, Marques JP, Bach Cuadra M. Improved susceptibility-weighted imaging for high contrast and resolution thalamic nuclei mapping at 7T. Magn Reson Med 2020; 84:1218-1234. [PMID: 32052486 DOI: 10.1002/mrm.28197] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2019] [Revised: 01/10/2020] [Accepted: 01/13/2020] [Indexed: 02/06/2023]
Abstract
PURPOSE The thalamus is an important brain structure and neurosurgical target, but its constituting nuclei are challenging to image non-invasively. Recently, susceptibility-weighted imaging (SWI) at ultra-high field has shown promising capabilities for thalamic nuclei mapping. In this work, several methodological improvements were explored to enhance SWI quality and contrast, and specifically its ability for thalamic imaging. METHODS High-resolution SWI was performed at 7T in healthy participants, and the following techniques were applied: (a) monitoring and retrospective correction of head motion and B0 perturbations using integrated MR navigators, (b) segmentation and removal of venous vessels on the SWI data using vessel enhancement filtering, and (c) contrast enhancement by tuning the parameters of the SWI phase-magnitude combination. The resulting improvements were evaluated with quantitative metrics of image quality, and by comparison to anatomo-histological thalamic atlases. RESULTS Even with sub-millimeter motion and natural breathing, motion and field correction produced clear improvements in both magnitude and phase data quality (76% and 41%, respectively). The improvements were stronger in cases of larger motion/field deviations, mitigating the dependence of image quality on subject performance. Optimizing the SWI phase-magnitude combination yielded substantial improvements in image contrast, particularly in the thalamus, well beyond previously reported SWI results. The atlas comparisons provided compelling evidence of anatomical correspondence between SWI features and several thalamic nuclei, for example, the ventral intermediate nucleus. Vein detection performed favorably inside the thalamus, and vein removal further improved visualization. CONCLUSION Altogether, the proposed developments substantially improve high-resolution SWI, particularly for thalamic nuclei imaging.
Collapse
Affiliation(s)
- João Jorge
- Medical Image Analysis Laboratory, Center for Biomedical Imaging (CIBM), University of Lausanne, Lausanne, Switzerland.,Laboratory for Functional and Metabolic Imaging, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Frédéric Gretsch
- Laboratory for Functional and Metabolic Imaging, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Elena Najdenovska
- Medical Image Analysis Laboratory, Center for Biomedical Imaging (CIBM), University of Lausanne, Lausanne, Switzerland.,Department of Radiology, Centre Hospitalier Universitaire Vaudois, University of Lausanne, Lausanne, Switzerland
| | - Constantin Tuleasca
- Department of Clinical Neurosciences, Neurosurgery Service and Gamma Knife Center, Centre Hospitalier Universitaire Vaudois, Lausanne, Switzerland.,Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland.,Faculty of Biology and Medicine (FBM), University of Lausanne (UNIL), Lausanne, Switzerland
| | - Marc Levivier
- Department of Clinical Neurosciences, Neurosurgery Service and Gamma Knife Center, Centre Hospitalier Universitaire Vaudois, Lausanne, Switzerland.,Faculty of Biology and Medicine (FBM), University of Lausanne (UNIL), Lausanne, Switzerland
| | - Philippe Maeder
- Department of Radiology, Centre Hospitalier Universitaire Vaudois, University of Lausanne, Lausanne, Switzerland
| | - Daniel Gallichan
- Cardiff University Brain Research Imaging Centre, School of Engineering, Cardiff University, Cardiff, UK
| | - José P Marques
- Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, The Netherlands
| | - Meritxell Bach Cuadra
- Medical Image Analysis Laboratory, Center for Biomedical Imaging (CIBM), University of Lausanne, Lausanne, Switzerland.,Department of Radiology, Centre Hospitalier Universitaire Vaudois, University of Lausanne, Lausanne, Switzerland.,Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| |
Collapse
|