Accelerated 3D whole-brain T1, T2, and proton density mapping: feasibility for clinical glioma MR imaging

Physics-guided self-supervised learning for retrospective T1 and T2 mapping from conventional weighted brain MRI: Technical developments and initial validation in glioblastoma

PURPOSE

To develop a self-supervised learning method to retrospectively estimate T1 and T2 values from clinical weighted MRI.

METHODS

A self-supervised learning approach was constructed to estimate T1, T2, and proton density maps from conventional T1- and T2-weighted images. MR physics models were employed to regenerate the weighted images from the network outputs, and the network was optimized based on loss calculated between the synthesized and input weighted images, alongside additional constraints based on prior information. The method was evaluated on healthy volunteer data, with conventional mapping as references. The reproducibility was examined on two 3.0T scanners. Performance in tumor characterization was inspected by applying the method to a public glioblastoma dataset.

RESULTS

For T1 and T2 estimation from three weighted images (T1 MPRAGE, T1 gradient echo sequences, and T2 turbo spin echo), the deep learning method achieved global voxel-wise error ≤9% in brain parenchyma and regional error ≤12.2% in six types of brain tissues. The regional measurements obtained from two scanners showed mean differences ≤2.4% and correlation coefficients >0.98, demonstrating excellent reproducibility. In the 50 glioblastoma patients, the retrospective quantification results were in line with literature reports from prospective methods, and the T2 values were found to be higher in tumor regions, with sensitivity of 0.90 and specificity of 0.92 in a voxel-wise classification task between normal and abnormal regions.

CONCLUSION

The self-supervised learning method is promising for retrospective T1 and T2 quantification from clinical MR images, with the potential to improve the availability of quantitative MRI and facilitate brain tumor characterization.

Magnetic resonance relaxometry in quantitative imaging of brain gliomas: A literature review

Magnetic resonance (MR) relaxometry is a quantitative imaging method that measures tissue relaxation properties. This review discusses the state of the art of clinical proton MR relaxometry for glial brain tumors. Current MR relaxometry technology also includes MR fingerprinting and synthetic MRI, which solve the inefficiencies and challenges of earlier techniques. Despite mixed results regarding its capability for brain tumor differential diagnosis, there is growing evidence that MR relaxometry can differentiate between gliomas and metastases and between glioma grades. Studies of the peritumoral zones have demonstrated their heterogeneity and possible directions of tumor infiltration. In addition, relaxometry offers T2* mapping that can define areas of tissue hypoxia not discriminated by perfusion assessment. Studies of tumor therapy response have demonstrated an association between survival and progression terms and dynamics of native and contrast-enhanced tumor relaxometric profiles. In conclusion, MR relaxometry is a promising technique for glial tumor diagnosis, particularly in association with neuropathological studies and other imaging techniques.

Automatic Brain Tissue and Lesion Segmentation and Multi-Parametric Mapping of Contrast-Enhancing Gliomas without the Injection of Contrast Agents: A Preliminary Study

This study aimed to develop a rapid, 1 mm3 isotropic resolution, whole-brain MRI technique for automatic lesion segmentation and multi-parametric mapping without using contrast by continuously applying balanced steady-state free precession with inversion pulses throughout incomplete inversion recovery in a single 6 min scan. Modified k-means clustering was performed for automatic brain tissue and lesion segmentation using distinct signal evolutions that contained mixed T1/T2/magnetization transfer properties. Multi-compartment modeling was used to derive quantitative multi-parametric maps for tissue characterization. Fourteen patients with contrast-enhancing gliomas were scanned with this sequence prior to the injection of a contrast agent, and their segmented lesions were compared to conventionally defined manual segmentations of T2-hyperintense and contrast-enhancing lesions. Simultaneous T1, T2, and macromolecular proton fraction maps were generated and compared to conventional 2D T1 and T2 mapping and myelination water fraction mapping acquired with MAGiC. The lesion volumes defined with the new method were comparable to the manual segmentations (r = 0.70, p < 0.01; t-test p > 0.05). The T1, T2, and macromolecular proton fraction mapping values of the whole brain were comparable to the reference values and could distinguish different brain tissues and lesion types (p < 0.05), including infiltrating tumor regions within the T2-lesion. Highly efficient, whole-brain, multi-contrast imaging facilitated automatic lesion segmentation and quantitative multi-parametric mapping without contrast, highlighting its potential value in the clinic when gadolinium is contraindicated.

A three-dimensional Magnetic Resonance Spin Tomography in Time-domain protocol for high-resolution multiparametric quantitative magnetic resonance imaging

Magnetic Resonance Spin TomogrAphy in Time-domain (MR-STAT) is a multiparametric quantitative MR framework, which allows for simultaneously acquiring quantitative tissue parameters such as T1, T2, and proton density from one single short scan. A typical two-dimensional (2D) MR-STAT acquisition uses a gradient-spoiled, gradient-echo sequence with a slowly varying RF flip-angle train and Cartesian readouts, and the quantitative tissue maps are reconstructed by an iterative, model-based optimization algorithm. In this work, we design a three-dimensional (3D) MR-STAT framework based on previous 2D work, in order to achieve better image signal-to-noise ratio, higher though-plane resolution, and better tissue characterization. Specifically, we design a 7-min, high-resolution 3D MR-STAT sequence, and the corresponding two-step reconstruction algorithm for the large-scale dataset. To reduce the long acquisition time, Cartesian undersampling strategies such as SENSE are adopted in our transient-state quantitative framework. To reduce the computational burden, a data-splitting scheme is designed for decoupling the 3D reconstruction problem into independent 2D reconstructions. The proposed 3D framework is validated by numerical simulations, phantom experiments, and in vivo experiments. High-quality knee quantitative maps with 0.8 × 0.8 × 1.5 mm3 resolution and bilateral lower leg maps with 1.6 mm isotropic resolution can be acquired using the proposed 7-min acquisition sequence and the 3-min-per-slice decoupled reconstruction algorithm. The proposed 3D MR-STAT framework could have wide clinical applications in the future.

Investigations on correlations between changes of optical tissue properties and NMR relaxation times

BACKGROUND

Accurate light dosimetry is a complex remaining challenge in interstitial photodynamic therapy (iPDT) for malignant gliomas. The light dosimetry should ideally be based on the tissue morphology and the individual optical tissue properties of each tissue type in the target region. First investigations are reported on using NMR information to estimate changes of individual optical tissue properties.

METHODS

Porcine brain tissue and optical tissue phantoms were investigated. To the porcine brain, supplements were added to simulate an edema or high blood content. The tissue phantoms were based on agar, Lipoveneous, ink, blood and gadobutrol (Gd-based MRI contrast agent). The concentrations of phantom ingredients and tissue additives are varied to compare concentration-dependent effects on optical and NMR properties. A 3-tesla whole-body MRI system was used to determine T1 and T2 relaxation times. Optical tissue properties, i.e., the spectrally resolved absorption and reduced scattering coefficient, were obtained using a single integrating sphere setup. The observed changes of NMR and optical properties were compared to each other.

RESULTS

By adjusting the NMR relaxation times and optical tissue properties of the tissue phantoms to literature values, recipes for human brain tumor, white matter and grey matter tissue phantoms were obtained that mimic these brain tissues simultaneously in both properties. For porcine brain tissue, it was observed that with increasing water concentration in the tissue, both NMR-relaxation times increased, while µa decreased and µs' increased at 635 nm. The addition of blood to porcine brain samples showed a constant T1, while T2 shortened and the absorption coefficient at 635 nm increased.

CONCLUSIONS

In this investigation, by changing sample contents, notable changes of both NMR relaxation times and optical tissue properties have been observed and their relations examined. The developed dual NMR/optical tissue phantoms can be used in iPDT research, clinical training and demonstrations.

Tailored magnetic resonance fingerprinting of post-operative pediatric brain tumor patients

PURPOSE

Brain and spinal cord tumors are the second most common cancer in children and account for one out of four cancers diagnosed. However, the long acquisition times associated with acquiring both data types prohibit using quantitative MR (qMR) in pediatric imaging protocols. This study aims to demonstrate the tailored magnetic resonance fingerprinting's (TMRF) ability to simultaneously provide quantitative maps (T1, T2) and multi-contrast qualitative images (T1 weighted, T1 FLAIR, T2 weighted) rapidly in pediatric brain tumor patients.

METHODS

In this work, we imaged five pediatric patients with brain tumors (resected/residual) using TMRF at 3 T. We compared the TMRF-derived T2 weighted images with those from the vendor-supplied sequence (as the gold standard, GS) for healthy and pathological tissue signal intensities. The relaxometric maps from TMRF were subjected to a region of interest (ROI) analysis to differentiate between healthy and pathological tissues. We performed the Wilcoxon rank sum test to check for significant differences between the two tissue types.

RESULTS

We found significant differences (p < 0.05) in both T1 and T2 ROI values between the two tissue types. A strong correlation was found between the TMRF-based T2 weighted and GS signal intensities for the healthy (correlation coefficient, r = 0.99) and pathological tissues (r = 0.88).

CONCLUSION

The TMRF implementation provides the two relaxometric maps and can potentially save ~2 min if it replaces the T2-weighted imaging in the current protocol.

Advanced MR Techniques for Preoperative Glioma Characterization: Part 1

Preoperative clinical magnetic resonance imaging (MRI) protocols for gliomas, brain tumors with dismal outcomes due to their infiltrative properties, still rely on conventional structural MRI, which does not deliver information on tumor genotype and is limited in the delineation of diffuse gliomas. The GliMR COST action wants to raise awareness about the state of the art of advanced MRI techniques in gliomas and their possible clinical translation or lack thereof. This review describes current methods, limits, and applications of advanced MRI for the preoperative assessment of glioma, summarizing the level of clinical validation of different techniques. In this first part, we discuss dynamic susceptibility contrast and dynamic contrast-enhanced MRI, arterial spin labeling, diffusion-weighted MRI, vessel imaging, and magnetic resonance fingerprinting. The second part of this review addresses magnetic resonance spectroscopy, chemical exchange saturation transfer, susceptibility-weighted imaging, MRI-PET, MR elastography, and MR-based radiomics applications. Evidence Level: 3 Technical Efficacy: Stage 2.

Evaluation of Gliomas with Magnetic Resonance Fingerprinting with PET Correlation-A Comparative Study

OBJECTIVES

Advanced MR imaging of brain tumors is still mainly based on qualitative imaging. PET imaging offers additive metabolic information, and MR fingerprinting (MRF) offers a novel approach to quantitative data acquisition. The purpose of this study was to evaluate the ability of MRF to predict tumor regions and grading in combination with PET.

METHODS

Seventeen patients with histologically verified infiltrating gliomas and available amino-acid PET data were enrolled. ROIs for solid tumor parts (SPo), perifocal edema (ED1), and normal-appearing white matter (NAWM) were selected on conventional MRI sequences and aligned to the MRF and PET images. The predictability of gliomas by region and grading as well as intermodal correlations were assessed.

RESULTS

For MRF, we calculated an overall predictability by region (SPo, ED1, and NAWM) for all of the MRF parameters of 76.5%, 47.1%, and 94.1%, respectively. The overall ability to distinguish low- from high-grade gliomas using MRF was 88.9% for LGG and 75% for HGG, with an accuracy of 82.4%, a ppV of 85.71%, and an npV of 80%. PET positivity was found in 13/17 patients for solid tumor parts, and in 3/17 patients for the edema region. However, there was no significant difference in region-specific MRF values between PET positive and PET negative patients.

CONCLUSIONS

MRF and PET provide quantitative measurements of the tumor tissue characteristics of gliomas, with good predictability. Nonetheless, the results are dissimilar, reflecting the different underlying mechanisms of each method.

Toward deep learning replacement of gadolinium in neuro-oncology: A review of contrast-enhanced synthetic MRI

Gadolinium-based contrast agents (GBCAs) have become a crucial part of MRI acquisitions in neuro-oncology for the detection, characterization and monitoring of brain tumors. However, contrast-enhanced (CE) acquisitions not only raise safety concerns, but also lead to patient discomfort, the need of more skilled manpower and cost increase. Recently, several proposed deep learning works intend to reduce, or even eliminate, the need of GBCAs. This study reviews the published works related to the synthesis of CE images from low-dose and/or their native -non CE- counterparts. The data, type of neural network, and number of input modalities for each method are summarized as well as the evaluation methods. Based on this analysis, we discuss the main issues that these methods need to overcome in order to become suitable for their clinical usage. We also hypothesize some future trends that research on this topic may follow.

Pre-contrast MAGiC in treated gliomas: a pilot study of quantitative MRI

Quantitative MR imaging is becoming more feasible to be used in clinical work since new approaches have been proposed in order to substantially accelerate the acquisition and due to the possibility of synthetically deriving weighted images from the parametric maps. However, their applicability has to be thoroughly validated in order to be included in clinical practice. In this pilot study, we acquired Magnetic Resonance Image Compilation scans to obtain T1, T2 and PD maps in 14 glioma patients. Abnormal tissue was segmented based on conventional images and using a deep learning segmentation technique to define regions of interest (ROIs). The quantitative T1, T2 and PD values inside ROIs were analyzed using the mean, the standard deviation, the skewness and the kurtosis and compared to the quantitative T1, T2 and PD values found in normal white matter. We found significant differences in pre-contrast T1 and T2 values between abnormal tissue and healthy tissue, as well as between T1w-enhancing and non-enhancing regions. ROC analysis was used to evaluate the potential of quantitative T1 and T2 values for voxel-wise classification of abnormal/normal tissue (AUC = 0.95) and of T1w enhancement/non-enhancement (AUC = 0.85). A cross-validated ROC analysis found high sensitivity (73%) and specificity (73%) with AUCs up to 0.68 on the a priori distinction between abnormal tissue with and without T1w-enhancement. These results suggest that normal tissue, abnormal tissue, and tissue with T1w-enhancement are distinguishable by their pre-contrast quantitative values but further investigation is needed.

Quantitative Relaxometry Metrics for Brain Metastases Compared to Normal Tissues: A Pilot MR Fingerprinting Study

The purpose of the present pilot study was to estimate T1 and T2 metric values derived simultaneously from a new, rapid Magnetic Resonance Fingerprinting (MRF) technique, as well as to assess their ability to characterize-brain metastases (BM) and normal-appearing brain tissues. Fourteen patients with BM underwent MRI, including prototype MRF, on a 3T scanner. In total, 108 measurements were analyzed: 42 from solid parts of BM's (21 each on T1 and T2 maps) and 66 from normal-appearing brain tissue (11 ROIs each on T1 and T2 maps for gray matter [GM], white matter [WM], and cerebrospinal fluid [CSF]). The BM's mean T1 and T2 values differed significantly from normal-appearing WM (p < 0.05). The mean T1 values from normal-appearing GM, WM, and CSF regions were 1205 ms, 840 ms, and 4233 ms, respectively. The mean T2 values were 108 ms, 78 ms, and 442 ms, respectively. The mean T1 and T2 values for untreated BM (n = 4) were 2035 ms and 168 ms, respectively. For treated BM (n = 17) the T1 and T2 values were 2163 ms and 141 ms, respectively. MRF technique appears to be a promising and rapid quantitative method for the characterization of free water content and tumor morphology in BMs.

Accuracy and repeatability of joint sparsity multi-component estimation in MR Fingerprinting

MR fingerprinting (MRF) is a promising method for quantitative characterization of tissues. Often, voxel-wise measurements are made, assuming a single tissue-type per voxel. Alternatively, the Sparsity Promoting Iterative Joint Non-negative least squares Multi-Component MRF method (SPIJN-MRF) facilitates tissue parameter estimation for identified components as well as partial volume segmentations. The aim of this paper was to evaluate the accuracy and repeatability of the SPIJN-MRF parameter estimations and partial volume segmentations. This was done (1) through numerical simulations based on the BrainWeb phantoms and (2) using in vivo acquired MRF data from 5 subjects that were scanned on the same week-day for 8 consecutive weeks. The partial volume segmentations of the SPIJN-MRF method were compared to those obtained by two conventional methods: SPM12 and FSL. SPIJN-MRF showed higher accuracy in simulations in comparison to FSL- and SPM12-based segmentations: Fuzzy Tanimoto Coefficients (FTC) comparing these segmentations and Brainweb references were higher than 0.95 for SPIJN-MRF in all the tissues and between 0.6 and 0.7 for SPM12 and FSL in white and gray matter and between 0.5 and 0.6 in CSF. For the in vivo MRF data, the estimated relaxation times were in line with literature and minimal variation was observed. Furthermore, the coefficient of variation (CoV) for estimated tissue volumes with SPIJN-MRF were 10.5% for the myelin water, 6.0% for the white matter, 5.6% for the gray matter, 4.6% for the CSF and 1.1% for the total brain volume. CoVs for CSF and total brain volume measured on the scanned data for SPIJN-MRF were in line with those obtained with SPM12 and FSL. The CoVs for white and gray matter volumes were distinctively higher for SPIJN-MRF than those measured with SPM12 and FSL. In conclusion, the use of SPIJN-MRF provides accurate and precise tissue relaxation parameter estimations taking into account intrinsic partial volume effects. It facilitates obtaining tissue fraction maps of prevalent tissues including myelin water which can be relevant for evaluating diseases affecting the white matter.

Synthetic MRI with Magnetic Resonance Spin TomogrAphy in Time-Domain (MR-STAT): Results from a Prospective Cross-Sectional Clinical Trial

BACKGROUND

Magnetic Resonance Spin TomogrAphy in Time-domain (MR-STAT) can reconstruct whole-brain multi-parametric quantitative maps (eg, T₁ , T₂ ) from a 5-minute MR acquisition. These quantitative maps can be leveraged for synthetization of clinical image contrasts.

PURPOSE

The objective was to assess image quality and overall diagnostic accuracy of synthetic MR-STAT contrasts compared to conventional contrast-weighted images.

STUDY TYPE

Prospective cross-sectional clinical trial.

POPULATION

Fifty participants with a median age of 45 years (range: 21-79 years) consisting of 10 healthy participants and 40 patients with neurological diseases (brain tumor, epilepsy, multiple sclerosis or stroke).

FIELD STRENGTH/SEQUENCE

3T/Conventional contrast-weighted imaging (T₁ /T₂ weighted, proton density [PD] weighted, and fluid-attenuated inversion recovery [FLAIR]) and a MR-STAT acquisition (2D Cartesian spoiled gradient echo with varying flip angle preceded by a non-selective inversion pulse).

ASSESSMENT

Quantitative T₁ , T₂ , and PD maps were computed from the MR-STAT acquisition, from which synthetic contrasts were generated. Three neuroradiologists blinded for image type and disease randomly and independently evaluated synthetic and conventional datasets for image quality and diagnostic accuracy, which was assessed by comparison with the clinically confirmed diagnosis.

STATISTICAL TESTS

Image quality and consequent acceptability for diagnostic use was assessed with a McNemar's test (one-sided α = 0.025). Wilcoxon signed rank test with a one-sided α = 0.025 and a margin of Δ = 0.5 on the 5-level Likert scale was used to assess non-inferiority.

RESULTS

All data sets were similar in acceptability for diagnostic use (≥3 Likert-scale) between techniques (T₁ w:P = 0.105, PDw:P = 1.000, FLAIR:P = 0.564). However, only the synthetic MR-STAT T₂ weighted images were significantly non-inferior to their conventional counterpart; all other synthetic datasets were inferior (T₁ w:P = 0.260, PDw:P = 1.000, FLAIR:P = 1.000). Moreover, true positive/negative rates were similar between techniques (conventional: 88%, MR-STAT: 84%).

DATA CONCLUSION

MR-STAT is a quantitative technique that may provide radiologists with clinically useful synthetic contrast images within substantially reduced scan time.

EVIDENCE LEVEL

1 Technical Efficacy: Stage 2.

Synthetic MRI improves radiomics-based glioblastoma survival prediction

Glioblastoma is an aggressive and fast-growing brain tumor with poor prognosis. Predicting the expected survival of patients with glioblastoma is a key task for efficient treatment and surgery planning. Survival predictions could be enhanced by means of a radiomic system. However, these systems demand high numbers of multicontrast images, the acquisitions of which are time consuming, giving rise to patient discomfort and low healthcare system efficiency. Synthetic MRI could favor deployment of radiomic systems in the clinic by allowing practitioners not only to reduce acquisition time, but also to retrospectively complete databases or to replace artifacted images. In this work we analyze the replacement of an actually acquired MR weighted image by a synthesized version to predict survival of glioblastoma patients with a radiomic system. Each synthesized version was realistically generated from two acquired images with a deep learning synthetic MRI approach based on a convolutional neural network. Specifically, two weighted images were considered for the replacement one at a time, a T2w and a FLAIR, which were synthesized from the pairs T1w and FLAIR, and T1w and T2w, respectively. Furthermore, a radiomic system for survival prediction, which can classify patients into two groups (survival >480 days and ≤ 480 days), was built. Results show that the radiomic system fed with the synthesized image achieves similar performance compared with using the acquired one, and better performance than a model that does not include this image. Hence, our results confirm that synthetic MRI does add to glioblastoma survival prediction within a radiomics-based approach.

Investigated diagnostic value of synthetic relaxometry, three-dimensional pseudo-continuous arterial spin labelling and diffusion-weighted imaging in the grading of glioma

BACKGROUND

To investigate the performance of synthetic relaxometry, three-dimensional pseudo-continuous arterial spin labelling (pCASL) and diffusion-weighted imaging (DWI) in differentiating high-grade gliomas (HGGs) from low-grade gliomas (LGGs) and to compare with the conventional MRI.

METHODS

Seventy-two patients with gliomas (including 27 LGGs and 45 HGGs) were studied using synthetic magnetic resonance imaging (sy-MRI), pCASL, and DWI with a 3.0 T MR scanner. T1 relaxometry (T1), T2 relaxometry (T2), as well as proton density (PD) from sy-MRI, cerebral blood flow (CBF) from pCASL, apparent diffusion coefficient (ADC) from DWI and enhancement quality (EQ), proportion enhancing (PE) from conventional contrast enhanced image based Visually-Accessible-Rembrandt-Images (VASARI) scoring system, were all analyzed by two radiologists. The Student's t-test, Mann-Whitney U test or Fisher's exact test was used to compare the parameters between LGGs and HGGs. The diagnostic performance of each parameter and their combination for glioma grading were analyzed.

RESULTS

Significant statistical differences in T1, PD, CBF, ADC, EQ and PE are observed between LGGs and HGGs (all P < 0.001). The ADC values have higher discrimination abilities compared with other univariable parameters, with the AUC of 0.905. AUC values for conventional contrast-enhanced method, EQ and PE from VASARI, and conventional contrast-free method, CBF + ADC, are 0.873 and 0.912 respectively. The combined T1, PD, CBF and ADC model had the best performance for differentiating LGGs and HGGs with AUC, sensitivity and specificity of 0.993, 95.5%, 100%, respectively.

CONCLUSIONS

Relaxometry parameters derived from synthetic MRI contributed to the discrimination of low-grade gliomas from high-grade gliomas. Proposed contrast-free approach combining T1, PD, CBF and ADC showed a strong discriminative power, and outperformed conventional approaches.