A survey of ficolin-3 activity in Systemic Lupus Erythematosus reveals a link to hematological disease manifestations and autoantibody profile.

The complement system plays a central role in the pathogenesis of Systemic Lupus Erythematosus (SLE), but most studies have focused on the classical pathway. Ficolin-3 is the main initiator of the lectin pathway of complement in humans, but its role in systemic autoimmune disease has not been conclusively determined. Here, we combined biochemical and genetic approaches to assess the contribution of ficolin-3 to SLE risk and disease manifestations. Ficolin-3 activity was measured by a functional assay in serum or plasma samples from Swedish SLE patients (n = 786) and controls matched for age and sex (n = 566). Genetic variants in an extended 300 kb genomic region spanning the FCN3 locus were analyzed for their association with ficolin-3 activity and SLE manifestations in a Swedish multicenter cohort (n = 985). Patients with ficolin-3 activity in the highest tertile showed a strong enrichment in an SLE cluster defined by anti-Sm/DNA/nucleosome antibodies (OR 3.0, p < 0.001) and had increased rates of hematological disease (OR 1.4, p = 0.078) and lymphopenia (OR = 1.6, p = 0.039). Genetic variants associated with low ficolin-3 activity mapped to an extended haplotype in high linkage disequilibrium upstream of the FCN3 gene. Patients carrying the lead genetic variant associated with low ficolin-3 activity had a lower frequency of hematological disease (OR 0.67, p = 0.018) and lymphopenia (OR 0.63, p = 0.031) and fewer autoantibodies (p = 0.0019). Loss-of-function variants in the FCN3 gene were not associated with SLE, but four (0.5 %) SLE patients developed acquired ficolin-3 deficiency where ficolin-3 activity in serum was depleted following diagnosis of SLE. Taken together, our results provide genetic and biochemical evidence that implicate the lectin pathway in hematological SLE manifestations. We also identify lectin pathway activation through ficolin-3 as a factor that contributes to the autoantibody response in SLE.

Distinct HLA associations with autoantibody-defined subgroups in idiopathic inflammatory myopathies.

BACKGROUND: In patients with idiopathic inflammatory myopathies (IIM), autoantibodies are associated with specific clinical phenotypes suggesting a pathogenic role of adaptive immunity. We explored if autoantibody profiles are associated with specific HLA genetic variants and clinical manifestations in IIM. METHODS: We included 1348 IIM patients and determined the occurrence of 14 myositis-specific or -associated autoantibodies. We used unsupervised cluster analysis to identify autoantibody-defined subgroups and logistic regression to estimate associations with clinical manifestations, HLA-DRB1, HLA-DQA1, HLA-DQB1 alleles, and amino acids imputed from genetic information of HLA class II and I molecules. FINDINGS: We identified eight subgroups with the following dominant autoantibodies: anti-Ro52, -U1RNP, -PM/Scl, -Mi2, -Jo1, -Jo1/Ro52, -TIF1γ or negative for all analysed autoantibodies. Associations with HLA-DRB1∗11, HLA-DRB1∗15, HLA-DQA1∗03, and HLA-DQB1∗03 were present in the anti-U1RNP-dominated subgroup. HLA-DRB1∗03, HLA-DQA1∗05, and HLA-DQB1∗02 alleles were overrepresented in the anti-PM/Scl and anti-Jo1/Ro52-dominated subgroups. HLA-DRB1∗16, HLA-DRB1∗07 alleles were most frequent in anti-Mi2 and HLA-DRB1∗01 and HLA-DRB1∗07 alleles in the anti-TIF1γ subgroup. The HLA-DRB1∗13, HLA-DQA1∗01 and HLA-DQB1∗06 alleles were overrepresented in the negative subgroup. Significant signals from variations in class I molecules were detected in the subgroups dominated by anti-Mi2, anti-Jo1/Ro52, anti-TIF1γ, and the negative subgroup. INTERPRETATION: Distinct HLA class II and I associations were observed for almost all autoantibody-defined subgroups. The associations support autoantibody profiles use for classifying IIM which would likely reflect underlying pathogenic mechanisms better than classifications based on clinical symptoms and/or histopathological features. FUNDING: See a detailed list of funding bodies in the Acknowledgements section at the end of the manuscript.

Single-cell profiling of muscle-infiltrating T cells in idiopathic inflammatory myopathies.

Idiopathic inflammatory myopathies (IIM) are rare autoimmune systemic diseases characterized by muscle weakness and the presence of muscle-infiltrating T cells. IIM represent a clinical challenge due to heterogeneity of symptoms and variability of response to immunosuppressive treatment. Here, we performed in-depth single-cell sequencing on muscle-infiltrating T cells and peripheral blood memory T cells in six patients with recently diagnosed IIM. We identified tissue resident memory T-cell (TRM ) signatures including the expression of HOBIT, XCL1 and CXCR6 in the muscle biopsies of all patients with IIM. Clonally expanded T-cell clones were mainly found among cytotoxic and TRM implying their role in the disease pathogenesis. Finally, identical expanded T-cell clones persisting at follow-up in the muscle tissue of two patients suggest their involvement in disease chronicity. Our study reveals a muscle tissue resident memory T-cell signature in patients with IIM and a transcriptomic map to identify novel therapeutic targets in IIM.

Genome-wide investigation of persistence to treatment with methotrexate in early rheumatoid arthritis.

OBJECTIVES: To investigate the influence of genetic factors on persistence to treatment of early rheumatoid arthritis (RA) with methotrexate (MTX) monotherapy. METHODS: We conducted a genome-wide association study (GWAS) in a sample of 3902 Swedish early RA patients initiating MTX in DMARD-monotherapy as their first ever DMARD. The outcome, short- and long-term persistence to this treatment, was defined as remaining on MTX at one and at three years, respectively, with no additional DMARDs added. As genetic predictors, we investigated individual SNPs, and a polygenic risk score (PRS) based on SNPs associated with RA risk. The SNP-based heritability of persistence was estimated overall and by RA serostatus. RESULTS: No individual SNP reached genome-wide significance (p < 5e-8), neither for persistence at one nor at three years. The RA PRS was not significantly associated with persistence at one (RR = 0.98 (0.96-1.01)) nor three years (RR = 0.96 (0.93-1.00)). The heritability for persistence was estimated to be 0.45 (0.15-0.75) at one year and 0.14 (0-0.40) at three years. Results in seropositive RA were comparable to those in the analysis of RA overall, while heritability estimates and PRS RRs were attenuated towards the null in seronegative RA. CONCLUSIONS: Despite being the largest GWAS on an MTX treatment outcome to date, no genome-wide significant associations were detected. The modest heritability observed, coupled with the broad spread of suggestively associated loci, indicate that genetic influence is of polygenic nature. Nevertheless, persistence to MTX monotherapy was lower in patients with a greater genetic disposition, per the PRS, towards RA.

Repurposing antidiabetic drugs for rheumatoid arthritis: results from a two-sample Mendelian randomization study.

Despite increasing therapeutic options to treat rheumatoid arthritis (RA), many patients fail to reach treatment targets. The use of antidiabetic drugs like thiazolidinediones has been associated with lower RA risk. We aimed to explore the repurposing potential of antidiabetic drugs in RA prevention by assessing associations between genetic variation in antidiabetic drug target genes and RA using Mendelian randomization (MR). A two-sample MR design was used to estimate the association between the antidiabetic drug and RA risk using summary statistics from genome-wide association studies (GWAS). We selected independent genetic variants from the gene(s) that encode the target protein(s) of the investigated antidiabetic drug as instruments. We extracted the associations of instruments with blood glucose concentration and RA from the UK Biobank and a GWAS meta-analysis of clinically diagnosed RA, respectively. The effect of genetic variation in the drug target(s) on RA risk was estimated by the Wald ratio test or inverse-variance weighted method. Insulin and its analogues, thiazolidinediones, and sulfonylureas had valid genetic instruments (n = 1, 1, and 2, respectively). Genetic variation in thiazolidinedione target (gene: PPARG) was inversely associated with RA risk (odds ratio [OR] 0.38 per 0.1mmol/L glucose lowering, 95% confidence interval [CI] 0.20-0.73). Corresponding ORs (95%CIs) were 0.83 (0.44-1.55) for genetic variation in the targets of insulin and its analogues (gene: INSR), and 1.12 (0.83, 1.49) 1.25 (0.78-2.00) for genetic variation in the sulfonylurea targets (gene: ABCC8 and KCNJ11). In conclusion, genetic variation in the thiazolidinedione target is associated with a lower RA risk. The underlying mechanisms warrant further exploration.

The human bone marrow plasma cell compartment in rheumatoid arthritis - Clonal relationships and anti-citrulline autoantibody producing cells.

A majority of circulating IgG is produced by plasma cells residing in the bone marrow (BM). Long-lived BM plasma cells constitute our humoral immune memory and are essential for infection-specific immunity. They may also provide a reservoir of potentially pathogenic autoantibodies, including rheumatoid arthritis (RA)-associated anti-citrullinated protein autoantibodies (ACPA). Here we investigated paired human BM plasma cell and peripheral blood (PB) B-cell repertoires in seropositive RA, four ACPA+ RA patients and one ACPA- using two different single-cell approaches, flow cytometry sorting, and transcriptomics, followed by recombinant antibody generation. Immunoglobulin (Ig) analysis of >900 paired heavy-light chains from BM plasma cells identified by either surface CD138 expression or transcriptome profiles (including gene expression of MZB1, JCHAIN and XBP1) demonstrated differences in IgG/A repertoires and N-linked glycosylation between patients. For three patients, we identified clonotypes shared between BM plasma cells and PB memory B cells. Notably, four individuals displayed plasma cells with identical heavy chains but different light chains, which may indicate receptor revision or clonal convergence. ACPA-producing BM plasma cells were identified in two ACPA+ patients. Three of 44 recombinantly expressed monoclonal antibodies from ACPA+ RA BM plasma cells were CCP2+, specifically binding to citrullinated peptides. Out of these, two clones reacted with citrullinated histone-4 and activated neutrophils. In conclusion, single-cell investigation of B-cell repertoires in RA bone marrow provided new understanding of human plasma cells clonal relationships and demonstrated pathogenically relevant disease-associated autoantibody expression in long-lived plasma cells.

Anti-phosphatidylserine/prothrombin antibodies and thrombosis associate positively with HLA-DRB1*13 and negatively with HLA-DRB1*03 in SLE.

OBJECTIVES: Emerging evidence demonstrates that aPS-PT associate with thrombotic events. Genetic predisposition, including HLA-DRB1 alleles, is known to contribute to the occurrence of conventional aPL [anti-ß2glycoprotein-I (anti-ß2GPI) and aCL]. We investigated associations between aPS-PT and HLA-DRB1* alleles and thrombosis in SLE. Conventional aPL were included for comparison. METHODS: We included 341 consecutive SLE patients, with information on general cardiovascular risk factors, including blood lipids, LA and thrombotic events. aPS/PT, anti-ß2GPI and aCL of IgA/G/M isotypes and LA were quantified. RESULTS: aPS/PT antibodies associated positively with HLA-DRB1*13 [odds ratio (OR) 2.7, P = 0.002], whereas anti-ß2GPI and aCL antibodies associated primarily with HLA-DRB1*04 (OR 2.5, P = 0.0005). These associations remained after adjustment for age, gender and other HLA-DRB1* alleles. HLA-DRB1*13, but not DRB1*04, remained as an independent risk factor for thrombosis and APS after adjustment for aPL and cardiovascular risk factors. The association between DRB1*13 and thrombosis was mediated by aPS-PT positivity. HLA-DRB1*03, on the other hand, associated negatively with thrombotic events as well as all aPL using both uni- and multivariate analyses. HLA-DRB1*03 had a thrombo-protective effect in aPL-positive patients. Additionally, HLA-DRB1*03 was associated with a favourable lipid profile regarding high-density lipoprotein and triglycerides. CONCLUSIONS: HLA-DRB1*13 confers risk for both aPS-PT and thrombotic events in lupus. The association between HLA-DRB1*13 and thrombosis is largely, but not totally, mediated through aPS-PT. HLA-DRB1*03 was negatively associated with aPL and positively with favourable lipid levels. Thus, HLA-DRB1*03 seems to identify a subgroup of SLE patients with reduced vascular risk.

Complement C4 Copy Number Variation is Linked to SSA/Ro and SSB/La Autoantibodies in Systemic Inflammatory Autoimmune Diseases.

OBJECTIVE: Copy number variation of the C4 complement components, C4A and C4B, has been associated with systemic inflammatory autoimmune diseases. This study was undertaken to investigate whether C4 copy number variation is connected to the autoimmune repertoire in systemic lupus erythematosus (SLE), primary Sjögren's syndrome (SS), or myositis. METHODS: Using targeted DNA sequencing, we determined the copy number and genetic variants of C4 in 2,290 well-characterized Scandinavian patients with SLE, primary SS, or myositis and 1,251 healthy controls. RESULTS: A prominent relationship was observed between C4A copy number and the presence of SSA/SSB autoantibodies, which was shared between the 3 diseases. The strongest association was detected in patients with autoantibodies against both SSA and SSB and 0 C4A copies when compared to healthy controls (odds ratio [OR] 18.0 [95% confidence interval (95% CI) 10.2-33.3]), whereas a weaker association was seen in patients without SSA/SSB autoantibodies (OR 3.1 [95% CI 1.7-5.5]). The copy number of C4 correlated positively with C4 plasma levels. Further, a common loss-of-function variant in C4A leading to reduced plasma C4 was more prevalent in SLE patients with a low copy number of C4A. Functionally, we showed that absence of C4A reduced the individuals' capacity to deposit C4b on immune complexes. CONCLUSION: We show that a low C4A copy number is more strongly associated with the autoantibody repertoire than with the clinically defined disease entities. These findings may have implications for understanding the etiopathogenetic mechanisms of systemic inflammatory autoimmune diseases and for patient stratification when taking the genetic profile into account.

Four Systemic Lupus Erythematosus Subgroups, Defined by Autoantibodies Status, Differ Regarding HLA-DRB1 Genotype Associations and Immunological and Clinical Manifestations.

OBJECTIVE: The heterogeneity of systemic lupus erythematosus (SLE) constitutes clinical and therapeutical challenges. We therefore studied whether unrecognized disease subgroups can be identified by using autoantibody profiling together with HLA-DRB1 alleles and immunological and clinical data. METHODS: An unsupervised cluster analysis was performed based on detection of 13 SLE-associated autoantibodies (double-stranded DNA, nucleosomes, ribosomal P, ribonucleoprotein [RNP] 68, RNPA, Smith [Sm], Sm/RNP, Sjögren's syndrome antigen A [SSA]/Ro52, SSA/Ro60, Sjögren's syndrome antigen B [SSB]/La, cardiolipin [CL]-Immunoglobulin G [IgG], CL-Immunoglobulin M [IgM], and ß2 glycoprotein I [ß2 GPI]-IgG) in 911 patients with SLE from two cohorts. We evaluated whether each SLE subgroup is associated with HLA-DRB1 alleles, clinical manifestations (n = 743), and cytokine levels in circulation (n = 446). RESULTS: Our analysis identified four subgroups among the patients with SLE. Subgroup 1 (29.3%) was dominated by anti-SSA/Ro60/Ro52/SSB autoantibodies and was strongly associated with HLA-DRB1*03 (odds ratio [OR] = 4.73; 95% confidence interval [CI] = 4.52-4.94). Discoid lesions were more common for this disease subgroup (OR = 1.71, 95% CI = 1.18-2.47). Subgroup 2 (28.7%) was dominated by anti-nucleosome/SmRNP/DNA/RNPA autoantibodies and associated with HLA-DRB1*15 (OR = 1.62, 95% CI = 1.41-1.84). Nephritis was most common in this subgroup (OR = 1.61, 95% CI = 1.14-2.26). Subgroup 3 (23.8%) was characterized by anti-ß2 GPI-IgG/anti-CL-IgG/IgM autoantibodies and a higher frequency of HLA-DRB1*04 compared with the other patients with SLE. Vascular events were more common in Subgroup 3 (OR = 1.74, 95% CI = 1.2-2.5). Subgroup 4 (18.2%) was negative for the investigated autoantibodies, and this subgroup was not associated with HLA-DRB1. Additionally, the levels of eight cytokines significantly differed among the disease subgroups. CONCLUSION: Our findings suggest that four fairly distinct subgroups can be identified on the basis of the autoantibody profile in SLE. These four SLE subgroups differ regarding associations with HLA-DRB1 alleles and immunological and clinical features, suggesting dissimilar disease pathways.

Understanding interactions between risk factors, and assessing the utility of the additive and multiplicative models through simulations.

Understanding the genetic background of complex diseases requires the expansion of studies beyond univariate associations. Therefore, it is important to use interaction assessments of risk factors in order to discover whether, and how genetic risk variants act together on disease development. The principle of interaction analysis is to explore the magnitude of the combined effect of risk factors on disease causation. In this study, we use simulations to investigate different scenarios of causation to show how the magnitude of the effect of two risk factors interact. We mainly focus on the two most commonly used interaction models, the additive and multiplicative risk scales, since there is often confusion regarding their use and interpretation. Our results show that the combined effect is multiplicative when two risk factors are involved in the same chain of events, an interaction called synergism. Synergism is often described as a deviation from additivity, which is a broader term. Our results also confirm that it is often relevant to estimate additive effect relationships, because they correspond to independent risk factors at low disease prevalence. Importantly, we evaluate the threshold of more than two required risk factors for disease causation, called the multifactorial threshold model. We found a simple mathematical relationship (square root) between the threshold and an additive-to-multiplicative linear effect scale (AMLES), where 0 corresponds to an additive effect and 1 to a multiplicative. We propose AMLES as a metric that could be used to test different effects relationships at the same time, given that it can simultaneously reveal additive, multiplicative and intermediate risk effects relationships. Finally, the utility of our simulation study was demonstrated using real data by analyzing and interpreting gene-gene interaction odds ratios from a rheumatoid arthritis case-control cohort.

The spectrum of association in HLA region with rheumatoid arthritis in a diverse Asian population: evidence from the MyEIRA case-control study.

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, PGWAS = 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, PGWAS = 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, PGWAS = 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.

A Gene-Environment Interaction Between Smoking and Gene polymorphisms Provides a High Risk of Two Subgroups of Sarcoidosis.

The influence and effect of cigarette smoking in sarcoidosis is unclear. Here, we evaluated gene-environment interaction between multiple genetic variants including HLA genes and smoking in sarcoidosis defined by two clinical phenotypes, Löfgren's syndrome (LS) and patients without Löfgren's syndrome (non-LS). To quantify smoking effects in sarcoidosis, we performed a gene-environment interaction study in a Swedish population-based case-control study consisting of 3,713 individuals. Cases and controls were classified according to their cigarette smoking status and genotypes by Immunochip platform. Gene-smoking interactions were quantified by an additive interaction model using a logistic regression adjusted by sex, age and first two principal components. The estimated attributable proportion (AP) was used to quantify the interaction effect. Assessment of smoking effects with inclusion of genetic information revealed 53 (in LS) and 34 (in non-LS) SNP-smoking additive interactions at false discovery rate (FDR) below 5%. The lead signals interacting with smoking were rs12132140 (AP = 0.56, 95% CI = 0.22-0.90), p = 1.28e-03) in FCRL1 for LS and rs61780312 (AP = 0.62, 95% CI = 0.28-0.90), p = 3e-04) in IL23R for non-LS. We further identified 16 genomic loci (in LS) and 13 (in non-LS) that interact with cigarette smoking. These findings suggest that sarcoidosis risk is modulated by smoking due to genetic susceptibility. Therefore, patients having certain gene variants, are at a higher risk for the disease. Consideration of individual's genetic predisposition is crucial to quantify effects of smoking in sarcoidosis.

Systematic approach demonstrates enrichment of multiple interactions between non-HLA risk variants and HLA-DRB1 risk alleles in rheumatoid arthritis.

OBJECTIVE: In anti-citrullinated protein antibody positive rheumatoid arthritis (ACPA-positive RA), a particular subset of HLA-DRB1 alleles, called shared epitope (SE) alleles, is a highly influential genetic risk factor. Here, we investigated whether non-HLA single nucleotide polymorphisms (SNP), conferring low disease risk on their own, interact with SE alleles more frequently than expected by chance and if such genetic interactions influence the HLA-DRB1 SE effect concerning risk to ACPA-positive RA. METHODS: We computed the attributable proportion (AP) due to additive interaction at genome-wide level for two independent ACPA-positive RA cohorts: the Swedish epidemiological investigation of rheumatoid arthritis (EIRA) and the North American rheumatoid arthritis consortium (NARAC). Then, we tested for differences in the AP p value distributions observed for two groups of SNPs, non-associated and associated with disease. We also evaluated whether the SNPs in interaction with HLA-DRB1 were cis-eQTLs in the SE alleles context in peripheral blood mononuclear cells from patients with ACPA-positive RA (SE-eQTLs). RESULTS: We found a strong enrichment of significant interactions (AP p<0.05) between the HLA-DRB1 SE alleles and the group of SNPs associated with ACPA-positive RA in both cohorts (Kolmogorov-Smirnov test D=0.35 for EIRA and D=0.25 for NARAC, p<2.2e-16 for both). Interestingly, 564 out of 1492 SNPs in consistent interaction for both cohorts were significant SE-eQTLs. Finally, we observed that the effect size of HLA-DRB1 SE alleles for disease decreases from 5.2 to 2.5 after removal of the risk alleles of the two top interacting SNPs (rs2476601 and rs10739581). CONCLUSION: Our data demonstrate that there are massive genetic interactions between the HLA-DRB1 SE alleles and non-HLA genetic variants in ACPA-positive RA.

EOMES-positive CD4+ T cells are increased in PTPN22 (1858T) risk allele carriers.

The presence of the PTPN22 risk allele (1858T) is associated with several autoimmune diseases including rheumatoid arthritis (RA). Despite a number of studies exploring the function of PTPN22 in T cells, the exact impact of the PTPN22 risk allele on T-cell function in humans is still unclear. In this study, using RNA sequencing, we show that, upon TCR-activation, naïve human CD4+ T cells homozygous for the PTPN22 risk allele overexpress a set of genes including CFLAR and 4-1BB, which are important for cytotoxic T-cell differentiation. Moreover, the protein expression of the T-box transcription factor Eomesodermin (EOMES) was increased in T cells from healthy donors homozygous for the PTPN22 risk allele and correlated with a decreased number of naïve CD4+ T cells. There was no difference in the frequency of other CD4+ T-cell subsets (Th1, Th17, Tfh, Treg). Finally, an accumulation of EOMES+ CD4+ T cells was observed in synovial fluid of RA patients with a more pronounced production of Perforin-1 in PTPN22 risk allele carriers. Altogether, we propose a novel mechanism of action of PTPN22 risk allele through the generation of cytotoxic CD4+ T cells and identify EOMES+ CD4+ T cells as a relevant T-cell subset in RA pathogenesis.

Integration of known DNA, RNA and protein biomarkers provides prediction of anti-TNF response in rheumatoid arthritis: results from the COMBINE study.

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

Development of autoantibodies against muscle-specific FHL1 in severe inflammatory myopathies.

Mutations of the gene encoding four-and-a-half LIM domain 1 (FHL1) are the causative factor of several X-linked hereditary myopathies that are collectively termed FHL1-related myopathies. These disorders are characterized by severe muscle dysfunction and damage. Here, we have shown that patients with idiopathic inflammatory myopathies (IIMs) develop autoimmunity to FHL1, which is a muscle-specific protein. Anti-FHL1 autoantibodies were detected in 25% of IIM patients, while patients with other autoimmune diseases or muscular dystrophies were largely anti-FHL1 negative. Anti-FHL1 reactivity was predictive for muscle atrophy, dysphagia, pronounced muscle fiber damage, and vasculitis. FHL1 showed an altered expression pattern, with focal accumulation in the muscle fibers of autoantibody-positive patients compared with a homogeneous expression in anti-FHL1-negative patients and healthy controls. We determined that FHL1 is a target of the cytotoxic protease granzyme B, indicating that the generation of FHL1 fragments may initiate FHL1 autoimmunity. Moreover, immunization of myositis-prone mice with FHL1 aggravated muscle weakness and increased mortality, suggesting a direct link between anti-FHL1 responses and muscle damage. Together, our findings provide evidence that FHL1 may be involved in the pathogenesis not only of genetic FHL1-related myopathies but also of autoimmune IIM. Importantly, these results indicate that anti-FHL1 autoantibodies in peripheral blood have promising potential as a biomarker to identify a subset of severe IIM.

Gene-gene interaction and RNA splicing profiles of MAP2K4 gene in rheumatoid arthritis.

We performed gene-gene interaction analysis, with HLA-DRB1 shared epitope (SE) alleles for 195 SNPs within immunologically important MAP2K, MAP3K and MAP4K gene families, in 2010 rheumatoid arthritis (RA) patients and 2280 healthy controls. We found a significant statistical interaction for rs10468473 with SE alleles in autoantibody-positive RA. Individuals heterozygous for rs10468473 demonstrated higher expression of total MAP2K4 mRNA in blood, compared to A-allele homozygous. We discovered a novel, putatively translated, "cassette exon" RNA splice form of MAP2K4, differentially expressed in peripheral blood mononuclear cells from 88 RA cases and controls. Within the group of RA patients, we observed a correlation of MAP2K4 isoform expression with carried SE alleles, autoantibody, and rheumatoid factor profiles. TNF-dependent modulation of isoform expression pattern was detected in the Jurkat cell line. Our data suggest a genetic interaction between MAP2K4 and HLA-DRB1, and the importance of rs10468473 and MAP2K4 splice variants in the development of autoantibody-positive RA.

A systemic sclerosis and systemic lupus erythematosus pan-meta-GWAS reveals new shared susceptibility loci.

Systemic sclerosis (SSc) and systemic lupus erythematosus (SLE) are two archetypal systemic autoimmune diseases which have been shown to share multiple genetic susceptibility loci. In order to gain insight into the genetic basis of these diseases, we performed a pan-meta-analysis of two genome-wide association studies (GWASs) together with a replication stage including additional SSc and SLE cohorts. This increased the sample size to a total of 21,109 (6835 cases and 14,274 controls). We selected for replication 19 SNPs from the GWAS data. We were able to validate KIAA0319L (P = 3.31 × 10(-11), OR = 1.49) as novel susceptibility loci for SSc and SLE. Furthermore, we also determined that the previously described SLE susceptibility loci PXK (P = 3.27 × 10(-11), OR = 1.20) and JAZF1 (P = 1.11 × 10(-8), OR = 1.13) are shared with SSc. Supporting these new discoveries, we observed that KIAA0319L was overexpressed in peripheral blood cells of SSc and SLE patients compared with healthy controls. With these, we add three (KIAA0319L, PXK and JAZF1) and one (KIAA0319L) new susceptibility loci for SSc and SLE, respectively, increasing significantly the knowledge of the genetic basis of autoimmunity.

Evidence of new risk genetic factor to systemic lupus erythematosus: the UBASH3A gene.

The ubiquitin associated and Src-homology 3 (SH3) domain containing A (UBASH3a) is a suppressor of T-cell receptor signaling, underscoring antigen presentation to T-cells as a critical shared mechanism of diseases pathogenesis. The aim of the present study was to determine whether the UBASH3a gene influence the susceptibility to systemic lupus erythematosus (SLE) in Caucasian populations. We evaluated five UBASH3a polymorphisms (rs2277798, rs2277800, rs9976767, rs13048049 and rs17114930), using TaqMan® allelic discrimination assays, in a discovery cohort that included 906 SLE patients and 1165 healthy controls from Spain. The SNPs that exhibit statistical significance difference were evaluated in a German replication cohort of 360 SLE patients and 379 healthy controls. The case-control analysis in the Spanish population showed a significant association between the rs9976767 and SLE (Pc = 9.9E-03 OR = 1.21 95%CI = 1.07-1.37) and a trend of association for the rs2277798 analysis (P = 0.09 OR = 0.9 95%CI = 0.79-1.02). The replication in a German cohort and the meta-analysis confirmed that the rs9976767 (Pc = 0.02; Pc = 2.4E-04, for German cohort and meta-analysis, respectively) and rs2277798 (Pc = 0.013; Pc = 4.7E-03, for German cohort and meta-analysis, respectively) UBASH3a variants are susceptibility factors for SLE. Finally, a conditional regression analysis suggested that the most likely genetic variation responsible for the association was the rs9976767 polymorphism. Our results suggest that UBASH3a gene plays a role in the susceptibility to SLE. Moreover, our study indicates that UBASH3a can be considered as a common genetic factor in autoimmune diseases.

New insight on the Xq28 association with systemic sclerosis.

OBJECTIVE: To evaluate whether the systemic sclerosis (SSc)-associated IRAK1 non-synonymous single-nucleotide polymorphism rs1059702 is responsible for the Xq28 association with SSc or whether there are other independent signals in the nearby methyl-CpG-binding protein 2 gene (MECP2). METHODS: We analysed a total of 3065 women with SSc and 2630 unaffected controls from five independent Caucasian cohorts. Four tag single-nucleotide polymorphisms of MECP2 (rs3027935, rs17435, rs5987201 and rs5945175) and the IRAK1 variant rs1059702 were genotyped using TaqMan predesigned assays. A meta-analysis including all cohorts was performed to test the overall effect of these Xq28 polymorphisms on SSc. RESULTS: IRAK1 rs1059702 and MECP2 rs17435 were associated specifically with diffuse cutaneous SSc (PFDR=4.12×10(-3), OR=1.27, 95% CI 1.09 to 1.47, and PFDR=5.26×10(-4), OR=1.30, 95% CI 1.14 to 1.48, respectively), but conditional logistic regression analysis showed that the association of IRAK1 rs1059702 with this subtype was explained by that of MECP2 rs17435. On the other hand, IRAK1 rs1059702 was consistently associated with presence of pulmonary fibrosis (PF), because statistical significance was observed when comparing SSc patients PF+ versus controls (PFDR=0.039, OR=1.30, 95% CI 1.07 to 1.58) and SSc patients PF+ versus SSc patients PF- (p=0.025, OR=1.26, 95% CI 1.03 to 1.55). CONCLUSIONS: Our data clearly suggest the existence of two independent signals within the Xq28 region, one located in IRAK1 related to PF and another in MECP2 related to diffuse cutaneous SSc, indicating that both genes may have an impact on the clinical outcome of the disease.