KRAS gene polymorphisms are associated with the risk of glioma: a two-center case-control study
Introduction
Glioma is the most common primary brain tumor in adults, accounting for more than two-thirds of all cases of primary malignant brain tumors (1). Approximate 200,000 patients worldwide are diagnosed with glioma annually, nearly 50% of the cases occurring between the ages of 40 and 60 (2-5). Studies have suggested that glioma originates from glial stem cells and progenitor cells (6). Gliomas can be classified into astrocytoma, oligodendroglioma, oligoastrocytoma, and ependymoma according to the origin of tissue cells (7,8). Mortality for glioma patients remains high despite advances in the treatment of gliomas such as radiation, chemotherapy, and surgery (9). The survival rates are about 14 months for high-grade glioma (III, IV) (10) and about 5 years for low-grade glioma (I, II) (11,12).
The formation and development of glioma is a complex process caused by abnormal genetic changes or epigenetic abnormalities (13). The etiology of glioma has been partially elucidated but remains unclear (14-17). With the progress of genomic association research and other technologies, genetic mutation has attracted more and more attention as the cause of glioma formation. Many genes have been identified to play a key role in the carcinogenesis of gliomas, such as MDM2, CDKN2A/ARF, RB, PTEN, MET, and EGFR (11,18). A group of miRNAs involved in glioma initiation by targeting RAS gene were also identified (19,20).
RAS gene which consists of KRAS, NRAS, and HRAS, is one of the most frequently activated oncogenes in human cancers. KRAS is the predominant mutated RAS gene in cancers (84% of all RAS missense mutations), followed by NRAS (12%), with HRAS rarely mutated (4%) (21). Ras has the activity of GTP enzyme, which can hydrolyze GTP to GTP (22,23). KRAS (Kirsten rat sarcoma 2 viral oncogene homolog) gene is a proto-oncogene that encodes a small GTPase transductor protein called KRAS. KRAS belongs to a group of small guanosine triphosphate (GTP) binding proteins, known as RAS superfamily or RAS-like GTPases (24). KRAS is of vital importance in the RAS gene because its mutation frequency is much higher than HRAS and NRAS (25). NRAS encodes the widely expressed small GTPase NRAS, which acts as a membrane-associated molecular switch in the RAS-MAPK pathway (26). The mutation frequency of NRAS ranks the second in the RAS gene. Q61K mutation is the most important one that is often associated with melanoma, thyroid cancer, and ovarian tumor (27). HRAS mutations, however, are the least, accounting for less than 1% of all RAS mutations. HRAS was the only oncoprotein that was significantly expressed in medullary thyroid tumor, head and neck cancer, and bladder cancer (28,29). More and more studies have shown that the RAS gene is associated with a variety of cancers. Pediatric low-grade gliomas are frequently driven by genetic alterations in the RAS/MAPK pathway (30). The genomic landscape of high-grade gliomas described from 118 pediatric patients revealed activation of the RTK/RAS/PI3K pathway in 67% of non-brain stem high-grade gliomas and 69% of diffuse intrinsic pontine gliomas (31). However, the relationship between RAS gene polymorphism and glioma has not been reported. Characterizing the role of RAS gene polymorphisms on the risk of glioma not only helps to demonstrate the specific role of RAS gene in specific cancer, but also sheds some light to the etiology of glioma. Therefore, in this study, we analyzed the relationship between NRAS and KRAS gene polymorphisms and glioma susceptibility.
We present the following article in accordance with the STROBE reporting checklist (available at http://dx.doi.org/10.21037/tp-20-359).
Methods
Study population
The current study was a two-center case-control study. The study population consisted of 191 glioma patients and 248 cancer-free controls. The subjects were all Han Chinese from southern China and were not of any blood relationship. Cases were defined by a confirmed diagnosis of glioma by histopathology. Controls were collected in the same geographical region (Guangzhou and Wenzhou) as the cases during the same period. All procedures performed in this study involving human participants were by the Declaration of Helsinki (as revised in 2013). This study was approved by the Institutional Review Committee of Guangzhou Women and Children Medical Center (No. 2016021650). Written consent was obtained from the participants’ parents or their legal guardian.
Genotyping
We selected four potential functional polymorphisms in KRAS and NRAS genes from the NCBI dbSNP database (http://www.ncbi.nlm.nih.gov/projects/SNP)and SNPinfo (http://snpinfo.niehs.nih.gov/snpfunc.htm). The genomic DNA kit (TianGen Biotech Co. Ltd., Beijing, China) was used for DNA extraction. Genomic DNA was extracted from the peripheral venous blood by subjects. Concentration of DNA was determined using the NanoDrop ND-1000 Spectrophotometer. Genotyping of NRAS and KRAS polymorphisms was performed using the TaqMan real-time PCR in a blind method by a 7900 Sequence Detection system (Thermo Fisher Scientific, Waltham, MA, USA) (32-34). The conditions of reactions were set as follow: pre-read stage at 60 °C for 30 seconds, holding stage at 95 °C 10 minutes, repeated 45 cycles each of denaturation at 95 °C for 15 seconds, annealing and extension at 60 °C for 1 minute. Genotyping was also performed on 10% of randomly selected samples to ensure genotyping accuracy. The results were 100% consistent with the original genotype.
Statistical analysis
A two-sided χ2 test was used to examine the demographic characteristics and genotype frequency distribution in patients and controls. The goodness-of-fit χ2 test was performed to evaluate the deviations from Hardy-Weinberg equilibrium (HWE) in selected polymorphisms in controls. To estimate the association between glioma risk and RAS polymorphisms, we calculated the odds ratios (ORs) and 95% confidence intervals (CIs) from the logistic regression analysis. SAS 9.1 (SAS Institute, Cary, NC, USA) was used to complete statistical analyses. The difference was considered statistically significant when P<0.05.
Results
Population characteristics
No significant differences were observed between the patients and the controls in terms of age and gender frequency (Table S1). The subtypes of glioma were astrocytic tumors (136; 71.2%), ependymoma (33; 17.28%), neuronal and mixed neuronal-glial tumors (14; 7.33%), and embryonal tumors (7; 3.66%). In these subjects, 1 case could not be recognized. According to the WHO criteria, 110 (57.59%), 38 (19.90%), 17 (8.90%), and 25 (13.09%) cases were classified into stage I, II, III, and IV, respectively, with an exception of 1 case (0.52%) which is classified into NA (not available) due to the lack of information.
Association between glioma risk and NRAS and KRAS polymorphisms
We successfully genotyped the NRAS gene SNP (rs2273267 A>T) and KRAS gene SNPs (rs12587 G>T, rs7973450 A>G, rs7312175 G>A) in 248 controls and 191 cases. The correlations of these SNPs with glioma risk are shown in Table 1. All the SNP genotype frequencies in controls follow HWE. The result indicated that the KRAS rs7312175 G>A (GA versus GG: adjusted OR =1.66, 95% CI: 1.05–2.64, P=0.030) variant alleles were associated with increased glioma susceptibility. The SNPs, NRAS rs2273267 A>T, KRAS rs12587 G>T, and KRAS rs7973450 A>G, could not impact the risk of glioma, respectively.
Full table
Stratification analysis
Stratification analysis was conducted to explore the association between the KRAS rs7212175 G>A polymorphism and glioma susceptibility by age, gender, subtypes, and clinical stages (Table 2). Compared with the GG genotype, the rs7212175 GA/AA genotypes were associated with an increased risk of glioma at age less than 60 months (adjusted OR =2.47, 95% CI: 1.30–4.69, P=0.006), males (adjusted OR =1.96, 95% CI: 1.06–3.64, P=0.033), astrocytic tumors (adjusted OR =1.84, 95% CI: 1.13–3.00, P=0.015), patients in clinical stages I (adjusted OR =1.83, 95% CI: 1.09–3.06, P=0.022), and patients in clinical stage I+II (adjusted OR =1.63, 95% CI: 1.01–2.63, P=0.045). The presence of two to three risk genotypes was associated with an enhanced glioma risk in age <60 months (adjusted OR =2.09, 95% CI: 1.21–3.62, P=0.009), astrocytic tumors (adjusted OR =1.62, 95% CI: 1.05–2.50, P=0.029), and patients in clinical stages I (adjusted OR =1.69, 95% CI: 1.07–2.68, P=0.025). We also examined the relationship between NRAS rs2273267 A>T polymorphism and glioma risk under the same subgroups yet failed to obtain any significant results (Table 3).
Full table
Full table
Discussion
Glioma is an extremely malignant tumor that accounts for the majority of intracranial tumors and has a high recurrence rate and mortality rate (35). To better understand the pathogenesis of glioma, improve the survival rate, and reduce the mortality rate in patients, we conducted this study to determine the relationship between KRAS and NRAS gene SNPs and glioma susceptibility. The results showed that KRAS gene rs7312175 G>A was significantly correlated with glioma susceptibility.
KRAS, located on the short arm of chromosome 12, encodes active guanosine diphosphate (GDP)/inactive guanosine triphosphate (GTP) binding protein, a membrane-bound protein involved in cell signal transduction (36). RAS protein plays a critical role in cell switch between the GTP binding state and GDP binding state circulation. Under normal circumstances RAS and GDP binding pathway is in a static state, whereas in the condition of external stimulation RAS and GTP binding signal is transmitted to the downstream pathway to regulate cell proliferation and differentiation (21,37,38). When KRAS mutation abnormally activates this protein, its GTPase activity is weakened, leading to increased binding of GTP to activated KRAS, which activates the effector signaling pathway and drives abnormal cell proliferation and tumor formation (1,39). More and more studies have been conducted to explore the role of KRAS gene SNPs on cancer. Chin et al. explored the effects of let-7 binding single nucleotide polymorphisms (SNPs) on target gene regulation and the risk of non-small cell carcinoma in a case-control study. The results identified a mutant gene at the let-7 complementary locus in KRAS 3’-UTR, which altered let-7 mediated KRAS expression regulation and significantly increased the risk of non-small cell lung cancer in moderate smokers (40). Liu et al. confirmed that KRAS gene SNP rs8720 CT increased mRNA and protein expression of KRAS and thus contribute to increased colorectal carcinoma (CRC) risk and reduced survival (41). Smits et al. showed that KRAS-LCS6 mutation existed in early CRC, which could indicate a good prognosis of patients (42).
Our previous study investigated the relationship between KRAS/NRAS gene polymorphisms and the risk of pediatric tumors. We found that KRAS/NRAS polymorphisms are irrelevant to hepatoblastoma susceptibility in Chinese population (23), while KRAS rs12587 G>T was significantly associated with the susceptibility to nephroblastoma (29). In this study, we are the first to investigate the relationship between KRAS/NRAS gene SNPs and glioma susceptibility. In the four SNPs analyzed, KRAS rs7312175 G>A is associated with the risk of glioma, while NRAS rs2273267 A>T, KRAS rs12587 G>T, KRAS rs7973450 A>G are not associated with glioma risk. Stratified analysis showed that individuals harboring the rs7312175 GA alleles were more likely to develop glioma in the following subgroups: children less than 60 months, male, tumorous subtypes were astrocytic tumors and clinical stages I.
Although this study verified the correlation between the KRAS gene SNPs and glioma risk to some extent and obtained positive results, there were still some deficiencies. Firstly, the sample size was small. Secondly, the subjects of this study were all Han people, which could not be applied to other non-Han people. Besides, only 4 gene polymorphisms were considered, and environmental exposure factors were not considered. Of note, in the current study, the HRAS gene SNPs were not included. More SNPs are warranted to be analyzed in the future.
In conclusion, we demonstrated that KRAS rs7312175 G>A polymorphism was significantly associated with the risk of glioma. This study has a certain significance for the better understanding of glioma. The SNP identified here may be used for the diagnosis and treatment of the disease to improve the survival rate of glioma patients.
Acknowledgments
Funding: This study was supported by grants from the National Natural Science Foundation of China (No. 81802346, 82002636) and China Postdoctoral Science Foundation (No. 2020T130132).
Footnote
Reporting Checklist: The authors have completed the STROBE reporting checklist. Available at http://dx.doi.org/10.21037/tp-20-359
Data Sharing Statement: Available at http://dx.doi.org/10.21037/tp-20-359
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at http://dx.doi.org/10.21037/tp-20-359). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. All procedures performed in this study involving human participants were by the Declaration of Helsinki (as revised in 2013). Our study design received approval from institutional review board of Guangzhou Women and Children’s Medical Center (No. 2016021650). Written informed consent was acquired from the parents or guardians of each participant.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Janes MR, Zhang J, Li LS, et al. Targeting KRAS Mutant Cancers with a Covalent G12C-Specific Inhibitor. Cell 2018;172:578-89.e17. [Crossref] [PubMed]
- Ostrom QT, Bauchet L, Davis FG, et al. The epidemiology of glioma in adults: a "state of the science" review. Neuro Oncol 2014;16:896-913. [Crossref] [PubMed]
- Bilmin K, Kujawska T, Grieb P. Sonodynamic Therapy for Gliomas. Perspectives and Prospects of Selective Sonosensitization of Glioma Cells. Cells 2019;8:1428. [Crossref] [PubMed]
- Kros JM, Mustafa DM, Dekker LJ, et al. Circulating glioma biomarkers. Neuro Oncol 2015;17:343-60. [PubMed]
- Chen WL, Chen HJ, Hou GQ, et al. LINC01198 promotes proliferation and temozolomide resistance in a NEDD4-1-dependent manner, repressing PTEN expression in glioma. Aging (Albany NY) 2019;11:6053-68. [Crossref] [PubMed]
- Weller M, Wick W, Aldape K, et al. Glioma. Nat Rev Dis Primers 2015;1:15017. [Crossref] [PubMed]
- Chen R, Smith-Cohn M, Cohen AL, et al. Glioma Subclassifications and Their Clinical Significance. Neurotherapeutics 2017;14:284-97. [Crossref] [PubMed]
- Xi J, Sun Q, Ma L, et al. Long non-coding RNAs in glioma progression. Cancer Lett 2018;419:203-9. [Crossref] [PubMed]
- Xu P, Tao X, Zhao C, et al. DTX3L is upregulated in glioma and is associated with glioma progression. Int J Mol Med 2017;40:491-8. [Crossref] [PubMed]
- Stylli SS, Luwor RB, Ware TM, et al. Mouse models of glioma. J Clin Neurosci 2015;22:619-26. [Crossref] [PubMed]
- Yamada D, Fujikawa K, Kawabe K, et al. RUNX2 Promotes Malignant Progression in Glioma. Neurochem Res 2018;43:2047-54. [Crossref] [PubMed]
- Wang S, Yao F, Lu X, et al. Temozolomide promotes immune escape of GBM cells via upregulating PD-L1. Am J Cancer Res 2019;9:1161-71. [PubMed]
- Tsankova NM, Canoll P. Advances in genetic and epigenetic analyses of gliomas: a neuropathological perspective. J Neurooncol 2014;119:481-90. [Crossref] [PubMed]
- Meng Q, Li S, Liu Y, et al. Circular RNA circSCAF11 Accelerates the Glioma Tumorigenesis through the miR-421/SP1/VEGFA Axis. Mol Ther Nucleic Acids 2019;17:669-77. [Crossref] [PubMed]
- Moradimotlagh A, Arefian E, Rezazadeh Valojerdi R, et al. MicroRNA-129 Inhibits Glioma Cell Growth by Targeting CDK4, CDK6, and MDM2. Mol Ther Nucleic Acids 2020;19:759-64. [Crossref] [PubMed]
- Cordier D, Gerber M, Brand S. Effects of two types of exercise training on psychological well-being, sleep, quality of life and physical fitness in patients with high-grade glioma (WHO III and IV): study protocol for a randomized controlled trial. Cancer Commun (Lond) 2019;39:46. [Crossref] [PubMed]
- Wang FY, Wang-Gou SY, Cao H, et al. Proteomics identifies EGF-like domain multiple 7 as a potential therapeutic target for epidermal growth factor receptor-positive glioma. Cancer Commun (Lond) 2020;40:518-30. [Crossref] [PubMed]
- Braganhol E, Wink MR, Lenz G, et al. Purinergic Signaling in Glioma Progression. Adv Exp Med Biol 2020;1202:87-108. [Crossref] [PubMed]
- Zhi T, Jiang K, Zhang C, et al. MicroRNA-1301 inhibits proliferation of human glioma cells by directly targeting N-Ras. Am J Cancer Res 2017;7:982-98. [PubMed]
- Li Y, Li Y, Ge P, et al. MiR-126 Regulates the ERK Pathway via Targeting KRAS to Inhibit the Glioma Cell Proliferation and Invasion. Mol Neurobiol 2017;54:137-45. [Crossref] [PubMed]
- Waters AM, Der CJ. KRAS: The Critical Driver and Therapeutic Target for Pancreatic Cancer. Cold Spring Harb Perspect Med 2018;8:a031435 [Crossref] [PubMed]
- Hallin J, Engstrom LD, Hargis L, et al. The KRAS(G12C) Inhibitor MRTX849 Provides Insight toward Therapeutic Susceptibility of KRAS-Mutant Cancers in Mouse Models and Patients. Cancer Discov 2020;10:54-71. [Crossref] [PubMed]
- Yang T, Wen Y, Li J, et al. NRAS and KRAS polymorphisms are not associated with hepatoblastoma susceptibility in Chinese children. Exp Hematol Oncol 2019;8:11. [Crossref] [PubMed]
- Simanshu DK, Nissley DV, McCormick F. RAS Proteins and Their Regulators in Human Disease. Cell 2017;170:17-33. [Crossref] [PubMed]
- Haigis KM. KRAS Alleles: The Devil Is in the Detail. Trends Cancer 2017;3:686-97. [Crossref] [PubMed]
- Cirstea IC, Kutsche K, Dvorsky R, et al. A restricted spectrum of NRAS mutations causes Noonan syndrome. Nat Genet 2010;42:27-9. [Crossref] [PubMed]
- Song Z, Liu F, Zhang J. Targeting NRAS(Q61K) mutant delays tumor growth and angiogenesis in non-small cell lung cancer. Am J Cancer Res 2017;7:831-44. [PubMed]
- Untch BR, Dos Anjos V, Garcia-Rendueles MER, et al. Tipifarnib Inhibits HRAS-Driven Dedifferentiated Thyroid Cancers. Cancer Res 2018;78:4642-57. [Crossref] [PubMed]
- Fu W, Zhuo Z, Hua RX, et al. Association of KRAS and NRAS gene polymorphisms with Wilms tumor risk: a four-center case-control study. Aging (Albany NY) 2019;11:1551-63. [Crossref] [PubMed]
- Ryall S, Zapotocky M, Fukuoka K, et al. Integrated Molecular and Clinical Analysis of 1,000 Pediatric Low-Grade Gliomas. Cancer Cell 2020;37:569-83.e5. [Crossref] [PubMed]
- Wu G, Diaz AK, Paugh BS, et al. The genomic landscape of diffuse intrinsic pontine glioma and pediatric non-brainstem high-grade glioma. Nat Genet 2014;46:444-50. [Crossref] [PubMed]
- Li Y, Zhuo ZJ, Zhou H, et al. H19 gene polymorphisms and neuroblastoma susceptibility in Chinese children: a six-center case-control study. J Cancer 2019;10:6358-63. [Crossref] [PubMed]
- Zhuo Z, Lu H, Zhu J, et al. METTL14 Gene Polymorphisms Confer Neuroblastoma Susceptibility: An Eight-Center Case-Control Study. Mol Ther Nucleic Acids 2020;22:17-26. [Crossref] [PubMed]
- Zhuo Z, Zhou C, Fang Y, et al. Correlation between the genetic variants of base excision repair (BER) pathway genes and neuroblastoma susceptibility in eastern Chinese children. Cancer Commun (Lond) 2020;40:641-6. [Crossref] [PubMed]
- Terrasi A, Bertolini I, Martelli C, et al. Specific V-ATPase expression sub-classifies IDHwt lower-grade gliomas and impacts glioma growth in vivo. EBioMedicine 2019;41:214-24. [Crossref] [PubMed]
- Ying HQ, Wang F, He BS, et al. The involvement of Kras gene 3'-UTR polymorphisms in risk of cancer and influence on patient response to anti-EGFR therapy in metastatic colorectal cancer: a meta-analysis. Onco Targets Ther 2014;7:1487-96. [PubMed]
- Herbst RS, Schlessinger J. Small molecule combats cancer-causing KRAS protein at last. Nature 2019;575:294-5. [Crossref] [PubMed]
- Yin C, Zhu B, Zhang T, et al. Pharmacological Targeting of STK19 Inhibits Oncogenic NRAS-Driven Melanomagenesis. Cell 2019;176:1113-27.e16. [Crossref] [PubMed]
- Ferrer I, Zugazagoitia J, Herbertz S, et al. KRAS-Mutant non-small cell lung cancer: From biology to therapy. Lung Cancer 2018;124:53-64. [Crossref] [PubMed]
- Chin LJ, Ratner E, Leng S, et al. A SNP in a let-7 microRNA complementary site in the KRAS 3' untranslated region increases non-small cell lung cancer risk. Cancer Res 2008;68:8535-40. [Crossref] [PubMed]
- Liu X, Yang G, Huang J, et al. KRAS SNPs are related to colorectal cancer susceptibility and survival in Chinese people. Biomark Med 2020;14:13-22. [Crossref] [PubMed]
- Smits KM, Paranjape T, Nallur S, et al. A let-7 microRNA SNP in the KRAS 3'UTR is prognostic in early-stage colorectal cancer. Clin Cancer Res 2011;17:7723-31. [Crossref] [PubMed]