Time-varying functional network information extracted from brief instances of spontaneous brain activity

Altered dynamic large-scale brain networks and combined machine learning in primary angle-closure glaucoma

Primary angle-closure glaucoma (PACG) is a severe and irreversible blinding eye disease characterized by progressive retinal ganglion cell death. However, prior research has predominantly focused on static brain activity changes, neglecting the exploration of how PACG impacts the dynamic characteristics of functional brain networks. This study enrolled forty-four patients diagnosed with PACG and forty-four age, gender, and education level-matched healthy controls (HCs). The study employed Independent Component Analysis (ICA) techniques to extract resting-state networks (RSNs) from resting-state functional magnetic resonance imaging (rs-fMRI) data. Subsequently, the RSNs was utilized as the basis for examining and comparing the functional connectivity variations within and between the two groups of resting-state networks. To further explore, a combination of sliding time window and k-means cluster analyses identified seven stable and repetitive dynamic functional network connectivity (dFNC) states. This approach facilitated the comparison of dynamic functional network connectivity and temporal metrics between PACG patients and HCs for each state. Subsequently, a support vector machine (SVM) model leveraging functional connectivity (FC) and FNC was applied to differentiate PACG patients from HCs. Our study underscores the presence of modified functional connectivity within large-scale brain networks and abnormalities in dynamic temporal metrics among PACG patients. By elucidating the impact of changes in large-scale brain networks on disease evolution, researchers may enhance the development of targeted therapies and interventions to preserve vision and cognitive function in PACG.

Real-time fMRI neurofeedback modulates induced hallucinations and underlying brain mechanisms

Hallucinations can occur in the healthy population, are clinically relevant and frequent symptoms in many neuropsychiatric conditions, and have been shown to mark disease progression in patients with neurodegenerative disorders where antipsychotic treatment remains challenging. Here, we combine MR-robotics capable of inducing a clinically-relevant hallucination, with real-time fMRI neurofeedback (fMRI-NF) to train healthy individuals to up-regulate a fronto-parietal brain network associated with the robotically-induced hallucination. Over three days, participants learned to modulate occurrences of and transition probabilities to this network, leading to heightened sensitivity to induced hallucinations after training. Moreover, participants who became sensitive and succeeded in fMRI-NF training, showed sustained and specific neural changes after training, characterized by increased hallucination network occurrences during induction and decreased hallucination network occurrences during a matched control condition. These data demonstrate that fMRI-NF modulates specific hallucination network dynamics and highlights the potential of fMRI-NF as a novel antipsychotic treatment in neurodegenerative disorders and schizophrenia.

Brain network dynamics in patients with single- and multiple-domain amnestic mild cognitive impairment

INTRODUCTION

Brain network dynamics have been extensively explored in patients with amnestic mild cognitive impairment (aMCI); however, differences in single- and multiple-domain aMCI (SD-aMCI and MD-aMCI) remain unclear.

METHODS

Using multicenter datasets, coactivation patterns (CAPs) were constructed and compared among normal control (NC), SD-aMCI, MD-aMCI, and Alzheimer's disease (AD) patients based on individual high-order cognitive network (HOCN) and primary sensory network (PSN) parcellations. Correlations between spatiotemporal characteristics and neuropsychological scores were analyzed.

RESULTS

Compared to NC, SD-aMCI showed temporal alterations in HOCN-dominant CAPs, while MD-aMCI showed alterations in PSN-dominant CAPs. In addition, transitions from SD-aMCI to AD may involve PSN, while MD-aMCI to AD involves both PSN and HOCN. Results were generally consistent across datasets from Chinese and White populations.

DISCUSSION

The HOCN and PSN are distinctively involved in aMCI subtypes and in the transformation between aMCI subtypes and AD, highlighting the necessity of aMCI subtype classification in AD studies.

HIGHLIGHTS

Individual functional network parcellations and coactivation pattern (CAP) analysis were performed to characterize spatiotemporal differences between single- and multiple-domain amnestic mild cognitive impairment (SD-aMCI and MD-aMCI), and between distinct aMCI subtypes and Alzheimer's disease (AD). The analysis of multicenter datasets converged on four pairs of recurrent CAPs, including primary sensory networks (PSN)-dominant CAPs, high-order cognitive networks (HOCN)-dominant CAPs, and PSN-HOCN-interacting CAPs. The HOCN and PSN are distinctively involved in aMCI subtypes and in the transformation between distinct aMCI subtypes and AD.

Alterations of resting-state network dynamics in Alzheimer's disease based on leading eigenvector dynamics analysis

Alzheimer's disease (AD) is a neurodegenerative disease, and mild cognitive impairment (MCI) is considered a transitional stage between healthy aging and dementia. Early detection of MCI can help slow down the progression of AD. At present, there are few studies exploring the characteristics of abnormal dynamic brain activity in AD. This article uses a method called leading eigenvector dynamics analysis (LEiDA) to study resting-state functional magnetic resonance imaging (rs-fMRI) data of AD, MCI, and cognitively normal (CN) participants. By identifying repetitive states of phase coherence, intergroup differences in brain dynamic activity indicators are examined, and the neurobehavioral scales were used to assess the relationship between abnormal dynamic activities and cognitive function. The results showed that in the indicators of occurrence probability and lifetime, the globally synchronized state of the patient group decreased. The activity state of the limbic regions significantly detected the difference between AD and the other two groups. Compared to CN, AD and MCI have varying degrees of increase in default and visual region activity states. In addition, in the analysis related to the cognitive scales, it was found that individuals with poorer cognitive abilities were less active in the globally synchronized state and more active in limbic region activity state and visual region activity state. Taken together, these findings reveal abnormal dynamic activity of resting-state networks in patients with AD and MCI, provide new insights into the dynamic analysis of brain networks, and contribute to a deeper understanding of abnormal spatial dynamic patterns in AD patients.NEW & NOTEWORTHY Alzheimer's disease (AD) is a neurodegenerative disease, but few studies have explored the characteristics of abnormal dynamic brain activity in AD patients. Here, our report reveals the abnormal dynamic activity of the patients' resting-state network, providing new insights into the dynamic analysis of brain networks and helping to gain a deeper understanding of the abnormal spatial dynamic patterns in AD patients.

Distinct Time-Resolved Brain-Wide Coactivations in Oxygenated and Deoxygenated Hemoglobin

OBJECTIVE

Human resting-state networks (RSNs) estimated from oxygenated (HbO) and deoxygenated hemoglobin (HbR) data exhibit strong similarities, while task-based studies show different dynamics in HbR and HbO responses. Such a discrepancy might be explained due to time-averaged estimations of RSNs. Our study investigated differences between HbO and HbR on time-resolved brain-wide coactivation patterns (CAPs).

METHODS

Diffuse optical tomography was reconstructed from resting-state whole-head functional near-infrared spectroscopy data of HbR and HbO in individual healthy participants. Time-averaged RSNs were obtained using the group-level independent component analysis. Time-resolved CAPs were estimated using a clustering approach on the time courses of all obtained RSNs. Characteristics of the RSNs and CAPs from HbR and HbO were compared.

RESULTS

Spatial patterns of HbR and HbO RSNs exhibited significant similarities. Meanwhile, HbR CAPs revealed much more organized spatial and dynamic characteristics than HbO CAPs. The entire set of HbR CAPs suggests a superstructure resulted from brain-wide neuronal dynamics, which is less evident in the set of HbO CAPs. These differences between HbO and HbR CAPs were consistently replicated in individual session data.

CONCLUSION

Our results suggest that human resting brain-wide neuronal activations are preserved better in time-resolved brain-wide patterns, i.e., CAPs, from HbR than those from HbO, while such a difference is lost between time-averaged HbR and HbO RSNs.

SIGNIFICANCE

Our results reveal, for the first time, HbR concentration fluctuations are more directly coupled with resting dynamics of brain-wide neuronal activations in human brains.

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.

A shifting role of thalamocortical connectivity in the emergence of cortical functional organization

The cortical patterning principle has been a long-standing question in neuroscience, yet how this translates to macroscale functional specialization in the human brain remains largely unknown. Here we examine age-dependent differences in resting-state thalamocortical connectivity to investigate its role in the emergence of large-scale functional networks during early life, using a primarily cross-sectional but also longitudinal approach. We show that thalamocortical connectivity during infancy reflects an early differentiation of sensorimotor networks and genetically influenced axonal projection. This pattern changes in childhood, when connectivity is established with the salience network, while decoupling externally and internally oriented functional systems. A developmental simulation using generative network models corroborated these findings, demonstrating that thalamic connectivity contributes to developing key features of the mature brain, such as functional segregation and the sensory-association axis, especially across 12-18 years of age. Our study suggests that the thalamus plays an important role in functional specialization during development, with potential implications for studying conditions with compromised internal and external processing.

Alterations in neural circuit dynamics between the limbic network and prefrontal/default mode network in patients with generalized anxiety disorder

BACKGROUND

Widespread functional alterations have been implicated in patients with generalized anxiety disorder (GAD). However, most studies have primarily focused on static brain network features in patients with GAD. The current research focused on exploring the dynamics within functional brain networks among individuals diagnosed with GAD.

METHODS

Seventy-five participants were divided into patients with GAD and healthy controls (HCs), and resting-state functional magnetic resonance imaging data were collected. The severity of symptoms was measured using the Hamilton Anxiety Scale and the Patient Health Questionnaire. Co-activation pattern (CAP) analysis, centered on the bed nucleus of the stria terminalis, was applied to explore network dynamics. The capability of these dynamic characteristics to distinguish between patients with GAD and HCs was evaluated using a support vector machine.

RESULTS

Patients with GAD exhibited disruptions in the limbic-prefrontal and limbic-default-mode network circuits. Particularly noteworthy was the marked reduction in dynamic indicators such as occurrence, EntriesFromBaseline, ExitsToBaseline, in-degree, out-degree, and resilience. Moreover, these decreased dynamic features effectively distinguished the GAD group from the HC in this study.

CONCLUSIONS

The current findings revealed the underlying brain networks associated with compromised emotion regulation in individuals with GAD. The dynamic reduction in connectivity between the limbic-default mode network and limbic-prefrontal networks could potentially act as a biomarker and therapeutic target for GAD in the future.

Estimation of static and dynamic functional connectivity in resting-state fMRI using zero-frequency resonator

Resting-state functional magnetic resonance imaging (rs-fMRI) is increasingly being used to infer the functional organization of the brain. Blood oxygen level-dependent (BOLD) features related to spontaneous neuronal activity, are yet to be clearly understood. Prior studies have hypothesized that rs-fMRI is spontaneous event-related and these events convey crucial information about the neuronal activity in estimating resting state functional connectivity (FC). Attempts have been made to extract these temporal events using a predetermined threshold. However, the thresholding methods in addition to being very sensitive to noise, may consider redundant events or exclude the low-valued inflection points. Here, we extract the event-related temporal onsets from the rs-fMRI time courses using a zero-frequency resonator (ZFR). The ZFR reflects the transient behavior of the BOLD events at its output. The conditional rate (CR) of the BOLD events occurring in a time course with respect to a seed time course is used to derive static FC. The temporal activity around the estimated events called high signal-to-noise ratio (SNR) segments are also obtained in the rs-fMRI time course and are then used to compute static and dynamic FCs during rest. Coactivation pattern (CAP) is the dynamic FC obtained using the high SNR segments driven by the ZFR. The static FC demonstrates that the ZFR-based CR distinguishes the coactivation and non-coactivation scores well in the distribution. CAP analysis demonstrated the stable and longer dwell time dominant resting state functional networks with high SNR segments driven by the ZFR. Static and dynamic FC analysis underpins that the ZFR-driven temporal onsets of BOLD events derive reliable and consistent FCs in the resting brain using a subset of the time points.

Human brain state dynamics reflect individual neuro-phenotypes

Neural activity and behavior vary within an individual (states) and between individuals (traits). However, the mapping of state-trait neural variation to behavior is not well understood. To address this gap, we quantify moment-to-moment changes in brain-wide co-activation patterns derived from resting-state functional magnetic resonance imaging. In healthy young adults, we identify reproducible spatio-temporal features of co-activation patterns at the single subject level. We demonstrate that a joint analysis of state-trait neural variations and feature reduction reveal general motifs of individual differences, encompassing state-specific and general neural features that exhibit day-to-day variability. The principal neural variations co-vary with the principal variations of behavioral phenotypes, highlighting cognitive function, emotion regulation, alcohol and substance use. Person-specific probability of occupying a particular co-activation pattern is reproducible and associated with neural and behavioral features. This combined analysis of state-trait variations holds promise for developing reproducible neuroimaging markers of individual life functional outcome.

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.

An investigation into the abnormal dynamic connection mechanism of generalized anxiety disorders based on non-homogeneous Markov models

BACKGROUND

The dynamic and hierarchical nature of the functional brain network. The neural dynamical systems tend to converge to multiple attractors (stable fixed points or dynamical states) in long run. Little is known about how the changes in this brain dynamic "long-term" behavior of the connectivity flow of brain network in generalized anxiety disorder (GAD).

METHODS

This study recruited 92 patients with GAD and 77 healthy controls (HC). We applied a reachable probability approach combining a Non-homogeneous Markov model with transition probability to quantify all possible connectivity flows and the hierarchical structure of brain functional systems at the dynamic level and the stationary probability vector (10-step transition probabilities) to describe the steady state of the system in the long run. A random forest algorithm was conducted to predict the severity of anxiety.

RESULTS

The dynamic functional patterns in distributed brain networks had larger possibility to converge in bilateral thalamus, posterior cingulate cortex (PCC), right superior occipital gyrus (SOG) and smaller possibility to converge in bilateral superior temporal gyrus (STG) and right parahippocampal gyrus (PHG) in patients with GAD compared to HC. The abnormal transition probability pattern could predict anxiety severity in patients with GAD.

LIMITATIONS

Small samples and subjects taking medications may have influenced our results. Future studies are expected to rule out the potential confounding effects.

CONCLUSION

Our results have revealed abnormal dynamic neural communication and integration in emotion regulation in patients with GAD, which give new insights to understand the dynamics of brain function of patients with GAD.

Modulation of habenular and nucleus accumbens functional connectivity by ketamine in major depression

INTRODUCTION

Major depressive disorder (MDD) is associated with dysfunctional reward processing, which involves functional circuitry of the habenula (Hb) and nucleus accumbens (NAc). Since ketamine elicits rapid antidepressant and antianhedonic effects in MDD, this study sought to investigate how serial ketamine infusion (SKI) treatment modulates static and dynamic functional connectivity (FC) in Hb and NAc functional networks.

METHODS

MDD participants (n = 58, mean age = 40.7 years, female = 28) received four ketamine infusions (0.5 mg/kg) 2-3 times weekly. Resting-state functional magnetic resonance imaging (fMRI) scans and clinical assessments were collected at baseline and 24 h post-SKI. Static FC (sFC) and dynamic FC variability (dFCv) were calculated from left and right Hb and NAc seeds to all other brain regions. Changes in FC pre-to-post SKI, and correlations with changes with mood and anhedonia were examined. Comparisons of FC between patients and healthy controls (HC) at baseline (n = 55, mean age = 32.6, female = 31), and between HC assessed twice (n = 16) were conducted as follow-up analyses.

RESULTS

Following SKI, significant increases in left Hb-bilateral visual cortex FC, decreases in left Hb-left inferior parietal cortex FC, and decreases in left NAc-right cerebellum FC occurred. Decreased dFCv between left Hb and right precuneus and visual cortex, and decreased dFCv between right NAc and right visual cortex both significantly correlated with improvements in mood ratings. Decreased FC between left Hb and bilateral visual/parietal cortices as well as increased FC between left NAc and right visual/parietal cortices both significantly correlated with improvements in anhedonia. No differences were observed between HC at baseline or over time.

CONCLUSION

Subanesthetic ketamine modulates functional pathways linking the Hb and NAc with visual, parietal, and cerebellar regions in MDD. Overlapping effects between Hb and NAc functional systems were associated with ketamine's therapeutic response.

Coactivation pattern analysis reveals altered whole-brain functional transient dynamics in autism spectrum disorder

Recent studies on autism spectrum disorder (ASD) have identified recurring states dominated by similar coactivation pattern (CAP) and revealed associations between dysfunction in seed-based large-scale brain networks and clinical symptoms. However, the presence of abnormalities in moment-to-moment whole-brain dynamics in ASD remains uncertain. In this study, we employed seed-free CAP analysis to identify transient brain activity configurations and investigate dynamic abnormalities in ASD. We utilized a substantial multisite resting-state fMRI dataset consisting of 354 individuals with ASD and 446 healthy controls (HCs, from HC groups and 2). CAP were generated from a subgroup of all HC subjects (HC group 1) through temporal K-means clustering, identifying four CAPs. These four CAPs exhibited either the activation or inhibition of the default mode network (DMN) and were grouped into two pairs with opposing spatial CAPs. CAPs for HC group 2 and ASD were identified by their spatial similarity to those for HC group 1. Compared with individuals in HC group 2, those with ASD spent more time in CAPs involving the ventral attention network but less time in CAPs related to executive control and the dorsal attention network. Support vector machine analysis demonstrated that the aberrant dynamic characteristics of CAPs achieved an accuracy of 74.87% in multisite classification. In addition, we used whole-brain dynamics to predict symptom severity in ASD. Our findings revealed whole-brain dynamic functional abnormalities in ASD from a single transient perspective, emphasizing the importance of the DMN in abnormal dynamic functional activity in ASD and suggesting that temporally dynamic techniques offer novel insights into time-varying neural processes.

Blink-related arousal network surges are shaped by cortical vigilance states

The vigilance state and the excitability of cortical networks impose wide-range effects on brain dynamics that arousal surges could promptly modify. We previously reported an association between spontaneous eye-blinks and BOLD activation in the brain arousal ascending network (AAN) and in thalamic nuclei based on 3T MR resting state brain images. Here we aimed to replicate our analyses using 7T MR images in a larger cohort of participants collected from the Human Connectome Project (HCP), which also contained simultaneous eye-tracking recordings, and to assess the interaction between the blink-associated arousal surges and the vigilance states. For this purpose, we compared blink associated BOLD activity under a vigilant versus a drowsy state, a classification made based on the pupillary data obtained during the fMRI scans. We conducted two main analyses: i) Cross-correlation analysis between the BOLD signal and blink events (eye blink time-series were convolved with the canonical and also with the temporal derivative of the Hemodynamic Response Function, HRF) within preselected regions of interests (ROIs) (i.e., brainstem AAN, thalamic and cerebellar nuclei) together with an exploratory voxel-wise analyses to assess the whole-brain, and ii) blink-event analysis of the BOLD signals to reveal the signal changes onset to the blinks in the preselected ROIs. Consistent with our prior findings on 3T MRI, we showed significant positive cross correlations between BOLD peaks in brainstem and thalamic nuclei that preceded or were overlapping with blink moments and that sharply decreased post-blink. Whole brain analysis revealed blink-related activation that was strongest in cerebellum, insula, lateral geniculate nucleus (LGN) and visual cortex. Drowsiness impacted HRF BOLD (enhancing it), time-to-peak (delaying it) and post-blink BOLD activity (accentuating decreases). Responses in the drowsy state could be related to the differences in the excitability of cortical, subcortical and cerebellar tissue, such that cerebellar and thalamic regions involved in visual attention processing were more responsive for the vigilant state, but AAN ROIs, as well as cerebellar and thalamic ROIs connected to pre-motor, frontal, temporal and DMN regions were less responsive. Such qualitative and quantitative differences in the blink related BOLD signal changes could reflect delayed cortical processing and the ineffectiveness of arousal surges during states of drowsiness. Future studies that manipulate arousal are needed to corroborate a mechanistic interaction of arousal surges with vigilance states and cortical excitability.

Spontaneous persistent activity and inactivity in vivo reveals differential cortico-entorhinal functional connectivity

Understanding the functional connectivity between brain regions and its emergent dynamics is a central challenge. Here we present a theory-experiment hybrid approach involving iteration between a minimal computational model and in vivo electrophysiological measurements. Our model not only predicted spontaneous persistent activity (SPA) during Up-Down-State oscillations, but also inactivity (SPI), which has never been reported. These were confirmed in vivo in the membrane potential of neurons, especially from layer 3 of the medial and lateral entorhinal cortices. The data was then used to constrain two free parameters, yielding a unique, experimentally determined model for each neuron. Analytic and computational analysis of the model generated a dozen quantitative predictions about network dynamics, which were all confirmed in vivo to high accuracy. Our technique predicted functional connectivity; e. g. the recurrent excitation is stronger in the medial than lateral entorhinal cortex. This too was confirmed with connectomics data. This technique uncovers how differential cortico-entorhinal dialogue generates SPA and SPI, which could form an energetically efficient working-memory substrate and influence the consolidation of memories during sleep. More broadly, our procedure can reveal the functional connectivity of large networks and a theory of their emergent dynamics.

The influence of accelerated brain aging on coactivation pattern dynamics in Parkinson's disease

Aging is widely acknowledged as the primary risk factor for brain degeneration, with Parkinson's disease (PD) tending to follow accelerated aging trajectories. We aim to investigate the impact of structural brain aging on the temporal dynamics of a large-scale functional network in PD. We enrolled 62 PD patients and 32 healthy controls (HCs). The level of brain aging was determined by calculating global and local brain age gap estimates (G-brainAGE and L-brainAGE) from structural images. The neural network activity of the whole brain was captured by identifying coactivation patterns (CAPs) from resting-state functional images. Intergroup differences were assessed using the general linear model. Subsequently, a spatial correlation analysis between the L-brainAGE difference map and CAPs was conducted to uncover the anatomical underpinnings of functional alterations. Compared to HCs (-3.73 years), G-brainAGE was significantly higher in PD patients (+1.93 years), who also exhibited widespread elevation in L-brainAGE. G-brainAGE was correlated with disease severity and duration. PD patients spent less time in CAPs involving activated default mode and the fronto-parietal network (DMN-FPN), as well as the sensorimotor and salience network (SMN-SN), and had a reduced transition frequency from other CAPs to the DMN-FPN and SMN-SN CAPs. Furthermore, the pattern of localized brain age acceleration showed spatial similarities with the SMN-SN CAP. Accelerated structural brain aging in PD adversely affects brain function, manifesting as dysregulated brain network dynamics. These findings provide insights into the neuropathological mechanisms underlying neurodegenerative diseases and imply the possibility of interventions for modifying PD progression by slowing the brain aging process.

Early altered directionality of resting brain network state transitions in the TgF344-AD rat model of Alzheimer's disease

Introduction

Alzheimer's disease (AD) is a progressive neurodegenerative disease resulting in memory loss and cognitive decline. Synaptic dysfunction is an early hallmark of the disease whose effects on whole-brain functional architecture can be identified using resting-state functional MRI (rsfMRI). Insights into mechanisms of early, whole-brain network alterations can help our understanding of the functional impact of AD's pathophysiology.

Methods

Here, we obtained rsfMRI data in the TgF344-AD rat model at the pre- and early-plaque stages. This model recapitulates the major pathological and behavioral hallmarks of AD. We used co-activation pattern (CAP) analysis to investigate if and how the dynamic organization of intrinsic brain functional networks states, undetectable by earlier methods, is altered at these early stages.

Results

We identified and characterized six intrinsic brain states as CAPs, their spatial and temporal features, and the transitions between the different states. At the pre-plaque stage, the TgF344-AD rats showed reduced co-activation of hub regions in the CAPs corresponding to the default mode-like and lateral cortical network. Default mode-like network activity segregated into two distinct brain states, with one state characterized by high co-activation of the basal forebrain. This basal forebrain co-activation was reduced in TgF344-AD animals mainly at the pre-plaque stage. Brain state transition probabilities were altered at the pre-plaque stage between states involving the default mode-like network, lateral cortical network, and basal forebrain regions. Additionally, while the directionality preference in the network-state transitions observed in the wild-type animals at the pre-plaque stage had diminished at the early-plaque stage, TgF344-AD animals continued to show directionality preference at both stages.

Discussion

Our study enhances the understanding of intrinsic brain state dynamics and how they are impacted at the early stages of AD, providing a nuanced characterization of the early, functional impact of the disease's neurodegenerative process.

Sex-specific resting state brain network dynamics in patients with major depressive disorder

Sex-specific neurobiological changes have been implicated in Major Depressive Disorder (MDD). Dysfunctions of the default mode network (DMN), salience network (SN) and frontoparietal network (FPN) are critical neural characteristics of MDD, however, the potential moderating role of sex on resting-state network dynamics in MDD has not been sufficiently evaluated. Thus, resting-state functional magnetic resonance imaging (fMRI) data were collected from 138 unmedicated patients with first-episode MDD (55 males) and 243 healthy controls (HCs; 106 males). Recurring functional network co-activation patterns (CAPs) were extracted, and time spent in each CAP (the total amount of volumes associated to a CAP), persistence (the average number of consecutive volumes linked to a CAP), and transitions across CAPs involving the SN, DMN and FPN were quantified. Relative to HCs, MDD patients exhibited greater persistence in a CAP involving activation of the DMN and deactivation of the FPN (DMN + FPN-). In addition, relative to the sex-matched HCs, the male MDD group spent more time in two CAPs involving the SN and DMN (i.e., DMN + SN- and DMN-SN + ) and transitioned more frequently from the DMN + FPN- CAP to the DMN + SN- CAP relative to the male HC group. Conversely, the female MDD group showed less persistence in the DMN + SN- CAP relative to the female HC group. Our findings highlight that the imbalance between SN and DMN could be a neurobiological marker supporting sex differences in MDD. Moreover, the dominance of the DMN accompanied by the deactivation of the FPN could be a sex-independent neurobiological correlate related to depression.

The co-activation patterns of multiple brain regions in Juvenile Myoclonic Epilepsy

Juvenile myoclonic epilepsy (JME) as an idiopathic generalized epilepsy has been studied by many advanced neuroimaging techniques to elucidate its neuroanatomical basis and pathophysiological mechanisms. In this paper, we used co-activation patterns (CAPs) to explore the differences of dynamic brain activity changes in resting state between JME patients and healthy controls. 27 cases JME patients and 27 cases healthy of fMRI data were collected. The structural image data of the subjects were analyzed by voxel-based morphological analysis, and the regions with gray matter volume atrophy and high voxel were selected as the regions of interest. Further, the mean disease duration was used as boundary to divide the patients' data into the below-average time and the above-average time groups, which were defined as patient disease duration groups. And these data were used to construct CAPs and to compare changes in brain dynamics. It was found that the number of patterns occurrences and the possibility of switching between patterns were smaller than those in the healthy control, which indicated patients with damage to brain regions. For the patient time control group, the number of patterns occurrences and the possibility of switching between patterns were similar, while there was linear regression between the three values and disease duration. Collectively, this study provides important evidence for revealing the key brain regions of JME by studying the transformation between CAPs. Future studies could investigate the effects of receiving treatment on patient dynamic brain activity.

Thalamic contributions to psychosis susceptibility: Evidence from co-activation patterns accounting for intra-seed spatial variability (μCAPs)

The temporal variability of the thalamus in functional networks may provide valuable insights into the pathophysiology of schizophrenia. To address the complexity of the role of the thalamic nuclei in psychosis, we introduced micro-co-activation patterns (μCAPs) and employed this method on the human genetic model of schizophrenia 22q11.2 deletion syndrome (22q11.2DS). Participants underwent resting-state functional MRI and a data-driven iterative process resulting in the identification of six whole-brain μCAPs with specific activity patterns within the thalamus. Unlike conventional methods, μCAPs extract dynamic spatial patterns that reveal partially overlapping and non-mutually exclusive functional subparts. Thus, the μCAPs method detects finer foci of activity within the initial seed region, retaining valuable and clinically relevant temporal and spatial information. We found that a μCAP showing co-activation of the mediodorsal thalamus with brain-wide cortical regions was expressed significantly less frequently in patients with 22q11.2DS, and its occurrence negatively correlated with the severity of positive psychotic symptoms. Additionally, activity within the auditory-visual cortex and their respective geniculate nuclei was expressed in two different μCAPs. One of these auditory-visual μCAPs co-activated with salience areas, while the other co-activated with the default mode network (DMN). A significant shift of occurrence from the salience+visuo-auditory-thalamus to the DMN + visuo-auditory-thalamus μCAP was observed in patients with 22q11.2DS. Thus, our findings support existing research on the gatekeeping role of the thalamus for sensory information in the pathophysiology of psychosis and revisit the evidence of geniculate nuclei hyperconnectivity with the audio-visual cortex in 22q11.2DS in the context of dynamic functional connectivity, seen here as the specific hyper-occurrence of these circuits with the task-negative brain networks.

The brain entropy dynamics in resting state

As a novel measure for irregularity and complexity of the spontaneous fluctuations of brain activities, brain entropy (BEN) has attracted much attention in resting-state functional magnetic resonance imaging (rs-fMRI) studies during the last decade. Previous studies have shown its associations with cognitive and mental functions. While most previous research assumes BEN is approximately stationary during scan sessions, the brain, even at its resting state, is a highly dynamic system. Such dynamics could be characterized by a series of reoccurring whole-brain patterns related to cognitive and mental processes. The present study aims to explore the time-varying feature of BEN and its potential links with general cognitive ability. We adopted a sliding window approach to derive the dynamical brain entropy (dBEN) of the whole-brain functional networks from the HCP (Human Connectome Project) rs-fMRI dataset that includes 812 young healthy adults. The dBEN was further clustered into 4 reoccurring BEN states by the k-means clustering method. The fraction window (FW) and mean dwell time (MDT) of one BEN state, characterized by the extremely low overall BEN, were found to be negatively correlated with general cognitive abilities (i.e., cognitive flexibility, inhibitory control, and processing speed). Another BEN state, characterized by intermediate overall BEN and low within-state BEN located in DMN, ECN, and part of SAN, its FW, and MDT were positively correlated with the above cognitive abilities. The results of our study advance our understanding of the underlying mechanism of BEN dynamics and provide a potential framework for future investigations in clinical populations.

Dynamic aberrances of substantia nigra-relevant coactivation patterns in first-episode treatment-naïve patients with schizophrenia

BACKGROUND

Although dopaminergic disturbances are well-known in schizophrenia, the understanding of dopamine-related brain dynamics remains limited. This study investigates the dynamic coactivation patterns (CAPs) associated with the substantia nigra (SN), a key dopaminergic nucleus, in first-episode treatment-naïve patients with schizophrenia (FES).

METHODS

Resting-state fMRI data were collected from 84 FES and 94 healthy controls (HCs). Frame-wise clustering was implemented to generate CAPs related to SN activation or deactivation. Connectome features of each CAP were derived using an edge-centric method. The occurrence for each CAP and the balance ratio for antagonistic CAPs were calculated and compared between two groups, and correlations between temporal dynamic metrics and symptom burdens were explored.

RESULTS

Functional reconfigurations in CAPs exhibited significant differences between the activation and deactivation states of SN. During SN activation, FES more frequently recruited a CAP characterized by activated default network, language network, control network, and the caudate, compared to HCs (F = 8.54, FDR-p = 0.030). Moreover, FES displayed a tilted balance towards a CAP featuring SN-coactivation with the control network, caudate, and thalamus, as opposed to its antagonistic CAP (F = 7.48, FDR-p = 0.030). During SN deactivation, FES exhibited increased recruitment of a CAP with activated visual and dorsal attention networks but decreased recruitment of its opposing CAP (F = 6.58, FDR-p = 0.034).

CONCLUSION

Our results suggest that neuroregulatory dysfunction in dopaminergic pathways involving SN potentially mediates aberrant time-varying functional reorganizations in schizophrenia. This finding enriches the dopamine hypothesis of schizophrenia from the perspective of brain dynamics.

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.

Mapping dynamic spatial patterns of brain function with spatial-wise attention

Objective: Using functional magnetic resonance imaging (fMRI) and deep learning to discover the spatial pattern of brain function, or functional brain networks (FBNs) has been attracted many reseachers. Most existing works focus on static FBNs or dynamic functional connectivity among fixed spatial network nodes, but ignore the potential dynamic/time-varying characteristics of the spatial networks themselves. And most of works based on the assumption of linearity and independence, that oversimplify the relationship between blood-oxygen level dependence signal changes and the heterogeneity of neuronal activity within voxels.Approach: To overcome these problems, we proposed a novel spatial-wise attention (SA) based method called Spatial and Channel-wise Attention Autoencoder (SCAAE) to discover the dynamic FBNs without the assumptions of linearity or independence. The core idea of SCAAE is to apply the SA to generate FBNs directly, relying solely on the spatial information present in fMRI volumes. Specifically, we trained the SCAAE in a self-supervised manner, using the autoencoder to guide the SA to focus on the activation regions. Experimental results show that the SA can generate multiple meaningful FBNs at each fMRI time point, which spatial similarity are close to the FBNs derived by known classical methods, such as independent component analysis.Main results: To validate the generalization of the method, we evaluate the approach on HCP-rest, HCP-task and ADHD-200 dataset. The results demonstrate that SA mechanism can be used to discover time-varying FBNs, and the identified dynamic FBNs over time clearly show the process of time-varying spatial patterns fading in and out.Significance: Thus we provide a novel method to understand human brain better. Code is available athttps://github.com/WhatAboutMyStar/SCAAE.

Spatiotemporal dynamics of hippocampal-cortical networks underlying the unique phenomenological properties of trauma-related intrusive memories

Trauma-related intrusive memories (TR-IMs) possess unique phenomenological properties that contribute to adverse post-traumatic outcomes, positioning them as critical intervention targets. However, transdiagnostic treatments for TR-IMs are scarce, as their underlying mechanisms have been investigated separate from their unique phenomenological properties. Extant models of more general episodic memory highlight dynamic hippocampal-cortical interactions that vary along the anterior-posterior axis of the hippocampus (HPC) to support different cognitive-affective and sensory-perceptual features of memory. Extending this work into the unique properties of TR-IMs, we conducted a study of eighty-four trauma-exposed adults who completed daily ecological momentary assessments of TR-IM properties followed by resting-state functional magnetic resonance imaging (rs-fMRI). Spatiotemporal dynamics of anterior and posterior hippocampal (a/pHPC)-cortical networks were assessed using co-activation pattern analysis to investigate their associations with different properties of TR-IMs. Emotional intensity of TR-IMs was inversely associated with the frequency and persistence of an aHPC-default mode network co-activation pattern. Conversely, sensory features of TR-IMs were associated with more frequent co-activation of the HPC with sensory cortices and the ventral attention network, and the reliving of TR-IMs in the "here-and-now" was associated with more persistent co-activation of the pHPC and the visual cortex. Notably, no associations were found between HPC-cortical network dynamics and conventional symptom measures, including TR-IM frequency or retrospective recall, underscoring the utility of ecological assessments of memory properties in identifying their neural substrates. These findings provide novel insights into the neural correlates of the unique features of TR-IMs that are critical for the development of individualized, transdiagnostic treatments for this pervasive, difficult-to-treat symptom.

State-transition dynamics of resting-state functional magnetic resonance imaging data: model comparison and test-to-retest analysis

Electroencephalogram (EEG) microstate analysis entails finding dynamics of quasi-stable and generally recurrent discrete states in multichannel EEG time series data and relating properties of the estimated state-transition dynamics to observables such as cognition and behavior. While microstate analysis has been widely employed to analyze EEG data, its use remains less prevalent in functional magnetic resonance imaging (fMRI) data, largely due to the slower timescale of such data. In the present study, we extend various data clustering methods used in EEG microstate analysis to resting-state fMRI data from healthy humans to extract their state-transition dynamics. We show that the quality of clustering is on par with that for various microstate analyses of EEG data. We then develop a method for examining test-retest reliability of the discrete-state transition dynamics between fMRI sessions and show that the within-participant test-retest reliability is higher than between-participant test-retest reliability for different indices of state-transition dynamics, different networks, and different data sets. This result suggests that state-transition dynamics analysis of fMRI data could discriminate between different individuals and is a promising tool for performing fingerprinting analysis of individuals.

Transient resting-state salience-limbic co-activation patterns in functional neurological disorders

BACKGROUND

Functional neurological disorders were historically regarded as the manifestation of a dynamic brain lesion which might be linked to trauma or stress, although this association has not yet been directly tested yet. Analysing large-scale brain network dynamics at rest in relation to stress biomarkers assessed by salivary cortisol and amylase could provide new insights into the pathophysiology of functional neurological symptoms.

METHODS

Case-control resting-state functional magnetic resonance imaging study of 79 patients with mixed functional neurological disorders (i.e., functional movement disorders, functional seizures, persistent perceptual-postural dizziness) and 74 age- and sex-matched healthy controls. Using a two-step hierarchical data-driven neuroimaging approach, static functional connectivity was first computed between 17 resting-state networks. Second, dynamic alterations in these networks were examined using co-activation pattern analysis. Using a partial least squares correlation analysis, the multivariate pattern of correlation between altered temporal characteristics and stress biomarkers as well as clinical scores were evaluated.

RESULTS

Compared to healthy controls, patients presented with functional aberrancies of the salience-limbic network connectivity. Thus, the insula and amygdala were selected as seed-regions for the subsequent analyses. Insular co-(de)activation patterns related to the salience network, the somatomotor network and the default mode network were detected, which patients entered more frequently than controls. Moreover, an insular co-(de)activation pattern with subcortical regions together with a wide-spread co-(de)activation with diverse cortical networks was detected, which patients entered less frequently than controls. In patients, dynamic alterations conjointly correlated with amylase measures and duration of symptoms.

CONCLUSION

The relationship between alterations in insular co-activation patterns, stress biomarkers and clinical data proposes inter-related mechanisms involved in stress regulation and functional (network) integration. In summary, altered functional brain network dynamics were identified in patients with functional neurological disorder supporting previously raised concepts of impaired attentional and interoceptive processing.

Controlling the human connectome with spatially diffuse input signals

The human brain is never at "rest"; its activity is constantly fluctuating over time, transitioning from one brain state-a whole-brain pattern of activity-to another. Network control theory offers a framework for understanding the effort - energy - associated with these transitions. One branch of control theory that is especially useful in this context is "optimal control", in which input signals are used to selectively drive the brain into a target state. Typically, these inputs are introduced independently to the nodes of the network (each input signal is associated with exactly one node). Though convenient, this input strategy ignores the continuity of cerebral cortex - geometrically, each region is connected to its spatial neighbors, allowing control signals, both exogenous and endogenous, to spread from their foci to nearby regions. Additionally, the spatial specificity of brain stimulation techniques is limited, such that the effects of a perturbation are measurable in tissue surrounding the stimulation site. Here, we adapt the network control model so that input signals have a spatial extent that decays exponentially from the input site. We show that this more realistic strategy takes advantage of spatial dependencies in structural connectivity and activity to reduce the energy (effort) associated with brain state transitions. We further leverage these dependencies to explore near-optimal control strategies such that, on a per-transition basis, the number of input signals required for a given control task is reduced, in some cases by two orders of magnitude. This approximation yields network-wide maps of input site density, which we compare to an existing database of functional, metabolic, genetic, and neurochemical maps, finding a close correspondence. Ultimately, not only do we propose a more efficient framework that is also more adherent to well-established brain organizational principles, but we also posit neurobiologically grounded bases for optimal control.

Spatial-temporal convolutional attention for discovering and characterizing functional brain networks in task fMRI

Functional brain networks (FBNs) are spatial patterns of brain function that play a critical role in understanding human brain function. There are many proposed methods for mapping the spatial patterns of brain function, however they oversimplify the underlying assumptions of brain function and have various limitations such as linearity and independence. Additionally, current methods fail to account for the dynamic nature of FBNs, which limits their effectiveness in accurately characterizing these networks. To address these limitations, we present a novel deep learning and spatial-wise attention based model called Spatial-Temporal Convolutional Attention (STCA) to accurately model dynamic FBNs. Specifically, we train STCA in a self-supervised manner by utilizing a Convolutional Autoencoder to guide the STCA module in assigning higher attention weights to regions of functional activity. To validate the reliability of the results, we evaluate our approach on the HCP-task motor behavior dataset, the experimental results demonstrate that the STCA derived FBNs have higher spatial similarity with the templates and that the spatial similarity between the templates and the FBNs derived by STCA fluctuates with the task design over time, suggesting that STCA can reflect the dynamic changes of brain function, providing a powerful tool to better understand human brain function. Code is available at https://github.com/SNNUBIAI/STCAE.

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.

Modular subgraphs in large-scale connectomes underpin spontaneous co-fluctuation events in mouse and human brains

Previous studies have adopted an edge-centric framework to study fine-scale network dynamics in human fMRI. To date, however, no studies have applied this framework to data collected from model organisms. Here, we analyze structural and functional imaging data from lightly anesthetized mice through an edge-centric lens. We find evidence of "bursty" dynamics and events - brief periods of high-amplitude network connectivity. Further, we show that on a per-frame basis events best explain static FC and can be divided into a series of hierarchically-related clusters. The co-fluctuation patterns associated with each cluster centroid link distinct anatomical areas and largely adhere to the boundaries of algorithmically detected functional brain systems. We then investigate the anatomical connectivity undergirding high-amplitude co-fluctuation patterns. We find that events induce modular bipartitions of the anatomical network of inter-areal axonal projections. Finally, we replicate these same findings in a human imaging dataset. In summary, this report recapitulates in a model organism many of the same phenomena observed in previously edge-centric analyses of human imaging data. However, unlike human subjects, the murine nervous system is amenable to invasive experimental perturbations. Thus, this study sets the stage for future investigation into the causal origins of fine-scale brain dynamics and high-amplitude co-fluctuations. Moreover, the cross-species consistency of the reported findings enhances the likelihood of future translation.

The neural correlates of arousal: Ventral posterolateral nucleus-global transient co-activation

Arousal and awareness are two components of consciousness whose neural mechanisms remain unclear. Spontaneous peaks of global (brain-wide) blood-oxygenation-level-dependent (BOLD) signal have been found to be sensitive to changes in arousal. By contrasting BOLD signals at different arousal levels, we find decreased activation of the ventral posterolateral nucleus (VPL) during transient peaks in the global signal in low arousal and awareness states (non-rapid eye movement sleep and anesthesia) compared to wakefulness and in eyes-closed compared to eyes-open conditions in healthy awake individuals. Intriguingly, VPL-global co-activation remains high in patients with unresponsive wakefulness syndrome (UWS), who exhibit high arousal without awareness, while it reduces in rapid eye movement sleep, a state characterized by low arousal but high awareness. Furthermore, lower co-activation is found in individuals during N3 sleep compared to patients with UWS. These results demonstrate that co-activation of VPL and global activity is critical to arousal but not to awareness.

Individual variability in neural representations of mind-wandering

Mind-wandering is a frequent, daily mental activity, experienced in unique ways in each person. Yet neuroimaging evidence relating mind-wandering to brain activity, for example in the default mode network (DMN), has relied on population-rather than individual-based inferences due to limited within-individual sampling. Here, three densely-sampled individuals each reported hundreds of mind-wandering episodes while undergoing multi-session functional magnetic resonance imaging. We found reliable associations between mind-wandering and DMN activation when estimating brain networks within individuals using precision functional mapping. However, the timing of spontaneous DMN activity relative to subjective reports, and the networks beyond DMN that were activated and deactivated during mind-wandering, were distinct across individuals. Connectome-based predictive modeling further revealed idiosyncratic, whole-brain functional connectivity patterns that consistently predicted mind-wandering within individuals but did not fully generalize across individuals. Predictive models of mind-wandering and attention that were derived from larger-scale neuroimaging datasets largely failed when applied to densely-sampled individuals, further highlighting the need for personalized models. Our work offers novel evidence for both conserved and variable neural representations of self-reported mind-wandering in different individuals. The previously-unrecognized inter-individual variations reported here underscore the broader scientific value and potential clinical utility of idiographic approaches to brain-experience associations.

Dynamical models reveal anatomically reliable attractor landscapes embedded in resting state brain networks

Analyses of functional connectivity (FC) in resting-state brain networks (RSNs) have generated many insights into cognition. However, the mechanistic underpinnings of FC and RSNs are still not well-understood. It remains debated whether resting state activity is best characterized as noise-driven fluctuations around a single stable state, or instead, as a nonlinear dynamical system with nontrivial attractors embedded in the RSNs. Here, we provide evidence for the latter, by constructing whole-brain dynamical systems models from individual resting-state fMRI (rfMRI) recordings, using the Mesoscale Individualized NeuroDynamic (MINDy) platform. The MINDy models consist of hundreds of neural masses representing brain parcels, connected by fully trainable, individualized weights. We found that our models manifested a diverse taxonomy of nontrivial attractor landscapes including multiple equilibria and limit cycles. However, when projected into anatomical space, these attractors mapped onto a limited set of canonical RSNs, including the default mode network (DMN) and frontoparietal control network (FPN), which were reliable at the individual level. Further, by creating convex combinations of models, bifurcations were induced that recapitulated the full spectrum of dynamics found via fitting. These findings suggest that the resting brain traverses a diverse set of dynamics, which generates several distinct but anatomically overlapping attractor landscapes. Treating rfMRI as a unimodal stationary process (i.e., conventional FC) may miss critical attractor properties and structure within the resting brain. Instead, these may be better captured through neural dynamical modeling and analytic approaches. The results provide new insights into the generative mechanisms and intrinsic spatiotemporal organization of brain networks.

Aberrant resting-state co-activation network dynamics in major depressive disorder

Major depressive disorder (MDD) is a globally prevalent and highly disabling disease characterized by dysfunction of large-scale brain networks. Previous studies have found that static functional connectivity is not sufficient to reflect the complicated and time-varying properties of the brain. The underlying dynamic interactions between brain functional networks of MDD remain largely unknown, and it is also unclear whether neuroimaging-based dynamic properties are sufficiently robust to discriminate individuals with MDD from healthy controls since the diagnosis of MDD mainly depends on symptom-based criteria evaluated by clinical observation. Resting-state functional magnetic resonance imaging (fMRI) data of 221 MDD patients and 215 healthy controls were shared by REST-meta-MDD consortium. We investigated the spatial-temporal dynamics of MDD using co-activation pattern analysis and made individual diagnoses using support vector machine (SVM). We found that MDD patients exhibited aberrant dynamic properties (such as dwell time, occurrence rate, transition probability, and entropy of Markov trajectories) in some transient networks including subcortical network (SCN), activated default mode network (DMN), de-activated SCN-cerebellum network, a joint network, activated attention network (ATN), and de-activated DMN-ATN, where some dynamic properties were indicative of depressive symptoms. The trajectories of other networks to deactivated DMN-ATN were more accessible in MDD patients. Subgroup analyses also showed subtle dynamic changes in first-episode drug-naïve (FEDN) MDD patients. Finally, SVM achieved preferable accuracies of 84.69%, 76.77%, and 88.10% in discriminating patients with MDD, FEDN MDD, and recurrent MDD from healthy controls with their dynamic metrics. Our findings reveal that MDD is characterized by aberrant dynamic fluctuations of brain network and the feasibility of discriminating MDD patients using dynamic properties, which provide novel insights into the neural mechanism of MDD.

On the variability of dynamic functional connectivity assessment methods

BACKGROUND

Dynamic functional connectivity (dFC) has become an important measure for understanding brain function and as a potential biomarker. However, various methodologies have been developed for assessing dFC, and it is unclear how the choice of method affects the results. In this work, we aimed to study the results variability of commonly used dFC methods.

METHODS

We implemented 7 dFC assessment methods in Python and used them to analyze the functional magnetic resonance imaging data of 395 subjects from the Human Connectome Project. We measured the similarity of dFC results yielded by different methods using several metrics to quantify overall, temporal, spatial, and intersubject similarity.

RESULTS

Our results showed a range of weak to strong similarity between the results of different methods, indicating considerable overall variability. Somewhat surprisingly, the observed variability in dFC estimates was found to be comparable to the expected functional connectivity variation over time, emphasizing the impact of methodological choices on the final results. Our findings revealed 3 distinct groups of methods with significant intergroup variability, each exhibiting distinct assumptions and advantages.

CONCLUSIONS

Overall, our findings shed light on the impact of dFC assessment analytical flexibility and highlight the need for multianalysis approaches and careful method selection to capture the full range of dFC variation. They also emphasize the importance of distinguishing neural-driven dFC variations from physiological confounds and developing validation frameworks under a known ground truth. To facilitate such investigations, we provide an open-source Python toolbox, PydFC, which facilitates multianalysis dFC assessment, with the goal of enhancing the reliability and interpretability of dFC studies.

Psychedelic concentrations of nitrous oxide reduce functional differentiation in frontoparietal and somatomotor cortical networks

Despite the longstanding use of nitrous oxide and descriptions of its psychological effects more than a century ago, there is a paucity of neurobiological investigation of associated psychedelic experiences. We measure the brain's functional geometry (through analysis of cortical gradients) and temporal dynamics (through analysis of co-activation patterns) using human resting-state functional magnetic resonance imaging data acquired before and during administration of 35% nitrous oxide. Both analyses demonstrate that nitrous oxide reduces functional differentiation in frontoparietal and somatomotor networks. Importantly, the subjective psychedelic experience induced by nitrous oxide is inversely correlated with the degree of functional differentiation. Thus, like classical psychedelics acting on serotonin receptors, nitrous oxide flattens the functional geometry of the cortex and disrupts temporal dynamics in association with psychoactive effects.

Characterizing systemic physiological effects on the blood oxygen level dependent signal of resting-state fMRI in time-frequency space using wavelets

Systemic physiological dynamics, such as heart rate variability (HRV) and respiration volume per time (RVT), are known to account for significant variance in the blood oxygen level dependent (BOLD) signal of resting-state functional magnetic resonance imaging (rsfMRI). However, synchrony between these cardiorespiratory changes and the BOLD signal could be due to neuronal (i.e., autonomic activity inducing changes in heart rate and respiration) or vascular (i.e., cardiorespiratory activity facilitating hemodynamic changes and thus the BOLD signal) effects and the contributions of these effects may differ spatially, temporally, and spectrally. In this study, we characterize these brain-body dynamics using a wavelet analysis in rapidly sampled rsfMRI data with simultaneous pulse oximetry and respiratory monitoring of the Human Connectome Project. Our time-frequency analysis across resting-state networks (RSNs) revealed differences in the coherence of the BOLD signal and heartbeat interval (HBI)/RVT dynamics across frequencies, with unique profiles per network. Somatomotor (SMN), visual (VN), and salience (VAN) networks demonstrated the greatest synchrony with both systemic physiological signals when compared to other networks; however, significant coherence was observed in all RSNs regardless of direct autonomic involvement. Our phase analysis revealed distinct frequency profiles of percentage of time with significant coherence between BOLD and systemic physiological signals for different phase offsets across RSNs, suggesting that the phase offset and temporal order of signals varies by frequency. Lastly, our analysis of temporal variability of coherence provides insight on potential influence of autonomic state on brain-body communication. Overall, the novel wavelet analysis enables an efficient characterization of the dynamic relationship between cardiorespiratory activity and the BOLD signal in spatial, temporal, and spectral dimensions to inform our understanding of autonomic states and improve our interpretation of the BOLD signal.

Ketamine treatment modulates habenular and nucleus accumbens static and dynamic functional connectivity in major depression

Dysfunctional reward processing in major depressive disorder (MDD) involves functional circuitry of the habenula (Hb) and nucleus accumbens (NAc). Ketamine elicits rapid antidepressant and alleviates anhedonia in MDD. To clarify how ketamine perturbs reward circuitry in MDD, we examined how serial ketamine infusions (SKI) modulate static and dynamic functional connectivity (FC) in Hb and NAc networks. MDD participants (n=58, mean age=40.7 years, female=28) received four ketamine infusions (0.5mg/kg) 2-3 times weekly. Resting-state fMRI scans and clinical assessments were collected at baseline and 24 hours post-SKI completion. Static FC (sFC) and dynamic FC variability (dFCv) were calculated from left and right Hb and NAc seeds to all other brain regions. Paired t-tests examined changes in FC pre-to-post SKI, and correlations were used to determine relationships between FC changes with mood and anhedonia. Following SKI, significant increases in left Hb-bilateral visual cortex FC, decreases in left Hb-left inferior parietal cortex FC, and decreases in left NAc-right cerebellum FC occurred. Decreased dFCv between left Hb and right precuneus and visual cortex, and decreased dFCv between right NAc and right visual cortex both significantly correlated with improvements in Hamilton Depression Rating Scale. Decreased FC between left Hb and bilateral visual/parietal cortices as well as increased FC between left NAc and right visual/parietal cortices both significantly correlated with improvements in anhedonia. Subanesthetic ketamine modulates functional pathways linking the Hb and NAc with visual, parietal, and cerebellar regions. Overlapping effects between Hb and NAc functional systems were associated with ketamine's therapeutic response.

Spatial and Spectral Components of the BOLD Global Signal in Rat Resting-State Functional MRI

PURPOSE

In resting-state fMRI (rs-fMRI), the global signal average captures widespread fluctuations related to unwanted sources of variance such as motion and respiration, as well as widespread neural activity; however, relative contributions of neural and non-neural sources to the global signal remain poorly understood. This study sought to tackle this problem through the comparison of the BOLD global signal to an adjacent non-brain tissue signal, where neural activity was absent, from the same rs-fMRI scan obtained from anesthetized rats. In this dataset, motion was minimal and ventilation was phase-locked to image acquisition to minimize respiratory fluctuations. Data were acquired using three different anesthetics: isoflurane, dexmedetomidine, and a combination of dexmedetomidine and light isoflurane.

METHODS

A power spectral density estimate, a voxel-wise spatial correlation via Pearson's correlation, and a co-activation pattern analysis were performed using the global signal and the non-brain tissue signal. Functional connectivity was calculated using Pearson's linear correlation on default mode network (DMN) regions.

RESULTS

We report differences in the spectral composition of the two signals and show spatial selectivity within DMN structures that show an increased correlation to the global signal and decreased intra-network connectivity after global signal regression. All of the observed differences between the global signal and the non-brain tissue signal were maintained across anesthetics.

CONCLUSION

These results show that the global signal is distinct from the noise contained in the tissue signal, as support for a neural contribution. This study provides a unique perspective to the contents of the global signal and their origins.

Large-scale momentary brain co-activation patterns are associated with hyperalgesia and mediate focal neurochemistry and cross-network functional connectivity in fibromyalgia

ABSTRACT

Fibromyalgia has been characterized by augmented cross-network functional communication between the brain's sensorimotor, default mode, and attentional (salience/ventral and dorsal) networks. However, the underlying mechanisms of these aberrant communication patterns are unknown. In this study, we sought to understand large-scale topographic patterns at instantaneous timepoints, known as co-activation patterns (CAPs). We found that a sustained pressure pain challenge temporally modulated the occurrence of CAPs. Using proton magnetic resonance spectroscopy, we found that greater basal excitatory over inhibitory neurotransmitter levels within the anterior insula orchestrated higher cross-network connectivity between the anterior insula and the default mode network through lower occurrence of a CAP encompassing the attentional networks during sustained pain. Moreover, we found that hyperalgesia in fibromyalgia was mediated through increased occurrence of a CAP encompassing the sensorimotor network during sustained pain. In conclusion, this study elucidates the role of momentary large-scale topographic brain patterns in shaping noxious information in patients with fibromyalgia, while laying the groundwork for using precise spatiotemporal dynamics of the brain for the potential development of therapeutics.

Structural and functional network analysis of twins using fMRI data

Similarities between twins have been widely demonstrated, underscoring the remarkable influence of genetics across numerous traits. In this study, we explore the genetic underpinnings of the human brain by examining MRI data from the Queensland Twin Imaging study. Specifically, this study seeks to compare brain structure and function between twins and unrelated subjects, with an emphasis on describing the effects of genetic factors. To achieve these goals, we employed the source-based morphometry method to extract intrinsic components and elucidate recognizable patterns. Our results show that twins exhibit a higher degree of similarity in gray and white matter density compared with unrelated individuals. In addition, four distinct states of brain activity were identified using coactivation patterns analysis. Furthermore, twins demonstrated a greater degree of similarity in the temporal and spatial features of each state compared with unrelated subjects. Taken together, these results support the hypothesis that twins show greater similarity in both brain structure and dynamic functional brain activity. Further exploration of these methods may advance our understanding of the complex interplay between genes, environment, and brain networks.

Living on the edge: network neuroscience beyond nodes

Network neuroscience has emphasized the connectional properties of neural elements - cells, populations, and regions. This has come at the expense of the anatomical and functional connections that link these elements to one another. A new perspective - namely one that emphasizes 'edges' - may prove fruitful in addressing outstanding questions in network neuroscience. We highlight one recently proposed 'edge-centric' method and review its current applications, merits, and limitations. We also seek to establish conceptual and mathematical links between this method and previously proposed approaches in the network science and neuroimaging literature. We conclude by presenting several avenues for future work to extend and refine existing edge-centric analysis.

Recent advances in the neuroscience of spontaneous and off-task thought: implications for mental health

People spend a remarkable 30-50% of awake life thinking about something other than what they are currently doing. These experiences of being "off-task" can be described as spontaneous thought when mental dynamics are relatively flexible. Here we review recent neuroscience developments in this area and consider implications for mental wellbeing and illness. We provide updated overviews of the roles of the default mode network and large-scale network dynamics, and we discuss emerging candidate mechanisms involving hippocampal memory (sharp-wave ripples, replay) and neuromodulatory (noradrenergic and serotonergic) systems. We explore how distinct brain states can be associated with or give rise to adaptive and maladaptive forms of thought linked to distinguishable mental health outcomes. We conclude by outlining new directions in the neuroscience of spontaneous and off-task thought that may clarify mechanisms, lead to personalized biomarkers, and facilitate therapy developments toward the goals of better understanding and improving mental health.

Frequency dependent whole-brain coactivation patterns analysis in Alzheimer's disease

Background

The brain in resting state has complex dynamic properties and shows frequency dependent characteristics. The frequency-dependent whole-brain dynamic changes of resting state across the scans have been ignored in Alzheimer's disease (AD).

Objective

Coactivation pattern (CAP) analysis can identify different brain states. This paper aimed to investigate the dynamic characteristics of frequency dependent whole-brain CAPs in AD.

Methods

We utilized a multiband CAP approach to model the state space and study brain dynamics in both AD and NC. The correlation between the dynamic characteristics and the subjects' clinical index was further analyzed.

Results

The results showed similar CAP patterns at different frequency bands, but the occurrence of patterns was different. In addition, CAPs associated with the default mode network (DMN) and the ventral/dorsal visual network (dorsal/ventral VN) were altered significantly between the AD and NC groups. This study also found the correlation between the altered dynamic characteristics of frequency dependent CAPs and the patients' clinical Mini-Mental State Examination assessment scale scores.

Conclusion

This study revealed that while similar CAP spatial patterns appear in different frequency bands, their dynamic characteristics in subbands vary. In addition, delineating subbands was more helpful in distinguishing AD from NC in terms of CAP.

Deconstructing the Mapper algorithm to extract richer topological and temporal features from functional neuroimaging data

Capturing and tracking large-scale brain activity dynamics holds the potential to deepen our understanding of cognition. Previously, tools from Topological Data Analysis, especially Mapper, have been successfully used to mine brain activity dynamics at the highest spatiotemporal resolutions. Even though it is a relatively established tool within the field of Topological Data Analysis, Mapper results are highly impacted by parameter selection. Given that non-invasive human neuroimaging data (e.g., from fMRI) is typically fraught with artifacts and no gold standards exist regarding "true" state transitions, we argue for a thorough examination of Mapper parameter choices to better reveal their impact. Using synthetic data (with known transition structure) and real fMRI data, we explore a variety of parameter choices for each Mapper step, thereby providing guidance and heuristics for the field. We also release our parameter-exploration toolbox as a software package to make it easier for scientists to investigate and apply Mapper on any dataset.

Trait- and state-like co-activation pattern dynamics in current and remitted major depressive disorder

BACKGROUND

Distinguishing between trait- and state-like neural alternations in major depressive disorder (MDD) may advance our understanding of this recurring disorder. We aimed to investigate dynamic functional connectivity alternations in unmedicated individuals with current or past MDD using co-activation pattern analyses.

METHODS

Resting-state functional magnetic resonance imaging data were acquired from individuals with first-episode current MDD (cMDD, n = 50), remitted MDD (rMDD, n = 44), and healthy controls (HCs, n = 64). Using a data-driven consensus clustering technique, four whole-brain states of spatial co-activation were identified and associated metrics (dominance, entries, transition frequency) were analyzed with respect to clinical characteristics.

RESULTS

Relative to rMDD and HC, cMDD showed increased dominance and entries of state 1 (primarily involving default mode network (DMN)), and decreased dominance of state 4 (mostly involving frontal-parietal network (FPN)). Among cMDD, state 1 entries correlated positively with trait rumination. Conversely, relative to cMDD and HC, individuals with rMDD were characterized by increased state 4 entries. Relative to HC, both MDD groups showed increased state 4-to-1 (FPN to DMN) transition frequency but reduction in state 3 (spanning visual attention, somatosensory, limbic networks), with the former metric specifically related to trait rumination.

LIMITATIONS

Further confirmation with longitudinal studies are required.

CONCLUSIONS

Regardless of symptoms, MDD was characterized by increased FPN-to-DMN transitions and reduced dominance of a hybrid network. State-related effect emerged in regions critically implicated in repetitive introspection and cognitive control. Asymptomatic individuals with past MDD were uniquely linked to increased FPN entries. Our findings identify trait-like brain network dynamics that might increase vulnerability to future MDD.

Estimating dynamic individual coactivation patterns based on densely sampled resting-state fMRI data and utilizing it for better subject identification

As a complex dynamic system, the brain exhibits spatially organized recurring patterns of activity over time. Coactivation patterns (CAPs), which analyzes data from each single frame, have been utilized to detect transient brain activity states recently. However, previous CAP analyses have been conducted at the group level, which might neglect meaningful individual differences. Here, we estimated individual CAP states at both subject- and scan-level based on a densely sampled dataset: Midnight Scan Club. We used differential identifiability, which measures the gap between intra- and inter-subject similarity, to evaluate individual differences. We found individual CAPs at the subject-level achieved the best fingerprinting ability by maintaining high intra-subject similarity and enlarging inter-subject differences, and brain regions of association networks mainly contributed to the identifiability. On the other hand, scan-level CAP states were unstable across scans for the same participant. Expectedly, we found subject-specific CAPs became more reliable and discriminative with more data (i.e., longer duration). As the acquisition time of each participant is limited in practice, our results recommend a data collection strategy that collects more scans with appropriate duration (e.g., 12 ~ 15 min/scan) to obtain more reliable subject-specific CAPs, when total acquisition time is fixed (e.g., 150 min). In summary, this work has constructed reliable subject-specific CAP states with meaningful individual differences, and recommended an appropriate data collection strategy, which can guide subsequent investigations into individualized brain dynamics.

Efficient Inference of Spatially-Varying Gaussian Markov Random Fields With Applications in Gene Regulatory Networks

In this paper, we study the problem of inferring spatially-varying Gaussian Markov random fields (SV-GMRF) where the goal is to learn a network of sparse, context-specific GMRFs representing network relationships between genes. An important application of SV-GMRFs is in inference of gene regulatory networks from spatially-resolved transcriptomics datasets. The current work on inference of SV-GMRFs are based on the regularized maximum likelihood estimation (MLE) and suffer from overwhelmingly high computational cost due to their highly nonlinear nature. To alleviate this challenge, we propose a simple and efficient optimization problem in lieu of MLE that comes equipped with strong statistical and computational guarantees. Our proposed optimization problem is extremely efficient in practice: we can solve instances of SV-GMRFs with more than 2 million variables in less than 2 minutes. We apply the developed framework to study how gene regulatory networks in Glioblastoma are spatially rewired within tissue, and identify prominent activity of the transcription factor HES4 and ribosomal proteins as characterizing the gene expression network in the tumor peri-vascular niche that is known to harbor treatment resistant stem cells.