Cerna Solutions

Nov 27, 2023

Abstract Breast cancer is a major global health concern, and recent researches have highlighted the critical roles of non-coding RNAs in both cancer and the immune system. The competing endogenous RNA hypothesis suggests that various types of RNA, including coding and non-coding RNAs, compete for microRNA targets, acting as molecular sponges. This study introduces the Pre_CLM_BCS pipeline to investigate the potential of long non-coding RNAs and circular RNAs as biomarkers in breast cancer subtypes. The pipeline identifies specific modules within each subtype that contain at least one long non-coding RNA or circular RNA exhibiting significantly distinct expression patterns when compared to other subtypes. The results reveal potential biomarker genes for each subtype, such as circ_001845, circ_001124, circ_003925, circ_000736, and circ_003996 for the basal-like subtype, circ_00306 and circ_00128 for the luminal B subtype, circ_000709 and NPHS1 for the normal-like subtype, CAMKV and circ_001855 for the luminal A subtype, and circ_00128 and circ_00173 for the HER2+ subtype. Additionally, certain long non-coding RNAs and circular RNAs, including RGS5-AS1, C6orf223, HHLA3-AS1, circ_000349, circ_003996, circ_003925, circ_002665, circ_001855, and DLEU1, are identified as potential regulators of T cell mechanisms, underscoring their importance in understanding breast cancer progression in various subtypes. This pipeline provides valuable insights into cancer and immune-related processes in breast cancer subtypes. Introduction Breast cancer continues to hold its position as the prevailing type of cancer, constituting a significant portion of diagnoses among women globally, with approximately one in three cases being attributed to this malignancy. Statistical data from 2021 reveal that an estimated 281,550 new cases were diagnosed; accounting for 14.8% of all cancer diagnoses, and tragically, breast cancer was responsible for 43,600 deaths, representing 7.2% of the total cancer-related fatalities 1 . Characterized by its immunohistochemical properties, hormone receptors, and specific protein involvement, breast cancer is classified into five primary molecular subtypes: basal-like, luminal B, normal-like, luminal A, and HER2+ 2 . These subtypes exhibit significant differences in characteristics 3 . Although substantial advancements have been made in clinical treatment, their efficacy could be further improved with a more precise understanding of the distinguishing features and molecular mechanisms associated with each breast cancer subtype 4 , 5 . This knowledge could significantly influence treatment outcomes and patient survival rates 6 . Non-coding RNAs (ncRNAs) play pivotal roles in the regulating of gene expression 7 . The ncRNAs include long ncRNAs (lncRNAs), microRNAs (miRNAs), and circular RNAs (circRNAs). MiRNAs are short, single-stranded ncRNA molecules. They primarily modulate gene expression by binding to target mRNAs leading to their degradation or repression of translation 8 . LncRNAs are a broad class of ncRNAs exceeding 200 nucleotides in length. They regulate mRNA production by influencing transcription of protein-coding genes 9 . Moreover, lncRNAs may affect tumorigenesis by controlling key cancer-related genes 10 . CircRNAs are a relatively new class of endogenous small ncRNAs. Unlike most RNAs, they lack both a 5′ cap and 3′ poly-A tail, typically resulting from splicing errors 11 . Advances in RNA sequencing technologies and bioinformatics have accelerated research on circRNAs, highlighting their significant role in gene expression regulation, primarily by acting as miRNA sponges 12 . Similar to miRNAs and lncRNAs, circRNAs have been linked to various complex human diseases 13 , 14 , including several types of cancers 15 . Their dysregulation and potential functional roles make circRNAs promising candidates for further exploration and potential therapeutic interventions in cancer research. While the roles of miRNAs, lncRNAs, and circRNAs in cancer pathology and physiology are increasingly being elucidated, their interaction within competing endogenous RNAs (ceRNA) networks promises new insights into their regulatory mechanisms 16 . Some studies suggest that ceRNAs are pivotal to many vital biological processes and could serve as potential diagnostic markers, prognostic indicators, or therapeutic targets 17 , 18 . The regulatory influence of miRNAs goes beyond their capacity to recognize target sites on various RNA molecules 19 . It has been suggested that miRNAs can facilitate regulatory crosstalk among different components of the transcriptome 20 . This intricate miRNA-mediated regulation is modulated by other RNA molecules, including mRNAs and ncRNAs, which have been found to interact with miRNAs and exhibit significant expression across diverse biological conditions. This mechanism introduces an additional layer of post-transcriptional gene regulation, providing a complementary perspective on the functional relevance of the vast number of transcribed, yet un-translated RNAs. The concept of ceRNAs further enhances this post-transcriptional regulatory landscape, where ncRNAs gain new significance. The interplay mediated by miRNAs among different types of RNA molecules has been observed in numerous contexts 21 , 22 . The ceRNAs hypothesis posits that a single miRNA can regulate multiple target RNAs, as long as they contain the specific miRNA response element. This hypothesis introduces the notion that RNAs can compete for a limited pool of miRNAs, serving as the foundation for constructing a ceRNA network 23 . Recently, the emergence of public databases has facilitated ceRNA network reconstruction, shedding light on the previously elusive mechanisms of tumorigenesis in cancers, such as breast cancer 24 . Although previous studies have separately analyzed lncRNA–miRNA–mRNA 25 and circRNA–miRNA–mRNA 26 ceRNA networks, no research to date has integrated both circRNAs and lncRNAs into the ceRNA network for the basal-like, luminal B, normal-like, luminal A, and HER2+ breast cancer subtypes. Examining ceRNA networks across these subtypes, and combining circRNAs and lncRNAs, can significantly enhance our understanding of immune-related tumorigenesis in various breast cancer subtypes. In this study, we present a novel pipeline, called Pre_CLM_BCS, designed to identify potential circRNAs, lncRNAs, and mRNAs as effective biomarkers for each breast cancer subtype. Our pipeline initiates by constructing a primary ceRNA network for each cancer subtype using gene expression data from patients. Subsequently, specific ceRNA networks are extracted for each breast cancer subtype from their respective primary ceRNA networks. Within these specific networks, modules containing at least one lncRNA or circRNA are identified and enriched with gene ontology terms and pathways. Finally, we employ survival analysis and differential gene expression to predict biomarkers from the enriched genes in the modules, which are further evaluated by support vector machine. The predicted biomarkers for each breast cancer subtype according to the pipelines are as follows: For the basal-like subtype: circ_001845, circ_001124, circ_003925, circ_000736, and circ_003996. For the luminal B subtype: circ_00306 and circ_00128. For the normal-like subtype: circ_000709 and nephrotic syndrome 1 (NPHS1 as an mRNA). For the luminal A subtype: cam kinase like vesicle associated (CAMKV as an mRNA) and circ_001855. For the HER2+ subtype: circ_00128 and circ_00173. Additionally, we find the following RNAs in modules that have regulatory effects on T cell mechanisms, further emphasizing the significance of these RNAs in understanding and addressing breast cancer progression across various subtypes: For the basal-like subtype: circ_003996, and circ_003925. For the luminal B subtype: chromosome 6 open reading frame 223 (C6orf223 as an lncRNA), HHLA3 Antisense RNA 1 (HHLA3-AS1 as an lncRNA), and circ_000349. For the normal-like subtype: deleted in lymphocytic leukemia 1 (DLEU1 as an lncRNA). For the luminal A subtype: circ_002665 and circ_001855. For the HER2+ subtype: regulator of G protein signaling 5 antisense RNA 1 (RGS5-AS1 as an lncRNA). Materials and methods In this section, our main goal is to introduce the proposed pipeline, Pre_CLM_BCS, for predicting suitable circRNAs, lncRNAs, and mRNAs as biomarkers for five breast cancer subtypes: basal-like, luminal B, normal-like, luminal A, and HER2+. Figure  1 visually outlines our pipeline, and the following subsections offer in-depth explanations for each step of the process. Figure 1 $$BCS = \left\{ {basal{-}like, luminal\;B, normal{-}like, luminal\;A, HER2 + } \right\}.$$ In this study, we extract miRNA, mRNA, and lncRNA gene expression data, along with clinical data for breast cancer samples from The Cancer Genome Atlas (TCGA) database (as of August 2015). We utilize the miRNA-seq dataset for miRNA gene expression and the RNA-SeqV2 level 3 dataset for mRNA and lncRNA expression profiles. To annotate the mRNA and lncRNA genes from the RNA-SeqV2 level 3 dataset, we use the BioMart tool available on the Ensembl website ( www.ensembl.org ) as a genome browser 27 . Concurrently, circRNA expression data for breast cancer subtypes is sourced based on the findings of Asha et al. 28 . Additionally, we collect clinical data, including immunohistochemical information, for breast cancer samples from TCGA database, in order to conduct survival analysis. Considering the conventional categorization of breast cancer into five subtypes, our study samples are categorized accordingly. This results in the creation of five distinct datasets for mRNA, lncRNA, circRNA, and miRNA. For each cancer subtype, we establish a dataset denoted as \(\Delta^{C}\), where \(C \in BCS\). Each dataset encapsulates the gene expression of mRNAs, lncRNAs, circRNAs, and miRNAs for samples within the respective breast cancer subtypes. The data preprocessing for each dataset \(\Delta^{C} \) comprises three main stages: 1. 2. 3. Detection of outliers. To identify outliers, our pipeline utilizes hierarchical clustering of the samples based on their expression profiles within the R environment (version 3.6) 29 . Primary ceRNA networks In the second step of the Pre_CLM_BCS pipeline, we construct a primary ceRNA network named \(P^{C} = \left( {N^{C} ,E^{C} } \right)\) for each breast cancer subtype \(C \in BCS\) using gene expression data from dataset \(\Delta^{C}\). Here, \(N^{C} \) represents the collection of mRNAs, lncRNAs, and circRNAs in the dataset as nodes. The set of edges (interactions) in this network is defined as \( E^{C}\), where each edge between two RNAs from set \(N^{C}\) indicates their interaction 30 . To identify these interactions, we utilize the Pearson correlation coefficient (PCC) for the gene expression dataset \(\Delta^{C}.\) Our approach for constructing a primary ceRNA network is rooted in the ceRNA regulatory model. According to this model, lncRNAs and circRNAs compete with miRNAs, resulting in opposing effects on mRNA expression 31 , 32 . Therefore, when the expression of an miRNA increases, there is a corresponding decrease in the levels of lncRNAs, mRNAs, and circRNAs. This competition weakens the positive correlation among miRNAs, lncRNAs, and mRNAs 33 , 34 . Consequently, a simultaneous increase in mRNA expression along with either lncRNA or circRNA implies the potential downregulation of a specific miRNA. To explore interactions as edges in the primary ceRNA network \(P^{C}\) among lncRNAs and mRNAs, we utilize the miRTarBase and TarBase databases to access experimentally verified mRNA–miRNA interactions. 35 , 36 We also extract experimentally verified mRNA–miRNA interactions from the miRTarBase database. To predict interactions as edges among the union of mRNAs and lncRNAs in the set \(N^{C}\), following the ceRNA regulatory model, we compute the PCC on the gene expression data of the RNAs, extracting their interactions from the miRTarBase and TarBase databases. In other words, assuming there is a documented interaction between miRNA A and both mRNA B and lncRNA C in the databases, we predict an edge between mRNA B and lncRNA C if the gene expression data shows a negative PCC between miRNA A and mRNA B, as well as between miRNA A and lncRNA C. Simultaneously, there should be a positive PCC between the gene expression data of mRNA B and lncRNA C. Therefore, in cases where lncRNAs and miRNAs share a common miRNA, as indicated in the mentioned databases, we forecast potential interaction between them. Due to the absence of experimental interaction data between circRNAs and miRNAs, we predict interactions among circRNAs or between circRNAs and mRNAs solely based on the PCC. In other words, to infer the interactions between circRNAs and mRNAs, we calculate the PCC for each pairing of circRNAs with miRNAs, mRNAs with miRNAs, and circRNAs with mRNAs. We predict an interaction between a circRNA and an mRNA when there is a positive PCC between their gene expression data, as well as a negative PCC between the gene expression data of the circRNA and a miRNA, as well as between the mRNA and that miRNA. For both the scenarios: negative PCC values between microRNAs with lncRNAs, mRNAs, and circRNAs, and positive PCC values among lncRNAs, mRNAs, and circRNAs pairings, a stringent threshold of an absolute PCC surpassing 0.4 and a p value beneath 0.05 is instated for statistical robustness. Meanwhile, we remove RNAs from the set \(N^{C}\) if they lack edges connecting them to other RNAs. In culmination, the visual representation of the primary ceRNA network \(P^{C}\) for each \(C \in CBCS\) can be achieved via Cytoscape 37 . Specific ceRNA networks In the third step of pipeline, we aim to establish a specific ceRNA network for each breast cancer subtype to delve into the distinctive regulatory mechanisms inherent to each subtype. To achieve this, we first extract a common ceRNA network from the primary ceRNA networks. Subsequently, we derive the specific ceRNA network for each subtype, C, by subtracting the common ceRNA network from the primary ceRNA network \(P^{C} = \left( {N^{C} ,E^{C} } \right)\). To achieve this, we start by constructing the common ceRNA network called \(cCENT = \left( {N^{cCENT} ,E^{cCENT} } \right)\) by amalgamating the primary ceRNA networks, \(P^{basal{-}like}\), \(P^{luminal\;B}\), \(P^{normal{-}like}\), \(P^{luminal\;A}\) and \(P^{HER2 + }\) where $$ \begin{aligned} N^{cCENT} & = N^{basal{-}like} \cap N^{luminal\;B} \cap N^{normal{-}like} \cap N^{luminal\;A} \cap { }N^{HER2 + } , \\ E^{cCENT} & = E^{basal{-}like} \cap E^{luminal\;B} \cap E^{normal{-}like} \cap E^{luminal\;A} \cap { }E^{HER2 + } . \\ \end{aligned} $$ Then, a specific ceRNA network named \(sCENT^{C} = \left( {N_{sCENT}^{C} ,E_{sCENT}^{C} } \right)\) is defined for each breast cancer subtype \(C \in BCS\) as bellow: $$ \begin{aligned} N_{sCENT}^{C} & = N^{C} - N^{cCENT} , \\ E_{sCENT}^{C} & = E^{C} - E^{cCENT} , \\ \end{aligned} $$ where \(N^{C}\) and \(E^{C}\) represent the number nodes (RNAs) and edges (interactions) among RNAs in the corresponding primary ceRNA network of breast cancer subtype C, respectively. Module detection In the fourth step of the Pre_CLM_BCS pipeline, we employ ClusterMaker Cytoscape application 38 to detect modules within the specific ceRNA networks. This application leverages the Markov cluster algorithm 39 . For each breast cancer subtype \(C \in BCS\), the specific ceRNA network \(sCENT^{C}\) is input into the application to identify modules. We select a collection of modules named \(LC - MODUL^{C}\) that contains one lncRNA or circRNA for further analysis in breast cancer subtype \(C\). Enrichment analysis In the fifth step of the pipeline, potential Gene Ontology (GO) terms of the modules in the set \(LC - MODUL^{C}\) of each breast cancer subtype, \(C \in BCS\), are determined through GO enrichment analysis conducted using the ToppFun web tool. GO terms with FDR-corrected p values less than 0.05 are considered significant 40 . Biomarker detection In the sixth step of our pipeline, we apply two methodologies to extract biomarkers from the set \(LC\_Module^{C}\), focusing on breast cancer subtype \( C\). In the first methodology, we utilize survival analysis to identify effective mRNAs in breast cancer subtype \(C\). In the second one, we employ differential gene expression data to predict biomarkers by comparing breast cancer subtype \(C\) to the other subtypes separately. Survival analysis For each breast cancer subtype, denoted as \(C \in BCS\), we calculate the union of mRNA genes from the modules listed in \(LC\_Module^{C}\). Subsequently, we explore the relationship between the expression of these mRNA genes and patients' overall survival by employing the Kaplan–Meier (KM) estimation and Log-rank test. TCGA survival data matching the mRNA genes are retrieved from TCGA database for survival analysis. A total of 88, 100, 76, 252, and 42 samples are collected for basal-like, luminal B, normal-like, luminal A, and HER2+ subtypes, respectively. For each gene’s survival analysis, patients are divided into two groups based on their expression values. Those with values exceeding the mean are categorized as the high-expression group, while those with values below the mean are classified as the low-expression group. We compute Kaplan–Meier curves and Log-rank test statistical results for both groups using the Survminer 41 and survival 42 packages in R. Differential gene expression analysis For the analysis of differentially expressed genes (DEG), we create the set of genes named \(G^{C}\) from \(LC\_Module^{C} \) for each breast cancer subtype \(C \in BCS\). Then, we conduct DEG analysis for each pair of subtypes \(C,C^{\prime} \in BCS\), where \(C \ne C^{\prime},\) as the between-subtype differentially expressed genes, denoted as \(BSDEG^{C,C^{\prime}}\). This analysis is performed using the edgeR package on the gene expression data of the genes in \(G^{C}\) 43 . Each gene with adjusted p values less than 0.05 (Benjamini–Hochberg) and an absolute logFC value greater than 0.5 is included in the set \(BSDEG^{C,C^{\prime}}\). Subsequently, we extract the set of biomarkers (\(Bio\_DEG^{C}\)) for breast subtype \(C\), as follows: $$ Bio\_DEG^{C} = \bigcap\limits_{\begin{subarray}{l} C^{\prime} \in BCS \\ C^{\prime} \ne C \end{subarray} } {BSDEG^{{C,C^{\prime}}} } . $$ To validate the predicted biomarkers in \(Bio\_DEG^{C} { }\) for each cancer subtype \(C\), we employ a support vector machine (SVM) classifier, and implement fivefold cross-validation. The validation procedures using the SVM classifier and fivefold cross-validation are performed individually for each subtype against all other subtypes. This process is repeated for each of the five subtypes. After running SVM on the breast cancer subtype \( C \in BCS\), we calculate the receiver operating characteristic (ROC) curve and the area under the curve (AUC) to evaluate the efficiency of the biomarkers, \(Bio\_DEG^{C}\), in distinguishing the breast cancer subtype from the others. The caret 44 , plotROC 45 , MLmetrics 46 , and ggplot 47 R packages are utilized to implement the SVM algorithm and to calculate and illustrate the ROC curves. Results In this section, we offer a step-by-step representation of the results obtained from executing the pipeline outlined in Fig. 1 for biomarker detection to distinguish breast cancer subtypes. Dataset Figure  2 illustrates the outcomes of the initial step in the Pre_CLM_BCS pipeline. During this phase, we extract gene expression of 30,480 genes, encompassing mRNA, miRNA, circRNA, and lncRNA. After the completion of the first and second data preprocessing stages within this step, we narrow down the initial gene pool to 14,056 genes. The figure also breaks down the retained genes by their respective types. Figure 2 The number of genes and samples for breast cancer subtypes after preprocessing: basal-like, luminal B, normal-like, luminal B, and HER2+. Note: The picture was created using biorender software (biorender.com). The validation results of the model on each breast cancer subtype demonstrate satisfactory performance. The average AUC values of fivefold for basal-like, HER2+, luminal A, luminal B, and normal-like are obtained 0.93, 0.92, 0.70, 0.82, and 0.86, respectively. Overall, these results indicate that the selected genes in the set \(Bio\_DEG^{C}\) for each subtype \(C\) serve as robust biomarker sets for distinguishing between breast cancer subtypes. Among the biomarkers predicted in this stage, circRNAs and lncRNAs emerge as promising candidates for breast cancer subtype prediction. Discussion In the realm of immuno-oncology, a field at the intersection of immunology and cancer research, the development of potent immunotherapy treatments holds promise for various cancer types 49 . Understanding the immune system’s intricacies is paramount to advancing these treatments. Breast cancer research has underscored the pivotal role of T cells in tumor progression, across various subtypes 50 , 51 , 52 . This section investigates T cell regulation and activation mechanisms within these subtypes, revealing at least one module in each subtype enriched in these processes. as follows: In the HER2+ subtype, module C in Fig. 9 comprises ZDHHC14 and CROT mRNA genes, along with the lncRNA RGS5-AS1. This module is associated with the positive activation of T cells, which is influenced by ZDHHC14 and CROT mRNA genes. These genes are subject to regulation by the lncRNA RGS5-AS1 through hsa-miR-20a-5p. Additionally, it’s worth noting that ZDHHC14 is known as a tumor suppressor gene with relevance across various cancer types 53 . In the luminal B subtype, modules G and H in Fig. 6 are pivotal in regulating T cell proliferation and activation, while module C is associated with nipple morphogenesis and influences the T cell apoptotic process. These modules reveal the presence of hsa-miR-149-3p and hsa-miR-133a-3p miRNAs, which are shared between lncRNA C6orf223 and the mRNA molecules in module G, and between lncRNA HHLA3-AS1 and the mRNA molecules in module H. This discovery underscores the potential of these miRNAs as influential genes involved in suppressing various oncogenic pathways 54 . As a result of our analysis, we suggest that biologists consider lncRNAs C6orf223 (in module G) and HHLA3-AS1 (in module H), as well as circ_000349 (in module C), as potential biomarkers for the luminal B subtype. In the basal-like subtype, module A (see Fig. 5 ) includes BCL2L12 mRNA linked to T Helper-2 (Th2) activation processes. This gene’s role in modulating immune responses suggests it as a potential therapeutic target 55 . CircRNAs, circ_003996, circ_003925, circ_004109, and circ_000736, in this module further regulate BCL2L12 via hsa-miR-539-5p. In the luminal A subtype, circ_002665 (module C in Fig. 8 ) and circ_001855 (module D in Fig. 8 ) appear to play essential roles in regulating T-cell activation mechanisms. These circRNAs are involved in the regulation of HOXB6 in Module C and HLA-F and CAMKV in Module D. Previous studies have discussed the roles of HOXB6 in breast cancer 56 , and the dysregulation of HLA-F and CAMKV has been reported in various cancers 57 , 58 , 59 . This analysis underscores the precise involvement of HOXB6, HLA-F, and CAMKV in tumor progression through T cell activation processes in the luminal A subtype. In the normal-like subtype, module B (Fig. 7 ), which is enriched in T Helper Cell Surface Molecules, includes DLEU1. DLEU1 has been identified as an up-regulated lncRNA in breast cancer tissues and cells, particularly in tumors with high malignancy 60 , 61 . This lncRNA interacts with NPHS1 mRNA. Based on our analysis, NPHS1 exhibits significant prognostic value and shows differential expression between normal-like subtype and the other four breast cancer subtypes. This underscores NPHS1 as a potential dependable immune molecular marker for normal-like subtype. Interestingly, three shared miRNAs, namely hsa-miR-378a-5p, hsa-miR-125a-3p, and hsa-miR-150-5p, are found to exist between DLEU1 and NPHS1. These miRNAs are known as anti-apoptotic agents in breast cancer, further highlighting their potential roles in immune regulation within the normal-like subtype 62 , 63 , 64 . These findings emphasize the complex interplay between various RNA molecules and immune mechanisms across breast cancer subtypes, offering insights into potential therapeutic targets and diagnostic markers. Conclusion In conclusion, breast cancer poses a significant global health challenge. Recent research highlights the critical roles of ncRNAs, specifically lncRNAs and circRNAs, in both cancer and immune system functions. The competing endogenous RNA hypothesis provides a valuable framework to understand the interplay between mRNAs and ncRNAs in the context of miRNA regulation. This study employed an innovative pipeline named Pre_CLM_BCS to investigate the roles of ncRNAs in various breast cancer subtypes. By constructing subtype-specific ceRNA networks, we unveiled distinct modules associated with immune and cancer-related processes within each subtype. Table 5 reveals the results of our pipeline, which has allowed us to identify several potential biomarkers specific to different breast cancer subtypes. Additionally, our analysis has unveiled certain ncRNAs that appear to play pivotal roles in the progression of their respective breast cancer subtypes, possibly through their regulatory effects on T cell mechanisms. Table 5 A compilation of potential biomarkers identified by our pipeline. In the realm of immuno-oncology, understanding T cell regulation and activation mechanisms is essential. This study revealed the involvement of specific RNA molecules in these processes across different breast cancer subtypes. These findings provide insights into potential therapeutic targets and diagnostic markers, shedding light on the intricate interplay between various RNA molecules and immune mechanisms in the context of breast cancer subtypes. This research contributes to the ongoing efforts to improve breast cancer diagnosis and treatment, ultimately reducing the global burden of this disease. While we have developed a pipeline to investigate the competitive regulatory mechanisms of circRNAs and lncRNAs within the context of miRNAs competition for the detection of breast cancer subtypes, it’s evident that each step of the pipeline has potential for improvement. In our future research, we aim to explore the implementation of graph neural networks (GNNs) 65 to enhance the construction of specific ceRNA networks in the third step of the pipeline. GNNs have shown promise in improving the accuracy of network analysis and prediction, and we anticipate that leveraging this technology will lead to the discovery of more effective biomarkers. Furthermore, we plan to investigate the integration of ceRNA networks with gene-protein signaling networks 66 , 67 , 68 . This combination promises to provide a more holistic view of the regulatory interactions within a biological system. By unifying coding and ncRNA interactions with protein signaling pathways, we can gain a deeper understanding of how these components mutually influence each other’s functions, ultimately shedding light on complex biological processes and diseases. This integrated approach holds the potential to identify key regulatory nodes, uncover novel biomarkers, and elucidate the intricate interplay between ceRNA networks and protein signaling. Data availability The study utilized miRNA-seq and RNA-SeqV2 level 3 data, which encompass mRNA and lncRNA expression data of breast cancer, sourced from TCGA Research Network ( http://cancergenome.nih.gov/ ) using the National Cancer Institute (NCI) Genomic Data Commons (GDC) resource ( https://gdc.cancer.gov/ ). CircRNA expression data for breast cancer was acquired based on the information provided by Asha et al. ( https://doi.org/10.18632/oncotarget.13134 ) 28 . References Siegel, R. L., Miller, K. D., Fuchs, H. E. & Jemal, A. Cancer statistics, 2021. CA Cancer J. Clin. 71, 7–33 (2021). Rights and permissions Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/ .