Genome‐wide screening and cohorts validation identifying novel lncRNAs as prognostic biomarkers for clear cell renal cell carcinoma

An Integrative Analysis Framework for Identifying the Prognostic Markers from Multidimensional RNA Data of Clear Cell Renal Cell Carcinoma

The altered regulatory status of long noncoding RNA (lncRNA), miRNA, and mRNA and their interactions play critical roles in tumor proliferation, metastasis, and progression, which ultimately influence cancer prognosis. However, there are limited studies of comprehensive identification of prognostic biomarkers from combined data sets of the three RNA types in the highly metastatic clear cell renal cell carcinoma (ccRCC). The current study employed an integrative analysis framework of functional genomics approaches and machine learning methods to the lncRNA, miRNA, and mRNA data and identified 16 RNAs (3 lncRNAs, 6 miRNAs, and 7 mRNAs) of prognostic value, with 9 of them novel. A 16 RNA-based score was established for prognosis prediction of ccRCC with significance (P < 0.0001). The area under the curve for the score model was 0.868 to 0.870 in the training cohort and 0.714 to 0.778 in the validation cohort. Construction of the lncRNA-miRNA-mRNA interaction network showed that the downstream mRNAs and upstream lncRNAs in the network initiated from the miRNA or lncRNA markers exhibit significant enrichment in functional classifications associated with cancer metastasis, proliferation, progression, or prognosis. The functional analysis provided clear support for the role of the RNA biomarkers in predicting cancer prognosis. This study provides promising biomarkers for predicting prognosis of ccRCC using multidimensional RNA data, and these findings are expected to facilitate potential clinical applications of the biomarkers.

lncRNA HIF1A-AS2: A potential oncogene in human cancers (Review)

Long non-coding RNAs (lncRNAs) are transcripts that are >200 nucleotides, but with no open reading frame. An increasing number of lncRNAs have been identified following the development of second-generation sequencing technologies, and they have since become a research hotspot. Functionally, they play a vital role in tumor progression, including in tumor proliferation, migration, invasion, apoptosis and acquisition of drug resistance. They regulate gene expression primarily through interaction with DNA, RNA and proteins at the epigenetic, transcriptional and post-transcriptional levels. Endogenous hypoxia-inducible factor 1α antisense RNA 2 (lncRNA HIF1A-AS2) is aberrantly expressed and involved the development/progression of various types of tumors, such as bladder cancer, glioblastoma, breast cancer and osteosarcoma. It plays a vital role in the proliferation, apoptosis, migration, invasion and epithelial-mesenchymal transformation of various tumor cells. This review summarizes the current body of knowledge on the biological functions and related molecular mechanisms of lncRNA HIF1A-AS2 in the development/progression of human tumors and other diseases.

Construction of Competitive Endogenous RNA Network and Verification of 3-Key LncRNA Signature Associated With Distant Metastasis and Poor Prognosis in Patients With Clear Cell Renal Cell Carcinoma

Clear cell renal cell carcinoma (ccRCC) is a common malignancy with high distant metastasis rate. Long non-coding RNAs (LncRNAs) are reported to be upregulated or downregulated in multiple cancers and play a crucial role in the metastasis of tumors or prognosis. Therefore, the purpose of our study is to construct a prognostic signature for ccRCC based on distant metastasis-related lncRNAs and explore the involved potential competitive endogenous RNA (ceRNA) network. The differentially expressed genes (DEGs) screened from the database of the cancer genome atlas (TCGA) were used to construct a co-expression network and identify the distant metastasis-related module by weighted gene co-expression network analysis (WGCNA). Key genes with metastatic and prognostic significance were identified through rigorous screening, including survival analysis, correlation analysis, and expression analyses in stage, grade, and distant metastasis, and were verified in the data set of gene expression omnibus (GEO) and the database from gene expression profiling interactive analysis (GEPIA). The potential upstream miRNAs and lncRNAs were predicted via five online databases and LncBase. Here, we constructed a ceRNA network of key genes that are significantly associated with the distant metastasis and prognosis of patients with ccRCC. The distant metastasis-related lncRNAs were used to construct a risk score model through the univariate, least absolute shrinkage selection operator (LASSO), and multivariate Cox regression analyses, and the patients were divided into high- and low-risk groups according to the median of the risk score. The Kaplan–Meier survival analysis demonstrated that mortality was significantly higher in the high-risk group than in the low-risk group. Considering the other clinical phenotype, the Cox regression analyses indicated that the lncRNAs model could function as an independent prognostic factor. Quantitative real-time (qRT)-PCR in the tissues and cells of ccRCC verified the high-expression level of three lncRNAs. Gene set enrichment analysis (GSEA) revealed that the lncRNA prognostic signature was mainly enriched in autophagy- and immune-related pathways, indicating that the autophagy and immune functions may play an important role in the distant metastasis of ccRCC. In summary, the constructed distant metastasis-related lncRNA signature could independently predict prognosis in patients with ccRCC, and the related ceRNA network provided a new sight on the potential mechanism of distant metastasis and a promising therapeutic target for ccRCC.

A panel of eight autophagy-related long non-coding RNAs is a good predictive parameter for clear cell renal cell carcinoma

Clear-cell renal cell carcinoma (ccRCC) carries a variable prognosis. Prognostic biomarkers can stratify patients according to risk, and can provide crucial information for clinical decision-making. We screened for an autophagy-related long non-coding lncRNA (lncRNA) signature to improve postoperative risk stratification in The Cancer Genome Atlas (TCGA) database. We confirmed this model in ICGC and SYSU cohorts as a significant and independent prognostic signature. Western blotting, autophagic-flux assay and transmission electron microscopy were used to verify that regulation of expression of 8 lncRNAs related to autophagy affected changes in autophagic flow in vitro. Our data suggest that 8-lncRNA signature related to autophagy is a promising prognostic tool in predicting the survival of patients with ccRCC. Combination of this signature with clinical and pathologic parameters could aid accurate risk assessment to guide clinical management, and this 8-lncRNAs signature related to autophagy may serve as a therapeutic target.

Integrative Analysis of DNA Methylation Identified 12 Signature Genes Specific to Metastatic ccRCC

Background: Abnormal epigenetic alterations can contribute to the development of human malignancies. Identification of these alterations for early screening and prognosis of clear cell renal cell carcinoma (ccRCC) has been a highly sought-after goal. Bioinformatic analysis of DNA methylation data provides broad prospects for discovery of epigenetic biomarkers. However, there is short of exploration of methylation-driven genes of ccRCC.

Methods: Gene expression data and DNA methylation data in metastatic ccRCC were sourced from the Gene Expression Omnibus (GEO) database. Differentially methylated genes (DMGs) at 5′-C-phosphate-G- 3′ (CpG) sites and differentially expressed genes (DEGs) were screened and the overlapping genes in DMGs and DEGs were then subject to gene set enrichment analysis. Next, the weighted gene co-expression network analysis (WGCNA) was used to search hub DMGs associated with ccRCC. Cox regression and ROC analyses were performed to screen potential biomarkers and develop a prognostic model based on the screened hub genes.

Results: Three hundred and fourteen overlapping DMGs were obtained from two independent GEO datasets. The turquoise module contained 79 hub DMGs, which represent the most significant module screened by WGCNA. Furthermore, a total of 12 hub genes (CETN3, DCAF7, GPX4, HNRNPA0, NUP54, SERPINB1, STARD5, TRIM52, C4orf3, C12orf51, and C17orf65) were identified in the TCGA database by multivariate Cox regression analyses. All the 12 genes were then used to generate the model for diagnosis and prognosis of ccRCC. ROC analysis showed that these genes exhibited good diagnostic efficiency for metastatic and non-metastatic ccRCC. Furthermore, the prognostic model with the 12 methylation-driven genes demonstrated a good prediction of 5-year survival rates for ccRCC patients.

Conclusion: Integrative analysis of DNA methylation data identified 12 signature genes, which could be used as epigenetic biomarkers for prognosis of metastatic ccRCC. This prognostic model has a good prediction of 5-year survival for ccRCC patients.

Long non-coding RNA LINC01234 regulates proliferation, migration and invasion via HIF-2α pathways in clear cell renal cell carcinoma cells

Long non-coding RNAs (lncRNAs) have been proved to have an important role in different malignancies including clear cell renal cell carcinoma (ccRCC). However, their role in disease progression is still not clear. The objective of the study was to identify lncRNA-based prognostic biomarkers and further to investigate the role of one lncRNA LINC01234 in progression of ccRCC cells. We found that six adverse prognostic lncRNA biomarkers including LINC01234 were identified in ccRCC patients by bioinformatic analysis using The Cancer Genome Atlas database. LINC01234 knockdown impaired cell proliferation, migration and invasion in vitro as compared to negative control. Furthermore, the epithelial-mesenchymal transition was inhibited after LINC01234 knockdown. Additionally, LINC01234 knockdown impaired hypoxia-inducible factor-2a (HIF-2α) pathways, including a suppression of the expression of HIF-2α, vascular endothelial growth factor A, epidermal growth factor receptor, c-Myc, Cyclin D1 and MET. Together, these datas showed that LINC01234 was likely to regulate the progression of ccRCC by HIF-2α pathways, and LINC01234 was both a promising prognostic biomarker and a potential therapeutic target for ccRCC.

Identification of a three-long noncoding RNA prognostic model involved competitive endogenous RNA in kidney renal clear cell carcinoma

Background

Long noncoding RNA (lncRNA) is generally identified as competing endogenous RNA (ceRNA) that plays a vital role in the pathogenesis of kidney renal clear cell carcinoma (KIRC), the most common subtype of renal cell carcinoma with poor prognosis and unclear pathogenesis. This study established a novel ceRNA network and thus identified a three-lncRNA prognostic model in KIRC patients.

Methods

Differentially expressed genes (DEGs) were screened out from The Cancer Genome Atlas (TCGA) database. The lncATLAS was applied to determine the differentially expressed lncRNAs (DElncRNAs) of the cytoplasm. The miRcode, miRDB, miRTarBase, and TargetScan databases were utilized to predict the interactions of DElncRNAs, DEmiRNAs, and DEmRNAs. Cytoscape was used to construct the ceRNA network. Then, a lncRNA prognostic model (LPM) was constructed based on ceRNA-related lncRNA that was significantly related to overall survival (OS), and its predictive ability was evaluated. Moreover, an LPM-based nomogram model was constructed. The significantly different expression of genes in the LPM was validated in an independent clinical cohort (N = 21) by quantitative RT-PCR.

Results

A novel ceRNA regulatory network, including 73 lncRNAs, 8 miRNAs, and 21 mRNAs was constructed. Functional enrichment analysis indicated that integral components of membrane and PI3K-Akt signaling pathway represented the most significant GO terms and pathway, respectively. The LPM established based on three lncRNAs (MIAT, LINC00460, and LINC00443) of great prognostic value from the ceRNA network was proven to be independent of conventional clinical parameters to differentiate patients with low or high risk of poor survival, with the AUC of 1-, 5- and 10-year OS were 0.723, 0.714 and 0.826 respectively. Furthermore, the nomogram showed a better predictive value in KIRC patients than individual prognostic parameters. The expression of MIAT and LINC00460 was significantly upregulated in the KIRC samples, while the expression of LINC00443 was significantly downregulated compared with the adjacent normal samples in the clinical cohort, TCGA, and GTEx.

Conclusion

This LPM based on three-lncRNA could serve as an independent prognostic factor with a tremendous predictive ability for KIRC patients, and the identified novel ceRNA network may provide insight into the prognostic biomarkers and therapeutic targets of KIRC.

Significance of PI3K signalling pathway in clear cell renal cell carcinoma in relation to VHL and HIF status

Renal cell carcinoma (RCC) includes diverse tumour types characterised by various genetic abnormalities. The genetic changes, like mutations, deletions and epigenetic alterations, play a crucial role in the modification of signalling networks, tumour pathogenesis and prognosis. The most prevalent RCC type, clear cell RCC (ccRCC), is asymptomatic in the early stages and has a poorer prognosis compared with the papillary and the chromophobe types RCCs. Generally, ccRCC is refractory to chemotherapy and radiation therapy. Loss of von Hippel-Lindau (VHL) gene and upregulation of hypoxia-inducible factors (HIF), the signature of most sporadic ccRCC, promote multiple growth factors. Hence, VHL/HIF and a variety of pathways, including phosphatase and TEnsin homolog on chromosome 10/phosphatidylinositol-3-kinase (PI3K)/AKT, are closely connected and contribute to the ontogeny of ccRCC. In the recent decade, multiple targeting agents have been developed based on blocking major signalling pathways directly or indirectly involved in ccRCC tumour progression, metastasis, angiogenesis and survival. However, most of these drugs have limitations; either metastatic ccRCC develops resistance to these agents, or despite blocking receptors, tumour cells use alternate signalling pathways. This review compiles the state of knowledge about the PI3K/AKT signalling pathway confined to ccRCC and its cross-talks with VHL/HIF pathway.

Construction and validation of a seven-gene signature for predicting overall survival in patients with kidney renal clear cell carcinoma via an integrated bioinformatics analysis

Kidney renal clear cell carcinoma (KIRC) remains a significant challenge worldwide because of its poor prognosis and high mortality rate, and accurate prognostic gene signatures are urgently required for individual therapy. This study aimed to construct and validate a seven-gene signature for predicting overall survival (OS) in patients with KIRC. The mRNA expression profile and clinical data of patients with KIRC were obtained from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC). Prognosis-associated genes were identified, and a prognostic gene signature was constructed. Then, the prognostic efficiency of the gene signature was assessed. The results obtained using data from the TCGA were validated using those from the ICGC and other online databases. Gene set enrichment analyses (GSEA) were performed to explore potential molecular mechanisms. A seven-gene signature (PODXL, SLC16A12, ZIC2, ATP2B3, KRT75, C20orf141, and CHGA) was constructed, and it was found to be effective in classifying KIRC patients into high- and low-risk groups, with significantly different survival based on the TCGA and ICGC validation data set. Cox regression analysis revealed that the seven-gene signature had an independent prognostic value. Then, we established a nomogram, including the seven-gene signature, which had a significant clinical net benefit. Interestingly, the seven-gene signature had a good performance in distinguishing KIRC from normal tissues. GSEA revealed that several oncological signatures and GO terms were enriched. This study developed a novel seven-gene signature and nomogram for predicting the OS of patients with KIRC, which may be helpful for clinicians in establishing individualized treatments.

Identification of a three-m6A related gene risk score model as a potential prognostic biomarker in clear cell renal cell carcinoma

Background

Clear cell renal cell carcinoma (ccRCC) is one of the most prevalent malignancies worldwide, N6-methyladenosine (m6A) has been shown to play important roles in regulating gene expression and phenotypes in both health and disease. Here, our purpose is to construct a m6A-regulrator-based risk score (RS) for prediction of the prognosis of ccRCC.

Methods

We used clinical and expression data of m6A related genes from The Cancer Genome Atlas (TCGA) dataset and the Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression analysis to develop an RS to predict survival of patients with ccRCC, and analyzed correlations between RS and other clinical indicators such as age, grade and stage. Validation of this RS was then engaged in another cohort, E-MTAB-1980 from the ArrayExpress dataset. Finally, we used quantitative real-time PCR to analyze the expression profile of genes consists of the RS.

Results

A three-gene RS including METTL3, METTL14 and HNRNPA2B1 which can predict overall survival (OS) of ccRCC patients from TCGA. After applying this RS into the validation cohort from Arrayexpress, we found that it successfully reproduced the result; furthermore, the results of PCR validation were in line with our analysis.

Conclusion

To sum up, our study has identified an RS composed of m6A related genes that may predict the prognosis of ccRCC patients, which might be helpful for future therapeutic strategies. Our results call for further experimental studies for validations.