Robustness of connectome harmonics to local gray matter and long-range white matter connectivity changes

Time-resolved coupling between connectome harmonics and subjective experience under the psychedelic DMT

Exploring the intricate relationship between brain's structure and function, and how this affects subjective experience is a fundamental pursuit in neuroscience. Psychedelic substances offer a unique insight into the influences of specific neurotransmitter systems on perception, cognition and consciousness. Specifically, their impact on brain function propagates across the structural connectome - a network of white matter pathways linking different regions. To comprehensively grasp the effects of psychedelic compounds on brain function, we used a theoretically rigorous framework known as connectome harmonic decomposition. This framework provides a robust method to characterize how brain function intricately depends on the organized network structure of the human connectome. We show that the connectome harmonic repertoire under DMT is reshaped in line with other reported psychedelic compounds - psilocybin, LSD and ketamine. Furthermore, we show that the repertoire entropy of connectome harmonics increases under DMT, as with those other psychedelics. Importantly, we demonstrate for the first time that measures of energy spectrum difference and repertoire entropy of connectome harmonics indexes the intensity of subjective experience of the participants in a time-resolved manner reflecting close coupling between connectome harmonics and subjective experience.

Connectome spectrum electromagnetic tomography: A method to reconstruct electrical brain source networks at high-spatial resolution

Connectome spectrum electromagnetic tomography (CSET) combines diffusion MRI-derived structural connectivity data with well-established graph signal processing tools to solve the M/EEG inverse problem. Using simulated EEG signals from fMRI responses, and two EEG datasets on visual-evoked potentials, we provide evidence supporting that (i) CSET captures realistic neurophysiological patterns with better accuracy than state-of-the-art methods, (ii) CSET can reconstruct brain responses more accurately and with more robustness to intrinsic noise in the EEG signal. These results demonstrate that CSET offers high spatio-temporal accuracy, enabling neuroscientists to extend their research beyond the current limitations of low sampling frequency in functional MRI and the poor spatial resolution of M/EEG.

Decomposing cortical activity through neuronal tracing connectome-eigenmodes in marmosets

Deciphering the complex relationship between neuroanatomical connections and functional activity in primate brains remains a daunting task, especially regarding the influence of monosynaptic connectivity on cortical activity. Here, we investigate the anatomical-functional relationship and decompose the neuronal-tracing connectome of marmoset brains into a series of eigenmodes using graph signal processing. These cellular connectome eigenmodes effectively constrain the cortical activity derived from resting-state functional MRI, and uncover a patterned cellular-functional decoupling. This pattern reveals a spatial gradient from coupled dorsal-posterior to decoupled ventral-anterior cortices, and recapitulates micro-structural profiles and macro-scale hierarchical cortical organization. Notably, these marmoset-derived eigenmodes may facilitate the inference of spontaneous cortical activity and functional connectivity of homologous areas in humans, highlighting the potential generalizing of the connectomic constraints across species. Collectively, our findings illuminate how neuronal-tracing connectome eigenmodes constrain cortical activity and improve our understanding of the brain's anatomical-functional relationship.

Mode-based morphometry: A multiscale approach to mapping human neuroanatomy

Voxel-based morphometry (VBM) and surface-based morphometry (SBM) are two widely used neuroimaging techniques for investigating brain anatomy. These techniques rely on statistical inferences at individual points (voxels or vertices), clusters of points, or a priori regions-of-interest. They are powerful tools for describing brain anatomy, but offer little insights into the generative processes that shape a particular set of findings. Moreover, they are restricted to a single spatial resolution scale, precluding the opportunity to distinguish anatomical variations that are expressed across multiple scales. Drawing on concepts from classical physics, here we develop an approach, called mode-based morphometry (MBM), that can describe any empirical map of anatomical variations in terms of the fundamental, resonant modes-eigenmodes-of brain anatomy, each tied to a specific spatial scale. Hence, MBM naturally yields a multiscale characterization of the empirical map, affording new opportunities for investigating the spatial frequency content of neuroanatomical variability. Using simulated and empirical data, we show that the validity and reliability of MBM are either comparable or superior to classical vertex-based SBM for capturing differences in cortical thickness maps between two experimental groups. Our approach thus offers a robust, accurate, and informative method for characterizing empirical maps of neuroanatomical variability that can be directly linked to a generative physical process.

Spatial cell-type enrichment predicts mouse brain connectivity

A fundamental neuroscience topic is the link between the brain's molecular, cellular, and cytoarchitectonic properties and structural connectivity. Recent studies relate inter-regional connectivity to gene expression, but the relationship to regional cell-type distributions remains understudied. Here, we utilize whole-brain mapping of neuronal and non-neuronal subtypes via the matrix inversion and subset selection algorithm to model inter-regional connectivity as a function of regional cell-type composition with machine learning. We deployed random forest algorithms for predicting connectivity from cell-type densities, demonstrating surprisingly strong prediction accuracy of cell types in general, and particular non-neuronal cells such as oligodendrocytes. We found evidence of a strong distance dependency in the cell connectivity relationship, with layer-specific excitatory neurons contributing the most for long-range connectivity, while vascular and astroglia were salient for short-range connections. Our results demonstrate a link between cell types and connectivity, providing a roadmap for examining this relationship in other species, including humans.

The coupling between the spatial and temporal scales of neural processes revealed by a joint time-vertex connectome spectral analysis

Brain oscillations are produced by the coordinated activity of large groups of neurons and different rhythms are thought to reflect different modes of information processing. These modes, in turn, are known to occur at different spatial scales. Nevertheless, how these rhythms support different spatial modes of information processing at the brain scale is not yet fully understood. Here we use "Joint Time-Vertex Spectral Analysis" to characterize the joint spectral content of brain activity both in time (temporal frequencies) and in space over the connectivity graph (spatial connectome harmonics). This method allows us to characterize the relationship between spatially localized or distributed neural processes on one side and their respective temporal frequency bands in source-reconstructed M/EEG signals. We explore this approach on two different datasets, an auditory steady-state response (ASSR) and a visual grating task. Our results suggest that different information processing mechanisms are carried out at different frequency bands: while spatially distributed activity (which may also be interpreted as integration) specifically occurs at low temporal frequencies (alpha and theta) and low graph spatial frequencies, localized electrical activity (i.e., segregation) is observed at high temporal frequencies (high and low gamma) over restricted high spatial graph frequencies. Crucially, the estimated contribution of the distributed and localized neural activity predicts performance in a behavioral task, demonstrating the neurophysiological relevance of the joint time-vertex spectral representation.

Geometric constraints on human brain function

The anatomy of the brain necessarily constrains its function, but precisely how remains unclear. The classical and dominant paradigm in neuroscience is that neuronal dynamics are driven by interactions between discrete, functionally specialized cell populations connected by a complex array of axonal fibres1-3. However, predictions from neural field theory, an established mathematical framework for modelling large-scale brain activity4-6, suggest that the geometry of the brain may represent a more fundamental constraint on dynamics than complex interregional connectivity7,8. Here, we confirm these theoretical predictions by analysing human magnetic resonance imaging data acquired under spontaneous and diverse task-evoked conditions. Specifically, we show that cortical and subcortical activity can be parsimoniously understood as resulting from excitations of fundamental, resonant modes of the brain's geometry (that is, its shape) rather than from modes of complex interregional connectivity, as classically assumed. We then use these geometric modes to show that task-evoked activations across over 10,000 brain maps are not confined to focal areas, as widely believed, but instead excite brain-wide modes with wavelengths spanning over 60 mm. Finally, we confirm predictions that the close link between geometry and function is explained by a dominant role for wave-like activity, showing that wave dynamics can reproduce numerous canonical spatiotemporal properties of spontaneous and evoked recordings. Our findings challenge prevailing views and identify a previously underappreciated role of geometry in shaping function, as predicted by a unifying and physically principled model of brain-wide dynamics.

Estimating Outlier-Immunized Common Harmonic Waves for Brain Network Analyses on the Stiefel Manifold

Since brain network organization is essentially governed by the harmonic waves derived from the Eigen-system of the underlying Laplacian matrix, discovering the harmonic-based alterations provides a new window to understand the pathogenic mechanism of Alzheimer's disease (AD) in a unified reference space. However, current reference (common harmonic waves) estimation studies over the individual harmonic waves are often sensitive to outliers, which are obtained by averaging the heterogenous individual brain networks. To address this challenge, we propose a novel manifold learning approach to identify a set of outlier-immunized common harmonic waves. The backbone of our framework is calculating the geometric median of all individual harmonic waves on the Stiefel manifold, instead of Fréchet mean, thus improving the robustness of learned common harmonic waves to the outliers. A manifold optimization scheme with theoretically guaranteed convergence is tailored to solve our method. The experimental results on synthetic data and real data demonstrate that the common harmonic waves learned by our approach are not only more robust to the outliers than the state-of-the-art methods, but also provide a putative imaging biomarker to predict the early stage of AD.

Mode-based morphometry: A multiscale approach to mapping human neuroanatomy

Voxel-based morphometry (VBM) and surface-based morphometry (SBM) are two widely used neuroimaging techniques for investigating brain anatomy. These techniques rely on statistical inferences at individual points (voxels or vertices), clusters of points, or a priori regions-of-interest. They are powerful tools for describing brain anatomy, but offer little insights into the generative processes that shape a particular set of findings. Moreover, they are restricted to a single spatial resolution scale, precluding the opportunity to distinguish anatomical variations that are expressed across multiple scales. Drawing on concepts from classical physics, here we develop an approach, called mode-based morphometry (MBM), that can describe any empirical map of anatomical variations in terms of the fundamental, resonant modes--eigenmodes--of brain anatomy, each tied to a specific spatial scale. Hence, MBM naturally yields a multiscale characterization of the empirical map, affording new opportunities for investigating the spatial frequency content of neuroanatomical variability. Using simulated and empirical data, we show that the validity and reliability of MBM are either comparable or superior to classical vertex-based SBM for capturing differences in cortical thickness maps between two experimental groups. Our approach thus offers a robust, accurate, and informative method for characterizing empirical maps of neuroanatomical variability that can be directly linked to a generative physical process.

Distributed harmonic patterns of structure-function dependence orchestrate human consciousness

A central question in neuroscience is how consciousness arises from the dynamic interplay of brain structure and function. Here we decompose functional MRI signals from pathological and pharmacologically-induced perturbations of consciousness into distributed patterns of structure-function dependence across scales: the harmonic modes of the human structural connectome. We show that structure-function coupling is a generalisable indicator of consciousness that is under bi-directional neuromodulatory control. We find increased structure-function coupling across scales during loss of consciousness, whether due to anaesthesia or brain injury, capable of discriminating between behaviourally indistinguishable sub-categories of brain-injured patients, tracking the presence of covert consciousness. The opposite harmonic signature characterises the altered state induced by LSD or ketamine, reflecting psychedelic-induced decoupling of brain function from structure and correlating with physiological and subjective scores. Overall, connectome harmonic decomposition reveals how neuromodulation and the network architecture of the human connectome jointly shape consciousness and distributed functional activation across scales.

Disruption in structural-functional network repertoire and time-resolved subcortical fronto-temporoparietal connectivity in disorders of consciousness

Understanding recovery of consciousness and elucidating its underlying mechanism is believed to be crucial in the field of basic neuroscience and medicine. Ideas such as the global neuronal workspace (GNW) and the mesocircuit theory hypothesize that failure of recovery in conscious states coincide with loss of connectivity between subcortical and frontoparietal areas, a loss of the repertoire of functional networks states and metastable brain activation. We adopted a time-resolved functional connectivity framework to explore these ideas and assessed the repertoire of functional network states as a potential marker of consciousness and its potential ability to tell apart patients in the unresponsive wakefulness syndrome (UWS) and minimally conscious state (MCS). In addition, the prediction of these functional network states by underlying hidden spatial patterns in the anatomical network, that is so-called eigenmodes, was supplemented as potential markers. By analysing time-resolved functional connectivity from functional MRI data, we demonstrated a reduction of metastability and functional network repertoire in UWS compared to MCS patients. This was expressed in terms of diminished dwell times and loss of nonstationarity in the default mode network and subcortical fronto-temporoparietal network in UWS compared to MCS patients. We further demonstrated that these findings co-occurred with a loss of dynamic interplay between structural eigenmodes and emerging time-resolved functional connectivity in UWS. These results are, amongst others, in support of the GNW theory and the mesocircuit hypothesis, underpinning the role of time-resolved thalamo-cortical connections and metastability in the recovery of consciousness.

Predicting Functional Connectivity From Observed and Latent Structural Connectivity via Eigenvalue Mapping

Understanding how complex dynamic activity propagates over a static structural network is an overarching question in the field of neuroscience. Previous work has demonstrated that linear graph-theoretic models perform as well as non-linear neural simulations in predicting functional connectivity with the added benefits of low dimensionality and a closed-form solution which make them far less computationally expensive. Here we show a simple model relating the eigenvalues of the structural connectivity and functional networks using the Gamma function, producing a reliable prediction of functional connectivity with a single model parameter. We also investigate the impact of local activity diffusion and long-range interhemispheric connectivity on the structure-function model and show an improvement in functional connectivity prediction when accounting for such latent variables which are often excluded from traditional diffusion tensor imaging (DTI) methods.

Edges in brain networks: Contributions to models of structure and function

Network models describe the brain as sets of nodes and edges that represent its distributed organization. So far, most discoveries in network neuroscience have prioritized insights that highlight distinct groupings and specialized functional contributions of network nodes. Importantly, these functional contributions are determined and expressed by the web of their interrelationships, formed by network edges. Here, we underscore the important contributions made by brain network edges for understanding distributed brain organization. Different types of edges represent different types of relationships, including connectivity and similarity among nodes. Adopting a specific definition of edges can fundamentally alter how we analyze and interpret a brain network. Furthermore, edges can associate into collectives and higher order arrangements, describe time series, and form edge communities that provide insights into brain network topology complementary to the traditional node-centric perspective. Focusing on the edges, and the higher order or dynamic information they can provide, discloses previously underappreciated aspects of structural and functional network organization.

Functional harmonics reveal multi-dimensional basis functions underlying cortical organization

The human brain consists of specialized areas that flexibly interact to form a multitude of functional networks. Complementary to this notion of modular organization, brain function has been shown to vary along a smooth continuum across the whole cortex. We demonstrate a mathematical framework that accounts for both of these perspectives: harmonic modes. We calculate the harmonic modes of the brain's functional connectivity graph, called "functional harmonics," revealing a multi-dimensional, frequency-ordered set of basis functions. Functional harmonics link characteristics of cortical organization across several spatial scales, capturing aspects of intra-areal organizational features (retinotopy, somatotopy), delineating brain areas, and explaining macroscopic functional networks as well as global cortical gradients. Furthermore, we show how the activity patterns elicited by seven different tasks are reconstructed from a very small subset of functional harmonics. Our results suggest that the principle of harmonicity, ubiquitous in nature, also underlies functional cortical organization in the human brain.