Tumor grade in soft-tissue sarcoma: Prediction with magnetic resonance imaging texture analysis.

To determine the value of 3T magnetic resonance imaging (MRI) texture analysis in differentiating high- from low-grade soft-tissue sarcoma.Forty-two patients with soft-tissue sarcomas who underwent 3T MRI were analyzed. Qualitative and texture analysis were performed on T1-, T2- and fat-suppressed contrast-enhanced (CE) T1-weighted images. Various features of qualitative and texture analysis were compared between high- and low-grade sarcoma. Areas under the receiver operating characteristic curves (AUC) were calculated for texture features. Multivariate logistic regression analysis was used to analyze the value of qualitative and texture analysis.There were 11 low- and 31 high-grade sarcomas. Among qualitative features, signal intensity on T1-weighted images, tumor margin on T2-weighted images, tumor margin on fat-suppressed CE T1-weighted images and peritumoral enhancement were significantly different between high- and low-grade sarcomas. Among texture features, T2 mean, T1 SD, CE T1 skewness, CE T1 mean, CE T1 difference variance and CE T1 contrast were significantly different between high- and low-grade sarcomas. The AUCs of the above texture features were > 0.7: T2 mean, .710 (95% confidence interval [CI] .543-.876); CE T1 mean, .768 (.590-.947); T1 SD, .730 (.554-.906); CE T1 skewness, .751 (.586-.916); CE T1 difference variance, .721 (.536-.907); and CE T1 contrast, .727 (.530-.924). The multivariate logistic regression model of both qualitative and texture features had numerically higher AUC than those of only qualitative or texture features.Texture analysis at 3T MRI may provide additional diagnostic value to the qualitative MRI imaging features for the differentiation of high- and low-grade sarcomas.

Introduction

Soft-tissue sarcomas are heterogeneous malignant neoplasms, and tumor histological grade, size, margin, stage, and patient age, among others, have been reported to be prognostic factors. Among these, histological grade is one of the most important.[ Histological grade is determined by pathological analysis of surgically obtained specimens and is based on mitotic count, differentiation, and necrosis.[ However, predicting tumor grade before surgery is critical because high-grade sarcomas usually require neoadjuvant chemotherapy.[

MRI is a tool for preoperatively diagnosing and grading soft-tissue sarcomas and determining tumor extent.[ There have been studies investigating the prediction of soft-tissue sarcoma grade using MRI features.[

Texture analysis is one such new technique for analyzing images as a part of radiomics, and can be used to quantitatively analyze tumoral heterogeneity, which arises from various causes, such as cellularity, angiogenesis, extravascular extracellular matrix or areas of necrosis, and has implications for tumor prognosis.[ Using various mathematical methods, texture analysis evaluates the gray-scale intensity and spatial relationship of pixels for assessing tumoral heterogeneity.[ Texture analysis has been studied to predict tumor stage and survival in various cancers,[ differentiate benign from malignant lesions in breast,[ and to predict treatment response in various other cancers.[ To our knowledge, however, there have been a small number of studies [ investigating whether texture analysis based on 3T preoperative MRI can predict the grade of soft-tissue sarcomas.

Thus, we hypothesized that texture analysis, which provides information about tumoral heterogeneity, may be useful in predicting the grade of soft-tissue sarcomas. The purpose of this study was to determine the value of 3T MRI texture analysis for differentiating high-grade from low-grade soft-tissue sarcomas.

Methods

Patient population

The institutional review board approved this retrospective study. It was performed in accordance with procedures complied with HIPAA guidelines. The requirement for written informed consent was waived by the institutional review board due to its retrospective nature. Between June 2010 and December 2017, 108 patients were referred for surveillance of soft-tissue sarcoma at the authors’ institute. Inclusion criteria were as follows: patients who underwent preoperative 3T MRI, and who underwent definite histopathologic confirmation including tumor grade after surgery. Patients treated with preoperative chemotherapy or radiation therapy, had poor image quality, small (< 1 cm) lesion, unavailable contrast-enhance (CE) images, diagnosed only by biopsy, inappropriate pathological findings from the surgical specimen for assessment of tumor grade, and software technical error, were excluded (Fig. 1). Forty-two patients with soft-tissue sarcomas, who underwent 3T MRI, including CE imaging, and were pathologically graded after surgery, were ultimately included in this study.

Figure 1

Flow diagram of the study.

MRI protocols

MR examinations were performed using a 3T MRI unit (MAGNETOM Verio; Siemens Healthineers, Erlangen, Germany). The MRI protocols included longitudinal fat-suppressed T2-weighted turbo spin-echo (TSE) imaging, axial T1-weighted TSE imaging, and axial T2-weighted TSE imaging with and without fat-suppression. After contrast material injection, longitudinal and axial fat-suppressed T1-weighted images were acquired. The acquisition parameters were as follows: field of view, 100 to 280 mm2; matrix size, 512 x 216; slice thickness, 3 to 10 mm; intersection gap, 0 mm; TR, 680 to 870/4000-5600 msec; TE, 11 to 21/63 to 83 msec; turbo fator, 3/13; and number of excitations, 1.

Qualitative analysis of conventional MRI

Conventional MRI was reviewed by a musculoskeletal radiologist (J.H.H, 3 years’ experience of musculoskeletal radiology). The reader was blinded to the histologic grade of the tumors. The following qualitative MRI features were assessed: size, tumor location, signal intensity, signal intensity heterogeneity and tumor margin on T1- and T2-weighted images, peritumoral high signal intensity on T2-weighted images, enhancement appearance, tumor margin on fat suppressed CE T1-weighted images, peritumoral enhancement, cortex extension and marrow extension.

MRI texture analysis

Axial T1- and T2-weighted images and fat-suppressed CE T1-weighted images were obtained from the picture archiving and communication system (PACS) and loaded into prototype software (Multiparametric Analysis, Siemens Healthineers, Erlangen, Germany) for co-registration, lesion segmentation, and texture analysis. A musculoskeletal radiologist (J.H.H, 3 years’ experience in musculoskeletal radiology), who was blinded to clinical information and the pathology report, reviewed the MR images on the PACS and drew regions of interest (ROI) of soft-tissue sarcomas using the post-processing workstation. On multiple slices, an ROI was drawn to include the entire tumor to the maximum extent possible on CE T1-weighted images. T1- and T2-weighted images were spatially registered with fat-suppressed CE T1-weighted images using deformable registration. Thus, ROIs were co-localized on all axial images. To avoid partial-volume effects, tumor margins were not included. After drawing an ROI, the volume of interest (VOI) in the tumor was automatically calculated using the workstation. And then, texture analysis was performed.

In our study, statistical based texture analysis was performed. The statistical based technique generates several subsets of texture features, such as first-order statistics, from the histogram of signal intensities in the segmented VOI and second-order statistics, also known as the gray-level co-occurrence matrix (GLCM). Texture analysis, including first- and second-order statistics, was performed using the Multiparametric Analysis prototype software. Mean intensity, standard deviation (SD), skewness, and kurtosis were derived from the image histogram, and difference entropy, difference variance, contrast, and entropy were derived from the GLCM.

Pathological analysis

One pathologist (C.K.J, 17 years’ experience in musculoskeletal pathology) assessed the tumor grade based on pathology findings from surgical specimens obtained from the institution. The pathologist was blinded to the MRI findings. Tumors were graded according to the French Federation of Comprehensive Cancer Centers system.[ Tumor grade was classified into 1 of 3 groups: grade 1, grade 2, and grade 3.

Statistical analysis

The pathological findings were used as the standard of reference. Tumor grades 2 and 3 were combined as high-grade for statistical analysis. Qualitative MRI features of low- and high-grade sarcomas were compared using Chi-squared test. Texture features were compared between high- and low-grade soft-tissue sarcomas using the Mann–Whitney U test. The areas under the receiver operating characteristic curves (AUC) for texture features were calculated to evaluate the diagnostic performance of each feature. Multivariate logistic regression analysis was conducted using features of qualitative and texture analysis which were significant in univariate analysis to analyze the diagnostic performance to predict tumor grade. The sensitivity and specificity of each multivariate model was calculated based on the cut-off value with highest Youden index. Statistical analysis was performed using SPSS version 19 (IBM Corporation, Armonk, NY); P < .05 was considered to be statistically significant.

Results

Patients

Forty-two patients (mean age, 59.1 years [range, 17–95 years]) were included in this study. There were 17 men (mean age, 56.3 years [range, 17–95 years]) and 25 women (mean age, 61 years [range, 23–84 years]). Among the 42 soft-tissue sarcomas, 11 lesions were grade 1, 8 were grade 2, and 23 were grade 3 on histopathology. The mean age of each group was 65.5 ± 11.2 and 56.8 ± 21.1 years for low- and high-grade sarcoma, respectively (P value = .206). The female proportion was 81.8% and 51.6% for low and high grade sarcoma, respectively (P value = .151). Table 1 shows the histopathologic types and grades of the included soft tissue sarcomas. The histological types of the included soft-tissue sarcomas included: undifferentiated pleomorphic sarcoma (n = 6), liposarcoma (n = 16, 9 low grade), myxofibrosarcoma (n = 5, 1 low grade), leiomyosarcoma (n = 3), malignant peripheral nerve sheath tumor (n = 3), synovial sarcoma (n = 2), rhabdomyosarcoma (n = 2), Ewing sarcoma (n = 1), angiosarcoma (n = 1), fibromyxoid sarcoma (n = 1, low grade), clear cell sarcoma (n = 1), and fibrosarcoma (n = 1).

Table 1

Histopathologic types and grades of soft tissue sarcomas.

The mean interval from MRI to surgical resection was 9.2 days (range, 1–37 days).

Qualitative MRI features are summarized in Table 2. High-grade sarcomas were more likely to have poorly-defined margin on T2- (23/31 [74.2%] vs 3/11 [27.3%], P = .011) and fat-suppressed CE T1-weighted images (25/31 [80.6%] vs 3/11 [27.3%], P = .002). Peritumoral enhancement was more common in high- than low-grade sarcoma (25/31 [80.6%] vs 4/11 [36.4%], P = .019). Lesion hyperintensity on T1-weighted images was more frequently encountered in low- than high-grade sarcoma (8/11 [72.7%] vs 10/31 [32.3%], P = .033).

Table 2

Qualitative imaging features of high- and low-grade soft tissue sarcomas.

The sensitivity and specificity of qualitative features which showed significant difference between high- and low-grade sarcomas were as follows: 67.7% and 72.7% for hyperintensity on T1-weighted images; 74.2% and 72.7% for poorly-defined margin on T2-weighted images; 80.6% and 72.7% for poorly-defined margin on CE T1-weighted images; and 80.6% and 63.6% for peritumoral enhancement. The sensitivity and specificity of poorly-defined margin on CE T1-weighted images was the highest.

Texture features based on histogram and GLCM are summarized in Table 3. Among the texture features based on histogram, T2 mean, CE T1 mean, T1 SD, and CE T1 skewness revealed significant differences: T2 mean, T1 SD and CE T1 skewness were significantly lower in high-grade than in low-grade sarcomas (655.0 ± 235.3 vs 850.0 ± 258.9, P = .041; 93.2 ± 65.6 vs 142.9 ± 66.3, P = .025; and .058 ± .679 vs .631 ± .570, P = .014, respectively). On the other hand, CE T1 mean values were significantly higher in high-grade than in low-grade sarcomas (649.2 ± 312.2 vs 385.9 ± 220.1, P = .009). In texture features based on GLCM, CE T1 difference variance and CE T1 contrast demonstrated significant differences, which were higher in high-grade than in low-grade sarcomas (.239 ± .086 vs .175 ± .064, P = .031; .430 ± .130 vs .354 ± .162, P = .027). Differences in other features did not reach statistical significance. Figures 2–4 show MRI findings and these texture features of soft tissue sarcoma in grade 1, 2, and 3.

Table 3

Comparison of texture features of T1, T2, and contrast-enhanced T1 images between low-grade and high-grade soft tissue sarcomas.

Figure 2

A 50-year-old woman with myxoid liposarcoma, grade 1. A mass with a lobulated margin in the medial side of knee is hypointense on axial T1-weighted imaging (a), heterogeneously hyperintense on axial T2-weighted imaging (b) and exhibits heterogeneous enhancement on axial fat-suppressed contrast-enhanced T1-weighted imaging (c). Whole-tumor T2 (upper) and contrast-enhanced (CE) T1 (lower) histograms (d) reveal negative T2 skewness (–0.104) and high CE T1 skewness (1.345). Texture features were high T1 standard deviation (108) and low CE T1 difference variance (0.176). These findings suggest low-grade sarcoma. (e) Grossly, the tumor is well circumscribed and shows gelatinous, tan-yellow cut surface. (f) In this microscopic image (hematoxylin-eosin stain; original magnification, x100), the tumor of low cellularity is composed of small spindled or ovoid non-lipogenic tumor cells with scant cytoplasm and lipoblasts in background of abundant myxoid stroma and arborizing capillary vasculature.

Figure 4

A 61-year-old woman with undifferentiated sarcoma, grade 3. Axial T1-weighted imaging (a) demonstrates a homogeneously hypointense mass within the posterior thigh. Axial T2-weighted imaging (b) shows relatively homogeneous hypointensity to intermediate signal. The mass shows homogeneous, intense contrast enhancement on axial fat-suppressed contrast-enhanced T1-weighted imaging (c). Whole-tumor T2 (upper) and contrast-enhanced (CE) T1 (lower) histograms (d) reveal high T2 skewness (1.233) and negative CE T1 skewness (-0.529). Texture features are high T2 mean (855.3), high T1 standard deviation (234), and high CE T1 difference variance (0.305). Texture features, except for T2 mean and T1 standard deviation, are compatible with high-grade sarcoma. (e) Grossly, an ill-defined tumor shows a lobulated, firm, tan-white cut surface. (f) In this microscopic image (hematoxylin-eosin stain; original magnification, x400), pleomorphic tumor cells are admixed with abundant chronic inflammatory cells and show frequent mitosis (arrows).

A 63-year-old man with leiomyosarcoma, grade 2. An intramuscular mass with circumscribed margin in the upper arm exhibits intermediate signal intensity on axial T1-weighted imaging (a), slight hyperintensity on axial T2-weighted imaging (b), and heterogeneous enhancement on axial fat-suppressed contrast-enhanced T1-weighted imaging (c). Whole-tumor T2 (upper) and contrast-enhanced (CE) T1 (lower) histograms (d) reveal high T2 skewness (0.402) and low CE T1 skewness (0.382). Texture features are low T1 standard deviation (49.2) and high CE T1 difference variance (0.275). These findings suggest high-grade sarcoma. (e) Grossly, the tumor is well circumscribed and shows firm to fleshy, whitish cut surface. (f) In this microscopic image (hematoxylin-eosin stain; original magnification, x400), spindle tumor cells are arranged in a fascicular architecture and show eosinophilic cytoplasm, elongated blunt-ended nuclei, and frequent mitosis (arrows).

We obtained receiver operating characteristic curves for texture features which demonstrated a significant difference between high- and low-grade sarcomas (Fig. 5). The diagnostic performance for assessment of tumor grade in soft-tissue sarcomas is summarized in Table 4. AUC of the above texture features that demonstrated significant differences were all > 0.7 (.710 [95% CI .543, .876] in T2 mean; .768 [95% CI .590, .947] in CE T1mean; .730 [95% CI .554, .906] in T1 SD; .751 [95% CI .586, .916] in CE T1 skewness; .721 [95% CI .536, .907] in CE T1 difference variance; .727 [95% CI .530, .924] in CE T1 contrast). In addition, we performed multivariate analysis for each texture features to investigate the impact of age and sex on their diagnostic performance. CE T1 mean, CE T1 skewness and CE T1 difference variance remained significant even after the adjustment for age and sex (Table 4).

Figure 5

Receiver operating characteristic curves for texture features, which demonstrate a significant difference between high- and low-grade sarcomas.

Table 4

Diagnostic performance of texture features for assessment of tumor grade in soft tissue sarcomas.

Multivariate models using qualitative and texture analysis

We constructed multivariate logistic regression models using features of qualitative and texture analysis which were significant in univariate analysis. The diagnostic performances of the models are summarized in Table 5. The multivariate model composed of only the texture features had numerically higher AUC than that composed of only the qualitative features. The model composed of both qualitative and texture features had even higher AUC than those models composed of only qualitative or texture features, although the difference was not statistically significant (Fig. 6). The multivariate model composed of both qualitative and texture features had higher sensitivity and specificity than the model composed of only qualitative features.

Table 5

Multivariate logistic regression models using qualitative and/or texture analysis.

Figure 6

Receiver operating characteristic curves for models using features of qualitative and/or texture analysis which were significant in univariate analysis.

Discussion

The results of our study demonstrated that texture analysis at 3T MRI may provide additional diagnostic value to the qualitative MRI features for the differentiation of high- and low-grade sarcomas. Among all texture features, T2 mean, CE T1 mean, T1 SD, CE T1 skewness, CE T1 difference variance, and CE T1 contrast showed significant difference between high- and low-grade sarcomas. The AUCs of these features were fairly high (> 0.7). In multivariate logistic regression analysis, the AUC of both texture and qualitative features was numerically higher than that of either qualitative or texture features alone.

The histological grade of soft-tissue sarcoma is one of the most important prognostic factors. Preoperative prediction of tumor grade is essential for planning a patient's treatment strategy.[ Preoperative biopsy can be used for diagnosing and grading of soft-tissue sarcomas. However, biopsy specimens cannot represent the entire tumor; grading of sarcoma using biopsy, therefore, may be limited. Yang et al[ found that preoperative tumor grades that were evaluated using biopsy were changed after surgical resection in 5 of 66 soft-tissue sarcomas. Texture analysis of whole tumor volume is a non-invasive tool for preoperatively predicting tumor grade.

Fernebro et al and Zhao et al[ revealed that MRI features may be used to differentiate high-grade from low-grade sarcomas. Zhao et al[ reported that high-grade tumors exhibited larger size, poorly defined tumor margin, more heterogeneous signal intensity, and peritumoral high signal intensity on T2-weighted images and peritumoral enhancement on T1-weighted CE images. Previous studies[ have reported that tumoral heterogeneity, tumor margin, and peripheral growth pattern are important factors in predicting tumor grade. We also found that high-grade tumors are more likely to have poorly-defined margin on T2- and fat-suppressed CE T1-weighted images and peritumoral enhancement. These results are consistent with previous studies.[ Low-grade sarcomas showed more frequently hyperintense signal intensity on T1-weighted images in this study. We suspect that this result is due to higher number of liposarcoma among low-grade sarcomas. Unlike previous studies, tumoral heterogeneity on T1-, T2-, and fat-suppressed CE T1-weighted images were not significantly different between high- and low-grade sarcoma. We assumed that the assessment of tumoral heterogeneity by visual and qualitative analysis might be limited. In our study, texture analysis was also used to perform quantitative analysis regarding tumoral heterogeneity and spatial arrangement. We found that texture features, which provide quantitative information about tumoral heterogeneity and growth pattern, may help to differentiate high-grade from low-grade soft-tissue sarcomas.

First-order features are calculated from the histogram of pixel intensity values and reflect tumor heterogeneity; SD, kurtosis, and skewness are such first-order features. Several studies have reported that these first-order features were useful in predicting tumor stage and in assessing treatment response in tumors.[ In our study, T2 mean, CE T1 mean, T1 SD and CE T1 skewness were significantly different between low-grade and high-grade sarcoma. CE T1 skewness was significantly lower in high-grade than in low-grade tumor. This was consistent previous studies[ that have used non-CE and CE computed tomography and MR imaging for other organs. Skewness means the asymmetry of the distribution. Lower skewness means more symmetric distribution. Features derived from histogram analysis reflect cellular distribution of tumor as well as distribution of contrast agent.[ It means that first order features are associated with tumor vascular permeability and then tumors with high vascular permeability may show less heterogeneity in first order texture analysis.[

Second-order features calculated by GLCMs may reflect the relationship between neighboring pixels. Such features have been used to differentiate benign from malignant breast lesions, and to stage tumors in the kidney and rectum.[ In our study, difference entropy, difference variance, contrast and entropy were used as second-order features. Among these features, CE T1 difference variance and CE T1 contrast exhibited significant differences, which were higher in high-grade than low-grade sarcomas. Contrast means degree of difference in the neighboring values, and difference variance describes the dispersion (with regard to the mean) of the difference gray-level distribution of the image. Our current study showed high-grade sarcoma have dissimilar and wider distribution in GLCM. These results suggest high-grade sarcoma are more spatially heterogeneous on CE T1-weighted images.

There are small number of recent studies which have predicted grade of soft-tissue sarcoma using texture analysis.[ They performed texture analysis based on only ADC maps or fat suppressed T2-weighted images or T1- and T2-weighted images. We used T1- and T2-weighted images and fat-suppressed CE T1-weighted images for texture analysis. Among total 24 texture features, CE T1 mean and CE T1 skewness showed the highest discriminatory power to differentiate high-grade from low-grade sarcoma. Previous researches revealed that contrast-enhanced CT or MR texture analysis was useful to differentiate benign from malignant tumor, stage the tumor and predict the survival.[ Our results also suggest that texture analysis based on CE T1 weighted images could be useful tool for prediction of tumor grade in soft tissue sarcoma.

To compare the overall diagnostic performance of qualitative and texture features, we constructed multivariate logistic models composed of qualitative, texture features and both. Although, the difference was not statistically significant, the multivariate model of texture features had higher AUC than that of qualitative features. Furthermore, the addition of texture features to the multivariate model of qualitative features increased both sensitivity and specificity.

Our study had several limitations, the first of which were its retrospective design, lack of a standardized MRI protocol, and relatively small sample size. Although selection bias was possible, we recruited consecutive patients who fulfilled the inclusion criteria. We performed texture analysis using dataset with different voxel size due to its retrospective nature. We used only some of the many variable texture features and, therefore, further studies using more diverse texture features may be needed.

In conclusion, our study demonstrated that several texture features are significantly different between high- and low-grade sarcomas and diagnostic performance of these texture features were fairly high. The addition of texture analysis to the conventional qualitative MRI analysis may provide additional diagnostic value for the prediction of soft tissue sarcoma grade.

Acknowledgments

The authors thank Siemens Healthineers (Robert Grimm, Yohan Son, and Mun Young Paek) for providing the Multiparametric Analysis prototype software.

Author contributions

Acquisition of data: Ji Hyun Hong, Won-Hee Jee, Chan-Kwon Jung, Yang-Guk Chung.

Analysis and interpretation of data: Ji Hyun Hong, Won-Hee Jee.

Critical review: Won-Hee Jee.

Drafting of manuscript: Ji Hyun Hong, Won-Hee Jee, Chan-Kwon Jung, Yang-Guk Chung.

Study conception and design: Ji Hyun Hong, Won-Hee Jee.