Next Article in Journal
DNA Methylation of Imprinted Genes KCNQ1, KCNQ1OT1, and PHLDA2 in Peripheral Blood Is Associated with the Risk of Breast Cancer
Next Article in Special Issue
Brain Functional Connectivity in Low- and High-Grade Gliomas: Differences in Network Dynamics Associated with Tumor Grade and Location
Previous Article in Journal
Comment on Colletti et al. Prevalence and Management of Cancer of the Rectal Stump after Total Colectomy and Rectal Sparing in Patients with Familial Polyposis: Results from a Registry-Based Study. Cancers 2022, 14, 298
Previous Article in Special Issue
Quantifying Tumor Heterogeneity via MRI Habitats to Characterize Microenvironmental Alterations in HER2+ Breast Cancer
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 14
Issue 11
10.3390/cancers14112651
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantitative Synthetic Magnetic Resonance Imaging for Brain Metastases: A Feasibility Study

by
Amaresha Shridhar Konar
1,†,
Akash Deelip Shah
2,†,
Ramesh Paudyal
1,
Maggie Fung
3,
Suchandrima Banerjee
3,
Abhay Dave
4,
Vaios Hatzoglou
2,‡ and
Amita Shukla-Dave
1,2,*,‡
1
Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York City, NY 10065, USA
2
Department of Radiology, Memorial Sloan Kettering Cancer Center, New York City, NY 10065, USA
3
General Electric Health Care, New York City, NY 10065, USA
4
Touro College of Osteopathic Medicine, New York City, NY 10027, USA
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
‡
These authors contributed equally to this work.
Cancers 2022, 14(11), 2651; https://0-doi-org.brum.beds.ac.uk/10.3390/cancers14112651
Submission received: 31 March 2022 / Revised: 16 May 2022 / Accepted: 20 May 2022 / Published: 27 May 2022
(This article belongs to the Special Issue Characterization of Tumor Physiology Using Magnetic Resonance Imaging (MRI))
Browse Figures
Versions Notes

Abstract

:

Simple Summary

This preliminary study aims to characterize brain metastases (BM) using T1 and T2 maps generated from newer, rapid, synthetic MRI (MAGnetic resonance image Compilation; MAGiC) in a clinical setting. In addition, synthetic MR could provide contrast images analogous to standard T1- and T2-weighted images. The reproducibility and repeatability of this method have been previously established for brain imaging. This study reports and analyzes the quantitative T1 and T2 values for 11 BM patients (17 BM lesions) with a total of 82 regions of interest (ROIs) delineated by an experienced neuroradiologist. The initial results, which need to be further validated in a larger patient cohort, demonstrated the ability of T1 and T2 metric values to characterize BMs and normal-appearing brain tissues. The T1 and T2 metrics could be potential surrogate biomarkers for BM free water content (cellularity) and tumor morphology, respectively.

Abstract

The present preliminary study aims to characterize brain metastases (BM) using T1 and T2 maps generated from newer, rapid, synthetic MRI (MAGnetic resonance image Compilation; MAGiC) in a clinical setting. We acquired synthetic MRI data from 11 BM patients on a 3T scanner. A multiple-dynamic multiple-echo (MDME) sequence was used for data acquisition and synthetic image reconstruction, including post-processing. MDME is a multi-contrast sequence that enables absolute quantification of physical tissue properties, including T1 and T2, independent of the scanner settings. In total, 82 regions of interest (ROIs) were analyzed, which were obtained from both normal-appearing brain tissue and BM lesions. The mean values obtained from the 48 normal-appearing brain tissue regions and 34 ROIs of BM lesions (T1 and T2) were analyzed using standard statistical methods. The mean T1 and T2 values were 1143 ms and 78 ms, respectively, for normal-appearing gray matter, 701 ms and 64 ms for white matter, and 4206 ms and 390 ms for cerebrospinal fluid. For untreated BMs, the mean T1 and T2 values were 1868 ms and 100 ms, respectively, and 2211 ms and 114 ms for the treated group. The quantitative T1 and T2 values generated from synthetic MRI can characterize BM and normal-appearing brain tissues.

1. Introduction

Magnetic resonance imaging (MRI) is one of the most widely employed imaging modalities for assessing brain tumors, owing to its noninvasive provision of high-resolution structural (T1- (spin–lattice) and T2- (spin–spin) weighted) images and functional (diffusion, perfusion, and metabolic) information [1]. Qualitative multiplanar T1 and T2 images (in particular post-Gadolinium T1-weighted) have been widely used for the anatomical location and morphological characteristics in brain tumors and metastases [2,3]. The varying signal intensities of T1 and T2 are based on the differences in the brain tissue’s intrinsic relaxation times associated with acquisition parameters (e.g., flip angle, echo time, and repetition time) and can be quantitatively mapped to measure pathophysiological changes [4,5]. The measurements of quantitative (q) relaxometry metrics may serve as surrogate biomarkers for free water content (cellularity) [6] and tumor morphology [7]. In clinical practice, standard quantitative T1 and T2 acquisition methods are time-consuming and not practically feasible in busy clinics [8,9,10,11,12]. Furthermore, the maps can be challenging to interpret due to low spatial resolution [13]. A single, time-efficient acquisition method that could simultaneously measure multiple tissue properties would account for these limitations and is therefore an unmet need in clinical brain imaging.
Rapid MRI acquisition methods have been developed to address the above technical limitations, showing promise in brain imaging [14,15,16]. Synthetic MRI is a newer, rapid method that simultaneously offers qualitative images and quantitative T1 and T2 maps in a clinically feasible timeframe [15,17,18]. The computational viability of data acquisition and processing to implement synthetic MRI was first tested for the brain [17]. The synthetic MR images were mathematically derived to display image contrast analogous to standard T1- and T2-weighted images [19]. The reproducibility and repeatability of this method have been vigorously tested and established for brain imaging [20]. A major vendor, General Electric Healthcare (GEHC, Waukesha, WI, USA), has developed synthetic MRI and termed the product MAGiC (MAGnetic resonance image Compilation), which is used in this study [21].
The quality and diagnostic accuracy of synthetic contrast-weighted images have been compared with standard images in the brain [21,22]. Tanenbaum et al. performed a prospective trial on 109 subjects with neuroimaging indications, in which 1526 MR images were read by seven blinded neuroradiologists to compare synthetic qualitative images with standard brain MRI [21]. They concluded that the synthetic MR imaging quality was similar to that of standard proton density (PD), short tau inversion recovery (STIR), and T1- and T2-weighted contrast views across neurologic conditions. While artifacts were more common in synthetic T2 fluid attenuated inversion recovery (FLAIR), these were readily recognizable, did not mimic pathology, and only necessitated additional standard T2 FLAIR to confirm diagnosis [21]. This study extensively focused on qualitative image analysis, not quantitative relaxometry maps [21]. Hagiwara et al. studied ten patients with a combined total of 167 brain metastases (BM) lesions and reported that the T1 inversion recovery (IR) qualitative images generated by the synthetic MRI method in BM patients created better contrast than synthetic T1-weighted or standard T1 IR imaging [23]. Furthermore, the detection of brain metastases was comparable among these qualitative images [23].
Recently, Blystad et al. demonstrated the promise of quantitative synthetic MRI to assess nonvisible peritumoral contrast enhancement in the peritumoral edema of malignant gliomas via the change (Δ) in R1 (R1 = 1/T1) maps [24]. In accordance, Müller et al. previously showed that standard quantitative T1 mapping could discover a subtle enhancement in glioblastoma patients, which is not visually detectable with T1-weighted subtraction images [25]. These findings have clinical implications, as an early reduction in the enhancing region under therapy predicts a favorable therapy response [24,25].
We will focus on BM as it is the most common intracranial tumor in adults. Lung cancer, breast cancer, and melanoma are the most common primary tumor sites [26]. The incidence of BM has risen due to the increased rates of diagnosis and increased survival of cancer patients [27,28]. Multiparametric MRI, including quantitative relaxometry metrics, may further improve detection and treatment selection [29,30,31]. A recently published study measured the magnetic resonance relaxation time using a multi-dynamic multi-echo (MDME) sequence at three time points (1, 10, and 20 min) after contrast injection on seven BM patients and showed that these measurements were time-dependent [32]. The present preliminary study aims to investigate the value of T1 and T2 metric values derived from the newer, rapid synthetic MRI method (MAGiC) and assess their ability to characterize untreated and treated BM, as well as normal-appearing brain tissue.

2. Materials and Methods

2.1. Phantom Selection

An MRI system phantom procured from CaliberMRI (Boulder, CO, USA) was co-developed by International Society for Magnetic Resonance in Medicine (ISMRM)/National Institute of Standards and Technology (NIST) to provide a standard phantom to perform qMRI studies, including and not limited to T1 and T2 relaxation times [33]. MRI system phantom consists of 14 vials in each T1 and T2 array with appropriate chemical composition to provide precise T1 and T2 values, including the range of T1 and T2 values generally seen in gray and white matter of the brain. The Nuclear Magnetic Resonance (NMR)-based reference T1 and T2 relaxation times were provided by NIST. This phantom has SI-traceable components and was monitored for its stability and accuracy [34].

2.2. Patient Selection

The Institutional Review Board (IRB) approved this Health Insurance Portability and Accountability Act (HIPAA)-compliant prospective study for participating patients. Written informed consent was obtained from all eligible 21 patients with brain tumors. This study focuses on 14 patients who had BM. In final analysis, we investigated the data acquired from 11 patients who had a total of 17 evaluable BM lesions. We excluded three patients due to small lesion size <0.5 mm. There is a subset of patient overlap with our companion paper on Brain Magnetic Resonance Fingerprinting (MRF) in this same journal issue. MRF and synthetic MRI are two different MRI relaxometry techniques. There is no overlap in image acquisition and data analysis methods between these two studies.
The median age for these 11 BM patients was 52 years (range, 25–61 years; five male and six female). All BM patients were enrolled between June 2019 and August 2021. The inclusion criteria were age ≥ 18 years, and clinical or radiological diagnosis of BM. Patient characteristics are illustrated in Table 1.
Patients with BM underwent standard MRI, including MAGiC sequence, irrespective of the treatment group (untreated (n = 3) or treated (n = 8)). This study focuses on testing the MAGiC sequence in BM patients. They were scanned at a single time point, not longitudinally. The therapy regimens for eight BM patients who underwent treatment were stereotactic radiosurgery (SRS) n = 3, focal radiation therapy (RT) n = 2, and whole-brain RT n = 3.

2.3. MRI Data Acquisition

ISMRM/NIST system phantom was scanned on GE (General Electric Healthcare, Waukesha, WI, USA) MRI system (Discovery 3.0 T MR750w) using an eight-channel brain array coil. The Gold Standard (GS) T1 measurements from the T1 arrays were acquired by the IR spin echo method with specific acquisition parameters as follows: inversion time (TI) = 50, 75, 100, 125, 150, 250, 500, 1000, 2000, 3000 ms; repetition time (TR) = 4500 ms; echo time (TE) = 7.34 ms; acquisition matrix 128 × 128; matrix reconstructed to 256 × 256; field of view (FOV) = 25 cm; slice thickness = 5 mm. The scan time for each TI measurement was approximately four minutes and the total scan time for GS T1 acquisition was around 40 min. The GS T2 measurements from the T2 array were obtained using a multiple single-echo spin echo method with the following acquisition parameters: TEs = 9, 12, 15, 20, 25, 30, 40, 45, 50, 60, 75, 80, 100, 120, 160 ms; TR = 5000 ms; acquisition matrix 128 × 128; matrix reconstructed to 256 × 256; FOV = 25 cm; slice thickness = 5 mm. The scan time of each TE measurement was approximately 21: 30 (min: sec); the total scan time was approximately five hours. The nonlinear least-squares curve fitting was performed in MATLAB (The MathWorks. Inc., Natick, MA, USA). MAGiC phantom data were acquired in a coronal plane with an FOV of 25 cm. Remaining acquisition parameters used for the MAGiC phantom data acquisition were the same as the patient data acquisition detailed below.
All patient MRI examinations were performed on a GE MRI system (Discovery 3.0 T MR750w) using a 16-channel head and neck coil for brain imaging. Images were prospectively acquired using a fixed set of scanning parameters closely approximating the current standard-of-care for the brain. The standard clinical MR acquisition parameters for brain imaging was multiplanar (axial, coronal, and sagittal). In this study, we used conventional axial T1-weighted (w) (both pre- and post-contrast images), T2w fat-suppressed, and T2w FLAIR images were acquired pre-contrast with slice thickness (ST) of 3 mm and FOV of 20–24 cm. The acquisition parameters were as follows: T1w imaging: TR = 2000 ms, TI = 1101 ms; TE = 25 ms; number of averages (NA) = 1; acquisition matrix 320 × 224; matrix reconstructed to 320 × 256; scan time ~2.43 min. T2w imaging with fat-suppressed fast spin echo: TR = 4796 ms; TE = 121 ms; NA = 1; acquisition matrix 256 × 256; matrix reconstructed to 320 × 256; scan time ~2.27 min. FLAIR imaging: TR = 9946 ms; TE = 127 ms; TI = 2375 ms; NA = 1, acquisition matrix 256 × 192; matrix reconstructed to 256 × 256, scan time ~5.13 min.
MAGiC brain data were acquired prior to contrast agent injection using multiple-dynamic multiple-echo (MDME) sequence, followed by synthetic image reconstruction and post-processing. MDME is a multi-contrast sequence and uses a multi-echo acquisition that enables absolute quantification of physical tissue properties, including T1, T2, and proton density (PD) independent of scanner settings. It uses multiple inversion times (TIs) to measure T1 relaxation. MDME parameters acquired in one scan are used in synthetic imaging to calculate pixel intensity and produce an appearance similar to standard MR images with modifiable TE, TR, and TI [35]. This provides quantitative (T1, T2, and PD maps) and qualitative images (T1, T2, T1 FLAIR, T2 FLAIR, PD, and STIR contrast images). Two-dimensional axial MAGiC brain data were acquired with four automatically calculated saturation delays, two echo times (TE) (23.4 ms and effective TE of 93.8 ms), and an auto TR of 4000 ms (ranging between 4000 ms–15,000 ms). The other parameters were FA = 90°; Echo Train Length (ETL) = 12; ST = 5 mm; FOV = 25 cm; acquisition matrix 320 × 256; matrix reconstructed to 256 × 256. Total scan time was approximately 6 min.

2.4. MRI Data Post-Processing

MAGiC is a combined package of image acquisition and a post-processing. MAGiC post-processing generates parametric maps of both the T1 and T2 relaxation times, as well as PD, using the data obtained by running the MDME sequence. Contrast-weighted T1, T2, T1 FLAIR, T2 FLAIR, STIR, phase-sensitive inversion recovery (PSIR), and double inversion recovery (DIR) images can be synthesized using the quantitative maps for any desired combination of TE, TR, and TI. The contrast-weighted image settings can be defined on the post-process screen or the MAGiC session application. This study focuses on the estimation of the quantitative T1 and T2 parametric maps for BM patients.

2.5. Regions of Interest Delineation

Regions of interest (ROI) were delineated by a neuroradiologist with over seven years of experience. To determine the extent of the tumor, anatomical T2w and contrast-enhanced T1w images were used. ROIs were drawn, excluding cystic and necrotic regions. Based on the literature, we modified the lesions’ size threshold for ROI delineation to ≥5 mm [36]. The size of the BM lesions ranged between 5 mm and 26 mm in this study. MAGiC-derived mean T1 and T2 metric values were obtained from the selected ROIs for normal and tumor tissues.
The neuroradiologist delineated 17 ROIs on the central slice of the evaluable BM lesions. ROIs were also delineated in normal-appearing gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF) for the eight patients (three untreated; five treated with SRS or focal RT). We did not assess the GM, WM, and CSF for the three patients who received whole-brain RT as there can be microscopic changes in the normal-appearing tissues after radiation exposure. MAGiC simultaneously generates T1 and T2 maps. Hence, ROIs were drawn on the T1 map and the same ROIs were used for T2 maps to obtain the mean value. In total, 82 ROIs were analyzed for all patients in this study.

2.6. Statistical Analysis

For the T1 and T2 values obtained using the phantom, the relative percentage difference between the method (phantom vendor provided (VP), GS acquisition, and the MAGiC) was calculated. The univariate analysis was performed on the mean metric value extracted from the ROIs using the Wilcoxon Rank Sum test (WRST) to find the difference between the untreated and treated groups. WRST is also called the Wilcoxon Mann–Whitney test and the Mann–Whitney U test to compare two independent samples. The T1 and T2 values from both untreated and treated groups were compared with the healthy-appearing brain tissues. In addition, T1 and T2 values (mean, median, and range) were reported from healthy-appearing GM, WM, and CSF. For WRST, the significance level was set at p ≤ 0.05.

3. Results

Table 2 reports the T1 and T2 values, as well as the relative percentage difference, between the three methods (VP, GS, and MAGiC). The percentage difference between the methods for T1 values showed a maximum of 6.6% between VP and GS, 16% between VP and MAGiC, and 23.9% between GS and MAGiC. Similarly, the percentage difference between the methods for T2 values showed a maximum of 16.8% between VP and GS, 18.4% between VP and MAGiC, and 30.5% between GS and MAGiC.
Figure 1 shows qualitative (standard clinically acquired T2w and post-contrast T1w images) and quantitative (T1 and T2 maps generated from MAGiC for normal-appearing GM, WM, and CSF) MRI data from a representative BM patient treated with focal RT. The mean T1 values for normal-appearing GM, WM, and CSF regions were 1098 ms, 767 ms, and 4289 ms, respectively. The mean T2 values were 75 ms, 65 ms, and 394 ms, respectively.
Figure 2 shows representative qualitative, multi-contrast-weighted images that were synthesized at the fixed combination of TE, TR, and TI from a BM patient treated with whole-brain RT. This figure illustrates that multiple diagnostic images can be obtained in a rapid, single MRI acquisition in a clinically feasible time.
Figure 3 shows the representative quantitative T1 and T2 maps from the same patient as Figure 2. The mean T1 and T2 values measured within the BM lesions were 1933 ms and 96 ms, respectively.
The box plot in Figure 4 exhibits T1 and T2 values obtained from normal-appearing GM, WM, and CSF regions using MAGiC-generated T1 and T2 maps. GM mean values for T1 and T2 were 1143 ms and 78 ms, respectively. Mean T1 and T2 values for WM were 701 ms and 64 ms, respectively. Mean T1 and T2 values for CSF regions were 4206 ms and 390 ms, respectively.
Table 3 shows the mean and range for T1 and T2 values estimated from MAGiC for BM lesions (untreated n = 3; treated n = 8). Mean T1 and T2 values measured within BM lesions for the untreated group were 1868 ms and 100 ms, respectively. Similarly, the treated group’s mean T1 and T2 values were 2211 ms and 114 ms, respectively. In comparison, the WRST showed no significant difference for the T1 and T2 values from BM lesions between the treated and untreated groups (p > 0.05). However, WRST performed for the mean T1 and T2 values between the BM and normal-appearing WM showed significant differences for both the untreated and treated groups (p < 0.05).
Figure 5 depicts a heat map of the range of T1 and T2 values from all the BM lesions, exhibiting a clear difference between BM and normal-appearing brain tissue (GM and WM). All 17 metastases from 11 patients were analyzed separately (patients were coded as “P” and metastatic lesions as “BM”) and four patients had more than one metastasis.

4. Discussion

Quantitative MRI using relaxometry reveals fundamental information and absolute numerical values of tissue [37]. However, the long acquisition times associated with quantitative T1 and T2 measurements limit its clinical applications [8,9,10,11,12]. Reducing acquisition scan times can facilitate diagnosis, staging, and treatment selection [38]. Riederer et al. illustrated that synthesized MR images can be generated using arbitrary pulse-sequence parameters [15,39]. Feasibility studies have previously reported that synthetic MRI-generated T1 and T2 maps were diagnostically acceptable in the normal brain [21,40,41,42]. The utility of this method has been further tested for extracranial organs, including the prostate, breast, and rectum [43,44,45,46]. The present study tested the robustness of the MAGiC method for the quantitative measurement of T1 and T2 values on a phantom and patients with BM. The accurate measurement of T1 and T2 values is critical in qMRI applications to probe tissue physiology. Measured Gold Standard (GS) T1 and T2 values would provide further evidence of the robustness of relaxometry mapping performed using the MAGiC method. GS acquisition is time-consuming and practically challenging to perform on patients. Therefore, we performed an additional experiment using the ISMRM/NIST MRI system phantom. We evaluated the robustness of MAGiC T1 and T2 measurements on the phantom by comparing vendor-provided data and data acquired with the GS method. We considered the T1 and T2 values from only seven vials that mimic brain tissue [9,47].
In our preliminary experience, synthetic MRI-generated T1 and T2 maps exhibited tumor heterogeneity for the first time and the metric values were able to differentiate between normal-appearing tissue and BM lesions. Using a 3T MRI and standard T1 and T2 mapping methods, Wansapura et al. reported mean T1 values of 1331 ms and 832 ms, and mean T2 values of 110 ms and 80 ms for normal-appearing GM and WM, respectively [9]. Our quantitative synthetic MRI results reported similar T1 and T2 values for normal-appearing GM and WM as in their published study [9].
Intrinsic T1 and T2 image contrast depends on the composition of tissues, including tissue water content, iron concentration, and the relative proportion of various macromolecules [48]. The quantification of T1 and T2 provides surrogate metrics of the underlying pathophysiology with tissue measurement of free water content and vascular morphology [12]. Quantitative T1 and T2 values may reflect an alteration in tissue composition and are surrogate biomarkers for tumor water content and morphology. These differences can explain variations in macromolecular tissue compositions. By adjusting the MR scan parameters for tissues, we can improve our ability to capture these changes. Quantifying relaxation times and proton density by the multi-echo acquisition of a saturation-recovery using a turbo spin echo readout sequence on a 1.5T MRI, Warntjes et al. reported T1 and T2 values of 561 ms and 73 ms, respectively, in the frontal white matter of the adult brain [35]. In the present study, the newer, rapid MAGiC method synthetically generated multiple contrast-weighted images at a fixed combination of TE, TR, and TI for all BM patients. As Tanenbaum et al. have exhaustively demonstrated the similar image quality of multiple synthetically (MAGiC) generated images to that of diagnostic MRI across neurologic conditions [21], we did not report here on qualitative diagnostic comparisons for brain imaging. Tanenbaum et al. also pointed out that, though artifacts were more common in synthetic T2 FLAIR, these were readily recognizable and did not mimic pathology, but they could necessitate additional standard T2 FLAIR to confirm diagnosis [21]. It has also been previously reported that BM detection was comparable among qualitative synthetic and standard MR images [21].
The implementation of advanced quantitative imaging techniques requires significant time, technical, and cost resources. Therefore, the quantitative characterization of tissue properties eludes effective clinical utility in day-to-day practice. Synthetic MRI-generated multiple contrast images, including qualitative and quantitative T1 and T2 mapping, represent an alternative solution for clinical diagnostic imaging due to their rapid acquisition and data quantification. In addition, these newly developed methods hold promise for the identification of subtle enhancements in the peritumoral area that are not visible on standard imaging. These findings could have clinical implications as the early reduction in enhancing regions under therapy predicts a favorable therapy response [24,25].
The focus of this study was to test the ability of quantitative T1 and T2 values to characterize BM lesions and normal-appearing brain tissues. This study had several limitations. The sample size was small for both treated and untreated groups. The results, therefore, need to be validated in studies with a larger BM patient cohort. The slightly higher T1 relaxation times observed for untreated and treated BMs could be due in part to pathophysiological variations caused by differences in the primary tumor sites. This needs further investigation. The measured T2 value for CSF using MAGiC was underestimated, which may be attributed to factors such as CSF motion, eddy current effects, inadequate sampling of the relaxation curve, partial volume, and therefore warrants further investigation. The role of both synthetic qualitative and quantitative images in the diagnostic workup requires multi-center testing before incorporation into rapid brain imaging protocols for clinical use. Future studies may also evaluate synthetically-generated multiple contrast images using MAGiC for BM patients.

5. Conclusions

In conclusion, the present preliminary results demonstrate that quantitative T1 and T2 values from synthetic MRI can differentiate between BM and normal-appearing brain tissues. Further evaluation in a larger patient cohort is required to characterize BM’s pathophysiological features and to assess the longitudinal treatment response.

Author Contributions

Conceptualization, A.S.K., A.D.S., V.H. and A.S.-D.; methodology, A.S.K., A.D.S., R.P., V.H. and A.S.-D.; validation, A.S.K., A.D.S., R.P., V.H. and A.S.-D.; formal analysis, A.S.K. and A.D.S.; investigation, A.S.K., A.D.S., R.P., A.D., V.H. and A.S.-D.; resources, A.S.K., A.D.S., V.H. and A.S.-D.; data curation, A.S.K., A.D.S., R.P., V.H. and A.S.-D.; writing—original draft preparation, A.S.K., A.D.S., R.P., M.F., S.B., A.D., V.H. and A.S.-D.; writing—review and editing, A.S.K., A.D.S., R.P., M.F., S.B., A.D., V.H. and A.S.-D.; visualization, A.S.K., A.D.S., R.P., V.H. and A.S.-D.; supervision V.H. and A.S.-D.; funding acquisition, MSKCC grant. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by NIH, P30 CA008748.

Institutional Review Board Statement

The Institutional Review Board approved the study (IRB protocol # 13-239, approved in 2013).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study will be provided upon reasonable request.

Acknowledgments

We are thankful to Cecile Berberat for editing the full manuscript, and to MR technologists for helping us in acquiring the datasets.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sawlani, V.; Patel, M.D.; Davies, N.; Flintham, R.; Wesolowski, R.; Ughratdar, I.; Pohl, U.; Nagaraju, S.; Petrik, V.; Kay, A.; et al. Multiparametric MRI: Practical approach and pictorial review of a useful tool in the evaluation of brain tumours and tumour-like lesions. Insights Imaging 2020, 11, 84. [Google Scholar] [CrossRef] [PubMed]
  2. Mehrabian, H.; Detsky, J.; Soliman, H.; Sahgal, A.; Stanisz, G.J. Advanced Magnetic Resonance Imaging Techniques in Management of Brain Metastases. Front. Oncol. 2019, 9, 440. [Google Scholar] [CrossRef] [PubMed]
  3. Tong, E.; McCullagh, K.L.; Iv, M. Advanced Imaging of Brain Metastases: From Augmenting Visualization and Improving Diagnosis to Evaluating Treatment Response. Front. Neurol. 2020, 11, 270. [Google Scholar] [CrossRef]
  4. Deoni, S.C.; Peters, T.M.; Rutt, B.K. Determination of optimal angles for variable nutation proton magnetic spin-lattice, T1, and spin-spin, T2, relaxation times measurement. Magn. Reson. Med. 2004, 51, 194–199. [Google Scholar] [CrossRef]
  5. Ngo, F.Q.; Bay, J.W.; Kurland, R.J.; Weinstein, M.A.; Hahn, J.F.; Glassner, B.J.; Woolley, C.A.; Dudley, A.W., Jr.; Ferrario, C.M.; Meaney, T.F. Magnetic resonance of brain tumors: Considerations of imaging contrast on the basis of relaxation measurements. Magn. Reson. Imaging 1985, 3, 145–155. [Google Scholar] [CrossRef]
  6. McSheehy, P.M.; Weidensteiner, C.; Cannet, C.; Ferretti, S.; Laurent, D.; Ruetz, S.; Stumm, M.; Allegrini, P.R. Quantified tumor T1 is a generic early-response imaging biomarker for chemotherapy reflecting cell viability. Clin. Cancer Res. 2010, 16, 212–225. [Google Scholar] [CrossRef] [Green Version]
  7. Kim, E.; Lee, E.; Plummer, C.; Gil, S.; Popel, A.S.; Pathak, A.P. Vasculature-specific MRI reveals differential anti-angiogenic effects of a biomimetic peptide in an orthotopic breast cancer model. Angiogenesis 2015, 18, 125–136. [Google Scholar] [CrossRef] [Green Version]
  8. Look, D.C.; Locker, D.R. Time saving in measurement of NMR and EPR relaxation times. Rev. Sci. Instrume 1970, 41, 250–251. [Google Scholar] [CrossRef] [Green Version]
  9. Wansapura, J.P.; Holland, S.K.; Dunn, R.S.; Ball, W.S., Jr. NMR relaxation times in the human brain at 3.0 tesla. J. Magn. Reson. Imaging 1999, 9, 531–538. [Google Scholar] [CrossRef]
  10. Gelman, N.; Ewing, J.R.; Gorell, J.M.; Spickler, E.M.; Solomon, E.G. Interregional variation of longitudinal relaxation rates in human brain at 3.0 T: Relation to estimated iron and water contents. Magn. Reson. Med. 2001, 45, 71–79. [Google Scholar] [CrossRef]
  11. Schmitt, P.; Griswold, M.A.; Jakob, P.M.; Kotas, M.; Gulani, V.; Flentje, M.; Haase, A. Inversion recovery TrueFISP: Quantification of T1, T2, and spin density. Magn. Reson. Med. 2004, 51, 661–667. [Google Scholar] [CrossRef] [PubMed]
  12. Deoni, S.C.; Rutt, B.K.; Peters, T.M. Rapid combined T1 and T2 mapping using gradient recalled acquisition in the steady state. Magn. Reson. Med. 2003, 49, 515–526. [Google Scholar] [CrossRef] [PubMed]
  13. Deoni, S.C.; Peters, T.M.; Rutt, B.K. High-resolution T1 and T2 mapping of the brain in a clinically acceptable time with DESPOT1 and DESPOT2. Magn. Reson. Med. 2005, 53, 237–241. [Google Scholar] [CrossRef]
  14. Riederer, S.J.; Lee, J.N.; Farzaneh, F.; Wang, H.Z.; Wright, R.C. Magnetic resonance image synthesis. Clinical implementation. Acta Radiol. Suppl. 1986, 369, 466–468. [Google Scholar] [PubMed]
  15. Riederer, S.J.; Suddarth, S.A.; Bobman, S.A.; Lee, J.N.; Wang, H.Z.; MacFall, J.R. Automated MR image synthesis: Feasibility studies. Radiology 1984, 153, 203–206. [Google Scholar] [CrossRef]
  16. Glad, I.K.; Sebastiani, G. A Bayesian approach to synthetic magnetic resonance imaging. Biometrika 1995, 82, 237–250. [Google Scholar] [CrossRef]
  17. Maitra, R.; Riddles, J.J. Synthetic magnetic resonance imaging revisited. IEEE Trans. Med. Imaging 2010, 29, 895–902. [Google Scholar] [CrossRef] [Green Version]
  18. Blystad, I.; Warntjes, J.B.; Smedby, O.; Landtblom, A.M.; Lundberg, P.; Larsson, E.M. Synthetic MRI of the brain in a clinical setting. Acta Radiol. 2012, 53, 1158–1163. [Google Scholar] [CrossRef]
  19. Andica, C.; Hagiwara, A.; Hori, M.; Kamagata, K.; Koshino, S.; Maekawa, T.; Suzuki, M.; Fujiwara, H.; Ikeno, M.; Shimizu, T.; et al. Review of synthetic MRI in pediatric brains: Basic principle of MR quantification, its features, clinical applications, and limitations. J. Neuroradiol. 2019, 46, 268–275. [Google Scholar] [CrossRef]
  20. Hagiwara, A.; Warntjes, M.; Hori, M.; Andica, C.; Nakazawa, M.; Kumamaru, K.K.; Abe, O.; Aoki, S. SyMRI of the Brain: Rapid Quantification of Relaxation Rates and Proton Density, With Synthetic MRI, Automatic Brain Segmentation, and Myelin Measurement. Investig. Radiol. 2017, 52, 647–657. [Google Scholar] [CrossRef]
  21. Tanenbaum, L.N.; Tsiouris, A.J.; Johnson, A.N.; Naidich, T.P.; DeLano, M.C.; Melhem, E.R.; Quarterman, P.; Parameswaran, S.X.; Shankaranarayanan, A.; Goyen, M.; et al. Synthetic MRI for Clinical Neuroimaging: Results of the Magnetic Resonance Image Compilation (MAGiC) Prospective, Multicenter, Multireader Trial. Am. J. Neuroradiol. 2017, 38, 1103–1110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Di Giuliano, F.; Minosse, S.; Picchi, E.; Marfia, G.A.; Da Ros, V.; Muto, M.; Muto, M.; Pistolese, C.A.; Laghi, A.; Garaci, F. Comparison between synthetic and conventional magnetic resonance imaging in patients with multiple sclerosis and controls. Magn. Reson. Mater. Phys. Biol. Med. 2020, 33, 549–557. [Google Scholar] [CrossRef] [PubMed]
  23. Hagiwara, A.; Hori, M.; Suzuki, M.; Andica, C.; Nakazawa, M.; Tsuruta, K.; Takano, N.; Sato, S.; Hamasaki, N.; Yoshida, M.; et al. Contrast-enhanced synthetic MRI for the detection of brain metastases. Acta Radiol. Open 2016, 5, 2058460115626757. [Google Scholar] [CrossRef] [PubMed]
  24. Blystad, I.; Warntjes, J.B.M.; Smedby, Ö.; Lundberg, P.; Larsson, E.M.; Tisell, A. Quantitative MRI using relaxometry in malignant gliomas detects contrast enhancement in peritumoral oedema. Sci. Rep. 2020, 10, 17986. [Google Scholar] [CrossRef] [PubMed]
  25. Müller, A.; Jurcoane, A.; Kebir, S.; Ditter, P.; Schrader, F.; Herrlinger, U.; Tzaridis, T.; Mädler, B.; Schild, H.H.; Glas, M.; et al. Quantitative T1-mapping detects cloudy-enhancing tumor compartments predicting outcome of patients with glioblastoma. Cancer Med. 2017, 6, 89–99. [Google Scholar] [CrossRef] [Green Version]
  26. Sacks, P.; Rahman, M. Epidemiology of Brain Metastases. Neurosurg. Clin. N. Am. 2020, 31, 481–488. [Google Scholar] [CrossRef]
  27. Kelly, P.J. Gliomas: Survival, origin and early detection. Surg. Neurol. Int. 2010, 1, 96. [Google Scholar] [CrossRef] [Green Version]
  28. Schouten, L.J.; Rutten, J.; Huveneers, H.A.; Twijnstra, A. Incidence of brain metastases in a cohort of patients with carcinoma of the breast, colon, kidney, and lung and melanoma. Cancer 2002, 94, 2698–2705. [Google Scholar] [CrossRef]
  29. Nowosielski, M.; Radbruch, A. The emerging role of advanced neuroimaging techniques for brain metastases. Chin. Clin. Oncol. 2015, 4, 23. [Google Scholar]
  30. Shah, A.D.; Shridhar Konar, A.; Paudyal, R.; Oh, J.H.; LoCastro, E.; Nuñez, D.A.; Swinburne, N.; Vachha, B.; Ulaner, G.A.; Young, R.J.; et al. Diffusion and Perfusion MRI Predicts Response Preceding and Shortly After Radiosurgery to Brain Metastases: A Pilot Study. J. Neuroimaging 2021, 31, 317–323. [Google Scholar] [CrossRef]
  31. Taunk, N.K.; Oh, J.H.; Shukla-Dave, A.; Beal, K.; Vachha, B.; Holodny, A.; Hatzoglou, V. Early posttreatment assessment of MRI perfusion biomarkers can predict long-term response of lung cancer brain metastases to stereotactic radiosurgery. Neuro-Oncology 2018, 20, 567–575. [Google Scholar] [CrossRef] [PubMed]
  32. Kang, K.M.; Choi, S.H.; Hwang, M.; Yoo, R.E.; Yun, T.J.; Kim, J.H.; Sohn, C.H. Application of Synthetic MRI for Direct Measurement of Magnetic Resonance Relaxation Time and Tumor Volume at Multiple Time Points after Contrast Administration: Preliminary Results in Patients with Brain Metastasis. Korean J. Radiol. 2018, 19, 783–791. [Google Scholar] [CrossRef] [PubMed]
  33. Russek, S.; Boss, M.; Jackson, E.; Jennings, D.; Evelhoch, J.; Gunter, J.; Sorensen, A. Characterization of NIST/ISMRM MRI system phantom. In Proceedings of the 20th Annual Meeting of ISMRM, Melbourne, Australia, 5–11 May 2012; p. 2456. [Google Scholar]
  34. Keenan, K.E.; Stupic, K.F.; Boss, M.A.; Russek, S.E.; Chenevert, T.L.; Prasad, P.V.; Reddick, W.E.; Zheng, J.; Hu, P.; Jackson, E.F. Comparison of T1 measurement using ISMRM/NIST system phantom. In Proceedings of the 24th Annual Meeting of ISMRM, Singapore, 7–13 May 2016; p. 3290. [Google Scholar]
  35. Warntjes, J.B.; Leinhard, O.D.; West, J.; Lundberg, P. Rapid magnetic resonance quantification on the brain: Optimization for clinical usage. Magn. Reson. Med. 2008, 60, 320–329. [Google Scholar] [CrossRef] [PubMed]
  36. Eisenhauer, E.A.; Therasse, P.; Bogaerts, J.; Schwartz, L.H.; Sargent, D.; Ford, R.; Dancey, J.; Arbuck, S.; Gwyther, S.; Mooney, M.; et al. New response evaluation criteria in solid tumours: Revised RECIST guideline (version 1.1). Eur. J. Cancer 2009, 45, 228–247. [Google Scholar] [CrossRef] [PubMed]
  37. Cheng, H.L.; Stikov, N.; Ghugre, N.R.; Wright, G.A. Practical medical applications of quantitative MR relaxometry. J. Magn. Reson. Imaging 2012, 36, 805–824. [Google Scholar] [CrossRef]
  38. Callaghan, M.F.; Mohammadi, S.; Weiskopf, N. Synthetic quantitative MRI through relaxometry modelling. NMR Biomed. 2016, 29, 1729–1738. [Google Scholar] [CrossRef] [Green Version]
  39. Bobman, S.; Riederer, S.; Lee, J.; Suddarth, S.; Wang, H.; MacFall, J. Synthesized MR images: Comparison with acquired images. Radiology 1985, 155, 731–738. [Google Scholar] [CrossRef]
  40. Granberg, T.; Uppman, M.; Hashim, F.; Cananau, C.; Nordin, L.; Shams, S.; Berglund, J.; Forslin, Y.; Aspelin, P.; Fredrikson, S. Clinical feasibility of synthetic MRI in multiple sclerosis: A diagnostic and volumetric validation study. Am. J. Neuroradiol. 2016, 37, 1023–1029. [Google Scholar] [CrossRef] [Green Version]
  41. Betts, A.M.; Leach, J.L.; Jones, B.V.; Zhang, B.; Serai, S. Brain imaging with synthetic MR in children: Clinical quality assessment. Neuroradiology 2016, 58, 1017–1026. [Google Scholar] [CrossRef]
  42. Lee, S.M.; Choi, Y.H.; Cheon, J.-E.; Kim, I.-O.; Cho, S.H.; Kim, W.H.; Kim, H.J.; Cho, H.-H.; You, S.-K.; Park, S.-H. Image quality at synthetic brain magnetic resonance imaging in children. Pediatr. Radiol. 2017, 47, 1638–1647. [Google Scholar] [CrossRef]
  43. Cui, Y.; Han, S.; Liu, M.; Wu, P.Y.; Zhang, W.; Zhang, J.; Li, C.; Chen, M. Diagnosis and Grading of Prostate Cancer by Relaxation Maps From Synthetic MRI. J. Magn. Reson. Imaging 2020, 52, 552–564. [Google Scholar] [CrossRef] [PubMed]
  44. Jung, Y.; Gho, S.M.; Back, S.N.; Ha, T.; Kang, D.K.; Kim, T.H. The feasibility of synthetic MRI in breast cancer patients: Comparison of T(2) relaxation time with multiecho spin echo T(2) mapping method. Br. J. Radiol. 2018, 92, 20180479. [Google Scholar] [CrossRef]
  45. Zhao, L.; Liang, M.; Wu, P.-y.; Yang, Y.; Zhang, H.; Zhao, X. A preliminary study of synthetic magnetic resonance imaging in rectal cancer: Imaging quality and preoperative assessment. Insights Imaging 2021, 12, 120. [Google Scholar] [CrossRef] [PubMed]
  46. Meng, T.; He, N.; He, H.; Liu, K.; Ke, L.; Liu, H.; Zhong, L.; Huang, C.; Yang, A.; Zhou, C.; et al. The diagnostic performance of quantitative mapping in breast cancer patients: A preliminary study using synthetic MRI. Cancer Imaging 2020, 20, 88. [Google Scholar] [CrossRef] [PubMed]
  47. Kumar, R.; Delshad, S.; Macey, P.M.; Woo, M.A.; Harper, R.M. Development of T2-relaxation values in regional brain sites during adolescence. Magn. Reson. Imaging 2011, 29, 185–193. [Google Scholar] [CrossRef] [Green Version]
  48. Bansal, R.; Hao, X.; Liu, F.; Xu, D.; Liu, J.; Peterson, B.S. The effects of changing water content, relaxation times, and tissue contrast on tissue segmentation and measures of cortical anatomy in MR images. Magn. Reson. Imaging 2013, 31, 1709–1730. [Google Scholar] [CrossRef] [Green Version]
Cancers 14 02651 g001 550
Figure 1. Representative MRI data from a 58-year-old male BM patient treated with focal radiation therapy (RT), showing normal-appearing brain tissue. (A,B) standard T2w and post-contrast T1w images depicting the anatomical structures, (C,D) T1 and T2 maps were estimated using MAGnetic resonance image Compilation (MAGiC) for a single slice exhibiting gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF).
Figure 1. Representative MRI data from a 58-year-old male BM patient treated with focal radiation therapy (RT), showing normal-appearing brain tissue. (A,B) standard T2w and post-contrast T1w images depicting the anatomical structures, (C,D) T1 and T2 maps were estimated using MAGnetic resonance image Compilation (MAGiC) for a single slice exhibiting gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF).
Cancers 14 02651 g001
Cancers 14 02651 g002 550
Figure 2. Representative MRI data of a 40-year-old male BM patient treated with whole-brain RT using standard and synthetically-generated (MAGiC) images. (A–C) T2w, T2w fluid attenuated inversion recovery (FLAIR) and T1w post-contrast images, respectively, from standard clinical imaging. (D–H) T2w, T1w, phase-sensitive inversion recovery (PSIR), T1w FLAIR, and short tau inversion recovery (STIR) images, respectively, from MAGiC.
Figure 2. Representative MRI data of a 40-year-old male BM patient treated with whole-brain RT using standard and synthetically-generated (MAGiC) images. (A–C) T2w, T2w fluid attenuated inversion recovery (FLAIR) and T1w post-contrast images, respectively, from standard clinical imaging. (D–H) T2w, T1w, phase-sensitive inversion recovery (PSIR), T1w FLAIR, and short tau inversion recovery (STIR) images, respectively, from MAGiC.
Cancers 14 02651 g002
Cancers 14 02651 g003 550
Figure 3. Representative MRI data from the same patient as Figure 2, showing quantitative images. (A and B) synthetically-generated T1 and T2 maps, respectively, from MAGiC.
Figure 3. Representative MRI data from the same patient as Figure 2, showing quantitative images. (A and B) synthetically-generated T1 and T2 maps, respectively, from MAGiC.
Cancers 14 02651 g003
Cancers 14 02651 g004 550
Figure 4. T1 and T2 values for normal-appearing GM, WM, and CSF using T1 and T2 maps generated from MAGiC methods. Boxes represent the interquartile range; whiskers represent the range of all values; the horizontal line within the box is the median value.
Figure 4. T1 and T2 values for normal-appearing GM, WM, and CSF using T1 and T2 maps generated from MAGiC methods. Boxes represent the interquartile range; whiskers represent the range of all values; the horizontal line within the box is the median value.
Cancers 14 02651 g004
Cancers 14 02651 g005 550
Figure 5. The heat map exhibits the T1 and T2 metric values of all the BM lesions and normal-appearing GM and WM.
Figure 5. The heat map exhibits the T1 and T2 metric values of all the BM lesions and normal-appearing GM and WM.
Cancers 14 02651 g005
Table
Table 1. Patient characteristics.
Table 1. Patient characteristics.
CharacteristicsValue
Total patients11
Total number of BM lesions17
Demographics
  Median age (y)52
  Age range (y)25–61
  Male/Female5/6
Location of primary tumor
  Lung5
  Colon1
  Melanoma2
  Other3
Untreated/Treated3/8
Table
Table 2. T1 and T2 values from three different methods.
Table 2. T1 and T2 values from three different methods.
Vial #T1 (in ms)Percent of Difference (in %)
VPGSMAGiCVP and GSVP and MAGiCGS and MAGiC
118381779.717193.26.53.4
213981350.911793.415.712.7
3998.3957.78524.114.711
4725.8678.26226.614.38.3
55094834535.1116.2
6367345.93275.710.95.5
7258.7242.13006.41623.9
Vial #T2 (in ms)Percent of Difference (in %)
VPGSMAGiCVP and GSVP and MAGiCGS and MAGiC
1645.8537.459116.88.510
2423.6357.441415.62.315.8
3286245.9287140.316.7
4184.8162.6186120.614.4
5134.1118.314111.85.119.2
694.481.610313.69.126.2
762.556.7749.318.430.5
VP: Phantom Vendor Provided; GS: Gold Standard; MAGiC: MAGnetic resonance image Compilation.
Table
Table 3. T1 and T2 values in patients from the untreated and treated groups.
Table 3. T1 and T2 values in patients from the untreated and treated groups.
Relaxometry ValuesUntreated BMTreated BM
T1
(ms)		Median (min, max)1845 (1583,2177)2311 (1654, 2558)
Mean ± SD1868 ± 2982211 ± 269
T2
(ms)		Median (min, max)97 (85, 119)104 (92, 154)
Mean ± SD100 ± 17114 ± 20
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Konar, A.S.; Shah, A.D.; Paudyal, R.; Fung, M.; Banerjee, S.; Dave, A.; Hatzoglou, V.; Shukla-Dave, A. Quantitative Synthetic Magnetic Resonance Imaging for Brain Metastases: A Feasibility Study. Cancers 2022, 14, 2651. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers14112651

AMA Style

Konar AS, Shah AD, Paudyal R, Fung M, Banerjee S, Dave A, Hatzoglou V, Shukla-Dave A. Quantitative Synthetic Magnetic Resonance Imaging for Brain Metastases: A Feasibility Study. Cancers. 2022; 14(11):2651. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers14112651

Chicago/Turabian Style

Konar, Amaresha Shridhar, Akash Deelip Shah, Ramesh Paudyal, Maggie Fung, Suchandrima Banerjee, Abhay Dave, Vaios Hatzoglou, and Amita Shukla-Dave. 2022. "Quantitative Synthetic Magnetic Resonance Imaging for Brain Metastases: A Feasibility Study" Cancers 14, no. 11: 2651. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers14112651

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop