An APOD-EFEMP1-COL5A2 fibroblast programme marks a matrix-remodelling, vascular and myeloid-associated niche linked to adverse gastric cancer outcomes

Planteamiento del Problema

El pronóstico del cáncer gástrico está fuertemente influenciado por el microambiente tumoral (TME), particularmente por la organización de la matriz extracelular (ECM), la arquitectura vascular y la localización de células inmunes por parte de los fibroblastos asociados al cáncer (CAFs). Sin embargo, la comprensión clínica actual es limitada porque las mediciones de abundancia estromal en bulk colapsan estados de fibroblastos biológicamente distintos en una sola señal. Existe una falta de marcadores compactos y reproducibles que conecten estados transcripcionales específicos de los fibroblastos con el contexto tisular y los resultados clínicos sin sobreajustarse a conjuntos de datos específicos o artefactos técnicos. El estudio tiene como objetivo definir un programa de fibroblastos específico y compacto que pueda rastrearse a través de plataformas de transcriptómica de célula única (single-cell), espacial y bulk para identificar un nicho vinculado a resultados adversos.

Metodología

El estudio emplea una jerarquía de evidencia multicapa y multiplataforma para evaluar un programa de tres genes: APOD-EFEMP1-COL5A2. El análisis distingue estrictamente entre unidades biológicas (células, spots, muestras, cohortes) para evitar la pseudoreplicación.

1. Atlas de Célula Única de Descubrimiento (GSE183904):

   • Integración de 40 muestras de cáncer gástrico (156,111 células controladas por calidad).
   • Se identificaron ocho compartimentos principales. Los fibroblastos fueron purificados utilizando puntuación UCell para genes de identidad de CAF (p. ej., DCN, LUM, FAP) y la exclusión de firmas de contaminación (epitelial, plasmática, inmune, endotelial).
   • Se definieron 2,415 CAFs de alta confianza a partir de 36 muestras.
   • Definición del Programa: Dentro de este conjunto de CAFs de alta confianza, las células se estratificaron en estados "Alto de Programa", "Medio" y "Bajo" basados en los cuartiles por muestra de la media del z-score de APOD, EFEMP1 y COL5A2.

2. Caracterización Transcriptómica:

   • Análisis de Pseudobulk Emparejado: Se agregaron los conteos brutos de 21 muestras que contenían estados de CAF alto y bajo emparejados. Se utilizó edgeR con un diseño emparejado para identificar genes diferencialmente expresados (DEGs) y realizar el enriquecimiento funcional (GO, KEGG, Hallmark).
   • Validación Independiente de Célula Única (GSE163558): Se recalculó la puntuación del programa dentro de un objeto de CAF independiente (255 células, 6 muestras) para evaluar la detectabilidad y la direccionalidad, observando las limitaciones de concentración de muestras.

3. Transcriptómica Espacial (GSE251950 y GSE246011):

   • Se analizaron 10 secciones de tejido (31,202 spots) para probar la colocalización espacial.
   • Se calculó un índice de expresión espacial de tres genes (media de los genes transformados en z) y se correlacionó con las puntuaciones de características de núcleo de CAF, endotelio, angiogénesis y mieloides dentro de cada sección.
   • Se realizó un análisis de enriquecimiento de vecindad mediante pruebas de permutación para evaluar la asociación local con características específicas del TME.

4. Análisis de Transcriptómica Bulk y Clínica (10 Cohortes):

   • Se puntuaron 10 cohortes independientes de bulk (incluyendo TCGA-STAD, GSE15459, GSE26253) utilizando la fórmula del programa.
   • Análisis de Supervivencia: Se estratificaron los pacientes por la mediana de la puntuación del programa. Se realizó una regresión de riesgos proporcionales de Cox para la Supervivencia Global (OS) y la Supervivencia Libre de Recurrencia (RFS). Se realizó un metaanálisis de efectos aleatorios para la OS.
   • Correlación con el TME: Se correlacionaron las puntuaciones del programa con las características estromales, de exclusión inmune y vasculares dentro de cada cohorte antes de resumirlas.

5. Inferencia de Mecanismo y Perturbación:

   • Red Regulatoria: Se clasificaron los factores de transcripción (TFs) basándose en puntuaciones integradas que vinculan la expresión del TF con la actividad del programa.
   • Análisis de Ligando-Receptor: Se utilizó una puntuación de producto de expresión (ligando emisor × receptor receptor) para priorizar las ramas de comunicación, enfocándose específicamente en las vías de matriz-integrina y MIF/CD74.
   • Perturbación Virtual: Se simularon reducciones (knockdowns) de los genes del programa (p. ej., COL5A2) y de los nodos de comunicación (p. ej., MIF/CD74) para identificar vulnerabilidades computacionales. El mapeo de fármacos se limitó a conceptos basados en el conocimiento debido a la falta de datos de IC50 clasificables.

Resultados Clave

1. Identificación y Localización del Programa

• El programa APOD-EFEMP1-COL5A2 se localiza preferencialmente en un compartimento de CAF de alta confianza y purificado.
• En el atlas de descubrimiento, se identificaron 616 CAFs de "Programa-Alto" en 34 muestras. Estas células no estaban confinadas a un solo grupo, sino que mostraron un enriquecimiento preferencial en grupos específicos (p. ej., Clúster 0, OR 3.51).
• Los CAFs de "Programa-Alto" exhibieron un estado transcripcional coherente caracterizado por la regulación positiva de genes de ECM, colágeno, adhesión focal y genes asociados a vasos sanguíneos.

2. Identidad Transcriptómica

• El análisis de pseudobulk emparejado (21 muestras) identificó 261 genes significativos (187 hacia arriba, 74 hacia abajo).
• Las vías reguladas positivamente clave incluyeron la organización de fibrillas de colágeno (NES 2.44) y la morfogénesis de vasos sanguíneos.
• El estado de "Programa-Alto" está fuertemente asociado con el remodelado de la matriz, la señalización de colágeno/integrina y las características vasculares.

3. Validación Independiente

• Célula Única: El programa fue detectable en el conjunto de datos independiente GSE163558, aunque las células de estado alto estaban altamente concentradas en pocas muestras (una muestra contribuyó con el 75% de las células de estado alto). Se observó consistencia direccional para los genes definitorios y las características de la ECM.
• Espacial: En GSE251950, el índice espacial de tres genes mostró correlaciones positivas con las características de núcleo de CAF, endotelio, angiogénesis y mieloides en las 10 secciones. El análisis de vecindad confirmó el enriquecimiento de estas características alrededor de los spots de alto índice.

4. Asociación Clínica

• Supervivencia: Las puntuaciones más altas del programa se asociaron significativamente con peores resultados.
  • TCGA-STAD y GSE15459 (OS): El metaanálisis arrojó un Hazard Ratio (HR) agrupado de 1.80 (9C 1.40–2.33, P = 6.24×10⁻⁶).
  • GSE26253 (RFS): Los tumores con programa-alto tuvieron una supervivencia libre de recurrencia más corta (HR 1.59, 9C 1.18–2.14, P = 0.0025).
• Fenotipo del TME: El programa se correlacionó fuertemente con las firmas de núcleo de CAF, remodelado de ECM/matriz y exclusión inmune estromal en las 10 cohortes de bulk.

5. Modelado Pronóstico

• El programa proporcionó información pronóstica modesta pero reproducible.
• La adición de la puntuación del programa a las variables clínicas (p. ej., edad) mejoró el C-index en 0.060 en TCGA-STAD y en 0.103 en GSE15459.
• El programa añadió un valor incremental más allá de las características genéricas de CAF (p. ej., puntuaciones de núcleo de CAF), con valores de delta C-index positivos a través de las cohortes.

6. Hipótesis Mecanísticas

• Rama Matriz-Integrina-Vascular: Los CAFs de "Programa-Alto" están vinculados computacionalmente con ligandos de matriz/colágeno (p. ej., COL1A1, COL5A2) y receptores de integrina (p. ej., ITGB1) en células endoteliales y de pericitos/músculo liso.
• Rama SP1-MIF/CD74-Mieloide: SP1 fue identificado como un principal candidato regulador aguas arriba. El programa está vinculado a un eje de comunicación MIF-CD74, con puntuaciones de producto de expresión altas para receptores mieloides/macrófagos.
• Perturbación Virtual: El knockdown de COL5A2 mostró la reversión más fuerte predicha del remodelado de la ECM y la señalización de colágeno-integrina. La simulación del bloqueo de MIF/CD74 priorizó los receptores mieloides/macrófagos.

Significancia y Reivindicaciones

El artículo afirma establecer un programa estromal reproducible y biológicamente coherente que vincula la biología de la matriz de los CAF, la organización tisular asociada a la vascularización y el contexto mieloide con resultados adversos en el cáncer gástrico.

• Límites de la Evidencia: El estudio establece explícitamente que no establece un biomarcador universal, un ensayo clínico listo para su uso, ni un mecanismo causal. Los hallazgos se describen como "hipótesis computacionales" y "evidencia de apoyo" en lugar de una prueba definitiva.
• Propuesta de Valor: La contribución primaria es un eje compacto de tres genes que sirve como un control estable para un fenotipo estromal más amplio y complejo. Logra conectar la localización de célula única, el contexto espacial y los resultados clínicos de bulk sin confundir distintas unidades biológicas.
• Potencial Traslacional: El estudio prioriza hipótesis específicas testeables para el seguimiento experimental:
  1. Validar la rama de comunicación matriz-integrina-vascular.
  2. Probar el eje regulador y de comunicación SP1-MIF/CD74.
  3. Investigar COL5A2 y MIF/CD74 como posibles nodos de intervención.
• Limitaciones: Los autores reconocen las limitaciones, incluyendo la naturaleza retrospectiva de los datos públicos, la concentración de muestras en los conjuntos de validación, los límites de resolución de nivel de spot en los datos espaciales y la naturaleza computacional de los resultados de perturbación y mapeo de fármacos. No se reclaman predicciones formales de eficacia de fármacos o respuesta al tratamiento.

En resumen, el artículo proporciona un marco computacional rigurooso y multicapa que define un nicho asociado a CAFs adversos en el cáncer gástrico, ofreciendo un conjunto enfocado de hipótesis para la futura validación experimental respecto al remodelado de la matriz, la organización vascular y las interacciones mieloides.