Epigenetic targets of Janus kinase inhibitors are linked to genetic risks of rheumatoid arthritis.

BACKGROUND: Current strategies that target cytokines (e.g., tumor necrosis factor (TNF)-α), or signaling molecules (e.g., Janus kinase (JAK)) have advanced the management for allergies and autoimmune diseases. Nevertheless, the molecular mechanism that underpins its clinical efficacy have largely remained elusive, especially in the local tissue environment. Here, we aimed to identify the genetic, epigenetic, and immunological targets of JAK inhibitors (JAKis), focusing on their effects on synovial fibroblasts (SFs), the major local effectors associated with destructive joint inflammation in rheumatoid arthritis (RA). METHODS: SFs were activated by cytokines related to inflammation in RA, and were treated with three types of JAKis or a TNF-α inhibitor (TNFi). Dynamic changes in transcriptome and chromatin accessibility were profiled across samples to identify drug targets. Furthermore, the putative targets were validated using luciferase assays and clustered regularly interspaced short palindromic repeat (CRISPR)-based genome editing. RESULTS: We found that both JAKis and the TNFi targeted the inflammatory module including IL6. Conversely, specific gene signatures that were preferentially inhibited by either of the drug classes were identified. Strikingly, RA risk enhancers for CD40 and TRAF1 were distinctively regulated by JAKis and the TNFi. We performed luciferase assays and CRISPR-based genome editing, and successfully fine-mapped the single causal variants in these loci, rs6074022-CD40 and rs7021049-TRAF1. CONCLUSIONS: JAKis and the TNFi had a direct impact on different RA risk enhancers, and we identified nucleotide-resolution targets for both drugs. Distinctive targets of clinically effective drugs could be useful for tailoring the application of these drugs and future design of more efficient treatment strategies.

Systematic review and meta-analysis for the 2024 update of the Japan College of Rheumatology clinical practice guidelines for the management of rheumatoid arthritis.

OBJECTIVES: To update evidence on the efficacy and safety of disease-modifying antirheumatic drugs (DMARDs) and provide information to the taskforce for the 2024 update of the Japan College of Rheumatology (JCR) clinical practice guidelines (CPG) for the management of rheumatoid arthritis (RA). METHODS: We searched various databases for randomised controlled trials on RA published until June 2022, with no language restriction. For each of the 15 clinical questions, 2 independent reviewers screened the articles, evaluated the core outcomes, and performed meta-analyses. RESULTS: Subcutaneous injection of methotrexate (MTX) showed similar efficacy to oral MTX in MTX-naïve RA patients. Ozoralizumab combined with MTX improved drug efficacy compared to the placebo in RA patients with inadequate response (IR) to csDMARD. Rituximab with and without concomitant csDMARDs showed similar efficacy to other bDMARDs in bDMARD-IR RA patients. Combined Janus kinase inhibitors and MTX achieved similar clinical responses and equal safety during a 4-year period compared to tumour necrosis factor inhibitors in MTX-IR RA patients. Biosimilars showed efficacy equivalent to that of the original bDMARDs in csDMARD-IR and bDMARD-IR RA patients. CONCLUSION: This systematic review provides latest evidence for the 2024 update of the JCR CPG for RA management.

Case report: Inflammatory sternoclavicular joint arthritis induced by an immune checkpoint inhibitor with remarkable responsiveness to infliximab.

This report describes the case of a 48-year-old woman who presented with sternoclavicular joint arthritis after administration of an immune checkpoint inhibitor (ICI), durvalumab, for small cell lung carcinoma. The onset of arthritis transpired 18 months after the commencement of the ICI therapeutic regimen and demonstrated resilience to glucocorticoid treatment. After excluding infectious aetiologies and metastatic involvement, the patient was diagnosed with ICI-induced arthritis (ICI-IA). Considering the articular implications akin to the SAPHO syndrome, the patient was treated with infliximab, resulting in complete resolution. This finding implies that biological DMARDs can serve as effective interventions for ICI-induced sternoclavicular joint arthritis. Given the heterogeneous nature of its pathogenesis, the selection of therapeutic agents may require customization based on the distinct clinical presentation of each individual case.

A 52-Year-Old With Painful Fingertips.

A patient had cold, cyanotic fingertips with small ulcerations. Laboratory testing showed leukocytosis, eosinophilia, and elevated D-dimer level; results of tests for antinuclear antibodies, antiphospholipid antibodies, antineutrophil cytoplasmic antibodies against proteinase 3 and myeloperoxidase, and HIV were negative. What is the diagnosis and what would you do next?

Immunomics analysis of rheumatoid arthritis identified precursor dendritic cells as a key cell subset of treatment resistance.

OBJECTIVES: Little is known about the immunology underlying variable treatment response in rheumatoid arthritis (RA). We performed large-scale transcriptome analyses of peripheral blood immune cell subsets to identify immune cells that predict treatment resistance. METHODS: We isolated 18 peripheral blood immune cell subsets of 55 patients with RA requiring addition of new treatment and 39 healthy controls, and performed RNA sequencing. Transcriptome changes in RA and treatment effects were systematically characterised. Association between immune cell gene modules and treatment resistance was evaluated. We validated predictive value of identified parameters for treatment resistance using quantitative PCR (qPCR) and mass cytometric analysis cohorts. We also characterised the identified population by synovial single cell RNA-sequencing analysis. RESULTS: Immune cells of patients with RA were characterised by enhanced interferon and IL6-JAK-STAT3 signalling that demonstrate partial normalisation after treatment. A gene expression module of plasmacytoid dendritic cells (pDC) reflecting the expansion of dendritic cell precursors (pre-DC) exhibited strongest association with treatment resistance. Type I interferon signalling was negatively correlated to pre-DC gene expression. qPCR and mass cytometric analysis in independent cohorts validated that the pre-DC associated gene expression and the proportion of pre-DC were significantly higher before treatment in treatment-resistant patients. A cluster of synovial DCs showed both features of pre-DC and pro-inflammatory conventional DC2s. CONCLUSIONS: An increase in pre-DC in peripheral blood predicted RA treatment resistance. Pre-DC could have pathophysiological relevance to RA treatment response.

Anti-synthetase Syndrome That Relapsed with Pulmonary Arterial Hypertension and Malignancy.

A 69-year-old man with a history of anti-synthetase antibody-positive polymyositis and interstitial lung disease (ILD) stable for more than 20 years suddenly developed pulmonary artery hypertension (PAH) with a mean PA pressure of 46 mmHg. At this stage, ILD was mild, but it became acutely exacerbated later, and high-dose corticosteroid and intravenous cyclophosphamide ameliorated both PAH and ILD. The tricuspid regurgitation pressure gradient decreased from 80 to 49 mmHg and ILD recovered almost completely. During a systemic examination, bone metastatic cancer of unknown origin was found. We herein report the relationship between anti-synthetase syndrome (ASS) and PAH as well as ASS and malignancy.

Transcriptome Profiling of Immune Cell Types in Peripheral Blood Reveals Common and Specific Pathways Involved in the Pathogenesis of Myositis-Specific Antibody-Positive Inflammatory Myopathies.

OBJECTIVE: Idiopathic inflammatory myopathies (IIM) demonstrate characteristic clinical phenotypes depending on the myositis-specific antibody (MSAs) present. We aimed to identify common or MSA-specific immunological pathways in different immune cell types from peripheral blood by transcriptome analysis. METHODS: We recruited 33 patients with IIM who were separated into the following groups: 15 patients with active disease at onset and 18 with inactive disease under treatment. All patients were positive for MSAs: anti-melanoma differentiation-associated gene 5 (MDA5) antibody (Ab) in 10 patients, anti-Mi-2 Ab in 7, and anti-aminoacyl-transfer RNA synthetase (ARS) Ab in 16. The patients were compared with 33 healthy controls. Twenty-four immune cell types sorted from peripheral blood were analyzed by flow cytometry, RNA sequencing, and differentially expressed gene analysis combined with pathway analysis. RESULTS: The frequencies of memory B cell types were significantly decreased in active patients, and the frequency of plasmablasts was prominently increased in active patients with anti-MDA5 Ab in comparison with healthy controls. The expression of type I interferon (IFN)-stimulated genes of all immune cell types was increased in the active, but not inactive, patients. Endoplasmic reticulum stress-related genes in all IIM memory B cells and oxidative phosphorylation-related genes in inactive IIM double negative B cells were also increased, suggesting prominent B cell activation in IIM. Furthermore, active patients with anti-MDA5 Ab, anti-Mi-2 Ab, or anti-ARS Ab were distinguished by IFN-stimulated and oxidative phosphorylation-related gene expression in plasmablasts. CONCLUSION: Unique gene expression patterns in patients with IIM with different disease activity levels and MSA types suggest different pathophysiologies. Especially, B cells may contribute to common and MSA-specific immunological pathways in IIM.

A case of aortoduodenal fistula caused by IgG4-related periaortitis.

An 86-year-old man who underwent endovascular aortic repair for impending rupture of an abdominal aortic aneurysm a year ago presented to our hospital because of fatigue and black stools. Multiple bacterial specimens were detected in blood cultures, and computed tomography following oral administration of gastrografin demonstrated gastrografin in the abdominal aorta. The diagnosis of aortic duodenal fistula was confirmed and emergency abdominal aortic replacement was performed. The pathological findings of the aorta included a large number of immunoglobulin G4 (IgG4)-positive plasma cells infiltrating all layers of the aortic wall, with particularly marked thickening of the adventitia. The serum IgG4 level was 241 mg/dl and IgG4-related periaortitis was diagnosed. Aortoduodenal fistula is a rare but fatal complication of IgG4-related periaortitis. Patients should be followed carefully after endovascular aortic repair for inflammatory abdominal aortic aneurysms.

Four cases of dermatomyositis with abnormally high anti-MDA-5 antibody titres and not high levels of serum ferritin.

Anti-melanoma differentiation-associated gene 5 (MDA-5) antibody-positive dermatomyositis is a fatal disease presenting with rapidly progressive interstitial lung disease. High ferritin levels are a well-known poor prognostic factor. A high anti-MDA-5 antibody titre was also recently identified as a poor prognostic factor. We encountered four cases that had extremely high anti-MDA-5 antibody titres without high levels of ferritin in the initial examination. All cases were female with ages ranging between 29 and 54 years (mean age, 44 years). In the initial examination, anti-MDA-5 antibody titres were 2060-3040 (normal range, <32 index), ferritin levels were 87-480 ng/ml (normal range, 2.6-129.4 ng/ml), KL-6 level was 186-1806 U/ml (normal range, <500 U/ml), and creatine kinase level was normal in all patients. One patient had respiratory distress on exertion. Computed Tomography (CT) images showed mild ground-glass attenuation/reticular shadows near the pleura in all patients. Three patients were treated with a combination of high-dose glucocorticoids, intermittent intravenous cyclophosphamide, and calcineurin inhibitors, and two required plasma exchange due to the worsening of lung lesion. In these patients, ferritin and KL-6 levels tended to elevate after the beginning of treatment. Very mild pulmonary lesions disappeared in one patient treated with moderate doses of a glucocorticoid and calcineurin inhibitor. All patients survived, and one required oxygen on exertion at discharge. The condition of patients with abnormally high anti-MDA-5 antibody titres may deteriorate even though ferritin levels were not high and lung shadows are minimal at presentation. Therefore, intensive treatment needs to be considered early in the course of the disease regardless of the serum ferritin level.

Splicing QTL analysis focusing on coding sequences reveals mechanisms for disease susceptibility loci.

Splicing quantitative trait loci (sQTLs) are one of the major causal mechanisms in genome-wide association study (GWAS) loci, but their role in disease pathogenesis is poorly understood. One reason is the complexity of alternative splicing events producing many unknown isoforms. Here, we propose two approaches, namely integration and selection, for this complexity by focusing on protein-structure of isoforms. First, we integrate isoforms with the same coding sequence (CDS) and identify 369-601 integrated-isoform ratio QTLs (i2-rQTLs), which altered protein-structure, in six immune subsets. Second, we select CDS incomplete isoforms annotated in GENCODE and identify 175-337 isoform-ratio QTL (i-rQTL). By comprehensive long-read capture RNA-sequencing among these incomplete isoforms, we reveal 29 full-length isoforms with unannotated CDSs associated with GWAS traits. Furthermore, we show that disease-causal sQTL genes can be identified by evaluating their trans-eQTL effects. Our approaches highlight the understudied role of protein-altering sQTLs and are broadly applicable to other tissues and diseases.

Transcriptome analysis of immune cells from Behçet's syndrome patients: the importance of IL-17-producing cells and antigen-presenting cells in the pathogenesis of Behçet's syndrome.

BACKGROUND: Behçet's syndrome (BS) is an immune-mediated disease characterized by recurrent oral ulcers, genital ulcers, uveitis, and skin symptoms. HLA-B51, as well as other genetic polymorphisms, has been reported to be associated with BS; however, the pathogenesis of BS and its relationship to genetic risk factors still remain unclear. To address these points, we performed immunophenotyping and transcriptome analysis of immune cells from BS patients and healthy donors. METHODS: ImmuNexUT is a comprehensive database consisting of RNA sequencing data and eQTL database of immune cell subsets from patients with immune-mediated diseases and healthy donors, and flow cytometry data and transcriptome data from 23 BS patients and 28 healthy donors from the ImmuNexUT study were utilized for this study. Differential gene expression analysis and weighted gene co-expression network analysis (WGCNA) were performed to identify genes associated with BS and clinical features of BS. eQTL database was used to assess the relationship between genetic risk factors of BS with those genes. RESULTS: The frequency of Th17 cells was increased in BS patients, and transcriptome analysis of Th17 cells suggested the activation of the NFκB pathway in Th17 cells of BS patients. Next, WGCNA was used to group genes into modules with similar expression patterns in each subset. Modules of antigen-presenting cells were associated with BS, and pathway analysis suggested the activation of antigen-presenting cells of BS patients. Further examination of genes in BS-associated modules indicated that the expression of YBX3, a member of a plasmacytoid dendritic cell (pDC) gene module associated with BS, is influenced by a BS risk polymorphism, rs2617170, in pDCs, suggesting that YBX3 may be a key molecule connecting genetic risk factors of BS with disease pathogenesis. Furthermore, pathway analysis of modules associated with HLA-B51 indicated that the association of IL-17-associated pathways in memory CD8+ T cells with HLA-B51; therefore, IL-17-producing CD8+ T cells, Tc17 cells, may play a critical role in BS. CONCLUSIONS: Various cells including CD4+ T cells, CD8+ T cells, and antigen-presenting cells are important in the pathogenesis of BS. Tc17 cells and YBX3 may be potential therapeutic targets in BS.

Immune cell multiomics analysis reveals contribution of oxidative phosphorylation to B-cell functions and organ damage of lupus.

OBJECTIVE: Systemic lupus erythematosus (SLE) is the prototypical systemic autoimmune disease. While the long-term prognosis has greatly improved, better long-term survival is still necessary. The type I interferon (IFN) signature, a prominent feature of SLE, is not an ideal therapeutic target or outcome predictor. To explore immunological pathways in SLE more precisely, we performed transcriptomic, epigenomic and genomic analyses using 19 immune cell subsets from peripheral blood. METHODS: We sorted 19 immune cell subsets and identified the mRNA expression profiles and genetic polymorphisms in 107 patients with SLE and 92 healthy controls. Combined differentially expressed genes and expression quantitative trait loci analysis was conducted to find key driver genes in SLE pathogenesis. RESULTS: We found transcriptomic, epigenetic and genetic importance of oxidative phosphorylation (OXPHOS)/mitochondrial dysfunction in SLE memory B cells. Particularly, we identified an OXPHOS-regulating gene, PRDX6 (peroxiredoxin 6), as a key driver in SLE B cells. Prdx6-deficient B cells showed upregulated mitochondrial respiration as well as antibody production. We revealed OXPHOS signature was associated with type I IFN signalling-related genes (ISRGs) signature in SLE memory B cells. Furthermore, the gene sets related to innate immune signalling among ISRGs presented correlation with OXPHOS and these two signatures showed associations with SLE organ damage as well as specific clinical phenotypes. CONCLUSION: This work elucidated the potential prognostic marker for SLE. Since OXPHOS consists of the electron transport chain, a functional unit in mitochondria, these findings suggest the importance of mitochondrial dysfunction as a key immunological pathway involved in SLE.

Dynamic landscape of immune cell-specific gene regulation in immune-mediated diseases.

Genetic studies have revealed many variant loci that are associated with immune-mediated diseases. To elucidate the disease pathogenesis, it is essential to understand the function of these variants, especially under disease-associated conditions. Here, we performed a large-scale immune cell gene-expression analysis, together with whole-genome sequence analysis. Our dataset consists of 28 distinct immune cell subsets from 337 patients diagnosed with 10 categories of immune-mediated diseases and 79 healthy volunteers. Our dataset captured distinctive gene-expression profiles across immune cell types and diseases. Expression quantitative trait loci (eQTL) analysis revealed dynamic variations of eQTL effects in the context of immunological conditions, as well as cell types. These cell-type-specific and context-dependent eQTLs showed significant enrichment in immune disease-associated genetic variants, and they implicated the disease-relevant cell types, genes, and environment. This atlas deepens our understanding of the immunogenetic functions of disease-associated variants under in vivo disease conditions.

Identifying the most influential gene expression profile in distinguishing ANCA-associated vasculitis from healthy controls.

OBJECTIVE: Previous gene expression analyses seeking genes specific to antineutrophil cytoplasmic antibody-associated vasculitis (AAV) have been limited due to crude cell separation and the use of microarrays. This study aims to identify AAV-specific gene expression profiles in a way that overcomes those limitations. METHODS: Blood samples were collected from 26 AAV patients and 28 healthy controls (HCs). Neutrophils were isolated by negative selection, whereas 19 subsets of peripheral blood mononuclear cells were sorted by fluorescence assisted cell sorting. RNA-sequencing was then conducted for each sample, and iterative weighted gene correlation network analysis (iterativeWGCNA) and random forest were consecutively applied to identify the most influential gene module in distinguishing AAV from HCs. Correlations of the identified module with clinical parameters were evaluated, and the biological role was assessed with hub gene identification and pathway analysis. Particularly, the module's association with neutrophil extracellular trap formation, NETosis, was analyzed. Finally, the module's overlap with GWAS-identified autoimmune disease genes (GADGs) was assessed for validation. RESULTS: A neutrophil module (Neu_M20) was ranked top in the random forest analysis among 255 modules created by iterativeWGCNA. Neu_M20 correlated with disease activity and neutrophil counts but not with the presence of antineutrophil cytoplasmic antibody. The module comprised pro-inflammatory genes, including those related to NETosis, supported by experimental evidence. The genes in the module significantly overlapped GADGs. CONCLUSION: We identified the distinct group of pro-inflammatory genes in neutrophils, which characterize AAV. Further investigations are warranted to confirm our findings as they could serve as novel therapeutic targets.

Contribution of a European-Prevalent Variant near CD83 and an East Asian-Prevalent Variant near IL17RB to Herpes Zoster Risk in Tofacitinib Treatment: Results of Genome-Wide Association Study Meta-Analyses.

OBJECTIVE: Tofacitinib is an oral JAK inhibitor for the treatment of rheumatoid arthritis (RA), psoriatic arthritis, and ulcerative colitis, and has been previously investigated for psoriasis (PsO). This meta-analysis of genome-wide association studies (GWAS) was performed to identify genetic factors associated with increased risk/faster onset of herpes zoster (HZ) in subjects with RA or PsO receiving tofacitinib treatment, and to determine potential mechanisms that could be attributed to the varying rates of HZ across ethnicities. METHODS: In an ethnicity/indication-specific, trans-ethnic, trans-population meta-analysis of GWAS in subjects with RA or PsO from phase II, phase III, and long-term extension studies of tofacitinib, 8 million genetic variants were evaluated for their potential association with time to an HZ event and incidence of an HZ event (case versus control) with tofacitinib treatment, using Cox proportional hazard and logistic regression analyses, respectively. RESULTS: In total, 5,246 subjects were included (3,168 with RA and 2,078 with PsO). After adjustment for age, baseline absolute lymphocyte count, genetically defined ethnicity, and concomitant methotrexate use (in RA subjects only), 4 loci were significantly associated with faster onset of HZ in European subjects (P < 5 × 10-8 ), including a single-nucleotide polymorphism (SNP) near CD83 (frequency of risk allele ~2% in European subjects versus ~0.1% in East Asian subjects). In the trans-ethnic, trans-population meta-analysis, the CD83 SNP remained significant. Four additional significant loci were identified in the meta-analysis, among which a SNP near IL17RB was associated with faster onset of HZ (meta-analysis hazard ratio 3.6 [95% confidence interval 2.40-5.44], P = 7.6 × 10-10 ; frequency of risk allele ~12% in East Asian subjects versus <0.2% in European subjects). CONCLUSION: Genetic analysis of tofacitinib-treated subjects with RA or PsO identified multiple loci associated with increased HZ risk. Prevalent variants near the immune-relevant genes CD83 and IL17RB in European and East Asian populations, respectively, may contribute to risk of HZ in tofacitinib-treated subjects.

Decreased peripheral blood memory B cells are associated with the presence of interstitial lung disease in rheumatoid arthritis: a case-control study.

OBJECTIVES: Interstitial lung disease sometimes occurs in rheumatoid arthritis patients. Although the underlying immunological mechanisms responsible for interstitial lung disease associated with rheumatoid arthritis have not yet been clarified, some reports have suggested possible roles of B cells. To examine the role of B-cell subsets in interstitial lung disease in rheumatoid arthritis patients, we analyzed peripheral blood B-cell subsets. METHODS: We analyzed the frequencies of the peripheral blood B-cell subsets by flow cytometry in rheumatoid arthritis patients with and without interstitial lung disease (n = 16 and 81, respectively) and in healthy donors (n = 110) by high-resolution computed tomography. RESULTS: Compared with healthy donors, rheumatoid arthritis patients showed statistically higher frequencies of naive B cells and lower frequencies of memory B cells. Moreover, the frequencies of memory B cells were lower in rheumatoid arthritis patients with interstitial lung disease than in those without. Multivariate analysis showed that the frequency of memory B cells, particularly switched memory B cells, was significantly decreased in rheumatoid arthritis patients with interstitial lung disease, even after adjusting for prednisolone dose. CONCLUSIONS: We suspect memory B cells play important roles in interstitial lung disease associated with rheumatoid arthritis.

Parsing multiomics landscape of activated synovial fibroblasts highlights drug targets linked to genetic risk of rheumatoid arthritis.

OBJECTIVES: Synovial fibroblasts (SFs) are one of the major components of the inflamed synovium in rheumatoid arthritis (RA). We aimed to gain insight into the pathogenic mechanisms of SFs through elucidating the genetic contribution to molecular regulatory networks under inflammatory condition. METHODS: SFs from RA and osteoarthritis (OA) patients (n=30 each) were stimulated with eight different cytokines (interferon (IFN)-α, IFN-γ, tumour necrosis factor-α, interleukin (IL)-1ß, IL-6/sIL-6R, IL-17, transforming growth factor-ß1, IL-18) or a combination of all 8 (8-mix). Peripheral blood mononuclear cells were fractioned into five immune cell subsets (CD4+ T cells, CD8+ T cells, B cells, natural killer (NK) cells, monocytes). Integrative analyses including mRNA expression, histone modifications (H3K27ac, H3K4me1, H3K4me3), three-dimensional (3D) genome architecture and genetic variations of single nucleotide polymorphisms (SNPs) were performed. RESULTS: Unstimulated RASFs differed markedly from OASFs in the transcriptome and epigenome. Meanwhile, most of the responses to stimulations were shared between the diseases. Activated SFs expressed pathogenic genes, including CD40 whose induction by IFN-γ was significantly affected by an RA risk SNP (rs6074022). On chromatin remodelling in activated SFs, RA risk loci were enriched in clusters of enhancers (super-enhancers; SEs) induced by synergistic proinflammatory cytokines. An RA risk SNP (rs28411362), located in an SE under synergistically acting cytokines, formed 3D contact with the promoter of metal-regulatory transcription factor-1 (MTF1) gene, whose binding motif showed significant enrichment in stimulation specific-SEs. Consistently, inhibition of MTF1 suppressed cytokine and chemokine production from SFs and ameliorated mice model of arthritis. CONCLUSIONS: Our findings established the dynamic landscape of activated SFs and yielded potential therapeutic targets associated with genetic risk of RA.

Integrated bulk and single-cell RNA-sequencing identified disease-relevant monocytes and a gene network module underlying systemic sclerosis.

OBJECTIVE: Immunological disturbances have been reported in systemic sclerosis (SSc). This study assessed the transcriptome disturbances in immune cell subsets in SSc and characterized a disease-related gene network module and immune cell cluster at single cell resolution. METHODS: Twenty-one Japanese SSc patients were enrolled and compared with 13 age- and sex-matched healthy controls (HC). Nineteen peripheral blood immune cell subsets were sorted by flow cytometry and bulk RNA-seq analysis was performed for each. Differential expression and pathway analyses were conducted. Iterative weighted gene correlation network analysis (iWGCNA) of each subset revealed clustered co-expressed gene network modules. Random forest analysis prioritized a disease-related gene module. Single cell RNA-seq analysis of 878 monocytes was integrated with bulk RNA-seq analysis and with a public database for single cell RNA-seq analysis of SSc patients. RESULTS: Inflammatory pathway genes were differentially expressed in widespread immune cell subsets of SSc. An inflammatory gene module from CD16+ monocytes, which included KLF10, PLAUR, JUNB and JUND, showed the greatest discrimination between SSc and HC. One of the clusters of SSc monocytes identified by single-cell RNA-seq analysis characteristically expressed these inflammatory co-expressed genes and was similar to lung infiltrating FCN1hi monocytes expressing IL1B. CONCLUSIONS: Our integrated analysis of bulk and single cell RNA-seq analysis identified an inflammatory gene module and a cluster of monocytes that are relevant to SSc pathophysiology. They could serve as candidate novel therapeutic targets in SSc.