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.

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.

Associations between vitamin D and autoimmune diseases: Mendelian randomization analysis.

OBJECTIVE: The VITAL trial of vitamin D supplementation suggested a possible protective effect for autoimmune diseases but uncertainties remain. We investigated potential causal effects of vitamin D on composite and individual autoimmune diseases using Mendelian randomization. METHODS: We used data from 332,984 participants of the UK Biobank of whom 23,089 had at least one autoimmune disease defined using ICD code and/or self-report. Diseases were further considered in mechanistic subgroups driven by "autoimmunity" (n = 12,774) or "autoinflammation" (n = 11,164), then individually. We selected variants within gene regions implicated in vitamin D biology to generate a weighted genetic score. We performed population-wide analysis using the ratio method, then examined non-linear effects across five quantiles based on 25-hydroxycholecalciferol levels. RESULTS: Genetically-predicted vitamin D was associated with lower risk of diseases in the autoinflammation group (OR 0.95 per 10 ng/ml increase in 25-hydroxycholecalciferol; 95%CI 0.91-0.99; p = 0.03) but not the autoimmunity group (OR 0.99; 95%CI 0.95-1.03; p = 0.64) or combined. When considering individual diseases, genetically-predicted vitamin D was associated with lower risk of psoriasis (OR 0.91; 95%CI 0.85-0.97; p = 0.005), the most common disease in the autoinflammation group, and suggestively with systemic lupus erythematosus (OR 0.84; 95%CI 0.69-1.02; p = 0.08); results were replicated using data from independent studies. We found no evidence for a plausible non-linear relationship between vitamin D and any outcome. CONCLUSIONS: We found genetic evidence to support a causal link between 25-hydroxycholecalciferol concentrations and psoriasis and systemic lupus erythematosus. These results have implications for potential disease prevention strategies, and the interpretation and design of vitamin D supplementation trials.

Lessons From Transcriptome Analysis of Autoimmune Diseases.

Various immune cell types, including monocytes, macrophages, and adaptive immune T and B cells, play major roles in inflammation in systemic autoimmune diseases. However, the precise contribution of these cells to autoimmunity remains elusive. Transcriptome analysis has added a new dimension to biology and medicine. It enables us to observe the dynamics of gene expression in different cell types in patients with diverse diseases as well as in healthy individuals, which cannot be achieved with genomic information alone. In this review, we summarize how transcriptome analysis has improved our understanding of the pathological roles of immune cells in autoimmune diseases with a focus on the ImmuNexUT database we reported. We will also discuss the common experimental and analytical design of transcriptome analyses. Recently, single-cell RNA-seq analysis has provided atlases of infiltrating immune cells, such as pro-inflammatory monocytes and macrophages, peripheral helper T cells, and age or autoimmune-associated B cells in various autoimmune disease lesions. With the integration of genomic data, expression quantitative trait locus (eQTL) analysis can help identify candidate causal genes and immune cells. Finally, we also mention how the information obtained from these analyses can be used practically to predict patient prognosis.

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.

Clinical Characteristics of Patients with Spondyloarthritis in Japan in Comparison with Other Regions of the World.

OBJECTIVE: To delineate clinical characteristics of patients with spondyloarthritis (SpA) in Japan in comparison to other areas of the world. METHODS: Using the ASAS-COMOSPA (Assessment of Spondyloarthritis international Society-COMOrbidities in SPondyloArthritis) data, an international cross-sectional observational study of patients with SpA, we analyzed information on demographics, disease characteristics, comorbidities, and risk factors. Patients were classified by region: Japan, other Asian countries (China, Singapore, South Korea, Taiwan), and non-Asian countries (Europe, the Americas, Africa). Patient characteristics, including diagnosis and treatment, were compared. RESULTS: Among 3984 patients included in the study, 161 were from centers in Japan, 933 from other Asian countries, and 2890 from other regions. Of patients with SpA in Japan, 42 (26.1%) had peripheral SpA, substantially more than in other countries. This trend was explained by the predominance of psoriatic arthritis (PsA) among Japanese patients with SpA. In contrast to the relatively low number in Japan, 54% of patients from other Asian countries had pure axial SpA (axSpA) without peripheral features. HLA-B27 testing, considered an integral part of the classification of axSpA, was performed in only 63.6% of Japanese patients with axSpA. More than half of Japanese patients with axSpA were classified using imaging criteria. CONCLUSION: In our study, there was a more substantial number of peripheral SpA cases observed in Japan compared to other parts of Asia and other regions of the world. Aside from ethnic differences, increasing recognition of PsA in Japan, as well as a potential underdiagnosis of axSpA due to the insufficient use of HLA-B27 testing, may partly explain regional discrepancies.