Fig. 5

Bioinformatics prediction of circNR3C2 and hsa-miR-513a-3p. a Flowchart showing the prediction procedures of circular RNAs and miRNAs possibly regulating HRD1 expression in breast cancer. b Volcano plot showing an overview of the differential expression of circular RNAs identified in GSE113230, in 3 TNBC tissues compared to adjacent normal tissues. The upregulated and downregulated circRNAs in red and blue respectively were characterized as fold change (FC) ≥ 2 or ≤ 0.5 with adjusted P value < 0.05. c Venn diagram showing the overlap between the predicted mRNA-binding miRNAs and circRNA-sponged miRNAs, namely the miRNA candidates predicted to bind to HRD1 mRNA and sponged by the top 10 downregulated circRNAs in TNBC. d Interaction network of the circRNA and miRNA candidates. Size of the representative shapes (circle and triangle) depended on the node degree analyzed by Cytoscape software. e Gene Ontology (GO) enrichment analysis showing the biological processes and molecular functions (especially GO:0061630 and GO:0016567) in which genes targeted by the miRNA candidates were significantly enriched. f Heatmap and hierarchical clustering analysis showing the enrichment of miRNA candidates in the indicated biological processes with corresponding statistical significance. GO:0043687, post-translational protein modification. GO:0006464, cellular protein modification process. GO:0007173, epidermal growth factor signaling pathway. GO:0044267, cellular protein metabolic process. The color key reflects the P value. g KEGG pathway enrichment analysis showing the biological pathways in which genes targeted by the miRNA candidates were significantly enriched. h Heatmap and hierarchical clustering analysis showing the enrichment of miRNA candidates in the indicated biological pathway with corresponding statistical significance. Hsa04350, TGF-beta signaling pathway. The color key reflects the P value. i Correlation between the expression of circNR3C2 and HRD1 in 30 TNBC tissues