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 insuffisant entre les optodes et le cuir chevelu, une rétention insuffisante des canaux ou des estimations de réseau incohérentes qui limitent l'interprétabilité. Comme illustré dans la Figure 3 et la Figure 4, l'ensemble 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, écart-type : 8,99 ; 10 femmes) et 19 témoins sains (TS ; âge moyen : 33,5 ans, écart-type : 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 révélé 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 de réseau 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 rang somme.
Prétraitement et construction du réseau
Le schéma expérimental et le flux de travail de prétraitement sont présentés dans Figure 1 et 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 des séries temporelles régionales à partir des nœuds du réseau 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 réalisation 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 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 réseau du mode par défaut (DMN) spécifique à une condition pour un participant représentatif souffrant de trouble dépressif majeur (MDD), accompagné du modèle de référence du réseau DMN dérivé de sujets contrôles sains. 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épendantes 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 DMN 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, le classement selon les différentes métriques de distance varie selon les conditions, reflétant l'information complémentaire capturée par chaque métrique sur les représentations du réseau.
Sélection de fréquence fondée sur un consensus
Pour intégrer les informations issues de quatre métriques, les rangs ont été regroupés pour chaque participant (Figure 5C). Une séparation nette des scores de rang cumulés entre les conditions indique une différenciation réussie des réponses spécifiques aux fréquences au sein du réseau, tandis que des scores de rang très similaires 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 de rang cumulé le plus faible 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 l'accord entre les métriques quantifié par le coefficient W de Kendall.
Parmi 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 sélectionnée variait d’un participant à l’autre, comme résumé dans le 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 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) Conception expérimentale intra-sujets. Chaque participant a complété toutes les conditions lors d'une seule séance, incluant un é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 proche infrarouge fonctionnelle (fNIRS) ont été enregistrées en continu (~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 de 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 au cuir chevelu (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 pour les mouvements grâce à la méthode de 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), notamment le cortex préfrontal médian (mPFC), le précunéus (Prec), le lobule pariétal inférieur (IPL) et le lobe temporal médial (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 les 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) 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 à 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 à travers 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 des données au niveau des participants (≥ 70 % des canaux totaux conservés). Les participants situés en dessous du seuil sont mis en évidence (HC13, HC14, HC15, HC18). (B) Qualité des données par canal à travers les participants. La carte thermique montre le pourcentage de 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 consulter une version agrandie de cette figure.

Figure 5. Sélection représentative fondée sur la distance pour 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 du mode par défaut pour 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 comprennent 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 classements pour chaque condition. La condition présentant le classement total le plus bas a été choisie 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 la 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 du témoin sain est indiquée pour quatre métriques basé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 proposée fondée sur l'agrégation des classements selon 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 des gains de stimulation obtenus chez l'ensemble des participants pour chaque métrique, ainsi que les valeurs correspondantes de W de Kendall.