Structural mediation of human brain activity revealed by white-matter interpolation of fMRI

Functional Connectivity Prediction With Deep Learning for Graph Transformation

Inferring resting-state functional connectivity (FC) from anatomical brain wiring, known as structural connectivity (SC), is of enormous significance in neuroscience for understanding biological neuronal networks and treating mental diseases. Both SC and FC are networks where the nodes are brain regions, and in SC, the edges are the physical fiber nerves among the nodes, while in FC, the edges are the nodes' coactivation relations. Despite the importance of SC and FC, until very recently, the rapidly growing research body on this topic has generally focused on either linear models or computational models that rely heavily on heuristics and simple assumptions regarding the mapping between FC and SC. However, the relationship between FC and SC is actually highly nonlinear and complex and contains considerable randomness; additional factors, such as the subject's age and health, can also significantly impact the SC-FC relationship and hence cannot be ignored. To address these challenges, here, we develop a novel SC-to-FC generative adversarial network (SF-GAN) framework for mapping SC to FC, along with additional metafeatures based on a newly proposed graph neural network-based generative model that is capable of learning the stochasticity. Specifically, a new graph-based conditional generative adversarial nets model is proposed, where edge convolution layers are leveraged to encode the graph patterns in the SC in the form of a graph representation. New edge deconvolution layers are then utilized to decode the representation back to FC. Additional metafeatures of subjects' profile information are integrated into the graph representation with newly designed sparse-regularized layers that can automatically select features that impact FC. Finally, we have also proposed new post hoc explainer of our SF-GAN, which can identify which subgraphs in SC strongly influence which subgraphs in FC by a new multilevel edge-correlation-guided graph clustering problem. The results of experiments conducted to test the new model confirm that it significantly outperforms existing state-of-the-art methods, with additional interpretability for identifying important metafeatures and subgraphs.

Improved Functionnectome by dissociating the contributions of white matter fiber classes to functional activation

Integrating the underlying brain circuit's structural and functional architecture is required to explore the functional organization of cognitive networks. In that regard, we recently introduced the Functionnectome. This structural-functional method combines an fMRI acquisition with tractography-derived white matter connectivity data to map cognitive processes onto the white matter. However, this multimodal integration faces three significant challenges: (1) the necessarily limited overlap between tractography streamlines and the grey matter, which may reduce the amount of functional signal associated with the related structural connectivity; (2) the scrambling effect of crossing fibers on functional signal, as a single voxel in such regions can be structurally connected to several cognitive networks with heterogeneous functional signals; and (3) the difficulty of interpretation of the resulting cognitive maps, as crossing and overlapping white matter tracts can obscure the organization of the studied network. In the present study, we tackled these problems by developing a streamline-extension procedure and dividing the white matter anatomical priors between association, commissural, and projection fibers. This approach significantly improved the characterization of the white matter involvement in the studied cognitive processes. The new Functionnectome priors produced are now readily available, and the analysis workflow highlighted here should also be generalizable to other structural-functional approaches. We improved the Functionnectome approach to better study the involvement of white matter in brain function by separating the analysis of the three classes of white matter fibers (association, commissural, and projection fibers). This step successfully clarified the activation maps and increased their statistical significance.

Estimating dynamic individual coactivation patterns based on densely sampled resting-state fMRI data and utilizing it for better subject identification

As a complex dynamic system, the brain exhibits spatially organized recurring patterns of activity over time. Coactivation patterns (CAPs), which analyzes data from each single frame, have been utilized to detect transient brain activity states recently. However, previous CAP analyses have been conducted at the group level, which might neglect meaningful individual differences. Here, we estimated individual CAP states at both subject- and scan-level based on a densely sampled dataset: Midnight Scan Club. We used differential identifiability, which measures the gap between intra- and inter-subject similarity, to evaluate individual differences. We found individual CAPs at the subject-level achieved the best fingerprinting ability by maintaining high intra-subject similarity and enlarging inter-subject differences, and brain regions of association networks mainly contributed to the identifiability. On the other hand, scan-level CAP states were unstable across scans for the same participant. Expectedly, we found subject-specific CAPs became more reliable and discriminative with more data (i.e., longer duration). As the acquisition time of each participant is limited in practice, our results recommend a data collection strategy that collects more scans with appropriate duration (e.g., 12 ~ 15 min/scan) to obtain more reliable subject-specific CAPs, when total acquisition time is fixed (e.g., 150 min). In summary, this work has constructed reliable subject-specific CAP states with meaningful individual differences, and recommended an appropriate data collection strategy, which can guide subsequent investigations into individualized brain dynamics.

Atlasing white matter and grey matter joint contributions to resting-state networks in the human brain

Over the past two decades, the study of resting-state functional magnetic resonance imaging has revealed that functional connectivity within and between networks is linked to cognitive states and pathologies. However, the white matter connections supporting this connectivity remain only partially described. We developed a method to jointly map the white and grey matter contributing to each resting-state network (RSN). Using the Human Connectome Project, we generated an atlas of 30 RSNs. The method also highlighted the overlap between networks, which revealed that most of the brain's white matter (89%) is shared between multiple RSNs, with 16% shared by at least 7 RSNs. These overlaps, especially the existence of regions shared by numerous networks, suggest that white matter lesions in these areas might strongly impact the communication within networks. We provide an atlas and an open-source software to explore the joint contribution of white and grey matter to RSNs and facilitate the study of the impact of white matter damage to these networks. In a first application of the software with clinical data, we were able to link stroke patients and impacted RSNs, showing that their symptoms aligned well with the estimated functions of the networks.

Joint connectivity matrix independent component analysis: Auto-linking of structural and functional connectivities

The study of human brain connectivity, including structural connectivity (SC) and functional connectivity (FC), provides insights into the neurophysiological mechanism of brain function and its relationship to human behavior and cognition. Both types of connectivity measurements provide crucial yet complementary information. However, integrating these two modalities into a single framework remains a challenge, because of the differences in their quantitative interdependencies as well as their anatomical representations due to distinctive imaging mechanisms. In this study, we introduced a new method, joint connectivity matrix independent component analysis (cmICA), which provides a data-driven parcellation and automated-linking of SC and FC information simultaneously using a joint analysis of functional magnetic resonance imaging (MRI) and diffusion-weighted MRI data. We showed that these two connectivity modalities produce common cortical segregation, though with various degrees of (dis)similarity. Moreover, we show conjoint FC networks and structural white matter tracts that directly link these cortical parcellations/sources, within one analysis. Overall, data-driven joint cmICA provides a new approach for integrating or fusing structural connectivity and FC systematically and conveniently, and provides an effective tool for connectivity-based multimodal data fusion in brain.

Frequency-specific coactivation patterns in resting-state and their alterations in schizophrenia: An fMRI study

The resting‐state human brain is a dynamic system that shows frequency‐dependent characteristics. Recent studies demonstrate that coactivation pattern (CAP) analysis can identify recurring brain states with similar coactivation configurations. However, it is unclear whether and how CAPs depend on the frequency bands. The current study investigated the spatial and temporal characteristics of CAPs in the four frequency sub‐bands from slow‐5 (0.01–0.027 Hz), slow‐4 (0.027–0.073 Hz), slow‐3 (0.073–0.198 Hz), to slow‐2 (0.198–0.25 Hz), in addition to the typical low‐frequency range (0.01–0.08 Hz). In the healthy subjects, six CAP states were obtained at each frequency band in line with our prior study. Similar spatial patterns with the typical range were observed in slow‐5, 4, and 3, but not in slow‐2. While the frequency increased, all CAP states displayed shorter persistence, which caused more between‐state transitions. Specifically, from slow‐5 to slow‐4, the coactivation not only changed significantly in distributed cortical networks, but also increased in the basal ganglia as well as the amygdala. Schizophrenia patients showed significant alteration in the persistence of CAPs of slow‐5. Using leave‐one‐pair‐out, hold‐out and resampling validations, the highest classification accuracy (84%) was achieved by slow‐4 among different frequency bands. In conclusion, our findings provide novel information about spatial and temporal characteristics of CAP states at different frequency bands, which contributes to a better understanding of the frequency aspect of biomarkers for schizophrenia and other disorders.

Post-stroke reorganization of transient brain activity characterizes deficits and recovery of cognitive functions

Functional magnetic resonance imaging (fMRI) has been widely employed to study stroke pathophysiology. In particular, analyses of fMRI signals at rest were directed at quantifying the impact of stroke on spatial features of brain networks. However, brain networks have intrinsic time features that were, so far, disregarded in these analyses. In consequence, standard fMRI analysis failed to capture temporal imbalance resulting from stroke lesions, hence restricting their ability to reveal the interdependent pathological changes in structural and temporal network features following stroke. Here, we longitudinally analyzed hemodynamic-informed transient activity in a large cohort of stroke patients (n = 103) to assess spatial and temporal changes of brain networks after stroke. Metrics extracted from the hemodynamic-informed transient activity were replicable within- and between-individuals in healthy participants, hence supporting their robustness and their clinical applicability. While large-scale spatial patterns of brain networks were preserved after stroke, their durations were altered, with stroke subjects exhibiting a varied pattern of longer and shorter network activations compared to healthy individuals. Specifically, patients showed a longer duration in the lateral precentral gyrus and anterior cingulum, and a shorter duration in the occipital lobe and in the cerebellum. These temporal alterations were associated with white matter damage in projection and association pathways. Furthermore, they were tied to deficits in specific behavioral domains as restoration of healthy brain dynamics paralleled recovery of cognitive functions (attention, language and spatial memory), but was not significantly correlated to motor recovery. These findings underscore the critical importance of network temporal properties in dissecting the pathophysiology of brain changes after stroke, thus shedding new light on the clinical potential of time-resolved methods for fMRI analysis.

Dynamics of Functional Network Organization Through Graph Mixture Learning

Understanding the organizational principles of human brain activity at the systems level remains a major challenge in network neuroscience. Here, we introduce a fully data-driven approach based on graph learning to extract meaningful repeating network patterns from regionally-averaged timecourses. We use the Graph Laplacian Mixture Model (GLMM), a generative model that treats functional data as a collection of signals expressed on multiple underlying graphs. By exploiting covariance between activity of brain regions, these graphs can be learned without resorting to structural information. To validate the proposed technique, we first apply it to task fMRI with a known experimental paradigm. The probability of each graph to occur at each time-point is found to be consistent with the task timing, while the spatial patterns associated to each epoch of the task are in line with previously established activation patterns using classical regression analysis. We further on apply the technique to resting state data, which leads to extracted graphs that correspond to well-known brain functional activation patterns. The GLMM allows to learn graphs entirely from the functional activity that, in practice, turn out to reveal high degrees of similarity to the structural connectome. The Default Mode Network (DMN) is always captured by the algorithm in the different tasks and resting state data. Therefore, we compare the states corresponding to this network within themselves and with structure. Overall, this method allows us to infer relevant functional brain networks without the need of structural connectome information. Moreover, we overcome the limitations of windowing the time sequences by feeding the GLMM with the whole functional signal and neglecting the focus on sub-portions of the signals.

Functionnectome as a framework to analyse the contribution of brain circuits to fMRI

In recent years, the field of functional neuroimaging has moved away from a pure localisationist approach of isolated functional brain regions to a more integrated view of these regions within functional networks. However, the methods used to investigate functional networks rely on local signals in grey matter and are limited in identifying anatomical circuitries supporting the interaction between brain regions. Mapping the brain circuits mediating the functional signal between brain regions would propel our understanding of the brain's functional signatures and dysfunctions. We developed a method to unravel the relationship between brain circuits and functions: The Functionnectome. The Functionnectome combines the functional signal from fMRI with white matter circuits' anatomy to unlock and chart the first maps of functional white matter. To showcase this method's versatility, we provide the first functional white matter maps revealing the joint contribution of connected areas to motor, working memory, and language functions. The Functionnectome comes with an open-source companion software and opens new avenues into studying functional networks by applying the method to already existing datasets and beyond task fMRI.

Large-scale functional network dynamics in human callosal agenesis: Increased subcortical involvement and preserved laterality

In the human brain, the corpus callosum is the major white-matter commissural tract enabling the transmission of sensory-motor, and higher level cognitive information between homotopic regions of the two cerebral hemispheres. Despite developmental absence (i.e., agenesis) of the corpus callosum (AgCC), functional connectivity is preserved, including interhemispheric connectivity. Subcortical structures have been hypothesised to provide alternative pathways to enable this preservation. To test this hypothesis, we used functional Magnetic Resonance Imaging (fMRI) recordings in children with AgCC and typically developing children, and a time-resolved approach to retrieve temporal characteristics of whole-brain functional networks. We observed an increased engagement of the cerebellum and amygdala/hippocampus networks in children with AgCC compared to typically developing children. There was little evidence that laterality of activation networks was affected in AgCC. Our findings support the hypothesis that subcortical structures play an essential role in the functional reconfiguration of the brain in the absence of a corpus callosum.

Diffusion-informed spatial smoothing of fMRI data in white matter using spectral graph filters

Brain activation mapping using functional magnetic resonance imaging (fMRI) has been extensively studied in brain gray matter (GM), whereas in large disregarded for probing white matter (WM). This unbalanced treatment has been in part due to controversies in relation to the nature of the blood oxygenation level-dependent (BOLD) contrast in WM and its detectability. However, an accumulating body of studies has provided solid evidence of the functional significance of the BOLD signal in WM and has revealed that it exhibits anisotropic spatio-temporal correlations and structure-specific fluctuations concomitant with those of the cortical BOLD signal. In this work, we present an anisotropic spatial filtering scheme for smoothing fMRI data in WM that accounts for known spatial constraints on the BOLD signal in WM. In particular, the spatial correlation structure of the BOLD signal in WM is highly anisotropic and closely linked to local axonal structure in terms of shape and orientation, suggesting that isotropic Gaussian filters conventionally used for smoothing fMRI data are inadequate for denoising the BOLD signal in WM. The fundamental element in the proposed method is a graph-based description of WM that encodes the underlying anisotropy observed across WM, derived from diffusion-weighted MRI data. Based on this representation, and leveraging graph signal processing principles, we design subject-specific spatial filters that adapt to a subject's unique WM structure at each position in the WM that they are applied at. We use the proposed filters to spatially smooth fMRI data in WM, as an alternative to the conventional practice of using isotropic Gaussian filters. We test the proposed filtering approach on two sets of simulated phantoms, showcasing its greater sensitivity and specificity for the detection of slender anisotropic activations, compared to that achieved with isotropic Gaussian filters. We also present WM activation mapping results on the Human Connectome Project's 100-unrelated subject dataset, across seven functional tasks, showing that the proposed method enables the detection of streamline-like activations within axonal bundles.

CHARACTERIZATION OF SPATIAL DYNAMICS OF FMRI DATA IN WHITE MATTER USING DIFFUSION-INFORMED WHITE MATTER HARMONICS

In this work, we leverage the Laplacian eigenbasis of voxel-wise white matter (WM) graphs derived from diffusion-weighted MRI data, dubbed WM harmonics, to characterize the spatial structure of WM fMRI data. Our motivation for such a characterization is based on studies that show WM fMRI data exhibit a spatial correlational anisotropy that coincides with underlying fiber patterns. By quantifying the energy content of WM fMRI data associated with subsets of WM harmonics across multiple spectral bands, we show that the data exhibits notable subtle spatial modulations under functional load that are not manifested during rest. WM harmonics provide a novel means to study the spatial dynamics of WM fMRI data, in such way that the analysis is informed by the underlying anatomical structure.

Dynamic functional networks in idiopathic normal pressure hydrocephalus: Alterations and reversibility by CSF tap test

Idiopathic Normal Pressure Hydrocephalus (iNPH)—the leading cause of reversible dementia in aging—is characterized by ventriculomegaly and gait, cognitive and urinary impairments. Despite its high prevalence estimated at 6% among the elderlies, iNPH remains underdiagnosed and undertreated due to the lack of iNPH‐specific diagnostic markers and limited understanding of pathophysiological mechanisms. INPH diagnosis is also complicated by the frequent occurrence of comorbidities, the most common one being Alzheimer's disease (AD). Here we investigate the resting‐state functional magnetic resonance imaging dynamics of 26 iNPH patients before and after a CSF tap test, and of 48 normal older adults. Alzheimer's pathology was evaluated by CSF biomarkers. We show that the interactions between the default mode, and the executive‐control, salience and attention networks are impaired in iNPH, explain gait and executive disturbances in patients, and are not driven by AD‐pathology. In particular, AD molecular biomarkers are associated with functional changes distinct from iNPH functional alterations. Finally, we demonstrate a partial normalization of brain dynamics 24 hr after a CSF tap test, indicating functional plasticity mechanisms. We conclude that functional changes involving the default mode cross‐network interactions reflect iNPH pathophysiological mechanisms and track treatment response, possibly contributing to iNPH differential diagnosis and better clinical management.

Regionconnect: Rapidly extracting standardized brain connectivity information in voxel-wise neuroimaging studies

Reporting white matter findings in voxel-wise neuroimaging studies typically lacks specificity in terms of brain connectivity. Therefore, the purpose of this work was to develop an approach for rapidly extracting standardized brain connectivity information for white matter regions with significant findings in voxel-wise neuroimaging studies. The new approach was named regionconnect and is based on precalculated average healthy adult brain connectivity information stored in standard space in a fashion that allows fast retrieval and integration. Towards this goal, the present work first generated and evaluated the white matter connectome of the IIT Human Brain Atlas v.5.0. It was demonstrated that the edges of the atlas connectome are representative of those of individual participants of the Human Connectome Project in terms of the spatial organization of streamlines and spatial patterns of track-density. Next, the new white matter connectome was used to develop multi-layer, connectivity-based labels for each white matter voxel of the atlas, consistent with the fact that each voxel may contain axons from multiple connections. The regionconnect algorithm was then developed to rapidly integrate information contained in the multi-layer labels across voxels of a white matter region and to generate a list of the most probable connections traversing that region. Usage of regionconnect does not require high angular resolution diffusion MRI or any MRI data. The regionconnect algorithm as well as the white matter tractogram and connectome, multi-layer, connectivity-based labels, and associated resources developed for the IIT Human Brain Atlas v.5.0 in this work are available at www.nitrc.org/projects/iit. An interactive, online version of regionconnect is also available at www.iit.edu/~mri.

Hemodynamic Response Function in Brain White Matter in a Resting State

The hemodynamic response function (HRF) characterizes temporal variations of blood oxygenation level-dependent (BOLD) signals. Although a variety of HRF models have been proposed for gray matter responses to functional demands, few studies have investigated HRF profiles in white matter particularly under resting conditions. In the present work we quantified the nature of the HRFs that are embedded in resting state BOLD signals in white matter, and which modulate the temporal fluctuations of baseline signals. We demonstrate that resting state HRFs in white matter could be derived by referencing to intrinsic avalanches in gray matter activities, and the derived white matter HRFs had reduced peak amplitudes and delayed peak times as compared with those in gray matter. Distributions of the time delays and correlation profiles in white matter depend on gray matter activities as well as white matter tract distributions, indicating that resting state BOLD signals in white matter encode neural activities associated with those of gray matter. This is the first investigation of derivations and characterizations of resting state HRFs in white matter and their relations to gray matter activities. Findings from this work have important implications for analysis of BOLD signals in the brain.