Identification of a Novel Diagnosis Model based on 5 Hub Genes for Chronic Thromboembolic Pulmonary Hypertension

Abstract

Background: As a type of precapillary pulmonary hypertension, chronic thromboembolic pulmonary hypertension (CTEPH) results from incomplete pulmonary embolism resolution. In this study, we aimed to determine biomarker genes for predicting the prognosis of CTEPH. Method: RNAseq of CTEPH was collected from the public database, namely Gene Expression Omnibus (GEO), including GSE84538 and GSE188938, which combined a dataset (GSE). Differentially expressed genes (DEG) or miRNA (DEM) were identified by limma package. Functional enrichment analysis was performed by the WebGestaltR package. Then, the miRNA-mRNA network was presented by Cytoscape, and the protein-protein interactions (PPI) network was constructed by STRING. MCODE was mined by mature MCODE algorithm. Immune infiltration analysis was conducted by ESTIMATER and ssGSEA analysis. A diagnosis model was established by SVM algorithm. Result: In the GSE dataset, CTEPH samples had a lower GOBP_RESPONSE_- TO_OXIDATIVE_STRESS score. A total of 628 DEGs and 31 DEMs were identified between CTEPH and normal samples. Afterward, DEGs were intersected with genes, which correlated with the GOBP_RESPONSE_TO_OXIDATIVE_STRESS score. A 26 DEMs-152 DEGs network was constructed, and a PPI network was established based on 152 DEGs to find 149 target genes. From the above 149 target genes, 3 modules were extracted to obtain 15 core targets. Finally, 5 hub genes were obtained by the intersection of 15 core targets and genes in MCODE2. A total of 5 hub genes were positively correlated with most immune cell scores as well as GOBP_RESPONSE_TO_OXIDATIVE_ STRESS. It was found that a diagnosis model based on 5 hub genes had a well diagnostic ability for CTEPH. Conclusion: We identified 5 hub genes associated with oxidative stress. It can be concluded that they may be beneficial in diagnosing CTEPH.