8q24 clear cell renal cell carcinoma germline variant is associated with VHL mutation status and clinical aggressiveness
BMC Urology volume 20, Article number: 173 (2020)
The four most commonly-mutated genes in clear cell renal cell carcinoma (ccRCC) tumors are BAP1, PBRM1, SETD2 and VHL. And, there are currently 14 known RCC germline variants that have been reproducibly shown to be associated with RCC risk. However, the association of germline genetics with tumor genetics and clinical aggressiveness are unknown.
We analyzed 420 ccRCC patients from The Cancer Genome Atlas. Molecular subtype was determined based on acquired mutations in BAP1, PBRM1, SETD2 and VHL. Aggressive subtype was defined clinically using Mayo SSIGN score and molecularly using the ccA/ccB gene expression subtype. Publically-available Hi-C data were used to link germline risk variants with candidate target genes.
The 8q24 variant rs35252396 was significantly associated with VHL mutation status (OR = 1.6, p = 0.0037) and SSIGN score (OR = 1.9, p = 0.00094), after adjusting for multiple comparisons. We observed that, while some germline variants have interactions with nearby genes, some variants demonstrate long-range interactions with target genes.
These data further demonstrate the link between rs35252396, HIF pathway and ccRCC clinical aggressiveness, providing a more comprehensive picture of how germline genetics and tumor genetics interact with respect to tumor development and progression.
The majority (> 90%) of kidney cancer is classified as renal cell carcinoma (RCC) and approximately 85% of RCCs are further classified as the clear cell subtype (ccRCC). The etiology of ccRCC has been extensively studied and smoking, obesity and hypertension are recognized environmental risk factors that increase the risk of developing ccRCC. Additionally, genome-wide association studies (GWAS) have to date identified 14 germline variants that are associated with risk of RCC [1,2,3,4,5]. The value of these germline genetic explorations notwithstanding, the functional impact of the germline variants associated with RCC and ccRCC specifically remains largely unknown. Furthermore, the association of germline genetics with tumor genetics and tumor aggressiveness are largely unknown. In some cancers investigators have reported that germline variants are associated with specific molecularly-defined tumor subtypes, and in some cases the association is large enough to suggest clinical relevance (e.g., rs55705857 has an odds ratio > 6 in IDH-mutated glioma) [6,7,8]. To date, similar analyses linking germline variants with tumor subtypes have not been performed for RCC or ccRCC specifically. The value of identifying associations of genetic variants with specific molecular subtypes of a tumor centers on the ability to provide evidence of the involvement of specific developmental pathways that can help inform the biology of ccRCC tumor evolution, progression and novel prevention efforts. Related to this, we now have access to a catalog of acquired tumor alterations that are commonly present in ccRCC [9, 10]. Motivated by the opportunity to combine data on germline genetics associated with ccRCC risk with specific acquired molecular alterations found in ccRCC tumor tissue, we advance the field by evaluating for the first time the association of known ccRCC germline variants with these acquired alterations in order to better understand ccRCC development and progression. Moreover, we also evaluated the association of these known germline variants with well-known and validated clinical measures of ccRCC aggressiveness. Finally, we leverage Hi-C data to identify candidate target genes for the 14 germline variants that are the focus of this investigation.
The Cancer Genome Atlas (TCGA) data
Raw genotyping data from Affymetrix 6.0 array were obtained from germline DNA for 420 ccRCC patients, as well as corresponding clinical and pathological data. Somatic mutations of the four most frequently-mutated genes in ccRCC (BAP1, PBRM1, SETD2 and VHL) were obtained from TCGA . The clear cell A and clear cell B (ccA/ccB) gene expression subtype classification has been reproducibly shown to be associated with outcome [11, 12] and was obtained from . The Mayo SSIGN score has also been reproducibly shown to be associated with outcome [13,14,15,16]. The Mayo SSIGN score is derived from an additive model that contains tumor stage, tumor size, tumor grade and presence of necrosis, and was calculated as described previously .
Quality control was performed on the genotying data for the 420 ccRCC TCGA patients, including 95% call rate (zero germline variants failed), 95% sample call rate (zero samples failed), Hardy–Weinberg equilibrium (821 variants failed with p < 0.000001), minor allele frequency (MAF; 216,942 variants had MAF < 0.05), sex check (100% concordance) and population stratification (all Caucasian). Genotyping data were phased and imputed using the Michigan Imputation Server with the Haplotype Reference Consortium (release 1) as the reference population. The imputation quality for each of the 14 variants is provided in Additional file 1: Table S1. An additive logistic regression model was used to assess the association between each of the 14 variants and subtype, with genotype coded as having 0, 1, or 2 copies of the minor allele for observed data and dosage was modeled as continuous for imputed data. A general linear model was used to assess the association between age of diagnosis and each of the 14 variants. To account for multiple testing, p values < 0.004 were considered statistically significant (0.05/14 = 0.004).
Associations between known RCC germline variants and genes within two-to-five megabase (Mb) were evaluated using Hi-C data via the HUGIn web browser . Additional genetic information, limited to what was available in HUGIn, was included such as frequently interacting regions (FIREs), topologically associating domain (TAD) boundary regions, and occupancy of histone marks H3K27ac, H3K4me1 and H3K4me3. Because RCC is thought to originate from the renal tubular epithelium, and mesodermal stem cells form the tubule of the kidney, analyses were performed using mesendoderm cells. Mesenchymal stem cells were also analyzed ; results were similar across the two cell lines. We also used Hi-C data from ccRCC cell line (VHL mutant) Caki2 (GSM2827127 and GSM2827128) , (VHL wild-type) embryonic kidney cell lines HEK293T (GSM1081530 and GSM1081531) and HEK293T RAD21cv that was treated with tobacco etch mosaic virus protease (GSM1081526 and GSM1081527) . RAD21 is a core subunit of cohesin complex, which is known to play a role in mediating chromosomal loops. In RAD21cv (a RAD21-EGFP variant) cells, RAD21cv replaced endogenous RAD21 and was incorporated into the cohesin complex. For Caki2 and HEK293 data, reads were mapped with Bowtie 2  and alignments from two replicates were combined. Chromatin interaction was identified with HOMER (https://homer.ucsd.edu/homer/interactions/) at 20-kb resolution, which takes into account the dependence of interaction frequency and linear distance along each chromosome. For each risk locus, we combined the virtual 4C plot generated by the HUGIn web browser with the interaction plot generated from internally-analyzed Hi-C data.
TCGA ccRCC cohort
Table 1 describes the 420 TCGA ccRCC patients that were analyzed. Of the 376 patients with available whole exome sequencing data, 150 (40%) had a VHL mutation, 34 (9%) BAP1 mutation, 119 (32%) PBRM1 mutation and 48 (13%) SETD2 mutation. We also subtyped patients according to disease aggressiveness using pathological indices defined by the Mayo SSIGN score  as well as molecularly according to the ccA/ccB gene expression subtype [11, 12]. Of the 355 patients that had available pathology data to calculate the Mayo SSIGN score, 81 (23%) were classified as aggressive (SSIGN score > 8). Using ccA/ccB to classify aggressiveness, of the 352 patients who had data, 168 (48%) were poor prognosis (ccB) subtype.
Association of RCC germline variants with frequently mutated genes
Using a case-case design, we evaluated the association of each of the 14 RCC germline variants with known ccRCC acquired alterations in BAP1, PBRM1, SETD2 and VHL (Table 2). We observed a statistically significant association after adjusting for multiple comparisons between the 8q24 variant rs35252396 and VHL mutation (OR = 1.60, p = 0.0037). While not significant after adjusting for multiple testing, we also observed a candidate association between EPAS1 variant rs7579899 and SETD2 mutation (OR = 1.87, p = 0.012) (Table 2).
Association of RCC germline variants with aggressive ccRCC
We observed a statistically significant association between the 8q24 variant rs35252396 and Mayo SSIGN score (OR = 1.92, p = 0.00094) (Table 2). However, we did not observe a statistically-significant association between the known germline variants and ccA/ccB gene expression subtype.
Association of RCC germline variants with age at diagnosis
We did not observe a statistically-significant association between the known germline variants and age at diagnosis (Additional file 1: Table S2).
Association of RCC germline variants with nearby genes
We evaluated the interaction of each of the known germline variants with putative target genes in mesendoderm cell lines and mesenchymal stem cells using publically-available Hi-C data. Hi-C identifies chromatin interactions to evaluate the three dimensional chromatin structures inside the nucleus, which may identify long-range interactions. Some germline variants demonstrated interactions with nearby genes: e.g., rs4381241 with FAF1 (518.5 kb away), rs57579899 with EPAS1 (16.8 kb away), rs12105918 with ZEB2/ZEB2-AS1 (~ 70 kb away), rs1800057 with ATM (50.2 kb away) and rs4903064 with DPF3 (81.4 kb away) (Fig. 1; Additional file 1: Figure S1). However, some of these variants showed additional long-range interactions that have not been reported previously: e.g., rs4381241 with CDKN2C and TTC39A, rs7579899 with PRKCE, rs12105918 with ARHGAP15, GTDC1 and TEX41, rs1800057 with CUL5, ACAT1, NPAT and EXPH5 and rs4903064 with RGS6. Other germline variants also demonstrated long-range interactions that have not been reported previously: e.g., rs10936602 with SEC62 and PHC3, rs67311347 with ENTPD3, CTNNB1 and ULK4, rs2241261 with PEBP4 and EGR3, rs74911261 with EXPH5, rs718314 with SSPN and ITPR2 and rs4765623 with FAM101A (Additional file 1: Figure S1). Additionally, rs35252396 demonstrated interactions with PCAT1 and PCAT2 in mesenchymal stem cells but not in mesendoderm cells (Fig. 2). Similarly, rs11813268 had interactions with OBFC1 in mesenchymal stem cells but not in mesendoderm cells. Finally, we failed to identify genes within ± 1 Mb whose promoters interacted with rs7105934 in either mesenchymal stem cells or mesendoderm cells. All identified interactions were further evaluated in independent cell lines: VHL-mutant ccRCC cell line (Caki) and VHL wild-type embryonic kidney cell lines (HEK293 and HEK293 RAD21cv). A large proportion of the interactions were also identified in ccRCC or embryonic kidney cell lines. For example, rs10936602 interaction with SEC62, rs2241261 interactions with PEBP4 and EGR3, rs1800057 and rs74911261 interactions with EXPH5, as well as rs718314 interaction with SSPN were confirmed in all three cell lines. Further, rs12105918 interaction with GTDC1, rs11813268 interaction with OBFC1 and rs4765623 interaction with FAM101A were only identified in the ccRCC cell line, while rs7579899 interaction with PRKCE, rs35252396 interaction with PCAT2, and rs4903064 interaction with DPF3 were only identified in the embryonic kidney cell lines.
There are currently 14 known germline variants that are associated with risk of RCC; however, it remains unclear how germline genetics modify the risk of developing RCC or their association with tumor aggressiveness. To date, associations between germlines genetics and tumor genetics is largely unknown for RCC. With respect to tumor aggressiveness, previous investigators did not observe significant associations between a polygenic risk model that was derived from these 14 germline variants and age at onset or tumor stage . However, the analyses were performed on overall RCC and not within relevant RCC subtypes. Herein, we not only evaluated ccRCC specifically, but also the association of these 14 variants with specific molecular (BAP1, PBRM1, SETD2 and VHL) and clinically-aggressive subtypes of ccRCC. Using a case-case analysis, we observed that the 8q24 germline variant rs35252396 was significantly associated with tumor VHL mutation status as well as with the Mayo SSIGN score. Notably, the Mayo SSIGN score has been reproducibly shown to be associated with clinical outcome [13,14,15,16]. The 8q24 germline variant is located within PVT1, a candidate oncogene that is thought to regulate MYC to promote tumor formation. Of note, the 8q24 variant also overlaps a DNase I hypersensitive site and H3K4me1 peak from fetal kidney , indicating its location within a regulatory region. Through functional laboratory studies, investigators recently demonstrated that the 8q24 germline variant affects HIF binding to a MYC enhancer . While the EPAS1 (rs7579899) and CCND1 (rs7105934) variants are also linked to the HIF pathway, we did not observe a statistically-significant association between these two variants and VHL mutation. Additionally, while it did not pass our multiple testing significance threshold, we also observed a candidate association between the EPAS1 germline variant rs57579899 and SETD2 tumor mutation (p = 0.012). We previously reported that loss of SETD2 activity was associated with greater risk of ccRCC death .
While previous investigators did not observe a significant association between a RCC-derived polygenic risk model and tumor stage , we observed a significant association between the 8q24 germline variant and the Mayo SSIGN score. The difference could be due to the fact that the original analysis  was performed for overall RCC whereas our analyses were performed within a more homogeneous subtype of RCC, particularly, ccRCC. The Mayo SSIGN score is derived from an additive model that contains tumor stage, tumor size, tumor grade and presence of necrosis and the model has reproducibly been shown to be associated with outcome in ccRCC, with higher SSIGN score being associated with poorer prognosis [13,14,15,16]. While we observed a significant association of the 8q24 germline variant and the Mayo SSIGN score, we did not observe a significant association between any of the 14 variants and the molecularly-defined ccA/ccB expression subtype that has been linked to ccRCC aggressiveness [11, 12, 17].
Previous studies have performed expression quantitative trait loci (eQTL) analyses to evaluate function of the 14 RCC germline variants [2, 5]. Herein, we used Hi-C data to identify candidate target genes underlying the association of each of the 14 germline variants with ccRCC risk. While some variants demonstrated interactions with nearby genes, we observed additional long-range interactions. Laboratory studies are necessary to further understand these observations.
This study has limitations. We only analyzed associations of germline variants with somatic mutations, and no other acquired molecular alterations such as copy number variation. Additionally, because there are limited GWAS data available on patients who also have tumor molecular data, we did not validate the observed associations between the 8q24 germline variant and VHL tumor mutation nor the association between the EPAS1 germline variant and SETD2 tumor mutation. Similarly, Hi-C data are currently limited and thus we were not able to validate the Hi-C results. As such, the observed associations reported herein require validation in an external cohort.
We identified a significant association between the 8q24 germline variant and the presence of VHL somatic mutation. Additionally, we demonstrated for the first time an association between the 8q24 germline variant and ccRCC clinical aggressiveness as measured by the Mayo SSIGN score. Importantly, we additionally defined candidate target genes underlying the association between each of the 14 germline variants and risk of ccRCC. Together, these results further elucidate genes and pathways associated with ccRCC development. Specifically, these data further demonstrate the link between rs35252396 in the 8q24 region, HIF pathway and clinical aggressiveness, providing a more comprehensive biological understanding of the development of VHL mutated ccRCC tumors. Future work should evaluate how the rs35252396 germline variant and VHL mutation interact to affect treatment outcome and prognosis.
Availability of data and materials
The data that support the findings of this study are available from TCGA but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of TCGA.
Clear cell A
Clear cell B
Clear cell renal cell carcinoma
Expression quantitative trait loci
Frequently interacting regions
Genome wide association study
Minor allele frequency
- p :
Renal cell carcinoma
Stage, sign, grade and necrosis
Topologically associating domain
The Cancer Genome Atlas
Gudmundsson J, Sulem P, Gudbjartsson DF, et al. A common variant at 8q24.21 is associated with renal cell cancer. Nat Commun. 2013;4:2776.
Henrion M, Frampton M, Scelo G, et al. Common variation at 2q22.3 (ZEB2) influences the risk of renal cancer. Hum Mol Genet. 2013;22:825–31.
Wu X, Scelo G, Purdue MP, et al. A genome-wide association study identifies a novel susceptibility locus for renal cell carcinoma on 12p11.23. Hum Mol Genet. 2012;21:456–62.
Purdue MP, Johansson M, Zelenika D, et al. Genome-wide association study of renal cell carcinoma identifies two susceptibility loci on 2p21 and 11q13.3. Nat Genet. 2011;43:60–5.
Scelo G, Purdue MP, Brown KM, et al. Genome-wide association study identifies multiple risk loci for renal cell carcinoma. Nat Commun. 2017;8:15724.
Eckel-Passow JE, Lachance DH, Molinaro AM, et al. Glioma groups based on 1p/19q, IDH, and TERT promoter mutations in tumors. N Engl J Med. 2015;372:2499–508.
Kuchenbaecker KB, McGuffog L, Barrowdale D, et al. Evaluation of polygenic risk scores for breast and ovarian cancer risk prediction in BRCA1 and BRCA2 mutation carriers. J Natl Cancer Inst. 2017;109:djw302.
Lecarpentier J, Silvestri V, Kuchenbaecker KB, et al. Prediction of breast and prostate cancer risks in male BRCA1 and BRCA2 mutation carriers using polygenic risk scores. J Clin Oncol. 2017;35:2240–50.
Tcga TCGA. Comprehensive molecular characterization of clear cell renal cell carcinoma. Nature. 2013;499:43–9.
Sato Y, Yoshizato T, Shiraishi Y, et al. Integrated molecular analysis of clear-cell renal cell carcinoma. Nat Genet. 2013;45:860–7.
Brannon AR, Reddy A, Seiler M, et al. Molecular stratification of clear cell renal cell carcinoma by consensus clustering reveals distinct subtypes and survival patterns. Genes Cancer. 2010;1:152–63.
Brooks SA, Brannon AR, Parker JS, et al. ClearCode34: a prognostic risk predictor for localized clear cell renal cell carcinoma. Eur Urol. 2014;66:77–84.
Frank I, Blute ML, Cheville JC, et al. An outcome prediction model for patients with clear cell renal cell carcinoma treated with radical nephrectomy based on tumor stage, size, grade and necrosis: the SSIGN score. J Urol. 2002;168:2395–400.
Zigeuner R, Hutterer G, Chromecki T, et al. External validation of the Mayo Clinic stage, size, grade, and necrosis (SSIGN) score for clear-cell renal cell carcinoma in a single European centre applying routine pathology. Eur Urol. 2010;57:102–9.
Ficarra V, Novara G, Galfano A, et al. The “Stage, Size, Grade and Necrosis” score is more accurate than the University of California Los Angeles Integrated Staging System for predicting cancer-specific survival in patients with clear cell renal cell carcinoma. BJU Int. 2009;103:165–70.
Parker WP, Cheville JC, Frank I, et al. Application of the Stage, Size, Grade, and Necrosis (SSIGN) Score for clear cell renal cell carcinoma in contemporary patients. Eur Urol. 2017;71:665–73.
Eckel-Passow JE, Igel DA, Serie DJ, et al. Assessing the clinical use of clear cell renal cell carcinoma molecular subtypes identified by RNA expression analysis. Urol Oncol. 2015;33(68):e17-23.
Schmitt AD, Hu M, Jung I, et al. A compendium of chromatin contact maps reveals spatially active regions in the human genome. Cell Rep. 2016;17:2042–59.
Myszczyszyn A, Czarnecka AM, Matak D, et al. The role of hypoxia and cancer stem cells in renal cell carcinoma pathogenesis. Stem Cell Rev. 2015;11:919–43.
ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489:57–74.
Zuin J, Dixon JR, van der Reijden MI, et al. Cohesin and CTCF differentially affect chromatin architecture and gene expression in human cells. Proc Natl Acad Sci USA. 2014;111:996–1001.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.
Kundaje A, Meuleman W, Ernst J, et al. Integrative analysis of 111 reference human epigenomes. Nature. 2015;518:317–30.
Grampp S, Platt JL, Lauer V, et al. Genetic variation at the 8q24.21 renal cancer susceptibility locus affects HIF binding to a MYC enhancer. Nat Commun. 2016;7:13183.
Ho TH, Kapur P, Joseph RW, et al. Loss of histone H3 lysine 36 trimethylation is associated with an increased risk of renal cell carcinoma-specific death. Mod Pathol. 2016;29:34–42.
The results published here are in whole or part based upon data generated by the TCGA Research Network: https://cancergenome.nih.gov/.
Ethics approval and consent to participate
Consent for publiation
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Hi-C interactions for each of the 14 known RCC germline variants. Chromosome 1: rs4381241. Supplementary Table S1. Information on the 4 known RCC germline variants. R-square denotes the imputation quality and FRQ denotes the frequency of the variants for each of the aquired alterations that were evaluated. Supplementary Table S2. Association of known RCC germline variants with age of diagnosis.
About this article
Cite this article
Eckel-Passow, J.E., Yan, H., Kosel, M.L. et al. 8q24 clear cell renal cell carcinoma germline variant is associated with VHL mutation status and clinical aggressiveness. BMC Urol 20, 173 (2020). https://doi.org/10.1186/s12894-020-00745-9