Identification of immune-related hub genes and analysis of infiltrated immune cells of idiopathic pulmonary artery hypertension

OBJECTIVES: Idiopathic pulmonary artery hypertension (IPAH) is a rare but life-threaten disease. However, the mechanism underlying IPAH is unclear. In this study, underlying mechanism, infiltration of immune cells, and immune-related hub genes of IPAH were analyzed via bioinformatics. METHODS: GSE15197, GSE48149, GSE113439, and GSE117261 were merged as lung dataset. Weighted gene correlation network analysis (WGCNA) was used to construct the co-expression gene networks of IPAH. Gene Ontology and pathway enrichment analysis were performed using DAVID, gene set enrichment analysis (GSEA), and gene set variation analysis (GSVA). Infiltration of immune cells in lung samples was analyzed using CIBERSORT. GSE22356 and GSE33463 were merged as peripheral blood mononuclear cells (PBMCs) dataset. Immune-related differentially expressed genes (IRDEGs) of lung and PBMCs dataset were analyzed. Based on the intersection between two sets of IRDEGs, hub genes were screened using machine learning algorithms and validated by RT-qPCR. Finally, competing endogenous RNA (ceRNA) networks of hub genes were constructed. RESULTS: The gray module was the most relevant module and genes in the module enriched in terms like inflammatory and immune responses. The results of GSEA and GSVA indicated that increasement in cytosolic calcium ion, and metabolism dysregulation play important roles in IPAH. The proportions of T cells CD4 memory resting and macrophage M1 were significantly greater in IPAH group, while the proportions of monocytes and neutrophils were significantly lower in IPAH group. IRDEGs of two datasets were analyzed and the intersection between two set of IRDEGs were identified as candidate hub genes. Predictive models for IPAH were constructed using data from PBMCs dataset with candidate hub genes as potential features via LASSO regression and XGBoost algorithm, respectively. CXCL10 and VIPR1 were identified as hub genes and ceRNA networks of CXCL10 was constructed. CONCLUSION: Inflammatory response, increasement in cytosolic calcium ion, and metabolism dysregulation play important roles in IPAH. T cells CD4 memory resting and macrophage M1 were significantly infiltrated in lung samples from patients with IPAH. IRDEGs of lung dataset and PBMCs dataset were analyzed, and CXCL10 and VIPR1 were identified as hub genes.