Fine-grain atlases of functional modes for fMRI analysis

Beyond the chronic pain stage: default mode network perturbation depends on years lived with back pain

ABSTRACT

Research has indicated that the default mode network (DMN) is perturbated in patients with chronic pain when compared with healthy controls, and this perturbation is correlated with the duration of pain during the chronic pain stage. It remains unclear whether DMN adaptations manifest during the subacute pain stage and progress over time because of the duration of pain experience, rather than being a specific correlate of the chronic pain stage. Furthermore, information regarding whether these adaptations are related to cognitive processes of adaptation is lacking. To this end, we examined the DMN in 31 patients with chronic back pain (CBP), 77 patients with subacute back pain (SBP), as well as 39 healthy pain-free controls (HC) applying a graph-theoretic network approach on functional resting-state magnetic resonance imaging. Beyond the comparison between groups, we used a linear analysis considering the years lived with pain (YLP) across all patients with back pain and additionally performed a mediation analysis of the role of cognitive pain coping. In line with previous studies, we found significant DMN perturbation in CBP compared with HC. However, this did not apply to the comparison of CBP with SBP. Instead, we observed a positive correlation between DMN perturbation and YLP. This was significantly mediated by coping attitudes towards pain. Default mode network perturbation may thus reflect neural adaptation processes to pain experience rather than a single correlate of the chronic pain stage and be modulated by cognitive adaption. This points to potentially underinvestigated significant adaptation processes that could enable more fine-grained patient stratification.

Gradients of Brain Organization: Smooth Sailing from Methods Development to User Community

Multimodal neuroimaging grants a powerful in vivo window into the structure and function of the human brain. Recent methodological and conceptual advances have enabled investigations of the interplay between large-scale spatial trends - or gradients - in brain structure and function, offering a framework to unify principles of brain organization across multiple scales. Strong community enthusiasm for these techniques has been instrumental in their widespread adoption and implementation to answer key questions in neuroscience. Following a brief review of current literature on this framework, this perspective paper will highlight how pragmatic steps aiming to make gradient methods more accessible to the community propelled these techniques to the forefront of neuroscientific inquiry. More specifically, we will emphasize how interest for gradient methods was catalyzed by data sharing, open-source software development, as well as the organization of dedicated workshops led by a diverse team of early career researchers. To this end, we argue that the growing excitement for brain gradients is the result of coordinated and consistent efforts to build an inclusive community and can serve as a case in point for future innovations and conceptual advances in neuroinformatics. We close this perspective paper by discussing challenges for the continuous refinement of neuroscientific theory, methodological innovation, and real-world translation to maintain our collective progress towards integrated models of brain organization.

Post-treatment alterations in white matter integrity in PTSD: Effects on symptoms and functional connectivity a secondary analysis of an RCT

Post-traumatic stress disorder (PTSD) has been linked to altered communication within the limbic system, including reduced structural connectivity in the uncinate fasciculus (UNC; i.e., decreased fractional anisotropy; FA) and reduced resting-state functional connectivity (RSFC) between the hippocampus and ventromedial prefrontal cortex (vmPFC). Previous research has demonstrated attenuation of PTSD symptoms and alterations in RSFC following exposure-based psychotherapy. However, the relationship between changes in structural and functional connectivity patterns and PTSD symptoms following treatment remains unclear. To investigate this, we conducted a secondary analysis of data from a randomized clinical trial of intensive exposure therapy, evaluating alterations in UNC FA, hippocampus-vmPFC RSFC, and PTSD symptoms before (pre-treatment), 7 days after (post-treatment), and 30 days after (follow-up) the completion of therapy. Our results showed that post-treatment changes in RSFC were positively correlated with post-treatment and follow-up changes in UNC FA and that post-treatment changes in UNC FA were positively correlated with post-treatment and follow-up changes in PTSD symptoms. These findings suggest that early changes in functional connectivity are associated with sustained changes in anatomical connectivity, which in turn are linked to reduced PTSD symptom severity.

Introducing HiBoP: a Unity-based visualization software for large iEEG datasets

BACKGROUND

Intracranial EEG data offer a unique spatio-temporal precision to investigate human brain functions. Large datasets have become recently accessible thanks to new iEEG data-sharing practices and tighter collaboration with clinicians. Yet, the complexity of such datasets poses new challenges, especially regarding the visualization and anatomical display of iEEG.

NEW METHOD

We introduce HiBoP, a multi-modal visualization software specifically designed for large groups of patients and multiple experiments. Its main features include the dynamic display of iEEG responses induced by tasks/stimulations, the definition of Regions and electrodes Of Interest, and the shift between group-level and individual-level 3D anatomo-functional data.

RESULTS

We provide a use-case with data from 36 patients to reveal the global cortical dynamics following tactile stimulation. We used HiBoP to visualize high-gamma responses [50-150 Hz], and define three major response components in primary somatosensory and premotor cortices and parietal operculum.

COMPARISON WITH EXISTING METHODS(S)

Several iEEG softwares are now publicly available with outstanding analysis features. Yet, most were developed in languages (Python/Matlab) chosen to facilitate the inclusion of new analysis by users, rather than the quality of the visualization. HiBoP represents a visualization tool developed with videogame standards (Unity/C#), and performs detailed anatomical analysis rapidly, across multiple conditions, patients, and modalities with an easy export toward third-party softwares.

CONCLUSION

HiBoP provides a user-friendly environment that greatly facilitates the exploration of large iEEG datasets, and helps users decipher subtle structure/function relationships.

Contrastive voxel clustering for multiscale modeling of brain network

Resting-state functional magnetic resonance imaging (fMRI) provides an efficient way to analyze the functional connectivity between brain regions. A comprehensive understanding of brain functionality requires a unified description of multi-scale layers of neural structure. However, existing brain network modeling methods often simplify this property by averaging Blood oxygen level dependent (BOLD) signals at the brain region level for fMRI-based analysis with the assumption that BOLD signals are homogeneous within each brain region, which ignores the heterogeneity of voxels within each Region of Interest (ROI). This study introduces a novel multi-stage self-supervised learning framework for multiscale brain network analysis, which effectively delineates brain functionality from voxel to ROIs and up to sample level. A Contrastive Voxel Clustering (CVC) module is proposed to simultaneously learn the voxel-level features and clustering assignments, which ensures the retention of informative clustering features at the finest voxel-level and concurrently preserves functional connectivity characteristics. Additionally, based on the extracted features and clustering assignments at the voxel level by CVC, a Brain ROI-based Graph Neural Network (BR-GNN) is built to extract functional connectivity features at the brain ROI-level and used for sample-level prediction, which integrates the functional clustering maps with the pre-established structural ROI maps and creates a more comprehensive and effective analytical tool. Experiments are performed on two datasets, which illustrate the effectiveness and generalization ability of the proposed method by analyzing voxel-level clustering results and brain ROIs-level functional characteristics. The proposed method provides a multiscale modeling framework for brain functional connectivity analysis, which will be further used for other brain disease identification. Code is available at https://github.com/yanliugroup/fmri-cvc.

fMRI-based spatio-temporal parcellations of the human brain

PURPOSE OF REVIEW

Human brain parcellation based on functional magnetic resonance imaging (fMRI) plays an essential role in neuroscience research. By segmenting vast and intricate fMRI data into functionally similar units, researchers can better decipher the brain's structure in both healthy and diseased states. This article reviews current methodologies and ideas in this field, while also outlining the obstacles and directions for future research.

RECENT FINDINGS

Traditional brain parcellation techniques, which often rely on cytoarchitectonic criteria, overlook the functional and temporal information accessible through fMRI. The adoption of machine learning techniques, notably deep learning, offers the potential to harness both spatial and temporal information for more nuanced brain segmentation. However, the search for a one-size-fits-all solution to brain segmentation is impractical, with the choice between group-level or individual-level models and the intended downstream analysis influencing the optimal parcellation strategy. Additionally, evaluating these models is complicated by our incomplete understanding of brain function and the absence of a definitive "ground truth".

SUMMARY

While recent methodological advancements have significantly enhanced our grasp of the brain's spatial and temporal dynamics, challenges persist in advancing fMRI-based spatio-temporal representations. Future efforts will likely focus on refining model evaluation and selection as well as developing methods that offer clear interpretability for clinical usage, thereby facilitating further breakthroughs in our comprehension of the brain.

Structurally-constrained encoding framework using a multi-voxel reduced-rank latent model for human natural vision

Objective. Voxel-wise visual encoding models based on convolutional neural networks (CNNs) have emerged as one of the prominent predictive tools of human brain activity via functional magnetic resonance imaging signals. While CNN-based models imitate the hierarchical structure of the human visual cortex to generate explainable features in response to natural visual stimuli, there is still a need for a brain-inspired model to predict brain responses accurately based on biomedical data.Approach. To bridge this gap, we propose a response prediction module called the Structurally Constrained Multi-Output (SCMO) module to include homologous correlations that arise between a group of voxels in a cortical region and predict more accurate responses.Main results. This module employs all the responses across a visual area to predict individual voxel-wise BOLD responses and therefore accounts for the population activity and collective behavior of voxels. Such a module can determine the relationships within each visual region by creating a structure matrix that represents the underlying voxel-to-voxel interactions. Moreover, since each response module in visual encoding tasks relies on the image features, we conducted experiments using two different feature extraction modules to assess the predictive performance of our proposed module. Specifically, we employed a recurrent CNN that integrates both feedforward and recurrent interactions, as well as the popular AlexNet model that utilizes feedforward connections.Significance.We demonstrate that the proposed framework provides a reliable predictive ability to generate brain responses across multiple areas, outperforming benchmark models in terms of stability and coherency of features.

Network mechanisms of ongoing brain activity's influence on conscious visual perception

Sensory inputs enter a constantly active brain, whose state is always changing from one moment to the next. Currently, little is known about how ongoing, spontaneous brain activity participates in online task processing. We employed 7 Tesla fMRI and a threshold-level visual perception task to probe the effects of prestimulus ongoing brain activity on perceptual decision-making and conscious recognition. Prestimulus activity originating from distributed brain regions, including visual cortices and regions of the default-mode and cingulo-opercular networks, exerted a diverse set of effects on the sensitivity and criterion of conscious recognition, and categorization performance. We further elucidate the mechanisms underlying these behavioral effects, revealing how prestimulus activity modulates multiple aspects of stimulus processing in highly specific and network-dependent manners. These findings reveal heretofore unknown network mechanisms underlying ongoing brain activity's influence on conscious perception, and may hold implications for understanding the precise roles of spontaneous activity in other brain functions.

Electroconvulsive therapy disrupts functional connectivity between hippocampus and posterior default mode network

BACKGROUND

The mechanisms underlying memory deficits after electroconvulsive therapy (ECT) remain unclear but altered functional interactions between hippocampus and neocortex may play a role.

OBJECTIVES

To test whether ECT reduces functional connectivity between hippocampus and posterior regions of the default mode network (DMN) and to examine whether altered hippocampal-neocortical functional connectivity correlates with memory impairment. A secondary aim was to explore if these connectivity changes are present 6 months after ECT.

METHODS

In-patients with severe depression (n = 35) received bitemporal ECT. Functional connectivity of the hippocampus was probed with resting-state fMRI before the first ECT-session, after the end of ECT, and at a six-month follow-up. Memory was assessed with the Verbal Learning Test - Delayed Recall. Seed-based connectivity analyses established connectivity of four hippocampal seeds, covering the anterior and posterior parts of the right and left hippocampus.

RESULTS

Compared to baseline, three of four hippocampal seeds became less connected to the core nodes of the posterior DMN in the week after ECT with Cohen's d ranging from -0.9 to -1.1. At the group level, patients showed post-ECT memory impairment, but individual changes in delayed recall were not correlated with the reduction in hippocampus-DMN connectivity. At six-month follow-up, no significant hippocampus-DMN reductions in connectivity were evident relative to pre-ECT, and memory scores had returned to baseline.

CONCLUSION

ECT leads to a temporary disruption of functional hippocampus-DMN connectivity in patients with severe depression, but the change in connectivity strength is not related to the individual memory impairment.

Systematic evaluation of fMRI data-processing pipelines for consistent functional connectomics

Functional interactions between brain regions can be viewed as a network, enabling neuroscientists to investigate brain function through network science. Here, we systematically evaluate 768 data-processing pipelines for network reconstruction from resting-state functional MRI, evaluating the effect of brain parcellation, connectivity definition, and global signal regression. Our criteria seek pipelines that minimise motion confounds and spurious test-retest discrepancies of network topology, while being sensitive to both inter-subject differences and experimental effects of interest. We reveal vast and systematic variability across pipelines' suitability for functional connectomics. Inappropriate choice of data-processing pipeline can produce results that are not only misleading, but systematically so, with the majority of pipelines failing at least one criterion. However, a set of optimal pipelines consistently satisfy all criteria across different datasets, spanning minutes, weeks, and months. We provide a full breakdown of each pipeline's performance across criteria and datasets, to inform future best practices in functional connectomics.

Systematic Review Between Resting-State fMRI and Task fMRI in Planning for Brain Tumour Surgery

As an alternative to task-based functional magnetic resonance imaging (T-fMRI), resting-state functional magnetic resonance imaging (Rs-fMRI) is suggested for preoperative mapping of patients with brain tumours, with an emphasis on treatment guidance and neurodegeneration prediction. A systematic review was conducted of 18 recent studies involving 1035 patients with brain tumours and Rs-fMRI protocols. This was accomplished by searching the electronic databases PubMed, Scopus, and Web of Science. For clinical benefit, we compared Rs-fMRI to standard T-fMRI and intraoperative direct cortical stimulation (DCS). The results of Rs-fMRI and T-fMRI were compared and their correlation with intraoperative DCS results was examined through a systematic review. Our exhaustive investigation demonstrated that Rs-fMRI is a dependable and sensitive preoperative mapping technique that detects neural networks in the brain with precision and identifies crucial functional regions in agreement with intraoperative DCS. Rs-fMRI comes in handy, especially in situations where T-fMRI proves to be difficult because of patient-specific factors. Additionally, our exhaustive investigation demonstrated that Rs-fMRI is a valuable tool in the preoperative screening and evaluation of brain tumours. Furthermore, its capability to assess brain function, forecast surgical results, and enhance decision-making may render it applicable in the clinical management of brain tumours.

Multivariate association between psychosocial environment, behaviors, and brain functional networks in adolescent depression

BACKGROUND

Adolescent depression shows high clinical heterogeneity. Brain functional networks serve as a powerful tool for investigating neural mechanisms underlying depression profiles. A key challenge is to characterize how variation in brain functional organization links to behavioral features and psychosocial environmental influences.

METHODS

We recruited 80 adolescents with major depressive disorder (MDD) and 42 healthy controls (HCs). First, we estimated the differences in functional connectivity of resting-state networks (RSN) between the two groups. Then, we used sparse canonical correlation analysis to characterize patterns of associations between RSN connectivity and symptoms, cognition, and psychosocial environmental factors in MDD adolescents. Clustering analysis was applied to stratify patients into homogenous subtypes according to these brain-behavior-environment associations.

RESULTS

MDD adolescents showed significantly hyperconnectivity between the ventral attention and cingulo-opercular networks compared with HCs. We identified one reliable pattern of covariation between RSN connectivity and clinical/environmental features in MDD adolescents. In this pattern, psychosocial factors, especially the interpersonal and family relationships, were major contributors to variation in connectivity of salience, cingulo-opercular, ventral attention, subcortical and somatosensory-motor networks. Based on this association, we categorized patients into two subgroups which showed different environment and symptoms characteristics, and distinct connectivity alterations. These differences were covered up when the patients were taken as a whole group.

CONCLUSION

This study identified the environmental exposures associated with specific functional networks in MDD youths. Our findings emphasize the importance of the psychosocial context in assessing brain function alterations in adolescent depression and have the potential to promote targeted treatment and precise prevention.

hyve, a compositional visualisation engine for brain imaging data

Visualisations facilitate the interpretation of geometrically structured data and results. However, heterogeneous geometries-such as volumes, surfaces, and networks-have traditionally mandated different software approaches. We introduce hyve, a Python library that uses a compositional functional framework to enable parametric implementation of custom visualisations for different brain geometries. Under this framework, users compose a reusable visualisation protocol from geometric primitives for representing data geometries, input primitives for common data formats and research objectives, and output primitives for producing interactive displays or configurable snapshots. hyve also writes documentation for user-constructed protocols, automates serial production of multiple visualisations, and includes an API for semantically organising an editable multi-panel figure. Through the seamless composition of input, output, and geometric primitives, hyve supports creating visualisations for a range of neuroimaging research objectives.

Replication of a neuroimaging biomarker for striatal dysfunction in psychosis

To bring biomarkers closer to clinical application, they should be generalizable, reliable, and maintain performance within the constraints of routine clinical conditions. The functional striatal abnormalities (FSA), is among the most advanced neuroimaging biomarkers in schizophrenia, trained to discriminate diagnosis, with post-hoc analyses indicating prognostic properties. Here, we attempt to replicate its diagnostic capabilities measured by the area under the curve (AUC) in receiver operator characteristic curves discriminating individuals with psychosis (n = 101) from healthy controls (n = 51) in the Human Connectome Project for Early Psychosis. We also measured the test-retest (run 1 vs 2) and phase encoding direction (i.e., AP vs PA) reliability with intraclass correlation coefficients (ICC). Additionally, we measured effects of scan length on classification accuracy (i.e., AUCs) and reliability (i.e., ICCs). Finally, we tested the prognostic capability of the FSA by the correlation between baseline scores and symptom improvement over 12 weeks of antipsychotic treatment in a separate cohort (n = 97). Similar analyses were conducted for the Yeo networks intrinsic connectivity as a reference. The FSA had good/excellent diagnostic discrimination (AUC = 75.4%, 95% CI = 67.0-83.3%; in non-affective psychosis AUC = 80.5%, 95% CI = 72.1-88.0%, and in affective psychosis AUC = 58.7%, 95% CI = 44.2-72.0%). Test-retest reliability ranged between ICC = 0.48 (95% CI = 0.35-0.59) and ICC = 0.22 (95% CI = 0.06-0.36), which was comparable to that of networks intrinsic connectivity. Phase encoding direction reliability for the FSA was ICC = 0.51 (95% CI = 0.42-0.59), generally lower than for networks intrinsic connectivity. By increasing scan length from 2 to 10 min, diagnostic classification of the FSA increased from AUC = 71.7% (95% CI = 63.1-80.3%) to 75.4% (95% CI = 67.0-83.3%) and phase encoding direction reliability from ICC = 0.29 (95% CI = 0.14-0.43) to ICC = 0.51 (95% CI = 0.42-0.59). FSA scores did not correlate with symptom improvement. These results reassure that the FSA is a generalizable diagnostic - but not prognostic - biomarker. Given the replicable results of the FSA as a diagnostic biomarker trained on case-control datasets, next the development of prognostic biomarkers should be on treatment-response data.

Frequency-specific segregation and integration of human cerebral cortex: An intrinsic functional atlas

The cognitive and behavioral functions of the human brain are supported by its frequency multiplexing mechanism. However, there is limited understanding of the dynamics of the functional network topology. This study aims to investigate the frequency-specific topology of the functional human brain using 7T rs-fMRI data. Frequency-specific parcellations were first performed, revealing frequency-dependent dynamics within the frontoparietal control, parietal memory, and visual networks. An intrinsic functional atlas containing 456 parcels was proposed and validated using stereo-EEG. Graph theory analysis suggested that, in addition to the task-positive vs. task-negative organization observed in static networks, there was a cognitive control system additionally from a frequency perspective. The reproducibility and plausibility of the identified hub sets were confirmed through 3T fMRI analysis, and their artificial removal had distinct effects on network topology. These results indicate a more intricate and subtle dynamics of the functional human brain and emphasize the significance of accurate topography.

Continuous evaluation of denoising strategies in resting-state fMRI connectivity using fMRIPrep and Nilearn

Reducing contributions from non-neuronal sources is a crucial step in functional magnetic resonance imaging (fMRI) connectivity analyses. Many viable strategies for denoising fMRI are used in the literature, and practitioners rely on denoising benchmarks for guidance in the selection of an appropriate choice for their study. However, fMRI denoising software is an ever-evolving field, and the benchmarks can quickly become obsolete as the techniques or implementations change. In this work, we present a denoising benchmark featuring a range of denoising strategies, datasets and evaluation metrics for connectivity analyses, based on the popular fMRIprep software. The benchmark prototypes an implementation of a reproducible framework, where the provided Jupyter Book enables readers to reproduce or modify the figures on the Neurolibre reproducible preprint server (https://neurolibre.org/). We demonstrate how such a reproducible benchmark can be used for continuous evaluation of research software, by comparing two versions of the fMRIprep. Most of the benchmark results were consistent with prior literature. Scrubbing, a technique which excludes time points with excessive motion, combined with global signal regression, is generally effective at noise removal. Scrubbing was generally effective, but is incompatible with statistical analyses requiring the continuous sampling of brain signal, for which a simpler strategy, using motion parameters, average activity in select brain compartments, and global signal regression, is preferred. Importantly, we found that certain denoising strategies behave inconsistently across datasets and/or versions of fMRIPrep, or had a different behavior than in previously published benchmarks. This work will hopefully provide useful guidelines for the fMRIprep users community, and highlight the importance of continuous evaluation of research methods.

New organizational principles and 3D cytoarchitectonic maps of the dorsolateral prefrontal cortex in the human brain

Areas of the dorsolateral prefrontal cortex (DLPFC) are part of the frontoparietal control, default mode, salience, and ventral attention networks. The DLPFC is involved in executive functions, like working memory, value encoding, attention, decision-making, and behavioral control. This functional heterogeneity is not reflected in existing neuroanatomical maps. For example, previous cytoarchitectonic studies have divided the DLPFC into two or four areas. Macroanatomical parcellations of this region rely on gyri and sulci, which are not congruent with cytoarchitectonic parcellations. Therefore, this study aimed to provide a microstructural analysis of the human DLPFC and 3D maps of cytoarchitectonic areas to help address the observed functional variability in studies of the DLPFC. We analyzed ten human post-mortem brains in serial cell-body stained brain sections and mapped areal boundaries using a statistical image analysis approach. Five new areas (i.e., SFG2, SFG3, SFG4, MFG4, and MFG5) were identified on the superior and middle frontal gyrus, i.e., regions corresponding to parts of Brodmann areas 9 and 46. Gray level index profiles were used to determine interregional cytoarchitectural differences. The five new areas were reconstructed in 3D, and probability maps were generated in commonly used reference spaces, considering the variability of areas in stereotaxic space. Hierarchical cluster analysis revealed a high degree of similarity within the identified DLPFC areas while neighboring areas (frontal pole, Broca's region, area 8, and motoric areas) were separable. Comparisons with functional imaging studies revealed specific functional profiles of the DLPFC areas. Our results indicate that the new areas do not follow a simple organizational gradient assumption in the DLPFC. Instead, they are more similar to those of the ventrolateral prefrontal cortex (Broca's areas 44, 45) and frontopolar areas (Fp1, Fp2) than to the more posterior areas. Within the DLPFC, the cytoarchitectonic similarities between areas do not seem to follow a simple anterior-to-posterior gradient either, but cluster along other principles. The new maps are part of the publicly available Julich Brain Atlas and provide a microstructural reference for existing and future imaging studies. Thus, our study represents a further step toward deciphering the structural-functional organization of the human prefrontal cortex.

Adolescents' pain-related ontogeny shares a neural basis with adults' chronic pain in basothalamo-cortical organization

During late adolescence, the brain undergoes ontogenic organization altering subcortical-cortical circuitry. This includes regions implicated in pain chronicity, and thus alterations in the adolescent ontogenic organization could predispose to pain chronicity in adulthood - however, evidence is lacking. Using resting-state functional magnetic resonance imaging from a large European longitudinal adolescent cohort and an adult cohort with and without chronic pain, we examined links between painful symptoms and brain connectivity. During late adolescence, thalamo-, caudate-, and red nucleus-cortical connectivity were positively and subthalamo-cortical connectivity negatively associated with painful symptoms. Thalamo-cortical connectivity, but also subthalamo-cortical connectivity, was increased in adults with chronic pain compared to healthy controls. Our results indicate a shared basis in basothalamo-cortical circuitries between adolescent painful symptomatology and adult pain chronicity, with the subthalamic pathway being differentially involved, potentially due to a hyperconnected thalamo-cortical pathway in chronic pain and ontogeny-driven organization. This can inform neuromodulation-based prevention and early intervention.

Should one go for individual- or group-level brain parcellations? A deep-phenotyping benchmark

The analysis and understanding of brain characteristics often require considering region-level information rather than voxel-sampled data. Subject-specific parcellations have been put forward in recent years, as they can adapt to individual brain organization and thus offer more accurate individual summaries than standard atlases. However, the price to pay for adaptability is the lack of group-level consistency of the data representation. Here, we investigate whether the good representations brought by individualized models are merely an effect of circular analysis, in which individual brain features are better represented by subject-specific summaries, or whether this carries over to new individuals, i.e., whether one can actually adapt an existing parcellation to new individuals and still obtain good summaries in these individuals. For this, we adapt a dictionary-learning method to produce brain parcellations. We use it on a deep-phenotyping dataset to assess quantitatively the patterns of activity obtained under naturalistic and controlled-task-based settings. We show that the benefits of individual parcellations are substantial, but that they vary a lot across brain systems.

Riemannian frameworks for the harmonization of resting-state functional MRI scans

Magnetic Resonance Imaging provides unprecedented images of the brain. Unfortunately, scanners and acquisition protocols can significantly impact MRI scans. The development of statistical methods able to reduce this variability without altering the relevant information in the scans, often coined harmonization methods, has been the topic of an increasing research effort supported by the recent growth of publicly available neuroimaging data sets and new possibilities for combining them to achieve greater statistical power. In this work, we focus on the challenges specifically raised by the harmonization of resting-state functional MRI scans. We propose to harmonize resting-state fMRI scans by reducing the impact of covariates such as scanner differences and scanning protocols on their associated functional connectomes and then propagating the changes back to the rs-fMRI time series. We use Riemannian geometric frameworks to preserve the mathematical properties of functional connectomes during their harmonization, and we demonstrate how state-of-the-art harmonization methods can be embedded within these frameworks to reduce covariates effects while preserving the relevant clinical information associated with aging or brain disorders. During our experiments, a large set of synthetic data was generated and processed to compare eighty variants of the proposed approach. The framework achieving the best harmonization was then applied to three low-dimensional data sets made of 712 sets of fMRI time series provided by the ABIDE consortium and two high-dimensional data sets obtained by processing 1527 rs-fMRI scans provided by the Human Connectome Project, the Framingham Heart Study and the Genetics of Brain Structure and Function study. These experiments established that our new framework could successfully harmonize low-dimensional connectomes and voxelwise functional time series and confirmed the need for preserving connectomes properties during their harmonization.

Brain decoding of the Human Connectome Project tasks in a dense individual fMRI dataset

Brain decoding aims to infer cognitive states from patterns of brain activity. Substantial inter-individual variations in functional brain organization challenge accurate decoding performed at the group level. In this paper, we tested whether accurate brain decoding models can be trained entirely at the individual level. We trained several classifiers on a dense individual functional magnetic resonance imaging (fMRI) dataset for which six participants completed the entire Human Connectome Project (HCP) task battery >13 times over ten separate fMRI sessions. We evaluated nine decoding methods, from Support Vector Machines (SVM) and Multi-Layer Perceptron (MLP) to Graph Convolutional Neural Networks (GCN). All decoders were trained to classify single fMRI volumes into 21 experimental conditions simultaneously, using ∼7 h of fMRI data per participant. The best prediction accuracies were achieved with GCN and MLP models, whose performance (57-67 % accuracy) approached state-of-the-art accuracy (76 %) with models trained at the group level on >1 K hours of data from the original HCP sample. Our SVM model also performed very well (54-62 % accuracy). Feature importance maps derived from MLP -our best-performing model- revealed informative features in regions relevant to particular cognitive domains, notably in the motor cortex. We also observed that inter-subject classification achieved substantially lower accuracy than subject-specific models, indicating that our decoders learned individual-specific features. This work demonstrates that densely-sampled neuroimaging datasets can be used to train accurate brain decoding models at the individual level. We expect this work to become a useful benchmark for techniques that improve model generalization across multiple subjects and acquisition conditions.

Differences between subclinical attention-deficit/hyperactivity and autistic traits in default mode, salience, and frontoparietal network connectivities in young adult Japanese

Attention-deficit hyperactivity disorder (ADHD) and autism spectrum disorder (ASD) are associated with attentional impairments, with both commonalities and differences in the nature of their attention deficits. This study aimed to investigate the neural correlates of ADHD and ASD traits in healthy individuals, focusing on the functional connectivity (FC) of attention-related large-scale brain networks (LSBNs). The participants were 61 healthy individuals (30 men; age, 21.9 ± 1.9 years). The Adult ADHD Self-Report Scale (ASRS) and Autism Spectrum Quotient (AQ) were administered as indicators of ADHD and ASD traits, respectively. Performance in the continuous performance test (CPT) was used as a behavioural measure of sustained attentional function. Functional magnetic resonance imaging scans were performed during the resting state (Rest) and auditory oddball task (Odd). Considering the critical role in attention processing, we focused our analyses on the default mode (DMN), frontoparietal (FPN), and salience (SN) networks. Region of interest (ROI)-to-ROI analyses (false discovery rate < 0.05) were performed to determine relationships between psychological measures with within-network FC (DMN, FPN, and SN) as well as with between-network FC (DMN-FPN, DMN-SN, and FPN-SN). ASRS scores, but not AQ scores, were correlated with less frequent commission errors and shorter reaction times in the CPT. During Odd, significant positive correlations with ASRS were demonstrated in multiple FCs within DMN, while significant positive correlations with AQ were demonstrated in multiple FCs within FPN. AQs were negatively correlated with FPN-SN FCs. During Rest, AQs were negatively and positively correlated with one FC within the SN and multiple FCs between the DMN and SN, respectively. These findings of the ROI-to-ROI analysis were only partially replicated in a split-half replication analysis, a replication analysis with open-access data sets, and a replication analysis with a structure-based atlas. The better CPT performance by individuals with subclinical ADHD traits suggests positive effects of these traits on sustained attention. Differential associations between LSBN FCs and ASD/ADHD traits corroborate the notion of differences in sustained and selective attention between clinical ADHD and ASD.

Lack of reproducibility of resting-state functional MRI findings in migraine with aura

BACKGROUND

Several studies have applied resting-state functional MRI to examine whether functional brain connectivity is altered in migraine with aura patients. These studies had multiple limitations, including small sample sizes, and reported conflicting results. Here, we performed a large, cross-sectional brain imaging study to reproduce previous findings.

METHODS

We recruited women aged 30-60 years from the nationwide Danish Twin Registry. Resting-state functional MRI of women with migraine with aura, their co-twins, and unrelated migraine-free twins was performed at a single centre. We carried out an extensive series of brain connectivity data analyses. Patients were compared to migraine-free controls and to co-twins.

RESULTS

Comparisons were based on data from 160 patients, 30 co-twins, and 136 controls. Patients were similar to controls with regard to age, and several lifestyle characteristics. We replicated clear effects of age on resting-state networks. In contrast, we failed to detect any differences, and to replicate previously reported differences, in functional connectivity between migraine patients with aura and non-migraine controls or their co-twins in any of the analyses.

CONCLUSION

Given the large sample size and the unbiased population-based design of our study, we conclude that women with migraine with aura have normal resting-state brain connectivity outside of migraine attacks.

A personalized cortical atlas for functional regions of interest

Advances in functional MRI (fMRI) allow mapping an individual's brain function in vivo. Task fMRI can localize domain-specific regions of cognitive processing or functional regions of interest (fROIs) within an individual. Moreover, data from resting state (no task) fMRI can be used to define an individual's connectome, which can characterize that individual's functional organization via connectivity-based parcellations. However, can connectivity-based parcellations alone predict an individual's fROIs? Here, we describe an approach to compute individualized rs-fROIs (i.e., regions that correspond to given fROI constructed using only resting state data) for motor control, working memory, high-level vision, and language comprehension. The rs-fROIs were computed and validated using a large sample of young adults (n = 1,018) with resting state and task fMRI from the Human Connectome Project. First, resting state parcellations were defined across a sequence of resolutions from broadscale to fine-grained networks in a training group of 500 individuals. Second, 21 rs-fROIs were defined from the training group by identifying the rs network that most closely matched task-defined fROIs across all individuals. Third, the selectivity of rs-fROIs was investigated in a training set of the remaining 518 individuals. All computed rs-fROIs were indeed selective for their preferred category. Critically, the rs-fROIs had higher selectivity than probabilistic atlas parcels for nearly all fROIs. In conclusion, we present a potential approach to define selective fROIs on an individual-level circumventing the need for multiple task-based localizers.NEW & NOTEWORTHY We compute individualized resting state parcels that identify an individual's own functional regions of interest (fROIs) for high-level vision, language comprehension, motor control, and working memory, using only their functional connectome. This approach demonstrates a rapid and powerful alternative for finding a large set of fROIs in an individual, using only their unique connectivity pattern, which does not require the costly acquisition of multiple fMRI localizer tasks.

Long term structural and functional neural changes following a single infusion of Ketamine in PTSD

NMDA receptor antagonists have a vital role in extinction, learning, and reconsolidation processes. During the reconsolidation window, memories are activated into a labile state and can be reconsolidated in an altered form. This concept might have significant clinical implications in treating PTSD. In this pilot study we tested the potential of a single infusion of ketamine, followed by brief exposure therapy, to enhance post-retrieval extinction of PTSD trauma memories. 27 individuals diagnosed with PTSD were randomly assigned to receive either ketamine (0.5 mg/kg 40 min; N = 14) or midazolam (0.045 mg/kg; N = 13) after retrieval of the traumatic memory. 24 h following infusion, participants received a four-day trauma-focused psychotherapy. Symptoms and brain activity were assessed before treatment, at the end of treatment, and at 30-day follow-up. Amygdala activation to trauma scripts (a major biomarker of fear response) served as the main study outcome. Although PTSD symptoms improved equally in both groups, post-treatment, ketamine recipients showed a lower amygdala (-0.33, sd = 0.13, 95%HDI [-0.56,-0.04]) and hippocampus (-0.3 (sd = 0.19), 95%HDI [-0.65, 0.04]; marginal effect) reactivation to trauma memories, compared to midazolam recipients. Post-retrieval ketamine administration was also associated with decreased connectivity between the amygdala and hippocampus (-0.28, sd = 0.11, 95%HDI [-0.46, -0.11]), with no change in amygdala-vmPFC connectivity. Moreover, reduction in fractional anisotropy in bi-lateral uncinate fasciculus was seen in the Ketamine recipients compared with the midazolam recipients (right: post-treatment: -0.01108, 95% HDI [-0.0184,-0.003]; follow-up: -0.0183, 95% HDI [-0.02719,-0.0107]; left: post-treatment: -0.019, 95% HDI [-0.028,-0.011]; follow-up: -0.017, 95% HDI [-0.026,-0.007]). Taken together it is possible that ketamine may enhance post-retrieval extinction of the original trauma memories in humans. These preliminary findings show promising direction toward the capacity to rewrite human traumatic memories and modulate the fear response for at least 30 days post-extinction. When combined with psychotherapy for PTSD, further investigation of ketamine dose, timing of administration, and frequency of administration, is warranted.

Homogeneous-Multiset-CCA-Based Brain Covariation and Contravariance Connectivity Network Modeling

Brain connectivity networks based on functional magnetic resonance imaging (fMRI) have expanded our understanding of brain functions in both healthy and diseased states. However, most current studies construct connectivity networks using averaged regional time courses with the strong assumption that the activities of voxels contained in each brain region are similar, ignoring their possible variations. Additionally, pairwise correlation analysis is often adopted with more attention to positive relationships, while joint interactions at the network level as well as anti-correlations are less investigated. In this paper, to provide a new strategy for regional activity representation and brain connectivity modeling, a novel homogeneous multiset canonical correlation analysis (HMCCA) model is proposed, which enforces sign constraints on the weights of voxels to guarantee homogeneity within each brain region. It is capable of obtaining regional representative signals and constructing covariation and contravariance networks simultaneously, at both group and subject levels. Validations on two sessions of fMRI data verified its reproducibility and reliability when dealing with brain connectivity networks. Further experiments on subjects with and without Parkinson's disease (PD) revealed significant alterations in brain connectivity patterns, which were further associated with clinical scores and demonstrated superior prediction ability, indicating its potential in clinical practice.

An atlas of trait associations with resting-state and task-evoked human brain functional organizations in the UK Biobank

Functional magnetic resonance imaging (fMRI) has been widely used to identify brain regions linked to critical functions, such as language and vision, and to detect tumors, strokes, brain injuries, and diseases. It is now known that large sample sizes are necessary for fMRI studies to detect small effect sizes and produce reproducible results. Here we report a systematic association analysis of 647 traits with imaging features extracted from resting-state and task-evoked fMRI data of more than 40,000 UK Biobank participants. We used a parcellation-based approach to generate 64,620 functional connectivity measures to reveal fine-grained details about cerebral cortex functional organizations. The difference between functional organizations at rest and during task was examined, and we have prioritized important brain regions and networks associated with a variety of human traits and clinical outcomes. For example, depression was most strongly associated with decreased connectivity in the somatomotor network. We have made our results publicly available and developed a browser framework to facilitate the exploration of brain function-trait association results (http://fmriatlas.org/).

Replication of a neuroimaging biomarker for striatal dysfunction in psychosis

To bring biomarkers closer to clinical application, they should be generalizable, reliable, and maintain performance within the constraints of routine clinical conditions. The functional striatal abnormalities (FSA), is among the most advanced neuroimaging biomarkers in schizophrenia, trained to discriminate diagnosis, with post-hoc analyses indicating prognostic properties. Here, we attempt to replicate its diagnostic capabilities measured by the area under the curve (AUC) in receiver operator characteristic curves discriminating individuals with psychosis (n=101) from healthy controls (n=51) in the Human Connectome Project for Early Psychosis. We also measured the test-retest (run 1 vs 2) and phase encoding direction (i.e., AP vs PA) reliability with intraclass correlation coefficients (ICC). Additionally, we measured effects of scan length on classification accuracy (i.e., AUCs) and reliability (i.e., ICCs). Finally, we tested the prognostic capability of the FSA by the correlation between baseline scores and symptom improvement over 12 weeks of antipsychotic treatment in a separate cohort (n=97). Similar analyses were conducted for the Yeo networks intrinsic connectivity as a reference. The FSA had good/excellent diagnostic discrimination (AUC=75.4%, 95%CI=67.0%-83.3%; in non-affective psychosis AUC=80.5%, 95%CI=72.1-88.0%, and in affective psychosis AUC=58.7%, 95%CI=44.2-72.0%). Test-retest reliability ranged between ICC=0.48 (95%CI=0.35-0.59) and ICC=0.22 (95%CI=0.06-0.36), which was comparable to that of networks intrinsic connectivity. Phase encoding direction reliability for the FSA was ICC=0.51 (95%CI=0.42-0.59), generally lower than for networks intrinsic connectivity. By increasing scan length from 2 to 10 minutes, diagnostic classification of the FSA increased from AUC=71.7% (95%CI=63.1%-80.3%) to 75.4% (95%CI=67.0%-83.3%) and phase encoding direction reliability from ICC=0.29 (95%CI=0.14-0.43) to ICC=0.51 (95%CI=0.42-0.59). FSA scores did not correlate with symptom improvement. These results reassure that the FSA is a generalizable diagnostic - but not prognostic - biomarker. Given the replicable results of the FSA as a diagnostic biomarker trained on case-control datasets, next the development of prognostic biomarkers should be on treatment-response data.

AtOM, an ontology model to standardize use of brain atlases in tools, workflows, and data infrastructures

Brain atlases are important reference resources for accurate anatomical description of neuroscience data. Open access, three-dimensional atlases serve as spatial frameworks for integrating experimental data and defining regions-of-interest in analytic workflows. However, naming conventions, parcellation criteria, area definitions, and underlying mapping methodologies differ considerably between atlases and across atlas versions. This lack of standardized description impedes use of atlases in analytic tools and registration of data to different atlases. To establish a machine-readable standard for representing brain atlases, we identified four fundamental atlas elements, defined their relations, and created an ontology model. Here we present our Atlas Ontology Model (AtOM) and exemplify its use by applying it to mouse, rat, and human brain atlases. We discuss how AtOM can facilitate atlas interoperability and data integration, thereby increasing compliance with the FAIR guiding principles. AtOM provides a standardized framework for communication and use of brain atlases to create, use, and refer to specific atlas elements and versions. We argue that AtOM will accelerate analysis, sharing, and reuse of neuroscience data.

Replication of a neuroimaging biomarker for striatal dysfunction in psychosis

To bring biomarkers closer to clinical application, they should be generalizable, reliable, and maintain performance within the constraints of routine clinical conditions. The functional striatal abnormalities (FSA), is among the most advanced neuroimaging biomarkers in schizophrenia, trained to discriminate diagnosis, with post-hoc analyses indicating prognostic properties. Here, we attempt to replicate its diagnostic capabilities measured by the area under the curve (AUC) in receiver operator characteristic curves discriminating individuals with psychosis (n=101) from healthy controls (n=51) in the Human Connectome Project for Early Psychosis. We also measured the test-retest (run 1 vs 2) and phase encoding direction (i.e., AP vs PA) reliability with intraclass correlation coefficients (ICC). Additionally, we measured effects of scan length on classification accuracy (i.e., AUCs) and reliability (i.e., ICCs). Finally, we tested the prognostic capability of the FSA by the correlation between baseline scores and symptom improvement over 12 weeks of antipsychotic treatment in a separate cohort (n=97). Similar analyses were conducted for the Yeo networks intrinsic connectivity as a reference. The FSA had good/excellent diagnostic discrimination (AUC=75.4%, 95%CI=67.0%-83.3%; in non-affective psychosis AUC=80.5%, 95%CI=72.1-88.0%, and in affective psychosis AUC=58.7%, 95%CI=44.2-72.0%). Test-retest reliability ranged between ICC=0.48 (95%CI=0.35-0.59) and ICC=0.22 (95%CI=0.06-0.36), which was comparable to that of networks intrinsic connectivity. Phase encoding direction reliability for the FSA was ICC=0.51 (95%CI=0.42-0.59), generally lower than for networks intrinsic connectivity. By increasing scan length from 2 to 10 minutes, diagnostic classification of the FSA increased from AUC=71.7% (95%CI=63.1%-80.3%) to 75.4% (95%CI=67.0%-83.3%) and phase encoding direction reliability from ICC=0.29 (95%CI=0.14-0.43) to ICC=0.51 (95%CI=0.42-0.59). FSA scores did not correlate with symptom improvement. These results reassure that the FSA is a generalizable diagnostic - but not prognostic - biomarker. Given the replicable results of the FSA as a diagnostic biomarker trained on case-control datasets, next the development of prognostic biomarkers should be on treatment-response data.

Large-scale brain networks and intra-axial tumor surgery: a narrative review of functional mapping techniques, critical needs, and scientific opportunities

In recent years, a paradigm shift in neuroscience has been occurring from "localizationism," or the idea that the brain is organized into separately functioning modules, toward "connectomics," or the idea that interconnected nodes form networks as the underlying substrates of behavior and thought. Accordingly, our understanding of mechanisms of neurological function, dysfunction, and recovery has evolved to include connections, disconnections, and reconnections. Brain tumors provide a unique opportunity to probe large-scale neural networks with focal and sometimes reversible lesions, allowing neuroscientists the unique opportunity to directly test newly formed hypotheses about underlying brain structural-functional relationships and network properties. Moreover, if a more complete model of neurological dysfunction is to be defined as a "disconnectome," potential avenues for recovery might be mapped through a "reconnectome." Such insight may open the door to novel therapeutic approaches where previous attempts have failed. In this review, we briefly delve into the most clinically relevant neural networks and brain mapping techniques, and we examine how they are being applied to modern neurosurgical brain tumor practices. We then explore how brain tumors might teach us more about mechanisms of global brain dysfunction and recovery through pre- and postoperative longitudinal connectomic and behavioral analyses.

Continuous Evaluation of Denoising Strategies in Resting-State fMRI Connectivity Using fMRIPrep and Nilearn

Reducing contributions from non-neuronal sources is a crucial step in functional magnetic resonance imaging (fMRI) connectivity analyses. Many viable strategies for denoising fMRI are used in the literature, and practitioners rely on denoising benchmarks for guidance in the selection of an appropriate choice for their study. However, fMRI denoising software is an ever-evolving field, and the benchmarks can quickly become obsolete as the techniques or implementations change. In this work, we present a denoising benchmark featuring a range of denoising strategies, datasets and evaluation metrics for connectivity analyses, based on the popular fMRIprep software. The benchmark is implemented in a fully reproducible framework, where the provided research objects enable readers to reproduce or modify core computations, as well as the figures of the article using the Jupyter Book project and the Neurolibre reproducible preprint server (https://neurolibre.org/). We demonstrate how such a reproducible benchmark can be used for continuous evaluation of research software, by comparing two versions of the fMRIprep software package. The majority of benchmark results were consistent with prior literature. Scrubbing, a technique which excludes time points with excessive motion, combined with global signal regression, is generally effective at noise removal. Scrubbing however disrupts the continuous sampling of brain images and is incompatible with some statistical analyses, e.g. auto-regressive modeling. In this case, a simple strategy using motion parameters, average activity in select brain compartments, and global signal regression should be preferred. Importantly, we found that certain denoising strategies behave inconsistently across datasets and/or versions of fMRIPrep, or had a different behavior than in previously published benchmarks. This work will hopefully provide useful guidelines for the fMRIprep users community, and highlight the importance of continuous evaluation of research methods. Our reproducible benchmark infrastructure will facilitate such continuous evaluation in the future, and may also be applied broadly to different tools or even research fields.

Brain functional connectivity mirrors genetic pleiotropy in psychiatric conditions

Pleiotropy occurs when a genetic variant influences more than one trait. This is a key property of the genomic architecture of psychiatric disorders and has been observed for rare and common genomic variants. It is reasonable to hypothesize that the microscale genetic overlap (pleiotropy) across psychiatric conditions and cognitive traits may lead to similar overlaps at the macroscale brain level such as large-scale brain functional networks. We took advantage of brain connectivity, measured by resting-state functional MRI to measure the effects of pleiotropy on large-scale brain networks, a putative step from genes to behaviour. We processed nine resting-state functional MRI datasets including 32 726 individuals and computed connectome-wide profiles of seven neuropsychiatric copy-number-variants, five polygenic scores, neuroticism and fluid intelligence as well as four idiopathic psychiatric conditions. Nine out of 19 pairs of conditions and traits showed significant functional connectivity correlations (rFunctional connectivity), which could be explained by previously published levels of genomic (rGenetic) and transcriptomic (rTranscriptomic) correlations with moderate to high concordance: rGenetic-rFunctional connectivity = 0.71 [0.40-0.87] and rTranscriptomic-rFunctional connectivity = 0.83 [0.52; 0.94]. Extending this analysis to functional connectivity profiles associated with rare and common genetic risk showed that 30 out of 136 pairs of connectivity profiles were correlated above chance. These similarities between genetic risks and psychiatric disorders at the connectivity level were mainly driven by the overconnectivity of the thalamus and the somatomotor networks. Our findings suggest a substantial genetic component for shared connectivity profiles across conditions and traits, opening avenues to delineate general mechanisms-amenable to intervention-across psychiatric conditions and genetic risks.

Identification of overlapping and interacting networks reveals intrinsic spatiotemporal organization of the human brain

The human brain is a complex network that exhibits dynamic fluctuations in activity across space and time. Depending on the analysis method, canonical brain networks identified from resting-state fMRI (rs-fMRI) are typically constrained to be either orthogonal or statistically independent in their spatial and/or temporal domains. We avoid imposing these potentially unnatural constraints through the combination of a temporal synchronization process ("BrainSync") and a three-way tensor decomposition method ("NASCAR") to jointly analyze rs-fMRI data from multiple subjects. The resulting set of interacting networks comprises minimally constrained spatiotemporal distributions, each representing one component of functionally coherent activity across the brain. We show that these networks can be clustered into six distinct functional categories and naturally form a representative functional network atlas for a healthy population. This functional network atlas could help explore group and individual differences in neurocognitive function, as we demonstrate in the context of ADHD and IQ prediction.

Comparison Between Gradients and Parcellations for Functional Connectivity Prediction of Behavior

Resting-state functional connectivity (RSFC) is widely used to predict behavioral measures. To predict behavioral measures, representing RSFC with parcellations and gradients are the two most popular approaches. Here, we compare parcellation and gradient approaches for RSFC-based prediction of a broad range of behavioral measures in the Human Connectome Project (HCP) and Adolescent Brain Cognitive Development (ABCD) datasets. Among the parcellation approaches, we consider group-average "hard" parcellations (Schaefer et al., 2018), individual-specific "hard" parcellations (Kong et al., 2021a), and an individual-specific "soft" parcellation (spatial independent component analysis with dual regression; Beckmann et al., 2009). For gradient approaches, we consider the well-known principal gradients (Margulies et al., 2016) and the local gradient approach that detects local RSFC changes (Laumann et al., 2015). Across two regression algorithms, individual-specific hard-parcellation performs the best in the HCP dataset, while the principal gradients, spatial independent component analysis and group-average "hard" parcellations exhibit similar performance. On the other hand, principal gradients and all parcellation approaches perform similarly in the ABCD dataset. Across both datasets, local gradients perform the worst. Finally, we find that the principal gradient approach requires at least 40 to 60 gradients to perform as well as parcellation approaches. While most principal gradient studies utilize a single gradient, our results suggest that incorporating higher order gradients can provide significant behaviorally relevant information. Future work will consider the inclusion of additional parcellation and gradient approaches for comparison.

Supervised Phenotype Discovery From Multimodal Brain Imaging

Data-driven discovery of image-derived phenotypes (IDPs) from large-scale multimodal brain imaging data has enormous potential for neuroscientific and clinical research by linking IDPs to subjects' demographic, behavioural, clinical and cognitive measures (i.e., non-imaging derived phenotypes or nIDPs). However, current approaches are primarily based on unsupervised approaches, without the use of information in nIDPs. In this paper, we proposed a semi-supervised, multimodal, and multi-task fusion approach, termed SuperBigFLICA, for IDP discovery, which simultaneously integrates information from multiple imaging modalities as well as multiple nIDPs. SuperBigFLICA is computationally efficient and largely avoids the need for parameter tuning. Using the UK Biobank brain imaging dataset with around 40,000 subjects and 47 modalities, along with more than 17,000 nIDPs, we showed that SuperBigFLICA enhances the prediction power of nIDPs, benchmarked against IDPs derived by conventional expert-knowledge and unsupervised-learning approaches (with average nIDP prediction accuracy improvements of up to 46%). It also enables the learning of generic imaging features that can predict new nIDPs. Further empirical analysis of the SuperBigFLICA algorithm demonstrates its robustness in different prediction tasks and the ability to derive biologically meaningful IDPs in predicting health outcomes and cognitive nIDPs, such as fluid intelligence and hypertension.

Disentangling self from pain: mindfulness meditation-induced pain relief is driven by thalamic-default mode network decoupling

ABSTRACT

For millenniums, mindfulness was believed to diminish pain by reducing the influence of self-appraisals of noxious sensations. Today, mindfulness meditation is a highly popular and effective pain therapy that is believed to engage multiple, nonplacebo-related mechanisms to attenuate pain. Recent evidence suggests that mindfulness meditation-induced pain relief is associated with the engagement of unique cortico-thalamo-cortical nociceptive filtering mechanisms. However, the functional neural connections supporting mindfulness meditation-based analgesia remain unknown. This mechanistically focused clinical trial combined functional magnetic resonance imaging with psychophysical pain testing (49°C stimulation and pain visual analogue scales) to identify the neural connectivity supporting the direct modulation of pain-related behavioral and neural responses by mindfulness meditation. We hypothesized that mindfulness meditation-based pain relief would be reflected by greater decoupling between brain mechanisms supporting appraisal (prefrontal) and nociceptive processing (thalamus). After baseline pain testing, 40 participants were randomized to a well-validated, 4-session mindfulness meditation or book-listening regimen. Functional magnetic resonance imaging and noxious heat (49°C; right calf) were combined during meditation to test study hypotheses. Mindfulness meditation significantly reduced behavioral and neural pain responses when compared to the controls. Preregistered (NCT03414138) whole-brain analyses revealed that mindfulness meditation-induced analgesia was moderated by greater thalamus-precuneus decoupling and ventromedial prefrontal deactivation, respectively, signifying a pain modulatory role across functionally distinct neural mechanisms supporting self-referential processing. Two separate preregistered seed-to-seed analyses found that mindfulness meditation-based pain relief was also associated with weaker contralateral thalamic connectivity with the prefrontal and primary somatosensory cortex, respectively. Thus, we propose that mindfulness meditation is associated with a novel self-referential nociceptive gating mechanism to reduce pain.

A quantified comparison of cortical atlases on the basis of trait morphometricity

BACKGROUND

Many different brain atlases exist that subdivide the human cortex into dozens or hundreds of regions-of-interest (ROIs). Inconsistency across studies using one or another cortical atlas may contribute to the replication crisis across the neurosciences.

METHODS

Here, we provide a quantitative comparison between seven popular cortical atlases (Yeo, Desikan-Killiany, Destrieux, Jülich-Brain, Gordon, Glasser, Schaefer) and vertex-wise measures (thickness, surface area, and volume), to determine which parcellation retains the most information in the analysis of behavioural traits (incl. age, sex, body mass index, and cognitive ability) in the UK Biobank sample (N∼40,000). We use linear mixed models to compare whole-brain morphometricity; the proportion of trait variance accounted for when using a given atlas.

RESULTS

Commonly-used atlases resulted in a considerable loss of information compared to vertex-wise representations of cortical structure. Morphometricity increased linearly as a function of the log-number of ROIs included in an atlas, indicating atlas-based analyses miss many true associations and yield limited prediction accuracy. Likelihood ratio tests revealed that low-dimensional atlases accounted for unique trait variance rather than variance common between atlases, suggesting that previous studies likely returned atlas-specific findings. Finally, we found that the commonly-used atlases yielded brain-behaviour associations on par with those obtained with random parcellations, where specific region boundaries were randomly generated.

DISCUSSION

Our findings motivate future structural neuroimaging studies to favour vertex-wise cortical representations over coarser atlases, or to consider repeating analyses across multiple atlases, should the use of low-dimensional atlases be necessary. The insights uncovered here imply that cortical atlas choices likely contribute to the lack of reproducibility in ROI-based studies.

Performance scaling for structural MRI surface parcellations: a machine learning analysis in the ABCD Study

The use of predefined parcellations on surface-based representations of the brain as a method for data reduction is common across neuroimaging studies. In particular, prediction-based studies typically employ parcellation-driven summaries of brain measures as input to predictive algorithms, but the choice of parcellation and its influence on performance is often ignored. Here we employed preprocessed structural magnetic resonance imaging (sMRI) data from the Adolescent Brain Cognitive Development Study® to examine the relationship between 220 parcellations and out-of-sample predictive performance across 45 phenotypic measures in a large sample of 9- to 10-year-old children (N = 9,432). Choice of machine learning (ML) pipeline and use of alternative multiple parcellation-based strategies were also assessed. Relative parcellation performance was dependent on the spatial resolution of the parcellation, with larger number of parcels (up to ~4,000) outperforming coarser parcellations, according to a power-law scaling of between 1/4 and 1/3. Performance was further influenced by the type of parcellation, ML pipeline, and general strategy, with existing literature-based parcellations, a support vector-based pipeline, and ensembling across multiple parcellations, respectively, as the highest performing. These findings highlight the choice of parcellation as an important influence on downstream predictive performance, showing in some cases that switching to a higher resolution parcellation can yield a relatively large boost to performance.

Meta-analysis of the functional neuroimaging literature with probabilistic logic programming

Inferring reliable brain-behavior associations requires synthesizing evidence from thousands of functional neuroimaging studies through meta-analysis. However, existing meta-analysis tools are limited to investigating simple neuroscience concepts and expressing a restricted range of questions. Here, we expand the scope of neuroimaging meta-analysis by designing NeuroLang: a domain-specific language to express and test hypotheses using probabilistic first-order logic programming. By leveraging formalisms found at the crossroads of artificial intelligence and knowledge representation, NeuroLang provides the expressivity to address a larger repertoire of hypotheses in a meta-analysis, while seamlessly modeling the uncertainty inherent to neuroimaging data. We demonstrate the language's capabilities in conducting comprehensive neuroimaging meta-analysis through use-case examples that address questions of structure-function associations. Specifically, we infer the specific functional roles of three canonical brain networks, support the role of the visual word-form area in visuospatial attention, and investigate the heterogeneous organization of the frontoparietal control network.

Interpreting mental state decoding with deep learning models

In mental state decoding, researchers aim to identify the set of mental states (e.g., experiencing happiness or fear) that can be reliably identified from the activity patterns of a brain region (or network). Deep learning (DL) models are highly promising for mental state decoding because of their unmatched ability to learn versatile representations of complex data. However, their widespread application in mental state decoding is hindered by their lack of interpretability, difficulties in applying them to small datasets, and in ensuring their reproducibility and robustness. We recommend approaching these challenges by leveraging recent advances in explainable artificial intelligence (XAI) and transfer learning, and also provide recommendations on how to improve the reproducibility and robustness of DL models in mental state decoding.

The Executive Functioning Paradox in Substance Use Disorders

Deficits in neurocognitive functioning are trait-like vulnerabilities that have been widely studied in persons with substance use disorders (SUD), but their role in the craving-use association and relapse vulnerability remains poorly understood. The main objectives of this study were to examine whether executive capacities moderate the magnitude of the craving-substance use relationship, and if this influence is correlated with the functional connectivity of cerebral networks, combining rsfMRI examinations and ecological momentary assessment (EMA). Eighty-six patients beginning outpatient treatment for alcohol, tobacco or cannabis addiction and 40 healthy controls completed neuropsychological tests followed by EMA to collect real-time data on craving. Fifty-four patients and 30 healthy controls also completed a resting-state fMRI before the EMA. Among the patients with SUD, better verbal fluency and resistance to interference capacities were associated with a greater propensity to use substances when the individual was experiencing craving. Preliminary rsfMRI results identified specific networks that interacted with executive performance capacities to influence the magnitude of the craving-use association. Individuals with better executive functioning may be more prone to relapse after craving episodes. Specifically, better resistance to interference and cognitive flexibility skills may reduce attention to distracting stimuli, leading to a greater awareness of craving and susceptibility to use substances.

Functional gradients in the human lateral prefrontal cortex revealed by a comprehensive coordinate-based meta-analysis

The lateral prefrontal cortex (LPFC) of humans enables flexible goal-directed behavior. However, its functional organization remains actively debated after decades of research. Moreover, recent efforts aiming to map the LPFC through meta-analysis are limited, either in scope or in the inferred specificity of structure-function associations. These limitations are in part due to the limited expressiveness of commonly-used data analysis tools, which restricts the breadth and complexity of questions that can be expressed in a meta-analysis. Here, we adopt NeuroLang, a novel approach to more expressive meta-analysis based on probabilistic first-order logic programming, to infer the organizing principles of the LPFC from 14,371 neuroimaging studies. Our findings reveal a rostrocaudal and a dorsoventral gradient, respectively explaining the most and second most variance in meta-analytic connectivity across the LPFC. Moreover, we identify a unimodal-to-transmodal spectrum of coactivation patterns along with a concrete-to-abstract axis of structure-function associations extending from caudal to rostral regions of the LPFC. Finally, we infer inter-hemispheric asymmetries along the principal rostrocaudal gradient, identifying hemisphere-specific associations with topics of language, memory, response inhibition, and sensory processing. Overall, this study provides a comprehensive meta-analytic mapping of the LPFC, grounding future hypothesis generation on a quantitative overview of past findings.

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.

Cytoarchitecture, intersubject variability, and 3D mapping of four new areas of the human anterior prefrontal cortex

The dorsolateral prefrontal cortex (DLPFC) plays a key role in cognitive control and executive functions, including working memory, attention, value encoding, decision making, monitoring, and controlling behavioral strategies. However, the relationships between this variety of functions and the underlying cortical areas, which specifically contribute to these functions, are not yet well-understood. Existing microstructural maps differ in the number, localization, and extent of areas of the DLPFC. Moreover, there is a considerable intersubject variability both in the sulcal pattern and in the microstructure of this region, which impedes comparison with functional neuroimaging studies. The aim of this study was to provide microstructural, cytoarchitectonic maps of the human anterior DLPFC in 3D space. Therefore, we analyzed 10 human post-mortem brains and mapped their borders using a well-established approach based on statistical image analysis. Four new areas (i.e., SFS1, SFS2, MFG1, and MFG2) were identified in serial, cell-body stained brain sections that occupy the anterior superior frontal sulcus and middle frontal gyrus, i.e., a region corresponding to parts of Brodmann areas 9 and 46. Differences between areas in cytoarchitecture were captured using gray level index profiles, reflecting changes in the volume fraction of cell bodies from the surface of the brain to the cortex-white matter border. A hierarchical cluster analysis of these profiles indicated that areas of the anterior DLPFC displayed higher cytoarchitectonic similarity between each other than to areas of the neighboring frontal pole (areas Fp1 and Fp2), Broca's region (areas 44 and 45) of the ventral prefrontal cortex, and posterior DLPFC areas (8d1, 8d2, 8v1, and 8v2). Area-specific, cytoarchitectonic differences were found between the brains of males and females. The individual areas were 3D-reconstructed, and probability maps were created in the MNI Colin27 and ICBM152casym reference spaces to take the variability of areas in stereotaxic space into account. The new maps contribute to Julich-Brain and are publicly available as a resource for studying neuroimaging data, helping to clarify the functional and organizational principles of the human prefrontal cortex.

Comprehensive decoding mental processes from Web repositories of functional brain images

Associating brain systems with mental processes requires statistical analysis of brain activity across many cognitive processes. These analyses typically face a difficult compromise between scope-from domain-specific to system-level analysis-and accuracy. Using all the functional Magnetic Resonance Imaging (fMRI) statistical maps of the largest data repository available, we trained machine-learning models that decode the cognitive concepts probed in unseen studies. For this, we leveraged two comprehensive resources: NeuroVault-an open repository of fMRI statistical maps with unconstrained annotations-and Cognitive Atlas-an ontology of cognition. We labeled NeuroVault images with Cognitive Atlas concepts occurring in their associated metadata. We trained neural networks to predict these cognitive labels on tens of thousands of brain images. Overcoming the heterogeneity, imbalance and noise in the training data, we successfully decoded more than 50 classes of mental processes on a large test set. This success demonstrates that image-based meta-analyses can be undertaken at scale and with minimal manual data curation. It enables broad reverse inferences, that is, concluding on mental processes given the observed brain activity.

Functional connectivity alterations between default mode network and occipital cortex in patients with obsessive-compulsive disorder (OCD)

Altered brain network connectivity is a potential biomarker for obsessive-compulsive disorder (OCD). A meta-analysis of resting-state MRI studies by Gürsel et al. (2018) described altered functional connectivity in OCD patients within and between the default mode network (DMN), the salience network (SN), and the frontoparietal network (FPN), as well as evidence for aberrant fronto-striatal circuitry. Here, we tested the replicability of these meta-analytic rsfMRI findings by measuring functional connectivity during resting-state fMRI in a new sample of OCD patients (n = 24) and matched controls (n = 33). We performed seed-to-voxel analyses using 30 seed regions from the prior meta-analysis. OCD patients showed reduced functional connectivity between the SN and the DMN compared to controls, replicating previous findings. We did not observe significant group differences of functional connectivity within the DMN, SN, nor FPN. Additionally, we observed reduced connectivity between the visual network to both the DMN and SN in OCD patients, in particular reduced functional connectivity between lateral parietal seeds and the left inferior lateral occipital pole. Furthermore, the right lateral parietal seed (associated with the DMN) was more strongly correlated with a cluster in the right lateral occipital cortex and precuneus (a region partly overlapping with the Dorsal Attentional Network (DAN)) in patients. Importantly, this latter finding was positively correlated to OCD symptom severity. Overall, our study partly replicated prior meta-analytic findings, highlighting hypoconnectivity between SN and DMN as a potential biomarker for OCD. Furthermore, we identified changes between the SN and the DMN with the visual network. This suggests that abnormal connectivity between cortex regions associated with abstract functions (transmodal regions such as the DMN), and cortex regions associated with constrained neural processing (unimodal regions such as the visual cortex), may be important in OCD.

The role of neural load effects in predicting individual differences in working memory function

Studies of working memory (WM) function have tended to adopt either a within-subject approach, focusing on effects of load manipulations, or a between-subjects approach, focusing on individual differences. This dichotomy extends to WM neuroimaging studies, with different neural correlates being identified for within- and between-subjects variation in WM. Here, we examined this issue in a systematic fashion, leveraging the large-sample Human Connectome Project dataset, to conduct a well-powered, whole-brain analysis of the N-back WM task. We first demonstrate the advantages of parcellation schemes for dimension reduction, in terms of load-related effect sizes. This parcel-based approach is then utilized to directly compare the relationship between load-related (within-subject) and behavioral individual differences (between-subject) effects through both correlational and predictive analyses. The results suggest a strong linkage of within-subject and between-subject variation, with larger load-effects linked to stronger brain-behavior correlations. In frontoparietal cortex no hemispheric biases were found towards one type of variation, but the Dorsal Attention Network did exhibit greater sensitivity to between over within-subjects variation, whereas in the Somatomotor network, the reverse pattern was observed. Cross-validated predictive modeling capitalizing on this tight relationship between the two effects indicated greater predictive power for load-activated than load-deactivated parcels, while also demonstrating that load-related effect size can serve as an effective guide to feature (i.e., parcel) selection, in maximizing predictive power while maintaining interpretability. Together, the findings demonstrate an important consistency across within- and between-subjects approaches to identifying the neural substrates of WM, which can be effectively harnessed to develop more powerful predictive models.

Parcellation-induced variation of empirical and simulated brain connectomes at group and subject levels

Recent developments of whole-brain models have demonstrated their potential when investigating resting-state brain activity. However, it has not been systematically investigated how alternating derivations of the empirical structural and functional connectivity, serving as the model input, from MRI data influence modeling results. Here, we study the influence from one major element: the brain parcellation scheme that reduces the dimensionality of brain networks by grouping thousands of voxels into a few hundred brain regions. We show graph-theoretical statistics derived from the empirical data and modeling results exhibiting a high heterogeneity across parcellations. Furthermore, the network properties of empirical brain connectomes explain the lion’s share of the variance in the modeling results with respect to the parcellation variation. Such a clear-cut relationship is not observed at the subject-resolved level per parcellation. Finally, the graph-theoretical statistics of the simulated connectome correlate with those of the empirical functional connectivity across parcellations. However, this relation is not one-to-one, and its precision can vary between models. Our results imply that network properties of both empirical connectomes can explain the goodness-of-fit of whole-brain models to empirical data at a global group level but not at a single-subject level, which provides further insights into the personalization of whole-brain models.

The structural and functional connectivities of the brain, which reflect the anatomical connections of axonal bundles and the amount of coactivation between brain regions, respectively, only weakly correlate with each other. In order to enhance and investigate this relationship, large-scale whole-brain dynamical models were involved in this branch of research. However, how the definitions of the brain regions parcellated according to a so-called brain atlas influence these models has so far not been systematically assessed. In this article, we show that this influence can be large, and link group-averaged, atlas-induced deviations to network properties extracted from both types of connectivity. Additionally, we demonstrate that the same association does not apply to subject-specific variations. These results may contribute to the further personalization of the whole-brain models.