Back to Authors

Fernández-Zapico, Martín E.

https://orcid.org/0000-0002-8089-3907

Mayo Clinic

Mayo Clinic

Citations 199

Longitudinal Analysis of Matched Patient Biospecimens Reveals Neural Reprogramming of Cancer-Associated Fibroblasts Following Chemotherapy in Pancreatic Cancer
Henstridge, Aylin Z. , Azorsa, Peter, Arya, Nandini , Elhossiny, Ahmed M. , Kadiyala, Padma , Goldwasser, Molly , Branch, Georgina , Sahai, Vaibhav , Peterson, Nicole , Machicado, Jorge D. , Kwon, Richard S. , Schulman, Allison R. , Wamsteker, Erik‐Jan , Philips, George, Menees, Stacy B. , Xia, Jonathan Y. , Hogenson, Tara L. , Rao, Arvind , Shi, Jiaqi , Frankel, Timothy L. , Bednar, Filip , Magliano, Marina Pasca di , Lyssiotis, Costas A. , Truty, Mark J. , Fernández-Zapico, Martín E. , Carpenter, Eileen S. , Azorsa, Peter A , Philips, George M , Pasca di Magliano, Marina , Lyssiotis, Costas A , Sahai, Vaibhav

Cancer Research | September 21, 2026 | Volume: No volume available | Issue: No issue available | Pages: No page available

10.1158/0008-5472.can-26-0823 Get It

Abstract Pancreatic ductal adenocarcinoma (PDAC), the most common subtype of pancreatic cancer, is a deadly disease with a complex tumor microenvironment (TME). Cytotoxic combination chemotherapy treatments are the mainstay of PDAC therapy, but patients rapidly develop chemoresistance, highlighting the need to elucidate how chemotherapy alters the TME and whether these changes drive chemoresistance. Here, we examined matched pre- and post-treatment tissue specimens using single-cell RNA-sequencing and found near-universal enrichment after treatment of axonal guidance genes in cancer-associated fibroblasts (CAFs). These neural CAFs were enriched near sites of perineural invasion, coinciding with regions of increased tumor cell proliferation, and were enriched in tumor areas distant from nerves after chemotherapy. Metastatic recurrent lesions had the highest prevalence of neural CAFs versus primary tumors and untreated metastases. Neural CAFs showed elevated non-canonical WNT mediators and axonal guidance genes, which complemented matching cognate binding partners in tumor epithelial cells, suggesting a role in tumor-stroma crosstalk. Additionally, in vitro studies revealed enrichment of axonal guidance genes in patient-derived CAF lines upon exposure to chemotherapy or tumor cells. Altogether, these findings implicate fibroblast-derived axonal-guidance genes in promoting PDAC invasion and point to a promising therapeutic target for this disease.
Spatial Configuration of Pancreatic Cancer Is Associated with Disease Recurrence after Neoadjuvant Therapy and Curative-Intent Resection
Cisneros, Luis , Toruner, Merih Deniz , Siddiqui, Zafar , Maraone, Andrea V. , Lin, Miranda , Xiao, Alexander , Thiels, Cornelius A. , Truty, Mark J. , George, Ben , Jazieh, Khalid , Scheel, Robert , Hartley, Christopher , Elhalaby, Rofyda , Dizona, Paul , Fernández-Zapico, Martín E. , Carr, Ryan M. , Cisneros, Luis H , Toruner, Merih D , Xiao, Alexander H , Thiels, Cornelius A , Scheel, Robb , McWilliams, Robert R , Maley, Carlo C. , Hartley, Christopher P , Elhalaby, Rofyda Essam , Dizona, Paul , Shi, Qian , Fernandez-Zapico, Martin E , Carr, Ryan M , Truty, Mark J , Garzaro, Andrea V. Maraone , Lin, Miranda , George, Ben , McWilliams, Robert R. , Shi, Qian

Clinical Cancer Research | July 15, 2026 | Volume: No volume available | Issue: No issue available | Pages: OF1-OF16

10.1158/1078-0432.ccr-25-4968 Get It

Abstract Purpose: Pancreatic ductal adenocarcinoma (PDAC) frequently recurs after neoadjuvant therapy (NT) and curative-intent resection. Although major pathologic response predicts favorable outcomes, most patients achieve only minor response with a heterogeneous recurrence risk. We asked whether the spatial organization of residual PDAC encodes clinically relevant biology beyond residual tumor burden. Experimental Design: In a retrospective cohort of 203 patients with resected PDAC treated with NT and restricted to minor pathologic response, routine hematoxylin and eosin (H&E) whole-slide images were segmented into cancer and stroma using an artificial intelligence–enabled pipeline. Spatial composition (patch density, edge density) and configuration (compactness/complexity, intermixing) were quantified and tested for associations with disease-free survival (DFS) using multivariable models adjusted for standard clinicopathologic factors. Results: Shorter DFS was associated with a fragmented, interface-rich tumor–stroma ecology featuring higher edge density and diversity and reduced homotypic aggregation, independent of clinicopathologic variables. Two spatial risk models were independently prognostic: (i) cancer mean shape index plus stromal shape index variability [adjusted hazard ratio (HR), 1.71; P = 0.003] and (ii) mean stromal patch area plus edge density (adjusted HR, 2.19; P = 0.002). Both models stratified outcomes in which pathologic response grading and residual cancer area did not. High-risk configurations were further associated with reduced intratumoral tumor-infiltrating lymphocyte (TIL) density and infiltration ratio, with relative TIL accumulation at the cancer periphery and within the stroma, consistent with an immune-excluded phenotype. Conclusions: Residual cancer–stroma topology quantified from standard H&E slides yields independent prognostic signals after NT in PDAC, provides a cellular immune correlate for spatial risk, and motivates prospective validation and spatially informed adjuvant strategies.
Histone lysine methyltransferases KMT2C and KMT2D join the all-star tumor suppressor team in gastrointestinal cancer
Ruiz, Nicole M. Peña, Fernandez-Zapico, Martin E., Peña Ruiz, Nicole M, Ruiz, Nicole M. Pena , Fernández-Zapico, Martín E.

Journal of Clinical Investigation | July 14, 2026 | Volume: 136 | Issue: 14 | Pages: No page available

10.1172/jci208392 Get It

Members of the type 2 histone lysine methyltransferase family (KMT2s) are key drivers of enhancer activation and are the most mutated group of epigenetic regulators in different cancer types. Within this family, KMT2C and KMT2D have the highest mutational incidence across various cancers. To evaluate their role in gastric cancer, Wang et al. developed a Pten deficiency–driven genetically engineered mouse model with inducible loss of Kmt2c and Kmt2d in gastric epithelial cells. Through extensive in vitro, in vivo, and in silico analyses, the authors revealed that the concomitant loss of Kmt2c and Kmt2d promotes gastric carcinogenesis while enhancing antigen presentation and sensitivity to immunotherapy and targeted approaches like mTOR inhibition, highlighting the tumor-suppressive roles of KMT2C/D in gastric cancer and uncovering a vulnerability for this dismal condition.
A Novel Emerin/Nesprin-3 Interaction Mediates Nuclear Size Changes Induced by Oncogenic KRAS in Pancreatic Cancer Cells
LaRue-Nolan, Kayla C. , Flores, Luis F., Vera, Renzo E., Sigafoos, Ashley N. , Garcia, Laura Hilario, Almada, Luciana L. , Tader, Brooke R , Rodriguez-Quevedo, Maria F. , Marks, David L. , Roux, Kyle J., Fernandez-Zapico, Martin E., Flores, Luis F. , Hilario Garcia, Laura, Vera, Renzo E., Hilario-Garcia, Laura, May, Danielle G. , Roux, Kyle J. , Fernández-Zapico, Martín E.

Molecular and Cellular Biology | May 23, 2026 | Volume: 46 | Issue: 7 | Pages: 746-765

10.1080/10985549.2026.2670710 Get It

Pancreatic tumor microenvironment reprogramming via alloantigen-expressing virotherapy elicits tumor rejection and improves immunotherapy response
Tesfay, Mulu Z. , Cios, Aleksandra , Cheng, Zetao , Ferdous, Khandoker Usran, Shelton, Randal S., Mustafa, Bahaa , Elliott, Natalie , Simões, Camila , Miousse, Isabelle R. , Siegel, Eric R. , Govindarajan, Rang, Cannon, Martin J. , Fernández-Zapico, Martín E. , Elliott, Natalie , Raupach, Elizabeth A. , Chabu, Chiswili Yves , Cios, Aleksandra , Urbaniak, Alicja Joanna , Bauer, Michael , Post, Steven R. , Chamcheu, Jean Christopher , Grdzelishvili, Valery Z. , Basnakian, Alexei G. , Moaven, Omeed , Borad, Mitesh J. , Nagalo, Bolni Marius

Journal for ImmunoTherapy of Cancer | May 12, 2026 | Volume: 14 | Issue: 5 | Pages: e013671

10.1136/jitc-2025-013671 Get It

Background Pancreatic ductal adenocarcinoma (PDAC), the most common malignant type of pancreatic cancer, is characterized by a dense desmoplastic stroma, low neoantigen burden, and a highly immunosuppressive tumor microenvironment (TME). Methods We present a novel strategy harnessing acute transplant rejection mechanisms by employing a recombinant oncolytic rVMG vector engineered to express murine alloantigens H-2Kb (rVMG-H-2Kb) or H-2Kk (rVMG-H-2Kk), thereby inducing tumor-specific antigenic mismatch responses. Results In vitro, rVMG-H-2Kb and rVMG-H-2Kk exhibited strong replication and cytolytic activity while inducing cell surface expression of both H-2Kk and endogenous H-2Kb, in addition to upregulation of antigen presentation genes (β2-microglobulin, Tap1, and Tapbp). In two immunocompetent PDAC models, intratumoral and systemic delivery of rVMG-H-2Kb and rVMG-H-2Kk delayed tumor progression, with rVMG-H-2Kk conferring a survival advantage. Multiplex immunohistochemistry and immunophenotyping revealed substantial TME remodeling, marked by increased effector T-cell infiltration, regulatory T-cell depletion, and reduced fibrosis. Spatial transcriptomics further showed compartment-specific immune activation and epithelial metabolic reprogramming, corroborating with enhanced tumor immunogenicity. Despite these effects, rVMG-H-2Kk also induced compensatory immunosuppressive pathways, including upregulation of antiviral response genes and immune checkpoint receptors such as programmed death-ligand-2. Importantly, combination therapy with rVMG-H-2Kk and murine checkpoint blockade (anti-programmed cell death protein-1 and anti-cytotoxic T-lymphocyte-associated protein 4) drastically improved survival compared to checkpoint blockade alone. Strikingly, surviving mice resisted tumor rechallenge, indicating the establishment of durable antitumor memory. Conclusion These findings establish rVMG-mediated alloantigen delivery as a novel immunotherapeutic platform capable of converting immune-cold tumors into immune-hot, sensitizing tumors to immune checkpoint inhibitors, and establishing durable antitumor immunity in PDAC.
Longitudinal Analysis of Matched Patient Biospecimens Reveals Neural Reprogramming of Cancer-Associated Fibroblasts Following Chemotherapy in Pancreatic Cancer
Henstridge, Aylin Z. , Arya, Nandini , Kadiyala, Padma , Branch, Georgina , Peterson, Nicole , Kwon, Richard S. , Schulman, Allison R. , Wamsteker, Erik, Philips, George, Menees, Stacy B. , Hogenson, Tara L., Rao, Arvind, Shi, Jiaqi , Bednar, Filip , Lyssiotis, Costas A. , Truty, Mark J. , Fernández-Zapico, Martín E. , Shi, Jiaqi , Bednar, Filip , Pasca di Magliano, Marina , Lyssiotis, Costas A , Carpenter, Eileen S , Elhossiny, Ahmed M. , Sahai, Vaibhav , Machicado, Jorge D. , Xia, Jonathan Y. , Hogenson, Tara L., Rao, Arvind, Frankel, Timothy L. , Magliano, Marina Pasca di , Carpenter, Eileen S.

No container title available | December 8, 2025 | Volume: No volume available | Issue: No issue available | Pages: No page available

10.64898/2025.12.01.691614 Get It

Abstract Pancreatic ductal adenocarcinoma (PDAC), the most common subtype of pancreatic cancer, is a deadly disease with a complex tumor microenvironment (TME). How chemotherapy alters the TME, and whether these changes drive chemoresistance, is poorly understood. We examined matched pre- and post-treatment tissue specimens and found near-universal enrichment of axonal guidance genes in cancer-associated fibroblasts (CAFs) after treatment. These CAFs were enriched near sites of perineural invasion, coinciding with regions of increased tumor cell proliferation, and were enriched in tumor areas distant from nerves after chemotherapy. Metastatic recurrence lesions had the highest prevalence of these CAFs versus primary tumors and untreated metastasis. These CAFs showed elevated non-canonical Wnt mediators and axonal guidance genes, which complemented matching cognate binding partners in tumor epithelial cells, suggesting a role in tumor-stroma crosstalk. Our findings implicate fibroblast-derived axonal-guidance genes in promoting PDAC invasion and point to a promising target for this disease. Statement of Significance Therapeutic resistance remains a major challenge in pancreatic cancer. Our findings shed light on the changes in the complex tumor microenvironment in response to chemotherapy and identify fibroblasts high in axonal-guidance genes that may drive tumor progression and chemoresistance, thus uncovering a potential avenue for targeting pancreatic cancer.
MC210807: Phase 1 clinical trial assessing the safety and efficacy of onvansertib, a novel, oral, PLK1 inhibitor in Relapsed/Refractory myeloproliferative chronic myelomonocytic leukemia (CMML).
Mangaonkar, Abhishek A. , Geyer, Susan M. , Burgess, Ashya, Rice, Kimberly, Asmus, E.J., Koenig, Patricia, Lasho, Terra, Finke, Christy M. , Reither, Rachel, Al‐Kali, Aref , Fernández-Zapico, Martín E. , Patnaik, Mrinal M. , McCullough, Kristen , Rice, Kimberly, Lasho, Terra, Carr, Ryan M , Fernandez, Jenna , Saliba, Antoine N.

Blood | December 5, 2025 | Volume: 146 | Issue: Supplement 1 | Pages: 5634-5634

10.1182/blood-2025-5634 Get It

Abstract BACKGROUND: DNA methyltransferase inhibition offers palliative benefit but does not alter disease biology. in chronic myelomonocytic leukemia (CMML, Merlevede J et al. Nat Comm 2016). This was conclusively demonstrated in the phase 3 DACOTA study, where decitabine did not improve event-free survival in comparison to hydroxyurea (median 12.1 versus 10.3 months, P=0.27) in advanced myeloproliferative CMML (MP-CMML). Hence, there is in urgent need to develop disease modifying therapies for this dismal condition. We and others have demonstrated a role for RAS mutations in MP-CMML through the activation of the KMT2A-PLK1 axis. Targeting of this axis using onvansertib, an inhibitor of PLK1, has demonstrated its pre-clinical efficacy in MP-CMML patient-derived organoid and xenograft models. (Carr RM et al. Nat Comm 2021). METHODS: We conducted a single center phase 1 study to assess the safety and efficacy of onvansertib in relapsed/refractory MP-CMML (NCT05549661). We used a Bayesian Optimal Interval (BOIN) study design to determine the maximum tolerated dose (MTD). The three dose cohorts included 6, 9 and 12 mg/m2 given orally from days 1 to 21 of each 28-day cycle. Key inclusion criteria included a WHO-defined diagnosis of MP-CMML with a white blood cell (WBC) count > 13,000/mm3 and adequate organ function (estimated GFR > 50 mL/min/m2). Response was evaluated per the MDS/MPN International Working Group Criteria (Savona et al. Blood 2015). Targeting efficacy was assessed by measuring phosphorylation of TCTP (a PLK1 substrate), specifically pTCTP/TCTP ratios by immunoblotting in peripheral blood mononuclear cell (PBMC) samples pre and post cycle 1 (Days 1, 7, 21) and end of cycle 3. Pharmacokinetic studies and translational correlative studies were conducted per protocol. RESULTS: Between April 2023 and July 2025, a total of 9 patients were enrolled. The median age at enrollment was 68 years (range: 60-83); all caucasian with 4 (44%) males. All patients had a diagnosis of CMML-1 with a median bone marrow (BM) blast percentage of 4 (0-9), median hemoglobin 10 (7.1-13.2) g/dL, median WBC count of 11.1 (3.1-72.7) x 109/L (4 patients had a WBC count < 13 x 109/L at enrollment due to hydroxyurea use but were confirmed to be MP-CMML at diagnosis), absolute monocyte count 2.1 (1.9-19.6) x 109/L, and platelet count of 36 (20-494) x 109/L. TET2 (89%), SRSF2 (78%), ASXL1 (56%), NRAS, KRAS and CBL (22% each) were the most frequent molecular abnormalities. Previous therapies included hydroxyurea (n=8) and decitabine (n=1). Four patients were enrolled at the 6mg/m2dose level, 2 at the 9 mg/m2 dose level and we currently have one patient at the 12 mg/m2 dose level. While there were 2 planned patients at the 6 mg/m2 dose level, due to a case of pericardial effusion that was possibly drug related, an additional two patients were enrolled at the same dose level, without any dose limiting toxicity (DLT). Targeting efficacy was confirmed by pTCTP immunoblotting in PBMC samples from 4 patients (6 and 9 mg/m2 dose levels) showing suppressed pTCTP/TCTP ratios. No confirmed DLT has been identified thus far, and no grade 3 or 4 toxicities attributable to the study drug have been reported. Possible or unrelated grade 3 adverse effects included infections (sepsis/cellulitis, n=2, 22%), acute kidney injury (n=1, 11%), pericardial effusion (n=1, 11%), and atrial fibrillation with rapid ventricular response (n=1, 11%). Median duration on therapy was 124 (range: 15-654) days with most common cause of study discontinuation being disease progression (33%). Response was assessed as per the MDS/MPN IWG criteria in patients who received at least 1 cycle (n=8) with 3 (38%) patients meeting criteria for either hematological or optimal marrow response. Specifically, there was one (14%) patient with >100% platelet response, one (14%) with normalization of WBC count (baseline WBC count >50 x 109/L), one (14%) optimal marrow response with BM blast% reduction from 8 to 2, facilitating candidacy for an allogeneic hematopoietic stem cell transplantation, and three patients (38%) remained as stable disease. CONCLUSION: In MP-CMML, PLK1 inhibition with onvansertib is relatively well tolerated and shows preliminary efficacy in approximately 40% of patients, with one patient having an optimal marrow response at the 9 mg/m2 dose. Future directions include a dose expansion phase and studying upfront combination with DNA methyltransferase inhibitors.
Abstract A015: Tumor landscape analysis: An ecologically informed framework to understand tumor microenvironments
Fernández-Zapico, Martín E. , Carr, Ryan M , Cisneros, Luis , Toruner, Merih Deniz , Maley, Carlo C.

Cancer Research | December 4, 2025 | Volume: 85 | Issue: 23_Supplement | Pages: A015-A015

10.1158/1538-7445.canevol25-a015 Get It

Abstract Background: The spatial organization of the tumor microenvironment (TME) encodes the selective pressures that drive cancer evolution, but many tools focus on local cell contacts or single-modality readouts. We developed Tumor Landscape Analysis (TLA), a modular framework that borrows from landscape ecology and spatial statistics to quantify heterogeneity, fragmentation, and diversity across scales, from single cells to whole-slide ecosystems. Methods: TLA ingests whole-cell or nuclear centroids and region-level segmentations to compute (i) point statistics such as nearest-neighbor distance, Ripley’s H, and Morisita-Horn colocalization; (ii) kernel-smoothed local spatial profiles including Getis-Ord Gi* hot and cold spots; and (iii) fragmentation, shape, and diversity metrics such as patch density, contagion, shape index, and Shannon entropy. We also introduce Local Microenvironments (LMEs), unsupervised ecological niches defined by joint cell-type abundance and spatial mixing. Pilot analyses were performed in ductal carcinoma in situ (DCIS; n=353 samples, 175 patients) with cell-level maps of epithelial, lymphocyte, and fibroblast populations, and in pancreatic ductal adenocarcinoma (PDAC; n=203 resected cases after neoadjuvant therapy) using cancer and stroma segmentations. Results: In DCIS, spatial relationships rather than raw abundance associated with recurrence risk. Shorter lymphocyte-to-fibroblast minimum nearest-neighbor distance and lower Ripley’s H, indicating greater intermixing, correlated with earlier recurrence independent of cell counts. This highlights immune–stromal proximity as a candidate eco-evolutionary driver of persistence. In PDAC, tissue-level fragmentation captured prognostic architecture after therapy. Two shape-based signatures emerged: (1) stromal uniformity, where lower standard deviation of stromal patch shape index associated with improved disease-free survival, consistent with a more constrained, less reactive stroma; and (2) tumor geometric complexity, where higher mean cancer patch shape index, reflecting more irregular or dendritic nests, associated with lower recurrence risk, whereas compact nests predicted worse outcomes. Across datasets, LME maps localized niche phenotypes, such as mixed-immune or segmented tumor habitats, that align with hypothesized selection landscapes and provide inputs for fragmentation analysis. Conclusions: TLA operationalizes an eco-evolutionary view of cancer by unifying interpretable, multiscale spatial metrics in a single workflow that works with any rasterized data input. By quantifying how habitats, edges, and mixing patterns are organized, TLA reveals niche structures and architectural motifs linked to clinical outcomes and complements foundation models and modality-specific pipelines with theory-grounded, clinically translatable descriptors. Ongoing work integrates TLA features with molecular data and outcomes to derive robust biomarkers and to inform adaptive, niche-targeted strategies in precision oncology. Citation Format: Ryan M. Carr, Luis Cisneros, Merih Deniz. Toruner, Martin E. Fernandez-Zapico, Carlo Maley. Tumor landscape analysis: An ecologically informed framework to understand tumor microenvironments [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Cancer Evolution: The Dynamics of Progression and Persistence; 2025 Dec 4-6; Albuquerque, NM. Philadelphia (PA): AACR; Cancer Res 2025;85(23_Suppl):Abstract nr A015.
Novel DAG‐Lactone‐Based Vav1 Inhibitors Show Anti‐proliferative Activity in Pancreatic Cancer Models
Bellomo, Ana , Törüner, Murat , Elhalem, Eleonora , Peach, Megan L. , Donadio, Lucía Gandolfi , Hogenson, Tara L. , Fernández-Zapico, Martín E. , Comin, María J. , Fernandez-Zapico, Martin E , Comin, María J , Toruner, Murat, Giordano, Luciana

ChemMedChem | October 20, 2025 | Volume: 20 | Issue: 23 | Pages: No page available

10.1002/cmdc.202500589 Get It

Targeting C1 domains is a promising strategy for modulating the activity of signaling proteins driving tumor growth and progression. While most small molecules developed to date have focused on typical C1 domains, the development of regulators targeting atypical C1 domains remains under‐explored. Herein, we describe the design and synthesis of novel cationic diacylglycerol DAG‐lactones to efficiently interact with the negatively charged residues present in the atypical C1 domain of Vav1, a guanine nucleotide exchange factor playing a critical role in tumor development, including pancreatic cancer. We evaluated the therapeutic potential of this new family of compounds using models from this dismal condition where Vav1 is aberrantly expressed. Treatment of cultured pancreatic tumor cells with sn‐1 cationic DAG‐lactones inhibited proliferation of Vav1‐expressing cells while Vav1‐negative cells showed no response. Additionally, we demonstrated that these compounds inhibited growth of patient‐derived organoids models of pancreatic cancer. These findings underscore the translational value of these cationic DAG‐lactones for pancreatic cancer patients expressing Vav1 and serve as foundation for future approaches targeting atypical C1 domain‐containing signaling proteins.