The present study showed a higher frequency of rs1333040 CT and TT genotypes in the CRC group and also this variant was associated with a higher risk of CRC in the co-dominant, dominant, and log-additive models. However, the ANRIL rs10757274, rs10757274, and rs1333048 polymorphisms showed no association with CRC in all genetic models. The CGAA haplotype of ANRIL rs1333040/ rs10757274/ rs4977574/rs1333048 polymorphisms was associated with the TGGC haplotype and lower risk of CRC. In addition, the rs1333040 polymorphism (CT and TT genotypes) was associated with higher TNM stages (III and IV in CRC patients in the dominant model. The ANRIL rs10757274 polymorphism was lower in CRC patients over the age of 50 years only in the dominant model. The ANRIL rs10757274 variant was associated with an increased risk of differentiated CRC tumors compared to poor or moderate in the dominant model.
ANRIL is a well-known lncRNA (also called CDKN2BAS) and encodes a large antisense non–coding RNA. ANRIL affects several biological processes and regulates the processes that can play key roles in tumor promotion and progression. Therefore, it is a potential candidate in various diseases like cancers. GWAS has introduced several genetic variants in INK4b-ARF-INK4a (such as ANRIL SNPs) and proposed the impacts of these variants on the susceptibility to the development of several diseases like cancers[12].
Important roles of ANRIL in the invasion and metastasis of TC cells were described by Zhao et al. They indicated that ANRIL could reduce the expression of p15INK4b tumor suppressor by blocking the TGF-β /Smad signaling pathway [19]. The results obtained from Liu et al.’s study indicated the up-regulation of ANRIL in association with poor prognosis and tumorigenesis of oral squamous cell carcinoma (OSCC)[20]. In Vitro declined proliferation, migration, and invasion of osteosarcoma cells in knockdown of ANRIL expression condition were found by Cheng et al.’s study demonstrating its role in cellular biology and oncogenesis [21]. Likewise, higher expression of ANRIL in hepatocellular carcinoma is reported by Huang et al. [22]. In many studies, ANRIL was introduced as an oncogene involved in several cancers, such as gastric cancer[23], breast cancer[24], lung cancer[25] and thyroid cancer[11].
Regarding the effects of genetic factors on cancer pathogenesis, several studies investigated the effects of genetic variants in different genes on the promotion and progression of cancer. Long non-coding RNAs like ANRIL were among the candidate genes [9, 11].
In a study on genome-wide association performed by Bahrami et al, the common genetic variants in a region on chromosome 9p21 were introduced as a risk factor for CRC. The genes in this locus were CDKN2B, CDKN2A, and 3' end of CDKN2BAS (ANRIL), and several genetic polymorphisms were recognized in this region. ANRIL showed post-transcriptional modulatory activity and disturbed the expression of other genes near itself [26]. In the large meta-analysis performed by Timofeeva et al., a novel locus for squamous cell carcinoma at 9p21, including the ANRIL rs1333040 variant was also introduced [27].
The effects of ANRIL polymorphism on various diseases have been studied and the relation between these variants and several cancers has been indicated. However, there is no published report on the possible effects of ANRIL variants and colorectal cancer. Kang et al. (2015) showed no association between ANRIL rs2151280 T/C polymorphism and esophageal cancer risk [28]. In their study, Deng et al. (2019) found that the ANRIL rs2151280 variant genotype (A > G) was related to the susceptibility to glioma [29]. In another study, Zhou et al. indicated that the ANRIL rs1537373 variant might affect transcription factor binding and control CDKN2B expression, associated with the predisposition to pancreatic cancer (PC)[30].
Omrani et al. performed the random forest (RF) algorithm to predict prostate cancer risk and showed that ANRIL rs1333048 AA and rs10757278 GG could predict the risk of prostate cancer in the Iranian population [31]. In addition, Gong et al showed no relation between ANRIL rs10120688 and rs1333049 polymorphisms and lung susceptibility; however, they indicated that ANRIL rs10120688 polymorphism was associated with response to platinum-based chemotherapy [32].
A meta-analysis was performed by Huang et al. (2018) on 17 studies investigating the effects of long non-coding RNA polymorphisms on cancer risk. The ANRIL rs4977574, rs1333048, and rs10757278 polymorphisms but not rs1333045 were accompanied by overall cancer risk[33]. The ANRIL rs17694493 GG genotype was associated with 4.93- 4.93-fold increased risk of gastric cancer as compared with the CC genotype in Petkevicius et al’s study [23]. Chen et al. indicated the ANRIL rs1333049 as a biomarker for better prognosis for lung cancer patients with platinum-based chemotherapy [25]. The possible effects of ANRIL rs1333045, rs4977574, rs1333048, and rs10757278 genotypes on prostate cancer were examined by Taheri et al.’s study performed in Iran and despite the results of the present study, they indicated the relation between rs4977574, rs1333048 and rs10757278 polymorphisms but not rs1333045 and PC risk [34]. In their study, Yao et al. evaluated the relation between ANRIL rs1333045, rs4977574, rs1333048, and rs10757278 polymorphisms and cervical cancer (CC) and they showed only the effect of rs4977574A allele and lower risk of CC[35]. Khorshidi et al. studied the relation between ANRIL rs1333045, rs1333048, rs4977574, and rs10757278 variants and breast cancer (BC) and found that there was a relation between rs1333045TT genotype and BC before multiple-testing correction. In addition, they showed a higher frequency of TCGA haplotype in BC patients[24].
In another study, Maruei-Milan et al. found no association between none of the ANRIL rs11333048, rs4977574, rs1333040, and rs10757274 SNPs and papillary thyroid cancer (PTC). However, the AAAC haplotype showed a decreased risk of PTC, and the CAAC, and also CAGT haplotypes were associated with an increased risk of cancer. In this study, the rs1333048 CC genotype was higher in patients with larger tumors in the recessive model. The rs4977574 AC genotype was associated with smaller tumor size in the over-dominant model. rs10757274 AA and rs1333040 CC genotypes were associated with a decreased risk of III-IV cancer stages[11].
In the present study, for the first time, we studied the effects of ANRIL polymorphisms on colorectal cancer. However, our study had some limitations, first, a small sample size, particularly the analysis of subgroups might affect the accuracy of the present results. Second, the results would have been more valuable if the study had been performed on cancerous tissue and healthy lateral tissue about ANRIL expression. Further studies are needed to assay these aspects