Genomic and Immunophenotypic Landscape of Acquired Resistance to PD-(L)1 Blockade in Non-Small-Cell Lung Cancer.

PURPOSE: Although immune checkpoint inhibitors (ICI) have extended survival in patients with non-small-cell lung cancer (NSCLC), acquired resistance (AR) to ICI frequently develops after an initial benefit. However, the mechanisms of AR to ICI in NSCLC are largely unknown. METHODS: Comprehensive tumor genomic profiling, machine learning-based assessment of tumor-infiltrating lymphocytes, multiplexed immunofluorescence, and/or HLA-I immunohistochemistry (IHC) were performed on matched pre- and post-ICI tumor biopsies from patients with NSCLC treated with ICI at the Dana-Farber Cancer Institute who developed AR to ICI. Two additional cohorts of patients with intervening chemotherapy or targeted therapies between biopsies were included as controls. RESULTS: We performed comprehensive genomic profiling and immunophenotypic characterization on samples from 82 patients with NSCLC and matched pre- and post-ICI biopsies and compared findings with a control cohort of patients with non-ICI intervening therapies between biopsies (chemotherapy, N = 32; targeted therapies, N = 89; both, N = 17). Putative resistance mutations were identified in 27.8% of immunotherapy-treated cases and included acquired loss-of-function mutations in STK11, B2M, APC, MTOR, KEAP1, and JAK1/2; these acquired alterations were not observed in the control groups. Immunophenotyping of matched pre- and post-ICI samples demonstrated significant decreases in intratumoral lymphocytes, CD3e+ and CD8a+ T cells, and PD-L1-PD1 engagement, as well as increased distance between tumor cells and CD8+PD-1+ T cells. There was a significant decrease in HLA class I expression in the immunotherapy cohort at the time of AR compared with the chemotherapy (P = .005) and the targeted therapy (P = .01) cohorts. CONCLUSION: These findings highlight the genomic and immunophenotypic heterogeneity of ICI resistance in NSCLC, which will need to be considered when developing novel therapeutic strategies aimed at overcoming resistance.

Qualification of a multiplexed tissue imaging assay and detection of novel patterns of HER2 heterogeneity in breast cancer.

Emerging data suggests that HER2 intratumoral heterogeneity (ITH) is associated with therapy resistance, highlighting the need for new strategies to assess HER2 ITH. A promising approach is leveraging multiplexed tissue analysis techniques such as cyclic immunofluorescence (CyCIF), which enable visualization and quantification of 10-60 antigens at single-cell resolution from individual tissue sections. In this study, we qualified a breast cancer-specific antibody panel, including HER2, ER, and PR, for multiplexed tissue imaging. We then compared the performance of these antibodies against established clinical standards using pixel-, cell- and tissue-level analyses, utilizing 866 tissue cores (representing 294 patients). To ensure reliability, the CyCIF antibodies were qualified against HER2 immunohistochemistry (IHC) and fluorescence in situ hybridization (FISH) data from the same samples. Our findings demonstrate the successful qualification of a breast cancer antibody panel for CyCIF, showing high concordance with established clinical antibodies. Subsequently, we employed the qualified antibodies, along with antibodies for CD45, CD68, PD-L1, p53, Ki67, pRB, and AR, to characterize 567 HER2+ invasive breast cancer samples from 189 patients. Through single-cell analysis, we identified four distinct cell clusters within HER2+ breast cancer exhibiting heterogeneous HER2 expression. Furthermore, these clusters displayed variations in ER, PR, p53, AR, and PD-L1 expression. To quantify the extent of heterogeneity, we calculated heterogeneity scores based on the diversity among these clusters. Our analysis revealed expression patterns that are relevant to breast cancer biology, with correlations to HER2 ITH and potential relevance to clinical outcomes.

Phase II Study of Eribulin plus Pembrolizumab in Metastatic Soft-tissue Sarcomas: Clinical Outcomes and Biological Correlates.

PURPOSE: Eribulin modulates the tumor-immune microenvironment via cGAS-STING signaling in preclinical models. This non-randomized phase II trial evaluated the combination of eribulin and pembrolizumab in patients with soft-tissue sarcomas (STS). PATIENTS AND METHODS: Patients enrolled in one of three cohorts: leiomyosarcoma (LMS), liposarcomas (LPS), or other STS that may benefit from PD-1 inhibitors, including undifferentiated pleomorphic sarcoma (UPS). Eribulin was administered at 1.4 mg/m2 i.v. (days 1 and 8) with fixed-dose pembrolizumab 200 mg i.v. (day 1) of each 21-day cycle, until progression, unacceptable toxicity, or completion of 2 years of treatment. The primary endpoint was the 12-week progression-free survival rate (PFS-12) in each cohort. Secondary endpoints included the objective response rate, median PFS, safety profile, and overall survival (OS). Pretreatment and on-treatment blood specimens were evaluated in patients who achieved durable disease control (DDC) or progression within 12 weeks [early progression (EP)]. Multiplexed immunofluorescence was performed on archival LPS samples from patients with DDC or EP. RESULTS: Fifty-seven patients enrolled (LMS, n = 19; LPS, n = 20; UPS/Other, n = 18). The PFS-12 was 36.8% (90% confidence interval: 22.5-60.4) for LMS, 69.6% (54.5-89.0) for LPS, and 52.6% (36.8-75.3) for UPS/Other cohorts. All 3 patients in the UPS/Other cohort with angiosarcoma achieved RECIST responses. Toxicity was manageable. Higher IFNα and IL4 serum levels were associated with clinical benefit. Immune aggregates expressing PD-1 and PD-L1 were observed in a patient that completed 2 years of treatment. CONCLUSIONS: The combination of eribulin and pembrolizumab demonstrated promising activity in LPS and angiosarcoma.

Quality Control for Single Cell Analysis of High-plex Tissue Profiles using CyLinter.

Tumors are complex assemblies of cellular and acellular structures patterned on spatial scales from microns to centimeters. Study of these assemblies has advanced dramatically with the introduction of high-plex spatial profiling. Image-based profiling methods reveal the intensities and spatial distributions of 20-100 proteins at subcellular resolution in 103-107 cells per specimen. Despite extensive work on methods for extracting single-cell data from these images, all tissue images contain artefacts such as folds, debris, antibody aggregates, optical aberrations and image processing errors that arise from imperfections in specimen preparation, data acquisition, image assembly, and feature extraction. We show that these artefacts dramatically impact single-cell data analysis, obscuring meaningful biological interpretation. We describe an interactive quality control software tool, CyLinter, that identifies and removes data associated with imaging artefacts. CyLinter greatly improves single-cell analysis, especially for archival specimens sectioned many years prior to data collection, such as those from clinical trials.

3D multiplexed tissue imaging reconstruction and optimized region of interest (ROI) selection through deep learning model of channels embedding.

Introduction: Tissue-based sampling and diagnosis are defined as the extraction of information from certain limited spaces and its diagnostic significance of a certain object. Pathologists deal with issues related to tumor heterogeneity since analyzing a single sample does not necessarily capture a representative depiction of cancer, and a tissue biopsy usually only presents a small fraction of the tumor. Many multiplex tissue imaging platforms (MTIs) make the assumption that tissue microarrays (TMAs) containing small core samples of 2-dimensional (2D) tissue sections are a good approximation of bulk tumors although tumors are not 2D. However, emerging whole slide imaging (WSI) or 3D tumor atlases that use MTIs like cyclic immunofluorescence (CyCIF) strongly challenge this assumption. In spite of the additional insight gathered by measuring the tumor microenvironment in WSI or 3D, it can be prohibitively expensive and time-consuming to process tens or hundreds of tissue sections with CyCIF. Even when resources are not limited, the criteria for region of interest (ROI) selection in tissues for downstream analysis remain largely qualitative and subjective as stratified sampling requires the knowledge of objects and evaluates their features. Despite the fact TMAs fail to adequately approximate whole tissue features, a theoretical subsampling of tissue exists that can best represent the tumor in the whole slide image. Methods: To address these challenges, we propose deep learning approaches to learn multi-modal image translation tasks from two aspects: 1) generative modeling approach to reconstruct 3D CyCIF representation and 2) co-embedding CyCIF image and Hematoxylin and Eosin (H&E) section to learn multi-modal mappings by a cross-domain translation for minimum representative ROI selection. Results and discussion: We demonstrate that generative modeling enables a 3D virtual CyCIF reconstruction of a colorectal cancer specimen given a small subset of the imaging data at training time. By co-embedding histology and MTI features, we propose a simple convex optimization for objective ROI selection. We demonstrate the potential application of ROI selection and the efficiency of its performance with respect to cellular heterogeneity.

Endocrine Therapy Synergizes with SMAC Mimetics to Potentiate Antigen Presentation and Tumor Regression in Hormone Receptor-Positive Breast Cancer.

Immunotherapies have yet to demonstrate significant efficacy in the treatment of hormone receptor-positive (HR+) breast cancer. Given that endocrine therapy (ET) is the primary approach for treating HR+ breast cancer, we investigated the effects of ET on the tumor immune microenvironment (TME) in HR+ breast cancer. Spatial proteomics of primary HR+ breast cancer samples obtained at baseline and after ET from patients enrolled in a neoadjuvant clinical trial (NCT02764541) indicated that ET upregulated ß2-microglobulin and influenced the TME in a manner that promotes enhanced immunogenicity. To gain a deeper understanding of the underlying mechanisms, the intrinsic effects of ET on cancer cells were explored, which revealed that ET plays a crucial role in facilitating the chromatin binding of RelA, a key component of the NF-κB complex. Consequently, heightened NF-κB signaling enhanced the response to interferon-gamma, leading to the upregulation of ß2-microglobulin and other antigen presentation-related genes. Further, modulation of NF-κB signaling using a SMAC mimetic in conjunction with ET augmented T-cell migration and enhanced MHC-I-specific T-cell-mediated cytotoxicity. Remarkably, the combination of ET and SMAC mimetics, which also blocks prosurvival effects of NF-κB signaling through the degradation of inhibitors of apoptosis proteins, elicited tumor regression through cell autonomous mechanisms, providing additional support for their combined use in HR+ breast cancer. SIGNIFICANCE: Adding SMAC mimetics to endocrine therapy enhances tumor regression in a cell autonomous manner while increasing tumor immunogenicity, indicating that this combination could be an effective treatment for HR+ patients with breast cancer.

Sequential apoptotic and multiplexed proteomic evaluation of single cancer cells.

A potential cause of cancer relapse is pretreatment chemoresistant subpopulations. Identifying targetable features of subpopulations that are poorly primed for therapy-induced cell death may improve cancer therapy. Here, we develop and validate real-time BH3 profiling, a live and functional single-cell measurement of pretreatment apoptotic sensitivity that occurs upstream of apoptotic protease activation. On the same single cells, we perform cyclic immunofluorescence, which enables multiplexed immunofluorescence of more than 30 proteins on the same cell. Using cultured cells and rapid ex vivo cultures of colon cancer patient-derived xenograft (PDX) models, we identify Bak as a univariate correlate of apoptotic priming, find that poorly primed subpopulations can correspond to specific stages of the cell cycle, and, in some PDX models, identify increased expression of Bcl-XL, Mcl-1, or Her2 in subpopulations that are poorly primed for apoptosis. Last, we generate and validate mathematical models of single-cell priming that describe how targetable proteins contribute to apoptotic priming.

High-plex immunofluorescence imaging and traditional histology of the same tissue section for discovering image-based biomarkers.

Precision medicine is critically dependent on better methods for diagnosing and staging disease and predicting drug response. Histopathology using hematoxylin and eosin (H&E)-stained tissue (not genomics) remains the primary diagnostic method in cancer. Recently developed highly multiplexed tissue imaging methods promise to enhance research studies and clinical practice with precise, spatially resolved single-cell data. Here, we describe the 'Orion' platform for collecting H&E and high-plex immunofluorescence images from the same cells in a whole-slide format suitable for diagnosis. Using a retrospective cohort of 74 colorectal cancer resections, we show that immunofluorescence and H&E images provide human experts and machine learning algorithms with complementary information that can be used to generate interpretable, multiplexed image-based models predictive of progression-free survival. Combining models of immune infiltration and tumor-intrinsic features achieves a 10- to 20-fold discrimination between rapid and slow (or no) progression, demonstrating the ability of multimodal tissue imaging to generate high-performance biomarkers.

Immune Profiling of Dermatologic Adverse Events from Checkpoint Blockade using Tissue Cyclic Immunofluorescence.

In this study, we demonstrate the utility of whole-slide CyCIF (tissue-based cyclic immunofluorescence) imaging for characterizing immune cell infiltrates in immune checkpoint inhibitor (ICI)-induced dermatologic adverse events (dAEs). We analyzed six cases of ICI-induced dAEs, including lichenoid, bullous pemphigoid, psoriasis, and eczematous eruptions, comparing immune profiling results obtained using both standard immunohistochemistry (IHC) and CyCIF. Our findings indicate that CyCIF provides more detailed and precise single-cell characterization of immune cell infiltrates than IHC, which relies on semi-quantitative scoring by pathologists. This pilot study highlights the potential of CyCIF to advance our understanding of the immune environment in dAEs by revealing tissue-level spatial patterns of immune cell infiltrates, allowing for more precise phenotypic distinctions and deeper exploration of disease mechanisms. By demonstrating that CyCIF can be performed on friable tissues, such as bullous pemphigoid, we provide a foundation for future studies to examine the drivers of specific dAEs using larger cohorts of phenotyped toxicity and suggest a broader role for highly multiplexed tissue imaging in phenotyping the immune mediated disease that they resemble.

Lymphocyte networks are dynamic cellular communities in the immunoregulatory landscape of lung adenocarcinoma.

Lymphocytes are key for immune surveillance of tumors, but our understanding of the spatial organization and physical interactions that facilitate lymphocyte anti-cancer functions is limited. We used multiplexed imaging, quantitative spatial analysis, and machine learning to create high-definition maps of lung tumors from a Kras/Trp53-mutant mouse model and human resections. Networks of interacting lymphocytes ("lymphonets") emerged as a distinctive feature of the anti-cancer immune response. Lymphonets nucleated from small T cell clusters and incorporated B cells with increasing size. CXCR3-mediated trafficking modulated lymphonet size and number, but T cell antigen expression directed intratumoral localization. Lymphonets preferentially harbored TCF1+ PD-1+ progenitor CD8+ T cells involved in responses to immune checkpoint blockade (ICB) therapy. Upon treatment of mice with ICB or an antigen-targeted vaccine, lymphonets retained progenitor and gained cytotoxic CD8+ T cell populations, likely via progenitor differentiation. These data show that lymphonets create a spatial environment supportive of CD8+ T cell anti-tumor responses.

Multiplexed 3D atlas of state transitions and immune interaction in colorectal cancer.

Advanced solid cancers are complex assemblies of tumor, immune, and stromal cells characterized by high intratumoral variation. We use highly multiplexed tissue imaging, 3D reconstruction, spatial statistics, and machine learning to identify cell types and states underlying morphological features of known diagnostic and prognostic significance in colorectal cancer. Quantitation of these features in high-plex marker space reveals recurrent transitions from one tumor morphology to the next, some of which are coincident with long-range gradients in the expression of oncogenes and epigenetic regulators. At the tumor invasive margin, where tumor, normal, and immune cells compete, T cell suppression involves multiple cell types and 3D imaging shows that seemingly localized 2D features such as tertiary lymphoid structures are commonly interconnected and have graded molecular properties. Thus, while cancer genetics emphasizes the importance of discrete changes in tumor state, whole-specimen imaging reveals large-scale morphological and molecular gradients analogous to those in developing tissues.

Computational multiplex panel reduction to maximize information retention in breast cancer tissue microarrays.

Recent state-of-the-art multiplex imaging techniques have expanded the depth of information that can be captured within a single tissue sample by allowing for panels with dozens of markers. Despite this increase in capacity, space on the panel is still limited due to technical artifacts, tissue loss, and long imaging acquisition time. As such, selecting which markers to include on a panel is important, since removing important markers will result in a loss of biologically relevant information, but identifying redundant markers will provide a room for other markers. To address this, we propose computational approaches to determine the amount of shared information between markers and select an optimally reduced panel that captures maximum amount of information with the fewest markers. Here we examine several panel selection approaches and evaluate them based on their ability to reconstruct the full panel images and information within breast cancer tissue microarray datasets using cyclic immunofluorescence as a proof of concept. We show that all methods perform adequately and can re-capture cell types using only 18 of 25 markers (72% of the original panel size). The correlation-based selection methods achieved the best single-cell marker mean intensity predictions with a Spearman correlation of 0.90 with the reduced panel. Using the proposed methods shown here, it is possible for researchers to design more efficient multiplex imaging panels that maximize the amount of information retained with the limited number of markers with respect to certain evaluation metrics and architecture biases.

Oncometabolite d-2HG alters T cell metabolism to impair CD8+ T cell function.

Gain-of-function mutations in isocitrate dehydrogenase (IDH) in human cancers result in the production of d-2-hydroxyglutarate (d-2HG), an oncometabolite that promotes tumorigenesis through epigenetic alterations. The cancer cell-intrinsic effects of d-2HG are well understood, but its tumor cell-nonautonomous roles remain poorly explored. We compared the oncometabolite d-2HG with its enantiomer, l-2HG, and found that tumor-derived d-2HG was taken up by CD8+ T cells and altered their metabolism and antitumor functions in an acute and reversible fashion. We identified the glycolytic enzyme lactate dehydrogenase (LDH) as a molecular target of d-2HG. d-2HG and inhibition of LDH drive a metabolic program and immune CD8+ T cell signature marked by decreased cytotoxicity and impaired interferon-γ signaling that was recapitulated in clinical samples from human patients with IDH1 mutant gliomas.

Single cell spatial analysis reveals the topology of immunomodulatory purinergic signaling in glioblastoma.

How the glioma immune microenvironment fosters tumorigenesis remains incompletely defined. Here, we use single-cell RNA-sequencing and multiplexed tissue-imaging to characterize the composition, spatial organization, and clinical significance of extracellular purinergic signaling in glioma. We show that microglia are the predominant source of CD39, while tumor cells principally express CD73. In glioblastoma, CD73 is associated with EGFR amplification, astrocyte-like differentiation, and increased adenosine, and is linked to hypoxia. Glioblastomas enriched for CD73 exhibit inflammatory microenvironments, suggesting that purinergic signaling regulates immune adaptation. Spatially-resolved single-cell analyses demonstrate a strong spatial correlation between tumor-CD73 and microglial-CD39, with proximity associated with poor outcomes. Similar spatial organization is present in pediatric high-grade gliomas including H3K27M-mutant diffuse midline glioma. These data reveal that purinergic signaling in gliomas is shaped by genotype, lineage, and functional state, and that core enzymes expressed by tumor and myeloid cells are organized to promote adenosine-rich microenvironments potentially amenable to therapeutic targeting.

A human breast atlas integrating single-cell proteomics and transcriptomics.

The breast is a dynamic organ whose response to physiological and pathophysiological conditions alters its disease susceptibility, yet the specific effects of these clinical variables on cell state remain poorly annotated. We present a unified, high-resolution breast atlas by integrating single-cell RNA-seq, mass cytometry, and cyclic immunofluorescence, encompassing a myriad of states. We define cell subtypes within the alveolar, hormone-sensing, and basal epithelial lineages, delineating associations of several subtypes with cancer risk factors, including age, parity, and BRCA2 germline mutation. Of particular interest is a subset of alveolar cells termed basal-luminal (BL) cells, which exhibit poor transcriptional lineage fidelity, accumulate with age, and carry a gene signature associated with basal-like breast cancer. We further utilize a medium-depletion approach to identify molecular factors regulating cell-subtype proportion in organoids. Together, these data are a rich resource to elucidate diverse mammary cell states.

Temporal and spatial topography of cell proliferation in cancer.

Proliferation is a fundamental trait of cancer cells, but its properties and spatial organization in tumours are poorly characterized. Here we use highly multiplexed tissue imaging to perform single-cell quantification of cell cycle regulators and then develop robust, multivariate, proliferation metrics. Across diverse cancers, proliferative architecture is organized at two spatial scales: large domains, and smaller niches enriched for specific immune lineages. Some tumour cells express cell cycle regulators in the (canonical) patterns expected of freely growing cells, a phenomenon we refer to as 'cell cycle coherence'. By contrast, the cell cycles of other tumour cell populations are skewed towards specific phases or exhibit non-canonical (incoherent) marker combinations. Coherence varies across space, with changes in oncogene activity and therapeutic intervention, and is associated with aggressive tumour behaviour. Thus, multivariate measures from high-plex tissue images capture clinically significant features of cancer proliferation, a fundamental step in enabling more precise use of anti-cancer therapies.

Scope2Screen: Focus+Context Techniques for Pathology Tumor Assessment in Multivariate Image Data.

Inspection of tissues using a light microscope is the primary method of diagnosing many diseases, notably cancer. Highly multiplexed tissue imaging builds on this foundation, enabling the collection of up to 60 channels of molecular information plus cell and tissue morphology using antibody staining. This provides unique insight into disease biology and promises to help with the design of patient-specific therapies. However, a substantial gap remains with respect to visualizing the resulting multivariate image data and effectively supporting pathology workflows in digital environments on screen. We, therefore, developed Scope2Screen, a scalable software system for focus+context exploration and annotation of whole-slide, high-plex, tissue images. Our approach scales to analyzing 100GB images of 109 or more pixels per channel, containing millions of individual cells. A multidisciplinary team of visualization experts, microscopists, and pathologists identified key image exploration and annotation tasks involving finding, magnifying, quantifying, and organizing regions of interest (ROIs) in an intuitive and cohesive manner. Building on a scope-to-screen metaphor, we present interactive lensing techniques that operate at single-cell and tissue levels. Lenses are equipped with task-specific functionality and descriptive statistics, making it possible to analyze image features, cell types, and spatial arrangements (neighborhoods) across image channels and scales. A fast sliding-window search guides users to regions similar to those under the lens; these regions can be analyzed and considered either separately or as part of a larger image collection. A novel snapshot method enables linked lens configurations and image statistics to be saved, restored, and shared with these regions. We validate our designs with domain experts and apply Scope2Screen in two case studies involving lung and colorectal cancers to discover cancer-relevant image features.

Narrative online guides for the interpretation of digital-pathology images and tissue-atlas data.

Multiplexed tissue imaging facilitates the diagnosis and understanding of complex disease traits. However, the analysis of such digital images heavily relies on the experience of anatomical pathologists for the review, annotation and description of tissue features. In addition, the wider use of data from tissue atlases in basic and translational research and in classrooms would benefit from software that facilitates the easy visualization and sharing of the images and the results of their analyses. In this Perspective, we describe the ecosystem of software available for the analysis of tissue images and discuss the need for interactive online guides that help histopathologists make complex images comprehensible to non-specialists. We illustrate this idea via a software interface (Minerva), accessible via web browsers, that integrates multi-omic and tissue-atlas features. We argue that such interactive narrative guides can effectively disseminate digital histology data and aid their interpretation.