A Riemannian approach to predicting brain function from the structural connectome

Elucidating genetic and molecular basis of altered higher-order brain structure-function coupling in major depressive disorder

Previous studies have shown that major depressive disorder (MDD) patients exhibit structural and functional impairments, but few studies have investigated changes in higher-order coupling between structure and function. Here, we systematically investigated the effect of MDD on higher-order coupling between structural connectivity (SC) and functional connectivity (FC). Each brain region was mapped into embedding vector by the node2vec algorithm. We used support vector machine (SVM) with the brain region embedding vector to distinguish MDD patients from health controls (HCs) and identify the most discriminative brain regions. Our study revealed that MDD patients had decreased higher-order coupling in connections between the most discriminative brain regions and local connections in rich-club organization and increased higher-order coupling in connections between the ventral attentional network and limbic network compared with HCs. Interestingly, transcriptome-neuroimaging association analysis demonstrated the correlations between regional rSC-FC coupling variations between MDD patients and HCs and α/β-hydrolase domain-containing 6 (ABHD6), β 1,3-N-acetylglucosaminyltransferase-9(β3GNT9), transmembrane protein 45B (TMEM45B), the correlation between regional dSC-FC coupling variations and retinoic acid early transcript 1E antisense RNA 1(RAET1E-AS1), and the correlations between regional iSC-FC coupling variations and ABHD6, β3GNT9, katanin-like 2 protein (KATNAL2). In addition, correlation analysis with neurotransmitter receptor/transporter maps found that the rSC-FC and iSC-FC coupling variations were both correlated with neuroendocrine transporter (NET) expression, and the dSC-FC coupling variations were correlated with metabotropic glutamate receptor 5 (mGluR5). Further mediation analysis explored the relationship between genes, neurotransmitter receptor/transporter and MDD related higher-order coupling variations. These findings indicate that specific genetic and molecular factors underpin the observed disparities in higher-order SC-FC coupling between MDD patients and HCs. Our study confirmed that higher-order coupling between SC and FC plays an important role in diagnosing MDD. The identification of new biological evidence for MDD etiology holds promise for the development of innovative antidepressant therapies.

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

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

Can structure predict function at individual level in the human connectome?

Several studies predicting Functional Connectivity (FC) from Structural Connectivity (SC) at individual level have been published in recent years, each promising increased performance and utility. We investigated three of these studies, analyzing whether the results truly represent a meaningful individual-level mapping from SC to FC. Using data from the Human Connectome Project shared accross the three studies, we constructed a predictor by averaging FC of training data and analyzed its performance in the same way. In each case, we found that group average FC is an equivalent or better predictor of individual FC than the predictive models in terms of raw prediction performance. Furthermore, we showed that additional analyses performed by the authors of the three studies, in which they attempt to show that their predicted FC has value beyond raw prediction performance, could also be reproduced using the group average FC predictor. This makes it unclear whether any of the three methods represent a meaningful individual-level predictive model. We conclude that either the methods are not appropriate for the data, that the sample size is too small, or that the data does not contain sufficient information to learn a mapping from SC to FC. We advise future individual-level studies to explicitly report results in comparison to the performance of the group average, and carefully demonstrate that their predictions contain meaningful individual-level information. Finally, we believe that investigating alternatives for the construction of SC and FC may improve the chances of developing a meaningful individual-level mapping from SC to FC.

Diffusion time-related structure-function coupling reveals differential association with inter-individual variations in body mass index

Body mass index (BMI) is an indicator of obesity, and recent neuroimaging studies have demonstrated that inter-individual variations in BMI are associated with altered brain structure and function. However, the mechanism underlying the alteration of structure-function correspondence according to BMI is under-investigated. In this study, we studied structural and functional connectivity derived from diffusion MRI tractography and inter-regional correlations of functional MRI time series, respectively. We combined the structural and functional connectivity information using the Riemannian optimization approach. First, the low-dimensional principal eigenvectors (i.e., gradients) of the structural connectivity were generated by applying diffusion map embedding with varying diffusion times. A transformation was identified so that the structural and functional embeddings share the same coordinate system, and subsequently, the functional connectivity matrix was simulated. Then, we generated gradients from the simulated functional connectivity matrix. We found the most apparent cortical hierarchical organization differentiating between low-level sensory and higher-order transmodal regions in the middle of the diffusion time, indicating that the hierarchical organization of the brain may reflect the intermediate mechanisms of mono- and polysynaptic communications. Associations between the functional gradients and BMI were strongest when the hierarchical structure was the most evident. Moreover, the gradient-BMI association map was related to the microstructural features, and the findings indicated that the BMI-related structure-function coupling was significantly associated with brain microstructure, particularly in higher-order transmodal areas. Finally, transcriptomic association analysis revealed the potential biological underpinnings specifying gene enrichment in the striatum, hypothalamus, and cortical cells. Our findings provide evidence that structure-function correspondence is strongly coupled with BMI when hierarchical organization is the most apparent and that the associations are related to the multiscale properties of the brain, leading to an advanced understanding of the neural mechanisms related to BMI.

Predicting resting-state brain functional connectivity from the structural connectome using the heat diffusion model: a multiple-timescale fusion method

Objective.Understanding the intricate relationship between structural connectivity (SC) and functional connectivity (FC) is pivotal for understanding the complexities of the human brain. To explore this relationship, the heat diffusion model (HDM) was utilized to predict FC from SC. However, previous studies using the HDM have typically predicted FC at a critical time scale in the heat kernel equation, overlooking the dynamic nature of the diffusion process and providing an incomplete representation of the predicted FC.Approach.In this study, we propose an alternative approach based on the HDM. First, we introduced a multiple-timescale fusion method to capture the dynamic features of the diffusion process. Additionally, to enhance the smoothness of the predicted FC values, we employed the Wavelet reconstruction method to maintain local consistency and remove noise. Moreover, to provide a more accurate representation of the relationship between SC and FC, we calculated the linear transformation between the smoothed FC and the empirical FC.Main results.We conducted extensive experiments in two independent datasets. By fusing different time scales in the diffusion process for predicting FC, the proposed method demonstrated higher predictive correlation compared with method considering only critical time points (Singlescale). Furthermore, compared with other existing methods, the proposed method achieved the highest predictive correlations of 0.6939±0.0079 and 0.7302±0.0117 on the two datasets respectively. We observed that the visual network at the network level and the parietal lobe at the lobe level exhibited the highest predictive correlations, indicating that the functional activity in these regions may be closely related to the direct diffusion of information between brain regions.Significance.The multiple-timescale fusion method proposed in this study provides insights into the dynamic aspects of the diffusion process, contributing to a deeper understanding of how brain structure gives rise to brain function.

Connectome-wide structure-function coupling models implicate polysynaptic alterations in autism

Autism spectrum disorder (ASD) is one of the most common neurodevelopmental diagnoses. Although incompletely understood, structural and functional network alterations are increasingly recognized to be at the core of the condition. We utilized multimodal imaging and connectivity modeling to study structure-function coupling in ASD and probed mono- and polysynaptic mechanisms on structurally-governed network function. We examined multimodal magnetic resonance imaging data in 80 ASD and 61 neurotypical controls from the Autism Brain Imaging Data Exchange (ABIDE) II initiative. We predicted intrinsic functional connectivity from structural connectivity data in each participant using a Riemannian optimization procedure that varies the times that simulated signals can unfold along tractography-derived personalized connectomes. In both ASD and neurotypical controls, we observed improved structure-function prediction at longer diffusion time scales, indicating better modeling of brain function when polysynaptic mechanisms are accounted for. Prediction accuracy differences (∆prediction accuracy) were marked in transmodal association systems, such as the default mode network, in both neurotypical controls and ASD. Differences were, however, lower in ASD in a polysynaptic regime at higher simulated diffusion times. We compared regional differences in ∆prediction accuracy between both groups to assess the impact of polysynaptic communication on structure-function coupling. This analysis revealed that between-group differences in ∆prediction accuracy followed a sensory-to-transmodal cortical hierarchy, with an increased gap between controls and ASD in transmodal compared to sensory/motor systems. Multivariate associative techniques revealed that structure-function differences reflected inter-individual differences in autistic symptoms and verbal as well as non-verbal intelligence. Our network modeling approach sheds light on atypical structure-function coupling in autism, and suggests that polysynaptic network mechanisms are implicated in the condition and that these can help explain its wide range of associated symptoms.

Eigenmode-based approach reveals a decline in brain structure-function liberality across the human lifespan

While brain function is supported and constrained by the underlying structure, the connectome-based link estimated by current approaches is either relatively moderate or accompanied by high model complexity, with the essential principles underlying structure-function coupling remaining elusive. Here, by proposing a mapping method based on network eigendecomposition, we present a concise and strong correspondence between structure and function. We show that the explanation of functional connectivity can be significantly improved by incorporating interactions between different structural eigenmodes. We also demonstrate the pronounced advantage of the present mapping in capturing individual-specific information with simple implementation. Applying our methodology to the human lifespan, we find that functional diversity decreases with age, with functional interactions increasingly dominated by the leading functional mode. We also find that structure-function liberality weakens with age, which is driven by the decreases in functional components that are less constrained by anatomy, while the magnitude of structure-aligned components is preserved. Overall, our work enhances the understanding of structure-function coupling from a collective, connectome-oriented perspective and promotes a more refined identification of functional portions relevant to human aging, holding great potential for mechanistic insights into individual differences associated with cognition, development, and neurological disorders.

BrainStat: A toolbox for brain-wide statistics and multimodal feature associations

Analysis and interpretation of neuroimaging datasets has become a multidisciplinary endeavor, relying not only on statistical methods, but increasingly on associations with respect to other brain-derived features such as gene expression, histological data, and functional as well as cognitive architectures. Here, we introduce BrainStat - a toolbox for (i) univariate and multivariate linear models in volumetric and surface-based brain imaging datasets, and (ii) multidomain feature association of results with respect to spatial maps of post-mortem gene expression and histology, task-based fMRI meta-analysis, as well as resting-state fMRI motifs across several common surface templates. The combination of statistics and feature associations into a turnkey toolbox streamlines analytical processes and accelerates cross-modal research. The toolbox is implemented in both Python and MATLAB, two widely used programming languages in the neuroimaging and neuroinformatics communities. BrainStat is openly available and complemented by an expandable documentation.

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

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