Bladder cancer stage-associated hub genes revealed by WGCNA co-expression network analysis.
Ontology highlight
ABSTRACT:
SUBMITTER: Di Y
PROVIDER: S-EPMC6350372 | biostudies-literature | 2019
REPOSITORIES: biostudies-literature
Publications
Bladder cancer stage-associated hub genes revealed by WGCNA co-expression network analysis.
Di Yu Y Chen Dongshan D Yu Wei W Yan Lei L
Hereditas 20190128
<h4>Background</h4>Bladder cancer was a malignant disease in patients, our research aimed at discovering the possible biomarkers for the diseases.<h4>Results</h4>The gene chip GSE31684, including 93samples, was downloaded from the GEO datasets and co-expression network was constructed by the data. Molecular complex detection(MCODE) was used to identify hub genes. The most significant cluster including 16 genes: <i>CDH11</i>, <i>COL3A1</i>, <i>COL6A3</i>, <i>COL5A1</i>, <i>AEBP1</i>, <i>COL1A2</i ...[more]
Similar Datasets
Project description:The present study aimed to identify the recurrence‑associated genes in colon cancer, which may provide theoretical evidence for the development of novel methods to prevent tumor recurrence. Colon cancer and matched normal samples microarray data (E‑GEOD‑39582) were downloaded from ArrayExpress. Genes with significant variation were identified, followed by the screening of differentially expressed genes (DEGs). Subsequently, the co‑expression network of DEGs was constructed using the weighted correlation network analysis (WGCNA) method, which was verified using the validation dataset. The significant modules associated with recurrence in the network were subsequently screened and verified in another independent dataset E‑GEOD‑33113. Function and pathway enrichment analyses were also conducted to determine the roles of selected genes. Survival analysis was performed to identify the association between these genes and survival. A total of 434 DEGs were identified in the colon samples, and stress‑associated endoplasmic reticulum protein family member 2 (SERP2) and long non‑coding RNA‑0219 (LINC0219) were determined to be the vital DEGs between all the three sub‑type groups with different clinical features. The brown module was identified to be the most significant module in the co‑expression network associated with the recurrence of colon cancer, which was verified in the E‑GEOD‑33113 dataset. Top 10 genes in the brown module, including EGF containing fibulin like extracellular matrix protein 2 (EFEMP2), fibrillin 1 (FBN1) and secreted protein acidic and cysteine rich (SPARC) were also associated with survival time of colon cancer patients. Further analysis revealed that the function of cell adhesion, biological adhesion, extracellular matrix (ECM) organization, pathways of ECM‑receptor interaction and focal adhesion were the significantly changed terms in colon cancer. In conclusion, SERP2, EFEMP2, FBN1, SPARC, and LINC0219 were revealed to be the recurrence‑associated molecular and prognostic indicators in colon cancer by WGCNA co‑expression network analysis.
| S-EPMC5865817 | biostudies-literature
Project description:BackgroundBladder cancer (BC) is among the most frequent cancers globally. Although substantial efforts have been put to understand its pathogenesis, its underlying molecular mechanisms have not been fully elucidated.MethodsThe robust rank aggregation approach was adopted to integrate 4 eligible bladder urothelial carcinoma microarray datasets from the Gene Expression Omnibus. Differentially expressed gene sets were identified between tumor samples and equivalent healthy samples. We constructed gene co-expression networks using weighted gene co-expression network to explore the alleged relationship between BC clinical characteristics and gene sets, as well as to identify hub genes. We also incorporated the weighted gene co-expression network and robust rank aggregation to screen differentially expressed genes.ResultsCDH11, COL6A3, EDNRA, and SERPINF1 were selected from the key module and validated. Based on the results, significant downregulation of the hub genes occurred during the early stages of BC. Moreover, receiver operating characteristics curves and Kaplan-Meier plots showed that the genes exhibited favorable diagnostic and prognostic value for BC. Based on gene set enrichment analysis for single hub gene, all the genes were closely linked to BC cell proliferation.ConclusionsThese results offer unique insight into the pathogenesis of BC and recognize CDH11, COL6A3, EDNRA, and SERPINF1 as potential biomarkers with diagnostic and prognostic roles in BC.
| S-EPMC9794320 | biostudies-literature
Project description:The aim of the present study was to identify hub genes and signaling pathways associated with bladder cancer (BC) utilizing centrality analysis and pathway enrichment analysis. The differentially expressed genes (DEGs) were screened from the ArrayExpress database between normal subjects and BC patients. Co-expression networks of BC were constructed using differentially co-expressed genes and links, and hub genes were investigated by degree centrality analysis of co-expression networks in BC. The enriched signaling pathways were investigated by Kyoto Encyclopedia of Genes and Genomes database analysis based on the DEGs. The hub gene expression in BC tissues was validated using reverse transcription-quantitative polymerase chain reaction (RT-qPCR) and western blotting. A total of 329 DEGs were screened, including 147 upregulated and 182 downregulated genes. The co-expression network constructed between BC and normal controls consisted of 182 nodes and 434 edges, and the two genes in each gene pair were differentially co-expressed genes. Centrality analysis of co-expression networks suggested that the top 5 hub genes with high degree included lectin, galactoside-binding, soluble, 4 (LGALS4), protein tyrosine phosphatase, receptor type N2 (PTPRN2), transmembrane protease, serine 11E (TMPRSS11E), tripartite motif containing 31 (TRIM31) and potassium voltage-gated channel subfamily D member 3 (KCND3). Pathway analysis revealed that the 329 DEGs were significantly enriched in 5 terms (cell cycle, DNA replication, oocyte meiosis, p53 signaling pathway and peroxisome proliferator-activated receptor signaling pathway). According to RT-qPCR and western blot analysis, 4/5 hub genes were significantly expressed, including LGALS4, PTPRN2, TMPRSS11E, TRIM31; however, KCND3 was not significantly expressed. In the present study, 5 hub genes were successfully identified (LGALS4, PTPRN2, TMPRSS11E, TRIM31 and KCND3) and 5 biological pathways that may be underlying biomarkers for early diagnosis and treatment associated with bladder cancer were revealed.
| S-EPMC5494668 | biostudies-other
Project description:Weighted gene coexpression network analysis (WGCNA) is a novel approach that can quickly analyze the relationships between genes and traits. In this study, the milk yield, lactose, fat, and protein of Holstein dairy cows were detected in a lactation cycle. Meanwhile, a total of 18 gene expression profiles were detected using mammary glands from six lactation stages (day 7 to calving, -7 d; day 30 post-calving, 30 d; day 90 post-calving, 90 d; day 180 post-calving, 180 d; day 270 post-calving, 270 d; day 315 post-calving, 315 d). On the basis of the 18 profiles, WGCNA identified for the first time 10 significant modules that may be related to lactation stage, milk yield, and the main milk composition content. Genes in the 10 significant modules were examined with gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The results revealed that the galactose metabolism pathway was a potential candidate for milk yield and milk lactose synthesis. In -7 d, ion transportation was more frequent and cell proliferation related terms became active. In late lactation, the suppressor of cytokine signaling 3 (SOCS3) might play a role in apoptosis. The sphingolipid signaling pathway was a potential candidate for milk fat synthesis. Dairy cows at 315 d were in a period of cell proliferation. Another notable phenomenon was that nonlactating dairy cows had a more regular circadian rhythm after a cycle of lactation. The results provide an important theoretical basis for the further molecular breeding of dairy cows.
| S-EPMC7911360 | biostudies-literature
Project description:Osteosarcoma (OS), the most common malignant bone tumor, accounts for the heavy healthy threat in the period of children and adolescents. OS occurrence usually correlates with early metastasis and high death rate. This study aimed to better understand the mechanism of OS metastasis.Based on Gene Expression Omnibus (GEO) database, we downloaded 4 expression profile data sets associated with OS metastasis, and selected differential expressed genes. Weighted gene co-expression network analysis (WGCNA) approach allowed us to investigate the most OS metastasis-correlated module. Gene Ontology functional and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were used to give annotation of selected OS metastasis-associated genes.We select 897 differential expressed genes from OS metastasis and OS non-metastasis groups. Based on these selected genes, WGCNA further explored 142 genes included in the most OS metastasis-correlated module. Gene Ontology functional and KEGG pathway enrichment analyses showed that significantly OS metastasis-associated genes were involved in pathway correlated with insulin-like growth factor binding.Our research figured out several potential molecules participating in metastasis process and factors acting as biomarker. With this study, we could better explore the mechanism of OS metastasis and further discover more therapy targets.
| S-EPMC6023727 | biostudies-literature
Project description:ObjectiveMature hair follicles represent an important stage of hair follicle development, which determines the stability of hair follicle structure and its ability to enter the hair cycle. Here, we used weighted gene co-expression network analysis (WGCNA) to identify hub genes of mature skin and hair follicles in Inner Mongolian cashmere goats.MethodsWe used transcriptome sequencing data for the skin of Inner Mongolian cashmere goats from fetal days 45-135 days, and divided the co expressed genes into different modules by WGCNA. Characteristic values were used to screen out modules that were highly expressed in mature skin follicles. Module hub genes were then selected based on the correlation coefficients between the gene and module eigenvalue, gene connectivity, and Gene Ontology (GO)/Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. The results were confirmed by quantitative polymerase chain reaction (qPCR).ResultsTen modules were successfully defined, of which one, with a total of 3166 genes, was selected as a specific module through sample and gene expression pattern analyses. A total of 584 candidate hub genes in the module were screened by the correlation coefficients between the genes and module eigenvalue and gene connectivity. Finally, GO/KEGG functional enrichment analyses detected WNT10A as a key gene in the development and maturation of skin hair follicles in fetal Inner Mongolian cashmere goats. qPCR showed that the expression trends of 13 genes from seven fetal skin samples were consistent with the sequencing results, indicating that the sequencing results were reliable.n.
| S-EPMC7755285 | biostudies-literature
Project description:BackgroundHomologous repair deficiency (HRD) causes double-strand break repair to be impeded, which is a common driver of carcinogenesis. However, the therapeutic and prognostic potential of HRD in invasive breast cancer (BRCA) has not been fully explored using comprehensive bioinformatics analysis.Materials and methodsHRD score was defined as the unweighted sum of LOH, TAI, and LST scores and obtained from the previous study according to Theo A et al. To characterize BRCA tumor microenvironment (TME) subtypes, "ConsensusClusterPlus" R package was used to conduct unsupervised clustering. The xCell algorithm was utilized for tumor composition analysis to estimate the TME in TCGA-BRCA. A WGCNA analysis was conducted to uncover the gene coexpression modules and hub genes in the HRD-related gene module of BRCA. The functional enrichment study was carried out using Metascape. A novel analysis pipeline, Genetic Perturbation Similarity Analysis (GPSA), was used to explore the single-gene perturbation closely related to HRD based on 3048 stable knockdown/knockout cell lines. The prognostic variables were evaluated using univariate COX analysis. Kaplan-Meier (KM) survival analysis was performed to assess the prognostic potential of HRD score. Receiver operator characteristic (ROC) curve was utilized to judge the diagnostic utility. Drug sensitivity was estimated through the R package "oncoPredict" and Genomics of Drug Sensitivity in Cancer (GDSC) database. XSum algorithm was performed to screen the candidate small molecule drugs based on the connectivity map (CMAP) database.ResultsLow HRD score suggested a better prognosis in BRCA patients. The tumor with low HRD score had considerably greater degree of infiltration of stromal cells and infiltration of immunocytes was significantly enhanced in the high HRD score group. Using WGCNA, ten co-expression modules were obtained. The turquoise module and 25 hub genes were identified as the most correlated with HRD in BRCA. Functional enrichment analysis revealed that the turquoise gene module was mainly concentrated in the "cell cycle" pathways. Candidate HRD-related gene signatures (MELK) were screened out through WGCNA and GPSA analysis pipeline and then validated on independent validation sets. A small molecule drug (Clofibrate) that has the potential to reverse the increase of high HRD score was screened out to improve oncological outcomes in BRCA. Molecular docking suggested MELK to be one of possible molecular targets in the Clofibrate treatment of BRCA.ConclusionBased on bioinformatic analysis, we fully explored the therapeutic and prognostic potential of HRD in BRCA. A novel HRD-related gene signature (MELK) were identified through the combination of WGCNA and GPSA analysis. In addition, we detailed the TME landscape in BRCA and identified four unique TME subtypes in group with high or low HRD score group. Moreover, Clofibrate were screened out to improve oncological outcomes in BRCA by reversing the increase of high HRD score. Thus, our study contributes to the development of personalized clinical management and treatment regimens of BRCA.
| S-EPMC10081085 | biostudies-literature
Project description:BackgroundMuscular invasive bladder cancer (MIBC) is a common malignant tumor in the world. Because of their heterogeneity in prognosis and response to treatment, biomarkers that can predict survival or help make treatment decisions in patients with MIBC are essential for individualized treatment.AimWe aimed to integrate bioinformatics research methods to identify a set of effective biomarkers capable of predicting, diagnosing, and treating MIBC. To provide a new theoretical basis for the diagnosis and treatment of bladder cancer.Methods and resultsGene expression profiles and clinical data of MIBC were obtained by downloading from the Cancer Genome Atlas database. A dataset of 129 MIBC cases and controls was included. 2084 up-regulated genes and 2961 down-regulated genes were identified by differentially expressed gene (DEG) analysis. Then, gene ontology analysis was performed to explore the biological functions of DEGs, respectively. The up-regulated DEGs are mainly enriched in epidermal cell differentiation, mitotic nuclear division, and so forth. They are also involved in the cell cycle, p53 signaling pathway, PPAR signaling pathway, and so forth. The weighted gene co-expression network analysis yielded five modules related to pathological stages and grading, of which blue and turquoise were the most relevant modules for MIBC. Next, Using Kaplan-Meier survival analysis to identify further hub genes, the screening criteria at p ≤ .05, we found CNKSR1, HIP1R, CFL2, TPM1, CSRP1, SYNM, POPDC2, PJA2, and RBBP8NL genes associated with the progression and prognosis of MIBC patients. Finally, immunohistochemistry experiments further confirmed that CNKSR1 plays a vital role in the tumorigenic context of MIBC.ConclusionThe research suggests that CNKSR1, POPDC2, and PJA2 may be novel biomarkers as therapeutic targets for MIBC, especially we used immunohistochemical further to validate CNKSR1 as a therapeutic target for MIBC which may help to improve the prognosis for MIBC.
| S-EPMC9458504 | biostudies-literature
Project description:Weighted gene co-expression network analysis (WGCNA) is a research method in systematic biology. It is widely used to identify gene modules related to target traits in multi-sample transcriptome data. In order to further explore the molecular mechanism of maize response to low-temperature stress at the seedling stage, B144 (cold stress tolerant) and Q319 (cold stress sensitive) provided by the Maize Research Institute of Heilongjiang Academy of Agricultural Sciences were used as experimental materials, and both inbred lines were treated with 5 °C for 0 h, 12 h, and 24 h, with the untreated material as a control. Eighteen leaf samples were used for transcriptome sequencing, with three biological replicates. Based on the above transcriptome data, co-expression networks of weighted genes associated with low-temperature-tolerance traits were constructed by WGCNA. Twelve gene modules significantly related to low-temperature tolerance at the seedling stage were obtained, and a number of hub genes involved in low-temperature stress regulation pathways were discovered from the four modules with the highest correlation with target traits. These results provide clues for further study on the molecular genetic mechanisms of low-temperature tolerance in maize at the seedling stage.
| S-EPMC10454290 | biostudies-literature
Project description:BackgroundAcute mountain sickness has become a heavily researched topic in recent years. However, the genetic mechanism and effects have not been elucidated. Our goal is to construct a gene co-expression network to identify the key modules and hub genes associated with high altitude hypoxia.ResultsThe GSE46480 dataset of rapidly transported healthy adults with acute mountain sickness was selected and analyzed by weighted gene co-expression network analysis (WGCNA) to construct a co-expression network. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the data set were carried out using Database for Annotation Visualization and Integrated Discovery (DAVID), and the hub genes were selected. We found that the turquoise module was most significantly correlated with acute mountain sickness. The functional enrichment analysis showed that the turquoise module was related to the apoptotic process, protein transport, and translation processes. The metabolic pathway analysis identified hsa03010:ribosome and hsa04144:endocytosis as the most important pathways in the turquoise module. Ten top 10 hub genes (MRPL3, PSMC6, AIMP1, HAT1, DPY30, ATP5L, COX7B, UQCRB, DPM1, and COMMD6) for acute mountain sickness were identified.ConclusionOne module and 10 hub genes were identified, which were related to acute mountain sickness. The reference provided by this module may help to elucidate the mechanism of acute mountain sickness. In addition, the hub genes may be used in the future as a biomarker and therapeutic target for accurate diagnosis and treatment.
| S-EPMC7164164 | biostudies-literature