Search

Genome-wide association studies and candidate gene studies have identified several genetic variants that increase glioma risk. The majority of these variants are non-coding and the mechanisms behind the increased risk in carriers are not known. In this study, we hypothesize that some of the established glioma risk variants induce aberrant DNA methylation in the developing tumor, either locally (gene-specific) or globally (genome-wide). In a pilot data set including 77 glioma patients, we used Illumina beadchip technology to analyze genetic variants in blood and DNA methylation in matched tumor samples. To validate our findings, we used data from the Cancer Genome Atlas, including 401 glioblastoma patients. Consensus clustering identified the glioma CpG island methylator phenotype (gCIMP) and two additional subgroups with distinct patterns of global DNA methylation. In the pilot dataset, gCIMP was associated with two genetic variants in CDKN2B-AS1, rs1412829 and rs4977756 (9p21.3, p = 8.1 x 10-7 and 4.8 x 10-5, respectively). The association was in the same direction in the TCGA dataset, although statistically significant only when combining individuals with AG and GG genotypes. We also investigated the relation between glioma risk variants and DNA methylation in the promoter region of genes located within 30 kb of each variant. One association in the pilot dataset, between the TERT risk variant rs2736100 and lower methylation of cg23827991 (in TERT; p = 0.001), was confirmed in the TCGA dataset (p = 0.001). In conclusion, we found an association between rs1412829 and rs4977756 (9p21.3, CDKN2B-AS1) and global DNA methylation pattern in glioma, for which a trend was seen also in the TCGA glioblastoma dataset. We also found an association between rs2736100 (in TERT) and levels of methylation at cg23827991 (localized in the same gene, 3.3 kbp downstream of the risk variant), which was validated in the TCGA dataset. Except for this one association, we did not find strong evidence for gene-specific DNA methylation mediated by glioma risk variants., Competing Interests: The authors have declared that no competing interests exist.

During the last years, genome wide association studies have discovered common germline genetic variants associated with specific glioma subtypes. We aimed to study the association between these germline risk variants and tumor phenotypes, including copy number aberrations and protein expression. A total of 91 glioma patients were included. Thirteen well known genetic risk variants in TERT, EGFR, CCDC26, CDKN2A, CDKN2B, PHLDB1, TP53, and RTEL1 were selected for investigation of possible correlations with the glioma somatic markers: EGFR amplification, 1p/19q codeletion and protein expression of p53, Ki-67, and mutated IDH1. The CDKN2A/B risk variant, rs4977756, and the CDKN2B risk variant, rs1412829 were inversely associated (p = 0.049 and p = 0.002, respectively) with absence of a mutated IDH1, i.e., the majority of patients homozygous for the risk allele showed no or low expression of mutated IDH1. The RTEL1 risk variant, rs6010620 was associated (p = 0.013) with not having 1p/19q codeletion, i.e., the majority of patients homozygous for the risk allele did not show 1p/19q codeletion. In addition, the EGFR risk variant rs17172430 and the CDKN2B risk variant rs1412829, both showed a trend for association (p = 0.055 and p = 0.051, respectively) with increased EGFR copy number, i.e., the majority of patients homozygote for the risk alleles showed chromosomal gain or amplification of EGFR. Our findings indicate that CDKN2A/B risk genotypes are associated with primary glioblastoma without IDH mutation, and that there is an inverse association between RTEL1 risk genotypes and 1p/19q codeletion, suggesting that these genetic variants have a molecular impact on the genesis of high graded brain tumors. Further experimental studies are needed to delineate the functional mechanism of the association between genotype and somatic genetic aberrations.

A number of gene variants have been associated with an increased risk of developing glioma. We hypothesized that the reported risk variants may be associated with tumor genomic instability. To explore potential correlations between germline risk variants and somatic genetic events, we analyzed matched tumor and blood samples from 95 glioma patients by means of SNP genotyping. The generated genotype data was used to calculate genome-wide allele-specific copy number profiles of the tumor samples. We compared the copy number profiles across samples and found two EGFR gene variants (rs17172430 and rs11979158) that were associated with homozygous deletion at the CDKN2A/B locus. One of the EGFR variants (rs17172430) was also associated with loss of heterozygosity at the EGFR locus. Our findings were confirmed in a separate dataset consisting of matched blood and tumor samples from 300 glioblastoma patients, compiled from publically available TCGA data. These results imply there is a functional effect of germline EGFR variants on tumor progression.