Global and structured waves of rs-fMRI signal identified as putative propagation of spontaneous neural activity

The contributions of resting-state functional-MRI studies to our understanding of male patients with obstructive sleep apnea: a systematic review

Objectives

Obstructive sleep apnea (OSA) is a condition marked by the recurrent partial or complete obstruction of the upper airway during sleep. This leads to intermittent pauses in breathing, fragmented sleep, and frequent awakenings throughout the night. Many of these symptoms are believed to be linked to brain damage; however, the fundamental neurological processes underlying these impairments remain largely unknown. This study investigates resting-state functional MRI (rs-fMRI) findings in male patients with OSA to better understand the specific mechanisms associated with this condition in this demographic.

Methods

The search was conducted in the PubMed and Google Scholar databases, encompassing literature from their inception to July 2024. Studies were identified based on predetermined inclusion and exclusion criteria and were evaluated by two independent reviewers.

Results

A total of 16 eligible original rs-fMRI studies on male patients with OSA were included in this review. The findings indicate that patients with OSA exhibit alterations in resting-state brain activity. These neural changes may help explain the effects of OSA on emotion, cognition, and quality of life. Additionally, these findings could be used in the future to evaluate treatment outcomes.

Conclusion

This study highlights significant changes in local brain activities, interested region related functional connectivity, and whole-brain functional connectivity networks in patients with OSA. These findings offer valuable insights into the neural alterations at the core of OSA and may serve as potential biomarkers for diagnosis and intervention.

Self-supervised graph contrastive learning with diffusion augmentation for functional MRI analysis and brain disorder detection

Resting-state functional magnetic resonance imaging (rs-fMRI) provides a non-invasive imaging technique to study patterns of brain activity, and is increasingly used to facilitate automated brain disorder analysis. Existing fMRI-based learning methods often rely on labeled data to construct learning models, while the data annotation process typically requires significant time and resource investment. Graph contrastive learning offers a promising solution to address the small labeled data issue, by augmenting fMRI time series for self-supervised learning. However, data augmentation strategies employed in these approaches may damage the original blood-oxygen-level-dependent (BOLD) signals, thus hindering subsequent fMRI feature extraction. In this paper, we propose a self-supervised graph contrastive learning framework with diffusion augmentation (GCDA) for functional MRI analysis. The GCDA consists of a pretext model and a task-specific model. In the pretext model, we first augment each brain functional connectivity network derived from fMRI through a graph diffusion augmentation (GDA) module, and then use two graph isomorphism networks with shared parameters to extract features in a self-supervised contrastive learning manner. The pretext model can be optimized without the need for labeled training data, while the GDA focuses on perturbing graph edges and nodes, thus preserving the integrity of original BOLD signals. The task-specific model involves fine-tuning the trained pretext model to adapt to downstream tasks. Experimental results on two rs-fMRI cohorts with a total of 1230 subjects demonstrate the effectiveness of our method compared with several state-of-the-arts.

Reduced visual and middle temporal gyrus activity correlates with years of exercise in athletes using resting-state fMRI

BACKGROUND AND PURPOSE

Different types of physical training can lead to changes in brain activity and function, and these changes can vary depending on the type of training. However, it remains unclear whether there are commonalities in how different types of training affect brain activity and function. The purpose of this study is to compare the brain activity states of professional athletes with those of ordinary university students and to explore the relationship between training and differences in brain activity states.

METHODS

This study primarily utilizes resting-state MRI and the degree centrality metric to investigate spontaneous brain activity in 86 high-level athletes with extensive training and 74 age- and gender-matched nonathletes. Additionally, a correlation analysis between brain activity in relevant regions and years of training was conducted.

RESULTS

The analysis revealed that, compared to nonathletes, high-level athletes exhibited reduced activity in the Calcarine (a visual area) and Middle Temporal Gyrus. Furthermore, changes in the activity of the Calcarine and Middle Temporal Gyrus were significantly correlated with the number of years of professional training.

CONCLUSIONS

The study results indicate that long-term physical training is associated with changes in brain activity in athletes, providing insights into the neural mechanisms underlying behavioral performance in professional athletes.

Resting-State fMRI: Emerging Concepts for Future Clinical Application

Resting-state functional magnetic resonance imaging (rsfMRI) has been developed as a method of investigating spontaneous neural activity. Based on its low-frequency signal synchronization, rsfMRI has made it possible to identify multiple macroscopic structures termed resting-state networks (RSNs) on a single scan of less than 10 minutes. It is easy to implement even in clinical practice, in which assigning tasks to patients can be challenging. These advantages have accelerated the adoption and growth of rsfMRI. Recently, studies on the global rsfMRI signal have attracted increasing attention. Because it primarily arises from physiological events, less attention has hitherto been paid to the global signal than to the local network (i.e., RSN) component. However, the global signal is not a mere nuisance or a subsidiary component. On the contrary, it is quantitatively the dominant component that accounts for most of the variance in the rsfMRI signal throughout the brain and provides rich information on local hemodynamics that can serve as an individual-level diagnostic biomarker. Moreover, spatiotemporal analyses of the global signal have revealed that it is closely and fundamentally associated with the organization of RSNs, thus challenging the basic assumptions made in conventional rsfMRI analyses and views on RSNs. This review introduces new concepts emerging from rsfMRI spatiotemporal analyses focusing on the global signal and discusses how they may contribute to future clinical medicine. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: Stage 1.

Resting-state networks representation of the global phenomena

Resting-state functional magnetic resonance imaging (rsfMRI) has been widely applied to investigate spontaneous neural activity, often based on its macroscopic organization that is termed resting-state networks (RSNs). Although the neurophysiological mechanisms underlying the RSN organization remain largely unknown, accumulating evidence points to a substantial contribution from the global signals to their structured synchronization. This study further explored the phenomenon by taking advantage of the inter- and intra-subject variations of the time delay and correlation coefficient of the signal timeseries in each region using the global mean signal as the reference signal. Consistent with the hypothesis based on the empirical and theoretical findings, the time lag and correlation, which have consistently been proven to represent local hemodynamic status, were shown to organize networks equivalent to RSNs. The results not only provide further evidence that the local hemodynamic status could be the direct source of the RSNs' spatial patterns but also explain how the regional variations in the hemodynamics, combined with the changes in the global events' power spectrum, lead to the observations. While the findings pose challenges to interpretations of rsfMRI studies, they further support the view that rsfMRI can offer detailed information related to global neurophysiological phenomena as well as local hemodynamics that would have great potential as biomarkers.

Development of top-down cortical propagations in youth

Hierarchical processing requires activity propagating between higher- and lower-order cortical areas. However, functional neuroimaging studies have chiefly quantified fluctuations within regions over time rather than propagations occurring over space. Here, we leverage advances in neuroimaging and computer vision to track cortical activity propagations in a large sample of youth (n = 388). We delineate cortical propagations that systematically ascend and descend a cortical hierarchy in all individuals in our developmental cohort, as well as in an independent dataset of densely sampled adults. Further, we demonstrate that top-down, descending hierarchical propagations become more prevalent with greater demands for cognitive control as well as with development in youth. These findings emphasize that hierarchical processing is reflected in the directionality of propagating cortical activity and suggest top-down propagations as a potential mechanism of neurocognitive maturation in youth.

Reproducibility of cerebral perfusion measurements using BOLD delay

BOLD delay is an emerging, noninvasive method for assessing cerebral perfusion that does not require the use of intravenous contrast agents and is thus particularly suited for longitudinal monitoring. In this study, we assess the reproducibility of BOLD delay using data from 136 subjects with normal cerebral perfusion scanned on two separate occasions with scanners, sequence parameters, and intervals between scans varying between subjects. The effects of various factors on the reproducibility of BOLD delay, defined here as the differences in BOLD delay values between the scanning sessions, were investigated using a linear mixed model. Reproducibility was additionally assessed using the intraclass correlation coefficient of BOLD delay between sessions. Reproducibility was highest in the posterior cerebral artery territory. The mean BOLD delay test-retest difference after accounting for the aforementioned factors was 1.2 s (95% CI = 1.0 to 1.4 s). Overall, BOLD delay shows good reproducibility, but care should be taken when interpreting longitudinal BOLD delay changes that are either very small or are located in certain brain regions.

Reliability and Sensitivity to Alterered Hemodynamics Measured with Resting-state fMRI Metrics: Comparison with 123I-IMP SPECT

Blood oxygenation level-dependent (BOLD) contrast is sensitive to local hemodynamic changes and thus is applicable to imaging perfusion or vascular reactivity. However, knowledge about its measurement characteristics compared to reference standard perfusion imaging is limited. This study longitudinally evaluated perfusion in patients with steno-occlusive disease using resting-state functional MRI (rsfMRI) acquired before and within nine days of anterior circulation revascularization in patients with large cerebral artery steno-occlusive diseases. The reliability and sensitivity to longitudinal changes of rsfMRI temporal correlation (Rc) and time delay (TDc) relative to the cerebellar signal were examined voxel-wise in comparison with single-photon emission CT (SPECT) cerebral blood flow (CBF) using the within-subject standard deviation (Sw) and intraclass correlation coefficients (ICCs). For statistical comparisons, the standard deviation (SD) of longitudinal changes within the cerebellum, the number of voxels with significant changes in the left middle cerebral artery territory ipsilateral to surgery, and their average changes relative to the cerebellar SD were evaluated. The test-retest reliability of the fMRI metrics was also similarly evaluated using the human connectome project (HCP) healthy young adult dataset. The test-retest time interval was 31 ± 18 days. Test-retest reliability was significantly higher for SPECT (cerebellar SD: -2.59 ± 0.20) than for fMRI metrics (cerebellar SD: Rc, -2.34 ± 0.24, p = 0.04; TDc, -2.19 ± 0.21, p = 0.003). Sensitivity to postoperative changes, which was evaluated as the number of voxels, was significantly higher for fMRI TDc (8.78 ± 0.72) than for Rc (7.42 ± 1.48, p = 0.03) or SPECT CBF (6.88 ± 0.67, p < 0.001). The ratio between the average Rc, TDc, and SPECT CBF changes within the left MCA target region and cerebellar SD was also significantly higher for fMRI TDc (1.21 ± 0.79) than Rc (0.48 ± 0.94, p = 0.006) or SPECT CBF (0.23 ± 0.57, p = 0.001). The measurement variability of time delay was also larger than that of temporal correlation in HCP data within the cerebellum (t = -8.7, p < 0.001) or in the whole-brain (t = -27.4, p < 0.001) gray matter. These data suggest that fMRI time delay is more sensitive to the hemodynamic changes than SPECT CBF, although the reliability is lower. The implication for fMRI connectivity studies is that temporal correlation can be significantly decreased due to altered hemodynamics, even in cases with normal CBF.

Selective blockade of rat brain T-type calcium channels provides insights on neurophysiological basis of arousal dependent resting state functional magnetic resonance imaging signals

A number of studies point to slow (0.1–2 Hz) brain rhythms as the basis for the resting-state functional magnetic resonance imaging (rsfMRI) signal. Slow waves exist in the absence of stimulation, propagate across the cortex, and are strongly modulated by vigilance similar to large portions of the rsfMRI signal. However, it is not clear if slow rhythms serve as the basis of all neural activity reflected in rsfMRI signals, or just the vigilance-dependent components. The rsfMRI data exhibit quasi-periodic patterns (QPPs) that appear to increase in strength with decreasing vigilance and propagate across the brain similar to slow rhythms. These QPPs can complicate the estimation of functional connectivity (FC) via rsfMRI, either by existing as unmodeled signal or by inducing additional wide-spread correlation between voxel-time courses of functionally connected brain regions. In this study, we examined the relationship between cortical slow rhythms and the rsfMRI signal, using a well-established pharmacological model of slow wave suppression. Suppression of cortical slow rhythms led to significant reduction in the amplitude of QPPs but increased rsfMRI measures of intrinsic FC in rats. The results suggest that cortical slow rhythms serve as the basis of only the vigilance-dependent components (e.g., QPPs) of rsfMRI signals. Further attenuation of these non-specific signals enhances delineation of brain functional networks.

The Spatiotemporal Dynamics of Cerebral Autoregulation in Functional Magnetic Resonance Imaging

The thigh-cuff release (TCR) maneuver is a physiological challenge that is widely used to assess dynamic cerebral autoregulation (dCA). It is often applied in conjunction with Transcranial Doppler ultrasound (TCD), which provides temporal information of the global flow response in the brain. This established method can only yield very limited insights into the regional variability of dCA, whereas functional MRI (fMRI) has the ability to reveal the spatial distribution of flow responses in the brain with high spatial resolution. The aim of this study was to use whole-brain blood-oxygenation-level-dependent (BOLD) fMRI to characterize the spatiotemporal dynamics of the flow response to the TCR challenge, and thus pave the way toward mapping dCA in the brain. We used a data driven approach to derive a novel basis set that was then used to provide a voxel-wise estimate of the TCR associated haemodynamic response function (HRF TCR ). We found that the HRF TCR evolves with a specific spatiotemporal pattern, with gray and white matter showing an asynchronous response, which likely reflects the anatomical structure of cerebral blood supply. Thus, we propose that TCR challenge fMRI is a promising method for mapping spatial variability in dCA, which will likely prove to be clinically advantageous.

Pre-processing of Sub-millimeter GE-BOLD fMRI Data for Laminar Applications

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.

Altered spontaneous brain activity patterns in patients with hyperthyroidism exophthalmos using amplitude of low-frequency fluctuation: a resting-state fMRI study

AIM

To assess changed spontaneous brain activity in hyperthyroid exophthalmos (HE) patients by the amplitude of the low-frequency fluctuation (ALFF) method, and to analyze the correlation between brain activity and ALFF values in these patients.

METHODS

Totally 18 HE and 18 hyperthyroid non-exophthalmos (HNE) patients were enrolled. The participants were tested by resting-state functional magnetic resonance imaging, and receiver operating characteristic (ROC) curves were generated to classify the ALFF values of the study population. Pearson's correlation analysis was utilized to evaluate the relationship between the ALFF values obtained from different brain areas and clinical manifestations.

RESULTS

Contrary to HNE patients, we observed lower ALFF values in the left calcarine fissure and surrounding cortex (LCFSC) in HE patients. In the ROC curve analysis of the LCFSC, the area under the curve reflected a high degree of accuracy. In addition, there was positive correlation between mean ALFF values of the LCFSC and the best-corrected visual acuity of the affected eyes.

CONCLUSION

The study displays abnormal brain activity in LCFSC in patients with HE, which might suggest pathological mechanism of visual impairment of HE patients.

Brain Activity Fluctuations Propagate as Waves Traversing the Cortical Hierarchy

The brain exhibits highly organized patterns of spontaneous activity as measured by resting-state functional magnetic resonance imaging (fMRI) fluctuations that are being widely used to assess the brain's functional connectivity. Some evidence suggests that spatiotemporally coherent waves are a core feature of spontaneous activity that shapes functional connectivity, although this has been difficult to establish using fMRI given the temporal constraints of the hemodynamic signal. Here, we investigated the structure of spontaneous waves in human fMRI and monkey electrocorticography. In both species, we found clear, repeatable, and directionally constrained activity waves coursed along a spatial axis approximately representing cortical hierarchical organization. These cortical propagations were closely associated with activity changes in distinct subcortical structures, particularly those related to arousal regulation, and modulated across different states of vigilance. The findings demonstrate a neural origin of spatiotemporal fMRI wave propagation at rest and link it to the principal gradient of resting-state fMRI connectivity.

Origin of the Time Lag Phenomenon and the Global Signal in Resting-State fMRI

The global mean signal of resting-state fMRI (rs-fMRI) shows a characteristic spatiotemporal pattern that is closely related to the pattern of vascular perfusion. Although being increasingly adopted in the mapping of the flow of neural activity, the mechanism that gives rise to the BOLD signal time lag remains controversial. In the present study, we compared the time lag of the global mean signal with those of the local network components obtained by applying temporal independent component analysis to the resting-state fMRI data, as well as by using simultaneous wide-field visual stimulation, and demonstrated that the time lag patterns are highly similar across all types of data. These results suggest that the time lag of the rs-fMRI signal reflects the local variance of the hemodynamic responses rather than the arrival or transit time of the stimulus, whether the trigger is neuronal or non-neuronal in origin as long as it is mediated by local hemodynamic responses. Examinations of the internal carotid artery signal further confirmed that the arterial signal is tightly inversely coupled with the global mean signal in accordance with previous studies, presumably reflecting the blood flow or blood pressure changes that are occurring almost simultaneously in the internal carotid artery and the cerebral pial/capillary arteries, within the low-frequency component in human rs-fMRI.

Lag Analysis of Fast fMRI Reveals Delayed Information Flow Between the Default Mode and Other Networks in Narcolepsy

Narcolepsy is a chronic neurological disease characterized by dysfunction of the hypocretin system in brain causing disruption in the wake-promoting system. In addition to sleep attacks and cataplexy, patients with narcolepsy commonly report cognitive symptoms while objective deficits in sustained attention and executive function have been observed. Prior resting-state functional magnetic resonance imaging (fMRI) studies in narcolepsy have reported decreased inter/intranetwork connectivity regarding the default mode network (DMN). Recently developed fast fMRI data acquisition allows more precise detection of brain signal propagation with a novel dynamic lag analysis. In this study, we used fast fMRI data to analyze dynamics of inter resting-state network (RSN) information signaling between narcolepsy type 1 patients (NT1, n = 23) and age- and sex-matched healthy controls (HC, n = 23). We investigated dynamic connectivity properties between positive and negative peaks and, furthermore, their anticorrelative (pos-neg) counterparts. The lag distributions were significantly (P < 0.005, familywise error rate corrected) altered in 24 RSN pairs in NT1. The DMN was involved in 83% of the altered RSN pairs. We conclude that narcolepsy type 1 is characterized with delayed and monotonic inter-RSN information flow especially involving anticorrelations, which are known to be characteristic behavior of the DMN regarding neurocognition.

Organization of Propagated Intrinsic Brain Activity in Individual Humans

Spontaneous infra-slow (<0.1 Hz) fluctuations in functional magnetic resonance imaging (fMRI) signals are temporally correlated within large-scale functional brain networks, motivating their use for mapping systems-level brain organization. However, recent electrophysiological and hemodynamic evidence suggest state-dependent propagation of infra-slow fluctuations, implying a functional role for ongoing infra-slow activity. Crucially, the study of infra-slow temporal lag structure has thus far been limited to large groups, as analyzing propagation delays requires extensive data averaging to overcome sampling variability. Here, we use resting-state fMRI data from 11 extensively-sampled individuals to characterize lag structure at the individual level. In addition to stable individual-specific features, we find spatiotemporal topographies in each subject similar to the group average. Notably, we find a set of early regions that are common to all individuals, are preferentially positioned proximal to multiple functional networks, and overlap with brain regions known to respond to diverse behavioral tasks-altogether consistent with a hypothesized ability to broadly influence cortical excitability. Our findings suggest that, like correlation structure, temporal lag structure is a fundamental organizational property of resting-state infra-slow activity.

Global vs. Network-Specific Regulations as the Source of Intrinsic Coactivations in Resting-State Networks

Spontaneous neural activities are endowed with specific patterning characterized by synchronizations within functionally relevant distant regions that are termed as resting-state networks (RSNs). Although the mechanisms that organize the large-scale neural systems are still largely unknown, recent studies have proposed a hypothesis that network-specific coactivations indeed emerge as the result of globally propagating neural activities with specific paths of transmission. However, the extent to which such a centralized global regulation, rather than network-specific control, contributes to the RSN synchronization remains unknown. In the present study, we investigated the contribution from each mechanism by directly identifying the global as well as local component of resting-state functional MRI (fMRI) data provided by human connectome project, using temporal independent component analysis (ICA). Based on the spatial distribution pattern, each ICA component was classified as global or local. Time lag mapping of each IC revealed several paths of global or semi-global propagations that are partially overlapping yet spatially distinct to each other. Consistent with previous studies, the time lag of global oscillation, although being less spatially homogenous than what was assumed to be, contributed to the RSN synchronization. However, an equivalent contribution was also shown on the part of the more locally confined activities that are independent to each other. While allowing the view that network-specific coactivation occurs as part of the sequences of global neural activities, these results further confirm an equally important role of the network-specific regulation for its coactivation, regardless of whether vascular artifacts contaminate the global component in fMRI measures.

Axial variation of deoxyhemoglobin density as a source of the low-frequency time lag structure in blood oxygenation level-dependent signals

Perfusion-related information is reportedly embedded in the low-frequency component of a blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal. The blood-propagation pattern through the cerebral vascular tree is detected as an interregional lag variation of spontaneous low-frequency oscillations (sLFOs). Mapping of this lag, or phase, has been implicitly treated as a projection of the vascular tree structure onto real space. While accumulating evidence supports the biological significance of this signal component, the physiological basis of the “perfusion lag structure,” a requirement for an integrative resting-state fMRI-signal model, is lacking. In this study, we conducted analyses furthering the hypothesis that the sLFO is not only largely of systemic origin, but also essentially intrinsic to blood, and hence behaves as a virtual tracer. By summing the small fluctuations of instantaneous phase differences between adjacent vascular regions, a velocity response to respiratory challenges was detected. Regarding the relationship to neurovascular coupling, the removal of the whole lag structure, which can be considered as an optimized global-signal regression, resulted in a reduction of inter-individual variance while preserving the fMRI response. Examination of the T2* and S₀, or non-BOLD, components of the fMRI signal revealed that the lag structure is deoxyhemoglobin dependent, while paradoxically presenting a signal-magnitude reduction in the venous side of the cerebral vasculature. These findings provide insight into the origin of BOLD sLFOs, suggesting that they are highly intrinsic to the circulating blood.

The Not-So-Global Blood Oxygen Level-Dependent Signal

Global signal regression is a controversial processing step for resting-state functional magnetic resonance imaging, partly because the source of the global blood oxygen level-dependent (BOLD) signal remains unclear. On the one hand, nuisance factors such as motion can readily introduce coherent BOLD changes across the whole brain. On the other hand, the global signal has been linked to neural activity and vigilance levels, suggesting that it contains important neurophysiological information and should not be discarded. Any widespread pattern of coordinated activity is likely to contribute appreciably to the global signal. Such patterns may include large-scale quasiperiodic spatiotemporal patterns, known also to be tied to performance on vigilance tasks. This uncertainty surrounding the separability of the global BOLD signal from concurrent neurological processes motivated an examination of the global BOLD signal's spatial distribution. The results clarify that although the global signal collects information from all tissue classes, a diverse subset of the BOLD signal's independent components contribute the most to the global signal. Further, the timing of each network's contribution to the global signal is not consistent across volunteers, confirming the independence of a constituent process that comprises the global signal.

Altered dynamics of brain segregation and integration in poststroke aphasia

Poststroke aphasia (PSA) results from direct effect of focal lesions and dysfunction of distributed language networks. However, how flexible the activity at specific nodes control global dynamics is currently unknown. In this study, we demonstrate that alterations in the regional activity may cause imbalances between segregation and integration in temporo-spatial pattern, and the transient dynamics are disrupted in PSA patients. Specifically, we applied dynamic framework to eyes-closed resting-state functional MRI data from PSA patients (n = 17), and age-, gender-, and education-matched healthy controls (HCs, n = 20). Subsequently, we calculated two basis brain organizational principles: "dynamic segregation," obtained from dynamic amplitude of low-frequency fluctuations (dALFF), which represent the specialized processing within interconnected brain regions; and "dynamic integration," obtained from dynamic functional connectivity, which measures the efficient communication between interconnected brain regions. We found that both measures were decreased in the PSA patients within the left frontal and temporal subregions compared to the HCs. PSA patients displayed increased flexibility of interaction between left temporo-frontal subregions and right temporo-parieto-frontal subnetworks. Furthermore, we found that dALFF in the pars triangularis of left inferior frontal gyrus was associated with aphasia quotient. These findings suggest that the reduced temporal flexibility of regional activity in language-relevant cortical regions in PSA is related to the disrupted organization of intrahemispheric networks, leading to a loss of the corresponding functions. By using dynamic framework, our results offer valuable information about the alterations in segregation and integration of spatiotemporal information across networks and illuminate how dysfunction in flexible activity may underlie language deficits in PSA.

On time delay estimation and sampling error in resting-state fMRI

Accumulating evidence indicates that resting-state functional magnetic resonance imaging (rsfMRI) signals correspond to propagating electrophysiological infra-slow activity (<0.1 Hz). Thus, pairwise correlations (zero-lag functional connectivity (FC)) and temporal delays among regional rsfMRI signals provide useful, complementary descriptions of spatiotemporal structure in infra-slow activity. However, the slow nature of fMRI signals implies that practical scan durations cannot provide sufficient independent temporal samples to stabilize either of these measures. Here, we examine factors affecting sampling variability in both time delay estimation (TDE) and FC. Although both TDE and FC accuracy are highly sensitive to data quantity, we use surrogate fMRI time series to study how the former is additionally related to the magnitude of a given pairwise correlation and, to a lesser extent, the temporal sampling rate. These contingencies are further explored in real data comprising 30-min rsfMRI scans, where sampling error (i.e., limited accuracy owing to insufficient data quantity) emerges as a significant but underappreciated challenge to FC and, even more so, to TDE. Exclusion of high-motion epochs exacerbates sampling error; thus, both sides of the bias-variance (or data quality-quantity) tradeoff associated with data exclusion should be considered when analyzing rsfMRI data. Finally, we present strategies for TDE in motion-corrupted data, for characterizing sampling error in TDE and FC, and for mitigating the influence of sampling error on lag-based analyses.

How networks communicate: propagation patterns in spontaneous brain activity

Initially regarded as 'noise', spontaneous (intrinsic) activity accounts for a large portion of the brain's metabolic cost. Moreover, it is now widely known that infra-slow (less than 0.1 Hz) spontaneous activity, measured using resting state functional magnetic resonance imaging of the blood oxygen level-dependent (BOLD) signal, is correlated within functionally defined resting state networks (RSNs). However, despite these advances, the temporal organization of spontaneous BOLD fluctuations has remained elusive. By studying temporal lags in the resting state BOLD signal, we have recently shown that spontaneous BOLD fluctuations consist of remarkably reproducible patterns of whole brain propagation. Embedded in these propagation patterns are unidirectional 'motifs' which, in turn, give rise to RSNs. Additionally, propagation patterns are markedly altered as a function of state, whether physiological or pathological. Understanding such propagation patterns will likely yield deeper insights into the role of spontaneous activity in brain function in health and disease.This article is part of the themed issue 'Interpreting blood oxygen level-dependent: a dialogue between cognitive and cellular neuroscience'.

Low-frequency hippocampal-cortical activity drives brain-wide resting-state functional MRI connectivity

The hippocampus, including the dorsal dentate gyrus (dDG), and cortex engage in bidirectional communication. We propose that low-frequency activity in hippocampal-cortical pathways contributes to brain-wide resting-state connectivity to integrate sensory information. Using optogenetic stimulation and brain-wide fMRI and resting-state fMRI (rsfMRI), we determined the large-scale effects of spatiotemporal-specific downstream propagation of hippocampal activity. Low-frequency (1 Hz), but not high-frequency (40 Hz), stimulation of dDG excitatory neurons evoked robust cortical and subcortical brain-wide fMRI responses. More importantly, it enhanced interhemispheric rsfMRI connectivity in various cortices and hippocampus. Subsequent local field potential recordings revealed an increase in slow oscillations in dorsal hippocampus and visual cortex, interhemispheric visual cortical connectivity, and hippocampal-cortical connectivity. Meanwhile, pharmacological inactivation of dDG neurons decreased interhemispheric rsfMRI connectivity. Functionally, visually evoked fMRI responses in visual regions also increased during and after low-frequency dDG stimulation. Together, our results indicate that low-frequency activity robustly propagates in the dorsal hippocampal-cortical pathway, drives interhemispheric cortical rsfMRI connectivity, and mediates visual processing.

Frequency-phase analysis of resting-state functional MRI

We describe an analysis method that characterizes the correlation between coupled time-series functions by their frequencies and phases. It provides a unified framework for simultaneous assessment of frequency and latency of a coupled time-series. The analysis is demonstrated on resting-state functional MRI data of 34 healthy subjects. Interactions between fMRI time-series are represented by cross-correlation (with time-lag) functions. A general linear model is used on the cross-correlation functions to obtain the frequencies and phase-differences of the original time-series. We define symmetric, antisymmetric and asymmetric cross-correlation functions that correspond respectively to in-phase, 90° out-of-phase and any phase difference between a pair of time-series, where the last two were never introduced before. Seed maps of the motor system were calculated to demonstrate the strength and capabilities of the analysis. Unique types of functional connections, their dominant frequencies and phase-differences have been identified. The relation between phase-differences and time-delays is shown. The phase-differences are speculated to inform transfer-time and/or to reflect a difference in the hemodynamic response between regions that are modulated by neurotransmitters concentration. The analysis can be used with any coupled functions in many disciplines including electrophysiology, EEG or MEG in neuroscience.

Combined spatiotemporal ICA (stICA) for continuous and dynamic lag structure analysis of MREG data

This study investigated lag structure in the resting-state fMRI by applying a novel independent component (ICA) method to magnetic resonance encephalography (MREG) data. Briefly, the spatial ICA (sICA) was used for defining the frontal and back nodes of the default mode network (DMN), and the temporal ICA (tICA), which is enabled by the high temporal resolution of MREG (TR=100ms), was used to separate both neuronal and physiological components of these two spatial map regions. Subsequently, lag structure was investigated between the frontal (DMNvmpf) and posterior (DMNpcc) DMN nodes using both conventional method with all-time points and a sliding-window approach. A rigorous noise exclusion criterion was applied for tICs to remove physiological pulsations, motion and system artefacts. All the de-noised tICs were used to calculate the null-distributions both for expected lag variability over time and over subjects. Lag analysis was done for the three highest correlating denoised tICA pairs. Mean time lag of 0.6s (± 0.5 std) and mean absolute correlation of 0.69 (± 0.08) between the highest correlating tICA pairs of DMN nodes was observed throughout the whole analyzed period. In dynamic 2min window analysis, there was large variability over subjects as ranging between 1-10sec. Directionality varied between these highly correlating sources an average 28.8% of the possible number of direction changes. The null models show highly consistent correlation and lag structure between DMN nodes both in continuous and dynamic analysis. The mean time lag of a null-model over time between all denoised DMN nodes was 0.0s and, thus the probability of having either DMNpcc or DMNvmpf as a preceding component is near equal. All the lag values of highest correlating tICA pairs over subjects lie within the standard deviation range of a null-model in whole time window analysis, supporting the earlier findings that there is a consistent temporal lag structure across groups of individuals. However, in dynamic analysis, there are lag values exceeding the threshold of significance of a null-model meaning that there might be biologically meaningful variation in this measure. Taken together the variability in lag and the presence of high activity peaks during strong connectivity indicate that individual avalanches may play an important role in defining dynamic independence in resting state connectivity within networks.

The Effects of Acetazolamide on the Evaluation of Cerebral Hemodynamics and Functional Connectivity Using Blood Oxygen Level-Dependent MR Imaging in Patients with Chronic Steno-Occlusive Disease of the Anterior Circulation

BACKGROUND AND PURPOSE

Measuring cerebrovascular reactivity with the use of vasodilatory stimuli, such as acetazolamide, is useful for chronic cerebrovascular steno-occlusive disease. The purpose of this study was to evaluate the effects of acetazolamide on the assessment of hemodynamic impairment and functional connectivity by using noninvasive resting-state blood oxygen level-dependent MR imaging.

MATERIALS AND METHODS

A 20-minute resting-state blood oxygen level-dependent MR imaging scan was acquired with infusion of acetazolamide starting at 5 minutes after scan initiation. A recently developed temporal-shift analysis technique was applied on blood oxygen level-dependent MR imaging data before and after acetazolamide infusion to identify regions with hemodynamic impairment, and the results were compared by using contrast agent-based DSC perfusion imaging as the reference standard. Functional connectivity was compared with and without correction on the signal by using information from temporal-shift analysis, before and after acetazolamide infusion.

RESULTS

Visually, temporal-shift analysis of blood oxygen level-dependent MR imaging data identified regions with compromised hemodynamics as defined by DSC, though performance deteriorated in patients with bilateral disease. The Dice similarity coefficient between temporal-shift and DSC maps was higher before (0.487 ± 0.150 by using the superior sagittal sinus signal as a reference for temporal-shift analysis) compared with after acetazolamide administration (0.384 ± 0.107) (P = .006, repeated-measures ANOVA). Functional connectivity analysis with temporal-shift correction identified brain network nodes that were otherwise missed. The accuracy of functional connectivity assessment decreased after acetazolamide administration (P = .015 for default mode network, repeated-measures ANOVA).

CONCLUSIONS

Temporal-shift analysis of blood oxygen level-dependent MR imaging can identify brain regions with hemodynamic compromise in relation to DSC among patients with chronic cerebrovascular disease. The use of acetazolamide reduces the accuracy of temporal-shift analysis and network connectivity evaluation.