9 HC-enriched pathways and 7 HC-enriched species were negatively correlated with ESR, CRP, RF, or anti-CCP, while 7 RA-enriched pathways were positively correlated with at least one of these clinical indices. The gene discussed is ESR1; the disease is rheumatoid arthritis.