BackgroundBreast cancer is a malignant tumour that seriously threatens women’s life and health and exhibits high inter-individual heterogeneity, emphasising the need for more in-depth research on its pathogenesis. While internal 7-methylguanosine (m7G) modifications affect RNA processing and function and are believed to be involved in human diseases, little is currently known about the role of m7G modification in breast cancer.Methods and ResultsWe elucidated the expression, copy number variation incidence and prognostic value of 24 m7G-related genes (m7GRGs) in breast cancer. Subsequently, based on the expression of these 24 m7GRGs, consensus clustering was used to divide tumour samples from the TCGA-BRCA dataset into four subtypes based on significant differences in their immune cell infiltration and stromal scores. Differentially expressed genes between subtypes were mainly enriched in immune-related pathways such as ‘Ribosome’, ‘TNF signalling pathway’ and ‘Salmonellainfection’. Support vector machines and multivariate Cox regression analysis were applied based on these 24 m7GRGs, and four m7GRGs—AGO2, EIF4E3, DPCS and EIF4E—were identified for constructing the prediction model. An ROC curve indicated that a nomogram model based on the risk model and clinical factors had strong ability to predict the prognosis of breast cancer. The prognoses of patients in the high- and low-TMB groups were significantly different (p= 0.03). Moreover, the four-gene signature was able to predict the response to chemotherapy.ConclusionsIn conclusion, we identified four different subtypes of breast cancer with significant differences in the immune microenvironment and pathways. We elucidated prognostic biomarkers associated with breast cancer and constructed a prognostic model involving four m7GRGs. In addition, we predicted the candidate drugs related to breast cancer based on the prognosis model.Supplementary InformationThe online version contains supplementary material available at 10.1186/s12885-023-11012-z.