Accurate subtyping of breast cancer is crucial for its diagnosis, management, and prognostication. This study aimed to determine the association of magnetic resonance imaging (MRI) breast features with the molecular subtype and aggressiveness of breast cancer in a multi-ethnic population. Treatment-naive patients with invasive breast carcinoma were included in this retrospective study. Breast MRI features were recorded based on the American College of Radiology-Breast Imaging Reporting and Data System (ACR-BIRADS) criteria, with tumour size, and apparent diffusion coefficient value (ADC). The statistical association was tested with Pearson Chi-Square Test of Independence for categorical data or the Kruskal-Wallis/ Mann Whitney U test for numerical data between the MRI features and molecular subtype, receptor status, tumour grade, lymphovascular infiltration (LVI) and axillary lymph node (ALN). Multinomial logistic regression was used to test the predictive likelihood of the significant features. The breast cancer subtypes were determined via immunohistochemistry (IHC) and dual-color dual-hapten in-situ hybridization (D-DISH). The expression statuses of ER, PR, and HER-2, LVI, and ALN were obtained from the histopathology report. The ER / PR / HER-2 was evaluated according to the American Society of Clinical Oncology / College of American Pathologists. The study included 194 patients; 41.8% (n = 81) Chinese, 40.7% (n = 79) Malay, and 17.5% (n = 34) Indian, involving 71.6%(n = 139) luminal-like, 12.9%(n = 25) HER-2 enriched, and 15.5%(n = 30) Triple-negative breast cancer (TNBC). TNBC was associated with rim enhancement (p = 0.002) and peritumoral oedema (p = 0.004). HER-2 enriched tumour was associated with larger tumour size (p = 0.041). Luminal-like cancer was associated with irregular shape (p = 0.005) with circumscribed margin (p = 0.003). Other associations were ER-negative tumour with circumscribed margin (p = 0.002) and PR-negative with round shape (p = 0.001). Tumour sizes were larger in ER-negative (p = 0.044) and PR-negative (p = 0.022). Rim enhancement was significantly associated with higher grade (p = 0.001), and moderate peritumoral oedema with positive axillary lymph node (p = 0.002). Certain MRI features can be applied to differentiate breast cancer molecular subtypes, receptor status and aggressiveness, even in a multi-ethnic population.