We also used chromVAR21 to compute transcription factor (TF) motif deviation in single cells by estimating the enrichment of TF binding motifs in open chromatin regions and examined the enrichment of TF motifs in immune cell types concerning diagnosis and identified NFKB1 TF motifs with increased enrichment with BD patients in monocytes (cluster 18, Fig. 1e, Supplementary Fig. 2b, d). This evidence concerns the gene TF and Behcet disease.