A Signature for Smoking Status of Coronary Heart Disease Patients through Weighted Gene Coexpression Network Analysis.
Ontology highlight
ABSTRACT: Background
Smoking is one of the risk factors of coronary heart disease (CHD), while its underlying mechanism is less well defined.Purpose
To identify and testify 6 key genes of CHD related to smoking through weighted gene coexpression network analysis (WGCNA), protein-protein interaction (PPI) network analysis, and pathway analysis.Methods
CHD patients' samples were first downloaded from Gene Expression Omnibus (GEO). Then, genes of interest were obtained after analysis of variance (ANOVA). Thereafter, 23 coexpressed modules that were determined after genes with similar expression were incorporated via WGCNA. The biological functions of genes in the modules were researched by enrichment analysis. Pearson correlation analysis and PPI network analysis were used to screen core genes related to smoking in CHD.Results
The violet module was the most significantly associated with smoking (r = -0.28, p = 0.006). Genes in this module mainly participated in biological functions related to the heart. Altogether, 6 smoking-related core genes were identified through bioinformatics analyses. Their expressions in animal models were detected through the animal experiment.Conclusion
This study identified 6 core genes to serve as underlying biomarkers for monitoring and predicting smoker's CHD risk.
SUBMITTER: Shang Z
PROVIDER: S-EPMC8791244 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundCoronary heart disease (CHD) is the most prevalent disease with an unelucidated pathogenetic mechanism and is mediated by complex molecular interactions of exosomes. Here, we aimed to identify differentially expressed exosome genes for the disease development and prognosis of CHD.MethodSix CHD samples and 32 normal samples were downloaded from the exoRbase database to identify the candidate genes in the CHD. The differentially expressed genes (DEGs) were identified. And then, weighted gene correlation network analysis (WGCNA) was used to investigate the modules in coexpressed genes between CHD samples and normal samples. DEGs and the module of the WGCNA were intersected to obtain the most relevant exosome genes. After that, the function enrichment analyses and protein-protein interaction network (PPI) were performed for the particular module using STRING and Cytoscape software. Finally, the CIBERSORT algorithm was used to analyze the immune infiltration of exosome genes between CHD samples and normal samples.ResultWe obtain a total of 715 overlapping exosome genes located at the intersection of the DEGs and key modules. The Gene Ontology enrichment of DEGs in the blue module included inflammatory response, neutrophil degranulation, and activation of CHD. In addition, protein-protein networks were constructed, and hub genes were identified, such as LYZ, CAMP, HP, ORM1, and LTF. The immune infiltration profiles varied significantly between normal controls and CHD. Finally, we found that mast cells activated and eosinophils had a positive correlation. B cell memory had a significant negative correlation with B cell naive. Besides, neutrophils and mast cells were significantly increased in CHD patients.ConclusionThe underlying mechanism may be related to neutrophil degranulation and the immune response. The hub genes and the difference in immune infiltration identified in the present study may provide new insights into the diagnostic and provide candidate targets for CHD.
| S-EPMC8536412 | biostudies-literature
Project description:The present study sought to identify potential hub genes and pathways of acute coronary syndrome (ACS). We downloaded the dataset (GSE56045) from the Gene Expression Omnibus (GEO) database and analyzed weighted gene coexpression networks (WGCNA). Gene Ontology annotation, Disease Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed using R software. The protein-protein interaction (PPI) network was constructed using Cytoscape, and the Molecular Complex Detection app was employed to identify significant modules and hub genes. The hub genes were then validated in other microarrays and patients by RT-PCR. Two modules were identified and associated with coronary artery disease (CAD) and included 219 genes. After function and PPI analyses, 24 genes were identified to be potentially associated with CAD. Linear correlation was performed to calculate the relationship between the gene expression levels and coronary artery calcification score and found that CCR7 (R = -0.081, P = 0.0065), CD2 (R = -0.075, P = 0.0012), CXCR5 (R = -0.065, P = 0.029) and IL7R (R = -0.06, P = 0.043) should be validated in other dataset. By comparing the gene expression levels in different groups in GSE23561, GSE34822, GSE59867, GSE60993 and GSE129935, only two genes (CCR7 and CXCR5) showed significance. The nomogram showed that CXCR5 showed the risk of ACS. Further analysis in chest patients found CXCR5 played a key role resulting in ACS. Our WGCNA analysis identified CXCR5 as a risk factor for ACS, and the potential pathogenesis may be associated with immune inflammation.
| S-EPMC7732301 | biostudies-literature
Project description:BackgroundEwing sarcoma (ES) is the second most common pediatric bone tumor with a high rate of metastasis, high recurrence, and low survival rate. Therefore, the identification of new biomarkers which can improve the prognosis of ES patients is urgently needed.MethodsHere, GSE17679 dataset was downloaded from GEO databases. WGCNA method was used to identify one module associating with OVS (overall vital survival) and event. cytoHubba was used to screen out 50 hub genes from the module genes. Then, GSE17679 dataset was randomly divided into train cohort and test cohort. Next, univariate Cox analysis, LASSO regression analysis, and multivariate Cox analysis were conducted on 50 hub genes combined with train cohort data to select pivotal genes. Finally, an optimal 7-gene-based risk assessment model was established, which was verified by test cohort, entire GSE17679, and two independent datasets (GSE63157 and TCGA-SARC).ResultsThe results of the functional enrichment analysis revealed that the OVS and event-associated module were mainly enriched in the protein transcription, cell proliferation, and cell-cycle control. And the train cohort was divided into high-risk and low-risk subgroups based on the median risk score; the results showed that the survival of the low-risk subgroup was significantly longer than high-risk. ROC analysis revealed that AUC values of 1, 3, and 5-year survival were 0.85, 0.94, and 0.88, and Kaplan-Meier analysis also revealed that P value < 0.0001, indicating that this model was accurate, which was also verified in the test, entire cohort, and two independent datasets (GSE63157 and TCGA-SARC). Then, we performed a comprehensive analysis (differential expression analysis, correlation analysis and survival analysis) of seven pivotal genes, and found that four genes (NCAPG, KIF4A, NUF2 and CDC20) plays a more crucial role in the prognosis of ES.ConclusionTaken together, this study established an optimal 7-gene-based risk assessment model and identified 4 potential therapeutic targets, to improve the prognosis of ES patients.
| S-EPMC8814720 | biostudies-literature
Project description:Intracranial aneurysm (IA) can cause fatal subarachnoid hemorrhage (SAH) after rupture, and identifying patients with unruptured IAs is essential for reducing SAH fatalities. The epithelial-mesenchymal transition (EMT) may be vital to IA progression. Here, identified key EMT-related genes in aneurysms and their pathogenic mechanisms via bioinformatic analysis. The GSE13353, GSE75436, and GSE54083 datasets from Gene Expression Omnibus were analyzed with limma to identify differentially expressed genes (DEGs) among unruptured aneurysms, ruptured aneurysms, and healthy samples. The results revealed that three EMT-related DEGs (ADIPOQ, WNT11, and CCL21) were shared among all groups. Coexpression modules and hub genes were identified via weighted gene co-expression network analysis, revealing two significant modules (red and green) and 14 EMT-related genes. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses suggested that cytokine interactions were closely related. Gene set enrichment analysis revealed that unruptured aneurysms were enriched for the terms "inflammatory response" and "vascular endothelial growth". Protein-protein interaction analysis identified seven key genes, which were evaluated with the GSE54083 dataset to determine their sensitivity and specificity. In the external validation set, we verified the differential expression of seven genes in unruptured aneurysms and normal samples. Together, these findings indicate that FN1, and SPARC may help distinguish normal patients from patients with asymptomatic IAs.
| S-EPMC8741966 | biostudies-literature
Project description:Systems-oriented genetic approaches that incorporate gene expression and genotype data are valuable in the quest for genetic regulatory loci underlying complex traits. Gene coexpression network analysis lends itself to identification of entire groups of differentially regulated genes-a highly relevant endeavor in finding the underpinnings of complex traits that are, by definition, polygenic in nature. Here we describe one such approach based on liver gene expression and genotype data from an F(2) mouse inter-cross utilizing weighted gene coexpression network analysis (WGCNA) of gene expression data to identify physiologically relevant modules. We describe two strategies: single-network analysis and differential network analysis. Single-network analysis reveals the presence of a physiologically interesting module that can be found in two distinct mouse crosses. Module quantitative trait loci (mQTLs) that perturb this module were discovered. In addition, we report a list of genetic drivers for this module. Differential network analysis reveals differences in connectivity and module structure between two networks based on the liver expression data of lean and obese mice. Functional annotation of these genes suggests a biological pathway involving epidermal growth factor (EGF). Our results demonstrate the utility of WGCNA in identifying genetic drivers and in finding genetic pathways represented by gene modules. These examples provide evidence that integration of network properties may well help chart the path across the gene-trait chasm.
| S-EPMC1998880 | biostudies-literature
Project description:AimsHeart failure (HF) is a chronic heart disease with a high incidence and mortality. Due to the regulatory complexity of gene coexpression networks, the underlying hub genes regulation in HF remain incompletely appreciated. We aimed to explore potential key modules and genes for HF using weighted gene coexpression network analysis (WGCNA).Methods and resultsThe expression profiles by high throughput sequencing of heart tissues samples from HF and non-HF samples were obtained from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) between HF and non-HF samples were firstly identified. Then, a coexpression network was constructed to identify key modules and potential hub genes. The biological functions of potential hub genes were analysed by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses. Finally, a protein-protein interaction (PPI) network was constructed using the STRING online tool. A total of 135 DEGs (133 up-regulated and 2 down-regulated DEGs) between HF and non-HF samples were identified in the GSE135055 and GSE123976 datasets. Moreover, a total of 38 modules were screened based on WGCNA in the GSE135055 dataset, and six potential hub genes (UCK2, ASB1, CCNI, CUX1, IRX6, and STX16) were screened from the key module by setting the gene significance over 0.2 and the module membership over 0.8. Furthermore, 78 potential hub genes were obtained by taking the intersection of the 135 DEGs and all genes in the key module, and enrichment analysis revealed that they were mainly involved in the MAPK and PI3K-AKT signalling pathways. Finally, in a PPI network constructed with the 78 potential hub genes, CUX1 and ASB1 were identified as hub genes in HF because they were also identified as potential hub genes in the WGCNA.ConclusionsTo the best of our knowledge, our study is the first to employ WGCNA to identify the key module and hub genes for HF. Our study identified a module and two genes that might play important roles in HF, which may provide potential biomarkers for the diagnosis of HF and improve our knowledge of the molecular mechanisms underlying HF.
| S-EPMC8934958 | biostudies-literature
Project description:Dilated cardiomyopathy (DCM), a heterogeneous cardiomyopathy, is a major cause of heart failure and heart transplant. Currently, immunotherapy is believed to be an effective treatment method for DCM. However, individual differences are so obvious that the clinical effect is not satisfactory. In order to find immune-related biomarkers of DCM to guide treatment and improve clinical efficacy, we downloaded a GSE120895 dataset from the Gene Expression Omnibus (GEO) database using CIBERSORT and WGCNA algorithms in RStudio and visualizing the protein-protein interaction (PPI) network for key modules by Cytoscape, and finally obtained six hub genes. A GSE17800 dataset was downloaded from the GEO dataset to verify the diagnostic values of hub genes, MYG1, FLOT1, and ATG13, which were excellent. Our study revealed unpublished potential immune mechanisms, biomarkers, and therapeutic targets of DCM.
| S-EPMC8950518 | biostudies-literature
Project description:BackgroundOvarian cancer is one of the most lethal diseases of women. The prognosis of ovarian cancer patients was closely correlated with immune cell expression and immune responses. Therefore, it is important to identify a robust prognostic signature, which correlates not only with prognoses but also with immune responses in ovarian cancer, thus, providing immune-related patient therapies.MethodsThe weighted gene coexpression network analysis (WGCNA) was used to identify candidate genes correlated with ovarian cancer prognoses. Univariate and multivariate Cox regression analyses were used to construct the prognostic signature. The Kaplan-Meier method was used to predict survival, and the immune-related bioinformatics analysis was performed using the R software. The relationship between the signature and clinical parameters was analyzed with the GraphPad Prism 7 and SPSS software.ResultsGene expression from The Cancer Genome Atlas dataset was used to perform the WGCNA analysis, and candidate prognostic-related genes in patients with ovarian cancer were identified. According to the Cox regression analysis, the prognostic signature was constructed, which divided patients into two groups. The high-risk group showed the least favorable prognosis. Three independent cohorts from the Gene Expression Omnibus (GEO) database were used for the validation studies. According to the immune analyses, the GEO database signatures were significantly correlated with the immune statuses of ovarian cancer patients. By analyzing the combination of the prognostic signature and total mutational burden (TMB), ovarian cancer patients were divided into four groups. In these groups, memory B cell, resting memory CD4 T cell, M2 macrophage, resting mast cell, and neutrophil were found with significant distinctions among these groups.ConclusionsThis novel signature predicted the prognosis of ovarian cancer patients precisely and independently and showed significant correlations with immune responses. Therefore, this prognostic signature could indicate targeted immunotherapies for ovarian cancer patients.
| S-EPMC7749778 | biostudies-literature
Project description:Symptomatic glycerol kinase deficiency (GKD) is associated with episodic metabolic and central nervous system deterioration. We report here the first application of Weighted Gene Co-Expression Network Analysis (WGCNA) to investigate a knockout (KO) murine model of a human genetic disease. WGCNA identified networks and key hub transcripts from liver mRNA of glycerol kinase (Gyk) KO and wild type (WT) mice. Day of life 1 (dol1) samples from KO mice contained a network module enriched for organic acid metabolism before Gyk KO mice develop organic acidemia and die on dol3-4 and the module containing Gyk was enriched with apoptotic genes. Roles for the highly connected Acot, Psat and Plk3 transcripts were confirmed in cell cultures and subsequently validated by causality testing. We provide evidence that GK may have an apoptotic moonlighting role that is lost in GKD. This systems biology strategy has improved our understanding of GKD pathogenesis and suggests possible treatments. Male WT and KO mouse pups were sacrificed on day of life (dol) 1 and each liver was harvested. Total RNA from 4 KO and 3 WT livers was isolated individually. Affymetrix mus 430 2.0 GeneChips were used to analyze differences in liver gene expression between KO and WT mice. Dol1 and 3 Gyk KO mice represent different disease states. Dol 1 was chosen because mice are phenotypically asymptomatic with respect to Glycerol Kinase Deficiency (GKD) and allowed us to look at alterations that occur before the overt disease state. Dol 3 mice are phenotypically symptomatic with respect to GKD.
2009-09-10 | E-GEOD-12748 | biostudies-arrayexpress
Project description:Symptomatic glycerol kinase deficiency (GKD) is associated with episodic metabolic and central nervous system deterioration. We report here the first application of Weighted Gene Co-Expression Network Analysis (WGCNA) to investigate a knockout (KO) murine model of a human genetic disease. WGCNA identified networks and key hub transcripts from liver mRNA of glycerol kinase (Gyk) KO and wild type (WT) mice. Day of life 1 (dol1) samples from KO mice contained a network module enriched for organic acid metabolism before Gyk KO mice develop organic acidemia and die on dol3-4 and the module containing Gyk was enriched with apoptotic genes. Roles for the highly connected Acot, Psat and Plk3 transcripts were confirmed in cell cultures and subsequently validated by causality testing. We provide evidence that GK may have an apoptotic moonlighting role that is lost in GKD. This systems biology strategy has improved our understanding of GKD pathogenesis and suggests possible treatments.
2009-09-10 | GSE12748 | GEO