Dimitrios A. Pappas

Genetics of rheumatoid arthritis contributes to biology and drug discovery

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological data sets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA). Here we performed a genome-wide association study meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ∼10 million single-nucleotide polymorphisms. We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 101 (refs 2, 3, 4). We devised an in silico pipeline using established bioinformatics methods based on functional annotation, cis-acting expression quantitative trait loci and pathway analyses—as well as novel methods based on genetic overlap with human primary immunodeficiency, haematological cancer somatic mutations and knockout mouse phenotypes—to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.

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.2×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.

Which spinal lesions are associated with new bone formation in patients with ankylosing spondylitis treated with anti-TNF agents? A long-term observational study using MRI and conventional radiography

Objective To study the relationship of spinal inflammation and fatty degeneration (FD) as detected by MRI and new bone formation seen on conventional radiographs (CRs) in ankylosing spondylitis (AS). Methods CRs at baseline, 2 years and 5 years and spinal MRIs at baseline and 2 years of 73 AS patients treated with infliximab in European AS Infliximab Cohort were available. Relative risks (RR) were calculated with a general linear model after adjustment for within-patient variation. Results In a total of 1466 vertebral edges (VEs) without baseline syndesmophytes, 61 syndesmophytes developed at 5 years, the majority of which (57.4%) had no corresponding detectable MRI lesions at baseline. VEs with both inflammation and FD at baseline had the highest risk (RR 3.3, p=0.009) for syndesmophyte formation at 5 years, followed by VEs that developed new FD or did not resolve FD at 2 years (RR=2.3, p=0.034), while inflammation at baseline with no FD at 2 years had the lowest risk for syndesmophyte formation at 5 years (RR=0.8). Of the VEs with inflammation at baseline, >70% resolved completely, 28.8% turned into FD after 2 years, but only 1 syndesmophyte developed within 5 years. Conclusions Parallel occurrence of inflammation and FD at baseline and development of FD without prior inflammation after 2 years were significantly associated with syndesmophyte formation after 5 years of anti-tumour necrosis factor (TNF) therapy. However, the sequence ‘inflammation–FD–new bone formation’ was rarely observed, an argument against the TNF-brake hypothesis. Whether an early suppression of inflammation leads to a decrease of the risk for new bone formation remains to be demonstrated.

Association of fine specificity and repertoire expansion of anticitrullinated peptide antibodies with rheumatoid arthritis associated interstitial lung disease

Background Interstitial lung disease (ILD) is associated with high morbidity and mortality in rheumatoid arthritis (RA). Citrullinated proteins are observed in RA lung tissues; however, the association of specific anticitrullinated peptide antibodies (ACPA) with ILD in RA is unknown. Methods RA patients underwent multidetector CT (MDCT) of the chest, from which ILD features and a semiquantitative ILD Score (ILDS; range 0–32) were assessed. Anti-CCP (CCP2) and levels of a panel of antibodies against 17 citrullinated and four non-citrullinated peptides were assessed from concurrent serum samples using a custom Bio-Plex bead array. High level ACPA was defined as ≥the group 75th percentile. Results Among the 177 RA patients studied, median levels of CCP2 and all specific ACPAs were 46–273% higher among RA patients with versus those without ILD (all p values <0.05), and higher levels correlated with higher ILDS. In contrast, levels of non-citrullinated protein antibodies were not higher in those with ILD. RA patients had a median of 2 high level ACPA reactivities (range 0–16), with each high level ACPA associated, on average, with a 0.10 unit higher ILDS (p=0.001). This association remained significant after adjusting for characteristics associated with ILD (age, gender, current and former smoking, Disease Activity Score for 28 joints, current prednisone and leflunomide use). More high level ACPA were observed in those with versus without pulmonary function restriction or impaired diffusion. Conclusions Our findings of a broader ACPA repertoire in RA ILD suggest a possible role for ACPA in the pathogenesis of ILD.

High-density genotyping of immune loci in Koreans and Europeans identifies eight new rheumatoid arthritis risk loci.

OBJECTIVE A highly polygenic aetiology and high degree of allele-sharing between ancestries have been well elucidated in genetic studies of rheumatoid arthritis. Recently, the high-density genotyping array Immunochip for immune disease loci identified 14 new rheumatoid arthritis risk loci among individuals of European ancestry. Here, we aimed to identify new rheumatoid arthritis risk loci using Korean-specific Immunochip data. METHODS We analysed Korean rheumatoid arthritis case-control samples using the Immunochip and genome-wide association studies (GWAS) array to search for new risk alleles of rheumatoid arthritis with anticitrullinated peptide antibodies. To increase power, we performed a meta-analysis of Korean data with previously published European Immunochip and GWAS data for a total sample size of 9299 Korean and 45,790 European case-control samples. RESULTS We identified eight new rheumatoid arthritis susceptibility loci (TNFSF4, LBH, EOMES, ETS1-FLI1, COG6, RAD51B, UBASH3A and SYNGR1) that passed a genome-wide significance threshold (p<5×10(-8)), with evidence for three independent risk alleles at 1q25/TNFSF4. The risk alleles from the seven new loci except for the TNFSF4 locus (monomorphic in Koreans), together with risk alleles from previously established RA risk loci, exhibited a high correlation of effect sizes between ancestries. Further, we refined the number of single nucleotide polymorphisms (SNPs) that represent potentially causal variants through a trans-ethnic comparison of densely genotyped SNPs. CONCLUSIONS This study demonstrates the advantage of dense-mapping and trans-ancestral analysis for identification of potentially causal SNPs. In addition, our findings support the importance of T cells in the pathogenesis and the fact of frequent overlap of risk loci among diverse autoimmune diseases.

Myocardial citrullination in rheumatoid arthritis: a correlative histopathologic study

IntroductionThe aim of this study was to explore the presence and localization of myocardial citrullination in samples from rheumatoid arthritis (RA) patients compared to rheumatic and non-rheumatic disease control groups.MethodsArchived myocardial samples obtained during autopsy from 1995 to 2009 were assembled into four groups: RA; scleroderma; fatal myocarditis; and non-rheumatic disease controls. Samples were examined by immunohistochemistry (IHC) for the presence and localization of citrullination and peptidyl arginine deiminase enzymes (PADs) by a single cardiovascular pathologist blinded to disease group and clinical characteristics.ResultsMyocardial samples from seventeen RA patients were compared with those from fourteen controls, five fatal myocarditis patients, and ten scleroderma patients. Strong citrullination staining was detected exclusively in the myocardial interstitium in each of the groups. However, average and peak anti-citrulline staining was 59% and 44% higher, respectively, for the RA group compared to the combined non-RA groups (P < 0.05 for both comparisons). Myocardial fibrosis did not differ between the groups. In contrast to citrullination, PADs 1 to 3 and 6 were detected in cardiomyocytes (primarily PADs 1 and 3), resident inflammatory cells (primarily PADs 2 and 4), and, to a smaller extent, in endothelial cells and vascular smooth muscle cells. PAD staining did not co-localize with anti-citrulline staining in the interstitium and did not vary by disease state.ConclusionsStaining for citrullination was higher in the myocardial interstitium of RA compared to other disease states, a finding that could link autoimmunity to the known increase in myocardial dysfunction and heart failure in RA.

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.

A weighted genetic risk score using all known susceptibility variants to estimate rheumatoid arthritis risk

Background There is currently great interest in the incorporation of genetic susceptibility loci into screening models to identify individuals at high risk of disease. Here, we present the first risk prediction model including all 46 known genetic loci associated with rheumatoid arthritis (RA). Methods A weighted genetic risk score (wGRS) was created using 45 RA non-human leucocyte antigen (HLA) susceptibility loci, imputed amino acids at HLA-DRB1 (11, 71 and 74), HLA-DPB1 (position 9) HLA-B (position 9) and gender. The wGRS was tested in 11 366 RA cases and 15 489 healthy controls. The risk of developing RA was estimated using logistic regression by dividing the wGRS into quintiles. The ability of the wGRS to discriminate between cases and controls was assessed by receiver operator characteristic analysis and discrimination improvement tests. Results Individuals in the highest risk group showed significantly increased odds of developing anti-cyclic citrullinated peptide-positive RA compared to the lowest risk group (OR 27.13, 95% CI 23.70 to 31.05). The wGRS was validated in an independent cohort that showed similar results (area under the curve 0.78, OR 18.00, 95% CI 13.67 to 23.71). Comparison of the full wGRS with a wGRS in which HLA amino acids were replaced by a HLA tag single-nucleotide polymorphism showed a significant loss of sensitivity and specificity. Conclusions Our study suggests that in RA, even when using all known genetic susceptibility variants, prediction performance remains modest; while this is insufficiently accurate for general population screening, it may prove of more use in targeted studies. Our study has also highlighted the importance of including HLA variation in risk prediction models.

Herpes Zoster Reactivation in Patients with Rheumatoid Arthritis: Analysis of Disease Characteristics and Disease-Modifying Antirheumatic Drugs

To identify the rheumatoid arthritis (RA) characteristics associated with increased herpes zoster (HZ) risk in the Corrona registry RA patients, and to evaluate the risk in initiators of tumor necrosis factor inhibitors (TNFi) or non‐TNFi biologic agents or (among those who were currently on or had been previously treated with methotrexate [MTX]) conventional synthetic disease‐modifying antirheumatic drugs (csDMARDs) other than MTX.

Brief Report: Identification of BACH2 and RAD51B as Rheumatoid Arthritis Susceptibility Loci in a Meta-Analysis of Genome-Wide Data.

Objective A recent high-density fine-mapping (ImmunoChip) study of genetic associations in rheumatoid arthritis (RA) identified 14 risk loci with validated genome-wide significance, as well as a number of loci showing associations suggestive of significance (P = 5 × 10−5 < 5 × 10−8), but these have yet to be replicated. The aim of this study was to determine whether these potentially significant loci are involved in the pathogenesis of RA, and to explore whether any of the loci are associated with a specific RA serotype. Methods A total of 16 single-nucleotide polymorphisms (SNPs) were selected for genotyping and association analyses in 2 independent validation cohorts, comprising 6,106 RA cases and 4,290 controls. A meta-analysis of the data from the original ImmunoChip discovery cohort and from both validation cohorts was carried out, for a combined total of 17,581 RA cases and 20,160 controls. In addition, stratified analysis of patient subsets, defined according to their anti–cyclic citrullinated peptide (anti-CCP) antibody status, was performed. Results A significant association with RA risk (P < 0.05) was replicated for 6 of the SNPs assessed in the validation cohorts. All SNPs in the validation study had odds ratios (ORs) for RA susceptibility in the same direction as those in the ImmunoChip discovery study. One SNP, rs72928038, mapping to an intron of BACH2, achieved genome-wide significance in the meta-analysis (P = 1.2 × 10−8, OR 1.12), and a second SNP, rs911263, mapping to an intron of RAD51B, was significantly associated in the anti-CCP–positive RA subgroup (P = 4 × 10−8, OR 0.89), confirming that both are RA susceptibility loci. Conclusion This study provides robust evidence for an association of RA susceptibility with genes involved in B cell differentiation (BACH2) and DNA repair (RAD51B). The finding that the RAD51B gene exhibited different associations based on serologic subtype adds to the expanding knowledge base in defining subgroups of RA.