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

Selective control of prefrontal neural timescales by parietal cortex

Intrinsic neural timescales quantify how long spontaneous neuronal activity patterns persist, reflecting dynamics of endogenous fluctuations. We measured intrinsic timescales of frontal eye field (FEF) neurons and examined their changes during posterior parietal cortex (PPC) inactivation. We observed two distinct classes of FEF neurons based on their intrinsic timescales: short-timescale neurons (∼25 ms) or long-timescale neurons (∼100 ms). Short-timescale neurons showed stronger transient visual responses, suggesting their role in rapid visual processing, whereas long-timescale neurons exhibited pronounced sustained attentional modulation, suggesting their role in maintaining stimulus-driven attention. During PPC inactivation, intrinsic timescales increased in both neuron types, with a significantly larger effect in short-timescale neurons. In addition, PPC inactivation reduced attentional modulation, particularly in long-timescale neurons. Our findings provide the first causal evidence linking intrinsic local neural timescales to long-range inter-area communications. These findings also suggest the presence of at least two distinct network motifs that support different neuronal dynamics and functional computations within the FEF.

Functionally specialized spectral organization of the resting human cortex

Ample studies across various neuroimaging modalities have suggested that the human cortex at rest is hierarchically organized along the spectral and functional axes. However, the relationship between the spectral and functional organizations of the human cortex remains largely unexplored. Here, we reveal the confluence of functional and spectral cortical organizations by examining the functional specialization in spectral gradients of the cortex. These spectral gradients, derived from functional magnetic resonance imaging data at rest using our temporal de-correlation method to enhance spectral resolution, demonstrate regional frequency biases. The grading of spectral gradients across the cortex - aligns with many existing brain maps - is found to be highly functionally specialized through discovered frequency-specific resting-state functional networks, functionally distinctive spectral profiles, and an intrinsic coordinate system that is functionally specialized. By demonstrating the functionally specialized spectral gradients of the cortex, we shed light on the close relation between functional and spectral organizations of the resting human cortex.

Sleep Modulates Neural Timescales and Spatiotemporal Integration in the Human Cortex

Spontaneous neural dynamics manifest across multiple temporal and spatial scales, which are thought to be intrinsic to brain areas and exhibit hierarchical organization across the cortex. In wake, a hierarchy of timescales is thought to naturally emerge from microstructural properties, gene expression, and recurrent connections. A fundamental question is timescales' organization and changes in sleep, where physiological needs are different. Here, we describe two measures of neural timescales, obtained from broadband activity and gamma power, which display complementary properties. We leveraged intracranial electroencephalography in 106 human epilepsy patients (48 females) to characterize timescale changes from wake to sleep across the cortical hierarchy. We show that both broadband and gamma timescales are globally longer in sleep than in wake. While broadband timescales increase along the sensorimotor-association axis, gamma ones decrease. During sleep, slow waves can explain the increase of broadband and gamma timescales, but only broadband ones show a positive association with slow-wave density across the cortex. Finally, we characterize spatial correlations and their relationship with timescales as a proxy for spatiotemporal integration, finding high integration at long distances in wake for broadband and at short distances in sleep for gamma timescales. Our results suggest that mesoscopic neural populations possess different timescales that are shaped by anatomy and are modulated by the sleep/wake cycle.

Longitudinal multimodal neuroimaging after traumatic brain injury

Traumatic brain injury is a major cause of long-term cognitive impairment, yet the mechanisms underlying recovery remain poorly understood. Neuroimaging methods such as diffusion MRI, functional MRI, and positron emission tomography (PET) provide insight into micro- and macro-scale changes post-TBI, but the relationships between regional cellular and functional alterations remain unclear. In this study, we conducted a longitudinal, multimodal neuroimaging analysis quantifying TBI-related pathologies in four biomarkers, namely flumazenil PET derived binding potential, dMRI-derived structural connectivity, and resting-state fMRI-derived functional connectivity and fractional amplitude of low-frequency fluctuations in individuals with mild-to-severe brain injury at the subacute (4-6 months post-injury) and chronic (1-year postinjury) stages. Brain injury related regional pathologies, and their changes over time, were correlated across the four biomarkers. Our results reveal complex, dynamic changes over time. We found that flumazenil-PET binding potential was significantly reduced in frontal and thalamic regions in brain injured subjects, consistent with neuronal loss, with partial recovery over time. Functional hyperconnectivity was observed in brain injured subjects initially but declined while remaining elevated compared to non-injured controls, whereas cortical structural hypoconnectivity persisted. Importantly, we observed that brain injury related alterations across MRI modalities became more strongly correlated with flumazenil-PET at the chronic stage. Regions with chronic reductions in flumazenil-PET binding also showed weaker structural node strength and lower amplitude of low frequency fluctuations, a relationship that was not found at the subacute stage. This observation could suggest a progressive convergence of structural and functional disruptions with neuronal loss over time. Additionally, regions with declining structural node strength also exhibited decreases in functional node strength, while these same regions showed increased amplitude of low frequency fluctuations over time. This pattern suggests that heightened intrinsic regional activity may serve as a compensatory mechanism in regions increasingly disconnected due to progressive axonal degradation. Altogether, these findings advance our understanding of how multimodal neuroimaging captures the evolving interplay between neuronal integrity, structural connectivity, and functional dynamics after brain injury. Clarifying these interrelationships could inform prognostic models and enhance knowledge of degenerative, compensatory, and recovery mechanisms in traumatic brain injury.

Network Mechanisms Underlying the Regional Diversity of Variance and Time Scales of the Brain's Spontaneous Activity Fluctuations

The brain's activity fluctuations have different temporal scales across the brain regions, with associative regions displaying slower timescales than sensory areas. This hierarchy of timescales has been shown to correlate with both structural brain connectivity and intrinsic regional properties. Here, using publicly available human resting-state fMRI and dMRI data, it was found that, while more structurally connected brain regions presented activity fluctuations with longer timescales, their activity fluctuations presented lower variance. The opposite relationships between the structural connectivity and the variance and temporal scales of resting-state fluctuations, respectively, were not trivially explained by simple network propagation principles. To understand these structure-function relationships, two commonly used whole-brain models were studied, namely, the Hopf and Wilson-Cowan models. These models use the brain's connectome to couple local nodes (representing brain regions) displaying noise-driven oscillations. The models show that the variance and temporal scales of activity fluctuations can oppositely relate to connectivity within specific parameter regions, even when all nodes have the same intrinsic dynamics-but also when intrinsic dynamics are constrained by the myelinization-related macroscopic gradient. These results show that, setting aside intrinsic regional differences, connectivity and network state are sufficient to explain the regional differences in fluctuations' scales. State dependence supports the vision that structure-function relationships can serve as biomarkers of altered brain states. Finally, the results indicate that the hierarchies of timescales and variances reflect a balance between stability and responsivity, with greater and faster responsiveness at the network periphery, while the network core ensures overall system robustness.

Structurally informed models of directed brain connectivity

Understanding how one brain region exerts influence over another in vivo is profoundly constrained by models used to infer or predict directed connectivity. Although such neural interactions rely on the anatomy of the brain, it remains unclear whether, at the macroscale, structural (or anatomical) connectivity provides useful constraints on models of directed connectivity. Here, we review the current state of research on this question, highlighting a key distinction between inference-based effective connectivity and prediction-based directed functional connectivity. We explore the methods via which structural connectivity has been integrated into directed connectivity models: through prior distributions, fixed parameters in state-space models and inputs to structure learning algorithms. Although the evidence suggests that integrating structural connectivity substantially improves directed connectivity models, assessments of reliability and out-of-sample validity are lacking. We conclude this Review with a strategy for future research that addresses current challenges and identifies opportunities for advancing the integration of structural and directed connectivity to ultimately improve understanding of the brain in health and disease.

Decoding cortical chronotopy-Comparing the influence of different cortical organizational schemes

The brain's diverse intrinsic timescales enable us to perceive stimuli with varying temporal persistency. This study aimed to uncover the cortical organizational schemes underlying these variations, revealing the neural architecture for processing a wide range of sensory experiences. We collected resting-state fMRI, task-fMRI, and diffusion-weighted imaging data from 47 individuals. Based on this data, we extracted six organizational schemes: (1) the structural Rich Club (RC) architecture, shown to synchronize the connectome; (2) the structural Diverse Club architecture, as an alternative to the RC based on the network's module structure; (3) the functional uni-to-multimodal gradient, reflected in a wide range of structural and functional features; and (4) the spatial posterior/lateral-to-anterior/medial gradient, established for hierarchical levels of cognitive control. Also, we explored the effects of (5) structural graph theoretical measures of centrality and (6) cytoarchitectural differences. Using Bayesian model comparison, we contrasted the impact of these organizational schemes on (1) intrinsic resting-state timescales and (2) inter-subject correlation (ISC) from a task involving hierarchically nested digit sequences. As expected, resting-state timescales were slower in structural network hubs, hierarchically higher areas defined by the functional and spatial gradients, and thicker cortical regions. ISC analysis demonstrated hints for the engagement of higher cortical areas with more temporally persistent stimuli. Finally, the model comparison identified the uni-to-multimodal gradient as the best organizational scheme for explaining the chronotopy in both task and rest. Future research should explore the microarchitectural features that shape this gradient, elucidating how our brain adapts and evolves across different modes of processing.

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.

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.

Behavioural stochastic resonance across the lifespan

Stochastic resonance (SR) is the phenomenon wherein the introduction of a suitable level of noise enhances the detection of subthreshold signals in non linear systems. It manifests across various physical and biological systems, including the human brain. Psychophysical experiments have confirmed the behavioural impact of stochastic resonance on auditory, somatic, and visual perception. Aging renders the brain more susceptible to noise, possibly causing differences in the  SR phenomenon between young and elderly individuals. This study investigates the impact of noise on motion detection accuracy throughout the lifespan, with 214 participants ranging in age from 18 to 82. Our objective was to determine the optimal noise level to induce an SR-like response in both young and old populations. Consistent with existing literature, our findings reveal a diminishing advantage with age, indicating that the efficacy of noise addition progressively diminishes. Additionally, as individuals age, peak performance is achieved with lower levels of noise. This study provides the first insight into how SR changes across the lifespan of healthy adults and establishes a foundation for understanding the pathological alterations in perceptual processes associated with aging.

Dendrites contribute to the gradient of intrinsic timescales encompassing cortical and subcortical brain networks

Introduction

Cytoarchitectonic studies have uncovered a correlation between higher levels of cortical hierarchy and reduced dendritic size. This hierarchical organization extends to the brain's timescales, revealing longer intrinsic timescales at higher hierarchical levels. However, estimating the contribution of single-neuron dendritic morphology to the hierarchy of timescales, which is typically characterized at a macroscopic level, remains challenging.

Method

Here we mapped the intrinsic timescales of six functional networks using functional magnetic resonance imaging (fMRI) data, and characterized the influence of neuronal dendritic size on intrinsic timescales of brain regions, utilizing a multicompartmental neuronal modeling approach based on digitally reconstructed neurons.

Results

The fMRI results revealed a hierarchy of intrinsic timescales encompassing both cortical and subcortical brain regions. The neuronal modeling indicated that neurons with larger dendritic structures exhibit shorter intrinsic timescales. Together these findings highlight the contribution of dendrites at the neuronal level to the hierarchy of intrinsic timescales at the whole-brain level.

Discussion

This study sheds light on the intricate relationship between neuronal structure, cytoarchitectonic maps, and the hierarchy of timescales in the brain.

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.