Discovering Genetic Factors for psoriasis through exhaustively searching for significant second order SNP-SNP interactions

In this paper, we aim at discovering genetic factors of psoriasis through searching for statistically significant SNP-SNP interactions exhaustively from two real psoriasis genome-wide association study datasets (phs000019.v1.p1 and phs000982.v1.p1) downloaded from the database of Genotypes and Phenotypes. To deal with the enormous search space, our search algorithm is accelerated with eight biological plausible interaction patterns and a pre-computed look-up table. After our search, we have discovered several SNPs having a stronger association to psoriasis when they are in combination with another SNP and these combinations may be non-linear interactions. Among the top 20 SNP-SNP interactions being found in terms of pairwise p-value and improvement metric value, we have discovered 27 novel potential psoriasis-associated SNPs where most of them are reported to be eQTLs of a number of known psoriasis-associated genes. On the other hand, we have inferred a gene network after selecting the top 10000 SNP-SNP interactions in terms of improvement metric value and we have discovered a novel long distance interaction between XXbac-BPG154L12.4 and RNU6-283P which is not a long distance haplotype and may be a new discovery. Finally, our experiments with the synthetic datasets have shown that our pre-computed look-up table technique can significantly speed up the search process.

[1]  James T. Elder,et al.  Fine mapping of eight psoriasis susceptibility loci , 2014, European Journal of Human Genetics.

[2]  K. P. O'Brien,et al.  HLA-Cw*0602 associates more strongly to psoriasis in the Swedish population than variants of the novel 6p21.3 gene PSORS1C3. , 2005, Acta dermato-venereologica.

[3]  Matti Pirinen,et al.  A genome-wide association study identifies new psoriasis susceptibility loci and an interaction between HLA-C and ERAP1 , 2010, Nature Genetics.

[4]  Jack T Stapleton,et al.  The Major Genetic Determinants of HIV-1 Control Affect HLA Class I Peptide Presentation , 2010, Science.

[5]  Yongxiang Yu,et al.  Investigation of 20 non‐HLA (human leucocyte antigen) psoriasis susceptibility loci in Chinese patients with psoriatic arthritis and psoriasis vulgaris , 2013, The British journal of dermatology.

[6]  Weixiong Zhang,et al.  Allele-Specific Network Reveals Combinatorial Interaction That Transcends Small Effects in Psoriasis GWAS , 2014, PLoS Comput. Biol..

[7]  S. Tsai,et al.  The genetic polymorphisms of POU5F1 gene are associated with psoriasis vulgaris in Chinese. , 2007, Journal of dermatological science.

[8]  X. Estivill,et al.  Coding haplotype analysis supports HCR as the putative susceptibility gene for psoriasis at the MHC PSORS1 locus. , 2002, Human molecular genetics.

[9]  Kristel Van Steen,et al.  A roadmap to multifactor dimensionality reduction methods , 2015, Briefings Bioinform..

[10]  James T. Elder,et al.  Multiple Loci within the Major Histocompatibility Complex Confer Risk of Psoriasis , 2009, PLoS genetics.

[11]  Scott M. Williams,et al.  Widespread epistasis regulates glucose homeostasis and gene expression , 2017, PLoS Genetics.

[12]  Jason H. Moore,et al.  BIOINFORMATICS REVIEW , 2005 .

[13]  S. Rosset,et al.  Measuring missing heritability: Inferring the contribution of common variants , 2014, Proceedings of the National Academy of Sciences.

[14]  Christian von Mering,et al.  STRING 8—a global view on proteins and their functional interactions in 630 organisms , 2008, Nucleic Acids Res..

[15]  Alan M. Kwong,et al.  Next-generation genotype imputation service and methods , 2016, Nature Genetics.

[16]  Scott F. Saccone,et al.  A Genome-Wide Association Study of Psoriasis and Psoriatic Arthritis Identifies New Disease Loci , 2008, PLoS genetics.

[17]  Gonçalo Abecasis,et al.  Deletion of the late cornified envelope LCE3B and LCE3C genes as a susceptibility factor for psoriasis , 2009, Nature Genetics.

[18]  I. Kockum,et al.  Genetic Association with ERAP1 in Psoriasis Is Confined to Disease Onset after Puberty and Not Dependent on HLA-C*06 , 2012, The Journal of investigative dermatology.

[19]  Xuejun Zhang,et al.  HLA‐DQA1 and DQB1 alleles are associated with genetic susceptibility to psoriasis vulgaris in Chinese Han , 2004, International journal of dermatology.

[20]  G. Rocheleau,et al.  A survey about methods dedicated to epistasis detection , 2015, Front. Genet..

[21]  Chris S. Haley,et al.  Detecting epistasis in human complex traits , 2014, Nature Reviews Genetics.

[22]  Cheng Soon Ong,et al.  GWIS - model-free, fast and exhaustive search for epistatic interactions in case-control GWAS , 2013, BMC Genomics.

[23]  Liangdan Sun,et al.  Epigenome-wide association data implicates DNA methylation-mediated genetic risk in psoriasis , 2016, Clinical Epigenetics.

[24]  Zhong Zhou,et al.  LncReg: a reference resource for lncRNA-associated regulatory networks , 2015, Database J. Biol. Databases Curation.

[25]  J. Shendure,et al.  A general framework for estimating the relative pathogenicity of human genetic variants , 2014, Nature Genetics.

[26]  Can Yang,et al.  GBOOST: a GPU-based tool for detecting gene-gene interactions in genome-wide case control studies , 2011, Bioinform..

[27]  Paul Weston,et al.  Interaction between ERAP1 and HLA-B27 in ankylosing spondylitis implicates peptide handling in the mechanism for HLA-B27 in disease susceptibility , 2011, Nature Genetics.

[28]  M. Weale,et al.  Conditional analysis identifies three novel major histocompatibility complex loci associated with psoriasis. , 2012, Human molecular genetics.

[29]  Lin He,et al.  SHEsisEpi, a GPU-enhanced genome-wide SNP-SNP interaction scanning algorithm, efficiently reveals the risk genetic epistasis in bipolar disorder , 2010, Cell Research.

[30]  A. Barton,et al.  One SNP at a Time: Moving beyond GWAS in Psoriasis. , 2016, The Journal of investigative dermatology.

[31]  H. Inoko,et al.  Major histocompatibility complex class I chain-related gene A in Thai psoriasis patients: MICA association as a part of human leukocyte antigen-B-Cw haplotypes. , 2004, Tissue antigens.

[32]  K. P. O'Brien,et al.  Polymorphisms in the SEEK1 and SPR1 genes on 6p21.3 associate with psoriasis in the Swedish population , 2003, Experimental dermatology.

[33]  Association of SEEK1 and psoriatic arthritis in two distinct Canadian populations. , 2005, Annals of the rheumatic diseases.

[34]  Attila Balint,et al.  Systematic analysis of complex genetic interactions , 2018, Science.

[35]  Q. Zou,et al.  An overview of SNP interactions in genome-wide association studies. , 2015, Briefings in functional genomics.

[36]  Jane Worthington,et al.  Genotypic variability-based genome-wide association study identifies non-additive loci HLA-C and IL12B for psoriasis , 2017, Journal of Human Genetics.

[37]  Wei Jiang,et al.  A network approach to exploring the functional basis of gene‐gene epistatic interactions in disease susceptibility , 2018, Bioinform..

[38]  James T. Elder,et al.  Genome-wide Association Analysis of Psoriatic Arthritis and Cutaneous Psoriasis Reveals Differences in Their Genetic Architecture. , 2015, American journal of human genetics.

[39]  Cisca Wijmenga,et al.  The MHC locus and genetic susceptibility to autoimmune and infectious diseases , 2017, Genome Biology.

[40]  Momiao Xiong,et al.  Psoriasis prediction from genome-wide SNP profiles , 2011, BMC dermatology.

[41]  T. G. Kim,et al.  MICA 5.1 allele is a susceptibility marker for psoriasis in the Korean population. , 2000, Tissue antigens.

[42]  Bertil Schmidt,et al.  GPU-accelerated exhaustive search for third-order epistatic interactions in case-control studies , 2015, J. Comput. Sci..

[43]  Chris S. Haley,et al.  EpiGPU: exhaustive pairwise epistasis scans parallelized on consumer level graphics cards , 2011, Bioinform..

[44]  H. Yu,et al.  Investigating the genetic association of HCP5, SPATA2, TNIP1, TNFAIP3 and COG6 with psoriasis in Chinese population , 2014, International journal of immunogenetics.

[45]  A. Skytthe,et al.  Heritability of psoriasis in a large twin sample , 2013, The British journal of dermatology.

[46]  Ellen T. Gelfand,et al.  The Genotype-Tissue Expression (GTEx) project , 2013, Nature Genetics.

[47]  M. Daly,et al.  Genome-wide association studies for common diseases and complex traits , 2005, Nature Reviews Genetics.

[48]  Bertil Schmidt,et al.  FPGA-based Acceleration of Detecting Statistical Epistasis in GWAS , 2014, ICCS.

[49]  Stefan Jenisch,et al.  Sequence and haplotype analysis supports HLA-C as the psoriasis susceptibility 1 gene. , 2006, American journal of human genetics.

[50]  Yadong Wang,et al.  LncRNA2Target: a database for differentially expressed genes after lncRNA knockdown or overexpression , 2014, Nucleic Acids Res..

[51]  H. Inoko,et al.  HLA Class I and II Alleles and Susceptibility to Generalized Pustular Psoriasis: Significant Associations with HLA‐Cw1 and HLA‐DQB1*0303 , 1998, The Journal of dermatology.

[52]  Momiao Xiong,et al.  A Novel Statistic for Genome-Wide Interaction Analysis , 2010, PLoS genetics.

[53]  G. Gandhi,et al.  Psoriasis-Associated Genetic Polymorphism in North Indian Population in the CCHCR1 Gene and in a Genomic Segment Flanking the HLA-C Region , 2011, Disease markers.

[54]  Mike Tyers,et al.  BioGRID: a general repository for interaction datasets , 2005, Nucleic Acids Res..

[55]  A. Bowcock,et al.  The immunogenetics of Psoriasis: A comprehensive review. , 2015, Journal of autoimmunity.

[56]  B. Cohen,et al.  Genetic Interactions Between Transcription Factors Cause Natural Variation in Yeast , 2009, Science.

[57]  D. Gladman,et al.  HLA antigens in psoriatic arthritis. , 1986, The Journal of rheumatology.

[58]  M. Ritchie.  Finding the epistasis needles in the genome-wide haystack. , 2015, Methods in molecular biology.

[59]  C. Carcassi,et al.  Psoriasis is associated with a SNP haplotype of the corneodesmosin gene (CDSN). , 2002, Tissue antigens.

[60]  K. Wuepper,et al.  Psoriasis vulgaris: a genetic approach. , 1990, The Journal of investigative dermatology.

[61]  Pui-Yan Kwok,et al.  Genomewide Scan Reveals Association of Psoriasis with IL-23 and NF-κB Pathways , 2008, Nature Genetics.

[62]  Liangdan Sun,et al.  Psoriasis Regression Analysis of MHC Loci Identifies Shared Genetic Variants with Vitiligo , 2011, PloS one.

[63]  A. Bowcock,et al.  The Molecular Revolution in Cutaneous Biology: The Era of Genome-Wide Association Studies and Statistical, Big Data, and Computational Topics. , 2017, The Journal of investigative dermatology.

[64]  Dipanwita Roy Chowdhury,et al.  Human protein reference database as a discovery resource for proteomics , 2004, Nucleic Acids Res..

[65]  J. Kere,et al.  Clinical associations of the risk alleles of HLA-Cw6 and CCHCR1*WWCC in psoriasis. , 2007, Acta dermato-venereologica.

[66]  S. Heath,et al.  ADAM33, a New Candidate for Psoriasis Susceptibility , 2007, PloS one.

[67]  Kenny Q. Ye,et al.  An integrated map of genetic variation from 1,092 human genomes , 2012, Nature.

[68]  A. Morris,et al.  Data quality control in genetic case-control association studies , 2010, Nature Protocols.

[69]  Bertil Schmidt,et al.  Hybrid CPU/GPU Acceleration of Detection of 2-SNP Epistatic Interactions in GWAS , 2014, Euro-Par.

[70]  H. Kreysel,et al.  HLA-linked complement polymorphisms (C2, BF) in psoriasis , 2004, Archives of Dermatological Research.

[71]  J. Ocampo-Candiani,et al.  Candidate gene polymorphisms and risk of psoriasis: A pilot study , 2016, Experimental and therapeutic medicine.

[72]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[73]  Feng Xu,et al.  Gene-gene Interaction Analysis by IAC (Interaction Analysis by Chi-Square) - A Novel Biological Constraint-based Interaction Analysis Framework , 2016, BIOINFORMATICS.

[74]  Xia Li,et al.  RAID: a comprehensive resource for human RNA-associated (RNA–RNA/RNA–protein) interaction , 2014, RNA.

[75]  Melissa A. Basford,et al.  Systematic comparison of phenome-wide association study of electronic medical record data and genome-wide association study data , 2013, Nature Biotechnology.

[76]  Jun Zhu,et al.  Development of GMDR-GPU for Gene-Gene Interaction Analysis and Its Application to WTCCC GWAS Data for Type 2 Diabetes , 2013, PloS one.

[77]  R. Segal,et al.  Polymorphism in MICA rather than HLA-B/C genes is associated with psoriatic arthritis in the Jewish population. , 2001, Human immunology.

[78]  Y. Chang,et al.  A study of candidate genes for psoriasis near HLA‐C in Chinese patients with psoriasis , 2003, The British journal of dermatology.

[79]  Guoliang Chen,et al.  Deletion of the late cornified envelope genes LCE3C and LCE3B is associated with psoriasis in a Chinese population. , 2011, The Journal of investigative dermatology.

[80]  S. Tsai,et al.  Psoriasis vulgaris in Chinese individuals is associated with PSORS1C3 and CDSN genes , 2006, The British journal of dermatology.

[81]  A. Barton,et al.  Confirmation of TNIP1 and IL23A as susceptibility loci for psoriatic arthritis , 2011 .

[82]  Yusuke Nakamura,et al.  Genome-wide association study identifies variations in 6p21.3 associated with nevirapine-induced rash. , 2011, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[83]  J. Reveille,et al.  Psoriasis and psoriatic arthritis associated with human immunodeficiency virus infection. , 1991, Rheumatic diseases clinics of North America.