Expression profiles and functional prediction of long non-coding RNAs LINC01133, ZEB1-AS1 and ABHD11-AS1 in the luminal subtype of breast cancer
Journal of Translational Medicine volume 19, Article number: 364 (2021)
Luminal breast cancer (BC) is the most frequent subtype accounting for more than 70% of BC. LncRNAs, a class of non-coding RNAs with more than 200 nucleotides, are involved in a variety of cellular processes and biological functions. Abberant expression is related to the development of various cancers, such as breast cancer. LINC01133, ZEB1-AS1, and ABHD11-AS1 were reported to be dysregulated in different cancers. However, their expression level in luminal BC remains poorly known. The aim of the present study was to evaluate the potential roles of these lncRNAs in BC, especially in luminal subtypes.
A comprehensive analysis was performed using the Lnc2Cancer database to identify novel cancer-associated lncRNA candidates. After conducting a literature review, three novel lncRNAs named LINC01133, ZEB1-AS1, and ABHD11-AS1 were chosen as target genes of the present study. Quantitative real‐time polymerase chain reaction (qRT-PCR) was used to evaluate the expression level of the mentioned lncRNAs in both luminal BC tissues and cell lines. Then, the correlation of the three mentioned lncRNAs expression with clinicopathological characteristics of the patients was studied. Moreover, several datasets were used to discover the potential roles and functions of LINC01133, ZEB1-AS1 and ABHD11-AS1 in luminal subtype of BC.
According to the qRT-PCR assay, the expression levels of LINC01133 and ZEB1-AS1 were decreased in luminal BC tissues and cell lines. On the other hand, ABHD11-AS1 was upregulated in the above-mentioned samples. The expression levels of LINC01133, ZEB1-AS1, and ABHD11-AS1 were not associated with any of the clinical features. Also, the results obtained from the bioinformatics analyses were consistent with qRT-PCR data. Functional annotation of the co-expressed genes with the target lncRNAs, protein–protein interactions and significantly enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways across luminal BC were also obtained using bioinformatics analysis.
Taken together, our findings disclosed the dysregulation of LINC01133, ZEB1-AS1, and ABHD11-AS1 in luminal BC. It was revealed that LINC01133 and ZEB1-AS1 expression was significantly downregulated in luminal BC tissues and cell lines, while ABHD11-AS1 was upregulated considerably in the mentioned tissues and cell lines. Also, bioinformatics and systems biology analyses have helped to identify the possible role of these lncRNAs in luminal BC. However, further analysis is needed to confirm the current findings.
Breast cancer (BC) is the most common cancer among women worldwide, accounting for 30% of all female cancers alone [1, 2]. It is also the second leading cause of cancer death with 11.6% of the total cancer deaths . The luminal subtype of BC, expressing both estrogen and progesterone receptors (ER+/PR+), accounts for a large percentage of this cancer (more than 70%) [4, 5]. Luminal type of BC is classified into two groups, luminal A and luminal B, according to human epidermal growth factor receptor 2 (HER2) status and levels of ki-67 . The leading choice for luminal BC treatment is endocrine therapy. However, resistance to endocrine therapy is a major challenge for clinicians . Also, similar treatment strategies can provide a variety of responses in luminal subtype patients . In general, despite the better prognosis, due to the heterogeneity and resistance to hormone therapy of luminal subtype, the treatment effectiveness is still limited . Therefore, identifying novel biomarkers and molecular mechanisms of BC is required to perform early diagnosis and personalized therapeutic strategies. In addition, the limited number of studies on the dysregulation of long non-coding RNAs (lncRNAs) in luminal BC, as well as the high frequency of this subtype, increase the importance of further research on this subgroup.
The GENCODE project results have indicated that only about 2% of the human genome encodes proteins, whereas the vast majority of it is transcribed into non-coding RNAs . LncRNAs, a class of non-coding RNAs with more than 200 nucleotides, have been shown to play essential roles in many different cellular processes, such as genetic imprinting, transcriptional responses, development, etc. [9, 10]. Therefore, the aberrant expression of lncRNAs can lead to the pathogenesis of various diseases including cancers . They have been identified as crucial molecules involved in cancer proliferation, invasion, and resistance to therapy [12, 13]. Some lncRNAs are expressed highly specific in different cancer types, and some of them are associated with the survival of patients. These two features make them ideal prognostic and diagnostic tools . Due to the development of lncRNA-based therapeutic approaches, determining the role of different lncRNAs in the tumorigenic process is of great importance. Hence, in the present study, we aimed to figure out the expression level and the underlying molecular mechanisms of LINC01133, ZEB1-AS1, and ABHD11-AS1 in luminal BC.
LINC01133, a novel lncRNA, is located on chromosome 1q23.2. It was reported that LINC01133 exerted a tumor suppressive role in gastric cancer , ovarian cancer , and esophageal squamous cell carcinoma  progression. It was also found that LINC01133 played oncogenic role in the development of some other cancers [18,19,20]. These dual roles of LINC01133 could attribute to the tissue-specific expression of this lncRNA. ZEB1-AS1 is associated with the progression and development of several cancers. Its aberrant expression was detected in many tumors , such as gastric cancer , hepatocellular carcinoma , and glioma . ABHD11-AS1, which is located at 7q11.23, is a newly discovered lncRNA. Previous studies have reported the upregulation of this lncRNA in different tumors, including papillary thyroid carcinoma  and pancreatic cancer . To the best of our knowledge, the present research is the first study that determines the dysregulation of ABHD11-AS1 in BC so far. Despite all the studies that have been done on these lncRNAs, their biological processes and potential mechanisms, specifically in luminal subtype of BC, are still unclear.
In the present study, we examined the expression level of three lncRNAs in luminal BC by qRT-PCR assay. Furthermore, we performed a correlation analysis between the expression of these lncRNAs and clinicopathological parameters of patients. Finally, to investigate the potential role and mechanism of the lncRNAs in luminal BC, different bioinformatics and systems biology analyses were accomplished using various databases.
Materials and methods
Identification of novel lncRNAs and performing literature mining
Lnc2Cancer database was used to discover novel lncRNAs that are associated with different types of cancers. Lnc2Cancer is a comprehensive database which provides lncRNA-cancer correlations between 165 different cancer types and 1614 lncRNAs [27, 28]. Then, to find lncRNAs that have been less or not studied in breast cancer, a comprehensive PubMed data mining was performed on the candidate lncRNAs. The following keywords were searched: “long non-coding RNA” or “lncRNA”, “breast cancer” or “breast carcinoma”, and “the candidate lncRNAs”. LncRNAs whose expression and role in luminal BC were not studied were selected. Also, among the selected lncRNAs, those that played an important role in several cancer-associated signaling pathways were included in the study.
Breast cancer tissue samples
To perform the present case–control study, luminal A and B BC tissue samples and paired adjacent non-cancerous tissues were obtained from 79 patients at Breast Cancer Research Biobank (BCRC-BB) . BCRC-BB follows every ethical guideline about storing and utilizing human biological samples for biobanks. Tissue specimens were snap-frozen in liquid nitrogen and stored at − 80 °C until performing tests. The inclusion criterion was that patients should not receive any local or systemic treatment before surgery. Clinicopathological parameters including age at diagnosis, tumor size, tumor subtype, grade, stage, pathology of tumors, and ER/PR/HER2 status were retrieved from the hospital records. Written consent was signed by all of the patients, and all the samples were approved by relevant pathologists. Also, the Ethics committee of Tehran University of Medical Sciences (TUMS) approved the present research (Code: IR.TUMS.MEDICINE.REC.1398.792).
Breast cancer cell lines and culture
Luminal cell lines, MCF7 and T47D, and a non-invasive breast epithelial cell line, MCF10A, were used for the current project. MCF7 and T47D were grown in DMEM (Sigma-Aldrich, St. Louis, MO, USA), 10% fetal bovine serum (Gibco, Carlsbad, CA, USA), 100 U/ml penicillin, and 100 μg/ml streptomycin (Sigma-Aldrich, St. Louis, MO, USA). MCF10A were maintained in DMEM supplemented with 5% fetal bovine serum, epidermal growth factor (EGF) (20 ng/ml), insulin (10 µg/ml), hydrocortisone (0.5 mg/ml), cholera toxin (100 ng/ml), penicillin (100 U/ml), and streptomycin (100 μg/ml). An incubator at 37 °C with 5% CO2 and 95% humidity was required for maintaining cells.
Isolation of total RNA and qRT-PCR
Total RNA isolation from breast tissues and cultured cell lines was performed by RiboEx™ (GeneAll), applying the manufacturer’s instructions. Then, complementary DNA (cDNA) synthesis was carried out by performing reverse transcription process of 1 µg of isolated RNA samples using 5× All-In-One RT MasterMix kit (Applied Biological Materials). Next, q-RT PCR was carried out by using AMPLIQON Real Q Plus 2× Master Mix Green Low ROX on LightCycler96 Roche system. Moreover, using serially diluted samples, the efficiencies for each pair of primers of the genes of interest were calculated. The expression levels of the genes of interest were normalized to the expression level of Beta-2-Microglobulin (β2M) as an internal control. The primer sequences are listed in Table 1.
Several web servers and databases were used to select the genes of the current study and their role in luminal BC (Table 2).
Gene correlation analysis
While lncRNAs have constituted a rapidly expanding field of research, only a small number of them has been identified functionally thus far. Analyzing the mechanisms of co-expressed genes with lncRNAs can be very useful in predicting the potential function of novel RNAs. For this purpose, the UCSC Xena Browser  was used to select a list of luminal samples. The following parameters were included: “Luminal A”, “Luminal B”, and “The TCGA Sample ID”. Then, the TCGA ID of the specimens was used to search the FPKM file of the selected samples in the Genomic Data Commons (GDC) database . Finally, a correlation analysis was performed to obtain the co-expressed genes with LINC01133, ZEB1-AS1, and ABHD11-AS1 across the selected samples. Pearson correlation coefficient for the co-expressed genes has been computed by a standard method (R ≥ 0.4).
Expression analysis among molecular subtypes of BRCA
The expression pattern of the lncRNAs LINC01133, ZEB1-AS1 and ABHD11-AS1 was further evaluated using TANRIC database  across TCGA Breast invasive carcinoma (BRCA) dataset. The analysis was performed among different BC subtypes of the PAM50 classification (Normal-like, Luminal A, Luminal B, Basal and HER2+). TANRIC database has characterized the expression pattern of several lncRNAs using TCGA and different independent datasets (> 8000 samples overall). In addition, the expression level of the mentioned lncRNAs was evaluated among luminal and also non-luminal cell lines using GENEVESTIGATOR database . GENEVESTIGATOR integrates thousands of RNA-Seq and microarrays data in order to evaluate gene expression among several cancers.
The International Cancer Genome Consortium (ICGC) portal  was used to identify genomic alterations of the selected lncRNAs across luminal BC. The mutational data was obtained from BRCA-EU project including ER+ HER2− tumor subtypes (luminal A and B groups). Also, figures of the lncRNAs were illustrated by trackViewer R Package  according to the location of the mutations.
Functional enrichment analysis
To investigate the associated pathways and functions of LIN01133, ZEB1-AS1 and ABHD11-AS1, functional enrichment analysis of the genes co-expressed with these three lncRNAs was performed using various databases. First, the Database for Annotation, Visualization and Integrated Discovery (DAVID) [35, 36] was used for Gene Ontology (GO) term enrichment analysis. To this end, the list of co-expressed genes with LINC01133, ZEB1-AS1 and ABHD11-AS1 was used which was obtained in the previous step. Then, to remove redundant GO terms based on similar measures, the list of GO terms obtained from DAVID was submitted to the REVIGO . At this level, Cytoscape software  was utilized for performing a detailed visualization of the networks specified by REViGO. The enrichr database  was also used to evaluate pathways in which three lncRNAs may be involved. Calculating and visualizing protein–protein interactions (PPI) of the co-expressed genes with LINC01133, ZEB1-AS1 and ABHD11-AS1 were done using the STRING database . Also, the hub genes of the lncRNAs co-expressed genes were explored by the degree method available in CytoHubba , a Cytoscape application.
IBM SPSS Statistics (Statistical Package for the Social Sciences, Version 24.0) was used for analyzing the data obtained from the qRT-PCR assay. The tumoral and adjacent non-tumoral tissues were extracted from paired sources. The differences between ∆Ct values of LINC01133, ZEB1-AS1, and ABHD11-AS1 expression levels in tumor and adjacent non-tumor tissue specimens were analyzed using paired sample t-test. For clinicopathological correlation analysis, the median expression of the lncRNAs was used for dividing 79 patients into two categories: relatively high and relatively low expression of the lncRNAs. Also, correlation analysis of LINC01133, ZEB1-AS1, ABHD11-AS1 and clinicopathological features was carried out by the chi-square and independent t-test. Moreover, Kaplan–Meier method and log-rank test were used to calculate the overall survival rates and compare differences between survival curves, respectively. The data is represented as mean ± standard deviation (SD). P-values of less than 0.05 were considered significant. The diagnostic value of each gene was calculated via the receiver operating characteristic (ROC) curve analysis. ROC curves were illustrated by IBM SPSS Statistics version 24.0.
Significantly dysregulated expression of LINC01133, ZEB1-AS1 and ABHD11-AS1 in luminal BC tissues and cell lines
According to the results of the qRT-PCR assay, LINC01133 and ZEB1-AS1 were significantly downregulated in luminal A and B BC tissues (pLINC01133 = 0.001) (pZEB1-AS1(T5,6) < 0.001) (pZEB1-AS1(T1-4) = 0.002) (Fig. 1A, C, E). Also, ABHD11-AS1 expression was upregulated considerably in luminal BC tissues (p < 0.001, Fig. 1G). As shown in Fig. 1B, D, F), LINC01133, ZEB1-AS1(T5,6), and ZEB1-AS1(T1-4) were downregulated in 64%, 80% and 71% of samples, respectively. Moreover, ABHD11-AS1 was upregulated in 70% of cases (Fig. 1H).
The results obtained from two luminal BC cell lines, T47D and MCF7, were consistent with the tissue expression data. The expression of LINC01133 and ZEB1-AS1 was decreased in the mentioned luminal cell lines compared with normal BC cell line (Fig. 2A–C). ABHD11-AS1 expression was upregulated significantly in luminal cell lines (Fig. 2D).
Correlation between LINC01133, ZEB1-AS1, and ABHD11-AS1 expression levels and the clinicopathological features
As shown in Tables 3 and 4, no significant correlation was detected between LINC01133, ZEB1-AS1(T1-4), ZEB1-AS1(T5,6), and ABHD11-AS1 expression with clinicopathological features. Also, according to Kaplan–Meier analysis, the expression of the lncRNAs was not significantly associated with overall survival of the patients (Additional file 1).
Evaluation of diagnostic value of LINC01133, ZEB1-AS1 and ABHD11-AS1 in luminal BC
The outcome of ROC analysis was defined based on the type of tissue (tumor or adjacent) and the predictor was characterized as gene expression level (upper or below the median expression). According to the ROC analysis of ZEB1-AS1(T5,6) shown in Fig. 3, the area under ROC curve is 0.729 across luminal A and B subtypes of BC (pZEB1-AS1(T5,6) < 0.001). On the other hand, LINC01133, ZEB1-AS1(T1-4), and ABHD11-AS1 were not suggested as good biomarkers for diagnosis of luminal subtypes of BC (pLINC01133 = 0.001/AUCLINC01133 = 0.661) (pZEB1-AS1(T1-4) = 0.005/AUCZEB1-AS1(T1-4) = 0.629) (pABHD11-AS1 = 0.009/AUCABHD11-AS1 = 0.375).
Expression pattern of lncRNAs LINC01133, ZEB1-AS1 and ABHD11-AS1 among PAM50-breast cancer subtype classifications
According to the differential expression analysis, lower expression of LINC01133 was observed in luminal A and B BC subtypes. The expression of ZEB1-AS1 in luminal B and HER2 subtypes was lower than other groups. Moreover, significant differences of the expression of these lncRNAs were observed among different PAM50 subtypes (pLINC01133 = 1.13e−32) (pZEB1-AS1 = 9.82e−5). On the other hand, luminal A and B had higher expression of ABHD11-AS1 followed by Normal-like, HER2, and Basal subtypes. The expression of ABHD11-AS1 was also significantly different among molecular subtypes of BC (pABHD11-AS1 = 2.79e−8) (Additional file 2).
Besides, results obtained from the GENEVESTIGATOR database across luminal BC cell lines were consistent with the experimental data of the present study which indicated that LINC01133 and ZEB1-AS1 were downregulated in luminal cell lines. Also, expression analysis of these lncRNAs in several non-luminal BC cell lines showed that their expression was also downregulated across these cell lines (Additional files 3, 4). However, as ABHD11-AS1 is a novel lncRNA, no information was available about this lncRNA across GENEVESTIGATOR database.
Genetic alterations of LINC01133, ZEB1-AS1 and ABHD11-AS1 across luminal BC
According to the results obtained from the ICGC data portal, LINC01133 mutations in ER+ HER2– BC take place most frequently in intronic, downstream, and upstream regions, respectively. Also, most LINC01133 mutations have been identified to be of the substitution type (Fig. 4A).
It was also found that ZEB1-AS1 mutations occur more frequently on introns and exons and are of the substitution type (Fig. 4B).
Most of ABHD11-AS1 mutations are substitutions which mostly take place in upstream and downstream regions (Fig. 4C).
Co-expression gene network and functional annotation analysis
A list containing 420 luminal samples (230 luminal A and 190 luminal B) was obtained using UCSC Xena Browser and the Genomic Data Commons (GDC) databases (Additional file 5). Results obtained from the co-expression analysis of LINC01133, ZEB1-AS1, and ABHD11-AS1 across the selected luminal A and B BC samples suggested that these lncRNAs are significantly co-expressed with 256, 149 and 126 genes, respectively (R ≥ 0.4) (Additional file 6).
Data retrieved from GO term enrichment analysis by DAVID included a list for categories of biological process (BP), cellular component (CC), and molecular function (MF) (Additional file 7). Also, results of summarizing and excluding redundant GO terms of the target lncRNAs of the present study, which were performed by REVIGO, were concordant with the data obtained from DAVID GO term enrichment analysis. The top 5 BP, CC, and MF terms related to each lncRNA are shown in Fig. 5. According to it, the genes co-expressed with LINC01133 are significantly involved in biological processes such as ion transport, response to external stimulus, ion transmembrane transport, positive regulation of phosphorylation, etc. It was also determined that molecular function of these genes considerably consists of ion/cation/metal binding, cation transmembrane transporter activity, etc. Each of the above activities takes place inside and/or in the vicinity of a cell. According to CC terms, it was predicted that the co-expressed genes with LINC01133 are considerably found in extracellular region/exosome, membrane-bounded vesicle and so on (p < 0.05) (Fig. 5A).
It has been found that the co-expressed genes with ZEB1-AS1 are considerably involved in the regulation of RNA metabolic process, regulation of nucleobase-containing compound metabolic process, nucleic acid-templated transcription, etc. As shown in Fig. 5B, the most significant molecular functions of these genes are signal transducer activity, kinase activity, phosphotransferase activity, etc. These genes are mostly enriched in the synapse and Golgi lumen (p < 0.05).
Moreover, the genes co-expressed with ABHD11-AS1 are considerably involved in biological processes such as transmembrane transport, reproductive process, ion transport, neurological system process, etc. They are also considerably involved in various molecular functions including signal transducer activity, transmembrane signaling receptor activity, molecular transducer activity, etc. Additionally, the GO cellular component terms suggested that the co-expressed genes with ABHD11-AS1 are considerably enriched in the intrinsic component of nuclear membrane part, apicolateral plasma membrane, etc. (p < 0.05) (Fig. 5C).
The results of gene set enrichment analysis by Enrichr are demonstrated in Fig. 6. According to KEGG, Reactome and WikiPathways analysis, the co-expressed genes with LINC01133 are mainly involved in sodium/proton exchanges, mucin-type O-glycan biosynthesis, and nuclear factor erythroid 2-related factor 2 (NRF2) pathway (Fig. 6A). Moreover, ZEB1-AS1 co-expressed genes are significantly enriched in tight junction interactions, estrogen biosynthesis, glycerophospholipid catabolism, and Wnt signaling pathway across luminal BC (Fig. 6B). The co-expressed genes with ABHD11-AS1 are mainly involved in the processes of Wnt ligand biogenesis, FGFR1 mutant receptor activation, mesodermal commitment pathway, and Wnt signaling (Fig. 6C).
Finally, the association analysis of the genes co-expressed with LINC01133 demonstrated that 42 out of 256 genes had a strong PPI (interaction score > 0.4) with each other (Fig. 7A). Also, according to PPI network analysis of the co-expressed genes with ZEB1-AS1 and ABHD11-AS1, 20 and 14 genes have been recognized with strong interactions with each other, respectively (Fig. 7B, C).
According to the CytoHubba, the genes KPRP, PPARG, ST6GALNAC1, and UGT1A1 were adopted as the top hub genes of the PPI network of LINC01133 co-expressed genes ranked by degree metric. Also, five (BAG4, ASH2L, DDHD2, LETM2, and WHSC1L1) and four (LACTB2, SPANXA1, SPANXA2, and SPANXC) genes were adopted as the hob genes related to ZEB1-AS1 and ABHD11-AS1 networks, respectively.
The luminal subtype of BC, which has shown better prognosis among the other subtypes, accounts for approximately two-thirds of this cancer [43, 44]. However, it is reported that the luminal B subtype has demonstrated poorer recurrence-free survival in adjuvant treatment categories in comparison to luminal A . Also, it was proved that a large percentage of ER+ patients with lymph node positive benefit less from adjuvant chemotherapy . Due to the heterogeneity of this cancer, the same therapies may demonstrate diverse outcomes in patients. So, distinguishing the underlying mechanisms and specific biomarkers can help the individualized treatment of especially high-risked luminal patients.
In this study, LINC01133 and ZEB1-AS1 were significantly downregulated in luminal A and B BC tissues compared to their adjacent non-tumoral tissues. Also, they were downregulated in luminal cell lines, namely T47D and MCF7, compared to a normal breast cell line (MCF10A). Moreover, ABHD11-AS1 was upregulated considerably in the mentioned tissues and cell lines.
LncRNAs have shown greater expression variability among different individuals and a greater degree of tissue-specificity, compared with coding-genes . The downregulation of LINC01133, ZEB1-AS1(T1-4), and ZEB1-AS1(T5,6) was detected in 64%, 71% and 80% of patients, respectively. Also, the upregulation of ABHD11-AS1 was found in a large proportion of patients (70% of cases). Song et al. reported that the downregulation of LINC01133 is considerably associated with the poor prognosis of BC patients . On the other hand, Luo et al. suggested that ZEB1‐AS1 promotes triple-negative breast cancer progression . So, this lncRNA may also play an important role in luminal BC patients, although further investigation is needed. Also, ZEB1-AS1(T5,6) may be used as a biomarker for luminal BC progression. However, more in-vitro tests are needed for more accurate conclusions.
Limited studies have been performed on the dysregulation of these lncRNAs and their roles in cancers so far. So, the combination of various datasets as well as the analysis of their co-expressed genes could provide better identification of functional roles of the mentioned lncRNAs in luminal BC. According to the result of TCGA data analysis available in TANRIC database, LINC01133, ZEB1-AS1, and ABHD11-AS1 expression were significantly different among PAM50 subtype classification. Besides, the expression pattern of LINC01133 and ZEB1-AS1 across luminal and non-luminal BC cell lines was consistent with the data from TANRIC database and also experimental data of the present study. According to the data obtained from the ICGC data portal, most of the mutations of LINC01133, ZEB1-AS1, and ABHD11-AS1 in ER+ HER2– patients are from substitution category. The effect of mutations and variations of lncRNAs have not been extensively studied as they constitute a novel class of RNAs. Ponjavic et al. demonstrated that functional lncRNAs show reduced single substitutions, deletions and insertions in their sequences. Although nucleotide substitutions occur mostly in protein-coding sequences compared with noncoding sequences, their effect on lncRNAs are so significant and should not be overlooked . Substitution mutations in lncRNAs may contribute to the pathogenesis of various cancers by disrupting the secondary structure of RNAs [18, 49]. The RNA secondary structure plays important roles in multiple cellular processes such as gene regulation and localization, splicing, stability and also translation . Also, mutations in the upstream regions of lncRNAs, which constitute a considerable percentage of LINC01133 and ABHD11-AS1 mutations, can have significant effects on them. For example, the presence of single-nucleotide polymorphisms in the promoter of a lncRNA may have an effect on its expression pattern . So, further studies are recommended to determine the effect of these mutations on the structure of the mentioned lncRNAs.
Remarkably, the protein-coding genes are significantly better annotated than lncRNAs . So, in this study, an attempt was made to determine the function of the target lncRNAs with the help of the known functional information of some coding and non-coding genes, using “guilt-by-association” principle. According to this principle, genes that are involved in some related and/or similar biological pathways may show similar expression patterns across different experimental conditions .
According to the “guilt-by-association” principle and the GO term enrichment analysis of the co-expressed genes with LINC01133, this lncRNA might be involved in ion transport, response to external stimulus, and positive regulation of phosphorylation, all of which are cancer‐associated biological processes. Also, ion/cation/metal binding and cation transmembrane transporter activity are known as the most crucial molecular functions in which LINC01133 is involved. LINC01133 might be mostly found in extracellular region/exosome and membrane-bounded vesicles.
According to the results obtained from the GO term enrichment analysis of ZEB1-AS1 co-expressed genes, this lncRNA is involved in the regulation of RNA metabolic process, regulation of nucleobase-containing compound metabolic process, nucleic acid-templated transcription, etc. Deregulation of all of these biological processes can lead to the development and progression of various cancers, including BC. In addition, ZEB1-AS1 might be plausibly involved in molecular functions including signal transducer activity, kinase activity and phosphotransferase activity, all of which can be linked to cancer progression. ZEB1-AS1 might be mostly enriched in synapse and Golgi lumen.
ABHD11-AS1 might be involved in transmembrane transport, reproductive process, ion transport, neurological system process, etc. ABHD11-AS1 also might be considerably involved in different molecular functions, like signal transducer activity, transmembrane signaling receptor activity, molecular transducer activity, etc. Dysregulation of all of the mentioned functions can lead to the development of various cancers. Moreover, ABHD11-AS1 is likely to be found in intrinsic component of nuclear membrane part, apicolateral plasma membrane, etc.
Significantly enriched KEGG, Reactome and WikiPathways of the co-expressed genes with LINC01133 include sodium/proton exchanges, mucin-type O-glycan biosynthesis, and NRF2 pathway. Abnormal glycoprotein structure of tumor cells can affect the survival, growth and metastasis of cancer cells . Nuclear receptors (NR) are known as regulators of physiological processes and play oncogenic or anti-oncogenic roles in cancerous cells . Growing evidence support the involvement of several NRs in the regulation of various pathways related to the initiation and development of BC . Overexpression of Nrf2 has been shown to increase the expression of glucose‐6‐phosphate dehydrogenase (G6PD) and Hypoxia‐inducing factor 1α (HIF‐1α) in BC cell lines, including MCF7. Also, overexpression of Nrf2 increases Notch1 expression via the G6PD/HIF-1α pathway. Notch signaling pathway regulates BC cell proliferation and migration by affecting the downstream gene, HES‐1, and the epithelial-to-mesenchymal transition (EMT) pathway, respectively . All of these evidences suggest that LINC01133 could play an essential role in BC.
Moreover, according to the pathway enrichment analysis, ZEB1-AS1 might be involved in tight junction interactions, estrogen biosynthesis, glycerophospholipid catabolism, and Wnt signaling pathway. All of these pathways can be associated with the progression and development of BC. Steroid hormones increase cell proliferation in BC . Also, alterations in tight junction complexes could facilitate BC initiation and progression by impairing their control over crucial cellular processes such as cell motility and polarity . Moreover, the Wnt pathway is significantly activated in breast tumors .
ABHD11-AS1 might also be mainly involved in Wnt ligand biogenesis, FGFR1 mutant receptor activation, mesodermal commitment pathway, and Wnt signaling pathway. The involvement of fibroblast growth factors (FGF) has been found in various cellular processes, including proliferation, drug resistance, anti-apoptosis and angiogenesis. Also, amplification of FGFR1 has been found in ER+ BC .
Data from STRING database showed that 42, 20, and 14 co-expressed genes with LINC01133, ZEB1-AS1 and ABHD11-AS1 have strong interactions with each other. These data confirmed the involvement of these lncRNAs and their co-expressed genes in similar biological pathways.
Taken together, the dysregulation of the three potential lncRNAs, LINC01133, ZEB1-AS1, and ABHD11-AS1 across luminal A and B subtypes of BC was reported in the present study. The bioinformatics analyses performed in this study helped us better identify the possible role of these lncRNAs in luminal BC. However, more experimental studies are needed to confirm these findings and to verify the exact roles of the mentioned lncRNAs.
All our findings have demonstrated the downregulation of LINC01133 and ZEB1-AS1 across luminal BC tissue specimens and cell lines compared to adjacent non-tumoral tissues and cell line. On the other hand, ABHD11-AS1 was upregulated in the mentioned BC tissues and cell lines.
Availability of data and materials
The datasets supporting the conclusions of this article are available in: [Lnc2Cancer database] at (https://www.bio-bigdata.com/lnc2cancer/). [UCSC Xena browser] at (http://xena.ucsc.edu/); [GDC database] at (https://portal.gdc.cancer.gov/). [TANRIC database] at (https://www.tanric.org); [GENEVESTIGATOR database] at (https://genevestigator.com); [International Cancer Genome Consortium] at (https://dcc.icgc.org). [DAVID database] at (https://david.ncifcrf.gov). [REVIGO database] at (https://revigo.irb.hr). [Cytoscape software] at (https://www.cytoscape.org). [STRING database] at (https://string-db.org). [Enrichr database] at (https://maayanlab.cloud/Enrichr/). Citation of all data is provided in the references list. The datasets supporting the conclusions of this article are also included within the article and its additional files.
Long non-coding RNAs
Quantitative real-time polymerase chain reaction
Kyoto Encyclopedia of Genes and Genomes
Human epidermal growth factor receptor 2
Breast Cancer Research Biobank
Tehran University of Medical Sciences
Epidermal growth factor
Polymerase chain reaction
Genomic Data Commons
Breast invasive carcinoma
The International Cancer Genome Consortium
The Database for Annotation, Visualization and Integrated Discovery
- IBM SPSS Statistics:
Statistical Package for the Social Sciences
Receiver operating characteristics
The International Cancer Genome Consortium
Nuclear factor erythroid 2-related factor 2
Hypoxia‐inducing factor 1α
Fibroblast growth factors
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70:7–30.
Zaidi Z, Dib HA. The worldwide female breast cancer incidence and survival, 2018. AACR. 2019. https://doi.org/10.1158/1538-7445.AM2019-4191.
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68:394–424.
Mathias C, Zambalde EP, Rask P, Gradia DF, de Oliveira JC. Long non-coding RNAs differential expression in breast cancer subtypes: what do we know? Clin Genet. 2019;95:558–68.
Makki J. Diversity of breast carcinoma: histological subtypes and clinical relevance. Clin Med Insights Pathol. 2015;8:CPath.S31563.
Abdel-Hafiz HA. Epigenetic mechanisms of tamoxifen resistance in luminal breast cancer. Diseases. 2017;5:16.
Yang F, Lv S-X, Lv L, Liu Y-H, Dong S-Y, Yao Z-H, Dai X-X, Zhang X-H, Wang O-C. Identification of lncRNA FAM83H-AS1 as a novel prognostic marker in luminal subtype breast cancer. OncoTargets Ther. 2016;9:7039.
Yang G, Shen T, Yi X, Zhang Z, Tang C, Wang L, Zhou Y, Zhou W. Crosstalk between long non-coding RNA s and Wnt/β-catenin signalling in cancer. J Cell Mol Med. 2018;22:2062–70.
Wang KC, Chang HY. Molecular mechanisms of long noncoding RNAs. Mol Cell. 2011;43:904–14.
Zarkou V, Galaras A, Giakountis A, Hatzis P. Crosstalk mechanisms between the WNT signaling pathway and long non-coding RNAs. Non-coding RNA Res. 2018;3:42–53.
Sanchez Calle A, Kawamura Y, Yamamoto Y, Takeshita F, Ochiya T. Emerging roles of long non-coding RNA in cancer. Cancer Sci. 2018;109:2093–100.
Sahu A, Singhal U, Chinnaiyan AM. Long noncoding RNAs in cancer: from function to translation. Trends Cancer. 2015;1:93–109.
Chandra Gupta S, Nandan TY. Potential of long non-coding RNAs in cancer patients: from biomarkers to therapeutic targets. Int J Cancer. 2017;140:1955–67.
Yarmishyn AA, Kurochkin IV. Long noncoding RNAs: a potential novel class of cancer biomarkers. Front Genet. 2015;6:145.
Yang X-Z, Cheng T-T, He Q-J, Lei Z-Y, Chi J, Tang Z, Liao Q-X, Zhang H, Zeng L-S, Cui S-Z. LINC01133 as ceRNA inhibits gastric cancer progression by sponging miR-106a-3p to regulate APC expression and the Wnt/β-catenin pathway. Mol Cancer. 2018;17:126.
Liu M, Shen C, Wang C. Long noncoding RNA LINC01133 confers tumor-suppressive functions in ovarian cancer by regulating leucine-rich repeat kinase 2 as an miR-205 sponge. Am J Pathol. 2019;189:2323–39.
Yang X-Z, He Q-J, Cheng T-T, Chi J, Lei Z-Y, Tang Z, Liao Q-X, Zhang H, Zeng L-S, Cui S-Z. Predictive value of LINC01133 for unfavorable prognosis was impacted by alcohol in esophageal squamous cell carcinoma. Cell Physiol Biochem. 2018;48:251–62.
Zheng YF, Zhang XY, Bu YZ. LINC01133 aggravates the progression of hepatocellular carcinoma by activating the PI3K/AKT pathway. J Cell Biochem. 2019;120:4172–9.
Huang C-S, Chu J, Zhu X-X, Li J-H, Huang X-T, Cai J-P, Zhao W, Yin X-Y. The C/EBPβ-LINC01133 axis promotes cell proliferation in pancreatic ductal adenocarcinoma through upregulation of CCNG1. Cancer Lett. 2018;421:63–72.
Zeng H-F, Qiu H-Y, Feng F-B. Long noncoding RNA LINC01133 functions as an miR-422a sponge to aggravate the tumorigenesis of human osteosarcoma. Oncol Res. 2018;26:335–43.
Ni X, Ding Y, Yuan H, Shao J, Yan Y, Guo R, Luan W, Xu M. Long non-coding RNA ZEB1-AS1 promotes colon adenocarcinoma malignant progression via miR-455-3p/PAK2 axis. Cell Prolif. 2020;53:e12723.
Zhang L-L, Zhang L-F, Guo X-H, Zhang D-Z, Yang F, Fan Y-Y. Downregulation of miR-335-5p by long noncoding RNA ZEB1-AS1 in gastric cancer promotes tumor proliferation and invasion. DNA Cell Biol. 2018;37:46–52.
Li T, Xie J, Shen C, Cheng D, Shi Y, Wu Z, Deng X, Chen H, Shen B, Peng C. Upregulation of long noncoding RNA ZEB1-AS1 promotes tumor metastasis and predicts poor prognosis in hepatocellular carcinoma. Oncogene. 2016;35:1575–84.
Lv Q-L, Hu L, Chen S-H, Sun B, Fu M-L, Qin C-Z, Qu Q, Wang G-H, He C-J, Zhou H-H. A long noncoding RNA ZEB1-AS1 promotes tumorigenesis and predicts poor prognosis in glioma. Int J Mol Sci. 2016;17:1431.
Wen J, Wang H, Dong T, Gan P, Fang H, Wu S, Li J, Zhang Y, Du R, Zhu Q. STAT3-induced upregulation of lncRNA ABHD11-AS1 promotes tumour progression in papillary thyroid carcinoma by regulating miR-1301-3p/STAT3 axis and PI3K/AKT signalling pathway. Cell Prolif. 2019;52:e12569.
Qiao X, Lv S, Qiao Y, Li Q, Ye B, Wang C, Miao L. Long noncoding RNA ABHD11-AS1 predicts the prognosis of pancreatic cancer patients and serves as a promoter by activating the PI3K-AKT pathway. Eur Rev Med Pharmacol Sci. 2018;22:8630–9.
Gao Y, Wang P, Wang Y, Ma X, Zhi H, Zhou D, Li X, Fang Y, Shen W, Xu Y. Lnc2Cancer v2.0: updated database of experimentally supported long non-coding RNAs in human cancers. Nucleic Acids Res. 2019;47:D1028–33.
Ning S, Zhang J, Wang P, Zhi H, Wang J, Liu Y, Gao Y, Guo M, Yue M, Wang L. Lnc2Cancer: a manually curated database of experimentally supported lncRNAs associated with various human cancers. Nucleic Acids Res. 2016;44:D980–5.
Majidzadeh-a K, Kaviani A, Esmaeili R, Farahmand L, Shojamoradi MH, Zare AA, Eini L, Abbasvandi F, Olfatbakhsh A, Moazen H. Iranian Breast Cancer Bio-Bank: the activity and challenging issues. Cell Tissue Bank. 2013;14:11–20.
Goldman MJ, Craft B, Hastie M, Repečka K, McDade F, Kamath A, Banerjee A, Luo Y, Rogers D, Brooks AN. Visualizing and interpreting cancer genomics data via the Xena platform. Nat Biotechnol. 2020;38:675–8.
Jensen MA, Ferretti V, Grossman RL, Staudt LM. The NCI Genomic Data Commons as an engine for precision medicine. Blood. 2017;130:453–9.
Li J, Han L, Roebuck P, Diao L, Liu L, Yuan Y, Weinstein JN, Liang H. TANRIC: an interactive open platform to explore the function of lncRNAs in cancer. Cancer Res. 2015;75:3728–37.
Hruz T, Laule O, Szabo G, Wessendorp F, Bleuler S, Oertle L, Widmayer P, Gruissem W, Zimmermann P. Genevestigator v3: a reference expression database for the meta-analysis of transcriptomes. Adv Bioinform. 2008. https://doi.org/10.1155/2008/420747.
Zhang J, Baran J, Cros A, Guberman JM, Haider S, Hsu J, Liang Y, Rivkin E, Wang J, Whitty B. International Cancer Genome Consortium Data Portal—a one-stop shop for cancer genomics data. Database. 2011. https://doi.org/10.1093/database/bar026.
Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13.
Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44.
Supek F, Bošnjak M, Škunca N, Šmuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS ONE. 2011;6:e21800.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Chin C-H, Chen S-H, Wu H-H, Ho C-W, Ko M-T, Lin C-Y. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;8:S11.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P. STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47:D607–13.
Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, Koplev S, Jenkins SL, Jagodnik KM, Lachmann A. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016;44:W90–7.
Ou J, Zhu LJ. trackViewer: a Bioconductor package for interactive and integrative visualization of multi-omics data. Nat Methods. 2019;16:453–4.
Chen X, Cong Y, Pan L, Jiang Y, Meng Q, Sun L, Pang H, Zhao Y, Dong X, Cai L. Luminal (Her2 negative) prognostic index and survival of breast cancer patients. Cancer Epidemiol. 2014;38:286–90.
Cheang MC, Chia SK, Voduc D, Gao D, Leung S, Snider J, Watson M, Davies S, Bernard PS, Parker JS. Ki67 index, HER2 status, and prognosis of patients with luminal B breast cancer. JNCI. 2009;101:736–50.
Kornienko AE, Dotter CP, Guenzl PM, Gisslinger H, Gisslinger B, Cleary C, Kralovics R, Pauler FM, Barlow DP. Long non-coding RNAs display higher natural expression variation than protein-coding genes in healthy humans. Genome Biol. 2016;17:14.
Song Z, Zhang X, Lin Y, Wei Y, Liang S, Dong C. LINC01133 inhibits breast cancer invasion and metastasis by negatively regulating SOX4 expression through EZH2. J Cell Mol Med. 2019;23:7554–65.
Luo N, Zhang K, Li X, Hu Y. ZEB1 induced-upregulation of long noncoding RNA ZEB1-AS1 facilitates the progression of triple negative breast cancer by binding with ELAVL1 to maintain the stability of ZEB1 mRNA. J Cell Biochem. 2020;121(10):4176–87.
Ponting CP, Oliver PL, Reik W. Evolution and functions of long noncoding RNAs. Cell. 2009;136:629–41.
Minotti L, Agnoletto C, Baldassari F, Corrà F, Volinia S. SNPs and somatic mutation on long non-coding RNA: new frontier in the cancer studies? High-throughput. 2018;7:34.
He F, Wei R, Zhou Z, Huang L, Wang Y, Tang J, Zou Y, Shi L, Gu X, Davis M. Integrative analysis of somatic mutations in non-coding regions altering RNA secondary structures in cancer genomes. Sci Rep. 2019;9:1–12.
Perron U, Provero P, Molineris I. In silico prediction of lncRNA function using tissue specific and evolutionary conserved expression. BMC Bioinform. 2017;18:144.
Wong DC, Sweetman C, Ford CM. Annotation of gene function in citrus using gene expression information and co-expression networks. BMC Plant Biol. 2014;14:186.
Brockhausen I. Mucin-type O-glycans in human colon and breast cancer: glycodynamics and functions. EMBO Rep. 2006;7:599–604.
Cheng HS, Lee JXT, Wahli W, Tan NS. Exploiting vulnerabilities of cancer by targeting nuclear receptors of stromal cells in tumor microenvironment. Mol Cancer. 2019;18:51.
Doan TB, Graham JD, Clarke CL. Emerging functional roles of nuclear receptors in breast cancer. J Mol Endocrinol. 2017;58:R169–90.
Zhang HS, Zhang ZG, Du GY, Sun HL, Liu HY, Zhou Z, Gou XM, Wu XH, Yu XY, Huang YH. Nrf2 promotes breast cancer cell migration via up-regulation of G6PD/HIF-1α/Notch1 axis. J Cell Mol Med. 2019;23:3451–63.
Capper CP, Rae JM, Auchus R. The metabolism, analysis, and targeting of steroid hormones in breast and prostate cancer. Horm Cancer. 2016;7:149–64.
Brennan K, Offiah G, McSherry EA, Hopkins AM. Tight junctions: a barrier to the initiation and progression of breast cancer? J Biomed Biotechnol. 2010. https://doi.org/10.1155/2010/46060.
Khramtsov AI, Khramtsova GF, Tretiakova M, Huo D, Olopade OI, Goss KH. Wnt/β-catenin pathway activation is enriched in basal-like breast cancers and predicts poor outcome. Am J Pathol. 2010;176:2911–20.
Katoh M, Nakagama H. FGF receptors: cancer biology and therapeutics. Med Res Rev. 2014;34:280–300.
The authors would like to appreciate Ms. Marzie Samimifar for proofreading the English language of the manuscript, Ms. Narges Jafarbeik-Iravani and Dr. Rasoul Abdollahzadeh for valuable assistance during this project. We also thank other members of Genetics department of Breast Cancer Research Center of Motamed Cancer Institute, ACECR, Tehran, Iran.
This project was supported by the Faculty of Medicine, Tehran University of Medical Sciences (TUMS). Our tissue samples were obtained from Breast Cancer Research Center of Motamed Cancer Institute (Tehran, Iran).
Ethics approval and consent to participate
The Ethics committee of Tehran University of Medical Sciences (TUMS) has approved the present research (Code: IR.TUMS.MEDICINE.REC.1398.792).
Consent for publication
The authors declare that there are no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The kaplan-meier survival functions of LINC01133, ZEB1-AS1(T1-4), ZEB1-AS1(T5,6), and ABHD11-AS1 across luminal BC samples.
Differential expression of LINC01133, ZEB1-AS1 and ABHD11-AS1 across different subtypes of BC, obtained by TANRIC database.
The expression level of LINC01133 across various BC cell lines, obtained by GENEVESTIGATOR database.
The expression level of ZEB1-AS1 across various BC cell lines, obtained by GENEVESTIGATOR database.
The list of 420 luminal samples obtained from UCSC Xena Browser and the Genomic Data Commons (GDC) database.
The list of co-expressed genes with LINC01133, ZEB1-AS1 and ABHD11-AS1 across the selected luminal A and B BC samples.
The list of biological process (BP), cellular component (CC), and molecular function (MF) categories, retrieved from DAVID database.
About this article
Cite this article
Mehrpour Layeghi, S., Arabpour, M., Shakoori, A. et al. Expression profiles and functional prediction of long non-coding RNAs LINC01133, ZEB1-AS1 and ABHD11-AS1 in the luminal subtype of breast cancer. J Transl Med 19, 364 (2021). https://doi.org/10.1186/s12967-021-03026-7
- Breast cancer
- Luminal A
- Luminal B
- Bioinformatics analyses