TYK2 Protein-Coding Variants Protect against Rheumatoid Arthritis and Autoimmunity, with No Evidence of Major Pleiotropic Effects on Non-Autoimmune Complex Traits

Despite the success of genome-wide association studies (GWAS) in detecting a large number of loci for complex phenotypes such as rheumatoid arthritis (RA) susceptibility, the lack of information on the causal genes leaves important challenges to interpret GWAS results in the context of the disease biology. Here, we genetically fine-map the RA risk locus at 19p13 to define causal variants, and explore the pleiotropic effects of these same variants in other complex traits. First, we combined Immunochip dense genotyping (n = 23,092 case/control samples), Exomechip genotyping (n = 18,409 case/control samples) and targeted exon-sequencing (n = 2,236 case/controls samples) to demonstrate that three protein-coding variants in TYK2 (tyrosine kinase 2) independently protect against RA: P1104A (rs34536443, OR = 0.66, P = 2.3x10-21), A928V (rs35018800, OR = 0.53, P = 1.2x10-9), and I684S (rs12720356, OR = 0.86, P = 4.6x10-7). Second, we show that the same three TYK2 variants protect against systemic lupus erythematosus (SLE, Pomnibus = 6x10-18), and provide suggestive evidence that two of the TYK2 variants (P1104A and A928V) may also protect against inflammatory bowel disease (IBD; Pomnibus = 0.005). Finally, in a phenome-wide association study (PheWAS) assessing >500 phenotypes using electronic medical records (EMR) in >29,000 subjects, we found no convincing evidence for association of P1104A and A928V with complex phenotypes other than autoimmune diseases such as RA, SLE and IBD. Together, our results demonstrate the role of TYK2 in the pathogenesis of RA, SLE and IBD, and provide supporting evidence for TYK2 as a promising drug target for the treatment of autoimmune diseases.

Joshua C. Denny | Lisa Bastarache | Jing Cui | David E. Larson | Elaine R. Mardis | Shawn N. Murphy | Helena Canhão | Robert M. Plenge | Yukinori Okada | Katherine P. Liao | Isaac Kohane | Michael F. Seldin | Soumya Raychaudhuri | Andrew Cagan | Peter K. Gregersen | Robert R. Graham | Tushar R. Bhangale | Elizabeth W. Karlson | John Bowes | Marieke J. H. Coenen | Vivian Gainer | Anne Barton | Barbara Franke | Roman Kosoy | Timothy Behrens | Alexandra Zhernakova | Jane Worthington | Xavier Mariette | Annette Lee | P. Gregersen | L. Kiemeney | B. Franke | L. Fulton | R. Fulton | E. Mardis | A. Zhernakova | Y. Okada | A. Barton | S. Eyre | Jane Worthington | D. Larson | I. Kohane | S. Raychaudhuri | C. Fronick | M. O'Laughlin | J. Denny | V. Gainer | S. Murphy | E. Stahl | J. Kremer | K. Liao | E. Karlson | R. Plenge | P. Tak | M. Seldin | J. Bowes | T. Behrens | Javier Martín | R. Graham | W. Ortmann | T. Bhangale | J. Greenberg | D. Pappas | D. Diogo | J. Cui | M. Coenen | X. Mariette | N. de Vries | L. Padyukov | L. Bastarache | H. Canhão | J. Crusius | A. Cagan | C. Richard-Miceli | J. Fonseca | M. Nurmohamed | F. Kurreeman | T. Mikuls | Tracie L. Deluca | R. Kosoy | M. Ransom | Lars Klareskog | Javier Martin | Steve Eyre | Leonid Padyukov | Paul P. Tak | Ward Ortmann | Annette Lee | Robert S. Fulton | Lambertus A. Kiemeney | Niek de Vries | Dorothée Diogo | Fina Kurreeman | Eli A. Stahl | J. Bart A. Crusius | Catrina C. Fronick | Jeffrey D. Greenberg | Joel M. Kremer | Dimitrios A. Pappas | João E. Fonseca | Michael Ransom | Michelle O'Laughlin | Lars Klareskog | Corrine Richard-Miceli | Michael T. Nurmohamed | Ted R. Mikuls | Lucinda L. Fulton | A. Barton | Michael Ransom

[1]  Judy H. Cho,et al.  Deep Resequencing of GWAS Loci Identifies Rare Variants in CARD9, IL23R and RNF186 That Are Associated with Ulcerative Colitis , 2013, PLoS genetics.

[2]  Y. Okada,et al.  Genome-wide association studies to advance our understanding of critical cell types and pathways in rheumatoid arthritis: recent findings and challenges , 2014, Current opinion in rheumatology.

[3]  P. Bork,et al.  A method and server for predicting damaging missense mutations , 2010, Nature Methods.

[4]  R. Plenge,et al.  JAK and STAT signaling molecules in immunoregulation and immune-mediated disease. , 2012, Immunity.

[5]  Brian T Fay,et al.  Associations of disease activity and treatments with mortality in men with rheumatoid arthritis: results from the VARA registry. , 2011, Rheumatology.

[6]  Isabelle Cleynen,et al.  Resequencing of positional candidates identifies low frequency IL23R coding variants protecting against inflammatory bowel disease , 2011, Nature Genetics.

[7]  Matti Pirinen,et al.  Identification of 15 new psoriasis susceptibility loci highlights the role of innate immunity , 2012 .

[8]  Marylyn D. Ritchie,et al.  PheWAS: demonstrating the feasibility of a phenome-wide scan to discover gene–disease associations , 2010, Bioinform..

[9]  M. Petri,et al.  Longitudinal expression of type I interferon responsive genes in systemic lupus erythematosus , 2009, Lupus.

[10]  J. Kremer,et al.  The CORRONA database. , 2006, Autoimmunity reviews.

[11]  D. Altshuler,et al.  Validating therapeutic targets through human genetics , 2013, Nature Reviews Drug Discovery.

[12]  M. Stephens,et al.  Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. , 2003, Genetics.

[13]  Miles Parkes,et al.  Genetic insights into common pathways and complex relationships among immune-mediated diseases , 2013, Nature Reviews Genetics.

[14]  J. Casanova,et al.  A patient with tyrosine kinase 2 deficiency without hyper-IgE syndrome. , 2012, The Journal of pediatrics.

[15]  V. Sexl,et al.  Tyrosine kinase 2 (TYK2) in cytokine signalling and host immunity. , 2011, Frontiers in bioscience.

[16]  Christian Gieger,et al.  Loci influencing lipid levels and coronary heart disease risk in 16 European population cohorts , 2009, Nature Genetics.

[17]  Shamil R Sunyaev,et al.  Pooled association tests for rare variants in exon-resequencing studies. , 2010, American journal of human genetics.

[18]  Inês Barroso,et al.  Genetic Variants Influencing Circulating Lipid Levels and Risk of Coronary Artery Disease , 2010, Arteriosclerosis, thrombosis, and vascular biology.

[19]  Christopher A. Miller,et al.  VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing. , 2012, Genome research.

[20]  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.

[21]  A. Syvänen,et al.  Association of NCF2, IKZF1, IRF8, IFIH1, and TYK2 with Systemic Lupus Erythematosus , 2011, PLoS genetics.

[22]  H. H. Kuper,et al.  The efficacy of anti-TNF in rheumatoid arthritis, a comparison between randomised controlled trials and clinical practice , 2007, Annals of the rheumatic diseases.

[23]  Philip Rosenstiel,et al.  Genome-wide association study for Crohn's disease in the Quebec Founder Population identifies multiple validated disease loci , 2007, Proceedings of the National Academy of Sciences.

[24]  T. Pincus,et al.  Databases of patients with early rheumatoid arthritis in the USA. , 2003, Clinical and experimental rheumatology.

[25]  K. Garber.  Pfizer's first-in-class JAK inhibitor pricey for rheumatoid arthritis market , 2013, Nature Biotechnology.

[26]  Suzette J. Bielinski,et al.  Use of diverse electronic medical record systems to identify genetic risk for type 2 diabetes within a genome-wide association study , 2012, J. Am. Medical Informatics Assoc..

[27]  Jonathan C. Cohen,et al.  Sequence variations in PCSK9, low LDL, and protection against coronary heart disease. , 2006, The New England journal of medicine.

[28]  A. Dunning,et al.  Beyond GWASs: illuminating the dark road from association to function. , 2013, American journal of human genetics.

[29]  C. Spettell,et al.  The comparative risk of serious infections among rheumatoid arthritis patients starting or switching biological agents , 2011, Annals of the rheumatic diseases.

[30]  L. Cardon,et al.  Use of genome-wide association studies for drug repositioning , 2012, Nature Biotechnology.

[31]  H. Hakonarson,et al.  ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data , 2010, Nucleic acids research.

[32]  H. Drummond,et al.  Haplotype-tagging analysis of common variants of the IL23R gene demonstrates gene-wide extent of association with IBD. , 2013, Inflammatory bowel diseases.

[33]  Lars Alfredsson,et al.  A Candidate Gene Approach Identifies the TRAF1/C5 Region as a Risk Factor for Rheumatoid Arthritis , 2007, PLoS Medicine.

[34]  M. Pangalos,et al.  Lessons learned from the fate of AstraZeneca's drug pipeline: a five-dimensional framework , 2014, Nature Reviews Drug Discovery.

[35]  Peter Szolovits,et al.  Association between low density lipoprotein and rheumatoid arthritis genetic factors with low density lipoprotein levels in rheumatoid arthritis and non-rheumatoid arthritis controls , 2013, Annals of the rheumatic diseases.

[36]  H. Canhão,et al.  Reuma.pt - the rheumatic diseases portuguese register. , 2011, Acta reumatologica portuguesa.

[37]  S. Browning,et al.  A Groupwise Association Test for Rare Mutations Using a Weighted Sum Statistic , 2009, PLoS genetics.

[38]  Pak Chung Sham,et al.  Genetic Power Calculator: design of linkage and association genetic mapping studies of complex traits , 2003, Bioinform..

[39]  D. Reich,et al.  Principal components analysis corrects for stratification in genome-wide association studies , 2006, Nature Genetics.

[40]  Annette Lee,et al.  CSK regulatory polymorphism is associated with systemic lupus erythematosus and influences B cell signaling and activation , 2012, Nature Genetics.

[41]  C. Carlson,et al.  Genetic variants associated with the white blood cell count in 13,923 subjects in the eMERGE Network , 2011, Human Genetics.

[42]  Tariq Ahmad,et al.  Genome-wide meta-analysis increases to 71 the number of confirmed Crohn's disease susceptibility loci , 2010, Nature Genetics.

[43]  T. Morio,et al.  Human tyrosine kinase 2 deficiency reveals its requisite roles in multiple cytokine signals involved in innate and acquired immunity. , 2006, Immunity.

[44]  J. Allison,et al.  Administrative codes combined with medical records based criteria accurately identified bacterial infections among rheumatoid arthritis patients. , 2009, Journal of clinical epidemiology.

[45]  J. Todd,et al.  Rare Variants of IFIH1, a Gene Implicated in Antiviral Responses, Protect Against Type 1 Diabetes , 2009, Science.

[46]  Wentian Li,et al.  Genome-Wide Association Scan Identifies Candidate Polymorphisms Associated with Differential Response to Anti-TNF Treatment in Rheumatoid Arthritis , 2008, Molecular medicine.

[47]  P. Gaffney,et al.  Two functional lupus-associated BLK promoter variants control cell-type- and developmental-stage-specific transcription. , 2014, American journal of human genetics.

[48]  Christian Gieger,et al.  Multiple Loci Are Associated with White Blood Cell Phenotypes , 2011, PLoS genetics.

[49]  D. Clayton,et al.  The imprinted DLK1-MEG3 gene region on chromosome 14q32.2 alters susceptibility to type 1 diabetes , 2009, Nature Genetics.

[50]  R. Collins,et al.  Common variants at 30 loci contribute to polygenic dyslipidemia , 2009, Nature Genetics.

[51]  M. Ban,et al.  A rare variant of the TYK2 gene is confirmed to be associated with multiple sclerosis , 2010, European Journal of Human Genetics.

[52]  Elodie Ey,et al.  Genetic and Functional Analyses of SHANK2 Mutations Suggest a Multiple Hit Model of Autism Spectrum Disorders , 2012, PLoS genetics.

[53]  Peter Szolovits,et al.  Genetic basis of autoantibody positive and negative rheumatoid arthritis risk in a multi-ethnic cohort derived from electronic health records. , 2011, American journal of human genetics.

[54]  Hongyu Zhao,et al.  Rare independent mutations in renal salt handling genes contribute to blood pressure variation , 2008, Nature Genetics.

[55]  R. Dixon,et al.  Enabling structure-based drug design of Tyk2 through co-crystallization with a stabilizing aminoindazole inhibitor , 2012, BMC Structural Biology.

[56]  Thomas Meitinger,et al.  Loss-of-function mutations in SLC30A8 protect against type 2 diabetes , 2014, Nature Genetics.

[57]  Ken Garber,et al.  Pfizer's JAK inhibitor sails through phase 3 in rheumatoid arthritis , 2011, Nature Biotechnology.

[58]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[59]  J. Fletcher,et al.  EMerging BiomARKers in Inflammatory Bowel Disease (EMBARK) Study Identifies Fecal Calprotectin, Serum MMP9, and Serum IL-22 as a Novel Combination of Biomarkers for Crohn's Disease Activity: Role of Cross-Sectional Imaging , 2013, The American Journal of Gastroenterology.

[60]  S. Purcell,et al.  Pleiotropy in complex traits: challenges and strategies , 2013, Nature Reviews Genetics.

[61]  Inês Barroso,et al.  Rare MTNR1B variants impairing melatonin receptor 1B function contribute to type 2 diabetes , 2012, Nature Genetics.

[62]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[63]  L. Rönnblom,et al.  Two Rare Disease-Associated Tyk2 Variants Are Catalytically Impaired but Signaling Competent , 2013, The Journal of Immunology.

[64]  J. Kere,et al.  Polymorphisms in the tyrosine kinase 2 and interferon regulatory factor 5 genes are associated with systemic lupus erythematosus. , 2005, American journal of human genetics.

[65]  P. Gregersen,et al.  The North American Rheumatoid Arthritis Consortium--bringing genetic analysis to bear on disease susceptibility, severity, and outcome. , 1998, Arthritis care and research : the official journal of the Arthritis Health Professions Association.

[66]  J. French,et al.  A phase IIb dose-ranging study of the oral JAK inhibitor tofacitinib (CP-690,550) versus placebo in combination with background methotrexate in patients with active rheumatoid arthritis and an inadequate response to methotrexate alone. , 2012, Arthritis and rheumatism.

[67]  X. Montalban,et al.  Tyrosine kinase 2 variant influences T lymphocyte polarization and multiple sclerosis susceptibility. , 2011, Brain : a journal of neurology.

[68]  S. Yusuf,et al.  Excess of Rare Variants in Non–Genome-Wide Association Study Candidate Genes in Patients With Hypertriglyceridemia , 2012, Circulation. Cardiovascular genetics.

[69]  Jun S. Liu,et al.  Genetics of rheumatoid arthritis contributes to biology and drug discovery , 2013 .

[70]  N. Nishimoto,et al.  Clinical value of blocking IL-6 receptor , 2009, Current opinion in rheumatology.

[71]  A. Koch.  The pathogenesis of rheumatoid arthritis. , 2007, American journal of orthopedics.

[72]  I. Kohane,et al.  Electronic medical records for discovery research in rheumatoid arthritis , 2010, Arthritis care & research.

[73]  P. Gregersen,et al.  MHC associations with clinical and autoantibody manifestations in European SLE , 2014, Genes and Immunity.

[74]  Steven J. Schrodi,et al.  A PheWAS approach in studying HLA-DRB1*1501 , 2013, Genes and Immunity.

[75]  Joshua M. Korn,et al.  Deep resequencing of GWAS loci identifies independent rare variants associated with inflammatory bowel disease , 2011, Nature Genetics.

[76]  J. Kere,et al.  Evidence for Genetic Association and Interaction Between the TYK2 and IRF5 Genes in Systemic Lupus Erythematosus , 2009, The Journal of Rheumatology.

[77]  Melissa A. Basford,et al.  Genome- and Phenome-Wide Analyses of Cardiac Conduction Identifies Markers of Arrhythmia Risk , 2013, Circulation.

[78]  D. Roden,et al.  Development of a Large‐Scale De‐Identified DNA Biobank to Enable Personalized Medicine , 2008, Clinical pharmacology and therapeutics.

[79]  Eli Stahl,et al.  High density genetic mapping identifies new susceptibility loci for rheumatoid arthritis , 2012, Nature Genetics.

[80]  Jennifer G. Robinson,et al.  The interleukin-6 receptor as a target for prevention of coronary heart disease: a mendelian randomisation analysis , 2012, The Lancet.

[81]  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.

[82]  B. Browning,et al.  Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. , 2007, American journal of human genetics.

[83]  David C. Wilson,et al.  Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease , 2012, Nature.

[84]  A. Barton,et al.  Genome-Wide Association Study of Genetic Predictors of Anti–Tumor Necrosis Factor Treatment Efficacy in Rheumatoid Arthritis Identifies Associations with Polymorphisms at Seven Loci , 2011, Arthritis and rheumatism.

[85]  S. Henikoff,et al.  Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm , 2009, Nature Protocols.

[86]  M. Carroll,et al.  Trends in lipids and lipoproteins in US adults, 1988-2010. , 2012, JAMA.

[87]  Erick R. Scott,et al.  Rare coding variants in the phospholipase D3 gene confer risk for Alzheimer’s disease , 2013, Nature.

[88]  P. Ravaud,et al.  A single tumour necrosis factor haplotype influences the response to adalimumab in rheumatoid arthritis , 2007, Annals of the rheumatic diseases.

[89]  James T. Elder,et al.  Rare and common variants in CARD14, encoding an epidermal regulator of NF-kappaB, in psoriasis. , 2012, American journal of human genetics.

[90]  D. Barford,et al.  KSR2 Mutations Are Associated with Obesity, Insulin Resistance, and Impaired Cellular Fuel Oxidation , 2013, Cell.

[91]  Sivakumar Gowrisankar,et al.  Rare variants in CFI, C3 and C9 are associated with high risk of advanced age-related macular degeneration , 2013, Nature Genetics.

[92]  S. Hebbring.  The challenges, advantages and future of phenome-wide association studies , 2014, Immunology.

[93]  Annette Lee,et al.  Analysis of Families in the Multiple Autoimmune Disease Genetics Consortium (madgc) Collection: the Ptpn22 620w Allele Associates with Multiple Autoimmune Phenotypes , 2022 .

[94]  Xihong Lin,et al.  Optimal tests for rare variant effects in sequencing association studies. , 2012, Biostatistics.

[95]  Christopher G. Chute,et al.  A genome- and phenome-wide association study to identify genetic variants influencing platelet count and volume and their pleiotropic effects , 2013, Human Genetics.