Next Article in Journal
In Vitro Setup for Determination of Nanoparticle-Mediated Magnetic Cell and Drug Accumulation in Tumor Spheroids under Flow Conditions
Previous Article in Journal
Comparison of Neurocognitive Functioning and Fine Motor Skills in Pediatric Cancer Survivors and Healthy Children
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 14
Issue 23
10.3390/cancers14235980
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Systematic Investigation of the Diagnostic and Prognostic Impact of LINC01087 in Human Cancers

by
Fatima Domenica Elisa De Palma
1,2,3,4,†,
Vincent Carbonnier
3,4,†,
Francesco Salvatore
2,5,
Guido Kroemer
3,4,6,
Jonathan G. Pol
3,4,*,‡ and
Maria Chiara Maiuri
1,3,4,*,‡
1
Department of Molecular Medicine and Medical Biotechnologies, University of Napoli Federico II, 80131 Napoli, Italy
2
CEINGE-Biotecnologie Avanzate Franco Salvatore, 80145 Napoli, Italy
3
Équipe Labellisée par la Ligue Contre le Cancer, Centre de Recherche des Cordeliers, Inserm U1138, Université Paris Cité, Sorbonne Université, Institut Universitaire de France, 75005 Paris, France
4
Metabolomics and Cell Biology Platforms, Gustave Roussy Cancer Campus, 94805 Villejuif, France
5
Centro Interuniversitario per Malattie Multigeniche e Multifattoriali e Loro Modelli Animali (Federico II, 80131, Napoli, Tor Vergata, Rome and “G. D’Annunzio”, Chieti-Pescara), 80131 Napoli, Italy
6
Department of Biology, Institut du Cancer Paris CARPEM, Hôpital Européen Georges Pompidou, 75004 Paris, France
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this paper as first authors.
‡
These authors contributed equally to this paper as senior authors.
Cancers 2022, 14(23), 5980; https://doi.org/10.3390/cancers14235980
Submission received: 27 September 2022 / Revised: 22 November 2022 / Accepted: 28 November 2022 / Published: 3 December 2022
(This article belongs to the Topic Non-coding RNAs in Cancer: Markers in Disease Progression and Therapy)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

LINC01087 is a recently described long non-coding RNA. Its potential implication within oncology is being increasingly acknowledged. Here, we explored the clinical interest of LINC01087 in the diagnosis and prognostication of multiple cancer types. A series of in-depth in silico analyses revealed LINC01087 as a potential diagnostic indicator of some hormone-related tumors, including breast cancer (BC), ovarian carcinoma, and testicular germ cell tumors, as well as other cancer histotypes, such as esophageal and stomach cancers. Moreover, LINC01087 appears as a prognostic indicator in BC and renal papillary cell carcinoma. Therefore, the quantitation of LINC01087 in tissue specimens might have a clinical utility to improve the management of these pathologies in the future.

Abstract

(1) Background: Long non-coding RNAs may constitute epigenetic biomarkers for the diagnosis, prognosis, and therapeutic response of a variety of tumors. In this context, we aimed at assessing the diagnostic and prognostic value of the recently described long intergenic non-coding RNA 01087 (LINC01087) in human cancers. (2) Methods: We studied the expression of LINC01087 across 30 oncological indications by interrogating public resources. Data extracted from the TCGA and GTEx databases were exploited to plot receiver operating characteristic curves (ROC) and determine the diagnostic performance of LINC01087. Survival data from TCGA and KM-Plotter directories allowed us to graph Kaplan–Meier curves and evaluate the prognostic value of LINC01087. To investigate the function of LINC01087, gene ontology (GO) annotation and Kyoto Encyclopedia of Gene and Genomes (KEGG) enrichment analyses were performed. Furthermore, interactions between LINC01087 and both miRNA and mRNA were studied by means of bioinformatics tools. (3) Results: LINC01087 was significantly deregulated in 7 out of 30 cancers, showing a predominant upregulation. Notably, it was overexpressed in breast (BC), esophageal (ESCA), and ovarian (OV) cancers, as well as lung squamous cell carcinoma (LUSC), stomach adenocarcinoma (STAD), and uterine carcinosarcoma (UCS). By contrast, LINC01087 displayed downregulation in testicular germ cell tumors (TGCT). ROC curve analyses identified LINC01087 as a potential diagnostic indicator in BC, ESCA, OV, STAD, and TGCT. Moreover, high and low expression of LINC01087 predicted a favorable prognosis in BC and papillary cell carcinoma, respectively. In silico analyses indicated that deregulation of LINC01087 in cancer was associated with a modulation of genes related to ion channel, transporter, and peptide receptor activity. (4) Conclusions: the quantification of an altered abundance of LINC01087 in tissue specimens might be clinically useful for the diagnosis and prognosis of some hormone-related tumors, including BC, OV, and TGCT, as well as other cancer types such as ESCA and STAD. Moreover, our study revealed the potential of LINC01087 (and perhaps other lncRNAs) to regulate neuroactive molecules in cancer.

Graphical Abstract

1. Introduction

Carcinogenesis is a dynamic process that inexorably and time-dependently perturbs the normal functions of an organism [1,2]. Besides environmental and lifestyle factors, genetic and epigenetic risk factors contribute dramatically to the transformation of normal cells into malignant cells [1]. Consequently, cancers are highly heterogeneous and complex. In this context, patients can effectively benefit from precision medicine, which may lead to more efficient diagnostic and therapeutic strategies.
Progress in high-throughput technologies contributed to the identification of new molecules, referred to as biomarkers, such as gene sequences, proteins, and metabolites, which are quantifiable in a variety of human samples (e.g., tissue biopsies, blood, feces, sperm, and urine) [3,4,5,6]. This allows biomarkers to precisely discriminate between healthy and disease status, as well as to predict clinical outcomes. Epigenetic modifications, such as aberrant DNA methylation (hypermethylation of tumor suppressor genes or hypomethylation of oncogenes), histone modifications, and dysregulated non-coding RNAs, have been explored as potential biomarkers in the diagnosis and follow-up of patients [7,8,9,10,11]. For instance, the methylation of the O(6)- methylguanine- DNA methyltransferase (MGMT) gene is recognized as a potential prognostic/predictive marker in glioblastoma [12]. Moreover, high DNA methylation of the gene coding for branched chain amino acid transaminase 1 (BCAT1) in plasma may distinguish between lung disease and healthy individuals [13]. Recently, Meyer and colleagues identified DNA methylation patterns that correlate with the response to neoadjuvant chemotherapy in triple-negative breast cancer (TNBC) [14]. Non-coding RNAs, including long (lncRNAs), short (miRNAs), and circular (circRNAs) RNAs, have been recognized as a promising type of early diagnostic and prognostic epigenetic biomarkers in human diseases, taking advantage of their stability to ribonucleases in solid tissues and body fluids [11,15,16,17,18,19,20,21]. LncRNAs consist of transcripts longer than 200 nucleotides, lacking open reading frames [22,23,24]. They can act as tumor suppressors or oncogenes, thereby contributing to cancer pathogenesis and progression when deregulated in tumors [25,26]. Moreover, several studies reported an association between altered lncRNA expression (up or down) and patient survival [25,27]. For instance, a high abundance of the Hox transcript antisense intergenic RNA (HOTAIR) is associated with an unfavorable prognosis in various malignancies, including bladder, gastric, and esophageal squamous cell carcinomas [28,29,30]. Inversely, a decreased expression of the maternally expressed gene 3 (MEG3) lncRNA predicts poor survival in breast cancer (BC) patients [31]. Additionally, in the clinical management of highly heterogeneous malignancies, such as BC or lung cancers, lncRNAs appeared to be highly efficient in discriminating between the different subtypes of organ-specific malignancies [32,33]. Notably, Li and colleagues identified three different molecular subtypes in glioma patients using a lncRNA expression signature [34]. Furthermore, a lncRNA-based molecular classification was described in BC. In particular, among the 4 clusters described, clusters I and II correlated with the basal-like and HER2-enriched subtypes (according to the Prediction Analysis of Microarray 50, PAM50, classification), respectively [35].
The long intergenic non-coding RNA 01087 (LINC01087) is a recently described lncRNA whose function remains poorly investigated. The contribution of LINC01087 to the pathogenesis of cancer, as well as its pivotal role as a biomarker, has been described only in a few oncological indications (i.e., BC, glioma, esophageal squamous cell carcinoma, and thyroid cancer) [32,36,37,38,39]. Recent evidence reported LINC01087 to be overexpressed in glioma patients with respect to healthy subjects [36,37]. Mechanistically, it has been demonstrated that LINC01087 participates in the development of glioma by targeting some microRNAs (miRNAs) (i.e., miR-384 and miR-1277-5p) and, therefore, altering the related signaling pathways [36,37]. For instance, it has been demonstrated in vitro that a high expression of LINC01087 reduces the level of Bcl2 by acting as a miR-384 sponge, thus promoting apoptosis in glioma cells [36]. In addition, its overexpression promotes proliferation, migration, and invasion of thyroid cancer cells through the upregulation of the gene protein phosphatase, Mg2+/Mn2+ dependent 1E (PPM1E) [38]. Furthermore, LINC01087 demonstrated a strong interest in the stratification and prognosis of BC patients [32]. Indeed, the abundance of LINC01087 can discriminate between the luminal and triple-negative (TNBC) subtypes of BC. Furthermore, a high tumor level of this lncRNA is an indicator of prolonged overall survival in BC patients [32].
Here, we shed light on the value of LINC01087 as an indicator for the diagnosis and prognostication of multiple cancer types. Using bioinformatics methods, we first assessed the level of expression of LINC01087 across several tumors. Then, we validated the diagnostic value of LINC01087 using receiver operating characteristics (ROC) curves analysis. We correlated its altered expression with patient survival and clinicopathological features. Finally, the potential functions of LINC01087 were identified by building co-expression networks and gene ontology (GO) enrichment analyses. Our results revealed that LINC01087 is mainly upregulated in multiple cancers, particularly in BC, esophageal, stomach, and ovarian carcinomas, as well as in testicular germ cell tumors. Additionally, LINC01087 could serve as a prognostic indicator in BC and papillary renal cell carcinoma. Interestingly, in silico analyses revealed a potential role for LINC01087 in carcinogenesis through the regulation of gene products that sense or exchange extracellular molecules.

2. Materials and Methods

2.1. Transcriptomics Data Analysis

RNA-sequencing (RNA-seq) and clinical data of 24 primary tumors and paired normal tissues were extracted from The Cancer Genome Atlas (TCGA) pan-cancer 2018 online repository using TCGAbiolinks package in R [23,24,40]. RNA-seq data of 27 types of healthy organs/tissues were downloaded from Genotype-Tissue Expression (GTEx) portal (http://www.gtexportal.org/home/datasets, 1 August 2022). To perform the differential expression analysis of LINC01087, we compared the transcriptomic data of primary tumors extracted from TCGA to paired normal tissues (also retrieved from TCGA) combined with healthy samples of the tissue of origin (extracted from GTEx). The TCGA and GTEx datasets originally comprised 32 and 35 tissues, respectively. However, we excluded several breast cancer and liver subtypes (total of 8 tissues) from both datasets, as well as some organs/tissues with no data (n = 3) from GTEx, resulting in 24 (TCGA) and 27 (GTEx) tissues, reflecting a total of 30 distinct tissues. Of note, only datasets with information about LINC01087 expression were included. Samples with some unassigned values (i.e., “NA”) were removed. Analysis of the differential expression of LINC01087 between normal and malignant tissues was performed using DeSeq2 R/Bioconductor in the 30 types of organs/tissues. LINC01087 was considered as significantly deregulated in cancer when log2 fold change (FC) absolute value ≥ 1.5 and p-value ≤ 0.05 (Wilcoxon test). For each dataset, the differential expression of LINC01087 was calculated based on the dataset median value used as a threshold. The association between the expression of LINC01087 in different tumor types and the clinicopathological parameters (i.e., age, gender, tumor stage, tumor size, tumor status, lymph node, and metastasis involvement) was calculated using Fisher’s exact test. For each dataset, a high versus low discrimination of LINC01087 expression was achieved, using the dataset median value as a threshold.

2.2. Survival Curve Analysis

To evaluate the effect of LINC01087 on the overall survival (OS) across patients affected with various cancer histotypes, we used the pan-cancer RNA-seq data from TCGA (https://portal.gdc.cancer.gov 15 September 2022) and the Kaplan–Meier Plotter (KM-Plotter) datasets (https://kmplot.com, 2 July 2022) [41,42]. In KM-Plotter, survival curves were generated by selecting the three following filters: median, JetSet best probe set, and all datasets. Survival R package v.3.4 and Survminer R package v.0.4.9 were applied to TCGA datasets in order to calculate the hazard ratio (HR), 95% confidence intervals (CI), and p-values. For each TCGA dataset, depending on the number of available samples, patients were subdivided into two groups according to the median expression level of LINC01087: high or low. A p-value ≤ 0.05 was considered statistically significant.

2.3. ROC Curve Analysis

ROC curves were generated and the area under the curve (AUC) was calculated using the R packages MASS and pROC [43,44]. All tests were performed in R statistical software (ver. 4). AUC ≥ 0.7 and p-value ≤ 0.05 were applied as a cut-off criterion to consider LINC01087 as a good diagnostic biomarker. For each dataset, a high versus low discrimination of LINC01087 expression was considered based on the dataset best cut-off value.

2.4. Network and Pathway Enrichment Analyses of LINC01087

For the identification of interactions between LINC01087 and mRNAs and miRNAs, we employed a dual strategy. In the first step, we performed a Spearman correlation test to identify mRNAs and miRNAs that correlated with LINC01087 expression in BC, ESCA, OV, STAD, and TGCT by extracting data from TCGA datasets. We defined significant associations as having a correlation coefficient (R) > 0.4 and a p-value < 0.05. Volcano plots and Venn diagrams were drawn using EnhancedVolcano and VennDiagram R packages, respectively. In the second step, we used the web tool DIANA-LncBASE v3 (https://diana.e-ce.uth.gr/lncbasev3, 14 October 2022) to identify candidate miRNA targets of LINC01087 [45]. Then, we took advantage of the miRTarBase (https://mirtarbase.cuhk.edu.cn/~miRTarBase/miRTarBase_2022/php/index.php, 20 October 2022) to predict mRNA targets of the miRNAs derived from both analyses [46]. Finally, the target mRNA candidates were compared with the mRNAs that correlated with LINC01087 expression in the initial analysis.
To investigate the biological pathways and processes, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) term enrichment analyses were performed using the clusterProfiler R package on genes significantly associated with LINC01087 expression (R > 0.1 p < 0.05) in BC, ESCA, STAD, and TGCT tumor tissues derived from TCGA (p-value were adjusted using the Benjamini–Hochberg procedure).

3. Results

3.1. Diagnostic Value of LINC01087

To assess the diagnostic potential of LINC01087, we measured its abundance in 30 types of tissue/organs extracted from the TCGA and GTEx databases (Table 1). In particular, we compared the transcriptomics data of primary tumors from TCGA to paired normal tissues (retrieved from TCGA) combined with healthy samples of the tissue of origin (extracted from GTEx).
LINC01087 demonstrated significant upregulation (|log2 FC| ≥ 1.5, p ≤ 0.05) in 6 out of the 30 cancers as compared to their normal (i.e., non-tumor or healthy) tissue counterparts (Figure 1 and Table 1).
Notably, LINC01087 displayed elevated levels in BC (log2FC = 2.2, p = 2.32e-35), esophageal carcinoma (ESCA, log2FC = 3.7, p = 4.08e-35), lung squamous cell carcinoma (LUSC, log2FC = 2.73, p = 0.003), ovarian cancer (OV, log2FC = 1.5, p = 1.11e-14), stomach adenocarcinoma (STAD, log2FC = 4, p = 2.49e-50), and in uterine carcinosarcoma (UCS, log2FC = 3.2, p = 0.0001) (Figure 1 and Table 1). In addition, a smaller (log2FC < 1.5) but significantly (p ≤ 0.05) altered expression of LINC01087 was detected in clear cell (ccRCC, log2FC = 0.8, p = 1.25e-15) and papillary (pRCC, log2FC = 0.9, p = 0.005) renal cell carcinomas, hepatocellular carcinoma (HCC, log2FC = 0.8, p = 4.28e-07), lung adenocarcinoma (LUAD, log2FC = 1.3, p = 0.002), rectum adenocarcinoma (READ, log2FC = 0.7, p = 0.001), and skin cutaneous melanoma (SKCM, log2FC = 0.5, p = 7.99e-23) (Table 1).
By contrast, LINC01087 demonstrated a strong downregulation in testicular germ cell tumors (TGCT, log2FC = −6.1, p = 6.59e-73) (Figure 1 and Table 1). Moreover, it was less (log2FC ≥ −1.5), but significantly (p ≤ 0.05) expressed in acute myeloid leukemia (LAML, log2FC = −1, p = 9.46e-08), brain lower grade glioma (LGG, log2FC = −0.3, p = 2.43e-04), pancreatic adenocarcinoma (PAAD, log2FC = −1.2, p = 6.86e-20), pheochromocytoma and paraganglioma (PGPG, log2FC = −0.7, p = 3.87e-07), and thyroid carcinoma (THCA, log2FC = −0.6, p = 3.13e-10) (Table 1).
Next, we measured the diagnostic performance of LINC01087 by plotting and interpreting ROC curves (Table 1). Values of area under the curve (AUC) superior to 0.7 designated LINC01087 as a good diagnostic biomarker in BC (AUC = 0.7, p = 1.16e-35), as well as esophageal carcinoma (AUC = 0.8, p = 2.04e-35), ovarian (AUC = 0.7, p = 5.56e-15) and stomach (AUC = 0.8, p = 1.25e-50) cancers, and in testicular germ cell tumors (AUC = 1, p = 3.30e-73) (Table 1). Furthermore, we observed lower (AUC < 0.7), but significant (p< 0.05), AUC values of LINC01087 in clear cell (ccRCC: AUC = 0.7, p = 6.24e-16) and papillary renal cell carcinomas (pRCC: AUC = 0.6, p = 0.002), acute myeloid leukemia (AML: AUC = 0.6, p = 4.73e-08), lower grade glioma (LGG: AUC = 0.5, p = 0.0001), hepatocellular carcinoma (HCC: AUC = 0.6, p = 2.14e-07), lung adenocarcinoma (LUAD: AUC = 0.5, p = 0.001), lung squamous cell carcinoma (LUSC: AUC = 0.6, p = 0.001), pheochromocytoma and paraganglioma (PCPG: AUC = 0.6, p = 1.93e-07), skin cutaneous melanoma (SKCM: AUC = 0.6, p = 3.99e-23), and uterine carcinosarcoma (UCS: AUC = 0.6, p = 6.69e-05) (Table 1).
Finally, we investigated whether the deregulation of LINC01087 could be associated with clinicopathological features (i.e., age, gender, tumor stage, tumor size, tumor status, lymph node, and metastasis involvement). To address this point, we focused on cancers in which LINC01087 showed both up- or down-regulation (|log2FC|> 1.5, p < 0.05) and diagnostic performance (AUC > 0.7, p < 0.05), when such information was available. These criteria restricted the analysis to BC, ESCA, OV, STAD, and TGCT (Table 2, Table S1). LINC01087 expression appeared to be significantly associated with the tumor status (p = 0.04) and tumor size (p = 0.01) in BC (Table 2). In TGCT, LINC01087 significantly correlated with the size of tumor (p = 0.05) (Table 2).
Altogether, these results imply that the perturbation of LINC01087 expression might play a critical role in carcinogenesis. Moreover, LINC01087 shows great potential as an indicator for the diagnosis of BC, as well as esophageal carcinoma and stomach adenocarcinomas, OV and TGCT. These findings require further experimentation to validate such interest.

3.2. Prognostic Value of LINC01087 in Cancer

Subsequently, we examined the prognostic value of LINC01087 by evaluating the relationship between the expression of LINC01087 and patient overall survival (OS). For this purpose, cohorts gathered in the TCGA directory were used for derivation, while inputs of the KM-Plotter database initially served as validation cohorts.
Low expression of LINC01087 denoted a favorable prognosis in pRCC (n = 412, p = 2.86E-04) (Table 1). Contrarily, the presence of LINC01087 correlated with a better OS in patients affected with BC (n = 1685, p = 0.01). Such prognostic interest was validated using KM-Plotter datasets (Table 1).
Further interrogation of the KM-Plotter database revealed that a low abundance of LINC01087 matched an extended OS in HCC, as well as in uterine corpus endometrial carcinoma (Table 1).

3.3. Impact of LINC01087 on Cancer-Related Pathways

To capture significant insights into LINC01087 functions in cancers, we performed a correlation-based expression analysis extracting data from TCGA datasets. In particular, we looked for genes that significantly correlated (Spearman correlation test: R ≥ |0.4|, p ≤ 0.05) with LINC01087 expression in each of the tissue types where LINC01087 demonstrated a strong diagnostic potential (log2 FC ≥ |1.5|, p ≤ 0.05, and AUC ≥ 0.7, p ≤ 0.05), namely BC, ESCA, OV, STAD, and TGCT (Analysis 1—Table 1, Figure S1, and Tables S2–S6). We found significant positive (R ≥ 0.4, p ≤ 0.05) correlations between the expression of LINC01087 and 574 transcripts in BC, 691 in ESCA, 141 in OV, 773 in STAD, and 636 in TGCT (Tables S2–S6). LINC01087 was associated with pseudogenes, microRNAs (miRNAs), and other lncRNAs, as well as with protein-coding genes. Out of the 10 transcripts that were the most significantly correlated with LINC01087 in BC, 6 were pseudogenes, 2 were lncRNAs, and 2 were protein-coding genes belonging to the POTE family (POTEH, POTEI). Similarly, in ESCA, among the top 10 genes associated with LINC01087 were pseudogenes and lncRNAs, as well as the neurofilament heavy polypeptide (NEFH) gene, whose deregulation had been associated with the development of several malignancies, such as BC [47]. Several genes were negatively and modestly (R ≤ −0.2, p ≤ 0.05) correlated with LINC01087; therefore, they were not considered for further analysis. To refine our analysis, we compared the correlating gene sets of every cancer type (R ≥ 0.4, p ≤ 0.05) for overlapping genes. We identified 2 transcripts that were positively associated with LINC01087 expression in all 5 cancers: the pseudogene POTE ankyrin domain family member K (POTEKP) and the uncharacterized lncRNA AC093838.1 (Figure 2A).
POTEKP showed the strongest positive correlation with LINC01087 in BC (R = 0.96), ESCA (R = 0.95), STAD (R = 0.93), and TGCT (R = 0.98), with the exception of OV, where it showed a moderate yet significant correlation (R = 0.61) (Tables S2–S6). In OV, AC093838.1 demonstrated the strongest association with LINC01087 (R = 0.91). The latter transcript also displayed a strong correlation with LINC01087 in the 4 other types of cancer (Tables S2–S6). In addition, the transcript of mediator complex subunit 15 pseudogene 4 (MED15P4) showed a strong association with LINC01087 in 3 out of the 5 tumor types (i.e., OV, STAD and TGCT), and a moderate correlation in ESCA (Tables S2–S6).
In the next step, we aimed at understanding the biological impact of LINC01087 in cancer. To this aim, we performed GO annotation and KEGG enrichment analyses on the set of genes correlated with LINC01087 expression in the 5 tumor types (n = 2458, |R| ≥ 0.1, p ≤ 0.05). We identified a total of 35 GO terms subcategorized as follows: 1 “biological process” (BP), 13 “cellular component” (CC) and 21 “molecular function” (MF) (Figure 2B, Table S7). Among them, the 3 largest GO terms referred to “channel activity”, “passive transmembrane transporter activity”, and “synaptic membrane”. Notably, these clusters encompassed genes encoding gamma-aminobutyric acid [GABA] receptors (i.e., GABRG1, GABRG2, GABBR2, and GABRA2,) and glutamate ionotropic receptors (i.e., GRIN3A, GRIA2, GRIA4, and GRIN1) (Figure 2B, Table S7). Of note, this annotation analysis revealed an over-representation of genes uncovered in ESCA, such as GABRA1, GRIN3A GRIA2, GRIA4, and GLRB (Figure 2B, Supplementary Table S7). GABRG2 was linked with both ESCA and STAD. AQP11, KCNH3, GABRA6, and CHRNA10 were exclusively associated with TGCT, while GJA8 and CLCN1 were correlating with LINC01087 only in OV (Figure 2B, Table S7). Concordant results were obtained using the KEGG annotation analysis (Figure 2C, Table S8).
In summary, this orthogonal investigation revealed that LINC01087 might impact cancer progression by modulating the ability of cells to interact with their environment.

3.4. LINC01087 Functions as Competing Endogenous RNA

LncRNAs participate in carcinogenesis by modulating gene expression via either direct binding to mRNAs, or indirectly by functioning as a competing endogenous RNA (ceRNA) in the cytoplasm. This latter mechanism is characterized by the binding of lncRNAs to miRNAs, thus preventing their interaction with target mRNAs [48]. LINC01087 displays a cytoplasmic distribution as predicted by lncLocator database (score = 0.76, Table S9). Therefore, we investigated the potential of LINC01087 to act as a miRNA sponge.
To construct a potential LINC01087-miRNA-mRNA network, we used a dual strategy. In a first step based on the TCGA datasets, we screened for miRNAs the expression of which correlated with that of LINC01087 in the 5 types of cancer (R ≥ |0.4|, p ≤ 0.05) (Analysis 2—Figure 3). Sixty-four miRNAs were positively correlated with the expression of LINC01087 (Figure 3, Table S10). The top 5 transcripts demonstrated a moderate correlation (0.64 ≤ R ≤ 0.72) and comprised miR-5692B, miR-4666A, miR-31199-1, miR-1252, and miR-7157 (Table S10). Then, the miRTarBase was interrogated, and 2396 potential mRNA targets of these 64 miRNAs were identified (Figure 3). Finally, we compared this list of candidates (n = 2396, Analysis 2) with the initial clusters of genes positively correlated with LINC01087 in the 5 tumor types from TCGA (n = 2458, Analysis 1). A total of 31 genes overlapped between the two datasets (Figure 3, Table S11).
In the second step, we used DIANA-LncBASE v3 to discover miRNAs interacting with LINC01087. A total of 20 miRNAs were identified as potential targets of LINC01087 (Analysis 3—Figure 3, Table S12). Most of them had already been involved in carcinogenesis, such as miR-19a-3p and miR-181-5p [49,50]. Furthermore, 4 of the 20 miRNAs belonged to the let-7 miRNA family (i.e., hsa-let-7a-5p, hsa-let-7d-5p, hsa-let-7e-5p, hsa-miR-98-5p). However, no overlapping miRNAs were observed between the DIANA- and TCGA-based analyses. Then, we executed the miRTarBase to detect mRNA targets of these 20 miRNAs. A total of 5532 mRNA targets were identified in Analysis 3. Sixty-eight of them were common with the list of 2458 genes identified in Analysis 1 (Figure 3, Table S13).
Next, we compared the 31 mRNA targets annotated in Analysis 2 with the 68 mRNA targets uncovered in Analysis 3 and identified 10 overlapping messengers representing putative targets of LINC01087: chromosome 8 open reading frame 37 (C8orf37), ephrin A5 (EFNA5), hook microtubule-tethering protein 3 (HOOK3), myozenin 3 (MYOZ3), pleomorphic adenoma gene 1 (PLAG1), DNA Polymerase Iota (POLI), solute carrier family 2member 3 (SLC2A3), small nuclear ribonucleoprotein polypeptide A’ (SNRPA1), as well as two isoleucine/glutamine (IQ) motif genes, namely Purkinje cell protein 4 like 1 (PCP4L1) and IQ motif containing G (IQCG) (Figure 3, Tables S11 and S13).
Ultimately, this study led to the construction of a core network of 10 mRNA targets and 28 miRNAs related to LINC01087 (Table 3).
In more detail, LINC01087 seemed to interfere with 2 miRNAs targeting the messenger of C8orf37 (namely has-miR-644a and has-miR-7-5p), EFNA5 (hsa-miR-1305, hsa-miR-92a-3p), HOOK3 (hsa-miR-4505, hsa-miR-181a-5p), MYOZ3 (hsa-miR-7975, hsa-miR-197-3p), and PCP4L1 (hsa-miR-4505, hsa-miR-423-5p), with 3 miRNAs targeting the mRNAs of SLC2A3 (hsa-miR-606, hsa-miR-34a-5p, hsa-miR-98-5p) and SNRPA1 (hsa-miR-3174, hsa-miR-34a-5p, hsa-miR-98-5p), and 4 miRNAs interfering with the transcripts of IQCG (hsa-miR-3671, hsa-miR-181a-5p, hsa-miR-181b-5p, hsa-miR-181d-5p), PLAG1 (hsa-miR-569, hsa-miR-181a-5p, hsa-miR-181b-5p, hsa-miR-98-5p), and POLI (hsa-miR-4255, hsa-miR-19a-3p, hsa-miR-19b-3p, hsa-miR-92a-3p). Interestingly, most of these mRNAs and miRNAs have been connected to tumor suppressive and/or oncogenic functions in multiple cancers [49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67].
Altogether, these results documented a constellation of both protein-coding and translation-interfering transcripts that are directly or indirectly regulated by LINC01087 and that might contribute to the initiation and progression of neoplasms.

4. Discussion

LncRNAs represent promising epigenetic diagnostic, prognostic, and therapeutic biomarkers [25]. No less than 23 trials (source: clinicaltrials.gov (accessed on 1 August 2022)) aim at evaluating their clinical relevance in a variety of malignancies, including bladder (NCT05270174), colorectal (NCT04729855), gastric (NCT05397548, NCT05334849), head and neck (NCT04946968), kidney (NCT04946266), liver (NCT05088811), and prostate (NCT05141383) cancers. In the present report, bioinformatics analyses uncovered the role of LINC01087 in the diagnosis and prognostication of different human cancers.
To examine the clinical interest of LINC01087, we first analyzed its expression across a large array of tumor types using TCGA and GTEx databases, then calculated the ROC scores to assess its diagnostic performance. Among the 30 cancers evaluated, LINC01087 displayed a significant modulation of its expression in 7 of them, mainly hormone-dependent. More precisely, LINC01087 showed an increased expression in BC, OC, and UCS, and an inversely reduced level in TGCT as compared to non-tumor/healthy control samples. Additionally, an elevation of LINC01087 expression was witnessed in LUSC, as well as in two digestive cancers, namely ESCA and STAD.
Our results were consistent with the literature regarding the abundance of LINC01087 in BC and ESCA [36,37,38,68]. Nevertheless, previous research reported a higher expression of LINC01087 in THCA patients than in normal individuals [38], which was not validated in this study. This discrepancy could originate from differences in the size of the cohorts (n = 1224 in the present work versus n = 30 in the former article by Yin et al.) as well as the detection method applied (RNA-seq versus RT-qPCR) [38].
Out of the 7 cancers harboring an aberrant expression of LINC01087, a diagnostic test revealed promising performances in 5 of them (i.e., ROC curve AUC > 0.7). In detail, the measurement of LINC01087 levels was an excellent (AUC = 1) indicator of malignancy in TGCT, and a fair (0.7 < AUC ≤ 0.8) diagnostic indicator in BC, ESCA, OV, and STAD. Of note, we reported that the measurement of the level of LINC01087 in BC also allowed us to further distinguish between the luminal and triple-negative subtypes [32].
Next, we evaluated the prognostic value of LINC01087 by comparing cancer patient OS according to its level of expression using the cohorts of TCGA and KM-plotter databases for derivation and validation, respectively. LINC01087 deregulation predicted patient outcome in two oncological indications. First, a downregulation of LINC01087 correlated with favorable survival in patients suffering from pRCC. Second, a high expression level of LINC01087 indicated a better prognosis in subjects affected with BC.
We proceeded by evaluating the potential functional significance of LINC01087 in the 5 cancer types in which it was most upregulated, namely, BC, ESCA, OV, STAD, and TGCT. At first sight, we screened for genes correlated with LINC01087 expression in the 5 tumor types. Such analysis identified 2 transcripts overlapping in the 5 cancer types, namely POTEKP and AC093838.1. These genes are both located on chromosome 2 but their role in human malignancies has not been elucidated.
Interestingly, GO annotation and KEGG enrichment analyses performed on the genes positively correlated with LINC01087 expression in BC, ESCA, OV, STAD, and TGCT, found a spectrum of genes and pathways associated with the regulation of ion channel activities and neurotransmission. This cluster consisted mainly of genes related to the receptors for GABA (i.e., GABRG1, GABRG2, GABBR2, GABRA2, GABRA6) and α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) (i.e., GRIN3A, GRIA2, GRIA4, and GRIN1). Based on this information, we speculate that the involvement of LINC010187 in tumorigenesis might be associated with the “neuro-cancer axis” [69]. Of note, neuronal cells and immune and cancer cells can synthesize and release neurotransmitters such as GABA [70,71,72,73]. GABA exerts its multiple functions through ionotropic GABA A receptors and metabolotropic GABA B receptors. The aberrant expression of GABA and GABA receptors has been associated with cancer development, metastasis, and suppression of immunosurveillance [71,74]. In this scenario, lncRNAs, as well as miRNAs, may regulate neuronal activity and target GABA receptors [59,75,76]. For instance, one recent study examined the role of the lncRNA Gm37494 in the context of osteoarthritis, focusing on its capability to target GABRA1 by binding to miR-181a-5p [76]. Therefore, it appears plausible, yet remains to be demonstrated, that LINC01087 (and probably other lncRNAs) regulate GABA receptor activity during the formation and progression of specific cancer types.
LncRNAs can regulate gene expression by interacting with DNA, RNA, and proteins [23]. Notably, they can either directly alter the transcription of neighboring and distant protein-coding genes, or indirectly, by targeting and thereby modulating the expression of miRNAs acting as ceRNAs. Therefore, we constructed a LINC01087-miRNA-mRNA network by interrogating the TCGA and DIANA-LncBASE datasets. Next, miRTarBase was used to discover possible interactions between miRNAs and mRNA targets in both datasets. Finally, we overlapped the mRNA targets with the cluster of genes correlated with LINC01087 expression. Our analysis yielded a list of 28 miRNAs and 10 mRNAs that constituted a core panel of putative targets of LINC01087.
Some mRNAs and miRNAs of the established LINC01087 network are related to malignant disease [51,52,53]. For instance, miR-3174 (targeting SNRPA1 mRNA) exhibits oncogenic activity, whereas miR-34a-5p and hsa-miR-644a (targeting SLC2A3 and C8orf37 mRNAs, respectively) possess tumor-suppressive properties [60,61,77,78,79]. Nevertheless, most of the miRNAs related to LINC01087 (i.e., miR-7-5p, miR-19a-3p, miR-19-b-3p, miR-92b-3p, hsa-miR-98-5p, miR-181a-5p, miR-181d-5p, miR-197-3p, miR-423-5p, miR-569, and hsa-miR-1305) have an ambiguous impact on carcinogenesis, acting as either tumor suppressor genes or oncogenes depending on the tumor type [58,63,65,66,67,80,81,82,83,84,85,86,87,88]. Among them, miR-92b-3p inhibited the proliferation and invasion of pancreatic cancer cells through the downregulation of GABRA3 [59]. In an analogous way, we found that LINC01087 expression was correlated with several GABRA genes, such as GABRA2 and GABRA6. Furthermore, some miRNAs of the LINC01087 network like miR-34a-5p, miR-197-3p, miR-423-5p, and hsa-miR-3174, were reported as targets of lncRNAs. For instance, miR-34a-5p can be trapped by LINC00665, HOTAIR, or else nuclear paraspeckle assembly transcript 1 (NEAT1) [62,64,89,90,91,92,93]. NEAT1 is also able to target miR-98-5p to promote glioma cancer progression [94]. Out of the 28 miRNAs, a few have been scarcely explored in cancer (e.g., miR-606, miR-3671, miR-4255, and miR-7975) [95,96,97]. For instance, miR-606 was mentioned in a single study to be upregulated in pancreatic cancer [98]. Similarly, one study reported a high expression level of miR-4255 in exosomes collected from patients with prostate cancer [99]. Moreover, Lin and collaborators demonstrated that the glycoprotein lipocalin 2 prevents oral cancer metastasis through miR-4505-mediated suppression of carbonic anhydrase IX expression [100].
Concerning the protein-encoding genes identified in the LINC01087 network, several of them were described as potential cancer-relevant biomarkers. Indeed, PLAG1 was reported as an immunohistochemical marker in pleomorphic adenoma [52]. Moreover, PLAG1 was found to be upregulated in OV, in which its silencing promoted chemosensitivity [50]. Another example concerns the role of HOOK3 as a predictor of poor prognosis in prostate cancer [51]. Furthermore, in accordance with our data, He and colleagues predicted HOOK3 as a target of miR-181-5p [101]. Similarly, upregulation of POLI predicted poor prognosis in esophageal squamous cell carcinoma (ESCC) patients [68,69,70]. Furthermore, the single amino acid polymorphism Thr706Ala in POLI was indicative of a higher risk of developing lung adenocarcinoma and squamous cell carcinoma in individuals of <61 years [102]. The aforementioned investigations support a possible role of the LINC01087-POLI interaction in cancer, in particular ESCA.
Here, we found a high expression of SLC2A3 in ESCA. SLC2A3, also known as glucose transporter 3 (GLUT3), is a member of the SLC2 family, mediating glucose transport across the plasma membrane. The oncogenic potential of SLC2A3 has been extensively investigated. An overexpression or copy number gain of SLC2A3 has been witnessed in liquid and solid tumors [103,104,105,106,107,108]. For instance, high levels of SLC2A3 were associated with unfavorable outcomes in CRC, gastric (GC), and head and neck squamous cell carcinoma (HNSCC) [105,106,107]. Additionally, SLC2A3 seems directly targeted by miRNAs along tumorigenesis [109,110,111]. For instance, miR-29c inhibited cell growth and glucose metabolism in prostate cancer by targeting SLC2A3 [110]. Moreover, a decreased expression of miR-106a increased the expression of SLC2A3, thereby contributing to progression and poor prognosis in glioblastoma patients [109]. In GC, cell glucose metabolism and growth are regulated by the miRNA-129-5p/SCL2A3 axis via the phosphatidylinositol 3-kinase/protein kinase B (PI3K-Akt) and mitogen-activated protein kinase (MAPK) signaling pathways [111]. Accordingly, we observed a relationship between LINC01087 and certain SLC membrane transporters (e.g., SLC28A1, SLC9A3, SLC4A9, and SLC17A8) in luminal and TNBC BC subtypes [32]. These similarities support an intimate link between the solute carrier family and LINC01087 in cancer.
EFNA5, a member of the ephrinA ligand family, is deregulated (up or down) in a variety of cancers [112,113,114,115,116]. A recent study explored the mRNA and protein expression, as well as DNA methylation, immune-related signature, and prognostic value of the EFNA family using a multi-omics integrative approach [56]. EFNA5 mRNA was found to be overexpressed in most tumors compared to normal tissues. Moreover, an association between the mRNA level of EFNA5 and DNA methylation of its promoter were documented for ccRCC, chRCC, pRCC, mesothelioma, prostate cancer, thymoma, and diffuse large B-cell lymphoma. At the protein level, EFNA5 was downregulated in glioblastoma multiforme and ccRCC, but highly expressed in endometrial carcinoma. In addition, a high abundance of EFNA5 in both TCGA and gene expression omnibus (GEO) databases correlated with poor OS in patients with advanced LUAD [56].
The gene SNRPA1, a spliceosome component responsible for processing pre-mRNA into mRNA, showed strong prognostic potential in prostate cancer, HCC, and ccRCC [55,117,118]. Its oncogenic role in CRC has also been demonstrated [119]. Here, we observed high expression of SNRPA1 in STAD, as well as a positive correlation between LINC01087 and SNRPA1. Interestingly, Chen and colleagues reported the interaction between the lncRNA RMRP (RNA component of mitochondrial RNA-processing endoribonuclease) and the protein SNRPA1 as a mediator of TP53 inactivation and as a promoter of cell proliferation in CRC [120].
IQCG is a protein-encoding gene involved in sperm axoneme assembly. Here, we reported an association between LINC01087, miR-181a-5p, and IQCG mRNAs in cancer. The link between miR-181a-5p and IQCG has previously been observed in endometrial carcinoma using an in silico approach [101].
Conversely, there is limited knowledge on C8orf37 and PCP4L1 in malignancies. While C8orf37 has been only associated with genetic disorders, the possible involvement of PCP4L1 in CRC and prostate cancer progression requires further confirmation [121,122,123,124].
Based on our results and previous observations, we speculate that the deregulation of LINC01087 might be related to the carcinogenesis of BC, ESCA, OV, STAD, and TGCT. However, the mechanisms behind LINC01087 involvement in the development/progression of these malignant diseases remain to be explored experimentally.
Since its discovery, interest in and functional knowledge of LINC01087 in the field of oncology have rapidly expanded. Further bench-to-bedside investigations are required to demonstrate the clinical relevance of LINC01087 as a diagnostic and prognostic biomarker.

5. Conclusions

In summary, the deregulation of LINC01087 plays a potential role in the diagnosis and prognosis of different tumor types. Notably, a high expression of LINC01087 exhibits a strong diagnostic value in breast, esophageal, stomach, and ovarian carcinomas. Moreover, a decreased expression of LINC01087 could be used for the diagnosis of malignant testicular germ tumors. LINC01087 also appears to predict the survival of patients suffering from breast and papillary kidney malignancies. Furthermore, the involvement of LINC01087 in cancer development might be related to its direct and indirect (miRNA-mediated) role in modulating biochemical exchanges between the intracellular milieu and the surrounding microenvironment. Future studies must shed light on the biological activity of LINC01087 and confirm its clinical utility as a diagnostic and/or prognostic biomarker. Ultimately, the quantitation of LINC01087 may ease clinical decision-making and improve the management of these oncological pathologies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers14235980/s1, Figure S1: Differentially expressed genes correlated with LINC01087 expression in several tissue types extracted from TCGA. Table S1: Non-significant associations between LINC01087 expression and clinicopathological characteristics in TCGA and GTEx data. Table S2: Genes significantly correlated with LINC01087 expression in breast cancer (TGCA); Table S3: Genes significantly correlated with LINC01087 expression in esophageal carcinoma (TGCA); Table S4: Genes significantly correlated with LINC01087 expression in ovarian cancer (TGCA); Table S5: Genes significantly correlated with LINC01087 expression in stomach adenocarcinoma (TGCA); Table S6: Genes significantly correlated with LINC01087 expression in testicular germ cell tumor carcinomas (TGCA); Table S7: GO annotation of genes positively associated with LINC01087 expression in BC, ESCA, OV, STAD and TGCT; Table S8: KEGG pathways associated with LINC01087 expression in BC, ESCA, OV, STAD and TGCT; Table S9: Prediction of subcellular localization of LINC0187 using “lncLocator” database; Table S10: List of the 64 miRNAs significantly correlated with LINC01087 expression in BC, ESCA, OV, STAD and TGCT in TCGA datasets (Analysis 2); Table S11: List of the 31 mRNA targets of LINC01087-related miRNAs that overlapped between Analysis 1 and Analysis 2; Table S12: List of the 20 miRNAs identified to interact with LINC01087 in the DIANA-LncBASE repository (Analysis 3); Table S13: List of the 68 mRNA targets of LINC01087-related miRNAs that overlapped between Analysis 1 and Analysis 3.

Author Contributions

Conceptualization: F.D.E.D.P. and M.C.M.; Data Curation: V.C.; Formal Analysis: F.D.E.D.P. and V.C.; Supervision: J.G.P. and M.C.M.; Validation: F.D.E.D.P. and J.G.P.; Writing-original draft: F.D.E.D.P., J.G.P. and G.K.; Writing-review and editing: F.D.E.D.P., V.C., J.G.P., F.S., G.K. and M.C.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets used and/or analysed during the current study were published previously as indicated in “Methods” Section. In particular, RNA-seq data were downloaded from the TCGA pancancer 2018 using TCGAbiolinks package in R and from Genotype-Tissue Expression (GTEx) portal (http://www.gtexportal.org/home/datasets, 1 August 2022). The TCGA clinical data were downloaded from TCGA pancancer 2018 using TCGAbiolinks package in R. Survival curves were generated extracting data from TCGA (https://portal.gdc.cancer.gov, 15 September 2022) and were plotted using the online survival analysis tool “Kaplan-Meier Plotter” (http://kmplot.com/analysis/, 2 July 2022) [41,42]. miRNA targets of LINC01087 were identified using both TCGA datasets and DIANA-LncBASE v3 tool (https://diana.e-ce.uth.gr/lncbasev3, 14 October 2022) [45]. mRNA-miRNA interactions were established using miRTarBase (https://mirtarbase.cuhk.edu.cn/~miRTarBase/miRTarBase_2019/php/index.php 20 October 2022) [46].

Acknowledgments

J.G.P. is supported by the SIRIC Cancer Research and Personalized Medicine (CARPEM); Multi-Organism Institute (ITMO) Aviesan Cancer (National Alliance for Life Sciences and Health), Institut National du Cancer (INCa), and Fondation pour la Recherche Médicale (FRM). M.C.M. is supported by CARPEM. GK is supported by the Ligue contre le Cancer (équipe labellisée); Agence National de la Recherche (ANR)—Projets blancs; AMMICa US23/CNRS UMS3655; Association pour la recherche sur le cancer (ARC); Cancéropôle Ile-de-France; Fondation pour la Recherche Médicale (FRM); a donation by Elior; Equipex Onco-Pheno-Screen; European Joint Programme on Rare Diseases (EJPRD); European Research Council (ICD-Cancer), European Union Horizon 2020 Projects Oncobiome and Crimson; Fondation Carrefour; Institut National du Cancer (INCa); Institut Universitaire de France; LabEx Immuno-Oncology (ANR-18-IDEX-0001); a Cancer Research ASPIRE Award from the Mark Foundation; the RHU Immunolife; Seerave Foundation; SIRIC Stratified Oncology Cell DNA Repair and Tumor Immune Elimination (SOCRATE); and SIRIC Cancer Research and Personalized Medicine (CARPEM). This study contributes to the IdEx Université de Paris ANR-18-IDEX-0001. The graphical abstract and Figure 3 were created with “BioRender.com”.

Conflicts of Interest

F.D.E.D.P., F.S., G.K., J.G.P. and M.C.M. are listed as co- inventors on an international patent (“Methods for diagnosing, prognosing and managing treatment of breast cancer”; WO 2021/170777 A1). V.C. declares no conflict of interest.

References

  1. Hanahan, D.; Weinberg, R.A. The Hallmarks of Cancer. Cell 2000, 100, 57–70. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Hanahan, D. Hallmarks of Cancer: New Dimensions. Cancer Discov. 2022, 12, 31–46. [Google Scholar] [CrossRef] [PubMed]
  3. Henry, N.L.; Hayes, D.F. Cancer Biomarkers. Mol. Oncol. 2012, 6, 140–146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Precone, V.; Del Monaco, V.; Esposito, M.V.; De Palma, F.D.E.; Ruocco, A.; Salvatore, F.; D’Argenio, V. Cracking the Code of Human Diseases Using Next-Generation Sequencing: Applications, Challenges, and Perspectives. Biomed Res. Int. 2015, 2015, 161648. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Olivier, M.; Asmis, R.; Hawkins, G.A.; Howard, T.D.; Cox, L.A. The Need for Multi-Omics Biomarker Signatures in Precision Medicine. Int. J. Mol. Sci. 2019, 20, 4781. [Google Scholar] [CrossRef] [Green Version]
  6. Srivastava, A.; Creek, D.J. Discovery and Validation of Clinical Biomarkers of Cancer: A Review Combining Metabolomics and Proteomics. Proteomics 2019, 19, e1700448. [Google Scholar] [CrossRef]
  7. Arechederra, M.; Recalde, M.; Gárate-Rascón, M.; Fernández-Barrena, M.G.; Ávila, M.A.; Berasain, C. Epigenetic Biomarkers for the Diagnosis and Treatment of Liver Disease. Cancers 2021, 13, 1265. [Google Scholar] [CrossRef]
  8. Parrella, P. The Value of Epigenetic Biomarkers in Breast Cancer. Biomark. Med. 2018, 12, 937–940. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Khan, S.A.; Reddy, D.; Gupta, S. Global Histone Post-Translational Modifications and Cancer: Biomarkers for Diagnosis, Prognosis and Treatment? World J. Biol. Chem. 2015, 6, 333–345. [Google Scholar] [CrossRef]
  10. Lu, H.; Liu, Y.; Wang, J.; Fu, S.; Wang, L.; Huang, C.; Li, J.; Xie, L.; Wang, D.; Li, D.; et al. Detection of Ovarian Cancer Using Plasma Cell-Free DNA Methylomes. Clin. Epigenetics 2022, 14, 74. [Google Scholar] [CrossRef]
  11. Anastasiadou, E.; Jacob, L.S.; Slack, F.J. Non-Coding RNA Networks in Cancer. Nat. Rev. Cancer 2018, 18, 5–18. [Google Scholar] [CrossRef] [PubMed]
  12. von Deimling, A.; Korshunov, A.; Hartmann, C. The next Generation of Glioma Biomarkers: MGMT Methylation, BRAF Fusions and IDH1 Mutations. Brain Pathol. 2011, 21, 74–87. [Google Scholar] [CrossRef] [PubMed]
  13. Palanca-Ballester, C.; Hervas, D.; Villalba, M.; Valdes-Sanchez, T.; Garcia, D.; Alcoriza-Balaguer, M.I.; Benet, M.; Martinez-Tomas, R.; Briones-Gomez, A.; Galbis-Caravajal, J.; et al. Translation of a Tissue Epigenetic Signature to Circulating Free DNA Suggests BCAT1 as a Potential Noninvasive Diagnostic Biomarker for Lung Cancer. Clin. Epigenetics 2022, 14, 116. [Google Scholar] [CrossRef] [PubMed]
  14. Meyer, B.; Clifton, S.; Locke, W.; Luu, P.-L.; Du, Q.; Lam, D.; Armstrong, N.J.; Kumar, B.; Deng, N.; Harvey, K.; et al. Identification of DNA Methylation Biomarkers with Potential to Predict Response to Neoadjuvant Chemotherapy in Triple-Negative Breast Cancer. Clin. Epigenetics 2021, 13, 226. [Google Scholar] [CrossRef]
  15. De Palma, F.D.E.; Luglio, G.; Tropeano, F.P.; Pagano, G.; D’Armiento, M.; Kroemer, G.; Maiuri, M.C.; De Palma, G.D. The Role of Micro-RNAs and Circulating Tumor Markers as Predictors of Response to Neoadjuvant Therapy in Locally Advanced Rectal Cancer. Int. J. Mol. Sci. 2020, 21, 7040. [Google Scholar] [CrossRef]
  16. Sarkar, D.; Diermeier, S.D. Circular RNAs: Potential Applications as Therapeutic Targets and Biomarkers in Breast Cancer. Noncoding RNA 2021, 7, 2. [Google Scholar] [CrossRef]
  17. De Palma, F.D.E.; Raia, V.; Kroemer, G.; Maiuri, M.C. The Multifaceted Roles of MicroRNAs in Cystic Fibrosis. Diagnostics 2020, 10, 1102. [Google Scholar] [CrossRef]
  18. De Palma, F.D.E.; D’Argenio, V.; Pol, J.; Kroemer, G.; Maiuri, M.C.; Salvatore, F. The Molecular Hallmarks of the Serrated Pathway in Colorectal Cancer. Cancers 2019, 11, 1017. [Google Scholar] [CrossRef] [Green Version]
  19. De Palma, F.D.E.; Salvatore, F.; Pol, J.G.; Kroemer, G.; Maiuri, M.C. Circular RNAs as Potential Biomarkers in Breast Cancer. Biomedicines 2022, 10, 725. [Google Scholar] [CrossRef]
  20. Cabús, L.; Lagarde, J.; Curado, J.; Lizano, E.; Pérez-Boza, J. Current Challenges and Best Practices for Cell-Free Long RNA Biomarker Discovery. Biomark. Res. 2022, 10, 62. [Google Scholar] [CrossRef]
  21. Di Cecilia, S.; Zhang, F.; Sancho, A.; Li, S.D.; Aguilo, F.; Sun, Y.; Rengasamy, M.; Zhang, W.; Vecchio, L.D.; Salvatore, F.; et al. RBM5-AS1 Is Critical for Self-Renewal of Colon Cancer Stem-like Cells. Cancer Res. 2016, 76, 5615–5627. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Uszczynska-Ratajczak, B.; Lagarde, J.; Frankish, A.; Guigó, R.; Johnson, R. Towards a Complete Map of the Human Long Non-Coding RNA Transcriptome. Nat. Rev. Genet. 2018, 19, 535–548. [Google Scholar] [CrossRef] [PubMed]
  23. Statello, L.; Guo, C.-J.; Chen, L.-L.; Huarte, M. Gene Regulation by Long Non-Coding RNAs and Its Biological Functions. Nat. Rev. Mol. Cell Biol. 2021, 22, 96–118. [Google Scholar] [CrossRef] [PubMed]
  24. Zampetaki, A.; Albrecht, A.; Steinhofel, K. Long Non-Coding RNA Structure and Function: Is There a Link? Front. Physiol. 2018, 9, 1201. [Google Scholar] [CrossRef]
  25. Qian, Y.; Shi, L.; Luo, Z. Long Non-Coding RNAs in Cancer: Implications for Diagnosis, Prognosis, and Therapy. Front. Med. 2020, 7, 612393. [Google Scholar] [CrossRef]
  26. Schmitt, A.M.; Chang, H.Y. Long Noncoding RNAs in Cancer Pathways. Cancer Cell 2016, 29, 452–463. [Google Scholar] [CrossRef] [Green Version]
  27. Hou, Z.; Zhao, W.; Zhou, J.; Shen, L.; Zhan, P.; Xu, C.; Chang, C.; Bi, H.; Zou, J.; Yao, X.; et al. A Long Noncoding RNA Sox2ot Regulates Lung Cancer Cell Proliferation and Is a Prognostic Indicator of Poor Survival. Int. J. Biochem. Cell Biol. 2014, 53, 380–388. [Google Scholar] [CrossRef]
  28. Kogo, R.; Shimamura, T.; Mimori, K.; Kawahara, K.; Imoto, S.; Sudo, T.; Tanaka, F.; Shibata, K.; Suzuki, A.; Komune, S.; et al. Long Noncoding RNA HOTAIR Regulates Polycomb-Dependent Chromatin Modification and Is Associated with Poor Prognosis in Colorectal Cancers. Cancer Res. 2011, 71, 6320–6326. [Google Scholar] [CrossRef] [Green Version]
  29. Okugawa, Y.; Toiyama, Y.; Hur, K.; Toden, S.; Saigusa, S.; Tanaka, K.; Inoue, Y.; Mohri, Y.; Kusunoki, M.; Boland, C.R.; et al. Metastasis-Associated Long Non-Coding RNA Drives Gastric Cancer Development and Promotes Peritoneal Metastasis. Carcinogenesis 2014, 35, 2731–2739. [Google Scholar] [CrossRef] [Green Version]
  30. Deng, H.-Y.; Wang, Y.-C.; Ni, P.-Z.; Lin, Y.-D.; Chen, L.-Q. Long Noncoding RNAs Are Novel Potential Prognostic Biomarkers for Esophageal Squamous Cell Carcinoma: An Overview. J. Thorac. Dis. 2016, 8, E653–E659. [Google Scholar] [CrossRef]
  31. Lu, C.; Wei, D.; Zhang, Y.; Wang, P.; Zhang, W. Long Non-Coding RNAs as Potential Diagnostic and Prognostic Biomarkers in Breast Cancer: Progress and Prospects. Front. Oncol. 2021, 11, 710538. [Google Scholar] [CrossRef] [PubMed]
  32. De Palma, F.D.E.; Del Monaco, V.; Pol, J.G.; Kremer, M.; D’Argenio, V.; Stoll, G.; Montanaro, D.; Uszczyńska-Ratajczak, B.; Klein, C.C.; Vlasova, A.; et al. The Abundance of the Long Intergenic Non-Coding RNA 01087 Differentiates between Luminal and Triple-Negative Breast Cancers and Predicts Patient Outcome. Pharmacol. Res. 2020, 161, 105249. [Google Scholar] [CrossRef] [PubMed]
  33. Flippot, R.; Malouf, G.G.; Su, X.; Mouawad, R.; Spano, J.-P.; Khayat, D. Cancer Subtypes Classification Using Long Non-Coding RNA. Oncotarget 2016, 7, 54082–54093. [Google Scholar] [CrossRef] [PubMed]
  34. Li, R.; Qian, J.; Wang, Y.; Zhang, J.; You, Y. Long Noncoding RNA Profiles Reveal Three Molecular Subtypes in Glioma. CNS Neurosci. Ther. 2014, 20, 339–343. [Google Scholar] [CrossRef]
  35. Su, X.; Malouf, G.G.; Chen, Y.; Zhang, J.; Yao, H.; Valero, V.; Weinstein, J.N.; Spano, J.-P.; Meric-Bernstam, F.; Khayat, D.; et al. Comprehensive Analysis of Long Non-Coding RNAs in Human Breast Cancer Clinical Subtypes. Oncotarget 2014, 5, 9864–9876. [Google Scholar] [CrossRef] [Green Version]
  36. Tian, Y.-H.; Jia, L.-W.; Liu, Z.-F.; Chen, Y.-H. LINC01087 Inhibits Glioma Cell Proliferation and Migration, and Increases Cell Apoptosis via MiR-384/Bcl-2 Axis. Aging 2021, 13, 20808–20819. [Google Scholar] [CrossRef]
  37. Chen, W.; Wang, F.; Zhang, J.; Li, C.; Hong, L. LINC01087 Indicates a Poor Prognosis of Glioma Patients with Preoperative MRI. Funct. Integr. Genom. 2022, 22, 55–64. [Google Scholar] [CrossRef]
  38. Yin, Y.; Huang, J.; Shi, H.; Huang, Y.; Huang, Z.; Song, M.; Yin, L. LINC01087 Promotes the Proliferation, Migration, and Invasion of Thyroid Cancer Cells by Upregulating PPM1E. J. Oncol. 2022, 2022, 7787378. [Google Scholar] [CrossRef]
  39. Zhang, J.; Xiao, F.; Qiang, G.; Zhang, Z.; Ma, Q.; Hao, Y.; Xing, H.; Liang, C. Novel LncRNA Panel as for Prognosis in Esophageal Squamous Cell Carcinoma Based on CeRNA Network Mechanism. Comput. Math. Methods Med. 2021, 2021, 8020879. [Google Scholar] [CrossRef]
  40. R: A Language and Environment for Statistical Computing. Available online: https://www.gbif.org/tool/81287/r-a-language-and-environment-for-statistical-computing (accessed on 3 June 2021).
  41. Nagy, Á.; Munkácsy, G.; Győrffy, B. Pancancer Survival Analysis of Cancer Hallmark Genes. Sci. Rep. 2021, 11, 6047. [Google Scholar] [CrossRef]
  42. Lánczky, A.; Győrffy, B. Web-Based Survival Analysis Tool Tailored for Medical Research (KMplot): Development and Implementation. J. Med. Internet Res. 2021, 23, e27633. [Google Scholar] [CrossRef] [PubMed]
  43. Robin, X.; Turck, N.; Hainard, A.; Tiberti, N.; Lisacek, F.; Sanchez, J.-C.; Müller, M. PROC: An Open-Source Package for R and S+ to Analyze and Compare ROC Curves. BMC Bioinform. 2011, 12, 77. [Google Scholar] [CrossRef] [PubMed]
  44. Venables, W.N.; Ripley, B.D. Modern Applied Statistics with S, Statistics and Computing, 4th ed.; Springer: New York, NY, USA, 2002; ISBN 978-0-387-95457-8. [Google Scholar]
  45. Karagkouni, D.; Paraskevopoulou, M.D.; Tastsoglou, S.; Skoufos, G.; Karavangeli, A.; Pierros, V.; Zacharopoulou, E.; Hatzigeorgiou, A.G. DIANA-LncBase v3: Indexing Experimentally Supported MiRNA Targets on Non-Coding Transcripts. Nucleic Acids Res. 2020, 48, D101–D110. [Google Scholar] [CrossRef] [PubMed]
  46. Huang, H.-Y.; Lin, Y.-C.-D.; Li, J.; Huang, K.-Y.; Shrestha, S.; Hong, H.-C.; Tang, Y.; Chen, Y.-G.; Jin, C.-N.; Yu, Y.; et al. MiRTarBase 2020: Updates to the Experimentally Validated MicroRNA–Target Interaction Database. Nucleic Acids Res. 2019, 48, D148–D154. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Calmon, M.F.; Jeschke, J.; Zhang, W.; Dhir, M.; Siebenkäs, C.; Herrera, A.; Tsai, H.-C.; O’Hagan, H.M.; Pappou, E.P.; Hooker, C.M.; et al. Epigenetic Silencing of Neurofilament Genes Promotes an Aggressive Phenotype in Breast Cancer. Epigenetics 2015, 10, 622–632. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Carlevaro-Fita, J.; Johnson, R. Global Positioning System: Understanding Long Noncoding RNAs through Subcellular Localization. Mol. Cell 2019, 73, 869–883. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Li, Y.; Kuscu, C.; Banach, A.; Zhang, Q.; Pulkoski-Gross, A.; Kim, D.; Liu, J.; Roth, E.; Li, E.; Shroyer, K.R.; et al. MiR-181a-5p Inhibits Cancer Cell Migration and Angiogenesis via Downregulation of Matrix Metalloproteinase-14. Cancer Res. 2015, 75, 2674–2685. [Google Scholar] [CrossRef] [Green Version]
  50. Ardizzone, A.; Calabrese, G.; Campolo, M.; Filippone, A.; Giuffrida, D.; Esposito, F.; Colarossi, C.; Cuzzocrea, S.; Esposito, E.; Paterniti, I. Role of MiRNA-19a in Cancer Diagnosis and Poor Prognosis. Int. J. Mol. Sci. 2021, 22, 4697. [Google Scholar] [CrossRef] [PubMed]
  51. Melling, N.; Harutyunyan, L.; Hube-Magg, C.; Kluth, M.; Simon, R.; Lebok, P.; Minner, S.; Tsourlakis, M.C.; Koop, C.; Graefen, M.; et al. High-Level HOOK3 Expression Is an Independent Predictor of Poor Prognosis Associated with Genomic Instability in Prostate Cancer. PLoS ONE 2015, 10, e0134614. [Google Scholar] [CrossRef] [Green Version]
  52. Huang, W.; Li, B.-R.; Feng, H. PLAG1 Silencing Promotes Cell Chemosensitivity in Ovarian Cancer via the IGF2 Signaling Pathway. Int. J. Mol. Med. 2020, 45, 703–714. [Google Scholar] [CrossRef] [PubMed]
  53. Xu, X.; Jaberi-Douraki, M.; Wallace, N.A. Predicting the Prognostic Value of POLI Expression in Different Cancers via a Machine Learning Approach. Int. J. Mol. Sci. 2022, 23, 8571. [Google Scholar] [CrossRef] [PubMed]
  54. Katabi, N.; Xu, B.; Jungbluth, A.A.; Zhang, L.; Shao, S.Y.; Lane, J.; Ghossein, R.; Antonescu, C.R. PLAG1 Immunohistochemistry Is a Sensitive Marker for Pleomorphic Adenoma: A Comparative Study with PLAG1 Genetic Abnormalities. Histopathology 2018, 72, 285. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Jiang, A.; Meng, J.; Gong, W.; Zhang, Z.; Gan, X.; Wang, J.; Wu, Z.; Liu, B.; Qu, L.; Wang, L. Elevated SNRPA1, as a Promising Predictor Reflecting Severe Clinical Outcome via Effecting Tumor Immunity for CcRCC, Is Related to Cell Invasion, Metastasis, and Sunitinib Sensitivity. Front. Immunol. 2022, 13, 842096. [Google Scholar] [CrossRef] [PubMed]
  56. Jiao, Z.; Feng, X.; Cui, Y.; Wang, L.; Gan, J.; Zhao, Y.; Meng, Q. Expression Characteristic, Immune Signature, and Prognosis Value of EFNA Family Identified by Multi-Omics Integrative Analysis in Pan-Cancer. BMC Cancer 2022, 22, 871. [Google Scholar] [CrossRef]
  57. Ebron, J.S.; Shankar, E.; Singh, J.; Sikand, K.; Weyman, C.M.; Gupta, S.; Lindner, D.J.; Liu, X.; Campbell, M.J.; Shukla, G.C. MiR-644a Disrupts Oncogenic Transformation and Warburg Effect by Direct Modulation of Multiple Genes of Tumor-Promoting Pathways. Cancer Res. 2019, 79, 1844–1856. [Google Scholar] [CrossRef] [Green Version]
  58. Heverhagen, A.E.; Legrand, N.; Wagner, V.; Fendrich, V.; Bartsch, D.K.; Slater, E.P. Overexpression of MicroRNA MiR-7-5p Is a Potential Biomarker in Neuroendocrine Neoplasms of the Small Intestine. Neuroendocrinology 2018, 106, 312–317. [Google Scholar] [CrossRef]
  59. Long, M.; Zhan, M.; Xu, S.; Yang, R.; Chen, W.; Zhang, S.; Shi, Y.; He, Q.; Mohan, M.; Liu, Q.; et al. MiR-92b-3p Acts as a Tumor Suppressor by Targeting Gabra3 in Pancreatic Cancer. Mol. Cancer 2017, 16, 167. [Google Scholar] [CrossRef] [Green Version]
  60. Wang, Q.; Yang, X.; Zhou, X.; Wu, B.; Zhu, D.; Jia, W.; Chu, J.; Wang, J.; Wu, J.; Kong, L. MiR-3174 Promotes Proliferation and Inhibits Apoptosis by Targeting FOXO1 in Hepatocellular Carcinoma. Biochem. Biophys. Res. Commun. 2020, 526, 889–897. [Google Scholar] [CrossRef]
  61. Wang, X.; Zhao, Y.; Lu, Q.; Fei, X.; Lu, C.; Li, C.; Chen, H. MiR-34a-5p Inhibits Proliferation, Migration, Invasion and Epithelial-Mesenchymal Transition in Esophageal Squamous Cell Carcinoma by Targeting LEF1 and Inactivation of the Hippo-YAP1/TAZ Signaling Pathway. J. Cancer 2020, 11, 3072–3081. [Google Scholar] [CrossRef]
  62. Yu, Z.; Zhao, H.; Feng, X.; Li, H.; Qiu, C.; Yi, X.; Tang, H.; Zhang, J. Long Non-Coding RNA FENDRR Acts as a MiR-423-5p Sponge to Suppress the Treg-Mediated Immune Escape of Hepatocellular Carcinoma Cells. Mol. Ther. Nucleic Acids 2019, 17, 516–529. [Google Scholar] [CrossRef]
  63. Zheng, Y.P.; Wu, L.; Gao, J.; Wang, Y. Tumor Suppressive Role of MiR-569 in Lung Cancer. Oncol. Lett. 2018, 15, 4087–4092. [Google Scholar] [CrossRef] [PubMed]
  64. Li, Z.; Hong, S.; Liu, Z. LncRNA LINC00641 Predicts Prognosis and Inhibits Bladder Cancer Progression through MiR-197-3p/KLF10/PTEN/PI3K/AKT Cascade. Biochem. Biophys. Res. Commun. 2018, 503, 1825–1829. [Google Scholar] [CrossRef] [PubMed]
  65. Lv, S.-Y.; Shan, T.-D.; Pan, X.-T.; Tian, Z.-B.; Liu, X.-S.; Liu, F.-G.; Sun, X.-G.; Xue, H.-G.; Li, X.-H.; Han, Y.; et al. The LncRNA ZEB1-AS1 Sponges MiR-181a-5p to Promote Colorectal Cancer Cell Proliferation by Regulating Wnt/β-Catenin Signaling. Cell Cycle 2018, 17, 1245–1254. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Cai, Y.; Hao, Y.; Ren, H.; Dang, Z.; Xu, H.; Xue, X.; Gao, Y. MiR-1305 Inhibits The Progression Of Non-Small Cell Lung Cancer By Regulating MDM2. Cancer Manag. Res. 2019, 11, 9529–9540. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Zhang, S.; Wan, H.; Zhang, X. LncRNA LHFPL3-AS1 Contributes to Tumorigenesis of Melanoma Stem Cells via the MiR-181a-5p/BCL2 Pathway. Cell Death Dis. 2020, 11, 950. [Google Scholar] [CrossRef]
  68. She, J.-K.; Fu, D.-N.; Zhen, D.; Gong, G.-H.; Zhang, B. LINC01087 Is Highly Expressed in Breast Cancer and Regulates the Malignant Behavior of Cancer Cells Through MiR-335-5p/Rock1. OncoTargets Ther. 2020, 13, 9771–9783. [Google Scholar] [CrossRef]
  69. Wang, H.; Zheng, Q.; Lu, Z.; Wang, L.; Ding, L.; Xia, L.; Zhang, H.; Wang, M.; Chen, Y.; Li, G. Role of the Nervous System in Cancers: A Review. Cell Death Discov. 2021, 7, 76. [Google Scholar] [CrossRef]
  70. Li, Y.-H.; Liu, Y.; Li, Y.-D.; Liu, Y.-H.; Li, F.; Ju, Q.; Xie, P.-L.; Li, G.-C. GABA Stimulates Human Hepatocellular Carcinoma Growth through Overexpressed GABAA Receptor Theta Subunit. World J. Gastroenterol. 2012, 18, 2704–2711. [Google Scholar] [CrossRef]
  71. Takehara, A.; Hosokawa, M.; Eguchi, H.; Ohigashi, H.; Ishikawa, O.; Nakamura, Y.; Nakagawa, H. γ-Aminobutyric Acid (GABA) Stimulates Pancreatic Cancer Growth through Overexpressing GABAA Receptor π Subunit. Cancer Res. 2007, 67, 9704–9712. [Google Scholar] [CrossRef] [Green Version]
  72. Watanabe, M.; Maemura, K.; Kanbara, K.; Tamayama, T.; Hayasaki, H. GABA and GABA Receptors in the Central Nervous System and Other Organs. In International Review of Cytology, A Survey of Cell Biology; Jeon, K.W., Ed.; Academic Press: Cambridge, MA, USA, 2002; Volume 213, pp. 1–47. [Google Scholar]
  73. Le, A.; Udupa, S.; Zhang, C. The Metabolic Interplay between Cancer and Other Diseases. Trends Cancer 2019, 5, 809–821. [Google Scholar] [CrossRef]
  74. Gumireddy, K.; Li, A.; Kossenkov, A.V.; Sakurai, M.; Yan, J.; Li, Y.; Xu, H.; Wang, J.; Zhang, P.J.; Zhang, L.; et al. The MRNA-Edited Form of GABRA3 Suppresses GABRA3-Mediated Akt Activation and Breast Cancer Metastasis. Nat. Commun. 2016, 7, 10715. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Keihani, S.; Kluever, V.; Mandad, S.; Bansal, V.; Rahman, R.; Fritsch, E.; Gomes, L.C.; Gärtner, A.; Kügler, S.; Urlaub, H.; et al. The Long Noncoding RNA NeuroLNC Regulates Presynaptic Activity by Interacting with the Neurodegeneration-Associated Protein TDP-43. Sci. Adv. 2019, 5, eaay2670. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Yuan, A.; Wu, P.; Zhong, Z.; He, Z.; Li, W. Long Non-Coding RNA Gm37494 Alleviates Osteoarthritis Chondrocyte Injury via the MicroRNA-181a-5p/GABRA1 Axis. J. Orthop. Surg. Res. 2022, 17, 304. [Google Scholar] [CrossRef] [PubMed]
  77. Zhang, J.-X.; Chen, Z.-H.; Xu, Y.; Chen, J.-W.; Weng, H.-W.; Yun, M.; Zheng, Z.-S.; Chen, C.; Wu, B.-L.; Li, E.-M.; et al. Downregulation of MicroRNA-644a Promotes Esophageal Squamous Cell Carcinoma Aggressiveness and Stem Cell–like Phenotype via Dysregulation of PITX2. Clin. Cancer Res. 2017, 23, 298–310. [Google Scholar] [CrossRef] [Green Version]
  78. Liang, W.; Liao, Y.; Li, Z.; Wang, Y.; Zheng, S.; Xu, X.; Ran, F.; Tang, B.; Wang, Z. MicroRNA-644a Promotes Apoptosis of Hepatocellular Carcinoma Cells by Downregulating the Expression of Heat Shock Factor 1. Cell Commun. Signal. 2018, 16, 30. [Google Scholar] [CrossRef] [Green Version]
  79. Wu, X.; Cheng, Y.-S.L.; Matthen, M.; Yoon, A.; Schwartz, G.K.; Bala, S.; Taylor, A.M.; Momen-Heravi, F. Down-Regulation of the Tumor Suppressor MiR-34a Contributes to Head and Neck Cancer by up-Regulating the MET Oncogene and Modulating Tumor Immune Evasion. J. Exp. Clin. Cancer Res. 2021, 40, 70. [Google Scholar] [CrossRef]
  80. Xiao, H. MiR-7-5p Suppresses Tumor Metastasis of Non-Small Cell Lung Cancer by Targeting NOVA2. Cell Mol. Biol. Lett. 2019, 24, 60. [Google Scholar] [CrossRef]
  81. Liu, Y.; Cheng, T.; Du, Y.; Hu, X.; Xia, W. LncRNA LUCAT1/MiR-181a-5p Axis Promotes Proliferation and Invasion of Breast Cancer via Targeting KLF6 and KLF15. BMC Mol. Cell Biol. 2020, 21, 69. [Google Scholar] [CrossRef]
  82. Lee, J.Y.; Ryu, D.; Lim, S.W.; Ryu, K.J.; Choi, M.E.; Yoon, S.E.; Kim, K.; Park, C.; Kim, S.J. Exosomal MiR-1305 in the Oncogenic Activity of Hypoxic Multiple Myeloma Cells: A Biomarker for Predicting Prognosis. J. Cancer 2021, 12, 2825–2834. [Google Scholar] [CrossRef]
  83. Wang, Y.; Yang, H.; Zhang, G.; Luo, C.; Zhang, S.; Luo, R.; Deng, B. Hsa-MiR-7-5p Suppresses Proliferation, Migration and Promotes Apoptosis in Hepatocellular Carcinoma Cell Lines by Inhibiting SPC24 Expression. Biochem. Biophys. Res. Commun. 2021, 561, 80–87. [Google Scholar] [CrossRef]
  84. Pan, Y.; Jin, K.; Xie, X.; Wang, K.; Zhang, H. MicroRNA-19a-3p Inhibits the Cellular Proliferation and Invasion of Non-small Cell Lung Cancer by Downregulating UBAP2L. Exp. Ther. Med. 2020, 20, 2252–2261. [Google Scholar] [CrossRef] [PubMed]
  85. Jiang, X.-M.; Yu, X.-N.; Liu, T.-T.; Zhu, H.-R.; Shi, X.; Bilegsaikhan, E.; Guo, H.-Y.; Song, G.-Q.; Weng, S.-Q.; Huang, X.-X.; et al. MicroRNA-19a-3p Promotes Tumor Metastasis and Chemoresistance through the PTEN/Akt Pathway in Hepatocellular Carcinoma. Biomed. Pharmacother. 2018, 105, 1147–1154. [Google Scholar] [CrossRef] [PubMed]
  86. Jiang, T.; Ye, L.; Han, Z.; Liu, Y.; Yang, Y.; Peng, Z.; Fan, J. MiR-19b-3p Promotes Colon Cancer Proliferation and Oxaliplatin-Based Chemoresistance by Targeting SMAD4: Validation by Bioinformatics and Experimental Analyses. J. Exp. Clin. Cancer Res. 2017, 36, 131. [Google Scholar] [CrossRef] [Green Version]
  87. Wei, Y.; Guo, S.; Tang, J.; Wen, J.; Wang, H.; Hu, X.; Gu, Q. MicroRNA-19b-3p Suppresses Gastric Cancer Development by Negatively Regulating Neuropilin-1. Cancer Cell Int. 2020, 20, 193. [Google Scholar] [CrossRef]
  88. Chaluvally-Raghavan, P.; Zhang, F.; Pradeep, S.; Hamilton, M.P.; Zhao, X.; Rupaimoole, R.; Moss, T.; Lu, Y.; Yu, S.; Pecot, C.V.; et al. Copy Number Gain of Hsa-MiR-569 at 3q26.2 Leads to Loss of TP53INP1 and Aggressiveness of Epithelial Cancers. Cancer Cell 2014, 26, 863–879. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Dai, Y.; Zhang, Y.; Hao, M.; Zhu, R. LINC00665 Functions as a Competitive Endogenous RNA to Regulate AGTR1 Expression by Sponging MiR-34a-5p in Glioma. Oncol. Rep. 2021, 45, 1202–1212. [Google Scholar] [CrossRef]
  90. Shao, T.; Hu, Y.; Tang, W.; Shen, H.; Yu, Z.; Gu, J. The Long Noncoding RNA HOTAIR Serves as a MicroRNA-34a-5p Sponge to Reduce Nucleus Pulposus Cell Apoptosis via a NOTCH1-Mediated Mechanism. Gene 2019, 715, 144029. [Google Scholar] [CrossRef] [PubMed]
  91. Ji, Y.; Wang, M.; Li, X.; Cui, F. The Long Noncoding RNA NEAT1 Targets MiR-34a-5p and Drives Nasopharyngeal Carcinoma Progression via Wnt/β-Catenin Signaling. Yonsei Med. J. 2019, 60, 336–345. [Google Scholar] [CrossRef]
  92. Gao, Y.; Luo, X.; Zhang, J. LincRNA-ROR Is Activated by H3K27 Acetylation and Induces EMT in Retinoblastoma by Acting as a Sponge of MiR-32 to Activate the Notch Signaling Pathway. Cancer Gene Ther. 2021, 28, 42–54. [Google Scholar] [CrossRef]
  93. Li, Q.-Y.; Shen, J.-Q.; Li, J.-H.; Dai, D.-F.; Saeed, M.; Li, C.-X. LINC00958/MiR-3174/PHF6 Axis Is Responsible for Triggering Proliferation, Migration and Invasion of Endometrial Cancer. Eur. Rev. Med. Pharm. Sci. 2021, 25, 6853–6861. [Google Scholar] [CrossRef]
  94. Li, Y.; Wang, X.; Zhao, Z.; Shang, J.; Li, G.; Zhang, R. LncRNA NEAT1 Promotes Glioma Cancer Progression via Regulation of MiR-98-5p/BZW1. Biosci. Rep. 2021, 41, BSR20200767. [Google Scholar] [CrossRef] [PubMed]
  95. Cilek, E.E.; Ozturk, H.; Dedeoglu, B.G. Construction of MiRNA-MiRNA Networks Revealing the Complexity of MiRNA-Mediated Mechanisms in Trastuzumab Treated Breast Cancer Cell Lines. PLoS ONE 2017, 12, e0185558. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Enokida, H.; Yoshino, H.; Matsushita, R.; Nakagawa, M. The Role of MicroRNAs in Bladder Cancer. Investig. Clin. Urol. 2016, 57, S60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. Koshizuka, K.; Hanazawa, T.; Fukumoto, I.; Kikkawa, N.; Okamoto, Y.; Seki, N. The MicroRNA Signatures: Aberrantly Expressed MicroRNAs in Head and Neck Squamous Cell Carcinoma. J. Hum. Genet. 2017, 62, 3–13. [Google Scholar] [CrossRef]
  98. Shams, R.; Saberi, S.; Zali, M.; Sadeghi, A.; Ghafouri-Fard, S.; Aghdaei, H.A. Identification of Potential MicroRNA Panels for Pancreatic Cancer Diagnosis Using Microarray Datasets and Bioinformatics Methods. Sci. Rep. 2020, 10, 7559. [Google Scholar] [CrossRef]
  99. Zhou, C.; Chen, Y.; He, X.; Zheng, Z.; Xue, D. Functional Implication of Exosomal MiR-217 and MiR-23b-3p in the Progression of Prostate Cancer. OncoTargets Ther. 2020, 13, 11595–11606. [Google Scholar] [CrossRef]
  100. Lin, C.-W.; Yang, W.-E.; Lee, W.-J.; Hua, K.-T.; Hsieh, F.-K.; Hsiao, M.; Chen, C.-C.; Chow, J.-M.; Chen, M.-K.; Yang, S.-F.; et al. Lipocalin 2 Prevents Oral Cancer Metastasis through Carbonic Anhydrase IX Inhibition and Is Associated with Favourable Prognosis. Carcinogenesis 2016, 37, 712–722. [Google Scholar] [CrossRef] [Green Version]
  101. He, S.; Zeng, S.; Zhou, Z.-W.; He, Z.-X.; Zhou, S.-F. Hsa-MicroRNA-181a Is a Regulator of a Number of Cancer Genes and a Biomarker for Endometrial Carcinoma in Patients: A Bioinformatic and Clinical Study and the Therapeutic Implication. Drug Des. Dev. Ther. 2015, 9, 1103–1175. [Google Scholar] [CrossRef] [Green Version]
  102. Sakiyama, T.; Kohno, T.; Mimaki, S.; Ohta, T.; Yanagitani, N.; Sobue, T.; Kunitoh, H.; Saito, R.; Shimizu, K.; Hirama, C.; et al. Association of Amino Acid Substitution Polymorphisms in DNA Repair Genes TP53, POLI, REV1 and LIG4 with Lung Cancer Risk. Int. J. Cancer 2005, 114, 730–737. [Google Scholar] [CrossRef]
  103. Hoff, A.M.; Kraggerud, S.M.; Alagaratnam, S.; Berg, K.C.G.; Johannessen, B.; Høland, M.; Nilsen, G.; Lingjærde, O.C.; Andrews, P.W.; Lothe, R.A.; et al. Frequent Copy Number Gains of SLC2A3 and ETV1 in Testicular Embryonal Carcinomas. Endocr. Relat. Cancer 2020, 27, 457–468. [Google Scholar] [CrossRef]
  104. Ziegler, G.C.; Almos, P.; McNeill, R.V.; Jansch, C.; Lesch, K.-P. Cellular Effects and Clinical Implications of SLC2A3 Copy Number Variation. J. Cell. Physiol. 2020, 235, 9021–9036. [Google Scholar] [CrossRef] [PubMed]
  105. Yao, X.; He, Z.; Qin, C.; Deng, X.; Bai, L.; Li, G.; Shi, J. SLC2A3 Promotes Macrophage Infiltration by Glycolysis Reprogramming in Gastric Cancer. Cancer Cell Int. 2020, 20, 503. [Google Scholar] [CrossRef] [PubMed]
  106. Gao, H.; Liang, J.; Duan, J.; Chen, L.; Li, H.; Zhen, T.; Zhang, F.; Dong, Y.; Shi, H.; Han, A. A Prognosis Marker SLC2A3 Correlates With EMT and Immune Signature in Colorectal Cancer. Front. Oncol. 2021, 11, 638099. [Google Scholar] [CrossRef]
  107. Chu, M.; Zheng, K.; Li, X.; Luo, Z.; Yang, X.; Wei, C. Comprehensive Analysis of the Role of SLC2A3 on Prognosis and Immune Infiltration in Head and Neck Squamous Cell Carcinoma. Anal. Cell. Pathol. 2022, 2022, e2371057. [Google Scholar] [CrossRef] [PubMed]
  108. Liu, J.; Hong, J.; Han, H.; Park, J.; Kim, D.; Park, H.; Ko, M.; Koh, Y.; Shin, D.-Y.; Yoon, S.-S. Decreased Vitamin C Uptake Mediated by SLC2A3 Promotes Leukaemia Progression and Impedes TET2 Restoration. Br. J. Cancer 2020, 122, 1445–1452. [Google Scholar] [CrossRef] [PubMed]
  109. Dai, D.-W.; Lu, Q.; Wang, L.-X.; Zhao, W.-Y.; Cao, Y.-Q.; Li, Y.-N.; Han, G.-S.; Liu, J.-M.; Yue, Z.-J. Decreased MiR-106a Inhibits Glioma Cell Glucose Uptake and Proliferation by Targeting SLC2A3 in GBM. BMC Cancer 2013, 13, 478. [Google Scholar] [CrossRef] [Green Version]
  110. Li, J.; Fu, F.; Wan, X.; Huang, S.; Wu, D.; Li, Y. Up-Regulated MiR-29c Inhibits Cell Proliferation and Glycolysis by Inhibiting SLC2A3 Expression in Prostate Cancer. Gene 2018, 665, 26–34. [Google Scholar] [CrossRef]
  111. Chen, D.; Wang, H.; Chen, J.; Li, Z.; Li, S.; Hu, Z.; Huang, S.; Zhao, Y.; He, X. MicroRNA-129-5p Regulates Glycolysis and Cell Proliferation by Targeting the Glucose Transporter SLC2A3 in Gastric Cancer Cells. Front Pharm. 2018, 9, 502. [Google Scholar] [CrossRef] [Green Version]
  112. Kania, A.; Klein, R. Mechanisms of Ephrin–Eph Signalling in Development, Physiology and Disease. Nat. Rev. Mol. Cell Biol. 2016, 17, 240–256. [Google Scholar] [CrossRef]
  113. Wang, T.-H.; Chang, J.-L.; Ho, J.-Y.; Wu, H.-C.; Chen, T.-C. EphrinA5 Suppresses Colon Cancer Development by Negatively Regulating Epidermal Growth Factor Receptor Stability. FEBS J. 2012, 279, 251–263. [Google Scholar] [CrossRef]
  114. Herath, N.I.; Spanevello, M.D.; Sabesan, S.; Newton, T.; Cummings, M.; Duffy, S.; Lincoln, D.; Boyle, G.; Parsons, P.G.; Boyd, A.W. Over-Expression of Eph and Ephrin Genes in Advanced Ovarian Cancer: Ephrin Gene Expression Correlates with Shortened Survival. BMC Cancer 2006, 6, 144. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  115. Giaginis, C.; Tsourouflis, G.; Zizi-Serbetzoglou, A.; Kouraklis, G.; Chatzopoulou, E.; Dimakopoulou, K.; Theocharis, S.E. Clinical Significance of Ephrin (Eph)-A1, -A2, -A4, -A5 and -A7 Receptors in Pancreatic Ductal Adenocarcinoma. Pathol. Oncol. Res. 2010, 16, 267–276. [Google Scholar] [CrossRef]
  116. Li, J.-J.; Liu, D.-P.; Liu, G.-T.; Xie, D. EphrinA5 Acts as a Tumor Suppressor in Glioma by Negative Regulation of Epidermal Growth Factor Receptor. Oncogene 2009, 28, 1759–1768. [Google Scholar] [CrossRef] [Green Version]
  117. Feng, J.; Guo, J.; Zhao, P.; Shen, J.; Chai, B.; Wang, J. MTOR Up-Regulation of SNRPA1 Contributes to Hepatocellular Carcinoma Development. Biosci. Rep. 2020, 40, BSR20193815. [Google Scholar] [CrossRef] [PubMed]
  118. Prognostic Value of Survival-Associated Splicing Factor SNRPA1 Overexpression and Its Potential Mechanism in Liver Cancer. Available online: https://www.fortunejournals.com/articles/prognostic-value-of-survivalassociated-splicing-factor-snrpa1-overexpression-and-its-potential-mechanism-in-liver-cancer.html (accessed on 1 August 2022).
  119. Zeng, Q.; Lei, F.; Chang, Y.; Gao, Z.; Wang, Y.; Gao, Q.; Niu, P.; Li, Q. An Oncogenic Gene, SNRPA1, Regulates PIK3R1, VEGFC, MKI67, CDK1 and Other Genes in Colorectal Cancer. Biomed Pharm. 2019, 117, 109076. [Google Scholar] [CrossRef] [PubMed]
  120. Chen, Y.; Hao, Q.; Wang, S.; Cao, M.; Huang, Y.; Weng, X.; Wang, J.; Zhang, Z.; He, X.; Lu, H.; et al. Inactivation of the Tumor Suppressor P53 by Long Noncoding RNA RMRP. Proc. Natl. Acad. Sci. USA 2021, 118, e2026813118. [Google Scholar] [CrossRef]
  121. Ding, B.; Yao, M.; Fan, W.; Lou, W. Whole-Transcriptome Analysis Reveals a Potential Hsa_circ_0001955/Hsa_circ_0000977-Mediated MiRNA-MRNA Regulatory Sub-Network in Colorectal Cancer. Aging 2020, 12, 5259–5279. [Google Scholar] [CrossRef]
  122. Siddappa, M.; Wani, S.A.; Long, M.D.; Leach, D.A.; Mathé, E.A.; Bevan, C.L.; Campbell, M.J. Identification of Transcription Factor Co-Regulators That Drive Prostate Cancer Progression. Sci. Rep. 2020, 10, 20332. [Google Scholar] [CrossRef]
  123. Estrada-Cuzcano, A.; Neveling, K.; Kohl, S.; Banin, E.; Rotenstreich, Y.; Sharon, D.; Falik-Zaccai, T.C.; Hipp, S.; Roepman, R.; Wissinger, B.; et al. Mutations in C8orf37, Encoding a Ciliary Protein, Are Associated with Autosomal-Recessive Retinal Dystrophies with Early Macular Involvement. Am. J. Hum. Genet. 2012, 90, 102–109. [Google Scholar] [CrossRef] [Green Version]
  124. Heon, E.; Kim, G.; Qin, S.; Garrison, J.E.; Tavares, E.; Vincent, A.; Nuangchamnong, N.; Scott, C.A.; Slusarski, D.C.; Sheffield, V.C. Mutations in C8ORF37 Cause Bardet Biedl Syndrome (BBS21). Hum. Mol. Genet. 2016, 25, 2283–2294. [Google Scholar] [CrossRef]
Cancers 14 05980 g001 550
Figure 1. Deregulation of LINC01087 expression in representative cancers. Boxplots depicting the expression level of LINC01087 in normal and malignant tissues derived from the TCGA and GTEx databases. Data show median, quartiles, and individual values that are expressed as gene raw count. Comparison was performed using Wilcoxon rank-sum test, *** p ≤ 0.001, **** p ≤ 0.0001. Abbreviations: BC, breast cancer (n = 1685); ESCA, esophageal carcinoma (n = 1619); GTEx, genotype-tissue expression; LUSC, lung squamous cell carcinoma (n = 551); OV, ovarian cancer (n = 561); STAD, stomach adenocarcinoma (n = 766); TCGA, the cancer genome atlas; TGCT, testicular germ cell tumors (n = 517); UCS, uterine carcinosarcoma (n = 199).
Figure 1. Deregulation of LINC01087 expression in representative cancers. Boxplots depicting the expression level of LINC01087 in normal and malignant tissues derived from the TCGA and GTEx databases. Data show median, quartiles, and individual values that are expressed as gene raw count. Comparison was performed using Wilcoxon rank-sum test, *** p ≤ 0.001, **** p ≤ 0.0001. Abbreviations: BC, breast cancer (n = 1685); ESCA, esophageal carcinoma (n = 1619); GTEx, genotype-tissue expression; LUSC, lung squamous cell carcinoma (n = 551); OV, ovarian cancer (n = 561); STAD, stomach adenocarcinoma (n = 766); TCGA, the cancer genome atlas; TGCT, testicular germ cell tumors (n = 517); UCS, uterine carcinosarcoma (n = 199).
Cancers 14 05980 g001
Cancers 14 05980 g002 550
Figure 2. LINC010187-related gene expression changes in representative cancers. (A) Venn diagram showing the overlapping genes positively associated with LINC01087 in cancers where it shows a strong diagnostic interest. The percentage of overlapping genes is reported in brackets. (B,C) Bubble maps depicting results of (B) GO term and (C) KEGG pathway enrichment analyses. The size of the circles indicates the number of genes in each set (count), while the color code informs about the adjusted p-values (p.adjust). Panel (B) illustrates the annotated clusters that included ≥ 4 genes. Panel (B) also indicates the “Biological process”, “Cellular component”, and “Molecular function” subcategories of GO hierarchy in turquoise, green, and fuchsia, respectively. See Supplementary Tables S7 and S8 for details. Abbreviations: BC, breast cancer; ESCA, esophageal carcinoma; GO, gene ontology; KEGG, Kyoto encyclopedia of genes and genomes; OV, ovarian cancer; STAD, stomach adenocarcinoma; TGCT, testicular germ cell tumors.
Figure 2. LINC010187-related gene expression changes in representative cancers. (A) Venn diagram showing the overlapping genes positively associated with LINC01087 in cancers where it shows a strong diagnostic interest. The percentage of overlapping genes is reported in brackets. (B,C) Bubble maps depicting results of (B) GO term and (C) KEGG pathway enrichment analyses. The size of the circles indicates the number of genes in each set (count), while the color code informs about the adjusted p-values (p.adjust). Panel (B) illustrates the annotated clusters that included ≥ 4 genes. Panel (B) also indicates the “Biological process”, “Cellular component”, and “Molecular function” subcategories of GO hierarchy in turquoise, green, and fuchsia, respectively. See Supplementary Tables S7 and S8 for details. Abbreviations: BC, breast cancer; ESCA, esophageal carcinoma; GO, gene ontology; KEGG, Kyoto encyclopedia of genes and genomes; OV, ovarian cancer; STAD, stomach adenocarcinoma; TGCT, testicular germ cell tumors.
Cancers 14 05980 g002
Cancers 14 05980 g003 550
Figure 3. Schematic representation of the strategy used to characterize the LINC01087-miRNA-mRNA network. A Spearman correlation test was applied to TCGA datasets to identify 2458 transcripts (Analysis 1) and 64 miRNAs (Analysis 2), the expression of which correlated with that of LINC01087 in BC, ESCA, OV, STAD, and TGCT. In parallel, the DIANA-LncBASE v3 toolbox allowed to discover 20 miRNAs that interacted with LINC01087 (Analysis 3). Next, the miRTarBase predicted 2396 mRNA targets of the miRNAs identified in Analysis 2, and 5532 mRNA targets of the miRNAs identified in Analysis 3. The expression of 31 among these mRNA targets from Analysis 2 and 68 mRNA targets from Analysis 3 had been correlated with that of LINC01087 in Analysis 1. Among them, 10 messenger RNAs were shared between the three analyses, hence representing the core of the genes regulated by LINC0187 and related miRNAs. The resulting LINC01087-miRNA-mRNA network is detailed in Table 3.
Figure 3. Schematic representation of the strategy used to characterize the LINC01087-miRNA-mRNA network. A Spearman correlation test was applied to TCGA datasets to identify 2458 transcripts (Analysis 1) and 64 miRNAs (Analysis 2), the expression of which correlated with that of LINC01087 in BC, ESCA, OV, STAD, and TGCT. In parallel, the DIANA-LncBASE v3 toolbox allowed to discover 20 miRNAs that interacted with LINC01087 (Analysis 3). Next, the miRTarBase predicted 2396 mRNA targets of the miRNAs identified in Analysis 2, and 5532 mRNA targets of the miRNAs identified in Analysis 3. The expression of 31 among these mRNA targets from Analysis 2 and 68 mRNA targets from Analysis 3 had been correlated with that of LINC01087 in Analysis 1. Among them, 10 messenger RNAs were shared between the three analyses, hence representing the core of the genes regulated by LINC0187 and related miRNAs. The resulting LINC01087-miRNA-mRNA network is detailed in Table 3.
Cancers 14 05980 g003
Table
Table 1. Diagnostic and prognostic value of LINC01087 in a variety of human cancers.
Table 1. Diagnostic and prognostic value of LINC01087 in a variety of human cancers.
DiagnosisPrognosis
TCGA and GTExTCGA KM Plotter
ExpressionClinical Interest (ROC)OS
Type of TumorAbbreviationN. of Samples (N, T)p-Valuelog2FCp-ValueAUCp-ValueHR
(95%CI) 		p-ValueHR
(95%CI)
Adrenocortical carcinomaACC337
(258, 79)		NS−0.170.030.6NS1.15
(0.53–2.5)		NANA
Bladder cancerBLCA452
(40, 412)		NS2.50.030.6NS1.01
(0.75–1.36)		NS1
(0.74–1.37)
Breast cancerBC1685
(572, 1113)		2.32E-352.21.16E-350.70.010.65
(0.47–0.9)		0.00630.64
(0.46–0.88)
Cervical squamous cell carcinoma and endocervical adenocarcinomaCESC328
(22, 306)		NS0.2NS0.5NS0.91
(0.57–1.44)		NS0.92
(0.56–1.51)
CholangiocarcinomaCHOL270
(235, 35)		NS2.5NS0.5NS1.24
(0.46–3.38)		NANA
Colon adenocarcinomaCOAD1300
(820, 480)		NS0.80.050.5NS0.97
(0.61–1.52)		NANA
Esophageal carcinomaESCA1619
(1456, 163)		4.08E-353.72.04E-350.8NS1.50
(0.92–2.45)		NS1.7
(0.89–3.23)
Glioblastoma multiformeGBM2816
(2647, 169)		NS−0.1NS0.5NS1.09
(0.76–1.57)		NANA
Lower grade gliomaLGG3174
(2642, 532)		2.43E-04−0.30.00010.5NS0.98
(0.63–1.52)		NANA
Head-neck squamous cell carcinomaHNSC710
(206, 504)		1.55E-110.17.77E-120.6NS0.92
(0.69–1.21)		NS0.8
(0.59–1.09)
Hepatocellular carcinomaHCC650
(276, 374)		4.28E-070.82.14E-070.6NS1.38
(0.96–1.98)		0.0371.47
(1.02–2.13)
Acute myeloid leukemiaAML1079
(929, 150)		9.46E-08−14.73E-080.6NS1.26
(0.82–1.93)		NANA
Lung adenocarcinomaLUAD1176
(637, 539)		0.0021.30.0010.5NS0.99
(0.73–1.35)		NS1.05
(0.77–1.44)
squamous cell carcinoma LUSC551
(49, 502)		0.0032.70.0010.6NS0.98
(0.75–1.29)		NS124
(0.94–1.63)
Skin cutaneous melanomaSKCM2282
(1810, 472)		7.99E-230.53.99E-230.6NS0.84
(0.64–1.12)		NANA
Ovarian cancerOV561
(180, 381)		1.11E-141.55.56E-150.7NS0.97
(0.74–1.26)		NS0.95
(0.73–1.24)
Pancreatic adenocarcinomaPAAD511
(332, 179)		6.86E-20−1.13.43E-200.7NS0.69
(0.45–1.06)		NS0.84
(0.53–1.34)
Pheochromocytoma and ParagangliomaPCPG445
(261,184)		3.87E-07−0.71.93E-070.6NS3.3e-09
(0–Inf)		NS0
(0–Inf)
Prostate adenocarcinomaPRAD798
(297, 501)		NS1.6NS0.5NS0.94
(0.25–3.50)		NANA
Rectum adenocarcinomaREAD956
(789,167)		0.0010.74.73E-040.6NS0.82
(0.37–1.82)		NS0.67
(0.29–1.53)
Chromophobe renal cell carcinomachRCC179
(114, 65)		NS0.1NS0.5NS0.45
(0.05–3.56)		NANA
Clear cell renal cell carcinomaccRCC702
(161, 541)		1.25E-150.86.24E-160.7NS1.20
(0.89–1.63)		NS1.19
(0.88–1.61)
Papillary renal cell carcinoma pRCC412
(121, 291)		0.0050.90.0020.62.86E-042.85
(1.58–5.17)		0.00252.46
(1.35–4.49)
SarcomaSARC265
(2, 263)		NS0.6NS0.6NS0.66
(0.43–1.02)		NS0.71
(0.46–1.12)
Stomach adenocarcinomaSTAD766
(391, 375)		2.49E-5041.25E-500.8NS1.02
(0.73–1.41)		NS1.05
(0.76–1.47)
Testicular germ cell tumorsTGCT517
(361, 156)		6.59E-73−6.13.30E-731.0NS0.76
(0.07–8.38)		NS0.85
(0.08–9.44)
ThymomaTHYM121
(2, 119)		NS5.20.0470.8NS0.91
(0.24–3.40)		NS1.1
(0.3–4.4)
Thyroid carcinomaTHCA1224
(712, 512)		3.13E-10−0.61.56E-100.6NS1.84
(0.67–5.08)		NS1.05
(0.36–3.02)
Uterine corpus endometrial carcinomaUCEC731
(177, 554)		NS1.1NS0.5NS1.48
(0.97–2.24)		0.00881.73
(1.14–2.62)
Uterine carcinosarcomaUCS199
(142, 57)		0.00013.26.69E-050.6NS1.57
(0.80–3.08)		NANA
All prognostic data extracted from the Kaplan–Meier plotter (KM-Plotter) database were generated using RNA-sequencing pan-cancer datasets. Significant diagnostic and prognostic interests of LINC01087 are highlighted in bold. Abbreviations: AUC, area under the curve; CI, confidence interval; FC, fold change; GTEx, genotype-tissue expression; HR, hazard ratio; KM, Kaplan–Meier; N, normal samples; NA, not available; NS, not significant; OS, overall survival; RFS, relapse-free survival; ROC, receiver operating characteristic; T, tumor samples; TCGA, the cancer genome atlas.
Table
Table 2. Significant Associations between LINC01087 Expression and Clinicopathological Characteristics in TCGA and GTEx data.
Table 2. Significant Associations between LINC01087 Expression and Clinicopathological Characteristics in TCGA and GTEx data.
BC
Clinicopathological Featuresn. of Total CasesLINC01087
Expression		p-Value
LowHigh
pTT1-T29294534760.04
T3-T418010377
Tumor size≤2 cm2821231590.01
>2 cm827433394
TGCT
Clinicopathological Featuresn. of Total CasesLINC01087
Expression		p-Value
Tumor size≤2 cm8041390.05
>2 cm584018
Abbreviations: BC, breast cancer; pN, lymph node stage; pM, metastasis stage; pT, tumor stage; TGCT, testicular germ cell tumors. Fisher’s test.
Table
Table 3. LINC01087-miRNA-mRNA network identified in the present study.
Table 3. LINC01087-miRNA-mRNA network identified in the present study.
lncRNA10 miRNA Targets Extracted
from Analysis 2		18 miRNA Targets Extracted
from Analysis 3		10 mRNA Targets Shared between
Analyses 1, 2, and 3
LINC01087hsa-miR-644ahsa-miR-7-5pC8orf37
hsa-miR-1305hsa-miR-92a-3pEFNA5
hsa-miR-4505hsa-miR-181a-5pHOOK3
hsa-miR-3671hsa-miR-181a-5p, hsa-miR-181b-5p, hsa-miR-181d-5pIQCG
hsa-miR-7975hsa-miR-197-3pMYOZ3
hsa-miR-4505hsa-miR-423-5pPCP4L1
hsa-miR-569hsa-miR-181a-5p, hsa-miR-181b-5p, hsa-miR-98-5pPLAG1
hsa-miR-4255hsa-miR-19a-3p, hsa-miR-19b-3p, hsa-miR-92a-3pPOLI
hsa-miR-606hsa-miR-34a-5p, hsa-miR-98-5pSLC2A3
hsa-miR-3174hsa-miR-34a-5p, hsa-miR-98-5pSNRPA1
In Analysis 2, miRNA targets were evidenced using TCGA; In Analysis 3, miRNA targets were evidenced using DIANA-LncBASE. In both Analyses 2 and 3, mRNA targets of LINC01087-related miRNAs were identified using miRTarBase.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

De Palma, F.D.E.; Carbonnier, V.; Salvatore, F.; Kroemer, G.; Pol, J.G.; Maiuri, M.C. Systematic Investigation of the Diagnostic and Prognostic Impact of LINC01087 in Human Cancers. Cancers 2022, 14, 5980. https://doi.org/10.3390/cancers14235980

AMA Style

De Palma FDE, Carbonnier V, Salvatore F, Kroemer G, Pol JG, Maiuri MC. Systematic Investigation of the Diagnostic and Prognostic Impact of LINC01087 in Human Cancers. Cancers. 2022; 14(23):5980. https://doi.org/10.3390/cancers14235980

Chicago/Turabian Style

De Palma, Fatima Domenica Elisa, Vincent Carbonnier, Francesco Salvatore, Guido Kroemer, Jonathan G. Pol, and Maria Chiara Maiuri. 2022. "Systematic Investigation of the Diagnostic and Prognostic Impact of LINC01087 in Human Cancers" Cancers 14, no. 23: 5980. https://doi.org/10.3390/cancers14235980

APA Style

De Palma, F. D. E., Carbonnier, V., Salvatore, F., Kroemer, G., Pol, J. G., & Maiuri, M. C. (2022). Systematic Investigation of the Diagnostic and Prognostic Impact of LINC01087 in Human Cancers. Cancers, 14(23), 5980. https://doi.org/10.3390/cancers14235980

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop