Common genetic polymorphisms in pre-microRNAs and risk of bladder cancer
World Journal of Surgical Oncologyvolume 13, Article number: 297 (2015)
At present, inconsistent association between single nucleotide polymorphism (SNP) in pre-miRNAs (hsa-mir-196a2 rs11614913 C/T, hsa-mir-499 rs3746444 A/G, and hsa-mir-146a rs2910164 C/G) and bladder cancer were obtained in limited studies. We performed a case–control study to test whether these three common polymorphisms are associated with bladder cancer. One hundred fifty-nine patients affected by bladder cancer and 298 unrelated healthy subjects were enrolled in the study.
Using polymerase chain reaction–restriction fragment length polymorphism assay (PCR–RFLP), genotypes of these three SNPs were determined, and their associations with bladder cancer, as well as with clinic pathological factors, and tumor progression were analyzed.
No association between bladder cancer risk and variant allele of hsa-mir-196a2 rs11614913 C/T, hsa-mir-499 rs3746444 A/G, or hsa-mir-146a rs2910164 C/G was observed. Heterozygous genotype (CT genotype) of rs11614913 was associated with a significantly decreased bladder cancer risk (P = 0.004, OR = 0.56, 95 % CI = 0.38–0.83). Further stratified analyses showed that rs2910164 is associated with the tumor stage in a recessive model and with metastasis in a dominant model (P = 0.012, OR = 0.20, 95 % CI = 0.05–0.72 and P = 0.04, OR = 2.63, 95 % CI = 1.03–6.67, respectively). No association between hsa-mir-499 rs3746444 A/G and bladder cancer was observed.
Our results suggested hsa-mir-196a2 rs11614913 C/T is associated with a significantly decreased risk of bladder cancer and hsa-mir-146a rs2910164 GG genotype is associated with clinical stage and metastasis in bladder cancer.
As the most common malignant tumor of the urinary tract, it is estimated that there were 74,690 new cases and 56,390 patients died from bladder cancer in United States in 2014 . Tobacco smoking and chemical carcinogens are the two well-established external environmental factors of bladder cancer . Besides, accumulating evidence indicates that inappropriate diet and lifestyle may have a relationship with bladder cancer . These factors could be an explanation for the rising morbidity of bladder cancer all over the world. It is worth noting that the majority of exposed individuals do not suffer from bladder cancer in their lifetime. This phenomenon suggested that other causes, such as the genetic factor may also play a role in bladder carcinogenesis.
MicroRNAs (miRNAs) are a class of ancient, non-coding, and single-stranded molecules that participate in transcriptional and translational regulation by pairing with target mRNAs . Although represents only 1–4 % in the human genome, a single miRNA can regulate approximately 200 target mRNAs . This means that miRNAs may play critical roles in a diverse range of biological processes, such as cell proliferation, tissue differentiation, embryonic development, and apoptosis [5, 6]. Mutation of miRNAs and dysregulation of miRNAs and their targets genes have been observed in various diseases and several reports have demonstrated that miRNA plays an important role during carcinogenesis [7, 8].
Single nucleotide polymorphism (SNP) is the most common type of genetic variation that usually appears in miRNA coding genes, miRNA biogenesis genes, or in the miRNA-binding sites of mRNAs. The variations may alter the expression and function of miRNA, which may relate to the tumorigenesis . Recently, two studies analyzed the relationship between common SNPs in pre-miRNAs and urinary bladder cancer . While one obtained the negative results, the other indicates that miR-146a rs2910164 plays an important role in bladder cancer. In addition, they have associated the C allele with a significantly decreased risk and recurrence of bladder cancer . To clarify the inconclusive results of previous studies, we explored the association between these three SNPs (rs2910164, rs11614913, and rs3746444) in pre-miRNAs (hsa-mir-146a, hsa-mir-196a2, and hsa-mir-499, respectively) and the risk of bladder cancer in a population of southwest China.
Patients and controls
All participants were form the Han population living in Sichuan Province of China and examined at the Department of urology, West China Hospital of Sichuan University. This study was approved by the Investigation and Ethics Committee of the West China Hospital of Sichuan University. Details were told to all participants and written informed consent was signed after the interview. The present study enrolled 159 unrelated patients with bladder cancer (124 men and 35 women, mean age: 64.09 ± 12.31) from the West China Hospital of Sichuan University between 2009 and 2012. The clinical diagnosis of bladder cancer was based on histopathologically confirmed specimens. Two hundred ninety-eight unrelated healthy individuals (140 men and 158 women, mean age: 60.21 ± 11.32) from a routine health survey were enrolled as controls. Patients who previously had cancer, radiotherapy and chemotherapy, or metastasized cancer from other or unknown origins were excluded. Control subjects were genetically unrelated individuals, and those with history of bladder cancer or other type of illness were excluded in this study.
Three SNPs, including hsa-mir-196a2 rs11614913 C/T, hsa-mir-499 rs3746444 A/G, and hsa-mir-146a rs2910164 C/G were genotyped by the PCR–RFLP method. The instruction manual details that 300 μl of EDTA-anticoagulated peripheral blood samples should be collected from the patients and that their genomic DNA be isolated by a DNA isolation kit from Bioteke (Peking, China).
PIRA PCR designer (http://primer1.soton.ac.uk/primer2.html) was used to establish primers. The primers are used for amplification of these three SNPs and their specific restriction enzyme was shown in Table 1. The PCR reactions mixture was consisted of 2.0 μl 10× PCR buffer, 1.5 mmol/L MgCl2, 0.2 mmol/L dNTPs, 0.4 μmol/L each primer, 100 ng of genomic DNA and 1 U of Taq DNA polymerase. Briefly, cycling parameters were summarized as follows: 94 °C for 2 min; 30 cycles at 94 °C for 30 s; 30 s with different temperature parameters: 62 °C for rs11614913, 66 °C for rs3746444 and 58 °C for rs2910164; 45 s at 72 °C and the last extension procedure at 72 °C for 7 min. PCR–RFLP assay was used to obtained the genotypes of the three SNPs: rs11614913 C/T, rs3746444 A/G, rs2910164 C/G .
Two researchers, who were irrelevant for this study, viewed the results independently. If there is no consensus, repeated assays would be required. Besides, 10 % duplicated samples were randomly selected to perform for the three SNPs, with a concordance rate of 100 %.
Genotypic frequencies of the observed SNPs were available by directly counting and the χ 2 test was used to evaluate the Hardy-Weinberg equilibrium. Data were carried out by using computer software SPSS (SPSS Inc., Chicago, IL, USA) and odds ratio (OR) and respective 95 % confidence intervals (95 % CI) were reported to compare the differences between the two groups. Genotypic association tests were conducted by using SNPstats, including codominant, dominant, recessive, overdominant, and log-additive genetic models .
Results and discussion
Genotyping for all the three SNPs were performed in 159 bladder cancer patients and 298 control subjects. As expected, the results showed genotype frequencies for these three SNPs in control subjects were consistent with those from the Hardy-Weinberg equilibrium. Allele frequencies of rs11614913 C/T, rs3746444 A/G, and rs2910164 C/G in bladder cancer patients and controls are summarized in Table 2. As shown in Table 2, no significant difference for the allele frequencies of these three SNPs was observed between bladder cancer patients and control subjects (P = 0.976 for rs11614913, P = 0.348 for rs3746444, and P = 0.417 for rs2910164). Results of genotypic association tests performed by SNPstats are shown in Table 3. As shown in Table 3, significant association was identified for rs11614913 in an overdominant model. Compared with the combined homozygous genotypes (TT/CC genotypes), a significant decrease in bladder cancer risk was associated with the heterozygous genotype (CT genotype) of rs11614913 (P = 0.004, OR = 0.56, 95 % CI = 0.38–0.83). No association between the rs3746444 and rs2910164 genotypes and bladder cancer risk was observed in codominant, dominant, recessive, overdominant, or log-additive genetic model. These results suggest that the mir-196a2 rs11614913 C/T, but not the mir-499 rs3746444 A/G or the mir146a2 rs2910164 C/G, was associated with bladder cancer in the Chinese population studied. This indicates that polymorphisms in pre-microRNAs may play a role in the pathogenesis of bladder cancer.
Further analysis stratified bladder cancer patients’ characteristics including age, gender, smoking status, tumor stage, tumor grade, recurrence, and metastasis. As shown in Table 4, rs2910164 was found to be associated with the tumor stage in a recessive model (P = 0.012, OR = 0.20, 95 % CI = 0.05–0.72). Table 5 shows that rs2910164 is associated with metastasis in a dominant model (P = 0.04, OR = 2.63, 95 % CI = 1.03–6.67). No association was observed between these three SNPs and bladder cancer patients’ age, gender, smoking status, tumor grade, or recurrence (data not shown).
In the present study, the association between SNPs in pre-miRNAs and bladder cancer, as well as with patients’ characteristics, such as age, gender, tumor grade, tumor stage, smoking status, recurrence, and metastasis were analyzed. Our results revealed that hsa-mir-196a2 rs11614913 C/T is associated with bladder cancer in an overdominant model. Similarly, the hsa-mir-146a rs2910164 is associated with the tumor stage in recessive model and with metastasis in dominant model, respectively. Hsa-mir-499 rs3746444, however, is not associated with bladder cancer or any of the patients’ characteristics.
The miRNAs were composed of 21 to 24 nucleotides which can provide sequence specificity to RNA silencing pathway by pairing with target mRNAs . Myriad evidence has indicated that miRNA is involved in diversity of biological process and many human genes are under its regulation, as well as some well-known tumor suppressor genes [5, 15]. Nowadays, the roles of altered miRNA expression in various diseases have been well established, including malignant cancer . Gottardo and his colleagues analyzed the expression profile of miRNAs in bladder cancer and found that some miRNAs are significantly deregulated in bladder cancer compared with normal bladder mucosa, suggesting the involvement of these genes in the development and progression of bladder cancer . Additionally, several additional miRNA-related SNPs have been demonstrated to be associated with bladder cancer risk, this means that miRNA-related SNPs are involved in mechanism of bladder carcinogenesis .
Emerging studies have reported the association between the hsa-miR-196a2 rs11614913 polymorphism and susceptibility to several types of cancers, including hepatocellular carcinoma, gastric cancer, and breast cancer [18–20]. MiR-196a is one of the mature products encoded by miR196a2 gene, and miR-196a* which contains the SNP rs11614913 is the other mature product . Hoffman et al. showed relative to empty vector control, the expression of mature miR-196a is increased in breast cancer cells both transfected with pre-miR-196a-C or miR-196a-T. Moreover, levels of the mature miR-196a were more than double in cells transfected with pre-miR-196a-C than those with miR-196a-T . This finding suggested that SNP rs11614913 can affect the processing of the pre-miRNA into its mature, functional form, and this polymorphism might be a risk factor of tumorigenesis. This present study’s findings reveal that rs11614913 is associated with bladder cancer in an overdominant model, and the heterozygous genotype (CT genotype) was associated with decrease bladder cancer risk.
The rs2910164 is a common G/C polymorphism which resided at position +60 relative to the first nucleotide of pre-miR-146a . Compared with the G allele, the Gibbs free was lower for rs2910164 C allele, suggesting that the secondary structure was less stable for the C allele. Genetic alterations might affect the expression and function of mature miR-146a . By now, the susceptibility for carcinogenesis caused by the rs2910164 has been demonstrated in a series of studies. Jazdzewski et al. reported that the rs2910164 CC genotype is associated with lower risk in papillary thyroid carcinoma . Studies in hepatocellular carcinoma, cervical squamous cell carcinoma, and bladder cancer also illustrate analogous results [11, 22, 12]. In vitro, Wang et al. found miR-146a can decrease bladder tumor cell proliferation and that the expression level of miR-146a is higher in rs2910164 C allele than G allele . All of these results demonstrated that the C allele may act in a protective role and is more likely to exercise an inhibiting effect on the tumorigenesis. But the studies conducted in the breast/ovarian cancer showed an antithetical conclusion . This diversity might be a result of organ-specific characteristics of rs2910164 and reflected the detailed mechanism of SNPs on cancer is not fully understood yet.
Invasive and metastatic nature was the main characteristic of malignance tumor, which was the primary cause of death for cancer patients . Although exact molecular mechanisms remain unclear, several studies indicate that constitutive nuclear factor-κB (NF-κB) activity could promote the invasion and migration of human cancers [24, 25]. Known as the target of miR-146a, interleukin 1 receptor-associated kinase 1 (IRAK-1) plays an important role in activation of NF-κB [26, 27]. It has been reported that tumor cells would cease to be capable of invasion and migration if there was a downregulation of IRAK-1 induced by miR-146a [26, 28, 29]. In addition, cell movement mode was also believed to be a primary event in invasion and metastasis . WASP family verprolin homologous protein 2 (WASF2) is involved in cell movement mode and several studies have confirmed WASF2 could promote tumor progression [30–33]. Recently, Yao found mir-146a can suppress the migration and invasion of gastric cancer cells by reduction of WASF2 expression . Taken together, the above results suggest that mir-146a can repress tumor invasions and metastases in various ways.
In the present study, we found that SNP rs2910164 is associated with the tumor stage in a recessive model, with GG genotype frequency being significantly lower in the group with invasive bladder cancer. Moreover, allele G carriers (CG/GG genotype) are associated with metastasis in a dominant model. As for bladder cancer, 70 % patients belong to superficial cases . Although the comprehensive treatment such as transurethral resection and postoperative intravesical chemotherapy are known as the gold standard nowadays, the progression rate within 5 years is still high . Our study demonstrates that miR-146a rs2910164 is related to tumor stage and metastasis of bladder cancer. These findings reveal more clues about these SNPs of pre-miRNAs in the carcinogenesis of bladder cancer and may provide more information about tumor progression that has not yet been observed during an initial diagnosis.
In conclusion, we found that hsa-mir-196a2 rs11614913 C/T is associated with a significantly decreased risk of bladder cancer in an overdominant model. In addition, hsa-mir-146a rs2910164 GG genotype is associated with clinical stage in a recessive model and CG/GG genotype is associated with metastasis in a dominant model. It is well known that the genetic polymorphism may execute a different effect among different racial or ethnic populations . Hence, future large-scale investigations in different population should be performed to validate our results among other races and ethnicities.
Siegel R, Ma J, Zou Z, Jemal A. Cancer statistics, 2014. CA Cancer J Clin. 2014;64:9–29.
Cohen SM, Shirai T, Steineck G. Epidemiology and etiology of premalignant and malignant urothelial changes. Scand J Urol Nephrol Suppl. 2000;105–15.
Silberstein JL, Parsons JK. Evidence-based principles of bladder cancer and diet. Urology. 2010;75:340–6.
Garzon R, Marcucci G, Croce CM. Targeting microRNAs in cancer: rationale, strategies and challenges. Nat Rev Drug Discov. 2010;9:775–89.
Esquela-Kerscher A, Slack FJ. Oncomirs - microRNAs with a role in cancer. Nat Rev Cancer. 2006;6:259–69.
Ambros V. The functions of animal microRNAs. Nature. 2004;431:350–5.
Kumar MS, Lu J, Mercer KL, Golub TR, Jacks T. Impaired microRNA processing enhances cellular transformation and tumorigenesis. Nat Genet. 2007;39:673–7.
Kwak PB, Iwasaki S, Tomari Y. The microRNA pathway and cancer. Cancer Sci. 2010;101:2309–15.
Ye Y, Wang KK, Gu J, Yang H, Lin J, Ajani JA, et al. Genetic variations in microRNA-related genes are novel susceptibility loci for esophageal cancer risk. Cancer Prev Res (Philadelphia, Pa). 2008;1:460–9.
Mittal RD, Gangwar R, George GP, Mittal T, Kapoor R. Investigative role of pre-microRNAs in bladder cancer patients: a case–control study in North India. DNA Cell Biol. 2011;30:401–6.
Wang M, Chu H, Li P, Yuan L, Fu G, Ma L, et al. Genetic variants in miRNAs predict bladder cancer risk and recurrence. Cancer Res. 2012;72:6173–82.
Zhou B, Wang K, Wang Y, Xi M, Zhang Z, Song Y, et al. Common genetic polymorphisms in pre-microRNAs and risk of cervical squamous cell carcinoma. Mol Carcinog. 2011;50:499–505.
Sole X, Guino E, Valls J, Iniesta R, Moreno V. SNPStats: a web tool for the analysis of association studies. Bioinformatics (Oxford, England). 2006;22:1928–9.
Lagos-Quintana M, Rauhut R, Lendeckel W, Tuschl T. Identification of novel genes coding for small expressed RNAs. Science (New York, NY). 2001;294:853–8.
Shen J, Ambrosone CB, DiCioccio RA, Odunsi K, Lele SB, Zhao H. A functional polymorphism in the miR-146a gene and age of familial breast/ovarian cancer diagnosis. Carcinogenesis. 2008;29:1963–6.
Gottardo F, Liu CG, Ferracin M, Calin GA, Fassan M, Bassi P, et al. Micro-RNA profiling in kidney and bladder cancers. Urol Oncol. 2007;25:387–92.
Yang H, Dinney CP, Ye Y, Zhu Y, Grossman HB, Wu X. Evaluation of genetic variants in microRNA-related genes and risk of bladder cancer. Cancer Res. 2008;68:2530–7.
Qi P, Dou TH, Geng L, Zhou FG, Gu X, Wang H, et al. Association of a variant in MIR 196A2 with susceptibility to hepatocellular carcinoma in male Chinese patients with chronic hepatitis B virus infection. Hum Immunol. 2010;71:621–6.
Peng S, Kuang Z, Sheng C, Zhang Y, Xu H, Cheng Q. Association of microRNA-196a-2 gene polymorphism with gastric cancer risk in a Chinese population. Dig Dis Sci. 2010;55:2288–93.
Hoffman AE, Zheng T, Yi C, Leaderer D, Weidhaas J, Slack F, et al. microRNA miR-196a-2 and breast cancer: a genetic and epigenetic association study and functional analysis. Cancer Res. 2009;69(14):5970–7.
Jazdzewski K, Murray EL, Franssila K, Jarzab B, Schoenberg DR, de la Chapelle A. Common SNP in pre-miR-146a decreases mature miR expression and predisposes to papillary thyroid carcinoma. Proc Natl Acad Sci U S A. 2008;105(20):7269–74.
Xu T, Zhu Y, Wei QK, Yuan Y, Zhou F, Ge YY, et al. A functional polymorphism in the miR-146a gene is associated with the risk for hepatocellular carcinoma. Carcinogenesis. 2008;29(11):2126–31.
Gupta GP, Massague J. Cancer metastasis: building a framework. Cell. 2006;127(4):679–95.
Karin M. Nuclear factor-kappaB in cancer development and progression. Nature. 2006;441(7092):431–6.
Van Waes C. Nuclear factor-kappaB in development, prevention, and therapy of cancer. Clinical Cancer Res. 2007;13(4):1076–82.
Bhaumik D, Scott GK, Schokrpur S, Patil CK, Campisi J, Benz CC. Expression of microRNA-146 suppresses NF-kappaB activity with reduction of metastatic potential in breast cancer cells. Oncogene. 2008;27(42):5643–7.
Liu G, Park YJ, Abraham E. Interleukin-1 receptor-associated kinase (IRAK) -1-mediated NF-kappaB activation requires cytosolic and nuclear activity. FASEB J. 2008;22(7):2285–96.
Chen G, Umelo IA, Lv S, Teugels E, Fostier K, Kronenberger P, et al. miR-146a Inhibits Cell Growth, Cell Migration and Induces Apoptosis in Non-Small Cell Lung Cancer Cells. PLoS ONE. 2013;8(3):e60317.
Li Y, Vandenboom 2nd TG, Wang Z, Kong D, Ali S, Philip PA, et al. miR-146a suppresses invasion of pancreatic cancer cells. Cancer Res. 2010;70(4):1486–95.
Wyckoff JB, Pinner SE, Gschmeissner S, Condeelis JS, Sahai E. ROCK- and myosin-dependent matrix deformation enables protease-independent tumor-cell invasion in vivo. Curr Biol. 2006;16(15):1515–23.
Semba S, Iwaya K, Matsubayashi J, Serizawa H, Kataba H, Hirano T, et al. Coexpression of actin-related protein 2 and Wiskott-Aldrich syndrome family verproline-homologous protein 2 in adenocarcinoma of the lung. Clin Cancer Res. 2006;12(8):2449–54.
Iwaya K, Norio K, Mukai K. Coexpression of Arp2 and WAVE2 predicts poor outcome in invasive breast carcinoma. Modern Pathol. 2007;20(3):339–43.
Iwaya K, Oikawa K, Semba S, Tsuchiya B, Mukai Y, Otsubo T, et al. Correlation between liver metastasis of the colocalization of actin-related protein 2 and 3 complex and WAVE2 in colorectal carcinoma. Cancer Sci. 2007;98(7):992–9.
Yao Q, Cao Z, Tu C, Zhao Y, Liu H, Zhang S. MicroRNA-146a acts as a metastasis suppressor in gastric cancer by targeting WASF2. Cancer Lett. 2013.
van Rhijn BW, Burger M, Lotan Y, Solsona E, Stief CG, Sylvester RJ, et al. Recurrence and progression of disease in non-muscle-invasive bladder cancer: from epidemiology to treatment strategy. Eur Urol. 2009;56(3):430–42.
Kurth KH, Denis L, Bouffioux C, Sylvester R, Debruyne FM, Pavone-Macaluso M, et al. Factors affecting recurrence and progression in superficial bladder tumours. Eur J Cancer (Oxford, England: 1990). 1995;31A(11):1840–6.
Wang J, Bi J, Liu X, Li K, Di J, Wang B. Has-miR-146a polymorphism (rs2910164) and cancer risk: a meta-analysis of 19 case–control studies. Mol Biol Rep. 2012;39(4):4571–9.
This work was supported by the National Natural Science Foundation of China (No. 81202023).
The authors declare that they have no competing interests.
DS gathered patients’ information and drafted the manuscript. WW carried out the genotyping of these SNPs. LX performed the statistical analysis. ZP conceived of the study and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
DS is a postdoctoral fellow of the Department of Urology, West China Hospital, Sichuan University, Chengdu, PR China.
WW is an associate research fellow of the Department of Pathology, West China Second University Hospital, Chengdu, PR China.
LX is a professor of the Department of Urology, West China Hospital, Sichuan University, Chengdu, PR China.
ZP is a professor of the Department of Urology, West China Hospital, Sichuan University, Chengdu, PR China.