Abstract Osteosarcoma exhibits profound heterogeneity that has long challenged efforts to understand its mechanisms and advance therapeutic progress. To unravel this complexity, we compiled the largest cross-species single-cell transcriptomic dataset, integrating 775,441 cells from human patients, dog patients, patient-derived xenografts, and mouse models. To our knowledge, this dataset represents the first multi-species, multi-technology, and multi-site (primary and metastatic) harmonization of single-cell data for any solid tumor, enabling a unified framework for interrogating inter- and intra-tumor heterogeneity across biological and evolutionary contexts. Through this work, we define subpopulations of osteosarcoma tumor cells that are conserved across tumors, species, and disease sites. These subpopulations span a continuum of differentiation states, from quiescent progenitor-like cells to more differentiated matrix-producing and inflammatory phenotypes, suggesting a conserved developmental hierarchy. Analysis of the stromal compartment revealed both established and previously unappreciated features of osteosarcoma, including the presence of bone-associated osteoclast-like macrophages in both primary and lung metastatic sites and enrichment of inflammatory and scar-associated macrophages in metastatic lung lesions, which we have previously implicated in metastatic progression. The resulting atlas provides a rich and unprecedented resource for exploration and discovery. Using this resource, we characterized tumor-host interactions occurring in primary and metastatic sites and compared them across species. This analysis revealed a striking number of matrix-derived signals within metastatic lung lesions, far exceeding those identified in primary bone lesions. For example, we found that tumor-derived fibronectin engages syndecans and integrin receptors on epithelial cells, inducing a pathogenic phenotype remarkably similar to that described in pulmonary fibrosis. Validation of these interactions using spatial transcriptomic data identifies distinct neighborhoods that support specific tumor cell subpopulations, with patterns conserved across samples. Collectively, this work establishes a transformative resource and conceptual framework for understanding tumor heterogeneity, evolution, and microenvironmental remodeling. It serves as a powerful platform for hypothesis generation, model fidelity assessment, and therapeutic discovery, guiding the next generation of translational advances in osteosarcoma biology. Citation Format: Yogesh Budhathoki, Matthew Cannon, Troy A. Mceachron, Anand G. Patel, Matthew Gust, Jaime F. Modiano, Dylan T. Ammons, Kathryn Cronise, Daniel Regan, Heather Gardner, Ryan D. Roberts. Single cell RNA-seq analysis of osteosarcoma reveals conserved and distinct ecosystems across sites and species 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 1428.
Building similarity graph...
Analyzing shared references across papers
Loading...
Yogesh Budhathoki
Matthew Cannon
Troy A. Mceachron
Cancer Research
University of Minnesota
National Cancer Institute
Tufts University
Building similarity graph...
Analyzing shared references across papers
Loading...
Budhathoki et al. (Fri,) studied this question.
www.synapsesocial.com/papers/69d1fceba79560c99a0a29dd — DOI: https://doi.org/10.1158/1538-7445.am2026-1428