Development and Validation of an MRI Radiomics-Based Signature to Predict Histological Grade in Patients with Invasive Breast Cancer

Abstract

Background

Histological grade is an important factor in the overall prognosis of patients with invasive breast cancer. Therefore, the non-invasive assessment of histological grade in breast cancer patients is an increasing concern for clinicians. We aimed to establish an MRI-based radiomics model for preoperative prediction of invasive breast cancer histological grade.

Methods

We enrolled 901 patients with invasive breast cancer and pre-operative MRI. Patients were randomly divided into the training cohort (n=630) and validation cohort (n=271) with a ratio of 7:3. A radiomics signature was constructed by extracting radiomics features from MRI images and was developed according to multivariate logistic regression analysis. The diagnostic performance of the radiomics model was assessed using receiver operating characteristic (ROC) curve analysis.

Results

Of the 901 patients, 618 (68.6%) were histological grade 1 or 2 while 283 (31.4%) were histological grade 3. The area under the ROC curve (AUC) of the radiomics model for the prediction of the histological grade was 0.761 (95% CI 0.728–0.794) in the training cohort and 0.722 (95% CI 0.664–0.777) in the validation cohort. The decision curve analysis (DCA) demonstrated the radiomics model’s clinical application value.

Conclusion

This is a preliminary study suggesting that the development of an MRI-based radiomics model can predict the histological grade of invasive breast cancer.

Keywords:

• breast cancer
• histologic grade
• magnetic resonance imaging
• radiomics
• signature

Introduction

Breast cancer is the most prevalent malignancy among women worldwide, accounting for almost 25% of all cancer cases in women. It is also the leading cause of cancer-related deaths in women worldwide.¹ Studies have shown that tumor histological grade is an independent prognostic factor in specific subgroups of breast cancer patients, including estrogen receptor positivity,² and lymph node metastasis.³ Invasive ductal carcinoma (IDC), the most common histological type of breast cancer, accounts for approximately 80% of all breast cancers. Therefore, accurate identification of tumor's histological grade in the IDC can provide a useful guide to the prognosis.

Magnetic resonance imaging (MRI) has been widely used in the study of breast cancer. Studies have shown that radiological analysis of MRI may help predict breast cancer subtypes,^4–6 neoadjuvant chemotherapy responses,^7,8 Ki-67 expression level,^9–11 pathological stage,¹² histological grade,^11,13,14 tumor malignancy,¹⁵ and breast cancer recurrence.^16,17 The eighth edition of the cancer staging system of the American Joint Committee on Cancer (AJCC), combined histological grade with anatomical staging to determine the clinical prognostic stage, where histological grade 3 luminal carcinoma is one stage higher than its anatomical stage.^14,18 Therefore, we wish to distinguish histological grade 3 tumors from grade 1 or 2 tumors. However, to our knowledge, we used the largest sample size with effective predictors to develop and validate an MRI-based radiomics model for preoperative prediction of invasive breast cancer histological grade.

Materials and Methods

Patients

Patients were collected from The Cancer Imaging Archive (TCIA; http://www.cancerimagingarchive.net/) datasets.¹⁹, 922 patients whose pre-operative MRI acquired from 1 January 2000 to 23 March 2014 with post-operative pathology confirmed invasive breast cancer. Patients with prior breast surgery, history of breast cancer, or neoadjuvant therapy before the MRI acquisition were excluded. 901 patients were eventually included in this study, due to the lack of histological grade information for 21 patients. Patients were randomly divided into the training cohort (n=630) and validation cohort (n=271) with a ratio of 7:3. The study flow diagram is shown in Figure 1.

Figure 1 The flow diagram of the study.

Imaging Acquisition, Segmentation, and Radiomics Feature Extraction

For the patients included in this study, cross-sectional breast MRI images were acquired in the prone position using a 1.5T or 3T scanner. Detailed scanner-related information and MR acquisition parameters are shown in the Tables S1 and S2, respectively. The following sequences were acquired by MRI: a non-fat-saturated T1-weighted sequence, a fat-saturated gradient-echo T1-weighted enhancement sequence, and four post-enhancement T1-weighted sequences typically acquired after intravenous contrast injection was used with a dose of 0.2 mL/kg body weight via peripheral veins.

A fuzzy C-means automatic segmentation²⁰ was used to obtain the tumor mask. The radiomics feature of breast cancer was automatically extracted from the N4-corrected²¹ T1-non-fat saturated (T1-NFS) images and first post-contrast sequences. A comprehensive set of 529 radiomics features were collected from TCIA datasets¹⁹ that have been shown to be effective predictors and could quantify characteristics of the breast, tumor, and fibroglandular tissue. The process of image segmentation and radiomics feature extraction can be available in the TCIA datasets publication.⁴

Construction of the Radiomics Model

To avoid overfitting the radiomics signature, features were further selected in three steps before the radiomics signature was constructed. First, features with ICC >0.75 within the training cohort were retained. Second, statistically, significant features were screened out using the univariate logistic analysis in the training cohort. Third, the most valuable features were selected using the least absolute shrinkage and selection operator (LASSO). A radiomics score (Rad-score) was calculated by using a formula based on the radiomics features.

The Rad-score was used to establish a radiomics signature through multivariate logistic regression. The diagnostic performance of the radiomics signature model in predicting histological grade was evaluated based on the area under the receiver operator characteristic (ROC) curve (AUC) in both the training cohort and testing cohort. To evaluate the clinical effectiveness of the radiomics signature, a decision curve analysis (DCA) was performed by calculating the net benefit of a threshold probability range across the training and testing cohorts.

Statistical Analysis

Statistical tests were performed using SPSS 25.0 software (IBM) and R programming language (ver. 3.5.1, http://www.r-project.org). We chose 26 features from the training cohort and tested them on the testing cohort to investigate the prognostic usefulness of imaging features.

Results

Clinical Features of the Patients

Nine hundred and one patients were enrolled in the study as shown in Figure 1. The clinical features of invasive breast cancer in the training and testing cohorts are shown in Table 1. Of the 901 patients, Histologic grade 1 was observed in 164 patients, grade 2 in 454 patients, and grade 3 in 283 patients, and 618 (68.6%) were histological grade 1 or 2 while 283 (31.4%) were histological grade 3.

Table 1 Clinical and Tumor Features in the Training and Validation Cohorts

Characteristics	Training Cohort (N = 618)	Validation Cohort (N = 283)
Median age (years)	53.28	49.86
Race
White	465	173
Black	108	93
Others*	36	8
Not available	9	9
Menopausal status
Pre	250	147
Post	356	132
Not available	12	4
Tumor staging (size)
T1	302	95
T2	246	145
T3	59	31
T4	8	12
Not available	3	0
ER status
Positive	536	137
Negative	82	146
PR status
Positive	486	103
Negative	132	180
HER2 status
Positive	89	66
Negative	529	217

Abbreviations: ER, estrogen receptor; PR, progesterone receptor; HER2, human epidermal growth factor; *Includes Asian, Native, Hispanic.

Construction of the Radiomics Model

A total of 529 radiomics features were extracted from the MRI images. After univariate correlation analysis, 143 radiomics features showed significant differences in predicting histological grade. These features were sequentially imported into LASSO (Figure 2) to obtain the most valuable features, resulting in 26 useful features. Finally, the radiomics signature was established by using the 26 features (Table 2).

Table 2 Radiomics Features Selection from the CTE in the Training Cohort

Radiomics Features	Coefficient
Volume_cu_mm_Tumor	−0.2952
Inf_mea_of_corr2_Tumor	−0.3053
Max_Enhancement_from_char_curv	−0.2096
BreastDensity_T1	−0.0853
Grouping_based_variance_of_washout_slope_3D_tumor_Group_1	0.9801
Grouping_based_variance_of_washout_slope_2D_tumorSlice_Group_2	−0.1994
Grouping_based_proportion_of_tumor_voxels_2D_tumorSlice_Group_1	−0.2213
Grouping_based_proportion_of_tumor_voxels_2D_tumorSlice_Group_2	0.0376
Grouping_based_variance_of_peak_enhancement_slope_3D_tissue_PostCon_Group_3	0.0850
Grouping_based_variance_of_washin_slope_3D_tissue_PostCon_Group_2	0.1022
Grouping_based_variance_of_washout_slope_3D_tissue_PostCon_Group_1	0.0532
Mean_norm_DHOG_max_timepoint_binsize_6_with_filling_Tumor	0.1349
SER_Washout_tumor_vol_cu_mm	−0.1489
SER_map_Correlation2_tumor	−0.0687
SER_map_Sum_of_Squares_variance_tumor	−0.1204
PE_map_Cluster_Prominence_tumor	0.1228
WashinRate_map_Correlation2_tumor	0.3493
WashinRate_map_Max_Probability_tumor	0.1759
WashinRate_map_information_measure_correlation1_tumor	−0.0017
WashinRate_map_skewness_tumor	−0.3948
SER_map_skewness_tissue_T1	0.2037
Peak_SER_tissue_PostCon	0.2229
SER_map_Cluster_Shade_tissue_PostCon	−0.3002
PE_map_kurtosis_tissue_PostCon	−0.2192
WashinRate_map_Autocorrelation_tissue_PostCon	0.1409
WashinRate_map_Cluster_Shade_tissue_PostCon	−0.0817

Figure 2 The LASSO regression model was used to select radiomics features. (A) LASSO coefficient profiles of the 26 radiomics features. A coefficient profile plot was generated versus the selected log (λ) value using ten-fold cross-validation, where optimal λ resulted in 8 features with nonzero coefficients. (B) The 26 radiomics features’ LASSO coefficient profiles. The log (λ) sequence was used to create a coefficient profile plot. Using 10-fold cross-validation, the dotted vertical line was drawn at the value chosen.

The radiomics model based on the training cohort and testing cohort is shown in Figure 3. Table 3 summarizes the discriminatory efficacies of the radiomics model. The AUC was 0.761 (95% CI 0.728–0.794) in the training cohort and 0.722 (95% CI 0.664–0.777) in the validation cohort for predicting histological grade. The DCA indicated that in most training and testing cohorts within reasonable threshold probabilities, the clinical radiomics nomogram added a greater overall net benefit in predicting histological grade in invasive breast cancer. The DCA for the radiomics model in the training and testing cohorts is shown in Figure 4.

Table 3 Performance of the Radiomics Model in the Training and Validation Cohorts

	Training Cohort (n=630)	Validation Cohort (n=271)
Radiomics model	AUC (95% CI) SEN SPE ACC	AUC (95% CI) SEN SPE ACC
	0.761 (0.728–0.794) 0.771 0.601 0.717	0.722 (0.664–0.777) 0.763 0.624 0.720

Abbreviations: AUC, area under the curve; SEN, sensitivity; SPE, specificity; ACC, accuracy; 95% CI, 95% confidence interval.

Figure 3 The receiver operator characteristic curves of the radiomics model. (A) The radiomics model of ROC curves in the training cohort. (B) The radiomics model of ROC curves in the validation cohort.

Figure 4 The decision curve analysis (DCA) of the radiomics model. (A) The radiomics model of DCA in the training cohort. (B) The radiomics model of DCA in the validation cohort.

Discussion

In the present study, we assessed the development and validation of a radiomics model to predict histological grade in invasive breast cancer. According to our findings, the radiomics signature could predict histological grade in both training and testing cohorts. At the same time, DCA demonstrated the radiomics model’s clinical application value.

Studies have shown that histological grade is a key factor in the overall prognosis of patients.^22–25 Grading of breast tumors requires histological evaluation of tissue samples obtained by guided biopsy or surgical excision. Breast cancers are graded according to the Scarff-Bloom-Richardson grading system.²⁶ This grading system relies on the frequency of cell mitosis (cytokinesis rate), Tubule formation (the percentage of tumors containing tubular structures), and the grade of nuclear pleomorphism.²⁷

Saha et al⁴ analyzed a total of 922 individuals with invasive breast cancer and pre-operative MRI to predict the molecular, genomic, and proliferation characteristics. A computer algorithm was used to extract 529 features of the tumor and surrounding tissue from the MRI images. Those features could be found in the published literature as well as those developed in their research.

Our patients and radiomics features were collected from this study. Song et al¹⁴ used multiparametric MRI to create machine learning-based prediction models that performed well in predicting Ki-67 and histologic grade in patients with luminal breast tumors. Scholars²⁸ reported that a total of 205 patients with contrast-enhanced spectral mammography examination invasive breast cancer were retrospectively enrolled and the radiomics model is a non-invasive predictive method that showed good application prospects in predicting histological grade. It is reported that the use of a combination of dynamic contrast-enhanced magnetic resonance imaging and T2-weighted imaging radiomic features could be used to predict histological grade in ductal breast carcinoma.²⁹ Although these studies had good diagnostic efficacy, the number of cases is relatively small. We enrolled 901 patients in this study and developed an MRI radiomics model with effective predictors to predict histological grade in patients with invasive breast cancer.

A total of 529 radiomics features extracted from the MRI images were enrolled in our study. After univariate correlation analysis, 143 radiomics features showed significant differences in predicting histological grade. These features were sequentially imported into LASSO to obtain the most valuable features, resulting in 26 useful features. Finally, the radiomics signature was established by using the 26 features.

Our study had several limitations. First, our 901 breast cancer patients were acquired from different machines. The variability of imaging acquisition parameters is a source of noise and future studies should focus on uniform imaging parameters for scanning. Second, this was a retrospective study. The multicenter, large-sample, and prospective cases need to be included in future studies.

We developed an MRI-based radiomics model to predict histological grade in invasive breast cancer. This radiomics signature can provide a preliminary method to predict histological grade of invasive breast cancer.

Ethical Approval and Consent to Participate

The study protocol was approved by the institutional review board of the First Affiliated Hospital of Wannan Medical College. Written informed consent was not required due to the database is publicly available.¹⁹ The personal information of patients was strictly protected. The study was carried out in accordance with the tenets of the Declaration of Helsinki 1964 and its later amendments.

Disclosure

The authors declare that they have no conflicts of interest in relation to this work.

Acknowledgments

We would like to thank the participants who contributed to our study.

Additional information

Funding

There is no funding to report.