- Research article
- Open Access
Genetic variants in MIR17HG affect the susceptibility and prognosis of glioma in a Chinese Han population
BMC Cancer volume 20, Article number: 976 (2020)
lncRNA MIR17HG was upregulated in glioma, and participated in promoting proliferation, migration and invasion of glioma. However, the role of MIR17HG polymorphisms in the occurrence and prognosis of glioma is still unclear.
In the study, 592 glioma patients and 502 control subjects were recruited. Agena MassARRAY platform was used to detect the genotype of MIR17HG polymorphisms. Logistic regression analysis was used to evaluate the relationship between MIR17HG single nucleotide polymorphisms (SNPs) and glioma risk by odds ratio (OR) and 95% confidence intervals (CIs). Kaplan–Meier curves, Cox hazards models were performed for assessing the role of these SNPs in glioma prognosis by hazard ratios (HR) and 95% CIs.
We found that rs7318578 (OR = 2.25, p = 3.18 × 10− 5) was significantly associated with glioma susceptibility in the overall participants. In the subgroup with age < 40 years, rs17735387 (OR = 1.53, p = 9.05 × 10− 3) and rs7336610 (OR = 1.35, p = 0.016) were related to the higher glioma susceptibility. More importantly, rs17735387 (HR = 0.82, log-rank p = 0.026) were associated with the longer survival of glioma patients. The GA genotype of rs17735387 had a better overall survival (HR = 0.75, log-rank p = 0.013) and progression free survival (HR = 0.73, log-rank p = 0.032) in patients with I-II glioma. We also found that rs72640334 was related to the poor prognosis (HR = 1.49, Log-rank p = 0.035) in female patients. In the subgroup of patients with age ≥ 40 years, rs17735387 was associated with a better prognosis (HR = 0.036, Log-rank p = 0.002).
Our study firstly reported that MIR17HG rs7318578 was a risk factor for glioma susceptibility and rs17735387 was associated with the longer survival of glioma among Chinese Han population, which might help to enhance the understanding of MIR17HG gene in gliomagenesis. In subsequent studies, we will continue to collect samples and follow up to further validate our findings and further explore the function of these MIR17HG SNPs in glioma in a larger sample size.
Glioma is the most frequent neoplasms originated from neuroglial stem or progenitor cells, accounting for 80% of primary malignant brain cancers with approximately 101,600 individuals diagnosed in China each year [1, 2]. Despite the efforts of diagnosis and therapeutics, the prognosis of glioma is still depressing. Until now, the aetiology of glioma remains unclear. However, environmental and occupational exposures have been identified to be associated with the occurrence and development of glioma, especially high-dosage ionizing radiation . In addition, genetic factors are also given a pivotal contribution to the occurrence and prognosis of glioma [4,5,6]. Several association studies have revealed that single nucleotide polymorphisms (SNPs) were associated with glioma risk and survival [7,8,9].
MIR17HG gene, located on chromosome 13q31.3, is the host gene of the microRNA 17–92 cluster. Functional studies have confirmed that the MIR17HG gene might be related to cell survival, proliferation, differentiation, and angiogenesis . LncRNA MIR17HG, also as a long noncoding RNA which regulating the expression of miRNA, played a carcinogenic effect in various cancers including rectal cancer, gastric cancer, and lung cancer [11,12,13]. A recent research showed that lncRNA MIR17HG was overexpressed in glioma, and lncRNA MIR17HG knockdown inhibited the proliferation, migration and invasion of glioma, suggesting that lncRNA MIR17HG might facilitate the malignant progress of glioma . Recently, increasing evidences indicated that genetic polymorphisms of MIR17HG were associated with the occurrence of multiple tumors, such as lymphoma, colorectal cancer, breast cancer [15,16,17]. However, the role of MIR17HG variants in glioma occurrence and prognosis is still unclear.
Here, we analyzed the association of selected MIR17HG SNPs and glioma susceptibility among the Chinese Han population, and examined the possible role of these polymorphisms in different glioma subgroups stratified by age, gender and grade. We also evaluated the influence of MIR17HG genetic variants on the survival of glioma patients.
This study recruited 592 glioma patients and 502 control subjects. All participants were genetically unrelated Chinese Han population. Glioma patients who diagnosed and confirmed by histopathology were enrolled from the department of Neurosurgery at Tangdu Hospital from February 2014 to March 2018. Patients with history of cancer and other systemic or complex diseases were excluded. Age- and gender-matched healthy controls were recruited from the physical examination center of the hospital. The controls were free from any cancer and any disease related to brain and central nervous system. Standardized questionnaires and medical records were used to collect demographic and clinical information. The follow-up information was obtained by telephone and return visit every 3 months; and the survival time, progress and outcome were recorded. After, approximately 5 mL blood samples were collected for further analysis. Our study was approved by the Ethics Committee of the Second Affiliated Hospital of Hainan Medical University and was in the Declaration of Helsinki. Written informed consent was obtained from each participant.
Genomic DNA was purified by a commercially available GoldMag DNA Purification Kit (GoldMag Co. Ltd., Xi′an City, China). NanoDrop 2000 (Thermo Scientifc, Waltham, MA, USA) was used to check DNA quality. Five MIR17HG SNPs (rs17735387, rs72640334, rs7318578, rs7336610, and rs75267932) were identified based on the NCBI dbSNP database, the 1000 Genomes Project data with minor allele frequencies (MAFs) > 5% in Chinese Han Beijing (CHB) population and Haploview software with a pairwise linkage disequilibrium (r2 > 0.80). MIR17HG polymorphisms were genotyped using Agena MassARRAY platform (Agena, San Diego, CA, U.S.A.) as previously described . The primers sequences were presented in Supplementary Table 1. Genotyping was in a blinded manner, and the call rate was ≥0.99. For quality control, 10% of blind and random samples were repeated genotyping, and the result was 100% reproducibility.
Statistical analysis were performed using SPSS version 18.0 (SPSS Inc., Chicago, IL, USA) and PLINK 2.1.7 package. The Chi square test or Student’s t-test was carried out to compare the differences in age and gender distributions between patients and controls, as appropriate. Hardy–Weinberg equilibrium (HWE) was performed for the controls using goodness-of-fit χ2 test. Logistic regression analysis was used to analyze the genetic effects of MIR17HG SNPs on the risk of glioma by calculating odds ratio (OR) and 95% confidence intervals (CIs) adjusted for age and sex. Multiple testing correction was performed by the false discovery (FDR). The overall survival (OS) and progression-free survival (PFS) of glioma patients were plotted by Kaplan–Meier survival curves. Univariate and multivariate Cox proportional hazards models were performed to assess the role of MIR17HG polymorphisms in the prognosis of glioma by calculating hazard ratio (HR) and 95% CIs. A two-tailed p value of < 0.05 was statistically significant.
The characteristics of patients and controls were presented in Table 1. The case group consisted of 592 glioma patients (40.53 ± 13.90 years, 55.1% males) and 502 healthy controls (40.46 ± 18.08 years, 54.8% males). The frequency distribution of age (p = 0.934) and sex (p = 0.924) between cases and controls were no statistical differences. Among the cases, there were 378 patients with WHO 2007 grade I + II and 214 patients with grade III + IV.
The genotyping results of MIR17HG variants
Five SNPs in MIR17HG were genotyped to determine the possible effect of MIR17HG variants on the risk or prognosis of glioma. The minor allele frequencies in patients and controls were displayed in Supplementary Table 2. The genotype frequencies of all the studied variants in the control group were in HWE (p > 0.05), and the genotyping rate exceeded 99.5%.
The correlation between MIR17HG variants and glioma risk
The genotype and allele frequencies of these SNPs in MIR17HG were displayed in Table 2. Compared with the control group, the frequencies of C allele (34.9% vs 28.9%) and CC genotype (19.7% vs 9.0%) of rs7318578 were higher in glioma patients. In details, rs7318578 C allele (OR = 1.32, 95% CI: 1.10–1.58, p = 2.63 × 10− 3) and CC genotype (OR = 2.25, 95% CI: 1.54–3.31, p = 3.18 × 10− 5) were related to the increased glioma susceptibility compared with the A allele and AA genotype, respectively, and the significance still existed after the FDR controlling procedure (FDR-p = 0.032 and FDR-p = 0.001 respectively). Moreover, rs7318578 variant showed a 1.26-fold increased risk of glioma under the additive model (OR = 1.26, 95% CI: 1.07–1.49, p = 6.23 × 10− 3). There was no association between other SNPs and the risk of glioma.
We further explored the association between glioma risk and MIR17HG SNPs by stratifying for age, sex and WHO grade. Among subjects of age ≥ 40 years, carriers with rs7318578 CC genotype showed a 2.46-fold increased the susceptibility to glioma compared with individuals with the AA genotype (OR = 2.46, 95% CI: 1.42–4.28, p = 1.41 × 10− 3, FDR-p = 0.035, Table 3). Additionally, rs17735387 was a risk factor for glioma occurrence: A vs G: OR = 1.53, 95% CI: 1.11–2.11, p = 9.05 × 10− 3; AA vs GG: OR = 3.27, 95% CI: 1.09–9.80, p = 0.034; GA + AA vs GG: OR = 1.57, 95% CI: 1.07–2.30, p = 0.021; additive: OR = 1.56, 95% CI: 1.12–2.18, p = 8.55 × 10− 3 at age < 40 years. MIR17HG rs7318578 C allele (OR = 1.37, 95% CI: 1.05–1.79, p = 0.020) and CC genotype (OR = 1.88, 95% CI: 1.08–3.28, p = 0.026) was associated with the increased risk of glioma in subjects aged younger 40 years. Results of multiple models showed that rs7336610 was associated with the high glioma susceptibility at age < 40 years (C vs T: OR = 1.35, 95% CI: 1.06–1.73, p = 0.016; TC vs TT: OR = 1.56, 95% CI: 1.02–2.39, p = 0.041; CC vs TT: OR = 1.72, 95% CI: 1.02–2.92, p = 0.044; TC + CC vs TT: OR = 1.61, 95% CI: 1.07–2.41, p = 0.022; additive: OR = 1.33, 95% CI: 1.02–1.73, p = 0.034).
Stratified by gender (Table 3), the significant association between rs7318578 and the glioma of risk was observed in males (CC vs AA: OR = 1.80, 95% CI: 1.10–2.95, p = 0.020) and females (CC vs AA: OR = 3.08, 95% CI: 1.67–5.67, p = 3.19 × 10− 4, FDR-p = 0.046 and additive: OR = 1.43, 95% CI: 1.11–1.84, p = 5.96 × 10− 3). Especially, the association under the allele model in females was still significant (C vs A: OR = 1.53, 95% CI: 1.16–2.01, p = 2.49 × 10− 3, FDR-p = 0.029).
In the stratified analysis by WHO grade, rs7336610 showed a genotype difference between patients with grade III-IV and patients with grade I-II, with OR from 1.31 to 1.72 (TC vs TT: OR = 1.58, 95% CI: 1.02–2.43, p = 0.039; CC vs TT: OR = 1.72, 95% CI: 1.04–2.86, p = 0.036; TC + CC vs TT: OR = 1.62, 95% CI: 1.07–2.45, p = 0.022; and additive: OR = 1.31, 95% CI: 1.02–1.68, p = 0.035), as shown in Table 4.
The correlation between MIR17HG variants and glioma prognosis
In this study, 592 patients had complete follow-up data. The detail information for the follow-up was as following: the median, min and max follow-up time were 11 months, 2 months and 8 months, respectively. The median time to events for OS and PFS were 11 months and 8 months, respectively; total number of events for OS and DFS were 527 patients and 523 patients, respectively.
Next, we investigated the correlation between MIR17HG variants and PFS or OS of glioma by Kaplan–Meier survival method, univariate and multivariate Cox proportional hazard model. Rs17735387 was related to the PFS of glioma (Log-rank p = 0.026), as shown in Fig. 1 and Table 5. Multivariate Cox proportional hazard mode adjusted for age, sex WHO grade, surgical method, use of radiotherapy and chemotherapy showed that carriers of rs17735387 GA genotype might present a longer PFS than patients with GG genotype (HR = 0.82, 95% CI: 0.68–0.99, p = 0.042; Table 6). No statistically significant association was found between other MIR17HG polymorphisms and the prognosis of glioma.
In patients with low-grade glioma (I-II), the Kaplan–Meier method (Table 5) revealed the association between MIR17HG rs17735387 and OS (Log-rank p = 0.032, Fig. 2a) or PFS (Log-rank p = 0.013, Fig. 2b). Univariate Cox proportional hazard model presented that the GA genotype of rs17735387 might had a better OS (HR = 0.77, p = 0.042) and PFS (HR = 0.75, p = 0.024) when compared with GG genotype among patients with I-II glioma (Table 6). Moreover, the multivariate Cox proportional hazard model also displayed that a better prognosis for glioma was also seen for rs17735387-GA genotype (OS: HR = 0.75, p = 0.024 and PFS: HR = 0.73, p = 0.016). However, no association between MIR17HG polymorphisms and the prognosis of glioma in high-grade glioma patients was found.
The age and sex stratified analyses were performed to assess the association between MIR17HG polymorphisms and the prognosis of glioma. In female patients, Kaplan–Meier method (Table 5) revealed the association of rs72640334 with OS (Log-rank p = 0.035, Fig. 2c) or PFS (Log-rank p = 0.049, Fig. 2d). The results of univariate Cox proportional hazard model showed that rs72640334 was related to the poor prognosis (OS, HR = 1.49, p = 0.027 and PFS, HR = 1.48, p = 0.034, Table 6). Kaplan–Meier method (Table 5) revealed the association between rs17735387 and OS (Log-rank p = 0.002, Fig. 2e) or PFS (Log-rank p = 0.002, Fig. 2f) among patients with age ≥ 40 years. In the subgroup of patients with age ≥ 40 years, GA genotype (multivariate: OS, HR = 0.77, p = 0.047) and AA (univariate: PFS, HR = 0.036, p = 0.007; multivariate: OS, HR = 0.46, p = 0.049 and PFS, HR = 0.45, p = 0.042, Table 6) genotype of rs17735387 were associated with a better prognosis.
This study explored the possible relationship between MIR17HG variants and the occurrence and prognosis of glioma in a Chinese Han population. Our data revealed that rs7318578, rs17735387 and rs7336610 polymorphisms were associated with the increased susceptibility to glioma. We also found that rs17735387 was related to a better prognosis of patients with glioma. To our knowledge, we firstly reported that MIR17HG polymorphisms might be related to glioma susceptibility and patients’ survival.
MIR17HG gene is also called c13orf25 and Oncomir-1, which encodes a polycistronic miR-17-92 cluster encompassed six miRNAs (miR-17, miR-18a, miR-19a, miR-20a, miR-19b-1, and miR-92a-1). The miR-17-92 cluster was deregulated in glioma, indicating that these miRNA played a key role of in gliomagenesis [19, 20]. Schulte JH et al. reported that miR-17-92 cluster amplification in neuroblastomas was associated with a poor prognosis . lncRNA MIR17HG was upregulated in glioma tissues and cell lines, and acted as competing endogenous RNA (ceRNA) to sponge miR-346/miR-425-5p in regulating the malignant of glioma . Yuze Cao et al. reported that lncRNA MIR17HG-mediated ceRNA network was identified as a potential prognostic biomarker for glioblastoma . Moreover, Xue Leng et al. observed that MIR17HG was highly expressed in glioma and participated in piR-DQ590027/ lncRNA MIR17HG/miR-153(miR-377)/FOXR2 pathway which involved in regulating the permeability of glioma-conditioned normal blood-brain barrier . These results suggested that lncRNA MIR17HG could be of pathogenic importance in the development and prognosis of glioma. Several previous studies have reported the effect of MIR17HG genetic polymorphisms on the risk of various disease including tumors [24, 25], but not in glioma.
Considering the importance of MIR17HG in the carcinogenic process of glioma, we hypothesized that MIR17HG polymorphisms might also are associated with glioma development. Here, we explored the relationship between five SNPs in MIR17HG and the risk and prognosis of glioma in a Chinese Han population. We found that rs7318578 might had a higher susceptibility to glioma. The incidence rates of glioma, that is, the rate of newly diagnosed tumor, are associated with increasing age and male gender . We further analyzed whether the genotypic effects of MIR17HG on the risk of glioma were dependent on age and sex. We found that rs7318578 was related to the increase risk of glioma in the subjects with age ≥ 40 years or in females. In addition, rs17735387 and rs7336610 also had a higher susceptibility to glioma in the subgroup aged < 40 years. These indicated that the effect of MIR17HG polymorphisms on glioma occurrence might present age and sex difference. More importantly, we found that rs17735387 was related to the better prognosis of patients with glioma, particularly in low-grade glioma. Previously, rs7336610 was reported to be associated with the risk of multiple myeloma and breast cancer, while rs17735387 had no relationship with the risk and prognosis of multiple myeloma [16, 24]. These results suggested that MIR17HG polymorphisms might have a different effect on the occurrence of different cancer types. However, our findings need further studies to confirm.
Inevitably, some limitations should not be ignored. First, all individuals including glioma patients and healthy controls were from the same hospital, therefore the selection bias cannot be ruled out. Second, due to the lack of data on environmental exposure and diet, the interaction between environment and genetics needs to be further explored in larger prospective studies. Third, the effect of these SNPs on miR-17-92 cluster or lncRNA MIR17HG was not assessed.
In conclusion, we reported that MIR17HG rs7318578 might be a risk factor for the susceptibility of glioma and rs17735387 was associated with the longer survival of glioma among Chinese Han population. Our study firstly provided evidence about the effect of MIR17HG polymorphisms on the risk and prognosis of glioma, which might help to enhance the understanding of MIR17HG gene in gliomagenesis. In subsequent studies, we will continue to collect samples and follow up to further validate our findings and further explore the function of these MIR17HG SNPs in glioma in a larger sample size.
Availability of data and materials
All the data regarding the findings are available within the manuscript. Anyone who is interested in the information should contact the corresponding author.
Minor allele frequencies
Weller M, Wick W, Aldape K, Brada M, Berger M, Pfister SM, et al. Glioma. Nat Rev Dis Primers. 2015;1:15017.
Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, et al. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66:115–32.
Ostrom QT, Gittleman H, Stetson L, Virk S, Barnholtz-Sloan JS. Epidemiology of intracranial Gliomas. Prog Neurol Surg. 2018;30:1–11.
Chen H, Chen Y, Zhao Y, Fan W, Zhou K, Liu Y, et al. Association of sequence variants on chromosomes 20, 11, and 5 (20q13.33, 11q23.3, and 5p15.33) with glioma susceptibility in a Chinese population. Am J Epidemiol. 2011;173:915–22.
Wang D, Hu E, Wu P, Yuan W, Xu S, Sun Z, et al. Genetic variant near TERC influencing the risk of gliomas with older age at diagnosis in a Chinese population. J Neuro-Oncol. 2015;124:57–64.
Kinnersley B, Houlston RS, Bondy ML. Genome-wide association studies in Glioma. Cancer Epidemiol Biomark Prev. 2018;27:418–28.
Chen Y, Wu Y, Huang X, Qu P, Li G, Jin T, et al. Leukocyte telomere length: a novel biomarker to predict the prognosis of glioma patients. J Cancer Res Clin Oncol. 2015;141:1739–47.
Li S, Jin T, Zhang J, Lou H, Yang B, Li Y, et al. Polymorphisms of TREH, IL4R and CCDC26 genes associated with risk of glioma. Cancer Epidemiol. 2012;36:283–7.
Li G, Jin T, Liang H, Zhang Z, He S, Tu Y, et al. RTEL1 tagging SNPs and haplotypes were associated with glioma development. Diagn Pathol. 2013;8:83.
Olive V, Jiang I, He L. mir-17-92, a cluster of miRNAs in the midst of the cancer network. Int J Biochem Cell Biol. 2010;42:1348–54.
Molinari C, Salvi S, Foca F, Teodorani N, Saragoni L, Puccetti M, et al. miR-17-92a-1 cluster host gene (MIR17HG) evaluation and response to neoadjuvant chemoradiotherapy in rectal cancer. OncoTargets Ther. 2016;9:2735–42.
Yuan J, Tan L, Yin Z, Zhu W, Tao K, Wang G, et al. MIR17HG-miR-18a/19a axis, regulated by interferon regulatory factor-1, promotes gastric cancer metastasis via Wnt/beta-catenin signalling. Cell Death Dis. 2019;10:454.
Mao JT, Xue B, Smoake J, Lu QY, Park H, Henning SM, et al. MicroRNA-19a/b mediates grape seed procyanidin extract-induced anti-neoplastic effects against lung cancer. J Nutr Biochem. 2016;34:118–25.
Cao S, Zheng J, Liu X, Liu Y, Ruan X, Ma J, et al. FXR1 promotes the malignant biological behavior of glioma cells via stabilizing MIR17HG. J Exp Clin Cancer Res. 2019;38:37.
Chen P, Bai Y, Li Y, Yuan Y, Cheng Y, Pang J, et al. Association between polymorphisms of MIR17HG and risk of colorectal cancer in the Chinese Han population. Mol Genet Genomic Med. 2019;7:e667.
Chacon-Cortes D, Smith RA, Lea RA, Youl PH, Griffiths LR. Association of microRNA 17-92 cluster host gene (MIR17HG) polymorphisms with breast cancer. Tumour Biol. 2015;36:5369–76.
Dobashi A, Tsuyama N, Asaka R, Togashi Y, Ueda K, Sakata S, et al. Frequent BCOR aberrations in extranodal NK/T-cell lymphoma, nasal type. Genes Chromosom Cancer. 2016;55:460–71.
Ellis JA, Ong B. The MassARRAY((R)) System for Targeted SNP Genotyping. Methods Mol Biol (Clifton, NJ). 2017;1492:77–94.
Sun G, SiMa G, Wu C, Fan Y, Tan Y, Wang Z. Decreased MiR-17 in glioma cells increased cell viability and migration by increasing the expression of Cyclin D1, p-Akt and Akt. PLoS One. 2018;13:e0190515.
Wang W, Zhang A, Hao Y, Wang G, Jia Z. The emerging role of miR-19 in glioma. J Cell Mol Med. 2018;22:4611–6.
Schulte JH, Marschall T, Martin M, Rosenstiel P, Mestdagh P, Schlierf S, et al. Deep sequencing reveals differential expression of microRNAs in favorable versus unfavorable neuroblastoma. Nucleic Acids Res. 2010;38:5919–28.
Cao Y, Wang P, Ning S, Xiao W, Xiao B, Li X. Identification of prognostic biomarkers in glioblastoma using a long non-coding RNA-mediated, competitive endogenous RNA network. Oncotarget. 2016;7:41737–47.
Leng X, Ma J, Liu Y, Shen S, Yu H, Zheng J, et al. Mechanism of piR-DQ590027/MIR17HG regulating the permeability of glioma conditioned normal BBB. J Exp Clin Cancer Res. 2018;37:246.
Wu H, Huang T, Ye Z, Fu X, Hu K, Yang X. Correlation of MicroRNA 17-92 cluster host gene (MIR17HG) polymorphisms with susceptibility and prognosis for multiple myeloma. Clin Lymphoma Myeloma Leuk. 2019;19:e359–e66.
Walther C, Mayrhofer M, Nilsson J, Hofvander J, Jonson T, Mandahl N, et al. Genetic heterogeneity in rhabdomyosarcoma revealed by SNP array analysis. Genes Chromosom Cancer. 2016;55:3–15.
Ostrom QT, Bauchet L, Davis FG, Deltour I, Fisher JL, Langer CE, et al. The epidemiology of glioma in adults: a “state of the science” review. Neuro-oncology. 2014;16:896–913.
We are grateful to the individuals for their participation in this study.
This work was supported by 2016 Cultivation Fund of Hainan Medical College (HN2016-06) and 2019 Hainan Provincial Higher Education Scientific Research Project (Hnky2019-42).
Ethics approval and consent to participate
Our research was approved by the Ethics Committee of the Second Affiliated Hospital of Hainan Medical University and was in the Declaration of Helsinki. Written informed consent was obtained from each participant.
Consent for publication
The authors declare that they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Feng, J., Ouyang, Y., Xu, D. et al. Genetic variants in MIR17HG affect the susceptibility and prognosis of glioma in a Chinese Han population. BMC Cancer 20, 976 (2020). https://doi.org/10.1186/s12885-020-07417-9
- Genetic variants