Longitudinal single-cell and spatial transcriptomics reveals intratumor heterogeneity, therapeutic response, and comparative value of canine marginal zone lymphoma.

UNLABELLED: Diffuse B-cell lymphomas are the most prevalent canine hematologic malignancies, and their clinical presentation resembles that of human B-cell non-Hodgkin lymphomas (NHLs). However, a lack of a clear subtyping framework perpetuates imprecise treatment approaches that fail to address mechanisms of therapy resistance. Here, we employed deep phenotyping via serial sampling and longitudinal transcriptomics to evaluate intratumoral composition and therapy response in a 6-year-old neutered male goldendoodle with stage 5A marginal zone lymphoma (MZL). Following sequential treatment consisting of a single IV dose of oncolytic vesicular stomatitis virotherapy (VSV) and standard CHOP chemotherapy 30 days later, we performed scRNAseq on seven serial lymph node biopsies and spatial sequencing on a resection taken 10 days post-VSV treatment. B cells (>90%) comprised three transcriptionally distinct and temporally stable subpopulations, with specific copy number profiles and robust spatial organization despite disrupted lymph node architecture. Analysis across subpopulations revealed downregulation of the quiescence program including KLF2 , suggesting similarity to KLF2 -deficient human MZL. While VSV successfully reached target B cells, downstream transcriptional responses were predominantly localized to the T-cell compartment. T cells (≤6%) showed increased proliferation and upregulation of cytotoxicity markers post-VSV administration, enriching for leukocyte-mediated cytotoxicity pathways consistent with an anti-viral immune response. Ultimately, these results provide an in vivo proof-of-concept for the treatment paradigm in canine lymphoma, while emphasizing the need for improved subtyping frameworks and multidimensional treatment strategies targeting intratumoral heterogeneity. SIGNIFICANCE: Using longitudinal single-cell and spatial transcriptomics, we mapped the molecular architecture and therapeutic response of canine marginal zone lymphoma in vivo . We identified conserved KLF2 downregulation linking canine and human MZL, alongside stable tumor heterogeneity and a distinct immune response to oncolytic virotherapy. These results demonstrate the utility of deep serial phenotyping to track therapy response, while validating naturally occurring canine cancers as comparative models for human therapeutic development.