Idiopathic pulmonary arterial hypertension (IPAH) is a life-threatening cardiopulmonary disease lacking specific diagnostic markers and targeted therapy, and its mechanism of development remains to be elucidated. The present study aimed to explore novel diagnostic biomarkers and therapeutic targets in IPAH by integrated bioinformatics analysis. Four eligible datasets (GSE117261, GSE15197, GSE53408, GSE48149) was firstly downloaded from GEO database and subsequently integrated by Robust rank aggregation (RRA) method to screen robust differentially expressed genes (DEGs). Then functional annotation of robust DEGs was performed by GO and KEGG enrichment analysis. The protein-protein interaction (PPI) network was constructed followed by using MCODE and CytoHubba plug-in to identify hub genes. Finally, 10 hub genes were screened including ENO1, TALDO1, TXNRD1, SHMT2, IDH1, TKT, PGD, CXCL10, CXCL9, and CCL5. The GSE113439 dataset was used as a validation cohort to appraise these hub genes and TXNRD1 was selected for verification at the protein level. The experiment results confirmed that serum TXNRD1 concentration was lower in IPAH patients and the level of TXNRD1 had great predictive efficiency (AUC:0.795) as well as presents negative correlation with mean pulmonary arterial pressure (mPAP) and pulmonary vascular resistance (PVR). Consistently, the expression of TXNRD1 was proved to be inhibited in animal and cellular model of PAH. In addition, GSEA analysis was performed to explore the functions of TXNRD1 and the results revealed that TXNRD1 was closely correlated with mTOR signaling pathway, MYC targets, and unfolded protein response. Finally, knockdown of TXNRD1 was shown to exacerbate proliferative disorder, migration and apoptosis resistance in PASMCs. In conclusion, our study demonstrates that TXNRD1 is a promising candidate biomarker for diagnosis of IPAH and plays an important role in PAH pathogenesis, although further research is necessary.