MKRN1/2 serve as tumor suppressors in renal clear cell carcinoma by regulating the expression of p53

Abstract

BACKGROUND:

Kidney renal clear cell carcinoma (KIRC) belongs to renal cell carcinoma which is a very aggressive malignant tumor with poor prognosis and high mortality. The MKRN family includes three members MKRN1, MKRN2 and MKRN3, which are closely related to cancers, and have been involved in many studies.

OBJECTIVE:

This study aimed to explore the roles of MKRN family in KIRC.

METHODS:

The expression of MKRNs was analyzed using the UALCAN database, prognostic analysis was performed with the GEPIA2 and Kaplan-Meier Plotter database, and correlation analysis was assessed by GEPIA2. The CCK-8 and colony formation assay were performed to detect cell proliferation, wound healing assays were performed to detect cell migration, cell cycles were detected by flow cytometry analysis, GST pull-down and co-immunoprecipitation assays were performed to detect the interaction of proteins, and the expression of MKRNs, p53 and other proteins were detect by immunoblotting analysis or quantitative PCR (qPCR).

RESULTS:

MKRN1 and MKRN2 were lowly expressed in KIRC samples compared to the corresponding normal tissues, and KIRC patients with high levels of MKRN1 and MKRN2 showed higher overall survival (OS) and disease free survival (DFS) rates. The overexpression of MKRN1 and MKRN2 inhibited the proliferation of human KIRC cells by arresting the cell cycles, but shows little effect on cells migration. The expression of MKRN1 and MKRN2 are correlated, and MKRN1 directly interacts with MKRN2. Moreover, both MKRN1 and MKRN2 were closely correlated with the expression of TP53 in KIRC tumor, and promoted the expression of p53 both at protein and mRNA levels.

CONCLUSIONS:

Our study suggests that MKRN1 and MKRN2 serve as tumor suppressors in KIRC, and act as promising therapeutic targets for KIRC treatment.

Keywords

Kidney renal clear cell carcinoma MKRN1 MKRN2 proliferation p53

Abbreviation

KIRC	Kidney renal clear cell carcinoma
OS	overall survival
DFS	disease free survival
MKRN	MAKORIN
Co-IP	Co-immunoprecipitation
qPCR	quantitative PCR
IB	immunoblotting
PI	propidiumiodide
EMT	epithelial-mesenchymal transition

1. Introduction

The MAKORIN (MKRN) gene family encodes proteins with a makorin-specific Cys-His arrangement, a RING zinc finger and several zinc-finger motifs [1, 2]. Three members, namely MKRN1, MKRN2 and MKRN3, are included in this family, which share highly homologous amino acid sequences, especially in the zinc finger motifs, suggesting that they may have similar functions [1]. All these MKRN family proteins can function as E3 ligase, and many substrates for them have been identified.

MKRN1 has been involved in regulating the length of telomere, and the degradation of cell cycle regulator p21 and tumor suppressor protein p53 [3, 4]. MKRN1 shows both the negative and positive effects on RNA polymerase II-dependent transcription, and was identified as a repressor for androgen receptor, c-Jun and retinoic acid (RA) receptor transcriptional activity [5]. MKRN1 has been shown to bind with and ubiquitylate adenomatous polyposis coli (APC) which is a tumor-suppressor function, increasing its proteasome degradation and favorably regulating Wnt/-catenin-mediated biological activities [6]. MKRN1 promotes the ubiquitination and proteasome degradation of Fas-associated protein with death domain (FADD) regulates death receptor-mediated apoptosis and necroptosis in HeLa cells [7].

MKRN2 promotes the degradation of inflammatory responsed p65, regulates the proliferation and migration of human neuroblastoma cells by ubiquitinating IGF2BP3, and inhibits the metastasis of non-small cell lung cancer by regulating the PI3K/Akt signaling pathway [8, 9, 10].

MKRN3 is frequently mutated in non-small cell lung cancers (NSCLCs), and genomic MKRN3 aberrations are significantly enriched in NSCLC samples. MKRN3 modulates cell proliferation through PABPC1 nonproteolytic ubiquitination and subsequently, PABPC1-mediated global protein synthesis [11], which suggested that MKRN3 act as a tumor suppressor. However, elevated expression of MKRN3 was observed in in squamous cell carcinoma of the head and neck [12].

Clear cell renal carcinoma (KIRC) is the most common subtype of renal cell carcinoma with a poor prognosis and few treatment options [13]. Due to distant metastases, the 5-year survival rate of KIRC patients is only 10–20% [14]. Many targets for the treatment of advanced KIRC have been discovered and used in clinical practice, but few patients have benefited [15], so other targets are urgently needed. The MKRN family members are closely related to cancers and were involved in many studies as above descrbed, but their roles in KIRC have not been illustrated. The aim of this study is to explore the roles of MKRN1, MKRN2 and MKRN3 in human KIRC.

2. Materials and methods

2.1 Analyses involving publicly available databases

2.1.1 University of alabama cancer database (UALCAN)

ULACAN (http://ualcan.path.uab.edu) is an easy-to-use interactive portal for in-depth analysis of TCGA gene expression data [16], and were used to assess the expression of MKRN family members in different type of cancers, especially the mRNA and protein levels of MKRN1 and MKRN2 in human kidney renal clear cell carcinoma (KIRC). Student’s $t$ -test was used to generate a $p$ -value for the expression analysis.

2.1.2 Gene expression profiling interactive analysis 2 (GEPIA2)

GEPIA2 (http://gepia2.cancer-pku.cn) is a powerful tool for analyzing gene expression based on RNA sequencing data from Cancer Genome Atlas (TCGA) and Tissue Genotype Expression Database (GTEx) [17]. In this study, GEPIA2 was used for severity analysis, prognostic value analysis, and correlation analysis. Student’s $t$ -test was used to generate a $P$ -value for the pathological stage analysis, and Pearson were used to generate a $P$ -value and correlation R-value for correlation analysis. For prognostic value analysis, hazard ratios with 95% confidence intervals and log-rank $P$ -values were also computed.

2.1.3 Kaplan-meier plotter

The Kaplan-Meier plotter (https://kmplot.com/analysis) [18] is also used to analyze the prognostic value of MKRNs in KIRC patients. Detailed information on the number of high risk cases and the median mRNA expression level were collected. The hazard ratio with 95% confidence interval and the log-rank $P$ -value are also calculated.

2.1.4 cBioPortal

cbiopportal (http://www.cbioportal.org) is a public database that provides visualization and multidimensional cancer genome data [19]. The genetic alterations of MKRNs in KIRC patients were analyzed using cBioPortal based on The Cancer Genome Atlas (TCGA) database.

2.2 Cell culture and transfection

The human embryonic kidney cell line HEK293T and human kidney renal clear cell carcinoma cell lines (786-O and OSRC-2) were purchased from the cell bank of Chinese Academy of Sciences (CAS). HEK293T cells were cultured in DMEM (high glucose), while786-O and OSRC-2 were cultured in RPMI1640, supplemented with 10% fetal bovine serum (FBS) and penicillin/streptomycin (all from Gibco, USA) in a 37 ${{}^{\circ}}$ C humidified atmosphere of 5% CO ${}_{2}$ . Plasmids used in this study were transfected into cells using Lipofectamine ${}^{\@setsize{\scriptsize}{8pt}{\viipt}{\@viipt}\textregistered}$ 2000 (Thermo Fisher, USA) according to the manufacturer’s instructions. Briefly, cells were seeded into 6-well-plates, transfections were done when cells to be 90% confluent. 3 $\mu$ g plasmids DNA and 6 $\mu$ l Lipofectamine ${}^{\@setsize{\scriptsize}{8pt}{\viipt}{\@viipt}\textregistered}$ 2000 Reagent were diluted in 150ul Opti-MEM ${}^{\@setsize{\scriptsize}{8pt}{\viipt}{\@viipt}\textregistered}$ Medium (Gibco, USA) respectively, incubated at room temperature for 5 min, and then mixed together for 20 min before added to cells. 24 hours after transfection, cells were screened by 2 $\mu$ g/ml puromycin (Beyotime, China) for seven days to established the stably transfected cell lines.

2.3 Plasmids construction

The plasmids containing MKRN1 and MKRN2 were amplified from the cDNA of human HEK293T cells using specific primers in Table S1, and inserted into pCDH, pCDNA3.0, pET28a or pGEX4T-1 vectors that digested with BamH I or Xho I, and recombination with amplified fragments using the ClonExpress II One Step Cloning Kit (Vazyme, China) according to the manufacturer’s instructions. Three shRNAs for MKRN1 and four shRNAs for MKRN2 were synthesized as oligos (Tingke, Shanghai, China), annealed and inserted into the pLKO.1 vector which was digested with restriction endonucleases EcoR I and Age I, the specific sequences for shRNAs seen in Table S2.

2.4 Cell proliferation assay

Five thousand 786-O or OSRC-2 cells that stably transfected with indicated plasmids that expressing genes or shRNAs were seeded into 96-well plates. Incubate the cells with Cell Counting Kit-8 (CCK8) solution (Beyotime, China) for 3 hours at the specified time, and the product was quantified by spectrophotometry at a wavelength of 450 nm in a microplate reader (Bio-Rad, USA). The experiments were conducted with eight replicates and repeated three times.

2.5 Colony formation assay

One thousand cells that stably transfected with pCDH, pCDH-MKRN1 or pCDH-MKRN2 were seeded into 6-well plates. Seven days later, Cells were fixed with 4% paraformaldehyde (Merck, Germany) at room temperature for 15 min, and stained with 0.2% crystal violet (Beyotime, China) seven days later. Images were obtained and the number of colonies was counted and calculated.

2.6 Cell cycle assay

786-O or OSRC-2 cells that stably transfected with pCDH, pCDH-MKRN1 or pCDH-MKRN2, were plated into 6-well plates and cultured overnight. Next day, the cells were collected and fixed with 70% ethanol for 1 hour, and then treated with 0.2% Triton X-100, 100 $\mu$ g/mL RNase A, and 50 $\mu$ g/mL propidiumiodide (PI) for 45 min and analyzed by flow cytometry as previously described [20].

2.7 Wound healing assay

786-O cells that stably transfected with pCDH, pCDH-MKRN1 or pCDH-MKRN2 were seeded into 6-well plates at a density of 5 $\times$ 10 ${}^{4}$ cells/ml. After 24 hours, draw a straight line with a sterile pipette tip to remove cells in the area. Images were taken using a microscope at 0 hours and 24 hours. Measure the degree of cell migration up to the center of the scratch and calculate the distance using the method previously described [21].

2.8 Recombinant protein purification

GST tagged MKRN1 or His6- tagged MKRN2 were expressed in BL21 E. coli cells. After isopropyl- $\beta$ -d-thiogalactopyranoside (IPTG, Sangon, China) induction, E. coli cells were pelleted by centrifuge, lysed in PBS buffer and incubated with glutathione or Ni ${}^{2+}$ beads (Sangon, China) to enrich the GST $-$ or His6-tagged proteins, followed by elution with 15 mM reduced L-glutathione or 500 mM imidazole (Sangon) that dissolved in PBS buffer, and then dialysis in PBS buffer supplemented with 15% glycerol before being aliquoted and preserved at $-$ 80 ${}^{\circ}$ C.

2.9 GST pull-down assay

Purified GST or GST-MKRN1 (50 $\mu$ g), MKRN2-His6 proteins (50 $\mu$ g) and Glutathione Sepharose4B (Sangon) were incubated at 4 ${}^{\circ}$ C overnight in 800 $\mu$ l GST pull-down buffer (20 mM Tris-Cl, 100 mM NaCl, 1 mM EDTA, 5 mM MgCl ${}_{2}$ , 1 mM DTT, 1% NP-40, 10 mg/mL of BSA, pH 7.4). These beads are then pelletized and washed three times with a pull-down buffer. The recovered immunoprecipitates were boiled in 2xSDS-PAGE protein loading buffer for 10 min, and then subjected to immunoblotting analysis with the indicated antibodies.

2.10 Co-immunoprecipitation (Co-IP) and immunoblotting

For Co-IP, cells that co-transfected with pCDNA3.0-MKRN1-HA and pCDNA3.0-MKRN2-Flag were lysed in 1 ml Co-IP buffer (50 mM Tris-HCl, 150 mM NaCl, 5 mM EDTA and 1% NP-40, pH 7.4) supplemented with a protease inhibitor cocktail (Roche, Switzerland). The cell lysate was then centrifuged and incubated with anti-Flag affinity gel (Sigma, USA) at 4 ${}^{\circ}$ C overnight. The beads are pelletized and washed 3 times with Co-IP buffer. The immunoprecipitates were enriched denatured at 99 ${}^{\circ}$ C for 10 min in 50 $\mu$ l 2X SDS-PAGE loading buffer. Input, immunoprecipitation and other cell lysates are processed by 10% SDS-PAGE and transferred to PVDF membrane (Bio-Rad, USA). The membrane was blocked with 10% skimmed milk at RT for one hours, and stained with the specified antibodies as follows: anti-MKRN1 (1:1000 dilution; ab72054, Abcam, UK), anti-MKRN2 (1:800 dilution; 12238-1-AP, Proteintech, China), anti-CCNA1 (1:1000 dilution; ab270940, Abcam), anti-CCNB1 (1:800 dilution; 28603-1-AP, Proteintech), anti-CCND1 (1:200 0 dilution; 60186-1-Ig, Proteintech), anti-CASPASE3 (1:500 dilution; ab32351, Abcam), anti-PARP(1:1000 dilution; 80174-1-RR, Proteintech), anti-Vimentin (1:10000 dilution; 10366-1-AP, Proteintech), anti-E-cadherin (1:1000 dilution; 20874-1-AP, Proteintech), anti-GAPDH (1:8000 dilution; 60004-1-Ig, Proteintech), anti-Flag (1:3000 dilution; 80010-1-RR, Proteintech), anti-HA (1:1000 dilution; 51064-2-AP, Proteintech), anti-GST (1:5000 dilution; HRP-66001, Proteintech), and anti-His (1:1000 dilution; 66005-1-Ig, Proteintech). Then the membranes were incubated with HRP-conjugated goat anti-mouse IgG (1:5000 dilution; SA00001-1, Proteintech) or goat anti-rabbit IgG (1:5000 dilution; SA00001-2, Proteintech) secondary antibodies, and the signals were visualized using a Tanon 5200 Imaging System (Tanon, China).

2.11 Quantitative PCR (qPCR)

Total RNAs were extracted from cells using total RNA kit (Tiangen, China). Complementary DNA (cDNA) was synthesized using ReverTra Ace RT Master mix (Toyobo, Japan). Quantitative PCR (qPCR) assay was performed to determine the relative abundances of TP53 and GAPDH mRNAs using specific primers in Table S3, stained by SYBR Green (Tiangen, China) on the CFX96 real-time PCR system (Bio-Rad, USA). The relative abundances of TP53 mRNAs were normalized to that of GAPDH mRNAs, using the 2 ${}^{-\Delta\Delta\text{Ct}}$ method as previously described [22]. All experiments were repeated three times.

2.12 Statistical analysis

Data of experiments are expressed as mean $\pm$ standard deviation (SD) and analyzed by student’s $t$ -test using GraphPad Prism 5 (GraphPad Software Inc., USA). ${}^{*}P<$ 0.05 was considered to be significantly different, ${}^{**}P<$ 0.01 was considered to be very significantly different.

Figure 1.

Low expression levels of MKRN1 and MKRN2 in human kidney renal clear cell carcinoma (KIRC). (A) The mRNA expressions of MKRN1, MKRN2 and MKRN3 in human kidney renal clear cell carcinoma (KIRC) and corresponding normal tissues were analyzed using the UALCAN database. ${}^{*}P<$ 0.05 was considered to be significant difference. (B) The proteins expressions of MKRN1 and MKRN2 in KIRC and corresponding normal tissues were analyzed using the UALCAN database. Z-values represent standard deviations from the median across samples for the given cancer type. Log2 Spectral count ratio values from CPTAC were first normalized within each sample profile, and then normalized across samples. (C) The correlations between the expression of MKRN1, MKRN2 or MKRN3 and tumor stage in KIRC patients were analyzed using the GEPIA2 database. The F value is a value on the F distribution that is a theoretical distribution, which is used in analysis of variance (ANOVA). The Pr ( $>$ F) value is the $p$ -value of the F-statistic.

Figure 2.

Prognostic values of MKRNs in KIRC. (A, B) The overall survival (OS) map (A) and disease free survival (DFS) map (B) of MKRNs in TCGA cancers analyzed using the GEPIA2 database. HR, hazard rate ratio $=$ treatment hazard rate/placebo hazard rate. HR, hazard rate ratio $=$ treatment hazard rate/placebo hazard rate. The hazard ratio is constant under the Cox proportional hazard model. (C) The expression levels of MKRNs and OS rates in KIRC patients were analyzed using the GEPIA2 database. The p (HR) is used to reject the null hypothesis that HR $=$ 1. (D) The expression levels of MKRNs and DFS rates in KIRC patients were analyzed using the GEPIA2 database.

3. Results

3.1 Low expression levels of MKRN1 and MKRN2 in human kidney renal clear cell carcinoma (KIRC)

A pan-cancer analysis of MKRN1, MKRN2 and MKRN3 in cancers was curated using the UALCAN database (Fig. S1A–C), and MKRN1 and MKRN2 were found lowly expressed in human kidney renal clear cell carcinoma (KIRC) samples compared to the corresponding normal tissues both at mRNA and protein levels (Fig. 1A and B). Since there is no data in the UALCAN CPTAC database, the protein expression of MKRN3 is not shown. The mRNA expression of MKRN3 is low in KIRC samples, and no obvious difference was observed in human KIRC tissues and the corresponding normal tissues (Fig. S1C and A).

Then the correlation between MKRNs (MKRN1, MKRN2 and MKRN3) mRNA expression and the pathological stage of KIRC patients were assessed using the GEPIA2 database. All three groups were highly variable, and the mRNA expression levels of MKRN1 and MKRN2 were gradually decreased with the progression of KIRC (Fig. 1C).

3.2 High levels of MKRN1 and MKRN2 predicts better survival in KIRC patients

To assess the value of differentially expressed MKRN family members in the progression of cancers, correlations between different MKRNs and clinical outcomes were analyzed using GEPIA2 (Fig. 2A and B). KIRC patients with high mRNA expression levels of MKRN1 and MKRN2 were significantly associated with better overall survival (OS) and disease free survival (DFS) rates, while this effect was not observed on MKRN3 (Fig. 2C and D). Furthermore, the Kaplan-Meier plotter were also used to analyze the prognostic values of MKRNs in patients with KIRC, high MKRN1 mRNA expression was significantly correlated with better OS and DFS rates (Fig. S2A and B). So, MKRN1 and MKRN2, but not MKRN3, were chosen for further study.

Figure 3.

MKRN1 and MKRN2 inhibited the proliferation of human KIRC cells. (A) Protein expression profiles of MKRN1 and MKRN2 in human KIRC cell lines and human embryonic kidney 293T (HEK293T) cells detected by immunoblotting analysis. HEK293T cells acted as normal control. (B, C) MKRN1 and MKRN2 inhibited the proliferation of human KIRC cells. 786-O and OSRC-2 cells that stably transfected with pCDH (empty vector), and pCDH-MKRN1 (B) or pCDH-MKRN2 (C) were seeded into 96-well plates and the cell viability were detected by the CCK-8 assay at indicated time points. Three independent experiments. (D, E) MKRN1 and MKRN2 inhibited the colony formation of human KIRC cells. 786-O and OSRC-2 cells that stably transfected with pCDH (empty vector), and pCDH-MKRN1 (D) or pCDH-MKRN2 (E) were seeded into 6-well plates, the cells colonies were stained with crystal violet, and photos were taken. The colony numbers were counted and calculated.Three samples each group.

Figure 4.

The expression of MKRN1 and MKRN2 are correlated, and MKRN1 interacts with MKRN2. (A) Correlation analysis of MKRN1 and MKRN2 mRNA levels in KIRC tumor, KIRC normal and kidney cortex using the GEPIA2 database. (B) MKRN1 interacted with MKRN2 in vitro. Recombinant GST-tagged MKRN1 and His6-tagged MKRN2 proteins were purified, and GST pull-down assay was performed before subjected to immunoblotting analysis. PD, pull-down. (C) MKRN1 formed a complex with MKRN2 in HEK293T cells. Cells were co-transfected with plasmids that expressed Flag tagged MKRN2 and HA tagged MKRN1, and the lysates were immunoprecipitated using anti-Flag affinity gels before subjected to immunoblotting analysis with indicated antibodies. Co-IP, Co-immunoprecipitation.

Figure 5.

MKRN1 and MKRN2 promoted the expression of p53. (A, B) Correlation analysis of MKRN1 (A) or MKRN2 (B) and TP53 mRNA levels in KIRC tumor, KIRC normal and kidney cortex using the GEPIA2 database. (C) MKRN1 and MKRN2 promoted the expression of p53 protein in a dose-dependent manner. OSRC-2 cells were transiently transfected with different amounts of pCDH (empty vectors), pCDH-MKRN1or pCDH-MKRN2, and subjected to immunoblotting analysis using indicated antibodies 48 hours later. (D) MKRN1 and MKRN2 synergistically promoted the expression of p53 protein. OSRC-2 cells were transiently transfected with pCDH-MKRN1 or pCDH-MKRN2, and subjected to immunoblotting analysis using indicated antibodies 48 hours later.

3.3 MKRN1 and MKRN2 inhibited the proliferation of human KIRC cells

The genetic alterations of MKRNs in KIRC patients were analyzed using the cBioPortal online tool. As showed in Fig. S3A and B, MKRNs were altered in 15 of 354 KIRC patients, accounting for 4.2%. MKRN1 were altered in 3 samples of 354 patients, including 2 gene mutations and 1 amplification; MKRN2 were altered in 9 samples of 354 patients, all were deep deletion; MKRN3 were altered in 3 samples of 354 patients, all are gene mutations. The results might suggest that gene mutations of MKRN1 and MKRN2 were seldom found in KIRC patients, so their expression changes become the focus of our further study.

The protein levels of MKRN1 and MKRN2 were lower in the two human KIRC cell lines (786-O and OSRC-2) compared to human embryonic kidney 293T (HEK293T) cells (Fig. 3A). To investigate the effect of MKRN1 and MKRN2 on KIRC cells, 786-O and OSRC-2 that stably transfected with MKRN1 or MKRN2 were established (Fig. S4A). As showed in Fig. 3B and C, significant cells growth arrest were observed in MKRN1 or MKRN2 overexpression cell lines (pCDH-MKRN1or pCDH-MKRN2) compared with that of control groups (pCDH), which was consistent with the results of colony formation assay (Fig. 3D and E). Flow cytometry analysis revealed that MKRN1 and MKRN2 inhibited the cell cycle progression of both786-O and OSRC-2 cells (Fig. S4B). Immunoblotting analysis found that MKRN1 and MKRN2 inhibited the expression of cell cycle proteins including CCNA1, CCNB1 and CCND1, but showed little effect on apoptosis protein markers (CASPASE3 and PARP) and epithelial-mesenchymal transition (EMT) protein markers (Vimentin and E-cadherin) (Fig. S4C). Wound healing assays were performed with 786-O cells that stably transfected with MKRN1 or MKRN2, and no obvious effect on the migration of human these cells were observed (Fig. S4D).

Furthermore, the shRNAs targeting MKRN1 and MKRN2 were designed and the knockdown efficiency was tested (Fig. S5A and C). The shMKRN1-1 and shMKRN1-3 for MKRN1, and shMKRN2-1 and shMKRN2-2 for MKRN2 were chosen for further study. The CCK-8 assay indicated that the ablation of MKRN1 and MKRN2 promoted the proliferation of 786-O and OSRC-2 cells (Fig. S5B and D).

3.4 The expression of MKRN1 and MKRN2 are correlated, and MKRN1 interacted with MKRN2

A correlation analysis was performed between MKRN1 and MKRN2 using the GEPIA2 database. The mRNA expression of MKRN1 and MKRN2 were closely correlated in KIRC tumor, KIRC normal and kidney cortex (Fig. 4A). Further study indicated that MKRN1 directly interacted with MKRN2 in vitro as detected by GST pull-down assay (Fig. 4B), and MKRN1 formed a complex with MKRN2 in HEK293T cells (Fig. 4C). These data suggested that MKRN1 and MKRN2 are synergistic in their biological functions.

3.5 MKRN1 and MKRN2 promoted the expression of p53

A correlation analysis was performed between MKRN1 or MKRN2 and TP53 using the GEPIA2 database. Both MKRN1 and MKRN2 were closely correlated with the mRNA expression of TP53 in KIRC tumor, KIRC normal and kidney cortex (Fig. 5A and B). MKRN1 or MKRN2 promoted the expression of p53 in a dose-dependent manner both at protein and mRNA levels (Figs 5C and S6A). We also found that MKRN1 and MKRN2 synergistically promoted the expression of p53 both at protein and mRNA levels (Figs 5D and S6B). Further study indicated that the ablation of MKRN1 or MKRN2 inhibited the expression of p53 (Fig. S6C and D). These results suggested that MKRN1 or MKRN2 regulated the expression of p53 mainly at transcription level.

4. Discussion

Kidney renal clear cell carcinoma (KIRC) is a subtype of renal cell carcinoma, which is a very aggressive malignant tumor shows poor prognosis and high mortality [15], so promising biomarkers or therapeutic targets for KIRC are urgently needed. MKRN1 and MKRN2 is lowly expressed in KIRC samples compared to the corresponding normal tissues (Fig. 1) and patients with high levels of MKRN1 and MKRN2 shows higher overall survival (OS) and disease free survival (DFS) rates (Fig. 2), which suggest that they might act as an important prognostic indicators, even good therapeutic targets for KIRC treatment.

MKRN1 and MKRN2 are widely expressed in human tissue [1, 23], and we found that they are extensively expressed across TCGA cancers (Fig. S1A and B), while the expression of MKRN3 is low in most tumors (Fig. S1C). The expression of MKRN3 is tissue-specific, mainly expressed in the brain and reproductive system [24], which may explain why MKRN3 is highly expressed in Glioblastoma (Fig. S1C). Due to its low expression, its role in KIRC is limited.

The expression of MKRN1 and MKRN2 are highly correlated in KIRC tumor and normal tissues, and they show a direct interaction (Fig. 4). One possibility is that MKRN1 and MKRN2 can form a heterodimer, and they indeed synergistically promoted the mRNA expression of p53 (Figs 5D and S6B).

Both MKRN1 and MKRN2 were closely correlated with the expression of TP53 gene which encodes the most famous tumor suppressor p53 (Fig. 5A and B). KIRC cells that transfected with MKRN1 or MKRN2 showed higher p53 protein and mRNA levels (Figs 5C and S6A), and we inferred that MKRN1 or MKRN2 regulated the expression of p53 mainly at transcription level. As previously reported, MKRN2 binds selectively to the 3’UTR of a diverse subset of mRNAs [25], one possibility is that MKRN2 might binds to the 3’ UTR of TP53 mRNA or the regulators of TP53 mRNA. These are several studies seem to be contrary to our results: MKRN1 destabilize both p53 and p21 through ubiquitination and proteasome-dependent degradation that promotes the proliferation of H1299 and U2OS cells [4], and the ubiquitination of p53 by E3 ligase MKRN2 promotes melanoma cell proliferation [26]. Whether MKRN1/2 acts as an oncogene or tumor suppressor gene might depend on the cancer types, and the detailed reason need to be further explored.

The von Hippel-Lindau tumor suppressor gene (VHL) is mutated or silenced in $>$ 50% of sporadic clear cell renal cell carcinomas [27]. The VHL gene product, VHL protein, is the substrate recognition module of an E3 ubiquitin ligase that targets the hypoxia-inducible factor (HIF) for destruction in the presence of oxygen [27, 28], so VHL/HIFs pathway plays an important in KIRC development. Far from now, no published articles for the relationship between MKRNs genes and VHL/HIFs pathway were found, and this may be a good direction for future study.

Thus, our study suggests that MKRN1 and MKRN2 are tumor suppressor genes, and might be potential therapeutic targets for KIRC treatment. There are also some limitations or weaknesses in our study, such as no in vivo experiments or animal studies were not involved, and an ongoing effort is to test the functions of MKRN1 and MKRN2 in animal models and patient samples of KIRC.

5. Conclusions

Based on bioinformatics analysis and validated by experiments, MKRN1 and MKRN2 serve as tumor suppressors in kidney renal clear cell carcinoma, and act as promising therapeutic targets for KIRC treatment.

Data availability

All the data generated or analyzed during this study are included in the published article, and the cell lines and plasmids could be offered on reasonable request.

Ethics approval and consent to participate

No animals and patients’ samples were involved in this study.

Consent for publication

Not applicable.

Author contributions

Conception: Chuanyin Li, Chengcheng Wang and Yonghui Tao.

Interpretation or analysis of data: Chuanyin Li, Chengcheng Wang, Yun Yang, Yanyan Luo and Shuting Huang.

Preparation of the manuscript: Chuanyin Li, Chengcheng Wang, Yun Yang and Yanyan Luo.

Revision for important intellectual content: Chuanyin Li, Chengcheng Wang, Yonghui Tao and Yun Yang.

Supervision: Chengcheng Wang and Chuanyin Li.

Funding

This work was funded by the National Natural Science Foundation of China (31900804, 31960179), the Guizhou Province Science and Technology Foundation ([2020]1Y086), the Shanghai Sailing Program (19YF1444700), the Interdisciplinary Program of Shanghai Jiao Tong University (ZH2018QNA64) and the China Postdoctoral Science Foundation (2018M642109).

Supplementary data

The supplementary files are available to download from http://dx.doi.org/10.3233/CBM-210559.

sj-pdf-1-cbm-10.3233_CBM-210559.pdf - Supplemental material

Supplemental material, sj-pdf-1-cbm-10.3233_CBM-210559.pdf

sj-docx-1-cbm-10.3233_CBM-210559.docx - Supplemental material

Supplemental material, sj-docx-1-cbm-10.3233_CBM-210559.docx

Footnotes

Acknowledgments

We thank all members of our laboratory for their support.

Conflict of interest

All authors declare that they have no conflict of interest.

References

Abreu

A.P.

Macedo

D.B.

Brito

V.N.

Kaiser

U.B.

and Latronico

A.C.

, A new pathway in the control of the initiation of puberty: The MKRN3 gene, J Mol Endocrinol 54 (2015), R131–139.

Jong

M.T.

Gray

T.A.

Glenn

C.C.

Saitoh

Driscoll

D.J.

and Nicholls

R.D.

, A novel imprinted gene, encoding a RING zinc-finger protein, and overlapping antisense transcript in the Prader-Willi syndrome critical region, Hum Mol Genet 8 (1999), 783–793.

Kim

J.H.

Park

S.M.

Kang

M.R.

S.Y.

Lee

T.H.

Muller

M.T.

and Chung

I.K.

, Ubiquitin ligase MKRN1 modulates telomere length homeostasis through a proteolysis of hTERT, Genes Dev 19 (2005), 776–781.

Lee

E.W.

Lee

M.S.

Camus

Ghim

Yang

M.R.

N.C.

Lane

D.P.

and Song

, Differential regulation of p53 and p21 by MKRN1 E3 ligase controls cell cycle arrest and apoptosis, EMBO J 28 (2009), 2100–2113.

Omwancha

Zhou

X.F.

Chen

S.Y.

Baslan

Fisher

C.J.

Zheng

Cai

and Shemshedini

, Makorin RING finger protein 1 (MKRN1) has negative and positive effects on RNA polymerase II-dependent transcription, Endocrine 29 (2006), 363–373.

Lee

H.K.

Lee

E.W.

Seo

Jeong

Lee

S.H.

Kim

S.Y.

Jho

E.H.

Choi

C.H.

Chung

J.Y.

and Song

, Ubiquitylation and degradation of adenomatous polyposis coli by MKRN1 enhances Wnt/beta-catenin signaling, Oncogene 37 (2018), 4273–4286.

Lee

E.W.

Kim

J.H.

Ahn

Y.H.

Seo

Jeong

Kim

S.J.

J.Y.

Park

K.M.

Lee

H.W.

Park

E.J.

Chun

K.H.

and ong

, Ubiquitination and degradation of the FADD adaptor protein regulate death receptor-mediated apoptosis and necroptosis, Nat Commun 3 (2012), 978.

Shin

Ito

Ichikawa

Tokunaga

Sakata-Sogawa

and Tanaka

, MKRN2 is a novel ubiquitin E3 ligase for the p65 subunit of NF-kappaB and negatively regulates inflammatory responses, Sci Rep 7 (2017), 46097.

Jia

Tang

Yang

Yuan

Han

Fang

Huang

and Geng

, Ubiquitination of IGF2BP3 by E3 ligase MKRN2 regulates the proliferation and migration of human neuroblastoma SHSY5Y cells, Biochem Biophys Res Commun 529 (2020), 43–50.

10.

Jiang

Ren

Wudu

Wang

Song

Jiang

and Qiu

, MKRN2 inhibits migration and invasion of non-small-cell lung cancer by negatively regulating the PI3K/Akt pathway, J Exp Clin Cancer Res 37 (2018), 189.

11.

Zheng

Tang

Zang

Y.S.

Zeng

Liu

Shen

Pang

Wang

Xie

Luo

Jia

Huang

Wei

Huang

Chen

Zhu

Zhang

Xiao

Zhang

Fletcher

J.A.

and Wang

, E3 ligase MKRN3 is a tumor suppressor regulating PABPC1 ubiquitination in non-small cell lung cancer, J Exp Med 218 (2021).

12.

Zhang

Liu

Wang

and Qiu

, Elevated expression of MKRN3 in squamous cell carcinoma of the head and neck and its clinical significance, Cancer Cell Int 21 (2021), 557.

13.

Grignon

D.J.

and Che

, Clear cell renal cell carcinoma, Clin Lab Med 25 (2005), 305–316.

14.

Hsieh

J.J.

Purdue

M.P.

Signoretti

Swanton

Albiges

Schmidinger

Heng

D.Y.

Larkin

and Ficarra

, Renal cell carcinoma, Nat Rev Dis Primers 3 (2017), 17009.

15.

Leng

Pan

and Zhang

, Overexpression of IRF3 predicts poor prognosis in clear cell renal cell carcinoma, Int J Gen Med 14 (2021), 5675–5692.

16.

Chandrashekar

D.S.

Bashel

Balasubramanya

S.A.H.

Creighton

C.J.

Ponce-Rodriguez

Chakravarthi

and Varambally

, UALCAN: A portal for facilitating tumor subgroup gene expression and survival analyses, Neoplasia 19 (2017), 649–658.

17.

Tang

Kang

Gao

and Zhang

, GEPIA:A web server for cancer and normal gene expression profiling and interactive analyses, Nucleic Acids Res 45 (2017), W98–W102.

18.

Nagy

Munkacsy

and Gyorffy

, Pancancer survival analysis of cancer hallmark genes, Sci Rep 11 (2021), 6047.

19.

Gao

Aksoy

B.A.

Dogrusoz

Dresdner

Gross

Sumer

S.O.

Sun

Jacobsen

Sinha

Larsson

Cerami

Sander

and Schultz

, Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal, Sci Signal 6 (2013), pl1.

20.

Liu

Pan

Ding

Yang

Chen

Zhang

and Nie

, UBQLN4 is an ATM substrate that stabilizes the anti-apoptotic proteins BCL2A1 and BCL2L10 in mesothelioma, Mol Oncol 15 (2021), 3738–3752.

21.

Gao

Zhou

Qiu

Peng

Qiu

Liang

Gao

and Luo

, Total flavonoids from Semen Cuscutae target MMP9 and promote invasion of EVT cells via Notch/AKT/MAPK signaling pathways, Sci Rep 8 (2018), 17342.

22.

Tao

Gao

Zhang

Zou

Ruan

Wang

G.L.

and Hu

, A CRISPR-based approach for targeted DNA demethylation, Cell Discov 2 (2016), 16009.

23.

Gray

T.A.

Hernandez

Carey

A.H.

Schaldach

M.A.

Smithwick

M.J.

Rus

Marshall Graves

J.A.

Stewart

C.L.

and Nicholls

R.D.

, The ancient source of a distinct gene family encoding proteins featuring RING and C(3)H zinc-finger motifs with abundant expression in developing brain and nervous system, Genomics 66 (2000), 76–86.

24.

Han

Zhang

Guo

Yang

Peng

Tian

Wang

Zhou

Han

Wang

and Hu

, MKRN3-mediated ubiquitination of Poly(A)-binding proteins modulates the stability and translation of GNRH1 mRNA in mammalian puberty, Nucleic Acids Res 49 (2021), 3796–3813.

25.

Wolf

E.J.

Miles

Lee

E.S.

Nabeel-Shah

Greenblatt

J.F.

Palazzo

A.F.

Tropepe

and Emili

, MKRN2 Physically Interacts with GLE1 to Regulate mRNA Export and Zebrafish Retinal Development, Cell Rep 31 (2020), 107693.

26.

Zhang

Cui

and Zheng

, Ubiquitination of P53 by E3 ligase MKRN2 promotes melanoma cell proliferation, Oncol Lett 19 (2020), 1975–1984.

27.

Kaelin

W.G.

, Jr., The von Hippel-Lindau tumor suppressor gene and kidney cancer, Clin Cancer Res 10 (2004), 6290S–5S.

28.

Kaelin

W.G.

, Jr., Von Hippel-Lindau disease: Insights into oxygen sensing, protein degradation, and cancer, J Clin Invest 132 (2022).