Unveiling the dynamics of B lymphocytes in systemic lupus erythematosus patients treated with belimumab through longitudinal single-cell RNA sequencing.

OBJECTIVES: Unraveling the mechanisms underlying treatment response for targeted therapeutics in systemic lupus erythematosus (SLE) patients is challenging due to the limited understanding of diverse responses of circulating immune cells, particularly B cells. We investigated B lymphocyte dynamics during anti-BAFF treatment, utilizing longitudinal single-cell transcriptome data. METHODS: We conducted single-cell RNA sequencing on PBMCs in four Korean SLE patients before and after belimumab treatment at the following time points: 2 weeks, 1, 3, 6, and 12 months. RESULTS: Analyzing over 73 000 PBMCs, we identified 8 distinct subsets of B cells and plasmablasts and analyzed dynamic changes within these cell subsets: initial declines in naive and transitional B cells followed by an increase at three months, contrasted by an initial increase and subsequent decrease in memory B cells by the third month. Meanwhile, plasmablasts exhibited a consistent decline throughout treatment. B cell activation pathways, specifically in naive and memory B cells, were downregulated during the third and sixth months. These findings were validated at the protein level throughout the first four weeks of treatment using flow cytometry. Comparative analysis with bulk transcriptome data from 22 Japanese SLE patients showed increased NR4A1 expression six months post-belimumab treatment, indicating its role in restricting self-reactive B cells, thereby contributing to the biological responses of anti-BAFF treatment. CONCLUSION: The observed B cell dynamics provided insights into the immunological mechanisms underlying the therapeutic effects of anti-BAFF in SLE patients. Furthermore, it underscores the need for research in predicting drug responses based on immune profiling.

Associations of upper respiratory mucosa microbiota with Rheumatoid arthritis, autoantibodies, and disease activity.

The lung is recognized as a site for initiating the formation of self-antigen and autoimmune responses in rheumatoid arthritis (RA). We aimed to investigate the association of upper respiratory microbiota with RA, autoantibody production, and disease activity. Forty-six patients with RA and 17 controls were examined. Nasopharyngeal swab samples were sequenced for microbiome profiling using the V3-V4 region of the 16S rRNA gene. The microbial diversity and relative abundance were compared between RA patients and controls. Correlation analyses were conducted to evaluate the relationship between microbial abundance and clinical markers such as autoantibodies and disease activity. Microbial diversity analysis revealed no major differences between RA patients and healthy controls. However, beta diversity analysis indicated a subtle distinction in microbial composition (unweighted UniFrac distance) between the two groups (P = 0.03), hinting at a minor subset of microbiota associated with disease status. Differential abundance analysis uncovered specific taxa at various taxonomic levels, including Saccharibacteria (TM7) [O-1] (PFDR = 2.53 × 10-2), TM7 [F-1] (PFDR = 5.20 × 10-3), Microbacterium (PFDR = 3.37 × 10-4), and Stenotrophomonas (PFDR = 2.57 × 10-3). The relative abundance of ten genera correlated significantly with anti-cyclic citrullinated peptide (anti-CCP) antibody levels (PFDR < 0.05) and 11 genera were significantly associated with disease activity markers, including ESR, CRP, DAS28-ESR, and DAS-CRP (PFDR < 0.05). In particular, Saccharibacteria TM7 [G-3] and Peptostreptococcaceae [XI] [G-1] were correlated with all disease activity biomarkers. Dysbiosis in the upper respiratory mucosa is associated with RA, anti-CCP antibody levels, and disease activity.

Association of HLA-DRB1 locus with treatment response to abatacept or TNF inhibitors in patients with seropositive rheumatoid arthritis.

The strongest genetic risk factor for rheumatoid arthritis (RA) has been known as HLA-DRB1 based on amino acid positions 11, 71, and 74. This study analyzed the association between specific HLA-DRB1 locus and treatment response to abatacept or TNF inhibitors (TNFi) in patients with seropositive RA. A total of 374 Korean RA patients were treated with abatacept (n = 110) or TNFi (n = 264). Associations between HLA-DRB1 and treatment response after 6 months were analyzed using multivariable logistic regression. Seropositive RA patients with HLA-DRB1 shared epitope (SE) had a favorable response to abatacept (OR = 3.67, P = 0.067) and an inversely associated response to TNFi (OR 0.57, P = 0.058) based on EULAR response criteria, but the difference was not statistically significant in comparison to those without SE. In analyses using amino acid positions of HLA-DRB1, a significant association was found between valine at amino acid position 11 of SE and good response to abatacept (OR = 6.46, P = 5.4 × 10-3). The VRA haplotype also showed a good response to abatacept (OR = 4.56, P = 0.013), but not to TNFi. Our results suggest that treatment response to abatacept or TNFi may differ depending on HLA-DRB1 locus in seropositive RA, providing valuable insights for selecting optimal therapy.