No abstract

Glioblastoma is the most lethal primary malignant brain tumor in adults. Simplified two-dimensional (2D) cell culture and neurospheres in vitro models fail to recapitulate the complexity of the tumor microenvironment, limiting its ability to predict therapeutic response. Three-dimensional (3D) scaffold-based models have emerged as a promising alternative for addressing these concerns. One such 3D system is gelatin methacrylate (GelMA) hydrogels, and we aimed to understand the suitability of using this system to mimic treatment-resistant glioblastoma cells that reside in specific niches. We characterized the phenotype of patient-derived glioma cells cultured in GelMA hydrogels (3D-GMH) for their tumorigenic properties using invasion and chemoresponse assays. In addition, we used integrated single-cell and spatial transcriptome analysis to compare cells cultured in 3D-GMH to neoplastic cells in vivo. Finally, we assessed tumor-immune cell interactions with a macrophage infiltration assay and a cytokine array. We show that the 3D-GMH system enriches treatment-resistant mesenchymal cells that are not represented in neurosphere cultures. Cells cultured in 3D-GMH resemble a mesenchymal-like cellular phenotype found in perivascular and hypoxic regions and recruit macrophages by secreting cytokines, a hallmark of the mesenchymal phenotype. Our 3D-GMH model effectively mimics the phenotype of glioma cells that are found in the perivascular and hypoxic niches of the glioblastoma core in situ, in contrast to the neurosphere cultures that enrich cells of the infiltrative edge of the tumor. This contrast highlights the need for due diligence in selecting an appropriate model when designing a study‘s objectives.

We have studied differentially regulated nuclear proteome of the clinical tissue specimens of glioblastoma (GBM, WHO Grade IV) and lower grades of gliomas (Grade II and III) using high resolution mass spectrometry- based quantitative proteomics approach. The results showed altered expression of many regulatory proteins from the nucleus such as DNA binding proteins, transcription and post transcriptional processing factors and also included enrichment of nuclear proteins that are targets of granzyme signaling – an immune surveillance pathway. Protein - protein interaction network analysis using integrated proteomics and transcriptomics data of transcription factors and proteins for cell invasion process (drawn from another GBM dataset) revealed YBX1, a ubiquitous RNA and DNA-binding protein and a transcription factor, as a key interactor of major cell invasion-associated proteins from GBM. To verify the regulatory link between them, the co-expression of YBX1 and six of the interacting proteins (EGFR, MAPK1, CD44, SOX2, TNC and MMP13) involved in cell invasion network was examined by immunohistochemistry on tissue micro arrays. Our analysis suggests YBX1 as a potential regulator of these key molecules involved in tumor invasion and thus as a promising target for development of new therapeutic strategies for GBM.

G protein-coupled receptor 56 (GPR56/ADGRG1) is an adhesion GPCR with an essential role in brain development and cancer. Elevated expression of GPR56 was observed in the clinical specimens of Glioblastoma (GBM), a highly invasive primary brain tumor. However, we found the expression to be variable across the specimens, presumably due to the intratumor heterogeneity of GBM. Therefore, we re-examined GPR56 expression in public domain spatial gene expression data and single-cell expression data for GBM, which revealed that GPR56 expression was high in cellular tumors, infiltrating tumor cells, and proliferating cells, low in microvascular proliferation and peri-necrotic areas of the tumor, especially in hypoxic mesenchymal-like cells. To gain a better understanding of the consequences of GPR56 downregulation in tumor cells and other molecular changes associated with it, we generated a sh-RNA-mediated GPR56 knockdown in the GBM cell line U373 and performed transcriptomics, proteomics, and phospho-proteomics analysis. Our analysis revealed enrichment of gene signatures, pathways, and phosphorylation of proteins potentially associated with mesenchymal (MES) transition in the tumor and concurrent increase in cell invasion and migration behavior of the GPR56 knockdown GBM cells. Interestingly, our analysis also showed elevated expression of Transglutaminase 2 (TG2) - a known interactor of GPR56, in the knockdown cells. The inverse expression of GPR56 and TG2 was also observed in intratumoral, spatial gene expression data for GBM and in GBM cell lines cultured in vitro under hypoxic conditions. Integrating all these observations, we propose a putative functional link between the inverse expression of the two proteins, the hypoxic niche and the mesenchymal status in the tumor. Hypoxia-induced downregulation of GPR56 and activation of TG2 may result in a network of molecular events that contribute to the mesenchymal transition of GBM cells, and we propose a putative model to explain this functional and regulatory relationship of the two proteins.

No abstract

Gliomas are heavily infiltrated with immune cells of myeloid origin. Past studies have shown that high-grade gliomas have a higher proportion of alternatively activated and suppressive myeloid cells when compared to low-grade gliomas, which correlate with poor prognosis. However, the differences in immune cell phenotypes within high-grade gliomas (between grade 3 and grade 4 or GBM) are relatively less explored, and a correlation of phenotypic characteristics between immune cells in the blood and high-grade tumors has not been performed. Additionally, myeloid cells of granulocytic origin present in gliomas remain poorly characterized. Herein, we address these questions through phenotypic characterizations of monocytes and neutrophils present in blood and tumors of individuals with glioblastoma (GBM, IDH-wild type) or grade 3 IDH-mutant gliomas. We observe that neutrophils are highly heterogeneous among individuals with glioma, and are different from healthy controls. We also show that CD163 expressing M2 monocytes are present in greater proportions in GBM tissue when compared to grade 3 IDH-mutant glioma tissue, and a larger proportion of granulocytic myeloid-derived suppressor cells are present in grade 3 IDH-mutant gliomas when compared to GBM. Finally, we demonstrate that the expression levels of CD86 and CD63 showed a high correlation between blood and tumor and suggest that these may be used as possible markers for prognosis.

Splice variants are known to be important in the pathophysiology of tumors, including the brain cancers. We applied a proteogenomics pipeline to identify splice variants in glioblastoma (GBM, grade IV glioma), a highly malignant brain tumor, using in-house generated mass spectrometric proteomic data and public domain RNASeq dataset. Our analysis led to the identification of a novel exon that maps to the long isoform of Neural cell adhesion molecule 1 (NCAM1), expressed on the surface of glial cells and neurons, important for cell adhesion and cell signaling. The presence of the novel exon is supported with the identification of five peptides spanning it. Additional peptides were also detected in sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) gel separated proteins from GBM patient tissue, underscoring the presence of the novel peptides in the intact brain protein. The novel exon was detected in the RNASeq dataset in 18 of 25 GBM samples and separately validated in additional 10 GBM tumor tissues using quantitative real-time-polymerase chain reaction (qRT-PCR). Both transcriptomic and proteomic data indicate downregulation of NCAM1, including the novel variant, in GBM. Domain analysis of the novel NCAM1 sequence indicates that the insertion of the novel exon contributes extra low-complexity region in the protein that may be important for protein-protein interactions and hence for cell signaling associated with tumor development. Taken together, the novel NCAM1 variant reported in this study exemplifies the importance of future multiomics research and systems biology applications in GBM.

We propose a novel strategy for discovering motifs from gene expression data. The gene expression data in our experiments comes from DNA Microarray analysis of the bacterium E. coli in response to recovery from nutrient starvation. We have annotated the data and identified the upregulated genes. Our interest is to find common regulatory motifs that are responsible for the upregulation of these specific genes. We assume that a common motif that a regulatory protein can bind to will be present in the upstream region of the upregulated genes and will not be present in the upstream regions of genes that showed a constant level of expression over time. Our objective is to find the common motifs that are present in at least some of the upstream sequences of upregulated genes and not present in the control set, which is the set of genes whose expression remained the same. Because it is possible that there could be several subsets of co-regulated genes under different control mechanisms among the co-expressed genes, we do not want to require motifs to be present in all upregulated sequences. Therefore, we propose a new algorithm for finding such motifs through stages of pre-processing, denoising, agglomerative clustering and consensus checking. Through this process, we have found some motifs that are good candidates for further validation.