Dysfunction of the magnocellular subdivision of the visual thalamus in developmental dyslexia.

Developmental dyslexia (DD) is one of the most common learning disorders, affecting millions of children and adults worldwide. To date, scientific research has attempted to explain DD primarily based on pathophysiological alterations in the cerebral cortex. In contrast, several decades ago, pioneering research on five post-mortem human brains suggested that a core characteristic of DD might be morphological alterations in a specific subdivision of the visual thalamus - the magnocellular LGN (M-LGN). However, due to considerable technical challenges in investigating LGN subdivisions non-invasively in humans, this finding was never confirmed in-vivo, and its relevance for DD pathology remained highly controversial. Here, we leveraged recent advances in high-resolution magnetic resonance imaging (MRI) at high field strength (7 Tesla) to investigate the M-LGN in DD in-vivo. Using a case-control design, we acquired data from a large sample of young adults with DD (n = 26; age 28 ± 7 years; 13 females) and matched control participants (n = 28; age 27 ± 6 years; 15 females). Each participant completed a comprehensive diagnostic behavioral test battery and participated in two MRI sessions, including three functional MRI experiments and one structural MRI acquisition. We measured blood-oxygen-level-dependent responses and longitudinal relaxation rates to compare both groups on LGN subdivision function and myelination. Based on previous research, we hypothesized that the M-LGN is altered in DD and that these alterations are associated with a key DD diagnostic score, i.e., rapid letter and number naming (RANln). The results showed aberrant responses of the M-LGN in DD compared to controls, which was reflected in a different functional lateralization of this subdivision between groups. These alterations were associated with RANln performance, specifically in male DD. We also found lateralization differences in the longitudinal relaxation rates of the M-LGN in DD relative to controls. Conversely, the other main subdivision of the LGN, the parvocellular LGN (P-LGN), showed comparable blood-oxygen-level-dependent responses and longitudinal relaxation rates between groups. The present study is the first to unequivocally show that M-LGN alterations are a hallmark of DD, affecting both the function and microstructure of this subdivision. It further provides a first functional interpretation of M-LGN alterations and a basis for a better understanding of sex-specific differences in DD with implications for prospective diagnostic and treatment strategies.

Distance from main arteries influences microstructural and functional brain tissue characteristics.

Given the substantial dependence of neurons on continuous supply of energy, the distribution of major cerebral arteries opens a question whether the distance from the main supply arteries constitutes a modulating factor for the microstructural and functional properties of brain tissue. To tackle this question, multimodal MRI acquisitions of 102 healthy volunteers over the full adult age span were utilised. Relaxation along a fictitious field in the rotating frame of rank n = 4 (RAFF4), adiabatic T1ρ, T2ρ,  and intracellular volume fraction (fICVF) derived from diffusion-weighted imaging were implemented to quantify microstructural (cellularity, myelin density, iron concentration) tissue characteristics and degree centrality and fractional amplitude of low-frequency fluctuations to probe for functional metrics. Inverse correlation of arterial distance with robust homogeneity was detected for T1ρ, T2ρ and RAFF4 for cortical grey matter and white matter, showing substantial complex microstructural differences between brain tissue close and farther from main arterial trunks. Albeit with wider variability, functional metrics pointed to increased connectivity and neuronal activity in areas farther from main arteries. Surprisingly, multiple of these microstructural and functional distance-based gradients diminished with higher age, pointing to uniformization of brain tissue with ageing. All in all, this pilot study provides a novel insight on brain regionalisation based on artery distance, which merits further investigation to validate its biological underpinnings.

Model-based super-resolution reconstruction for pseudo-continuous Arterial Spin Labeling.

Arterial spin labeling (ASL) is a promising, non-invasive perfusion magnetic resonance imaging technique for quantifying cerebral blood flow (CBF). Unfortunately, ASL suffers from an inherently low signal-to-noise ratio (SNR) and spatial resolution, undermining its potential. Increasing spatial resolution without significantly sacrificing SNR or scan time represents a critical challenge towards routine clinical use. In this work, we propose a model-based super-resolution reconstruction (SRR) method with joint motion estimation that breaks the traditional SNR/resolution/scan-time trade-off. From a set of differently oriented 2D multi-slice pseudo-continuous ASL images with a low through-plane resolution, 3D-isotropic, high resolution, quantitative CBF maps are estimated using a Bayesian approach. Experiments on both synthetic whole brain phantom data, and on in vivo brain data, show that the proposed SRR Bayesian estimation framework outperforms state-of-the-art ASL quantification.

High-fidelity intravoxel incoherent motion parameter mapping using locally low-rank and subspace modeling.

PURPOSE: Intravoxel incoherent motion (IVIM) is a quantitative magnetic resonance imaging (MRI) method used to quantify perfusion properties of tissue non-invasively without contrast. However, clinical applications are limited by unreliable parameter estimates, particularly for the perfusion fraction (f) and pseudodiffusion coefficient (D*). This study aims to develop a high-fidelity reconstruction for reliable estimation of IVIM parameters. The proposed method is versatile and amenable to various acquisition schemes and fitting methods. METHODS: To address current challenges with IVIM, we adapted several advanced reconstruction techniques. We used a low-rank approximation of IVIM images and temporal subspace modeling to constrain the magnetization dynamics of the bi-exponential diffusion signal decay. In addition, motion-induced phase variations were corrected between diffusion directions and b-values, facilitating the use of high SNR real-valued diffusion data. The proposed method was evaluated in simulations and in vivo brain acquisitions in six healthy subjects and six individuals with a history of SARS-CoV-2 infection and compared with the conventionally reconstructed magnitude data. Following reconstruction, IVIM parameters were estimated voxel-wise. RESULTS: Our proposed method reduced noise contamination in simulations, resulting in a 60%, 58.9%, and 83.9% reduction in the NRMSE for D, f, and D*, respectively, compared to the conventional reconstruction. In vivo, anisotropic properties of D, f, and D* were preserved with the proposed method, highlighting microvascular differences in gray matter between individuals with a history of COVID-19 and those without (p = 0.0210), which wasn't observed with the conventional reconstruction. CONCLUSION: The proposed method yielded a more reliable estimation of IVIM parameters with less noise than the conventional reconstruction. Further, the proposed method preserved anisotropic properties of IVIM parameter estimates and demonstrated differences in microvascular perfusion in COVID-affected subjects, which weren't observed with conventional reconstruction methods.

Clarifying Human White Matter.

Progress in magnetic resonance imaging (MRI) now makes it possible to identify the major white matter tracts in the living human brain. These tracts are important because they carry many of the signals communicated between different brain regions. MRI methods coupled with biophysical modeling can measure the tissue properties and structural features of the tracts that impact our ability to think, feel, and perceive. This review describes the fundamental ideas of the MRI methods used to identify the major white matter tracts in the living human brain.

Rapid quantitative magnetization transfer imaging: Utilizing the hybrid state and the generalized Bloch model.

PURPOSE: To explore efficient encoding schemes for quantitative magnetization transfer (qMT) imaging with few constraints on model parameters. THEORY AND METHODS: We combine two recently proposed models in a Bloch-McConnell equation: the dynamics of the free spin pool are confined to the hybrid state, and the dynamics of the semi-solid spin pool are described by the generalized Bloch model. We numerically optimize the flip angles and durations of a train of radio frequency pulses to enhance the encoding of three qMT parameters while accounting for all eight parameters of the two-pool model. We sparsely sample each time frame along this spin dynamics with a three-dimensional radial koosh-ball trajectory, reconstruct the data with subspace modeling, and fit the qMT model with a neural network for computational efficiency. RESULTS: We extracted qMT parameter maps of the whole brain with an effective resolution of 1.24 mm from a 12.6-min scan. In lesions of multiple sclerosis subjects, we observe a decreased size of the semi-solid spin pool and longer relaxation times, consistent with previous reports. CONCLUSION: The encoding power of the hybrid state, combined with regularized image reconstruction, and the accuracy of the generalized Bloch model provide an excellent basis for efficient quantitative magnetization transfer imaging with few constraints on model parameters.

Minimization of eddy current artifacts in sequences with periodic dynamics.

PURPOSE: To minimize eddy current artifacts in periodic pulse sequences with balanced gradient moments as, for example, used for quantitative MRI. THEORY AND METHODS: Eddy current artifacts in balanced sequences result from large jumps in k-space. In quantitative MRI, one often samples some spin dynamics repeatedly while acquiring different parts of k-space. We swap individual k-space lines between different repetitions in order to minimize jumps in temporal succession without changing the overall trajectory. This reordering can be formulated as a traveling salesman problem and we tackle the discrete optimization with a simulated annealing algorithm. RESULTS: Compared to the default ordering, we observe a substantial reduction of artifacts in the reconstructed images and the derived quantitative parameter maps. Comparing two variants of our algorithm, one that resembles the pairing approach originally proposed by Bieri et al., and one that minimizes all k-space jumps equally, we observe slightly lower artifact levels in the latter. CONCLUSION: The proposed reordering scheme effectively reduces eddy current artifacts in sequences with balanced gradient moments. In contrast to previous approaches, we capitalize on the periodicity of the sampled signal dynamics, enabling both efficient k-space sampling and minimizing artifacts caused by eddy currents.

Improving quantitative MRI using self-supervised deep learning with model reinforcement: Demonstration for rapid T1 mapping.

PURPOSE: This paper proposes a novel self-supervised learning framework that uses model reinforcement, REference-free LAtent map eXtraction with MOdel REinforcement (RELAX-MORE), for accelerated quantitative MRI (qMRI) reconstruction. The proposed method uses an optimization algorithm to unroll an iterative model-based qMRI reconstruction into a deep learning framework, enabling accelerated MR parameter maps that are highly accurate and robust. METHODS: Unlike conventional deep learning methods which require large amounts of training data, RELAX-MORE is a subject-specific method that can be trained on single-subject data through self-supervised learning, making it accessible and practically applicable to many qMRI studies. Using quantitative T 1 $$ {\mathrm{T}}_1 $$ mapping as an example, the proposed method was applied to the brain, knee and phantom data. RESULTS: The proposed method generates high-quality MR parameter maps that correct for image artifacts, removes noise, and recovers image features in regions of imperfect image conditions. Compared with other state-of-the-art conventional and deep learning methods, RELAX-MORE significantly improves efficiency, accuracy, robustness, and generalizability for rapid MR parameter mapping. CONCLUSION: This work demonstrates the feasibility of a new self-supervised learning method for rapid MR parameter mapping, that is readily adaptable to the clinical translation of qMRI.

High-resolution myelin-water fraction and quantitative relaxation mapping using 3D ViSTa-MR fingerprinting.

PURPOSE: This study aims to develop a high-resolution whole-brain multi-parametric quantitative MRI approach for simultaneous mapping of myelin-water fraction (MWF), T1, T2, and proton-density (PD), all within a clinically feasible scan time. METHODS: We developed 3D visualization of short transverse relaxation time component (ViSTa)-MRF, which combined ViSTa technique with MR fingerprinting (MRF), to achieve high-fidelity whole-brain MWF and T1/T2/PD mapping on a clinical 3T scanner. To achieve fast acquisition and memory-efficient reconstruction, the ViSTa-MRF sequence leverages an optimized 3D tiny-golden-angle-shuffling spiral-projection acquisition and joint spatial-temporal subspace reconstruction with optimized preconditioning algorithm. With the proposed ViSTa-MRF approach, high-fidelity direct MWF mapping was achieved without a need for multicompartment fitting that could introduce bias and/or noise from additional assumptions or priors. RESULTS: The in vivo results demonstrate the effectiveness of the proposed acquisition and reconstruction framework to provide fast multi-parametric mapping with high SNR and good quality. The in vivo results of 1 mm- and 0.66 mm-isotropic resolution datasets indicate that the MWF values measured by the proposed method are consistent with standard ViSTa results that are 30× slower with lower SNR. Furthermore, we applied the proposed method to enable 5-min whole-brain 1 mm-iso assessment of MWF and T1/T2/PD mappings for infant brain development and for post-mortem brain samples. CONCLUSIONS: In this work, we have developed a 3D ViSTa-MRF technique that enables the acquisition of whole-brain MWF, quantitative T1, T2, and PD maps at 1 and 0.66 mm isotropic resolution in 5 and 15 min, respectively. This advancement allows for quantitative investigations of myelination changes in the brain.

GPU-accelerated Bloch simulations and MR-STAT reconstructions using the Julia programming language.

PURPOSE: MR-STAT is a relatively new multiparametric quantitative MRI technique in which quantitative paramater maps are obtained by solving a large-scale nonlinear optimization problem. Managing reconstruction times is one of the main challenges of MR-STAT. In this work we leverage GPU hardware to reduce MR-STAT reconstruction times. A highly optimized, GPU-compatible Bloch simulation toolbox is developed as part of this work that can be utilized for other quantitative MRI techniques as well. METHODS: The Julia programming language was used to develop a flexible yet highly performant and GPU-compatible Bloch simulation toolbox called BlochSimulators.jl. The runtime performance of the toolbox is benchmarked against other Bloch simulation toolboxes. Furthermore, a (partially matrix-free) modification of a previously presented (matrix-free) MR-STAT reconstruction algorithm is proposed and implemented using the Julia language on GPU hardware. The proposed algorithm is combined with BlochSimulators.jl and the resulting MR-STAT reconstruction times on GPU hardware are compared to previously presented MR-STAT reconstruction times. RESULTS: The BlochSimulators.jl package demonstrates superior runtime performance on both CPU and GPU hardware when compared to other existing Bloch simulation toolboxes. The GPU-accelerated partially matrix-free MR-STAT reconstruction algorithm, which relies on BlochSimulators.jl, allows for reconstructions of 68 seconds per two-dimensional (2D slice). CONCLUSION: By combining the proposed Bloch simulation toolbox and the partially matrix-free reconstruction algorithm, 2D MR-STAT reconstructions can be performed in the order of one minute on a modern GPU card. The Bloch simulation toolbox can be utilized for other quantitative MRI techniques as well, for example for online dictionary generation for MR Fingerprinting.

Joint MAPLE: Accelerated joint T1 and T 2 * $$ {{\mathrm{T}}_2}

PURPOSE: Quantitative MRI finds important applications in clinical and research studies. However, it is encoding intensive and may suffer from prohibitively long scan times. Accelerated MR parameter mapping techniques have been developed to help address these challenges. Here, an accelerated joint T1, T 2 * $$ {{\mathrm{T}}_2}^{\ast } $$ , frequency and proton density mapping technique with scan-specific self-supervised network reconstruction is proposed to synergistically combine parallel imaging, model-based, and deep learning approaches to speed up parameter mapping. METHODS: Proposed framework, Joint MAPLE, includes parallel imaging, signal modeling, and data consistency blocks which are optimized jointly in a combined loss function. A scan-specific self-supervised reconstruction is embedded into the framework, which takes advantage of multi-contrast data from a multi-echo, multi-flip angle, gradient echo acquisition. RESULTS: In comparison with parallel reconstruction techniques powered by low-rank methods, emerging scan specific networks, and model-based T 2 * $$ {{\mathrm{T}}_2}^{\ast } $$ estimation approaches, the proposed framework reduces the reconstruction error in parameter maps by approximately two-fold on average at acceleration rates as high as R = 16 with uniform sampling. It can outperform evaluated parallel reconstruction techniques up to four-fold on average in the presence of challenging sub-sampling masks. It is observed that Joint MAPLE performs well at extreme acceleration rates of R = 25 and R = 36 with error values less than 20%. CONCLUSION: Joint MAPLE enables higher fidelity parameter estimation at high acceleration rates by synergistically combining parallel imaging and model-based parameter mapping and exploiting multi-echo, multi-flip angle datasets. Utilizing a scan-specific self-supervised reconstruction obviates the need for large data sets for training while improving the parameter estimation ability.

Tensor denoising of quantitative multi-parameter mapping.

PURPOSE: Quantitative multi-parameter mapping (MPM) provides maps of physical quantities representing physiologically meaningful tissue characteristics, which allows to investigate microstructure-function relationships reflecting normal or pathologic processes in the brain. However, the achievable spatial resolution and stability of MPM for basic research or clinical applications is severely constrained by SNR limits of the MR acquisition process, resulting in relatively long acquisition times. To increase SNR, we denoise MPM acquisitions using principal component analysis along tensors exploiting the Marchenko-Pastur law (tMPPCA). METHODS: tMPPCA denoising was applied to three sets of MPM raw data before the quantification of maps of proton density, magnetization transfer saturation, R1, and R2*. The regional SNR gain for high-resolution MPM was investigated as well as reproducibility gains for clinically optimized protocols with moderate and high acceleration factors at different image resolutions. RESULTS: Substantial noise reduction in raw data was achieved, resulting in reduced noise for quantitative mapping up to sixfold without introducing bias of mean values (below 1%). Scan-rescan fluctuations were reduced up to threefold. Denoising allowed to decrease the voxel volume fourfold at the same scan time or reduce the scan time twofold at same voxel volume without loss of sensitivity. CONCLUSIONS: tMPPCA denoising can (a) improve of fine spatial and temporal patterns, (b) considerably reduce scan time for clinical applications, or (c) increase resolution to potentially push cutting-edge MPM protocols from the upper to the lower limit of the mesoscopic scale.

Bias-reduced neural networks for parameter estimation in quantitative MRI.

PURPOSE: To develop neural network (NN)-based quantitative MRI parameter estimators with minimal bias and a variance close to the Cramér-Rao bound. THEORY AND METHODS: We generalize the mean squared error loss to control the bias and variance of the NN's estimates, which involves averaging over multiple noise realizations of the same measurements during training. Bias and variance properties of the resulting NNs are studied for two neuroimaging applications. RESULTS: In simulations, the proposed strategy reduces the estimates' bias throughout parameter space and achieves a variance close to the Cramér-Rao bound. In vivo, we observe good concordance between parameter maps estimated with the proposed NNs and traditional estimators, such as nonlinear least-squares fitting, while state-of-the-art NNs show larger deviations. CONCLUSION: The proposed NNs have greatly reduced bias compared to those trained using the mean squared error and offer significantly improved computational efficiency over traditional estimators with comparable or better accuracy.

Repeat it without me: Crowdsourcing the T1 mapping common ground via the ISMRM reproducibility challenge.

PURPOSE: T1 mapping is a widely used quantitative MRI technique, but its tissue-specific values remain inconsistent across protocols, sites, and vendors. The ISMRM Reproducible Research and Quantitative MR study groups jointly launched a challenge to assess the reproducibility of a well-established inversion-recovery T1 mapping technique, using acquisition details from a seminal T1 mapping paper on a standardized phantom and in human brains. METHODS: The challenge used the acquisition protocol from Barral et al. (2010). Researchers collected T1 mapping data on the ISMRM/NIST phantom and/or in human brains. Data submission, pipeline development, and analysis were conducted using open-source platforms. Intersubmission and intrasubmission comparisons were performed. RESULTS: Eighteen submissions (39 phantom and 56 human datasets) on scanners by three MRI vendors were collected at 3 T (except one, at 0.35 T). The mean coefficient of variation was 6.1% for intersubmission phantom measurements, and 2.9% for intrasubmission measurements. For humans, the intersubmission/intrasubmission coefficient of variation was 5.9/3.2% in the genu and 16/6.9% in the cortex. An interactive dashboard for data visualization was also developed: https://rrsg2020.dashboards.neurolibre.org. CONCLUSION: The T1 intersubmission variability was twice as high as the intrasubmission variability in both phantoms and human brains, indicating that the acquisition details in the original paper were insufficient to reproduce a quantitative MRI protocol. This study reports the inherent uncertainty in T1 measures across independent research groups, bringing us one step closer to a practical clinical baseline of T1 variations in vivo.

Fast and High-Resolution T2 Mapping Based on Echo Merging Plus k-t Undersampling with Reduced Refocusing Flip Angles (TEMPURA) as Methods for Human Renal MRI.

PURPOSE: To develop a highly accelerated multi-echo spin-echo method, TEMPURA, for reducing the acquisition time and/or increasing spatial resolution for kidney T2 mapping. METHODS: TEMPURA merges several adjacent echoes into one k-space by either combining independent echoes or sharing one echo between k-spaces. The combined k-space is reconstructed based on compressed sensing theory. Reduced flip angles are used for the refocusing pulses, and the extended phase graph algorithm is used to correct the effects of indirect echoes. Two sequences were developed: a fast breath-hold sequence; and a high-resolution sequence. The performance was evaluated prospectively on a phantom, 16 healthy subjects, and two patients with different types of renal tumors. RESULTS: The fast TEMPURA method reduced the acquisition time from 3-5 min to one breath-hold (18 s). Phantom measurements showed that fast TEMPURA had a mean absolute percentage error (MAPE) of 8.2%, which was comparable to a standardized respiratory-triggered sequence (7.4%), but much lower than a sequence accelerated by purely k-t undersampling (21.8%). High-resolution TEMPURA reduced the in-plane voxel size from 3 × 3 to 1 × 1 mm2, resulting in improved visualization of the detailed anatomical structure. In vivo T2 measurements demonstrated good agreement (fast: MAPE = 1.3%-2.5%; high-resolution: MAPE = 2.8%-3.3%) and high correlation coefficients (fast: R = 0.85-0.98; high-resolution: 0.82-0.96) with the standardized method, outperforming k-t undersampling alone (MAPE = 3.3-4.5%, R = 0.57-0.59). CONCLUSION: TEMPURA provides fast and high-resolution renal T2 measurements. It has the potential to improve clinical throughput and delineate intratumoral heterogeneity and tissue habitats at unprecedented spatial resolution.

Zero-DeepSub: Zero-shot deep subspace reconstruction for rapid multiparametric quantitative MRI using 3D-QALAS.

PURPOSE: To develop and evaluate methods for (1) reconstructing 3D-quantification using an interleaved Look-Locker acquisition sequence with T2 preparation pulse (3D-QALAS) time-series images using a low-rank subspace method, which enables accurate and rapid T1 and T2 mapping, and (2) improving the fidelity of subspace QALAS by combining scan-specific deep-learning-based reconstruction and subspace modeling. THEORY AND METHODS: A low-rank subspace method for 3D-QALAS (i.e., subspace QALAS) and zero-shot deep-learning subspace method (i.e., Zero-DeepSub) were proposed for rapid and high fidelity T1 and T2 mapping and time-resolved imaging using 3D-QALAS. Using an ISMRM/NIST system phantom, the accuracy and reproducibility of the T1 and T2 maps estimated using the proposed methods were evaluated by comparing them with reference techniques. The reconstruction performance of the proposed subspace QALAS using Zero-DeepSub was evaluated in vivo and compared with conventional QALAS at high reduction factors of up to nine-fold. RESULTS: Phantom experiments showed that subspace QALAS had good linearity with respect to the reference methods while reducing biases and improving precision compared to conventional QALAS, especially for T2 maps. Moreover, in vivo results demonstrated that subspace QALAS had better g-factor maps and could reduce voxel blurring, noise, and artifacts compared to conventional QALAS and showed robust performance at up to nine-fold acceleration with Zero-DeepSub, which enabled whole-brain T1, T2, and PD mapping at 1 mm isotropic resolution within 2 min of scan time. CONCLUSION: The proposed subspace QALAS along with Zero-DeepSub enabled high fidelity and rapid whole-brain multiparametric quantification and time-resolved imaging.

DeepEMC-T2 mapping: Deep learning-enabled T2 mapping based on echo modulation curve modeling.

PURPOSE: Echo modulation curve (EMC) modeling enables accurate quantification of T2 relaxation times in multi-echo spin-echo (MESE) imaging. The standard EMC-T2 mapping framework, however, requires sufficient echoes and cumbersome pixel-wise dictionary-matching steps. This work proposes a deep learning version of EMC-T2 mapping, called DeepEMC-T2 mapping, to efficiently estimate accurate T2 maps from fewer echoes. METHODS: DeepEMC-T2 mapping was developed using a modified U-Net to estimate both T2 and proton density (PD) maps directly from MESE images. The network implements several new features to improve the accuracy of T2/PD estimation. A total of 67 MESE datasets acquired in axial orientation were used for network training and evaluation. An additional 57 datasets acquired in coronal orientation with different scan parameters were used to evaluate the generalizability of the framework. The performance of DeepEMC-T2 mapping was evaluated in seven experiments. RESULTS: Compared to the reference, DeepEMC-T2 mapping achieved T2 estimation errors from 1% to 11% and PD estimation errors from 0.4% to 1.5% with ten/seven/five/three echoes, which are more accurate than standard EMC-T2 mapping. By incorporating datasets acquired with different scan parameters and orientations for joint training, DeepEMC-T2 exhibits robust generalizability across varying imaging protocols. Increasing the echo spacing and including longer echoes improve the accuracy of parameter estimation. The new features proposed in DeepEMC-T2 mapping all enabled more accurate T2 estimation. CONCLUSIONS: DeepEMC-T2 mapping enables simplified, efficient, and accurate T2 quantification directly from MESE images without dictionary matching. Accurate T2 estimation from fewer echoes allows for increased volumetric coverage and/or higher slice resolution without prolonging total scan times.

Self-supervised learning for denoising of multidimensional MRI data.

PURPOSE: To develop a fast denoising framework for high-dimensional MRI data based on a self-supervised learning scheme, which does not require ground truth clean image. THEORY AND METHODS: Quantitative MRI faces limitations in SNR, because the variation of signal amplitude in a large set of images is the key mechanism for quantification. In addition, the complex non-linear signal models make the fitting process vulnerable to noise. To address these issues, we propose a fast deep-learning framework for denoising, which efficiently exploits the redundancy in multidimensional MRI data. A self-supervised model was designed to use only noisy images for training, bypassing the challenge of clean data paucity in clinical practice. For validation, we used two different datasets of simulated magnetization transfer contrast MR fingerprinting (MTC-MRF) dataset and in vivo DWI image dataset to show the generalizability. RESULTS: The proposed method drastically improved denoising performance in the presence of mild-to-severe noise regardless of noise distributions compared to previous methods of the BM3D, tMPPCA, and Patch2self. The improvements were even pronounced in the following quantification results from the denoised images. CONCLUSION: The proposed MD-S2S (Multidimensional-Self2Self) denoising technique could be further applied to various multi-dimensional MRI data and improve the quantification accuracy of tissue parameter maps.

Measurement of T1ρ dispersion with compressed sensing and magnetization prepared radial balanced steady-state free precession in spontaneous human osteoarthritis.

PURPOSE: To assess the potential for accelerating continuous-wave (CW) T1ρ dispersion measurement with compressed sensing approach via studying the effect that the data reduction has on the ability to detect differences between intact and degenerated articular cartilage with different spin-lock amplitudes and to assess quantitative bias due to acceleration. METHODS: Osteochondral plugs (n = 27, 4 mm diameter) from femur (n = 14) and tibia (n = 13) regions from human cadaver knee joints were obtained from commercial biobank (Science Care, USA) under Ethical permission 134/2015. MRI of specimens was performed at 9.4T with magnetization prepared radial balanced SSFP (bSSFP) readout sequence, and the CWT1ρ relaxation time maps were computed from the measured data. The relaxation time maps were evaluated in the cartilage zones for different acceleration factors. For reference, Osteoarthritis Research Society International (OARSI) grading and biomechanical measurements were performed and correlated with the MRI findings. RESULTS: Four-fold acceleration of CWT1ρ dispersion measurement by compressed sensing approach was feasible without meaningful loss in the sensitivity to osteoarthritic (OA) changes within the articular cartilage. Differences were significant between intact and OA groups in the superficial and transitional zones, and CWT1ρ correlated moderately with the reference measurements (0.3 < r < 0.7). CONCLUSION: CWT1ρ was able to differentiate between intact and OA cartilage even with four-fold acceleration. This indicates that acceleration of CWT1ρ dispersion measurement by compressed sensing approach is feasible with negligible loss in the sensitivity to osteoarthritic changes in articular cartilage.

Mitigating slice cross-talk in multi-slice multi-echo spin echo T2 mapping.

PURPOSE: To investigate whether a T2 inter-slice variation could occur when a multi-slice multi-echo spin echo (MESE) sequence is used for image acquisition and to propose an enhanced method for reconstructing T2 maps that can effectively address and correct these variations. METHODS: Bloch simulations were performed accounting for the direct saturation effect to evaluate magnetization changes in multi-slice 2D MESE sequence. Experimental phantom scans were performed to validate these simulations. An improved version of the dictionary-based reconstruction approach was proposed, enabling the creation of a multi-slice dictionary of echo modulation curves (EMC). The corresponding method has been assayed considering inter-slice T2 variation with phantoms and in lower leg. RESULTS: Experimental and numerical study illustrate that direct saturation leads to a bias of EMCs. This bias during the T2 maps reconstructions using original single-slice EMC-dictionary method led to inter-slice T2 variation of 2.03% in average coefficient of variation (CV) in agarose phantoms, and up to 2.8% in vivo (for TR = 2 s, slice gap = 0%). A reduction of CV was observed when increasing the gap up to 100% (0.36% in phantoms, and up to 1.5% in vivo) or increasing TR up to 4 s (0.76% in phantoms, and up to 1.9% in vivo). Matching the multi-slice experimental data with multi-slice dictionaries provided a reduced CV of 0.54% in phantoms and up to 2.3% in vivo. CONCLUSION: T2 values quantified from multi-slice MESE images using single-slice dictionaries are biased. A dedicated multi-slice EMC method providing the appropriate dictionaries can reduce the inter-slice T2 variation.