Integration of single-cell and bulk RNA sequencing revealed immune heterogeneity and its association with disease activity in rheumatoid arthritis patients.

Rheumatoid arthritis (RA) is a chronic, inflammatory, systemic autoimmune disease characterized by cartilage, bone damage, synovial inflammation, hyperplasia, autoantibody production, and systemic features. To obtain an overall profile of the immune environment in RA patients and its association with clinical features, we performed single-cell transcriptome and T-cell receptor sequencing of mononuclear cells from peripheral blood (PBMC) and synovial fluid (SF) from RA patients, integrated with two large cohorts with bulk RNA sequencing for further validation and investigation. Dendritic cells (DCs) exhibited relatively high functional heterogeneity and tissue specificity in relation to both antigen presentation and proinflammatory functions. Peripheral helper T cells (TPHs) are likely to originate from synovial tissue, undergo activation and exhaustion, and are subsequently released into the peripheral blood. Notably, among all immune cell types, TPHs were found to have the most intense associations with disease activity. In addition, CD8 effector T cells could be clustered into two groups with different cytokine expressions and play distinct roles in RA development. By integrating single-cell data with bulk sequencing from two large cohorts, we identified interactions among TPHs, CD8 cells, CD16 monocytes, and DCs that strongly contribute to the proinflammatory local environment in RA joints. Of note, the swollen 28-joint counts exhibited a more pronounced association with this immune environment compared to other disease activity indexes. The immune environment alternated significantly from PBMCs to SF, which indicated that a series of immune cells was involved in proinflammatory responses in the local joints of RA patients. By integrating single-cell data with two large cohorts, we have uncovered associations between specific immune cell populations and clinical features. This integration provides a rapid and precise methodology for assessing local immune activation, offering valuable insights into the pathophysiological mechanisms at play in RA.