Efficient pulse sequence design framework for high-dimensional MR fingerprinting scans using systematic error index

PURPOSE

For effective optimization of MR fingerprinting (MRF) pulse sequences, estimating and minimizing errors from actual scan conditions are crucial. Although virtual-scan simulations offer an approximation to these errors, their computational demands become expensive for high-dimensional MRF frameworks, where interactions between more than two tissue properties are considered. This complexity makes sequence optimization impractical. We introduce a new mathematical model, the systematic error index (SEI), to address the scalability challenges for high-dimensional MRF sequence design.

METHODS

By eliminating the need to perform dictionary matching, the SEI model approximates quantification errors with low computational costs. The SEI model was validated in comparison with virtual-scan simulations. The SEI model was further applied to optimize three high-dimensional MRF sequences that quantify two to four tissue properties. The optimized scans were examined in simulations and healthy subjects.

RESULTS

The proposed SEI model closely approximated the virtual-scan simulation outcomes while achieving hundred- to thousand-times acceleration in the computational speed. In both simulation and in vivo experiments, the optimized MRF sequences yield higher measurement accuracy with fewer undersampling artifacts at shorter scan times than the heuristically designed sequences.

CONCLUSION

We developed an efficient method for estimating real-world errors in MRF scans with high computational efficiency. Our results illustrate that the SEI model could approximate errors both qualitatively and quantitatively. We also proved the practicality of the SEI model of optimizing sequences for high-dimensional MRF frameworks with manageable computational power. The optimized high-dimensional MRF scans exhibited enhanced robustness against undersampling and system imperfections with faster scan times.

Beyond the Conventional Structural MRI: Clinical Application of Deep Learning Image Reconstruction and Synthetic MRI of the Brain

ABSTRACT

Recent technological advancements have revolutionized routine brain magnetic resonance imaging (MRI) sequences, offering enhanced diagnostic capabilities in intracranial disease evaluation. This review explores 2 pivotal breakthrough areas: deep learning reconstruction (DLR) and quantitative MRI techniques beyond conventional structural imaging. DLR using deep neural networks facilitates accelerated imaging with improved signal-to-noise ratio and spatial resolution, enhancing image quality with short scan times. DLR focuses on supervised learning applied to clinical implementation and applications. Quantitative MRI techniques, exemplified by 2D multidynamic multiecho, 3D quantification using interleaved Look-Locker acquisition sequences with T2 preparation pulses, and magnetic resonance fingerprinting, enable precise calculation of brain-tissue parameters and further advance diagnostic accuracy and efficiency. Potential DLR instabilities and quantification and bias limitations will be discussed. This review underscores the synergistic potential of DLR and quantitative MRI, offering prospects for improved brain imaging beyond conventional methods.

Multiparametric Characterization of Focal Cortical Dysplasia Using 3D MR Fingerprinting

OBJECTIVES

To develop a multiparametric machine-learning (ML) framework using high-resolution 3 dimensional (3D) magnetic resonance (MR) fingerprinting (MRF) data for quantitative characterization of focal cortical dysplasia (FCD).

MATERIALS

We included 119 subjects, 33 patients with focal epilepsy and histopathologically confirmed FCD, 60 age- and gender-matched healthy controls (HCs), and 26 disease controls (DCs). Subjects underwent whole-brain 3 Tesla MRF acquisition, the reconstruction of which generated T1 and T2 relaxometry maps. A 3D region of interest was manually created for each lesion, and z-score normalization using HC data was performed. We conducted 2D classification with ensemble models using MRF T1 and T2 mean and standard deviation from gray matter and white matter for FCD versus controls. Subtype classification additionally incorporated entropy and uniformity of MRF metrics, as well as morphometric features from the morphometric analysis program (MAP). We translated 2D results to individual probabilities using the percentage of slices above an adaptive threshold. These probabilities and clinical variables were input into a support vector machine for individual-level classification. Fivefold cross-validation was performed and performance metrics were reported using receiver-operating-characteristic-curve analyses.

RESULTS

FCD versus HC classification yielded mean sensitivity, specificity, and accuracy of 0.945, 0.980, and 0.962, respectively; FCD versus DC classification achieved 0.918, 0.965, and 0.939. In comparison, visual review of the clinical magnetic resonance imaging (MRI) detected 48% (16/33) of the lesions by official radiology report. In the subgroup where both clinical MRI and MAP were negative, the MRF-ML models correctly distinguished FCD patients from HCs and DCs in 98.3% of cross-validation trials for the magnetic resonance imaging negative group and MAP negative group. Type II versus non-type-II classification exhibited mean sensitivity, specificity, and accuracy of 0.835, 0.823, and 0.83, respectively; type IIa versus IIb classification showed 0.85, 0.9, and 0.87. In comparison, the transmantle sign was present in 58% (7/12) of the IIb cases.

INTERPRETATION

The MRF-ML framework presented in this study demonstrated strong efficacy in noninvasively classifying FCD from normal cortex and distinguishing FCD subtypes. ANN NEUROL 2024.

Brain tumor-associated epilepsies in adulthood: Current state of diagnostic and individual treatment options

Brain tumors are one of the most frequent causes of structural epilepsy and set a major burden on treatment costs and the social integrity of patients. Although promising oncological treatment strategies are already available, epileptological treatment is often intractable and requires lifelong epileptological care. Therefore, treatment strategies must be adapted to age-related needs, and specific aspects of late-onset epilepsy (LOE) must be considered. The practical implementation of individual decisions from tumor boards and the current state of the art in scientific knowledge about pathological mechanisms, modern diagnostic procedures and biomarkers, and patient-individualized treatment options into practical epileptological disease management is a prerequisite. This narrative review focuses on the current work progress regarding pathogenesis, diagnosis, and therapy. Exemplarily, interdisciplinary approaches for optimized individualized therapy will be discussed, emphasizing the combination of neurological-epileptological and oncological perspectives.

Multiparametric quantification of T1 and T2 relaxation time of bone metastasis in comparison with red or fatty bone marrow using magnetic resonance fingerprinting

OBJECTIVES

To assess the T1 and T2 values of bone marrow lesions in spine and pelvis derived from magnetic resonance fingerprinting (MRF) and to evaluate the differences in values among bone metastasis, red marrow and fatty marrow.

METHODS

Sixty patients who underwent lumbar spine and pelvic MRI with magnetic resonance fingerprinting were retrospectively included. Among eligible patients, those with bone metastasis, benign red marrow deposition and normal fatty marrow were identified. Two radiologists independently measured the T1 and T2 values from metastatic bone lesions, fatty marrow, and red marrow deposition on three-dimensional-magnetic resonance fingerprinting. Intergroup comparison and interobserver agreement were analyzed.

RESULTS

T1 relaxation time was significantly higher in osteoblastic metastasis than in red marrow (1674.6 ± 436.3 vs 858.7 ± 319.5, p < .001). Intraclass correlation coefficients for T1 and T2 values were 0.96 (p < 0.001) and 0.83 (p < 0.001), respectively. T2 relaxation time of osteoblastic metastasis and red marrow deposition had no evidence of a difference (osteoblastic metastasis, 57.9 ± 25.0 vs red marrow, 58.0 ± 34.4, p = 0.45), as were the average T2 values of osteolytic metastasis and red marrow deposition (osteolytic metastasis, 45.3 ± 15.1 vs red marrow, 58.0 ± 34.4, p = 0.63).

CONCLUSIONS

We report the feasibility of three-dimensional-magnetic resonance fingerprinting based quantification of bone marrow to differentiate bone metastasis from red marrow. Simultaneous T1 and T2 quantification of metastasis and red marrow deposition was possible in spine and pelvis and showed significant different values with excellent inter-reader agreement.

ADVANCE IN KNOWLEDGE

T1 values from three-dimensional-magnetic resonance fingerprinting might be a useful quantifier for evaluating bone marrow lesions.

Combining magnetic resonance fingerprinting with voxel-based morphometric analysis to reduce false positives for focal cortical dysplasia detection

OBJECTIVE

We aim to improve focal cortical dysplasia (FCD) detection by combining high-resolution, three-dimensional (3D) magnetic resonance fingerprinting (MRF) with voxel-based morphometric magnetic resonance imaging (MRI) analysis.

METHODS

We included 37 patients with pharmacoresistant focal epilepsy and FCD (10 IIa, 15 IIb, 10 mild Malformation of Cortical Development [mMCD], and 2 mMCD with oligodendroglial hyperplasia and epilepsy [MOGHE]). Fifty-nine healthy controls (HCs) were also included. 3D lesion labels were manually created. Whole-brain MRF scans were obtained with 1 mm3 isotropic resolution, from which quantitative T1 and T2 maps were reconstructed. Voxel-based MRI postprocessing, implemented with the morphometric analysis program (MAP18), was performed for FCD detection using clinical T1w images, outputting clusters with voxel-wise lesion probabilities. Average MRF T1 and T2 were calculated in each cluster from MAP18 output for gray matter (GM) and white matter (WM) separately. Normalized MRF T1 and T2 were calculated by z-scores using HCs. Clusters that overlapped with the lesion labels were considered true positives (TPs); clusters with no overlap were considered false positives (FPs). Two-sample t-tests were performed to compare MRF measures between TP/FP clusters. A neural network model was trained using MRF values and cluster volume to distinguish TP/FP clusters. Ten-fold cross-validation was used to evaluate model performance at the cluster level. Leave-one-patient-out cross-validation was used to evaluate performance at the patient level.

RESULTS

MRF metrics were significantly higher in TP than FP clusters, including GM T1, normalized WM T1, and normalized WM T2. The neural network model with normalized MRF measures and cluster volume as input achieved mean area under the curve (AUC) of .83, sensitivity of 82.1%, and specificity of 71.7%. This model showed superior performance over direct thresholding of MAP18 FCD probability map at both the cluster and patient levels, eliminating ≥75% FP clusters in 30% of patients and ≥50% of FP clusters in 91% of patients.

SIGNIFICANCE

This pilot study suggests the efficacy of MRF for reducing FPs in FCD detection, due to its quantitative values reflecting in vivo pathological changes. © 2024 International League Against Epilepsy.

Motion Robust MR Fingerprinting Scan to Image Neonates With Prenatal Opioid Exposure

PURPOSE

To explore whether MR fingerprinting (MRF) scans provide motion-robust and quantitative brain tissue measurements for non-sedated infants with prenatal opioid exposure (POE).

STUDY TYPE

Prospective.

POPULATION

13 infants with POE (3 male; 12 newborns (age 7-65 days) and 1 infant aged 9-months).

FIELD STRENGTH/SEQUENCE

3T, 3D T1-weighted MPRAGE, 3D T2-weighted TSE and MRF sequences.

ASSESSMENT

The image quality of MRF and MRI was assessed in a fully crossed, multiple-reader, multiple-case study. Sixteen image quality features in three types-image artifacts, structure and myelination visualization-were ranked by four neuroradiologists (8, 7, 5, and 8 years of experience respectively), using a 3-point scale. MRF T1 and T2 values in 8 white matter brain regions were compared between babies younger than 1 month and babies between 1 and 2 months.

STATISTICAL TESTS

Generalized estimating equations model to test the significance of differences of regional T1 and T2 values of babies under 1 month and those older. MRI and MRF image quality was assessed using Gwet's second order auto-correlation coefficient (AC2) with confidence levels. The Cochran-Mantel-Haenszel test was used to assess the difference in proportions between MRF and MRI for all features and stratified by the type of features. A P value <0.05 was considered statistically significant.

RESULTS

The MRF of two infants were excluded in T1 and T2 value analysis due to severe motion artifact but were included in the image quality assessment. In infants under 1 month of age (N = 6), the T1 and T2 values were significantly higher compared to those between 1 and 2 months of age (N = 4). MRF images showed significantly higher image quality ratings in all three feature types compared to MRI images.

CONCLUSIONS

MR Fingerprinting scans have potential to be a motion-robust and efficient method for nonsedated infants.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 1.

DTI-MR fingerprinting for rapid high-resolution whole-brain T1 , T2 , proton density, ADC, and fractional anisotropy mapping

PURPOSE

This study aims to develop a high-efficiency and high-resolution 3D imaging approach for simultaneous mapping of multiple key tissue parameters for routine brain imaging, including T1 , T2 , proton density (PD), ADC, and fractional anisotropy (FA). The proposed method is intended for pushing routine clinical brain imaging from weighted imaging to quantitative imaging and can also be particularly useful for diffusion-relaxometry studies, which typically suffer from lengthy acquisition time.

METHODS

To address challenges associated with diffusion weighting, such as shot-to-shot phase variation and low SNR, we integrated several innovative data acquisition and reconstruction techniques. Specifically, we used M1-compensated diffusion gradients, cardiac gating, and navigators to mitigate phase variations caused by cardiac motion. We also introduced a data-driven pre-pulse gradient to cancel out eddy currents induced by diffusion gradients. Additionally, to enhance image quality within a limited acquisition time, we proposed a data-sharing joint reconstruction approach coupled with a corresponding sequence design.

RESULTS

The phantom and in vivo studies indicated that the T1 and T2 values measured by the proposed method are consistent with a conventional MR fingerprinting sequence and the diffusion results (including diffusivity, ADC, and FA) are consistent with the spin-echo EPI DWI sequence.

CONCLUSION

The proposed method can achieve whole-brain T1 , T2 , diffusivity, ADC, and FA maps at 1-mm isotropic resolution within 10 min, providing a powerful tool for investigating the microstructural properties of brain tissue, with potential applications in clinical and research settings.

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.

Time-efficient, high-resolution 3T whole-brain relaxometry using 3D-QALAS with wave-CAIPI readouts

PURPOSE

Volumetric, high-resolution, quantitative mapping of brain-tissue relaxation properties is hindered by long acquisition times and SNR challenges. This study combines time-efficient wave-controlled aliasing in parallel imaging (wave-CAIPI) readouts with the 3D quantification using an interleaved Look-Locker acquisition sequence with a T2 preparation pulse (3D-QALAS), enabling full-brain quantitative T1 , T2 , and proton density (PD) maps at 1.15-mm3 isotropic voxels in 3 min.

METHODS

Wave-CAIPI readouts were embedded in the standard 3D-QALAS encoding scheme, enabling full-brain quantitative parameter maps (T1 , T2 , and PD) at acceleration factors of R = 3 × 2 with minimum SNR loss due to g-factor penalties. The quantitative parameter maps were estimated using a dictionary-based mapping algorithm incorporating inversion efficiency and B1 -field inhomogeneity effects. The parameter maps using the accelerated protocol were quantitatively compared with those obtained from the conventional 3D-QALAS sequence using GRAPPA acceleration of R = 2 in the ISMRM/NIST phantom, and in 10 healthy volunteers.

RESULTS

When tested in both the ISMRM/NIST phantom and 10 healthy volunteers, the quantitative maps using the accelerated protocol showed excellent agreement against those obtained from conventional 3D-QALAS at RGRAPPA  = 2.

CONCLUSION

Three-dimensional QALAS enhanced with wave-CAIPI readouts enables time-efficient, full-brain quantitative T1 , T2 , and PD mapping at 1.15 mm3 in 3 min at R = 3 × 2 acceleration. The quantitative maps obtained from the accelerated wave-CAIPI 3D-QALAS protocol showed very similar values to those from the standard 3D-QALAS (R = 2) protocol, alluding to the robustness and reliability of the proposed method.

Magnetic resonance coherence pathway unraveling

Efficiently acquiring multi-contrast magnetic resonance imaging data is crucial for patient comfort and clinical throughput. Developing scan acceleration methods tailored for specific applications drastically improves the value of an MRI examination. Here, we propose a novel method to control the aliasing of simultaneously acquired images of multiple spin echo coherence pathways with the goal of producing high quality multi-contrast images from a single acquisition. Modulating the radiofrequency phase of several pulses applied in brief succession also uniquely modulates the phase of spin echo coherence pathways. A method, termed magnetic resonance coherence pathway unraveling (MR-CPU), to control the aliasing of simultaneously acquired coherence pathway images is developed here along with parallel imaging-based reconstruction methods to separate them. MR-CPU was validated in phantom experiments and tested in vivo. High levels of correlation between reference pathway images and MR-CPU-derived coherence pathway images were found from the phantom experiments. Minimal artifacts arising from the separation of the overlapped coherence pathway images were observed in vivo. MR-CPU provides a novel mechanism through which to acquire and separate multiple overlapped coherence pathway images, thus adding to the diagnostic potential of an MRI exam without the penalty of additional scan time.

Fast deep learning reconstruction techniques for preclinical magnetic resonance fingerprinting

We propose a deep learning (DL) model and a hyperparameter optimization strategy to reconstruct T1 and T2 maps acquired with the magnetic resonance fingerprinting (MRF) methodology. We applied two different MRF sequence routines to acquire images of ex vivo rat brain phantoms using a 7-T preclinical scanner. Subsequently, the DL model was trained using experimental data, completely excluding the use of any theoretical MRI signal simulator. The best combination of the DL parameters was implemented by an automatic hyperparameter optimization strategy, whose key aspect is to include all the parameters to the fit, allowing the simultaneous optimization of the neural network architecture, the structure of the DL model, and the supervised learning algorithm. By comparing the reconstruction performances of the DL technique with those achieved from the traditional dictionary-based method on an independent dataset, the DL approach was shown to reduce the mean percentage relative error by a factor of 3 for T1 and by a factor of 2 for T2 , and to improve the computational time by at least a factor of 37. Furthermore, the proposed DL method enables maintaining comparable reconstruction performance, even with a lower number of MRF images and a reduced k-space sampling percentage, with respect to the dictionary-based method. Our results suggest that the proposed DL methodology may offer an improvement in reconstruction accuracy, as well as speeding up MRF for preclinical, and in prospective clinical, investigations.

Improving motion robustness of 3D MR fingerprinting with a fat navigator

PURPOSE

To develop a 3D MR fingerprinting (MRF) method in combination with fat navigators to improve its motion robustness for neuroimaging.

METHODS

A rapid fat navigator was developed using the stack-of-spirals acquisition and non-Cartesian spiral GRAPPA. The fat navigator module was implemented in the 3D MRF sequence with high scan efficiency. The developed method was first validated in phantoms and five healthy subjects with intentional head motion. The method was further applied to infants with neonatal opioid withdrawal symptoms. The 3D MRF scans with fat navigators acquired with and without acceleration along the partition-encoding direction were both examined in the study.

RESULTS

Both phantom and in vivo results demonstrated that the added fat navigator modules did not influence the quantification accuracy in MRF. In combination with non-Cartesian spiral GRAPPA, a rapid fat navigator sampling with whole-brain coverage was achieved in ˜0.5 s at 3T, reducing its sensitivity to potential motion. Based on the motion waveforms extracted from fat navigators, the motion robustness of the 3D MRF was largely improved. With the proposed method, the motion-corrupted MRF datasets yielded T1 and T2 maps with significantly reduced artifacts and high correlations with measurements from the reference motion-free MRF scans.

CONCLUSION

We developed a 3D MRF method coupled with rapid fat navigators to improve its motion robustness for quantitative neuroimaging. Our results demonstrate that (1) accurate tissue quantification was preserved with the fat navigator modules and (2) the motion robustness for quantitative tissue mapping was largely improved with the developed method.

SSL-QALAS: Self-Supervised Learning for rapid multiparameter estimation in quantitative MRI using 3D-QALAS

PURPOSE

To develop and evaluate a method for rapid estimation of multiparametric T1 , T2 , proton density, and inversion efficiency maps from 3D-quantification using an interleaved Look-Locker acquisition sequence with T2 preparation pulse (3D-QALAS) measurements using self-supervised learning (SSL) without the need for an external dictionary.

METHODS

An SSL-based QALAS mapping method (SSL-QALAS) was developed for rapid and dictionary-free estimation of multiparametric maps from 3D-QALAS measurements. The accuracy of the reconstructed quantitative maps using dictionary matching and SSL-QALAS was evaluated by comparing the estimated T1 and T2 values with those obtained from the reference methods on an International Society for Magnetic Resonance in Medicine/National Institute of Standards and Technology phantom. The SSL-QALAS and the dictionary-matching methods were also compared in vivo, and generalizability was evaluated by comparing the scan-specific, pre-trained, and transfer learning models.

RESULTS

Phantom experiments showed that both the dictionary-matching and SSL-QALAS methods produced T1 and T2 estimates that had a strong linear agreement with the reference values in the International Society for Magnetic Resonance in Medicine/National Institute of Standards and Technology phantom. Further, SSL-QALAS showed similar performance with dictionary matching in reconstructing the T1 , T2 , proton density, and inversion efficiency maps on in vivo data. Rapid reconstruction of multiparametric maps was enabled by inferring the data using a pre-trained SSL-QALAS model within 10 s. Fast scan-specific tuning was also demonstrated by fine-tuning the pre-trained model with the target subject's data within 15 min.

CONCLUSION

The proposed SSL-QALAS method enabled rapid reconstruction of multiparametric maps from 3D-QALAS measurements without an external dictionary or labeled ground-truth training data.

Multi-parametric MRI for radiotherapy simulation

Magnetic resonance imaging (MRI) has become an important imaging modality in the field of radiotherapy (RT) in the past decade, especially with the development of various novel MRI and image-guidance techniques. In this review article, we will describe recent developments and discuss the applications of multi-parametric MRI (mpMRI) in RT simulation. In this review, mpMRI refers to a general and loose definition which includes various multi-contrast MRI techniques. Specifically, we will focus on the implementation, challenges, and future directions of mpMRI techniques for RT simulation.

Magnetic Resonance Fingerprinting: A Review of Clinical Applications

ABSTRACT

Magnetic resonance fingerprinting (MRF) is an approach to quantitative magnetic resonance imaging that allows for efficient simultaneous measurements of multiple tissue properties, which are then used to create accurate and reproducible quantitative maps of these properties. As the technique has gained popularity, the extent of preclinical and clinical applications has vastly increased. The goal of this review is to provide an overview of currently investigated preclinical and clinical applications of MRF, as well as future directions. Topics covered include MRF in neuroimaging, neurovascular, prostate, liver, kidney, breast, abdominal quantitative imaging, cardiac, and musculoskeletal applications.

T1 relaxation time of ISMRM/NIST T1 phantom spheres at 7 T

T1 relaxation times of the 14 T1 phantom spheres that make up the standard International Society for Magnetic Resonance in Medicine (ISMRM)/National Institute of Standards and Technology (NIST) system phantom are reported at 7 T. T1 values of six of the 14 T1 spheres at 7 T (with T1 > 270 ms) have been reported previously, but, to the best of our knowledge, not all of the T1s of the 14 T1 spheres at 7 T have been reported before. Given the increasing number of 7-T MRI systems in clinical settings and the increasing need for T1 phantoms that cover a wide range of T1 relaxation times to evaluate rapid T1 mapping techniques at 7 T, it is of high interest to obtain accurate T1 values for all the ISMRM/NIST T1 spheres at 7 T. In this work, T1 relaxation time was measured on a 7-T MRI scanner using an inversion-recovery spin-echo pulse sequence and derived by curve fitting to a signal equation that exhibits insensitivity to B 1 + inhomogeneity. Day-to-day reproducibility was within 0.4% and differences between two different RF coils within 1.5%. T1s of a subset of the 14 spheres were also measured by NMR at 7 T for comparison, and the T1 results were consistent between the MRI and NMR measurements. T1 measurements performed at 3 T on the same 14 spheres using the same sequence and fitting method yielded good agreement (mean percentage difference of -0.4%) with the reference T1 values available from the NIST, reflecting the accuracy of the reported technique despite being without the standard phantom housing. We found that the T1 values of all 14 NiCl₂ spheres are consistently lower at 7 T than at 3 T. Although our results were well reproduced, this study represents initial work to quantify the 7-T T1 values of all 14 NIST T1 spheres outside of the standard housing and does not warrant reproducibility of the ISMRM/NIST system phantom as a whole. A future study to assess the T1 values of a version of the ISMRM/NIST system phantom that fits inside typical commercial coils at 7 T will be very helpful. Nonetheless, the details on our acquisition and curve-fitting methods reported here allow the T1 measurements to be reproduced elsewhere. The T1 values of all 14 spheres reported here will be valuable for the development of quantitative MR fingerprinting and rapid T1 mapping for a large variety of research projects, not only in neuroimaging but also in body MRI, musculoskeletal MRI, and gadolinium contrast-enhanced MRI, each of which is concerned with much shortened T1.

Normative quantitative relaxation atlases for characterization of cortical regions using magnetic resonance fingerprinting

Quantitative magnetic resonance (MR) has been used to study cyto- and myelo-architecture of the human brain non-invasively. However, analyzing brain cortex using high-resolution quantitative MR acquisition can be challenging to perform using 3T clinical scanners. MR fingerprinting (MRF) is a highly efficient and clinically feasible quantitative MR technique that simultaneously provides T1 and T2 relaxation maps. Using 3D MRF from 40 healthy subjects (mean age = 25.6 ± 4.3 years) scanned on 3T magnetic resonance imaging, we generated whole-brain gyral-based normative MR relaxation atlases and investigated cortical-region-based T1 and T2 variations. Gender and age dependency of T1 and T2 variations were additionally analyzed. The coefficient of variation of T1 and T2 for each cortical-region was 3.5% and 7.3%, respectively, supporting low variability of MRF measurements across subjects. Significant differences in T1 and T2 were identified among 34 brain regions (P < 0.001), lower in the precentral, postcentral, paracentral lobule, transverse temporal, lateral occipital, and cingulate areas, which contain sensorimotor, auditory, visual, and limbic functions. Significant correlations were identified between age and T1 and T2 values. This study established whole-brain MRF T1 and T2 atlases of healthy subjects using a clinical 3T scanner, which can provide a quantitative and region-specific baseline for future brain studies and pathology detection.

Novel 3D magnetic resonance fingerprinting radiomics in adult brain tumors: a feasibility study

OBJECTIVES

To test the feasibility of using 3D MRF maps with radiomics analysis and machine learning in the characterization of adult brain intra-axial neoplasms.

METHODS

3D MRF acquisition was performed on 78 patients with newly diagnosed brain tumors including 33 glioblastomas (grade IV), 6 grade III gliomas, 12 grade II gliomas, and 27 patients with brain metastases. Regions of enhancing tumor, non-enhancing tumor, and peritumoral edema were segmented and radiomics analysis with gray-level co-occurrence matrices and gray-level run-length matrices was performed. Statistical analysis was performed to identify features capable of differentiating tumors based on type, grade, and isocitrate dehydrogenase (IDH1) status. Receiver operating curve analysis was performed and the area under the curve (AUC) was calculated for tumor classification and grading. For gliomas, Kaplan-Meier analysis for overall survival was performed using MRF T1 features from enhancing tumor region.

RESULTS

Multiple MRF T1 and T2 features from enhancing tumor region were capable of differentiating glioblastomas from brain metastases. Although no differences were identified between grade 2 and grade 3 gliomas, differentiation between grade 2 and grade 4 gliomas as well as between grade 3 and grade 4 gliomas was achieved. MRF radiomics features were also able to differentiate IDH1 mutant from the wild-type gliomas. Radiomics T1 features for enhancing tumor region in gliomas correlated to overall survival (p < 0.05).

CONCLUSION

Radiomics analysis of 3D MRF maps allows differentiating glioblastomas from metastases and is capable of differentiating glioblastomas from metastases and characterizing gliomas based on grade, IDH1 status, and survival.

KEY POINTS

• 3D MRF data analysis using radiomics offers novel tissue characterization of brain tumors. • 3D MRF with radiomics offers glioma characterization based on grade, IDH1 status, and overall patient survival.

Evaluating whole-brain tissue-property changes in MRI-negative pharmacoresistant focal epilepsies using MR fingerprinting

OBJECTIVE

We aim to quantify whole-brain tissue-property changes in patients with magnetic resonance imaging (MRI)-negative pharmacoresistant focal epilepsy by three-dimensional (3D) magnetic resonance fingerprinting (MRF).

METHODS

We included 30 patients with pharmacoresistant focal epilepsy and negative MRI by official radiology report, as well as 40 age- and gender-matched healthy controls (HCs). MRF scans were obtained with 1 mm³ isotropic resolution. Quantitative T1 and T2 relaxometry maps were reconstructed from MRF and registered to the Montreal Neurological Institute (MNI) space. A two-sample t test was performed in Functional Magnetic Resonance Imaging of the Brain (FMRIB) Software Library (FSL) to evaluate significant abnormalities in patients comparing to HCs, with correction by the threshold-free cluster enhancement (TFCE) method. Subgroups analyses were performed for extra-temporal epilepsy/temporal epilepsy (ETLE/TLE), and for those with/without subtle abnormalities detected by morphometric analysis program (MAP), to investigate each subgroup's pattern of MRF changes. Correlation analyses were performed between the mean MRF values in each significant cluster and seizure-related clinical variables.

RESULTS

Compared to HCs, patients exhibited significant group-level T1 increase ipsilateral to the epileptic origin, in the mesial temporal gray matter (GM) and white matter (WM), temporal pole GM, orbitofrontal GM, hippocampus, and amygdala, with scattered clusters in the neocortical temporal and insular GM. No significant T2 changes were detected. The ETLE subgroup showed a T1-increase pattern similar to the overall cohort, with additional involvement of the ipsilateral anterior cingulate GM. The subgroup of MAP+ patients also showed a T1-increase pattern similar to the overall cohort, with additional cluster in the ipsilateral lateral orbitofrontal GM. Higher T1 was associated with younger seizure-onset age, longer epilepsy duration, and higher seizure frequency.

SIGNIFICANCE

MRF revealed group-level T1 increase in limbic/paralimbic structures ipsilateral to the epileptic origin, in patients with pharmacoresistant focal epilepsy and no apparent lesions on MRI, suggesting that these regions may be commonly affected by seizures in the epileptic brain. The significant association between T1 increase and higher seizure burden may reflect progressive tissue damage.

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.

Motion-corrected 3D-EPTI with efficient 4D navigator acquisition for fast and robust whole-brain quantitative imaging

PURPOSE

To develop a motion estimation and correction method for motion-robust three-dimensional (3D) quantitative imaging with 3D-echo-planar time-resolved imaging.

THEORY AND METHODS

The 3D-echo-planar time-resolved imaging technique was designed with additional four-dimensional navigator acquisition (x-y-z-echoes) to achieve fast and motion-robust quantitative imaging of the human brain. The four-dimensional-navigator is inserted into the relaxation-recovery deadtime of the sequence in every pulse TR (∼2 s) to avoid extra scan time, and to provide continuous tracking of the 3D head motion and B0 -inhomogeneity changes. By using an optimized spatiotemporal encoding combined with a partial-Fourier scheme, the navigator acquires a large central k-t data block for accurate motion estimation using only four small-flip-angle excitations and readouts, resulting in negligible signal-recovery reduction to the 3D-echo-planar time-resolved imaging acquisition. By incorporating the estimated motion and B0 -inhomogeneity changes into the reconstruction, multi-contrast images can be recovered with reduced motion artifacts.

RESULTS

Simulation shows the cost to the SNR efficiency from the added navigator acquisitions is <1%. Both simulation and in vivo retrospective experiments were conducted, that demonstrate the four-dimensional navigator provided accurate estimation of the 3D motion and B0 -inhomogeneity changes, allowing effective reduction of image artifacts in quantitative maps. Finally, in vivo prospective undersampling acquisition was performed with and without head motion, in which the motion corrupted data after correction show close image quality and consistent quantifications to the motion-free scan, providing reliable quantitative measurements even with head motion.

CONCLUSION

The proposed four-dimensional navigator acquisition provides reliable tracking of the head motion and B0 change with negligible SNR cost, equips the 3D-echo-planar time-resolved imaging technique for motion-robust and efficient quantitative imaging.

Black Line Sign in Focal Cortical Dysplasia IIB: A 7T MRI and Electroclinicopathologic Study

BACKGROUND AND OBJECTIVES

We aim to provide detailed imaging-electroclinicopathologic characterization of the black line sign, a novel MRI marker for focal cortical dysplasia (FCD) IIB.

METHODS

7T T2*-weighted gradient-echo (T2*w-GRE) images were retrospectively reviewed in a consecutive cohort of patients with medically intractable epilepsy with pathology-proven FCD II, for the occurrence of the black line sign. We examined the overlap between the black line region and the seizure-onset zone (SOZ) defined by intracranial EEG (ICEEG) and additionally assessed whether complete inclusion of the black line region in the surgical resection was associated with postoperative seizure freedom. The histopathologic specimen was aligned with the MRI to investigate the pathologic underpinning of the black line sign. Region-of-interest-based quantitative MRI (qMRI) analysis on the 7T T1 map was performed in the black line region, entire lesional gray matter (GM), and contralateral/ipsilateral normal gray and white matter (WM).

RESULTS

We included 20 patients with FCD II (14 IIB and 6 IIA). The black line sign was identified in 12/14 (85.7%) of FCD IIB and 0/6 of FCD IIA on 7T T2*w-GRE. The black line region was highly concordant with the ICEEG-defined SOZ (5/7 complete and 2/7 partial overlap). Seizure freedom was seen in 8/8 patients whose black line region was completely included in the surgical resection; in the 2 patients whose resection did not completely include the black line region, both had recurring seizures. Inclusion of the black line region in the surgical resection was significantly associated with seizure freedom (p = 0.02). QMRI analyses showed that the T1 mean value of the black line region was significantly different from the WM (p < 0.001), but similar to the GM. Well-matched histopathologic slices in one case revealed accumulated dysmorphic neurons and balloon cells in the black line region.

DISCUSSION

The black line sign may serve as a noninvasive marker for FCD IIB. Both MRI-pathology and qMRI analyses suggest that the black line region was an abnormal GM component within the FCD. Being highly concordant with ICEEG-defined SOZ and significantly associated with seizure freedom when included in resection, the black line sign may contribute to the planning of ICEEG/surgery of patients with medically intractable epilepsy with FCD IIB.

CLASSIFICATION OF EVIDENCE

This study provides Class II evidence that in individuals with intractable focal epilepsy undergoing resection who have a 7T MRI with adequate image quality, the presence of the black line sign may suggest FCD IIB, be concordant with SOZ from ICEEG, and be associated with more seizure freedom if fully included in resection.

Characterizing Thalamic and Basal Ganglia Nuclei in Medically Intractable Focal Epilepsy by MR Fingerprinting

OBJECTIVES

Magnetic resonance fingerprinting (MRF) is a novel, quantitative and noninvasive technique to measure brain tissue properties. We aim to use MRF for characterizing normal-appearing thalamic and basal ganglia nuclei in the epileptic brain.

METHODS

A 3D MRF protocol (1mm³ isotropic resolution) was acquired from 48 patients with unilateral medically refractory focal epilepsy and 39 healthy controls (HCs). Whole-brain T1 and T2 maps (containing T1 and T2 relaxation times) were reconstructed for each subject. Ten subcortical nuclei in the thalamus and basal ganglia were segmented as regions of interest (ROIs), within which the mean T1 and T2 values, as well as their coefficient of variation (CV) were compared between the patients and HCs at group level. Subgroup and correlation analyses were performed to examine the relationship between significant MRF measures and various clinical characteristics. Using significantly abnormal MRF measures from the group-level analyses, support vector machine (SVM) and logistic regression machine learning models were built and tested with 5-fold and 10-fold cross-validations, to separate patients from HCs, and to separate patients with left-sided and right-sided epilepsy, at individual level.

RESULTS

MRF revealed increased T1 mean value in the ipsilateral thalamus and nucleus accumbens; increased T1 CV in bilateral thalamus, bilateral pallidum, and ipsilateral caudate; and increased T2 CV in the ipsilateral thalamus in patients compared to HCs (P<0.05, FDR corrected). The SVM classifier produced 78.2% average accuracy to separate individual patients from HCs, with AUC of 0.83. The logistic regression classifier produced 67.4% average accuracy to separate patients with left-sided and right-sided epilepsy, with AUC of 0.72.

SIGNIFICANCE

MRF revealed bilateral tissue-property changes in the normal-appearing thalamus and basal ganglia, with ipsilateral predominance and thalamic preference, suggesting subcortical involvement/impairment in patients with medically intractable focal epilepsy. The individual-level performance of the MRF-based machine-learning models suggests potential opportunities for predicting lateralization.

Learning-based optimization of acquisition schedule for magnetization transfer contrast MR fingerprinting

Magnetization transfer contrast MR fingerprinting (MTC-MRF) is a novel quantitative imaging method that simultaneously quantifies free bulk water and semisolid macromolecule parameters using pseudo-randomized scan parameters. To improve acquisition efficiency and reconstruction accuracy, the optimization of MRF sequence design has been of recent interest in the MRF field, but has been challenging due to the large number of degrees of freedom to be optimized in the sequence. Herein, we propose a framework for learning-based optimization of the acquisition schedule (LOAS), which optimizes RF saturation-encoded MRF acquisitions with a minimal number of scan parameters for tissue parameter determination. In a supervised learning framework, scan parameters were subsequently updated to minimize a predefined loss function that can directly represent tissue quantification errors. We evaluated the performance of the proposed approach with a numerical phantom and in in vivo experiments. For validation, MRF images were synthesized using the tissue parameters estimated from a fully connected neural network framework and compared with references. Our results showed that LOAS outperformed existing indirect optimization methods with regard to quantification accuracy and acquisition efficiency. The proposed LOAS method could be a powerful optimization tool in the design of MRF pulse sequences.

Blip up-down acquisition for spin- and gradient-echo imaging (BUDA-SAGE) with self-supervised denoising enables efficient T2 , T2 *, para- and dia-magnetic susceptibility mapping

PURPOSE

To rapidly obtain high resolution T₂ , T₂ *, and quantitative susceptibility mapping (QSM) source separation maps with whole-brain coverage and high geometric fidelity.

METHODS

We propose Blip Up-Down Acquisition for Spin And Gradient Echo imaging (BUDA-SAGE), an efficient EPI sequence for quantitative mapping. The acquisition includes multiple T₂ *-, T₂ '-, and T₂ -weighted contrasts. We alternate the phase-encoding polarities across the interleaved shots in this multi-shot navigator-free acquisition. A field map estimated from interim reconstructions was incorporated into the joint multi-shot EPI reconstruction with a structured low rank constraint to eliminate distortion. A self-supervised neural network (NN), MR-Self2Self (MR-S2S), was used to perform denoising to boost SNR. Using Slider encoding allowed us to reach 1 mm isotropic resolution by performing super-resolution reconstruction on volumes acquired with 2 mm slice thickness. Quantitative T₂ (=1/R₂ ) and T₂ * (=1/R₂ *) maps were obtained using Bloch dictionary matching on the reconstructed echoes. QSM was estimated using nonlinear dipole inversion on the gradient echoes. Starting from the estimated R₂ /R₂ * maps, R₂ ' information was derived and used in source separation QSM reconstruction, which provided additional para- and dia-magnetic susceptibility maps.

RESULTS

In vivo results demonstrate the ability of BUDA-SAGE to provide whole-brain, distortion-free, high-resolution, multi-contrast images and quantitative T₂ /T₂ * maps, as well as yielding para- and dia-magnetic susceptibility maps. Estimated quantitative maps showed comparable values to conventional mapping methods in phantom and in vivo measurements.

CONCLUSION

BUDA-SAGE acquisition with self-supervised denoising and Slider encoding enables rapid, distortion-free, whole-brain T₂ /T₂ * mapping at 1 mm isotropic resolution under 90 s.

3D Echo Planar Time-resolved Imaging (3D-EPTI) for ultrafast multi-parametric quantitative MRI

Multi-parametric quantitative MRI has shown great potential to improve the sensitivity and specificity of clinical diagnosis and to enhance our understanding of complex brain processes, but suffers from long scan time especially at high spatial resolution. To address this longstanding challenge, we introduce a novel approach, termed 3D Echo Planar Time-resolved Imaging (3D-EPTI), which significantly increases the acceleration capacity of MRI sampling, and provides high acquisition efficiency for multi-parametric MRI. This is achieved by exploiting the spatiotemporal correlation of MRI data at multiple timescales through new encoding strategies within and between efficient continuous readouts. Specifically, an optimized spatiotemporal CAIPI encoding within the readouts combined with a radial-block sampling strategy across the readouts enables an acceleration rate of 800 fold in the k-t space. A subspace reconstruction was employed to resolve thousands of high-quality multi-contrast images. We have demonstrated the ability of 3D-EPTI to provide robust and repeatable whole-brain simultaneous T1, T2, T2*, PD and B1+ mapping at high isotropic resolution within minutes (e.g., 1-mm isotropic resolution in 3 minutes), and to enable submillimeter multi-parametric imaging to study detailed brain structures.

Rigid real-time prospective motion-corrected three-dimensional multiparametric mapping of the human brain

PURPOSE

To develop a rigid real-time prospective motion-corrected multiparametric mapping technique and to test the performance of quantitative estimates.

METHODS

Motion tracking and correction were performed by integrating single-shot spiral navigators into a multiparametric imaging technique, three-dimensional quantification using an interleaved Look-Locker acquisition sequence with a T2 preparation pulse (3D-QALAS). The spiral navigator was optimized, and quantitative measurements were validated using a standard system phantom. The effect of motion correction on whole-brain T1 and T2 mapping under different types of head motion during the scan was evaluated in 10 healthy volunteers. Finally, six patients with Parkinson's disease, which is known to be associated with a high prevalence of motion artifacts, were scanned to evaluate the effectiveness of our method in the real world.

RESULTS

The phantom study demonstrated that the proposed motion correction method did not introduce quantitative bias. Improved parametric map quality and repeatability were shown in volunteer experiments with both in-plane and through-plane motions, comparable to the no-motion ground truth. In real-life validation in patients, the approach showed improved parametric map quality compared to images obtained without motion correction.

CONCLUSIONS

Real-time prospective motion-corrected multiparametric relaxometry based on 3D-QALAS provided robust and repeatable whole-brain multiparametric mapping.

Rapid MR relaxometry using deep learning: An overview of current techniques and emerging trends

Quantitative mapping of MR tissue parameters such as the spin-lattice relaxation time (T₁ ), the spin-spin relaxation time (T₂ ), and the spin-lattice relaxation in the rotating frame (T_1ρ ), referred to as MR relaxometry in general, has demonstrated improved assessment in a wide range of clinical applications. Compared with conventional contrast-weighted (eg T₁ -, T₂ -, or T_1ρ -weighted) MRI, MR relaxometry provides increased sensitivity to pathologies and delivers important information that can be more specific to tissue composition and microenvironment. The rise of deep learning in the past several years has been revolutionizing many aspects of MRI research, including image reconstruction, image analysis, and disease diagnosis and prognosis. Although deep learning has also shown great potential for MR relaxometry and quantitative MRI in general, this research direction has been much less explored to date. The goal of this paper is to discuss the applications of deep learning for rapid MR relaxometry and to review emerging deep-learning-based techniques that can be applied to improve MR relaxometry in terms of imaging speed, image quality, and quantification robustness. The paper is comprised of an introduction and four more sections. Section 2 describes a summary of the imaging models of quantitative MR relaxometry. In Section 3, we review existing "classical" methods for accelerating MR relaxometry, including state-of-the-art spatiotemporal acceleration techniques, model-based reconstruction methods, and efficient parameter generation approaches. Section 4 then presents how deep learning can be used to improve MR relaxometry and how it is linked to conventional techniques. The final section concludes the review by discussing the promise and existing challenges of deep learning for rapid MR relaxometry and potential solutions to address these challenges.

Using magnetic resonance fingerprinting to characterize periventricular nodular heterotopias in pharmacoresistant epilepsy

OBJECTIVE

We aimed to use a novel magnetic resonance fingerprinting (MRF) technique to examine in vivo tissue property characteristics of periventricular nodular heterotopia (PVNH). These characteristics were further correlated with stereotactic-electroencephalographic (SEEG) ictal onset findings.

METHODS

We included five patients with PVNH who had SEEG-guided surgery and at least 1 year of seizure freedom or substantial seizure reduction. High-resolution MRF scans were acquired at 3 T, generating three-dimensional quantitative T₁ and T₂  maps. We assessed the differences between T₁ and T₂  values from the voxels in the nodules located in the SEEG-defined seizure onset zone (SOZ) and non-SOZ, on -individual and group levels. Receiver operating characteristic analyses were performed to obtain the optimal classification performance. Quantification of SEEG ictal onset signals from the nodules was performed by calculating power spectrum density (PSD). The association between PSD and T₁ /T₂  values was further assessed at different frequency bands.

RESULTS

Individual-level analysis showed T₁ was significantly higher in SOZ voxels than non-SOZ voxels (p < .05), with an average 73% classification accuracy. Group-level analysis also showed higher T₁ was significantly associated with SOZ voxels (p < .001). At the optimal cutoff (normalized T₁ of 1.1), a 76% accuracy for classifying SOZ nodules from non-SOZ nodules was achieved. T₁  values were significantly associated with ictal onset PSD at the ultraslow, θ, β, γ, and ripple bands (p < .05). T₂  values were significantly associated with PSD only at the ultraslow band (p < .05).

SIGNIFICANCE

Quantitative MRF measures, especially T₁ , can provide additional noninvasive information to separate nodules in SOZ and non-SOZ. The T₁ and T₂ tissue property changes carry electrophysiological underpinnings relevant to the epilepsy, as shown by their significant positive associations with power changes during the SEEG seizure onset. The use of MRF as a supplementary noninvasive tool may improve presurgical evaluation for patients with PVNH and pharmacoresistant epilepsy.

Radiomics with 3-dimensional magnetic resonance fingerprinting: influence of dictionary design on repeatability and reproducibility of radiomic features

Objectives

We aimed to investigate the influence of magnetic resonance fingerprinting (MRF) dictionary design on radiomic features using in vivo human brain scans.

Methods

Scan-rescans of three-dimensional MRF and conventional T1-weighted imaging were performed on 21 healthy volunteers (9 males and 12 females; mean age, 41.3 ± 14.6 years; age range, 22–72 years). Five patients with multiple sclerosis (3 males and 2 females; mean age, 41.2 ± 7.3 years; age range, 32–53 years) were also included. MRF data were reconstructed using various dictionaries with different step sizes. First- and second-order radiomic features were extracted from each dataset. Intra-dictionary repeatability and inter-dictionary reproducibility were evaluated using intraclass correlation coefficients (ICCs). Features with ICCs > 0.90 were considered acceptable. Relative changes were calculated to assess inter-dictionary biases.

Results

The overall scan-rescan ICCs of MRF-based radiomics ranged from 0.86 to 0.95, depending on dictionary step size. No significant differences were observed in the overall scan-rescan repeatability of MRF-based radiomic features and conventional T1-weighted imaging (p = 1.00). Intra-dictionary repeatability was insensitive to dictionary step size differences. MRF-based radiomic features varied among dictionaries (overall ICC for inter-dictionary reproducibility, 0.62–0.99), especially when step sizes were large. First-order and gray level co-occurrence matrix features were the most reproducible feature classes among different step size dictionaries. T1 map-derived radiomic features provided higher repeatability and reproducibility among dictionaries than those obtained with T2 maps.

Conclusion

MRF-based radiomic features are highly repeatable in various dictionary step sizes. Caution is warranted when performing MRF-based radiomics using datasets containing maps generated from different dictionaries.

Key Points

• MRF-based radiomic features are highly repeatable in various dictionary step sizes.

• Use of different MRF dictionaries may result in variable radiomic features, even when the same MRF acquisition data are used.

• Caution is needed when performing radiomic analysis using data reconstructed from different dictionaries.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00330-022-08555-3.

Quantification of the Intrinsic T1 and T2 of Heschl's Gyri with MR Fingerprinting

PURPOSE

The human primary auditory cortex is located in the Heschl's gyrus (HG). To assess the intrinsic MR property in the gray matter of the HG (GM-HG) with T1 and T2 values using a commercially available MR fingerprinting (MRF) technique.

METHODS

The subjects were 10 healthy volunteers (with 20 HGs; mean age, 31.5 years old; range, 25-53 years old). Coronal T1 and T2 maps were obtained with commercially available MRF using a 3-Tesla MR system. Two radiologists measured the T1 and T2 values of the GM-HG, the GM in the superior temporal gyrus (GM-STG), and the GM in the middle temporal gyrus (GM-MTG) by drawing a ROI on coronal maps.

RESULTS

For both radiologists, the mean T1 and T2 values of the GM-HG were significantly lower than those in the GM-STG or GM-MTG (P < 0.01). The interobserver reliability using the intraclass correlation coefficients (ICC) (2,1) showed strong agreement for the measurement of the T1 and T2 values (ICCs =⃥ 0.80 and 0.78 for T1 and T2 values, respectively).

CONCLUSION

The T1 and T2 values on MRF for the GM-HG were lower than those for the GM-STG and GM-MTG, likely reflecting a higher myelin content and iron deposition in the GM-HG. Quantitative measurements using the MRF can clarify cortical properties with high reliability, which may indicate that MRF mapping provides new insights into the structure of the human cortical GM.

Simultaneous relaxometry and morphometry of human brain structures with 3D magnetic resonance fingerprinting: a multicenter, multiplatform, multifield-strength study

Relaxation times and morphological information are fundamental magnetic resonance imaging-derived metrics of the human brain that reflect the status of the underlying tissue. Magnetic resonance fingerprinting (MRF) enables simultaneous acquisition of T1 and T2 maps inherently aligned to the anatomy, allowing whole-brain relaxometry and morphometry in a single scan. In this study, we revealed the feasibility of 3D MRF for simultaneous brain structure-wise morphometry and relaxometry. Comprehensive test-retest scan analyses using five 1.5-T and three 3.0-T systems from a single vendor including different scanner types across 3 institutions demonstrated that 3D MRF-derived morphological information and relaxation times are highly repeatable at both 1.5 T and 3.0 T. Regional cortical thickness and subcortical volume values showed high agreement and low bias across different field strengths. The ability to acquire a set of regional T1, T2, thickness, and volume measurements of neuroanatomical structures with high repeatability and reproducibility facilitates the ability of longitudinal multicenter imaging studies to quantitatively monitor changes associated with underlying pathologies, disease progression, and treatments.

MRI of focal cortical dysplasia

Focal cortical dysplasia (FCD) are histopathologically categorized in ILAE type I to III. Mild malformations of cortical development (mMCD) including those with oligodendroglial hyperplasia (MOGHE) are to be integrated into this classification yet. Only FCD type II have distinctive MRI and molecular genetics alterations so far. Subtle FCD including FCD type II located in the depth of a sulcus are often overlooked requiring the use of dedicated sequences (MP2RAGE, FLAWS, EDGE) and/or voxel (VBM)- or surface-based (SBM) postprocessing. The added value of 7 Tesla MRI has to be proven yet.

Optimized multi-axis spiral projection MR fingerprinting with subspace reconstruction for rapid whole-brain high-isotropic-resolution quantitative imaging

PURPOSE

To improve image quality and accelerate the acquisition of 3D MR fingerprinting (MRF).

METHODS

Building on the multi-axis spiral-projection MRF technique, a subspace reconstruction with locally low-rank constraint and a modified spiral-projection spatiotemporal encoding scheme called tiny golden-angle shuffling were implemented for rapid whole-brain high-resolution quantitative mapping. Reconstruction parameters such as the locally low-rank regularization parameter and the subspace rank were tuned using retrospective in vivo data and simulated examinations. B₀ inhomogeneity correction using multifrequency interpolation was incorporated into the subspace reconstruction to further improve the image quality by mitigating blurring caused by off-resonance effect.

RESULTS

The proposed MRF acquisition and reconstruction framework yields high-quality 1-mm isotropic whole-brain quantitative maps in 2 min at better quality compared with 6-min acquisitions of prior approaches. The proposed method was validated to not induce bias in T₁ and T₂ mapping. High-quality whole-brain MRF data were also obtained at 0.66-mm isotropic resolution in 4 min using the proposed technique, where the increased resolution was shown to improve visualization of subtle brain structures.

CONCLUSIONS

The proposed tiny golden-angle shuffling, MRF with optimized spiral-projection trajectory and subspace reconstruction enables high-resolution quantitative mapping in ultrafast acquisition time.

Three-dimensional high-resolution T1 and T2 mapping of whole macaque brain at 9.4 T using magnetic resonance fingerprinting

PURPOSE

Quantitative T₁ and T₂ mapping in non-human primates with whole-brain coverage is challenged by the requirement of sub-millimeter resolution and the inhomogeneity of the transmit magnetic field (B₁ ⁺ ) covering a large field of view. The goal of the current study is to develop a magnetic resonance fingerprinting (MRF) method for simultaneous T₁ and T₂ mapping of the entire macaque brain within feasible scan time.

METHODS

A three-dimensional (3D) MRF sequence with both inversion- and T₂ -preparation modules was developed and evaluated on a 9.4 T preclinical scanner. Data acquisition used a 3D stack-of-spirals trajectory, with undersampling along both the in-plane and the through-plane directions. The effect of B₁ ⁺ inhomogeneity was accounted for by matching the acquired fingerprint to a dictionary simulated with the B₁ ⁺ factors measured from a separate scan. In vitro and ex vivo studies were performed to evaluate the accuracy and the undersampling capacity of the MRF method. The application of the MRF method for in vivo, brain-wide T₁ and T₂ mapping was demonstrated on macaques at 4, 6, and 12 years of age.

RESULTS

The MRF method enabled highly repeatable T₁ and T₂ mapping at high spatial resolution (0.35 × 0.35 × 1 mm³ ) with an acceleration factor of 24. In vivo studies showed significant age-related T₂ reduction in deep gray nuclei including the globus pallidus, the putamen, and the caudate nucleus.

CONCLUSIONS

This study demonstrates the first MRF study for brain-wide, multi-parametric quantification in non-human primates with sub-millimeter resolution.

Magnetic Resonance Fingerprinting for Preoperative Meningioma Consistency Prediction

RATIONALE AND OBJECTIVES

Preoperative meningioma consistency prediction is highly beneficial for surgical planning and prognostication. We aimed to use magnetic resonance fingerprinting (MRF)-derived T1 and T2 values to preoperatively predict meningioma consistency.

MATERIALS AND METHODS

A total of 51 patients with meningiomas were enrolled in this study. MRF, T1-weighted imaging, T2-weighted imaging, and diffusion-weighted imaging were performed in all patients before surgery using a 3T MRI scanner. MRF-derived T1 and T2 values, T1-weightd and T2-weighted signal intensities, as well as apparent diffusion coefficient value yield from diffusion-weighted imaging were compared between the soft, moderate and hard meningiomas. Receiver operating characteristic curve analyses were used to determine the diagnostic performance of T1, T2 value, and a combination of T1 and T2 values.

RESULTS

After Bonferroni corrections, quantitative T1 and T2 values yielded from MRF were significantly different between the soft, moderate and hard meningiomas (all p < 0.05). T2 signal intensity was significantly different between the soft and hard, soft and moderate meningiomas (both p < 0.05), whereas was not significantly different between the moderate and hard meningiomas. However, T1 signal intensity and apparent diffusion coefficient value had no significant differences between the soft, moderate and hard meningiomas (all p > 0.05). The combination of T1 and T2 values had greater areas under receiver operating characteristic curve curves compared to individual T1 or T2 value.

CONCLUSION

MRF may help to preoperatively differentiate between the soft, moderate and hard meningiomas and may be useful in guiding the surgical planning.

Magnetic Resonance Fingerprinting of the Pediatric Brain

Magnetic resonance fingerprinting (MRF) is increasingly being used to evaluate brain development and differentiate normal and pathologic tissues in children. MRF can provide reliable and accurate intrinsic tissue properties, such as T1 and T2 relaxation times. MRF is a powerful tool in evaluating brain disease in pediatric population. MRF is a new quantitative MR imaging technique for rapid and simultaneous quantification of multiple tissue properties.

Magnetic resonance fingerprinting residual signals can disassociate human grey matter regions

The importance of accurate structural discrimination of the human grey matter regions has motivated the development of observer-independent reproducible methods that account for inter-individual architectonic variations. We introduce a non-invasive statistical residual analysis framework, employing unique tissue-specific magnetic resonance fingerprinting (MRF) signals after adjusting for the effect of T₁ and T₂^* MR relaxometry parameters (here termed MRF residuals). A 7 T Siemens MR scanner was used to acquire MRF signals, quantitative transmit magnetic field (B₁⁺) maps and T₁-weighted anatomical images of eleven cortical areas (5L, 5M, 5Ci, 7A, 7P, 7PC, hIP3, BA2, BA4a, BA4p and BA6) from six female participants. MRF residual signal for each voxel was calculated as the difference between the actual and best matching MRF signal evolutions from a precomputed MRF dictionary covering a range of T₁, T₂^* and B₁⁺ values. To compare MRF residuals between regions of interest, normalised autocorrelation was used as a shape-based statistical signal characterisation method and the Euclidean distance between autocorrelation profiles of residuals was used to measure the interareal dissimilarity. In the eleven cortical areas in both cerebral hemispheres of six participants, the proposed MRF residual analysis consistently showed interareal dissimilarity profiles that concorded with histological studies, indicating that MRF residuals potentially contain tissue microstructural information. MRF residual signals provide additional area-specific information that is complementary to the MR relaxometry-based (T₁, T₂^*) information used previously for distinguishing microstructural differences between human cerebral cortex regions in vivo. The proposed approach led to more accurate identification of structural variations across cortical areas of interest.

Technical overview of magnetic resonance fingerprinting and its applications in radiation therapy

Magnetic resonance fingerprinting (MRF) is an emerging imaging technique for rapid and simultaneous quantification of multiple tissue properties. The technique has been developed for quantitative imaging of different organs. The obtained quantitative measures have the potential to improve multiple steps of a typical radiotherapy workflow and potentially further improve integration of magnetic resonance imaging guided clinical decision making. In this review paper, we first provide a technical overview of the MRF method from data acquisition to postprocessing, along with recent development in advanced reconstruction methods. We further discuss critical aspects that could influence its usage in radiation therapy, such as accuracy and precision, repeatability and reproducibility, geometric distortion, and motion robustness. Finally, future directions for MRF application in radiation therapy are discussed.

Automated design of pulse sequences for magnetic resonance fingerprinting using physics-inspired optimization

Magnetic resonance fingerprinting (MRF) is a method to extract quantitative tissue properties such as [Formula: see text] and [Formula: see text] relaxation rates from arbitrary pulse sequences using conventional MRI hardware. MRF pulse sequences have thousands of tunable parameters, which can be chosen to maximize precision and minimize scan time. Here, we perform de novo automated design of MRF pulse sequences by applying physics-inspired optimization heuristics. Our experimental data suggest that systematic errors dominate over random errors in MRF scans under clinically relevant conditions of high undersampling. Thus, in contrast to prior optimization efforts, which focused on statistical error models, we use a cost function based on explicit first-principles simulation of systematic errors arising from Fourier undersampling and phase variation. The resulting pulse sequences display features qualitatively different from previously used MRF pulse sequences and achieve fourfold shorter scan time than prior human-designed sequences of equivalent precision in [Formula: see text] and [Formula: see text] Furthermore, the optimization algorithm has discovered the existence of MRF pulse sequences with intrinsic robustness against shading artifacts due to phase variation.

Repeatability of MR fingerprinting in normal cervix and utility in cervical carcinoma

Background

Magnetic resonance fingerprinting (MRF) is a fast-imaging acquisition technique that generates quantitative and co-registered parametric maps. The aim of this feasibility study was to evaluate the agreement between MRF and phantom reference values, scan-rescan repeatability of MRF in normal cervix, and its ability to distinguish cervical carcinoma (CC) from normal cervical tissues.

Methods

An International Society of Magnetic Resonance in Medicine/National Institute of Standards and Technology (ISMRM/NIST) phantom was scanned using MRF 15 times over 65 days. Agreement between MRF and phantom reference T1 and T2 values was assessed by linear regression. Healthy volunteers and patients with suspected CC were prospectively recruited. MRF was repeated twice for healthy volunteers (MRF1 and MRF2). Volumes of interest of normal cervical tissues and CC were delineated on T1 and T2 maps. MRF scan-rescan repeatability was evaluated by Bland-Altman plots, within-subject coefficients of variation (wCV), and intraclass correlation coefficients (ICC). T1 and T2 values were compared between CC and normal cervical tissues using Mann-Whitney U test. Receiver operating characteristic (ROC) analysis was performed to evaluate diagnostic efficiency.

Results

Strong correlations were observed between MRF and phantom (R²=0.999 for T1, 0.981 for T2). Twelve healthy volunteers (28.7±5.1 years) and 28 patients with CC (54.6±15.2 years) were recruited for the in-vivo experiments. Repeatability of MRF parameters were wCV <3% for T1, <5% for T2 and ICC ≥0.92 for T1, ≥0.94 for T2. T1 value of CC (1,529±112 ms) was higher than normal mucosa [MRF1: 1,430±129 ms, MRF2: 1,440±130 ms; P=0.031, area under the curve (AUC) ≥0.717] and normal stroma (MRF1: 1,258±101 ms, MRF2: 1,276±105 ms; P<0.001, AUC ≥0.946). T2 value of CC (69±9 ms) was lower than normal mucosa (MRF1: 88±16 ms, MRF2: 87±13 ms; P<0.001, AUC ≥0.854), but was not different from normal stroma (P=0.919).

Conclusions

Excellent agreement was observed between MRF and phantom reference values. MRF exhibited excellent scan-rescan repeatability in normal cervix with potential value in differentiating CC from normal cervical tissues.

Multimodal Image Integration for Epilepsy Presurgical Evaluation: A Clinical Workflow

Multimodal image integration (MMII) is a promising tool to help delineate the epileptogenic zone (EZ) in patients with medically intractable focal epilepsies undergoing presurgical evaluation. We report here the detailed methodology of MMII and an overview of the utility of MMII at the Cleveland Clinic Epilepsy Center from 2014 to 2018, exemplified by illustrative cases. The image integration was performed using the Curry platform (Compumedics Neuroscan^™, Charlotte, NC, USA), including all available diagnostic modalities such as Magnetic resonance imaging (MRI), Positron Emission Tomography (PET), single-photon emission computed tomography (SPECT) and Magnetoencephalography (MEG), with additional capability of trajectory planning for intracranial EEG (ICEEG), particularly stereo-EEG (SEEG), as well as surgical resection planning. In the 5-year time span, 467 patients underwent MMII; of them, 98 patients (21%) had a history of prior neurosurgery and recurring seizures. Of the 467 patients, 425 patients underwent ICEEG implantation with further CT co-registration to identify the electrode locations. A total of 351 patients eventually underwent surgery after MMII, including 197 patients (56%) with non-lesional MRI and 223 patients (64%) with extra-temporal lobe epilepsy. Among 269 patients with 1-year post-operative follow up, 134 patients (50%) had remained completely seizure-free. The most common histopathological finding is focal cortical dysplasia. Our study illustrates the usefulness of MMII to enhance SEEG electrode trajectory planning, assist non-invasive/invasive data interpretation, plan resection strategy, and re-evaluate surgical failures. Information presented by MMII is essential to the understanding of the anatomo-functional-electro-clinical correlations in individual cases, which leads to the ultimate success of presurgical evaluation of patients with medically intractable focal epilepsies.

Improving surgical outcome with electric source imaging and high field magnetic resonance imaging

While most patients with focal epilepsy present with clear structural abnormalities on standard, 1.5 or 3 T MRI, some patients are MRI-negative. For those, quantitative MRI techniques, such as volumetry, voxel-based morphometry, and relaxation time measurements can aid in finding the epileptogenic focus. High-field MRI, just recently approved for clinical use by the FDA, increases the resolution and, in several publications, was shown to improve the detection of focal cortical dysplasias and mild cortical malformations. For those cases without any tissue abnormality in neuroimaging, even at 7 T, scalp EEG alone is insufficient to delimitate the epileptogenic zone. They may benefit from the use of high-density EEG, in which the increased number of electrodes helps improve spatial sampling. The spatial resolution of even low-density EEG can benefit from electric source imaging techniques, which map the source of the recorded abnormal activity, such as interictal epileptiform discharges, focal slowing, and ictal rhythm. These EEG techniques help localize the irritative, functional deficit, and seizure-onset zone, to better estimate the epileptogenic zone. Combining those technologies allows several drug-resistant cases to be submitted to surgery, increasing the odds of seizure freedom and providing a must needed hope for patients with epilepsy.

Magnetic resonance fingerprinting: an overview

Magnetic resonance fingerprinting (MRF) is an evolving quantitative MRI framework consisting of unique data acquisition, processing, visualization, and interpretation steps. MRF is capable of simultaneously producing multiple high-resolution property maps including T1, T2, M0, ADC, and T2* measurements. While a relatively new technology, MRF has undergone rapid development for a variety of clinical applications from brain tumor characterization and epilepsy imaging to characterization of prostate cancer, cardiac imaging, among others. This paper will provide a brief overview of current state of MRF technology including highlights of technical and clinical advances. We will conclude with a brief discussion of the challenges that need to be overcome to establish MRF as a quantitative imaging biomarker.

Feasibility of Novel Three-Dimensional Magnetic Resonance Fingerprinting of the Prostate Gland: Phantom and Clinical Studies

Objective

To evaluate the feasibility of a new three-dimensional (3D) MR fingerprinting (MRF) technique for the prostate gland by conducting phantom and clinical studies.

Materials and Methods

The new 3D MRF technique used in this study enables quick data acquisition and has a high resolution. For the phantom study, the MRF T1 and T2 values in an in-house phantom were compared with those of gold-standard mapping methods using linear regression analysis. For the clinical study, we evaluated 90 patients who underwent prostate imaging with MRF for suspected prostate cancer between September 2019 and February 2020. The mean T1 and T2 values were compared in the peripheral zone, transition zone, and focal lesions using paired t tests. The differences in the T1 and T2 values according to cancer aggressiveness were evaluated using one-way analysis of variance.

Results

In the phantom study, the MRF T1 and T2 values showed a perfect correlation with the gold-standard T1 and T2 values (R > 0.99). In the clinical study, the T1 and T2 values in the peripheral zone were significantly higher than those in the transitional zone (p < 0.001, both). The T1 and T2 values in prostate cancer were significantly lower than those in the peripheral and transitional zones. The higher the grade of cancer, the lower the T2 values.

Conclusion

The T1 and T2 values obtained from the 3D MRF showed a perfect correlation with the gold standard values in the phantom study. Differences in the T1 and T2 values among the different zones of the prostate gland were identified using 3D MRF in patients.

Whole-brain 3D MR fingerprinting brain imaging: clinical validation and feasibility to patients with meningioma

Purpose

MR fingerprinting (MRF) is a MR technique that allows assessment of tissue relaxation times. The purpose of this study is to evaluate the clinical application of this technique in patients with meningioma.

Materials and methods

A whole-brain 3D isotropic 1mm³ acquisition under a 3.0T field strength was used to obtain MRF T₁ and T₂-based relaxometry values in 4:38 s. The accuracy of values was quantified by scanning a quantitative MR relaxometry phantom. In vivo evaluation was performed by applying the sequence to 20 subjects with 25 meningiomas. Regions of interest included the meningioma, caudate head, centrum semiovale, contralateral white matter and thalamus. For both phantom and subjects, mean values of both T₁ and T₂ estimates were obtained. Statistical significance of differences in mean values between the meningioma and other brain structures was tested using a Friedman’s ANOVA test.

Results

MR fingerprinting phantom data demonstrated a linear relationship between measured and reference relaxometry estimates for both T₁ (r² = 0.99) and T₂ (r² = 0.97). MRF T₁ relaxation times were longer in meningioma (mean ± SD 1429 ± 202 ms) compared to thalamus (mean ± SD 1054 ± 58 ms; p = 0.004), centrum semiovale (mean ± SD 825 ± 42 ms; p < 0.001) and contralateral white matter (mean ± SD 799 ± 40 ms; p < 0.001). MRF T₂ relaxation times were longer for meningioma (mean ± SD 69 ± 27 ms) as compared to thalamus (mean ± SD 27 ± 3 ms; p < 0.001), caudate head (mean ± SD 39 ± 5 ms; p < 0.001) and contralateral white matter (mean ± SD 35 ± 4 ms; p < 0.001)

Conclusions

Phantom measurements indicate that the proposed 3D-MRF sequence relaxometry estimations are valid and reproducible. For in vivo, entire brain coverage was obtained in clinically feasible time and allows quantitative assessment of meningioma in clinical practice.

Investigation of the effect of acquisition schemes on time-resolved magnetic resonance fingerprinting

Purpose.This study aims to investigate the feasibility of different acquisition methods for time-resolved magnetic resonance fingerprinting (TR-MRF) in computer simulation.Methods.An extended cardiac-torso (XCAT) phantom is used to generate abdominal T1, T2, and proton density maps for MRF simulation. The simulated MRF technique consists of an IR-FISP MRF sequence with spiral trajectory acquisition. MRF maps were simulated with different numbers of repetitions from 1 to 15. Three different methods were used to generate TR-MRF maps: (1) continuous acquisition without delay between MRF repetitions; (2) continuous acquisition with 5 s delay between MRF repetitions; (3) triggered acquisition with variable delay between MRF repetitions to allow the next acquisition to start at different respiration phase. After the generation of TR-MRF maps, the image quality indexes including the absolute T1 and T2 values, signal-to-noise-ratio (SNR), tumor-to-liver contrast-to-noise ratio, error in the amplitude of diaphragm motion and tumor volume error were used to evaluate the reconstructed parameter maps. Three volunteers were recruited to test the feasibility of the selected acquisition method.Results.Dynamic MR parametric maps using three different acquisition methods were estimated. The overall and liver T1 value error, liver SNR in T1 and T2 maps, and tumor SNR from T1 maps from triggered method is statistically significantly better than the other two methods (p-value < 0.05). The other image quality indexes have no significant difference between the triggered method and the other two continuous acquisition methods. All image quality indexes exhibit no significant difference between the acquisition methods with 0 s and 5 s delay. The triggered method was successfully performed in three healthy volunteers.Conclusion.TR-MRF technique was investigated using three different acquisition methods in computer simulation where the triggered method showed better performance than the other two methods. The triggered method has been tested successfully in healthy volunteers.

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

Purpose

Advanced MRI-based biomarkers offer comprehensive and quantitative information for the evaluation and characterization of brain tumors. In this study, we report initial clinical experience in routine glioma imaging with a novel, fully 3D multiparametric quantitative transient-state imaging (QTI) method for tissue characterization based on T1 and T2 values.

Methods

To demonstrate the viability of the proposed 3D QTI technique, nine glioma patients (grade II–IV), with a variety of disease states and treatment histories, were included in this study. First, we investigated the feasibility of 3D QTI (6:25 min scan time) for its use in clinical routine imaging, focusing on image reconstruction, parameter estimation, and contrast-weighted image synthesis. Second, for an initial assessment of 3D QTI-based quantitative MR biomarkers, we performed a ROI-based analysis to characterize T1 and T2 components in tumor and peritumoral tissue.

Results

The 3D acquisition combined with a compressed sensing reconstruction and neural network-based parameter inference produced parametric maps with high isotropic resolution (1.125 × 1.125 × 1.125 mm³ voxel size) and whole-brain coverage (22.5 × 22.5 × 22.5 cm³ FOV), enabling the synthesis of clinically relevant T1-weighted, T2-weighted, and FLAIR contrasts without any extra scan time. Our study revealed increased T1 and T2 values in tumor and peritumoral regions compared to contralateral white matter, good agreement with healthy volunteer data, and high inter-subject consistency.

Conclusion

3D QTI demonstrated comprehensive tissue assessment of tumor substructures captured in T1 and T2 parameters. Aiming for fast acquisition of quantitative MR biomarkers, 3D QTI has potential to improve disease characterization in brain tumor patients under tight clinical time-constraints.

Automated fusion of multimodal imaging data for identifying epileptogenic lesions in patients with inconclusive magnetic resonance imaging

Many methods applied to data acquired by various imaging modalities have been evaluated for their benefit in localizing lesions in magnetic resonance (MR) negative epilepsy patients. No approach has proven to be a stand-alone method with sufficiently high sensitivity and specificity. The presented study addresses the potential benefit of the automated fusion of results of individual methods in presurgical evaluation. We collected electrophysiological, MR, and nuclear imaging data from 137 patients with pharmacoresistant MR-negative/inconclusive focal epilepsy. A subgroup of 32 patients underwent surgical treatment with known postsurgical outcomes and histopathology. We employed a Gaussian mixture model to reveal several classes of gray matter tissue. Classes specific to epileptogenic tissue were identified and validated using the surgery subgroup divided into two disjoint sets. We evaluated the classification accuracy of the proposed method at a voxel-wise level and assessed the effect of individual methods. The training of the classifier resulted in six classes of gray matter tissue. We found a subset of two classes specific to tissue located in resected areas. The average classification accuracy (i.e., the probability of correct classification) was significantly higher than the level of chance in the training group (0.73) and even better in the validation surgery subgroup (0.82). Nuclear imaging, diffusion-weighted imaging, and source localization of interictal epileptic discharges were the strongest methods for classification accuracy. We showed that the automatic fusion of results can identify brain areas that show epileptogenic gray matter tissue features. The method might enhance the presurgical evaluations of MR-negative epilepsy patients.