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

Sintesi Tecnica: "Stabile ma confuso: la regione cerebrale, non la diagnosi, guida un sottotipo molecolare che supera la validazione nel cervello psichiatrico post-mortem"

Problema
La sottotipizzazione molecolare non supervisionata dei trascrittomi cerebrali è sempre più utilizzata per risolvere l'eterogeneità delle diagnosi psichiatriche (es. autismo, schizofrenia). Tuttavia, gli algoritmi di clustering restituiranno sempre una partizione, indipendentemente dal fatto che tale partizione rifletta una biologia riproducibile o semplicemente artefatti tecnici/anatomici. Il documento identifica tre pericoli acuti negli studi post-mortem psichiatrici:

1. Dimensioni del campione ridotte (decine o poche centinaia di campioni) che limitano la potenza statistica.
2. Sovrapposizione dei soggetti tra coorti nominalmente indipendenti, che porta a false affermazioni di replicazione.
3. Dominanza anatomica, in cui la variazione trascrizionale per regione cerebrale supera di gran lunga la variazione per diagnosi. Di conseguenza, il clustering di dati multi-regione aggregati spesso recupera l'anatomia piuttosto che i sottotipi di malattia.

Gli autori sostengono che le attuali pratiche di validazione spesso non riescono a distinguere tra un cluster "stabile" (struttura riproducibile) e un cluster "valido" (biologicamente rilevante per la malattia), specialmente quando i confondenti come la regione cerebrale non vengono esplicitamente testati.

Metodologia
Gli autori hanno rianalizzato tre coorti post-mortem pubbliche (384 campioni totali) utilizzando una disciplina "orientata alla validazione" che assembla controlli stabiliti e aggiunge espliciti test di confondimento. I dataset includevano:

• Disturbo dello Spettro Autistico (ASD): GSE28521 (microarray, corteccia frontale/temporale/cervelletto) e GSE64018 (RNA-seq, giro temporale superiore).
• Schizofrenia, Disturbo Bipolare, Depressione Maggiore: GSE80655 (RNA-seq, cingolato anteriore, DLPFC, nucleo accumbens).

Pipeline Analitica:

1. Riduzione delle Caratteristiche: L'espressione genica è stata riassunta in punteggi a livello di pathway tramite single-sample enrichment (ssGSEA o media dei valori z per gene).
2. Clustering: I campioni sono stati raggruppati utilizzando il Gaussian Mixture Modeling (GMM) con selezione dei componenti tramite il Criterio di Informazione Bayesiana (BIC).
3. Controlli di Validazione: Ogni partizione candidata è stata sottoposta a quattro gate:
   • Nulla per permutazione delle etichette: Test se i cluster sono recuperabili da etichette permutate.
   • Nulla per set di geni casuali: Test se la partizione è specifica per pathway curati o deriva da set di geni casuali.
   • Stabilità del bootstrap: Misurazione dell'Indice Rand Corretto (ARI) tra le assegnazioni originali e quelle ricampionate (soglia $\ge$ 0,80).
   • Indipendenza dall'ascendenza: Controllo dell'associazione con i componenti principali dell'ascendenza.
4. Test dei Confondenti: È stata eseguita un'associazione esplicita tramite tabella di contingenza ($\chi^2$, Cramér's V) tra le partizioni e la regione cerebrale, la sovrapposizione dei donatori e la diagnosi clinica.
5. Analisi Non Circolare: Per evitare il "double dipping", i segnali di malattia sono stati valutati tramite contrasti caso-controllo basati sulla diagnosi all'interno delle regioni, utilizzando un pannello fisso di geni marcatori canonici per il tipo cellulare (es. PVALB, SST, GAD1) anziché caratteristiche derivate dal clustering stesso.

Risultati Chiave

1. Fallimento del sottotipo "GABAergico" dell'Autismo:
   Un ipotetico sottogruppo ASD caratterizzato dalla downregulation dei pathway GABAergici ha fallito la validazione della stabilità del bootstrap (ARI = 0,71, al di sotto della soglia di 0,80). Inoltre, la presunta replicazione in una seconda coorte (GSE64018) è stata invalidata perché le due coorti condividevano 10 soggetti su 24, il che significa che l'accordo era una corroborazione cross-platform sugli stessi cervelli, non una replicazione indipendente.

2. La partizione della Schizofrenia che "passa la validazione" è un classificatore di regione:
   Nella coorte della schizofrenia (GSE80655), una partizione con $k=3$ ha superato tutti i gate statistici (Bootstrap ARI = 0,923; ha superato i controlli del set di geni casuali e dell'ascendenza). Tuttavia, il test dei confondenti ha rivelato che questa partizione era guidata interamente dalla regione cerebrale, non dalla diagnosi:

   • Il sottotipo 0 e 2 erano prevalentemente nucleo accumbens (nAcc).
   • Il sottotipo 1 era esclusivamente corticale (AnCg e DLPFC).
   • L'associazione tra sottotipo e regione era schiacciante ($p = 4,3 \times 10^{-26}$, Cramér's V = 0,67), mentre l'associazione con la diagnosi non era significativa ($p = 0,41$).
   • Limitare l'analisi a una singola regione (DLPFC) ha causato il calo della stabilità al di sotto della soglia di validazione (ARI = 0,68), confermando che la "stabilità" era un artefatto dell'asse regionale.

3. Natura Sistematica dell'Artefatto:
   Questo artefatto guidato dalla regione era sistematico in tutte le diagnosi della coorte (Bipolare, Depressione Maggiore) ed era invariante rispetto al numero di cluster ($k$). La coorte aggregata delle quattro diagnosi mostrava la struttura regionale più forte ma la stabilità più debole.

4. Recupero della Biologia Nota tramite Analisi Non Circolare:
   Sebbene il framework di clustering abbia mancato i segnali di malattia, un'analisi convenzionale caso-controllo non circolare ha recuperato la biologia degli interneuroni stabilita:

   • ASD: Significativa deplezione dei marcatori della parvalbumina (PVALB) nella corteccia.
   • Schizofrenia: Significativa deplezione della somatostatina (SST) e degli enzimi di sintesi GABA (GAD1/GAD2) nelle regioni corticali.
   • Doppia Dissociazione: Un'interazione con segno opposto e significativa ha confermato che la deplezione di PVALB è più pronunciata nell'ASD, mentre la deplezione di SST è più pronunciata nella schizofrenia.

5. Fallimento del Benchmark:
   Un benchmark di calibrazione compagno (un modello a soglia adattiva che dichiarava $R^2 = 0,889$) non è riuscito a riprodursi. Gli autori hanno tracciato questo fallimento in un difetto metrico: l'Indice Rand Corretto (ARI) restituiva un valore spurio di 1,0 su input degeneri (dataset con un singolo vero cluster), gonfiando artificialmente la correlazione.

Significato e Rivendicazioni
Il documento non pretende di scoprire nuovi sottotipi molecolari; piuttosto, offre un "caso di studio documentato" dimostrando come una rivendicazione di sottotipo possa superare rigorosi controlli di stabilità statistica pur essendo interamente guidata da un confondente noto (regione cerebrale).

Gli autori concludono che:

• Stabilità $\neq$ Validità: Un'alta stabilità del bootstrap certifica la riproducibilità ma non garantisce che la partizione rifletta il segnale biologico inteso (malattia vs anatomia).
• Il Test dei Confondenti è Obbligatorio: Testare le partizioni contro la regione cerebrale (e altri confondenti) deve essere una componente di routine della batteria di validazione per qualsiasi analisi di sottotipizzazione che aggreghi regioni o coorti.
• Disciplina sopra la Scoperta: In coorti multi-regione e aggregate con campioni piccoli, la scoperta esplorativa è soggetta a fallimenti. Gli autori raccomandano una checklist specifica: riportare i controlli nulli, testare contro i confondenti, garantire la replicazione con soggetti disgiunti, evitare l'inferenza circolare e pubblicare le partizioni che falliscono i controlli.

Lo studio sottolinea che, sebbene la sottotipizzazione molecolare possa avere successo (citando Bowen et al., 2019, come esempio di successo con controllo della regione), i modelli di fallimento specifici qui documentati — guidati dalla dominanza anatomica e dai default di pipeline non esaminati — sono criticità fondamentali nelle attuali genomiche psichiatriche.