A Cellular Ground Truth to Develop MRI Signatures in Glioma Models by Correlative Light Sheet Microscopy and Atlas-Based Coregistration.

Glioblastoma is the most common malignant primary brain tumor with poor overall survival. Magnetic resonance imaging (MRI) is the main imaging modality for glioblastoma but has inherent shortcomings. The molecular and cellular basis of MR signals is incompletely understood. We established a ground truth-based image analysis platform to coregister MRI and light sheet microscopy (LSM) data to each other and to an anatomic reference atlas for quantification of 20 predefined anatomic subregions. Our pipeline also includes a segmentation and quantification approach for single myeloid cells in entire LSM datasets. This method was applied to three preclinical glioma models in male and female mice (GL261, U87MG, and S24), which exhibit different key features of the human glioma. Multiparametric MR data including T2-weighted sequences, diffusion tensor imaging, T2 and T2* relaxometry were acquired. Following tissue clearing, LSM focused on the analysis of tumor cell density, microvasculature, and innate immune cell infiltration. Correlated analysis revealed differences in quantitative MRI metrics between the tumor-bearing and the contralateral hemisphere. LSM identified tumor subregions that differed in their MRI characteristics, indicating tumor heterogeneity. Interestingly, MRI signatures, defined as unique combinations of different MRI parameters, differed greatly between the models. The direct correlation of MRI and LSM allows an in-depth characterization of preclinical glioma and can be used to decipher the structural, cellular, and, likely, molecular basis of tumoral MRI biomarkers. Our approach may be applied in other preclinical brain tumor or neurologic disease models, and the derived MRI signatures could ultimately inform image interpretation in a clinical setting.SIGNIFICANCE STATEMENT We established a histologic ground truth-based approach for MR image analyses and tested this method in three preclinical glioma models exhibiting different features of glioblastoma. Coregistration of light sheet microscopy to MRI allowed for an evaluation of quantitative MRI data in histologically distinct tumor subregions. Coregistration to a mouse brain atlas enabled a regional comparison of MRI parameters with a histologically informed interpretation of the results. Our approach is transferable to other preclinical models of brain tumors and further neurologic disorders. The method can be used to decipher the structural, cellular, and molecular basis of MRI signal characteristics. Ultimately, information derived from such analyses could strengthen the neuroradiological evaluation of glioblastoma as they enhance the interpretation of MRI data.

Sphingosine-1-Phosphate Recruits Macrophages and Microglia and Induces a Pro-Tumorigenic Phenotype That Favors Glioma Progression.

Glioblastoma is the most aggressive brain tumor in adults. Treatment failure is predominantly caused by its high invasiveness and its ability to induce a supportive microenvironment. As part of this, a major role for tumor-associated macrophages/microglia (TAMs) in glioblastoma development was recognized. Phospholipids are important players in various fundamental biological processes, including tumor-stroma crosstalk, and the bioactive lipid sphingosine-1-phosphate (S1P) has been linked to glioblastoma cell proliferation, invasion, and survival. Despite the urgent need for better therapeutic approaches, novel strategies targeting sphingolipids in glioblastoma are still poorly explored. Here, we showed that higher amounts of S1P secreted by glioma cells are responsible for an active recruitment of TAMs, mediated by S1P receptor (S1PR) signaling through the modulation of Rac1/RhoA. This resulted in increased infiltration of TAMs in the tumor, which, in turn, triggered their pro-tumorigenic phenotype through the inhibition of NFkB-mediated inflammation. Gene set enrichment analyses showed that such an anti-inflammatory microenvironment correlated with shorter survival of glioblastoma patients. Inhibition of S1P restored a pro-inflammatory phenotype in TAMs and resulted in increased survival of tumor-bearing mice. Taken together, our results establish a crucial role for S1P in fine-tuning the crosstalk between glioma and infiltrating TAMs, thus pointing to the S1P-S1PR axis as an attractive target for glioma treatment.

Brain tumor classification of virtual NMR voxels based on realistic blood vessel-induced spin dephasing using support vector machines.

Remodeling of tissue microvasculature commonly promotes neoplastic growth; however, there is no imaging modality in oncology yet that noninvasively quantifies microvascular changes in clinical routine. Although blood capillaries cannot be resolved in typical magnetic resonance imaging (MRI) measurements, their geometry and distribution influence the integral nuclear magnetic resonance (NMR) signal from each macroscopic MRI voxel. We have numerically simulated the expected transverse relaxation in NMR voxels with different dimensions based on the realistic microvasculature in healthy and tumor-bearing mouse brains (U87 and GL261 glioblastoma). The 3D capillary structure in entire, undissected brains was acquired using light sheet fluorescence microscopy to produce large datasets of the highly resolved cerebrovasculature. Using this data, we trained support vector machines to classify virtual NMR voxels with different dimensions based on the simulated spin dephasing accountable to field inhomogeneities caused by the underlying vasculature. In prediction tests with previously blinded virtual voxels from healthy brain tissue and GL261 tumors, stable classification accuracies above 95% were reached. Our results indicate that high classification accuracies can be stably attained with achievable training set sizes and that larger MRI voxels facilitated increasingly successful classifications, even with small training datasets. We were able to prove that, theoretically, the transverse relaxation process can be harnessed to learn endogenous contrasts for single voxel tissue type classifications on tailored MRI acquisitions. If translatable to experimental MRI, this may augment diagnostic imaging in oncology with automated voxel-by-voxel signal interpretation to detect vascular pathologies.

Large-scale characterization of the microvascular geometry in development and disease by tissue clearing and quantitative ultramicroscopy.

Three-dimensional assessment of optically cleared, entire organs and organisms has recently become possible by tissue clearing and selective plane illumination microscopy ("ultramicroscopy"). Resulting datasets can be highly complex, encompass over a thousand images with millions of objects and data of several gigabytes per acquisition. This constitutes a major challenge for quantitative analysis. We have developed post-processing tools to quantify millions of microvessels and their distribution in three-dimensional datasets from ultramicroscopy and demonstrate the capabilities of our pipeline within entire mouse brains and embryos. Using our developed acquisition, segmentation, and analysis platform, we quantify physiological vascular networks in development and the healthy brain. We compare various geometric vessel parameters (e.g. vessel density, radius, tortuosity) in the embryonic spinal cord and brain as well as in different brain regions (basal ganglia, corpus callosum, cortex). White matter tract structures (corpus callosum, spinal cord) showed lower microvascular branch densities and longer vessel branch length compared to grey matter (cortex, basal ganglia). Furthermore, we assess tumor neoangiogenesis in a mouse glioma model to compare tumor core and tumor border. The developed methodology allows rapid quantification of three-dimensional datasets by semi-automated segmentation of fluorescently labeled objects with conventional computer hardware. Our approach can aid preclinical investigations and paves the way towards "quantitative ultramicroscopy".

A multi-omics analysis reveals the unfolded protein response regulon and stress-induced resistance to folate-based antimetabolites.

Stress response pathways are critical for cellular homeostasis, promoting survival through adaptive changes in gene expression and metabolism. They play key roles in numerous diseases and are implicated in cancer progression and chemoresistance. However, the underlying mechanisms are only poorly understood. We have employed a multi-omics approach to monitor changes to gene expression after induction of a stress response pathway, the unfolded protein response (UPR), probing in parallel the transcriptome, the proteome, and changes to translation. Stringent filtering reveals the induction of 267 genes, many of which have not previously been implicated in stress response pathways. We experimentally demonstrate that UPR-mediated translational control induces the expression of enzymes involved in a pathway that diverts intermediate metabolites from glycolysis to fuel mitochondrial one-carbon metabolism. Concomitantly, the cells become resistant to the folate-based antimetabolites Methotrexate and Pemetrexed, establishing a direct link between UPR-driven changes to gene expression and resistance to pharmacological treatment.

Gibbs point field model quantifies disorder in microvasculature of U87-glioblastoma.

Microvascular proliferation in glioblastoma multiforme is a biological key mechanism to facilitate tumor growth and infiltration and a main target for treatment interventions. The vascular architecture can be obtained by Single Plane Illumination Microscopy (SPIM) to evaluate vascular heterogeneity in tumorous tissue. We make use of the Gibbs point field model to quantify the order of regularity in capillary distributions found in the U87 glioblastoma model in a murine model and to compare tumorous and healthy brain tissue. A single model parameter Γ was assigned that is linked to tissue-specific vascular topology through Monte-Carlo simulations. Distributions of the model parameter Γ differ significantly between glioblastoma tissue with mean 〈ΓG〉=2.1±0.4, as compared to healthy brain tissue with mean 〈ΓH〉=4.9±0.4, suggesting that the average Γ-value allows for tissue differentiation. These results may be used for diagnostic magnetic resonance imaging, where it has been shown recently that Γ is linked to tissue-inherent relaxation parameters.

PERK-mediated expression of peptidylglycine α-amidating monooxygenase supports angiogenesis in glioblastoma.

PKR-like kinase (PERK) plays a significant role in inducing angiogenesis in various cancer types including glioblastoma. By proteomics analysis of the conditioned medium from a glioblastoma cell line treated with a PERK inhibitor, we showed that peptidylglycine α-amidating monooxygenase (PAM) expression is regulated by PERK under hypoxic conditions. Moreover, PERK activation via CCT020312 (a PERK selective activator) increased the cleavage and thus the generation of PAM cleaved cytosolic domain (PAM sfCD) that acts as a signaling molecule from the cytoplasm to the nuclei. PERK was also found to interact with PAM, suggesting a possible involvement in the generation of PAM sfCD. Knockdown of PERK or PAM reduced the formation of tubes by HUVECs in vitro. Furthermore, in vivo data highlighted the importance of PAM in the growth of glioblastoma with reduction of PAM expression in engrafted tumor significantly increasing the survival in mice. In summary, our data revealed PAM as a potential target for antiangiogenic therapy in glioblastoma.

Glioblastoma multiforme restructures the topological connectivity of cerebrovascular networks.

Glioblastoma multiforme alters healthy tissue vasculature by inducing angiogenesis and vascular remodeling. To fully comprehend the structural and functional properties of the resulting vascular network, it needs to be studied collectively by considering both geometric and topological properties. Utilizing Single Plane Illumination Microscopy (SPIM), the detailed capillary structure in entire healthy and tumor-bearing mouse brains could be resolved in three dimensions. At the scale of the smallest capillaries, the entire vascular systems of bulk U87- and GL261-glioblastoma xenografts, their respective cores, and healthy brain hemispheres were modeled as complex networks and quantified with fundamental topological measures. All individual vessel segments were further quantified geometrically and modular clusters were uncovered and characterized as meta-networks, facilitating an analysis of large-scale connectivity. An inclusive comparison of large tissue sections revealed that geometric properties of individual vessels were altered in glioblastoma in a relatively subtle way, with high intra- and inter-tumor heterogeneity, compared to the impact on the vessel connectivity. A network topology analysis revealed a clear decomposition of large modular structures and hierarchical network organization, while preserving most fundamental topological classifications, in both tumor models with distinct growth patterns. These results augment our understanding of cerebrovascular networks and offer a topological assessment of glioma-induced vascular remodeling. The findings may help understand the emergence of hypoxia and necrosis, and prove valuable for therapeutic interventions such as radiation or antiangiogenic therapy.

A sensitive and simple targeted proteomics approach to quantify transcription factor and membrane proteins of the unfolded protein response pathway in glioblastoma cells.

Many cellular events are driven by changes in protein expression, measurable by mass spectrometry or antibody-based assays. However, using conventional technology, the analysis of transcription factor or membrane receptor expression is often limited by an insufficient sensitivity and specificity. To overcome this limitation, we have developed a high-resolution targeted proteomics strategy, which allows quantification down to the lower attomol range in a straightforward way without any prior enrichment or fractionation approaches. The method applies isotope-labeled peptide standards for quantification of the protein of interest. As proof of principle, we applied the improved workflow to proteins of the unfolded protein response (UPR), a signaling pathway of great clinical importance, and could for the first time detect and quantify all major UPR receptors, transducers and effectors that are not readily detectable via antibody-based-, SRM- or conventional PRM assays. As transcription and translation is central to the regulation of UPR, quantification and determination of protein copy numbers in the cell is important for our understanding of the signaling process as well as how pharmacologic modulation of these pathways impacts on the signaling. These questions can be answered using our newly established workflow as exemplified in an experiment using UPR perturbation in a glioblastoma cell lines.