$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Le protocole décrit a été utilisé pour collecter des données TMS-EEG à partir du cortex préfrontal dorsolatéral gauche (DLPFC) d’un volontaire sain. L’étude a été approuvée par le Comité d’éthique de coordination de l’hôpital universitaire d’Helsinki et réalisée conformément à la Déclaration d’Helsinki. Le sujet a signé un formulaire de consentement éclairé. Pour plus de détails sur l’acquisition des données, reportez-vous à la Table des matériaux.
Le DLPFC gauche est un domaine d’intérêt commun, utilisé en pratique clinique pour le traitement répétitif par SMT du trouble dépressif majeur (TDM)65. De plus en plus de littérature souligne l’importance de la connectivité fonctionnelle du cortex cingulaire antérieur sous-génital (sgACC) avec le DLPFC dans la physiopathologie du TDM 3,66,67. En particulier, il existe des indications que le traitement par TMS administré sur les zones du DLPFC gauche montrant une forte anti-corrélation avec le sgACC gauche pourrait avoir une efficacité plus élevée 38,68,69,70. À la suite de ces résultats, on suppose que les données TMS-EEG enregistrées à partir des cibles fonctionnellement anticorrélées avec l’activité BOLD de sgACC seront les plus informatives pour le développement ultérieur de biomarqueurs.
Le masque d’IRM structurel pour le DLPFC a été défini comme une combinaison des régions 9a, 9p, 9-46d, 46, a9-46v et p9-46v du HCPMMP36, englobant les zones de Brodmann 9 et 46. Le masque sgACC était la région 25, équivalente à l’aire de Brodmann 25. Conformément au protocole décrit, des cartes de corrélation graine-voxel ont été dérivées à l’aide de données IRMf et sélectionnées uniquement des voxels appartenant au retour d’intérêt DLPFC. Les cartes ont ensuite été seuillées pour n’afficher que les voxels avec une corrélation négative ; La valeur absolue de la corrélation a été prise pour simplifier la visualisation. Les grappes résultantes ont été superposées sur l’IRM structurelle dans la neuronavigation pendant l’expérience TMS-EEG (Figure 2A).
Les cibles initiales de cartographie TMS-EEG ont été sélectionnées dans les régions qui présentaient la plus forte anticorrélation fonctionnelle avec le sgACC. Ces cibles ont ensuite été affinées de manière itérative par la procédure de cartographie TMS-EEG afin d’identifier les sites de stimulation les plus prometteurs en fonction de la qualité de la TEP et de leur proximité avec les plaques corticales présentant la plus forte anticorrélation avec le sgACC. L’ensemble finalisé des sites de stimulation potentiels a ensuite été analysé à l’aide d’une tractographie en temps réel. Comme il n’existe pas de voie directe connue entre le sgACC et le DLPFC, des voies indirectes de la substance blanche ont dû être prises en compte. Dans les cas où des voies directes ou indirectes spécifiques sont connues a priori, celles-ci doivent être prioritaires. Pour le sgACC, le réseau cérébral à grande échelle affecté par MDD3 a été pris en compte lors de la sélection des cibles basée sur la tractographie. La connectivité structurelle de la cible TMS a également été évaluée, les zones présentant une large connectivité à d’autres régions du cerveau, telles que le cortex préfrontal ventromédian, l’aire de Broca ou le lobe pariétal, étant considérées comme prioritaires.
À la suite du protocole de cartographie, la cible, l’orientation et l’intensité de la stimulation ont été choisies en fonction du fait que les TEP résultants étaient supérieurs à 6 μV et peu affectés par les artefacts. Chaque itération d’optimisation impliquait l’enregistrement de 20 impulsions avec un prétraitement minimal. La cible optimisée finale a ensuite été utilisée pour collecter 300 impulsions pour l’analyse des données.
Le pipeline de prétraitement minimal, conçu pour simuler l’utilisation du rt-TEP en temps réel, impliquait les étapes suivantes :
1. Division des données en époques
2. Correction de la ligne de base
3. Suppression de l’artefact TMS en remplaçant les données de -2 à 10 ms autour de l’impulsion TMS par des zéros
4. Suppression des canaux défectueux
5. Référençage à la référence moyenne
6. Un filtre passe-bas à 80 Hz et un filtre coupe-bande à 48-52 Hz ont été appliqués lorsque nécessaire pour traiter le bruit de ligne et le bruit à haute fréquence.
L’ensemble du pipeline de prétraitement a suivi les procédures décrites dans Mutanen et al., 202471. Pour plus de détails sur le prétraitement TMS-EEG, voir Hernandez-Pavon et al., 202272, Rogasch et al., 201773, Mutanen et al., 201874, et Mutanen et al., 201675. En résumé, le processus comprenait :
1. Division des données en époques
2. Correction de la ligne de base
3. Suppression des artefacts TMS entre -2 et 10 ms autour de l’impulsion TMS, avec interpolation cubique utilisant 5 ms de données avant et après l’intervalle supprimé
4. Suppression des canaux défectueux (4 au total : Fpz, F1, FT10, TP9, PO4)
5. Élimination des mauvais procès
6. Élimination des dérives avec un détrend robuste
7. ICA pour l’élimination des artefacts oculaires
8. Correction de la ligne de base
9. Application de l’algorithme SOUND74 pour supprimer le bruit extracrânien, suivie d’un réréférencement à la référence moyenne
10. Algorithme SSP-SIR75 pour supprimer les artefacts musculaires évoqués par TMS
11. Filtrage passe-bas à 80 Hz et filtre coupe-bande à 48-52 Hz
12. Couper les extrémités de la fenêtre temporelle pour supprimer les éventuels effets de bord
13. Suppression supplémentaire des essais défectueux pour s’assurer que les essais résiduels à fort bruit sont exclus (22 rejets sur 300 au total)
La procédure expérimentale a commencé par l’identification du seuil du moteur au repos. Le RMT du sujet était de 41 % MSO.
Au cours de la procédure expérimentale, ce sont principalement les électrodes situées dans la région d’intérêt (AF3, F1, F3, F5, FC3) (Figure 2B) afin d’évaluer la qualité des TEP résultants. D’autres électrodes ont fait l’objet d’une surveillance continue pour détecter une faible impédance, et l’impédance a été abaissée au besoin, conformément à la procédure décrite à l’étape 3.1.10. Les premières réponses de la TMS dans la fenêtre de temps de 10 à 60 ms, qui reflètent l’activation directe du cortex par le champ électrique induit, ont été utilisées pour évaluer la qualité des TEP produits.
La procédure de cartographie (voir Figure 3) a commencé avec l’un des groupes d’IRMf (Figure 3A) et une étude de la connectivité structurelle de la zone cible (Figure 3B), qui a principalement montré des lignes de courant vers la zone homologue dans l’hémisphère controlatéral et le pôle frontal.
Pour la stimulation, la bobine a été orientée à 45 degrés par rapport à la ligne médiane. L’intensité de la stimulation a été réglée à 49 % MSO, ce qui correspond à 96 V/m estimé au maximum du champ E au point chaud et à 120 % RMT. Dans la figure 3C, une réponse potentielle de type TEP est observée dans les électrodes F3, F1 et FC3. L’électrode AF3 présente un artefact de résonance de grande amplitude, tandis que F5 est affecté par un petit artefact musculaire reconnaissable par le pic de haute fréquence et de haute amplitude juste après l’impulsion TMS 76,77. Pour minimiser l’artefact de sonnerie pendant l’enregistrement, un capuchon de filet appuyant sur les électrodes pour réduire le mouvement ou une fine couche de mousse sous la bobine peut être utilisé (voir Hernandez-Pavon et al., 202328). Bien que ces artefacts puissent souvent être éliminés par filtrage, puisque seul AF3 a été contaminé, il a été exclu de la référence moyenne pour éviter un filtrage inutile, ce qui a donné les données présentées à la figure 3D.
En se concentrant sur l’électrode F3 (Figure 3F), le signal semble largement insensible à l’artefact musculaire observé chez F5. La forme d’onde caractéristique de grande amplitude récupère environ 15 ms, et une déviation de 5 microvolts entre 22 ms et 32 ms est probablement une véritable réponse corticale à la TMS. Le filtrage du signal (Fig. 3E, G) permet de confirmer que l’amplitude n’est pas affectée par le bruit dans les autres canaux.
Pour déterminer si l’artefact musculaire dans le canal F5 pouvait être réduit par la rotation de la bobine, tous les autres paramètres de stimulation ont été maintenus constants (Figure 4), et la bobine a d’abord été tournée vers une orientation postérieure-antérieure (Figure 4A), où le maximum du champ E a été estimé à 80 V/m, puis dans une direction latérale-médiale (Figure 4E), avec un champ E de 101 V/m maximum. Une forte augmentation de l’activation musculaire a été observée avec l’orientation de la bobine postérieure-antérieure (Figure 4B, C), qui a été signalée comme inconfortable par le sujet. L’orientation latérale-médiale a produit un signal similaire à celui montré sur la figure 3, mais avec une amplitude plus grande, entraînant une réponse TEP précoce de 12 μV à l’électrode F3 (Figure 4F,G), ce qui permet une réduction de l’intensité de la stimulation. L’absence d’artefacts musculaires excessifs, de décomposition ou de bourdonnement, ainsi que la taille de la réponse TEP précoce, font de cette combinaison de paramètres de stimulation un candidat prometteur pour la collecte de données.
La figure 5 montre l’effet de l’intensité de la stimulation sur les réponses TEP dans une cible à faible excitabilité. Le deuxième cluster d’IRMf (Figure 5A) et sa connectivité structurelle (Figure 5B) ont ensuite été étudiés. Comme pour la première cible, les connexions structurelles semblent se limiter aux lobes frontaux. Il convient de noter que, bien que la distance entre les deux groupes d’IRMf soit de 28 mm, l’intensité de stimulation utilisée pour les cibles précédentes, 49 % MSO (77 V/m, 120 % RMT), n’a produit aucune réponse TMS discernable, soulignant l’importance de la sélection de l’intensité de stimulation basée sur TEP. La figure 5B et C montre la réponse TMS à 55 % MSO (89 V/m, 134 % RMT). Il est important de noter que la grande déviation observée après 10 ms n’est pas une véritable réponse TMS, mais plutôt une continuation de la récupération de l’artefact musculaire. Par conséquent, seuls les signaux survenant plus de 25 ms après l’impulsion doivent être pris en compte, ce qui donne une amplitude de 4 μV entre 26 ms et 60 ms. L’intensité a encore été augmentée à 60 % MSO (97 V/m, 146 % RMT) dans le but d’obtenir des réponses distinctes plus tôt (figures 5F, G). Par rapport à l’intensité MSO de 55 %, une augmentation attendue de l’amplitude de l’artefact musculaire a été observée, mais il y avait peu ou pas d’augmentation de l’amplitude du TEP. Sur la base de la forme d’onde observée et de la nécessité d’une intensité de stimulation plus élevée, cette cible serait considérée comme présentant une excitabilité corticale plus faible par rapport aux cibles précédentes, ce qui la rendrait moins informative pour la collecte de données. Cependant, une cible aux propriétés similaires pourrait être préférable pour les protocoles de traitement.
Au total, 16 combinaisons distinctes de paramètres de stimulation ont été étudiées à l’aide d’une procédure de cartographie de 20 essais. La durée totale de l’intervention, à l’exclusion de la préparation de l’EEG et de l’EMG (qui a duré environ 30 min), a été de 2 h 35 min. La cible finale de stimulation peut être vue sur la figure 6A. Bien qu’elle soit très proche de la cible initiale (Figure 3A), cette cible présente une connectivité structurelle beaucoup plus répandue (Figure 6B), fournissant potentiellement plus d’informations sur la propagation du signal du point chaud vers d’autres régions corticales. L’intensité de la stimulation a été maintenue à 49 % MSO (102 V/m, 120 % RMT).
Le prétraitement minimal des 20 premières impulsions collectées à partir de la cible de stimulation finale est illustré à la figure 6C, D. Une résonance de grande amplitude a contaminé l’électrode F1, qui a ensuite été rejetée. Cependant, le bruit résiduel est resté présent dans d’autres électrodes même après ce rejet. Malgré la présence de l’artefact de résonance, l’absence de l’artefact musculaire dans la région d’intérêt a permis de considérer les latences après 16 ms comme un véritable signal neuronal. Après filtrage, l’amplitude de la composante précoce entre 17 ms et 35 ms était de 9 μV.
Les figures 6E, F montrent le résultat de l’ensemble du pipeline de prétraitement appliqué à l’ensemble de données de 300 impulsions. Le TEP résultant ressemble beaucoup à la forme d’onde dérivée de 20 impulsions (Figure 6D), démontrant la pertinence de la surveillance en temps réel des TEP. La réponse précoce entre 20 ms et 40 ms est de 6 μV, démontrant la réduction attendue de l’amplitude due au prétraitement.

Figure 2 : Neuronavigation et configuration des électrodes. (A) Modèle de tête basé sur l’IRM 3D du sujet superposé à une connectivité dérivée de l’IRMf. (B) Schéma de placement des électrodes EEG, avec les électrodes d’intérêt sous la bobine de stimulation marquées en orange. Notez que les emplacements des électrodes sont numérisés sur le modèle de tête 3D (étape 3.7.4). Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 3 : Un exemple de traitement minimal des données TMS-EEG. (A) L’écran de neuronavigation indique le site de stimulation et l’orientation de la bobine. La flèche rouge indique la direction la plus forte de l’impulsion biphasique, et la flèche bleue indique la direction la plus faible. L’intensité de la stimulation est exprimée en pourcentage de la puissance maximale du stimulateur (MSO) et du seuil moteur au repos (RMT). (B) Tractographie en temps réel au lieu de stimulation. (C) Données TMS-EEG brutes provenant d’électrodes situées sous la bobine de stimulation. (D) Données TMS-EEG après rejet d’un canal bruyant, les canaux rejetés étant remplacés par des zéros. (E) Données TMS-EEG après filtrage. (F) Vue agrandie de l’électrode F3 après rejet du mauvais canal. (G) Vue agrandie de l’électrode F3 après filtrage. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 4 : Effet de l’orientation de la stimulation sur l’artefact musculaire. (A) Cible de stimulation avec orientation de la bobine postérieure-antérieure et intensité correspondante. (B) Données brutes pour l’orientation postérieure-antérieure, montrant la contamination par des artefacts musculaires. (C) Données de l’électrode F3 avec un gros artefact musculaire. (D) Tractographie en temps réel correspondant à l’emplacement de la stimulation. (E) Cible de stimulation avec orientation latérale-médiale et intensité correspondante. (F) Données brutes pour l’orientation latérale-médiale, montrant des TEP propres. (G) TEP de l’électrode F3, présentant une réponse précoce importante et sans artefact. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 5 : Effet de l’intensité de la stimulation sur les réponses TEP dans une cible à faible excitabilité. (A) Emplacement et intensité de la cible de stimulation. (B) Données brutes enregistrées à 55 % de MSO. C) Données d’électrode F3, sans composants TEP précoces distincts. (D) Tractographie en temps réel correspondant à l’emplacement de la stimulation. (E) L’intensité de la stimulation a augmenté tout en maintenant la cible constante. (F) Données brutes enregistrées à 60 % de MSO. (G) Données des électrodes F3 montrant des artefacts musculaires dus à une intensité de stimulation plus élevée. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 6 : Effet des pipelines de prétraitement minimaux et complets sur les données TMS-EEG. (A) Cible et intensité de stimulation finale. (B) Tractographie en temps réel pour la localisation de la stimulation. (C) La moyenne des données brutes a été calculée sur 20 essais. (D) TEP moyenne de 20 essais à l’électrode F3. (E) Les TEP entièrement prétraités ont fait l’objet d’une moyenne de plus de 300 essais. (F) TEP moyen entièrement prétraité de 300 essais à l’électrode F3. Veuillez cliquer ici pour voir une version agrandie de cette figure.