Sherman, Mark E.

Abstract High-grade serous ovarian cancer (HGSOC), the most prevalent and lethal ovarian cancer subtype, is typically diagnosed at an advanced stage due to the absence of effective early detection strategies. Mounting evidence indicates that HGSOC originates from the fallopian tube (FT), particularly from secretory epithelial cells in the fimbria. However, incomplete understanding of the cellular and molecular origins of HGSOC continues to hinder early interception efforts. To address this, we established a clinically and genetically annotated living FT organoid biobank at Mayo Clinic, representing over 200 patient donors. Using optimized enrichment protocols, we isolated FT epithelial stem/progenitor cells from tissue and Tao brushings, generating organoids independent of anatomical site or laterality that recapitulate native FT epithelial architecture and lineage differentiation. We applied integrated multi-omic analyses, including single-cell RNA sequencing, single-nucleus RNA/ATAC sequencing, and bulk RNA-seq, across fresh cells and organoids. Integration with public datasets produced the largest single-cell atlas of normal and high-risk FTs to date. This analysis defined high-specificity gene signatures for secretory and multiciliated lineages that outperform canonical lineage markers (e.g. PAX8, FOXJ1), identified lineage-defining transcription factors and regulatory networks, and showed concordant protein expression in the human FT epithelium. Notably, subcluster-level analysis revealed a rare hybrid epithelial-mesenchymal population transcriptionally and proteomically aligned with the mesenchymal-like subtype of HGSOC. Spatial immunoprofiling demonstrated that these cells express multiple defining markers and reside above the basement membrane, intercalated between epithelial cells and lack luminal contact, reminiscent of the histologically described peg cell population. Comparative analyses with fetal mesonephric tissues uncovered developmental parallels, and this population was stably maintained in 3D organoid cultures, providing a robust ex-vivo model to study early transformation. Together, our findings identify a rare, multipotent epithelial population with mesenchymal features as a candidate origin for mesenchymal-like HGSOC and establish the FT organoid biobank as a powerful platform for mechanistic and translational studies in ovarian carcinogenesis. Citation Format: Megan L. Ritting, Wenmei Yang, Syed Mohammed Musheer Aalam, Hui Zhao, Liang Feng, Jianning Song, Mihai G. Dumbrava, Wazim M. Ismail, Kenneth Schaufelberger, S. John Weroha, Scott H. Kaufmann, Alexandre Gaspar-Maia, Mark E. Sherman, Jamie N. Bakkum-Gamez, Nagarajan Kannan. A rare multipotent peg-like epithelial cell is a candidate cell-of-origin for high-grade serous ovarian 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 7636.

ABSTRACT The normal adult male breast has not been characterized at single-cell resolution, leaving the cellular basis of male breast cancer (MBC) biology undefined. Here we present an integrated single-cell RNA sequencing atlas of the adult human breast comprising 174,471 cells from 17 donors (3 male, 14 female), including 18,117 male-derived cells. This revealed that the male breast retains all three epithelial populations, basal (BC), luminal progenitor (LP), and luminal committed cells (LC), but with an increase in LC at the expense of BC and LP across all three male donors. Male LC were distinguished from female by elevated ESR1 and PGR mRNA, enrichment of RNA processing and ribosome biogenesis programs, reduced inflammatory cytokine and growth factor signaling, elevated estradiol gene set enrichment scores, and higher inferred activity of developmental patterning transcription factors. This pattern was observed across differential expression, gene ontology, ligand profiling, and regulon-based analyses, and was not restricted to sex chromosome-linked gene expression. This is consistent with the near-universal estrogen receptor (ER) positivity that characterizes MBC clinically. This atlas provides the first cellular and transcriptional reference for the normal male breast and a resource for investigating sex differences in mammary biology, germline susceptibility variant interpretation, and modeling breast malignancies.

PURPOSE To compare mortality risk among long-term survivors of high-grade serous ovarian cancer (HGSOC) with that of the general US population and to examine how mortality and recurrence trends have evolved under changing treatment approaches. MATERIALS AND METHODS This is a retrospective cohort study including patients with histologically confirmed HGSOC treated at a single referral cancer center from 1991 to 2022. Mortality risk, expressed as standardized mortality ratios (SMRs, ratio of observed deaths in the study population to expected deaths in an age-matched US population), was assessed at diagnosis and yearly on patients without events (disease persistence, recurrence, or death from any cause). RESULTS A total of 2,074 consecutive patients were included. The median follow-up was 12.5 years (IQR, 11.7-13.8). Most patients were stage III (1,396, 69.2%), 1,299 (62.8%) underwent primary cytoreductive surgery, and 1,688 (87.4%) received platinum-taxane as first-line therapy. At diagnosis, the mortality rate was 7.40 times higher than that in the general population (95% CI, 7.04 to 7.77). By 7 event-free years, the 95% CI for the SMR included 1 for the first time, and at 10 years, the SMR was 1.05 (95% CI, 0.72 to 1.49), indicating no statistically significant excess mortality compared with the general population. In those event-free at 10 years, the 5-year cumulative incidence of ovarian cancer–related deaths (8.2% [95% CI, 4.2 to 16.0]) was comparable with that of non–ovarian cancer–related deaths (6.3% [95% CI, 2.9 to 13.8]). CONCLUSION In this cohort, mortality risk steadily declined for patients who remained event-free, approaching that of the general population by 7 years after diagnosis. These results emphasize the importance of individualized, long-term follow-up to address both recurrence risk and survivorship needs.

Abstract Uptake of tamoxifen is limited to 5-15% of all eligible women due to concerns regarding medication side effects combined with the lack of ability to predict which women will benefit from the therapies. We used a unique retrospective sample of n=92 high-risk women with paired benign biopsies obtained before (PreTam) and during tamoxifen (OnTam) to discover tissue gene expression biomarkers of biologic response, defined as change (Δ) in lobular Ki67. Ki67 was digitally quantified in normal lobules from PreTam and OnTam benign tissue samples and samples with PreTam Ki67 = 0 were excluded. Data were analyzed from 51 women with both PreTam and OnTam Ki67 (median age = 49 years and median BMI = 25.5 kg/m2 at PreTam biopsy). OnTam biopsies were collected a median of 11.5 months after PreTam biopsy (median 6.5 months after starting tamoxifen; dose 20 mg/day in 86%). Median Ki67 significantly decreased overall (Δ Ki67 = -0.29%, p=0.003), but the effect was restricted to women with baseline Ki67 ≥ 0.5% (n=38, Δ Ki67 = -0.78%, p=0.0002); no significant change occurred with baseline Ki67 < 0.5 (Δ Ki67 = 0.19%, p=0.23). Δ Ki67 was inversely proportional to baseline Ki67. RNA-seq was successful for 171 benign tissue samples (including 53 PreTam-OnTam pairs). After edgeR normalization (>14,600 genes), PreTam vs OnTam comparison across 33 pairs with baseline Ki67 ≥ 0.5% yielded 66 differentially expressed genes (fold-change > 2 and p<0.01). At p<0.01, ranked gene set enrichment analysis revealed 27 gene sets enriched for reduced expression with Tam (with top-ranked gene set related to steroid hormone biosynthesis) and 59 gene sets enriched for increased expression with Tam (with top-ranking gene sets related to B cell and T cell receptor signaling, regulation of cell cycle, and other immune-related pathways). To identify predictors of response, we correlated baseline PreTam gene expression with Δ Ki67 in 37 women with available baseline PreTam RNAseq and Δ Ki67 data. Partial correlation analysis (via ‘ppcor’ package, adjusted for library qubit) revealed 73 genes with significant positive correlation (rho > 0.50) and 14 genes with significant negative correlation (rho < -0.50) with Δ Ki67 (p<0.01). At p<0.01, ranked gene set enrichment analysis yielded 67 gene sets enriched with negative correlation (i.e. increased expression associated with a reduction in Ki67), including pathways related to ribosome biogenesis, DNA recombination, and nucleic acid metabolic processes; there were 62 gene sets enriched for positive correlation (i.e. decreased expression associated with a reduction in Ki67), including pathways such as VEGF and ERBB signaling. Our data support (i) baseline Ki67 ≥ 0.5 % and (ii) an 87-gene baseline signature as complementary biomarkers of tamoxifen response. Prospective validation of this combined assay could target tamoxifen to the women most likely to benefit, potentially overcoming a major barrier to chemoprevention uptake. Citation Format: S. J. Winham, C. Wang, Y. Liu, B. M. McCauley, D. J. Tabrizi, D. L. Gehling, J. L. Fischer, L. R. Seymour, M. H. Solanki, M. E. Sherman, K. R. Lohani, T. L. Hoskin, D. C. Radisky, A. C. Degnim. Gene expression signatures of beneficial response to tamoxifen for breast cancer prevention [abstract]. In: Proceedings of the San Antonio Breast Cancer Symposium 2025; 2025 Dec 9-12; San Antonio, TX. Philadelphia (PA): AACR; Clin Cancer Res 2026;32(4 Suppl):Abstract nr PS3-01-18.

Abstract Benign breast disease (BBD) confers elevated risk for breast cancer, yet current tools lack precision in identifying lesions likely to progress. Immune mechanisms may influence progression, but the spatial and molecular features of immune responses in BBD versus invasive cancer remain poorly characterized. We applied Digital Spatial Proteomics (DSP) to evaluate immune-related protein expression across 1,540 regions of interest (ROIs) from breast tumors (n=58, Carolina Breast Cancer Study [CBCS]), a UNC BBD cohort (n=20), and a Mayo BBD cohort (n=177). Using Consensus Clustering in CBCS and UNC BBD, we identified two stable immune phenotypes: “High Immune,” enriched in tumors and high-risk benign lesions, and “Quiet Immune,” predominant in low-risk BBD. ROIs were classified by immune marker expression, and patients were labeled as “High” if ≥50% of their ROIs were High Immune. In CBCS, High Immune profiles were associated with ER-negative status (OR=2.92 [0.99-8.65] participant-level; OR=3.28 [1.34-8.04] ROI-level). In BBD, High Immune profiles trended toward higher odds of high-risk lesions (OR=1.00 [0.15-6.70] participant-level; OR=1.68 [0.57-4.97] ROI-level), though this pattern was not replicated in the Mayo cohort. BBD samples exhibited substantial spatial heterogeneity in immune expression (40% within-patient ROI agreement vs. 74% in tumors), potentially limiting biomarker reliability. Nonetheless, four immune proteins (CD44, CD80, CD14, CD66b) were significantly associated with future cancer development in the Mayo cohort, though these findings require validation due to within-person variability. Our findings highlight the challenges posed by spatial heterogeneity in assessing immune biomarkers in BBD and suggest that immune activation may emerge later in disease progression. These results raise important questions about the timing and consistency of immune responses in the natural history of breast cancer. Citation Format: A. Zhang, A. Cozzo, L. Olsson, E. Kirk, X. Gao, S. Nyante, B. Calhoun, R. Vierkant, M. Sherman, M. Troester. Immune Heterogeneity in Breast Cancer Natural History: Benign Breast Disease to Breast Cancer [abstract]. In: Proceedings of the San Antonio Breast Cancer Symposium 2025; 2025 Dec 9-12; San Antonio, TX. Philadelphia (PA): AACR; Clin Cancer Res 2026;32(4 Suppl):Abstract nr PS2-05-23.

Abstract Background: Artificial intelligence (AI) algorithms based on deep learning approaches show promise in improving breast cancer (BC) detection on mammography and may also improve prediction of future BC risk compared with clinical risk prediction models. Historical clinical risk prediction models underperform among women with BBD; however, AI models for BC detection on mammography have not been tested among women at elevated BC risk due to benign breast disease (BBD). We evaluated if an AI cancer detection algorithm can aid BC risk prediction among women with BBD. Methods: We examined women without a prior BC who had an initial BBD biopsy at Mayo Clinic, Rochester, between 2002 and 2013. Incident BC occurring after BBD was identified using the Mayo Tumor Registry and supplemented by follow-up questionnaires. Only women with a full field digital mammogram at least 2 years prior to BC diagnosis were eligible. A breast pathologist evaluated BBD according to increasing severity, as non-proliferative (NP), proliferative disease without atypia (PDWA) or atypical hyperplasia (AH). AI malignancy scores derived from the Transpara detection algorithm (1-10) and volumetric percent density (VPD) from Volpara (per one standard deviation, SD) were assessed from mammograms close to BBD diagnosis. Cox proportional hazards regression models were applied to estimate hazard ratios (HRs) with 95% Confidence Intervals (CIs), adjusted for age and BMI. C-statistics were estimated to assess the contribution of each factor to BC risk prediction. Likelihood ratio tests (LRTs) and bootstrapping methods were used to compare model performance with the addition of AI malignancy score. Results: The BBD cohort with mammography consisted of 3,125 women followed for a median 12.9 years, with 250 incident BC. Of these, 221 BC occurred at least two years after the mammogram and were used in analyses. As expected, increased BC risk was associated with BBD severity [for AH HR=3.18 (95%CI: 2.04, 4.95) and for PDWA HR=1.59 (95%CI: 1.17, 2.14) compared to NP] and with higher VPD [HR=1.24 per SD (95% CI: 1.05, 1.45)] (Table ). The AI-malignancy score alone was also associated with BC risk [HR=1.16 per 1 unit score (95%CI: 1.09,1.23)] and showed similar discriminatory accuracy [C-statistic=0.626 (95%CI: 0.586, 0.665)] as the model with BBD severity and VPD combined [C-statistic=0.627, (95%CI: 0.585, 0.669)] (Table ). In models with BBD severity and VPD, the AI-malignancy score was an independent risk factor for BC [HR=1.14 (95%CI: 1.07, 1.21); PLRT<0.001] but only achieved a marginally significant improvement in discriminatory accuracy [C-statistic=0.651 (0.609, 0.693), ΔC-statistic=0.024 (95% CI: 0.000, 0.047)]. Conclusion: In this preliminary study, an AI BC detection algorithm achieved similar performance in predicting future BC risk compared to established risk factors. Citation Format: C. M. Vachon, S. Winham, M. Jensen, D. Hursh, L. Pacheco-Spann, A. Norman, J. Fischer, S. Schrup, L. Seymour, D. Gehling, S. Nyante, M. Troester, N. Karssemeijer, R. Vierkant, D. Radisky, C. Scott, S. I. Maimone, A. Degnim, M. Sherman. An AI algorithm for breast cancer detection improves future BC risk prediction among women with benign breast disease [abstract]. In: Proceedings of the San Antonio Breast Cancer Symposium 2025; 2025 Dec 9-12; San Antonio, TX. Philadelphia (PA): AACR; Clin Cancer Res 2026;32(4 Suppl):Abstract nr PD4-02.

Abstract Background The incidence of early onset breast cancers (BCs) has increased, paralleling rising trends in delayed childbearing. We hypothesize that a distinct postpartum BC (PPBC) subtype, identifiable by time since last birth (TSLB) and biomarker expression, contributes to this trend. Methods We applied GeoMx Digital Spatial Profiling (DSP) to measure associations between TSLB and 71 proteins in 640 BCs from women aged ≤40 years included in the Young Women’s Breast Cancer Study. We analyzed data using univariable linear regression and multivariable sliced inverse regression to account for higher order interactions among biomarkers. Results In keratin-rich segments, progesterone receptor (PR) (P = 1.00 × 10−4) and PTEN (P = 1.00 × 10−3) were associated with longer TSLB; multivariable analyses revealed positive associations for GZMB (P = 4.00 × 10−4), SMA (P = 1.00 × 10−4), and NF-1 (P = 1.00 × 10−3). In keratin-poor segments, univariable significant positive associations were found for PR (P = 2.00 × 10−4) and PTEN (P = 1.00 × 10−3), whereas CD20 (P = 3.00 × 10−4) and CTLA4 (P = 4.00 × 10−5) were negatively associated; multivariable significant associations were found for fibronectin (P = 3.00 × 10−5) and pan-Akt (P = 1.00 × 10−3). Associations persisted after adjustment for multiple comparisons and BC molecular subtypes. Associations including for PR, PTEN, and CD20 were strongest among women with the shortest TSLB. OPAL multiplex immunofluorescence assays for PR, PTEN, CD20, SMA, and CTLA4 replicated several DSP findings, particularly when stratified by subtype and with compartment matching. In TCGA, RNA species linked to proteins associated with TSLB correlated strongly with a T-cell exhaustion signature previously linked to poor prognosis among premenopausal women. Conclusion These data support PPBC as a biologically coherent phenotype defined by TSLB and biomarker profile, with potential implications for prevention and therapy.

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 High-grade serous ovarian cancer (HGSOC), the most prevalent and lethal ovarian cancer subtype, is often diagnosed at an advanced stage due to the lack of effective early detection methods. Current evidence suggests that HGSOC originates from the fallopian tube (FT), specifically from secretory epithelial cells in the fimbriated region, based on the detection of early serous lesions and clonal genetic relationships linking them to HGSOC tumors. However, efforts to improve early diagnosis and interception are hindered by an incomplete understanding of HGSOC’s cellular and molecular origin. To address this, we established a diverse, clinically and genetically annotated living FT organoid biobank at the Mayo Clinic, encompassing specimens from over 200 patient donors. Using robust, multi-step protocols, we enriched FT epithelial stem/progenitor cells from tissue samples and Tao brushings, generating bi-lineage differentiated organoids independent of anatomical site (ampulla, fimbria) or laterality (right, left). To gain deeper insight into FT epithelial lineage dynamics, we applied multi-omic modalities—including single-cell RNA sequencing (scRNA-seq), single-nucleus RNA and ATAC sequencing (snMultiome), and bulk RNA-seq—to fresh cells and pooled organoids, integrating findings with public FT datasets to generate the largest known integrated single-cell atlas of normal and high-risk FTs. This enabled the generation and validation of gene expression-derived signatures for secretory and multiciliated lineages, identification of key lineage-defining transcription factors and their regulatory networks, and validation using immunohistochemistry (IHC) data from the Human Protein Atlas. Subcluster-level analysis of epithelial cells revealed transcriptionally distinct populations with unique biological functions, deepening insight into epithelial heterogeneity and FT lineage differentiation. Further, analysis of FT epithelial subclusters identified a rare population with transcriptomic and proteomic association with the mesenchymal-like HGSOC subtype. These cells retained epithelial identity while expressing mesenchymal and ECM-related genes, forming a hybrid phenotype. Spatial profiling via IHC localized them to a rare, basal, peg-like position within the FT epithelium. Developmental analysis showed shared transcriptional features with fetal mesonephric cells, suggesting both developmental plasticity and potential vulnerability to transformation. Pooled organoid scRNA-seq confirmed that they are preserved in culture, indicating they are stably maintained ex-vivo. Together, these findings suggest that this cell cluster represents a rare, multipotent epithelial population with mesenchymal features, offering a candidate cellular origin for the mesenchyme-like subtype of HGSOC. In conclusion, this living FT organoid biobank serves as a powerful, clinically annotated platform for dissecting epithelial lineage hierarchies, identifying rare cell populations, and uncovering mechanisms that may underlie subtype-specific ovarian carcinogenesis. Citation Format: Megan Ritting, Syed Mohammed Musheer Aalam, Mihai Dumbrava, Wazim Mohammed Ismail, Alexandre Gaspar-Maia, Mark Sherman, Jamie Bakkum-Gamez, Nagarajan Kannan. Tubal organoid biobanking and multi-omics framework uncovers rare basal peg cells as a source of epithelial stemness and a candidate for high-grade serous ovarian cancer origin [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Advances in Ovarian Cancer Research; 2025 Sep 19-21; Denver, CO. Philadelphia (PA): AACR; Cancer Res 2025;85(18_Suppl):Abstract nr A078.