Nakkireddy, Sumanth Reddy

Abstract Background: Understanding the spatiotemporal dynamics of the tumor microenvironment (TME) during treatment is crucial for improving outcomes in pancreatic ductal adenocarcinoma (PDAC). This study presents a pre-planned spatiotemporal analysis from a phase II trial evaluating neoadjuvant mFOLFIRINOX in resectable PDAC, aimed at elucidating treatment resistance mechanisms and identifying predictive biomarkers. Methods: This multi-institutional, single-arm phase II trial enrolled treatment-naïve patients with histologically confirmed resectable PDAC (NCT05066802) who underwent six cycles of neoadjuvant mFOLFIRINOX followed by surgery. The primary endpoint was the R0 resection rate. Matched pre-treatment tumor biopsies and post-treatment surgical specimens were analyzed from selected responders (R) and non-responders (NR) using high-resolution spatial transcriptomics (10X Genomics Visium). This approach enabled comprehensive mapping of TME evolution and treatment-induced changes. Results: Between May 2020 and February 2022, 27 patients were enrolled. Of these, 24 patients underwent surgery after neoadjuvant treatment (resection rate: 88.9%), and 23 achieved R0 resection (85.2%), meeting the primary endpoint. However, 20 resected patients (83.3%) recurred. Spatiotemporal analysis was performed on matched samples from four R and four NR. In pre-treatment samples from NR group, clusters dominated by fibroblasts (enriched for ECM, EMT, and CAF phenotypes) and B cell subtypes were distinct compared to R group, which exhibited elevated antibody responses and robust T cell and macrophage activity. NR samples showed increased IgM expression but lacked corresponding Fc receptors, suggesting limited immune effector function. Post-treatment NR surgical specimens displayed elevated expression of MUC5AC and CEACAM6 in malignant ductal cells, alongside genes promoting tumor growth. Neoadjuvant mFOLFIRINOX increased fibroblast, while reducing immune cell components. NR tumors formed distinct "resistance niches" after chemotherapy, marked by ECM production, CAF activation, and RTK signaling. Conversely, R tumors showed enhanced epithelial organization with increased type-I hemidesmosome complexes facilitating epithelial adhesion to the basement membrane. Conclusion: Three months of neoadjuvant mFOLFIRINOX was feasible and tolerable for resectable PDAC but associated with high recurrence rate, emphasizing the need for better patient selection. Spatiotemporal analysis revealed critical differences in TME dynamics according to treatment response, including the emergence of resistance niches and distinct pre-treatment spatial signatures. These findings support the integration of spatial biology into future pancreatic cancer trials to improve stratification and therapy. Citation Format: Choong-kun Lee, Changjin Hong, Ho Kyoung Hwang, Min Hwan Kim, Jean Rene Clemenceau, Seock-Jin Chung, Isabel Barnfather, Sumanth Reddy Nakkireddy, Jeong Hwan Park, Ji-Youn Sung, Minji Kim, Yooho Choi, Inyeop Jang, Sunho Park, Sung Hyun Kim, Chang Moo Kang, Hye Jin Choi, Tae Hyun Hwang. Spatiotemporal evolution of the tumor microenvironment predicts treatment response in phase II trial of neoadjuvant chemotherapy for resectable pancreatic cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 771.

Abstract Patients with head and neck squamous cell carcinomas (HNSCC) often have poor outcomes due to suboptimal risk management and treatment strategies; yet integrating novel prognostic biomarkers into clinical practice is challenging. Here, we report the presence of multinucleated giant cells (MGC)—a type of macrophages—in tumors from patients with HNSCC, which are associated with a favorable prognosis in treatment-naive and preoperative chemotherapy–treated patients. Importantly, MGC density increased in tumors following preoperative therapy, suggesting a role of these cells in the antitumoral response. To enable clinical translation of MGC density as a prognostic marker, we developed a deep-learning model to automate its quantification on routinely stained pathological whole slide images. Finally, we used spatial transcriptomic and proteomic approaches to describe the MGC-related tumor microenvironment and observed an increase in central memory CD4 T cells. We defined an MGC-specific signature resembling to TREM2-expressing mononuclear tumor-associated macrophages, which colocalized in keratin tumor niches. Significance: Novel individual biomarkers are needed to guide therapeutic decisions for patients with head and neck cancer. We report for the first time, granulomas of TREM2-expressing multinucleated giant macrophages in keratin-rich tumor niches, as a biomarker of favorable prognosis and developed a deep-learning model to automate its quantification on routinely stained pathological slides.

Abstract Background: Analyzing interactions and the microanatomical distribution of the tumor immune microenvironment (TIME) is essential for a comprehensive understanding of tumor progression. Imaging Mass Cytometry (IMC) is a high-dimensional tissue imaging system that enables the deeper and multiparametric in situ exploration of tumor microenvironments at a single-cell level. This study describes a pipeline for analyzing 39-plex IMC images of human HPV(+) oropharyngeal squamous cell carcinoma (OPSCC) tumor sections, elucidating immune cells' nature, functions, and interactions with tumor cells. The IMC data analysis provides valuable information for clinical studies that could be used for the identification of prognostic biomarkers and mechanisms of resistance to current immunotherapies. Methods: A retrospective study was conducted on (N=20) patients (10 progressors and 10 non-progressors) with HPV(+) OPSCC. The dataset comprises paired hematoxylin and eosin (H&E) images with IMC data from multiple Regions of Interest (ROIs) (N=82), enabling the in-situ profiling of 39 proteins in primary tumor tissue. Cell segmentation was followed by employing original images and their segmentation masks for cell phenotype identification using a probabilistic neural network method to spatially profile cancer cells, macrophages, NK-M cells, B-cells, and distinct T-cell populations. We identified protein expression at single-cell resolution and investigated the impact of multicellular spatial organization on the disease progression. The spatial neighborhood interaction score of tumor cells with various other cells was computed for all the ROIs across the two groups. Results: The markers showed varying expression levels of intensity distribution, likely due to differences in tissue sample composition. Specifically, CD3, CD4, CD8, CD19, CD20, and CD68 showed significantly higher intensity in non-progressors. Interestingly, FoxP3 also had a significantly higher intensity level in non-progressors. Cell phenotyping emphasized that non-progressors generally have increased levels of NK-M cell markers. Spatial analysis revealed a significant increase in interaction between CD3 cells and tumor cells in non-progressors (p-value=0.07, Mann-Whitney U test). Furthermore, a significant rise in the interaction between B-cells and FoxP3+ cells was observed in non-progressors (p-value<0.05, Mann-Whitney U test) suggesting a coordinated suppression of tumor growth. Conclusion: We developed a deep learning-based pipeline to analyze spatial relationships in the TIME of IMC data, demonstrating its clinical relevance in predicting disease progression. Our findings emphasize the importance of spatial relationships in the TME for response and suggest that proximity between either CD3+ T-cells to cancer cells or B-cells to FoxP3+ cells are candidate biomarkers. Citation Format: Sumanth Reddy Nakkireddy, Routman David, Kathleen Bartemes, Daniel Ma, Tae Hyun Hwang, Kathryn Van Abel, Chadi Nimeh Abdel-Halim. AI-enhanced spatial proteomics reveals B-cells & FoxP3+ cells interactions in HPV(+) oropharyngeal squamous cell carcinoma Progression [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 5488.

Abstract Traditional cancer research methodologies predominantly focus on 2D tumor modeling, including spatial profiling and AI-driven analysis at the 2D level. To advance beyond these limitations, we implemented Holotomography (HT), a 3D, label-free optical imaging method. This technique, leveraging refractive index contrasts, allows for enhanced visualization of tumor tissues, organoids, and cells in three dimensions. Using serial H&E sections and unstained slides, we created detailed 3D tumor models in gastric cancer, providing insights into tumor, immune, and stromal composition, as well as cell-cell interactions and extracellular matrix characteristics. AI algorithms further enhanced our 3D tumor modeling, offering novel perspectives on the initiation and progression of tumor, immune, and stromal components at the subcellular level. Analyzing gastric tumors treated with immune checkpoint inhibitors before, during, and after treatment revealed dynamic changes in spatial composition over time, offering unique insights into the evolution of the 3D tumor ecosystem. Additionally, by integrating spatial transcriptomic and in situ molecular imaging data, we developed a novel Multimodal 3D Tumor Model. This model not only uncovers subcellular-level molecular characteristics but also elucidates cellular interactions during treatment, suggesting potential biomarkers and therapeutic approaches. Lastly, by constructing 3D cellular interaction networks incorporating true 3D interactions, we provided an in-depth view of cellular communication within the tumor microenvironment. This pioneering approach, representing the first Multimodal 3D Tumor Model at the subcellular level, offers novel insights into gastric cancer's spatial evolution, potentially revolutionizing our understanding and treatment of gastric cancer. Citation Format: Yoon Ho Choi, Minji Kim, Katya Redekop, Inyeop Jang, Jean R. Clemenceau, Isabel Barnfather, Sumanth Reddy Nakkireddy, Changjin Hong, Sunho Park, Ji-Youn Sung, Jeong Hwan Park, Shitara Kohei, Aoki Yu, Sung Hak Lee, Sam Wang, Jaeho Cheong, Tae Hyun Hwang. AI-driven multimodal 3D tumor modeling with spatial molecular insights for predicting ICI response in gastric cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7400.

Abstract Although the majority of HPV (+) oropharyngeal squamous cell carcinoma (OPSCC) patients exhibit a favorable prognosis, approximately 20% face recurrent or metastatic disease, posing management challenges. Therefore, the identification of robust prognostic markers for risk stratification is essential. Our study focused on the comprehensive characterization of intratumor heterogeneity within the tumor immune microenvironment (TME) in both primary tumors and metastatic lymph nodes. Utilizing computational approaches, we integrated hematoxylin and eosin (H&E) and 11 adjacent immunohistochemistry (IHC)-stained slides to investigate the cellular composition and functional characteristics across different regions within the TME. Based on these detailed immune characteristics, we classified patients into specific immune subtypes. Our integrative TME analysis indicated that patients with high tumor-infiltrating lymphocytes (TIL), increased CD8 + levels, and reduced CD163 + cell counts within their primary tumors were likely to have a more favorable progression-free survival outcome. The integration of H&E and IHC data enhances the functional characterization of immune phenotypes of the TME with biological interpretability and improves patient stratification in HPV(+) OPSCC.

Abstract Background: Recent studies support that high levels of tumor-infiltrating lymphocytes (TILs) were associated with better prognosis in HPV(+) oropharyngeal squamous cell carcinoma (OPSCC). While these studies provide a potential role for TILs as a prognostic biomarker, the analyses were relied on manual quantification performed by pathologists, resulting in inter-observer variability. Deep Learning (DL)-based whole slide hematoxylin and eosin (H&E) image analyses may overcome these challenges. However, this analysis only provides the presence of cell types without consideration for the functional roles of each cell within the tumor immune microenvironment (TME). In this work, we develop a computational pipeline to integrate paired H&E and immunohistochemistry (IHC) images to functionally characterize TILs and investigate their prognostic utility. Methods: We analyzed 88 patients: 42 stage I, 39 stage II, and 7 stage III. Our data contain both H&E and IHCs examining FoxP3, CD3, PD-L1, CD20, etc. on serial sections of the tissue. In-house DL-based H&E analysis used to identify TILs, tumors, and stroma in each tumor, then performed registration between adjacent H&E and IHC images from the same tissue. The patients were then classified into three basic immune phenotypes: immune inflamed (IN; high TILs in the tumor region), immune excluded (EX; TILs are mostly localized in stroma), and immune desert (ID; few/no TILs) based on TIL enrichment in the TME. To functionally characterize TILs, we quantified protein expression from the adjacent IHCs. For example, we further classified each patient into different subtypes based on enriched protein expression (e.g., FoxP3 high IN, FoxP3 low IN). We used the Kaplan-Meier method and multivariate Cox analysis to evaluate the prognostic value of different subtypes enriched with different proteins to predict disease-free survival (DFS). Results: The IN group with 43 patients was significantly associated with good prognosis. Interestingly, further stratification of the IN subgroup based on Foxp3 quantification on TIL regions (i.e., high FoxP3 IN and low FoxP3 IN) showed that high protein expression of FoxP3 in TILs in the IN subgroup is significantly associated with a better prognosis compared to other immune subgroups (HR, 0.16; p-value, 0.003). Multivariate analysis, including other clinical covariates showed that the immune subtypes associated with high FoxP3 are independently associated with DFS. These results demonstrate that DL-based integrative IHC and H&E image analysis could be used to identify subgroups with distinct clinical outcomes. Furthermore, our results reveal unknown roles for Foxp3 expression in the TILs in HPV(+) OPSCC as a prognostic biomarker, a finding which should be evaluated in a larger cohort. Citation Format: Sumanth Reddy Nakkireddy, Inyeop Jang, Minji Kim, Linda X. Yin, Michael Rivera, Joaquin J. Garcia MD, Kathleen R. Bartemes, David M. Routman, Eric J. Moore, Daniel J. Ma, Kathryn M. Van Abel, Tae Hyun Hwang. Integrative spatial analysis of paired IHC and H&E images identifies Foxp3 enriched tumor-infiltrating lymphocytes associated with disease-free survival in human papillomavirus (HPV)-related oropharyngeal squamous cell carcinoma. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 4638.