Resting-state fMRI Signal Research Articles

Sex-specific age-related differences in cerebrospinal fluid clearance assessed by resting-state functional magnetic resonance imaging

Cerebrospinal fluid (CSF) flow may assist the clearance of brain wastes, such as amyloid-β (Aβ) and tau, and thus play an important role in aging and dementias. However, a lack of non-invasive tools to assess the CSF dynamics-related clearance in humans hindered the understanding of the relevant changes in healthy aging. The global infra-slow (<0.1 Hz) brain activity measured by the global mean resting-state fMRI signal (gBOLD) was recently found to be coupled by large CSF movements. This coupling has been found to correlate with various pathologies of Alzheimer's disease (AD), particularly Aβ pathology, linking it to waste clearance. Using resting-state fMRI data from a group of 719 healthy aging participants, we examined the sex-specific differences of the gBOLD-CSF coupling over a wide age range between 36-100 years of age. We found that this coupling index remains stable before around age 55 and then starts to decline afterward, particularly in females. Menopause may contribute to the accelerated decline in females.

Disparity in temporal and spatial relationships between resting-state electrophysiological and fMRI signals

Resting-state brain networks (RSNs) have been widely applied in health and disease, but the interpretation of RSNs in terms of the underlying neural activity is unclear. To address this fundamental question, we conducted simultaneous recordings of whole-brain resting-state functional magnetic resonance imaging (rsfMRI) and electrophysiology signals in two separate brain regions of rats. Our data reveal that for both recording sites, spatial maps derived from band-specific local field potential (LFP) power can account for up to 90% of the spatial variability in RSNs derived from rsfMRI signals. Surprisingly, the time series of LFP band power can only explain to a maximum of 35% of the temporal variance of the local rsfMRI time course from the same site. In addition, regressing out time series of LFP power from rsfMRI signals has minimal impact on the spatial patterns of rsfMRI-based RSNs. This disparity in the spatial and temporal relationships between resting-state electrophysiology and rsfMRI signals suggests that electrophysiological activity alone does not fully explain the effects observed in the rsfMRI signal, implying the existence of an rsfMRI component contributed by ‘electrophysiology-invisible’ signals. These findings offer a novel perspective on our understanding of RSN interpretation.

Disparity in temporal and spatial relationships between resting-state electrophysiological and fMRI signals.

Resting-state brain networks (RSNs) have been widely applied in health and disease, but the interpretation of RSNs in terms of the underlying neural activity is unclear. To address this fundamental question, we conducted simultaneous recordings of whole-brain resting-state functional magnetic resonance imaging (rsfMRI) and electrophysiology signals in two separate brain regions of rats. Our data reveal that for both recording sites, spatial maps derived from band-specific local field potential (LFP) power can account for up to 90% of the spatial variability in RSNs derived from rsfMRI signals. Surprisingly, the time series of LFP band power can only explain to a maximum of 35% of the temporal variance of the local rsfMRI time course from the same site. In addition, regressing out time series of LFP power from rsfMRI signals has minimal impact on the spatial patterns of rsfMRI-based RSNs. This disparity in the spatial and temporal relationships between resting-state electrophysiology and rsfMRI signals suggests that electrophysiological activity alone does not fully explain the effects observed in the rsfMRI signal, implying the existence of an rsfMRI component contributed by 'electrophysiology-invisible' signals. These findings offer a novel perspective on our understanding of RSN interpretation.

Widespread Autonomic Physiological Coupling Across the Brain-Body Axis.

The brain is closely attuned to visceral signals from the body's internal environment, as evidenced by the numerous associations between neural, hemodynamic, and peripheral physiological signals. We show that these brain-body co-fluctuations can be captured by a single spatiotemporal pattern. Across several independent samples, as well as single-echo and multi-echo fMRI data acquisition sequences, we identify widespread co-fluctuations in the low-frequency range (0.01 - 0.1 Hz) between resting-state global fMRI signals, neural activity, and a host of autonomic signals spanning cardiovascular, pulmonary, exocrine and smooth muscle systems. The same brain-body co-fluctuations observed at rest are elicited by arousal induced by cued deep breathing and intermittent sensory stimuli, as well as spontaneous phasic EEG events during sleep. Further, we show that the spatial structure of global fMRI signals is maintained under experimental suppression of end-tidal carbon dioxide (PETCO2) variations, suggesting that respiratory-driven fluctuations in arterial CO2 accompanying arousal cannot explain the origin of these signals in the brain. These findings establish the global fMRI signal as a significant component of the arousal response governed by the autonomic nervous system.

Disparity in temporal and spatial relationships between resting-state electrophysiological and fMRI signals.

Resting-state brain networks (RSNs) have been widely applied in health and disease, but the interpretation of RSNs in terms of the underlying neural activity is unclear. To address this fundamental question, we conducted simultaneous recordings of whole-brain resting-state functional magnetic resonance imaging (rsfMRI) and electrophysiology signals in two separate brain regions of rats. Our data reveal that for both recording sites, spatial maps derived from band-specific local field potential (LFP) power can account for up to 90% of the spatial variability in RSNs derived from rsfMRI signals. Surprisingly, the time series of LFP band power can only explain to a maximum of 35% of the temporal variance of the local rsfMRI time course from the same site. In addition, regressing out time series of LFP power from rsfMRI signals has minimal impact on the spatial patterns of rsfMRI-based RSNs. This disparity in the spatial and temporal relationships between resting-state electrophysiology and rsfMRI signals suggests that electrophysiological activity alone does not fully explain the effects observed in the rsfMRI signal, implying the existence of an rsfMRI component contributed by "electrophysiology-invisible" signals. These findings offer a novel perspective on our understanding of RSN interpretation.

Disparity in temporal and spatial relationships between resting-state electrophysiological and fMRI signals.

Resting-state brain networks (RSNs) have been widely applied in health and disease, but the interpretation of RSNs in terms of the underlying neural activity is unclear. To address this fundamental question, we conducted simultaneous recordings of whole-brain resting-state functional magnetic resonance imaging (rsfMRI) and electrophysiology signals in two separate brain regions of rats. Our data reveal that for both recording sites, spatial maps derived from band-specific local field potential (LFP) power can account for up to 90% of the spatial variability in RSNs derived from rsfMRI signals. Surprisingly, the time series of LFP band power can only explain to a maximum of 35% of the temporal variance of the local rsfMRI time course from the same site. In addition, regressing out time series of LFP power from rsfMRI signals has minimal impact on the spatial patterns of rsfMRI-based RSNs. This disparity in the spatial and temporal relationships between resting-state electrophysiology and rsfMRI signals suggests that electrophysiological activity alone does not fully explain the effects observed in the rsfMRI signal, implying the existence of an rsfMRI component contributed by "electrophysiology-invisible" signals. These findings offer a novel perspective on our understanding of RSN interpretation.

Gradient synchronization for multivariate functional data, with application to brain connectivity.

Quantifying the association between components of multivariate random curves is of general interest and is a ubiquitous and basic problem that can be addressed with functional data analysis. An important application is the problem of assessing functional connectivity based on functional magnetic resonance imaging (fMRI), where one aims to determine the similarity of fMRI time courses that are recorded on anatomically separated brain regions. In the functional brain connectivity literature, the static temporal Pearson correlation has been the prevailing measure for functional connectivity. However, recent research has revealed temporally changing patterns of functional connectivity, leading to the study of dynamic functional connectivity. This motivates new similarity measures for pairs of random curves that reflect the dynamic features of functional similarity. Specifically, we introduce gradient synchronization measures in a general setting. These similarity measures are based on the concordance and discordance of the gradients between paired smooth random functions. Asymptotic normality of the proposed estimates is obtained under regularity conditions. We illustrate the proposed synchronization measures via simulations and an application to resting-state fMRI signals from the Alzheimer's Disease Neuroimaging Initiative and they are found to improve discrimination between subjects with different disease status.

Uncovering Hemispheric Asymmetry and Directed Oscillatory Brain-Heart Interplay in Anxiety Processing: An fMRI Study.

Brain-heart interactions (BHI) are critical for generating and processing emotions, including anxiety. Understanding specific neural correlates would be instrumental for greater comprehension and potential therapeutic interventions of anxiety disorders. While prior work has implicated the pontine structure as a central processor in cardiac regulation in anxiety, the distributed nature of anxiety processing across the cortex remains elusive. To address this, we performed a whole-brain-heart analysis using the full frequency directed transfer function to study resting-state spectral differences in BHI between high and low anxiety groups undergoing fMRI scans. Our findings revealed a hemispheric asymmetry in low-frequency interplay (0.05 Hz - 0.15 Hz) characterized by ascending BHI to the left insula and descending BHI from the right insula. Furthermore, we provide evidence supporting the "pacemaker hypothesis", highlighting the pons' function in regulating cardiac activity. Higher frequency interplay (0.2 Hz - 0.4Hz) demonstrate a preference for ascending interactions, particularly towards ventral prefrontal cortical activity in high anxiety groups, suggesting the heart's role in triggering a cognitive response to regulate anxiety. These findings highlight the impact of anxiety on BHI, contributing to a better understanding of its effect on the resting-state fMRI signal, with further implications for potential therapeutic interventions in treating anxiety disorders.

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

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.

Connectivity gradients in spontaneous brain activity at multiple frequency bands.

The intrinsic organizational structure of the brain is reflected in spontaneous brain oscillations. Its functional integration and segregation hierarchy have been discovered in space by leveraging gradient approaches to low-frequency functional connectivity. This hierarchy of brain oscillations has not yet been fully understood, since previous studies have mainly concentrated on the brain oscillations from a single limited frequency range (~ 0.01-0.1Hz). In this work, we extended the frequency range and performed gradient analysis across multiple frequency bands of fast resting-state fMRI signals from the Human Connectome Project and condensed a frequency-rank cortical map of the highest gradient. We found that the coarse skeletons of the functional organization hierarchy are generalizable across the multiple frequency bands. Beyond that, the highest integration levels of connectivity vary in the frequency domain across different large-scale brain networks. These findings are replicated in another independent dataset and demonstrated that different brain networks can integrate information at varying rates, indicating the significance of examining the intrinsic architecture of spontaneous brain activity from the perspective of multiple frequency bands.

Critical scaling of whole-brain resting-state dynamics

Scale invariance is a characteristic of neural activity. How this property emerges from neural interactions remains a fundamental question. Here, we studied the relation between scale-invariant brain dynamics and structural connectivity by analyzing human resting-state (rs-) fMRI signals, together with diffusion MRI (dMRI) connectivity and its approximation as an exponentially decaying function of the distance between brain regions. We analyzed the rs-fMRI dynamics using functional connectivity and a recently proposed phenomenological renormalization group (PRG) method that tracks the change of collective activity after successive coarse-graining at different scales. We found that brain dynamics display power-law correlations and power-law scaling as a function of PRG coarse-graining based on functional or structural connectivity. Moreover, we modeled the brain activity using a network of spins interacting through large-scale connectivity and presenting a phase transition between ordered and disordered phases. Within this simple model, we found that the observed scaling features were likely to emerge from critical dynamics and connections exponentially decaying with distance. In conclusion, our study tests the PRG method using large-scale brain activity and theoretical models and suggests that scaling of rs-fMRI activity relates to criticality.

Generating dynamic carbon-dioxide traces from respiration-belt recordings: Feasibility using neural networks and application in functional magnetic resonance imaging.

In the context of functional magnetic resonance imaging (fMRI), carbon dioxide (CO2) is a well-known vasodilator that has been widely used to monitor and interrogate vascular physiology. Moreover, spontaneous fluctuations in end-tidal carbon dioxide (PETCO2) reflects changes in arterial CO2 and has been demonstrated as the largest physiological noise source for denoising the low-frequency range of the resting-state fMRI (rs-fMRI) signal. However, the majority of rs-fMRI studies do not involve CO2 recordings, and most often only heart rate and respiration are recorded. While the intrinsic link between these latter metrics and CO2 led to suggested possible analytical models, they have not been widely applied. In this proof-of-concept study, we propose a deep-learning (DL) approach to reconstruct CO2 and PETCO2 data from respiration waveforms in the resting state. We demonstrate that the one-to-one mapping between respiration and CO2 recordings can be well predicted using fully convolutional networks (FCNs), achieving a Pearson correlation coefficient (r) of 0.946 ± 0.056 with the ground truth CO2. Moreover, dynamic PETCO2 can be successfully derived from the predicted CO2, achieving r of 0.512 ± 0.269 with the ground truth. Importantly, the FCN-based methods outperform previously proposed analytical methods. In addition, we provide guidelines for quality assurance of respiration recordings for the purposes of CO2 prediction. Our results demonstrate that dynamic CO2 can be obtained from respiration-volume using neural networks, complementing the still few reports in DL of physiological fMRI signals, and paving the way for further research in DL based bio-signal processing.

Spontaneous Activity Patterns in Human Attention Networks Code for Hand Movements.

Recent evidence suggests that, in the absence of any task, spontaneous brain activity patterns and connectivity in the visual and motor cortex code for natural stimuli and actions, respectively. These "resting-state" activity patterns may underlie the maintenance and consolidation (replay) of information states coding for ecological stimuli and behaviors. In this study, we examine whether replay patterns occur in resting-state activity in association cortex grouped into high-order cognitive networks not directly processing sensory inputs or motor outputs. Fifteen participants (7 females) performed four hand movements during an fMRI study. Three movements were ecological. The fourth movement as control was less ecological. Before and after the task scans, we acquired resting-state fMRI scans. The analysis examined whether multivertex task activation patterns for the four movements computed at the cortical surface in different brain networks resembled spontaneous activity patterns measured at rest. For each movement, we computed a vector of r values indicating the strength of the similarity between the mean task activation pattern and frame-by-frame resting-state patterns. We computed a cumulative distribution function of r 2 values and used the 90th percentile cutoff value for comparison. In the dorsal attention network, resting-state patterns were more likely to match task patterns for the ecological movements than the control movement. In contrast, rest-task pattern correlation was more likely for less ecological movement in the ventral attention network. These findings show that spontaneous activity patterns in human attention networks code for hand movements.SIGNIFICANCE STATEMENT fMRI indirectly measures neural activity noninvasively. Resting-state (spontaneous) fMRI signals measured in the absence of any task resemble signals evoked by task performance both in topography and inter-regional (functional) connectivity. However, the function of spontaneous brain activity is unknown. We recently showed that spatial activity patterns evoked by visual and motor tasks in visual and motor cortex, respectively, occur at rest in the absence of any stimulus or response. Here we show that activity patterns related to hand movements replay at rest in frontoparietal regions of the human attention system. These findings show that spontaneous activity in the human cortex may mediate the maintenance and consolidation of information states coding for ecological stimuli and behaviors.

Connectome-based predictive modeling for functional recovery of acute ischemic stroke

Patients of acute ischemic stroke possess considerable chance of recovery of various levels in the first several weeks after stroke onset. Prognosis of functional recovery is important for decision-making in poststroke patient care and placement. Poststroke functional recovery has conventionally been based on demographic and clinical variables such as age, gender, and severity of stroke impairment. On the other hand, the concept of connectome has become a basis of interpreting the functional impairment and recovery of stroke patients. In this research, the connectome-based predictive modeling was used to provide predictive models for prognosing poststroke functional recovery. Predictive models were developed to use the brain connectivity at stroke onset to predict functional assessment scores at one or three months later, or to use the brain connectivity one-month poststroke to predict functional assessment scores at three months after stroke onset. The brain connectivity was computed from the resting-state fMRI signals. The functional assessment scores used in this research included modified Rankin Scale (mRS) and Barthel Index (BI). This research found significant models that used the brain connectivity at onset to predict the mRS one-month poststroke and to predict the BI three-month poststroke for patients with supratentorial infarction, as well as predictive models that used the brain connectivity one-month poststroke to predict the mRS three-month poststroke for patients with supratentorial infarction in the right hemisphere. The connectome-based predictive modeling could provide clinical value in prognosis of acute ischemic stroke.

Automatic diagnosis of late-life depression by 3D convolutional neural networks and cross-sample Entropy analysis from resting-state fMRI

Resting-state fMRI has been widely used in investigating the pathophysiology of late-life depression (LLD). Unlike the conventional linear approach, cross-sample entropy (CSE) analysis shows the nonlinear property in fMRI signals between brain regions. Moreover, recent advances in deep learning, such as convolutional neural networks (CNNs), provide a timely application for understanding LLD. Accurate and prompt diagnosis is essential in LLD; hence, this study aimed to combine CNN and CSE analysis to discriminate LLD patients and non-depressed comparison older adults based on brain resting-state fMRI signals. Seventy-seven older adults, including 49 patients and 28 comparison older adults, were included for fMRI scans. Three-dimensional CSEs with volumes corresponding to 90 seed regions of interest of each participant were developed and fed into models for disease classification and depression severity prediction. We obtained a diagnostic accuracy > 85% in the superior frontal gyrus (left dorsolateral and right orbital parts), left insula, and right middle occipital gyrus. With a mean root-mean-square error (RMSE) of 2.41, three separate models were required to predict depressive symptoms in the severe, moderate, and mild depression groups. The CSE volumes in the left inferior parietal lobule, left parahippocampal gyrus, and left postcentral gyrus performed best in each respective model. Combined complexity analysis and deep learning algorithms can classify patients with LLD from comparison older adults and predict symptom severity based on fMRI data. Such application can be utilized in precision medicine for disease detection and symptom monitoring in LLD.

What lies underneath: Precise classification of brain states using time-dependent topological structure of dynamics.

The self-organising global dynamics underlying brain states emerge from complex recursive nonlinear interactions between interconnected brain regions. Until now, most efforts of capturing the causal mechanistic generating principles have supposed underlying stationarity, being unable to describe the non-stationarity of brain dynamics, i.e. time-dependent changes. Here, we present a novel framework able to characterise brain states with high specificity, precisely by modelling the time-dependent dynamics. Through describing a topological structure associated to the brain state at each moment in time (its attractor or ‘information structure’), we are able to classify different brain states by using the statistics across time of these structures hitherto hidden in the neuroimaging dynamics. Proving the strong potential of this framework, we were able to classify resting-state BOLD fMRI signals from two classes of post-comatose patients (minimally conscious state and unresponsive wakefulness syndrome) compared with healthy controls with very high precision.

Dynamic spectral signatures of mirror movements in the sensorimotor functional connectivity network of patients with Kallmann syndrome.

In Kallmann syndrome (KS), the peculiar phenomenon of bimanual synkinesis or mirror movement (MM) has been associated with a spectral shift, from lower to higher frequencies, of the resting-state fMRI signal of the large-scale sensorimotor brain network (SMN). To possibly determine whether a similar frequency specificity exists across different functional connectivity SMN states, and to capture spontaneous transitions between them, we investigated the dynamic spectral changes of the SMN functional connectivity in KS patients with and without MM symptom. Brain MRI data were acquired at 3 Tesla in 39 KS patients (32 without MM, KSMM-, seven with MM, KSMM+) and 26 age- and sex-matched healthy control (HC) individuals. The imaging protocol included 20-min rs-fMRI scans enabling detailed spectro-temporal analyses of large-scale functional connectivity brain networks. Group independent component analysis was used to extract the SMN. A sliding window approach was used to extract the dynamic spectral power of the SMN functional connectivity within the canonical physiological frequency range of slow rs-fMRI signal fluctuations (0.01–0.25 Hz). K-means clustering was used to determine (and count) the most recurrent dynamic states of the SMN and detect the number of transitions between them. Two most recurrent states were identified, for which the spectral power peaked at a relatively lower (state 1) and higher (state 2) frequency. Compared to KS patients without MM and HC subjects, the SMN of KS patients with MM displayed significantly larger spectral power changes in the slow 3 canonical sub-band (0.073–0.198 Hz) and significantly fewer transitions between state 1 (less recurrent) and state 2 (more recurrent). These findings demonstrate that the presence of MM in KS patients is associated with reduced spontaneous transitions of the SMN between dynamic functional connectivity states and a higher recurrence and an increased spectral power change of the high-frequency state. These results provide novel information about the large-scale brain functional dynamics that could help to understand the pathologic mechanisms of bimanual synkinesis in KS syndrome and, potentially, other neurological disorders where MM may also occur.

Covariance-based vs. correlation-based functional connectivity dissociates healthy aging from Alzheimer disease

Prior studies of aging and Alzheimer disease have evaluated resting state functional connectivity (FC) using either seed-based correlation (SBC) or independent component analysis (ICA), with a focus on particular functional systems. SBC and ICA both are insensitive to differences in signal amplitude. At the same time, accumulating evidence indicates that the amplitude of spontaneous BOLD signal fluctuations is physiologically meaningful. We systematically compared covariance-based FC, which is sensitive to amplitude, vs. correlation-based FC, which is not, in affected individuals and controls drawn from two cohorts of participants including autosomal dominant Alzheimer disease (ADAD), late onset Alzheimer disease (LOAD), and age-matched controls. Functional connectivity was computed over 222 regions of interest and group differences were evaluated in terms of components projected onto a space of lower dimension. Our principal observations are: (1) Aging is associated with global loss of resting state fMRI signal amplitude that is approximately uniform across resting state networks. (2) Thus, covariance FC measures decrease with age whereas correlation FC is relatively preserved in healthy aging. (3) In contrast, symptomatic ADAD and LOAD both lead to loss of spontaneous activity amplitude as well as severely degraded correlation structure. These results demonstrate a double dissociation between age vs. Alzheimer disease and the amplitude vs. correlation structure of resting state BOLD signals. Modeling results suggest that the AD-associated loss of correlation structure is attributable to a relative increase in the fraction of locally restricted as opposed to widely shared variance.

Frequency-Specific Analysis of the Dynamic Reconfiguration of the Brain in Patients with Schizophrenia.

The analysis of resting-state fMRI signals usually focuses on the low-frequency range/band (0.01–0.1 Hz), which does not cover all aspects of brain activity. Studies have shown that distinct frequency bands can capture unique fluctuations in brain activity, with high-frequency signals (>0.1 Hz) providing valuable information for the diagnosis of schizophrenia. We hypothesized that it is meaningful to study the dynamic reconfiguration of schizophrenia through different frequencies. Therefore, this study used resting-state functional magnetic resonance (RS-fMRI) data from 42 schizophrenia and 40 normal controls to investigate dynamic network reconfiguration in multiple frequency bands (0.01–0.25 Hz, 0.01–0.027 Hz, 0.027–0.073 Hz, 0.073–0.198 Hz, 0.198–0.25 Hz). Based on the time-varying dynamic network constructed for each frequency band, we compared the dynamic reconfiguration of schizophrenia and normal controls by calculating the recruitment and integration. The experimental results showed that the differences between schizophrenia and normal controls are observed in the full frequency, which is more significant in slow3. In addition, as visual network, attention network, and default mode network differ a lot from each other, they can show a high degree of connectivity, which indicates that the functional network of schizophrenia is affected by the abnormal brain state in these areas. These shreds of evidence provide a new perspective and promote the current understanding of the characteristics of dynamic brain networks in schizophrenia.

Gaining insight into the neural basis of resting-state fMRI signal

The blood oxygenation level-dependent (BOLD)-based resting-state functional magnetic resonance imaging (rsfMRI) has been widely used as a non-invasive tool to map brain-wide connectivity architecture. However, the neural basis underpinning the resting-state BOLD signal remains elusive. In this study, we combined simultaneous calcium-based fiber photometry with rsfMRI in awake animals to examine the relationship of the BOLD signal and spiking activity at the resting state. We observed robust couplings between calcium and BOLD signals in the dorsal hippocampus as well as other distributed areas in the default mode network (DMN), suggesting that the calcium measurement can reliably predict the rsfMRI signal. In addition, using the calcium signal recorded as the ground truth, we assessed the impacts of different rsfMRI data preprocessing pipelines on functional connectivity mapping. Overall, our results provide important evidence suggesting that spiking activity measured by the calcium signal plays a key role in the neural mechanism of resting-state BOLD signal.