Gene Signatures and Cancer-Immune Phenotypes Based on m(6)A Regulators in Breast Cancer

The N(6)-methyladenosine (m(6)A) has been considered as a new layer of epitranscriptomic regulation on mRNA processing, stability, and translation. However, potential roles of m(6)A RNA methylation modification in tumor immune microenvironment (TIME) of breast cancer are yet fully understood. In this study, we comprehensively evaluated the genetic variations and transcript expressions of 15 m(6)A regulators in 1,079 breast cancer samples from the Cancer Genome Atlas (TCGA) database. We validated major regulators had significantly differential mRNA and protein expression in tumor tissue compared to normal tissues from 39 pairs of clinical breast cancer samples with different molecular subtypes, and especially high expression of m(6)A readers YTHDF1 and YTHDF3 predicted poor survival. Two clusters of breast cancer patients identified by the 15 m(6)A regulators’ pattern showed distinct overall survival, immune activation status, and immune cell infiltration, and clinical samples confirmed the diversity of lymphocytic infiltration. The profiles of these two clusters accorded with that of two classical cancer-immune phenotypes, immune-excluded and immune-inflamed phenotypes, it suggested that m(6)A regulators-based patterns might serve as crucial mediators of TIME in breast cancer. Moreover, the m(6)A phenotype-related gene signatures could also be survival predictor in breast cancer. Therefore, comprehensive evaluation of tumor m(6)A modification pattern will contribute to enhance our understanding of the characterization of immune cell infiltration in the tumor microenvironment and promote the responsiveness of breast cancer to immunotherapy.