HAVOC: Small-scale histomic mapping of cancer biodiversity across large tissue distances using deep neural networks.

Intratumoral heterogeneity can wreak havoc on current precision medicine strategies because of challenges in sufficient sampling of geographically separated areas of biodiversity distributed across centimeter-scale tumor distances. To address this gap, we developed a deep learning pipeline that leverages histomorphologic fingerprints of tissue to create "Histomic Atlases of Variation Of Cancers" (HAVOC). Using a number of objective molecular readouts, we demonstrate that HAVOC can define regional cancer boundaries with distinct biology. Using larger tumor specimens, we show that HAVOC can map biodiversity even across multiple tissue sections. By guiding profiling of 19 partitions across six high-grade gliomas, HAVOC revealed that distinct differentiation states can often coexist and be regionally distributed within these tumors. Last, to highlight generalizability, we benchmark HAVOC on additional tumor types. Together, we establish HAVOC as a versatile tool to generate small-scale maps of tissue heterogeneity and guide regional deployment of molecular resources to relevant biodiverse niches.

Niche deconvolution of the glioblastoma proteome reveals a distinct infiltrative phenotype within the proneural transcriptomic subgroup.

Glioblastoma is often subdivided into three transcriptional subtypes (classical, proneural, mesenchymal) based on bulk RNA signatures that correlate with distinct genetic and clinical features. Potential cellular-level differences of these subgroups, such as the relative proportions of glioblastoma's hallmark histopathologic features (e.g. brain infiltration, microvascular proliferation), may provide insight into their distinct phenotypes but are, however, not well understood. Here we leverage machine learning and reference proteomic profiles derived from micro-dissected samples of these major histomorphologic glioblastoma features to deconvolute and estimate niche proportions in an independent proteogenomically-characterized cohort. This approach revealed a strong association of the proneural transcriptional subtype with a diffusely infiltrating phenotype. Similarly, enrichment of a microvascular proliferation proteomic signature was seen within the mesenchymal subtype. This study is the first to link differences in the cellular pathology signatures and transcriptional profiles of glioblastoma, providing potential new insights into the genetic drivers and poor treatment response of specific subsets of glioblastomas.

The Brain Protein Atlas: A conglomerate of proteomics datasets of human neural tissue.

The human brain represents one of the most complex biological structures with significant spatiotemporal molecular plasticity occurring through early development, learning, aging, and disease. While much progress has been made in mapping its transcriptional architecture, more downstream phenotypic readouts are relatively scarce due to limitations with tissue heterogeneity and accessibility, as well as an inability to amplify protein species prior to global -OMICS analysis. To address some of these barriers, our group has recently focused on using mass-spectrometry workflows compatible with small amounts of formalin-fixed paraffin-embedded tissue samples. This has enabled exploration into spatiotemporal proteomic signatures of the brain and disease across otherwise inaccessible neurodevelopmental timepoints and anatomical niches. Given the similar theme and approaches, we introduce an integrated online portal, "The Brain Protein Atlas (BPA)" (www.brainproteinatlas.org), representing a public resource that allows users to access and explore these amalgamated datasets. Specifically, this portal contains a growing set of peer-reviewed mass-spectrometry-based proteomic datasets, including spatiotemporal profiles of human cerebral development, diffuse gliomas, clinically aggressive meningiomas, and a detailed anatomic atlas of glioblastoma. One barrier to entry in mass spectrometry-based proteomics data analysis is the steep learning curve required to extract biologically relevant data. BPA, therefore, includes several built-in analytical tools to generate relevant plots (e.g., volcano plots, heatmaps, boxplots, and scatter plots) and evaluate the spatiotemporal patterns of proteins of interest. Future iterations aim to expand available datasets, including those generated by the community at large, and analytical tools for exploration. Ultimately, BPA aims to improve knowledge dissemination of proteomic information across the neuroscience community in hopes of accelerating the biological understanding of the brain and various maladies.

Topographic mapping of the glioblastoma proteome reveals a triple-axis model of intra-tumoral heterogeneity.

Glioblastoma is an aggressive form of brain cancer with well-established patterns of intra-tumoral heterogeneity implicated in treatment resistance and progression. While regional and single cell transcriptomic variations of glioblastoma have been recently resolved, downstream phenotype-level proteomic programs have yet to be assigned across glioblastoma's hallmark histomorphologic niches. Here, we leverage mass spectrometry to spatially align abundance levels of 4,794 proteins to distinct histologic patterns across 20 patients and propose diverse molecular programs operational within these regional tumor compartments. Using machine learning, we overlay concordant transcriptional information, and define two distinct proteogenomic programs, MYC- and KRAS-axis hereon, that cooperate with hypoxia to produce a tri-dimensional model of intra-tumoral heterogeneity. Moreover, we highlight differential drug sensitivities and relative chemoresistance in glioblastoma cell lines with enhanced KRAS programs. Importantly, these pharmacological differences are less pronounced in transcriptional glioblastoma subgroups suggesting that this model may provide insights for targeting heterogeneity and overcoming therapy resistance.

Defining Protein Pattern Differences Among Molecular Subtypes of Diffuse Gliomas Using Mass Spectrometry.

Molecular characterization of diffuse gliomas has thus far largely focused on genomic and transcriptomic interrogations. Here, we utilized mass spectrometry and overlay protein-level information onto genomically defined cohorts of diffuse gliomas to improve our downstream molecular understanding of these lethal malignancies. Bulk and macrodissected tissues were utilized to quantitate 5,496 unique proteins over three glioma cohorts subclassified largely based on their IDH and 1p19q codeletion status (IDH wild type (IDHwt), n = 7; IDH mutated (IDHmt), 1p19q non-codeleted, n = 7; IDH mutated, 1p19q-codeleted, n = 10). Clustering analysis highlighted proteome and systems-level pathway differences in gliomas according to IDH and 1p19q-codeletion status, including 287 differentially abundant proteins in macrodissection-enriched tumor specimens. IDHwt tumors were enriched for proteins involved in invasiveness and epithelial to mesenchymal transition (EMT), while IDHmt gliomas had increased abundances of proteins involved in mRNA splicing. Finally, these abundance changes were compared with IDH-matched GBM stem-like cells (GSCs) to better pinpoint protein patterns enriched in putative cellular drivers of gliomas. Using this integrative approach, we outline specific proteins involved in chloride transport (e.g. chloride intracellular channel 1, CLIC1) and EMT (e.g. procollagen-lysine, 2-oxoglutarate 5-dioxygenase 3, PLOD3, and serpin peptidase inhibitor clade H member 1, SERPINH1) that showed concordant IDH-status-dependent abundance differences in both primary tissue and purified GSC cultures. Given the downstream position proteins occupy in driving biology and phenotype, understanding the proteomic patterns operational in distinct glioma subtypes could help propose more specific, personalized, and effective targets for the management of patients with these aggressive malignancies.