The MRI-based radiomic approach developed in this work, built on a comprehensive BC dataset including MRI sequences and molecular outcomes, can efficiently characterize BC heterogeneity according to the most examined biomarkers (ER, PR, HER2, and Ki67). This evidence concerns the gene ESR1 and breast cancer.