Abstract
The genome-wide association studies (GWAS) have enabled us in identifying different breast cancer (BC) susceptibility loci. However, majority of these are non-coding variants with no annotated biological function. We investigated such 78 noncoding genome wide associated SNPs of BC and further expanded the list to 2,162 variants with strong linkage-disequilibrium (LD, r
≥0.8). Using multiple publically available algorithms such as CADD, GWAVA, and FATHAMM, we classified all these variants into deleterious, damaging, or benign categories. Out of total 2,241 variants, 23 (1.02%) variants were extreme deleterious (rank 1), 70 (3.12%) variants were deleterious (rank 2), and 1,937 (86.43%) variants were benign (rank 3). The results show 14% of lead or associated variants are under strong negative selection (GERP++ RS ≥2), and ∼22% are under balancing selection (Tajima's D score >2) in CEU population of 1KGP-the regions being positively selected (GERP++ RS <0) in mammalian evolution. The expression quantitative trait loci of highest deleteriously ranked genes were tested on relevant adipose and breast tissues, the results of which were extended for protein expression on breast tissues. From the concordance analysis of ranking system of GWAVA, CADD, and FATHMM, eQTL and protein expression, we identified the deleterious SNPs localized in STXBP4 and ZNF404 genes which might play a role in BC development by dysregulating its gene expression. This simple approach will be easier to implement and to prioritize large scale GWAS data for variety of diseases and link to the potentially unrecognized functional roles of genes. J. Cell. Biochem. 118: 4296-4307, 2017. © 2017 Wiley Periodicals, Inc.