A computationally informed distinction of interoception and exteroception.

While interoception is of major neuroscientific interest, its precise definition and delineation from exteroception continue to be debated. Here, we propose a functional distinction between interoception and exteroception based on computational concepts of sensor-effector loops. Under this view, the classification of sensory inputs as serving interoception or exteroception depends on the sensor-effector loop they feed into, for the control of either bodily (physiological and biochemical) or environmental states. We explain the utility of this perspective by examining the perception of skin temperature, one of the most challenging cases for distinguishing between interoception and exteroception. Specifically, we propose conceptualising thermoception as inference about the thermal state of the body (including the skin), which is directly coupled to thermoregulatory processes. This functional view emphasises the coupling to regulation (control) as a defining property of perception (inference) and connects the definition of interoception to contemporary computational theories of brain-body interactions.

An observational treatment study of metacognition in anxious-depression.

Prior studies have found metacognitive biases are linked to a transdiagnostic dimension of anxious-depression, manifesting as reduced confidence in performance. However, previous work has been cross-sectional and so it is unclear if under-confidence is a trait-like marker of anxious-depression vulnerability, or if it resolves when anxious-depression improves. Data were collected as part of a large-scale transdiagnostic, four-week observational study of individuals initiating internet-based cognitive behavioural therapy (iCBT) or antidepressant medication. Self-reported clinical questionnaires and perceptual task performance were gathered to assess anxious-depression and metacognitive bias at baseline and 4-week follow-up. Primary analyses were conducted for individuals who received iCBT (n=649), with comparisons between smaller samples that received antidepressant medication (n=82) and a control group receiving no intervention (n=88). Prior to receiving treatment, anxious-depression severity was associated with under-confidence in performance in the iCBT arm, replicating previous work. From baseline to follow-up, levels of anxious-depression were significantly reduced, and this was accompanied by a significant increase in metacognitive confidence in the iCBT arm (ß=0.17, SE=0.02, p<0.001). These changes were correlated (r(647)=-0.12, p=0.002); those with the greatest reductions in anxious-depression levels had the largest increase in confidence. While the three-way interaction effect of group and time on confidence was not significant (F(2, 1632)=0.60, p=0.550), confidence increased in the antidepressant group (ß=0.31, SE = 0.08, p<0.001), but not among controls (ß=0.11, SE = 0.07, p=0.103). Metacognitive biases in anxious-depression are state-dependent; when symptoms improve with treatment, so does confidence in performance. Our results suggest this is not specific to the type of intervention.

Interoceptive and metacognitive facets of fatigue in multiple sclerosis.

Numerous disorders are characterised by fatigue as a highly disabling symptom. Fatigue plays a particularly important clinical role in multiple sclerosis (MS) where it exerts a profound impact on quality of life. Recent concepts of fatigue grounded in computational theories of brain-body interactions emphasise the role of interoception and metacognition in the pathogenesis of fatigue. So far, however, for MS, empirical data on interoception and metacognition are scarce. This study examined interoception and (exteroceptive) metacognition in a sample of 71 persons with a diagnosis of MS. Interoception was assessed by prespecified subscales of a standard questionnaire (Multidimensional Assessment of Interoceptive Awareness [MAIA]), while metacognition was investigated with computational models of choice and confidence data from a visual discrimination paradigm. Additionally, autonomic function was examined by several physiological measurements. Several hypotheses were tested based on a preregistered analysis plan. In brief, we found the predicted association of interoceptive awareness with fatigue (but not with exteroceptive metacognition) and an association of autonomic function with exteroceptive metacognition (but not with fatigue). Furthermore, machine learning (elastic net regression) showed that individual fatigue scores could be predicted out-of-sample from our measurements, with questionnaire-based measures of interoceptive awareness and sleep quality as key predictors. Our results support theoretical concepts of interoception as an important factor for fatigue and demonstrate the general feasibility of predicting individual levels of fatigue from simple questionnaire-based measures of interoception and sleep.

Predicting future depressive episodes from resting-state fMRI with generative embedding.

After a first episode of major depressive disorder (MDD), there is substantial risk for a long-term remitting-relapsing course. Prevention and early interventions are thus critically important. Various studies have examined the feasibility of detecting at-risk individuals based on out-of-sample predictions about the future occurrence of depression. However, functional magnetic resonance imaging (fMRI) has received very little attention for this purpose so far. Here, we explored the utility of generative models (i.e. different dynamic causal models, DCMs) as well as functional connectivity (FC) for predicting future episodes of depression in never-depressed adults, using a large dataset (N = 906) of task-free ("resting state") fMRI data from the UK Biobank (UKB). Connectivity analyses were conducted using timeseries from pre-computed spatially independent components of different dimensionalities. Over a three-year period, 50% of selected participants showed indications of at least one depressive episode, while the other 50% did not. Using nested cross-validation for training and a held-out test set (80/20 split), we systematically examined the combination of 8 connectivity feature sets and 17 classifiers. We found that a generative embedding procedure based on combining regression DCM (rDCM) with a support vector machine (SVM) enabled the best predictions, both on the training set (0.63 accuracy, 0.66 area under the curve, AUC) and the test set (0.62 accuracy, 0.64 AUC; p < 0.001). However, on the test set, rDCM was only slightly superior to predictions based on FC (0.59 accuracy, 0.61 AUC). Interpreting model predictions based on SHAP (SHapley Additive exPlanations) values suggested that the most predictive connections were widely distributed and not confined to specific networks. Overall, our analyses suggest (i) ways of improving future fMRI-based generative embedding approaches for the early detection of individuals at-risk for depression and that (ii) achieving accuracies of clinical utility may require combination of fMRI with other data modalities.

Using standardized patients for undergraduate clinical skills training in an introductory course to psychiatry.

BACKGROUND: The goal of this study was to assess the value and acceptance of Standardized or Simulated Patients (SPs) for training clinically inexperienced undergraduate medical students in psychiatric history taking, psychopathological assessment, and communication with psychiatric patients. METHODS: As part of a newly developed introductory course to psychiatry, pairs of 3rd year medical students conducted psychiatric assessments of SPs, including history and psychopathological state, under the supervision of a clinical lecturer. Prior to the assessment, students attended introductory lectures to communication in psychiatry and psychopathology but were clinically inexperienced. After the interview, the students' summary of their findings was discussed with other students and the lecturer. Students, lecturers, and actors were invited to a survey after the course. Questions for the students included self-reports about perceived learning success and authenticity of the interviews. RESULTS: 41 students, 6 actors and 8 lecturers completed the survey (response rates of 48%, 50%, and 100%, respectively). The survey results indicated that, despite their lack of clinical experience, students learned how to conduct a psychiatric interview, communicate in a non-judgmental and empathetic manner, take a psychiatric history and perform a psychopathological examination. SPs were perceived as authentic. The survey results suggested that this setting allowed for an enjoyable, non-distressful and motivating learning experience within a restricted time frame of just two afternoons. CONCLUSION: The results indicated that the SP approach presented is useful for teaching clinical skills in psychiatry to students with limited previous clinical experience and knowledge of psychiatry. We argue that SPs can be used to teach practical psychiatric skills already during an early phase of the curriculum. Limitations of our study include a limited sample size, a temporal gap between the course and the survey, reliance on self-reports, and lack of comparison to alternative interventions.

Auditory mismatch responses are differentially sensitive to changes in muscarinic acetylcholine versus dopamine receptor function.

The auditory mismatch negativity (MMN) has been proposed as a biomarker of NMDA receptor (NMDAR) dysfunction in schizophrenia. Such dysfunction may be caused by aberrant interactions of different neuromodulators with NMDARs, which could explain clinical heterogeneity among patients. In two studies (N = 81 each), we used a double-blind placebo-controlled between-subject design to systematically test whether auditory mismatch responses under varying levels of environmental stability are sensitive to diminishing and enhancing cholinergic vs. dopaminergic function. We found a significant drug × mismatch interaction: while the muscarinic acetylcholine receptor antagonist biperiden delayed and topographically shifted mismatch responses, particularly during high stability, this effect could not be detected for amisulpride, a dopamine D2/D3 receptor antagonist. Neither galantamine nor levodopa, which elevate acetylcholine and dopamine levels, respectively, exerted significant effects on MMN. This differential MMN sensitivity to muscarinic versus dopaminergic receptor function may prove useful for developing tests that predict individual treatment responses in schizophrenia.

Advances in spiral fMRI: A high-resolution dataset.

We present data collected for the research article "Advances in Spiral fMRI: A High-resolution Study with Single-shot Acquisition" (Kasper et al. 2022). All data was acquired on a 7T ultra-high field MR system (Philips Achieva), equipped with a concurrent magnetic field monitoring setup based on 16 NMR probes. For task-based fMRI, a visual quarterfield stimulation paradigm was employed, alongside physiological monitoring via peripheral recordings. This data collection contains different datasets pertaining to different purposes: (1) Measured magnetic field dynamics (k0, spiral k-space trajectories, 2nd order spherical harmonics, concomitant fields) during ultra-high field fMRI sessions from six subjects, as well as concurrent temperature curves of the gradient coil, to explore MR system and subject-induced variability of field fluctuations and assess the impact of potential correction methods. (2) MR Raw Data, i.e., coil and concurrent encoding magnetic field (trajectory) data, of a single subject, as well as nominal spiral gradient waveforms, precomputed B0 and coil sensitivity maps, to enable testing of alternative image reconstruction approaches for spiral fMRI data. (3) Reconstructed image time series of the same subject alongside behavioral and physiological logfiles, to reproduce the fMRI preprocessing and analysis, as well as figures presented in the research article related to this article, using the published analysis code repository. All data is provided in standardized formats for the respective research area. In particular, ISMRMRD (HDF5) is used to store raw coil data and spiral trajectories, as well as measured field dynamics. Likewise, the NIfTI format is used for all imaging data (anatomical MRI and spiral fMRI, B0 and coil sensitivity maps).

Advances in spiral fMRI: A high-resolution study with single-shot acquisition.

Spiral fMRI has been put forward as a viable alternative to rectilinear echo-planar imaging, in particular due to its enhanced average k-space speed and thus high acquisition efficiency. This renders spirals attractive for contemporary fMRI applications that require high spatiotemporal resolution, such as laminar or columnar fMRI. However, in practice, spiral fMRI is typically hampered by its reduced robustness and ensuing blurring artifacts, which arise from imperfections in both static and dynamic magnetic fields. Recently, these limitations have been overcome by the concerted application of an expanded signal model that accounts for such field imperfections, and its inversion by iterative image reconstruction. In the challenging ultra-high field environment of 7 Tesla, where field inhomogeneity effects are aggravated, both multi-shot and single-shot 2D spiral imaging at sub-millimeter resolution was demonstrated with high depiction quality and anatomical congruency. In this work, we further these advances towards a time series application of spiral readouts, namely, single-shot spiral BOLD fMRI at 0.8 mm in-plane resolution. We demonstrate that high-resolution spiral fMRI at 7 T is not only feasible, but delivers both excellent image quality, BOLD sensitivity, and spatial specificity of the activation maps, with little artifactual blurring. Furthermore, we show the versatility of the approach with a combined in/out spiral readout at a more typical resolution (1.5 mm), where the high acquisition efficiency allows to acquire two images per shot for improved sensitivity by echo combination.

Optogenetic activation of striatal D1R and D2R cells differentially engages downstream connected areas beyond the basal ganglia.

The basal ganglia (BG) are a group of subcortical nuclei responsible for motor and executive function. Central to BG function are striatal cells expressing D1 (D1R) and D2 (D2R) dopamine receptors. D1R and D2R cells are considered functional antagonists that facilitate voluntary movements and inhibit competing motor patterns, respectively. However, whether they maintain a uniform function across the striatum and what influence they exert outside the BG is unclear. Here, we address these questions by combining optogenetic activation of D1R and D2R cells in the mouse ventrolateral caudoputamen with fMRI. Striatal D1R/D2R stimulation evokes distinct activity within the BG-thalamocortical network and differentially engages cerebellar and prefrontal regions. Computational modeling of effective connectivity confirms that changes in D1R/D2R output drive functional relationships between these regions. Our results suggest a complex functional organization of striatal D1R/D2R cells and hint toward an interconnected fronto-BG-cerebellar network modulated by striatal D1R and D2R cells.

Inference on homeostatic belief precision.

Interoception and homeostatic/allostatic control are intertwined branches of closed-loop brain-body interactions (BBI). Given their importance in mental and psychosomatic disorders, establishing computational assays of BBI represents a clinically important but methodologically challenging endeavor. This technical note presents a novel approach, derived from a generic computational model of homeostatic/allostatic control that underpins (meta)cognitive theories of affective and psychosomatic disorders. This model views homeostatic setpoints as probability distributions ("homeostatic beliefs") whose parameters determine regulatory efforts and change dynamically under allostatic predictions. In particular, changes in homeostatic belief precision, triggered by anticipated threats to homeostasis, are thought to alter cerebral regulation of bodily states. Here, we present statistical procedures for inferring homeostatic belief precision from measured bodily states and/or regulatory (action) signals. We analyze the inference problem, derive two alternative estimators of homeostatic belief precision, and apply our method to simulated data. Our proposed approach may prove useful for assessing BBI in individual subjects.

Technical note: A fast and robust integrator of delay differential equations in DCM for electrophysiological data.

Dynamic causal models (DCMs) of electrophysiological data allow, in principle, for inference on hidden, bulk synaptic function in neural circuits. The directed influences between the neuronal elements of modeled circuits are subject to delays due to the finite transmission speed of axonal connections. Ordinary differential equations are therefore not adequate to capture the ensuing circuit dynamics, and delay differential equations (DDEs) are required instead. Previous work has illustrated that the integration of DDEs in DCMs benefits from sophisticated integration schemes in order to ensure rigorous parameter estimation and correct model identification. However, integration schemes that have been proposed for DCMs either emphasize speed (at the possible expense of accuracy) or robustness (but with computational costs that are problematic in practice). In this technical note, we propose an alternative integration scheme that overcomes these shortcomings and offers high computational efficiency while correctly preserving the nature of delayed effects. This integration scheme is available as open-source code in the Translational Algorithms for Psychiatry-Advancing Science (TAPAS) toolbox and can be easily integrated into existing software (SPM) for the analysis of DCMs for electrophysiological data. While this paper focuses on its application to the convolution-based formalism of DCMs, the new integration scheme can be equally applied to more advanced formulations of DCMs (e.g. conductance based models). Our method provides a new option for electrophysiological DCMs that offers the speed required for scientific projects, but also the accuracy required for rigorous translational applications, e.g. in computational psychiatry.

Model-based prediction of muscarinic receptor function from auditory mismatch negativity responses.

Drugs affecting neuromodulation, for example by dopamine or acetylcholine, take centre stage among therapeutic strategies in psychiatry. These neuromodulators can change both neuronal gain and synaptic plasticity and therefore affect electrophysiological measures. An important goal for clinical diagnostics is to exploit this effect in the reverse direction, i.e., to infer the status of specific neuromodulatory systems from electrophysiological measures. In this study, we provide proof-of-concept that the functional status of cholinergic (specifically muscarinic) receptors can be inferred from electrophysiological data using generative (dynamic causal) models. To this end, we used epidural EEG recordings over two auditory cortical regions during a mismatch negativity (MMN) paradigm in rats. All animals were treated, across sessions, with muscarinic receptor agonists and antagonists at different doses. Together with a placebo condition, this resulted in five levels of muscarinic receptor status. Using a dynamic causal model - embodying a small network of coupled cortical microcircuits - we estimated synaptic parameters and their change across pharmacological conditions. The ensuing parameter estimates associated with (the neuromodulation of) synaptic efficacy showed both graded muscarinic effects and predictive validity between agonistic and antagonistic pharmacological conditions. This finding illustrates the potential utility of generative models of electrophysiological data as computational assays of muscarinic function. In application to EEG data of patients from heterogeneous spectrum diseases, e.g. schizophrenia, such models might help identify subgroups of patients that respond differentially to cholinergic treatments. SIGNIFICANCE STATEMENT: In psychiatry, the vast majority of pharmacological treatments affect actions of neuromodulatory transmitters, e.g. dopamine or acetylcholine. As treatment is largely trial-and-error based, one of the goals for computational psychiatry is to construct mathematical models that can serve as "computational assays" and infer the status of specific neuromodulatory systems in individual patients. This translational neuromodeling strategy has great promise for electrophysiological data in particular but requires careful validation. The present study demonstrates that the functional status of cholinergic (muscarinic) receptors can be inferred from electrophysiological data using dynamic causal models of neural circuits. While accuracy needs to be enhanced and our results must be replicated in larger samples, our current results provide proof-of-concept for computational assays of muscarinic function using EEG.

A Hilbert-based method for processing respiratory timeseries.

In this technical note, we introduce a new method for estimating changes in respiratory volume per unit time (RVT) from respiratory bellows recordings. By using techniques from the electrophysiological literature, in particular the Hilbert transform, we show how we can better characterise breathing rhythms, with the goal of improving physiological noise correction in functional magnetic resonance imaging (fMRI). Specifically, our approach leads to a representation with higher time resolution and better captures atypical breathing events than current peak-based RVT estimators. Finally, we demonstrate that this leads to an increase in the amount of respiration-related variance removed from fMRI data when used as part of a typical preprocessing pipeline. Our implementation is publicly available as part of the PhysIO package, which is distributed as part of the open-source TAPAS toolbox (https://translationalneuromodeling.org/tapas).

Hemodynamic modeling of long-term aspirin effects on blood oxygenated level dependent responses at 7 Tesla in patients at cardiovascular risk.

Aspirin is considered a potential confound for functional magnetic resonance imaging (fMRI) studies. This is because aspirin affects the synthesis of prostaglandin, a vasoactive mediator centrally involved in neurovascular coupling, a process underlying blood oxygenated level dependent (BOLD) responses. Aspirin-induced changes in BOLD signal are a potential confound for fMRI studies of at-risk individuals or patients (e.g. with cardiovascular conditions or stroke) who receive low-dose aspirin prophylactically and are compared to healthy controls without aspirin. To examine the severity of this potential confound, we combined high field (7 Tesla) MRI during a simple hand movement task with a biophysically informed hemodynamic model. We compared elderly individuals receiving aspirin for primary or secondary prophylactic purposes versus age-matched volunteers without aspirin medication, testing for putative differences in BOLD responses. Specifically, we fitted hemodynamic models to BOLD responses from 14 regions activated by the task and examined whether model parameter estimates were significantly altered by aspirin. While our analyses indicate that hemodynamics differed across regions, consistent with the known regional variability of BOLD responses, we neither found a significant main effect of aspirin (i.e., an average effect across brain regions) nor an expected drug × region interaction. While our sample size is not sufficiently large to rule out small-to-medium global effects of aspirin, we had adequate statistical power for detecting the expected interaction. Altogether, our analysis suggests that patients with cardiovascular risk receiving low-dose aspirin for primary or secondary prophylactic purposes do not show strongly altered BOLD signals when compared to healthy controls without aspirin.

Hierarchical Bayesian models of social inference for probing persecutory delusional ideation.

While persecutory delusions (PDs) have been linked to fallacies of reasoning and social inference, computational characterizations of delusional tendencies are rare. Here, we examined 151 individuals from the general population on opposite ends of the PD spectrum (Paranoia Checklist [PCL]). Participants made trial-wise predictions in a probabilistic lottery, guided by advice from a more informed human and a nonsocial cue. Additionally, 2 frames differentially emphasized causes of invalid advice: (a) the adviser's possible intentions (dispositional frame) or (b) the rules of the game (situational frame). We applied computational modeling to examine possible reasons for group differences in behavior. Comparing different models, we found that a hierarchical Bayesian model (hierarchical Gaussian filter) explained participants' responses better than other learning models. Model parameters determining participants' belief updates about the adviser's fidelity and the contribution of prior beliefs about fidelity to trial-wise decisions, respectively, showed significant Group × Frame interactions: High PCL scorers held more rigid beliefs about the adviser's fidelity across both experimental frames and relied less on advice in situational frames than low scorers. These results suggest that PD tendencies are associated with rigid beliefs and prevent adaptive use of social information in "safe" contexts. This supports previous proposals of a link between PD and aberrant social inference. (PsycInfo Database Record (c) 2020 APA, all rights reserved).

Remote, Automated, and MRI-Compatible Administration of Interoceptive Inspiratory Resistive Loading.

Research on how humans perceive sensory inputs from their bodies ("interoception") has been rapidly gaining momentum, with interest across a host of disciplines from physiology through to psychiatry. However, studying interoceptive processes is not without significant challenges, and many methods utilized to access internal states have been largely devoted to capturing and relating naturally occurring variations in interoceptive signals (such as heartbeats) to measures of how the brain processes these signals. An alternative procedure involves the controlled perturbation of specific interoceptive axes. This is challenging because it requires non-invasive interventions that can be repeated many times within a subject and that are potent but safe. Here we present an effective methodology for instigating these perturbations within the breathing domain. We describe a custom-built circuitry that is capable of delivering inspiratory resistive loads automatically and precisely. Importantly, our approach is compatible with magnetic resonance imaging (MRI) environments, allowing for the administration of complicated experimental designs in neuroimaging as increasingly required within developing fields such as computational psychiatry/psychosomatics. We describe the experimental setup for both the control and monitoring of the inspiratory resistive loads, and demonstrate its possible utilities within different study designs. This methodology represents an important step forward from the previously utilized, manually controlled resistive loading setups, which present significant experimental burdens with prolonged and/or complicated sequences of breathing stimuli.

Volatility Estimates Increase Choice Switching and Relate to Prefrontal Activity in Schizophrenia.

BACKGROUND: Reward-based decision making is impaired in patients with schizophrenia (PSZ), as reflected by increased choice switching. The underlying cognitive and motivational processes as well as associated neural signatures remain unknown. Reinforcement learning and hierarchical Bayesian learning account for choice switching in different ways. We hypothesized that enhanced choice switching, as seen in PSZ during reward-based decision making, relates to higher-order beliefs about environmental volatility, and we examined the associated neural activity. METHODS: In total, 46 medicated PSZ and 43 healthy control subjects performed a reward-based decision-making task requiring flexible responses to changing action-outcome contingencies during functional magnetic resonance imaging. Detailed computational modeling of choice data was performed, including reinforcement learning and the hierarchical Gaussian filter. Trajectories of learning from computational modeling informed the analysis of functional magnetic resonance imaging data. RESULTS: A 3-level hierarchical Gaussian filter accounted best for the observed choice data. This model revealed a heightened initial belief about environmental volatility and a stronger influence of volatility on lower-level learning of action-outcome contingencies in PSZ as compared with healthy control subjects. This was replicated in an independent sample of nonmedicated PSZ. Beliefs about environmental volatility were reflected by higher activity in dorsolateral prefrontal cortex of PSZ as compared with healthy control subjects. CONCLUSIONS: Our study suggests that PSZ inferred the environment as overly volatile, which may explain increased choice switching. In PSZ, activity in dorsolateral prefrontal cortex was more strongly related to beliefs about environmental volatility. Our computational phenotyping approach may provide useful information to dissect clinical heterogeneity and could improve prediction of outcome.

Inflexible social inference in individuals with subclinical persecutory delusional tendencies.

It has been suspected that abnormalities in social inference (e.g., learning others' intentions) play a key role in the formation of persecutory delusions (PD). In this study, we examined the association between subclinical PD and social inference, testing the prediction that proneness to PD is related to altered social inference and beliefs about others' intentions. We included 151 participants scoring on opposite ends of Freeman's Paranoia Checklist (PCL). The participants performed a probabilistic advice-taking task with a dynamically changing social context (volatility) under one of two experimental frames. These frames differentially emphasised possible reasons behind unhelpful advice: (i) the adviser's possible intentions (dispositional frame) or (ii) the rules of the game (situational frame). Our design was thus 2 × 2 factorial (high vs. low delusional tendencies, dispositional vs. situational frame). We found significant group-by-frame interactions, indicating that in the situational frame high PCL scorers took advice less into account than low scorers. Additionally, high PCL scorers believed more frequently that incorrect advice was delivered intentionally and that such misleading behaviour was directed towards them personally. Overall, our results suggest that social inference in individuals with subclinical PD tendencies is shaped by negative prior beliefs about the intentions of others and is thus less sensitive to the attributional framing of adviser-related information. These findings may help future attempts of identifying individuals at risk for developing psychosis and understanding persecutory delusions in psychosis.