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.3 x 10(-21)), A928V (rs35018800, OR = 0.53, P = 1.2 x 10(-9)), and I684S (rs12720356, OR = 0.86, P = 4.6 x 10(-7)). Second, we show that the same three TYK2 variants protect against systemic lupus erythematosus (SLE, Pomnibus = 6 x 10(-18)), and provide suggestive evidence that two of the TYK2 variants (P1104A and A928V) may also protect against inflammatory bowel disease (IBD; P(omnibus) = 0.005). Finally, in a phenome-wide association study (PheWAS) assessing \u3e500 phenotypes using electronic medical records (EMR) in \u3e29,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

Distribution and consequences of VKORC1 polymorphisms in Germany

Runge, M., Von Keyserlingk, M., Braune, S., Freise, J., Eiler, T., Plenge-Bönig, A., Becker, D., Pelz, H.-J., Esther, A., Rost, S., Müller, C.R

Integration of Sequence Data from a Consanguineous Family with Genetic Data from an Outbred Population Identifies PLB1 as a Candidate Rheumatoid Arthritis Risk Gene

Integrating genetic data from families with highly penetrant forms of disease together with genetic data from outbred populations represents a promising strategy to uncover the complete frequency spectrum of risk alleles for complex traits such as rheumatoid arthritis (RA). Here, we demonstrate that rare, low-frequency and common alleles at one gene locus, phospholipase B1 (PLB1), might contribute to risk of RA in a 4-generation consanguineous pedigree (Middle Eastern ancestry) and also in unrelated individuals from the general population (European ancestry). Through identity-by-descent (IBD) mapping and whole-exome sequencing, we identified a non-synonymous c.2263G. C (p.G755R) mutation at the PLB1 gene on 2q23, which significantly co-segregated with RA in family members with a dominant mode of inheritance (P = 0.009). We further evaluated PLB1 variants and risk of RA using a GWAS meta-analysis of 8,875 RA cases and 29,367 controls of European ancestry. We identified significant contributions of two independent non-coding variants near PLB1 with risk of RA (rs116018341 [MAF = 0.042] and rs116541814 [MAF = 0.021], combined P = 3.26 x 10(-6)). Finally, we performed deep exon sequencing of PLB1 in 1,088 RA cases and 1,088 controls (European ancestry), and identified suggestive dispersion of rare protein-coding variant frequencies between cases and controls (P = 0.049 for C-alpha test and P = 0.055 for SKAT). Together, these data suggest that PLB1 is a candidate risk gene for RA. Future studies to characterize the full spectrum of genetic risk in the PLB1 genetic locus are warranted

Visualization of Shared Genomic Regions and Meiotic Recombination in High-Density SNP Data

A fundamental goal of single nucleotide polymorphism (SNP) genotyping is to determine the sharing of alleles between individuals across genomic loci. Such analyses have diverse applications in defining the relatedness of individuals (including unexpected relationships in nominally unrelated individuals, or consanguinity within pedigrees), analyzing meiotic crossovers, and identifying a broad range of chromosomal anomalies such as hemizygous deletions and uniparental disomy, and analyzing population structure.We present SNPduo, a command-line and web accessible tool for analyzing and visualizing the relatedness of any two individuals using identity by state. Using identity by state does not require prior knowledge of allele frequencies or pedigree information, and is more computationally tractable and is less affected by population stratification than calculating identity by descent probabilities. The web implementation visualizes shared genomic regions, and generates UCSC viewable tracks. The command-line version requires pedigree information for compatibility with existing software and determining specified relationships even though pedigrees are not required for IBS calculation, generates no visual output, is written in portable C++, and is well-suited to analyzing large datasets. We demonstrate how the SNPduo web tool identifies meiotic crossover positions in siblings, and confirm our findings by visualizing meiotic recombination in synthetic three-generation pedigrees. We applied SNPduo to 210 nominally unrelated Phase I / II HapMap samples and, consistent with previous findings, identified six undeclared pairs of related individuals. We further analyzed identity by state in 2,883 individuals from multiplex families with autism and identified a series of anomalies including related parents, an individual with mosaic loss of chromosome 18, an individual with maternal heterodisomy of chromosome 16, and unexplained replicate samples.SNPduo provides the ability to explore and visualize SNP data to characterize the relatedness between individuals. It is compatible with, but distinct from, other established analysis software such as PLINK, and performs favorably in benchmarking studies for the analyses of genetic relatedness

Psychiatric co-morbidity is associated with increased risk of surgery in Crohn's disease

Psychiatric co-morbidity, in particular major depression and anxiety, is common in patients with Crohn's disease (CD) and ulcerative colitis (UC). Prior studies examining this may be confounded by the co-existence of functional bowel symptoms. Limited data exist examining an association between depression or anxiety and disease-specific endpoints such as bowel surgery.National Institutes of Health (U.S.) (NIH U54-LM008748)American Gastroenterological AssociationNational Institutes of Health (U.S.) (NIH K08 AR060257)Beth Isreal Deaconess Medical Center (Katherine Swan Ginsburg Fund)National Institutes of Health (U.S.) (NIH R01-AR056768)National Institutes of Health (U.S.) (NIH U01-GM092691)National Institutes of Health (U.S.) (NIH R01-AR059648)Burroughs Wellcome Fund (Career Award for Medical Scientists)National Institutes of Health (U.S.) (NIH K24 AR052403)National Institutes of Health (U.S.) (NIH P60 AR047782)National Institutes of Health (U.S.) (NIH R01 AR049880

Single nucleotide polymorphisms at the TRAF1/C5 locus are associated with rheumatoid arthritis in a Han Chinese population

<p>Abstract</p> <p>Background</p> <p>Genetic variants in <it>TRAF1C5 </it>and <it>PTPN22 </it>genes have been shown to be significantly associated with arthritis rheumatoid in Caucasian populations. This study investigated the association between single nucleotide polymorphisms (SNPs) in <it>TRAF1/C5 </it>and <it>PTPN22 </it>genes and rheumatoid arthritis (RA) in a Han Chinese population. We genotyped SNPs rs3761847 and rs7021206 at the <it>TRAF1/C5 </it>locus and rs2476601 SNP in the <it>PTPN22 </it>gene in a Han Chinese cohort composed of 576 patients with RA and 689 controls. The concentrations of anti-cyclic citrullinated peptide antibodies (CCP) and rheumatoid factor (RF) were determined for all affected patients. The difference between the cases and the controls was compared using <it>χ</it>²analysis.</p> <p>Results</p> <p>Significant differences in SNPs rs3761847 and rs7021206 at <it>TRAF1/C5 </it>were observed between the case and control groups in this cohort; the allelic p-value was 0.0018 with an odds ratio of 1.28 for rs3761847 and 0.005 with an odds ratio of 1.27 for rs7021206. This significant association between rs3761847 and RA was independent of the concentrations of anti-CCP and RF. No polymorphism of rs2476601 was observed in this cohort.</p> <p>Conclusions</p> <p>We first demonstrated that genetic variants at the <it>TRAF1/C5 </it>locus are significantly associated with RA in Han Chinese, suggesting that <it>TRAF1/C5 </it>may play a role in the development of RA in this population, which expands the pathogenesis role of <it>TRAF1/C5 </it>in a different ethnicity.</p

Balancing selection is common in the extended MHC region but most alleles with opposite risk profile for autoimmune diseases are neutrally evolving

<p>Abstract</p> <p>Background</p> <p>Several susceptibility genetic variants for autoimmune diseases have been identified. A subset of these polymorphisms displays an opposite risk profile in different autoimmune conditions. This observation open interesting questions on the evolutionary forces shaping the frequency of these alleles in human populations.</p> <p>We aimed at testing the hypothesis whereby balancing selection has shaped the frequency of opposite risk alleles.</p> <p>Results</p> <p>Since balancing selection signatures are expected to extend over short genomic portions, we focused our analyses on 11 regions carrying putative functional polymorphisms that may represent the disease variants (and the selection targets). No exceptional nucleotide diversity was observed for <it>ZSCAN23</it>, <it>HLA-DMB</it>, <it>VARS2</it>, <it>PTPN22</it>, <it>BAT3</it>, <it>C6orf47</it>, and <it>IL10</it>; summary statistics were consistent with evolutionary neutrality for these gene regions. Conversely, <it>CDSN/PSORS1C1</it>, <it>TRIM10/TRIM40</it>, <it>BTNL2</it>, and <it>TAP2 </it>showed extremely high nucleotide diversity and most tests rejected neutrality, suggesting the action of balancing selection. For <it>TAP2 </it>and <it>BTNL2 </it>these signatures are not secondary to linkage disequilibrium with HLA class II genes. Nonetheless, with the exception of variants in <it>TRIM40 </it>and <it>CDSN</it>, our data suggest that opposite risk SNPs are not selection targets but rather have accumulated as neutral variants.</p> <p>Conclusion</p> <p>Data herein indicate that balancing selection is common within the extended MHC region and involves several non-HLA loci. Yet, the evolutionary history of most SNPs with an opposite effect for autoimmune diseases is consistent with evolutionary neutrality. We suggest that variants with an opposite effect on autoimmune diseases should not be considered a distinct class of disease alleles from the evolutionary perspective and, in a few cases, the opposite effect on distinct diseases may derive from complex haplotype structures in regions with high genetic diversity.</p

Defining the role of the MHC in autoimmunity: a review and pooled analysis

Abstract: The major histocompatibility complex (MHC) is one of the most extensively studied regions in the human genome because of the association of variants at this locus with autoimmune, infectious, and inflammatory diseases. However, identification of causal variants within the MHC for the majority of these diseases has remained difficult due to the great variability and extensive linkage disequilibrium (LD) that exists among alleles throughout this locus, coupled with inadequate study design whereby only a limited subset of about 20 from a total of approximately 250 genes have been studied in small cohorts of predominantly European origin. We have performed a review and pooled analysis of the past 30 years of research on the role of the MHC in six genetically complex disease traits -multiple sclerosis (MS), type 1 diabetes (T1D), systemic lupus erythematosus (SLE), ulcerative colitis (UC), Crohn&apos;s disease (CD), and rheumatoid arthritis (RA) -in order to consolidate and evaluate the current literature regarding MHC genetics in these common autoimmune and inflammatory diseases. We corroborate established MHC disease associations and identify predisposing variants that previously have not been appreciated. Furthermore, we find a number of interesting commonalities and differences across diseases that implicate both general and disease-specific pathogenetic mechanisms in autoimmunity

Genome-Wide Association Study and Gene Expression Analysis Identifies CD84 as a Predictor of Response to Etanercept Therapy in Rheumatoid Arthritis

Anti-tumor necrosis factor alpha (anti-TNF) biologic therapy is a widely used treatment for rheumatoid arthritis (RA). It is unknown why some RA patients fail to respond adequately to anti-TNF therapy, which limits the development of clinical biomarkers to predict response or new drugs to target refractory cases. To understand the biological basis of response to anti-TNF therapy, we conducted a genome-wide association study (GWAS) meta-analysis of more than 2 million common variants in 2,706 RA patients from 13 different collections. Patients were treated with one of three anti-TNF medications: etanercept (n = 733), infliximab (n = 894), or adalimumab (n = 1,071). We identified a SNP (rs6427528) at the 1q23 locus that was associated with change in disease activity score (ΔDAS) in the etanercept subset of patients (P = 8×10-8), but not in the infliximab or adalimumab subsets (P>0.05). The SNP is predicted to disrupt transcription factor binding site motifs in the 3′ UTR of an immune-related gene, CD84, and the allele associated with better response to etanercept was associated with higher CD84 gene expression in peripheral blood mononuclear cells (P = 1×10-11 in 228 non-RA patients and P = 0.004 in 132 RA patients). Consistent with the genetic findings, higher CD84 gene expression correlated with lower cross-sectional DAS (P = 0.02, n = 210) and showed a non-significant trend for better ΔDAS in a subset of RA patients with gene expression data (n = 31, etanercept-treated). A small, multi-ethnic replication showed a non-significant trend towards an association among etanercept-treated RA patients of Portuguese ancestry (n = 139, P = 0.4), but no association among patients of Japanese ancestry (n = 151, P = 0.8). Our study demonstrates that an allele associated with response to etanercept therapy is also associated with CD84 gene expression, and further that CD84 expression correlates with disease activity. These findings support a model in which CD84 genotypes and/or expression may serve as a useful biomarker for response to etanercept treatment in RA patients of European ancestry. © 2013 Cui et al