Open Access

Detection of dysregulated competing endogenous RNA modules associated with clear cell kidney carcinoma

  • Authors:
    • Hong Wang
    • Dahua Xu
    • Huiying Huang
    • Ying Cui
    • Chunhua Li
    • Chunrui Zhang
    • Shengnan Guo
    • Lining Zhang
    • Xiaomu Xu
    • Jiankai Xu
    • Jianping Lu
    • Liqiang Wang
    • Kongning Li
  • View Affiliations

  • Published online on: June 19, 2018     https://doi.org/10.3892/mmr.2018.9189
  • Pages: 1963-1972
  • Copyright: © Wang et al. This is an open access article distributed under the terms of Creative Commons Attribution License.

Metrics: Total Views: 0 (Spandidos Publications: | PMC Statistics: )
Total PDF Downloads: 0 (Spandidos Publications: | PMC Statistics: )


Abstract

Recent evidence has suggested that competitive endogenous RNAs (ceRNAs) are important regulatory molecules in clear cell kidney carcinoma (KIRC) and their dysregulation may contribute to cancer pathogenesis. However, the critical roles of dysregulated ceRNAs in KIRC remain unknown. In the present study, a KIRC dysregulated ceRNA‑ceRNA network (KDCCNet) was constructed based on the ʻceRNA hypothesisʼ by integrating microRNA regulation and expression profiles in cancerous and normal tissues. Two dysregulated patterns of ceRNAs interaction (gain and loss) exist in KDCCNet. The two modules, which are 95% loss interactions and 97% gain interactions, were demonstrated to be able to distinguish normal samples from cancer samples. Two long non‑coding (lnc)‑RNAs (glucuronidase β pseudogene 11 and LIFR antisense RNA 1) demonstrated significant associations with KIRC prognosis. The present study of the KDCCNet revealed a novel biological mechanism for KIRC and provides novel lncRNAs as candidate prognostic biomarkers.

Introduction

Human genome sequence data indicates that >98% of human genome transcripts are non-coding RNAs (ncRNAs), while protein-coding genes account for a small proportion of RNAs (1). Advances in the sensitivity of experimental assays have revealed that ncRNAs are involved in a variety of key biological processes (2). ncRNAs are divided into two groups based on length: Long non-coding RNAs (lncRNAs) and small non-coding RNAs. LncRNAs are a class of pervasively transcribed RNA molecules with lengths of >200 nucleotides. A previous study indicated that lncRNAs serve a critical role in regulating cellular processes (3). MicroRNAs (miRNAs) are small single-stranded RNAs of 18–25 nucleotides in length and are among the most important types of small ncRNAs. miRNAs can repress their targets by binding to miRNA response elements (MREs) in the 3′ untranslated region of mRNAs, causing translational repression or mRNA degradation. Notably, each miRNA regulates numerous RNA targets, including mRNA but also lncRNAs; the majority of RNA molecules harbor several MREs and are repressed by different miRNAs. As a result of this targeted multiplicity, lncRNAs and mRNA may compete for limited pools of miRNAs, therefore acting as competitive endogenous RNAs (ceRNAs) (4).

A number of studies have suggested that lncRNAs serve a critical role in regulating cancer processes through a ceRNA-associated mechanism (57). Various studies have investigated the ceRNA mechanism in specific types of cancer by constructing ceRNA networks (810). However, ceRNA regulation is required for the functioning of normal and cancerous cells. Therefore, detecting dysregulated ceRNA interactions by comparing ceRNA networks in normal and cancer cells would increase the understanding of the initiation and propagation of cancer.

Clear cell kidney carcinoma (KIRC) is the most common type of renal cell carcinoma, accounting for 70–80% of all renal cell carcinoma cases (11). The incidence and fatality rate of KIRC have increased over the past two decades (12). Although a number of studies have examined KIRC, few biomarkers have been identified (1315). Therefore, further studies of the key molecular mechanisms involved in KIRC and biomarkers for the early detection of patients with KIRC are urgently required. The dysregulated ceRNA network may be useful for investigating the pathogenesis and identifying novel biomarkers of KIRC.

In the present study, to establish an activated dysregulated ceRNA-ceRNA network, the regulatory association of lncRNAs, miRNAs and mRNAs were integrated based on the ‘ceRNA hypothesis’. Next, the expression patterns of lncRNAs, miRNAs and mRNAs were investigated in KIRC and normal tissues. By comparing ceRNA interactions in KIRC and normal tissues, an activated KIRC dysregulated ceRNA-ceRNA network (KDCCNet) was established. Topological feature analysis of KDCCNet demonstrated that lncRNAs are the key components of the KDCCNet. In the KDCCNet, there are two types of dysregulated patterns of ceRNAs interactions, these being gain and loss interactions. A total of two typical modules were discovered, which could distinguish between normal and cancer samples. The results of the present study provide insight into the underlying mechanisms of dysregulated ceRNA of KIRC and a global view of KDCCNet.

Materials and methods

Data resources

KIRC expression profiles of mRNAs, miRNAs and lncRNAs were obtained from The Cancer Genome Atlas (TCGA, http://cancergenome.nih.gov/) database release 1.0 including 517 cancer and 71 adjacent normal samples. All of the expression profiles were log2 transformed.

An independent dataset was downloaded from the GEO database (www.ncbi.nlm.nih.gov/geo/; accession number GSE36895), which is supported by the Affymetrix Human genome U133 plus 2.0 array (Affymetrix; Thermo Fisher Scientific, Inc., Waltham, MA, USA) (16). This dataset contains 29 tumor samples and 23 adjacent normal samples.

In order to investigate human miRNA-target regulatory associations, miRNA target information were collected from 7 databases, which includes miRanda 2010 (17), miRcode 11 (18), miRTarBase 6.0 (19), starBase V2.0 (20), TargetScan 7.1 (21), miRecords 2013 (22) and lnCeDB 2014 (23). Following redundancy analysis and standardization using R version 3.3.3, the databases were processed and 1,698,634 miRNA-mRNA interactions and 613,523 miRNA-lncRNA interactions were obtained. These interactions contained 18,013 mRNAs, 12,154 lncRNAs and 2,586 miRNAs. In addition, 93 KIRC-associated miRNAs were collected from the Human microRNA Disease Database V2.0 (24).

Construction of the ceRNA network

In order to identify ceRNA interactions, the miRNA target interactions were considered as well as the expression associations of the miRNA-mRNA, miRNA-lncRNA and lncRNA-mRNA. In this case two principles, explained below, were followed to identify ceRNA interactions.

Firstly, a hypergeometric test was used to compute the significance of shared miRNAs for each possible lncRNA-mRNA pair. The P-value was calculated according to:

P=1-∑i=or(Ki)(N-KM-i)(NM)

Where N represents the total number of human miRNAs; K represents the total number of miRNAs targeting mRNAs; M represents the total number of miRNAs targeting lncRNAs; and r represents the number of common miRNAs between lncRNA-mRNA pairs. All P-values were subjected to false discovery rate (FDR) correction and pairs with FDR<0.01 were considered as candidate ceRNA interaction pairs.

The Pearson correlation coefficient (R) was computed for each candidate ceRNA pair mentioned above, to identify the active ceRNA pairs in KIRC and normal tissues. Candidate ceRNA pairs with R>0 and P-adjusted<0.01 were identified as ceRNA interactions.

Following assembling all identified ceRNA pairs, the ceRNA network for KIRC and normal tissue was separately generated.

Identification of dysregulated ceRNA interactions

To identify dysregulated ceRNA interactions, the altered R was estimated for the expression of the ceRNA pair in the KIRC samples compared with normal samples. The altered R was estimated as follows:

ΔR=CorrKIRC(A,B)-Corrnormal(A,B)

CorrKIRC (A, B) is the R of a ceRNA pair in the KIRC samples and Corrnormal (A, B) is the R of certain ceRNA pairs in the normal samples.

Based on the assumption that the ΔR of a certain ceRNA pair is significant, the cancer/normal labels 1,000 times were randomly permuted. The significant P-value of each ΔR was given as follows:

P=∑i=1NSiN

When Si =1, it represents that the value of random ΔR for a certain ceRNA pair was greater than the real one, otherwise Si=0.

Identification of KIRC-associated dysregulated ceRNA modules

The edges of the dysregulated ceRNA network were weighted by ΔR. Based on value of ΔR, dysregulated ceRNAs interactions were divided into two patterns, gain ceRNAs interactions (ΔR>0) and loss ceRNAs interactions (ΔR<0). Then the Markov Cluster Algorithm clustering computational method was applied to identify the dysregulated ceRNA modules in GraphWeb (25) (biit.cs.ut.ee/graphweb/; Inflate parameter=1.8, count of nodes >10).

The classification performance of dysregulated modules through unsupervised hierarchical cluster analysis were visualized and 1- R was used as a distance measure between two samples.

Differential expression analysis

Differential expression analysis was performed to identify the genes associated with KIRC. This was achieved by comparing the RNA expression between KIRC samples and normal samples. Differential expression was detected by fold change (FC) and using the limma method (FDR<0.01, logFC>0.5).

Functional enrichment analysis

The gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis was performed using an online analysis tool (26) (Gene set enrichment analysis; software.broadinstitute.org/gsea/index.jsp; FDR<0.05).

Survival analysis

The univariate Cox regression was performed in order to evaluate the association between survival and the expression level of each lncRNA. Then a risk score was obtained by combining the Cox regression coefficients and lncRNA expression for each of the patients. The median risk score was used as the cut-off point and the patients were divided into high-risk and low-risk groups. The differences in the survival times were analyzed through the use of a log rank test.

Statistical analysis

Student's t-test was used to compare the difference of degree distribution between mRNA and lncRNA in KDCCNet (27). And the differences of the expression level, the transcript length and the number of exons of dysregulated and non-dysregulated lncRNAs were also compared based on the t-test. P<0.05 was considered to indicate a statically significant difference. The statistical analyses were performed using R version 3.3.3.

Results

Construction and global properties of KDCCNet

Regulatory interactions of ceRNAs are present in normal and cancer cells. In the present study, a two-stage analysis method was introduced for identifying ceRNA regulatory interactions in normal and KIRC tissues. As described in the Materials and methods, 585 and 3,646 ceRNA regulatory interactions were identified in normal tissue and KIRC tissue, respectively. Out of the ceRNA regulatory interactions in normal tissue, there were 67 lncRNAs and 415 mRNAs, while there were 127 lncRNAs and 1,706 mRNAs ceRNA regulatory interactions in KIRC tissue.

KIRC transformation from normal cells occurs depending upon the dynamics of multiple molecular and physiological alterations, including regulatory interactions of dysregulated ceRNAs (28). To evaluate the roles of dysregulated ceRNA in KIRC, a KDCCNet was constructed using a computational method based on the expression data of 517 KIRC patients and 71 adjacent normal pairs from TCGA. As described above, the KDCCNet was obtained (Fig. 1A). The KDCCNet contained 862 mRNAs, 123 lncRNAs and 1,195 edges. A total of two dysregulated patterns of ceRNAs interaction (gain and loss) exist in KDCCNet. In KDCCNet, 695 gains and 500 losses were detected (Table I). Therefore, certain normal ceRNA regulations may disappear or be reduced, while novel ceRNA regulations may be generated. These results suggest that dysregulated ceRNAs may serve an important role in KIRC progression.

Table I.

Statistics of the nodes and edges in the competitive endogenous RNAs dysregulated network.

Table I.

Statistics of the nodes and edges in the competitive endogenous RNAs dysregulated network.

TypemRNAlncRNAEdges
Dysregulated network8621231,195
Gain520  75695
Loss367  61500

[i] Lnc, long non-coding.

By topological feature analysis, the topological features of the KDCCNet were revealed. Examination of the degree of distribution of KDCCNet revealed the power law distribution (R=0.792). In the KDCCNet, the majority of ceRNAs had few interacting dysregulated ceRNA partners, while a small subset of ceRNAs regarded as hubs possessed a relatively large number of interacting dysregulated ceRNA partners (Fig. 1B). Furthermore, it was demonstrated that the characteristic path length of KDCCNet was significantly increased when compared with random linked networks (P<0.001; Fig. 1C).

LncRNAs as topological key components involved in KIRC progression

In the KDCCNet, there are two types of dysregulated ceRNAs: LncRNAs and mRNAs. To analyze the roles of dysregulated ceRNA in more detail, the contribution degree of dysregulated ceRNAs in the two groups was investigated. It was demonstrated that the degree distribution of lncRNAs was significantly increased compared with the mRNA in KDCCNet (t-test, P=3.33×10−13; Fig. 1B). This result indicated that lncRNAs covered the majority of hubs and lncRNAs were topological key nodes in the KDCCNet associated with KIRC progression. In addition, a recent study demonstrated that lncRNAs were key regulators of the ceRNA network (29).

As dysregulated lncRNAs serve more important roles in KIRC progression, it was hypothesized that there are significant differences between dysregulated lncRNAs and non-dysregulated lncRNAs. To further analyze dysregulated lncRNAs, the characteristics of dysregulated lncRNAs and non-dysregulated lncRNAs were compared. As dysregulated lncRNAs were expected to contribute to the KIRC aberrant transcriptome, the expression levels of dysregulated lncRNAs and non-dysregulated lncRNAs were investigated first. The results demonstrated that the expression level of dysregulated lncRNAs was significantly increased compared with non-dysregulated lncRNAs (P=1.74×10−115, t-test; Fig. 1D). Secondly, the distribution of lengths for dysregulated lncRNAs and non-dysregulated lncRNAs on chromosomes was analyzed. It was revealed that the length of dysregulated lncRNAs on chromosomes was significantly increased compared with non-dysregulated lncRNAs (P=1.97×10−29, t-test; Fig. 1D). Finally, the counts of exons of dysregulated lncRNAs and non-dysregulated lncRNAs were compared, which also demonstrated that the counts of exons of dysregulated lncRNAs was significantly increased compared with non-dysregulated lncRNAs (P=1.45×10−12, t-test; Fig. 1D).

These results demonstrated that lncRNAs are topological key nodes of the KDCCNet associated with KIRC progression. Dysregulated lncRNAs exhibited increased expression in KIRC tissues, longer chromosome lengths and larger numbers of exons. This makes dysregulated lncRNAs more suitable as molecular sponges of miRNAs in KIRC.

Module analysis of the KDCCNet

It is generally thought that cancer occurs through the interactive activity of a number of genes. Functionally coherent modules can be regarded as biomarkers for distinguishing cancer from normal tissue (30). Therefore, the Markov Cluster Algorithm clustering computational method was applied to identify dysregulated ceRNA modules in the KDCCNet. As a result, 18 dysregulated ceRNA modules were obtained (Inflate=1.8, count of ceRNAs >10). The modules were enriched with known KIRC miRNAs in the Human microRNA Disease Database. An unsupervised hierarchical cluster analysis was conducted to visualize the classification performance of the modules. Modules that could distinguish cancer from normal tissue in TCGA and validation dataset were selected (area under the curve>0.65). Finally, five dysregulated ceRNA modules were obtained.

The distribution of the two patterns of ceRNA interactions was analyzed and two typical modules (95% loss interactions and 97% gain interactions) were observed (Fig. 2A). They were defined as loss ceRNA interactions module (LC module) and gain ceRNA interactions module (GC module), respectively (Table II).

Table II.

A detailed description of the LC and GC modules.

Table II.

A detailed description of the LC and GC modules.

ModuleNodesEdgesLostGainΔRmRNAlncRNAmiRNA
LC module37390.950.050.7334343  7
GC module30290.030.970.740129112

[i] Lnc, long non-coding; GC, competitive endogenous RNA gains; LC, competitive endogenous RNA losses; ΔR, change in the Pearson correlation coefficient; miRNA, microRNA.

Typical dysregulated ceRNA modules contribute to KIRC

Firstly, unsupervised hierarchical cluster analysis demonstrated that the two typical modules as aforementioned could distinguish the normal samples from the cancer samples (Fig. 2B and C). The differentially expressed ceRNAs covered 67 and 50% in the LC and GC modules, respectively (Fig. 2A). Therefore, the two typical modules contributed to KIRC progression. In addition, the dysregulated ceRNAs of the two typical modules and their miRNA regulators were distributed on different chromosomes (Fig. 3).

Secondly, the GO and KEGG functional enrichment analysis of the LC and GC modules was investigated. Through functional enrichment analysis, it was demonstrated that the LC module was significantly enriched in vitamin D metabolic processes (FDR q-value=2.70×10−4; data not shown). Vitamin D metabolic processes have been demonstrated to inhibit the cell cycle and induce apoptosis in various types of cancer, including KIRC (31). This suggested that the LC module may contribute to KIRC progression by influencing vitamin D metabolic processes via loss dysregulated ceRNA interactions. Additionally, functional enrichment analysis revealed that the GC module was significantly enriched in the renal cell carcinoma pathway (FDR q-value=6.23×10−3; data not shown), which is associated with KIRC.

ceRNA contributes to cancer progression by influencing the expression of mRNAs in cancer associated GO biological processes or the KEGG pathway (32). It was hypothesized that if the two ceRNAs within the pair were differentially expressed simultaneously and at least one ceRNA was in KIRC-associated pathways, then they may exhibit more important effects on the dysregulated ceRNA module than the other ceRNAs which have no significant difference in expression or are independent of KIRC. It was demonstrated that two loss ceRNA interactions [glucuronidase β pseudogene 11 (GUSBP11) vs. cytochrome P450 family 27 subfamily B member 1 (CYP27B1), and RP11-752D24.2 vs. CYP27B1; Fig. 4A) composed of three differentially expressed ceRNAs within the LC module, influenced the vitamin D metabolic processes via CYP27B1 and one gain ceRNA interaction [LIFR antisense RNA 1 (LIFR-AS1) vs. vascular endothelial growth factor A (VEGFA); Fig. 4B] composed of two differentially expressed ceRNAs within the GC module, influenced the renal cell carcinoma pathway via VEGFA. Blomberg et al (33) demonstrated that differential expressed CYP27B1 was associated with KIRC, however the mechanism is unknown. The results of the present study further complemented the suggestion that CYP27B1 may be upregulated by ceRNAs (GUSBP11 and RP11-752D24.2) in KIRC. The above indicated that the three-dysregulated ceRNA interaction pairs may serve a key role via affecting the KIRC-associated pathways.

As lncRNAs are the key components for the typical modules, attention was paid to lncRNAs of the three ceRNA pairs mentioned above. Based on survival analysis, two lncRNAs were demonstrated to be associated with the survival of patients with KIRC (GUSBP11, P=1.78×10−13; LIFR-AS1, P=1.32×10−3). In fact, GUSBP11 has been discovered to be associated with the survival of patients with head and neck squamous cell carcinoma (34). Also, LIFR-AS1 was predicted to mediate the inhibitory action of leukemia inhibitory factor, a cytokine that is involved in embryonic uterine development and is associated with the tumor size of the uterine fibroid (35). The results of the present study established the link between KIRC, and GUSBP11 and LIFR-AS1 via dysregulated ceRNA analysis.

In addition, miRNAs mediate competitive interactions, therefore expression pattern alterations of miRNAs and ceRNAs may cause dysregulation of ceRNA interactions. The correlations between miRNAs and ceRNAs in GC and LC module are shown in Fig. 5 and miRNA regulators of the dysregulation modules were evaluated based on the Human microRNA Disease Database (HMDD). It was found that 4 miRNAs regulators of the LC module (miRNA-141, miRNA-203, miRNA-9 and miRNA-218) and 2 miRNA (miRNA-138 and miRNA-146b) regulators of the GC module were involved in KIRC from HMDD. For example, miR-203 inhibits renal cancer cell proliferation, migration and invasion (36). miR-203 was also a regulator of the LC module, and therefore it was hypothesized that this dysregulated ceRNA interaction may be involved in the mechanism of KIRC.

Discussion

In the present study, a multidimensional integration strategy was introduced to construct the KDCCNet. Topological analysis identified that the dysregulated lncRNAs were the key components of the KDCCNet. A total of two dysregulated patterns of ceRNAs interaction exist in the KDCCNet: i) Gain ceRNAs interactions: A lncRNA-mRNA pair demonstrated no competitive interactions in normal samples but did in the KIRC samples; ii) Loss ceRNAs interactions: A lncRNA-mRNA pair demonstrated competing activities in normal samples, though in not cancer samples. Gain of KIRC ceRNA interactions and also loss of normal ceRNA interactions occurred during the initiation and progression of KIRC. The two patterns contributed to KIRC progression. Next, two KIRC-associated dysregulated ceRNA modules were identified, the LC and GC modules, which can distinguish cancer from normal samples in TCGA and the validation dataset. Furthermore, two dysregulated lncRNAs (GUSBP11 and LIFR-AS1) in the two modules demonstrated associations with the survival of KIRC patients. In addition, certain KIRC associated genes/miRNAs were observed in the modules. For example, overproduction of VEGFA in the GC module occurs early during the development of VHL−/− renal cell carcinomas [the patients who have biallelic inactivating mutations of the von Hippel-Lindau (VHL) gene] and therefore reduces the selection pressure to activate ‘collateral’ angiogenic pathways (37,38). miRNA-203 in the LC module could inhibit renal cancer cell proliferation, migration and invasion (36). The present study of the KDCCNet revealed a novel biological mechanism for KIRC and provided a novel link between the two dysregulated lncRNAs (GUSBP11 and LIFR-AS1) and KIRC.

Although genome-wide analysis was performed to investigate the potential mechanism of dysregulated ceRNA interactions in KIRC progression and identified KIRC-associated dysregulated ceRNA modules, there were certain limitations to the present study. The understanding of the functional roles of lncRNAs is limited and further studies should examine the remaining regulatory elements, including transcription factors, DNA methylation and copy number variation, which may reveal a more detailed molecular mechanism of KIRC progression.

In conclusion, a KIRC-associated dysregulated ceRNA network was developed and a method for examining the competitive interactions of RNAs that are potentially dysregulated in KIRC progression. These results increase the understanding of lncRNA characteristics and provide a foundation for leveraging publicly available genomic data to study the functions and mechanisms of dysregulated ceRNAs in KIRC.

Acknowledgements

Not applicable.

Funding

The present study was supported in part by the Natural Science Foundation of Heilongjiang Province (grant no. B201302), National Natural Science Foundation of China (grant no. 31301094), Education Department Foundation of Heilongjiang Province (grant no. 12531227) and Health Department Foundation of Heilongjiang Province (grant no. 2012–798).

Availability of data and materials

The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.

Authors' contributions

KL and LW conceived and designed the experiments. HW, DX, HH, YC, CL, CZ, SG, LZ, XX, JX and JL collected and analyzed the data. KL and LW wrote the manuscript. All authors read and approved the final manuscript.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare they have no competing interests.

Glossary

Abbreviations

Abbreviations:

KIRC

clear cell kidney carcinoma

lncRNAs

long non-coding RNAs

ceRNA

competing endogenous RNAs

KDCCNet

KIRC dysregulated ceRNA-ceRNA network

References

1 

Ponting CP and Belgard Grant T: Transcribed dark matter: Meaning or myth? Hum Mol Genet. 19:R162–R168. 2010. View Article : Google Scholar : PubMed/NCBI

2 

Parker BJ, Moltke I, Roth A, Washietl S, Wen J, Kellis M, Breaker R and Pedersen JS: New families of human regulatory RNA structures identified by comparative analysis of vertebrate genomes. Genome Res. 21:1929–1943. 2011. View Article : Google Scholar : PubMed/NCBI

3 

Huarte M: The emerging role of lncRNAs in cancer. Nat Med. 21:1253–1261. 2015. View Article : Google Scholar : PubMed/NCBI

4 

Salmena L, Poliseno L, Tay Y, Kats L and Pandolfi PP: A ceRNA hypothesis: The rosetta stone of a hidden RNA language? Cell. 146:353–358. 2011. View Article : Google Scholar : PubMed/NCBI

5 

Liu XH, Sun M, Nie FQ, Ge YB, Zhang EB, Yin DD, Kong R, Xia R, Lu KH, Li JH, et al: Lnc RNA HOTAIR functions as a competing endogenous RNA to regulate HER2 expression by sponging miR-331-3p in gastric cancer. Mol Cancer. 13:922014. View Article : Google Scholar : PubMed/NCBI

6 

Xia T, Liao Q, Jiang X, Shao Y, Xiao B, Xi Y and Guo J: Long noncoding RNA associated-competing endogenous RNAs in gastric cancer. Sci Rep. 4:60882014. View Article : Google Scholar : PubMed/NCBI

7 

Tay Y, Rinn J and Pandolfi PP: The multilayered complexity of ceRNA crosstalk and competition. Nature. 505:344–352. 2014. View Article : Google Scholar : PubMed/NCBI

8 

Zhou X, Liu J and Wang W: Construction and investigation of breast-cancer-specific ceRNA network based on the mRNA and miRNA expression data. IET Syst Biol. 8:96–103. 2014. View Article : Google Scholar : PubMed/NCBI

9 

Xu J, Li Y, Lu J, Pan T, Ding N, Wang Z, Shao T, Zhang J, Wang L and Li X: The mRNA related ceRNA-ceRNA landscape and significance across 20 major cancer types. Nucleic Acids Res. 43:8169–8182. 2015. View Article : Google Scholar : PubMed/NCBI

10 

Poliseno L and Pandolfi PP: PTEN ceRNA networks in human cancer. Methods. 77–78:41–50. 2015. View Article : Google Scholar

11 

Zhang L, Xul B, Chen S, Lu K, Liu C, Wang Y, Zhao Y, Zhang X, Liu D and Chen M: The complex roles of microRNAs in the metastasis of renal cell carcinoma. J Nanosci Nanotechnol. 13:3195–3203. 2013. View Article : Google Scholar : PubMed/NCBI

12 

Siegel RL, Miller KD and Jemal A: Cancer statistics, 2017. CA Cancer J Clin. 67:7–30. 2017. View Article : Google Scholar : PubMed/NCBI

13 

Chen G, Wang Y, Wang L and Xu W: Identifying prognostic biomarkers based on aberrant DNA methylation in kidney renal clear cell carcinoma. Oncotarget. 8:5268–5280. 2017.PubMed/NCBI

14 

Zhan Y, Guo W, Zhang Y, Wang Q, Xu XJ and Zhu L: A five-gene signature predicts prognosis in patients with kidney renal clear cell carcinoma. Comput Math Methods Med. 2015:8427842015. View Article : Google Scholar : PubMed/NCBI

15 

Linehan WM: Genetic basis of kidney cancer: Role of genomics for the development of disease-based therapeutics. Genome Res. 22:2089–2100. 2012. View Article : Google Scholar : PubMed/NCBI

16 

Peña-Llopis S, Vega-Rubín-de-Celis S, Liao A, Leng N, Pavía-Jiménez A, Wang S, Yamasaki T, Zhrebker L, Sivanand S, Spence P, et al: BAP1 loss defines a new class of renal cell carcinoma. Nat Genet. 44:751–759. 2012. View Article : Google Scholar : PubMed/NCBI

17 

Enright AJ, John B, Gaul U, Tuschl T, Sander C and Marks DS: MicroRNA targets in Drosophila. Genome Biol. 5:R12003. View Article : Google Scholar : PubMed/NCBI

18 

Jeggari A, Marks DS and Larsson E: miRcode: A map of putative microrna target sites in the long non-coding transcriptome. Bioinformatics. 28:2062–2063. 2012. View Article : Google Scholar : PubMed/NCBI

19 

Hsu SD, Lin FM, Wu WY, Liang C, Huang WC, Chan WL, Tsai WT, Chen GZ, Lee CJ, Chiu CM, et al: MiRTarBase: A database curates experimentally validated microRNA-target interactions. Nucleic Acids Res. 39:(Database Issue). D163–D169. 2011. View Article : Google Scholar : PubMed/NCBI

20 

Li JH, Liu S, Zhou H, Qu LH and Yang JH: StarBase v2.0: Decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 42:(Database Issue). D92–D97. 2014. View Article : Google Scholar : PubMed/NCBI

21 

Agarwal V, Bell GW, Nam JW and Bartel DP: Predicting effective microRNA target sites in mammalian mRNAs. Elife. 4:2015. View Article : Google Scholar

22 

Xiao F, Zuo Z, Cai G, Kang S, Gao X and Li T: miRecords: An integrated resource for microRNA-target interactions. Nucleic Acids Res. 37:(Database Issue). D105–D110. 2009. View Article : Google Scholar : PubMed/NCBI

23 

Das S, Ghosal S, Sen R and Chakrabarti J: lnCeDB: Database of human long noncoding RNA acting as competing endogenous RNA. PLoS One. 9:e989652014. View Article : Google Scholar : PubMed/NCBI

24 

Li Y, Qiu C, Tu J, Geng B, Yang J, Jiang T and Cui Q: HMDD v2.0: A database for experimentally supported human microRNA and disease associations. Nucleic Acids Res. 42:(Database Issue). D1070–D1074. 2014. View Article : Google Scholar : PubMed/NCBI

25 

Reimand J, Tooming L, Peterson H, Adler P and Vilo J: GraphWeb: Mining heterogeneous biological networks for gene modules with functional significance. Nucleic Acids Res. 36:(Web Server Issue). W452–W459. 2008. View Article : Google Scholar : PubMed/NCBI

26 

Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES and Mesirov JP: Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 102:15545–15550. 2005. View Article : Google Scholar : PubMed/NCBI

27 

Statistician William Sealy Gosset (Student): The probable error of a mean. Biometrika. 6:1–25. 1908. View Article : Google Scholar

28 

Xiao H, Tang K, Liu P, Chen K, Hu J, Zeng J, Xiao W, Yu G, Yao W, Zhou H, et al: LncRNA MALAT1 functions as a competing endogenous RNA to regulate ZEB2 expression by sponging miR-200s in clear cell kidney carcinoma. Oncotarget. 6:38005–38015. 2015. View Article : Google Scholar : PubMed/NCBI

29 

Wang LK, Chen XF, He DD, Li Y and Fu J: Dissection of functional lncRNAs in Alzheimer's disease by construction and analysis of lncRNA-mRNA networks based on competitive endogenous RNAs. Biochem Biophys Res Commun. 485:569–576. 2017. View Article : Google Scholar : PubMed/NCBI

30 

Zhang Y, Xuan J, Clarke R and Ressom HW: Module-based breast cancer classification. Int J Data Min Bioinform. 7:284–302. 2013. View Article : Google Scholar : PubMed/NCBI

31 

Lee J and Park SH: Tumor-suppressive activity of 1,25-dihydroxyvitamin D3 against kidney cancer cells via up-regulation of FOXO3. Biosci Biotechnol Biochem. 80:1947–1953. 2016. View Article : Google Scholar : PubMed/NCBI

32 

Zhou M, Diao Z, Yue X, Chen Y, Zhao H, Cheng L and Sun J: Construction and analysis of dysregulated lncRNA-associated ceRNA network identified novel lncRNA biomarkers for early diagnosis of human pancreatic cancer. Oncotarget. 7:56383–56394. 2016.PubMed/NCBI

33 

Jensen Blomberg M, Andersen CB, Nielsen JE, Bagi P, Jørgensen A, Juul A and Leffers H: Expression of the vitamin D receptor, 25-hydroxylases, 1alpha-hydroxylase and 24-hydroxylase in the human kidney and renal clear cell cancer. J Steroid Biochem Mol Biol. 121:376–382. 2010. View Article : Google Scholar : PubMed/NCBI

34 

Cao W, Liu JN, Liu Z, Wang X, Han ZG, Ji T, Chen WT and Zou X: A three-lncRNA signature derived from the Atlas of ncRNA in cancer (TANRIC) database predicts the survival of patients with head and neck squamous cell carcinoma. Oral Oncol. 65:94–101. 2017. View Article : Google Scholar : PubMed/NCBI

35 

Aissani B, Zhang K and Wiener H: Genetic determinants of uterine fibroid size in the multiethnic NIEHS uterine fibroid study. Int J Mol Epidemiol Genet. 6:9–19. 2015.PubMed/NCBI

36 

Xu M, Gu M, Zhang K, Zhou J, Wang Z and Da J: miR-203 inhibition of renal cancer cell proliferation, migration and invasion by targeting of FGF2. Diagn Pathol. 10:242015. View Article : Google Scholar : PubMed/NCBI

37 

Kaelin WG Jr: The von Hippel-Lindau tumor suppressor protein and clear cell renal carcinoma. Clin Cancer Res. 13:S680–S684. 2007. View Article : Google Scholar

38 

Kim WY and Kaelin WG: Role of VHL gene mutation in human cancer. J Clin Oncol. 22:4991–5004. 2004. View Article : Google Scholar : PubMed/NCBI

Related Articles

Journal Cover

August-2018
Volume 18 Issue 2

Print ISSN: 1791-2997
Online ISSN:1791-3004

Sign up for eToc alerts

Recommend to Library

Copy and paste a formatted citation
x
Spandidos Publications style
Wang H, Xu D, Huang H, Cui Y, Li C, Zhang C, Guo S, Zhang L, Xu X, Xu J, Xu J, et al: Detection of dysregulated competing endogenous RNA modules associated with clear cell kidney carcinoma. Mol Med Rep 18: 1963-1972, 2018
APA
Wang, H., Xu, D., Huang, H., Cui, Y., Li, C., Zhang, C. ... Li, K. (2018). Detection of dysregulated competing endogenous RNA modules associated with clear cell kidney carcinoma. Molecular Medicine Reports, 18, 1963-1972. https://doi.org/10.3892/mmr.2018.9189
MLA
Wang, H., Xu, D., Huang, H., Cui, Y., Li, C., Zhang, C., Guo, S., Zhang, L., Xu, X., Xu, J., Lu, J., Wang, L., Li, K."Detection of dysregulated competing endogenous RNA modules associated with clear cell kidney carcinoma". Molecular Medicine Reports 18.2 (2018): 1963-1972.
Chicago
Wang, H., Xu, D., Huang, H., Cui, Y., Li, C., Zhang, C., Guo, S., Zhang, L., Xu, X., Xu, J., Lu, J., Wang, L., Li, K."Detection of dysregulated competing endogenous RNA modules associated with clear cell kidney carcinoma". Molecular Medicine Reports 18, no. 2 (2018): 1963-1972. https://doi.org/10.3892/mmr.2018.9189