Analysis of CD4+ T-helper-associated hub gene signature and immune dysregulation via RNA-sequencing data in a mouse tail model of lymphedema

Gland Surg. 2023 Sep 25;12(9):1141-1157. doi: 10.21037/gs-23-48. Epub 2023 Sep 20.

Abstract

Background: T-helper cells play an essential role in the progression of lymphedema. This study aimed to explore the biological significance of T-helper cell-associated genes (THAGs) in a mouse tail model of lymphedema by RNA-sequencing (RNA-seq) data.

Methods: The expression profiles of a murine model of secondary lymphedema were obtained from European Nucleotide Archive (ENA) database. Differentially expressed genes (DEGs) were screened and the enrichment analysis of DEGs was conducted. THAGs were constructed by crossing the T-helper-related gene sets obtained from Molecular Signatures Database with DEGs. Protein-protein interaction (PPI) network analysis was utilized to establish T-helper-associated hub genes (THAHGs). Single-sample gene set enrichment analysis (ssGSEA) was employed to decipher differences in immune cell infiltration. The correlation between THAHGs and immune infiltration was calculated by Pearson correlation analysis. Receiver operating characteristic (ROC) curves of THAHGs were drawn to evaluate their diagnostic properties. Additionally, potential drugs and upstream transcription factors (TFs) were predicted based on THAHGs.

Results: Enrichment analysis showed that lymphedematous tissue presented higher activation of biological process (BP) of T-helper 1 (Th1), T-helper 2 (Th2), T-helper 17 (Th17). The immune infiltration analysis further calculated that the relative immune abundance of follicular B cells, memory B cells, M1 macrophage, and CD4+ Tm cells was significantly elevated while the relative immune abundance of neutrophils and plasma cells were down-regulated in lymphedema. We established a list of THAHGs consisting of eight hub genes, compassing Cd4, Foxp3, Irf4, Ccr6, Il12rb1, Batf, Il1b, Cd74. THAHGs were shown to be significantly interrelated and related to immune infiltration by Pearson correlation analysis. ROC curves showcased that the area under curve (AUC) values of THAHGs were larger than 0.70. Gata3 was the most potential TF and thalidomide might be the immunoregulatory drug for lymphedema based on THAHGs.

Conclusions: Biological pathways associated with T-helpers were significantly enriched in mouse lymphedema tissue. The relative immune infiltration abundance of M1 macrophage, CD4+ Tm cells, and T-helper cells was higher in the lymphedema group. Besides, we identified the THAHGs containing eight genes, namely, Cd4, Foxp3, Irf4, Ccr6, Il12rb1, Batf, Il1b, and Cd74. The THAHGs were closely correlated with immune infiltration results and with good diagnostic properties.

Keywords: CD4; Foxp3; Lymphedema; T-helper 17 (Th17); T-helper cell.