Distinct transcriptome architectures underlying lupus establishment and exacerbation.

Systemic lupus erythematosus (SLE) is a complex autoimmune disease involving multiple immune cells. To elucidate SLE pathogenesis, it is essential to understand the dysregulated gene expression pattern linked to various clinical statuses with a high cellular resolution. Here, we conducted a large-scale transcriptome study with 6,386 RNA sequencing data covering 27 immune cell types from 136 SLE and 89 healthy donors. We profiled two distinct cell-type-specific transcriptomic signatures: disease-state and disease-activity signatures, reflecting disease establishment and exacerbation, respectively. We then identified candidate biological processes unique to each signature. This study suggested the clinical value of disease-activity signatures, which were associated with organ involvement and therapeutic responses. However, disease-activity signatures were less enriched around SLE risk variants than disease-state signatures, suggesting that current genetic studies may not well capture clinically vital biology. Together, we identified comprehensive gene signatures of SLE, which will provide essential foundations for future genomic and genetic studies.

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.

Machine learning-driven immunophenotypic stratification of mixed connective tissue disease corroborating the clinical heterogeneity.

OBJECTIVES: To stratify patients with mixed connective tissue disease (MCTD) based on their immunophenotype. METHODS: We analyzed the immunophenotype and transcriptome of 24 immune cell subsets from patients with MCTD, systemic lupus erythematosus (SLE), idiopathic inflammatory myopathy (IIM), and systemic sclerosis (SSc) from our functional genome database, ImmuNexUT (https://www.immunexut.org/). MCTD patients were stratified by employing machine learning models including Random Forest, trained by immunophenotyping data from SLE, IIM, and SSc patients. Transcriptomes were analyzed with gene set variation analysis (GSVA) and clinical features of MCTD subgroups were compared. RESULTS: This study included 215 patients, including 22 patients with MCTD. Machine learning models, constructed to classify SLE, IIM, and SSc patients based on immunophenotyping, were applied to MCTD patients, resulting in 16 classified as SLE-immunophenotype and 6 as non-SLE-immunophenotype. Among MCTD patients, patients with the SLE-immunophenotype had higher proportions of Th1 cells [2.85% (interquartile range (IQR) 1.54-3.91) vs 1.33% (IQR 0.99-1.74) p= 0.027] and plasmablasts [6.35% (IQR 4.17-17.49) vs 2.00% (IQR 1.20-2.80) p= 0.010]. Notably, the number of SLE-related symptoms was higher in patients with the SLE-immunophenotype [2.0 (IQR 1.0-2.0) vs 1.0 (IQR1.0-1.0) p= 0.038]. Moreover, GSVA scores of interferon-α and -γ responses were significantly higher in patients with the SLE-immunophenotype in central memory CD8+ T cells, while hedgehog signalling was higher in non-SLE-immunophenotype patients in 5 cell subsets. CONCLUSION: This study describes the stratification of MCTD patients based on immunophenotyping, suggesting the presence of distinct immunological processes behind the clinical subtypes of MCTD.

Multimodal repertoire analysis unveils B cell biology in immune-mediated diseases.

OBJECTIVES: Despite the involvement of B cells in the pathogenesis of immune-mediated diseases (IMDs), biological mechanisms underlying their function are scarcely understood. To overcome this gap, here we constructed and investigated a large-scale repertoire catalogue of five B cell subsets of patients with IMDs. METHODS: We mapped B cell receptor regions from RNA sequencing data of sorted B cell subsets. Our dataset consisted of 595 donors under IMDs and health. We characterised the repertoire features from various aspects, including their association with immune cell transcriptomes and clinical features and their response to belimumab treatment. RESULTS: Heavy-chain complementarity-determining region 3 (CDR-H3) length among naïve B cells was shortened among autoimmune diseases. Strong negative correlation between interferon signature strength and CDR-H3 length was observed in naïve B cells and suggested the role for interferon in premature B cell development. VDJ gene usage was skewed especially in plasmablasts and unswitched-memory B cells of patients with systemic lupus erythematosus (SLE). We developed a scoring system to quantify this skewing, and it positively correlated with peripheral helper T cell transcriptomic signatures and negatively correlated with the amount of somatic hyper mutations in plasmablasts, suggesting the association of extrafollicular pathway. Further, this skewing led to high usage of IGHV4-34 gene with 9G4 idiotypes in unswitched-memory B cells, which showed a prominent positive correlation with disease activity in SLE. Gene usage skewing in unswitched-memory B cells was ameliorated after belimumab treatment. CONCLUSIONS: Our multimodal repertoire analysis enabled us the system-level understanding of B cell abnormality in diseases.

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.

Functional evaluation of rare OASL variants by analysis of SLE patient-derived iPSCs.

BACKGROUND: Systemic lupus erythematosus (SLE) is a chronic systemic autoimmune disease characterized by genetic heterogeneity and an interferon (IFN) signature. The overall landscapes of the heritability of SLE remains unclear. OBJECTIVES: To identify and elucidate the biological functions of rare variants underlying SLE, we conducted analyses of patient-derived induced pluripotent stem cells (iPSCs) in combination with genetic analysis. METHODS: Two familial SLE patient- and two healthy donor (HD)-derived iPSCs were established. Type 1 IFN-secreting dendritic cells (DCs) were differentiated from iPSCs. Genetic analyses of SLE-iPSCs, and 117 SLE patients and 107 HDs in the ImmuNexUT database were performed independently. Genome editing of the variants on iPSCs was performed with the CRISPR/Cas9 system. RESULTS: Type 1 IFN secretion was significantly increased in DCs differentiated from SLE-iPSCs compared to HD-iPSCs. Genetic analyses revealed a rare variant in the 2'-5'-Oligoadenylate Synthetase Like (OASL) shared between SLE-iPSCs and another independent SLE patient, and significant accumulation of OASL variants among SLE patients (HD 0.93%, SLE 6.84%, OR 8.387) in the database. Genome editing of mutated OASL 202Q to wild-type 202 R or wild-type OASL 202 R to mutated 202Q resulted in reduced or enhanced Type 1 IFN secretion of DCs. Three other OASL variants (R60W, T261S and A447V) accumulated in SLE patients had also capacities to enhance Type 1 IFN secretion in response to dsRNA. CONCLUSIONS: We established a patient-derived iPSC-based strategy to investigate the linkage of genotype and phenotype in autoimmune diseases. Detailed case-based investigations using patient-derived iPSCs provide information to unveil the heritability of the pathogenesis of autoimmune diseases.

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.

Control of naive and effector CD4 T cell receptor repertoires by rheumatoid-arthritis-risk HLA alleles.

OBJECTIVE: Human Leukocyte Antigen (HLA) alleles regulate susceptibility to rheumatoid arthritis (RA) and immune-mediated diseases. This study aims to elucidate the impact of HLA alleles to T cell subsets. METHODS: We performed genome-wide and HLA allele association analysis for T cell receptor (TCR) beta chain repertoire in 13 purified T cell subsets from the ImmuNexUT database, consisting of 407 donors with ten immune-mediated diseases and healthy controls. RESULTS: HLA class II alleles were associated with TRBV gene usage and the public clones of CD4 T cells, while HLA class I alleles were associated with CD8 T cells. RA-risk and immune-mediated diseases-risk HLA alleles were associated with TRBV gene usage of naive and effector CD4 T cell subsets and public clones accumulating in Th17. Clonal diversity was independent of HLA alleles and was correlated with transcriptome changes that reflect TCR signaling. CONCLUSION: This study revealed in vivo evidence that both HLA alleles and environmental factors shape naive and effector TCR repertoires in RA and immune-mediated diseases patients.

Dysregulation of the gene signature of effector regulatory T cells in the early phase of systemic sclerosis.

OBJECTIVES: We evaluated flow-cytometric and transcriptome features of peripheral blood immune cells from early-phase (disease duration <5 years) SSc in comparison with late-phase SSc. METHODS: Fifty Japanese patients with SSc (12 early SSc cases and 38 late SSc cases) and 50 age- and sex-matched healthy controls were enrolled. A comparison of flow-cytometric subset proportions and RNA-sequencing of 24 peripheral blood immune cell subsets was performed. We evaluated differentially expressed genes (DEGs), characterized the co-expressed gene modules, and estimated the composition of subpopulations by deconvolution based on single-cell RNA-sequencing data. As a disease control, idiopathic inflammatory myositis (IIM) patients were also evaluated. RESULTS: Analysing the data from early and late SSc, fraction II effector regulatory T cell (Fr. II eTreg) genes showed a remarkable differential gene expression, enriched for genes related to oxidative phosphorylation. Although the flow-cytometric proportion of Fr. II eTregs was not changed in early SSc, deconvolution indicated expansion of the activated subpopulation. Co-expressed gene modules of Fr. II eTregs demonstrated enrichment of the DEGs of early SSc and correlation with the proportion of the activated subpopulation. These results suggested that DEGs in Fr. II eTregs from patients with early SSc were closely associated with the increased proportion of the activated subpopulation. Similar dysregulation of Fr. II eTregs was also observed in data from patients with early IIM. CONCLUSIONS: RNA-seq of immune cells indicated the dysregulation of Fr. II eTregs in early SSc with increased proportion of the activated subpopulation.

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.

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.

Immature platelet levels correlate with disease activity and predict treatment response of thrombocytopenia in lupus patients.

OBJECTIVE: The immature platelet fraction (IPF) represents recently produced platelets in bone marrow and this parameter is increased in patient with primary immune thrombocytopenia (ITP). We investigated the associations between IPF and absolute immature platelet count (AIPC), and clinical parameters in systemic lupus erythematosus (SLE), which has more complex pathological mechanisms than in primary ITP. METHODS: Patients with SLE were retrospectively reviewed at the University of Tokyo Hospital from May, 2012 to January, 2021. The correlations between clinical parameters and the number of immature platelets were assessed with Spearman's rank correlation coefficients. A multiple logistic regression model was used to identify the independent clinical parameters for IPF and AIPC. The difference in the distribution of time for a complete response (CR) after prednisolone (PSL) administration was also evaluated by log-rank test. RESULTS: A total of 282 SLE patients were enrolled, and 12.41% of those patients showed thrombocytopenia. IPF correlated with clinical parameters such as platelet count (r = -0.58), AIPC (r = 0.64) and systemic lupus erythematosus disease activity index 2000 (SLEDAI-2K) (r = 0.24). SLEDAI-2K [odds ratio (OR) (per unit increase), 1.07; 95% CI, 1.013 - 1.13] and thrombocytopenia (OR, 32.23; 95% CI, 11.072 - 93.80) were independent clinical parameters to account for IPF increase. IPF correlated with the number of bone marrow megakaryocytes (n = 19, r = 0.57). Notably, the probability of CR in response to PSL in AIPC-high patients was higher than in AIPC-low patients (hazard ratio, 4.62; 95% CI, 1.07 - 20.02). CONCLUSION: IPF correlated with disease activity of SLE and represented platelet production in the bone marrow, whereas AIPC predicted a rapid response to steroids in thrombocytopenic patients with SLE.

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.

A gene module associated with dysregulated TCR signaling pathways in CD4+ T cell subsets in rheumatoid arthritis.

We analyzed the transcriptome of detailed CD4+ T cell subsets including them after abatacept treatment, and examined the difference among CD4+ T cell subsets and identified gene sets that are closely associated disease activity and abatacept treatment. Seven CD4+ T cell subsets (naive, Th1, Th17, Th1/17, nonTh1/17, Tfh and Treg) were sorted from PBMCs taken from 10 RA patients and 10 healthy controls, and three RA patients donated samples before and 6 months after abatacept treatment. Paired-end RNA sequencing was performed using HiSeq 2500. A total of 149 samples except for 12 outliers were analyzed. Overview of expression pattern of RA revealed that administration of abatacept exerts a large shift toward the expression pattern of HC. Most of differentially expressed gene (DEG) upregulated in RA (n = 1776) were downregulated with abatacept treatment (n = 1349). Inversely, most of DEG downregulated in RA (n = 1860) were upregulated with abatacept treatment (n = 1294). This DEG-based analysis revealed shared pathway changes in RA CD4+ T cell subsets. Knowledge-based pathway analysis revealed the upregulation of activation-related pathways in RA that was substantially ameliorated by abatacept. Weighted gene co-expression network analysis (WGCNA) evaluated CD4+ T cells collectively and identified a gene module that consisted of 227 genes and was correlated with DAS28-CRP (Spearman's rho = 0.46, p = 4 × 10-9) and abatacept administration (Spearman's rho = -0.91, p = 5 × 10-57). The most highly connected 30 genes of this module included ZAP70 and JAK3, and pathway analysis of this module revealed dysregulation of the TCR signaling pathway network, which was ameliorated by abatacept.

Polymorphic lymphoproliferative disorders in patients with rheumatoid arthritis are associated with a better clinical outcome.

OBJECTIVES: The characteristics of lymphoproliferative disorders (LPD) in patients with rheumatoid arthritis (RA) remain unclear. Therefore, we retrospectively analyzed the clinical characteristics of these patients in our department. METHODS: Twenty RA patients who developed LPD between April 2003 and August 2016 in our department were analyzed. RESULTS: All of the RA patients who developed LPD had been treated with methotrexate (MTX). The median weekly and total dosages of MTX were 6.8 mg/week and 2530 mg, respectively. The median duration of MTX administration was eight years. Nineteen patients (95%) achieved complete remission (CR) and 15 (75%) achieved CR with MTX cessation alone. Based on the pathological findings, we divided MTX-associated LPD patients into two groups (n = 16); polymorphic LPD (31%) and other groups. CR with MTX cessation alone was achieved in 5 (100%) and 6 (54.5%) patients in the polymorphic LPD and other groups, respectively (p = .12). Moreover, the duration from the cessation of MTX to CR was significantly shorter in the polymorphic LPD group than in the other group (5.3 months vs 12.6 months, p = .01, respectively). CONCLUSION: Polymorphic LPD, which was the most frequent pathological diagnosis in this cohort, was associated with a higher incidence of CR and a significantly shorter duration to CR.

Enhanced gut homing receptor expression of unswitched memory B cells in rheumatoid arthritis.

****************************************************************************.

A case of refractory polyarteritis nodosa successfully treated with rituximab.

A 59-year-old man who presented with continuous fever, livedo reticularis, and left leg ischemia with multiple tibial artery stenosis and renal artery aneurysm, as demonstrated by arteriography, was diagnosed with polyarteritis nodosa (PAN) 6 years ago. Although he frequently relapsed in spite of intensive immunosuppressive therapies, the disease activity of PAN was controlled with repeated rituximab (RTX) therapies and steroid doses were tapered safely. Peripheral CD19+ B-cells disappeared soon after the 1st administration of RTX. Although CD19+ B-cells remained absent, 3.1% of CD3+CD20+ T-cells were observed in the peripheral blood prior to the 2nd administration of RTX. Recent studies have suggested the pathogenic role of CD3+CD20+ T-cells in autoimmune diseases in the context of RTX therapy; therefore, their roles in the pathogenesis of PAN also need to be considered.

Efficacy of intensive immunosuppression in exacerbated rheumatoid arthritis-associated interstitial lung disease.

OBJECTIVES: Acute or subacute exacerbations are recognized as a severe complication of rheumatoid arthritis-associated interstitial lung disease (RA-ILD). Nevertheless, the role of intensive immunosuppression in RA-ILD remains elusive. We attempted to evaluate the clinical characteristics and efficacy of immunosuppressive treatment in exacerbated RA-ILD. METHODS: Clinical data, including respiratory function, imaging, treatment, and prognosis, were retrospectively collected for 17 patients with RA-ILD who required hospitalization at the University of Tokyo Hospital due to an acute exacerbation (12 patients) or subacute exacerbation (5 patients). RESULTS: Patients with RA-ILD demonstrated a significantly higher titers of anticyclic citrullinated peptide antibodies compared with RA patients in Japanese Ninja registry, suggesting the role of adaptive immunity. Immunosuppressive treatment suppressed the deterioration of pulmonary functions with improved ground grass opacity and consolidation. In particular, in patients with less fibrosis on computed tomography (CT) images showed a better response to treatment. Although five patients treated with combination therapy, including cyclophosphamide, showed a severely decreased lung volume, these intensive therapies provided a good prognosis without fatalities for the average observation period of 474 days. CONCLUSIONS: Immunosuppressive therapy is effective for exacerbations of RA-ILD. For severe cases with low respiratory function, intensive therapy, including cyclophosphamide, has a potential to improve the prognosis.