Impact Factor 4.848 | CiteScore 3.5
More on impact ›

Original Research ARTICLE

Front. Oncol., 27 October 2020 | https://doi.org/10.3389/fonc.2020.541849

Radiomics Nomogram of DCE-MRI for the Prediction of Axillary Lymph Node Metastasis in Breast Cancer

Ning Mao1†, Yi Dai2†, Fan Lin1, Heng Ma1, Shaofeng Duan3, Haizhu Xie1, Wenlei Zhao1* and Nan Hong4*
  • 1Department of Radiology, Yantai Yuhuangding Hospital, Qingdao University, Yantai, China
  • 2Department of Radiology, Peking University Shenzhen Hospital, Shenzhen, China
  • 3Precision Health Institution, GE Healthcare, China, Shanghai, China
  • 4Department of Radiology, Peking University People’s Hospital, Beijing, China

Purpose: This study aimed to establish and validate a radiomics nomogram based on dynamic contrast-enhanced (DCE)-MRI for predicting axillary lymph node (ALN) metastasis in breast cancer.

Method: This retrospective study included 296 patients with breast cancer who underwent DCE-MRI examinations between July 2017 and June 2018. A total of 396 radiomics features were extracted from primary tumor. In addition, the least absolute shrinkage and selection operator (LASSO) algorithm was used to select the features. Radiomics signature and independent risk factors were incorporated to build a radiomics nomogram model. Calibration and receiver operator characteristic (ROC) curves were used to confirm the performance of the nomogram in the training and validation sets. The clinical usefulness of the nomogram was evaluated by decision curve analysis (DCA).

Results: The radiomics signature consisted of three ALN-status-related features, and the nomogram model included the radiomics signature and the MR-reported lymph node (LN) status. The model showed good calibration and discrimination with areas under the ROC curve (AUC) of 0.92 [95% confidence interval (CI), 0.87–0.97] in the training set and 0.90 (95% CI, 0.85–0.95) in the validation set. In the MR-reported LN-negative (cN0) subgroup, the nomogram model also exhibited favorable discriminatory ability (AUC, 0.79; 95% CI, 0.70–0.87). DCA findings indicated that the nomogram model was clinically useful.

Conclusions: The MRI-based radiomics nomogram model could be used to preoperatively predict the ALN metastasis of breast cancer.

Introduction

Breast cancer is a malignant tumor that endangers women’s health and quality of life. Axillary lymph node (ALN) is the first station of breast lymphatic drainage, which collects approximately 75% of breast lymph. Thus, ALN is the most easily metastasized site of breast cancer. ALN status is an important factor affecting the treatment of patients with breast cancer and is assess by the gold standards ALN dissection and sentinel lymph node (LN) biopsy. However, ALN dissection is invasive and has many complications, such as lymphedema, and sentinel LN biopsy is also invasive (1). Therefore, a non-invasive prediction tool for preoperative LN status is needed.

MRI has been widely used in breast examination because of its good soft tissue contrast, high sensitivity, and high negative predictive rate (2). Although this technique is superior to digital mammography and ultrasonography, its efficacy in identifying malignant nodes is unsatisfactory (24).

Radiomics can extract massive image features; transform medical images into high-dimensional and exploitable data; and use artificial intelligence to combine medical images, genes, and huge clinical data to establish a model that supports clinical decision-making and quantify tumor heterogeneity (59). This method has good clinical prospects (911). The combined analysis of multiple features including clinical ones is the most promising approach, especially for the clinical management of tumors (1216). Furthermore, nomograms, which allow the investigation of multiple features in parallel, transform complex regression equation into visual graphs (1719).

This study aimed to develop and validate a radiomics nomogram model based on dynamic contrast-enhanced (DCE)-MRI and clinical risk factors to determine its potential in predicting ALN metastasis in patients with breast cancer.

Materials and Methods

Patients

This retrospective study was approved by the Institutional Review Board. Inclusion criteria were as follows: (a) patients with breast cancer confirmed by histopathological examination, (b) available clinical information, and (c) surgery conducted after MR scanning. Exclusion criteria were as follows: (a) patients who underwent preoperative neoadjuvant chemotherapy or radiotherapy, (b) patients who underwent biopsy prior to MR scanning, and (c) patients with other tumors and (d) non-mass lesions without delineate boundaries. The patients were divided into two independent sets, namely, training (200 patients) and validation sets (96 patients).

Clinical data were obtained through the medical record systems. All images were reviewed by two radiologists with at least 10 years of experience in imaging diagnosis, and the largest diameter of the tumors, apparent diffusion coefficient (ADC) value, enhanced features and the short diameter of the largest LN were recorded. MRI-reported LN status refers to the imaging-based diagnosed LN status according to the radiologist. T2WI and DCE-MRI series were used for ALN diagnosis. A patient’s LN status was classified as positive (cN+) if one or more ALNs found on MR images met any one of the following MRI features: 1) visible ALN >10 mm in a short diameter, 2) ratio of the longest to shortest axes < 1.6, 3) eccentric cortical thickening, and 4) loss of fatty hilum. Those who did not met the above mentioned criteria and met the above criteria but showed no difference in terms of number, size, or shape compared with the contralateral ALN, the LN status was assumed to be negative (cN0) (2022). LN was classified as positive when at least one of the four criteria was satisfied. Agreement from MRI-reported LN-status analyzed from two observers was compared using κ statistic, in which 0 < κ ≤ 0.4 indicates poor agreement, 0.4 < κ < 0.75 indicates good agreement, and 0.75 ≤ κ < 1 indicates excellent agreement. All disagreements were resolved through consultation.

Pathological Evaluation

Pathology is the gold standard for LN metastasis. Radionuclide and methylene blue were used as tracers to ensure that all sentinel LNs were removed. The patients were injected with radionuclide 2−3 h prior to surgery. After anesthetization, methylene blue was injected into the patient’s breast, which was then gently rubbed to allow the dye to further spread along the lymph vessels. Radionuclide detector was used to identify LNs labeled with nuclide during surgery. Stained LNs were also searched along the blue-stained lymph-vessels from top to bottom, inside to outside, and toward the axilla. All LN specimens were fixed by 4% neutral formaldehyde, embedded in paraffin, sectioned in 4-um thickness, sequentially sectioned, and stained by hematoxylin and eosin stain. Finally, the morphology of LN tissues was observed by two pathologists under BX53 electron microscope, and the tumor cells were confirmed as LN metastasis. If the pathological result of LN biopsy was inconsistent with that of surgery, then the latter was used as the standard.

Histopathological information, such as histological grade, estrogen receptor, progesterone receptor, human epidermal growth factor receptor type 2, and Ki-67, was obtained from the medical record system. Threshold values were ≤ 1% for the estrogen receptor and progesterone receptor levels and ≤ 20% for Ki-67 (23).

MR Image Acquisition

Figure 1A presents the study flowchart. All images were obtained on a 3.0T MRI system (GE Discovery 750W) using an eight-channel breast-dedicated coil in prone position. The MRI sequences included axial T1-weighted imaging, axial T2-weighted imaging, DCE-MRI, and sagittal contrast-enhanced imaging.

FIGURE 1
www.frontiersin.org

Figure 1 Study flowchart (A) and radiomics workflow (B).

The scanning parameters were as follows: ① axial T1WI (TR = 460 ms, TE = 6.3 ms, slice thickness = 5 mm, slice spacing = 1 mm); ② axial fat suppression T2WI (TR = 5210 ms, TE = 84.7 ms, slice thickness = 5 mm, slice spacing = 1 mm); ③ axial DWI (SE-EPI sequence, TR = 2496 ms, TE = 71.9 ms, slice thickness = 5 mm, slice spacing = 1 mm, B = 0/800 s/mm2); ④ DCE scanning was performed on T1 fat suppression. The contrast medium was GD-DTPA, dose = 0.2 mmol/kg, TR = 5.7 ms, TE = 1.7 ms, slice thickness = 2 mm, slice spacing = 0 mm, FOV = 36 cm × 36cm, matrix = 288 × 320, phases = 8, and total time=6 min; ⑤ and sagittal contrast-enhanced imaging was performed after DCE (TR = 6.7 ms, TE = 1.7 ms, slice thickness = 2 mm, slice spacing = 0 mm, FOV = 28 cm × 28 cm, matrix = 200 × 256. Scan ranges for breast MRI were as follows: in the prone position, the bilateral breasts naturally hung over the center of the breast coil; in the horizontal axis position, the bilateral breasts were located in the center of the FOV, and the range included the entire bilateral breasts and bilateral axillary regions; and in the sagittal position, the positioning line was parallel to the long axis of the breast. All DICOM data were exported from Picture Archiving and Communication Systems.

Image Preprocessing

Image preprocessing was necessary prior to feature extraction. This process consisted of three steps, namely, standardization of the gray value of the region of interest (ROI), discretization of the gray level, and image resampling (2426).

Image Segmentation and Radiomics Feature Extraction

Figure 1B presents the radiomics workflow. DCE-MRI (the peak enhanced phase of the multiphase contrast-enhanced MRI selected in accordance with time intensity curve) was selected for radiomics analysis, and the primary tumor was manually segmented using 3D-ROI by two trained radiologists with at least 10 years of experience in breast imaging and who were also blinded to the LN status and pathologic results. Twenty-four breast lesions were randomly selected to calculate the intra- and inter-observer agreement of the feature extraction. First, the two radiologists extracted the radiomics features. After 2 weeks, reader 1 used the same method to extract the radiomics features. Inter- and intra-correlation coefficients (ICCs) were calculated to assess the reproducibility of the radiomics features, and ICCs > 0.80 were considered as good agreement. The remaining image segmentation was performed by reader 1. Image segmentation and radiomics feature extraction were performed on Artificial Intelligence Kit software (version 3.2.0; GE Healthcare, Shanghai, China).

Feature Selection and Radiomics Signature Building

The features with high repeatability (ICC > 0.80) were selected. Feature selection was then performed using the LASSO logistic regression method in the training set. This method is suitable for high-dimensional data (27). Radiomics score reflecting the risk of ALN metastasis was calculated for each patient by using a linear combination of selected features weighted by their respective coefficients. Receiver operator characteristic (ROC) curves were used to assess the performance of the radiomics signature in the two sets.

Construction of Radiomics Nomogram

Clinical factors included age, tumor size, tumor margin, and MRI-reported LN status. Risk factors were determined by univariate and multivariate logistic regressions. Collinearity was assessed by variance inflation factor (VIF). Likelihood ratio test with backward step-down selection was applied for logistic regression. A nomogram was established in the training set on the basis of multivariable logistic regression.

Assessment of Nomogram Performance

ROC curves were used to assess the predictive performance of the radiomics nomogram in the training set, and calibration curves were employed to evaluate the agreement between the observed and predicted results. Good agreement between the true state of ALN and the predicted probability based on radiomics nomogram was achieved when the calibration curves were close to the diagonal line. Hosmer–Lemeshow test was used to determine the goodness of fit of the radiomics nomogram.

Validation of Radiomics Nomogram

The radiomics nomogram was validated using the validation set with the same formula in the training set. ROC and calibration curves were used to assess the predictive performance of the radiomics nomogram.

Clinical Use

The clinical usefulness of the nomogram was assessed using decision curve analysis (DCA) in the validation set. The ROC curve was used to calculate the area under the ROC curve (AUC). However, ROC only considers the specificity and sensitivity of the method, and DCA determines the clinical practicability of radiomics nomograms by quantifying the net benefits under different threshold probabilities in the validation set. The calculation formula is as follows:

net benefit treated=TPnFPn(Pt1Pt),

where TP and FP are the true positive count and the false positive count, respectively; and n is the number of subjects; and Pt is the threshold probability.

Statistical Analysis

All statistical tests were conducted in R3.5.1. Chi-square or Fisher’s exact test was used to compare the differences in categorical variables, and a two-sample t test was applied to compare the differences in age and tumor size. LASSO logistic regression was used to select the most discriminating features and build the radiomics signature via 10-fold cross validation based on the minimum criteria. The radiomics signature was calculated by combining the features weighted by their coefficients. Clinical factors were used to construct the clinical model by using multivariable stepwise-backward logistic regression, and the clinical nomogram was provided. VIFs were accessed to exclude multi-collinearity, and the combined nomogram was built similarly to clinics, except for the combination of clinical factors and radiomics signature. ROC analysis, calibration curve, and DCA analysis were employed to evaluate the performance of the nomograms. DeLong’s test was used to compare the differences of ROC curves. In addition, “glmnet,” “glm,” “rms,” “pROC,” “Calibration Curves,” and “Decision Curve” packages were used. P < 0.05 indicates statistically significant difference.

Results

Patient Characteristics

Table 1 exhibits the patients’ characteristics in the training and validation sets. The kappa value obtained in agreement of observation is 0.85, indicating a good agreement between two observers in MRI-reported LN status classification. Molecular subtype was detected in both groups, which showed no significant differences between the metastatic and non-metastatic groups in terms of age, tumor size, ADC value, enhanced features and histological grade (p > 0.05). The proportions of ALN metastasis in the training and validation sets were 47.1% and 47.3%, respectively. These results justified their use as training and validation sets.

TABLE 1
www.frontiersin.org

Table 1 Patient characteristics in the training and validation sets.

Feature Selection, Radiomics Signature Building, and Validation

A total of 396 radiomics features were extracted from each MR image and divided into six groups, namely, histogram, form factor matrix, gray-level co-occurrence matrix (GLCM), gray-level size zone matrix (GLSZM), Haralick matrix, and run length matrix (RLM). The ICCs ranged from 0.863 to 0.982 and from 0.832 to 0.935 in the intra- and inter-observers, respectively. Three LN state-related features with non-zero coefficients, namely, GLCMEnergy_AllDirection_offset7, LargeAreaEmphasis, and Correlation_AllDirection_offset7_SD were selected from the LASSO model in the training set (Figures 2A, B). The calculation formula is as follows:

Radscore = -0.292+(0.501×Correlation_AllDirection_offset7_SD)+ (0.0551× GLCMEnergy_AllDirection_offset7)+(0.0821× LargeAreaEmphasis)
FIGURE 2
www.frontiersin.org

Figure 2 LASSO algorithm for radiomics feature selection and the predictive performance of the radiomics. (A) Mean square error path using 10-fold cross validation; (B) LASSO coefficient profiles of the radiomics features; and (C, D) ROC curves of the radiomics signature and nomogram in the training and validation sets.

Significant difference was observed in the radiomics scores between LN-negative and LN-positive patients in the two sets (P < 0.01). The radiomics signature yielded AUCs of 0.78 (95% CI, 0.73–0.83) in the training set and 0.79 (95% CI, 0.73–0.85) in the validation set (Figures 2C, D).

Development and Validation of the Radiomics Nomogram Model

The radiomics signature and the MRI-reported LN status were identified as risk factors of LN metastasis in breast cancer (Table 2). The MRI-reported LN status was a qualitative feature that could be easily obtained. No collinearity was observed because the VIF of the predictor ranged from 1.10 to 1.25. The nomogram model included the radiomics signature and the MRI-reported LN status (Figure 3A). In the calibration curve in Figures 3B, C, the gray line represents perfect prediction, and the dotted line represents the calibration curve of the radiomics nomogram. The calibration curve and the nonsignificant Hosmer–Lemeshow test showed good agreement between the true state of ALN and the predicted probability based on radiomics (P = 0.663). The radiomics nomogram yielded AUCs of 0.92 (95% CI, 0.87–0.97) in the training set and 0.90 (95% CI, 0.85–0.95) in the validation set (Figures 2C, D). Significant difference was observed between the differences of ROC curves in the two sets (P < 0.001). The nomogram model yielded an AUC of 0.79 (95% CI, 0.70–0.87) in the cN0 subgroup (Figure 4). The results of DCA are shown in Figure 5. When the threshold probability ranged from 0.1 to 1.0 in the validation set, the radiomics nomogram to predict LN metastasis provides more net benefit than the “treat all” or “treat none” scheme. Therefore, our nomogram excellently performed in discrimination, calibration, and clinical use.

TABLE 2
www.frontiersin.org

Table 2 Risk factors for ALN metastasis in breast cancer.

FIGURE 3
www.frontiersin.org

Figure 3 Radiomics nomogram with radiomics signature and LN status (A) and calibration curves of the radiomics nomogram in the training (B) and validation (C) sets. Calibration curves indicate that the predicted probability has a good agreement with the actual state of axillary lymph node.

FIGURE 4
www.frontiersin.org

Figure 4 ROC curves of the nomogram in the cN0 subgroup in the validation set.

FIGURE 5
www.frontiersin.org

Figure 5 DCA of the radiomics nomogram. The y axis represents the net benefits, while the x axis represents the threshold probability. The red line represents the radiomics nomogram. The blue line represents the assumption that all patients were included in the lymph node metastasis group. The black line represents the assumption that all patients were included in the non-lymph node metastasis group.

Discussion

LN metastasis is a negative prognostic factor of breast cancer (28, 29). Thus, non-invasive LN assessment tools are promising. In this study, a radiomics nomogram model based on MRI was developed to predict the pretreatment of ALN metastasis in breast cancer and was validated using an independent dataset. This nomogram model was composed of radiomics signature and MR-reported LN status with AUCs of 0.92 in the training set and 0.90 in the validation set. LN metastasis has been predicted on the basis of clinical information or radiomics features only (2, 17, 30, 31). This research combined clinical information with radiomics features and used visualization nomogram to predict LN metastasis.

LN status has a certain diagnostic performance in differentiating ALN metastases (22). In this study, the MRI-reported LN status remarkably differs between the metastatic and non-metastatic groups. Moreover, univariate and multivariate logistic regression models have identified the MRI-reported LN status as an independent predictor of ALN metastasis. Therefore, this status was used as a predictor of the model. Previous study (32) used dynamic gadopentetate dimeglumine (Gd) enhanced MRI to evaluate axilla status in patients with breast cancer, and used ROC curves to compare enhancement indices and nodal area with histopathology of excised nodes, with AUCs from 0.77 to 0.88. Our results showed that compared with conventional MRI, radiomics nomogram had higher AUC of 0.90.

The application of radiomics nomogram provides a new approach for establishing a LN metastasis prediction model by using multiple characteristics. We previously used CESM-based radiomics signature and CESM-reported LN status to construct a radiomics nomogram to predict axillary LN metastasis, yielding an AUC of 0.79 in external validation cohort (33). Qiu et al. used 21 texture features derived from ultrasound imaging and ultrasound-reported LN status to predict LN metastasis in breast cancers, with an AUC of 0.759 in validation set (34). In our present study, the proposed MRI-based radiomics nomogram showed better performance than CESM-based and ultrasound-based radiomics nomogram, which may be used as an individualized model to visualize the risk of ALN metastasis by doctors and patients, and may meet the requirements for the development of precision medicine (35).

Tan et al. (36) not only used radiomics signature to predict LN status but also incorporated molecular subtype and PR status in nomogram. Other previous studies also used clinic-pathological characteristics to establish models in predicting LN metastasis of breast cancer patients, such as lymphovascular invasion and serum miRNA expression (37, 38), which might have a limited clinical implication, because characteristics such as molecular subtype, PR status, lymphvascular invasion, and miRNA was usually obtained by biopsy or other examinations, which to some extent limited the clinical application of these prediction models. However, the proposed radiomics nomogram only incorporated the MR-reported LN status and radiomics signature, which could be obtained by a non-invasive way before surgery, with an acceptable performance in LN metastasis predicting.

The discrimination and calibration performance of radiomics nomograms does not represent their clinical usefulness. Thus, whether this technique could improve patient outcome was assessed using DCA. Within the threshold probability range of 0.1–1.0, the radiomics nomogram provided more net benefits than the “treat all” or “treat none” scheme.

The proposed nomogram model showed good discriminating performance in cN0 patients who are difficult to diagnose by using traditional methods.

This study offered other notable advantages. Prior to feature extraction, some preprocessing techniques were applied to improve feature discrimination, and ICCs were used to evaluate the reproducibility of the radiomics feature extraction. These methods improved the reliability of this study.

This study has several limitations. First, the patients were enrolled from a single institution with a limited number. Despite the promising prospect, a large sample size and a multicenter study are warranted to prove the robustness of the proposed nomogram. Second, image segmentation was conducted manually. Although ICCs exhibited good reproducibility in feature extraction, the automated method for image segmentation provides stability (39, 40). Third, the methodology was limited by its statistical robustness, which could be overcome only through the true-blinded testing of the hypothesis. Future studies should adopt a double-blinded prospective design. Fourth, this study was performed retrospectively. In the future, the authors aim to collaborate with surgical colleagues and develop a prospective study to validate the proposed nomogram. Finally, the radiomics features were not extracted from the LNs.

In summary, the radiomics nomogram combined with MRI-based radiomics and clinical risk factors exhibited good predictive performance, calibration, and clinical utility in identifying ALN metastasis in patients with breast cancer. MRI-based radiomics could serve as a potential tool to help clinicians generate optimal clinical decisions and avoid overtreatment for patients with breast cancer.

Data Availability Statement

All datasets generated for this study are included in the article/supplementary material.

Ethics Statement

The studies involving human participants were reviewed and approved by the Institutional Review Board of Yantai Yuhuangding Hospital, Affiliated Hospital of Qingdao University. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

NM and YD created the study design. HX, HM, and FL collected the data. NM and SD processed the data. NH and WZ conducted data analysis. NM and YD wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by the National Natural Science Foundation of China (82001775).

Conflict of Interest

SD was employed by the company GE Healthcare.

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.

Acknowledgments

We thank surgical oncologist Kun Zhang, pathologist Ping Yang, and Weiyi Chen.

References

1. Kootstra J, Hoekstra-Weebers JE, Rietman H, de Vries J, Baas P, Geertzen JH, et al. Quality of life after sentinel lymph node biopsy or axillary lymph node dissection in stage I/II breast cancer patients: a prospective longitudinal study. Ann Surg Oncol (2008) 15:2533–41. doi: 10.1245/s10434-008-9996-9

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Kvistad KA, Rydland J, Smethurst HB, Lundgren S, Fjosne HE, Haraldseth O. Axillary lymph node metastases in breast cancer: preoperative detection with dynamic contrast-enhanced MRI. Eur Radiol (2000) 10:1464–71. doi: 10.1007/s003300000370

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Zhao J, Zhang J, Zhu QL, Jiang YX, Sun Q, Zhou YD, et al. The value of contrast-enhanced ultrasound for sentinel lymph node identification and characterisation in pre-operative breast cancer patients: A prospective study. Eur Radiol (2018) 28:1654–61. doi: 10.1007/s00330-017-5089-0

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Fusco R, Sansone M, Granata V, Di Bonito M, Avino F, Catalano O, et al. Use of Quantitative Morphological and Functional Features for Assessment of Axillary Lymph Node in Breast Dynamic Contrast-Enhanced Magnetic Resonance Imaging. BioMed Res Int (2018) 2018:2610801. doi: 10.1155/2018/2610801

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Gillies RJ, Kinahan PE, Hricak H. Radiomics: Images Are More than Pictures, They Are Data. Radiology (2016) 278:563–77. doi: 10.1148/radiol.2015151169

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Lambin P, Rios-Velazquez E, Leijenaar R, Carvalho S, van Stiphout RG, Granton P, et al. Radiomics: extracting more information from medical images using advanced feature analysis. Eur J Cancer (2012) 48:441–6. doi: 10.1016/j.ejca.2011.11.036

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho S, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun (2014) 5:4006. doi: 10.1038/ncomms5006

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Mao N, Wang Q, Liu M, Dong J, Xiao C, Sun N, et al. Computerized Image Analysis to Differentiate Benign and Malignant Breast Tumors on Magnetic Resonance Diffusion Weighted Image: A Preliminary Study. J Comput Assist Tomogr (2019) 43:93–7. doi: 10.1097/rct.0000000000000793

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Mao N, Yin P, Wang Q, Liu M, Dong J, Zhang X, et al. Added Value of Radiomics on Mammography for Breast Cancer Diagnosis: A Feasibility Study. J Am Coll Radiol (2019) 16:485–91. doi: 10.1016/j.jacr.2018.09.041

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Yin P, Mao N, Liu X, Sun C, Wang S, Chen L, et al. Can clinical radiomics nomogram based on 3D multiparametric MRI features and clinical characteristics estimate early recurrence of pelvic chondrosarcoma? J Magn Reson Imaging (2020) 51:435–45. doi: 10.1002/jmri.26834

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Yin P, Mao N, Zhao C, Wu J, Sun C, Chen L, et al. Comparison of radiomics machine-learning classifiers and feature selection for differentiation of sacral chordoma and sacral giant cell tumour based on 3D computed tomography features. Eur Radiol (2019) 29:1841–7. doi: 10.1007/s00330-018-5730-6

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Birkhahn M, Mitra AP, Cote RJ. Molecular markers for bladder cancer: the road to a multimarker approach. Expert Rev Anticancer Ther (2007) 7:1717–27. doi: 10.1586/14737140.7.12.1717

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Vellinga TT, Kranenburg O, Frenkel N, Ubink I, Marvin D, Govaert K, et al. Lymphangiogenic Gene Expression Is Associated With Lymph Node Recurrence and Poor Prognosis After Partial Hepatectomy for Colorectal Liver Metastasis. Ann Surg (2017) 266:765–71. doi: 10.1097/sla.0000000000002430

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Gorelik E, Landsittel DP, Marrangoni AM, Modugno F, Velikokhatnaya L, Winans MT, et al. Multiplexed immunobead-based cytokine profiling for early detection of ovarian cancer. Cancer Epidemiol Biomarkers Prev (2005) 14:981–7. doi: 10.1158/1055-9965.epi-04-0404

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Paik S, Shak S, Tang G, Kim C, Baker J, Cronin M, et al. A multigene assay to predict recurrence of tamoxifen-treated, node-negative breast cancer. N Engl J Med (2004) 351:2817–26. doi: 10.1056/NEJMoa041588

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Sparano JA, Gray RJ, Makower DF, Pritchard KI, Albain KS, Hayes DF, et al. Prospective Validation of a 21-Gene Expression Assay in Breast Cancer. N Engl J Med (2015) 373:2005–14. doi: 10.1056/NEJMoa1510764

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Dihge L, Bendahl PO, Ryden L. Nomograms for preoperative prediction of axillary nodal status in breast cancer. Br J Surg (2017) 104:1494–505. doi: 10.1002/bjs.10583

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Huang YQ, Liang CH, He L, Tian J, Liang CS, Chen X, et al. Development and Validation of a Radiomics Nomogram for Preoperative Prediction of Lymph Node Metastasis in Colorectal Cancer. J Clin Oncol (2016) 34:2157–64. doi: 10.1200/jco.2015.65.9128

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Wu S, Zheng J, Li Y, Yu H, Shi S, Xie W, et al. A Radiomics Nomogram for the Preoperative Prediction of Lymph Node Metastasis in Bladder Cancer. Clin Cancer Res (2017) 23:6904–11. doi: 10.1158/1078-0432.ccr-17-1510

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Valente SA, Levine GM, Silverstein MJ, Rayhanabad JA, Weng-Grumley JG, Ji L, et al. Accuracy of predicting axillary lymph node positivity by physical examination, mammography, ultrasonography, and magnetic resonance imaging. Ann Surg Oncol (2012) 19:1825–30. doi: 10.1245/s10434-011-2200-7

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Mortellaro VE, Marshall J, Singer L, Hochwald SN, Chang M, Copeland EM, et al. Magnetic resonance imaging for axillary staging in patients with breast cancer. J Magn Reson Imaging (2009) 30:309–12. doi: 10.1002/jmri.21802

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Yoshimura G, Sakurai T, Oura S, Suzuma T, Tamaki T, Umemura T, et al. Evaluation of Axillary Lymph Node Status in Breast Cancer with MRI. Breast Cancer (1999) 6:249–58. doi: 10.1007/bf02967179

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Goldhirsch A, Winer EP, Coates AS, Gelber RD, Piccart-Gebhart M, Thürlimann B, et al. Personalizing the treatment of women with early breast cancer: highlights of the St Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2013. Ann Oncol (2013) 24:2206–23. doi: 10.1093/annonc/mdt303

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Collewet G, Strzelecki M, Mariette F. Influence of MRI acquisition protocols and image intensity normalization methods on texture classification. Magn Reson Imaging (2004) 22:81–91. doi: 10.1016/j.mri.2003.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Gibbs P, Turnbull LW. Textural analysis of contrast-enhanced MR images of the breast. Magn Reson Med (2003) 50:92–8. doi: 10.1002/mrm.10496

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Depeursinge A, Foncubierta-Rodriguez A, Van De Ville D, Muller H. Three-dimensional solid texture analysis in biomedical imaging: review and opportunities. Med Image Anal (2014) 18:176–96. doi: 10.1016/j.media.2013.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Sauerbrei W, Royston P, Binder H. Selection of important variables and determination of functional form for continuous predictors in multivariable model building. Stat Med (2007) 26:5512–28. doi: 10.1002/sim.3148

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Luini A, Gatti G, Ballardini B, Zurrida S, Galimberti V, Veronesi P, et al. Development of axillary surgery in breast cancer. Ann Oncol (2005) 16:259–62. doi: 10.1093/annonc/mdi060

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Cianfrocca M, Goldstein LJ. Prognostic and predictive factors in early-stage breast cancer. Oncologist (2004) 9:606–16. doi: 10.1634/theoncologist.9-6-606

PubMed Abstract | CrossRef Full Text | Google Scholar

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

31. Sodano C, Clauser P, Dietzel M, Kapetas P, Pinker K, Helbich TH, et al. Clinical relevance of total choline (tCho) quantification in suspicious lesions on multiparametric breast MRI. Eur Radiol (2020) 30(6):3371–82. doi: 10.1007/s00330-020-06678-z

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Murray AD, Staff RT, Redpath TW, Gilbert FJ, Ah-See AK, Brookes JA, et al. Dynamic contrast enhanced MRI of the axilla in women with breast cancer: comparison with pathology of excised nodes. Br J Radiol (2002) 75:220–8. doi: 10.1259/bjr.75.891.750220

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Mao N, Yin P, Li Q, Wang Q, Liu M, Ma H, et al. Radiomics nomogram of contrast-enhanced spectral mammography for prediction of axillary lymph node metastasis in breast cancer: a multicenter study. Eur Radiol (2020). doi: 10.1007/s00330-020-07016-z

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Qiu X, Jiang Y, Zhao Q, Yan C, Huang M, Jiang T. Could Ultrasound-Based Radiomics Noninvasively Predict Axillary Lymph Node Metastasis in Breast Cancer? J Ultrasound Med (2020) 39(10):1897–905. doi: 10.1002/jum.15294

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Balachandran VP, Gonen M, Smith JJ, DeMatteo RP. Nomograms in oncology: more than meets the eye. Lancet Oncol (2015) 16:e173–80. doi: 10.1016/s1470-2045(14)71116-7

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Tan H, Wu Y, Bao F, Zhou J, Wan J, Tian J, et al. Mammography-based radiomics nomogram: a potential biomarker to predict axillary lymph node metastasis in breast cancer. Br J Radiol (2020) 93:20191019. doi: 10.1259/bjr.20191019

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Bevilacqua JL, Kattan MW, Fey JV, Cody HS, Borgen PI, Van Zee KJ. Doctor, what are my chances of having a positive sentinel node? A validated nomogram for risk estimation. J Clin Oncol (2007) 25:3670–9. doi: 10.1200/jco.2006.08.8013

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Shiino S, Matsuzaki J, Shimomura A, Kawauchi J, Takizawa S, Sakamoto H, et al. Serum miRNA-based Prediction of Axillary Lymph Node Metastasis in Breast Cancer. Clin Cancer Res (2019) 25:1817–27. doi: 10.1158/1078-0432.Ccr-18-1414

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Jung SC, Choi SH, Yeom JA, Kim JH, Ryoo I, Kim SC, et al. Cerebral blood volume analysis in glioblastomas using dynamic susceptibility contrast-enhanced perfusion MRI: a comparison of manual and semiautomatic segmentation methods. PloS One (2013) 8:e69323. doi: 10.1371/journal.pone.0069323

PubMed Abstract | CrossRef Full Text | Google Scholar

40. de Hoop B, Gietema H, van Ginneken B, Zanen P, Groenewegen G, Prokop M. A comparison of six software packages for evaluation of solid lung nodules using semi-automated volumetry: what is the minimum increase in size to detect growth in repeated CT examinations. Eur Radiol (2009) 19:800–8. doi: 10.1007/s00330-008-1229-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: breast cancer, lymphatic metastasis, radiomics, nomogram, magnetic resonance imaging

Citation: Mao N, Dai Y, Lin F, Ma H, Duan S, Xie H, Zhao W and Hong N (2020) Radiomics Nomogram of DCE-MRI for the Prediction of Axillary Lymph Node Metastasis in Breast Cancer. Front. Oncol. 10:541849. doi: 10.3389/fonc.2020.541849

Received: 10 March 2020; Accepted: 29 September 2020;
Published: 27 October 2020.

Edited by:

Bo Gao, Affiliated Hospital of Guizhou Medical University, China

Reviewed by:

Ying-Shi Sun, Peking University Cancer Hospital, China
Haitao Yang, Chongqing Medical University, China

Copyright © 2020 Mao, Dai, Lin, Ma, Duan, Xie, Zhao and Hong. 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: Nan Hong, hongnan@bjmu.edu.cn; Wenlei Zhao, 271517561@qq.com

These authors have contributed equally to this work