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.

A Cancer Cell Program Promotes T Cell Exclusion and Resistance to Checkpoint Blockade.

Immune checkpoint inhibitors (ICIs) produce durable responses in some melanoma patients, but many patients derive no clinical benefit, and the molecular underpinnings of such resistance remain elusive. Here, we leveraged single-cell RNA sequencing (scRNA-seq) from 33 melanoma tumors and computational analyses to interrogate malignant cell states that promote immune evasion. We identified a resistance program expressed by malignant cells that is associated with T cell exclusion and immune evasion. The program is expressed prior to immunotherapy, characterizes cold niches in situ, and predicts clinical responses to anti-PD-1 therapy in an independent cohort of 112 melanoma patients. CDK4/6-inhibition represses this program in individual malignant cells, induces senescence, and reduces melanoma tumor outgrowth in mouse models in vivo when given in combination with immunotherapy. Our study provides a high-resolution landscape of ICI-resistant cell states, identifies clinically predictive signatures, and suggests new therapeutic strategies to overcome immunotherapy resistance.

MCMICRO: a scalable, modular image-processing pipeline for multiplexed tissue imaging.

Highly multiplexed tissue imaging makes detailed molecular analysis of single cells possible in a preserved spatial context. However, reproducible analysis of large multichannel images poses a substantial computational challenge. Here, we describe a modular and open-source computational pipeline, MCMICRO, for performing the sequential steps needed to transform whole-slide images into single-cell data. We demonstrate the use of MCMICRO on tissue and tumor images acquired using multiple imaging platforms, thereby providing a solid foundation for the continued development of tissue imaging software.

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.

Fibroblast-tumor cell signaling limits HER2 kinase therapy response via activation of MTOR and antiapoptotic pathways.

Despite the implementation of multiple HER2-targeted therapies, patients with advanced HER2+ breast cancer ultimately develop drug resistance. Stromal fibroblasts represent an abundant cell type in the tumor microenvironment and have been linked to poor outcomes and drug resistance. Here, we show that fibroblasts counteract the cytotoxic effects of HER2 kinase-targeted therapy in a subset of HER2+ breast cancer cell lines and allow cancer cells to proliferate in the presence of the HER2 kinase inhibitor lapatinib. Fibroblasts from primary breast tumors, normal breast tissue, and lung tissue have similar protective effects on tumor cells via paracrine factors. This fibroblast-mediated reduction in drug sensitivity involves increased expression of antiapoptotic proteins and sustained activation of the PI3K/AKT/MTOR pathway, despite inhibition of the HER2 and the RAS-ERK pathways in tumor cells. HER2 therapy sensitivity is restored in the fibroblast cocultures by combination treatment with inhibitors of MTOR or the antiapoptotic proteins BCL-XL and MCL-1. Expression of activated AKT in tumor cells recapitulates the effects of fibroblasts resulting in sustained MTOR signaling and poor lapatinib response. Lapatinib sensitivity was not altered by fibroblasts in tumor cells that exhibited sustained MTOR signaling due to a strong gain-of-function PI3KCA mutation. These findings indicate that in addition to tumor cell-intrinsic mechanisms that cause constitutive PI3K/AKT/MTOR pathway activation, secreted factors from fibroblasts can maintain this pathway in the context of HER2 inhibition. Our integrated proteomic-phenotypic approach presents a strategy for the discovery of protective mechanisms in fibroblast-rich tumors and the design of rational combination therapies to restore drug sensitivity.

Dosage of Dyrk1a shifts cells within a p21-cyclin D1 signaling map to control the decision to enter the cell cycle.

Mammalian cells have a remarkable capacity to compensate for heterozygous gene loss or extra gene copies. One exception is Down syndrome (DS), where a third copy of chromosome 21 mediates neurogenesis defects and lowers the frequency of solid tumors. Here we combine live-cell imaging and single-cell analysis to show that increased dosage of chromosome 21-localized Dyrk1a steeply increases G1 cell cycle duration through direct phosphorylation and degradation of cyclin D1 (CycD1). DS-derived fibroblasts showed analogous cell cycle changes that were reversed by Dyrk1a inhibition. Furthermore, reducing Dyrk1a activity increased CycD1 expression to force a bifurcation, with one subpopulation of cells accelerating proliferation and the other arresting proliferation by costabilizing CycD1 and the CDK inhibitor p21. Thus, dosage of Dyrk1a repositions cells within a p21-CycD1 signaling map, directing each cell to either proliferate or to follow two distinct cell cycle exit pathways characterized by high or low CycD1 and p21 levels.

NEK8 links the ATR-regulated replication stress response and S phase CDK activity to renal ciliopathies.

Renal ciliopathies are a leading cause of kidney failure, but their exact etiology is poorly understood. NEK8/NPHP9 is a ciliary kinase associated with two renal ciliopathies in humans and mice, nephronophthisis (NPHP) and polycystic kidney disease. Here, we identify NEK8 as a key effector of the ATR-mediated replication stress response. Cells lacking NEK8 form spontaneous DNA double-strand breaks (DSBs) that further accumulate when replication forks stall, and they exhibit reduced fork rates, unscheduled origin firing, and increased replication fork collapse. NEK8 suppresses DSB formation by limiting cyclin A-associated CDK activity. Strikingly, a mutation in NEK8 that is associated with renal ciliopathies affects its genome maintenance functions. Moreover, kidneys of NEK8 mutant mice accumulate DNA damage, and loss of NEK8 or replication stress similarly disrupts renal cell architecture in a 3D-culture system. Thus, NEK8 is a critical component of the DNA damage response that links replication stress with cystic kidney disorders.

Inferring reaction network structure from single-cell, multiplex data, using toric systems theory.

The goal of many single-cell studies on eukaryotic cells is to gain insight into the biochemical reactions that control cell fate and state. In this paper we introduce the concept of Effective Stoichiometric Spaces (ESS) to guide the reconstruction of biochemical networks from multiplexed, fixed time-point, single-cell data. In contrast to methods based solely on statistical models of data, the ESS method leverages the power of the geometric theory of toric varieties to begin unraveling the structure of chemical reaction networks (CRN). This application of toric theory enables a data-driven mapping of covariance relationships in single-cell measurements into stoichiometric information, one in which each cell subpopulation has its associated ESS interpreted in terms of CRN theory. In the development of ESS we reframe certain aspects of the theory of CRN to better match data analysis. As an application of our approach we process cytomery- and image-based single-cell datasets and identify differences in cells treated with kinase inhibitors. Our approach is directly applicable to data acquired using readily accessible experimental methods such as Fluorescence Activated Cell Sorting (FACS) and multiplex immunofluorescence.

A two-dimensional ERK-AKT signaling code for an NGF-triggered cell-fate decision.

Growth factors activate Ras, PI3K, and other signaling pathways. It is not well understood how these signals are translated by individual cells into a decision to proliferate or differentiate. Here, using single-cell image analysis of nerve growth factor (NGF)-stimulated PC12 cells, we identified a two-dimensional phospho-ERK (pERK)-phospho-AKT (pAKT) response map with a curved boundary that separates differentiating from proliferating cells. The boundary position remained invariant when different stimuli were used or upstream signaling components perturbed. We further identified Rasa2 as a negative feedback regulator that links PI3K to Ras, placing the stochastically distributed pERK-pAKT signals close to the decision boundary. This allows for uniform NGF stimuli to create a subpopulation of cells that differentiates with each cycle of proliferation. Thus, by linking a complex signaling system to a simpler intermediate response map, cells gain unique integration and control capabilities to balance cell number expansion with differentiation.

Determinants of drug-target interactions at the single cell level.

The physiochemical determinants of drug-target interactions in the microenvironment of the cell are complex and generally not defined by simple diffusion and intrinsic chemical reactivity. Non-specific interactions of drugs and macromolecules in cells are rarely considered formally in assessing pharmacodynamics. Here, we demonstrate that non-specific interactions lead to very slow incorporation kinetics of DNA binding drugs. We observe a rate of drug incorporation in cell nuclei three orders of magnitude slower than in vitro due to anomalous drug diffusion within cells. This slow diffusion, however, has an advantageous consequence: it leads to virtually irreversible binding of the drug to specific DNA targets in cells. We show that non-specific interactions drive slow drug diffusion manifesting as slow reaction front propagation. We study the effect of non-specific interactions in different cellular compartments by permeabilization of plasma and nuclear membranes in order to pinpoint differential compartment effects on variability in intracellular drug kinetics. These results provide the basis for a comprehensive model of the determinants of intracellular diffusion of small-molecule drugs, their target-seeking trajectories, and the consequences of these processes on the apparent kinetics of drug-target interactions.

SHPRH and HLTF act in a damage-specific manner to coordinate different forms of postreplication repair and prevent mutagenesis.

Postreplication repair (PRR) pathways play important roles in restarting stalled replication forks and regulating mutagenesis. In yeast, Rad5-mediated damage avoidance and Rad18-mediated translesion synthesis (TLS) are two forms of PRR. Two Rad5-related proteins, SHPRH and HLTF, have been identified in mammalian cells, but their specific roles in PRR are unclear. Here, we show that HLTF and SHPRH suppress mutagenesis in a damage-specific manner, preventing mutations induced by UV and MMS, respectively. Following UV, HLTF enhances PCNA monoubiquitination and recruitment of TLS polymerase η, while also inhibiting SHPRH function. In contrast, MMS promotes the degradation of HLTF and the interactions of SHPRH with Rad18 and polymerase κ. Our data suggest not only that cells differentially utilize HLTF and SHPRH for different forms of DNA damage, but also, surprisingly, that HLTF and SHPRH may coordinate the two main branches of PRR to choose the proper bypass mechanism for minimizing mutagenesis.

Adaptive resistance of melanoma cells to RAF inhibition via reversible induction of a slowly dividing de-differentiated state.

Treatment of BRAF-mutant melanomas with MAP kinase pathway inhibitors is paradigmatic of the promise of precision cancer therapy but also highlights problems with drug resistance that limit patient benefit. We use live-cell imaging, single-cell analysis, and molecular profiling to show that exposure of tumor cells to RAF/MEK inhibitors elicits a heterogeneous response in which some cells die, some arrest, and the remainder adapt to drug. Drug-adapted cells up-regulate markers of the neural crest (e.g., NGFR), a melanocyte precursor, and grow slowly. This phenotype is transiently stable, reverting to the drug-naïve state within 9 days of drug withdrawal. Transcriptional profiling of cell lines and human tumors implicates a c-Jun/ECM/FAK/Src cascade in de-differentiation in about one-third of cell lines studied; drug-induced changes in c-Jun and NGFR levels are also observed in xenograft and human tumors. Drugs targeting the c-Jun/ECM/FAK/Src cascade as well as BET bromodomain inhibitors increase the maximum effect (Emax) of RAF/MEK kinase inhibitors by promoting cell killing. Thus, analysis of reversible drug resistance at a single-cell level identifies signaling pathways and inhibitory drugs missed by assays that focus on cell populations.

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.

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.

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.

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.

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.

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.

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.