Search

ABSTRACTStudies investigating the immunogenicity of additional COVID-19 vaccine doses in immunosuppressed patients with inflammatory rheumatic diseases (IRD) are still limited. The objective was to explore the antibody response including response to omicron virus subvariants (sBA.1 and sBS.2) after third and fourth COVID-19 vaccine doses in Swedish IRD patients treated with immunomodulating drugs compared to controls. Antibody levels to spike wild-type antigens (full-length protein and S1) and the omicron variants sBA.1 and sBA.2 (full-length proteins) were measured. A positive response was defined as having antibody levels over cut-off or ≥fourfold increase in post-vaccination levels for both antigens. Patients with arthritis, vasculitis, and other autoimmune diseases (n = 414), and controls (n = 61) receiving biologic/targeted synthetic disease-modifying anti-rheumatic drugs (DMARDs) with or without conventional synthetic DMARDs participated. Of these, blood samples were available for 370 patients and 52 controls after three doses, and 65 patients and 15 controls after four doses. Treatment groups after three vaccine doses were rituximab (n = 133), abatacept (n = 22), IL6r inhibitors (n = 71), JAnus Kinase inhibitors (JAK-inhibitors) (n = 56), tumor necrosis factor inhibitor (TNF-inhibitors) (n = 61), IL12/23/17 inhibitors (n = 27), and controls (n = 52). The percentage of responders after three and four vaccine doses was lower in rituximab-treated patients (59% and 57%) compared to controls (100%) (P < 0.001). After three doses, the percentage of responders in all other groups was 100%, including response to omicron sBA.1 and sBA.2. In rituximab-treated patients, higher baseline immunoglobulin G (IgG) and longer time-period between rituximab and vaccination predicted better response. In this Swedish nationwide study including IRD patients three and four COVID-19 vaccine doses were immunogenic in patients treated with IL6r inhibitors, TNF-inhibitors, JAK-inhibitors, and IL12/23/17-inhibitors but not in rituximab. As >50% of rituximab patients responded to vaccines including omicron subvariants, these patients should be prioritized for additional vaccine doses.IMPORTANCEResults from this study provide further evidence that additional doses of COVID-19 vaccines are immunogenic and result in satisfactory antibody response in a majority of patients with inflammatory rheumatic diseases (IRD) receiving potent immunomodulating treatments such as biological or targeted disease-modifying anti-rheumatic drugs (DMARDs) given as monotherapy or combined with traditional DMARDs. We observed that rituximab treatment, both as monotherapy and combined with csDMARDs, impaired antibody response, and only roughly 50% of patients developed a satisfactory antibody response including response to omicron subvariants after the third vaccine. In addition, higher IgG levels at the last rituximab course before the third vaccine dose and a longer time after the last rituximab treatment increased the chance of a satisfactory antibody response. These results indicate that rituximab-treated patients should be prioritized for additional vaccine doses.CLINICAL TRIALSEudraCT (European Union Drug Regulating Authorities Clinical Trials Database) with number 2021-000880-63.

Abstract Background Our study aims to investigate an intrinsic link underlying sex hormone-binding globulin (SHBG) and rheumatoid arthritis (RA), which remains inconclusive in observational settings. Methods Summary statistics were collected from the largest GWAS(s) on SHBG adjusted for BMI (SHBGadjBMI; N overall = 368,929; N men = 180,094; N women = 188,908), crude SHBG (N overall = 370,125; N men = 180,726; N women = 189,473), and RA (N case = 22,350; N control = 74,823). A genome-wide cross-trait design was performed to quantify global and local genetic correlation, identify pleiotropic loci, and infer a causal relationship. Results Among the overall population, a significant global genetic correlation was observed for SHBGadjBMI and RA ( $$r_{{\text{g}}}$$ r g = 0.11, P = 1.0 × 10−4) which was further supported by local signal (1q25.2). A total of 18 independent pleiotropic SNPs were identified, of which three were highly likely causal variants and four were found to have effects on both traits through gene expression mediation. A putative causal association of SHBGadjBMI on RA was demonstrated (OR = 1.20, 95% CI = 1.01–1.43) without evidence of reverse causality (OR = 0.999, 95% CI = 0.997–1.000). Sex-specific analyses revealed distinct shared genetic regions (men: 1q32.1-q32.2 and 5p13.1; women: 1q25.2 and 22q11.21-q11.22) and diverse pleiotropic SNPs (16 in men and 18 in women, nearly half were sex-specific) underlying SHBGadjBMI and RA, demonstrating biological disparities between sexes. Replacing SHBGadjBMI with crude SHBG, a largely similar yet less significant pattern of results was observed. Conclusion Our cross-trait analysis suggests an intrinsic, as well as a sex-specific, link underlying SHBG and RA, providing novel insights into disease etiology.

Abstract Rheumatoid arthritis (RA) is an autoimmune disease characterized by systemic inflammation and is mediated by multiple immune cell types. In this work, we aimed to determine the relevance of changes in cell proportions in peripheral blood mononuclear cells (PBMCs) during the development of disease and following treatment. Samples from healthy blood donors, newly diagnosed RA patients, and established RA patients that had an inadequate response to MTX and were about to start tumor necrosis factor inhibitors (TNFi) treatment were collected before and after 3 months of treatment. We used in parallel a computational deconvolution approach based on RNA expression and flow cytometry to determine the relative cell-type frequencies. Cell-type frequencies from deconvolution of gene expression indicate that monocytes (both classical and non-classical) and CD4+ cells (Th1 and Th2) were increased in RA patients compared to controls, while NK cells and B cells (naïve and mature) were significantly decreased in RA patients. Treatment with MTX caused a decrease in B cells (memory and plasma cell), and a decrease in CD4 Th cells (Th1 and Th17), while treatment with TNFi resulted in a significant increase in the population of B cells. Characterization of the RNA expression patterns found that most of the differentially expressed genes in RA subjects after treatment can be explained by changes in cell frequencies (98% and 74% respectively for MTX and TNFi).

Abstract Genotypes causing pregnancy loss and perinatal mortality are depleted among living individuals and are therefore difficult to find. To explore genetic causes of recessive lethality, we searched for sequence variants with deficit of homozygosity among 1.52 million individuals from six European populations. In this study, we identified 25 genes harboring protein-altering sequence variants with a strong deficit of homozygosity (10% or less of predicted homozygotes). Sequence variants in 12 of the genes cause Mendelian disease under a recessive mode of inheritance, two under a dominant mode, but variants in the remaining 11 have not been reported to cause disease. Sequence variants with a strong deficit of homozygosity are over-represented among genes essential for growth of human cell lines and genes orthologous to mouse genes known to affect viability. The function of these genes gives insight into the genetics of intrauterine lethality. We also identified 1077 genes with homozygous predicted loss-of-function genotypes not previously described, bringing the total set of genes completely knocked out in humans to 4785.

Abstract Objectives To investigate the association between psychosocial vulnerability, defined as either low social support or low decision latitude at work, and disease remission at 3, 12, and 60 months in patients with rheumatoid arthritis (RA). Methods This cohort study included all patients enrolled in both the Swedish Epidemiological Investigation of Rheumatoid Arthritis (EIRA) 1996–2015 and the Swedish Rheumatology Quality Register (SRQ, n = 2820). Information on social support and decision latitude at work at RA diagnosis were identified from the EIRA questionnaire. Indexes for levels of social support and decision latitude at work, respectively, were calculated based on the questionnaire. Low social support and low decision latitude at work, respectively, were identified by a score in the lowest quartile and compared with the three other quartiles (not low). Disease-activity parameters were retrieved from SRQ at 3, 12, and 60 months. The associations between social support or decision latitude at work, respectively, and Disease Activity Score 28 joint count with C-reactive protein (DAS28-CRP) remission were analysed using logistic regression models adjusted for age, sex, smoking habits, alcohol habits, symptom duration, and educational level. Results Having low social support (n = 591) was not associated with DAS28-CRP remission at 3 (OR 0.93, 95% CI 0.74–1.16), 12 (OR 0.96, 95%CI 0.75–1.23), or 60 (OR 0.89, 95%CI 0.72–1.10) months compared to not low social support (n = 2209). No association was observed for low (n = 212) versus not low (n = 635) decision latitude at work and DAS28-CRP remission at 3 (OR 0.84, 95%CI 0.54–1.31), 12 (OR 0.81, 95%CI 0.56–1.16), or 60 (OR 1.37, 95%CI 0.94–2.01) months. Conclusion In a country with general access to healthcare, psychosocial vulnerability does not influence the likelihood of achieving remission in early RA.

Abstract Background Studies on the association between coffee, a modifiable lifestyle factor, and rheumatoid arthritis (RA), a chronic autoimmune disease primarily affecting the joints, have been conflicting. The aim of the present study was to study the association between coffee consumption and risk of RA in the context of different lifestyle factors. Methods We included 2184 cases (72% women, mean age 55 years) newly diagnosed with RA during 2005–2018 in Sweden and 4201 controls matched on age, sex, and residential area. Data on coffee consumption was collected through a food frequency questionnaire and categorized into < 2 (reference), 2–< 4, 4–< 6, and ≥ 6 cups/day. We calculated odds ratios (OR) with 95% confidence intervals (CI) for coffee consumption and risk of RA, in a crude model (taking matching factors into account), and then adjusted first for smoking and further for BMI, educational level, alcohol consumption, and physical activity. We also stratified analyses on sex, smoking, rheumatoid factor, and anti-CCP2 status. Results In the crude model, high coffee consumption was associated with increased risk of RA (OR = 1.50, 95% CI 1.20–1.88 for ≥ 6 cups/day compared to < 2 cups). After adjusting for smoking, the OR decreased and was no longer statistically significant (OR = 1.16, 95% CI 0.92–1.46) and decreased further in the full model (OR = 1.14 95% CI 0.89–1.45). This pattern held true in all strata. Conclusion The findings from this large, population-based case-control study did not support a significant association between coffee consumption and risk of RA as a whole nor within different subgroups.

Objective The study objective was to identify subgroups of patients with rheumatoid arthritis (RA) based on their health status 3 years after diagnosis and to assess potential associations to clinical presentation at diagnosis. Methods This observational study included patients with RA with 3‐year follow‐up data from the Swedish Epidemiological Investigation of RA study, collected from 2011 to 2018. Hierarchical agglomerative cluster analysis, based on symptoms of pain, fatigue, sleep quality, mood disturbances, and overall health‐related quality of life (HRQoL), was used to identify subgroups 3 years after diagnosis. Modified Poisson regression was used to estimate risk ratios (RRs) and 95% confidence intervals (CIs) for the associations between the subgroups and patient characteristics at diagnosis. Results A total of 1055 individuals constituted the study population, of whom 1011 had complete data on the clustering variables and were therefore eligible for analysis (73% women, median age 58 years). The following three clusters were identified: cluster 1 (466 patients with good health status), cluster 2 (398 patients in an intermediate group), and cluster 3 (147 patients with high levels of pain and fatigue together with markedly impaired HRQoL). Cluster 3 was associated to higher baseline pain (RR: 3.71 [95% CI: 2.14‐6.41]), global health (RR: 6.60 [95% CI: 3.53‐12.33]), and the Stanford Health Assessment Questionnaire (RR: 4.40 [95% CI: 2.46‐7.87]), compared with cluster 1 (highest compared with lowest quartiles). An inverse association was seen for baseline swollen joint count (RR: 0.51 [95% CI: 0.34‐0.85]). Conclusion A subgroup of patients with RA experience high levels of pain, fatigue, and psychosocial distress 3 years after diagnosis. This subgroup already displayed pronounced pain and functional disabilities at diagnosis.

BackgroundMethotrexate (MTX) is the first line treatment for rheumatoid arthritis (RA), but failure of satisfying treatment response occurs in a significant proportion of patients. Here we present a longitudinal multi-omics study aimed at detecting molecular and cellular processes in peripheral blood associated with a successful methotrexate treatment of rheumatoid arthritis.MethodsEighty newly diagnosed patients with RA underwent clinical assessment and donated blood before initiation of MTX, and 3 months into treatment. Flow cytometry was used to describe cell types and presence of activation markers in peripheral blood, the expression of 51 proteins was measured in serum or plasma, and RNA sequencing was performed in peripheral blood mononuclear cells (PBMC). Response to treatment after 3 months was determined using the EULAR response criteria. We assessed the changes in biological phenotypes during treatment, and whether these changes differed between responders and non-responders with regression analysis. By using measurements from baseline, we also tried to find biomarkers of future MTX response or, alternatively, to predict MTX response.ResultsAmong the MTX responders, (Good or Moderate according to EULAR treatment response classification, n = 60, 75%), we observed changes in 29 partly overlapping cell types proportions, levels of 13 proteins and expression of 38 genes during treatment. These changes were in most cases suppressions that were stronger among responders compared to non-responders. Within responders to treatment, we observed a suppression of FOXP3 gene expression, reduction of immunoglobulin gene expression and suppression of genes involved in cell proliferation. The proportion of many HLA-DR expressing T-cell populations were suppressed in all patients irrespective of clinical response, and the proportion of many IL21R+ T-cells were reduced exclusively in non-responders. Using only the baseline measurements we could not detect any biomarkers or prediction models that could predict response to MTX.ConclusionWe conclude that a deep molecular and cellular phenotyping of peripheral blood cells in RA patients treated with methotrexate can reveal previously not recognized differences between responders and non-responders during 3 months of treatment with MTX. This may contribute to the understanding of MTX mode of action and explain non-responsiveness to MTX therapy.

Abstract We aimed to search for common features in the autoreactive T cell receptor (TCR) repertoire in patients with rheumatoid arthritis (RA), focusing on the newly identified candidate antigen citrullinated Tenascin C (cit-TNC). Mononuclear cells from peripheral blood or synovial fluid of eight RA-patients positive for the RA-associated HLA-DRB1*04:01 allele were in-vitro cultured with recently identified citrullinated peptides from Tenascin C. Antigen-specific T cells were isolated using peptide-HLA tetramer staining and subsequently single-cell sequenced for paired alpha/beta TCR analyses by bioinformatic tools. TCRs were re-expressed for further studies of antigen-specificity and T cell responses. Autoreactive T cell lines could be grown out from both peripheral blood and synovial fluid. We demonstrate the feasibility of retrieving true autoreactive TCR sequences by validating antigen-specificity in T cell lines with re-expressed TCRs. One of the Tenascin C peptides, cit-TNC22, gave the most robust T cell responses including biased TCR gene usage patterns. The shared TCR-beta chain signature among the cit-TNC22-specific TCRs was evident in blood and synovial fluid of different patients. The identification of common elements in the autoreactive TCR repertoire gives promise to the possibility of both immune monitoring of the autoimmune components in RA and of future antigen- or TCR-targeted specific intervention in subsets of patients.

Abstract Background Tea is a popular beverage around the world and has properties that can affect the immune system. The association between tea consumption and the risk of rheumatoid arthritis (RA), a chronic autoimmune disease primarily affecting the joints, is not well studied and results are conflicting. Methods We collected data on tea consumption for 2237 incident RA cases diagnosed 2005–2018 and 4661 controls matched on age, sex, and residential area. Tea consumption was classified into no (0 cups/day), irregular (

Abstract Background Hormonal reproductive factors have been suggested to play an important role in the etiology of rheumatoid arthritis (RA), an autoimmune inflammatory disorder affecting primarily women. We conducted a two-sample Mendelian randomization (MR) study examining three relevant exposures, age at menarche (AAM), age at natural menopause (ANM), and age at first birth (AFB) with the risk of RA. Methods We collected summary statistics from the hitherto largest GWAS conducted in AAM (N = 329,345), ANM (N = 69,360), AFB (N = 251,151), and RA (N case = 14,361, N control = 43,923), all of European ancestry. We constructed strong instruments using hundreds of exposure-associated genetic variants and estimated causal relationship through different MR approaches including an inverse-variance weighted method, an MR-Egger regression and a weighted median method. We conducted a multivariable MR to control for pleiotropic effect acting in particular through obesity and socioeconomic status. We also performed important sensitivity analyses to verify model assumptions. Results We did not find any evidence in support for a causal association between genetically predicted reproductive factors and risk of RA (ORper-SD increment in AAM = 1.06 [0.98–1.15]; ORper-SD increment in ANM = 1.05 [0.98–1.11], OR per-SD increment in AFB = 0.85 [0.65–1.10]). Results remained consistent after removing palindromic SNPs (ORper-SD increment in AAM = 1.06 [0.97–1.15], ORper-SD increment in ANM = 1.05 [0.98–1.13], ORper-SD increment in AFB = 0.81 [0.61–1.07]) or excluding SNPs associated with potential confounding traits (ORper-SD increment in AAM = 1.03 [0.94–1.12], ORper-SD increment in ANM = 1.04 [0.95–1.14]). No outlying instrument was identified through the leave-one-out analysis. Conclusions Our MR study does not convincingly support a casual effect of reproductive factors, as reflected by age at menarche, age at menopause, and age at first birth, on the development of RA. Despite the largely augmented set of instruments we used, these instruments only explained a modest proportion of phenotypic variance of exposures. Our knowledge regarding this topic is still insufficient and future studies with larger sample size should be designed to replicate or dispute our findings.

Background: Why the adaptive immune system turns against citrullinated antigens in rheumatoid arthritis (RA) and whether anti-citrullinated protein antibodies (ACPAs) contribute to pathogenesis are questions that have triggered intense research, but still are not fully answered. Neutrophils may be crucial in this context, both as sources of citrullinated antigens and also as targets of ACPAs. To better understand how ACPAs and neutrophils contribute to RA, we studied the reactivity of a broad spectrum of RA patient-derived ACPA clones to activated or resting neutrophils, and we also compared neutrophil binding using polyclonal ACPAs from different patients. Methods: Neutrophils were activated by Ca2+ ionophore, PMA, nigericin, zymosan or IL-8, and ACPA binding was studied using flow cytometry and confocal microscopy. The roles of PAD2 and PAD4 were studied using PAD-deficient mice or the PAD4 inhibitor BMS-P5. Results: ACPAs broadly targeted NET-like structures, but did not bind to intact cells or influence NETosis. We observed high clonal diversity in ACPA binding to neutrophil-derived antigens. PAD2 was dispensable, but most ACPA clones required PAD4 for neutrophil binding. Using ACPA preparations from different patients, we observed high patient-to-patient variability in targeting neutrophil-derived antigens and similarly in another cellular effect of ACPAs, the stimulation of osteoclast differentiation. Conclusions: Neutrophils can be important sources of citrullinated antigens under conditions that lead to PAD4 activation, NETosis and the extrusion of intracellular material. A substantial clonal diversity in targeting neutrophils and a high variability among individuals in neutrophil binding and osteoclast stimulation suggest that ACPAs may influence RA-related symptoms with high patient-to-patient variability.

Abstract Background Fine-mapping of human leukocyte antigen (HLA) region for rheumatoid arthritis (RA) risk factors has identified several HLA alleles and its corresponding amino acid residues as independent signals (i.e., HLA-A, HLA-B, HLA-DPB1, and HLA-DQA1 genes), in addition to the well-established genetic factor in HLA-DRB1 gene. However, this was mainly performed in the Caucasian and East Asian populations, and data from different Asian regions is less represented. We aimed to evaluate whether there are independent RA risk variants in both anti-citrullinated protein antibody (ACPA)-positive and ACPA-negative RA patients from the multi-ethnic Malaysian population, using the fine-mapping of HLA region strategy. Methods We imputed the classical HLA alleles, amino acids, and haplotypes using the Immunochip genotyping data of 1260 RA cases (i.e., 530 Malays, 259 Chinese, 412 Indians, and 59 mixed ethnicities) and 1571 controls (i.e., 981 Malays, 205 Chinese, 297 Indians, and 87 mixed ethnicities) from the Malaysian Epidemiological Investigation of Rheumatoid Arthritis (MyEIRA) population-based case-control study. Stepwise logistic regression was performed to identify the independent genetic risk factors for RA within the HLA region. Results We confirmed that the HLA-DRB1 amino acid at position 11 with valine residue conferred the strongest risk effect for ACPA-positive RA (OR = 4.26, 95% CI = 3.30–5.49, P GWAS = 7.22 × 10−29) in the Malays. Our study also revealed that HLA-DRB1 amino acid at position 96 with histidine residue was negatively associated with the risk of developing ACPA-positive RA in the Indians (OR = 0.48, 95% CI = 0.37–0.62, P GWAS = 2.58 × 10−08). Interestingly, we observed that HLA-DQB1*03:02 allele was inversely related to the risk of developing ACPA-positive RA in the Malays (OR = 0.17, 95% CI = 0.09–0.30, P GWAS = 1.60 × 10−09). No association was observed between the HLA variants and risk of developing ACPA-negative RA in any of the three major ethnic groups in Malaysia. Conclusions Our results demonstrate that the RA-associated genetic factors in the multi-ethnic Malaysian population are similar to those in the Caucasian population, despite significant differences in the genetic architecture of HLA region across populations. A novel and distinct independent association between the HLA-DQB1*03:02 allele and ACPA-positive RA was observed in the Malays. In common with the Caucasian population, there is little risk from HLA region for ACPA-negative RA.

Based on the epidemiological link between periodontitis and rheumatoid arthritis (RA), and the unique feature of the periodontal bacterium Porphyromonas gingivalis to citrullinate proteins, it has been suggested that production of anti-citrullinated protein antibodies (ACPA), which are present in a majority of RA patients, may be triggered in the gum mucosa. To address this hypothesis, we investigated the antibody response to a citrullinated P. gingivalis peptide in relation to the autoimmune ACPA response in early RA, and examined citrulline-reactivity in monoclonal antibodies derived from human gingival B cells. Antibodies to a citrullinated peptide derived from P. gingivalis (denoted CPP3) and human citrullinated peptides were analyzed by multiplex array in 2,807 RA patients and 372 controls; associations with RA risk factors and clinical features were examined. B cells from inflamed gingival tissue were single-cell sorted, and immunoglobulin (Ig) genes were amplified, sequenced, cloned and expressed (n=63) as recombinant monoclonal antibodies, and assayed for citrulline-reactivities by enzyme-linked immunosorbent assay. Additionally, affinity-purified polyclonal anti-cyclic-citrullinated peptide (CCP2) IgG, and monoclonal antibodies derived from RA blood and synovial fluid B cells (n=175), were screened for CPP3-reactivity. Elevated anti-CPP3 antibody levels were detected in RA (11%), mainly CCP2+ RA, compared to controls (2%), p

IntroductionDigital diagnostic decision support tools promise to accelerate diagnosis and increase health care efficiency in rheumatology. Rheumatic? is an online tool developed by specialists in rheumatology and general medicine together with patients and patient organizations. It calculates a risk score for several rheumatic diseases. We ran a pilot study retrospectively testing Rheumatic? for its ability to differentiate symptoms from existing or emerging immune-mediated rheumatic diseases from other rheumatic and musculoskeletal complaints and disorders in patients visiting rheumatology clinics.Materials and MethodsThe performance of Rheumatic? was tested using in three university rheumatology centers: (A) patients at Risk for RA (Karolinska Institutet, n = 50 individuals with musculoskeletal complaints and anti-citrullinated protein antibody positivity) (B) patients with early joint swelling [dataset B (Erlangen) n = 52]. (C) Patients with early arthritis where the clinician considered it likely to be of auto-immune origin [dataset C (Leiden) n = 73]. In dataset A we tested whether Rheumatic? could predict the development of arthritis. In dataset B and C we tested whether Rheumatic? could predict the development of an immune-mediated rheumatic diseases. We examined the discriminative power of the total score with the Wilcoxon rank test and the area-under-the-receiver-operating-characteristic curve (AUC-ROC). Next, we calculated the test characteristics for these patients passing the first or second expert-based Rheumatic? scoring threshold.ResultsThe total test scores differentiated between: (A) Individuals developing arthritis or not, median 245 vs. 163, P < 0.0001, AUC-ROC = 75.3; (B) patients with an immune-mediated arthritic disease or not median 191 vs. 107, P < 0.0001, AUC-ROC = 79.0; but less patients with an immune-mediated arthritic disease or not amongst those where the clinician already considered an immune mediated disease most likely (median 262 vs. 212, P < 0.0001, AUC-ROC = 53.6). Threshold-1 (advising to visit primary care doctor) was highly specific in dataset A and B (0.72, 0.87, and 0.23, respectively) and sensitive (0.67, 0.61, and 0.67). Threshold-2 (advising to visit rheumatologic care) was very specific in all three centers but not very sensitive: specificity of 1.0, 0.96, and 0.91, sensitivity 0.05, 0.07, 0.14 in dataset A, B, and C, respectively.ConclusionRheumatic? is a web-based patient-centered multilingual diagnostic tool capable of differentiating immune-mediated rheumatic conditions from other musculoskeletal problems. The current scoring system needs to be further optimized.

BackgroundThe accumulation of risk for the development of rheumatoid arthritis (RA) is regarded as a continuum that may start with interacting environmental and genetic factors, proceed with the initiation of autoimmunity, and result in the formation of autoantibodies such as anti-citrullinated peptide antibodies (ACPA). In parallel, at-risk individuals may be asymptomatic or experience joint pain (arthralgia) that is itself non-specific or clinically suspicious for evolving RA, even in the absence of overt arthritis. Optimal strategies for the management of people at-risk of RA, both for symptom control and to delay or prevent progression to classifiable disease, remain poorly understood.MethodsTo help address this, groups of stakeholders from academia, clinical rheumatology, industry and patient research partners have collaborated to advance understanding, define and study different phases of the at-risk state. In this current report we describe different European initiatives in the field and the successful effort to build a European Registry of at-risk people to facilitate observational and interventional research.ResultsWe outline similarities and differences between cohorts of at-risk individuals at institutions spanning several countries, and how to best combine them within the new database. Over the past 2 years, besides building the technical infrastructure, we have agreed on a core set of variables that all partners should strive to collect for harmonization purposes.ConclusionWe emphasize to address this process from different angles and touch on the biologic, epidemiologic, analytic, and regulatory aspects of collaborative studies within a meta-database of people at-risk of RA.

HLA-DRB1 alleles have been associated with several autoimmune diseases. For anti-citrullinated protein antibody positive rheumatoid arthritis (RA), HLA-DRB1 shared epitope (SE) alleles are the major genetic risk factors. In order to study the genetic regulation of major histocompatibility complex (MHC) Class II gene expression in immune cells, we investigated transcriptomic profiles of a variety of immune cells from healthy individuals carrying different HLA-DRB1 alleles. Sequencing libraries from peripheral blood mononuclear cells, CD4+ T cells, CD8+ T cells, and CD14+ monocytes of 32 genetically pre-selected healthy female individuals were generated, sequenced and reads were aligned to the standard reference. For the MHC region, reads were mapped to available MHC reference haplotypes and AltHapAlignR was used to estimate gene expression. Using this method, HLA-DRB and HLA-DQ were found to be differentially expressed in different immune cells of healthy individuals as well as in whole blood samples of RA patients carrying HLA-DRB1 SE-positive versus SE-negative alleles. In contrast, no genes outside the MHC region were differentially expressed between individuals carrying HLA-DRB1 SE-positive and SE-negative alleles, thus HLA-DRB1 SE alleles have a strong cis effect on gene expression. Altogether, our findings suggest that immune effects associated with different allelic forms of HLA-DR and HLA-DQ may be associated not only with differences in the structure of these proteins, but also with differences in their expression levels.

Seropositive rheumatoid arthritis (RA) is characterized by the presence of rheumatoid factor (RF) and anti-citrullinated protein autoantibodies (ACPA) with different fine-specificities. Yet, other serum anti-modified protein autoantibodies (AMPA), e.g. anti-carbamylated (Carb), -acetylated (KAc), and malondialdehyde acetaldehyde (MAA) modified protein antibodies, have been described. In this comprehensive study, we analyze 30 different IgG and IgA AMPA reactivities to Cit, Carb, KAc, and MAA antigens detected by ELISA and autoantigen arrays in N=1985 newly diagnosed RA patients. Association with patient characteristics such as smoking and disease activity were explored. Carb and KAc reactivities by different assays were primarily seen in patients also positive for anti-citrulline reactivity. Modified vimentin (mod-Vim) peptides were used for direct comparison of different AMPA reactivities, revealing that IgA AMPA recognizing mod-Vim was mainly detected in subsets of patients with high IgG anti-Cit-Vim levels and a history of smoking. IgG reactivity to acetylation was mainly detected in a subset of patients with Cit and Carb reactivity. Anti-acetylated histone reactivity was RA-specific and associated with high anti-CCP2 IgG levels, multiple ACPA fine-specificities, and smoking status. This reactivity was also found to be present in CCP2+ RA-risk individuals without arthritis. Our data further demonstrate that IgG autoreactivity to MAA was increased in RA compared to controls with highest levels in CCP2+ RA, but was not RA-specific, and showed low correlation with other AMPA. Anti-MAA was instead associated with disease activity and was not significantly increased in CCP2+ individuals at risk of RA. Notably, RA patients could be subdivided into four different subsets based on their AMPA IgG and IgA reactivity profiles. Our serology results were complemented by screening of monoclonal antibodies derived from single B cells from RA patients for the same antigens as the RA cohort. Certain CCP2+ clones had Carb or Carb+KAc+ multireactivity, while such reactivities were not found in CCP2- clones. We conclude that autoantibodies exhibiting different patterns of ACPA fine-specificities as well as Carb and KAc reactivity are present in RA and may be derived from multireactive B-cell clones. Carb and KAc could be considered reactivities within the “Cit-umbrella” similar to ACPA fine-specificities, while MAA reactivity is distinctly different.

Recent advances in single-cell sequencing technologies enable the generation of large-scale data sets of paired TCR sequences from patients with autoimmune disease. Methods to validate and characterize patient-derived TCR data are needed, as well as relevant model systems that can support the development of antigen-specific tolerance inducing drugs. We have generated a pipeline to allow streamlined generation of ‘artificial’ T cells in a robust and reasonably high throughput manner for in vitro and in vivo studies of antigen-specific and patient-derived immune responses. Hereby chimeric (mouse-human) TCR alpha and beta constructs are re-expressed in three different formats for further studies: (i) transiently in HEK cells for peptide-HLA tetramer validation experiments, (ii) stably in the TCR-negative 58 ​T cell line for functional readouts such as IL-2 production and NFAT-signaling, and lastly (iii) in human HLA-transgenic mice for studies of autoimmune disease and therapeutic interventions. As a proof of concept, we have used human HLA-DRB1∗04:01 restricted TCR sequences specific for a type I diabetes-associated GAD peptide, and an influenza-derived HA peptide. We show that the same chimeric TCR constructs can be used in each of the described assays facilitating sequential validation and prioritization steps leading to humanized animal models.

Objective Rheumatoid arthritis (RA) is a chronic autoimmune inflammatory disorder with a global prevalence of approximately 0.5–1%. Patients with RA are at an increased risk of developing comorbidities (eg, cardiovascular disease, pulmonary disease, diabetes and depression). Despite this, there are limited recommendations for the management and implementation of associated comorbidities. This study aimed to identify good practice interventions in the care of RA and associated comorbidities.Methods A combination of primary research (180+ interviews with specialists across 12 European rheumatology centres) and secondary research (literature review of existing publications and guidelines/recommendations) were used to identify challenges in management and corresponding good practice interventions. Findings were prioritised and reviewed by a group of 18 rheumatology experts including rheumatologists, comorbidity experts, a patient representative and a highly specialised nurse.Results Challenges throughout the patient pathway (including delays in diagnosis and referral, shortage of rheumatologists, limited awareness of primary care professionals) and 18 good practice interventions were identified in the study. The expert group segmented and prioritised interventions according to three distinct stages of the disease: (1) suspected RA, (2) recent diagnosis of RA and (3) established RA. Examples of good practice interventions included enabling self-management (self-monitoring and disease management support, for example, lifestyle adaptations); early arthritis clinic; rapid access to care (online referral, triage, ultrasound-guided diagnosis); dedicated comorbidity specialists; enhanced communication with primary care (hotline, education sessions); and integrating patient registries into daily clinical practice.Conclusion Learning from implementation of good practice interventions in centres across Europe provides an opportunity to more widely improved care for patients with RA and associated comorbidities.

Abstract Background A relationship between rheumatoid arthritis (RA) and periodontitis has been suggested from findings that individuals with RA are prone to have advanced periodontitis and vice versa. In search of possible common pathogenetic features of these two diseases, we investigated the presence of citrullinated proteins and expression of endogenous peptidylarginine deiminases (PAD2 and PAD4), in periodontal tissue of individuals with periodontitis and healthy controls, in relation to the periodontal pathogens Porphyromonas gingivalis (P. gingivalis) and Aggregatibacter actinomycetemcomitans (A. actinomycetemcomitans), producing leukotoxin as virulence factor. These two oral bacteria have been suggested to be linked to anti-citrullinated protein antibodies in patients with RA. Methods Gingival tissue biopsies were obtained from 15 patients with periodontitis and 15 individuals without periodontal disease. Presence of CD3-positive lymphocytes, citrullinated proteins, PAD2, PAD4, P. gingivalis as well as A. actinomycetemcomitans and Mannheimia haemolytica produced leukotoxins were analysed by immunohistochemistry, followed by triple-blind semi-quantitative analysis. Mann–Whitney and Fisher’s exact tests were used to analyse differences between groups. PADI2 and PADI4 mRNA levels were assessed by RT-qPCR and analysed using Wilcoxon signed rank test. Results Increased staining of citrullinated proteins was observed in gingival connective tissue from subjects with periodontitis (80%, 12/15) compared to healthy gingival tissue (27%, 4/15), whereas no differences were observed in gingival epithelium. There was also an increased staining of the citrullinating enzymes PAD2 and PAD4 in gingival connective tissue of patients with periodontitis whereas similar levels of PAD2 and PAD4 were observed in the gingival epithelium of the two groups. Similarly, the mRNA levels of PADI2 and PADI4 were also increased in the gingival tissue of patients with periodontitis compared to healthy controls. Furthermore, presence of P. gingivalis and leukotoxins was comparable in both epithelium and connective tissue, from the different investigated individuals with and without periodontitis, and there were no correlations between the presence of periodontal pathogens and the expression of citrullinated proteins or PAD enzymes. Conclusion Chronic gingival inflammation is associated with increased local citrullination and PAD2 and PAD4 expression in periodontitis. The increased citrullination and PAD2 and PAD4 expression in periodontitis were, however, independent of the presence of periodontal pathogen P. gingivalis and A. actinomycetemcomitans leukotoxin.

Editorial summary Genomic medicine has considerable potential to provide novel diagnostic and therapeutic solutions for patients who have molecularly complex diseases and who are not responding to existing therapies. To bridge the gap between genomic medicine and clinical practice, integration of various data types, resources, and joint international initiatives will be required.

Abstract Background Methotrexate (MTX) is the standard first-line therapy in rheumatoid arthritis (RA) with variable clinical efficacy that is difficult to predict. The glycosylation status of immunoglobulin G (IgG) is altered in RA and influenced by MTX treatment. We aimed to further investigate if IgG glycosylation in untreated early RA can predict therapeutic response to MTX. Methods We used a shotgun proteomic approach to screen for the Fc glycopeptides in the serum of 12 control subjects and 59 untreated patients with early RA prior to and following MTX initiation. MTX treatment response was defined according to the European League Against Rheumatism at a median of 14 weeks (range 13–15) after treatment initiation. Seropositive patients were defined as those testing positive for anticitrullinated protein antibodies and/or rheumatoid factor at baseline (n = 44). Data analysis was performed using uni- and multivariate statistics. Results We could confirm a low abundance of galactosylated glycans in untreated patients with early RA compared with control subjects that was partially restored by MTX treatment. This was more evident among future nonresponders than among responders to MTX treatment. Results were further validated and confirmed by multivariate statistical analysis of the baseline Fc glycan, proteomic, and clinical data. We found that the ratio between the main agalactosylated (FA2) and main mono- and di-galactosylated Fc glycans (FA2G1 and FA2G2) of IgG1 ranked as the most prominent factor distinguishing responders from nonresponders. A low baseline ratio of FA2/[FA2G1 + FA2G2]-IgG1 was associated with nonresponse (OR 5.3 [1.6–17.0]) and was able to discriminate future nonresponders from responders to MTX therapy with a sensitivity of 70% (95% CI 46–88%) and a specificity of 69% (95% CI 52–83%). For seropositive patients (n = 44), this trend was improved with a sensitivity of 73% (95% CI 45–92%) for nonresponse and a specificity of 79% (95% CI 60–92%). Conclusions We show that the FA2/[FA2G1 + FA2G2] of IgG1 is a biomarker candidate that is significantly associated with nonresponding patients and has potential value for prediction of MTX clinical response.

Abstract Background Multiple factors, including interactions between genetic and environmental risks, are important in susceptibility to rheumatoid arthritis (RA). However, the underlying mechanism is not fully understood. This study was undertaken to evaluate whether DNA methylation can mediate the interaction between genotype and smoking in the development of anti-citrullinated peptide antibody (ACPA)-positive RA. Methods We investigated the gene-smoking interactions in DNA methylation using 393 individuals from the Epidemiological Investigation of Rheumatoid Arthritis (EIRA). The interaction between rs6933349 and smoking in the risk of developing ACPA-positive RA was further evaluated in a larger portion of the EIRA (1119 controls and 944 ACPA-positive patients with RA), and in the Malaysian Epidemiological Investigation of Rheumatoid Arthritis (MyEIRA) (1556 controls and 792 ACPA-positive patients with RA). Finally, mediation analysis was performed to investigate whether DNA methylation of cg21325723 mediates this gene-environment interaction on the risk of developing of ACPA-positive RA. Results We identified and replicated one significant gene-environment interaction between rs6933349 and smoking in DNA methylation of cg21325723. This gene-smoking interaction is a novel interaction in the risk of developing ACPA-positive in both Caucasian (multiplicative P value = 0.056; additive P value = 0.016) and Asian populations (multiplicative P value = 0.035; additive P value = 0.00027), and it is mediated through DNA methylation of cg21325723. Conclusions We showed that DNA methylation of cg21325723 can mediate the gene-environment interaction between rs6933349 and smoking, impacting the risk of developing ACPA-positive RA, thus being a potential regulator that integrates both internal genetic and external environmental risk factors.

Objectives The relationship between alcohol consumption and risk for rheumatoid arthritis (RA) is incompletely understood. We aimed to determine the influence of alcohol on anticitrullinated protein antibody (ACPA) positive and ACPA-negative RA and investigate potential interactions between alcohol consumption, smoking and the presence of human leucocyte antigen (HLA)-DRB1-shared epitope (SE).Methods A Swedish population-based case–control study with incident cases of RA was used (3353 cases, 2836 matched controls). Subjects with different HLA-DRB1-SE status, smoking and alcohol consumption were compared regarding risk of ACPA-positive and ACPA-negative RA, by calculating OR with 95% CI employing logistic regression. Interaction on the additive scale between alcohol, HLA-DRB1-SE and smoking was estimated by calculating the attributable proportion (AP) due to interaction.Results Compared with non-drinking, low and moderate alcohol consumption was dose dependently associated with a reduced risk of ACPA-positive and ACPA-negative RA. Independent of smoking habits, non-drinking and the presence of HLA-DRB1-SE interacted to increase the risk of ACPA-positive RA. Among HLA-DRB1-SE positive subjects, there was also a significant interaction between non-drinking and smoking with regard to risk for ACPA-positive RA. A three-way interaction was observed between alcohol, smoking and HLA-DRB1-SE with regard to risk for ACPA-positive RA (AP 0.7, 95% CI 0.6 to 0.8) that remained significant when the influence from the two-way interactions was removed (AP 0.4, 95% CI 0.2 to 0.6).Conclusions Our findings emphasize the need to investigate complex interactions between several environmental and genetic factors in order to better understand the etiology of RA. Whereas of great interest in an aetiological perspective, the finding of a protective role of alcohol on risk for RA must, however, be interpreted with caution in a clinical and public health perspective.

Objectives Previous studies on rheumatoid arthritis (RA) and body mass index (BMI) have yielded diverging results. We aimed to clarify the influence of BMI on the risk of developing anticitrullinated peptide antibody (ACPA)-positive and ACPA-negative RA by taking into consideration gender, smoking habits and human leukocyte antigen (HLA-DRB1) shared epitope (SE) status.Methods The present report is based on a Swedish population-based, case–control study with incident cases of RA (3572 cases, 5772 matched controls). Using logistic regression models, overweight/obese subjects were compared with normal weight subjects regarding risk of developing RA, by calculating ORs with 95% CIs.Results We observed diverging results for women and men. Among women, the risk of both ACPA-positive and ACPA-negative RA increased with increasing BMI, whereas an inverse association was observed among men for ACPA-positive RA. The results were similar regardless if RA onset before or after the age of 55 years was considered. When the analyses were stratified by smoking habits, the influence of BMI on RA risk was mainly restricted to smokers. Among women, a significant interaction was observed between smoking and overweight/obesity with regard to both subsets of RA. No interaction was observed between HLA-DRB1 SE and overweight/obesity with regard to RA risk.Conclusions The interaction between smoking and obesity regarding risk for RA in women warrants efforts to reduce these risk factors in those at risk for RA. The sex differences concerning the influence of obesity on RA risk merit further studies to verify these results and understand underlying mechanisms.

Rheumatoid arthritis (RA) associated anti-citrullinated protein autoantibodies (ACPA) target a wide range of modified proteins. Citrullination occurs during physiological processes such as apoptosis, yet little is known about the interaction of ACPA with nuclear antigens or apoptotic cells. Since uncleared apoptotic cells and neutrophil extracellular trap (NET) products have been postulated to be central sources of autoantigen and immunostimulation in autoimmune disease, we sought to characterize the anti-nuclear and anti-neutrophil reactivities of ACPA. Serology showed that a subset of anti-CCP2 seropositive RA patients had high reactivity to full-length citrullinated histones. In contrast, seronegative RA patients displayed elevated IgG reactivity to native histone compared to controls, but no citrulline-specific reactivity. Screening of 10 single B-cell derived monoclonal ACPA from RA patients revealed that four ACPA exhibited strong binding to apoptotic cells and three of these had anti-nuclear (ANA) autoantibody reactivity. Modified histones were confirmed to be the primary targets of this anti-nuclear ACPA subset following immunoprecipitation from apoptotic cell lysates. Monoclonal ACPA were also screened for reactivities against stimulated murine and human neutrophils, and all the nuclear-reactive monoclonal ACPA bound to NETs. Intriguingly, one ACPA mAb displayed a contrasting cytoplasmic perinuclear neutrophil binding and may represent a different NET-reactive ACPA subset. Notably, studies of CRISPR-Cas9 PAD4 KO cells and cells from PAD KO mice showed that the cytoplasmic NET-binding was fully dependent on PAD4, whilst nuclear- and histone-mediated NET reactivity was largely PAD-independent. Our further analysis revealed that the nuclear binding could be explained by consensus-motif driven ACPA cross-reactivity to acetylated histones. Specific acetylated histone peptides targeted by the monoclonal antibodies were identified and the anti-modified protein autoantibody (AMPA) profile of the ACPA was found to correlate with the functional activity of the antibodies. In conclusion, when investigating monoclonal ACPA, we could group ACPA into distinct subsets based on their nuclear binding-patterns and acetylation-mediated binding to apoptotic cells, neutrophils, and NETs. Differential anti-modified protein reactivities of RA-autoantibody subsets could have an important functional impact and provide insights in RA pathogenesis.

Abstract Background Twin studies are powerful models to elucidate epigenetic modifications resulting from gene–environment interactions. Yet, commonly a limited number of clinical twin samples are available, leading to an underpowered situation afflicted with false positives and hampered by low sensitivity. We investigated genome-wide DNA methylation data from two small sets of monozygotic twins representing different phases during the progression of rheumatoid arthritis (RA) to find novel genes for further research. Methods We implemented a robust statistical methodology aimed at investigating a small number of samples to identify differential methylation utilizing the comprehensive CHARM platform with whole blood cell DNA from two sets of twin pairs discordant either for ACPA (antibodies to citrullinated protein antigens)-positive RA versus ACPA-negative healthy or for ACPA-positive healthy (a pre-RA stage) versus ACPA-negative healthy. To deconvolute cell type-dependent differential methylation, we assayed the methylation patterns of sorted cells and used computational algorithms to resolve the relative contributions of different cell types and used them as covariates. Results To identify methylation biomarkers, five healthy twin pairs discordant for ACPAs were profiled, revealing a single differentially methylated region (DMR). Seven twin pairs discordant for ACPA-positive RA revealed six significant DMRs. After deconvolution of cell type proportions, profiling of the healthy ACPA discordant twin-set revealed 17 genome-wide significant DMRs. When methylation profiles of ACPA-positive RA twin pairs were adjusted for cell type, the analysis disclosed one significant DMR, associated with the EXOSC1 gene. Additionally, the results from our methodology suggest a temporal connection of the protocadherine beta-14 gene to ACPA-positivity with clinical RA. Conclusions Our biostatistical methodology, optimized for a low-sample twin design, revealed non-genetically linked genes associated with two distinct phases of RA. Functional evidence is still lacking but the results reinforce further study of epigenetic modifications influencing the progression of RA. Our study design and methodology may prove generally useful in twin studies.

Rheumatoid arthritis patients respond differently to anti-TNF treatment. Using community-based challenge, the authors show that currently available data does not reveal meaningful genetic predictors of response to anti-TNF therapy, thus confirming clinical observations.

Abstract OBJECTIVE: In rheumatoid arthritis (RA) several recent efforts have sought to discover means of predicting which patients would benefit from treatment. However, results have been discrepant with few successful replications. Our objective was to build a biobank with DNA, RNA and protein measurements to test the claim that the current state-of-the-art precision medicine will benefit RA patients. METHODS: We collected 451 blood samples from 61 healthy individuals and 185 RA patients initiating treatment, before treatment initiation and at a 3 month follow-up time. All samples were subjected to high-throughput RNA sequencing, DNA genotyping, extensive proteomics and flow cytometry measurements, as well as comprehensive clinical phenotyping. Literature review identified 2 proteins, 52 single-nucleotide polymorphisms (SNPs) and 72 gene-expression biomarkers that had previously been proposed as predictors of Tumor Necrosis Factor (TNF) inhibitor response (ΔDAS28-CRP), RESULTS: From these published TNFi biomarkers we found that 2 protein, 2 SNP and 8 mRNA biomarkers could be replicated in the 59 TNF initiating patients. Combining these replicated biomarkers into a single signature we found that we could explain 51% of the variation in ΔDAS28-CRP. This corresponds to a sensitivity of 0.73 and specificity of 0.78 for the prediction of three month ΔDAS28-CRP better than −1.2. CONCLUSIONS: The COMBINE biobank is currently the largest collection of multi-omics data from RA patients with high potential for discovery and replication. Taking advantage of this we surveyed the current state-of-the-art of drug-response stratification in RA, and identified a small set of previously published biomarkers available in peripheral blood which predicts clinical response to TNF blockade in this independent cohort.

BackgroundCollagen XI (CXI) is a heterotrimeric molecule with triple helical structure in which the α3(XI) chain is identical to the α1(II) chain of collagen II (CII), but with extensive posttranslational modifications. CXI molecules are intermingled in the cartilage collagen fibers, which are mainly composed of CII. One of the alpha chains in CXI is shared with CII and contains the immunodominant T cell epitope, but it is unclear whether there are shared B cell epitopes as the antibodies tend to recognize the triple helical structures.MethodsMice expressing the susceptible immune response gene Aq were immunized with CII or CXI. Serum antibody responses were measured, monoclonal antibodies were isolated and analyzed for specificity to CII, CXI, and triple helical collagen peptides using bead-based multiplex immunoassays, enzyme-linked immunosorbent assays, and Western blots. Arthritogenicity of the antibodies was investigated by passive transfer experiments.ResultsImmunization with CII or CXI leads to a strong T and B cell response, including a cross-reactive response to both collagen types. Immunization with CII leads to severe arthritis in mice, with a response toward CXI at the chronic stage, whereas CXI immunization induces very mild arthritis only. A series of monoclonal antibodies to CXI were isolated and of these, the L10D9 antibody bound to both CXI and CII equally strong, with a specific binding for the D3 epitope region of α3(XI) or α1(II) chain. The L10D9 antibody binds cartilage in vivo and induced severe arthritis. In contrast, the L5F3 antibody only showed weak binding and L7D8 antibody has no binding to cartilage and did not induce arthritis. The arthritogenic L10D9 antibody bound to an epitope shared with CII, the triple helical D3 epitope. Antibody levels to the shared D3 epitope were elevated in the sera from mice with arthritis as well as in rheumatoid arthritis.ConclusionCXI is immunologically not exposed in healthy cartilage but contains T and B cell epitopes cross-reactive with CII, which could be activated in both mouse and human arthritis and could evoke an arthritogenic response.

A large portion of common variant loci associated with genetic risk for schizophrenia reside within noncoding sequence of unknown function. Here, we demonstrate promoter and enhancer enrichment in schizophrenia variants associated with expression quantitative trait loci (eQTL). The enrichment is greater when functional annotations derived from the human brain are used relative to peripheral tissues. Regulatory trait concordance analysis ranked genes within schizophrenia genome-wide significant loci for a potential functional role, based on colocalization of a risk SNP, eQTL, and regulatory element sequence. We identified potential physical interactions of noncontiguous proximal and distal regulatory elements. This was verified in prefrontal cortex and -induced pluripotent stem cell–derived neurons for the L-type calcium channel (CACNA1C) risk locus. Our findings point to a functional link between schizophrenia-associated noncoding SNPs and 3D genome architecture associated with chromosomal loopings and transcriptional regulation in the brain.

INTRODUCTION:The possible hypothesis of a link between periodontitis and rheumatoid arthritis (RA), specifically anti-citrullinated protein antibody (ACPA) positive RA, prompted us to investigate the prevalence of periodontitis in the Swedish Epidemiological Investigation of RA (EIRA), a well-characterised population-based RA case-control cohort. METHODS:Periodontal status of 2,740 RA cases and 3,942 matched controls was retrieved through linking EIRA with the National Dental Health Registry (DHR), where dental diagnostic- and treatment codes on the adult Swedish population have been registered. Dental records from 100 cases and controls were reviewed to validate the periodontal diagnostic codes in DHR. RESULTS:The reviewed dental records confirmed 90% of the periodontitis diagnoses in DHR among RA cases, and 88% among controls. We found the positive predictive value of periodontitis diagnoses in the DHR to be 89% (95% CI 78 to 95%) with a sensitivity of 77% (95% CI: 65 to 86%). In total, 86% of EIRA participants were identified in DHR. The risk for periodontitis increased by age and current smoking status in both cases as well as controls. No significant differences in prevalence of periodontal disease in terms of gingivitis, periodontitis, peri-implantitis or increased risk for periodontitis or peri-implantitis were observed between RA cases and controls. In addition, there was no difference on the basis of seropositivity, ACPA or rheumatoid factor (RF), among patients with RA. CONCLUSIONS:Our data verify that smoking and ageing are risk factors for periodontitis, both in RA and controls. We found no evidence of an increased prevalence of periodontitis in patients with established RA compared to healthy controls, and no differences based on ACPA or RF status among RA subjects.

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

The Fc-glycan profile of IgG1 anti-citrullinated peptide antibodies (ACPA) in rheumatoid arthritis (RA) patients has recently been reported to be different from non-ACPA IgG1, a phenomenon which likely plays a role in RA pathogenesis. Herein we investigate the Fc-glycosylation pattern of all ACPA-IgG isotypes and simultaneously investigate in detail the IgG protein-chain sequence repertoire. IgG from serum or plasma (S/P, n = 14) and synovial fluid (SF, n = 4) from 18 ACPA-positive RA-patients was enriched using Protein G columns followed by ACPA-purification on cyclic citrullinated peptide-2 (CCP2)-coupled columns. Paired ACPA (anti-CCP2 eluted IgG) and IgG flow through (FT) fractions were analyzed by LC-MS/MS-proteomics. IgG peptides, isotypes and corresponding Fc-glycopeptides were quantified and interrogated using uni- and multivariate statistics. The Fc-glycans from the IgG4 peptide EEQFNSTYR was validated using protein A column purification. Relative to FT-IgG4, the ACPA-IgG4 Fc-glycan-profile contained lower amounts (p = 0.002) of the agalacto and asialylated core-fucosylated biantennary form (FA2) and higher content (p = 0.001) of sialylated glycans. Novel differences in the Fc-glycan-profile of ACPA-IgG1 compared to FT-IgG1 were observed in the distribution of bisected forms (n = 5, p = 0.0001, decrease) and mono-antennnary forms (n = 3, p = 0.02, increase). Our study also confirmed higher abundance of FA2 (p = 0.002) and lower abundance of afucosylated forms (n = 4, p = 0.001) in ACPA-IgG1 relative to FT-IgG1 as well as lower content of IgG2 (p = 0.0000001) and elevated content of IgG4 (p = 0.004) in ACPA compared to FT. One λ-variable peptide sequence was significantly increased in ACPA (p = 0.0001). In conclusion, the Fc-glycan profile of both ACPA-IgG1 and ACPA-IgG4 are distinct. Given that IgG1 and IgG4 have different Fc-receptor and complement binding affinities, this phenomenon likely affects ACPA effector- and immune-regulatory functions in an IgG isotype-specific manner. These findings further highlight the importance of antibody characterization in relation to functional in vivo and in vitro studies.

Translational medicine is becoming increasingly dependent upon data generated from health care, clinical research, and molecular investigations. This increasing rate of production and diversity in data has brought about several challenges, including the need to integrate fragmented databases, enable secondary use of patient clinical data from health care in clinical research, and to create information systems that clinicians and biomedical researchers can readily use. Our case study effectively integrates requirements from the clinical and biomedical researcher perspectives in a translational medicine setting. Our three principal achievements are (a) a design of a user-friendly web-based system for management and integration of clinical and molecular databases, while adhering to proper de-identification and security measures; (b) providing a real-world test of the system functionalities using clinical cohorts; and (c) system integration with a clinical decision support system to demonstrate system interoperability. We engaged two active clinical cohorts, 747 psoriasis patients and 2001 rheumatoid arthritis patients, to demonstrate efficient query possibilities across the data sources, enable cohort stratification, extract variation in antibody patterns, study biomarker predictors of treatment response in RA patients, and to explore metabolic profiles of psoriasis patients. Finally, we demonstrated system interoperability by enabling integration with an established clinical decision support system in health care. To assure the usefulness and usability of the system, we followed two approaches. First, we created a graphical user interface supporting all user interactions. Secondly we carried out a system performance evaluation study where we measured the average response time in seconds for active users, http errors, and kilobits per second received and sent. The maximum response time was found to be 0.12 seconds; no server or client errors of any kind were detected. In conclusion, the system can readily be used by clinicians and biomedical researchers in a translational medicine setting.

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.

ObjectivesGenome-wide association studies have facilitated the identification of over 30 susceptibility loci for rheumatoid arthritis (RA). However, evidence for a number of potential susceptibility genes have not so far reached genome-wide significance in studies of Caucasian RA.MethodsA cohort of 4286 RA patients from across Europe and 5642 population matched controls were genotyped for 25 SNPs, then combined in a meta-analysis with previously published data.ResultsSignificant evidence of association was detected for nine SNPs within the European samples. When meta-analysed with previously published data, 21 SNPs were associated with RA susceptibility. Although SNPs in the PTPN2 gene were previously reported to be associated with RA in both Japanese and European populations, we show genome-wide evidence for a different SNP within this gene associated with RA susceptibility in an independent European population (rs7234029, P = 4.4×10(-9)).ConclusionsThis study provides further genome-wide evidence for the association of the PTPN2 locus (encoding the T cell protein tyrosine phosphastase) with Caucasian RA susceptibility. This finding adds to the growing evidence for PTPN2 being a pan-autoimmune susceptibility gene.

Recent work has shown that much of the missing heritability of complex traits can be resolved by estimates of heritability explained by all genotyped SNPs. However, it is currently unknown how much heritability is missing due to poor tagging or additional causal variants at known GWAS loci. Here, we use variance components to quantify the heritability explained by all SNPs at known GWAS loci in nine diseases from WTCCC1 and WTCCC2. After accounting for expectation, we observed all SNPs at known GWAS loci to explain 1.29 x more heritability than GWAS-associated SNPs on average (P=3.3 x 10⁻⁵). For some diseases, this increase was individually significant: 2.07 x for Multiple Sclerosis (MS) (P=6.5 x 10⁻⁹) and 1.48 x for Crohn's Disease (CD) (P = 1.3 x 10⁻³); all analyses of autoimmune diseases excluded the well-studied MHC region. Additionally, we found that GWAS loci from other related traits also explained significant heritability. The union of all autoimmune disease loci explained 7.15 x more MS heritability than known MS SNPs (P < 1.0 x 10⁻¹⁶ and 2.20 x more CD heritability than known CD SNPs (P = 6.1 x 10⁻⁹), with an analogous increase for all autoimmune diseases analyzed. We also observed significant increases in an analysis of > 20,000 Rheumatoid Arthritis (RA) samples typed on ImmunoChip, with 2.37 x more heritability from all SNPs at GWAS loci (P = 2.3 x 10⁻⁶) and 5.33 x more heritability from all autoimmune disease loci (P < 1 x 10⁻¹⁶ compared to known RA SNPs (including those identified in this cohort). Our methods adjust for LD between SNPs, which can bias standard estimates of heritability from SNPs even if all causal variants are typed. By comparing adjusted estimates, we hypothesize that the genome-wide distribution of causal variants is enriched for low-frequency alleles, but that causal variants at known GWAS loci are skewed towards common alleles. These findings have important ramifications for fine-mapping study design and our understanding of complex disease architecture.

Nature Communications 7: Article number: 12460 (2016); Published: 23 August 2016; Updated: 10 October 2016. The HTML version of this Article incorrectly duplicated the authors S. Louis Bridges, Lindsey Criswell, Larry Moreland, Lars Klareskog, Saedis Saevarsdottir, Leonid Padyukov, Peter K. Gregersen, Stephen Friend, Robert Plenge, Gustavo Stolovitzky, Baldo Oliva, Yuanfang Guan and Lara M.

Genome-wide association (GWA) studies have identified numerous, replicable, genetic associations between common single nucleotide polymorphisms (SNPs) and risk of common autoimmune and inflammatory (immune-mediated) diseases, some of which are shared between two diseases. Along with epidemiological and clinical evidence, this suggests that some genetic risk factors may be shared across diseases-as is the case with alleles in the Major Histocompatibility Locus. In this work we evaluate the extent of this sharing for 107 immune disease-risk SNPs in seven diseases: celiac disease, Crohn's disease, multiple sclerosis, psoriasis, rheumatoid arthritis, systemic lupus erythematosus, and type 1 diabetes. We have developed a novel statistic for Cross Phenotype Meta-Analysis (CPMA) which detects association of a SNP to multiple, but not necessarily all, phenotypes. With it, we find evidence that 47/107 (44%) immune-mediated disease risk SNPs are associated to multiple-but not all-immune-mediated diseases (SNP-wise P(CPMA)

Non-Hodgkin lymphoma (NHL) represents a diverse group of hematological malignancies, of which follicular lymphoma (FL) is a prevalent subtype. A previous genome-wide association study has established a marker, rs10484561 in the human leukocyte antigen (HLA) class II region on 6p21.32 associated with increased FL risk. Here, in a three-stage genome-wide association study, starting with a genome-wide scan of 379 FL cases and 791 controls followed by validation in 1,049 cases and 5,790 controls, we identified a second independent FL-associated locus on 6p21.32, rs2647012 (OR(combined) = 0.64, P(combined) = 2 × 10(-21)) located 962 bp away from rs10484561 (r(2)

Epidemiology and candidate gene studies indicate a shared genetic basis for celiac disease (CD) and rheumatoid arthritis (RA), but the extent of this sharing has not been systematically explored. Previous studies demonstrate that 6 of the established non-HLA CD and RA risk loci (out of 26 loci for each disease) are shared between both diseases. We hypothesized that there are additional shared risk alleles and that combining genome-wide association study (GWAS) data from each disease would increase power to identify these shared risk alleles. We performed a meta-analysis of two published GWAS on CD (4,533 cases and 10,750 controls) and RA (5,539 cases and 17,231 controls). After genotyping the top associated SNPs in 2,169 CD cases and 2,255 controls, and 2,845 RA cases and 4,944 controls, 8 additional SNPs demonstrated P

BACKGROUND: To investigate the associations between HLA-DRB1 shared epitope (SE) alleles and rheumatoid arthritis in subsets of rheumatoid arthritis defined by autoantibodies in three Asian populations from Malaysia. METHODS: 1,079 rheumatoid arthritis patients and 1,470 healthy controls were included in the study. Levels of antibodies to citrullinated proteins (ACPA) and rheumatoid factors were assessed and the PCR-SSO method was used for HLA-DRB1 genotyping. RESULTS: The proportion of ACPA positivity among Malay, Chinese and Indian rheumatoid arthritis patients were 62.9%, 65.2% and 68.6%, respectively. An increased frequency of SE alleles was observed in ACPA-positive rheumatoid arthritis among the three Asian ethnic groups. HLA-DRB1*10 was highly associated with rheumatoid arthritis susceptibility in these Asian populations. HLA-DRB1*0405 was significantly associated with susceptibility to rheumatoid arthritis in Malays and Chinese, but not in Indians. HLA-DRB1*01 did not show any independent effect as a risk factor for rheumatoid arthritis in this study and HLA-DRB1*1202 was protective in Malays and Chinese. There was no association between SE alleles and ACPA- negative rheumatoid arthritis in any of the three Asian ethnic groups. CONCLUSION: The HLA-DRB1 SE alleles increase the risk of ACPA-positive rheumatoid arthritis in all three Asian populations from Malaysia.

OBJECTIVE: Glutathione-S-transferases (GSTs) play an important role in tobacco smoke detoxification, interestingly approximately 50% of individuals in most human populations lack the gene GSTM1 due to copy number variation (CNV). We aimed to investigate GSTM1 CNV in Rheumatoid Arthritis (RA) in relation to smoking and HLA-DRB1 shared epitope; the two best known risk factors for RA and in addition, to perform subanalyses in patients where relations between variations in GSTM1 and RA have previously been described. METHODS: qPCR was performed using TaqMan Copy Number assays (Applied Biosystems) for 2426 incident RA cases and 1257 controls from the Swedish EIRA. Odds ratio (OR) together with 95% confidence intervals (CI) was calculated and used as a measure of the relative risk of developing RA. RESULTS: No association between RA and GSTM1 CNV was observed when analyzing whole EIRA. However, ≥1 copy of GSTM1 appears to be a significant risk factor for autoantibody positive RA in non-smoking females ≥60 years (OR: 2.00 95% CI: 1.07-3.74), a population where such relationships have previously been described. Our data further suggest a protective effect of GSTM1 in ACPA-negative smoking men (OR: 0.56 95% CI: 0.35-0.90). CONCLUSION: We assessed the exact number of GSTM1 gene copies in relation to development and severity of RA. Our data provide support for the notion that variations in copy numbers of GSTM1 may influence risk in certain subsets of RA, but do not support a role for GSTM1 CNV as a factor that more generally modifies the influence of smoking on RA.

Patterns of genetic diversity have previously been shown to mirror geography on a global scale and within continents and individual countries. Using genome-wide SNP data on 5174 Swedes with extensive geographical coverage, we analyzed the genetic structure of the Swedish population. We observed strong differences between the far northern counties and the remaining counties. The population of Dalarna county, in north middle Sweden, which borders southern Norway, also appears to differ markedly from other counties, possibly due to this county having more individuals with remote Finnish or Norwegian ancestry than other counties. An analysis of genetic differentiation (based on pairwise F(st)) indicated that the population of Sweden's southernmost counties are genetically closer to the HapMap CEU samples of Northern European ancestry than to the populations of Sweden's northernmost counties. In a comparison of extended homozygous segments, we detected a clear divide between southern and northern Sweden with small differences between the southern counties and considerably more segments in northern Sweden. Both the increased degree of homozygosity in the north and the large genetic differences between the south and the north may have arisen due to a small population in the north and the vast geographical distances between towns and villages in the north, in contrast to the more densely settled southern parts of Sweden. Our findings have implications for future genome-wide association studies (GWAS) with respect to the matching of cases and controls and the need for within-county matching. We have shown that genetic differences within a single country may be substantial, even when viewed on a European scale. Thus, population stratification needs to be accounted for, even within a country like Sweden, which is often perceived to be relatively homogenous and a favourable resource for genetic mapping, otherwise inferences based on genetic data may lead to false conclusions.