Research Paper Volume 10, Issue 6 pp 1402—1414

The role of microRNA-608 polymorphism on the susceptibility and survival of cancer: a meta-analysis

Zhi-Ming Dai1, , Jian-Rui Lv1, , Kang Liu2, , Xiao-Ming Lei1, , Wei Li1, , Gang Wu1, , Xing-Han Liu2, , Yu-Xiao Zhu2, , Qian Hao2, , Zhi-Jun Dai2, ,

  • 1 Department of Anesthesiology, Second Affiliated Hospital of Xi’an Jiaotong University, Xi’an, Shaanxi 710004, China
  • 2 Department of Oncology, Second Affiliated Hospital of Xi’an Jiaotong University, Xi’an, Shaanxi 710004, China
* Equal contribution

Received: January 17, 2018       Accepted: June 10, 2018       Published: June 16, 2018      

https://doi.org/10.18632/aging.101476
How to Cite

Copyright: © 2018 Dai et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

The role of rs4919510 polymorphism in microRNA-608 (miR-608) and cancer susceptibility and prognosis remain controversial and debatable. We conducted a meta-analysis of twenty-four eligible publications on the association of rs4919510 polymorphism with cancer risk and/or prognosis. Odds ratios, hazard ratios, and 95% confidence interval were used to investigate the association between this polymorphism and susceptibility, overall survival, and recurrence-free survival of cancer. Overall, eighteen case-control studies and nine cohort studies evaluated the susceptibility and prognostic value of rs4919510 polymorphism in cancer, respectively. Pooled analysis showed that rs4919510 polymorphism was not associated with cancer risk in all five genetic models. When stratifying by different cancer sites, rs4919510 polymorphism was detected to have a significant association with a decreased risk of colorectal cancer in homozygous model (P = 0.006) and recessive model (P = 0.001), subgroup analysis also emerged a weakened correlation between rs4919510 polymorphism and an increased risk of papillary thyroid cancer in heterozygote model (P = 0.04). Furthermore, the prognosis of rs4919510 variant in cancer patients showed that rs4919510 GG genotype was significant association with poor recurrence-free survival in homozygous models (P = 0.04). The meta-analysis suggested that the microRNA-608 rs4919510 polymorphism maybe associate with a significantly decreased risk for colorectal cancer. Further investigations on larger populations are required to evaluate and confirm this relationship.

Introduction

MicroRNA (miRNA) is a growing family of naturally occurring, comprised of ~22 nucleotide, endogenous, single stranded non-coding RNA that is found in both prokaryotes and eukaryotes [1]. MiRNAs regulates the stability or translational efficiency of the targeted messenger RNAs (mRNAs) by targeting the 3' untranslated region (3´UTR) and altering the translation of genes at the post-transcription level [2]. At present, more than 1000 miRNA genes have been discovered in the human genome [3,4]. Numerous studies have indicated that miRNAs post-transcriptionally regulate a variety of biological processes, including cell proliferation, differentiation, apoptosis, and stress response. Thus, it is conceivable that the miRNA biogenesis pathway may have effects on cancer susceptibility, prognosis, and treatment response [5,6]. Single nucleotide polymorphisms (SNPs) are inherited genetic variations, which can affect gene expression or protein function. SNPs in miRNA genes may change the property of miRNAs by altering pri-miRNA/pre-miRNA processing, or by affecting the miRNA-mRNA interactions [7]. Therefore, SNPs in miRNAs may be considered as biomarkers for the diagnosis and/or prognosis of cancer.

The microRNA-608 (miR-608) gene lies within an intron of SEMA4G on the human chromosome 10q24 locus. Due to the loss of heterozygosity, the 10q24 locus has been reported to be correlated with certain types of human cancers, including brain, colorectal, prostate, and breast cancer [8-11]. The rs4919510 polymorphism, located in the mature miR-608 sequence, is a C/G single-nucleotide variation. The rs4919510 variation was considered to influence miRNA activities by binding to mir608 target sites within CD4 antigen, growth hormone receptor (GHR), retinoic X receptor beta (RXRB), and tumour protein p53 (TP53) genes with lower free energies than the wild type allele [12].

It has been hypothesised that rs4919510 polymorphism in miR-608 plays a part in carcinogenesis and clinic outcomes of cancer. Recently, several studies have revealed that the rs4919510 variation in the miR-608 gene could influence the susceptibility and/or prognosis of certain types of cancer, including colorectal [10,13-15], breast [11,16-19], thyroid [20,21], gastric [22,23], lung [24,25], head and neck cancers [26-28], and neuroblastoma [29], as well as esophageal squamous cell carcinoma [30], and hepatocellular carcinomas [31,32]. However, the results of the previous studies are inconsistent. In addition, an earlier meta-analysis found that rs4919510 polymorphism in miR-608 was significantly associated with decreased cancer risk in recessive model [33]. In recent years, there was a lot of new literature published. In this study, a phenotype meta-analysis was performed to evaluate the effects of miR-608 rs4919510 genetic polymorphism on susceptibility and prognosis of different types of cancer.

Results

Characteristics of eligible studies

Through the primary literature retrieval in Pubmed, Embase, and CNKI, 90 studies were identified for cancer susceptibility and/or prognosis assessment based on rs4919510 polymorphism of miR-608 gene. After review of titles and abstracts, we identified 39 potential articles eligible to be included for the evaluation. After retrieving the literature, two overlapping studies were found [25,34], and we excluded the earlier study. Other two of the selected studies did not meet the requirements of Hardy-Weinberg Equilibrium (HWE) [35,36], and hence twenty-four articles were assessed for eligibility to be included in our meta-analysis. The flow chart for the studies selection process is shown in Fig. 1. There were two separate groups from different ethnicities (African-Americans and Caucasians) in Ryan et al’s study [10]. Thus, we treated them as separate studies. Another report by Zheng et al. [27] also included two independent studies of prognosis. Finally, eighteen independent case-control studies from seventeen published articles were evaluated for the association of rs4919510 polymorphism with cancer risk [10,11,13,15,17-23,25,26,28-31] and nine studies from seven published articles were evaluated for the prognostic value of rs4919510 polymorphism in cancer [10,14-16,24,27,32]. All eligible studies were published between 2012 and 2018 and carried out in the USA, China, Chile, Iran, and Lithuania. The studies presented data for several types of cancer including colorectal, gastric, lung, breast, thyroid, neuroblastoma and head and neck cancers as well as hepatocellular, and oesophageal squamous cell carcinomas. The characteristics of the eligible studies are presented in Table 1 and Table 2.

Preferred reporting items for systematic reviews and meta-analyses flow diagram of the literature review process for microRNA-608 rs4919510 and cancer.

Figure 1. Preferred reporting items for systematic reviews and meta-analyses flow diagram of the literature review process for microRNA-608 rs4919510 and cancer.

Table 1. Characteristics of studies evaluating cancer susceptibility included in the meta-analysis.

First authorYearCountryCancer typeEthnicityStudy designGenotyping methodSource of controlCase/Control     Cases    ControlsP of HWENOS
GGCGCCGGCGCC
Ryan2012USACRCAACC*TaqmanPop94/1851248342895620.399
CaucasianCCTaqmanPop145/248748908711690.879
Huang2012ChinaBCAsianCCSNPstreamPop1138/19343815451926409143540.388
Wang2014ChinaHCCAsianCCMassArrayPop816/7202414151602273611320.588
Kupcinskas12014LithuaniaGCCaucasianCCRT-PCRPop363/350258825013862510.118
Kupcinskas22014LithuaniaCRCCaucasianCCRT-PCRPop192/42674713812963180.167
Wei2015ChinaPTCAsianCCMassArrayPop828/10312664281303265032020.758
Qiu2015ChinaNPCAsianCCTaqManPop906/10722554601912545322860.838
Zhang2015ChinaESCCAsianCCSNaPshotPop738/8822173841372914401510.488
Dong2015ChinaPTCAsianCCMassArrayHosp369/751136186472793701020.248
Jiang2016ChinaGCAsianCCMassArrayHosp898/9922784511652964832100.628
Miao2016ChinaHNCAsianCCIlluminaHosp576/15521772851145097622780.818
Ying2016ChinaCRCAsianCCMassArrayPop805/6182326904232505123130.157
Dai2016ChinaBCAsianCCMassArrayPop560/5831572961071832871130.988
Hashemi2016IranBCCaucasianCCPCR-RFLPPop160/1920201400431490.086
Morales2016ChileBCCaucasianCCTaqmanPop440/80740174226663104310.338
Yin2016ChinaLCAsianCCTaqManHosp575/6081582941231972691150.848
He2018ChinaNeuroblastomaAsianCCTaqmanPop393/812127190762274051790.958
CRC: colorectal cancer;BC: breast cancer; HCC:hepatocellular carcinoma;GC: gastric cancer; PTC:papillary thyroid cancer;NPC: nasopharyngeal carcinoma; ESCC: esophageal squamous cell carcinoma; LC: lung cancer; HNC: head and neck cancer; AA: African-American; CC*: case-control Pop: population based; Hosp: hospital based; HWE: Hardy-Weinberg equilibrium; NOS Newcastle-Ottawa scale.

Table 2. Characteristics of studies evaluating cancer prognosis included in the meta-analysis.

AuthorCountryType of tumorEthnicitySurvival analysisMean/Median age (SD/range)StageCasesGenotypesHR(95%CI)Median follow-up (months)
OS (Univariable)OS (Multivariable)RFS (Univariable)RFS (Multivariable)
Ryan 2012aUSACRCAAOS64.7 (11.7)I-IV94CG vs. CC0.70 (0.39-1.27)0.72 (0.40-1.30)--50.5
GG vs. CC0.38 (0.13-1.13)0.36 (0.12-1.07)--
CG/GG vs. CC0.77 (0.53-1.12)---
Ryan 2012bUSACRCCaucasianOS64.7 (11.7)I-IV145CG vs. CC1.25 (0.74-2.10)1.41 (0.81-2.44)--50.5
GG vs. CC3.54 (1.38-9.12)2.95 (1.13-7.71)--
CG/GG vs. CC1.24 (0.87-1.77)---
Xing 2012ChinaCRCChineseOS, RFS59.4 (22-90)0-IV408GC vs. CC0.99 (0.60-1.64)-0.85 (0.57-1.27)-23
GG vs. CC1.11 (0.65-1.88)-1.21 (0.80-1.81)-
CG/GG vs. CC1.03 (0.64-1.66)-1.03(0.71-1.49)-
Zheng 2013aChinaNPCChineseRFS54 (21-82)I-IV873CG vs. CC--1.41 (0.89-2.04)1.45 (0.92-2.18)37.2
GG vs. CC--2.02 (1.29-3.16)2.05 (1.35-3.21)
GC/GG vs. CC--1.61 (1.13-2.29)1.68 (1.18-2.61)
Zheng 2013bChinaNPCChineseRFS52 (19-87)I-IV828CG vs. CC--1.55 (1.03-2.15)1.62 (1.11-2.19)34.8
GG vs. CC--2.21 (1.49-3.24)2.24 (1.45-3.38)
GC/GG vs. CC--1.65 (1.17-2.39)1.71 (1.23-2.31)
Jiao 2014ChinaBCChineseOS51 (29-83)I-III196CG vs. CC1.42 (0.65-3.11)---114
GG vs. CC1.32 (0.72-2.41)---
CG/GG vs. CC1.34 (0.76-2.38)---
CC/CG vs. GG1.21 (0.61-2.40)---
Ma 2016ChinaHCCChineseOS50 (20-81)I-IV361GG vs. CC1.67 (1.11-2.49)1.61 (1.07-2.40)--53
CG/GG vs. CC1.47 (1.03-2.09)1.40 (0.98-2.00)--
CC/CG vs. GG0.76 (0.57-1.01)0.76 (0.57-1.01)--
Xia 2016ChinaNSCLCChineseOS57.12 (11.47)I-IV584CG vs. GG1.07 (0.84-1.37)0.98 (0.76-1.25)-
GG vs. CC1.00 (0.74-1.36)0.99 (0.73-1.34)
Ying 2016ChinaCRCChineseOS, RFS60.80 (13.03)0-III1218CG vs. CC0.89 (0.70-1.13)0.89 (0.70-1.14)0.92 (0.76-1.11)0.93 (0.76-1.13)26
GG vs. CC1.03 (0.75-1.41)1.10 (0.80-1.52)1.04 (0.81-1.34)1.08 (0.83-1.40)
CG/GG vs. CC0.92 (0.73-1.15)0.94 (0.75-1.18)0.94 (0.79-1.14)0.97 (0.80-1.16)
CC/CG vs. GG0.90 (0.68-1.19)0.85 (0.64-1.12)0.91 (0.72-1.15)0.88 (0.70-1.10)
G vs. C0.99 (0.85-1.16)1.02 (0.87-1.19)1.00 (0.89-1.13)1.02 (0.90-1.16)
CRC: colorectal cancer; NPC: nasopharyngeal carcinoma; BC: breast cancer; HCC:hepatocellular carcinoma; NSCLC: non-small cell lung cancer; AA: African-American; OS: overall survival; RFS: recurrence-free survival.

The association between rs4919510 polymorphism and cancer susceptibility

Overall, there were eighteen studies encompassing 10,345 cases and 14,160 controls included in our meta-analysis. The associations of the rs4919510 polymorphism in miR-608 with the risk of different types of cancers are shown in Table 3. No statistically significant association was observed in the overall population under all comparison models. In the subgroup analysis of ethnicity, a significant association with increased cancer risk was found for Asians under heterozygote comparison (OR = 1.08, 95% CI = 1.01-1.17, Fig. 2). When the analysis was stratified by different cancer types, rs4919510 polymorphism was detected significant association with a decreased colorectal cancer risk in two genetic models (GG vs. CC: OR= 0.74, 95% CI = 0.60-0.92, Fig. 3; GG vs. CC/CG: OR= 0.73, 95% CI = 0.61-0.88). Moreover, significant associations were observed between rs4919510 polymorphism and papillary thyroid cancer in heterozygote model (OR = 1.25, 95% CI = 1.01-1.54).

Table 3. Meta-analysis results for the overall and subgroup analyses of rs4919510 polymorphism and cancer susceptibility.

ComparisonsG vs. CGG vs. CCGG vs. CC/CGGG/CG vs. CCCG vs. CC
OR (95% CI)P -valueHeterogeneity
(I2, P-value)
OR (95% CI)P -valueHeterogeneity (I2, P-value)OR (95% CI)P -valueHeterogeneity (I2, P-value)OR (95% CI)P -valueHeterogeneity (I2, P-value)OR (95% CI)P -valueHeterogeneity (I2, P-value)
Overall1.01
(0.94-1.07)
0.8462%, 0.0011.03
(0.90-1.18)
0.6461%, 0.0010.97
(0.88-1.07)
0.5456%, 0.0031.05
(0.97-1.14)
0.2539%, 0.051.07
(1.00-1.14)
0.054%, 0.41
Ethnicity
Caucasian1.08
(0.89-1.32)
0.4551%, 0.091.36
(0.98-1.87)
0.060, 0.641.31
(0.96-1.80)
0.090, 0.611.05
(0.83-1.32)
0.6950%, 0.091.03
(0.88-1.20)
0.7545%, 0.12
Asian0.99
(0.93-1.07)
0.8767%, 0.0011.00
(0.86-1.16)
0.9969%, 0.0010.95
(0.86-1.05)
0.3164%, 0.0021.05
(0.96-1.16)
0.2943%, 0.061.08
(1.01-1.17)
0.030, 0.55
Tumor type
BC0.98
(0.85-1.13)
0.7856%, 0.081.06
(0.89-1.25)
0.530, 0.580.99
(0.87-1.12)
0.840, 0.390.99
(0.80-1.22)
0.9159%, 0.080.99
(0.80-1.24)
0.9650%, 0.08
CRC0.99
(0.80-1.22)
0.9150%, 0.080.74
(0.60-0.92)
0.00628%, 0.250.73
(0.61-0.88)
0.00121%, 0.280.97
(0.84-1.12)
0.6416%, 0.321.04
(0.89-1.21)
0.640, 0.73
GC1.11
(0.99-1.25)
0.080, 0.401.27
(1.00-1.62)
0.0538%, 0.201.11
(0.92-1.34)
0.2963%, 0.101.18
(0.98-1.42)
0.090, 0.851.13
(0.93-1.38)
0.220, 0.50
PTC1.07
(0.96-1.19)
0.240, 0.491.20
(0.95-1.50)
0.120, 0.471.01
(0.87-1.19)
0.850, 0.801.23
(1.00-1.51)
0.050, 0.401.25
(1.01-1.54)
0.040, 0.42
HNC1.06
(0.81-1.40)
0.6689%, 0.0031.13
(0.65-1.98)
0.660, 0.471.07
(0.77-1.48)
0.3480%, 0.031.10
(0.72-1.68)
0.6586%, 0.0091.09
(0.78-1.54)
0.6176%, 0.04
OR: odds ratio; CI: confidence interval; BC: breast cancer; CRC: colorectal cancer; GC: gastric cancer; PTC: papillary thyroid cancer;HNC: head and neck cancer.
Stratified analysis based on ethnicity for the association between microRNA-608 rs4919510 polymorphism and cancer risk using a heterozygous model (CG vs. CC). The squares and horizontal lines correspond to the study specific OR and 95% CI. The area of the squares reflects the weight (inverse of the variance). The diamond represents the summary OR and 95% CI.

Figure 2. Stratified analysis based on ethnicity for the association between microRNA-608 rs4919510 polymorphism and cancer risk using a heterozygous model (CG vs. CC). The squares and horizontal lines correspond to the study specific OR and 95% CI. The area of the squares reflects the weight (inverse of the variance). The diamond represents the summary OR and 95% CI.

Stratified analysis based on the different cancer sites for the association between microRNA-608 rs4919510 polymorphism and cancer risk using homozygote model (GG vs. CC). The squares and horizontal lines correspond to the study specific OR and 95% CI. The area of the squares reflects the weight (inverse of the variance). The diamond represents the summary OR and 95% CI.

Figure 3. Stratified analysis based on the different cancer sites for the association between microRNA-608 rs4919510 polymorphism and cancer risk using homozygote model (GG vs. CC). The squares and horizontal lines correspond to the study specific OR and 95% CI. The area of the squares reflects the weight (inverse of the variance). The diamond represents the summary OR and 95% CI.

The association between rs4919510 polymorphism and prognostic significance of cancer

Through retrieval in databases, nine studies reporting the association between rs4919510 polymorphism and the prognostic value of cancer were available. Survival data of the studies included our meta-analysis are summarised in Table 4. No significant association was observed between rs4919510 polymorphism and overall survival of cancer under univariate analysis or multivariate analysis. A significant result existing in homozygote comparison model indicates the GG allele maybe a significantly risk role in recurrence-free survival of cancer (GG vs CC: HR = 1.51, 95% CI = 1.02-2.23). However, the negative association was observed under multivariate analysis in homozygote model.

Table 4. Meta-analysis of the association between rs4919510 polymorphisms and prognosis of cancer.

Comparisons  ModelHRUnivariableHRMultivariable
HR (95% CI)p-valueHR (95% CI)p-value
OSGG vs. CC1.20 (0.90-1.60)0.221.18 (0.82-1.71)0.38
GG vs. CC/CG1.05 (0.91-1.21)0.511.12 (0.76-1.65)0.57
CG vs. CC0.93 (0.79-1.09)0.360.92 (0.75-1.14)0.45
RFSGG vs. CC1.51 (1.02-2.23)0.041.67 (1.01-2.76)0.05
GG vs. CC/CG1.24 (0.91-1.70)0.171.38 (0.91-2.09)0.13
CG vs. CC1.10 (0.83-1.47)0.491.26 (0.85-1.86)0.25
OS: overall survival; RFS: recurrence-free survival; HR: hazard ratio; CI: confidential interval; HRMultivariable, HR adjusted by gender, age, histological type, clinical stage, receipt of chemotherapy.

Publication bias

In our meta-analysis, both Begg’s funnel plot and Egger’s test were conducted to assess the publication bias. As shown in Fig. 4, the results revealed no publication bias for eligible studies under heterozygous model. The symmetry was then confirmed in Egger’s linear regression test (P = 0.21).

Begg’s funnel plot of publication bias on the relationships between rs4919510 polymorphism and cancer susceptibility under heterozygous model.

Figure 4. Begg’s funnel plot of publication bias on the relationships between rs4919510 polymorphism and cancer susceptibility under heterozygous model.

Heterogeneity and sensitivity analysis

Significant heterogeneities in the data of microRNA-608 rs4919510 polymorphism were observed in the overall meta-analysis as well as subgroup analysis, the results are shown in Table 3. Sources of heterogeneity were further studied by sensitivity analysis, the result showed that no individual study influenced the pooled OR value (Fig. 5), indicating that the results of our meta-analysis are statistically credible and robust.

Sensitivity analysis of association between microRNA-608 rs4919510 polymorphism and cancer risk.

Figure 5. Sensitivity analysis of association between microRNA-608 rs4919510 polymorphism and cancer risk.

Discussion

Aberrant miRNA function and expression could influence multiple physiological processes and participate in the occurrence, development, and prognosis of many human diseases including malignant tumours [7,37]. Furthermore, SNPs or mutations in miRNAs may directly exert their effects by binding to target mRNAs [38]. The rs4919510 polymorphism was located in the sequence of mature miR-608 and its G/C variation was common in several populations. A previous report found that the variant miR-608 was predicted to bind to mir608 target sites within immunity and defense, DNA repair, cell growth-related, and cell death-related genes with lower free energies than that of the wild type allele, for example, CD4, GHR, RXRB, and TP53 [12]. Therefore, it was assumed that the rs4919510 polymorphism possibly elicits different biological activities of these target mRNAs and is associated with the pathogenesis and prognosis of cancer as a predictive biomarker.

Our meta-analysis results showed that no significant association was observed between rs4919510 variant in miR-608 and risk of cancer in the overall population based on seventeen eligible case-control studies. However, an association was found between rs4919510 polymorphism and cancer risk in an Asian population based on the heterozygote model (P = 0.03), meanwhile, similar positive results were observed in PTC patient (P = 0.04). Therefore, rs4919510 polymorphism might be potential indicators of PTC risk.

A previous study found that rs4919510 polymorphism was not significantly associated with BC risk, but variant genotypes (GG/CG) influenced HER2-positive BC risk [11]. The results our own group proved rs4919510 alleles of miR-608 were not associated with the presence of BC a Northern Chinese population [17]. Negative results have also been reported in Chile populations [19]. Conversely, another research results revealed that rs4919510 GC genotype significantly decreased the risk of BC compared to CC genotype in an Iran population [18]. In a published pooled analysis, Wang et al [39] reported that there is no relationship between rs4919510 polymorphism and BC risk, our findings were partially in line with results from this meta-analysis.

Ryan et al. [10] and Kupcinskashas et al.’s [13] studies indicated no significant associations between rs4919510 polymorphism and CRC risk in Caucasian and African American populations. Opposite results were proved by Ying et al. who found that rs4919510 GG genotype was significantly associated with a decreased susceptibility to colorectal cancer in homozygous and recessive genetic models [15]. Remarkably, our meta-analysis results showed that rs4919510 GG genotype was significantly associated with a decreased susceptibility to CRC in homozygous and recessive genetic models. The results suggested that the presence of GG models might play a protect role in tumorigenesis in colorectal cancer. The strength of his conclusion should be confirmed on the base of more studies and sample size.

Qiu et al. [26] showed that carriers of rs4919510 CC genotypes had a reduced risk of nasopharyngeal carcinoma versus individuals who carried GG/CG genotypes in Chinese population. However, negative results were also reported in Chinese populations by Miao et al. [28]. Moreover, null significant association between this polymorphism and gastric cancer risk in Kupcinskas et al. [22] and Jiang et al.’s [23] studies. The other studies by Wang et al. [31] and Zhang et al. [30] showed that rs4919510 polymorphism was not associated with HCC and ESCC susceptibility respectively.

To date, there are three meta-analyses concerning the correlation of rs4919510 genetic susceptibility to cancer. Hu et al.’s study [40] indicated that no significant associations between rs4919510 polymorphism and overall risk of cancer or the risk of specific types of cancer. Another pooled analysis observed that rs4919510 polymorphism was significantly associated with elevated cancer risk in Chinese [41], however, the meta-analysis included some published article that the control population was not conformed the requirements of Hardy-Weinberg Equilibrium. Liu et al. [33] showed that rs4919510 variant was significantly associated with decreased cancer risk based on recessive model. Our study included more case-control studies than those two previous researches and the results challenged of their findings and provided stronger evidence.

Polymorphism of rs4919510 in miR-608 has been little explored for its effect on clinical outcomes for different types of cancer. Lin et al. [42] found that rs4919510 variation was associated with an increased renal cell carcinoma risk. Xing et al. [14] reported that the rs4919510 C allele was associated with a reduced risk of colorectal cancer and had a greater RFS and OS in a Chinese population. Interestingly, a study on colorectal cancer by Ryan et al. [10] showed that the rs4919510 GG genotype was associated with an increased risk of death in Caucasians and with a reduced risk of death in African Americans. However, negative results were reported in Asian populations in different types of cancer [15,16,24]. Furthermore, Zheng et al. [27] found that rs4919510 polymorphism was not associated with RFS but with OS of nasopharyngeal carcinoma in Asians.

As far as we know, this is the first meta-analysis regarding the rs4919510 polymorphism in miR-608 and its clinical prognosis significance in cancer. The results indicated that rs4919510 variant was significantly associated with the RFS of cancer under homozygote models (P = 0.04) in Asians by univariable analysis, however, which was no longer remained significant after adjustment for age, gender, histology and clinic stage and so on in multivariable model. No significant association was observed between this polymorphism and the OS of cancer. This might result from the different types of cancers, differences in clinical trials including different ethical backgrounds, different baseline characteristics, or various follow-up time periods used in each study. Furthermore, the exact functional mechanism by which rs4919510 polymorphism in miR-608 may affect cancer prognosis, there are several plausible hypotheses such as the polymorphism altering the transcription of target genes and reducing the sensitivity of the cancer cell to chemotherapeutic drugs [7,43].

Our study provides added evidence that the rs4919510 gene polymorphism may be associated with cancer susceptibility and prognosis. However, some potential limitations of our meta-analysis should be acknowledged. First, the sample size was relatively small in some studies, which may reduce the statistical power of the analysis and needs to be interpreted with caution. Second, in our meta-analysis, some cancers only had one study included, such as hepatocellular carcinoma, esophageal squamous cell carcinoma, lung cancer and neuroblastoma, which may lead to heterogeneity in quantitative analysis. Third, we perform stratified analysis by population, most of which are Chinese, and type of cancer, most of which are BC and CRC, rather than other baseline characteristics. Fourth, there was significant heterogeneity in some comparison models in our meta-analysis. Even with the use the appropriate meta-analytical techniques with the random-effects model, the heterogeneity might affect the precision of the overall results.

In conclusion, this meta-analysis suggests that the rs4919510 polymorphism in miR-608 is associated with cancer risk and prognosis. Further well-designed, prospective investigations including rs4919510 polymorphism and cancer susceptibility and prognosis with larger sample sizes are required in order to validate our findings.

Methods

Literature search and selection criteria

This meta-analysis followed the Preferred Reporting Items for Systematic Reviews and Meta-analyses (PRISMA) criteria [44]. We searched PubMed, Embase, and Chinese National Knowledge Infrastructure (CNKI) databases for studies published prior to March 2018 (last search: March 30, 2018). The search was conducted with and without medical subject heading (MeSH) terms for ‘‘MicroRNA-608/rs4919510’’, ‘‘polymorphism or variant or mutation or genotype’’, and ‘‘cancer’’. The languages were limited to English and the subjects were human. The reference lists of searched literatures were read manually to complete our investigation, if necessary. It is necessary to select appraisal quality tools to evaluate the selected literature [45]. The Newcastle-Ottawa Scale (NOS) was used for assessing the quality of case-control study, the results were showed in Table 1. The NOS score of all articles are not less than 6 scores that mean each selected literature was a high-quality study.

Studies that fit into our meta-analysis fell into one of the following categories: (1) evaluation of the association of rs4919510 polymorphism and cancer susceptibility and/or prognosis; (2) provision of sufficient data to estimate the odds ratio (OR) or hazard ratio (HR) and 95% confidence intervals (CI) according to rs4919510 polymorphism; (3) case-control design; (4) controls without any evidence of malignant disease. The following were the exclusion criteria: (1) repeats studies, reviews, and abstracts; (2) the study did not have a control group; (3) design based on family; (4) the control population was not conformed the requirements of Hardy-Weinberg Equilibrium. The diagnosis of each cancer in eligible study was made based on the histopathological evidence and those with secondary tumour or metastasized cancer from other organs were excluded. Two investigators (Z-MD and J-RL) retrieved the literature and extracted data from each eligible study independently.

Data extraction and synthesis

For each eligible study investigating the association between rs4919510 polymorphism and cancer risk, the following information was extracted by the two independent researchers: first author, year of publication, country, tumour type, ethnicity, source of control, size of the study population, genotyping methods, as well as Hardy-Weinberg value of controls. For studies examining the association between rs4919510 polymorphism and its prognostic value in cancer, the following information was collected from the eligible articles: first author, year of publication, origin of the study population, type of cancer, method of survival analysis, HRs with 95% CIs, as well as the follow-up time. Different ethnic groups were categorised as Caucasian, Asian, and African. Disagreements were resolved by discussion among all authors.

Statistical analysis

ORs and 95% CIs were used to measure cancer risk associated with the rs4919510 polymorphism. Five different ORs were calculated: (1) GG versus CC (homozygous carriers), (2) CG versus CC (heterozygous carriers), (3) GG+ CG versus CC (dominant model), (4) GG versus CC + CG (recessive model), and (5) G versus C (allele comparison). The Hardy-Weinberg equilibrium (HWE) for the controls was assessed by the Chi-square test in each study. HRs with their 95% CIs was combined to evaluate the effective value of the rs4919510 polymorphism on the prognosis of cancer. The significance of the pooled ORs/HRs was determined by the Z test. If the study did not report the HRs, the Engauge Digitizer version 4.1 was used to read the Kaplane-Meier curves to estimate the HRs and the 95% CIs.

Statistical heterogeneity between studies was measured by the Chi-square-based Q statistic and evaluated with the I2 test; higher I2 values means higher levels of heterogeneity (I2 = 75-100%: extreme heterogeneity; I2 = 50-75%: large heterogeneity; I2 = 25-50%: moderate heterogeneity; I2 < 25%: no heterogeneity) [46]. When the p value was greater than 0.10, the fixed-effects model was used, otherwise the random-effects model was used. The Egger’s test and Begg’s funnel plot were adopted to assess publication bias. Sensitivity analysis determined whether the individual data in fact have a major effect on the results of the review. All analysis was performed by the Review Manager version 5.3 (The Cochrane Collaboration, Oxford, United Kingdom) and the Stata version 14.0 (Stata Corp, College Station, TX).

Abbreviations

OR: odds ratio; HR: hazard ratio, CI: confidence interva; RFS: recurrence-free surviva; HWE: Hardy-Weinberg equilibriu; PRISMA: Preferred Reporting Items for Systematic Reviews and Meta-analyse; CNKI: Chinese National Knowledge Infrastructur; MeSH: medical subject headin; NOS: Newcastle-Ottawa Scale.

Author Contributions

ZJD, JRL and ZMD designed the research, analyzed the data, wrote the paper. XML, ZMD and KL reviewed drafts of the paper, and participated in its design and coordination. WL and GW contributed reagents/materials/analysis tools; DZM, XHL and YXZ prepared figures and/or tables. KL and QH reviewed drafts of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Funding

This study was supported by National Natural Science Foundation, China (No. 81471670); China Postdoctoral Science Foundation (No. 2014M560791; 2015T81037); the Fundamental Research Funds for the Central Universities, China (No. 2014qngz-04) and Science and Technology Plan of Innovation Project, Shaanxi province, China (No. 2015KTCL03-06).

References

  • 1. Lee RC, Feinbaum RL, Ambros V. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993; 75:843–54. https://doi.org/10.1016/0092-8674(93)90529-Y [PubMed]
  • 2. Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB. Prediction of mammalian microRNA targets. Cell. 2003; 115:787–98. https://doi.org/10.1016/S0092-8674(03)01018-3 [PubMed]
  • 3. Zamore PD, Haley B. Ribo-gnome: the big world of small RNAs. Science. 2005; 309:1519–24. https://doi.org/10.1126/science.1111444 [PubMed]
  • 4. Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009; 136:215–33. https://doi.org/10.1016/j.cell.2009.01.002 [PubMed]
  • 5. Ohtsuka M, Ling H, Doki Y, Mori M, Calin GA. MicroRNA Processing and Human Cancer. J Clin Med. 2015; 4:1651–67. https://doi.org/10.3390/jcm4081651 [PubMed]
  • 6. Iorio MV, Casalini P, Piovan C, Braccioli L, Tagliabue E. Breast cancer and microRNAs: therapeutic impact. Breast. 2011 (Suppl 3); 20:S63–70. https://doi.org/10.1016/S0960-9776(11)70297-1 [PubMed]
  • 7. Ryan BM, Robles AI, Harris CC. Genetic variation in microRNA networks: the implications for cancer research. Nat Rev Cancer. 2010; 10:389–402. https://doi.org/10.1038/nrc2867 [PubMed]
  • 8. Daido S, Takao S, Tamiya T, Ono Y, Terada K, Ito S, Ouchida M, Date I, Ohmoto T, Shimizu K. Loss of heterozygosity on chromosome 10q associated with malignancy and prognosis in astrocytic tumors, and discovery of novel loss regions. Oncol Rep. 2004; 12:789–95. [PubMed]
  • 9. Kim JH, Dhanasekaran SM, Mehra R, Tomlins SA, Gu W, Yu J, Kumar-Sinha C, Cao X, Dash A, Wang L, Ghosh D, Shedden K, Montie JE, et al. Integrative analysis of genomic aberrations associated with prostate cancer progression. Cancer Res. 2007; 67:8229–39. https://doi.org/10.1158/0008-5472.CAN-07-1297 [PubMed]
  • 10. Ryan BM, McClary AC, Valeri N, Robinson D, Paone A, Bowman ED, Robles AI, Croce C, Harris CC. rs4919510 in hsa-mir-608 is associated with outcome but not risk of colorectal cancer. PLoS One. 2012; 7:e36306. https://doi.org/10.1371/journal.pone.0036306 [PubMed]
  • 11. Huang AJ, Yu KD, Li J, Fan L, Shao ZM. Polymorphism rs4919510:C>G in mature sequence of human microRNA-608 contributes to the risk of HER2-positive breast cancer but not other subtypes. PLoS One. 2012; 7:e35252. https://doi.org/10.1371/journal.pone.0035252 [PubMed]
  • 12. Landi D, Gemignani F, Barale R, Landi S. A catalog of polymorphisms falling in microRNA-binding regions of cancer genes. DNA Cell Biol. 2008; 27:35–43. https://doi.org/10.1089/dna.2007.0650 [PubMed]
  • 13. Kupcinskas J, Bruzaite I, Juzenas S, Gyvyte U, Jonaitis L, Kiudelis G, Skieceviciene J, Leja M, Pauzas H, Tamelis A, Pavalkis D, Kupcinskas L. Lack of association between miR-27a, miR-146a, miR-196a-2, miR-492 and miR-608 gene polymorphisms and colorectal cancer. Sci Rep. 2014; 4:5993. https://doi.org/10.1038/srep05993 [PubMed]
  • 14. Xing J, Wan S, Zhou F, Qu F, Li B, Myers RE, Fu X, Palazzo JP, He X, Chen Z, Yang H. Genetic polymorphisms in pre-microRNA genes as prognostic markers of colorectal cancer. Cancer Epidemiol Biomarkers Prev. 2012; 21:217–27. https://doi.org/10.1158/1055-9965.EPI-11-0624 [PubMed]
  • 15. Ying HQ, Peng HX, He BS, Pan YQ, Wang F, Sun HL, Liu X, Chen J, Lin K, Wang SK. 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. https://doi.org/10.18632/oncotarget.12422 [PubMed]
  • 16. Jiao L, Zhang J, Dong Y, Duan B, Yu H, Sheng H, Huang J, Gao H. Association between miR-125a rs12976445 and survival in breast cancer patients. Am J Transl Res. 2014; 6:869–75. [PubMed]
  • 17. Dai ZM, Kang HF, Zhang WG, Li HB, Zhang SQ, Ma XB, Lin S, Wang M, Feng YJ, Liu K, Liu XH, Xu P, Dai ZJ. The Associations of Single Nucleotide Polymorphisms in miR196a2, miR-499, and miR-608 With Breast Cancer Susceptibility: A STROBE-Compliant Observational Study. Medicine (Baltimore). 2016; 95:e2826. https://doi.org/10.1097/MD.0000000000002826 [PubMed]
  • 18. Hashemi M, Sanaei S, Rezaei M, Bahari G, Hashemi SM, Mashhadi MA, Taheri M, Ghavami S. miR-608 rs4919510 C>G polymorphism decreased the risk of breast cancer in an Iranian subpopulation. Exp Oncol. 2016; 38:57–59. [PubMed]
  • 19. Morales S, Gulppi F, Gonzalez-Hormazabal P, Fernandez-Ramires R, Bravo T, Reyes JM, Gomez F, Waugh E, Jara L. Association of single nucleotide polymorphisms in Pre-miR-27a, Pre-miR-196a2, Pre-miR-423, miR-608 and Pre-miR-618 with breast cancer susceptibility in a South American population. BMC Genet. 2016; 17:109. https://doi.org/10.1186/s12863-016-0415-0 [PubMed]
  • 20. Wei WJ, Wang YL, Li DS, Wang Y, Wang XF, Zhu YX, Pan XD, Wang ZY, Wu Y, Jin L, Wang JC, Ji QH. Association study of single nucleotide polymorphisms in mature microRNAs and the risk of thyroid tumor in a Chinese population. Endocrine. 2015; 49:436–44. https://doi.org/10.1007/s12020-014-0467-8 [PubMed]
  • 21. Dong G, Zhang R, Xu J, Guo Y. Association between microRNA polymorphisms and papillary thyroid cancer susceptibility. Int J Clin Exp Pathol. 2015; 8:13450–57. [PubMed]
  • 22. Kupcinskas J, Wex T, Link A, Leja M, Bruzaite I, Steponaitiene R, Juzenas S, Gyvyte U, Ivanauskas A, Ancans G, Petrenkiene V, Skieceviciene J, Kupcinskas L, Malfertheiner P. Gene polymorphisms of micrornas in Helicobacter pylori-induced high risk atrophic gastritis and gastric cancer. PLoS One. 2014; 9:e87467. https://doi.org/10.1371/journal.pone.0087467 [PubMed]
  • 23. Jiang J, Jia ZF, Cao DH, Wu YH, Sun ZW, Cao XY. Association of the miR-146a rs2910164 polymorphism with gastric cancer susceptibility and prognosis. Future Oncol. 2016; 12:2215–26. https://doi.org/10.2217/fon-2016-0224 [PubMed]
  • 24. Lingzi X, Zhihua Y, Xuelian L, Yangwu R, Haibo Z, Yuxia Z, Baosen Z. Genetic variants in microRNAs predict non-small cell lung cancer prognosis in Chinese female population in a prospective cohort study. Oncotarget. 2016; 7:83101–14. https://doi.org/10.18632/oncotarget.13072 [PubMed]
  • 25. Yin Z, Cui Z, Ren Y, Xia L, Wang Q, Zhang Y, He Q, Zhou B. Association between polymorphisms in pre-miRNA genes and risk of lung cancer in a Chinese non-smoking female population. Lung Cancer. 2016; 94:15–21. https://doi.org/10.1016/j.lungcan.2016.01.013 [PubMed]
  • 26. Qiu F, Yang L, Zhang L, Yang X, Yang R, Fang W, Wu D, Chen J, Xie C, Huang D, Zhou Y, Lu J. Polymorphism in mature microRNA-608 sequence is associated with an increased risk of nasopharyngeal carcinoma. Gene. 2015; 565:180–86. https://doi.org/10.1016/j.gene.2015.04.008 [PubMed]
  • 27. Zheng J, Deng J, Xiao M, Yang L, Zhang L, You Y, Hu M, Li N, Wu H, Li W, Lu J, Zhou Y. A sequence polymorphism in miR-608 predicts recurrence after radiotherapy for nasopharyngeal carcinoma. Cancer Res. 2013; 73:5151–62. https://doi.org/10.1158/0008-5472.CAN-13-0395 [PubMed]
  • 28. Miao L, Wang L, Zhu L, Du J, Zhu X, Niu Y, Wang R, Hu Z, Chen N, Shen H, Ma H. Association of microRNA polymorphisms with the risk of head and neck squamous cell carcinoma in a Chinese population: a case-control study. Chin J Cancer. 2016; 35:77. https://doi.org/10.1186/s40880-016-0136-9 [PubMed]
  • 29. He J, Zou Y, Liu X, Zhu J, Zhang J, Zhang R, Yang T, Xia H. Association of Common Genetic Variants in Pre-microRNAs and Neuroblastoma Susceptibility: A Two-Center Study in Chinese Children. Mol Ther Nucleic Acids. 2018; 11:1–8. https://doi.org/10.1016/j.omtn.2018.01.003 [PubMed]
  • 30. Zhang P, Wang J, Lu T, Wang X, Zheng Y, Guo S, Yang Y, Wang M, Kolluri VK, Qiu L, Shen F, Fan L, Li J, et al. miR-449b rs10061133 and miR-4293 rs12220909 polymorphisms are associated with decreased esophageal squamous cell carcinoma in a Chinese population. Tumour Biol. 2015; 36:8789–95. https://doi.org/10.1007/s13277-015-3422-2 [PubMed]
  • 31. Wang R, Zhang J, Ma Y, Chen L, Guo S, Zhang X, Ma Y, Wu L, Pei X, Liu S, Wang J, Hu H, Liu J. Association study of miR‑149 rs2292832 and miR‑608 rs4919510 and the risk of hepatocellular carcinoma in a large‑scale population. Mol Med Rep. 2014; 10:2736–44. https://doi.org/10.3892/mmr.2014.2536 [PubMed]
  • 32. Ma XP, Yu G, Chen X, Xiao Q, Shi Z, Zhang LY, Chen H, Zhang P, Ding DL, Huang HX, Saiyin H, Chen TY, Lu PX, et al. MiR-608 rs4919510 is associated with prognosis of hepatocellular carcinoma. Tumour Biol. 2016; 37:9931–42. https://doi.org/10.1007/s13277-016-4897-1 [PubMed]
  • 33. Liu H, Zhou Y, Liu Q, Xiao G, Wang B, Li W, Ye D, Yu S. Association of miR-608 rs4919510 polymorphism and cancer risk: a meta-analysis based on 13,664 subjects. Oncotarget. 2017; 8:37023–31. https://doi.org/10.18632/oncotarget.9509 [PubMed]
  • 34. Yin Z, Cui Z, Guan P, Li X, Wu W, Ren Y, He Q, Zhou B. Interaction between Polymorphisms in Pre-MiRNA Genes and Cooking Oil Fume Exposure on the Risk of Lung Cancer in Chinese Non-Smoking Female Population. PLoS One. 2015; 10:e0128572. https://doi.org/10.1371/journal.pone.0128572 [PubMed]
  • 35. Li D, Zhu G, Di H, Li H, Liu X, Zhao M, Zhang Z, Yang Y. Associations between genetic variants located in mature microRNAs and risk of lung cancer. Oncotarget. 2016; 7:41715–24. https://doi.org/10.18632/oncotarget.9566 [PubMed]
  • 36. Ranjbar R, Chaleshi V, Aghdaei HA, Morovvati S. Investigating the association between miR-608 rs4919510 and miR-149 rs2292832 with Colorectal Cancer in Iranian Population. MicroRNA. 2018; 07. https://doi.org/10.2174/2211536607666180206145540 [PubMed]
  • 37. Duan S, Mi S, Zhang W, Dolan ME. Comprehensive analysis of the impact of SNPs and CNVs on human microRNAs and their regulatory genes. RNA Biol. 2009; 6:412–25. https://doi.org/10.4161/rna.6.4.8830 [PubMed]
  • 38. Negrini M, Nicoloso MS, Calin GA. MicroRNAs and cancer--new paradigms in molecular oncology. Curr Opin Cell Biol. 2009; 21:470–79. https://doi.org/10.1016/j.ceb.2009.03.002 [PubMed]
  • 39. Wang J, Kong X, Xing Z, Wang X, Zhai J, Fang Y, Gao J. A meta-analysis: is there any association between MiR-608 rs4919510 polymorphism and breast cancer risks? PLoS One. 2017; 12:e0183012. https://doi.org/10.1371/journal.pone.0183012 [PubMed]
  • 40. Hu Y, Yu CY, Wang JL, Guan J, Chen HY, Fang JY. MicroRNA sequence polymorphisms and the risk of different types of cancer. Sci Rep. 2014; 4:3648. https://doi.org/10.1038/srep03648 [PubMed]
  • 41. Wu S, Yuan W, Shen Y, Lu X, Li Y, Tian T, Jiang L, Zhuang X, Wu J, Chu M. The miR-608 rs4919510 polymorphism may modify cancer susceptibility based on type. Tumour Biol. 2017; 39:1010428317703819. https://doi.org/10.1177/1010428317703819 [PubMed]
  • 42. Lin J, Horikawa Y, Tamboli P, Clague J, Wood CG, Wu X. Genetic variations in microRNA-related genes are associated with survival and recurrence in patients with renal cell carcinoma. Carcinogenesis. 2010; 31:1805–12. https://doi.org/10.1093/carcin/bgq168 [PubMed]
  • 43. Sanoff HK, Sargent DJ, Green EM, McLeod HL, Goldberg RM. Racial differences in advanced colorectal cancer outcomes and pharmacogenetics: a subgroup analysis of a large randomized clinical trial. J Clin Oncol. 2009; 27:4109–15. https://doi.org/10.1200/JCO.2009.21.9527 [PubMed]
  • 44. Liberati A, Altman DG, Tetzlaff J, Mulrow C, Gøtzsche PC, Ioannidis JP, Clarke M, Devereaux PJ, Kleijnen J, Moher D. The PRISMA statement for reporting systematic reviews and meta-analyses of studies that evaluate health care interventions: explanation and elaboration. J Clin Epidemiol. 2009; 62:e1–34. https://doi.org/10.1016/j.jclinepi.2009.06.006 [PubMed]
  • 45. Zeng X, Zhang Y, Kwong JS, Zhang C, Li S, Sun F, Niu Y, Du L. The methodological quality assessment tools for preclinical and clinical studies, systematic review and meta-analysis, and clinical practice guideline: a systematic review. J Evid Based Med. 2015; 8:2–10. https://doi.org/10.1111/jebm.12141 [PubMed]
  • 46. Higgins JP, Thompson SG. Quantifying heterogeneity in a meta-analysis. Stat Med. 2002; 21:1539–58. https://doi.org/10.1002/sim.1186 [PubMed]