WJCDWorld Journal of Cardiovascular Diseases2164-5329Scientific Research Publishing10.4236/wjcd.2024.143010WJCD-132067ArticlesMedicine&Healthcare Identification and Validation of Vascular-Associated Biomarkers for the Prognosis and Potential Pathogenesis of Hypertension Using Comprehensive Bioinformatics Methods XiangguangChang1LeiGuo2LiyingZou1YazhaoMa3JilinFeng4Department of Laboratory, Yangbi Yi Autonomous County People’s Hospital, Dali, ChinaDepartment of Infectious Diseases, Yangbi Yi Autonomous County People’s Hospital, Dali, ChinaMedical Community Office, Yangbi Yi Autonomous County People’s Hospital, Dali, ChinaDepartment of Obstetrics and Gynecology, Yangbi Yi Autonomous County People’s Hospital, Dali, China2603202414031151286, February 202424, March 2024 27, March 2024© Copyright 2014 by authors and Scientific Research Publishing Inc. 2014This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

Background: Hypertension, also known as increased blood pressure, is a phenomenon in which blood flows in blood vessels and causes persistently higher-than-normal pressure on the vessel wall. The identification of novel prognostic and pathogenesis biomarkers plays a key role in the management of hypertension. Methods: The GSE7483 and GSE75815 datasets from the gene expression omnibus (GEO) database were used to identify the genes associated with hypertension that were differentially expressed genes (DEGs). The functional role of the DEGs was elucidated by gene body (GO) enrichment analysis. In addition, we performed an immune infiltration assay and GSEA on the DEGs of hypertensive patients and verified the expression of novel DEGs in the blood of hypertensive patients by RT-qPCR. Results: A total of 267 DEGs were identified from the GEO database. GO analysis revealed that these genes were associated mainly with biological processes such as fibroblast proliferation, cell structural organization, extracellular matrix organization, vasculature development regulation, and angiogenesis. We identified five possible biomarkers, Ecm1, Sparc, Sphk1, Thbsl, and Mecp2, which correlate with vascular development and angiogenesis characteristic of hypertension by bioinformatics, and explored the clinical expression levels of these genes by RT-qPCR, and found that Sparc, Sphk1, and Thbs1 showed significant up-regulation, in agreement with the results of the bioinformatics analysis. Conclusion: Our study suggested that Sparc, Sphk1 and Thbs1 may be potential novel biomarkers for the diagnosis, treatment and prognosis of hypertension and that they are involved in the regulation of vascular development and angiogenesis in hypertension.

Hypertension Biomarkers Differentially Expressed Genes Vascular Development and Angiogenesis Bioinformatics Analysis1. Introduction

Hypertension is a major risk factor for cardiovascular diseases and end-stage kidney damage and currently affects approximately 40% of the global adult population [1] . Hypertension is a risk factor for the development of cardiovascular diseases such as coronary artery disease (CAD), left ventricular hypertrophy, valvular heart disease, atrial fibrillation, stroke and arrhythmias, without timely diagnosis and appropriate treatment, it may lead to disease or death [2] [3] . It is estimated that 1.4 billion people worldwide suffer from high blood pressure, but only 14% of patients have their condition under control [4] , and the prevention and control situation is very grim. Antihypertensive drugs are currently the main means of treating hypertension, but there are significant individual differences in the efficacy of these drugs, and there are many adverse reactions [5] . Existing studies have shown that differences in individualized medication used caused by genetic polymorphisms are the main reason for the low efficacy and greater side effects of current antihypertensive drugs [6] .

Bioinformatics is a multidisciplinary field that applies information technology and computer science principles to the study of the molecular aspects of biology and immunology [7] . Over the past decade, single-cell RNA sequencing (SCS) has emerged as a powerful tool for characterizing different cellular functions and states [8] [9] . In recent years, with the cross-fusion of computer science and medical fields, an increasing number of algorithms have been used in the diagnosis, prediction and prognosis of clinical diseases [10] . Weighted gene coexpression network analysis (WGCNA) has been widely used to identify central genes in clinically important modules by constructing gene coexpression networks [11] . Microarray data and high-throughput sequencing technology are widely used in gene sequencing for screening differentially expressed genes and have shown enormous potential in medical research because of their high diagnostic performance for various diseases [12] [13] . However, to date, a multibiomarker-based diagnostic model has not been developed for the diagnosis of hypertension.

In this study, we used bioinformatics methods to identify and validate genes associated with hypertension. This study explored the biological functions of genes associated with hypertension through the use of public databases and constructed a prediction model and novel biomarker interaction network. The findings can aid in the prediction and treatment of hypertension risk. The data of hypertension patients were downloaded from the GEO dataset and analyzed using R software to identify DEGs. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genomes (KEGG) pathway analysis methods were used to analyze the overlapping DEGs in the dataset. Using a variety of bioinformatics methods, novel biomarkers associated with the prognosis and potential pathogenesis of hypertension were screened.

2. Materials and Methodology2.1. Data Sources

The hypertension-related datasets GSE7483 and GSE75815 were downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/geo/). The GSE7483 dataset utilizes the GPL890 platform and comprises both treated and untreated Wistar Kyoto (WKY) rats, as well as treated and untreated spontaneously hypertensive rats (SHR). From this dataset, we specifically isolated 12 hypertensive samples (untreated SHR rats) and 12 control samples (untreated WKY rats). Regarding the GSE75815 dataset, analysis was performed utilizing the GPL6887 platform. This dataset investigates an animal model of AngII-induced hypertension utilizing Illumina WG-6 v2.0 microarrays. Specifically, the thoracic aorta, abdominal aorta, and mesenteric arteries of C57BL/6J mice infused with AngII (490 ng/kg/min) for 2 consecutive weeks were analyzed. From this dataset, our focus was on the transcriptome of hypertensive mice and control mice. We selected 3 hypertensive samples (thoracic aorta of C57BL/6J wild-type hypertensive mice infused with AngII (490 ng/kg/min) for 2 weeks) and 3 control samples (thoracic aorta of C57BL/6J wild-type hypertensive mice infused with an AngII carrier for 2 weeks).

2.2. Differential Gene Analysis

R software was used for statistical analysis and visualization of the data. The data were downloaded from the GEO database through the GEO query package. Probes with one probe corresponding to multiple molecules were removed. When probes corresponding to the same molecule were encountered, only probes with the maximum signal value were retained. The DESeq2 package was used to normalise the counts. The limma package was used for difference analysis between two groups, and data visualization was performed. The threshold for identifying DEGs was set to P < 0.05. A volcano map and heatmap were generated to display the differential gene expression results. Overlap analysis of the different gene sets in this paper was performed using Venny 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/index.html).

2.3. Gene Ontology (GO) Analysis

GO enrichment analysis was performed using R software. GO annotation analysis (molecular function, cellular component, and biological process) of the DEGs was performed using the cluster Profiler package, loading species genomes using the Annotation Hub package, and the database for GO analysis was called from https://www.geneontology.org/. P < 0.05 was considered to indicate statistical significance.

2.4. Immune Infiltration Assay

Immune infiltration analysis was performed using R software, and the gene expression data of each sample were compared with those of immune cells through the mMCP-counter package to assess the infiltration level of immune cells.

2.5. Gene Set Enrichment Analysis (GSEA)

R software was used for GSEA conversion of rat genes to mouse genes using the biomaRT software package and the cluster profiler package was used to analyze the risk signature-related biological pathways and functions of the DEGs. The GMT files used were obtained from https://www.gsea-msigdb.org/gsea/downloads.jsp, and the “gene set database” was selected as follows: “M2: curated gene sets (CP: Canonical pathways)”. Based on the gene expression profiles and phenotype grouping, according to the description of clusterprofiler (https://bioconductor.org/packages/release/bioc/html/clusterProfiler.html), the minimum gene set was set to 5, the maximum gene set was 5000, and one thousand repetitions were used for sampling. A P < 0.05 and an FDR < 0.25 were considered significant, and visual analysis were performed with the enrichplot package.

2.6. Clinical Sample Collection

A total of 102 blood samples were collected from healthy controls and hypertensive patients at Yangbi Yi Autonomous County People’s Hospital. The study was conducted with the informed consent of the patients and with the approval of the Ethics Committee of Yangbi Yi Autonomous County People’s Hospital.

2.7. RNA Extraction and RT-qPCR

Total RNA was extracted using a magnetic bead method and a total RNA extraction kit (G3611-50T, Servicebio, China) and synthesized using a SweScript RT I first-strand cDNA synthesis kit (G3331-50, Servicebio, China) according to the manufacturer’s instructions. cDNA was synthesized, and RT‒qPCR was performed using SYBR Green qPCR Master Mix (G3320-01, Servicebio, China). The RT-qPCR amplification reaction program was set as follows: amplification at 95˚C for 10 min, amplification at 95˚C for 15 s, amplification at 60˚C for 30 s, and 40 cycles. Using GAPDH as an internal reference gene, the relative expression levels were calculated by the 2−ΔΔct method. The sequence of primers is shown in Table 1.

3. Results3.1. DEGs That Were Significantly Differentially Expressed

A total of 2357 DEGs that were differentially expressed according to the GSE7483 dataset were obtained. A total of 1181 genes were upregulated, and 1176 genes were downregulated (Figure 1(a)). A total of 1416 DEGs were generated from the GSE75815 database, of which 735 were upregulated and 681 were downregulated (Figure 1(b)), additionally, heatmaps were drawn (Figure 1(c) and Figure 1(d)). Additionally, 267 overlapping DEGs were identified in the 2 datasets using Venny 2.1. (Figure 1(e)).

3.2. Enrichment Analysis

The results of GO enrichment analysis showed that the biological processes involved mainly included chemotactic positive regulation, actin filament organization, actin filament assembly, fibroblast proliferation, extracellular structure organization, extracellular matrix organization, positive regulation of response to external stimuli, regulation of vascular system development, regulation of angiogenesis, and chemotactic regulation. The main molecular functions include cell adhesion molecule binding, collagen binding, carbohydrate binding, sulfur binding, heparin binding, growth factor binding, extracellular matrix binding, extracellular matrix structural components, integrin binding, glycosaminoglycan

Primer sequences
GenesSequence (F: Forward primer, R: Reversed primer)
Ecm1F: 5’-AGAGCACTTTCAUGUGTTGG-3’
R: 5’-GAGTCAGGGTGATCCATGG-3’
Mecp2F: 5’-TGCTGGGAAGTATGATGTG-3’
R: 5’-CACTTTAGAGCGAAAGGCT-3’
SparcF: 5’-CCTGGTCACCCTGTATGAG-3’
R: 5’-TCATGGATCTTCTTCACCCG-3’
Sphk1F: 5’-TATGCTGGCTATGAGCAGG-3’
R: 5’-AGAGACAGCAGGTTCATGG-3’
Thbs1F: 5’-CGATAGCCTCAACAACCGA-3’
R: 5’-GCCATCCTGTTTAAUTCTC-3’
GAPDHF: 5’-TCAAGATCATCAGCAATGCC-3’
R: 5’-CGATACCAUAGTTGTCATGGA-3’

binding, etc. The cellular components mainly included collagen-containing extracellular matrix, cell-matrix junctions, β-catenin-TCF complexes, membrane domains of cell adhesion protein complexes, membrane microdomains of integrin complexes, platelet α-granules, basement membrane, etc. (Figures 2(a)-(f)). Through GO enrichment analysis, we found that the DEGs were significantly associated mainly with angiogenesis, vascular development, and intercellular junctions, and the key genes are listed in Table 2.

Key genes at the intersection of DEGs
NamelogFC
GSE7483GSE75815
Ecm10.4350.478
Mecp2−0.249−0.445
Sparc0.2530.916
Sphk10.5611.22
Thbs10.2050.595
3.3. Immune Infiltration Assay

To determine whether hypertensive patients have different immune infiltrates, we compared immune infiltration among hypertensive patients. As shown in Figure 3(a) and Figure 3(b), the immune cells (lymphatics, fibroblasts, T cells, B cells, vessels, endothelial cells, memory B cells, mast cells, granulocytes, NK cells, monocytes, cd8 T cells, neutrophils, macrophages and eosinophils) of the hypertension group and the normal group were significantly different, indicating that hypertension is significantly associated with immune infiltration. The GSEA technique was used to investigate the changes in patients with hypertension and the possible underlying mechanisms. The results showed that multiple signaling pathways were significantly enriched, but the enrichment modes were different. In the hypertension group, genes associated with keratinocyte, immunoregulatory interactions between a lymphoid cell and a nonlymphoid cell, degradation of the extracellular matrix, initial triggering of complement, degradation of the extracellular matrix, and cell surface interactions at the vascular wall were highly aggregated (Figure 3(c)). We selected five differential DEGs (Ecm1, Sparc, Sphk1, Thbsl and Mecp2) that are highly correlated with vascular development and angiogenesis for subsequent experimental validation.

3.4. RT-qPCR Validation of the Five DEGs in the Blood of Patients with Clinical Hypertension

To further validate the expression of the five novel key DEGs enriched above in the blood of hypertensive patients, we detected their expression in 51 hypertensive blood samples and 51 non-hypertensive blood samples. The RT‒qPCR results for Sparc (P < 0.05), Sphk1 (P < 0.05), and Thbs1 (P < 0.05) in the blood of hypertensive patients were consistent with the results of the database analysis, revealing a significant upregulation trend. In addition, RT-qPCR results showed that there was no significant difference in the expression of Ecm1 and Mecp2 in the blood of hypertensive patients and non-hypertensive blood (Figure 4).

4. Discussion

In the field of vascular biology, animal models play a critical role in the discovery of the basic mechanisms of vascular health and the development of treatments for various diseases [14] . The development and formation of blood vessels is a complex process that involves the regulation of various cell types and molecular signaling pathways and includes multiple steps, such as the proliferation and migration of vascular endothelial cells, tube formation, and vascular differentiation and maturation [15] [16] [17] . Studies have shown that vascular endothelial growth factor (VEGF) and its receptor (VEGFR) play key roles in tumor angiogenesis [18] . In addition, gestational hypertension is associated with changes in vascular function during pregnancy, and these changes may affect the future cardiovascular health of the mother and the child [19] . Clinical noninvasive methods for studying human endothelial function and vascular tension have proven that patients with hypertension exhibit impaired endothelium-dependent vasodilation, enhanced vasoreactivity, and increased contractility [20] . In addition, some studies have shown that promoting the development and formation of blood vessels can ameliorate the development and progression of hypertension [21] . Therefore, regulating angiogenesis may provide a new therapeutic strategy for the treatment of hypertension. In this study, we investigated the expression levels of the Ecm1, Mecp2, Sparc, Sphk1, and Thbs1 genes, which are associated with angiogenesis and angiogenesis in hypertensive patients, and proposed that Sparc, Sphk1 and Thbs1 may be significant novel targets.

SphK1 is a key rate-limiting enzyme and intracellular signal transduction enzyme involved in sphingolipid metabolism, and it is involved in the regulation of various physiological processes, such as Ca2+ homeostasis, cell survival, migration, and inflammation [22] [23] . Under inflammatory conditions, abnormally expressed SphK1 is involved in the regulation of the inflammatory response and the management of various immune cell functions. Abnormal inflammation can activate SphK1, destroying its activity and balance [24] [25] . Moreover, SphK1 triggers proinflammatory signals, regulates the chemotactic response of immune cells and activates vascular endothelial cells (VECs) [26] . Our studies showed that the expression of SphK1 in the blood of patients with hypertension was significantly upregulated. Further bioinformatics, KEGG, enrichment, and MF analyses confirmed that SphK1 was involved in the regulation of angiogenesis and angiogenesis inhibition in patients with hypertension.

Thrombospondin-1 (Thbs1) is associated with various cardiovascular diseases. The association between Thbs1 and endothelial dysfunction suggests that Thbs1 plays an important role in hypertension [27] . In addition, Thbs1 is an important endogenous inhibitor of angiogenesis that contributes to microcirculation and capillary density by regulating VEGF signaling [28] . Our results showed that Thbs1 levels increase under hypertensive conditions, which prompted us to further investigate the molecular importance of Thbs1 in the development of hypertension. Using bioinformatics methods, we used the GSE7483 and GSE75815 datasets to determine the Thbs1 level in patients with hypertension. Regulatory mechanisms of angiogenesis and angiogenesis. These results confirmed that Thbs1 obstructed the development and production of blood vessels in cells, thereby aggravating the increase in human blood pressure. Moreover, GO enrichment analysis and KEGG pathway analysis revealed that Ecm1, Mecp2, Sparc, Sphk1, and Thbs1 were involved in the processes of cellular immune infiltration and angiogenesis.

Hypertension is a common chronic disease that affects a large population worldwide. A growing body of evidence shows a strong link between the immune system and high blood pressure [29] [30] . In the pathological process of hypertension, immune cells infiltrate key organs, such as the blood vessel walls, heart, and kidney. These infiltrating immune cells, such as T lymphocytes and macrophages, can secrete a variety of cytokines and growth factors, which cause vasoconstriction, smooth muscle cell proliferation and inflammation and subsequently induce hypertension [31] [32] [33] . However, the association between immune infiltration and hypertension has not been fully elucidated. In the present study, the Ecm1, Mecp2, Sparc, Sphk1, and Thbs1 genes were found to be significantly correlated with the level of immune cell infiltration in patients with hypertension. Furthermore, the correlation of this feature with the response to immunotherapy was evaluated. The analysis of the significantly differentially expressed DEGs showed that Ecm1, Mecp2, Sparc, Sphk1, and Thbs1 were related to autophagy. Autophagy plays an important role in the development of cardiovascular diseases such as atherosclerosis, cardiac ischemia/reperfusion, cardiomyopathy, heart failure and hypertension, and generated reactive oxygen species interact to prevent hypertension [34] [35] . This finding is consistent with the results obtained from our DEGs dataset.

5. Summary

In summary, this study performed bioinformatics analysis of data from the GEO database and established a model of angiogenesis-related genes and the prognosis of patients with hypertension. Sparc, Sphk1 and Thbs1 could be used to predict the prognosis and immune response in patients with hypertension. The genes we identified can be used in the clinically stratified management of patients and can also lay the foundation for future studies of angiogenesis as a therapeutic target for hypertension. This study has several limitations. First, although new biomarkers associated with hypertension have been predicted, their mechanism of action remains unclear, and further study is needed. Second, the results need to be verified by experimental studies.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Chang, X.G., Guo, L., Zou, L.Y., Ma, Y.Z. and Feng, J.L. (2024) Identification and Validation of Vascular-Associated Biomarkers for the Prognosis and Potential Pathogenesis of Hypertension Using Comprehensive Bioinformatics Methods. World Journal of Cardiovascular Diseases, 14, 115-128. https://doi.org/10.4236/wjcd.2024.143010

NOTESReferencesLiu, Y., Kong, X., Jiang, Y., et al. (2022) Association of AGTR1 A1166C and CYP2C9*3 Gene Polymorphisms with the Antihypertensive Effect of Valsartan. International Journal of Hypertension, 2022, Article ID: 7677252. https://doi.org/10.1155/2022/7677252Kjeldsen, S. (2018) Hypertension and Cardiovascular Risk: General Aspects. Pharmacological Research, 129, 95-99. https://doi.org/10.1016/j.phrs.2017.11.003Lim, S., Lee, K., Han, T., et al. (2024) Antihypertensive Effect of Milk Fermented by Lactiplantibacillus plantarum K79 on Spontaneously Hypertensive Rats. Food Science of Animal Resources, 44, 178-188. https://doi.org/10.5851/kosfa.2023.e70World Health Organization (2021) WHO Guidelines Approved by the Guidelines Review Committee. Guideline for the Pharmacological Treatment of Hypertension in Adults. Geneva.Kim, B., Lee, Y., Shin, J., et al. (2024) Blood Pressure and Variability Responses to the Downtitration of Antihypertensive Drugs. Journal of Hypertension. https://doi.org/10.1097/HJH.0000000000003668Wang, Z., Liu, Y., Liu, J., et al. (2011) HSG/Mfn2 Gene Polymorphism and Essential Hypertension: A Case-Control Association Study in Chinese. Journal of Atherosclerosis and Thrombosis, 18, 24-31. https://doi.org/10.5551/jat.5611Holtstrater, C., Schrors, B., Bukur, T. and Lower, M. (2020) Bioinformatics for Cancer Immunotherapy. In: Boegel, S., Ed., Bioinformatics for Cancer Immunotherapy, Humana, New York, 1-9. https://doi.org/10.1007/978-1-0716-0327-7_1Ye, F., Huang, W. and Guo, G. (2017) Studying Hematopoiesis Using Single-Cell Technologies. Journal of Hematology & Oncology, 10, Article No. 27. https://doi.org/10.1186/s13045-017-0401-7Treppner, M., Salas-Bastos, A., Hess, M., et al. (2021) Synthetic Single Cell RNA Sequencing Data from Small Pilot Studies Using Deep Generative Models. Scientific Reports, 11, Article No. 9403. https://doi.org/10.1038/s41598-021-88875-4Li, D., Hao, X. and Song, Y. (2018) Identification of the Key MicroRNAs and the MiRNA-MRNA Regulatory Pathways in Prostate Cancer by Bioinformatics Methods. BioMed Research International, 2018, Article ID: 6204128. https://doi.org/10.1155/2018/6204128Zhu, Y., Yang, X. and Zu, Y. (2022) Integrated Analysis of WGCNA and Machine Learning Identified Diagnostic Biomarkers in Dilated Cardiomyopathy with Heart Failure. Frontiers in Cell and Developmental Biology, 10, Article 1089915. https://doi.org/10.3389/fcell.2022.1089915Guo, Q. (2023) Bioinformatics Analysis of the Diversity of Gut Microbiota and Different Microbiota on Insulin Resistance in Diabetes Mellitus Patients. Heliyon, 9, e22117. https://doi.org/10.1016/j.heliyon.2023.e22117Tian, Y., Yang, J., Lan, M., et al. (2020) Construction and Analysis of a Joint Diagnosis Model of Random Forest and Artificial Neural Network for Heart Failure. Aging, 12, 26221-26235. https://doi.org/10.18632/aging.202405Cleaver, O. (2021) Mouse Models of Vascular Development and Disease. Current Opinion in Hematology, 28, 179-188. https://doi.org/10.1097/MOH.0000000000000649Ahir, B., Engelhard, H. and Lakka, S. (2020) Tumor Development and Angiogenesis in Adult Brain Tumor: Glioblastoma. Molecular Neurobiology, 57, 2461-2478. https://doi.org/10.1007/s12035-020-01892-8Langen, U., Ayloo, S. and Gu, C. (2019) Development and Cell Biology of the Blood-Brain Barrier. Annual Review of Cell and Developmental Biology, 35, 591-613. https://doi.org/10.1146/annurev-cellbio-100617-062608Patel, S., Nilsson, M., Le, X., et al. (2023) Molecular Mechanisms and Future Implications of VEGF/VEGFR in Cancer Therapy. Clinical Cancer Research, 29, 30-39. https://doi.org/10.1158/1078-0432.CCR-22-1366Hicklin, D. and Ellis, L. (2005) Role of the Vascular Endothelial Growth Factor Pathway in Tumor Growth and Angiogenesis. Journal of Clinical Oncology, 23, 1011-1027. https://doi.org/10.1200/JCO.2005.06.081Pasqua, T., Pagliaro, P., Rocca, C., et al. (2018) Role of NLRP-3 Inflammasome in Hypertension: A Potential Therapeutic Target. Current Pharmaceutical Biotechnology, 19, 708-714. https://doi.org/10.2174/1389201019666180808162011Flammer, A., Anderson, T., Celermajer, D., et al. (2012) The Assessment of Endothelial Function: From Research Into Clinical Practice. Circulation, 126, 753-767. https://doi.org/10.1161/CIRCULATIONAHA.112.093245Touyz, R., Alves-Lopes, R., Rios, F., et al. (2018) Vascular Smooth Muscle Contraction in Hypertension. Cardiovascular Research, 114, 529-539. https://doi.org/10.1093/cvr/cvy023Pyne, N.J. and Pyne, S. (2020) Recent Advances in the Role of Sphingosine 1-Phosphate in Cancer. FEBS Letters, 594, 3583-3601. https://doi.org/10.1002/1873-3468.13933Mcginley, M.P. and Cohen, J.A. (2021) Sphingosine 1-Phosphate Receptor Modulators in Multiple Sclerosis and Other Conditions. The Lancet, 398, 1184-1194. https://doi.org/10.1016/S0140-6736(21)00244-0Obinata, H. and Hla, T. (2019) Sphingosine 1-Phosphate and Inflammation. International Immunology, 31, 617-625. https://doi.org/10.1093/intimm/dxz037Schneider, G. (2020) S1P Signaling in the Tumor Microenvironment. Advances in Experimental Medicine and Biology, 1223, 129-153. https://doi.org/10.1007/978-3-030-35582-1_7Wang, Y., Wu, H., Deng, R., et al. (2021) Geniposide Downregulates the VEGF/SphK1/S1P Pathway and Alleviates Angiogenesis in Rheumatoid Arthritis in Vivo and in Vitro. Phytotherapy Research, 35, 4347-4362. https://doi.org/10.1002/ptr.7130Sharma, K., Chanana, N., Mohammad, G., et al. (2021) Hypertensive Patients Exhibit Enhanced Thrombospondin-1 Levels at High-Altitude. Life, 11, Article 893. https://doi.org/10.3390/life11090893Kuebler, W.M. (2017) What Mediates the Effects of Thrombospondin-1 in Pulmonary Hypertension? New Evidence for a Dual-Pronged Role of CD47. Cardiovascular Research, 113, 3-5. https://doi.org/10.1093/cvr/cvw232Zhong, C., Si, Y., Yang, H., et al. (2023) Identification of Monocyte-Associated Pathways Participated in the Pathogenesis of Pulmonary Arterial Hypertension Based on Omics-Data. Pulmonary Circulation, 13, e12319. https://doi.org/10.1002/pul2.12319Carnevale, L., Perrotta, M., Mastroiacovo, F., et al. (2024) Advanced Magnetic Resonance Imaging (MRI) to Define the Microvascular Injury Driven by Neuroinflammation in the Brain of a Mouse Model of Hypertension. Hypertension, 81, 636-647. https://doi.org/10.1161/HYPERTENSIONAHA.123.21940Rodriguez-Iturbe, B., Pons, H. and Johnson, R. (2017) Role of the Immune System in Hypertension. Physiological Reviews, 97, 1127-1164. https://doi.org/10.1152/physrev.00031.2016Ferreira, N., Tostes, R., Paradis, P., et al. (2021) Aldosterone, Inflammation, Immune System, and Hypertension. American Journal of Hypertension, 34, 15-27. https://doi.org/10.1093/ajh/hpaa137Gan, L., Ye, D., Feng, Y., et al. (2023) Immune Cells and Hypertension. Immunologic Research, 72, 1-13. https://doi.org/10.1007/s12026-023-09414-zMei, Y., Thompson, M., Cohen, R., et al. (2015) Autophagy and Oxidative Stress in Cardiovascular Diseases. Biochimica et Biophysica Acta, 1852, 243-251. https://doi.org/10.1016/j.bbadis.2014.05.005Mao, M., Song, S., Li, X., et al. (2023) Advances in Epigenetic Modifications of Autophagic Process in Pulmonary Hypertension. Frontiers in Immunology, 14, Article 1206406. https://doi.org/10.3389/fimmu.2023.1206406