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

Résumé technique : « Stable mais confondu : la région cérébrale, et non le diagnostic, pilote un sous-type moléculaire validé dans le cerveau psychiatrique postmortem »

Énoncé du problème
Le sous-typage moléculaire non supervisé des transcriptomes cérébraux est de plus en plus utilisé pour résoudre l'hétérogénéité des diagnostics psychiatriques (ex. : autisme, schizophrénie). Cependant, les algorithmes de partitionnement renverront toujours une partition, qu'elle reflète une biologie reproductible ou simplement des artefacts techniques ou anatomiques. L'article identifie trois dangers aigus dans les études psychiatriques postmortem :

1. Petites tailles d'échantillons (dizaines à quelques centaines) qui limitent la puissance statistique.
2. Chevauchement de sujets entre des cohortes nominalement indépendantes, menant à de fausses affirmations de réplication.
3. Dominance anatomique, où la variation transcriptomique par région cérébrale excède largement la variation par diagnostic. Par conséquent, le partitionnement de données multi-régions regroupées récupère souvent l'anatomie plutôt que des sous-types de maladies.

Les auteurs soutiennent que les pratiques de validation actuelles échouent souvent à distinguer un groupe « stable » (structure reproductible) d'un groupe « valide » (biologiquement pertinent pour la maladie), particulièrement lorsque les facteurs de confusion comme la région cérébrale ne sont pas explicitement testés.

Méthodologie
Les auteurs ont réanalysé trois cohortes postmortem publiques (384 échantillons au total) en utilisant une discipline de « validation d'abord » qui assemble des contrôles établis et ajoute des tests de confusion explicites. Les ensembles de données comprenaient :

• Trouble du spectre de l'autisme (TSA) : GSE28521 (microarray, cortex frontal/temporal/cervelet) et GSE64018 (RNA-seq, cortex temporal supérieur).
• Schizophrénie, trouble bipolaire, dépression majeure : GSE80655 (RNA-seq, cortex cingulaire antérieur, DLPFC, noyau accumbens).

Pipeline analytique :

1. Réduction des caractéristiques : L'expression génique a été résumée en scores de niveau de voie (pathway) via l'enrichissement à échantillon unique (ssGSEA ou moyenne des scores z par gène).
2. Partitionnement (Clustering) : Les échantillons ont été regroupés à l'aide d'une modélisation de mélange gaussien (GMM) avec sélection des composantes via le critère d'information bayésien (BIC).
3. Contrôles de validation : Chaque partition candidate a été soumise à quatre portes :
   • Null par permutation de labels : Tester si les groupes sont récupérables à partir de labels permutés.
   • Null par ensemble de gènes aléatoires : Tester si la partition est spécifique à des voies curatées ou si elle provient d'ensembles de gènes aléatoires.
   • Stabilité par bootstrap : Mesurer l'indice de Rand ajusté (ARI) entre les assignations originales et rééchantillonnées (seuil $\ge$ 0,80).
   • Indépendance de l'ascendance : Vérifier l'association avec les composantes principales de l'ascendance.
4. Test de confusion : Une association explicite par tableau de contingence ($\chi^2$, Cramér's V) a été réalisée entre les partitions et la région cérébrale, le chevauchement des donneurs et le diagnostic clinique.
5. Analyse non circulaire : Pour éviter le « double dipping », les signaux de la maladie ont été évalués via des contrastes cas-témoins basés sur le diagnostic au sein des régions, en utilisant un panel fixe de gènes marqueurs de types cellulaires canoniques (ex. : PVALB, SST, GAD1) plutôt que des caractéristiques dérivées du partitionnement lui-même.

Résultats clés

1. Échec du sous-type « GABAergique » de l'autisme :
   Un sous-groupe putatif de TSA caractérisé par une régulation à la baisse des voies GABAergiques a échoué à la validation de stabilité par bootstrap (ARI = 0,71, en dessous du seuil de 0,80). De plus, la réplication apparente dans une seconde cohorte (GSE64018) a été invalidée car les deux cohortes partageaient 10 sujets sur 24, ce qui signifie que l'accord était une corroboration inter-plateforme sur les mêmes cerveaux, et non une réplication indépendante.

2. La partition de la schizophrénie « validée » est un classificateur de région :
   Dans la cohorte de la schizophrenie (GSE80655), une partition $k=3$ a passé toutes les portes statistiques (Bootstrap ARI = 0,923 ; a passé les tests de gènes aléatoires et d'ascendance). Cependant, le test de confusion a révélé que cette partition était entièrement pilotée par la région cérébrale, et non par le diagnostic :

   • Les sous-types 0 et 2 étaient principalement issus du noyau accumbens (nAcc).
   • Le sous-type 1 était exclusivement cortical (AnCg et DLPFC).
   • L'association entre le sous-type et la région était écrasante ($p = 4,3 \times 10^{-26}$, Cramér's V = 0,67), tandis que l'association avec le diagnostic était non significative ($p = 0,41$).
   • La restriction de l'analyse à une seule région (DLPFC) a fait chuter la stabilité sous le seuil de validation (ARI = 0,68), confirmant que la « stabilité » était un artefact de l'axe de la région.

3. Nature systématique de l'artéfact :
   Cet artéfact piloté par la région était systématique à travers tous les diagnostics de la cohorte (Bipolaire, Dépression majeure) et invariant vis-à-vis du nombre de clusters ($k$). La cohorte regroupée des quatre diagnostics présentait la structure régionale la plus forte mais la plus faible stabilité.

4. Récupération de la biologie connue via l'analyse non circulaire :
   Alors que le cadre de partitionnement a manqué les signaux de la maladie, une analyse cas-témoins conventionnelle et non circulaire a récupéré la biologie des interneurones établie :

   • Autisme : Déplétion significative des marqueurs de la parvalbumine (PVALB) dans le cortex.
   • Schizophrénie : Déplétion significative de la somatostatine (SST) et des enzymes de synthèse du GABA (GAD1/GAD2) dans les régions corticales.
   • Double dissociation : Une interaction de signe opposé significative a confirmé que la déplétion de PVALB est plus prononcée dans l'autisme, tandis que la déplétion de SST est plus prononcée dans la schizophrénie.

5. Échec du benchmark :
   Un benchmark de calibration compagnon (un modèle à seuil adaptatif revendiquant un $R^2 = 0,889$) a échoué à se reproduire. Les auteurs ont attribué cela à un défaut de métrique : l'indice de Rand ajusté (ARI) renvoyait une valeur fallacieuse de 1,0 sur des entrées dégénérées (ensembles de données avec un seul vrai cluster), gonflant artificiellement la corrélation.

Signification et revendications
L'article ne prétend pas découvrir de nouveaux sous-types moléculaires ; il propose plutôt un « cas d'étude documenté » démontrant comment une revendication de sous-type peut passer des tests rigoureux de stabilité statistique tout en étant entièrement pilotée par un facteur de confusion connu (la région cérébrale).

Les auteurs concluent que :

• Stabilité $\neq$ Validité : Une haute stabilité par bootstrap certifie la reproductibilité mais ne garantit pas que la partition reflète le signal biologique visé (maladie vs anatomie).
• Le test de confusion est obligatoire : Tester les partitions contre la région cérébrale (et d'autres facteurs de confusion) doit être une composante routinière de la batterie de validation pour toute analyse de sous-typage regroupant des régions ou des cohortes.
• Discipline plutôt que découverte : Dans les petites cohortes multi-régions regroupées, la découverte exploratoire est sujette à l'échec. Les auteurs préconisent une liste de contrôle spécifique : rapporter les contrôles nuls, tester contre les facteurs de confusion, assurer une réplication avec des sujets disjoints, éviter l'inférence circulaire et publier les partitions qui échouent aux contrôles.

L'étude souligne que, bien que le sous-typage moléculaire puisse réussir (citant Bowen et al., 2019, comme un exemple de succès contrôlé par la région), les modes de défaillance spécifiques documentés ici — pilotés par la dominance anatomique et les défauts par défaut des pipelines — sont des pièges critiques dans la génomique psychiatrique actuelle.