Most of the HC-enriched pathways and species were negatively correlated with ESR, CRP, RF, and anti-CCP, while the RA-enriched pathways and species were positively correlated with these clinical indices. This evidence concerns the gene ESR1 and rheumatoid arthritis.