Research Paper Volume 11, Issue 2 pp 615—633

MiR-223-3p promotes cell proliferation and metastasis by downregulating SLC4A4 in clear cell renal cell carcinoma

Wen Xiao1,2,3, , Xuegang Wang1,2,3, , Tao Wang1,2,3, , Jinchun Xing1,2,3, ,

  • 1 Department of Urology, The First Affiliated Hospital of Xiamen University, Xiamen, Fujian, China
  • 2 Center of Diagnosis and Treatment of Urinary System Diseases, The First Affiliated Hospital of Xiamen University, Xiamen, Fujian, China
  • 3 The Key Laboratory of Urinary Tract Tumors and Calculi of Xiamen City, The First Affiliated Hospital of Xiamen University, Xiamen, Fujian, China

Received: October 7, 2018       Accepted: January 5, 2019       Published: January 22, 2019      

https://doi.org/10.18632/aging.101763
How to Cite

Copyright: © 2019 Xiao et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

MicroRNAs (miRNAs) are known to affect the occurrence and progression of cancer. We therefore evaluated the involvement of miR-223-3p in renal cell cancer. MiR-223-3p was highly expressed in clear cell renal cell cancer tissues. Clear cell renal cell cancer patients with higher miR-223-3p expression had higher tumor stages and grades and poorer prognoses. In renal cancer cells, overexpression of miR-223-3p enhanced cell proliferation and metastasis, while inhibition of miR-223-3p reduced the malignant capacity of the cells. MiR-223-3p was found to bind directly to solute carrier family 4, member 4 (SLC4A4) mRNA, thereby reducing SLC4A4 mRNA and protein expression. SLC4A4 overexpression restrained cell proliferation and metastasis by suppressing Kirsten rat sarcoma viral oncogene (KRAS) expression in renal cancer cells. SLC4A4 expression correlated negatively with miR-223-3p expression in patient samples. Given that miR-223-3p suppressed the SLC4A4/KRAS axis, miR-223-3p gene therapy could be an effective treatment for renal cancer.

Introduction

Kidney cancer accounts for 3.77% of adult malignancies, and is a common lethal urological malignancy. In 2018, about 65,340 new kidney cancer cases and 14,970 deaths were estimated in the United States [1], and about 403,262 new kidney cancer cases (2.2% of all cancers) and 175,098 deaths (1.8% of all cancers) were reported worldwide [2]. Renal cell carcinoma (RCC) accounts for more than 90% of kidney cancer cases. The disease originates from the renal epithelium and encompasses more than 10 histological and molecular subtypes [3]. The incidence of RCC has increased gradually in the past two decades, and up to 30% of patients with a primary diagnosis of RCC will experience cancer metastasis [4, 5].

Clear cell RCC (ccRCC) is the most common RCC subtype, and has the highest rate of mortality, invasion and metastasis [6]. Given the extreme lethality of ccRCC, early diagnosis is critical for its treatment. Thus, it would be significant to determine the molecular mechanisms responsible for the genesis and development of ccRCC.

MicroRNAs (miRNAs), first discovered in 1993 [7], are small, non-coding RNA molecules that contribute to multiple cellular processes, including cancer cell proliferation [8, 9], apoptosis [10], metastasis [11, 12], radiation sensitivity [13] and chemosensitivity [14]. MiRNAs may function as oncogenes (oncomiRNAs) or tumor suppressors by downregulating the transcription of various genes [15, 16]. Mature miRNAs bind to the mRNA of target genes to degrade them and inhibit gene translation [1719]. Recently, several miRNAs have been identified as cancer biomarkers associated with disease progression, including miR-31-5p, miR-144-3p and miR-129-3p [2022].

MiR-223-3p is frequently upregulated in tumors from different organs and tissues, such as hepatocellular carcinoma [23], lung cancer [24], gastric cancer [25] and prostate cancer [26]. Nevertheless, to the best of our knowledge, the involvement of miR-223-3p in ccRCC has not yet been investigated. Therefore, in the present study, we systematically examined the function of miR-223-3p in ccRCC, including its expression in ccRCC tissues, its association with the pathological tumor (T) stage, grade and prognosis, and its effects on RCC cell proliferation, migration and invasion in vitro. We also determined its target mRNA molecule and its downstream effects on gene expression in RCC.

Results

MiR-223-3p is significantly upregulated and associated with a poor prognosis in ccRCC patients in TCGA-KIRC

Gene expression data from ccRCC patients were obtained from The Cancer Genome Atlas (TCGA) ccRCC dataset (TCGA-KIRC). Detailed clinicopathological information on these patients is presented in Table 1. The relative expression of miR-223-3p, shown in log2 (FPKM+1) form, ranged from 4.85 to 10.62 units in normal tissues and from 3.59 to 10.92 units in tumor tissues. MiR-223-3p expression was significantly higher in ccRCC tissues than in the corresponding non-cancerous tissues. Similar results were obtained from 71 paired ccRCC tissues and corresponding non-cancerous tissues (Figure 1A). Receiver operator characteristic (ROC) curve analysis demonstrated that miR-223-3p could sufficiently discriminate ccRCC tissues from normal tissues (area under the curve [AUC]: 0.7609, 95% confidence interval: 0.6964 - 0.8254) and from paired normal tissues (AUC: 0.7337, 95% confidence interval: 0.6494 - 0.8180) (Figure 1B). We also screened the database for a classic oncomiRNA (miR-21) and several miRNAs from recently published articles (miR-452, miR-543 and miR-708). We found that miR-223-3p was less sensitive than miR-21, but more sensitive than miR-708 (Supplementary Figure 1).

Table 1. Correlations of miR-223-3p and SLC4A4 levels in tissues and clinicopathological parameters of ccRCC patients

Correlation of miR-223-3p and SLC4A4 expression with clinicopathological parameters.

miR-223-3p expressionSLC4A4 mRNA expression
ParameterLow (n=250)High (n=250)P valueLow (n=250)High (n=250)P value
Age(years)0.2830.795
<=60132119130121
>60118131120129
Sex0.2590.090
Male170157173154
female80937796
T stage0.0410.354
T1+T2169146152163
T3+T4811049887
N stage0.4170.173
N0+ NX245241240246
N159104
M stage0.0130.137
M0+ MX222201205218
M128494532
Grade0.0020.472
G1+G213196109118
G3+G4119154141132
TNM stage0.0180.202
I+II162135141156
III+IV8811510994
MiR-223-3p expression is upregulated in ccRCC and predicts a poor prognosis

Figure 1. MiR-223-3p expression is upregulated in ccRCC and predicts a poor prognosis. MiR-223-3p levels in 71 normal tissues and 506 ccRCC tissues were downloaded from the dataset of TCGA-KIRC. The levels of miR-223-3p were compared according to different clinicopathological parameters: (A) cancer versus para-cancer, and cancer versus paired para-cancer. (B) The ROC curve displays that miR-223-3p could effectively distinguish ccRCC from para-cancer and paired para-cancer tissues, with AUCs of 0.7609, CI 95% (0.6964 to 0.8254) (p < 0.0001) and 0.7337, CI 95% (0.6494 to 0.8180) (p < 0.0001), respectively. (CI) MiR-223-3p expression was analyzed in ccRCC patients according to (C) lymph node metastasis, (D) distant metastasis, (E) T stage, (F) TNM stage, (G) grade, (H) OS and (I) DFS. Data are shown as the mean ± SEM. * p < 0.05; ** p < 0.01; *** p < 0.001; ### p < 0.001.

T-test analysis demonstrated that miR-223-3p expression was higher in patients with lymph node metastases than in those without (Figure 1C). High miR-223-3p expression was also observed in patients with distant metastases (Figure 1D). MiR-223-3p expression was significantly higher in T stages III and IV than in T stages I and II, and was greater in T stage IV than in T stage III (Figure 1E). Higher miR-223-3p levels were associated with more advanced pathological TNM stages (T stage, lymphatic metastasis and distant metastasis grouping) and grades in ccRCC patients (Figure 1F and 1G).

Kaplan-Meier survival analysis was used to determine the correlation of miR-223-3p expression with the overall survival (OS) and disease-free survival (DFS) of ccRCC patients. Patients from TCGA-KIRC were divided into two groups based on the median miR-223-3p level. Patients with higher miR-223-3p expression (> 6.7 units) exhibited shorter OS (Figure 1H, p < 0.0001) and DFS (Figure 1I, p = 0.0346). The expression of miR-223-3p was then compared between patients with longer and shorter OS and DFS (OS-Good and OS-Poor, DFS-Good and DFS-Poor). High miR-223-3p expression indicated a poor prognosis (Supplementary Figure 2). Univariate and multivariate survival analyses demonstrated that miR-223-3p expression was an independent prognostic factor for OS in ccRCC patients (p = 0.000) (Table 2).

Table 2. Univariate and multivariate analyses of miR-223-3p and SLC4A4 mRNA level and patient overall survival

Univariate analysisMultivariate analysisc
VariableHRa95%CIbPHRa95% CIbP
Overall survival (n = 500)
Age (years)
>60 (n = 249)1.8041.322-2.4620.000*1.6831.231-2.3000.001*
≤60 (n = 251)
Sex
Male (n = 327)0.9170.671-1.2540.588
Female (n = 173)
T stage
T3 or T4 (n = 185)2.9302.155-3.9830.000*1.6281.137-2.3330.008*
T1 or T2 (n = 315)
N stage
N1 (n = 14)3.5001.842-6.6480.000*
N0 or NX (n =486)
M stage
M1 (n = 77)4.2013.069-5.7490.000*2.7061.886-3.8830.000*
M0 or MX (n = 423)
Grade
G3 or G4 (n = 273)2.6021.840-3.6790.000*1.5841.092-2.2980.015*
G1 or G2 (n = 227)
miR-223-3p
High (n = 250)2.1301.553-2.9210.000*1.8881.371-2.5990.000*
Low (n = 250)
SLC4A4
High (n = 250)0.5580.409-0.7610.000*0.6420.470-0.8780.005*
Low (n = 250)

These data indicate that miR-223-3p expression is upregulated and is significantly associated with various clinicopathological parameters (such as T stage, lymphatic metastasis and distant metastasis) in ccRCC. High miR-223-3p expression predicted a poor prognosis and thus could be a potential prognostic indicator for ccRCC.

MiR-223-3p promotes renal cancer cell proliferation, migration and invasion in vitro

We next evaluated miR-223-3p levels in ccRCC patients from The First Affiliated Hospital of Xiamen University. The results confirmed our observations in TCGA-KIRC (Figure 2A). We also determined the miR-223-3p levels of three renal cancer cell lines (A498, Caki-1 and ACHN) [27], and found that ACHN cells had the lowest expression, while Caki-1 cells had the highest (Figure 2B).

MiR-223-3p promotes RCC cell proliferation, migration and invasion in vitro

Figure 2. MiR-223-3p promotes RCC cell proliferation, migration and invasion in vitro. MiR-223-3p levels were measured in (A) renal cancer tissue samples and (B) cell lines. The transfection effects of (C) miR-223-3p mimics and (D) miR-223-3p inhibitors were measured in 786-O and Caki-1 cells. The effects of (E) miR-223-3p mimics and (F) miR-223-3p inhibitors on cell proliferation were measured by the CCK-8 assay in 786-O and Caki-1 cells. Representative pictures display the (G) migration and (H) invasion of 786-O cells treated with miR-223-3p inhibitors. Data are shown as the mean ± SEM. * p < 0.05; ** p < 0.01; *** p < 0.001.

Then, we used gain- or loss-of-function assays to evaluate the effects of miR-223-3p on the proliferation, migration and invasion of 786-O (RCC) cells and Caki-1 cells. The effects of miR-223-3p mimics and inhibitors are shown in Figure 2C and 2D. MiR-223-3p overexpression significantly enhanced the proliferation of 786-O and Caki-1 cells (Figure 2E), while miR-223-3p knockdown significantly reduced their proliferation (Figure 2F). Knockdown of miR-223-3p significantly repressed cell migration and invasion compared to the negative control (NC) in Transwell assays (Figure 2G and 2H). These data indicate that miR-223-3p promotes the proliferation, migration and invasion of renal cancer cells.

Bioinformatic analysis of miR-223-3p target genes in ccRCC

To identify the potential targets of miR-223-3p, we used two prediction software programs (TargetScan and miRDB) to score predicted molecules. We then selected the top 20 molecules scored by each program, and identified five molecules of interest (Figure 3A): solute carrier family 4, member 4 (SLC4A4), Ras homolog family member B (RHOB), interleukin 6 signal transducer (IL6ST), LIM domain only 2 (LMO2) and F-box and WD repeat domain containing 7 (FBXW7). To narrow the field of prospective candidates, we performed a correlation analysis in TCGA-KIRC. A heat map revealed the correlations between miR-223-3p and the five target molecules (Figure 3B). As shown in Figure 3C, only SLC4A4 and RHOB levels exhibited significant negative correlations with miR-223-3p levels. Then, we compared the mRNA levels of the five target molecules between cancerous and non-cancerous tissues in TCGA-KIRC. SLC4A4 expression was significantly lower in ccRCC tissues than in the corresponding non-cancerous tissues, while RHOB expression was higher in ccRCC tissues (Figure 3D). In addition, miR-223-3p expression exhibited a significant negative correlation with SLC4A4 expression in ccRCC patient samples from our university (Figure 3E).

Bioinformatic analysis of miR-223-3p target genes in ccRCC

Figure 3. Bioinformatic analysis of miR-223-3p target genes in ccRCC. (A) Bioinformatic prediction of the top 20 mRNA targets of miR-223-3p in TargetScan and miRDB. (B) Heat map depicting the expression of miR-223-3p and five target genes in samples from TCGA-KIRC. (C) Correlation analysis of the expression of miR-223-3p and five target genes in cancer samples from TCGA-KIRC. (D) Relative mRNA expression of five target genes in cancer samples from TCGA-KIRC. (E) A qRT-PCR analysis demonstrated the negative correlation between SLC4A4 and miR-223-3p expression in ccRCC tissues (R = -0.437, p = 0.037). Data are shown as the mean ± SEM. * p < 0.05; ** p < 0.01; *** p < 0.001.

MiR-223-3p directly binds to SLC4A4

To determine whether miR-223-3p directly binds to SLC4A4, we evaluated the changes in protein and mRNA levels by Western blotting and qRT-PCR, respectively, following cellular treatment with miR-223-3p mimics and inhibitors. In 786-O and Caki-1 cells, miR-223-3p mimics reduced both the protein and mRNA levels of SLC4A4 (Figure 4A and 4B), while miR-223-3p inhibitors had the opposite effects (Figure 4C and 4D).

SLC4A4 is a direct target of miR-223-3p

Figure 4. SLC4A4 is a direct target of miR-223-3p. (A) Western blotting and (B) qRT-PCR analysis of SLC4A4 expression in 786-O and Caki-1 cells transfected with miR-223-3p mimics versus the corresponding NC. (C) Western blotting and (D) qRT-PCR analysis of SLC4A4 expression in 786-O and Caki-1 cells transfected with miR-223-3p inhibitors versus the corresponding NC. (E) The predicted binding sites for miR-223-3p in the SLC4A4 3′-UTR. The red nucleotides are the seed-pairing target sites of miR-223-3p. (F) Luciferase reporter assays demonstrate that the reporter activity of 786-O and Caki-1 cells decreased by approximately 50% upon co-transfection of the wild-type SLC4A4 3′-UTR reporter construct and miR-223-3p mimics. Data are shown as the mean ± SEM. * p < 0.05; ** p < 0.01; *** p < 0.001.

Then, we determined the effect of miR-223-3p on the 3′-untranslated region (UTR) of SLC4A4 using a luciferase reporter assay. Luciferase reporter constructs containing either wild-type or mutated SLC4A4 binding sequences upstream of the firefly luciferase gene were generated (Figure 4E). Caki-1 and 786-O cells were co-transfected with the reporter vectors and mimics or mimic controls. Luciferase activity was significantly reduced after miR-223-3p mimic co-transfection with WT vectors (Figure 4F). These results suggest that SLC4A4 is a direct target of miR-223-3p.

SLC4A4 is significantly downregulated and associated with a poor prognosis in ccRCC patients in TCGA-KIRC

As SLC4A4 was found to be a direct target of miR-223-3p, we investigated SLC4A4 mRNA levels in TCGA-KIRC. The relative expression of SLC4A4 in log2 (FPKM+1) form ranged from 4.85 to 10.62 units in normal tissues and from 3.59 to 10.92 units in tumor tissues. The expression of SLC4A4 was significantly lower in ccRCC tissues than in non-cancerous tissues (Figure 5A). To confirm the results from TCGA-KIRC, we examined three additional datasets in the Oncomine database (Figure 5B). Low SLC4A4 expression was detected in patients with distant metastases (Figure 5C). SLC4A4 expression was significantly lower in T stage IV than in T stages I, II and III (Figure 5D). Lower SLC4A4 levels were associated with more advanced pathological TNM stages and grades in ccRCC patients (Figure 5E and 5F). Patients with lower SLC4A4 expression exhibited shorter OS (Figure 5G, t-test, p < 0.0001) and DFS (Figure 5H, t-test, p = 0.005). Univariate and multivariate survival analyses indicated that SLC4A4 expression was an independent prognostic factor for OS and DFS in ccRCC patients (Tables 2 and 3).

SLC4A4 expression is downregulated in ccRCC and predicts a poor prognosis

Figure 5. SLC4A4 expression is downregulated in ccRCC and predicts a poor prognosis. SLC4A4 mRNA levels in 72 normal tissues and 533 ccRCC tissues were downloaded from the dataset of TCGA-KIRC. (A) SLC4A4 mRNA levels were lower in cancer tissues than in para-cancer tissues. (B) SLC4A4 levels in three additional ccRCC datasets. (CH) SLC4A4 levels were compared in ccRCC patients according to the following clinicopathological parameters: (C) distant metastasis, (D) T stage, (E) TNM stage, (F) grade, (G) OS and (H) DFS. Data are shown as the mean ± SEM. * p < 0.05; ** p < 0.01; *** p < 0.001; ## p < 0.01; ### p < 0.001.

Table 3. Univariate and multivariate analyses of miR-223-3p and SLC4A4 mRNA level and patient disease–free survival

Univariate analysisMultivariate analysisc
VariableHRa 95%CIbPHRa95% CIbP
Disease–free survival (n = 407)
Age (years)
>60 (n = 184)1.4581.020-2.0830.038*
≤60 (n = 223)
Sex
Male (n = 272)1.3450.904-2.0010.143
Female (n = 135)
T stage
T3 or T4 (n = 140)4.3332.989-6.2800.000*2.1161.387-3.2260.001*
T1 or T2 (n = 267)
N stage
N1 (n = 12)5.8482.935-11.6560.000*2.8501.396-5.8180.004*
N0 or NX (n =395)
M stage
M1 (n = 52)8.1245.579-11.8310.000*4.9013.222-7.4550.000*
M0 or MX (n = 355)
Grade
G3 or G4 (n = 210)3.2052.120-4.8440.000*2.1951.425-3.3820.000*
G1 or G2 (n = 197)
miR-223-3p
High (n = 250)1.3420.938-1.9210.108
Low (n = 250)
SLC4A4
High (n = 250)0.5860.407-0.8440.004*0.6070.419-0.8800.008*
Low (n = 250)

SLC4A4 is downregulated and involved in the biological pathogenesis of ccRCC

As SLC4A4 was downregulated in ccRCC patients in TCGA-KIRC and was an independent prognostic factor for ccRCC, we tested SLC4A4 expression in ccRCC patient samples from our hospital by quantitative real-time polymerase chain reaction (qRT-PCR) and immunohistochemistry analyses. SLC4A4 was downregulated in ccRCC tissues compared to normal tissues (Figure 6A and 6B). Immunohistochemistry data from the Human Protein Atlas also demonstrated that SLC4A4 was downregulated in ccRCC tissues compared to normal tissues (Figure 6C). Gene set enrichment analysis demonstrated that SLC4A4 expression was associated with the epithelial-mesenchymal transition and Kirsten rat sarcoma viral oncogene (KRAS) signaling (Figure 6D).

SLC4A4 overexpression inhibits cell proliferation, migration and invasion in vitro

Figure 6. SLC4A4 overexpression inhibits cell proliferation, migration and invasion in vitro. SLC4A4 was downregulated in ccRCC tissues, as shown by (A) qRT-PCR and (B) immunohistochemistry analysis. (C) SLC4A4 was downregulated in ccRCC tissues, as shown by immunohistochemistry analysis from the Human Protein Atlas. (D) Gene set enrichment analysis comparing low and high SLC4A4 expression groups in the database of TCGA. Enrichment curves are shown for activated gene sets related to the epithelial-mesenchymal transition and KRAS signaling. (E) and (F) SLC4A4 mRNA and protein levels were successfully upregulated in 786-O and Caki-1 cells, and SLC4A4 overexpression inhibited KRAS expression. (GI) SLC4A4 overexpression inhibited (G) proliferation, (H) migration and (I) invasion in Caki-1 cells. Data are shown as the mean ± SEM. * p < 0.05; ** p < 0.01; *** p < 0.001.

Then, we evaluated the function of SLC4A4 in ccRCC by overexpressing SLC4A4 in RCC cell lines. The mRNA and protein levels of SLC4A4 were confirmed to be upregulated in 786-O and Caki-1 cells transfected with plasmid vectors expressing SLC4A4 (Figure 6E and 6F). Overexpression of SLC4A4 significantly repressed the proliferation of 786-O cells compared to the NC group (Figure 6G). Furthermore, SLC4A4 vector transfection markedly inhibited the migration and invasion of Caki-1 cells in Transwell assays (Figure 6H and 6I). These results demonstrate that SLC4A4 may function as a tumor suppressor in ccRCC cells.

SLC4A4 attenuates the function of miR-223-3p in renal cancer cells

To determine whether SLC4A4 could reverse the oncogenic effects of miR-223-3p, we co-transfected renal cancer cells with miR-223-3p mimics and SLC4A4. The mRNA levels of miR-223-3p and SLC4A4 are shown in Figure 7A and 7B, while the protein levels of SLC4A4 and KRAS are shown in Figure 7C. Overexpression of SLC4A4 impaired the ability of miR-223-3p mimics to induce KRAS expression (Figure 7C) and to enhance 786-O and Caki-1 cell proliferation (Figure 7D). The enhancement of Caki-1 cell migration and invasion by miR-223-3p was also reversed by co-transfection with SLC4A4 (Figure 7E and 7F). These data demonstrate that SLC4A4 partially alleviates the oncogenic effects of miR-223-3p in RCC.

SLC4A4 overexpression attenuates the effects of miR-223-3p mimics on cellular processes

Figure 7. SLC4A4 overexpression attenuates the effects of miR-223-3p mimics on cellular processes. (A) MiR-223-3p and (B) SLC4A4 levels were measured by qRT-PCR in 786-O and Caki-1 cells. Overexpression of SLC4A4 impaired the effects of miR-223-3p mimics on (C) KRAS expression, (D) cell proliferation, (E) migration and (F) invasion in Caki-1 cells. Data are shown as the mean ± SEM. * p < 0.05; ** p < 0.01; *** p < 0.001.

Discussion

Aberrant miRNA expression is known to be associated with cancer cell proliferation, invasion and metastasis, and miR-223-3p expression is dysregulated in many tumors. Some studies have revealed that miR-223-3p functions as an anti-tumor gene by inhibiting the metastasis and progression of osteosarcoma [28] and glioblastoma [29]. However, other studies have indicated that miR-223-3p is an oncomiRNA that is upregulated in lung cancer, gastric carcinoma and prostate cancer tissues compared with the corresponding non-cancerous tissues [2426]. However, the expression and function of miR-223-3p in ccRCC remained unknown.

In this study, miR-223-3p was found to be significantly upregulated and associated with various clinicopathological parameters in ccRCC. Higher miR-223-3p expression was associated with worse OS and DFS in ccRCC patients. Moreover, miR-223-3p expression was greater in renal cancer cells than in normal renal cells. Overexpression of miR-223-3p promoted the proliferation and motility of ccRCC cells. Thus, we identified miR-223-3p as a key oncomiRNA in renal cancer.

Then, we searched for candidate genes using two publicly available algorithms (miRDB and TargetScan), and identified five potential target genes of miR-223-3p. Only SLC4A4 and RHOB expression correlated negatively with miR-223-3p expression. However, in TCGA-KIRC, SLC4A4 expression was significantly lower in ccRCC tissues than in non-cancerous tissues, while RHOB expression was higher in ccRCC tissues. Moreover, clinical data analysis revealed a negative association between miR-223-3p and SLC4A4 levels in human ccRCC samples.

The kidney expresses the solute carrier (SLC) transporters at higher levels than any other organ [30]. SLC-family proteins have different effects on cancer development, depending on the substrates they transport [3133]. SLC4A4 has been reported to promote colon and breast cancer growth and migration. SLC4A4 mRNA expression depends exclusively on hypoxia inducible factor 1 subunit alpha [34], a tumor suppressor gene in ccRCC [35]. However, the function of SLC4A4 and the miRNAs regulating its expression in ccRCC have not been previously investigated.

We analyzed the expression and prognostic effects of SLC4A4 in the database of TCGA-KIRC. SLC4A4 was significantly downregulated in ccRCC tissues, and low SLC4A4 expression was associated with a poor prognosis. Gene set enrichment analysis demonstrated that the epithelial-mesenchymal transition and KRAS signaling were significantly enriched in response to low SLC4A4 expression in patients with ccRCC. SLC4A4 overexpression impaired RCC cell proliferation, migration and invasion in vitro. Furthermore, forced expression of SLC4A4 reversed the effects of miR-223-3p on cell growth and migration in Cell Counting kit-8 (CCK-8) and Transwell assays by suppressing KRAS signaling.

In the present study, miR-223-3p and SLC4A4 were identified as important independent biomarkers that could be used to predict the clinical outcomes of ccRCC. MiR-223-3p promoted cell invasion, migration, growth and proliferation in RCC by directly binding to SLC4A4. This is the first report indicating that SLC4A4 is a potential target gene of miR-223-3p in ccRCC. These results suggest that the miR-223-3p/SLC4A4 axis could be an ideal prognostic predictor and therapeutic candidate for human ccRCC.

Materials and Methods

Patient samples

Surgical specimens were obtained from 24 patients between 2017 and 2018 at the Department of Urology of The First Affiliated Hospital of Xiamen University. Resected tissues were frozen in liquid nitrogen and stored at -80°C for further experiments. Informed consent was obtained from all patients, and the experimental procedures were approved by the Institutional Review Board of Xiamen University.

RNA extraction and qRT-PCR

Tissue RNA was extracted with the TRIzol reagent (Thermo, Massachusetts, USA) according to the manufacturer’s instructions. The purity and concentration of each RNA sample was tested with a NanoDrop 2000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). A qPCR analysis was performed according to the manufacturer’s instructions (LightCycler 480II; Roche, Basel, Switzerland) with SYBR Green mix (Thermo, Massachusetts, USA). The relative levels of miR-223-3p and SLC4A4 were calculated by: 2-ΔCt (ΔCt = Ctgene–Ctnormalizer). MiRNA primers were purchased from RiboBio (Guangzhou, China). The primers for SLC4A4 and GAPDH were purchased from GENEWIZ (GENEWIZ, Suzhou, China): SLC4A4, Forward 5’-TTCACGGAACTGGATGAGCT-3’, Reverse 5’-ACTGTGGGAGAGAAGAAGCC-3’; GAPDH, Forward 5’-GAGTCAACGGATTTGGTCGT-3’, Reverse 5’-GACAAGCTTCCCGTTCTCAG-3’.

Cell culture and transient transfection

HK2, 786-O, ACHN and Caki-1 cells were purchased from The American Type Culture Collection (USA). Cancer cells were cultivated in high-glucose Dulbecco’s modified Eagle’s medium (DMEM; Wuhan Boster Biological Technology, Ltd., Wuhan, China) containing 10% fetal bovine serum (Gibco; Thermo Fisher Scientific, Inc., Waltham, MA, USA) at 37°C in a 5% CO2 incubator, as previously reported [36]. Human miR-223-3p mimics, miR-223-3p inhibitors and the corresponding controls were synthesized by RiboBio. Plasmid vectors expressing SLC4A4 or the NC were constructed by Genechem (Shanghai, China). Lipofectamine 2000 reagent (Invitrogen, Carlsbad, CA, USA) was used to transfect cells with miR-223-3p mimics, miR-223-3p inhibitors, the SLC4A4 plasmid and the corresponding controls in six-well plates according to the recommendations of the manufacturer.

Cell proliferation analysis

Caki-1 and 786-O cells were first transfected with miR-223-3p mimics, miR-223-3p inhibitors, SLC4A4 or the NC. Then, the cells were added to a 96-well plate at a density of 3x103 cells/well. The cell proliferation rate (optical density value) was determined with a CCK-8 assay (Dojindo Molecular Technologies, Inc, Rockville, MD, USA) according to the manufacturer’s protocol. Two hours after the addition of 10 µL of CCK-8 solution to each well, the optical density value was measured at 450 nm. Three independent experiments were conducted for each assay.

Migration and invasion assays

Migration assays were performed as previously described [37]. Cells were homogenized in the absence of plasma for 24 h, and then cultured at 1x105 cells/well in 24-well Transwell plates with polycarbonate membrane inserts (Corning Inc., Corning, NY, USA). For the invasion assays, the membrane was coated with Matrigel (Thermo Fisher Scientific), and cells were cultured at 2x105 cells/well. After incubation for 24 hours, the cells were fixed with 100% methanol and stained with 0.05% crystal violet. Five random fields were counted, and three independent experiments were conducted for each assay.

Luciferase assays

Wild-type and mutant SLC4A4 3′-UTR reporters were purchased from RiboBio. Tumor cells transfected with miR-223-3p mimics or miR-NC in 24-well plates were co-transfected with 500 ng of a luciferase reporter with Lipofectamine 2000 reagent (Thermo Fisher Scientific). All dual luciferase assays were performed with a Dual Luciferase Assay kit (Promega, Madison, WI, USA). Normalized Renilla-luciferase values were calculated relative to the control, according to the manufacturer’s protocol.

Western blotting

Tissues and cells were pyrolysed in a protein lysis system containing radioimmunoprecipitation assay buffer (Wuhan Boster Biological Technology, Ltd.), a protease inhibitor cocktail (Roche Diagnostics, Indianapolis, IN, USA) and PMSF (Wuhan Boster Biological Technology, Ltd.). Protein concentrations were measured with a bicinchoninic acid kit (Beyotime Institute of Biotechnology) at 562 nm. Thirty micrograms of each protein were subjected to sodium dodecyl sulfate polyacrylamide gel electrophoresis and transferred to polyvinylidene fluoride membranes (EMD Millipore, Bedford, MA, USA) for 90 min. After the transfer, the membranes were blocked in phosphate-buffered saline with 5% non-fat milk for one hour and then incubated with antibodies against SLC4A4 (1:1000; A5332; ABclonal Biotech Co., Ltd.), KRAS (1:1000; A11059; ABclonal Biotech Co., Ltd.) and GAPDH (1:2000; BM3876; Wuhan Boster Biological Technology, Ltd.) at 4°C overnight. The next day, the membranes were washed and then incubated with secondary antibodies (1:5000; BA1020; Wuhan Boster Biological Technology, Ltd.) at room temperature (25°C) for two hours. Finally, the membranes were washed, and protein levels were detected with a ChemiDoc-XRS+ (Bio-Rad Laboratories, Inc., Hercules, CA, USA).

Immunohistochemistry

Adjacent normal tissues and ccRCC tissues were fixed in formalin, dehydrated and embedded. Tissue sections were then incubated with a primary rabbit SLC4A4 polyclonal antibody (1:100) at 4°C overnight. The sections were washed three times with phosphate-buffered saline and then incubated with a goat anti-rabbit secondary antibody (1:200; GB23303; Servicebio, Inc., Woburn, MA, USA) for two hours at room temperature.

Bioinformatic analysis

RNA-seq data from ccRCC patients (TCGA-KIRC) and the patients’ clinical information were obtained from the Xena Functional Genomics Explorer (https://xenabrowser.net/heatmap/) of the University of California Santa Cruz [38] and the Oncomine database (https://www.oncomine.org). Detailed clinical information on staging can be found at (https://www.cancer.org/cancer/kidney-cancer/detection-diagnosis-staging/staging.html). Protein immunohistochemistry data were obtained from the Human Protein Atlas (https://www.proteinatlas.org/). Candidate miR-223-3p target genes were identified through two publicly available algorithms, miRDB (http://mirdb.org/miRDB) and TargetScan (http://www.targetscan.org). Gene set enrichment analysis was used to determine the SLC4A4 pathway involved in the pathogenesis of ccRCC in patients from TCGA-KIRC (http://www.broadinstitute.org/gsea) [39]. Following the performance of 1,000 permutations, p values < 0.05 and a false discovery rate < 25% were required for the enriched gene set analysis to be considered significantly enriched.

Statistical analysis

The RNA levels of paired samples were analyzed with a paired-sample t-test, while those of unpaired samples were analyzed by one-way ANOVA for four groups (#) or by the t-test for two groups (*). The AUC and ROC curve were used to distinguish the clinical classifications. A Kaplan-Meier curve was used to evaluate the survival rate according to miR-223-3p and SLC4A4 expression with a log-rank test. Long and short OS and DFS were discriminated based on good (≥ 5 years alive or disease-free) and poor (death or recurrence or progression in ≤ 2 years) prognoses [40]. The prognostic significance of miR-223-3p and SLC4A4 in ccRCC were analyzed by univariate and multivariate Cox proportional hazard regression analyses. P values < 0.05 were considered to be statistically significant. All statistical analyses were conducted with SPSS Statistics 22.0 (IBM SPSS, Chicago, IL, USA).

Supplementary Materials

Supplementary Figures

Author Contributions

Jinchun Xing and Wen Xiao designed the study. Wen Xiao and Xuegang Wang performed the experiments. Wen Xiao drafted the manuscript. Xuegang Wang helped to draft the manuscript. Tao Wang and Jinchun Xing conducted the statistical analyses. Jinchun Xing revised the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We wish to thank all our colleagues in the Department of Urology, The First Affiliated Hospital of Xiamen University.

Conflicts of Interest

The authors declare no conflicts of interest pertaining to this manuscript.

Funding

This study was supported by a grant from the National Natural Scientific Foundation of China (Grant no. 81602218).

References

  • 1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018.CA: a cancer journal for clinicians.2018;68:7-30. https://doi.org/10.3322/caac.21442
  • 2. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries.CA: a cancer journal for clinicians.2018; https://doi.org/10.3322/caac.21492 [PubMed]
  • 3. Hsieh JJ, Purdue MP, Signoretti S, Swanton C, Albiges L, Schmidinger M, Heng DY, Larkin J, Ficarra V. Renal cell carcinoma.Nature reviews Disease primers.2017;3:17009 https://doi.org/10.1038/nrdp.2017.9 [PubMed]
  • 4. Linehan WM. The genetic basis of kidney cancer: implications for management and use of targeted therapeutic approaches.European urology.2012;61:896-898. https://doi.org/10.1016/j.eururo.2012.02.022 [PubMed]
  • 5. Bhatt JR and Finelli A. Landmarks in the diagnosis and treatment of renal cell carcinoma.Nature reviews Urology.2014;11:517-525. https://doi.org/10.1038/nrurol.2014.194 [PubMed]
  • 6. Rini BI, Campbell SC, Escudier B. Renal cell carcinoma.Lancet.2009;373:1119-1132. https://doi.org/10.1016/S0140-6736(09)60229-4 [PubMed]
  • 7. Lee RC, Feinbaum RL, Ambros V. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14.Cell.1993;75:843-854. https://doi.org/10.1016/0092-8674(93)90529-Y [PubMed]
  • 8. Hiyoshi Y, Kamohara H, Karashima R, Sato N, Imamura Y, Nagai Y, Yoshida N, Toyama E, Hayashi N, Watanabe M, Baba H. MicroRNA-21 regulates the proliferation and invasion in esophageal squamous cell carcinoma.Clinical cancer research : an official journal of the American Association for Cancer Research.2009;15:1915-1922. https://doi.org/10.1158/1078-0432.CCR-08-2545 [PubMed]
  • 9. Gao Y, Ma X, Yao Y, Li H, Fan Y, Zhang Y, Zhao C, Wang L, Ma M, Lei Z, Zhang X. miR-155 regulates the proliferation and invasion of clear cell renal cell carcinoma cells by targeting E2F2.Oncotarget.2016;7:20324-20337. https://doi.org/10.18632/oncotarget.7951 [PubMed]
  • 10. Chen H, Liu C, Jiang H, Gao Y, Xu M, Wang J, Liu S, Fu Y, Sun X, Xu J, Zhang J, Dai L. Regulatory Role of miRNA-375 in Expression of BMP15/GDF9 Receptors and its Effect on Proliferation and Apoptosis of Bovine Cumulus Cells.Cellular physiology and biochemistry : international journal of experimental cellular physiology, biochemistry, and pharmacology.2017;41:439-450. https://doi.org/10.1159/000456597 [PubMed]
  • 11. Ma L, Teruya-Feldstein J, Weinberg RA. Tumour invasion and metastasis initiated by microRNA-10b in breast cancer.Nature.2007;449:682-688. https://doi.org/10.1038/nature06174 [PubMed]
  • 12. Zhou Q, Zheng X, Chen L, Xu B, Yang X, Jiang J, Wu C. Smad2/3/4 Pathway Contributes to TGF-beta-Induced MiRNA-181b Expression to Promote Gastric Cancer Metastasis by Targeting Timp3.Cellular physiology and biochemistry : international journal of experimental cellular physiology, biochemistry, and pharmacology.2016;39:453-466. https://doi.org/10.1159/000445638 [PubMed]
  • 13. Huang X, Taeb S, Jahangiri S, Emmenegger U, Tran E, Bruce J, Mesci A, Korpela E, Vesprini D, Wong CS, Bristow RG, Liu FF, Liu SK. miRNA-95 mediates radioresistance in tumors by targeting the sphingolipid phosphatase SGPP1.Cancer research.2013;73:6972-6986. https://doi.org/10.1158/0008-5472.CAN-13-1657 [PubMed]
  • 14. Chen J, Xu Y, Tao L, Pan Y, Zhang K, Wang R, Chen LB, Chu X. MiRNA-26a Contributes to the Acquisition of Malignant Behaviors of Doctaxel-Resistant Lung Adenocarcinoma Cells through Targeting EZH2.Cellular physiology and biochemistry : international journal of experimental cellular physiology, biochemistry, and pharmacology.2017;41:583-597. https://doi.org/10.1159/000457879 [PubMed]
  • 15. Lu YC, Cheng AJ, Lee LY, You GR, Li YL, Chen HY, Chang JT. MiR-520b as a novel molecular target for suppressing stemness phenotype of head-neck cancer by inhibiting CD44.Scientific reports.2017;7:2042 https://doi.org/10.1038/s41598-017-02058-8 [PubMed]
  • 16. Rupaimoole R and Slack FJ. MicroRNA therapeutics: towards a new era for the management of cancer and other diseases.Nature reviews Drug discovery.2017;16:203-222. https://doi.org/10.1038/nrd.2016.246 [PubMed]
  • 17. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function.Cell.2004;116:281-297. https://doi.org/10.1016/S0092-8674(04)00045-5 [PubMed]
  • 18. Krol J, Loedige I, Filipowicz W. The widespread regulation of microRNA biogenesis, function and decay.Nature reviews Genetics.2010;11:597-610. https://doi.org/10.1038/nrg2843 [PubMed]
  • 19. Ha M and Kim VN. Regulation of microRNA biogenesis.Nature reviews Molecular cell biology.2014;15:509-524. https://doi.org/10.1038/nrm3838 [PubMed]
  • 20. Bi J, Liu H, Cai Z, Dong W, Jiang N, Yang M, Huang J, Lin T. Circ-BPTF promotes bladder cancer progression and recurrence through the miR-31-5p/RAB27A axis.Aging.2018;10:1964-1976. https://doi.org/10.18632/aging.101520 [PubMed]
  • 21. Lou N, Ruan AM, Qiu B, Bao L, Xu YC, Zhao Y, Sun RL, Zhang ST, Xu GH, Ruan HL, Yuan CF, Han WW, Shi HC, et al. miR-144-3p as a novel plasma diagnostic biomarker for clear cell renal cell carcinoma.Urologic oncology.2017;35:36 e37-36 e14. https://doi.org/10.1016/j.urolonc.2016.07.012 [PubMed]
  • 22. Chen X, Ruan A, Wang X, Han W, Wang R, Lou N, Ruan H, Qiu B, Yang H, Zhang X. miR-129-3p, as a diagnostic and prognostic biomarker for renal cell carcinoma, attenuates cell migration and invasion via downregulating multiple metastasis-related genes.Journal of cancer research and clinical oncology.2014;140:1295-1304. https://doi.org/10.1007/s00432-014-1690-7 [PubMed]
  • 23. Oksuz Z, Serin MS, Kaplan E, Dogen A, Tezcan S, Aslan G, Emekdas G, Sezgin O, Altintas E, Tiftik EN. Serum microRNAs; miR-30c-5p, miR-223-3p, miR-302c-3p and miR-17-5p could be used as novel non-invasive biomarkers for HCV-positive cirrhosis and hepatocellular carcinoma.Molecular biology reports.2015;42:713-720. https://doi.org/10.1007/s11033-014-3819-9 [PubMed]
  • 24. Liu C, Yang Z, Deng Z, Zhou Y, Gong Q, Zhao R, Chen T. Upregulated lncRNA ADAMTS9-AS2 suppresses progression of lung cancer through inhibition of miR-223-3p and promotion of TGFBR3.IUBMB life.2018;70:536-546. https://doi.org/10.1002/iub.1752 [PubMed]
  • 25. Yang F, Xu Y, Liu C, Ma C, Zou S, Xu X, Jia J, Liu Z. NF-kappaB/miR-223-3p/ARID1A axis is involved in Helicobacter pylori CagA-induced gastric carcinogenesis and progression.Cell death & disease.2018;9:12 https://doi.org/10.1038/s41419-017-0020-9 [PubMed]
  • 26. Wei Y, Yang J, Yi L, Wang Y, Dong Z, Liu Z, Ou-yang S, Wu H, Zhong Z, Yin Z, Zhou K, Gao Y, Yan B, et al. MiR-223-3p targeting SEPT6 promotes the biological behavior of prostate cancer.Scientific reports.2014;4:7546 https://doi.org/10.1038/srep07546 [PubMed]
  • 27. Brodaczewska KK, Szczylik C, Fiedorowicz M, Porta C, Czarnecka AM. Choosing the right cell line for renal cell cancer research.Molecular cancer.2016;15:83 https://doi.org/10.1186/s12943-016-0565-8 [PubMed]
  • 28. Ji Q, Xu X, Song Q, Xu Y, Tai Y, Goodman SB, Bi W, Xu M, Jiao S, Maloney WJ, Wang Y. miR-223-3p Inhibits Human Osteosarcoma Metastasis and Progression by Directly Targeting CDH6.Molecular therapy : the journal of the American Society of Gene Therapy.2018;26:1299-1312. https://doi.org/10.1016/j.ymthe.2018.03.009 [PubMed]
  • 29. Ding Q, Shen L, Nie X, Lu B, Pan X, Su Z, Yan A, Yan R, Zhou Y, Li L, Xu J. MiR-223-3p overexpression inhibits cell proliferation and migration by regulating inflammation-associated cytokines in glioblastomas.Pathology, research and practice.2018;214:1330-1339. https://doi.org/10.1016/j.prp.2018.05.012 [PubMed]
  • 30. Zhang Y, Zhang Y, Sun K, Meng Z, Chen L. The SLC transporter in nutrient and metabolic sensing, regulation, and drug development.Journal of molecular cell biology.2018; https://doi.org/10.1093/jmcb/mjy052 [PubMed]
  • 31. Bhutia YD, Babu E, Ramachandran S, Yang S, Thangaraju M, Ganapathy V. SLC transporters as a novel class of tumour suppressors: identity, function and molecular mechanisms.The Biochemical journal.2016;473:1113-1124. https://doi.org/10.1042/BJ20150751 [PubMed]
  • 32. Wu C, Li D, Jia W, Hu Z, Zhou Y, Yu D, Tong T, Wang M, Lin D, Qiao Y, Zhou Y, Chang J, Zhai K, et al. Genome-wide association study identifies common variants in SLC39A6 associated with length of survival in esophageal squamous-cell carcinoma.Nature genetics.2013;45:632-638. https://doi.org/10.1038/ng.2638 [PubMed]
  • 33. Chen L, Zhou S, Qin J, Hu H, Ma H, Liu B, Wang X, Ma J, Ye S, Zhong C, Zhou G, Liang C. Combination of SLC administration and Tregs depletion is an attractive strategy for targeting hepatocellular carcinoma.Molecular cancer.2013;12:153 https://doi.org/10.1186/1476-4598-12-153 [PubMed]
  • 34. Parks SK and Pouyssegur J. The Na(+)/HCO3(-) Co-Transporter SLC4A4 Plays a Role in Growth and Migration of Colon and Breast Cancer Cells.Journal of cellular physiology.2015;230:1954-1963. https://doi.org/10.1002/jcp.24930 [PubMed]
  • 35. Shen C, Beroukhim R, Schumacher SE, Zhou J, Chang M, Signoretti S, Kaelin WG Jr. Genetic and functional studies implicate HIF1alpha as a 14q kidney cancer suppressor gene.Cancer discovery.2011;1:222-235. https://doi.org/10.1158/2159-8290.CD-11-0098 [PubMed]
  • 36. Xiao W, Lou N, Ruan H, Bao L, Xiong Z, Yuan C, Tong J, Xu G, Zhou Y, Qu Y, Hu W, Gao Y, Ru Z, et al. Mir-144-3p Promotes Cell Proliferation, Metastasis, Sunitinib Resistance in Clear Cell Renal Cell Carcinoma by Downregulating ARID1A.Cellular physiology and biochemistry : international journal of experimental cellular physiology, biochemistry, and pharmacology.2017;43:2420-2433. https://doi.org/10.1159/000484395 [PubMed]
  • 37. Liu R, Qin X, Ji C, Zeng W, Yang Y, Tan W. Pygopus 2 promotes kidney cancer OS-RC-2 cells proliferation and invasion in vitro and in vivo.Asian J Urol.2015;2:151-157. https://doi.org/10.1016/j.ajur.2015.06.009 [PubMed]
  • 38. Goldman M, Craft B, Swatloski T, Cline M, Morozova O, Diekhans M, Haussler D, Zhu J. The UCSC Cancer Genomics Browser: update 2015.Nucleic acids research.2015;43:D812-817. https://doi.org/10.1093/nar/gku1073 [PubMed]
  • 39. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles.Proceedings of the National Academy of Sciences of the United States of America.2005;102:15545-15550. https://doi.org/10.1073/pnas.0506580102 [PubMed]
  • 40. Zheng J, Wang L, Peng Z, Yang Y, Feng D, He J. Low level of PDZ domain containing 1 (PDZK1) predicts poor clinical outcome in patients with clear cell renal cell carcinoma.EBioMedicine.2017;15:62-72. https://doi.org/10.1016/j.ebiom.2016.12.003 [PubMed]