Identification of biomarkers and pathways in hypertensive nephropathy based on the ceRNA regulatory network

doi:10.21203/rs.2.16949/v1

Download PDF

Research article

Identification of biomarkers and pathways in hypertensive nephropathy based on the ceRNA regulatory network

https://doi.org/10.21203/rs.2.16949/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 11 Nov, 2020

Read the published version in BMC Nephrology →

You are reading this older preprint version

Read the latest preprint version →

Background Hypertensive nephropathy (HTN) is a kind of renal injury caused by chronic hypertension, which seriously affect people’s life. The purpose of this study was to identify the potential biomarkers of HTN and understand its possible mechanisms.

Methods The dataset numbered GSE28260 related to hypertensive and normotensive was downloaded from NCBI Gene Expression Omnibus. Then, the differentially expressed RNAs (DERs) were screened using R limma package, and functional analyses of DE-mRNA were performed by DAVID. Afterwards, a ceRNA network was established and KEGG pathway was analyzed based on the Gene Set Enrichment Analysis (GSEA) database. Finally, a ceRNA regulatory network directly associated with HTN was proposed.

Results A total of 947 DERs were identified, including 900 DE-mRNAs, 20 DE-lncRNAs and 27 DE-miRNAs. Based on these DE-mRNAs, they were involved in biological processes such as fatty acid beta-oxidation, IRE1-mediated unfolded protein response, and transmembrane transport, and many KEGG pathways like glycine, serine and threonine metabolism, carbon metabolism. Subsequently, lncRNAs KCTD21-AS1, LINC00470 and SNHG14 were found to be hub nodes in the ceRNA regulatory network. KEGG analysis showed that insulin signaling pathway, glycine, serine and threonine metabolism, pathways in cancer, lysosome, and apoptosis was associated with hypertensive. Finally, insulin signaling pathway was screened to directly associate with HTN and was regulated by mRNAs PPP1R3C, PPKAR2B and AKT3, miRNA has-miR-107, and lncRNAs SNHG14, TUG1, ZNF252P-AS1 and MIR503HG.

Conclusions Insulin signaling pathway was directly associated with HTN, and miRNA has-miR-107 and lncRNAs SNHG14, TUG1, ZNF252P-AS1 and MIR503HG were the biomarkers of HTN. These results would improve our understanding of the occurrence and development of HTN.

Urology & Nephrology

Hypertensive nephropathy

differentially expressed RNAs

ceRNA network directly associated with HTN

functional analysis

Hypertension, which is influenced by genetic and environmental factors, is a risk factor of cardiovascular and cerebrovascular disease. The kidney is not only an important organ causing hypertension, but also the target organ of hypertensive damage [1]. Hypertensive nephropathy (HTN) is a kind of renal injury caused by chronic hypertension [2]. When hypertension lasts for 5–10 years, renal arteriolosclerosis, thickening of tube walls, and narrowing of bureaucratic cavities can occur, which will trigger substantial ischemic injury to the kidney [3]. Although the antihypertensive agent such as cilnidipine and avosentan are commonly used for the treatment of HNT, the effect of clinical treatment is still not ideal [4]. Due to the increasing occurrence and mortality rate of HTN, molecular mechanisms of hypertensive effect on kidney are urgent to be revealed, which is helpful for the early diagnosis and targeted control therapy of HTN.

Previously, a number of biomarkers have been identified to reveal the potential mechanisms of HTN, such as β2-microglobulin (β2-MG), transforming growth factor-β (TGF-β) and periostin [5, 6]. β2-MG, a small globulin, was a good indicator for the diagnosis of early HTN [7]. The high level of β2-MG in serum indicated the damaged filtration function and increased filtration load of glomerular [5]. TGF-β was also reported to be associated with HTN and its high level affected the progression of HTN [8]. Furthermore, periostin, related to proteinuria, plasma creatinine and renal flow, has been identified as the key markers of progression and regression in HTN [6]. Although some protein biomarkers of HTN were identified, there were little researches on gene markers of HTN.

Competing endogenous RNAs (ceRNAs) can regulate gene expression by competing for the binding microRNAs (miRNAs). A ceRNA regulatory network displayed a novel level of long non-coding RNA (lncRNA), microRNA (miRNA) and Mrna [9]. Increasing evidence showed that ceRNAs play important roles in many biological processes, and its imbalance may lead to the initiation and development of human disease [10–12]. By constructed ceRNA regulatory network of gastric cancer, Liu et al. [13] indicated that deleted in lymphocytic leukemia 2 (DLEU2) and DX11 antisense RNA 1 (DDX11-AS1) were hub lncRNAs in gastric cancer and act as potential ceRNAs to sponge miRNA. LncRNA Arid2-IR was one of the most highly up-regulated lncRNAs in mouse model of nonimmune-related kidney diseases and might play a functional role in renal inflammation [14]. In addition, Chen et al. [2] showed that has-miR–335–5p and has-miR–26b–5p, which modulated the majority of differentially expressed RNAs (DERs), were the two most outstanding miRNAs in HTN, by analyzing the miRNA-DERs regulatory network. Although some studies have indicated potential biomarkers associated with HTN, single biomarkers have limitations to explore the mechanism of hypertension on kidney.

In this study, the DERs were screened between hypertensive and normotensive, and then a ceRNA regulatory network was constructed comprised of lncRNA, miRNA and mRNA. Subsequently, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of genes existed in the ceRNA regulatory network were performed based on the Gene Set Enrichment Analysis (GSEA). Finally, based on the results of KEGG, a ceRNA regulatory network directly associated with HTN was proposed. These results would improve our understanding of the occurrence and development of HTN.

Data pre-processing and annotation of RNA

We downloaded the dataset numbered GSE28260 from the NCBI Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/) on July 26, 2019, which included two cubes: gene expression profile GSE28345 and miRNA expression profile GSE28283 [15]. Each cube contained 8 samples, including 4 samples of normotensive and 4 samples of hypertensive. The renal cortical tissue expression profiles in normotensive and hypertensive were originated from Affymetrix Human Gene 1.0 ST Array and Agilent–021827 Human miRNA Microarray (V3), respectively.

The data of GSE28345 from Affymetrix platform were pre-treated with expression value background correction and normalization by R 3.4.1 oligo package (Version 3.6, http://www.bioconductor.org/packages/release/bioc/html/oligo.html) [16], including conversion of original data format, background correction (MAS method), supplementation of missing value by median method, and data standardization by quantiles method. Afterwards, based on that the one gene can correspond to multiple probes, the data of GSE28283 from Agilent platform was first commented and then calculated the average value of these multiple expression. Subsequently, the data was carried on a logarithm analysis by R 3.4.1 limma package (Version 3.34.0, https://bioconductor.org/packages/release/bioc/html/limma.html) [17].

Selection of DERs and functional analysis

The lncRNAs, miRNAs and mRNAs detected in expression profile were annotated by HUGO Gene Nomenclature Committee (HGNC) database which including 4409 lncRNAs, 1914 miRNAs and 19221 protein coding gene (http://www.genenames.org/) [18]. Then the DERs between hypertensive and normotensive groups were screened using R 3.4.1 limma package Version 3.34.0 [17]. False discovery rate (FDR) value < 0.05 and |log₂ fold change (FC)|＞0.5 were set as threshold to identify DERs.

Based on the screened DERs, the pheatmap (Version 1.0.8, https://cran.r-project.org/package = pheatmap) [19] of R 3.4.1 was utilized for bidirectional hierarchical clustering [20] according to the Euclidean distance [21] and was displayed with heat map. Afterwards, Gene ontology (GO) terms and KEGG pathway enrichment analysis of differentially expressed mRNAs (DE-mRNAs) were performed by DAVID (Version 6.8, https://david.ncifcrf.gov/) [22, 23]. The cut-off of the significantly enriched GO terms and pathways was P value < 0.05.

Construction of ceRNA regulatory network and KEGG analysis

We first downloaded the connection pairs of lncRNA-miRNA from DIANA-LncBase (Version 2, http://carolina.imis.athena-innovation.gr/diana_tools/web/index.php) database [24]. Then the differentially expressed regulatory relationship between lncRNA and miRNA was selected according to the miRNA target gene score (miTG-score) > 0.6. The lncRNA-miRNA network was established with the connection pairs of the opposite expression and was visualized by Cytoscape 3.6.1 (http://www.cytoscape.org/) [25]. Afterwards, starBase (Version 2.0, http://starbase.sysu.edu.cn/) database was used to predict the targeted genes regulated by miRNA connected with lncRNA, and miRNA-mRNA network was constructed and visualized. The starBase database provide the predictive information of the targeted genes by five algorithms (targetScan, picTar, RNA22, PITA and miRanda).

Using the above obtained lncRNA-miRNA pairs and mRNA-miRNA pairs, a ceRNA regulatory network was proposed by combining lncRNA-miRNA-mRNA pairs and visualized by Cytoscape 3.6.1. Finally, the KEGG pathway-enrichment analysis of genes included in the ceRNA regulatory network was performed based on the GSEA (http://software.broadinstitute.org/gsea/index.jsp) [26]. GSEA software is used to evaluate the distribution trend of a predefined set of genes in a gene list that is sequenced by phenotypic relevance, thus determining their contribution to the phenotype.

Establishment of ceRNA regulatory network directly associated with HTN

“Hypertensive renal” as a keyword was searched in Comparative Toxicogenomics Database 2019 update (http://ctd.mdibl.org/) database [27]. The KEGG pathways related to HTN were found and compared with the above KEGG pathways in the ceRNA regulatory network. RNAs associated with HTN were identified and a ceRNA regulatory network directly associated with HTN was constructed with these RNAs.

Screening of DERs

After annotation by the HGNC database, a total of 459 lncRNAs, 851 miRNAs and 17962 mRNAs were identified. The DERs between hypertensive and normotensive were screened using limma package, and 947 DERs were filtered based on the thresholds of FDR < 0.05 and |log₂FC| > 1 (Fig. 1A).. These DERs contained 900 DE-mRNAs (352 down-regulated mRNAs and 548 up-regulated mRNAs), 20 DE-lncRNAs (16 down-regulated lncRNAs and 4 up-regulated lncRNAs) and 27 DE-miRNAs (12 down-regulated miRNAs and 15 up-regulated miRNAs). Based on the expression level of selected DERs, the bidirectional hierarchical clustering was performed. The results showed that the expression value of DERs could well differentiate the hypertensive and normotensive and indicated that the screened DERs were characteristically expressive (Fig. 1B)..

GO terms and KEGG pathways of screened DE-mRNAs

The functional analyses, including GO terms and KEGG pathways, were performed on the above 900 DE-mRNAs, and a total of 21 biological process and 18 KEGG pathways were obtained (Fig. 2).. As shown in Fig. 2A, in the involved biological process, the significant enrichment terms of DE-mRNAs were associated with fatty acid beta-oxidation, proteolysis involved in cellular protein catabolic process, IRE1-mediated unfolded protein response, and transmembrane transport. In addition, the DE-mRNAs also played an important role in glycine, serine and threonine metabolism, carbon metabolism, peroxisome, lysosome, metabolic pathways, biosynthesis of antibiotics, arginine and proline metabolism and citrate cycle (TCA cycle) based on the analysis of KEGG (Fig. 2B)..

Construction of ceRNA regulatory network

Before establishing a ceRNA regulatory network, co-expression networks of lncRNA-miRNA and miRNA-mRNA were constructed. The lncRNA-miRNA co-expression network showed that a total of 19 connection pairs were obtained, including 18 nodes, 6 miRNAs (3 down-regulated miRNAs and 3 up-regulated miRNAs) and 12 lncRNAs (8 down-regulated lncRNAs and 4 up-regulated lncRNAs) (Fig. 3A).. Afterwards, the results of miRNA-mRNA network displayed that there were 232 nodes in the network and the network contained 6 miRNAs (3 down-regulated miRNAs and 3 up-regulated miRNAs) and 226 mRNAs (68 down-regulated mRNAs and 158 up-regulated mRNAs) (Fig. 3B)..

Based on the results of lncRNA-miRNA and miRNA-mRNA, a ceRNA regulatory network consisting of lncRNA-miRNA-mRNA was proposed (Fig. 3C).. In the ceRNA regulatory network, 244 nodes were obtained. In addition, 12 lncRNAs (8 down-regulation and 4 up-regulation), 6 miRNAs (3 down-regulation and 3 up-regulation) and 226 mRNAs (68 down-regulation and 158 up-regulation) were identified. Among them, the 8 down-regulated lncRNAs contained LINC00184, small nucleolar RNA host gene 5 (SNHG5),, Jpx transcript Xist activator (JPX),, urothelial cancer associated 1 (UCA1),, small nucleolar RNA host gene 14 (SNHG14),, MIR503 host gene (MIR503HG),, taurine up-regulated 1 (TUG1),, ZNF252P antisense RNA 1 (ZNF252P-AS1),, and the 4 up-regulated lncRNAs included family with sequence similarity 95 member B1 (FAM95B1),, LINC00470, KCTD21 antisense RNA 1 (KCTD21-AS1),, small nucleolar RNA host gene 22 (SNHG22).. The 6 miRNAs included hsa-miR–548c–5p, hsa-miR–369–3p, has-miR–1290, has-miR–532–5p, has-miR–28–5p and has-miR–107. Furthermore, in this ceRNA regulatory network, lncRNAs KCTD21-AS1, LINC00470 and SNHG14 were hub nodes, which targeted more miRNAs and mRNAs.

KEGG enrichment-pathway analysis of mRNA in the ceRNA regulatory network based on GSEA

After the establishment of ceRNA regulatory network, the KEGG pathways were analyzed by GSEA. In the GSEA results, there are three key statistical values, including enrichment score (ES), normalized enrichment score (NES) and nominal P value. Based on the above three key statistical values, five significant KEGG signaling pathways (insulin signaling pathway, glycine, serine and threonine metabolism, pathways in cancer, lysosome, and apoptosis) were identified with the cut-off of P < 0.05 (Table 1).. Obviously, mRNA protein kinase cAMP-dependent type II regulatory subunit beta (PRKAR2B) and AKT serine/threonine kinase 3 (AKT3) were both related to the pathways of insulin signaling pathway and apoptosis. AKT3, KIT proto-oncogene, receptor tyrosine kinase (KIT) and retinoic acid receptor beta (RARB) were also involved in pathways in cancer. Furthermore, mRNA protein phosphatase 1 regulatory subunit 3C (PPP1R3C) was also essential for insulin signaling pathway.

Establishment of ceRNA regulatory network directly associated with HTN

In order to better understand the potential biomarkers directly associated with HTN, the KEGG pathways related to HTN were obtained according to the CTD database. After compared with the above KEGG pathways of mRNAs, one pathway directly associated with HTN (insulin signaling pathway) was screened and contained three important mRNAs (PPP1R3C, PPKAR2B and AKT3).. Based on the three genes, a ceRNA regulatory network directly associated with HTN was established. As shown in Fig. 4, insulin signaling pathway was regulated by mRNAs PPP1R3C, PPKAR2B and AKT3, miRNA has-miR–107, and lncRNAs SNHG14, TUG1, ZNF252P-AS1 and MIR503HG.

Recently, HTN is an increasingly common kidney disease in the patients with hypertension, which seriously endangers people’s health. However, the possible mechanism of hypertension affecting HTN remains unclear. In this study, we identified 947 DERs, including 900 DE-mRNAs, 20 DE-lncRNAs and 27 DE-miRNAs. Based on the selected DE-mRNAs, it was found that these mRNAs were involved in biological processes of fatty acid beta-oxidation, proteolysis involved in cellular protein catabolic process, IRE1-mediated unfolded protein response, and transmembrane transport, and many KEGG pathways like glycine, serine and threonine metabolism, carbon metabolism, metabolic pathways, biosynthesis of antibiotics, and citrate cycle (TCA cycle). Afterwards, a ceRNA regulatory network was established with 12 lncRNAs, 6 miRNAs and 226 mRNAs. It is clear that lncRNAs KCTD21-AS1, LINC00470 and SNHG14 were hub nodes in the ceRNA regulatory network, which targeted more miRNAs and mRNAs. Then KEGG analysis showed that mRNAs included in this network played important roles in insulin signaling pathway, glycine, serine and threonine metabolism, pathways in cancer, lysosome, and apoptosis. Finally, a ceRNA regulatory network directly associated with HTN was proposed. Insulin signaling pathway was screened to directly associate with HTN and was modulated by mRNAs PPP1R3C, PPKAR2B and AKT3, miRNA has-miR–107, and lncRNAs SNHG14, TUG1, ZNF252P-AS1 and MIR503HG.

At the beginning, 947 DERs were screened out and were found to play roles in many biological processes and signaling pathways. In this study, the screened DERs between hypertensive and normotensive were significantly associated with fatty acid beta-oxidation, IRE1-mediated unfolded protein response, and transmembrane transport. It’s reported that fatty-acid is an effective energy source that lead to the generation of ATP, and involved in the lipid metabolism of cancer cells, such as breast cancer cells [28, 29]. Plate et al. [30] indicated that the activation of IRE1-mediated unfolded protein response could influence endoplasmic reticulum quality control, and secretory proteostasis for destabilized, disease-relevant proteins. Furthermore, the damage of transmembrane transport would result in the imbalance of salt and water, and abnormality of blood pressure, which finally caused renal disease [31]. Therefore, fatty acid beta-oxidation, IRE1-mediated unfolded protein response, and transmembrane transport might play important roles in HTN. In addition, in our study, these DERs were also participated in many KEGG pathways like glycine, serine and threonine metabolism, carbon metabolism, biosynthesis of antibiotics, and citrate cycle (TCA cycle). Recent studies have shown that glycine, serine and threonine metabolism, and carbon metabolism affected the antioxidant ability of cells and their hyperactivation would drive the occurrence of cancer [32, 33]. However, the pathways of biosynthesis of antibiotics, and citrate cycle (TCA cycle) were not clear in HTN.

To understand the relationship between target genes and functions, a ceRNA regulatory network was established composing with 12 lncRNAs, 6 miRNAs and 226 mRNAs. In this network, lncRNAs KCTD21-AS1, LINC00470 and SNHG14 were found to be hub nodes. KCTD21-AS1, which promoted protein degradation and reduced cellular signaling, was associated with many diseases, including breast cancer [34] and obesity [35]. A study by Liu et al. indicated that LINC00470 was up-regulated in glioblastoma cell and its high level expression was an unfavorable prognosis marker for astrocytoma patients [36]. Moreover, SNHG14, as a key lncRNA, facilitated the migration and invasion of clear cell renal cell carcinoma via sponging miR–203 and accelerate Neural Wiskott-Aldrich syndrome protein [37]. Therefore, we supposed that lncRNAs KCTD21-AS1, LINC00470 and SNHG14 may be related to HTN.

In the end, a ceRNA regulatory network directly associated with HTN was proposed to find the biomarkers directly related to HTN. Our study showed that insulin signaling pathway was screened to directly associate with HTN. Insulin signaling pathway was reported to be damaged in the glomerulus by high glucose and promote apoptotic environment [38]. Therefore, insulin signaling pathway is possible to be a target for the prevention and treatment of HTN. In addition, based on the ceRNA regulatory network directly associated with HTN, insulin signaling pathway was modulated by mRNAs PPP1R3C, PPKAR2B and AKT3, miRNA hsa-miR–107, and lncRNAs SNHG14, TUG1, ZNF252P-AS1 and MIR503HG. MiRNA hsa-miR–107, as the hub in this network, was reported to be overexpressed in gastric carcinoma and promote tumor growth and survival [39]. Many researchers have shown that PPP1R3C, PPKAR2B and AKT3 may play important roles in cervical cancer [40], cardiovascular events [41] and prostate cancer [42]. In addition, lncRNAs TUG1, an important regulator of cancers, could facilitate proliferation and suppressed apoptosis by regulating miR–132–3p in osteosarcoma cells [43]. MIR503HG could inhibit migration and invasion of cells via miR–103/OLFM4 axis in triple negative breast cancer. However, the studies with regard to ZNF252P-AS1were few. Therefore, combined with our results, lncRNAs SNHG14, TUG1, ZNF252P-AS1 and MIR503HG were the biomarkers of HTN and HTN may be affected by insulin signaling pathway regulated by has-miR–107.

However, there were some limitations of our study. Firstly, the sample size is not big enough and the research was lack of validation of experiment. Furthermore, studies of these biomarkers needed to be further investigated on the clinical practice.

We identified 947 DERs, including 900 DE-mRNAs, 20 DE-lncRNAs and 27 DE-miRNAs. These DE-mRNAs were involved in the fatty-acid beta-oxidation, IRE1-mediated unfolded protein response, and transmembrane transport, and KEGG pathways like glycine, serine and threonine metabolism. Based on the ceRNA regulatory network, KCTD21-AS1, LINC00470 and SNHG14 were hub nodes, which targeted more miRNAs and mRNAs. Finally, the ceRNA regulatory network directly associated with HTN revealed that insulin signaling pathway was directly related to HTN and found the pathway was regulated by 3 mRNAs (PPP1R3C, PPKAR2B and AKT3),, 1 miRNA (has-miR–107) and 3 lncRNAs (SNHG14, TUG1, ZNF252P-AS1 and MIR503HG).. These results indicated that SNHG14, TUG1, ZNF252P-AS1 and MIR503HG regulated by has-miR–107 may be the biomarkers closely associated with HTN and would improve our understanding of the occurrence and development of HTN.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Competing interests

The authors declare that they have no competing interests.

Funding

This study was provided by Chinese National Natural Science Foundation (Project No. 81804032).

Authors’ contributions

Lei Kang and Zhen Wang designed the research and revised the manuscript. Yingxia Yang collected the data. Zhongjie Liu analyzed and interpreted the data. Zhen Wang obtained the funding and drafted the manuscript. All authors read and approved the final manuscript.

Acknowledgements

Not applicable.

1.Chen J, Xu J, Li Y, Zhang J, Chen H, Lu J, Wang Z, Zhao X, Xu K, Li X et al: Competing endogenous RNA network analysis identifies critical genes among the different breast cancer subtypes. Oncotarget 2017, 8(6):10171–10184.

2.Chen Z, Wu H, Wang G, Feng Y: Identification of potential candidate genes for hypertensive nephropathy based on gene expression profile. BMC nephrology 2016, 17(1):149.

3.Tan YG, Chan DWK, Yap FKP, Yap TL: Hypertensive urgency in nephrogenic diabetes insipidus with concomitant Hinman syndrome. BMJ case reports 2019, 12(7).

4.Wang XC, Liu CH, Chen YJ, Wu Y, Yang LS, Liu HM, Liao HL: Clinical and pathological analysis of the kidney in patients with hypertensive nephropathy. Experimental and therapeutic medicine 2013, 6(5):1243–1246.

5.Cao J, Hou R, Lu J, Zhang K, Zhao C, Jiang H, Feng Y, Wang Y: The predictive value of beta2-MG and TGF-beta for elderly hypertensive nephropathy. Experimental and therapeutic medicine 2019, 17(4):3065–3070.

6.Guerrot D, Dussaule JC, Mael-Ainin M, Xu-Dubois YC, Rondeau E, Chatziantoniou C, Placier S: Identification of periostin as a critical marker of progression/reversal of hypertensive nephropathy. PloS one 2012, 7(3):e31974.

7.Kurashige T, Takahashi T, Yamazaki Y, Nagano Y, Kondo K, Nakamura T, Yamawaki T, Tsuburaya R, Hayashi YK, Nonaka I et al: Elevated urinary beta2 microglobulin in the first identified Japanese family afflicted by X-linked myopathy with excessive autophagy. Neuromuscular disorders: NMD 2013, 23(11):911–916.

8.Rouse RL, Stewart SR, Thompson KL, Zhang J: Kidney injury biomarkers in hypertensive, diabetic, and nephropathy rat models treated with contrast media. Toxicologic pathology 2013, 41(4):662–680.

9.Franco-Zorrilla JM, Valli A, Todesco M, Mateos I, Puga MI, Rubio-Somoza I, Leyva A, Weigel D, Garcia JA, Paz-Ares J: Target mimicry provides a new mechanism for regulation of microRNA activity. Nat Genet 2007, 39(8):1033–1037.

10.Liu XH, Sun M, Nie FQ, Ge YB, Zhang EB, Yin DD, Kong R, Xia R, Lu KH, Li JH et al: Lnc RNA HOTAIR functions as a competing endogenous RNA to regulate HER2 expression by sponging miR–331–3p in gastric cancer. Molecular cancer 2014, 13:92.

11.Wang J, Liu X, Wu H, Ni P, Gu Z, Qiao Y, Chen N, Sun F, Fan Q: CREB up-regulates long non-coding RNA, HULC expression through interaction with microRNA–372 in liver cancer. Nucleic Acids Res 2010, 38(16):5366–5383.

12.Zhou X, Gao Q, Wang J, Zhang X, Liu K, Duan Z: Linc-RNA-RoR acts as a “sponge” against mediation of the differentiation of endometrial cancer stem cells by microRNA–145. Gynecologic oncology 2014, 133(2):333–339.

13.Liu H, Zhang Z, Wu N, Guo H, Zhang H, Fan D, Nie Y, Liu Y: Integrative Analysis of Dysregulated lncRNA-Associated ceRNA Network Reveals Functional lncRNAs in Gastric Cancer. Genes 2018, 9(6).

14.Zhou Q, Huang XR, Yu J, Yu X, Lan HY: Long Noncoding RNA Arid2-IR Is a Novel Therapeutic Target for Renal Inflammation. Molecular therapy: the journal of the American Society of Gene Therapy 2015, 23(6):1034–1043.

15.Marques FZ, Campain AE, Tomaszewski M, Zukowska-Szczechowska E, Yang YH, Charchar FJ, Morris BJ: Gene expression profiling reveals renin mRNA overexpression in human hypertensive kidneys and a role for microRNAs. Hypertension 2011, 58(6):1093–1098.

16.Parrish RS, Spencer HJ, 3rd: Effect of normalization on significance testing for oligonucleotide microarrays. Journal of biopharmaceutical statistics 2004, 14(3):575–589.

17.Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK: limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res 2015, 43(7):e47.

18.Braschi B, Denny P, Gray K, Jones T, Seal R, Tweedie S, Yates B, Bruford E: Genenames.org: the HGNC and VGNC resources in 2019. Nucleic Acids Res 2019, 47(D1):D786-D792.

19.Wang L, Cao C, Ma Q, Zeng Q, Wang H, Cheng Z, Zhu G, Qi J, Ma H, Nian H et al: RNA-seq analyses of multiple meristems of soybean: novel and alternative transcripts, evolutionary and functional implications. BMC plant biology 2014, 14:169.

20.Bien J, Tibshirani R: Hierarchical Clustering With Prototypes via Minimax Linkage. Journal of the American Statistical Association 2011, 106(495):1075–1084.

21.Szekely GJ, Rizzo ML: Hierarchical clustering via joint between-within distances: Extending Ward’s minimum variance method. Journal of Classification 2005, 22(2):151–183.

22.Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature Protocols 2009, 4(1):44–57.

23.Huang DW, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Research 2009, 37(1):1–13.

24.Paraskevopoulou MD, Vlachos IS, Karagkouni D, Georgakilas G, Kanellos I, Vergoulis T, Zagganas K, Tsanakas P, Floros E, Dalamagas T et al: DIANA-LncBase v2: indexing microRNA targets on non-coding transcripts. Nucleic Acids Res 2016, 44(D1):D231–238.

25.Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 2003, 13(11):2498–2504.

26.Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES et al: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proceedings of the National Academy of Sciences of the United States of America 2005, 102(43):15545–15550.

27.Davis AP, Grondin CJ, Johnson RJ, Sciaky D, McMorran R, Wiegers J, Wiegers TC, Mattingly CJ: The Comparative Toxicogenomics Database: update 2019. Nucleic Acids Res 2019, 47(D1):D948-D954.

28.Monaco ME: Fatty acid metabolism in breast cancer subtypes. Oncotarget 2017, 8(17):29487–29500.

29.Schafer ZT, Grassian AR, Song L, Jiang Z, Gerhart-Hines Z, Irie HY, Gao S, Puigserver P, Brugge JS: Antioxidant and oncogene rescue of metabolic defects caused by loss of matrix attachment. Nature 2009, 461(7260):109–113.

30.Plate L, Wiseman RL: Regulating Secretory Proteostasis through the Unfolded Protein Response: From Function to Therapy. Trends in cell biology 2017, 27(10):722–737.

31.Carlstrom M, Wilcox CS, Arendshorst WJ: Renal autoregulation in health and disease. Physiological reviews 2015, 95(2):405–511.

32.Amelio I, Cutruzzola F, Antonov A, Agostini M, Melino G: Serine and glycine metabolism in cancer. Trends in biochemical sciences 2014, 39(4):191–198.

33.Mehrmohamadi M, Liu X, Shestov AA, Locasale JW: Characterization of the usage of the serine metabolic network in human cancer. Cell reports 2014, 9(4):1507–1519.

34.Faryna M, Konermann C, Aulmann S, Bermejo JL, Brugger M, Diederichs S, Rom J, Weichenhan D, Claus R, Rehli M et al: Genome-wide methylation screen in low-grade breast cancer identifies novel epigenetically altered genes as potential biomarkers for tumor diagnosis. FASEB journal: official publication of the Federation of American Societies for Experimental Biology 2012, 26(12):4937–4950.

35.Yoganathan P, Karunakaran S, Ho MM, Clee SM: Nutritional regulation of genome-wide association obesity genes in a tissue-dependent manner. Nutrition & metabolism 2012, 9(1):65.

36.Liu C, Zhang Y, She X, Fan L, Li P, Feng J, Fu H, Liu Q, Zhao C, Sun Y et al: A cytoplasmic long noncoding RNA LINC00470 as a new AKT activator to mediate glioblastoma cell autophagy. Journal of hematology & oncology 2018, 11(1):77.

37.Liu G, Ye Z, Zhao X, Ji Z: SP1-induced up-regulation of lncRNA SNHG14 as a ceRNA promotes migration and invasion of clear cell renal cell carcinoma by regulating N-WASP. American journal of cancer research 2017, 7(12):2515–2525.

38.Katsoulieris EN, Drossopoulou GI, Kotsopoulou ES, Vlahakos DV, Lianos EA, Tsilibary EC: High Glucose Impairs Insulin Signaling in the Glomerulus: An In Vitro and Ex Vivo Approach. PloS one 2016, 11(7):e0158873.

39.Ren W, Zhang X, Li W, Feng Q, Feng H, Tong Y, Rong H, Wang W, Zhang D, Zhang Z et al: Exosomal miRNA–107 induces myeloid-derived suppressor cell expansion in gastric cancer. Cancer management and research 2019, 11:4023–4040.

40.Dai F, Chen G, Wang Y, Zhang L, Long Y, Yuan M, Yang D, Liu S, Cheng Y: Identification of candidate biomarkers correlated with the diagnosis and prognosis of cervical cancer via integrated bioinformatics analysis. Onco Targets Ther 2019, 12:4517–4532.

41.Voora D, Cyr D, Lucas J, Chi JT, Dungan J, McCaffrey TA, Katz R, Newby LK, Kraus WE, Becker RC et al: Aspirin exposure reveals novel genes associated with platelet function and cardiovascular events. Journal of the American College of Cardiology 2013, 62(14):1267–1276.

42.Lin HP, Lin CY, Huo C, Jan YJ, Tseng JC, Jiang SS, Kuo YY, Chen SC, Wang CT, Chan TM et al: AKT3 promotes prostate cancer proliferation cells through regulation of Akt, B-Raf, and TSC1/TSC2. Oncotarget 2015, 6(29):27097–27112.

43.Li G, Liu K, Du X: Long Non-Coding RNA TUG1 Promotes Proliferation and Inhibits Apoptosis of Osteosarcoma Cells by Sponging miR–132–3p and Upregulating SOX4 Expression. Yonsei medical journal 2018, 59(2):226–235.

Table 1 KEGG enrichment-pathway analysis of mRNA in the ceRNA regulatory network based on GSEA.

Name	Size	ES	NES	NOM P-value	Gene
KEGG _ insulin _ signaling _ pathway	3	-0.8969	-1.9100	1.193E-03	PPP1R3C, PRKAR2B, AKT3
KEGG _ glycine _ serine _ and _ threonine _ metabolism	5	0.7330	1.6968	6.270E-03	MAOB, SHMT1, GATM, PSAT1, PHGDH
KEGG _ pathways _ in _ cancer	3	-0.6399	-1.5624	2.632E-02	AKT3, KIT, RARB
KEGG _ lysosome	5	0.6561	1.5132	3.030E-02	NAPSA, SLC17A5, TPP1, ARSB, CTSB
KEGG _ apoptosis	2	-0.6784	-1.4150	4.599E-02	PRKAR2B, AKT3

Download PDF

Journal Publication

published 11 Nov, 2020

Read the published version in BMC Nephrology →

Editorial decision: Major revision
29 May, 2020
Review #2 received at journal
03 Apr, 2020
Reviewer #2 agreed at journal
19 Mar, 2020
Review #1 received at journal
12 Dec, 2019
Reviewer #1 agreed at journal
23 Nov, 2019
Reviewers invited by journal
22 Nov, 2019
Submission checks completed at journal
04 Nov, 2019
Editor assigned by journal
10 Oct, 2019
Editor invited by journal
09 Oct, 2019
First submitted to journal
08 Oct, 2019

You are reading this older preprint version

Read the latest preprint version →

Identification of biomarkers and pathways in hypertensive nephropathy based on the ceRNA regulatory network

Status:

Journal Publication

Version 1

Abstract

Figures

Background

method

Data pre-processing and annotation of RNA

Selection of DERs and functional analysis

Construction of ceRNA regulatory network and KEGG analysis

Establishment of ceRNA regulatory network directly associated with HTN

Results

Screening of DERs

GO terms and KEGG pathways of screened DE-mRNAs

Construction of ceRNA regulatory network

Establishment of ceRNA regulatory network directly associated with HTN

discussion

conclusions

declarations

Ethics approval and consent to participate

Consent for publication

Availability of data and materials

Competing interests

Funding

Authors’ contributions

Acknowledgements

references

tables

Status:

Journal Publication

Version 1