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

Definição do Problema

O prognóstico do câncer gástrico é fortemente influenciado pelo microambiente tumoral (TME), particularmente pela organização da matriz extracelular (ECM), arquitetura vascular e localização de células imunes por fibroblastos associados ao câncer (CAFs). No entanto, o entendimento clínico atual é limitado porque as medições de abundância estromal em bulk colapsam estados de fibroblastos biologicamente distintos em um único sinal. Há uma carência de marcadores compactos e reprodutíveis que conectem estados transcricionais específicos de fibroblastos ao contexto tecidual e aos desfechos clínicos sem sofrer overfitting para conjuntos de dados específicos ou artefatos técnicos. O estudo visa definir um programa de fibroblastos específico e compacto que possa ser rastreado através de plataformas de transcriptômica de célula única (single-cell), espacial e bulk para identificar um nicho ligado a desfechos adversos.

Metodologia

O estudo emprega uma hierarquia de evidências multicamadas e multiplataforma para avaliar um programa de três genes: APOD-EFEMP1-COL5A2. A análise distingue estritamente entre unidades biológicas (células, spots, amostras, coortes) para evitar pseudorreplicação.

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

   • Integrou 40 amostras de câncer gástrico (156.111 células controladas por qualidade).
   • Identificou oito compartimentos principais. Os fibroblastos foram purificados usando a pontuação UCell para genes de identidade de CAF (ex: DCN, LUM, FAP) e exclusão de assinaturas de contaminação (epitelial, plasmócito, imune, endotelial).
   • Definiu 2.415 CAFs de alta confiança a partir de 36 amostras.
   • Definição do Programa: Dentro deste conjunto de CAFs de alta confiança, as células foram estratificadas em estados "Alto", "Médio" e "Baixo" baseados nos quartis por amostra do z-score médio de APOD, EFEMP1 e COL5A2.

2. Caracterização Transcricional:

   • Análise de Pseudobulk Pareada: Agregou contagens brutas de 21 amostras contendo estados de CAF alto e baixo pareados. Utilizou edgeR com um delineamento pareado para identificar genes diferencialmente expressos (DEGs) e realizar enriquecimento funcional (GO, KEGG, Hallmark).
   • Validação Independente de Célula Única (GSE163558): Recalculou a pontuação do programa dentro de um objeto de CAF separado (255 células, 6 amostras) para avaliar a detectabilidade e a direcionalidade, observando as limitações de concentração de amostras.

3. Transcriptômica Espacial (GSE251950 & GSE246011):

   • Analisou 10 seções de tecido (30.202 spots) para testar a colocalização espacial.
   • Calculou um índice de expressão espacial de três genes (média dos genes transformados em z) e o correlacionou com escores de características de núcleo de CAF, endotélio, angiogênese e mieloide dentro de cada seção.
   • Realizou análise de enriquecimento de vizinhança usando testes de permutação para avaliar a associação local com características específicas do TME.

4. Transcriptômica de Bulk e Análise Clínica (10 Coortes):

   • Pontuou 10 coortes independentes de bulk (incluindo TCGA-STAD, GSE15459, GSE26253) usando a fórmula do programa.
   • Análise de Sobrevivência: Estratificou pacientes pelo valor de mediana da pontuação do programa. Realizou regressão de riscos proporcionais de Cox para Sobrevivência Global (OS) e Sobrevivência Livre de Recorrência (RFS). Realizou meta-análise de efeitos aleatórios para OS.
   • Correlação com o TME: Correlacionou as pontuações do programa com características estromais, de exclusão imune e vasculares dentro de cada coorte antes de resumir.

5. Inferência de Mecanismo e Perturbação:

   • Rede Regulatória: Classificou fatores de transcrição (TFs) baseando-se em escores integrados que ligam a expressão do TF à atividade do programa.
   • Análise de Ligante-Receptor: Utilizou um escore de produto de expressão (ligante emissor × receptor receptor) para priorizar ramos de comunicação, focando especificamente nas vias de matriz-integrina e MIF/CD74.
   • Perturbação Virtual: Simulou nocaute (knockdown) de genes do programa (ex: COL5A2) e nós de comunicação (ex: MIF/CD74) para identificar vulnerabilidades computacionais. O mapeamento de drogas foi limitado a conceitos baseados em conhecimento devido à falta de dados de IC50 classificáveis.

Principais Resultados

1. Identificação e Localização do Programa

• O programa APOD-EFEMP1-COL5A2 localiza-se preferencialmente em um compartimento de CAF de alta confiança purificado.
• No atlas de descoberta, 2.415 CAFs de alta confiança foram identificados em 34 amostras. Estas células não estavam confinadas a um único cluster, mas mostraram enriquecimento preferencial em clusters específicos (ex: Cluster 0, OR 3.51).
• Os CAFs do estado "Programa-Alto" exibiram um estado transcricional coerente caracterizado pela regulação positiva de genes de ECM, colágeno, adesão focal e genes associados a vasos sanguíneos.

2. Identidade Transcricional

• A análise de pseudobulk pareada (21 amostras) identificou 261 genes significativos (187 para cima, 74 para baixo).
• As principais vias reguladas positivamente incluíram Organização de fibrila de colágeno (NES 2.44) e Morfogênese de vasos sanguíneos.
• O estado de programa-alto está fortemente associado à remodelagem de matriz, sinalização de colágeno/integrina e características vasculares.

3. Validação Independente

• Célula Única: O programa foi detectável no conjunto de dados independente GSE163558, embora as células de estado alto estivessem altamente concentradas em poucas amostras (uma amostra contribuiu com 75% das células de estado alto). A consistência direcional foi observada para os genes definidores e características de ECM.
• Espacial: Em GSE251950, o índice espacial de três genes mostrou correlações positivas com características de núcleo de CAF, endotélio, angiogênese e mieloide em todas as 10 seções. A análise de vizinhança confirmou o enriquecimento desses recursos ao redor dos spots de alto índice.

4. Associação Clínica

• Sobrevivência: Pontuações mais altas do programa foram significativamente associadas a desfechos piores.
  • TCGA-STAD & GSE15459 (OS): A meta-análise resultou em um Hazard Ratio (HR) agrupado de 1.80 (95% CI 1.40–2.33, P = 6.24×10⁻⁶).
  • GSE26253 (RFS): Tumores com programa-alto tiveram menor sobrevida livre de recorrência (HR 1.59, 95% CI 1.18–2.14, P = 0.0025).
• Fenótipo do TME: O programa correlacionou-se fortemente com assinaturas de núcleo de CAF, ECM/remodelagem de matriz e exclusão imune estromal em todas as 10 coortes de bulk.

5. Modelagem Prognóstica

• O programa forneceu informações prognósticas modestas, porém reprodutíveis.
• Adicionar a pontuação do programa a variáveis clínicas (ex: idade) melhorou o C-index em 0.060 no TCGA-STAD e em 0.103 no GSE15459.
• O programa adicionou valor incremental além de características genéricas de CAF (ex: escores de núcleo de CAF), com valores positivos de delta C-index através das coortes.

6. Hipóteses Mecanísticas

• Ramo Matriz-Integrina-Vascular: CAFs de Programa-Alto estão computacionalmente ligados a ligantes de matriz/colágeno (ex: COL1A1, COL5A2) e receptores de integrina (ex: ITGB1) em células endoteliais e pericitos/células musculares lisas.
• Ramo SP1-MIF/CD74-Mieloide: SP1 foi identificado como um principal candidato regulador a montante (upstream). O programa está ligado a um eixo de comunicação MIF-CD74, com altos escores de produto de expressão para receptores mieloides/macrófagos.
• Perturbação Virtual: O nocaute de COL5A2 mostrou a reversão mais forte prevista da remodelagem de ECM e da sinalização de colágeno-integrina. A simulação de bloqueio de MIF/CD74 priorizou receptores mieloides/macrófagos.

Significância e Alegações

O artigo afirma estabelecer um programa estromal biologicamente coerente e reprodutível que liga a biologia de matriz de CAF, a organização tecidual associada a vasos e o contexto mieloide a desfechos adversos no câncer gástrico.

• Limites de Evidência: O estudo declara explicitamente que não estabelece um biomarcador universal, um ensaio clínico pronto ou um mecanismo causal. Os achados são descritos como "hipóteses computacionais" e "evidências de suporte", em vez de prova definitiva.
• Proposta de Valor: A principal contribuição é um eixo compacto de três genes que serve como um controle estável para um fenótipo estromal mais amplo e complexo. Ele consegue unir a localização de célula única, o contexto espacial e os desfechos clínicos de bulk sem confundir unidades biológicas distintas.
• Potencial Translacional: O estudo prioriza hipóteses específicas testáveis para acompanhamento experimental:
  1. Validar o ramo de comunicação matriz-integrina-vascular.
  2. Testar o eixo regulatório e de comunicação SP1-MIF/CD74-mieloide.
  3. Investigar COL5A2 e MIF/CD74 como potenciais nós de intervenção.
• Limitações: Os autores reconhecem limitações, incluindo a natureza retrospectiva dos dados públicos, a concentração de amostras nos conjuntos de validação, os limites de resolução de nível de spot nos dados espaciais e a natureza computacional das perturbações e do mapeamento de drogas. Não são feitas previsões formais de eficácia de drogas ou resposta ao tratamento.

Em resumo, o artigo fornece um rigoroso arcabouço computacional multicamadas que define um nicho associado a CAFs adversos no câncer gástrico, oferecendo um conjunto focado de hipóteses para futura validação experimental em relação à remodelagem de matriz, organização vascular e interações mieloides.