GNLY and rheumatoid arthritis: The RA dataset contained 10,211 cells divided into 18 clusters (Figure 2A), which annotated as follows: T cells characterized by high expression of CD3D and CD3E (clusters 0,1,2,3,9,10,11 and12); monocytes with high expression of STXBP2 and FCN1 (cluster 4 and 7), NK cells with GNLY and NKG7 as marker genes (clusters 5and15); B cells marked by MS4A1 and CD79B (clusters 6,8and13); plasmablasts characterized by high expression of IGJ and CD27 (cluster 14); dendritic cells (DCs) (cluster 16) and pDCs (cluster 17) marked by CD1C, ENHO and PTPRA, MAP1A, respectively (Figure 2B).