Abstract
The Cancer Genome Atlas (TCGA) and other large-scale genomic data pipelines have been integral to the current understanding of the molecular events underlying renal cell carcinoma (RCC). These data networks have focused mostly on primary RCC, which often demonstrates indolent behavior. However, metastatic disease is the major cause of mortality associated with RCC and data sets examining metastatic tumors are sparse. Therefore, a more comprehensive analysis of gene expression and DNA methylome profiling of metastatic RCC in addition to primary RCC and normal kidney was performed. Integrative analysis of the methylome and transcriptome identified over 30 RCC-specific genes whose mRNA expression inversely correlated with promoter methylation, including several known targets of hypoxia inducible factors. Notably, genes encoding several metabolism-related proteins were identified as differentially regulated via methylation including hexokinase 2, aldolase C, stearoyl-CoA desaturase, and estrogen-related receptor-γ (ESRRG), which has a known role in the regulation of nuclear-encoded mitochondrial metabolism genes. Several gene expression changes could portend prognosis in the TCGA cohort. Mechanistically, ESRRG loss occurs via DNA methylation and histone repressive silencing mediated by the polycomb repressor complex 2. Restoration of ESRRG in RCC lines suppresses migratory and invasive phenotypes independently of its canonical role in mitochondrial metabolism.
Implications: Collectively, these data provide significant insight into the biology of aggressive RCC and demonstrate a novel role for DNA methylation in the promotion of HIF signaling and invasive phenotypes in renal cancer.
Introduction
Insight into the biologic basis of renal cancer has greatly benefitted from large-scale data sets such as The Cancer Genome Atlas (TCGA). These data have confirmed the high prevalence of tumor-initiating events such as alterations in the VHL tumor suppressor gene, which lead to stabilization of hypoxia inducible factors (HIF1α and HIF2α; refs. 1–3). Analyses of these datasets have also indicated the importance of epigenetics to renal cell carcinoma (RCC) malignancy, including the alteration of PBRM1, which encodes a component of the SWI/SNF chromatin remodeling complex and BAP1, a gene encoding a known histone deubiquitinase (1, 4). Although TCGA analysis has been invaluable to our understanding of primary ccRCC, far less is known about metastatic ccRCC.
Development of metastasis is the major cause of mortality associated with ccRCC. Patients with metastatic ccRCC have a poor prognosis, and remain largely incurable despite recent advances with combined T-cell checkpoint blockade (5–8). To date, studies of metastatic ccRCC using patient-derived samples has been limited. One recent pan-cancer analysis of metastasis included RCC; however, only eight of the 500 tumor samples were RCC (9). Two recent studies focused on pharmacogenes and kinases with altered expression in metastatic RCC (10, 11). However, integrative studies performed in an unbiased manner are lacking.
DNA methylation plays a significant role in the regulation of gene expression. Probably the most well-characterized alteration is gene promoter methylation that silences the corresponding gene expression. In RCC, VHL inactivation through promoter methylation has been identified (12). Although other genes have been described as methylated in RCC through cell line based studies, a recent analysis indicated that the methylomes of cultured RCC cells vastly differs from that of patient samples (11). The TCGA analysis of ccRCC (referred to as KIRC) demonstrated that increased promoter methylation of gene correlated with higher stage and higher grade tumors (1). Recent studies by Sato and colleagues performed an analysis of over 100 primary ccRCC and identified that a DNA hypermethylator phenotype defined a subgroup of patients with worsened overall survival (13). Although the TCGA analysis provides large-scale data sets, data do not include metastatic tumor tissues. Because gene signatures related to DNA methylation might differ between primary and metastatic RCC, we undertook a combined analysis of the transcriptome and methylome of metastatic as well as primary ccRCC and normal kidney to gain insight into the biological basis of aggressive ccRCC.
We performed concurrent transcriptomic and methylome analyses via array-based methodologies followed by an integrative analysis of these two data sets. These data represent the first such combined analysis for ccRCC that includes normal kidney, primary RCC and metastatic RCC tissues. Our analysis indicates a role for promoter methylation, both hyper- and hypomethylation, in the regulation of genes involved in cellular metabolism. We identified several genes silenced via promoter methylation and whose loss was correlated with worsened prognosis through subsequent analysis of the TCGA cohort. Notably, our data demonstrate a role for promoter hypomethylation in the amplification of HIF signaling in aggressive RCC. Our data also implicate novel genes whose deregulation is associated with tumor progression and invasive phenotypes.
Materials and Methods
Patient samples
Samples for array (normal = 9, primary = 9, metastasis n = 26) based studies were acquired from the NCI Cooperative Human Tissue Network (CHTN) and were not patient-matched. Patient-matched samples (n = 11/group) used for validation RT-qPCR studies were obtained from patients undergoing clinically indicated procedures at UAB Hospital and have been previously described (14). Further detailed information of these patient-matched samples is provided in Supplementary Table S4. All studies were performed in accord with the institutional IRB.
Cell culture
RCC cell lines were purchased from ATCC except RXF-393 (NCI) and RCC4 (P. Ratcliffe, Oxford). RCC cells RCC4, RCC10, A498, 786-0, and HK2 proximal tubular epithelial cells were grown in DMEM with 10% FBS with penicillin–streptomycin (100 U/mL) in 5% CO2 at 37°C. CAKI-1 cells were cultured in MEM medium containing 10% FBS with antibiotics. Cells were used within 10 passages of the initial stock with periodic testing for Mycoplasma infection. No further characterization was performed.
GSK126 treatment
CAKI-1 cells were seeded on six-well plates and treated with GSK126 (Cayman Chemical) at the indicated dose for 72 hours.
Plasmids and virus infections
Sequence-verified shRNA lentiviral plasmid vectors for the human EZH2 gene and a nonspecific shRNA control (shCTRL) constructs have been previously described (15). Lentiviral transduction was performed as previously described (16). PPARGC1A and ESRRG cDNAs were transduced in RCC cells using adenovirus-mediated gene transfer (Vector Biolabs) for 48 hours.
Gene expression profiling
Gene expression analysis was performed using the human HT-12 BeadChip and iScan system from Illumina, Inc. Total RNA was converted to cDNA by reverse transcription, followed by second-strand synthesis to generate double-stranded cDNA. After purification, the cDNA was converted to biotin-labeled cRNA, hybridized to a human HT-12 BeadChip, and stained with strepavidin-Cy3 for visualization. The human HT-12 BeadChips contain sequences representing approximately 46,000 curated and putative genes and ESTs. Quality standards for hybridization, labeling, staining, background signal, and basal level of housekeeping gene expression for each chip were verified. After scanning the probe array, the resulting image was analyzed using the GenomeStudio software (Illumina, Inc.). Raw data obtained from Illumina's beadstudio were preprocessed using a bioconductor package “beadarray” (17) in the open source software R v3.2.2 [https://cran.r-project.org/]. Raw signal values were quantile normalized and log2 transformed. Probes with no-match or of poor quality were then filtered out, leaving 34,476 probes for further analysis. Differential expression analysis was performed using the limma package (18). Only probes with an absolute fold change of two or more and a P-value of 0.05 or less were considered differentially expressed. Probes were annotated to their respective genes using the Bioconductor annotation package “illuminaHumanv4.db” (19). The gene expression profiling data have been deposited at Gene Expression Omnibus (GEO; #GSE105288).
Whole genome methylation analysis
Genomic DNA was isolated from 25 mg of fresh frozen tissues using the Invitrogen PureLink Genomic DNA Kit (Thermo Fisher Scientific). DNA samples were bisulfite converted with the Zymo EZ DNA Methylation Kit (Zymo Research). Genome-wide methylation profiling was accomplished with the Illumina Infinium HumanMethylation450 BeadChip Kit, which interrogates DNA methylation at 485,755 CpG sites. The BeadChip arrays were processed according to the manufacturer's instructions and scanned using an Illumina iScan system (Illumina). The Illumina GenomeStudio Methylation Module software was used to analyze the BeadChip image data. The raw intensity data were loaded to a biocondutor package “minfi” (20). The raw data were normalized using the subset-quantile within array normalization (SWAN) method and probes with a detection P value ≥0.01 in at least one sample were excluded from further analysis. Using the getBeta function of minfi, β values were calculated for each probe and M values were calculated using the formula: log2 ([methylated_probe+100]/[unmethylated_probe+100]). Differentially methylated probes were identified using M values determined via the “limma” package (18). Probes with an absolute fold change of 2 or more and P value of <0.05 were chosen as differentially methylated.
In silico validation
Using the Oncomine database (21), the expression level of 35 RCC associated genes was checked in public RCC microarray data sets (22–24).
In addition, UALCAN (25) was used to study the expression level of these genes in 72 normal kidney and 533 primary RCC samples and their expression effect on patient survival. UALCAN uses “level 3 RNA-seq data” from TCGA [https://tcga-data.nci.nih.gov/].
The promoter CpG methylation status of 35 ccRCC-associated genes was validated by KIRC TCGA 450 K methylation data. The raw data for 160 kidney normal and 325 KIRC primary tumors were downloaded from the TCGA portal and processed using the same method as above.
qRT-PCR
Total RNA from tissues was isolated using the RNeasy Mini Kit (Qiagen). Total RNA was extracted from cultured cells using Trizol reagent (Ambion). cDNA was generated using the high-capacity cDNA Reverse Transcription Kit (Applied Biosystems). Real-time PCR (RT-PCR) was performed using TaqMan gene expression reagents mixed with TaqMan primers and indicated cDNAs in QuantStudio 6K Flex Real-Time PCR System (Applied Biosystem). Target gene mRNA expression was normalized to either ribosomal protein (RPLPO) or TATA binding protein (TBP). Normalized CT values were calculated via the ΔΔCt analysis. The following probes were used: ALDOA (Hs00605108_g1), ALDOB (Hs01551887_m1), ALDOC(Hs00902799_g1), ESRRA (Hs01067166_g1), ESRRB (Hs01584024_m1),ESRRG (Hs00976243_m1), ACO2 (Hs00426616_g1), OGDH (Hs01081865_m1), SUCLG1 (Hs00388749_m1), MDH2 (Hs00938918_m1), TBP (Hs00427620_m1), and RPLP0 (Hs99999902_m1).
Immunoblotting analysis
Immunoblotting was performed as described previously (16). Antibodies used included EZH2 (Cell Signaling Technology; #5246), ESRRG (Proteintech; #14017-1-AP), α-tubulin (Sigma-Aldrich; # T5168), β-actin (Abcam; #Ab8226).
Proteomics analysis
RCC cells and human kidney tissues (n = 5/ group) were lysed in T-per lysis buffer with protease inhibitor. Following the BCA protein assay, 40 μg of each lysate was separated onto a 10% SDS-PAGE gels and the band around the 40 kDa was excised and digested with trypsin overnight prior to LC/MS analysis.
Wound healing and migration assay
RCC cells were seeded in six-well plates in growth medium and transfected with the indicated adenovirus for 48 hours. The cells were mechanically disrupted with a pipette tip, and cell migration was observed with an EVOS FL imaging system (Invitrogen). For the transwell migration assay, RCC cells were seeded in serum-free medium onto the upper chamber of the transwell (Corning) for 24 hours. The migrated cells were fixed and quantified by crystal violet staining.
Matrigel-based invasion assay
Cells were infected with the indicated adenovirus before seeding in the upper chambers of BD BioCoat Matrigel Matrix (Corning). After 16 hours, cells on the upper chamber (non-invading cells) were removed with a cotton swab, whereas cells on the lower side of the chamber (invading cells) were stained with Diff-Quik Stain Kit (Siemens). Invading cells were observed and counted under a light microscope.
Methylated DNA immunoprecipitation
Genomic DNA was extracted from RCC cells and human kidney tissues using the Wizard genomic DNA Isolation Kit (Promega). DNA immunoprecipitations were performed with 2 μg of genomic DNA and either 5-methylcytosine (5-mC) rabbit mAb or rabbit mAb IgG using SimpleDIP methylated DNA IP (MeDIP) Kit (Cell Signaling Technology). The enriched DNA was quantified by real-time PCR (Applied Biosystem). The human ESRRG primers were specifically designed for heavily methylated regions in the ESRRG promoter (5′-GGAGTTGCTGCAGTTTCAAAG-3′; 5′-CGTTCCGAGTCAGCTCCTA-3′). The enrichment of immunoprecipitated DNA was represented as signal relative to the total 10% of input DNA.
Results
Methylome/ transcriptome analysis of normal kidney, primary tumor, and metastasis
To obtain novel insight into the biology of aggressive kidney cancer, we obtained normal kidney, primary ccRCC, and metastatic ccRCC tumor deposits from the Cooperative Human Tissue Network (CHTN). DNA and RNA were extracted from tissue and analyzed via methylome and transcriptomic array analysis, respectively. Sites of metastasis for DNA/RNA array analysis are from sites including adrenal, pancreas, lung, and liver (Supplementary Table S1). Both primary and metastatic tumor were included to assess if alterations in the transcriptome/methylome observed were part of a continuum that paralleled tumor progression. Principal component analysis (PCA) of transcriptomic data demonstrated that metastatic tumor samples clustered together (Supplementary Fig. S1). The pipeline for our integrative analysis is demonstrated in Figure 1. In summary, we first identified genes demonstrating RCC-specific expression alterations. In tandem, we identified CpG sites (i.e., probes) with differential methylation levels in RCC (either primary and/or metastatic tumor). We focused specifically on CpG sites in the expected region of the gene promoter which included probes assigned to the 200 base pairs of the transcription start site (TSS200), within 1,500 base of the TSS (TSS1500), 5′ untranslated regions (5′UTR), and the first exon. These analyses identified 35 RCC specific genes changes demonstrating a statistically significant inverse correlation between mRNA expression status and DNA promoter methylation status (Fig. 2A). In general, both primary and metastatic tumor demonstrated similar directionality with respect to gene expression changes and methylation. However, metastatic tumor often demonstrated more prominent changes compared with primary tumor (Fig. 2A; Supplementary Table S2). For example, HK2 (hexokinase 2) and SLC16A3, which encodes for monocarboxylate transporter 4 (MCT-4), had increased expression in metastasis and a corresponding decrease in promoter methylation.
Analysis pipeline. Schematic representation of analysis flow to identify clear cell RCC associated genes based on gene expression and DNA methylation profile.
Transcriptome and methylome alterations in biologically aggressive kidney cancer. A, Heatmaps showing gene expression and DNA methylation profile of 35 RCC associated genes. Fold change value for each gene is represented by the probe with highest value. Fold change cut-off of 2 and P-value limit of 0.05 were considered to obtain differentially expressed or differentially methylated probes. B, Survival analysis of down-regulated genes SLC25A25, ESRRG, TEMEM213, MYH10, CLCNKB, and ATP6V0A4 using TCGA data. The KM plots were obtained from UALCAN (http://ualcan.path.uab.edu), which uses level 3 RNA-seq data and survival information from TCGA kidney renal clear cell carcinoma (KIRC) patients (n = 533). The significance of difference between survival curves was obtained by log-rank test.
External validation of transcriptome/methylome analyses
We next analyzed gene expression changes in publicly available data sets that mainly analyzed primary ccRCC. Similar gene expression changes were found in the TCGA data set, which performed transcriptomic analysis via RNA sequencing (Supplementary Figs. S2 and S3). The promoter methylation status of the 35 RCC specific genes was also consistent with KIRC TCGA 450 K methylation data (Supplementary Table S3). Further validation was obtained via analysis of publicly available microarray data sets utilizing the Oncomine analysis platform (Supplementary Figs. S4, S5, and S6). Given the more prominent changes in gene expression in metastatic tumors for several genes identified in this analysis, we analyzed TCGA RNAseq data to assess whether gene expression levels in primary tumor had prognostic significance. For several genes downregulated in metastatic RCC, lower expression levels were associated with worsened survival in TCGA analysis (Fig. 2B). In contrast, only one of the upregulated genes in metastatic RCC (LGALS1) demonstrated worsened prognosis with higher expression in the TCGA RNAseq data set (Supplementary Fig. S7).
HIF target genes are promoter hypomethylated with increased expression
Several HIF target genes were hypomethylated with increased expression in metastatic RCC, including ADM (adrenomedullin), TNFAIP6 (TNF alpha induced protein 6), SLC16A3 (monocarboxylate transporter 4), CAV1 (caveolin-1), and LGALS1 (galectin-1; refs. 26–30). A similar pattern was found for the HIF-responsive glycolysis genes HK2 (hexokinase 2) and ALDOC (aldolase C). Our data demonstrate that promoter methylation is an additional layer of regulation in the expression of HIF-regulated genes in RCC. Adolases catalyze the interconversion between fructose 1,6-bisphospate and the intermediates dihydroxyacetone phosphate and glyceraldehyde 3-phosphate. Prior studies have demonstrated alteration in the protein expression of ALDO isoforms in primary RCC (31, 32). In our study, mRNA expression of ALDOA and ALDOC was higher in a panel of RCC lines relative to untransformed HK2 renal epithelial cells (Figs. 3A and B). ALDOB mRNA levels were extremely low (data not shown). We also examined aldolase isoform expression using a proteomics-based approach as commercial antibodies demonstrated significant cross-reactivity between ALDO isoforms (data not shown). RCC lines demonstrated increased ALDOA and ALDOC protein levels relative to untransformed renal epithelial RPTEC lines (Fig. 3C). All cultured lines had very low to undetectable levels of ALDOB relative ALDOA and ALDOC (data not shown). We next utilized proteomics to assess the relative abundance of each ALDO isoform in RCC tissues relative to normal kidney. We identified increased ALDOA and ALDOC protein levels in tumor relative to normal kidney (Fig. 3D). ALDOB protein levels were much lower in tumor than normal kidney. We next determined ALDO isoform mRNA expression as a function of tumor progression by assessing transcript levels in normal kidney, primary tumor, and metastatic tumor deposits in patient-matched samples (Fig. 3E), identifying a trend toward increased transcript levels of both ALDOA and ALDOC with tumor progression.
The differential distribution of aldolases in RCC cells. A and B, Relative mRNA expressions of ALDOA and ALDOC in RCC cells. Total RNA was isolated from RCC cells and renal proximal tubular epithelial cells (HK2; n = 4/group). Transcript levels were normalized to those of RPLPO. Values are means ± SEM and asterisks indicate the differences relative to HK2 cells (*, P < 0.05; **, P < 0.01). C, Relative protein levels of aldolases in RCC cells compared with primary renal proximal tubule epithelial cells (RPTEC) by LC/MS proteomic analysis. D, Relative protein levels of aldolases in normal and tumor tissues by LC-MS proteomic analysis (n = 5/group). Asterisks indicate difference relative to normal kidney tissue (*, P < 0.05; **, P < 0.01). E, Relative mRNA expressions of aldolase genes in patient-matched normal kidney, primary and metastasis were measured by RT-qPCR (n = 11/group). Values are means ± SEM and asterisks indicate the differences relative to normal kidney tissue (*, P < 0.05; **, P < 0.01).
Characterization of genes silenced in RCC via promoter methylation
We identified promoter hypermethylation and silencing of several genes previously reported to be reduced in RCC including ATP1A1 (ATPase Na+/K+ transporting subunit alpha 1), ATP6V0A4 (ATPase H+ transporting V0 subunit A4), KCNJ1 (potassium voltage-gated channel subfamily J member 1), CLCNKB (chloride voltage-gated channel Kb), PLXNB1 (plexin B1), and TMEM213 (transmembrane protein 213; refs. 33–37). A novel finding was reduced expression of the gene encoding the estrogen receptor-related (ERR)-γ (ESRRG). The ERRs consist of three members (α, β, and γ) encoded by ESRRA, ESRRB, and ESRRG, respectively. The ERRs belong to the superfamily of nuclear hormone receptors. The ERRs most prominent role is the transcriptional regulation of mitochondrial metabolism. We validated the loss of ESRRG expression in patient-matched primary and metastatic tumor relative to normal kidney (Fig. 4A). Data on these samples including sites of metastasis and prior systemic therapy is provided in Supplementary Table S4. Although the greatest fold reduction in ESRRG expression was found in primary tumor relative to normal kidney, a further decline in expression was noted in metastatic tumor deposits. Consistent with these data, we demonstrated significantly reduced ESRRG expression in RCC lines relative to normal kidney (Fig. 4B). We next performed immunoblotting in another sample cohort and identified reduced ERR-γ protein in primary tumor relative to normal kidney (Fig. 4C). We also assessed the mRNA expression of the other ERRs in normal and tumor and identified modest reductions in ESRRA in tumor (Fig. 4D). Similar to ESRRG, we identified marked reductions in the expression of ESRRB (Fig. 4D). Similar changes in mRNA expression of the ERRs were identified in analyses of TCGA data (data not shown). We next examined the prognostic significance of the various ERR isoforms. As previously noted, low ESRRG expression was associated with worsened overall survival in the TCGA data. In contrast, no significant associations with survival were observed based on expression levels for either ESRRA or ESRRB (Fig. 4E).
Suppression of ESRRG in RCC. A, Relative mRNA expression of ESRRG in patient-matched normal kidney, primary tumor, and metastatic tumor (n = 11/group). Asterisks indicate differences relative to normal kidney tissue (**, P < 0.01). B, ESRRG mRNA expression in a panel of RCC cell lines relative to normal kidney. Transcript levels were normalized to those of TATA-binding protein (TBP). C, Western blot analysis of ESRRG in normal (N) and tumor (T). D, Relative mRNA expression of ESRRA and ESRRB in patient-matched normal kidney, primary tumor, and metastatic tumor (n = 11/group). Asterisks indicate differences relative to normal kidney tissue (**, P < 0.01). E, Kaplan–Meier survival curve analysis of patients from the TCGA data set on ccRCC based on expression of ESRRA and ESRRB.
Mechanism of ESRRG silencing in RCC
We sought further insight into the silencing of ESRRG in renal cancer given that its expression had prognostic significance. Analysis using the UCSC genome browser identified to two CpG islands near the 5′ region of the ESRRG gene (Fig. 5A). We identified multiple CpG sites within a CpG island situated right at the transcription start site whose methylation levels inversely correlated with ESRRG mRNA expression (Fig. 5A, CpG island #1). We therefore performed ChIP studies using an antibody-specific to 5-mC. ChIP qPCR utilizing primers to this CpG island demonstrated strong enrichment of 5-mC in this region in three RCC cell lines and tumors from patients (Fig. 5B–D). Recent studies of the methylation landscape of ccRCC indicate that many DNA hypermethylated sites are enriched for targets of the polycomb repressor complex 2 (PRC2; ref. 13). PRC2 lays down the repressive histone mark H3K27me3. Prior studies have demonstrated coordination between DNA methylation and histone repressive marks in the silencing of gene expression. The catalytic subunit, enhancer of zeste homolog 2 (EZH2) is required for PRC2 activity. We therefore assessed the role of PRC2 in ESRRG silencing in RCC cells via a knockdown approach. We confirmed EZH2 knockdown via stable transduction with shRNA (Fig. 5E). Knockdown of EZH2 with two nonoverlapping shRNAs led to increased ESRRG expression in CAKI-1 cells (Fig. 5F). We further evaluated the effect of EZH2 inhibition on ESRRG expression using a pharmacologic inhibitor (GSK126). GSK126 reduced H3K27me3 levels consistent with EZH2/PRC2 inhibition (Fig. 5G). Similar to EZH2 knockdown, GSK126 treatment resulted in the increased expression of ESRRG in CAKI-1 cells (Fig. 5H). Collectively, these data demonstrate coordination between DNA and histone methylation in the silencing of ESRRG in RCC.
Silencing of ESRRG gene by CpG methylation in ccRCC. A, A cluster of CpG loci in the 5′ region of the ESRRG gene. The marked CpG loci are heavily methylated in kidney tumors. B and C, Analysis of methylated DNA enrichment by ChIP-qPCR. Genomic DNA was extracted with either (B) CAKI-1 cells, or (C) RCC4 and RXF-393 cells (n = 3/group). D, ChIP was performed on kidney tumor samples with either rabbit IgG or 5-mC antibodies (n = 3/group). The enriched DNA was quantified by real-time PCR with primer sets targeting highly methylated regions of the ESRRG promoter. E, Immunoblot analysis of EZH2 knockdown in CAKI-1 cells transduced with shCTRL or two nonoverlapping EZH2-targeting shRNAs (shRNA EZH2-1 and shRNA2-2). F, CAKI-1 cells were stably transduced with shEZH2-1, shEZH2-2, or control shRNA (ShCTRL). The mRNA levels of ESRRG were analyzed by qRT-PCR (n = 3 per group). G, Immunoblot analysis of H3K27me3 in CAKI-1 cells treated with GSK126 at indicated doses for 72 hours. H, Relative mRNA expression of ESRRG in CAKI-1 cells treated with GSK126 at indicated doses for 72 hours (n = 3/group, three independent experiments).
ESRRG suppresses migratory and invasive phenotypes in RCC
Based on these data, we next assessed the biologic significance of ESRRG silencing in RCC. We transiently expressed ERR-γ via adenovirus. Adenovirus with a GFP construct was used as a control. We first assessed ESRRG mRNA levels in transduced cells. RT-qPCR analysis demonstrated that ESRRG expression was comparable to ESRRG transcript levels in normal kidney (Fig. 6A). We confirmed ERR-γ expression at the protein level via immunoblotting (Fig. 6B). We next assessed for effects on in vitro tumor phenotypes. In both CAKI-1 and RCC4 cells, ERR-γ suppressed wound healing as determined by scratch assay (Fig. 6C and D). RCC cell migration as determined by a transwell assay was significantly reduced by the expression of ERR-γ (Fig. 6E and F). We next assessed cell invasiveness using a Matrigel Boyden chamber assay. In both RCC lines tested, ERR-γ expression led to a significant reduction in cell invasiveness (Fig. 6G and H). As noted previously, ERRs including ERR-γ have well-established roles in mitochondrial metabolism. ERR-γ is known to promote the transcription of nuclear encoded genes related to mitochondrial metabolism including genes encoded tricarboxylic acid (TCA) cycle enzymes. ERR-γ has also been demonstrated to have a role in mesenchymal-to epithelial transition (MET) in line with our invasion assay results (38). We next wanted to assess if ERR-γ's effects on cell invasiveness were related to its role in mitochondrial metabolism. Notably, we did not observe significant effects on the expression of TCA cycle enzymes, known targets of ERRs (Fig. 6I). In contrast, transduction with peroxisome proliferator-activated receptor gamma coactivator 1-α (PGC-1α), a known transcriptional activator of nuclear encoded mitochondrial metabolism genes, induced mRNA expression of several TCA cycle enzymes (Fig. 6I). Collectively, these data demonstrate that ERR-γ suppresses invasive phenotypes independent of its canonical role in mitochondrial metabolism.
ESRRG reverts migratory and invasive behaviors in RCC cells. A, RCC cells were transduced with either Ad-GFP or Ad-ESRRG for 48 hours. Relative mRNA quantification of ESRRG was analyzed by qRT-PCR (n = 3, two independent experiments). B, Western blot analysis of ESRRG in nuclear lysates from RCC cells transduced with either Ad-GFP or Ad-ESRRG. C and D, Representative images of wound healing over time and the measurement of gap distance in RCC cells following ESRRG expression (n = 3/group, three independent experiments). E and F, Representative images and quantification of Transwell migratory assay in RCC cells (n = 3/group, three independent experiments). G and H, Boyden invasion chamber with Matrigel insert was used to assess invasiveness of RCC cells following ESRRG expression (n = 3/group, three independent experiments). All data are presented as mean ± SEM. I, Relative mRNA expression of TCA cycle enzymes in RCC cells transduced with the indicated adenovirus. After 48 hours, mRNA was extracted from RCC cells and analyzed by RT-qPCR for expression of TCA cycle enzymes (n = 3/group).
Discussion
Similar to many other malignancies, RCC has excellent outcomes when organ confined. In contrast, patients with metastatic disease face a poor prognosis. Despite the difference in prognosis between localized and metastatic disease, datasets examining metastatic tumor deposits in the context of RCC are generally lacking. Although prior studies from our group and others have performed transcriptomic analyses of primary and metastatic sites, molecular analyses using epigenetic profiling platforms were not performed (14, 39, 40). A recent study analyzed the transcriptomes and methylomes of both primary and metastatic RCC samples with a focus on pharmacogenes; however, normal samples were not included in the analysis (11). Thus, this study presents a unique data set with insight into the biology of aggressive renal cancer performed in an unbiased manner. The integration with methyl array data provides additional insight into the contribution of epigenetics to observed changes in gene expression.
A notable finding from this study is the presence of multiple alterations in metabolism related genes in metastatic RCC. These include the glycolysis genes HK2 and ALDOC. However, several other HIF responsive genes were identified as hypomethylated in RCC, including ADM, TNFAIP6, and SLC16A3. HIF signaling is activated in the setting of VHL alterations, the most common tumor-initiating event in ccRCC. Hence, the increased expression in primary RCC is consistent with prior studies. However, further increases within metastatic tumor deposits with corresponding promoter hypomethylation indicate an epigenetic basis for the amplification of HIF signaling with tumor progression. To the best of our knowledge, these are the first data from patient-derived samples to implicate a role for methylation in the regulation of HIF target genes in the context of ccRCC. These data are consistent with recent experimental studies demonstrating further increases in HIF signaling beyond the tumor initiating event of VHL loss. Studies by Gao and colleagues demonstrate that inactivation of PBRM1, a gene commonly mutated in ccRCC, promotes HIF signaling in VHL null ccRCC (41). Experimental studies of RCC by Vanharanta and colleagues demonstrate that DNA demethylation promotes the expression of the HIF target gene CYTIP in an in vivo model of ccRCC metastasis (42). Given the large number of HIF target genes, a major challenge in the field has been to identify which of these genes are contributory to the malignant phenotype. Adding another level of complexity in renal cancer is the findings that HIF-1α is tumor suppressive whereas HIF-2α is tumor promoting (43, 44). Hence, the identification of HIF targets genes that are epigenetically amplified through mechanisms such as DNA promoter hypomethylation or chromatin remodeling (e.g., PBRM1 mutation) may reveal those HIF target genes critical to tumor development and/or progression.
We also demonstrate clear evidence for cooperative effects between both DNA and histone repressive methylation in the silencing of ESRRG, which encodes for ERR-γ. The ERRs are known to physically interact with PGC-1α in the promotion on mitochondrial metabolism related genes (45, 46). Our RT-qPCR studies identified markedly reduced expression of ESRRB and ESRRG in metastatic RCC with modest reductions in ESRRA. However, these are the first data to identify the role for histone and DNA repressive methylation of ESRRG in cancer. Our data are in line with prior studies in ccRCC demonstrating that DNA hypermethylated regions are enriched for targets of the PRC2 complex which lays down the repressive H3K27me3 mark (13, 47). Treatment of RCC cells with azacytidine did not result in increased ESRRG expression (data not shown). However, our studies with EZH2 knockdown and EZH2 inhibition suggest that hypermethylated/silenced genes may be reactivated via PRC2 inhibition. EZH2 has been shown to be overexpressed in RCC and other cancer such as prostate, bladder, and breast cancer among others (48–51). Our data are relevant as EZH2 inhibitors are currently in clinical trials.
Loss of ESRRG has previously been implicated in other malignancies including prostate cancer where it impacts proliferation (52, 53). ESRRG has also been implicated in breast cancer. For example, high ESRRG expression is associated with a better prognosis in patients with breast cancer (54). However, data on its role on tumor biology are conflicting. Expression of ERR- γ in MCF-7 breast cancer cells has been shown to promote proliferation (55). In contrast, ERR-γ expression in MDA-MB-231 breast carcinoma cells suppresses invasion (38). We would like to point out that ERR-γ is highly expressed in the kidney in contrast to prostate and breast tissues. Moreover, we expressed ERR-γ in RCC cells at levels comparable to normal kidney to maintain physiologic relevance. Our data demonstrating that ERR-γ re-expression in RCC cells suppresses invasive phenotypes is consistent with a role for this transcription factor in mesenchymal to epithelial transition as previously reported. Further supportive evidence is provided by the fact that lower ESRRG expression correlates with worsened survival. No correlation with survival was found with the genes encoding the other ERR isoforms, ESRRA or ESRRB. Furthermore, despite the strong effects on invasive phenotypes in our studies, we did not see effects on the expression on mitochondrial metabolism genes under our experimental conditions. Hence, our data would support a role for ERR- γ independent of activities known to be redundant among the ERR isoforms (e.g., transcriptional regulation of nuclear encoded mitochondrial metabolism genes).
In addition, we identified several other genes that were hypermethylated and repressed that have been shown to suppress in vitro tumor phenotypes in RCC cells including PLXNB1, ATP1A, and KCNJ1 (33, 35, 37). Our study provides mechanistic insight into the silencing of these genes in kidney cancer. Moreover, they provide further rationale to assess the role of these genes in renal carcinogenesis.
We would like to acknowledge the limitations of this study. Our initial transcriptomic analysis did not utilize patient-matched samples. Moreover, our initial analysis used relatively few normal tissues and primary RCCs. Nevertheless, we were able to validated expression changes of ESRRG in a patient-matched set of tissues. Moreover, many of the changes observed in both our transcriptomic and methylome analyses were consistent with publically available data sets including TCGA data. In addition, there were some discordancy between expression changes and correlation with survival. TMEM178, PLXNB1, TSPAN8 demonstrated reduced expression in RCC but lower expression correlated with improved survival in the TCGA cohort. In contrast, DYSF expression was increased in tumor but higher expression correlated with improved outcomes in the TCGA cohort. One potential explanation is that these factors may have different roles in primary RCC as compared with metastasis.
We also note that the methyl array analysis was performed with standard bisulfite sequencing, which cannot resolve between 5-mC and 5-hydroxymethylcytosine (5-hmC). The conversion of 5-mC to 5-hmC is catalyzed by the ten-eleven translocation enzymes (TETs 1-3). This has traditionally been thought of as the initial step of converting of 5-mC to unmethylated cytosine. However, recent studies indicate that 5-hmC may actually be a stable mark that can accumulate in various genomic regions including enhancer elements and promote gene transcription (56, 57). In contrast to 5-mC, 5-hmC levels tend to have lower levels in promoter regions (58). Based on these data, we therefore focused on methylation in gene promoters that anticorrelated with gene expression. Studies from our lab and others have demonstrated loss of 5-hmC in ccRCC (16, 59). Hence, this may represent an alternate mechanism for gene repression. Future studies should consider these alternate methyl markers and their contribution to gene expression.
In summary, our studies provide critical insight into the biology of aggressive kidney cancers and provide valuable global profiling data of metastatic RCC. We demonstrate an epigenetic basis for the propagation of HIF-related pathways, such as glycolysis, beyond tumor initiating events such as VHL loss. Furthermore, our study is the first demonstrate a role for DNA/repressive methylation in the silencing of ESRRG in cancer and that this process promotes invasive phenotypes in RCC cells. Taken together, our study provides compelling evidence that ccRCC is an epigenetically-driven malignancy. Given the emergence of pharmacologic agents that can target the epigenetic modifiers, our data provide rationale for the consideration of epigenetic-based approaches to advanced kidney cancer.
Disclosure of Potential Conflicts of Interest
Guru Sonpavde reports receiving a commercial research grant from Bayer, Amgen-Onyx, Merck, Pfizer, Boehringer-Ingelheim, Pfizer, and Celgene; has received speakers bureau honoraria from OncLive, Physicians' Education Resource, Research to Practice, Clinical Care Options, and UpToDate; and is a member of Bayer, Sanofi, EMD Serono, Astellas Pharma, Pfizer, Novartis, Eisai, Janssen, Amgen, AstraZeneca, Merck, and Genentech. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: H.-Y. Nam, S. Varambally, S. Sudarshan
Development of methodology: H.-Y. Nam, C.B. Livi, S. Sudarshan
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): S. Shelar, E.-Y. Kho, G. Sonpavde, G. Naik, P. Ghatalia
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): H.-Y. Nam, D.S. Chandrashekar, A. Kundu, S. Shelar, E.-Y. Kho, G. Sonpavde, C.B. Livi, S. Varambally, S. Sudarshan
Writing, review, and/or revision of the manuscript: H.-Y. Nam, G. Sonpavde, G. Naik, S. Varambally, S. Sudarshan
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): G. Sonpavde, S. Sudarshan
Study supervision: S. Sudarshan
Acknowledgments
Array data were generated by the UT Health San Antonio Cancer Center Genomics Shared Resource which is supported by NIH-NCI P30 CA054174 and UT Health San Antonio. Proteomics studies were formed at the UAB Cancer Center Mass Spectrometry/Proteomics (MSP) Shared Facility. We thank Dr. Gregory M. Pavela (University of Alabama, Birmingham) for his editing of the manuscript. Support for this study was provided by Department of Veterans Affairs Merit Award BX002930 (to S. Sudarshan), NIH R01CA200653 (to S. Sudarshan), and NIH R01CA157845 (to S. Varambally). The research described was supported in part by the UAB Comprehensive Cancer Center (P30CA013148).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Footnotes
Note: Supplementary data for this article are available at Molecular Cancer Research Online (http://mcr.aacrjournals.org/).
Current address for C.B. Livi: Agilent Technologies, Santa Clara, CA.
- Received November 1, 2017.
- Revision received April 20, 2018.
- Accepted August 3, 2018.
- Published first August 21, 2018.
- ©2018 American Association for Cancer Research.