Probing glioblastoma and its microenvironment using singlenucleus and singlecell sequencing Singlecell (scSeq) and singlenucleus sequencing (snSeq) are powerful tools to investigate cancer genomics at single cell resolution. Multiple studies have recently illuminated intratumoral heterogeneity in glioblastoma, however, the majority focused on molecular complexity of tumor cells, without considering unexplored host cell types that contribute to the microenvironment around tumor. To address the glioblastoma microenvironment composition and potential tumorhost interactions, we performed deep coverage sequencing of freshly resected primary GBM patient tissue without implementing any tumor enrichment strategies. The sequencing resulted in 902 cells and 1186 nuclei, respectively, passing quality control and with low mitochondrial gene percentage. We customized reference transcriptome by listing gene transcript loci as exons to take into account immature RNA, which greatly improved the alignment rate for singlenucleus data. We applied Cell Ranger pipelines (Version 3.0.2) and Seurat package (Version 2.3.1) and discovered 10 clusters in both scSeq and snSeq. Pathway analysis of each cluster signature in scSeq data along with known GBM microenvironment cell signatures revealed glioma tumor population along with surrounding microgliamacrophages, astrocytes, pericytes, oligodendrocytes, T cells and endothelial cells. The analysis of snSeq was able to capture the majority of cell types from patient tissues (tumor and microenvironment cells), but interestingly presented different cell type composition in microenvironment cell types such as microgliamacrophages. Integrating singlecell and singlenucleus transcriptomic data using canonical correlation analysis facilitated a comparison of snSeq and scSeq, contrasting depiction for certain cell types (e.g. NKX62 gene in Oligodendrocytes). Differential analysis of pathways between tumor and microenvironment cells unveiled potentially rewired pathways such as double strand break repair pathway. Our results demonstrate the cellular diversity of brain tumor microenvironment and lay a foundation to further investigate the individual tumor and host cell transcriptomes that are influenced not only by their cell identity but also by their interaction with surrounding microenvironment.

Geostatistical visualization of ecological interactions in tumors Recent advances in our understanding of cancer progression have highlighted the roles played by molecular heterogeneity and by the tumor microenvironment in driving drug resistance and metastasis. The coupling of singlecell measurement technologies with algorithms, such as  t sne and SPADE, have enabled deep investigation of tumor heterogeneity. However, such techniques only capture molecular heterogeneity and do not enable the quantification nor visualization of intercellular interactions. They additionally do not allow the visualization of ecological niches that are critical to understanding tumor behavior. Novel computational tools to quantify and visualize spatial patterns in the tumor microenvironment are critically needed. Here, we take a tumor ecology perspective to examine how predation, mutualism, commensalism, and parasitism may impact tumor development and spatial patterning. We additionally quantify local spatial heterogeneity and the emergent global spatial behavior of the models using geostatistics. By visualizing emergent spatial patterns we demonstrate the potential utility of a geostatistical analysis in differentiating amongst cellcell interactions in the tumor microenvironment. These studies introduce both an ecological framework for characterizing intercellular interactions in cancer and a novel way of quantifying and visualizing spatial patterns in cancer.