Identification and Validation of Genes Related to Macrophage Polarization and Cell Death Modes Under Mycobacterium tuberculosis Infection.

Abstract

To investigate the correlation between M1/M2 macrophages (M1/M2 Mφ) and cell death mode under Mycobacterium tuberculosis (Mtb) infection. Raw gene expression profiles were collected from the Gene Expression Omnibus (GEO) database. Genes related to different cell death modes were collected from the KEGG, FerrDb and GSEA databases. The differentially expressed genes (DEGs) of the gene expression profiles were identified using the limma package in R. The intersection genes of M1/M2 Mφ with different cell death modes were obtained by the VennDiagram package. Hub genes were obtained by constructing the protein-protein interactions (PPI) network and Receiver Operating Characteristic (ROC) curve analysis. The expression of cell death modes marker genes and Hub genes were verified by Western Blot and Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR). Bioinformatics analysis was performed to screen Hub genes of Mtb-infected M1 Mφ and different cell death modes, naming NFKB1, TNF, CFLAR, TBK1, IL6, RELA, SOCS1, AIM2; Hub genes of Mtb-infected M2 Mφ and different cell death modes, naming TNF, BIRC3, MAP1LC3C, DEPTOR, UVRAG, SOCS1. Combined with experimental validation, M1 Mφ under Mtb infection showed higher expression of death (including apoptosis, autophagy, ferroptosis, and pyroptosis) genes compared to M2 Mφ and genes such as NFKB1, TNF, CFLAR, TBK1, IL6, RELA, AIM2, BIRC3, DEPTOR show differential expression. NFKB1, TNF, CFLAR, TBK1, IL6, RELA, AIM2 in Mtb-infected M1 Mφ, and TNF, BIRC3, DEPTOR in Mtb-infected M2 Mφ might be used as potential diagnostic targets for TB. At early stage of Mtb infection, apoptosis, autophagy, ferroptosis, and pyroptosis occurred more significantly in M1 Mφ than that in M2 Mφ, which may contribute to the transition of Mtb-infected Mφ from M1-dominant to M2-dominant and contribute to the immune escape mechanisms of Mtb.