Identification of potential gene signatures associated with osteosarcoma by integrated bioinformatics analysis.
Ontology highlight
ABSTRACT: Background
Osteosarcoma (OS) is the most primary malignant bone cancer in children and adolescents with a high mortality rate. This work aims to screen novel potential gene signatures associated with OS by integrated microarray analysis of the Gene Expression Omnibus (GEO) database.Material and methods
The OS microarray datasets were searched and downloaded from GEO database to identify differentially expressed genes (DEGs) between OS and normal samples. Afterwards, the functional enrichment analysis, protein-protein interaction (PPI) network analysis and transcription factor (TF)-target gene regulatory network were applied to uncover the biological function of DEGs. Finally, two published OS datasets (GSE39262 and GSE126209) were obtained from GEO database for evaluating the expression level and diagnostic values of key genes.Results
In total 1,059 DEGs (569 up-regulated DEGs and 490 down-regulated DEGs) between OS and normal samples were screened. Functional analysis showed that these DEGs were markedly enriched in 214 GO terms and 54 KEGG pathways such as pathways in cancer. Five genes (CAMP, METTL7A, TCN1, LTF and CXCL12) acted as hub genes in PPI network. Besides, METTL7A, CYP4F3, TCN1, LTF and NETO2 were key genes in TF-gene network. Moreover, Pax-6 regulated four key genes (TCN1, CYP4F3, NETO2 and CXCL12). The expression levels of four genes (METTL7A, TCN1, CXCL12 and NETO2) in GSE39262 set were consistent with our integration analysis. The expression levels of two genes (CXCL12 and NETO2) in GSE126209 set were consistent with our integration analysis. ROC analysis of GSE39262 set revealed that CYP4F3, CXCL12, METTL7A, TCN1 and NETO2 had good diagnostic values for OS patients. ROC analysis of GSE126209 set revealed that CXCL12, METTL7A, TCN1 and NETO2 had good diagnostic values for OS patients.
SUBMITTER: Jia Y
PROVIDER: S-EPMC8164836 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundThis study hoped to explore the potential biomarkers and associated metabolites during osteosarcoma (OS) progression based on bioinformatics integrated analysis.MethodsGene expression profiles of GSE28424, including 19 human OS cell lines (OS group) and 4 human normal long bone tissue samples (control group), were downloaded. The differentially expressed genes (DEGs) in OS vs. control were investigated. The enrichment investigation was performed based on DEGs, followed by protein-protein interaction network analysis. Then, the feature genes associated with OS were explored, followed by survival analysis to reveal prognostic genes. The qRT-PCR assay was performed to test the expression of these genes. Finally, the OS-associated metabolites and disease-metabolic network were further investigated.ResultsTotally, 357 DEGs were revealed between the OS vs. control groups. These DEGs, such as CXCL12, were mainly involved in functions like leukocyte migration. Then, totally, 38 feature genes were explored, of which 8 genes showed significant associations with the survival of patients. High expression of CXCL12, CEBPA, SPARCL1, CAT, TUBA1A, and ALDH1A1 was associated with longer survival time, while high expression of CFLAR and STC2 was associated with poor survival. Finally, a disease-metabolic network was constructed with 25 nodes including two disease-associated metabolites cyclophosphamide and bisphenol A (BPA). BPA showed interactions with multiple prognosis-related genes, such as CXCL12 and STC2.ConclusionWe identified 8 prognosis-related genes in OS. CXCL12 might participate in OS progression via leukocyte migration function. BPA might be an important metabolite interacting with multiple prognosis-related genes.
| S-EPMC8256509 | biostudies-literature
Project description:The purpose of this work was to extract key players such as mRNAs and long non-coding RNA (lncRNAs) in the etiopathogenesis of osteosarcoma (OS). The sequencing analyses (mRNAs and lncRNAs) of OS were conducted followed by differentially expressed mRNAs and lncRNAs (DEmRNAs and DElncRNAs) identification between U-2OS cells with has-miR-590-5p overexpression and negative control cells. Following this, the co-expression and functional enrichment analyses of DEmRNAs and DElncRNAs were carried out. Also, the miRNAs-DElncRNAs-DEmRNAs regulatory network was constructed with DElncRNAs-miRNAs and DElncRNAs-DEmRNAs pairs after the target gene analysis of miRNA. In addition, the ceRNA-has-miR-590-5p was further extracted based on the has-miR-590-5p-DElncRNAs and DElncRNAs-DEmRNAs interactions. Finally, the results of the bioinformatics analysis was verified by reverse-transcription polymerase chain reaction (RT-PCR). Totally, 980 DEmRNAs (539 up-regulated DEmRNAs and 441 down-regulated DEmRNAs) and 682 DElncRNAs (352 up-regulated DElncRNAs and 330 down-regulated DElncRNAs) were extracted between cells with hsa-miR-590-5p overexpression and normal cells. The functional analyses suggested that up-regulated genes were significantly enriched in several GO terms such as signal transduction and cytokine-cytokine receptor interaction pathway while down-regulated genes (SCUBE3, HIST1H4E and EDIL3) were associated with calcium ion binding, cell surface function and nucleosome assembly. Additionally, the miRNAs-DEmRNAs-DEmRNAs network represented 220 pairs among 41 miRNAs, 38 DElncRNAs and 61 DEmRNAs. Furthermore, the ceRNA-hsa-miR-590-5p network consisted of 70 interaction pairs including hsa-miR-590-5p-SCUBE3-CTB-113D17.1, hsa-miR-590-5p-EDIL3-CTB-113D17.1 and hsa-miR-590-5p-HIST1H4E-CTB-113D17.1) among hsa-miR-590-5p, 30 DEmRNAs and 4 down-regulated DElncRNAs. Meanwhile, the RT-PCR results incidated that compared with the blank (KB) and negative control (NC) group, the mRNA expression of SCUBE3, HIST1H4E, and EDIL3 were significantly descreased in mimics group (P value <0.05). The lncRNA CTB-113D17.1 might implicate with OS development probably via serving as a hsa-miR-590-5p sponge to regulate gene targets (SCUBE3, EDIL3 and HIST1H4E), which will facilitate the deep understandings of OS progression.
| S-EPMC7488886 | biostudies-literature
Project description:BackgroundIgA nephropathy (IgAN) is the most common form of primary glomerulonephritis worldwide, and its diagnosis depends mainly on renal biopsy. However, there is no specific treatment for IgAN. Moreover, its causes and underlying molecular events require further exploration.MethodsThe expression profiles of GSE64306 and GSE93798 were downloaded from the Gene Expression Omnibus (GEO) database and used to identify the differential expression of miRNAs and genes, respectively. The StarBase and TransmiR databases were employed to predict target genes and transcription factors of the differentially expressed miRNAs (DE-miRNAs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were conducted to predict biological functions. A comprehensive analysis of the miRNA-mRNA regulatory network was constructed, and protein-protein interaction (PPI) networks and hub genes were identified. CIBERSORT was used to examine the immune cells in IgAN, and correlation analyses were performed between the hub genes and infiltrating immune cells.ResultsFour downregulated miRNAs and 16 upregulated miRNAs were identified. Forty-five and twelve target genes were identified for the upregulated and downregulated DE-miRNAs, respectively. CDKN1A, CDC23, EGR1, HIF1A, and TRIM28 were the hub genes with the highest degrees of connectivity. CIBERSORT revealed increases in the numbers of activated NK cells, M1 and M2 macrophages, CD4 naive T cells, and regulatory T cells in IgAN. Additionally, HIF1A, CDC23, TRIM28, and CDKN1A in IgAN patients were associated with immune cell infiltration.ConclusionsA potential miRNA-mRNA regulatory network contributing to IgAN onset and progression was successfully established. The results of the present study may facilitate the diagnosis and treatment of IgAN by targeting established miRNA-mRNA interaction networks. Infiltrating immune cells may play significant roles in IgAN pathogenesis. Future studies on these immune cells may help guide immunotherapy for IgAN patients.
| S-EPMC8620631 | biostudies-literature
Project description:BackgroundDiabetic kidney disease (DKD) is a common complication of diabetes that is clinically characterized by progressive albuminuria due to glomerular destruction. The etiology of DKD is multifactorial, and numerous studies have demonstrated that cellular senescence plays a significant role in its pathogenesis, but the specific mechanism requires further investigation.MethodsThis study utilized 5 datasets comprising 144 renal samples from the Gene Expression Omnibus (GEO) database. We obtained cellular senescence-related pathways from the Molecular Signatures Database and evaluated the activity of senescence pathways in DKD patients using the Gene Set Enrichment Analysis (GSEA) algorithm. Furthermore, we identified module genes related to cellular senescence pathways through Weighted Gene Co-Expression Network Analysis (WGCNA) algorithm and used machine learning algorithms to screen for hub genes related to senescence. Subsequently, we constructed a cellular senescence-related signature (SRS) risk score based on hub genes using the Least Absolute Shrinkage and Selection Operator (LASSO), and verified mRNA levels of hub genes by RT-PCR in vivo. Finally, we validated the relationship between the SRS risk score and kidney function, as well as their association with mitochondrial function and immune infiltration.ResultsThe activity of cellular senescence-related pathways was found to be elevated among DKD patients. Based on 5 hub genes (LIMA1, ZFP36, FOS, IGFBP6, CKB), a cellular senescence-related signature (SRS) was constructed and validated as a risk factor for renal function decline in DKD patients. Notably, patients with high SRS risk scores exhibited extensive inhibition of mitochondrial pathways and upregulation of immune cell infiltration.ConclusionCollectively, our findings demonstrated that cellular senescence is involved in the process of DKD, providing a novel strategy for treating DKD.
| S-EPMC10313062 | biostudies-literature
Project description:BackgroundHepatocellular carcinoma (HCC) is a type of primary liver tumor with poor prognosis and high mortality, and its molecular mechanism remains incompletely understood. This study aimed to use bioinformatics technology to identify differentially expressed genes (DEGs) in HCC pathogenesis, hoping to identify novel biomarkers or potential therapeutic targets for HCC research.MethodsThe bioinformatics analysis of our research mostly involved the following two datasets: Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA). First, we screened DEGs based on the R packages (limma and edgeR). Using the DAVID database, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of DEGs were carried out. Next, the protein-protein interaction (PPI) network of the DEGs was built in the STRING database. Then, hub genes were screened through the cytoHubba plug-in, followed by verification using the GEPIA and Oncomine databases. We demonstrated differences in levels of the protein in hub genes using the Human Protein Atlas (HPA) database. Finally, the hub genes prognostic values were analyzed by the GEPIA database. Additionally, using the Comparative Toxicogenomics Database (CTD), we constructed the drug-gene interaction network.ResultsWe ended up with 763 DEGs, including 247 upregulated and 516 downregulated DEGs, that were mainly enriched in the epoxygenase P450 pathway, oxidation-reduction process, and metabolism-related pathways. Through the constructed PPI network, it can be concluded that the P53 signaling pathway and the cell cycle are the most obvious in module analysis. From the PPI, we filtered out eight hub genes, and these genes were significantly upregulated in HCC samples, findings consistent with the expression validation results. Additionally, survival analysis showed that high level gene expression of CDC20, CDK1, MAD2L1, BUB1, BUB1B, CCNB1, and CCNA2 were connected with the poor overall survival of HCC patients. Toxicogenomics analysis showed that only topotecan, oxaliplatin, and azathioprine could reduce the gene expression levels of all seven hub genes.ConclusionThe present study screened out the key genes and pathways that were related to HCC pathogenesis, which could provide new insight for the future molecularly targeted therapy and prognosis evaluation of HCC.
| S-EPMC7985067 | biostudies-literature
Project description:BackgroundGastric cancer (GC) is one of the common gastrointestinal malignancy worldwide and exhibits a poor prognosis. Increasing studies have indicated that microRNAs play critical roles in the cancer progression and have shown great potential as useful biomarkers. The search for potential diagnostic and prognostic biomarkers of gastric cancer (GC) with integrated bioinformatics analyses has been undertaken in previous studies.MethodsIn this study, the robust rank aggregation (RRA) method was used to perform an integrated analysis of differentially expressed miRNAs (DEMs) from five microarray datasets in the Gene Expression Omnibus (GEO) database to find robust biomarkers for GC. Ultimately, seven miRNAs were filtered from fourteen primary miRNAs using the validation set of The Cancer Genome Atlas (TCGA) database. Based on these results, diagnostic and survival analyses were performed, and logistic regression and Cox regression were used to determine the clinicopathological characteristics of the DEM expression and overall survival.ResultsNine eligible miRNA datasets related to GC were selected from the GEO database for integrated analysis in this study. Diagnostic analysis implied that these miRNAs could be regarded as promising candidate diagnostic biomarkers in GC tissues, but whether the results of the tissue analysis are consistent with those of peripheral blood analysis requires further validation. The logistic regression indicated that the ectopic expression of these DEMs was relevant to the histological type, anatomical region, and pathological grade of GC. However, the survival and Cox regression analyses suggested that the poor prognosis of GC patients was not strongly dependent on the ectopic expression of the seven miRNAs, but rather, a poor prognosis was associated with age, metastasis, and histological grade.ConclusionsBased on the results presented in this study it can be concluded that these miRNAs (miR-455-3p, miR-135b-5p, let-7a-3p, miR-195-5p, miR-204-5p, miR-149-5p, and miR-143-3p) might be potential biomarkers for the early diagnosis of GC patients, but this finding should be regarded with caution. A large-scale, prospective, and multicenter cohort study should be performed.
| S-EPMC8246217 | biostudies-literature
Project description:Selecting differentially expressed genes (DEGs) based on integrated bioinformatics analyses has been used in previous studies to explore potential biomarkers in gastric cancer (GC) with microarray and RNA sequencing data. However, the genes obtained may be inaccurate because of noisy data and errors, as well as insufficient clinical sample sizes. Thus, we aimed to find robust and strong DEGs with prognostic value for GC, where the robust rank aggregation method was employed to select significant DEGs from eight Gene Expression Omnibus data sets with a total of 140 up-regulated and 206 down-regulated genes. Network data mining was then used to screen hub genes, and 11 genes were filtered using Fisher's exact test. Based on these results, we built a prognostic signature with seven genes (FBN1, MMP1, PLAU, SPARC, COL1A2, COL2A1 and ATP4A) using stepwise multivariate Cox proportional hazard regression. According to the risk score for each patient, we found that high-risk group patients had significantly worse survival results compared with those in the low-risk group (log-rank test P-value < 0.001). This seven-gene signature was then validated with an external data set. Thus, we established a signature based on seven DEGs with prognostic value for GC patients using multi-steps bioinformatics methods, which may provide novel insights and potential biomarkers for prognosis, as well as possibly serving as new therapeutic targets in clinical applications.
| S-EPMC6201382 | biostudies-literature
Project description:Hepatocellular carcinoma (HCC) is a common yet deadly form of malignant cancer. However, the specific mechanisms involved in HCC diagnosis have not yet fully elucidated. Herein, we screened four publically available Gene Expression Omnibus (GEO) expression profiles (GSE14520, GSE29721, GSE45267 and GSE60502), and used them to identify 409 differentially expressed genes (DEGs), including 142 and 267 up- and down-regulated genes, respectively. The DAVID database was used to look for functionally enriched pathways among DEGs, and the STRING database and Cytoscape platform were used to generate a protein-protein interaction (PPI) network for these DEGs. The cytoHubba plug-in was utilized to detect 185 hub genes, and three key clustering modules were constructed with the MCODE plug-in. Gene functional enrichment analyses of these three key clustering modules were further performed, and nine core genes including BIRC5, DLGAP5, DTL, FEN1, KIAA0101, KIF4A, MCM2, MKI67, and RFC4, were identified in the most critical cluster. Subsequently, the hierarchical clustering and expression of core genes in TCGA liver cancer tissues were analyzed using the UCSC Cancer Genomics Browser, and whether elevated core gene expression was linked to a poor prognosis in HCC patients was assessed using the GEPIA database. The PPI of the nine core genes revealed an interaction between FEN1, MCM2, RFC4, and BIRC5. Furthermore, the expression of FEN1 was positively correlated with that of three other core genes in TCGA liver cancer tissues. FEN1 expression in HCC and other tumor types was assessed with the FIREBROWSE and ONCOMINE databases, and results were verified in HCC samples and hepatoma cells. FEN1 levels were also positively correlated with tumor size, distant metastasis and vascular invasion. In conclusion, we identified nine core genes associated with HCC development, offering novel insight into HCC progression. In particular, the aberrantly elevated FEN1 may represent a potential biomarker for HCC diagnosis and treatment.
| S-EPMC6733258 | biostudies-literature
Project description:BackgroundApoptosis is associated with the pathogenesis of Mycobacterium tuberculosis infection. This study aims to identify apoptosis-related genes as biomarkers for differentiating active tuberculosis (ATB) from latent tuberculosis infection (LTBI).MethodsThe tuberculosis (TB) datasets (GSE19491, GSE62525, and GSE28623) were downloaded from the Gene Expression Omnibus (GEO) database. The diagnostic biomarkers differentiating ATB from LTBI were identified by combining the data of protein-protein interaction network, differentially expressed gene, Weighted Gene Co-Expression Network Analysis (WGCNA), and receiver operating characteristic (ROC) analyses. Machine learning algorithms were employed to validate the diagnostic ability of the biomarkers. Enrichment analysis for biomarkers was established, and potential drugs were predicted. The association between biomarkers and N6-methyladenosine (m6A) or 5-methylated cytosine (m5C) regulators was evaluated.ResultsSix biomarkers including CASP1, TNFSF10, CASP4, CASP5, IFI16, and ATF3 were obtained for differentiating ATB from LTBI. They showed strong diagnostic performances, with area under ROC (AUC) values > 0.7. Enrichment analysis demonstrated that the biomarkers were involved in immune and inflammation responses. Furthermore, 24 drugs, including progesterone and emricasan, were predicted. The correlation analysis revealed that biomarkers were positively correlated with most m6A or m5C regulators.ConclusionThe six ARGs can serve as effective biomarkers differentiating ATB from LTBI and provide insight into the pathogenesis of Mycobacterium tuberculosis infection.
| S-EPMC10834671 | biostudies-literature
Project description:PurposeSmoking is recognized as a risk factor for Chronic Obstructive Pulmonary Disease (COPD), yet only 20-25% of smokers eventually develop COPD. Since its molecular pathogenesis remains unclear, there is an important need to further understand genetic differences between smokers with COPD and healthy smokers, screen out high-risk and susceptible groups among smokers, and find effective therapeutic targets.MethodsBioinformatics tools were used to screen biomarkers that were significantly associated with COPD smokers and healthy smokers. qRT-PCR and Western blotting analysis were used to detect hub gene expression in CSE-treated BEAS-2B cells and lung tissue of COPD mouse models.ResultsOur study identified 132 DEGs. The GO and KEGG analyses suggested that the ECM-receptor interaction, MAPK signaling pathway, Chemokine signaling pathway, PI3K-Akt signaling pathway, extracellular matrix organization and collagen fibril organization were associated with the occurrence and development of COPD. In addition, WGCNA analysis of GSE1650 showed that the brown module was most correlated with COPD. The intersection between the brown module and DEGs was used to identify 9 HUB genes (COL14A1, SULF1, MOXD1, CXCL12, CHRNA1, COMP, POU2AF1, MMP11, THBS2) that showed consistent expression and upregulation. Both the mRNA and protein expression levels of the Hub genes (except that of MMP11) were significantly upregulated in tobacco smoke exposed mouse emphysema models and CSE treated BEAS-2B cells.ConclusionOur results suggest that COL14A1, SULF1, MOXD1, CXCL12, CHRNA1, COMP, POU2AF1, and THBS2 may be potentially useful biomarkers for identifying smokers with a risk of developing COPD. The GO and KEGG functional enrichment analyses further confirmed the significant role played by ECM in the pathogenesis of COPD. The results of this study may provide further insights into the pathogenetic mechanisms involved in COPD.
| S-EPMC8901430 | biostudies-literature