Contrast-Enhanced Spectral Mammography (CESM) is a recently introduced mammographic method with characteristics particularly suitable for breast cancer radiomic analysis. This work aims to evaluate radiomic features for predicting histological outcome and two cancer molecular subtypes, namely Human Epidermal growth factor Receptor 2 (HER2)-positive and triple-negative. From 52 patients, 68 lesions were identified and confirmed on histological examination. Radiomic analysis was performed on regions of interest (ROIs) selected from both low-energy (LE) and ReCombined (RC) CESM images. Fourteen statistical features were extracted from each ROI. Expression of estrogen receptor (ER) was significantly correlated with variation coefficient and variation range calculated on both LE and RC images; progesterone receptor (PR) with skewness index calculated on LE images; and Ki67 with variation coefficient, variation range, entropy and relative smoothness indices calculated on RC images. HER2 was significantly associated with relative smoothness calculated on LE images, and grading tumor with variation coefficient, entropy and relative smoothness calculated on RC images. Encouraging results for differentiation between ER+/ER−, PR+/PR−, HER2+/HER2−, Ki67+/Ki67−, High-Grade/Low-Grade and TN/NTN were obtained. Specifically, the highest performances were obtained for discriminating HER2+/HER2− (90.87%), ER+/ER− (83.79%) and Ki67+/Ki67− (84.80%). Our results suggest an interesting role for radiomics in CESM to predict histological outcomes and particular tumors’ molecular subtype.

In recent years personalized medicine reached an increasing importance, especially in the design of oncological therapies. In particular, the development of patients’ profiling strategies suggests the possibility of promising rewards. In this work, we present an explainable artificial intelligence (XAI) framework based on an adaptive dimensional reduction which (i) outlines the most important clinical features for oncological patients’ profiling and (ii), based on these features, determines the profile, i.e., the cluster a patient belongs to. For these purposes, we collected a cohort of 267 breast cancer patients. The adopted dimensional reduction method determines the relevant subspace where distances among patients are used by a hierarchical clustering procedure to identify the corresponding optimal categories. Our results demonstrate how the molecular subtype is the most important feature for clustering. Then, we assessed the robustness of current therapies and guidelines; our findings show a striking correspondence between available patients’ profiles determined in an unsupervised way and either molecular subtypes or therapies chosen according to guidelines, which guarantees the interpretability characterizing explainable approaches to machine learning techniques. Accordingly, our work suggests the possibility to design data-driven therapies to emphasize the differences observed among the patients.

Cancer treatment planning benefits from an accurate early prediction of the treatment efficacy. The goal of this study is to give an early prediction of three-year Breast Cancer Recurrence (BCR) for patients who underwent neoadjuvant chemotherapy. We addressed the task from a new perspective based on transfer learning applied to pre-treatment and early-treatment DCE-MRI scans. Firstly, low-level features were automatically extracted from MR images using a pre-trained Convolutional Neural Network (CNN) architecture without human intervention. Subsequently, the prediction model was built with an optimal subset of CNN features and evaluated on two sets of patients from I-SPY1 TRIAL and BREAST-MRI-NACT-Pilot public databases: a fine-tuning dataset (70 not recurrent and 26 recurrent cases), which was primarily used to find the optimal subset of CNN features, and an independent test (45 not recurrent and 17 recurrent cases), whose patients had not been involved in the feature selection process. The best results were achieved when the optimal CNN features were augmented by four clinical variables (age, ER, PgR, HER2+), reaching an accuracy of 91.7% and 85.2%, a sensitivity of 80.8% and 84.6%, a specificity of 95.7% and 85.4%, and an AUC value of 0.93 and 0.83 on the fine-tuning dataset and the independent test, respectively. Finally, the CNN features extracted from pre-treatment and early-treatment exams were revealed to be strong predictors of BCR.

To date, some artificial intelligence (AI) methods have exploited Dynamic Contrast-Enhanced Magnetic Resonance Imaging (DCE-MRI) to identify finer tumor properties as potential earlier indicators of pathological Complete Response (pCR) in breast cancer patients undergoing neoadjuvant chemotherapy (NAC). However, they work either for sagittal or axial MRI protocols. More flexible AI tools, to be used easily in clinical practice across various institutions in accordance with its own imaging acquisition protocol, are required. Here, we addressed this topic by developing an AI method based on deep learning in giving an early prediction of pCR at various DCE-MRI protocols (axial and sagittal). Sagittal DCE-MRIs refer to 151 patients (42 pCR; 109 non-pCR) from the public I-SPY1 TRIAL database (DB); axial DCE-MRIs are related to 74 patients (22 pCR; 52 non-pCR) from a private DB provided by Istituto Tumori “Giovanni Paolo II” in Bari (Italy). By merging the features extracted from baseline MRIs with some pre-treatment clinical variables, accuracies of 84.4% and 77.3% and AUC values of 80.3% and 78.0% were achieved on the independent tests related to the public DB and the private DB, respectively. Overall, the presented method has shown to be robust regardless of the specific MRI protocol.

In the absence of lymph node abnormalities detectable on clinical examination or imaging, the guidelines provide for the dissection of the first axillary draining lymph nodes during surgery. It is not always possible to arrive at surgery without diagnostic doubts, and machine learning algorithms can support clinical decisions. The web calculator CancerMath (CM) allows you to estimate the probability of having positive lymph nodes valued on the basis of tumor size, age, histologic type, grading, expression of estrogen receptor, and progesterone receptor. We collected 993 patients referred to our institute with clinically negative results characterized by sentinel lymph node status, prognostic factors defined by CM, and also human epidermal growth factor receptor 2 (HER2) and Ki-67. Area Under the Curve (AUC) values obtained by the online CM application were comparable with those obtained after training its algorithm on our database. Nevertheless, by training the CM model on our dataset and using the same feature, we reached a sensitivity median value of 72%, whereas the online one was equal to 46%, despite a specificity reduction. We found that the addition of the prognostic factors Her2 and Ki67 could help improve performances on the classification of particular types of patients with the aim of reducing as much as possible the false positives that lead to axillary dissection. As showed by our experimental results, it is not particularly suitable for use as a support instrument for the prediction of metastatic lymph nodes on clinically negative patients.

Nuclei identification is a fundamental task in many areas of biomedical image analysis related to computational pathology applications. Nowadays, deep learning is the primary approach by which to segment the nuclei, but accuracy is closely linked to the amount of histological ground truth data for training. In addition, it is known that most of the hematoxylin and eosin (H&E)-stained microscopy nuclei images contain complex and irregular visual characteristics. Moreover, conventional semantic segmentation architectures grounded on convolutional neural networks (CNNs) are unable to recognize distinct overlapping and clustered nuclei. To overcome these problems, we present an innovative method based on gradient-weighted class activation mapping (Grad-CAM) saliency maps for image segmentation. The proposed solution is comprised of two steps. The first is the semantic segmentation obtained by the use of a CNN; then, the detection step is based on the calculation of local maxima of the Grad-CAM analysis evaluated on the nucleus class, allowing us to determine the positions of the nuclei centroids. This approach, which we denote as NDG-CAM, has performance in line with state-of-the-art methods, especially in isolating the different nuclei instances, and can be generalized for different organs and tissues. Experimental results demonstrated a precision of 0.833, recall of 0.815 and a Dice coefficient of 0.824 on the publicly available validation set. When used in combined mode with instance segmentation architectures such as Mask R-CNN, the method manages to surpass state-of-the-art approaches, with precision of 0.838, recall of 0.934 and a Dice coefficient of 0.884. Furthermore, performance on the external, locally collected validation set, with a Dice coefficient of 0.914 for the combined model, shows the generalization capability of the implemented pipeline, which has the ability to detect nuclei not only related to tumor or normal epithelium but also to other cytotypes.