Glioblastoma radiomics: can genomic and molecular characteristics correlate with imaging response patterns?

PURPOSE: For glioblastoma (GBM), imaging response (IR) or pseudoprogression (PSP) is frequently observed after chemoradiation and may connote a favorable prognosis. With tumors categorized by the Cancer Genome Atlas Project (mesenchymal, classical, neural, and proneural) and by methylguanine-methyltransferase (MGMT) methylation status, we attempted to determine if certain genomic or molecular subtypes of GBM were specifically associated with IR or PSP. METHODS: Patients with GBM treated at two institutions were reviewed. Kaplan-Meier method was used to estimate overall survival (OS) and progression-free survival (PFS). Mantel-cox test determined effect of IR and PSP on OS and PFS. Fisher's exact test was utilized to correlate IR and PSP with genomic subtypes and MGMT status. RESULTS: Eighty-two patients with GBM were reviewed. The median OS and PFS were 17.9 months and 8.9 months. IR was observed in 28 (40%) and was associated with improved OS (median 29.4 vs 14.5 months p < 0.01) and PFS (median 17.7 vs 5.5 months, p < 0.01). PSP was observed in 14 (19.2%) and trended towards improved PFS (15.0 vs 7.7 months p = 0.08). Tumors with a proneural component had a higher rate of IR compared to those without a proneural component (IR 60% vs 28%; p = 0.03). MGMT methylation was associated with IR (58% vs 24%, p = 0.032), but not PSP (34%, p = 0.10). CONCLUSION: IR is associated with improved OS and PFS. The proneural subtype and MGMT methylated tumors had higher rates of IR.

Agent-based computational modeling of glioblastoma predicts that stromal density is central to oncolytic virus efficacy.

Oncolytic viruses (OVs) are emerging cancer immunotherapy. Despite notable successes in the treatment of some tumors, OV therapy for central nervous system cancers has failed to show efficacy. We used an ex vivo tumor model developed from human glioblastoma tissue to evaluate the infiltration of herpes simplex OV rQNestin (oHSV-1) into glioblastoma tumors. We next leveraged our data to develop a computational, model of glioblastoma dynamics that accounts for cellular interactions within the tumor. Using our computational model, we found that low stromal density was highly predictive of oHSV-1 therapeutic success, suggesting that the efficacy of oHSV-1 in glioblastoma may be determined by stromal-to-tumor cell regional density. We validated these findings in heterogenous patient samples from brain metastatic adenocarcinoma. Our integrated modeling strategy can be applied to suggest mechanisms of therapeutic responses for central nervous system cancers and to facilitate the successful translation of OVs into the clinic.

Gelatin methacrylate hydrogels culture model for glioblastoma cells enriches for mesenchymal-like state and models interactions with immune cells.

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.

An anatomic transcriptional atlas of human glioblastoma.

Glioblastoma is an aggressive brain tumor that carries a poor prognosis. The tumor's molecular and cellular landscapes are complex, and their relationships to histologic features routinely used for diagnosis are unclear. We present the Ivy Glioblastoma Atlas, an anatomically based transcriptional atlas of human glioblastoma that aligns individual histologic features with genomic alterations and gene expression patterns, thus assigning molecular information to the most important morphologic hallmarks of the tumor. The atlas and its clinical and genomic database are freely accessible online data resources that will serve as a valuable platform for future investigations of glioblastoma pathogenesis, diagnosis, and treatment.

Cytoarchitecture of the spinal cord of the postnatal (P4) mouse.

Interpretation of the new wealth of gene expression and molecular mechanisms in the developing mouse spinal cord requires an accurate anatomical base on which data can be mapped. Therefore, we have assembled a spinal cord atlas of the P4 mouse to facilitate direct comparison with the adult specimens and to contribute to studies of the development of the mouse spinal cord. This study presents the anatomy of the spinal cord of the P4 C57Bl/6J mouse using Nissl and acetyl cholinesterase-stained sections. It includes a detailed map of the laminar organization of selected spinal cord segments and a description of named cell groups of the spinal cord such as the central cervical (CeCv), lateral spinal nucleus, lateral cervical, and dorsal nuclei. The motor neuron groups have also been identified according to the muscle groups they are likely to supply. General features of Rexed's laminae of the P4 spinal cord showed similarities to that of the adult (P56). However, certain differences were observed with regard to the extent of laminae and location of certain cell groups, such as the dorsal nucleus having a more dispersed structure and a more ventral and medial position or the CeCv being located in the medial part of lamina 5 in contrast to the adult where it is located in lamina 7. Motor neuron pools appeared to be more tightly packed in the P4 spinal cord. The dorsal horn was relatively larger and there was more white matter in the P56 spinal cord.

Adult spinal cord radial glia display a unique progenitor phenotype.

Radial glia (RG) are primarily embryonic neuroglial progenitors that express Brain Lipid Binding Protein (Blbp a.k.a. Fabp7) and Glial Fibrillary Acidic Protein (Gfap). We used these transcripts to demarcate the distribution of spinal cord radial glia (SCRG) and screen for SCRG gene expression in the Allen Spinal Cord Atlas (ASCA). We reveal that neonatal and adult SCRG are anchored in a non-ventricular niche at the spinal cord (SC) pial boundary, and express a "signature" subset of 122 genes, many of which are shared with "classic" neural stem cells (NSCs) of the subventricular zone (SVZ) and SC central canal (CC). A core expressed gene set shared between SCRG and progenitors of the SVZ and CC is particularly enriched in genes associated with human disease. Visualizing SCRG in a Fabp7-EGFP reporter mouse reveals an extensive population of SCRG that extend processes around the SC boundary and inwardly (through) the SC white matter (WM), whose abundance increases in a gradient from cervical to lumbar SC. Confocal analysis of multiple NSC-enriched proteins reveals that postnatal SCRG are a discrete and heterogeneous potential progenitor population that become activated by multiple SC lesions, and that CC progenitors are also more heterogeneous than previously appreciated. Gene ontology analysis highlights potentially unique regulatory pathways that may be further manipulated in SCRG to enhance repair in the context of injury and SC disease.

Genomic anatomy of the hippocampus.

Availability of genome-scale in situ hybridization data allows systematic analysis of genetic neuroanatomical architecture. Within the hippocampus, electrophysiology and lesion and imaging studies demonstrate functional heterogeneity along the septotemporal axis, although precise underlying circuitry and molecular substrates remain uncharacterized. Application of unbiased statistical component analyses to genome-scale hippocampal gene expression data revealed robust septotemporal molecular heterogeneity, leading to the identification of a large cohort of genes with robust regionalized hippocampal expression. Manual mapping of heterogeneous CA3 pyramidal neuron expression patterns demonstrates an unexpectedly complex molecular parcellation into a relatively coherent set of nine expression domains in the septal/temporal and proximal/distal axes with reciprocal, nonoverlapping boundaries. Unique combinatorial profiles of adhesion molecules within these domains suggest corresponding differential connectivity, which is demonstrated for CA3 projections to the lateral septum using retrograde labeling. This complex, discrete molecular architecture provides a novel paradigm for predicting functional differentiation across the full septotemporal extent of the hippocampus.