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

Énoncé du problème

Le pronostic du cancer gastrique est fortement influencé par le microenvironnement tumoral (MET), en particulier par l'organisation de la matrice extracellulaire (MEC), l'architecture vasculaire et la localisation des cellules immunitaires par les fibroblastes associés au cancer (CAF). Cependant, la compréhension clinique actuelle est limitée car les mesures d'abondance stromale globale fusionnent des états de fibroblastes biologiquement distincts en un signal unique. Il existe un manque de marqueurs compacts et reproductibles capables de relier des états transcriptionnels spécifiques de fibroblastes au contexte tissulaire et aux résultats cliniques sans surapprentissage (overfitting) sur des jeux de données spécifiques ou des artefacts techniques. L'étude vise à définir un programme de fibroblastes spécifique et compact qui puisse être tracé à travers les plateformes transcriptomiques de cellule unique (single-cell), spatiales et globales (bulk) afin d'identifier une niche liée à des issues défavorables.

Méthodologie

L'étude emploie une hiérarchie de preuves multicouche et multi-plateformes pour évaluer un programme de trois gènes : APOD-EFEMP1-COL5A2. L'analyse distingue strictement les unités biologiques (cellules, spots, échantillons, cohortes) pour éviter la pseudo-réplication.

1. Atlas de cellule unique de découverte (GSE183904) :

   • Intégration de 40 échantillons de cancer gastrique (156 111 cellules contrôlées par la qualité).
   • Identification de huit compartiments majeurs. Les fibroblastes ont été purifiés en utilisant le score UCell pour les gènes d'identité des CAF (ex. : DCN, LUM, FAP) et l'exclusion des signatures de contamination (épithéliales, plasmocytaires, immunitaires, endothéliales).
   • Définition de 2 415 CAFs de haute confiance issus de 36 échantillons.
   • Définition du programme : Au sein de cet ensemble de CAFs de haute confiance, les cellules ont été stratifiées en états « Programme-Haut », « Moyen » et « Bas » basés sur les quartiles par échantillon du z-score moyen de APOD, EFEMP1 et COL5A2.

2. Caractérisation transcriptionnelle :

   • Analyse de pseudobulk appariée : Agrégation des comptes bruts de 21 échantillons contenant des états de CAF élevés et bas appariés. Utilisation d'edgeR avec un plan apparié pour identifier les gènes différentiellement exprimés (DEG) et réaliser l'enrichissement fonctionnel (GO, KEGG, Hallmark).
   • Validation par cellule unique indépendante (GSE163558) : Recalcul du score du programme dans un objet CAF séparé (255 cellules, 6 échantillons) pour évaluer la détectabilité et la directionnalité, en notant les limitations de concentration d'échantillons.

3. Transcriptomique spatiale (GSE251950 & GSE246011) :

   • Analyse de 10 sections de tissus (31 202 spots) pour tester la colocalisation spatiale.
   • Calcul d'un indice d'expression spatiale à trois gènes (moyenne des gènes transformés en z) et corrélation avec les scores de caractéristiques des noyaux de CAF (CAF-core), endothéliales, de l'angiogenèse et myéloïdes au sein de chaque section.
   • Réalisation d'une analyse d'enrichissement de voisinage par tests de permutation pour évaluer l'association locale avec des caractures spécifiques du MET.

4. Transcriptomique globale et analyse clinique (10 cohortes) :

   • Calcul du score pour 10 cohortes globales indépendantes (incluant TCGA-STAD, GSE15459, GSE26253) en utilisant la formule du programme.
   • Analyse de survie : Stratification des patients par la médiane du score du programme. Réalisation de régressions de Cox à risques proportionnels pour la survie globale (OS) et la survie sans récidive (RFS). Réalisation d'une méta-analyse à effets aléatoires pour l'OS.
   • Corrélation avec le MET : Corrélation des scores du programme avec les caractéristiques stromales, d'exclusion immunitaire et vasculaires au sein de chaque cohorte avant la synthèse.

5. Inférence de mécanisme et de perturbation :

   • Réseau de régulation : Classement des facteurs de transcription (TF) basés sur des scores intégrés reliant l'expression des TF à l'activité du programme.
   • Analyse Ligand-Récepteur : Utilisation d'un score de produit d'expression (ligand émetteur × récepteur récepteur) pour prioriser les branches de communication, en se concentrant spécifiquement sur les voies matrice-intégrine et MIF/CD74.
   • Perturbation virtuelle : Simulation de l'inhibition (knockdown) des gènes du programme (ex. : COL5A2) et des nœuds de communication (ex. : MIF/CD74) pour identifier les vulnérabilités computationnelles. La cartographie des médicaments a été limitée aux concepts basés sur la connaissance en raison de l'absence de données d'IC50 classables.

Résultats clés

1. Identification et localisation du programme

• Le programme APOD-EFEMP1-COL5A2 se localise préférentiellement dans un compartiment de CAFs de haute confiance purifié.
• Dans l'atlas de découverte, 616 CAFs « Programme-Haut » ont été identifiés à travers 34 échantillons. Ces cellules n'étaient pas confinées à un seul cluster mais montraient un enrichissement préférentiel dans des clusters spécifiques (ex. : Cluster 0, OR 3,51).
• Les CAFs « Programme-Haut » présentaient un état transcriptionnel cohérent caractérisé par une régulation positive des gènes de la matrice extracellulaire, du collagène, des adhérences focales et des gènes associés aux vaisseaux.

2. Identité transcriptionnelle

• L'analyse de pseudobulk appariée (21 échantillons) a identifié 261 gènes significatifs (187 en hausse, 74 en baisse).
• Les voies clés régulées à la hausse incluaient l'organisation des fibrilles de collagène (NES 2,44) et la morphogenèse des vaisseaux sanguins.
• L'état « Programme-Haut » est fortement associé au remodelage de la matrice, à la signalisation collagène/intégrine et aux caractéristiques vasculaires.

3. Validation indépendante

• Cellule unique : Le programme était détectable dans le jeu de données indépendant GSE163558, bien que les cellules à haut niveau soient très concentrées dans quelques échantillons (un échantillon a contribué à 75 % des cellules à haut niveau). Une cohérence directionnelle a été observée pour les gènes définissants et les caractéristiques de la MEC.
• Spatial : Dans GSE251950, l'indice spatial à trois gènes a montré des corrélations positives avec les caractéristiques des noyaux de CAF (CAF-core), endothéliales, de l'angiogenèse et myéloïdes à travers les 10 sections. L'analyse de voisinage a confirmé l'enrichissement de ces caractéristiques autour des spots à haut indice.

4. Association clinique

• Survie : Des scores de programme plus élevés étaient significativement associés à des issues plus défavorables.
  • TCGA-STAD & GSE15459 (OS) : La méta-analyse a produit un Hazard Ratio (HR) combiné de 1,80 (IC 95 % 1,40–2,33, P = 6,24×10⁻⁶).
  • GSE26253 (RFS) : Les tumeurs à haut programme présentaient une survie sans récidive plus courte (HR 1,59, IC 95 % 1,18–2,14, P = 0,0025).
• Phénotype du MET : Le programme est corrélé fortement avec les signatures de noyau de CAF (CAF-core), de remodelage de la MEC/matrice et d'exclusion immunitaire stromale dans les 10 cohortes globales.

5. Modélisation pronostique

• Le programme a fourni des informations pronostiques modestes mais reproductibles.
• L'ajout du score du programme aux variables cliniques (ex. : âge) a amélioré le C-index de 0,060 dans TCGA-STAD et de 0,103 dans GSE15459.
• Le programme apporte une valeur incrémentale au-delà des caractéristiques génériques des CAF (ex. : scores de noyau de CAF), avec des valeurs delta C-index positives à travers les cohortes.

6. Hypothèses mécanistiques

• Branche Matrice-Intégrine-Vasculaire : Les CAFs « Programme-Haut » sont informatiquement liés aux ligands de matrice/collagène (ex. : COL1A1, COL5A2) et aux récepteurs d'intégrine (ex. : ITGB1) sur les cellules endothéliales et les cellules péricytes/muscle lisse.
• Branche SP1-MIF/CD74-Myéloïde : SP1 a été identifié comme un régulateur en amont candidat de premier rang. Le programme est lié à un axe de communication MIF-CD74, avec des scores de produit d'expression élevés pour les récepteurs myéloïdes/macrophages.
• Perturbation virtuelle : L'inhibition de COL5A2 a montré la plus forte inversion prédite du remodelage de la MEC et de la signalisation collagène-intégrine. La simulation du blocage de MIF/CD74 a priorisé les récepteurs myéloïdes/macrophages.

Signification et revendications

L'article affirme établir un programme stromal reproductible et biologiquement cohérent qui lie la biologie de la matrice des CAF, l'organisation tissulaire associée aux vaisseaux et le contexte myéloïde à des issues défavorables du cancer gastrique.

• Limites de preuve : L'étude stipule explicitement qu'elle n'établit pas un biomarqueur universel, un test clinique prêt à l'emploi, ou un mécanisme causal. Les résultats sont décrits comme des « hypothèses computationnelles » et des « preuves de soutien » plutôt que des preuves définitives.
• Proposition de valeur : La contribution primaire est un axe compact de trois gènes servant de levier stable pour un phénotype stromal plus large et complexe. Elle réussit à relier la localisation en cellule unique, le contexte spatial et les résultats cliniques globaux sans confondre des unités biologiques distinctes.
• Potentiel translationnel : L'étude priorise des hypothèses testables spécifiques pour un suivi expérimental :
  1. Valider la branche de communication matrice-intégrine-vasculaire.
  2. Tester l'axe de régulation et de communication SP1-MIF/CD74-myéloïde.
  3. Étudier COL5A2 et MIF/CD74 comme potentiels nœuds d'intervention.
• Limitations : Les auteurs reconnaissent les limitations, notamment la nature rétrospective des données publiques, la concentration des échantillons dans les ensembles de validation, la limite de résolution des spots dans les données spatiales, et la nature computationnelle des résultats de perturbation et de cartographie des médicaments. Aucune prédiction d'efficacité médicamenteuse ou de réponse au traitement n'est revendiquée.

En résumé, l'article fournit un cadre computationnel rigoureux et multicouche qui définit une niche associée aux CAFs défavorables dans le cancer gastrique, offrant un ensemble focalisé d'hypothèses pour une validation expérimentale future concernant le remodelage de la matrice, l'organisation vasculaire et les interactions myéloïdes.