Genetic variants in m6A modification genes are associated with colorectal cancer risk.

The N6-methyladenosine (m6A) modification plays important regulatory roles in gene expression, cancer occurrence and metastasis. Herein, we aimed to explore the association between genetic variants in m6A modification genes and susceptibility to colorectal cancer. We used logistic regression models to investigate the associations between candidate single-nucleotide polymorphisms (SNPs) in 20 m6A modification genes and colorectal cancer risk. The false discovery rate (FDR) method was used for multiple comparisons. Dual luciferase assays and RNA m6A quantifications were applied to assess transcriptional activity and measure m6A levels, respectively. We found that SND1 rs118049207 was significantly associated with colorectal cancer risk in a Nanjing population (odds ratio (OR) = 1.69, 95% confidence interval (95% CI) = 1.31-2.18, P = 6.51 × 10-6). This finding was further replicated in an independent Beijing population (OR = 1.36, 95% CI = 1.04-1.79, P = 2.41 × 10-2) and in a combined analysis (OR = 1.52, 95% CI = 1.27-1.84, P = 8.75 × 10-6). Stratification and interaction analyses showed that SND1 rs118049207 multiplicatively interacted with the sex and drinking status of the patients to enhance their colorectal cancer risk (P = 1.56 × 10-3 and 1.41 × 10-2, respectively). Furthermore, rs118049207 served as an intronic enhancer on SND1 driven by DMRT3. SND1 mRNA expression was markedly increased in colorectal tumour tissues compared with adjacent normal tissues. The colorimetric m6A quantification strategy revealed that SND1 could alter m6A levels in colorectal cancer cell lines. Our findings indicated that genetic variants in m6A modification genes might be promising predictors of colorectal cancer risk.