Abstract
Haibing Xiao1,2, Lin Bao2, Wen Xiao2, Hailong Ruan2, Zhengshuai Song2, Yan Qu1, Ke Chen2, Xiaoping Zhang2 and Hongmei Yang1
1Department of Pathogenic Biology, School of Basic Medicine, Tongji Medical College, Huazhong University of Science and Technology, Wuhan 430030, China
2Department of Urology, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan 430022, China
Correspondence to:
Hongmei Yang, email: [email protected]
Xiaoping Zhang, email: [email protected]
Keywords: Lucat1; EZH2; p57
Received: May 02, 2017 Accepted: July 29, 2017 Published: September 23, 2017
ABSTRACT
Background: Many long intergenic noncoding RNAs (lincRNAs) are encoded in the human genome. However, their biological functions, molecular mechanisms and prognostic values associated with clear cell renal cell carcinoma (ccRCC) have yet to be elucidated.
Methods: We screened the lncRNAs’ profile in ccRCC from The Cancer Genome Atlas (TCGA) database, and selected Lucat1 for further study. MTS, colony formation assay and transwell assay were performed to examine the effect of Lucat1 on proliferation and metastasis of ccRCC. The Chip and Rip assay was performed to verify that Lucat1 can bind to polycomb PRC2 complex and suppress p57 expression.
Results: In this study, we found that lncRNA Lucat1 expression was significantly up regulated in tumor tissues compared to matched adjacent non-tumor tissues. The Lucat1 expression level was also associated with grade, the clinical pathological stage and the survival time. Functional assays showed that Lucat1 can promote renal cancer cell proliferation in vitro and in vivo. Further analysis showed that Lucat1 can bind to polycomb PRC2 complex and suppress p57 expression.
Conclusions: Taken together, our results suggest that Lucat1, as a regulator of proliferation, may serve as a candidate prognostic biomarker and target for novel therapies in human ccRCC.
INTRODUCTION
The rate of renal cell carcinoma (RCC) has been rising throughout the world. Renal cell carcinoma (RCC) accounts for nearly 5% of adult malignancies with about 63,920 new cases and 13,860 deaths estimated in 2014 in the United States [1]. In China, about 66,800 new cases and 23,400 deaths were estimated for 2014 [2]. Clear cell renal cell carcinoma (ccRCC) is the most common subtype and represents approximately 75% of all renal tumors. Metastasis is common in ccRCC, and approximately one-third of ccRCC patients have metastasis at the time of diagnosis despite the wide use of ultrasound and computed tomography [3]. In addition, renal cancer patients respond poorly to radiation treatment and conventional chemotherapy [4]. Hence, a better understanding of the mechanisms involved in the pathogenesis of ccRCC and more effective therapeutic approaches are crucial [5].
Long non-coding RNA (lncRNA) is a heterogeneous class of transcription RNA with a minimum length of 200 bases and limited protein-coding potential [6, 7]. LncRNAs are involved in multilevel regulation of gene expression, including transcriptional regulation by recruiting chromatin-modifying complexes [8, 9] and post-transcriptional regulation by interacting with miRNAs, mRNAs, or proteins [10, 11].
Recently, numerous lncRNAs have been identified to have a direct role in recruiting PRC2. PRC2, a methyltransferase that is composed of enhancer of zeste homolog 2 (EZH2), suppressor of zeste 12 (SUZ12) and embryonic ectoderm development (EED), can catalyze the di- and trimethylation of lysine residue 27 of histone 3 (H3K27me3), thus modulating gene expression. Ahmad M. Khalil etal have found that nearly 20% of lincRNAs expressed in various cell types are bound by PRC2 [12]. These lncRNAs epigenetically regulate gene expression through binding to PRC2 in various biological processes, especially during cancer, such as HOTAIR, TUG1, MALAT1, PINT etc. [13–20].
Lucat1, also named SCAL1, is induced by cigarette smoke and elevated in lung cancer cell lines [21] and may contribute to cisplatin resistance in high grade serous ovarian cancer [22]. An abstract implies that Lucat1 may decrease the expression of p21 and p57 in human non-small cell lung cancer [23]. However, the biological functions of Lucat1 in the control of ccRCC tumorigenesis remain largely unknown.
In this study, the TCGA database was used to search for lncRNA gene expression profiles in ccRCC. We identified Lucat1 as a new candidate lncRNA that promotes the development of ccRCC. Our data indicated that Lucat1 has higher expression in renal cancer cell lines and renal cancer tissues. We also found that Lucat1 is capable of facilitating cell growth, migration and invasion through epigenetically suppressing p57 in renal cancer cell lines.
RESULTS
Lucat1 is upregulated in ccRCC tissues and cell lines and indicates a poor prognosis
The TCGA database was used to search for differentially expressed lncRNAs between ccRCC tissues and normal tissues. Based on ccRCC RNA-seq data, Lucat1 was found overexpressed in ccRCC tissues compared with normal (Figure 1A). To further support this conclusion, we examined the expression of Lucat1 expression in 45 renal cancer tissues and their corresponding noncancerous tissues from Union Hosptial and obtained the same result (Figure 1B). In comparison, we also found Lucat1 was ubiquitously expressed at higher levels in a panel of 5 human clear cell renal cell carcinoma lines than immortalized human proximal renal tubule epithelial cell line HK-2 (Figure 1C). Moreover, to assess the clinical significance of Lucat1, we evaluated the correlation between its level and clinic-pathological parameters. Results revealed that Lucat1 levels were remarkably correlated with grade, TNM stage and metastasis in ccRCC (Figure 1D, 1E and 1F) (Table 1). Nevertheless, Lucat1 levels were not associated with other clinical characteristics, including gender (p = 0.115), age (p = 0.108) and recurrence (p = 0.225). Additionally, multivariate cox regression analysis revealed that high Lucat1 expression, age, TNM stage, grade, metastasis are independent predictors of OS in ccRCC patients (Table 2). Kaplan-Meier analysis indicated that high Lucat1 expression was related to a poorer overall survival (log-rank test, P <0.001, Figure 1G) and disease free survival (log-rank test, P <0.001, Figure 1H). Taken together, these results confirmed that high Lucat1 expression was related to poor prognosis, and upregulated expression of Lucat1 might be crucial in ccRCC tumorigenesis and progression.
Figure 1: Expression of Lucat1 in ccRCC cell lines, tissues and its clinical parameters. (A) Lucat1 was found to be highly over-expressed in ccRCC tissues compared with normal tissues in the TCGA RNA-seq data (P < 0.001). (B) Relative expression of Lucat1 in 45 pairs of ccRCC tumor tissues and their corresponding adjacent non-cancerous tissues. (C) Real-time PCR analysis of Lucat1 expression in immortalized human renal tubule epithelial cell line HK-2 and indicated renal carcinoma cell lines. (D, E, F) High Lucat1 expression was signifcantly correlated with the TNM grade, stage and metastasis. (G, H) High Lucat1 expression contributed to a significant poorer oval-all survival and disease-free survival at the TCGA database (n=480).*, t-test p<0.05; **, t-test p<0.001; ##, ANOVA, p<0.001.
Table 1: The characteristic of Lucat1 in clear cell renal cell carcinoma
Characteristic |
| Total (n = 448) | lncRNA-Lucat1 | p Value | |
---|---|---|---|---|---|
Low (224) | High (224) | ||||
Gender | Male | 287 | 135 | 152 |
|
| Female | 161 | 89 | 72 | 0.115 |
Age | <=60 | 218 | 118 | 100 |
|
| >60 | 230 | 106 | 124 | 0.108 |
T | T1&T2 | 276 | 156 | 120 |
|
| T3&T4 | 172 | 68 | 104 | 0.001 |
N | N0 | 218 | 108 | 110 |
|
| N1 | 16 | 3 | 13 | 0.017a,c |
| Nx | 214 | 113 | 101 | 0.031 |
M | M0 | 377 | 199 | 178 |
|
| M1 | 71 | 25 | 46 | 0.009 |
Stage | 1,2 | 260 | 148 | 112 |
|
| 3,4 | 188 | 76 | 112 | 0.001 |
Grade | 1,2 | 198 | 111 | 87 |
|
| 3,4 | 245 | 109 | 136 | 0.015b |
| X | 5 | 4 | 1 | 0.015c |
Recurrence | No | 349 | 180 | 169 |
|
| Yes | 99 | 44 | 55 | 0.225 |
a: N0, N1 and Lucat1 expression, b: Grade1,2 Grade3,4 and Lucat1 expression, c: Fisher Two-tailed Fisher’s exact test.
Table 2: Univariate and multivariate analyses of clinicopathological factors for overall survival
Risk factors | Univariate analysis | Multivariate analysis | ||||
---|---|---|---|---|---|---|
HR | p-value | 95% CI | HR | p-value | 95% CI | |
Lucat1 expression | 1.155 | 0.000** | 1.082-1.233 | 1.124 | 0.002** | 1.043-1.211 |
Age | 1.695 | 0.002** | 1.220-2.354 | 1.488 | 0.018* | 1.070-2.070 |
Grade | 2.421 | 0.000** | 1.692-3.463 | 1.489 | 0.040* | 1.017-2.179 |
Stage | 4.308 | 0.000** | 3.046-6.095 | 2.466 | 0.000** | 1.624-3.744 |
M | 4.765 | 0.000** | 3.424-6.632 | 2.299 | 0.000** | 1.564-3.380 |
T | 3.597 | 0.000** | 2.588-5.001 |
|
|
|
N | 2.821 | 0.002** | 1.484-5.365 |
|
|
|
Recurrence | 2.407 | 0.000** | 1.745-3.319 |
|
|
|
Gender | 1.022 | 0.896 | 0.737-1.418 |
|
|
|
Knockdown of Lucat1 inhibited cell proliferation and metastasis in vitro
To explore the role of Lucat1 in renal cancer cells, we stably inhibited Lucat1 in two ccRCC cell lines ACHN and 786-O with lenti-viruses carrying shRNA for Lucat1 and a control nonspecific shRNA (LacZ) (Figure 2A). MTS (Figure 2B, 2C) assay and Colony formation assay (Figure 2D, 2E) showed that knockdown Lucat1 inhibited cell proliferation in ACHN, 786-O and OS-RC-2 cells. Further assay of transwell showed that knockdown Lucat1 suppressed renal cancer cell migration and invasion (Figure 2F, 2G, 2H and 2I). Moreover, knockdown Lucat1 expression resulted in an accumulation of ACHN and 786-O cells in the G1 phase of the cell cycle (Figure 2J, 2K). In ccRCC cell lines OS-RC-2, we get the similar result (Supplementary Figure 1). We also found that up-regulation of Lucat1 can promote ccRCC proliferation in vitro (Suppementary Figure S1).
Figure 2: Knockdown of Lucat1 inhibited cell proliferation and metastasis in vitro. (A) The effciency of Lucat1 silencing in short hairpin RNA-stably transduced renal cancer cell lines ACHN and 786-O. Relative gene expression was determined using the comparative delta-delta CT method (2-ΔΔCt). (B, C) MTS assays revealed cell growth curves of indicated cells. (D, E) Representative micrographs of crystal violet-stained cell colonies analyzed by clongenic formation and its relative count. (F, G, H, I) Migration and invasion assay for renal cancer cells. Representative photographs were taken at × 200 magnifcation; number of migrated cells was quantified in ten random images from each treatment group. (J, K) Flow cytometric determination of proportion of indicated cells in distinct cell cycle phases.
Lucat1 inhibits p57 expression via directly binding with PRC2 complex in ccRCC cells
Lucat1 is a long intergenic noncoding RNA, and it is estimated that up to 20% lincRNAs physically associate with PRC2 [12]. Thus, we hypothesized that Lucat1 might function in such a manner in ccRCC. As for PRC2 complex exerts its epigenetic functions in the nucleus, we first investigated the subcellular localization of Lucat1 in ccRCC and found that a considerable increase in Lucat1 expression in the nucleus versus the cytosol (Figure 3A, 3B, 3C), suggesting that Lucat1 may have a major regulatory function at the transcriptional level through PRC2. Furthermore, we performed RIP assay and the results showed that Lucat1 could directly bind with enhancer of zeste homolog 2 (EZH2, the catalytic subunit of the PRC2) and SUZ12, another member of the PRC2 complex in ACHN and 786-O cells (Figure 3D, 3E, 3F). We also found that another long non-coding RNA Hotair could bind with EZH2 and SUZ12 as previously reported. Previous studies showed that the expression of PRC2 target genes, such as the cell-cycle regulation genes CDKN1A (p21CIP1), CDKN1B (p27KIP1), CDKN1C (p57Kip2), CDKN2A (p16INK4a), CDKN2B (p15INK4b) were up-regulated at the mRNA and protein levels [14, 29]. We then graphed a correlation between Lucat1 and CDKN1A (p21CIP1), CDKN1B (p27KIP1), CDKN1C (p57Kip2), CDKN2A (p16INK4a), CDKN2B (p15INK4b) in TCGA data, and found only CDKN1C had a negative correlation with Lucat1 (Supplementray Figure 2). And then we performed a correlation between EZH2 and CDKN1A (p21CIP1), CDKN1B (p27KIP1), CDKN1C (p57Kip2), CDKN2A (p16INK4a), CDKN2B (p15INK4b) in TCGA data and obtained the same result (Supplementray Figure 2). So we mainly focused on Lucat1, EHZ2 and p57. To address whether Lucat1 is involved in transcriptional repression through the enrichment of PRC2 to p57 promoters, we conducted chromatin immunoprecipitation (ChIP) analysis by Lucat1 knockdown. ChIP assays demonstrated that knockdown of Lucat1 decreased the binding of EZH2 and H3K27me3 levels across the p57 promoters (Figure 3G, 3H). We next explored the role of Lucat1 on the expression of p57. The results disclosed that Lucat1 knockdown increased the expression of p57 at the mRNA and protein levels (Figure 3I, 3J).
Figure 3: Lucat is associated with PRC2 in ccRCC. (A, B and C) Lucat1 nuclear localization, as identified using qRT-PCR in fractionated ACHN and 786-O cells. After nuclear and cytosolic separation, RNA expression levels were measured by qRT-PCR. GAPDH was used as a cytosolic marker, and U6 was used as a nuclear marker. (D) ACHN and 768-O cells transfected with psi-Flag (control), or Flag-EZH2 and Flag-SUZ12 plasmids. (E, F) RIP experiments were performed with Flag antibody, and the coprecipitated RNA was subjected to qRT-PCR for Lucat1 and Hotair. The fold enrichment of Lucat1 and Hotair in RIPs is relative to its matching IgG control RIP. (G, H) ChIP-qPCR of H3K27me3 and EZH2 of the promoter region of the p57 locus after between ACHN and 786-O cell transfected with sh-Lacz or sh-Lucat1. Antibody enrichment was quantified relative to the amount of input DNA. Antibody directed against IgG was used as a negative control. (I) qRT-PCR detected the expression of p57 after transfection with sh-Lacz or sh-Lucat1. (J) Western blot assays detected the expression of p57 after transfection with sh-Lacz or sh-Lucat1.
Silencing of p57 is partly involved in the oncogenic function of Lucat1
We conducted rescue assays to identify whether p57 is mediated in the Lucat1 induced ccRCC cell growth. We stably inhibited p57 in two ccRCC cell lines ACHN and 786-O with lenti-viruses carrying shRNA for p57 and a control nonspecifc shRNA (LacZ) (Figure 4A, 4B). The effects of Lucat1 expression on endogenous p57 protein were monitored. It showed that sh-Lucat1 can promote the expression of p57, whereas sh-p57 can relieve the promotion of p57 by Lucat1 in renal cancer cell lines ACHN and 768-O (Figure 4C and 4D). Colony formation assay and MTS demonstrated that the knockdown of Lucat1 inhibited cell growth was partly reversed by sh-p57 treatment (Figure 4E, 4F, 4G, 4H and Supplementary Figure 3). Further study showed that knockdown of Lucat1 lead G0-1 cycle arrest was partly reversed by down-regulation of p57 (Figure 4I and 4J).
Figure 4: Silencing of p57 is partly involved in the oncogenic function of Lucat1. (A, B) The effciency of p57 silencing in short hairpin RNA-stably transduced renal cancer cell lines ACHN and 786-O by qRT-PCR and western blot. (C, D) Western blot assays detected the expression of p57 after transfection with sh-Lucat1 and/or sh-P57 in ACHN and 768-O. (E, F) Colony forming assays were performed to determine the cell viability in ACHN and 786-O cell lines. (G, H) MTS assays were performed to determine the cell viability in ACHN and 786-O cell lines. (I, J) Flow cytometric determination of proportion of indicated cells in distinct cell cycle phases.
These findings indicated that Lucat1 exerts oncogenic effects in ccRCC cells partly through repressing p57 expression.
Knockout of EZH2 is involved in the oncogenic function of Lucat1
CDKN1C (p57) is a direct target of EZH2 and is suppressed by epigenetic mechanisms in breast cancer, ovarian cancer, non-small-cell lung cancer etc. [30, 31]. But the mechanism in ccRCC remains unclear. We applied the efficient CRISPR/Cas genome editing system [32] targeting exon 2 of EZH2 to knock out this gene in 786-O cells [26]. Consistent with previous reports, EZH2-null clones exhibited a significant increase in p57 and HOXA9 (also a direct target of EZH2) level comparing to controls (Figure 5A, 5B). We then stably transfected sh-Lacz and sh-Lucat1 in 786-O EZH2 knockout cell lines and found that stably transfected sh-Lucat1 can’t clearly inhibit the expression of p57 (Figure 5C). Colony assay, MTS and cycle test also showed that sh-Lucat1 did not exist obvious effects on proliferation in 786-O EZH2 knockout cell lines (Figure 5D, 5E, 5F and Supplementary Figure 3). These findings indicated that Lucat1 exertion of oncogenic effects in ccRCC cells may be accomplished mainly through interaction with EZH2.
Figure 5: Lucat1 have its function through EZH2. (A) Schematic Model of knockout EZH2 with CRISPR/Cas9 system. (B) Western blot for EZH2, p57 and HOXA9 in 786-O EZH2 knockout cell lines and its control cell lines. (C) Western blot assays detected the expression of p57 after transfection with sh-Lucat1 and/or sh-P57 in 768-O and 786-O knockout cell lines. (D, E, F) Colony forming assays, MTS assays and flow cytometric were performed to determine the cell viability in 786-O and 786-O EZH2 knockout cell lines. (G, H, I) Relative expression of EZH2, P57 in 45 pairs of ccRCC tumor tissues and their corresponding adjacent non-cancerous tissues. (J, K, L) The correlation of Lucat1, EZH2 and P57 in renal cancer tissue.
Next, we studied the clinic-pathological relevance of the Lucat1, its relative molecule p57, and EZH2 expression in 45 matched normal and ccRCC clinical tissue samples. We found significantly elevated expression of EZH2, Lucat1 and a lower expression of p57 in the tumor samples (Figure 5G, 5H and 5I). Bivariate correlation analysis showed that expression of Lucat1 and EZH2 were significantly correlated with p57 transcript level of ccRCC tissues (Figure 5J, 5K), but there was low correlation between Lucat1 and EZH2 (Figure 5L).
Knockdown of Lucat1 represses tumor growth in vivo
To further determine whether Lucat1 affects tumorigenesis in vivo, sh-Lacz/sh-Lucat1 transfected ACHN cells were inoculated into nude mice. Up to 49 days, there was a dramatic decrease in tumor volume and weight in the sh-Lucat1 group compared with sh-LacZ group (Figure 6A, 6B and 6C). Moreover, we also found that the tumors developed from sh-Lacz cells showed stronger Ki-67 expression than tumors formed from sh-Lucat1 and that tumors that developed from sh-Lucat1 cells showed a stronger p57 expression than tumors formed in the control, as detected by IHC analysis (Figure 6D). In the fresh subcutaneous tumor, we detected the expression of p57 by Western Blot. There was higher expression of p57 in sh-Lucat1 group than the sh-Lacz group (Figure 6E).
Figure 6: Lucat1 induced proliferation in vivo. (A) Photographs of tumors excised 7 weeks after inoculation of stably transfected cells ACHN into the armpit of nude mice. (B) Mean tumor volume measured by caliper on the indicated days. (C) Tumor weight of each nude mouse at the end of 7 weeks. (D) IHC staining for H&E, ki-67 and p57 in slices of sectioned implanted tumors formed by indicated cells. Original magnifcation was ×200. (E) Western blotting confirmed protein expression of p57 in indicated tumors. (F) Schematic Model of Long noncoding RNA Lucat1 acts as an oncogene in clear-cell renal cell carcinoma by binding to polycomb PRC2 complex and repressing p57.
These results suggest that the level of Lucat1 expression is significantly associated with the proliferation capacity of ccRCC in vivo and may through decrease the expression of p57 (Figure 6F).
DISCUSSION
Recently, many studies have shown that lncRNAs are frequently in dysregulation in various tumors and have multiple functions in a wide range of biological processes, such as proliferation, apoptosis, cell cycle arrest or cell migration and invasion [5, 33]. Given that Lucat1is induced by cigarette smoke and elevated in lung cancer cell lines and may contribute to cisplatin resistance in high grade serous ovarian cancer, we speculated that Lucat1 might involve in the ccRCC progression. With the aid of high-throughput techniques, especially the TCGA data, we identified and characterized Lucat1 as an oncogenic lncRNA in ccRCC. The results revealed that Lucat1 was upregulated in ccRCC tumorous tissues and markedly correlated with poor prognosis shorter overall survival and progression free survival. What’s more, silencing Lucat1 also impaired cell proliferation, migration and invasion, and cell cycle arrest in vitro, and inhibited tumorigenesis of ccRCC cells in vivo.
Many lncRNAs affect this aspect of cell transformation, as they interface with the epigenetic machinery influencing gene expression. Several can associate with the polycomb repressive complex 2 (PRC2), a chromatin repressor complex that catalyzes H3K27 trimethylation and which is tightly linked to the aberrant proliferation of cancer cells. Some lncRNAs affect gene expression through post-transcriptional processes, including splicing turnover, exportation or translation of mRNAs, as well as the stability and post-translational modification of proteins [33]. P57 is an inhibitor for cyclin-dependent kinase, and is deemed as a candidate for a tumor-suppressive gene that has been embroiled in numerous of cancers [31, 34, 35]. Bishoy A et.al examined the expression of p16, p21, p27, p53 and p57 in a large cohort of patients in ccRCC treated with extirpative therapy, and found that p57 was the most common altered marker, altered in 60% of 452 patients [35]. The TCGA data also showed that p57 exhibited a lower expression in tumor than in adjacent normal tissue in urothelial bladder cancer, breast cancer, colon and rectal adenocarcinoma, chromophobe renal cell carcinoma, ccRCC, lung adenocarcinoma etc. In addition, we also found that knockdown of Lucat1 increased the expression of p57 in ccRCC cells. Furthermore, we found that knockdown of Lucat1 inhibited cell proliferation, migration and invasion, and cell cycle arrest, while knockdown of p57 reversed the negative role of cell proliferation, and cell cycle arrest in ACHN and 786-O cell lines, which indicated p57 was a novel Lucat1 target, and Lucat1 could function as oncogene through suppressing p57 expression in ccRCC cells. Many lncRNAs affect the expression of p57 through recruiting EZH2 [14, 23, 36-38].
We applied the efficient CRISPR/Cas genome editing system targeting exon 2 of EZH2 to knock out this gene in 786-O cells [26] and found a higher expression of p57 in the EZH2-knockout cell lines. Moreover, we also found that in EZH2-knockout cell lines, Lucat1 cannot affect the expression of p57 and the proliferation of ccRCC, which further verified that Lucat1 decrease the expression of p57 through EZH2.
We have tried to detect the genes that can regulate the expression of Lucat1 in ccRCC, and found that the PTEN mutant group exhibited a higher expression of Lucat1 than the PTEN wide-type group from TCGA (supplementary Figure 4). Further analysis showed that Lucat1 has the strongest positive correlation with miR-21 (Supplemental Figure 4, Supplementary Table 2), a clearly reported target for PTEN. However, down-regulation of PTEN or miR-21 exhibited no obvious effects on the expression of Lucat1 (data not shown). Further investigation for genes relative to the regulation of the expression of Lucat1 is needed. Thus, it is reasonable for us to conclude that Lucat1 promotes ccRCC cell proliferation through recruitment of the Polycomb PRC2 complex and inhibition of p57.
In summary, our results provide strong evidence that Lucat1 is up-regulated in ccRCC tumor tissues and cell lines, and is associated with poor prognosis in ccRCC patients. Lucat1 acts as an oncogene in clear-cell renal cell carcinoma by binding to polycomb PRC2 complex and suppressing p57. This lncRNA may serve as a target for novel therapies in ccRCC.
MATERIALS AND METHODS
Clinical sample preparation
A total of 45 paired clear cell renal cell carcinoma and corresponding noncancerous tissues were obtained sequentially from patients undergoing radical nephrectomy during the period of 2010–2016 in Union hospital. Corresponding noncancerous tissues were acquired at least 5cm away from the tumor site. The study protocol was approved by the ethics committee of Huazhong University of Science and Technology and Union Hospital, and a written informed consent form was obtained from all participants involved in this study.
Cell culture
ACHN, 786-O and HK-2 cells were maintained in Dulbecco’s modified Eagle’s medium, supplemented with 10% fetal bovine serum and 2 mmol/L l-glutamine in a humidified atmosphere of 5% CO2 maintained at 37 °C. OS-RC-2 and CaKi-1 cells were cultured in RPMI-1640 supplemented with 10% fetal bovine serum and 2 mmol/L l-glutamine [24].
TCGA analysis
The RNA-seq data of 448 tumours and 67 matched normal samples were downloaded from The Cancer Genome Atlas (TCGA) Data Portal (http://ibl.mdanderson.org/tanric/_design/basic/query.html) [25].
Oligonucleotide, lentivirus, plasmid and shRNA
Oligos corresponding to the target sequences were annealed and cloned into the AgeI and EcoRI sites of the plko.l plasmid (Addgene).
Quantitative real-time PCR (RT-qPCR)
All RNA from tissues and cells were extracted with the TRIzol reagent (Invitrogen, Carlsbad, CA) according to the manufacturer’s protocol with modifcation. cDNAs were synthesized using Rever Ace qPCR RT Kit (TOYOBO). Real-time PCR was performed using SYBR Green Realtime PCR Master Mix (Roche) and the ABI ViiA7 qPCR System (Applied Biosystems). Chromatin immunoprecipitation (ChIP) was performed to investigate whether EZH2 and H3K27 binding to p57 promoter. ChIP assays were performed as described previously [26]. RNA binding protein immunoprecipitation (RIP) experiments were performed using Magna RIP Kit (Millipore, Catalog No.17-701) according to the manufacturer’s instructions and a previously published RIP-Chip protocol [27]. RIP assays were carried out as described previously [26]. All primers were in Supplementary Table 1.
Colony formation, cell proliferation assay, cell cycle analysis, cell migration and invasion
Colony formation was measured nearly two weeks after seeding 1000 cells per well in 6-well plates. Cell proliferation was estimated using the MTS according to manufacturer instructions. Fluorescence-activated cell-sorting (FACS) (BD, USA) analysis was done using propidiumiodide (PI) stains for cell-cycle analysis according to the manufacturer’s protocol with three replications [28]. Migration and invasion assays were performed using uncoated and Matrigel™ coated Transwell® inserts according to manufacturer instructions. 105 and 1.5*105 cells for ACHN migration and invasion. 3*104 and 5*104 cells for 786-O migration and invasion. All experiments were triplicated [26].
Xenograft subcutaneously
Tumorgenesis in nude mice was determined as described previously [5]. Five mice each were injected subcutaneously with prepared cells at a single site. Tumor onset was measured with calipers at the site of injection weekly at different times on the same day. All experiments were approved by the Animal Care and Use Committee of Tongji Medical College of Huazhong University of Science and Technology.
Immunohistochemical staining
As previously described [28].
Statistical analysis
All statistical analyses were carried out using SPSS 18.0 statistical software. Continuous data were compared using Student’s 2-tailed t-test. Data are represented as mean ± SEM. In all cases, p < 0.05 was considered statistically significant. *p < 0.05;**p < 0.01.
Abbreviations
ccRCC: clear cell renal cell carcinoma, Lucat1: lung cancer associated transcript 1, lncRNA: long non-coding RNA. PRC2: polycomb repressive complex 2.
Author contributions
Conception and design: HBX, XPZ, HMY. Development of the methodology: XHB. Acquisition of data: LB, WX. Analysis and interpretation of data: HLR, ZSS. Writing revision of the manuscript: HBX. Administrative, technical, and material support: YQ. Study supervision: KC, XPZ and HMY. All authors read and approved the final manuscript.
AVAILABILITY OF DATA AND MATERIAL
The datasets used and/or analyzed are available from the corresponding author on reasonable request.
CONFLICTS OF INTEREST
The authors declare that they have no competing interests.
FUNDING
This work was supported by grants from the National Natural Science Foundation of China (No. 81372760, 81072095, 81672528, 31401094, 31070142, 81272560.) and the National High Technology Research and Development Program of China (863 Program; 2012AA021101).
REFERENCES
1. Siegel R, Ma J, Zou Z, Jemal A. Cancer statistics, 2014. CA Cancer J Clin. 2014; 64: 9-29. doi: 10.3322/caac.21208.
2. Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, Jemal A, Yu XQ, He J. Cancer statistics in China, 2015. CA Cancer J Clin. 2016; 66: 115-32. doi: 10.3322/caac.21338.
3. Wu Y, Wang YQ, Weng WW, Zhang QY, Yang XQ, Gan HL, Yang YS, Zhang PP, Sun MH, Xu MD, Wang CF. A serum-circulating long noncoding RNA signature can discriminate between patients with clear cell renal cell carcinoma and healthy controls. Oncogenesis. 2016; 5: e192. doi: 10.1038/oncsis.2015.48.
4. Linehan WM. Genetic basis of kidney cancer: role of genomics for the development of disease-based therapeutics. Genome Res. 2012; 22: 2089-100. doi: 10.1101/gr.131110.111.
5. Xiao H, Tang K, Liu P, Chen K, Hu J, Zeng J, Xiao W, Yu G, Yao W, Zhou H, Li H, Pan Y, Li A, et al. LncRNA MALAT1 functions as a competing endogenous RNA to regulate ZEB2 expression by sponging miR-200s in clear cell kidney carcinoma. Oncotarget. 2015; 6:38005-15. doi: 10.18632/oncotarget.5357.
6. Cech TR, Steitz JA. The noncoding RNA revolution-trashing old rules to forge new ones. Cell. 2014; 157: 77-94. doi: 10.1016/j.cell.2014.03.008.
7. Guttman M, Russell P, Ingolia NT, Weissman JS, Lander ES. Ribosome profiling provides evidence that large noncoding RNAs do not encode proteins. Cell. 2013; 154: 240-51. doi: 10.1016/j.cell.2013.06.009.
8. Rinn JL, Kertesz M, Wang JK, Squazzo SL, Xu X, Brugmann SA, Goodnough LH, Helms JA, Farnham PJ, Segal E, Chang HY. Functional demarcation of active and silent chromatin domains in human HOX loci by noncoding RNAs. Cell. 2007; 129: 1311-23. doi: 10.1016/j.cell.2007.05.022.
9. Yang L, Lin C, Jin C, Yang JC, Tanasa B, Li W, Merkurjev D, Ohgi KA, Meng D, Zhang J, Evans CP, Rosenfeld MG. lncRNA-dependent mechanisms of androgen-receptor-regulated gene activation programs. Nature. 2013; 500: 598-602. doi: 10.1038/nature12451.
10. Poliseno L, Salmena L, Zhang J, Carver B, Haveman WJ, Pandolfi PP. A coding-independent function of gene and pseudogene mRNAs regulates tumour biology. Nature. 2010; 465: 1033-8. doi: 10.1038/nature09144.
11. Karreth FA, Reschke M, Ruocco A, Ng C, Chapuy B, Leopold V, Sjoberg M, Keane TM, Verma A, Ala U, Tay Y, Wu D, Seitzer N, et al. The BRAF Pseudogene Functions as a Competitive Endogenous RNA and Induces Lymphoma In Vivo. Cell. 2015; 161: 319-32. doi: 10.1016/j.cell.2015.02.043.
12. Khalil AM, Guttman M, Huarte M, Garber M, Raj A, Rivea Morales D, Thomas K, Presser A, Bernstein BE, van Oudenaarden A, Regev A, Lander ES, Rinn JL. Many human large intergenic noncoding RNAs associate with chromatin-modifying complexes and affect gene expression. Proc Natl Acad Sci U S A. 2009; 106: 11667-72. doi: 10.1073/pnas.0904715106.
13. Gupta RA, Shah N, Wang KC, Kim J, Horlings HM, Wong DJ, Tsai MC, Hung T, Argani P, Rinn JL, Wang Y, Brzoska P, Kong B, et al. Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature. 2010; 464: 1071-6. doi: 10.1038/nature08975.
14. Zhang E, He X, Yin D, Han L, Qiu M, Xu T, Xia R, Xu L, Yin R, De W. Increased expression of long noncoding RNA TUG1 predicts a poor prognosis of gastric cancer and regulates cell proliferation by epigenetically silencing of p57. Cell Death Dis. 2016; 7: e2109. doi: 10.1038/cddis.2015.356.
15. Hirata H, Hinoda Y, Shahryari V, Deng G, Nakajima K, Tabatabai ZL, Ishii N, Dahiya R. Long noncoding RNA MALAT1 promotes aggressive renal cell carcinoma through Ezh2 and interacts with miR-205. Cancer Res. 2015; 75: 1322-31. doi: 10.1158/0008-5472.CAN-14-2931.
16. Marin-Bejar O, Marchese FP, Athie A, Sanchez Y, Gonzalez J, Segura V, Huang L, Moreno I, Navarro A, Monzo M, Garcia-Foncillas J, Rinn JL, Guo S, et al. Pint lincRNA connects the p53 pathway with epigenetic silencing by the Polycomb repressive complex 2. Genome Biol. 2013; 14: R104. doi: 10.1186/gb-2013-14-9-r104.
17. Huo Y, Li Q, Wang X, Jiao X, Zheng J, Li Z, Pan X. MALAT1 predicts poor survival in osteosarcoma patients and promotes cell metastasis through associating with EZH2. Oncotarget. 2017; 8:46993-47006. doi: 10.18632/oncotarget.16551.
18. Nie F, Yu X, Huang M, Wang Y, Xie M, Ma H, Wang Z, De W, Sun M. Long noncoding RNA ZFAS1 promotes gastric cancer cells proliferation by epigenetically repressing KLF2 and NKD2 expression. Oncotarget. 2016; 8:38227-38238. doi: 10.18632/oncotarget.9611.
19. Zang C, Nie FQ, Wang Q, Sun M, Li W, He J, Zhang M, Lu KH. Long non-coding RNA LINC01133 represses KLF2, P21 and E-cadherin transcription through binding with EZH2, LSD1 in non small cell lung cancer. Oncotarget. 2016; 7: 11696-707. doi: 10.18632/oncotarget.7077.
20. Zhu XX, Yan YW, Chen D, Ai CZ, Lu X, Xu SS, Jiang S, Zhong GS, Chen DB, Jiang YZ. Long non-coding RNA HoxA-AS3 interacts with EZH2 to regulate lineage commitment of mesenchymal stem cells. Oncotarget. 2016; 7: 63561-70. doi: 10.18632/oncotarget.11538.
21. Thai P, Statt S, Chen CH, Liang E, Campbell C, Wu R. Characterization of a novel long noncoding RNA, SCAL1, induced by cigarette smoke and elevated in lung cancer cell lines. Am J Respir Cell Mol Biol. 2013; 49: 204-11. doi: 10.1165/rcmb.2013-0159RC.
22. Zheng ZG, Xu H, Suo SS, Xu XL, Ni MW, Gu LH, Chen W, Wang LY, Zhao Y, Tian B, Hua YJ. The Essential Role of H19 Contributing to Cisplatin Resistance by Regulating Glutathione Metabolism in High-Grade Serous Ovarian Cancer. Sci Rep. 2016; 6: 26093. doi: 10.1038/srep26093.
23. Renhua G, Yue S, Shidai J, Jing F, Xiyi L. 165P: Long noncoding RNA LUCAT1 is associated with poor prognosis in human non-small cell lung cancer and affects cell proliferation via regulating p21 and p57 expression. J Thorac Oncol. 2016; 11: S129. doi: 10.1016/S1556-0864(16)30275-1.
24. Zhang S, Zhong G, He W, Yu H, Huang J, Lin T. LncRNA Up-Regulated in Nonmuscle Invasive Bladder Cancer Facilitates Tumor Growth and Acts as a Negative Prognostic Factor of Recurrence. J Urol. 2016; 196:1270-8. doi: 10.1016/j.juro.2016.05.107.
25. Yang X, Karuturi RK, Sun F, Aau M, Yu K, Shao R, Miller LD, Tan PB, Yu Q. CDKN1C (p57) is a direct target of EZH2 and suppressed by multiple epigenetic mechanisms in breast cancer cells. PLoS One. 2009; 4: e5011. doi: 10.1371/journal.pone.0005011.
26. Guo J, Cai J, Yu L, Tang H, Chen C, Wang Z. EZH2 regulates expression of p57 and contributes to progression of ovarian cancer in vitro and in vivo. Cancer Sci. 2011; 102: 530-9. doi: 10.1111/j.1349-7006.2010.01836.x.
27. Shalem O, Sanjana NE, Zhang F. High-throughput functional genomics using CRISPR-Cas9. Nat Rev Genet. 2015; 16: 299-311. doi: 10.1038/nrg3899.
28. Chen K, Xiao H, Zeng J, Yu G, Zhou H, Huang C, Yao W, Xiao W, Hu J, Guan W, Wu L, Huang J, Huang Q, et al. Alternative splicing of EZH2 pre-mRNA by SF3B3 contributes to the tumorigenic potential of renal cancer. Clin Cancer Res. 2016; 23:3428-3441. doi: 10.1158/1078-0432.CCR-16-2020.
29. Huarte M. The emerging role of lncRNAs in cancer. Nat Med. 2015; 21: 1253-61. doi: 10.1038/nm.3981.
30. Guo H, Lv Y, Tian T, Hu TH, Wang WJ, Sui X, Jiang L, Ruan ZP, Nan KJ. Downregulation of p57 accelerates the growth and invasion of hepatocellular carcinoma. Carcinogenesis. 2011; 32: 1897-904. doi: 10.1093/carcin/bgr220.
31. Gayed BA, Youssef RF, Bagrodia A, Kapur P, Darwish OM, Krabbe LM, Sagalowsky A, Lotan Y, Margulis V. Prognostic role of cell cycle and proliferative biomarkers in patients with clear cell renal cell carcinoma. J Urol. 2013; 190: 1662-7. doi: 10.1016/j.juro.2013.06.037.
32. Sun CC, Li SJ, Li G, Hua RX, Zhou XH, Li DJ. Long Intergenic Noncoding RNA 00511 Acts as an Oncogene in Non-small-cell Lung Cancer by Binding to EZH2 and Suppressing p57. Mol Ther Nucleic Acids. 2016; 5: e385. doi: 10.1038/mtna.2016.94.
33. Aguilar R, Bustos FJ, Saez M, Rojas A, Allende ML, van Wijnen AJ, van Zundert B, Montecino M. Polycomb PRC2 complex mediates epigenetic silencing of a critical osteogenic master regulator in the hippocampus. Biochim Biophys Acta. 2016; 1859: 1043-55. doi: 10.1016/j.bbagrm.2016.05.009.
34. Zhang E, Yin D, Han L, He X, Si X, Chen W, Xia R, Xu T, Gu D, De W, Guo R, Xu Z, Chen J. E2F1-induced upregulation of long noncoding RNA LINC00668 predicts a poor prognosis of gastric cancer and promotes cell proliferation through epigenetically silencing of CKIs. Oncotarget. 2016; 7: 23212-26. doi: 10.18632/oncotarget.6745.
35. Xiao H, Xiao W, Cao J, Li H, Guan W, Guo X, Chen K, Zheng T, Ye Z, Wang J, Xu H. miR-206 functions as a novel cell cycle regulator and tumor suppressor in clear-cell renal cell carcinoma. Cancer Lett. 2016; 374:107-116. doi: 10.1016/j.canlet.2016.01.032.
36. Li J, Han L, Roebuck P, Diao L, Liu L, Yuan Y, Weinstein JN, Liang H. TANRIC: An interactive open platform to explore the function of lncRNAs in cancer. Cancer Res. 2015; 75:3728-37. doi: 10.1158/0008-5472.CAN-15-0273.
37. Keene JD, Komisarow JM, Friedersdorf MB. RIP-Chip: the isolation and identification of mRNAs, microRNAs and protein components of ribonucleoprotein complexes from cell extracts. Nat Protoc. 2006; 1: 302-7. doi: 10.1038/nprot.2006.47.
38. Xiao H, Zeng J, Li H, Chen K, Yu G, Hu J, Tang K, Zhou H, Huang Q, Li A, Li Y, Ye Z, Wang J, et al. MiR-1 downregulation correlates with poor survival in clear cell renal cell carcinoma where it interferes with cell cycle regulation and metastasis. Oncotarget. 2015; 6: 13201-15. doi: 10.18632/oncotarget.3915.