Dataset Information

Identification of hub genes and regulatory factors of glioblastoma multiforme subgroups by RNA-seq data analysis.

ABSTRACT: Glioblastoma multiforme (GBM) is the most common malignant brain tumor. This study aimed to identify the hub genes and regulatory factors of GBM subgroups by RNA sequencing (RNA-seq) data analysis, in order to explore the possible mechanisms responsbile for the progression of GBM. The dataset RNASeqV2 was downloaded by TCGA-Assembler, containing 169 GBM and 5 normal samples. Gene expression was calculated by the reads per kilobase per million reads measurement, and nor malized with tag count comparison. Following subgroup classification by the non-negative matrix factorization, the differentially expressed genes (DEGs) were screened in 4 GBM subgroups using the method of significance analysis of microarrays. Functional enrichment analysis was performed by DAVID, and the protein-protein interaction (PPI) network was constructed based on the HPRD database. The subgroup-related microRNAs (miRNAs or miRs), transcription factors (TFs) and small molecule drugs were predicted with pre-defined criteria. A cohort of 19,515 DEGs between the GBM and control samples was screened, which were predominantly enriched in cell cycle- and immunoreaction-related pathways. In the PPI network, lymphocyte cytosolic protein 2 (LCP2), breast cancer 1 (BRCA1), specificity protein 1 (Sp1) and chromodomain-helicase-DNA-binding protein 3 (CHD3) were the hub nodes in subgroups 1-4, respectively. Paired box 5 (PAX5), adipocyte protein 2 (aP2), E2F transcription factor 1 (E2F1) and cAMP-response element-binding protein-1 (CREB1) were the specific TFs in subgroups 1-4, respectively. miR‑147b, miR‑770-5p, miR‑220a and miR‑1247 were the particular miRNAs in subgroups 1-4, respectively. Natalizumab was the predicted small molecule drug in subgroup 2. In conclusion, the molecular regulatory mechanisms of GBM pathogenesis were distinct in the different subgroups. Several crucial genes, TFs, miRNAs and small molecules in the different GBM subgroups were identified, which may be used as potential markers. However, further experimental validations may be required.

SUBMITTER: Li Y

PROVIDER: S-EPMC5029949 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundGlioblastoma multiforme (GBM) is a highly lethal, stage IV brain tumour with a prevalence of approximately 2 per 10,000 people globally. The cell surface proteins or surfaceome serve as information gateway in many oncogenic signalling pathways and are important in modulating cancer phenotypes. Dysregulation in surfaceome expression and activity have been shown to promote tumorigenesis. The expression of GBM surfaceome is a case in point; OMICS screening in a cell-based system identified that this sub-proteome is largely perturbed in GBM. Additionally, since these cell surface proteins have 'direct' access to drugs, they are appealing targets for cancer therapy. However, a comprehensive GBM surfaceome landscape has not been fully defined yet. Thus, this study aimed to define GBM-associated surfaceome genes and identify key cell-surface genes that could potentially be developed as novel GBM biomarkers for therapeutic purposes.MethodsWe integrated the RNA-Seq data from TCGA GBM (n = 166) and GTEx normal brain cortex (n = 408) databases to identify the significantly dysregulated surfaceome in GBM. This was followed by an integrative analysis that combines transcriptomics, proteomics and protein-protein interaction network data to prioritize the high-confidence GBM surfaceome signature.ResultsOf the 2381 significantly dysregulated genes in GBM, 395 genes were classified as surfaceome. Via the integrative analysis, we identified 6 high-confidence GBM molecular signature, HLA-DRA, CD44, SLC1A5, EGFR, ITGB2, PTPRJ, which were significantly upregulated in GBM. The expression of these genes was validated in an independent transcriptomics database, which confirmed their upregulated expression in GBM. Importantly, high expression of CD44, PTPRJ and HLA-DRA is significantly associated with poor disease-free survival. Last, using the Drugbank database, we identified several clinically-approved drugs targeting the GBM molecular signature suggesting potential drug repurposing.ConclusionsIn summary, we identified and highlighted the key GBM surface-enriched repertoires that could be biologically relevant in supporting GBM pathogenesis. These genes could be further interrogated experimentally in future studies that could lead to efficient diagnostic/prognostic markers or potential treatment options for GBM.

Project description:Glioblastoma (GBM) represents a profoundly aggressive and heterogeneous brain neoplasm linked to a bleak prognosis. Hypoxia, a common feature in GBM, has been linked to tumor progression and therapy resistance. In this study, we aimed to identify hypoxia-related differentially expressed genes (DEGs) and construct a prognostic signature for GBM patients using multi-omics analysis. Patient cohorts were collected from publicly available databases, including the Gene Expression Omnibus (GEO), the Chinese Glioma Genome Atlas (CGGA), and The Cancer Genome Atlas-Glioblastoma Multiforme (TCGA-GBM), to facilitate a comprehensive analysis. Hypoxia-related genes (HRGs) were obtained from the Molecular Signatures Database (MSigDB). Differential expression analysis revealed 41 hypoxia-related DEGs in GBM patients. A consensus clustering approach, utilizing these DEGs' expression patterns, identified four distinct clusters, with cluster 1 showing significantly better overall survival. Machine learning techniques, including univariate Cox regression and LASSO regression, delineated a prognostic signature comprising six genes (ANXA1, CALD1, CP, IGFBP2, IGFBP5, and LOX). Multivariate Cox regression analysis substantiated the prognostic significance of a set of three optimal signature genes (CP, IGFBP2, and LOX). Using the hypoxia-related prognostic signature, patients were classified into high- and low-risk categories. Survival analysis demonstrated that the high-risk group exhibited inferior overall survival rates in comparison to the low-risk group. The prognostic signature showed good predictive performance, as indicated by the area under the curve (AUC) values for one-, three-, and five-year overall survival. Furthermore, functional enrichment analysis of the DEGs identified biological processes and pathways associated with hypoxia, providing insights into the underlying mechanisms of GBM. Delving into the tumor immune microenvironment, our analysis revealed correlations relating the hypoxia-related prognostic signature to the infiltration of immune cells in GBM. Overall, our study highlights the potential of a hypoxia-related prognostic signature as a valuable resource for forecasting the survival outcome of GBM patients. The multi-omics approach integrating bulk sequencing, single-cell analysis, and immune microenvironment assessment enhances our understanding of the intricate biology characterizing GBM, thereby potentially informing the tailored design of therapeutic interventions.

Project description:In addition to direct targeting and repressing mRNAs, recent studies reported that microRNAs (miRNAs) can bridge up an alternative layer of post-transcriptional gene regulatory networks. The competing endogenous RNA (ceRNA) regulation depicts the scenario where pairs of genes (ceRNAs) sharing, fully or partially, common binding miRNAs (miRNA program) can establish coexpression through competition for a limited pool of the miRNA program. While the dynamics of ceRNA regulation among cellular conditions have been verified based on in silico and in vitro experiments, comprehensive investigation into the strength of ceRNA regulation in human datasets remains largely unexplored. Furthermore, pan-cancer analysis of ceRNA regulation, to our knowledge, has not been systematically investigated.In the present study we explored optimal conditions for ceRNA regulation, investigated functions governed by ceRNA regulation, and evaluated pan-cancer effects. We started by investigating how essential factors, such as the size of miRNA programs, the number of miRNA program binding sites, and expression levels of miRNA programs and ceRNAs affect the ceRNA regulation capacity in tumors derived from glioblastoma multiforme patients captured by The Cancer Genome Atlas (TCGA). We demonstrated that increased numbers of common targeting miRNAs as well as the abundance of binding sites enhance ceRNA regulation and strengthen coexpression of ceRNA pairs. Also, our investigation revealed that the strength of ceRNA regulation is dependent on expression levels of both miRNA programs and ceRNAs. Through functional annotation analysis, our results indicated that ceRNA regulation is highly associated with essential cellular functions and diseases including cancer. Furthermore, the highly intertwined ceRNA regulatory relationship enables constitutive and effective intra-function regulation of genes in diverse types of cancer.Using gene and microRNA expression datasets from TCGA, we successfully quantified the optimal conditions for ceRNA regulation, which hinge on four essential parameters of ceRNAs. Our analysis suggests optimized ceRNA regulation is related to disease pathways and essential cellular functions. Furthermore, although the strength of ceRNA regulation is dynamic among cancers, its governing functions are stably maintained. The findings of this report contribute to better understanding of ceRNA dynamics and its crucial roles in cancers.

Dataset Information

Identification of hub genes and regulatory factors of glioblastoma multiforme subgroups by RNA-seq data analysis.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets