Identification of Signal Pathways and Hub Genes of Pulmonary Arterial Hypertension by Bioinformatic Analysis

Wei, Rui-Qi; Zhang, Wen-Mei; Liang, Zhe; Piao, Chunmei; Zhu, Guangfa

doi:https://doi.org/10.1155/2022/1394088

Canadian Respiratory Journal

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Disclosure Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 1394088 | https://doi.org/10.1155/2022/1394088

Identification of Signal Pathways and Hub Genes of Pulmonary Arterial Hypertension by Bioinformatic Analysis

Rui-Qi Wei,¹Wen-Mei Zhang,²Zhe Liang,³Chunmei Piao,⁴and Guangfa Zhu²

Academic Editor: Alexandru Corlateanu

Received26 May 2022

Accepted19 Jul 2022

Published29 Aug 2022

Abstract

Pulmonary arterial hypertension (PAH) is a progressive and complex pulmonary vascular disease with poor prognosis. The aim of this study was to provide a new understanding of the pathogenesis of disease and potential treatment targets for patients with PAH based on multiple-microarray analysis.Two microarray datasets (GSE53408 and GSE113439) downloaded from the Gene Expression Omnibus (GEO) database were analysed. All the raw data were processed by R, and differentially expressed genes (DEGs) were screened out by the “limma” package. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed and visualized by R and Cytoscape software. Protein-protein interactions (PPI) of DEGs were analysed based on the NetworkAnalyst online tool. A total of 442 upregulated DEGs and 84 downregulated DEGs were identified. GO enrichment analysis showed that these DEGs were mainly enriched in mitotic nuclear division, organelle fission, chromosome segregation, nuclear division, and sister chromatid segregation. Significant KEGG pathway enrichment included ribosome biogenesis in eukaryotes, RNA transport, proteoglycans in cancer, dilated cardiomyopathy, rheumatoid arthritis, vascular smooth muscle contraction, focal adhesion, regulation of the actin cytoskeleton, and hypertrophic cardiomyopathy. The PPI network identified 10 hub genes including HSP90AA1, CDC5L, MDM2, LRRK2, CFTR, IQGAP1, CAND1, TOP2A, DDX21, and HIF1A. We elucidated potential biomarkers and therapeutic targets for PAH by bioinformatic analysis, which provides a theoretical basis for future study.

1. Introduction

According to the latest clinical guidelines, pulmonary hypertension (PH) is divided into five major categories, the first of which is pulmonary arterial hypertension (PAH). PAH is a progressive and fatal disorder accompanied by nonspecific clinical manifestations such as chest pain, dyspnea, and cyanosis [1]. The main pathophysiological changes in PAH are increases in pulmonary arterial pressure and pulmonary vascular resistance (PVR) because of vascular proliferation and adverse remodeling in the distal pulmonary arteriole, giving rise to right heart failure and even death [1, 2]. Although basic discoveries and pivotal clinical trials have led to the development of medications in recent years [3], patients with PAH still have a poor prognosis. According to a previous study of Chinese people, the 1- and 3-year survival rates are only 68% and 39%, respectively, and even though survival rates for PAH have improved in China in the targeted therapy era, the mortality is still high [4]. The existing medications, including prostacyclin analogues and receptor agonists, phosphodiesterase type 5 inhibitors, and endothelin receptor antagonists, are aimed at correcting the imbalance of vasoactive factors in PAH [5, 6]. There is a shortage of new drugs to address other pathological mechanisms such as vascular remodeling [6]. Therefore, it is important to elucidate the molecular mechanisms and to identify novel drug targets for potential therapy.

Microarrays are powerful and effective tools to detect the gene expression differences between health and disease conditions. Microarray technology has been used for more than two decades and is a well-established and mature biochemical technique [7]. A large amount of microarray data is freely available from the NCBI-Gene Expression Omnibus (GEO), which is an open access database. The objective of this study was to use transcriptomic microarray data analysis to find genes associated with the onset of pulmonary hypertension, and to provide a research basis for identifying new potential therapeutic targets. In this study, we searched for and screened the transcriptome chip data of PAH from GEO to find differentially expressed genes (DEGs) between PAH and controls using bioinformatics analyses. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed to detect statistically significant signaling pathways of DEGs. Also, a protein-protein interaction (PPI) network was established to identify hub genes. This research will provide a direction for the discovery of diagnostic biomarkers and targeted therapy strategies related to PAH.

2. Materials and Methods

2.1. Microarray Data Source and Screening

The GEO is a public database that stores microarray, next-generation sequencing, and other high-throughput sequencing data (https://www.ncbi.nlm.nih.gov/). Using the search terms “pulmonary hypertension” and “microarray,” and after filtering by “Homo sapiens” and “Expression profiling by array,” we identified 15 transcriptomic studies. Then, only human transcriptomic studies from lung tissue and with the same GEO platform were included. We ultimately selected two GEO series (GSE) studies (GSE53408, GSE113439) for further analysis. The study information is shown in Table 1.

In the two original studies, lung tissue samples were from PAH patients undergoing lung transplantation and healthy lung tissue was derived from normal tissue adjacent to cancer during lung cancer resection. First, total RNA was extracted by the Trizol method and reverse transcribed into cDNA. Second, cRNA was obtained by in vitro transcription using cDNA as a template and labeled by biotin. Third, GeneChips were hybridized with cRNA, washed, and stained in an Affymetrix Fluidics Station 400. Finally, GeneChips were scanned using a GeneArray scanner. The two GSE studies used two different scanners.

2.2. Data Processing and Quality Control

GSE53408 [8] contains 11 normal and 12 severe PAH lung tissue samples, while GSE113439 [9] contains 11 normal and 15 PAH samples (one CTEPAH sample (GSM3106338) was removed because it does not belong to WHO Group 1 pulmonary hypertension). We downloaded the raw CEL files of GSE53408 and GSE113439 from the GEO database. R software (version 4.0.2) was used to process the raw data. The “affy” package [10] was used to read the CEL files and extract Affymetrix GeneChip probe level data. After correcting background, normalizing, and summarizing by using the Robust Multichip Average (RMA) algorithm [11], we got the probe expression matrix. Finally, we converted probe IDs into gene symbols (international standard names of genes) using the GPL6244 annotation package (“hugene10sttranscriptcluster.db” [12]). Probe IDs without corresponding gene symbols were deleted, and the maximum value was regarded as the gene expression value when multiple probes corresponded to one gene. We further carried out quality control for the gene expression matrix obtained above using cluster analysis and principal component analysis (PCA).

2.3. Screening and Integration of DEGs

We used the “limma” package [13] for DEG analysis. Gene expression values have been log2 transformed and normalized after RMA algorithm processing and were fitted to a linear model by using weighted least squares. Differential expression analysis was performed between PAH and controls. “Limma” package calculated the log2 fold change (log2FC) and false discovery rate (FDR) adjusted value for each gene. FDR uses the Benjamini–Hochberg method to correct the value in order to control the number of false positives for multiple tests. Adjusted value <0.05 and |log2FC| ≥ 1.0 were considered as DEGs. Considering the large batch effect between GSE53408 and GSE113439, we analysed and obtained the DEGs for each dataset. Ultimate DEGs were integrated by using the Venn diagram tool (https://bioinformatics.psb.ugent.be/webtools/Venn/).

2.4. GO and KEGG Pathway Enrichment Analysis

GO enrichment analysis of DEGs was performed using the “enrichGO” function of the “clusterProfiler” package [14]. Genes were divided into three categories in function including biological process (BP), molecular function (MF), and cellular components (CC). GO terms were considered significantly enriched for the DEGs while an adjusted value <0.05. The Benjamini and Hochberg method was used to adjust the raw values. The KEGG (https://www.genome.jp/kegg/) pathway enrichment analysis of integrated DEGs in our study was performed using the DAVID webtools (https://david.ncifcrf.gov/). value <0.05 was considered to be significant statistically. The visualization of analytical results was achieved by using the R (4.0.2 version) and Cytoscape software (3.8.0 version).

2.5. Tissue-Specific PPI Network Analysis and Hub Gene Identification

Different gene lists were inputted into NetworkAnalyst for lung tissue-specific PPI network analyses [15]. NetworkAnalyst is a bioinformatic online analysis tool with power functions and friendly operation interface, which can be accessed in https://www.networkanalyst.ca/. PPI data from human tissues were obtained from the DifferentialNet database (https://netbio.bgu.ac.il/diffnet/) [16]. Hub genes were further identified by a selected minimum connected network. Each node represents a protein encoded by a gene while connections between nodes represent the interaction of proteins. The nodes of maximum connectivity were considered as core proteins or hub genes with important biological regulatory functions in PAH.

3. Results

3.1. Microarray Data Information

To identify DEGs in PAH compared with healthy individuals, two gene expression profiles were selected in our study. They were GSE53408 and GSE113439. Raw data were downloaded and processed by the “affy” package. The gene expression values were normalized by the RMA integrated preprocessing algorithm. Boxplots of data before and after normalization are shown in Figure 1(a). The results of cluster analysis and PCA indicate that the PAH group and the control group data can be readily distinguished, as shown in Figure 1(b). Therefore, subsequent genic difference analysis could be carried out.

(a)

(b)

(c)

(d)

Figure 1

Normalization of gene expression and cluster analysis and PCA analysis for each sample. (a) Normalization of gene expression for GSE53408 data. (b) Normalization of gene expression for GSE113439 data. Abscissa shows sample lists, and ordinates show the gene expression value (log2 transformed). Blue bars represent the data before normalization, and red bars represent the normalized data. (c) Cluster analysis and PCA analysis for each sample of GSE53408 data. (d) Cluster analysis and PCA analysis for each sample of GSE113439 data. The PAH group data and the normal group data are readily distinguished. Abbreviations: IPAH, idiopathic pulmonary arterial hypertension; CTDPAH, connective tissue disease-associated PAH; CHDPAH, congenital heart disease-associated PAH.

3.2. DEG Analysis in PAH

DEGs between PAH and normal lung tissues were evaluated by the “limma” package. There were 635 DEGs in GSE53408, including 525 upregulated genes and 110 downregulated genes (Figure 2(a)). And 564 DEGs were screened in GSE113439, including 468 upregulated genes and 96 downregulated genes (Figure 2(a)). And the top 100 DEGs are displayed in the form of heatmaps, as shown in Figure 2(b). Taking into account the original experimental processes such as different model Genechip scanners, the two sets of data have a large batch effect, so a Venn diagram was constructed to obtain the final DEGs, as shown in Figure 2(c). Finally, 526 DEGs were identified in total, with 442 upregulated genes and 84 downregulated genes (Table 2). The top twenty upregulated and downregulated DEGs sorted by average log2FC size of the two datasets are shown in Table S1.

(a)

(b)

(c)

Figure 2

Volcano plots, cluster heatmaps, and Venn diagram of DEGs for GSE53408 and GSE113439 datasets. (a) Volcano plots of GSE53408 and GSE113439 datasets. Red dots represent upregulated genes based on adjusted value <0.05 and log2FC ≥ 1. Blue dots represent downregulated genes based on adjusted value<0.05 and log2FC ≤ −1. The black points represent genes with no significant difference. FC is the fold change. (b) Cluster heatmaps of the top 100 significant genes of GSE53408 and GSE113439 data. Red indicates that the gene expression is upregulated, blue indicates that the gene expression is downregulated, and white indicates no significant changes in gene expression. Abbreviations: DEGs, differentially expressed genes. (c) Venn diagram of DEGs, left panel: the upregulated overlapping DEGs in the two datasets are 442; right panel: the downregulated overlapping DEGs in the two datasets are 84.

3.3. GO Term Enrichment Analysis of DEGs

GO enrichment analysis for integrated DEGs was performed using the “clusterProfiler” package. The enriched GO terms were divided into three ontologies including biological process (BP), molecular function (MF), and cell component (CC). Significant results of GO enrichment analysis are exhibited in Figure 3(a). In BP, DEGs were generally enriched in mitotic nuclear division, organelle fission, chromosome segregation, nuclear division, and sister chromatid segregation. In MF, DEGs were chiefly enriched in ATPase activity, helicase activity, DNA-dependent ATPase activity, catalytic activity acting on DNA, and DNA helicase activity. In CC, DEGs were mainly enriched in condensed chromosome, spindle, chromosomal region, mitotic spindle, and microtubule (Table S2 and Figure S1).

(a)

(b)

Figure 3

Barplots of GO enrichment analysis (top 20) and Dotplot of KEGG pathway enrichment analysis. (a) Y-axes show GO enrichment significance items of DEGs in three different functional groups (molecular function, biological processes, and cell composition), and x-axes represent the gene counts enriched in each GO term. Different bar colors represent different adjust values. Abbreviations: DEGs, differentially expressed genes; GO, gene ontology. (b) Dotplot of KEGG pathway enrichment analysis. Each dot represents a KEGG pathway, and the size of each dot indicates the number of genes in the enrichment pathway.

3.4. KEGG Pathway Enrichment Analysis of DEGs

After inputting the DEGs list into DAVID online tool, the most significantly enriched KEGG pathways of the DEGs were obtained (Table S3). The signaling pathways of DEGs were eventually enriched in ribosome biogenesis in eukaryotes, RNA transport, proteoglycans in cancer, dilated cardiomyopathy, rheumatoid arthritis, vascular smooth muscle contraction, focal adhesion, regulation of actin cytoskeleton, and hypertrophic cardiomyopathy (Figure 3(b)). Cytoscape software was used to calculate the network topological characteristics and draw each node and edge. Different pathway terms and genes are displayed by ovals in different colors, as shown in Figure 4(a).

(a)

(b)

3.5. Tissue-Specific PPI Network Analysis and Hub Gene Identification

Protein-protein interactions among the DEGs were predicted by the NetworkAnalyst tool. Human lung tissue-specific networks were constructed based on the DifferentialNet database. We chose a minimum connected network for identifying hub genes in PAH. This PPI network contains 1030 nodes and 3364 edges (Figure 4(b)). The top ten hub genes evaluated according to the degree of connectivity in the PPI network were HSP90AA1, CDC5L, MDM2, LRRK2, CFTR, IQGAP1, CAND1, TOP2A, DDX21, and HIF1A (Table 3). All of these identified hub genes were upregulated in PAH patients.

4. Discussion

Pulmonary arterial hypertension is a chronic and progressive disease with few typical manifestations and thus could be delayed in diagnosis and treatment [2]. Unfortunately, the fundamental pathogenesis of PAH preliminarily remains unclear warranting the necessity to explore basic mechanisms to delineate specific targets promising for therapeutic intervention.

In this study, we downloaded two microarray datasets, GSE53408 and GSE113439, from the GEO database and separately obtained the DEGs by bioinformatics analysis. Then the upregulated and downregulated genes from two datasets were intersected by a Venn diagram webtool to obtain the final integrated DEGs. We identified 526 DEGs, including 442 upregulated genes and 84 downregulated genes. The integrated DEGs were further subjected to GO and KEGG pathway enrichment analysis and tissue-specific PPI network analysis. There has been a number of PAH studies using microarray analyses [17–19], and the signaling pathway or hub genes identified by each study varied widely. This was mainly due to different microarray platforms or different data processing procedures. In Li’s study [18], they analysed the DEGs using the GEO2R online tool. Although the GEO2R tool is simple and convenient, the average gene expression value varies from one sample to another in the raw data, so it is hard to ensure the accuracy of the analysis results. Our study started by analyzing raw CEL format microarray data and normalized it using the RMA algorithm to keep the gene expression levels of each sample at the same level, which ensured the reliability of downstream analysis. Luo et al. [17] combined the disease group and control samples of the two datasets and used the “sva” R package to remove batch effects for DEG analysis that may cause some deviation considering the use of different GeneArray scanners between GSE53408 and GSE113439. In this study, we analysed GSE53408 and GSE113439 and merged the DEGs of the two datasets, which may be more reasonable. In addition, one CTEPH sample in GSE113439 was not selected in our study because we only paid attention to the WHO group 1 PH.

The DEGs associated with GO BP terms mainly consisted of mitotic nuclear division, organelle fission, chromosome segregation, nuclear division, and sister chromatid segregation, which all participated in the process of cell division and proliferation. In PAH progression, excessive proliferation of pulmonary artery smooth muscle cells plays an important role, which causes vascular medial thickening and pulmonary vascular remodeling [20]. DEGs associated with GO MF terms such as ATPase activity, helicase activity, DNA-dependent ATPase activity, catalytic activity acting on DNA, and DNA helicase activity indicated that DNA replication was active, which was consistent with the results of BP. The enriched KEGG pathways of DEGs such as vascular smooth muscle contraction, focal adhesion, and regulation of the actin cytoskeleton may be associated with the progression of PAH. As we all know, pulmonary vascular smooth muscle contraction is a critical pathological mechanism of PAH, especially in the state induced by hypoxia. Hypoxia promoted pulmonary vasoconstriction by interacting with a variety of protein kinases, calcium channels, and potassium channels [21]. Blocking these signal transduction pathways may provide novel therapeutic approaches for PAH depending on further research. The focal adhesion pathway was also identified in PAH by Zhao and Austin’s study [22]. Focal adhesion dynamics is regulated by focal adhesion kinase (FAK) family signaling via integrins [23]. Several studies have demonstrated the association of FAK with vascular diseases including PAH [23, 24]. Jia et al. [24] found that activation of FAK could stimulate pulmonary arterial smooth muscle cell proliferation in vitro. Paulin et al. [25] confirmed that FAK inhibitor improved hemodynamics, vascular remodeling, and right ventricular hypertrophy in monocrotaline induced PAH rat model. Based on the above research, FAK inhibition may open a new therapeutic way for PAH. Finally, regulation of the actin cytoskeleton was found to be involved in the process of smooth muscle cell contraction, proliferation, and migration [26]. Abelson tyrosine kinase (c-Abl) could modulate restructuring of actin cytoskeletal [27]. There was some research [28, 29] suggesting that c-Abl inhibitor imatinib attenuated the symptoms of patients with PAH.

Network analysis is a powerful tool that helps to identify hub genes and to discover novel targets for therapy [30]. We constructed a PPI network of proteins encoded by DEGs and screened 10 hub genes according to their connectivity degree: HSP90AA1, CDC5L, MDM2, LRRK2, CFTR, IQGAP1, CAND1, TOP2A, DDX21, and HIF1A. We searched for the relationship between each gene symbol and pulmonary hypertension in PubMed and found that four genes involved in PAH had been studied by previous researchers. Shen et al. [31] reported that MDM2-mediated ubiquitination of ACE2 contributed to the development of PAH, and the MDM2 inhibitor JNJ-165 reversed SU5416/hypoxia-induced PAH in mice model. CFTR is expressed in PASMCs and modulates pulmonary arterial tone [32]. In hypoxia-induced PAH mice, deficiency of CFTR signally attenuated right ventricular systolic pressure, pulmonary vessel medial wall thickness, and muscularization [33]. Grigolo et al. [34] found that antitopoisomerase II (TOP2A) autoantibodies were related to systemic sclerosis associated with pulmonary hypertension. Luo et al. [35] reported that the CD146-HIF-1α axis in PASMC derived vascular remodeling and PAH. Modulation of CD146-HIF-1α may be a new antiremodeling therapy strategy. The remaining six hub genes, HSP90AA1, CDC5L, LRRK2, IQGAP1, CAND1, and DDX21, were not found to have a direct connection with PAH previously. HSP90AA1, as a molecular chaperone, assisted the assembly of intracellular molecules and protein folding, and participated in multiple processes of fibroblast activation [36, 37]. Upregulation of CDC5L could promote cell cycle progression and proliferation [38]. Hongge et al. [39] found that LRRK2 was expressed in endothelial cells and regulated the procedure of monocyte adhesion, which was involved in vascular inflammation and immune response. IQGAP1, as a scaffold protein, maintained the integrity of the endothelium barrier [40]. Gharib et al. [41] discovered that IQGAP1 was a novel candidate biomarker and played a crucial role in hypoxia-induced pulmonary hypertension and vascular remodeling by global gene annotation analysis.

There are several limitations in our study. At first, microarray technology can only detect known genes, and due to discrepancies in chip probes designed by different companies, the conclusions obtained above might vary. Then, considering the relatively small sample size in the two public datasets (48 samples in total), there may be some bias on screened DEGs. Finally, all of our results were based on the forecasting of bioinformatic analysis, which should be further confirmed by in vitro and in vivo experiments.

5. Conclusion

In summary, we analysed the DEGs between patients and healthy individuals aiming to assist the understanding of the pathogenesis of PAH. Vascular smooth muscle contraction, focal adhesion, and regulation of the actin cytoskeleton pathways were found to be involved in PAH in this study. And ten hub genes, including HSP90AA1, CDC5L, MDM2, LRRK2, CFTR, IQGAP1, CAND1, TOP2A, DDX21, and HIF1A may be new therapeutic targets for patients with PAH to be further validated.

Data Availability

All data relevant to the study are included in the article and are available from the corresponding author on reasonable request.

Disclosure

Chunmei Piao and Guangfa Zhu are the co-first authors.

Conflicts of Interest

The authors have no conflicts of interest to declare.

Authors’ Contributions

Chunmei Piao and Guangfa Zhu contributed equally to this work. Rui-Qi Wei, Chunmei Piao, and Guangfa Zhu contributed to the study design, analysis, and writing of the article. Wen-Mei Zhang and Zhe Liang contributed to data analysis. All authors read and approved the manuscript.

Acknowledgments

The authors thank the contributors of the GEO (https://www.ncbi.nlm.nih.gov/geo/) database for sharing their data via open access. The authors also thank Accdon for language editing and retouching services. This work was supported by the National Natural Science Foundation of China (Nos. 82170409, 81800222, 81770468, and 81970067).

Supplementary Materials

Figure S1: genes enriched in each GO term. (a) Biological process of DEGs; (b) molecular function of DEGs; (c) cell component of DEGs. Brown dots represent GO terms, and black dots represent genes. Table S1: top 20 upregulated and downregulated DEGs according to average gene expression fold change from two datasets. Table S2: GO enriched analysis of DEGs between PAH and normal tissues (top five terms in BP, MF, and CC of DEGs according to value). Table S3: KEGG enriched analysis of DEGs between PAH and normal tissues. (Supplementary Materials)

References

M. Humbert, C. Guignabert, S. Bonnet et al., “Pathology and pathobiology of pulmonary hypertension: state of the art and research perspectives,” European Respiratory Journal, vol. 53, no. 1, 2019.
View at: Google Scholar
D. Kim and M. P. George, “Pulmonary hypertension,” Medical Clinics of North America, vol. 103, no. 3, pp. 413–423, 2019.
View at: Google Scholar
J. R. Klinger, C. G. Elliott, D. J. Levine et al., “Therapy for pulmonary arterial hypertension in adults: update of the chest guideline and expert panel report,” Chest, vol. 155, no. 3, pp. 565–586, 2019.
View at: Publisher Site | Google Scholar
R. Zhang, L. Z. Dai, W. P. Xie et al., “Survival of Chinese patients with pulmonary arterial hypertension in the modern treatment era,” Chest, vol. 140, no. 2, pp. 301–309, 2011.
View at: Publisher Site | Google Scholar
N. Galie, M. Humbert, J. L. Vachiery et al., “Esc/ers guidelines for the diagnosis and treatment of pulmonary hypertension: the joint task force for the diagnosis and treatment of pulmonary hypertension of the European society of cardiology (esc) and the European respiratory society (ers): endorsed by: association for European paediatric and congenital cardiology (aepc), international society for heart and lung transplantation (ishlt),” European Heart Journal, vol. 37, no. 1, pp. 67–119, 2015.
View at: Publisher Site | Google Scholar
O. Sitbon, M. Gomberg-Maitland, J. Granton et al., “Clinical trial design and new therapies for pulmonary arterial hypertension,” European Respiratory Journal, vol. 53, no. 1, p. 1801908, 2019.
View at: Publisher Site | Google Scholar
J. Hoffmann, J. Wilhelm, A. Olschewski, and G. Kwapiszewska, “Microarray analysis in pulmonary hypertension,” European Respiratory Journal, vol. 48, no. 1, pp. 229–241, 2016.
View at: Publisher Site | Google Scholar
Y. D. Zhao, H. Z. H. Yun, J. Peng et al., “De novo synthesize of bile acids in pulmonary arterial hypertension lung,” Metabolomics, vol. 10, no. 6, pp. 1169–1175, 2014.
View at: Publisher Site | Google Scholar
M. Mura, M. J. Cecchini, M. Joseph, and J. T. Granton, “Osteopontin lung gene expression is a marker of disease severity in pulmonary arterial hypertension,” Respirology, vol. 24, no. 11, pp. 1104–1110, 2019.
View at: Publisher Site | Google Scholar
L. Gautier, L. Cope, B. M. Bolstad, and R. A. Irizarry, “Affy--Analysis of Affymetrix Genechip data at the probe level,” Bioinformatics, vol. 20, no. 3, pp. 307–315, 2004.
View at: Publisher Site | Google Scholar
R. A. Irizarry, B. Hobbs, F. Collin et al., “Exploration, normalization, and summaries of high density oligonucleotide array probe level data,” Biostatistics, vol. 4, no. 2, pp. 249–264, 2003.
View at: Publisher Site | Google Scholar
J. W. MacDonald, “Hugene10sttranscriptcluster.Db: Affymetrix Hugene10 annotation data (chip Hugene10sttranscriptcluster),” 2017, https://bioconductor.org/packages/release/data/annotation/html/hugene10sttranscriptcluster.db.html.
View at: Google Scholar
M. E. Ritchie, B. Phipson, D. Wu et al., “Limma powers differential expression analyses for rna-sequencing and microarray studies,” Nucleic Acids Research, vol. 43, no. 7, p. e47, 2015.
View at: Publisher Site | Google Scholar
G. Yu, L. G. Wang, Y. Han, and Q. Y. He, “Clusterprofiler: an R package for comparing biological themes among gene clusters,” OMICS: A Journal of Integrative Biology, vol. 16, no. 5, pp. 284–287, 2012.
View at: Publisher Site | Google Scholar
J. Xia, M. J. Benner, and R. E. Hancock, “Networkanalyst—Integrative approaches for protein-protein interaction network analysis and visual exploration,” Nucleic Acids Research, vol. 42, pp. W167–W174, 2014.
View at: Google Scholar
O. Basha, R. Shpringer, C. M. Argov, and E. Yeger-Lotem, “The differentialnet database of differential protein-protein interactions in human tissues,” Nucleic Acids Research, vol. 46, no. D1, pp. D522–D526, 2018.
View at: Publisher Site | Google Scholar
J. Luo, H. Li, Z. Liu et al., “Integrative analyses of gene expression profile reveal potential crucial roles of mitotic cell cycle and microtubule cytoskeleton in pulmonary artery hypertension,” BMC Medical Genomics, vol. 13, no. 1, p. 86, 2020.
View at: Publisher Site | Google Scholar
Q. Li, L. Meng, and D. Liu, “Screening and identification of therapeutic targets for pulmonary arterial hypertension through microarray technology,” Frontiers in Genetics, vol. 11, p. 782, 2020.
View at: Publisher Site | Google Scholar
E. Zhao, H. Xie, and Y. Zhang, “Identification of differentially expressed genes associated with idiopathic pulmonary arterial hypertension by integrated bioinformatics approaches,” Journal of Computational Biology, vol. 28, no. 1, pp. 79–88, 2021.
View at: Publisher Site | Google Scholar
C. Morii, H. Y. Tanaka, Y. Izushi et al., “3d in vitro model of vascular medial thickening in pulmonary arterial hypertension,” Frontiers in Bioengineering and Biotechnology, vol. 8, p. 482, 2020.
View at: Publisher Site | Google Scholar
N. Sommer, I. Strielkov, O. Pak, and N. Weissmann, “Oxygen sensing and signal transduction in hypoxic pulmonary vasoconstriction,” European Respiratory Journal, vol. 47, no. 1, pp. 288–303, 2016.
View at: Publisher Site | Google Scholar
M. Zhao, E. D. Austin, A. R. Hemnes, J. E. Loyd, and Z. Zhao, “An evidence-based knowledgebase of pulmonary arterial hypertension to identify genes and pathways relevant to pathogenesis,” Molecular BioSystems, vol. 10, no. 4, pp. 732–740, 2014.
View at: Publisher Site | Google Scholar
J. M. Murphy, K. Jeong, and S. T. S. Lim, “Fak family kinases in vascular diseases,” International Journal of Molecular Sciences, vol. 21, no. 10, p. 3630, 2020.
View at: Publisher Site | Google Scholar
D. Jia, Q. Zhu, H. Liu et al., “Osteoprotegerin disruption attenuates HySu-induced pulmonary hypertension through integrin α v β₃/FAK/AKT pathway suppression,” Circulation: Cardiovascular Genetics, vol. 10, no. 1, Article ID e001591, 2017.
View at: Publisher Site | Google Scholar
R. Paulin, J. Meloche, A. Courboulin et al., “Targeting cell motility in pulmonary arterial hypertension,” European Respiratory Journal, vol. 43, no. 2, pp. 531–544, 2014.
View at: Publisher Site | Google Scholar
D. D. Tang, “The dynamic actin cytoskeleton in smooth muscle,” Advances in Pharmacology, vol. 81, pp. 1–38, 2018.
View at: Google Scholar
R. A. Cleary, R. Wang, T. Wang, and D. D. Tang, “Role of abl in airway hyperresponsiveness and airway remodeling,” Respiratory Research, vol. 14, no. 1, p. 105, 2013.
View at: Publisher Site | Google Scholar
K. C. Patterson, A. Weissmann, T. Ahmadi, and H. W. Farber, “Imatinib mesylate in the treatment of refractory idiopathic pulmonary arterial hypertension,” Annals of Internal Medicine, vol. 145, no. 2, pp. 152-153, 2006.
View at: Publisher Site | Google Scholar
H. A. Ghofrani, W. Seeger, and F. Grimminger, “Imatinib for the treatment of pulmonary arterial hypertension,” New England Journal of Medicine, vol. 353, no. 13, pp. 1412-1413, 2005.
View at: Publisher Site | Google Scholar
K. Ozturk, M. Dow, D. E. Carlin, R. Bejar, and H. Carter, “The emerging potential for network analysis to inform precision cancer medicine,” Journal of Molecular Biology, vol. 430, no. 18, pp. 2875–2899, 2018.
View at: Publisher Site | Google Scholar
H. Shen, J. Zhang, C. Wang et al., “Mdm2-Mediated ubiquitination of angiotensin-converting enzyme 2 contributes to the development of pulmonary arterial hypertension,” Circulation, vol. 142, no. 12, pp. 1190–1204, 2020.
View at: Publisher Site | Google Scholar
R. Robert, J. P. Savineau, C. Norez, F. Becq, and C. Guibert, “Expression and function of cystic fibrosis transmembrane conductance regulator in rat intrapulmonary arteries,” European Respiratory Journal, vol. 30, no. 5, pp. 857–864, 2007.
View at: Publisher Site | Google Scholar
C. Tabeling, H. Yu, L. Wang et al., “Cftr and sphingolipids mediate hypoxic pulmonary vasoconstriction,” Proceedings of the National Academy of Sciences of the U S A, vol. 112, no. 13, pp. E1614–E1623, 2015.
View at: Publisher Site | Google Scholar
B. Grigolo, I. Mazzetti, R. Meliconi et al., “Anti-topoisomerase ii alpha autoantibodies in systemic sclerosis-association with pulmonary hypertension and hla-B35,” Clinical and Experimental Immunology, vol. 121, no. 3, pp. 539–543, 2001.
View at: Publisher Site | Google Scholar
Y. Luo, X. Teng, L. Zhang et al., “CD146-HIF-1α hypoxic reprogramming drives vascular remodeling and pulmonary arterial hypertension,” Nature Communications, vol. 10, no. 1, p. 3551, 2019.
View at: Publisher Site | Google Scholar
V. Sontake, Y. Wang, R. K. Kasam et al., “Hsp90 regulation of fibroblast activation in pulmonary fibrosis,” JCI Insight, vol. 2, no. 4, Article ID e91454, 2017.
View at: Publisher Site | Google Scholar
P. S. Bellaye, C. Shimbori, T. Yanagihara et al., “Synergistic role of HSP90α and HSP90β to promote myofibroblast persistence in lung fibrosis,” European Respiratory Journal, vol. 51, no. 2, Article ID 1700386, 2018.
View at: Publisher Site | Google Scholar
H. Liu, D. Guo, Y. Sha et al., “Anxa7 promotes the cell cycle, proliferation and cell adhesion-mediated drug resistance of multiple myeloma cells by up-regulating Cdc5l,” Aging (Albany NY), vol. 12, no. 11, pp. 11100–11115, 2020.
View at: Publisher Site | Google Scholar
L. Hongge, G. Kexin, M. Xiaojie, X. Nian, and H. Jinsha, “The role of Lrrk2 in the regulation of monocyte adhesion to endothelial cells,” Journal of Molecular Neuroscience, vol. 55, no. 1, pp. 233–239, 2015.
View at: Publisher Site | Google Scholar
D. Mehta, K. Ravindran, and W. M. Kuebler, “Novel regulators of endothelial barrier function,” American Journal of Physiology—Lung Cellular and Molecular Physiology, vol. 307, no. 12, pp. L924–L935, 2014.
View at: Publisher Site | Google Scholar
S. A. Gharib, D. L. Luchtel, D. K. Madtes, and R. W. Glenny, “Global gene annotation analysis and transcriptional profiling identify key biological modules in hypoxic pulmonary hypertension,” Physiological Genomics, vol. 22, no. 1, pp. 14–23, 2005.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Rui-Qi Wei et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

625

Downloads

805

Citations