Polymorphism (rs16917496) at the miR-502 Binding Site of the Lysine Methyltransferase 5A (SET8) and Its Correlation with Colorectal Cancer in Iranians.

BACKGROUND: One of the gene expression regulatory mechanisms is mediated by small noncoding RNAs called microRNA (miRNA). They interact with a recognition sequence located mostly in 3'-untranslated regions (3'-UTRs) of mRNAs. Polymorphisms in miRNAs recognition sequences could affect gene expression which in turn may alter disease susceptibility. SET8, a member of the SET domain-containing methyltransferase, acts in a variety of biological processes such as genomic stability. Here, we report correlation of rs16917496 polymorphism, located in the recognition sequence of miR-502 within 3'-UTR of SET8, with colorectal cancer (CRC) in Iranians.
MATERIALS AND METHODS: One hundred and seventy CRC patients and 170 noncancer counterparts were recruited in this case-control study. Genotyping of rs16917496 was performed using polymerase chain reaction-restriction fragment length polymorphism method.
RESULTS: There was no significant association of rs16917496 with CRC in population under study (P value for genotype and allele distribution were >0.05). However, stratification analysis based on smoking status revealed that TT+TC genotypes of SET8 rs16917496 are strongly associated with increased risk of CRC (odds ratio: 5.8, 95% confidence interval: 1.37-24.34, P - 0.005) in smoker subgroup.
CONCLUSION: Correlation of rs16917496 T allele with CRC in smokers is emphasizing the importance of individuals' genotype in the recruitment of adverse health hazards of smoking more profoundly for certain people compared to others.

Introduction

Colorectal cancer (CRC) is the second most common cancer in women and the third prevalent one in men in developed countries.[123] According to the National Institute of Cancer Statistics as of 2013, 45 out of 100,000 people diagnosed with CRC in the world.[45] Although the incidence of CRC is lower in Asian than in Western countries, recent studies have shown increasing rates of CRC in developing countries including Iran.[467]

Genetic syndromes such as lynch and familial adenomatous polyposis, together with positive family history constituting up to 30% of CRC susceptibility; however, common genetic variants present in individuals’ genome considered as the main genetic risk of CRC.[8] Recent genome-wide association studies determined multiple CRC-related single nucleotide polymorphisms (SNPs).[910111213]

MicroRNAs (miRNAs) constitute a novel discovered a class of small noncoding RNAs (about 22 nucleotides long) that play important roles in posttranscriptional gene regulation by binding to the 3’-untranslated regions (3’-UTRs) of specific mRNAs resulting in mRNA cleavage or translational repression.[4141516] SNPs located in miRNAs-binding sites contribute to the disruption of miRNA recognition elements (MREs) or create new sites which eventually can lead to gene expression dysregulation and disease susceptibility.[21718]

The previous studies have suggested that the SNP rs16917496 (CC genotype), which is located within the miR-502 binding site of SET8 (also known as SETD8 and PR-SET7) 3’UTR, is correlated with reduced SET8 protein levels.[192021] SET8 is known as a histone H4-Lys-20 specific methyltransferase and has been recently reported that modulates p53 activity by particularly monomethylation p53 at lysine 382 (p53K382me1).[2223] SET8 is involved in vital cellular processes such as transcriptional regulation, genomic stability, heterochromatin formation, DNA replication, and cell cycle arrest.[242526] Therefore, it is proposed that deregulation of SET8 may affect individual's cancer risk. Considering regulatory effects of SET8 on p53 activity, in conditions such as DNA damage, with decreased expression of SET8, the proapoptotic and checkpoint activation functions of p53 would be enhanced consequently. However, CC genotype related to SNP rs16917496 could disrupt this regulation.[222728]

According to the available reports, this polymorphism is associated with breast, ovarian, lymphoblastic leukemia and cervical cancer susceptibility, and clinical outcome of hepatocellular carcinoma (HCC).[192027282930] However, CRC susceptibility in relation to SNP rs16917496 has not worked out yet.

Here, we report the results of our case–control study on the correlation of SNP rs16917496 with sporadic CRC in the Iranian population. In addition, other factors such as gender, age, body mass index (BMI), physical activity, smoking status, and nonsteroidal anti-inflammatory drugs (NSAIDs) consumption were also analyzed in the present study.

Materials and Methods

Study population

For this case–control study, a total of 340 individuals were selected among patients referred to the colonoscopy centers of the Isfahan University of Medical Sciences Hospitals, from mid-2013 to mid-2015. One hundred and seventy CRC cases were selected on the basis of diagnosis by colonoscopy and then followed their pathology report for the ultimate confirmation of their colonoscopy-based CRC diagnosis. Equal numbers of individuals who underwent colonoscopy procedure but diagnosed with no malignancy or other adverse colonic diseases were constituted our normal control group. All participants filled and signed an informed consent form according to the Isfahan University of Medical Sciences Ethical Committee instructions. To eliminate any established genetic risk factor for CRC, patients were individually interviewed to select cases with sporadic CRC without positive history of familial cancers. In addition, all participants were asked to fill up a structured questionnaire to register the parameters known to influence the CRC susceptibility risk including gender, age, BMI, physical activity, smoking status, and NSAIDs consumption (for NSAIDs consumption any person who use aspirin or nonaspirin NSAIDs occasional or daily for long term [>5 years] were considered as an NSAIDs consumer or regular consumer).

Single nucleotide polymorphism selection

Figure 1 summarizes the course of SNP selection in SET8 gene. As it is evident, SNPs located at 3’UTR of the target gene were considered for this study. Minor allele frequency higher than 0.05 and being located at MREs of miRNAs were the other two criteria considered for SNP selection. SNP rs16917496 was well fitted into our criteria and hence selected for this study.

Figure 1

Single nucleotide polymorphism selection algorithm

Single nucleotide polymorphism genotyping

Genomic DNA was extracted from peripheral blood using PrimePrep Genomic DNA Isolation Kit (Genet Bio, Korea). The quality and quantity of the extracted DNA were assessed by agarose gel electrophoresis and spectrophotometry. SNP genotyping performed by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method. PCR product size, restriction enzyme used, fragments subsequent to digestion, and primer sequences used are given in Table 1. PCR was performed in a total volume of 25 μl, including 2.5 μl of 10 × PCR buffer, 1 mM MgCl2, 200 μM each dNTPs, 0.2 μM each primer, 100 ng of genomic DNA, and 1.5 U of Taq DNA polymerase (provided from CinnaGen company, Iran). Thermal cycling condition included an initial denaturation of 94°C for 5 min and then 34 cycles of 94°C, 62°C, and 72°C, all for 30 s and a final extension at 72°C for 5 min. RFLP was carried out using SwaI restriction endonuclease enzyme (provided from Thermo Fisher Scientific Company, USA). Subsequently, digested products were differentiated on agarose gel to determine the genotype of each sample. Fragment sizes of 390 and 360 bp were indicative of TT homozygous genotype; a single 750 bp fragment displayed the presence of CC homozygous and three fragments of 390, 360, and 750 bp indicated TC genotype [Figure 2]. Confirmation of the RFLP-based genotyping was done by direct sequencing of randomly selected samples with different genotypes (10% of samples).

Table 1

Primer sequences and characteristic of polymerase chain reaction product and digestion reaction

Figure 2

Agarose gel electrophoresis analysis of restriction fragment length polymorphism-polymerase chain reaction products. Lines 1 and 8 are TT genotype, 2–7 and 10 are TC genotype, and 9 is CC genotype

Statistics

All statistical analyses were performed using SPSS version 22 (SPSS, Chicago, IL, USA). Genotype frequencies in cases and controls were tested for Hardy–Weinberg equilibrium using Chi-square test. Associations between SET8 (rs16917496) polymorphism with susceptibility to sporadic CRC were examined by logistic regression analysis. Age and sex were used as confounders. Correlation between this polymorphism and CRC was evaluated using odds ratios (ORs) and 95% confidence intervals (CIs). The significance level was set at P < 0.05. The difference in demographic and lifestyle characteristics distribution such as age, gender, smoking status, BMI, and NSAIDs consumption assessed by Pearson Chi-square test for categorical variables and t-test for continuous variables. Mann–Whitney test were used to compare physical activity between CRC and control groups.

Results

Demographic and lifestyle characteristics

Our study was intended to evaluate the relationship between rs16917496 SNP located at 3’UTR of SET8 gene and CRC susceptibility in a subset of the Iranian population. In this case–control study, 170 confirmed CRC patients (90 males and 80 females) and 170 (66 males and 104 females) unaffected individuals were selected among those referred to colonoscopy units of our University hospitals. Mean age in case and control groups was 57.72 ± 12.29 and 49.05 ± 13.79 years, respectively. Adjusted data (by age and gender) demonstrated that there was a considerable difference between case and control groups regarding physical activity (P < 0.001). Moreover, we found that healthy subjects in control group had more NSAIDs consumption compared with CRC patients group (P - 0.004). However, there were no statistically significant differences between patients and controls in terms of BMI and smoking status (P - 0.19 and P - 0.23, respectively). Characteristics of the study population are summarized in Table 2.

Table 2

Baseline characteristics of colorectal cancer patients and controls in study

Genotype and allele distribution

The genotype distribution of rs16917496 polymorphism in cases and controls was in agreement with Hardy–Weinberg equilibrium. The frequency of TT, CC, and TC genotype in control group was 40.6%, 17.6%, and 41.8%, respectively and this frequency in case group was 34.1%, 18.8%, and 47.1%, respectively. In addition, the frequencies of T and C alleles were 61.5% and 38.5% in controls, and 57.6% and 42.4% in cases, respectively. Statistical analysis revealed that no significant correlation exists between SNP rs16917496 and risk of CRC in Iranian population (P value for genotype and allele distribution, respectively: 0.32 and 0.15) [Table 3]. Because C is the ancestral allele for this SNP, we compared other two genotypes against CC genotype. It is determined that TT + TC genotype (variant genotypes) has no either pathologic or protective effect on sporadic CRC (OR: 1.08, 95% CI: 0.62–1.88). In addition, association between these genotypes (TT+TC) compared to CC with CRC risk was assessed in subgroups of participants stratified by age (under 55 and over 55) and BMI (under 25 and over 25 kg/m2). Our data indicated that there was no significant difference in two age groups (OR<55: 1.1, 95% CI: 0.51–2.35 and OR≥55: 1.14, 95% CI: 0.49–2.66) with genotype distribution. In addition, data for BMI stratification demonstrated that there was no correlation between BMI and the target genotypes (OR≤25: 0.78, 95% CI: 0.33–1.86 and OR>25:1.34, 95% CI: 0.64–2.82).

Table 3

Association between genotypes and allele frequency with colorectal cancer risk

Further analyses were carried out to evaluate the association between TT+TC compared to CC and CRC stratified by NSAIDs consuming and smoking status. Stratification based on irregular and regular NSAID subgroups displayed no correlation between these subgroups and SNP genotypes regarding to CRC risk (ORirregular: 1.08, 95% CI: 0.59–2.00 and ORregular: 1.14, 95% CI: 0.28–4.52). TT+TC genotypes were determined to be associated with increased risk of CRC (OR: 5.8, 95% CI: 1.37–24.34) in smoker subgroup. While in nonsmoker subgroup, no significant correlation was determined (OR: 1.61, 95% CI: 0.86–3.02) [Table 4].

Table 4

Stratification analysis of rs16917496 genotype frequency in sporadic colorectal cancer and control groups

Discussion

SNPs are the most frequent type of variation in the human genome.[31] SNPs in the miRNA-binding sites at 3’UTR region can modulate miRNA-binding by disrupting existing recognition elements (MREs) or creating new recognition sites.[1718] Altered mRNA-miRNA interaction would affect the amount of particular gene translation products available in the cell at a given time point. However, the outcome of such deregulation would manifest as disease susceptibility such as malignancies.[143233]

The previous studies demonstrated that SNP rs16917496 (CC genotype and C allele) which is located within the SET8 miR-502 binding site in its 3’UTR, is correlated with reduced SET8 protein levels.[2728] Here, we report an association between SNP rs16917496 genotypes and sporadic CRC risk in Iranian population. SET8 encodes a histone H4 methyltransferase that is essential for various biological processes, such as transcriptional regulation, heterochromatin formation, DNA replication, cell-cycle arrest, and maintenance of genome integrity.[2334] SET8 methyltransferase activity modifies p53 protein function by monomethylation of p53 at lysine 382.[22] In cellular stresses including DNA damage, p53-dependent cell cycle arrest or apoptosis activated for which depletion of cellular SET8 is an essential prerequisite as seen that p53K382me1 levels decrease with DNA damage.[2227] Given, the importance of SET8-P53 interaction in cell-cycle and apoptosis control, deregulation of SET8 under influence of SNP rs16917496 T>C within the miR-502 binding site in the 3’-UTR of the SET8 gene would adversely affect this vital axis in the cell.[20]

Our data are indicative of no significant association between allele and genotype frequency of this polymorphism and sporadic CRC risk in Iranian population. Reviewing previous studies regarding rs16917496 T>C shows contradictory results. In Chinese population, CC and CC+CT genotypes are reported to be associated with longer survival of HCC, nonsmall-cell lung carcinoma (NSCLC) (CC genotype), and SCLC (CC+CT genotype) patients.[202735] In addition, the other study in china indicated that CC genotype is associated with decreased risk of epithelial ovarian cancer and TT genotype is associated with increased risk of NSCLC.[2936] On the other hand, a few studies revealed that CC genotype was associated with increased risk of breast cancer among the premenopausal women in the Chinese population, and CT+TT genotype is associated with decreased risk of acute lymphoblastic leukemia in comparison with CC genotype in the Iranian population.[1928] Despite these inconsistent results, all analysis on SET8 expression indicated that C allele is associated with reduced SET8 protein level.[2027]

Our results revealed no correlation of this polymorphism with sporadic CRC in Iranian population. In stratification analysis, neither CC genotype nor combined TT and TC genotypes, did not alter the risk of CRC among individuals under and over 55 years (P - 0.40 and P - 0.38, respectively), irregular and regular NSAIDs consumers (P - 0.39 and P - 0.42, respectively), and BMI under 25 and over 25 kg/m2 (P - 0.27 and P - 0.22, respectively). However, in the subgroup of smokers, combined TT+TC genotype was associated with increased risk of sporadic CRC. However, there may be other players taking part in SET8 regulatory pathway. For instance, we can mention SCF (β-TRCP) which is marking SET8 for ubiquitination and degradation in a casein kinase I-mediated manner.[37]

Our finding regarding rs16917496 T allele correlation with CRC in smokers compared to nonsmokers, would be an additional proof of cigarette health hazard. However, our findings regarding genotype-smoking interaction highlight this important point that certain genotypes are more readily accommodate adverse health effect of smoking compared to the other genotypes. Finally, it is the first evaluation of SNP rs16917496 located within the SET8 miR-502 binding site in relation with CRC. However, the reasons for the contradictory results may be owing to differences in ethnic background, and the different forms of CRC entities (familial vs. sporadic) and low-sample size that needs to more investigation in other population and more sample size.

Financial support and sponsorship

The study was supported by Isfahan University of Medical Sciences.

Conflicts of interest

There are no conflicts of interest.