Integrated multi-omics analysis indicates MARCO and ZBTB20 as candidate biomarkers in rheumatoid arthritis and tuberculosis.

BACKGROUND

Rheumatoid arthritis (RA) patients have a clearly increased possibility to develop tuberculosis (TB). However, the common molecular background, useful diagnostic molecules, and treatment-related targets of this comorbidity are still not sufficiently explained. We attempted to search the shared genetic factors and candidate drug targets for RA and TB by an integrated multi-omics way.

METHODS

Bulk transcriptome data sets for RA and TB were obtained from Gene Expression Omnibus (GEO). Differential expression analysis and weighted gene co-expression network analysis (WGCNA) were first performed to collect overlapped disease-related genes. Least Absolute Shrinkage and Selection Operator (LASSO) regression was then applied to select hub genes. External data sets, single-cell RNA sequencing (scRNA-seq), and quantitative real-time PCR (qRT-PCR) using clinical samples were used for validation. We also carried out immune infiltration analysis, gene set enrichment analyses, molecular docking, and 100 ns molecular dynamics (MD) simulation, in order to describe possible function and treatment value of these genes.

RESULTS

A total of 312 co-expressed genes common to RA and TB were obtained, and these genes were mainly gathered in immune and inflammatory processes. MARCO and ZBTB20 were kept as hub genes and showed acceptable diagnostic ability in RA (AUC = 0.839 and 0.795, respectively) and TB (AUC = 0.849 and 0.915, respectively). Single-cell results suggested their expression in T cells, B cells, and myeloid cell lineages, tightly correlated with immune infiltration and inflammatory signaling. Computational drug-repurposing analysis identified Tipifarnib as a candidate compound, and molecular docking and molecular dynamics simulations predicted potentially stable interactions with MARCO and ZBTB20. In the clinical cohort, MARCO expression was significantly higher in the RA-TB group than in healthy controls, whereas ZBTB20 did not differ significantly among the four groups.

CONCLUSIONS

Our study suggests that MARCO and ZBTB20 may act as common diagnostic biomarkers for RA and TB. The clinical findings provide preliminary support for MARCO as a marker associated with RA-TB comorbidity, whereas the translational relevance of ZBTB20 requires further validation. These findings provide a molecular framework for further investigation of RA-TB comorbidity and may inform the future development of tissue-specific biomarkers and targeted interventions.