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

Probleemstelling

De prognose van maagkanker wordt sterk beïnvloed door de tumor micro-omgeving (TME), in het bijzonder de organisatie van de extracellulaire matrix (ECM), de vasculaire architectuur en de lokalisatie van immuuncellen door kanker-geassocieerde fibroblasten (CAFs). De huidige klinische kennis is echter beperkt omdat metingen van de bulk stroma-abundantie biologisch onderscheidbare fibroblast-toestanden samenvoegen tot één enkel signaal. Er is een gebrek aan compacte, reproduceerbare markers die specifieke transcriptionele fibroblast-toestanden verbinden met de weefselcontext en klinische uitkomsten zonder te overfitten op specifieke datasets of technische artefacten. De studie beoogt een specifiek, compact fibroblast-programma te definiëren dat traceerbaar is over single-cell, ruimtelijke en bulk transcriptomische platforms om een niche te identificeren die gekoppeld is aan ongunstige uitkomsten.

Methodologie

De studie maakt gebruik van een meerlagige, cross-platform bewijshierarchie om een APOD-EFEMP1-COL5A2 drie-genen programma te evalueren. De analyse maakt een strikt onderscheid tussen biologische eenheden (cellen, spots, monsters, cohorten) om pseudo-replicatie te voorkomen.

1. Ontdekkings Single-Cell Atlas (GSE183904):

   • Geïntegreerde 40 maagkanker monsters (156.111 kwaliteitsgecontroleerde cellen).
   • Identificeerde acht belangrijke compartimenten. Fibroblasten werden gezuiverd met behulp van UCell scoring voor CAF-identiteit genen (bijv. DCN, LUM, FAP) en uitsluiting van contaminatie-signatures (epitheel, plasma, immuun, endotheel).
   • Definieerde 2.415 high-confidence CAFs uit 36 monsters.
   • Programmadefinities: Binnen deze high-confidence CAF-set werden cellen gestratificeerd in "Programme-High", "Middle", en "Low" toestanden op basis van monster-specifieke kwartielen van de gemiddelde z-score van APOD, EFEMP1, en COL5A2.

2. Transcriptionele Karakterisering:

   • Gepaard Pseudobulk Analyse: Geaggregeerde ruwe counts van 21 monsters die gepaard gingen met hoge en lage CAF-toestanden. Gebruikte edgeR met een gepaard ontwerp om differentieel tot expressie gebrachte genen (DEGs) te identificeren en functionele verrijking (GO, KEGG, Hallmark) uit te voeren.
   • Onafhankelijke Single-Cell Validatie (GSE163558): Herrekende de programme score binnen een apart CAF-object (255 cellen, 6 monsters) om detecteerbaarheid en directionaliteit te beoordelen, waarbij werd opgemerkt dat de concentratie van monsters een beperking vormde.

3. Ruimtelijke Transcriptomics (GSE251950 & GSE246011):

   • Analyseerde 10 weefselsecties (31.202 spots) om ruimtelijke colocalisatie te testen.
   • Berekende een drie-genen ruimtelijke expressie-index (gemiddelde van z-getransformeerde genen) en correleerde dit met CAF-core, endotheel, angiogenese en myeloïde feature scores binnen elke sectie.
   • Voerde buurtverrijking-analyse uit met permutatietesten om de lokale associatie met specifieke TME-features te beoordelen.

4. Bulk Transcriptomics & Klinische Analyse (10 Cohorten):

   • Scorde 10 onafhankelijke bulk cohorten (inclusief TCGA-STAD, GSE15459, GSE26253) met de formule van het programma.
   • Overlevingsanalyse: Stratificeerde patiënten op basis van de mediaan programme score. Voerde Cox proportional-hazards regressie uit voor Overall Survival (OS) en Recurrence-Free Survival (RFS). Voerde een random-effects meta-analyse uit voor OS.
   • TME Correlatie: Correleerde programme scores met stroma, immuun-exclusie en vasculaire features binnen elk cohort voordat deze werden samengevat.

5. Mechanisme & Perturbatie Inferentie:

   • Regulatoir Netwerk: Rangschikte transcriptiefactoren (TFs) op basis van geïntegreerde scores die de TF-expressie koppelen aan de programme activiteit.
   • Ligand-Receptor Analyse: Gebruikte een expressie-product score (zender ligand × ontvanger receptor) om communicatievertakkingen te prioriteren, met specifieke focus op matrix-integrine en MIF/CD74 pathways.
   • Virtuele Perturbatie: Simuleerde knockdowns van programme genen (bijv. COL5A2) en communicatieknooppunten (bijv. MIF/CD74) om computationele kwetsbaarheden te identificeren. Drug mapping was beperkt tot kennisgebaseerde concepten vanwege het ontbreken van rankbare IC50-data.

Belangrijkste Resultaten

1. Identificatie en Lokalisatie van het Programma

• Het APOD-EFEMP1-COL5A2 programma lokaliseert zich bij voorkeur in een gezuiverd, high-confidence CAF-compartiment.
• In de ontdekkingsatlas werden 616 Programme-High CAFs geïdentificeerd over 34 monsters. Deze cellen waren niet beperkt tot één enkele cluster maar vertoonden een voorkeursverrijking in specifieke clusters (bijv. Cluster 0, OR 3.51).
• Programme-High CAFs vertoonden een coherente transcriptionele staat gekenmerkt door upregulatie van ECM, collageen, focal adhesion en vasculaire-geassocieerde genen.

2. Transcriptionele Identiteit

• Gepaarde pseudobulk analyse (21 monsters) identificeerde 261 significante genen (187 up, 74 down).
• Belangrijke geüpreguleerde pathways waren onder andere Collagen fibril organization (NES 2.44) en Blood vessel morphogenesis.
• De programme-high staat is sterk geassocieerd met matrix remodeling, collageen/integrine signalering en vasculaire kenmerken.

3. Onafhankelijke Validatie

• Single-Cell: Het programma was detecteerbaar in de onafhankelijke GSE163558 dataset, hoewel de high-state cellen sterk geconcentreerd waren in enkele monsters (één monster droeg 75% bij aan de high-state cellen). Er werd een directionele consistentie waargenomen voor de definiërende genen en ECM-features.
• Spatial: In GSE251950 vertoonde de drie-genen ruimtelijke index positieve correlaties met CAF-core, endotheel, angiogenese en myeloïde features over alle 10 secties. Buurtanalyse bevestigde de verrijking van deze features rondom high-index spots.

4. Klinische Associatie

• Overleving: Hogere programme scores waren significant geassocieerd met slechtere uitkomsten.
  • TCGA-STAD & GSE15459 (OS): De meta-analyse leverde een gepoolde Hazard Ratio (HR) op van 1.80 (95% CI 1.40–2.33, P = 6.24×10⁻⁶).
  • GSE26253 (RFS): Programme-high tumoren hadden een kortere recurrence-free survival (HR 1.59, 95% CI 1.18–2.14, P = 0.0025).
• TME Fenotype: Het programma correleerde sterk met CAF-core, ECM/matrix remodeling en stromale immuun-exclusie signatures in alle 10 de bulk cohorten.

5. Prognostische Modellering

• Het programma bood bescheiden maar reproduceerbare prognostische informatie.
• Het toevoegen van de programme score aan klinische variabelen (bijv. leeftijd) verbeterde de C-index met 0.060 in TCGA-STAD en met 0.103 in GSE15459.
• Het programma voegde incrementele waarde toe bovenop generieke CAF-features (bijv. CAF-core scores), met positieve delta C-index waarden over de cohorten heen.

6. Mechanistische Hypothesen

• Matrix-Integrine-Vasculaire Tak: Programme-High CAFs zijn computationeel gelinkt aan matrix/collageen liganden (bijv. COL1A1, COL5A2) en integrine receptoren (bijv. ITGB1) op endotheel- en pericyt/gladde spiercellen.
• SP1-MIF/CD74-Myeloïde Tak: SP1 werd geïdentificeerd als een top kandidaat upstream regulator. Het programma is gelinkt aan een MIF-CD74 communicatie-as, met hoge expressie-product scores voor myeloïde/macrofaag ontvangers.
• Virtuele Perturbatie: COL5A2 knockdown toonde de sterkste voorspelde omkering van ECM remodeling en collageen-integrine signalering. MIF/CD74 blokkade simulatie gaf prioriteit aan myeloïde/macrofaag ontvangers.

Betekenis en Claims

Het artikel claimt een reproduceerbaar, biologisch coherent stroma programma vast te stellen dat de CAF matrix biologie, vasculaire-geassocieerde weefselorganisatie en de myeloïde context verbindt met ongunstige uitkomsten bij maagkanker.

• Bewijsgrenzen: De studie stelt expliciet dat zij geen universele biomarker, een klaar voor klinische toepassing zijnde assay of een causaal mechanisme vaststelt. De bevindingen worden beschreven als "computationele hypothesen" en "ondersteunend bewijs" in plaats van definitief bewijs.
• Waardepropositie: De primaire bijdrage is een compacte, drie-genen as die dient als een stabiel handvat voor een breder, complex stroma fenotype. Het slaagt erin om single-cell lokalisatie, ruimtelijke context en bulk klinische uitkomsten te overbruggen zonder verschillende biologische eenheden te verwarren.
• Translationeel Potentieel: De studie geeft prioriteit aan specifieke testbare hypothesen voor experimentele opvolging:
  1. Validatie van de matrix-integrine-vasculaire communicatie tak.
  2. Testen van de SP1-MIF/CD74-myeloïde regulatoire en communicatieve as.
  3. Onderzoek naar COL5A2 en MIF/CD74 als potentiële interventie knooppunten.
• Beperkingen: De auteurs erkennen beperkingen zoals het retrospectieve karakter van publieke data, de concentratie van monsters in validatiesets, de resolutiebeperkingen op spot-niveau in ruimtelijke data, en het computationele karakter van de perturbatie en drug-mapping resultaten. Er worden geen formele claims gemaakt over effectiviteit van medicijnen of voorspellingen van behandelrespons.

Samenvattend biedt het artikel een rigoureus, meerlagig computationeel kader dat een ongunstig CAF-geassocieerd niche definieert in maagkanker, en biedt het een gefocuste set hypothesen voor toekomstige experimentele validatie met betrekking tot matrix remodeling, vasculaire organisatie en myeloïde interacties.