scholarly journals Identification and Verification of Five Potential Biomarkers Related to Skin and Thermal Injury Using Weighted Gene Co-Expression Network Analysis

2022 ◽  
Vol 12 ◽  
Author(s):  
Ronghua Yang ◽  
Zhengguang Wang ◽  
Jiehua Li ◽  
Xiaobing Pi ◽  
Xiaoxiang Wang ◽  
...  

Background: Burn injury is a life-threatening disease that does not have ideal biomarkers. Therefore, this study first applied weighted gene co-expression network analysis (WGCNA) and differentially expressed gene (DEG) screening methods to identify pivotal genes and diagnostic biomarkers associated with the skin burn process.Methods: After obtaining transcriptomic datasets of burn patient skin and normal skin from Gene Expression Omnibus (GEO) and performing differential analysis and functional enrichment, WGCNA was used to identify hub gene modules associated with burn skin processes in the burn patient peripheral blood sample dataset and determine the correlation between modules and clinical features. Enrichment analysis was performed to identify the functions and pathways of key module genes. Differential analysis, WGCNA, protein-protein interaction analysis, and enrichment analysis were utilized to screen for hub genes. Hub genes were validated in two other GEO datasets, tested by immunohistochemistry for hub gene expression in burn patients, and receiver operating characteristic curve analysis was performed. Finally, we constructed the specific drug activity, transcription factors, and microRNA regulatory network of the five hub genes.Results: A total of 1,373 DEGs in GSE8056 were obtained, and the top 5 upregulated genes were S100A12, CXCL8, CXCL5, MMP3, and MMP1, whereas the top 5 downregulated genes were SCGB1D2, SCGB2A2, DCD, TSPAN8, and KRT25. DEGs were significantly enriched in the immunity, epidermal development, and skin development processes. In WGCNA, the yellow module was identified as the most closely associated module with tissue damage during the burn process, and the five hub genes (ANXA3, MCEMP1, MMP9, S100A12, and TCN1) were identified as the key genes for burn injury status, which consistently showed high expression in burn patient blood samples in the GSE37069 and GSE13902 datasets. Furthermore, we verified using immunohistochemistry that these five novel hub genes were also significantly elevated in burn patient skin. In addition, MCEMP1, MMP9, and S100A12 showed perfect diagnostic performance in the receiver operating characteristic analysis.Conclusion: In conclusion, we analyzed the changes in genetic processes in the skin during burns and used them to identify five potential novel diagnostic markers in blood samples from burn patients, which are important for burn patient diagnosis. In particular, MCEMP1, MMP9, and S100A12 are three key blood biomarkers that can be used to identify skin damage in burn patients.

PeerJ ◽  
2020 ◽  
Vol 8 ◽  
pp. e10269
Author(s):  
Lingna Zhou ◽  
Liya Ding ◽  
Yuqi Gong ◽  
Jing Zhao ◽  
Gong Xin ◽  
...  

Background Host response diffuse large B-cell lymphoma (HR DLBCL) shares features of histologically defined T-cell/histiocyte-rich B-cell lymphoma, including fewer genetic abnormalities, frequent splenic and bone marrow involvement, and younger age at presentation. HR DLBCL is inherently less responsive to the standard treatment for DLBCL. Moreover, the mechanism of infiltration of HR DLBCL with preexisting abundant T-cells and dendritic cells is unknown, and their associated underlying immune responses incompletely defined. Here, hub genes and pathogenesis associated with HR DLBCL were explored to reveal molecular mechanisms and treatment targets. Methods Differentially expressed genes were identified in three datasets (GSE25638, GSE44337, GSE56315). The expression profile of the genes in the GSE53786 dataset was used to constructed a co-expression network. Protein-protein interactions analysis in the modules of interest identified candidate hub genes. Then screening of real hub genes was carried out by survival analysis within the GSE53786 and GSE10846 datasets. Expression of hub genes was validated in the Gene expression profiling interactive analysis, Oncomine databases and human tissue specimens. Functional enrichment analysis and Gene set enrichment analysis were utilized to investigate the potential mechanisms. Tumor Immune Estimation Resource and The Cancer Genome Atlas were used to mine the association of the hub gene with tumor immunity, potential upstream regulators were predicted using bioinformatics tools. Results A total of 274 common differentially expressed genes were identified. Within the key module, we identified CXCL10 as a real hub gene. The validation of upregulated expression level of CXCL10 was consistent with our study. CXCL10 might have a regulatory effect on tumor immunity. The predicted miRNA (hsa-mir-6849-3p) and transcription factor (IRF9) might regulate gene expression in the hub module.


2021 ◽  
Vol 2021 ◽  
pp. 1-10
Author(s):  
Tiancheng Zhang ◽  
Guihua Rao ◽  
Xiwen Gao

Background. Tuberculosis (TB) is a serious chronic bacterial infection caused by Mycobacterium tuberculosis (MTB). It is one of the deadliest diseases in the world and a heavy burden for people all over the world. However, the hub genes involved in the host response remain largely unclear. Methods. The data set GSE11199 was studied to clarify the potential gene network and signal transduction pathway in TB. The subjects were divided into latent tuberculosis and pulmonary tuberculosis, and the distribution of differentially expressed genes (DEGs) was analyzed between them using GEO2R. We verified the enriched process and pathway of DEGs by making use of the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO). The construction of protein-protein interaction (PPI) network of DEGs was achieved through making use of the Search Tool for the Retrieval of Interacting Genes (STRING), aiming at identifying hub genes. Then, the hub gene expression level in latent and pulmonary tuberculosis was verified by a boxplot. Finally, through making use of Gene Set Enrichment Analysis (GSEA), we further analyzed the pathways related to DEGs in the data set GSE11199 to show the changing pattern between latent and pulmonary tuberculosis. Results. We identified 98 DEGs in total in the data set GSE11199, 91 genes upregulated and 7 genes downregulated included. The enrichment of GO and KEGG pathways demonstrated that upregulated DEGs were mainly abundant in cytokine-mediated signaling pathway, response to interferon-gamma, endoplasmic reticulum lumen, beta-galactosidase activity, measles, JAK-STAT signaling pathway, cytokine-cytokine receptor interaction, etc. Based on the PPI network, we obtained 4 hub genes with a higher degree, namely, CTLA4, GZMB, GZMA, and PRF1. The box plot showed that these 4 hub gene expression levels in the pulmonary tuberculosis group were higher than those in the latent group. Finally, through Gene Set Enrichment Analysis (GSEA), it was concluded that DEGs were largely associated with proteasome and primary immunodeficiency. Conclusions. This study reveals the coordination of pathogenic genes during TB infection and offers the diagnosis of TB a promising genome. These hub genes also provide new directions for the development of latent molecular targets for TB treatment.


Author(s):  
Xitong Yang ◽  
Pengyu Wang ◽  
Shanquan Yan ◽  
Guangming Wang

AbstractStroke is a sudden cerebrovascular circulatory disorder with high morbidity, disability, mortality, and recurrence rate, but its pathogenesis and key genes are still unclear. In this study, bioinformatics was used to deeply analyze the pathogenesis of stroke and related key genes, so as to study the potential pathogenesis of stroke and provide guidance for clinical treatment. Gene Expression profiles of GSE58294 and GSE16561 were obtained from Gene Expression Omnibus (GEO), the differentially expressed genes (DEGs) were identified between IS and normal control group. The different expression genes (DEGs) between IS and normal control group were screened with the GEO2R online tool. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the DEGs were performed. Using the Database for Annotation, Visualization and Integrated Discovery (DAVID) and gene set enrichment analysis (GSEA), the function and pathway enrichment analysis of DEGS were performed. Then, a protein–protein interaction (PPI) network was constructed via the Search Tool for the Retrieval of Interacting Genes (STRING) database. Cytoscape with CytoHubba were used to identify the hub genes. Finally, NetworkAnalyst was used to construct the targeted microRNAs (miRNAs) of the hub genes. A total of 85 DEGs were screened out in this study, including 65 upward genes and 20 downward genes. In addition, 3 KEGG pathways, cytokine − cytokine receptor interaction, hematopoietic cell lineage, B cell receptor signaling pathway, were significantly enriched using a database for labeling, visualization, and synthetic discovery. In combination with the results of the PPI network and CytoHubba, 10 hub genes including CEACAM8, CD19, MMP9, ARG1, CKAP4, CCR7, MGAM, CD79A, CD79B, and CLEC4D were selected. Combined with DEG-miRNAs visualization, 5 miRNAs, including hsa-mir-146a-5p, hsa-mir-7-5p, hsa-mir-335-5p, and hsa-mir-27a- 3p, were predicted as possibly the key miRNAs. Our findings will contribute to identification of potential biomarkers and novel strategies for the treatment of ischemic stroke, and provide a new strategy for clinical therapy.


BMC Genomics ◽  
2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Ruoting Lin ◽  
Conor E. Fogarty ◽  
Bowei Ma ◽  
Hejie Li ◽  
Guoying Ni ◽  
...  

Abstract Background Papillary thyroid carcinoma (PTC) is the most common thyroid cancer. While many patients survive, a portion of PTC cases display high aggressiveness and even develop into refractory differentiated thyroid carcinoma. This may be alleviated by developing a novel model to predict the risk of recurrence. Ferroptosis is an iron-dependent form of regulated cell death (RCD) driven by lethal accumulation of lipid peroxides, is regulated by a set of genes and shows a variety of metabolic changes. To elucidate whether ferroptosis occurs in PTC, we analyse the gene expression profiles of the disease and established a new model for the correlation. Methods The thyroid carcinoma (THCA) datasets were downloaded from The Cancer Genome Atlas (TCGA), UCSC Xena and MisgDB, and included 502 tumour samples and 56 normal samples. A total of 60 ferroptosis related genes were summarised from MisgDB database. Gene set enrichment analysis (GSEA) and Gene set variation analysis (GSVA) were used to analyse pathways potentially involving PTC subtypes. Single sample GSEA (ssGSEA) algorithm was used to analyse the proportion of 28 types of immune cells in the tumour immune infiltration microenvironment in THCA and the hclust algorithm was used to conduct immune typing according to the proportion of immune cells. Spearman correlation analysis was performed on the ferroptosis gene expression and the correlation between immune infiltrating cells proportion. We established the WGCNA to identify genes modules that are highly correlated with the microenvironment of immune invasion. DEseq2 algorithm was further used for differential analysis of sequencing data to analyse the functions and pathways potentially involving hub genes. GO and KEGG enrichment analysis was performed using Clusterprofiler to explore the clinical efficacy of hub genes. Univariate Cox analysis was performed for hub genes combined with clinical prognostic data, and the results was included for lasso regression and constructed the risk regression model. ROC curve and survival curve were used for evaluating the model. Univariate Cox analysis and multivariate Cox analysis were performed in combination with the clinical data of THCA and the risk score value, the clinical efficacy of the model was further evaluated. Results We identify two subtypes in PTC based on the expression of ferroptosis related genes, with the proportion of cluster 1 significantly higher than cluster 2 in ferroptosis signature genes that are positively associated. The mutations of Braf and Nras are detected as the major mutations of cluster 1 and 2, respectively. Subsequent analyses of TME immune cells infiltration indicated cluster 1 is remarkably richer than cluster 2. The risk score of THCA is in good performance evaluated by ROC curve and survival curve, in conjunction with univariate Cox analysis and multivariate Cox analysis results based on the clinical data shows that the risk score of the proposed model could be used as an independent prognostic indicator to predict the prognosis of patients with papillary thyroid cancer. Conclusions Our study finds seven crucial genes, including Ac008063.2, Apoe, Bcl3, Acap3, Alox5ap, Atxn2l and B2m, and regulation of apoptosis by parathyroid hormone-related proteins significantly associated with ferroptosis and immune cells in PTC, and we construct the risk score model which can be used as an independent prognostic index to predict the prognosis of patients with PTC.


BMC Cancer ◽  
2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Baojie Wu ◽  
Shuyi Xi

Abstract Background This study aimed to explore and identify key genes and signaling pathways that contribute to the progression of cervical cancer to improve prognosis. Methods Three gene expression profiles (GSE63514, GSE64217 and GSE138080) were screened and downloaded from the Gene Expression Omnibus database (GEO). Differentially expressed genes (DEGs) were screened using the GEO2R and Venn diagram tools. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed. Gene set enrichment analysis (GSEA) was performed to analyze the three gene expression profiles. Moreover, a protein–protein interaction (PPI) network of the DEGs was constructed, and functional enrichment analysis was performed. On this basis, hub genes from critical PPI subnetworks were explored with Cytoscape software. The expression of these genes in tumors was verified, and survival analysis of potential prognostic genes from critical subnetworks was conducted. Functional annotation, multiple gene comparison and dimensionality reduction in candidate genes indicated the clinical significance of potential targets. Results A total of 476 DEGs were screened: 253 upregulated genes and 223 downregulated genes. DEGs were enriched in 22 biological processes, 16 cellular components and 9 molecular functions in precancerous lesions and cervical cancer. DEGs were mainly enriched in 10 KEGG pathways. Through intersection analysis and data mining, 3 key KEGG pathways and related core genes were revealed by GSEA. Moreover, a PPI network of 476 DEGs was constructed, hub genes from 12 critical subnetworks were explored, and a total of 14 potential molecular targets were obtained. Conclusions These findings promote the understanding of the molecular mechanism of and clinically related molecular targets for cervical cancer.


2021 ◽  
Author(s):  
Li Guoquan ◽  
Du Junwei ◽  
He Qi ◽  
Fu Xinghao ◽  
Ji Feihong ◽  
...  

Abstract BackgroundHashimoto's thyroiditis (HT), also known as chronic lymphocytic thyroiditis, is a common autoimmune disease, which mainly occurs in women. The early manifestation was hyperthyroidism, however, hypothyroidism may occur if HT was not controlled for a long time. Numerous studies have shown that multiple factors, including genetic, environmental, and autoimmune factors, were involved in the pathogenesis of the disease, but the exact mechanisms were not yet clear. The aim of this study was to identify differentially expressed genes (DEGs) by comprehensive analysis and to provide specific insights into HT. MethodsTwo gene expression profiles (GSE6339, GSE138198) about HT were downloaded from the Gene Expression Omnibus (GEO) database. The DEGs were assessed between the HT and normal groups using the GEO2R. The DEGs were then sent to the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. The hub genes were discovered using Cytoscape and CytoHubba. Finally, NetworkAnalyst was utilized to create the hub genes' targeted microRNAs (miRNAs). ResultsA total of 62 DEGs were discovered, including 60 up-regulated and 2 down-regulated DEGs. The signaling pathways were mainly engaged in cytokine interaction and cytotoxicity, and the DEGs were mostly enriched in immunological and inflammatory responses. IL2RA, CXCL9, IL10RA, CCL3, CCL4, CCL2, STAT1, CD4, CSF1R, and ITGAX were chosen as hub genes based on the results of the protein-protein interaction (PPI) network and CytoHubba. Five miRNAs, including mir-24-3p, mir-223-3p, mir-155-5p, mir-34a-5p, mir-26b-5p, and mir-6499-3p, were suggested as likely important miRNAs in HT. ConclusionsThese hub genes, pathways and miRNAs contribute to a better understanding of the pathophysiology of HT and offer potential treatment options for HT.


Hereditas ◽  
2021 ◽  
Vol 158 (1) ◽  
Author(s):  
Yun Tang ◽  
Xiaobo Yang ◽  
Huaqing Shu ◽  
Yuan Yu ◽  
Shangwen Pan ◽  
...  

Abstract Background Sepsis and septic shock are life-threatening diseases with high mortality rate in intensive care unit (ICU). Acute kidney injury (AKI) is a common complication of sepsis, and its occurrence is a poor prognostic sign to septic patients. We analyzed co-differentially expressed genes (co-DEGs) to explore relationships between septic shock and AKI and reveal potential biomarkers and therapeutic targets of septic-shock-associated AKI (SSAKI). Methods Two gene expression datasets (GSE30718 and GSE57065) were downloaded from the Gene Expression Omnibus (GEO). The GSE57065 dataset included 28 septic shock patients and 25 healthy volunteers and blood samples were collected within 0.5, 24 and 48 h after shock. Specimens of GSE30718 were collected from 26 patients with AKI and 11 control patents. AKI-DEGs and septic-shock-DEGs were identified using the two datasets. Subsequently, Gene Ontology (GO) functional analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, and protein-protein interaction (PPI) network analysis were performed to elucidate molecular mechanisms of DEGs. We also evaluated co-DEGs and corresponding predicted miRNAs involved in septic shock and AKI. Results We identified 62 DEGs in AKI specimens and 888, 870, and 717 DEGs in septic shock blood samples within 0.5, 24 and 48 h, respectively. The hub genes of EGF and OLFM4 may be involved in AKI and QPCT, CKAP4, PRKCQ, PLAC8, PRC1, BCL9L, ATP11B, KLHL2, LDLRAP1, NDUFAF1, IFIT2, CSF1R, HGF, NRN1, GZMB, and STAT4 may be associated with septic shock. Besides, co-DEGs of VMP1, SLPI, PTX3, TIMP1, OLFM4, LCN2, and S100A9 coupled with corresponding predicted miRNAs, especially miR-29b-3p, miR-152-3p, and miR-223-3p may be regarded as promising targets for the diagnosis and treatment of SSAKI in the future. Conclusions Septic shock and AKI are related and VMP1, SLPI, PTX3, TIMP1, OLFM4, LCN2, and S100A9 genes are significantly associated with novel biomarkers involved in the occurrence and development of SSAKI.


2020 ◽  
Author(s):  
Xi Pan ◽  
Jian-Hao Liu

Abstract Background Nasopharyngeal carcinoma (NPC) is a heterogeneous carcinoma that the underlying molecular mechanisms involved in the tumor initiation, progression, and migration are largely unclear. The purpose of the present study was to identify key biomarkers and small-molecule drugs for NPC screening, diagnosis, and therapy via gene expression profile analysis. Methods Raw microarray data of NPC were retrieved from the Gene Expression Omnibus (GEO) database and analyzed to screen out the potential differentially expressed genes (DEGs). The key modules associated with histology grade and tumor stage was identified by using weighted correlation network analysis (WGCNA). Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of genes in the key module were performed to identify potential mechanisms. Candidate hub genes were obtained, which based on the criteria of module membership (MM) and high connectivity. Then we used receiver operating characteristic (ROC) curve to evaluate the diagnostic value of hub genes. The Connectivity map database was further used to screen out small-molecule drugs of hub genes. Results A total of 430 DEGs were identified based on two GEO datasets. The green gene module was considered as key module for the tumor stage of NPC via WGCNA analysis. The results of functional enrichment analysis revealed that genes in the green module were enriched in regulation of cell cycle, p53 signaling pathway, cell part morphogenesis. Furthermore, four DEGs-related hub genes in the green module were considered as the final hub genes. Then ROC revealed that the final four hub genes presented with high areas under the curve, suggesting these hub genes may be diagnostic biomarkers for NPC. Meanwhile, we screened out several small-molecule drugs that have provided potentially therapeutic goals for NPC. Conclusions Our research identified four potential prognostic biomarkers and several candidate small-molecule drugs for NPC, which may contribute to the new insights for NPC therapy.


2020 ◽  
Vol 20 (1) ◽  
Author(s):  
Ming Chen ◽  
Junkai Zeng ◽  
Yeqing Yang ◽  
Buling Wu

Abstract Background Pulpitis is an inflammatory disease, the grade of which is classified according to the level of inflammation. Traditional methods of evaluating the status of dental pulp tissue in clinical practice have limitations. The rapid and accurate diagnosis of pulpitis is essential for determining the appropriate treatment. By integrating different datasets from the Gene Expression Omnibus (GEO) database, we analysed a merged expression matrix of pulpitis, aiming to identify biological pathways and diagnostic biomarkers of pulpitis. Methods By integrating two datasets (GSE77459 and GSE92681) in the GEO database using the sva and limma packages of R, differentially expressed genes (DEGs) of pulpitis were identified. Then, the DEGs were analysed to identify biological pathways of dental pulp inflammation with Gene Ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis and Gene Set Enrichment Analysis (GSEA). Protein–protein interaction (PPI) networks and modules were constructed to identify hub genes with the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and Cytoscape. Results A total of 470 DEGs comprising 394 upregulated and 76 downregulated genes were found in pulpitis tissue. GO analysis revealed that the DEGs were enriched in biological processes related to inflammation, and the enriched pathways in the KEGG pathway analysis were cytokine-cytokine receptor interaction, chemokine signalling pathway and NF-κB signalling pathway. The GSEA results provided further functional annotations, including complement system, IL6/JAK/STAT3 signalling pathway and inflammatory response pathways. According to the degrees of nodes in the PPI network, 10 hub genes were identified, and 8 diagnostic biomarker candidates were screened: PTPRC, CD86, CCL2, IL6, TLR8, MMP9, CXCL8 and ICAM1. Conclusions With bioinformatics analysis of merged datasets, biomarker candidates of pulpitis were screened and the findings may be as reference to develop a new method of pulpitis diagnosis.


TH Open ◽  
2020 ◽  
Vol 04 (04) ◽  
pp. e403-e412
Author(s):  
Aastha Mishra ◽  
Shankar Chanchal ◽  
Mohammad Z. Ashraf

AbstractSevere novel corona virus disease 2019 (COVID-19) infection is associated with a considerable activation of coagulation pathways, endothelial damage, and subsequent thrombotic microvascular injuries. These consistent observations may have serious implications for the treatment and management of this highly pathogenic disease. As a consequence, the anticoagulant therapeutic strategies, such as low molecular weight heparin, have shown some encouraging results. Cytokine burst leading to sepsis which is one of the primary reasons for acute respiratory distress syndrome (ARDS) drive that could be worsened with the accumulation of coagulation factors in the lungs of COVID-19 patients. However, the obscurity of this syndrome remains a hurdle in making decisive treatment choices. Therefore, an attempt to characterize shared biological mechanisms between ARDS and thrombosis using comprehensive transcriptomics meta-analysis is made. We conducted an integrated gene expression meta-analysis of two independently publicly available datasets of ARDS and venous thromboembolism (VTE). Datasets GSE76293 and GSE19151 derived from National Centre for Biotechnology Information–Gene Expression Omnibus (NCBI-GEO) database were used for ARDS and VTE, respectively. Integrative meta-analysis of expression data (INMEX) tool preprocessed the datasets and effect size combination with random effect modeling was used for obtaining differentially expressed genes (DEGs). Network construction was done for hub genes and pathway enrichment analysis. Our meta-analysis identified a total of 1,878 significant DEGs among the datasets, which when subjected to enrichment analysis suggested inflammation–coagulation–hypoxemia convolutions in COVID-19 pathogenesis. The top hub genes of our study such as tumor protein 53 (TP53), lysine acetyltransferase 2B (KAT2B), DExH-box helicase 9 (DHX9), REL-associated protein (RELA), RING-box protein 1 (RBX1), and proteasome 20S subunit beta 2 (PSMB2) gave insights into the genes known to be participating in the host–virus interactions that could pave the way to understand the various strategies deployed by the virus to improve its replication and spreading.


Sign in / Sign up

Export Citation Format

Share Document