Abstract
Histone demethylase JMJD4 is a burgeoning tumor marker, which has been proven to be associated with colon cancer, but the role it plays in kidney cancer has not yet been investigated. In the present study, we evaluated whether JMJD4 can be a prognostic marker of patients with clear cell renal cell carcinoma (ccRCC) using data from public platform and in vitro experiments. Our results revealed that the expression of JMJD4 is higher in cancerous tissue than in normal tissues (). High expression of JMJD4 is associated with a poor overall survival (OS) of ccRCC as compared with low expression of JMJD4 (). JMJD4 showed significant relevance with M stage (), gender (), OS (0.018), disease-specific survival (DSS) (0.007), and percussion free interval (PFI) (0.041). Univariate and multivariate Cox analyses demonstrated that high JMJD4 expression had independent predictive value for OS in ccRCC patients (, , and ). Besides, in vitro experiments confirmed that high expression of JMJD4 can significantly promote the invasion ability (), cloning ability (), and proliferation () of renal cell carcinoma. In summary, high JMJD4 expression may be a prognostic marker in patients with kidney cancer.
1. Introduction
Renal cell carcinoma (RCC) is a common malignant tumor of the genitourinary system. It can be mainly divided into three subtypes, clear cell renal cell carcinoma (ccRCC), chromophobe, and papillary renal cell carcinoma, which comprise 65–70%, 15–20%, and 5–7% of all RCCs [1, 2]. RCC accounts for nearly 3% of all adult malignant tumors, and its incidence and related mortality have been increasing globally in the past decades [3, 4]. Usually, some patients have been prompted to have advanced RCC at the time of diagnosis due to the insidious onset of RCC. This tricky feature leads to poor prognosis and high mortality of patients with RCC. In advanced RCC, increased resistance to targeted therapies which may be mediated by gene disorders, angiogenesis, and decreased intake of tyrosine kinase inhibitors by cancer cells is common. In view of the hidden onset of RCC and the high resistance rate of targeted drug therapy for advanced RCC, it is urgent to explore effective biomolecular markers to contribute to the early diagnosis and prognostic information for RCC.
The modification of histones is one of the epigenetic mechanisms that regulate gene expression without changing the genotype [5, 6]. Hitherto, histone modifications that have been widely studied are histone methylation and acetylation modifications, which are closely related to the occurrence of many cancers [7, 8]. Histone demethylase is one of the most important modification enzymes involved in the regulation of chromatin function by catalyzing the removal of the methyl group on the N-terminal lysine residue of histone [9]. Many histone demethylases have been proven to promote the genesis and development of cancer, such as KDM3A, KDM4A, KDM4D, and JMJD6. Study has confirmed that KDM3A, which is highly expressed in colorectal cancer, is conducive to metastasis of colorectal cancer [10]. KDM3A has also been found to be increased in pancreatic tumors and promote cancer genesis by regulating expression of DCLK1 [11]. KDM4A could function as an oncogene in lung cancer by mediating expression of Myc through the Wnt/β-catenin signaling pathway [12]. In our previous study, we found that KDM4D inhibition can reduce the proliferation and development of kidney cancer. In addition, JMJD6 has been confirmed to be related to the progression of kidney cancer recently [13]. JMJD4 is a newly discovered histone demethylase which is homologous to JMJD6. However, the relationship between JMJD4 and tumor, especially kidney cancer, has not been reported yet.
This study is aimed at evaluating the prognostic value of JMJD4 in kidney cancer. We also conducted in vitro experiments to evaluate the effect of JMJD4 on ccRCC cells. Our results demonstrated that JMJD4 may be a cancer-promoting gene, and high expression of JMJD4 portends a poor prognosis in these patients.
2. Materials and Methods
2.1. Data Acquisition
The gene expression data and clinical data of ccRCC patients (Workflow Type: FPKM) were acquired from the online TCGA website (https://portal.gdc.cancer.gov/). Immunohistochemistry staining of ccRCC tissues was downloaded from The Human Protein Atlas (https://www.proteinatlas.org/). Patients with missing data, unknown staging status, and unknown survival status were excluded. Finally, a total of 539 patients with ccRCC with clinical characteristics information were obtained.
2.2. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes Pathway (KEGG) Enrichment Analysis
Differentially expressed genes in the low and high JMJD4 data set were identified and categorized into specific signaling pathways by GO (https://www.geneontology.org) and the KEGG database (http://www.genome.ad.jp/kegg/). The significant GO and KEGG pathways were identified by Fisher’s exact test. Adjusted value < 0.05 and value < 0.2 were considered as statistically significant.
2.3. Cell Culture
Human kidney cancer cell line (Caki-1), purchased from American Tissue Culture Collection (ATCC; Manassas, VA, USA), were cultured in McCoy’s 5A (HyClone, Logan, UT, USA) medium at 37°C, 5% CO2 in a humidified environment. All media are supplemented with 10% fetal bovine serum (Hyclone). Then, cells were transfected with pcDNA-3.1(+)-JMJD4 plasmid and control vector plasmid.
2.4. Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) Analysis
Caki-1 cells were harvested, and total RNA was extracted. Then, cDNA was synthesized by EntiLink™ 1st Strand cDNA Synthesis Kit (ELK Biotechnology, Wuhan, China) according to the manufacturer’s instruction. Quantitative real-time PCR was performed using the StepOne™ Real-Time PCR (Life technologies, Wuhan, China) with the following primer sequences: JMJD4 sense primer, 5-AGATGGTGTTTGTGCCCAGT-3; JMJD4 antisense primer, 5-CATGTTGGCCAGGTTGAAGC-3; β-actin sense primer, 5-CATCATCCCTGCCTCTACTGG-3; and β-actin antisense primer, 5-GTGGGTGTCGCTGTTGAAGTC-3.
2.5. Colony Formation Assay
Caki-1 cells were digested and resuspended in medium and seeded on a 6-well plate at a density of 500 cells per well. The culture was maintained in McCoy’s 5A medium and terminated when visible clones appeared. Next, the cells were harvested and fixed with 4% paraformaldehyde for 15 minutes. The cells were then stained with crystal violet and dried at room temperature. Finally, the number of clones which contains more than 50 cells was counted under the microscope.
2.6. Tumor Invasion Assays
Caki-1 cells were trypsinized, resuspended in serum-free culture medium, and counted. The upper chamber membranes of Transwell were evenly covered with Matrigel (BD science, New Jersey, USA) and incubated for 1 h before the cells were seeded. After 24 hours of conventional culture, cells in the upper chamber were gently wiped off by a cotton swab. Each chamber was fixed with methanol, stained with crystal violet, washed with PBS, and dried at room temperature. Six fields were randomly selected to take pictures under a microscope.
2.7. Wound-Healing Assays
Cells were harvested and seeded on six-well plates with per well. To simulate a wound, the cell surface was scratched by a 10 μl sterile plastic pipette tip with a straightedge when cells reached a confluent of 85% to 90%. Then, the cells were cultured in serum-free medium and incubated at 37°C. Overnight, wound images were collected under a microscope. The distance of the wound was measured and recorded, and the migration distance of cells in each group was calculated.
2.8. Cell Viability Assay
The proliferation of kidney cancer cells was determined by CCK-8 assays (Dojindo Molecular Technologies, Kumamoto, Japan). Generally, renal cancer cells were seeded into 96-well plates and incubated for 0 h, 24 h, 48 h, and 72 h. Then, the medium in each well was discarded and replaced by 100 μl of serum-free medium mixed with 10 μl of CCK-8 under dark condition. After incubated for 2 h, OD values at the absorbance of 470 nm were recorded.
2.9. Western Blot Analysis
Caki-1 cells were lysed by RIPA lysis buffer (Thermo Fisher Scientific Inc., MA, USA) mixed with protease and phosphatase inhibitors. Then, proteins were separated by SDS-PAGE and transferred to PVDF membranes. After transfer, membranes were put into 5% bovine serum albumin and sealed at 25°C for 2 hours. Next, membranes were washed thrice with TBST on a shaker for 5 minutes each time and incubated with TBST diluted primary JMJD4 antibody (AP1030a, Abcepta, San Diego, CA, USA) and anti-β-actin antibody (ab8227, Abcam, Cambridge, MA, USA) at 4°C. Overnight, membranes were washed by TBST and incubated with horseradish peroxidase- (HRP-) labeled secondary antibody (Jackson 1 : 2000) for 2 h at 25°C. Finally, membranes were reacted with electrochemical luminescence regents (Thermo Fisher Scientific Inc., MA, USA) and photographed in dark room.
2.10. Statistical Analysis
Clinical and RNA-seq data were organized and analyzed through R software 3.6.3 (Auckland, New Zealand) according to previous studies [14]. The box plots are used to analyze the expression level of JMJD4 between normal tissue and cancerous tissue among patients with various kinds of cancers including ccRCC. The relationship between JMJD4 expression and clinical characteristics of ccRCC patients was assessed by Wilcoxon signed rank test and Logistic regression. Univariate and multivariate Cox analyses were used to determine whether JMJD is an independent prognosis factor. Overall survival (OS), disease-specific survival (DSS), and progression-free interval (PFI) curves were plotted by the Kaplan–Meier method and compared with the log rank test. The diagnostic value of JMJD4 was assessed by plotting the received operating characteristic (ROC) curve and calculating the area under the curve. Significances were computed by Graphpad Prism 7.0 (Graphpad Software, Inc., CA, USA).
3. Results
3.1. Clinical Characteristics of ccRCC Patients
Totally, clinical data of 539 patients with ccRCC were downloaded from the official site of TCGA and analyzed. The characteristics of 539 patients are shown in Table 1. There are 353 male participants (65.5%) and 186 female participants (34.5%) with a median age of patients is 61. In terms of Fuhrman grade, 14 patients were grade I (2.6%), 235 patients were grade II (44.3%), 207 patients were grade III (39%), and 75 patients were grade IV (14.1%).
3.2. Expression of JMJD4 in ccRCC Tissues
JMJD4 was significantly higher in ccRCC tissues () than in normal tissues both in paired and unpaired box plot analysis (Figures 1(a) and 1(b)). The expression of JMJD4 was also higher in bladder cancer, breast cancer, cholangiocarcinoma, colon cancer, renal cancer, lung cancer, etc., than in normal tissue ( or , Figure 1(c)).

(a)

(b)

(c)
3.3. Correlation of High Expression of JMJD4 with Clinical Characteristics
The relationships between JMJD4 and clinical characteristics of ccRCC patients are displayed in Table 2. JMJD4 showed significant relevance with M stage (), gender (), OS (0.018), DSS (0.007), and PFI (0.041). As shown in Table 3, Logistic regression analysis confirmed that high expression of JMJD4 was significantly associated with M stage (M1 vs. M0: , , and ), pathological stage (stage III and stage IV vs. stage I and stage II, 95%, , and ), and gender (male vs. female, , , and ).
3.4. High Expression of JMJD Predicts a Worse Clinical Outcome
Online immunohistochemical staining images of JMJD4 in ccRCC patients are shown in Figure 2(a). We performed Kaplan–Meier survival analysis to assess the prognostic value of JMJD4, and results demonstrated that high JMJD4 expression was associated with poor OS (), DSS (0.005), and PFI (0.02) as shown in Figures 2(b)–2(d). Subgroup analysis by different clinical characteristics demonstrated that high JMJD4 expression was significantly associated with poor OS in patients with advance Fuhrman grade (), female patients (), patients from T1 to T3 stage (), patients no greater than 60 years old (), and patients with low or normal hemoglobin () as shown in Figures 3(a)–3(e). High level of JMJD4 is also connected with poor PFI in patients from T1 to T3 stage () and poor DSS in patients no greater than 60 years old () and female patients ().

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)
Univariate Cox regression analysis manifested that high JMJD4 expression was significantly correlated with poor OS (, , and ). Multivariate cox regression further validated that JMJD4 has independent prognostic value in patients with ccRCC (Table 4 and Figure 4, , , and ).

3.5. Diagnostic Value of JMJD4 in ccRCC Patients
The area under the ROC was 0.717, which indicates a high diagnostic value of JMJD for ccRCC patients (Figure 5(a)). The 1-, 3-, and 5-year survival probability of ccRCC patients are shown in Figure 5(b) by constructing a nomogram with the combination of JMJD4 and other clinical characteristics.

(a)

(b)
3.6. Genes Enriched in Different Signaling Pathways
We analyzed the RNA-seq data from TCGA database. Volcano map of differentially expressed genes is shown in Figure 6(a), and more genes were downregulated. The result of KEGG pathway enrichment is shown in Figure 6(b). KEGG analysis has shown that three most enriched pathways were estrogen signaling pathway, Staphylococcus aureus infection, and neuroactive ligand-receptor pathway. There are 31 genes enriched in tumor-related pathways as shown in Figure 6(c), suggesting that JMJD4 is closely related to cancer development. GO analysis and KEGG pathway enrichment were performed to investigate genes correlated with JMJD4 associated with ccRCC. Pathways that meet value < 0.05 and value < 0.2 are considered to be statistically significant. In biological process (BP) of GO analysis, genes were mainly associated with cornification, keratinization, keratinocyte differentiation, epidermal cell differentiation, and skin development (Figure 7(a)). Significant pathways in cellular component (CC) included keratin filament, extracellular region, intermediate filament, intermediate filament cytoskeleton, and extracellular region part as shown in Figure 7(b). Most enriched genes in molecular function (MF) include transmembrane transporter activity, channel activity, passive transmembrane transporter activity, transporter activity, and ion transmembrane transporter activity (Figure 7(c)). Genes enriched in pathway mentioned above include DSG3, KRT4, KRT5, KRT6, KRT14, KRT15, KRT16, and IVL.

(a)

(b)

(c)

(a)

(b)

(c)
3.7. JMJD4 Promoted Renal Cancer Cell Proliferation, Invasion, and Colony Formation
To exploit the role of JMJD4 in renal cell proliferation and development, plasmid containing pcDNA-3.1(+)-JMJD4 was transfected into the Caki-1 cell line. The Transwell invasion assay and colony formation assay were then used to unravel the effect of JMJD4 on cell invasion and colony formation, and the results showed that JMJD4 significantly facilitated the invasion and colony formation ability of Caki-1 cells (, Figure 8(a)). The efficiency of overexpression was verified by qRT-PCR and western blot (Figures 8(b) and 8(c)). As a result, western blot analysis and qRT-PCR showed that the expression of JMJD4 was significantly upregulated after transfection of pcDNA-3.1(+)-JMJD4. To investigate the effect that JMJD4 exerts on renal cell proliferation, we performed CCK-8 assays. The results showed that overexpression of JMJD4 in Caki-1 cells promoted the renal cell growth (, Figure 8(d)).

(a)

(b)

(c)

(d)
4. Discussion
JumonjiC (JmjC) domain-containing histone demethylase (JMJD) is a protein family that plays an important role in cancer growth and development by regulating gene transcription and chromatin structure. Among the members of family, JMJD6 was reported to be tumorigenic factor which involves in various kinds of cancers such as prostate cancer [15], breast cancer [16, 17], melanoma [18], colon cancer [19], and lung cancer [20]. JMJD4 is a paralog of JMJD6 which shares 34% sequence with JMJD6 [21]. Knockout of the JMJD6 gene in mice will significantly reduce the proliferation of NIH3T fibroblasts, while knocking out JMJD4 will cause the same effect by RNA interference [22]. Study has confirmed that JMJD4-knockout embryos are born at a normal Mendelian ratio which proved that JMJD4 is unessential in embryogenesis [23]. Currently, only one study has reported that JMJD4 expression could be a prognostic indicator for patients with colon cancer [24]. The interaction between JMJD4 and other tumors still needs more research to reveal. Therefore, we investigate the relationship between JMJD4 and renal cancer.
In this study, we first use bioinformatics analysis to assess the role of JMJD4 in ccRCC from online public database such as TCGA and TIMER. The results of box plot analysis demonstrate that high expression of JMJD4 is a common feature of multiple human cancers, including ccRCC. And Logistic regression analysis indicated that high JMJD4 is closely related to M stage, gender, and prognosis in ccRCC patients. Then, to evaluate the prognostic value of JMJD4 in ccRCC patients, we performed Kaplan–Meier survival analysis, and results showed that patients with high JMJD4 expression had poor OS, DSS, and PFI. Moreover, Cox regression analysis confirmed that high JMJD4 expression has independent prognostic value for OS in ccRCC patients. Other clinicopathological characteristics are also related to the poor OS of ccRCC, including advanced Fuhrman grade, M stage, T stage, and age, as shown in Figure 4. Next, we performed ROC analysis and compute the AUC. The results suggest that high JMJD4 expression has certain diagnostic value in ccRCC patients (Figure 5(a)). Besides, to help practitioners predict the risk of individual death and optimize treatment options, we also constructed a prognostic nomogram (Figure 5(b)).
To further evaluate the mechanism of high expression of JMJD4 in promoting the progression of renal cancer, we performed GO and KEGG analysis. We found that the high JMJD4 expression was associated with collecting duct acid secretion, neuroactive ligand-receptor interaction, and fat digestion and absorption pathway in KEGG analysis. GO analysis shows that keratinization, keratinocyte differentiation, and epidermal cell differentiation were most enriched genes in BP. Skin development extracellular region, intermediate filament, and intermediate filament cytoskeleton were most enriched pathways in CC, while transmembrane transporter activity, channel activity, and passive transmembrane transporter activity were most enriched pathways in MF. The main genes enriched in these pathways include DSG3, KRT5, KRT6, and KRT14. DSG3 was reported to promote the growth and invasion of cancer cells [25]. KRT5 deficiency can prevent the migration of ovarian cancer cells, which prompts that it may act as an oncogene [26]. KRT6 is associated with notch1 signaling and contributes the progression of renal cancer [27]. KRT14 was proved to be a key regulator in metastasis of breast cancer [28]. The above data shows that the genes related to JMJD4 are mainly cancer-promoting genes, indicating that JMJD4 may promote the invasion and progression of cancer. Furthermore, to clarify the effect of JMJD4 on kidney cancer, we transfected the pcDNA-3.1(+)-JMJD4 plasmid into Caki-1 cells and verified the overexpression by qRT-PCR and western blot (Figures 8(a) and 8(b)). Then, we performed Transwell invasion experiments and plate colony formation experiments. The results showed that JMJD4 overexpression significantly promoted the invasiveness and colony formation ability of Caki-1 cells (Figures 8(c), ). We also conducted the CCK-8 experiment, and the results demonstrate that JMJD4 can significantly promote the proliferation of Caki-1 cells (Figures 8(d), ). The above in vitro data prove that JMJD4 acts as a cancer-promoting gene in renal cancer cells. One limitation of our research is that we did not conduct in vivo experiments for verification. Subsequent experiments can be performed on animal experiments to verify the cancer-promoting effect of JMJD4 in vivo.
In summary, our study is the first research that focuses on the role of JMJD4 played in ccRCC. We validated the predictive value and diagnostic value of JMJD4 in ccRCC patients, as well as the effect of JMJD4 on promoting the invasion and proliferation of ccRCC cells. Our results showed that high expression of JMJD can be a potential prognostic marker for ccRCC patients.
Data Availability
The datasets generated during the current study are available from the corresponding author on reasonable request.
Conflicts of Interest
The authors declare no conflict of interest.
Authors’ Contributions
Hao Yan and Yewei Bao contributed equally to this work.
Acknowledgments
This research was granted by the Shanghai Science and Technology Development Foundation (17JC1400904).
Supplementary Materials
Supplementary 1. Supplemental material 1: details of the enriched genes in KEGG analysis.
Supplementary 2. Supplemental material 2: details of the enriched genes in biological process of GO analysis.
Supplementary 3. Supplemental material 3: details of the enriched genes in molecular function of GO analysis.
Supplementary 4. Supplemental material 4: details of the enriched genes in cellular component of GO analysis.