Downregulation of CLDN7 due to promoter hypermethylation is associated with human clear cell renal cell carcinoma progression and poor prognosis
Journal of Experimental & Clinical Cancer Research volume 37, Article number: 276 (2018)
Metastasis is the primary cause of death in renal cell carcinoma (RCC). Loss of cell-to-cell adhesion, including tight junctions (TJs) is the initial step in the process of metastasis. Claudin-7 (CLDN7) is a major component of TJs. However, the clinical significance and its regulation of kidney tumorigenesis remain poorly understood.
A total of 120 fresh clear cell RCC (ccRCC) specimens and 144 primary RCC and adjacent nonmalignant renal paraffin specimens were obtained from Department of Urology, Peking University First Hospital. Expression of CLDN7 in ccRCC tissues and cell lines were determined using bioinformatic data mining, quantitative real-time PCR (qRT-PCR), Western blotting and immunostaining. The clinical significance of CLDN7 expression and promoter DNA methylation status was analyzed in ccRCC patients from Peking University First Hospital and The Cancer Genome Atlas. Additionally, the methylation specific-PCR, bisulfite genomic sequencing and demethylation analysis of CLDN7 were performed. Biological functions of CLDN7 were investigated by examining cell proliferation using MTS assays and EdU incorporation assays, cell migration by in vitro wound healing assays and transwell migration assays, cell invasion by transwell invasion assays, and cell apoptosis by flow cytometry. Mouse model experiments were performed to confirm the effects of CLDN7 on tumor growth and metastasis in vivo. The molecular mechanism of CLDN7 function was investigated using gene-set enrichment analysis (GSEA) and high-throughput cDNA sequencing (RNA-Seq) and confirmed by qRT-PCR, Western blot and immunostaining in vitro and in vivo.
Our findings revealed that CLDN7 is frequently downregulated via hypermethylation of its promoter in ccRCC. CLDN7 can help predict aggressive tumor status and poor prognosis in ccRCC patients. Interestingly, hypermethylation of the CLDN7 promoter was related to advanced ccRCC status and poor prognosis. Moreover, overexpression of CLDN7 induced cell apoptosis, suppressed proliferation, migration and invasion abilities of ccRCC cells both in vitro and in vivo. Additionally, GSEA and RNA-Seq results showed that CLDN7 had negative effects in cancer-associated signaling pathways and (epithelial-mesenchymal transition) EMT-related pathways. These results were validated by qRT-PCR, Western blot and immunostaining.
We have demonstrated a previously undescribed role of CLDN7 as a ccRCC suppressor and suggest that loss of CLDN7 potentiates EMT and tumor progression. CLDN7 may serve as a functional tumor suppressor in tumor progression and a potential biomarker and target in patients with ccRCC.
Renal cell carcinoma (RCC) accounts for up to 85% to 90% of all kidney cancers [1, 2]. According to the most recent pathological classification by the International Society of Urological Pathology (ISUP), RCC mainly includes clear cell (ccRCC), papillary (pRCC) and chromophobe (chRCC) subtypes , with ccRCC being the most common subtype . The 5-year survival rate of patients with kidney cancer increased approximately 30% in recent years. This improvement is largely due to low-stage and low-grade tumors being incidentally detected by improved early-detection techniques [2, 5, 6]. However, one third of patients with kidney cancer present with metastatic disease , and 20%~ 30% of the patients who undergo curative surgery relapse in distant sites during follow-up . Although there has been considerable progress in the systemic treatment of metastatic renal cell carcinoma in the past 20 years , metastatic RCC remains an incurable condition for the majority of patients.
Metastasis is the primary cause of death for RCC. The initial step in the process of metastasis is the loss of cell-to-cell adhesion in the neoplastic epithelium . Tight junctions (TJs) are the most apical components of the epithelial cell junctional complex and provide a form of cell–cell adhesion . Claudin-7 (CLDN7) is a major component of TJs in epithelial cells. In RCC, CLDN7 is a marker of renal tumors originating from the distal nephron marker, such as chRCC and oncocytoma . The expression of CLDN7 is significantly lower in ccRCC than non-ccRCC [10,11,12]. However, the clinical significance and molecular mechanisms of downregulation of CLDN7 in ccRCC remain unknown.
In this study, we confirmed that downregulation of CLDN7 due to hypermethylation may help predict aggressive tumor status and poor prognosis in ccRCC patients. Using combined data from human patients and in vitro and in vivo analyses in ccRCC cells manipulated for CLDN7 overexpression, we previously demonstrated a tumor suppressing role of CLDN7 and that it induced epithelial features in ccRCC cells. Bioinformatics and RNA sequencing analysis showed that genes that were influenced by CLDN7 overexpression were mostly enriched in cancer- and EMT-related pathways, indicating that CLDN7 may have an important role in the development and progression of ccRCC.
For DNA and RNA extraction, a total of 120 patients who had been pathologically diagnosed with ccRCC were included, and detailed clinicopathological features of the samples are summarized in Additional file 1: Table S1. A total of 120 paired ccRCC tissues and adjacent normal kidney tissues were immediately snap-frozen in liquid nitrogen following surgical resection. For immunohistochemical (IHC) staining, 144 patients were selected who had RCC and had undergone radical or partial nephrectomy between May 1, 2012 and August 30, 2012 at Peking University First Hospital. All the RCC paraffin specimens were pathologically diagnosed as RCC, including 129 ccRCC, 7 chRCC, 7 pRCC and 1 collecting duct RCC. This study was approved by ethics committee of the Peking University First Hospital (Beijing, China). Written informed consent was also obtained from all patients.
Bioinformatic data mining
We downloaded TCGA RCC RNA-Seq gene expression data, TCGA ccRCC Methylation 450 K data and clinical data from UCSC Xena (http://xena.ucsc.edu/). A total of 129 normal kidney, 534 ccRCC, 66 chRCC and 291 pRCC tissues had CLDN7 mRNA expression, and 319 ccRCC tissues had CLDN7 promoter methylation data. All 534 ccRCC tumors had clinical data for clinical correlation and survival analyses. The heat map and the correlation between CLDN7 expression and methylation status were further analyzed in the same patient cohort and verified using UCSC Xena data. The relationship between CLDN7 mRNA expression and promoter methylation status was searched in the ccRCC cohort in TCGA database using cBioportal (http://cbioportal.org). Additionally, gene-set enrichment analysis (GSEA) was performed to compare differences in molecular pathways in cell processes between the low CLDN7 and high CLDN7 groups on the data from the ccRCC dataset of TCGA.
Cell culture and transfection
The cell lines (HEK-293, Caki-2, Caki-1, OS-RC-2, 786-O, 769-P, ACHN and A498) were obtained from the American Type Culture Collection (Rockville, MD, USA). Cell lines were cultured according to conditions specified by the provider. The kidney cancer cell line KETR-3 was purchased from KeyGEN BioTECH Co., Ltd. (Jiangsu, China). For overexpression of CLDN7, recombinant pGC-LV-GV287-GFP vectors with the CLDN7 mRNA (NM_001185022.1) or with a scrambled control sequence (Control) were constructed by the Genechem Company (Genechem Co. Ltd., Shanghai, China). All the viral vectors contained GFP as a marker to track lentivirus-mediated target gene expression by fluorescence microscopy. Briefly, pGC-LV-GV287-GFP-CLDN7 or a control vector mixed with pHelper1.0 and pHelper2.0 were cotransfected into HEK-293 T cells with PEI (Sigma-Aldrich). Caki-1 and A498 cells were infected via lentiviruses according to the MOI value (the number of lentiviruses per number of cells) recommended by the manufacturer. Caki-1 or A498 CLDN7 and Caki-1 or A498 CTRL cells were sorted using a flow cytometer sorter (BD FACS AriaTM SORP, New Jersey, USA). The Luciferase-pcDNA3 was a gift from William Kaelin (Addgene plasmid # 18964) . Lentiviruses were produced using a three-vector system: Luciferase-pcDNA3: viral packaging (psPAX2): viral envelope (pMD2G) at a 4:3:1 ratio in HEK-293 T cells. Caki-1 CLDN7-Luc and Caki-1 CTRL-Luc cells were selected with neomycin (200 μg/mL). The transfection efficiency was validated (Additional file 2: Figure S1).
Quantitative real-time PCR (qRT-PCR) and reverse transcription PCR (RT-PCR)
Total RNA was extracted from the tissue samples or the transfected cells using the TRIzol reagent (Invitrogen; Thermo Fisher Scientific Inc.), according to the manufacturer’s instructions. cDNA was generated using reverse transcription (TansGEN, Beijing, China). qRT-PCR was performed using the ABI PRISM 7000 Fluorescent Quantitative PCR System (Applied Biosystems, Foster City, CA, USA), according to the manufacturer’s instructions, and normalized to GAPDH. RT-PCR was performed by electrophoresis on a 1.5% agarose gel. All experiments were repeated at least three times. The detailed primer sequences included in this study are shown in Additional file 3: Table S2.
Immunohistochemistry (IHC) and Western blot analysis
The immunohistochemistry (IHC) and IHC scoring were carried out according to protocols that have been described previously . Protein lysates were prepared by homogenization in 1% NP40 containing 1 mM PMSF and 20 μg protein was separated by SDS-PAGE. The immunoreactive bands were visualized using an Immobilon™ Western Kit (Millipore, Billerica, MA) using the SYNGENE G: BOX imaging system (Frederick, USA). Antibodies specific to CLDN7 (ab27487), BCL-2 (ab32124), PARP1 (ab32064), Caspase-3 (ab13847), E-cadherin (CDH1, ab76055), N-cadherin (CDH2, ab18203), Vimentin (ab92547) and TGFB1 (ab25121) were purchased from Abcam (Hong Kong, China). GAPDH (TA309157) and Ki-67 (TA500265) were purchased from ZSGB-BIO, Beijing, China. Cleaved-Caspase3 (Asp175) (9661S) was purchased from Cell Signaling Technology, MA, USA. The antibodies for the IHC and Western blot assays were diluted with phosphate buffered solution (PBS) or PBS plus Tween-20 (PBST), according to the manufacturer’s instructions.
Methylation-specific PCR (MSP) and bisulfite genomic sequencing PCR (BGS)
Genomic DNA (1 μg) was denatured using 0.3 M NaOH for 10 min at 37 °C. The samples were then incubated at 50 °C for 16 h after adding hydroquinone (Sigma-Aldrich, St. Louis, Missouri, USA) and sodium bisulfate (Sigma-Aldrich). Genomic DNA was analyzed via MSP using primer sets located within a CpG-rich area in the CLDN7 promoter. PCR samples were then resolved by electrophoresis on a 1.5% agarose gel. For the BGS assay, DNA was purified, and a CpG-rich promoter region was amplified by PCR. The PCR products were purified and cloned into a PCR 2.1-TA cloning vector (Invitrogen). A minimum of six positive clones from each product were selected for sequencing. The detailed primer sequences are shown in Additional file 3: Table S2.
Caki-1 and A498 cells were seeded in six-well plates at a concentration of 1 × 105 cells per well, grown for 24 h, and then treated with 5 μM 5-Aza-2′-deoxycytidine (5-Aza-dC, A, Sigma-Aldrich) for 4 days. Cells were cultured with or without 100 nM Trichostatin A (TSA, T, Sigma-Aldrich) for the final 24 h. RNA was isolated for RT-PCR analysis and DNA was extracted for CLDN7 MSP.
The metabolic activity of Caki-1 cells was assessed using CellTiter 96™ AQueous Nonradioactive Cell Proliferation Assay, according to the manufacturer’s instructions (Promega, Madison, WI, USA). The optical density of the wells was measured at 450 nm using a Multiscan microplate spectrophotometer (Thermo LabSystems, Milford, MA, USA).
Ethynyl-2-deoxyuridine (EdU) incorporation assay
Cell proliferation was also determined by Ethynyl-2-deoxyuridine incorporation assay using an EdU Apollo DNA in vitro kit (RIBOBIO, Guangzhou, China), according to the manufacturer’s instructions. Experiments were repeated at least three times.
Flow cytometry analysis assay
Cell apoptosis was assayed by staining with Annexin V-APC and PI (KeyGEN BioTECH) following manufacturer’s instructions and detected by a flow cytometer (FACSCalibur, Becton Dickinson, New Jersey, USA).
Wound healing assays
Cell migration was determined via a wound-healing assay. Briefly, approximately 3 × 105 cells were seeded in 12-well plates at equal densities and grown to 80%~ 90% confluency. Artificial gaps were generated by a 200 μL sterile pipette tip after transfection. Wounded areas were marked and photographed with a microscope (Leica DM IL, Leica Microsystems, Germany) equipped with a digital camera (Leica DFC300FX).
Transwell migratory and invasive assays
For the transwell migration assay, 1000 cells were plated into the upper chambers (24-well insert, pore size 8 μm, Corning) with 100 μL serum-free PRIM-1640. The lower chambers were filled with 500 μL PRIM-1640 containing 10% fetal bovine serum. Two days later, cells under the surface of the lower chamber were washed with PBS and stained with 0.5% crystal violet for 30 min.
For the invasion assay, 2000 cells were seeded on transwells (24-well insert, pore size 8 μm, Corning) coated with 60 μL Matrigel (1:3 dilution in PBS, Product #354234, Corning Inc., NY, USA). The culture conditions were the same as described for the transwell migration assay. After 72 h, adherent cells on the lower surface were stained with 0.5% crystal violet. The number of cells on the lower surface was photographed with a microscope (Leica DM IL, Leica Microsystems) equipped with a digital camera (Leica DFC300FX).
Mouse model experiments
Animal experiments were conducted in accordance with the NIH Guidelines for the Care and Use of Laboratory Animals with the approval of the Review Board of Peking University First Hospital, Beijing. Mice were maintained under pathogen-free conditions with regulated temperature and humidity levels. Mice were randomly assigned to cages in groups of 5 and fed ad libitum under controlled light/dark cycles.
Twenty 5-week-old male BALB/c nude mice were purchased from Vitalriver, Beijing, China. Approximately 5 × 106 Caki-1 CLDN7 or Caki-1 CTRL cells suspended in 100 μL Hank’s Balanced Salt Solution (Thermo Fisher Scientific Inc., USA) were mixed with Matrigel (1:1, Product #354234, Corning Inc., NY, USA). Then, 200 μL tumor cells were subcutaneously implanted on the right flank of 6-week BALB/c nude mice using a 28-gauge needle (Thermo Fisher Scientific Inc., USA). Tumor size was measured every fourth day and calculated using the formula: (length × width2)/2.
For the metastasis experiment, eight 5-week-old male B-NDG mice (NOD- Prkdcscid IL2rgtm1/Bcgen) that lacked mature T cells, B cells, and natural killer (NK) cells, were purchased from BIOCYTOGEN, Beijing, China. Approximately 1 × 106 Caki-1 CTRL-Luc or Caki-1 CLDN7-Luc cells were suspended in 200 μL PBS and injected into the lateral tail veins of each unanesthetized B-NDG mouse at six-weeks-old. Twenty days after injection, mice were anesthetized with isoflurane (YIPIN Pharmaceutical CO., LTD, Hebei, China). Ten minutes after D-Luciferin, sodium salt (150 mg/kg) was injected intraperitoneally, and cancer cells were detected with an in vivo imaging system, Xenogen IVIS (PerkinElmer, MA, USA). The total flux in photons per second were calculated and reported for each mouse’s lung and liver region using Living Image 4.3.1 (PerkinElmer/Caliper).
High-throughput cDNA sequencing (RNA-Seq)
The RNA-Seq experiments were performed by Novogene (Beijing, China). The RNA-seq library was prepared for sequencing using standard Illumina protocols. Briefly, total RNAs from Caki-1 CTRL and CLDN7 cells were isolated using TRIzol reagent (Invitrogen) and treated with RNase-free DNase I (New England Biolabs, MA, USA), to remove any contaminating genomic DNA. RNA extraction was performed using Dynabeads oligo(dT) (Invitrogen Dynal). Double-stranded complementary DNAs were synthesized from 1 μg of total RNA using Superscript II reverse transcriptase (Invitrogen) and random hexamer primers. Escherichia coli RNase H (New England Biolabs) were added to remove RNA complementary to the cDNA. The cDNA was then fragmented by nebulization, and the standard Illumina protocol was followed thereafter to create the mRNA-seq library. The libraries were sequenced on an Illumina HiSeq 2000 platform. Sequencing reads were aligned to the human genome (hg19) using the TopHat program (v2.1.1) set to the default parameters. Total read counts for each protein-coding gene were extracted using HTSeq (version 0.6.0) and then loaded into R package DESeq2 to calculate the differentially expressed genes with a cut-off fold change of ≥1.5 and an FDR < 0.05. Gene expression was calculated using the RPKM (reads per kilobase transcriptome per million reads) method. Experiments were repeated three times.
All data were expressed as means ± SD from at least three separate experiments. All data were analyzed using SPSS 20.0 statistical software (IBM. Chicago, IL, USA). The comparison of immunostaining intensities between the two groups was analyzed using a Mann-Whitney U test. The clinicopathological correlation with gene expression of the patients was analyzed using a chi-square test. Survival curves for patients were plotted using the Kaplan-Meier method, with log-rank tests for statistical significance. Uni- and multivariable Cox regression analyses were used to test the prognostic relevance between clinicopathological and immunohistochemical data. In case of multiple tests, a one-way ANOVA followed by Bonferroni-Holm procedure was applied. Other statistically significant differences were determined using Student’s T-test. Statistical significance was defined as a two-tailed p < 0.05. *p < 0.05, **p < 0.01.
CLDN7 downregulation was associated with poor prognostic features
Our previous microarray chip assay of 6 ccRCC tissues and paired adjacent normal kidney tissues showed that CLDN7 was significantly downregulated in ccRCC tissue compared with normal kidney tissue (Fig. 1A a). To determine the expression of CLDN7, we analyzed the RNA-Seq data of RCC cases from TCGA. CLDN7 mRNA expression was decreased significantly in ccRCC tissues, but not in chRCC and pRCC tissue, compared to adjacent normal kidney tissues (Fig. 1A b, c, d). Next, we performed qRT-PCR in the RCC cell lines and 12 paired human ccRCC tissues. CLDN7 mRNA expression was significantly lower in the RCC cells and ccRCC tissues than in the human embryonic kidney cells (HEK-293) and adjacent normal kidney tissues (Fig. 1b). Western blots produced similar results, with CLDN7 protein being downregulated in the RCC cell lines (Fig. 1c). Immunostaining of the 144 RCC tissues and adjacent normal kidney tissues also showed that intensity of immunostaining was significantly weaker in the ccRCC tissues compared with their adjacent normal tissues (Fig. 1d). Moreover, CLDN7 staining intensity was lower in ccRCC tissue than in non-ccRCC tissue. The representative immunostaining of CLDN7 in normal kidney tissues, ccRCC, pRCC and chRCC tissues are shown in Fig. 1e.
To explore the correlation between CLDN7 expression and clinicopathological features, we divided the TCGA 534 ccRCC patients into different groups based on the mean value (10.625 FPKM) of CLDN7 mRNA expression of the adjacent normal kidney tissues. The results showed that CLDN7 mRNA level was strongly correlated with gender, pathologic M, pathologic T, pathologic stage and overall survival (Additional file 4: Table S3).
We further analyzed the CLDN7 immunostaining intensity with clinicopathological parameters in 129 ccRCC patients from Peking University First Hospital. The results showed that the CLDN7 protein level was strongly correlated with age, pathologic T, histologic grade and risk of recurrence, metastasis and death after surgery (Table 1). Collectively, CLDN7 was downregulated in ccRCC patients and was significantly associated with tumor progression features.
To explore the potential prognostic significance of CLDN7, we evaluated 534 ccRCC patients from TCGA using a Kaplan-Meier analysis with a log-rank test. As shown in Fig. 1f, the low CLDN7 mRNA expression group had a poorer prognosis than the higher CLDN7 group. Moreover, in 144 RCC patients from Peking University First Hospital, immunostaining results confirmed that ccRCC patients with low CLDN7 protein expression levels experienced recurrent or metastatic tumors earlier after surgery, compared to those with higher expression levels (Fig. 1g). Univariate and multivariate Cox regression analyses indicated that higher ages, pathologic T stage and elevated creatinine levels after surgery were unfavorable prognostic factors in ccRCC patients, except for high protein expression of CLDN7, which was a favorable prognostic factor (HR = 0.227, 95% CI = 0.093–0.551, p = 0.001) (Table 2). Therefore, decreased CLDN7 expression indicates poor prognosis in ccRCC patients.
CLDN7 downregulation by promotor methylation in ccRCC was associated with adverse pathologic results and poor prognosis
Epigenetic alterations, such as promoter CpG methylation, could mediate the activation of oncogene and inactivate tumor suppressor genes, thus contributing to tumorigenesis . CLDN7 promoter hypermethylation has been reported in several cancers [16,17,18,19]. In the TCGA ccRCC dataset, we generated a comparable heatmap in which we found that when some regions of the CLDN7 promoter were hypermethylated, CLDN7 expression was lower (Fig. 2a). Moreover, we observed that CLDN7 methylation status was negatively related to its mRNA and protein expression in the TCGA ccRCC dataset (Fig. 2b). Collectively, CLDN7 downregulation in ccRCC was significantly associated with promoter hypermethylation.
Furthermore, we analyzed the 14 methylation regions of CLDN7 in the TCGA datasets (UCSC Xena). The results suggested that most of the methylation regions of CLDN7 (8/14) were significantly hypermethylated in ccRCC tissues compared with adjacent normal tissues (Fig. 2c and Additional file 5: Figure S2). Among the four most hypermethylated CLDN7 promoter regions (cg00072720, cg11941546, cg14034852 and cg05490983), we found that the methylation status of cg00072720 was negatively correlated with CLDN7 mRNA expression and had the largest negative correlation coefficient (Fig. 2d and Additional file 6: Figure S3). Taken together, we chose to use cg00072720 in the following experiments.
To explore the clinical application of CLDN7 methylation status, we performed methylation-specific PCR analyses in 108 ccRCC tissues from Peking University First Hospital. The relationship between CLDN7 methylation status and clinicopathologic features are summarized in Table 3. Hypermethylated CLDN7 was significantly correlated with advanced age, pathologic T and histologic G, which agrees with our bioinformatic analysis of the TCGA ccRCC Methylation 450 K dataset (n = 319, Additional file 7: Table S4). Interestingly, the CLDN7 promoter hypermethylated group had a poor prognosis (Additional file 8: Figure S4). Therefore, hypermethylation of CLDN7 promoter was associated with downregulation of CLDN7 and a poor prognosis in ccRCC patients.
Promoter hypermethylation contributed to downregulation of CLDN7 in ccRCC
As shown in Fig. 2, promoter hypermethylation was associated with low CLDN7 expression in ccRCC patients. We designed primers for the cg00072720 region that was mapped to the CLDN7 promoter for experimental validation. The results suggest that RCC cells had low mRNA levels of CLDN7, which is in accordance with their high hypermethylation (Figs. 1b and 3a). The CLDN7 promoter was hypermethylated in both Caki-1 and A498 cells. Treatment with 5-Aza-2′-deoxycytidine (5-aza-dC, 5-Aza, A) and Trichostatin A (TSA, T) induced demethylation of the CLDN7 promoter and increased CLDN7 mRNA expression (Fig. 3b). CLDN7 hypermethylation was also found in 34/108 (31.5%) ccRCC tissues (Fig. 3c and Table 3). Since MSP analyzed only a few CpG sites in the promoter, we sequenced the CLDN7 promoter region described above. Sequencing of sodium-bisulfite-treated DNA from two ccRCC patients (No. 83 and 85) showed hypermethylation of CpG sites in ccRCC tumor tissues (T) compared with adjacent normal kidney tissues (N) (Fig. 3d).
CLDN7 overexpression attenuated tumor proliferation and induced tumor apoptosis in metastatic ccRCC cells
Lower expression and higher methylation of CLDN7 were observed in Caki-1 and A498 cells. Therefore, we generated a GFP labeled, lentiviral vector-expressing control and CLDN7-expressing Caki-1 and A498 cells for a functional study. The effect of CLDN7 on the growth of Caki-1 and A498 cells was determined using MTS assays and an EdU incorporation assay. The results suggested that overexpression of CLDN7 inhibited cell proliferation in Caki-1 and A498 cells (Fig. 4a, b).
Cell proliferative alteration is usually caused by changes in the cell cycle or apoptosis. Flow cytometry was performed to evaluate the early apoptosis in the cells. Overexpression of CLDN7 induced significant apoptosis in Caki-1 and A498 cells (Fig. 4c). Consistent with the weak proliferative ability and high apoptotic rate observed in the CLDN7 overexpressed Caki-1 and A498 cells, a xenograft experiment in mice found that CLDN7 overexpressed tumors grew slower than those in a control group (Fig. 4d).
CLDN7 overexpression significantly decreased cell migratory and invasive ability
The migratory and invasive abilities of CLDN7 overexpressed Caki-1 and A498 cells were significantly decreased, compared with control cells. The results showed that migratory distances of the CLDN7 overexpressed group were significantly inhibited (Fig. 5a). A transwell migratory and invasive assay also found that migratory and invasive cells were clearly attenuated in Caki-1 and A498 cells that overexpressed CLDN7, compared with control cells (Fig. 5b). To confirm the inhibitory effect of CLDN7 on cell migration and invasion, an in vivo model of metastasis was established. Metastasis experiments indicated that overexpression of CLDN7 inhibited metastasis to liver and lung (Fig. 5c). The number and size of metastatic lesions formed by Caki-1 CLDN7-Luc cells were significantly fewer and smaller than the metastatic lesions formed by Caki-1 Control-Luc cells (Fig. 5c).
Mechanisms of tumor suppression by CLDN7 in ccRCC
To understand the tumor suppressive role of CLDN7 in ccRCC, gene-set enrichment analysis on the data from the database of TCGA was used to compare the differences in cell processes between the low- and high- CLDN7 groups. Interestingly, the results suggest that several pathways relating to cancer and TGF-β, WNT, Focal adhesion and NOTCH pathways relating to epithelial-mesenchymal transition (EMT) and tumor progression were all decreased in ccRCC with high CLDN7 mRNA expression (Additional file 9: Figure S5 and Additional file 10: Table S5).
To confirm these results, we performed mRNA sequencing in Caki-1 CLDN7 and control cells, which suggested that many genes are differently expressed when CLDN7 is overexpressed. These differently expressed genes could separate CLDN7 overexpressed Caki-1 cells from control cells (Fig. 6a). Furthermore, a Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis found that genes influenced by CLDN7 overexpression were mostly enriched in pathways involved with cancer and cell adhesion molecules (CAMs), extracellular matrix (ECM)-receptor and Hedgehog signaling pathways that relate to EMT and tumor progression (Fig. 6b). These genes screened by RNA-Seq were validated by qRT-PCR. The results showed that oncogenes (BCL2, HIF1A, GLI-1, ITGB-1 and AR) were significantly downregulated, while tumor suppressor gene p21 was significantly upregulated in CLDN7 overexpressed Caki-1 cells. Additionally, the epithelial marker, E-cadherin, was increased while mesenchymal markers, N-cadherin and Vimentin were decreased. This finding suggests that there is a mesenchymal-epithelial transition (MET) in CLDN7 overexpressed cells. Interestingly, the EMT inducers, TGFB1 and TWIST1 were both downregulated by CLDN7 (Fig. 6c). Cell pro-apoptotic markers, such as cleaved-PARP1 and cleaved-Caspase3, were upregulated by CLDN7 in Caki-1 and A498 cells, and the anti-apoptotic marker BCL2 was downregulated (Fig. 6d). Moreover, BCL2 and the cell proliferation marker Ki-67 were suppressed by CLDN7 while cleaved-Caspase3 was upregulated in xenograft tissues (Fig. 6d). To validate the mechanism of CLDN7 in cell adhesion and metastasis, EMT markers were evaluated in cultured cells and cell derived xenografts. It was found that CLDN7 overexpression promoted the expression of E-cadherin in vitro and in vivo. In contrast, expression of TGFB1, N-cadherin and Vimentin were all downregulated by CLDN7 (Fig. 6e). Taken together, CLDN7 suppressed cell growth and metastasis by inducing cell apoptosis and inhibiting EMT.
Previous studies have reported that CLDN7 is dysregulated in various cancers [10,11,12, 18,19,20,21,22,23]. During our study, we confirmed that there is a downregulation of CLDN7 in ccRCC. Bioinformatic Data Mining found that DNA hypermethylation in the promoter of CLDN7, and MSP and BGS results confirmed the CLDN7 promoter hypermethylation in ccRCC tissues. Analysis of follow up data, clinical features and CLDN7 expression and methylation data, demonstrated that the lower expression and higher methylation status of CLDN7 were significantly associated with tumor progression and poor prognosis. Interestingly, in vitro and in vivo assays both found that CLDN7 overexpression significantly inhibited cell proliferation and induced apoptosis. Cell migratory and invasive abilities were also suppressed by CLDN7. To explore the molecular mechanism of the tumor suppressive function of CLDN7 in ccRCC, GSEA was performed to evaluate the different gene expression profiles between low- and high-CLDN7 expression groups of ccRCC patients. We found that cancer pathways and EMT-related pathways both decreased in ccRCC patients with high CLDN7 expression . Furthermore, RNA-Seq found that different expressed genes influenced by CLDN7 overexpression were enriched in pathways related to cancer and EMT, which was confirmed by qRT-PCR, Western blot and IHC staining in vitro and in vivo. Taken together, it is suggested that CLDN7 may have a fundamental role in tumor progression in ccRCC, by downregulating genes in pathways relating to cancer and EMT at the transcriptional level.
Previous study has shown that TGFB1 exposure decreased expression of CLDN7 and diminished epithelial barrier function, however, CLDN7 overexpression resulted in protection from TGFB1-mediated barrier dysfunction . Many transcription repressors, including TWIST1, ZEB-1, Musashi-2, and Snail, promote EMT and can bind to E-box motifs, thus suppressing CLDN7 expression [25,26,27,28,29,30]. CLDN7, as a member of tight junctions, is often considered to be downstream of EMT. Interestingly, a recent study reported that overexpression of CLDN7 in colon cancer induced epithelial features and suppressed EMT through upregulation of Rab25 then decreased expression of p-Src and mitogen-activated protein kinase–extracellular signal–regulated kinase 1/2 . However, the present study found that in CLDN7 overexpressed ccRCC cells, Rab25 did not increased significantly. Notably, we discovered that CLDN7 downregulated TGF-beta signaling pathway significantly. Additionally, we found that the EMT inducer, TGFB1, was decreased by CLDN7. This result was validated in vitro and in vivo. Therefore, we propose a possible vicious circle between the loss of CLDN7 and upregulation of TGFB1 in ccRCC carcinogenesis and development. Further investigation of this cycle is part of the ongoing studies in our laboratory.
Finally, given the suppressive role of CLDN7 during the process of EMT induced by TGFB1 and Musashi-2 [24, 27], it is a novel strategy to inhibit tumor progression by increasing CLDN7 expression. Studies have demonstrated that CLDN7 overexpression inhibited human colon and lung cancer invasion though EMT and MAPK pathways [19, 31]. In addition, the EMT process itself has been shown to influence cellular resistance to a number of drugs . And CLDN7 was revealed to increase chemosensitivity through the activation of caspase pathway in lung cancer . In this study, we found that CLDN7 downregulation was associated with poorer prognosis and CLDN7 overexpression inhibited EMT-related pathways in ccRCC. Therefore, we propose that CLDN7 may serve as a biomarker or even a therapeutic target for ccRCC.
In conclusion, this is the first report to show that downregulation of CLDN7 is closely associated with metastatic features and poor prognosis in ccRCC patients. We also proposed the mechanism that promoter hypermethylation contributes to downregulation of CLDN7 expression. Moreover, CLDN7 promoter hypermethylation is also related to adverse pathologic results and poor prognosis. CLDN7 overexpression could inhibit cell growth and metastasis of ccRCC cells in vitro and in vivo. Notably, we have demonstrated a previously undescribed and important role of CLDN7 in cancer-related and EMT-related pathways in ccRCC. Further investigations are required to confirm the potential clinical application of CLDN7 and elucidate the detailed molecular mechanisms between CLDN7 and TGFB1 in ccRCC.
Bisulfite genomic sequencing PCR
Clear cell renal cell carcinoma
Chromophobe renal cell carcinoma
- E-cadherin, CDH1:
International Society of Urological Pathology
- N-cadherin, CDH2:
Polymerase Chain Reaction
Papillary renal cell carcinoma
Quantitative real-time PCR
Renal cell carcinoma
High-throughput cDNA sequencing
Reverse transcription PCR
The Cancer Genome Atlas
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2017. CA Cancer J Clin. 2017;67:7–30.
Choueiri TK, Motzer RJ. Systemic therapy for metastatic renal-cell carcinoma. N Engl J Med. 2017;376:354–66.
Srigley JR, Delahunt B, Eble JN, Egevad L, Epstein JI, Grignon D, Hes O, Moch H, Montironi R, Tickoo SK, et al. The International Society of Urological Pathology (ISUP) Vancouver classification of renal neoplasia. Am J Surg Pathol. 2013;37:1469–89.
Nickerson ML, Jaeger E, Shi Y, Durocher JA, Mahurkar S, Zaridze D, Matveev V, Janout V, Kollarova H, Bencko V, et al. Improved identification of von Hippel-Lindau gene alterations in clear cell renal tumors. Clin Cancer Res. 2008;14:4726–34.
Surveillance, Epidemiology, and End Results Program. SEER stat fact sheets: kidney and renal pelvis cancer. Bethesda: National Cancer Institute. http://seer.cancer.gov/statfacts/html/kidrp.html. Accessed 20 Aug 2018.
Klatte T, Pantuck AJ, Kleid MD, Belldegrun AS. Understanding the natural biology of kidney cancer: implications for targeted cancer therapy. Rev Urol. 2007;9:47–56.
Zisman A, Pantuck AJ, Wieder J, Chao DH, Dorey F, Said JW, deKernion JB, Figlin RA, Belldegrun AS. Risk group assessment and clinical outcome algorithm to predict the natural history of patients with surgically resected renal cell carcinoma. J Clin Oncol. 2002;20:4559–66.
Hirohashi S, Kanai Y. Cell adhesion system and human cancer morphogenesis. Cancer Sci. 2003;94:575–81.
Tsukita S, Furuse M, Itoh M. Multifunctional strands in tight junctions. Nat Rev Mol Cell Biol. 2001;2:285–93.
Choi YD, Kim KS, Ryu S, Park Y, Cho NH, Rha SH, Jang JJ, Ro JY, Juhng SW, Choi C. Claudin-7 is highly expressed in chromophobe renal cell carcinoma and renal oncocytoma. J Korean Med Sci. 2007;22:305–10.
Hornsby CD, Cohen C, Amin MB, Picken MM, Lawson D, Yin-Goen Q, Young AN. Claudin-7 immunohistochemistry in renal tumors: a candidate marker for chromophobe renal cell carcinoma identified by gene expression profiling. Arch Pathol Lab Med. 2007;131:1541–6.
Lechpammer M, Resnick MB, Sabo E, Yakirevich E, Greaves WO, Sciandra KT, Tavares R, Noble LC, DeLellis RA, Wang LJ. The diagnostic and prognostic utility of claudin expression in renal cell neoplasms. Mod Pathol. 2008;21:1320–9.
Safran M, Kim WY, O'Connell F, Flippin L, Gunzler V, Horner JW, Depinho RA, Kaelin WG Jr. Mouse model for noninvasive imaging of HIF prolyl hydroxylase activity: assessment of an oral agent that stimulates erythropoietin production. Proc Natl Acad Sci U S A. 2006;103:105–10.
Li Y, Chen D, Li Y, Jin L, Liu J, Su Z, Qi Z, Shi M, Jiang Z, Ni L, et al. Oncogenic cAMP responsive element binding protein 1 is overexpressed upon loss of tumor suppressive miR-10b-5p and miR-363-3p in renal cancer. Oncol Rep. 2016;35:1967–78.
Sharma S, Kelly TK, Jones PA. Epigenetics in cancer. Carcinogenesis. 2010;31:27–36.
Kominsky SL, Argani P, Korz D, Evron E, Raman V, Garrett E, Rein A, Sauter G, Kallioniemi OP, Sukumar S. Loss of the tight junction protein claudin-7 correlates with histological grade in both ductal carcinoma in situ and invasive ductal carcinoma of the breast. Oncogene. 2003;22:2021–33.
Lourenco SV, Coutinho-Camillo CM, Buim ME, de Carvalho AC, Lessa RC, Pereira CM, Vettore AL, Carvalho AL, Fregnani JH, Kowalski LP, Soares FA. Claudin-7 down-regulation is an important feature in oral squamous cell carcinoma. Histopathology. 2010;57:689–98.
Nakayama F, Semba S, Usami Y, Chiba H, Sawada N, Yokozaki H. Hypermethylation-modulated downregulation of claudin-7 expression promotes the progression of colorectal carcinoma. Pathobiology. 2008;75:177–85.
Bhat AA, Pope JL, Smith JJ, Ahmad R, Chen X, Washington MK, Beauchamp RD, Singh AB, Dhawan P. Claudin-7 expression induces mesenchymal to epithelial transformation (MET) to inhibit colon tumorigenesis. Oncogene. 2015;34:4570–80.
Singh AB, Dhawan P. Claudins and cancer: fall of the soldiers entrusted to protect the gate and keep the barrier intact. Semin Cell Dev Biol. 2015;42:58–65.
Lee JW, Lee SJ, Seo J, Song SY, Ahn G, Park CS, Lee JH, Kim BG, Bae DS. Increased expressions of claudin-1 and claudin-7 during the progression of cervical neoplasia. Gynecol Oncol. 2005;97:53–9.
Akimoto T, Takasawa A, Murata M, Kojima Y, Takasawa K, Nojima M, Aoyama T, Hiratsuka Y, Ono Y, Tanaka S, et al. Analysis of the expression and localization of tight junction transmembrane proteins, claudin-1, −4, −7, occludin and JAM-A, in human cervical adenocarcinoma. Histol Histopathol. 2016;31:921–31.
Ding L, Wang L, Sui L, Zhao H, Xu X, Li T, Wang X, Li W, Zhou P, Kong L. Claudin-7 indirectly regulates the integrin/FAK signaling pathway in human colon cancer tissue. J Hum Genet. 2016;61:711–20.
Nguyen N, Fernando SD, Biette KA, Hammer JA, Capocelli KE, Kitzenberg DA, Glover LE, Colgan SP, Furuta GT, Masterson JC. TGF-β1 alters esophageal epithelial barrier function by attenuation of Claudin-7 in eosinophilic esophagitis. Mucosal Immunol. 2018;11:415–26.
Katoh M. Epithelial-mesenchymal transition in gastric cancer (review). Int J Oncol. 2005;27:1677–83.
Larsen JE, Nathan V, Osborne JK, Farrow RK, Deb D, Sullivan JP, Dospoy PD, Augustyn A, Hight SK, Sato M, et al. ZEB1 drives epithelial-to-mesenchymal transition in lung cancer. J Clin Invest. 2016;126:3219–35.
Kudinov AE, Deneka A, Nikonova AS, Beck TN, Ahn YH, Liu X, Martinez CF, Schultz FA, Reynolds S, Yang DH, et al. Musashi-2 (MSI2) supports TGF-beta signaling and inhibits claudins to promote non-small cell lung cancer (NSCLC) metastasis. Proc Natl Acad Sci U S A. 2016;113:6955–60.
Wroblewski LE, Piazuelo MB, Chaturvedi R, Schumacher M, Aihara E, Feng R, Noto JM, Delgado A, Israel DA, Zavros Y, et al. Helicobacter pylori targets cancer-associated apical-junctional constituents in gastroids and gastric epithelial cells. Gut. 2015;64:720–30.
Lu Y, Yu L, Yang M, Jin X, Liu Z, Zhang X, Wang L, Lin D, Liu Y, Wang M, Quan C. The effects of shRNA-mediated gene silencing of transcription factor SNAI1 on the biological phenotypes of breast cancer cell line MCF-7. Mol Cell Biochem. 2014;388:113–21.
Mironchik Y, Winnard PT, Vesuna F, Kato Y, Wildes F, Pathak AP, Kominsky S, Artemov D, Bhujwalla Z, Van Diest P, et al. Twist overexpression induces in vivo angiogenesis and correlates with chromosomal instability in breast Cancer. Cancer Res. 2005;65:10801–9.
Lu Z, Ding L, Hong H, Hoggard J, Lu Q, Chen YH. Claudin-7 inhibits human lung cancer cell migration and invasion through ERK/MAPK signaling pathway. Exp Cell Res. 2011;317:1935–46.
Singh A, Settleman J. EMT, cancer stem cells and drug resistance: an emerging axis of evil in the war on cancer. Oncogene. 2010;29:4741–51.
Hoggard J, Fan J, Lu Z, Lu Q, Sutton L, Chen YH. Claudin-7 increases chemosensitivity to cisplatin through the upregulation of caspase pathway in human NCI-H522 lung cancer cells. Cancer Sci. 2013;104:611–8.
We thank the staff at the Department of Urology, Peking University Hospital, Beijing, China, and the Institute of Urology, Peking University, Beijing 100034, China, for technical support.
This study was funded by the National Key R&D Program of China (2016YFC0902601), National Natural Science Foundation of China (No. 81602253, No. 81672546, No. 81372746 and No. 81772703) and Natural Science Foundation of Beijing (7172219).
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. TCGA Kidney Clear Cell Carcinoma, Papillary Cell Carcinoma and Chromophobe CLDN7 mRNA expression data, methylation beta value and clinical data were downloaded from UCSC Xena (https://xenabrowser.net/heatmap/).
Ethics approval and consent to participate
All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. This study was approved by the Institutional Review Board of the Peking University First Hospital. Animal experiments were conducted in accordance with the NIH Guidelines for the Care and Use of Laboratory Animals with the approval of the Review Board of Peking University First Hospital, Beijing.
Consent for publication
Written informed consent for publication of their clinical details and/or clinical images was obtained from the patient/parent/guardian/ relative of the patient. A copy of the consent form is available for review by the Editor of this journal.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. The clinicopathological features of 120 ccRCC patients from Peking University First Hospital. (DOCX 14 kb)
Figure S1. Transfection validation. (JPG 2066 kb)
Table S2. The primers designed for this study. (DOCX 14 kb)
Table S3. Correlation between CLDN7 expression and 653 clinicopathological features in 534 ccRCC patients from TCGA. (DOCX 15 kb)
Figure S2. Six DNA methylation sites of CLDN7 promoter that not significantly hypermethylated in TCGA ccRCC dataset. (JPG 1171 kb)
Figure S3. The six DNA methylation regions of CLDN7 promoter that are negatively correlated with CLDN7 mRNA expression in TCGA ccRCC dataset. (JPG 1164 kb)
Table S4. Correlation between CLDN7 promoter DNA methylation site (cg00072720) and clinicopathological features in 319 ccRCC patients from TCGA. (DOCX 15 kb)
Figure S4. The CLDN7 promoter DNA methylation site, cg00072720, was associated with poor overall survival time while in hypermethylated status. (JPG 470 kb)
Figure S5. Gene-set enrichment analysis is used to identify the pathways in two different CLDN7 mRNA level groups. (JPG 2095 kb)
Table S5. Gene-set enrichment analysis between high- and low- CLDN7 group in Kidney clear cell carcinoma (KIRC) cohort from TCGA (532 cases). (DOCX 17 kb)
About this article
Cite this article
Li, Y., Gong, Y., Ning, X. et al. Downregulation of CLDN7 due to promoter hypermethylation is associated with human clear cell renal cell carcinoma progression and poor prognosis. J Exp Clin Cancer Res 37, 276 (2018). https://doi.org/10.1186/s13046-018-0924-y