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.

Deep brain stimulation has state-dependent effects on motor connectivity in Parkinson's disease.

Subthalamic nucleus deep brain stimulation is an effective treatment for advanced Parkinson's disease; however, its therapeutic mechanism is unclear. Previous modelling of functional MRI data has suggested that deep brain stimulation has modulatory effects on a number of basal ganglia pathways. This work uses an enhanced data collection protocol to collect rare functional MRI data in patients with subthalamic nucleus deep brain stimulation. Eleven patients with Parkinson's disease and subthalamic nucleus deep brain stimulation underwent functional MRI at rest and during a movement task; once with active deep brain stimulation, and once with deep brain stimulation switched off. Dynamic causal modelling and Bayesian model selection were first used to compare a series of plausible biophysical models of the cortico-basal ganglia circuit that could explain the functional MRI activity at rest in an attempt to reproduce and extend the findings from our previous work. General linear modelling of the movement task functional MRI data revealed deep brain stimulation-associated signal increases in the primary motor and cerebellar cortices. Given the significance of the cerebellum in voluntary movement, we then built a more complete model of the motor system by including cerebellar-basal ganglia interactions, and compared the modulatory effects deep brain stimulation had on different circuit components during the movement task and again using the resting state data. Consistent with previous results from our independent cohort, model comparison found that the rest data were best explained by deep brain stimulation-induced increased (effective) connectivity of the cortico-striatal, thalamo-cortical and direct pathway and reduced coupling of subthalamic nucleus afferent and efferent connections. No changes in cerebellar connectivity were identified at rest. In contrast, during the movement task, there was functional recruitment of subcortical-cerebellar pathways, which were additionally modulated by deep brain stimulation, as well as modulation of local (intrinsic) cortical and cerebellar circuits. This work provides in vivo evidence for the modulatory effects of subthalamic nucleus deep brain stimulation on effective connectivity within the cortico-basal ganglia loops at rest, as well as further modulations in the cortico-cerebellar motor system during voluntary movement. We propose that deep brain stimulation has both behaviour-independent effects on basal ganglia connectivity, as well as behaviour-dependent modulatory effects.

Analysis of family-wise error rates in statistical parametric mapping using random field theory.

This technical report revisits the analysis of family-wise error rates in statistical parametric mapping-using random field theory-reported in (Eklund et al. []: arXiv 1511.01863). Contrary to the understandable spin that these sorts of analyses attract, a review of their results suggests that they endorse the use of parametric assumptions-and random field theory-in the analysis of functional neuroimaging data. We briefly rehearse the advantages parametric analyses offer over nonparametric alternatives and then unpack the implications of (Eklund et al. []: arXiv 1511.01863) for parametric procedures. Hum Brain Mapp, 40:2052-2054, 2019. © 2017 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc.

BIDS apps: Improving ease of use, accessibility, and reproducibility of neuroimaging data analysis methods.

The rate of progress in human neurosciences is limited by the inability to easily apply a wide range of analysis methods to the plethora of different datasets acquired in labs around the world. In this work, we introduce a framework for creating, testing, versioning and archiving portable applications for analyzing neuroimaging data organized and described in compliance with the Brain Imaging Data Structure (BIDS). The portability of these applications (BIDS Apps) is achieved by using container technologies that encapsulate all binary and other dependencies in one convenient package. BIDS Apps run on all three major operating systems with no need for complex setup and configuration and thanks to the comprehensiveness of the BIDS standard they require little manual user input. Previous containerized data processing solutions were limited to single user environments and not compatible with most multi-tenant High Performance Computing systems. BIDS Apps overcome this limitation by taking advantage of the Singularity container technology. As a proof of concept, this work is accompanied by 22 ready to use BIDS Apps, packaging a diverse set of commonly used neuroimaging algorithms.

Physiological basis of vascular autocalibration (VasA): Comparison to hypercapnia calibration methods.

PURPOSE: The statistical power of functional MRI (fMRI) group studies is significantly hampered by high intersubject spatial and magnitude variance. We recently presented a vascular autocalibration method (VasA) to account for vascularization differences between subjects and hence improve the sensitivity in group studies. Here, we validate the novel calibration method by means of direct comparisons of VasA with more established measures of baseline venous blood volume (and indirectly vascular reactivity), the M-value. METHODS: Seven healthy volunteers participated in two 7 T (T) fMRI experiments to compare M-values with VasA estimates: (i) a hypercapnia experiment to estimate voxelwise M-value maps, and (ii) an fMRI experiment using visual stimulation to estimate voxelwise VasA maps. RESULTS: We show that VasA and M-value calibration maps show the same spatial profile, providing strong evidence that VasA is driven by local variations in vascular reactivity as reflected in the M-value. CONCLUSION: The agreement of vascular reactivity maps obtained with VasA when compared with M-value maps confirms empirically the hypothesis that the VasA method is an adequate tool to account for variations in fMRI response amplitudes caused by vascular reactivity differences in healthy volunteers. VasA can therefore directly account for them and increase the statistical power of group studies. The VasA toolbox is available as a statistical parametric mapping (SPM) toolbox, facilitating its general application. Magn Reson Med, 2016. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine.

Evaluation of 2D multiband EPI imaging for high-resolution, whole-brain, task-based fMRI studies at 3T: Sensitivity and slice leakage artifacts.

Functional magnetic resonance imaging (fMRI) studies that require high-resolution whole-brain coverage have long scan times that are primarily driven by the large number of thin slices acquired. Two-dimensional multiband echo-planar imaging (EPI) sequences accelerate the data acquisition along the slice direction and therefore represent an attractive approach to such studies by improving the temporal resolution without sacrificing spatial resolution. In this work, a 2D multiband EPI sequence was optimized for 1.5mm isotropic whole-brain acquisitions at 3T with 10 healthy volunteers imaged while performing simultaneous visual and motor tasks. The performance of the sequence was evaluated in terms of BOLD sensitivity and false-positive activation at multiband (MB) factors of 1, 2, 4, and 6, combined with in-plane GRAPPA acceleration of 2× (GRAPPA 2), and the two reconstruction approaches of Slice-GRAPPA and Split Slice-GRAPPA. Sensitivity results demonstrate significant gains in temporal signal-to-noise ratio (tSNR) and t-score statistics for MB 2, 4, and 6 compared to MB 1. The MB factor for optimal sensitivity varied depending on anatomical location and reconstruction method. When using Slice-GRAPPA reconstruction, evidence of false-positive activation due to signal leakage between simultaneously excited slices was seen in one instance, 35 instances, and 70 instances over the ten volunteers for the respective accelerations of MB 2×GRAPPA 2, MB 4×GRAPPA 2, and MB 6×GRAPPA 2. The use of Split Slice-GRAPPA reconstruction suppressed the prevalence of false positives significantly, to 1 instance, 5 instances, and 5 instances for the same respective acceleration factors. Imaging protocols using an acceleration factor of MB 2×GRAPPA 2 can be confidently used for high-resolution whole-brain imaging to improve BOLD sensitivity with very low probability for false-positive activation due to slice leakage. Imaging protocols using higher acceleration factors (MB 3 or MB 4×GRAPPA 2) can likely provide even greater gains in sensitivity but should be carefully optimized to minimize the possibility of false activations.

Vascular autorescaling of fMRI (VasA fMRI) improves sensitivity of population studies: A pilot study.

The blood oxygenation level-dependent (BOLD) signal is widely used for functional magnetic resonance imaging (fMRI) of brain function in health and disease. The statistical power of fMRI group studies is significantly hampered by high inter-subject variance due to differences in baseline vascular physiology. Several methods have been proposed to account for physiological vascularization differences between subjects and hence improve the sensitivity in group studies. However, these methods require the acquisition of additional reference scans (such as a full resting-state fMRI session or ASL-based calibrated BOLD). We present a vascular autorescaling (VasA) method, which does not require any additional reference scans. VasA is based on the observation that slow oscillations (<0.1Hz) in arterial blood CO2 levels occur naturally due to changes in respiration patterns. These oscillations yield fMRI signal changes whose amplitudes reflect the blood oxygenation levels and underlying local vascularization and vascular responsivity. VasA estimates proxies of the amplitude of these CO2-driven oscillations directly from the residuals of task-related fMRI data without the need for reference scans. The estimates are used to scale the amplitude of task-related fMRI responses, to account for vascular differences. The VasA maps compared well to cerebrovascular reactivity (CVR) maps and cerebral blood volume maps based on vascular space occupancy (VASO) measurements in four volunteers, speaking to the physiological vascular basis of VasA. VasA was validated in a wide variety of tasks in 138 volunteers. VasA increased t-scores by up to 30% in specific brain areas such as the visual cortex. The number of activated voxels was increased by up to 200% in brain areas such as the orbital frontal cortex while still controlling the nominal false-positive rate. VasA fMRI outperformed previously proposed rescaling approaches based on resting-state fMRI data and can be readily applied to any task-related fMRI data set, even retrospectively.

Resting state functional MRI in Parkinson's disease: the impact of deep brain stimulation on 'effective' connectivity.

Depleted of dopamine, the dynamics of the parkinsonian brain impact on both 'action' and 'resting' motor behaviour. Deep brain stimulation has become an established means of managing these symptoms, although its mechanisms of action remain unclear. Non-invasive characterizations of induced brain responses, and the effective connectivity underlying them, generally appeals to dynamic causal modelling of neuroimaging data. When the brain is at rest, however, this sort of characterization has been limited to correlations (functional connectivity). In this work, we model the 'effective' connectivity underlying low frequency blood oxygen level-dependent fluctuations in the resting Parkinsonian motor network-disclosing the distributed effects of deep brain stimulation on cortico-subcortical connections. Specifically, we show that subthalamic nucleus deep brain stimulation modulates all the major components of the motor cortico-striato-thalamo-cortical loop, including the cortico-striatal, thalamo-cortical, direct and indirect basal ganglia pathways, and the hyperdirect subthalamic nucleus projections. The strength of effective subthalamic nucleus afferents and efferents were reduced by stimulation, whereas cortico-striatal, thalamo-cortical and direct pathways were strengthened. Remarkably, regression analysis revealed that the hyperdirect, direct, and basal ganglia afferents to the subthalamic nucleus predicted clinical status and therapeutic response to deep brain stimulation; however, suppression of the sensitivity of the subthalamic nucleus to its hyperdirect afferents by deep brain stimulation may subvert the clinical efficacy of deep brain stimulation. Our findings highlight the distributed effects of stimulation on the resting motor network and provide a framework for analysing effective connectivity in resting state functional MRI with strong a priori hypotheses.

The Past, Present, and Future of the Brain Imaging Data Structure (BIDS).

The Brain Imaging Data Structure (BIDS) is a community-driven standard for the organization of data and metadata from a growing range of neuroscience modalities. This paper is meant as a history of how the standard has developed and grown over time. We outline the principles behind the project, the mechanisms by which it has been extended, and some of the challenges being addressed as it evolves. We also discuss the lessons learned through the project, with the aim of enabling researchers in other domains to learn from the success of BIDS.

Convolution models for induced electromagnetic responses.

In Kilner et al. [Kilner, J.M., Kiebel, S.J., Friston, K.J., 2005. Applications of random field theory to electrophysiology. Neurosci. Lett. 374, 174-178.] we described a fairly general analysis of induced responses-in electromagnetic brain signals-using the summary statistic approach and statistical parametric mapping. This involves localising induced responses-in peristimulus time and frequency-by testing for effects in time-frequency images that summarise the response of each subject to each trial type. Conventionally, these time-frequency summaries are estimated using post-hoc averaging of epoched data. However, post-hoc averaging of this sort fails when the induced responses overlap or when there are multiple response components that have variable timing within each trial (for example stimulus and response components associated with different reaction times). In these situations, it is advantageous to estimate response components using a convolution model of the sort that is standard in the analysis of fMRI time series. In this paper, we describe one such approach, based upon ordinary least squares deconvolution of induced responses to input functions encoding the onset of different components within each trial. There are a number of fundamental advantages to this approach: for example; (i) one can disambiguate induced responses to stimulus onsets and variably timed responses; (ii) one can test for the modulation of induced responses-over peristimulus time and frequency-by parametric experimental factors and (iii) one can gracefully handle confounds-such as slow drifts in power-by including them in the model. In what follows, we consider optimal forms for convolution models of induced responses, in terms of impulse response basis function sets and illustrate the utility of deconvolution estimators using simulated and real MEG data.

Set-level threshold-free tests on the intrinsic volumes of SPMs.

Conventionally, set-level inference on statistical parametric maps (SPMs) is based on the topological features of an excursion set above some threshold-for example, the number of clusters or Euler characteristic. The expected Euler characteristic-under the null hypothesis-can be predicted from an intrinsic measure or volume of the SPM, such as the resel counts or the Lipschitz-Killing curvatures (LKC). We propose a new approach that performs a null hypothesis omnibus test on an SPM, by testing whether its intrinsic volume (described by LKC coefficients) is different from the volume of the underlying residual fields: intuitively, whether the number of peaks in the statistical field (testing for signal) and the residual fields (noise) are consistent or not. Crucially, this new test requires no arbitrary feature-defining threshold but is nevertheless sensitive to distributed or spatially extended patterns. We show the similarities between our approach and conventional topological inference-in terms of false positive rate control and sensitivity to treatment effects-in two and three dimensional simulations. The test consistently improves on classical approaches for moderate (>20) degrees of freedom. We also demonstrate the application to real data and illustrate the comparison of the expected and observed Euler characteristics over the complete threshold range.

Genome-wide epigenetic and mRNA-expression profiling followed by CRISPR/Cas9-mediated gene-disruptions corroborate the MIR141/MIR200C-ZEB1/ZEB2-FGFR1 axis in acquired EMT-associated EGFR TKI-resistance in NSCLC cells.

Background: Epithelial-mesenchymal-transition (EMT) is an epigenetic-based mechanism contributing to the acquired treatment resistance against receptor tyrosine kinase inhibitors (TKIs) in non-small cell lung cancer (NSCLC) cells harboring epidermal growth factor receptor (EGFR)-mutations. Delineating the exact epigenetic and gene-expression alterations in EMT-associated EGFR TKI-resistance (EMT-E-TKI-R) is vital for improved diagnosis and treatment of NSCLC patients. Methods: We characterized genome-wide changes in mRNA-expression, DNA-methylation and the histone-modification H3K36me3 in EGFR-mutated NSCLC HCC827 cells in result of acquired EMT-E-TKI-R. CRISPR/Cas9 was used to functional examine key findings from the omics analyses. Results: Acquired EMT-E-TKI-R was analyzed with three omics approaches. RNA-sequencing identified 2,233 and 1,972 up- and down-regulated genes, respectively, and among these were established EMT-markers. DNA-methylation EPIC array analyses identified 14,163 and 7,999 hyper- and hypo-methylated, respectively, differential methylated positions of which several were present in EMT-markers. Finally, H3K36me3 chromatin immunoprecipitation (ChIP)-sequencing detected 2,873 and 3,836 genes with enrichment and depletion, respectively, and among these were established EMT-markers. Correlation analyses showed that EMT-E-TKI-R mRNA-expression changes correlated better with H3K36me3 changes than with DNA-methylation changes. Moreover, the omics data supported the involvement of the MIR141/MIR200C-ZEB1/ZEB2-FGFR1 signaling axis for acquired EMT-E-TKI-R. CRISPR/Cas9-mediated analyses corroborated the importance of ZEB1 in acquired EMT-E-TKI-R, MIR200C and MIR141 to be in an EMT-E-TKI-R-associated auto-regulatory loop with ZEB1, and FGFR1 to mediate cell survival in EMT-E-TKI-R. Conclusions: The current study describes the synchronous genome-wide changes in mRNA-expression, DNA-methylation, and H3K36me3 in NSCLC EMT-E-TKI-R. The omics approaches revealed potential novel diagnostic markers and treatment targets. Besides, the study consolidates the functional impact of the MIR141/MIR200C-ZEB1/ZEB2-FGFR1-signaling axis in NSCLC EMT-E-TKI-R.

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.

The problem of low variance voxels in statistical parametric mapping; a new hat avoids a 'haircut'.

Statistical parametric mapping (SPM) locates significant clusters based on a ratio of signal to noise (a 'contrast' of the parameters divided by its standard error) meaning that very low noise regions, for example outside the brain, can attain artefactually high statistical values. Similarly, the commonly applied preprocessing step of Gaussian spatial smoothing can shift the peak statistical significance away from the peak of the contrast and towards regions of lower variance. These problems have previously been identified in positron emission tomography (PET) (Reimold et al., 2006) and voxel-based morphometry (VBM) (Acosta-Cabronero et al., 2008), but can also appear in functional magnetic resonance imaging (fMRI) studies. Additionally, for source-reconstructed magneto- and electro-encephalography (M/EEG), the problems are particularly severe because sparsity-favouring priors constrain meaningfully large signal and variance to a small set of compactly supported regions within the brain. (Acosta-Cabronero et al., 2008) suggested adding noise to background voxels (the 'haircut'), effectively increasing their noise variance, but at the cost of contaminating neighbouring regions with the added noise once smoothed. Following theory and simulations, we propose to modify--directly and solely--the noise variance estimate, and investigate this solution on real imaging data from a range of modalities.

Dopamine-dependent prediction errors underpin reward-seeking behaviour in humans.

Theories of instrumental learning are centred on understanding how success and failure are used to improve future decisions. These theories highlight a central role for reward prediction errors in updating the values associated with available actions. In animals, substantial evidence indicates that the neurotransmitter dopamine might have a key function in this type of learning, through its ability to modulate cortico-striatal synaptic efficacy. However, no direct evidence links dopamine, striatal activity and behavioural choice in humans. Here we show that, during instrumental learning, the magnitude of reward prediction error expressed in the striatum is modulated by the administration of drugs enhancing (3,4-dihydroxy-L-phenylalanine; L-DOPA) or reducing (haloperidol) dopaminergic function. Accordingly, subjects treated with L-DOPA have a greater propensity to choose the most rewarding action relative to subjects treated with haloperidol. Furthermore, incorporating the magnitude of the prediction errors into a standard action-value learning algorithm accurately reproduced subjects' behavioural choices under the different drug conditions. We conclude that dopamine-dependent modulation of striatal activity can account for how the human brain uses reward prediction errors to improve future decisions.

Dynamic causal modelling of COVID-19 and its mitigations.

This technical report describes the dynamic causal modelling of mitigated epidemiological outcomes during the COVID-9 coronavirus outbreak in 2020. Dynamic causal modelling is a form of complex system modelling, which uses 'real world' timeseries to estimate the parameters of an underlying state space model using variational Bayesian procedures. Its key contribution-in an epidemiological setting-is to embed conventional models within a larger model of sociobehavioural responses-in a way that allows for (relatively assumption-free) forecasting. One advantage of using variational Bayes is that one can progressively optimise the model via Bayesian model selection: generally, the most likely models become more expressive as more data becomes available. This report summarises the model (on 6-Nov-20), eight months after the inception of dynamic causal modelling for COVID-19. This model-and its subsequent updates-is used to provide nowcasts and forecasts of latent behavioural and epidemiological variables as an open science resource. The current report describes the underlying model structure and the rationale for the variational procedures that underwrite Bayesian model selection.

Predicting language lateralization from gray matter.

It has long been predicted that the degree to which language is lateralized to the left or right hemisphere might be reflected in the underlying brain anatomy. We investigated this relationship on a voxel-by-voxel basis across the whole brain using structural and functional magnetic resonance images from 86 healthy participants. Structural images were converted to gray matter probability images, and language activation was assessed during naming and semantic decision. All images were spatially normalized to the same symmetrical template, and lateralization images were generated by subtracting right from left hemisphere signal at each voxel. We show that the degree to which language was left or right lateralized was positively correlated with the degree to which gray matter density was lateralized. Post hoc analyses revealed a general relationship between gray matter probability and blood oxygenation level-dependent signal. This is the first demonstration that structural brain scans can be used to predict language lateralization on a voxel-by-voxel basis in the normal healthy brain.

A parametric empirical Bayesian framework for fMRI-constrained MEG/EEG source reconstruction.

We describe an asymmetric approach to fMRI and MEG/EEG fusion in which fMRI data are treated as empirical priors on electromagnetic sources, such that their influence depends on the MEG/EEG data, by virtue of maximizing the model evidence. This is important if the causes of the MEG/EEG signals differ from those of the fMRI signal. Furthermore, each suprathreshold fMRI cluster is treated as a separate prior, which is important if fMRI data reflect neural activity arising at different times within the EEG/MEG data. We present methodological considerations when mapping from a 3D fMRI Statistical Parametric Map to a 2D cortical surface and thence to the covariance components used within our Parametric Empirical Bayesian framework. Our previous introduction of a canonical (inverse-normalized) cortical mesh also allows deployment of fMRI priors that live in a template space; for example, from a group analysis of different individuals. We evaluate the ensuing scheme with MEG and EEG data recorded simultaneously from 12 participants, using the same face-processing paradigm under which independent fMRI data were obtained. Because the fMRI priors become part of the generative model, we use the model evidence to compare (i) multiple versus single, (ii) valid versus invalid, (iii) binary versus continuous, and (iv) variance versus covariance fMRI priors. For these data, multiple, valid, binary, and variance fMRI priors proved best for a standard Minimum Norm inversion. Interestingly, however, inversion using Multiple Sparse Priors benefited little from additional fMRI priors, suggesting that they already provide a sufficiently flexible generative model.

Dynamic causal modelling of COVID-19.

This technical report describes a dynamic causal model of the spread of coronavirus through a population. The model is based upon ensemble or population dynamics that generate outcomes, like new cases and deaths over time. The purpose of this model is to quantify the uncertainty that attends predictions of relevant outcomes. By assuming suitable conditional dependencies, one can model the effects of interventions (e.g., social distancing) and differences among populations (e.g., herd immunity) to predict what might happen in different circumstances. Technically, this model leverages state-of-the-art variational (Bayesian) model inversion and comparison procedures, originally developed to characterise the responses of neuronal ensembles to perturbations. Here, this modelling is applied to epidemiological populations-to illustrate the kind of inferences that are supported and how the model per se can be optimised given timeseries data. Although the purpose of this paper is to describe a modelling protocol, the results illustrate some interesting perspectives on the current pandemic; for example, the nonlinear effects of herd immunity that speak to a self-organised mitigation process.

Second waves, social distancing, and the spread of COVID-19 across the USA.

We recently described a dynamic causal model of a COVID-19 outbreak within a single region. Here, we combine several of these (epidemic) models to create a (pandemic) model of viral spread among regions. Our focus is on a second wave of new cases that may result from loss of immunity-and the exchange of people between regions-and how mortality rates can be ameliorated under different strategic responses. In particular, we consider hard or soft social distancing strategies predicated on national (Federal) or regional (State) estimates of the prevalence of infection in the population. The modelling is demonstrated using timeseries of new cases and deaths from the United States to estimate the parameters of a factorial (compartmental) epidemiological model of each State and, crucially, coupling between States. Using Bayesian model reduction, we identify the effective connectivity between States that best explains the initial phases of the outbreak in the United States. Using the ensuing posterior parameter estimates, we then evaluate the likely outcomes of different policies in terms of mortality, working days lost due to lockdown and demands upon critical care. The provisional results of this modelling suggest that social distancing and loss of immunity are the two key factors that underwrite a return to endemic equilibrium.