Abstract
BACKGROUND: Recent studies have identified a complex relationship between methylation patterns and the development of various cancers. Breast cancer (BC) is the second leading cause of cancer mortality among women. Approximately 5-20% of BC patients are at risk of BC brain metastases (BCBM). Although 5-methylcytosine (m5C) has been identified as an important regulatory modifier, its distribution in BCBM is not well understood. This study aimed to investigate the distribution of m5C in BCBM. MATERIALS AND METHODS: Samples from BCBM (231-BR cells) and BC (MDA-MB-231 cells) groups were subjected to a comprehensive analysis of the m5C methylation in long non-coding RNA (lncRNA) using methylated RNA immunoprecipitation next-generation sequencing (MeRIP-seq). The expression levels of methylated genes in BC and adjacent tissues were verified through quantitative real-time polymerase chain reaction (RT-qPCR). Enrichment pathway analyses were through Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) to predict the potential functions of m5C in BCBM. RESULTS: The MeRIP-seq analysis identified 23,934 m5C peaks in BCBM and 21,236 m5C in BC. A total of 9,480 annotated genes (BCBM) and 8,481 annotated genes (BC) were mapped. Notably, 1,819 methylation sites in lncRNA were upregulated in BCBM, whereas 2,415 methylation sites were upregulated in BC. Significant m5C hypermethylated lncRNAs included ENST00000477316, ENST00000478098 and uc002gtt.1, whereas hypomethylated lncRNAs included ENST00000600912, ENST00000493668, ENST00000544651 and ENST00000464989. These results were verified by qPCR and MeRIP-qPCR in BC and BCBM. Considering the strong association between m5C RNA methylation regulators and lncRNA, we examined the expression levels of 13 m5C RNA methylation regulators and observed significant differences between BC tissues and adjacent normal tissues. In addition, the interaction between regulators of altered expression and the differentially expressed genes in vitro was analyzed. The GO and KEGG pathways analyses revealed that genes significantly associated with m5C sites in lncRNA were linked to the BCBM signaling pathways. CONCLUSION: This uncovered significant variations in the levels and distribution of m5C in BCBM compared to BC. The findings provide a new theoretical understanding of the mechanisms of BCBM.