Correlation signature of the macroscopic states of the gene regulatory network in cancer.
Ontology highlight
ABSTRACT:
SUBMITTER: Slavov N
PROVIDER: S-EPMC2657377 | biostudies-literature | 2009 Mar
REPOSITORIES: biostudies-literature
Publications
Correlation signature of the macroscopic states of the gene regulatory network in cancer.
Slavov Nikolai N Dawson Kenneth A KA
Proceedings of the National Academy of Sciences of the United States of America 20090225 11
Although cancer types differ substantially, many cancers share common gene expression signatures. Consistent with this observation, we find convergent and representative distributions and correlation vectors that are distinct in cancer and noncancer ensembles. These differences originate in many genes, but comparatively few genes account for the major differences. We identify genes with different combinatorial regulation in cancer and noncancer as indicated by significant differences in their co ...[more]
Similar Datasets
Project description:Endometrial hyperplasia (EH) is a precursor for endometrial cancer (EC). However, biomarkers for the progression from EH to EC and standard prognostic biomarkers for EC have not been identified. In this study, we aimed to identify key genes with prognostic significance for the progression from EH to EC. Weighted-gene correlation network analysis (WGCNA) was used to identify hub genes utilizing microarray data (GSE106191) downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were identified from the Uterine Corpus Endometrial Carcinoma (UCEC) dataset of The Cancer Genome Atlas database. The Limma-Voom R package was applied to detect differentially expressed genes (DEGs; mRNAs) between cancer and normal samples. Genes with |log2 (fold change [FC])| > 1.0 and p < 0.05 were considered as DEGs. Univariate and multivariate Cox regression and survival analyses were performed to identify potential prognostic genes using hub genes overlapping in the two datasets. All analyses were conducted using R Bioconductor and related packages. Through WGCNA and overlapping genes in hub modules with DEGs in the UCEC dataset, we identified 42 hub genes. The results of the univariate and multivariate Cox regression analyses revealed that four hub genes, BUB1B, NDC80, TPX2, and TTK, were independently associated with the prognosis of EC (Hazard ratio [95% confidence interval]: 0.591 [0.382-0.912], p = 0.017; 0.605 [0.371-0.986], p = 0.044; 1.678 [1.132-2.488], p = 0.01; 2.428 [1.372-4.29], p = 0.02, respectively). A nomogram was established with a risk score calculated using the four genes' coefficients in the multivariate analysis, and tumor grade and stage had a favorable predictive value for the prognosis of EC. The survival analysis showed that the high-risk group had an unfavorable prognosis compared with the low-risk group (p < 0.0001). The receiver operating characteristic curves also indicated that the risk model had a potential predictive value of prognosis with area under the curve 0.807 at 2 years, 0.783 at 3 years, and 0.786 at 5 years. We established a four-gene signature with prognostic significance in EC using WGCNA and established a nomogram to predict the prognosis of EC.
| S-EPMC8488359 | biostudies-literature
Project description:Background and objectivesSepsis is a leading cause of mortality and morbidity in the intensive care unit. Regulatory mechanisms underlying the disease progression and prognosis are largely unknown. The study aimed to identify master regulators of mortality-related modules, providing potential therapeutic target for further translational experiments.MethodsThe dataset GSE65682 from the Gene Expression Omnibus (GEO) database was utilized for bioinformatic analysis. Consensus weighted gene co-expression netwoek analysis (WGCNA) was performed to identify modules of sepsis. The module most significantly associated with mortality were further analyzed for the identification of master regulators of transcription factors and miRNA.ResultsA total number of 682 subjects with various causes of sepsis were included for consensus WGCNA analysis, which identified 27 modules. The network was well preserved among different causes of sepsis. Two modules designated as black and light yellow module were found to be associated with mortality outcome. Key regulators of the black and light yellow modules were the transcription factor CEBPB (normalized enrichment score = 5.53) and ETV6 (NES = 6), respectively. The top 5 miRNA regulated the most number of genes were hsa-miR-335-5p (n = 59), hsa-miR-26b-5p (n = 57), hsa-miR-16-5p (n = 44), hsa-miR-17-5p (n = 42), and hsa-miR-124-3p (n = 38). Clustering analysis in 2-dimension space derived from manifold learning identified two subclasses of sepsis, which showed significant association with survival in Cox proportional hazard model (p = 0.018).ConclusionsThe present study showed that the black and light-yellow modules were significantly associated with mortality outcome. Master regulators of the module included transcription factor CEBPB and ETV6. miRNA-target interactions identified significantly enriched miRNA.
| S-EPMC7545567 | biostudies-literature
Project description:Gene expression profile or transcriptome can represent cellular states, thus understanding gene regulation mechanisms can help understand how cells respond to external stress. Interaction between transcription factor (TF) and target gene (TG) is one of the representative regulatory mechanisms in cells. In this paper, we present a novel computational method to construct condition-specific transcriptional networks from transcriptome data. Regulatory interaction between TFs and TGs is very complex, specifically multiple-to-multiple relations. Experimental data from TF Chromatin Immunoprecipitation sequencing is useful but produces one-to-multiple relations between TF and TGs. On the other hand, co-expression networks of genes can be useful for constructing condition transcriptional networks, but there are many false positive relations in co-expression networks. In this paper, we propose a novel method to construct a condition-specific and combinatorial transcriptional network, applying kernel canonical correlation analysis (kernel CCA) to identify multiple-to-multiple TF-TG relations in certain biological condition. Kernel CCA is a well-established statistical method for computing the correlation of a group of features vs. another group of features. We, therefore, employed kernel CCA to embed TFs and TGs into a new space where the correlation of TFs and TGs are reflected. To demonstrate the usefulness of our network construction method, we used the blood transcriptome data for the investigation on the response to high fat diet in a human and an arabidopsis data set for the investigation on the response to cold/heat stress. Our method detected not only important regulatory interactions reported in previous studies but also novel TF-TG relations where a module of TF is regulating a module of TGs upon specific stress.
| S-EPMC8172963 | biostudies-literature
Project description:BackgroundIdiopathic pulmonary fibrosis (IPF) is a severe lung disease characterised by extensive pathological changes. The objective for this study was to identify the gene network and regulators underlying disease pathology in IPF and its association with lung function.MethodsLung Tissue Research Consortium dataset with 262 IPF and control subjects (GSE47460) was randomly divided into two non-overlapping groups for cross-validated differential gene expression analysis. Consensus weighted gene coexpression network analysis identified overlapping coexpressed gene modules between both IPF groups. Modules were correlated with lung function (diffusion capacity, DLCO; forced expiratory volume in 1 s, FEV1; forced vital capacity, FVC) and enrichment analyses used to identify biological function and transcription factors. Module correlation with miRNA data (GSE72967) identified associated regulators. Clinical relevance in IPF was assessed in a peripheral blood gene expression dataset (GSE93606) to identify modules related to survival.ResultsCorrelation network analysis identified 16 modules in IPF. Upregulated modules were associated with cilia, DNA replication and repair, contractile fibres, B-cell and unfolded protein response, and extracellular matrix. Downregulated modules were associated with blood vessels, T-cell and interferon responses, leucocyte activation and degranulation, surfactant metabolism, and cellular metabolic and catabolic processes. Lung function correlated with nine modules (eight with DLCO, five with FVC). Intermodular network of transcription factors and miRNA showed clustering of fibrosis, immune response and contractile modules. The cilia-associated module was able to predict survival (p=0.0097) in an independent peripheral blood IPF cohort.ConclusionsWe identified a correlation gene expression network with associated regulators in IPF that provides novel insight into the pathological process of this disease.
| S-EPMC6467239 | biostudies-literature
Project description:BackgroundSciatic nerve injury (SNI), which frequently occurs under the traumatic hip and hip fracture dislocation, induces serious complications such as motor and sensory loss, muscle atrophy, or even disabling. The present work aimed to determine the regulating factors and gene network related to the SNI pathology.MethodsSciatic nerve injury dataset GSE18803 with 24 samples was divided into adult group and neonate group. Weighted gene co-expression network analysis (WGCNA) was carried out to identify modules associated with SNI in the two groups. Moreover, differentially expressed genes (DEGs) were determined from every group, separately. Subsequently, co-expression network and protein-protein interaction (PPI) network were overlapped to identify hub genes, while functional enrichment and Reactome analysis were used for a comprehensive analysis of potential pathways. GSE30165 was used as the test set for investigating the hub gene involvement within SNI. Gene set enrichment analysis (GSEA) was performed separately using difference between samples and gene expression level as phenotype label to further prove SNI-related signaling pathways. In addition, immune infiltration analysis was accomplished by CIBERSORT. Finally, Drug-Gene Interaction database (DGIdb) was employed for predicting the possible therapeutic agents.Results14 SNI status modules and 97 DEGs were identified in adult group, while 15 modules and 21 DEGs in neonate group. A total of 12 hub genes was overlapping from co-expression and PPI network. After the results from both test and training sets were overlapped, we verified that the ten real hub genes showed remarkably up-regulation within SNI. According to functional enrichment of hub genes, the above genes participated in the immune effector process, inflammatory responses, the antigen processing and presentation, and the phagocytosis. GSEA also supported that gene sets with the highest significance were mostly related to the cytokine-cytokine receptor interaction. Analysis of hub genes possible related signaling pathways using gene expression level as phenotype label revealed an enrichment involved in Lysosome, Chemokine signaling pathway, and Neurotrophin signaling pathway. Immune infiltration analysis showed that Macrophages M2 and Regulatory T cells may participate in the development of SNI. At last, 25 drugs were screened from DGIdb to improve SNI treatment.ConclusionsThe gene expression network is determined in the present work based on the related regulating factors within SNI, which sheds more light on SNI pathology and offers the possible biomarkers and therapeutic targets in subsequent research.
| S-EPMC8522103 | biostudies-literature
Project description:Inferring regulatory relationships among many genes based on their temporal variation in transcript abundance has been a popular research topic. Due to the nature of microarray experiments, classical tools for time series analysis lose power since the number of variables far exceeds the number of the samples. In this paper, we describe some of the existing multivariate inference techniques that are applicable to hundreds of variables and show the potential challenges for small-sample, large-scale data. We propose a directed partial correlation (DPC) method as an efficient and effective solution to regulatory network inference using these data. Specifically for genomic data, the proposed method is designed to deal with large-scale datasets. It combines the efficiency of partial correlation for setting up network topology by testing conditional independence, and the concept of Granger causality to assess topology change with induced interruptions. The idea is that when a transcription factor is induced artificially within a gene network, the disruption of the network by the induction signifies a genes role in transcriptional regulation. The benchmarking results using GeneNetWeaver, the simulator for the DREAM challenges, provide strong evidence of the outstanding performance of the proposed DPC method. When applied to real biological data, the inferred starch metabolism network in Arabidopsis reveals many biologically meaningful network modules worthy of further investigation. These results collectively suggest DPC is a versatile tool for genomics research. The R package DPC is available for download (http://code.google.com/p/dpcnet/).
| S-EPMC3071805 | biostudies-literature
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:IntroductionThis bioinformatic study confirmed a new miRNA-mRNA regulatory network and a prognostic signature in endometrial cancer (EC).Materials and methodsWe downloaded RNA-seq and miRNA-seq data of EC from the TCGA database, then used EdegR package to screen differentially expressed miRNAs and mRNAs (DE-miRNAs and DE-mRNAs). Then, we constructed a regulatory network of EC-associated miRNAs and hub genes by Cytoscape, and determined the expression of unexplored miRNAs in EC tissues and normal adjacent tissues by quantitative Real-Time PCR (qRT-PCR). A prognostic signature model and a predictive nomogram were constructed. Finally, we explored the association between the prognostic model and the immune cell infiltration.ResultsA total of 11,531 DE-mRNAs and 236 DE-miRNAs, as well as 275 and 118 candidate DEGs for upregulated and downregulated DE-miRNAs were screened out. The miRNA-mRNA network included 5 downregulated and 13 upregulated DE-miRNAs. qRT-PCR proved that the expression levels of miRNA-18a-5p, miRNA-18b-5p, miRNA-449c-5p and miRNA-1224-5p and their target genes (NR3C1, CTGF, MYC, and TNS1) were consistent with our predictions. Univariate and multivariate Cox proportional hazards regression analyses of the hub genes revealed a significant prognostic value of NR3C1, EZH2, AND GATA4, and these genes were closely related to eight types of immune infiltration cells.ConclusionWe identified three genes as candidate biomarkers for EC, which may provide a theoretical basis for targeted therapy.
| S-EPMC8039850 | biostudies-literature
Project description:Changes in iron regulation characterize the malignant state. However, the pathways that effect these changes and their specific impact on prognosis remain poorly understood. We capitalized on publicly available microarray datasets comprising 674 breast cancer cases to systematically investigate how expression of genes related to iron metabolism is linked to breast cancer prognosis. Of 61 genes involved in iron regulation, 49% were statistically significantly associated with distant metastasis-free survival. Cases were divided into test and training cohorts, and the supervised principal component method was used to stratify cases into risk groups. Optimal risk stratification was achieved with a model comprising 16 genes, which we term the iron regulatory gene signature (IRGS). Multivariable analysis revealed that the IRGS contributes information not captured by conventional prognostic indicators (HR = 1.61; 95% confidence interval: 1.16-2.24; P = 0.004). The IRGS successfully stratified homogeneously treated patients, including ER+ patients treated with tamoxifen monotherapy, both with (P = 0.006) and without (P = 0.03) lymph node metastases. To test whether multiple pathways were embedded within the IRGS, we evaluated the performance of two gene dyads with known roles in iron biology in ER+ patients treated with tamoxifen monotherapy (n = 371). For both dyads, gene combinations that minimized intracellular iron content [anti-import: TFRC(Low)/HFE(High); or pro-export: SLC40A1 (ferroportin)(High)/HAMP(Low)] were associated with favorable prognosis (P < 0.005). Although the clinical utility of the IRGS will require further evaluation, its ability to both identify high-risk patients within traditionally low-risk groups and low-risk patients within high-risk groups has the potential to affect therapeutic decision making.
| S-EPMC3206152 | biostudies-literature
Project description:In this study, we infer the breast cancer gene regulatory network from gene expression data. This network is obtained from the application of the BC3Net inference algorithm to a large-scale gene expression data set consisting of 351 patient samples. In order to elucidate the functional relevance of the inferred network, we are performing a Gene Ontology (GO) analysis for its structural components. Our analysis reveals that most significant GO-terms we find for the breast cancer network represent functional modules of biological processes that are described by known cancer hallmarks, including translation, immune response, cell cycle, organelle fission, mitosis, cell adhesion, RNA processing, RNA splicing and response to wounding. Furthermore, by using a curated list of census cancer genes, we find an enrichment in these functional modules. Finally, we study cooperative effects of chromosomes based on information of interacting genes in the beast cancer network. We find that chromosome 21 is most coactive with other chromosomes. To our knowledge this is the first study investigating the genome-scale breast cancer network.
| S-EPMC3909882 | biostudies-literature