This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Chronic kidney disease (CKD) is a worldwide public health problem with increasing incidence, poor prognosis, and high treatment cost. As a global health problem, CKD affects 10-16% of adults in Asia, Europe, and the United States and can progress to kidney failure [1]. The development of disease-related modules and genes is becoming increasingly popular. These methods are extremely useful in aiding the clinical search for diagnostic and therapeutic indicators. CKD is a complex disease related to genetic and environmental risk elements [2–4]. In response to the growing need to identify patients with CKD at an early stage and improve risk stratification for progression to end-stage renal disease, numerous studies have been conducted in large numbers of patients to investigate new and existing kidney disease biomarkers.
A number of different methodologies have been used, ranging from candidate single-gene studies to genome-wide multiomics analyses, to identify potential drug candidates. It is hoped that new methodologies will be developed to discover new CKD biomarkers, which will help us better understand the biology of kidney disease by using genetic, epigenetic, and transcriptome investigations [5]. CKD is an important issue given the increasing number of such patients worldwide. CKD is characterized by glomerular filtration rate (GFR) of less than 60 mL/min/1.73 m2 and signs of renal injury lasting at least 3 months. Reduced estimated GFR (eGFR) and severity of proteinuria independently predicted end-stage renal disease and mortality in patients with CKD [6, 7]. There is an urgent need to identify new biomarkers in patients with CKD to better detect people at high risk of rapid decline in kidney function, so that effective therapies can be used to curb disease progression [8, 9].
Weighted gene coexpression network analysis (WGCNA), as a method to screen disease-related modules, is the most common and useful method for discovering the link between genes and clinical characteristics [10]. Complex disorders including glioblastoma multiforme, cardiovascular disease in patients with diabetes, and Sjogren’s syndrome have been studied using WGCNA in prior studies [11–13]. WGCNA analysis provides an alternative approach for exploring genetic biomarkers that predict prognosis for CKD. The expressing data can be used to construct significantly correlated genes and their coexpression modules. In addition, these modules can further analyze modular characteristic genes (ME) and intramodular hub genes [14, 15]. Therefore, WGCNA may be a valuable tool for a comprehensive understanding of CKD-related genomic changes [10]. However, there was still little researches on CKD. In this study, WGCNA was performed to identify and screen related genes, so as to further explore the possible mechanisms of the critical genes and provide candidate markers for the diagnosis of CKD.
2. Materials and Methods
2.1. Microarray Data Source
The information was from the GEO datasets for patients with CKD extract (https://www.ncbi.nlm.nih.gov/gds/). The keywords “Chronic kidney disease” and “Homo sapiens” were applied as queries to search ckD-related datasets from GEO datasets. The GEO dataset met the following criteria: (1) the dataset contains CKD specimens and normal specimens; (2) each sample was assigned a group label; (iii) platform type is limited to “microarray”; (4) each probe has an available gene symbol or GeneBank ID; (5) the number of samples in the dataset was greater than 10. Finally, two datasets including GSE66494 and GSE62792 were selected as analysis datasets, including 65 CKD patients and 14 healthy controls.
2.2. Identification of Differentially Expressed Genes (DEGs)
In order to find differences in gene expression between CKD samples and healthy controls, GEO2R was used [16]. The log-fold changes in expressions and adjusted
2.3. Clustering and Enrichment Analysis
Clustering and enrichment analyses were carried out using GO (Gene Ontology) analysis. Pathway analysis is described in the Kyoto Encyclopedia of Genes and Genomes (KEGG). Enrichment analysis uses DAVID (https://david. http://ncifcrf.gov/tools. jsp), including biological process, molecular function, cells, and KEGG analysis. In addition, only FDR of GO or KEGG terms less than 0.05 was considered significant. We carried out visualization of the top 10 GO terms and the top 10 KEGG pathways.
2.4. Weighted Gene Coexpression Network Analysis
We compile and organize data files containing gene expression and phenotypes in standard formats. Firstly, to validate the accuracy of the study, we performed sample cluster analysis to verify the association of all data in the training queue. As a mean of ensuring that gene interaction followed a scale-free distribution, a study known as soft threshold selection analysis was employed. Besides, dynamic tree cutting algorithm was applied to identify modules via hierarchical clustering. Then, the protein expression abundance was clustered in R software (https://http://www.r-project.org/) to construct a weighted gene network. Subsequently, the correlation and correlation coefficient between the expressing spectrum and groups were calculated. We further identify important modules associated with traits. The dissimilarity degrees of MES in the module tree were calculated, and some modules (dissimilarity degree of
2.5. Construction of an Interactive Network
The molecular network was based on the interaction of text mining, experiment, database, coexpression, neighborhood, gene fusion, cooccurrence, and so on. Molecular networks and subnetworks were optimized for social networks constructed by STRING (https://http://string-db.org/cgi/input.pl) and Cytoscape software (version 3.7.2). Cytoscape software was used to depict the results from STRING as PPI networks. The PPI network was cleansed of nodes that were not connected to any other nodes. Degree centrality was used to determine the PPI network’s hub genes. The subnetwork was analyzed and scored using Cytoscape’s MCODE plug-in.
2.6. Statistical Analysis
Statistical analysis was performed by SPSS software (SPSS Inc., Chicago, IL, USA) and R.4.1.1 (R Core Team, Massachusetts, USA). Data visualization was performed using PRISM. The differential expression threshold was set to 1.5.
3. Results
3.1. Differential Expression Analysis in CKD Patients
First, we screened gene expression profiling datasets from CKD patients and healthy individuals in the GEO database. A total of renal biopsy specimens of CKD patients were selected for microarray analysis data, and datasets GSE66494 and GSE62792 were selected as analysis datasets, including data of 65 CKD patients and 14 healthy controls. Data were obtained from the platform Agilent-014850 Whole Human Genome Microarray 4x44K G4112F (Probe Name version) and stored in the GPL6480 platform. We calculated DEGs between CKD patients and healthy individuals (Figures 1(a) and 1(b)). We found that a total of 882 DEGs were screened out in the above two datasets (Figure 1(c)). We found that a total of 882 DEGs were screened in the above two datasets, and their gene expression profiles were listed in the heatmap (Figure 1(d)). Additionally, DEGs were observed to significantly distinguish CKD patients from healthy controls in principal component analysis (Figure 1(e)).
[figure(s) omitted; refer to PDF]
3.2. Clustering and Enrichment Analysis
In order to further learn the function of CKD-related genes, we conducted GO function and KEGG pathway enrichment analysis of related genes. GO enrichment results showed that, for BP, these related molecules mainly play physiological functions including cell-cell adhesion, in utero embryonic development, VEGF receptor signaling pathway, negative regulation of transcription, nuclear-transcribed mRNA poly (A), mitotic cytokinesis, intracellular signal transduction, positive regulation of transcription, mitotic metaphase plate congression, and multicellular organism growth (Figure 2(a)). The results showed that for CC, these related molecules were mainly located in the intracellular, organelle, intracellular part, intracellular organelle, membrane-bounded organelle, cytoplasm, cytoplasmic part, intracellular membrane, organelle part, intracellular organelle part (Figure 2(a)). The results showed that, for MF, the main physiological functions of these related molecules included protein binding, binding, poly (A) RNA binding, nucleotide binding, nucleoside phosphate binding, small molecule binding, heterocyclic compound binding, organic cyclic compound binding, RNA binding, and carbohydrate derivative binding (Figure 2(a)). In order to understand the enrichment of the pathway, KEGG analysis was used to analyze the pathway, and it was found that these molecules were mainly involved in 14 pathways including protein processing in endoplasmic reticulum, pathways in cancer, proteoglycans in cancer, TNF signaling pathway, platelet activation, aldosterone synthesis and secretion, adrenergic signaling in cardiomyocytes, focal adhesion, N-glycan biosynthesis, and HTLV-I infection (Figure 2(b)).
[figure(s) omitted; refer to PDF]
3.3. Construction of Coexpression Networks
WGCNA was used to identify disease-related modules in which genes exhibited coordinated expression patterns, which greatly improved the chances of identifying hub genes. In order to construct gene coexpression network, GSE66494 and GSE62792 data were used for cluster analysis. A total of 79 samples, including 10891 gene expression data, were used to construct hierarchical clustering trees. The analysis results showed that no obvious outlier samples were found. Thus, the analysis program retained all samples to construct the weighted coexpression network (Figure 3(a)). The dynamic mixed shearing methods were used to merge the modules with high similarity of feature genes, and 10 gene modules with different colors were finally obtained, among which the gray module was the gene without coexpression (Figure 3(b)). In addition, we also conducted hierarchical clustering of these gene modules, and these characteristic gene modules could be grouped into two categories. Disease was considered as the main clinical feature for correlation analysis of different gene modules. MEblack and MEpink were negatively correlated with the disease, while MEyellow, MEturquiose, MEblue, MEred, MEpurple, MEmagenta, MEbrown, and MEgreen were positively correlated with the disease. However, the results showed that MEturquiose, MEred, and MEblue in group were significantly correlated with disease (Figure 3(c)). Then, the MEturquiose, MEblue, and MEred were significantly positively correlated modules in the CKD and healthy groups, and it turned out that there was a complex network of connections between these molecules (Figures 3(d)–3(f)).
[figure(s) omitted; refer to PDF]
3.4. Interaction Network Analysis
To understand the interaction network status of CKD-related molecules, the interaction network was constructed through STRING. It turned out that there was a complex network of connections between these molecules in DEGs (Figure 4(a)). We analyzed and extracted key subnetworks through the Cytoscape plugin MCODE. The results showed that there were 3 subnetworks with high scores, including 14 key genes in total (Figures 4(b)–4(d)). Most importantly, a subnetwork containing 6 nodes was found with a score of 12.08, including CCAAT enhancer binding protein zeta (CEBPZ), interferon gamma inducible protein 16 (IFI16), Ly1 antibody reactive (LYAR), biogenesis of ribosomes BRX1 (BRIX1), BMS1 ribosome biogenesis factor (BMS1), and DEAD-box helicase 18 (DDX18). We used these molecules as candidate markers to distinguish CKD from healthy people.
[figure(s) omitted; refer to PDF]
3.5. Validation of Candidate Markers
First, we aggregated the details, including degree, MCODE_score in Cytoscape, fold changes, and
[figure(s) omitted; refer to PDF]
4. Discussion
CKD has a wide range of underlying causes, including both hereditary and environmental factors, and is considered to be a global public health problem, with the adjusted prevalence of CKD in the European adult population ranging from 3.3% to 17.3% [17]. CKD is more common among the elderly and is associated with a higher risk of cardiovascular disease (CVD). One of the most common causes of kidney CKD is diabetes mellitus (DM). In addition, few biomarkers have been found in clinical practice, although other diagnostic biomarkers for CKD have been studied [18]. Therefore, improved early detection and treatment of CKD necessitated new molecular biomarkers.
IFI16 has been shown to influence ribosome biogenesis and has been identified as a promoter associated with stem cell-like properties in colorectal cancer [19]. LYAR can enhance the stem cell-like properties of breast cancer and lead to poor prognosis of breast cancer, which is expected to be a potential biomarker for breast cancer treatment [20]. BRIX1 is a new potential target in psoriasis and diffuse superficial actinic sweat keratosis [21]. BMS1 is an RNA- and DNA-binding protein involved in nucleolar processing of 7S to 5.8SrRNA. When exposed to cytotoxic agents, the nucleolar localization of PUF-A redistributes into the nuclear cytoplasm. DDX18 is a risk site associated with DTC susceptibility [22]. CEBPZ was also involved in cell growth and differentiation, especially hematopoietic differentiation [23]. DDX18 was identified to be associated with stroke, and serum PDCD11-AB levels may serve as a potential biomarker for TRANSIENT ischemic attack [24, 25].
The proliferation of omics-related biomarker studies over the past decade reflects the need for new, effective, noninvasive tools that can identify people at risk for CKD and help target kidney disease management [26]. CKD was a growing public health problem with high morbidity and mortality. New biomarkers were developed to improve risk stratification and clinical decision-making and to guide the enrichment of patients with CKD in clinical trials. Despite tremendous efforts, only a few biomarkers have so far found large-scale clinical application. Although our study has identified multiple putative biomarkers, these were mainly from small, single-center studies. The utility of such biomarkers needed to be confirmed in different populations and in larger cohorts. In addition, no in vivo or in vitro studies were performed. Both the above deficiencies are also the main aspects of our further research.
5. Conclusions
We identified 6 hub genes in CKD, which ere demonstrated in the validation dataset. They could be used these molecules as candidate markers to distinguish CKD from healthy people. Our study indicated that 6 hub genes may play key roles in the occurrence and development of CKD.
Authors’ Contributions
Jia Wang and Yuan Yin contributed equally to this work.
Acknowledgments
This work was supported by the Traditional Chinese Medicine Scientific Research Project of Shanghai Hongkou District Health Commission in 2020 (HKQ-ZYY-2020-33).
[1] S. I. Hallan, J. Coresh, B. C. Astor, A. Asberg, N. R. Powe, S. Romundstad, H. A. Hallan, S. Lydersen, J. Holmen, "International comparison of the relationship of chronic kidney disease prevalence and ESRD risk," Journal of the American Society of Nephrology, vol. 17 no. 8, pp. 2275-2284, DOI: 10.1681/ASN.2005121273, 2006.
[2] N. H. Arar, V. S. Voruganti, S. D. Nath, F. Thameem, R. Bauer, S. A. Cole, J. Blangero, J. W. Mac Cluer, A. G. Comuzzie, H. E. Abboud, "A genome-wide search for linkage to chronic kidney disease in a community-based sample: the SAFHS," Nephrology, Dialysis, Transplantation, vol. 23 no. 10, pp. 3184-3191, DOI: 10.1093/ndt/gfn215, 2008.
[3] K. Nabrdalik, J. Gumprecht, P. Adamczyk, S. Gorczynska-Kosiorz, J. Zywiec, W. Grzeszczak, "Association of rs 1800471 polymorphism of TGFB1 gene with chronic kidney disease occurrence and progression and hypertension appearance," Archives of Medical Science, vol. 9 no. 2, pp. 230-237, DOI: 10.5114/aoms.2013.34418, 2013.
[4] C. M. O'Seaghdha, C. S. Fox, "Genome-wide association studies of chronic kidney disease: what have we learned?," Nature Reviews. Nephrology, vol. 8 no. 2, pp. 89-99, DOI: 10.1038/nrneph.2011.189, 2011.
[5] J. Yang, F. H. J. Claas, M. Eikmans, "Genome-wide association studies in kidney transplantation: advantages and constraints," Transplant Immunology, vol. 49,DOI: 10.1016/j.trim.2018.04.004, 2018.
[6] B. C. Astor, K. Matsushita, R. T. Gansevoort, M. van der Velde, M. Woodward, A. S. Levey, P. E. Jong, J. Coresh, C Chronic Kidney Disease Prognosis, M. El-Nahas, K. U. Eckardt, B. L. Kasiske, J. Wright, L. Appel, T. Greene, A. Levin, O. Djurdjev, D. C. Wheeler, M. J. Landray, J. N. Townend, J. Emberson, L. E. Clark, A. Macleod, A. Marks, T. Ali, N. Fluck, G. Prescott, D. H. Smith, J. R. Weinstein, E. S. Johnson, M. L. Thorp, J. F. Wetzels, P. J. Blankestijn, A. D. van Zuilen, V. Menon, M. Sarnak, G. Beck, F. Kronenberg, B. Kollerits, M. Froissart, B. Stengel, M. Metzger, G. Remuzzi, P. Ruggenenti, A. Perna, H. J. Heerspink, B. Brenner, D. de Zeeuw, P. Rossing, H. H. Parving, P. Auguste, K. Veldhuis, Y. Wang, L. Camarata, B. Thomas, T. Manley, "Lower estimated glomerular filtration rate and higher albuminuria are associated with mortality and end-stage renal disease. A collaborative meta- analysis of kidney disease population cohorts," Kidney International, vol. 79 no. 12, pp. 1331-1340, DOI: 10.1038/ki.2010.550, 2011.
[7] F National Kidney, "K/DOQI clinical practice guidelines for chronic kidney disease: evaluation, classification, and stratification," American Journal of Kidney Diseases, vol. 39, 2002.
[8] R. G. Fassett, S. K. Venuthurupalli, G. C. Gobe, J. S. Coombes, M. A. Cooper, W. E. Hoy, "Biomarkers in chronic kidney disease: a review," Kidney International, vol. 80 no. 8, pp. 806-821, DOI: 10.1038/ki.2011.198, 2011.
[9] A. S. Levey, J. Coresh, "Chronic kidney disease," Lancet, vol. 379 no. 9811, pp. 165-180, DOI: 10.1016/S0140-6736(11)60178-5, 2012.
[10] P. Langfelder, S. Horvath, "WGCNA: an R package for weighted correlation network analysis," BMC Bioinformatics, vol. 9 no. 1,DOI: 10.1186/1471-2105-9-559, 2008.
[11] W. Liang, F. Sun, Y. Zhao, L. Shan, H. Lou, "Identification of susceptibility modules and genes for cardiovascular disease in diabetic patients using WGCNA analysis," Journal Diabetes Research, vol. 2020, article 4178639,DOI: 10.1155/2020/4178639, 2020.
[12] Q. Yang, R. Wang, B. Wei, C. Peng, L. Wang, G. Hu, D. Kong, C. Du, "Candidate biomarkers and molecular mechanism investigation for glioblastoma multiforme utilizing WGCNA," BioMed Research International, vol. 2018,DOI: 10.1155/2018/4246703, 2018.
[13] Q. Yao, Z. Song, B. Wang, Q. Qin, J. A. Zhang, "Identifying key genes and functionally enriched pathways in Sjögren’s syndrome by weighted gene co-expression network analysis," Frontiers in Genetics, vol. 10,DOI: 10.3389/fgene.2019.01142, 2019.
[14] Z. Wang, Z. Lyu, L. Pan, G. Zeng, P. Randhawa, "Defining housekeeping genes suitable for RNA-seq analysis of the human allograft kidney biopsy tissue," BMC Medical Genomics, vol. 12 no. 1,DOI: 10.1186/s12920-019-0538-z, 2019.
[15] L. Xia, X. Su, J. Shen, Q. Meng, J. Yan, C. Zhang, Y. Chen, H. Wang, M. Xu, "ANLN functions as a key candidate gene in cervical cancer as determined by integrated bioinformatic analysis," Cancer Management and Research, vol. 10, pp. 663-670, DOI: 10.2147/CMAR.S162813, 2018.
[16] T. Barrett, S. E. Wilhite, P. Ledoux, C. Evangelista, I. F. Kim, M. Tomashevsky, K. A. Marshall, K. H. Phillippy, P. M. Sherman, M. Holko, A. Yefanov, H. Lee, N. Zhang, C. L. Robertson, N. Serova, S. Davis, A. Soboleva, "NCBI GEO: archive for functional genomics data sets--update," Nucleic Acids Research, vol. 41, pp. D991-D995, DOI: 10.1093/nar/gks1193, 2013.
[17] A. S. Levey, R. Atkins, J. Coresh, E. P. Cohen, A. J. Collins, K. U. Eckardt, M. E. Nahas, B. L. Jaber, M. Jadoul, A. Levin, N. R. Powe, J. Rossert, D. C. Wheeler, N. Lameire, G. Eknoyan, "Chronic kidney disease as a global public health problem: approaches and initiatives - a position statement from kidney disease improving global outcomes," Kidney International, vol. 72 no. 3, pp. 247-259, DOI: 10.1038/sj.ki.5002343, 2007.
[18] K. Bruck, V. S. Stel, G. Gambaro, S. Hallan, H. Volzke, J. Arnlov, M. Kastarinen, I. Guessous, J. Vinhas, B. Stengel, H. Brenner, J. Chudek, S. Romundstad, C. Tomson, A. O. Gonzalez, A. K. Bello, J. Ferrieres, L. Palmieri, G. Browne, V. Capuano, W. Van Biesen, C. Zoccali, R. Gansevoort, G. Navis, D. Rothenbacher, P. M. Ferraro, D. Nitsch, C. Wanner, K. J. Jager, "CKD prevalence varies across the European general population," Journal of the American Society of Nephrology, vol. 27 no. 7, pp. 2135-2147, DOI: 10.1681/ASN.2015050542, 2016.
[19] Y. Zou, J. Zhang, L. Zhang, X. Yan, "Interferon-induced protein 16 expression in colorectal cancer and its correlation with proliferation and immune signature markers," Oncology Letters, vol. 22 no. 3,DOI: 10.3892/ol.2021.12948, 2021.
[20] S. Li, J. Ma, A. Zheng, X. Song, S. Chen, F. Jin, "DEAD-box helicase 27 enhances stem cell-like properties with poor prognosis in breast cancer," Journal of Translational Medicine, vol. 19 no. 1,DOI: 10.1186/s12967-021-03011-0, 2021.
[21] F. Abdallah, C. Pichon, "Evidence on the direct correlation between miR-31 and IL-22 axis in IMQ- induced psoriasis," Experimental Dermatology, vol. 28 no. 11, pp. 1336-1340, DOI: 10.1111/exd.14001, 2019.
[22] G. Figlioli, R. Elisei, C. Romei, O. Melaiu, M. Cipollini, F. Bambi, B. Chen, A. Kohler, A. Cristaudo, K. Hemminki, F. Gemignani, A. Forsti, S. Landi, "A comprehensive meta-analysis of case-control association studies to evaluate polymorphisms associated with the risk of differentiated thyroid carcinoma," Cancer Epidemiology, Biomarkers & Prevention, vol. 25 no. 4, pp. 700-713, DOI: 10.1158/1055-9965.EPI-15-0652, 2016.
[23] D. Chen, Y. Qin, M. Dai, L. Li, H. Liu, Y. Zhou, C. Qiu, Y. Chen, Y. Jiang, "BGN and COL11A1 regulatory network analysis in colorectal cancer (CRC) reveals that BGN influences CRC cell biological functions and interacts with miR-6828-5p," Cancer Management and Research, vol. 12, pp. 13051-13069, DOI: 10.2147/CMAR.S277261, 2020.
[24] D. F. Merlo, S. Agramunt, L. Anna, H. Besselink, M. Botsivali, N. J. Brady, M. Ceppi, L. Chatzi, B. Chen, I. Decordier, P. B. Farmer, S. Fleming, V. Fontana, A. Försti, E. Fthenou, F. Gallo, P. Georgiadis, H. Gmuender, R. W. Godschalk, B. Granum, L. J. Hardie, K. Hemminki, K. Hochstenbach, L. E. Knudsen, M. Kogevinas, K. Kovács, S. A. Kyrtopoulos, M. Løvik, J. K. Nielsen, U. C. Nygaard, M. Pedersen, P. Rydberg, B. Schoket, D. Segerbäck, R. Singh, J. Sunyer, M. Törnqvist, H. van Loveren, F. J. van Schooten, K. V. Loock, H. von Stedingk, J. Wright, J. C. Kleinjans, M. Kirsch-Volders, J. H. M. van Delft, the NewGeneris Consortium, "Micronuclei in cord blood lymphocytes and associations with biomarkers of exposure to carcinogens and hormonally active factors, gene polymorphisms, and gene expression: the New Generis cohort," Environmental Health Perspectives, vol. 122 no. 2, pp. 193-200, DOI: 10.1289/ehp.1206324, 2014.
[25] Y. Yoshida, H. Wang, T. Hiwasa, T. Machida, E. Kobayashi, S. Mine, G. Tomiyoshi, R. Nakamura, N. Shinmen, H. Kuroda, H. Takizawa, K. Kashiwado, I. Kamitsukasa, H. Shin, T. Wada, A. Aotsuka, E. Nishi, M. Ohno, M. Takemoto, K. Yokote, S. Takahashi, J. Matsushima, X. M. Zhang, M. Takiguchi, Y. Iwadate, "Elevation of autoantibody level against PDCD11 in patients with transient ischemic attack," Oncotarget, vol. 9 no. 10, pp. 8836-8848, DOI: 10.18632/oncotarget.23653, 2018.
[26] A. Levin, P. E. Stevens, "Summary of KDIGO 2012 CKD guideline: behind the scenes, need for guidance, and a framework for moving forward," Kidney International, vol. 85 no. 2014, pp. 49-61, DOI: 10.1038/ki.2013.444, .
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright © 2022 Jia Wang et al. This is an open access article distributed under the Creative Commons Attribution License (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License. https://creativecommons.org/licenses/by/4.0/
Abstract
Chronic kidney disease (CKD) is an ongoing deterioration of renal function that often progresses to end-stage renal disease. In this study, we aimed to screen and identify potential key genes for CKD using the weighted gene coexpression network (WGCNA) analysis tool. Gene expression data related to CKD were screened from GEO database, and expression datasets of GSE66494 and GSE62792 were obtained. After discrete analysis of samples, WGCNA analysis was performed to construct gene coexpression module, and the correlation between the module and disease was calculated. The modules with a significant correlation with the disease were selected for Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Then, the interaction network of related molecules was constructed, and the high score subnetwork was selected, and the candidate key molecules were identified. A total of 882 DEGs were identified in the screening datasets. A subnetwork containing 6 nodes was found with a high score of 12.08, including CEBPZ, IFI16, LYAR, BRIX1, BMS1, and DDX18. DEGs could significantly differentiate CKD and healthy individuals in principal component analysis. In addition, the MEturquiose, MEred, and MEblue in group were significantly correlated with disease in WGCNA. These 6 hub genes were found to significantly discriminate between CKD and healthy controls in the validation dataset, suggesting that they could use these molecules as candidate markers to distinguish CKD from healthy people. Overall, our study indicated that 6 hub genes may play key roles in the occurrence and development of CKD.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer