Single-section spatial hypoxia-cytotoxic associations do not consistently reproduce across breast cancer patients

Technische Samenvatting: Ruimtelijke Hypoxie–Cytotoxische Associaties in Enkele Secties Reproduceren Niet Consistent Over Breaschtkankerpatiënten

Probleemstelling
Ruimtelijke transcriptomics maakt het mogelijk om gelokaliseerde tumor–immuunrelaties binnen intact weefsel te bestuderen. Er bestaat echter een kritieke statistische uitdaging: de duizenden meetpunten (spots) binnen een enkele weefselsectie zijn ruimtelijk afhankelijke observaties van één enkel biologisch specimen, en geen onafhankelijke biologische replica's. Het behandelen van deze spots als onafhankelijke bewijslast voor patiëntniveau-claims vormt "pseudoreplicatie". Een sterke, statistisch significante associatie die op het niveau van de spots wordt waargenomen binnen één sectie, weerspiegelt daarom mogelijk geen consistent biologisch effect over verschillende patiënten heen. Deze studie adresseert de kloof tussen ruimtelijke associaties binnen een sectie en de reproduceerbaarheid op patiëntniveau in borstkanker.

Methodologie
De auteurs hanteerden een twee-fasen secundaire-analyse-ontwerp met behulp van publieke borstkankerdatasets:

1. Discovery-fase (Ontdekkingsfase):

   • Dataset: Een enkele 10x Genomics Visium borstkankersectie (3.798 spots).
   • Genprogramma's: Zes vooraf gespecificeerde transcriptionele programma's werden berekend: tumor epitheel, hypoxie (6 genen), cytotoxische lymfocyt (8 genen), myeloïde, fibroblast/ECM, en EMT-regulerend.
   • Analyse: De auteurs construeerden een gesymmetriseerde zes-naaste-buren grafiek om ruimtelijke buurten te definiëren. Ze beoordeelden de ruimtelijke autocorrelatie (Moran's I) en evalueerden de associatie tussen hypoxie-scores en de gemiddelde cytotoxische scores van naburige spots met behulp van Spearman-correlatie en rank-biserial effectgroottes (het vergelijken van hoog- versus laag-hypoxische spots).

2. Externe Validatiefase:

   • Dataset: Een onafhankelijke HER2-positieve borstkanker-cohort (Andersson et al.) bestaande uit 36 secties van acht verschillende patiënten (13.619 totalen spots).
   • Protocol: Elke sectie werd onafhankelijk verwerkt (library-size normalisatie, log-transformatie, gen-gewijze standaardisatie). De directionele associatie tussen hypoxie en buurt-cytotoxiciteit werd voor elke sectie berekend.
   • Aggregatie: Sectie-niveau correlaties werden op patiëntniveau geaggregeerd met behulp van de mediaan.
   • Statistische Toetsing: Een eenzijdige exacte Wilcoxon signed-rank test bepaalde of de patiënt-niveau mediane correlaties significant minder dan nul waren.
   • Gevoeligheidsanalyses: De auteurs testten de robuustheid door het variëren van buurt-groottes (4, 8, 10 buren), het gebruik van een effector-alleen cytotoxische signatuur, het toepassen van Kendall-correlatie, en het substitueren van de compacte hypoxie-signatuur met de 200-gen MSigDB Hallmark hypoxie set.

Belangrijkste Resultaten

• Discovery Sectie: Er werd een duidelijke inverse associatie waargenomen. Hypoxie-gerelateerde transcriptie was negatief geassocieerd met de naburige cytotoxische genactiviteit (Spearman $\rho = -0,202$). Hoog-hypoxische spots hadden significant lagere buurt-cytotoxische scores dan laag-hypoxische spots (rank-biserial effect = $-0,286$).
• Validatie op Patiëntniveau: De inverse associatie reproduceerde zich niet consistent over de externe cohort.
  • Slechts 19 van de 36 secties en 5 van de 8 patiënten vertoonden negatieve associaties.
  • De mediane patiënt-niveau correlatie was $-0,043$, wat niet statistisch onderscheidbaar was van nul (eenzijdige exacte Wilcoxon $P = 0,473$).
  • Patiënt-niveau heterogeniteit was aanzienlijk; bijvoorbeeld vertoonde Patiënt B uniform positieve sectie-niveau associaties, terwijl Patiënt A uniform negatieve associaties vertoonde.
• Gevoeligheidsanalyses:
  • Het variëren van buurt-groottes en het gebruik van Kendall-correlatie leverde zwakke, heterogene patiënt-effecten op (medianen variërend van $-0,025$ tot $-0,032$), waarvan er geen significant waren.
  • Het gebruik van de MSigDB Hallmark hypoxie-signatuur draaide de trend om, wat een positieve mediane patiënt-correlatie produceerde ($0,047$), wat aangeeft dat de geobserveerde richting sterk afhankelijk is van de specifieke transcriptionele definitie van "hypoxie".

Belangrijkste Bijdragen

• Demonstratie van de Risico's van Pseudoreplicatie: De studie levert empirisch bewijs dat significante resultaten op de spot-niveau binnen een enkele ruimtelijke transcriptomics-sectie geen garantie bieden voor biologische consistentie op patiëntniveau.
• Methodologisch Kader: Het vestigt een rigoureuze twee-fasen workflow voor ruimtelijke transcriptomics die de patiënt behandelt als de eenheid van replicatie, waardoor het samenvoegen van ruimtelijk afhankelijke data om de statistische kracht te inflateren wordt voorkomen.
• Karakterisering van Heterogeniteit: De resultaten benadrukken dat ruimtelijke tumor-immuunrelaties (specifiek hypoxie-cytotoxiciteit) zeer heterogeen zijn tussen patiënten en gevoelig zijn voor de keuze van gen-signaturen en buurt-definities.

Betekenis en Claims
Het artikel claimt dat een sterke associatie op sectie-niveau tussen hypoxie en cytotoxische exclusie geen algemeen geldende regel is voor borstkankerpatiënten. In plaats daarvan ondersteunt het bewijs ruimtelijke context-afhankelijkheid in plaats van een universeel hypoxie-geassocieerd cytotoxisch-exclusiepatroon.

De auteurs waarschuwen tegen het interpreteren van ruimtelijke associaties op de spot-niveau van een enkele sectie als patiënt-niveau biologische effecten zonder externe validatie. Zij stellen dat patiënt-heterogeniteit moet worden beschouwd als een biologisch resultaat in plaats van een technische hinderlaag. De studie concludeert dat toekomstig onderzoek naar ruimtelijke transcriptomics prioriteit moet geven aan het toevoegen van onafhankelijke patiënten en het vooraf specificeren van patiënt-niveau estimanden, in plaats van simpelweg het aantal spots binnen een enkele sectie te vergroten. De bevindingen suggereren dat alleen op transcriptoom gebaseerde ruimtelijke associaties mechanismen zoals metabole competitie of verschillen in cel abundantie of weefselarchitectuur niet kunnen onderscheiden zonder verdere validatie.