<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">Health</journal-id><journal-title-group><journal-title>Health</journal-title></journal-title-group><issn pub-type="epub">1949-4998</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/health.2024.163015</article-id><article-id pub-id-type="publisher-id">Health-131783</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Biomedical&amp;Life Sciences</subject><subject> Medicine&amp;Healthcare</subject></subj-group></article-categories><title-group><article-title>
 
 
  Identification and Validation of Novel Biomarkers Related to the Calcium Metabolism Pathway in Hypertension Patients Based on Comprehensive Bioinformatics Methods
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Xiangguang</surname><given-names>Chang</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Lei</surname><given-names>Guo</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Liying</surname><given-names>Zou</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Yazhao</surname><given-names>Ma</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Jilin</surname><given-names>Feng</given-names></name><xref ref-type="aff" rid="aff4"><sup>4</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Laboratory, Yangbi Yi Autonomous County People’s Hospital, Dali, China</addr-line></aff><aff id="aff3"><addr-line>Department of Infectious Diseases, Yangbi Yi Autonomous County People’s Hospital, Dali, China</addr-line></aff><aff id="aff4"><addr-line>Medical Community Office, Yangbi Yi Autonomous County People’s Hospital, Dali, China</addr-line></aff><aff id="aff2"><addr-line>Department of Obstetrics and Gynecology, Yangbi Yi Autonomous County People’s Hospital, Dali, China</addr-line></aff><pub-date pub-type="epub"><day>18</day><month>03</month><year>2024</year></pub-date><volume>16</volume><issue>03</issue><fpage>173</fpage><lpage>186</lpage><history><date date-type="received"><day>6,</day>	<month>February</month>	<year>2024</year></date><date date-type="rev-recd"><day>15,</day>	<month>March</month>	<year>2024</year>	</date><date date-type="accepted"><day>18,</day>	<month>March</month>	<year>2024</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  <b>Background:</b> Hypertension is a universal risk factor for cardiovascular diseases and is thus the leading cause of death worldwide. The identification of 
  novel prognostic and pathogenesis biomarkers plays a key role in disease 
  management. <b>Methods:</b> The GSE145854 and GSE164494 datasets were downloa
  ded from the Gene Expression Omnibus (GEO) database and used for s
  creening and validating hypertension signature genes, respectively. Gene Ontology (GO) enrichment analysis was performed on the differentially expressed genes (DEGs) related to calcium ion metabolism in patients with hypertension. The core genes related to immune infiltration were analyzed and screened, and the activity of the signature genes and related pathways was quantified using gene set enrichment analysis (GSEA). The infiltration of immune cells in the blood samples was analyzed, and the DEGs that were abnormally expressed in the clinical blood samples of patients with hypertension were verified via RT
  -
  qPCR. <b>Results:</b> A total of 176 DEGs were screened. GO showed that DEGs was involved in the regulation of calcium ion metabolism in biological processes (BP), actin mediated cell contraction, negative regulation of cell movement, and calcium ion transmembrane transport, and in the regulation of protease activity in molecular functions (MF). KEGG analysis revealed that the DEGs were involved mainly in the cGMP
  -
  PKG signalin
  g pathway, ubiquitin-protein transferase, tight junction-associated pro
  teins, 
  and the regulation of myocardial cells. MF analysis revealed the
   
  immune infiltration function of the cells. RT
  -
  qPCR revealed that the expression of Cacna1d, Serpine1, Slc8a3, and Trpc4 was up
   
  regulated in hypertension,
   
  the expression of Myoz2 and Slc25a23 was down
   
  regulated. <b>Conclusion:</b> Cacna1d, Serpine1, Slc8a3, Trpc4, Myoz2 and Slc25a23 may be involved in the regulation of calcium metabolism pathways and play key roles in hypertension. These differentially expressed calcium metabolism-related genes may serve as prognostic markers of hypertension.
 
</p></abstract><kwd-group><kwd>Hypertension</kwd><kwd> Biomarkers</kwd><kwd> Differentially Expressed Genes</kwd><kwd> Ca&lt;sup&gt;2+&lt;/sup&gt; Metabo-lism</kwd><kwd> Bioinformatics Analysis</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Hypertension is a worldwide health problem and one of the common causes of death. High blood pressure is often an underlying and progressive process, so it is often difficult to diagnose in a timely manner [<xref ref-type="bibr" rid="scirp.131783-ref1">1</xref>] . Its pathogenesis is complex and includes genetic susceptibility and the influence of environmental factors [<xref ref-type="bibr" rid="scirp.131783-ref2">2</xref>] . In addition, studies have shown that hypertension is associated with the development of cardiovascular diseases, which increase the risk of cardiovascular diseases and their complications and may also lead to cognitive and functional decline in elderly individuals [<xref ref-type="bibr" rid="scirp.131783-ref3">3</xref>] . However, by lowering blood pressure to a target level, or “controlling” high blood pressure with lifestyle changes and medication, the risk of adverse outcomes from high blood pressure can be greatly reduced [<xref ref-type="bibr" rid="scirp.131783-ref4">4</xref>] . Although drug treatment has been used for a long time, in many cases, patients’ blood pressure still cannot reach the ideal range [<xref ref-type="bibr" rid="scirp.131783-ref5">5</xref>] .</p><p>In recent years, the rapid development of bioinformatics technology has provided important support for the identification of biomarkers and the identification of hub genes in the module of clinical significance [<xref ref-type="bibr" rid="scirp.131783-ref6">6</xref>] . Through comprehensive analysis of high-throughput data, bioinformatics and machine learning analyses can aid in understanding the molecular mechanisms of this disease and screening key genes [<xref ref-type="bibr" rid="scirp.131783-ref7">7</xref>] . These methods have played key roles in transcriptomic, metabolomic and proteomic studies, facilitating the search for disease driver genes, discovery of drug targets, classification of tumor subtypes, discovery of molecular markers, drug sensitivity screening, etc., and querying and predicting protein and RNA interactions [<xref ref-type="bibr" rid="scirp.131783-ref8">8</xref>] [<xref ref-type="bibr" rid="scirp.131783-ref9">9</xref>] [<xref ref-type="bibr" rid="scirp.131783-ref10">10</xref>] .</p><p>This study aimed to provide basic insights into the etiology and related molecular basis of hypertension by analyzing the differentially expressed genes (DEGs) through comprehensive bioinformatics. Our study aimed to elucidate the role of calcium ion metabolism-related genes in the pathogenesis of hypertension and to create a strategy for potential drug prediction.</p></sec><sec id="s2"><title>2. Materials and Methodology</title><sec id="s2_1"><title>2.1. Sample Collection</title><p>We collected blood samples from 51 non-hypertensive people (Age: 9 - 99; Male: n = 25; Female: n = 26) and 51 hypertensive patients (Age: 43 - 88; Male: n = 27; Female: n = 24) in 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.</p></sec><sec id="s2_2"><title>2.2. RT-qPCR</title><p>According to the manufacturer’s instructions, we used a magnetic bead method and total RNA extraction kit (G3611-50T; Servicebio, China) to extract total RNA from blood and a SweScript RT I first-strand cDNA synthesis kit (G3331-50; Servicebio, China) for cDNA synthesis. In subsequent experiments, we used SYBR Green qPCR Master Mix (G3320-01; Servicebio, China) for RT-qPCR. The RT-qPCR amplification reaction program was as follows. First, the amplification reaction was performed at 95˚C for 10 min, followed by 40 cycles of amplification at 95˚C for 15 s and amplification at 60˚C. The amplification reaction was performed at 60˚C for 30 s. For the quantitative analysis of gene expression levels, we used GAPDH as an internal reference gene and the 2<sup>−ΔΔCt</sup> method to calculate the relative expression level. The sequences of primers used were as follows: (<xref ref-type="table" rid="table1">Table 1</xref>).</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Primer sequences</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Genes</th><th align="center" valign="middle" >Sequence (F: Forward primer, R: Reversed primer)</th></tr></thead><tr><td align="center" valign="middle"  rowspan="2"  >Cacna1d</td><td align="center" valign="middle" >F: 5'-AAATCCAAACTCAGCCGAC-3'</td></tr><tr><td align="center" valign="middle" >R: 5'-AAACGTGACAGACTTCACG-3'</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Cbarp</td><td align="center" valign="middle" >F: 5'-GCCCCAUTTGTCACCAT-3'</td></tr><tr><td align="center" valign="middle" >R: 5'-CTGATCTCTGCGAAGTCAGTG-3'</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Ccl19</td><td align="center" valign="middle" >F: 5'-CCTGCTGTAGTGTTCACCA-3'</td></tr><tr><td align="center" valign="middle" >R: 5'-TCTGGATGATGCGTTCTACC-3'</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Serpine 1</td><td align="center" valign="middle" >F: 5'-TTCAAGTTGATGACAGGGC-3'</td></tr><tr><td align="center" valign="middle" >R: 5'-CTCATCCTTGTTCCATGGC-3'</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Slc8a3</td><td align="center" valign="middle" >F: 5'-GATGGGAAAGCCAGTATTGG-3'</td></tr><tr><td align="center" valign="middle" >R: 5'-TCCACCGTAGTCTTGAACTC-3'</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Trpc4</td><td align="center" valign="middle" >F: 5'-CATTTGTAAGTACAGTGCCC-3'</td></tr><tr><td align="center" valign="middle" >R: 5'-CAAATAAAGCCTCTGCCAC-3'</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Gt2i</td><td align="center" valign="middle" >F: 5'-TGCAAAGGAAAGGATTCGT-3'</td></tr><tr><td align="center" valign="middle" >R: 5'-CTTCCTTTACTCCTGAAGCTG-3'</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Myoz2</td><td align="center" valign="middle" >F: 5'-TGCTCCAGGATATTCTGGAC-3'</td></tr><tr><td align="center" valign="middle" >R: 5'-GATTGATAGTACTTAGGGACAGC-3'</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Nos1ap</td><td align="center" valign="middle" >F: 5'-GGATACGGTATGAGTTAUGC-3'</td></tr><tr><td align="center" valign="middle" >R: 5'-TCACTTTCACTCCATCCAC-3'</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Slc25a23</td><td align="center" valign="middle" >F: 5'-TCTACGAGACTCTGAAGUACTG-3'</td></tr><tr><td align="center" valign="middle" >R: 5'-CTGGATATGGTACCGCAGG-3'</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Spink 1</td><td align="center" valign="middle" >F: 5'-TTTCTTCTCAGTGCCTTGG-3'</td></tr><tr><td align="center" valign="middle" >R: 5'-CATTGTAACATTTGGCCTCTC-3'</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >GAPDH</td><td align="center" valign="middle" >F: 5'-TCAAGATCATCAGCAATGCC-3'</td></tr><tr><td align="center" valign="middle" >R: 5'-CGATACCAUAGTTGTCATGGA-3'</td></tr></tbody></table></table-wrap></sec><sec id="s2_3"><title>2.3. Determination of DEGs</title><p>When using R software for statistical analysis and visualization, we used the GEOquery package to download the hypertension-related datasets GSE145854 and GSE164494 from the GEO database (https://www.ncbi.nlm.nih.gov/geo/). GSE145854 using the GPL17117 platform, we describe the transcriptome sequencing results of Wistar Kyoto (WKY) Rats and Spontaneously Hypertensive Rats (SHR) thoracic aorta, where we used 3 SHR samples as hypertensive samples and 3 WKY samples as normal controls. GSE164494 using the GPL21273 platform and the Ang II-induced hypertensive mouse model, we used 5 wild-type mice treated with 500 ng/kg/min AngII infusion as hypertensive samples and 5 wild-type mice treated with saline infusion as normal controls. During data processing, we excluded the case in which one probe corresponded to multiple molecules and when multiple probes corresponding to the same molecule were encountered, we retained only the probe with the maximum signal value. Next, we analysed the differences between the two samples using the limma package and normalised them using the DESeq2 package. The threshold for identifying DEGs was set at a P &lt; 0.05. We created volcano plots and heatmaps to visualize the expression patterns of the DEGs. To further understand the relationships between different gene sets, we performed overlap analysis using the Venny 2.1 tool (https://bioinfogp.cnb.csic.es/tools/venny/index.html). This approach helped to reveal common and unique gene sets that were biologically important in our study.</p></sec><sec id="s2_4"><title>2.4. Gene Ontology (GO) Analysis</title><p>We used R software for GO enrichment analysis, and the cluster Profiler package was used for GO annotation analysis of the DEGs, which included Molecular Function (MF), Cellular Component (CC), and Biological Process (BP) terms. In the KEGG pathway enrichment analysis, we also used the cluster Profiler package, and a significance level of P &lt; 0.05 was used as the standard for statistical significance.</p></sec><sec id="s2_5"><title>2.5. Gene Set Enrichment Analysis (GSEA)</title><p>Genome enrichment analysis (GSEA) was performed using R software to identify biological pathways and functions associated with phenotype groupings. Specifically, the gene expression profile and phenotype grouping data were organized into an appropriate format, the risk characteristic-related biological pathway and functional analysis of the DEGs were performed using the cluster Profiler package. To investigate the association of differential DEGs with relevant signalling pathways and biological processes in hypertension, the gene set database was selected as “M2: curated gene sets (CP: Canonical pathways)”. The minimum gene set was 5, the maximum gene set was 5000, the resampling was 1,000 times, the significance thresholds were set to P &lt; 0.05 and FDR &lt; 0.25, and the enrichplot package was used to visualize the GSEA results.</p></sec><sec id="s2_6"><title>2.6. Immune Infiltration Assay</title><p>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.</p></sec></sec><sec id="s3"><title>3. Results</title><sec id="s3_1"><title>3.1. DEGs that were Significantly Differentially Expressed</title><p>According to the search terms and retrieval conditions, the GSE145854 and GSE164494 datasets were downloaded from the GEO database. Among them, GSE145854 had 570 upregulated DEGs and 425 downregulated DEGs (<xref ref-type="fig" rid="fig1">Figure 1</xref>(A)), and GSE164494 had 2278 upregulated DEGs and 2301 downregulated DEGs (<xref ref-type="fig" rid="fig1">Figure 1</xref>(B)). In addition, 176 overlapping DEGs were identified in the 2 datasets using Venny 2.1 (<xref ref-type="fig" rid="fig1">Figure 1</xref>(C)).</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Key genes involved in calcium ion metabolism</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >name</th><th align="center" valign="middle"  colspan="2"  >logFC</th></tr></thead><tr><td align="center" valign="middle" >GSE164494</td><td align="center" valign="middle" >GSE145854</td></tr><tr><td align="center" valign="middle" >Cacna1d</td><td align="center" valign="middle" >2.767629029</td><td align="center" valign="middle" >1.067715</td></tr><tr><td align="center" valign="middle" >Cbarp</td><td align="center" valign="middle" >1.038753192</td><td align="center" valign="middle" >0.396428</td></tr><tr><td align="center" valign="middle" >Ccl19</td><td align="center" valign="middle" >3.583316657</td><td align="center" valign="middle" >0.475507333</td></tr><tr><td align="center" valign="middle" >Gt2i</td><td align="center" valign="middle" >−1.855870609</td><td align="center" valign="middle" >−0.307862667</td></tr><tr><td align="center" valign="middle" >Myoz2</td><td align="center" valign="middle" >−1.654615041</td><td align="center" valign="middle" >−0.877018667</td></tr><tr><td align="center" valign="middle" >Nos1ap</td><td align="center" valign="middle" >−1.097176652</td><td align="center" valign="middle" >−0.553486333</td></tr><tr><td align="center" valign="middle" >Serpine 1</td><td align="center" valign="middle" >1.432787915</td><td align="center" valign="middle" >0.951215</td></tr><tr><td align="center" valign="middle" >Slc25a23</td><td align="center" valign="middle" >−2.905241285</td><td align="center" valign="middle" >−0.714465</td></tr><tr><td align="center" valign="middle" >Slc8a3</td><td align="center" valign="middle" >2.344115167</td><td align="center" valign="middle" >0.775163333</td></tr><tr><td align="center" valign="middle" >Spink 1</td><td align="center" valign="middle" >−2.262386725</td><td align="center" valign="middle" >−0.563688333</td></tr><tr><td align="center" valign="middle" >Trpv6</td><td align="center" valign="middle" >−1.099406098</td><td align="center" valign="middle" >0.303186</td></tr></tbody></table></table-wrap></sec><sec id="s3_2"><title>3.2. Enrichment Analysis</title><p>To demonstrate the hypertension-related functional annotation and pathway enrichment, GO function and KEGG pathway analyses were performed on the DEGs. The GO annotation of DEGs consisted of three parts: biological process (BP), cellular component (CC) and molecular function (MF). The functional enrichment of the DEGs was analyzed. The DEGs were involved mainly in the regulation of calcium ion transmembrane transport, the immune response, cell physiological functions, cell migration, and cytokine production (Figures 2(A)-(D)). The results of KEGG analysis showed that these genes were mainly associated with the interactions of the cGMP-PKG signaling pathway, ubiquitin-protein transferase, tight junction-related proteins and the regulation of myocardial cells (<xref ref-type="fig" rid="fig2">Figure 2</xref>(E)). The intersection DEGs were significantly associated with calcium ion metabolism and were simultaneously enriched in other metabolic pathways. The related key genes are listed in the table below (<xref ref-type="table" rid="table2">Table 2</xref>).</p></sec><sec id="s3_3"><title>3.3. Immune Infiltration Assay</title><p>The GSEA technique was used to investigate the changes and possible underlying mechanisms in hypertensive patients and healthy volunteers. In hypertensive patients, Role of Phospholipids in Phagocytosis, Mitotic Prometaphase, Factors Involved in Megakaryocyte Development and Platelet Production, Role of Lat2 Ntal Lab on Calcium Mobilization, Fceri Mediated Mapk Activation, Chromatin Modifying Enzymes pathway was significantly enriched, but the enrichment patterns were different between the two groups (<xref ref-type="fig" rid="fig3">Figure 3</xref>(A)). To determine whether hypertensive patients have different immune infiltrates, we compared immune infiltration among hypertensive patients. Immune infiltration assay is mainly based on analysis of tissue sample transcriptome sequencing data or gene chip data. Different cells have some signature markers, as do immune cells, so immune cell infiltration can be assessed and quantified based on the amount of expression of these marker genes in tissues. The results of our analysis show that the expression of monocytes, T cells, mast cells, B cells, endothelial cells, basophils, memory B cells, lymphatic vessels, fibroblasts, cd8 T cells, eosinophils, neutrophils, NK cells, macrophages, and vessels were significantly different between the control group and the hypertensive patients (<xref ref-type="fig" rid="fig3">Figure 3</xref>(B)). These results suggest that this feature model has a potential role in predicting the immune response to immunotherapy in hypertensive patients.</p></sec><sec id="s3_4"><title>3.4. RT-qPCR validation of 11 DEGs in the Blood of Patients with Clinical Hypertension</title><p>To further verify the expression levels of 11 novel DEGs in the blood of hypertensive patients. RT-qPCR analysis was performed on 51 blood samples from hypertensive patients and non-hypertensive patients. Cacna1d, Cbarp, Ccl19, Serpine1, Slc8a3, and Trpc4 tended to be upregulated, while Cacna1d (P &lt; 0.05), Serpine1 (P &lt; 0.05), Slc8a3 (P &lt; 0.05), and Trpc4 (P &lt; 0.05) were significantly upregulated. Gtf2i, Myoz2, Nos1ap, Slc25a23, and Spink1 tended to be downregulated, and Myoz2 (P &lt; 0.05) and Slc25a23 (P &lt; 0.05) had significant trends (<xref ref-type="fig" rid="fig4">Figure 4</xref>).</p></sec></sec><sec id="s4"><title>4. Discussion</title><p>According to the Global Health Observatory, 1.13 billion people worldwide have hypertension, which increases the incidence of diseases such as left ventricular hypertrophy, coronary heart disease, heart failure, atrial fibrillation and peripheral arterial disease [<xref ref-type="bibr" rid="scirp.131783-ref11">11</xref>] [<xref ref-type="bibr" rid="scirp.131783-ref12">12</xref>] . Many pathophysiological factors are known to be associated with the pathogenesis of hypertension, including structural and functional abnormalities as well as molecular and cellular mechanisms of cardiovascular changes [<xref ref-type="bibr" rid="scirp.131783-ref13">13</xref>] [<xref ref-type="bibr" rid="scirp.131783-ref14">14</xref>] . In addition, impaired vasodilation [<xref ref-type="bibr" rid="scirp.131783-ref15">15</xref>] , impaired Ca<sup>2+</sup> signaling [<xref ref-type="bibr" rid="scirp.131783-ref16">16</xref>] , oxidative stress [<xref ref-type="bibr" rid="scirp.131783-ref17">17</xref>] , proinflammatory cytokines [<xref ref-type="bibr" rid="scirp.131783-ref18">18</xref>] and fibrotic growth factor [<xref ref-type="bibr" rid="scirp.131783-ref19">19</xref>] are considered to play important roles in this process. Ca<sup>2+</sup> serve as universal secondary messengers and participate in various biological processes, such as proliferation, cell death, migration and immune response [<xref ref-type="bibr" rid="scirp.131783-ref20">20</xref>] . There is now growing evidence that calcium ion homeostasis is an important driving factor in the occurrence and development of cancers and affects the treatment response of cancer patients [<xref ref-type="bibr" rid="scirp.131783-ref21">21</xref>] . With the understanding of cellular calcium homeostasis mechanisms and specific calcium signaling-targeted therapies (e.g., cardiovascular disease [<xref ref-type="bibr" rid="scirp.131783-ref22">22</xref>] and neuronal diseases [<xref ref-type="bibr" rid="scirp.131783-ref23">23</xref>] ), Ca<sup>2+</sup> signaling has become a very attractive target for the development of novel disease treatment drugs [<xref ref-type="bibr" rid="scirp.131783-ref24">24</xref>] . The researchers used bioinformatics technology to analyze the molecules that regulate calcium metabolism mediating hypertension, and preliminary screening analysis found that core targets such as Cacna1d, Cbarp, Ccl19, Serpine1, Slc8a3, Trpc4, Gtf2i, Myoz2, Nos1ap, Slc25a23, and Spink1 may regulate calcium metabolism and hypertension.</p><p>Bioinformatics methods have been used to analyze gene sequencing data of various diseases to identify DEGs and perform various analyses. An increasing number of powerful databases and powerful online tools have been established to help us reposition known complex disease mechanisms [<xref ref-type="bibr" rid="scirp.131783-ref25">25</xref>] . In a growing body of studies, microarray technology is being used to search for DEGs and their molecular functions (MFs), biological processes (BPs) or cellular components (CCs), and related regulatory pathways in specific disease states [<xref ref-type="bibr" rid="scirp.131783-ref26">26</xref>] . For example, Habib Rahman and his colleagues used bioinformatics and machine learning methods to identify novel factors that improve the identification and characterization of glioblastoma tumors and their progression [<xref ref-type="bibr" rid="scirp.131783-ref27">27</xref>] . With neighborhood-based benchmarking and multilayer network topology techniques, they also identified novel putative biomarkers showing how T2D interacts [<xref ref-type="bibr" rid="scirp.131783-ref28">28</xref>] . In our study, key genes were identified via comprehensive bioinformatics analysis techniques, such as GO and KEGG enrichment analysis and MF analysis, in the GEO datasets GSE145854 and GSE164494. Relying on the GEO database, the specific infiltration of immune cells in the blood of hypertensive patients and healthy patients was compared to elucidate the pathogenesis and potential therapeutic targets of hypertension.</p><p>The results of this study revealed 11 common DEGs between the GSE145854 dataset and the GSE164494 dataset, including 6 upregulated DEGs, i.e., Cacna1d, Cbarp, Ccl19, Serpine1, Slc8a3, and Trpc4, and 5 downregulated DEGs, i.e., Gtf2i, Myoz2, Nos1ap, Slc25a23 and Spink1. A recent genome-wide association study (GWAS) of Nos1ap, Slc25a23, and Spink1 revealed that the rs9810888 polymorphism of the calcium voltage-gated channel subunit α1 D gene (Cacna1d) was associated with blood pressure in Chinese adults [<xref ref-type="bibr" rid="scirp.131783-ref29">29</xref>] . Stanton AM et al. also showed that the Cacna1d allele is associated with elevated blood pressure and blood pressure salt sensitivity and that obstruction of Cacna1d expression can relieve the upregulation of blood pressure [<xref ref-type="bibr" rid="scirp.131783-ref30">30</xref>] . Cacna1d is an important gene in the calcium channel pathway. It plays an important role in the regulation of cellular calcium and iron levels and is a key signal for the contraction and dilation of vascular smooth muscle [<xref ref-type="bibr" rid="scirp.131783-ref31">31</xref>] . Our results also indicated that Cacna1d was involved in the regulation of calcium ion metabolism and the process of hypertension.</p><p>In summary, this study conducted bioinformatics analysis based on GEO data set, and we found that genes Cacna1d, Serpine1, Slc8a3, Trpc4, Myoz2 and Slc25a23 were the most significant markers of calcium ion metabolism and hypertension. In addition, down-regulating the expression of Cacna1d, Serpine1, Slc8a3, Trpc4, or promoting the expression of Myoz2 and Slc25a23 are molecules that regulate calcium ion metabolism and predict the occurrence of hypertension, and can also be considered as potential targets to prevent the concentration of calcium ions in the blood vessel wall and calcium ion metabolism, so as to prevent the occurrence of vascular calcification and vascular diseases. Although it is currently unclear how the immune system affects vascular remodeling and calcium ion metabolism, it can be speculated that the immune system has a long way to go in this pathological process. The study has the problems of insufficient data and result bias, so further experiments are needed to verify and improve the study results. Further experiments on immune cells can help to identify the specific targets of the immune system involved in vascular remodeling and calcium ion metabolism, facilitating the development of immunomodulatory therapy for patients with hypertension. Moreover, it is necessary to explore the specific underlying mechanism through animal and cell experiments.</p></sec><sec id="s5"><title>Conflicts of Interest</title><p>The authors declare no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s6"><title>Cite this paper</title><p>Chang, X.G., Guo, L., Zou, L.Y., Ma, Y.Z. and Feng, J.L. (2024) Identification and Validation of Novel Biomarkers Related to the Calcium Metabolism Pathway in Hypertension Patients Based on Comprehensive Bioinformatics Methods. Health, 16, 173-186. https://doi.org/10.4236/health.2024.163015</p></sec><sec id="s7"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.131783-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Smereczańska, M., Domian, N., M&amp;#322;ynarczyk, M., et al. (2023) Evaluation of the Expression and Localization of the Multifunctional Protein CacyBP/SIP and Elements of the MAPK Signaling Pathway in the Adrenal Glands of Rats with Primary and Secondary Hypertension. International Journal of Molecular Sciences, 25, Article No. 84. &lt;/br&gt;https://doi.org/10.3390/ijms25010084</mixed-citation></ref><ref id="scirp.131783-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Gharipour, M., Nezafati, P., Sadeghian, L., et al. (2022) Precision Medicine and Metabolic Syndrome. ARYA Atherosclerosis, 18, 1-10.</mixed-citation></ref><ref id="scirp.131783-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Buford, T. (2016) Hypertension and Aging. Aging Research Reviews, 26, 96-111. &lt;/br&gt;https://doi.org/10.1016/j.arr.2016.01.007</mixed-citation></ref><ref id="scirp.131783-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Zhou, D., Xi, B., Zhao, M., et al. (2018) Uncontrolled Hypertension Increases Risk of All-Cause and Cardiovascular Disease Mortality in US Adults: The NHANES III Linked Mortality Study. Scientific Reports, 8, Article No. 9418. &lt;/br&gt;https://doi.org/10.1038/s41598-018-27377-2</mixed-citation></ref><ref id="scirp.131783-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Deussen, A. and Kopaliani, I. (2023) Targeting Inflammation in Hypertension. Current Opinion in Nephrology and Hypertension, 32, 111-117. &lt;/br&gt;https://doi.org/10.1097/MNH.0000000000000862</mixed-citation></ref><ref id="scirp.131783-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Zhu, 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 ID: 1089915. &lt;/br&gt;https://doi.org/10.3389/fcell.2022.1089915</mixed-citation></ref><ref id="scirp.131783-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Orlov, Y., Anashkina, A., Klimontov, V., et al. (2021) Medical Genetics, Genomics and Bioinformatics Aid in Understanding Molecular Mechanisms of Human Diseases. International Journal of Molecular Sciences, 22, Article No. 9962. &lt;/br&gt;https://doi.org/10.3390/ijms22189962</mixed-citation></ref><ref id="scirp.131783-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Picard, M., Scott-Boyer, M., Bodein, A., et al. (2021) Integration Strategies of Multiomics Data for Machine Learning Analysis. Computational and Structural Biotechnology Journal, 19, 3735-3746. &lt;/br&gt;https://doi.org/10.1016/j.csbj.2021.06.030</mixed-citation></ref><ref id="scirp.131783-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Zheng, T., Zheng, Z., Zhou, H., et al. (2023) The Multifaceted Roles of COL4A4 in Lung Adenocarcinoma: An Integrated Bioinformatics and Experimental Study. Computers in Biology and Medicine, 170, Article ID: 107896. &lt;/br&gt;https://doi.org/10.1016/j.compbiomed.2023.107896</mixed-citation></ref><ref id="scirp.131783-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Zhou, X., Chen, Y., Zhang, Z., et al. (2024) Identification of Differentially Expressed Genes, Signaling Pathways and Immune Infiltration in Postmenopausal Osteoporosis by Integrated Bioinformatics Analysis. Heliyon, 10, E23794. &lt;/br&gt;https://doi.org/10.1016/j.heliyon.2023.e23794</mixed-citation></ref><ref id="scirp.131783-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Klinik Für Innere Medizin III A E N, Universit&amp;#228;tsklinikum Jena, Jena, Germany (2019) Recent Advances in Hypertension Research. Acta Physiologica (Oxford, England), 226, E13295.</mixed-citation></ref><ref id="scirp.131783-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Wei, Y., George, N., Chang, C., et al. (2017) Assessing Sex Differences in the Risk of Cardiovascular Disease and Mortality per Increment in Systolic Blood Pressure: A Systematic Review and Meta-Analysis of Follow-Up Studies in the United States. PLOS ONE, 12, E0170218. &lt;/br&gt;https://doi.org/10.1371/journal.pone.0170218</mixed-citation></ref><ref id="scirp.131783-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Eid, A., El-Yazbi, A., Zouein, F., et al. (2018) Inositol 1,4,5-Trisphosphate Receptors in Hypertension. Frontiers in Physiology, 9, Article No. 1018. &lt;/br&gt;https://doi.org/10.3389/fphys.2018.01018</mixed-citation></ref><ref id="scirp.131783-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Ma, J., Li, Y., Yang, X., et al. (2023) Signaling Pathways in Vascular Function and Hypertension: Molecular Mechanisms and Therapeutic Interventions. Signal Transduction and Targeted Therapy, 8, Article No. 168. &lt;/br&gt;https://doi.org/10.1038/s41392-023-01430-7</mixed-citation></ref><ref id="scirp.131783-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Pan, X., Shao, Y., Wu, F., et al. (2018) FGF21 Prevents Angiotensin II-Induced Hypertension and Vascular Dysfunction by Activation of ACE2/Angiotensin-(1-7) Axis in Mice. Cell Metabolism, 27, 1323-1337.E5. &lt;/br&gt;https://doi.org/10.1016/j.cmet.2018.04.002</mixed-citation></ref><ref id="scirp.131783-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Wilson, C., Zhang, X., Buckley, C., et al. (2019) Increased Vascular Contractility in Hypertension Results from Impaired Endothelial Calcium Signaling. Hypertension (Dallas, Tex: 1979), 74, 1200-1214. &lt;/br&gt;https://doi.org/10.1161/HYPERTENSIONAHA.119.13791</mixed-citation></ref><ref id="scirp.131783-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Zhang, Z., Zhao, L., Zhou, X., et al. (2022) Role of Inflammation, Immunity, and Oxidative Stress in Hypertension: New Insights and Potential Therapeutic Targets. Frontiers in Immunology, 13, Article ID: 1098725. &lt;/br&gt;https://doi.org/10.3389/fimmu.2022.1098725</mixed-citation></ref><ref id="scirp.131783-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Pioli, M. and De Faria, A. (2019) Pro-Inflammatory Cytokines and Resistant Hypertension: Potential for Novel Treatments? Current Hypertension Reports, 21, Article No. 95. &lt;/br&gt;https://doi.org/10.1007/s11906-019-1003-2</mixed-citation></ref><ref id="scirp.131783-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Teixeira, D., Peruchetti, D., Souza, M., et al. (2020) A High Salt Diet Induces Tubular Damage Associated with a Pro-Inflammatory and Pro-Fibrotic Response in a Hypertension-Independent Manner. Biochimica et Biophysica Acta Molecular Basis of Disease, 1866, Article ID: 165907. &lt;/br&gt;https://doi.org/10.1016/j.bbadis.2020.165907</mixed-citation></ref><ref id="scirp.131783-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Giorgi, C., Danese, A., Missiroli, S., et al. (2018) Calcium Dynamics as a Machine for Decoding Signals. Trends in Cell Biology, 28, 258-273. &lt;/br&gt;https://doi.org/10.1016/j.tcb.2018.01.002</mixed-citation></ref><ref id="scirp.131783-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Marchi, S., Giorgi, C., Galluzzi, L., et al. (2020) Ca Fluxes and Cancer. Molecular Cell, 78, 1055-1069. &lt;/br&gt;https://doi.org/10.1016/j.molcel.2020.04.017</mixed-citation></ref><ref id="scirp.131783-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Lei, M., Wu, L., Terrar, D., et al. (2018) Modernized Classification of Cardiac Antiarrhythmic Drugs. Circulation, 138, 1879-1896. &lt;/br&gt;https://doi.org/10.1161/CIRCULATIONAHA.118.035455</mixed-citation></ref><ref id="scirp.131783-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Tiscione, S., Casas, M., Horvath, J., et al. (2021) IP3R-Driven Increases in Mitochondrial Ca2+ Promote Neuronal Death in NPC Disease. Proceedings of the National Academy of Sciences of the United States of America, 118, e2110629118. &lt;/br&gt;https://doi.org/10.1073/pnas.2110629118</mixed-citation></ref><ref id="scirp.131783-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Cui, C., Merritt, R., Fu, L., et al. (2017) Targeting Calcium Signaling in Cancer Therapy. Acta Pharmaceutica Sinica B, 7, 3-17. &lt;/br&gt;https://doi.org/10.1016/j.apsb.2016.11.001</mixed-citation></ref><ref id="scirp.131783-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Stearman, R., Bui, Q., Speyer, G., et al. (2019) Systems Analysis of the Human Pulmonary Arterial Hypertension Lung Transcriptome. American Journal of Respiratory Cell and Molecular Biology, 60, 637-649. &lt;/br&gt;https://doi.org/10.1165/rcmb.2018-0368OC</mixed-citation></ref><ref id="scirp.131783-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Saygin, D., Tabib, T., Bittar, H., et al. (2020) Transcriptional Profiling of Lung Cell Populations in Idiopathic Pulmonary Arterial Hypertension. Pulmonary Circulation, 10, 1-15. &lt;/br&gt;https://doi.org/10.1177/2045894020908782</mixed-citation></ref><ref id="scirp.131783-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Rahman, M., Rana, H., Peng, S., et al. (2021) Bioinformatics and Machine Learning Methodologies to Identify the Effects of Central Nervous System Disorders on Glioblastoma Progression. Briefings in Bioinformatics, 22, bbaa365. &lt;/br&gt;https://doi.org/10.1093/bib/bbaa365</mixed-citation></ref><ref id="scirp.131783-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Rahman, M., Peng, S., Hu, X., et al. (2020) A Network-Based Bioinformatics Approach to Identify Molecular Biomarkers for Type 2 Diabetes That Are Linked to the Progression of Neurological Diseases. International Journal of Environmental Research and Public Health, 17, Article No. 1035. &lt;/br&gt;https://doi.org/10.3390/ijerph17031035</mixed-citation></ref><ref id="scirp.131783-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Lu, X., Wang, L., Lin, X., et al. (2015) Genome-Wide Association Study in Chinese Identifies Novel Loci for Blood Pressure and Hypertension. Human Molecular Genetics, 24, 865-874. &lt;/br&gt;https://doi.org/10.1093/hmg/ddu478</mixed-citation></ref><ref id="scirp.131783-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">Stanton, A., Heydarpour, M., Williams, J., et al. (2023) CACNA1D Gene Polymorphisms Associate with Increased Blood Pressure and Salt Sensitivity of Blood Pressure in White Individuals. Hypertension (Dallas, Tex: 1979), 80, 2665-2673. &lt;/br&gt;https://doi.org/10.1161/HYPERTENSIONAHA.123.21229</mixed-citation></ref><ref id="scirp.131783-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">Pande, J., Mallhi, K., Sawh, A., et al. (2006) Aortic Smooth Muscle and Endothelial Plasma Membrane Ca2+ Pump Isoforms Are Inhibited Differently by the Extracellular Inhibitor Caloxin 1b1. American Journal of Physiology Cell Physiology, 290, C1341-C1349. &lt;/br&gt;https://doi.org/10.1152/ajpcell.00573.2005</mixed-citation></ref></ref-list></back></article>