59
Views
3
CrossRef citations to date
0
Altmetric
Original Research

An intronic genetic variation of MGMT affects enhancer activity and is associated with glioma susceptibility

, , , , &
Pages 3995-4003 | Published online: 27 Sep 2018

Abstract

Purpose

O6-methylguanine-DNA methyltransferase (MGMT) plays a crucial role in repairing damaged DNA caused by alkylating agents. A number of cancer susceptibility loci have been recognized as enhancer variants. This study aimed to explore the significance of enhancer variants of MGMT in glioma susceptibility.

Patients and methods

A retrospective case-control study consisting of 150 glioma patients and 327 controls was conducted to test whether enhancer variants of MGMT are associated with glioma susceptibility. Genotypes were determined by Sequenom MassARRAY technology. Associations were estimated by logistic regression. Biochemical assays were used to examine the function of glioma susceptibility locus.

Results

We found that the A allele of rs10764901, an intronic variant of MGMT, was associated with a significantly decreased risk of glioma. The rs10764901 AA genotype carriers had an OR of 0.49 (95% CI, 0.24–0.98; P=0.045) compared with the rs10764901 GG genotype. When the rs10764901 AG and AA genotypes were pooled for analysis, a significantly decreased risk of glioma was also found (OR, 0.63; 95% CI, 0.43–0.93; P=0.021). Functional analyses showed that the rs10764901 A allele drove a lower luciferase expression and had higher transcription factor binding affinity than the G allele.

Conclusion

An enhancer variant of MGMT rs10764901 affects the regulatory activity of enhancer by altering the binding affinity of transcription factors and is associated with glioma susceptibility.

Introduction

O6-methylguanine-DNA methyltransferase (MGMT) is an important DNA repair enzyme that plays a crucial role in protection against the carcinogenic effects of alkylating agents.Citation1 MGMT transfers mutagenic and cytotoxic alkyl groups from DNA to its active site and then loses activity irreversibly.Citation2 Therefore, a certain amount of MGMT is required to maintain its DNA repair activity. Previous studies have shown that aberrant expression of MGMT plays an important role in carcinogenesis.Citation3Citation5 The increased frequency of oncogene mutations might be the pivotal carcinogenic mechanism of MGMT inactivation.Citation6,Citation7 On the other hand, the MGMT level was also reported to be associated with response to alkylating drugs in cancer.Citation8,Citation9 Promoter methylation, which is the most frequent mechanism of inactivation of MGMT in glioma,Citation10 has been well established as a major prognostic factor and predictor of chemotherapy response for glioma.Citation11,Citation12 However, the expression level of MGMT is not exactly matched to the promoter methylation status,Citation13,Citation14 which indicates that other mechanisms also play a non-negligible role in MGMT expression regulation.

Some studies have revealed that several genetic variations of MGMT affect gene expression and then play an important role in carcinogenesis and progression of several cancers.Citation15Citation17 For example, the T allele of a promoter genetic variation rs16906252 was found to reduce expression of MGMT and be associated with increased risk of MGMT-methylated colorectal cancer and improved survival of patients with glioma.Citation15,Citation17 In recent years, some novel significant genetic susceptibility factors of cancers have been identified by unbiased genome-wide association studies (GWAS). It is worthy to note that most of them are located in non-coding regions of genome.Citation18 Furthermore, a number of non-coding genetic variations identified by GWAS occur in regions that act as enhancers and have functional effects on enhancer activity.Citation18Citation20

In order to survey the role of enhancer variants of MGMT in glioma, we investigated their associations with glioma susceptibility and explored their functional relevance. Our data demonstrate that an intronic variant of MGMT rs10764901 affects the regulatory activity of an enhancer by altering the binding affinity of transcription factors and is associated with glioma susceptibility.

Patients and methods

Selection of candidate genetic variations

Several tracks of Ensembl Genome Browser (http://www.ensembl.org/index.html, Release 75) were used to select candidate genetic variations.Citation21 First, we used the VISTA track to explore potential enhancers of MGMT.Citation22 Not only the gene body of MGMT but also its 100 kb upstream and downstream flanking regions were analyzed. Then, the Genomic Evolutionary Rate Profiling (GERP) conservation scores based on the 37 eutherian mammals were used to measure the evolutionary conservation of sequences of the potential enhancers of MGMT.Citation23 In order to screen the common genetic variations in our study population, the 1000 Genomes Project East Asian (1 KG EAS) common short variants (single-nucleotide polymorphisms [SNPs] and indels) track, which displays the variants with minor allelic frequency (MAF) >1% genotyped in East Asian individuals by the 1000 Genomes project, was used to analyze the variants located in the potential enhancers of MGMT.Citation24 Because of the relatively small sample size, only the common variants with MAF >10% in Southern Han Chinese (CHS) population were selected for genotyping in this study.

Study subjects

This study recruited 150 incident glioma patients and 327 healthy controls. A portion of patients and controls was enrolled in our previous studies on glioma.Citation25,Citation26 In the present study, we extended the sample size of glioma patients from 138 to 150. All subjects were unrelated CHS living in Fujian, Zhejiang, and surrounding provinces. Patients with histopathologically confirmed glioma were recruited from January 2010 to October 2014 at the First Affiliated Hospital, Fujian Medical University (Fuzhou, n=122) and the First Affiliated Hospital, College of Medicine, Zhejiang University (Hangzhou, n=28). Healthy controls were cancer-free individuals recruited during the same period. At recruitment, written informed consent was obtained from each subject, and the detailed information on demographic and clinical characteristics were collected. This study was approved by the Institutional Review Board of the First Affiliated Hospital, Fujian Medical University.

Genotype analysis

Genomic DNA samples of 138 glioma patients and 327 healthy controls enrolled in our previous studies were isolated from the peripheral blood lymphocytes using a commercial Tiangen TIANamp Genomic DNA kit (Tiangen Biotech., Beijing, China).Citation25,Citation26 To increase the sample size of cases, 8.0% (12/150) of genomic DNA samples of cases were isolated from paraffin-embedded normal tissue adjacent to cancer specimens using KAPA Express Extract Kits (KAPA Biosys-tems, Wilmington, MA, USA). Genotypes of the candidate enhancer variants were determined by the Sequenom Mas-sARRAY iPLEX platform (Sequenom Inc., San Diego, CA, USA). The information of primers is shown in . To control the data quality, genotyping was performed without knowledge of the case/control status of the subjects. And a 5% random sample was tested twice by different persons, and the reproducibility was 100%.

Functional analyses of genetic variation

In the present study, we performed functional studies only for rs10764901 which was identified to be associated with glioma susceptibility significantly.

For luciferase reporter gene assays, two 1,366-bp DNA fragments corresponding to the potential enhancer (hs589) and containing rs10764901 G or A allele were amplified by polymerase chain reaction using same DNA sample template (primers are available upon request) and subcloned into the pGL3-promoter vector (Promega Corporation, Fitchburg, WI, USA). The resultant plasmid that contains rs10764901 G allele was designated as P-G, while the other one was designated as P-A. The two constructs were identical except for the different allele at rs10764901 polymorphic site. They were restriction mapped and sequenced to confirm their authenticity. U251 cell was used for luciferase assays. The constructed reporter plasmid or the blank pGL3-promoter plasmid was co-electrotransfected with pRL-TK (Promega Corporation) to U251 cell, using Celetrix Electroporator (Celetrix Biotech., Manassas, VA, USA). Three independent transfection experiments were done, and each was performed in triplicate. The luciferase activity was analyzed by a Dual-Luciferase Reporter Assay System (Promega Corporation).

For electrophoretic mobility shift assays (EMSA), synthetic double-stranded and 3’ biotin-labeled oligonucleotides corresponding to rs10764901 G (Probe-G) or A (Probe-A) sequences and U251 cell nuclear extract were incubated at 25°C for 20 minutes using the Light Shift Chemiluminescent EMSA Kit (Pierce, Rockford, IL, USA). For competition assays, non-labeled oligonucleotides at 150-fold molar excess were added to the reaction mixture before the addition of biotin-labeled probes. The reaction mixture was separated on 7% polyacrylamide gel electrophoresis, and the products were detected by Stabilized Streptavidin-Horseradish Peroxidase Conjugate (Pierce). Densitometric analyses of EMSA images were performed by ImageJ software (version 1.50i, NIH). The information of probes is shown in .

Statistical analysis

The differences in sex and age between cases and controls were examined by two-sided chi-squared test and Wilcoxon rank sum test, respectively. The associations of the candidate enhancer variants with glioma susceptibility were estimated by ORs and their 95% CIs computed by multivariate logistic regression model and adjusted for sex and age. Bonferroni correction was used to adjust the P-value for multiple testing. Meanwhile, the false-positive report probability (FPRP) was also calculated to assess the significance as described by Wacholder et al.Citation27 Because of the relatively small sample size, we set the significance value of FPRP at 0.50 and the prior probability at 0.10 in this study. The statistical power to detect an OR of 1.50 (or its reciprocal, 0.67) was used to calculate FPRP value. We also performed the subgroup analyses by sex, age, and WHO grade. A t-test was used to examine the differences in luciferase activity between P-G, P-A, and the blank pGL3-promoter. The differences in densitometries of DNA-protein complexes formed with Probe-A and Probe-G were also analyzed by t-test. These statistical analyses were implemented in Statistic Analysis System software (version 9.4, SAS Institute, Cary, NC, USA) and STATA statistical software (version 14.0; StataCorp LP, College Station, TX, USA). P-value of <0.05 was used as the criterion of statistical significance, and all statistical tests were two sided.

Results

We first used Ensembl Genome Browser to choose candidate variants located in the potential enhancers of MGMT. As shown in , four potential enhancers (hs656, hs696, hs331, and hs589) located in the intron of MGMT were identified by VISTA track. The potential enhancers hs656 and hs331 were excluded from further analyses because of the absence of common variants with MAF >10% in CHS population. The remaining two potential enhancers hs696 and hs589 were shown to have extreme evolutionary sequence conservation by the GERP conservation scores track. Two common variants rs577227 (MAF=0.39) and rs10764901 (MAF=0.41) were found in the potential enhancers hs696 and hs589, respectively. They were both selected as candidate enhancer variants for further genotyping.

Figure 1 Selection of candidate enhancer variants of MGMT.

Notes: The VISTA track of Ensembl Genome Browser displays the potential enhancers of MGMT. The GERP scores represent the evolutionary conservation of sequences. The 1 KG EAS common short variants (SNPs and indels) track displays the variants with frequency >1% in East Asian individuals.

Abbreviations: MGMT, O6-methylguanine-DNA methyltransferase; SNP, single-nucleotide polymorphism. GERP, Genomic Evolutionary Rate Profiling; 1KG EAS, 1000 Genomes Project East Asian.

Figure 1 Selection of candidate enhancer variants of MGMT.Notes: The VISTA track of Ensembl Genome Browser displays the potential enhancers of MGMT. The GERP scores represent the evolutionary conservation of sequences. The 1 KG EAS common short variants (SNPs and indels) track displays the variants with frequency >1% in East Asian individuals.Abbreviations: MGMT, O6-methylguanine-DNA methyltransferase; SNP, single-nucleotide polymorphism. GERP, Genomic Evolutionary Rate Profiling; 1KG EAS, 1000 Genomes Project East Asian.

The detailed baseline characteristics of study subjects are presented in . There was no significant differences of age and sex distributions between glioma patients and controls, which suggested that the frequency matching of cases and controls was adequate. Of the 150 patients, 80 (53.3%) patients were classified into astrocytic tumors, 37 (24.7%) were classified into oligodendroglial tumors, 16 (10.7%) were classified into oligoastrocytic tumors, and 17 (11.3%) were other histological types. The detailed tumor WHO classification data were successfully collected from 136 (90.7%) patients. Among them, seven (4.7%) patients were WHO grade I, 45 (30.0%) were WHO grade II, 37 (24.7%) were WHO grade III, and 47 (31.3%) were WHO grade IV. There are no clear grade records for the remaining 14 (9.3%) patients. Maximal safe resection or subtotal resection were performed on all patients. Meanwhile, 101 (67.3%) and 112 (74.7%) patients also underwent radiotherapy- and alkylate-based chemotherapy, respectively.

Table 1 Selected characteristics of glioma patients and controls

The genotype distributions of rs10764901 and rs577227 in controls were in agreement with that expected under Hardy-Weinberg equilibrium. The multivariate logistic regression was used to investigate the associations between the candidate enhancer variants and glioma susceptibility with adjustment for sex and age (). As a result, a significant association between rs10764901 and glioma susceptibility was found. Subjects carrying rs10764901 AG or AA genotype had an adjusted OR of 0.67 (95% CI, 0.44–1.02; P=0.059) or 0.49 (95% CI, 0.24–0.98; P=0.045), respectively, compared with individuals carrying GG genotype. These results suggest that rs10764901 A allele may be a protective allele and acts in an allele dose-dependent manner (trend test; P=0.018). A more significant decreased risk of glioma was seen when groups of rs10764901 AG and AA genotype were pooled for analysis (adjusted OR, 0.63; 95% CI, 0.43–0.93; P=0.021). The association for rs10764901 was still significant after Bonferroni correction for two testings. Moreover, the estimated FPRP value for rs10764901 was 0.332 which is below the prespeci-fied FPRP value of 0.50. It also suggests that our finding is noteworthy. We also examined whether genotype frequencies of rs10764901 differed between blood DNA and adjacent normal tissue DNA. No significant difference was found (). The significant association between rs10764901 and glioma susceptibility was also found when we repeated the association analyses without the 12 cases with adjacent normal tissue DNA (). These results verified that there was no significant heterogeneity in genotype frequencies for rs10764901 between different DNA sources. The other candidate enhancer variant rs577227 was not associated with glioma susceptibility significantly. Compared with subjects carrying rs577227 TT genotype, the adjusted ORs for TC or CC genotype were 1.33 (95% CI, 0.87–2.03; P=0.193) or 1.32 (95% CI, 0.69–2.53; P=0.410), respectively. When we combined rs577227 TC and CC genotype for analysis, the adjusted OR was 1.34 (95% CI, 0.89–2.01; P=0.156).

Table 2 Genotype frequencies of rs10764901 and rs577227 in glioma patients and controls and their contributions to the risk of glioma

We then carried out subgroup analyses for rs10764901 based on sex, age, and WHO grade in a dominant model (). The protective effect of rs10764901 A allele was observed across all subgroups, especially in males and patients with WHO grade III glioma. Significantly decreased risk of glioma was observed in males (adjusted OR, 0.53; 95% CI, 0.32–0.86; P=0.011). Similarly, it was found that rs10764901 was significantly associated with deceased risk of WHO grade III glioma (adjusted OR, 0.49; 95% CI, 0.25–0.99; P=0.046). No significant association was found in other subgroups although the trend is similar in all subgroups. It might be due to the small sample sizes within a stratum.

Figure 2 Subgroup analyses for rs10764901 based on sex, age, and WHO grade in a dominant model.

Note: The central black dot represents the OR, and the horizontal line indicates the 95% CI.

Figure 2 Subgroup analyses for rs10764901 based on sex, age, and WHO grade in a dominant model.Note: The central black dot represents the OR, and the horizontal line indicates the 95% CI.

To investigate the regulatory activity of the potential enhancer hs589 and whether rs10764901 has a functional effect on the regulatory activity of hs589, we constructed two luciferase reporter gene plasmids (P-G and P-A) encompassing hs589 with rs10764901 G or A allele, respectively, and then transiently transfected into U251 cells. As shown in , about six- or threefold higher reporter gene expression was observed for the plasmid P-G or P-A, respectively, than the pGL3-promoter plasmid (P-G vs pGL3-promoter, P<0.001; P-A vs pGL3-promoter, P<0.001). These results suggest that hs589 might be a real functional enhancer element. Then, we compared the transcriptional activity of rs10764901 G and A alleles. The P-G construct drove about twofold higher luciferase expression than the P-A construct (P=0.002). This finding indicates that rs10764901 G>A change affects the enhancer activity of hs589.

Figure 3 Results of luciferase reporter gene and EMSA.

Notes: (A) Luciferase expression of constructs containing the potential enhancer (hs589) with rs10764901 G allele (P-G) or A allele (P-A) in U251 cells. P-G, P-A, or pGL3-promoter construct was co-electrotransfected with pRL-TK to standardize transfection efficiency, respectively. Data shown are the means ± SE from three independent transfection experiments, each performed in triplicate. (B) EMSA with biotin-labeled oligonucleotides containing the rs10764901 G allele (Probe-G) or A allele (Probe-A) and nuclear extracts from U251 cells. Lanes 1 and 5 show mobilities of the labeled oligonucleotides without nuclear extracts; lanes 2 and 6 show mobilities of the labeled oligonucleotides with nuclear extracts in the presence of non-labeled G or non-labeled A competitors; lanes 3 and 7 show mobilities of the labeled oligonucleotides with nuclear extracts in the presence of non-labeled mut competitor; lanes 4 and 8 show mobilities of the labeled oligonucleotides with nuclear extracts in the absence of competitor.

Abbreviation: EMSA, electrophoretic mobility shift assay.

Figure 3 Results of luciferase reporter gene and EMSA.Notes: (A) Luciferase expression of constructs containing the potential enhancer (hs589) with rs10764901 G allele (P-G) or A allele (P-A) in U251 cells. P-G, P-A, or pGL3-promoter construct was co-electrotransfected with pRL-TK to standardize transfection efficiency, respectively. Data shown are the means ± SE from three independent transfection experiments, each performed in triplicate. (B) EMSA with biotin-labeled oligonucleotides containing the rs10764901 G allele (Probe-G) or A allele (Probe-A) and nuclear extracts from U251 cells. Lanes 1 and 5 show mobilities of the labeled oligonucleotides without nuclear extracts; lanes 2 and 6 show mobilities of the labeled oligonucleotides with nuclear extracts in the presence of non-labeled G or non-labeled A competitors; lanes 3 and 7 show mobilities of the labeled oligonucleotides with nuclear extracts in the presence of non-labeled mut competitor; lanes 4 and 8 show mobilities of the labeled oligonucleotides with nuclear extracts in the absence of competitor.Abbreviation: EMSA, electrophoretic mobility shift assay.

To explore whether the differences in enhancer activity between rs10764901 G and A alleles are due to their ability to bind transcriptional activators or inhibitors, we conducted EMSA using U251 cell nuclear extract. As shown in , U251 cell nuclear extract was able to bind both Probe-G and Probe-A, and the binding patterns were similar. However, more DNA-protein complex was formed when Probe-A was incubated with U251 cell nuclear extract (Lane 4) compared to Probe-G (Lane 8) under the same experimental conditions, with the relative densitometries (mean ± SE) from three independent experiments being 1.08±0.06vs 0.73±0.07 (P=0.016). It indicates that the protein binding affinity of Probe-A was higher than that of Probe-G. In competition assays, 150-fold excess of non-labeled A or non-labeled G probe could completely eliminate the DNA-protein interaction (Lane 2 and Lane 6, respectively). Meanwhile, the DNA–protein interaction was partly eliminated by 150-fold excess of non-labeled mut probe (Lane 3 and Lane 7). These results indicate that the binding is sequence specific.

Discussion

In the present study, we demonstrated that an intronic variant of MGMT rs10764901 is significantly associated with glioma susceptibility by a retrospective case-control study involving 150 glioma patients and 327 controls. The rs10764901 A allele carriers were found to have a significantly decreased risk of glioma. Functional analyses revealed that the rs10764901 A allele had higher transcription factor binding affinity than the G allele and drove a lower luciferase expression.

MGMT is well known as a unique DNA repair gene that reverses DNA alkylation damage alone.Citation1 It is widely considered to be a tumor suppressor gene because some studies have found that overexpression of MGMT is associated with decreased risk of several cancers.Citation3,Citation4 However, we found an interesting result that rs10764901 A allele, which drove a lower gene expression, significantly reduced the risk of glioma. It indicates that MGMT may function as an oncogene under certain circumstances. Consistent with our results, a study on colorectal cancer found that the T allele of MGMT promoter variant rs16906252 is associated with constitutively reduced gene expression and reduces the risk of MGMT-unmethylated colorectal cancer.Citation17 These results indicate that MGMT overexpression may be a risk factor for MGMT-unmethylated cancer. Interestingly, a converse result that rs16906252 T allele is associated with an elevated risk of MGMT-methylated colorectal cancer was observed.Citation17 These clues indicate that the role of MGMT in carcinogenesis is far more complex, and it might be completely opposite in cancer subclassifications with different methylation status of MGMT. Moreover, it was also reported that, in the context of different methylation status of MGMT, the prognostic effect of TERT gene is dichotomous,Citation28 which supports the hypothesis that the mechanism of carcinogenesis is different for cancer subclassifications with different methylation status of MGMT. This postulation can partly explain the inconclusive results reported by previous published studies on MGMT genetic variations.Citation29 The studies enrolling more patients with MGMT methylation might reach the opposite conclusion to those recruiting more MGMT-unmethylated patients.

Unfortunately, detailed information on methylation status of MGMT of glioma patients enrolled in this study is scant. Thus, we have no hard evidence to indicate that the MGMT-unmethylated glioma constituent ratio is higher than MGMT-methylated glioma in this study, which might partly explain why a deleterious role of MGMT overexpression in glioma carcinogenesis was inferred from our study. However, it is worthy to note that nearly two thirds of patients enrolled in this study are males. MGMT unmethylation has been found to be significantly higher in male patients with glioma.Citation30,Citation31 Therefore, it is plausible to deduce that MGMT-unmethylated glioma is in the majority in our study.

Still some other clues indicate that the physiological role of MGMT is far more than DNA repair. For instance, MGMT is expressed ubiquitously and its expression level far exceed its need to repair DNA damage.Citation1 The exact role of phosphorylated MGMT is still an enigma.Citation1 All of these and our findings indicate a broad function of MGMT needed further investigation.

A major limitation of this study is the relatively small sample size which might limit the statistical power. In addition, the lack of information on methylation status of MGMT results in the absence of subgroup analyses based on methylation status. Therefore, we could not explore directly whether rs10764901 has a converse role in developing glioma with different methylation status of MGMT.

Conclusion

We identified an enhancer variant of MGMT rs10764901 as a novel susceptibility locus for glioma for the first time. Our data revealed that rs10764901 A allele drove a lower gene expression and significantly reduced the risk of glioma. These results indicate that the underlying role of MGMT in glioma carcinogenesis needs further investigation. In addition, they warn us of the need for cancer risk stratification in study on genetic predisposition of cancer.

Acknowledgments

This study was supported by Joint Funds for the Innovation of Science and Technology, Fujian Province, China (grant no. 2016Y9016); Medical Elite Cultivation Program of Fujian Provincial Health and Family Planning Commission, China (grant no. 2016-ZQN-44); and National Natural Science Foundation, China (grant no. 81301772). The sponsors had no role in study design, data collection, data analysis, data interpretation, or writing of the report.

Supplementary material

Table S1 Information of primers for Sequenom MassARRAY iPLEX assays

Table S2 Information of probes for electrophoretic mobility shift assays

Table S3 Genotype frequencies of rs10764901 genotyped using blood DNA or adjacent normal tissue DNA in glioma patients

Table S4 Genotype frequencies of rs10764901 genotyped using blood DNA in glioma patients and controls and their contributions to the risk of glioma

Author contributions

LH and XS designed this manuscript; LH, WX, LD, DY, and SZ conducted the study and collected glioma tissues and blood samples as well as corresponding clinical characteristics of cases and controls; LH and XS analyzed the data and wrote this manuscript. All authors contributed toward data analysis, drafting, and critically revising the paper and agree to be accountable for all aspects of the work.

Disclosure

The authors report no conflicts of interest in this work.

References

  • GersonSLMGMT: its role in cancer aetiology and cancer therapeuticsNat Rev Cancer20044429630715057289
  • WibleyJEPeggAEMoodyPCCrystal structure of the human O(6)-alkylguanine-DNA alkyltransferaseNucleic Acids Res200028239340110606635
  • LiuLAllayEDumencoLLGersonSLRapid repair of O6-methylgua-nine-DNA adducts protects transgenic mice from N-methylnitrosourea-induced thymic lymphomasCancer Res19945417464846528062258
  • ZhouZQManguinoDKewittKSpontaneous hepatocellular carcinoma is reduced in transgenic mice overexpressing human O6-methylguanine-DNA methyltransferaseProc Natl Acad Sci U S A20019822125661257111606727
  • KawateHSakumiKTsuzukiTSeparation of killing and tumori-genic effects of an alkylating agent in mice defective in two of the DNA repair genesProc Natl Acad Sci U S A1998959511651209560238
  • NakamuraMWatanabeTYonekawaYKleihuesPOhgakiHPromoter methylation of the DNA repair gene MGMT in astrocytomas is frequently associated with G:C --> A:T mutations of the TP53 tumor suppressor geneCarcinogenesis200122101715171911577014
  • EstellerMToyotaMSanchez-CespedesMInactivation of the DNA repair gene O6-methylguanine-DNA methyltransferase by promoter hypermethylation is associated with G to A mutations in K-ras in colorectal tumorigenesisCancer Res20006092368237110811111
  • BrentTPHoughtonPJHoughtonJAO6-Alkylguanine-DNA alkyltransferase activity correlates with the therapeutic response of human rhabdomyosarcoma xenografts to 1-(2-chloroethyl)-3-(trans-4-methylcyclohexyl)-1-nitrosoureaProc Natl Acad Sci U S A1985829298529893857628
  • ScholdSCBrentTPvon HofeEO6-alkylguanine-DNA alk-yltransferase and sensitivity to procarbazine in human brain-tumor xenograftsJ Neurosurg19897045735772926498
  • EstellerMHamiltonSRBurgerPCBaylinSBHermanJGInactivation of the DNA repair gene O6-methylguanine-DNA methyltransferase by promoter hypermethylation is a common event in primary human neoplasiaCancer Res199959479379710029064
  • EstellerMGarcia-FoncillasJAndionEInactivation of the DNA-repair gene MGMT and the clinical response of gliomas to alkylating agentsN Engl J Med2000343191350135411070098
  • HegiMEDiserensACGorliaTMGMT gene silencing and benefit from temozolomide in glioblastomaN Engl J Med200535210997100315758010
  • UnoMOba-ShinjoSMCamargoAACorrelation of MGMT promoter methylation status with gene and protein expression levels in glioblastomaClinics201166101747175522012047
  • HsuCYLinSCHoHLExclusion of histiocytes/endothelial cells and using endothelial cells as internal reference are crucial for interpretation of MGMT immunohistochemistry in glioblastomaAm J Surg Pathol201337226427123282970
  • McDonaldKLRapkinsRWOlivierJThe T genotype of the MGMT C>T (rs16906252) enhancer single-nucleotide polymorphism (SNP) is associated with promoter methylation and longer survival in glioblastoma patientsEur J Cancer201349236036822975219
  • HsuCYHoHLLinSCHoTDHoDMThe MGMT promoter single-nucleotide polymorphism rs1625649 had prognostic impact on patients with MGMT methylated glioblastomaPLoS One20171210e018643029036186
  • Kuroiwa-TrzmielinaJWangFRapkinsRWSNP rs16906252C>T Is an Expression and Methylation Quantitative Trait Locus Associated with an Increased Risk of Developing MGMT-Methylated Colorectal CancerClin Cancer Res201622246266627727267851
  • KhuranaEFuYChakravartyDDemichelisFRubinMAGersteinMRole of non-coding sequence variants in cancerNat Rev Genet20161729310826781813
  • OldridgeDAWoodACWeichert-LeaheyNGenetic predisposition to neuroblastoma mediated by a LMO1 super-enhancer polymorphismNature2015528758241842126560027
  • PomerantzMMAhmadiyehNJiaLThe 8q24 cancer risk variant rs6983267 shows long-range interaction with MYC in colorectal cancerNat Genet200941888288419561607
  • FlicekPAmodeMRBarrellDEnsembl 2014Nucleic Acids Res201442Database issueD749D75524316576
  • ViselAMinovitskySDubchakIPennacchioLAVISTA Enhancer Browser – a database of tissue-specific human enhancersNucleic Acids Res200735Database issueD88D9217130149
  • PatenBHerreroJBealKFitzgeraldSBirneyEEnredoBEEnredo and Pecan: genome-wide mammalian consistency-based multiple alignment with paralogsGenome Res200818111814182818849524
  • 1000 Genomes Project ConsortiumAutonABrooksLDDurbinRMA global reference for human genetic variationNature20155267571687426432245
  • HuangLDaiLXuWZhangSYanDShiXIdentification of expression quantitative trait loci of MTOR associated with the progression of gliomaOncol Lett201815166567129387238
  • HuangLXuWYanDDaiLShiXIdentification of expression quantitative trait loci of RPTOR for susceptibility to gliomaTumour Biol20163722305231126361958
  • WacholderSChanockSGarcia-ClosasMEl GhormliLRothmanNAssessing the probability that a positive report is false: an approach for molecular epidemiology studiesJ Natl Cancer Inst200496643444215026468
  • ChenRSmith-CohnMCohenALColmanHGlioma Subclassifications and Their Clinical SignificanceNeurotherapeutics201714228429728281173
  • LiuYSheteSHoskingFRobertsonLHoulstonRBondyMGenetic advances in glioma: susceptibility genes and networksCurr Opin Genet Dev201020323924420211558
  • ZawlikIVaccarellaSKitaDMittelbronnMFranceschiSOhgakiHPromoter methylation and polymorphisms of the MGMT gene in glioblastomas: a population-based studyNeuroepidemiology2009321212918997474
  • FranceschiETosoniAMinichilloSThe Prognostic Roles of Gender and O6-Methylguanine-DNA Methyltransferase Methylation Status in Glioblastoma Patients: The Female PowerWorld Neurosurg2018112e342e34729337169