Identification of Three Significant Genes Associated with Immune Cells Infiltration in Dysfunctional Adipose Tissue-Induced Insulin-Resistance of Obese Patients via Comprehensive Bioinformatics Analysis

Abstract

Background. Low-grade chronic inflammation in dysfunctional adipose tissue links obesity with insulin resistance through the activation of tissue-infiltrating immune cells. Numerous studies have reported on the pathogenesis of insulin-resistance. However, few studies focused on genes from genomic database. In this study, we would like to explore the correlation of genes and immune cells infiltration in adipose tissue via comprehensive bioinformatics analyses and experimental validation in mice and human adipose tissue. Methods. Gene Expression Omnibus (GEO) datasets (GSE27951, GSE55200, and GSE26637) of insulin-resistant individuals or type 2 diabetes patients and normal controls were downloaded to get differently expressed genes (DEGs), and GO and KEGG pathway analyses were performed. Subsequently, we integrated DEGs from three datasets and constructed commonly expressed DEGs’ PPI net-works across datasets. Center regulating module of DEGs and hub genes were screened through MCODE and cytoHubba in Cytoscape. Three most significant hub genes were further analyzed by GSEA analysis. Moreover, we verified the predicted hub genes by performing RT qPCR analysis in animals and human samples. Besides, the relative fraction of 22 immune cell types in adipose tissue was detected by using the deconvolution algorithm of CIBERSORT (Cell Type Identification by Estimating Relative Subsets of RNA Transcripts). Furthermore, based on the significantly changed types of immune cells, we performed correlation analysis between hub genes and immune cells. And, we performed immunohistochemistry and immunofluorescence analysis to verify that the hub genes were associated with adipose tissue macrophages (ATM). Results. Thirty DEGs were commonly expressed across three datasets, most of which were upregulated. DEGs mainly participated in the process of multiple immune cells’ infiltration. In protein-protein interaction network, we identified CSF1R, C1QC, and TYROBP as hub genes. GSEA analysis suggested high expression of the three hub genes was correlated with immune cells functional pathway’s activation. Immune cell infiltration and correlation analysis revealed that there were significant positive correlations between TYROBP and M0 macrophages, CSF1R and M0 macrophages, Plasma cells, and CD8 T cells. Finally, hub genes were associated with ATMs infiltration by experimental verification. Conclusions. This article revealed that CSF1R, C1QC, and TYROBP were potential hub genes associated with immune cells’ infiltration and the function of proinflammation, especially adipose tissue macrophages, in the progression of obesity-induced diabetes or insulin-resistance.