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.

Chromosomal alterations among age-related haematopoietic clones in Japan.

The extent to which the biology of oncogenesis and ageing are shaped by factors that distinguish human populations is unknown. Haematopoietic clones with acquired mutations become common with advancing age and can lead to blood cancers1-10. Here we describe shared and population-specific patterns of genomic mutations and clonal selection in haematopoietic cells on the basis of 33,250 autosomal mosaic chromosomal alterations that we detected in 179,417 Japanese participants in the BioBank Japan cohort and compared with analogous data from the UK Biobank. In this long-lived Japanese population, mosaic chromosomal alterations were detected in more than 35.0% (s.e.m., 1.4%) of individuals older than 90 years, which suggests that such clones trend towards inevitability with advancing age. Japanese and European individuals exhibited key differences in the genomic locations of mutations in their respective haematopoietic clones; these differences predicted the relative rates of chronic lymphocytic leukaemia (which is more common among European individuals) and T cell leukaemia (which is more common among Japanese individuals) in these populations. Three different mutational precursors of chronic lymphocytic leukaemia (including trisomy 12, loss of chromosomes 13q and 13q, and copy-neutral loss of heterozygosity) were between two and six times less common among Japanese individuals, which suggests that the Japanese and European populations differ in selective pressures on clones long before the development of clinically apparent chronic lymphocytic leukaemia. Japanese and British populations also exhibited very different rates of clones that arose from B and T cell lineages, which predicted the relative rates of B and T cell cancers in these populations. We identified six previously undescribed loci at which inherited variants predispose to mosaic chromosomal alterations that duplicate or remove the inherited risk alleles, including large-effect rare variants at NBN, MRE11 and CTU2 (odds ratio, 28-91). We suggest that selective pressures on clones are modulated by factors that are specific to human populations. Further genomic characterization of clonal selection and cancer in populations from around the world is therefore warranted.

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.

Six-month safety and effectiveness of tofacitinib in patients with rheumatoid arthritis in Japan: Interim analysis of post-marketing surveillance.

OBJECTIVES: We evaluated the real-world safety/effectiveness of tofacitinib, an oral Janus kinase inhibitor for the treatment of rheumatoid arthritis (RA), in patients with RA in Japan registered in a post-marketing surveillance study. METHODS: This interim analysis included data from July 2013 to December 2018. Adverse events (AEs), serious AEs (SAEs), Simplified Disease Activity Index (SDAI)/Clinical Disease Activity Index (CDAI)/Disease Activity Score in 28 joints, erythrocyte sedimentation rate [DAS28-4(ESR)] scores, and rates of SDAI/CDAI/DAS28-4(ESR)-defined remission and low disease activity were analysed using 6 months of data. Risk factors for serious infections were assessed by multivariable analyses. RESULTS: Safety and disease activity were evaluated in 6866 and 6649 patients, respectively. Overall, 32.73%/7.37% of patients reported AEs/SAEs. Clinically important AEs with tofacitinib included serious infections/infestations [3.13% of patients; incidence rate (IR; patients with events) 6.91/100 patient-years (PY)], herpes zoster (3.63%; IR 8.02/100 PY), and malignancies (0.68%; IR 1.45/100 PY). SDAI/CDAI/DAS28-4(ESR) scores and remission/low disease activity rates improved over 6 months. Male sex, older age, Steinbrocker's stage IV, history of infection, and diabetes mellitus at baseline were independent risk factors for serious infection. CONCLUSIONS: In patients with RA receiving tofacitinib in Japan, safety was consistent with the reported profile, and disease activity improved over 6 months. STUDY IDENTIFIER: NCT01932372.

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.

Combined plasma metabolomic and transcriptomic analysis identify histidine as a biomarker and potential contributor in SLE pathogenesis.

OBJECTIVES: To investigate metabolite alterations in the plasma of SLE patients to identify novel biomarkers and provide insight into SLE pathogenesis. METHODS: Patients with SLE (n = 41, discovery cohort and n = 37, replication cohort), healthy controls (n = 30 and n = 29) and patients with RA (n = 19, disease control) were recruited. Metabolic profiles of the plasma samples were analysed using liquid chromatography-time-of-flight mass spectrometry and capillary electrophoresis-time-of-flight mass spectrometry. Transcriptome data was analysed using RNA-sequencing for 18 immune cell subsets. The importance of histidine (His) in plasmablast differentiation was investigated by using mouse splenic B cells. RESULTS: We demonstrate that a specific amino acid combination including His can effectively distinguish between SLE patients and healthy controls. Random forest and partial least squares-discriminant analysis identified His as an effective classifier for SLE patients. A decrease in His plasma levels correlated with damage accrual independent of prednisolone dosage and type I IFN signature. The oxidative phosphorylation signature in plasmablasts negatively correlated with His levels. We also showed that plasmablast differentiation induced by innate immune signals was dependent on His. CONCLUSIONS: Plasma His levels are a potential biomarker for SLE patients and are associated with damage accrual. Our data suggest the importance of His as a pathogenic metabolite in SLE pathogenesis.

[Myelodysplastic syndrome with der (1;7)(q10;p10) complicated with eosinophilia and organizing pneumonia].

The unbalanced translocation der (1;7)(q10;p10) is a characteristic cytogenetic abnormality observed in myelodysplastic syndrome (MDS). A 63-year-old man presented to our hospital with fever and lung disease. The chromosomal analysis of bone marrow cells showed 46, XY, +1, der (1;7)(q10;p10) in all four metaphases. The patient was diagnosed with MDS. Bronchoscope examination revealed organizing pneumonia. The patient's eosinophil count rose to 39% after 30 days. His fever and dyspnea worsened, and a skin rash (systemic erythema) appeared simultaneously. Therefore, the patient was commenced on azacitidine and corticosteroids. Although treatment with both drugs could control disease progression transiently, the WT-1 value and the percentage of myeloblasts in the patient's bone marrow increased. Therefore, the patient received hematopoietic stem cell transplantation from his haplo-identical donor daughter. Some reports have demonstrated that patients with MDS with der (1;7)(q10;p10) have better prognosis than those with other abnormalities, such as -7/7q-. However, reported cases with severe complications show very poor prognosis. MDS with der (1;7)(q10;p10) complicated by eosinophilia and organizing pneumonia have not been reported, and its prognosis is expected to be very poor. Our case suggests that such cases might quickly require hematopoietic stem cell transplantation before the disease worsens.

IMPACT: Genomic Annotation of Cell-State-Specific Regulatory Elements Inferred from the Epigenome of Bound Transcription Factors.

Despite significant progress in annotating the genome with experimental methods, much of the regulatory noncoding genome remains poorly defined. Here we assert that regulatory elements may be characterized by leveraging local epigenomic signatures where specific transcription factors (TFs) are bound. To link these two features, we introduce IMPACT, a genome annotation strategy that identifies regulatory elements defined by cell-state-specific TF binding profiles, learned from 515 chromatin and sequence annotations. We validate IMPACT using multiple compelling applications. First, IMPACT distinguishes between bound and unbound TF motif sites with high accuracy (average AUPRC 0.81, SE 0.07; across 8 tested TFs) and outperforms state-of-the-art TF binding prediction methods, MocapG, MocapS, and Virtual ChIP-seq. Second, in eight tested cell types, RNA polymerase II IMPACT annotations capture more cis-eQTL variation than sequence-based annotations, such as promoters and TSS windows (25% average increase in enrichment). Third, integration with rheumatoid arthritis (RA) summary statistics from European (N = 38,242) and East Asian (N = 22,515) populations revealed that the top 5% of CD4+ Treg IMPACT regulatory elements capture 85.7% of RA h2, the most comprehensive explanation for RA h2 to date. In comparison, the average RA h2 captured by compared CD4+ T histone marks is 42.3% and by CD4+ T specifically expressed gene sets is 36.4%. Lastly, we find that IMPACT may be used in many different cell types to identify complex trait associated regulatory elements.

Pannexin1 channel-dependent secretome from apoptotic tumor cells shapes immune-escape microenvironment.

Apoptotic cell death is a critical step in organism development and tissue homeostasis. Apoptotic cells affect immune cell activities in normal tissues. It is not clear whether similar cell death machinery causes tumor environments to evade anti-tumor immune responses. Here, using a mouse transplant model, we found a large number of tumor cells undergoing intrinsic apoptosis in tumors derived from the 4T1 breast cancer cell line, where neutrophils significantly accumulated. Interestingly, these apoptotic 4T1 tumor cells directly induced neutrophil extracellular traps (NETs) in a pannexin 1 (Panx1) channel-dependent manner, and knockdown of Panx1 in 4T1 cells led to a reduction in tumor size. Spermidine released through Panx1 from apoptotic 4T1 cells induced NETs in bone marrow-derived neutrophils in vitro. In addition, inhibition of spermidine synthesis suppressed tumor growth in the mouse transplant model. Collectively, our data suggested a new immune-escape mechanism for tumors by Panx1-mediated secretome from intrinsic apoptotic cells, which may provide a new therapeutic target for cancer.

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.

Multi-trait and cross-population genome-wide association studies across autoimmune and allergic diseases identify shared and distinct genetic component.

OBJECTIVES: Autoimmune and allergic diseases are outcomes of the dysregulation of the immune system. Our study aimed to elucidate differences or shared components in genetic backgrounds between autoimmune and allergic diseases. METHODS: We estimated genetic correlation and performed multi-trait and cross-population genome-wide association study (GWAS) meta-analysis of six immune-related diseases: rheumatoid arthritis, Graves' disease, type 1 diabetes for autoimmune diseases and asthma, atopic dermatitis and pollinosis for allergic diseases. By integrating large-scale biobank resources (Biobank Japan and UK biobank), our study included 105 721 cases and 433 663 controls. Newly identified variants were evaluated in 21 778 cases and 712 767 controls for two additional autoimmune diseases: psoriasis and systemic lupus erythematosus. We performed enrichment analyses of cell types and biological pathways to highlight shared and distinct perspectives. RESULTS: Autoimmune and allergic diseases were not only mutually classified based on genetic backgrounds but also they had multiple positive genetic correlations beyond the classifications. Multi-trait GWAS meta-analysis newly identified six allergic disease-associated loci. We identified four loci shared between the six autoimmune and allergic diseases (rs10803431 at PRDM2, OR=1.07, p=2.3×10-8, rs2053062 at G3BP1, OR=0.90, p=2.9×10-8, rs2210366 at HBS1L, OR=1.07, p=2.5×10-8 in Japanese and rs4529910 at POU2AF1, OR=0.96, p=1.9×10-10 across ancestries). Associations of rs10803431 and rs4529910 were confirmed at the two additional autoimmune diseases. Enrichment analysis demonstrated link to T cells, natural killer cells and various cytokine signals, including innate immune pathways. CONCLUSION: Our multi-trait and cross-population study should elucidate complex pathogenesis shared components across autoimmune and allergic diseases.

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.

Functional genetics for studying the human immune system.

Although small numbers of immune-mediated diseases are inherited due to rare genetic mutations, most are multifactorial diseases caused by multiple elements including genetic and environmental factors. In the case of autoimmune diseases, many disease-susceptibility genes, including several in the major histocompatibility gene complex, have been reported, and over the past 10 years, genome-wide association studies (GWAS) have been used to analyze disease-susceptibility loci in representative diseases. Furthermore, many disease-susceptibility variants have been found to be related to gene expression levels. The expression of genes involved in disease pathogenesis is often cell-type-specific, and this is closely related to epigenome alterations. Genomic information is present even before the onset of a disease and has a clear causal relationship to the disease (i.e. the outcome). Therefore, it is important to establish functional genetics in human immunology to understand the pathogenesis of diseases using these pieces of information. We can then apply these results to drug discovery. Here, we will review these issues, especially focusing on autoimmune diseases, and discuss current and future directions of human immune system research.

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.

Identification of U11snRNA as an endogenous agonist of TLR7-mediated immune pathogenesis.

The activation of innate immune receptors by pathogen-associated molecular patterns (PAMPs) is central to host defense against infections. On the other hand, these receptors are also activated by immunogenic damage-associated molecular patterns (DAMPs), typically released from dying cells, and the activation can evoke chronic inflammatory or autoimmune disorders. One of the best known receptors involved in the immune pathogenesis is Toll-like receptor 7 (TLR7), which recognizes RNA with single-stranded structure. However, the causative DAMP RNA(s) in the pathogenesis has yet to be identified. Here, we first developed a chemical compound, termed KN69, that suppresses autoimmunity in several established mouse models. A subsequent search for KN69-binding partners led to the identification of U11 small nuclear RNA (U11snRNA) as a candidate DAMP RNA involved in TLR7-induced autoimmunity. We then showed that U11snRNA robustly activated the TLR7 pathway in vitro and induced arthritis disease in vivo. We also found a correlation between high serum level of U11snRNA and autoimmune diseases in human subjects and established mouse models. Finally, by revealing the structural basis for U11snRNA's ability to activate TLR7, we developed more potent TLR7 agonists and TLR7 antagonists, which may offer new therapeutic approaches for autoimmunity or other immune-driven diseases. Thus, our study has revealed a hitherto unknown immune function of U11snRNA, providing insight into TLR7-mediated autoimmunity and its potential for further therapeutic applications.

Amino acid signatures of HLA Class-I and II molecules are strongly associated with SLE susceptibility and autoantibody production in Eastern Asians.

Human leukocyte antigen (HLA) is a key genetic factor conferring risk of systemic lupus erythematosus (SLE), but precise independent localization of HLA effects is extremely challenging. As a result, the contribution of specific HLA alleles and amino-acid residues to the overall risk of SLE and to risk of specific autoantibodies are far from completely understood. Here, we dissected (a) overall SLE association signals across HLA, (b) HLA-peptide interaction, and (c) residue-autoantibody association. Classical alleles, SNPs, and amino-acid residues of eight HLA genes were imputed across 4,915 SLE cases and 13,513 controls from Eastern Asia. We performed association followed by conditional analysis across HLA, assessing both overall SLE risk and risk of autoantibody production. DR15 alleles HLA-DRB1*15:01 (P = 1.4x10-27, odds ratio (OR) = 1.57) and HLA-DQB1*06:02 (P = 7.4x10-23, OR = 1.55) formed the most significant haplotype (OR = 2.33). Conditioned protein-residue signals were stronger than allele signals and mapped predominantly to HLA-DRB1 residue 13 (P = 2.2x10-75) and its proxy position 11 (P = 1.1x10-67), followed by HLA-DRB1-37 (P = 4.5x10-24). After conditioning on HLA-DRB1, novel associations at HLA-A-70 (P = 1.4x10-8), HLA-DPB1-35 (P = 9.0x10-16), HLA-DQB1-37 (P = 2.7x10-14), and HLA-B-9 (P = 6.5x10-15) emerged. Together, these seven residues increased the proportion of explained heritability due to HLA to 2.6%. Risk residues for both overall disease and hallmark autoantibodies (i.e., nRNP: DRB1-11, P = 2.0x10-14; DRB1-13, P = 2.9x10-13; DRB1-30, P = 3.9x10-14) localized to the peptide-binding groove of HLA-DRB1. Enrichment for specific amino-acid characteristics in the peptide-binding groove correlated with overall SLE risk and with autoantibody presence. Risk residues were in primarily negatively charged side-chains, in contrast with rheumatoid arthritis. We identified novel SLE signals in HLA Class I loci (HLA-A, HLA-B), and localized primary Class II signals to five residues in HLA-DRB1, HLA-DPB1, and HLA-DQB1. These findings provide insights about the mechanisms by which the risk residues interact with each other to produce autoantibodies and are involved in SLE pathophysiology.