Qualité des données et faisabilité de l'acquisition
L'application des critères prédéfinis de qualité des données a conduit à l'exclusion de canaux et de participants présentant une qualité de signal insuffisante. Les résultats satisfaisants se caractérisent par une forte rétention des canaux dans toutes les conditions, des enregistrements hémodynamiques stables et des schémas de connectivité interprétables permettant une estimation fiable des caractéristiques du réseau. Les résultats médiocres ou sous-optimaux incluent une faible qualité du signal, des artefacts de mouvement excessifs, un couplage inadéquat entre les optodes et le cuir chevelu, une rétention insuffisante des canaux ou des estimations de réseau incohérentes limitant l'interprétabilité. Comme illustré dans la Figure 3 et la Figure 4, le jeu de données final après évaluation de la qualité comprenait 19 participants souffrant de trouble dépressif majeur (TDM ; âge moyen : 31,21 ans, ÉT : 8,99 ; 10 femmes) et 19 témoins sains (TS ; âge moyen : 33,5 ans, ÉT : 13,51 ; 7 femmes). Sept participants (3 TDM et 4 TS) ont été exclus en raison d'une qualité de données insuffisante, principalement liée à des artefacts de mouvement ou à un couplage inadéquat entre les optodes et le cuir chevelu.
Les cartes de qualité par canal ont en outre mis en évidence des distributions spatiales hétérogènes de la qualité du signal entre les participants, avec des régions localisées de qualité réduite chez certains individus. Ces résultats soulignent l'importance d'appliquer des seuils prédéfinis basés sur l'indice de couplage au cuir chevelu et sur le pic du spectre de puissance avant toute analyse ultérieure. Ces procédures de contrôle de qualité garantissent que les mesures ultérieures des réseaux sont calculées à partir de données fiables et constituent un élément essentiel du protocole. Les résultats représentatifs présentés ici dépendent directement de la réussite du criblage de la qualité des données, du prétraitement normalisé, de l'extraction des séries temporelles basée sur les régions d'intérêt (ROI) et de l'application de la procédure de sélection prédéfinie par le rang et la somme.
Prétraitement et construction du réseau
Le schéma expérimental et le flux de travail de prétraitement sont illustrés dans la Figure 1 et la Figure 2. Après l'étape de contrôle de qualité (Figure 2), les signaux ont été correctement convertis en concentration d'hémoglobine et corrigés des artefacts de mouvement ainsi que du bruit physiologique systémique. L'achèvement réussi de ces étapes de prétraitement se manifeste par la conservation de canaux de haute qualité, des séries temporelles régionales stables et des estimations de connectivité interprétables, tandis que les résultats médiocres incluent une perte de signal excessive, des artefacts non résolus ou des mesures de réseau peu fiables. L'extraction de séries temporelles régionales à partir des nœuds du réseau en mode par défaut a permis de construire des matrices de connectivité fonctionnelle pour chaque condition. Ces résultats représentatifs dépendent directement de la bonne exécution du contrôle de qualité des données17, du prétraitement normalisé18, de la définition des régions d'intérêt (ROI) et de l'application de la procédure de classement prédéfinie. Les matrices de connectivité fonctionnelle du réseau en mode par défaut ainsi obtenues constituent la base des comparaisons ultérieures fondées sur la distance.
Résultats représentatifs au niveau du réseau
Figure 5A montre un exemple de profil de connectivité fonctionnelle du RMD spécifique à une condition pour un participant représentatif atteint de TED, aux côtés du modèle de référence du réseau RMD dérivé de sujets témoins en bonne santé. Ces matrices permettent une visualisation directe des schémas de connectivité interrégionale selon les conditions et aident à distinguer un résultat de protocole satisfaisant d’estimations de réseau moins interprétables ou de qualité médiocre. L’inspection visuelle révèle des variations dépendant de la condition dans l’organisation du réseau par rapport au modèle, illustrant comment la fréquence de stimulation peut moduler la structure du réseau RMD au cours d’une seule séance.
Les mesures de distance (Figure 5B) quantifient ces différences selon plusieurs métriques (connectivité, valeurs propres, vecteurs propres et contrôlabilité moyenne). Pour chaque métrique, les conditions peuvent être classées selon leur proximité par rapport au modèle de référence, des valeurs de distance plus faibles indiquant une similarité plus grande. Notamment, l'ordre obtenu à partir des différentes métriques de distance varie selon les conditions, reflétant ainsi l'information complémentaire captée par chaque métrique concernant les représentations du réseau.
Sélection de fréquence fondée sur un consensus
Afin d'intégrer les informations provenant de quatre métriques, les classements ont été regroupés pour chaque participant (Figure 5C). Une séparation nette des scores cumulés de classement entre les conditions indique une différenciation réussie des réponses spécifiques à la fréquence au sein du réseau, tandis que des scores de classement très proches entre les conditions suggèrent une sensibilité limitée du protocole ou une faible distinction entre les effets de stimulation. La condition présentant le score cumulé de classement le plus bas a été considérée comme la fréquence de stimulation candidate représentative pour ce participant. Tableau 1 résume ces sélections pour l'ensemble des participants, ainsi que les valeurs correspondantes de gain de stimulation et le degré d'accord entre les métriques quantifié par le coefficient W de Kendall.
Dans le groupe TDM, 12 des 19 participants ont présenté le rang cumulatif le plus bas pour une fréquence de stimulation active plutôt que pour la condition placebo. La fréquence de stimulation choisie variait selon les participants, comme résumé dans Tableau 1. L'accord entre les métriques, quantifié à l'aide du coefficient W de Kendall, différait également selon les individus, reflétant une variabilité dans la cohérence des classements. Des valeurs plus élevées du coefficient W de Kendall indiquent une concordance plus forte entre les métriques et une plus grande confiance dans le fait que plusieurs mesures de réseau soutiennent le même classement des conditions, tandis que des valeurs plus faibles indiquent un accord plus faible et une hétérogénéité plus marquée entre les métriques. Ainsi, la sélection finale doit être interprétée comme une intégration fondée sur les données de caractéristiques de réseau complémentaires, plutôt que comme une décision basée sur une seule métrique.
Ces résultats représentatifs et fondés sur la faisabilité illustrent les types de sorties au niveau du réseau générées par le protocole et démontrent que le flux de travail peut être mis en œuvre pour des comparaisons intra-sujets entre différentes conditions de stimulation. Les résultats visent à illustrer l'acquisition du signal, le prétraitement, l'extraction des caractéristiques et la cartographie paramétrique basée sur les rangs, plutôt qu'à constituer une preuve d'avantage thérapeutique, d'efficacité clinique ou de validation des performances, comme le montrent les Figures 1–5 et le Tableau 1.

Figure 1. Conception expérimentale et configuration de la sonde fNIRS. (A) Plan expérimental intra-sujets. Chaque participant a accompli toutes les conditions lors d'une seule séance, incluant l'état de repos de base (RS0), une stimulation factice, quatre fréquences de stimulation active (2, 10, 25 et 40 Hz) et un état de repos post-stimulation (RS1). L'ordre des stimulations actives était randomisé entre les participants. Chaque bloc durait 6 min et était séparé par des intervalles d'environ 2 min sans stimulation. Les données de spectroscopie fonctionnelle dans l'infrarouge proche (fNIRS) ont été enregistrées en continu (environ 54 min au total) et segmentées par condition pour l'analyse. (B) Disposition de la sonde fNIRS montrant 32 sources (rouges) et 30 détecteurs (bleus), formant 89 canaux à longue séparation et 8 canaux à courte séparation (entourés). Les canaux à courte séparation ont été utilisés pour éliminer le bruit physiologique superficiel. (C) Profil de sensibilité corticale de la disposition de la sonde par rapport au réseau neuronal par défaut (DMN), basé sur l'atlas d'Étiquetage Anatomique Automatisé (AAL), illustrant la sensibilité de mesure dépendante de la profondeur. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 2. Pipeline de prétraitement et d'extraction des caractéristiques pour l'analyse de réseau basée sur l'fNIRS.
Les données brutes d'intensité lumineuse ont été segmentées en blocs spécifiques aux conditions (RS0, faux traitement, fréquences de stimulation et RS1). La qualité des données a été évaluée à l'aide de l'indice de couplage cutané (SCI > 0,8) et du pic du spectre de puissance (PSP > 0,1), et seuls les canaux de haute qualité ont été conservés. Les signaux ont été convertis en densité optique, puis corrigés des mouvements grâce à la réparation par distribution des dérivées temporelles (TDDR) et à un filtrage par ondelettes. Les variations de concentration en hémoglobine ont été calculées, et une régression par canaux à courte distance a été appliquée afin de réduire le bruit physiologique systémique. Pour l'analyse du réseau, des séries temporelles ont été extraites de régions d'intérêt (ROIs) appartenant au réseau par défaut (DMN), incluant le cortex préfrontal médian (mPFC), le précunéus (Prec), le lobule pariétal inférieur (IPL) et le lobe temporal médian (MTL). La connectivité fonctionnelle a été calculée entre les ROIs. Les caractéristiques extraites comprenaient les valeurs de connectivité, les valeurs propres, les vecteurs propres et la contrôlabilité moyenne. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 3. Évaluation de la qualité des données chez des participants souffrant de trouble dépressif majeur (TDM).
(A) Nombre de canaux satisfaisant aux critères de qualité pendant la condition de base (RS0) pour chaque participant atteint de TDM. Les barres indiquent tous les canaux (gris), les canaux à longue séparation (bleu) et les canaux à courte séparation (orange) répondant aux critères de qualité (≥70 % du temps d’enregistrement avec SCI > 0,8 et PSP > 0,1). Les lignes en pointillés indiquent le nombre maximal de canaux à longue séparation disponibles (89), de canaux à courte séparation (8) et le seuil de qualité acceptable des données au niveau des participants (≥70 % du total des canaux conservés). Les participants situés en dessous du seuil sont mis en évidence (D13, D18, D19). (B) Qualité des données par canal chez les participants. La carte thermique montre le pourcentage du temps d’enregistrement satisfaisant aux critères de qualité pour chaque canal et chaque participant pendant RS0. Une teinte plus foncée indique une qualité de données plus faible. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 4. Évaluation de la qualité des données chez les participants témoins sains. (A) Nombre de canaux satisfaisant aux critères de qualité durant la condition de base (RS0) pour chaque participant témoin sain. Les barres indiquent tous les canaux (gris), les canaux à longue séparation (bleu) et les canaux à courte séparation (orange) satisfaisant aux critères de qualité (≥70 % du temps d’enregistrement avec SCI > 0,8 et PSP > 0,1). Les lignes en pointillés indiquent le nombre maximal de canaux disponibles et le seuil de qualité acceptable au niveau des participants (≥ 70 % des canaux totaux conservés). Les participants en dessous du seuil sont mis en évidence (HC13, HC14, HC15, HC18). (B) Qualité des données par canal, croisée entre participants. La carte thermique montre le pourcentage du temps d’enregistrement satisfaisant aux critères de qualité pour chaque canal et chaque participant durant RS0. Une teinte plus foncée indique une qualité de données plus faible. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 5. Sélection représentative, basée sur la distance, de la fréquence de stimulation chez un participant unique.
(A) Matrices de connectivité fonctionnelle montrant les interactions par paires entre quatre régions du réseau par défaut chez un participant représentatif atteint de trouble dépressif majeur (D11), dans différentes conditions (RS0, faux traitement et fréquences de stimulation). Les régions incluent le cortex préfrontal médian (mPFC), le lobule pariétal inférieur (IPL), le précunéus (Prec) et le lobe temporal médial (MTL). La matrice centrale représente le modèle de sujet témoin sain dérivé des données RS0 au niveau du groupe. (B) Mesures de distance entre chaque condition et le modèle de sujet témoin sain selon quatre métriques : distance de connectivité (D_conn), distance des valeurs propres (D_eig), distance des vecteurs propres (D_evec) et distance de la contrôlabilité moyenne (D_AC). Les barres représentent les valeurs de distance ; des valeurs plus faibles indiquent une plus grande similarité avec le modèle. Les classements sont indiqués pour chaque condition (classement 1 = distance la plus faible). (C) Classement consensuel entre les métriques, obtenu en additionnant les rangs pour chaque condition. La condition présentant le rang total le plus bas a été sélectionnée comme fréquence candidate de stimulation pour ce participant. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.
| ID du participant | Connectivité | Valeur propre | Vecteur propre | Contrôlabilité moyenne | Consensus | Kendall-W |
| D01 | 40HZ (G=1.00) | 40HZ (G=1.00) | 40HZ (G=1.00) | 40HZ (G=0.29) | 40HZ | 0.85 |
| D02 | Sham (G=0.18) | Sham (G=0.00) | Sham (G=0.11) | Sham (G=-0.09) | Sham | 0.89 |
| D03 | 10HZ (G=0.20) | Sham (G=0.05) | 2HZ (G=0.56) | Sham (G=0.63) | 40HZ | 0.34 |
| D07 | Sham (G=0.29) | 25HZ (G=0.06) | 2HZ (G=0.33) | 25HZ (G=0.45) | 25HZ | 0.36 |
| D10 | Sham (G=1.00) | Sham (G=1.00) | Sham (G=1.00) | Sham (G=1.00) | Sham | 0.59 |
| D11 | 10HZ (G=0.71) | 10HZ (G=0.66) | 2HZ (G=0.73) | 40HZ (G=1.00) | 10HZ | 0.34 |
| D15 | 40HZ (G=0.58) | 40HZ (G=0.53) | 40HZ (G=0.21) | 40HZ (G=-0.27) | 40HZ | 0.96 |
| D16 | 25HZ (G=1.00) | 10HZ (G=1.00) | 25HZ (G=0.90) | 25HZ (G=0.40) | 25HZ | 0.64 |
| D17 | 2HZ (G=0.51) | 2HZ (G=0.75) | 25HZ (G=0.73) | 10HZ (G=0.82) | 10HZ | 0.21 |
| D20 | 25HZ (G=-0.08) | 25HZ (G=0.04) | 25HZ (G=-0.27) | Sham (G=0.31) | 25HZ | 0.84 |
| D21 | Sham (G=0.37) | 2HZ (G=0.08) | Sham (G=0.27) | 40HZ (G=-0.28) | Sham | 0.41 |
| D22 | 2HZ (G=-0.06) | 10HZ (G=0.19) | 2HZ (G=0.75) | 2HZ (G=1.00) | 2HZ | 0.47 |
| D23 | Sham (G=0.53) | 25HZ (G=-0.14) | Sham (G=0.71) | 40HZ (G=0.23) | Sham | 0.31 |
| D24 | 25HZ (G=0.59) | 25HZ (G=0.18) | 25HZ (G=0.18) | Sham (G=0.34) | 25HZ | 0.5 |
| D25 | 25HZ (G=1.00) | Sham (G=-0.01) | 10HZ (G=1.00) | Sham (G=0.19) | 25HZ | 0.16 |
| D26 | 10HZ (G=0.52) | Sham (G=0.01) | 10HZ (G=1.00) | 10HZ (G=0.71) | 10HZ | 0.93 |
| D27 | 10HZ (G=-0.11) | 10HZ (G=0.20) | 2HZ (G=-0.20) | 40HZ (G=1.00) | 2HZ | 0.05 |
| D28 | 10HZ (G=0.51) | 10HZ (G=0.04) | 10HZ (G=1.00) | 10HZ (G=-0.02) | 10HZ | 0.81 |
| D29 | 10HZ (G=0.81) | 10HZ (G=0.54) | Sham (G=1.00) | 10HZ (G=1.00) | 10HZ | 0.85 |
| Gain | 0.50, 0.52 | 0.32, 0.18 | 0.58, 0.73 | 0.46, 0.40 | | 0.55, 0.50 |
| (Moyenne, Médiane) |
Tableau 1 : Sélection de fréquence spécifique à chaque participant selon les métriques de distance.
Pour chaque participant, la condition de stimulation associée à la distance minimale par rapport au modèle témoin sain est indiquée pour quatre métriques fondées sur les réseaux : distance de connectivité (D_conn), distance des valeurs propres (D_eig), distance des vecteurs propres (D_evec) et distance de la contrôlabilité moyenne (D_AC). Les valeurs entre parenthèses (G) indiquent le gain de stimulation. La colonne Consensus indique la fréquence de stimulation candidate déterminée par agrégation des classements sur les quatre métriques. Le coefficient W de Kendall quantifie l'accord entre les métriques pour chaque participant. La dernière ligne résume les valeurs moyennes et médianes du gain de stimulation obtenues sur l'ensemble des participants pour chaque métrique, ainsi que les valeurs correspondantes de W de Kendall.