Novel approach to functional SNPs discovery from genome-wide data reveals promising variants for colon cancer risk

In the majority of colorectal cancer (CRC) cases, the genetic basis of predisposition remains unexplained. The goal of the study was to assess the regulatory SNPs (rSNPs) in the human genome and to reveal СRC drivers based on the available chromatin immunoprecipitation sequencing (ChIP-Seq, ChIA-PET) and transcriptional profiling (RNA-Seq) data. We combined positional (locations within genome regulatory elements) and functional (associated with allele-specific binding and expression) criteria followed by an analysis using genome-wide association studies (GWAS) and minor allele frequency (MAF) datasets. DeSeq2 analysis through 70 CRC patients reinforced the regulatory potential. rSNPs (1,476) that were associated with significant (P < 0.01) allele-specific events resulting in thirty that exhibited a link with CRC according to the MAF and 27, with a risk of malignancy in general according to GWAS. Selected rSNPs may modify the expression of genes for tumor suppressors and the regulators of signaling pathways, including noncoding RNAs. However, the rSNPs from the most represented group affect the expression of genes related to splicing. Our findings strongly suggest that the identified variants might contribute to CRC susceptibility, which indicates that aberrant splicing is one of the key mechanisms for unraveling disease etiopathogenesis and provides useful inputs for interpreting how genotypic variation corresponds to phenotypic outcome.