Integrative genomic analysis of N6-methyladenosine-single nucleotide polymorphisms (m⁶A-SNPs) associated with breast...

Due to the important role of N6-methyladenosine (m⁶A) in breast cancer, single nucleotide polymorphisms (SNPs) in genes with m⁶A modification may also be involved in breast cancer pathogenesis. In this study, we used a public genome-wide association study dataset to identify m⁶A-SNPs associated with breast cancer and to further explore their potential functions. We found 113 m⁶A-SNPs associated with breast cancer that reached the genome-wide suggestive threshold (5.0E-05), and 86 m⁶A-SNPs had eQTL signals. Only six genes were differentially expressed between controls and breast cancer cases in GEO datasets (GSE15852, GSE115144, and GSE109169), and the SNPs rs4829 and rs9610915 were located next to the m⁶A modification sites in the 3'UTRs of TOM1L1 and MAFF, respectively. In addition, we found that polyadenylate-binding protein cytoplasmic 1 might have a potential interaction with rs4829 (TOM1L1) and rs9610915 (MAFF). In summary, these findings indicated that the SNPs rs4829 and rs9610915 are potentially associated with breast cancer because they had eQTL signals, altered gene expression, and were located next to the m⁶A modification sites in the 3'UTRs of their coding genes. However, further studies are still needed to clarify how genetic variation affects the epigenetic modification, m⁶A, and its subsequent functions in the pathogenesis of breast cancer.