Dataset Information

Integration of RNA-Seq and proteomics data identifies glioblastoma multiforme surfaceome signature.

ABSTRACT:

Background

Glioblastoma 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.

Methods

We 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.

Results

Of 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.

Conclusions

In 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.

SUBMITTER: Syafruddin SE

PROVIDER: S-EPMC8306276 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

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:ObjectiveAlthough glucagon-secreting α-cells and insulin-secreting β-cells have opposing functions in regulating plasma glucose levels, the two cell types share a common developmental origin and exhibit overlapping transcriptomes and epigenomes. Notably, destruction of β-cells can stimulate repopulation via transdifferentiation of α-cells, at least in mice, suggesting plasticity between these cell fates. Furthermore, dysfunction of both α- and β-cells contributes to the pathophysiology of type 1 and type 2 diabetes, and β-cell de-differentiation has been proposed to contribute to type 2 diabetes. Our objective was to delineate the molecular properties that maintain islet cell type specification yet allow for cellular plasticity. We hypothesized that correlating cell type-specific transcriptomes with an atlas of open chromatin will identify novel genes and transcriptional regulatory elements such as enhancers involved in α- and β-cell specification and plasticity.MethodsWe sorted human α- and β-cells and performed the "Assay for Transposase-Accessible Chromatin with high throughput sequencing" (ATAC-seq) and mRNA-seq, followed by integrative analysis to identify cell type-selective gene regulatory regions.ResultsWe identified numerous transcripts with either α-cell- or β-cell-selective expression and discovered the cell type-selective open chromatin regions that correlate with these gene activation patterns. We confirmed cell type-selective expression on the protein level for two of the top hits from our screen. The "group specific protein" (GC; or vitamin D binding protein) was restricted to α-cells, while CHODL (chondrolectin) immunoreactivity was only present in β-cells. Furthermore, α-cell- and β-cell-selective ATAC-seq peaks were identified to overlap with known binding sites for islet transcription factors, as well as with single nucleotide polymorphisms (SNPs) previously identified as risk loci for type 2 diabetes.ConclusionsWe have determined the genetic landscape of human α- and β-cells based on chromatin accessibility and transcript levels, which allowed for detection of novel α- and β-cell signature genes not previously known to be expressed in islets. Using fine-mapping of open chromatin, we have identified thousands of potential cis-regulatory elements that operate in an endocrine cell type-specific fashion.

Project description: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.

Project description:BackgroundCoordinated efforts to collect large-scale data sets provide a basis for systems level understanding of complex diseases. In order to translate these fragmented and heterogeneous data sets into knowledge and medical benefits, advanced computational methods for data analysis, integration and visualization are needed.MethodsWe introduce a novel data integration framework, Anduril, for translating fragmented large-scale data into testable predictions. The Anduril framework allows rapid integration of heterogeneous data with state-of-the-art computational methods and existing knowledge in bio-databases. Anduril automatically generates thorough summary reports and a website that shows the most relevant features of each gene at a glance, allows sorting of data based on different parameters, and provides direct links to more detailed data on genes, transcripts or genomic regions. Anduril is open-source; all methods and documentation are freely available.ResultsWe have integrated multidimensional molecular and clinical data from 338 subjects having glioblastoma multiforme, one of the deadliest and most poorly understood cancers, using Anduril. The central objective of our approach is to identify genetic loci and genes that have significant survival effect. Our results suggest several novel genetic alterations linked to glioblastoma multiforme progression and, more specifically, reveal Moesin as a novel glioblastoma multiforme-associated gene that has a strong survival effect and whose depletion in vitro significantly inhibited cell proliferation. All analysis results are available as a comprehensive website.ConclusionsOur results demonstrate that integrated analysis and visualization of multidimensional and heterogeneous data by Anduril enables drawing conclusions on functional consequences of large-scale molecular data. Many of the identified genetic loci and genes having significant survival effect have not been reported earlier in the context of glioblastoma multiforme. Thus, in addition to generally applicable novel methodology, our results provide several glioblastoma multiforme candidate genes for further studies.Anduril is available at http://csbi.ltdk.helsinki.fi/anduril/The glioblastoma multiforme analysis results are available at http://csbi.ltdk.helsinki.fi/anduril/tcga-gbm/