# The Construction and Analysis of a ceRNA Network Related to Salt-Sensitivity Hypertensives

Authors: Xiu-Juan Liu; Hong-Lin Yin; Yan Li; Hao Hao; Yang Liu; Quan-Lin Zhao
Journal: BioMed Research International
 (2022)

---

## Abstract

Background. Salt-sensitivity hypertensives (SSH) are an independent risk factor for cardiovascular disease. However, the mechanism of SSH is not clear. This study is aimed at constructing a competing endogenous RNA (ceRNA) network related to SSH. Methods. Data sets were collected from the Gene Expression Omnibus database (GEO) to extract data on salt sensitivity RNA of patients with or without hypertensives in GSE135111. Firstly, we analyzed differentially expressed genes (DEGs, log2FC≥0.5 and P<0.05) and differentially expressed lncRNAs (DELs, log2FC ≥1 and P<0.05) between SSH and salt-sensitive normotension (SSN). Then, the gene ontology (GO), KEGG pathway enrichment analysis, and PPI network construction of DEGs were performed, and the hub genes in the PPI network by cytoHubba (12 methods) were screened out. Finally, a ceRNA network was constructed based on lncRNA-miRNA-mRNA pairs and hub genes. Results. 163 DEGs and 65 DELs were screened out. The GO and KEGG pathway analyses of DEGs were mainly enriched in metabolism (e.g., insulin secretion and cellular response to glucagon stimulus and peptidyl-tyrosine dephosphorylation,) and plasma membrane signaling (e.g., cell adhesion and chemical synaptic transmission and integral component of membrane). Additionally, a ceRNA network, including 1 mRNA (EGLN3), 2 miRNAs (hsa-miR-17-5p and hsa-miR-20b-5p), and 1 lncRNA (C1orf143) was successfully constructed. Conclusions. In conclusion, the proposed ceRNA network may help elucidate the regulatory mechanism by which lncRNAs function as ceRNAs and contribute to the pathogenesis of SSH. Importantly, candidate lncRNAs, miRNAs, and mRNAs can be further evaluated as a potential therapeutic targets for SSH.


---

## Body

## 1. Background

Hypertension, one of the most common diseases in humans, is one of the most well-known major risk factors for cardiovascular disease (CVD) and stroke. Accumulating study confirmed a strong relationship between sodium intake and blood pressure. Studies have confirmed that reducing the intake of salt in the diet can effectively reduce blood pressure. However, individuals have different responses to dietary salt. Some people will significantly increase their blood pressure after increasing their dietary salt intake, whereas another part of the population does not show significant changes in blood pressure. This phenomenon is called salt sensitivity [1]. Salt sensitivity is related to multiple factors, including the physiological environment, genetics, and demographic factors. Common demographic factors usually include sex, race, and age. However, the mechanism of salt-sensitivity hypertensives (SSH) is not clear. The mechanism of SSH has become the focus of many researchers.SSH can be defined as elevated blood pressure caused by relatively high salt intake [2]. The mechanism behind the rise in blood pressure caused by increased salt intake is complex. Current studies have found that it is related to an increase in blood volume caused by osmotic pressure, impaired endothelial function, and an imbalance in the regulation of nitric oxide and endothelin [3]. Furthermore, abnormal activation of the renin-angiotensin-aldosterone system, enhanced sympathetic nervous system, insulin resistance, and renal mechanism are also contribute [4]. Genetics, nutritional, and environmental factors are also involved in the development of salt-sensitive hypertension [5]. However, the BP of salt sensitivity differs among individuals. The blood pressure of salt-sensitive people can be divided into SSH and salt-sensitive normotension (SSN). Although both are sensitive to salt, their blood pressure is different. Its internal mechanism has yet to be studied while few people pay attention to their differences.Herein, we collect hypertension expression data for these three RNAs from the GEO database and construct a competing endogenous RNA (ceRNA) network related to SSH based on the base of DEGs (differentially expressed genes) and DELs (differentially expressed lncRNAs). This is the first study to investigate the ceRNA network related to SSH. After the analysis of DEGs and DELs, we found that some RNAs were tightly related to SSH. The DEGs and pathways predicted in our study may reveal the potential molecular mechanism of SSH. The workflow of this study is shown in Figure1.Figure 1
Workflow of this study. Notes: workflow of the study.


## 2. Methods

### 2.1. Data Collection and Processing

First, the GSE135111 microarray data, including 5 SSH blood samples and 5 SSN blood samples, was downloaded from the GEO database [2]. The microarray data were normalized by R software such as probe identification transformation and log2 transformation. Finally, the normalization data were used for further analysis.


### 2.2. Difference Analysis

In this study, the Limma package [6] was used to identify DEGs (log2FC≥0.5 and P value < 0.05), DELs (log2FC≥1 and P value < 0.05), and DECs (log2FC≥1 and P value <0.05). Then, DEGs and DELs were used in the next analysis.


### 2.3. GO and KEGG Pathway Enrichment Analyses of DEGs

Gene ontology (GO) is widely used to annotate genes, gene products, and sequences. KEGG is a comprehensive database for the biological interpretation of genome sequences and other high-throughput data. To represent the characteristics of DEGs, the enrichment analysis of the GO and KEGG pathways of DEGs was performed using a clusterProfiler package with the following criterion:P value <0.05.


### 2.4. Construction of the PPI Network and Screening of Hub Genes

To select hub genes related to SSH herein, the DEGs were mapped into PPIs (protein-protein interaction) by the STRING database, and a combined score of >0.4 was set as a threshold value. In addition, nodes with higher degrees of interaction from the PPI network were considered as hub nodes. As everyone knows, cytoHubba [7] is a tool for screening hub genes in the Cytoscape software [8]. Therefore, in this study, the hub gene modules (top 50 genes) were screened out by all 12 methods in Cytoscape software. Only the remaining overlapping genes in all 12 methods were selected as hub genes related to SSH.


### 2.5. Construction of lncRNA-miRNA-mRNA Pairs

As we know, lncRNA-miRNA pairs and miRNA-mRNA pairs can form lncRNA-miRNA-mRNA pairs. miRNA can bind to targeted mRNA to promote the degradation of mRNA, while lncRNA can bind to targeted miRNA to inhibit the degradation of mRNA. Herein, we used ggalluvial R package [9] to construct lncRNA-miRNA-mRNA pairs through miRcode (version 11; http://www.mircode.org/mircode/), miRDB (version 7.0; http://mirdb.org/), miRTarBase (http://mirtarbase.mbc.nctu.edu.tw/index.html), and TargetScan (version 7.2; http://targetscan.org/vert_72/) on the base of DEGs and DELs. miRcode provides “whole transcriptome” human microRNA target predictions based on the comprehensive GENCODE gene annotation, including >10,000 long noncoding RNA genes. Coding genes are also covered, including atypical regions such as 5′UTRs and CDS. miRDB is an online database for miRNA target prediction and functional annotations. All miRDB targets were predicted by a bioinformatics tool, MirTarget, which was developed by analyzing thousands of miRNA-target interactions from high-throughput sequencing experiments. miRTarBase is a database of experimentally validated microRNA targets. TargetScan predicts biological targets of miRNAs by searching for the presence of conserved 8mer, 7mer, and 6mer sites that match the seed region of each miRNA.Firstly, we predicted the lncRNA-miRNA pairs through the miRcode database on the base of DELs. The target genes for these miRNA signatures were then obtained using the miRDB, miRTarBase, and TargetScan databases. Genes present in all three databases were regarded as the target genes for these miRNAs. Comparing predicted target genes with DEGs, only the remaining overlapping genes and their interaction pairs were used for constructing the lncRNA-miRNA-mRNA pairs.


### 2.6. Construction of the ceRNA Network

In this study, the overlapping genes of hub genes and lncRNA-miRNA-mRNA pairs were used as potential key genes related to SSH. In addition, we construct a ceRNA network on hub genes and lncRNA-miRNA-mRNA pairs.


## 2.1. Data Collection and Processing

First, the GSE135111 microarray data, including 5 SSH blood samples and 5 SSN blood samples, was downloaded from the GEO database [2]. The microarray data were normalized by R software such as probe identification transformation and log2 transformation. Finally, the normalization data were used for further analysis.


## 2.2. Difference Analysis

In this study, the Limma package [6] was used to identify DEGs (log2FC≥0.5 and P value < 0.05), DELs (log2FC≥1 and P value < 0.05), and DECs (log2FC≥1 and P value <0.05). Then, DEGs and DELs were used in the next analysis.


## 2.3. GO and KEGG Pathway Enrichment Analyses of DEGs

Gene ontology (GO) is widely used to annotate genes, gene products, and sequences. KEGG is a comprehensive database for the biological interpretation of genome sequences and other high-throughput data. To represent the characteristics of DEGs, the enrichment analysis of the GO and KEGG pathways of DEGs was performed using a clusterProfiler package with the following criterion:P value <0.05.


## 2.4. Construction of the PPI Network and Screening of Hub Genes

To select hub genes related to SSH herein, the DEGs were mapped into PPIs (protein-protein interaction) by the STRING database, and a combined score of >0.4 was set as a threshold value. In addition, nodes with higher degrees of interaction from the PPI network were considered as hub nodes. As everyone knows, cytoHubba [7] is a tool for screening hub genes in the Cytoscape software [8]. Therefore, in this study, the hub gene modules (top 50 genes) were screened out by all 12 methods in Cytoscape software. Only the remaining overlapping genes in all 12 methods were selected as hub genes related to SSH.


## 2.5. Construction of lncRNA-miRNA-mRNA Pairs

As we know, lncRNA-miRNA pairs and miRNA-mRNA pairs can form lncRNA-miRNA-mRNA pairs. miRNA can bind to targeted mRNA to promote the degradation of mRNA, while lncRNA can bind to targeted miRNA to inhibit the degradation of mRNA. Herein, we used ggalluvial R package [9] to construct lncRNA-miRNA-mRNA pairs through miRcode (version 11; http://www.mircode.org/mircode/), miRDB (version 7.0; http://mirdb.org/), miRTarBase (http://mirtarbase.mbc.nctu.edu.tw/index.html), and TargetScan (version 7.2; http://targetscan.org/vert_72/) on the base of DEGs and DELs. miRcode provides “whole transcriptome” human microRNA target predictions based on the comprehensive GENCODE gene annotation, including >10,000 long noncoding RNA genes. Coding genes are also covered, including atypical regions such as 5′UTRs and CDS. miRDB is an online database for miRNA target prediction and functional annotations. All miRDB targets were predicted by a bioinformatics tool, MirTarget, which was developed by analyzing thousands of miRNA-target interactions from high-throughput sequencing experiments. miRTarBase is a database of experimentally validated microRNA targets. TargetScan predicts biological targets of miRNAs by searching for the presence of conserved 8mer, 7mer, and 6mer sites that match the seed region of each miRNA.Firstly, we predicted the lncRNA-miRNA pairs through the miRcode database on the base of DELs. The target genes for these miRNA signatures were then obtained using the miRDB, miRTarBase, and TargetScan databases. Genes present in all three databases were regarded as the target genes for these miRNAs. Comparing predicted target genes with DEGs, only the remaining overlapping genes and their interaction pairs were used for constructing the lncRNA-miRNA-mRNA pairs.


## 2.6. Construction of the ceRNA Network

In this study, the overlapping genes of hub genes and lncRNA-miRNA-mRNA pairs were used as potential key genes related to SSH. In addition, we construct a ceRNA network on hub genes and lncRNA-miRNA-mRNA pairs.


## 3. Results

### 3.1. Identification of DEGs and DELs

After preprocessing and data integration, 163 DEGs (86 upregulated and 77 downregulated), 65 DELs (34 upregulated and 31 downregulated), and 4 DECs (1 upregulated and 3 downregulated) were screened (Figures2(a) and 2(b)). Herein, the heat map of the lncRNAs, circRNAs, and mRNAs showed that the SSH clustered separately from the paired SSN (Figure 2(c)).Figure 2
Difference analysis of mRNAs and lncRNAs. Notes: (a) expression of mRNAs between two sets of samples (SSH vs. SSN). Blue represents the downregulated genes and red represents the upregulated genes, and the names of the 10 genes with the lowestP value have been indicated. (b) Expression of lncRNAs between two sets of samples (SSH vs. SSN). Blue represents the downregulated lncRNAs and red represents the upregulated lncRNAs, and the names of the 10 lncRNAs with the lowest P value have been indicated. (c) Hierarchical clustering heat map of DEGs (top 57), DELs (65), and DECs (4). Firebrick indicates that the relative expression of mRNAs/lncRNAs/circRNAs were upregulated, navy indicates that the relative expression of mRNAs/lncRNAs/circRNAs were downregulated.
(a)(b)(c)


### 3.2. Enrichment Analysis of DEGs

As shown in Table1 and Figure 3(a). GO and KEGG pathway enrichment analyses of DEGs were mainly enriched in metabolism (e.g., insulin secretion and cellular response to glucagon stimulus and peptidyl-tyrosine dephosphorylation) and plasma membrane signaling (e.g., cell adhesion and chemical synaptic transmission and integral component of membrane).Table 1
GO and KEGG pathway enrichment analyses of DEGs.
CategoryTermCount/10-log10 (P value)BPCell adhesion0.91.613392753Chemical synaptic transmission0.61.434237533Cellular response to glucagon stimulus0.31.428932795Peptidyl-tyrosine dephosphorylation0.41.39799905Adenylate cyclase-inhibiting G-protein coupled receptor signaling pathway0.31.30207393CCPlasma membrane4.41.791142931Intrinsic component of plasma membrane0.31.711353914Collagen trimer0.41.469118021Integral component of membrane5.11.467740938Transport vesicle membrane0.31.463466518Integral component of plasma membrane1.81.359780188MFTransmembrane receptor protein tyrosine phosphatase activity0.32.104180723Protein tyrosine phosphatase activity0.41.343534176PATHWAYGlutamatergic synapse0.51.867759405Insulin secretion0.41.500852037Tight junction0.41.475102357Figure 3
Enrichment analysis of DEGs and PPI network construction. Notes: (a) GO and KEGG pathway enrichment analyses of DEGs. (b) PPI network. Red represents upregulated DEGs. Navy blue represents downregulated DEGs, and yellow represents other genes associated with DEGs.
(a)(b)


### 3.3. PPI Network Construction and Screening of Hub Genes

In this study, these DEGs demonstrated significant interactions. A total of 59 nodes of the 163 DEGs were mapped in the PPI network (Figure3(b)). Four genes (7 ≥Degree ≥4), including ADCY2, TAS2R38, TNFSF11, and ADCY6, were located in the center of the PPI network. Additionally, 31 hub genes were selected, including EGLN3, TNFSF11, and DPPA4, and were screened (Table 2 and Figure 4(a)). Furthermore, the top 10 hub genes (radiality method) are shown in Figure 4(b).Table 2
The first 50 hub genes provided by MCC method.
MCCBetweennessBottleNeckClosenessClusteringCoefficientDegreeDMNCEcCentricityEPCMNCRadialityStressADCY6PSMC4TNFSF11ADCY6HTR1FADCY6HTR1FGCGADCY6ADCY6TNFSF11PSMC4ADCY2TNFSF11PSMC4TNFSF11NTPCRADCY2TAS2R38TNFSF11ADCY2ADCY2PSMC4TNFSF11TAS2R38CSF2ADCY2PSMC4ADCY2TAS2R38ADCY6PSMC4TAS2R38HTR1FCSF2CSF2HTR1FGTSE1GRIN2AADCY2TAS2R38TNFSF11ADCY2DPPA4HTR1FTAS2R38GTSE1ADCY6TNFSF11ADCY6KCNT1CSF2GCGKCNT1GCGCSF2GCGGCGEGLN3GTSE1KCNT1OBBP2PTPRBGTSE1ADCY6OBBP2NTPCREGLN3NTPCRNTPCROCSTAMPGCGOBBP2EGLN3CSF2TAS2R38REP15CSF2REP15GTSE1TNFSF11REP15SPIBOBBP2CSF2GCGADCY6GCGCLDN24GCGCLDN24OCSTAMPPSMC4CLDN24OBBP2NTSR1GCGNTSR1GCGOBBP2CLDN22HTR1FCLDN22SPIBCSF2CLDN22SDC2EGLN3PTPRBC7orf72GTSE1EGLN3NPPAPTPRBNPPAADCY6NTSR1NPPAKIF14ADCY2PSMC4KIF14DPPA4HTR1FTNNI3PSMC4TNNI3ADCY2PTPRRTNNI3GAS2L3TAS2R38GTSE1KCNT1OBBP2NTPCRETV5GTSE1ETV5NTSR1GTSE1ETV5C7orf72C7orf72DPPA4TAS2R38NTSR1KIF14FUT1DPPA4FUT1PTPRBAKAP5FUT1ADCY6KIF14FUT1PLCXD2EGLN3OCSTAMPRASIP1FUT1RASIP1OBBP2C4BPARASIP1ADCY2PLCXD2GRIN2AADCY2FUT1SPIBCYYR1GRIN2ACYYR1GAS2L3EGLN3CYYR1GCGKCNT1C7orf72GRIN2AC7orf72NTSR1SLC16A14C7orf72SLC16A14C7orf72OBBP2SLC16A14FCRL1GRIN2ANTSR1PTPRBSPRY1C7orf72ZNRF3NTSR1ZNRF3SDC2KCNT1ZNRF3ZNF577PTPRBSPRY1FUT1CMTM6PTPRRPTPRKSPRY1PTPRKKIF14PLCXD2PTPRKTAS2R38FUT1CMTM6DAOTAS2R38AKAP5GRIN2ACMTM6GRIN2APTPRRKIF14GRIN2AHTR1FDAOEGLN3SPRY1DAOSDC2FAM124AEGLN3FAM124AAKAP5PTPRBFAM124AARHGAP44SPRY1DAOCMTM6KIF14GAS2L3GTSF1DAOGTSF1HTR1FOCSTAMPGTSF1NTSR1CMTM6KIF14DPPA4PLCXD2PLCXD2KCNT1KIF14KCNT1TAS2R38SPIBKCNT1NTPCRDPPA4PLCXD2REP15REP15FCRL1TCL1APLCXD2TCL1APLCXD2C7orf72TCL1APTPRRREP15NTPCRCLDN24CLDN24ZNF577OBBP2NTPCROBBP2NTPCRGRIN2AOBBP2AKAP5CLDN24REP15CLDN22CLDN22C4BPAFCRL1REP15FCRL1GRIN2ASDC2FCRL1CCDC175CLDN22CLDN24NPPANPPAKCNT1PTPRRCLDN24PTPRRKCNT1FUT1PTPRRC4BPANPPACLDN22TNNI3TNNI3ARHGAP44GAS2L3CLDN22GAS2L3FCRL1SPRY1GAS2L3PLCXD2TNNI3NPPAETV5ETV5CCDC175CCDC175NPPACCDC175CCDC175GAS2L3CCDC175FAM124AETV5TNNI3RASIP1RASIP1PTPRBC7orf72TNNI3C7orf72ARHGAP44CMTM6C7orf72KCNT1RASIP1ETV5CYYR1CYYR1GRIN2AARHGAP44ETV5ARHGAP44ZNF577DPPA4ARHGAP44GRIN2ACYYR1RASIP1SLC16A14SLC16A14FUT1KCNMB3RASIP1KCNMB3CLDN24ZNF577KCNMB3FUT1SLC16A14CYYR1ZNRF3ZNRF3FAM124AC4BPACYYR1C4BPACLDN22DAOC4BPADAOZNRF3SLC16A14PTPRKPTPRKDAOAPOC1SLC16A14APOC1NPPAFCRL1APOC1KCNMB3PTPRKZNRF3FAM124AFAM124ASPRY1SLC17A4ZNRF3SLC17A4TNNI3ARHGAP44SLC17A4RASIP1FAM124APTPRKGTSF1GTSF1CMTM6ZNF577PTPRKZNF577CYYR1SIGLEC14ZNF577PTPRBGTSF1FAM124ATCL1ATCL1ADPPA4SDC2FAM124ASDC2SLC16A14KCNMB3SDC2HOGA1TCL1AGTSF1FCRL1FCRL1KCNMB3CSF2GTSF1CSF2ZNRF3FAM124ACSF2SPRY1FCRL1TCL1APTPRRPTPRRSIGLEC14AKAP5TCL1AAKAP5PTPRKCCDC175AKAP5CMTM6PTPRRFCRL1GAS2L3GAS2L3RASIP1NTSR1FCRL1NTSR1FAM124ARASIP1NTSR1SIGLEC14GAS2L3PTPRRCCDC175CCDC175HOGA1SPRY1PTPRRSPRY1GTSF1ETV5SPRY1REP15CCDC175GAS2L3ARHGAP44ARHGAP44REP15PTPRBGAS2L3PTPRBTCL1AREP15PTPRBETV5ARHGAP44CCDC175KCNMB3KCNMB3ETV5SIGLEC14CCDC175SIGLEC14C4BPAGTSF1SIGLEC14DPPA4KCNMB3ARHGAP44C4BPAC4BPAGTSF1CMTM6ARHGAP44CMTM6APOC1DNMT3BCMTM6GTSF1C4BPAKCNMB3APOC1APOC1TCL1ATNFSF11KCNMB3TNFSF11SLC17A4TCL1ATNFSF11TCL1AAPOC1C4BPASLC17A4SLC17A4DNMT3BPSMC4C4BPAPSMC4SPRY1HOGA1PSMC4DNMT3BSLC17A4APOC1ZNF577ZNF577CLDN24EGLN3APOC1EGLN3SIGLEC14AQP6EGLN3CLDN24ZNF577SLC17A4SDC2SDC2CLDN22GTSE1SLC17A4GTSE1CMTM6PVRL3GTSE1CLDN22SDC2ZNF577AKAP5AKAP5NPPAHOGA1ZNF577HOGA1RPL22L1ZNRF3HOGA1NPPAAKAP5SDC2HTR1FHTR1FTNNI3DAOSDC2DAOMRPL15PTPRKDAOTNNI3HTR1FAKAP5SIGLEC14SIGLEC14CYYR1RPL22L1AKAP5RPL22L1DNMT3BCLDN24RPL22L1CYYR1SIGLEC14Figure 4
Screening of hub genes. Notes: (a) upset plot. The horizontal bar chart shows the number of elements in each collection. The bar chart above represents the number of elements for each intersection and the point represent the relevant set involved. (b) Top 10 hub genes. The redder the color, the stronger their interaction.
(a)(b)


### 3.4. Construction of lncRNA-miRNA-mRNA Pairs

Herein, lncRNA-miRNA-mRNA pairs based on DEGs and DELs, including 2 mRNAs (EGLN3 and TNFRSF21), 2 miRNAs (hsa-miR-17-5p and hsa-miR-20b-5p), and 1 lncRNA (C1orf143) was constructed. As shown in Figure 5(a).Figure 5
ceRNA network construction. Notes: (a) Sankey diagram for the lncRNA-miRNA-mRNA pairs. Each rectangle represents a mRNA, miRNA, or lncRNA, and the connection degree of each gene is visualized based on the size of the rectangle. (b) Venn diagram of the hub gens and the lncRNA-miRNA-mRNA pairs. The overlapping genes in both represent the potential key genes related to SSH. (c) A ceRNA network. Red represents lncRNAs. Navy blue represents miRNAs. Violet represents mRNAs.
(a)(b)(c)


### 3.5. Construction of the ceRNA Network

In this study, an overlapping gene (EGLN3) between hub genes and lncRNA-miRNA-mRNA pairs was screened out and was identified as a potential key gene related to SSH (Figure 5(b)). Then, based on the potential key genes, a ceRNA network was successfully constructed, including 1 lncRNA (C1orf143), 2 miRNAs (hsa-miR-17-5p and hsa-miR-20b-5p), and 1 mRNA (EGLN3).


## 3.1. Identification of DEGs and DELs

After preprocessing and data integration, 163 DEGs (86 upregulated and 77 downregulated), 65 DELs (34 upregulated and 31 downregulated), and 4 DECs (1 upregulated and 3 downregulated) were screened (Figures2(a) and 2(b)). Herein, the heat map of the lncRNAs, circRNAs, and mRNAs showed that the SSH clustered separately from the paired SSN (Figure 2(c)).Figure 2
Difference analysis of mRNAs and lncRNAs. Notes: (a) expression of mRNAs between two sets of samples (SSH vs. SSN). Blue represents the downregulated genes and red represents the upregulated genes, and the names of the 10 genes with the lowestP value have been indicated. (b) Expression of lncRNAs between two sets of samples (SSH vs. SSN). Blue represents the downregulated lncRNAs and red represents the upregulated lncRNAs, and the names of the 10 lncRNAs with the lowest P value have been indicated. (c) Hierarchical clustering heat map of DEGs (top 57), DELs (65), and DECs (4). Firebrick indicates that the relative expression of mRNAs/lncRNAs/circRNAs were upregulated, navy indicates that the relative expression of mRNAs/lncRNAs/circRNAs were downregulated.
(a)(b)(c)


## 3.2. Enrichment Analysis of DEGs

As shown in Table1 and Figure 3(a). GO and KEGG pathway enrichment analyses of DEGs were mainly enriched in metabolism (e.g., insulin secretion and cellular response to glucagon stimulus and peptidyl-tyrosine dephosphorylation) and plasma membrane signaling (e.g., cell adhesion and chemical synaptic transmission and integral component of membrane).Table 1
GO and KEGG pathway enrichment analyses of DEGs.
CategoryTermCount/10-log10 (P value)BPCell adhesion0.91.613392753Chemical synaptic transmission0.61.434237533Cellular response to glucagon stimulus0.31.428932795Peptidyl-tyrosine dephosphorylation0.41.39799905Adenylate cyclase-inhibiting G-protein coupled receptor signaling pathway0.31.30207393CCPlasma membrane4.41.791142931Intrinsic component of plasma membrane0.31.711353914Collagen trimer0.41.469118021Integral component of membrane5.11.467740938Transport vesicle membrane0.31.463466518Integral component of plasma membrane1.81.359780188MFTransmembrane receptor protein tyrosine phosphatase activity0.32.104180723Protein tyrosine phosphatase activity0.41.343534176PATHWAYGlutamatergic synapse0.51.867759405Insulin secretion0.41.500852037Tight junction0.41.475102357Figure 3
Enrichment analysis of DEGs and PPI network construction. Notes: (a) GO and KEGG pathway enrichment analyses of DEGs. (b) PPI network. Red represents upregulated DEGs. Navy blue represents downregulated DEGs, and yellow represents other genes associated with DEGs.
(a)(b)


## 3.3. PPI Network Construction and Screening of Hub Genes

In this study, these DEGs demonstrated significant interactions. A total of 59 nodes of the 163 DEGs were mapped in the PPI network (Figure3(b)). Four genes (7 ≥Degree ≥4), including ADCY2, TAS2R38, TNFSF11, and ADCY6, were located in the center of the PPI network. Additionally, 31 hub genes were selected, including EGLN3, TNFSF11, and DPPA4, and were screened (Table 2 and Figure 4(a)). Furthermore, the top 10 hub genes (radiality method) are shown in Figure 4(b).Table 2
The first 50 hub genes provided by MCC method.
MCCBetweennessBottleNeckClosenessClusteringCoefficientDegreeDMNCEcCentricityEPCMNCRadialityStressADCY6PSMC4TNFSF11ADCY6HTR1FADCY6HTR1FGCGADCY6ADCY6TNFSF11PSMC4ADCY2TNFSF11PSMC4TNFSF11NTPCRADCY2TAS2R38TNFSF11ADCY2ADCY2PSMC4TNFSF11TAS2R38CSF2ADCY2PSMC4ADCY2TAS2R38ADCY6PSMC4TAS2R38HTR1FCSF2CSF2HTR1FGTSE1GRIN2AADCY2TAS2R38TNFSF11ADCY2DPPA4HTR1FTAS2R38GTSE1ADCY6TNFSF11ADCY6KCNT1CSF2GCGKCNT1GCGCSF2GCGGCGEGLN3GTSE1KCNT1OBBP2PTPRBGTSE1ADCY6OBBP2NTPCREGLN3NTPCRNTPCROCSTAMPGCGOBBP2EGLN3CSF2TAS2R38REP15CSF2REP15GTSE1TNFSF11REP15SPIBOBBP2CSF2GCGADCY6GCGCLDN24GCGCLDN24OCSTAMPPSMC4CLDN24OBBP2NTSR1GCGNTSR1GCGOBBP2CLDN22HTR1FCLDN22SPIBCSF2CLDN22SDC2EGLN3PTPRBC7orf72GTSE1EGLN3NPPAPTPRBNPPAADCY6NTSR1NPPAKIF14ADCY2PSMC4KIF14DPPA4HTR1FTNNI3PSMC4TNNI3ADCY2PTPRRTNNI3GAS2L3TAS2R38GTSE1KCNT1OBBP2NTPCRETV5GTSE1ETV5NTSR1GTSE1ETV5C7orf72C7orf72DPPA4TAS2R38NTSR1KIF14FUT1DPPA4FUT1PTPRBAKAP5FUT1ADCY6KIF14FUT1PLCXD2EGLN3OCSTAMPRASIP1FUT1RASIP1OBBP2C4BPARASIP1ADCY2PLCXD2GRIN2AADCY2FUT1SPIBCYYR1GRIN2ACYYR1GAS2L3EGLN3CYYR1GCGKCNT1C7orf72GRIN2AC7orf72NTSR1SLC16A14C7orf72SLC16A14C7orf72OBBP2SLC16A14FCRL1GRIN2ANTSR1PTPRBSPRY1C7orf72ZNRF3NTSR1ZNRF3SDC2KCNT1ZNRF3ZNF577PTPRBSPRY1FUT1CMTM6PTPRRPTPRKSPRY1PTPRKKIF14PLCXD2PTPRKTAS2R38FUT1CMTM6DAOTAS2R38AKAP5GRIN2ACMTM6GRIN2APTPRRKIF14GRIN2AHTR1FDAOEGLN3SPRY1DAOSDC2FAM124AEGLN3FAM124AAKAP5PTPRBFAM124AARHGAP44SPRY1DAOCMTM6KIF14GAS2L3GTSF1DAOGTSF1HTR1FOCSTAMPGTSF1NTSR1CMTM6KIF14DPPA4PLCXD2PLCXD2KCNT1KIF14KCNT1TAS2R38SPIBKCNT1NTPCRDPPA4PLCXD2REP15REP15FCRL1TCL1APLCXD2TCL1APLCXD2C7orf72TCL1APTPRRREP15NTPCRCLDN24CLDN24ZNF577OBBP2NTPCROBBP2NTPCRGRIN2AOBBP2AKAP5CLDN24REP15CLDN22CLDN22C4BPAFCRL1REP15FCRL1GRIN2ASDC2FCRL1CCDC175CLDN22CLDN24NPPANPPAKCNT1PTPRRCLDN24PTPRRKCNT1FUT1PTPRRC4BPANPPACLDN22TNNI3TNNI3ARHGAP44GAS2L3CLDN22GAS2L3FCRL1SPRY1GAS2L3PLCXD2TNNI3NPPAETV5ETV5CCDC175CCDC175NPPACCDC175CCDC175GAS2L3CCDC175FAM124AETV5TNNI3RASIP1RASIP1PTPRBC7orf72TNNI3C7orf72ARHGAP44CMTM6C7orf72KCNT1RASIP1ETV5CYYR1CYYR1GRIN2AARHGAP44ETV5ARHGAP44ZNF577DPPA4ARHGAP44GRIN2ACYYR1RASIP1SLC16A14SLC16A14FUT1KCNMB3RASIP1KCNMB3CLDN24ZNF577KCNMB3FUT1SLC16A14CYYR1ZNRF3ZNRF3FAM124AC4BPACYYR1C4BPACLDN22DAOC4BPADAOZNRF3SLC16A14PTPRKPTPRKDAOAPOC1SLC16A14APOC1NPPAFCRL1APOC1KCNMB3PTPRKZNRF3FAM124AFAM124ASPRY1SLC17A4ZNRF3SLC17A4TNNI3ARHGAP44SLC17A4RASIP1FAM124APTPRKGTSF1GTSF1CMTM6ZNF577PTPRKZNF577CYYR1SIGLEC14ZNF577PTPRBGTSF1FAM124ATCL1ATCL1ADPPA4SDC2FAM124ASDC2SLC16A14KCNMB3SDC2HOGA1TCL1AGTSF1FCRL1FCRL1KCNMB3CSF2GTSF1CSF2ZNRF3FAM124ACSF2SPRY1FCRL1TCL1APTPRRPTPRRSIGLEC14AKAP5TCL1AAKAP5PTPRKCCDC175AKAP5CMTM6PTPRRFCRL1GAS2L3GAS2L3RASIP1NTSR1FCRL1NTSR1FAM124ARASIP1NTSR1SIGLEC14GAS2L3PTPRRCCDC175CCDC175HOGA1SPRY1PTPRRSPRY1GTSF1ETV5SPRY1REP15CCDC175GAS2L3ARHGAP44ARHGAP44REP15PTPRBGAS2L3PTPRBTCL1AREP15PTPRBETV5ARHGAP44CCDC175KCNMB3KCNMB3ETV5SIGLEC14CCDC175SIGLEC14C4BPAGTSF1SIGLEC14DPPA4KCNMB3ARHGAP44C4BPAC4BPAGTSF1CMTM6ARHGAP44CMTM6APOC1DNMT3BCMTM6GTSF1C4BPAKCNMB3APOC1APOC1TCL1ATNFSF11KCNMB3TNFSF11SLC17A4TCL1ATNFSF11TCL1AAPOC1C4BPASLC17A4SLC17A4DNMT3BPSMC4C4BPAPSMC4SPRY1HOGA1PSMC4DNMT3BSLC17A4APOC1ZNF577ZNF577CLDN24EGLN3APOC1EGLN3SIGLEC14AQP6EGLN3CLDN24ZNF577SLC17A4SDC2SDC2CLDN22GTSE1SLC17A4GTSE1CMTM6PVRL3GTSE1CLDN22SDC2ZNF577AKAP5AKAP5NPPAHOGA1ZNF577HOGA1RPL22L1ZNRF3HOGA1NPPAAKAP5SDC2HTR1FHTR1FTNNI3DAOSDC2DAOMRPL15PTPRKDAOTNNI3HTR1FAKAP5SIGLEC14SIGLEC14CYYR1RPL22L1AKAP5RPL22L1DNMT3BCLDN24RPL22L1CYYR1SIGLEC14Figure 4
Screening of hub genes. Notes: (a) upset plot. The horizontal bar chart shows the number of elements in each collection. The bar chart above represents the number of elements for each intersection and the point represent the relevant set involved. (b) Top 10 hub genes. The redder the color, the stronger their interaction.
(a)(b)


## 3.4. Construction of lncRNA-miRNA-mRNA Pairs

Herein, lncRNA-miRNA-mRNA pairs based on DEGs and DELs, including 2 mRNAs (EGLN3 and TNFRSF21), 2 miRNAs (hsa-miR-17-5p and hsa-miR-20b-5p), and 1 lncRNA (C1orf143) was constructed. As shown in Figure 5(a).Figure 5
ceRNA network construction. Notes: (a) Sankey diagram for the lncRNA-miRNA-mRNA pairs. Each rectangle represents a mRNA, miRNA, or lncRNA, and the connection degree of each gene is visualized based on the size of the rectangle. (b) Venn diagram of the hub gens and the lncRNA-miRNA-mRNA pairs. The overlapping genes in both represent the potential key genes related to SSH. (c) A ceRNA network. Red represents lncRNAs. Navy blue represents miRNAs. Violet represents mRNAs.
(a)(b)(c)


## 3.5. Construction of the ceRNA Network

In this study, an overlapping gene (EGLN3) between hub genes and lncRNA-miRNA-mRNA pairs was screened out and was identified as a potential key gene related to SSH (Figure 5(b)). Then, based on the potential key genes, a ceRNA network was successfully constructed, including 1 lncRNA (C1orf143), 2 miRNAs (hsa-miR-17-5p and hsa-miR-20b-5p), and 1 mRNA (EGLN3).


## 4. Discussion

It is generally known that salt sensitivity of blood pressure (SSBP) is an independent risk factor for cardiovascular disease. Although they are both sensitive to salt, their blood pressure is different. The blood pressure of salt-sensitive people can be divided into SSH and salt-sensitive normotension (SSN) [10]. As for SSH, salt sensitivity is related to multiple factors, including physiological environment, genetics, and demographic factors. Common demographic factors generally include gender, race, and age [11]. However, the pathogenic mechanisms of SSBP are still uncertain. Therefore, the mechanism of salt-sensitivity hypertensives (SSH) has become the focus of many researchers.Nowadays, it has become a popular method for evaluating DEGs through gene expression analysis to explore the causes of diseases [12]. These different genome-wide expression profiling techniques make them more valuable on account of complementary results. In this study, we obtained the analysis of the expression of mRNAs and lncRNAs of whole blood samples from patients with SSH and SSN with active disease and to explored potential RNAs related to SSH by the bioinformatic analysis. Then, we selected 163 DEGs and 65 DELs as our subsequent research object. Furthermore, an enrichment analysis of the GO, KEGG pathway, and the construction of DEG PPI networks were performed, and we found that these 163 genes may participate in the process of SSH through metabolism (e.g., insulin secretion and cell response to glucagon stimulus and tyrosine dephosphorylation) and plasma membrane signaling (e.g., cell adhesion and chemical synaptic transmission and integral component of membrane). This is consistent with previous studies.In metabolism, such as previous studies have demonstrated that single nucleotide polymorphisms (SNPs) of the sodium-bicarbonate cotransporter gene (SLC4A5) are associated with hypertension [13]. Furusho et al. found that mutations of with-no-lysine kinase 1 (WNK1) could lead to abnormally increased salt reabsorption and salt-sensitive hypertension [14]. Based on the epidemiological baseline survey, fasting blood glucose was found to be an independent and dose-dependent related factor of blood pressure salt sensitivity [15]. At the same time, it is well known that the ion channel is closely related to SSH. It has been found that there are many pathways mediated phosphorylation of ion channels in the body, including PKA-mediated phosphorylation, PKC-mediated phosphorylation, PI-3 K/PKB and PI-3 K/SGK3 pathway phosphorylation, MAPKs pathway-mediated phosphorylation, and Src channel regulation. However, there are still some signaling pathways that have not yet been reported in the results of our study such as tyrosine dephosphorylation. Tyrosine dephosphorylation is regulated by protein tyrosine phosphatase (PTP). Studies have found that the PTP gene mainly expresses 112 PTP in the human genome. It is not only an indispensable specific regulator in cell signal transduction but also a key signal molecule. It plays a key role in the regulation of different physiological events and is related to many diseases such as metabolism, cardiovascular diseases, cancer, and autoimmune diseases [15].The composition of its catalytic sites determines the sensitivity of dephosphorylation [16]. The same PTP can also be divided into two categories: the receptor tyrosine phosphatase and the nonreceptor tyrosine phosphatase. It is worth noting that the K+ channel is closely related to SSH. At the same time, previous articles have reported that K+ channels contain multiple tyrosine residues [17]. Some protein enzymes can act on channel protein synthesis and transport to change the phosphorylation state of channel proteins and regulate channel functions. Meanwhile, our results suggest that PTP may be involved in the occurrence of SSH. This gives us great hints. However, whether PTP participates in SSH through the dephosphorylation of K+ channel protein remains to be verified by experiments.In plasma membrane signaling, at present, many types of research are thorough. Such as G-protein coupled receptor signaling pathway and insulin secretion [18]. Note that we constructed an SSH-related ceRNA network, including 1 mRNA (EGLN3), 2 miRNAs (hsa-miR-17-5p and hsa-miR-20b-5p), and 1 lncRNA (C1orf143). Among the ceRNA regulatory network, 1 lncRNA (C1orf143) is associated with two miRNAs (hsa-miR-17-5p and hsa-miR-20b-5p), and we proposed that possible competition for C1orf143 binding to hsa-miR-17-5p and hsa-miR-20b-5p influences the downstream regulation of EGLN3. However, this requires experimental verification.EGLN3 (also named Prolyl-hydroxylase 3, PHD3, HPH1, and SM-20) belongs to the EGLN family of prolyl hydroxylases and can catalyze hydroxylation. Studies have shown that EGLN3 is involved in the metabolism and angiogenesis of oxygen from tumor cells, which in turn affects tumor cell proliferation [19]. However, there is no report about EGLN3 participating in SSH. It has been reported that EGLN3 is regulated by miRNA to play a role in different biological processes. However, so far, there has been no research on the EGLN3 ceRNA network in SSH. What is more, lncRNAs play a significant role in the development of the disease according to many studies. In our study, we found that C1orf143 could lncRNAs play important biological roles by regulating gene expression (a ceRNA network) in SSH. It is worth noting that there is no research report on C1orf143. This is a new gene. Our research results indicate that C1orf143 may play an important role in SSH. This gives us great hints. However, whether C1orf143 is related to SSH still needs further experimental verification.On the whole, based on comprehensive bioinformatics analysis of multiple cohort datasets of SSH and SSN patients, 163 DEGs were identified. The enrichment analysis of DEGs involving related molecules or pathways may deepen our understanding of SSH. Additionally, the SSH-related ceRNA network, including 1 mRNA (EGLN3), 2 miRNAs (hsa-miR-17-5p and hsa-miR-20b-5p), and 1 lncRNA (C1orf143) was successfully constructed. However, some limitations still existed in this research. The main method of our research was bioinformatics technology, a useful tool to understand interactions, pathways, and networks.In conclusion, these findings could enrich DEG expression profile between SSH and SSN and provide novel information on the occurrence of SSH. Other scientific researchers are expected to verify this at the genetic level.



---

*Source: 8258351-2022-10-14.xml*