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

Resumen Técnico: El rastreo de linaje mitocondrial espacial revela un nicho premetastásico y un cambio de destino programado por el microambiente en el osteosarcoma

Planteamiento del problema
La progresión del osteosarcoma se caracteriza por dinámicas espaciotemporales complejas; sin embargo, las relaciones de linaje que gobiernan las decisiones de destino de las células tumorales aún se comprenden poco. Si bien la secuenciación de ARN de célula única (scRNA-seq) ha revelado la heterogeneidad intratumoral, esta disocia las células de su arquitectura tisular nativa, borrando las coordenadas espaciales cruciales para comprender los vecindarios celulares. Por el contrario, la transcriptómica espacial (ST) preserva la arquitectura pero captura solo instantáneas estáticas, incapaz de distinguir si las células adyacentes comparten un origen clonal común o si han convergido en estados similares mediante evolución convergente. Además, los métodos de rastreo de linaje existentes suelen depender de la manipulación transgénica (por ejemplo, el barcoding por CRISPR), lo que los hace inaplicables a especímenes clínicos humanos primarios. Existe una necesidad crítica de integrar espacio, tiempo e historia de linaje utilizando trazadores endógenos para determinar si los estados prometastásicos surgen de una plasticidad microambiental transitoria o de un bloqueo clonal determinista.

Metodología
Los autores integraron tres dimensiones analíticas —espacio, tiempo y linaje— en un único espécimen de osteosarcoma humano en etapa tardía y no metastásica, libre de tratamiento:

1. Transcriptómica de célula única y espacial: Se realizó scRNA-seq en 11,203 células para definir las poblaciones celulares e inferir variaciones en el número de copias (CNVs) para distinguir células malignas de no malignas. La transcriptómica espacial (10x Visium) se realizó en la misma muestra para mapear la expresión génica sobre la arquitectura tisular.
2. Reconstrucción pseudotemporal: Los autores aplicaron cell2fate, un modelo bayesiano de velocidad de ARN, para inferir la dinámica de transcripción y los módulos génicos, reconstruyendo la progresión temporal de los estados tumorales y la remodelación estromal.
3. Rastreo de linaje endógeno: Para resolver la ascendencia clonal sin manipulación genética, los autores utilizaron variantes naturales de ADN mitocondrial (mtDNA) somático. Enriquecieron los transcritos mitocondriales de las librerías de cDNA de Visium utilizando el protocolo MAESTER.
4. Análisis filogenético y espacial:
   • Coherencia espacial: Se aplicó un filtro de autocorrelación espacial basado en la I de Moran para confirmar la coherencia clonal espacial y filtrar artefactos derivados de la profundidad.
   • Reconstrucción de linaje: Se utilizó un marco de trabajo MitoDrift, empleando un modelo de deriva de Wright-Fisher acoplado con un árbol de Markov oculto (WF-HMT), para reconstruir árboles filogenéticos y modelar la segregación estocástica mitocondrial.
   • Modelado de transición de destino: Se ajustó un modelo de Markov de tiempo continuo y estado discreto (ARD) para estimar las tasas de transición direccionales entre estados celulares (COL3A1⁺, ALPL⁺, THY1⁺ y otros).
   • Análisis de interacción: La comunicación célula-célula se analizó mediante CellChat y SpaCET para identificar centros de señalización dentro del contexto espacial.

Resultados clave

• Heterogeneidad tumoral y perfilado de CNV: La scRNA-seq identificó 12 poblaciones distintas, con las células tumorales (66.5%) dominando el paisaje. Las subpoblaciones malignas incluyeron células ALPL⁺ (osteoblásticas), THY1⁺ (mesenquimales/tipo célula madre) y COL3A1⁺ (progenitoras/remodeladoras de matriz). El perfilado de CNV confirmó el origen maligno de estas subpoblaciones, mostrando las células THY1⁺ una ganancia única en el cromosoma 16.
• Bifurcación de destino escalonada: La reconstrucción pseudotemporal reveló una cascada escalonada donde los progenitores COL3A1⁺ se bifurcan en dos linajes distintos: un linaje osteoblástico ALPL⁺ y un linaje mesenquimal THY1⁺. Esta bifurcación es temporalmente coincidente con el pico de actividad de los módulos de remodelación impulsados por mieloides.
• Programas transcripcionales divergentes:
  • La transición al linaje ALPL⁺ es impulsada por un programa de EMT (transición epitelio-mesénquima) y un cambio metabólico hacia la fosforilación oxidativa.
  • La transición al linaje THY1⁺ se caracteriza por una reprogramación inflamatoria mediada por NF-κB, la regulación a la baja del procesamiento de antígenos y un reajuste metabólico, distinto de la vía de la ALPL⁺ impulsada por la EMT.
• Descubrimiento del nicho THY1_Endo: La transcriptómica espacial identificó un nicho específico prometastásico donde las células tumorales THY1⁺ colocalizan con células endoteliales disfuncionales PLVAP⁺. Este nicho THY1_Endo está enriquecido con señalización de PTN (pleiotrofina) y NOTCH, formando un centro de señalización espacialmente confinado. El análisis de interacción célula-célula confirmó que las células THY1⁺ participan en las interacciones de señalización más fuertes con las células endoteliales y subconjuntos específicos de macrófagos, una relación no observada con otros subtipos tumorales.
• Acoplamiento linaje-estado y bloqueo de destino direccional:
  • Coherencia espacial: Las variantes mitocondriales mostraron una autocorrelación espacial significativa (Z = 10.07, p ≈ 4×10⁻²⁴), confirmando que las frecuencias alélicas de variantes de mtDNA (VAFs) retienen información clonal suficiente a la resolución de Visium.
  • Irreversibilidad: La reconstrucción filogenética mediante el modelo WF-HMT reveló una asimetría sorprendente en las tasas de transición. La tasa de COL3A1⁺ a THY1⁺ fue de 724.08, comparada con solo 4.09 para la transición inversa (una asimetría de ~177 veces). En contraste, la tasa de transición de COL3A1⁺ a ALPL⁺ fue menor (0.47).
  • Independencia clonal: Los bajos valores de K de Blomberg (<0.1) indicaron que los estados THY1⁺ y ALPL⁺ surgieron múltiples veces de forma independiente a partir de diferentes progenitores COL3A1⁺, pero una vez establecidos, el destino THY1⁺ es esencialmente irreversible.
  • Firmas mutacionales: Se identificaron firmas de mutaciones mitocondriales específicas de clados en clones de alta pureza de THY1⁺ y ALPL⁺, lo que sugiere que el cebado mutacional precede al bloqueo fenotípico.

Significación y afirmaciones
El artículo propone un modelo de relevo "tiempo-espacio-linaje" para la progresión del osteosarcoma. Los autores afirman que el destino mesenquimal THY1⁺, metastásicamente competente, no está gobernado únicamente por programas celulares autónomos, sino que requiere una estricta licencia espacial mediante la colocalización directa con endotelios disfuncionales para formar el nicho THY1_Endo.

Las afirmaciones clave incluyen:

1. Validación metodológica: El estudio demuestra que las variantes de ADN mitocondrial somático pueden servir como poderosos trazadores de linaje endógeno en tumores humanos primarios, resolviendo la historia clonal y las decisiones de destino sin manipulación transgénica.
2. Perspectiva mecánica: La investigación establece que el estado THY1⁺ representa un diferenciación terminal irreversible y clonalmente bloqueada, impulsada por una transición direccional desde progenitores COL3A1⁺, en lugar de un estado de plasticidad transitoria.
3. Dependencia del nicho: La competencia metastásica de las células THY1⁺ es contingente a un entorno de señalización espacialmente restringido (PTN/NOTCH) proporcionado por el nicho THY1_Endo.
4. Implicación clínica: El nicho THY1_Endo se propone como un objetivo espacial tangible. Los autores sugieren que las estrategias terapéuticas deberían apuntar a interrumpir tanto el compromiso de linaje (el destino irreversible) como el soporte del nicho (la señalización espacial) para interceptar la cascada metastásica.

El estudio concluye que la integración de la velocidad de ARN, la transcriptómica espacial y el rastreo de linaje mitocondrial proporciona un marco multidimensional esencial para comprender y potencialmente interceptar la evolución del cáncer, cerrando la brecha entre los estados de célula única, la organización espacial, la dinámica temporal y la ascendencia clonal.