Abstract

The objective of this study was to conduct a meta-analysis to systematically summarize and investigate the association of miRNA-124 rs531564, miRNA-218 rs11134527, miRNA-146a rs2910164, miRNA-196a2 rs11614913, and miRNA-499 rs3746444 polymorphisms with cervical cancer. A systematic review was performed to identify relevant studies using Embase and PubMed databases. A chi-square-based -test combined with the inconsistency index () was used to check the heterogeneity between studies. A total of six case-control studies on rs2910164 and rs11614913, 4 studies on rs3746444 and rs11134527, and three studies on rs531564 were included. No evidence of association was found between miR-146a rs2910164, miR-196a2 rs11614913, miRNA-499 rs3746444, and miR-218 rs11134527 polymorphisms and cervical cancer risk in all the genetic models. The miR-124 rs531564 polymorphism was associated with a statistically increased risk of cervical cancer in a homozygote model (CC vs. GG: , 95% CI: 1.40-5.91, ), dominant model (GC/CC vs. GG: , 95% CI: 1.07-1.80, ), and recessive model (CC vs. GC/GG: , 95% CI: 1.58-3.23, ). However, this finding should be interpreted with caution for limited samples and heterogeneity. Large-scale and well-designed studies are needed to validate our result.

1. Introduction

Cervical cancer (CC) is the fourth most common cancer in females, with 570,000 new cases and 311,000 deaths estimated for 2018 worldwide [1]. Cervical cancer is a multistep process involving the transformation of the normal cervical epithelium to cervical intraepithelial neoplasia that is subsequently transformed to cervical cancer [2]. Due to cervical cytology screening, cervical cancer can be detected at an early stage, greatly reducing the incidence and mortality. However, cervical cancer is still one of the deadliest female-specific cancers due to its tendency to metastasize and recur after treatment. It is well established that the persistence of Human Papillomavirus (HPV) infection is the main cause of cervical cancer and is indeed deemed as a necessary cause for the disease. However, evidence suggests that genetic risk factors also play a crucial role in the pathogenesis of cervical cancer [35].

MicroRNAs (miRNAs) are short nonprotein-coding small RNAs of approximately 23 nucleotides that mediate posttranscriptional regulation through base pairing to partially complementary sites of the target mRNA and promote translational repression or messenger RNA degradation [6, 7]. Increasing evidences show that single nucleotide polymorphisms (SNPs) or variation in miRNA sequence could potentially alter various biological processes by influencing target selection of miRNAs and is closely associated with various human diseases [8]. The relationship between genetic mutations in microRNA and cancer susceptibility has been getting increasing attention. Several recent studies have indicated that miRNA-related SNPs can remarkably alter the biogenesis and function of the corresponding miRNAs and may function as oncogenic or antioncogenic molecules in human cancers [9, 10].

Recently, though a number of studies have been conducted to investigate the association of miR-124 G>C (rs531564), miR-218 A>G (rs11134527), miR-146a G>C (rs2910164), miR-196a2 T>C (rs11614913), and miR-499 A>G (rs3746444) polymorphisms with cervical cancer [2, 1120], the results were inconclusive. Therefore, a meta-analysis was performed to evaluate the association between these five SNPs and the risk of cervical cancer.

2. Methods and Materials

2.1. Publication Search

The PubMed and Embase databases were searched by two independent investigators covering all papers published up to December 27, 2020. The following terms were used: “cervical or CIN”, “microRNA or miRNA or microRNAs”, and “polymorphism or polymorphisms”. The references of the retrieved articles and review articles were also searched manually for further relevant studies. All eligible studies met the following criteria: (1) full-text study, (2) case-control design, (3) analysis of the relationship between miRNA polymorphisms and the risk of cervical cancer, (4) studies focusing on human beings, and (5) sufficient data to calculate risk estimates: the odds ratio (OR) with 95% confidence interval (CI) and a value. The selection of the study was completed independently by two investigators (Liu and Dong) according to the inclusion criteria.

2.2. Data Extraction

Each publication was sought for the following information: the first author’s surname, the publication year, country of origin, ethnicity, genotyping method, the number of cases and controls, and Hardy-Weinberg equilibrium (HWE). According to the amount of the sample, we defined sample size as the large sample size (≥200) and small sample size (<200).

Two reviewers (Liu and Dong) extracted the eligible research data repeatedly. Disagreement between two reviewers was discussed with another reviewer (Wang) until a consensus is reached.

When Hardy-Weinberg equilibrium (HWE) in the control was not reported, the test was used to evaluate the HWE of polymorphism in the control group. If , it is considered deviated from HWE.

2.3. Methodological Quality Assessment

The quality of the selected studies was assessed according to the methodological quality assessment scale which was mentioned in prior meta-analysis [22] (see Table S1).

In this scale, five items were evaluated, including (1) the representativeness of the subjects, (2) the ascertainment of cervical cancer, (3) the control of the genotyping methods, (4) HWE in controls, and (5) assessment of association. The quality score ranged from 0 to 9, and a high score indicated good quality of the study. Studies with a score less than 5 were removed from the subsequent analyses.

2.4. Statistics

We conducted our meta-analysis based on the PRISMA checklists and followed the guideline [23]. The strength of association was assessed between miRNA polymorphisms and the risk of cervical cancer by crude OR corresponding to 95% CI [24]. Pooled ORs were performed for allelic comparison (miR-146a: G versus C, miR-196a: T versus C, miR-499: A versus G, miR-218: A versus G, and miR-124: G versus C), heterozygote model (miR-146a: GC versus GG, miR-196a: TC versus TT, miR-499: AG versus AA, miR-218: AG versus AA, and miR-124: GC versus GG), homozygote model (miR-146a: CC versus GG, miR-196a: CC versus TT, miR-499: GG versus AA, miR-218: GG versus AA, and miR-124: CC versus GG), dominant model (miR-146a: CC+GC versus GG, miR-196a: TC+CC versus TT, miR-499: AG+GG versus AA, miR-218: AG+GG versus AA, and miR-124: GC+CC versus GG), and recessive model (miR-146a: CC versus CG+GG, miR-196a: CC versus TC+TT, miR-499: GG versus AG+AA, miR-218: GG versus AG+AA, and miR-124: CC versus GC+GG), respectively.

A chi-square-based -test combined with the inconsistency index () was used to check the heterogeneity between studies [25]. The random-effects model was used when heterogeneity tests yielded significant results ( or ) [26, 27]. Besides, subgroup analyses were stratified by sample size and HWE. Publication bias was assessed by visual inspection of funnel plots and using Egger’s test. Sensitivity analyses were performed to identify individual studies’ effect on pooled results and test the reliability of results. value of Egger’s test <0.05 was considered representative of statistically significant publication bias. All statistical analyses were carried out using the STATA 12.0 software (Stata Corp., College Station, TX, USA).

3. Results

3.1. Study Characteristics

A total of 247 studies were acquired from Embase and PubMed databases (Embase: 98 and PubMed: 149). The detailed screening process is shown in Figure 1. After reading the titles and the abstracts, 235 articles were excluded, of which 61 were duplicate ones, 149 had no relation to this topic, 22 were reviews or meta-analyses, and 3 were conference abstracts. Finally, 12 eligible case-control studies were included in our meta-analysis (Table 1) [1021].

In six studies, genotype frequencies of rs531564, rs11134527, rs2910164, rs11614913, and rs3746444 were presented separately, so each of them was considered as a separate study in this meta-analysis. Therefore, all the six included studies containing 2291 cases and 2850 controls for rs2910164, six studies containing 2291 cases and 2850 controls for rs11614913, four studies involving 1554 cases and 1962 controls for rs3746444, four studies containing 3167 cases and 3132 controls for rs11134527, and three studies containing 975 cases and 950 controls for rs531564 were finally analyzed in our meta-analysis. For these five miRNA in this meta-analysis, the subjects in all included studies were of the Asian population.

The TaqMan assay and polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) were used for genotyping to determine the SNPs in the included studies. HWE of genotype distribution in the controls was tested in all studies, and two articles were not consistent with HWE. And according to the methodology quality assessment, all the studies were above a score of 6.0 and recruited into the following analyses.

3.2. Association between miRNA-146a rs2910164 Polymorphism and Cervical Cancer Susceptibility

The association strength between miRNA-146a rs2910164 polymorphism and cervical cancer risk is shown in Table 2. As shown in Table 2, significant association was not identified in any genetic model (C vs. G: , 95% CI 0.90-1.30, , Figure 2; GC vs. GG: , 95% CI 0.76-1.15, ; CC vs. GG: OR =0.85, 95% CI 0.57-1.26, ; GC/GG vs. CC: , 95% CI 0.70-1.15, ; and GG vs. GC/CC: , 95% CI 0.68-1.19, ).

Next, subgroup analysis was conducted according to sample size and HWE. In the large sample size and HWE group, a significant association was also not identified in any genetic model. And the heterogeneity is still high. When we delete three articles (Shizhi Wan, Shruti SrivaStava, and Guange Chen), heterogeneity is significantly reduced (Table 3). A significantly increased risk of cervical cancer susceptibility was identified in allelic comparison (C vs. G: , 95% CI: 1.19-1.56, ), while there is a significantly decreased risk of cervical cancer susceptibility in heterozygote comparison (GC vs. GG: , 95% CI: 0.60-0.97, ), homozygote model (CC vs. GG: , 95% CI: 0.40-0.70, ), dominant model (GC/CC vs. GG: , 95% CI: 0.54-0.85, ), and recessive model (CC vs. GC/GG: , 95% CI: 0.53-0.82, ).

3.3. Association between miR-196a2 rs11614913 Polymorphism and Cervical Cancer Susceptibility

Significant association was not identified in any genetic model (T vs. C: , 95% CI: 0.67-1.04, , Figure 3; TC vs. TT: , 95% CI: 0.90-1.58, ; CC vs. TT: , 95% CI: 0.94-2.10, ; TC/CC vs. TT: , 95% CI: 0.93-1.75, ; and CC vs. TC/TT: , 95% CI: 0.97-1.29, ) when all eligible studies were pooled.

Next, subgroup analysis was conducted according to sample size and HWE. In large sample size, heterogeneity is significantly reduced. A significantly increased risk of cervical cancer susceptibility was identified in homozygote comparison (CC vs. TT: , 95% CI: 1.00-1.41, ). In the HWE group, heterogeneity is also significantly reduced. But significant association was not identified in any genetic model.

3.4. Association between miR-499 rs3746444 Polymorphism and Cervical Cancer Susceptibility

Significant association was not identified in any genetic model (G vs. A: , 95% CI: 0.50-1.15, , Figure 4; AG vs. AA: , 95% CI: 0.82-2.20, ; GG versus AA: , 95% CI: 0.61-3.30, ; AG/GG vs. AA: , 95% CI: 0.83-2.42, ; and GG vs. AG/AA: , 95% CI: 0.68-2.35; ).

Next, subgroup analysis was conducted according to sample size and HWE. In the large sample size and HWE group, a significant association was also not identified in any genetic model. The heterogeneity in the HWE group is still high, and in a large sample size, the heterogeneity is high except the homozygote model and recessive model. When we delete the two articles (Bin Zhou and Shruti SrivaStava), heterogeneity is significantly reduced in all models (Table 3). However, significant association was also not identified in any genetic model.

3.5. Association between miR-218 rs11134527 Polymorphism and Cervical Cancer Susceptibility

Significant association was not identified in any genetic model (G vs. A: , 95% CI: 0.83-1.15, ; AG vs. AA: , 95% CI: 0.89-1.29, , Figure 5; GG vs. AA: , 95% CI: 0.73-1.40, ; AG/GG vs. AA: , 95% CI: 0.86-1.31, ; and GG vs. AG/AA: , 95% CI: 0.76-1.23, ).

The sample size of the 4 included articles is all large size and subgroup analysis was conducted according to HWE. In the HWE group, heterogeneity is also significantly reduced. But significant association was not identified in any genetic model. When we delete two articles (Li Chuanyin and Guange Chen), heterogeneity is significantly reduced in all models (Table 3). We observed a significantly increased risk of cervical cancer susceptibility in allelic comparison (G versus A: , 95% CI: 1.02-1.21, ), while there is a significantly decreased risk of cervical cancer susceptibility in the homozygote model (GG vs. AA: , 95% CI: 0.66-0.94, ) and recessive model (GG vs. AG/AA: , 95% CI: 0.68-0.94, ). Significant association was not identified in heterozygote comparison (AG vs. AA: , 95% CI: 0.85-1.10; ) and dominant model (GG/AG vs. AA: , 95% CI: 0.81-1.04, ).

3.6. Association between miR-124 rs531564 Polymorphism and Cervical Cancer Susceptibility

We observed a significantly increased risk of cervical cancer susceptibility in the homozygote model (CC vs. GG: , 95% CI: 1.40-5.91, , Figure 6), dominant model (GC/CC vs. GG: , 95% CI: 1.07-1.80, ), and recessive model (CC vs. GC/GG: , 95% CI: 1.58-3.23, ), while there is a significantly decreased risk of cervical cancer susceptibility in allelic comparison (C vs. G: , 95% CI: 0.40-0.83, ). Significant association was not identified in heterozygote comparison (GC vs.GG: , 95% CI: 0.97-1.66, ).

Next, subgroup analysis was conducted according to sample size. In a large sample size, a significantly increased risk of cervical cancer susceptibility could be identified in the homozygote model (CC vs. GG: , 95% CI: 1.04-12.63, ) and recessive model (CC vs. GC/GG: , 95% CI: 1.51-3.31, ), while there is a significantly decreased risk of cervical cancer susceptibility in allelic comparison (C vs. G: , 95% CI: 0.32-0.66; ). Significant association was not identified in heterozygote comparison (GC vs.GG: , 95% CI: 0.97-1.66, ) and dominant model (GC/CC vs. GG: , 95% CI: 0.90-10.85, ).

3.7. Sources of Heterogeneity

Since significant heterogeneity was found for miR-146 (rs2910164), miR-196a2 (rs11614913), miR-499 (rs3746444), and miR-218(rs11134527), we conducted a subgroup analysis by sample size, and the results showed that the sample size is the source of heterogeneity of miR-196a2 (rs11614913) in each genetic model. However, for other SNPs, heterogeneity was still high in large sample size. However, when we delete some articles, we found that the heterogeneity is significantly reduced.

3.8. Sensitivity Analysis

The influence of each study on the pooled ORs in this meta-analysis was examined by excluding one study at a time using STATA 12.0 software. The results showed that no significant alteration in the pooled ORs was found in any of the genetic models for the five SNPs.

3.9. Publication Bias

Due to the limited number of studies included (<10), no funnel plot or Egger’s test was performed.

4. Discussion

In this meta-analysis, the relationship between these five SNPs in microRNAs (rs2910164, rs11614913, rs3746444, rs11134527, and rs531564) and the risk of cervical cancer was evaluated from 11 published studies. We demonstrated that no evidence was found for the association between miRNA-146a rs2910164, miR-196a2 rs11614913, miRNA-499 rs3746444, and miR-218 rs11134527 and cervical cancer risk in any genetic models. The miR-124 rs531564 polymorphism was associated with a statistically increased risk of cervical cancer in the homozygote model, dominant model, and recessive model. However, due to limited samples and heterogeneity, this finding should be interpreted with caution.

Recent studies have demonstrated that miRNAs are abnormally expressed in cervical cancer and play a role in tumor generation [28]. miR-146a (rs2910164) involves a G>C nucleotide substitution, which changes the G:U pair to C:U [29]. Different locations may cause the function of miR-146a to be different, and the binding affinity to the 3-UTR region of the target gene may also be different [30, 31]. Furthermore, it is reported that miR-146a can restrict the NF-κB signaling pathway mediated by IRAK1 and TRAF6, which is an important signaling factor involved in the occurrence of cancers [32, 33]. Previous studies have shown that miR-124 is closely related to the occurrence and development of cervical cancer [34, 35]. Evidence has been shown that miR-124 rs531564 polymorphism affects miRNA processing, and the G allele affects the circular secondary structure of pre-miR-124 and changes the expression of other miRNAs [36].

In 2013, Yamamoto et al. found that the expression of miR-218 in cervical cancer tissues was significantly lower than that of adjacent noncancerous tissues. The restoration of miR-218 significantly inhibited the proliferation, migration, and invasion of cancer cells [37]. Laminin-5 β3 (LAMB3) has been verified as a transcriptional target of miR-218 [38], and the expression of LAMB3 is increased in the presence of the HPV-16 E6 oncoprotein, and this effect is mediated by miR-218 [19]. The miR-196a2 is located on the human chromosome 12q13.13. In 2016, Torruella-Loran et al. reported that rs11614913 in miR-196a2 has the function of regulating the expression of several genes involved in cancer [39]. Results of a functional study showed that the mutation reduced the efficiency of miRNA precursor processing into its mature form and reduced the ability to regulate target genes [40]. miR-499 is one of the microRNAs that play an important role in posttranscriptional regulation by regulating multiple genes and signal transduction pathways. The miR-499 (rs3746444) polymorphism leads to mismatches from A:U to G:U in the stem structure of the miR-499 precursor. The presence of this mismatch will affect Sox6 and Rod1 genes, which are important for the occurrence of cancers [41, 42].

Significant heterogeneity was found in this meta-analysis for association of miR-218 rs11134527, miR-146a rs2910164, miR-196a2 rs11614913, miR-499 rs3746444, and cervical cancer. Therefore, we conducted a subgroup analysis based on sample size and HWE. When we restricted the sample size to the large sample size group and the HWE group, heterogeneity is significantly reduced for miR-196a2 rs11614913. However, the small sample size group of miR-196a2 includes only two studies, and the large sample size group includes three studies. Therefore, further research is needed to confirm these results. Significant heterogeneity was also found for the association of miRNA-146a rs2910164, miR-499 rs3746444, miR-218 rs11134527, and cervical cancer in the large sample size group and the HWE group. However, when we delete some researches in each analysis, heterogeneity is significantly reduced, suggesting that these researches contributed to the source of heterogeneity. However, owing to the limited articles included, further evaluation of the results is still needed.

In the sensitivity analysis, no significant changes were found after omitting one study at a time, indicating that our meta-analysis results are relatively stable and credible.

Our meta-analysis has several advantages. First, this is the first meta-analysis of the association between five miRNA polymorphisms and the susceptibility to cervical cancer. In addition, according to the methodological quality assessment, all included studies are of high quality. Moreover, there are no restrictions in the literature search, so the selection bias is well controlled. Our meta-analysis has some limitations. Firstly, the number of included studies is too few. Secondly, heterogeneity was detected in miR-146a rs2910164, miR-196a2 rs11614913, miR-499 rs3746444, and miR-218 rs11134527. After deleting some studies in each analysis, heterogeneity is significantly reduced. The pooled ORs became significant without evidence of heterogeneity. Thirdly, the included studies were all Asian ethnicity; therefore, it is uncertain whether these results can be generalized to other populations. Fourth, publication bias might also exist, because only published studies were included in this meta-analysis, and results that are not statistically significant usually have less chance of publication.

5. Conclusion

In summary, our results suggested that miR-124 rs531564 polymorphism is significantly associated with increased risk of cervical cancer, while miR-146a rs2910164, miR-196a2 rs11614913, miR-499 rs3746444, and miR-218 rs2292832 polymorphisms may not be associated with the susceptibility of cervical cancer. However, since few studies have been included, there is not enough data to fully confirm the association between cervical cancer and the five miRNA polymorphisms, so the results should be interpreted with caution.

Data Availability

The data supporting this meta-analysis are from previously reported studies which have been cited.

Conflicts of Interest

The authors declare that there are no conflicts of interest associated with the manuscript.

Acknowledgments

We would like to acknowledge the Department of Obstetrics and Gynecology, Xuanwu Hospital, Capital Medical University, for helping me at work.

Supplementary Materials

Table S1: scale for quality assessment. (Supplementary Materials)