Identification of circRNA-miRNA-mRNA network in luminal breast cancers by integrated analysis of microarray datasets

Introduction: Circular RNAs (circRNAs) regulatory network is important in human cancer. We, therefore, mapped the regulatory networks driven by circRNA in luminal-subtype breast cancer. Methods: Breast cancer-related microarray datasets from GEO database were analyzed for the differentially expressed circRNAs, miRNAs, and mRNAs. The potential downstream RNAs were collected using Circular RNA Interactome or Targetscan database. Protein-protein interaction (PPI) analysis was performed for the filtered genes to identify hub genes. The functions were annotated by the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. CircRNA-miRNA-mRNA networks were mapped using Cytoscape software. Hsa_circ_0086735-miR-1296-5p-STAT1 axis was used for verification. The expression levels of hsa_circ_0086735, miR-1296-5p, and STAT1 mRNA were confirmed by qRT-PCR in luminal-subtype tissues and cell lines. The interactions among them were verified by Luciferase reporter assay and RNA pull-down assay. Cell proliferation and apoptosis were assayed. Overall and distant metastasis-free survival was analyzed. Results: A total of 70 genes were finally targeted and enriched in multi-process and multi-pathway. Networks containing 96 circRNA-miRNA-mRNA axes were constructed. Hsa_circ_0086735 and STAT1 mRNA was upregulated in luminal breast cancer, while miR-1296-5p was downregulated. Hsa_circ_0086735-miR-1296-5p-STAT1 axis promotes breast cancer progression and contributes to tamoxifen resistance. High hsa_circ_0086735 was associated with poor overall and distant metastasis-free survival. Discussion: This study identified the hsa_circ_0086735-miR-1296-5p-STAT1 as an important regulatory axis in luminal-subtype breast cancer, aiding to determine potential therapeutic targets.