1.
Introduction
Multiple sclerosis (MS) is a chronic immune disease of the brain characterized by formation of lesions in the central nervous system and myelin loss [1]. In MS, there is a loss of oligodendrocytes, resulting in thinning or complete loss of myelin. Over the years, the process of remyelination will take place, and the inability of oligodendrocytes to rebuild the myelin sheath will result in the development of scars like lesions or plaque around the damaged axons [1].
Magnetic resonance imaging (MRI) is the preferred imaging modality for diagnosing MS. Visual assessment of these plaques using contrast-enhanced MRI can be challenging due to subjective sensitivity evaluation in early stages. Although MRI can be used to assess the lesion changes in MS, there is a need to use advanced MRI methods to investigate new MS plaques to allow early diagnosis. Automated assessment of the mean lesion volume, myelin-to-white-matter (WM) ratio and WM fraction (WMF) and comparing them to data obtained from a healthy reference population will provide a reliable and reproducible method that is comparable to the visual assessment [2]. Also, the use of deep learning algorithms to automate the quantification of plaque number, along with the mean lesion volume of myelin loss, can improve the sensitivity and accuracy of MS diagnosis [3].
Synthetic MRI NEURO (SyMRI version: 11.3.6; SyntheticMR, Linköping, Sweden) is a new innovation in the market which has transformed the subjective way of reporting in MRI to quantitative reporting. This software helps to measure the accurate tissue properties of the different tissues in the brain. SyMRI works on the multi-delay multi-echo acquisition (MDME) sequence, which acquires the relaxation properties of the patient in an acquisition time of 6 minutes [4].
The basic concept of the sequence is to work on the relaxation properties of the particular tissues, so it also provides relaxation maps with a color coding scheme, and T1, T2, R1, R2 and PD relaxations maps are provided by the software [5]. These relaxation times can be used for early detection of pathologies, or to see infiltration of any pathology to the adjacent normal tissues. Relaxation maps provided by the software constitute the raw data which is post-processed, and synthetic contrast weighted images are formed. A total of 12 different contrast images are provided, which include the T1w, T2w, FLAIR, proton density weighted (PDw), Diffusion inversion recovery (DIR), Phase-sensitive inversion recovery (PSIR), Short-TI Inversion Recovery (STIR) and other contrast images; with the suppression of WM, gray matter (GM) or fat can be generated in less than 10 seconds [4]. The images are synthetic, so parameters like repetition time (TR), echo time (TE) and inversion time (TI) can also be adjusted as per the requirements.
Brain tissue segmentation i.e., myelin, GM, WM and cerebrospinal fluid (CSF) segmentation, is also provided by the software. The myelin segmentation provided by the software is the first of its kind, and it can help better the visualization of the myelin. Also, automatically quantified reports for the segmentation are provided by the software, which gives information about the accurate volume of particular material in the brain, which includes the intracranial volume (ICV), brain parenchymal volume (BPV), brain parenchymal fraction (BPF) [4].
Quantitative values for each slice are also available in the software; in addition to this, if more information is required on a particular area that can also be acquired by freehand drawing the region of interest (ROI) for that particular region. Another unique feature is the reference curves provided by the software, which helps to compare the particular patient with a normal patient of the same age; it also tells if the patient has normal values or if their values are high or low [4].
The software comes with an easy interface and a lot other features, such as distance measurement, contrast enhancement, freehand ROI, automated ROI and windowing tools [4].
SyMRI Neuro is the only clinical USA Food and Drug Administration (FDA)-approved product available in the market which provides myelin segmentation and myelin estimation, and it easily provides us the myelin details.
This study was planned to be conducted in two phases. In Phase 1, we conducted a pilot study to evaluate the feasibility of SyMRI to detect the myelin loss and WMF loss in MS patients with white-matter hyperintensities (WMH) as compared to the non-MS patients with WMHs. In Phase 2, our aim is to analyze the number of plaques on the conventional MRI and synthetic generated T1w FLAIR images, to compare them with the number of plaques generated by artificial intelligence (AI) outputs and to determine how the mean lesion volume and lesion count can benefit the radiologists in day-to-day clinical practice.
2.
Materials and methods
2.1. Study Population
A retrospective case-control study was conducted in a private diagnostic center to quantify the myelin loss in patients with MS. MRI images were acquired from 15 MS patients with WMHs (mean age = 47.9 years; 7 males, 8 females) (herein referred to as the test group), and from 15 non-MS patients with WMHs (mean age = 43.8 Years; 8 males, 7 females) (herein referred to as the control group). The study was conducted in accordance with the Declaration of Helsinki, 1996. Ethics approval was not required because a waiver was obtained due to the use of retrospective data. Data were collected between March 2021 and March 2022. The patients in this study were scanned for internal research purposes. Written informed consent was obtained from all the patients.
2.2. Image Acquisition
All patients underwent MRI scans on a 3T MRI scanner (Discovery MR750w; GE Healthcare; Milwaukee, WI, USA) with a 48-channel head coil. All subjects underwent a fast multi-delay multi-echo acquisition (MDME) sequence called MAGiC, which is a customized version of SyntheticMR's SyMRI® IMAGE software marketed by GE Healthcare under a license agreement. Quantitative SyMRI was performed by applying an MDME sequence using a 2D axial pulse sequence with different combinations of TEs and saturation delay times (TR = 4000 ms, TE1 = 22 ms, TE2 = 89 ms, slice thickness =4 mm, interslice gap = 1 mm, flip angle = 90°, matrix = 320 × 256 and rFOV = 0.8). At our diagnostic center, two TEs and four delay times were used to generate eight real images and eight imaginary ones, which were then used to quantify the longitudinal T1 and transverse T2 relaxation times and the PD. The data acquired at each section were used to produce T1, T2 and PD maps, which were then used to calculate the synthetic magnetic resonance images. The total image acquisition time was 6 minutes. All subjects also underwent conventional diffusion-weighted imaging, T1w FLAIR and T2w imaging. Coronal T2 FLAIR images were obtained using the following parameters: TE = 120.4 ms; TR = 9000 ms; TI = 2248 ms; slice thickness = 5 mm; interslice gap = 5.5 mm.
2.3. Data Analysis
White-matter lesions (WMLs) that appear hyperintense/bright on T1w FLAIR images were included in the analysis.
SyMRI image analysis was done using SyMRI post-processing software (SyMRI version: 11.3.6; SyntheticMR, Linköping, Sweden) that automatically and simultaneously generates different synthetic weighted contrasts, i.e., T1w, T2w, PDw, T1w FLAIR, T2 FLAIR, DIR and PSIR images. Also, the SyMRI imaging data were used to automatically segment the total brain volume, including the WM, GM, CSF volume, total ICV and the myelin (MyC) partial maps, along with the R1 and R2 relaxation maps.
The MyC-to-WM ratio was calculated by using SPSS software to quantify the signal intensities of the normal appearing white matter in the test group and control group, and their mean values were recorded. All measurements were performed in a randomized and blinded manner. Also, age-stratified reference curves for the given population were obtained and compared with the normative value, together with data plotting, curve fitting and confidence interval calculation. For each tissue type, data were plotted against the patient age.
Brain tissue measurements were normalized to ICV, yielding the WM fraction (WMF = WM volume/ICV) and GM fraction (GMF = GM volume/ICV).
We hypothesized that there would be myelin loss in the test group, and that this may be revealed with the quantitative synthetic MRI results for myelin quantification. Hence, we used the quantitative synthetic MRI results to investigate the association between myelin partial maps and WM volume.
2.4. Statistical Analysis
The data were evaluated for normality by using histograms and a Kolmogorov–Smirnov test; the variables are reported as mean ± SD.
Differences in the means of the WM, GM and CSF volumes, myelin partial fraction map, WM-to-myelin ratio and WMF were assessed using a Mann-Whitney U test.
Statistical analysis was performed using SPSS v22 (IBM Released 2013, Armonk, NY, USA: IBM). The significance level for the tests was set at p < 0.05 (two-sided).
3.
Results
3.1. Demographics of the Participants
Table 1 compares the demographic data of the ages, MyC-to-WM fractions, total myelin volumes and total WM volumes between the test group and control group. The mean age of the test group was 47.93 ± 19.44 years, and that of the control group was 43.8 ± 8.03 years. No significant differences in age between the test and control groups was noted (P = 0.819).
3.2. WMF Analysis
Non-parametric independent sample t-tests were used to assess the differences in WMF, WM volume, GMF and ICV between the test group and the control group. The results showed that the WMF was significantly lower in the test group than in the control group (38.8% vs 33.2%, p < 0.001). Interestingly, we noted that the mean WM volume was significantly different between the control and the test group (511.64 ± 91.55 vs 441.6 ± 90.55, p = 0.012). Also, there were no significant differences in GMF and ICV between the test group and the control group (see Table 2).
3.3. Myelin Fraction Analysis
Differences in the means of myelin volume fraction and MyC-to-WM fraction were assessed using the Mann-Whitney U nonparametric t-test. We noted a significant difference in the mean myelin volume between the test group and the control group (138.29 ± 29.28 vs.158.66 ± 32.31, p = 0.044).
3.4. Grey Matter Volume Fraction Analysis
GMF analysis revealed that there were no significant differences in the mean GM volume fraction (i.e., the GMF) or BPV between the test and the control group.
4.
Discussion and conclusion
In this study, we aimed to find the benefit and use of SyMRI software for automatic segmentation of the different types of brain tissue volumes in MS patients. We found that the myelin volume, WM volume and WMF were significantly lower in the test group as compared to the control group. There was no significant difference in the GM volume between both the groups. It is well established in theory that MS primarily affects the WM [1], and we found that to be true in our study as well.
To our knowledge, this is the first study to estimate the WM content in patients with MS by using synthetic MRI in an Indian population. A previous multi-center, multi-reader study has indicated that there is comparable accuracy between synthetic MRI data and conventional MRI on neuroradiology reports [6]. Another study indicated that the contrast ratio is better for synthetic versus conventional DIRs; also, the number of plaques reported for synthetic MRI is larger than that for conventional MRI [3]. Another strength of this study is that we took into account the individual brain sizes in this study by comparing the WMF (with respect to the ICV) and the WM volume. Synthetic MRI has also been found to be faster than conventional MRI, thus saving valuable time for the radiographer and the patient.
However, through this study, we have found that the 6-min MAGiC sequence for automatic segmentation of the brain volume and generation of different contrast images can decrease the workload of the radiologist and the infrastructure requirements for undergoing an MRI scan, as well as increase the efficiency of the system as a whole.
The main limitation of this study is the small sample size. Also, we did not perform any ROI analysis on the scans to look for the exact location of WM loss. Synthetic T2FLAIR is more prone to artifacts than conventional T2FLAIR, which is specifically seen in the posterior limb of the internal capsule [6], or next to the CSF [3]. Few studies have shown that GM might be heavily affected by MS as well [7],[8], but that was not the case with this study.
Although MRI is a sensitive method for detecting the WMLs, we did not assess the correlation between the number of plaques and disease severity.
Though the healthcare industry has been relatively slow in adapting AI into their workflow, deep learning (DL) tools are commercially available for clinical settings. However, adoption of these software packages in everyday clinical practice is limited. For Phase 2 of the study, we are specifically interested in comparing the myelin loss in MS patients to that of non-MS patients, and we would also like to explore the added benefits of using these DL tools for assessment of the quantitative outputs in a clinical setting. These DL tools, along with the SyMRI software, can be useful to quantify the WM, myelin loss and number of plaques in patients with MS. We will also correlate the number of plaques and the amount of WM loss with disease severity or the reproducibility of synthetic MRIs with conventional MRIs.
To conclude, this study shows that SyMRI can be used as an alternative to conventional MRI to generate the different contrasts. Thus, the measurements of WM volume, WMF and myelin volume can be assessed using SyMRI with good reliability and reproducibility.