Shi, Ji

Background Tumor-infiltrating lymphocytes (TILs) are established determinants of clinical outcomes in breast cancer. 1–6 While stromal TILs (sTILs) are commonly evaluated due to easier scoring and higher interobserver reproducibility, intraepithelial TILs may provide additional insights into anti-tumor immunity. 7 We aimed to comprehensively characterize immunotype-associated cellular landscapes in triple-negative breast cancer (TNBC) and identify features associated with clinical outcomes. Methods 212 TNBC tumors were classified into three immunotypes based on immune infiltrate location and density: immune desert (ID), immune excluded (IE), and immune inflamed (IN). 8 Single-cell spatial transcriptomics was performed on 75 tumors to define the cellular composition and functional states of the tumor microenvironment across immunotypes. Findings were validated and integrated with bulk RNA sequencing data and clinical outcome datasets from FinXX, CALGB-40603, I-SPY 2, and Real-World Clinico-Genomic Data (RWCGD) cohorts. Results Patients with IN tumors had significantly improved outcomes compared with ID tumors. Despite high sTILs, IE tumors showed poor outcomes similar to ID tumors. Single-cell spatial analysis revealed that ID and IE tumors exhibited reduced major histocompatibility complex I/II expression and fewer tumor-resident plasmacytoid dendritic cells (pDCs) compared with IN tumors, which were enriched for interferon-alpha (IFNα) and interferon-gamma (IFNγ) responses. High IFN response scores were associated with favorable outcomes across multiple therapy types in independent datasets. Deconvolution of RWCGD bulk RNA-seq data confirmed that pDC abundance correlated with improved outcomes specifically in hormone receptor–negative subtypes. Conclusions Our study highlights pDCs and IFN signaling as hallmarks of effective anti-tumor immunity in TNBC. Immunotype-based profiling underscores the limited prognostic value of sTILs alone in immune-excluded tumors and supports pDCs and IFN pathways as potential biomarkers for prognosis and therapeutic development in TNBC.

Abstract Background: Most colorectal cancer (CRC) arises from polyps and is mainly prevented by polypectomy. The most important polyps to manage with colonoscopy are those with highest CRC risk- namely, advanced polyps (> 1cm, villous histology or high-grade dysplasia (HGD). Yet, 48% of advanced polyps recur within 1 to 3 years of removal, and up to 5% of advanced polyps under surveillance still progress to CRC. We performed this pilot study to investigate molecular and microenvironment heterogeneities of polyps and how they might impact their clinical behavior. Methods: GeoMx spatial transcriptomics was performed on FFPE tissues from three different polyp outcome phenotypes (POPs) including the polyp that does not recur (POP-NR), that recurs following polypectomy but cured by colonoscopy (POP-R) or the polyp despite polypectomy develops CRC at the polypectomy(ies) site (POP-CRC). Normal colon and polyp with low or HGD from the index polyp were assessed from 6 patients with POP-NR, 9 with POP-R and 12 with POP-CRC with a minimum of 2 follow up colonoscopies at 3-year intervals. Epithelium was identified as PanCK positive, and stroma identified as PanCK negative and positive nuclear staining. Cell type deconvolution was implemented using single-cell adult human intestinal tract catalogue (Elmentaite et al., https://www.gutcellatlas.org/). Differential gene expression (DEG), cell type abundance and functional enrichment analyses was performed using linear mixed effects model, and observations with p-value < 0.05 and log2FC >= |1| are reported. Results: The greatest number of DEGs were identified in stroma of the index POP-CRC compared to POP-NR polyps (n = 17 down and 11 upregulated genes), followed by epithelium of the index POP-CRC vs POP-NR polyps (n = 10 down and one upregulated gene(s)). Between the index POP-CRC and POP-R, epithelium showed downregulation of 7 genes, and upregulation of no genes. In stroma two genes were down and none upregulated. Four genes were downregulated in stroma of the index POP-R vs POP-NR polyp and 3 in epithelium, and one gene was upregulated in stroma and another in epithelium. MZT2B was upregulated in stroma of both the index POP-R (log2FC = 1.09, p-value = 0.0014) and POP-CRC (log2FC = 1.29, p-value = 0.0001) and has been implicated as a prognostic marker associated with worse prognosis in certain cancers. IgM and IgA plasma cells, proximal progenitors, MMP9+ inflammatory macrophages and myofibroblasts were found to be downregulated in epithelium of index POP-CRC and POP-R compared to POP-NR polyps; while microfold cells were downregulated in stroma of index POP-CRC and POP-R compared to POP-NR polyps. Conclusions: The stromal microenvironment and less so, epithelial features present in the index polyp differ based on a polyp’s future clinical behavior. The polyp-immune interaction warrants further study as a potential prevention target against polyp progression. Citation Format: Mrunal Dehankar, Lisa A. Boardman, Daniel O'Brien, E. Aubrey Thompson, Jennifer M. Kachergus, Ji Shi, Alexej Abyzov, Milovan Suvakov, Rondell P. Graham, Chen Wang. Stromal immune composition significantly contributes to adenomatous polyp recurrence or progression to cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2026; Part 1 (Regular Abstracts); 2026 Apr 17-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2026;86(7 Suppl):Abstract nr 4018.

Background: Discovering risk biomarkers in small benign breast disease (BBD) biopsies is constrained by scarce tissue and microanatomic heterogeneity of terminal duct lobular units (TDLUs). We tested whether tissue-sparing tissue microarray (TMA)–based Digital Spatial Profiling (DSP) can deliver reproducible, biologically coherent protein measurements across benign lobules and breast cancers (BCs), and how well DSP aligns with standard immunoassays. Methods: We performed a pilot using tissues from the Mayo Clinic BBD cohort using TMAs representing four contexts: terminal duct lobular units (TDLUs) from BBD biopsies preceding BC and matched BBD-controls, subsequent BCs, and BC-associated TDLUs. We profiled 79 proteins by DSP (37 retained after QC) and benchmarked against chromogenic IHC and OPAL immunofluorescence. Reproducibility was evaluated using intraclass correlation coefficients (ICCs), cross-platform agreement (weighted kappa), marker correlations, and mixed-effects models with false-discovery-rate (FDR) control. Results: We analyzed 368 BBD-TDLU cores (88 cases; 88 controls), 204 BC cores and 110 BC-associated TDLU cores. ICCs were highest in BC tissues, and lower in BC-associated TDLUs and BBD-TDLUs. Agreement was slight–to-fair in TDLUs but moderate (ER/PR) to substantial (BCL2) in BC. DSP recapitulated expected immunologic correlations (CD45 with T-cell, B-cell, and macrophage markers) and tissue-type gradients (BC > BC-associated TDLUs > BBD-TDLUs). Exploratory case–control differences in BBD-TDLUs did not persist after FDR control. Conclusions: TMA-based DSP is feasible in archival breast tissues and yields biologically coherent, cross-platform-benchmarked profiles that are particularly robust in BC, while conserving scarce TDLUS and clarifying current limits of single-marker risk stratification from benign lobules. These data provide a foundation for refined sampling and expanded panels in future TDLU-focused studies.

Abstract Parkinson’s disease (PD) is the fastest growing neurodegenerative disease in the world1. Gastrointestinal (GI) dysfunction can occur decades before motor impairments and in up to 80% of individuals living with PD2,3,4. We investigated peripheral relationships that may underlie environmental contributions to mechanisms along the gut-blood axis that contribute to PD pathogenesis. Single-cell multiomic spatial molecular imaging (SMI) of colonic tissue localized inflammatory injury within epithelial cells that appear to be associated with iron mishandling in both inflammatory bowel disease (IBD) and PD biosamples. We found that both the single-cell SMI of RNA and protein revealed parallel cross-modal dysregulation in the gut epithelium, in both IBD and PD biosamples. These data are accompanied by plasma (PD) and stool (IBD) protein depletion of CCL22. Our findings suggest iron mishandling along the gut barrier likely contributes to systemic inflammation, which may be the catalyst that primes circulating immune cells to body-first PD pathogenesis.

Abstract Parkinson’s disease (PD) is the fastest growing neurodegenerative disease in the world 1 . Gastrointestinal (GI) dysfunction can occur decades before motor impairments and in up to 80% of individuals living with PD 2,3,4 . We investigated peripheral relationships that may underlie mechanisms along the gut-blood axis that contribute to PD pathogenesis. Single-cell multiomic spatial molecular imaging (SMI) of colonic tissue localized inflammatory injury within epithelial cells that appear to be associated with iron mishandling in both inflammatory bowel disease (IBD) and PD biosamples. We found that both the single-cell SMI of RNA and protein revealed parallel cross-modal dysregulation in the gut epithelium, in both IBD and PD biosamples. These data are accompanied by plasma (PD) and stool (IBD) protein depletion of CCL22. Our findings suggest iron mishandling along the gut barrier likely contributes to systemic inflammation, which may be the catalyst that primes circulating immune cells to body-first PD pathogenesis.

Abstract Brain white matter tracts undergo structural and functional changes linked to late-life cognitive decline, but the cellular and molecular contributions to their selective vulnerability are not well defined. In naturally aged mice, we demonstrate that senescent and disease-associated microglia (DAM) phenotypes converge in hippocampus-adjacent white matter. Through gold-standard gene expression and immunolabeling combined with high-dimensional spatial mapping, we identified microglial cell fates in aged white matter characterized by aberrant morphology, microenvironment reorganization, and expression of senescence and DAM markers, including galectin 3 (GAL3/Lgals3), B-cell lymphoma 2 (Bcl2), and cyclin dependent kinase inhibitors, including Cdkn2a/p16ink4a. Pharmacogenetic or pharmacological targeting of p16ink4a or BCL2 reduced white matter GAL3+ DAM abundance and rejuvenated microglial fimbria organization. Our results demonstrate dynamic changes in microglial identity in aged white matter that can be reverted by senotherapeutic intervention to promote homeostatic maintenance in the aged brain.

Abstract Cerebral malaria is the deadliest complication that can arise from Plasmodium infection. CD8 T-cell engagement of brain vasculature is a putative mechanism of neuropathology in cerebral malaria. To define contributions of brain endothelial cell major histocompatibility complex (MHC) class I antigen-presentation to CD8 T cells in establishing cerebral malaria pathology, we developed novel H-2Kb LoxP and H-2Db LoxP mice crossed with Cdh5-Cre mice to achieve targeted deletion of discrete class I molecules, specifically from brain endothelium. This strategy allowed us to avoid off-target effects on iron homeostasis and class I-like molecules, which are known to perturb Plasmodium infection. This is the first endothelial-specific ablation of individual class-I molecules enabling us to interrogate these molecular interactions. In these studies, we interrogated human and mouse transcriptomics data to compare antigen presentation capacity during cerebral malaria. Using the Plasmodium berghei ANKA model of experimental cerebral malaria (ECM), we observed that H-2Kb and H-2Db class I molecules regulate distinct patterns of disease onset, CD8 T-cell infiltration, targeted cell death and regional blood–brain barrier disruption. Strikingly, ablation of either molecule from brain endothelial cells resulted in reduced CD8 T-cell activation, attenuated T-cell interaction with brain vasculature, lessened targeted cell death, preserved blood–brain barrier integrity and prevention of ECM and the death of the animal. We were able to show that these events were brain-specific through the use of parabiosis and created the novel technique of dual small animal MRI to simultaneously scan conjoined parabionts during infection. These data demonstrate that interactions of CD8 T cells with discrete MHC class I molecules on brain endothelium differentially regulate development of ECM neuropathology. Therefore, targeting MHC class I interactions therapeutically may hold potential for treatment of cases of severe malaria.

Abstract Purpose: Young women’s breast cancers (YBCs), defined herein as breast cancer (BC) diagnosed at age ≤ 40 years, are etiologically and biologically distinct from postmenopausal BCs, demonstrating a higher proportion of estrogen receptor (ER) negative subtypes and a worse prognosis. Studies have not assessed whether risk factors such as recent parity and body mass index (BMI) are linked to protein expression profiles in YBCs. Accordingly, we examined protein expression assessed with GeoMx® Digital Spatial Profiling (DSP) in 640 YBCs prepared as tissue microarrays (TMAs) in the Young Women’s BC Study (YWBCS). Methods: The YWBCS is a multi-institutional cohort of 1,297 women with incident BC diagnosed at age <40 years identified through medical record review, with risk factor data collected via questionnaires. TMAs including 1,505 cores of 640 YBCs that were collected prior to therapy and were suitable for molecular analysis were evaluated for 72 markers with NanoString GeoMx DSP. Associations of log-transformed biomarker expression with BMI and time from last live birth to BC diagnosis were carried out using linear mixed modeling approaches, accounting for multiple tumor tissue cores per individual. Results: Mean age at BC was 35.5 years (range 17-40). Tumor subtypes based on clinical immunostaining were the following: luminal B 46%; luminal A 36%; triple negative 12% and HER2+ 6%. Of 72 biomarkers tested, 53 (74%) were evaluable after QC and normalization. GeoMx ER and PR protein levels were higher in luminal A and luminal B YBCs than in HER2+ and TNBC (p=2x10-94 and 2x10-52, respectively), supporting the accuracy of the GeoMx assay. Compared to nulliparas or those with most recent birth > 3 years before BC, those who gave birth ≤ 3 years prior to BC had tumors with higher expression of PR (8x10-6) and ER (p=0.006), and lower expression of BCL-2 (p=0.006). In contrast, women with last birth ≥ 6 years prior to BC had higher expression of CD8 (p=0.006) and PD-1 (p=0.001) than nulliparas and those with more recent births. Compared to women with BMI between 18.5 and 29.9, those with higher or lower BMIs had higher expression of CD34 (p=0.0003), CD4 (p=0.0009) and CD14 (p=0.001). BCL-2 expression was higher among women with BMI < 18.5 compared to those with higher values (p=0.005). Women with BMI of 30+ had higher β2-microglobulin levels than those with lower BMIs (p=0.003). Inverse dose-response associations (with higher BMI related to lower levels) were observed for expression of ER (p=0.01), pan cytokeratin (p=0.004) and four proteins in the p13K/AKT signaling pathway (p ≤ 0.01 for each). Conclusions: Preliminarily, we identify differences in hormonal and immune markers in YBCs in relation to timing of births and BMI, suggesting that these factors may contribute to etiologic heterogeneity. Multivariate modeling is ongoing to incorporate molecular subtype and other risk factors. Citation Format: Robert A. Vierkant, Laura M. Pacheco-Spann, Stacey J. Winham, Jennifer M. Kachergus, Ji Shi, Craig Snow, Fergus J. Couch, Janet E. Olson, Chen Wang, Derek C. Radisky, Amy C. Degnim, E. Aubrey Thompson, Laura C. Collins, Kathryn J. Ruddy, Ann H. Partridge, Mark E. Sherman. Associations of parity and body mass index with NanoString digital spatial protein profiling in early onset breast cancer [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 1000.

Abstract Background: Biopsy diagnoses of benign breast disease (BBD) confer a 1.5- to 4-fold increased risk of developing breast cancer (BC) compared with women without BBD. Previously, we reported that decreased numbers of specific immune cell types in lobules of BBD biopsies predicted increased BC risk, suggesting the promise of future tissue biomarker studies to define BC risk markers among BBD patients. Thus, we applied protein-based GeoMx® Digital Spatial Profiling (DSP) to BBD biopsies preceding BC (cases) and to BBD biopsies from cancer-free patients (controls) to identify possible BC risk markers, which we then evaluated in subsequent BC tissues and in surrounding normal lobules of cases. Methods: Archived pathology slides of BBD biopsies were reviewed and used to guide preparation of TMAs containing 1.0-mm diameter FFPE cores of lobules from an age- and cohort-period-matched set of 91 cases and 88 controls from the Mayo Clinic BBD Cohort. For patients who later developed BC, we prepared TMAs of BC tissue and surrounding mapped normal lobules. We applied GeoMx® DSP (immune and canonical signaling proteins) to both sets of TMAs. Following QC and data normalization, associations of case status with log-transformed biomarker expression in lobules of BBD biopsies were carried out using linear mixed modeling approaches, accounting for multiple ROIs per individual. Biomarkers significantly associated with case status (p<0.05) were further examined in BCs and adjacent normal lobules, using similar approaches. Results: The mean age at BBD biopsy was 52 years, and at BC diagnosis of cases, 61.4 years (mean time from BBD to BC was 10.2 years). A family history of BC was more frequent among cases (70% versus 43%; chi-square p=0.002). Of 72 biomarkers tested, 46 (64.4%) were evaluable after QC and normalization and 5 were associated with BC risk after adjustment for family history of BC: BCL2 (p=0.005), STING (p=0.006), CD44 (p=0.02), S100 protein (p=0.03) and pan-AKT (p=0.05); each showed higher levels in lobules of BBD biopsies of controls than cases. Three unique patterns appeared when examining these biomarkers across tissue type within cases: for BCL2 (p=5 x 10-9) and STING (p=2 x 10-19), levels were high in both BC and lobules surrounding BC but low in preceding BBD; for GAPDH (p=4 x 10-53) and pan-AKT (p=2 x 10-33), levels were high in BC, low in preceding BBD, and moderate in lobules surrounding BC; and for CD44 (p=2 x 10-6) and S100B (p=2 x 10-49), levels were low in BC, high in lobules surrounding BC and moderate in preceding BBD. Conclusions: Using a novel TMA of lobules in combination with DSP, we preliminarily identified immune-based and PI3 kinase-related protein biomarkers in BBD biopsies associated with BC risk. In case-only analyses, these markers demonstrated complex differences between lobules in BBD biopsies, subsequent BCs and adjacent normal lobules. Citation Format: Robert Alan Vierkant, Jodi M. Carter, Stacey J. Winham, Chen Wang, Jennifer M. Kachergus, Ji Shi, Raymond M. Moore, Bryan M. McCauley, Laura M. Pacheco-Spann, E A. Thompson, Derek C. Radisky, Amy C. Degnim, Mark E. Sherman. Towards prediction of breast cancer risk in benign biopsies with high-plex GeoMx spatial protein profiling [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 2203.