Analysis of Immune Intratumor Heterogeneity Highlights Immunoregulatory and Coinhibitory Lymphocytes as Hallmarks of Recurrence in Stage I Non-Small Cell Lung Cancer.

Our understanding of the molecular mechanisms underlying postsurgical recurrence of non-small cell lung cancer (NSCLC) is rudimentary. Molecular and T cell repertoire intratumor heterogeneity (ITH) have been reported to be associated with postsurgical relapse; however, how ITH at the cellular level impacts survival is largely unknown. Here we report the analysis of 2880 multispectral images representing 14.2% to 27% of tumor areas from 33 patients with stage I NSCLC, including 17 cases (relapsed within 3 years after surgery) and 16 controls (without recurrence ≥5 years after surgery) using multiplex immunofluorescence. Spatial analysis was conducted to quantify the minimum distance between different cell types and immune cell infiltration around malignant cells. Immune ITH was defined as the variance of immune cells from 3 intratumor regions. We found that tumors from patients having relapsed display different immune biology compared with nonrecurrent tumors, with a higher percentage of tumor cells and macrophages expressing PD-L1 (P =.031 and P =.024, respectively), along with an increase in regulatory T cells (Treg) (P =.018), antigen-experienced T cells (P =.025), and effector-memory T cells (P =.041). Spatial analysis revealed that a higher level of infiltration of PD-L1+ macrophages (CD68+PD-L1+) or antigen-experienced cytotoxic T cells (CD3+CD8+PD-1+) in the tumor was associated with poor overall survival (P =.021 and P =.006, respectively). A higher degree of Treg ITH was associated with inferior recurrence-free survival regardless of tumor mutational burden (P =.022), neoantigen burden (P =.021), genomic ITH (P =.012) and T cell repertoire ITH (P =.001). Using multiregion multiplex immunofluorescence, we characterized ITH at the immune cell level along with whole exome and T cell repertoire sequencing from the same tumor regions. This approach highlights the role of immunoregulatory and coinhibitory signals as well as their spatial distribution and ITH that define the hallmarks of tumor relapse of stage I NSCLC.

Towards the characterization of the tumor microenvironment through dictionary learning-based interpretable classification of multiplexed immunofluorescence images.

Objective.Histology image analysis is a crucial diagnostic step in staging and treatment planning, especially for cancerous lesions. With the increasing adoption of computational methods for image analysis, significant strides are being made to improve the performance metrics of image segmentation and classification frameworks. However, many developed frameworks effectively function as black boxes, granting minimal context to the decision-making process. Thus, there is a need to develop methods that offer reasonable discriminatory power and a biologically-informed intuition to the decision-making process.Approach.In this study, we utilized and modified a discriminative feature-based dictionary learning (DFDL) paradigm to generate a classification framework that allows for discrimination between two distinct clinical histologies. This framework allows us (i) to discriminate between 2 clinically distinct diseases or histologies and (ii) provides interpretable group-specific representative dictionary image patches, or 'atoms', generated during classifier training. This implementation is performed on multiplexed immunofluorescence images from two separate patient cohorts- a pancreatic cohort consisting of cancerous and non-cancerous tissues and a metastatic non-small cell lung cancer (mNSCLC) cohort of responders and non-responders to an immunotherapeutic treatment regimen. The analysis was done at both the image-level and subject-level. Five cell types were selected, namely, epithelial cells, cytotoxic lymphocytes, antigen presenting cells, HelperT cells, and T-regulatory cells, as our phenotypes of interest.Results.We showed that DFDL had significant discriminant capabilities for both the pancreatic pathologies cohort (subject-level AUC-0.8878) and the mNSCLC immunotherapy response cohort (subject-level AUC-0.7221). The secondary analysis also showed that more than 50% of the obtained dictionary atoms from the classifier contained biologically relevant information.Significance.Our method shows that the generated dictionary features can help distinguish patients presenting two different histologies with strong sensitivity and specificity metrics. These features allow for an additional layer of model interpretability, a highly desirable element in clinical applications for identifying novel biological phenomena.

GaWRDenMap: a quantitative framework to study the local variation in cell-cell interactions in pancreatic disease subtypes.

Spatial pattern modelling concepts are being increasingly used in capturing disease heterogeneity. Quantification of heterogeneity in the tumor microenvironment is extremely important in pancreatic ductal adenocarcinoma (PDAC), which has been shown to co-occur with other pancreatic diseases and neoplasms with certain attributes that make visual discrimination difficult. In this paper, we propose the GaWRDenMap framework, that utilizes the concepts of geographically weighted regression (GWR) and a density function-based classification model, and apply it to a cohort of multiplex immunofluorescence images from patients belonging to six different pancreatic diseases. We used an internal cohort of 228 patients comprised of 34 Chronic Pancreatitis (CP), 71 PDAC, 70 intraductal papillary mucinous neoplasm (IPMN), 16 mucinous cystic neoplasm (MCN), 29 pancreatic intraductal neoplasia (PanIN) and 8 IPMN-associated PDAC patients. We utilized GWR to model the relationship between epithelial cells and immune cells on a spatial grid. The GWR model estimates were used to generate density signatures which were used in subsequent pairwise classification models to distinguish between any two pairs of disease groups. Image-level, as well as subject-level analysis, were performed. When applied to this dataset, our classification model showed significant discrimination ability in multiple pairwise comparisons, in comparison to commonly used abundance-based metrics, like the Morisita-Horn index. The model was able to best discriminate between CP and PDAC at both the subject- and image-levels. It was also able to reasonably discriminate between PDAC and IPMN. These results point to a potential difference in the spatial arrangement of epithelial and immune cells between CP, PDAC and IPMN, that could be of high diagnostic significance. Further validation on a more comprehensive dataset would be warranted.

Central nervous system immune interactome is a function of cancer lineage, tumor microenvironment, and STAT3 expression.

BACKGROUNDImmune cell profiling of primary and metastatic CNS tumors has been focused on the tumor, not the tumor microenvironment (TME), or has been analyzed via biopsies.METHODSEn bloc resections of gliomas (n = 10) and lung metastases (n = 10) were analyzed via tissue segmentation and high-dimension Opal 7-color multiplex imaging. Single-cell RNA analyses were used to infer immune cell functionality.RESULTSWithin gliomas, T cells were localized in the infiltrating edge and perivascular space of tumors, while residing mostly in the stroma of metastatic tumors. CD163+ macrophages were evident throughout the TME of metastatic tumors, whereas in gliomas, CD68+, CD11c+CD68+, and CD11c+CD68+CD163+ cell subtypes were commonly observed. In lung metastases, T cells interacted with CD163+ macrophages as dyads and clusters at the brain-tumor interface and within the tumor itself and as clusters within the necrotic core. In contrast, gliomas typically lacked dyad and cluster interactions, except for T cell CD68+ cell dyads within the tumor. Analysis of transcriptomic data in glioblastomas revealed that innate immune cells expressed both proinflammatory and immunosuppressive gene signatures.CONCLUSIONOur results show that immunosuppressive macrophages are abundant within the TME and that the immune cell interactome between cancer lineages is distinct. Further, these data provide information for evaluating the role of different immune cell populations in brain tumor growth and therapeutic responses.FUNDINGThis study was supported by the NIH (NS120547), a Developmental research project award (P50CA221747), ReMission Alliance, institutional funding from Northwestern University and the Lurie Comprehensive Cancer Center, and gifts from the Mosky family and Perry McKay. Performed in the Flow Cytometry & Cellular Imaging Core Facility at MD Anderson Cancer Center, this study received support in part from the NIH (CA016672) and the National Cancer Institute (NCI) Research Specialist award 1 (R50 CA243707). Additional support was provided by CCSG Bioinformatics Shared Resource 5 (P30 CA046592), a gift from Agilent Technologies, a Research Scholar Grant from the American Cancer Society (RSG-16-005-01), a Precision Health Investigator Award from University of Michigan (U-M) Precision Health, the NCI (R37-CA214955), startup institutional research funds from U-M, and a Biomedical Informatics & Data Science Training Grant (T32GM141746).