Timescales of spontaneous fMRI fluctuations relate to structural connectivity in the brain

A network control theory pipeline for studying the dynamics of the structural connectome

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains the dynamics of a system. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter the dynamics of a system in a desired way. An interesting development for NCT in the neuroscience field is its application to study behavior and mental health symptoms. To date, NCT has been validated to study different aspects of the human structural connectome. NCT outputs can be monitored throughout developmental stages to study the effects of connectome topology on neural dynamics and, separately, to test the coherence of empirical datasets with brain function and stimulation. Here, we provide a comprehensive pipeline for applying NCT to structural connectomes by following two procedures. The main procedure focuses on computing the control energy associated with the transitions between specific neural activity states. The second procedure focuses on computing average controllability, which indexes nodes' general capacity to control the dynamics of the system. We provide recommendations for comparing NCT outputs against null network models, and we further support this approach with a Python-based software package called 'network control theory for python'. The procedures in this protocol are appropriate for users with a background in network neuroscience and experience in dynamical systems theory.

Extracting interpretable signatures of whole-brain dynamics through systematic comparison

The brain's complex distributed dynamics are typically quantified using a limited set of manually selected statistical properties, leaving the possibility that alternative dynamical properties may outperform those reported for a given application. Here, we address this limitation by systematically comparing diverse, interpretable features of both intra-regional activity and inter-regional functional coupling from resting-state functional magnetic resonance imaging (rs-fMRI) data, demonstrating our method using case-control comparisons of four neuropsychiatric disorders. Our findings generally support the use of linear time-series analysis techniques for rs-fMRI case-control analyses, while also identifying new ways to quantify informative dynamical fMRI structures. While simple statistical representations of fMRI dynamics performed surprisingly well (e.g., properties within a single brain region), combining intra-regional properties with inter-regional coupling generally improved performance, underscoring the distributed, multifaceted changes to fMRI dynamics in neuropsychiatric disorders. The comprehensive, data-driven method introduced here enables systematic identification and interpretation of quantitative dynamical signatures of multivariate time-series data, with applicability beyond neuroimaging to diverse scientific problems involving complex time-varying systems.

Beyond network connectivity: A classification approach to brain age prediction with resting-state fMRI

The brain is a complex, dynamic organ that shows differences in the same subject at various periods. Understanding how brain activity changes across age as a function of the brain networks has been greatly abetted by fMRI. Canonical analysis consists of determining how alterations in connectivity patterns (CPs) of certain regions are affected. An alternative approach is taken here by not considering connectivity but rather features computed from recordings at the regions of interest (ROIs). Using machine learning (ML) we assess how neural signals are altered by and prospectively predictive of age and sex via a methodology that is novel in drawing upon pairwise classification across six decades of subjects' chronological ages. ML is used to answer the equally important questions of what properties of the computed features are most predictive as well as which brain networks are most affected by aging. It was found that there is decreased differentiation among the neural signals of older subjects that are separated in age by the same number of years as younger subjects. Furthermore, the burstiness of the signals change at different rates between males and females. The findings provide insight into brain aging via an ROI-based analysis, the consideration of several feature groups, and a novel classification-based ML pipeline. There is also a contribution to understanding the effects of data aggregated from different recording centers on the conclusions of fMRI studies.

The significance of structural rich club hubs for the processing of hierarchical stimuli

The brain's structural network follows a hierarchy that is described as rich club (RC) organization, with RC hubs forming the well-interconnected top of this hierarchy. In this study, we tested whether RC hubs are involved in the processing of hierarchically higher structures in stimulus sequences. Moreover, we explored the role of previously suggested cortical gradients along anterior-posterior and medial-lateral axes throughout the frontal cortex. To this end, we conducted a functional magnetic resonance imaging (fMRI) experiment and presented participants with blocks of digit sequences that were structured on different hierarchically nested levels. We additionally collected diffusion weighted imaging data of the same subjects to identify RC hubs. This classification then served as the basis for a region of interest analysis of the fMRI data. Moreover, we determined structural network centrality measures in areas that were found as activation clusters in the whole-brain fMRI analysis. Our findings support the previously found anterior and medial shift for processing hierarchically higher structures of stimuli. Additionally, we found that the processing of hierarchically higher structures of the stimulus structure engages RC hubs more than for lower levels. Areas involved in the functional processing of hierarchically higher structures were also more likely to be part of the structural RC and were furthermore more central to the structural network. In summary, our results highlight the potential role of the structural RC organization in shaping the cortical processing hierarchy.

Evidence for Multiscale Multiplexed Representation of Visual Features in EEG

Distinct neural processes such as sensory and memory processes are often encoded over distinct timescales of neural activations. Animal studies have shown that this multiscale coding strategy is also implemented for individual components of a single process, such as individual features of a multifeature stimulus in sensory coding. However, the generalizability of this encoding strategy to the human brain has remained unclear. We asked if individual features of visual stimuli were encoded over distinct timescales. We applied a multiscale time-resolved decoding method to electroencephalography (EEG) collected from human subjects presented with grating visual stimuli to estimate the timescale of individual stimulus features. We observed that the orientation and color of the stimuli were encoded in shorter timescales, whereas spatial frequency and the contrast of the same stimuli were encoded in longer timescales. The stimulus features appeared in temporally overlapping windows along the trial supporting a multiplexed coding strategy. These results provide evidence for a multiplexed, multiscale coding strategy in the human visual system.

Cortical timescales and the modular organization of structural and functional brain networks

Recent years have seen growing interest in characterizing the properties of regional brain dynamics and their relationship to other features of brain structure and function. In particular, multiple studies have observed regional differences in the "timescale" over which activity fluctuates during periods of quiet rest. In the cerebral cortex, these timescales have been associated with both local circuit properties as well as patterns of inter-regional connectivity, including the extent to which each region exhibits widespread connectivity to other brain areas. In the current study, we build on prior observations of an association between connectivity and dynamics in the cerebral cortex by investigating the relationship between BOLD fMRI timescales and the modular organization of structural and functional brain networks. We characterize network community structure across multiple scales and find that longer timescales are associated with greater within-community functional connectivity and diverse structural connectivity. We also replicate prior observations of a positive correlation between timescales and structural connectivity degree. Finally, we find evidence for preferential functional connectivity between cortical areas with similar timescales. We replicate these findings in an independent dataset. These results contribute to our understanding of functional brain organization and structure-function relationships in the human brain, and support the notion that regional differences in cortical dynamics may in part reflect the topological role of each region within macroscale brain networks.

Towards a biologically annotated brain connectome

The brain is a network of interleaved neural circuits. In modern connectomics, brain connectivity is typically encoded as a network of nodes and edges, abstracting away the rich biological detail of local neuronal populations. Yet biological annotations for network nodes - such as gene expression, cytoarchitecture, neurotransmitter receptors or intrinsic dynamics - can be readily measured and overlaid on network models. Here we review how connectomes can be represented and analysed as annotated networks. Annotated connectomes allow us to reconceptualize architectural features of networks and to relate the connection patterns of brain regions to their underlying biology. Emerging work demonstrates that annotated connectomes help to make more veridical models of brain network formation, neural dynamics and disease propagation. Finally, annotations can be used to infer entirely new inter-regional relationships and to construct new types of network that complement existing connectome representations. In summary, biologically annotated connectomes offer a compelling way to study neural wiring in concert with local biological features.

Decreased intrinsic neural timescales in obsessive compulsive disorder and two distinct subtypes revealed by heterogeneity through discriminative analysis

BACKGROUND

OCD is featured as the destruction of information storage and processing. The cognition of neurobiological and clinical heterogeneity is in suspense and poorly studied.

METHODS

Ninety-nine patients and matched HCs(n = 104) were recruited and underwent resting-state functional MRI scans. We applied INT to evaluate altered local neural dynamics representing the ability of information integration. Moreover, considering OCD was a highly heterogeneous disorder, we investigated putative OCD subtypes from INT using a novel semi-supervised machine learning, named HYDRA.

RESULTS

Compared with HCs, patients with OCD showed decreased INTs in extensive brain regions, including bilateral cerebellum and precuneus, STG/MTG and PCC, hippocampus in DMN; right IFG/MFG/SFG, SPL and bilateral angular gyrus in CEN and insula, SMA in SN. Moreover, many other regions involved in visual processing also had disrupted dynamics of local neural organization, consisting of bilateral CUN, LING and fusiform gyrus and occipital lobe. HYDRA divided patients into two distinct neuroanatomical subtypes from INT. Subtype 1 showed decreased INTs in distributed networks, while subtype 2 presented increased in several common regions which were also found to be decreased in subtype 1, such as STG, IPL, postcentral gyrus and left insula, supramarginal gyrus.

CONCLUSION

This study showed distinct abnormalities from the perspective of dynamics of local neural organization in OCD. Such alteration and dimensional approach may provide a new insight into the prior traditional cognition of this disorder and to some extent do favor of more precise diagnosis and treatment response in the future.

Neuromodulation of striatal D1 cells shapes BOLD fluctuations in anatomically connected thalamic and cortical regions

Understanding how the brain's macroscale dynamics are shaped by underlying microscale mechanisms is a key problem in neuroscience. In animal models, we can now investigate this relationship in unprecedented detail by directly manipulating cellular-level properties while measuring the whole-brain response using resting-state fMRI. Here, we focused on understanding how blood-oxygen-level-dependent (BOLD) dynamics, measured within a structurally well-defined striato-thalamo-cortical circuit in mice, are shaped by chemogenetically exciting or inhibiting D1 medium spiny neurons (MSNs) of the right dorsomedial caudate putamen (CPdm). We characterize changes in both the BOLD dynamics of individual cortical and subcortical brain areas, and patterns of inter-regional coupling (functional connectivity) between pairs of areas. Using a classification approach based on a large and diverse set of time-series properties, we found that CPdm neuromodulation alters BOLD dynamics within thalamic subregions that project back to dorsomedial striatum. In the cortex, changes in local dynamics were strongest in unimodal regions (which process information from a single sensory modality) and weakened along a hierarchical gradient towards transmodal regions. In contrast, a decrease in functional connectivity was observed only for cortico-striatal connections after D1 excitation. Our results show that targeted cellular-level manipulations affect local BOLD dynamics at the macroscale, such as by making BOLD dynamics more predictable over time by increasing its self-correlation structure. This contributes to ongoing attempts to understand the influence of structure-function relationships in shaping inter-regional communication at subcortical and cortical levels.

Neurophysiological signatures of cortical micro-architecture

Systematic spatial variation in micro-architecture is observed across the cortex. These micro-architectural gradients are reflected in neural activity, which can be captured by neurophysiological time-series. How spontaneous neurophysiological dynamics are organized across the cortex and how they arise from heterogeneous cortical micro-architecture remains unknown. Here we extensively profile regional neurophysiological dynamics across the human brain by estimating over 6800 time-series features from the resting state magnetoencephalography (MEG) signal. We then map regional time-series profiles to a comprehensive multi-modal, multi-scale atlas of cortical micro-architecture, including microstructure, metabolism, neurotransmitter receptors, cell types and laminar differentiation. We find that the dominant axis of neurophysiological dynamics reflects characteristics of power spectrum density and linear correlation structure of the signal, emphasizing the importance of conventional features of electromagnetic dynamics while identifying additional informative features that have traditionally received less attention. Moreover, spatial variation in neurophysiological dynamics is co-localized with multiple micro-architectural features, including gene expression gradients, intracortical myelin, neurotransmitter receptors and transporters, and oxygen and glucose metabolism. Collectively, this work opens new avenues for studying the anatomical basis of neural activity.

Resting-state fMRI signals contain spectral signatures of local hemodynamic response timing

Functional magnetic resonance imaging (fMRI) has proven to be a powerful tool for noninvasively measuring human brain activity; yet, thus far, fMRI has been relatively limited in its temporal resolution. A key challenge is understanding the relationship between neural activity and the blood-oxygenation-level-dependent (BOLD) signal obtained from fMRI, generally modeled by the hemodynamic response function (HRF). The timing of the HRF varies across the brain and individuals, confounding our ability to make inferences about the timing of the underlying neural processes. Here, we show that resting-state fMRI signals contain information about HRF temporal dynamics that can be leveraged to understand and characterize variations in HRF timing across both cortical and subcortical regions. We found that the frequency spectrum of resting-state fMRI signals significantly differs between voxels with fast versus slow HRFs in human visual cortex. These spectral differences extended to subcortex as well, revealing significantly faster hemodynamic timing in the lateral geniculate nucleus of the thalamus. Ultimately, our results demonstrate that the temporal properties of the HRF impact the spectral content of resting-state fMRI signals and enable voxel-wise characterization of relative hemodynamic response timing. Furthermore, our results show that caution should be used in studies of resting-state fMRI spectral properties, because differences in fMRI frequency content can arise from purely vascular origins. This finding provides new insight into the temporal properties of fMRI signals across voxels, which is crucial for accurate fMRI analyses, and enhances the ability of fast fMRI to identify and track fast neural dynamics.

Functional brain networks reflect spatial and temporal autocorrelation

High-throughput experimental methods in neuroscience have led to an explosion of techniques for measuring complex interactions and multi-dimensional patterns. However, whether sophisticated measures of emergent phenomena can be traced back to simpler, low-dimensional statistics is largely unknown. To explore this question, we examined resting-state functional magnetic resonance imaging (rs-fMRI) data using complex topology measures from network neuroscience. Here we show that spatial and temporal autocorrelation are reliable statistics that explain numerous measures of network topology. Surrogate time series with subject-matched spatial and temporal autocorrelation capture nearly all reliable individual and regional variation in these topology measures. Network topology changes during aging are driven by spatial autocorrelation, and multiple serotonergic drugs causally induce the same topographic change in temporal autocorrelation. This reductionistic interpretation of widely used complexity measures may help link them to neurobiology.

Intrinsic timescales in the visual cortex change with selective attention and reflect spatial connectivity

Intrinsic timescales characterize dynamics of endogenous fluctuations in neural activity. Variation of intrinsic timescales across the neocortex reflects functional specialization of cortical areas, but less is known about how intrinsic timescales change during cognitive tasks. We measured intrinsic timescales of local spiking activity within columns of area V4 in male monkeys performing spatial attention tasks. The ongoing spiking activity unfolded across at least two distinct timescales, fast and slow. The slow timescale increased when monkeys attended to the receptive fields location and correlated with reaction times. By evaluating predictions of several network models, we found that spatiotemporal correlations in V4 activity were best explained by the model in which multiple timescales arise from recurrent interactions shaped by spatially arranged connectivity, and attentional modulation of timescales results from an increase in the efficacy of recurrent interactions. Our results suggest that multiple timescales may arise from the spatial connectivity in the visual cortex and flexibly change with the cognitive state due to dynamic effective interactions between neurons.

Resting-state fMRI signals contain spectral signatures of local hemodynamic response timing

Functional magnetic resonance imaging (fMRI) has proven to be a powerful tool for noninvasively measuring human brain activity; yet, thus far, fMRI has been relatively limited in its temporal resolution. A key challenge is understanding the relationship between neural activity and the blood-oxygenation-level-dependent (BOLD) signal obtained from fMRI, generally modeled by the hemodynamic response function (HRF). The timing of the HRF varies across the brain and individuals, confounding our ability to make inferences about the timing of the underlying neural processes. Here we show that resting-state fMRI signals contain information about HRF temporal dynamics that can be leveraged to understand and characterize variations in HRF timing across both cortical and subcortical regions. We found that the frequency spectrum of resting-state fMRI signals significantly differs between voxels with fast versus slow HRFs in human visual cortex. These spectral differences extended to subcortex as well, revealing significantly faster hemodynamic timing in the lateral geniculate nucleus of the thalamus. Ultimately, our results demonstrate that the temporal properties of the HRF impact the spectral content of resting-state fMRI signals and enable voxel-wise characterization of relative hemodynamic response timing. Furthermore, our results show that caution should be used in studies of resting-state fMRI spectral properties, as differences can arise from purely vascular origins. This finding provides new insight into the temporal properties of fMRI signals across voxels, which is crucial for accurate fMRI analyses, and enhances the ability of fast fMRI to identify and track fast neural dynamics.

Neurophysiological signatures of cortical micro-architecture

Systematic spatial variation in micro-architecture is observed across the cortex. These micro-architectural gradients are reflected in neural activity, which can be captured by neurophysiological time-series. How spontaneous neurophysiological dynamics are organized across the cortex and how they arise from heterogeneous cortical micro-architecture remains unknown. Here we extensively profile regional neurophysiological dynamics across the human brain by estimating over 6 800 timeseries features from the resting state magnetoencephalography (MEG) signal. We then map regional time-series profiles to a comprehensive multi-modal, multi-scale atlas of cortical micro-architecture, including microstructure, metabolism, neurotransmitter receptors, cell types and laminar differentiation. We find that the dominant axis of neurophysiological dynamics reflects characteristics of power spectrum density and linear correlation structure of the signal, emphasizing the importance of conventional features of electromagnetic dynamics while identifying additional informative features that have traditionally received less attention. Moreover, spatial variation in neurophysiological dynamics is colocalized with multiple micro-architectural features, including genomic gradients, intracortical myelin, neurotransmitter receptors and transporters, and oxygen and glucose metabolism. Collectively, this work opens new avenues for studying the anatomical basis of neural activity.

Disrupted relationship between intrinsic neural timescales and alpha peak frequency during unconscious states - A high-density EEG study

Our brain processes the different timescales of our environment's temporal input stochastics. Is such a temporal input processing mechanism key for consciousness? To address this research question, we calculated measures of input processing on shorter (alpha peak frequency, APF) and longer (autocorrelation window, ACW) timescales on resting-state high-density EEG (256 channels) recordings and compared them across different consciousness levels (awake/conscious, ketamine and sevoflurane anaesthesia, unresponsive wakefulness, minimally conscious state). We replicate and extend previous findings of: (i) significantly longer ACW values, consistently over all states of unconsciousness, as measured with ACW-0 (an unprecedented longer version of the well-know ACW-50); (ii) significantly slower APF values, as measured with frequency sliding, in all four unconscious states. Most importantly, we report a highly significant correlation of ACW-0 and APF in the conscious state, while their relationship is disrupted in the unconscious states. In sum, we demonstrate the relevance of the brain's capacity for input processing on shorter (APF) and longer (ACW) timescales - including their relationship - for consciousness. Albeit indirectly, e.g., through the analysis of electrophysiological activity at rest, this supports the mechanism of temporo-spatial alignment to the environment's temporal input stochastics, through relating different neural timescales, as one key predisposing factor of consciousness.

Modelling brain dynamics by Boolean networks

Understanding the relationship between brain architecture and brain function is a central issue in neuroscience. We modeled realistic spatio-temporal patterns of brain activity on a human connectome with a Boolean networks model with the aim of computationally replicating certain cognitive functions as they emerge from the standardization of many fMRI studies, identified as patterns of human brain activity. Results from the analysis of simulation data, carried out for different parameters and initial conditions identified many possible paths in the space of parameters of these network models, with normal (ordered asymptotically constant patterns), chaotic (oscillating or disordered) but also highly organized configurations, with countless spatial–temporal patterns. We interpreted these results as routes to chaos, permanence of the systems in regimes of complexity, and ordered stationary behavior, associating these dynamics to cognitive processes. The most important result of this work is the study of emergent neural circuits, i.e., configurations of areas that synchronize over time, both locally and globally, determining the emergence of computational analogues of cognitive processes, which may or may not be similar to the functioning of biological brain. Furthermore, results put in evidence the creation of how the brain creates structures of remote communication. These structures have hierarchical organization, where each level allows for the emergence of brain organizations which behave at the next superior level. Taken together these results allow the interplay of dynamical and topological roots of the multifaceted brain dynamics to be understood.

Changes in dynamic transitions between integrated and segregated states underlie visual hallucinations in Parkinson's disease

Hallucinations are a core feature of psychosis and common in Parkinson's. Their transient, unexpected nature suggests a change in dynamic brain states, but underlying causes are unknown. Here, we examine temporal dynamics and underlying structural connectivity in Parkinson's-hallucinations using a combination of functional and structural MRI, network control theory, neurotransmitter density and genetic analyses. We show that Parkinson's-hallucinators spent more time in a predominantly Segregated functional state with fewer between-state transitions. The transition from integrated-to-segregated state had lower energy cost in Parkinson's-hallucinators; and was therefore potentially preferable. The regional energy needed for this transition was correlated with regional neurotransmitter density and gene expression for serotoninergic, GABAergic, noradrenergic and cholinergic, but not dopaminergic, receptors. We show how the combination of neurochemistry and brain structure jointly shape functional brain dynamics leading to hallucinations and highlight potential therapeutic targets by linking these changes to neurotransmitter systems involved in early sensory and complex visual processing.

Priming cardiovascular exercise improves complex motor skill learning by affecting the trajectory of learning-related brain plasticity

In recent years, mounting evidence from animal models and studies in humans has accumulated for the role of cardiovascular exercise (CE) in improving motor performance and learning. Both CE and motor learning may induce highly dynamic structural and functional brain changes, but how both processes interact to boost learning is presently unclear. Here, we hypothesized that subjects receiving CE would show a different pattern of learning-related brain plasticity compared to non-CE controls, which in turn associates with improved motor learning. To address this issue, we paired CE and motor learning sequentially in a randomized controlled trial with healthy human participants. Specifically, we compared the effects of a 2-week CE intervention against a non-CE control group on subsequent learning of a challenging dynamic balancing task (DBT) over 6 consecutive weeks. Structural and functional MRI measurements were conducted at regular 2-week time intervals to investigate dynamic brain changes during the experiment. The trajectory of learning-related changes in white matter microstructure beneath parieto-occipital and primary sensorimotor areas of the right hemisphere differed between the CE vs. non-CE groups, and these changes correlated with improved learning of the CE group. While group differences in sensorimotor white matter were already present immediately after CE and persisted during DBT learning, parieto-occipital effects gradually emerged during motor learning. Finally, we found that spontaneous neural activity at rest in gray matter spatially adjacent to white matter findings was also altered, therefore indicating a meaningful link between structural and functional plasticity. Collectively, these findings may lead to a better understanding of the neural mechanisms mediating the CE-learning link within the brain.

Intrinsic neural timescales: temporal integration and segregation

We are continuously bombarded by external inputs of various timescales from the environment. How does the brain process this multitude of timescales? Recent resting state studies show a hierarchy of intrinsic neural timescales (INT) with a shorter duration in unimodal regions (e.g., visual cortex and auditory cortex) and a longer duration in transmodal regions (e.g., default mode network). This unimodal-transmodal hierarchy is present across acquisition modalities [electroencephalogram (EEG)/magnetoencephalogram (MEG) and fMRI] and can be found in different species and during a variety of different task states. Together, this suggests that the hierarchy of INT is central to the temporal integration (combining successive stimuli) and segregation (separating successive stimuli) of external inputs from the environment, leading to temporal segmentation and prediction in perception and cognition.

Reduced intrinsic neural timescales in schizophrenia along posterior parietal and occipital areas

We computed intrinsic neural timescales (INT) based on resting-state functional magnetic resonance imaging (rsfMRI) data of healthy controls (HC) and patients with schizophrenia spectrum disorder (SZ) from three independently collected samples. Five clusters showed decreased INT in SZ compared to HC in all three samples: right occipital fusiform gyrus (rOFG), left superior occipital gyrus (lSOG), right superior occipital gyrus (rSOG), left lateral occipital cortex (lLOC) and right postcentral gyrus (rPG). In other words, it appears that sensory information in visual and posterior parietal areas is stored for reduced lengths of time in SZ compared to HC. Finally, we found that symptom severity appears to modulate INT of these areas in SZ.

Modular origins of high-amplitude cofluctuations in fine-scale functional connectivity dynamics

The topology of structural brain networks shapes brain dynamics, including the correlation structure of brain activity (functional connectivity) as estimated from functional neuroimaging data. Empirical studies have shown that functional connectivity fluctuates over time, exhibiting patterns that vary in the spatial arrangement of correlations among segregated functional systems. Recently, an exact decomposition of functional connectivity into frame-wise contributions has revealed fine-scale dynamics that are punctuated by brief and intermittent episodes (events) of high-amplitude cofluctuations involving large sets of brain regions. Their origin is currently unclear. Here, we demonstrate that similar episodes readily appear in silico using computational simulations of whole-brain dynamics. As in empirical data, simulated events contribute disproportionately to long-time functional connectivity, involve recurrence of patterned cofluctuations, and can be clustered into distinct families. Importantly, comparison of event-related patterns of cofluctuations to underlying patterns of structural connectivity reveals that modular organization present in the coupling matrix shapes patterns of event-related cofluctuations. Our work suggests that brief, intermittent events in functional dynamics are partly shaped by modular organization of structural connectivity.

High-level cognition during story listening is reflected in high-order dynamic correlations in neural activity patterns

Our thoughts arise from coordinated patterns of interactions between brain structures that change with our ongoing experiences. High-order dynamic correlations in neural activity patterns reflect different subgraphs of the brain's functional connectome that display homologous lower-level dynamic correlations. Here we test the hypothesis that high-level cognition is reflected in high-order dynamic correlations in brain activity patterns. We develop an approach to estimating high-order dynamic correlations in timeseries data, and we apply the approach to neuroimaging data collected as human participants either listen to a ten-minute story or listen to a temporally scrambled version of the story. We train across-participant pattern classifiers to decode (in held-out data) when in the session each neural activity snapshot was collected. We find that classifiers trained to decode from high-order dynamic correlations yield the best performance on data collected as participants listened to the (unscrambled) story. By contrast, classifiers trained to decode data from scrambled versions of the story yielded the best performance when they were trained using first-order dynamic correlations or non-correlational activity patterns. We suggest that as our thoughts become more complex, they are reflected in higher-order patterns of dynamic network interactions throughout the brain.

Network Controllability in Transmodal Cortex Predicts Positive Psychosis Spectrum Symptoms

BACKGROUND

The psychosis spectrum (PS) is associated with structural dysconnectivity concentrated in transmodal cortex. However, understanding of this pathophysiology has been limited by an overreliance on examining direct interregional connectivity. Using network control theory, we measured variation in both direct and indirect connectivity to a region to gain new insights into the pathophysiology of the PS.

METHODS

We used psychosis symptom data and structural connectivity in 1068 individuals from the Philadelphia Neurodevelopmental Cohort. Applying a network control theory metric called average controllability, we estimated each brain region's capacity to leverage its direct and indirect structural connections to control linear brain dynamics. Using nonlinear regression, we determined the accuracy with which average controllability could predict PS symptoms in out-of-sample testing. We also examined the predictive performance of regional strength, which indexes only direct connections to a region, as well as several graph-theoretic measures of centrality that index indirect connectivity. Finally, we assessed how the prediction performance for PS symptoms varied over the functional hierarchy spanning unimodal to transmodal cortex.

RESULTS

Average controllability outperformed all other connectivity features at predicting positive PS symptoms and was the only feature to yield above-chance predictive performance. Improved prediction for average controllability was concentrated in transmodal cortex, whereas prediction performance for strength was uniform across the cortex, suggesting that indexing indirect connections through average controllability is crucial in association cortex.

CONCLUSIONS

Examining interindividual variation in direct and indirect structural connections to transmodal cortex is crucial for accurate prediction of positive PS symptoms.

The brain and its time: intrinsic neural timescales are key for input processing

We process and integrate multiple timescales into one meaningful whole. Recent evidence suggests that the brain displays a complex multiscale temporal organization. Different regions exhibit different timescales as described by the concept of intrinsic neural timescales (INT); however, their function and neural mechanisms remains unclear. We review recent literature on INT and propose that they are key for input processing. Specifically, they are shared across different species, i.e., input sharing. This suggests a role of INT in encoding inputs through matching the inputs' stochastics with the ongoing temporal statistics of the brain's neural activity, i.e., input encoding. Following simulation and empirical data, we point out input integration versus segregation and input sampling as key temporal mechanisms of input processing. This deeply grounds the brain within its environmental and evolutionary context. It carries major implications in understanding mental features and psychiatric disorders, as well as going beyond the brain in integrating timescales into artificial intelligence.

Dynamical consequences of regional heterogeneity in the brain's transcriptional landscape

Brain regions vary in their molecular and cellular composition, but how this heterogeneity shapes neuronal dynamics is unclear. Here, we investigate the dynamical consequences of regional heterogeneity using a biophysical model of whole-brain functional magnetic resonance imaging (MRI) dynamics in humans. We show that models in which transcriptional variations in excitatory and inhibitory receptor (E:I) gene expression constrain regional heterogeneity more accurately reproduce the spatiotemporal structure of empirical functional connectivity estimates than do models constrained by global gene expression profiles or MRI-derived estimates of myeloarchitecture. We further show that regional transcriptional heterogeneity is essential for yielding both ignition-like dynamics, which are thought to support conscious processing, and a wide variance of regional-activity time scales, which supports a broad dynamical range. We thus identify a key role for E:I heterogeneity in generating complex neuronal dynamics and demonstrate the viability of using transcriptomic data to constrain models of large-scale brain function.

Topographic gradients of intrinsic dynamics across neocortex

The intrinsic dynamics of neuronal populations are shaped by both microscale attributes and macroscale connectome architecture. Here we comprehensively characterize the rich temporal patterns of neural activity throughout the human brain. Applying massive temporal feature extraction to regional haemodynamic activity, we systematically estimate over 6000 statistical properties of individual brain regions' time-series across the neocortex. We identify two robust spatial gradients of intrinsic dynamics, one spanning a ventromedial-dorsolateral axis and dominated by measures of signal autocorrelation, and the other spanning a unimodal-transmodal axis and dominated by measures of dynamic range. These gradients reflect spatial patterns of gene expression, intracortical myelin and cortical thickness, as well as structural and functional network embedding. Importantly, these gradients are correlated with patterns of meta-analytic functional activation, differentiating cognitive versus affective processing and sensory versus higher-order cognitive processing. Altogether, these findings demonstrate a link between microscale and macroscale architecture, intrinsic dynamics, and cognition.

Individual differences in haemoglobin concentration influence bold fMRI functional connectivity and its correlation with cognition

Resting-state connectivity measures the temporal coherence of the spontaneous neural activity of spatially distinct regions, and is commonly measured using BOLD-fMRI. The BOLD response follows neuronal activity, when changes in the relative concentration of oxygenated and deoxygenated haemoglobin cause fluctuations in the MRI T2* signal. Since the BOLD signal detects changes in relative concentrations of oxy/deoxy-haemoglobin, individual differences in haemoglobin levels may influence the BOLD signal-to-noise ratio in a manner independent of the degree of neural activity. In this study, we examined whether group differences in haemoglobin may confound measures of functional connectivity. We investigated whether relationships between measures of functional connectivity and cognitive performance could be influenced by individual variability in haemoglobin. Finally, we mapped the neuroanatomical distribution of the influence of haemoglobin on functional connectivity to determine where group differences in functional connectivity are manifest. In a cohort of 518 healthy elderly subjects (259 men), each sex group was median-split into two groups with high and low haemoglobin concentration. Significant differences were obtained in functional connectivity between the high and low haemoglobin groups for both men and women (Cohen's d 0.17 and 0.03 for men and women respectively). The haemoglobin connectome in males showed a widespread systematic increase in functional connectivity correlation values, whilst the female connectome showed predominantly parietal and subcortical increases and temporo-parietal decreases. Despite the haemoglobin groups having no differences in cognitive measures, significant differences in the linear relationships between cognitive performance and functional connectivity were obtained for all 5 cognitive tests in males, and 4 out of 5 tests in females. Our findings confirm that individual variability in haemoglobin levels that give rise to group differences are an important confounding variable in BOLD-fMRI-based studies of functional connectivity. Controlling for haemoglobin variability as a potentially confounding variable is crucial to ensure the reproducibility of human brain connectome studies, especially in studies that compare groups of individuals, compare sexes, or examine connectivity-cognition relationships.