Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 13 January 2023
Sec. Cancer Imaging and Image-directed Interventions
This article is part of the Research Topic Methods in Cancer Imaging and Image-directed Interventions View all 17 articles

Magnetic resonance imaging radiomics to differentiate ovarian sex cord-stromal tumors and primary epithelial ovarian cancers

Meiying ChengMeiying Cheng1Shifang TanShifang Tan1Tian RenTian Ren2Zitao ZhuZitao Zhu3Kaiyu WangKaiyu Wang4Lingjie ZhangLingjie Zhang1Lingsong MengLingsong Meng1Xuhong YangXuhong Yang5Teng PanTeng Pan6Zhexuan YangZhexuan Yang1Xin Zhao*Xin Zhao1*
  • 1Department of Radiology, Third Affiliated Hospital of Zhengzhou University, Zhengzhou, Henan, China
  • 2Department of Information, Third Affiliated Hospital of Zhengzhou University, Zhengzhou, Henan, China
  • 3Medical College, Wuhan University, Wuhan, China
  • 4Magnetic resonance imaging (MRI) Research, GE Healthcare (China), Beijing, China
  • 5Department of Research, Huiying Medical Technology Co., Ltd., Beijing, China
  • 6Department of Research, Third Affiliated Hospital of Zhengzhou University, Zhengzhou, Beijing, China

Objective: To evaluate the diagnostic ability of magnetic resonance imaging (MRI) based radiomics and traditional characteristics to differentiate between Ovarian sex cord-stromal tumors (SCSTs) and epithelial ovarian cancers (EOCs).

Methods: We consecutively included a total of 148 patients with 173 tumors (81 SCSTs in 73 patients and 92 EOCs in 75 patients), who were randomly divided into development and testing cohorts at a ratio of 8:2. Radiomics features were extracted from each tumor, 5-fold cross-validation was conducted for the selection of stable features based on development cohort, and we built radiomics model based on these selected features. Univariate and multivariate analyses were used to identify the independent predictors in clinical features and conventional MR parameters for differentiating SCSTs and EOCs. And nomogram was used to visualized the ultimately predictive models. All models were constructed based on the logistic regression (LR) classifier. The performance of each model was evaluated by the receiver operating characteristic (ROC) curve. Calibration and decision curves analysis (DCA) were used to evaluate the performance of models.

Results: The final radiomics model was constructed by nine radiomics features, which exhibited superior predictive ability with AUCs of 0.915 (95%CI: 0.869-0.962) and 0.867 (95%CI: 0.732-1.000) in the development and testing cohorts, respectively. The mixed model which combining the radiomics signatures and traditional parameters achieved the best performance, with AUCs of 0.934 (95%CI: 0.892-0.976) and 0.875 (95%CI: 0.743-1.000) in the development and testing cohorts, respectively.

Conclusion: We believe that the radiomics approach could be a more objective and accurate way to distinguish between SCSTs and EOCs, and the mixed model developed in our study could provide a comprehensive, effective method for clinicians to develop an appropriate management strategy.

Introduction

Ovarian sex cord-stromal tumors (SCSTs) are rare nonepithelial neoplasms, and represent about 7% of all primary ovarian tumors (1, 2). According to the classification of ovarian tumors by World Health Organization (WHO) (2020), SCSTs are divided into the following three clinicopathologic subcategories: pure stromal tumors, pure sex cord tumors, and mixed sex cord-stromal tumors (3). Fibromas, thecomas, and granulosa cell tumors account for the majority of SCSTs.

Morphologically, SCSTs usually present as solid masses (4), resembling malignant tumors. Although some of the SCSTs, such as fibromas, have a few specific characteristics on conventional imaging features, diagnostic dilemmas may often arise if the tumor shows increased cellularity, or due to necrosis, hemorrhage, edema, or cystic degeneration (2, 4). Other types are even more confusing on account of morphological complexity. Moreover, the rarity of SCSTs contributes to a low degree of suspicion, which makes it susceptible to misdiagnosis as the more common epithelial ovarian cancer (EOC). Clinically, SCSTs most commonly present at early stages (I) and are primarily surgically treated with an overall favorable prognosis (2, 5), while epithelial tumors usually present at advanced stages (III or IV) and are treated with chemotherapy and surgical debulking (2). Hence, a more accurate and objective assessment method to identify SCSTs from ovarian cancers is imperative.

Magnetic resonance imaging (MRI) has been widely used to detect and evaluate adnexal lesions, especially for the indeterminate adnexal masses on ultrasonography (6, 7). The high soft-tissue resolution and ability to characterize the composition of different fluid types allow it to characterize the lesion types more accurately (7). However, all the evaluation by MRI requires the subjective interpretation of radiologists, in addition to the complexity and overlap in imaging characteristics of different diseases, it is still a challenge to differentiate the types of tumors with MRI alone (6).

Radiomics is an emerging method for postprocessing any type of medical image and generating new quantification metrics which have proven to provide important insights into tumor biology, shifting radiology from the traditional visual analysis to a more objective and automated analysis (6). Recently, radiomics studies on adnexal tumors have demonstrated some encouraging advances, including clinical outcomes in ovarian cancer (6, 8, 9), tumor category (9, 10), and subtype differentiation (9, 11). However, there are few studies involving SCSTs. In this study, we developed an MR image-based radiomics diagnostic model to prove that SCSTs and EOCs are separable.

Patients and methods

Patient population

This retrospective pilot study was approved by our Institutional Review Board with a waiver of informed consent. From January 2017 to December 2021, we retrospectively retrieved all MRI pelvic examinations referring to ovarian or adnexal lesions from our Institutional Picture Archiving and Communication System (PACS) and obtained 1867 results. Every report was analyzed by the researchers. Then we identified patients who met the following inclusion criteria (1): pathologically confirmed as SCSTs or EOCs; (2) MRI was performed within 1 month before pelvic or laparoscopic surgery at our institution; (3) no chemotherapy, radiotherapy, or previous gynecological operation prior to MRI examination. The exclusion criteria were: (1) poor MR image quality, so that the focus cannot be clearly observed or delineated; (2) multiple nodules that the primary focus cannot be identified due to a mutual fusion and extensive adhesion pattern or any other reason. Finally, we consecutively reviewed a total of 148 patients with 173 tumors (81 SCSTs in 73 patients and 92 EOCs in 75 patients) as the primary cohort. Details of histopathology are presented in Table 1. Then the primary dataset was randomly split into the development and testing cohorts with a fixed ratio of 8:2 in each category, resulting in 137 tumors (64 SCSTs and 73 EOCs) for the development cohort and 36 tumors (17 SCSTs and 19 EOCs) for testing cohort. A flowchart of the patient selection process is shown in Figure 1.

TABLE 1
www.frontiersin.org

Table 1 Histopathological types of the selected samples.

FIGURE 1
www.frontiersin.org

Figure 1 Flowchart for the patient selection process in this study.

Clinical characteristics such as patient age, menstrual status, endocrine level, cancer antigen 125 (CA125), and risk of ovarian malignancy algorithm (ROMA), were obtained from the hospital information system.

MR image acquisition and interpretation

MR examinations were performed on the 3.0 T system (SIGNA Pioneer, GE Healthcare, and Skyra, Siemens Healthcare) with the phased-array abdominal coil. The conventional MR sequences included T1 weighted imaging (T1WI) in the axial plane, T2 weighted imaging (T2WI) in the axial and sagittal plane, fat-suppressed T2WI (FS-T2WI) in the axial and coronal plane, diffusion-weighted imaging (DWI) with the b value of 1000 s/mm2 in the axial plane, and multiphase contrast-enhanced fat-suppressed T1WI in the axial, sagittal, and coronal plane. Detailed information about the acquisition parameters is presented in Table S1.

Two radiologists (reader 1 and reader 2, with 7 and 15 years of experience in MRI, respectively) independently recorded the conventional imaging features while blinded to the histological results. The recorded features include (1) Maximum diameter (MD) (measured at the lesion slice with the maximum diameter of the tumor in the three-dimensional measurements); (2) Visibility of hemorrhagic component; (3) Solid and cystic components (predominantly cystic, cystic-solid, and predominantly solid corresponding to less than 1/3, 1/3-2/3, and more than 2/3 solid component, respectively); (4) Signal intensity (SI) of the solid components on T2WI (Hypo-, iso-, or hyperintense was relative to the external myometrium; A few purely cystic lesions were not recorded); (5) Apparent diffusion coefficient (ADC) value (measured manually on the DWI derived ADC maps; regions of interest were placed at target areas of the tumor, and areas such as necrosis, hemorrhage, vascular structures were avoided as much as possible; three measurements were obtained and averaged). Discrepancies were resolved by a consensus, or a third radiologist (reader 3 with more than 20 years of experience in gynecologic imaging) would serve as an arbitrator.

Tumor segmentation

In order to reduce the discrepancies related to various scanning parameters and eliminate the internal dependence of image radiomics features on voxel size, all images were set to a fixed voxel size (1mm×1mm×1mm) with linear interpolation algorithm.

The 3D segmentation was performed on open-source software (ITK-SNAP version 3.8.0, http://www.itksnap.org). The volume of interest (VOI) for each lesion on each slice was manually delineated on FS-T2WI by a reader 1. Then, reader 3 confirmed all the regions. To examine the reproducibility of extracted radiomics features and to obtain more robust radiomics features, the intra-class correlation coefficient (ICC) was used to assess the intra- and inter-observer agreement of VOI delineation. So, four weeks later, 30 patients (15 SCSTs and 15 EOCs) were randomly selected. Reader 1 then re-delineated the VOI for intra-consistency testing, in the meantime, reader 2 outlined the VOI according to the same procedure to test inter-consistency. The feature with an ICC > 0.75 was selected for further analysis (12). ICC can be obtained from the following equation:

ICC=(MSRMSE)MSR+(MSCMSEn),

where MSR is mean square for rows, MSC is the mean square for columns, MSE is mean square for error and n represents the number of subjects.

Feature extraction

The extraction of radiomics features was conducted in the Radcloud software (Huiying Medical Technology Co., Ltd, Beijing, China). A total of 1,409 features were extracted from the FS-T2WI sequences using the pyradiomics function package (https://pyradiomics.readthedocs.io/). These features could be divided into the following three categories. First, first-order statistical features, such as peak value, mean, and variance, quantitatively describe the voxel intensity distribution of the lesion area in MR images through common basic indicators. Second, two-dimensional morphological features, describe the two-dimensional shape and size of the lesion. Third, texture features, such as Gray Level Co-occurrence Matrix (GLCM), Gray Level Run Length Matrix (GLCM), and Gray Level Size Zone Matrix (GLSZM), quantify the heterogeneity of the lesion texture. In addition, several filters, such as exponential, logarithm, square, square root, and wavelet (including wavelet-LHL, wavelet-LHH, wavelet-HLL, wavelet-LLH, wavelet-HLH, wavelet-HHH, wavelet-HHL, and wavelet-LLL) filters, were applied to calculate the first-order statistical features and texture features of the transformed image.

Feature selection

Before the selection of radiomics features, normalization processing was performed for all extracted features, and the features were normalized to the normal distribution by mean and variance scaling.

We run a 5-fold cross validation to select features, where each fold we did feature selection and model building on 80% of the development data (the training cohort), and evaluated on the remaining 20% (the validation cohort). In each fold, we implement the dimensionality reduction process in three steps. Features with a variance value of >0.8 were first selected. Then SelectKBest was applied to select the features with a p-value less than 0.05. Finally, the Lassolars algorithm was used to screen the optimal radiomics features, among which the top 10 features with the best contribution would be selected. The least absolute shrinkage and selection operator (LASSO) is a regression analysis method that can perform both variable selection and regularization to improve the identification accuracy and interpretability of the model. Regression algorithm least angle regression (LARS) provides variables through the linear combination of high-dimensional data. It is related to positive stepwise regression. Lassolars algorithm used in this paper is a combination of Lars algorithm and lasso model, which can automatically select the optimal parameters λ and have better performance than LASSO algorithm alone (13). Finally, we counted the frequency of each feature selected in the 5 folds, the features that get selected three or more times repeatedly were considered stable features, and different models were built using various combinations of those selected features.

Classifier modeling

Univariate analysis and multivariate logistic regression analysis were successively performed based on development data to screen the clinical and conventional MR features. Based on the selected radiomics features, clinical features, and conventional MR parameters, five prediction models (clinical model, conventional MR model, traditional model, radiomics model and mixed model) were constructed by using the LR classifier. LR classifier includes classification, function establishment, solving optimal model parameters through optimization iteration, and verifying the model performance. In addition, the nomogram was also constructed to visualize the results of the logistic regression. The nomogram develops scoring criteria based on the magnitude of the regression coefficients for all independent variables. It scores each value level of each independent variable, giving an overall score and finally calculating the probability of disease risk for each patient through a conversion function between the score and the outcome probability. The calibration curve was drawn to assess the agreement between the predicted results and actual presence. The detailed process of radiomics analysis is presented in Figure 2.

FIGURE 2
www.frontiersin.org

Figure 2 Process of radiomics analysis in this study.

The performance of the model was estimated by using receiver operating characteristic (ROC) curves and confusion matrix analysis, and the area under the curve (AUC), sensitivity, specificity, and accuracy were calculated. The DeLong test was used to compare the performance of different models. The decision curve analysis (DCA) was used to evaluate the clinical utility of the prediction models by quantifying the net benefits at different threshold probabilities in the dataset.

Statistical analysis

Statistical analysis was performed using SPSS version 25.0 (IBM). Quantitative variables are shown as mean ± standard deviation. Categorical variables were assessed by Chi-square tests or Fisher exact test, and differences in continuous variables were assessed by t-test or Mann-Whitney U test. P< 0.05 was considered statistically significant. Model construction was executed using R software 3.5.3 (https://cran.r-project.org/) and Python 2.7 software (https://www.python.org/). The packages of “pyradiomics” (https://pyradiomics.readthedocs.io/), “scikitlearn” (https://scikit-learn.org/), and “matplotlib” (https://matplotlib.org/) were used for feature selection, model building, and plotting in this study.

Results

Patients and tumor characteristics

The comparisons of the clinical data and MR parameters between SCSTs and EOCs groups in the primary and development cohorts are summarized in Table 2. Age, ROMA index, serum CA125, ADC value, MD, SI on T2WI, solid and cystic components showed significant differences between the SCSTs and EOCs groups, while no significant difference was observed in menstrual status, endocrine level, and hemorrhagic component. Among the variables with significant differences, only ROMA index (P< 0.001), ADC value (P = 0.004), solid and cystic components (P = 0.043) remained as independent predictors on the multivariate logistic regression analysis.

TABLE 2
www.frontiersin.org

Table 2 The clinical data and MR parameters of the primary and development cohorts.

Feature extraction and selection

Of all the extracted radiomics features, the median ICC was 0.90, and 1277 of 1409 features (91%) were robust and were selected for subsequent analysis, with ICC > 0.75. Lassolars algorithms on feature selection for each fold are shown in Figure S1. The selected radiomic features and the corresponding coefficients in each fold are shown in Figure 3 and Table S2. We did model building in every fold, the AUCs for the 5-fold cross validation are reported in Table S3 and Figure S2, and the mean AUCs were 0.883 ± 0.018 in the training cohort and 0.849 ± 0.021 in the validation cohort. Finally, there were three features got selected in 3 of the 5 folds, three features got selected in 4 of the 5 folds, and three features got selected in every fold. The frequencies of the radiomics features are summarized in Table 3.

FIGURE 3
www.frontiersin.org

Figure 3 The selected radiomic features in each fold. (A) 6 radiomic features selected in fold 1. (B) 9 radiomic features selected in fold 2. (C) 10 radiomic features selected in fold 3. (D) 10 radiomic features selected in fold 4. (E) 10 radiomic features selected in fold 5.

TABLE 3
www.frontiersin.org

Table 3 Details of the selected features in the developed radiomics models.

Construction and performance of the prediction models

Based on the selected clinical variable ROMA index, a clinical model was established. Based on parameter ADC, solid and cystic components, a conventional MR model was established. Then, based on the combination of the above clinical factors and conventional MR parameters, a traditional model was established. The AUCs of the clinical model, conventional MR model, and traditional model were 0.729 (95%CI: 0.645-0.812), 0.737 (95%CI: 0.654-0.819), 0.755 (95%CI: 0.677-0.834) in the development cohort, respectively, and 0.680 (95%CI: 0.498-0.862), 0.693 (95%CI: 0.515-0.872), 0.735 (95%CI: 0.569-0.902) in the testing cohort, respectively (Figure 4 and Table 4).

FIGURE 4
www.frontiersin.org

Figure 4 ROC curves of conventional MR model, clinical model, traditional model, radiomics model and mixed model within the development (A) and testing (B) cohorts.

TABLE 4
www.frontiersin.org

Table 4 Predictive performances of the final models in the development and testing cohorts.

Three different radiomics models (Models 1-3) were developed using the following combinations of the radiomics features: 3 features (got selected in every fold), 6 features (got selected 4 or more times repeatedly) and 9 features (got selected 3 or more times repeatedly). Among them, Model 3 showed the best performance, which was determined as the final radiomics model, with the AUCs of 0.915 (95% CI: 0.869-0.962) and 0.867 (95% CI: 0.732-1.000) in the development and testing cohorts, respectively (Figure 4 and Table 4). However, no significant differences were observed between Model 3 and the other two models. The AUCs of the three models and the comparisons in terms of diagnostic performance among them were shown in Table 5. Rad-score based on 9 features were weighted by their respective coefficients, the calculation formula for the Rad-score is provided in the Supplementary Materials page 7.

TABLE 5
www.frontiersin.org

Table 5 Predictive performances of the three developed radiomics models.

Finally, we established a mixed model based on the Rad-score, clinical characteristics (ROMA), and conventional MR parameters (ADC, solid and cystic components). The AUCs of the mixed model were 0.934 (95% CI: 0.892-0.976) and 0.875 (95% CI: 0.743-1.000) in the development and testing cohorts, respectively (Figure 4 and Table 4).

The difference between the traditional model and radiomics model was statistically significant in the testing cohort (traditional model AUC vs. radiomics model AUC: 0.735 vs. 0.867, Delong test P< 0.001); however, there was no evidence of a difference in the testing set between the radiomics model compared with the mixed model (radiomics model AUC vs. mixed model AUC: 0.867 vs. 0.875, Delong test P = 0.561).

A radiomics nomogram was constructed by using the selected variables from multivariate logistic regression and Rad-score to provide a visualized outcome measure (Figure 5A). The total score for this nomogram was calculated using the formula: Nomo-score = -1.5402 + 11.4118 × Rad-score -0.0335 × ROMA + 1.3928 × ADC + 0.3532 × Components. The calibration curves demonstrated good diagnostic consistency between the predictions of the radiomics nomogram and the actual observations of the samples (Figure 5B).

FIGURE 5
www.frontiersin.org

Figure 5 (A) Radiomics nomogram with Rad-score, two conventional MR findings and a clinical factor, including components, ADC, and ROMA. (B) Calibration curves of radiomics nomogram. The diagonal line represented the perfect prediction of the radiomics nomogram. The black solid line represented the calibration curve of nomogram in the testing cohorts. The calibration curves were close to the diagonal line, which indicated good prediction performance of the nomogram.

DCA revealed that the radiomics model and the mixed model provided a better net benefit than the traditional model across the majority of the range of reasonable threshold probabilities (Figure 6).

FIGURE 6
www.frontiersin.org

Figure 6 Decision curve analysis (DCA) for the three models within the development (A) and testing (B) cohorts. The net benefit versus the threshold probability was plotted. The x-axis represented the threshold probability, while the y-axis represented the net benefits. The sensitivity and specificity of the model were calculated at each threshold to determine the net benefit. The DCAs showed that the net benefits of the mixed model (black line) and the radiomics model (red line) were superior to the benefits of the traditional model (blue line) with the threshold probability range from 0 to 1.

Discussion

SCSTs are composed of a heterogeneous group with complicated image manifestation. Although a few cases are mainly cystic, most SCSTs are predominantly solid masses (2, 4, 5). In MRI, the presence of solid tissue in an adnexal lesion is the primary cause of increased risk stratification (7, 14), so the particular lesion that requires to be differentiated is EOC. With the optimization of various diagnostic schemes, multiple methods have been gradually combined to evaluate adnexal tumors (15). In the current study, we assessed the ability of the clinical model, conventional MR model, traditional model, and the radiomics model to distinguish SCSTs and EOCs. As shown in the result, radiomics model yielded prominent predictive performance, significantly higher than that of the traditional model. The mixed model stood out among all the models.

Generally, SCSTs are considered to be clinically different from epithelial tumors to some extent, they occur across a wide age range and some patients may have clinical signs of hormone production (5, 16). In our study, only 8 in 73 patients with SCSTs showed elevated hormone levels, mainly presenting as a slight increase in prolactin. No significant difference was observed in menstrual status or endocrine level between the two groups. In previous studies (17, 18), serum CA125 and ROMA index have been proven to be effective diagnostic markers for EOCs, which could be used to differentiate between benign and malignant lesions. The present results were consistent with the previous studies, revealing that serum CA125 and ROMA index in the SCSTs group were significantly lower than those in EOCs, owing to the blunt nature of SCSTs. And in the multivariate logistic regression analysis, ROMA index was an independent predictor. However, the diagnostic efficacy of the clinical model was not satisfactory (AUCs = 0.729 and 0.680 in the development and testing cohorts, respectively), suggesting that clinical characteristics can only provide limited information for the differential diagnosis.

Conventional MR images could provide abundant information regarding the pathological features of tumors. Yin et al. (19) reported their results of 36 thecomas/fibrothecomas and 40 malignant pelvic solid tumors via conventional MRI and DWI examinations. They found that signal intensity on T2WI, capsule, and the lowest ADC value were important indicators in discriminating thecomas/fibrothecomas from malignant pelvic solid tumors. Their study lacked other types of sex cord tumors, and they didn’t build predictive models. In our study, ADC value, solid and cystic components were independent predictors for differentiating EOCs and SCSTs. The ADC value of SCSTs was relatively higher than that of EOCs, and SCSTs were more likely to be dominated by solid masses. Conventional MR model showed that the discrimination ability of MR parameters was also limited. Moreover, after negotiation, most conventional MR parameters required subjective interpretation by the radiologist, resulting in a poor inter-observer agreement in the assessment of these features.

In this study, we used different combinations of the selected radiomics features to establish three radiomics models, all of which achieved good prediction efficiency. When combining the most efficient group of radiomics features and the traditional parameters, we found that the diagnostic performance of the final radiomics model was comparable to that of the mixed model. To date, in a few studies, radiomics have been applied to patients with adnexal tumors, such as evaluating the ability of texture features to characterize the histopathological classification of ovarian cancer (911) and to predict prognosis (8, 9, 20). However, SCSTs have been rarely discussed. To our knowledge, this study was the first one to establish an MR-based radiomics model focusing on the differentiation of SCSTs from others. The results are encouraging, and show great potential to improve the prediction accuracy for ovarian tumors and derive predictive imaging biomarkers. Misdiagnosis could thus be further avoided.

Radiomics extracts high-throughput quantitative data from medical images, which is helpful for disease diagnosis, staging, management, and prognostication. In the current study, we ran a 5-fold cross validation setup during the feature selection, where each fold we did feature selection and model building to provide stable features. Nine radiomics features were finally selected, mainly including first-order statistical features and Gray Level Size Zone Matrix (GLSZM). The features contain a variety of traceable image information. First-order statistics are based on histograms of the original image and describe the distribution of voxel intensities within the image region (21, 22). The GLSZM features quantify the Second-Order joint probabilities of images (21). For example, Zone Entropy (ZE) of GLSZM measures the randomness in the distribution of zone sizes and gray levels, a higher value indicates more heterogeneity in the texture patterns; High Gray Level Zone Emphasis of GLSZM measures the distribution of the higher gray-level values, with a higher value indicating a greater proportion of higher gray-level values and size zones in the image. In our study, GLSZM is a quite important set of features for differentiating SCSTs and EOCs, the intrinsic definition of these features imply that they may capture the presence of necrotic, edematous or cellular regions within the tumors (22).

To date, there is no consensus on which or how many MRI sequences should be used to establish a radiomics model. Some studies (11, 15, 23) have included several individual sequences at the same time or combined multiple sequences to establish a radiomics model for evaluation, suggesting that the multi-sequence combination model may have better performance. Only FS-T2WI sequences were used in our radiomic analysis, considering that FS-T2WI is an important and common sequence in the conventional pelvic MRI scanning protocol with the highest spatial resolution which would improve the visualization of ROIs. Sufficient predictive performance has been achieved for differentiating SCSTs and EOCs. Further research with more sequences such as DWI, ADC, and T1WI+C, would be carried out to improve our radiomics classifier.

Several limitations should be noted. First, for a radiomics study, the sample size was relatively small, so the results might be biased. Second, this study is a single-center study without external verification, so the reproducibility and generalizability of the models need to be further tested. In the future, a multicenter study with a larger dataset size should be conducted to perform an optimal radiomics analysis. Third, the majority of patients in the EOC group were in advanced stages, resulting in a demonstrated advantage of the radiomics model. Future research should include more indeterminate adnexal masses.

In conclusion, by comparing various models, we found the MR-based radiomics model achieved excellent prediction performance for differentiating SCSTs and EOCs. The mixed model which combining the radiomics features and traditional parameters achieve a performance comparable to the radiomics model. Therefore, we believe that the radiomics approach could be a more objective and accurate way for discriminating SCSTs and EOCs. Meanwhile, the mixed model developed in our study could provide a comprehensive, effective manner for clinicians to diagnose and develop appropriate management strategies.

Data availability statement

The datasets presented in this article are not readily available because of institutional restriction. Requests to access the datasets should be directed to Xin Zhao, E-mail:zhaoxinxct@vip.163.com.

Ethics statement

The studies involving human participants were reviewed and approved by The Third Affiliated Hospital of Zhengzhou University. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author contributions

Conception and design: MC. Administrative support: XZ. Provision of study materials or patients: ST, ZZ, and TP. Collection and assembly of data: TR, LZ, ZY, and LM. Data analysis and interpretation: MC and XY. Writing - review and editing: MC, ST, and KW. Final approval of manuscript: All authors. All authors contributed to the article and approved the submitted version.

Funding

This research was funded by Henan Province Medical Science and Technology Research Plan Joint Construction Project (Grant No. LHGJ20190380).

Conflict of interest

KW was employed by General Electric GE Healthcare. XY was employed by Huiying Medical Technology.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.1073983/full#supplementary-material

References

1. Horta M, Cunha TM. Sex cord-stromal tumors of the ovary: a comprehensive review and update for radiologists. Diagn Interv Radiol (2015) 21(4):277–86. doi: 10.5152/dir.2015.34414

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Javadi S, Ganeshan DM, Jensen CT, Iyer RB, Bhosale PR. Comprehensive review of imaging features of sex cord-stromal tumors of the ovary. Abdom Radiol (NY) (2021) 46(4):1519–29. doi: 10.1007/s00261-021-02998-w

PubMed Abstract | CrossRef Full Text | Google Scholar

3. McCluggage WG, Lax SF, Longacre TA, Malpica A, Soslow RA. WHO classification of tumors: Female genital tumors. 5th edition. Lyon: Internal Agency for Research on Cancer (IARC) (2020).

Google Scholar

4. Zhao SH, Li HM, Qiang JW, Wang DB, Fan H. The value of MRI for differentiating benign from malignant sex cord-stromal tumors of the ovary: emphasis on diffusion-weighted MR imaging. J Ovarian Res (2018) 11(1):73. doi: 10.1186/s13048-018-0444-6

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Al Harbi R, McNeish IA, El-Bahrawy M. Ovarian sex cord-stromal tumors: an update on clinical features, molecular changes, and management. Int J Gynecol Cancer (2021) 31(2):161–8. doi: 10.1136/ijgc-2020-002018

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Nougaret S, McCague C, Tibermacine H, Vargas HA, Rizzo S, Sala E. Radiomics and radiogenomics in ovarian cancer: a literature review. Abdom Radiol (NY) (2021) 46(6):2308–22. doi: 10.1007/s00261-020-02820-z

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Sadowski EA, Thomassin-Naggara I, Rockall A, Maturen KE, Forstner R, Jha P, et al. O-RADS MRI risk stratification system: guide for assessing adnexal lesions from the ACR O-RADS committee. Radiology (2022) 303(1):35–47. doi: 10.1148/radiol.204371

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Vargas HA, Veeraraghavan H, Micco M, Nougaret S, Lakhman Y, Meier AA, et al. A novel representation of inter-site tumor heterogeneity from pre-treatment computed tomography textures classifies ovarian cancers by clinical outcome. Eur Radiol (2017) 27(9):3991–4001. doi: 10.1007/s00330-017-4779-y

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Zhang H, Mao Y, Chen X, Wu G, Liu X, Zhang P, et al. Magnetic resonance imaging radiomics in categorizing ovarian masses and predicting clinical outcome: a preliminary study. Eur Radiol (2019) 29(7):3358–71. doi: 10.1007/s00330-019-06124-9

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Song XL, Ren JL, Zhao D, Wang L, Ren H, Niu J. Radiomics derived from dynamic contrast-enhanced MRI pharmacokinetic protocol features: the value of precision diagnosis ovarian neoplasms. Eur Radiol (2021) 31(1):368–78. doi: 10.1007/s00330-020-07112-0

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Jian J, Li Y, Pickhardt PJ, Xia W, He Z, Zhang R, et al. MR image-based radiomics to differentiate type iota and type IotaIota epithelial ovarian cancers. Eur Radiol (2021) 31(1):403–10. doi: 10.1007/s00330-020-07091-2

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Zhang Q, Peng Y, Liu W, Bai J, Zheng J, Yang X, et al. Radiomics based on multimodal MRI for the differential diagnosis of benign and malignant breast lesions. J Magn Reson Imaging (2020) 52(2):596–607. doi: 10.1002/jmri.27098

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Lathwal A, Arora C, Raghava GPS. Prediction of risk scores for colorectal cancer patients from the concentration of proteins involved in mitochondrial apoptotic pathway. PloS One (2019) 14(9):e0217527. doi: 10.1371/journal

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Imaoka I, Wada A, Kaji Y, Hayashi T, Hayashi M, Matsuo M, et al. Developing an MR imaging strategy for diagnosis of ovarian masses. Radiographics (2006) 26(5):1431–48. doi: 10.1148/rg.265045206

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Qian L, Ren J, Liu A, Gao Y, Hao F, Zhao L, et al. MR imaging of epithelial ovarian cancer: a combined model to predict histologic subtypes. Eur Radiol (2020) 30(11):5815–25. doi: 10.1007/s00330-020-06993-5

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Hanley KZ, Mosunjac MB. Practical review of ovarian sex cord-stromal tumors. Surg Pathol Clin (2019) 12(2):587–620. doi: 10.1016/j.path.2019.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Romagnolo C, Leon AE, Fabricio ASC, Taborelli M, Polesel J, Del Pup L, et al. HE4, CA125 and risk of ovarian malignancy algorithm (ROMA) as diagnostic tools for ovarian cancer in patients with a pelvic mass: an Italian multicenter study. Gynecol Oncol (2016) 141(2):303–11. doi: 10.1016/j.ygyno.2016.01.016

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Suri A, Perumal V, Ammalli P, Suryan V, Bansal SK. Diagnostic measures comparison for ovarian malignancy risk in epithelial ovarian cancer patients: a meta-analysis. Sci Rep (2021) 11(1):17308. doi: 10.1038/s41598-021-96552-9

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Yin B, Li W, Cui Y, Chu C, Ding M, Chen J, et al. Value of diffusion-weighted imaging combined with conventional magnetic resonance imaging in the diagnosis of thecomas/fibrothecomas and their differential diagnosis with malignant pelvic solid tumors. World J Surg Oncol (2016) 14(1):5. doi: 10.1186/s12957-015-0760-x

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Wei W, Liu Z, Rong Y, Zhou B, Bai Y, Wei W, et al. A computed tomography-based radiomic prognostic marker of advanced high-grade serous ovarian cancer recurrence: A multicenter study. Front Oncol (2019) 9:255. doi: 10.3389/fonc.2019.00255

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Li W, Feng J, Zhu D, Xiao Z, Liu J, Fang Y, et al. Nomogram model based on radiomics signatures and age to assist in the diagnosis of knee osteoarthritis. Exp Gerontol (2022) 171:112031. doi: 10.1016/j.exger.2022.112031

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Yang Y, Yan LF, Zhang X, Nan HY, Hu YC, Han Y, et al. Optimizing texture retrieving model for multimodal MR image-based support vector machine for classifying glioma. J Magn Reson Imaging (2019) 49(5):1263–74. doi: 10.1002/jmri.26524

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Li C, Wang H, Chen Y, Zhu C, Gao Y, Wang X, et al. Nomograms of combining MRI multisequences radiomics and clinical factors for differentiating high-grade from low-grade serous ovarian carcinoma. Front Oncol (2022) 12:816982. doi: 10.3389/fonc.2022.816982

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: ovarian sex cord-stromal tumor, epithelial ovarian cancer, magnetic resonance imaging, radiomics, prediction model

Citation: Cheng M, Tan S, Ren T, Zhu Z, Wang K, Zhang L, Meng L, Yang X, Pan T, Yang Z and Zhao X (2023) Magnetic resonance imaging radiomics to differentiate ovarian sex cord-stromal tumors and primary epithelial ovarian cancers. Front. Oncol. 12:1073983. doi: 10.3389/fonc.2022.1073983

Received: 19 October 2022; Accepted: 28 December 2022;
Published: 13 January 2023.

Edited by:

Jasper Nijkamp, Aarhus University, Denmark

Reviewed by:

Wei Wei, Xi’an Polytechnic University, China
Niu Jinliang, Second Hospital of Shanxi Medical University, China

Copyright © 2023 Cheng, Tan, Ren, Zhu, Wang, Zhang, Meng, Yang, Pan, Yang and Zhao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xin Zhao, zhaoxinxct@vip.163.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.