Resting-state fMRI signals contain spectral signatures of local hemodynamic response timing

  1. Sydney M Bailes
  2. Daniel EP Gomez
  3. Beverly Setzer
  4. Laura D Lewis  Is a corresponding author
  1. Boston University, United States
  2. Massachusetts Institute of Technology, United States

Abstract

Functional magnetic resonance imaging (fMRI) has proven to be a powerful tool for noninvasively measuring human brain activity; yet, thus far, fMRI has been relatively limited in its temporal resolution. A key challenge is understanding the relationship between neural activity and the blood-oxygenation-level-dependent (BOLD) signal obtained from fMRI, generally modeled by the hemodynamic response function (HRF). The timing of the HRF varies across the brain and individuals, confounding our ability to make inferences about the timing of the underlying neural processes. Here we show that resting-state fMRI signals contain information about HRF temporal dynamics that can be leveraged to understand and characterize variations in HRF timing across both cortical and subcortical regions. We found that the frequency spectrum of resting-state fMRI signals significantly differs between voxels with fast versus slow HRFs in human visual cortex. These spectral differences extended to subcortex as well, revealing significantly faster hemodynamic timing in the lateral geniculate nucleus of the thalamus. Ultimately, our results demonstrate that the temporal properties of the HRF impact the spectral content of resting-state fMRI signals and enable voxel-wise characterization of relative hemodynamic response timing. Furthermore, our results show that caution should be used in studies of resting-state fMRI spectral properties, because differences in fMRI frequency content can arise from purely vascular origins. This finding provides new insight into the temporal properties of fMRI signals across voxels, which is crucial for accurate fMRI analyses, and enhances the ability of fast fMRI to identify and track fast neural dynamics.

Data availability

The data used in this paper has been deposited on OpenNeuro (https://doi.org/10.18112/openneuro.ds004645.v1.0.0).

The following data sets were generated

Article and author information

Author details

  1. Sydney M Bailes

    Department of Biomedical Engineering, Boston University, Boston, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1712-9233
  2. Daniel EP Gomez

    Department of Biomedical Engineering, Boston University, Boston, United States
    Competing interests
    The authors declare that no competing interests exist.
  3. Beverly Setzer

    Department of Biomedical Engineering, Boston University, Boston, United States
    Competing interests
    The authors declare that no competing interests exist.
  4. Laura D Lewis

    Institute for Medical Engineering and Science, Massachusetts Institute of Technology, Cambridge, United States
    For correspondence
    ldlewis@mit.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-4003-0277

Funding

National Institutes of Health (R00-MH111748)

  • Laura D Lewis

National Institutes of Health (U19-NS123717)

  • Laura D Lewis

National Institutes of Health (R01-AG070135)

  • Laura D Lewis

Searle Scholars Program

  • Laura D Lewis

Pew Charitable Trusts (Pew Biomedical Scholars Program)

  • Laura D Lewis

Sloan School of Management, Massachusetts Institute of Technology (Sloan Fellowship)

  • Laura D Lewis

One Mind (One Mind Rising Star Award)

  • Laura D Lewis

National Institutes of Health (T32-GM008764)

  • Sydney M Bailes

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Ethics

Human subjects: All experimental procedures were approved by the Massachusetts General Hospital Institutional Review Board and all subjects provided informed consent. (protocol number: 2014P001068),

Copyright

© 2023, Bailes et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,263
    views
  • 157
    downloads
  • 6
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Sydney M Bailes
  2. Daniel EP Gomez
  3. Beverly Setzer
  4. Laura D Lewis
(2023)
Resting-state fMRI signals contain spectral signatures of local hemodynamic response timing
eLife 12:e86453.
https://doi.org/10.7554/eLife.86453

Share this article

https://doi.org/10.7554/eLife.86453

Further reading

    1. Neuroscience
    Yiting Li, Wenqu Yin ... Baoming Li
    Research Article

    Time estimation is an essential prerequisite underlying various cognitive functions. Previous studies identified ‘sequential firing’ and ‘activity ramps’ as the primary neuron activity patterns in the medial frontal cortex (mPFC) that could convey information regarding time. However, the relationship between these patterns and the timing behavior has not been fully understood. In this study, we utilized in vivo calcium imaging of mPFC in rats performing a timing task. We observed cells that showed selective activation at trial start, end, or during the timing interval. By aligning long-term time-lapse datasets, we discovered that sequential patterns of time coding were stable over weeks, while cells coding for trial start or end showed constant dynamism. Furthermore, with a novel behavior design that allowed the animal to determine individual trial interval, we were able to demonstrate that real-time adjustment in the sequence procession speed closely tracked the trial-to-trial interval variations. And errors in the rats’ timing behavior can be primarily attributed to the premature ending of the time sequence. Together, our data suggest that sequential activity maybe a stable neural substrate that represents time under physiological conditions. Furthermore, our results imply the existence of a unique cell type in the mPFC that participates in the time-related sequences. Future characterization of this cell type could provide important insights in the neural mechanism of timing and related cognitive functions.

    1. Neuroscience
    Rossella Conti, Céline Auger
    Research Article

    Granule cells of the cerebellum make up to 175,000 excitatory synapses on a single Purkinje cell, encoding the wide variety of information from the mossy fibre inputs into the cerebellar cortex. The granule cell axon is made of an ascending portion and a long parallel fibre extending at right angles, an architecture suggesting that synapses formed by the two segments of the axon could encode different information. There are controversial indications that ascending axon (AA) and parallel fibre (PF) synapse properties and modalities of plasticity are different. We tested the hypothesis that AA and PF synapses encode different information, and that the association of these distinct inputs to Purkinje cells might be relevant to the circuit and trigger plasticity, similar to the coincident activation of PF and climbing fibre inputs. Here, by recording synaptic currents in Purkinje cells from either proximal or distal granule cells (mostly AA and PF synapses, respectively), we describe a new form of associative plasticity between these two distinct granule cell inputs. We show for the first time that synchronous AA and PF repetitive train stimulation, with inhibition intact, triggers long-term potentiation (LTP) at AA synapses specifically. Furthermore, the timing of the presentation of the two inputs controls the outcome of plasticity and induction requires NMDAR and mGluR1 activation. The long length of the PFs allows us to preferentially activate the two inputs independently, and despite a lack of morphological reconstruction of the connections, these observations reinforce the suggestion that AA and PF synapses have different coding capabilities and plasticity that is associative, enabling effective association of information transmitted via granule cells.