Genetics of Systemic Sclerosis: An Update

Systemic sclerosis (SSc) is an autoimmune disease characterized by vasculopathy, immune cell activation, and fibrosis of the skin and internal organs. Over the past few years, a role for genetics in the susceptibility for SSc has been established. This review aims to provide an update on the progress made in the past year or so within the field of SSc genetics research. This year has been of particular interest due to the publication of a large genome-wide association study, further investigations into gene–gene interactions, and the tendency to validate genetic results in functional models

The systemic lupus erythematosus IRF5 risk haplotype is associated with systemic sclerosis

Systemic sclerosis (SSc) is a fibrotic autoimmune disease in which the genetic component plays an important role. One of the strongest SSc association signals outside the human leukocyte antigen (HLA) region corresponds to interferon (IFN) regulatory factor 5 (IRF5), a major regulator of the type I IFN pathway. In this study we aimed to evaluate whether three different haplotypic blocks within this locus, which have been shown to alter the protein function influencing systemic lupus erythematosus (SLE) susceptibility, are involved in SSc susceptibility and clinical phenotypes. For that purpose, we genotyped one representative single-nucleotide polymorphism (SNP) of each block (rs10488631, rs2004640, and rs4728142) in a total of 3,361 SSc patients and 4,012 unaffected controls of Caucasian origin from Spain, Germany, The Netherlands, Italy and United Kingdom. A meta-analysis of the allele frequencies was performed to analyse the overall effect of these IRF5 genetic variants on SSc. Allelic combination and dependency tests were also carried out. The three SNPs showed strong associations with the global disease (rs4728142: P = 1.34×10<sup>−8</sup>, OR = 1.22, CI 95% = 1.14–1.30; rs2004640: P = 4.60×10<sup>−7</sup>, OR = 0.84, CI 95% = 0.78–0.90; rs10488631: P = 7.53×10<sup>−20</sup>, OR = 1.63, CI 95% = 1.47–1.81). However, the association of rs2004640 with SSc was not independent of rs4728142 (conditioned P = 0.598). The haplotype containing the risk alleles (rs4728142*A-rs2004640*T-rs10488631*C: P = 9.04×10<sup>−22</sup>, OR = 1.75, CI 95% = 1.56–1.97) better explained the observed association (likelihood P-value = 1.48×10<sup>−4</sup>), suggesting an additive effect of the three haplotypic blocks. No statistical significance was observed in the comparisons amongst SSc patients with and without the main clinical characteristics. Our data clearly indicate that the SLE risk haplotype also influences SSc predisposition, and that this association is not sub-phenotype-specific

Genome-Wide Scan Identifies TNIP1, PSORS1C1, and RHOB as Novel Risk Loci for Systemic Sclerosis

Systemic sclerosis (SSc) is an orphan, complex, inflammatory disease affecting the immune system and connective tissue. SSc stands out as a severely incapacitating and life-threatening inflammatory rheumatic disease, with a largely unknown pathogenesis. We have designed a two-stage genome-wide association study of SSc using case-control samples from France, Italy, Germany, and Northern Europe. The initial genome-wide scan was conducted in a French post quality-control sample of 564 cases and 1,776 controls, using almost 500 K SNPs. Two SNPs from the MHC region, together with the 6 loci outside MHC having at least one SNP with a P<10−5 were selected for follow-up analysis. These markers were genotyped in a post-QC replication sample of 1,682 SSc cases and 3,926 controls. The three top SNPs are in strong linkage disequilibrium and located on 6p21, in the HLA-DQB1 gene: rs9275224, P = 9.18×10−8, OR = 0.69, 95% CI [0.60–0.79]; rs6457617, P = 1.14×10−7 and rs9275245, P = 1.39×10−7. Within the MHC region, the next most associated SNP (rs3130573, P = 1.86×10−5, OR = 1.36 [1.18–1.56]) is located in the PSORS1C1 gene. Outside the MHC region, our GWAS analysis revealed 7 top SNPs (P<10−5) that spanned 6 independent genomic regions. Follow-up of the 17 top SNPs in an independent sample of 1,682 SSc and 3,926 controls showed associations at PSORS1C1 (overall P = 5.70×10−10, OR:1.25), TNIP1 (P = 4.68×10−9, OR:1.31), and RHOB loci (P = 3.17×10−6, OR:1.21). Because of its biological relevance, and previous reports of genetic association at this locus with connective tissue disorders, we investigated TNIP1 expression. A markedly reduced expression of the TNIP1 gene and also its protein product were observed both in lesional skin tissue and in cultured dermal fibroblasts from SSc patients. Furthermore, TNIP1 showed in vitro inhibitory effects on inflammatory cytokine-induced collagen production. The genetic signal of association with TNIP1 variants, together with tissular and cellular investigations, suggests that this pathway has a critical role in regulating autoimmunity and SSc pathogenesis

Associations of baseline use of biologic or targeted synthetic DMARDs with COVID-19 severity in rheumatoid arthritis: Results from the COVID-19 Global Rheumatology Alliance physician registry.

OBJECTIVE: To investigate baseline use of biologic or targeted synthetic (b/ts) disease-modifying antirheumatic drugs (DMARDs) and COVID-19 outcomes in rheumatoid arthritis (RA). METHODS: We analysed the COVID-19 Global Rheumatology Alliance physician registry (from 24 March 2020 to 12 April 2021). We investigated b/tsDMARD use for RA at the clinical onset of COVID-19 (baseline): abatacept (ABA), rituximab (RTX), Janus kinase inhibitors (JAKi), interleukin 6 inhibitors (IL-6i) or tumour necrosis factor inhibitors (TNFi, reference group). The ordinal COVID-19 severity outcome was (1) no hospitalisation, (2) hospitalisation without oxygen, (3) hospitalisation with oxygen/ventilation or (4) death. We used ordinal logistic regression to estimate the OR (odds of being one level higher on the ordinal outcome) for each drug class compared with TNFi, adjusting for potential baseline confounders. RESULTS: Of 2869 people with RA (mean age 56.7 years, 80.8% female) on b/tsDMARD at the onset of COVID-19, there were 237 on ABA, 364 on RTX, 317 on IL-6i, 563 on JAKi and 1388 on TNFi. Overall, 613 (21%) were hospitalised and 157 (5.5%) died. RTX (OR 4.15, 95% CI 3.16 to 5.44) and JAKi (OR 2.06, 95% CI 1.60 to 2.65) were each associated with worse COVID-19 severity compared with TNFi. There were no associations between ABA or IL6i and COVID-19 severity. CONCLUSIONS: People with RA treated with RTX or JAKi had worse COVID-19 severity than those on TNFi. The strong association of RTX and JAKi use with poor COVID-19 outcomes highlights prioritisation of risk mitigation strategies for these people

Integration of sequence data from a consanguineous family with genetic data from an outbred population identifies PLB1 as a candidate rheumatoid arthritis risk gene

Integrating genetic data from families with highly penetrant forms of disease together with genetic data from outbred populations represents a promising strategy to uncover the complete frequency spectrum of risk alleles for complex traits such as rheumatoid arthritis (RA). Here, we demonstrate that rare, low-frequency and common alleles at one gene locus, phospholipase B1 (PLB1), might contribute to risk of RA in a 4-generation consanguineous pedigree (Middle Eastern ancestry) and also in unrelated individuals from the general population (European ancestry). Through identity-by-descent (IBD) mapping and whole-exome sequencing, we identified a non-synonymous c.2263G>C (p.G755R) mutation at the PLB1 gene on 2q23, which significantly co-segregated with RA in family members with a dominant mode of inheritance (P = 0.009). We further evaluated PLB1 variants and risk of RA using a GWAS meta-analysis of 8,875 RA cases and 29,367 controls of European ancestry. We identified significant contributions of two independent non-coding variants near PLB1 with risk of RA (rs116018341 [MAF = 0.042] and rs116541814 [MAF = 0.021], combined P = 3.2×10-6). Finally, we performed deep exon sequencing of PLB1 in 1,088 RA cases and 1,088 controls (European ancestry), and identified suggestive dispersion of rare protein-coding variant frequencies between cases and controls (P = 0.049 for C-alpha test and P = 0.055 for SKAT). Together, these data suggest that PLB1 is a candidate risk gene for RA. Future studies to characterize the full spectrum of genetic risk in the PLB1 genetic locus are warranted. © 2014 Plenge et al

Genetics of rheumatoid arthritis contributes to biology and drug discovery

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

Guidelines for the use and interpretation of assays for monitoring autophagy (4th edition)1.

In 2008, we published the first set of guidelines for standardizing research in autophagy. Since then, this topic has received increasing attention, and many scientists have entered the field. Our knowledge base and relevant new technologies have also been expanding. Thus, it is important to formulate on a regular basis updated guidelines for monitoring autophagy in different organisms. Despite numerous reviews, there continues to be confusion regarding acceptable methods to evaluate autophagy, especially in multicellular eukaryotes. Here, we present a set of guidelines for investigators to select and interpret methods to examine autophagy and related processes, and for reviewers to provide realistic and reasonable critiques of reports that are focused on these processes. These guidelines are not meant to be a dogmatic set of rules, because the appropriateness of any assay largely depends on the question being asked and the system being used. Moreover, no individual assay is perfect for every situation, calling for the use of multiple techniques to properly monitor autophagy in each experimental setting. Finally, several core components of the autophagy machinery have been implicated in distinct autophagic processes (canonical and noncanonical autophagy), implying that genetic approaches to block autophagy should rely on targeting two or more autophagy-related genes that ideally participate in distinct steps of the pathway. Along similar lines, because multiple proteins involved in autophagy also regulate other cellular pathways including apoptosis, not all of them can be used as a specific marker for bona fide autophagic responses. Here, we critically discuss current methods of assessing autophagy and the information they can, or cannot, provide. Our ultimate goal is to encourage intellectual and technical innovation in the field

Value of ultrasonography as a marker of early response to abatacept in patients with rheumatoid arthritis and an inadequate response to methotrexate: results from the APPRAISE study

Objectives: To study the responsiveness of a combined power Doppler and greyscale ultrasound (PDUS) score for assessing synovitis in biologic-naïve patients with rheumatoid arthritis (RA) starting abatacept plus methotrexate (MTX). Methods: In this open-label, multicentre, single-arm study, patients with RA (MTX inadequate responders) received intravenous abatacept (∼10 mg/kg) plus MTX for 24 weeks. A composite PDUS synovitis score, developed by the Outcome Measures in Rheumatology–European League Against Rheumatism (OMERACT–EULAR)-Ultrasound Task Force, was used to evaluate individual joints. The maximal score of each joint was added into a Global OMERACT–EULAR Synovitis Score (GLOESS) for bilateral metacarpophalangeal joints (MCPs) 2–5 (primary objective). The value of GLOESS containing other joint sets was explored, along with clinical efficacy. Results: Eighty-nine patients completed the 24-week treatment period. The earliest PDUS sign of improvement in synovitis was at week 1 (mean change in GLOESS (MCPs 2–5): −0.7 (95% CIs −1.2 to −0.1)), with continuous improvement to week 24. Early improvement was observed in the component scores (power Doppler signal at week 1, synovial hyperplasia at week 2, joint effusion at week 4). Comparable changes were observed for 22 paired joints and minimal joint subsets. Mean Disease Activity Score 28 (C reactive protein) was significantly reduced from weeks 1 to 24, reaching clinical meaningful improvement (change ≥1.2) at week 8. Conclusions: In this first international prospective study, the composite PDUS score is responsive to abatacept. GLOESS demonstrated the rapid onset of action of abatacept, regardless of the number of joints examined. Ultrasound is an objective tool to monitor patients with RA under treatment. Trial registration number: NCT00767325