The ceRNA Network Has Potential Prognostic Value in Clear Cell Renal Cell Carcinoma: A Study Based on TCGA Database

Liu, Haosheng; Zhu, Zhaowen; Fang, Jianxiong; Liu, Tianqi; Zhang, Zhenhui; Zhao, Chao; Pu, Xiaoyong; Liu, Jiumin

doi:https://doi.org/10.1155/2020/4830847

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 4830847 | https://doi.org/10.1155/2020/4830847

The ceRNA Network Has Potential Prognostic Value in Clear Cell Renal Cell Carcinoma: A Study Based on TCGA Database

Haosheng Liu,¹Zhaowen Zhu,²Jianxiong Fang,³Tianqi Liu,⁴Zhenhui Zhang,³Chao Zhao,¹Xiaoyong Pu,⁵and Jiumin Liu⁵

Academic Editor: Kui Li

Received29 Nov 2019

Accepted25 Feb 2020

Published26 Jun 2020

Abstract

Clear cell renal cell carcinoma (ccRCC) is a very common cancer in urology. Many evidences suggest that complex changed pathways take a nonnegligible part in the occurrence and development of ccRCC. Nevertheless, the underlying mechanism is not clear. In this study, the expression data between ccRCC and normal tissue samples in TCGA database were compared to distinguish differentially expressed genes (DEGs: mRNAs, miRNAs, and lncRNAs). Afterwards, we used GO enrichment and KEGG pathway enrichment analyses to explore the functions of these DEGs. We also found the correlation between three RNAs and created a competing endogenous RNA (ceRNA) network. Moreover, we used univariate Cox regression analysis to select DEGs that are connected with overall survival (OS) of ccRCC patients. We found 1652 mRNAs, 1534 lncRNAs, and 173 miRNAs that were distinguished in ccRCC compared with normal tissues. According to GO analysis, the maladjusted mRNAs are mainly concentrated in immune cell activation and kidney development, while according to KEGG, they are mainly concentrated in pathways related to cancer. A total of 5 mRNAs, 1 miRNA, and 4 lncRNAs were connected with patients’ OS. In this article, a network of lncRNA-miRNA-mRNA was established; it is expected to be able to indicate possible molecular mechanisms for initial of ccRCC and provide a new viewpoint for diagnosis of ccRCC.

1. Introduction

RCC accounts for about 3% of human cancers; in the last 20 years, morbidity has increased about 2% per year both worldwide and in Europe, causing about 99200 RCC patients and 39100 RCC-related deaths in EU in 2018 [1]. RCC is the commonest kidney essence lesion, accounting for approximately 90% of all renal malignances. Many patients with renal mass have no symptoms until the terminal stage. At present, >60% of RCC patients are discovered by ultrasound (US) or CT performed for other disease. Flank pain, visible haematuria, and palpable abdominal mass which are classical signs of RCC are less now [2]. Therefore, an in-depth study of the biomarkers of tumorigenesis and progression in RCC is badly needed to gain the solutions and realize diagnosis goals.

In the human genome, protein-coding genes account for less than 3 percent; however, more than 80% of our genomes have no protein-coding capacity. Such transcription is called noncoding RNAs (ncRNAs), which are divided into long noncoding RNAs (lncRNAs) and small noncoding RNAs (sncRNAs). More and more evidences showed that lncRNAs are connected with tumorigenesis and development through a variety of biological processes, including transcriptional regulation, self-renewal, and carcinogenesis [3]. The rapid development of RNA-Seq technology allows people to discover new lncRNAs related to urologic malignancies [4]. The article analyzed the differentially expressed lncRNAs and their function in ccRCC.

microRNAs are sncRNAs with regulatory function [5]. After binding to the 3-UTR of mRNA, it can inhibit the process of transcription and even induce mRNA degradation. However, miRNAs are also influenced by other ncRNAs, for example, lncRNAs and circular RNAs (circRNAs), which serve as sponges for miRNAs [6, 7]. Currently, dysregulated miRNAs are believed to be connected with the occurrence of tumors including gastrointestinal cancer [8], colon cancer [9], prostate cancer [10], and colorectal cancer [11]. Thus, miRNAs may be a good biomarker for cancer diagnosis.

Competitive endogenous RNA theory indicated that RNAs which have miRNA binding sites struggle for posttranscriptional control. This theory has got wide interest as a whole function of lncRNAs and circRNAs and a substitution function of mRNAs. Moreover, the hypothesis assumes that specific RNAs can regulate the expression of miRNA target gene by inhibiting miRNA activity [12]. However, there are few studies on ceRNA of ccRCC.

We explored the interaction among three differentially expressed RNAs from The Cancer Genome Atlas (TCGA) database. In addition, these RNAs were annotated for possible biological functions. Then, we established a ceRNA network in ccRCC. Besides, we used K-M and receiver operating characteristic analysis to distinguish prognostic DEGs for predicting OS of ccRCC. Our results provide new insights into diagnosis strategies for ccRCC.

2. Materials and Methods

2.1. Data Download and Processing

We downloaded ccRCC patients’ transcriptome profiling and clinical data from TCGA database (https://tcga-data.nci.nih.gov/tcga/), which was inputted on the Illumina HiSeq RNA-Seq platform. The exclusion criteria were (i) samples without clinical data and (ii) samples without complete information of the stage and overall survival period. Finally, a total of 534 cases contain 472 clear cell renal cell carcinoma tissues and 62 adjacent nontumor renal tissues were included in our study. We annotated the names of RNAs through Homo sapiens.GRCH38.84.chr.gtf.gz downloaded from Ensembl Genome Browser 97 (http://www.ensembl.org). Both RNA profile data and clinical characteristics of ccRCC are publicly available and in open-access platforms. Therefore, approval by local ethics committee was not needed.

2.2. Analysis of Differentially Expressed Genes

We used the edgeR package to make comparisons between tumor tissues and normal specimens to select the DEGs. DEmRNAs and DElncRNAs were distinguished through the threshold of a and an adj. value < 0.05, and the DEmiRNAs was identified using the threshold of a and an adj. value < 0.05.

2.3. Cluster Analysis of the DEGs

R software was used for cluster analysis based on the expression value of each sample. With R package “pheatmap,” hierarchical clustering analysis was conducted. We showed the results through a clustergram. The column means the samples, and the row shows the gene expression level.

2.4. Functional Enrichment Analysis

We analyzed the DEGs’ functional enrichment with the online bioinformatics tools KOBAS 3.0 (http://kobas.cbi.pku.edu.cn/index.php) and DAVID (https://david.ncifcrf.gov/, version 6.8). KEGG pathways and GO terms were considered to be statistically significant when value < 0.05.

2.5. Predict the Target of miRNAs

We used miRDB (http://www.mirdb.org/), TargetScan (http://www.targetscan.org/), and miRTarBase (http://amp.pharm.mssm.edu/Harmonizome/resource/MiRTarBase) to predict the possible target genes of DEmiRNAs. The overlapped target genes were searched by Venn overlap analysis. miRNAs’ target lncRNAs were predicted through the miRcode dataset (http://www.mircode.org/).

2.6. Construction of the ceRNA Network

We identified DEGs and the relationships between miRNA-mRNA and miRNA-lncRNA. Based upon these results, Cytoscape 3.6.1 (version 3.6.1, San Diego CA) was used to establish the ceRNA network.

2.7. Construction of the Protein-Protein Interaction (PPI) Network

To explore the interaction between mRNAs involved in this ceRNA network, a PPI network was constructed by the Search Tool for the Retrieval of Interacting Genes (http://string.embl.de/) with a composite score greater than 0.4 as the cut-off criterion [13].

2.8. Survival and Receiver Operating Characteristic Analyses

We used a univariate Cox regression model to analyze the connection between ccRCC patients’ OS and ceRNA network DEGs. value < 0.05 was regarded as statistically significant. Then, those DEGs were presented with the survival curve. The cut-off point divided all samples into high and low expression groups. In addition, sensitivity and specificity were calculated through ROC curves and AUC values.

2.9. Statistical Analysis

We analyzed the data through R and Cytoscape software. Data were expressed as . The differences between two groups were analyzed through fold change and Student’s -test. The survival curves were constructed by K-M method, and log-rank tests were used to test survival differences. When , we thought the difference are significant.

3. Results

3.1. Clustering Analysis of DEGs

The number of DEGs was 1652 (mRNA), 173 (miRNA), and 1534 (lncRNA), respectively. Hierarchical clustering analysis was used to analyze all the DEGs. (Figure 1). Besides, we showed the top 10 upregulated and downregulated DEGs in ccRCC (Tables 1–3).

(a)

(b)

(c)

3.2. GO Enrichment Analysis of DEmRNAs

We used DAVID database to analyze these dysregulated genes and perform enrichment analysis to determine their respective functions. The upregulated genes were mainly enriched in T cell activation and regulation of leukocyte activation (Figures 2(a) and 2(c)). The downregulated genes were mostly enriched in renal system development and urogenital system development (Figures 2(b) and 2(d)).

(a)

(b)

(c)

(d)

3.3. KEGG Pathway Enrichment Analysis of DEmRNAs

The KEGG pathway enrichment method was also used to analyze the DEmRNAs. These DEmRNAs had the strongest correlation with neuroactive ligand-receptor interaction, cytokine-cytokine receptor interaction, metabolism of xenobiotics by cytochrome P450, cell adhesion molecules (CAMs), and complement and coagulation cascades (Figure 3). Based on these results, we believed that the DEmRNAs play an indispensable part in the progression of ccRCC.

(a)

(b)

3.4. ceRNA Network

There were 44 DElncRNAs that may bind to 25 DEmiRNAs. We searched for the DEmiRNAs’ target mRNAs through using miRDB, miRTarBase, and TargetScan. Based on the Venn intersection analysis, 33 DEmRNAs were identified. At last, we established a ceRNA network (Figure 4).

3.5. Protein-Protein Interaction Network

To clarify the relationships between the 33 DEmRNAs involved in this ceRNA network, a PPI network was constructed using the STRING database containing 32 nodes and 15 edges (Figure 5). KEGG pathway enrichment analysis of these mRNAs was performed with KOBAS 3.0 with as the cut-off criterion. The results showed that these mRNAs were mainly enriched in pathways related to cancer (Table 4).

3.6. Analysis of the Survival-Associated DEGs

It is important to screen DEGs that could predict the prognosis of ccRCC. We used a univariate Cox proportional hazards regression model to analyze ccRCC patients’ OS. There were 4 lncRNAs, 1 miRNA, and 5 mRNAs screened out: all of the lncRNAs, miR-144, and NETO2 positively influenced ccRCC patients’ OS, whereas OS was negatively connected with NOD2, PAPPA, PCDH, and SPI2 (Figure 6). Compared with normal tissues, AC011383.1, PSORS1C3, miR-144, NETO2, NOD2, and SPI2 levels were increased in ccRCC tissues while ALDH1L1-AS2, DNAJC3-AS1, PAPPA, and PCDH9 levels were decreased () (Figure 7). Receiver operating characteristic analysis was performed on these DEGs; their AUC value ranged from 0.736 to 0.976 (Figure 8, Table 5).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

4. Discussion

The abnormal expression of protein-coding and protein-noncoding transcription is one feature of the cancer transcriptome [14]. The prognosis of ccRCC has been a big challenge due to the difficultly of diagnosis. So, it is significant to find effective biomarkers of ccRCC. In these years, dysregulated mRNAs [15], miRNAs [16], and lncRNAs [17] have been reported in ccRCC. There have been many reports of these RNAs on kidney tumor alone, but the interaction among mRNA, miRNA, and lncRNA is still unclear. Many of the genes with known ceRNA interactors have been found to be associated with liver cancer, leukaemias, lymphomas, and so on [18]. It is necessary to study the ceRNA network in ccRCC systematically.

In this study, 4 prognostic DElncRNAs (AC011383.1, ALDH1L1-AS2, DNAJC3-AS1, and PSORS1C3) were included in the ceRNA network. They may be independent prognostic factors of ccRCC patients’ OS. Liang et al. reported that DNAJC3-AS1 plays a positive role in osteosarcoma evolvement by regulating DNAJC3. Moreover, it is a possible marker and treatment point for osteosarcoma [19].

Besides, we identified one independent prognostic DEmiRNA (miR-144). The upregulated miR-144 could enhance ccRCC malignancy and resistance to sunitinib by regulating ARID1A [20]. Nevertheless, Liu et al. believed that miR-144 inhibited cancer cell proliferation and metastasis in renal cell carcinoma [21]. The results in this study show that miR-144 has a positive impaction on ccRCC patients’ OS.

We found 5 mRNAs that are connected with ccRCC patients’ prognosis: NETO2, NOD2, PAPPA, PCDH9, and SPI1. Oparina et al. report that NETO2’s expression level was increased in 90% ccRCC and 50% of papillary renal cancers. It is a possible biomarker in kidney cancer [22]. Mey et al. found that NOD2’s expression level is higher in tumor tissues compared to normal tissues [23]. Dalgin et al. use microarray gene expression profiling to identify specific renal cell carcinoma markers. They identified 158 genes that dysregulated in tumor tissues; these genes are related to proteolysis and cell adhesion, including PAPPA [24].

However, there is no research about PCDH9 and SPI1 in ccRCC. In the present work, we found the two novel mRNAs as key predictors of ccRCC prognosis. PCDH9 is a member of the protocadherin family. It is associated with a kind of tumors, for example, melanoma [25], ovarian cancer [26], and medulloblastoma [27]. Ren et al. proved that PCDH9 is a tumor-inhibiting gene and has prognostic value in prostate cancer [28]. SPI1 is a member of the ETS family. It is relevant to multiple malignancies, for example, papillary thyroid carcinoma ³⁰, pediatric T cell acute lymphoblastic leukemia [29], and breast cancer [30]. SPI1 is closely related to clinical features such as grade, metastasis, and stage, which means that SPI1 may be a potential prognostic biomarker of ccRCC.

A PPI network was built to illustrate the relationship between the mRNAs involved in the ceRNA network. Some of these mRNAs are closely connected with each other, suggesting that their interactions may play an important role in the development of ccRCC, which further confirms the significant role of this ceRNA network in ccRCC.

However, this study has some limitations. First, although we established the ceRNA network, we have not demonstrated other regulatory models. Besides, other studies may draw a different conclusion due to the differences in inclusion and exclusion criteria. Finally, there are few experimental data explaining the mechanisms of ceRNA, and we need further experiments to illustrate the role of ceRNA in ccRCC.

5. Conclusion

This article illustrated that lncRNAs, miRNAs, and mRNAs involved in the ceRNA network may be possible biomarkers. They are expected to predict the survival rate in ccRCC patients. Nevertheless, we need more experiments to validate these RNAs’ biological function.

Data Availability

Both RNA profiles data and clinical characteristics of ccRCC are publicly available and in open access platforms.

Conflicts of Interest

The authors declare that there is no conflict of interests.

Authors’ Contributions

Haosheng Liu and Zhaowen Zhu contributed equally to this work and are co-first authors.

Acknowledgments

The work described in this paper was fully supported by the National Natural Science Foundation of China (81570691), Guangdong Natural Science Foundation (2019A1515012019, 2017A030313715), Science and Technology Program of Guangdong Province (2017ZC0275).

References

J. Ferlay, M. Colombet, I. Soerjomataram et al., “Cancer incidence and mortality patterns in Europe: estimates for 40 countries and 25 major cancers in 2018,” European Journal of Cancer, vol. 103, pp. 356–387, 2018.
View at: Publisher Site | Google Scholar
A. Thorstenson, M. Bergman, A. H. Scherman-Plogell et al., “Tumour characteristics and surgical treatment of renal cell carcinoma in Sweden 2005–2010: a population-based study from the national Swedish kidney cancer register,” Scandinavian Journal of Urology, vol. 48, no. 3, pp. 231–238, 2014.
View at: Publisher Site | Google Scholar
T. Hung, Y. Wang, M. F. Lin et al., “Extensive and coordinated transcription of noncoding RNAs within cell-cycle promoters,” Nature Genetics, vol. 43, no. 7, pp. 621–629, 2011.
View at: Publisher Site | Google Scholar
E. S. Martens-Uzunova, R. Böttcher, C. M. Croce, G. Jenster, T. Visakorpi, and G. A. Calin, “Long noncoding RNA in prostate, bladder, and kidney cancer,” European Urology, vol. 65, no. 6, pp. 1140–1151, 2014.
View at: Publisher Site | Google Scholar
D. P. Bartel, “Metazoan microRNAs,” Cell, vol. 173, no. 1, pp. 20–51, 2018.
View at: Publisher Site | Google Scholar
C. Braicu, A. A. Zimta, D. Gulei, A. Olariu, and I. Berindan-Neagoe, “Comprehensive analysis of circular RNAs in pathological states: biogenesis, cellular regulation, and therapeutic relevance,” Cellular and Molecular Life Sciences, vol. 76, no. 8, pp. 1559–1577, 2019.
View at: Publisher Site | Google Scholar
M. Seles, G. C. Hutterer, T. Kiesslich et al., “Current insights into long non-coding RNAs in renal cell carcinoma,” International Journal of Molecular Sciences, vol. 17, no. 4, p. 573, 2016.
View at: Publisher Site | Google Scholar
M. Konno, J. Koseki, A. Asai et al., “Distinct methylation levels of mature microRNAs in gastrointestinal cancers,” Nature Communications, vol. 10, no. 1, p. 3888, 2019.
View at: Publisher Site | Google Scholar
L. Min, S. Zhu, L. Chen et al., “Evaluation of circulating small extracellular vesicles derived miRNAs as biomarkers of early colon cancer: a comparison with plasma total miRNAs,” Journal Extracellular Vesicles, vol. 8, no. 1, p. 1643670, 2019.
View at: Publisher Site | Google Scholar
J. Jeon, E. Olkhov-Mitsel, H. Xie et al., “Temporal stability and prognostic biomarker potential of the prostate cancer urine miRNA transcriptome,” JNCI: Journal of the National Cancer Institute, vol. 112, no. 3, pp. 247–255, 2020.
View at: Publisher Site | Google Scholar
Y. Okugawa, Y. Toiyama, K. Hur et al., “Circulating miR-203 derived from metastatic tissues promotes myopenia in colorectal cancer patients,” Journal of Cachexia Sarcopenia and Muscle, vol. 10, no. 3, pp. 536–548, 2019.
View at: Publisher Site | Google Scholar
S. Anfossi, A. Babayan, K. Pantel, and G. A. Calin, “Clinical utility of circulating non-coding RNAs -- an update,” Nature Reviews. Clinical Oncology, vol. 15, no. 9, pp. 541–563, 2018.
View at: Publisher Site | Google Scholar
J. Wang, C. Zhang, W. He, and X. Gou, “Construction and comprehensive analysis of dysregulated long non-coding RNA-associated competing endogenous RNA network in clear cell renal cell carcinoma,” Journal of Cellular Biochemistry, vol. 120, no. 2, pp. 2576–2593, 2019.
View at: Publisher Site | Google Scholar
F. A. Karreth and P. P. Pandolfi, “ceRNA cross-talk in cancer: when ce-bling rivalries go awry,” Cancer Discovery, vol. 3, no. 10, pp. 1113–1121, 2013.
View at: Publisher Site | Google Scholar
M. R. Morris and F. Latif, “The epigenetic landscape of renal cancer,” Nature Reviews Nephrology, vol. 13, no. 1, pp. 47–60, 2017.
View at: Publisher Site | Google Scholar
E. A. Braga, M. V. Fridman, V. I. Loginov, A. A. Dmitriev, and S. G. Morozov, “Molecular mechanisms in clear cell renal cell carcinoma: role of miRNAs and hypermethylated miRNA genes in crucial oncogenic pathways and processes,” Frontiers in Genetics, vol. 10, 2019.
View at: Publisher Site | Google Scholar
X. Liu, Y. Hao, W. Yu et al., “Long non-coding RNA emergence during renal cell carcinoma tumorigenesis,” Cellular Physiology and Biochemistry, vol. 47, no. 2, pp. 735–746, 2018.
View at: Publisher Site | Google Scholar
Y. Tay, J. Rinn, and P. P. Pandolfi, “The multilayered complexity of ceRNA crosstalk and competition,” Nature, vol. 505, no. 7483, pp. 344–352, 2014.
View at: Publisher Site | Google Scholar
R. Liang, Z. Liu, Z. Chen et al., “Long noncoding RNA DNAJC3-AS1 promotes osteosarcoma progression via its sense-cognate gene DNAJC3,” Cancer Medicine, vol. 8, no. 2, pp. 761–772, 2019.
View at: Publisher Site | Google Scholar
W. Xiao, N. Lou, H. Ruan et al., “Mir-144-3p promotes cell proliferation, metastasis, sunitinib resistance in clear cell renal cell carcinoma by downregulating ARID1A,” Cellular Physiology and Biochemistry, vol. 43, no. 6, pp. 2420–2433, 2017.
View at: Publisher Site | Google Scholar
F. Liu, N. Chen, R. Xiao, W. Wang, and Z. Pan, “MiR-144-3p serves as a tumor suppressor for renal cell carcinoma and inhibits its invasion and metastasis by targeting MAP3K8,” Biochemical and Biophysical Research Communications, vol. 480, no. 1, pp. 87–93, 2016.
View at: Publisher Site | Google Scholar
N. Y. Oparina, A. F. Sadritdinova, A. V. Snezhkina et al., “Increase in NETO2 gene expression is a potential molecular genetic marker in renal and lung cancers,” Russian Journal of Genetics, vol. 48, no. 5, pp. 506–512, 2012.
View at: Publisher Site | Google Scholar
L. Mey, M. Jung, F. Roos et al., “NOD1 and NOD2 of the innate immune system is differently expressed in human clear cell renal cell carcinoma, corresponding healthy renal tissue, its vasculature and primary isolated renal tubular epithelial cells,” Journal of Cancer Research and Clinical Oncology, vol. 145, no. 6, pp. 1405–1416, 2019.
View at: Publisher Site | Google Scholar
G. S. Dalgin, D. T. Holloway, L. S. Liou, and C. DeLisi, “Identification and characterization of renal cell carcinoma gene markers,” Cancer Inform, vol. 3, pp. 65–92, 2007.
View at: Google Scholar
H. Z. Yang, J. Zhang, J. Zeng et al., “Resveratrol inhibits the proliferation of melanoma cells by modulating cell cycle,” International Journal of Food Sciences and Nutrition, vol. 71, no. 1, pp. 84–93, 2020.
View at: Publisher Site | Google Scholar
N. Izycka, K. Sterzynska, R. Januchowski, and E. Nowak-Markwitz, “Semaphorin 3A (SEMA3A), protocadherin 9 (PCdh9), and S100 calcium binding protein A3 (S100A3) as potential biomarkers of carcinogenesis and chemoresistance of different neoplasms, including ovarian cancer - review of literature,” Ginekologia Polska, vol. 90, no. 4, pp. 223–227, 2019.
View at: Publisher Site | Google Scholar
C. J. Robbins, M. J. Bou-Dargham, K. Sanchez, M. C. Rosen, and Q. A. Sang, “Decoding somatic driver gene mutations and affected signaling pathways in human medulloblastoma subgroups,” Journal of Cancer, vol. 9, no. 24, pp. 4596–4610, 2018.
View at: Publisher Site | Google Scholar
S. Ren, G. H. Wei, D. Liu et al., “Whole-genome and transcriptome sequencing of prostate cancer identify new genetic alterations driving disease progression,” European Urology, vol. 73, no. 3, pp. 322–339, 2018.
View at: Publisher Site | Google Scholar
H. Y. Wu, Y. Wei, and S. L. Pan, “Down-regulation and clinical significance of miR-7-2-3p in papillary thyroid carcinoma with multiple detecting methods,” IET Systems Biology, vol. 13, no. 5, pp. 225–233, 2019.
View at: Publisher Site | Google Scholar
S. Kimura, “Genetic and epigenetic landscape of pediatric T-cell acute lymphoblastic leukemia,” Rinshō Ketsueki, vol. 60, no. 5, pp. 459–467, 2019.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Haosheng Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

464

Downloads

1174

Citations