Comprehensive analysis of lncRNA-mediated ceRNA network in renal cell carcinoma based on GEO database

Abstract

Renal cell carcinoma (RCC) ranks among the leading causes of cancer-related mortality. Despite extensive research, the precise etiology and progression of RCC remain incompletely elucidated. Long noncoding RNA (lncRNA) has been identified as competitive endogenous RNA (ceRNA) capable of binding to microRNA (miRNA) sites, thereby modulating the expression of messenger RNAs (mRNA) and target genes. This regulatory network is known to exert a pivotal influence on cancer initiation and progression. However, the specific role and functional significance of the lncRNA-miRNA-mRNA ceRNA network in RCC remain poorly understood. The RCC transcriptome data was obtained from the gene expression omnibus database. The identification of differentially expressed long noncoding RNAs (DElncRNAs), differentially expressed miRNAs, and differentially expressed mRNAs (DEmRNAs) between RCC and corresponding paracancer tissues was performed using the “Limma” package in R 4.3.1 software. We employed a weighted gene co-expression network analysis to identify the key DElncRNAs that are most relevant to RCC. Subsequently, we utilized the encyclopedia of RNA interactomes database to predict the interactions between these DElncRNAs and miRNAs, and the miRDB database to predict the interactions between miRNAs and mRNAs. Therefore, key DElncRNAs were obtained to verify the expression of their related genes in the The Cancer Genome Atlas database and to analyze the prognosis. The construction of RCC-specific lncRNA-miRNA-mRNA ceRNA network was carried out using Cytoscape 3.7.0. A total of 286 DElncRNAs, 56 differentially expressed miRNAs, and 2065 DEmRNAs were identified in RCC. Seven key DElncRNAs (GAS6 antisense RNA 1, myocardial infarction associated transcript, long intergenic nonprotein coding RNA 921, MMP25 antisense RNA 1, Chromosome 22 Open Reading Frame 34, MIR34A host gene, MIR4435-2 host gene) were identified using weighted gene co-expression network analysis and encyclopedia of RNA interactomes databases. Subsequently, a network diagram comprising 217 nodes and 463 edges was constructed based on these key DElncRNAs. The functional analysis of DEmRNAs in the ceRNA network was conducted using Kyoto Encyclopedia of Genes and Genomes and gene ontology. We constructed RCC-specific ceRNA networks and identified the crucial lncRNAs associated with RCC using bioinformatics analysis, which will help us further understand the pathogenesis of this disease.