Bulk transcriptomic classifiers stratify pancreatic ductal adenocarcinoma (PDAC) into classical and basal-like subtypes with prognostic and therapeutic relevance, yet increasing evidence indicates that these epithelial programs frequently coexist within individual tumors. How these intermediate states affect the local tumor microenvironment remains poorly defined. Here, we integrate multiplexed ion beam imaging (MIBI) with bulk RNA sequencing to resolve epithelial subtype identity at the level of spatially contiguous cancer nests and quantify their associated microenvironments. Across 47 primary tumor samples from 34 patients, we identified classical, intermediate, and basal cancer cell states at single-cell resolution and delineated discrete cancer nests with mixed or dominant subtype compositions. Distance-resolved spatial analysis reveals that basal-rich cancer nests are surrounded by locally immunosuppressive microenvironments characterized by reduced expression of MHC class II and co-stimulatory molecules in proximal myeloid cells, independent of myeloid abundance. These regions are enriched in fibroblast-dominated neighborhoods and distinct cell-cell interaction architectures. Using EcoTyper analysis of two independent bulk RNA-seq cohorts, including an OHSU discovery cohort (N = 277 patients) and TCGA as a validation cohort (N = 147 patients), we identified poor-prognosis tumor ecotypes enriched for basal epithelial states that similarly exhibited depleted myeloid antigen presentation signatures, linking spatial niche phenotypes to transcriptional ecotypes and patient outcomes. Together, these findings demonstrate that epithelial subtype programs in PDAC are organized at the level of spatially defined cancer nests and that basal cancer programs reside within localized niches of myeloid antigen presentation dysfunction, linking intratumoral architecture to immune suppression and clinical prognosis.

Understanding the role of tertiary lymphoid structures (TLS) is crucial in non-small cell lung cancer (NSCLC), as they are associated with patient prognosis and treatment outcomes. Specific cellular ecosystems that originate anti-tumor activity or predict immunotherapy response remain poorly characterized. To this end, we developed a high-resolution, multimodal spatial atlas jointly profiling transcriptomics, proteomics, and histology to characterize TLS maturation in NSCLC alongside secondary lymph organs as a baseline. Using this atlas, we proposed a pathologist-in-the-loop framework that combines a variational graph autoencoder (VGAE) with diffusion pseudotime to refine human expert annotations and characterize TLS maturation. These spatial molecular representations were extended to H&E whole-slide images via a vision transformer-based foundation model. Next, we resolved cellular composition, spatial organization, and cell-cell interactions within these data and defined two divergent spatial ecosystems. Clinical evidence suggests that these ecosystems are associated with distinct patient outcomes: a mature germinal center niche with favorable prognoses, and a tumor-macrophage-fibroblast niche with unfavorable prognoses. In summary, our work decodes key components of TLS heterogeneity, identifies hallmark spatial patterns involved in NSCLC adaptive immunity, and provides a framework for translating spatial omics insights into clinical applications. ### Competing Interest Statement The authors have declared no competing interest.

Immune checkpoint inhibition (ICI) benefits only a subset of patients with metastatic triple-negative breast cancer and determinants of response remain unclear. We assembled a longitudinal cohort of 103 female patients from the phase 2 TONIC trial, with samples spanning primary tumors, pretreatment metastases and on-treatment metastases during nivolumab therapy. We profiled 37 proteins in 270 tumors using highly multiplexed imaging and developed SpaceCat, an open-source pipeline that extracts more than 800 imaging features per sample, including cell density, diversity, spatial interactions and functional marker expression. Metastatic but not primary tumors contained features predictive of outcome. Spatial metrics such as immune diversity and T cell infiltration at tumor borders were most informative, while ratios of T cells to cancer cells and PDL1 on myeloid cells were also associated with response. Multivariate models stratified patients with the highest performance on treatment (area under the curve = 0.90). Bulk RNA-seq confirmed the predictive value of on-treatment samples. These findings highlight the value of longitudinal profiling to resolve evolving tumor microenvironment dynamics driving ICI response.

Abstract Tertiary lymphoid structures (TLS) predict benefit from immune checkpoint inhibitors (CPIs), yet mature, germinal-center-rich TLS are infrequent in solid tumors by histological review. Here, using 38-plex MIBI spatial proteomics across 165 lymphoid structures from 14 NSCLC resections, we establish a continuum of TLS maturity using high dimensional compositional, spatial, and molecular features. We demonstrate that histologically-defined lymphoid aggregates (LA) comprise a heterogeneous class of structures, which span this continuum of maturity. We identify a subset of lymphoid aggregates that harbor follicular dendritic cell networks, T follicular helper cells, and activated B cell states characteristic of mature TLS, yet are not readily distinguished from other LA structures in our histological review. We developed a novel digital pathology classifier to identify mature LAs in CPI trials, and demonstrate in a retrospective analysis of Atezolizumab in advanced NSCLC that the inclusion of mature LAs greatly expands the biomarker-eligible population while maintaining strong predicted benefit. Together, these data redefine the biological spectrum of tumor-associated lymphoid aggregates and provide a framework for implementing maturity-informed TLS biomarker strategies.

Immune checkpoint inhibition (ICI) has fundamentally changed cancer treatment. However, only a minority of patients with metastatic triple negative breast cancer (TNBC) benefit from ICI, and the determinants of response remain largely unknown. To better understand the factors influencing patient outcome, we assembled a longitudinal cohort with tissue from multiple timepoints, including primary tumor, pre-treatment metastatic tumor, and on-treatment metastatic tumor from 117 patients treated with ICI (nivolumab) in the phase II TONIC trial. We used highly multiplexed imaging to quantify the subcellular localization of 37 proteins in each tumor. To extract meaningful information from the imaging data, we developed SpaceCat, a computational pipeline that quantifies features from imaging data such as cell density, cell diversity, spatial structure, and functional marker expression. We applied SpaceCat to 678 images from 294 tumors, generating more than 800 distinct features per tumor. Spatial features were more predictive of patient outcome, including features like the degree of mixing between cancer and immune cells, the diversity of the neighboring immune cells surrounding cancer cells, and the degree of T cell infiltration at the tumor border. Non-spatial features, including the ratio between T cell subsets and cancer cells and PD-L1 levels on myeloid cells, were also associated with patient outcome. Surprisingly, we did not identify robust predictors of response in the primary tumors. In contrast, the metastatic tumors had numerous features which predicted response. Some of these features, such as the cellular diversity at the tumor border, were shared across timepoints, but many of the features, such as T cell infiltration at the tumor border, were predictive of response at only a single timepoint. We trained multivariate models on all of the features in the dataset, finding that we could accurately predict patient outcome from the pre-treatment metastatic tumors, with improved performance using the on-treatment tumors. We validated our findings in matched bulk RNA-seq data, finding the most informative features from the on-treatment samples. Our study highlights the importance of profiling sequential tumor biopsies to understand the evolution of the tumor microenvironment, elucidating the temporal and spatial dynamics underlying patient responses and underscoring the need for further research on the prognostic role of metastatic tissue and its utility in stratifying patients for ICI.

Recent advancements in transcriptomics and proteomics have opened the possibility for spatially resolved molecular characterization of tissue architecture with the promise of enabling a deeper understanding of tissue biology in either homeostasis or disease. The wealth of data generated by these technologies has recently driven the development of a wide range of computational methods. These methods have the requirement of advanced coding fluency to be applied and integrated across the full spatial omics analysis process, thus presenting a hurdle for widespread adoption by the biology research community. To address this, we introduce SPEX (Spatial Expression Explorer), a web-based analysis platform that employs modular analysis pipeline design, accessible through a user-friendly interface. SPEX's infrastructure allows for streamlined access to open-source image data management systems, analysis modules, and fully integrated data visualization solutions. Analysis modules include essential steps covering image processing, single-cell analysis, and spatial analysis. We demonstrate SPEX's ability to facilitate the discovery of biological insights in spatially resolved omics datasets from healthy tissue to tumor samples.

Spatial profiling of tissues promises to elucidate tumor-microenvironment interactions and generate prognostic and predictive biomarkers. We analyzed single-cell spatial data from 3 multiplex imaging technologies: cyclic immunofluorescence (CycIF) data we generated from 102 patients with breast cancer with clinical follow-up as well as publicly available mass cytometry and multiplex ion-beam imaging datasets. Similar single-cell phenotyping results across imaging platforms enabled combined analysis of epithelial phenotypes to delineate prognostic subtypes among patients who are estrogen-receptor+ (ER+). We utilized discovery and validation cohorts to identify biomarkers with prognostic value. Increased lymphocyte infiltration was independently associated with longer survival in triple-negative (TN) and high-proliferation ER+ breast tumors. An assessment of 10 spatial analysis methods revealed robust spatial biomarkers. In ER+ disease, quiescent stromal cells close to tumor were abundant in tumors with good prognoses, while tumor cell neighborhoods containing mixed fibroblast phenotypes were enriched in poor-prognosis tumors. In TN disease, macrophage/tumor and B/T lymphocyte neighbors were enriched, and lymphocytes were dispersed in good-prognosis tumors, while tumor cell neighborhoods containing vimentin+ fibroblasts were enriched in poor-prognosis tumors. In conclusion, we generated comparable single-cell spatial proteomic data from several clinical cohorts to enable prognostic spatial biomarker identification and validation.

Despite recent advances, the adoption of computer vision methods into clinical and commercial applications has been hampered by the limited availability of accurate ground truth tissue annotations required to train robust supervised models. Generating such ground truth can be accelerated by annotating tissue molecularly using immunofluorescence (IF) staining and mapping these annotations to a post-IF hematoxylin and eosin (H&E) (terminal H&E) stain. Mapping the annotations between IF and terminal H&E increases both the scale and accuracy by which ground truth could be generated. However, discrepancies between terminal H&E and conventional H&E caused by IF tissue processing have limited this implementation. We sought to overcome this challenge and achieve compatibility between these parallel modalities using synthetic image generation, in which a cycle-consistent generative adversarial network was applied to transfer the appearance of conventional H&E such that it emulates terminal H&E. These synthetic emulations allowed us to train a deep learning model for the segmentation of epithelium in terminal H&E that could be validated against the IF staining of epithelial-based cytokeratins. The combination of this segmentation model with the cycle-consistent generative adversarial network stain transfer model enabled performative epithelium segmentation in conventional H&E images. The approach demonstrates that the training of accurate segmentation models for the breadth of conventional H&E data can be executed free of human expert annotations by leveraging molecular annotation strategies such as IF, so long as the tissue impacts of the molecular annotation protocol are captured by generative models that can be deployed prior to the segmentation process.

PDF file - 954K, Supplemental Figures: Figure S1. Chemical structure of the Akt inhibitor GDC-0068. Figure S2. Effect of GDC-0068 on cell cycle and apoptotic response on PC-3 and MCF7-neo/HER2 cell lines. Figure S3. Western blot analysis of indicated proteins in isogenic MCF10A cells with and without PTEN knockout in the presence of 20 or 0.2 ng/ml EGF. Figure S4. Representative immunohistochemistry staining images of cleaved caspase-3 in TOV-21G.x1 tumors. Figure S5. Combination effects between GDC-0068 and chemotherapeutic agents in vitro.

Supplemental Figures S1-S9. S1: Kinase inhibitor screen and response of OP449 and INK128. S2: CI values for DT1154 combined with INK128 and viability curves. S3: Kmeans clusters and signaling in response to PP2A activators and INK128. S4: Oncoprint and survival curves of patients with alterations in the MYC/mTOR pathway. S5: pT58 MYC expression levels in PDA cells treated with DT1154 and INK128. S6: Signaling with B56a knockdown or T58A MYC overexpression in response to DT1154 and INK128 treatment. S7: Signaling showing decreased expression of MYC increases the efficacy of INK128. S8: Treatment of KPC and KPCM/+ cell lines with BEZ235, Dasatinib, and Sunitinib. S9: Endpoint analysis of long- and short-term in vivo treatment with DT1154 and INK128.

XLSX file - 74K, Table S1. Enzymatic potency, selectivity and cellular potency of GDC-0068 Table S2. GDC-0068 IC50 profile on cell viability and genetic background of cancer cell lines Table S3. Maximum percent body weight changes of mice treated with GDC-0068 single agent and in combination with chemotherapeutic agents

Multiplex ion beam imaging (MIBI) and imaging mass cytometry (IMC) enable highly multiplexed antibody (40+) staining of frozen or formalin fixed, paraffin-embedded (FFPE) human or murine tissues through detection of metal ions liberated from primary antibodies by time-of-flight mass spectrometry (TOF). These methods make detection of more than 50 targets theoretically possible while maintaining spatial orientation. As such, they are ideal tools to identify the multiple immune, epithelial, and stromal cell subsets in the tumor microenvironment and to characterize spatial relationships and tumor-immune status in either murine models or human samples. This chapter summarizes methods for antibody conjugation and validation, staining, and preliminary data collection using IMC or MIBI in both human and mouse pancreatic adenocarcinoma samples. These protocols are intended to facilitate use of these complex platforms in not only tissue-based tumor immunology studies but also tissue-based oncology or immunology studies more broadly.

Results of Kinase Inhibitor Screen including compound names, IC50s, and fold change compared to vehicle.

Results of the Phosphokinase Array with phosphorylation sites, fold change relative to vehicle, and Kmeans clustering.

Recent advancements in transcriptomics and proteomics have opened the possibility for spatially resolved molecular characterization of tissue architecture with the promise of enabling a deeper understanding of tissue biology in either homeostasis or disease. The wealth of data generated by these technologies has recently driven the development of computational pipelines that, nevertheless, have the requirement of coding fluency to be applied. To remove this hurdle, we present SPEX (Spatial Expression Explorer), a comprehensive image analysis software implemented as a userfriendly web-based application with modules that can be put together by the user as pipelines conveniently through a graphical user interface. SPEX’s infrastructure allows for streamlined access to open source image data management systems and analysis modules for cell segmentation, cell phenotyping, cell-cell co-occurrence and spatially informed omics analyses. We demonstrate SPEX’s ability to facilitate the discovery of biological insights in spatially resolved omics datasets from healthy tissue to tumor samples.

Background. DCIS consists of a molecularly heterogeneous group of premalignant lesions, with variable risk of invasive progression. Understanding biomarkers for invasive progression could help individualize treatment recommendations based upon tumor biology. As part of the NCI Human Tumor Atlas Network (HTAN), we conducted comprehensive genomic analyses on two large DCIS case-control cohorts. Methods. We performed smart3-seq and low-pass whole genome sequencing on two independent, retrospective, longitudinally sampled DCIS case-control cohorts. TBCRC 038 was a multicenter cohort diagnosed with DCIS between 1998 and 2016 at one of the Translational Breast Cancer Research sites; the RAHBT (Resource of Archival Human Breast Tissue) cohort included women identified through the St. Louis Breast Tissue Repository, and the Women’s Health Repository diagnosed between 1997 and 2001. We studied the spectrum of molecular changes present and sought genomic predictors of subsequent ipsilateral breast events (iBEs: DCIS recurrence or invasive progression) in both DCIS epithelium and stroma in formalin fixed paraffin embedded tissue. We generated de novo tumor and stroma-centric subtypes for DCIS that represents fundamental transcriptomic organization. Copy number analysis was performed using low-pass DNA sequencing. Non-negative matrix factorization (NMF) was applied to the RNA expression of all coding genes to identify clusters. A negative-binomial regression model was used to identify differentially expressed genes. Results. We analyzed 677 DCIS samples from 481 patients with 7.1 years median follow-up. In TBCRC samples, we identified three clusters via NMF in TBCRC referred to as ER low, quiescent, and ER high. The ER-low cluster had significantly higher levels of ERBB2 and lower levels of ESR1 compared to quiescent and ER-high clusters. Quiescent cluster lesions were less proliferative and less metabolically active than ER high and ER low subtypes. These findings were replicated in the RAHBT cohort. Focusing on the stromal component of DCIS from laser capture microdissection in RAHBT samples, we identified four distinct DCIS-associated stromal clusters. A “normal-like” stromal cluster with ECM organization and PI3K-AKT signaling; a “collagen-rich” stromal cluster; a “desmoplastic” stromal cluster with high fibroblast and total myeloid abundance, mostly associated with macrophages and myeloid dendritic cells (mDC); and an “immune-dense” stromal cluster. Further, we compared differentially expressed genes in patients with or without subsequent iBEs within 5 years of diagnosis. Hypothesizing that the resulting 812 DE genes (DESeq2) represent multiple routes to subsequent iBEs, we leveraged NMF to identify paths to progression. In both TBCRC and RAHBT cohorts, poor outcome groups exhibited increased ER, MYC signaling, and oxidative phosphorylation, supporting that these pathways are important for DCIS recurrence and progression. Conclusion. Comprehensive genomic profiling in two independent DCIS cohorts with longitudinal outcomes shows distinct DCIS stromal expression patterns and immune cell composition. RNA expression profiles reveal underlying tumor biology that is associated with later iBEs in both cohorts. These studies provide new insight into DCIS biology and will guide the design of diagnostic strategies to prevent invasive progression. Citation Format: Siri H Strand, Belén Rivero-Gutiérrez, Kathleen E Houlahan, Jose A Seoane, Lorraine M King, Tyler Risom, Lunden Simpson, Sujay Vennam, Aziz Khan, Timothy Hardman, Bryan E Harmon, Fergus J Couch, Kristalyn Gallagher, Mark Kilgore, Shi Wei, Angela DeMichele, Tari King, Priscilla F McAuliffe, Julie Nangia, Joanna Lee, Jennifer Tseng, Anna Maria Storniolo, Alastair Thompson, Gaorav Gupta, Robyn Burns, Deborah J Veis, Katherine DeSchryver, Chunfang Zhu, Magdalena Matusiak, Jason Wang, Shirley X Zhu, Jen Tappenden, Daisy Yi Ding, Dadong Zhang, Jingqin Luo, Shu Jiang, Sushama Varma, Cody Straub, Sucheta Srivastava, Christina Curtis, Rob Tibshirani, Robert Michael Angelo, Allison Hall, Kouros Owzar, Kornelia Polyak, Carlo Maley, Jeffrey R Marks, Graham A Colditz, E Shelley Hwang, Robert B West. The Breast PreCancer Atlas DCIS genomic signatures define biology and correlate with clinical outcomes: An analysis of TBCRC 038 and RAHBT cohorts [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr GS4-07.

Ductal carcinoma in situ (DCIS) is the most common precursor of invasive breast cancer (IBC), with variable propensity for progression. We perform multiscale, integrated molecular profiling of DCIS with clinical outcomes by analyzing 774 DCIS samples from 542 patients with 7.3 years median follow-up from the Translational Breast Cancer Research Consortium 038 study and the Resource of Archival Breast Tissue cohorts. We identify 812 genes associated with ipsilateral recurrence within 5 years from treatment and develop a classifier that predicts DCIS or IBC recurrence in both cohorts. Pathways associated with recurrence include proliferation, immune response, and metabolism. Distinct stromal expression patterns and immune cell compositions are identified. Our multiscale approach employed in situ methods to generate a spatially resolved atlas of breast precancers, where complementary modalities can be directly compared and correlated with conventional pathology findings, disease states, and clinical outcome.

Next-generation tools for multiplexed imaging have driven a new wave of innovation in understanding how single-cell function and tissue structure are interrelated. In previous work, we developed multiplexed ion beam imaging by time of flight, a highly multiplexed platform that uses secondary ion mass spectrometry to image dozens of antibodies tagged with metal reporters. As instrument throughput has increased, the breadth and depth of imaging data have increased as well. To extract meaningful information from these data, we have developed tools for cell identification, cell classification, and spatial analysis. In this review, we discuss these tools and provide examples of their application in various contexts, including ductal carcinoma in situ, tuberculosis, and Alzheimer's disease. We hope the synergy between multiplexed imaging and automated image analysis will drive a new era in anatomic pathology and personalized medicine wherein quantitative spatial signatures are used routinely for more accurate diagnosis, prognosis, and therapeutic selection.