Long-term soil warming changes the profile of primary metabolites in fine roots of Norway spruce in a temperate montane forest.

Climate warming poses major threats to temperate forests, but the response of tree root metabolism has largely remained unclear. We examined the impact of long-term soil warming (>14 years, +4°C) on the fine root metabolome across three seasons for 2 years in an old spruce forest, using a liquid chromatography-mass spectrometry platform for primary metabolite analysis. A total of 44 primary metabolites were identified in roots (19 amino acids, 12 organic acids and 13 sugars). Warming increased the concentration of total amino acids and of total sugars by 15% and 21%, respectively, but not organic acids. We found that soil warming and sampling date, along with their interaction, directly influenced the primary metabolite profiles. Specifically, in warming plots, concentrations of arginine, glycine, lysine, threonine, tryptophan, mannose, ribose, fructose, glucose and oxaloacetic acid increased by 51.4%, 19.9%, 21.5%, 19.3%, 22.1%, 23.0%, 38.0%, 40.7%, 19.8% and 16.7%, respectively. Rather than being driven by single compounds, changes in metabolite profiles reflected a general up- or downregulation of most metabolic pathway network. This emphasises the importance of metabolomics approaches in investigating root metabolic pathways and understanding the effects of climate change on tree root metabolism.

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.

Does long-term soil warming affect microbial element limitation? A test by short-term assays of microbial growth responses to labile C, N and P additions.

Increasing global temperatures have been reported to accelerate soil carbon (C) cycling, but also to promote nitrogen (N) and phosphorus (P) dynamics in terrestrial ecosystems. However, warming can differentially affect ecosystem C, N and P dynamics, potentially intensifying elemental imbalances between soil resources, plants and soil microorganisms. Here, we investigated the effect of long-term soil warming on microbial resource limitation, based on measurements of microbial growth (18 O incorporation into DNA) and respiration after C, N and P amendments. Soil samples were taken from two soil depths (0-10, 10-20 cm) in control and warmed (>14 years warming, +4°C) plots in the Achenkirch soil warming experiment. Soils were amended with combinations of glucose-C, inorganic/organic N and inorganic/organic P in a full factorial design, followed by incubation at their respective mean field temperatures for 24 h. Soil microbes were generally C-limited, exhibiting 1.8-fold to 8.8-fold increases in microbial growth upon C addition. Warming consistently caused soil microorganisms to shift from being predominately C limited to become C-P co-limited. This P limitation possibly was due to increased abiotic P immobilization in warmed soils. Microbes further showed stronger growth stimulation under combined glucose and inorganic nutrient amendments compared to organic nutrient additions. This may be related to a prolonged lag phase in organic N (glucosamine) mineralization and utilization compared to glucose. Soil respiration strongly positively responded to all kinds of glucose-C amendments, while responses of microbial growth were less pronounced in many of these treatments. This highlights that respiration-though easy and cheap to measure-is not a good substitute of growth when assessing microbial element limitation. Overall, we demonstrate a significant shift in microbial element limitation in warmed soils, from C to C-P co-limitation, with strong repercussions on the linkage between soil C, N and P cycles under long-term warming.

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.

Long-term soil warming alters fine root dynamics and morphology, and their ectomycorrhizal fungal community in a temperate forest soil.

Climate warming is predicted to affect temperate forests severely, but the response of fine roots, key to plant nutrition, water uptake, soil carbon, and nutrient cycling is unclear. Understanding how fine roots will respond to increasing temperature is a prerequisite for predicting the functioning of forests in a warmer climate. We studied the response of fine roots and their ectomycorrhizal (EcM) fungal and root-associated bacterial communities to soil warming by 4°C in a mixed spruce-beech forest in the Austrian Limestone Alps after 8 and 14 years of soil warming, respectively. Fine root biomass (FRB) and fine root production were 17% and 128% higher in the warmed plots, respectively, after 14 years. The increase in FRB (13%) was not significant after 8 years of treatment, whereas specific root length, specific root area, and root tip density were significantly higher in warmed plots at both sampling occasions. Soil warming did not affect EcM exploration types and diversity, but changed their community composition, with an increase in the relative abundance of Cenoccocum at 0-10 cm soil depth, a drought-stress-tolerant genus, and an increase in short- and long-distance exploration types like Sebacina and Boletus at 10-20 cm soil depth. Warming increased the root-associated bacterial diversity but did not affect their community composition. Soil warming did not affect nutrient concentrations of fine roots, though we found indications of limited soil phosphorus (P) and potassium (K) availability. Our findings suggest that, in the studied ecosystem, global warming could persistently increase soil carbon inputs due to accelerated fine root growth and turnover, and could simultaneously alter fine root morphology and EcM fungal community composition toward improved nutrient foraging.

Conductance-based dynamic causal modeling: A mathematical review of its application to cross-power spectral densities.

Dynamic Causal Modeling (DCM) is a Bayesian framework for inferring on hidden (latent) neuronal states, based on measurements of brain activity. Since its introduction in 2003 for functional magnetic resonance imaging data, DCM has been extended to electrophysiological data, and several variants have been developed. Their biophysically motivated formulations make these models promising candidates for providing a mechanistic understanding of human brain dynamics, both in health and disease. However, due to their complexity and reliance on concepts from several fields, fully understanding the mathematical and conceptual basis behind certain variants of DCM can be challenging. At the same time, a solid theoretical knowledge of the models is crucial to avoid pitfalls in the application of these models and interpretation of their results. In this paper, we focus on one of the most advanced formulations of DCM, i.e. conductance-based DCM for cross-spectral densities, whose components are described across multiple technical papers. The aim of the present article is to provide an accessible exposition of the mathematical background, together with an illustration of the model's behavior. To this end, we include step-by-step derivations of the model equations, point to important aspects in the software implementation of those models, and use simulations to provide an intuitive understanding of the type of responses that can be generated and the role that specific parameters play in the model. Furthermore, all code utilized for our simulations is made publicly available alongside the manuscript to allow readers an easy hands-on experience with conductance-based DCM.

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.

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).

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.

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.

Regression dynamic causal modeling for resting-state fMRI.

"Resting-state" functional magnetic resonance imaging (rs-fMRI) is widely used to study brain connectivity. So far, researchers have been restricted to measures of functional connectivity that are computationally efficient but undirected, or to effective connectivity estimates that are directed but limited to small networks. Here, we show that a method recently developed for task-fMRI-regression dynamic causal modeling (rDCM)-extends to rs-fMRI and offers both directional estimates and scalability to whole-brain networks. First, simulations demonstrate that rDCM faithfully recovers parameter values over a wide range of signal-to-noise ratios and repetition times. Second, we test construct validity of rDCM in relation to an established model of effective connectivity, spectral DCM. Using rs-fMRI data from nearly 200 healthy participants, rDCM produces biologically plausible results consistent with estimates by spectral DCM. Importantly, rDCM is computationally highly efficient, reconstructing whole-brain networks (>200 areas) within minutes on standard hardware. This opens promising new avenues for connectomics.

Timing of repetition suppression of event-related potentials to unattended objects.

Current theories of object perception emphasize the automatic nature of perceptual inference. Repetition suppression (RS), the successive decrease of brain responses to repeated stimuli, is thought to reflect the optimization of perceptual inference through neural plasticity. While functional imaging studies revealed brain regions that show suppressed responses to the repeated presentation of an object, little is known about the intra-trial time course of repetition effects to everyday objects. Here, we used event-related potentials (ERPs) to task-irrelevant line-drawn objects, while participants engaged in a distractor task. We quantified changes in ERPs over repetitions using three general linear models that modeled RS by an exponential, linear, or categorical "change detection" function in each subject. Our aim was to select the model with highest evidence and determine the within-trial time-course and scalp distribution of repetition effects using that model. Model comparison revealed the superiority of the exponential model indicating that repetition effects are observable for trials beyond the first repetition. Model parameter estimates revealed a sequence of RS effects in three time windows (86-140, 322-360, and 400-446 ms) and with occipital, temporoparietal, and frontotemporal distribution, respectively. An interval of repetition enhancement (RE) was also observed (320-340 ms) over occipitotemporal sensors. Our results show that automatic processing of task-irrelevant objects involves multiple intervals of RS with distinct scalp topographies. These sequential intervals of RS and RE might reflect the short-term plasticity required for optimization of perceptual inference and the associated changes in prediction errors and predictions, respectively, over stimulus repetitions during automatic object processing.

Visual Mismatch and Predictive Coding: A Computational Single-Trial ERP Study.

Predictive coding (PC) posits that the brain uses a generative model to infer the environmental causes of its sensory data and uses precision-weighted prediction errors (pwPEs) to continuously update this model. While supported by much circumstantial evidence, experimental tests grounded in formal trial-by-trial predictions are rare. One partial exception is event-related potential (ERP) studies of the auditory mismatch negativity (MMN), where computational models have found signatures of pwPEs and related model-updating processes. Here, we tested this hypothesis in the visual domain, examining possible links between visual mismatch responses and pwPEs. We used a novel visual "roving standard" paradigm to elicit mismatch responses in humans (of both sexes) by unexpected changes in either color or emotional expression of faces. Using a hierarchical Bayesian model, we simulated pwPE trajectories of a Bayes-optimal observer and used these to conduct a comprehensive trial-by-trial analysis across the time × sensor space. We found significant modulation of brain activity by both color and emotion pwPEs. The scalp distribution and timing of these single-trial pwPE responses were in agreement with visual mismatch responses obtained by traditional averaging and subtraction (deviant-minus-standard) approaches. Finally, we compared the Bayesian model to a more classical change model of MMN. Model comparison revealed that trial-wise pwPEs explained the observed mismatch responses better than categorical change detection. Our results suggest that visual mismatch responses reflect trial-wise pwPEs, as postulated by PC. These findings go beyond classical ERP analyses of visual mismatch and illustrate the utility of computational analyses for studying automatic perceptual processes.SIGNIFICANCE STATEMENT Human perception is thought to rely on a predictive model of the environment that is updated via precision-weighted prediction errors (pwPEs) when events violate expectations. This "predictive coding" view is supported by studies of the auditory mismatch negativity brain potential. However, it is less well known whether visual perception of mismatch relies on similar processes. Here we combined computational modeling and electroencephalography to test whether visual mismatch responses reflected trial-by-trial pwPEs. Applying a Bayesian model to series of face stimuli that violated expectations about color or emotional expression, we found significant modulation of brain activity by both color and emotion pwPEs. A categorical change detection model performed less convincingly. Our findings support the predictive coding interpretation of visual mismatch responses.

Feature-specific prediction errors for visual mismatch.

Predictive coding (PC) theory posits that our brain employs a predictive model of the environment to infer the causes of its sensory inputs. A fundamental but untested prediction of this theory is that the same stimulus should elicit distinct precision weighted prediction errors (pwPEs) when different (feature-specific) predictions are violated, even in the absence of attention. Here, we tested this hypothesis using functional magnetic resonance imaging (fMRI) and a multi-feature roving visual mismatch paradigm where rare changes in either color (red, green), or emotional expression (happy, fearful) of faces elicited pwPE responses in human participants. Using a computational model of learning and inference, we simulated pwPE and prediction trajectories of a Bayes-optimal observer and used these to analyze changes in blood oxygen level dependent (BOLD) responses to changes in color and emotional expression of faces while participants engaged in a distractor task. Controlling for visual attention by eye-tracking, we found pwPE responses to unexpected color changes in the fusiform gyrus. Conversely, unexpected changes of facial emotions elicited pwPE responses in cortico-thalamo-cerebellar structures associated with emotion and theory of mind processing. Predictions pertaining to emotions activated fusiform, occipital and temporal areas. Our results are consistent with a general role of PC across perception, from low-level to complex and socially relevant object features, and suggest that monitoring of the social environment occurs continuously and automatically, even in the absence of attention.

Dynamic causal modelling revisited.

This paper revisits the dynamic causal modelling of fMRI timeseries by replacing the usual (Taylor) approximation to neuronal dynamics with a neural mass model of the canonical microcircuit. This provides a generative or dynamic causal model of laminar specific responses that can generate haemodynamic and electrophysiological measurements. In principle, this allows the fusion of haemodynamic and (event related or induced) electrophysiological responses. Furthermore, it enables Bayesian model comparison of competing hypotheses about physiologically plausible synaptic effects; for example, does attentional modulation act on superficial or deep pyramidal cells - or both? In this technical note, we describe the resulting dynamic causal model and provide an illustrative application to the attention to visual motion dataset used in previous papers. Our focus here is on how to answer long-standing questions in fMRI; for example, do haemodynamic responses reflect extrinsic (afferent) input from distant cortical regions, or do they reflect intrinsic (recurrent) neuronal activity? To what extent do inhibitory interneurons contribute to neurovascular coupling? What is the relationship between haemodynamic responses and the frequency of induced neuronal activity? This paper does not pretend to answer these questions; rather it shows how they can be addressed using neural mass models of fMRI timeseries.

Switch costs in inhibitory control and voluntary behaviour: A computational study of the antisaccade task.

An integral aspect of human cognition is the ability to inhibit stimulus-driven, habitual responses, in favour of complex, voluntary actions. In addition, humans can also alternate between different tasks. This comes at the cost of degraded performance when compared to repeating the same task, a phenomenon called the "task-switch cost." While task switching and inhibitory control have been studied extensively, the interaction between them has received relatively little attention. Here, we used the SERIA model, a computational model of antisaccade behaviour, to draw a bridge between them. We investigated task switching in two versions of the mixed antisaccade task, in which participants are cued to saccade either in the same or in the opposite direction to a peripheral stimulus. SERIA revealed that stopping a habitual action leads to increased inhibitory control that persists onto the next trial, independently of the upcoming trial type. Moreover, switching between tasks induces slower and less accurate voluntary responses compared to repeat trials. However, this only occurs when participants lack the time to prepare the correct response. Altogether, SERIA demonstrates that there is a reconfiguration cost associated with switching between voluntary actions. In addition, the enhanced inhibition that follows antisaccade but not prosaccade trials explains asymmetric switch costs. In conclusion, SERIA offers a novel model of task switching that unifies previous theoretical accounts by distinguishing between inhibitory control and voluntary action generation and could help explain similar phenomena in paradigms beyond the antisaccade task.

Inhibition failures and late errors in the antisaccade task: influence of cue delay.

In the antisaccade task participants are required to saccade in the opposite direction of a peripheral visual cue (PVC). This paradigm is often used to investigate inhibition of reflexive responses as well as voluntary response generation. However, it is not clear to what extent different versions of this task probe the same underlying processes. Here, we explored with the Stochastic Early Reaction, Inhibition, and late Action (SERIA) model how the delay between task cue and PVC affects reaction time (RT) and error rate (ER) when pro- and antisaccade trials are randomly interleaved. Specifically, we contrasted a condition in which the task cue was presented before the PVC with a condition in which the PVC served also as task cue. Summary statistics indicate that ERs and RTs are reduced and contextual effects largely removed when the task is signaled before the PVC appears. The SERIA model accounts for RT and ER in both conditions and better so than other candidate models. Modeling demonstrates that voluntary pro- and antisaccades are frequent in both conditions. Moreover, early task cue presentation results in better control of reflexive saccades, leading to fewer fast antisaccade errors and more rapid correct prosaccades. Finally, high-latency errors are shown to be prevalent in both conditions. In summary, SERIA provides an explanation for the differences in the delayed and nondelayed antisaccade task. NEW & NOTEWORTHY In this article, we use a computational model to study the mixed antisaccade task. We contrast two conditions in which the task cue is presented either before or concurrently with the saccadic target. Modeling provides a highly accurate account of participants' behavior and demonstrates that a significant number of prosaccades are voluntary actions. Moreover, we provide a detailed quantitative analysis of the types of error that occur in pro- and antisaccade trials.

The Stochastic Early Reaction, Inhibition, and late Action (SERIA) model for antisaccades.

The antisaccade task is a classic paradigm used to study the voluntary control of eye movements. It requires participants to suppress a reactive eye movement to a visual target and to concurrently initiate a saccade in the opposite direction. Although several models have been proposed to explain error rates and reaction times in this task, no formal model comparison has yet been performed. Here, we describe a Bayesian modeling approach to the antisaccade task that allows us to formally compare different models on the basis of their evidence. First, we provide a formal likelihood function of actions (pro- and antisaccades) and reaction times based on previously published models. Second, we introduce the Stochastic Early Reaction, Inhibition, and late Action model (SERIA), a novel model postulating two different mechanisms that interact in the antisaccade task: an early GO/NO-GO race decision process and a late GO/GO decision process. Third, we apply these models to a data set from an experiment with three mixed blocks of pro- and antisaccade trials. Bayesian model comparison demonstrates that the SERIA model explains the data better than competing models that do not incorporate a late decision process. Moreover, we show that the early decision process postulated by the SERIA model is, to a large extent, insensitive to the cue presented in a single trial. Finally, we use parameter estimates to demonstrate that changes in reaction time and error rate due to the probability of a trial type (pro- or antisaccade) are best explained by faster or slower inhibition and the probability of generating late voluntary prosaccades.

A hemodynamic model for layered BOLD signals.

High-resolution blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) at the sub-millimeter scale has become feasible with recent advances in MR technology. In principle, this would enable the study of layered cortical circuits, one of the fundaments of cortical computation. However, the spatial layout of cortical blood supply may become an important confound at such high resolution. In particular, venous blood draining back to the cortical surface perpendicularly to the layered structure is expected to influence the measured responses in different layers. Here, we present an extension of a hemodynamic model commonly used for analyzing fMRI data (in dynamic causal models or biophysical network models) that accounts for such blood draining effects by coupling local hemodynamics across layers. We illustrate the properties of the model and its inversion by a series of simulations and show that it successfully captures layered fMRI data obtained during a simple visual experiment. We conclude that for future studies of the dynamics of layered neuronal circuits with high-resolution fMRI, it will be pivotal to include effects of blood draining, particularly when trying to infer on the layer-specific connections in cortex--a theme of key relevance for brain disorders like schizophrenia and for theories of brain function such as predictive coding.

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.