1
|
Polimeni JR, Lewis LD. Imaging faster neural dynamics with fast fMRI: A need for updated models of the hemodynamic response. Prog Neurobiol 2021; 207:102174. [PMID: 34525404 PMCID: PMC8688322 DOI: 10.1016/j.pneurobio.2021.102174] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2020] [Revised: 07/30/2021] [Accepted: 09/08/2021] [Indexed: 12/20/2022]
Abstract
Fast fMRI enables the detection of neural dynamics over timescales of hundreds of milliseconds, suggesting it may provide a new avenue for studying subsecond neural processes in the human brain. The magnitudes of these fast fMRI dynamics are far greater than predicted by canonical models of the hemodynamic response. Several studies have established nonlinear properties of the hemodynamic response that have significant implications for fast fMRI. We first review nonlinear properties of the hemodynamic response function that may underlie fast fMRI signals. We then illustrate the breakdown of canonical hemodynamic response models in the context of fast neural dynamics. We will then argue that the canonical hemodynamic response function is not likely to reflect the BOLD response to neuronal activity driven by sparse or naturalistic stimuli or perhaps to spontaneous neuronal fluctuations in the resting state. These properties suggest that fast fMRI is capable of tracking surprisingly fast neuronal dynamics, and we discuss the neuroscientific questions that could be addressed using this approach.
Collapse
Affiliation(s)
- Jonathan R Polimeni
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, USA; Department of Radiology, Harvard Medical School, Boston, MA, USA; Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, USA.
| | - Laura D Lewis
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, USA; Department of Biomedical Engineering, Boston University, Boston, MA, USA.
| |
Collapse
|
2
|
Wang L, Li C, Chen D, Lv X, Go R, Wu J, Yan T. Hemodynamic response varies across tactile stimuli with different temporal structures. Hum Brain Mapp 2020; 42:587-597. [PMID: 33169898 PMCID: PMC7814760 DOI: 10.1002/hbm.25243] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2020] [Revised: 10/01/2020] [Accepted: 10/03/2020] [Indexed: 11/23/2022] Open
Abstract
Tactile stimuli can be distinguished based on their temporal features (e.g., duration, local frequency, and number of pulses), which are fundamental for vibrotactile frequency perception. Characterizing how the hemodynamic response changes in shape across experimental conditions is important for designing and interpreting fMRI studies on tactile information processing. In this study, we focused on periodic tactile stimuli with different temporal structures and explored the hemodynamic response function (HRF) induced by these stimuli. We found that HRFs were stimulus‐dependent in tactile‐related brain areas. Continuous stimuli induced a greater area of activation and a stronger and narrower hemodynamic response than intermittent stimuli with the same duration. The magnitude of the HRF increased with increasing stimulus duration. By normalizing the characteristics into topographic matrix, nonlinearity was obvious. These results suggested that stimulation patterns and duration within a cycle may be key characters for distinguishing different stimuli. We conclude that different temporal structures of tactile stimuli induced different HRFs, which are essential for vibrotactile perception and should be considered in fMRI experimental designs and analyses.
Collapse
Affiliation(s)
- Luyao Wang
- School of Mechatronical Engineering, Beijing Institute of Technology, Beijing, China.,Beijing Advanced Innovation Center for Intelligent Robots and Systems, Beijing Institute of Technology, Beijing, China
| | - Chunlin Li
- School of Biomedical Engineering, Capital Medical University, Beijing, China
| | - Duanduan Chen
- School of Life Science, Beijing Institute of Technology, Beijing, China
| | - Xiaoyu Lv
- School of Mechatronical Engineering, Beijing Institute of Technology, Beijing, China.,Beijing Advanced Innovation Center for Intelligent Robots and Systems, Beijing Institute of Technology, Beijing, China
| | - Ritsu Go
- School of Mechatronical Engineering, Beijing Institute of Technology, Beijing, China.,Beijing Advanced Innovation Center for Intelligent Robots and Systems, Beijing Institute of Technology, Beijing, China
| | - Jinglong Wu
- School of Mechatronical Engineering, Beijing Institute of Technology, Beijing, China.,Beijing Advanced Innovation Center for Intelligent Robots and Systems, Beijing Institute of Technology, Beijing, China.,Graduate School of Interdisciplinary Science and Engineering in Health Systems, Okayama University, Okayama, Japan
| | - Tianyi Yan
- School of Life Science, Beijing Institute of Technology, Beijing, China.,Beijing Advanced Innovation Center for Intelligent Robots and Systems, Beijing Institute of Technology, Beijing, China
| |
Collapse
|
3
|
Aso T, Urayama SI, Fukuyama H, Le Bihan D. Comparison of diffusion-weighted fMRI and BOLD fMRI responses in a verbal working memory task. Neuroimage 2013; 67:25-32. [DOI: 10.1016/j.neuroimage.2012.11.005] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2012] [Revised: 10/04/2012] [Accepted: 11/02/2012] [Indexed: 11/25/2022] Open
|
4
|
Marxen M, Cassidy RJ, Dawson TL, Ross B, Graham SJ. Transient and sustained components of the sensorimotor BOLD response in fMRI. Magn Reson Imaging 2012; 30:837-47. [PMID: 22495237 DOI: 10.1016/j.mri.2012.02.007] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2011] [Revised: 11/18/2011] [Accepted: 02/14/2012] [Indexed: 11/15/2022]
Abstract
Blood oxygenation level-dependent (BOLD) signal time courses in functional magnetic resonance imaging are estimated within the framework of general linear modeling by convolving an input function, that represents neural activity, with a canonical hemodynamic response function (HRF). Here we investigate the performance of different neural input functions and latency-optimized HRFs for modeling BOLD signals in response to vibrotactile somatosensory stimuli of variable durations (0.5, 1, 4, 7 s) in 14 young, healthy adults who were required to make button press responses at each stimulus cessation. Informed by electrophysiology and the behavioral task, three nested models with an increasing number of parameters were considered: a boxcar; boxcar and offset transient; and onset transient, boxcar and offset transient (TBT). The TBT model provided the best fit of the group-averaged BOLD time courses based on χ(2) and F statistics. Only the TBT model was capable of fitting the bimodal shape of the BOLD response to the 7-s stimulus and the relative peak amplitudes for all stimulus lengths in key somatosensory and motor areas. This suggests that the TBT model provides a more comprehensive description of brain sensorimotor responses in this experiment than provided by the simple boxcar model. Work comparing the activation maps obtained with the TBT model with magnetoencephalography data is under way.
Collapse
Affiliation(s)
- Michael Marxen
- Department of Psychiatry and Psychotherapy, Technische Universität Dresden, Würzburger Straße 35, Dresden, Germany.
| | | | | | | | | |
Collapse
|
5
|
Ramkumar P, Parkkonen L, Hari R. Oscillatory response function: towards a parametric model of rhythmic brain activity. Hum Brain Mapp 2010; 31:820-34. [PMID: 19957269 DOI: 10.1002/hbm.20907] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
Rhythmic brain activity, measured by magnetoencephalography (MEG), is modulated during stimulation and task performance. Here, we introduce an oscillatory response function (ORF) to predict the dynamic suppression-rebound modulation of brain rhythms during a stimulus sequence. We derived a class of parametric models for the ORF in a generalized convolution framework. The model parameters were estimated from MEG data acquired from 10 subjects during bilateral tactile stimulation of fingers (stimulus rates of 4 Hz and 10 Hz in blocks of 0.5, 1, 2, and 4 s). The envelopes of the 17-23 Hz rhythmic activity, computed for sensors above the rolandic region, correlated 25%-43% better with the envelopes predicted by the models than by the stimulus time course (boxcar). A linear model with separate convolution kernels for onset and offset responses gave the best prediction. We studied the generalizability of this model with data from 5 different subjects during a separate bilateral tactile sequence by first identifying neural sources of the 17-23 Hz activity using cortically constrained minimum norm estimates. Both the model and the boxcar predicted strongest modulation in the primary motor cortex. For short-duration stimulus blocks, the model predicted the envelope of the cortical currents 20% better than the boxcar did. These results suggest that ORFs could concisely describe brain rhythms during different stimuli, tasks, and pathologies.
Collapse
Affiliation(s)
- Pavan Ramkumar
- Brain Research Unit, Low Temperature Laboratory, Helsinki University of Technology, 02015 TKK, Finland.
| | | | | |
Collapse
|
6
|
Affiliation(s)
- Riitta Hari
- Brain Research Unit, Low Temperature Laboratory, Aalto University School of Science and Technology, AALTO, Espoo, Finland.
| | | | | |
Collapse
|
7
|
Nangini C, Hlushchuk Y, Hari R. Predicting stimulus-rate sensitivity of human somatosensory fMRI signals with MEG. Hum Brain Mapp 2009; 30:1824-32. [PMID: 19378275 DOI: 10.1002/hbm.20787] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
Abstract
With increasing stimulus rate (SR), cortical EEG and MEG responses typically decrease in amplitude whereas BOLD fMRI signals increase. To address this discrepancy, we predicted BOLD responses with squared MEG waveforms using a recently proposed energy-density model. Tactile stimuli were delivered to finger tips at SRs of 1, 4, or 10 Hz in successive 25-s blocks, and brain signals were detected from area 3b of the primary somatosensory cortex of nine healthy adults using a 306-channel whole-scalp neuromagnetometer and a 3-T fMRI magnet. The main MEG deflections decreased in amplitude as a function of SR, whereas the BOLD signals increased from 1- to 4-Hz SR, with no further change at 10 Hz. MEG energy densities, obtained over the whole stimulus train and convolved with different hemodynamic response functions, predicted both the shape and amplitude of the BOLD signals well, and incorporation of nonlinear terms into the model did not offer any further advantage. Thus, squared MEG waveforms obtained over the entire stimulus train provided an appropriate estimate of area 3b neuronal activity associated with the BOLD signal.
Collapse
Affiliation(s)
- Cathy Nangini
- Brain Research Unit, Low Temperature Laboratory, Helsinki University of Technology, TKK, Espoo, Finland.
| | | | | |
Collapse
|
8
|
Nangini C, Tam F, Graham SJ. A novel method for integrating MEG and BOLD fMRI signals with the linear convolution model in human primary somatosensory cortex. Hum Brain Mapp 2008; 29:97-106. [PMID: 17290370 PMCID: PMC4896808 DOI: 10.1002/hbm.20361] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
Characterizing the neurovascular coupling between hemodynamic signals and their neural origins is crucial to functional neuroimaging research, even more so as new methods become available for integrating results from different functional neuroimaging modalities. We present a novel method to relate magnetoencephalography (MEG) and BOLD fMRI data from primary somatosensory cortex within the context of the linear convolution model. This model, which relates neural activity to BOLD signal change, has been widely used to predict BOLD signals but typically lacks experimentally derived measurements of neural activity. In this study, an fMRI experiment is performed using variable-duration (< or =1 s) vibrotactile stimuli applied at 22 Hz, analogous to a previously published MEG study (Nangini et al., [2006]: Neuroimage 33:252-262), testing whether MEG source waveforms from the previous study can inform the convolution model and improve BOLD signal estimates across all stimulus durations. The typical formulation of the convolution model in which the input is given by the stimulus profile is referred to as Model 1. Model 2 is based on an energy argument relating metabolic demand to the postsynaptic currents largely responsible for the MEG current dipoles, and uses the energy density of the estimated MEG source waveforms as input to the convolution model. It is shown that Model 2 improves the BOLD signal estimates compared to Model 1 under the experimental conditions implemented, suggesting that MEG energy density can be a useful index of hemodynamic activity.
Collapse
Affiliation(s)
- Cathy Nangini
- Department of Medical Biophysics, University of Toronto, Ontario, Canada.
| | | | | |
Collapse
|