BackgroundPredicting breast cancer molecular subtypes can help guide individualized clinical treatment of patients who need the rational preoperative treatment. This study aimed to investigate the efficacy of preoperative prediction of breast cancer molecular subtypes by contrast-enhanced mammography (CEM) radiomic features. MethodsThis retrospective two-centre study included women with breast cancer who underwent CEM preoperatively between August 2016 and May 2022. We included 356 patients with 386 lesions, which were grouped into training (n=162), internal test (n=160) and external test sets (n=64). Radiomics features were extracted from low-energy (LE) images and recombined (RC) images and selected. Three dichotomous tasks were established according to postoperative immunohistochemical results: Luminal vs. non-Luminal, human epidermal growth factor receptor (HER2)-enriched vs. non-HER2-enriched, and triple-negative breast cancer (TNBC) vs. non-TNBC. For each dichotomous task, the LE, RC and LE+RC radiomics models were built by the support vector machine (SVM) classifier. The prediction performance of the models was assessed by the area under the receiver operating characteristic curve (AUC). Then, the accuracy, sensitivity, specificity, positive predictive value, and negative predictive value were calculated for the models. DeLong’s test was utilized to compare the AUCs. ResultsRadiomics models based on CEM are valuable for predicting breast cancer molecular subtypes. The LE+RC model achieved the best performance in the test set. The LE+RC model predicted Luminal, HER2-enriched, and TNBC subtypes with AUCs of 0.93, 0.89, and 0.87 in the internal test set and 0.82, 0.83 and 0.69 in the external test set, respectively. In addition, the LE model performed more satisfactorily than the RC model. ConclusionCEM radiomics features can effectively predict breast cancer molecular subtypes preoperatively, and the LE+RC model has the best predictive performance.