Insula-cortico-subcortical networks predict interoceptive awareness and stress resilience

BACKGROUND

Interoception, the neural sensing of visceral signals, and interoceptive awareness (IA), the conscious perception of interoception, are crucial for life survival functions and mental health. Resilience, the capacity to overcome adversity, has been associated with reduced interoceptive disturbances. Here, we sought evidence for our Insula Modular Active Control (IMAC) model that suggest that the insula, a brain region specialized in the processing of interoceptive information, realizes IA and contributes to resilience and mental health via cortico-subcortical connections.

METHODS

64 healthy participants (32 females; ages 18-34 years) answered questionnaires that assess IA and resilience. Mental health was evaluated with the Beck Depression Inventory II that assesses depressive mood. Participants also underwent a 15 minute resting-state functional resonance imaging session. Pearson correlations and mediation analyses were used to investigate the relationship between IA and resilience and their contributions to depressive mood. We then performed insula seed-based functional connectivity analyzes to identify insula networks involved in IA, resilience and depressive mood.

RESULTS

We first demonstrated that resilience mediates the relationship between IA and depressive mood. Second, shared and distinct intra-insula, insula-cortical and insula-subcortical networks were associated with IA, resilience and also predicted the degree of experienced depressive mood. Third, while resilience was associated with stronger insula-precuneus, insula-cerebellum and insula-prefrontal networks, IA was linked with stronger intra-insula, insula-striatum and insula-motor networks.

CONCLUSIONS

Our findings help understand the roles of insula-cortico-subcortical networks in IA and resilience. These results also highlight the potential use of insula networks as biomarkers for depression prediction.

oFVSD: a Python package of optimized forward variable selection decoder for high-dimensional neuroimaging data

The complexity and high dimensionality of neuroimaging data pose problems for decoding information with machine learning (ML) models because the number of features is often much larger than the number of observations. Feature selection is one of the crucial steps for determining meaningful target features in decoding; however, optimizing the feature selection from such high-dimensional neuroimaging data has been challenging using conventional ML models. Here, we introduce an efficient and high-performance decoding package incorporating a forward variable selection (FVS) algorithm and hyper-parameter optimization that automatically identifies the best feature pairs for both classification and regression models, where a total of 18 ML models are implemented by default. First, the FVS algorithm evaluates the goodness-of-fit across different models using the k-fold cross-validation step that identifies the best subset of features based on a predefined criterion for each model. Next, the hyperparameters of each ML model are optimized at each forward iteration. Final outputs highlight an optimized number of selected features (brain regions of interest) for each model with its accuracy. Furthermore, the toolbox can be executed in a parallel environment for efficient computation on a typical personal computer. With the optimized forward variable selection decoder (oFVSD) pipeline, we verified the effectiveness of decoding sex classification and age range regression on 1,113 structural magnetic resonance imaging (MRI) datasets. Compared to ML models without the FVS algorithm and with the Boruta algorithm as a variable selection counterpart, we demonstrate that the oFVSD significantly outperformed across all of the ML models over the counterpart models without FVS (approximately 0.20 increase in correlation coefficient, r, with regression models and 8% increase in classification models on average) and with Boruta variable selection algorithm (approximately 0.07 improvement in regression and 4% in classification models). Furthermore, we confirmed the use of parallel computation considerably reduced the computational burden for the high-dimensional MRI data. Altogether, the oFVSD toolbox efficiently and effectively improves the performance of both classification and regression ML models, providing a use case example on MRI datasets. With its flexibility, oFVSD has the potential for many other modalities in neuroimaging. This open-source and freely available Python package makes it a valuable toolbox for research communities seeking improved decoding accuracy.

An insula hierarchical network architecture for active interoceptive inference

In the brain, the insular cortex receives a vast amount of interoceptive information, ascending through deep brain structures, from multiple visceral organs. The unique hierarchical and modular architecture of the insula suggests specialization for processing interoceptive afferents. Yet, the biological significance of the insula's neuroanatomical architecture, in relation to deep brain structures, remains obscure. In this opinion piece, we propose the Insula Hierarchical Modular Adaptive Interoception Control (IMAC) model to suggest that insula modules (granular, dysgranular and agranular), forming parallel networks with the prefrontal cortex and striatum, are specialized to form higher order interoceptive representations. These interoceptive representations are recruited in a context-dependent manner to support habitual, model-based and exploratory control of visceral organs and physiological processes. We discuss how insula interoceptive representations may give rise to conscious feelings that best explain lower order deep brain interoceptive representations, and how the insula may serve to defend the body and mind against pathological depression.

Frequency drift in MR spectroscopy at 3T

PURPOSE

Heating of gradient coils and passive shim components is a common cause of instability in the B0 field, especially when gradient intensive sequences are used. The aim of the study was to set a benchmark for typical drift encountered during MR spectroscopy (MRS) to assess the need for real-time field-frequency locking on MRI scanners by comparing field drift data from a large number of sites.

METHOD

A standardized protocol was developed for 80 participating sites using 99 3T MR scanners from 3 major vendors. Phantom water signals were acquired before and after an EPI sequence. The protocol consisted of: minimal preparatory imaging; a short pre-fMRI PRESS; a ten-minute fMRI acquisition; and a long post-fMRI PRESS acquisition. Both pre- and post-fMRI PRESS were non-water suppressed. Real-time frequency stabilization/adjustment was switched off when appropriate. Sixty scanners repeated the protocol for a second dataset. In addition, a three-hour post-fMRI MRS acquisition was performed at one site to observe change of gradient temperature and drift rate. Spectral analysis was performed using MATLAB. Frequency drift in pre-fMRI PRESS data were compared with the first 5:20 minutes and the full 30:00 minutes of data after fMRI. Median (interquartile range) drifts were measured and showed in violin plot. Paired t-tests were performed to compare frequency drift pre- and post-fMRI. A simulated in vivo spectrum was generated using FID-A to visualize the effect of the observed frequency drifts. The simulated spectrum was convolved with the frequency trace for the most extreme cases. Impacts of frequency drifts on NAA and GABA were also simulated as a function of linear drift. Data from the repeated protocol were compared with the corresponding first dataset using Pearson's and intraclass correlation coefficients (ICC).

RESULTS

Of the data collected from 99 scanners, 4 were excluded due to various reasons. Thus, data from 95 scanners were ultimately analyzed. For the first 5:20 min (64 transients), median (interquartile range) drift was 0.44 (1.29) Hz before fMRI and 0.83 (1.29) Hz after. This increased to 3.15 (4.02) Hz for the full 30 min (360 transients) run. Average drift rates were 0.29 Hz/min before fMRI and 0.43 Hz/min after. Paired t-tests indicated that drift increased after fMRI, as expected (p < 0.05). Simulated spectra convolved with the frequency drift showed that the intensity of the NAA singlet was reduced by up to 26%, 44 % and 18% for GE, Philips and Siemens scanners after fMRI, respectively. ICCs indicated good agreement between datasets acquired on separate days. The single site long acquisition showed drift rate was reduced to 0.03 Hz/min approximately three hours after fMRI.

DISCUSSION

This study analyzed frequency drift data from 95 3T MRI scanners. Median levels of drift were relatively low (5-min average under 1 Hz), but the most extreme cases suffered from higher levels of drift. The extent of drift varied across scanners which both linear and nonlinear drifts were observed.

Association of the oxytocin receptor gene with attitudinal trust: role of amygdala volume

Previous studies have shown that genetic variations in rs53576, a common variant of the oxytocin receptor gene (OXTR) resulting from a single nucleotide polymorphism involving an adenine (A)/guanine (G) transition, are associated with attitudinal trust in men. However, the pathway from gene to behaviour has not been elucidated. We conducted the present study to determine whether amygdala volume mediates the association between OXTR rs53576 genotypes and attitudinal trust. Our results revealed that the left amygdala volume was significantly smaller in GG men than in AA and AG men, whereas it was significantly smaller in AA and AG women than in GG women. In addition, the left amygdala volume was negatively associated with attitudinal trust in men, whereas there was no such association in women. We also found a significant mediation effect of the left amygdala volume on the association between OXTR rs53576 genotypes and attitudinal trust in men. The results of our study suggest that the left amygdala volume plays a pivotal role in the association between OXTR rs53576 genotypes and attitudinal trust in men.

Behavioural Differences and Neural Substrates of Altruistic and Spiteful Punishment

Altruistic punishment following social norm violations promotes human cooperation. However, experimental evidence indicates that some forms of punishment are spiteful rather than altruistic. Using two types of punishment games and seven non-strategic games, we identified strong behavioural differences between altruistic and spiteful punishers. Altruistic punishers who rejected unfair offers in the ultimatum game and punished norm violators in the third-party punishment game behaved pro-socially in various non-strategic games. Spiteful punishers who rejected unfair offers in the ultimatum game but did not punish norm violators in the third-party punishment game behaved selfishly in non-strategic games. In addition, the left caudate nucleus was larger in spiteful punishers than in altruistic punishers. These findings are in contrast to the previous assumption that altruistic punishers derive pleasure from enforcement of fairness norms, and suggest that spiteful punishers derive pleasure from seeing the target experience negative consequences.