Quasi-periodic patterns (QPP): large-scale dynamics in resting state fMRI that correlate with local infraslow electrical activity

Brain states as wave-like motifs

There is ample evidence of wave-like activity in the brain at multiple scales and levels. This emerging literature supports the broader adoption of a wave perspective of brain activity. Specifically, a brain state can be described as a set of recurring, sequential patterns of propagating brain activity, namely a wave. We examine a collective body of experimental work investigating wave-like properties. Based on these works, we consider brain states as waves using a scale-agnostic framework across time and space. Emphasis is placed on the sequentiality and periodicity associated with brain activity. We conclude by discussing the implications, prospects, and experimental opportunities of this framework.

Variation in the Distribution of Large-scale Spatiotemporal Patterns of Activity Across Brain States

A few large-scale spatiotemporal patterns of brain activity (quasiperiodic patterns or QPPs) account for most of the spatial structure observed in resting state functional magnetic resonance imaging (rs-fMRI). The QPPs capture well-known features such as the evolution of the global signal and the alternating dominance of the default mode and task positive networks. These widespread patterns of activity have plausible ties to neuromodulatory input that mediates changes in nonlocalized processes, including arousal and attention. To determine whether QPPs exhibit variations across brain conditions, the relative magnitude and distribution of the three strongest QPPs were examined in two scenarios. First, in data from the Human Connectome Project, the relative incidence and magnitude of the QPPs was examined over the course of the scan, under the hypothesis that increasing drowsiness would shift the expression of the QPPs over time. Second, using rs-fMRI in rats obtained with a novel approach that minimizes noise, the relative incidence and magnitude of the QPPs was examined under three different anesthetic conditions expected to create distinct types of brain activity. The results indicate that both the distribution of QPPs and their magnitude changes with brain state, evidence of the sensitivity of these large-scale patterns to widespread changes linked to alterations in brain conditions.

Rhythmicity of neuronal oscillations delineates their cortical and spectral architecture

Neuronal oscillations are commonly analyzed with power spectral methods that quantify signal amplitude, but not rhythmicity or 'oscillatoriness' per se. Here we introduce a new approach, the phase-autocorrelation function (pACF), for the direct quantification of rhythmicity. We applied pACF to human intracerebral stereoelectroencephalography (SEEG) and magnetoencephalography (MEG) data and uncovered a spectrally and anatomically fine-grained cortical architecture in the rhythmicity of single- and multi-frequency neuronal oscillations. Evidencing the functional significance of rhythmicity, we found it to be a prerequisite for long-range synchronization in resting-state networks and to be dynamically modulated during event-related processing. We also extended the pACF approach to measure 'burstiness' of oscillatory processes and characterized regions with stable and bursty oscillations. These findings show that rhythmicity is double-dissociable from amplitude and constitutes a functionally relevant and dynamic characteristic of neuronal oscillations.

Global spatiotemporal synchronizing structures of spontaneous neural activities in different cell types

Increasing evidence has revealed the large-scale nonstationary synchronizations as traveling waves in spontaneous neural activity. However, the interplay of various cell types in fine-tuning these spatiotemporal patters remains unclear. Here, we performed comprehensive exploration of spatiotemporal synchronizing structures across different cell types, states (awake, anesthesia, motion) and developmental axis in male mice. We found traveling waves in glutamatergic neurons exhibited greater variety than those in GABAergic neurons. Moreover, the synchronizing structures of GABAergic neurons converged toward those of glutamatergic neurons during development, but the evolution of waves exhibited varying timelines for different sub-type interneurons. Functional connectivity arises from both standing and traveling waves, and negative connections can be elucidated by the spatial propagation of waves. In addition, some traveling waves were correlated with the spatial distribution of gene expression. Our findings offer further insights into the neural underpinnings of traveling waves, functional connectivity, and resting-state networks, with cell-type specificity and developmental perspectives.

Altered resting-state functional connectivity and dynamic network properties in cognitive impairment: an independent component and dominant-coactivation pattern analyses study

Introduction

Cognitive impairment (CI) due to Alzheimer's disease (AD) encompasses a decline in cognitive abilities and can significantly impact an individual's quality of life. Early detection and intervention are crucial in managing CI, both in the preclinical and prodromal stages of AD prior to dementia.

Methods

In this preliminary study, we investigated differences in resting-state functional connectivity and dynamic network properties between 23 individual with CI due to AD based on clinical assessment and 15 healthy controls (HC) using Independent Component Analysis (ICA) and Dominant-Coactivation Pattern (d-CAP) analysis. The cognitive status of the two groups was also compared, and correlations between cognitive scores and d-CAP switching probability were examined.

Results

Results showed comparable numbers of d-CAPs in the Default Mode Network (DMN), Executive Control Network (ECN), and Frontoparietal Network (FPN) between HC and CI groups. However, the Visual Network (VN) exhibited fewer d-CAPs in the CI group, suggesting altered dynamic properties of this network for the CI group. Additionally, ICA revealed significant connectivity differences for all networks. Spatial maps and effect size analyses indicated increased coactivation and more synchronized activity within the DMN in HC compared to CI. Furthermore, reduced switching probabilities were observed for the CI group in DMN, VN, and FPN networks, indicating less dynamic and flexible functional interactions.

Discussion

The findings highlight altered connectivity patterns within the DMN, VN, ECN, and FPN, suggesting the involvement of multiple functional networks in CI. Understanding these brain processes may contribute to developing targeted diagnostic and therapeutic strategies for CI due to AD.

Time-varying functional connectivity predicts fluctuations in sustained attention in a serial tapping task

The mechanisms for how large-scale brain networks contribute to sustained attention are unknown. Attention fluctuates from moment to moment, and this continuous change is consistent with dynamic changes in functional connectivity between brain networks involved in the internal and external allocation of attention. In this study, we investigated how brain network activity varied across different levels of attentional focus (i.e., "zones"). Participants performed a finger-tapping task, and guided by previous research, in-the-zone performance or state was identified by low reaction time variability and out-of-the-zone as the inverse. In-the-zone sessions tended to occur earlier in the session than out-of-the-zone blocks. This is unsurprising given the way attention fluctuates over time. Employing a novel method of time-varying functional connectivity, called the quasi-periodic pattern analysis (i.e., reliable, network-level low-frequency fluctuations), we found that the activity between the default mode network (DMN) and task positive network (TPN) is significantly more anti-correlated during in-the-zone states versus out-of-the-zone states. Furthermore, it is the frontoparietal control network (FPCN) switch that differentiates the two zone states. Activity in the dorsal attention network (DAN) and DMN were desynchronized across both zone states. During out-of-the-zone periods, FPCN synchronized with DMN, while during in-the-zone periods, FPCN switched to synchronized with DAN. In contrast, the ventral attention network (VAN) synchronized more closely with DMN during in-the-zone periods compared with out-of-the-zone periods. These findings demonstrate that time-varying functional connectivity of low frequency fluctuations across different brain networks varies with fluctuations in sustained attention or other processes that change over time.

Intrinsic forebrain arousal dynamics governs sensory stimulus encoding

The neural encoding of sensory stimuli is subject to the brain's internal circuit dynamics. Recent work has demonstrated that the resting brain exhibits widespread, coordinated activity that plays out over multisecond timescales in the form of quasi-periodic spiking cascades. Here we demonstrate that these intrinsic dynamics persist during the presentation of visual stimuli and markedly influence the efficacy of feature encoding in the visual cortex. During periods of passive viewing, the sensory encoding of visual stimuli was determined by quasi-periodic cascade cycle evolving over several seconds. During this cycle, high efficiency encoding occurred during peak arousal states, alternating in time with hippocampal ripples, which were most frequent in low arousal states. However, during bouts of active locomotion, these arousal dynamics were abolished: the brain remained in a state in which visual coding efficiency remained high and ripples were absent. We hypothesize that the brain's observed dynamics during awake, passive viewing reflect an adaptive cycle of alternating exteroceptive sensory sampling and internal mnemonic function.

QPPLab: A generally applicable software package for detecting, analyzing, and visualizing large-scale quasiperiodic spatiotemporal patterns (QPPs) of brain activity

One prominent brain dynamic process detected in functional neuroimaging data is large-scale quasi-periodic patterns (QPPs) which display spatiotemporal propagations along brain cortical gradients. QPP associates with the infraslow neural activity related to attention and arousal fluctuations and has been identified in both resting and task-evoked brains across various species. Several QPP detection and analysis tools were developed for distinct applications with study-specific parameter methods. This MATLAB package provides a simplified and user-friendly generally applicable toolbox for detecting, analyzing, and visualizing QPPs from fMRI timeseries of the brain. This paper describes the software functions and presents its ease of use on any brain datasets.

Metadata

[Table: see text].

The interaction between random and systematic visual stimulation and infraslow quasiperiodic spatiotemporal patterns of whole brain activity

One prominent feature of the infraslow BOLD signal during rest or task is quasi-periodic spatiotemporal pattern (QPP) of signal changes that involves an alternation of activity in key functional networks and propagation of activity across brain areas, and that is known to tie to the infraslow neural activity involved in attention and arousal fluctuations. This ongoing whole-brain pattern of activity might potentially modify the response to incoming stimuli or be modified itself by the induced neural activity. To investigate this, we presented checkerboard sequences flashing at 6Hz to subjects. This is a salient visual stimulus that is known to produce a strong response in visual processing regions. Two different visual stimulation sequences were employed, a systematic stimulation sequence in which the visual stimulus appeared every 20.3 secs and a random stimulation sequence in which the visual stimulus occurred randomly every 14~62.3 secs. Three central observations emerged. First, the two different stimulation conditions affect the QPP waveform in different aspects, i.e., systematic stimulation has greater effects on its phase and random stimulation has greater effects on its magnitude. Second, the QPP was more frequent in the systematic condition with significantly shorter intervals between consecutive QPPs compared to the random condition. Third, the BOLD signal response to the visual stimulus across both conditions was swamped by the QPP at the stimulus onset. These results provide novel insights into the relationship between intrinsic patterns and stimulated brain activity.

Manifold learning for fMRI time-varying functional connectivity

Whole-brain functional connectivity (FC) measured with functional MRI (fMRI) evolves over time in meaningful ways at temporal scales going from years (e.g., development) to seconds [e.g., within-scan time-varying FC (tvFC)]. Yet, our ability to explore tvFC is severely constrained by its large dimensionality (several thousands). To overcome this difficulty, researchers often seek to generate low dimensional representations (e.g., 2D and 3D scatter plots) hoping those will retain important aspects of the data (e.g., relationships to behavior and disease progression). Limited prior empirical work suggests that manifold learning techniques (MLTs)-namely those seeking to infer a low dimensional non-linear surface (i.e., the manifold) where most of the data lies-are good candidates for accomplishing this task. Here we explore this possibility in detail. First, we discuss why one should expect tvFC data to lie on a low dimensional manifold. Second, we estimate what is the intrinsic dimension (ID; i.e., minimum number of latent dimensions) of tvFC data manifolds. Third, we describe the inner workings of three state-of-the-art MLTs: Laplacian Eigenmaps (LEs), T-distributed Stochastic Neighbor Embedding (T-SNE), and Uniform Manifold Approximation and Projection (UMAP). For each method, we empirically evaluate its ability to generate neuro-biologically meaningful representations of tvFC data, as well as their robustness against hyper-parameter selection. Our results show that tvFC data has an ID that ranges between 4 and 26, and that ID varies significantly between rest and task states. We also show how all three methods can effectively capture subject identity and task being performed: UMAP and T-SNE can capture these two levels of detail concurrently, but LE could only capture one at a time. We observed substantial variability in embedding quality across MLTs, and within-MLT as a function of hyper-parameter selection. To help alleviate this issue, we provide heuristics that can inform future studies. Finally, we also demonstrate the importance of feature normalization when combining data across subjects and the role that temporal autocorrelation plays in the application of MLTs to tvFC data. Overall, we conclude that while MLTs can be useful to generate summary views of labeled tvFC data, their application to unlabeled data such as resting-state remains challenging.

Longitudinal investigation of changes in resting-state co-activation patterns and their predictive ability in the zQ175 DN mouse model of Huntington's disease

Huntington's disease (HD) is a neurodegenerative disorder caused by expanded (≥ 40) glutamine-encoding CAG repeats in the huntingtin gene, which leads to dysfunction and death of predominantly striatal and cortical neurons. While the genetic profile and clinical signs and symptoms of the disease are better known, changes in the functional architecture of the brain, especially before the clinical expression becomes apparent, are not fully and consistently characterized. In this study, we sought to uncover functional changes in the brain in the heterozygous (HET) zQ175 delta-neo (DN) mouse model at 3, 6, and 10 months of age, using resting-state functional magnetic resonance imaging (RS-fMRI). This mouse model shows molecular, cellular and circuitry alterations that worsen through age. Motor function disturbances are manifested in this model at 6 and 10 months of age. Specifically, we investigated, longitudinally, changes in co-activation patterns (CAPs) that are the transient states of brain activity constituting the resting-state networks (RSNs). Most robust changes in the temporal properties of CAPs occurred at the 10-months time point; the durations of two anti-correlated CAPs, characterized by simultaneous co-activation of default-mode like network (DMLN) and co-deactivation of lateral-cortical network (LCN) and vice-versa, were reduced in the zQ175 DN HET animals compared to the wild-type mice. Changes in the spatial properties, measured in terms of activation levels of different brain regions, during CAPs were found at all three ages and became progressively more pronounced at 6-, and 10 months of age. We then assessed the cross-validated predictive power of CAP metrics to distinguish HET animals from controls. Spatial properties of CAPs performed significantly better than the chance level at all three ages with 80% classification accuracy at 6 and 10 months of age.

Interacting spiral wave patterns underlie complex brain dynamics and are related to cognitive processing

The large-scale activity of the human brain exhibits rich and complex patterns, but the spatiotemporal dynamics of these patterns and their functional roles in cognition remain unclear. Here by characterizing moment-by-moment fluctuations of human cortical functional magnetic resonance imaging signals, we show that spiral-like, rotational wave patterns (brain spirals) are widespread during both resting and cognitive task states. These brain spirals propagate across the cortex while rotating around their phase singularity centres, giving rise to spatiotemporal activity dynamics with non-stationary features. The properties of these brain spirals, such as their rotational directions and locations, are task relevant and can be used to classify different cognitive tasks. We also demonstrate that multiple, interacting brain spirals are involved in coordinating the correlated activations and de-activations of distributed functional regions; this mechanism enables flexible reconfiguration of task-driven activity flow between bottom-up and top-down directions during cognitive processing. Our findings suggest that brain spirals organize complex spatiotemporal dynamics of the human brain and have functional correlates to cognitive processing.

The interaction between random and systematic visual stimulation and infraslow quasiperiodic spatiotemporal patterns of whole brain activity

One prominent feature of the infraslow BOLD signal during rest or task is quasi-periodic spatiotemporal pattern (QPP) of signal changes that involves an alternation of activity in key functional networks and propagation of activity across brain areas, and that is known to tie to the infraslow neural activity involved in attention and arousal fluctuations. This ongoing whole-brain pattern of activity might potentially modify the response to incoming stimuli or be modified itself by the induced neural activity. To investigate this, we presented checkerboard sequences flashing at 6Hz to subjects. This is a salient visual stimulus that is known to produce a strong response in visual processing regions. Two different visual stimulation sequences were employed, a systematic stimulation sequence in which the visual stimulus appeared every 20.3 secs and a random stimulation sequence in which the visual stimulus occurred randomly every 14∼62.3 secs. Three central observations emerged. First, the two different stimulation conditions affect the QPP waveform in different aspects, i.e., systematic stimulation has greater effects on its phase and random stimulation has greater effects on its magnitude. Second, the QPP was more frequent in the systematic condition with significantly shorter intervals between consecutive QPPs compared to the random condition. Third, the BOLD signal response to the visual stimulus across both conditions was swamped by the QPP at the stimulus onset. These results provide novel insights into the relationship between intrinsic patterns and stimulated brain activity.

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.

Sex-specific age-related changes in glymphatic function assessed by resting-state functional magnetic resonance imaging

The glymphatic system that clears out brain wastes, such as amyloid-β (Aβ) and tau, through cerebrospinal fluid (CSF) flow may play an important role in aging and dementias. However, a lack of non-invasive tools to assess the glymphatic function in humans hindered the understanding of the glymphatic 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 used to measure the glymphatic process and found to correlate with various pathologies of Alzheimer's disease (AD), including Aβ pathology. Using resting-state fMRI data from a large group of 719 healthy aging participants, we examined the sex-specific changes of the gBOLD-CSF coupling, as a measure of glymphatic function, over a wide age range between 36-100 years old. 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.

Dynamic 3D imaging of cerebral blood flow in awake mice using self-supervised-learning-enhanced optical coherence Doppler tomography

Cerebral blood flow (CBF) is widely used to assess brain function. However, most preclinical CBF studies have been performed under anesthesia, which confounds findings. High spatiotemporal-resolution CBF imaging of awake animals is challenging due to motion artifacts and background noise, particularly for Doppler-based flow imaging. Here, we report ultrahigh-resolution optical coherence Doppler tomography (µODT) for 3D imaging of CBF velocity (CBFv) dynamics in awake mice by developing self-supervised deep-learning for effective image denoising and motion-artifact removal. We compare cortical CBFv in awake vs. anesthetized mice and their dynamic responses in arteriolar, venular and capillary networks to acute cocaine (1 mg/kg, i.v.), a highly addictive drug associated with neurovascular toxicity. Compared with awake, isoflurane (2-2.5%) induces vasodilation and increases CBFv within 2-4 min, whereas dexmedetomidine (0.025 mg/kg, i.p.) does not change vessel diameters nor flow. Acute cocaine decreases CBFv to the same extent in dexmedetomidine and awake states, whereas decreases are larger under isoflurane, suggesting that isoflurane-induced vasodilation might have facilitated detection of cocaine-induced vasoconstriction. Awake mice after chronic cocaine show severe vasoconstriction, CBFv decreases and vascular adaptations with extended diving arteriolar/venular vessels that prioritize blood supply to deeper cortical capillaries. The 3D imaging platform we present provides a powerful tool to study dynamic changes in vessel diameters and morphology alongside CBFv networks in the brain of awake animals that can advance our understanding of the effects of drugs and disease conditions (ischemia, tumors, wound healing).

Intrinsic macroscale oscillatory modes driving long range functional connectivity in female rat brains detected by ultrafast fMRI

Spontaneous fluctuations in functional magnetic resonance imaging (fMRI) signals correlate across distant brain areas, shaping functionally relevant intrinsic networks. However, the generative mechanism of fMRI signal correlations, and in particular the link with locally-detected ultra-slow oscillations, are not fully understood. To investigate this link, we record ultrafast ultrahigh field fMRI signals (9.4 Tesla, temporal resolution = 38 milliseconds) from female rats across three anesthesia conditions. Power at frequencies extending up to 0.3 Hz is detected consistently across rat brains and is modulated by anesthesia level. Principal component analysis reveals a repertoire of modes, in which transient oscillations organize with fixed phase relationships across distinct cortical and subcortical structures. Oscillatory modes are found to vary between conditions, resonating at faster frequencies under medetomidine sedation and reducing both in number, frequency, and duration with the addition of isoflurane. Peaking in power within clear anatomical boundaries, these oscillatory modes point to an emergent systemic property. This work provides additional insight into the origin of oscillations detected in fMRI and the organizing principles underpinning spontaneous long-range functional connectivity.

Triple network activation causes tinnitus in patients with sudden sensorineural hearing loss: A model-based volume-entropy analysis

Tinnitus can be defined as the conscious perception of phantom sounds in the absence of corresponding external auditory signals. Tinnitus can develop in the setting of sudden sensorineural hearing loss (SSNHL), but the underlying mechanism is largely unknown. Using electroencephalography, we investigated differences in afferent node capacity between 15 SSNHL patients without tinnitus (NT) and 30 SSNHL patients with tinnitus (T). Where the T group showed increased afferent node capacity in regions constituting a "triple brain network" [default mode network (DMN), central executive network (CEN), and salience network (SN)], the NT group showed increased information flow in regions implicated in temporal auditory processing and noise-canceling pathways. Our results demonstrate that when all components of the triple network are activated due to sudden-onset auditory deprivation, tinnitus ensues. By contrast, auditory processing-associated and tinnitus-suppressing networks are highly activated in the NT group, to overcome the activation of the triple network and effectively suppress the generation of tinnitus.

Altered basal forebrain function during whole-brain network activity at pre- and early-plaque stages of Alzheimer's disease in TgF344-AD rats

BACKGROUND

Imbalanced synaptic transmission appears to be an early driver in Alzheimer's disease (AD) leading to brain network alterations. Early detection of altered synaptic transmission and insight into mechanisms causing early synaptic alterations would be valuable treatment strategies. This study aimed to investigate how whole-brain networks are influenced at pre- and early-plague stages of AD and if these manifestations are associated with concomitant cellular and synaptic deficits.  METHODS: To this end, we used an established AD rat model (TgF344-AD) and employed resting state functional MRI and quasi-periodic pattern (QPP) analysis, a method to detect recurrent spatiotemporal motifs of brain activity, in parallel with state-of-the-art immunohistochemistry in selected brain regions.

RESULTS

At the pre-plaque stage, QPPs in TgF344-AD rats showed decreased activity of the basal forebrain (BFB) and the default mode-like network. Histological analyses revealed increased astrocyte abundance restricted to the BFB, in the absence of amyloid plaques, tauopathy, and alterations in a number of cholinergic, gaba-ergic, and glutamatergic synapses. During the early-plaque stage, when mild amyloid-beta (Aβ) accumulation was observed in the cortex and hippocampus, QPPs in the TgF344-AD rats normalized suggesting the activation of compensatory mechanisms during this early disease progression period. Interestingly, astrogliosis observed in the BFB at the pre-plaque stage was absent at the early-plaque stage. Moreover, altered excitatory/inhibitory balance was observed in cortical regions belonging to the default mode-like network. In wild-type rats, at both time points, peak activity in the BFB preceded peak activity in other brain regions-indicating its modulatory role during QPPs. However, this pattern was eliminated in TgF344-AD suggesting that alterations in BFB-directed neuromodulation have a pronounced impact in network function in AD.

CONCLUSIONS

This study demonstrates the value of rsfMRI and advanced network analysis methods to detect early alterations in BFB function in AD, which could aid early diagnosis and intervention in AD. Restoring the global synaptic transmission, possibly by modulating astrogliosis in the BFB, might be a promising therapeutic strategy to restore brain network function and delay the onset of symptoms in AD.

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 enigma of infra-slow fluctuations in the human EEG

Spontaneous Infra-Slow Fluctuations (ISFs) of the human EEG (EEG-ISFs) were discovered 60 years ago when appropriate amplifiers for their recordings were designed. To avoid skin-related artifacts the recording of EEG-ISFs required puncturing the skin under the electrode. In the beginning of the 21st century the interest in EEG-ISFs was renewed with the appearance of commercially available DC-coupled amplified and by observation of ISFs of the blood oxygen level–dependent (BOLD) functional magnetic resonance imaging signal at a similar frequency. The independent components of irregular EEG-ISFs were shown to correlate with BOLD signals which in turn were driven by changes in arousal level measured by galvanic skin response (GSR), pupil size and HRV. There is no consensus regarding the temporal organization of EEG-ISFs: some studies emphasize the absence of peaks on EEG-ISFs spectra, some studies report prominent oscillations with frequency around 0.1 or 0.02 Hz, while some emphasize multiple discrete infraslow oscillations. No studies used parameters of EEG-ISFs as neuromarkers to discriminate psychiatric patients from healthy controls. Finally, a set of working hypotheses is suggested that must be tested in future research to solve the enigma of EEG-ISFs.

A parsimonious description of global functional brain organization in three spatiotemporal patterns

Resting-state functional magnetic resonance imaging (MRI) has yielded seemingly disparate insights into large-scale organization of the human brain. The brain's large-scale organization can be divided into two broad categories: zero-lag representations of functional connectivity structure and time-lag representations of traveling wave or propagation structure. In this study, we sought to unify observed phenomena across these two categories in the form of three low-frequency spatiotemporal patterns composed of a mixture of standing and traveling wave dynamics. We showed that a range of empirical phenomena, including functional connectivity gradients, the task-positive/task-negative anti-correlation pattern, the global signal, time-lag propagation patterns, the quasiperiodic pattern and the functional connectome network structure, are manifestations of these three spatiotemporal patterns. These patterns account for much of the global spatial structure that underlies functional connectivity analyses and unifies phenomena in resting-state functional MRI previously thought distinct.

Spatiotemporal patterns of spontaneous brain activity: a mini-review

The brain exists in a state of constant activity in the absence of any external sensory input. The spatiotemporal patterns of this spontaneous brain activity have been studied using various recording and imaging techniques. This has enabled considerable progress to be made in elucidating the cellular and network mechanisms that are involved in the observed spatiotemporal dynamics. This mini-review outlines different spatiotemporal dynamic patterns that have been identified in four commonly used modalities: electrophysiological recordings, optical imaging, functional magnetic resonance imaging, and electroencephalography. Signal sources for each modality, possible sources of the observed dynamics, and future directions are also discussed.

NoLiTiA: An Open-Source Toolbox for Non-linear Time Series Analysis

In many scientific fields including neuroscience, climatology or physics, complex relationships can be described most parsimoniously by non-linear mechanics. Despite their relevance, many neuroscientists still apply linear estimates in order to evaluate complex interactions. This is partially due to the lack of a comprehensive compilation of non-linear methods. Available packages mostly specialize in only one aspect of non-linear time-series analysis and most often require some coding proficiency to use. Here, we introduce NoLiTiA, a free open-source MATLAB toolbox for non-linear time series analysis. In comparison to other currently available non-linear packages, NoLiTiA offers (1) an implementation of a broad range of classic and recently developed methods, (2) an implementation of newly proposed spatially and time-resolved recurrence amplitude analysis and (3) an intuitive environment accessible even to users with little coding experience due to a graphical user interface and batch-editor. The core methodology derives from three distinct fields of complex systems theory, including dynamical systems theory, recurrence quantification analysis and information theory. Besides established methodology including estimation of dynamic invariants like Lyapunov exponents and entropy-based measures, such as active information storage, we include recent developments of quantifying time-resolved aperiodic oscillations. In general, the toolbox will make non-linear methods accessible to the broad neuroscientific community engaged in time series processing.

Comparison of Resting-State Functional MRI Methods for Characterizing Brain Dynamics

Resting-state functional MRI (fMRI) exhibits time-varying patterns of functional connectivity. Several different analysis approaches have been developed for examining these resting-state dynamics including sliding window connectivity (SWC), phase synchrony (PS), co-activation pattern (CAP), and quasi-periodic patterns (QPP). Each of these approaches can be used to generate patterns of activity or inter-areal coordination which vary across time. The individual frames can then be clustered to produce temporal groupings commonly referred to as “brain states.” Several recent publications have investigated brain state alterations in clinical populations, typically using a single method for quantifying frame-wise functional connectivity. This study directly compares the results of k-means clustering in conjunction with three of these resting-state dynamics methods (SWC, CAP, and PS) and quantifies the brain state dynamics across several metrics using high resolution data from the human connectome project. Additionally, these three dynamics methods are compared by examining how the brain state characterizations vary during the repeated sequences of brain states identified by a fourth dynamic analysis method, QPP. The results indicate that the SWC, PS, and CAP methods differ in the clusters and trajectories they produce. A clear illustration of these differences is given by how each one results in a very different clustering profile for the 24s sequences explicitly identified by the QPP algorithm. PS clustering is sensitive to QPPs with the mid-point of most QPP sequences grouped into the same single cluster. CAPs are also highly sensitive to QPPs, separating each phase of the QPP sequences into different sets of clusters. SWC (60s window) is less sensitive to QPPs. While the QPPs are slightly more likely to occur during specific SWC clusters, the SWC clustering does not vary during the 24s QPP sequences, the goal of this work is to improve both the practical and theoretical understanding of different resting-state dynamics methods, thereby enabling investigators to better conceptualize and implement these tools for characterizing functional brain networks.

Functional Connectivity of the Brain Across Rodents and Humans

Resting-state functional magnetic resonance imaging (rs-fMRI), which measures the spontaneous fluctuations in the blood oxygen level-dependent (BOLD) signal, is increasingly utilized for the investigation of the brain's physiological and pathological functional activity. Rodents, as a typical animal model in neuroscience, play an important role in the studies that examine the neuronal processes that underpin the spontaneous fluctuations in the BOLD signal and the functional connectivity that results. Translating this knowledge from rodents to humans requires a basic knowledge of the similarities and differences across species in terms of both the BOLD signal fluctuations and the resulting functional connectivity. This review begins by examining similarities and differences in anatomical features, acquisition parameters, and preprocessing techniques, as factors that contribute to functional connectivity. Homologous functional networks are compared across species, and aspects of the BOLD fluctuations such as the topography of the global signal and the relationship between structural and functional connectivity are examined. Time-varying features of functional connectivity, obtained by sliding windowed approaches, quasi-periodic patterns, and coactivation patterns, are compared across species. Applications demonstrating the use of rs-fMRI as a translational tool for cross-species analysis are discussed, with an emphasis on neurological and psychiatric disorders. Finally, open questions are presented to encapsulate the future direction of the field.

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.

Aberrant functional brain network dynamics in patients with functional constipation

The aberrant static functional connectivity of brain network has been widely investigated in patients with functional constipation (FCon). However, the dynamics of brain functional connectivity in FCon patients remained unknown. This study aimed to detect the brain dynamics of functional connectivity states and network topological organizations of FCon patients and investigate the correlations of the aberrant brain dynamics with symptom severity. Eighty-three FCon patients and 80 healthy subjects (HS) were included in data analysis. The spatial group independent component analysis, sliding-window approach, k-means clustering, and graph-theoretic analysis were applied to investigate the dynamic temporal properties and coupling patterns of functional connectivity states, as well as the time-variation of network topological organizations in FCon patients. Four reoccurring functional connectivity states were identified in k-means clustering analysis. Compared to HS, FCon patients manifested the lower occurrence rate and mean dwell time in the state with a complex connection between default mode network and cognitive control network, as well as the aberrant anterior insula-cortical coupling patterns in this state, which were significantly correlated with the symptom severity. The graph-theoretic analysis demonstrated that FCon patients had higher sample entropy at the nodal efficiency of anterior insula than HS. The current findings provided dynamic perspectives for understanding the brain connectome of FCon and laid the foundation for the potential treatment of FCon based on brain connectomics.

MRI magnetic compatible electrical neural interface: From materials to application

Neural electrical interfaces are important tools for local neural stimulation and recording, which potentially have wide application in the diagnosis and treatment of neural diseases, as well as in the transmission of neural activity for brain-computer interface (BCI) systems. At the same time, magnetic resonance imaging (MRI) is one of the effective and non-invasive techniques for recording whole-brain signals, providing details of brain structures and also activation pattern maps. Simultaneous recording of extracellular neural signals and MRI combines two expressions of the same neural activity and is believed to be of great importance for the understanding of brain function. However, this combination makes requests on the magnetic and electronic performance of neural interface devices. MRI-compatibility refers here to a technological approach to simultaneous MRI and electrode recording or stimulation without artifacts in imaging. Trade-offs between materials magnetic susceptibility selection and electrical function should be considered. Herein, prominent trends in selecting materials of suitable magnetic properties are analyzed and material design, function and application of neural interfaces are outlined together with the remaining challenge to fabricate MRI-compatible neural interface.

Imaging the temporal dynamics of brain states with highly sampled fMRI

The spontaneous dynamics of the brain modulate its function from moment to moment, shaping neural computation and cognition. Functional MRI (fMRI), while classically used as a tool for spatial localization, is increasingly being used to identify the temporal dynamics of brain activity. fMRI analyses focused on the temporal domain have revealed important new information about the dynamics underlying states such as arousal, attention, and sleep. Dense temporal sampling – either by using fast fMRI acquisition, or multiple repeated scan sessions within individuals – can further enrich the information present in these studies. This review focuses on recent developments in using fMRI to identify dynamics across brain states, particularly vigilance and sleep states, and the potential for highly temporally sampled fMRI to answer these questions.

An isotropic EPI database and analytical pipelines for rat brain resting-state fMRI

Resting-state functional magnetic resonance imaging (fMRI) has drastically expanded the scope of brain research by advancing our knowledge about the topologies, dynamics, and interspecies translatability of functional brain networks. Several databases have been developed and shared in accordance with recent key initiatives in the rodent fMRI community to enhance the transparency, reproducibility, and interpretability of data acquired at various sites. Despite these pioneering efforts, one notable challenge preventing efficient standardization in the field is the customary choice of anisotropic echo planar imaging (EPI) schemes with limited spatial coverage. Imaging with anisotropic resolution and/or reduced brain coverage has significant shortcomings including reduced registration accuracy and increased deviation in brain feature detection. Here we proposed a high-spatial-resolution (0.4 mm), isotropic, whole-brain EPI protocol for the rat brain using a horizontal slicing scheme that can maintain a functionally relevant repetition time (TR), avoid high gradient duty cycles, and offer unequivocal whole-brain coverage. Using this protocol, we acquired resting-state EPI fMRI data from 87 healthy rats under the widely used dexmedetomidine sedation supplemented with low-dose isoflurane on a 9.4 T MRI system. We developed an EPI template that closely approximates the Paxinos and Watson's rat brain coordinate system and demonstrated its ability to improve the accuracy of group-level approaches and streamline fMRI data pre-processing. Using this database, we employed a multi-scale dictionary-learning approach to identify reliable spatiotemporal features representing rat brain intrinsic activity. Subsequently, we performed k-means clustering on those features to obtain spatially discrete, functional regions of interest (ROIs). Using Euclidean-based hierarchical clustering and modularity-based partitioning, we identified the topological organizations of the rat brain. Additionally, the identified group-level FC network appeared robust across strains and sexes. The "triple-network" commonly adapted in human fMRI were resembled in the rat brain. Through this work, we disseminate raw and pre-processed isotropic EPI data, a rat brain EPI template, as well as identified functional ROIs and networks in standardized rat brain coordinates. We also make our analytical pipelines and scripts publicly available, with the hope of facilitating rat brain resting-state fMRI study standardization.

Contribution of animal models toward understanding resting state functional connectivity

Functional connectivity, which reflects the spatial and temporal organization of intrinsic activity throughout the brain, is one of the most studied measures in human neuroimaging research. The noninvasive acquisition of resting state functional magnetic resonance imaging (rs-fMRI) allows the characterization of features designated as functional networks, functional connectivity gradients, and time-varying activity patterns that provide insight into the intrinsic functional organization of the brain and potential alterations related to brain dysfunction. Functional connectivity, hence, captures dimensions of the brain's activity that have enormous potential for both clinical and preclinical research. However, the mechanisms underlying functional connectivity have yet to be fully characterized, hindering interpretation of rs-fMRI studies. As in other branches of neuroscience, the identification of the neurophysiological processes that contribute to functional connectivity largely depends on research conducted on laboratory animals, which provide a platform where specific, multi-dimensional investigations that involve invasive measurements can be carried out. These highly controlled experiments facilitate the interpretation of the temporal correlations observed across the brain. Indeed, information obtained from animal experimentation to date is the basis for our current understanding of the underlying basis for functional brain connectivity. This review presents a compendium of some of the most critical advances in the field based on the efforts made by the animal neuroimaging community.

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.

Behavioural relevance of spontaneous, transient brain network interactions in fMRI

How spontaneously fluctuating functional magnetic resonance imaging (fMRI) signals in different brain regions relate to behaviour has been an open question for decades. Correlations in these signals, known as functional connectivity, can be averaged over several minutes of data to provide a stable representation of the functional network architecture for an individual. However, associations between these stable features and behavioural traits have been shown to be dominated by individual differences in anatomy. Here, using kernel learning tools, we propose methods to assess and compare the relation between time-varying functional connectivity, time-averaged functional connectivity, structural brain data, and non-imaging subject behavioural traits. We applied these methods to Human Connectome Project resting-state fMRI data to show that time-varying fMRI functional connectivity, detected at time-scales of a few seconds, has associations with some behavioural traits that are not dominated by anatomy. Despite time-averaged functional connectivity accounting for the largest proportion of variability in the fMRI signal between individuals, we found that some aspects of intelligence could only be explained by time-varying functional connectivity. The finding that time-varying fMRI functional connectivity has a unique relationship to population behavioural variability suggests that it might reflect transient neuronal communication fluctuating around a stable neural architecture.

Propagating patterns of intrinsic activity along macroscale gradients coordinate functional connections across the whole brain

The intrinsic activity of the human brain, observed with resting-state fMRI (rsfMRI) and functional connectivity, exhibits macroscale spatial organization such as functional networks and gradients. Dynamic analysis techniques have shown that functional connectivity is a mere summary of time-varying patterns with distinct spatial and temporal characteristics. A better understanding of these patterns might provide insight into aspects of the brain's intrinsic activity that cannot be inferred by functional connectivity or the spatial maps derived from it, such as functional networks and gradients. Here, we describe three spatiotemporal patterns of coordinated activity across the whole brain obtained by averaging similar ~20-second-long segments of rsfMRI timeseries. In each of these patterns, activity propagates along a particular macroscale functional gradient, simultaneously across the cerebral cortex and in most other brain regions. In some regions, like the thalamus, the propagation suggests previously-undescribed gradients. The coordinated activity across areas is consistent with known tract-based connections, and nuanced differences in the timing of peak activity between regions point to plausible driving mechanisms. The magnitude of correlation within and particularly between functional networks is remarkably diminished when these patterns are regressed from the rsfMRI timeseries, a quantitative demonstration of their significant role in functional connectivity. Taken together, our results suggest that a few recurring patterns of propagating intrinsic activity along macroscale gradients give rise to and coordinate functional connections across the whole brain.

Dynamic time warping outperforms Pearson correlation in detecting atypical functional connectivity in autism spectrum disorders

Resting state fMRI (rsfMRI) is frequently used to study brain function, including in clinical populations. Similarity of blood-oxygen-level-dependent (BOLD) fluctuations during rsfMRI between brain regions is thought to reflect intrinsic functional connectivity (FC), potentially due to history of coactivation. To quantify similarity, studies have almost exclusively relied on Pearson correlation, which assumes linearity and can therefore underestimate FC if the hemodynamic response function differs regionally or there is BOLD signal lag between timeseries. Here we show in three cohorts of children, adolescents and adults, with and without autism spectrum disorders (ASDs), that measuring the similarity of BOLD signal fluctuations using non-linear dynamic time warping (DTW) is more robust to global signal regression (GSR), has higher test-retest reliability and is more sensitive to task-related changes in FC. Additionally, when comparing FC between individuals with ASDs and typical controls, more group differences are detected using DTW. DTW estimates are also more related to ASD symptom severity and executive function, while Pearson correlation estimates of FC are more strongly associated with respiration during rsfMRI. Together these findings suggest that non-linear methods such as DTW improve estimation of resting state FC, particularly when studying clinical populations whose hemodynamics or neurovascular coupling may be altered compared to typical controls.

Dynamic Brain Connectivity in Resting State Functional MR Imaging

Dynamic functional connectivity adds another dimension to resting-state functional MR imaging analysis. In recent years, dynamic functional connectivity has been increasingly used in resting-state functional MR imaging, and several studies have demonstrated that dynamic functional connectivity patterns correlate with different physiologic and pathologic brain states. In fact, evidence suggests that dynamic functional connectivity is a more sensitive marker than static functional connectivity; therefore, it might be a promising tool to add to clinical functional neuroimaging. This article provides a broad overview of dynamic functional connectivity and reviews its general principles, techniques, and potential clinical applications.

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.

Relationship Between Basic Properties of BOLD Fluctuations and Calculated Metrics of Complexity in the Human Connectome Project

Resting state functional MRI (rs-fMRI) creates a rich four-dimensional data set that can be analyzed in a variety of ways. As more researchers come to view the brain as a complex dynamical system, tools are increasingly being drawn from other fields to characterize the complexity of the brain's activity. However, given that the signal measured with rs-fMRI arises from the hemodynamic response to neural activity, the extent to which complexity metrics reflect neural complexity as compared to signal properties related to image quality remains unknown. To provide some insight into this question, correlation dimension, approximate entropy and Lyapunov exponent were calculated for different rs-fMRI scans from the same subject to examine their reliability. The metrics of complexity were then compared to several properties of the rs-fMRI signal from each brain area to determine if basic signal features could explain differences in the complexity metrics. Differences in complexity across brain areas were highly reliable and were closely linked to differences in the frequency profiles of the rs-fMRI signal. The spatial distributions of the complexity and frequency metrics suggest that they are both influenced by location-dependent signal properties that can obscure changes related to neural activity.

Ultra-slow Oscillations in fMRI and Resting-State Connectivity: Neuronal and Vascular Contributions and Technical Confounds

Ultra-slow, ∼0.1-Hz variations in the oxygenation level of brain blood are widely used as an fMRI-based surrogate of "resting-state" neuronal activity. The temporal correlations among these fluctuations across the brain are interpreted as "functional connections" for maps and neurological diagnostics. Ultra-slow variations in oxygenation follow a cascade. First, they closely track changes in arteriole diameter. Second, interpretable functional connections arise when the ultra-slow changes in amplitude of γ-band neuronal oscillations, which are shared across even far-flung but synaptically connected brain regions, entrain the ∼0.1-Hz vasomotor oscillation in diameter of local arterioles. Significant confounds to estimates of functional connectivity arise from residual vasomotor activity as well as arteriole dynamics driven by self-generated movements and subcortical common modulatory inputs. Last, methodological limitations of fMRI can lead to spurious functional connections. The neuronal generator of ultra-slow variations in γ-band amplitude, including that associated with self-generated movements, remains an open issue.

Bring the Noise: Reconceptualizing Spontaneous Neural Activity

Definitions of what constitutes the 'signal of interest' in neuroscience can be controversial, due in part to continuously evolving notions regarding the significance of spontaneous neural activity. This review highlights how the challenge of separating brain signal from noise has led to new conceptualizations of brain functional organization at both the micro- and macroscopic level. Recent debates in the functional neuroimaging community surrounding artifact removal processes have revived earlier discussions surrounding how to appropriately isolate and measure neuronal signals against a background of noise from various sources. Insights from electrophysiological studies and computational modeling can inform current theory and data analytic practices in human functional neuroimaging, given that signal and noise may be inextricably linked in the brain.

Single-scale time-dependent window-sizes in sliding-window dynamic functional connectivity analysis: A validation study

During the past ten years, dynamic functional connectivity (FC) has been extensively studied using the sliding-window method. A fixed window-size is usually selected heuristically, since no consensus exists yet on choice of the optimal window-size. Furthermore, without a known ground-truth, the validity of the computed dynamic FC remains unclear and questionable. In this study, we computed single-scale time-dependent (SSTD) window-sizes for the sliding-window method. SSTD window-sizes were based on the frequency content at every time point of a time series and were computed without any prior information. Therefore, they were time-dependent and data-driven. Using simulated sinusoidal time series with frequency shifts, we demonstrated that SSTD window-sizes captured the time-dependent period (inverse of frequency) information at every time point. We further validated the dynamic FC values computed with SSTD window-sizes with both a classification analysis using fMRI data with a low sampling rate and a regression analysis using fMRI data with a high sampling rate. Specifically, we achieved both a higher classification accuracy in predicting cognitive impairment status in fighters and a larger explained behavioral variance in healthy young adults when using dynamic FC matrices computed with SSTD window-sizes as features, as compared to using dynamic FC matrices computed with the conventional fixed window-sizes. Overall, our study computed and validated SSTD window-sizes in the sliding-window method for dynamic FC analysis. Our results demonstrate that dynamic FC matrices computed with SSTD window-sizes can capture more temporal dynamic information related to behavior and cognitive function.

Aberrant static and dynamic functional patterns of frontoparietal control network in antipsychotic-naïve first-episode psychosis subjects

Psychotic disorders are disabling clinical syndromes characterized by widespread alterations in cortical information processing. Disruption of frontoparietal network (FPN) connectivity has emerged as a common footprint across the psychosis spectrum. Our goal was to characterize the static and dynamic resting‐state functional connectivity (FC) of the FPN in antipsychotic‐naïve first‐episode psychosis (FEP) subjects. We compared the static FC of the FPN in 40 FEP and 40 healthy control (HC) subjects, matched on age, sex, and socioeconomic status. To study the dynamic FC, we measured quasiperiodic patterns (QPPs) that consist of infraslow spatioemporal patterns embedded in the blood oxygen level‐dependent signal that repeats over time, exhibiting alternation of high and low activity. Relative to HC, we found functional hypoconnectivity between the right middle frontal gyrus and the right middle temporal gyrus, as well as the left inferior temporal gyrus and the left inferior parietal gyrus in FEP (p < .05, false discovery rate corrected). The correlation of the QPP with all functional scans was significantly stronger for FEP compared to HC, suggesting a greater impact of the QPPs to intrinsic brain activity in psychotic population. Regressing the QPP from the functional scans erased all significant group differences in static FC, suggesting that abnormal connectivity in FEP could result from altered QPP. Our study supports that alterations of cortical information processing are not a function of psychotic chronicity or antipsychotic medication exposure and may be regarded as trait specific. In addition, static connectivity abnormality may be partly related to altered brain network temporal dynamics.

Deep Temporal Organization of fMRI Phase Synchrony Modes Promotes Large-Scale Disconnection in Schizophrenia

Itinerant dynamics of the brain generates transient and recurrent spatiotemporal patterns in neuroimaging data. Characterizing metastable functional connectivity (FC) - particularly at rest and using functional magnetic resonance imaging (fMRI) - has shaped the field of dynamic functional connectivity (DFC). Mainstream DFC research relies on (sliding window) correlations to identify recurrent FC patterns. Recently, functional relevance of the instantaneous phase synchrony (IPS) of fMRI signals has been revealed using imaging studies and computational models. In the present paper, we identify the repertoire of whole-brain inter-network IPS states at rest. Moreover, we uncover a hierarchy in the temporal organization of IPS modes. We hypothesize that connectivity disorder in schizophrenia (SZ) is related to the (deep) temporal arrangement of large-scale IPS modes. Hence, we analyze resting-state fMRI data from 68 healthy controls (HC) and 51 SZ patients. Seven resting-state networks (and their sub-components) are identified using spatial independent component analysis. IPS is computed between subject-specific network time courses, using analytic signals. The resultant phase coupling patterns, across time and subjects, are clustered into eight IPS states. Statistical tests show that the relative expression and mean lifetime of certain IPS states have been altered in SZ. Namely, patients spend (45%) less time in a globally coherent state and a subcortical-centered state, and (40%) more time in states reflecting anticoupling within the cognitive control network, compared to the HC. Moreover, the transition profile (between states) reveals a deep temporal structure, shaping two metastates with distinct phase synchrony profiles. A metastate is a collection of states such that within-metastate transitions are more probable than across. Remarkably, metastate occupation balance is altered in SZ, in favor of the less synchronous metastate that promotes disconnection across networks. Furthermore, the trajectory of IPS patterns is less efficient, less smooth, and more restricted in SZ subjects, compared to the HC. Finally, a regression analysis confirms the diagnostic value of the defined IPS measures for SZ identification, highlighting the distinctive role of metastate proportion. Our results suggest that the proposed IPS features may be used for classification studies and for characterizing phase synchrony modes in other (clinical) populations.

Questions and controversies in the study of time-varying functional connectivity in resting fMRI

The brain is a complex, multiscale dynamical system composed of many interacting regions. Knowledge of the spatiotemporal organization of these interactions is critical for establishing a solid understanding of the brain's functional architecture and the relationship between neural dynamics and cognition in health and disease. The possibility of studying these dynamics through careful analysis of neuroimaging data has catalyzed substantial interest in methods that estimate time-resolved fluctuations in functional connectivity (often referred to as "dynamic" or time-varying functional connectivity; TVFC). At the same time, debates have emerged regarding the application of TVFC analyses to resting fMRI data, and about the statistical validity, physiological origins, and cognitive and behavioral relevance of resting TVFC. These and other unresolved issues complicate interpretation of resting TVFC findings and limit the insights that can be gained from this promising new research area. This article brings together scientists with a variety of perspectives on resting TVFC to review the current literature in light of these issues. We introduce core concepts, define key terms, summarize controversies and open questions, and present a forward-looking perspective on how resting TVFC analyses can be rigorously and productively applied to investigate a wide range of questions in cognitive and systems neuroscience.