Research article

Association between miR-124-1 rs531564 polymorphism and risk of cancer: An updated meta-analysis of case-control studies

Abdolkarim Moazeni-Roodi1, Mohammad Hashemi2[*]

1Department of Clinical Biochemistry, Iranshahr University of Medical Sciences, Iranshahr, Iran

2Department of Clinical Biochemistry, School of Medicine, Zahedan University of Medical Sciences, Zahedan, Iran

EXCLI J 2018;17:Doc608

 

Abstract

Many studies examined the association between miR-124-1 rs531564 polymorphism and the risk of some human cancers, but the findings remain controversial. This update meta-analysis aimed to validate the association between rs531564 polymorphism of miR-124-1 and cancer risk. Eligible studies including 6,502 cancer cases and 7,213 controls were documented by searching Web of Science, PubMed, Scopus, and Google scholar databases. Pooled odds ratios (ORs) with 95 % confidence intervals (CIs) were estimated to quantitatively evaluate the association between rs531564 variant and cancer risk. The results indicated that rs531564 variant significantly decreased the risk of cancer in homozygous codominant (OR=0.54, 95 % CI=0.43-0.69, p<0.00001, GG vs CC), dominant (OR=0.84, 95 % CI=0.72-0.99, p=0.03, CG+GG vs CC), recessive (OR=0.65, 95 % CI=0.54-0.78, p<0.00001, GG vs CG+CC), and allele (OR=0.84, 95 % CI=0.73-0.96, p=0.008, G vs C) genetic model. Stratified analysis by cancer type revealed that rs531564 variant was associated with gastric cancer, cervical cancer, esophageal squamous cell carcinoma and colorectal cancer risk. In summary, the findings of this meta-analysis support an association between miR-124-1 rs531564 polymorphism and cancer risk. Larger and well-designed studies are required to estimate this association in detail.

Keywords: miR-124-1, rs531564, polymorphism, cancer, meta-analysis

Introduction

Cancer is one of the main leading cause of morbidity and mortality worldwide (Global Burden of Disease Cancer Collaboration, 2015[11]). In 2016, approximately 17.2 million new cancer cases and 8.9 million deaths occurred worldwide (Global Burden of Disease Cancer Collaboration, 2018[10]). It has been proposed that the complex interaction of various genetic loci and diverse environmental factors play a role in cancer development (Borek, 1993[2]; Lichtenstein et al., 2000[18]). Despite physical disparities, all human populations are 99 % genetically identical, and the remaining 1 % genetic variations is responsible for human diversity (International HapMap Consortium, 2007[14]; Ryan et al., 2010[25]). Single-nucleotide polymorphisms (SNPs) contribute to phenotypic differences both within and among populations (Omrani et al., 2014[23]).

MicroRNA (miRNA) are a class of noncoding RNAs consisting of 18-25 nucleotides in length that bind to the 3`untranslated region (3'UTR) of target mRNAs to regulate gene expression (Ryan et al., 2010[25]). Variations in miRNAs genes, including pri-miRNAs, pre-miRNAs, and mature miRNAs, impact on miRNAs biogenesis, processing, target binding, and expression level of mature miRNAs (Mishra et al., 2008[22]). Preceding studies have shown that miRNAs play a crucial role in various tumor-associated biological processes, including proliferation, metastasis, apoptosis and differentiation (He et al., 2013[12]; Liu et al., 2013[19]; Ge et al., 2016[9]).

In human, miR-124 is encoded by three miRNA genes including MIR124-1 (8p23.1), MIR124-2 (8q12.3), and MIR124-3 (20q13.33). A functional polymorphism rs531564 located in the pri-miRNA region of the miR124-1 gene affect the expression levels of the mature miR-124 (Qi et al., 2012[24]). To date, several epidemiological studies inspected the association between miR-124-1 rs531564 polymorphism and the risk of various cancer including gastric cancer (Zhou et al., 2012[37]; Asgarpour et al., 2017[1]; Singh et al., 2017[28]), cervical cancer (Wu and Zhang, 2014[31]; Xiong et al., 2014[32]; Chuanyin et al., 2017[5]), breast cancer (Ma et al., 2013[20]; Ying et al., 2016[34]; Danesh et al., 2018[6]), renal cell carcinoma (Liang et al., 2017[17]), osteosarcoma (Shi et al., 2016[27]), esophageal squamous cell carcinoma (ESCC) (Yin et al., 2013[33]; Zhang et al., 2014[36]; Wu et al., 2018[31]), colorectal cancer (Gao et al., 2015[8]), but the findings are still controversial. Therefore, we performed an updated meta-analysis to find out the impact of rs11134527 polymorphism on cancer risk.

Methods

Literature search

A comprehensive search in Web of Science, PubMed, Scopus, and Google Scholar databases was conducted for all articles describing an association between miR-124-1 rs531564 polymorphism and cancer risk published up to June 08, 2018. The search strategy was “cancer OR carcinoma, tumor OR neoplasms”, AND “miR-124-1 OR microRNA-124-1 OR miRNA-124-1” AND “polymorphism OR mutation OR variant OR rs531564”. Relevant studies included the meta-analysis if they met the following inclusion criteria: 1) Original case-control studies that evaluated the miR-124-1 polymorphism and cancer risk; 2) studies provided necessary information of the genotype frequencies of miR-124-1 rs531564 variant in both cases and controls. The exclusion criteria were: 1) conference abstract, case reports, reviews, duplication data; 2) insufficient genotype information provided.

Data extraction

The authors independently searched the literatures, extracted the relevant data and finally discussed disagreement. The following data were recorded from each study including the first author's name, publication year, country, ethnicity of participants, cancer type, genotyping methods of miR-124-1 rs531564 polymorphism, number of genotypes in case-control groups and result of the HWE test (Table 1(Tab. 1), References in Table 1: Asgarpour et al., 2017[1]; Chuanyin et al., 2017[5]; Danesh et al., 2018[6]; Gao et al., 2015[8]; Liang et al., 2017[17]; Ma et al., 2013[20]; Shi et al., 2016[27]; Singh et al., 2017[28]; Wu et al., 2014[31]; Wu et al., 2018[30]; Xiong et al., 2014[32]; Yin et al., 2013[33]; Ying et al., 2016[34]; Zhang et al., 2014[36]; Zhou et al., 2012[37]).

Statistical analysis

Meta-analysis was achieved by Revman 5.3 software (Version 5.3. Copenhagen: The Nordic Cochrane Centre, The Cochrane Collaboration, 2014) and STATA 14.1 software (Stata Corporation, College Station, TX, USA). Hardy-Weinberg equilibrium (HWE) for each study was calculated by the χ2 test. The association between miR-124-1 rs531564 polymorphism and cancer risk was assessed by pooled odds ratios (ORs) and their 95 % confidence intervals (CIs) for codominant (CG vs CC and GG vs CC), dominant (CG+GG vs CC), recessive (GG vs CG+CC), overdominant (CG vs CC+GG) and the allelic (G vs C) genetic inheritance models. The significance of the pooled OR was assessed by the Z-test, and P<0.05 was considered to be statistically significant. The choice of using fixed or random effects model was determined by the results of the between-study heterogeneity test, which was measured using the Q test and I2 statistic. If the test result was I2 ≥ 50 % or PQ < 0.1, indicating the presence of heterogeneity, the random effect model was selected; otherwise, the fixed-effects model was used.

Begg's funnel plot was conducted under all inheritance models to evaluate the publication bias and the asymmetric plots implied potential publication bias. The degree of asymmetry was tested using Egger's test and p < 0.05 was considered significant publication bias.

Sensitivity analysis was performed to evaluate the stability of the studies on the pooled ORs. A single study in the analysis was neglected each time to calculate the outcomes again.

Results

Study characteristics

Totally 15 case-control studies including 6,502 cancer cases and 7,213 controls were included in the meta-analyses. Table 1(Tab. 1) shows the characteristics and relevant data of the included studies.

Main analysis results

In the current meta-analysis of 15 eligible studies, the findings support an association between miR-124-1 rs531564 polymorphism and cancer risk. The rs531564 variant significantly decreased the risk of cancer in homozygous codominant (OR=0.54, 95 % CI=0.43-0.69, p<0.00001, GG vs CC), dominant (OR=0.84, 95 % CI=0.72-0.99, p=0.03, CG+GG vs CC), recessive (OR=0.65, 95 % CI=0.54-0.78, p<0.00001, GG vs CG+CC), and allele (OR=0.84, 95 % CI=0.73-0.96, p=0.008, G vs C) inheritance model (Figure 1(Fig. 1) and Table 2(Tab. 2); References in Figure 1: Asgarpour et al., 2017[1]; Chuanyin et al., 2017[5]; Danesh et al., 2018[6]; Gao et al., 2015[8]; Liang et al., 2017[17]; Ma et al., 2013[20]; Shi et al., 2016[27]; Singh et al., 2017[28]; Wu et al., 2014[31]; Wu et al., 2018[30]; Xiong et al., 2014[32]; Yin et al., 2013[33]; Ying et al., 2016[34]; Zhang et al., 2014[36]; Zhou et al., 2012[37]).

Subgroup analysis by cancer type

Stratified analysis of miR-124-1 rs531564 polymorphism was achieved by cancer type (Table 3(Tab. 3)). The data implied that rs531564 variant increased the risk of gastric cancer in overdominant (OR=1.27, 95 % CI=1.02-1.58, p=0.03, CG vs GG+CC) inheritance model. The rs531564 variant was associated with significantly decrease in risk of cervical cancer in codominant, dominant, recessive and allele inheritance model (Table 3(Tab. 3)). Furthermore, the variant significantly decreased the risk of ESCC in recessive and allele models. The rs531564 variant was significantly associated with protection against colorectal in recessive model (Table 3(Tab. 3)). No significant association was found between rs531564 variant and breast cancer risk.

Heterogeneity and publication bias

Heterogeneity among the findings included in the meta-analysis is shown in Table 2(Tab. 2). The data showed heterogeneity existed between studies.

The potential publication bias was estimated using a Begg's funnel plot and Egger's test. Neither Begg's funnel plot nor Egger's test detected any obvious evidence of publication bias in analyses for all genetic models (Table 2(Tab. 2) and Figure 2(Fig. 2)).

Sensitivity analysis

We executed sensitivity analysis to evaluate the effect of a specific study on the overall estimate. The relevant pooled ORs showed no significant change appeared when each study was ignored, one at a time, from the overall meta-analysis in homozygous codominant, dominant, recessive, and allele models (Figure 3(Fig. 3), References in Figure 3: Asgarpour et al., 2017[1]; Chuanyin et al., 2017[5]; Danesh et al., 2018[6]; Gao et al., 2015[8]; Liang et al., 2017[17]; Ma et al., 2013[20]; Shi et al., 2016[27]; Singh et al., 2017[28]; Wu et al., 2014[31]; Wu et al., 2018[30]; Xiong et al., 2014[32]; Yin et al., 2013[33]; Ying et al., 2016[34]; Zhang et al., 2014[36]; Zhou et al., 2012[37]). This indicates that the results of this meta-analysis are relatively stable and reliable.

Discussion

It has been well known that miRNAs is involved in carcinogenesis as tumor suppressor gene or oncogene (Calin et al., 2004[4]; Ryan et al., 2010[25]). Dysregulation of miRNAs contributes to the initiation and progression of human malignancies (Shen et al., 2015[26]; He et al., 2018[13]; Skjefstad et al., 2018[29]). It has been documented that miR-124 is a tumor suppressor miRNA in many cancers (Jin et al., 2017[15]; Yuan et al., 2017[35]; Cai et al., 2018[3]; Ma et al., 2018[21]).

An increasing number of studies have focused on associations between miR-124-1 rs531564 polymorphism and various cancer susceptibility (Zhou et al., 2012[37]; Ma et al., 2013[20]; Yin et al., 2013[33]; Wu and Zhang, 2014[31]; Xiong et al., 2014[32]; Zhang et al., 2014[36]; Gao et al., 2015[8]; Shi et al., 2016[27]; Ying et al., 2016[34]; Asgarpour et al., 2017[1]; Chuanyin et al., 2017[5]; Liang et al., 2017[17]; Singh et al., 2017[28]; Danesh et al., 2018[6]; Wu et al., 2018[30]), but the findings were inconsistent. We performed an updated meta-analysis of 15 case-control studies to find out the impact of rs531564 polymorphism of miR-124-1 gene on overall cancer risk.

The findings of our meta-analysis showed that the rs531564 polymorphism was significantly associated with protection against cancer in homozygous codominant, dominant, recessive and allele inheritance model. Our findings are in agreement with the results of two meta-analyses regarding the impact of miR-124-1 rs531564 variant on cancer risk (Fang et al., 2015[7]; Li et al., 2015[16]). The meta-analysis performed by Li et al. (2015[16]) enrolled 4 case-control studies and the findings revealed that rs531564 polymorphism significantly reduced cancer risk. The other study conducted by Fang et al. (2015[7]) with 5 case-control studies also showed that the pri-miR1241 rs531564 polymorphism significantly reduced cancer risk.

We performed stratified analysis by cancer type and the findings revealed that the rs531564 variant was significantly associated with gastric cancer, cervical cancer ESCC, colorectal and breast cancer risk.

There are several limitations in our meta-analysis that should be addressed. First, heterogeneity existed between some studies. It can be supposed that the heterogeneity probably derived from difference of ethnicity, source of control, status and cancer type. Second, the languages of the studies were limited to English. Third, we did not evaluate potential gene-environmental interactions due to lack of relevant data across the included studies. Finally, all subjects are of Asian descent, so our meta-analysis was limited to Asian population and lack of other ethnicities.

In conclusion, our meta-analysis proposed that miR-124-1 rs531564 polymorphism may be an important protective factor for cancer in Asians. Additional well designed case-control studies with larger sample sizes are required to validate the findings.

Conflict of interest

The authors have declared that no competing interests exist.

 

References

1. Asgarpour J, Mehrzad J, Yazdi SAT. MicroRNA binding site polymorphisms are associated with the development of gastric cancer. J Cell Mol Res. 2017;9:73-7.
2. Borek C. Molecular mechanisms in cancer induction and prevention. Environ Health Perspect. 1993;101 (Suppl 3):237-45.
3. Cai WL, Huang WD, Li B, Chen TR, Li ZX, Zhao CL, et al. MicroRNA-124 inhibits bone metastasis of breast cancer by repressing Interleukin-11. Mol Cancer. 2018;17:9.
4. Calin GA, Sevignani C, Dumitru CD, Hyslop T, Noch E, Yendamuri S, et al. Human microRNA genes are frequently located at fragile sites and genomic regions involved in cancers. Proc Natl Acad Sci U S A. 2004;101:2999-3004.
5. Chuanyin L, Xiaona W, Zhiling Y, Yu Z, Shuyuan L, Jie Y, et al. The association between polymorphisms in microRNA genes and cervical cancer in a Chinese Han population. Oncotarget. 2017;8:87914-27.
6. Danesh H, Hashemi M, Bizhani F, Hashemi SM, Bahari G. Association study of miR-100, miR-124-1, miR-218-2, miR-301b, miR-605, and miR-4293 polymorphisms and the risk of breast cancer in a sample of Iranian population. Gene. 2018;647:73-8.
7. Fang C, Zeng H, Li A, Xu X, Long X. Association of the pri-miR-124-1 rs531564 polymorphism with cancer risk: A meta-analysis. Mol Clin Oncol. 2015;3:892-6.
8. Gao XR, Wang HP, Zhang SL, Wang MX, Zhu ZS. Pri-miR-124 rs531564 polymorphism and colorectal cancer risk. Sci Rep. 2015;5:14818.
9. Ge L, Zheng B, Li M, Niu L, Li Z. MicroRNA-497 suppresses osteosarcoma tumor growth in vitro and in vivo. Oncol Lett. 2016;11:2207-12.
10. Global Burden of Disease Cancer Collaboration, Fitzmaurice C, Akinyemiju TF, Al Lami FH, Alam T, Alizadeh-Navaei R, et al. Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 29 cancer groups, 1990 to 2016: a systematic analysis for the Global Burden of Disease Study. JAMA Oncol. 2018;epub. doi: 10.1001/jamaoncol.2018.2706.
11. Global Burden of Disease Cancer Collaboration, Fitzmaurice C, Dicker D, Pain A, Hamavid H, Moradi-Lakeh M, et al. The global burden of cancer 2013. JAMA Oncol. 2015;1:505-27.
12. He J, Zhang W, Zhou Q, Zhao T, Song Y, Chai L, et al. Low-expression of microRNA-107 inhibits cell apoptosis in glioma by upregulation of SALL4. Int J Biochem Cell Biol. 2013;45:1962-73.
13. He RQ, Yang X, Liang L, Chen G, Ma J. MicroRNA-124-3p expression and its prospective functional pathways in hepatocellular carcinoma: A quantitative polymerase chain reaction, gene expression omnibus and bioinformatics study. Oncol Lett. 2018;15:5517-32.
14. International HapMap Consortium, Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, et al. A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007;449:851-61.
15. Jin H, Li Q, Cao F, Wang SN, Wang RT, Wang Y, et al. MiR-124 inhibits lung tumorigenesis induced by K-ras mutation and NNK. Mol Ther Nucleic Acids. 2017;9:145-54.
16. Li WJ, Wang Y, Gong Y, Tu C, Feng TB, Qi CJ. MicroRNA-124 rs531564 polymorphism and cancer risk: a meta-analysis. Asian Pac J Cancer Prev. 2015;16:7905-9.
17. Liang G, Liu Q, Sun J, Liu W, Shen W. The rs20541 G> A polymorphism in the interleukin 13 gene is associated with a decreased risk of renal cell carcinoma in a Chinese population. Int J Clin Exp Pathol. 2017;10:708-16.
18. Lichtenstein P, Holm NV, Verkasalo PK, Iliadou A, Kaprio J, Koskenvuo M, et al. Environmental and heritable factors in the causation of cancer - analyses of cohorts of twins from Sweden, Denmark, and Finland. N Engl J Med. 2000;343:78-85.
19. Liu Z, Gersbach E, Zhang X, Xu X, Dong R, Lee P, et al. MiR-106a represses the Rb tumor suppressor p130 to regulate cellular proliferation and differentiation in high-grade serous ovarian carcinoma. Mol Cancer Res. 2013;11:1314-25.
20. Ma F, Zhang P, Lin D, Yu D, Yuan P, Wang J, et al. There is no association between microRNA gene polymorphisms and risk of triple negative breast cancer in a Chinese Han population. PLoS One. 2013;8:e60195.
21. Ma J, Weng L, Wang Z, Jia Y, Liu B, Wu S, et al. MiR-124 induces autophagy-related cell death in cholangiocarcinoma cells through direct targeting of the EZH2-STAT3 signaling axis. Exp Cell Res. 2018;366:103-13.
22. Mishra PJ, Mishra PJ, Banerjee D, Bertino JR. MiRSNPs or MiR-polymorphisms, new players in microRNA mediated regulation of the cell: Introducing microRNA pharmacogenomics. Cell Cycle. 2008;7:853-8.
23. Omrani M, Hashemi M, Eskandari-Nasab E, Hasani SS, Mashhadi MA, Arbabi F, et al. Hsa-mir-499 rs3746444 gene polymorphism is associated with susceptibility to breast cancer in an Iranian population. Biomark Med. 2014;8:259-67.
24. Qi L, Hu Y, Zhan Y, Wang J, Wang BB, Xia HF, et al. A SNP site in pri-miR-124 changes mature miR-124 expression but no contribution to Alzheimer's disease in a Mongolian population. Neurosci Lett. 2012;515:1-6.
25. Ryan BM, Robles AI, Harris CC. Genetic variation in microRNA networks: the implications for cancer research. Nat Rev Cancer. 2010;10:389-402.
26. Shen B, Zhang Y, Yu S, Yuan Y, Zhong Y, Lu J, et al. MicroRNA-339, an epigenetic modulating target is involved in human gastric carcinogenesis through targeting NOVA1. FEBS Lett. 2015;589:3205-11.
27. Shi ZW, Wang JL, Zhao N, Guan Y, He W. Single nucleotide polymorphism of hsa-miR-124a affects risk and prognosis of osteosarcoma. Cancer Biomark. 2016;17:249-57.
28. Singh DK, Zhang WB, Xu Y, Yin J, Gu HY, Jiang PC. Hsa-miR-34b/c rs4938723 T>C, pri-miR-124-1 rs531564 C>G, pre-miR-125a rs12975333 G>T and hsa-miR-423 rs6505162 C>A polymorphisms and the risk of gastric cardia adenocarcinoma. Int J Clin Exp Med. 2017;10:14919-26.
29. Skjefstad K, Johannessen C, Grindstad T, Kilvaer T, Paulsen EE, Pedersen M, et al. A gender specific improved survival related to stromal miR-143 and miR-145 expression in non-small cell lung cancer. Sci Rep. 2018;8:8549.
30. Wu F, Li M, You W, Ji Y, Cui X, Hu J, et al. A genetic variant in miR-124 decreased the susceptibility to esophageal squamous cell carcinoma in a Chinese Kazakh population. Genet Test Mol Biomark. 2018;22:29-34.
31. Wu H, Zhang J. MiR-124 rs531564 polymorphism influences genetic susceptibility to cervical cancer. Int J Clin Exp Med. 2014;7:5847-51.
32. Xiong X, Cheng J, Liu X, Tang S, Luo X. [Correlation analysis between miR-124 rs531564 polymorphisms and susceptibility to cervical cancer]. Nan Fang Yi Ke Da Xue Xue Bao. 2014;34:210-3.
33. Yin J, Wang X, Zheng L, Shi Y, Wang L, Shao A, et al. Hsa-miR-34b/c rs4938723 T>C and hsa-miR-423 rs6505162 C>A polymorphisms are associated with the risk of esophageal cancer in a Chinese population. PLoS ONE. 2013;8:e80570.
34. Ying HQ, Peng HX, He BS, Pan YQ, Wang F, Sun HL, et al. MiR-608, pre-miR-124-1 and pre-miR26a-1 polymorphisms modify susceptibility and recurrence-free survival in surgically resected CRC individuals. Oncotarget. 2016;7:75865-73.
35. Yuan L, Li S, Zhou Q, Wang D, Zou D, Shu J, et al. MiR-124 inhibits invasion and induces apoptosis of ovarian cancer cells by targeting programmed cell death 6. Oncol Lett. 2017;14:7311-7.
36. Zhang J, Huang X, Xiao J, Yang Y, Zhou Y, Wang X, et al. Pri-miR-124 rs531564 and pri-miR-34b/c rs4938723 polymorphisms are associated with decreased risk of esophageal squamous cell carcinoma in Chinese populations. PLoS ONE. 2014;9:e100055.
37. Zhou Y, Du WD, Chen G, Ruan J, Xu S, Zhou FS, et al. Association analysis of genetic variants in microRNA networks and gastric cancer risk in a Chinese Han population. J Cancer Res Clin Oncol. 2012;138:939-45.
 
 

Figure 1: The forest plot for relationship between miR-124-1 rs531564 polymorphism and cancer susceptibility for CG vs CC (A), GG vs CC (B), CG+GG vs CC (C), GG vs CG+CC (D), CG vs CC+GG (E), and G vs C (F)

Figure 2: The funnel plot for the test of publication bias. The funnel plot for CG vs CC (A), GG vs CC (B), CG+GG vs CC (C), GG vs CG+CC (D), CG vs CC+GG (E), and G vs C (F)

Figure 3: Sensitivity analyses for studies on miR-124-1 rs531564 polymorphism using different genetic models. Sensitivity analyses for CG vs CC (A), GG vs CC (B), CG+GG vs CC (C), GG vs CG+CC (D), CG vs CC+GG (E), and G vs C (F)

 

Table 1: Characteristics of all studies included in the meta-analysis

Table 2: The pooled ORs and 95 % CIs for the association between miR-124-1 rs531564 polymorphism and cancer susceptibility.

Table 3: Stratified analysis of miR-124-1 rs531564 variant on susceptibility to cancer

[*] Corresponding Author:

Professor of Clinical Biochemistry PhD Mohammad Hashemi, Dept. of Clinical Biochemistry, School of Medicine, Zahedan University of Medical Sciences, Zahedan, Iran; Tel: +98 541 3235122, eMail: mhd.hashemi@gmail.com, hashemim@zaums.ac.ir