Assessing the correlation with DAS28 for each gene individually, the most positively correlated gene was S100A8 with mean R = 0.28 (95% CI [0.19, 0.37]) and most anti-correlated gene HSP90AB1 with mean r = -0.23 (95% CI [-0.32, -0.14]) (Supplementary Figures 7, 8A). The RA score performed better than any single gene, positively correlated with DAS28 where the average correlation was 0.33 with 95% CI [0.24, 0.41] (Supplementary Figures 8B, C), suggesting this score could be helpful as a disease activity biomarker. Here, S100A8 is linked to rheumatoid arthritis.