Bioinformatics analysis of the immune cell infiltration characteristics and correlation with crucial diagnostic markers in pulmonary arterial hypertension

Abstract

Abstract Background Pulmonary arterial hypertension (PAH) is a pathophysiological syndrome, characterized by pulmonary vascular remodeling. Immunity and inflammation are progressively recognized properties of PAH, which are crucial for the initiation and maintenance of pulmonary vascular remodeling. This study explored immune cell infiltration characteristics and potential biomarkers of PAH using comprehensive bioinformatics analysis. Methods Microarray data of GSE117261, GSE113439 and GSE53408 datasets were downloaded from Gene Expression Omnibus database. The differentially expressed genes (DEGs) were identified in GSE117261 dataset. The proportions of infiltrated immune cells were evaluated by CIBERSORT algorithm. Feature genes of PAH were selected by least absolute shrinkage and selection operator (LASSO) regression analysis and validated by fivefold cross-validation, random forest and logistic regression. The GSE113439 and GSE53408 datasets were used as validation sets and logistic regression and receiver operating characteristic (ROC) curve analysis were performed to evaluate the prediction value of PAH. The PAH-associated module was identified by weighted gene association network analysis (WGCNA). The intersection of genes in the modules screened and DEGs was used to construct protein–protein interaction (PPI) network and the core genes were selected. After the intersection of feature genes and core genes, the hub genes were identified. The correlation between hub genes and immune cell infiltration was analyzed by Pearson correlation analysis. The expression level of LTBP1 in the lungs of monocrotaline-induced PAH rats was determined by Western blotting. The localization of LTBP1 and CD4 in lungs of PAH was assayed by immunofluorescence. Results A total of 419 DEGs were identified, including 223 upregulated genes and 196 downregulated genes. Functional enrichment analysis revealed that a significant enrichment in inflammation, immune response, and transforming growth factor β (TGFβ) signaling pathway. CIBERSORT analysis showed that ten significantly different types of immune cells were identified between PAH and control. Resting memory CD4+ T cells, CD8+ T cells, γδ T cells, M1 macrophages, and resting mast cells in the lungs of PAH patients were significantly higher than control. Seventeen feature genes were identified by LASSO regression for PAH prediction. WGCNA identified 15 co-expression modules. PPI network was constructed and 100 core genes were obtained. Complement C3b/C4b receptor 1 (CR1), thioredoxin reductase 1 (TXNRD1), latent TGFβ binding protein 1 (LTBP1), and toll-like receptor 1 (TLR1) were identified as hub genes and LTBP1 has the highest diagnostic efficacy for PAH (AUC = 0.968). Pearson correlation analysis showed that LTBP1 was positively correlated with resting memory CD4+ T cells, but negatively correlated with monocytes and neutrophils. Western blotting showed that the protein level of LTBP1 was increased in the lungs of monocrotaline-induced PAH rats. Immunofluorescence of lung tissues from rats with PAH showed increased expression of LTBP1 in pulmonary arteries as compared to control and LTBP1 was partly colocalized with CD4+ cells in the lungs. Conclusion LTBP1 was correlated with immune cell infiltration and identified as the critical diagnostic maker for PAH.