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

Resumo Técnico: Associações Espaciais Hipóxia–Citotóxica de Seção Única Não se Reproduzem Consistentemente entre Pacientes com Câncer de Mama

Problema
A transcriptômica espacial permite o estudo de relações localizadas entre tumor e sistema imune em tecidos intactos. No entanto, existe um desafio estatístico crítico: os milhares de pontos de medição dentro de uma única seção de tecido são observações espacialmente dependentes de um único espécime biológico, não réplicas biológicas independentes. Tratar esses pontos como evidência independente para alegações ao nível do paciente constitui "pseudoreplicação". Consequentemente, uma associação forte e estatisticamente significativa observada ao nível do ponto dentro de uma única seção pode não refletir um efeito biológico consistente entre diferentes pacientes. Este estudo aborda a lacuna entre as associações espaciais dentro da seção e a reprodutibilidade ao nível do paciente no câncer de mama.

Metodologia
Os autores empregaram um design de análise secundária de dois estágios usando conjuntos de dados públicos de câncer de mama:

1. Fase de Descoberta:

   • Conjunto de Dados: Uma única seção de câncer de mama 10x Genomics Visium (3.798 pontos).
   • Programas Gênicos: Seis programas transcricionais pré-especificados foram calculados: epitelial tumoral, hipóxia (6 genes), linfócito citotóxico (8 genes), mieloide, fibroblasto/ECM e regulador de EMT.
   • Análise: Os autores construíram um grafo de seis vizinhos mais próximos simetrizado para definir vizinhanças espaciais. Eles avaliaram a autocorrelação espacial (I de Moran) e avaliaram a associação entre os escores de hipóxia e a média dos escores citotóxicos de pontos vizinhos usando correlação de Spearman e tamanhos de efeito bisserial de rank (comparando pontos de alta vs. baixa hipóxia).

2. Fase de Validação Externa:

   • Conjunto de Dados: Uma coorte independente de câncer de mama HER2-positivo (Andersson et al.) compreendendo 36 seções de oito pacientes distintos (13.169 pontos totais).
   • Protocolo: Cada seção foi processada independentemente (normalização de tamanho de biblioteca, transformação logarítmica, padronização gene a gene). A associação direcional entre hipóxia e citotoxicidade de vizinhança foi calculada para cada seção.
   • Agregação: As correlações ao nível da seção foram agregadas ao nível do paciente usando a mediana.
   • Teste Estatístico: Um teste de Wilcoxon signed-rank unilateral determinou se as medianas de correlação ao nível do paciente eram significativamente menores que zero.
   • Análises de Sensibilidade: Os autores testaram a robustez variando os tamanhos de vizinhança (4, 8, 10 vizinhos), usando uma assinatura citotóxica apenas de efetores, usando correlação de Kendall e substituindo a assinatura compacta de hipóxia pelo conjunto de 200 genes MSigDB Hallmark hypoxia.

Resultos Principais

• Seção de Descoberta: Uma clara associação inversa foi observada. A transcrição relacionada à hipóxia foi negativamente associada à atividade gênica citotóxica vizinha (Spearman $\rho = -0,202$). Pontos de alta hipóxia apresentaram escores citotóxicos de vizinhança significativamente menores do que pontos de baixa hipóxia (efeito bisserial de rank = $-0,286$).
• Validação ao Nível do Paciente: A associação inversa não se reproduziu consistentemente na coorte externa.
  • Apenas 19 de 36 seções e 5 de 8 pacientes mostraram associações negativas.
  • A mediana da correlação ao nível do paciente foi de $-0,043$, o que não foi estatisticamente distinguível de zero (Wilcoxon signed-rank exato unilateral $P = 0,473$).
  • A heterogeneidade ao nível do paciente foi substancial; por exemplo, o Paciente B apresentou associações de nível de seção uniformemente positivas, enquanto o Paciente A apresentou associações uniformemente negativas.
• Análises de Sensibilidade:
  • Variar os tamanhos de vizinhança e usar a correlação de Kendall resultou em efeitos de pacientes fracos e heterogêneos (medianas variando de $-0,025$a$-0,032$), nenhum dos quais foi significativo.
  • O uso da assinatura MSigDB Hallmark hypoxia inverteu a tendência, produzindo uma mediana de correlação de paciente positiva ($0,047$), indicando que a direção observada é altamente dependente da definição transcricional específica de "hipóxia".

Principais Contribuições

• Demonstração dos Riscos de Pseudoreplicação: O estudo fornece evidência empírica de que a significância ao nível do ponto em uma única seção de transcriptômica espacial não garante a consistência biológica ao nível do paciente.
• Estrutura Metodológica: Estabelece um fluxo de trabalho rigoroso de dois estágios para a transcriptômica espacial que trata o paciente como a unidade de replicação, evitando o agrupamento de dados espacialmente dependentes para inflar o poder estatístico.
• Caracterização da Heterogeneidade: Os resultados destacam que as relações espaciais tumor-imune (especificamente hipóxia-citotoxicidade) são altamente heterogêneas entre pacientes e sensíveis à escolha das assinaturas gênicas e definições de vizinhança.

Significância e Alegações
O artigo afirma que uma forte associação de seção única entre hipóxia e exclusão citotóxica não é uma regra generalizável entre pacientes com câncer de mama. Em vez disso, a evidência apoia a dependência do contexto espacial em vez de um padrão universal de exclusão citotóxica associado à hipóxia.

Os autores alertam contra a interpretação de associações espaciais de nível de ponto de uma única seção como efeitos biológicos ao nível do paciente sem validação externa. Eles argumentam que a heterogeneidade do paciente deve ser tratada como um resultado biológico e não como um incômodo técnico. O estudo conclui que a pesquisa futura em transcriptômica espacial deve priorizar a adição de pacientes independentes e a pré-especificação de estimandos ao nível do paciente, em vez de simplesmente aumentar o número de pontos dentro de uma única seção. Os achados sugerem que associações espaciais baseadas apenas no transcriptoma não podem distinguir mecanismos como competição metabólica de diferenças na abundância celular ou arquitetura tecidual sem validação adicional.