8q24 Clear Cell Renal Cell Carcinoma Germline Variant is Associated with VHL Mutation Status and Clinical Aggressiveness

doi:10.21203/rs.3.rs-48012/v2

Background

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.

Methods

We analyzed 420 ccRCC patients from The Cancer Genome Atlas (TCGA). 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.

Results

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.

Conclusions

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.

Urology & Nephrology

kidney

GWAS

subtype

Hi-C

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-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-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 [9]. 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 [9]. The Mayo SSIGN score has also been reproducibly shown to be associated with outcome [13-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 [17].

Statistical Methods

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 Supplementary 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).

Hi-C Analysis

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 [18]. 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 [19]; results were similar across the two cell lines. We also used Hi-C data from ccRCC cell line (VHL mutant) Caki2 (GSM2827127 and GSM2827128) [20], (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) [21]. 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 [22] and alignments from two replicates were combined. Chromatin interaction was identified with HOMER (http://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 [13] 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 (Supplementary 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) (Figure 1; Supplementary 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 (Supplementary Figure S1). Additionally, rs35252396 demonstrated interactions with PCAT1 and PCAT2 in mesenchymal stem cells but not in mesendoderm cells (Figure 2). Similarly, rs11813268 had interactions with OBFC1 in mesenchymal stem cells but not in mesendoderm cells. Finally, we failed to identify genes within ±1Mb 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 [5]. 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-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 [23], 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 [24]. 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 [25].

While previous investigators did not observe a significant association between a RCC-derived polygenic risk model and tumor stage [5], 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 [5] 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-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.

ccA clear cell A

ccB clear cell B

ccRCC clear cell renal cell carcinoma

eQTL expression quantitative trait loci

FIREs frequently interacting regions

GWAS genome wide association study

kb kilobase

MAF minor allele frequency

OR odds ratio

p pvalue

RCC renal cell carcinoma

SSIGN stage, sign, grade and necrosis

TAD topologically associating domain

TCGA The Cancer Genome Atlas

- Ethics approval and consent to participate

Not applicable.

- Consent to publish

Not applicable.

- 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.

- Competing interests

The authors have no competing interests.

- Funding

Not applicable.

- Authors' Contributions

The study was designed by JEEP and AP. Statistical and bioinformatics analyses were performed by HY, MK, DS and PD. Results were interpreted by all authors and all authors contributed to the final version of the paper.

- Acknowledgements

The results published here are in whole or part based upon data generated by the TCGA Research Network: http://cancergenome.nih.gov/.

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-831.
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-462.
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-65.
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-2508.
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.
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-2250.
TCGA TCGA. Comprehensive molecular characterization of clear cell renal cell carcinoma. Nature 2013; 499: 43-49.
Sato Y, Yoshizato T, Shiraishi Y, et al. Integrated molecular analysis of clear-cell renal cell carcinoma. Nat Genet 2013; 45: 860-867.
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-163.
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-2400.
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-109.
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-170.
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-673.
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-2059.
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-943.
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 U S A 2014; 111: 996-1001.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods 2012; 9: 357-359.
Kundaje A, Meuleman W, Ernst J, et al. Integrative analysis of 111 reference human epigenomes. Nature 2015; 518: 317-330.
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.

Table 1: Description of 420 TCGA ccRCC patients.
		N (%)
Sex
Female		137 (32.6%)
Male		283 (67.4%)

Max Tumor Size
N		376
Mean (SD)		6.5 (3.5)
Median		5.5
Q1, Q3		4.0, 8.5
Range		(1.1-25.0)

Stage
Missing		44
Stage I		176 (46.8%)
Stage II		135 (35.9%)
Stage IV		65 (17.3%)

grade
Missing		44
G1		7 (1.9%)
G2		156 (41.5%)
G3		152 (40.4%)
G4		60 (16.0%)
GX		1 (0.3%)

Percent Necrosis
Missing		64
0		189 (53.1%)
2		4 (1.1%)
3		2 (0.6%)
5		42 (11.8%)
8		1 (0.3%)
10		23 (6.5%)
15		13 (3.7%)
20		6 (1.7%)
25		4 (1.1%)
30		72 (20.2%)

SSIGN Group
Missing		65
Low Risk (0-3)		159 (44.8%)
Intermediate Risk (4-7)		115 (32.4%)
High Risk (8+)		81 (22.8%)

BAP1 Mutated
Missing		44
No		342 (91.0%)
Yes		34 (9.0%)

PBRM1 Mutated
Missing		44
No		257 (68.4%)
Yes		119 (31.6%)

SETD2 Mutated
Missing		44
No		328 (87.2%)
Yes		48 (12.8%)

VHL Mutated
Missing		44
No		226 (60.1%)
Yes		150 (39.9%)

ccA/ccB Expression Subtype
Missing		68
ccB		168 (47.7%)
ccA		184 (52.3%)

SupplementaryMaterials20200908.docx

8q24 Clear Cell Renal Cell Carcinoma Germline Variant is Associated with VHL Mutation Status and Clinical Aggressiveness

Status:

Journal Publication

Version 2

Abstract

Figures

Background

Methods

Results

Discussion

Conclusion

Abbreviations

Declarations

References

Tables

Supplementary Files

Status:

Journal Publication

Version 2