Impact Factor 4.137 | CiteScore 4.28
More on impact ›

Original Research ARTICLE

Front. Oncol., 09 October 2019 | https://doi.org/10.3389/fonc.2019.01007

Multiplanar MRI-Based Predictive Model for Preoperative Assessment of Lymph Node Metastasis in Endometrial Cancer

Xiaojuan Xu1, Hailin Li2,3,4, Siwen Wang2,3, Mengjie Fang2,3, Lianzhen Zhong2,3, Wenwen Fan1, Di Dong2,3*, Jie Tian2,5* and Xinming Zhao1*
  • 1Department of Diagnostic Imaging, National Cancer Center/National Clinical Research Center for Cancer/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
  • 2CAS Key Lab of Molecular Imaging, Institute of Automation, Chinese Academy of Sciences, Beijing, China
  • 3School of Artificial Intelligence, University of Chinese Academy of Sciences, Beijing, China
  • 4School of Automation, Harbin University of Science and Technology, Harbin, China
  • 5Beijing Advanced Innovation Center for Big Data-Based Precision Medicine, School of Medicine, Beihang University, Beijing, China

Introduction: Assessment of lymph node metastasis (LNM) is crucial for treatment decision and prognosis prediction for endometrial cancer (EC). However, the sensitivity of the routinely used magnetic resonance imaging (MRI) is low in assessing normal-sized LNM (diameter, 0–0.8 cm). We aimed to develop a predictive model based on magnetic resonance (MR) images and clinical parameters to predict LNM in normal-sized lymph nodes (LNs).

Materials and Methods: A total of 200 retrospective patients were enrolled and divided into a training cohort (n = 140) and a test cohort (n = 60). All patients underwent preoperative MRI and had pathological result of LNM status. In total, 4,179 radiomic features were extracted. Four models including a clinical model, a radiomic model, and two combined models were built. Area under the receiver operating characteristic (ROC) curves (AUC) and calibration curves were used to assess these models. Subgroup analysis was performed according to LN size. All patients underwent surgical staging and had pathological results.

Results: All of the four models showed predictive ability in LNM. One of the combined models, ModelCR1, consisting of radiomic features, LN size, and cancer antigen 125, showed the best discrimination ability on the training cohort [AUC, 0.892; 95% confidence interval [CI], 0.834–0.951] and test cohort (AUC, 0.883; 95% CI, 0.786–0.980). The subgroup analysis showed that this model also indicated good predictive ability in normal-sized LNs (0.3–0.8 cm group, accuracy = 0.846; <0.3 cm group, accuracy = 0.849). Furthermore, compared with the routinely preoperative MR report, the sensitivity and accuracy of this model had a great improvement.

Conclusions: A predictive model was proposed based on MR radiomic features and clinical parameters for LNM in EC. The model had a good discrimination ability, especially for normal-sized LNs.

Introduction

Endometrial cancer (EC) is the most common gynecological malignancy in industrialized countries (1, 2). In China, EC is the second most common malignancy of the female genital tract with patients steadily increasing, especially in high urbanization areas (3). Lymph node metastasis (LNM) is an important risk factor for EC prognosis. Systematic lymphadenectomy is routinely performed according to International Federation of Gynecology and Obstetrics (FIGO). However, there is long-term controversy regarding whether it is necessary for low-risk or stage IA disease (4), as the incidence of LNM is very low in these patients (5). In addition, indiscriminative lymphadenectomy may lead to overtreatment and increase in post-operative complications, including chronic lymphedema, lymphocysts, infection, and nerve/vascular injuries (6).

Several histopathological findings, such as histological subtype, depth of myometrial invasion (DMI), primary tumor diameter (PTD), lymphovascular space invasion, and tumor grade, are known to be risk factors for LNM (4, 5), and researchers proposed various risk-classification models (4, 7). However, most of them are only available post-operatively. Sentinel lymph nodes mapping was proposed to evaluate LNM intraoperatively (8), but the technological dependence on experienced surgeons and relatively high false-negative rates limited its clinical application. Accurate preoperative and non-invasive evaluation of LNM is crucial, which can provide valuable information for prognosis prediction and treatment decision, especially in determining the extent of lymphadenectomy.

Magnetic resonance imaging (MRI) is a routinely used imaging modality for preoperative evaluation of EC. It plays an important role in assessing DMI (9), but its value for LNM assessment remains unsatisfactory, with reported sensitivities of 25–50% (10, 11). Radiomics, as a novel data mining technique, could extract high-dimensional quantitative features from medical images and select reliable features for the establishment of prediction models that could be used in computer-assisted decision support. Some recent researches showed that radiomics had the potential to evaluate therapeutic effects, predict the recurrence and metastasis, predict survival time (1214), and aid the differential diagnosis of cancers (15). Currently, radiomic investigations in preoperative prediction of LNM showed encouraging achievement (1618). However, to our knowledge, there is no literature that has determined whether a radiomics-based study would render superior prediction of metastasis in different size groups of LNs, and there has been no study on EC.

The purpose of this study was to investigate the efficacy of multiplanar enhanced MRI-based radiomics for preoperative prediction of metastasis in normal-sized (diameter 0–0.8 cm on MRI) LNs in EC patients.

Materials and Methods

Study Design and Participants

This retrospective study with anonymous data was approved by the Ethics Committee of our hospital, and the informed consent requirement was waived.

Two hundred consecutive patients with EC who had been treated between January 2011 and December 2017 were enrolled. Figure 1 shows the patient recruitment pathway. Patients were divided into two independent cohorts: 140 patients treated between January 2011 and March 2016 in the training cohort, and 60 patients treated between April 2016 and December 2017 in the test cohort.

FIGURE 1
www.frontiersin.org

Figure 1. Recruitment pathway for patients in this study.

As shown in Table 1, clinical parameters including age, blood serum cancer antigen 125 (CA125) level, preoperative histological type, and differentiation were derived from medical records.

TABLE 1
www.frontiersin.org

Table 1. Baseline characteristics of the training and test cohorts.

All MR imaging data were reviewed together by two board-certified radiologists (reader 1 and reader 2) specialized in gynecological tumor imaging with 6 and 20 years of experience. The PTD, DMI, involvement of the cervix, cornua, adnexa, parametrium, and LN status including the size and positive or negative were recorded. Maximal short-axis diameter of LN was measured on delayed phase of dynamic contrast enhanced (DCE) sequence at axial–sectional images (see details below). Patients with pelvic LN > 8 mm or abdominal LN > 10 mm, or with non-homogeneous enhancement and central necrosis on DCE images were regarded as MR report LN-positive (19). The consistency between the two radiologists was assessed by calculating the Cohen's kappa coefficients. Any disagreement was resolved by consultation. Note that LN status was defined by case.

MR Image Acquisition, Region of Interest Segmentation, and Radiomic Feature Extraction

Before receiving standard FIGO surgical staging, all patients underwent pelvic DCE MRI on two 3.0-T MR scanners (Signa HDxt and Discovery MR750, GE Medical Systems) with 8-channel phased array body coils. Two non-enhanced and one enhanced sequence were obtained and collected for analysis. Detailed scanning parameters are listed in Table 2.

TABLE 2
www.frontiersin.org

Table 2. Detailed acquired parameters in two MR scanners.

Tumor volume of interest (VOI), covering the whole tumor volume on each MR image, were manually segmented by reader 1 using ITK-SNAP software (www.itksnap.org, version 3.6.0). Radiomic feature extraction was performed with algorithms implemented in Python (www.python.org, version 3.6.5) (20). Three-dimensional radiomic features were extracted from the corresponding VOIs, including first-order statistics, shape-based, and texture features. More information about the radiomic feature extraction methodology can be found in Supplementary Method 1.

Surgery and Histopathologic Work-Up

All patients underwent FIGO surgical staging, and accepted template systematic lymphadenectomy. All lymph node specimens were processed and evaluated according to a standard protocol. Histologic analysis of each template lymph node dissection specimen included the following parameters: total number of histologically detected lymph nodes and number of positive nodes in each region as follows: external iliac, internal iliac and obturator, and common iliac. Note that the histopathologic LN status was still considered by case level in our analysis. The 2014 World Health Organization (WHO) classification (21) and the 2009 revised FIGO staging criteria for EC (22) were used for histological diagnosis, grading, and pathological staging.

Data Analysis

Feature Selection and Model Construction

Stability analysis of radiomic features between inter-/intra-reader segmentations was firstly carried out. Thirty patients were randomly chosen, and all of their images were segmented separately by the two radiologists, thereinto, reader 1 then re-segmented these images 1 week later. The intraclass correlation coefficients (ICCs) are usually adopted to assess the stability of radiomic features extracted from VOIs delineated by different readers or segmented by the same reader at different times. The radiomic features with ICC >0.75 were retained since they had good agreement between different segmentations.

Then, stability analysis between different versions of MR scanners on radiomic features was carried out. With all the patients randomly assigned to two MR scanners, Mann–Whitney U test was used to find out whether a radiomic feature showed statistical difference between different versions of MR scanners in the training cohort. We removed the radiomic features that had significant differences in the two versions of MR scanners, which would improve the generalization capability of our classifier.

Figure 2 shows the workflow of model development and decision-making process for model selection. Four models were constructed, including a clinical model with only clinical parameters (ModelC), a radiomic model with only radiomic features (ModelR), and two combined model (ModelCR1 and ModelCR2). After model evaluation, the final model was selected to be visualized as a clinical useful preoperative nomogram. The detailed construction processes of the four models were as follows.

FIGURE 2
www.frontiersin.org

Figure 2. Flow diagram of radiomic model construction. (A) MR images segmentation. The tumor region in each MRI slice was manually segmented, and then the whole tumor volume was reconstructed in order to extract 3D radiomic feature. (B) Radiomic feature extraction. Three types of radiomic features were extracted from tumor volume. (C) Feature selection process including stability, univariate analysis, and multivariate analysis. The construction of ModelC starts from univariate analysis. (D) Clinical application. After evaluating the four models, an optimal model was selected to plot nomogram for clinical computer-assisted decision support. MRI, magnetic resonance image; 3D, three-dimensional.

ModelC

The original feature set of the ModelC consisted of all of the 10 clinical parameters, including age, CA125, tumor pathologic type and differentiation by biopsy, tumor long-axis diameter, DMI, the ratio of tumor infiltration depth to myometrium depth, LN size, and adnexa or other organ involvement, which were all observed on MR images. After feature selection via Mann–Whitney U test and the least absolute shrinkage and selection operator (LASSO) method, six features were retained to fit the ModelC. The logistic regression model was constructed to examine the ability of the clinical parameters in classifying LNM.

ModelR

The stable original feature set of the ModelR consisted of 3,040 radiomic features that were dimensionally reduced by stability analysis. Univariate analysis (Mann–Whitney U test and chi-square test) was performed to evaluate the difference in LNM status, and the distribution of the p-values for the radiomic features is shown in Supplementary Figure 1. The retained significant features were then ranked by minimum redundancy/maximum relevance (mRMR), and the most redundant features were removed. The LASSO method was used to prevent overfitting. Two radiomic features were finally selected to fit the ModelR. Similarly, logistic regression model was constructed to examine the ability of the radiomic features in classifying LNM.

ModelCR1

All of the clinical parameters and retained 3,040 radiomic features formed the stable feature set. Univariate analysis was performed on this feature set. After removing features with p-values > 0.05, we computed the mRMR ranking for the feature set based on the concordance index (23, 24), and the top 5% features were retained. Then, the recursive feature elimination (RFE) method was performed to further select the LNM-related features. The RFE algorithm repeatedly constructed the model and removed the features, depending on the root mean square error of the model by a cross-validation in the training cohort.

ModelCR2

The stable original feature set of the ModelCR2 consisted of stable radiomic features and clinical parameters except LN size. Univariate analysis was first performed on this feature set. Then, in the multivariable analysis, mRMR and the LASSO method were performed successively. Logistic regression model was constructed to examine the classification ability of the combination of radiomic features and clinical primary lesion information.

Assessment and Validation of Model Performance

The receiver operating characteristic (ROC) curves were plotted to assess the performance of the four models in both cohorts. Area under ROC curve (AUC) was calculated for quantitative comparison. The model with the highest AUC was selected as the final model. Delong test was used to compare AUCs between the training cohort and test cohort, and a p-value > 0.05 indicated that there was no significant difference in AUCs, which ensured that the model had an enough low risk of over-fitting. Calibration curve was plotted to evaluate the agreement between prediction result and gold standard.

In previous research, good effects were gained in predicting the metastasis of an EC-LN larger than 1 cm (25). However, there is no study that ever focused on metastasis prediction on different sized LNs. So, we carried out a subgroup analysis on LN size. Patients were divided into three subgroups according to the LN size measured on MRI, including enlarged LNs with diameter larger than 0.8 cm (>0.8 cm), normal-sized LNs with diameter between 0.3 and 0.8 cm (0.3–0.8 cm), and normal-sized LNs with diameter smaller than 0.3 cm (<0.3 cm). F-score (F1=2Recall×PrecisionRecall+Precision) was calculated in these subgroups, assuming that recall (equivalently, sensitivity, TPTP+FN) and precision (equivalently, PPV, positive predictive value, TPTP+FP) are of equal importance, where TP, FN, and FP represent true positive, false negative, and false positive, respectively. The higher F-score synthetically reflects higher sensitivity and higher PPV.

Clinical Utility of the Final Model

In order to determine the clinical significance of the final model, decision curves were plotted by quantifying the net benefits in the training and test cohort. For the convenience of clinical application, a visualized preoperative nomogram was developed based on the formula exported by the logistic regression of the final model.

Statistical Analysis

In this study, statistical analysis programs were completed by R software (version 3.5.0; https://www.r-project.org). All statistical hypothesis tests were two-sided, and p-values < 0.05 were considered significant.

Results

Patient Characteristics

The clinical and pathological characteristics in the two cohorts are shown in Tables 1, 3, respectively. Pathologically LN positive (pN+) patients formed 37.14% (52/140) and 25.00% (15/60) of the training and test cohorts, respectively, and there was no significant difference between them (p-value = 0.133, χ2 test). The clinical parameters age and CA125 had no differences between the two cohorts (p-value = 0.077 and 0.539 respectively, Mann–Whitney U-test). In total, sensitivity and specificity were 28.36% (19/67) and 93.98% (125/133) according to the MR report LN status after consensus within two radiologists in our study. Also, the judgments by two radiologists on MRI were basically stable (sensitivity was 0.642 and 0.552, and the specificity was 0.917 and 0.940, respectively). The inconsistency of judgment was resolved by consultation. The Cohen's kappa coefficients to test consistency of the main MR indicators evaluated by the two radiologists are listed in Supplementary Table 1.

TABLE 3
www.frontiersin.org

Table 3. Pathological characteristics of the patients in our study.

Feature Selection and Model Construction

In total, 1,393 radiomic features were extracted from each of the three MR scanning sequences. Then, 4,179 radiomic features were reduced to 3,040 by stability analysis.

In ModelC, six clinical parameters were selected including CA125, tumor differentiation by biopsy, DMI, the ratio of tumor infiltration depth to myometrium depth, LN size, and adnexa involvement, which were all observed on MR. In ModelR, two radiomic features were selected including correlation and HGLE. In ModelCR1, four risk factors including two clinical parameters (CA125 and LN size) and two radiomic features (correlation and HGLE) were used to build the prediction model (Figure 3) (13). The two radiomic features were extracted from the delayed phase of the 3D-Iso-LAVA and sagittal T2WI FSE, respectively. In ModelCR2, the LN size was removed and the same other three indicators (CA125, correlation, and HGLE) were selected. The detailed calculation formulas for ModelCR1 and ModelCR2 were given in Supplementary Method 2.

FIGURE 3
www.frontiersin.org

Figure 3. (A) LASSO coefficient profiles of the clinical parameters in ModelC. According to the 1 standard error of the minimum criteria (the 1-SE criteria), the dotted line was plotted at the selected log(λ) (−2.914) via 10-fold cross-validation. (B) LASSO coefficient profiles of the radiomic features in ModelR. A log(λ) value of −1.750 was chosen (10-fold cross-validation, 1-SE criteria). (C) Feature selection using the RFE method in ModelCR1. The rank of feature importance was obtained using the random forest method; RFE built the model continuously by eliminating the lower ranking feature. The RMSE was used to select the optimal feature set in a 10-fold cross-validation. (D) LASSO coefficient profiles of the combined feature set in ModelCR2. A log(λ) value of −1.983 was chosen (10-fold cross-validation, 1-SE criteria). RFE, recursive feature elimination; RMSE, root mean square error.

Assessment of Predictive Models

ModelCR1 showed a significant ability in detecting pN+ with an AUC of 0.892 [95% confidence interval [CI]: 0.834–0.951] in the training cohort and an AUC of 0.883 (95% CI: 0.786–0.980) in the test cohort. Nomogram (Figure 4C) was established for ModelCR1. The p values calculated from Delong tests were 0.875, 0.8416, 0.7008, and 0.5865 for ModelCR1, ModelCR2, ModelR, and ModelC, respectively, indicating that there were no significant differences in AUCs between the training cohort and test cohort for each model. Performances of the four models in the training and test cohort are shown in Figures 4A,B. Based on the threshold determined by Youden's index in the training cohort, we used net reclassification index (NRI) to analyze the improvement brought by ModelCR1 compared with other models. The results showed that ModelCR1 outperformed ModelR (NRI = 0.306, P < 0.001), ModelC (NRI = 0.134, P = 0.010), and ModelCR2 (NRI = 0.090, P = 0.077). Meanwhile, ModelCR1 also significantly surpassed MR reports by radiologists (NRI = 0.489, P = 0.006). Besides, the calibration curves were plotted in both cohorts for further performance evaluation of ModelCR1 (Figures 5A,B). Calibration curves show good fitness for probability of LNM (Hosmer–Lemeshow test, p-value = 0.961 in the training cohort, 0.803 in the test cohort). Figure 5C shows patients' risk scores calculated from ModelCR1, intuitively indicating its high classification ability.

FIGURE 4
www.frontiersin.org

Figure 4. (A) ROC of the four models in training cohort. (B) ROC of the four models in test cohort. (C) Preoperative nomogram of ModelCR1. ROC, receiver operating characteristic.

FIGURE 5
www.frontiersin.org

Figure 5. (A) The calibration curve in the training cohort. (B) The calibration curve in the test cohort. (C) Patient risk score output by ModelCR1, while red bars show scores for those who were pathologically LN(–).

As shown in Figure 6A, in the subgroup of enlarged LNs, ModelCR1 achieved the highest sensitivity of 0.970, equal to that predicted by MR report. In the subgroup of normal-sized LNs (0.3–0.8 cm), ModelCR1 displayed the highest accuracy of 0.846 and a sensitivity of 0.647, which far surpassed the MR report (accuracy, 0.785; sensitivity, 0.235). In the subgroup of normal-sized LNs (<0.3 cm), ModelCR1 showed the best accuracy of 0.849 and a moderate sensitivity of 0.471, however, still greatly outperforming the MR report (accuracy, 0.817; sensitivity, 0.000). Meanwhile, F-score and accuracy in three subgroups are shown in Figures 6B,C, respectively. The highest F-score and most powerful accuracy of ModelCR1 were reflected among the five predictive models.

FIGURE 6
www.frontiersin.org

Figure 6. (A) Heatmap of showing the classifiers' performance in different LN size subgroups. A deeper blue indicates a larger value. (B) Line chart of F-score, (C) histogram of accuracy in three subgroups. Red, blue, green, yellow, and purple lines and boxes, respectively, represent ModelCR1, MR report, ModelCR2, ModelR, and ModelC.

Discussion

In the present study, we developed four predictive models based on multiplanar DCE MR images and clinical parameters for LNM in EC patients. ModelCR1, which consisted of radiomic features, LN size, and CA125, showed the best discrimination ability, especially in patients with normal-sized LNs (diameter, 0–0.8 cm on MRI) and the sensitivity was greatly improved compared with the routine MR reports. The high F-scores indicated that while the sensitivity increased significantly, the PPV remained high.

A non-invasive and convenient preoperative assessment for LNM is crucial for EC treatment decision and prognosis prediction. Patients' data from preoperative procedures such as MRI, biopsy, and CA125 have been studied to assess LNM in recent years. MRI still remained the cornerstone in LN assessment in EC, showing satisfactory specificity but relatively low sensitivity. The combination of relative apparent diffusion coefficient (ADC) value and LN size was reported to result in a significant increase in sensitivity from 25 to 83% compared with conventional MRI (26); however, there have been conflicting reports in the literature regarding the detection of LNM at diffusion-weighted MR imaging (DWI). Nakai et al. (27) used 1.5-T MRI to evaluate nodal ADC values in gynecologic malignancies and were unable to differentiate benign from malignant LNs. Wang et al. (28) proposed a tumor biomarker predictive method by combining human epididymis secretory protein 4 (HE4) and CA125, achieving a high sensitivity of 94.1% but a low specificity of 30.7%. Notably, there is no clearly defined HE4 cutoff value for EC at present. In our study, we incorporated CA125 in our models, which was more generally accepted than HE4. Kang et al. (29) developed a low-risk prediction model for LNM based on MRI and serum CA125 data in endometrioid-type EC patients, and obtained sensitivity and specificity of 84.9 and 55.5%, respectively. Here, three MRI parameters including DMI, LN enlargement, and extension beyond uterine corpus were identified to be independent risk factors for LNM. In our study, we obtained CA125 and MR report LN size as risk factors for EC LNM prediction, which was similar to that result, and showed good discrimination ability on both cohorts and different LN size subgroups, especially for those normal-sized LNs, which previous researches had not yet focused on.

In our study, we collected and analyzed all available preoperative clinical parameters and established four prediction models. We aimed to determine the prediction efficiency of different models compared with the MR report in different sized LN subgroups. MRI uses several common morphological criteria in differentiating benign from malignant nodes (30) but nodal size still remains the commonly accepted standard. Low sensitivity is a recognized limitation when nodal size criteria are used on cross-sectional imaging, especially for normal-sized nodes due to limited spatial resolution. In this study, with node size gradually decreasing, the MR report and ModelC showed a decreasing sensitivity, whereas the ModelR and ModelCR2 were more stable because of the high sensitivity in each sized LN subgroup (Supplementary Figure 2). The performance of the above classifiers confirmed our thoughts: When the LN was normal sized on MRI, combining LN size in classifiers could improve prediction accuracy but greatly reduce sensitivity. It is already accepted that normal-sized LNs may also contain metastases (31). The results of the MR report rely too much on LN size so that when LN size is normal on MRI, the sensitivity becomes very low. The concept remains the same when LN size is enlarged (>0.8 cm), then the specificity becomes very low. This can be due to the fact that it is usually difficult to differentiate enlarged nodes because of benign pathology, such as infection, granulomatous disease, and reactive hyperplasia vs. malignant disease (30). LN size was not a significant predictor for the fusion model based on radiomics, although it could improve the predictive accuracy. This may indicate that there was enough information contained in the primary tumor region that could detect LNM. It is feasible to predict LNM status without the dependence on information of LNs.

Although ModelCR1 showed a slightly lower sensitivity than ModelR and ModelCR2 in the normal-sized group, its accuracy in each group is the highest, and the F-score with normal-sized LNs is greatly improved. Therefore, it was proposed as the optimal prediction model. To our knowledge, this is the first subgroup analysis on different sized LNs with preoperative nomogram study in EC.

Due to a variety of MR scanner parameter settings and scanner models, it is difficult to guarantee different scanners with exactly the same imaging quality, thus making it difficult to ensure the stability of radiomic features. By eliminating the radiomic features sensitive to scanner models and parameter settings in the training cohort, the radiomic model generalization ability can be improved. The radiomic texture features (correlation and HGLE) selected in ModelCR1, ModelCR2, and ModelR reflected two kinds of heterogeneity of VOI with a Pearson correlation coefficient of 0.095. Correlation shows the linear dependency of gray-level values to the corresponding voxels in the gray-level cooccurence matrix (GLCM) of MRI. HGLE is a measure of the proportion of areas with higher gray values in the tumor. These two radiomic features indicated that the extent of heterogeneity of tumor is associated with LNM. The more heterogeneous the tumor, the higher the risk of LNM.

The limitations of the present study include two aspects. First, there was no external validation. Multicenter investigation with a larger dataset was needed to further validate the generalization ability of our model. Second, genomic information was not yet incorporated into our models. A combination of gene marker panels and radiomic features will be promising in evaluation of EC.

In conclusion, our study presented a predictive model based on multiplanar contrast enhanced MR images and incorporated both the radiomic features and clinical parameters, which showed good predictive accuracy for preoperative LNM in EC, especially in patients with normal-sized LNs.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

Ethics Statement

The studies involving human participants were reviewed and approved by National Cancer Center/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

XX, HL, XZ, DD, JT, MF, and SW: conception and design. XX, HL, WF, and LZ: collection and assembly of data. HL, XX, DD, XZ, and MF: data analysis and interpretation. All authors: manuscript writing and final approval of manuscript.

Funding

This work was supported by the Beijing Hope Run Special Fund of Cancer Foundation of China (grant LC2016B01), the National Natural Science Foundation of China (81971776, 81771924, 81501616, 81671851, 81527805, and 61671449), National Key R&D Program of China (2017YFA0205200, 2017YFC1308700, 2017YFC1308701, 2017YFC1309100, and 2016YFC0103803), the Beijing Municipal Science and Technology Commission (Z171100000117023, Z161100002616022), the Beijing Natural Science Foundation (L182061), and the Youth Innovation Promotion Association CAS (2017175).

Conflict of Interest

The 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.

Supplementary Material

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

Abbreviations

LNM, lymph node metastasis; EC, endometrial cancer; LN, lymph node; AUC, area under the curve; MRI, magnetic resonance imaging; IGO, International Federation of Gynecology and Obstetrics; DMI, depth of myometrial invasion; A125, cancer antigen 125; PTD, primary tumor diameter; WHO, World Health Organization; DCE, dynamic contrast enhanced; VOI, volume of interest; ICC, intraclass correlation coefficient; mRMR, minimum redundancy/maximum relevance; RFE, recursive feature elimination; HGLE, high gray-level emphasis; LASSO, least absolute shrinkage and selection operator; ROC, receiver operating characteristic; PV, positive predictive value; TP, true positive; FN, false negative; FP, false positive; pN+, pathologically LN positive; 3D-iso-LAVA-XV, three-dimensional liver acquisition with volume acceleration DCE with isotropy scanning; CI, confidence interval; ADC, apparent diffusion coefficient; DWI, diffusion-weighted MR imaging; HE4, human epididymis secretory protein 4; GLCM, gray-level cooccurrence matrix.

References

1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018. CA Cancer J Clin. (2018) 68:7–30. doi: 10.3322/caac.21442

CrossRef Full Text | Google Scholar

2. Ghanem AI, Khan NT, Mahan M, Ibrahim A, Buekers T, Elshaikh MA. The impact of lymphadenectomy on survival endpoints in women with early stage uterine endometrioid carcinoma: a matched analysis. Eur J Obstet Gynecol Reprod Biol. (2017) 210:225–30. doi: 10.1016/j.ejogrb.2016.12.033

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Chen W, Zheng R, Zhang S, Zeng H, Zuo T, Xia C, et al. Cancer incidence and mortality in China in 2013: an analysis based on urbanization level. Chin J Cancer Res. (2017) 29:1–10. doi: 10.21147/j.issn.1000-9604.2017.01.01

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Korkmaz V, Meydanli MM, Yalçin I, Sari ME, Sahin H, Kocaman E, et al. Comparison of three different risk-stratification models for predicting lymph node involvement in endometrioid endometrial cancer clinically confined to the uterus. J Gynecol Oncol. (2017) 28:e78. doi: 10.3802/jgo.2017.28.e78

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Papathemelis T, Scharl S, Kronberger K, Gerken M, Scharl A, Pauer A, et al. Survival benefit of pelvic and paraaortic lymphadenectomy in high-grade endometrial carcinoma: a retrospective population-based cohort analysis. J Cancer Res Clin Oncol. (2017) 143:2555–62. doi: 10.1007/s00432-017-2508-1

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Biglia N, Librino A, Ottino MC, Panuccio E, Daniele A, Chahin A. Lower limb lymphedema and neurological complications after lymphadenectomy for gynecological cancer. Int J Gynecol Cancer. (2015) 25:521–5. doi: 10.1097/IGC.0000000000000341

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Colombo N, Preti E, Landoni F, Carinelli S, Colombo A, Marini C, et al. ESMO Guidelines Working Group. Endometrial cancer: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann Oncol. (2013) 24(Suppl. 6):vi33–8. doi: 10.1093/annonc/mdt353

CrossRef Full Text | Google Scholar

8. Khoury-Collado F, St Clair C, Abu-Rustum NR. Sentinel lymph node mapping in endometrial cancer: an update. Oncologist. (2016) 21:461–6. doi: 10.1634/theoncologist.2015-0473

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Taufiq M, Masroor I, Hussain Z. Diagnostic accuracy of diffusion weighted magnetic resonance imaging in the detection of myometrial invasion in endometrial carcinoma. J Coll Physicians Surg Pak. (2016) 26:13–7.

PubMed Abstract | Google Scholar

10. Kim HJ, Cho A, Yun M, Kim YT, Kang WJ. Comparison of FDG PET/CT and MRI in lymph node staging of endometrial cancer. Ann Nucl Med. (2016) 30:104–13. doi: 10.1007/s12149-015-1037-8

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Teng F, Zhang YF, Wang YM, Yu J, Lang X, Tian WY, et al. Contrast-enhanced MRI in preoperative assessment of myometrial and cervical invasion, and lymph node metastasis: diagnostic value and error analysis in endometrial carcinoma. Acta Obstet Gynecol Scand. (2015) 94:266–73. doi: 10.1111/aogs.12570

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Dong D, Tang L, Li ZY, Fang MJ, Gao JB, Shan XH, et al. Development and validation of an individualized nomogram to identify occult peritoneal metastasis in patients with advanced gastric cancer. Ann Oncol. (2019) 30:431–8. doi: 10.1093/annonc/mdz001

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Liang W, Xu L, Yang P, Zhang L, Wan D, Huang Q, et al. Novel nomogram for preoperative prediction of early recurrence in intrahepatic cholangiocarcinoma. Front Oncol. (2018) 8:360. doi: 10.3389/fonc.2018.00360

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Song J, Shi J, Dong D, Fang M, Zhong W, Wang K, et al. A new approach to predict progression-free survival in stage IV EGFR-mutant NSCLC patients with EGFR-TKI therapy. Clin Cancer Res. (2018) 24:3583–92. doi: 10.1158/1078-0432.CCR-17-2507

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Li L, Wang K, Ma X, Liu Z, Wang S, Du J, et al. Radiomic analysis of multiparametric magnetic resonance imaging for differentiating skull base chordoma and chondrosarcoma. Eur J Radiol. (2019) 118:81–7. doi: 10.1016/j.ejrad.2019.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Lu W, Zhong L, Dong D, Fang M, Dai Q, Leng S, et al. Radiomic analysis for preoperative prediction of cervical lymph node metastasis in patients with papillary thyroid carcinoma. Eur J Radiol. (2019) 118:231–8. doi: 10.1016/j.ejrad.2019.07.018

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Dong Y, Feng Q, Yang W, Lu Z, Deng C, Zhang L, et al. Preoperative prediction of sentinel lymph node metastasis in breast cancer based on radiomics of T2-weighted fat-suppression and diffusion-weighted MRI. Eur Radiol. (2018) 28:582–91. doi: 10.1007/s00330-017-5005-7

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Shen C, Liu Z, Wang Z, Guo J, Zhang H, Wang Y, et al. Building CT radiomics based nomogram for preoperative esophageal cancer patients lymph node metastasis prediction. Transl Oncol. (2018) 11:815–24. doi: 10.1016/j.tranon.2018.04.005

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Meissnitzer M, Forstner R. MRI of endometrium cancer—How we do it. Cancer Imaging. (2016) 16:11. doi: 10.1186/s40644-016-0069-1

CrossRef Full Text | Google Scholar

20. van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, et al. Computational radiomics system to decode the radiographic phenotype. Cancer Res. (2017) 77:e104–7. doi: 10.1158/0008-5472.CAN-17-0339

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Kurman RJ, Carcangiu ML, Herrington CS, Young RH. WHO Classification of Tumors of Female Reproductive Organs [M]. 4th edn. Lyon: IARC Press (2014).

PubMed Abstract | Google Scholar

22. Adapted from FIGO Committee on Gynecologic cancers. FIGO classification of cancer of the vulva, cervix, and corpus uteri. Int J Gynecol Obstetr. (2014) 115:97–8. doi: 10.1016/j.ijgo.2014.02.003

CrossRef Full Text | Google Scholar

23. Harrel FE Jr, Lee KL, Mark DB. Tutorial in biostatistics: multivariable prognostic models: issues in developing models, evaluating assumptions and adequacy, and measuring and reducing error. Stat Med. (1996) 15:361–87. doi: 10.1002/(SICI)1097-0258(19960229)15:4<361::AID-SIM168>3.0.CO;2-4

CrossRef Full Text | Google Scholar

24. Pencina MJ, D'Agostino RB. Overall C as a measure of discrimination in survival analysis: model specific population value and confidence interval estimation. Stat Med. (2004) 23:2109–23. doi: 10.1002/sim.1802

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Koplay M, Dogan NU, Erdogan H, Sivri M, Erol C, Nayman A, et al. Diagnostic efficacy of diffusion-weighted MRI for pre-operative assessment of myometrial and cervical invasion and pelvic lymph node metastasis in endometrial carcinoma. J Med Imaging Radiat Oncol. (2014) 58:538–46. doi: 10.1111/1754-9485.12209

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Lin G, Ho KC, Wang JJ, Ng KK, Wai YY, Chen YT, et al. Detection of lymph node metastasis in cervical and uterine cancers by diffusion-weighted magnetic resonance imaging at 3T. J Magn Reson Imaging. (2008) 28:128–35. doi: 10.1002/jmri.21412

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Nakai G, Matsuki M, Inada Y, Tatsugami F, Tanikake M, Narabayashi I, et al. Detection and evaluation of pelvic lymph nodes in patients with gynecologic malignancies using body diffusion-weighted magnetic resonance imaging. J Comput Assist Tomogr. (2008) 32:764–8. doi: 10.1097/RCT.0b013e318153fd43

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Wang Y, Han C, Teng F, Bai Z, Tian W, Xue F. Predictive value of serum HE4 and CA125 concentrations for lymphatic metastasis of endometrial cancer. Int J Gynaecol Obstet. (2017) 136:58–63. doi: 10.1002/ijgo.12010

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Kang S, Nam JH, Bae DS, Kim JW, Kim MH, Chen X, et al. Preoperative assessment of lymph node metastasis in endometrial cancer: a Korean Gynecologic Oncology Group study. Cancer. (2017) 123:263–72. doi: 10.1002/cncr.30349

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Lai G, Rockall AG. Lymph node imaging in gynecologic malignancy. Semin Ultrasound CT MR. (2010) 31:363–76. doi: 10.1053/j.sult.2010.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Thoeny HC, Froehlich JM, Triantafyllou M, Huesler J, Bains LJ, Vermathen P, et al. Metastases in normal-sized pelvic lymph nodes: detection with diffusion-weighted MR imaging. Radiology. (2014) 273:125–35. doi: 10.1148/radiol.14132921

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: endometrial cancer, lymph node, metastasis, magnetic resonance imaging, radiomics

Citation: Xu X, Li H, Wang S, Fang M, Zhong L, Fan W, Dong D, Tian J and Zhao X (2019) Multiplanar MRI-Based Predictive Model for Preoperative Assessment of Lymph Node Metastasis in Endometrial Cancer. Front. Oncol. 9:1007. doi: 10.3389/fonc.2019.01007

Received: 28 July 2019; Accepted: 18 September 2019;
Published: 09 October 2019.

Edited by:

Chunxiao Guo, University of Texas MD Anderson Cancer Center, United States

Reviewed by:

Zaiyi Liu, Guangdong Provincial People's Hospital, China
Anna Colleen Crouch, University of Texas MD Anderson Cancer Center, United States

Copyright © 2019 Xu, Li, Wang, Fang, Zhong, Fan, Dong, Tian 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: Di Dong, di.dong@ia.ac.cn; Jie Tian, jie.tian@ia.ac.cn; Xinming Zhao, xinmingzh@sina.com

These authors have contributed equally to this work