Results of the 2020 fastMRI Challenge for Machine Learning MR Image Reconstruction.

Accelerating MRI scans is one of the principal outstanding problems in the MRI research community. Towards this goal, we hosted the second fastMRI competition targeted towards reconstructing MR images with subsampled k-space data. We provided participants with data from 7,299 clinical brain scans (de-identified via a HIPAA-compliant procedure by NYU Langone Health), holding back the fully-sampled data from 894 of these scans for challenge evaluation purposes. In contrast to the 2019 challenge, we focused our radiologist evaluations on pathological assessment in brain images. We also debuted a new Transfer track that required participants to submit models evaluated on MRI scanners from outside the training set. We received 19 submissions from eight different groups. Results showed one team scoring best in both SSIM scores and qualitative radiologist evaluations. We also performed analysis on alternative metrics to mitigate the effects of background noise and collected feedback from the participants to inform future challenges. Lastly, we identify common failure modes across the submissions, highlighting areas of need for future research in the MRI reconstruction community.

Trading off spatio-temporal properties in 3D high-speed fMRI using interleaved stack-of-spirals trajectories.

PURPOSE: Highly undersampled acquisitions have been proposed to push the limits of temporal resolution in functional MRI. This contribution is aimed at identifying parameter sets that let the user trade-off between ultra-high temporal resolution and spatial signal quality by varying the sampling densities. The proposed method maintains the synergies of a temporal resolution that enables direct filtering of physiological artifacts for highest statistical power, and 3D read-outs with optimal use of encoding capabilities of multi-coil arrays for efficient sampling and high signal-to-noise ratio (SNR). METHODS: One- to four-shot interleaved spherical stack-of-spiral trajectories with repetition times from 96 to 352 ms at a nominal resolution of 3 mm using different sampling densities were compared for image quality and temporal SNR (tSNR). The one- and three-shot trajectories were employed in a resting state study for functional characterization. RESULTS: Compared to a previously described single-shot trajectory, denser sampled trajectories of the same type are shown to be less prone to blurring and off-resonance vulnerability that appear in addition to the variable density artifacts of the point spread function. While the multi-shot trajectories lead to a decrease in tSNR efficiency, the high SNR due to the 3D read-out, combined with notable increases in image quality, leads to superior overall results of the three-shot interleaved stack of spirals. A resting state analysis of 15 subjects shows significantly improved functional sensitivity in areas of high off-resonance gradients. CONCLUSION: Mild variable-density sampling leads to excellent tSNR behavior and no increased off-resonance vulnerability, and is suggested unless maximum temporal resolution is sought.

Improving the sensitivity of spin-echo fMRI at 3T by highly accelerated acquisitions.

PURPOSE: Spin-echo (SE) functional MRI (fMRI) can be highly advantageous compared to gradient-echo (GE) fMRI with respect to magnetic field-inhomogeneity artifacts. However, at 3T, the majority of blood oxygenation level-dependent (BOLD) fMRI experiments are performed using T2∗ -weighted GE sequences because of their superior sensitivity compared to SE-fMRI. The presented SE implementation of a highly accelerated GE pulse sequence therefore aims to improve the sensitivity of SE-fMRI while profiting from a reduction of susceptibility-induced signal dropout. METHODS: Spin-echo MR encephalography (SE-MREG) is compared with the more conventionally used spin-echo echo-planar imaging (SE-EPI) and spin-echo simultaneous multislice (SE-SMS) at 3T in terms of capability to detect neuronal activations and resting-state functional connectivity. For activation analysis, healthy subjects underwent consecutive SE-MREG (pulse repetition time [TR] = 0.25 seconds), SE-SMS (TR = 1.3 seconds), and SE-EPI (TR = 4.4 seconds) scans in pseudorandomized order applied to a visual block design paradigm for generation of t-statistics maps. For the investigation of functional connectivity, additional resting-state data were acquired for 5 minutes and a seed-based correlation analysis using Stanford's FIND (Functional Imaging in Neuropsychiatric Disorders) atlas was performed. RESULTS: The increased sampling rate of SE-MREG relative to SE-SMS and SE-EPI improves the sensitivity to detect BOLD activation by 33% and 54%, respectively, and increases the capability to extract resting-state networks. Compared with a brain region that is not affected by magnetic field inhomogeneities, SE-MREG shows 2.5 times higher relative signal strength than GE-MREG in mesial temporal structures. CONCLUSION: SE-MREG offers a viable possibility for whole-brain fMRI with consideration of brain regions that are affected by strong susceptibility-induced magnetic field gradients.

15 Years MR-encephalography.

OBJECTIVE: This review article gives an account of the development of the MR-encephalography (MREG) method, which started as a mere 'Gedankenexperiment' in 2005 and gradually developed into a method for ultrafast measurement of physiological activities in the brain. After going through different approaches covering k-space with radial, rosette, and concentric shell trajectories we have settled on a stack-of-spiral trajectory, which allows full brain coverage with (nominal) 3 mm isotropic resolution in 100 ms. The very high acceleration factor is facilitated by the near-isotropic k-space coverage, which allows high acceleration in all three spatial dimensions. METHODS: The methodological section covers the basic sequence design as well as recent advances in image reconstruction including the targeted reconstruction, which allows real-time feedback applications, and-most recently-the time-domain principal component reconstruction (tPCR), which applies a principal component analysis of the acquired time domain data as a sparsifying transformation to improve reconstruction speed as well as quality. APPLICATIONS: Although the BOLD-response is rather slow, the high speed acquisition of MREG allows separation of BOLD-effects from cardiac and breathing related pulsatility. The increased sensitivity enables direct detection of the dynamic variability of resting state networks as well as localization of single interictal events in epilepsy patients. A separate and highly intriguing application is aimed at the investigation of the glymphatic system by assessment of the spatiotemporal patterns of cardiac and breathing related pulsatility. DISCUSSION: MREG has been developed to push the speed limits of fMRI. Compared to multiband-EPI this allows considerably faster acquisition at the cost of reduced image quality and spatial resolution.

The potential of MR-Encephalography for BCI/Neurofeedback applications with high temporal resolution.

Real-time functional magnetic resonance imaging (rt-fMRI) enables the update of various brain-activity measures during an ongoing experiment as soon as a new brain volume is acquired. However, the recorded Blood-oxygen-level dependent (BOLD) signal also contains physiological artifacts such as breathing and heartbeat, which potentially cause misleading false positive effects especially problematic in brain-computer interface (BCI) and neurofeedback (NF) setups. The low temporal resolution of echo planar imaging (EPI) sequences (which is in the range of seconds) prevents a proper separation of these artifacts from the BOLD signal. MR-Encephalography (MREG) has been shown to provide the high temporal resolution required to unalias and correct for physiological fluctuations and leads to increased specificity and sensitivity for mapping task-based activation and functional connectivity as well as for detecting dynamic changes in connectivity over time. By comparing a simultaneous multislice echo planar imaging (SMS-EPI) sequence and an MREG sequence using the same nominal spatial resolution in an offline analysis for three different experimental fMRI paradigms (perception of house and face stimuli, motor imagery, Stroop task), the potential of this novel technique for future BCI and NF applications was investigated. First, adapted general linear model pre-whitening which accounts for the high temporal resolution in MREG was implemented to calculate proper statistical results and be able to compare these with the SMS-EPI sequence. Furthermore, the respiration- and cardiac pulsation-related signals were successfully separated from the MREG signal using independent component analysis which were then included as regressors for a GLM analysis. Only the MREG sequence allowed to clearly separate cardiac pulsation and respiration components from the signal time course. It could be shown that these components highly correlate with the recorded respiration and cardiac pulsation signals using a respiratory belt and fingertip pulse plethysmograph. Temporal signal-to-noise ratios of SMS-EPI and MREG were comparable. Functional connectivity analysis using partial correlation showed a reduced standard error in MREG compared to SMS-EPI. Also, direct time course comparisons by down-sampling the MREG signal to the SMS-EPI temporal resolution showed lower variance in MREG. In general, we show that the higher temporal resolution is beneficial for fMRI time course modeling and this aspect can be exploited in offline application but also, is especially attractive, for real-time BCI and NF applications.

Targeted partial reconstruction for real-time fMRI with arbitrary trajectories.

PURPOSE: A partial image reconstruction formalism is introduced for the targeted extraction of real-time feedback from arbitrary trajectories when full image reconstruction in real time is computationally too demanding. METHODS: Explicit calculation and storage of linear combinations of lines of the reconstruction matrix by an incomplete basis change in spatial coordinates lead to translation of the expensive full reconstruction from a frame-wise application to a region of interest (ROI)-wise application. This step is independent from signal data and can be executed before the experiment. Subsequently, the results of the sum over fully reconstructed voxels can be evaluated directly. Data from a high-speed fMRI acquisition was used to investigate the targeted partial reconstruction of a functional ROI atlas, incorporating an intravolume dephasing correction. The same data and ROIs were used for a comparison of the time series obtained with those obtained from already existing methods for compartment-wise reconstruction. To examine real-time feasibility, the reconstruction was implemented and tested for online reconstruction performance. RESULTS: The reconstruction yields results that are virtually identical to the standard reconstruction (i.e., the magnitude sums over the ROIs), with negligible discrepancies even after termination of the conjugate gradient algorithm at a feasible number of iterations. Notably, more discrepancies arise with existing compartment-wise reconstructions. The online real-time implementation evaluated 1 ROI within 2.8 ms in the case of a highly parallel 3D whole brain acquisition. CONCLUSION: The high reconstruction fidelity and speed are satisfying for the exemplary application of real-time functional feedback using a highly parallel 3D whole brain acquisition.

Transition from spatial coherence to incoherence in coupled chaotic systems.

We investigate the spatio-temporal dynamics of coupled chaotic systems with nonlocal interactions, where each element is coupled to its nearest neighbors within a finite range. Depending upon the coupling strength and coupling radius, we find characteristic spatial patterns such as wavelike profiles and study the transition from coherence to incoherence leading to spatial chaos. We analyze the origin of this transition based on numerical simulations and support the results by theoretical derivations, identifying a critical coupling strength and a scaling relation of the coherent profiles. To demonstrate the universality of our findings, we consider time-discrete as well as time-continuous chaotic models realized as a logistic map and a Rössler or Lorenz system, respectively. Thereby, we establish the coherence-incoherence transition in networks of coupled identical oscillators.