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.

Autoantibodies to nuclear valosin-containing protein-like protein: systemic sclerosis-specific antibodies revealed by in vitro human proteome.

OBJECTIVES: To identify and characterize undescribed systemic sclerosis (SSc)-specific autoantibodies targeting nucleolar antigens and to assess their clinical significance. METHODS: We conducted proteome-wide autoantibody screening (PWAS) against serum samples from SSc patients with nucleolar patterned anti-nuclear antibodies (NUC-ANAs) of specific antibodies (Abs) unknown, utilizing wet protein arrays fabricated from in vitro human proteome. Controls included SSc patients with already-known SSc-specific autoantibodies, patients with other connective tissue diseases and healthy subjects. The selection of nucleolar antigens was performed by database search in the Human Protein Atlas. The presence of autoantibodies was certified by immunoblots and immunoprecipitations. Indirect immunofluorescence assays on HEp-2 cells were also conducted. Clinical assessment was conducted by retrospective review of electronic medical records. RESULTS: PWAS identified three candidate autoantibodies, including anti-nuclear valosin-containing protein-like (NVL) Ab. Additional measurements in disease controls revealed that only anti-NVL Abs are exclusively detected in SSc. Detection of anti-NVL Abs was reproduced by conventional assays such as immunoblotting and immunoprecipitation. Indirect immunofluorescence assays demonstrated homogeneous nucleolar patterns. Anti-NVL Ab-positive cases were characterized by significantly low prevalence of diffuse skin sclerosis and interstitial lung disease, compared with SSc cases with NUC-ANAs other than anti-NVL Abs, such as anti-U3-RNP and anti-Th/To Abs. CONCLUSION: Anti-NVL Ab is an SSc-specific autoantibody associated with a unique combination of clinical features, including limited skin sclerosis and lack of lung involvement.

Lincoln Sign in Synovitis, Acne, Pustulosis, Hyperostosis, and Osteomyelitis Syndrome.

Synovitis, acne, pustulosis, hyperostosis, and osteitis (SAPHO) syndrome is characterized by skin lesions and osteomyelitis. Although the pathogenesis of SAPHO syndrome remains unclear, chronic focal infection may play a role in the development of excessive immune responses.

Association Between Proton Pump Inhibitors and the Risk of Intestinal Behçet Disease.

OBJECTIVE: This study aimed to investigate the association between proton pump inhibitor (PPI) use and the incidence of intestinal Behçet disease (BD) in patients with BD. METHODS: A retrospective cohort study was conducted at The University of Tokyo Hospital, including patients with BD diagnosed between April 2005 and November 2023. Cox models and Kaplan-Meier analyses were used to evaluate hazard ratios (HRs) and cumulative incidence of intestinal BD, respectively. Secondary analyses were performed to assess the duration and dose-response relationship of PPI use. RESULTS: Among 194 patients with BD, 25.3% developed intestinal BD during a mean follow-up of 12 years. PPI users had a significantly higher incidence of intestinal BD compared to nonusers (adjusted HR 2.48, 95% CI 1.38-4.47, P = 0.002), with a confirmed duration/dose-dependent relationship. The cumulative incidence of intestinal BD was markedly elevated in PPI users (log-rank P < 0.001). The result was similar to that in the propensity score-matched cohort. CONCLUSION: This study demonstrates a significant association between PPI use and increased incidence of intestinal BD in patients with BD. Caution in prescribing PPIs for patients with BD is warranted due to the potential risk of severe complications associated with intestinal BD.

Upper gastrointestinal involvement of Behçet's disease in Japan: endoscopic findings and clinical features.

AIM: Behçet's disease (BD) can involve any gastrointestinal (GI) tract site. We analyzed the characteristics, risk factors, and treatment responses to upper GI (UGI) involvement in patients with BD. METHODS: This retrospective cohort study analyzed UGI findings in 101 patients with BD who underwent endoscopy between April 2005 and December 2022 at the University of Tokyo Hospital. The patients were divided into two groups based on the presence or absence of UGI findings. Patient backgrounds, clinical symptoms, colonoscopy (CS) findings, and blood test findings were compared between the groups. RESULTS: In total, 18.8% (19/101) of the patients had UGI lesions. The prevalence rates in the esophagus, stomach, and duodenum were 6.9%, 6.9%, and 8.9%, respectively. Of these 19 patients, BD treatment were intensified in 10 (52.6%) patients after esophagogastroduodenoscopy (EGD), and all showed improvement in symptoms or endoscopic findings. In the multivariate analysis, symptoms (OR: 37.1, P < 0.001), CRP > 1 mg/dL (OR: 11.0, P = 0.01), and CS findings (OR: 5.16, P = 0.04) were independent predictors of UGI involvement in BD patients. The prediction model for UGI involvement using these three factors was highly accurate, with an AUC of 0.899 on the ROC curve. In the subgroup analysis of intestinal BD, symptoms (OR: 12.8, P = 0.01) and ESR > 20 mm/h (OR: 11.5, P = 0.007) were independent predictors. CONCLUSIONS: EGD should be conducted in BD patients with high CRP, GI symptoms, and lower GI involvement, which leads to better management of BD in terms of improving symptoms and endoscopic findings.

Baricitinib for anti-melanoma differentiation-associated protein 5 antibody-positive dermatomyositis-associated interstitial lung disease: a case series and literature review on Janus kinase inhibitors for the disease.

Anti-melanoma differentiation-associated protein 5 antibody-positive dermatomyositis (anti-MDA5-DM) is frequently complicated by progressive interstitial lung disease (ILD), the prognosis of which is poor, and management is a major challenge. We treated three patients with anti-MDA5-DM-associated ILD (anti-MDA5-DM-ILD) using the Janus kinase (JAK) inhibitor, baricitinib, which improved lung opacities and saved two patients. We reviewed 6 patients with anti-MDA5-DM-ILD who had been treated with tofacitinib at our institution. Five of the patients survived, although discontinuation of tofacitinib due to complications was frequently observed. In addition, a literature search of patients with anti-MDA5-DM-ILD who were treated with JAK inhibitors yielded 21 articles involving 79 cases. All patients except one were treated with tofacitinib, and the survival rate was 75.9%. Although not statistically confirmed, the deceased patients tended to be older and had higher ferritin levels. A total of 92 complications were observed, 11 of which resulted in JAK inhibitor discontinuation. Cytomegalovirus reactivation comprised a substantial percentage of all complications and of those patients who required JAK inhibitor discontinuation. Five cases with fatal infective complications were also observed. While tofacitinib has been proposed to be a therapeutic option for anti-MDA5-DM-ILD, other JAK inhibitors, including baricitinib, are a treatment option. Further investigation is warranted to optimize treatment of anti-MDA5-DM-ILD.

Characteristics of anti-melanoma differentiation associated gene 5 antibody-positive dermatomyositis with thrombotic microangiopathy.

OBJECTIVES: Anti-melanoma differentiation associated gene 5 antibody (anti-MDA5 Ab)-positive dermatomyositis (DM) is a representative of rapidly progressive interstitial pneumonia. However, its association with thrombotic microangiopathy (TMA), characterized by thrombocytopenia, haemolytic anaemia, and organ dysfunction, has not been defined. This study aimed to elucidate the characteristics of anti-MDA5 Ab-positive DM accompanied by TMA. METHODS: We reviewed our hospital records from November 2009 to September 2022. We included patients in accordance with the 2017 European League Against Rheumatism/American College of Rheumatology classification criteria and the criteria of Bohan and Peter. TMA was diagnosed according to the criteria for transplantation-associated TMA proposed by the International Working Group. RESULTS: This study enrolled a total of 26 anti-MDA5 Ab-positive DM patients, four of whom developed TMA. The patients with TMA had an increased urine protein/creatinine ratio. In addition, these four of them showed significantly elevated levels of ferritin and anti-MDA5 Ab titers and were considered to have high disease activity; yet, all of them survived. CONCLUSIONS: Our study indicated that anti-MDA5 Ab-positive DM patients with hyperferritinemia, a high anti-MDA5 Ab titer, and an increased urine protein/creatinine ratio should be carefully managed, bearing in mind a complication of TMA.

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.

Significance of anti-transcobalamin receptor antibodies in cutaneous arteritis revealed by proteome-wide autoantibody screening.

Cutaneous arteritis (CA) is a single-organ vasculitis that exclusively affects the small to medium-sized arteries of the skin. Diagnosis depends on a histological investigation with skin biopsy, which could be burdensome for both patients and clinicians. Moreover, the pathogenesis of CA remains unstudied, and treatment has not yet been established. Herein, we applied our proteome-wide autoantibody screening method to explore autoantibodies in the serum of CA patients. As a result, anti-transcobalamin receptor (TCblR) antibodies (Abs) were specifically detected in 24% of CA patients. Patients with positive anti-TCblR Abs were spared from peripheral neuropathy compared to those with negative anti-TCblR Abs, showing characteristics as CA confined to the skin. In addition, we revealed that anti-TCblR Abs trigger the autocrine loop of interleukin-6 mediated by tripartite motif-containing protein 21 in human endothelial cells and induce periarterial inflammation in murine skin. Furthermore, we demonstrated that methylcobalamin, a ligand of TCblR, ameliorates inflammation caused by anti-TCblR Abs both in vitro and in vivo. Collectively, our investigation unveils the pathologic significance of anti-TCblR Abs in CA and their potential as a diagnostic marker and a pathophysiology-oriented therapeutic target.

Long-term risks of malignancy in myositis-specific antibody-positive idiopathic inflammatory myopathy.

To identify the long-term risks of malignancy in patients with myositis-specific antibody (MSA)-positive idiopathic inflammatory myopathy (IIM). This retrospective cohort study included 216 IIM patients (aged > 18 years). Of these, 109 patients were positive for antibodies against anti-aminoacyl-tRNA synthetase (ARS), melanoma differentiation-associated gene 5 (MDA5), Mi-2, and transcriptional intermediary factor 1-γ (TIF1-γ). Age- and sex-matched standardized incidence ratios (SIRs) were calculated to compare the incidence of malignancy in IIM patients to that of the general population. The malignancy-free survival rate was estimated by Kaplan-Meier methods. Our study included 109 patients, 64 with anti-ARS, 28 with anti-MDA5, 9 with anti-Mi-2, and 8 with anti-TIF1-γ antibodies; 16 and 5 patients were diagnosed with a malignancy within 3 years before or after and within 4 to 10 years after their IIM onset, respectively. The SIRs of malignancy within 3 years of IIM onset for each MSA were calculated as follows: 2.12 (95% confidence interval [CI] 0.98-4.35) for anti-ARS, 1.87 (95% CI 0.48-4.97) for anti-MDA5, 2.11 (95% CI 0.11-13.69) for anti-Mi-2, and 9.30 (95% CI 2.98-25.58) for anti-TIF1-γ antibodies. The SIR at 4 to 10 years after IIM onset in patients with an anti-MDA5 antibody was 4.62 (95% CI 1.19-14.72); other MSAs did not have statistically significant SIRs. The long-term SIR of malignancy in patients with an anti-MDA5 antibody was 4.62 (95% CI 1.19-14.72), and the SIR among patients with an anti-TIF1-γ antibody within 3 years of IIM onset was 9.30 (95% CI 2.98-25.58). Screening for malignancies in patients with late phase of IIM and an anti-MDA5 antibody may be beneficial.

A diagnostic score for eosinophilic granulomatosis with polyangiitis among eosinophilic disorders.

BACKGROUND: Eosinophilic granulomatosis with polyangiitis (EGPA) is a form of systemic vasculitis with eosinophilic inflammation. However, existing classification criteria are all designed to classify EGPA among vasculitis and there is no established method distinguishing EGPA from other eosinophilic disorders. The aim of the present study was to propose a scoring system to differentiate EGPA among eosinophilic disorders. METHODS: Non-supervised hierarchical clustering using Ward's method and principal component analysis (PCA) were performed for 19 clinical parameters of 58 patients with eosinophilia-related diseases at a tertiary university hospital. The newly proposed scoring system was externally validated in 40 patients at another tertiary institution. RESULTS: Two distinct clusters were identified, and clinical features including peripheral neuropathy, asthma, skin involvement, lung involvement, rheumatoid factor (RF) positivity, myeloperoxidase (MPO)-anti-neutrophil cytoplasmic antibody (ANCA) positivity, IgE elevation, C-reactive protein (CRP) elevation, and vasculitis pathological findings were predominantly observed in one of these clusters (p < 0.05). Ten features defining the cluster with a high rate of vasculitis were weighted by PCA to create the E-CASE (EGPA classification among systemic eosinophilia) scoring system, on a 16-point scale. Based on the distribution of scores in the primary cohort, we defined an E-CASE score ≥12 as positive, ≤ 8 as negative, and 9-11 as undeterminable. The sensitivity and specificity of the E-CASE score in the validation cohort were 93.3% and 100%, respectively. CONCLUSIONS: We developed and verified a novel scoring system for differentiating EGPA from other types of eosinophilic disorders.

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.

Autotaxin is a potential link between genetic risk factors and immunological disturbances of plasmacytoid dendritic cells in systematic lupus erythematosus.

BACKGROUND: The importance of autotaxin, an enzyme that catalyzes lysophospholipid production, has recently been recognized in various diseases, including cancer and autoimmune diseases. Herein, we examined the role of autotaxin in systemic lupus erythematosus (SLE), utilizing data from ImmuNexUT, a comprehensive database consisting of transcriptome data and expression quantitative trait locus (eQTL) data of immune cells from patients with immune-mediated disorders. METHODS: Serum autotaxin concentrations in patients with SLE and healthy controls (HCs) were compared. The transcriptome data of patients with SLE and age- and sex-matched HCs were obtained from ImmuNexUT. The expression of ENPP2, the gene encoding autotaxin, was examined in peripheral blood immune cells. Next, weighted gene correlation network analysis (WGCNA) was performed to identify genes with expression patterns similar to ENPP2. The ImmuNexUT eQTL database and public epigenomic databases were used to infer the relationship between autotaxin and pathogenesis of SLE. RESULTS: Autotaxin levels were elevated in the serum of patients with SLE compared to HCs. Furthermore, the expression of ENPP2 was higher in plasmacytoid dendritic cells (pDCs) than in other immune cell subsets, and its expression was elevated in pDCs of patients with SLE compared to HCs. In WGCNA, ENPP2 belonged to a module that correlated with disease activity. This module was enriched in interferon-associated genes and included genes whose expression was influenced by single-nucleotide polymorphisms associated with SLE, suggesting that it is a key module connecting genetic risk factors of SLE with disease pathogenesis. Analysis utilizing the ImmuNexUT eQTL database and public epigenomic databases suggested that the increased expression of ENPP2 in pDCs from patients with SLE may be caused by increased expression of interferon-associated genes and increased binding of STAT3 complexes to the regulatory region of ENPP2. CONCLUSIONS: Autotaxin may play a critical role in connecting genetic risk factors of SLE to disease pathogenesis in pDCs.

ANCA-associated vasculitis with protein-losing enteropathy is characterized by hypocomplementemia.

Protein-losing enteropathy (PLE) has been reported to be associated with various systemic autoimmune diseases. However, reports regarding PLE in ANCA-associated vasculitis (AAV) patients are limited. We herein aimed to describe the clinical characteristics of AAV with PLE. We conducted a retrospective chart review of patients who were diagnosed with AAV and who began treatment at the University of Tokyo Hospital between June 2003 and June 2020. Among 68 AAV patients, there were four patients (5.9%) with PLE, consisting of two patients with MPA, one patient with GPA, and one patient with EGPA. Clinical courses were described, and their data were compared with AAV patients without PLE. Demographic characteristics, disease activity, and the pattern of organ involvement were similar between patients with PLE and without PLE. Patients with PLE had hypocomplementemia more frequently than the patients without PLE (CH50 75.0% vs 1.8%, p < 0.001, C3 50.0% vs 1.8%, p = 0.01, C4 75.0% vs 3.5%, p = 0.001). Although hypoalbuminemia improved with immunosuppressive therapy for AAV, the improvement in hypoalbuminemia was slow in most cases. We also performed a systematic review on PLE associated with vasculitis. Thirteen reports were included, and Henoch-Schonlein Purpura patients with PLE also tended to have hypocomplementemia. In conclusion, PLE is a rare complication of AAV and complement system may associate with the mechanism of PLE.

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.