Pseudo Test-Retest Evaluation of Millimeter-Resolution Whole-Brain Dynamic Contrast-enhanced MRI in Patients with High-Grade Glioma.

Background Advances in sub-Nyquist-sampled dynamic contrast-enhanced (DCE) MRI enable monitoring of brain tumors with millimeter resolution and whole-brain coverage. Such undersampled quantitative methods need careful characterization regarding achievable test-retest reproducibility. Purpose To demonstrate a fully automated high-resolution whole-brain DCE MRI pipeline with 30-fold sparse undersampling and estimate its reproducibility on the basis of reference regions of stable tissue types during multiple posttreatment time points by using longitudinal clinical images of high-grade glioma. Materials and Methods Two methods for sub-Nyquist-sampled DCE MRI were extended with automatic estimation of vascular input functions. Continuously acquired three-dimensional k-space data with ramped-up flip angles were partitioned to yield high-resolution, whole-brain tracer kinetic parameter maps with matched precontrast-agent T1 and M0 maps. Reproducibility was estimated in a retrospective study in participants with high-grade glioma, who underwent three consecutive standard-of-care examinations between December 2016 and April 2019. Coefficients of variation and reproducibility coefficients were reported for histogram statistics of the tracer kinetic parameters plasma volume fraction and volume transfer constant (Ktrans) on five healthy tissue types. Results The images from 13 participants (mean age ± standard deviation, 61 years ± 10; nine women) with high-grade glioma were evaluated. In healthy tissues, the protocol achieved a coefficient of variation less than 57% for median Ktrans, if Ktrans was estimated consecutively. The maximum reproducibility coefficient for median Ktrans was estimated to be at 0.06 min-1 for large or low-enhancing tissues and to be as high as 0.48 min-1 in smaller or strongly enhancing tissues. Conclusion A fully automated, sparsely sampled DCE MRI reconstruction with patient-specific vascular input function offered high spatial and temporal resolution and whole-brain coverage; in healthy tissues, the protocol estimated median volume transfer constant with maximum reproducibility coefficient of 0.06 min-1 in large, low-enhancing tissue regions and maximum reproducibility coefficient of less than 0.48 min-1 in smaller or more strongly enhancing tissue regions. Published under a CC BY 4.0 license. Online supplemental material is available for this article. See also the editorial by Lenkinski in this issue.

Sparse precontrast T1 mapping for high-resolution whole-brain DCE-MRI.

PURPOSE: To develop and evaluate an efficient precontrast T1 mapping technique suitable for quantitative high-resolution whole-brain dynamic contrast-enhanced-magnetic resonance imaging (DCE-MRI). METHODS: Variable flip angle (VFA) T1 mapping was considered that provides 1 × 1 × 2 mm3 resolution to match a recent high-resolution whole-brain DCE-MRI protocol. Seven FAs were logarithmically spaced from 1.5° to 15°. T1 and M0 maps were estimated using model-based reconstruction. This approach was evaluated using an anatomically realistic brain tumor digital reference object (DRO) with noise-mimicking 3T neuroimaging and fully sampled data acquired from one healthy volunteer. Methods were also applied on fourfold prospectively undersampled VFA data from 13 patients with high-grade gliomas. RESULTS: T1 -mapping precision decreased with undersampling factor R, althoughwhereas bias remained small before a critical R. In the noiseless DRO, T1 bias was <25 ms in white matter (WM) and <11 ms in brain tumor (BT). T1 standard deviation (SD) was <119.5 ms in WM (coefficient of variation [COV] ~11.0%) and <253.2 ms in BT (COV ~12.7%). In the noisy DRO, T1 bias was <50 ms in WM and <30 ms in BT. For R ≤ 10, T1 SD was <107.1 ms in WM (COV ~9.9%) and <240.9 ms in BT (COV ~12.1%). In the healthy subject, T1 bias was <30 ms for R ≤ 16. At R = 4, T1 SD was 171.4 ms (COV ~13.0%). In the prospective brain tumor study, T1 values were consistent with literature values in WM and BT. CONCLUSION: High-resolution whole-brain VFA T1 mapping is feasible with sparse sampling, supporting its use for quantitative DCE-MRI.

Impact of (k,t) sampling on DCE MRI tracer kinetic parameter estimation in digital reference objects.

PURPOSE: To evaluate the impact of (k,t) data sampling on the variance of tracer-kinetic parameter (TK) estimation in high-resolution whole-brain dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) using digital reference objects. We study this in the context of TK model constraints, and in the absence of other constraints. METHODS: Three anatomically and physiologically realistic brain-tumor digital reference objects were generated. Data sampling strategies included uniform and variable density; zone-based, lattice, pseudo-random, and pseudo-radial; with 50-time frames and 4-fold to 25-fold undersampling. In all cases, we assume a fully sampled first time frame, and prior knowledge of the arterial input function. TK parameters were estimated by indirect estimation (i.e., image-time-series reconstruction followed by model fitting), and direct estimation from the under-sampled data. We evaluated methods based on the Cramér-Rao bound and Monte-Carlo simulations, over the range of signal-to-noise ratio (SNR) seen in clinical brain DCE-MRI. RESULTS: Lattice-based sampling provided the lowest SDs, followed by pseudo-random, pseudo-radial, and zone-based. This ranking was consistent for the Patlak and extended Tofts model. Pseudo-random sampling resulted in 19% higher averaged SD compared to lattice-based sampling. Zone-based sampling resulted in substantially higher SD at undersampling factors above 10. CRB analysis showed only a small difference between uniform and variable density for both lattice-based and pseudo-random sampling up to undersampling factors of 25. CONCLUSION: Lattice sampling provided the lowest SDs, although the differences between sampling schemes were not substantial at low undersampling factors. The differences between lattice-based and pseudo-random sampling strategies with both uniform and variable density were within the range of error induced by other sources, at up to 25-fold undersampling.

Deblurring for spiral real-time MRI using convolutional neural networks.

PURPOSE: To develop and evaluate a fast and effective method for deblurring spiral real-time MRI (RT-MRI) using convolutional neural networks. METHODS: We demonstrate a 3-layer residual convolutional neural networks to correct image domain off-resonance artifacts in speech production spiral RT-MRI without the knowledge of field maps. The architecture is motivated by the traditional deblurring approaches. Spatially varying off-resonance blur is synthetically generated by using discrete object approximation and field maps with data augmentation from a large database of 2D human speech production RT-MRI. The effect of off-resonance range, shift-invariance of blur, and readout durations on deblurring performance are investigated. The proposed method is validated using synthetic and real data with longer readouts, quantitatively using image quality metrics and qualitatively via visual inspection, and with a comparison to conventional deblurring methods. RESULTS: Deblurring performance was found superior to a current autocalibrated method for in vivo data and only slightly worse than an ideal reconstruction with perfect knowledge of the field map for synthetic test data. Convolutional neural networks deblurring made it possible to visualize articulator boundaries with readouts up to 8 ms at 1.5 T, which is 3-fold longer than the current standard practice. The computation time was 12.3 ± 2.2 ms per frame, enabling low-latency processing for RT-MRI applications. CONCLUSION: Convolutional neural networks deblurring is a practical, efficient, and field map-free approach for the deblurring of spiral RT-MRI. In the context of speech production imaging, this can enable 1.7-fold improvement in scan efficiency and the use of spiral readouts at higher field strengths such as 3 T.

Joint arterial input function and tracer kinetic parameter estimation from undersampled dynamic contrast-enhanced MRI using a model consistency constraint.

PURPOSE: To develop and evaluate a model-based reconstruction framework for joint arterial input function (AIF) and kinetic parameter estimation from undersampled brain tumor dynamic contrast-enhanced MRI (DCE-MRI) data. METHODS: The proposed method poses the tracer-kinetic (TK) model as a model consistency constraint, enabling the flexible inclusion of different TK models and TK solvers, and the joint estimation of the AIF. The proposed method is evaluated using an anatomic realistic digital reference object (DRO), and nine retrospectively down-sampled brain tumor DCE-MRI datasets. We also demonstrate application to 30-fold prospectively undersampled brain tumor DCE-MRI. RESULTS: In DRO studies with up to 60-fold undersampling, the proposed method provided TK maps with low error that were comparable to fully sampled data and were demonstrated to be compatible with a third-party TK solver. In retrospective undersampling studies, this method provided patient-specific AIF with normalized root mean-squared-error (normalized by the 90th percentile value) less than 8% at up to 100-fold undersampling. In the 30-fold undersampled prospective study, the proposed method provided high-resolution whole-brain TK maps and patient-specific AIF. CONCLUSION: The proposed model-based DCE-MRI reconstruction enables the use of different TK solvers with a model consistency constraint and enables joint estimation of patient-specific AIF. TK maps and patient-specific AIF with high fidelity can be reconstructed at up to 100-fold undersampling in k,t-space. Magn Reson Med 79:2804-2815, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

A multispeaker dataset of raw and reconstructed speech production real-time MRI video and 3D volumetric images.

Real-time magnetic resonance imaging (RT-MRI) of human speech production is enabling significant advances in speech science, linguistics, bio-inspired speech technology development, and clinical applications. Easy access to RT-MRI is however limited, and comprehensive datasets with broad access are needed to catalyze research across numerous domains. The imaging of the rapidly moving articulators and dynamic airway shaping during speech demands high spatio-temporal resolution and robust reconstruction methods. Further, while reconstructed images have been published, to-date there is no open dataset providing raw multi-coil RT-MRI data from an optimized speech production experimental setup. Such datasets could enable new and improved methods for dynamic image reconstruction, artifact correction, feature extraction, and direct extraction of linguistically-relevant biomarkers. The present dataset offers a unique corpus of 2D sagittal-view RT-MRI videos along with synchronized audio for 75 participants performing linguistically motivated speech tasks, alongside the corresponding public domain raw RT-MRI data. The dataset also includes 3D volumetric vocal tract MRI during sustained speech sounds and high-resolution static anatomical T2-weighted upper airway MRI for each participant.

Efficient DCE-MRI Parameter and Uncertainty Estimation Using a Neural Network.

Quantitative DCE-MRI provides voxel-wise estimates of tracer-kinetic parameters that are valuable in the assessment of health and disease. These maps suffer from many known sources of variability. This variability is expensive to compute using current methods, and is typically not reported. Here, we demonstrate a novel approach for simultaneous estimation of tracer-kinetic parameters and their uncertainty due to intrinsic characteristics of the tracer-kinetic model, with very low computation time. We train and use a neural network to estimate the approximate joint posterior distribution of tracer-kinetic parameters. Uncertainties are estimated for each voxel and are specific to the patient, exam, and lesion. We demonstrate the methods' ability to produce accurate tracer-kinetic maps. We compare predicted parameter ranges with uncertainties introduced by noise and by differences in post-processing in a digital reference object. The predicted parameter ranges correlate well with tracer-kinetic parameter ranges observed across different noise realizations and regression algorithms. We also demonstrate the value of this approach to differentiate significant from insignificant changes in brain tumor pharmacokinetics over time. This is achieved by enforcing consistency in resolving model singularities in the applied tracer-kinetic model.

Tracer kinetic models as temporal constraints during brain tumor DCE-MRI reconstruction.

PURPOSE: To apply tracer kinetic models as temporal constraints during reconstruction of under-sampled brain tumor dynamic contrast enhanced (DCE) magnetic resonance imaging (MRI). METHODS: A library of concentration vs time profiles is simulated for a range of physiological kinetic parameters. The library is reduced to a dictionary of temporal bases, where each profile is approximated by a sparse linear combination of the bases. Image reconstruction is formulated as estimation of concentration profiles and sparse model coefficients with a fixed sparsity level. Simulations are performed to evaluate modeling error, and error statistics in kinetic parameter estimation in presence of noise. Retrospective under-sampling experiments are performed on a brain tumor DCE digital reference object (DRO), and 12 brain tumor in-vivo 3T datasets. The performances of the proposed under-sampled reconstruction scheme and an existing compressed sensing-based temporal finite-difference (tFD) under-sampled reconstruction were compared against the fully sampled inverse Fourier Transform-based reconstruction. RESULTS: Simulations demonstrate that sparsity levels of 2 and 3 model the library profiles from the Patlak and extended Tofts-Kety (ETK) models, respectively. Noise sensitivity analysis showed equivalent kinetic parameter estimation error statistics from noisy concentration profiles, and model approximated profiles. DRO-based experiments showed good fidelity in recovery of kinetic maps from 20-fold under-sampled data. In-vivo experiments demonstrated reduced bias and uncertainty in kinetic mapping with the proposed approach compared to tFD at under-sampled reduction factors >= 20. CONCLUSIONS: Tracer kinetic models can be applied as temporal constraints during brain tumor DCE-MRI reconstruction. The proposed under-sampled scheme resulted in model parameter estimates less biased with respect to conventional fully sampled DCE MRI reconstructions and parameter estimation. The approach is flexible, can use nonlinear kinetic models, and does not require tuning of regularization parameters.

Radiofrequency transmit calibration: A multi-center evaluation of vendor-provided radiofrequency transmit mapping methods.

PURPOSE: To determine the accuracy and test-retest repeatability of fast radiofrequency (RF) transmit measurement approaches used in Dynamic Contrast Enhanced Magnetic Resonance Imaging (DCE-MRI). Spatial variation in the transmitted RF field introduces bias and increased variance in quantitative DCE-MRI metrics including tracer kinetic parameter maps. If unaccounted for, these errors can dominate all other sources of bias and variance. The amount and pattern of variation depend on scanner-specific hardware and software. METHODS: Human tissue mimicking torso and brain phantoms were constructed. RF transmit maps were measured and compared across eight different commercial scanners, from three major vendors, and three clinical sites. Vendor-recommended rapid methods for RF mapping were compared to a slower reference method. Imaging was repeated at all sites after 2 months. Ranges and magnitude of RF inhomogeneity were compared scanner-wise at two time points. Limits of Agreement of vendor-recommended methods and double-angle reference method were assessed. RESULTS: At 3 T, B1 + inhomogeneity spans across 35% in the head and 120% in the torso. Fast vendor provided methods are within 30% agreement with the reference double angle method for both the head and the torso phantom. CONCLUSIONS: If unaccounted for, B1 + inhomogeneity can severely impact tracer-kinetic parameter estimation. Depending on the scanner, fast vendor provided B1 + mapping sequences allow unbiased and reproducible measurements of B1 + inhomogeneity to correct for this source of bias.