scholarly journals Integrated analysis of lymphocyte infiltration-associated lncRNA for ovarian cancer via TCGA, GTEx and GEO datasets

PeerJ ◽  
2020 ◽  
Vol 8 ◽  
pp. e8961 ◽  
Author(s):  
Meijing Wu ◽  
Xiaobin Shang ◽  
Yue Sun ◽  
Jing Wu ◽  
Guoyan Liu

Background Abnormal expression of long non-coding RNAs (lncRNA) play a significant role in the incidence and progression of high-grade serous ovarian cancer (HGSOC), which is a leading cause of mortality among gynecologic malignant tumor patients. In this study, our aim is to identify lncRNA-associated competing endogenous RNA (ceRNA ) axes that could define more reliable prognostic parameters of HGSOC, and to investigate the lncRNAs’ potential mechanism of in lymphocyte infiltration. Methods The RNA-seq and miRNA expression profiles were downloaded from The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression (GTEx) database; while for obtaining the differentially expressed lncRNAs (DELs), miRNAs (DEMs), and genes (DEGs), we used edgeR, limma and DESeq2. After validating the RNA, miRNA and gene expressions, using integrated three RNA expression profiles (GSE18520, GSE27651, GSE54388) and miRNA profile (GSE47841) from the Gene Expression Omnibus (GEO) database, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway analyses through ClusterProfiler. The prognostic value of these genes was determined with Kaplan–Meier survival analysis and Cox regression analysis. The ceRNA network was constructed using Cytoscape. The correlation between lncRNAs in ceRNA network and immune infiltrating cells was analyzed by using Tumor IMmune Estimation Resource (TIMER), and gene markers of tumor-infiltrating immune cells were identified using Spearman’s correlation after removing the influence of tumor purity. Results A total of 33 DELs (25 upregulated and eight downregulated), 134 DEMs (76 upregulated and 58 downregulated), and 1,612 DEGs (949 upregulated and 663 downregulated) were detected that could be positively correlated with overall survival (OS) of HGSOC. With the 1,612 analyzed genes, we constructed a ceRNA network, which indicated a pre-dominant involvement of the immune-related pathways. Furthermore, our data revealed that LINC00665 influenced the infiltration level of macrophages and dendritic cells (DCs). On the other hand, FTX and LINC00665, which may play their possible roles through the ceRNA axis, demonstrated a potential to inhibit Tregs and prevent T-cell exhaustion. Conclusion We defined several prognostic biomarkers for the incidence and progression of HGSOC and constructed a network for ceRNA axes; among which three were indicated to have a positive correlation with lymphocyte infiltration, namely: FTX-hsa-miR-150-5p-STK11, LINC00665-hsa-miR449b-5p-VAV3 and LINC00665-hsa-miR449b-5p-RRAGD.


2019 ◽  
Vol 2019 ◽  
pp. 1-9 ◽  
Author(s):  
Dongkai Zhou ◽  
Bingqiang Gao ◽  
Qifan Yang ◽  
Yang Kong ◽  
Weilin Wang

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.



2020 ◽  
Vol 11 ◽  
Author(s):  
Jian-Rong Sun ◽  
Chen-Fan Kong ◽  
Kun-Min Xiao ◽  
Jia-Lu Yang ◽  
Xiang-Ke Qu ◽  
...  

Hepatocellular carcinoma (HCC) is one of the most common types of malignancy and is associated with high mortality. Prior research suggests that long non-coding RNAs (lncRNAs) play a crucial role in the development of HCC. Therefore, it is necessary to identify lncRNA-associated therapeutic biomarkers to improve the accuracy of HCC prognosis. Transcriptomic data of HCC obtained from The Cancer Genome Atlas (TCGA) database were used in the present study. Differentially expressed RNAs (DERNAs), including 74 lncRNAs, 16 miRNAs, and 35 mRNAs, were identified using bioinformatics analysis. The DERNAs were subsequently used to reconstruct a competing endogenous RNA (ceRNA) network. A lncRNA signature was revealed using Cox regression analysis, including LINC00200, MIR137HG, LINC00462, AP002478.1, and HTR2A-AS1. Kaplan-Meier plot demonstrated that the lncRNA signature is highly accurate in discriminating high- and low-risk patients (P &lt; 0.05). The area under curve (AUC) value exceeded 0.7 in both training and validation cohort, suggesting a high prognostic potential of the signature. Furthermore, multivariate Cox regression analysis indicated that both the TNM stage and the lncRNA signature could serve as independent prognostic factors for HCC (P &lt; 0.05). Then, a nomogram comprising the TNM stage and the lncRNA signature was determined to raise the accuracy in predicting the survival of HCC patients. In the present study, we have introduced a ceRNA network that could contribute to provide a new insight into the identification of potential regulation mechanisms for the development of HCC. The five-lncRNA signature could serve as a reliable biosignature for HCC prognosis, while the nomogram possesses strong potential in clinical applications.



PeerJ ◽  
2020 ◽  
Vol 8 ◽  
pp. e10437
Author(s):  
Xinnan Zhao ◽  
Miao He

Background Ovarian cancer (OC) is a highly malignant disease with a poor prognosis and high recurrence rate. At present, there is no accurate strategy to predict the prognosis and recurrence of OC. The aim of this study was to identify gene-based signatures to predict OC prognosis and recurrence. Methods mRNA expression profiles and corresponding clinical information regarding OC were collected from The Cancer Genome Atlas (TCGA) database. Gene set enrichment analysis (GSEA) and LASSO analysis were performed, and Kaplan–Meier curves, time-dependent ROC curves, and nomograms were constructed using R software and GraphPad Prism7. Results We first identified several key signalling pathways that affected ovarian tumorigenesis by GSEA. We then established a nine-gene-based signature for overall survival (OS) and a five-gene-based-signature for relapse-free survival (RFS) using LASSO Cox regression analysis of the TCGA dataset and validated the prognostic value of these signatures in independent GEO datasets. We also confirmed that these signatures were independent risk factors for OS and RFS by multivariate Cox analysis. Time-dependent ROC analysis showed that the AUC values for OS and RFS were 0.640, 0.663, 0.758, and 0.891, and 0.638, 0.722, 0.813, and 0.972 at 1, 3, 5, and 10 years, respectively. The results of the nomogram analysis demonstrated that combining two signatures with the TNM staging system and tumour status yielded better predictive ability. Conclusion In conclusion, the two-gene-based signatures established in this study may serve as novel and independent prognostic indicators for OS and RFS.



2021 ◽  
Author(s):  
Yuan Xu ◽  
Guofu Lin ◽  
Yifei Liu ◽  
Xianbin Lin ◽  
Hai Lin ◽  
...  

Abstract Background: Accumulating evidence indicates that long non-coding RNAs (lncRNAs) are involving in the tumorigenesis and metastasis of lung cancer. The aim of the study is to systematically characterize the lncRNA-associated competing endogenous RNA (ceRNA) network and identify key lncRNAs in the development of stage I lung adenocarcinoma (LUAD). Methods: Totally, 1,955 DEmRNAs, 165 DEmiRNAs and 1,107 DElncRNAs were obtained in 10 paired normal and LUAD tissues. And a total of 8,912 paired lncRNA-miRNA-mRNA network was constructed. Using the Cancer Genome Atlas (TCGA) dataset, the module of ME turquoise was revealed to be most relevant to the progression of LUAD though Weighted Gene Co-expression Network Analysis (WGCNA). Results: Of the lncRNAs identified, LINC00639, RP4-676L2.1 and FENDRR were in ceRNA network established by our RNA-sequencing dataset. Using univariate Cox regression analysis, FENDRR was a risk factor of progression free survival (PFS) of stage I LUAD patients (HRs=1.69, 95%CI 1.07-2.68, P< .050). Subsequently, differential expression of FENDRR in paired normal and LUAD tissues was detected significant by real-time quantitative (qRT-PCR) (P <0.001). Conclusions: This study, for the first time, deciphered the regulatory role of FENDRR/miR-6815-5p axis in the progression of early-stage LUAD, which is needed to be established in vitro and in vivo.



2021 ◽  
Author(s):  
Zhuoqi Li ◽  
Jing Zhou ◽  
Liankun Gu ◽  
Baozhen Zhang

Abstract Colorectal cancer (CRC) is one of the most common and deadly malignant carcinomas. Many long noncoding RNAs (lncRNA) have been reported to play an important role in the tumorigenesis of CRC by interacting with miRNAs and influencing the expression of some mRNAs through a competing endogenous RNA (ceRNA) network. Pseudogenes are one kind of lncRNA and can act as RNA sponges for miRNAs and regulate gene expression via ceRNA networks, but there are few studies about pseudogenes in CRC. In this study, total of 31 differentially expressed (DE) pseudogenes, 17 DE miRNAs and 152 DE mRNAs were identified by analyzing the expression profiles of colon adenocarcinoma (COAD) obtained from The Cancer Genome Atlas (TCGA). And a ceRNA network was constructed based on these RNAs. Kaplan–Meier analysis showed that 7 pseudogenes, 4 miRNAs and 30 mRNAs were significantly associated with overall survival. Then multivariate Cox regression analysis on the ceRNA-related DE pseudogenes was performed and a 5-pseudogene signature with the greatest prognostic value for CRC was identified. What’s more, the results were validated by the Gene Expression Omnibus (GEO) database, and quantitative real‐time PCR (qRT‐PCR) in 113 pairs of CRC tissues. In conclusion, this study provides a pseudogene-associated ceRNA network and 7 prognostic pseudogene biomarkers, and a 5-pseudogene prognostic risk signature that may be useful to predict the survival of CRC patients.



2020 ◽  
Vol 40 (8) ◽  
Author(s):  
Sihan Chen ◽  
Guodong Cao ◽  
Wei Wu ◽  
Yida Lu ◽  
Xiaobo He ◽  
...  

Abstract Colon adenocarcinoma (COAD) is a malignant gastrointestinal tumor, often occurring in the left colon, which is regulated by glycolysis-related processes. In past studies, multiple genes that influence the prognosis for survival have been discovered through bioinformatics analysis. However, the prediction of disease prognosis using a single gene is not an accurate method. In the present study, a mechanistic model was established to achieve better prediction for the prognosis of COAD. COAD-related data downloaded from The Cancer Genome Atlas (TCGA) were correlated with the glycolysis process using gene set enrichment analysis (GSEA) to determine the glycolysis-related genes that regulate COAD. Using COX regression analysis, glycolysis-related genes associated with the prognosis of COAD were identified, and the genes screened to establish a predictive model. The risk scores of this model were correlated with relevant clinical data to obtain a connection diagram between the model and survival rate, tumor characteristic data, etc. Finally, genes in the model were correlated with cells in the tumor microenvironment, finding that they affected specific immune cells in the model. Seven genes related to glycolysis were identified (PPARGC1A, DLAT, 6PC2, P4HA1, STC2, ANKZF1, and GPC1), which affect the prognosis of patients with COAD and constitute the model for prediction of survival of COAD patients.



2021 ◽  
Vol 12 ◽  
Author(s):  
Min Zhou ◽  
Shasha Hong ◽  
Bingshu Li ◽  
Cheng Liu ◽  
Ming Hu ◽  
...  

Background: DNA methylation affects the development, progression, and prognosis of various cancers. This study aimed to identify DNA methylated-differentially expressed genes (DEGs) and develop a methylation-driven gene model to evaluate the prognosis of ovarian cancer (OC).Methods: DNA methylation and mRNA expression profiles of OC patients were downloaded from The Cancer Genome Atlas, Genotype-Tissue Expression, and Gene Expression Omnibus databases. We used the R package MethylMix to identify DNA methylation-regulated DEGs and built a prognostic signature using LASSO Cox regression. A quantitative nomogram was then drawn based on the risk score and clinicopathological features.Results: We identified 56 methylation-related DEGs and constructed a prognostic risk signature with four genes according to the LASSO Cox regression algorithm. A higher risk score not only predicted poor prognosis, but also was an independent poor prognostic indicator, which was validated by receiver operating characteristic (ROC) curves and the validation cohort. A nomogram consisting of the risk score, age, FIGO stage, and tumor status was generated to predict 3- and 5-year overall survival (OS) in the training cohort. The joint survival analysis of DNA methylation and mRNA expression demonstrated that the two genes may serve as independent prognostic biomarkers for OS in OC.Conclusion: The established qualitative risk score model was found to be robust for evaluating individualized prognosis of OC and in guiding therapy.



2020 ◽  
Author(s):  
Dai Zhang ◽  
Si Yang ◽  
Yiche Li ◽  
Meng Wang ◽  
Jia Yao ◽  
...  

Abstract Background: Ovarian cancer (OV) is deemed as the most lethal gynecological cancer in women. The aim of this study was construct an effective gene prognostic model for OV patients.Methods: The expression profiles of glycolysis-related genes (GRGs) and clinical data of patients with OV were extracted from The Cancer Genome Atlas (TCGA) database. Univariate, multivariate, and least absolute shrinkage and selection operator Cox regression analyses were conducted, and a prognostic signature based on GRGs was constructed. The predictive ability of the signature was analyzed in training and test sets.Results: Based on nine GRGs (ISG20, CITED2, PYGB, IRS2, ANGPTL4, TGFBI, LHX9, PC, and DDIT4), a gene risk signature was identified to predict the outcome of patients with OV. The signature showed a good prognostic ability for OV, particularly high-grade OV, in the TCGA dataset, with areas under the curve of 0.709, 0.762, and 0.808 for 3-, 5- and 10-year survival, respectively. Similar results were found in the test sets, and the signature was also an independent prognostic factor. Moreover, a nomogram combining the prediction model and clinical factors was constructed.Conclusion: Our study established a nine-GRG risk model and a nomogram to better perform on OV patients’ survival prediction. The risk model represents a promising and independent prognostic predictor for OV patients. Moreover, our study of GRGs could offer guidances for underlying mechanisms explorations in the future.



2020 ◽  
Author(s):  
Gaochen Lan ◽  
Xiaoling Yu ◽  
Yanna Zhao ◽  
Jinjian Lan ◽  
Wan Li ◽  
...  

Abstract Background: Breast cancer is the most common malignant disease among women. At present, more and more attention has been paid to long non-coding RNAs (lncRNAs) in the field of breast cancer research. We aimed to investigate the expression profiles of lncRNAs and construct a prognostic lncRNA for predicting the overall survival (OS) of breast cancer.Methods: The expression profiles of lncRNAs and clinical data with breast cancer were obtained from The Cancer Genome Atlas (TCGA). Differentially expressed lncRNAs were screened out by R package (limma). The survival probability was estimated by the Kaplan‑Meier Test. The Cox Regression Model was performed for univariate and multivariate analysis. The risk score (RS) was established on the basis of the lncRNAs’ expression level (exp) multiplied regression coefficient (β) from the multivariate cox regression analysis with the following formula: RS=exp a1 * β a1 + exp a2 * β a2 +……+ exp an * β an. Functional enrichment analysis was performed by Metascape.Results: A total of 3404 differentially expressed lncRNAs were identified. Among them, CYTOR, MIR4458HG and MAPT-AS1 were significantly associated with the survival of breast cancer. Finally, The RS could predict OS of breast cancer (RS=exp CYTOR * β CYTOR + exp MIR4458HG * β MIR4458HG + exp MAPT-AS1 * β MAPT-AS1). Moreover, it was confirmed that the three-lncRNA signature could be an independent prognostic biomarker for breast cancer (HR=3.040, P=0.000).Conclusions: This study established a three-lncRNA signature, which might be a novel prognostic biomarker for breast cancer.



2021 ◽  
Author(s):  
Dan Bai ◽  
Huhu Feng ◽  
Jiajun Yang ◽  
Haitao Shi ◽  
Fanpu Ji ◽  
...  

Abstract Early diagnosis and prognosis rely on the successful identification of biomarkers and understanding the mechanisms. Through pan-cancer analysis amongst three types of gastrointestinal tumors, we constructed competitive endogenous RNA (ceRNA) networks, differentially expressed set of genes were distinguished, validated, and analyzed, their relevance to survival elucidated with immune infiltration profiles. Shared genes in esophageal, gastric and colon cancers were found significantly enriched in the processes of cell cycle, cell differentiation, DNA replication, synaptic transmission, the cyclic guanosine monophosphate protein-dependent protein kinase (cGMP-PKG) signaling pathway, and the glutamate receptor and other functions. Principal component analysis of the ceRNA network suggested the expression patterns of identified genes. Using Cox regression analysis of mRNAs, miRNAs and lncRNAs in the ceRNA network, genes including hsa-mir-196b, hsa-mir-584, PPP1R12B, SYNM, PDE2A, ALDH6A1 and MIR22HG were found significantly survival-related. Receiver operating characteristic (ROC) curve analysis and decision curve analysis (DCA) showed that the identified genes were not related to the survival in thyroid nor breast cancers, but effective for the prognosis of gastrointestinal tumors. These results could provide new theoretical and experimental clues, unravel the mechanisms, assisting molecular diagnosis and prognosis of gastrointestinal tumors.



Sign in / Sign up

Export Citation Format

Share Document