Neurodesk: an accessible, flexible and portable data analysis environment for reproducible neuroimaging.

Neuroimaging research requires purpose-built analysis software, which is challenging to install and may produce different results across computing environments. The community-oriented, open-source Neurodesk platform ( https://www.neurodesk.org/ ) harnesses a comprehensive and growing suite of neuroimaging software containers. Neurodesk includes a browser-accessible virtual desktop, command-line interface and computational notebook compatibility, allowing for accessible, flexible, portable and fully reproducible neuroimaging analysis on personal workstations, high-performance computers and the cloud.

Web-based processing of physiological noise in fMRI: addition of the PhysIO toolbox to CBRAIN.

Neuroimaging research requires sophisticated tools for analyzing complex data, but efficiently leveraging these tools can be a major challenge, especially on large datasets. CBRAIN is a web-based platform designed to simplify the use and accessibility of neuroimaging research tools for large-scale, collaborative studies. In this paper, we describe how CBRAIN's unique features and infrastructure were leveraged to integrate TAPAS PhysIO, an open-source MATLAB toolbox for physiological noise modeling in fMRI data. This case study highlights three key elements of CBRAIN's infrastructure that enable streamlined, multimodal tool integration: a user-friendly GUI, a Brain Imaging Data Structure (BIDS) data-entry schema, and convenient in-browser visualization of results. By incorporating PhysIO into CBRAIN, we achieved significant improvements in the speed, ease of use, and scalability of physiological preprocessing. Researchers now have access to a uniform and intuitive interface for analyzing data, which facilitates remote and collaborative evaluation of results. With these improvements, CBRAIN aims to become an essential open-science tool for integrative neuroimaging research, supporting FAIR principles and enabling efficient workflows for complex analysis pipelines.

Neurodesk: An accessible, flexible, and portable data analysis environment for reproducible neuroimaging.

Neuroimaging data analysis often requires purpose-built software, which can be challenging to install and may produce different results across computing environments. Beyond being a roadblock to neuroscientists, these issues of accessibility and portability can hamper the reproducibility of neuroimaging data analysis pipelines. Here, we introduce the Neurodesk platform, which harnesses software containers to support a comprehensive and growing suite of neuroimaging software (https://www.neurodesk.org/). Neurodesk includes a browser-accessible virtual desktop environment and a command line interface, mediating access to containerized neuroimaging software libraries on various computing platforms, including personal and high-performance computers, cloud computing and Jupyter Notebooks. This community-oriented, open-source platform enables a paradigm shift for neuroimaging data analysis, allowing for accessible, flexible, fully reproducible, and portable data analysis pipelines.

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

Mono-planar T-Hex: Speed and flexibility for high-resolution 3D imaging.

PURPOSE: The aim of this work is the reconciliation of high spatial and temporal resolution for MRI. For this purpose, a novel sampling strategy for 3D encoding is proposed, which provides flexible k-space segmentation along with uniform sampling density and benign filtering effects related to signal decay. METHODS: For time-critical MRI applications such as functional MRI (fMRI), 3D k-space is usually sampled by stacking together 2D trajectories such as echo planar imaging (EPI) or spiral readouts, where each shot covers one k-space plane. For very high temporal and medium to low spatial resolution, tilted hexagonal sampling (T-Hex) was recently proposed, which allows the acquisition of a larger k-space volume per excitation than can be covered with a planar readout. Here, T-Hex is described in a modified version where it instead acquires a smaller k-space volume per shot for use with medium temporal and high spatial resolution. RESULTS: Mono-planar T-Hex sampling provides flexibility in the choice of speed, signal-to-noise ratio (SNR), and contrast for rapid MRI acquisitions. For use with a conventional gradient system, it offers the greatest benefit in a regime of high in-plane resolution <1 mm. The sampling scheme is combined with spirals for high sampling speed as well as with more conventional EPI trajectories. CONCLUSION: Mono-planar T-Hex sampling combines fast 3D encoding with SNR efficiency and favorable depiction characteristics regarding noise amplification and filtering effects from T2∗ decay, thereby providing flexibility in the choice of imaging parameters. It is attractive both for high-resolution time series such as fMRI and for applications that require rapid anatomical imaging.

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.

Feasibility of spiral fMRI based on an LTI gradient model.

Spiral imaging is very well suited for functional MRI, however its use has been limited by the fact that artifacts caused by gradient imperfections and B0 inhomogeneity are more difficult to correct compared to EPI. Effective correction requires accurate knowledge of the traversed k-space trajectory. With the goal of making spiral fMRI more accessible, we have evaluated image reconstruction using trajectories predicted by the gradient impulse response function (GIRF), which can be determined in a one-time calibration step. GIRF-predicted reconstruction was tested for high-resolution (0.8 mm) fMRI at 7T. Image quality and functional results of the reconstructions using GIRF-prediction were compared to reconstructions using the nominal trajectory and concurrent field monitoring. The reconstructions using nominal spiral trajectories contain substantial artifacts and the activation maps contain misplaced activation. Image artifacts are substantially reduced when using the GIRF-predicted reconstruction, and the activation maps for the GIRF-predicted and monitored reconstructions largely overlap. The GIRF reconstruction provides a large increase in the spatial specificity of the activation compared to the nominal reconstruction. The GIRF-reconstruction generates image quality and fMRI results similar to using a concurrently monitored trajectory. The presented approach does not prolong or complicate the fMRI acquisition. Using GIRF-predicted trajectories has the potential to enable high-quality spiral fMRI in situations where concurrent trajectory monitoring is not available.

TAPAS: An Open-Source Software Package for Translational Neuromodeling and Computational Psychiatry.

Psychiatry faces fundamental challenges with regard to mechanistically guided differential diagnosis, as well as prediction of clinical trajectories and treatment response of individual patients. This has motivated the genesis of two closely intertwined fields: (i) Translational Neuromodeling (TN), which develops "computational assays" for inferring patient-specific disease processes from neuroimaging, electrophysiological, and behavioral data; and (ii) Computational Psychiatry (CP), with the goal of incorporating computational assays into clinical decision making in everyday practice. In order to serve as objective and reliable tools for clinical routine, computational assays require end-to-end pipelines from raw data (input) to clinically useful information (output). While these are yet to be established in clinical practice, individual components of this general end-to-end pipeline are being developed and made openly available for community use. In this paper, we present the Translational Algorithms for Psychiatry-Advancing Science (TAPAS) software package, an open-source collection of building blocks for computational assays in psychiatry. Collectively, the tools in TAPAS presently cover several important aspects of the desired end-to-end pipeline, including: (i) tailored experimental designs and optimization of measurement strategy prior to data acquisition, (ii) quality control during data acquisition, and (iii) artifact correction, statistical inference, and clinical application after data acquisition. Here, we review the different tools within TAPAS and illustrate how these may help provide a deeper understanding of neural and cognitive mechanisms of disease, with the ultimate goal of establishing automatized pipelines for predictions about individual patients. We hope that the openly available tools in TAPAS will contribute to the further development of TN/CP and facilitate the translation of advances in computational neuroscience into clinically relevant computational assays.

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.

Whole-brain estimates of directed connectivity for human connectomics.

Connectomics is essential for understanding large-scale brain networks but requires that individual connection estimates are neurobiologically interpretable. In particular, a principle of brain organization is that reciprocal connections between cortical areas are functionally asymmetric. This is a challenge for fMRI-based connectomics in humans where only undirected functional connectivity estimates are routinely available. By contrast, whole-brain estimates of effective (directed) connectivity are computationally challenging, and emerging methods require empirical validation. Here, using a motor task at 7T, we demonstrate that a novel generative model can infer known connectivity features in a whole-brain network (>200 regions, >40,000 connections) highly efficiently. Furthermore, graph-theoretical analyses of directed connectivity estimates identify functional roles of motor areas more accurately than undirected functional connectivity estimates. These results, which can be achieved in an entirely unsupervised manner, demonstrate the feasibility of inferring directed connections in whole-brain networks and open new avenues for human connectomics.

T-Hex: Tilted hexagonal grids for rapid 3D imaging.

PURPOSE: The purpose of this work is to devise and demonstrate an encoding strategy for 3D MRI that reconciles high speed with flexible segmentation, uniform k-space density, and benign T2∗ effects. METHODS: Fast sampling of a 3D k-space is typically accomplished by 2D readouts per shot using EPI trains or spiral readouts. Tilted hexagonal (T-Hex) sampling is a way of acquiring more k-space volume per excitation while maintaining uniform sampling density and a smooth T2∗ filter. The k-space volume covered per shot is controlled by the tilting angle. Image reconstruction is performed with a 3D extension of the iterative SENSE approach, incorporating actual field dynamics and static off-resonance. T-Hex imaging is compared with established 3D schemes in terms of speed and noise performance. RESULTS: Tilted hexagonal acquisition is found to achieve greater imaging speed than known alternatives, particularly in combination with spiral trajectories. The interplay of the proposed 3D trajectories, array detection, and off-resonance is successfully addressed by iterative inversion of the full signal model. Enhanced coverage per shot is of greatest utility for high speed in an intermediate resolution regime of 1 to 4 mm. T-Hex EPI combines the benefits of extended coverage per shot with increased robustness against off-resonance effects. CONCLUSION: Sampling of tilted hexagonal grids is a feasible means of gaining 3D imaging speed with near-optimal SNR efficiency and benign depiction properties. It is a particularly promising technique for time-resolved applications such as fMRI.

CG-SENSE revisited: Results from the first ISMRM reproducibility challenge.

PURPOSE: The aim of this work is to shed light on the issue of reproducibility in MR image reconstruction in the context of a challenge. Participants had to recreate the results of "Advances in sensitivity encoding with arbitrary k-space trajectories" by Pruessmann et al. METHODS: The task of the challenge was to reconstruct radially acquired multicoil k-space data (brain/heart) following the method in the original paper, reproducing its key figures. Results were compared to consolidated reference implementations created after the challenge, accounting for the two most common programming languages used in the submissions (Matlab/Python). RESULTS: Visually, differences between submissions were small. Pixel-wise differences originated from image orientation, assumed field-of-view, or resolution. The reference implementations were in good agreement, both visually and in terms of image similarity metrics. DISCUSSION AND CONCLUSION: While the description level of the published algorithm enabled participants to reproduce CG-SENSE in general, details of the implementation varied, for example, density compensation or Tikhonov regularization. Implicit assumptions about the data lead to further differences, emphasizing the importance of sufficient metadata accompanying open datasets. Defining reproducibility quantitatively turned out to be nontrivial for this image reconstruction challenge, in the absence of ground-truth results. Typical similarity measures like NMSE of SSIM were misled by image intensity scaling and outlier pixels. Thus, to facilitate reproducibility, researchers are encouraged to publish code and data alongside the original paper. Future methodological papers on MR image reconstruction might benefit from the consolidated reference implementations of CG-SENSE presented here, as a benchmark for methods comparison.

Cholinergic and dopaminergic effects on prediction error and uncertainty responses during sensory associative learning.

Navigating the physical world requires learning probabilistic associations between sensory events and their change in time (volatility). Bayesian accounts of this learning process rest on hierarchical prediction errors (PEs) that are weighted by estimates of uncertainty (or its inverse, precision). In a previous fMRI study we found that low-level precision-weighted PEs about visual outcomes (that update beliefs about associations) activated the putative dopaminergic midbrain; by contrast, precision-weighted PEs about cue-outcome associations (that update beliefs about volatility) activated the cholinergic basal forebrain. These findings suggested selective dopaminergic and cholinergic influences on precision-weighted PEs at different hierarchical levels. Here, we tested this hypothesis, repeating our fMRI study under pharmacological manipulations in healthy participants. Specifically, we performed two pharmacological fMRI studies with a between-subject double-blind placebo-controlled design: study 1 used antagonists of dopaminergic (amisulpride) and muscarinic (biperiden) receptors, study 2 used enhancing drugs of dopaminergic (levodopa) and cholinergic (galantamine) modulation. Pooled across all pharmacological conditions of study 1 and study 2, respectively, we found that low-level precision-weighted PEs activated the midbrain and high-level precision-weighted PEs the basal forebrain as in our previous study. However, we found pharmacological effects on brain activity associated with these computational quantities only when splitting the precision-weighted PEs into their PE and precision components: in a brainstem region putatively containing cholinergic (pedunculopontine and laterodorsal tegmental) nuclei, biperiden (compared to placebo) enhanced low-level PE responses and attenuated high-level PE activity, while amisulpride reduced high-level PE responses. Additionally, in the putative dopaminergic midbrain, galantamine compared to placebo enhanced low-level PE responses (in a body-weight dependent manner) and amisulpride enhanced high-level precision activity. Task behaviour was not affected by any of the drugs. These results do not support our hypothesis of a clear-cut dichotomy between different hierarchical inference levels and neurotransmitter systems, but suggest a more complex interaction between these neuromodulatory systems and hierarchical Bayesian quantities. However, our present results may have been affected by confounds inherent to pharmacological fMRI. We discuss these confounds and outline improved experimental tests for the future.

The relationship between resting-state functional connectivity, antidepressant discontinuation and depression relapse.

The risk of relapsing into depression after stopping antidepressants is high, but no established predictors exist. Resting-state functional magnetic resonance imaging (rsfMRI) measures may help predict relapse and identify the mechanisms by which relapses occur. rsfMRI data were acquired from healthy controls and from patients with remitted major depressive disorder on antidepressants. Patients were assessed a second time either before or after discontinuation of the antidepressant, and followed up for six months to assess relapse. A seed-based functional connectivity analysis was conducted focusing on the left subgenual anterior cingulate cortex and left posterior cingulate cortex. Seeds in the amygdala and dorsolateral prefrontal cortex were explored. 44 healthy controls (age: 33.8 (10.5), 73% female) and 84 patients (age: 34.23 (10.8), 80% female) were included in the analysis. 29 patients went on to relapse and 38 remained well. The seed-based analysis showed that discontinuation resulted in an increased functional connectivity between the right dorsolateral prefrontal cortex and the parietal cortex in non-relapsers. In an exploratory analysis, this functional connectivity predicted relapse risk with a balanced accuracy of 0.86. Further seed-based analyses, however, failed to reveal differences in functional connectivity between patients and controls, between relapsers and non-relapsers before discontinuation and changes due to discontinuation independent of relapse. In conclusion, changes in the connectivity between the dorsolateral prefrontal cortex and the posterior default mode network were associated with and predictive of relapse after open-label antidepressant discontinuation. This finding requires replication in a larger dataset.

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

Neural arbitration between social and individual learning systems.

Decision making requires integrating knowledge gathered from personal experiences with advice from others. The neural underpinnings of the process of arbitrating between information sources has not been fully elucidated. In this study, we formalized arbitration as the relative precision of predictions, afforded by each learning system, using hierarchical Bayesian modeling. In a probabilistic learning task, participants predicted the outcome of a lottery using recommendations from a more informed advisor and/or self-sampled outcomes. Decision confidence, as measured by the number of points participants wagered on their predictions, varied with our definition of arbitration as a ratio of precisions. Functional neuroimaging demonstrated that arbitration signals were independent of decision confidence and involved modality-specific brain regions. Arbitrating in favor of self-gathered information activated the dorsolateral prefrontal cortex and the midbrain, whereas arbitrating in favor of social information engaged the ventromedial prefrontal cortex and the amygdala. These findings indicate that relative precision captures arbitration between social and individual learning systems at both behavioral and neural levels.

A Reconfigurable Platform for Magnetic Resonance Data Acquisition and Processing.

Developments in magnetic resonance imaging (MRI) in the last decades show a trend towards a growing number of array coils and an increasing use of a wide variety of sensors. Associated cabling and safety issues have been addressed by moving data acquisition closer to the coil. However, with the increasing number of radio-frequency (RF) channels and trend towards higher acquisition duty-cycles, the data amount is growing, which poses challenges for throughput and data handling. As it is becoming a limitation, early compression and preprocessing is becoming ever more important. Additionally, sensors deliver diverse data, which require distinct and often low-latency processing for run-time updates of scanner operation. To address these challenges, we propose the transition to reconfigurable hardware with an application tailored assembly of interfaces and real-time processing resources. We present an integrated solution based on a system-on-chip (SoC), which offers sufficient throughput and hardware-based parallel processing power for very challenging applications. It is equipped with fiber-optical modules serving as versatile interfaces for modular systems with in-field operation. We demonstrate the utility of the platform on the example of concurrent imaging and field sensing with hardware-based coil compression and trajectory extraction. The preprocessed data are then used in expanded encoding model based image reconstruction of single-shot and segmented spirals as used in time-series and anatomical imaging respectively.