Integrative analysis of ceRNA network reveals functional lncRNAs associated with independent recurrent prognosis in colon adenocarcinoma.
Ontology highlight
ABSTRACT: Background
Long non-coding RNAs (lncRNAs), acting as competing endogenous RNA (ceRNA) have been reported to regulate the expression of targeted genes by sponging miRNA in colon adenocarcinoma (COAD).Methods
However, their potential implications for recurrence free survival prognosis and functional roles remains largely unclear in COAD. In this study, we downloaded the TCGA dataset (training dataset) and GSE39582 (validation dataset) of COAD patients with prognostic information.Results
A total of 411 differentially expressed genes (DElncRNAs: 12 downregulated and 43 upregulated), 18 DE miRNAs (9 downregulated and 9 upregulated) and 338 DEmRNAs (113 downregulated and 225 upregulated) were identified in recurrence samples compared with non-recurrence samples with the thresholds of FDR < 0.05 and |log2FC|> 0.263. Based on six signature lncRNAs (LINC00899, LINC01503, PRKAG2-AS1, RAD21-AS1, SRRM2-AS1 and USP30-AS1), the risk score (RS) system was constructed. Two prognostic clinical features, including pathologic stage and RS model status were screened for building the nomogram survival model. Moreover, a recurrent-specific ceRNA network was successfully constructed with 2 signature lncRNAs, 4 miRNAs and 113 mRNAs. Furthermore, we further manifested that SRRM2-AS1 predicted a poor prognosis in COAD patients. Furthermore, knockdown of SRRM2-AS1 significantly suppressed cell proliferation, migration, invasion and EMT markers in HT-29 and SW1116 cells.Conclusion
These identified novel lncRNA signature and ceRNA network associated with recurrence prognosis might provide promising therapeutic targets for COAD patients.
SUBMITTER: Mao Y
PROVIDER: S-EPMC8259330 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Intrahepatic cholangiocarcinoma (ICC) is the second most common lethal liver cancer worldwide. Currently, despite the latest developments in genomics and transcriptomics for ICC in recent years, the molecular pathogenesis promoting ICC remains elusive, especially in regulatory mechanisms of long noncoding RNAs (lncRNAs), which acts as competing endogenous RNA (ceRNA). In order to elucidate the molecular mechanism of functional lncRNA, expression profiles of lncRNAs, microRNAs (miRNAs), and messenger RNAs (mRNAs) were obtained from The Cancer Genome Atlas (TCGA) database and an integrative analysis of the ICC-associated ceRNA network was performed. Moreover, gene oncology enrichment analyses for the genes in the ceRNA network were implemented and novel prognostic biomarker lncRNA molecules were identified. In total, 6,738 differentially expressed mRNAs (DEmRNAs), 2,768 lncRNAs (DElncRNAs), and 173 miRNAs (DEmiRNAs) were identified in tumor tissues and adjacent nontumor ICC tissues with the thresholds of adjusted P < 0.01 and |logFC| > 2. An ICC-specific ceRNA network was successfully constructed with 30 miRNAs, 16 lncRNAs, and 80 mRNAs. Gene oncology enrichment analyses revealed that they were associated with the adaptive immune response, T cell selection and positive regulation of GTPase activity categories. Among the ceRNA networks, DElncRNAs ARHGEF26-AS1 and MIAT were found to be hub genes in underexpressed and overexpressed networks, respectively. Notably, univariate Cox regression analysis indicated that DElncRNAs HULC significantly correlated with overall survival (OS) in ICC patients (P value?<?0.05), and an additional survival analysis for HULC was reconfirmed in an independent ICC cohort from the Gene Expression Omnibus (GEO) database. These findings contribute to a more comprehensive understanding of the ICC-specific ceRNA network and provide novel strategies for subsequent functional studies of lncRNAs in ICC.
| S-EPMC6899321 | biostudies-literature
Project description:BackgroundThe survival prognosis is the hallmark of cancer progression. Here, we aimed to develop a recurrence-related gene signature to predict the prognosis of colon adenocarcinoma (COAD).MethodsThe proteomic data from the Clinical Proteomic Tumor Analysis Consortium (CPTAC) and genomic data from the cancer genomic maps [The Cancer Genome Atlas (TCGA)] dataset were analyzed to identify co-differentially expressed genes (cDEGs) between recurrence samples and non-recurrence samples in COAD using limma package. Functional enrichment analysis, including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway was conducted. Univariate and multivariate Cox regressions were applied to identify the independent prognostic feature cDEGs and establish the signature whose performance was evaluated by Kaplan-Meier curve, receiver operating characteristic (ROC), Harrell's concordance index (C-index), and calibration curve. The area under the receiver operating characteristic (ROC) curve (AUROC) and a nomogram were calculated to assess the predictive accuracy. GSE17538 and GSE39582 were used for external validation. Quantitative real-time PCR and Western blot analysis were carried out to validate our findings.ResultsWe identified 86 cDEGs in recurrence samples compared with non-recurrence samples. These genes were primarily enriched in the regulation of carbon metabolic process, fructose and mannose metabolism, and extracellular exosome. Then, an eight-gene-based signature (CA12, HBB, NCF1, KBTBD11, MMAA, DMBT1, AHNAK2, and FBLN2) was developed to separate patients into high- and low-risk groups. Patients in the low-risk group had significantly better prognosis than those in the high-risk group. Four prognostic clinical features, including pathological M, N, T, and RS model status, were screened for building the nomogram survival model. The PCR and Western blot analysis results suggested that CA12 and AHNAK2 were significantly upregulated, while MMAA and DMBT1 were downregulated in the tumor sample compared with adjacent tissues, and in non-recurrent samples compared with non-recurrent samples in COAD.ConclusionThese identified recurrence-related gene signatures might provide an effective prognostic predictor and promising therapeutic targets for COAD patients.
| S-EPMC9281446 | biostudies-literature
Project description:Mounting evidence suggests that long noncoding RNAs (lncRNAs) play important roles in the regulation of gene expression by acting as competing endogenous RNA (ceRNA). However, the regulatory mechanisms of lncRNA as ceRNA in gastric cancer (GC) are not fully understood. Here, we first constructed a dysregulated lncRNA-associated ceRNA network by integrating analysis of gene expression profiles of lncRNAs, microRNAs (miRNAs), and messenger RNAs (mRNAs). Then, we determined three lncRNAs (RP5-1120P11, DLEU2, and DDX11-AS1) as hub lncRNAs, in which associated ceRNA subnetworks were involved in cell cycle-related processes and cancer-related pathways. Furthermore, we confirmed that the two lncRNAs (DLEU2 and DDX11-AS1) were significantly upregulated in GC tissues, promote GC cell proliferation, and negatively regulate miRNA expression, respectively. The hub lncRNAs (DLEU2 and DDX11-AS1) could have oncogenic functions, and act as potential ceRNAs to sponge miRNA. Our findings not only provide novel insights on ceRNA regulation in GC, but can also provide opportunities for the functional characterization of lncRNAs in future studies.
| S-EPMC6027299 | biostudies-literature
Project description:BackgroundLong noncoding RNAs (lncRNAs) play important roles in competing endogenous RNA (ceRNA) networks involved in the development and progression of various cancers, including muscle-invasive bladder cancer (MIBC).PurposeThis study aims to construct the lncRNA-associated ceRNA network and identify lncRNA signatures correlated with the clinical features of MIBC tissue samples from The Cancer Genome Atlas (TGCA) database.MethodsThe differential expression profiles of MIBC associated lncRNAs, miRNAs and mRNAs were obtained from TCGA. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed to determine the principal functions of significantly dysregulated mRNAs. The dysregulated lncRNA-associated ceRNA network of MIBC was constructed based on the bioinformatics data, and the correlations between lncRNA expression and clinical features were analyzed using a weighted gene coexpression network analysis (WGCNA). Six cancer specific lncRNAs from the ceRNA network were randomly selected to detect their expression in 32 paired MIBC tissue samples and 5 bladder cancer cell lines using quantitative real-time polymerase chain reaction (qRT-PCR).ResultsThe ceRNA network was constructed with 30 lncRNAs, 13 miRNAs and 32 mRNAs. Seventeen lncRNAs in the ceRNA network correlated with certain clinical features, and only 1 lncRNA (MIR137HG) correlated with the overall survival (OS) of patients with MIBC (log-rank test P<0.05). GO and KEGG analyses revealed roles for the potential mRNA targets of MIR137HG in epithelial cell differentiation and the peroxisome proliferator-activated receptor (PPAR) and tumor necrosis factor (TNF) signaling pathways. The expression data from TCGA were highly consistent with the verification results of the MIBC tissue samples and bladder cancer cell lines.ConclusionThese findings improve our understanding of the regulatory mechanism of the lncRNA-miRNA-mRNA ceRNA network and reveal potential lncRNAs as prognostic biomarkers of MIBC.
| S-EPMC6614857 | biostudies-literature
Project description:Atopic dermatitis (AD) is a prevalent inflammatory skin disease characterized by its chronic nature and relapse. Ample evidence suggests that non-coding RNAs play a major role in AD pathogenesis. However, the mechanism remains unknown, particularly in AD recurrence. Dynamic morphological and cytokine changes were measured throughout the whole course of an FITC-induced AD recurrence murine model. Microarray assay and integrative analysis were performed to comprehensively explore long non-coding RNA (lncRNA), messenger RNA (mRNA), and microRNA (miRNA) networks. Our results showed that an AD recurrence model was established. Overall, 5766 lncRNAs, 4025 mRNAs, and 202 miRNAs changed after elicitation, whereas, 419 lncRNAs, 349 mRNAs, and more notably, only 23 miRNAs, were dysregulated in the remission phase. Gene ontology (GO) and KEGG pathway enrichment analyses were used to investigate the potential functions of the dysregulated genes. The altered regulation of seven miRNAs and seven lncRNAs were validated in different stages of the model. The competing endogenous RNA (ceRNA) network inferred that lncRNA humanlincRNA0490+ could compete for miR-155-5p binding, through which it might affect Pki? expression. Altogether, our findings have provided a novel perspective on the potential roles of non-coding RNAs in AD, and suggest that specific non-coding RNAs could be new therapeutic targets against AD recurrence.
| S-EPMC6214145 | biostudies-literature
Project description:Long noncoding RNAs (lncRNAs) have important biological functions as competing endogenous RNAs (ceRNAs) in tumors, yet the functions and regulatory mechanisms of lncRNA-related ceRNAs in gastric cancer have not been fully elucidated. In this study, we constructed a lncRNA-miRNA-mRNA ceRNA network and identified potential lncRNA biomarkers in gastric cancer. Basing on the RNA profiles downloaded from The Cancer Genome Atlas (TCGA) platform, the gastric cancer-specific differentially expressed lncRNAs, miRNAs, and mRNAs were screened for constructing a ceRNA network using bioinformatic tools. The enrichment analysis of the biological processes in Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes pathways was performed on the ceRNA-related DEmRNAs. According to the modularization of protein-protein interaction (PPI) network, we extracted a ceRNA subnetwork and analyzed the correlation between the expression of the lncRNAs involved and specific clinical features of patients. Next, the expression of highly up-regulated in liver cancer (HULC) and RP11-314B1.2 showed significant changes in several pathological processes involved in gastric cancer, and nine lncRNAs were found to be correlated with the overall survival of patients with gastric cancer. Through the univariate and multivariate Cox regression analyses, two lncRNAs (LINC00106 and RP11-999E24.3) were identified and utilized to establish a risk score model for assessing the prognosis of patients. The analysis results were also partially verified using quantitative real-time PCR. The findings from this study indicate that HULC, RP11-314B1.2, LINC00106, and RP11-999E24.3 could be considered as potential therapeutic targets or prognostic biomarkers in gastric cancer, and provide a new perspective for cancer pathogenesis research.
| S-EPMC7050084 | biostudies-literature
Project description:The altered expression of long non-coding RNAs (lncRNAs) has been implicated in the development and human diseases. However, functional roles and regulatory mechanisms of lncRNA as competing endogenous RNAs (ceRNAs) in osteoporosis and their potential clinical implication for osteoporosis risk are largely unexplored. In this study, we performed integrated analysis for paired expression profiles and regulatory relationships of dysregulated lncRNAs, mRNAs, and miRNAs based on "ceRNA hypothesis," and constructed an osteoporosis-related dysregulated miRNA-mediated lncRNA-mRNA ceRNA network (DysCeNet) composed of 105 nodes (including eight miRNAs, 24 mRNAs, and 73 lncRNAs) and 515 edges. Functional analysis suggested that the DysCeNet was involved in known osteoporosis or bone metabolism-related biological processes and pathways. Then, we performed random forest-based feature selection for 73 lncRNAs with ceRNA activity and identified 25 of 73 lncRNAs as potential diagnostic biomarkers. A random forest-based classifier composed of 25 lncRNA biomarkers (RF-25lncRNA) was developed for predicting osteoporosis risk. Performance evaluation with the leave-one-out cross-validation (LOOCV) procedure showed that the RF-25lncRNA achieved a good performance in distinguishing high- and low-bone mineral density (BMD) subjects in different osteoporosis datasets. Our study for the first time revealed a global view of lncRNA-associated ceRNA regulation in osteoporosis and provided novel lncRNAs with ceRNA activity as candidate epigenetic diagnostic biomarkers for early detection of osteoporosis risk.
| S-EPMC7136400 | biostudies-literature
Project description:BackgroundLung adenocarcinoma (LUAD) is the most common subtype of nonsmall-cell lung cancer (NSCLC) and has a high incidence rate and mortality. The survival of LUAD patients has increased with the development of targeted therapeutics, but the prognosis of these patients is still poor. Long noncoding RNAs (lncRNAs) play an important role in the occurrence and development of LUAD. The purpose of this study was to identify novel abnormally regulated lncRNA-microRNA (miRNA)-messenger RNA (mRNA) competing endogenous RNA (ceRNA) networks that may suggest new therapeutic targets for LUAD or relate to LUAD prognosis.MethodsWe used the SBC human ceRNA array V1.0 to screen for differentially expressed (DE) lncRNAs and mRNAs in four paired LUAD samples. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to annotate the DE lncRNAs and mRNAs. R bioinformatics packages, The Cancer Genome Atlas (TCGA) LUAD database, and Kaplan-Meier (KM) survival analysis tools were used to validate the microarray data and construct the lncRNA-miRNA-mRNA ceRNA regulatory network. Then, quantitative real-time PCR (qRT-PCR) was used to validate the DE lncRNAs in 7 LUAD cell lines.ResultsA total of 2819 DE lncRNAs and 2396 DE mRNAs (P < 0.05 and fold change ≥ 2 or ≤ 0.5) were identified in four paired LUAD tissue samples. In total, 255 of the DE lncRNAs were also identified in TCGA. The GO and KEGG analysis results suggested that the DE genes were most enriched in angiogenesis and cell proliferation, and were closely related to human cancers. Moreover, the differential expression of ENST00000609697, ENST00000602992, and NR_024321 was consistent with the microarray data, as determined by qRT-PCR validation in 7 LUAD cell lines; however, only ENST00000609697 was associated with the overall survival of LUAD patients (log-rank P = 0.029). Finally, through analysis of ENST00000609697 target genes, we identified the ENST00000609697-hsa-miR-6791-5p-RASL12 ceRNA network, which may play a tumor-suppressive role in LUAD.ConclusionENST00000609697 was abnormally expressed in LUAD. Furthermore, downregulation of ENST00000609697 and its target gene RASL12 was associated with poor prognosis in LUAD. The ENST00000609697-hsa-miR-6791-5p-RASL12 axis may play a tumor-suppressive role. These results suggest new potential prognostic and therapeutic biomarkers for LUAD.
| S-EPMC8611860 | biostudies-literature
Project description:The aberrant expression of long noncoding RNAs (lncRNAs) has drawn increasing attention in the field of hepatocellular carcinoma (HCC) biology. In the present study, we obtained the expression profiles of lncRNAs, microRNAs (miRNAs), and messenger RNAs (mRNAs) in 371 HCC tissues and 50 normal tissues from The Cancer Genome Atlas (TCGA) and identified hepatocarcinogenesis-specific differentially expressed genes (DEGs, log fold change ≥ 2, FDR < 0.01), including 753 lncRNAs, 97 miRNAs, and 1,535 mRNAs. Because the specific functions of lncRNAs are closely related to their intracellular localizations and because the cytoplasm is the main location for competitive endogenous RNA (ceRNA) action, we analyzed not only the interactions among these DEGs but also the distributions of lncRNAs (cytoplasmic, nuclear or both). Then, an HCC-associated deregulated ceRNA network consisting of 37 lncRNAs, 10 miRNAs, and 26 mRNAs was constructed after excluding those lncRNAs located only in the nucleus. Survival analysis of this network demonstrated that 15 lncRNAs, 3 miRNAs, and 16 mRNAs were significantly correlated with the overall survival of HCC patients (p < 0.01). Through multivariate Cox regression and lasso analysis, a risk score system based on 13 lncRNAs was constructed, which showed good discrimination and predictive ability for HCC patient survival time. This ceRNA network-construction approach, based on lncRNA distribution, not only narrowed the scope of target lncRNAs but also provided specific candidate molecular biomarkers for evaluating the prognosis of HCC, which will help expand our understanding of the ceRNA mechanisms involved in the early development of HCC.
| S-EPMC6618076 | biostudies-literature
Project description:Background Cancer metabolism is largely altered compared to normal cells. This study aims to explore critical metabolism pathways in colon adenocarcinoma (COAD), and reveal the possible mechanism of their role in cancer progression. Methods Expression data and sequencing data of COAD samples were obtained from The Cancer Genome Atlas and Gene Expression Omnibus databases. The expression profiles between tumor and normal samples were compared to identify differential metabolism pathways through single sample gene set enrichment analysis. Results Fatty acid synthesis was identified as a key metabolism pathway in COAD. Based on fatty acid-related lncRNAs, two molecular subtypes (C1 and C2) were defined. C2 subtype with worse prognosis had higher immune infiltration and higher expression of immune checkpoints. Five transcription factors (TFs) including FOS, JUN, HIF1A, STAT3 and STAT2 were highly expressed in C2 subtype. Five fatty acid-related lncRNAs were identified to be biomarkers for predicting COAD prognosis. Finally, further experients showed that knockdown of lncRNA PAXIP1-AS1 decreased the triglyceride content and the fatty acid synthase and acetyl-CoA carboxylase 1 expressions, which suggested that lncRNA PAXIP1-AS1 plays an important role in fatty acid metabolism of COAD. Conclusions This study demonstrated that fatty acid synthesis was greatly altered in COAD. Fatty acid-related lncRNAs were speculated to be involved in cancer progression through associating with TFs. The five screened TFs may serve as new drug targets for treating COAD. Supplementary Information The online version contains supplementary material available at 10.1186/s13062-022-00332-y.
| S-EPMC9331591 | biostudies-literature