Abstract
Recent advances in human genetic studies have opened new avenues for the identification of susceptibility genes for many complex genetic disorders, especially in the field of rare cancers such as glioma. Glioma is one of the least understood human tumors and the etiology for glioma is barely known. Hundreds of single-nucleotide polymorphisms (SNPs) are found to be related to the risk of glioma in previous studies. This study is committed to investigate the role of heredity in this disorder. To examine and validate how common variants contribute to glioma susceptibility in the Han Chinese population, we evaluated 12 tagging SNPs in a case–control study in the Chinese Han population from Xi'an city of China (301 cases and 302 controls). Overall, two protective alleles and one risk allele for glioma were found by genetic model analyses. In dominant model, the allele “T” of rs6947203 in the RPA3 gene acts as a protective allele [odds ratio (OR), 0.59; 95% confidence interval (CI), 0.22–0.90; p=0.014]. In recessive model, the allele “C” of rs1042522 in the TP53 gene acts as a risk allele (OR, 1.65; 95% CI, 1.05–2.59; p=0.0314). In additive model, the allele “G” of rs4140805 in the RPA3 gene (OR, 0.73; 95% CI, 0.53–0.99; p=0.0437) and the allele “T” of rs6947203 in the RPA3 gene (OR, 0.62; 95% CI, 0.42–0.92; p=0.0177) both act as protective alleles. We also observed a haplotype of “CC” in the TP53 gene with an increased risk of 34% of developing glioma (p=0.0306). Our results, combined with previous studies, ascertain the potential role of the TP53 gene to glioma onset.
Introduction
Glioma is the most common brain tumor in humans and has a very poor prognosis with only 3% of glioblastoma patients surviving 5 years after diagnosis. 1 Over the past years, several relevant biomarkers with a potential clinical interest have been identified in gliomas using various techniques, such as karyotype analysis, microsatellite analysis, fluorescent in situ hybridization, and array-based chromosome comparative genomic hybridization technology. 2 Current evidence suggests that the heredity risk is due to the coinheritance of multiple low-risk variants. 3 Genomic alterations described in gliomas, such as EGFR, CDKN2A, TP53, RB1, and PTEN, have been found in other common epithelial tumors. 4,5 Several genes affecting cell cycle and DNA repair functions have been proposed to play a role in glioma pathogenesis and progression. 6,7 Recent studies have emphasized on identifying the relationship between glioma susceptibility and single-nucleotide polymorphism (SNP). 8
In this study, we investigate and validate the potential associations between glioma risk and 12 tagging SNPs (tSNPs) that were previously reported to be associated with glioma susceptibility. Our data shed new light on the association between common SNPs and glioma susceptibility in the Chinese population.
Materials and Methods
Study participants
We recruited a total of recently diagnosed and histologically confirmed 350 glioma patients from September 2009 to November 2011 into an ongoing molecular epidemiological study at the Department of Neurosurgery of the Tangdu Hospital affiliated with the Fourth Military Medical University in Xi'an city, China. All glioma cases were previously healthy. Age, gender, or disease stage restrictions were ignored for case recruitment. In our study, we excluded 49 cases that had incomplete clinical information. Finally, 301 glioma cases were successfully genotyped.
A random sample of unrelated healthy individuals from July 2011 to October 2011 were recruited from the medical center at Tangdu Hospital for genetic association research of human complex diseases, such as liver, renal cancer, and glioma. The recruitment and exclusion standard were used. Generally, all the subjects were healthy and had no chronic diseases involving vital organs (lung, heart, kidney, brain, and liver), and thus we reduced the known environmental and therapeutic factors that influence the variation of human complex diseases. We recruited a total of 302 healthy controls in this study. All the participants were restricted to Han Chinese who lived in Xi'an city and its surrounding areas.
Clinical data and demographic
A standardized questionnaire including information such as region, age, gender, alcohol use, ethnicity, education status, and family history of cancer were used to collect necessary data through in-person interview. For patients, we consulted treating physicians or medical chart review to collect related information. We also tested the α-fetoprotein and plasma carcinoembryonic antigen to ensure the quality of controls.
SNP selection and genotyping
Twelve tSNPs with minor allele frequency (MAF) >5% in HapMap Asian population in nine genes (CCNH, CHAF1A, ERCC1, PARP1, RPA3, TP53, XRCC3, XRCC4, and XRCC6) previously reported to be associated with 3,9 –13 glioma risk were successfully genotyped. Genomic DNA from whole blood was extracted using phenol–chloroform method and their concentration was measured by spectrometry (DU530 UV/VIS spectrophotometer; Beckman Instruments, Fullerton, CA). We used Sequenom MassARRAY Assay Design 3.0 Software to design Multiplexed SNP MassEXTEND assay. 14 We performed Sequenom MassARRAY RS1000 to genotype SNPs using the standard protocol recommended by the manufacturer. Sequenom Typer 4.0 Software was used to perform data management and analyses. 14,15
Statistical analysis
Microsoft Excel and SPSS 16.0 statistical package (SPSS, Chicago, IL) were used for data management. All p-values in this study were two-sided and p=0.05 was considered the cutoff value of statistical significance. An exact test was used to test each tSNP frequency's departure from Hardy–Weinberg equilibrium in control subjects. χ2 test was used to calculate the genotype frequencies of cases and controls. 16 We tested odds ratios (ORs) and 95% confidence intervals (CIs) using unconditional logistic regression analysis with adjustment for age and gender. 17 The possibility of gender differences as a source of population substructure was evaluated by a genotype test for each tSNP in men and women, and the number of significant results at the 5% level was compared with the number expected by χ2 test. 16 We did not detect ethnic specificity because all participants' ethnicity was limited to Han Chinese.
We used three models (allele model, dominant model, and recessive model) to perform the associations of certain tSNP contributed to the risk of glioma. ORs and 95% CIs were calculated using unconditional logistic regression analysis with adjustment for age and gender. 17
We also used the Haploview software package (version 4.2) to estimate the pairwise linkage disequilibrium (LD) between markers and partition haplotype blocks. 18
Finally, we used the SHEsis software platform (
Results
We included a total of 603 participants, with 301 patients (157 men, 144 women; median age at diagnosis 41.5 years) and 302 controls (155 men, 145 women; median age 42.3 years) for genetic association analyses. Basic characteristics of the cases such as gender, age, and pathology are listed in Table 1. The primers of the 12 selected tSNPs are shown in Table 2, which were designed by Sequenom MassARRAY Assay Design 3.0 Software. 14 Twelve tSNPs were successfully genotyped in all cases and controls. Average call rate of tSNPs was 99.57% (range from 98.51% to 100%) (listed in Table 3). The relationships between the tSNPs and glioma risk were listed in Table 4. We found no association between tSNPs and glioma onset at 5% level.
SNP, single-nucleotide polymorphism; PCR, polymorphism chain reaction; UEP, unextended mini-sequencing primer.
HWE, Hardy–Weinberg equilibrium.
OR, odds ratio; CI, confidence interval.
We further analyzed the association of tSNPs and glioma risk using logistic test including dominant model, recessive model, and additive model. The MA of each tSNP was assumed to be a risk factor and their frequencies (MAF) were listed in Table 5. We found two protective alleles and one risk allele, the allele “T” of rs6947203 in the gene RPA3 acts as a protective allele in dominant (OR, 0.59; 95% CI, 0.22–0.90; p=0.014) and additive (OR, 0.62; 95% CI, 0.42–0.92; p=0.0177) model. The allele “G” of rs4140805 in the gene RPA3 (OR, 0.73; 95% CI, 0.53–0.99; p=0.0437) was found to play a protective role in additive model. The allele “C” of rs1042522 in the gene TP53 (OR, 1.65; 95% CI, 1.05–2.59; p=0.0314) acts as a risk allele in recessive model.
means the p-values ≤0.05 and have statistical significance.
MAF, minor allele frequency.
The relationship between the TP53 haplotypes and glioma risk are listed in Table 6. Haplotype “CC” in the TP53 gene was found to be associated with risk of glioma (OR, 1.3384; 95% CI, 1.0273–1.7437; Fisher's p=0.0307; Pearson's p=0.0306).
p value≤0.05 indicates statistical significance.
Discussion
In the current study, we successfully genotyped 12 tSNPs in nine genes (CCNH, CHAF1A, ERCC1, PARP1, RPA3, TP53, XRCC3, XRCC4, and XRCC6) in the Han Chinese population and indentified 2 protective tSNPs in the RPA3 gene (allele “T” of rs6947203 and allele “G” of rs4140805) and 1 risk tSNP in the TP53 gene (allele “C” of rs1042522) to be associated with glioma. We also observed a haplotype of “CC” of the TP53 gene with an increased risk of 34% of developing glioma.
TP53 is located in 17p13.1, encodes tumor protein TP53, which responds to diverse cellular stresses to regulate target genes that induce cell cycle arrest, apoptosis, senescence, DNA repair, or changes in metabolism. TP53 protein is expressed at low level in normal cells and at a high level in a variety of transformed cell lines, where it is believed to contribute to transformation and malignancy. TP53 is a DNA-binding protein containing transcription activation, DNA-binding, and oligomerization domains. It is postulated to bind to a TP53-binding site and activate expression of downstream genes that inhibit growth and/or invasion, and thus function as a tumor suppressor. Mutants of TP53 that frequently occur in a number of different human cancers fail to bind the consensus DNA binding site, and hence cause the loss of tumor suppressor activity. Alterations of this gene occur not only as somatic mutations in human malignancies, but also as germline mutations in some cancer-prone families with Li–Fraumeni syndrome. Multiple TP53 variants due to alternative promoters and multiple alternative splicing have been found. These variants encode distinct isoforms, which can regulate TP53 transcriptional activity. Our result demonstrates that mutation in TP53 plays a risk role in glioma onset; this is supported by the role for TP53 inactivation in astrocytoma initiation and is consistent with the high frequency of TP53 mutations in Grade II human astrocytomas. 20,21 Another study in Indians showed higher Arg/Arg (rs1042522 G/C) genotype in gliomas compared with normal population (38% vs. 13%). 22 Previous study in New York also found that for Ex4+119 C>G SNP (rs1042522), women with the heterozygous genotype (G/C) had a 32% increase in breast cancer risk. 23 Ours and the previous study analysis of SNP suggest that TP53 genotypes, haplotypes, and locus–locus interactions are associated with glioma risk. These findings indicate TP53 gene plays an important role in glioma onset.
As to gene RPA3, which is located in 7p22, two protective tSNPs (allele “T” of rs6947203 and allele “G” of rs4140805) have been found for this gene in our study; however, we have not found any previous study that showed the relationship between these two tSNPs and glioma onset, in fact, study based on this gene is rare. The protein encoded by this gene has 8 exons and 121 amino acids, it is a single-stranded DNA-binding protein that functions in many aspects of DNA metabolism and has a central role in DNA replication, playing an essential function in both initiation and elongation. 24 –28
Our study is the first to find SNPs in the RPA3 gene are associated to glioma onset. Further study should be concentrated on the exact mechanism of how tSNPs affect glioma onset and whether it has population differentiation.
Some limitations cannot be ignored in our study. First, population admixture, which may cause type-I error (false positive) for association study. In our study, all the samples we used were from the same hospital to avoid two or more definite selection bias and they are neutralized because they did not differ in geographical distributions or genotype frequencies. The race of all participants was limited to Han Chinese who lived in Xi'an city or nearby, and so substantial population admixture can be ignored in our study.
Second, sample size (301 glioma patients and 302 control subjects) was not relatively large enough for association studies. We could not find any gender-specific significant tSNPs because our sample size was not large enough and could not be subgrouped by gender (data not shown). We selected SNPs whose frequencies were higher than 5% in HapMap Asian populations to affirm the statistical power. We also planned a haplotype-based study to investigate whether certain haplotypes were associated with glioma onset.
In conclusion, our study provides new evidences for the relationship between tSNPs and glioma onset, which may shed light on the etiology of glioma and the tSNPs we found can be applied in clinical diagnosis and therapy.
Footnotes
Acknowledgments
The authors are grateful to all the patients and individuals for their participation. We would also like to thank the clinicians and other hospital staff who contributed to the blood sample and data collection for this study.
Disclosure Statement
This work is supported by the National 863 High-Technology Research and Development Program (No. 2012AA02A519).
