Uncovering therapeutic targets for macrophage-mediated T cell suppression and PD-L1 therapy sensitization.

Tumor-associated macrophages (TAMs) and other myelomonocytic cells are implicated in regulating responsiveness to immunotherapies, including immune checkpoint inhibitors (ICIs) targeting the PD-1/PD-L1 axis. We have developed an ex vivo high-throughput approach to discover modulators of macrophage-mediated T cell suppression, which can improve clinical outcomes of ICIs. We screened 1,430 Food and Drug Administration (FDA)-approved small-molecule drugs using a co-culture assay employing bone-marrow-derived macrophages (BMDMs) and splenic-derived T cells. This identified 57 compounds that disrupted macrophage-mediated T cell suppression. Seven compounds exerted prominent synergistic T cell expansion activity when combined with αPD-L1. These include four COX1/2 inhibitors and two myeloid cell signaling inhibitors. We demonstrate that the use of cyclooxygenase (COX)1/2 inhibitors in combination with αPD-L1 decreases tumor growth kinetics and enhances overall survival in triple-negative breast cancer (TNBC) tumor models in a CD8+ T cell-dependent manner. Altogether, we present a rationalized approach for identifying compounds that synergize with ICI to potentially enhance therapeutic outcomes for patients with solid tumors.

Computational methods and biomarker discovery strategies for spatial proteomics: a review in immuno-oncology.

Advancements in imaging technologies have revolutionized our ability to deeply profile pathological tissue architectures, generating large volumes of imaging data with unparalleled spatial resolution. This type of data collection, namely, spatial proteomics, offers invaluable insights into various human diseases. Simultaneously, computational algorithms have evolved to manage the increasing dimensionality of spatial proteomics inherent in this progress. Numerous imaging-based computational frameworks, such as computational pathology, have been proposed for research and clinical applications. However, the development of these fields demands diverse domain expertise, creating barriers to their integration and further application. This review seeks to bridge this divide by presenting a comprehensive guideline. We consolidate prevailing computational methods and outline a roadmap from image processing to data-driven, statistics-informed biomarker discovery. Additionally, we explore future perspectives as the field moves toward interfacing with other quantitative domains, holding significant promise for precision care in immuno-oncology.

PI3Kγ inhibition circumvents inflammation and vascular leak in SARS-CoV-2 and other infections.

Virulent infectious agents such as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and methicillin-resistant Staphylococcus aureus (MRSA) induce tissue damage that recruits neutrophils, monocyte, and macrophages, leading to T cell exhaustion, fibrosis, vascular leak, epithelial cell depletion, and fatal organ damage. Neutrophils, monocytes, and macrophages recruited to pathogen-infected lungs, including SARS-CoV-2-infected lungs, express phosphatidylinositol 3-kinase gamma (PI3Kγ), a signaling protein that coordinates both granulocyte and monocyte trafficking to diseased tissues and immune-suppressive, profibrotic transcription in myeloid cells. PI3Kγ deletion and inhibition with the clinical PI3Kγ inhibitor eganelisib promoted survival in models of infectious diseases, including SARS-CoV-2 and MRSA, by suppressing inflammation, vascular leak, organ damage, and cytokine storm. These results demonstrate essential roles for PI3Kγ in inflammatory lung disease and support the potential use of PI3Kγ inhibitors to suppress inflammation in severe infectious diseases.

Metabolic abnormalities in the bone marrow cells of young offspring born to mothers with obesity.

BACKGROUND/OBJECTIVES: Intrauterine metabolic reprogramming occurs in mothers with obesity during gestation, putting the offspring at high risk of developing obesity and associated metabolic disorders even before birth. We have generated a mouse model of maternal high-fat diet-induced obesity that recapitulates the metabolic changes seen in humans born to women with obesity. METHODS: Here, we profiled and compared the metabolic characteristics of bone marrow cells of newly weaned 3-week-old offspring of dams fed either a high-fat (Off-HFD) or a regular diet (Off-RD). We utilized a state-of-the-art flow cytometry, and targeted metabolomics approach coupled with a Seahorse metabolic analyzer. RESULTS: We revealed significant metabolic perturbation in the offspring of HFD-fed vs. RD-fed dams, including utilization of glucose primarily via oxidative phosphorylation. We also show a reduction in levels of amino acids, a phenomenon previously linked to bone marrow aging. Using flow cytometry, we found changes in the immune complexity of bone marrow cells and identified a unique B cell population expressing CD19 and CD11b in the bone marrow of three-week-old offspring of high-fat diet-fed mothers. Our data also revealed increased expression of Cyclooxygenase-2 (COX-2) on myeloid CD11b, and on CD11bhi B cells. CONCLUSIONS: Altogether, we demonstrate that the offspring of mothers with obesity show metabolic and immune changes in the bone marrow at a very young age and prior to any symptomatic metabolic disease.

Machine Learning Links T-cell Function and Spatial Localization to Neoadjuvant Immunotherapy and Clinical Outcome in Pancreatic Cancer.

Tumor molecular data sets are becoming increasingly complex, making it nearly impossible for humans alone to effectively analyze them. Here, we demonstrate the power of using machine learning (ML) to analyze a single-cell, spatial, and highly multiplexed proteomic data set from human pancreatic cancer and reveal underlying biological mechanisms that may contribute to clinical outcomes. We designed a multiplex immunohistochemistry antibody panel to compare T-cell functionality and spatial localization in resected tumors from treatment-naïve patients with localized pancreatic ductal adenocarcinoma (PDAC) with resected tumors from a second cohort of patients treated with neoadjuvant agonistic CD40 (anti-CD40) monoclonal antibody therapy. In total, nearly 2.5 million cells from 306 tissue regions collected from 29 patients across both cohorts were assayed, and over 1,000 tumor microenvironment (TME) features were quantified. We then trained ML models to accurately predict anti-CD40 treatment status and disease-free survival (DFS) following anti-CD40 therapy based on TME features. Through downstream interpretation of the ML models' predictions, we found anti-CD40 therapy reduced canonical aspects of T-cell exhaustion within the TME, as compared with treatment-naïve TMEs. Using automated clustering approaches, we found improved DFS following anti-CD40 therapy correlated with an increased presence of CD44+CD4+ Th1 cells located specifically within cellular neighborhoods characterized by increased T-cell proliferation, antigen experience, and cytotoxicity in immune aggregates. Overall, our results demonstrate the utility of ML in molecular cancer immunology applications, highlight the impact of anti-CD40 therapy on T cells within the TME, and identify potential candidate biomarkers of DFS for anti-CD40-treated patients with PDAC.

Metabolic abnormalities in the bone marrow cells of young offspring born to obese mothers.

Intrauterine metabolic reprogramming occurs in obese mothers during gestation, putting the offspring at high risk of developing obesity and associated metabolic disorders even before birth. We have generated a mouse model of maternal high-fat diet-induced obesity that recapitulates the metabolic changes seen in humans. Here, we profiled and compared the metabolic characteristics of bone marrow cells of newly weaned 3-week-old offspring of dams fed either a high-fat (Off-HFD) or a regular diet (Off-RD). We utilized a state-of-the-art targeted metabolomics approach coupled with a Seahorse metabolic analyzer. We revealed significant metabolic perturbation in the offspring of HFD-fed vs. RD-fed dams, including utilization of glucose primarily via oxidative phosphorylation, and reduction in levels of amino acids, a phenomenon previously linked to aging. Furthermore, in the bone marrow of three-week-old offspring of high-fat diet-fed mothers, we identified a unique B cell population expressing CD19 and CD11b, and found increased expression of Cyclooxygenase-2 (COX-2) on myeloid CD11b, and on CD11b hi B cells, with all the populations being significantly more abundant in offspring of dams fed HFD but not a regular diet. Altogether, we demonstrate that the offspring of obese mothers show metabolic and immune changes in the bone marrow at a very young age and prior to any symptomatic metabolic disease.

Machine learning links T cell function and spatial localization to neoadjuvant immunotherapy and clinical outcome in pancreatic cancer.

Tumor molecular datasets are becoming increasingly complex, making it nearly impossible for humans alone to effectively analyze them. Here, we demonstrate the power of using machine learning to analyze a single-cell, spatial, and highly multiplexed proteomic dataset from human pancreatic cancer and reveal underlying biological mechanisms that may contribute to clinical outcome. A novel multiplex immunohistochemistry antibody panel was used to audit T cell functionality and spatial localization in resected tumors from treatment-naive patients with localized pancreatic ductal adenocarcinoma (PDAC) compared to a second cohort of patients treated with neoadjuvant agonistic CD40 (αCD40) monoclonal antibody therapy. In total, nearly 2.5 million cells from 306 tissue regions collected from 29 patients across both treatment cohorts were assayed, and more than 1,000 tumor microenvironment (TME) features were quantified. We then trained machine learning models to accurately predict αCD40 treatment status and disease-free survival (DFS) following αCD40 therapy based upon TME features. Through downstream interpretation of the machine learning models' predictions, we found αCD40 therapy to reduce canonical aspects of T cell exhaustion within the TME, as compared to treatment-naive TMEs. Using automated clustering approaches, we found improved DFS following αCD40 therapy to correlate with the increased presence of CD44+ CD4+ Th1 cells located specifically within cellular spatial neighborhoods characterized by increased T cell proliferation, antigen-experience, and cytotoxicity in immune aggregates. Overall, our results demonstrate the utility of machine learning in molecular cancer immunology applications, highlight the impact of αCD40 therapy on T cells within the TME, and identify potential candidate biomarkers of DFS for αCD40-treated patients with PDAC.

Dual-modality imaging of immunofluorescence and imaging mass cytometry for whole-slide imaging and accurate segmentation.

Imaging mass cytometry (IMC) is a powerful technique capable of detecting over 30 markers on a single slide. It has been increasingly used for single-cell-based spatial phenotyping in a wide range of samples. However, it only acquires a rectangle field of view (FOV) with a relatively small size and low image resolution, which hinders downstream analysis. Here, we reported a highly practical dual-modality imaging method that combines high-resolution immunofluorescence (IF) and high-dimensional IMC on the same tissue slide. Our computational pipeline uses the whole-slide image (WSI) of IF as a spatial reference and integrates small-FOV IMC into a WSI of IMC. The high-resolution IF images enable accurate single-cell segmentation to extract robust high-dimensional IMC features for downstream analysis. We applied this method in esophageal adenocarcinoma of different stages, identified the single-cell pathology landscape via reconstruction of WSI IMC images, and demonstrated the advantage of the dual-modality imaging strategy.

MYC Deregulation and PTEN Loss Model Tumor and Stromal Heterogeneity of Aggressive Triple-Negative Breast Cancer.

Triple-negative breast cancer (TNBC) patients have a poor prognosis and few treatment options. Mouse models of TNBC are important for development of new therapies, however, few mouse models represent the complexity of TNBC. Here, we develop a female TNBC murine model by mimicking two common TNBC mutations with high co-occurrence: amplification of the oncogene MYC and deletion of the tumor suppressor PTEN. This Myc;Ptenfl model develops heterogeneous triple-negative mammary tumors that display histological and molecular features commonly found in human TNBC. Our research involves deep molecular and spatial analyses on Myc;Ptenfl tumors including bulk and single-cell RNA-sequencing, and multiplex tissue-imaging. Through comparison with human TNBC, we demonstrate that this genetic mouse model develops mammary tumors with differential survival and therapeutic responses that closely resemble the inter- and intra-tumoral and microenvironmental heterogeneity of human TNBC, providing a pre-clinical tool for assessing the spectrum of patient TNBC biology and drug response.

Digitize your Biology! Modeling multicellular systems through interpretable cell behavior.

Cells are fundamental units of life, constantly interacting and evolving as dynamical systems. While recent spatial multi-omics can quantitate individual cells' characteristics and regulatory programs, forecasting their evolution ultimately requires mathematical modeling. We develop a conceptual framework-a cell behavior hypothesis grammar-that uses natural language statements (cell rules) to create mathematical models. This allows us to systematically integrate biological knowledge and multi-omics data to make them computable. We can then perform virtual "thought experiments" that challenge and extend our understanding of multicellular systems, and ultimately generate new testable hypotheses. In this paper, we motivate and describe the grammar, provide a reference implementation, and demonstrate its potential through a series of examples in tumor biology and immunotherapy. Altogether, this approach provides a bridge between biological, clinical, and systems biology researchers for mathematical modeling of biological systems at scale, allowing the community to extrapolate from single-cell characterization to emergent multicellular behavior.

An open protocol for modeling T Cell Clonotype repertoires using TCRß CDR3 sequences.

T cell receptor repertoires can be profiled using next generation sequencing (NGS) to measure and monitor adaptive dynamical changes in response to disease and other perturbations. Genomic DNA-based bulk sequencing is cost-effective but necessitates multiplex target amplification using multiple primer pairs with highly variable amplification efficiencies. Here, we utilize an equimolar primer mixture and propose a single statistical normalization step that efficiently corrects for amplification bias post sequencing. Using samples analyzed by both our open protocol and a commercial solution, we show high concordance between bulk clonality metrics. This approach is an inexpensive and open-source alternative to commercial solutions.

Challenges and opportunities for modeling aging and cancer.

Age is among the main risk factors for cancer, and any cancer study in adults is faced with an aging tissue and organism. Yet, pre-clinical studies are carried out using young mice and are not able to address the impact of aging and associated comorbidities on disease biology and treatment outcomes. Here, we discuss the limitations of current mouse cancer models and suggest strategies for developing novel models to address these major gaps in knowledge and experimental approaches.

Supervised learning of high-confidence phenotypic subpopulations from single-cell data.

Accurately identifying phenotype-relevant cell subsets from heterogeneous cell populations is crucial for delineating the underlying mechanisms driving biological or clinical phenotypes. Here, by deploying a learning with rejection strategy, we developed a novel supervised learning framework called PENCIL to identify subpopulations associated with categorical or continuous phenotypes from single-cell data. By embedding a feature selection function into this flexible framework, for the first time, we were able to select informative features and identify cell subpopulations simultaneously, which enables the accurate identification of phenotypic subpopulations otherwise missed by methods incapable of concurrent gene selection. Furthermore, the regression mode of PENCIL presents a novel ability for supervised phenotypic trajectory learning of subpopulations from single-cell data. We conducted comprehensive simulations to evaluate PENCILs versatility in simultaneous gene selection, subpopulation identification and phenotypic trajectory prediction. PENCIL is fast and scalable to analyze 1 million cells within 1 hour. Using the classification mode, PENCIL detected T-cell subpopulations associated with melanoma immunotherapy outcomes. Moreover, when applied to scRNA-seq of a mantle cell lymphoma patient with drug treatment across multiple time points, the regression mode of PENCIL revealed a transcriptional treatment response trajectory. Collectively, our work introduces a scalable and flexible infrastructure to accurately identify phenotype-associated subpopulations from single-cell data.

Dual-modality imaging of immunofluorescence and imaging mass cytometry for whole slide imaging with accurate single-cell segmentation.

Imaging mass cytometry (IMC) is a powerful multiplexed tissue imaging technology that allows simultaneous detection of more than 30 makers on a single slide. It has been increasingly used for singlecell-based spatial phenotyping in a wide range of samples. However, it only acquires a small, rectangle field of view (FOV) with a low image resolution that hinders downstream analysis. Here, we reported a highly practical dual-modality imaging method that combines high-resolution immunofluorescence (IF) and high-dimensional IMC on the same tissue slide. Our computational pipeline uses the whole slide image (WSI) of IF as a spatial reference and integrates small FOVs IMC into a WSI of IMC. The high-resolution IF images enable accurate single-cell segmentation to extract robust high-dimensional IMC features for downstream analysis. We applied this method in esophageal adenocarcinoma of different stages, identified the single-cell pathology landscape via reconstruction of WSI IMC images, and demonstrated the advantage of the dual-modality imaging strategy. Motivation: Highly multiplexed tissue imaging allows visualization of the spatially resolved expression of multiple proteins at the single-cell level. Although imaging mass cytometry (IMC) using metal isotope-conjugated antibodies has a significant advantage of low background signal and absence of autofluorescence or batch effect, it has a low resolution that hampers accurate cell segmentation and results in inaccurate feature extraction. In addition, IMC only acquires mm 2 -sized rectangle regions, which limits its application and efficiency when studying larger clinical samples with non-rectangle shapes. To maximize the research output of IMC, we developed the dual-modality imaging method based on a highly practical and technical improvement requiring no extra specialized equipment or agents and proposed a comprehensive computational pipeline that combines IF and IMC. The proposed method greatly improves the accuracy of cell segmentation and downstream analysis and is able to obtain whole slide image IMC to capture the comprehensive cellular landscape of large tissue sections.

T cell egress via lymphatic vessels is tuned by antigen encounter and limits tumor control.

Antigen-specific CD8+ T cell accumulation in tumors is a prerequisite for effective immunotherapy, and yet the mechanisms of lymphocyte transit are not well defined. Here we show that tumor-associated lymphatic vessels control T cell exit from tumors via the chemokine CXCL12, and intratumoral antigen encounter tunes CXCR4 expression by effector CD8+ T cells. Only high-affinity antigen downregulates CXCR4 and upregulates the CXCL12 decoy receptor, ACKR3, thereby reducing CXCL12 sensitivity and promoting T cell retention. A diverse repertoire of functional tumor-specific CD8+ T cells, therefore, exit the tumor, which limits the pool of CD8+ T cells available to exert tumor control. CXCR4 inhibition or loss of lymphatic-specific CXCL12 boosts T cell retention and enhances tumor control. These data indicate that strategies to limit T cell egress might be an approach to boost the quantity and quality of intratumoral T cells and thereby response to immunotherapy.

An open protocol for modeling T Cell Clonotype repertoires using TCRß CDR3 sequences.

T cell receptor repertoires can be profiled using next generation sequencing (NGS) to measure and monitor adaptive dynamical changes in response to disease and other perturbations. Genomic DNA-based bulk sequencing is cost-effective but necessitates multiplex target amplification using multiple primer pairs with highly variable amplification efficiencies. Here, we utilize an equimolar primer mixture and propose a single statistical normalization step that efficiently corrects for amplification bias post sequencing. Using samples analyzed by both our open protocol and a commercial solution, we show high concordance between bulk clonality metrics. This approach is an inexpensive and open-source alternative to commercial solutions.

Alignment, segmentation and neighborhood analysis in cyclic immunohistochemistry data using CASSATT.

Cyclic immunohistochemistry (cycIHC) uses sequential rounds of colorimetric immunostaining and imaging for quantitative mapping of location and number of cells of interest. Additionally, cycIHC benefits from the speed and simplicity of brightfield microscopy, making the collection of entire tissue sections and slides possible at a trivial cost compared to other high dimensional imaging modalities. However, large cycIHC datasets currently require an expert data scientist to concatenate separate open-source tools for each step of image pre-processing, registration, and segmentation, or the use of proprietary software. Here, we present a unified and user-friendly pipeline for processing, aligning, and analyzing cycIHC data - Cyclic Analysis of Single-Cell Subsets and Tissue Territories (CASSATT). CASSATT registers scanned slide images across all rounds of staining, segments individual nuclei, and measures marker expression on each detected cell. Beyond straightforward single cell data analysis outputs, CASSATT explores the spatial relationships between cell populations. By calculating the log odds of interaction frequencies between cell populations within tissues and tissue regions, this pipeline helps users identify populations of cells that interact-or do not interact-at frequencies that are greater than those occurring by chance. It also identifies specific neighborhoods of cells based on the assortment of neighboring cell types that surround each cell in the sample. The presence and location of these neighborhoods can be compared across slides or within distinct regions within a tissue. CASSATT is a fully open source workflow tool developed to process cycIHC data and will allow greater utilization of this powerful staining technique.

Modulation of myeloid and T cells in vivo by Bruton's tyrosine kinase inhibitor ibrutinib in patients with metastatic pancreatic ductal adenocarcinoma.

BACKGROUND: In preclinical studies of pancreatic ductal adenocarcinoma (PDAC), ibrutinib improved the antitumor efficacy of the standard of care chemotherapy. This led to a phase 1b clinical trial to determine the safety, tolerability, and immunologic effects of ibrutinib treatment in patients with advanced PDAC. METHODS: Previously untreated patients with PDAC were enrolled in a phase 1b clinical trial (ClinicalTrials.gov) to determine the safety, toxicity, and maximal tolerated dose of ibrutinib when administered with the standard regimen of gemcitabine and nab-paclitaxel. To study the immune response to ibrutinib alone, the trial included an immune response arm where patients were administered with ibrutinib daily for a week followed by ibrutinib combined with gemcitabine and nab-paclitaxel. Endoscopic ultrasonography-guided primary PDAC tumor biopsies and blood were collected before and after ibrutinib monotherapy. Changes in abundance and functional state of immune cells in the blood was evaluated by mass cytometry by time of flight and statistical scaffold analysis, while that in the local tumor microenvironment (TME) were assessed by multiplex immunohistochemistry. Changes in B-cell receptor and T-cell receptor repertoire were assessed by sequencing and analysis of clonality. RESULTS: In the blood, ibrutinib monotherapy significantly increased the frequencies of activated inducible T cell costimulator+(ICOS+) CD4+ T cells and monocytes. Within the TME, ibrutinib monotherapy led to a trend in decreased B-cell abundance but increased interleukin-10+ B-cell frequency. Monotherapy also led to a trend in increased mature CD208+dendritic cell density, increased late effector (programmed cell death protein 1 (PD-1-) eomesodermin (EOMES+)) CD8+ T-cell frequency, with a concomitantly decreased dysfunctional (PD-1+ EOMES+) CD8+ T-cell frequency. When ibrutinib was combined with chemotherapy, most of these immune changes were not observed. Patients with partial clinical responses had more diverse T and B cell receptor repertoires prior to therapy initiation. CONCLUSION: Ibrutinib monotherapy skewed the immune landscape both in the circulation and TME towards activated T cells, monocytes and DCs. These effects were not observed when combining ibrutinib with standard of care chemotherapy. Future studies may focus on other therapeutic combinations that augment the immunomodulatory effects of ibrutinib in solid tumors. TRIAL REGISTRATION NUMBER: NCT02562898.