The ovarian reserve as measured by the anti-Müllerian hormone is not diminished in patients with rheumatoid arthritis compared to the healthy population.

OBJECTIVES: Rheumatoid arthritis (RA) is the most prevalent chronic inflammatory arthritis, affecting 0.5-1% worldwide population and predominates in females. Altered fertility has been reported due to a decrease in ovarian reserve secondary to sustained inflammation. The anti-Müllerian Hormone (AMH) is currently the most reliable biomarker of ovarian reserve. However, few and contradictory studies have been reported to analyse the relationship between fertility in RA female patients and AMH. The aim of present study is to determine the AMH serum concentrations in a long-standing RA patient group and control group. We also sought to determine the correlation between AMH serum levels and disease activity measured by different parameters and the effect of biological DMARDs. METHODS: Serum AMH levels were measured in 60 women with long-standing RA aged 20-50 y.o. and compared to 59 healthy women. AMH was assessed by an electrochemiluminescence immunoassay method (ECLIA, Roche Diagnostics) and a large data set of clinical and molecular data was annotated. Demographic parameters, RA disease activity measured by DAS28 score and inflammatory biomarkers such as ESR, CRP, lymphocyte CD4+, CD8+, NK cells, IL-10 and IL-6 were determined. A comprehensive gynaecological self-administered questionnaire was given. Serum AMH levels were age-correlated. Differences between groups were calculated using Student's t-test or Mann-Whitney U-test for continuous variables and Fisher's exact test for categorical variables. Multivariate analysis was conducted by the partial correlation coefficient. Linear regression analysis was performed to study the effect of different variables on proportional AMH change. p-values <0.005 were considered significant. RESULTS: The median age was similar in AR and control groups (37.4±6.23 vs. 37.3±6.27 p=0.937). Mean disease duration was 8.37±5.36 years. The number of previous treatments was <3 in 71.7% of patients and ≥3 in 28.3%. Disease activity measured by DAS28 was 2.89± 1.54. The age-adjusted mean serum concentration of AMH was 1.27 ng/ml [IQR 0.42; 2.24] in RA patients and 1.31 ng/ml [IQR 0.46; 3.09] in controls (p=0.608). Neither disease activity (p=0.862), nor current or previous bDMARDs treatments (p=0.871) were associated with AMH levels. However, a negative linear correlation was observed between AMH and IL-10 levels (p=0.033). CONCLUSIONS: Our study shows that ovarian reserve determined by AMH serum levels is not reduced in rheumatoid arthritis patients compared with healthy controls. In our series, AMH levels were not affected by disease activity, however, a significant correlation was observed between AMH and IL-10 levels. These results support the role of cytokines profile in the female reproductive system and will focus further investigations in this critical area, mainly once biological DMARDs have been recommended in RA pregnant patients.

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.

Epigenome-wide association study of rheumatoid arthritis identifies differentially methylated loci in B cells.

Epigenetic regulation of immune cell types could be critical for the development and maintenance of autoimmune diseases like rheumatoid arthritis (RA). B cells are highly relevant in RA, since patients express autoantibodies and depleting this cell type is a successful therapeutic approach. Epigenetic variation, such as DNA methylation, may mediate the pathogenic activity of B cells. In this study, we performed an epigenome-wide association study (EWAS) for RA with three different replication cohorts, to identify disease-specific alterations in DNA methylation in B cells. CpG methylation in isolated B lymphocytes was assayed on the Illumina HumanMethylation450 BeadChip in a discovery cohort of RA patients (N = 50) and controls (N = 75). Differential methylation was observed in 64 CpG sites (q < 0.05). Six biological pathways were also differentially methylated in RA B cells. Analysis in an independent cohort of patients (N = 15) and controls (N = 15) validated the association of 10 CpG sites located on 8 genes CD1C, TNFSF10, PARVG, NID1, DHRS12, ITPK1, ACSF3 and TNFRSF13C, and 2 intergenic regions. Differential methylation at the CBL signaling pathway was replicated. Using an additional case-control cohort (N = 24), the association between RA risk and CpGs cg18972751 at CD1C (P = 2.26 × 10-9) and cg03055671 at TNFSF10 (P = 1.67 × 10-8) genes was further validated. Differential methylation at genes CD1C, TNFSF10, PARVG, NID1, DHRS12, ITPK1, ACSF3, TNFRSF13C and intergenic region chr10p12.31 was replicated in a cohort of systemic lupus erythematosus (SLE) patients (N = 47) and controls (N = 56). Our results highlight genes that may drive the pathogenic activity of B cells in RA and suggest shared methylation patterns with SLE.

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.

A genome-wide association study identifies a novel locus at 6q22.1 associated with ulcerative colitis.

The genetic analysis of ulcerative colitis (UC) has provided new insights into the etiology of this prevalent inflammatory bowel disease. However, most of the heritability of UC (>70%) has still not been characterized. To identify new risk loci for UC we have performed the first genome-wide association study (GWAS) in a Southern European population and undertaken a meta-analysis study combining the newly genotyped 825 UC patients and 1525 healthy controls from Spain with the six previously published GWAS comprising 6687 cases and 19 718 controls from Northern-European ancestry. We identified a novel locus with genome-wide significance at 6q22.1 [rs2858829, P = 8.97 × 10(-9), odds ratio (OR) (95% confidence interval, CI] = 1.12 (1.08-1.16)] that was validated with genotype data from a replication cohort of the same Southern European ancestry consisting in 1073 cases and 1279 controls [combined P = 7.59 × 10(-10), OR (95% CI) = 1.12 (1.08-1.16)]. Furthermore, we confirmed the association of 33 reported associations with UC and we nominally validated the GWAS results of nine new risk loci (P < 0.05, same direction of effect). SNP rs2858829 lies in an intergenic region and is a strong cis-eQTL for FAM26F gene, a gene that is shown to be selectively upregulated in UC colonic mucosa with active inflammation. Our results provide new insight into the genetic risk background of UC, confirming that there is a genetic risk component that differentiates from Crohn's Disease, the other major form of inflammatory bowel disease.

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.

A genome-wide association study identifies SLC8A3 as a susceptibility locus for ACPA-positive rheumatoid arthritis.

OBJECTIVE: RA patients with serum ACPA have a strong and specific genetic background. The objective of the study was to identify new susceptibility genes for ACPA-positive RA using a genome-wide association approach. METHODS: A total of 924 ACPA-positive RA patients with joint damage in hands and/or feet, and 1524 healthy controls were genotyped in 582 591 single-nucleotide polymorphisms (SNPs) in the discovery phase. In the validation phase, the most significant SNPs in the genome-wide association study representing new candidate loci for RA were tested in an independent cohort of 863 ACPA-positive patients with joint damage and 1152 healthy controls. All individuals from the discovery and validation cohorts were Caucasian and of Southern European ancestry. RESULTS: In the discovery phase, 60 loci not previously associated with RA risk showed evidence for association at P < 5×10(-4) and were tested for replication in the validation cohort. A total of 12 loci were replicated at the nominal level (P < 0.05, same direction of effect as in the discovery phase). When combining the discovery and validation cohorts, an intronic SNP in the Solute Carrier family 8 gene (SLC8A3) was found to be associated with ACPA-positive RA at a genome-wide level of significance RA [odds ratio (95% CI): 1.42 (1.25, 1.6), Pcombined = 3.19×10(-8)]. CONCLUSIONS: SLC8A3 was identified as a new risk locus for ACPA-positive RA. This study demonstrates the advantage of analysing relevant subsets of RA patients to identify new genetic risk variants.

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.

Risk variants for psoriasis vulgaris in a large case-control collection and association with clinical subphenotypes.

Recent genome-wide association studies (GWASs) have identified >20 new loci associated with the susceptibility to psoriasis vulgaris (PsV) risk. We investigated the association of PsV and its main clinical subphenotypes with 32 loci having previous genome-wide evidence of association with PsV (P < 5e-8) or strong GWAS evidence (P < 5e-5 in discovery and P < 0.05 in replication sample) in a large cohort of PsV patients (n = 2005) and controls (n = 1497). We provide the first independent replication for COG6 (P = 0.00079) and SERPINB8 (P = 0.048) loci with PsV. In those patients having developed psoriatic arthritis (n = 955), we found, for the first time, a strong association with IFIH1 (P = 0.013). Analyses of clinically relevant PsV subtypes yielded a significant association of severity of cutaneous disease with variation at LCE3D locus (P = 0.0005) in PsV and nail involvement with IL1RN in purely cutaneous psoriasis (PsC, P = 0.007). In an exploratory analysis of epistasis, we replicated the previously described HLA-C-ERAP1 interaction with PsC. Our findings show that common genetic variants associated with a complex phenotype like PsV influence different subphenotypes of high clinical relevance.

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.

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.

Genome-wide transcriptional analysis of T cell activation reveals differential gene expression associated with psoriasis.

BACKGROUND: Psoriasis is a chronic autoimmune disease in which T cells have a predominant role in initiating and perpetuating the chronic inflammation in skin. However, the mechanisms that regulate T cell activation in psoriasis are still incompletely understood. The objective of the present study was to characterize the main genetic pathways associated with T cell activation in psoriasis. RESULTS: Gene expression profiles from in vitro activated T cells were obtained from 17 psoriasis patients and 7 healthy controls using Illumina HT-12 v4 microarrays. From a total of 47,321 analyzed transcripts, 42 genes were found to be differentially expressed between psoriasis and controls (FDR p-value < 0.1, absolute fold-change > 1.2). Using an independent cohort of 8 patients and 8 healthy controls we validated the overexpression of SPATS2L (p-value =0.0009) and KLF6 (p-value =0.0012) genes in activated T cells from psoriasis patients. Using weighted correlation analysis we identified SPATS2L and KLF6 coexpression networks, which were also significantly associated with psoriasis (p-value < 0.05). Gene Ontology analysis allowed the identification of several biological processes associated with each coexpression network. Finally, using Gene Set Enrichment Analysis over the global T cell transcriptome we also found additional genetic pathways strongly associated with psoriasis (p-value < 0.0001). CONCLUSIONS: This study has identified two new genes, SPATS2L and KLF6, strongly associated with T cell activation in psoriasis. Functional analyses of the gene expression profiles also revealed new biological processes and genetic pathways associated with psoriasis. The results of this study provide an important insight into the biology of this common chronic inflammatory disease.

Impact of the use of musculoskeletal ultrasound by rheumatologists in patients with shoulder and hand complaints compared with traditional clinical care.

OBJECTIVES: To compare the routine use of musculoskeletal ultrasonography (MSUS) with traditional clinical care in daily practice at shoulder and hand level. METHODS: An observational study was performed in four rheumatology departments. Within each department, 2 rheumatologists were selected; one rheumatologist used MSUS, and the other followed traditional rheumatology care. Consecutive patients with nontraumatic pain, hand numbness or disability, or pain and/or limitations in the shoulder were selected. We collected information regarding the clinical and MSUS diagnoses, changes in diagnosis and treatment following MSUS, local injections, the rheumatologist's satisfaction and the use of health care resources. A descriptive analysis was performed. RESULTS: A total of 168 patients were analysed, with 104 and 64 patients in the MSUS and traditional care groups, respectively. MSUS led to a diagnosis and therapeutic change in 53 (52%) and 55 patients (54%), respectively. The rate of local injection was 47% in the MSUS group (73% unexpected, 61% performed using US) compared with 21% in the traditional group (p=0.001). According to the rheumatologists, MSUS was useful in 72 cases (71%) and extremely useful in 20 cases (20%), and the rheumatologists reported a higher satisfaction with their patient evaluations (p<0.001). The MSUS group required fewer additional tests (38% vs. 81%, respectively, p<0.001), fewer medical visits (46% vs. 84%, p<0.001), and lower direct costs (11 vs. 30 euros, p<0.001) than the traditional care group. CONCLUSIONS: Compared with traditional care, the routine use of MSUS in rheumatology practice at hand and shoulder level can lead to important improvements in care, thereby reducing the number of additional tests and medical visits.

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.

Targeting of the CD80/86 proinflammatory axis as a therapeutic strategy to prevent severe COVID-19.

An excessive immune response known as cytokine storm is the hallmark of severe COVID-19. The cause of this cytokine rampage is yet not known. Based on recent epidemiological evidence, we hypothesized that CD80/86 signaling is essential for this hyperinflammation, and that blocking this proinflammatory axis could be an effective therapeutic approach to protect against severe COVID-19. Here we provide exploratory evidence that abatacept, a drug that blocks CD80/86 co-stimulation, produces changes at the systemic level that are highly antagonistic of the proinflammatory processes elicited by COVID-19. Using RNA-seq from blood samples from a longitudinal cohort of n = 38 rheumatic patients treated with abatacept, we determined the immunological processes that are significantly regulated by this treatment. We then analyzed available blood RNA-seq from two COVID19 patient cohorts, a very early cohort from the epicenter of the pandemic in China (n = 3 COVID-19 cases and n = 3 controls), and a recent and larger cohort from the USA (n = 49 severe and n = 51 mild COVD-19 patients). We found a highly significant antagonism between SARS-CoV-2 infection and COVID-19 severity with the systemic response to abatacept. Analysis of previous single-cell RNA-seq data from bronchoalveolar lavage fluid from mild and severe COVID-19 patients and controls, reinforce the implication of the CD80/86 proinflammatory axis. Our functional results further support abatacept as a candidate therapeutic approach to prevent severe COVID-19.

Non-invasive vagus nerve stimulation for rheumatoid arthritis: a proof-of-concept study.

BACKGROUND: Vagus nerve stimulation delivered with an implanted device has been shown to improve rheumatoid arthritis severity. We aimed to investigate the safety and efficacy of non-invasive stimulation of the auricular branch of the vagus nerve for the treatment of patients with moderately to severely active rheumatoid arthritis. METHODS: This prospective, multicentre, open-label, single-arm proof-of-concept study enrolled patients aged 18-80 years with active rheumatoid arthritis who had an inadequate response to conventional synthetic disease-modifying antirheumatic drugs (DMARDs) and up to one biological DMARD. Biological DMARDs were stopped at least 4 weeks before enrolment and concomitant use was not allowed during the study. All eligible participants were assigned to use a non-invasive, wearable vagus nerve stimulation device for up to 30 min per day, which delivered pulses of 20 kHz. Follow-up visits occurred at week 1, week 2, week 4, week 8, and week 12 after the baseline visit. The primary endpoint was the mean change in Disease Activity Score of 28 joints with C-reactive protein (DAS28-CRP) at week 12 compared with baseline. Secondary endpoints included the mean change in the Health Assessment Questionnaire-Disability Index (HAQ-DI), the proportion of patients with a minimal clinically important difference of 0·22 on HAQ-DI, the proportion achieving American College of Rheumatology (ACR) 20, ACR50, and ACR70 response, and safety analysis. This study is registered with ClinicalTrials.gov (NCT04116866). FINDINGS: Of 35 patients screened for eligibility, 30 (86%) were enrolled at six centres in Spain between Dec 27, 2018, and Oct 24, 2019, of whom 27 (90%) completed the week 12 visit. The mean change in DAS28-CRP at 12 weeks was -1·4 (95%CI -1·9 to -0·9; p<0·0001) from a mean baseline of 5·3 (SD 1·0). 11 (37%) of 30 patients reached DAS28-CRP of 3·2 or less, and seven (23%) patients reached DAS28-CRP of less than 2·6 at week 12. The mean HAQ-DI change was -0·5 (95%CI -0·7 to -0·2; p<0·0001) from a mean baseline of 1·6 (SD 0·7), and 17 (57%) patients reached a minimal clinically important difference of 0·22 or more. ACR20 responses were reached by 16 (53%) patients, ACR50 responses by 10 (33%) patients, and ACR70 by five (17%) patients. Four adverse events were reported, none of which were serious and all of which resolved without intervention. INTERPRETATION: Use of the device was well tolerated, and patients had clinically meaningful reductions in DAS28-CRP. This was an uncontrolled, open-label study, and the results must be interpreted in this context. Further evaluation in larger, controlled studies is needed to confirm whether this non-invasive approach might offer an alternative treatment for rheumatoid arthritis. FUNDING: Nesos.

Incidence of COVID-19 in a cohort of adult and paediatric patients with rheumatic diseases treated with targeted biologic and synthetic disease-modifying anti-rheumatic drugs.

Objectives: To investigate the incidence of COVID-19 in a cohort of adult and paediatric patients with rheumatic diseases receiving targeted biologic and synthetic disease modifying anti-rheumatic drugs (tDMARDs) and to explore the possible effect of these treatments in the clinical expression of COVID-19. Methods: A cross-sectional study comprising of a telephone survey and electronic health records review was performed including all adult and paediatric patients with rheumatic diseases treated with tDMARDs in a large rheumatology tertiary centre in Barcelona, Spain. Demographics, disease activity, COVID-19 related symptoms and contact history data were obtained from the start of the 2020 pandemic. Cumulative incidence of confirmed cases (SARS-CoV-2 positive PCR test) was compared to the population estimates for the same city districts from a governmental COVID-19 health database. Suspected cases were defined following WHO criteria and compared to those without compatible symptoms. Results: 959 patients with rheumatic diseases treated with tDMARDs were included. We identified 11 confirmed SARS-CoV-2 positive cases in the adult cohort and no confirmed positive cases in the paediatric cohort. COVID-19 incidence rates of the rheumatic patient cohort were very similar to that of the general population [(0.48% (95% CI 0.09 to 0.87%)] and [0.58% (95% CI 0.56 to 0.60%)], respectively. We found significant differences in tDMARDs proportions between the suspected and non-suspected cases (p=0.002). Conclusion: Adult and paediatric patients with rheumatic diseases on tDMARDs do not seem to present a higher risk of COVID-19 or a more severe disease outcome when compared to general population.

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.

Genome-wide association study meta-analysis identifies five new loci for systemic lupus erythematosus.

BACKGROUND: Systemic lupus erythematosus (SLE) is a common systemic autoimmune disease with a complex genetic inheritance. Genome-wide association studies (GWAS) have significantly increased the number of significant loci associated with SLE risk. To date, however, established loci account for less than 30% of the disease heritability and additional risk variants have yet to be identified. Here we performed a GWAS followed by a meta-analysis to identify new genome-wide significant loci for SLE. METHODS: We genotyped a cohort of 907 patients with SLE (cases) and 1524 healthy controls from Spain and performed imputation using the 1000 Genomes reference data. We tested for association using logistic regression with correction for the principal components of variation. Meta-analysis of the association results was subsequently performed on 7,110,321 variants using genetic data from a large cohort of 4036 patients with SLE and 6959 controls of Northern European ancestry. Genetic association was also tested at the pathway level after removing the effect of known risk loci using PASCAL software. RESULTS: We identified five new loci associated with SLE at the genome-wide level of significance (p < 5 × 10- 8): GRB2, SMYD3, ST8SIA4, LAT2 and ARHGAP27. Pathway analysis revealed several biological processes significantly associated with SLE risk: B cell receptor signaling (p = 5.28 × 10- 6), CTLA4 co-stimulation during T cell activation (p = 3.06 × 10- 5), interleukin-4 signaling (p = 3.97 × 10- 5) and cell surface interactions at the vascular wall (p = 4.63 × 10- 5). CONCLUSIONS: Our results identify five novel loci for SLE susceptibility, and biologic pathways associated via multiple low-effect-size loci.