A functional variant in GREM1 confers risk for colorectal cancer by disrupting a hsa-miR-185-3p binding site

The transforming growth factor beta (TGF-β) pathway has been implicated in carcinogenesis of intestinal canal. Except for common variants indentified by genome-wide association studies, variants with lower frequency can also explain a part of the disease heritability, especially those in gene regulatory regions. In this study, we searched for colorectal cancer (CRC) related functional low-frequency variants (minor allele frequency 1-5%) in untranslated regions (UTR) involved in the TGF-β signaling using a next-generation sequencing based approach. A case-control study including 1,841 CRC cases and 1,837 controls was performed to identify CRC associated variants and biological experiments were applied to further explore the potential functions of the significant variants. Three low-frequency UTR variants were selected as our candidates and subsequent association analyses showed that a low-frequency variant rs12915554 in the 3’ UTR of GREM1 was significantly associated with CRC risk (Additive model: OR=1.43, 95%CI: 1.04-1.95, P=0.026). Functional annotations suggested that rs12915554 variation increased the expression of GREM1 by perturbing a hsa-miR-185-3p binding site. Moreover, higher expression level of GREM1 was investigated in colon tumor tissues compared with adjacent normal tissues using TCGA data. In conclusion, low-frequency UTR variant rs12915554 in the gene GREM1 was in relation to CRC susceptibility in a Chinese population and this variation might promote CRC development through enhancing GREM1 expression in a miRNA-mediated posttranscriptional manner.