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.

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.

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.

Phase-based masking for quantitative susceptibility mapping of the human brain at 9.4T.

PURPOSE: To develop improved tissue masks for QSM. METHODS: Masks including voxels at the brain surface were automatically generated from the magnitude alone (MM) or combined with test functions from the first (PG) or second (PB) derivative of the sign of the wrapped phase. Phase images at 3T and 9.4T were simulated at different TEs and used to generate a mask, PItoh , with between-voxel phase differences less than π. MM, PG, and PB were compared with PItoh . QSM were generated from 3D multi-echo gradient-echo data acquired at 9.4T (21 subjects aged: 20-56y), and from the QSM2016 challenge 3T data using different masks, unwrapping, background removal, and dipole inversion algorithms. QSM contrast was quantified using age-based iron concentrations. RESULTS: Close to air cavities, phase wraps became denser with increasing field and echo time, yielding increased values of the test functions. Compared with PItoh , PB had the highest Dice coefficient, while PG had the lowest and MM the highest percentage of voxels outside PItoh. Artifacts observed in QSM at 9.4T with MM were mitigated by stronger background filters but yielded a reduced QSM contrast. With PB, QSM contrast was greater and artifacts diminished. Similar results were obtained with challenge data, evidencing larger effects of mask close to air cavities. CONCLUSION: Automatic, phase-based masking founded on the second derivative of the sign of the wrapped phase, including cortical voxels at the brain surface, was able to mitigate artifacts and restore QSM contrast across cortical and subcortical brain regions.

QSMxT: Robust masking and artifact reduction for quantitative susceptibility mapping.

PURPOSE: Quantitative susceptibility mapping (QSM) estimates the spatial distribution of tissue magnetic susceptibilities from the phase of a gradient-echo signal. QSM algorithms require a signal mask to delineate regions with reliable phase for subsequent susceptibility estimation. Existing masking techniques used in QSM have limitations that introduce artifacts, exclude anatomical detail, and rely on parameter tuning and anatomical priors that narrow their application. Here, a robust masking and reconstruction procedure is presented to overcome these limitations and enable automated QSM processing. Moreover, this method is integrated within an open-source software framework: QSMxT. METHODS: A robust masking technique that automatically separates reliable from less reliable phase regions was developed and combined with a two-pass reconstruction procedure that operates on the separated sources before combination, extracting more information and suppressing streaking artifacts. RESULTS: Compared with standard masking and reconstruction procedures, the two-pass inversion reduces streaking artifacts caused by unreliable phase and high dynamic ranges of susceptibility sources. It is also robust across a range of acquisitions at 3 T in volunteers and phantoms, at 7 T in tumor patients, and in an in silico head phantom, with significant artifact and error reductions, greater anatomical detail, and minimal parameter tuning. CONCLUSION: The two-pass masking and reconstruction procedure separates reliable from less reliable phase regions, enabling a more accurate QSM reconstruction that mitigates artifacts, operates without anatomical priors, and requires minimal parameter tuning. The technique and its integration within QSMxT makes QSM processing more accessible and robust to streaking artifacts.

Feasibility of Hepatic Fat Quantification Using Proton Density Fat Fraction by Multi-Echo Chemical-Shift-Encoded MRI at 7T.

Fat fraction quantification and assessment of its distribution in the hepatic tissue become more important with the growing epidemic of obesity, and the increasing prevalence of diabetes mellitus type 2 and non-alcoholic fatty liver disease. At 3Tesla, the multi-echo, chemical-shift-encoded magnetic resonance imaging (CSE-MRI)-based acquisition allows the measurement of proton density fat-fraction (PDFF) even in clinical protocols. Further improvements in SNR can be achieved by the use of phased array coils and increased static magnetic field. The purpose of the study is to evaluate the feasibility of PDFF imaging using a multi-echo CSE-MRI technique at ultra-high magnetic field (7Tesla). Thirteen volunteers (M/F) with a broad range of age, body mass index, and hepatic PDFF were measured at 3 and 7T by multi-gradient-echo MRI and single-voxel spectroscopy MRS. All measurements were performed in breath-hold (exhalation); the MRI protocols were optimized for a short measurement time, thus minimizing motion-related problems. 7T data were processed off-line using Matlab® (MRI:multi-gradient-echo) and jMRUI (MRS), respectively. For quantitative validation of the PDFF results, a similar protocol was performed at 3T, including on-line data processing provided by the system manufacturer, and correlation analyses between 7 and 3T data were performed off-line. The multi-echo CSE-MRI measurements at 7T with a phased-array coil configuration and an optimal post-processing yielded liver volume coverage ranging from 30 to 90% for high- and low-BMI subjects, respectively. PDFFs ranged between 1 and 20%. We found significant correlations between 7T MRI and -MRS measurements (R2 ≅ 0.97; p < 0.005), and between MRI-PDFF at 7T and 3T fields (R2 ≅ 0.94; p < 0.005) in the evaluated volumes. Based on the measurements and analyses performed, the multi-echo CSE-MRI method using a 32-channel coil at 7T showed its aptitude for MRI-based quantitation of PDFF in the investigated volumes. The results are the first step toward qMRI of the whole liver at 7T with further improvements in hardware.

Improved susceptibility weighted imaging at ultra-high field using bipolar multi-echo acquisition and optimized image processing: CLEAR-SWI.

PURPOSE: Susceptibility Weighted Imaging (SWI) has become established in the clinical investigation of stroke, microbleeds, tumor vascularization, calcification and iron deposition, but suffers from a number of shortcomings and artefacts. The goal of this study was to reduce the sensitivity of SWI to strong B1 and B0 inhomogeneities at ultra-high field to generate homogeneous images with increased contrast and free of common artefacts. All steps in SWI processing have been addressed - coil combination, phase unwrapping, image combination over echoes, phase filtering and homogeneity correction - and applied to an efficient bipolar multi-echo acquisition to substantially improve the quality of SWI. PRINCIPAL RESULTS: Our findings regarding the optimal individual processing steps lead us to propose a Contrast-weighted, Laplace-unwrapped, bipolar multi-Echo, ASPIRE-combined, homogeneous, improved Resolution SWI, or CLEAR-SWI. CLEAR-SWI was compared to two other multi-echo SWI methods and standard, single-echo SWI with the same acquisition time at 7 T in 10 healthy volunteers and with single-echo SWI in 13 patients with brain tumors. CLEAR-SWI had improved contrast-to-noise and homogeneity, reduced signal dropout and was not compromised by the artefacts which affected standard SWI in 10 out of 13 cases close to tumors (as assessed by expert raters), as well as generating T2* maps and phase images which can be used for Quantitative Susceptibility Mapping. In a comparison with other multi-echo SWI methods, CLEAR-SWI had the fewest artefacts, highest SNR and generally higher contrast-to-noise. MAJOR CONCLUSIONS: CLEAR-SWI eliminates the artefacts common in standard, single-echo SWI, reduces signal dropouts and improves image homogeneity and contrast-to-noise. Applied clinically, in a study of brain tumor patients, CLEAR-SWI was free of the artefacts which affected standard, single-echo SWI.

The traveling heads 2.0: Multicenter reproducibility of quantitative imaging methods at 7 Tesla.

OBJECT: This study evaluates inter-site and intra-site reproducibility at ten different 7 T sites for quantitative brain imaging. MATERIAL AND METHODS: Two subjects - termed the "traveling heads" - were imaged at ten different 7 T sites with a harmonized quantitative brain MR imaging protocol. In conjunction with the system calibration, MP2RAGE, QSM, CEST and multi-parametric mapping/relaxometry were examined. RESULTS: Quantitative measurements with MP2RAGE showed very high reproducibility across sites and subjects, and errors were in concordance with previous results and other field strengths. QSM had high inter-site reproducibility for relevant subcortical volumes. CEST imaging revealed systematic differences between the sites, but reproducibility was comparable to results in the literature. Relaxometry had also very high agreement between sites, but due to the high sensitivity, differences caused by different applications of the B1 calibration of the two RF coil types used were observed. CONCLUSION: Our results show that quantitative brain imaging can be performed with high reproducibility at 7 T and with similar reliability as found at 3 T for multicenter studies of the supratentorial brain.

Phase unwrapping with a rapid opensource minimum spanning tree algorithm (ROMEO).

PURPOSE: To develop a rapid and accurate MRI phase-unwrapping technique for challenging phase topographies encountered at high magnetic fields, around metal implants, or postoperative cavities, which is sufficiently fast to be applied to large-group studies including Quantitative Susceptibility Mapping and functional MRI (with phase-based distortion correction). METHODS: The proposed path-following phase-unwrapping algorithm, ROMEO, estimates the coherence of the signal both in space-using MRI magnitude and phase information-and over time, assuming approximately linear temporal phase evolution. This information is combined to form a quality map that guides the unwrapping along a 3D path through the object using a computationally efficient minimum spanning tree algorithm. ROMEO was tested against the two most commonly used exact phase-unwrapping methods, PRELUDE and BEST PATH, in simulated topographies and at several field strengths: in 3T and 7T in vivo human head images and 9.4T ex vivo rat head images. RESULTS: ROMEO was more reliable than PRELUDE and BEST PATH, yielding unwrapping results with excellent temporal stability for multi-echo or multi-time-point data. It does not require image masking and delivers results within seconds, even in large, highly wrapped multi-echo data sets (eg, 9 seconds for a 7T head data set with 31 echoes and a 208 × 208 × 96 matrix size). CONCLUSION: Overall, ROMEO was both faster and more accurate than PRELUDE and BEST PATH, delivering exact results within seconds, which is well below typical image acquisition times, enabling potential on-console application.

Real-time motion and retrospective coil sensitivity correction for CEST using volumetric navigators (vNavs) at 7T.

PURPOSE: To explore the impact of temporal motion-induced coil sensitivity changes on CEST-MRI at 7T and its correction using interleaved volumetric EPI navigators, which are applied for real-time motion correction. METHODS: Five healthy volunteers were scanned via CEST. A 4-fold correction pipeline allowed the mitigation of (1) motion, (2) motion-induced coil sensitivity variations, ΔB1- , (3) motion-induced static magnetic field inhomogeneities, ΔB0 , and (4) spatially varying transmit RF field fluctuations, ΔB1+ . Four CEST measurements were performed per session. For the first 2, motion correction was turned OFF and then ON in absence of voluntary motion, whereas in the other 2 controlled head rotations were performed. During post-processing ΔB1- was removed additionally for the motion-corrected cases, resulting in a total of 6 scenarios to be compared. In all cases, retrospective ∆B0 and - ΔB1+ corrections were performed to compute artifact-free magnetization transfer ratio maps with asymmetric analysis (MTRasym ). RESULTS: Dynamic ΔB1- correction successfully mitigated signal deviations caused by head motion. In 2 frontal lobe regions of volunteer 4, induced relative signal errors of 10.9% and 3.9% were reduced to 1.1% and 1.0% after correction. In the right frontal lobe, the motion-corrected MTRasym contrast deviated 0.92%, 1.21%, and 2.97% relative to the static case for Δω = 1, 2, 3 ± 0.25 ppm. The additional application of ΔB1- correction reduced these deviations to 0.10%, 0.14%, and 0.42%. The fully corrected MTRasym values were highly consistent between measurements with and without intended head rotations. CONCLUSION: Temporal ΔB1- cause significant CEST quantification bias. The presented correction pipeline including the proposed retrospective ΔB1- correction significantly reduced motion-related artifacts on CEST-MRI.

Intra-session and inter-subject variability of 3D-FID-MRSI using single-echo volumetric EPI navigators at 3T.

PURPOSE: In this study, we demonstrate the first combination of 3D FID proton MRSI and spatial encoding via concentric-ring trajectories (CRTs) at 3T. FID-MRSI has many benefits including high detection sensitivity, in particular for J-coupled metabolites (e.g., glutamate/glutamine). This makes it highly attractive, not only for clinical, but also for, potentially, functional MRSI. However, this requires excellent reliability and temporal stability. We have, therefore, augmented this 3D-FID-MRSI sequence with single-echo, imaging-based volumetric navigators (se-vNavs) for real-time motion/shim-correction (SHMOCO), which is 2× quicker than the original double-echo navigators (de-vNavs), hence allowing more efficient integration also in short-TR sequences. METHODS: The tracking accuracy (position and B0 -field) of our proposed se-vNavs was compared to the original de-vNavs in phantoms (rest and translation) and in vivo (voluntary head rotation). Finally, the intra-session stability of a 5:40 min 3D-FID-MRSI scan was evaluated with SHMOCO and no correction (NOCO) in 5 resting subjects. Intra/inter-subject coefficients of variation (CV) and intra-class correlations (ICC) over the whole 3D volume and in selected regions of interest ROI were assessed. RESULTS: Phantom and in vivo scans showed highly consistent tracking performance for se-vNavs compared to the original de-vNavs, but lower frequency drift. Up to ~30% better intra-subject CVs were obtained for SHMOCO (P < 0.05), with values of 9.3/6.9/6.5/7.8% over the full VOI for Glx/tNAA/tCho/m-Ins ratios to tCr. ICCs were good-to-high (91% for Glx/tCr in motor cortex), whereas the inter-subject variability was ~11-19%. CONCLUSION: Real-time motion/shim corrected 3D-FID-MRSI with time-efficient CRT-sampling at 3T allows reliable, high-resolution metabolic imaging that is fast enough for clinical use and even, potentially, for functional MRSI.

The influence of spatial resolution on the spectral quality and quantification accuracy of whole-brain MRSI at 1.5T, 3T, 7T, and 9.4T.

PURPOSE: Inhomogeneities in the static magnetic field (B0 ) deteriorate MRSI data quality by lowering the spectral resolution and SNR. MRSI with low spatial resolution is also prone to lipid bleeding. These problems are increasingly problematic at ultra-high fields. An approach to tackling these challenges independent of B0 -shim hardware is to increase the spatial resolution. Therefore, we investigated the effect of improved spatial resolution on spectral quality and quantification at 4 field strengths. METHODS: Whole-brain MRSI data was simulated for 3 spatial resolutions and 4 B0 s based on experimentally acquired MRI data and simulated free induction decay signals of metabolites and lipids. To compare the spectral quality and quantification, we derived SNR normalized to the voxel size (nSNR), linewidth and metabolite concentration ratios, their Cramer-Rao-lower-bounds (CRLBs), and the absolute percentage error (APE) of estimated concentrations compared to the gold standard for the whole-brain and 8 brain regions. RESULTS: At 7T, we found up to a 3.4-fold improved nSNR (in the frontal lobe) and a 2.8-fold reduced linewidth (in the temporal lobe) for 1 cm3 versus 0.25 cm3 resolution. This effect was much more pronounced at higher and less homogenous B0 (1.6-fold improved nSNR and 1.8-fold improved linewidth in the parietal lobe at 3T). This had direct implications for quantification: the volume of reliably quantified spectra increased with resolution by 1.2-fold and 1.5-fold (when thresholded by CRLBs or APE, respectively). CONCLUSION: MRSI data quality benefits from increased spatial resolution particularly at higher B0 , and leads to more reliable metabolite quantification. In conjunction with the development of better B0 shimming hardware, this will enable robust whole-brain MRSI at ultra-high field.

Computationally Efficient Combination of Multi-channel Phase Data From Multi-echo Acquisitions (ASPIRE).

PURPOSE: To develop a simple method for combining multi-echo phase information from a number of coils in an array that requires no volume coil or additional scans and yields signal-to-noise ratio-optimal images that reflect only ΔB0-related phase. THEORY AND METHODS: Two SNR optimal coil combination methods were developed which retrieve the ΔB0-related phase by determining the coil-dependent phase offsets. The first variant, MCPC-3D-S, requires the unwrapping of one phase image; the second variant, ASPIRE, allows unwrapping to be avoided if two echoes j and k satisfy the echo time relation m⋅TEk=(m+1)⋅TEj, where m is an integer, making this a particularly fast and robust approach. Both developed methods constitute improvements over a prior method, MCPC-3D, in terms of robustness and computational expense. RESULTS: In the brain at 7 T, phase matching and contrast-to-noise ratio were higher with MCPC-3D-S and ASPIRE than with phase difference reconstruction, and similar to the reference coil-dependent Roemer combination. Unlike the Roemer and virtual reference coil methods, the proposed approaches also eliminated all non- ΔB0-related phase. CONCLUSION: MCPC-3D-S is an improvement over prior multi-echo methods, which is useful if the ASPIRE echo time condition cannot be fulfilled. ASPIRE is a particularly fast and robust approach that runs on the scanner's reconstructor in a small fraction of the acquisition time. Magn Reson Med 79:2996-3006, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Combining two open source tools for neural computation (BioPatRec and Netlab) improves movement classification for prosthetic control.

BACKGROUND: Controlling a myoelectric prosthesis for upper limbs is increasingly challenging for the user as more electrodes and joints become available. Motion classification based on pattern recognition with a multi-electrode array allows multiple joints to be controlled simultaneously. Previous pattern recognition studies are difficult to compare, because individual research groups use their own data sets. To resolve this shortcoming and to facilitate comparisons, open access data sets were analysed using components of BioPatRec and Netlab pattern recognition models. METHODS: Performances of the artificial neural networks, linear models, and training program components were compared. Evaluation took place within the BioPatRec environment, a Matlab-based open source platform that provides feature extraction, processing and motion classification algorithms for prosthetic control. The algorithms were applied to myoelectric signals for individual and simultaneous classification of movements, with the aim of finding the best performing algorithm and network model. Evaluation criteria included classification accuracy and training time. RESULTS: Results in both the linear and the artificial neural network models demonstrated that Netlab's implementation using scaled conjugate training algorithm reached significantly higher accuracies than BioPatRec. CONCLUSIONS: It is concluded that the best movement classification performance would be achieved through integrating Netlab training algorithms in the BioPatRec environment so that future prosthesis training can be shortened and control made more reliable. Netlab was therefore included into the newest release of BioPatRec (v4.0).

Correcting dynamic distortions in 7T echo planar imaging using a jittered echo time sequence.

PURPOSE: To develop a distortion correction method for echo planar imaging (EPI) that is able to measure dynamic changes in B0 . THEORY AND METHODS: The approach we propose is based on single-echo EPI with a jittering of the echo time between two values for alternate time points. Field maps are calculated between phase images from adjacent volumes and are used to remove distortion from corresponding magnitude images. The performance of our approach was optimized using an analytical model and by comparison with field maps from dual-echo EPI. The method was tested in functional MRI experiments at 7T with motor tasks and compared with the conventional static approach. RESULTS: Unwarping using our method was accurate even for head rotations up to 8.2°, where the static approach introduced errors up to 8.2 mm. Jittering the echo time between 19 and 25 ms had no measurable effect on blood oxygenation level-dependent (BOLD) sensitivity. Our approach reduced the distortions in activated regions to <1 mm and repositioned active voxels correctly. CONCLUSION: This method yields accurate distortion correction in the presence of motion. No reduction in BOLD sensitivity was observed. As such, it is suitable for application in a wide range of functional MRI experiments. Magn Reson Med 76:1388-1399, 2016. © 2015 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.