Edge-centric functional network representations of human cerebral cortex reveal overlapping system-level architecture

Smoking affects global and regional brain entropy in depression patients regardless of depression: Preliminary findings

OBJECTIVE

This study examines the effect of smoking on global and regional brain entropy in patients with Major Depressive Disorder (MDD), aiming to elucidate the relationship between smoking habits and brain network complexity in depression.

METHODS

The study enrolled 24 MDD patients, divided into smokers and non-smokers, from Alanya Alaaddin Keykubat University and Istanbul Medipol University. Resting-state fMRI data were acquired and processed. The complexity of neuronal activity was assessed using dispersion entropy, with statistical significance determined by a suite of tests including Kolmogorov-Smirnov, Student's t-test, and Mann-Whitney U test.

RESULTS

The smoking cohort exhibited higher global brain entropy compared to the non-smoking group (p = 0.033), with significant differences in various brain networks, indicating that smoking may alter global brain activity and network dynamics in individuals with MDD.

CONCLUSION

The study provides evidence that smoking is associated with increased brain entropy in MDD patients, suggesting that chronic smoking may influence cognitive and emotional networks. This underscores the importance of considering smoking history in the treatment and prognosis of MDD. The findings call for further research to understand the mechanistic links between smoking, brain entropy, and depression.

Network state dynamics underpin basal craving in a transdiagnostic population

Emerging fMRI methods quantifying brain dynamics present an opportunity to capture how fluctuations in brain responses give rise to individual variations in affective and motivation states. Although the experience and regulation of affective states affect psychopathology, their underlying time-varying brain responses remain unclear. Here, we present a novel framework to identify network states matched to an affective experience and examine how the dynamic engagement of these network states contributes to this experience. We apply this framework to investigate network state dynamics underlying basal craving, an affective experience with important clinical implications. In a transdiagnostic sample of healthy controls and individuals diagnosed with or at risk for craving-related disorders (total N = 252), we utilized connectome-based predictive modeling (CPM) to identify brain networks predictive of basal craving. An edge-centric timeseries approach was leveraged to quantify the moment-to-moment engagement of the craving-positive and craving-negative subnetworks during independent scan runs. We found that dynamic markers of network engagement, namely more persistence in a craving-positive network state and less dwelling in a craving-negative network state, characterized individuals with higher craving. We replicated the latter results in a separate dataset, incorporating distinct participants (N = 173) and experimental stimuli. The associations between basal craving and network state dynamics were consistently observed even when craving-predictive networks were defined in the replication dataset. These robust findings suggest that network state dynamics underpin individual differences in basal craving. Our framework additionally presents a new avenue to explore how the moment-to-moment engagement of behaviorally meaningful network states supports our affective experiences.

fMRI signals in white matter rewire gray matter community organization

Human brain gray matter (GM) has usually been clustered into multiple functional networks. The white matter (WM) fiber bundles are known to interconnect these networks simultaneously, engaging in numerous cognitive functions. However, the exact interconnections between GM and WM are still unclear, whether functional signals in WM rewires GM community organization remains to be explored. In this study, we divided brain functional connections into three types by using edge-centric method, including intra-GM, intra-WM and GM-WM connections, and calculated the edge community evaluation indexes for quantifying GM community engagement. The results showed that the involvement of WM significantly enhanced community entropy in the heteromodal system, while the sensory-attention system remained barely changed. In addition, delta community entropy showed a significant correlation with clinical cognitive scale. Our results suggested that WM rewired GM community organization, enhancing the community engagement of brain regions in the heteromodal system. This involvement was observed to be disrupted in disease groups. Our study revealed that considering the functional signals of GM and WM simultaneously could better understand the brain's functional organization.

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.

Abnormal Dynamic Reconstruction of Overlapping Communities in Schizophrenia Patients

OBJECTIVE

This study aims to explore the changes in dynamic overlapping communities in the brains of schizophrenia (SZ) patients and further investigate the dynamic restructuring patterns of overlapping communities in SZ patients.

MATERIALS AND METHODS

A total of 43 SZ patients and 49 normal controls (NC) were selected for resting-state functional MRI (rs-fMRI) scans. Dynamic functional connectivity analysis was conducted separately on SZ patients and NC using rs-fMRI and Jackknife Correlation techniques to construct dynamic brain network models. Based on these models, a dynamic overlapping community detection method was utilized to explore the abnormal overlapping community structure in SZ patients using evaluation metrics such as the structural stability of overlapping communities, nodes' functional diversity, and activity level of overlapping communities.

RESULTS

The stability of communities in SZ patients showed a decreasing trend. The changes in the overlapping community structure of SZ patients may be related to a decrease in the diversity of overlapping node functions. Additionally, compared to the NC group, the activity level of overlapping communities of SZ patients was significantly reduced.

CONCLUSION

The structure or organization of the brain functional network in SZ patients is abnormal or disrupted, and the activity of the brain network in information processing and transmission is weakened in SZ patients.

Macroscale brain states support the control of semantic cognition

A crucial aim in neuroscience is to understand how the human brain adapts to varying cognitive demands. This study investigates network reconfiguration during controlled semantic retrieval in differing contexts. We analyze brain responses to two semantic tasks of varying difficulty - global association and feature matching judgments - which are contrasted with non-semantic tasks on the cortical surface and within a whole-brain state space. Demanding semantic association tasks elicit activation in anterior prefrontal and temporal regions, while challenging semantic feature matching and non-semantic tasks predominantly activate posterior regions. Task difficulty also modulates activation along different dimensions of functional organization, suggesting different mechanisms of cognitive control. More demanding semantic association judgments engage cognitive control and default mode networks together, while feature matching and non-semantic tasks are skewed towards cognitive control networks. These findings highlight the brain's dynamic ability to tailor its networks to support diverse neurocognitive states, enriching our understanding of controlled cognition.

Brain Network Interconnectivity Dynamics Explain Metacognitive Differences in Listening Behavior

Complex auditory scenes pose a challenge to attentive listening, rendering listeners slower and more uncertain in their perceptual decisions. How can we explain such behaviors from the dynamics of cortical networks that pertain to the control of listening behavior? We here follow up on the hypothesis that human adaptive perception in challenging listening situations is supported by modular reconfiguration of auditory-control networks in a sample of N = 40 participants (13 males) who underwent resting-state and task functional magnetic resonance imaging (fMRI). Individual titration of a spatial selective auditory attention task maintained an average accuracy of ∼70% but yielded considerable interindividual differences in listeners' response speed and reported confidence in their own perceptual decisions. Whole-brain network modularity increased from rest to task by reconfiguring auditory, cinguloopercular, and dorsal attention networks. Specifically, interconnectivity between the auditory network and cinguloopercular network decreased during the task relative to the resting state. Additionally, interconnectivity between the dorsal attention network and cinguloopercular network increased. These interconnectivity dynamics were predictive of individual differences in response confidence, the degree of which was more pronounced after incorrect judgments. Our findings uncover the behavioral relevance of functional cross talk between auditory and attentional-control networks during metacognitive assessment of one's own perception in challenging listening situations and suggest two functionally dissociable cortical networked systems that shape the considerable metacognitive differences between individuals in adaptive listening behavior.

A synergistic workspace for human consciousness revealed by Integrated Information Decomposition

How is the information-processing architecture of the human brain organised, and how does its organisation support consciousness? Here, we combine network science and a rigorous information-theoretic notion of synergy to delineate a 'synergistic global workspace', comprising gateway regions that gather synergistic information from specialised modules across the human brain. This information is then integrated within the workspace and widely distributed via broadcaster regions. Through functional MRI analysis, we show that gateway regions of the synergistic workspace correspond to the human brain's default mode network, whereas broadcasters coincide with the executive control network. We find that loss of consciousness due to general anaesthesia or disorders of consciousness corresponds to diminished ability of the synergistic workspace to integrate information, which is restored upon recovery. Thus, loss of consciousness coincides with a breakdown of information integration within the synergistic workspace of the human brain. This work contributes to conceptual and empirical reconciliation between two prominent scientific theories of consciousness, the Global Neuronal Workspace and Integrated Information Theory, while also advancing our understanding of how the human brain supports consciousness through the synergistic integration of information.

Lack of evidence for predictive utility from resting state fMRI data for individual exposure-based cognitive behavioral therapy outcomes: A machine learning study in two large multi-site samples in anxiety disorders

Data-based predictions of individual Cognitive Behavioral Therapy (CBT) treatment response are a fundamental step towards precision medicine. Past studies demonstrated only moderate prediction accuracy (i.e. ability to discriminate between responders and non-responders of a given treatment) when using clinical routine data such as demographic and questionnaire data, while neuroimaging data achieved superior prediction accuracy. However, these studies may be considerably biased due to very limited sample sizes and bias-prone methodology. Adequately powered and cross-validated samples are a prerequisite to evaluate predictive performance and to identify the most promising predictors. We therefore analyzed resting state functional magnet resonance imaging (rs-fMRI) data from two large clinical trials to test whether functional neuroimaging data continues to provide good prediction accuracy in much larger samples. Data came from two distinct German multicenter studies on exposure-based CBT for anxiety disorders, the Protect-AD and SpiderVR studies. We separately and independently preprocessed baseline rs-fMRI data from n = 220 patients (Protect-AD) and n = 190 patients (SpiderVR) and extracted a variety of features, including ROI-to-ROI and edge-functional connectivity, sliding-windows, and graph measures. Including these features in sophisticated machine learning pipelines, we found that predictions of individual outcomes never significantly differed from chance level, even when conducting a range of exploratory post-hoc analyses. Moreover, resting state data never provided prediction accuracy beyond the sociodemographic and clinical data. The analyses were independent of each other in terms of selecting methods to process resting state data for prediction input as well as in the used parameters of the machine learning pipelines, corroborating the external validity of the results. These similar findings in two independent studies, analyzed separately, urge caution regarding the interpretation of promising prediction results based on neuroimaging data from small samples and emphasizes that some of the prediction accuracies from previous studies may result from overestimation due to homogeneous data and weak cross-validation schemes. The promise of resting-state neuroimaging data to play an important role in the prediction of CBT treatment outcomes in patients with anxiety disorders remains yet to be delivered.

Reliability of variability and complexity measures for task and task-free BOLD fMRI

Brain activity continuously fluctuates over time, even if the brain is in controlled (e.g., experimentally induced) states. Recent years have seen an increasing interest in understanding the complexity of these temporal variations, for example with respect to developmental changes in brain function or between-person differences in healthy and clinical populations. However, the psychometric reliability of brain signal variability and complexity measures-which is an important precondition for robust individual differences as well as longitudinal research-is not yet sufficiently studied. We examined reliability (split-half correlations) and test-retest correlations for task-free (resting-state) BOLD fMRI as well as split-half correlations for seven functional task data sets from the Human Connectome Project to evaluate their reliability. We observed good to excellent split-half reliability for temporal variability measures derived from rest and task fMRI activation time series (standard deviation, mean absolute successive difference, mean squared successive difference), and moderate test-retest correlations for the same variability measures under rest conditions. Brain signal complexity estimates (several entropy and dimensionality measures) showed moderate to good reliabilities under both, rest and task activation conditions. We calculated the same measures also for time-resolved (dynamic) functional connectivity time series and observed moderate to good reliabilities for variability measures, but poor reliabilities for complexity measures derived from functional connectivity time series. Global (i.e., mean across cortical regions) measures tended to show higher reliability than region-specific variability or complexity estimates. Larger subcortical regions showed similar reliability as cortical regions, but small regions showed lower reliability, especially for complexity measures. Lastly, we also show that reliability scores are only minorly dependent on differences in scan length and replicate our results across different parcellation and denoising strategies. These results suggest that the variability and complexity of BOLD activation time series are robust measures well-suited for individual differences research. Temporal variability of global functional connectivity over time provides an important novel approach to robustly quantifying the dynamics of brain function. PRACTITIONER POINTS: Variability and complexity measures of BOLD activation show good split-half reliability and moderate test-retest reliability. Measures of variability of global functional connectivity over time can robustly quantify neural dynamics. Length of fMRI data has only a minor effect on reliability.

Unravelling consciousness and brain function through the lens of time, space, and information

Disentangling how cognitive functions emerge from the interplay of brain dynamics and network architecture is among the major challenges that neuroscientists face. Pharmacological and pathological perturbations of consciousness provide a lens to investigate these complex challenges. Here, we review how recent advances about consciousness and the brain's functional organisation have been driven by a common denominator: decomposing brain function into fundamental constituents of time, space, and information. Whereas unconsciousness increases structure-function coupling across scales, psychedelics may decouple brain function from structure. Convergent effects also emerge: anaesthetics, psychedelics, and disorders of consciousness can exhibit similar reconfigurations of the brain's unimodal-transmodal functional axis. Decomposition approaches reveal the potential to translate discoveries across species, with computational modelling providing a path towards mechanistic integration.

EEG decoding with spatiotemporal convolutional neural network for visualization and closed-loop control of sensorimotor activities: A simultaneous EEG-fMRI study

Closed-loop neurofeedback training utilizes neural signals such as scalp electroencephalograms (EEG) to manipulate specific neural activities and the associated behavioral performance. A spatiotemporal filter for high-density whole-head scalp EEG using a convolutional neural network can overcome the ambiguity of the signaling source because each EEG signal includes information on the remote regions. We simultaneously acquired EEG and functional magnetic resonance images in humans during the brain-computer interface (BCI) based neurofeedback training and compared the reconstructed and modeled hemodynamic responses of the sensorimotor network. Filters constructed with a convolutional neural network captured activities in the targeted network with spatial precision and specificity superior to those of the EEG signals preprocessed with standard pipelines used in BCI-based neurofeedback paradigms. The middle layers of the trained model were examined to characterize the neuronal oscillatory features that contributed to the reconstruction. Analysis of the layers for spatial convolution revealed the contribution of distributed cortical circuitries to reconstruction, including the frontoparietal and sensorimotor areas, and those of temporal convolution layers that successfully reconstructed the hemodynamic response function. Employing a spatiotemporal filter and leveraging the electrophysiological signatures of the sensorimotor excitability identified in our middle layer analysis would contribute to the development of a further effective neurofeedback intervention.

Tau follows principal axes of functional and structural brain organization in Alzheimer's disease

Alzheimer's disease (AD) is a brain network disorder where pathological proteins accumulate through networks and drive cognitive decline. Yet, the role of network connectivity in facilitating this accumulation remains unclear. Using in-vivo multimodal imaging, we show that the distribution of tau and reactive microglia in humans follows spatial patterns of connectivity variation, the so-called gradients of brain organization. Notably, less distinct connectivity patterns ("gradient contraction") are associated with cognitive decline in regions with greater tau, suggesting an interaction between reduced network differentiation and tau on cognition. Furthermore, by modeling tau in subject-specific gradient space, we demonstrate that tau accumulation in the frontoparietal and temporo-occipital cortices is associated with greater baseline tau within their functionally and structurally connected hubs, respectively. Our work unveils a role for both functional and structural brain organization in pathology accumulation in AD, and supports subject-specific gradient space as a promising tool to map disease progression.

Systematic evaluation of fMRI data-processing pipelines for consistent functional connectomics

Functional interactions between brain regions can be viewed as a network, enabling neuroscientists to investigate brain function through network science. Here, we systematically evaluate 768 data-processing pipelines for network reconstruction from resting-state functional MRI, evaluating the effect of brain parcellation, connectivity definition, and global signal regression. Our criteria seek pipelines that minimise motion confounds and spurious test-retest discrepancies of network topology, while being sensitive to both inter-subject differences and experimental effects of interest. We reveal vast and systematic variability across pipelines' suitability for functional connectomics. Inappropriate choice of data-processing pipeline can produce results that are not only misleading, but systematically so, with the majority of pipelines failing at least one criterion. However, a set of optimal pipelines consistently satisfy all criteria across different datasets, spanning minutes, weeks, and months. We provide a full breakdown of each pipeline's performance across criteria and datasets, to inform future best practices in functional connectomics.

YAP in development and disease: Navigating the regulatory landscape from retina to brain

The distinctive role of Yes-associated protein (YAP) in the nervous system has attracted widespread attention. This comprehensive review strategically uses the retina as a vantage point, embarking on an extensive exploration of YAP's multifaceted impact from the retina to the brain in development and pathology. Initially, we explore the crucial roles of YAP in embryonic and cerebral development. Our focus then shifts to retinal development, examining in detail YAP's regulatory influence on the development of retinal pigment epithelium (RPE) and retinal progenitor cells (RPCs), and its significant effects on the hierarchical structure and functionality of the retina. We also investigate the essential contributions of YAP in maintaining retinal homeostasis, highlighting its precise regulation of retinal cell proliferation and survival. In terms of retinal-related diseases, we explore the epigenetic connections and pathophysiological regulation of YAP in diabetic retinopathy (DR), glaucoma, and proliferative vitreoretinopathy (PVR). Lastly, we broaden our exploration from the retina to the brain, emphasizing the research paradigm of "retina: a window to the brain." Special focus is given to the emerging studies on YAP in brain disorders such as Alzheimer's disease (AD) and Parkinson's disease (PD), underlining its potential therapeutic value in neurodegenerative disorders and neuroinflammation.

The Brain's Topographical Organization Shapes Dynamic Interaction Patterns That Support Flexible Behavior Based on Rules and Long-Term Knowledge

Adaptive behavior relies both on specific rules that vary across situations and stable long-term knowledge gained from experience. The frontoparietal control network (FPCN) is implicated in the brain's ability to balance these different influences on action. Here, we investigate how the topographical organization of the cortex supports behavioral flexibility within the FPCN. Functional properties of this network might reflect its juxtaposition between the dorsal attention network (DAN) and the default mode network (DMN), two large-scale systems implicated in top-down attention and memory-guided cognition, respectively. Our study tests whether subnetworks of FPCN are topographically proximal to the DAN and the DMN, respectively, and how these topographical differences relate to functional differences: the proximity of each subnetwork is anticipated to play a pivotal role in generating distinct cognitive modes relevant to working memory and long-term memory. We show that FPCN subsystems share multiple anatomical and functional similarities with their neighboring systems (DAN and DMN) and that this topographical architecture supports distinct interaction patterns that give rise to different patterns of functional behavior. The FPCN acts as a unified system when long-term knowledge supports behavior but becomes segregated into discrete subsystems with different patterns of interaction when long-term memory is less relevant. In this way, our study suggests that the topographical organization of the FPCN and the connections it forms with distant regions of cortex are important influences on how this system supports flexible behavior.

Inferring brain causal and temporal-lag networks for recognizing abnormal patterns of dementia

Brain functional network analysis has become a popular method to explore the laws of brain organization and identify biomarkers of neurological diseases. However, it is still a challenging task to construct an ideal brain network due to the limited understanding of the human brain. Existing methods often ignore the impact of temporal-lag on the results of brain network modeling, which may lead to some unreliable conclusions. To overcome this issue, we propose a novel brain functional network estimation method, which can simultaneously infer the causal mechanisms and temporal-lag values among brain regions. Specifically, our method converts the lag learning into an instantaneous effect estimation problem, and further embeds the search objectives into a deep neural network model as parameters to be learned. To verify the effectiveness of the proposed estimation method, we perform experiments on the Alzheimer's Disease Neuroimaging Initiative (ADNI) database by comparing the proposed model with several existing methods, including correlation-based and causality-based methods. The experimental results show that our brain networks constructed by the proposed estimation method can not only achieve promising classification performance, but also exhibit some characteristics of physiological mechanisms. Our approach provides a new perspective for understanding the pathogenesis of brain diseases. The source code is released at https://github.com/NJUSTxiazw/CTLN.

Unveiling the Core Functional Networks of Cognition: An Ontology-Guided Machine Learning Approach

Deciphering the functional architecture that underpins diverse cognitive functions is fundamental quest in neuroscience. In this study, we employed an innovative machine learning framework that integrated cognitive ontology with functional connectivity analysis to identify brain networks essential for cognition. We identified a core assembly of functional connectomes, primarily located within the association cortex, which showed superior predictive performance compared to two conventional methods widely employed in previous research across various cognitive domains. Our approach achieved a mean prediction accuracy of 0.13 across 16 cognitive tasks, including working memory, reading comprehension, and sustained attention, outperforming the traditional methods' accuracy of 0.08. In contrast, our method showed limited predictive power for sensory, motor, and emotional functions, with a mean prediction accuracy of 0.03 across 9 relevant tasks, slightly lower than the traditional methods' accuracy of 0.04. These cognitive connectomes were further characterized by distinctive patterns of resting-state functional connectivity, structural connectivity via white matter tracts, and gene expression, highlighting their neurogenetic underpinnings. Our findings reveal a domain-general functional network fingerprint that pivotal to cognition, offering a novel computational approach to explore the neural foundations of cognitive abilities.

Edge-Based General Linear Models Capture Moment-to-Moment Fluctuations in Attention

Although we must prioritize the processing of task-relevant information to navigate life, our ability to do so fluctuates across time. Previous work has identified fMRI functional connectivity (FC) networks that predict an individual's ability to sustain attention and vary with attentional state from 1 min to the next. However, traditional dynamic FC approaches typically lack the temporal precision to capture moment-to-moment network fluctuations. Recently, researchers have "unfurled" traditional FC matrices in "edge cofluctuation time series" which measure timepoint-by-timepoint cofluctuations between regions. Here we apply event-based and parametric fMRI analyses to edge time series to capture moment-to-moment fluctuations in networks related to attention. In two independent fMRI datasets examining young adults of both sexes in which participants performed a sustained attention task, we identified a reliable set of edges that rapidly deflects in response to rare task events. Another set of edges varies with continuous fluctuations in attention and overlaps with a previously defined set of edges associated with individual differences in sustained attention. Demonstrating that edge-based analyses are not simply redundant with traditional regions-of-interest-based approaches, up to one-third of reliably deflected edges were not predicted from univariate activity patterns alone. These results reveal the large potential in combining traditional fMRI analyses with edge time series to identify rapid reconfigurations in networks across the brain.

Similarity in evoked responses does not imply similarity in macroscopic network states

It is commonplace in neuroscience to assume that if two tasks activate the same brain areas in the same way, then they are recruiting the same underlying networks. Yet computational theory has shown that the same pattern of activity can emerge from many different underlying network representations. Here we evaluated whether similarity in activation necessarily implies similarity in network architecture by comparing region-wise activation patterns and functional correlation profiles from a large sample of healthy subjects (N = 242). Participants performed two executive control tasks known to recruit nearly identical brain areas, the color-word Stroop task and the Multi-Source Interference Task (MSIT). Using a measure of instantaneous functional correlations, based on edge time series, we estimated the task-related networks that differed between incongruent and congruent conditions. We found that the two tasks were much more different in their network profiles than in their evoked activity patterns at different analytical levels, as well as for a wide range of methodological pipelines. Our results reject the notion that having the same activation patterns means two tasks engage the same underlying representations, suggesting that task representations should be independently evaluated at both node and edge (connectivity) levels.

Network-wise analysis of movie-specific information in dynamic functional connectivity using COBE

Movie watching during functional magnetic resonance imaging has emerged as a promising tool to measure the complex behavior of the brain in response to a stimulus similar to real-life situations. It has been observed that presenting a movie (sequence of events) as a stimulus will lead to a unique time course of dynamic functional connectivity related to movie stimuli that can be compared across the participants. We assume that the observed dynamic functional connectivity across subjects can be divided into following 2 components: (i) specific to a movie stimulus (depicting group-level behavior in functional connectivity) and (ii) individual-specific behavior (not necessarily common across the subjects). In this work, using the dynamic time warping distance measure, we have shown the extent of similarity between the temporal sequences of functional connectivity while the underlying movie stimuli were same and different. Further, the temporal sequence of functional connectivity patterns related to a movie is enhanced by suppressing the subject-specific components of dynamic functional connectivity using common and orthogonal basis extraction. Quantitative analysis using the F-ratio measure reveals significant differences in dynamic functional connectivity within the somatomotor network and default mode network, as well as between the occipital network and somatomotor networks.

Higher-order functional connectivity analysis of resting-state functional magnetic resonance imaging data using multivariate cumulants

Blood-level oxygenation-dependent (BOLD) functional magnetic resonance imaging (fMRI) is the most common modality to study functional connectivity in the human brain. Most research to date has focused on connectivity between pairs of brain regions. However, attention has recently turned towards connectivity involving more than two regions, that is, higher-order connectivity. It is not yet clear how higher-order connectivity can best be quantified. The measures that are currently in use cannot distinguish between pairwise (i.e., second-order) and higher-order connectivity. We show that genuine higher-order connectivity can be quantified by using multivariate cumulants. We explore the use of multivariate cumulants for quantifying higher-order connectivity and the performance of block bootstrapping for statistical inference. In particular, we formulate a generative model for fMRI signals exhibiting higher-order connectivity and use it to assess bias, standard errors, and detection probabilities. Application to resting-state fMRI data from the Human Connectome Project demonstrates that spontaneous fMRI signals are organized into higher-order networks that are distinct from second-order resting-state networks. Application to a clinical cohort of patients with multiple sclerosis further demonstrates that cumulants can be used to classify disease groups and explain behavioral variability. Hence, we present a novel framework to reliably estimate genuine higher-order connectivity in fMRI data which can be used for constructing hyperedges, and finally, which can readily be applied to fMRI data from populations with neuropsychiatric disease or cognitive neuroscientific experiments.

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.

Group-common and individual-specific effects of structure-function coupling in human brain networks with graph neural networks

The human cerebral cortex is organized into functionally segregated but synchronized regions bridged by the structural connectivity of white matter pathways. While structure-function coupling has been implicated in cognitive development and neuropsychiatric disorders, studies yield inconsistent findings. The extent to which the structure-function coupling reflects reliable individual differences or primarily group-common characteristics remains unclear, at both the global and regional brain levels. By leveraging two independent, high-quality datasets, we found that the graph neural network accurately predicted unseen individuals' functional connectivity from structural connectivity, reflecting a strong structure-function coupling. This coupling was primarily driven by network topology and was substantially stronger than that of the linear models. Moreover, we observed that structure-function coupling was dominated by group-common effects, with subtle yet significant individual-specific effects. The regional group and individual effects of coupling were hierarchically organized across the cortex along a sensorimotor-association axis, with lower group and higher individual effects in association cortices. These findings emphasize the importance of considering both group and individual effects in understanding cortical structure-function coupling, suggesting insights into interpreting individual differences of the coupling and informing connectivity-guided therapeutics.

Personalized functional network mapping for autism spectrum disorder and attention-deficit/hyperactivity disorder

Autism spectrum disorder (ASD) and Attention-deficit/hyperactivity disorder (ADHD) are two typical neurodevelopmental disorders that have a long-term impact on physical and mental health. ASD is usually comorbid with ADHD and thus shares highly overlapping clinical symptoms. Delineating the shared and distinct neurophysiological profiles is important to uncover the neurobiological mechanisms to guide better therapy. In this study, we aimed to establish the behaviors, functional connectome, and network properties differences between ASD, ADHD-Combined, and ADHD-Inattentive using resting-state functional magnetic resonance imaging. We used the non-negative matrix fraction method to define personalized large-scale functional networks for each participant. The individual large-scale functional network connectivity (FNC) and graph-theory-based complex network analyses were executed and identified shared and disorder-specific differences in FNCs and network attributes. In addition, edge-wise functional connectivity analysis revealed abnormal edge co-fluctuation amplitude and number of transitions among different groups. Taken together, our study revealed disorder-specific and -shared regional and edge-wise functional connectivity and network differences for ASD and ADHD using an individual-level functional network mapping approach, which provides new evidence for the brain functional abnormalities in ASD and ADHD and facilitates understanding the neurobiological basis for both disorders.

Edge time series components of functional connectivity and cognitive function in Alzheimer's disease

Understanding the interrelationships of brain function as measured by resting-state magnetic resonance imaging and neuropsychological/behavioral measures in Alzheimer's disease is key for advancement of neuroimaging analysis methods in clinical research. The edge time-series framework recently developed in the field of network neuroscience, in combination with other network science methods, allows for investigations of brain-behavior relationships that are not possible with conventional functional connectivity methods. Data from the Indiana Alzheimer's Disease Research Center sample (53 cognitively normal control, 47 subjective cognitive decline, 32 mild cognitive impairment, and 20 Alzheimer's disease participants) were used to investigate relationships between functional connectivity components, each derived from a subset of time points based on co-fluctuation of regional signals, and measures of domain-specific neuropsychological functions. Multiple relationships were identified with the component approach that were not found with conventional functional connectivity. These involved attentional, limbic, frontoparietal, and default mode systems and their interactions, which were shown to couple with cognitive, executive, language, and attention neuropsychological domains. Additionally, overlapping results were obtained with two different statistical strategies (network contingency correlation analysis and network-based statistics correlation). Results demonstrate that connectivity components derived from edge time-series based on co-fluctuation reveal disease-relevant relationships not observed with conventional static functional connectivity.

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.

Multimodal measures of spontaneous brain activity reveal both common and divergent patterns of cortical functional organization

Large-scale functional networks have been characterized in both rodent and human brains, typically by analyzing fMRI-BOLD signals. However, the relationship between fMRI-BOLD and underlying neural activity is complex and incompletely understood, which poses challenges to interpreting network organization obtained using this technique. Additionally, most work has assumed a disjoint functional network organization (i.e., brain regions belong to one and only one network). Here, we employ wide-field Ca2+ imaging simultaneously with fMRI-BOLD in mice expressing GCaMP6f in excitatory neurons. We determine cortical networks discovered by each modality using a mixed-membership algorithm to test the hypothesis that functional networks exhibit overlapping organization. We find that there is considerable network overlap (both modalities) in addition to disjoint organization. Our results show that multiple BOLD networks are detected via Ca2+ signals, and networks determined by low-frequency Ca2+ signals are only modestly more similar to BOLD networks. In addition, the principal gradient of functional connectivity is nearly identical for BOLD and Ca2+ signals. Despite similarities, important differences are also detected across modalities, such as in measures of functional connectivity strength and diversity. In conclusion, Ca2+ imaging uncovers overlapping functional cortical organization in the mouse that reflects several, but not all, properties observed with fMRI-BOLD signals.

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.

Whole-brain multivariate hemodynamic deconvolution for functional MRI with stability selection

Conventionally, analysis of functional MRI (fMRI) data relies on available information about the experimental paradigm to establish hypothesized models of brain activity. However, this information can be inaccurate, incomplete or unavailable in multiple scenarios such as resting-state, naturalistic paradigms or clinical conditions. In these cases, blind estimates of neuronal-related activity can be obtained with paradigm-free analysis methods such as hemodynamic deconvolution. Yet, current formulations of the hemodynamic deconvolution problem have three important limitations: (1) their efficacy strongly depends on the appropriate selection of regularization parameters, (2) being univariate, they do not take advantage of the information present across the brain, and (3) they do not provide any measure of statistical certainty associated with each detected event. Here we propose a novel approach that addresses all these limitations. Specifically, we introduce multivariate sparse paradigm free mapping (Mv-SPFM), a novel hemodynamic deconvolution algorithm that operates at the whole brain level and adds spatial information via a mixed-norm regularization term over all voxels. Additionally, Mv-SPFM employs a stability selection procedure that removes the need to select regularization parameters and also lets us obtain an estimate of the true probability of having a neuronal-related BOLD event at each voxel and time-point based on the area under the curve (AUC) of the stability paths. Besides, we present a formulation tailored for multi-echo fMRI acquisitions (MvME-SPFM), which allows us to better isolate fluctuations of BOLD origin on the basis of their linear dependence with the echo time (TE) and to assign physiologically interpretable units (i.e., changes in the apparent transverse relaxation ΔR2∗) to the resulting deconvolved events. Remarkably, we demonstrate that Mv-SPFM achieves comparable performance even when using a single-echo formulation. We demonstrate that this algorithm outperforms existing state-of-the-art deconvolution approaches, and shows higher spatial and temporal agreement with the activation maps and BOLD signals obtained with a standard model-based linear regression approach, even at the level of individual neuronal events. Furthermore, we show that by employing stability selection, the performance of the algorithm depends less on the selection of temporal and spatial regularization parameters λ and ρ. Consequently, the proposed algorithm provides more reliable estimates of neuronal-related activity, here in terms of ΔR2∗, for the study of the dynamics of brain activity when no information about the timings of the BOLD events is available. This algorithm will be made publicly available as part of the splora Python package.

Functional hierarchies in brain dynamics characterized by signal reversibility in ferret cortex

Brain signal irreversibility has been shown to be a promising approach to study neural dynamics. Nevertheless, the relation with cortical hierarchy and the influence of different electrophysiological features is not completely understood. In this study, we recorded local field potentials (LFPs) during spontaneous behavior, including awake and sleep periods, using custom micro-electrocorticographic (μECoG) arrays implanted in ferrets. In contrast to humans, ferrets remain less time in each state across the sleep-wake cycle. We deployed a diverse set of metrics in order to measure the levels of complexity of the different behavioral states. In particular, brain irreversibility, which is a signature of non-equilibrium dynamics, captured by the arrow of time of the signal, revealed the hierarchical organization of the ferret's cortex. We found different signatures of irreversibility and functional hierarchy of large-scale dynamics in three different brain states (active awake, quiet awake, and deep sleep), showing a lower level of irreversibility in the deep sleep stage, compared to the other. Irreversibility also allowed us to disentangle the influence of different cortical areas and frequency bands in this process, showing a predominance of the parietal cortex and the theta band. Furthermore, when inspecting the embedded dynamic through a Hidden Markov Model, the deep sleep stage was revealed to have a lower switching rate and lower entropy production. These results suggest functional hierarchies in organization that can be revealed through thermodynamic features and information theory metrics.

Connectomes for 40,000 UK Biobank participants: A multi-modal, multi-scale brain network resource

We mapped functional and structural brain networks for more than 40,000 UK Biobank participants. Structural connectivity was estimated with tractography and diffusion MRI. Resting-state functional MRI was used to infer regional functional connectivity. We provide high-quality structural and functional connectomes for multiple parcellation granularities, several alternative measures of interregional connectivity, and a variety of common data pre-processing techniques, yielding more than one million connectomes in total and requiring more than 200,000 h of compute time. For a single subject, we provide 28 out-of-the-box versions of structural and functional brain networks, allowing users to select, e.g., the parcellation and connectivity measure that best suit their research goals. Furthermore, we provide code and intermediate data for the time-efficient reconstruction of more than 1000 different versions of a subject's connectome based on an array of methodological choices. All connectomes are available via the UK Biobank data-sharing platform and our connectome mapping pipelines are openly available. In this report, we describe our connectome resource in detail for users, outline key considerations in developing an efficient pipeline to map an unprecedented number of connectomes, and report on the quality control procedures that were completed to ensure connectome reliability and accuracy. We demonstrate that our structural and functional connectivity matrices meet a number of quality control checks and replicate previously established findings in network neuroscience. We envisage that our resource will enable new studies of the human connectome in health, disease, and aging at an unprecedented scale.

Abnormal Thalamo-Cortical Interactions in Overlapping Communities of Migraine: An Edge Functional Connectivity Study

OBJECTIVE

Migraine has been demonstrated to exhibit abnormal functional connectivity of large-scale brain networks, which is closely associated with its pathophysiology and has not yet been explored by edge functional connectivity. We used an edge-centric approach combined with motif analysis to evaluate higher-order communication patterns of brain networks in migraine.

METHODS

We investigated edge-centric metrics in 108 interictal migraine patients and 71 healthy controls. We parcellated the brain into networks using independent component analysis. We applied edge graph construction, k-means clustering, community overlap detection, graph-theory-based evaluations, and clinical correlation analysis. We conducted motif analysis to explore the interactions among regions, and a classification model to test the specificity of edge-centric results.

RESULTS

The normalized entropy of lateral thalamus was significantly increased in migraine, which was positively correlated with the baseline headache duration, and negatively correlated with headache duration reduction following preventive medications at 3-month follow-up. Network-wise entropy of the sensorimotor network was significantly elevated in migraine. The community similarity between lateral thalamus and postcentral gyrus was enhanced in migraine. Migraine patients showed overrepresented L-shape and diverse motifs, and underrepresented forked motifs with lateral thalamus serving as the reference node. Furthermore, migraine patients presented with overrepresented L-shape triads, where the postcentral gyrus shared different edges with the lateral thalamus. The classification model showed that entropy of the lateral thalamus had the highest discriminative power, with an area under the curve of 0.86.

INTERPRETATION

Our findings indicated an abnormal higher-order thalamo-cortical communication pattern in migraine patients. The thalamo-cortical-somatosensory disturbance of concerted working may potentially lead to aberrant information flow and deficit pain processing of migraine. ANN NEUROL 2023;94:1168-1181.

The virtual aging brain: Causal inference supports interhemispheric dedifferentiation in healthy aging

The mechanisms of cognitive decline and its variability during healthy aging are not fully understood, but have been associated with reorganization of white matter tracts and functional brain networks. Here, we built a brain network modeling framework to infer the causal link between structural connectivity and functional architecture and the consequent cognitive decline in aging. By applying in-silico interhemispheric degradation of structural connectivity, we reproduced the process of functional dedifferentiation during aging. Thereby, we found the global modulation of brain dynamics by structural connectivity to increase with age, which was steeper in older adults with poor cognitive performance. We validated our causal hypothesis via a deep-learning Bayesian approach. Our results might be the first mechanistic demonstration of dedifferentiation during aging leading to cognitive decline.

Edge Time Series Components of Functional Connectivity and Cognitive Function in Alzheimer's Disease

Understanding the interrelationships of brain function as measured by resting-state magnetic resonance imaging and neuropsychological/behavioral measures in Alzheimer's disease is key for advancement of neuroimaging analysis methods in clinical research. The edge time-series framework recently developed in the field of network neuroscience, in combination with other network science methods, allows for investigations of brain-behavior relationships that are not possible with conventional functional connectivity methods. Data from the Indiana Alzheimer's Disease Research Center sample (53 cognitively normal control, 47 subjective cognitive decline, 32 mild cognitive impairment, and 20 Alzheimer's disease participants) were used to investigate relationships between functional connectivity components, each derived from a subset of time points based on co-fluctuation of regional signals, and measures of domain-specific neuropsychological functions. Multiple relationships were identified with the component approach that were not found with conventional functional connectivity. These involved attentional, limbic, frontoparietal, and default mode systems and their interactions, which were shown to couple with cognitive, executive, language, and attention neuropsychological domains. Additionally, overlapping results were obtained with two different statistical strategies (network contingency correlation analysis and network-based statistics correlation). Results demonstrate that connectivity components derived from edge time-series based on co-fluctuation reveal disease-relevant relationships not observed with conventional static functional connectivity.

CAMBA framework: Unveiling the brain asymmetry alterations and longitudinal changes after stroke using resting-state EEG

Hemispheric asymmetry or lateralization is a fundamental principle of brain organization. However, it is poorly understood to what extent the brain asymmetries across different levels of functional organizations are evident in health or altered in brain diseases. Here, we propose a framework that integrates three degrees of brain interactions (isolated nodes, node-node, and edge-edge) into a unified analysis pipeline to capture the sliding window-based asymmetry dynamics at both the node and hemisphere levels. We apply this framework to resting-state EEG in healthy and stroke populations and investigate the stroke-induced abnormal alterations in brain asymmetries and longitudinal asymmetry changes during poststroke rehabilitation. We observe that the mean asymmetry in patients was abnormally enhanced across different frequency bands and levels of brain interactions, with these abnormal patterns strongly associated with the side of the stroke lesion. Compared to healthy controls, patients displayed significant alterations in asymmetry fluctuations, disrupting and reconfiguring the balance of inter-hemispheric integration and segregation. Additionally, analyses reveal that specific abnormal asymmetry metrics in patients tend to move towards those observed in healthy controls after short-term brain-computer interface rehabilitation. Furthermore, preliminary evidence suggests that baseline clinical and asymmetry features can predict poststroke improvements in the Fugl-Meyer assessment of the lower extremity (mean absolute error of about 2). Overall, these findings advance our understanding of hemispheric asymmetry. Our framework offers new insights into the mechanisms underlying brain alterations and recovery after a brain lesion, may help identify prognostic biomarkers, and can be easily extended to different functional modalities.

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.

A new causal centrality measure reveals the prominent role of subcortical structures in the causal architecture of the extended default mode network

Network representation has been an incredibly useful concept for understanding the behavior of complex systems in social sciences, biology, neuroscience, and beyond. Network science is mathematically founded on graph theory, where nodal importance is gauged using measures of centrality. Notably, recent work suggests that the topological centrality of a node should not be over-interpreted as its dynamical or causal importance in the network. Hence, identifying the influential nodes in dynamic causal models (DCM) remains an open question. This paper introduces causal centrality for DCM, a dynamics-sensitive and causally-founded centrality measure based on the notion of intervention in graphical models. Operationally, this measure simplifies to an identifiable expression using Bayesian model reduction. As a proof of concept, the average DCM of the extended default mode network (eDMN) was computed in 74 healthy subjects. Next, causal centralities of different regions were computed for this causal graph, and compared against several graph-theoretical centralities. The results showed that the subcortical structures of the eDMN were more causally central than the cortical regions, even though the graph-theoretical centralities unanimously favored the latter. Importantly, model comparison revealed that only the pattern of causal centrality was causally relevant. These results are consistent with the crucial role of the subcortical structures in the neuromodulatory systems of the brain, and highlight their contribution to the organization of large-scale networks. Potential applications of causal centrality-to study causal models of other neurotypical and pathological functional networks-are discussed, and some future lines of research are outlined.

Confounding Effects on the Performance of Machine Learning Analysis of Static Functional Connectivity Computed from rs-fMRI Multi-site Data

Resting-state functional magnetic resonance imaging (rs-fMRI) is a non-invasive imaging technique widely used in neuroscience to understand the functional connectivity of the human brain. While rs-fMRI multi-site data can help to understand the inner working of the brain, the data acquisition and processing of this data has many challenges. One of the challenges is the variability of the data associated with different acquisitions sites, and different MRI machines vendors. Other factors such as population heterogeneity among different sites, with variables such as age and gender of the subjects, must also be considered. Given that most of the machine-learning models are developed using these rs-fMRI multi-site data sets, the intrinsic confounding effects can adversely affect the generalizability and reliability of these computational methods, as well as the imposition of upper limits on the classification scores. This work aims to identify the phenotypic and imaging variables producing the confounding effects, as well as to control these effects. Our goal is to maximize the classification scores obtained from the machine learning analysis of the Autism Brain Imaging Data Exchange (ABIDE) rs-fMRI multi-site data. To achieve this goal, we propose novel methods of stratification to produce homogeneous sub-samples of the 17 ABIDE sites, as well as the generation of new features from the static functional connectivity values, using multiple linear regression models, ComBat harmonization models, and normalization methods. The main results obtained with our statistical models and methods are an accuracy of 76.4%, sensitivity of 82.9%, and specificity of 77.0%, which are 8.8%, 20.5%, and 7.5% above the baseline classification scores obtained from the machine learning analysis of the static functional connectivity computed from the ABIDE rs-fMRI multi-site data.

Controversies and progress on standardization of large-scale brain network nomenclature

Progress in scientific disciplines is accompanied by standardization of terminology. Network neuroscience, at the level of macroscale organization of the brain, is beginning to confront the challenges associated with developing a taxonomy of its fundamental explanatory constructs. The Workgroup for HArmonized Taxonomy of NETworks (WHATNET) was formed in 2020 as an Organization for Human Brain Mapping (OHBM)-endorsed best practices committee to provide recommendations on points of consensus, identify open questions, and highlight areas of ongoing debate in the service of moving the field toward standardized reporting of network neuroscience results. The committee conducted a survey to catalog current practices in large-scale brain network nomenclature. A few well-known network names (e.g., default mode network) dominated responses to the survey, and a number of illuminating points of disagreement emerged. We summarize survey results and provide initial considerations and recommendations from the workgroup. This perspective piece includes a selective review of challenges to this enterprise, including (1) network scale, resolution, and hierarchies; (2) interindividual variability of networks; (3) dynamics and nonstationarity of networks; (4) consideration of network affiliations of subcortical structures; and (5) consideration of multimodal information. We close with minimal reporting guidelines for the cognitive and network neuroscience communities to adopt.

High-amplitude network co-fluctuations linked to variation in hormone concentrations over the menstrual cycle

Many studies have shown that the human endocrine system modulates brain function, reporting associations between fluctuations in hormone concentrations and brain connectivity. However, how hormonal fluctuations impact fast changes in brain network organization over short timescales remains unknown. Here, we leverage a recently proposed framework for modeling co-fluctuations between the activity of pairs of brain regions at a framewise timescale. In previous studies we showed that time points corresponding to high-amplitude co-fluctuations disproportionately contributed to the time-averaged functional connectivity pattern and that these co-fluctuation patterns could be clustered into a low-dimensional set of recurring "states." Here, we assessed the relationship between these network states and quotidian variation in hormone concentrations. Specifically, we were interested in whether the frequency with which network states occurred was related to hormone concentration. We addressed this question using a dense-sampling dataset (N = 1 brain). In this dataset, a single individual was sampled over the course of two endocrine states: a natural menstrual cycle and while the subject underwent selective progesterone suppression via oral hormonal contraceptives. During each cycle, the subject underwent 30 daily resting-state fMRI scans and blood draws. Our analysis of the imaging data revealed two repeating network states. We found that the frequency with which state 1 occurred in scan sessions was significantly correlated with follicle-stimulating and luteinizing hormone concentrations. We also constructed representative networks for each scan session using only "event frames"-those time points when an event was determined to have occurred. We found that the weights of specific subsets of functional connections were robustly correlated with fluctuations in the concentration of not only luteinizing and follicle-stimulating hormones, but also progesterone and estradiol.

Hierarchical organization of spontaneous co-fluctuations in densely sampled individuals using fMRI

Edge time series decompose functional connectivity into its framewise contributions. Previous studies have focused on characterizing the properties of high-amplitude frames (time points when the global co-fluctuation amplitude takes on its largest value), including their cluster structure. Less is known about middle- and low-amplitude co-fluctuations (peaks in co-fluctuation time series but of lower amplitude). Here, we directly address those questions, using data from two dense-sampling studies: the MyConnectome project and Midnight Scan Club. We develop a hierarchical clustering algorithm to group peak co-fluctuations of all magnitudes into nested and multiscale clusters based on their pairwise concordance. At a coarse scale, we find evidence of three large clusters that, collectively, engage virtually all canonical brain systems. At finer scales, however, each cluster is dissolved, giving way to increasingly refined patterns of co-fluctuations involving specific sets of brain systems. We also find an increase in global co-fluctuation magnitude with hierarchical scale. Finally, we comment on the amount of data needed to estimate co-fluctuation pattern clusters and implications for brain-behavior studies. Collectively, the findings reported here fill several gaps in current knowledge concerning the heterogeneity and richness of co-fluctuation patterns as estimated with edge time series while providing some practical guidance for future studies.

Individual differences in time-varying and stationary brain connectivity during movie watching from childhood to early adulthood: Age, sex, and behavioral associations

Spatially remote brain regions exhibit dynamic functional interactions across various task conditions. While time-varying functional connectivity during movie watching shows sensitivity to movie content, stationary functional connectivity remains relatively stable across videos. These findings suggest that dynamic and stationary functional interactions may represent different aspects of brain function. However, the relationship between individual differences in time-varying and stationary connectivity and behavioral phenotypes remains elusive. To address this gap, we analyzed an open-access functional MRI dataset comprising participants aged 5-22 years, who watched two cartoon movie clips. We calculated regional brain activity, time-varying connectivity, and stationary connectivity, examining associations with age, sex, and behavioral assessments. Model comparison revealed that time-varying connectivity was more sensitive to age and sex effects compared with stationary connectivity. The preferred age models exhibited quadratic log age or quadratic age effects, indicative of inverted-U shaped developmental patterns. In addition, females showed higher consistency in regional brain activity and time-varying connectivity than males. However, in terms of behavioral predictions, only stationary connectivity demonstrated the ability to predict full-scale intelligence quotient. These findings suggest that individual differences in time-varying and stationary connectivity may capture distinct aspects of behavioral phenotypes.

Optimizing network neuroscience computation of individual differences in human spontaneous brain activity for test-retest reliability

A rapidly emerging application of network neuroscience in neuroimaging studies has provided useful tools to understand individual differences in intrinsic brain function by mapping spontaneous brain activity, namely intrinsic functional network neuroscience (ifNN). However, the variability of methodologies applied across the ifNN studies-with respect to node definition, edge construction, and graph measurements-makes it difficult to directly compare findings and also challenging for end users to select the optimal strategies for mapping individual differences in brain networks. Here, we aim to provide a benchmark for best ifNN practices by systematically comparing the measurement reliability of individual differences under different ifNN analytical strategies using the test-retest design of the Human Connectome Project. The results uncovered four essential principles to guide ifNN studies: (1) use a whole brain parcellation to define network nodes, including subcortical and cerebellar regions; (2) construct functional networks using spontaneous brain activity in multiple slow bands; and (3) optimize topological economy of networks at individual level; and (4) characterize information flow with specific metrics of integration and segregation. We built an interactive online resource of reliability assessments for future ifNN (https://ibraindata.com/research/ifNN).

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.

Leading basic modes of spontaneous activity drive individual functional connectivity organization in the resting human brain

Spontaneous activity of the human brain provides a window to explore intrinsic principles of functional organization. However, most studies have focused on interregional functional connectivity. The principles underlying rich repertoires of instantaneous activity remain largely unknown. We apply a recently proposed eigen-microstate analysis to three resting-state functional MRI datasets to identify basic modes that represent fundamental activity patterns that coexist over time. We identify five leading basic modes that dominate activity fluctuations. Each mode exhibits a distinct functional system-dependent coactivation pattern and corresponds to specific cognitive profiles. In particular, the spatial pattern of the first leading basis mode shows the separation of activity between the default-mode and primary and attention regions. Based on theoretical modelling, we further reconstruct individual functional connectivity as the weighted superposition of coactivation patterns corresponding to these leading basic modes. Moreover, these leading basic modes capture sleep deprivation-induced changes in brain activity and interregional connectivity, primarily involving the default-mode and task-positive regions. Our findings reveal a dominant set of basic modes of spontaneous activity that reflect multiplexed interregional coordination and drive conventional functional connectivity, furthering the understanding of the functional significance of spontaneous brain activity.

Co-evolving dynamics and topology in a coupled oscillator model of resting brain function

Dynamic models of ongoing BOLD fMRI brain dynamics and models of communication strategies have been two important approaches to understanding how brain network structure constrains function. However, dynamic models have yet to widely incorporate one of the most important insights from communication models: the brain may not use all of its connections in the same way or at the same time. Here we present a variation of a phase delayed Kuramoto coupled oscillator model that dynamically limits communication between nodes on each time step. An active subgraph of the empirically derived anatomical brain network is chosen in accordance with the local dynamic state on every time step, thus coupling dynamics and network structure in a novel way. We analyze this model with respect to its fit to empirical time-averaged functional connectivity, finding that, with the addition of only one parameter, it significantly outperforms standard Kuramoto models with phase delays. We also perform analyses on the novel time series of active edges it produces, demonstrating a slowly evolving topology moving through intermittent episodes of integration and segregation. We hope to demonstrate that the exploration of novel modeling mechanisms and the investigation of dynamics of networks in addition to dynamics on networks may advance our understanding of the relationship between brain structure and function.

Enhancing precision in human neuroscience

Human neuroscience has always been pushing the boundary of what is measurable. During the last decade, concerns about statistical power and replicability - in science in general, but also specifically in human neuroscience - have fueled an extensive debate. One important insight from this discourse is the need for larger samples, which naturally increases statistical power. An alternative is to increase the precision of measurements, which is the focus of this review. This option is often overlooked, even though statistical power benefits from increasing precision as much as from increasing sample size. Nonetheless, precision has always been at the heart of good scientific practice in human neuroscience, with researchers relying on lab traditions or rules of thumb to ensure sufficient precision for their studies. In this review, we encourage a more systematic approach to precision. We start by introducing measurement precision and its importance for well-powered studies in human neuroscience. Then, determinants for precision in a range of neuroscientific methods (MRI, M/EEG, EDA, Eye-Tracking, and Endocrinology) are elaborated. We end by discussing how a more systematic evaluation of precision and the application of respective insights can lead to an increase in reproducibility in human neuroscience.

Low-frequency repetitive transcranial magnetic stimulation alters the individual functional dynamical landscape

Repetitive transcranial magnetic stimulation (rTMS) is a noninvasive approach to modulate brain activity and behavior in humans. Still, how individual resting-state brain dynamics after rTMS evolves across different functional configurations is rarely studied. Here, using resting state fMRI data from healthy subjects, we aimed to examine the effects of rTMS to individual large-scale brain dynamics. Using Topological Data Analysis based Mapper approach, we construct the precise dynamic mapping (PDM) for each participant. To reveal the relationship between PDM and canonical functional representation of the resting brain, we annotated the graph using relative activation proportion of a set of large-scale resting-state networks (RSNs) and assigned the single brain volume to corresponding RSN-dominant or a hub state (not any RSN was dominant). Our results show that (i) low-frequency rTMS could induce changed temporal evolution of brain states; (ii) rTMS didn't alter the hub-periphery configurations underlined resting-state brain dynamics; and (iii) the rTMS effects on brain dynamics differ across the left frontal and occipital lobe. In conclusion, low-frequency rTMS significantly alters the individual temporo-spatial dynamics, and our finding further suggested a potential target-dependent alteration of brain dynamics. This work provides a new perspective to comprehend the heterogeneous effect of rTMS.