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

Problematica

La prognosi del cancro gastrico è fortemente influenzata dal microambiente tumorale (TME), in particolare dall'organizzazione della matrice extracellulare (ECM), dall'architettura vascolare e dalla localizzazione delle cellule immunitarie da parte dei fibroblasti associati al cancro (CAF). Tuttavia, l'attuale comprensione clinica è limitata perché le misurazioni dell'abbondanza stromale "bulk" collassano stati fibroblastici biologicamente distinti in un unico segnale. Manca la presenza di marcatori compatti e riproducibili che colleghino specifici stati trascrizionali dei fibroblasti al contesto tissutale e agli esiti clinici senza incorrere in overfitting su specifici dataset o artefatti tecnici. Lo studio mira a definire un programma fibroblastico specifico e compatto che possa essere tracciato attraverso piattaforme di trascrittomica single-cell, spaziale e bulk per identificare una nicchia legata a esiti avversi.

Metodologia

Lo studio impiega una gerarchia di evidenze multi-livello e cross-platform per valutare un programma a tre geni APOD-EFEMP1-COL5A2. L'analisi distingue rigorosamente tra unità biologiche (cellule, spot, campioni, coorti) per evitare la pseudo-replicazione.

1. Atlante Single-Cell di Scoperta (GSE183904):

   • Integrazione di 40 campioni di cancro gastrico (156.111 cellule controllate per qualità).
   • Identificazione di otto compartimenti principali. I fibroblasti sono stati purificati utilizzando lo scoring UCell per i geni di identità CAF (es. DCN, LUM, FAP) ed escludendo le firme di contaminazione (epiteliali, plasmatiche, immunitarie, endoteliali).
   • Definizione di 2.415 CAF ad alta confidenza da 36 campioni.
   • Definizione del Programma: All'interno di questo set di CAF ad alta confidenza, le cellule sono state stratificate in stati "High-Programme", "Middle" e "Low" basandosi sui quartili per campione del mean z-score di APOD, EFEMP1 e COL5A2.

2. Caratterizzazione Trascrizionale:

   • Analisi Pseudobulk Accoppiata: Aggregazione dei conteggi grezzi da 21 campioni contenenti coppie di stati CAF high e low. È stato utilizzato edgeR con un design accoppiato per identificare geni differenzialmente espressi (DEG) e per eseguire l'arricchimento funzionale (GO, KEGG, Hallmark).
   • Validazione Single-Cell Indipendente (GSE163558): Ricalcolo dello score del programma all'interno di un separato oggetto CAF (255 cellule, 6 campioni) per valutare la rilevabilità e la direzionalità, annotando i limiti di concentrazione del campione.

3. Trascrittomica Spaziale (GSE251950 & GSE246011):

   • Analisi di 10 sezioni tissutali (31.202 spot) per testare la colocalizzazione spaziale.
   • Calcolo di un indice di espressione spaziale a tre geni (media dei geni z-trasformati) e correlazione con i feature score di CAF-core, endotelio, angiogenesi e mieloide all'interno di ogni sezione.
   • Esecuzione di un'analisi di arricchimento del vicinato tramite test di permutazione per valutare l'associazione locale con specifiche caratteristiche del TME.

4. Analisi Trascrittomica Bulk e Clinica (10 Coorti):

   • Scoring di 10 coorti bulk indipendenti (incluse TCGA-STAD, GSE15459, GSE26253) utilizzando la formula del programma.
   • Analisi della Sopravvivenza: Stratificazione dei pazienti per mediana del punteggio del programma. Eseguita regressione di Cox a rischi proporzionali per la Sopravvivenza Globale (OS) e la Sopravvivenza Libera da Recidiva (RFS). Eseguita meta-analisi a effetti casuali per l'OS.
   • Correlazione TME: Correlazione dei punteggi del programma con le caratteristiche stromali, di esclusione immunitaria e vascolari all'interno di ogni coorte prima della sintesi.

5. Inferenza Meccanicistica e di Perturbazione:

   • Rete Regolatoria: Classificazione dei fattori di trascrizione (TF) basata su punteggi integrati che collegano l'espressione dei TF all'attività del programma.
   • Analisi Ligando-Recettore: Utilizzo di uno score di prodotto di espressione (ligando emittente × recettore ricevente) per dare priorità alle branche di comunicazione, focalizzandosi specificamente sulle vie matrice-integrina e MIF/CD74.
   • Perturbazione Virtuale: Simulazione di knockdown dei geni del programma (es. COL5A2) e dei nodi di comunicazione (es. MIF/CD74) per identificare vulnerabilità computazionali. La mappatura dei farmaci è stata limitata a concetti basati sulla conoscenza a causa della mancanza di dati IC50 classificabili.

Risultati Chiave

1. Identificazione e Localizzazione del Programma

• Il programma APOD-EFEMP1-COL5A2 si localizza preferenzialmente in un compartimento CAF purificato ad alta confidenza.
• Nell'atlante di scoperta, sono stati identificati 616 CAF "Programme-High" in 34 campioni. Queste cellule non erano confinate a un singolo cluster ma mostravano un arricchimento preferenziale in cluster specifici (es. Cluster 0, OR 3.51).
• I CAF "Programme-High" esibivano uno stato trascrizionale coerente caratterizzato dall'upregulation di geni relativi a ECM, collagene, adesione focale e geni associati ai vasi.

2. Identità Trascrizionale

• L'analisi pseudobulk accoppiata (21 campioni) ha identificato 261 geni significativi (187 up, 74 down).
• Le principali vie upregulated includevano l'Organizzazione delle fibrille di collagene (NES 2.44) e la Morfogenesi dei vasi sanguigni.
• Lo stato "Programme-High" è fortemente associato al rimodellamento della matrice, alla segnalazione di collagene/integrina e a caratteristiche vascolari.

3. Validazione Indipendente

• Single-Cell: Il programma era rilevabile nel dataset indipendente GSE163558, sebbene le cellule ad alto stato fossero altamente concentrate in pochi campioni (un campione ha contribuito con il 75% delle cellule ad alto stato). È stata osservata una coerenza direzionale per i geni definitori e le caratteristiche ECM.
• Spaziale: In GSE251950, l'indice spaziale a tre geni ha mostrato correlazioni positive con le caratteristiche di CAF-core, endotelio, angiogenesi e mieloide in tutte le 10 sezioni. L'analisi del vicinato ha confermato l'arricchimento di queste caratteristiche attorno agli spot ad alto indice.

4. Associazione Clinica

• Sopravvivenza: Punteggi del programma più elevati sono stati significativamente associati a esiti peggiori.
  • TCGA-STAD & GSE15459 (OS): La meta-analisi ha prodotto un Hazard Ratio (HR) combinato di 1.80 (95% CI 1.40–2.33, P = 6.24×10⁻⁶).
  • GSE26253 (RFS): I tumori con programma-alto hanno presentato una sopravvivenza libera da recidiva più breve (HR 1.59, 95% CI 1.18–2.14, P = 0.0025).
• Fenotipo TME: Il programma correla fortemente con le firme di CAF-core, rimodellamento ECM/matrice e esclusione immunitaria stromale in tutte le 10 coorti bulk.

5. Modellazione Prognostica

• Il programma ha fornito informazioni prognostiche modeste ma riproducibili.
• L'aggiunta dello score del programma alle variabili cliniche (es. età) ha migliorato il C-index di 0.060 in TCGA-STAD e di 0.103 in GSE15459.
• Il programma ha aggiunto un valore incrementale rispetto alle generiche caratteristiche CAF (es. score CAF-core), con valori positivi di delta C-index attraverso le coorti.

6. Ipotesi Meccanicistiche

• Branca Matrice-Integrina-Vascolare: I CAF "Programme-High" sono computazionalmente collegati a ligandi della matrice/collagene (es. COL1A1, COL5A2) e recettori dell'integrina (es. ITGB1) su cellule endoteliali e periciti/cellule muscolari lisce.
• Branca SP1-MIF/CD74-Mieloide: SP1 è stato identificato come un top candidato regolatore a monte. Il programma è collegato a un asse di comunicazione MIF-CD74, con alti score di prodotto di espressione per i riceventi mieloidi/macrofagici.
• Perturbazione Virtuale: Il knockdown di COL5A2 ha mostrato la più forte inversione prevista del rimodellamento dell'ECM e della segnalazione collagene-integrina. La simulazione del blocco di MIF/CD74 ha dato priorità ai riceventi mieloidi/macrofagici.

Significatività e Claim

Il paper sostiene di aver stabilito un programma stromale biologicamente coerente e riproducibile che collega la biologia della matrice dei CAF, l'organizzazione tissutale associata ai vasi e il contesto mieloide agli esiti avversi del cancro gastrico.

• Confini dell'Evidenza: Lo studio dichiara esplicitamente di non stabilire un biomarcatore universale, un saggio clinico pronto all'uso o un meccanismo causale. I risultati sono descritti come "ipotesi computazionali" e "evidenze di supporto" piuttosto che prove definitive.
• Proposta di Valore: Il contributo primario è un asse compatto di tre geni che funge da stabile maniglia per un fenotipo stromale più ampio e complesso. Esso riesce a collegare la localizzazione single-cell, il contesto spaziale e gli esiti clinici bulk senza confondere distinte unità biologiche.
• Potenziale Traslazionale: Lo studio dà priorità a specifiche ipotesi testabili per il follow-up sperimentale:
  1. Validare la branca di comunicazione matrice-integrina-vascolare.
  2. Testare l'asse regolatorio e di comunicazione SP1-MIF/CD74-mieloide.
  3. Indagare COL5A2 e MIF/CD74 come potenziali nodi di intervento.
• Limitazioni: Gli autori riconoscono le limitazioni, tra cui la natura retrospettiva dei dati pubblici, la concentrazione dei campioni nei set di validazione, i limiti di risoluzione degli spot nei dati spaziali e la natura computazionale dei risultati di perturbazione e mappatura dei farmaci. Non vengono rivendicate previsioni formali di efficacia dei farmaci o di risposta al trattamento.

In sintesi, il paper fornisce un rigoroso framework computazionale multi-livello che definisce una nicchia associata a CAF con esiti avversi nel cancro gastrico, offrendo un set focalizzato di ipotesi per la futura validazione sperimentale riguardanti il rimodellamento della matrice, l'organizzazione vascolare e le interazioni mieloidi.