Abstract
Brain metastasis in basal-like breast cancer poses a significant challenge in cancer management due to its aggressive nature and limited treatment options. This study conducted a comprehensive analysis to explore the potential role of circular RNAs (circRNAs) as members of endogenous networks in developing breast cancer brain metastasis. Here, we utilized RNA sequencing data from primary breast cancer and brain metastasis tissue with basal-like subtype (n = 11). After quality controlling and preprocessing of fastq files, gene expression of mRNA and circRNAs were extracted from matched samples and normalized. Then, we employed the weighted gene co-expression network analysis approach to identify brain metastasis-associated circRNA modules ( Spearman Correlation > 0.5 , P - value < 0.05 ). Moreover, we found five protein-coding genes of PHLDA1, SLC12A2, MMP2, RGP1, and MAP2K6, significantly upregulated in brain metastatic tissues compared to primary breast cancer ( FDR < 0.05 ). These genes were enriched in the "GnRH signaling pathway" and "Fluid shear stress and atherosclerosis" pathways ( FDR < 0.05 ). Next, to explore the potential interactions between circRNAs and protein-coding genes, we reconstructed a competing endogenous RNA (ceRNA) network using mutual miRNAs between the circRNA module and upregulated mRNAs. Notably, we could detect two axes of circ_0087558/miR-604/MMP2 and MMP2/miR-1248/Circ_0087558/miR-643/MAP2K6 in ceRNA network. In conclusion, the identified circRNA-miRNA-mRNA axes might be therapeutic targets or diagnostic biomarkers for this challenging subtype of breast cancer. However, due to the small number of samples, further experimental validations are essential.