Polymorphisms of the pri-miR-34b/c promoter and TP53 codon 72 are associated with risk of colorectal cancer
- Authors:
- Published online on: December 13, 2013 https://doi.org/10.3892/or.2013.2926
- Pages: 995-1002
Abstract
Introduction
Colorectal cancer (CRC) is one of the most common types of cancer (the third and fourth most common cancer in women and men, respectively) in the world and the second leading cause of cancer-related mortality in developed countries. More than 1.2 million cases are diagnosed globally each year, with ~600,000 deaths (1,2). In order to find new diagnostic and therapeutic tools to reduce CRC-related morbidity, it is key to understand the etiology and biology of CRC. Progression to CRC is caused by an accumulation of various genetic and epigenetic alterations, leading to transformation from a normal tissue to a malignant, and potentially metastatic, tumor. CRC develops through two main genetic pathways that are characterized by different forms of genomic instability, the chromosomal instability (CIN, 85%) and microsatellite instability (MSI, 15%) pathways (3–5). Many tumor-suppressor genes and oncogenes have been described, and the discovery of new tumor markers, including those for CRC, continues at a rapid pace.
A new group of biomarkers, microRNAs (miRs), has recently been established. miRs are 20–25-nucleotide non-coding RNAs that negatively regulate gene transcription at the transcriptional or post-transcriptional level by interacting with the 3′ untranslated regions (UTRs) of specific messenger RNAs (mRNAs) and are key modulators in the control of biological processes, such as cell development, differentiation, proliferation and apoptosis (6–9). Estimates based on bioinformatics and microarray analyses suggest that >30% of all genes are regulated by multiple miRs (10). miR-34s form an evolutionarily conserved miR family, with miR-34a, miR-34b and miR-34c occurring in vertebrates. The miR-34a and miR-34b/c loci are regulated directly by interaction with TP53, which induces apoptosis, cell cycle arrest and senescence (11–13). In addition, reduced miR-34a expression is a frequent feature of both pancreatic tumors and neuroblastomas (14,15) and reduced miR-34b/c expression has been observed in non-small cell lung cancer (16). CpG methylation of miR-34b/c has been found in CRC (2,17), oral squamous cell carcinoma (18), and malignant melanoma, where it correlated with metastatic potential (19).
The downregulation of miR-34 family members in cancer suggests that these miRs function as tumor-suppressor genes, suggesting a possible role as prognostic markers (20). Several mechanisms regulating miR expression, including gene amplification, deletion, epigenetic alterations and single-nucleotide substitution, have been implicated, but not demonstrated (20,21). Although single nucleotide polymorphisms (SNPs) in miRs are not considered functionally important, nucleotide variations in primary (pri)- or precursor (pre)-miRs may affect miR processing and modify miR expression (22). Recently, studies reported that a potentially functional SNP, rs4938723, in the promoter region of pri-miR-34b/c may contribute to susceptibility to hepatocellular carcinoma (23), CRC (24), endometrial cancer (25) and survival in breast cancer (26). However, there are few reports on the relationship between SNPs in the miR-34b/c promoter and risk and prognostic significance in CRC patients.
In the present study, we investigated whether the SNPs rs4938723 (T>C) in the promoter region of miR-34b/c and Arg72Pro (G>C) in codon 72 of TP53 are independently or complementarily associated with the risk and clinical outcomes of CRC and whether the combined effect of these SNPs and metabolic risks (diabetes and hypertension) is related to progression of CRC, which is known to be associated with metabolic disease (27,28) in the Korean population.
Materials and methods
Patients and clinical samples
From June 2000 to January 2009, blood samples were collected from 545 patients diagnosed with CRC at CHA Bundang Medical Center of CHA University in South Korea. We retrospectively obtained information concerning the age; gender; underlying conditions [hypertension (HTN), diabetes mellitus (DM), body mass index (BMI), smoking and alcohol consumption]; tumor size, stage and site; time to progression; and time to mortality. We estimated homocysteine and folic acid levels. The American Joint Committee on Cancer: Classification and Stage Groupings, 7th edition was used for tumor assessment. The cancer-free control group consisted of 428 individuals who were randomly selected from participants in a health-screening program to exclude those with a history of cancer and other medical diseases. All study subjects provided written consent and all were ethnic Koreans. The subjects’ recruitment was approved by the Institutional Review Board of CHA Bundang Medical Center.
Genotyping
DNA was extracted from leukocytes using a G-DEX™ II Genomic DNA Extraction kit (iNtRON Biotechnology, Seongnam, Korea), according to the manufacturer’s instructions. The SNPs miR-34b/c rs4938723 and TP53 Arg72Pro rs1042522 were genotyped by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assays. Primer sequences used for amplification of rs4938723 were: (forward) 5′-CCT CTG GGA ACC TTC TTT GAC CAA T-3′ and (reverse) 5′-TGA GAT CAA GGC CAT ACC ATT CAA GA-3′. Primer sequences used for amplification of TP53 Arg72Pro were: (forward) 5′-TTG CCG TCC CAA GCA ATG GAT GA-3′ and (reverse) 5′-TCT GGG AAG GGA CAG AAG ATG AC-3′. For each of the polymorphisms, 30% of the PCR assays were randomly chosen and repeated, followed by DNA sequencing to validate the RFLP findings. Sequencing was performed using an ABI 3730xl DNA Analyzer (Applied Biosystems, Foster City, CA, USA). The concordance of the quality control samples was 100%.
Statistical analysis
The genotypes for each SNP were analyzed as a three-group categorical variable (reference model) and were also grouped according to the dominant and recessive model. Odds ratios (ORs), adjusted odds ratios (AORs) and 95% confidence intervals (CIs) were used to calculate the strength of association. To analyze baseline characteristics, we used a Chi-square test for categorical data when comparing patient and control baseline data. Multivariate analysis was performed to select independent risk factors for CRC among genotypes and clinical variables using logistic regression analysis. The overall survival was compared using the Kaplan-Meier method and potential variables were verified by multivariate analysis using a Cox regression model. All tests were two-tailed, and a P-value <0.05 was deemed to indicate a statistically significant difference. Analyses were performed using GraphPad Prism 4.0 (GraphPad Software, San Diego, CA, USA) and MedCalc version 11.1.1.0 (MedCalc Software, Ostend, Belgium). The distribution of allele frequencies for the rs4938723 (T>C) and TP53 Arg72Pro (G>C) gene polymorphisms were calculated by Chi-square test to determine whether the observed genotype distributions conformed to the expected Hardy-Weinberg equilibrium (29).
Results
Population
Baseline characteristics of patients with CRC and controls are shown in Table I. The mean age was 62 years in CRC patients and 61 years in controls. Of the CRC patients, 302 (55.4%) were male. Approximately 61.5% of CRC cases had HTN and 33.6% had DM; these were significantly higher values than observed among controls (34.1 and 6.1%, respectively) (P<0.001). However, the number of smokers in the control group exceeded that for CRC patients. Tumors occurred most frequently in the rectum and proximal colon. With regard to tumor-node-metastasis (TNM) staging, the majority of patients (80.7%) had stage II–III disease.
Variant genotypes of TP53 codon 72 are related to reduced CRC risk
Genotype frequencies of miR-34b/c rs4938723 and TP53 Arg72Pro rs1042522 in CRC patients and controls are shown in Table II. In multivariate analysis, the variant genotypes of TP53 Arg72Pro GC and GC/CC were associated with a significantly decreased risk of CRC compared with the wild-type GG genotype (AOR = 0.727, 95% CI = 0.550–0.960 for GC; AOR = 0.735, 95% CI = 0.565–0.958 for GC/CC). However, no overall association was observed between miR-34b/c rs4938723 and CRC risk in our study population. The observed genotype frequencies for miR-34b/c rs4938723 (T>C) and TP53 Arg72Pro (G>C) polymorphisms in the cases and controls were all as expected for Hardy-Weinberg equilibrium (P>0.05).
Combined genotype effects of rs4938723 and TP53 (TT/GC and CC/GG) are associated with decreased CRC risk in colon cancer patients
Combined genotype analyses were conducted to evaluate the combined effects of the two polymorphisms on the risk of CRC (Table III). Nine combined genotypes were estimated from the two polymorphisms in CRC patients. In the multivariate analysis, combined genotypes TT/GC and CC/GG were associated with significantly decreased CRC risk when compared with the wild-type TT/GG genotype (AOR = 0.628, 95% CI = 0.422–0.934 for TT/GC; AOR = 0.381, 95% CI = 0.183–0.793 for CC/GG, respectively). This association was observed only in patients with colon, not rectal, cancer and, moreover, was observed only in proximal colon cancer patients (data not shown).
Table IIICombined genotype frequencies of miR-34b/c rs4938723 and TP53 Arg72Pro in CRC patients and controls. |
SNP rs4938723 with DM is associated with increased CRC risk
Table IV shows CRC risk by combined genetic-environmental effects (HTN, DM, homocysteine and folic acid). TP53 Arg72Pro GG and all genotypes of rs4938723 with HTN were significantly associated with increased risk of CRC. All DM patients, in particular, showed strong positive association with CRC, regardless of genotype. The polymorphism rs4938723 with DM was associated with a significantly increased CRC risk compared with wild-type TT and TP53 Arg72Pro CC with DM significantly decreased the risk when compared with wild-type GG. No differences were observed between homocysteine and folic acid levels for any genotype (data not shown). A homocysteine level >11.7 μmol/l was associated with increased risk of CRC in rs4938723 TC and TC/CC, but with decreased risk for TP53 Arg72Pro SNPs. Folic acid level <4.45 ng/ml was associated with an increased CRC risk, but the rs4938723 and TP53 Arg72Pro polymorphisms decreased the risk of CRC when compared with wild-type (Table IV).
Polymorphisms rs4938723 and TP53 Arg72Pro show a trend toward, but are not significantly associated with, improved survival
The 3-year survival rate was estimated in each miR-34b/c rs4938723 and TP53 Arg72Pro patient group (Table V). In the multivariate analysis, overall survival of the variant miR-34b/c rs4938723 and TP53 Arg72Pro genotypes was more evident in subjects carrying polymorphisms than in wild-type. However, we did not find a significant association between miR-34b/c rs4938723 and TP53 Arg72Pro genotypes and survival of CRC (Fig. 1). Stratified analyses (Table VI) showed that the risk reduction effect of the variant TP53 Arg72Pro GC/CC genotypes was significant in subjects without HTN (AOR = 0.635, 95% CI = 0.418–0.966). However, no significant interactions were observed between miR-34b/c rs4938723 TC/CC and other clinical features.
Table VMultivariate survival analysis according to miR-34b/c rs4938723 and TP53 Arg72Pro polymorphisms. |
Discussion
In the present study on CRC in a Korean population, we investigated the correlation between the SNPs rs4938723, in the promoter region of pri-miR-34b/c and TP53 Arg72Pro and the risk of CRC. We found that the TP53 Arg72Pro polymorphism was significantly associated with a decreased risk of CRC. The combined genotypes rs4938723 CC and TP53 GG had a tendency to protect individuals against CRC. This finding suggests that SNPs in the promoter regions of miRs may have important roles in the etiology of CRC, providing novel biomarkers for malignancies. This tendency was shown for cancer of the colon, particularly the proximal colon, but not for rectal cancer. Further investigation of the relationship between the site of cancer and genotype is required. Recently, studies reported that TP53 influences the expression of miR-34b/c in a one-way relationship. However, rs4938723 TT and TP53 Arg72Pro GC were also associated with decreased CRC risk in the present study, which suggests that TP53 and miR-34b/c may interact bidirectional and supports the concept that CRC is a complex disease involving multiple genes.
We also found that HTN or DM was associated with a significantly increased risk of CRC, regardless of genotype (Table IV). SNP rs4938723 with DM was associated with increased CRC risk when compared with wild-type, but TP53 Arg72Pro CC with DM was associated with decreased risk when compared with wild-type. Moreover, decreased levels of homocysteine and folic acid showed a positive correlation with the risk of CRC (Table VII). At folic acid levels of <4.45 ng/ml, both the rs4938723 and TP53 codon 72 SNPs were related to a decreased risk of CRC. At homocysteine levels >11.70 μmol/l, rs4938723 was positively correlated with CRC risk, but the TP53 codon 72 SNP was inversely correlated. These findings suggest that genetic factors and environmental factors, such as metabolic diseases, influence each other in CRC and may implicate miRs as pathophysiologic linkers/modulators between metabolic syndrome and cancer.
Rs4938723 is located within the CpG island of pri-miR-34b/c (423-bp upstream of the transcription start site), making it a potential binding site for GATA-X transcription factors. According to the web-based SNP analysis tool TFSEARCH 1.3, GATA family members bind to promoters of many genes and directly activate or suppress expression of target genes and may be involved in carcinogenesis (30). Therefore, SNP rs4938723 (T>C) may affect the expression level of miR-34b/c. Ectopic miR-34b/c caused cell cycle arrest in the G1 phase (31) and miR-34b/c inhibited proliferation and colony formation in soft agar (11). The relationship between miR expression in tumors and prognosis, with regard to both the treatment response rate and survival, is of increasing interest. This polymorphism has been reported to be associated with susceptibility to hepatocellular carcinoma (23), CRC (24), survival of breast cancer (26), renal cell carcinoma (32), non-small cell lung (33) and oral cancer (18).
Arg72Pro is the most widely investigated of the variations in the TP53 gene. The 72Arg allele induces apoptosis more efficiently than the 72Pro allele (34). It has been reported that homozygosity for Pro of TP53 Arg72Pro is potentially a risk factor for cancer of the lung, esophagus, stomach, breast, nasopharynx, urothelium and prostate (35–37). In CRC, meta-analysis was performed to estimate the effect of the TP53 Arg72Pro polymorphism and CRC risk, and no significant association was identified (38). In contrast to that result, we found a negative relationship between the TP53 Arg72Pro polymorphism and CRC risk in the Korean population.
Our study has several limitations. First, the enrolled patients were selected at a single institution in Korea and the sample size may limit the statistical power, especially for evaluating the connection between SNPs and CRC risk and survival in detailed genotype subgroups. Second, as the study was hospital-based, the cases and controls could not be representative of the general population. This possible selection bias could not be avoided. Finally, the mechanisms underlying the effects of genetic polymorphisms on the levels of pre/mature miRs and the identity of the miR target genes remain unknown and should be determined in future studies.
In conclusion, although we identified no significant prognostic value of SNPs for CRC, we found that SNPs rs4938723 and TP53 Arg72Pro show a trend toward improved survival and that the TP53 Arg72Pro CC genotype and dominant model (GC + CC) significantly decrease the risk of CRC. Additionally, the combined genotype rs4938723 and TP53 Arg72Pro (TT/GC and CC/GG) was significantly associated with decreased CRC risk in the Korean population. Also, SNP rs4938723 with DM was more closely related to increased CRC risk than the wild-type genotype. Future studies in multiethnic populations are warranted to validate our results and to define the functional effects of these SNPs on CRC.
Acknowledgements
This study was partially supported by a National Research Foundation (NRF) of Korea Grant funded by the Korean Government (2009-0075784 and 2012R1A1A2007033), and a Priority Research Centers Program Grant, administered by the NRF and funded by the Ministry of Education (2009-0093821), Republic of Korea.
References
Siegel R, Naishadham D and Jemal A: Cancer statistics, 2012. CA Cancer J Cin. 62:10–29. 2012. View Article : Google Scholar | |
Hrašovec S and Glavač D: MicroRNAs as novel biomarkers in colorectal cancer. Front Genet. 3:1802012. | |
Cunningham JM, Christensen ER, Tester DJ, et al: Hypermethylation of the hMLH1 promoter in colon cancer with microsatellite instability. Cancer Res. 58:3455–3460. 1998.PubMed/NCBI | |
Toyota M, Ahuja N, Ohe-Toyota M, Herman JG, Baylin SB and Issa JP: CpG island methylator phenotype in colorectal cancer. Proc Natl Acad Sci USA. 96:8681–8686. 1999. View Article : Google Scholar : PubMed/NCBI | |
Shen L, Toyota M, Kondo Y, et al: Integrated genetic and epigenetic analysis identifies three different subclasses of colon cancer. Proc Natl Acad Sci USA. 104:18654–18659. 2007. View Article : Google Scholar : PubMed/NCBI | |
Iorio MV and Croce CM: MicroRNA dysregulation in cancer: diagnostics, monitoring and therapeutics. A comprehensive review. EMBO Mol Med. 4:143–159. 2012. View Article : Google Scholar : PubMed/NCBI | |
Fazi F and Nervi C: MicroRNA: basic mechanisms and transcriptional regulatory networks for cell fate determination. Cardiovasc Res. 79:553–561. 2008. View Article : Google Scholar : PubMed/NCBI | |
Meister G and Tuschl T: Mechanisms of gene silencing by double-stranded RNA. Nature. 431:343–349. 2004. View Article : Google Scholar : PubMed/NCBI | |
Pillai RS, Bhattacharyya SN and Filipowicz W: Repression of protein synthesis by miRNAs: how many mechanisms? Trends Cell Biol. 17:118–126. 2007. View Article : Google Scholar : PubMed/NCBI | |
Lim LP, Lau NC, Garrett-Engele P, et al: Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature. 433:769–773. 2005. View Article : Google Scholar : PubMed/NCBI | |
Corney DC, Flesken-Nikitin A, Godwin AK, Wang W and Nikitin AY: MicroRNA-34b and MicroRNA-34c are targets of p53 and cooperate in control of cell proliferation and adhesion-independent growth. Cancer Res. 67:8433–8438. 2007. View Article : Google Scholar | |
He L, He X, Lim LP, et al: A microRNA component of the p53 tumour suppressor network. Nature. 447:1130–1134. 2007. View Article : Google Scholar : PubMed/NCBI | |
Tarasov V, Jung P, Verdoodt B, et al: Differential regulation of microRNAs by p53 revealed by massively parallel sequencing: miR-34a is a p53 target that induces apoptosis and G1-arrest. Cell Cycle. 6:1586–1593. 2007. View Article : Google Scholar : PubMed/NCBI | |
Chang TC, Wentzel EA, Kent OA, et al: Transactivation of miR-34a by p53 broadly influences gene expression and promotes apoptosis. Mol Cell. 26:745–752. 2007. View Article : Google Scholar : PubMed/NCBI | |
Welch C, Chen Y and Stallings RL: MicroRNA-34a functions as a potential tumor suppressor by inducing apoptosis in neuroblastoma cells. Oncogene. 26:5017–5022. 2007. View Article : Google Scholar : PubMed/NCBI | |
Bommer GT, Gerin I, Feng Y, et al: p53-mediated activation of miRNA34 candidate tumor-suppressor genes. Curr Biol. 17:1298–1307. 2007. View Article : Google Scholar : PubMed/NCBI | |
Toyota M, Suzuki H, Sasaki Y, Maruyama R, Imai K, Shinomura Y and Tokino T: Epigenetic silencing of microRNA-34b/c and B-cell translocation gene 4 is associated with CpG island methylation in colorectal cancer. Cancer Res. 68:4123–4132. 2008.PubMed/NCBI | |
Kozaki K, Imoto I, Mogi S, Omura K and Inazawa J: Exploration of tumor-suppressive microRNAs silenced by DNA hypermethylation in oral cancer. Cancer Res. 68:2094–2105. 2008. View Article : Google Scholar : PubMed/NCBI | |
Lujambio A, Calin GA, Villanueva A, et al: A microRNA DNA methylation signature for human cancer metastasis. Proc Natl Acad Sci USA. 105:13556–13561. 2008. View Article : Google Scholar : PubMed/NCBI | |
Kong YW, Ferland-McCollough D, Jackson TJ and Bushell M: microRNAs in cancer management. Lancet Ooncol. 13:e249–e258. 2012. View Article : Google Scholar : PubMed/NCBI | |
Lujambio A, Ropero S, Ballestar E, et al: Genetic unmasking of an epigenetically silenced microRNA in human cancer cells. Cancer Res. 67:1424–1429. 2007. View Article : Google Scholar : PubMed/NCBI | |
Iwai N and Naraba H: Polymorphisms in human pre-miRNAs. Biochem Biophys Res Commun. 331:1439–1444. 2005. View Article : Google Scholar : PubMed/NCBI | |
Xu Y, Liu L, Liu J, et al: A potentially functional polymorphism in the promoter region of miR-34b/c is associated with an increased risk for primary hepatocellular carcinoma. Int J Cancer. 128:412–417. 2011. View Article : Google Scholar : PubMed/NCBI | |
Gao LB, Li LJ, Pan XM, Li ZH, Liang WB, Bai P, Zhu YH and Zhang L: A genetic variant in the promoter region of miR-34b/c is associated with a reduced risk of colorectal cancer. Biol Chem. 394:415–420. 2013.PubMed/NCBI | |
Hiroki E, Suzuki F, Akahira J, Nagase S, Ito K, Sugawara J, Miki Y, Suzuki T, Sasano H and Yaegashi N: MicroRNA-34b functions as a potential tumor suppressor in endometrial serous adenocarcinoma. Int J Cancer. 131:E395–E404. 2012. View Article : Google Scholar : PubMed/NCBI | |
Bensen JT, Tse CK, Nyante SJ, Barnholtz-Sloan JS, Cole SR and Millikan RC: Association of germline microRNA SNPs in pre-miRNA flanking region and breast cancer risk and survival: the Carolina Breast Cancer Study. Cancer Causes Control. 24:1099–1109. 2013. View Article : Google Scholar : PubMed/NCBI | |
Ahmed RL, Schmitz KH, Anderson KE, Rosamond WD and Folsom AR: The metabolic syndrome and risk of incident colorectal cancer. Cancer. 107:28–36. 2006. View Article : Google Scholar : PubMed/NCBI | |
Moghaddam AA, Woodward M and Huxley R: Obesity and risk of colorectal cancer: a meta-analysis of 31 studies with 70,000 events. Cancer Epidemiol Biomarkers Prev. 16:2533–2547. 2007. View Article : Google Scholar : PubMed/NCBI | |
Wittke-Thompson JK, Pluzhnikov A and Cox NJ: Rational inferences about departures from Hardy-Weinberg equilibrium. Am J Hum Genet. 76:967–986. 2005. View Article : Google Scholar : PubMed/NCBI | |
Chou J, Provot S and Werb Z: GATA3 in development and cancer differentiation: cells GATA have it! J Cell Physiol. 222:42–49. 2010. | |
Xiang Y, Fan S, Cao J, Huang S and Zhang LP: Association of the microRNA-499 variants with susceptibility to hepatocellular carcinoma in a Chinese population. Mol Biol Rep. 39:7019–7023. 2012. View Article : Google Scholar : PubMed/NCBI | |
Lin J, Horikawa Y, Tamboli P, Clague J, Wood CG and Wu X: Genetic variations in microRNA-related genes are associated with survival and recurrence in patients with renal cell carcinoma. Carcinogenesis. 31:1805–1812. 2010. View Article : Google Scholar : PubMed/NCBI | |
Hu Z, Chen J, Tian T, et al: Genetic variants of miRNA sequences and non-small cell lung cancer survival. J Clin Invest. 118:2600–2608. 2008.PubMed/NCBI | |
Dumont P, Leu JI, Della Pietra AC III, George DL and Murphy M: The codon 72 polymorphic variants of p53 have markedly different apoptotic potential. Nat Genet. 33:357–365. 2003. View Article : Google Scholar : PubMed/NCBI | |
Irarrázabal CE, Rojas C, Aracena R, Márquez C and Gil L: Chilean pilot study on the risk of lung cancer associated with codon 72 polymorphism in the gene of protein p53. Toxicol Lett. 144:69–76. 2003.PubMed/NCBI | |
Kuroda Y, Tsukino H, Nakao H, Imai H and Katoh T: p53 codon 72 polymorphism and urothelial cancer risk. Cancer Lett. 189:77–83. 2003. View Article : Google Scholar | |
Suzuki K, Matsui H, Ohtake N, et al: A p53 codon 72 polymorphism associated with prostate cancer development and progression in Japanese. J Biomed Sci. 10:430–435. 2003. View Article : Google Scholar : PubMed/NCBI | |
Dahabreh IJ, Linardou H, Bouzika P, Varvarigou V and Murray S: TP53 Arg72Pro polymorphism and colorectal cancer risk: a systematic review and meta-analysis. Cancer Epidemiol Biomarkers Prev. 19:1840–1847. 2010. View Article : Google Scholar |