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: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:Salvia miltiorrhiza, known as Danshen, has attracted worldwide interest for its substantial effects on coronary heart disease (CHD). Danshensu (DSS) is one of the main active ingredients of Danshen on CHD. Although it has been proven to have a good clinical effect on CHD, the action mechanisms remain elusive. In the current study, a coexpression network-based approach was used to illustrate the beneficial properties of DSS in the context of CHD. By integrating the gene expression profile data and protein-protein interactions (PPIs) data, two coexpression protein interaction networks (CePIN) in a CHD state (CHD CePIN) and a non-CHD state (non-CHD CePIN) were generated. Then, shared nodes and unique nodes in CHD CePIN were attained by conducting a comparison between CHD CePIN and non-CHD CePIN. By calculating the topological parameters of each shared node and unique node in the networks, and comparing the differentially expressed genes, target proteins involved in disease regulation were attained. Then, Gene Ontology (GO) enrichment was utilized to identify biological processes associated to target proteins. Consequently, it turned out that the treatment of CHD with DSS may be partly attributed to the regulation of immunization and blood circulation. Also, it indicated that sodium/hydrogen exchanger 3 (SLC9A3), Prostaglandin G/H synthase 2 (PTGS2), Oxidized low-density lipoprotein receptor 1 (OLR1), and fibrinogen gamma chain (FGG) may be potential therapeutic targets for CHD. In summary, this study provided a novel coexpression protein interaction network approach to provide an explanation of the mechanisms of DSS on CHD and identify key proteins which maybe the potential therapeutic targets for CHD.
| S-EPMC5486119 | 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:Heart development is a complex process requiring dynamic transcriptional regulation. Disturbance of this process will lead to severe developmental defects such as congenital heart disease/defect (CHD). CHD is a group of complex disorder with high genetic heterogeneity, common pathways associated with CHD remains largely unknown. In the manuscript, we focused on the tissue specific genes in human fetal heart samples to explore such pathways. We used the RNA microarray dataset of human fetal tissues from ENCODE project to identify genes with heart tissue specific expression. A transcriptional network was constructed for these genes based on the Pearson correlation coefficients of their expression levels. Function, selective constraints and disease associations of these genes were then examined. Our analysis identified a network consisted of 316 genes with human fetal heart specific expression. The network was highly co-regulated and showed evolutionary conserved tissue expression pattern in tetrapod. Genes in this network are enriched in CHD specific genes and disease mutations. Using the transcriptomic data, we discovered a highly concerted gene network that might reflect a common pathway associated with the etiology of CHD. Such analysis should be helpful for disease associated gene identification in clinical studies.
| S-EPMC5402266 | 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: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
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: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:BackgroundTo date, vascular dementia (VaD) lacks effective treatment in clinical practice. There is also growing evidence that VaD may be closely related to the immune response. The development of high-throughput technology, and the recently discovered group of new mediators called competitive endogenous RNAs (ceRNA), provides a unique opportunity to study the immunomodulation of VaD.MethodsIn this study, we used gene expression profiles in the Gene Expression Omnibus (GEO) database to obtain immune-related gene coexpression modules through a weighted gene coexpression network analysis (WGCNA) and gene enrichment analysis. We extracted and merged long non-coding RNA (lncRNA) and microRNA (miRNA) expressions from the GEO database and mapped them with related databases. Subsequently, we used Cytoscape to construct a lncRNA-miRNA-mRNA network, and then we performed an enrichment analysis on the mRNAs in the network to determine their regulatory function. Subsequently, we used an ImmuCellAI immune infiltration analysis and constructed a ceRNA sub-network of related immune cells. Finally, we conducted a gene set enrichment analysis (GSEA) to determine the potential regulatory pathways of the key factors.ResultsAs a result, we identified the blue module as the key module of immunity and constructed the related CeRNA network. Immune infiltration analysis showed that natural killer T (NKT) cells may be the key immune cells of VaD. Using a Pearson correlation analysis, we identified that B4GALT1, PPP1R3B, MICB, HHAT, DSC2, DNA2, SCARA3, and lncRNA NEAT1 may be the key factors of VaD. Our subsequent GSEA analysis showed that lncRNA NEAT1 may be regulated by NK cells and toll-like receptors.ConclusionsOur research provides new therapeutic targets for vascular dementia from the immunological perspective for the first time, including B4GALT1, PPP1R3B, MICB, HHAT, DSC2, DNA2, SCARA3, and lncRNA NEAT1, and our research hopes to provide new treatment options for VaD.
| S-EPMC8184445 | biostudies-literature
Project description:BackgroundPreeclampsia (PE) is a multisystemic syndrome which has short- and long-term risk to mothers and children and has pluralistic etiology.ObjectiveThis study is aimed at constructing a competitive endogenous RNA (ceRNA) network for pathways most related to PE using a data mining strategy based on weighted gene coexpression network analysis (WGCNA).MethodsWe focused on pathways involving hypoxia, angiogenesis, and epithelial mesenchymal transition according to the gene set variation analysis (GSVA) scores. The gene sets of these three pathways were enriched by gene set enrichment analysis (GSEA). WGCNA was used to study the underlying molecular mechanisms of the three pathways in the pathogenesis of PE by analyzing the relationship among pathways and genes. The soft threshold power (β) and topological overlap matrix allowed us to obtain 15 modules, among which the red module was chosen for the downstream analysis. We chose 10 hub genes that satisfied ∣log2Fold Change | >2 and had a higher degree of connectivity within the module. These candidate genes were subsequently confirmed to have higher gene significance and module membership in the red module. Coexpression networks were established for the hub genes to unfold the connection between the genes in the red module and PE. Finally, ceRNA networks were constructed to further clarify the underlying molecular mechanism involved in the occurrence of PE. 56 circRNAs, 17 lncRNAs, and 20 miRNAs participated in the regulation of the hub genes. Coagulation factor II thrombin receptor (F2R) and lumican (LUM) were considered the most relevant genes, and ceRNA networks of them were constructed.ConclusionThe microarray data mining process based on bioinformatics methods constructed lncRNA and miRNA networks for ten hub genes that were closely related to PE and focused on ceRNAs of F2R and LUM finally. The results of our study may provide insight into the mechanisms underlying PE occurrence.
| S-EPMC8759911 | biostudies-literature