Dataset Information

Identification of potential biomarkers related to glioma survival by gene expression profile analysis.

ABSTRACT:

Background

Recent studies have proposed several gene signatures as biomarkers for different grades of gliomas from various perspectives. However, most of these genes can only be used appropriately for patients with specific grades of gliomas.

Methods

In this study, we aimed to identify survival-relevant genes shared between glioblastoma multiforme (GBM) and lower-grade glioma (LGG), which could be used as potential biomarkers to classify patients into different risk groups. Cox proportional hazard regression model (Cox model) was used to extract relative genes, and effectiveness of genes was estimated against random forest regression. Finally, risk models were constructed with logistic regression.

Results

We identified 104 key genes that were shared between GBM and LGG, which could be significantly correlated with patients' survival based on next-generation sequencing data obtained from The Cancer Genome Atlas for gene expression analysis. The effectiveness of these genes in the survival prediction of GBM and LGG was evaluated, and the average receiver operating characteristic curve (ROC) area under the curve values ranged from 0.7 to 0.8. Gene set enrichment analysis revealed that these genes were involved in eight significant pathways and 23 molecular functions. Moreover, the expressions of ten (CTSZ, EFEMP2, ITGA5, KDELR2, MDK, MICALL2, MAP 2?K3, PLAUR, SERPINE1, and SOCS3) of these genes were significantly higher in GBM than in LGG, and comparing their expression levels to those of the proposed control genes (TBP, IPO8, and SDHA) could have the potential capability to classify patients into high- and low- risk groups, which differ significantly in the overall survival. Signatures of candidate genes were validated, by multiple microarray datasets from Gene Expression Omnibus, to increase the robustness of using these potential prognostic factors. In both the GBM and LGG cohort study, most of the patients in the high-risk group had the IDH1 wild-type gene, and those in the low-risk group had IDH1 mutations. Moreover, most of the high-risk patients with LGG possessed a 1p/19q-noncodeletion.

Conclusion

In this study, we identified survival relevant genes which were shared between GBM and LGG, and those enabled to classify patients into high- and low-risk groups based on expression level analysis. Both the risk groups could be correlated with the well-known genetic variants, thus suggesting their potential prognostic value in clinical application.

SUBMITTER: Hsu JB

PROVIDER: S-EPMC7402580 | biostudies-literature | 2019 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Identification of potential biomarkers related to glioma survival by gene expression profile analysis.

Hsu Justin Bo-Kai JB Chang Tzu-Hao TH Lee Gilbert Aaron GA Lee Tzong-Yi TY Chen Cheng-Yu CY

BMC medical genomics 20190320 Suppl 7

<h4>Background</h4>Recent studies have proposed several gene signatures as biomarkers for different grades of gliomas from various perspectives. However, most of these genes can only be used appropriately for patients with specific grades of gliomas.<h4>Methods</h4>In this study, we aimed to identify survival-relevant genes shared between glioblastoma multiforme (GBM) and lower-grade glioma (LGG), which could be used as potential biomarkers to classify patients into different risk groups. Cox pr ...[more]

PMID: 30894197

Similar Datasets

Project description:PurposeIncreasing evidence suggests myopia is not a simple refractive error, many other factors might also be involved. Here, we assessed myopic and normal corneas' gene expression profiles to identify possible diagnostic and therapeutic biomarkers for myopia.Materials and methodsWe obtained the expression profile of ten patients and seven normal control samples from the GSE112155 and GSE151631 datasets based on the Gene Expression Omnibus (GEO) database. We used the "limma" R package to determine the differentially expressed genes (DEGs) between myopic and normal corneas. Weighted gene co-expression network analysis (WGCNA) was used to identify critical co-expressed modules related to myopia, and enrichment analyses were used to annotate the function of genes encompassed in the compulsory module. We also validated these findings in two external datasets (GSE24641 and GSE136701).ResultsWe identified that the DEGs were significantly enriched in ultraviolet (UV) response, TNF-α signaling via NFκB, Angiogenesis, Myogenesis pathways, etc. We used 2095 genes to construct the co-expression gene modules and found five interesting modules because the eigengene expression of these modules was significantly differentially expressed between myopic and normal corneas. Notably, the enrichment analysis found that the genes encompassed in lightgreen module were significantly enriched in immune-related pathways. These findings were proved by subsequent analysis based on Xcell software. We found the component of B cells, CD4+ memory T cells, CD8+ central memory T cells, plasmacytoid dendritic cells, T helper 2 (Th2) cells, regulatory T cells (Tregs), etc. were significantly increased in myopic corneas, while CD8+ T cells, CD4+ T central memory cells, natural killer T (NKT) cells, and T helper 1 (Th1) cells were significantly decreased.ConclusionOur findings identified some markers that might detect diagnosis and treatment for myopia from cornea aspect. Future studies are warranted to verify the functional role of immune-related pathways in cornea during the pathogenesis or progression of myopia.

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.

Project description:Background: Burn injury is a life-threatening disease that does not have ideal biomarkers. Therefore, this study first applied weighted gene co-expression network analysis (WGCNA) and differentially expressed gene (DEG) screening methods to identify pivotal genes and diagnostic biomarkers associated with the skin burn process. Methods: After obtaining transcriptomic datasets of burn patient skin and normal skin from Gene Expression Omnibus (GEO) and performing differential analysis and functional enrichment, WGCNA was used to identify hub gene modules associated with burn skin processes in the burn patient peripheral blood sample dataset and determine the correlation between modules and clinical features. Enrichment analysis was performed to identify the functions and pathways of key module genes. Differential analysis, WGCNA, protein-protein interaction analysis, and enrichment analysis were utilized to screen for hub genes. Hub genes were validated in two other GEO datasets, tested by immunohistochemistry for hub gene expression in burn patients, and receiver operating characteristic curve analysis was performed. Finally, we constructed the specific drug activity, transcription factors, and microRNA regulatory network of the five hub genes. Results: A total of 1,373 DEGs in GSE8056 were obtained, and the top 5 upregulated genes were S100A12, CXCL8, CXCL5, MMP3, and MMP1, whereas the top 5 downregulated genes were SCGB1D2, SCGB2A2, DCD, TSPAN8, and KRT25. DEGs were significantly enriched in the immunity, epidermal development, and skin development processes. In WGCNA, the yellow module was identified as the most closely associated module with tissue damage during the burn process, and the five hub genes (ANXA3, MCEMP1, MMP9, S100A12, and TCN1) were identified as the key genes for burn injury status, which consistently showed high expression in burn patient blood samples in the GSE37069 and GSE13902 datasets. Furthermore, we verified using immunohistochemistry that these five novel hub genes were also significantly elevated in burn patient skin. In addition, MCEMP1, MMP9, and S100A12 showed perfect diagnostic performance in the receiver operating characteristic analysis. Conclusion: In conclusion, we analyzed the changes in genetic processes in the skin during burns and used them to identify five potential novel diagnostic markers in blood samples from burn patients, which are important for burn patient diagnosis. In particular, MCEMP1, MMP9, and S100A12 are three key blood biomarkers that can be used to identify skin damage in burn patients.

Project description:BackgroundEndometriosis is an inflammatory gynecological disease leading to deep pelvic pain, dyspareunia, and infertility. The pathophysiology of endometriosis is complex and depends on a variety of biological processes and pathways. Therefore, there is an urgent need to identify reliable biomarkers for early detection and accurate diagnosis to predict clinical outcomes and aid in the early intervention of endometriosis. We screened transcription factor- (TF-) immune-related gene (IRG) regulatory networks as potential biomarkers to reveal new molecular subgroups for the early diagnosis of endometriosis.MethodsTo explore potential therapeutic targets for endometriosis, the Gene Expression Omnibus (GEO), Immunology Database and Analysis Portal (ImmPort), and TF databases were used to obtain data related to the recognition of differentially expressed genes (DEGs), differentially expressed IRGs (DEIRGs), and differentially expressed TFs (DETFs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on the DETFs and DEIRGs. Then, DETFs and DEIRGs were further validated in the external datasets of GSE51981 and GSE1230103. Then, we used quantitative real-time polymerase chain reaction (qRT-PCR) to verify the hub genes. Simultaneously, the Pearson correlation analysis and protein-protein interaction (PPI) analyses were used to indicate the potential mechanisms of TF-IRGs at the molecular level and obtain hub IRGs. Finally, the receiver operating characteristic (ROC) curve analysis was used to assess the diagnostic value of the hub IRGs.ResultsWe screened a total of 94 DETFs and 121 DEIRGs in endometriosis. Most downregulated DETFs showed decreased expression in the endometria of moderate/severe endometriosis patients. The top-ranked upregulated DEIRGs were upregulated in the endometra of infertile women. Functional analysis showed that DETFs and DEIRGs may be involved in the biological behaviors and pathways of endometriosis. The TF-IRG PPI network was successfully constructed. Compared with the control group, high C3, VCAM1, ITGB2, and C3AR1 expression had statistical significance in endometriosis among the hub DEIRGs. They also showed higher sensitivity and specificity by ROC analysis for the diagnosis of endometriosis. Finally, compared with controls, C3 and VCAM1 were highly expressed in endometriosis tissue samples. In addition, they also showed high specificity and sensitivity for diagnosing endometriosis.ConclusionOverall, we discovered the TF-IRG regulatory network and analyzed 4 hub IRGs that were closely related to endometriosis, which contributes to the diagnosis of endometriosis. Additionally, we verified that DETFs or DEIRGs were associated with the clinicopathological features of endometriosis, and external datasets also confirmed the hub IRGs. Finally, C3 and VCAM1 were highly expressed in endometriosis tissue samples compared with controls and may be potential biomarkers of endometriosis, which are helpful for the early diagnosis of endometriosis.

Project description:PurposeGlaucoma is the second commonest cause of blindness. We assessed the gene expression profile of astrocytes in the optic nerve head to identify possible prognostic biomarkers for glaucoma.MethodA total of 20 patient and nine normal control subject samples were derived from the GSE9944 (six normal samples and 13 patient samples) and GSE2378 (three normal samples and seven patient samples) datasets, screened by microarray-tested optic nerve head tissues, were obtained from the Gene Expression Omnibus (GEO) database. We used a weighted gene coexpression network analysis (WGCNA) to identify coexpressed gene modules. We also performed a functional enrichment analysis and least absolute shrinkage and selection operator (LASSO) regression analysis. Genes expression was represented by boxplots, functional geneset enrichment analyses (GSEA) were used to profile the expression patterns of all the key genes. Then the key genes were validated by the external dataset.ResultsA total 8,606 genes and 19 human optic nerve head samples taken from glaucoma patients in the GSE9944 were compared with normal control samples to construct the co-expression gene modules. After selecting the most common clinical traits of glaucoma, their association with gene expression was established, which sorted two modules showing greatest correlations. One with the correlation coefficient is 0.56 (P = 0.01) and the other with the correlation coefficient is -0.56 (P = 0.01). Hub genes of these modules were identified using scatterplots of gene significance versus module membership. A functional enrichment analysis showed that the former module was mainly enriched in genes involved in cellular inflammation and injury, whereas the latter was mainly enriched in genes involved in tissue homeostasis and physiological processes. This suggests that genes in the green-yellow module may play critical roles in the onset and development of glaucoma. A LASSO regression analysis identified three hub genes: Recombinant Bone Morphogenetic Protein 1 gene (BMP1), Duchenne muscular dystrophy gene (DMD) and mitogens induced GTP-binding protein gene (GEM). The expression levels of the three genes in the glaucoma group were significantly lower than those in the normal group. GSEA further illuminated that BMP1, DMD and GEM participated in the occurrence and development of some important metabolic progresses. Using the GSE2378 dataset, we confirmed the high validity of the model, with an area under the receiver operator characteristic curve of 85%.ConclusionWe identified several key genes, including BMP1, DMD and GEM, that may be involved in the pathogenesis of glaucoma. Our results may help to determine the prognosis of glaucoma and/or to design gene- or molecule-targeted drugs.

Dataset Information

Identification of potential biomarkers related to glioma survival by gene expression profile analysis.

Background

Methods

Results

Conclusion

Publications

Identification of potential biomarkers related to glioma survival by gene expression profile analysis.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets