scRNAseq Analysis of Osteosarcoma Reveals Conserved and Distinct Ecosystems Across Sites and Species

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 dataset 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. ]]>
Inclusion Criteria: Subjects must meet one of the following criteria: Patients of any age group diagnosed with a malignant, intermediate or low grade solid tumors Patients with newly diagnosed solid tumors who have imaging suggestive of a tumor with malignant potential who are undergoing tumor biopsy or resectionPatients previously diagnosed with a solid tumor that are undergoing a subsequent biopsy or resection of primary, metastatic or recurrent diseasePatients with a diagnosis of a solid tumor undergoing treatment. Patients with a diagnosis of a solid tumor who have completed treatment within the last 5 yearsSubjects must be able to understand and sign the informed consent form directly or through an authorized representativeAll existing enrollments on IRB#11-00478 will be rolled over to this study with continued specimen and data collection.There were no exclusion criteria.]]>