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

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.

Manifold learning for fMRI time-varying functional connectivity

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

Intermediately synchronised brain states optimise trade-off between subject specificity and predictive capacity

Functional connectivity (FC) refers to the statistical dependencies between activity of distinct brain areas. To study temporal fluctuations in FC within the duration of a functional magnetic resonance imaging (fMRI) scanning session, researchers have proposed the computation of an edge time series (ETS) and their derivatives. Evidence suggests that FC is driven by a few time points of high-amplitude co-fluctuation (HACF) in the ETS, which may also contribute disproportionately to interindividual differences. However, it remains unclear to what degree different time points actually contribute to brain-behaviour associations. Here, we systematically evaluate this question by assessing the predictive utility of FC estimates at different levels of co-fluctuation using machine learning (ML) approaches. We demonstrate that time points of lower and intermediate co-fluctuation levels provide overall highest subject specificity as well as highest predictive capacity of individual-level phenotypes.

Edge-based general linear models capture high-frequency 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 one minute to the next. However, traditional dynamic FC approaches typically lack the temporal precision to capture moment-by-moment network fluctuations. Recently, researchers have 'unfurled' traditional FC matrices in 'edge cofluctuation time series' which measure time point-by-time point cofluctuations between regions. Here we apply event-based and parametric fMRI analyses to edge time series to capture high-frequency fluctuations in networks related to attention. In two independent fMRI datasets 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.

Temporal Mapper: Transition networks in simulated and real neural dynamics

Characterizing large-scale dynamic organization of the brain relies on both data-driven and mechanistic modeling, which demands a low versus high level of prior knowledge and assumptions about how constituents of the brain interact. However, the conceptual translation between the two is not straightforward. The present work aims to provide a bridge between data-driven and mechanistic modeling. We conceptualize brain dynamics as a complex landscape that is continuously modulated by internal and external changes. The modulation can induce transitions between one stable brain state (attractor) to another. Here, we provide a novel method-Temporal Mapper-built upon established tools from the field of topological data analysis to retrieve the network of attractor transitions from time series data alone. For theoretical validation, we use a biophysical network model to induce transitions in a controlled manner, which provides simulated time series equipped with a ground-truth attractor transition network. Our approach reconstructs the ground-truth transition network from simulated time series data better than existing time-varying approaches. For empirical relevance, we apply our approach to fMRI data gathered during a continuous multitask experiment. We found that occupancy of the high-degree nodes and cycles of the transition network was significantly associated with subjects' behavioral performance. Taken together, we provide an important first step toward integrating data-driven and mechanistic modeling of brain dynamics.

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 to 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.

Homotopic local-global parcellation of the human cerebral cortex from resting-state functional connectivity

Resting-state fMRI is commonly used to derive brain parcellations, which are widely used for dimensionality reduction and interpreting human neuroscience studies. We previously developed a model that integrates local and global approaches for estimating areal-level cortical parcellations. The resulting local-global parcellations are often referred to as the Schaefer parcellations. However, the lack of homotopic correspondence between left and right Schaefer parcels has limited their use for brain lateralization studies. Here, we extend our previous model to derive homotopic areal-level parcellations. Using resting-fMRI and task-fMRI across diverse scanners, acquisition protocols, preprocessing and demographics, we show that the resulting homotopic parcellations are as homogeneous as the Schaefer parcellations, while being more homogeneous than five publicly available parcellations. Furthermore, weaker correlations between homotopic parcels are associated with greater lateralization in resting network organization, as well as lateralization in language and motor task activation. Finally, the homotopic parcellations agree with the boundaries of a number of cortical areas estimated from histology and visuotopic fMRI, while capturing sub-areal (e.g., somatotopic and visuotopic) features. Overall, these results suggest that the homotopic local-global parcellations represent neurobiologically meaningful subdivisions of the human cerebral cortex and will be a useful resource for future studies. Multi-resolution parcellations estimated from 1479 participants are publicly available (https://github.com/ThomasYeoLab/CBIG/tree/master/stable_projects/brain_parcellation/Yan2023_homotopic).

Hierarchical overlapping modular structure in the human cerebral cortex improves individual identification

The idea that brain networks have a hierarchical modular organization is pervasive. Increasing evidence suggests that brain modules overlap. However, little is known about the hierarchical overlapping modular structure in the brain. In this study, we developed a framework to uncover brain hierarchical overlapping modular structures based on a nested-spectral partition algorithm and an edge-centric network model. Overlap degree between brain modules is symmetrical across hemispheres, with highest overlap observed in the control and salience/ventral attention networks. Furthermore, brain edges are clustered into two groups: intrasystem and intersystem edges, to form hierarchical overlapping modules. At different levels, modules are self-similar in the degree of overlap. Additionally, the brain's hierarchical structure contains more individual identifiable information than a single-level structure, particularly in the control and salience/ventral attention networks. Our results offer pathways for future studies aimed at relating the organization of hierarchical overlapping modules to brain cognitive behavior and disorders.

Edge-centric functional network analyses reveal disrupted network configuration in autism spectrum disorder

BACKGROUND

Neuroscientific evidence suggests that the pathological symptoms associated with autism spectrum disorders (ASD) are not confined to a single brain region but involve networks of the brain on a larger spatial scale. Analyzing diagrams of edge-edge interactions could provide important perspectives on the organization and function of complex systems.

METHODS

Resting-state fMRI data from 238 ASD patients and 311 healthy controls (HCs) were included in the current study. We used the thalamus as the mediating node to calculate the edge functional connectivity (eFC) of the brain network and compared the ASD subjects and HCs.

RESULTS

Compared with the HCs, the ASD subjects exhibited abnormalities in the central node thalamus and four brain regions (amygdala, nucleus accumbens, pallidum and hippocampus), as well as in the eFC formed by the inferior frontal gyrus (IFG) (or middle temporal gyrus (MTG)). In addition, ASD subjects showed variable characteristics of the eFC between nodes in different networks.

CONCLUSIONS

The changes in these brain regions may be due to the disturbance in the reward system, which leads to coherence in the instantaneous comovement of the functional connections formed by these brain regions in ASD. This notion also reveals a functional network feature between the cortical and subcortical regions in ASD.

Edge-centric functional network predicts risk propensity in economic decision-making: evidence from a resting-state fMRI study

Despite node-centric studies revealing an association between resting-state functional connectivity and individual risk propensity, the prediction of future risk decisions remains undetermined. Herein, we applied a recently emerging edge-centric method, the edge community similarity network (ECSN), to alternatively describe the community structure of resting-state brain activity and to probe its contribution to predicting risk propensity during gambling. Results demonstrated that inter-individual variability of risk decisions correlates with the inter-subnetwork couplings spanning the visual network (VN) and default mode network (DMN), cingulo-opercular task control network, and sensory/somatomotor hand network (SSHN). Particularly, participants who have higher community similarity of these subnetworks during the resting state tend to choose riskier and higher yielding bets. And in contrast to low-risk propensity participants, those who behave high-risky show stronger couplings spanning the VN and SSHN/DMN. Eventually, based on the resting-state ECSN properties, the risk rate during the gambling task is effectively predicted by the multivariable linear regression model at the individual level. These findings provide new insights into the neural substrates of the inter-individual variability in risk propensity and new neuroimaging metrics to predict individual risk decisions in advance.

Multivariate information theory uncovers synergistic subsystems of the human cerebral cortex

One of the most well-established tools for modeling the brain is the functional connectivity network, which is constructed from pairs of interacting brain regions. While powerful, the network model is limited by the restriction that only pairwise dependencies are considered and potentially higher-order structures are missed. Here, we explore how multivariate information theory reveals higher-order dependencies in the human brain. We begin with a mathematical analysis of the O-information, showing analytically and numerically how it is related to previously established information theoretic measures of complexity. We then apply the O-information to brain data, showing that synergistic subsystems are widespread in the human brain. Highly synergistic subsystems typically sit between canonical functional networks, and may serve an integrative role. We then use simulated annealing to find maximally synergistic subsystems, finding that such systems typically comprise ≈10 brain regions, recruited from multiple canonical brain systems. Though ubiquitous, highly synergistic subsystems are invisible when considering pairwise functional connectivity, suggesting that higher-order dependencies form a kind of shadow structure that has been unrecognized by established network-based analyses. We assert that higher-order interactions in the brain represent an under-explored space that, accessible with tools of multivariate information theory, may offer novel scientific insights.

Brain-wide network analysis of resting-state neuromagnetic data

The present study performed a brain-wide network analysis of resting-state magnetoencephalograms recorded from 53 healthy participants to visualize elaborate brain maps of phase- and amplitude-derived graph-theory metrics at different frequencies. To achieve this, we conducted a vertex-wise computation of threshold-independent graph metrics by combining proportional thresholding and a conjunction analysis and applied them to a correlation analysis of age and brain networks. Source power showed a frequency-dependent cortical distribution. Threshold-independent graph metrics derived from phase- and amplitude-based connectivity showed similar or different distributions depending on frequency. Vertex-wise age-brain correlation maps revealed that source power at the beta band and the amplitude-based degree at the alpha band changed with age in local regions. The present results indicate that a brain-wide analysis of neuromagnetic data has the potential to reveal neurophysiological network features in the human brain in a resting state.

Connectional Hierarchy in Human Brain Revealed by Individual Variability of Functional Network Edges

The human cerebral cortex is connected by intricate inter-areal wiring at the macroscale. The cortical hierarchy from primary sensorimotor to higher-order association areas is a unifying organizational principle across various neurobiological properties; however, previous studies have not clarified whether the connections between cortical regions exhibit a similar hierarchical pattern. Here, we identify a connectional hierarchy indexed by inter-individual variability of functional connectivity edges, which continuously progresses along a hierarchical gradient from within-network connections to between-network edges connecting sensorimotor and association networks. We found that this connectional hierarchy of variability aligns with both hemodynamic and electromagnetic connectivity strength and is constrained by structural connectivity strength. Moreover, the patterning of connectional hierarchy is related to inter-regional similarity in transcriptional and neurotransmitter receptor profiles. Using the Neurosynth cognitive atlas and cortical vulnerability maps in 13 brain disorders, we found that the connectional hierarchy of variability is associated with similarity networks of cognitive relevance and that of disorder vulnerability. Finally, we found that the prominence of this hierarchical gradient of connectivity variability declines during youth. Together, our results reveal a novel hierarchal organizational principle at the connectional level that links multimodal and multiscale human connectomes to individual variability in functional connectivity.

The Entangled Brain

The Entangled Brain (Pessoa, L., 2002. MIT Press) promotes the idea that we need to understand the brain as a complex, entangled system. Why does the complex systems perspective, one that entails emergent properties, matter for brain science? In fact, many neuroscientists consider these ideas a distraction. We discuss three principles of brain organization that inform the question of the interactional complexity of the brain: (1) massive combinatorial anatomical connectivity; (2) highly distributed functional coordination; and (3) networks/circuits as functional units. To motivate the challenges of mapping structure and function, we discuss neural circuits illustrating the high anatomical and functional interactional complexity typical in the brain. We discuss potential avenues for testing for network-level properties, including those relying on distributed computations across multiple regions. We discuss implications for brain science, including the need to characterize decentralized and heterarchical anatomical-functional organization. The view advocated has important implications for causation, too, because traditional accounts of causality provide poor candidates for explanation in interactionally complex systems like the brain given the distributed, mutual, and reciprocal nature of the interactions. Ultimately, to make progress understanding how the brain supports complex mental functions, we need to dissolve boundaries within the brain-those suggested to be associated with perception, cognition, action, emotion, motivation-as well as outside the brain, as we bring down the walls between biology, psychology, mathematics, computer science, philosophy, and so on.

Brain network integration, segregation and quasi-periodic activation and deactivation during tasks and rest

Previous studies have shown that a re-organization of the brain's functional connectome expressed in terms of network integration and segregation may play a pivotal role for brain function. However, it has been proven difficult to fully capture both processes independently in a single methodological framework. In this study, by starting from pair-wise assessments of instantaneous phase synchronization and community membership, we assemble spatiotemporally flexible networks that reflect changes in integration/segregation that occur at a spectrum of spatial as well as temporal scales. This is achieved by iteratively assembling smaller networks into larger units under the constraint that the smaller units should be internally integrated, i.e. belong to the same community. The assembled subnetworks can be partly overlapping and differ in size across time. Our results show that subnetwork integration and segregation occur simultaneously in the brain. During task performance, global changes in synchronization between networks arise that are tied to the underlying temporal design of the experiment. We show that a hallmark property of the dynamics of the brain's functional connectome is a presence of quasi-periodic patterns of network activation and deactivation, which during task performance becomes intertwined with the underlying temporal structure of the experimental paradigm. Additionally, we show that the degree of network integration throughout a n-back working memory task is correlated to performance.

Uncovering individual differences in fine-scale dynamics of functional connectivity

Functional connectivity (FC) profiles contain subject-specific features that are conserved across time and have potential to capture brain-behavior relationships. Most prior work has focused on spatial features (nodes and systems) of these FC fingerprints, computed over entire imaging sessions. We propose a method for temporally filtering FC, which allows selecting specific moments in time while also maintaining the spatial pattern of node-based activity. To this end, we leverage a recently proposed decomposition of FC into edge time series (eTS). We systematically analyze functional magnetic resonance imaging frames to define features that enhance identifiability across multiple fingerprinting metrics, similarity metrics, and data sets. Results show that these metrics characteristically vary with eTS cofluctuation amplitude, similarity of frames within a run, transition velocity, and expression of functional systems. We further show that data-driven optimization of features that maximize fingerprinting metrics isolates multiple spatial patterns of system expression at specific moments in time. Selecting just 10% of the data can yield stronger fingerprints than are obtained from the full data set. Our findings support the idea that FC fingerprints are differentially expressed across time and suggest that multiple distinct fingerprints can be identified when spatial and temporal characteristics are considered simultaneously.

A spatiotemporal complexity architecture of human brain activity

The human brain operates in large-scale functional networks. These networks are an expression of temporally correlated activity across brain regions, but how global network properties relate to the neural dynamics of individual regions remains incompletely understood. Here, we show that the brain's network architecture is tightly linked to critical episodes of neural regularity, visible as spontaneous "complexity drops" in functional magnetic resonance imaging signals. These episodes closely explain functional connectivity strength between regions, subserve the propagation of neural activity patterns, and reflect interindividual differences in age and behavior. Furthermore, complexity drops define neural activity states that dynamically shape the connectivity strength, topological configuration, and hierarchy of brain networks and comprehensively explain known structure-function relationships within the brain. These findings delineate a principled complexity architecture of neural activity-a human "complexome" that underpins the brain's functional network organization.

Considering dynamic nature of the brain: the clinical importance of connectivity variability in machine learning classification and prediction

The brain connectivity of resting-state fMRI (rs-fMRI) represents an intrinsic state of brain architecture, and it has been used as a useful neural marker for detecting psychiatric conditions as well as for predicting psychosocial characteristics. However, most studies using brain connectivity have focused more on the strength of functional connectivity over time (static-FC) but less attention to temporal characteristics of connectivity changes (FC-variability). The primary goal of the current study was to investigate the effectiveness of using the FC-variability in classifying an individual's pathological characteristics from others and predicting psychosocial characteristics. In addition, the current study aimed to prove that benefits of the FC-variability are reliable across various analysis procedures. To this end, three open public large resting-state fMRI datasets including individuals with Autism Spectrum Disorder (ABIDE; N = 1249), Schizophrenia disorder (COBRE; N = 145), and typical development (NKI; N = 672) were utilized for the machine learning (ML) classification and prediction based on their static-FC and the FC-variability metrics. To confirm the robustness of FC-variability utility, we benchmarked the ML classification and prediction with various brain parcellations and sliding window parameters. As a result, we found that the ML performances were significantly improved when the ML included FC-variability features in classifying pathological populations from controls (e.g., individuals with autism spectrum disorder vs. typical development) and predicting psychiatric severity (e.g., score of autism diagnostic observation schedule), regardless of parcellation selection and sliding window size. Additionally, the ML performance deterioration was significantly prevented with FC-variability features when excessive features were inputted into the ML models, yielding more reliable results. In conclusion, the current finding proved the usefulness of the FC-variability and its reliability.

Angular gyrus: an anatomical case study for association cortex

The angular gyrus is associated with a spectrum of higher order cognitive functions. This mini-review undertakes a broad survey of putative neuroanatomical substrates, guided by the premise that area-specific specializations derive from a combination of extrinsic connections and intrinsic area properties. Three levels of spatial resolution are discussed: cellular, supracellular connectivity, and synaptic micro-scale, with examples necessarily drawn mainly from experimental work with nonhuman primates. A significant factor in the functional specialization of the human parietal cortex is the pronounced enlargement. In addition to "more" cells, synapses, and connections, however, the heterogeneity itself can be considered an important property. Multiple anatomical features support the idea of overlapping and temporally dynamic membership in several brain wide subnetworks, but how these features operate in the context of higher cognitive functions remains for continued investigations.

Ventral striatal subregional dysfunction in late-life grief: Relationships with yearning and depressive symptoms

Bereaved older adults experiencing high grief in the first year after an attachment loss is at increased risk for prolonged grief disorder (PGD) via unknown mechanisms. Yearning, a core grief symptom, is linked to the ventral striatal (VS) brain function, but the role of this neuronal system in late-life grief is poorly understood. As a first step, we examined the VS subregional abnormalities associated with multidimensional symptoms in bereaved elders during the first year post-loss. Sixty-five bereaved elders completed clinical assessments within 13 months post-loss. Ventral caudate (VCau) and nucleus accumbens (NAcc) functional connectivity (FC) was assessed using seed-based resting-state functional MRI. VCau and NAcc FC differences between high (inventory of complicated grief [ICG] score≥30; n = 35) and low (ICG score<30; n = 30) grief, and the relationships between ventral striatal subregional FC and clinical measures (yearning and depressive symptoms) were assessed after covariate adjustments (α < 0.05; 3dClustSim corrected). Relative to low grief participants, those with high grief showed higher FC between VCau and the medial prefrontal, orbitofrontal, and subgenual cingulate cortices. VCau FC abnormalities positively correlated with yearning (r2 = 0.24, p < 0.001). In contrast, FC between VCau and temporoparietal junction negatively correlated with depressive symptoms, a commonly co-occurring symptom (r2 = 0.37, p < 0.001). The FC between NAcc and insula/striatum positively correlated with yearning (r2 = 0.35, p < 0.001); no other NAcc FC findings were seen in the full sample. In women, higher FC between the NAcc and bilateral posterior cingulate, precuneus, and visual areas were found in those with high, relative to low grief symptoms. Distinct VS subregional abnormalities associate with yearning and depressive symptoms in bereaved elders. Whether ventral striatal dysfunction correlates with PGD development and/or worsening depression remains to be elucidated.

Leveraging edge-centric networks complements existing network-level inference for functional connectomes

The human connectome is modular with distinct brain regions clustering together to form large-scale communities, or networks. This concept has recently been leveraged in novel inferencing procedures by averaging the edge-level statistics within networks to induce more powerful inferencing at the network level. However, these networks are constructed based on the similarity between pairs of nodes. Emerging work has described novel edge-centric networks, which instead use the similarity between pairs of edges to construct networks. In this work, we use these edge-centric networks in a network-level inferencing procedure and compare this novel method to traditional inferential procedures and the network-level procedure using node-centric networks. We use data from the Human Connectome Project, the Healthy Brain Network, and the Philadelphia Neurodevelopmental Cohort and use a resampling technique with various sample sizes (n=40, 80, 120) to probe the power and specificity of each method. Across datasets and sample sizes, using the edge-centric networks outperforms using node-centric networks for inference as well as edge-level FDR correction and NBS. Additionally, the edge-centric networks were found to be more consistent in clustering effect sizes of similar values as compared to node-centric networks, although node-centric networks often had a lower average within-network effect size variability. Together, these findings suggest that using edge-centric networks for network-level inference can procure relatively powerful results while remaining similarly accurate to the underlying edge-level effects across networks, complementing previous inferential methods.

Evidence of familial confounding of the association between cannabis use and cerebellar-cortical functional connectivity using a twin study

Cerebellar-cortical resting-state functional connectivity (rsFC) has been reported to be altered in cannabis users. However, this association may be due to genetic and environmental confounding rather than a causal relationship between cannabis use and changes in rsFC. In this co-twin control study, linear mixed models were used to assess relationships between the number of lifetime cannabis uses (NLCU) and age of cannabis onset (ACO) with cerebellar-cortical rsFC. The rsFC with seven functional networks was evaluated in 147 monozygotic and 82 dizygotic twin pairs. Importantly, the use of genetically informed models in this twin sample facilitated examining whether shared genetic or environmental effects underlie crude associations between cannabis measures and connectivity. Individual-level phenotypic analyses (i.e., accounting for twin-pair non-independence) showed that individuals in the full sample with earlier ACO and higher NLCU had lower cerebellar rsFC within the VA, DA, and FP networks. Yet, there were no significant differences in cerebellar-cortical rsFC between monozygotic twins who were discordant for cannabis measures. These findings suggest shared genetic or environmental confounds contribute to associations between cannabis use and altered cerebellar-cortical rsFC, rather than unique causal impacts of cannabis use on cerebellar-cortical rsFC.

It's about time: Linking dynamical systems with human neuroimaging to understand the brain

Most human neuroscience research to date has focused on statistical approaches that describe stationary patterns of localized neural activity or blood flow. While these patterns are often interpreted in light of dynamic, information-processing concepts, the static, local, and inferential nature of the statistical approach makes it challenging to directly link neuroimaging results to plausible underlying neural mechanisms. Here, we argue that dynamical systems theory provides the crucial mechanistic framework for characterizing both the brain's time-varying quality and its partial stability in the face of perturbations, and hence, that this perspective can have a profound impact on the interpretation of human neuroimaging results and their relationship with behavior. After briefly reviewing some key terminology, we identify three key ways in which neuroimaging analyses can embrace a dynamical systems perspective: by shifting from a local to a more global perspective, by focusing on dynamics instead of static snapshots of neural activity, and by embracing modeling approaches that map neural dynamics using "forward" models. Through this approach, we envisage ample opportunities for neuroimaging researchers to enrich their understanding of the dynamic neural mechanisms that support a wide array of brain functions, both in health and in the setting of psychopathology.

The Brain Connectome for Chinese Reading

Chinese, as a logographic language, fundamentally differs from alphabetic languages like English. Previous neuroimaging studies have mainly focused on alphabetic languages, while the exploration of Chinese reading is still an emerging and fast-growing research field. Recently, a growing number of neuroimaging studies have explored the neural circuit of Chinese reading. Here, we summarize previous research on Chinese reading from a connectomic perspective. Converging evidence indicates that the left middle frontal gyrus is a specialized hub region that connects the ventral with dorsal pathways for Chinese reading. Notably, the orthography-to-phonology and orthography-to-semantics mapping, mainly processed in the ventral pathway, are more specific during Chinese reading. Besides, in addition to the left-lateralized language-related regions, reading pathways in the right hemisphere also play an important role in Chinese reading. Throughout, we comprehensively review prior findings and emphasize several challenging issues to be explored in future work.

Dynamic and stationary brain connectivity during movie watching as revealed by functional MRI

Spatially remote brain regions show synchronized activity as typically revealed by correlated functional MRI (fMRI) signals. An emerging line of research has focused on the temporal fluctuations of connectivity; however, its relationships with stationary connectivity have not been clearly illustrated. We examined dynamic and stationary connectivity when the participants watched four different movie clips. We calculated point-by-point multiplication between two regional time series to estimate the time-resolved dynamic connectivity, and estimated the inter-individual consistency of the dynamic connectivity time series. Widespread consistent dynamic connectivity was observed for each movie clip, which also showed differences between the clips. For example, a cartoon movie clip, Wall-E, showed more consistent of dynamic connectivity with the posterior cingulate cortex and supramarginal gyrus, while a court drama clip, A Few Good Men, showed more consistent of dynamic connectivity with the auditory cortex and temporoparietal junction, which might suggest the involvement of specific brain processing for different movie contents. In contrast, the stationary connectivity as measured by the correlations between regional time series was highly similar among the movie clips, and showed fewer statistically significant differences. The patterns of consistent dynamic connectivity could be used to classify different movie clips with higher accuracy than the stationary connectivity and regional activity. These results support the functional significance of dynamic connectivity in reflecting functional brain changes, which could provide more functionally relevant information than stationary connectivity.

Subcortical-cortical dynamical states of the human brain and their breakdown in stroke

The mechanisms controlling dynamical patterns in spontaneous brain activity are poorly understood. Here, we provide evidence that cortical dynamics in the ultra-slow frequency range (<0.01–0.1 Hz) requires intact cortical-subcortical communication. Using functional magnetic resonance imaging (fMRI) at rest, we identify Dynamic Functional States (DFSs), transient but recurrent clusters of cortical and subcortical regions synchronizing at ultra-slow frequencies. We observe that shifts in cortical clusters are temporally coincident with shifts in subcortical clusters, with cortical regions flexibly synchronizing with either limbic regions (hippocampus/amygdala), or subcortical nuclei (thalamus/basal ganglia). Focal lesions induced by stroke, especially those damaging white matter connections between basal ganglia/thalamus and cortex, provoke anomalies in the fraction times, dwell times, and transitions between DFSs, causing a bias toward abnormal network integration. Dynamical anomalies observed 2 weeks after stroke recover in time and contribute to explaining neurological impairment and long-term outcome.

Favaretto et al. show that the brain rapidly alternates between transient connectivity patterns, with cortical regions flexibly synchronizing with two groups of subcortical regions, and that this dynamic is abnormal in stroke patients.

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

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

Precision dynamical mapping using topological data analysis reveals a hub-like transition state at rest

In the absence of external stimuli, neural activity continuously evolves from one configuration to another. Whether these transitions or explorations follow some underlying arrangement or lack a predictable ordered plan remains to be determined. Here, using fMRI data from highly sampled individuals (~5 hours of resting-state data per individual), we aimed to reveal the rules that govern transitions in brain activity at rest. Our Topological Data Analysis based Mapper approach characterized a highly visited transition state of the brain that acts as a switch between different neural configurations to organize the spontaneous brain activity. Further, while the transition state was characterized by a uniform representation of canonical resting-state networks (RSNs), the periphery of the landscape was dominated by a subject-specific combination of RSNs. Altogether, we revealed rules or principles that organize spontaneous brain activity using a precision dynamics approach.

Improving power in functional magnetic resonance imaging by moving beyond cluster-level inference

Inference in neuroimaging typically occurs at the level of focal brain areas or circuits. Yet, increasingly, well-powered studies paint a much richer picture of broad-scale effects distributed throughout the brain, suggesting that many focal reports may only reflect the tip of the iceberg of underlying effects. How focal versus broad-scale perspectives influence the inferences we make has not yet been comprehensively evaluated using real data. Here, we compare sensitivity and specificity across procedures representing multiple levels of inference using an empirical benchmarking procedure that resamples task-based connectomes from the Human Connectome Project dataset (∼1,000 subjects, 7 tasks, 3 resampling group sizes, 7 inferential procedures). Only broad-scale (network and whole brain) procedures obtained the traditional 80% statistical power level to detect an average effect, reflecting >20% more statistical power than focal (edge and cluster) procedures. Power also increased substantially for false discovery rate- compared with familywise error rate-controlling procedures. The downsides are fairly limited; the loss in specificity for broad-scale and FDR procedures was relatively modest compared to the gains in power. Furthermore, the broad-scale methods we introduce are simple, fast, and easy to use, providing a straightforward starting point for researchers. This also points to the promise of more sophisticated broad-scale methods for not only functional connectivity but also related fields, including task-based activation. Altogether, this work demonstrates that shifting the scale of inference and choosing FDR control are both immediately attainable and can help remedy the issues with statistical power plaguing typical studies in the field.

Altered brain network topology in children with auditory processing disorder: A resting-state multi-echo fMRI study

Children with auditory processing disorder (APD) experience hearing difficulties, particularly in the presence of competing sounds, despite having normal audiograms. There is considerable debate on whether APD symptoms originate from bottom-up (e.g., auditory sensory processing) and/or top-down processing (e.g., cognitive, language, memory). A related issue is that little is known about whether functional brain network topology is altered in APD. Therefore, we used resting-state functional magnetic resonance imaging data to investigate the functional brain network organization of 57 children from 8 to 14 years old, diagnosed with APD (n = 28) and without hearing difficulties (healthy control, HC; n = 29). We applied complex network analysis using graph theory to assess the whole-brain integration and segregation of functional networks and brain hub architecture. Our results showed children with APD and HC have similar global network properties -i.e., an average of all brain regions- and modular organization. Still, the APD group showed different hub architecture in default mode-ventral attention, somatomotor and frontoparietal-dorsal attention modules. At the nodal level -i.e., single-brain regions-, we observed decreased participation coefficient (PC - a measure quantifying the diversity of between-network connectivity) in auditory cortical regions in APD, including bilateral superior temporal gyrus and left middle temporal gyrus. Beyond auditory regions, PC was also decreased in APD in bilateral posterior temporo-occipital cortices, left intraparietal sulcus, and right posterior insular cortex. Correlation analysis suggested a positive association between PC in the left parahippocampal gyrus and the listening-in-spatialized-noise -sentences task where APD children were engaged in auditory perception. In conclusion, our findings provide evidence of altered brain network organization in children with APD, specific to auditory networks, and shed new light on the neural systems underlying children's listening difficulties.

Closing the mechanistic gap: the value of microarchitecture in understanding cognitive networks

Cognitive neuroscience aims to provide biologically relevant accounts of cognition. Contemporary research linking spatial patterns of neural activity to psychological constructs describes 'where' hypothesised functions occur, but not 'how' these regions contribute to cognition. Technological, empirical, and conceptual advances allow this mechanistic gap to be closed by embedding patterns of functional activity in macro- and microscale descriptions of brain organisation. Recent work on the default mode network (DMN) and the multiple demand network (MDN), for example, highlights a microarchitectural landscape that may explain how activity in these networks integrates varied information, thus providing an anatomical foundation that will help to explain how these networks contribute to many different cognitive states. This perspective highlights emerging insights into how microarchitecture can constrain network accounts of human cognition.

Structural Brain Network Abnormalities in Parkinson’s Disease With Freezing of Gait

Objective

Diffusion tensor imaging (DTI) studies have investigated white matter (WM) integrity abnormalities in Parkinson’s disease (PD). However, little is known about the topological changes in the brain network. This study aims to reveal these changes by comparing PD without freezing of gait (FOG) (PD FOG–), PD with FOG (PD FOG+), and healthy control (HC).

Methods

21 PD FOG+, 34 PD FOG-, and 23 HC were recruited, and DTI images were acquired. The graph theoretical analysis and network-based statistical method were used to calculate the topological parameters and assess connections.

Results

PD FOG+ showed a decreased normalized clustering coefficient, small-worldness, clustering coefficient, and increased local network efficiency compared with HCs. PD FOG+ showed decreased centrality, degree centrality, and nodal efficiency in the striatum, frontal gyrus, and supplementary motor area (SMA). PD FOG+ showed decreased connections in the frontal gyrus, cingulate gyrus, and caudate nucleus (CAU). The between centrality of the left SMA and left CAU was negatively correlated with FOG questionnaire scores.

Conclusion

This study demonstrates that PD FOG+ exhibits disruption of global and local topological organization in structural brain networks, and the disrupted topological organization can be potential biomarkers in PD FOG+. These new findings may provide increasing insight into the pathophysiological mechanism of PD FOG+.

Brain network analysis between Parkinson's Disease and Health Control based on edge functional connectivity

Parkinson's Disease (PD) is the second largest neurodegenerative disease. Brain functional connectivity (FC) studies for PD were useful. In this study, we employed a novel brain network construction method, edge functional connectivity (eFC), to explore FC differences between healthy control (HC) subjects and PD patients. The data used in this study included 34 HCs and 47 PDs from Huashan Hospital, Fudan University, China. Resting state functional magnetic resonance imaging (rsfMRI) and clinical information were selected. Firstly, we constructed eFC brain network and calculated network matrix for the HC and PD groups. Then, we compared brain network matrix between eFC and the traditional nodal functional connectivity (nFC) method. Receiver operating characteristic curve (ROC) analysis was applied to validate the efficiency of the eFC brain network. The results showed that both nFC and eFC brain networks could identify significantly different characteristics between the HC and PD groups. Important hubs were mainly concentrated in visual network, sensorimotor network, subcortex and cerebellum. In addition, new hubs in basal ganglia and cerebellum regions were found in eFC. Furthermore, eFC achieved better classification results (AUC=0.985) than nFC (AUC=0.861) in discriminating PD from CN subjects.

Edge-based network analysis reveals frequency-specific network dynamics in aberrant anxiogenic processing in rats

Uncovering interactions between edges of brain networks can reveal the organizational principle of the networks and also their dysregulations underlying aberrant behaviours such as in neuropsychiatric diseases. In this study, we looked into the applicability of edge-based network analysis in uncovering possible network mechanisms of aberrant anxiogenic processing. Utilizing a rat model of prodromal Parkinson's disease we examined how a dorsomedial striatum-tied associative network (DSAN) may mediate context-based anxiogenic behaviour. Following dopamine depletion in the dorsomedial striatum, an exaggerated bottom-up signalling (posterior parietal-hippocampal-retrosplenial to anterior prefrontal-cingulate-amygdala regions) and gradient specific to the theta frequency in this network was observed. This change was accompanied by increased anxiety behaviour of the animals. By employing an edge-based approach in correlating informational flow (phase transfer entropy) with functional connectivity of all edges of this network, we further explore how the abnormal bottom-up signalling might be explained by alterations to the informational flow-connectivity motifs in the network. Our results demonstrate usage of edge-based network analysis in revealing concurrent informational processing and functional organization dynamics across multiple pathways in a brain network. This approach in unveiling network abnormalities and its impact on behavioural outcomes would be useful in probing the network basis of neuropsychiatric conditions.

Temporally and anatomically specific contributions of the human amygdala to threat and safety learning

Neural plasticity in subareas of the rodent amygdala is widely known to be essential for Pavlovian threat conditioning and safety learning. However, less consistent results have been observed in human neuroimaging studies. Here, we identify and test three important factors that may contribute to these discrepancies: the temporal profile of amygdala response in threat conditioning, the anatomical specificity of amygdala responses during threat conditioning and safety learning, and insufficient power to identify these responses. We combined data across multiple studies using a well-validated human threat conditioning paradigm to examine amygdala involvement during threat conditioning and safety learning. In 601 humans, we show that two amygdala subregions tracked the conditioned stimulus with aversive shock during early conditioning while only one demonstrated delayed responding to a stimulus not paired with shock. Our findings identify cross-species similarities in temporal- and anatomical-specific amygdala contributions to threat and safety learning, affirm human amygdala involvement in associative learning and highlight important factors for future associative learning research in humans.

Unbalanced functional connectivity at rest affects the ERP correlates of affective prediction in high intolerance of uncertainty individuals: A high density EEG investigation

In a recent study we outlined the link between Intolerance of Uncertainty (IU) and the neural correlates of affective predictions, as constructed by the brain (generation stage) to prepare to relevant stimuli (implementation stage), and update predictive models according to incoming stimuli (updating stage). In this study we further explored whether the brain's functional organization at rest can modulate neural activity elicited within an emotional S1-S2 paradigm as a function of IU and uncertainty of S1-S2 contingencies. We computed resting state functional connectivity (RS-FC) from a 3-min resting period recorded with high density EEG, and we tested whether RS graph theory nodal measures (i.e., strength, clustering coefficient, betweenness centrality) predicted in-task ERP modulation as a function of IU. We found that RS-FC differently predicted in-task ERPs within the generation and updating stages. Higher IU levels were associated to altered RS-FC patterns within both domain-specific (i.e., right superior temporal sulcus) and domain-general regions (i.e., right orbitofrontal cortex), predictive of a reduced modulation of in-task ERPs in the generation and updating stages. This is presumably ascribable to an unbalancing between synchronization and integration within these regions, which may disrupt the exchange of information between top-down and bottom-up pathways. This altered RS-FC pattern may in turn result in the construction of less efficient affective predictions and a reduced ability to deal with contextual uncertainty in individuals high in IU.

Interacting brains revisited: A cross-brain network neuroscience perspective

Elucidating the neural basis of social behavior is a long‐standing challenge in neuroscience. Such endeavors are driven by attempts to extend the isolated perspective on the human brain by considering interacting persons' brain activities, but a theoretical and computational framework for this purpose is still in its infancy. Here, we posit a comprehensive framework based on bipartite graphs for interbrain networks and address whether they provide meaningful insights into the neural underpinnings of social interactions. First, we show that the nodal density of such graphs exhibits nonrandom properties. While the current hyperscanning analyses mostly rely on global metrics, we encode the regions' roles via matrix decomposition to obtain an interpretable network representation yielding both global and local insights. With Bayesian modeling, we reveal how synchrony patterns seeded in specific brain regions contribute to global effects. Beyond inferential inquiries, we demonstrate that graph representations can be used to predict individual social characteristics, outperforming functional connectivity estimators for this purpose. In the future, this may provide a means of characterizing individual variations in social behavior or identifying biomarkers for social interaction and disorders.

Neural connectome features of procrastination: Current progress and future direction

Procrastination refers to an irrationally delay for intended courses of action despite of anticipating a negative consequence due to this delay. Previous studies tried to reveal the neural substrates of procrastination in terms of connectome-based biomarkers. Based on this, we proposed a unified triple brain network model for procrastination and pinpointed out what challenges we are facing in understanding neural mechanism of procrastination. Specifically, based on neuroanatomical features, the unified triple brain network model proposed that connectome-based underpinning of procrastination could be ascribed to the abnormalities of self-control network (i.e., dorsolateral prefrontal cortex, DLPFC), emotion-regulation network (i.e., orbital frontal cortex, OFC), and episodic prospection network (i.e., para-hippocampus cortex, PHC). Moreover, based on the brain functional features, procrastination had been attributed to disruptive neural circuits on FPN (frontoparietal network)-SCN (subcortical network) and FPN-SAN (salience network), which led us to hypothesize the crucial roles of interplay between these networks on procrastination in unified triple brain network model. Despite of these findings, poor interpretability and computational model limited further understanding for procrastination from theoretical and neural perspectives. On balance, the current study provided an overview to show current progress on the connectome-based biomarkers for procrastination, and proposed the integrative neurocognitive model of procrastination.

Time-resolved structure-function coupling in brain networks

The relationship between structural and functional connectivity in the brain is a key question in systems neuroscience. Modern accounts assume a single global structure-function relationship that persists over time. Here we study structure-function coupling from a dynamic perspective, and show that it is regionally heterogeneous. We use a temporal unwrapping procedure to identify moment-to-moment co-fluctuations in neural activity, and reconstruct time-resolved structure-function coupling patterns. We find that patterns of dynamic structure-function coupling are region-specific. We observe stable coupling in unimodal and transmodal cortex, and dynamic coupling in intermediate regions, particularly in insular cortex (salience network) and frontal eye fields (dorsal attention network). Finally, we show that the variability of a region’s structure-function coupling is related to the distribution of its connection lengths. Collectively, our findings provide a way to study structure-function relationships from a dynamic perspective.

Temporal unwrapping analysis of diffusion weighted MRI connectivity and functional MRI scans reveals that the coupling between structure and function in the human brain is regionally heterogeneous and provides a framework to evaluate these relationships from a dynamic perspective.

Edge-centric analysis of stroke patients: An alternative approach for biomarkers of lesion recovery

Most neuroimaging studies of post-stroke recovery rely on analyses derived from standard node-centric functional connectivity to map the distributed effects in stroke patients. Here, given the importance of nonlocal and diffuse damage, we use an edge-centric approach to functional connectivity in order to provide an alternative description of the effects of this disorder. These techniques allow for the rendering of metrics such as normalized entropy, which describes the diversity of edge communities at each node. Moreover, the approach enables the identification of high amplitude co-fluctuations in fMRI time series. We found that normalized entropy is associated with stroke lesion severity and continually increases across the time of patients' recovery. Furthermore, high amplitude co-fluctuations not only relate to the lesion severity but are also associated with patients' level of recovery. The current study is the first edge-centric application for a clinical population in a longitudinal dataset and demonstrates how a different perspective for functional data analysis can further characterize topographic modulations of brain dynamics.

A mathematical perspective on edge-centric brain functional connectivity

Edge time series are increasingly used in brain imaging to study the node functional connectivity (nFC) dynamics at the finest temporal resolution while avoiding sliding windows. Here, we lay the mathematical foundations for the edge-centric analysis of neuroimaging time series, explaining why a few high-amplitude cofluctuations drive the nFC across datasets. Our exposition also constitutes a critique of the existing edge-centric studies, showing that their main findings can be derived from the nFC under a static null hypothesis that disregards temporal correlations. Testing the analytic predictions on functional MRI data from the Human Connectome Project confirms that the nFC can explain most variation in the edge FC matrix, the edge communities, the large cofluctuations, and the corresponding spatial patterns. We encourage the use of dynamic measures in future research, which exploit the temporal structure of the edge time series and cannot be replicated by static null models.

Dissociable multi-scale patterns of development in personalized brain networks

The brain is organized into networks at multiple resolutions, or scales, yet studies of functional network development typically focus on a single scale. Here, we derive personalized functional networks across 29 scales in a large sample of youths (n = 693, ages 8-23 years) to identify multi-scale patterns of network re-organization related to neurocognitive development. We found that developmental shifts in inter-network coupling reflect and strengthen a functional hierarchy of cortical organization. Furthermore, we observed that scale-dependent effects were present in lower-order, unimodal networks, but not higher-order, transmodal networks. Finally, we found that network maturation had clear behavioral relevance: the development of coupling in unimodal and transmodal networks are dissociably related to the emergence of executive function. These results suggest that the development of functional brain networks align with and refine a hierarchy linked to cognition.

Neural alterations in opioid-exposed infants revealed by edge-centric brain functional networks

Prenatal opioid exposure has been linked to adverse effects spanning multiple neurodevelopmental domains, including cognition, motor development, attention, and vision. However, the neural basis of these abnormalities is largely unknown. A total of 49 infants, including 21 opioid-exposed and 28 controls, were enrolled and underwent MRI (43 ± 6 days old) after birth, including resting state functional MRI. Edge-centric functional networks based on dynamic functional connections were constructed, and machine-learning methods were employed to identify neural features distinguishing opioid-exposed infants from unexposed controls. An accuracy of 73.6% (sensitivity 76.25% and specificity 69.33%) was achieved using 10 times 10-fold cross-validation, which substantially outperformed those obtained using conventional static functional connections (accuracy 56.9%). More importantly, we identified that prenatal opioid exposure preferentially affects inter- rather than intra-network dynamic functional connections, particularly with the visual, subcortical, and default mode networks. Consistent results at the brain regional and connection levels were also observed, where the brain regions and connections associated with visual and higher order cognitive functions played pivotal roles in distinguishing opioid-exposed infants from controls. Our findings support the clinical phenotype of infants exposed to opioids in utero and may potentially explain the higher rates of visual and emotional problems observed in this population. Finally, our findings suggested that edge-centric networks could better capture the neural differences between opioid-exposed infants and controls by abstracting the intrinsic co-fluctuation along edges, which may provide a promising tool for future studies focusing on investigating the effects of prenatal opioid exposure on neurodevelopment.

Cortico-subcortical interactions in overlapping communities of edge functional connectivity

Both cortical and subcortical regions can be functionally organized into networks. Regions of the basal ganglia are extensively interconnected with the cortex via reciprocal connections that relay and modulate cortical function. Here we employ an edge-centric approach, which computes co-fluctuations among region pairs in a network to investigate the role and interaction of subcortical regions with cortical systems. By clustering edges into communities, we show that cortical systems and subcortical regions couple via multiple edge communities, with hippocampus and amygdala having a distinct pattern from striatum and thalamus. We show that the edge community structure of cortical networks is highly similar to one obtained from cortical nodes when the subcortex is present in the network. Additionally, we show that the edge community profile of both cortical and subcortical nodes can be estimates solely from cortico-subcortical interactions. Finally, we used a motif analysis focusing on edge community triads where a subcortical region coupled to two cortical regions and found that two community triads where one community couples the subcortex to the cortex were overrepresented. In summary, our results show organized coupling of the subcortex to the cortex that may play a role in cortical organization of primary sensorimotor/attention and heteromodal systems and puts forth the motif analysis of edge community triads as a promising method for investigation of communication patterns in networks.

Statelets: Capturing recurrent transient variations in dynamic functional network connectivity

Dynamic functional network connectivity (dFNC) analysis is a widely used approach for capturing brain activation patterns, connectivity states, and network organization. However, a typical sliding window plus clustering (SWC) approach for analyzing dFNC models the system through a fixed sequence of connectivity states. SWC assumes connectivity patterns span throughout the brain, but they are relatively spatially constrained and temporally short‐lived in practice. Thus, SWC is neither designed to capture transient dynamic changes nor heterogeneity across subjects/time. We propose a state‐space time series summarization framework called “statelets” to address these shortcomings. It models functional connectivity dynamics at fine‐grained timescales, adapting time series motifs to changes in connectivity strength, and constructs a concise yet informative representation of the original data that conveys easily comprehensible information about the phenotypes. We leverage the earth mover distance in a nonstandard way to handle scale differences and utilize kernel density estimation to build a probability density profile for local motifs. We apply the framework to study dFNC of patients with schizophrenia (SZ) and healthy control (HC). Results demonstrate SZ subjects exhibit reduced modularity in their brain network organization relative to HC. Statelets in the HC group show an increased recurrence across the dFNC time‐course compared to the SZ. Analyzing the consistency of the connections across time reveals significant differences within visual, sensorimotor, and default mode regions where HC subjects show higher consistency than SZ. The introduced approach also enables handling dynamic information in cross‐modal and multimodal applications to study healthy and disordered brains.

Predicting Brain Amyloid-β PET Grades with Graph Convolutional Networks Based on Functional MRI and Multi-Level Functional Connectivity

BACKGROUND

The detection of amyloid-β (Aβ) deposition in the brain provides crucial evidence in the clinical diagnosis of Alzheimer's disease (AD). However, the current positron emission tomography (PET)-based brain Aβ examination suffers from the problems of coarse visual inspection (in many cases, with 2-class stratification) and high scanning cost.

OBJECTIVE

1) To characterize the non-binary Aβ deposition levels in the AD continuum based on clustering of PET data, and 2) to explore the feasibility of predicting individual Aβ deposition grades with non-invasive functional magnetic resonance imaging (fMRI).

METHODS

1) Individual whole-brain Aβ-PET images from the OASIS-3 dataset (N = 258) were grouped into three clusters (grades) with t-SNE and k-means. The demographical data as well as global and regional standard uptake value ratios (SUVRs) were compared among the three clusters with Chi-square tests or ANOVA tests. 2) From resting-state fMRI, both conventional functional connectivity (FC) and high-order FC networks were constructed and the topological architectures of the two networks were jointly learned with graph convolutional networks (GCNs) to predict the Aβ-PET grades for each individual.

RESULTS

We found three clearly separated clusters, indicating three Aβ-PET grades. There were significant differences in gender, age, cognitive ability, APOE type, as well as global and regional SUVRs among the three grades we found. The prediction of Aβ-PET grades with GCNs on FC for the 258 participants in the AD continuum reached a satisfactory averaged accuracy (78.8%) in the two-class classification tasks.

CONCLUSION

The results demonstrated the feasibility of using deep learning on a non-invasive brain functional imaging technique to approximate PET-based Aβ deposition grading.

Linking Multi-Layer Dynamical GCN With Style-Based Recalibration CNN for EEG-Based Emotion Recognition

Electroencephalography (EEG)-based emotion computing has become one of the research hotspots of human-computer interaction (HCI). However, it is difficult to effectively learn the interactions between brain regions in emotional states by using traditional convolutional neural networks because there is information transmission between neurons, which constitutes the brain network structure. In this paper, we proposed a novel model combining graph convolutional network and convolutional neural network, namely MDGCN-SRCNN, aiming to fully extract features of channel connectivity in different receptive fields and deep layer abstract features to distinguish different emotions. Particularly, we add style-based recalibration module to CNN to extract deep layer features, which can better select features that are highly related to emotion. We conducted two individual experiments on SEED data set and SEED-IV data set, respectively, and the experiments proved the effectiveness of MDGCN-SRCNN model. The recognition accuracy on SEED and SEED-IV is 95.08 and 85.52%, respectively. Our model has better performance than other state-of-art methods. In addition, by visualizing the distribution of different layers features, we prove that the combination of shallow layer and deep layer features can effectively improve the recognition performance. Finally, we verified the important brain regions and the connection relationships between channels for emotion generation by analyzing the connection weights between channels after model learning.

On co-activation pattern analysis and non-stationarity of resting brain activity

The non-stationarity of resting-state brain activity has received increasing attention in recent years. Functional connectivity (FC) analysis with short sliding windows and coactivation pattern (CAP) analysis are two widely used methods for assessing the dynamic characteristics of brain activity observed with functional magnetic resonance imaging (fMRI). However, the statistical nature of the dynamics captured by these techniques needs to be verified. In this study, we found that the results of CAP analysis were similar for real fMRI data and simulated stationary data with matching covariance structures and spectral contents. We also found that, for both the real and simulated data, CAPs were clustered into spatially heterogeneous modules. Moreover, for each of the modules in the real data, a spatially similar module was found in the simulated data. The present results suggest that care needs to be taken when interpreting observations drawn from CAP analysis as it does not necessarily reflect non-stationarity or a mixture of states in resting brain activity.