Our analysis found that 34 upregulated genes from DMARDs-naïve RA patients were predicted to be activated by MyD88 (z-score = 5.384; p-value = 1.23E-16), which we referred to as MyD88-dependent genes (Figure 2C) (Interaction confidence = 0.9; PPI enrichment p-value < 1.1E-13). This evidence concerns the gene MYD88 and rheumatoid arthritis.