Identifying subtypes and developing prognostic models based on N6-methyladenosine and immune microenvironment related genes in breast cancer

Abstract

Breast cancer (BC) is the most prevalent cancer in women globally. The tumor microenvironment (TME), comprising epithelial tumor cells and stromal elements, is vital for breast tumor development. N6-methyladenosine (m6A) modification plays a key role in RNA metabolism, influencing its various aspects such as stability and translation. There is a notable link between m6A methylation and immune cells in the TME, although this relationship is complex and not fully deciphered. In this research, BC expression and clinicopathological data from TCGA were scrutinized to assess expression profiles, mutations, and CNVs of 31 m6A genes and immune microenvironment-related genes, examining their correlations, functions, and prognostic impacts. Lasso and Cox regression identified prognostic genes for constructing a nomogram. Single-cell analyses mapped the distribution and patterns of these genes in BC cell development. We investigated associations between gene-derived risk scores and factors like immune infiltration, TME, checkpoints, TMB, CSC indices, and drug response. As a complement to computational analyses, in vitro experiments were conducted to confirm these expression patterns. We included 31 m6A regulatory genes and discovered a correlation between these genes and the extent of immune cell infiltration. Subsequently, a 7-gene risk score was generated, encompassing HSPA2, TAP1, ULBP2, CXCL1, RBP1, STC2, and FLT3. It was observed that the low-risk group exhibited better overall survival (OS) in BC, with higher immune scores but lower tumor mutational burden (TMB) and cancer stem cell (CSC) indices, as well as lower IC50 values for commonly used drugs. To enhance clinical applicability, age and stage were incorporated into the risk score, and a more comprehensive nomogram was constructed to predict OS. This nomogram was validated and demonstrated good predictive performance, with area under the curve (AUC) values for 1-year, 3-year, and 5-year OS being 0.848, 0.807, and 0.759, respectively. Our findings highlight the profound impact of prognostic-related genes on BC immune response and prognostic outcomes, suggesting that modulation of the m6A-immune pathway could offer new avenues for personalized BC treatment and potentially improve clinical outcomes.