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.

Three-dimensional EPI with shot-selective CAIPIRIHANA for rapid high-resolution quantitative susceptibility mapping at 3 T.

PURPOSE: QSM provides insight into healthy brain aging and neuropathologies such as multiple sclerosis (MS), traumatic brain injuries, brain tumors, and neurodegenerative diseases. Phase data for QSM are usually acquired from 3D gradient-echo (3D GRE) scans with long acquisition times that are detrimental to patient comfort and susceptible to patient motion. This is particularly true for scans requiring whole-brain coverage and submillimeter resolutions. In this work, we use a multishot 3D echo plannar imaging (3D EPI) sequence with shot-selective 2D CAIPIRIHANA to acquire high-resolution, whole-brain data for QSM with minimal distortion and blurring. METHODS: To test clinical viability, the 3D EPI sequence was used to image a cohort of MS patients at 1-mm isotropic resolution at 3 T. Additionally, 3D EPI data of healthy subjects were acquired at 1-mm, 0.78-mm, and 0.65-mm isotropic resolution with varying echo train lengths (ETLs) and compared with a reference 3D GRE acquisition. RESULTS: The appearance of the susceptibility maps and the susceptibility values for segmented regions of interest were comparable between 3D EPI and 3D GRE acquisitions for both healthy and MS participants. Additionally, all lesions visible in the MS patients on the 3D GRE susceptibility maps were also visible on the 3D EPI susceptibility maps. The interplay among acquisition time, resolution, echo train length, and the effect of distortion on the calculated susceptibility maps was investigated. CONCLUSION: We demonstrate that the 3D EPI sequence is capable of rapidly acquiring submillimeter resolutions and providing high-quality, clinically relevant susceptibility maps.

Modelling the depth-dependent VASO and BOLD responses in human primary visual cortex.

Functional magnetic resonance imaging (fMRI) using a blood-oxygenation-level-dependent (BOLD) contrast is a common method for studying human brain function noninvasively. Gradient-echo (GRE) BOLD is highly sensitive to the blood oxygenation change in blood vessels; however, the spatial signal specificity can be degraded due to signal leakage from activated lower layers to superficial layers in depth-dependent (also called laminar or layer-specific) fMRI. Alternatively, physiological variables such as cerebral blood volume using the VAscular-Space-Occupancy (VASO) contrast have shown higher spatial specificity compared to BOLD. To better understand the physiological mechanisms such as blood volume and oxygenation changes and to interpret the measured depth-dependent responses, models are needed which reflect vascular properties at this scale. For this purpose, we extended and modified the "cortical vascular model" previously developed to predict layer-specific BOLD signal changes in human primary visual cortex to also predict a layer-specific VASO response. To evaluate the model, we compared the predictions with experimental results of simultaneous VASO and BOLD measurements in a group of healthy participants. Fitting the model to our experimental data provided an estimate of CBV change in different vascular compartments upon neural activity. We found that stimulus-evoked CBV change mainly occurs in small arterioles, capillaries, and intracortical arteries and that the contribution from venules and ICVs is smaller. Our results confirm that VASO is less susceptible to large vessel effects compared to BOLD, as blood volume changes in intracortical arteries did not substantially affect the resulting depth-dependent VASO profiles, whereas depth-dependent BOLD profiles showed a bias towards signal contributions from intracortical veins.

Unsupervised physiological noise correction of functional magnetic resonance imaging data using phase and magnitude information (PREPAIR).

Of the sources of noise affecting blood oxygen level-dependent functional magnetic resonance imaging (fMRI), respiration and cardiac fluctuations are responsible for the largest part of the variance, particularly at high and ultrahigh field. Existing approaches to removing physiological noise either use external recordings, which can be unwieldy and unreliable, or attempt to identify physiological noise from the magnitude fMRI data. Data-driven approaches are limited by sensitivity, temporal aliasing, and the need for user interaction. In the light of the sensitivity of the phase of the MR signal to local changes in the field stemming from physiological processes, we have developed an unsupervised physiological noise correction method using the information carried in the phase and the magnitude of echo-planar imaging data. Our technique, Physiological Regressor Estimation from Phase and mAgnItude, sub-tR (PREPAIR) derives time series signals sampled at the slice TR from both phase and magnitude images. It allows physiological noise to be captured without aliasing, and efficiently removes other sources of signal fluctuations not related to physiology, prior to regressor estimation. We demonstrate that the physiological signal time courses identified with PREPAIR agree well with those from external devices and retrieve challenging cardiac dynamics. The removal of physiological noise was as effective as that achieved with the most used approach based on external recordings, RETROICOR. In comparison with widely used recording-free physiological noise correction tools-PESTICA and FIX, both performed in unsupervised mode-PREPAIR removed significantly more respiratory and cardiac noise than PESTICA, and achieved a larger increase in temporal signal-to-noise-ratio at both 3 and 7 T.

Improved dynamic distortion correction for fMRI using single-echo EPI and a readout-reversed first image (REFILL).

The boundaries between tissues with different magnetic susceptibilities generate inhomogeneities in the main magnetic field which change over time due to motion, respiration and system instabilities. The dynamically changing field can be measured from the phase of the fMRI data and corrected. However, methods for doing so need multi-echo data, time-consuming reference scans and/or involve error-prone processing steps, such as phase unwrapping, which are difficult to implement robustly on the MRI host. The improved dynamic distortion correction method we propose is based on the phase of the single-echo EPI data acquired for fMRI, phase offsets calculated from a triple-echo, bipolar reference scan of circa 3-10 s duration using a method which avoids the need for phase unwrapping and an additional correction derived from one EPI volume in which the readout direction is reversed. This Reverse-Encoded First Image and Low resoLution reference scan (REFILL) approach is shown to accurately measure B0 as it changes due to shim, motion and respiration, even with large dynamic changes to the field at 7 T, where it led to a > 20% increase in time-series signal to noise ratio compared to data corrected with the classic static approach. fMRI results from REFILL-corrected data were free of stimulus-correlated distortion artefacts seen when data were corrected with static field mapping. The method is insensitive to shim changes and eddy current differences between the reference scan and the fMRI time series, and employs calculation steps that are simple and robust, allowing most data processing to be performed in real time on the scanner image reconstruction computer. These improvements make it feasible to routinely perform dynamic distortion correction in fMRI.

Mesoscopic in vivo human T2* dataset acquired using quantitative MRI at 7 Tesla.

Mesoscopic (0.1-0.5 mm) interrogation of the living human brain is critical for advancing neuroscience and bridging the resolution gap with animal models. Despite the variety of MRI contrasts measured in recent years at the mesoscopic scale, in vivo quantitative imaging of T2* has not been performed. Here we provide a dataset containing empirical T2* measurements acquired at 0.35 × 0.35 × 0.35 mm3 voxel resolution using 7 Tesla MRI. To demonstrate unique features and high quality of this dataset, we generate flat map visualizations that reveal fine-scale cortical substructures such as layers and vessels, and we report quantitative depth-dependent T2* (as well as R2*) values in primary visual cortex and auditory cortex that are highly consistent across subjects. This dataset is freely available at https://doi.org/10.17605/OSF.IO/N5BJ7, and may prove useful for anatomical investigations of the human brain, as well as for improving our understanding of the basis of the T2*-weighted (f)MRI signal.

Bayesian population receptive field modeling in human somatosensory cortex.

Somatosensation is fundamental to our ability to sense our body and interact with the world. Our body is continuously sampling the environment using a variety of receptors tuned to different features, and this information is routed up to primary somatosensory cortex. Strikingly, the spatial organization of the peripheral receptors in the body are well maintained, with the resulting representation of the body in the brain being referred to as the somatosensory homunculus. Recent years have seen considerable advancements in the field of high-resolution fMRI, which have enabled an increasingly detailed examination of the organization and properties of this homunculus. Here we combined advanced imaging techniques at ultra-high field (7T) with a recently developed Bayesian population receptive field (pRF) modeling framework to examine pRF properties in primary somatosensory cortex. In each subject, vibrotactile stimulation of the fingertips (i.e., the peripheral mechanoreceptors) modulated the fMRI response along the post-central gyrus and these signals were used to estimate pRFs. We found the pRF center location estimates to be in accord with previous work as well as evidence of other properties in line with the underlying neurobiology. Specifically, as expected from the known properties of cortical magnification, we find a larger representation of the index finger compared to the other stimulated digits (middle, index, little). We also show evidence that the little finger is marked by the largest pRF sizes, and that pRF size increases from anterior to posterior regions of S1. The ability to estimate somatosensory pRFs in humans provides an unprecedented opportunity to examine the neural mechanisms underlying somatosensation and is critical for studying how the brain, body, and environment interact to inform perception and action.

On the analysis of rapidly sampled fMRI data.

Recent advances in parallel imaging and simultaneous multi-slice techniques have permitted whole-brain fMRI acquisitions at sub-second sampling intervals, without significantly sacrificing the spatial coverage and resolution. Apart from probing brain function at finer temporal scales, faster sampling rates may potentially lead to enhanced functional sensitivity, owing possibly to both cleaner neural representations (due to less aliased physiological noise) and additional statistical benefits (due to more degrees of freedom for a fixed scan duration). Accompanying these intriguing aspects of fast acquisitions, however, confusion has also arisen regarding (1) how to preprocess/analyze these fast fMRI data, and (2) what exactly is the extent of benefits with fast acquisitions, i.e., how fast is fast enough for a specific research aim? The first question is motivated by the altered spectral distribution and noise characteristics at short sampling intervals, while the second question seeks to reconcile the complicated trade-offs between the functional contrast-to-noise ratio and the effective degrees of freedom. Although there have been recent efforts to empirically approach different aspects of these two questions, in this work we discuss, from a theoretical perspective accompanied by some illustrative, proof-of-concept experimental in vivo human fMRI data, a few considerations that are rarely mentioned, yet are important for both preprocessing and optimizing statistical inferences for studies that employ acquisitions with sub-second sampling intervals. Several summary recommendations include concerns regarding advisability of relying on low-pass filtering to de-noise physiological contributions, employment of statistical models with sufficient complexity to account for the substantially increased serial correlation, and cautions regarding using rapid sampling to enhance functional sensitivity given that different analysis models may associate with distinct trade-offs between contrast-to-noise ratios and the effective degrees of freedom. As an example, we demonstrate that as TR shortens, the intrinsic differences in how noise is accommodated in general linear models and Pearson correlation analyses (assuming Gaussian distributed stochastic signals and noise) can result in quite different outcomes, either gaining or losing statistical power.

Serial correlations in single-subject fMRI with sub-second TR.

When performing statistical analysis of single-subject fMRI data, serial correlations need to be taken into account to allow for valid inference. Otherwise, the variability in the parameter estimates might be under-estimated resulting in increased false-positive rates. Serial correlations in fMRI data are commonly characterized in terms of a first-order autoregressive (AR) process and then removed via pre-whitening. The required noise model for the pre-whitening depends on a number of parameters, particularly the repetition time (TR). Here we investigate how the sub-second temporal resolution provided by simultaneous multislice (SMS) imaging changes the noise structure in fMRI time series. We fit a higher-order AR model and then estimate the optimal AR model order for a sequence with a TR of less than 600 ms providing whole brain coverage. We show that physiological noise modelling successfully reduces the required AR model order, but remaining serial correlations necessitate an advanced noise model. We conclude that commonly used noise models, such as the AR(1) model, are inadequate for modelling serial correlations in fMRI using sub-second TRs. Rather, physiological noise modelling in combination with advanced pre-whitening schemes enable valid inference in single-subject analysis using fast fMRI sequences.

Using multi-echo simultaneous multi-slice (SMS) EPI to improve functional MRI of the subcortical nuclei of the basal ganglia at ultra-high field (7T).

The nuclei of the basal ganglia pose a special problem for functional MRI, especially at ultra-high field, because T2* variations between different regions result in suboptimal BOLD sensitivity when using gradient-echo echo-planar imaging (EPI). Specifically, the iron-rich lentiform nucleus of the basal ganglia, including the putamen and globus pallidus, suffers from substantial signal loss when imaging is performed using conventional single-echo EPI with echo times optimized for the cortex. Multi-echo EPI acquires several echoes at different echo times for every imaging slice, allowing images to be reconstructed with a weighting of echo times that is optimized individually for each voxel according to the underlying tissue or T2* properties. Here we show that multi-echo simultaneous multi-slice (SMS) EPI can improve functional activation of iron-rich subcortical regions while maintaining sensitivity within cortical areas. Functional imaging during a motor task known to elicit strong activations in the cortex and the subcortex (basal ganglia) was performed to compare the performance of multi-echo SMS EPI to single-echo SMS EPI. Notably within both the caudate nucleus and putamen of the basal ganglia, multi-echo SMS EPI yielded higher tSNR (an average 84% increase) and CNR (an average 58% increase), an approximate 3-fold increase in supra-threshold voxels, and higher t-values (an average 39% increase). The degree of improvement in the group level t-statistics was negatively correlated to the underlying T2* of the voxels, such that the shorter the T2*, as in the iron-rich nuclei of the basal ganglia, the higher the improvement of t-values in the activated region.

Measuring the effects of attention to individual fingertips in somatosensory cortex using ultra-high field (7T) fMRI.

Attention to sensory information has been shown to modulate the neuronal processing of that information. For example, visuospatial attention acts by modulating responses at retinotopically appropriate regions of visual cortex (Puckett and DeYoe, 2015; Tootell et al. 1998). Much less, however, is known about the neuronal processing associated with attending to other modalities of sensory information. One reason for this is that visual cortex is relatively large, and therefore easier to access non-invasively in humans using tools such as functional magnetic resonance imaging (fMRI). With high-resolution fMRI, however, it is now possible to access smaller cortical areas such as primary somatosensory cortex (Martuzzi et al., 2014; Sanchez-Panchuelo et al., 2010; Schweisfurth et al. 2014; Schweizer et al. 2008). Here, we combined a novel experimental design and high-resolution fMRI at ultra-high field (7T) to measure the effects of attention to tactile stimulation in primary somatosensory cortex, S1. We find that attention modulates somatotopically appropriate regions of S1, and importantly, that this modulation can be measured at the level of the cortical representation of individual fingertips.

Analysis and correction of field fluctuations in fMRI data using field monitoring.

This work investigates the role of magnetic field fluctuations as a confound in fMRI. In standard fMRI experiments with single-shot EPI acquisition at 3 Tesla the uniform and gradient components of the magnetic field were recorded with NMR field sensors. By principal component analysis it is found that differences of field evolution between the EPI readouts are explainable by few components relating to slow and within-shot field dynamics of hardware and physiological origin. The impact of fluctuating field components is studied by selective data correction and assessment of its influence on image fluctuation and SFNR. Physiological field fluctuations, attributed to breathing, were found to be small relative to those of hardware origin. The dominant confounds were hardware-related and attributable to magnet drift and thermal changes. In raw image time series, field fluctuation caused significant SFNR loss, reflected by a 67% gain upon correction. Large part of this correction can be accomplished by traditional image realignment, which addresses slow and spatially uniform field changes. With realignment, explicit field correction increased the SFNR on the order of 6%. In conclusion, field fluctuations are a relevant confound in fMRI and can be addressed effectively by retrospective data correction. Based on the physics involved it is anticipated that the advantage of full field correction increases with field strength, with non-Cartesian readouts, and upon phase-sensitive BOLD analysis.

Monitoring, analysis, and correction of magnetic field fluctuations in echo planar imaging time series.

PURPOSE: To assess the utility of concurrent magnetic field monitoring for observing and correcting for variations in k-space trajectories and global background fields that occur in single-shot echo planar imaging (EPI) time series as typically used in functional MRI (fMRI). METHODS: Field monitoring was performed using an array of NMR field probes operated concurrently with series of single-shot EPI acquisitions from a static phantom. The observed fluctuations in field evolution were analyzed in terms of their temporal and spatial behavior at the field level as well as at the level of reconstructed image series. The potential to correct for such fluctuations was assessed by accounting for them upon image reconstruction. An indication of the number and relative magnitude of underlying effects was obtained via principal component analysis. RESULTS: Trajectory and global field variations were found to induce substantial image fluctuations. Global field fluctuations induced standard deviations in image intensity up to 31%. Fluctuations in the trajectory induced ghosting artifacts with standard deviations up to 2%. Concurrent magnetic field monitoring reduced the fluctuations in the EPI time series to a maximum of 1.2%. CONCLUSION: Concurrent magnetic field monitoring holds the potential to improve the net sensitivity of fMRI by reducing signal fluctuations unrelated to brain activity.

VesselBoost: A Python Toolbox for Small Blood Vessel Segmentation in Human Magnetic Resonance Angiography Data.

Magnetic resonance angiography (MRA) performed at ultra-high magnetic field provides a unique opportunity to study the arteries of the living human brain at the mesoscopic level. From this, we can gain new insights into the brain's blood supply and vascular disease affecting small vessels. However, for quantitative characterization and precise representation of human angioarchitecture to, for example, inform blood-flow simulations, detailed segmentations of the smallest vessels are required. Given the success of deep learning-based methods in many segmentation tasks, we here explore their application to high-resolution MRA data, and address the difficulty of obtaining large data sets of correctly and comprehensively labelled data. We introduce VesselBoost, a vessel segmentation package, which utilizes deep learning and imperfect training labels for accurate vasculature segmentation. Combined with an innovative data augmentation technique, which leverages the resemblance of vascular structures, VesselBoost enables detailed vascular segmentations.

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.

Imaging of the pial arterial vasculature of the human brain in vivo using high-resolution 7T time-of-flight angiography.

The pial arterial vasculature of the human brain is the only blood supply to the neocortex, but quantitative data on the morphology and topology of these mesoscopic arteries (diameter 50-300 µm) remains scarce. Because it is commonly assumed that blood flow velocities in these vessels are prohibitively slow, non-invasive time-of-flight magnetic resonance angiography (TOF-MRA)-which is well suited to high 3D imaging resolutions-has not been applied to imaging the pial arteries. Here, we provide a theoretical framework that outlines how TOF-MRA can visualize small pial arteries in vivo, by employing extremely small voxels at the size of individual vessels. We then provide evidence for this theory by imaging the pial arteries at 140 µm isotropic resolution using a 7 Tesla (T) magnetic resonance imaging (MRI) scanner and prospective motion correction, and show that pial arteries one voxel width in diameter can be detected. We conclude that imaging pial arteries is not limited by slow blood flow, but instead by achievable image resolution. This study represents the first targeted, comprehensive account of imaging pial arteries in vivo in the human brain. This ultra-high-resolution angiography will enable the characterization of pial vascular anatomy across the brain to investigate patterns of blood supply and relationships between vascular and functional architecture.

New acquisition techniques and their prospects for the achievable resolution of fMRI.

This work reviews recent advances in technologies for functional magnetic resonance imaging (fMRI) of the human brain and highlights the push for higher functional specificity based on increased spatial resolution and specific MR contrasts to reveal previously undetectable functional properties of small-scale cortical structures. We discuss how the combination of MR hardware, advanced acquisition techniques and various MR contrast mechanisms have enabled recent progress in functional neuroimaging. However, these advanced fMRI practices have only been applied to a handful of neuroscience questions to date, with the majority of the neuroscience community still using conventional imaging techniques. We thus discuss upcoming challenges and possibilities for fMRI technology development in human neuroscience. We hope that readers interested in functional brain imaging acquire an understanding of current and novel developments and potential future applications, even if they don't have a background in MR physics or engineering. We summarize the capabilities of standard fMRI acquisition schemes with pointers to relevant literature and comprehensive reviews and introduce more recent developments.