A unified approach for characterizing static/dynamic connectivity frequency profiles using filter banks

A survey of brain functional network extraction methods using fMRI data

Functional network (FN) analyses play a pivotal role in uncovering insights into brain function and understanding the pathophysiology of various brain disorders. This paper focuses on classical and advanced methods for deriving brain FNs from functional magnetic resonance imaging (fMRI) data. We systematically review their foundational principles, advantages, shortcomings, and interrelations, encompassing both static and dynamic FN extraction approaches. In the context of static FN extraction, we present hypothesis-driven methods such as region of interest (ROI)-based approaches as well as data-driven methods including matrix decomposition, clustering, and deep learning. For dynamic FN extraction, both window-based and windowless methods are surveyed with respect to the estimation of time-varying FN and the subsequent computation of FN states. We also discuss the scope of application of the various methods and avenues for future improvements.

The overlap across psychotic disorders: A functional network connectivity analysis

Functional network connectivity (FNC) has previously been shown to distinguish patient groups from healthy controls (HC). However, the overlap across psychiatric disorders such as schizophrenia (SZ), bipolar (BP), and schizoaffective disorder (SAD) is not evident yet. This study focuses on studying the overlap across these three psychotic disorders in both dynamic and static FNC (dFNC/sFNC). We used resting-state fMRI, demographics, and clinical information from the Bipolar-Schizophrenia Network on Intermediate Phenotypes cohort (BSNIP). The data includes three groups of patients with schizophrenia (SZ, N = 181), bipolar (BP, N = 163), and schizoaffective (SAD, N = 130) and HC (N = 238) groups. After estimating each individual's dFNC, we group them into three distinct states. We evaluated two dFNC features, including occupancy rate (OCR) and distance travelled over time. Finally, the extracted features, including both sFNC and dFNC, are tested statistically across patients and HC groups. In addition, we explored the link between the clinical scores and the extracted features. We evaluated the connectivity patterns and their overlap among SZ, BP, and SAD disorders (false discovery rate or FDR corrected p < 0.05). Results showed dFNC captured unique information about overlap across disorders where all disorder groups showed similar pattern of activity in state 2. Moreover, the results showed similar patterns between SZ and SAD in state 1 which was different than BP. Finally, the distance travelled feature of SZ (average R = 0.245, p < 0.01) and combined distance travelled from all disorders was predictive of the PANSS symptoms scores (average R = 0.147, p < 0.01).

Studying time-resolved functional connectivity via communication theory: on the complementary nature of phase synchronization and sliding window Pearson correlation

Time-resolved functional connectivity (trFC) assesses the time-resolved coupling between brain regions using functional magnetic resonance imaging (fMRI) data. This study aims to compare two techniques used to estimate trFC, to investigate their similarities and differences when applied to fMRI data. These techniques are the sliding window Pearson correlation (SWPC), an amplitude-based approach, and phase synchronization (PS), a phase-based technique. To accomplish our objective, we used resting-state fMRI data from the Human Connectome Project (HCP) with 827 subjects (repetition time: 0.7s) and the Function Biomedical Informatics Research Network (fBIRN) with 311 subjects (repetition time: 2s), which included 151 schizophrenia patients and 160 controls. Our simulations reveal distinct strengths in two connectivity methods: SWPC captures high-magnitude, low-frequency connectivity, while PS detects low-magnitude, high-frequency connectivity. Stronger correlations between SWPC and PS align with pronounced fMRI oscillations. For fMRI data, higher correlations between SWPC and PS occur with matched frequencies and smaller SWPC window sizes (~30s), but larger windows (~88s) sacrifice clinically relevant information. Both methods identify a schizophrenia-associated brain network state but show different patterns: SWPC highlights low anti-correlations between visual, subcortical, auditory, and sensory-motor networks, while PS shows reduced positive synchronization among these networks. In sum, our findings underscore the complementary nature of SWPC and PS, elucidating their respective strengths and limitations without implying the superiority of one over the other.

Inflammation affects dynamic functional network connectivity pattern changes via plasma NFL in cognitive impairment patients

BACKGROUND

Plasma neurofilament light chain (NFL) is a biomarker of inflammation and neurodegenerative diseases such as Alzheimer's disease (AD). However, the underlying neural mechanisms by which NFL affects cognitive function remain unclear. In this study, we investigated the effects of inflammation on cognitive integrity in patients with cognitive impairment through the functional interaction of plasma NFL with large-scale brain networks.

METHODS

This study included 29 cognitively normal, 55 LowNFL patients, and 55 HighNFL patients. Group independent component analysis (ICA) was applied to the resting-state fMRI data, and 40 independent components (IC) were extracted for the whole brain. Next, the dynamic functional network connectivity (dFNC) of each subject was estimated using the sliding-window method and k-means clustering, and five dynamic functional states were identified. Finally, we applied mediation analysis to investigate the relationship between plasma NFL and dFNC indicators and cognitive scales.

RESULTS

The present study explored the dynamics of whole-brain FNC in controls and LowNFL and HighNFL patients and highlighted the temporal properties of dFNC states in relation to psychological scales. A potential mechanism for the association between dFNC indicators and NFL levels in cognitively impaired patients.

CONCLUSIONS

Our findings suggested the decreased ability of information processing and communication in the HighNFL group, which helps us to understand their abnormal cognitive functions clinically. Characteristic changes in the inflammation-coupled dynamic brain network may provide alternative biomarkers for the assessment of disease severity in cognitive impairment patients.

Dynamic functional connectivity in schizophrenia and bipolar disorder: A review of the evidence and associations with psychopathological features

Alterations of functional network connectivity have been implicated in the pathophysiology of schizophrenia (SCZ) and bipolar disorder (BD). Recent studies also suggest that the temporal dynamics of functional connectivity (dFC) can be altered in these disorders. Here, we summarized the existing literature on dFC in SCZ and BD, and their association with psychopathological and cognitive features. We systematically searched PubMed, Web of Science, and Scopus for studies investigating dFC in SCZ and BD and identified 77 studies. Our findings support a general model of dysconnectivity of dFC in SCZ, whereas a heterogeneous picture arose in BD. Although dFC alterations are more severe and widespread in SCZ compared to BD, dysfunctions of a triple network system underlying goal-directed behavior and sensory-motor networks were present in both disorders. Furthermore, in SCZ, positive and negative symptoms were associated with abnormal dFC. Implications for understanding the pathophysiology of disorders, the role of neurotransmitters, and treatments on dFC are discussed. The lack of standards for dFC metrics, replication studies, and the use of small samples represent major limitations for the field.

Alterations in grey matter structure linked to frequency-specific cortico-subcortical connectivity in schizophrenia via multimodal data fusion

Schizophrenia (SZ) is a complex psychiatric disorder that is currently defined by symptomatic and behavioral, rather than biological, criteria. Neuroimaging is an appealing avenue for SZ biomarker development, as several neuroimaging-based studies comparing individuals with SZ to healthy controls (HC) have shown measurable group differences in brain structure, as well as functional brain alterations in both static and dynamic functional network connectivity (sFNC and dFNC, respectively). The recently proposed filter-banked connectivity (FBC) method extends the standard dFNC sliding-window approach to estimate FNC within an arbitrary number of distinct frequency bands. The initial implementation used a set of filters spanning the full connectivity spectral range, providing a unified approach to examine both sFNC and dFNC in a single analysis. Initial FBC results found that individuals with SZ spend more time in a less structured, more disconnected low-frequency (i.e., static) FNC state than HC, as well as preferential SZ occupancy in high-frequency connectivity states, suggesting a frequency-specific component underpinning the functional dysconnectivity observed in SZ. Building on these findings, we sought to link such frequency-specific patterns of FNC to covarying data-driven structural brain networks in the context of SZ. Specifically, we employ a multi-set canonical correlation analysis + joint independent components analysis (mCCA + jICA) data fusion framework to study the connection between grey matter volume (GMV) maps and FBC states across the full connectivity frequency spectrum. Our multimodal analysis identified two joint sources that captured co-varying patterns of frequency-specific functional connectivity and alterations in GMV with significant group differences in loading parameters between the SZ group and HC. The first joint source linked frequency-modulated connections between the subcortical and sensorimotor networks and GMV alterations in the frontal and temporal lobes, while the second joint source identified a relationship between low-frequency cerebellar-sensorimotor connectivity and structural changes in both the cerebellum and motor cortex. Together, these results show a strong connection between cortico-subcortical functional connectivity at both high and low frequencies and alterations in cortical GMV that may be relevant to the pathogenesis and pathophysiology of SZ.

Phase and amplitude, two sides of functional connectivity

The dynamics of the human brain can be captured by estimating time-resolved functional network connectivity (trFNC). The most used method for estimating trFNC is sliding window Pearson correlation (SWPC). Methods based on instantaneous phase synchrony, which uses phase information for estimating trFNC are being increasingly used. These two approaches are similar under specific assumptions. Prior works have focused on which of these approaches is the best. Some works argue that SWPC can capture amplitude information and therefore we believe that instantaneous phase synchrony methods and SWPC capture different aspects of connectivity since phase synchrony methods work with the phase of the signal. Here we show that these two approaches result in different time-resolved information and therefore should be viewed as complimentary views of connectivity.

A method for estimating and characterizing explicitly nonlinear dynamic functional network connectivity in resting-state fMRI data

The past 10 years have seen an explosion of approaches that focus on the study of time-resolved change in functional connectivity (FC). FC characterization among networks at a whole-brain level is frequently termed functional network connectivity (FNC). Time-resolved or dynamic functional network connectivity (dFNC) focuses on the estimation of transient, recurring, whole-brain patterns of FNC. While most approaches in this area have attempted to capture dynamic linear correlation, we are particularly interested in whether explicitly nonlinear relationships, above and beyond linear, are present and contain unique information. This study thus proposes an approach to assess explicitly nonlinear dynamic functional network connectivity (EN dFNC) derived from the relationship among independent component analysis time courses. Linear relationships were removed at each time point to evaluate, typically ignored, explicitly nonlinear dFNC using normalized mutual information (NMI). Simulations showed the proposed method estimated explicitly nonlinearity over time, even within relatively short windows of data. We then, applied our approach on 151 schizophrenia patients, and 163 healthy controls fMRI data and found three unique, highly structured, mostly long-range, functional states that also showed significant group differences. In particular, explicitly nonlinear relationships tend to be more widespread than linear ones. Results also highlighted a state with long range connections to the visual domain, which were significantly reduced in schizophrenia. Overall, this work suggests that quantifying EN dFNC may provide a complementary and potentially valuable tool for studying brain function by exposing relevant variation that is typically ignored.

A method for capturing dynamic spectral coupling in resting fMRI reveals domain-specific patterns in schizophrenia

Introduction

Resting-state functional magnetic resonance imaging (rs-fMRI) is a powerful tool for assessing functional brain connectivity. Recent studies have focused on shorter-term connectivity and dynamics in the resting state. However, most of the prior work evaluates changes in time-series correlations. In this study, we propose a framework that focuses on time-resolved spectral coupling (assessed via the correlation between power spectra of the windowed time courses) among different brain circuits determined via independent component analysis (ICA).

Methods

Motivated by earlier work suggesting significant spectral differences in people with schizophrenia, we developed an approach to evaluate time-resolved spectral coupling (trSC). To do this, we first calculated the correlation between the power spectra of windowed time-courses pairs of brain components. Then, we subgrouped each correlation map into four subgroups based on the connectivity strength utilizing quartiles and clustering techniques. Lastly, we examined clinical group differences by regression analysis for each averaged count and average cluster size matrices in each quartile. We evaluated the method by applying it to resting-state data collected from 151 (114 males, 37 females) people with schizophrenia (SZ) and 163 (117 males, 46 females) healthy controls (HC).

Results

Our proposed approach enables us to observe the change of connectivity strength within each quartile for different subgroups. People with schizophrenia showed highly modularized and significant differences in multiple network domains, whereas males and females showed less modular differences. Both cell count and average cluster size analysis for subgroups indicate a higher connectivity rate in the fourth quartile for the visual network in the control group. This indicates increased trSC in visual networks in the controls. In other words, this shows that the visual networks in people with schizophrenia have less mutually consistent spectra. It is also the case that the visual networks are less spectrally correlated on short timescales with networks of all other functional domains.

Conclusions

The results of this study reveal significant differences in the degree to which spectral power profiles are coupled over time. Importantly, there are significant but distinct differences both between males and females and between people with schizophrenia and controls. We observed a more significant coupling rate in the visual network for the healthy controls and males in the upper quartile. Fluctuations over time are complex, and focusing on only time-resolved coupling among time-courses is likely to miss important information. Also, people with schizophrenia are known to have impairments in visual processing but the underlying reasons for the impairment are still unknown. Therefore, the trSC approach can be a useful tool to explore the reasons for the impairments.

Edge-centric analysis of time-varying functional brain networks with applications in autism spectrum disorder

The interaction between brain regions changes over time, which can be characterized using time-varying functional connectivity (tvFC). The common approach to estimate tvFC uses sliding windows and offers limited temporal resolution. An alternative method is to use the recently proposed edge-centric approach, which enables the tracking of moment-to-moment changes in co-fluctuation patterns between pairs of brain regions. Here, we first examined the dynamic features of edge time series and compared them to those in the sliding window tvFC (sw-tvFC). Then, we used edge time series to compare subjects with autism spectrum disorder (ASD) and healthy controls (CN). Our results indicate that relative to sw-tvFC, edge time series captured rapid and bursty network-level fluctuations that synchronize across subjects during movie-watching. The results from the second part of the study suggested that the magnitude of peak amplitude in the collective co-fluctuations of brain regions (estimated as root sum square (RSS) of edge time series) is similar in CN and ASD. However, the trough-to-trough duration in RSS signal is greater in ASD, compared to CN. Furthermore, an edge-wise comparison of high-amplitude co-fluctuations showed that the within-network edges exhibited greater magnitude fluctuations in CN. Our findings suggest that high-amplitude co-fluctuations captured by edge time series provide details about the disruption of functional brain dynamics that could potentially be used in developing new biomarkers of mental disorders.

A Unified Framework for Modularizing and Comparing Time-Resolved Functional Connectivity Methods

Functional connectivity is a widely used measure for finding the relationships between functional entities of the brain. Recently, more focus has been put on the methods that aim to estimate these relationships in a time-resolved fashion. However, the similarities and differences between these methods are not always clear and can result in unfair and incorrect comparisons. Here, we present a framework that provides a unified, systematic view for some of the more well-known methods. Using the proposed unified framework, we explain different methodologies using a unified language and show how they are similar and different conceptually. We give examples how this framework exposes important assumptions made by various methods, which can help clarify differences in results and facilitate reproducibility. We also show how such a framework will enable us to develop methods that improve upon previous methods.

Multi-spatial-scale dynamic interactions between functional sources reveal sex-specific changes in schizophrenia

We introduce an extension of independent component analysis (ICA), called multiscale ICA, and design an approach to capture dynamic functional source interactions within and between multiple spatial scales. Multiscale ICA estimates functional sources at multiple spatial scales without imposing direct constraints on the size of functional sources, overcomes the limitation of using fixed anatomical locations, and eliminates the need for model-order selection in ICA analysis. We leveraged this approach to study sex-specific and sex-common connectivity patterns in schizophrenia. Results show dynamic reconfiguration and interaction within and between multi-spatial scales. Sex-specific differences occur (a) within the subcortical domain, (b) between the somatomotor and cerebellum domains, and (c) between the temporal domain and several others, including the subcortical, visual, and default mode domains. Most of the sex-specific differences belong to between-spatial-scale functional interactions and are associated with a dynamic state with strong functional interactions between the visual, somatomotor, and temporal domains and their anticorrelation patterns with the rest of the brain. We observed significant correlations between multi-spatial-scale functional interactions and symptom scores, highlighting the importance of multiscale analyses to identify potential biomarkers for schizophrenia. As such, we recommend such analyses as an important option for future functional connectivity studies.

Multiframe Evolving Dynamic Functional Connectivity (EVOdFNC): A Method for Constructing and Investigating Functional Brain Motifs

The study of brain network connectivity as a time-varying property began relatively recently and, to date, has remained primarily concerned with capturing a handful of discrete static states that characterize connectivity as measured on a timescale shorter than that of the full scan. Capturing group-level representations of temporally evolving patterns of connectivity is a challenging and important next step in fully leveraging the information available in large resting state functional magnetic resonance imaging (rs-fMRI) studies. We introduce a flexible, extensible data-driven framework for the stable identification of group-level multiframe (movie-style) dynamic functional network connectivity (dFNC) states. Our approach employs uniform manifold approximation and embedding (UMAP) to produce a continuity-preserving planar embedding of high-dimensional time-varying measurements of whole-brain functional network connectivity. Planar linear exemplars summarizing dominant dynamic trends across the population are computed from local linear approximations to the two-dimensional 2D embedded trajectories. A high-dimensional representation of each 2D exemplar segment is obtained by averaging the dFNC observations corresponding to the n planar nearest neighbors of τ evenly spaced points along the 2D line segment representation (where n is the UMAP number-of-neighbors parameter and τ is the temporal duration of trajectory segments being approximated). Each of the 2D exemplars thus “lifts” to a multiframe high-dimensional dFNC trajectory of length τ. The collection of high-dimensional temporally evolving dFNC representations (EVOdFNCs) derived in this manner are employed as dynamic basis objects with which to characterize observed high-dimensional dFNC trajectories, which are then expressed as weighted combination of these basis objects. Our approach yields new insights into anomalous patterns of fluidly varying whole-brain connectivity that are significantly associated with schizophrenia as a broad diagnosis as well as with certain symptoms of this serious disorder. Importantly, we show that relative to conventional hidden Markov modeling with single-frame unvarying dFNC summary states, EVOdFNCs are more sensitive to positive symptoms of schizophrenia including hallucinations and delusions, suggesting that a more dynamic characterization is needed to help illuminate such a complex brain disorder.

Dynamic expression of brain functional systems disclosed by fine-scale analysis of edge time series

Functional connectivity (FC) describes the statistical dependence between neuronal populations or brain regions in resting-state fMRI studies and is commonly estimated as the Pearson correlation of time courses. Clustering or community detection reveals densely coupled sets of regions constituting resting-state networks or functional systems. These systems manifest most clearly when FC is sampled over longer epochs but appear to fluctuate on shorter timescales. Here, we propose a new approach to reveal temporal fluctuations in neuronal time series. Unwrapping FC signal correlations yields pairwise co-fluctuation time series, one for each node pair or edge, and allows tracking of fine-scale dynamics across the network. Co-fluctuations partition the network, at each time step, into exactly two communities. Sampled over time, the overlay of these bipartitions, a binary decomposition of the original time series, very closely approximates functional connectivity. Bipartitions exhibit characteristic spatiotemporal patterns that are reproducible across participants and imaging runs, capture individual differences, and disclose fine-scale temporal expression of functional systems. Our findings document that functional systems appear transiently and intermittently, and that FC results from the overlay of many variable instances of system expression. Potential applications of this decomposition of functional connectivity into a set of binary patterns are discussed.

Brain Density Clustering Analysis: A New Approach to Brain Functional Dynamics

BACKGROUND

A number of studies in recent years have explored whole-brain dynamic connectivity using pairwise approaches. There has been less focus on trying to analyze brain dynamics in higher dimensions over time.

METHODS

We introduce a new approach that analyzes time series trajectories to identify high traffic nodes in a high dimensional space. First, functional magnetic resonance imaging (fMRI) data are decomposed using spatial ICA to a set of maps and their associated time series. Next, density is calculated for each time point and high-density points are clustered to identify a small set of high traffic nodes. We validated our method using simulations and then implemented it on a real data set.

RESULTS

We present a novel approach that captures dynamics within a high dimensional space and also does not use any windowing in contrast to many existing approaches. The approach enables one to characterize and study the time series in a potentially high dimensional space, rather than looking at each component pair separately. Our results show that schizophrenia patients have a lower dynamism compared to healthy controls. In addition, we find patients spend more time in nodes associated with the default mode network and less time in components strongly correlated with auditory and sensorimotor regions. Interestingly, we also found that subjects oscillate between state pairs that show opposite spatial maps, suggesting an oscillatory pattern.

CONCLUSION

Our proposed method provides a novel approach to analyze the data in its native high dimensional space and can possibly provide new information that is undetectable using other methods.

Alzheimer's Disease Projection From Normal to Mild Dementia Reflected in Functional Network Connectivity: A Longitudinal Study

Background

Alzheimer's disease (AD) is the most common age-related problem and progresses in different stages, including mild cognitive impairment (early stage), mild dementia (middle-stage), and severe dementia (late-stage). Recent studies showed changes in functional network connectivity obtained from resting-state functional magnetic resonance imaging (rs-fMRI) during the transition from healthy aging to AD. By assuming that the brain interaction is static during the scanning time, most prior studies are focused on static functional or functional network connectivity (sFNC). Dynamic functional network connectivity (dFNC) explores temporal patterns of functional connectivity and provides additional information to its static counterpart.

Method

We used longitudinal rs-fMRI from 1385 scans (from 910 subjects) at different stages of AD (from normal to very mild AD or vmAD). We used group-independent component analysis (group-ICA) and extracted 53 maximally independent components (ICs) for the whole brain. Next, we used a sliding-window approach to estimate dFNC from the extracted 53 ICs, then group them into 3 different brain states using a clustering method. Then, we estimated a hidden Markov model (HMM) and the occupancy rate (OCR) for each subject. Finally, we investigated the link between the clinical rate of each subject with state-specific FNC, OCR, and HMM.

Results

All states showed significant disruption during progression normal brain to vmAD one. Specifically, we found that subcortical network, auditory network, visual network, sensorimotor network, and cerebellar network connectivity decrease in vmAD compared with those of a healthy brain. We also found reorganized patterns (i.e., both increases and decreases) in the cognitive control network and default mode network connectivity by progression from normal to mild dementia. Similarly, we found a reorganized pattern of between-network connectivity when the brain transits from normal to mild dementia. However, the connectivity between visual and sensorimotor network connectivity decreases in vmAD compared with that of a healthy brain. Finally, we found a normal brain spends more time in a state with higher connectivity between visual and sensorimotor networks.

Conclusion

Our results showed the temporal and spatial pattern of whole-brain FNC differentiates AD form healthy control and suggested substantial disruptions across multiple dynamic states. In more detail, our results suggested that the sensory network is affected more than other brain network, and default mode network is one of the last brain networks get affected by AD In addition, abnormal patterns of whole-brain dFNC were identified in the early stage of AD, and some abnormalities were correlated with the clinical score.

Tools of the trade: estimating time-varying connectivity patterns from fMRI data

Given the dynamic nature of the brain, there has always been a motivation to move beyond 'static' functional connectivity, which characterizes functional interactions over an extended period of time. Progress in data acquisition and advances in analytical neuroimaging methods now allow us to assess the whole brain's dynamic functional connectivity (dFC) and its network-based analog, dynamic functional network connectivity at the macroscale (mm) using fMRI. This has resulted in the rapid growth of analytical approaches, some of which are very complex, requiring technical expertise that could daunt researchers and neuroscientists. Meanwhile, making real progress toward understanding the association between brain dynamism and brain disorders can only be achieved through research conducted by domain experts, such as neuroscientists and psychiatrists. This article aims to provide a gentle introduction to the application of dFC. We first explain what dFC is and the circumstances under which it can be used. Next, we review two major categories of analytical approaches to capture dFC. We discuss caveats and considerations in dFC analysis. Finally, we walk readers through an openly accessible toolbox to capture dFC properties and briefly review some of the dynamic metrics calculated using this toolbox.