Meta-analysis of Hsa-mir-499 polymorphism (rs3746444) for cancer risk: evidence from 31 case-control studies

Background MicroRNAs (miRNAs) are a family of endogenous, small and non-coding RNAs that regulate gene expression negatively at the post-transcriptional level by suppressing translation or degrading target mRNAs, and are involved in diverse biological and pathological processes. Single nucleotide polymorphisms (SNPs) which are located in the miRNA-coding genes may participate in the process of development and diseases by altering the expression of mature miRNA. Recent studies investigating the association between hsa-mir-499 polymorphism (rs3746444) and cancer risk have yielded conflicting results. Methods In this meta-analysis, we conducted a search of case–control studies on the associations of SNP rs3746444 with susceptibility to cancer in electronic databases. A total of 31 studies involving 12799 cases and 14507 controls were retrieved and the strength of the association was estimated by pooled odds ratios (ORs) and 95% confidence intervals (CIs). Hardy-Weinberg equilibrium (HWE) was assessed by the goodness-of-fit chi-square test in controls. Subgroup analyses were done by racial descent and cancer type. Publication bias of literatures was evaluated by visual inspection of funnel plots and the linear regression asymmetry test by Egger et al. Sensitivity analysis was conducted by excluding one study at a time to examine the influence of individual data set on the pooled ORs. Results Overall, significant association between rs3746444 polymorphism and susceptibility to cancer was identified in TC versus TT and TC/CC versus TT (dominant) models. In the stratified analyses, increased risks were found in Asians, but not in Caucasians in all comparison models tested. Moreover, significant association with an increased risk was found in Chinese population. Also, much higher significant association with increased cancer risks were found in Iranian population. In different cancer types, a decreased risk was found in esophageal cancer. Conclusion Our meta-analysis suggested that hsa-mir-499 rs3746444 T > C polymorphism is associated with the risk of cancer in Asians, mainly in Iranian and Chinese population. However, rs3746444 T > C polymorphism is negatively associated with the risk of esophageal cancer. Electronic supplementary material The online version of this article (doi:10.1186/s12881-014-0126-1) contains supplementary material, which is available to authorized users.


Background
In early 1990s, microRNAs (miRNAs) were first discovered through analysis of developmental timing mutants in C. elegans by several groups, simultaneously [1,2]. In the last decade, the study of miRNA biology has attracted remarkable attention, resulting in rapid advances, and revealing miRNAs as key gene regulators in diverse biological pathways [3]. miRNAs participate in regulation of stem cell functions [4], development [5], drug resistance [6], metabolism and metabolic disorders [7], cardiovascular [8] and malignant diseases [9,10].
miRNAs are hairpin-derived RNAs~20-24 nucleotides long, which post-transcriptionally repress the expression of target genes usually by binding to the 3′-untranslated region (3′ UTR) of messenger RNA (mRNA) in a broad range of organisms in both normal physiological and disease contexts [11,12]. The majority of human miRNA loci are located within intronic regions and are transcribed by RNA polymerase II as part of their hosting transcription units [13]. Typically, following RNA polymerase II-mediated transcription, the long primary precursors (pri-miRNAs) are cleaved by the nuclear RNase III (Drosha) to release~70 nt pre-miRNAs [14]. The resulting transcripts adopt a stem-loop structure and are then exported to the cytoplasm where they are subsequently processed by another RNase III (Dicer) to generate mature double-stranded~22 nt miRNAs [15]. Subsequently, one strand of this duplex is incorporated into an Argonautecontaining RNA-induced silencing complex (RISC), resulting in the translational repression and/or degradation of their target mRNAs [16].
Recent evidences support that single nucleotide polymorphisms (SNPs), occurring in DNA sequences of miRNA-coding genes or in miRNA-binding site in mRNAs, contribute to the biogenesis and functions of miRNAs [17]. Gain/loss-of-function of miRNA polymorphisms may result in enhancing of the combination of the miRNA to the targets or losing control of the mRNAs, which may be associated with diseases [18].
Reports on the associations between SNPs in pre-miRNA or SNPs in binding sites in 3′ UTR of miRNA targeting mRNAs and human diseases have provided new insights into the molecular mechanisms of pathophysiological processes in humans. rs4846049 (G > T) of MTHFR gene is associated with increased risk for coronary heart diseases, and the potentially pathogenetic mechanism may be SNP-modified posttranscriptional gene regulation by miRNA-149 to MTHFR [19]. MiR-196a binding-site SNP regulates RAP1A expression contributing to esophageal squamous cell carcinoma risk and metastasis [20]. A miRNA-binding SNP (1010A/G) located within 3′-UTR of HOXB5 is associated with gene expression and may be a promising prognostic factor for bladder cancer [21]. The A to G base change of rs999885 may provide a protective effect against chronic HBV infection but an increased risk for HCC in HBV persistent carriers by altering the expression of the miR-106b-25 cluster [22]. miR-200b/200c/429binding site polymorphism in the 3′-UTR of AP-2α gene is associated with cisplatin resistance, suggesting that SNP (rs1045385) A > C variation may be a potential prognostic marker for cisplatin treatment [23].
The results of these observations remain controversial and inconclusive. In the present study, we conducted a meta-analysis in order to derive more precise and more comprehensive estimation of the associations between the SNP hsa-mir-499 rs3746444 T > C and susceptibility to cancers to quantify the potential between-study heterogeneity.

Study selection
We performed a publication search in PubMed, EMBASE, ISI Web of Science, The Cochrane Library, ScienceDirect, EBSCO, Ovid, Wiley Online Library, and HighWire databases with the following search terms: (miR-499 OR rs3746444) AND (cancer), by two independent investigators (Chen Chen and Shenglan Yang, last search update: June 29, 2014). Hand searches were also performed to identify additional articles in the reference lists of included articles not retrieved by initial electronic search. Publication date and publication language were not restricted in our search. All studies matching the inclusion criteria were retrieved for further examination and data extraction. All of the investigators have received training in literature search, statistics and evidence-based medicine.

Inclusion and exclusion criteria
Studies included in current meta-analysis had to meet all the following criteria: (1) evaluated the associations between the hsa-mir-499 rs3746444 polymorphism and cancer risk; (2) studied on human beings; (3) diseases were confirmed by histology, imaging or pathology; (4) a case-control design; (5) detailed genotype data were provided for the calculation of odds ratio (OR) and 95% confidence interval (CIs); (6) if serial studies of the same population from the same group were reported, the latest study was included. Studies were excluded when they represented duplicates of previous publications, or were meta -analyses, meeting abstracts, letters, reviews, or editorial articles.

Data extraction
Two investigators (Chen Chen and Shenglan Yang) independently extracted data from the included studies using a standard protocol and data-collection form according to the inclusion criteria listed above, and reached consensus on all items. Data extracted from eligible studies included the first author's name, year of publication, country of origin, ethnicity, cancer type, genotyping method, total numbers of cases and controls, and genotype frequencies of cases and controls. For study including subjects of different countries of origin group from same ethnicity, we combined them together. The ethnic descents were categorized as Caucasian or Asian. If different results were generated, the two authors would check the data and have had a discussion to come to an agreement. Two senior investigators (Yan Wang and Dao Wen Wang) were invited to the discussion if disagreement still existed.

Statistical analysis
For each study, the departure of frequencies of hsa-mir-499 polymorphism from expectation under Hardy-Weinberg equilibrium (HWE) was assessed by the goodness-of-fit chi-square test in controls. P < 0.05 was considered representative of a departure from HWE.
OR corresponding to 95% CI was used to assess the strength of association between hsa-miR-499 rs3746444 T > C polymorphism and susceptibility to cancer risk. The significance of the pooled OR was determined by the Z-test, and P < 0.05 was considered as statistically significant. Pooled ORs were calculated for allele frequency comparison (C versus T, TC versus TT, CC versus TT, TC/CC versus TT (dominant) and CC versus TC/TT (recessive)), respectively. Subgroup analyses were done by racial descent and cancer type. Statistical heterogeneity among the studies was estimated using chi-square-based Q-test, a P value greater than 0.1 indicates no significant heterogeneity and the pooled OR was estimated by the fixed-effects model (the Mantel-Haenszel method); otherwise, the random-effects model (the DerSimonian and Laird method) was employed [46]. Publication bias of literatures was evaluated by funnel plots and the linear regression asymmetry test by Egger et al. [47]. An asymmetric plot suggests a possible publication bias and the P value of Egger's test less than 0.05 was considered representative of statistically significant publication bias. Sensitivity analysis was conducted by deleting one study at a time to examine the influence of individual data set on the pooled ORs.
All of the statistical tests were performed with Review Manage (v.4.2; Oxford, England) and STATA software version 11.0 (STATA Corporation, College Station, TX, USA). All the P values are two-sided.

Meta-analysis results
The association between hsa-mir-499 rs3746444 polymorphism and susceptibility to cancer was analyzed in 31 independent studies. The results are shown in Table 2 and Additional file 2: Figures S1, S2 and S3. Significant association between rs3746444 polymorphism and susceptibility to cancer was identified in TC versus TT and TC/CC versus TT (dominant) models, when all the eligible studies were pooled (C versus T: OR = 1.08, 95% CI 0.99-1.17, P = 0.08; TC versus TT: OR = 1.11, 95% CI 1.00-1. 23

Publication bias
Begg's funnel plot and Egger's test were performed to assess the publication bias of included studies. Symmetrical funnel plots were obtained in all comparison models ( Figure 2). Further, Egger's test confirmed the absence of publication bias in all genetic models (P > 0.05).

Sensitivity analysis
Exclusion of one study at a time was performed to examine the influence of the individual data set to the pooled ORs. Analyses were performed using STATA by means of the metainf program. Results showed that the pooled OR estimates were similar with those of the excluded studies (Figure 3).

Discussion
In the present study, hsa-mir-499 rs3746444 polymorphism and susceptibility to cancer was evaluated in different genetic models in 12799 cancer cases and 14507 controls. Significant association between rs3746444 polymorphism and susceptibility to cancer was identified in TC versus TT and TC/CC versus TT (dominant) models. Subgroup analysis absed on ethnicity revealed significantly increased     Many studies investigated the role of SNPs presented in precursor and mature miRNAs, and their influences on susceptibility and progression of various diseases. miRNA-associated SNPs can have direct or indirect effects: 1) direct effects are produced by the SNPs in the pri-miRNA, pre-miRNA or mature miRNA that impair or enhance miRNA processing or function; 2) indirect effects are derived from the SNPs in miRNA promoters that affect transcription and SNPs in an mRNA that create or destroy a target site [71].
Genome-wide patterns of human polymorphisms in miRNAs and miRNA target sites in 3′ UTRs of mRNAs revealed that only~10% of human pre-miRNAs have documented SNPs, and <1% of mature miRNAs have SNPs in the functional seed region [72]. hsa-mir-499 rs3746444 polymorphism is located in the seed region of the mature miR-499 sequence. This T > C polymorphism-resulted mismatch may affect target mRNA expression. Several case-control studies have investigated the association between rs3746444 polymorphism and risks of various cancers. However, the results were controversial.  Different from past studies, we analyzed the specific association among different ethnic groups throughout Asia. In this meta-analysis, we first found that rs3746444 polymorphism was associated with increased cancer risk mainly in Asian population, especially in Chinese and Iranians. It is important to notice that, although rs3746444 polymorphism contributes to cancer risk in Chinese and Iranians, the genetic models are different. The genetic variations between Chinese Han population and Arabian population may contribute to the differences. Interestingly, we also found that rs3746444 polymorphism showed a decreased risk in esophageal cancer. This suggests that potentially different mechanisms may underlie tumorgenesis in different pathological backgrounds and the environment variations.
It has been suggested that miRNAs may influence gene expression of~30% of protein-coding genes by binding to the mRNA incompletely, and a miRNA may affect multiple target mRNAs involved in same pathophysiological progress [3,73]. This means that rs3746444 polymorphism may contribute to single gene function in certain disease model at molecular level. When taken together in epidemiological studies of populations, various effects of rs3746444 polymorphism on different genes may result in different associations with diseases at phenotype level.
Our meta-analysis pooled the largest numbers of cases and controls from included studies, which significantly increases the statistical power but some limitations should be considered. Firstly, lack of the consideration of combined genetic factors together with environmental exposures, while a more precise analysis needs to be conducted if more detailed data are available. Secondly, although sensitivity analysis showed that each study can not affect the overall results, some heterogeneity was evident in some of the comparisons. This may be due to different ethnicities and different pathology. And expanding the sample size may yield a better result. Thirdly, some unpublished data (such as negative report, not written in English and lack of studies in African population) may potentially influence the results of our meta-analysis. At last, different genotyping strategies may contribute to the bias in the analysis.
Recently several meta-analyses systematically reviewed the potential association of rs3746444 polymorphism with susceptibility to diseases [40][41][42][43][44][45][74][75][76][77]. However, the results were controversial across these meta-analyses. Such discrepancies in the results of these systematic reviews may due to different search periods, databases searched, inclusion and exclusion criteria. Our efforts to search the maximum range of databases yielded the most comprehensive studies in the similar time period especially the most recent publications [53,59]. Further, diagnosis of cancer is largely relied on confirmation by histology, imaging or pathology. Some of the other studies omitted this fundamental condition. Most importantly, deviations from Hardy-Weinberg equilibrium can inflate the chance of a false-positive association, and sensitivity analysis should be performed, pooling with and without studies not in HWE, to test the robustness of the results [78,79]. Only few of other studies performed sensitivity analysis, which may limit the reliability of their conclusions. Although some studies claimed that their studies were in HWE, a careful examination of the studies revealed that the data did not reach HWE. In conclusion, in our study, relevant literatures selected from broad databases with stringent standards would be expected to provide a more reliable conclusion.

Conclusions
In summary, though with limitations, our meta-analysis suggested that hsa-mir-499 rs3746444 T > C polymorphism is associated with the risk of cancer, especially in Asians, mainly in Chinese and Iranians. However, rs3746444 T > C polymorphism is negatively associated with the risk of esophageal cancer. Larger studies from different ethnic groups and studying different types of cancers with detailed information are needed to further clarify the association between rs3746444 polymorphism and cancer risk.