MicroRNAs (miRNAs) are single-stranded RNAs of 20–30 nt that mediate post-transcriptional gene silencing. Over the past decade, it has been reported that miRNAs play key roles in the initiation and progression of human cancer. For instance, the let-7 family of miRNAs can target oncogenes, such as MYC and RAS family members, to inhibit tumor growth (Johnson et al., 2005; Kim et al., 2009). Osada & Takahashi (2011) reported that the miR-17–92 cluster interacts with E2F1 and MYC to promote tumor development in lung cancer. Although individual miRNAs can have either tumorigenic or tumor-suppressive functions, Lu et al. (2005) analyzed 217 mammalian miRNAs from 334 samples and found that miRNAs were globally down-regulated in cancers. Thomson et al. (2006) indicated that this widespread down-regulation of miRNAs is due to a failure in a Drosha processing step, suggesting that miRNA biogenesis may be impaired in cancer.
The majority of the miRNA biogenesis pathway is shown in Fig. 1. In the nucleus, microRNA genes are transcribed as primary miRNAs (pri-miRNAs) by RNA polymerase II (Pol II). These long pri-miRNAs are then cleaved by the double-stranded RNaseIII enzyme Drosha and its essential cofactor DGCR8 (Denli et al., 2004; Lee et al., 2003), generating pre-miRNAs (∼65 nt). Pre-miRNAs are then exported from the nucleus to the cytoplasm by Exportin-5 (XOP5) and its cofactor Ran (Yi et al., 2003). Following nuclear export, the pre-miRNA loop is cleaved by the RNaseIII enzyme Dicer to generate miRNA duplexes (Park et al., 2011). Dicer, TRBP or PACT, and Argonaute proteins Ago1–4 regulate the processing of pre-miRNA and the assembly of the RNA-induced silencing complex (RISC) (Chendrimada et al., 2005). Together with members of the GW182 family of proteins, one strand of the RNA duplex remains anchored to the Ago protein as the mature miRNA while the other one degraded (Gregory et al., 2005).
As mentioned in Fig. 1, the Drosha, DGCR8, Exportin-5, Ran, and Dicer proteins are crucial components in the miRNA biogenesis pathway. Dysregulation of these miRNA-processing machinery components may disturb the product, leading to disease thus this study was focused on this five core components in miRNA biogenesis pathway and explore its link with cancer risk. Indeed, the expression levels of Drosha and Dicer are down-regulated in ovarian cancer (Kobel, Gilks & Huntsman, 2009) and neuroblastomas (Lin et al., 2010), while DGCR8 is up-regulated in esophageal cancer (Sugito et al., 2006), bladder cancer (Catto et al., 2009), and prostate cancer. Exportin-5 is also down-regulated in bladder cancer (Catto et al., 2009). There are several studies implicating single-nucleotide polymorphisms (SNPs) present in those five miRNA biogenesis pathway genes were link with cancer risk, e.g., rs10719 T>C in DROSHA; rs417309 G>A and rs1640299 T>G in DGCR8; rs11077 in XPO5; rs14035 C>T and rs3803012 A>G in RAN; and rs1057035 T>C, rs3742330 A>G, and rs13078 T>A in DICER1. However, the conclusions of previous studies remain inconsistent. A comprehensive review and meta-analysis of SNPs in miRNA biogenesis pathway components is needed. In this study, a meta-analysis was conducted to evaluate the association between SNPs in five genes (DROSHA, DGCR8, XPO5, RAN, and DICER1) involved in the canonical microRNA biogenesis pathway and human cancer risk.
Materials and Methods
Identification of eligible studies
Two reviewers (He and Zhu) searched the electronic databases PubMed, Google Scholar, CBM, and CNKI for studies containing at least one of the five gene names and its aliases (“DROSHA ETOHI2 HSA242976 RANSE3L RN3 RNASE3L RNASEN,” “DGCR8 C22orf12 DGCRK6 Gy1 pahsa,” “XPO5 exp5,” “RAN ARA24 Gsp1 TC4,” or “DICER DCR1 Dicer1e HERNA K12H4 8-LIKE MNG1 RMSE2”). The following terms were also used in each search: cancer, carcinoma, neoplasm, tumor, gene name and its aliases, gene polymorphism, allele, and variation. The search was limited to articles published in English or Chinese through Jun 21, 2016.The reference lists of pertinent articles were checked as well. Articles were included in the study if (1) they reported a correlation between a SNP from DROSHA(rs10719), DGCR8(rs417309 or rs1640299), XPO5(rs11077), RAN(rs14035), or DICER1(rs1057035, rs3742330, or rs13078) and cancer risk; (2) they presented a case-control study; and (3) genotype frequencies were available. Articles were excluded if (1) they were review articles or focused on animals; (2) they lacked data necessary for a systematic review; or (3) the control subjects exhibited a departure from Hardy–Weinberg equilibrium (HWE). In the case where there was more than one publication using the same samples, the article with the most recent publication date was selected. The two reviewer reached consensus on each study.
The following information were extracted by two reviewers (He and Zhao) independently and disagreement were solved by discussion: first author’s name, publication year, country of origin, ethnicity, cancer types, genotype frequencies in case and control groups, HWE results, and source of control samples. Quality assessment of articles was conducted using the Newcastle-Ottawa Quality Assessment Scale (NOS) (Stang, 2010). If more than one type of cancer or multistage research was involved in a single article, data were extracted as independent articles.
Hardy–Weinberg equilibrium (HWE) was evaluated in the controls of each SNP using the Chi-square test. Odds ratios (OR) with 95% confidence intervals (CI) were calculated by the Z test to evaluate the association between the SNPs and cancer susceptibility under a dominant model, P < 0.05 was considered statistically significant. The heterogeneity between articles was examined using a Q-test and I2 index (Higgins & Thompson, 2002). Random-effects models (Dersimonian & Nan, 1986) were used if heterogeneity between articles was reported (P < 0.10, I2 > 50%), otherwise fixed-effects models (Mantel & Haenszel, 1959) were applied. Sensitivity analyses were performed with a leave-one-out method by removing each article one at a time and repeating the analysis. Publication bias was evaluated by Egger’s test and Begg’s test, with a P-value > 0.10 considered evidence for no potential publication bias (Begg & Mazumdar, 1994; Egger et al., 1997). All P value were two sided and all analyses were performed by Stata statistical software (Version 12.0; StataCorp, College Station, Texas USA).
A total of 832 relevant articles were identified based on our search strategy. Article selection procedures for all SNPs from DROSHA, DGCR8, XPO5, RAN, and DICER1 are shown in a flow chart (Fig. 2). As a result, 22 case-control articles including nine SNPs (rs10719, rs417309, rs1640299, rs11077, rs14035, rs3803012, rs1057035, rs3742330, and rs13078) met the inclusion criteria for this study. Characteristics of the eligible studies are summarized in Table 1. Among the applicable articles, 15 articles were about studies of subjects who were of Asian descent, and seven articles were about studies of subjects who were of Caucasian descent. The analyzed articles for each SNP are shown in Table S1.
|Author||Year||Cancer type||Country||Ethnicity||Controls||Case||Control||Method||Polymorphism site||nos|
|Horikawa et al.||2008||renal cell carcinoma||American||Caucasian||PB||279||278||SNPlex||rs417309, rs13078, rs3742330, rs10719, rs14035||8|
|Yang et al.||2008||bladder cancer||American||Caucasian||HB||746||746||SNPlex||rs417309, rs1640299, rs13078, rs3742330, rs11077||7|
|Ye et al.||2008||esophageal cancer||American||Caucasian||HB||346||346||SNPlex||rs417309, rs1640299, rs13078, rs3742330, rs14035, rs11077||7|
|Kim et al.||2010||lung cancer||Korea||Asian||HB||100||100||rs417309, rs1640299, rs13078, rs3742330, rs10719, rs14035||7|
|Sung et al.||2011||breast cancer||Korea||Asian||HB||559||567||TaqMan||rs1057035, rs11077||7|
|Ma et al.||2012||head and neck cancer||China||Asian||HB||397||900||TaqMan||rs1057035, rs3803012||7|
|Chen et al.||2013||cervical carcinoma||China||Asian||HB||1486||1549||TaqMan||rs1057035, rs3803012||7|
|Jiang et al.||2013||breast cancer||China||Asian||HB||878||900||TaqMan||rs417309, rs1640299, rs1057035, rs13078, rs10719, rs3803012||7|
|Liu et al.||2013||hepatocellular carcinoma||China||Asian||HB||1300||2688||TaqMan||rs417309, rs164029, rs1057035, rs3803012||7|
|Slaby et al.||2013||colorectal cancer||Czech||Caucasian||HB||197||202||TLDA||rs1057035||7|
|Yuan et al.||2013||bladder cancer||China||Asian||HB||685||730||TaqMan||rs1057035, rs13078, rs3742330, rs10719||7|
|Roy et al.||2014||oral cancer||India||Asian||HB||451||452||Taqman||rs14035||7|
|Cho et al.||2015||Colorectal Cancer||Korean||Asian||HB||408||400||PCR-RFLP||rs3742330, rs10719, rs14035, rs11077||7|
|Martin-Guerrero et al.||2015||Lymphocytic Leukemia||Spanish||Caucasian||HB||123||391||Taqman||rs417309, rs1640299, rs1057035,rs13078, rs10719, rs14035||7|
|Xie et al. (2015)||2015||gastric cancer||China||Asian||HB||137||142||rs3742330, rs14035, rs11077||7|
|Zhao et al. (2015)||2015||colorectal cancer||China||Asian||HB||163||142||rs3742330, rs14035, rs11077||7|
|Zu et al. (2013)||2013||lung cancer||China||Asian||HB||600||600||TaqMan||rs1057035||7|
|Zheng Liang (2013)||2013||esophageal cancer||China||Asian||HB||380||380||rs3742330||7|
|Zhang (2012)||2012||gastric cancer||China||Asian||HB||1674||1852||TaqMan||rs3803012, rs1057035||7|
|Buas et al. (2015)||2015||esophageal cancer||Europe||Caucasian||HB||600||600||TaqMan||rs14035, rs11077||7|
|Gutierrez-Camino et al. (2014)||2014||Lymphocytic Leukemia||Spanish||Caucasian||HB||213||387||rs1640299||7|
|Gene(locus)||Position||studies||Method||Cases /controls||WW vs.WM+MMa||Pb||I2||WW vs. WM+MMa||Pc||I2||WW vs.WM+MMa||Pd||I2|
|ORb(95% CI)||ORc(95% CI)||ORd(95% CI)|
A significant association was observed between cancer risk and the DGCR8 rs417309 G>A polymorphism in the overall pooled analysis (GA + AA vs. GG: OR 1.22, 95% CI [1.04–1.42]) and in the Asian samples (GA + AA vs. GG: OR 1.44, 95% CI [1.13–1.83]). For the rs1057035 T>C variant in DICER1, a significant (13%) increase in cancer risk was found in the overall pooled analysis (AA vs. AC + CC: OR 1.13, 95% CI [1.05–1.22]). Subgroup analysis determined by sample ethnicity revealed that genotypes containing the major allele increased the risk of cancer in those of Asian descent (AA vs. AC + CC: OR 1.13, 95% CI [1.04–1.22]).
Sensitivity analysis was conducted to assess the robustness of the results. Each individual article was omitted to measure its effect on the pooled ORs (Fig. S3). The sensitivity analysis forest plot showed that no individual article dramatically affected the pooled OR for any SNP.
Funnel plots did not show obvious asymmetry for any locus (Fig. 4). The results of Egger’s test and Begg’s test revealed that rs10719 in DROSHA presented publication bias (Begg’s test: Z = 1.47, P = 0.14; Egger’s test: P = 0.09, 95% CI [−0.87–7.07]), suggesting that the results for rs10719 should be treated with caution. For the other SNPs, neither test indicated potential publication bias (Table 3).
In this study, we searched for articles indicating an association between SNPs in miRNA biogenesis pathway genes and human cancer risk. A total of 21 articles were identified, and 9 SNPs were evaluated for association with cancer susceptibility. The results demonstrated that the GG genotype of rs417309 in DGCR8 was significantly rarer among cases compared with controls in the overall pooled analysis and the TT genotype of rs1057035 in DICER1 was associated with a 13% increase in cancer risk. In addition significant associations were found in the Asian population but not in the Caucasian population, which suggested a possible ethnic difference in the genetic background and the environment. Significant association were not observed in DROSHA (rs10719 T>C), XPO5 (rs11077 A>G) and RAN (rs14035 C>T, rs3803012A>G). Genome-wide association studies (GWAS) involving associations between esophageal adenocarcinoma risk and RAN(rs14035 C>T) or XOP5(rs11077 A>G) were included in this meta-analysis (Buas et al., 2015). The results of our meta-analysis were consistent with the findings of these GWAS publications (OR 1.10, 95% CI [0.76–1.59] compared with OR 1.17, 95% CI [0.87–1.57] in our meta-analysis for RAN rs14035 C>T and OR 0.94, 95% CI [0.85–1.04] compared with OR 0.93, 95% CI [0.87–1.01] in our meta-analysis for XOP5 rs11077 A>G).
|Category||Studies||Begg’s test||Egger’s test|
DROSHA and its essential cofactor DGCR8 are critical protein that executes the initial step in microRNA processing (Lee et al., 2003). Dysregulation of DROSHA and DGCR8 has been observed in many cancers such as epithelial skin cancer (Sand et al., 2010), breast cancer (Chen et al., 2013) and ovarian cancer (Guo et al., 2015). Pre-miRNAs are exported into the cytoplasm to generate mature miRNAs through XOP5 and it cofactor RAN (Lin & Gregory, 2015), some heterozygous XPO5 variations were found in colon, gastric and endometrial cancer (Melo et al., 2010). The variation of XPO5 impairs pre-miRNA exportation and lead to defect in miRNAs biogenesis. Genetic and epigenetic association studies reported that the genetic variation of XPO5 was associated with the risk of breast cancer (Leaderer et al., 2011). DICER1 is another enzymes that plays critical role in the cleavage of pre-miRNAs into their mature form (Lin & Gregory, 2015). The variation of DICER1 were correlated with cancer risk through affecting cell proliferation and apoptosis (Dedes et al., 2011) all those component are key enzymes in miRNAs mature process. Since the global impairment of mature miRNAs is emerging as a common feature of human tumors (Melo et al., 2010), and given the critical functions of Drosha, DGCR8, Exportin-5, Ran, and Dicer in miRNA biogenesis, it is logical to presume that genetic polymorphisms in these genes may affect the processing of miRNAs and, thus, cancer risk. As a matter of fact, several studies have shown significantly association about single nucleotide polymorphism of MicroRNA biogenesis pathway gene with different cancer risk. Cho et al. (2015) reported that RAN rs14035 CT heterozygotes and XPO5 rs11077 AA carriers experienced reduced risk of colorectal cancer in Korean population. Leaderer et al. (2011) performed both genetic and epigenetic association studies of XPO5 in a case-control study of breast cancer and found XPO5 rs11544382 was associated with cancer risk (the analysis of XPO5 rs11544382 with cancer risk was not performed due to lack of availed data). Chen et al. (2013) analyzed 1,486 cervical cases and 1,549 cancer-free controls in Chinese population and found that a single nucleotide polymorphisms in DICER rs1057035 and RAN rs3803012 was associated with cervical cancer risk. A case-control and further research conducted by Jiang et al. (2013) suggested that DGCR8 rs417309 G>A might affect breast cancer risk through interrupting the miRNA binding. However, conclusions of miRNA biogenesis pathway genes polymorphisms and cancer risk remain inconsistent, which may cause by the heterogeneity of the cancer subtype, sample size and the ethnicity of patients. Yu, Kuang & Yin (2015) conducted a meta-analysis synthesizes sevens studies with multi-type cancer and found the C allele of the DICER rs1057035 polymorphism probably decreases cancer risk (We included more studies in our studies and found the same results in DICER rs1057035). This study based on the association of a single SNP or gene with cancer risk, which may weaken gene’s biological value. All five miRNA biogenesis pathway were evaluated in our analysis, which may provide a more comprehensive view of the true system. To the best of our knowledge, this study represents the first meta-analysis to focus on the association between human cancer risk and SNPs across the entire miRNAs biogenesis pathway.
Sensitivity analysis indicated that the results of our study were robust. Egger’s test and Begg’s test revealed that publication bias has little influence on the results, with the exception of rs10719 in DROSHA, suggesting that the results for rs10719 should be considered with caution. Despite these results, some limitations still apply to this meta-analysis. First, all results were based on unadjusted estimates because few adjusted ORs were reported in the included articles. Biases caused by confounding variables such as sex, age, and smoking may therefore persist. Secondly, heterogeneity among different cancers may cause real effects to be hidden when pooling all cancer types. A stratified analysis based on cancer type was not performed due to the lack of an adequate number of publications for a single cancer type. Third, although the results of the publication bias analysis indicated no bias for any locus except rs10719, it is possible that there actually was publication bias but that it was not identified by the Egger’s and Begg’s tests due to the low power given to the small number of studies (Begg & Mazumdar, 1994; Egger et al., 1997). Finally, we applied a traditional fixed model and a random model to pooled multiple-ethnicity samples; neither of these approaches is ideal for capturing the heterogeneity of effects observed across different ethnic populations (Li & Keating, 2014).
In summary, this meta-analysis suggests a potential role of the miRNA biogenesis genes DGCR8 (rs417309 G>A) and DICER1 (rs1057035 T>C) in cancer risk. The identified markers can thus potentially be used as biomarkers for cancer diagnosis. However, future well-designed studies, including those focused on a single cancer type with larger sample sizes, are still needed to confirm these observed associations due to the limitations mentioned above.