Identification of co-expression modules and pathways correlated with osteosarcoma and its metastasis.
Ontology highlight
ABSTRACT: Background
Osteosarcoma is the most common bone tumor that occurs in children.Methods
To identify co-expression modules and pathways correlated with osteosarcoma and its clinical characteristics, we performed weighted gene co-expression network analysis (WGCNA) on RNA-seq data of osteosarcoma with 52 samples. Then we performed pathway enrichment analysis on genes from significant modules.Results
A total of 5471 genes were included in WGCNA, and 16 modules were identified. Module-trait analysis identified that a module involved in microtubule bundle formation, drug metabolism-cytochrome P450, and IL-17 signaling pathway was negatively correlated with osteosarcoma and positively correlated with metastasis; a module involved in DNA replication was positively correlated with osteosarcoma; a module involved in cell junction was positively correlated with metastasis; and a module involved in heparin binding negatively correlated with osteosarcoma. Moreover, expression levels in four of the top ten differentially expressed genes were validated in another independent dataset.Conclusions
Our analysis might provide insight for molecular mechanisms of osteosarcoma.
SUBMITTER: Wang JS
PROVIDER: S-EPMC6408756 | biostudies-literature | 2019 Mar
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BACKGROUNDOsteosarcoma is a type of bone cancer casting huge threat to the human health worldwide. Previously, gene expression analyses were performed to identify biomarkers for cancer; however, systemic co-expression analysis for osteosarcoma is still in need. The aim of this study was to construct a gene co-expression network that predicts clusters of candidate genes associated with the pathogenesis of osteosarcoma.METHODSHere, we extracted the large scale of datasets from the GEO database. With systematical approaches, we identified the co-expression modules by using weighted gene co-expression network analysis (WGCNA) and investigated the functional enrichments of important modules at GO and KEGG terms.RESULTSFirst, seven co-expression modules, which contain different genes, were conducted for 2228 genes in the 22 human osteosarcoma samples. Then, correlation study showed that the hub genes between pairwise modules displayed great differences. Lastly, functional enrichments of the co-expression modules showed that the module 5 enriched in immune response, antigen processing, and presentation, which is in consistence with GO result. Therefore, we speculated that the module 5 may play a key role in the pathogenesis of osteosarcoma.CONCLUSIONSHere, we speculated that genes of the module 5 were the essential genes that were associated to human osteosarcoma. Together, our findings not only provided outline of co-expression gene modules for human osteosarcoma, but also promoted the understanding of these modules at functional aspects.
| S-EPMC5932805 | biostudies-literature
Project description:BackgroundOsteosarcoma (OS) is a primary malignant bone tumor. Patients with different immune characteristics respond differently to chemotherapy and have a lower chance of survival. The potential pathogenesis and therapeutic targets of OS must be investigated further.MethodsOS expression profile data and clinical information were downloaded from the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) and the Gene Expression Omnibus (GEO) databases. The immune-related gene set was obtained from the ImmPort database, and the immune-related gene expression profiles were used for non-negative matrix factorization (NMF) cluster analysis of patients in the 2 databases to find the best clustering number. In the TARGET database, OS patients were classified into low-risk and high-risk groups based on the differences in their survival rates. Weighted correlation network analysis (WGCNA) was applied to the low-risk and high-risk groups to determine the module with the lowest conservatism in order to differentiate the prognosis of the 2 groups.ResultsA total of 500 key genes associated with poor prognosis were identified. Gene Ontology (GO) enrichment analysis revealed that the biological processes of these genes were primarily focused on the regulation of small guanosine triphosphatase (GTPase) mediated signal transduction, collagen-containing extracellular matrix, and Rho GTPase binding. A random survival forest identified EPHB3, TEAD1, and KRR1P1 as key genes. Their expression level was linked to overall survival. We discovered that the core genes were associated to immune cell infiltration. Simultaneously, paired survival analysis of two genes revealed differences in survival. We also reverse-predicted the main genes and developed their competitive endogenous RNA (ceRNA) network. Finally, utilizing the CellMiner database, we observed that the genes TEAD1 and EPHB3 were connected to drug sensitivity.ConclusionsIn this study, we identified the modules and key genes related to the poor prognosis of OS patients by using WGCNA, and verified their impact on the prognosis of OS patients and their role in the immune microenvironment of OS. In addition, targeted gene related antitumor drugs were screened out. The discoveries may lead to novel molecular targets and treatment methods for OS patients.KeywordsOsteosarcoma (OS); weighted gene co-expression network analysis (WGCNA); gene.
| S-EPMC9011312 | biostudies-literature
Project description:Background: Osteosarcoma (OS) is the most prevalent malignant primary bone tumor in children. Selectin P ligand gene (SELPLG) has been studied in several cancers. Our research aimed to explore the role of SELPLG in OS. Methods: All OS patient data was obtained from TARGET and GEO databases. Differential expression analyses were conducted in limma package of R. Functional analyses included GO and KEGG enrichment analyses. Immune cell infiltration analysis was done in CIBERSORT software. The overall survival was calculated using survival and survminer package of R. Results: Significantly lower SELPLG expression was observed in metastatic OS samples compared with non-metastatic OS samples, both in TARGET and in GSE21257. Low SELPLG expression was an independent undesirable prognostic factor for OS patients, in both TARGET and GEO datasets. Totally 62 differentially expressed gene (DEG) overlaps were found between high SELPLG vs. low SELPLG and non-metastatic vs. metastatic OS samples, affecting metastases and thereby influencing the prognosis, which were significantly enriched in 40 GO and six KEGG terms. Five types of immune cells were significantly differentially infiltrated between high and low SELPLG expression OS patients. Conclusion: SELPLG is closely correlated with metastases and prognosis of OS patients. The OS patients with low SELPLG expression have relatively poorer prognosis and SELPLG is a potential prognostic biomarker for OS.
| S-EPMC8825369 | biostudies-literature
Project description:Osteosarcoma (OS) is the most common histological type of primary bone cancer. The present study was designed to identify the key genes and signaling pathways involved in the metastasis of OS. Microarray data of GSE39055 were downloaded from the Gene Expression Omnibus database, which included 19 OS biopsy specimens before metastasis (control group) and 18 OS biopsy specimens after metastasis (case group). After the differentially expressed genes (DEGs) were identified using the Linear Models for Microarray Analysis package, hierarchical clustering analysis and unsupervised clustering analysis were performed separately, using orange software and the self-organization map method. Based upon the Database for Annotation, Visualization and Integrated Discovery tool and Cytoscape software, enrichment analysis and protein-protein interaction (PPI) network analysis were conducted, respectively. After function deviation scores were calculated for the significantly enriched terms, hierarchical clustering analysis was performed using Cluster 3.0 software. Furthermore, logistic regression analysis was used to identify the terms that were significantly different. Those terms that were significantly different were validated using other independent datasets. There were 840 DEGs in the case group. There were various interactions in the PPI network [including intercellular adhesion molecule-1 (ICAM1), transforming growth factor β1 (TGFB1), TGFB1-platelet-derived growth factor subunit B (PDGFB) and PDGFB-platelet‑derived growth factor receptor-β (PDGFRB)]. Regulation of cell migration, nucleotide excision repair, the Wnt signaling pathway and cell migration were identified as the terms that were significantly different. ICAM1, PDGFB, PDGFRB and TGFB1 were identified to be enriched in cell migration and regulation of cell migration. Nucleotide excision repair and the Wnt signaling pathway were the metastasis-associated pathways of OS. In addition, ICAM1, PDGFB, PDGFRB and TGFB1, which were involved in cell migration and regulation of cell migration may affect the metastasis of OS.
| S-EPMC5819903 | 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:Acute myocardial infarction (AMI) seriously threatens human life. In this study we aimed to systemically analyze the function of key gene modules in human platelets in AMI. We used weighted gene co-expression network analysis (WGCNA) to construct a co-expression module, and analyzed the relationship between potential modules and clinical characteristics based on platelet RNA-seq RPKM count reads of 16 ST-segment elevation myocardial infarction (STEMI) patients and 16 non-STEMI (NSTEMI) patients provided by the GEO database. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed with the DAVID tool. Hub genes were calculated by the Cytohubba package. A total of 3653 genes was selected to construct the co-expression modules. A significant correlation between BMI and the module with color of sky-blue in STEMI. In NSTEMI, there was a significant correlation between the sky blue module and CAD, the Salmon module and HT, and the Cyan module and HT. In STEMI, the Hub genes were mainly enriched in functions related to cell membrane signal transduction including Aqp1, Armcx1, Gsta4, Hist3h2a and Il17re. In NSTEMI, the Hub genes are related mainly to energy metabolism in the sky-blue module including Olr1, Nap1l3, Gfer, Dohh, Crispld1 and Ccdc8b; they are mainly related to extracellular space and calcium binding in the Cyan module, including Clec12b, Chd4, Asgr1, Armcx4, Chid1 and Alkbh7. The hub genes in the Salmon module include Ell3, Aldh1b1, Cavin4, Cabp4, Eif1ay and Dus3l. Our results provide a framework for co-expression gene modules in STEMI and NSTEMI patients, and identify key targets as biomarkers for patients with different subtypes of AMI.
| S-EPMC8129354 | biostudies-literature
Project description:Osteosarcoma is the most common malignant bone tumor in adolescents and young adults, and identifying biomarkers for prognosis and therapy is necessary. Bone morphogenetic protein receptor 2 (BMPR2) is involved in various cellular functions, including cell adhesion, proliferation and invasion, inflammation, apoptosis and metastatic spread. However, the correlation between BMPR2 expression levels and prognosis and tumor-infiltrating immune cells in osteosarcoma is not well understood. In the present study, the expression level of BMPR2 was investigated using the Oncomine and R2 databases. The association between the expression level of BMPR2 and the clinical prognosis of patients with cancer was analyzed using the R2 database. The relationship between the expression level of BMPR2 and immune cell infiltration in the stroma of osteosarcoma was assessed using the Tumor Immune Estimation Resource (TIMER) and CIBERSORT. The correlations between BMPR2 expression level and infiltrated immune cell gene marker sets in osteosarcoma were validated in the TIMER and R2 databases. Analysis of a cohort of patients with osteosarcoma revealed that BMPR2 expression was significantly higher in osteosarcoma compared with in normal tissue and was correlated with poor prognosis. M0 macrophages, M2 macrophages, resting mast, γ δ T and CD8+ T cells were the top five immune cells with the highest degrees of infiltration in osteosarcoma. In addition, BMPR2 expression level showed a significant negative correlation with the gene markers of CD8+ T cells, monocytes and M2 macrophages. Low levels of infiltrating CD8+ T cells, monocytes or M2 macrophages in osteosarcoma was significantly associated with poor survival. These data suggested that CD8+ T cells, monocytes and M2 macrophages play significant roles in the establishment of the immune microenvironment of osteosarcoma. High BMPR2 expression was associated with poor prognosis and low infiltration of CD8+ T cells, monocytes and M2 macrophages in osteosarcoma. Hence, BMPR2 can be considered a biomarker of the immune infiltration, metastasis and prognosis of osteosarcoma.
| S-EPMC7988701 | biostudies-literature
Project description:The rapid accumulation of cancer-related data owing to high-throughput technologies has provided unprecedented choices to understand the progression of cancer and discover functional networks in multiple cancers. Establishment of co-expression networks will help us to discover the systemic properties of carcinogenesis features and regulatory mechanisms of multiple cancers. Here, we proposed a computational workflow to identify differentially co-expressed gene modules across 8 cancer types by using combined gene differential expression analysis methods and a higher-order generalized singular value decomposition. Four co-expression modules were identified; and oncogenes and tumor suppressors were significantly enriched in these modules. Functional enrichment analysis demonstrated the significantly enriched pathways in these modules, including ECM-receptor interaction, focal adhesion and PI3K-Akt signaling pathway. The top-ranked miRNAs (mir-199, mir-29, mir-200) and transcription factors (FOXO4, E2A, NFAT, and MAZ) were identified, which play an important role in deregulating cellular energetics; and regulating angiogenesis and cancer immune system. The clinical significance of the co-expressed gene clusters was assessed by evaluating their predictability of cancer patients' survival. The predictive power of different clusters and subclusters was demonstrated. Our results will be valuable in cancer-related gene function annotation and for the evaluation of cancer patients' prognosis.
| S-EPMC5762563 | biostudies-literature
Project description:Spinocerebellar ataxia type 2 (SCA2) is an autosomal dominant neurodegenerative disease caused by CAG repeat expansion in the ATXN2 gene. The repeat resides in an encoded region of the gene resulting in polyglutamine (polyQ) expansion which has been assumed to result in gain of function, predominantly, for the ATXN2 protein. We evaluated temporal cerebellar expression profiles by RNA sequencing of ATXN2Q127 mice versus wild-type (WT) littermates. ATXN2Q127 mice are characterized by a progressive motor phenotype onset, and have progressive cerebellar molecular and neurophysiological (Purkinje cell firing frequency) phenotypes. Our analysis revealed previously uncharacterized early and progressive abnormal patterning of cerebellar gene expression. Weighted Gene Coexpression Network Analysis revealed four gene modules that were significantly correlated with disease status, composed primarily of genes associated with GTPase signaling, calcium signaling and cell death. Of these genes, few overlapped with differentially expressed cerebellar genes that we identified in Atxn2-/- knockout mice versus WT littermates, suggesting that loss-of-function is not a significant component of disease pathology. We conclude that SCA2 is a disease characterized by gain of function for ATXN2.
| S-EPMC5886232 | biostudies-literature
Project description:Breast cancer (BC) is the most common leading cause of cancer-related death in women worldwide. Gene expression profiling analysis for human BCs has been studied previously. However, co-expression analysis for BC cell lines is still devoid to date. The aim of the study was to identify key pathways and hub genes that may serve as a biomarker for BC and uncover potential molecular mechanism using weighted correlation network analysis. We analyzed microarray data of BC cell lines (GSE 48213) listed in the Gene Expression Omnibus database. Gene co-expression networks were used to construct and explore the biological function in hub modules using the weighted correlation network analysis algorithm method. Meanwhile, Gene ontology and KEGG pathway analysis were performed using Cytoscape plug-in ClueGo. The network of the key module was also constructed using Cytoscape. A total of 5000 genes were selected, 28 modules of co-expressed genes were identified from the gene co-expression network, one of which was found to be significantly associated with a subtype of BC lines. Functional enrichment analysis revealed that the brown module was mainly involved in the pathway of the autophagy, spliceosome, and mitophagy, the black module was mainly enriched in the pathway of colorectal cancer and pancreatic cancer, and genes in midnightblue module played critical roles in ribosome and regulation of lipolysis in adipocytes pathway. Three hub genes CBR3, SF3B6, and RHPN1 may play an important role in the development and malignancy of the disease. The findings of the present study could improve our understanding of the molecular pathogenesis of breast cancer.
| S-EPMC6639467 | biostudies-literature