- Research
- Open access
- Published:
Identification of Ferroptosis-related potential biomarkers and immunocyte characteristics in Chronic Thromboembolic Pulmonary Hypertension via bioinformatics analysis
BMC Cardiovascular Disorders volume 23, Article number: 504 (2023)
Abstract
Background
Chronic Thromboembolic Pulmonary Hypertension (CTEPH) is a form of pulmonary hypertension with a high mortality rate. A new type of iron-mediated cell death is Ferroptosis, which is characterized by the accumulation of lethal iron ions and lipid peroxidation leading to mitochondrial atrophy and increased mitochondrial membrane density. Now, there is a lack of Ferroptosis-related biomarkers (FRBs) associated with pathogenic process of CTEPH.
Methods
The differentially expressed genes (DEGs) of CTEPH were obtained by GEO2R. Genes related to Ferroptosis were obtained from FerrDb database. The intersection of Ferroptosis and DEGs results in FRBs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed in Database for Annotation, Visualization and Integrated Discovery (DAVID) database. The optimal potential biomarkers for CTEPH were analyzed by least absolute shrinkage and selection operator (LASSO) and support vector machine-recursive feature elimination (SVM-RFE) machine learning. The four hub genes were verified from the Gene Expression Omnibus (GEO) dataset GSE188938. Immune infiltration was analyzed by CIBERSORT. SPSS software was used to analyze the Spearman rank correlation between FRBs identified and infiltration-related immune cells, and p < 0.05 was considered as statistically significant.
Results
In this study, potential genetic biomarkers associated with Ferroptosis in CTEPH were investigated and explored their role in immune infiltration. In total, we identified 17 differentially expressed Ferroptosis-associated genes by GEOquery package. The key FRBs including ARRDC3, HMOX1, BRD4, and YWHAE were screened using Lasso and SVM-RFE machine learning methods.Through gene set GSE188938 verification, only upregulation of gene ARRDC3 showed statistical difference. In addition, immune infiltration analysis using the CIBERSORT algorithm revealed the infiltration of Eosinophils and Neutrophils in CTEPH samples was less than that in the control group. And correlation analysis revealed that ARRDC3 was positively correlated with T cells follicular helper (r = 0.554, p = 0.017) and negatively correlated with Neutrophils (r = -0.47, p = 0.049).
Conclusions
In conclusion, ARRDC3 upregulation with different immune cell infiltration were involved in the development of CTEPH. ARRDC3 might a potential Ferroptosis-related biomarker for CTEPH treatment. This study provided a new insight into pathogenesis CTEPH.
Introduction
Chronic thromboembolic pulmonary hypertension (CTEPH) is a rare progressive pulmonary vascular disease that usually results from acute pulmonary embolism. CTEPH is characterized by persistent thromboembolic pulmonary artery obstruction and the presence of variable small vessel lesions [1]. In recent years, studies have shown that between 3 and 30 people per 1 million develop CTEPH [2]. The early symptoms of CTEPH are usually not obvious. The mortality of patients with advanced CTEPH is high, and the recent 5-year survival rate is less than 40% [3]. The pathogenesis of CTEPH is believed to be caused by multiple factors, including hypercoagulability, fibrinolysis disease, inflammation, pulmonary microvascular hyperplasia [4]. CTEPH is difficult to diagnose in early stages, and biomarkers to predict its severity and clinical outcome are still lacking. As a result, it is crucial to explore more valuable genetic biomarkers for the potential etiology and therapeutic targets of CTEPH.
Ferroptosis is characterized by mitochondrial atrophy and increased mitochondrial membrane density, the accumulation of iron and lipid reactive oxygen species and the involvement of a unique set of genes. Divalent iron ions in the unstable iron pool of the cell generate hydroxyl radicals (lipid ROS) through the Fenton reaction, causing peroxidation of polyunsaturated fatty acids. This process is accompanied by overloading or inactivation of the lipid peroxide scavenging system, and eventually the accumulation of toxic lipid metabolites leads to cell death [5]. The pulmonary circulation has the physiological basis for Ferroptosis because it is rich in blood oxygen supply, polyunsaturated fatty acids and iron ions. Although Ferroptosis was first formally proposed in tumors, more and more studies have shown that Ferroptosis is involved in the development of chronic obstructive pulmonary disease(COPD), acute lung injury, pulmonary fibrosis and other lung diseases [6]. Ferroptosis is involved in the damage of pulmonary vessels through oxidative stress and lipid peroxidation [7]. Recently, a growing number of studies have focused on discovering biomarkers of Ferroptosis at multiple levels, including morphology, biochemistry, protein, and genes [8]. However, it remains a great challenge to explore the reliable gene biomarkers and specific regulatory details associated with Ferroptosis in CTEPH.
Bioinformatics analysis is a new cross analysis combining molecular biology and information technology. The aim is to efficiently and extensively obtain differential gene expression of diseases by using gene chips, and provide signaling pathways for exploring the specific mechanisms of diseases or prediction for drug therapy [9,10,11,12]. Immune cells play an important role in the pathological process of CTEPH by participating in the release of inflammatory factors, cytotoxic effects, thrombosis and so on [13]. In the present study, we applied the CIBERSORT tool to the assessment of immune infiltration in the pathological process of CTEPH. Besides, we screened and verified Ferroptosis-related genes as CTEPH biomarkers, analyzed their role in the pulmonary vascular immune microenvironment by bioinformatics analysis. And, the pathway analysis was carried out to provide specific pathway guidance for the subsequent experimental research.
Materials and methods
Data download
We used the “GEOquery” package of R software (version 4.1.3) to download the CTEPH expression profile datasets GSE130391 from the Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/) database (Supplementary Material 1) [14]. Details of the GSE130391 dataset: Samples surgically extracted from CTEPH pulmonary arteries were compared to post-transplant pulmonary arteries from IPAH and failed donor pulmonary arteries from controls. The goal was to determine if altered gene expression in CTEPH associated with development of chronic thrombi. The CTEPH group (n = 14) and the Control group (n = 4) were selected for subsequent analysis.
Differential expression analysis
We first extracted 676 differentially expressed genes (DEGs) between CTEPH and Control from GSE130391 database using GEO2R (Supplementary Material 2). DEGs with p < 0.05 and |log2FC|> 1 were considered significant. Genes related to Ferroptosis were obtained from FerrDb database [15]. The intersection genes of Ferroptosis and DEGs were considered to be the most important genes. Venn diagram was drawn by Venn diagram intersection (bioinformatics.psb.ugent.be/webtools/Venn/).
Functional enrichment
GO and KEGG analysis of 17 genes were performed in Database for Annotation, Visualization and Integrated Discovery (DAVID) [16]. Metascape (http://metas.cape.org/) was used to analyse Reactome pathway enrichment and Immunologic Signatures enrichment [17].The KEGG formal permission had been obtained from Kanehisa laboratories, which copyrighted KEGG pathway database [18,19,20]. p < 0.05 was considered significant enrichment. Gene set enrichment analysis (GSEA) was performed on the gene expression matrix through Xiangtao Tool [21]. p < 0.05 was considered significant enrichment.
Identification of optimal potential biomarkers for CTEPH
Glmnet package of R was applied for the LASSO algorithm. And the average misjudgment rate was compared by its tenfold cross validation in addition. SVM-RFE is a machine learning method based on support vector machine (SVM), which deletes the feature vectors generated by SVM to find the best variables. “e1071” package was used to identify the diagnostic value of these biomarkers in CTEPH [22]. The area under the curve (AUC) was used to evaluate the diagnostic ability of the best genetic biomarkers by SVM-RFE. The intersection genes of Lasso and SVM-RFE were considered as the best diagnostic gene biomarker for CTEPH. The code for Lasso and SVM-RFE machine learning algorithms has been shown in Supplementary Material 1.
Verification of optimal potential biomarkers for CTEPH
The CTEPH-related dataset GSE188938 was obtained from the GEO database. R software was used to obtain gene expression profiles(Control, n = 5;CTEPH, n = 7). In the study of GSE188938, Professor Xu sought to identify key genes and pathways that contribute to the diagnosis and treatment of CTEPH. It provided the research direction for further understanding the molecular mechanism of CTEPH [23]. SPSS software was used for difference analysis between Control group and CTEPH group.
Immune infiltration analysis
We uploaded GSE130391 gene expression matrix data to CIBERSORT (https://cibersort.stanford.edu/), and obtained 22 kinds of immune cell infiltration information [24]. p < 0.05 of the samples were identified significantly.“corrplot” package and “ggplot2” package were used to draw a correlation heatmap and boxplot. The proportion of 22 types of infiltrating immune cell types in each tissue from the GSE130391 dataset were showed in Supplementary Material 5.
Correlation analysis between diagnostic biomarkers and infiltration-related immune cells
SPSS software was used to analyze the Spearman rank correlation between FRBs identified and the level of related immune cells. P < 0.05 was considered statistically significant [25]. We also uesed "ggplot2" package in R software to diplay the correlation.
Statistical analysis
The data were analyzed by one-way Analysis of Variance(ANOVA) to determine differences between the Control group and the CTEPH group. For all statistical tests, p < 0.05 was considered as statistically significant.
Results
Identification of Ferroptosis related genes and DEGs in CTEPH
There were 676 differentially expressed genes (DEGs) between CTEPH and Control from GSE130391 database. PCA diagram showed the expression of these genes (Fig. 1A). 259 genes related to Ferroptosis were obtained from FerrDb database (Supplementary Material 3). There were 17 intersection genes of Ferroptosis and DEGs (Fig. 1B, Table 1). The clustering heatmap showed the expression pattern of 17 Ferroptosis related biomarkers (FRBs) among samples (Fig. 1C). The correlation between these 17 genes were shown in Fig. 1D.
Functional enrichment of 17 Ferroptosis-related DEGs
GO enrichment analysis indicated that the above 17 Ferroptosis related biomarkers were significantly related to regulation of inflammatory response, Cell growth cycle, cell death, iron ion homeostasis, mitophagy and so on (Fig. 2A). KEGG analysis showed 17 FRBs significantly enriched in Leishmaniasis, Pathways in cancer, Cell cycle, etc. (Fig. 2B). Reactome pathway analysis indicated that Cellular responses to stress and Signaling by Rho GTPases were enriched (Fig. 2C). Surprisingly, the 17 FRBs were also obviously enriched in many immune-related signatures (Fig. 2D). It was showed the top five GSEA enrichment gene sets involved in inflammation, CD5 cells, natural killer cells and B lymphocyte in Fig. 2E. Collectively, These evidences indicated that 17 intersection genes mentioned may play an important role in the pathogenesis of CTEPH by participating in the regulation of immune cells and cytokines.
Screening and verification of potential biomarkers for CTEPH
We screened out 5 genes from 17 FRBs as the most suitable potential biomarkers for CTEPH by using LASSO algorithm (Fig. 3A-B). 13 potential biomarkers were selected from DEGs by SVM-RFE algorithm. AUC = 0.788 indicated that SVM-RFE had a good gene selection effect (Fig. 3C). The intersection potential biomarkers obtained by the two algorithms were considered the most accurate potential biomarkers, including ARRDC3, HMOX1, BRD4 and YWHAE (Fig. 3D).
Verification of four hub genes from the GEO dataset GSE188938
In the GEO dataset GSE188938, ARRDC3 was significantly upregulated in CTEPH group compared with control group (p = 0.0051, Fig. 4A). However, we found that the gene expression of HMOX1, BRD4 and YWHAE were not statistically significant (Fig. 4B-D).These results indicated that ARRDC3 may be the most important hub gene of CTEPH and a potential therapeutic target. The expression of ARRDC3, HMOX1, BRD4 and YWHAE in GSE188938 were showed in Supplementary Material 4.
Immune infiltration analysis
The boxplot of the immune cell infiltration difference showed a reduction in Eosinophils and Neutrophils at CTEPH compared with Control samples (Fig. 5A). In correlation heatmap, Pearson correlation analysis of the 22 types of immune cells reveals that Eosinophils is negatively correlated with Neutrophils (Fig. 5B). As a result, Eosinophils and Neutrophils may be involved in the pathological regulation of CTEPH.
Correlation analysis between ARRDC3 and infiltration-related immune cells
Correlation analysis revealed that ARRDC3 was positively correlated with T cells follicular helper (r = 0.554, p = 0.017, Fig. 6A-B). ARRDC3 was negatively correlated with Neutrophils (r = -0.47, p = 0.049, Fig. 6C, Supplementary Material 5).
Discussion
Despite the advancements in CTEPH treatment during the past decades, pulmonary endarterectomy (PEA) treatment has not been completely satisfactory and there are still some patients with recurrent attacks. With the development of sequencing technology, it is increasingly expected to uncover the biological effects of genes related to CTEPH. Ferroptosis is believed to involve in the development of vascular inflammation and thrombosis. However, few studies have focused on the genes and underlying regulatory details associated with iron deficiency and the immune infiltration landscape, which have profound implications for patients with CTEPH. In this study, a total of 17 differentially expressed genes related to Ferroptosis in CTEPH patients were identified in the metadata cohort. Among them, four genes including ARRDC3, HMOX1, BRD4, and YWHAE are considered as key genes for Ferroptosis in CTEPH development by the integration of LASSO and SVM-RFE machine learning methods. Only the upregulation of ARRDC3 was found to have a statistical difference through gene set GSE188938 verification. In addition, analysis of immune infiltration using the CIBERSORT algorithm revealed less infiltration of Eosinophils and Neutrophils in CTEPH samples than in controls. ARRDC3 was found to be positively correlated with T cells follicular helper and negatively correlated with Neutrophils according to the correlation analysis.
Firstly, We downloaded the CTEPH expression profile dataset GSE130391 from the GEO database and identified a total of 676 DEGs on CTEPH. There were 17 intersection genes of Ferroptosis and DEGs, indicating that Ferroptosis is likely to be involved in the pathogenesis of CTEPH through these 17 genes. GO enrichment analysis revealed that 17 intersection genes were mainly associated with regulation of inflammatory response, regulation of cell cycle and cell death. In the results of KEGG enrichment analysis, there were significantly enriched in pathways of cancer, cell cycle, chemical carcinogenesis-reactive oxygen species. Evidences have confirmed that immune response participated in CTEPH development [26] and chronic inflammatory diseases including inflammatory bowel disease, osteomyelitis and anticardiolipin antibody syndrome and chronic venous ulcers would increase the risk of CTEPH [27]. CTEPH pulmonary vascular remodeling could be attributed partly to the abnormal gene regulation of cell cycle and cell death via apoptosis and Ferroptosis in pulmonary artery endothelial cells and smooth muscle cells [5, 28]. These studies mentioned above were consistent with our GO and KEGG results. In the present study, Immunologic Signatures enrichment analysis also indicated that Ferroptosis-related genes were involved in CTEPH pathogenesis by regulating autophagy, immune cells, cytokines and multiple kinases. In addition, the main gene sets enriched in GSEA analysis were NFLAMMATORY RESPONSE, CD5 TARGETS, NATURAL KILLER DIFFERENTIATION, CELL CYCLE MITOTIC, HADDAD B LYMPHOCYTE PROGENITOR. All above results suggested that Ferroptosis and immune inflammatory response might play important roles in CTEPH progression.
Next, four genes including ARRDC3, HMOX1, BRD4, and YWHAE were identified as key hub genes in CTEPH combined with Ferroptosis by LASSO and SVM-RFE algorithms which are mainly used to filter the most characteristic variables and construct the best classification models [26]. ARRDC3 is associated with inflammation, cellular invasion, etc., which promotes H. pylori-associated gastritis by regulating protease-activated receptor 1 and inhibits invasive metastasis of breast cancer cells by regulating G protein-coupled receptor lysosomes [29]. HMOX1 could regulate the proliferation and apoptosis of human cardiomyocytes and play an anti-oxidative stress role in apoptosis induced by high glucose [30, 31]. BRD4 inhibition prevented diabetic cardiomyopathy induced by high-fat diet through activating mediated mitochondrial autophagy [32]. YWHAE promoted proliferation and metastasis of breast cancer cells, and regulated proliferation and cell cycle protein D1 through RAF/MAPK and Hippo pathways [33]. The above studies indicated that ARRDC3, HMOX1, BRD4, and YWHAE play a role in promoting inflammation, cell proliferation and metastasis, and abnormal cell apoptosis. The pathogenesis of CTEPH is closely related to chronic inflammation, abnormal cell proliferation and abnormal cell cycle regulation [34]. In this study, our results indicated that ARRDC3, HMOX1, BRD4, and YWHAE were involved in CTEPH progress possibly through regulating inflammatory, vascular cell proliferation, smooth muscle cell invasion and apoptosis. Through gene set GSE188938 verification, only ARRDC3 upregulation showed statistical difference suggesting ARRDC3 might be the most important hub gene of CTEPH and a potential therapeutic target for CTEPH.
Additionally, we screened the proportion of 22 immune cells in CTEPH using CIBERSORT analysis, and found that there were significant differences in the cells infiltration of Eosinophils and Neutrophils compared with the control sample. We also analyzed the correlation between ARRDC3 and immune cells infiltration. ARRDC3 was positively correlated with T cells follicular helper and negatively correlated with Neutrophils. Eosinophils and Neutrophils were believed to be involved in pulmonary artery atherosclerosis and thrombotic lesions leaded to CTEPH [35, 36]. Ferroptosis-inducing agents induced Ferroptosis-like cell death of Eosinophils [37]. Neutrophil extracellular traps mediated m6A modification and regulated sepsis-associated acute lung injury by activating Ferroptosis in alveolar epithelial cells [38]. These studies were similar with our results that both Eosinophils and Neutrophils were involved in Ferroptosis processes of CTEPH. Our study found that Ferroptosis-related hub gene ARRDC3 was related to immune cell infiltration including T cells follicular helper and Neutrophils, and further experiments are necessary to verify its exact effect and underlying regulatory mechanism.
The results of bioinformatics analysis provided a new clue and viewpoint for our further research on CTEPH. Undeniably, there were some limitations in the current study, such as the sample size of the selected database, heterogeneity of the clinical parameters of CTEPH patients and lack of grading of disease severity. To minimize sample heterogeneity, it is necessary to include more samples and improve clinical data. Further experimental verification are warranted to confirm the effects and functional correlation of the identified genes and pathways, and provide more accurate characterization of immune cell populations in CTEPH in vivo and vitro studies.
In conclusion, ARRDC3 upregulation with different immune cell infiltration were involved in the development of CTEPH and ARRDC3 might a potential and promising Ferroptosis-related biomarker for CTEPH treatment. Our study might provide new insights into the prevention, monitoring and potential therapeutic intervention for CTEPH.
Availability of data and materials
The data used to support the findings of this study are available from Supplementary Material.
Abbreviations
- CTEPH:
-
Chronic thromboembolic pulmonary hypertension
- FRBs:
-
Ferroptosis related biomarkers
- DEGs:
-
Differentially expressed genes
- GO:
-
Gene ontology
- KEGG:
-
Kyoto encyclopedia of genes and genomes
- GSEA:
-
Gene set enrichment analysis
- GEO:
-
Gene expression omnibus
- LASSO:
-
Least absolute shrinkage and selection operator
- SVM-RFE:
-
Support vector machine-recursive feature elimination
- COPD:
-
Chronic obstructive pulmonary disease
- DAVID:
-
Database for annotation, visualization and integrated discovery
References
Teerapuncharoen K, Bag R. Chronic Thromboembolic Pulmonary Hypertension. Lung. 2022;200(3):283–99.
Martinez C, Wallenhorst C, Teal S, Cohen AT, Peacock AJ. Incidence and risk factors of chronic thromboembolic pulmonary hypertension following venous thromboembolism, a population-based cohort study in England. Pulm Circ. 2018;8:767756354.
Yandrapalli S, Tariq S, Kumar J, Aronow WS, Malekan R, Frishman WH, et al. Chronic thromboembolic pulmonary hypertension: epidemiology, diagnosis, and management. Cardiol Rev. 2018;26:62–72.
Mullin CJ, Klinger JR. Chronic thromboembolic pulmonary Hypertension. Heart Fail Clin. 2018;14:339–51.
Jiang X, Stockwell BR, Conrad M. Ferroptosis: mechanisms, biology and role in disease. Nat Rev Mol Cell Biol. 2021;22:266–82.
Tao N, Li K, Liu J. Molecular mechanisms of ferroptosis and its role in pulmonary disease. Oxid Med Cell Longev. 2020;2020:9547127.
Mao H, Zhao Y, Li H, Lei L. Ferroptosis as an emerging target in inflammatory diseases. Prog Biophys Mol Biol. 2020;155:20–8.
Tang D, Chen X, Kang R, Kroemer G. Ferroptosis: molecular mechanisms and health implications. Cell Res. 2021;31:107–25.
Wang S, Xia P, Zhang L, Yu L, Liu H, Meng Q, et al. Systematical identification of breast cancer-related circular RNA modules for deciphering circRNA functions based on the non-negative matrix factorization algorithm. Int J Mol Sci. 2019;20:919.
Ceylan H. Identification of hub genes associated with obesity-induced hepatocellular carcinoma risk based on integrated bioinformatics analysis. Med Oncol. 2021;38:63.
Ma L, Lu H, Chen R, Wu M, Jin Y, Zhang J, et al. Identification of key genes and potential new biomarkers for ovarian aging: a study based on RNA-Sequencing data. Front Genet. 2020;11:590660.
Ceylan H. Integrated bioinformatics analysis to identify alternative therapeutic targets for Alzheimer’s disease: insights from a synaptic machinery perspective. J Mol Neurosci. 2022;72:273–86.
Koudstaal T, Boomars KA, Kool M. Pulmonary Arterial Hypertension and Chronic Thromboembolic Pulmonary Hypertension: an immunological perspective. J Clin Med. 2020;9:561.
Halliday SJ, Matthews DT, Talati MH, Austin ED, Su YR, Absi TS, et al. A multifaceted investigation into molecular associations of chronic thromboembolic pulmonary hypertension pathogenesis. JRSM Cardiovasc Dis. 2020;9:405416946.
Zhou N, Bao J. FerrDb: a manually curated resource for regulators and markers of ferroptosis and ferroptosis-disease associations. Database (Oxford). 2019;10:1523.
Dennis GJ, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, et al. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4(5):3.
Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10:1523.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. PRotein Sci. 2019;28:1947–51.
Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51:587–92.
Luo X, Xiang T, Huang H, Ye L, Huang Y, Wu Y. Identification of significant immune-related genes for epilepsy via bioinformatics analysis. Ann Transl Med. 2021;9:1161.
Liang Y, Lin F, Huang Y. Identification of biomarkers associated with diagnosis of osteoarthritis patients based on bioinformatics and machine learning. J Immunol Res. 2022;2022:5600190.
Xu W, Deng M, Meng X, Sun X, Tao X, Wang D, et al. The alterations in molecular markers and signaling pathways in chronic thromboembolic pulmonary hypertension, a study with transcriptome sequencing and bioinformatic analysis. Front Cardiovasc Med. 2022;9:961305.
Kawada JI, Takeuchi S, Imai H, Okumura T, Horiba K, Suzuki T, et al. Immune cell infiltration landscapes in pediatric acute myocarditis analyzed by Cibersort. J Cardiol. 2021;77:174–8.
Yang J, Fan Y, Liu S. ATF3 as a potential diagnostic marker of early-stage osteoarthritis and its correlation with immune infiltration through bioinformatics analysis. Bone Joint Res. 2022;11:679–89.
Xi Q, Liu Z, Song Y, Gan H, Huang Z, Luo Q, et al. Proteomic analyses of endarterectomized tissues from patients with chronic thromboembolic pulmonary hypertension. Cardiology. 2020;145:48–52.
Humbert M, Farber HW, Ghofrani HA, Benza RL, Busse D, Meier C, et al. Risk assessment in pulmonary arterial hypertension and chronic thromboembolic pulmonary hypertension. Eur Respir J. 2019;53:1802004.
Deng C, Wu D, Yang M, Chen Y, Ding H, Zhong Z, et al. The role of tissue factor and autophagy in pulmonary vascular remodeling in a rat model for chronic thromboembolic pulmonary hypertension. Respir Res. 2016;17:65.
Liu YG, Teng YS, Shan ZG, Cheng P, Hao CJ, Lv YP, et al. Arrestin domain containing 3 promotes Helicobacter pylori-associated gastritis by regulating protease-activated receptor 1. Jci Insight. 2020;5:e135849.
Chen C, Shen H, Huang Q, Li Q. The Circular RNA CDR1as Regulates the Proliferation and Apoptosis of Human Cardiomyocytes Through the miR-135a/HMOX1 and miR-135b/HMOX1 Axes. Genet Test Mol Biomarkers. 2020;24:537–48.
Zhang X, Guan T, Yang B, Chi Z, Wan Q, Gu HF. SLC30A7 has anti-oxidant stress effects in high glucose-induced apoptosis via the NFE2L2/HMOX1 signal transduction pathway. Diabetes Res Clin Pract. 2021;172:108445.
Hao K, Jiang W, Zhou M, Li H, Chen Y, Jiang F, et al. Targeting BRD4 prevents acute gouty arthritis by regulating pyroptosis. Int J Biol Sci. 2020;16:3163–73.
Mu J, Zhang D, Tian Y, Xie Z, Zou MH. BRD4 inhibition by JQ1 prevents high-fat diet-induced diabetic cardiomyopathy by activating PINK1/Parkin-mediated mitophagy in vivo. J Mol Cell Cardiol. 2020;149:1–14.
Yan L, Li X, Liu Z, Zhao Z, Luo Q, Zhao Q, et al. Research progress on the pathogenesis of CTEPH. Heart Fail Rev. 2019;24:1031–40.
Miao R, Dong X, Gong J, Wang Y, Guo X, Li Y, et al. Possible immune regulation mechanisms for the progression of chronic thromboembolic pulmonary hypertension. Thromb Res. 2021;198:122–31.
Quarck R, Wynants M, Verbeken E, Meyns B, Delcroix M. Contribution of inflammation and impaired angiogenesis to the pathobiology of chronic thromboembolic pulmonary hypertension. Eur Respir J. 2015;46:431–43.
Wu Y, Chen H, Xuan N, Zhou L, Wu Y, Zhu C, et al. Induction of ferroptosis-like cell death of eosinophils exerts synergistic effects with glucocorticoids in allergic airway inflammation. Thorax. 2020;75:918–27.
Zhang H, Liu J, Zhou Y, Qu M, Wang Y, Guo K, et al. Neutrophil extracellular traps mediate m(6)A modification and regulates sepsis-associated acute lung injury by activating ferroptosis in alveolar epithelial cells. Int J Biol Sci. 2022;18:3337–57.
Acknowledgements
We thank professor Halliday SJ for generous sharing GSE130391 dataset from GEO.
Funding
This work was supported by the research grant from Guangzhou Medical University Undergraduate Science and Technology Innovation Project (2022A043).
Author information
Authors and Affiliations
Contributions
LJP was responsible for the design. LSF and ZYZ performed execution and data curation of the project. LWH carried through supervision and processed writing-original draft. All authors commented on and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study analyzed existing data sets and did not require ethical review.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1:
Supplementary Material 1. Data of GSE130391.
Additional file 2:
Supplementary Material 2. 676 DEGs between CTEPH and Control group.
Additional file 3:
Supplementary Material 3. 259 Ferroptosis-related genes.
Additional file 4:
Supplementary Material 4. The expression of ARRDC3, HMOX1, BRD4 and YWHAE in GSE188938.
Additional file 5:
Supplementary Material 5. The expression of 22 immune cells in the sample and Correlation analysis.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Lin, J., Lin, S., Zhang, Y. et al. Identification of Ferroptosis-related potential biomarkers and immunocyte characteristics in Chronic Thromboembolic Pulmonary Hypertension via bioinformatics analysis. BMC Cardiovasc Disord 23, 504 (2023). https://doi.org/10.1186/s12872-023-03511-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12872-023-03511-5