Radiomic analysis reveals DCE-MRI features for prediction of molecular subtypes of breast cancer

Abstract

The purpose of this study was to investigate the role of features derived from breast dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) and to incorporated clinical information to predict the molecular subtypes of breast cancer. In particular, 60 breast cancers with the following four molecular subtypes were analyzed: luminal A, luminal B, human epidermal growth factor receptor-2 (HER2)-over-expressing and basal-like. The breast region was segmented and the suspicious tumor was depicted on sequentially scanned MR images from each case. In total, 90 features were obtained, including 88 imaging features related to morphology and texture as well as dynamic features from tumor and background parenchymal enhancement (BPE) and 2 clinical information-based parameters, namely, age and menopausal status. An evolutionary algorithm was used to select an optimal subset of features for classification. Using these features, we trained a multi-class logistic regression classifier that calculated the area under the receiver operating characteristic curve (AUC). The results of a prediction model using 24 selected features showed high overall classification performance, with an AUC value of 0.869. The predictive model discriminated among the luminal A, luminal B, HER2 and basal-like subtypes, with AUC values of 0.867, 0.786, 0.888 and 0.923, respectively. An additional independent dataset with 36 patients was utilized to validate the results. A similar classification analysis of the validation dataset showed an AUC of 0.872 using 15 image features, 10 of which were identical to those from the first cohort. We identified clinical information and 3D imaging features from DCE-MRI as candidate biomarkers for discriminating among four molecular subtypes of breast cancer.

About the Authors:

Ming Fan

Affiliation: Institute of Biomedical Engineering and Instrumentation, Hangzhou Dianzi University, Hangzhou, China

Hui Li

Affiliation: Institute of Biomedical Engineering and Instrumentation, Hangzhou Dianzi University, Hangzhou, China

Shijian Wang

Affiliation: Institute of Biomedical Engineering and Instrumentation, Hangzhou Dianzi University, Hangzhou, China

Bin Zheng

Affiliations Institute of Biomedical Engineering and Instrumentation, Hangzhou Dianzi University, Hangzhou, China, School of Electrical and Computer Engineering, University of Oklahoma, Norman, Oklahoma, United States of America

Juan Zhang

* E-mail: [email protected] (JZ); [email protected] (LL)

Affiliation: Zhejiang Cancer Hospital, Zhejiang Hangzhou, China

Lihua Li

* E-mail: [email protected] (JZ); [email protected] (LL)

Affiliation: Institute of Biomedical Engineering and Instrumentation, Hangzhou Dianzi University, Hangzhou, China

Introduction

Breast cancer is the most common disease that affects women globally. Breast cancer subtypes, which are molecular indicators that are indicative of outcome, can be used to guide targeted therapy. A classification system categorizes invasive breast carcinomas into four distinct molecular subtypes, luminal A, luminal B, human epidermal growth factor receptor-2 (HER2)-over-expressing and basal-like, which provide a clinically useful basis for therapeutic selection[1]. The various molecular subtypes have distinct recurrence and survival rates. In particular, the luminal A subtype shows the best survival among the four molecular subtypes [2], while the worst survival is observed for the HER2 and basal-like subtypes [3]. Patients with the luminal B subtype exhibit higher proliferation and have a poorer prognosis than those with the luminal A subtype [4]. Moreover, the basal-like and HER2 subtypes of breast cancer are more sensitive to preoperative chemotherapy [5,6] and are more likely to undergo a pathologic complete response [7] but have a worse prognosis than luminal and normal-like cancers[8,9]. Finally, the basal-like molecular subtype is significantly associated with nodal involvement, which is important for the prognosis of early-stage breast cancer[10].

Molecular subtyping is thus beneficial for the diagnosis and individualized treatment of breast cancer. However, the determination of subtypes by genetic analysis is invasive and expensive, requiring specialized equipment and technical expertise. Immunohistochemical surrogate biomarkers of estrogen receptor (ER), progesterone receptor (PR), and HER2 status have been used to define molecular subtypes. Therefore, demand exists for alternative means of classifying breast cancers into distinct molecular subtypes. Magnetic resonance imaging (MRI) has played an evolving role in the diagnosis and treatment of breast cancer, and studies have highlighted the value of dynamic contrast-enhanced MRI (DCE-MRI) in reflecting the anatomic and functional properties of tumors and facilitating treatment [11–13]. Imaging surrogates are also important for developing personalized treatment and determining tumor heterogeneity when assessing cancer evolution [14,15]. The morphologic and kinetic features extracted from DCE-MRI, including tumor shape, size, spiculation and contrast enhancement, have been shown to be associated with different histological types, tumor grades, and microvessel distributions [11]. In our previous work, we demonstrated that bilateral asymmetry of kinetic features between the left and right breasts [16] and dynamic features in non-tumor breast stroma [17,18] are highly accurate in discriminating benign from malignant tumors.

Over the past few years, pilot studies have attempted to correlate the tumor features extracted from MRI with the molecular subtypes of breast cancer [19,20]; this filed is referred to as “radiogenomics”, while the efforts to discover and use quantitative features are termed “radiomics” [19–21]. Detailed reviews of radiogenomics and radiomics have recently been published elsewhere [22–24]. Many studies have analyzed the associations between qualitative and semi-quantitative MRI features and molecular subtypes [14,19,25] or triple-negative breast cancers [27–29]. Koo et al correlated DCE MRI perfusion parameters (K^trans, k_ep and v_e) with prognostic factors and the triple negative subtype of breast cancer [30]. Semi-automatically extracted imaging features have also been computed to determine associations between quantitative MRI features and breast cancer subtypes [26,31,32]. The results demonstrated a moderate correlation between two dynamic imaging features and the luminal A and B subtypes of breast cancer [32]. An earlier study of 48 breast cancer patients also identified imaging features associated with the luminal B subtype [31]. That study showed the enhancement dynamics of the tumor along with enhancement of fibroglandular breast tissue, namely background parenchymal enhancement (BPE). A recent study identified 44 MRI features as well as 6 clinical and pathologic features per tumor for classifying molecular subtypes [33]. However, the predictive value of 3D volumetric features from DCE-MRI for the prediction of the four molecular subtypes of breast cancer has not been fully examined.

The purpose of the present study was to investigate the use of features extracted from DCE-MRI and clinical information for the prediction of the molecular subtypes of breast cancer. To our knowledge, our study is the first to classify the four distinct molecular subtypes of breast cancer using 3D volumetric imaging features. In particular, we performed high-throughput extraction of a large number of quantitative and minable imaging features, assuming that the subsets of these features convey prognostic and predictive information. Based on our previous work, we extract imaging features, including the shape, texture, statistics, dynamics and bilateral asymmetry of both the breast lesion and the non-tumor background parenchymal region. To assess the stability of these features for prediction, an additional dataset was used for validation. The successful application of this technique would allow MRI examinations to provide valuable information related to the prognosis of breast cancer.

Materials and methods

Patient population

This is a retrospective study approved by the Internal Research Review and Ethical Committee of the Zhejiang Cancer Hospital. All patient information was anonymized before the data were stored in our research dataset. In this study, two cohorts were included to assess and validate the prediction results using the same imaging protocols. For the first cohort, retrospective data from 116 pretreatment breast MRI evaluations performed at Zhejiang Cancer Hospital between January and July 2011 were collected. The first cohort includes 116 pretreatment breast DCE-MRI evaluations acquired from 116 previously underwent breast cancer diagnosis and treatment at this hospital. Forty patients with no pathologic examination or incomplete pathology data and those without a complete MRI sequence were excluded. Patients who underwent breast cancer treatment (e.g., chemotherapy or radiation therapy) before MRI were also excluded (n = 16). Ultimately, 60 patients were analyzed.

A reproducibility cohort of patients was also studied to assess variability. We collected 57 pretreatment breast MRI evaluations, excluding patients who did not receive pathologic examination or had incomplete pathology data (n = 19), as well as those who underwent chemotherapy or radiation therapy before MRI (n = 2). The final dataset included 36 breast cancer patients aged from 29 to 59 years, with a mean age of 47.06 years.

A pathology report was collected from the initial breast biopsy or surgery, and the ER, PR, and HER2 status of each patient with invasive cancer was determined using streptavidin-peroxidase (SP) for immunohistochemistry (IHC) analyses. The ER and/or PR status was determined as positive when at least 1% of the tumor cell nuclei showed staining for the ER or PR [34]. A sample was scored as HER2+ according to the American Society of Clinical Oncology (ASCO)/College of American Pathologists (CAP) guideline recommendations for HER2 testing in breast cancer [35]. The HER2 status was specifically considered positive when the IHC staining intensity score was greater than or equal to three. A HER2 score of 2+, as determined by IHC, with confirmation of gene amplification by fluorescence in situ hybridization (FISH), was also considered to indicate a positive HER2 status. Based on the receptor status, the molecular subtype classification was determined using criteria reported by prior investigators. Four molecular subtypes have specifically been determined: luminal A: ER and/or PR positive, HER2 negative; luminal B: ER and/or PR positive, HER2 positive; HER2: ER and PR negative, HER2 positive; and basal-like: ER, PR, and HER2 negative. Luminal breast cancers were assigned to the luminal B classification if the expression of Ki-67 was at least 14% [1]. The characteristics of the patients in the reproducibility cohort are shown in S1 Table. The statistical test showed no significant differences in clinical characteristics between the subtypes.

MR image acquisition

The same DCE-MRI scanning method was applied to each woman in our cohort following the current DCE-MRI examination protocol adopted in the Radiology Department of Zhejiang Cancer Hospital. All patients were scanned in the prone position using a Siemens MRI scanner (MAGNETOM Espree-Pink 1.5 T) and a dedicated eight-channel double-breast coil. A precontrast (baseline) series of fat-suppressed T1-weighted and fat-saturated T2-weighted 3D image scans was first acquired, followed by a DCE series after the intravenous injection of a contrast agent at a dose of 0.2 mmol per kilogram body weight and a saline flush of 20 ml at the same flow rate of 2 ml/s. Three postcontrast series (S-0, S-1, S-2 and S-3) of 3D MRI scans and data acquisitions were sequentially performed at time intervals of 1.5 min, 2.5 min and 8.5 min. The parameters were as follows: repetition time (TR)/echo time (TE), 4.4/1.6 ms; flip angle, 12°; acquisition matrix, 512 × 512; slice thickness, 1.2 mm in three MRI scan series (S-0 to S-2) and 0.8 mm in S-3; and spatial resolution, 0.625 × 0.625 × 1.2 mm³ in S-0, S-1 and S-2 and 0.625 × 0.625 × 1.2 mm³ in S-3. The first three image scan series all had the same number of image slices (88) and the same in-depth resolution, whereas the third series (S-3), with 160 image slices, showed high in-depth resolution. Therefore, only the images acquired from the first three MRI scan series (S-0, S-1, and S-2) were selected and used in the imaging feature computation and analysis.

Computer-based DCE-MRI analysis

Two radiologists with 20 years of experience first annotated the center location of the suspicious breast tumor in each case. The marked tumor center location was initially used as an initial segmentation “seed”, and a volumetric breast tumor boundary contour was automatically segmented [36]. In this method, the segmentation was first conducted with a spatial fuzzy C-means (FCM) algorithm and was then refined using a Markov random field (MRF)-based approach, the parameters of which were adaptively adjusted using segmentation results from contiguous slices. The visual examination and manual correction (if needed) of the segmentation results were performed by our investigators. In this study, manual correction of tumor boundary segmentation was necessary in less than 20% of cases. A computerized scheme was also applied to automatically segment the breast region depicted on all breast MR image slices and to, register sequential MRI scans (S-0, S-1, and S-2) following the same procedure as in our previous study [17]. In brief, the following steps were used for image registration: 1) performing initial image filtering using a Sobel filter, 2) searching for a maximum information window on each image slice, 3) performing an information matching process based on the correlation coefficients of two matched windows, and 4) registering S-1 and S-2 images with the S-0 image through a linearly shifting process. Fig 1 shows an example of CAD-segmented breast regions depicted on the DCE-MRI image slice before contrast agent injection (S-0). The left and right breast areas depicted on each breast MR image slice were segmented to exclude the voxels located behind the chest wall, and a breast tumor was also identified. From the segmented DCE-MRI, we identified the BPE in the abnormal breast exhibiting a breast tumor and the contralateral normal breast in cancer patients.

[Figure omitted. See PDF.]

Fig 1. An example of a breast DCE-MRI slice.

(a) and its CAD-segmented breast regions depicted on the image slice before contrast agent injection (S-0) in which the non-breast regions behind the chest wall are removed.

https://doi.org/10.1371/journal.pone.0171683.g001

From the registered MRI sequential scans, we applied a computerized scheme to extract 88 imaging features. The computerized scheme calculated 5 morphologic features, 45 texture features, 15 first-order statistical features, 29 dynamic features from the breast lesion and BPE, and 9 bilateral differences in background parenchymal area. Table 1 summarizes these DCE-MRI features and clinical information.

[Figure omitted. See PDF.]

Table 1. Summary of extracted DCE-MR image features.

https://doi.org/10.1371/journal.pone.0171683.t001

Morphologic features.

To assess tumor surface irregularity, the intra-tumor variability was measured for each tumor. For this purpose, the radial length was calculated for the distance between a given point on the surface of the tumor and the tumor centroid, and the mean and variance of the radial length were determined. The roughness, compactness, diameter and volume of the breast lesion were also obtained.

First-order statistics.

To evaluate tumor heterogeneity in terms of the data distribution, histograms were computed using the voxel intensities inside the 3D lesion from which statistical measures of the shape of the histograms were derived. These measures included 1) kurtosis, 2) skewness, 3) variance, 4) minimum, 5) maximum and 6) average values.

Texture features.

The second-order texture features were assessed using the gray level co-occurrence matrix (GLCM) based on 3D analysis to reflect tumor spatial complexity [37]. Volumetric texture features such as contrast, correlation, energy, homogeneity and entropy were included for both the lesion and the background parenchyma.

Dynamic features.

Dynamic features were evaluated to assess histopathology and angiogenesis in breast cancer. These features were calculated from the relative average difference in all registered pixel values between image series using the following equation:where S represents the regions of interest, including the lesion area and background parenchyma in the normal or abnormal breast; I_T(i) represents the value of the ith-matched pixel voxel in the Tth image slice; and N_j is the total number of voxels detected inside the 3D volume. The statistical features for the dynamic enhancement of in the breast lesion and background parenchyma were also calculated, including the variances, means, and maximum and minimum values. The background parenchymal features included features that were in the breast lesion and in lesion-free breast areas [16]. The computerized scheme sorted and ranked voxels based on the computed BPE values. Experiments and data analysis were conducted to determine the optimal threshold (p%) to select the percentage of voxels in the top BPE value ranking list to maximize the diagnostic value of using each BPE feature[16]. The contrast enhancement features in the background parenchyma region were computed from the top 1% of voxels in the breast exhibiting a breast lesion and the other, negative (lesion-free) breast.

Bilateral difference in BPE.

To evaluate breast tissue asymmetry, an important image phenotype that is highly related to the abnormal biological processes leading to cancer development[38], the bilateral asymmetry (p% = 1%) of characteristic kinetic features was also assessed between the abnormal breast exhibiting a breast tumor and the contralateral normal breast. These features were obtained by calculating the differences in the maximum, mean, and variances of between the left and right breasts.

Statistical analysis

All analyses were performed using a publicly available data mining and machine-learning software platform, Weka [39], to conduct feature selection and validation tasks. We employed a multi-class version of the logistic regression implemented in Weka for discrimination of the four molecular subtypes through a one-against-others strategy. Clinical information, namely, age and menopausal status, was also included in the prediction model. The Kruskal—Wallis test was used to test the statistical significance of features selected across the four subtypes. An evolutionary algorithm (EA)-based optimization method, i.e., “EvolutionarySearch”, was used to search for optimal feature subsets for classification. The mutation probability and cross over probability were set at 0.01 and 0.6, respectively. The EA chromosome population in each generation was 500, and the maximum number of generations was 200. To evaluate the performance of the classifier, a receiver-operating characteristic (ROC) analysis was performed, and the area under the ROC curve (AUC) was computed. The EA chromosome that achieved the highest AUC was selected to establish the optimal feature pool and build the optimal classifier.

To avoid overfitting of classifiers to the available data, the leave-one-out cross-validation (LOOCV) test was employed. Specifically, at each iteration of the LOOCV process, one sample was used for testing, and the other sample was used for training. This procedure was repeated for each sample. In each evaluation, we performed a Wrapper Subset Evaluator (WSE) feature selection using all training cases to search for optimal features from the optimal feature subset pool generated from EA. Using these features, a multi-class logistic classifier was trained on the basis of the training set and was tested on the one independent left-out testing case. The importance of the predictive imaging features in the classifier was evaluated by the selection frequencies of features over all of the LOOCV loops. The sensitivity and specificity were calculated using a threshold criterion that was determined as the value that would maximize the average sensitivity and specificity using the best-threshold function in R software, version 3.3.0 (R Development Core Team 2015).

Results

Molecular subtype classification in the first cohort

The characteristics of the included patients in the first cohort are shown in Table 2. Sixty women with breast cancer who fulfilled the selection criteria were included in the study sample. All of the patients were of Han Chinese origin, and the median age was 48 years (range: 32–78 years). Among the patients, 35 were premenopausal, and two had a family history of breast cancer. The distribution of the invasive breast cancers by tumor type was as follows: invasive ductal carcinoma (n = 51), intraductal carcinoma (n = 4), mucinous carcinoma (n = 2), carcinoma spongiosum (n = 2) and poorly differentiated adenocarcinoma (n = 1). The subjects included 34 luminal A, 8 luminal B, 7 HER2 and 11 basal-like breast cancer patients. Differences in the categorical variables (menopausal status, family history and tumor type) between distinct molecular subtypes were examined using Fisher’s exact test, and group differences in continuous measures were analyzed using a parametric two-sample independent t-test. The statistical tests showed no significant associations between the molecular subtypes and the pathologic conditions.

[Figure omitted. See PDF.]

Table 2. Characteristics of the molecular subtypes of patients.

https://doi.org/10.1371/journal.pone.0171683.t002

To evaluate the effectiveness of DCE-MRI features for the prediction of the molecular subtypes of breast cancer, molecular status was classified using the logistic regression model. The multi-class classifier based on 24 selected features produced an overall AUC value of 0.869. For the luminal A subtype, the classifier using all features computed an AUC value of 0.867 with a 95% confidence interval (CI) of [0.826, 0.908] and an 88.2% classification sensitivity level at 76.9% specificity. Regarding the luminal B subtype, the computed AUC was 0.786 with a 95% CI of [0.585, 0.987], and the optimal sensitivity and specificity were 86.5% and 62.5%, respectively. For the HER2 molecular subtype, the classifier had an AUC value of 0.888±0.083, a sensitivity of 81.1% and a specificity of 100%. Finally, for the basal-like subtype, the computed AUC, sensitivity and specificity were 0.923±0.037, 81.1% and 100%, respectively.

The Kruskal—Wallis test was performed to evaluate differences in image features among the four subtypes. However, no significant differences among the subtypes were revealed using the Kruskal—Wallis test with a p<0.05. We only observed a trend for an association (p = 0.06) of enhancement between S-2 and S-1 images from the contralateral normal breast. Higher values for this feature were observed in the HER2 subtype than the other subtypes. The importance of each feature representing the frequencies of each feature during the LOOCV loops is shown in Table 3. The top seven features with greater than 50% frequencies in the prediction model included two first-order statistical features (skewness and kurtosis), dynamic features in BPE and in breast lesions, two bilateral asymmetry features and one clinical factor (age). For example, Fig 2 shows that luminal A tumors had lower values of skewness and kurtosis features compared with luminal B tumors.

[Figure omitted. See PDF.]

Fig 2. Representative imaging features differentiate tumors with different molecular subtypes.

A) Postcontrast MRI of luminal A and luminal B breast cancers. B) Postcontrast MRI of luminal A and luminal B breast cancers and corresponding kurtosis density using a kernel smoothing function.

https://doi.org/10.1371/journal.pone.0171683.g002

[Figure omitted. See PDF.]

Table 3. List of 24 features that were selected in leave-one-out training and testing cycles to test 60 cases.

https://doi.org/10.1371/journal.pone.0171683.t003

Fig 3 shows and compares ROC curves representing the classification performance levels of applying the 15 dynamic features and the 8 combined texture and first-order statistical features to classify the luminal A, luminal B, HER2 and basal-like tumor molecular subtypes. For the luminal A subtype, the AUC value using the 24 features was significantly increased compared with using dynamic features (p = 0.003) or combined texture and statistical features (p = 0.033). For the luminal B subtype, using the combination of features, the computed AUC was significantly increased to 0.786 compared with using dynamic features (p = 0.04). Regarding the HER2 molecular subtype, the AUC value using all features was 0.888±0.083, which was significantly higher than that using dynamic features (p = 0.001) or combined texture and statistical features (p = 0.045). Finally, for the basal-like subtype, the computed AUC using all features was significantly increased to 0.923 compared with using dynamic features (p = 0.003) or combined texture and statistical features (p = 0.0003). The results indicate that the radiomic features had higher discriminative power in predicting molecular subtypes of breast cancer compared with using only dynamic, texture or statistical features.

[Figure omitted. See PDF.]

Fig 3. ROC curve in the first cohort for classifying the four molecular subtypes of breast cancer.

The classifiers based on dynamic features, morphologic features, first-order statistical features and clinical information are shown. Features are combined to classify between (a) luminal A and non-luminal A tumors; (b) luminal B and non-luminal B tumors; (c) HER2-positive and non-HER2-positive tumors; and (d) basal-like and non-basal-like tumors.

https://doi.org/10.1371/journal.pone.0171683.g003

Molecular subtype classification in the reproducibility cohort

Data from the 36 patients in the reproducibility cohort are displayed in S1 Table. The validation dataset included equal numbers (n = 9) of the four subtypes of patients. Among these patients, there were 34 invasive ductal carcinomas, 1 intraductal carcinoma and 1 poorly differentiated adenocarcinoma. The statistical tests showed no significant associations between the molecular subtypes and the pathologic conditions.

The same procedures of feature selection and classification were also performed for the validation cohort. The multi-class classifier based on the 15 selected features yielded an overall AUC value of 0.872. For the luminal A subtype, the classifier using all features computed an AUC value of 0.905 with a 95% CI of [0.7753, 0.9825], a sensitivity of 87.5% and a specificity of 92.9%.. Regarding the luminal B subtype, the computed AUC was 0.835 with a 95% CI of [0.7159, 0.9755], and a 59.3% classification sensitivity level at 100.0% specificity, respectively. For the HER2 molecular subtype, the classifier had an AUC value of 0.947±0.053, and the optimal sensitivity and specificity were 88.9% and 88.9%, respectively. Finally, for the basal like subtype, the computed AUC, sensitivity and specificity were 0.802±0.172, 88.9% and 63.0%, respectively. S1 Fig shows the classification performance by applying 15 image features to classify the molecular subtypes of luminal A, luminal B, HER2 and basal-like tumors, respectively. The Kruskal—Wallis test showed a statistically significant difference in bilateral differences (F83) between subtypes with p = 0.007, as shown in S2 Fig. A higher value for this feature was observed in the basal-like subtype than other subtypes. The importance of each feature representing the frequencies of each feature during the LOOCV loops is shown in Table 4. Among the 15 features, 10 (bold font in Table 4) were identical to those from the first cohort, including 3 texture, 4 BPE, 1 bilateral asymmetry and 2 dynamic features in breast lesions. When the features selected in the first cohort was utilized, multi-class classifier generated an overall AUC value of 0.630 for classifying the four subtypes. After feature selection, the classifier yielded an overall AUC of 0.730 and AUCs of 0.607, 0.807, 0.835, 0.671 for luminal A, luminal B, HER2 and basal-like subtypes, respectively.

[Figure omitted. See PDF.]

Table 4. List of 15 image features that were selected in leave-one-out training and testing cycles to test 36 cases in the validation dataset.

https://doi.org/10.1371/journal.pone.0171683.t004

Discussion

In this paper, we have presented results showing that features derived from DCE-MRI can be viewed as potential biomarkers for predicting the four molecular subtypes of breast cancer. Radiomic analyses were performed, and the features of morphologic, first-order statistic, texture, bilateral asymmetry and dynamic features in both breast lesions and BPE were extracted and combined for the prediction of molecular subtypes of breast cancer. A multi-class classifier was trained and tested in all the cases, and the AUCs were 0.867, 0.786, 0.888 and 0.923 for the prediction of the luminal A, luminal B, HER2 and basal-like subtypes, respectively, using 24 features. An additional independent validation cohort was also included in our study, and this cohort showed a similar performance with an overall AUC of 0.872. To our knowledge, our study is the first to identify full 3D volumetric DCE-MRI features for the prediction of the four molecular subtypes of breast cancer.

We observed imaging features that predict the distinct molecular subtypes of patients. In our experiments, low kurtosis and skewness were observed for the luminal A subtype, and these features exhibited high importance in the prediction model, which is in line with a previous study that highlighted skewness as a predictor for discriminating molecular subtypes of breast cancer[33]. Kurtosis and skewness features have been identified as biomarkers of tumor heterogeneity [40], and higher values of these biomarkers have been associated with treatment failure [41], while lower values indicate a response to treatment[42]. These findings were consistent with the clinical characteristics of the luminal A subtype, which is less aggressive and has the best prognosis, a high survival rate and a low recurrence rate compared to other subtypes[43–45]. We also observed that the features of dynamic enhancement of the breast lesion and BPE and bilateral asymmetry were among the features of the model. The differences in BPE between the normal and abnormal breasts may reflect the aggressiveness of the breast tumor, which is one of the main features of the HER2 subtype. We also observed that among all the subtypes, HER2 tumors had the highest enhancement values in the normal breasts, while luminal A and luminal B tumors had the lowest values. Compared to the luminal subtypes, the typically more aggressive HER2 subtype is associated with increased neoangiogenesis, which may be responsible for the higher enhancement, even in the normal breast.

Related studies by Mazurowski et al. evaluated 23 imaging features, including lesion, texture and dynamic features [31]. They found that the enhancement dynamics of the tumor and background parenchyma were correlated with the luminal B subtype in a 48-patient cohort. A subsequent study from the same group found associations between DCE-MRI features and the luminal A and B subtypes [32]. Similar to previous studies [31,32], our study also identified imaging biomarkers of BPE for prediction. However, these studies did not include the features of skewness, kurtosis and dynamic bilateral asymmetry, which were evaluated as important features for the prediction of subtypes in our study. The predictive nature of skewness and kurtosis is consistent with recently published research for the prediction of three molecular subtypes, namely, luminal-like, HER2-over-expressing and basal-like [33]. The authors examined morphological, histogram-based first-order, and GLCM-based texture features; they obtained an overall accuracy of 71.2% for prediction using imaging features, which was increased to 83.4% when nuclear and histologic grades were incorporated into the prediction model. However, the authors did not examine dynamic or BPE features for prediction, which have been shown as important image biomarkers in prediction. They employed feature evaluation and selection using the ratio of between-group to within-group sums of squares (BW-ratio) test, while our study used a global optimization scheme for radiogenomic feature selection with an AUC value as the fitness criterion. Moreover, the above studies acquired images with 1.5 T and 3.0 T MR systems with different protocols, which may have affected the imaging features and the association results.

In this study, an additional validation cohort was also included to assess variability. The prediction performance of this cohort was similar to that of the first cohort in terms of AUC value, and most of the selected features in the prediction model were the same (10 out of 15). Moreover, the classifier from the validation dataset used similar or identical features as the first dataset, such as skewness, bilateral asymmetric features, and BPE features, which indicates high reproducibility for these image features. However, the features most frequently selected in leave-one-out selection in the first cohort were not among the most selected features in the validation dataset. Therefore, further analyses with more samples should be conducted to validate the radiomic features for classifying subtypes of breast cancer.

This is a preliminary study with several limitations. For example, the statistical power of the study was limited by the relatively small sample size, which was less than in other radiogenomics studies with hundreds of features [46]. Other limitations are that only age and menopausal status were included as clinical information-based parameters in our prediction model; features such as nuclear grade, histologic grade and lymph node status were not used because relevant data were not available for the cohort. Another limitation is that we did not include semi-quantitative or pharmacokinetic analyses of signal intensity time-course data, which is still an important clinical MRI technique of choice for breast cancer diagnosis [47,48]. Despite these limitations, our results are encouraging, as they identify a model using imaging features and menopausal status to predict the luminal A, luminal B, HER2 and basal-like subtypes of breast cancer.

In the present study, imaging features were extracted from DCE-MRI scans of breast cancers. The experimental results demonstrated that the molecular subtype could be predicted based on imaging features from DCE-MRI. Future studies with larger sample sizes examining more radiomic features should be conducted to refine our findings.

Supporting information

[Figure omitted. See PDF.]

S1 Fig. ROC curve in the reproducibility dataset for discriminating the four molecular subtypes of breast cancer.

The classifiers based on dynamic features, morphologic features and first-orde r statistic features are shown. Features are combined to classify between (a) luminal A and non-luminal A tumors; (b) luminal B and non-luminal B tumors; (c) HER2-positive and non-HER2-positive tumors; and (d) basal-like and non-basal-like tumors.

https://doi.org/10.1371/journal.pone.0171683.s001

(TIFF)

S2 Fig. Bilateral asymmetric features between abnormal and normal breast differentiated tumors with four molecular subtypes.

https://doi.org/10.1371/journal.pone.0171683.s002

(TIFF)

S1 Table. Characteristics of the molecular subtypes of patients.

https://doi.org/10.1371/journal.pone.0171683.s003

(DOCX)

Acknowledgments

This work was supported in part by financial support from the National Natural Science Foundation of China (61271063 and 61401131), the National Basic Research Program of China (973 Program) (2013CB329502), and the Natural Science Foundation of Zhejiang Province of China (LZ15F010001 and LQ14F010011).

Author Contributions

1. Conceptualization: MF BZ LL.

2. Data curation: JZ.

3. Formal analysis: HL SW.

4. Funding acquisition: MF LL.

5. Investigation: LL.

6. Methodology: MF HL SW.

7. Project administration: LL.

8. Resources: JZ LL.

9. Writing – original draft: MF LL.

10. Writing – review & editing: MF LL.

Citation: Fan M, Li H, Wang S, Zheng B, Zhang J, Li L (2017) Radiomic analysis reveals DCE-MRI features for prediction of molecular subtypes of breast cancer. PLoS ONE 12(2): e0171683. https://doi.org/10.1371/journal.pone.0171683

References

1. Goldhirsch A, Wood WC, Coates AS, Gelber RD, Thurlimann B, Senn HJ. Strategies for subtypes—dealing with the diversity of breast cancer: highlights of the St. Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2011. Ann Oncol. 2011;22(8): 1736–47. pmid:21709140

2. Park S, Koo JS, Kim MS, Park HS, Lee JS, Kim SI, et al. Characteristics and outcomes according to molecular subtypes of breast cancer as classified by a panel of four biomarkers using immunohistochemistry. Breast. 2012;21(1): 50–7. pmid:21865043

3. Engstrom MJ, Opdahl S, Hagen AI, Romundstad PR, Akslen LA, Haugen OA, et al. Molecular subtypes, histopathological grade and survival in a historic cohort of breast cancer patients. Breast Cancer Res Treat. 2013;140(3): 463–73. pmid:23901018

4. Cheang MC, Chia SK, Voduc D, Gao D, Leung S, Snider J, et al. Ki67 index, HER2 status, and prognosis of patients with luminal B breast cancer. J Natl Cancer Inst. 2009;101(10): 736–50. pmid:19436038

5. Rouzier R, Perou CM, Symmans WF, Ibrahim N, Cristofanilli M, Anderson K, et al. Breast cancer molecular subtypes respond differently to preoperative chemotherapy. Clin Cancer Res. 2005;11(16): 5678–85. pmid:16115903

6. Carey LA, Dees EC, Sawyer L, Gatti L, Moore DT, Collichio F, et al. The triple negative paradox: primary tumor chemosensitivity of breast cancer subtypes. Clin Cancer Res. 2007;13(8): 2329–34. pmid:17438091

7. Yoo C, Ahn JH, Jung KH, Kim SB, Kim HH, Shin HJ, et al. Impact of immunohistochemistry-based molecular subtype on chemosensitivity and survival in patients with breast cancer following neoadjuvant chemotherapy. J Breast Cancer. 2012;15(2): 203–10. pmid:22807938

8. Sotiriou C, Pusztai L. Gene-expression signatures in breast cancer. N Engl J Med. 2009;360(8): 790–800. pmid:19228622

9. Dent R, Trudeau M, Pritchard KI, Hanna WM, Kahn HK, Sawka CA, et al. Triple-negative breast cancer: clinical features and patterns of recurrence. Clin Cancer Res. 2007;13(15 Pt 1): 4429–34. pmid:17671126

10. Crabb SJ, Cheang MC, Leung S, Immonen T, Nielsen TO, Huntsman DD, et al. Basal breast cancer molecular subtype predicts for lower incidence of axillary lymph node metastases in primary breast cancer. Clin Breast Cancer. 2008;8(3): 249–56. pmid:18650155

11. Hylton N. Dynamic contrast-enhanced magnetic resonance imaging as an imaging biomarker. J Clin Oncol. 2006;24(20): 3293–8. pmid:16829653

12. Zhao B, Schwartz LH, Larson SM. Imaging surrogates of tumor response to therapy: anatomic and functional biomarkers. J Nucl Med. 2009;50(2): 239–49. pmid:19164218

13. Parikh J, Selmi M, Charles-Edwards G, Glendenning J, Ganeshan B, Verma H, et al. Changes in primary breast cancer heterogeneity may augment midtreatment MR imaging assessment of response to neoadjuvant chemotherapy. Radiology. 2014;272(1): 100–12. pmid:24654970

14. Yamaguchi K, Abe H, Newstead GM, Egashira R, Nakazono T, Imaizumi T, et al. Intratumoral heterogeneity of the distribution of kinetic parameters in breast cancer: comparison based on the molecular subtypes of invasive breast cancer. Breast Cancer. 2015;22(5): 496–502. pmid:24402638

15. Gatenby RA, Grove O, Gillies RJ. Quantitative imaging in cancer evolution and ecology. Radiology. 2013;269(1): 8–15. pmid:24062559

16. Yang Q, Li L, Zhang J, Shao G, Zhang C, Zheng B. Computer-aided diagnosis of breast DCE-MRI images using bilateral asymmetry of contrast enhancement between two breasts. J Digit Imaging. 2014;27(1): 152–60. pmid:24043592

17. Yang Q, Li L, Zhang J, Shao G, Zheng B. A computerized global MR image feature analysis scheme to assist diagnosis of breast cancer: a preliminary assessment. Eur J Radiol. 2014;83(7): 1086–91. pmid:24743001

18. Yang Q, Li L, Zhang J, Shao G, Zheng B. A new quantitative image analysis method for improving breast cancer diagnosis using DCE-MRI examinations. Med Phys. 2015;42(1): 103–9. pmid:25563251

19. Yamamoto S, Maki DD, Korn RL, Kuo MD. Radiogenomic analysis of breast cancer using MRI: a preliminary study to define the landscape. AJR Am J Roentgenol. 2012;199(3): 654–63. pmid:22915408

20. Sutton EJ, Oh JH, Dashevsky BZ, Veeraraghavan H, Apte AP, Thakur SB, et al. Breast cancer subtype intertumor heterogeneity: MRI-based features predict results of a genomic assay. J Magn Reson Imaging. 2015;42(5): 1398–406. pmid:25850931

21. Kumar V, Gu Y, Basu S, Berglund A, Eschrich SA, Schabath MB, et al. Radiomics: the process and the challenges. Magn Reson Imaging. 2012;30(9): 1234–48. pmid:22898692

22. Mazurowski MA. Radiogenomics: what it is and why it is important. J Am Coll Radiol. 2015;12(8): 862–6. pmid:26250979

23. Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures, they are data. Radiology. 2016;278(2): 563–77. pmid:26579733

24. Grimm LJ. Breast MRI radiogenomics: Current status and research implications. J Magn Reson Imaging. 2016;43(6): 1269–78. pmid:26663695

25. Blaschke E, Abe H. MRI phenotype of breast cancer: kinetic assessment for molecular subtypes. J Magn Reson Imaging. 2015;42(4): 920–4. pmid:25758675

26. Bae MS, Seo M, Kim KG, Park IA, Moon WK. Quantitative MRI morphology of invasive breast cancer: correlation with immunohistochemical biomarkers and subtypes. Acta Radiol. 2015;56(3): 269–75. pmid:24558165

27. Costantini M, Belli P, Distefano D, Bufi E, Matteo MD, Rinaldi P, et al. Magnetic resonance imaging features in triple-negative breast cancer: comparison with luminal and HER2-overexpressing tumors. Clin Breast Cancer. 2012;12(5): 331–9. pmid:23040001

28. Sung JS, Jochelson MS, Brennan S, Joo S, Wen YH, Moskowitz C, et al. MR imaging features of triple-negative breast cancers. Breast J. 2013;19(6): 643–9. pmid:24015869

29. Uematsu T, Kasami M, Yuen S. Triple-negative breast cancer: correlation between MR imaging and pathologic findings. Radiology. 2009;250(3): 638–47. pmid:19244039

30. Koo HR, Cho N, Song IC, Kim H, Chang JM, Yi A, et al. Correlation of perfusion parameters on dynamic contrast-enhanced MRI with prognostic factors and subtypes of breast cancers. J Magn Reson Imaging. 2012;36(1): 145–51. pmid:22392859

31. Mazurowski MA, Zhang J, Grimm LJ, Yoon SC, Silber JI. Radiogenomic analysis of breast cancer: luminal B molecular subtype is associated with enhancement dynamics at MR imaging. Radiology. 2014;273(2): 365–72. pmid:25028781

32. Grimm LJ, Zhang J, Mazurowski MA. Computational approach to radiogenomics of breast cancer: luminal A and luminal B molecular subtypes are associated with imaging features on routine breast MRI extracted using computer vision algorithms. J Magn Reson Imaging. 2015;42(4): 902–7. pmid:25777181

33. Sutton EJ, Dashevsky BZ, Oh JH, Veeraraghavan H, Apte AP, Thakur SB, et al. Breast cancer molecular subtype classifier that incorporates MRI features. J Magn Reson Imaging. 2016.

34. Hammond ME, Hayes DF, Wolff AC, Mangu PB, Temin S. American society of clinical oncology/college of american pathologists guideline recommendations for immunohistochemical testing of estrogen and progesterone receptors in breast cancer. J Oncol Pract. 2010;6(4): 195–7. pmid:21037871

35. Wolff AC, Hammond ME, Hicks DG, Dowsett M, McShane LM, Allison KH, et al. Recommendations for human epidermal growth factor receptor 2 testing in breast cancer: American Society of Clinical Oncology/College of American Pathologists clinical practice guideline update. J Clin Oncol. 2013;31(31): 3997–4013. pmid:24101045

36. Zhang C, Li L, editors. 3D segmentation of masses in DCE-MRI images using FCM and adaptive MRF. Medical Imaging 2014: Image Processing, February 16, 2014–February 18, 2014; 2014; San Diego, CA, United states: SPIE.

37. Chen W, Giger ML, Li H, Bick U, Newstead GM. Volumetric texture analysis of breast lesions on contrast-enhanced magnetic resonance images. Magn Reson Med. 2007;58(3): 562–71. pmid:17763361

38. Scutt D, Lancaster GA, Manning JT. Breast asymmetry and predisposition to breast cancer. Breast Cancer Res. 2006;8(2): R14. pmid:16563179

39. Witten I, Frank E, Hall M Data mining: practical machine learing tools and techniques, 3rd ed. Amsterdam: Elsevier, http://www.cs.waikao.ac.nz/ml/weka/.

40. Just N. Improving tumour heterogeneity MRI assessment with histograms. Br J Cancer. 2014;111(12): 2205–13. pmid:25268373

41. King AD, Chow KK, Yu KH, Mo FK, Yeung DK, Yuan J, et al. Head and neck squamous cell carcinoma: diagnostic performance of diffusion-weighted MR imaging for the prediction of treatment response. Radiology. 2013;266(2): 531–8. pmid:23151830

42. Foroutan P, Kreahling JM, Morse DL, Grove O, Lloyd MC, Reed D, et al. Diffusion MRI and novel texture analysis in osteosarcoma xenotransplants predicts response to anti-checkpoint therapy. PLoS One. 2013;8(12): e82875. pmid:24358232

43. Voduc KD, Cheang MC, Tyldesley S, Gelmon K, Nielsen TO, Kennecke H. Breast cancer subtypes and the risk of local and regional relapse. J Clin Oncol. 2010;28(10): 1684–91. pmid:20194857

44. Arvold ND, Taghian AG, Niemierko A, Abi Raad RF, Sreedhara M, Nguyen PL, et al. Age, breast cancer subtype approximation, and local recurrence after breast-conserving therapy. J Clin Oncol. 2011;29(29): 3885–91. pmid:21900114

45. Metzger-Filho O, Sun Z, Viale G, Price KN, Crivellari D, Snyder RD, et al. Patterns of recurrence and outcome according to breast cancer subtypes in lymph node-negative disease: results from international Breast Cancer Study Group Trials VIII and IX. J Clin Oncol. 2013;31(25): 3083–90. pmid:23897954

46. Parmar C, Leijenaar RT, Grossmann P, Rios Velazquez E, Bussink J, Rietveld D, et al. Radiomic feature clusters and prognostic signatures specific for Lung and Head & Neck cancer. Sci Rep. 2015;5(11044. pmid:26251068

47. Fusco R, Sansone M, Filice S, Granata V, Catalano O, Amato DM, et al. Integration of DCE-MRI and DW-MRI Quantitative Parameters for Breast Lesion Classification. Biomed Res Int. 2015;2015(237863. pmid:26339597

48. Fusco R, Sansone M, Filice S, Carone G, Amato DM, Sansone C, et al. Pattern Recognition Approaches for Breast Cancer DCE-MRI Classification: A Systematic Review. J Med Biol Eng. 2016;36(4): 449–59. pmid:27656117