Abstract
Celiac disease (CD) is a gluten intolerance disorder with known genetic contribution. The recent fine mapping and genome-wide association studies (GWAS) have identified up to 57 non-HLA CD susceptibility SNPs, majority of which are non-coding variants lacking any functional annotation. Therefore, we adopted multidimensional computational approach for uncovering the plausible mechanisms through which these GWAS SNPs are connected to CD pathogenesis. At initial phase, we identified that 25 (43.85%) out of 57 CD-SNPs lies in evolutionarily constrained genetic element regions. In follow-up phases, through computational (CADD, GWAVA, and FATHMM algorithms) deleterious intensity measurements, we have discovered that 42 (3.94%) out of 1065 variants (57 CD-lead and 1008-linked SNPs; r(2)0.8) are differentially deleterious in nature to CD. Further functional scrutinization of these CD variants by public domain eQTL mapping, gene expression, knockout mouse model, and pathway analyses revealed that deleterious SNPs of CCR2 gene influences its expression levels and may also elicit a cascade of T-cell-mediated immunological events leading to intestinal gluten intolerance in genetically susceptible individuals. This study demonstrates the utility of integrated in silico analysis of annotations, gene expression, and pathways in prioritizing the potential complex disease variants from large-scale open source genomic data. J. Cell. Biochem. 118: 2193-2207, 2017. (c) 2017 Wiley Periodicals, Inc.