Quantitative, qualitative and spatial analysis of lymphocyte infiltration in periampullary and pancreatic adenocarcinoma

Spatial distance between tumor and lymphocyte can predict the survival of patients with resectable lung adenocarcinoma

Background and objective

Spatial interaction between tumor-infiltrating lymphocytes (TILs) and tumor cells is valuable in predicting the effectiveness of immune response and prognosis amongst patients with lung adenocarcinoma (LUAD). Recent evidence suggests that the spatial distance between tumor cells and lymphocytes also influences the immune responses, but the distance analysis based on Hematoxylin and Eosin (H&E) -stained whole-slide images (WSIs) remains insufficient. To address this issue, we aim to explore the relationship between distance and prognosis prediction of patients with LUAD in this study.

Methods

We recruited patients with resectable LUAD from three independent cohorts in this multi-center study. We proposed a simple but effective deep learning-driven workflow to automatically segment different cell types in the tumor region using the HoVer-Net model, and quantified the spatial distance (DIST) between tumor cells and lymphocytes based on H&E-stained WSIs. The association of DIST with disease-free survival (DFS) was explored in the discovery set (D1, n = 276) and the two validation sets (V1, n = 139; V2, n = 115).

Results

In multivariable analysis, the low DIST group was associated with significantly better DFS in the discovery set (D1, HR, 0.61; 95 % CI, 0.40-0.94; p = 0.027) and the two validation sets (V1, HR, 0.54; 95 % CI, 0.32-0.91; p = 0.022; V2, HR, 0.44; 95 % CI, 0.24-0.81; p = 0.009). By integrating the DIST with clinicopathological factors, the integrated model (full model) had better discrimination for DFS in the discovery set (C-index, D1, 0.745 vs. 0.723) and the two validation sets (V1, 0.621 vs. 0.596; V2, 0.671 vs. 0.650). Furthermore, the computerized DIST was associated with immune phenotypes such as immune-desert and inflamed phenotypes.

Conclusions

The integration of DIST with clinicopathological factors could improve the stratification performance of patients with resectable LUAD, was beneficial for the prognosis prediction of LUAD patients, and was also expected to assist physicians in individualized treatment.

Longitudinal molecular profiling elucidates immunometabolism dynamics in breast cancer

Although metabolic reprogramming within tumor cells and tumor microenvironment (TME) is well described in breast cancer, little is known about how the interplay of immune state and cancer metabolism evolves during treatment. Here, we characterize the immunometabolic profiles of tumor tissue samples longitudinally collected from individuals with breast cancer before, during and after neoadjuvant chemotherapy (NAC) using proteomics, genomics and histopathology. We show that the pre-, on-treatment and dynamic changes of the immune state, tumor metabolic proteins and tumor cell gene expression profiling-based metabolic phenotype are associated with treatment response. Single-cell/nucleus RNA sequencing revealed distinct tumor and immune cell states in metabolism between cold and hot tumors. Potential drivers of NAC based on above analyses were validated in vitro. In summary, the study shows that the interaction of tumor-intrinsic metabolic states and TME is associated with treatment outcome, supporting the concept of targeting tumor metabolism for immunoregulation.

Comprehensive analysis and identification of prognostic biomarkers and immunotherapeutic targets in the NADPH oxidase family (and its regulatory subunits) in pancreatic ductal adenocarcinoma

PURPOSE

This study aimed to evaluate the role of NADPH in pancreatic ductal adenocarcinoma using bioinformatic analyses and experimental validations.

METHODS

We compared the expression levels, performed GO and KEGG analysis of NADPH oxidase family and its regulatory subunits, and determined the survival of patients with pancreatic ductal adenocarcinoma by GEPIA, David and KM plotter. The relationship between their expression with immune infiltration levels, phagocytotic/NK cell immune checkpoints, recruitment-related molecules were detected by Timer 2.0 and TISIDB, respectively. Subsequently, their correlation with NK cell infiltration level was verified by immunohistochemistry.

RESULTS

The expression of some members of the NADPH oxidase family and its regulatory subunits was significantly increased in pancreatic ductal adenocarcinoma tissues compared to that in normal tissues and was positively correlated with natural killer (NK) cell infiltration. Furthermore, the NADPH oxidase family and its regulatory subunits were associated with survival and immune status in patients with pancreatic ductal adenocarcinoma, including chemokines, immune checkpoints, and immune infiltration levels of NK cells, monocytes, and myeloid-derived suppressor cells.

CONCLUSIONS

These results suggest the NADPH oxidase family and its regulatory subunits might serve as indicators for predicting the responsiveness to immunotherapy and outcome of patients with pancreatic ductal adenocarcinoma, providing a new perspective or strategy for immunotherapy in pancreatic ductal adenocarcinoma.

B cell immune profiles in dysbiotic vermiform appendixes of pancreatic cancer patients

Pancreatic ductal adenocarcinoma (PDAC) remains one of the deadliest solid tumors and is resistant to immunotherapy. B cells play an essential role in PDAC progression and immune responses, both locally and systemically. Moreover, increasing evidence suggests that microbial compositions inside the tumor, as well as in the oral cavity and the gut, are important factors in shaping the PDAC immune landscape. However, the gut-associated lymphoid tissue (GALT) has not previously been explored in PDAC patients. In this study, we analyzed healthy vermiform appendix (VA) from 20 patients with PDAC and 32 patients with colon diseases by gene expression immune profiling, flow cytometry analysis, and microbiome sequencing. We show that the VA GALT of PDAC patients exhibits markers of increased inflammation and cytotoxic cell activity. In contrast, B cell function is decreased in PDAC VA GALT based on gene expression profiling; B cells express significantly fewer MHC class II surface receptors, whereas plasma cells express the immune checkpoint molecule HLA-G. Additionally, the vermiform appendix microbiome of PDAC patients is enriched with Klebsiella pneumoniae, Bifidobacterium animalis, and Adlercreutzia equolifaciens, while certain commensals are depleted. Our findings may suggest impaired B cell function within the GALT of PDAC patients, which could potentially be linked to microbial dysbiosis. Additional investigations are imperative to validate our observations and explore these potential targets of future therapies.

Tertiary lymphoid structures in pancreatic cancer: a new target for immunotherapy

Pancreatic cancer (PC) is extremely malignant and shows limited response to available immunotherapies due to the hypoxic and immunosuppressive nature of its tumor microenvironment (TME). The aggregation of immune cells (B cells, T cells, dendritic cells, etc.), which is induced in various chronic inflammatory settings such as infection, inflammation, and tumors, is known as the tertiary lymphoid structure (TLS). Several studies have shown that TLSs can be found in both intra- and peritumor tissues of PC. The role of TLSs in peritumor tissues in tumors remains unclear, though intratumoral TLSs are known to play an active role in a variety of tumors, including PC. The formation of intratumoral TLSs in PC is associated with a good prognosis. In addition, TLSs can be used as an indicator to assess the effectiveness of treatment. Targeted induction of TLS formation may become a new avenue of immunotherapy for PC. This review summarizes the formation, characteristics, relevant clinical outcomes, and clinical applications of TLSs in the pancreatic TME. We aim to provide new ideas for future immunotherapy of PC.

Role of Immune Microenvironment in Pancreatic Ductal Adenocarcinoma: Could It Be Considered a Predictor of Prognosis?

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) is characterized by a highly immunosuppressive tumor microenvironment (TME). The aim of this study is to determine the potential significant TME immune markers of long-term survival.

METHODS

We retrospectively included patients with a diagnosis of resectable PDAC having undergone upfront surgery. Immunohistochemical (IHC) staining using tissue microarray for PD-L1, CD3, CD4, CD8, FOXP3, CD20, iNOS and CD163 was performed in order to characterize the TME. The primary endpoint was long-term survival, defined as the Overall Survival > 24 months from surgery.

RESULTS

A total of 38 consecutive patients were included, and 14 (36%) of them were long-term survivors. Long-term survivors showed a higher density of CD8+ lymphocytes intra- and peri-acinar (p = 0.08), and a higher CD8/FOXP3 intra- and peri-tumoral ratio (p = 0.05). A low density of intra- and peri-tumoral FOXP3 infiltration is a good predictor of long-term survival (p = 0.04). A significant association of the low density of intra- and peri-tumoral tumor-associated macrophages (TAMs) iNOS+ with long-term survival was detected (p = 0.04).

CONCLUSIONS

Despite the retrospective nature and small sample size, our study showed that the high infiltration of CD8+ lymphocytes and low infiltration of FOXP3+ and TAMs iNOS+ are predictors of good prognosis. A preoperative assessment of these potential immune markers could be useful and determinant in the staging process and in PDAC management.

mosna reveals different types of cellular interactions predictive of response to immunotherapies in cancer

Single-cell spatially resolved proteomic or transcriptomic methods offer the opportunity to discover cell types interactions of biological or clinical importance. To extract relevant information from these data, we present mosna, a Python package to analyze spatially resolved experiments and discover patterns of cellular spatial organization. It includes the detection of preferential interactions between specific cell types and the discovery of cellular niches. We exemplify the proposed analysis pipeline on spatially resolved proteomic data from cancer patient samples annotated with clinical response to immunotherapy, and we show that mosna can identify a number of features describing cellular composition and spatial distribution that can provide biological hypotheses regarding factors that affect response to therapies.

B cells in the tumor microenvironment: Multi-faceted organizers, regulators, and effectors of anti-tumor immunity

Our understanding of tumor-infiltrating lymphocytes (TILs) is rapidly expanding beyond T cell-centric perspectives to include B cells and plasma cells, collectively referred to as TIL-Bs. In many cancers, TIL-Bs carry strong prognostic significance and are emerging as key predictors of response to immune checkpoint inhibitors. TIL-Bs can perform multiple functions, including antigen presentation and antibody production, which allow them to focus immune responses on cognate antigen to support both T cell responses and innate mechanisms involving complement, macrophages, and natural killer cells. In the stroma of the most immunologically "hot" tumors, TIL-Bs are prominent components of tertiary lymphoid structures, which resemble lymph nodes structurally and functionally. Additionally, TIL-Bs participate in a variety of other lympho-myeloid aggregates and engage in dynamic interactions with the tumor stroma. Here, we summarize our current understanding of TIL-Bs in human cancer, highlighting the compelling therapeutic opportunities offered by their unique tumor recognition and effector mechanisms.

An immune score reflecting pro- and anti-tumoural balance of tumour microenvironment has major prognostic impact and predicts immunotherapy response in solid cancers

BACKGROUND

Cancer immunity is based on the interaction of a multitude of cells in the spatial context of the tumour tissue. Clinically relevant immune signatures are therefore anticipated to fundamentally improve the accuracy in predicting disease progression.

METHODS

Through a multiplex in situ analysis we evaluated 15 immune cell classes in 1481 tumour samples. Single-cell and bulk RNAseq data sets were used for functional analysis and validation of prognostic and predictive associations.

FINDINGS

By combining the prognostic information of anti-tumoural CD8⁺ lymphocytes and tumour supportive CD68⁺CD163⁺ macrophages in colorectal cancer we generated a signature of immune activation (SIA). The prognostic impact of SIA was independent of conventional parameters and comparable with the state-of-art immune score. The SIA was also associated with patient survival in oesophageal adenocarcinoma, bladder cancer, lung adenocarcinoma and melanoma, but not in endometrial, ovarian and squamous cell lung carcinoma. We identified CD68⁺CD163⁺ macrophages as the major producers of complement C1q, which could serve as a surrogate marker of this macrophage subset. Consequently, the RNA-based version of SIA (ratio of CD8A to C1QA) was predictive for survival in independent RNAseq data sets from these six cancer types. Finally, the CD8A/C1QA mRNA ratio was also predictive for the response to checkpoint inhibitor therapy.

INTERPRETATION

Our findings extend current concepts to procure prognostic information from the tumour immune microenvironment and provide an immune activation signature with high clinical potential in common human cancer types.

FUNDING

Swedish Cancer Society, Lions Cancer Foundation, Selanders Foundation, P.O. Zetterling Foundation, U-CAN supported by SRA CancerUU, Uppsala University and Region Uppsala.

The prognostic impact of tumor-infiltrating B lymphocytes in patients with solid malignancies: A systematic review and meta-analysis

This study reviewed the prognostic effect of tumor-infiltrating B lymphocytes (TIBLs) on solid malignancies, to determine the potential role of TIBLs in predicting cancer patient's prognosis and their response to immunotherapy. A total of 45 original papers involving 11,099 individual patients were included in this meta-analysis covering 7 kinds of cancer. The pooled results suggested that high levels of TIBLs were correlated with favorable OS in lung, esophageal, gastric, colorectal, liver, and breast cancer; improved RFS in lung cancer; and improved DFS in gastrointestinal neoplasms. Additionally, TIBLs were significantly correlated with negative lymphatic invasion in gastric cancer, small tumor size in hepatocellular carcinoma, and negative distant metastasis in colorectal cancer. Additionally, TIBLs were reported as a discriminative feature of patients treated with immunotherapy with improved survival. We concluded that TIBLs play a favorable prognostic role among the common solid malignancie, providing theoretical evidence for further prognosis prediction for solid tumors.

Expression of cartilage oligomeric matrix protein in periampullary adenocarcinoma is associated with pancreatobiliary-type morphology, higher levels of fibrosis and immune cell exclusion

Cartilage oligomeric matrix protein (COMP) is an emerging regulator of tumor progression. The aim of this study was to evaluate the expression of COMP in periampullary adenocarcinoma with respect to prognostic value for survival and relapse, levels of fibrosis and infiltrating immune cells. COMP expression was evaluated using immunohistochemistry in primary tumors and subsets of paired lymph node metastases in tissue microarrays including 175 patients with periampullary adenocarcinoma. Collagen content was assessed with Sirius Red-Fast Green staining. High COMP levels were detected in cancer cells and in stroma, in 46% and 57% of the patients, respectively. High COMP expression was strongly associated with more aggressive pancreatobiliary-type (PB-type) compared to intestinal-type tumors (p < .0001). Importantly, high expression of COMP correlated with the exclusion of cytotoxic T-cells from the cancer cell compartment of the tumors, particularly in PB-type tumors. Higher levels of fibrosis measured by the density of collagen fibers correlated with high COMP levels in both cancer cells and stroma. This in turn could lead to exclusion of cytotoxic T-cells from accessing the cancer cells, a recognized immunotherapy resistance mechanism. Targeting COMP could therefore be considered as a novel therapeutic strategy in PB-type periampullary adenocarcinoma.

Therapeutic Targets and Prognostic Biomarkers Among CXC Chemokines in Pancreatic Ductal Adenocarcinoma Microenvironment

OBJECTIVES

Pancreatic ductal adenocarcinoma (PDAC) is characterized by occult onset, rapid progression, and poor prognosis. CXC chemokines play an important role in tumor microenvironment and development. However, the potential mechanistic values of CXC chemokines as clinical biomarkers and therapeutic targets in PDAC have not been fully clarified.

METHODS

The altered expression, interaction network, and clinical data of CXC chemokines in patients with PDAC were analyzed by the data from the Gene Expression Omnibus and the Tumor Cancer Genome Atlas.

RESULTS

CXCL5 transcriptional level was significantly elevated in PDAC tissues. A significant correlation was found between the expression of CXC1/3/5/8 and the pathological stage of PDAC patients. The PDAC patients with low transcriptional levels of CXCL5/9/10/11/17 were associated with a significantly better prognosis. The functions of differentially expressed CXC chemokines are primarily related to the chemokine signaling pathway, cytokine-cytokine receptor interaction, and viral protein interaction with cytokine and cytokine receptor. RELA, NFKB1, and SP1 are key transcription factors for CXC chemokines, and the SRC family of tyrosine kinases, mitogen-activated protein kinases, CDK5, PRKCQ, ROCK1, ITK, IKBKE, JAK3, and NTRK2 are CXC chemokine targets.

CONCLUSIONS

The results indicated that CXC chemokines might serve as therapeutic targets and prognostic biomarkers in PDAC.

Branching Copy-Number Evolution and Parallel Immune Profiles across the Regional Tumor Space of Resected Pancreatic Cancer

Pancreatic ductal adenocarcinoma (PDAC) remains a highly lethal disease. The only option for curative treatment is resection of the tumor followed by standard adjuvant chemotherapy. Yet, early relapse due to chemoresistance is almost inevitable. Herein, we delineated the genetic intratumor heterogeneity in resected PDAC, with the aim to identify evolutionary patterns that may be associated with overall survival (OS) following treatment with curative intent. Potential relationships with the adjacent immune microenvironment were also examined. The genetic and immune landscapes of the regional tumor space were analyzed in nine patients with resected PDAC. Targeted deep sequencing and genome wide SNP array were followed by clonal deconvolution and phylogenetic analysis. A mathematical complexity score was developed to calculate the network extent of each phylogeny. Spatial variation in abundancy and tumor nest infiltration of immune cells was analyzed by double IHC staining. Copy-number heterogeneity was denoted as the major contributing factor to the branching architectures of the produced phylogenetic trees. Increased tree complexity was significantly inversely associated with OS, and larger regional maximum aberrations (higher treetops) were associated with increased PD-L1 expression on tumor cells. Contrastingly, an FREM1 gene amplification, found in one patient, coincided with a particularly vigorous immune response. Findings from this limited case series suggest that complex evolutionary patterns may be associated with a shorter survival in surgically treated patients with PDAC. Some hypothesis-generating associations with the surrounding immune microenvironment were also detected.

IMPLICATIONS

Evolutionary copy-number patterns may be associated with survival in patients with resected PDAC.

NOX4 regulates TGFβ-induced proliferation and self-renewal in glioblastoma stem cells

Glioblastoma (GBM) is the most aggressive and common glioma subtype, with a median survival of 15 months after diagnosis. Current treatments have limited therapeutic efficacy; thus, more effective approaches are needed. The glioblastoma tumoural mass is characterised by a small cellular subpopulation – glioblastoma stem cells (GSCs) – that has been held responsible for glioblastoma initiation, cell invasion, proliferation, relapse and resistance to chemo‐ and radiotherapy. Targeted therapies against GSCs are crucial, as is understanding the molecular mechanisms that govern the GSCs. Transforming growth factor β (TGFβ) signalling and reactive oxygen species (ROS) production are known to govern and regulate cancer stem cell biology. Among the differentially expressed genes regulated by TGFβ in a transcriptomic analysis of two different patient‐derived GSCs, we found NADPH oxidase 4 (NOX4) as one of the top upregulated genes. Interestingly, when patient tissues were analysed, NOX4 expression was found to be higher in GSCs versus differentiated cells. A functional analysis of the role of NOX4 downstream of TGFβ in several patient‐derived GSCs showed that TGFβ does indeed induce NOX4 expression and increases ROS production in a NOX4‐dependent manner. NOX4 downstream of TGFβ regulates GSC proliferation, and NOX4 expression is necessary for TGFβ‐induced expression of stem cell markers and of the transcription factor nuclear factor erythroid 2‐related factor 2 (NRF2), which in turn controls the cell’s antioxidant and metabolic responses. Interestingly, overexpression of NOX4 recapitulates the effects induced by TGFβ in GSCs: enhanced proliferation, stemness and NRF2 expression. In conclusion, this work functionally establishes NOX4 as a key mediator of GSC biology.

Spatial Distribution and Predictive Significance of Dendritic Cells and Macrophages in Esophageal Cancer Treated With Combined Chemoradiotherapy and PD-1 Blockade

Background

The first clinical study (NCT03671265) of first-line chemoradiotherapy combined with PD-1 blockade showed promising treatment outcomes in locally advanced esophageal squamous cell carcinoma (ESCC). However, partial patients did not respond to the combination treatment. The roles of dendritic cells (DCs) and macrophages in this combination treatment remain poorly understood.

Methods

We performed multiplexed immunofluorescence method to identify CD11c⁺ DCs, CD68⁺ macrophages, and their PD-L1^- or PD-L1⁺ subpopulations in paired tumor biopsies (n = 36) collected at baseline and during the combination treatment (after radiation, 40 Gy) from the phase Ib trial (NCT03671265). We applied whole exome sequencing in the baseline tumor biopsies (n = 14) to estimate tumor mutation burden (TMB). We dynamically investigated the spatial distribution of DCs and macrophages under chemoradiotherapy combined with PD-1 blockade, and evaluated the association between their spatial distribution and combination outcome, and TMB.

Results

The results showed that high percentages of PD-L1^- DCs and macrophages in the baseline tumor compartment, but not in the stromal compartment, predicted improved OS and PFS. Chemoradiotherapy combined with PD-1 blockade promoted DCs and macrophages to migrate closer to tumor cells. During combination treatment, PD-L1^- tumor cells were nearest to PD-L1^- DCs and macrophages, while PD-L1⁺ tumor cells were next to PD-L1⁺ DCs and macrophages. High TMB was closely associated with a shorter distance from tumor cells to DCs and macrophages. Shorter distance between PD-L1⁺ tumor cells and PD-L1⁺ DCs or PD-L1^- macrophages during the combination was correlated with better OS. Shorter distance between PD-L1^- tumor cells and PD-L1^- macrophages during combination was associated with both longer OS and PFS.

Conclusions

PD-L1^- or PD-L1⁺ DCs and macrophages exhibit distinct spatial distribution in ESCC. The close distance between tumor cells and these antigen-presenting cells (APCs) is critical to the clinical outcome in chemoradiotherapy combined with PD-1 blockade in ESCC patients. Our results highlight the predictive potential of spatial patterns of APCs in chemoradiotherapy combined with immunotherapy and reveal the underlying mechanism of APCs participating in chemoradiotherapy-induced antitumor immune response in ESCC.

Spatial Immunology in Liver Metastases from Colorectal Carcinoma according to the Histologic Growth Pattern

Simple Summary

In the era of immunotherapy, the tumor microenvironment (TME) has attracted special interest. However, colorectal liver metastases (CRC-LM) present histological peculiarities that could affect the interaction of immune and tumor cells such as fibrotic encapsulation and dense intratumoral stroma. We explored the spatial distribution of lymphocytic infiltrates in CRC-LM in the context of the histologic growth patterns using multispectral digital pathology providing data on three different scenarios, tumor periphery, invasive margin, and central tumoral areas. Our results illustrate a similar poor cell density of CD8⁺ cells between different metastases subtypes in intratumoral regions. However, in encapsulated metastases, cytotoxic cells reach the tumor cells while remaining retained in stromal areas in non-encapsulating metastases. Some aspects are still unresolved, such as understanding the reason why most lymphocytes are largely retained in the capsule.

Abstract

Colorectal cancer liver metastases (CRC-LM) present differential histologic growth patterns (HGP) that determine the interaction between immune and tumor cells. We explored the spatial distribution of lymphocytic infiltrates in CRC-LM in the context of the HGP using multispectral digital pathology. We did not find statistically significant differences of immune cell densities in the central regions of desmoplastic (_dHGP) and non-desmoplastic (_ndHGP) metastases. The spatial evaluation reported that _dHGP-metastases displayed higher infiltration by CD8⁺ and CD20⁺ cells in peripheral regions as well as CD4⁺ and CD45RO⁺ cells in _ndHGP-metastases. However, the reactive stroma regions at the invasive margin (IM) of _ndHGP-metastases displayed higher density of CD4⁺, CD20⁺, and CD45RO⁺ cells. The antitumor status of the TIL infiltrates measured as CD8/CD4 reported higher values in the IM of encapsulated metastases up to 400 μm towards the tumor center (p < 0.05). Remarkably, the IM of _dHGP-metastases was characterized by higher infiltration of CD8⁺ cells in the epithelial compartment parameter assessed with the ratio CD8^epithelial/CD8^stromal, suggesting anti-tumoral activity in the encapsulating lesions. Taking together, the amount of CD8⁺ cells is comparable in the IM of both HGP metastases types. However, in _dHGP-metastases some cytotoxic cells reach the tumor nests while remaining retained in the stromal areas in _ndHGP-metastases.

Interplay between copy number alterations and immune profiles in the early breast cancer Scandinavian Breast Group 2004-1 randomized phase II trial: results from a feasibility study

Emerging data indicate that genomic alterations can shape immune cell composition in early breast cancer. However, there is a need for complementary imaging and sequencing methods for the quantitative assessment of combined somatic copy number alteration (SCNA) and immune profiling in pathological samples. Here, we tested the feasibility of three approaches-CUTseq, for high-throughput low-input SCNA profiling, multiplexed fluorescent immunohistochemistry (mfIHC) and digital-image analysis (DIA) for quantitative immuno-profiling- in archival formalin-fixed paraffin-embedded (FFPE) tissue samples from patients enrolled in the randomized SBG-2004-1 phase II trial. CUTseq was able to reproducibly identify amplification and deletion events with a resolution of 100 kb using only 6 ng of DNA extracted from FFPE tissue and pooling together 77 samples into the same sequencing library. In the same samples, mfIHC revealed that CD4 + T-cells and CD68 + macrophages were the most abundant immune cells and they mostly expressed PD-L1 and PD-1. Combined analysis showed that the SCNA burden was inversely associated with lymphocytic infiltration. Our results set the basis for further applications of CUTseq, mfIHC and DIA to larger cohorts of early breast cancer patients.

The Immune Landscape of Colorectal Cancer

Simple Summary

We sought to provide a detailed overview of the immune landscape of colorectal cancer in the largest study to date in terms of patient numbers and analyzed immune cell types. We applied a multiplex in situ staining method in combination with an advanced scanning and image analysis pipeline akin to flow cytometry, and analyzed 5968 individual multi-layer images of tissue defining in a total of 39,078,450 cells. We considered the location of immune cells with respect to the stroma, and tumor cell compartment and tumor regions in the central part or the invasive margin. To the best of our knowledge, this study is the first comprehensive spatial description of the immune landscape in colorectal cancer using a large population-based cohort and a multiplex immune cell identification.

Abstract

While the clinical importance of CD8+ and CD3+ cells in colorectal cancer (CRC) is well established, the impact of other immune cell subsets is less well described. We sought to provide a detailed overview of the immune landscape of CRC in the largest study to date in terms of patient numbers and in situ analyzed immune cell types. Tissue microarrays from 536 patients were stained using multiplexed immunofluorescence panels, and fifteen immune cell subclasses, representing adaptive and innate immunity, were analyzed. Overall, therapy-naïve CRC patients clustered into an ‘inflamed’ and a ‘desert’ group. Most T cell subsets and M2 macrophages were enriched in the right colon (p-values 0.046–0.004), while pDC cells were in the rectum (p = 0.008). Elderly patients had higher infiltration of M2 macrophages (p = 0.024). CD8+ cells were linked to improved survival in colon cancer stages I-III (q = 0.014), while CD4+ cells had the strongest impact on overall survival in metastatic CRC (q = 0.031). Finally, we demonstrated repopulation of the immune infiltrate in rectal tumors post radiation, following an initial radiation-induced depletion. This study provides a detailed analysis of the in situ immune landscape of CRC paving the way for better diagnostics and providing hints to better target the immune microenvironment.

Discordance of PD-L1 Expression at the Protein and RNA Levels in Early Breast Cancer

Simple Summary

Despite the increasing use of checkpoint inhibitors for early and metastatic breast cancer, Programmed Death Ligand 1 (PD-L1) remains the only validated albeit imperfect predictive biomarker. Significant discordance in PD-L1 protein expression depending on the antibody used has been demonstrated, while the weak correlation and discordant prognostic information between protein and gene expression underscore its biologic heterogeneity. In this study, we use material from two patient cohorts of early breast cancer and multiple methodologies (immunohistochemistry, RNA fluorescent in situ hybridization, immunofluorescence, bulk gene expression, and multiplex fluorescent immunohistochemistry) to demonstrate the significant discordance in PD-L1 expression among various methods and between different areas of the same tumor, which hints toward the presence of spatial, intratumoral and biological heterogeneity.

Abstract

We aimed to assess if the discrepant prognostic information between Programmed Death Ligand 1 (PD-L1) protein versus mRNA expression in early breast cancer (BC) could be attributed to heterogeneity in its expression. PD-L1 protein and mRNA expression in BC tissue microarrays from two clinical patient cohorts were evaluated (105 patients; cohort 1: untreated; cohort 2: neoadjuvant chemotherapy-treated). Immunohistochemistry (IHC) with SP142, SP263 was performed. PD-L1 mRNA was evaluated using bulk gene expression and RNA-FISH RNAscope^®, the latter scored in a semi-quantitative manner and combined with immunofluorescence (IF) staining for the simultaneous detection of PD-L1 protein expression. PD-L1 expression was assessed in cores as a whole and in two regions of interest (ROI) from the same core. The cell origin of PD-L1 expression was evaluated using multiplex fluorescent IHC. IHC PD-L1 expression between SP142 and SP263 was concordant in 86.7% of cores (p < 0.001). PD-L1 IF/IHC was weakly correlated with spatial mRNA expression (concordance 54.6–71.2%). PD-L1 was mostly expressed by lymphocytes intra-tumorally, while its stromal expression was mostly observed in macrophages. Our results demonstrate only moderate concordance between the various methods of assessing PD-L1 expression at the protein and mRNA levels, which may be attributed to both analytical performance and spatial heterogeneity.

Tysserand - Fast and accurate reconstruction of spatial networks from bioimages

SUMMARY

Networks provide a powerful framework to analyze spatial omics experiments. However, we lack tools that integrate several methods to easily reconstruct networks for further analyses with dedicated libraries. In addition, choosing the appropriate method and parameters can be challenging.We propose tysserand, a Python library to reconstruct spatial networks from spatially resolved omics experiments. It is intended as a common tool to which the bioinformatics community can add new methods to reconstruct networks, choose appropriate parameters, clean resulting networks and pipe data to other libraries.

AVAILABILITY AND IMPLEMENTATION

tysserand software and tutorials with a Jupyter notebook to reproduce the results are available at https://github.com/VeraPancaldiLab/tysserand.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Spatial architecture of the immune microenvironment orchestrates tumor immunity and therapeutic response

Tumors are not only aggregates of malignant cells but also well-organized complex ecosystems. The immunological components within tumors, termed the tumor immune microenvironment (TIME), have long been shown to be strongly related to tumor development, recurrence and metastasis. However, conventional studies that underestimate the potential value of the spatial architecture of the TIME are unable to completely elucidate its complexity. As innovative high-flux and high-dimensional technologies emerge, researchers can more feasibly and accurately detect and depict the spatial architecture of the TIME. These findings have improved our understanding of the complexity and role of the TIME in tumor biology. In this review, we first epitomized some representative emerging technologies in the study of the spatial architecture of the TIME and categorized the description methods used to characterize these structures. Then, we determined the functions of the spatial architecture of the TIME in tumor biology and the effects of the gradient of extracellular nonspecific chemicals (ENSCs) on the TIME. We also discussed the potential clinical value of our understanding of the spatial architectures of the TIME, as well as current limitations and future prospects in this novel field. This review will bring spatial architectures of the TIME, an emerging dimension of tumor ecosystem research, to the attention of more researchers and promote its application in tumor research and clinical practice.

Spatial clustering of CD68+ tumor associated macrophages with tumor cells is associated with worse overall survival in metastatic clear cell renal cell carcinoma

Immune infiltration is typically quantified using cellular density, not accounting for cellular clustering. Tumor-associated macrophages (TAM) activate oncogenic signaling through paracrine interactions with tumor cells, which may be better reflected by local cellular clustering than global density metrics. Using multiplex immunohistochemistry and digital pathologic analysis we quantified cellular density and cellular clustering for myeloid cell markers in 129 regions of interest from 55 samples from 35 patients with metastatic ccRCC. CD68+ cells were found to be clustered with tumor cells and dispersed from stromal cells, while CD163+ and CD206+ cells were found to be clustered with stromal cells and dispersed from tumor cells. CD68+ density was not associated with OS, while high tumor/CD68+ cell clustering was associated with significantly worse OS. These novel findings would not have been identified if immune infiltrate was assessed using cellular density alone, highlighting the importance of including spatial analysis in studies of immune cell infiltration of tumors.

Significance: Increased clustering of CD68+ TAMs and tumor cells was associated with worse overall survival for patients with metastatic ccRCC. This effect would not have been identified if immune infiltrate was assessed using cell density alone, highlighting the importance of including spatial analysis in studies of immune cell infiltration of tumors.

The prognostic impact of the tumour stroma fraction: A machine learning-based analysis in 16 human solid tumour types

Background

The development of a reactive tumour stroma is a hallmark of tumour progression and pronounced tumour stroma is generally considered to be associated with clinical aggressiveness. The variability between tumour types regarding stroma fraction, and its prognosis associations, have not been systematically analysed.

Methods

Using an objective machine-learning method we quantified the tumour stroma in 16 solid cancer types from 2732 patients, representing retrospective tissue collections of surgically resected primary tumours. Image analysis performed tissue segmentation into stromal and epithelial compartment based on pan-cytokeratin staining and autofluorescence patterns.

Findings

The stroma fraction was highly variable within and across the tumour types, with kidney cancer showing the lowest and pancreato-biliary type periampullary cancer showing the highest stroma proportion (median 19% and 73% respectively). Adjusted Cox regression models revealed both positive (pancreato-biliary type periampullary cancer and oestrogen negative breast cancer, HR(95%CI)=0.56(0.34-0.92) and HR(95%CI)=0.41(0.17-0.98) respectively) and negative (intestinal type periampullary cancer, HR(95%CI)=3.59(1.49-8.62)) associations of the tumour stroma fraction with survival.

Interpretation

Our study provides an objective quantification of the tumour stroma fraction across major types of solid cancer. Findings strongly argue against the commonly promoted view of a general associations between high stroma abundance and poor prognosis. The results also suggest that full exploitation of the prognostic potential of tumour stroma requires analyses that go beyond determination of stroma abundance.

Funding

The Swedish Cancer Society, The Lions Cancer Foundation Uppsala, The Swedish Government Grant for Clinical Research, The Mrs Berta Kamprad Foundation, Sweden, Sellanders foundation, P.O.Zetterling Foundation, and The Sjöberg Foundation, Sweden.

Stromal Protein-Mediated Immune Regulation in Digestive Cancers

Simple Summary

Solid cancers are surrounded by a network of non-cancerous cells comprising different cell types, including fibroblasts, and acellular protein structures. This entire network is called the tumor microenvironment (TME) and it provides a physical barrier to the tumor shielding it from infiltrating immune cells, such as lymphocytes, or therapeutic agents. In addition, the TME has been shown to dampen efficient immune responses of infiltrated immune cells, which are key in eliminating cancer cells from the organism. In this review, we will discuss how TME proteins in particular are involved in this dampening effect, known as immunosuppression. We will focus on three different types of digestive cancers: pancreatic cancer, colorectal cancer, and gastric cancer. Moreover, we will discuss current therapeutic approaches using TME proteins as targets to reverse their immunosuppressive effects.

Abstract

The stromal tumor microenvironment (TME) consists of immune cells, vascular and neural structures, cancer-associated fibroblasts (CAFs), as well as extracellular matrix (ECM), and favors immune escape mechanisms promoting the initiation and progression of digestive cancers. Numerous ECM proteins released by stromal and tumor cells are crucial in providing physical rigidity to the TME, though they are also key regulators of the immune response against cancer cells by interacting directly with immune cells or engaging with immune regulatory molecules. Here, we discuss current knowledge of stromal proteins in digestive cancers including pancreatic cancer, colorectal cancer, and gastric cancer, focusing on their functions in inhibiting tumor immunity and enabling drug resistance. Moreover, we will discuss the implication of stromal proteins as therapeutic targets to unleash efficient immunotherapy-based treatments.

Prognostic value of preoperative glucose to lymphocyte ratio in patients with resected pancreatic cancer

BACKGROUND

Inflammatory factors and fasting blood glucose were verified to be associated with the prognosis of pancreatic ductal adenocarcinoma. The goal of this study is to confirm the prognostic role of preoperative blood glucose to lymphocyte ratio for patients with resected pancreatic ductal adenocarcinoma.

METHODS

A total of 259 pancreatic ductal adenocarcinoma patients were enrolled and randomly divided into training cohort and validation cohort. The training cohort was used to generate an optimal cutoff value and the validation cohort was used to further validate the model.

RESULTS

A total of 259 patients were incorporated in this study and randomly divided into the training cohort (n = 130, 1/2 of 259) and the validation cohort (129, 1/2 of 259). The optimal cutoff value of glucose to lymphocyte ratio was calculated to be 3.47 for overall survival. Cox regression analysis found that preoperative blood glucose to lymphocyte ratio was independent risk factor (p = 0.040) for overall survival. Prognostic values of glucose to lymphocyte ratio on overall survival were observed in younger male patients with pancreatic body and tail cancer, American Joint Committee on Cancer 8th N1 stage, without microvascular and peripancreatic fat invasion, and Carbohydrate antigen 19-9 higher than 200 U/ml. A prognostic prediction model of overall survival was designed and presented in nomogram.

CONCLUSION

Preoperative blood glucose to lymphocyte ratio is an independent biomarker to predict the overall survival for pancreatic ductal adenocarcinoma patients who underwent curative resection.

Topographical Distribution and Spatial Interactions of Innate and Semi-Innate Immune Cells in Pancreatic and Other Periampullary Adenocarcinoma

Background

The clinical management of pancreatic and other periampullary neoplasms remains challenging. In contrast to other cancer types, immunotherapies are largely ineffective, and the reason for the deprived immune response and the immune inhibiting cellular composition is only fragmentarily understood. The aim of this study was to comprehensively map the abundance, topographic distribution and spatial interaction of innate and innate-like immune cells in the tumor microenvironment of periampullary adenocarcinoma.

Methods

Multiplexed immunofluorescent imaging was performed on tissue microarrays with tumors from a consecutive cohort of 175 patients with resected periampullary adenocarcinoma. To obtain a detailed spatial analysis of immune cell infiltration, two multiplex immune panels including antibodies against CD3, NKp46, CD56, CD68, CD163 and CD1a, CD208, CD123, CD15, CD68 and pan-cytokeratin were applied.

Results

The infiltration of natural killer (NK) and NK-like T (NKT) cells was lower in malignant compared to benign tissue. NKT cells were more abundant in intestinal type compared to pancreatobiliary type tumors, and were associated with more favorable clinicopathological features and a prolonged survival. The interaction of NKp46⁺ NKT cells with macrophages was also associated with a prolonged survival.

Conclusions

This study provides a comprehensive map of the innate immune landscape in periampullary adenocarcinoma. NK cells, and even more so NKT cells, are revealed to be central players in the local immune response in a clinically relevant context.

CLEC3B as a potential diagnostic and prognostic biomarker in lung cancer and association with the immune microenvironment

Background

Lung cancer is the leading cause of cancer-related mortality globally. Discovering effective biomarkers for early diagnosis and prognosis is important to reduce the mortality rate and ensure efficient therapy for lung cancer patients. C-type lectin domain family 3 member B (CLEC3B) has been reported in various cancers, but its correlation with lung cancer remains elusive.

Methods

The GEO, TCGA and Oncomine databases were analyzed to examine the expression of CLEC3B in lung cancer. The CLEC3B mRNA levels in 15 patient tissue samples were detected by real-time PCR and the CLEC3B protein levels in 34 patient tissue samples were detected by immunohistochemistry. A Chi-square test was performed to analyze the correlation of CLEC3B expression and clinicopathological factors. The diagnostic value of CLEC3B was revealed by receiver operating characteristic (ROC) curves. Univariate and multivariate Cox proportional hazards regression models and Kaplan–Meier plots were used to evaluate the prognostic value of CLEC3B in lung cancer. The TIMER database was used to evaluate the correlation of CLEC3B and immune infiltration. Gene set enrichment analysis revealed tumor‐associated biological processes related to CLEC3B.

Results

CLEC3B is significantly downregulated in lung cancer patients compared with nontumor controls according to database analysis and patient tissue sample detection (p < 0.001). Specifically, CLEC3B is significantly downregulated in stage IA lung cancer patients (p < 0.001) and has a high diagnostic accuracy (area under the receiver operating characteristic curve > 0.9). Moreover, low expression of CLEC3B is related to poor progression-free survival (HR = 0.60, 95% CI 0.49–0.74, p = 8.3e−07) and overall survival (HR = 0.66, 95% CI 0.58–0.75, p = 2.1e−10), indicating it as a risk factor for lung cancer. Multivariate analysis value showed that low expression of CLEC3B may be an independent risk factor for disease‐free survival in lung cancer patients (HR = 0.655, 95% CI 0.430–0.996, Cox p = 0.048). In addition, we also investigated the potential role of CLEC3B in tumor-immune interactions and found that CLEC3B might be associated with the immune infiltration and immune activation of lung cancer, especially in squamous cell carcinoma.

Conclusions

Our findings indicate that CLEC3B expression is downregulated in lung cancer and reveal the diagnostic and prognostic potential of CLEC3B in lung cancer and its potential as an immune-related therapeutic target in lung cancer.