# Core autophagy genes and immune infiltration characteristics in rheumatoid arthritis: A bioinformatics study

**Authors:** Zining Peng, Qian Deng, Yuanbo Huang, Fanyu Meng, Yuan Long, Yuanyuan Wei, Weitian Yan, Xiaoyu Zhang, Jiangyun Peng, Zhaofu Li, Nian Liu

PMC · DOI: 10.1371/journal.pone.0326168 · 2025-07-11

## TL;DR

This study identifies 11 autophagy-related genes linked to rheumatoid arthritis and finds they correlate with immune cell activity in joints, suggesting new treatment possibilities.

## Contribution

Systematic identification of 11 RA-autophagy core genes with immune infiltration correlations and preclinical validation in a CIA rat model.

## Key findings

- 11 core autophagy genes (e.g., IFNG, FOXO1) were identified as potential RA therapeutic targets.
- RA synovium showed 14 differentially abundant immune cell types compared to healthy controls.
- Animal experiments confirmed significant changes in autophagy and mRNA levels of core genes.

## Abstract

This study aims to identify autophagy-related biomarkers in rheumatoid arthritis (RA) synovium, analyze their immune infiltration characteristics, and validate therapeutic potential through multi-level experimental approaches.

We used public datasets to obtain synovial tissue genes of healthy people and RA patients, screened differentially expressed genes (DEGs) of RA, and intersected with the human autophagy gene database (HADb) to obtain RA autophagy genes. GO and KEGG enrichment analysis and single-gene genome enrichment analysis were performed. The diagnostic value of RA core autophagy genes in the validation set was screened and verified. The immune cell correlation analysis of RA autophagy core genes was performed to obtain the correlation between single disease autophagy core genes and immune cells. Finally, we prepared CIA rat models to verify the autophagy protein P62, Beclin-1 and the 11 core genes associated with RA-autophagy.

A total of 1098 RA DEGs were obtained. Intersecting with 222 autophagy genes obtained from the HADb database yielded 27 RA autophagy genes. GO analysis of RA autophagy genes showed 307 biological mechanisms. KEGG enrichment analysis obtained 86 signaling pathways, including FoxO, Necroptosis and other pathways related to RA autophagy. GSEA analysis found that the control group had a higher correlation with adipokine signaling pathways and others. And 11 RA autophagy-related core genes (IFNG, EGFR, MYC, CXCR4, MAPK8, CASP1, TNFSF10, CTSB, FAS, FOXO1, FOXO3) were obtained by screening the PPI network, and there were differences in expression in the training set (P < 0.001). External validation set verification showed diagnostic efficacy. Analysis of immune infiltration in RA autophagy-related genes revealed 14 immune cell types differentially abundant in synovial tissues of RA patients vs. normal controls. Significant correlations exist between autophagy genes and immune subsets. Finally, animal experiments showed that joint autophagy was enhanced (P < 0.001), and the mRNA of 11 RA-autophagy core genes had significant changes (P < 0.001).

We systematically identified 11 autophagy-related core genes as potential therapeutic targets for RA. Our CIA model validation provides preclinical evidence supporting their translational potential. These genes showed significant correlations with 14 synovial immune cell subtypes, may serve as novel therapeutic targets by modulating immune infiltration and inflammatory pathways. Future investigations should focus on elucidating the mechanistic basis of the observed gene-immune cell interactions in both autophagic and immune pathways to facilitate the development of precision therapies.

## Linked entities

- **Genes:** IFNG (interferon gamma) [NCBI Gene 3458], EGFR (epidermal growth factor receptor) [NCBI Gene 1956], MYC (MYC proto-oncogene, bHLH transcription factor) [NCBI Gene 4609], CXCR4 (C-X-C motif chemokine receptor 4) [NCBI Gene 7852], MAPK8 (mitogen-activated protein kinase 8) [NCBI Gene 5599], CASP1 (caspase 1) [NCBI Gene 834], TNFSF10 (TNF superfamily member 10) [NCBI Gene 8743], CTSB (cathepsin B) [NCBI Gene 1508], FAS (Fas cell surface death receptor) [NCBI Gene 355], FOXO1 (forkhead box O1) [NCBI Gene 2308], FOXO3 (forkhead box O3) [NCBI Gene 2309]
- **Proteins:** GTF2H1 (general transcription factor IIH subunit 1), BECN1 (beclin 1)
- **Diseases:** rheumatoid arthritis (MONDO:0008383)
- **Species:** Mus musculus (taxon 10090)

## Full-text entities

- **Genes:** NUP62 (nucleoporin 62) [NCBI Gene 23636] {aka IBSN, SNDI, p62}, CXCR4 (C-X-C motif chemokine receptor 4) [NCBI Gene 7852] {aka CD184, D2S201E, FB22, HM89, HSY3RR, LCR1}, FOXO3 (forkhead box O3) [NCBI Gene 2309] {aka AF6q21, FKHRL1, FKHRL1P2, FOXO2, FOXO3A}, MAPK8 (mitogen-activated protein kinase 8) [NCBI Gene 5599] {aka JNK, JNK-46, JNK1, JNK1A2, JNK21B1/2, PRKM8}, IFNG (interferon gamma) [NCBI Gene 3458] {aka IFG, IFI, IMD69}, CTSB (cathepsin B) [NCBI Gene 1508] {aka APPS, CPSB, KWE, RECEUP}, CASP1 (caspase 1) [NCBI Gene 834] {aka ICE, IL1BC, P45}, EGFR (epidermal growth factor receptor) [NCBI Gene 1956] {aka ERBB, ERBB1, ERRP, HER1, NISBD2, NNCIS}, MYC (MYC proto-oncogene, bHLH transcription factor) [NCBI Gene 4609] {aka MRTL, MYCC, bHLHe39, c-Myc}, TNFSF10 (TNF superfamily member 10) [NCBI Gene 8743] {aka APO2L, Apo-2L, CD253, TANCR, TL2, TNLG6A}, FAS (Fas cell surface death receptor) [NCBI Gene 355] {aka ALPS1A, APO-1, APT1, CD95, FAS1, FASTM}, BECN1 (beclin 1) [NCBI Gene 8678] {aka ATG6, VPS30, beclin1}, FOXO1 (forkhead box O1) [NCBI Gene 2308] {aka FKH1, FKHR, FOXO1A}
- **Diseases:** RA (MESH:D001172), inflammatory (MESH:D007249)
- **Species:** Homo sapiens (human, species) [taxon 9606], Rattus norvegicus (brown rat, species) [taxon 10116]

## Figures

13 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12250634/full.md

---
Source: https://tomesphere.com/paper/PMC12250634