Abstract
Systemic sclerosis (SSc) is the most common connective tissue disease causing pulmonary hypertension (PAH). However, the cause and potential immune molecular events associated with PAH are still unclear. Therefore, it is particularly essential to analyze the changes in SSc-PAH-related immune cells and their immune-related genes. Three microarray datasets (GSE22356, GSE33463, and GSE19617) were obtained by the Gene Expression Omnibus (GEO). Compared with SSc, we found neutrophils have a statistically higher abundance, while T-cell CD4 naive and T-cell CD4 memory resting have a statistically lower abundance in peripheral blood mononuclear cells (PBMCs). Moreover, the results of Gene Set Enrichment Analysis (GSEA) showed there is a differential enrichment of multiple pathways between SSc and SSc-PAH. By combining differentiated expressed genes (DEGs) and immune-related genes (IRGs), fifteen IRGs were selected. In addition, we also analyzed the first five rich Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and the most abundant Gene Ontology (GO)-molecular functional terms. Furthermore, interleukin-7 receptor (IL-7R), tyrosine-protein kinase (LCK), histone deacetylase 1 (HDAC1), and epidermal growth factor receptor (EGFR) genes were identified as hub genes via protein-protein interaction (PPI) network analysis. The Comparative Toxic Genomics Database (CTD) analysis result showed that LCK, HDAC1, and EGFR have a higher score with SSc. Coexpression network analysis confirmed that IL-7R, LCK, and HDAC1 are key genes related to immune regulation in SSc without PAH and are involved in T-cell immune regulation. Subsequently, using GSE22356 and GSE33463 as the test sets and GSE19617 as the verification set, it was verified that the mRNA expression levels of the three central genes of SSc-PAH were significantly lower than those of the SSc without PAH samples. Consistent with previous predictions, the expressions of IL-7R, LCK, and HDAC1 are positively correlated with the numbers of T-cell CD4 naive and T-cell CD4 memory, while the expressions of IL-7R and LCK are negatively correlated with the numbers of neutrophils in the peripheral blood. Therefore, this evidence may suggest that these three immune-related genes: IL-7R, LCK, and HDAC1, may be highly related to the immunological changes in SSc-PAH. These three molecules can reduce T cells in SSc-PAH PBMCs through the regulation of T-cell activation, which suggests that these three molecules may be involved in the development of SSc-PAH. Meanwhile, the low expression of IL-7R, LCK, and HDAC1 detected in the peripheral blood of SSc may indicate the possibility of PAH and hopefully become a biomarker for the early detection of SSc-PAH. Finally, 49 target miRNAs of 3 specifically expressed hub genes were obtained, and 49 mRNA-miRNA pairs were identified, which provided directions for our further research.