To develop a simple and convenient method based on multiparametric magnetic resonance imaging (MRI) and clinical features to non-invasively predict tumour-infiltrating lymphocytes (TILs) in breast cancer (BC) and to explore the relationship between TIL levels and disease-free survival (DFS).A total of 172 BC patients were enrolled between November 2017 and June 2021 in this retrospective study. The patients were divided into high (≥10%) and low (<10%) TIL groups. Clinicopathological data were collected. MRI features were reviewed by two radiologists. Predictors associated with TILs were determined by using multivariable logistic regression analyses. Kaplan-Meier survival curves based on TIL levels were used to estimate DFS.A total of 102 patients with low TILs and 70 patients with high TILs were included in the study. Tumour size (odds ratio [OR], 1.040; 95% confidence interval [CI]: 1.006, 1.075; p=0.020), apparent diffusion coefficient (ADC; OR, 1.003; 95% CI: 1.001, 1.005; p=0.015), clinical axillary lymph node status (CALNS; OR, 3.222; 95% CI: 1.372,7.568; p=0.007), and enhancement pattern (OR, 0.284; 95% CI: 0.143, 0.563; p<0.001) were independently associated with TIL levels. These features were used in the ALSE model (where A is ADC, L is CALNS, S is size, and E is enhancement pattern). High TILs were associated with better DFS (p=0.016).The ALSE model derived from multiparametric MRI and clinical features could non-invasively predict TIL levels in BC, and high TILs were associated with longer DFS, especially in human epidermal growth factor receptor 2 (HER2)-positive BC and triple-negative BC (TNBC).