Dataset Information

Identification of prognostic genes in kidney renal clear cell carcinoma by RNA?seq data analysis.

ABSTRACT: The present study aimed to analyze RNA-seq data of kidney renal clear cell carcinoma (KIRC) to identify prognostic genes. RNA?seq data were downloaded from The Cancer Genome Atlas. Feature genes with a coefficient of variation (CV) >0.5 were selected using the genefilter package in R. Gene co?expression networks were constructed with the WGCNA package. Cox regression analysis was performed using the survive package. Furthermore, a functional enrichment analysis was conducted using Database for Annotation, Visualization and Integrated Discovery tools. A total of 533 KIRC samples were collected, from which 6,758 feature genes with a CV >0.5 were obtained for further analysis. The KIRC samples were divided into two sets: The training set (n=319 samples) and the validation set (n=214 samples). Subsequently, gene co?expression networks were constructed for the two sets. A total of 12 modules were identified, and the green module was significantly associated with survival time. Genes from the green module were revealed to be implicated in the cell cycle and p53 signaling pathway. In addition, a total of 11 hub genes were revealed, and 10 of them (CCNA2, CDC20, CDCA8, GTSE1, KIF23, KIF2C, KIF4A, MELK, TOP2A and TPX2) were validated as possessing prognostic value, as determined by conducting a survival analysis on another gene expression dataset. In conclusion, a total of 10 prognostic genes were identified in KIRC. These findings may help to advance the understanding of this disease, and may also provide potential biomarkers for therapeutic development.

SUBMITTER: Gu Y

PROVIDER: S-EPMC5364979 | biostudies-literature | 2017 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Identification of prognostic genes in kidney renal clear cell carcinoma by RNA‑seq data analysis.

Gu Yanqin Y Lu Linfeng L Wu Lingfeng L Chen Hao H Zhu Wei W He Yi Y

Molecular medicine reports 20170213 4

The present study aimed to analyze RNA-seq data of kidney renal clear cell carcinoma (KIRC) to identify prognostic genes. RNA‑seq data were downloaded from The Cancer Genome Atlas. Feature genes with a coefficient of variation (CV) >0.5 were selected using the genefilter package in R. Gene co‑expression networks were constructed with the WGCNA package. Cox regression analysis was performed using the survive package. Furthermore, a functional enrichment analysis was conducted using Database for A ...[more]

PMID: 28260099

Similar Datasets

Project description:BackgroundLong 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.MethodsDifferentially 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.ResultsA 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.ConclusionThis 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.

Project description:Kidney Renal Clear Cell Carcinoma (KIRC) is one of fatal genitourinary diseases and accounts for most malignant kidney tumours. KIRC has been shown resistance to radiotherapy and chemotherapy. Like many types of cancers, there is no curative treatment for metastatic KIRC. Using advanced sequencing technologies, The Cancer Genome Atlas (TCGA) project of NIH/NCI-NHGRI has produced large-scale sequencing data, which provide unprecedented opportunities to reveal new molecular mechanisms of cancer. We combined differentially expressed genes, pathways and network analyses to gain new insights into the underlying molecular mechanisms of the disease development.Followed by the experimental design for obtaining significant genes and pathways, comprehensive analysis of 537 KIRC patients' sequencing data provided by TCGA was performed. Differentially expressed genes were obtained from the RNA-Seq data. Pathway and network analyses were performed. We identified 186 differentially expressed genes with significant p-value and large fold changes (P < 0.01, |log(FC)| > 5). The study not only confirmed a number of identified differentially expressed genes in literature reports, but also provided new findings. We performed hierarchical clustering analysis utilizing the whole genome-wide gene expressions and differentially expressed genes that were identified in this study. We revealed distinct groups of differentially expressed genes that can aid to the identification of subtypes of the cancer. The hierarchical clustering analysis based on gene expression profile and differentially expressed genes suggested four subtypes of the cancer. We found enriched distinct Gene Ontology (GO) terms associated with these groups of genes. Based on these findings, we built a support vector machine based supervised-learning classifier to predict unknown samples, and the classifier achieved high accuracy and robust classification results. In addition, we identified a number of pathways (P < 0.04) that were significantly influenced by the disease. We found that some of the identified pathways have been implicated in cancers from literatures, while others have not been reported in the cancer before. The network analysis leads to the identification of significantly disrupted pathways and associated genes involved in the disease development. Furthermore, this study can provide a viable alternative in identifying effective drug targets.Our study identified a set of differentially expressed genes and pathways in kidney renal clear cell carcinoma, and represents a comprehensive computational approach to analysis large-scale next-generation sequencing data. The pathway and network analyses suggested that information from distinctly expressed genes can be utilized in the identification of aberrant upstream regulators. Identification of distinctly expressed genes and altered pathways are important in effective biomarker identification for early cancer diagnosis and treatment planning. Combining differentially expressed genes with pathway and network analyses using intelligent computational approaches provide an unprecedented opportunity to identify upstream disease causal genes and effective drug targets.

Project description:Background: Kidney renal clear cell carcinoma (KIRC) is the malignancy originated from the renal epithelium, with a high rate of distant metastasis. Aberrant alternative splicing (AS) of pre-mRNA are widely reported to be involved in the tumorigenesis and metastasis of multiple cancers. The aim of this study is to explore the mechanism of alternative splicing events (ASEs) underlying tumorigenesis and metastasis of KIRC. Methods: RNA-seq of 537 KIRC samples downloaded from the TCGA database and ASEs data from the TCGASpliceSeq database were used to identify ASEs in patients with KIRC. The univariate and Lasso regression analysis were used to screen the most significant overall survival-related ASEs (OS-SEs). Based on those, the OS-SEs model was proposed. The interaction network of OS-SEs and splicing factors (SFs) with absolute value of correlation coefficient value >0.750 was constructed by Pearson correlation analysis. The OS-SEs significantly related to distant metastasis and clinical stage were identified by non-parametric test, and those were also integrated into co-expression analysis with prognosis-related Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways identified by Gene Set Variation Analysis (GSVA). ASEs with significance were selected for multiple online database validation. Results: A total of prognostic 6,081 overall survival-related ASEs (OS-SEs) were identified by univariate Cox regression analysis and a prediction model was constructed based on 5 OS-SEs screened by Lasso regression with the Area Under Curve of 0.788. Its risk score was also illustrated to be an independent predictor, which the good reliability of the model. Among 390 identified candidate SFs, DExD-Box Helicase 39B (DDX39B) was significantly correlated with OS and metastasis. After external database validation, Retained Intron of Ras Homolog Family Member T2 (RHOT2) and T-Cell Immune Regulator 1 (TCIRG1) were identified. In the co-expression analysis, overlapped co-expression signal pathways for RHOT2 and TCIRG1 were sphingolipid metabolism and N-glycan biosynthesis. Conclusions: Based on the results of comprehensive bioinformatic analysis, we proposed that aberrant DDX39B regulated RHOT2-32938-RI and TCIRG1-17288-RI might be associated with the tumorigenesis, metastasis, and poor prognosis of KIRC via sphingolipid metabolism or N-glycan biosynthesis pathway.

Dataset Information

Identification of prognostic genes in kidney renal clear cell carcinoma by RNA?seq data analysis.

Publications

Identification of prognostic genes in kidney renal clear cell carcinoma by RNA‑seq data analysis.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets