Seven-chain adaptive immune receptor repertoire analysis in rheumatoid arthritis reveals novel features associated with disease and clinically relevant phenotypes.

BACKGROUND: In rheumatoid arthritis (RA), the activation of T and B cell clones specific for self-antigens leads to the chronic inflammation of the synovium. Here, we perform an in-depth quantitative analysis of the seven chains that comprise the adaptive immune receptor repertoire (AIRR) in RA. RESULTS: In comparison to controls, we show that RA patients have multiple and strong differences in the B cell receptor repertoire including reduced diversity as well as altered isotype, chain, and segment frequencies. We demonstrate that therapeutic tumor necrosis factor inhibition partially restores this alteration but find a profound difference in the underlying biochemical reactivities between responders and non-responders. Combining the AIRR with HLA typing, we identify the specific T cell receptor repertoire associated with disease risk variants. Integrating these features, we further develop a molecular classifier that shows the utility of the AIRR as a diagnostic tool. CONCLUSIONS: Simultaneous sequencing of the seven chains of the human AIRR reveals novel features associated with the disease and clinically relevant phenotypes, including response to therapy. These findings show the unique potential of AIRR to address precision medicine in immune-related diseases.

Prevalence and clinical characteristics of patients with rheumatoid arthritis with interstitial lung disease using unstructured healthcare data and machine learning.

OBJECTIVES: Real-world data regarding rheumatoid arthritis (RA) and its association with interstitial lung disease (ILD) is still scarce. This study aimed to estimate the prevalence of RA and ILD in patients with RA (RAILD) in Spain, and to compare clinical characteristics of patients with RA with and without ILD using natural language processing (NLP) on electronic health records (EHR). METHODS: Observational case-control, retrospective and multicentre study based on the secondary use of unstructured clinical data from patients with adult RA and RAILD from nine hospitals between 2014 and 2019. NLP was used to extract unstructured clinical information from EHR and standardise it into a SNOMED-CT terminology. Prevalence of RA and RAILD were calculated, and a descriptive analysis was performed. Characteristics between patients with RAILD and RA patients without ILD (RAnonILD) were compared. RESULTS: From a source population of 3 176 165 patients and 64 241 683 EHRs, 13 958 patients with RA were identified. Of those, 5.1% patients additionally had ILD (RAILD). The overall age-adjusted prevalence of RA and RAILD were 0.53% and 0.02%, respectively. The most common ILD subtype was usual interstitial pneumonia (29.3%). When comparing RAILD versus RAnonILD patients, RAILD patients were older and had more comorbidities, notably concerning infections (33.6% vs 16.5%, p<0.001), malignancies (15.9% vs 8.5%, p<0.001) and cardiovascular disease (25.8% vs 13.9%, p<0.001) than RAnonILD. RAILD patients also had higher inflammatory burden reflected in more pharmacological prescriptions and higher inflammatory parameters and presented a higher in-hospital mortality with a higher risk of death (HR 2.32; 95% CI 1.59 to 2.81, p<0.001). CONCLUSIONS: We found an estimated age-adjusted prevalence of RA and RAILD by analysing real-world data through NLP. RAILD patients were more vulnerable at the time of inclusion with higher comorbidity and inflammatory burden than RAnonILD, which correlated with higher mortality.

Longitudinal analysis of blood DNA methylation identifies mechanisms of response to tumor necrosis factor inhibitor therapy in rheumatoid arthritis.

BACKGROUND: Rheumatoid arthritis (RA) is a chronic, immune-mediated inflammatory disease of the joints that has been associated with variation in the peripheral blood methylome. In this study, we aim to identify epigenetic variation that is associated with the response to tumor necrosis factor inhibitor (TNFi) therapy. METHODS: Peripheral blood genome-wide DNA methylation profiles were analyzed in a discovery cohort of 62 RA patients at baseline and at week 12 of TNFi therapy. DNA methylation of individual CpG sites and enrichment of biological pathways were evaluated for their association with drug response. Using a novel cell deconvolution approach, altered DNA methylation associated with TNFi response was also tested in the six main immune cell types in blood. Validation of the results was performed in an independent longitudinal cohort of 60 RA patients. FINDINGS: Treatment with TNFi was associated with significant longitudinal peripheral blood methylation changes in biological pathways related to RA (FDR<0.05). 139 biological functions were modified by therapy, with methylation levels changing systematically towards a signature similar to that of healthy controls. Differences in the methylation profile of T cell activation and differentiation, GTPase-mediated signaling, and actin filament organization pathways were associated with the clinical response to therapy. Cell type deconvolution analysis identified CpG sites in CD4+T, NK, neutrophils and monocytes that were significantly associated with the response to TNFi. INTERPRETATION: Our results show that treatment with TNFi restores homeostatic blood methylation in RA. The clinical response to TNFi is associated to methylation variation in specific biological pathways, and it involves cells from both the innate and adaptive immune systems. FUNDING: The Instituto de Salud Carlos III.

Interactions between rheumatoid arthritis antibodies are associated with the response to anti-tumor necrosis factor therapy.

BACKGROUND: Blocking of the Tumor Necrosis Factor (TNF) activity is a successful therapeutic approach for 50-60% of rheumatoid arthritis (RA) patients. However, there are yet no biomarkers to stratify patients for anti-TNF therapy. Rheumatoid factor (RF) and anti-cyclic-citrullinated antibodies (anti-CCP) have been evaluated as biomarkers of response but the results have shown limited consistency. Anti-carbamylated protein (anti-CarP) and anti-peptidylarginine deiminase type 4 (anti-PAD4) antibodies have been much less studied. Despite being linked to common immune processes, the interaction between these markers has not been evaluated yet. Our aim was to analyze the interaction between these four antibodies in relation to the response to anti-TNF therapy. METHODS: For this objective, a prospective cohort of n = 80 RA patients starting anti-TNF therapy was recruited. Serum determinations at baseline were performed for RF, anti-CCP, anti-CarP and anti-PAD4 antibodies using enzyme-linked immunosorbent assays (ELISA). The clinical response to anti-TNF therapy was determined at week 12 using the change in DAS28 score. Association was performed using multivariate linear regression adjusting for baseline DAS28, sex and age. RESULTS: The interaction between pairs of antibodies was tested by the addition of an interaction term. We found two highly significant antibody interactions associated with treatment response: anti-CarP with anti-PAD4 (p = 0.0062), and anti-CCP with RF (p = 0.00068). The latter antibody interaction was replicated in an independent retrospective cohort of RA patients (n = 199, p = 0.04). CONCLUSIONS: The results of this study suggest that antibody interaction effects are important factors in the response to anti-TNF therapy in RA.

A Combined Transcriptomic and Genomic Analysis Identifies a Gene Signature Associated With the Response to Anti-TNF Therapy in Rheumatoid Arthritis.

Background: Rheumatoid arthritis (RA) is the most frequent autoimmune disease involving the joints. Although anti-TNF therapies have proven effective in the management of RA, approximately one third of patients do not show a significant clinical response. The objective of this study was to identify new genetic variation associated with the clinical response to anti-TNF therapy in RA. Methods: We performed a sequential multi-omic analysis integrating different sources of molecular information. First, we extracted the RNA from synovial biopsies of 11 RA patients starting anti-TNF therapy to identify gene coexpression modules (GCMs) in the RA synovium. Second, we analyzed the transcriptomic association between each GCM and the clinical response to anti-TNF therapy. The clinical response was determined at week 14 using the EULAR criteria. Third, we analyzed the association between the GCMs and anti-TNF response at the genetic level. For this objective, we used genome-wide data from a cohort of 348 anti-TNF treated patients from Spain. The GCMs that were significantly associated with the anti-TNF response were then tested for validation in an independent cohort of 2,706 anti-TNF treated patients. Finally, the functional implication of the validated GCMs was evaluated via pathway and cell type epigenetic enrichment analyses. Results: A total of 149 GCMs were identified in the RA synovium. From these, 13 GCMs were found to be significantly associated with anti-TNF response (P < 0.05). At the genetic level, we detected two of the 13 GCMs to be significantly associated with the response to adalimumab (P = 0.0015) and infliximab (P = 0.021) in the Spain cohort. Using the independent cohort of RA patients, we replicated the association of the GCM associated with the response to adalimumab (P = 0.0019). The validated module was found to be significantly enriched for genes involved in the nucleotide metabolism (P = 2.41e-5) and epigenetic marks from immune cells, including CD4+ regulatory T cells (P = 0.041). Conclusions: These findings show the existence of a drug-specific genetic basis for anti-TNF response, thereby supporting treatment stratification in the search for response biomarkers in RA.

Genetic association between CD96 locus and immunogenicity to anti-TNF therapy in Crohn's disease.

The production of antibodies to anti-tumor necrosis factor alpha (TNF) agents is one of the main causes of treatment failure in Crohn's disease (CD). To date, however, the contribution of genetics to anti-TNF immunogenicity in CD is still unknown. The objective of the present study was to identify genetic variation associated with anti-TNF immunogenicity in CD. We performed a two-stage genome-wide association study in a cohort of 96 and 123 adalimumab-treated patients, respectively. In the discovery stage, we identified a genome-wide significant association between the CD96 locus and the production of antibodies to anti-TNF treatment (P = 1.88e-09). This association was validated in the replication stage (P < 0.05). The risk allele for anti-TNF immunogenicity was found to be also associated with a lack of response to anti-TNF therapy (P = 0.019). These findings represent an important step toward the understanding of the immunogenicity-based mechanisms that underlie anti-TNF response in CD.

Genome-wide pathway analysis identifies VEGF pathway association with oral ulceration in systemic lupus erythematosus.

BACKGROUND: Systemic lupus erythematosus (SLE) is a genetically complex rheumatic disease characterized by heterogeneous clinical manifestations of unknown etiology. Recent studies have suggested the existence of a genetic basis for SLE heterogeneity. The objective of the present study was to identify new genetic variation associated with the clinically relevant phenotypes in SLE. METHODS: A two-stage pathway-based approach was used to identify the genetic variation associated with the main clinical phenotypes in SLE. In the discovery stage, 482 SLE patients were genotyped using Illumina Human Quad610 microarrays. Association between 798 reference genetic pathways from the Molecular Signatures Database and 11 SLE phenotypes was tested using the set-based method implemented in PLINK software. Pathways significantly associated after multiple test correction were subsequently tested for replication in an independent cohort of 425 SLE patients. Using an in silico approach, we analyzed the functional effects of common SLE therapies on the replicated genetic pathways. The association of known SLE risk variants with the development of the clinical phenotypes was also analyzed. RESULTS: In the discovery stage, we found a significant association between the vascular endothelial growth factor (VEGF) pathway and oral ulceration (P value for false discovery rate (P FDR) < 0.05), and between the negative regulation signaling pathway of retinoic acid inducible gene-I/melanoma differentiation associated gene 5 and the production of antinuclear antibodies (P FDR < 0.05). In the replication stage, we validated the association between the VEGF pathway and oral ulceration. Therapies commonly used to treat mucocutaneous phenotypes in SLE were found to strongly influence VEGF pathway gene expression (P = 4.60e-4 to 5.38e-14). Analysis of known SLE risk loci identified a strong association between PTPN22 and the risk of hematologic disorder and with the development of antinuclear antibodies. CONCLUSIONS: The present study has identified VEGF genetic pathway association with the risk of oral ulceration in SLE. New therapies targeting the VEGF pathway could be more effective in reducing the severity of this phenotype. These findings represent a first step towards the understanding of the genetic basis of phenotype heterogeneity in SLE.

A deletion at ADAMTS9-MAGI1 locus is associated with psoriatic arthritis risk.

OBJECTIVE: Copy number variants (CNVs) have been associated with the risk to develop multiple autoimmune diseases. Our objective was to identify CNVs associated with the risk to develop psoriatic arthritis (PsA) using a genome-wide analysis approach. METHODS: A total of 835 patients with PsA and 1498 healthy controls were genotyped for CNVs using the Illumina HumanHap610 BeadChip genotyping platform. Genomic CNVs were characterised using CNstream analysis software and analysed for association using the χ(2) test. The most significant genomic CNV associations with PsA risk were independently tested in a validation sample of 1133 patients with PsA and 1831 healthy controls. In order to test for the specificity of the variants with PsA aetiology, we also analysed the association to a cohort of 822 patients with purely cutaneous psoriasis (PsC). RESULTS: A total of 165 common CNVs were identified in the genome-wide analysis. We found a highly significant association of an intergenic deletion between ADAMTS9 and MAGI1 genes on chromosome 3p14.1 (p=0.00014). Using the independent patient and control cohort, we validated the association between ADAMTS9-MAGI1 deletion and PsA risk (p=0.032). Using next-generation sequencing, we characterised the 26 kb associated deletion. Finally, analysing the PsC cohort we found a lower frequency of the deletion compared with the PsA cohort (p=0.0088) and a similar frequency to that of healthy controls (p>0.3). CONCLUSIONS: The present genome-wide scan for CNVs associated with PsA risk has identified a new deletion associated with disease risk and which is also differential from PsC risk.

Variation at FCGR2A and functionally related genes is associated with the response to anti-TNF therapy in rheumatoid arthritis.

OBJECTIVE: Anti-TNF therapies have been highly efficacious in the management of rheumatoid arthritis (RA), but 25-30% of patients do not show a significant clinical response. There is increasing evidence that genetic variation at the Fc receptor FCGR2A is associated with the response to anti-TNF therapy. We aimed to validate this genetic association in a patient cohort from the Spanish population, and also to identify new genes functionally related to FCGR2A that are also associated with anti-TNF response. METHODS: A total of 348 RA patients treated with an anti-TNF therapy were included and genotyped for FCGR2A polymorphism rs1081274. Response to therapy was determined at 12 weeks, and was tested for association globally and independently for each anti-TNF drug (infliximab, etanercept and adalimumab). Using gene expression profiles from macrophages obtained from synovial fluid of RA patients, we searched for genes highly correlated with FCGR2A expression. Tag SNPs were selected from each candidate gene and tested for association with the response to therapy. RESULTS: We found a significant association between FCGR2A and the response to adalimumab (P=0.022). Analyzing the subset of anti-CCP positive RA patients (78%), we also found a significant association between FCGR2A and the response to infliximab (P=0.035). DHX32 and RGS12 were the most consistently correlated genes with FCGR2A expression in RA synovial fluid macrophages (P<0.001). We found a significant association between the genetic variation at DHX32 (rs12356233, corrected P=0.019) and a nominally significant association between RGS12 and the response to adalimumab (rs4690093, uncorrected P=0.040). In the anti-CCP positive group of patients, we also found a nominally significant association between RGS12 and the response to infliximab (rs2857859, uncorrected P=0.042). CONCLUSIONS: In the present study we have validated the FCGR2A association in an independent population, and we have identified new genes associated with the response to anti-TNF therapy in RA.

Identification of risk loci for Crohn's disease phenotypes using a genome-wide association study.

BACKGROUND & AIMS: Crohn's disease is a highly heterogeneous inflammatory bowel disease comprising multiple clinical phenotypes. Genome-wide association studies (GWASs) have associated a large number of loci with disease risk but have not associated any specific genetic variants with clinical phenotypes. We performed a GWAS of clinical phenotypes in Crohn's disease. METHODS: We genotyped 576,818 single-nucleotide polymorphisms in a well-characterized cohort of 1090 Crohn's disease patients of European ancestry. We assessed their association with 17 phenotypes of Crohn's disease (based on disease location, disease behavior, disease course, age at onset, and extraintestinal manifestations). A total of 57 markers with strong associations to Crohn's disease phenotypes (P < 2 × 10(-4)) were subsequently analyzed in an independent replication cohort of 1296 patients of European ancestry. RESULTS: We replicated the association of 4 loci with different Crohn's disease phenotypes. Variants in MAGI1, CLCA2, 2q24.1, and LY75 loci were associated with a complicated stricturing disease course (Pcombined = 2.01 × 10(-8)), disease location (Pcombined = 1.3 × 10(-6)), mild disease course (Pcombined = 5.94 × 10(-7)), and erythema nodosum (Pcombined = 2.27 × 10(-6)), respectively. CONCLUSIONS: In a GWAS, we associated 4 loci with clinical phenotypes of Crohn's disease. These findings indicate a genetic basis for the clinical heterogeneity observed for this inflammatory bowel disease.

PDE3A-SLCO1C1 locus is associated with response to anti-tumor necrosis factor therapy in psoriatic arthritis.

Aim: Variation at PDE3A-SLCO1C1 locus has been recently associated with the response to anti-TNF therapy in rheumatoid arthritis. We undertook the present study to determine whether PDE3A-SLCO1C1 is also associated with the response to anti-TNF therapy in psoriatic arthritis. Patients & methods: Genomic DNA was obtained from 81 psoriatic arthritis patients that had been treated with anti-TNF therapy. PDE3A-SLCO1C1 SNP rs3794271 was genotyped using Taqman realt-time PCR. The clinical response to anti-TNF therapy was measured as the change from baseline in the level of disease activity according to the DAS28 score. Results: A significant association between rs3794271 and anti-TNF response in psoriatic arthritis was found (beta = -0.71; p = 0.0036). Conclusion: PDE3A-SLCO1C1 locus is also associated with response to anti-TNF therapy in psoriatic arthritis. Original submitted 12 May 2014; Revision submitted 18 August 2014.

GWAS replication study confirms the association of PDE3A-SLCO1C1 with anti-TNF therapy response in rheumatoid arthritis.

AIM: The present study was undertaken to replicate the association of candidate genes for anti-TNF response in rheumatoid arthritis. Candidate genes were selected from a recent genome-wide association study on anti-TNF response performed in a population from Denmark. MATERIALS & METHODS: Genomic DNA was obtained from 315 Spanish rheumatoid arthritis patients having received an anti-TNF agent as their first biological therapy. SNPs from NR2FR2, MAP2K6, CBLN2 and PDE3A-SLCO1C1 candidate loci were genotyped. RESULTS: The PDE3A-SLCO1C1 locus rs3794271 SNP showed a highly significant association with anti-TNF treatment response (p = 1.74 × 10⁻5). Combining the statistical evidence from the Spanish and Danish rheumatoid arthritis cohorts, the associated rs3794271 SNP reached a genome-wide significance level of association (p = 3.3 × 10⁻¹°). CONCLUSION: The present findings establish the PDE3A-SLCO1C1 locus as a strong genetic marker of anti-TNF therapy response.

Baseline comorbidities in patients with rheumatoid arthritis who have been prescribed biological therapy: a case control study.

AIM: To determine whether rheumatoid arthritis (RA) patients who have been prescribed biological agents exhibit a different comorbidity burden than RA patients who take disease-modifying antirheumatic drugs (DMARDs) alone, and to understand the association between comorbidity and other variables, as well as the association between comorbidity and multimorbidity. METHODS: This observational case-control study included 114 RA patients treated with biological agents and a control group comprising 163 sex- and age-matched RA patients treated with DMARDs only. Current and previous data regarding the patients' disease activity, comorbidities, and treatments were collected. The data were analysed using bivariate and multivariate regression models. RESULTS: The patients who were prescribed biological agents exhibited poorer disease control, received more DMARDs and steroids, and underwent more total joint arthroplasties compared with the patients in the control group. However, the risk factors for cardiovascular disease and the comorbidity frequency were similar between cases and controls. The most prevalent comorbidities were hypertension, obesity, and respiratory, thyroid, and upper gastrointestinal disorders. The incidence of cardiovascular disease was low, and only 29% of the patients exhibited multimorbidities. A bivariate association of age, late diagnosis, joint replacements and a high score on the health assessment questionnaire score (HAQ) with comorbidity was observed. There were also correlations between the Charlson index and age, joint reconstructive surgery, disease activity (DAS28), and HAQ score. However, when binary logarithmic regression models were applied, only patient age remained significantly associated with comorbidity and multimorbidity [hazard ratio, 1.08; 95% confidence interval, 1.05-1.12; p<0.0005]. CONCLUSION: RA patients taking biological drugs have a comorbidity burden equivalent to those treated with DMARDs alone. Age is the main predictive factor of comorbidity in these patients.

A genome-wide association study on a southern European population identifies a new Crohn's disease susceptibility locus at RBX1-EP300.

OBJECTIVE: Genome-wide association studies (GWAS) have identified multiple risk loci for Crohn's disease (CD). However, the cumulative risk exerted by these loci is low, and the likelihood that additional, as-yet undiscovered loci contribute to the risk of CD is very high. We performed a GWAS on a southern European population to identify new CD risk loci. DESIGN: We genotyped 620 901 genome markers on 1341 CD patients and 1518 controls from Spain. The top association signals representing new candidate risk loci were subsequently analysed in an independent replication cohort of 1365 CD patients and 1396 controls. RESULTS: We identified a genome-wide significant association on chromosome 22q13.2 in the intergenic region between the RBX1 and EP300 genes (single nucleotide polymorphism rs4820425, OR 1.27, 95% CI 1.17 to 1.38, p=3.42E-8). We also found suggestive evidence for the association of the IFNGR2 (21q22.11), FOXP2 (7q31), MACROD2 (20p12.1) and AIF1 (6p21.3) loci with CD risk. CONCLUSIONS: In this GWAS performed on a southern European cohort, we have identified a new risk locus for CD between RBX1 and EP300. This study demonstrates that using populations of different ancestry is a useful strategy to identify new risk loci for CD.

Effectiveness, predictive response factors, and safety of anti-tumor necrosis factor (TNF) therapies in anti-TNF-naive rheumatoid arthritis.

OBJECTIVE: To evaluate the effectiveness and safety of anti-tumor necrosis factor (anti-TNF) therapies in rheumatoid arthritis (RA), and to identify the factors involved in this response. METHODS: Dynamic prospective cohort study of patients with RA treated with anti-TNF under clinical practice conditions. Effectiveness was evaluated using Disease Activity Score (DAS) 28, European League Against Rheumatism (EULAR) response, Health Assessment Questionnaire (HAQ), and time to treatment failure. Prior adherence was evaluated retrospectively and safety was evaluated by adverse events (AE). The analysis was restricted to anti-TNF-naive patients. RESULTS: The study included 161 patients treated for RA during 6 years (60 infliximab, 79 etanercept, and 22 adalimumab). At 6 months, 15% reached a good EULAR response and 38% a moderate response. A mean decrease of -1.5 (p < 0.0001) was observed in the DAS28 and of -0.34 in the HAQ (p < 0.0001); however, women showed poorer progress in terms of DAS and HAQ. In the first year, 64.3% did not experience treatment failure and this figure was 50.5% after 2 years. In one-third, glucocorticoids were withdrawn and in the remainder the dose was reduced by 50%. Adherence to treatment, selection of etanercept, and intensification of infliximab were associated with a lower probability of premature failure in the multivariate model. AE were similar to other those in studies and no outstanding differences in safety were found between the 3 anti-TNF therapies. CONCLUSIONS: Anti-TNF treatments are effective and safe, reducing the activity of the disease, disability, and the need for corticosteroids. Patients who displayed good adherence prior to the anti-TNF treatment and were treated with etanercept or with increasing doses of infliximab had the best chance of displaying a response.