Brain parcellation driven by dynamic functional connectivity better capture intrinsic network dynamics

Static and dynamic brain functional connectivity patterns in patients with unilateral moderate-to-severe asymptomatic carotid stenosis

Background and purpose

Asymptomatic carotid stenosis (ACS) is an independent risk factor for ischemic stroke and vascular cognitive impairment, affecting cognitive function across multiple domains. This study aimed to explore differences in static and dynamic intrinsic functional connectivity and temporal dynamics between patients with ACS and those without carotid stenosis.

Methods

We recruited 30 patients with unilateral moderate-to-severe (stenosis ≥ 50%) ACS and 30 demographically-matched healthy controls. All participants underwent neuropsychological testing and 3.0T brain MRI scans. Resting-state functional MRI (rs-fMRI) was used to calculate both static and dynamic functional connectivity. Dynamic independent component analysis (dICA) was employed to extract independent circuits/networks and to detect time-frequency modulation at the circuit level. Further imaging-behavior associations identified static and dynamic functional connectivity patterns that reflect cognitive decline.

Results

ACS patients showed altered functional connectivity in multiple brain regions and networks compared to controls. Increased connectivity was observed in the inferior parietal lobule, frontal lobe, and temporal lobe. dICA further revealed changes in the temporal frequency of connectivity in the salience network. Significant differences in the temporal variability of connectivity were found in the fronto-parietal network, dorsal attention network, sensory-motor network, language network, and visual network. The temporal parameters of these brain networks were also related to overall cognition and memory.

Conclusions

These results suggest that ACS involves not only changes in the static large-scale brain network connectivity but also dynamic temporal variations, which parallel overall cognition and memory recall.

Functional anatomy and topographical organization of the frontotemporal arcuate fasciculus

Traditionally, the frontotemporal arcuate fasciculus (AF) is viewed as a single entity in anatomo-clinical models. However, it is unclear if distinct cortical origin and termination patterns within this bundle correspond to specific language functions. We use track-weighted dynamic functional connectivity, a hybrid imaging technique, to study the AF structure and function in two distinct datasets of healthy subjects. Here we show that the AF can be subdivided based on dynamic changes in functional connectivity at the streamline endpoints. An unsupervised parcellation algorithm reveals spatially segregated subunits, which are then functionally quantified through meta-analysis. This approach identifies three distinct clusters within the AF - ventral, middle, and dorsal frontotemporal AF - each linked to different frontal and temporal termination regions and likely involved in various language production and comprehension aspects. Our findings may have relevant implications for the understanding of the functional anatomy of the AF as well as its contribution to linguistic and non-linguistic functions.

Track-Weighted Dynamic Functional Connectivity Profiles and Topographic Organization of the Human Pulvinar

The human pulvinar is considered a prototypical associative thalamic nucleus as it represents a key node in several cortico-subcortical networks. Through this extensive connectivity to widespread brain areas, it has been suggested that the pulvinar may play a central role in modulating cortical oscillatory dynamics of complex cognitive and executive functions. Additionally, derangements of pulvinar activity are involved in different neuropsychiatric conditions including Lewy-body disease, Alzheimer's disease, and schizophrenia. Anatomical investigations in nonhuman primates have demonstrated a topographical organization of cortico-pulvinar connectivity along its dorsoventral and rostrocaudal axes; this specific organization shows only partial overlap with the traditional subdivision into subnuclei (anterior, lateral, medial, and inferior) and is thought to coordinate information processing within specific brain networks. However, despite its relevance in mediating higher-order cognitive functions, such a structural and functional organization of the pulvinar in the human brain remains poorly understood. Track-weighted dynamic functional connectivity (tw-dFC) is a recently developed technique that combines structural and dynamic functional connectivity, allowing the identification of white matter pathways underlying the fluctuations observed in functional connectivity between brain regions over time. Herein, we applied a data-driven parcellation approach to reveal topographically organized connectivity clusters within the human pulvinar complex, in two large cohorts of healthy human subjects. Unsupervised clustering of tw-dFC time series within the pulvinar complex revealed dorsomedial, dorsolateral, ventral anterior, and ventral posterior connectivity clusters. Each of these clusters shows functional coupling to specific, widespread cortico-subcortical white matter brain networks. Altogether, our findings represent a relevant step towards a better understanding of pulvinar anatomy and function, and a detailed characterization of his role in healthy and pathological conditions.

Prediction of individual performance and verbal intelligence scores from resting-state fMRI in children and adolescents

The neuroimaging basis of intelligence remains elusive; however, there is a growing body of research employing connectome-based predictive modeling to estimate individual intelligence scores, aiming to identify the optimal set of neuroimaging features for accurately predicting an individual's cognitive abilities. Compared to adults, the disparities in cognitive performance among children and adolescents are more likely to captivate individuals' interest and attention. Limited research has been dedicated to exploring neuroimaging markers of intelligence specifically in the pediatric population. In this study, we utilized resting-state functional magnetic resonance imaging (fMRI) and intelligence quotient (IQ) scores of 170 healthy children and adolescents obtained from a public database to identify brain functional connectivity markers associated with individual intellectual behavior. Initially, we extracted and summarized relevant resting-state features from whole-brain or functional network connectivity that were most pertinent to IQ scores. Subsequently, these features were employed to establish prediction models for both performance and verbal IQ scores. Within a 10-fold cross-validation framework, our findings revealed that prediction models based on whole-brain functional connectivity effectively predicted performance IQ scores( R = 0.35 , P = 2.2 × 10 - 4 ) but not verbal IQ scores( R = 0.12 , P = 0.20 ). Results of prediction models based on brain functional network connectivity further demonstrated the exceptional predictive ability of the default mode network (DMN) and fronto-parietal task control network (FTPN) for performance IQ scores ( R = 0.71 , P = 2.2 × 10 - 18 ). The above findings have also been validated using an independent dataset. Our findings suggest that the performance IQ of children and adolescents primarily relies on the connectivity of brain regions associated with DMN and FTPN. Moreover, variations in intellectual performance during childhood and adolescences are closely linked to alterations in brain functional network connectivity.

Graph Learning for Cortical Parcellation from Tensor Decompositions of Resting-State fMRI

Cortical parcellation has long been a cornerstone in the field of neuroscience, enabling the cerebral cortex to be partitioned into distinct, non-overlapping regions that facilitate the interpretation and comparison of complex neuroscientific data. In recent years, these parcellations have frequently been based on the use of resting-state fMRI (rsfMRI) data. In parallel, methods such as independent components analysis have long been used to identify large-scale functional networks with significant spatial overlap between networks. Despite the fact that both forms of decomposition make use of the same spontaneous brain activity measured with rsfMRI, a gap persists in establishing a clear relationship between disjoint cortical parcellations and brain-wide networks. To address this, we introduce a novel parcellation framework that integrates NASCAR, a three-dimensional tensor decomposition method that identifies a series of functional brain networks, with state-of-the-art graph representation learning to produce cortical parcellations that represent near-homogeneous functional regions that are consistent with these brain networks. Further, through the use of the tensor decomposition, we avoid the limitations of traditional approaches that assume statistical independence or orthogonality in defining the underlying networks. Our findings demonstrate that these parcellations are comparable or superior to established atlases in terms of homogeneity of the functional connectivity across parcels, task contrast alignment, and architectonic map alignment. Our methodological pipeline is highly automated, allowing for rapid adaptation to new datasets and the generation of custom parcellations in just minutes, a significant advancement over methods that require extensive manual input. We describe this integrated approach, which we refer to as Untamed, as a tool for use in the fields of cognitive and clinical neuroscientific research. Parcellations created from the Human Connectome Project dataset using Untamed, along with the code to generate atlases with custom parcel numbers, are publicly available at https://untamed-atlas.github.io.

Identifying canonical and replicable multi-scale intrinsic connectivity networks in 100k+ resting-state fMRI datasets

Despite the known benefits of data-driven approaches, the lack of approaches for identifying functional neuroimaging patterns that capture both individual variations and inter-subject correspondence limits the clinical utility of rsfMRI and its application to single-subject analyses. Here, using rsfMRI data from over 100k individuals across private and public datasets, we identify replicable multi-spatial-scale canonical intrinsic connectivity network (ICN) templates via the use of multi-model-order independent component analysis (ICA). We also study the feasibility of estimating subject-specific ICNs via spatially constrained ICA. The results show that the subject-level ICN estimations vary as a function of the ICN itself, the data length, and the spatial resolution. In general, large-scale ICNs require less data to achieve specific levels of (within- and between-subject) spatial similarity with their templates. Importantly, increasing data length can reduce an ICN's subject-level specificity, suggesting longer scans may not always be desirable. We also find a positive linear relationship between data length and spatial smoothness (possibly due to averaging over intrinsic dynamics), suggesting studies examining optimized data length should consider spatial smoothness. Finally, consistency in spatial similarity between ICNs estimated using the full data and subsets across different data lengths suggests lower within-subject spatial similarity in shorter data is not wholly defined by lower reliability in ICN estimates, but may be an indication of meaningful brain dynamics which average out as data length increases.

Enhancing precision in human neuroscience

Human neuroscience has always been pushing the boundary of what is measurable. During the last decade, concerns about statistical power and replicability - in science in general, but also specifically in human neuroscience - have fueled an extensive debate. One important insight from this discourse is the need for larger samples, which naturally increases statistical power. An alternative is to increase the precision of measurements, which is the focus of this review. This option is often overlooked, even though statistical power benefits from increasing precision as much as from increasing sample size. Nonetheless, precision has always been at the heart of good scientific practice in human neuroscience, with researchers relying on lab traditions or rules of thumb to ensure sufficient precision for their studies. In this review, we encourage a more systematic approach to precision. We start by introducing measurement precision and its importance for well-powered studies in human neuroscience. Then, determinants for precision in a range of neuroscientific methods (MRI, M/EEG, EDA, Eye-Tracking, and Endocrinology) are elaborated. We end by discussing how a more systematic evaluation of precision and the application of respective insights can lead to an increase in reproducibility in human neuroscience.

Considering dynamic nature of the brain: the clinical importance of connectivity variability in machine learning classification and prediction

The brain connectivity of resting-state fMRI (rs-fMRI) represents an intrinsic state of brain architecture, and it has been used as a useful neural marker for detecting psychiatric conditions as well as for predicting psychosocial characteristics. However, most studies using brain connectivity have focused more on the strength of functional connectivity over time (static-FC) but less attention to temporal characteristics of connectivity changes (FC-variability). The primary goal of the current study was to investigate the effectiveness of using the FC-variability in classifying an individual's pathological characteristics from others and predicting psychosocial characteristics. In addition, the current study aimed to prove that benefits of the FC-variability are reliable across various analysis procedures. To this end, three open public large resting-state fMRI datasets including individuals with Autism Spectrum Disorder (ABIDE; N = 1249), Schizophrenia disorder (COBRE; N = 145), and typical development (NKI; N = 672) were utilized for the machine learning (ML) classification and prediction based on their static-FC and the FC-variability metrics. To confirm the robustness of FC-variability utility, we benchmarked the ML classification and prediction with various brain parcellations and sliding window parameters. As a result, we found that the ML performances were significantly improved when the ML included FC-variability features in classifying pathological populations from controls (e.g., individuals with autism spectrum disorder vs. typical development) and predicting psychiatric severity (e.g., score of autism diagnostic observation schedule), regardless of parcellation selection and sliding window size. Additionally, the ML performance deterioration was significantly prevented with FC-variability features when excessive features were inputted into the ML models, yielding more reliable results. In conclusion, the current finding proved the usefulness of the FC-variability and its reliability.

Neural activity changes in first-episode, drug-naïve patients with major depressive disorder after transcutaneous auricular vagus nerve stimulation treatment: A resting-state fMRI study

Introduction

Major depressive disorder (MDD) is a disease with prominent individual, medical, and economic impacts. Drug therapy and other treatment methods (such as Electroconvulsive therapy) may induce treatment-resistance and have associated side effects including loss of memory, decrease of reaction time, and residual symptoms. Transcutaneous auricular vagus nerve stimulation (taVNS) is a novel and non-invasive treatment approach which stimulates brain structures with no side-effects. However, it remains little understood whether and how the neural activation is modulated by taVNS in MDD patients. Herein, we used the regional homogeneity (ReHo) to investigate the brain activity in first-episode, drug-naïve MDD patients after taVNS treatment.

Materials and methods

Twenty-two first-episode, drug-naïve MDD patients were enrolled in the study. These patients received the first taVNS treatment at the baseline time, and underwent resting-state MRI scanning twice, before and after taVNS. All the patients then received taVNS treatments for 4 weeks. The severity of depression was assessed by the 17-item Hamilton Depression Rating Scale (HAMD) at the baseline time and after 4-week’s treatment. Pearson analysis was used to assess the correlation between alterations of ReHo and changes of the HAMD scores. Two patients were excluded due to excessive head movement, two patients lack clinical data in the fourth week, thus, imaging analysis was performed in 20 patients, while correlation analysis between clinical and imaging data was performed in only 18 patients.

Results

There were significant differences in the ReHo values in first-episode, drug-naïve MDD patients between pre- or post- taVNS. The primary finding is that the patients exhibited a significantly lower ReHo in the left/right median cingulate cortex, the left precentral gyrus, the left postcentral gyrus, the right calcarine cortex, the left supplementary motor area, the left paracentral lobule, and the right lingual gyrus. Pearson analysis revealed a positive correlation between changes of ReHo in the right median cingulate cortex/the left supplementary motor area and changes of HAMD scores after taVNS.

Conclusion

The decreased ReHo were found after taVNS. The sensorimotor, limbic and visual-related brain regions may play an important role in understanding the underlying neural mechanisms and be the target brain regions in the further therapy.

Parcellating the human brain using resting-state dynamic functional connectivity

Brain cartography has expanded substantially over the past decade. In this regard, resting-state functional connectivity (FC) plays a key role in identifying the locations of putative functional borders. However, scant attention has been paid to the dynamic nature of functional interactions in the human brain. Indeed, FC is typically assumed to be stationary across time, which may obscure potential or subtle functional boundaries, particularly in regions with high flexibility and adaptability. In this study, we developed a dynamic FC (dFC)-based parcellation framework, established a new functional human brain atlas termed D-BFA (DFC-based Brain Functional Atlas), and verified its neurophysiological plausibility by stereo-EEG data. As the first dFC-based whole-brain atlas, the proposed D-BFA delineates finer functional boundaries that cannot be captured by static FC, and is further supported by good correspondence with cytoarchitectonic areas and task activation maps. Moreover, the D-BFA reveals the spatial distribution of dynamic variability across the brain and generates more homogenous parcels compared with most alternative parcellations. Our results demonstrate the superiority and practicability of dFC in brain parcellation, providing a new template to exploit brain topographic organization from a dynamic perspective. The D-BFA will be publicly available for download at https://github.com/sliderplm/D-BFA-618.

White matter substrates of functional connectivity dynamics in the human brain

The contribution of structural connectivity to functional connectivity dynamics is still far from being elucidated. Herein, we applied track-weighted dynamic functional connectivity (tw-dFC), a model integrating structural, functional, and dynamic connectivity, on high quality diffusion weighted imaging and resting-state fMRI data from two independent repositories. The tw-dFC maps were analyzed using independent component analysis, aiming at identifying spatially independent white matter components which support dynamic changes in functional connectivity. Each component consisted of a spatial map of white matter bundles that show consistent fluctuations in functional connectivity at their endpoints, and a time course representative of such functional activity. These components show high intra-subject, inter-subject, and inter-cohort reproducibility. We provided also converging evidence that functional information about white matter activity derived by this method can capture biologically meaningful features of brain connectivity organization, as well as predict higher-order cognitive performance.

A novel 5D brain parcellation approach based on spatio-temporal encoding of resting fMRI data from deep residual learning

OBJECTIVE

Brain parcellation is an essential aspect of computational neuroimaging research and deals with segmenting the brain into (possibly overlapping) sub-regions employed to study brain anatomy or function. In the context of functional parcellation, brain organization which is often measured via temporal metrics such as coherence, is highly dynamic. This dynamic aspect is ignored in most research, which typically applies anatomically based, fixed regions for each individual, and can produce misleading results.

METHODS

In this work, we propose a novel spatio-temporal-network (5D) brain parcellation scheme utilizing a deep residual network to predict the probability of each voxel belonging to a brain network at each point in time.

RESULTS

We trained 53 4D brain networks and evaluate the ability of these networks to capture spatial and temporal dynamics as well as to show sensitivity to individual or group-level variation (in our case with age).

CONCLUSION

The proposed system generates informative spatio-temporal networks that vary not only across individuals but also over time and space.

SIGNIFICANCE

The dynamic 5D nature of the developed approach provides a powerful framework that expands on existing work and has potential to identify novel and typically ignored findings when studying the healthy and disordered brain.

Revealing the Relevant Spatiotemporal Scale Underlying Whole-Brain Dynamics

The brain rapidly processes and adapts to new information by dynamically transitioning between whole-brain functional networks. In this whole-brain modeling study we investigate the relevance of spatiotemporal scale in whole-brain functional networks. This is achieved through estimating brain parcellations at different spatial scales (100-900 regions) and time series at different temporal scales (from milliseconds to seconds) generated by a whole-brain model fitted to fMRI data. We quantify the richness of the dynamic repertoire at each spatiotemporal scale by computing the entropy of transitions between whole-brain functional networks. The results show that the optimal relevant spatial scale is around 300 regions and a temporal scale of around 150 ms. Overall, this study provides much needed evidence for the relevant spatiotemporal scales and recommendations for analyses of brain dynamics.

A detailed analysis of anatomical plausibility of crossed and uncrossed streamline rendition of the dentato-rubro-thalamic tract (DRT(T)) in a commercial stereotactic planning system

Background

An increasing number of neurosurgeons use display of the dentato-rubro-thalamic tract (DRT) based on diffusion weighted imaging (dMRI) as basis for their routine planning of stimulation or lesioning approaches in stereotactic tremor surgery. An evaluation of the anatomical validity of the display of the DRT with respect to modern stereotactic planning systems and across different tracking environments has not been performed.

Methods

Distinct dMRI and anatomical magnetic resonance imaging (MRI) data of high and low quality from 9 subjects were used. Six subjects had repeated MRI scans and therefore entered the analysis twice. Standardized DICOM structure templates for volume of interest definition were applied in native space for all investigations. For tracking BrainLab Elements (BrainLab, Munich, Germany), two tensor deterministic tracking (FT2), MRtrix IFOD2 (https://www.mrtrix.org), and a global tracking (GT) approach were used to compare the display of the uncrossed (DRTu) and crossed (DRTx) fiber structure after transformation into MNI space. The resulting streamlines were investigated for congruence, reproducibility, anatomical validity, and penetration of anatomical way point structures.

Results

In general, the DRTu can be depicted with good quality (as judged by waypoints). FT2 (surgical) and GT (neuroscientific) show high congruence. While GT shows partly reproducible results for DRTx, the crossed pathway cannot be reliably reconstructed with the other (iFOD2 and FT2) algorithms.

Conclusion

Since a direct anatomical comparison is difficult in the individual subjects, we chose a comparison with two research tracking environments as the best possible “ground truth.” FT2 is useful especially because of its manual editing possibilities of cutting erroneous fibers on the single subject level. An uncertainty of 2 mm as mean displacement of DRTu is expectable and should be respected when using this approach for surgical planning. Tractographic renditions of the DRTx on the single subject level seem to be still illusive.

Hidden Markov Modeling Reveals Prolonged "Baseline" State and Shortened Antagonistic State across the Adult Lifespan

The brain networks undergo functional reorganization across the whole lifespan, but the dynamic patterns behind the reorganization remain largely unclear. This study models the dynamics of spontaneous activity of large-scale networks using hidden Markov model (HMM), and investigates how it changes with age on two adult lifespan datasets of 176/157 subjects (aged 20-80 years). Results for both datasets showed that 1) older adults tended to spend less time on a state where default mode network (DMN) and attentional networks show antagonistic activity, 2) older adults spent more time on a "baseline" state with moderate-level activation of all networks, accompanied with lower transition probabilities from this state to the others and higher transition probabilities from the others to this state, and 3) HMM exhibited higher sensitivity in uncovering the age effects compared with temporal clustering method. Our results suggest that the aging brain is characterized by the shortening of the antagonistic instances between DMN and attention systems, as well as the prolongation of the inactive period of all networks, which might reflect the shift of the dynamical working point near criticality in older adults.

Brain parcellation driven by dynamic functional connectivity better capture intrinsic network dynamics

Until now, dynamic functional connectivity (dFC) based on functional magnetic resonance imaging is typically estimated on a set of predefined regions of interest (ROIs) derived from an anatomical or static functional atlas which follows an implicit assumption of functional homogeneity within ROIs underlying temporal fluctuation of functional coupling, potentially leading to biases or underestimation of brain network dynamics. Here, we presented a novel computational method based on dynamic functional connectivity degree (dFCD) to derive meaningful brain parcellations that can capture functional homogeneous regions in temporal variance of functional connectivity. Several spatially distributed but functionally meaningful areas that are well consistent with known intrinsic connectivity networks were identified through independent component analysis (ICA) of time‐varying dFCD maps. Furthermore, a systematical comparison with commonly used brain atlases, including the Anatomical Automatic Labeling template, static ICA‐driven parcellation and random parcellation, demonstrated that the ROI‐definition strategy based on the proposed dFC‐driven parcellation could better capture the interindividual variability in dFC and predict observed individual cognitive performance (e.g., fluid intelligence, cognitive flexibility, and sustained attention) based on chronnectome. Together, our findings shed new light on the functional organization of resting brains at the timescale of seconds and emphasized the significance of a dFC‐driven and voxel‐wise functional homogeneous parcellation for network dynamics analyses in neuroscience.