Search

Purpose: Lung metastasis is responsible for nearly all deaths caused by osteosarcoma, the most common pediatric bone tumor. How malignant bone cells coerce the lung microenvironment to support metastatic growth is unclear. The purpose of this study is to identify metastasis-specific therapeutic vulnerabilities by delineating the cellular and molecular mechanisms underlying osteosarcoma lung metastatic niche formation., Experimental Design: Using single-cell transcriptomics (scRNA-seq), we characterized genome- and tissue-wide molecular changes induced within lung tissues by disseminated osteosarcoma cells in both immunocompetent murine models of metastasis and patient samples. We confirmed transcriptomic findings at the protein level and determined spatial relationships with multi-parameter immunofluorescence and spatial transcriptomics. Based on these findings, we evaluated the ability of nintedanib, a kinase inhibitor used to treat patients with pulmonary fibrosis, to impair metastasis progression in both immunocompetent murine osteosarcoma and immunodeficient human xenograft models. Single-nucleus and spatial transcriptomics was used to perform molecular pharmacodynamic studies that define the effects of nintedanib on tumor and non-tumor cells within the metastatic microenvironment., Results: Osteosarcoma cells induced acute alveolar epithelial injury upon lung dissemination. scRNA-seq demonstrated that the surrounding lung stroma adopts a chronic, non-resolving wound-healing phenotype similar to that seen in other models of lung injury. Accordingly, metastasis-associated lung demonstrated marked fibrosis, likely due to the accumulation of pathogenic, pro-fibrotic, partially differentiated epithelial intermediates and macrophages. Our data demonstrated that nintedanib prevented metastatic progression in multiple murine and human xenograft models by inhibiting osteosarcoma-induced fibrosis., Conclusions: Fibrosis represents a targetable vulnerability to block the progression of osteosarcoma lung metastasis. Our data support a model wherein interactions between osteosarcoma cells and epithelial cells create a pro-metastatic niche by inducing tumor deposition of extracellular matrix proteins such as fibronectin that is disrupted by the anti-fibrotic TKI nintedanib. Our data shed light on the non-cell autonomous effects of TKIs on metastasis and provide a roadmap for using single-cell and spatial transcriptomics to define the mechanism of action of TKI on metastases in animal models.

Fusion-positive rhabdomyosarcoma is an aggressive pediatric cancer molecularly characterized by arrested myogenesis. The defining genetic driver, PAX3::FOXO1, functions as a chimeric gain-of-function transcription factor. An incomplete understanding of PAX3::FOXO1's in vivo epigenetic mechanisms has hindered therapeutic development. Here, we establish a PAX3::FOXO1 zebrafish injection model and semi-automated ChIP-seq normalization strategy to evaluate how PAX3::FOXO1 initially interfaces with chromatin in a developmental context. We investigated PAX3::FOXO1's recognition of chromatin and subsequent transcriptional consequences. We find that PAX3::FOXO1 interacts with inaccessible chromatin through partial/homeobox motif recognition consistent with pioneering activity. However, PAX3::FOXO1-genome binding through a composite paired-box/homeobox motif alters chromatin accessibility and redistributes H3K27ac to activate neural transcriptional programs. We uncover neural signatures that are highly representative of clinical rhabdomyosarcoma gene expression programs that are enriched following chemotherapy. Overall, we identify partial/homeobox motif recognition as a new mode for PAX3::FOXO1 pioneer function and identify neural signatures as a potentially critical PAX3::FOXO1 tumor initiation event.

Osteosarcoma is an aggressive malignancy characterized by high genomic complexity. Identification of few recurrent mutations in protein coding genes suggests that somatic copy-number aberrations (SCNA) are the genetic drivers of disease. Models around genomic instability conflict-it is unclear whether osteosarcomas result from pervasive ongoing clonal evolution with continuous optimization of the fitness landscape or an early catastrophic event followed by stable maintenance of an abnormal genome. We address this question by investigating SCNAs in >12,000 tumor cells obtained from human osteosarcomas using single-cell DNA sequencing, with a degree of precision and accuracy not possible when inferring single-cell states using bulk sequencing. Using the CHISEL algorithm, we inferred allele- and haplotype-specific SCNAs from this whole-genome single-cell DNA sequencing data. Surprisingly, despite extensive structural complexity, these tumors exhibit a high degree of cell-cell homogeneity with little subclonal diversification. Longitudinal analysis of patient samples obtained at distant therapeutic timepoints (diagnosis, relapse) demonstrated remarkable conservation of SCNA profiles over tumor evolution. Phylogenetic analysis suggests that the majority of SCNAs were acquired early in the oncogenic process, with relatively few structure-altering events arising in response to therapy or during adaptation to growth in metastatic tissues. These data further support the emerging hypothesis that early catastrophic events, rather than sustained genomic instability, give rise to structural complexity, which is then preserved over long periods of tumor developmental time., Significance: Chromosomally complex tumors are often described as genomically unstable. However, determining whether complexity arises from remote time-limited events that give rise to structural alterations or a progressive accumulation of structural events in persistently unstable tumors has implications for diagnosis, biomarker assessment, mechanisms of treatment resistance, and represents a conceptual advance in our understanding of intratumoral heterogeneity and tumor evolution., (© 2023 The Authors; Published by the American Association for Cancer Research.)

DNA methylation is traditionally thought to be established during early development and to remain mostly unchanged thereafter in healthy tissues, although recent studies have shown that this epigenetic mark can be more dynamic. Epigenetic changes occur in the liver after birth, but the timing and underlying biological processes leading to DNA methylation changes are not well understood. We hypothesized that this epigenetic reprogramming was the result of terminal differentiation of hepatocyte precursors. Using genomic approaches, we characterized the DNA methylation patterns in mouse liver from E18.5 until adulthood to determine if the timing of the DNA methylation change overlaps with hepatocyte terminal differentiation, and to examine the genomic context of these changes and identify the regulatory elements involved. Out of 271,325 CpGs analyzed throughout the genome, 214,709 CpGs changed DNA methylation by more than 5% ( e.g. , from 5 to 10% methylation) between E18.5 and 9 wk of age, and 18,863 CpGs changed DNA methylation by more than 30%. Genome-scale data from six time points between E18.5 and P20 show that DNA methylation changes coincided with the terminal differentiation of hepatoblasts into hepatocytes. We also showed that epigenetic reprogramming occurred primarily in intergenic enhancer regions while gene promoters were less affected. Our data suggest that normal postnatal hepatic development and maturation involves extensive epigenetic remodeling of the genome, and that enhancers play a key role in controlling the transition from hepatoblasts to fully differentiated hepatocytes. Our study provides a solid foundation to support future research aimed at further revealing the role of epigenetics in stem cell biology., (Copyright © 2016 Cannon et al.)