Stable but confounded: brain region, not diagnosis, drives a validation-passing molecular subtype in postmortem psychiatric brain

Resumo Técnico: "Estável, mas confundido: região cerebral, não o diagnóstico, impulsiona um subtipo molecular que passa na validação em cérebro psiquiátrico postmortem"

Problema
A subtipagem molecular não supervisionada de transcriptomas cerebrais é cada vez mais utilizada para resolver a heterogeneidade de diagnósticos psiquiátricos (ex: autismo, esquizofrenia). No entanto, algoritmos de agrupamento (clustering) sempre retornarão uma partição, independentemente de essa partição refletir uma biologia reproduzível ou meramente artefatos técnicos/anatômicos. O artigo identifica três perigos agudos em estudos psiquiátricos postmortem:

1. Tamanhos de amostra pequenos (dezenas a poucas centenas) que limitam o poder estatístico.
2. Sobreposição de sujeitos entre coortes nominalmente independentes, levando a falsas alegações de replicação.
3. Dominância anatômica, onde a variação transcriptômica por região cerebral excede vastamente a variação por diagnóstico. Consequentemente, o agrupamento de dados multi-regiões combinados frequentemente recupera a anatomia em vez de subtipos de doenças.

Os autores argumentam que as práticas de validação atuais frequentemente falham em distinguir entre um cluster "estável" (estrutura reproduzível) e um cluster "válido" (biologicamente relevante para a doença), particularmente quando confundidores como a região cerebral não são explicitamente testados.

Metodologia
Os autores reanalisaram três coortes públicas postmortem (384 amostras no total) usando uma disciplina de "validação primeiro", que reúne controles estabelecidos e adiciona testes explícitos de confundidores. Os conjuntos de dados incluíram:

• Transtorno do Espectro Autista (TEA): GSE28521 (microarray, córtex frontal/temporal/cerebelo) e GSE64018 (RNA-seq, giro temporal superior).
• Esquizofrenia, Transtorno Bipolar, Depressão Maior: GSE80655 (RNA-seq, cingulado anterior, DLPFC, núcleo accumbens).

Pipeline Analítico:

1. Redução de Características: A expressão gênica foi resumida em escores de nível de via (pathway-level scores) usando enriquecimento de amostra única (ssGSEA ou média de z-scores por gene).
2. Agrupamento (Clustering): As amostras foram agrupadas usando Modelagem de Mistura Gaussiana (GMM) com seleção de componentes via Critério de Informação Bayesiano (BIC).
3. Controles de Validação: Cada partição candidata foi submetida a quatro etapas:
   • Nulo de permutação de rótulos: Testando se os clusters são recuperáveis a partir de rótulos permutados.
   • Nulo de conjunto de genes aleatórios: Testando se a partição é específica para vias curadas ou se surge de conjuntos de genes aleatórios.
   • Estabilidade de Bootstrap: Medindo o Índice Rand Ajustado (ARI) entre as atribuições originais e reamostradas (limiar $\ge$ 0,80).
   • Independência de ancestralidade: Verificando a associação com componentes principais de ancestralidade.
4. Teste de Confundidores: Foi realizada a associação explícita por tabela de contingência ($\chi^2$, V de Cramér) entre as partições e a região cerebral, sobreposição de doadores e diagnóstico clínico.
5. Análise Não Circular: Para evitar o "double dipping" (uso duplo de dados), os sinais de doença foram avaliados via contrastes caso-controle baseados no diagnóstico dentro das regiões, utilizando um painel fixo de genes marcadores de tipos celulares canônicos (ex: PVALB, SST, GAD1), em vez de características derivadas do próprio agrupamento.

Resultos Principais

1. Falha do Subtipo "GABAérgico" do Autismo:
   Um suposto subgrupo de TEA caracterizado pela regulação negativa da via GABAérgica falhou na validação de estabilidade de bootstrap (ARI = 0,71, abaixo do limiar de 0,80). Além disso, a aparente replicação em uma segunda coorte (GSE64018) foi invalidada porque as duas coortes compartilhavam 10 de 24 sujeitos, significando que o acordo era uma corroboração entre plataformas sobre os mesmos cérebros, e não uma replicação independente.

2. A Partição de Esquizofrenia que "Passa na Validação" é um Classificador de Região:
   Na coorte de esquizofrenia (GSE80655), uma partição $k=3$ passou por todos os portões estatísticos (Bootstrap ARI = 0,923; passou nos testes de conjunto de genes aleatórios e ancestralidade). No entanto, o teste de confundidores revelou que esta partição era impulsionada inteiramente pela região cerebral, não pelo diagnóstico:

   • O Subtipo 0 e 2 eram predominantemente núcleo accumbens (nAcc).
   • O Subtipo 1 era exclusivamente cortical (AnCg e DLPFC).
   • A associação entre o subtipo e a região foi esmagadora ($p = 4,3 \times 10^{-26}$, V de Cramér = 0,67), enquanto a associação com o diagnóstico não foi significativa ($p = 0,41$).
   • Restringir a análise a uma única região (DLPFC) fez com que a estabilidade caísse abaixo do limiar de validação (ARI = 0,68), confirmando que a "estabilidade" era um artefato do eixo da região.

3. Natureza Sistemática do Artefato:
   Este artefato impulsionado pela região foi sistemático em todos os diagnósticos da coorte (Bipolar, Depressão Maior) e invariante ao número de clusters ($k$). A coorte combinada das quatro diagnósticos apresentou a estrutura regional mais forte, mas a estabilidade mais fraca.

4. Recuperação de Biologia Conhecida via Análise Não Circular:
   Embora a estrutura de agrupamento tenha perdido os sinais de doença, uma análise convencional de caso-controle não circular recuperou a biologia de interneurônios estabelecida:

   • TEA: Depleção significativa de marcadores de parvalbumina (PVALB) no córtex.
   • Esquizofrenia: Depleção significativa de marcadores de somatostatina (SST) e enzimas de síntese de GABA (GAD1/GAD2) em regiões corticais.
   • Dupla Dissociação: Uma interação de sinal oposto confirmou que a depleção de PVALB é mais pronunciada no TEA, enquanto a depleção de SST é mais pronunciada na esquizofrenia.

5. Falha de Benchmark:
   Um benchmark de calibração companheiro (um modelo de limiar adaptativo que alegava $R^2 = 0,889$) falhou em reproduzir os resultados. Os autores rastrearam isso a um defeito métrico: o Índice Rand Ajustado (ARI) retornou um valor espúrio de 1,0 em entradas degeneradas (datasets com um único cluster verdadeiro), inflando artificialmente a correlação.

Significância e Alegações
O artigo não alega descobrir novos subtipos moleculares; em vez disso, oferece um "estudo de caso aplicado" demonstrando como uma alegação de subtipo pode passar em rigorosos testes de estabilidade estatística, mas ser inteiramente impulsionada por um confundidor conhecido (região cerebral).

Os autores concluem que:

• Estabilidade $\neq$ Validade: Alta estabilidade de bootstrap certifica a reprodutibilidade, mas não garante que a partição reflita o sinal biológico pretendido (doença vs. anatomia).
• Teste de Confundidores é Obrigatório: Testar as partições contra a região cerebral (e outros confundidores) deve ser uma parte rotineira da bateria de validação para qualquer análise de subtipagem que combine regiões ou coortes.
• Disciplina sobre Descoberta: Em pequenas coortes combinadas de múltiplas regiões, a descoberta exploratória é propensa ao erro. Os autores defendem um checklist específico: relatar controles nulos, testar contra confundidores, garantir replicação com sujeitos disjuntos, evitar inferência circular e publicar partições que falham nos controles.

O estudo enfatiza que, embora a subtipagem molecular possa ter sucesso (citando Bowen et al., 2019, como um exemplo de sucesso controlado por região), os modos de falha específicos documentados aqui — impulsionados pela dominância anatômica e pelos padrões padrão (defaults) não examinados dos pipelines — são armadilhas críticas na genômica psiquiátrica atual.