Dataset Information

Identification and analysis of key genes in osteosarcoma using bioinformatics.

ABSTRACT: Osteosarcoma (OS) is an invasive malignant neoplasm of the bones. The present study identified and analyzed key genes associated with OS. Expression profiling of the dataset GSE49003, which included 6 metastatic and 6 non-metastatic OS cell lines and was obtained from the Gene Expression Omnibus, was performed. Following data preprocessing, the differentially expressed genes (DEGs) were selected using the limma package in R. Subsequently, bidirectional hierarchical clustering using the pheatmap package in R and an unpaired Students' t-test was performed for the DEGs. Based on the Search Tool for the Retrieval of Interacting Genes database and Cytoscape software, a protein-protein interaction (PPI) network for the DEGs was constructed. Using Database for Annotation, Visualization and Integrated Discovery software and the Kyoto Encyclopedia of Genes and Genomes Orthology Based Annotation System server, functional and pathway enrichment analyses were performed for the DEGs corresponding to the proteins of the network. In addition, the transcription factors (TFs) and CpG islands of the gene promoter were searched for using the TRANSFAC database and CpG Island Searcher software, respectively. A total of 323 DEGs were identified between the metastatic and non-metastatic samples. In the PPI network, upregulated epidermal growth factor receptor (EGFR) exhibits a high degree and was therefore highly interconnected with other proteins. Enrichment analysis revealed that EGFR was enriched in cytoskeleton organization, organic substance response and the signaling pathway of focal adhesion. The TFs early growth response 1, nuclear factor-?B complex subunits, peroxisome proliferator activated receptor ?, signal transducer and activator of transcription 3 and MYC proto-oncogene were identified in the EGFR promoter region. Furthermore, multiple CpG islands, starting from the 400 bp of the EGFR promoter sequence, were predicted. Methylated modification of the CpG islands in the EGFR promoter may help to regulate EGFR expression. The TFs identified in the EGFR promoter may function in the progression of OS.

SUBMITTER: Diao C

PROVIDER: S-EPMC5778824 | biostudies-literature | 2018 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Identification and analysis of key genes in osteosarcoma using bioinformatics.

Diao Chunyu C Xi Yong Y Xiao Tao T

Oncology letters 20171219 3

Osteosarcoma (OS) is an invasive malignant neoplasm of the bones. The present study identified and analyzed key genes associated with OS. Expression profiling of the dataset GSE49003, which included 6 metastatic and 6 non-metastatic OS cell lines and was obtained from the Gene Expression Omnibus, was performed. Following data preprocessing, the differentially expressed genes (DEGs) were selected using the limma package in R. Subsequently, bidirectional hierarchical clustering using the pheatmap ...[more]

PMID: 29435005

Similar Datasets

Project description:BackgroundOsteosarcoma (OS) metastasis is the most common cause of cancer-related mortality, however, no sufficient clinical biomarkers have been identified. In this study, we identified five genes to help predict metastasis at diagnosis.MethodsWe performed weighted gene co-expression network analysis (WGCNA) to identify the most relevant gene modules associated with OS metastasis. An important machine learning algorithm, the support vector machine (SVM), was employed to predict key genes for classifying the OS metastasis phenotype. Finally, we investigated the clinical significance of key genes and their enriched pathways.ResultsEighteen modules were identified in WGCNA, among which the pink, red, brown, blue, and turquoise modules demonstrated good preservation. In the five modules, the brown and red modules were highly correlated with OS metastasis. Genes in the two modules closely interacted in protein-protein interaction networks and were therefore chosen for further analysis. Genes in the two modules were primarily enriched in the biological processes associated with tumorigenesis and development. Furthermore, 65 differentially expressed genes were identified as common hub genes in both WGCNA and protein-protein interaction networks. SVM classifiers with the maximum area under the curve were based on 30 and 15 genes in the brown and red modules, respectively. The clinical significance of the 45 hub genes was analyzed. Of the 45 genes, 17 were found to be significantly correlated with survival time. Finally, 5/17 genes, including ADAP2 (P = 0.0094), LCP2 (P = 0.013), ARHGAP25 (P = 0.0049), CD53 (P = 0.016), and TLR7 (P = 0.04) were significantly correlated with the metastatic phenotype. In vitro verification, western blotting, wound healing analyses, transwell invasion assays, proliferation assays, and colony formation assays indicated that ARHGAP25 promoted OS cell migration, invasion, proliferation, and epithelial-mesenchymal transition.ConclusionWe identified five genes, namely ADAP2, LCP2, ARHGAP25, CD53, and TLR7, as candidate biomarkers for the prediction of OS metastasis; ARHGAP25 inhibits MG63 OS cell growth, migration, and invasion in vitro, indicating that ARHGAP25 can serve as a promising specific and prognostic biomarker for OS metastasis.

Project description:Osteosarcoma, which arises from bone tissue, is considered to be one of the most common types of cancer in children and teenagers. As the etiology of osteosarcoma has not been fully elucidated, the overall prognosis for patients is generally poor. In recent years, the development of bioinformatical technology has allowed researchers to identify numerous molecular biological characteristics associated with the prognosis of osteosarcoma using online databases. In the present study, Gene Expression Omnibus (GEO) database was used and three microarray datasets were obtained. The GEO2R web tool was utilized and differentially expressed genes (DEGs) in osteosarcoma tissue were identified. Venn analysis was performed to determine the intersection of the DEG profiles. DEGs were analyzed by Gene Ontology function and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis. Protein-protein interactions (PPIs) between these DEGs were analyzed using the Search Tool for the Retrieval of Interacting Genes database, and the PPI network was then visualized using Cytoscape software. The top ten genes were identified based on measurement of degree, density of maximum neighborhood component, maximal clique centrality and mononuclear cell counts in the PPI network, and five overlapping genes [origin recognition complex subunit 6 (ORC6), IGF-binding protein 5 (IGFBP5), minichromosome maintenance 10 replication initiation factor (MCM10), MET proto-oncogene, receptor tyrosine kinase (MET) and centromere protein F (CENPF)] were identified. Additionally, three module networks were analyzed by Molecular Complex Detection (MCODE), and six key genes [ORC6, MCM10, DEP domain containing 1 (DEPDC1), CENPF, TIMELESS interacting protein (TIPIN) and shugoshin 1 (SGOL1)] were screened. Combined with the results from Cytoscape and MCODE, eight hub genes (ORC6, MCM10, DEPDC1, CENPF, TIPIN, SGOL1, MET and IGFBP5) were obtained. Furthermore, Kaplan-Meier plotter survival analysis was used to evaluate the prognostic value of these eight hub genes in patients with osteosarcoma. Oncomine and GEPIA databases were applied to further confirm the expression levels of hub genes in tissue. Finally, the functional roles of the core gene CENPF were investigated using Cell Counting Kit-8, wound healing and Transwell assays, which indicated that CENPF knockdown inhibited the proliferation, migration and invasion of osteosarcoma cells. These results provided potential prognostic markers, as well as a basis for further investigation of the mechanism underlying osteosarcoma.

Project description:Gastric cancer (GC) is one of the most common types of cancer worldwide. Patients must be identified at an early stage of tumor progression for treatment to be effective. The aim of the present study was to identify potential biomarkers with diagnostic value in patients with GC. To examine potential therapeutic targets for GC, four Gene Expression Omnibus (GEO) datasets were downloaded and screened for differentially expressed genes (DEGs). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were subsequently performed to study the function and pathway enrichment of the identified DEGs. A protein-protein interaction (PPI) network was constructed. The CytoHubba plugin of Cytoscape was used to calculate the degree of connectivity of proteins in the PPI network, and the two genes with the highest degree of connectivity were selected for further analysis. Additionally, the two DEGs with the largest and smallest log Fold Change values were selected. These six key genes were further examined using Oncomine and the Kaplan-Meier plotter platform. A total of 99 upregulated and 172 downregulated genes common to all four GEO datasets were screened. The DEGs were primarily enriched in the Biological Process terms: 'extracellular matrix organization', 'collagen catabolic process' and 'cell adhesion'. These three KEGG pathways were significantly enriched in the categories: 'ECM-receptor interaction', 'protein digestion and absorption', and 'focal adhesion'. Based on Oncomine, expression of ATP4A and ATP4B were downregulated in GC, whereas expression of the other genes were all upregulated. The Kaplan-Meier plotter platform confirmed that upregulated expression of the identified key genes was significantly associated with worse overall survival of patients with GC. The results of the present study suggest that FN1, COL1A1, INHBA and CST1 may be potential biomarkers and therapeutic targets for GC. Additional studies are required to explore the potential value of ATP4A and ATP4B in the treatment of GC.

Project description:Nasopharyngeal carcinoma (NPC) is one of the most common malignancies in the head and neck. The aim of the current study was to identify the key pathways and genes involved in NPC through bioinformatics analysis and to identify potential molecular mechanisms underlying NPC proliferation and progression. Three gene expression profiles (GSE12452, GSE34573 and GSE64634) were downloaded from the Gene Expression Omnibus database. A total of 76 samples were analyzed, of which 59 were NPC samples and 17 were normal samples. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were subsequently conducted. The protein-protein interaction (PPI) network of the differentially expressed genes (DEGs) was constructed using Cytoscape software. Analysis of GSE12452, GSE34573 and GSE64634 datasets identified 1,301 (553 upregulated and 748 downregulated), 1,232 (348 upregulated and 884 downregulated) and 1,218 (555 upregulated and 663 downregulated) DEGs, respectively. Using Venn diagram analysis, 268 DEGs (59 upregulated and 209 downregulated) that intersected all three datasets, were selected for further analysis. The results of GO analysis revealed that upregulated DEGs were significantly enriched in biological processes, including 'cell adhesion', 'cell division', 'mitosis' and 'mitotic cell cycle'. The downregulated DEGs were mainly enriched in 'microtubule-based movement', 'cilium movement', 'cilium axoneme assembly' and 'epithelial cell differentiation'. The KEGG pathway analysis results revealed that the upregulated DEGs were highly associated with several pathways, including 'extracellular matrix-receptor interaction', 'human papillomavirus infection', 'arrhythmogenic right ventricular cardiomyopathy' and 'focal adhesion', whereas the downregulated DEGs were enriched in 'metabolic pathways', 'Huntington's disease', 'fluid shear stress and atherosclerosis' and 'chemical carcinogenesis'. On the basis of the PPI network of the DEGs, the following top 10 hub genes were identified: Dynein axonemal light intermediate chain 1, dynein axonemal intermediate chain 2, calmodulin 1, coiled-coil domain containing 114, dynein axonemal heavy chain 5, radial spoke head 9 homolog, radial spoke head component 4A, NDC80 kinetochore complex component, thymidylate synthetase and coiled-coil domain containing 39. In conclusion, by performing a comprehensive bioinformatics analysis of DEGs, putative targets that could be used to elucidate the molecular mechanisms underlying NPC were identified.

Dataset Information

Identification and analysis of key genes in osteosarcoma using bioinformatics.

Publications

Identification and analysis of key genes in osteosarcoma using bioinformatics.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets