Spatial mitochondrial lineage tracing uncovers a premetastatic niche and microenvironment programmed fate switching in osteosarcoma

Technical Summary: Spatial Mitochondrial Lineage Tracing Uncovers a Premetastatic Niche and Microenvironment-Programmed Fate Switching in Osteosarcoma

Problem Statement
Osteosarcoma progression is characterized by complex spatiotemporal dynamics, yet the lineage relationships governing tumor cell fate decisions remain poorly understood. While single-cell RNA sequencing (scRNA-seq) has revealed intratumoral heterogeneity, it dissociates cells from their native tissue architecture, erasing spatial coordinates crucial for understanding cellular neighborhoods. Conversely, spatial transcriptomics (ST) preserves architecture but captures only static snapshots, unable to distinguish whether adjacent cells share a common clonal origin or have converged on similar states via convergent evolution. Furthermore, existing lineage tracing methods often rely on transgenic manipulation (e.g., CRISPR barcoding), rendering them inapplicable to primary human clinical specimens. There is a critical need to integrate space, time, and lineage history using endogenous tracers to determine whether pro-metastatic states arise from transient microenvironmental plasticity or deterministic clonal locking.

Methodology
The authors integrated three analytical dimensions—space, time, and lineage—on a single treatment-naïve, late-stage, non-metastatic human osteosarcoma specimen:

1. Single-Cell and Spatial Transcriptomics: scRNA-seq was performed on 11,203 cells to define cellular populations and infer Copy Number Variations (CNVs) to distinguish malignant from non-malignant cells. Spatial transcriptomics (10x Visium) was performed on the same sample to map gene expression onto tissue architecture.
2. Pseudotemporal Reconstruction: The authors applied cell2fate, a Bayesian RNA velocity model, to infer transcriptional dynamics and gene modules, reconstructing the temporal progression of tumor states and stromal remodeling.
3. Endogenous Lineage Tracing: To resolve clonal ancestry without genetic manipulation, the authors utilized naturally occurring somatic mitochondrial DNA (mtDNA) variants. They enriched mitochondrial transcripts from Visium cDNA libraries using the MAESTER protocol.
4. Spatial and Phylogenetic Analysis:
   • Spatial Coherence: A Moran's I-based spatial autocorrelation filter was applied to mitochondrial variant allele frequencies (VAFs) to confirm clonal spatial coherence and filter out depth-driven artifacts.
   • Lineage Reconstruction: A MitoDrift framework, employing a Wright-Fisher drift model coupled with a hidden Markov tree (WF-HMT), was used to reconstruct phylogenetic trees and model stochastic mitochondrial segregation.
   • Fate Transition Modeling: A discrete-state continuous-time Markov model (ARD) was fitted to estimate directional transition rates between cell states (COL3A1⁺, ALPL⁺, THY1⁺, and others).
   • Interaction Analysis: Cell-cell communication was analyzed using CellChat and SpaCET to identify signaling hubs within the spatial context.

Key Results

• Tumor Heterogeneity and CNV Profiling: scRNA-seq identified 12 distinct populations, with tumor cells (66.5%) dominating the landscape. Malignant subpopulations included ALPL⁺ (osteoblastic), THY1⁺ (mesenchymal/stem-like), and COL3A1⁺ (progenitor/matrix-remodeling) cells. CNV profiling confirmed the malignant origin of these subpopulations, with THY1⁺ cells exhibiting a unique gain on chromosome 16.
• Stepwise Fate Bifurcation: Pseudotemporal reconstruction revealed a stepwise cascade where COL3A1⁺ progenitors bifurcate into two distinct lineages: an ALPL⁺ osteoblastic lineage and a THY1⁺ mesenchymal lineage. This bifurcation is temporally coincident with the peak activity of myeloid-driven remodeling modules.
• Divergent Transcriptional Programs:
  • The transition to the ALPL⁺ lineage is driven by an EMT program and a metabolic switch toward oxidative phosphorylation.
  • The transition to the THY1⁺ lineage is characterized by NF-κB-mediated inflammatory reprogramming, antigen processing downregulation, and metabolic rewiring, distinct from the EMT-driven ALPL⁺ path.
• Discovery of the THY1_Endo Niche: Spatial transcriptomics identified a specific pro-metastatic niche where THY1⁺ tumor cells co-localize with PLVAP⁺ dysfunctional endothelial cells. This THY1_Endo niche is enriched for PTN (pleiotrophin) and NOTCH signaling, forming a spatially confined signaling hub. Cell-cell interaction analysis confirmed that THY1⁺ cells engage in the strongest signaling interactions with endothelial cells and specific macrophage subsets, a relationship not observed with other tumor subtypes.
• Lineage-State Coupling and Directional Fate Locking:
  • Spatial Coherence: Mitochondrial variants showed significant spatial autocorrelation (Z = 10.07, p ≈ 4×10⁻²⁴), confirming that mtDNA VAFs retain sufficient clonal information at Visium resolution.
  • Irreversibility: Phylogenetic reconstruction using the WF-HMT model revealed a striking asymmetry in transition rates. The rate from COL3A1⁺ to THY1⁺ was 724.08, compared to only 4.09 for the reverse transition (a ~177-fold asymmetry). In contrast, the COL3A1⁺ to ALPL⁺ transition rate was lower (0.47).
  • Clonal Independence: Low Blomberg's K values (<0.1) indicated that THY1⁺ and ALPL⁺ states arose multiple times independently from different COL3A1⁺ progenitors, yet once established, the THY1⁺ fate is essentially irreversible.
  • Mutation Signatures: Clade-specific mitochondrial mutations were identified in high-purity THY1⁺ and ALPL⁺ clones, suggesting mutational priming precedes phenotypic locking.

Significance and Claims
The paper proposes a "time-space-lineage" relay model for osteosarcoma progression. The authors claim that the metastatic THY1⁺ mesenchymal fate is not governed solely by cell-autonomous programs but requires strict spatial licensing via direct co-localization with dysfunctional endothelia to form the THY1_Endo niche.

Key claims include:

1. Methodological Validation: The study demonstrates that somatic mitochondrial variants can serve as powerful endogenous lineage tracers in primary human tumors, resolving clonal history and fate decisions without transgenic manipulation.
2. Mechanistic Insight: The research establishes that the THY1⁺ state represents an irreversible, clonally locked terminal differentiation driven by a directional transition from COL3A1⁺ progenitors, rather than a transient plastic state.
3. Niche Dependency: The metastatic competence of THY1⁺ cells is contingent upon a spatially restricted signaling environment (PTN/NOTCH) provided by the THY1_Endo niche.
4. Clinical Implication: The THY1_Endo niche is proposed as a tangible spatial target. The authors suggest that therapeutic strategies should aim to disrupt both the lineage commitment (the irreversible fate) and the niche support (the spatial signaling) to intercept the metastatic cascade.

The study concludes that integrating RNA velocity, spatial transcriptomics, and mitochondrial lineage tracing provides a multi-dimensional framework essential for understanding and potentially intercepting cancer evolution, bridging the gap between single-cell states, spatial organization, temporal dynamics, and clonal ancestry.