Abstract
In this study, we identified a key enhancer RNA (eRNA) region in breast cancer (BRCA) by applying an integrated analysis method. Reported eRNA region and genes affected by them were selected as presumed target pairs. Kaplan-Meier (KM) survival and correlation analyses were performed to screen valuable eRNA region. Based on the KM value and its correlation with the paired target genes, we carefully selected ELOVL2-AS1 as a potential key eRNA region in BRCA. Subsequently, we analyzed the expression of ELOVL2-AS1 and ELOVL2 in four BRCA subtypes and in different BRCA cell lines. The expression of ELOVL2-AS1 and ELOVL2 in triple negative breast cancer (TNBC) was significantly lower than those in Luminal A. After that, we analyzed the function of genes that are positively correlated with ELOVL2-AS1. We found that the co-expression gene mainly related to cilia and cilia characteristics of TNBC is significantly weaker than that of Luminal A. Considering the stronger invasion and metastasis of TNBC (compared with Luminal A) and the close relationship between decreased cilia and metastasis, we overexpressed ELOVL2-AS1 in TNBC and observed its effect on cell migration. The results show that it can inhibit the migration of TNBC. Finally, we analyzed the assay for transposase-accessible chromatin sequencing data, chromatin interaction analysis with paired-end tag sequencing data, and chromatin immunoprecipitation sequencing data and identified the chromatin interaction between ELOVL2-AS1 and ELOVL2, suggesting a direct regulatory interaction.