Effects of spatial smoothing on group-level differences in functional brain networks

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.

Inter-regional correlation estimators for functional magnetic resonance imaging

Functional magnetic resonance imaging (fMRI) functional connectivity between brain regions is often computed using parcellations defined by functional or structural atlases. Typically, some kind of voxel averaging is performed to obtain a single temporal correlation estimate per region pair. However, several estimators can be defined for this task, with various assumptions and degrees of robustness to local noise, global noise, and region size. In this paper, we systematically present and study the properties of 9 different functional connectivity estimators taking into account the spatial structure of fMRI data, based on a simple fMRI data spatial model. These include 3 existing estimators and 6 novel estimators. We demonstrate the empirical properties of the estimators using synthetic, animal, and human data, in terms of graph structure, repeatability and reproducibility, discriminability, dependence on region size, as well as local and global noise robustness. We prove analytically the link between regional intra-correlation and inter-region correlation, and show that the choice of estimator has a strong influence on inter-correlation values. Some estimators, including the commonly used correlation of averages (ca), are positively biased, and have more dependence to region size and intra-correlation than robust alternatives, resulting in spatially-dependent bias. We define the new local correlation of averages estimator with better theoretical guarantees, lower bias, significantly lower dependence on region size (Spearman correlation 0.40 vs 0.55, paired t-test T=27.2, p=1.1e-47), at negligible cost to discriminative power, compared to the ca estimator. The difference in connectivity pattern between the estimators is not distributed uniformly throughout the brain, but rather shows a clear ventral-dorsal gradient, suggesting that region size and intra-correlation plays an important role in shaping functional networks defined using the ca estimator, and leading to non-trivial differences in their connectivity structure. We provide an open source R package and equivalent Python implementation to facilitate the use of the new estimators, together with preprocessed rat time-series.

Effects of Hyperdirect Pathway Theta Burst Transcranial Magnetic Stimulation on Inhibitory Control, Craving, and Smoking in Adults With Nicotine Dependence: A Double-Blind, Randomized Crossover Trial

BACKGROUND

Nicotine dependence is associated with dysregulated hyperdirect pathway (HDP)-mediated inhibitory control (IC). However, there are currently no evidence-based treatments that have been shown to target the HDP to improve IC and reduce cigarette cravings and smoking.

METHODS

Following a baseline nonstimulation control session, this study (N = 37; female: n = 17) used a double-blind, randomized crossover design to examine the behavioral and neural effects of intermittent theta burst stimulation (iTBS) and continuous TBS (cTBS) to the right inferior frontal gyrus (rIFG)-a key cortical node of the HDP. Associations between treatment effects were also explored.

RESULTS

At baseline, HDP IC task-state functional connectivity was positively associated with IC task performance, which confirmed the association between HDP circuit function and IC. Compared with iTBS, rIFG cTBS improved IC task performance. Compared with the baseline nonstimulation control session, both TBS conditions reduced cigarette craving and smoking; however, although craving and smoking were lower for cTBS, no differences were found between the two active conditions. In addition, although HDP IC task-state functional connectivity was greater following cTBS than iTBS, there was no significant difference between conditions. Finally, cTBS-induced improvement in IC task performance was associated with reduced craving, and cTBS-induced reduction in craving was associated with reduced smoking.

CONCLUSIONS

These findings warrant further investigation into the effects of rIFG cTBS for increasing IC and reducing craving and smoking among individuals with nicotine dependence. Future sham-controlled cTBS studies may help further elucidate the mechanisms by which rIFG cTBS mediates IC and smoking behavior.

A new causal centrality measure reveals the prominent role of subcortical structures in the causal architecture of the extended default mode network

Network representation has been an incredibly useful concept for understanding the behavior of complex systems in social sciences, biology, neuroscience, and beyond. Network science is mathematically founded on graph theory, where nodal importance is gauged using measures of centrality. Notably, recent work suggests that the topological centrality of a node should not be over-interpreted as its dynamical or causal importance in the network. Hence, identifying the influential nodes in dynamic causal models (DCM) remains an open question. This paper introduces causal centrality for DCM, a dynamics-sensitive and causally-founded centrality measure based on the notion of intervention in graphical models. Operationally, this measure simplifies to an identifiable expression using Bayesian model reduction. As a proof of concept, the average DCM of the extended default mode network (eDMN) was computed in 74 healthy subjects. Next, causal centralities of different regions were computed for this causal graph, and compared against several graph-theoretical centralities. The results showed that the subcortical structures of the eDMN were more causally central than the cortical regions, even though the graph-theoretical centralities unanimously favored the latter. Importantly, model comparison revealed that only the pattern of causal centrality was causally relevant. These results are consistent with the crucial role of the subcortical structures in the neuromodulatory systems of the brain, and highlight their contribution to the organization of large-scale networks. Potential applications of causal centrality-to study causal models of other neurotypical and pathological functional networks-are discussed, and some future lines of research are outlined.

Functional magnetic particle imaging (fMPI) of cerebrovascular changes in the rat brain during hypercapnia

Objective.Non-invasive functional brain imaging modalities are limited in number, each with its own complex trade-offs between sensitivity, spatial and temporal resolution, and the directness with which the measured signals reflect neuronal activation. Magnetic particle imaging (MPI) directly maps the cerebral blood volume (CBV), and its high sensitivity derives from the nonlinear magnetization of the superparamagnetic iron oxide nanoparticle (SPION) tracer confined to the blood pool. Our work evaluates functional MPI (fMPI) as a new hemodynamic functional imaging modality by mapping the CBV response in a rodent model where CBV is modulated by hypercapnic breathing manipulation.Approach.The rodent fMPI time-series data were acquired with a mechanically rotating field-free line MPI scanner capable of 5 s temporal resolution and 3 mm spatial resolution. The rat's CBV was modulated for 30 min with alternating 5 min hyper-/hypocapnic states, and processed using conventional fMRI tools. We compare our results to fMRI responses undergoing similar hypercapnia protocols found in the literature, and reinforce this comparison in a study of one rat with 9.4T BOLD fMRI using the identical protocol.Main results.The initial image in the time-series showed mean resting brain voxel SNR values, averaged across rats, of 99.9 following the first 10 mg kg-1SPION injection and 134 following the second. The time-series fit a conventional General Linear Model with a 15%-40% CBV change and a peak pixel CNR between 12 and 29, 2-6× higher than found in fMRI.Significance.This work introduces a functional modality with high sensitivity, although currently limited spatial and temporal resolution. With future clinical-scale development, a large increase in sensitivity could supplement other modalities and help transition functional brain imaging from a neuroscience tool focusing on population averages to a clinically relevant modality capable of detecting differences in individual patients.

Spatial Smoothing Effect on Group-Level Functional Connectivity during Resting and Task-Based fMRI

Spatial smoothing is a preprocessing step applied to neuroimaging data to enhance data quality by reducing noise and artifacts. However, selecting an appropriate smoothing kernel size can be challenging as it can lead to undesired alterations in final images and functional connectivity networks. However, there is no sufficient information about the effects of the Gaussian kernel size on group-level results for different cases yet. This study investigates the influence of kernel size on functional connectivity networks and network parameters in whole-brain rs-fMRI and tb-fMRI analyses of healthy adults. The analysis includes {0, 2, 4, 6, 8, 10} mm kernels, commonly used in practical analyses, covering all major brain networks. Graph theoretical measures such as betweenness centrality, global/local efficiency, clustering coefficient, and average path length are examined for each kernel. Additionally, principal component analysis (PCA) and independent component analysis (ICA) parameters, namely kurtosis and skewness, are evaluated for the functional images. The findings demonstrate that kernel size directly affects node connections, resulting in modifications to functional network structures and PCA/ICA parameters. However, network metrics exhibit greater resilience to these changes.

Identifying reproducible resting state networks and functional connectivity alterations following chronic restraint stress in anaesthetized rats

Background

Resting-state functional MRI (rs-fMRI) in rodent models have the potential to bridge invasive experiments and observational human studies, increasing our understanding of functional alterations in the brains of patients with depression. A major limitation in current rodent rs-fMRI studies is that there has been no consensus on healthy baseline resting-state networks (RSNs) that are reproducible in rodents. Therefore, the present study aimed to construct reproducible RSNs in a large dataset of healthy rats and then evaluate functional connectivity changes within and between these RSNs following a chronic restraint stress (CRS) model within the same animals.

Methods

A combined MRI dataset of 109 Sprague Dawley rats at baseline and after two weeks of CRS, collected during four separate experiments conducted by our lab in 2019 and 2020, was re-analysed. The mICA and gRAICAR toolbox were first applied to detect optimal and reproducible ICA components and then a hierarchical clustering algorithm (FSLNets) was applied to construct reproducible RSNs. Ridge-regularized partial correlation (FSLNets) was used to evaluate the changes in the direct connection between and within identified networks in the same animals following CRS.

Results

Four large-scale networks in anesthetised rats were identified: the DMN-like, spatial attention-limbic, corpus striatum, and autonomic network, which are homologous across species. CRS decreased the anticorrelation between DMN-like and autonomic network. CRS decreased the correlation between amygdala and a functional complex (nucleus accumbens and ventral pallidum) in the right hemisphere within the corpus striatum network. However, a high individual variability in the functional connectivity before and after CRS within RSNs was observed.

Conclusion

The functional connectivity changes detected in rodents following CRS differ from reported functional connectivity alterations in patients with depression. A simple interpretation of this difference is that the rodent response to CRS does not reflect the complexity of depression as it is experienced by humans. Nonetheless, the high inter-subject variability of functional connectivity within networks suggests that rats demonstrate different neural phenotypes, like humans. Therefore, future efforts in classifying neural phenotypes in rodents might improve the sensitivity and translational impact of models used to address aetiology and treatment of psychiatric conditions including depression.

Amygdala-based functional connectivity mediates the relationship between thought control ability and trait anxiety

Thought control ability (TCA) refers to the ability to exclude unwanted thoughts. There has been consistent evidence on the protective effect of TCA on anxiety, that higher TCA is associated with lower anxiety. However, the underlying neural mechanism remains unclear. In this study, with a large sample (N = 495), we investigated how seed-based resting-state functional connectivity (RSFC) mediates the relationship between TCA and anxiety. Our behaviour results replicated previous findings that TCA is negatively associated with trait anxiety after controlling for gender, age, and depression. More importantly, the RSFC results revealed that TCA is negatively associated with the left amygdala - left frontal pole (LA-LFP), left amygdala - left inferior temporal gyrus (LA-LITG), and left hippocampus - left inferior frontal gyrus (LH-LIFG) connectivity. In addition, a mediation analysis demonstrated that the LA-LFP and LA-LITG connectivity in particular mediated the influence of TCA on trait anxiety. Overall, our study extends previous research by revealing the neural bases underlying the protective effect of TCA on anxiety and pinpointing specific mediating RSFC pathways. Future studies could explore whether targeted TCA training (behavioural or neural) can help alleviate anxiety.

Accelerated cortical thinning precedes and predicts conversion to psychosis: The NAPLS3 longitudinal study of youth at clinical high-risk

Progressive grey matter loss has been demonstrated among clinical high-risk (CHR) individuals who convert to psychosis, but it is unknown whether these changes occur prior to psychosis onset. Identifying illness-related neurobiological mechanisms that occur prior to conversion is essential for targeted early intervention. Among participants in the third wave of the North American Prodrome Longitudinal Study (NAPLS3), this report investigated if steeper cortical thinning was observable prior to psychosis onset among CHR individuals who ultimately converted (CHR-C) and assessed the shortest possible time interval in which rates of cortical thinning differ between CHR-C, CHR non-converters (CHR-NC), and health controls (HC). 338 CHR-NC, 42 CHR-C, and 62 HC participants (age 19.3±4.2, 44.8% female, 52.5% racial/ethnic minority) completed up to 5 MRI scans across 8 months. Accelerated thinning among CHR-C compared to CHR-NC and HC was observed in multiple prefrontal, temporal, and parietal cortical regions. CHR-NC also exhibited accelerated cortical thinning compared to HC in several of these areas. Greater percent decrease in cortical thickness was observed among CHR-C compared to other groups across 2.9±1.8 months, on average, in several cortical areas. ROC analyses discriminating CHR-C from CHR-NC by percent thickness change in a left hemisphere region of interest, scanner, age, age2, and sex had an AUC of 0.74, with model predictive power driven primarily by percent thickness change. Findings indicate that accelerated cortical thinning precedes psychosis onset and differentiates CHR-C from CHR-NC and HC across short time intervals. Mechanisms underlying cortical thinning may provide novel treatment targets prior to psychosis onset.

Greater functional connectivity between the ventral frontal cortex and occipital cortex in herpes zoster patients than post-herpetic neuralgia patients

OBJECTIVE

This study aimed to compare whole brain network between herpes zoster (HZ) patients and post-herpetic neuralgia (PHN) patients, as well as to investigate the associations between whole brain network changes and pain intensity and the accuracy of classifying between different types of pain.

METHODS

PHN patients (n = 50) and HZ patients (n = 50) and healthy controls (HCs) (n = 50) underwent resting-state functional magnetic resonance imaging (rs-fMRI). Functional connectivity and global and local graph theory metrics were calculated by using Dosenbach-160 atlas. The relationship between neuroimaging indicators and clinical scales was evaluated using correlation analysis, and receiver operating characteristic (ROC) curves evaluated the feasibility of classifying PHN and HZ patients using specific neuroimaging indicators.

RESULTS

(1) 10 greater average connectivities were found in HZ group among the default mode, frontoparietal, cingulo-opercular, sensorimotor, occipital networks (ONs), and cerebellum (p < 0.001). (2) HZ patients exhibited higher global efficiency than those in the PHN and HCs (t = 2.178, p = 0.038). (3) Multiple linear regression analyses indicated that functional connectivity between the ventral frontal cortex in the cingulo-opercular network and the occipital gyrus in the ON influenced the visual analog score pain scores (β = 4.273; p = 0.004).

CONCLUSION

The variation of functional connectivity between ventral frontal cortex in the cingulo-opercular network and occipital gyrus in the ON may be a robust neuroimaging marker of the transition from HZ to PHN patients.

ADVANCES IN KNOWLEDGE

Whole-brain network analysis may be effective in distinguishing HZ and PHN patients and predicting pain intensity.

Reduced hippocampal-cortical connectivity during memory suppression predicts the ability to forget unwanted memories

The ability to suppress unwelcome memories is important for productivity and well-being. Successful memory suppression is associated with hippocampal deactivations and a concomitant disruption of this region's functionality. Much of the previous neuroimaging literature exploring such suppression-related hippocampal modulations has focused on the region's negative coupling with the prefrontal cortex. Task-based changes in functional connectivity between the hippocampus and other brain regions still need further exploration. In the present study, we utilize psychophysiological interactions and seed connectome-based predictive modeling to investigate the relationship between the hippocampus and the rest of the brain as 134 participants attempted to suppress unwanted memories during the Think/No-Think task. The results show that during retrieval suppression, the right hippocampus exhibited decreased functional connectivity with visual cortical areas (lingual and cuneus gyrus), left nucleus accumbens and the brain-stem that predicted superior forgetting of unwanted memories on later memory tests. Validation tests verified that prediction performance was not an artifact of head motion or prediction method and that the negative features remained consistent across different brain parcellations. These findings suggest that systemic memory suppression involves more than the modulation of hippocampal activity-it alters functional connectivity patterns between the hippocampus and visual cortex, leading to successful forgetting.

Selecting Multiple Node Statistics Jointly from Functional Connectivity Networks for Brain Disorders Identification

Functional connectivity networks (FCN) analysis is instructive for the diagnosis of brain diseases, such as mild cognitive impairment (MCI) and major depressive disorder (MDD) at their early stages. As the critical step of FCN analysis, feature representation provides the basis for finding potential biomarkers of brain diseases. In previous studies, different node statistics (e.g. local efficiency and local clustering coefficients) are usually extracted from FCNs as features for the diagnosis/classification task, which can specifically locate disease-related regions on the node level, so as to help us understand the neurodevelopmental roots of brain disorders. However, each node statistic is proposed only considering a kind of specific network property, which has one-sidedness and limitations. As a result, it is incomplete to represent a node with only one statistic. To resolve this issue, we put forward a novel scheme to select multiple node statistics jointly from the estimated FCNs for automated classification, called multiple node statistics feature selection (MNSFS). Specifically, we first extract multiple statistics from the same nodes and assign each kind of statistic into a group. Then, sparse group least absolute shrinkage and selection operator (sgLASSO) is used to select groups (nodes) and statistics in the groups towards a better classification performance. Such a technique enables us to simultaneously locate the discriminative brain regions, as well as the specific statistics associated with these brain regions, making the classification results more interpretable. We conducted our scheme on two public databases for identifying subjects with MCI and MDD from normal controls. Experimental results show that the proposed scheme achieves superior classification accuracy and features interpreted on the benchmark datasets.

Can we have a second helping? A preregistered direct replication study on the neurobiological mechanisms underlying self-control

Self‐control is of vital importance for human wellbeing. Hare et al. (2009) were among the first to provide empirical evidence on the neural correlates of self‐control. This seminal study profoundly impacted theory and empirical work across multiple fields. To solidify the empirical evidence supporting self‐control theory, we conducted a preregistered replication of this work. Further, we tested the robustness of the findings across analytic strategies. Participants underwent functional magnetic resonance imaging while rating 50 food items on healthiness and tastiness and making choices about food consumption. We closely replicated the original analysis pipeline and supplemented it with additional exploratory analyses to follow‐up on unexpected findings and to test the sensitivity of results to key analytical choices. Our replication data provide support for the notion that decisions are associated with a value signal in ventromedial prefrontal cortex (vmPFC), which integrates relevant choice attributes to inform a final decision. We found that vmPFC activity was correlated with goal values regardless of the amount of self‐control and it correlated with both taste and health in self‐controllers but only taste in non‐self‐controllers. We did not find strong support for the hypothesized role of left dorsolateral prefrontal cortex (dlPFC) in self‐control. The absence of statistically significant group differences in dlPFC activity during successful self‐control in our sample contrasts with the notion that dlPFC involvement is required in order to effectively integrate longer‐term goals into subjective value judgments. Exploratory analyses highlight the sensitivity of results (in terms of effect size) to the analytical strategy, for instance, concerning the approach to region‐of‐interest analysis.

Threat imminence reveals links among unfolding of anticipatory physiological response, cortical-subcortical intrinsic functional connectivity, and anxiety

Excessive expression of fear responses in anticipation of threat occurs in anxiety, but understanding of underlying pathophysiological mechanisms is limited. Animal research indicates that threat-anticipatory defensive responses are dynamically organized by threat imminence and rely on conserved circuitry. Insight from basic neuroscience research in animals on threat imminence could guide mechanistic research in humans mapping abnormal function in this circuitry to aberrant defensive responses in pathological anxiety. 50 pediatric anxiety patients and healthy-comparisons (33 females) completed an instructed threat-anticipation task whereby cues signaled delivery of painful (threat) or non-painful (safety) thermal stimulation. Temporal changes in skin-conductance indexed anxiety effects on anticipatory responding as function of threat imminence. Multivariate network analyses of resting-state functional connectivity data from a subsample were used to identify intrinsic-function correlates of anticipatory-response dynamics, within a specific, distributed network derived from translational research on defensive responding. By considering threat imminence, analyses revealed specific anxiety effects. Importantly, pathological anxiety was associated with excessive deployment of anticipatory physiological response as threat, but not safety, outcomes became more imminent. Magnitude of increase in threat-anticipatory physiological responses corresponded with magnitude of intrinsic connectivity within a cortical-subcortical circuit. Moreover, more severe anxiety was associated with stronger associations between anticipatory physiological responding and connectivity that ventromedial prefrontal cortex showed with hippocampus and basolateral amygdala, regions implicated in animal models of anxiety. These findings link basic and clinical research, highlighting variations in intrinsic function in conserved defensive circuitry as a potential pathophysiological mechanism in anxiety.

Principles and open questions in functional brain network reconstruction

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network representation involves often covert theoretical assumptions and methodological choices which affect the way networks are reconstructed from experimental data, and ultimately the resulting network properties and their interpretation. Here, we review some fundamental conceptual underpinnings and technical issues associated with brain network reconstruction, and discuss how their mutual influence concurs in clarifying the organization of brain function.

Age differences in functional network reconfiguration with working memory training

Demanding cognitive functions like working memory (WM) depend on functional brain networks being able to communicate efficiently while also maintaining some degree of modularity. Evidence suggests that aging can disrupt this balance between integration and modularity. In this study, we examined how cognitive training affects the integration and modularity of functional networks in older and younger adults. Twenty three younger and 23 older adults participated in 10 days of verbal WM training, leading to performance gains in both age groups. Older adults exhibited lower modularity overall and a greater decrement when switching from rest to task, compared to younger adults. Interestingly, younger but not older adults showed increased task-related modularity with training. Furthermore, whereas training increased efficiency within, and decreased participation of, the default-mode network for younger adults, it enhanced efficiency within a task-specific salience/sensorimotor network for older adults. Finally, training increased segregation of the default-mode from frontoparietal/salience and visual networks in younger adults, while it diffusely increased between-network connectivity in older adults. Thus, while younger adults increase network segregation with training, suggesting more automated processing, older adults persist in, and potentially amplify, a more integrated and costly global workspace, suggesting different age-related trajectories in functional network reorganization with WM training.