$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Pour enregistrer simultanément à partir de jusqu’à 12 régions cérébrales SBN, nous avons construit un réseau multifibre personnalisable suivant le protocole ci-dessus. Après la récupération de l’opération, les souris ont été enregistrées lors d’interactions sociales (Figure 5). Pour toutes les expériences, les souris mâles du sujet CD1 étaient autorisées à interagir avec des souris Balb/c mâles dans leur cage domestique. Nous avons utilisé des souris Vgat-Cre pour exprimer jRCaMP1b40 dans des neurones GABAergiques de manière dépendante de Cre dans 11 régions cérébrales SBN, GCaMP6f39 dans des cellules non GABAergiques de manière excluantCre 45 dans les mêmes régions cérébrales, et GRAB-DA1m33 dans le noyau accumbens42 (Tableau 1). Une confirmation histologique du ciblage des fibres se trouve dans nos publicationsantérieures 42,43. Ici, nous avons cherché à tester la performance de différentes méthodes de correction des artefacts de mouvement. Plus précisément, nous avons comparé l’efficacité de l’utilisation de : (1) un filtre médian seul ; (2) Régression lasso avec régularisation L1 utilisant uniquement des signaux isosbestiques à fibres égales (approche à prédicteur unique)11 et (3) Régression lasso utilisant des signaux isosbestiques de toutes les fibres comme prédicteurs (approche multi-prédictrice)41,42 (Figure 4A). De manière cruciale, les méthodes basées sur Lasso corrigent le signal brut en utilisant le signal fluorescent indépendant de l’activité généré en réponse à la LED à fréquence isosestique (415 nm), tandis que le lissage du filtre médian ne prend pas en compte le signal indépendant de l’activité.
Nous avons évalué la correction de mouvement à l’aide de deux mesures différentes : le bruit absolu à haute fréquence et le jitter. Le bruit haute fréquence a été défini comme l’écart-type de la première dérivée du signal, qui reflète l’ampleur des fluctuations rapides du signal. Le jitter a été défini comme le coefficient de variation de la dérivée du signal, qui reflète la magnitude de la variabilité du signal d’une trame à l’autre par rapport à la moyenne de la dérivée du signal. Nous avons constaté que les deux approches basées sur le lasso (prédicteur unique et multi-prédictor) surpassaient nettement le filtre médian en termes de réduction du bruit et du jitter à haute fréquence, la correction multi-prédicteur surpassant légèrement la correction à un seul prédicteur (Figure 4C–D ; Bruit haute fréquence, n=414 pour 23 signaux répartis sur 18 sessions d’enregistrement : test de Kruskal-Wallis, H = 501,9420, p < 0,001 ; tests post-hoc de Mann-Whitney U avec correction de Bonferroni : régression lissée vs régression à site unique ***p < 0,001, régression lissée vs multi-sites***p < 0,001, Régression monosite vs multisite **p = 0,0073 ; Jitter : test de Kruskal-Wallis, H = 495,6685, p < 0,001 ; tests post-hoc de Mann-Whitney U avec correction de Bonferroni : régression lissée vs régression à site unique ***p < 0,001, régression lissée vs multi-site ***p < 0,001, Régression à site unique vs multisite ***p = 0,0006). Plus précisément, l’approche multi-prédicteur (utilisant toutes les isosestiques de fibre) a montré des bénéfices pour les signaux en canal rouge (RCaMP) : 5 des 12 signaux rouges ont montré un bruit résiduel significativement inférieur à la correction à un seul prédicteur (n = 18 enregistrements/région). Pour le jitter, le BNST (test t, t = 5,42, ***p < 0,001), le POA (test t, t = 2,68, *p = 0,011), l’AH (test t, t = 2,66, *p = 0,012) et le PAG (test U de Mann-Whitney, U = 279, ***p < 0,001) ont montré une supériorité multi-prédictrice. Pour le bruit à haute fréquence, le BNST (test U de Mann-Whitney, U = 285, ***p < 0,001), POA (test U de Mann-Whitney, U = 292, ***p < 0,001), AH (test t, t = 7,10, ***p < 0,001), PAG (test t, t = 13,72, ***p < 0,001) et LHb (test t, t = 5,47, ***p < 0,001) ont montré une supériorité multi-prédictif. Les signaux à canal vert (GCaMP) n’ont montré aucune différence significative entre les deux approches Lasso dans les 12 régions, que ce soit pour les métriques de gigue ou de bruit haute fréquence (toutes p>0,05, pour la liste complète des tests statistiques voir Tableau 3). Cela suggère que pour la photométrie bicolore, l’intégration d’informations isosbestiques transversales améliore la correction du mouvement pour les signaux rouges, probablement parce que les isosbestiques à canal vert fournissent des informations de mouvement plus nettes. Ces données démontrent que la correction des artefacts de mouvement utilisant des méthodes de régression basées sur l’isosbestique surpasse largement les approches de filtrage simples, la régression multi-canal offrant des avantages spécifiques pour les enregistrements bicolores.

Figure 1. Enregistrement à partir du Social Behavior Network utilisant la photométrie multisite à fibre bicolore. (A) Représentation schématique de la manière de construire un réseau multifibre. (B) Représentation schématique du flux de travail chirurgical. (C) Représentation schématique d’expériences d’imagerie calcique bicolore. (D) Exemple d’utilisation d’un signal isosbestique pour calculer un signal corrigé de mouvement. Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Figure 2. Images des interventions clés pour la fabrication et la chirurgie. (A) Cordon de brassage fabriqué prêt à être utilisé pour l’enregistrement. (B) Réseau multifibre placé dans un support stéréotaxique personnalisé avec la mèche de 0,25 mm x 5,5 mm placée au-dessus d’un trou traversant de fibre. (C) Réseau multifibre avec tous les trous traversants de la fibre dégagés et écrous hexagonals collés. (D) Fendre une fibre à l’aide d’un graveur de fibre. (E) Réseau multifibres avec toutes les fibres collées et placées sur une rondelle de polissage personnalisée. (F) Vue agrandie de la façon d’aligner et de visser le connecteur du cordon de brassage à la matrice multifibre. (G) Réseau multifibre entièrement construit avec des goupilles de baguette collées en place. (H) Souris fixée à un cadre stéréotaxique avec incision sur le cuir chevelu réalisée pour exposer le crâne. (I) Réseau multifibre fabriqué attaché à un manipulateur stéréotaxique pour manœuvrer le réseau multifibre pendant la chirurgie. (J) Vue agrandie de l’ensemble multifibre en train d’être descendu dans le cerveau. (K) Réseau multifibre entièrement abaissé et fixé au crâne avec du ciment dentaire. (L) Réseau multifibre sur un animal qui est solidement connecté au cordon de brassage avant l’enregistrement. (M) Vue dézoomée de l’animal en cage avec un cordon de brassage connecté avant et/ou pendant l’enregistrement. Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Figure 3. Séparation spectrale. Exemples de traces montrant une diaphonie spectrale minimale entre indicateurs calciques. Les signaux GCaMP6f (bleus) montrent une activation robuste lors de l’excitation à 470 nm mais une réponse négligeable à une excitation à 560 nm. Les signaux jRCaMP1b (rouge) présentent le motif inverse avec une forte activation à partir de 560 nm mais une réponse minimale à une illumination de 470 nm. (A) Enregistrer simultanément les signaux GCaMP6f et jRCaMP1b avec un éclairage uniquement de 470 nm. (B) Identique à (a), mais avec un éclairage uniquement à 560 nm. Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Figure 4. Comparaison des différentes méthodes pour le traitement du signal multifibre. (A) Schéma représentant la correction multi-site basée sur la régression du signal photométrique. Cette approche implique des étapes de traitement séquentielles par lesquelles tous les signaux isosbestiques disponibles à 415 nm sont mappés linéairement à chaque signal dépendant du calcium. Les composantes du mouvement prédite sont soustraites de ce signal pour produire la lecture finale de l’activité. Signaux ΔF/F à score z provenant de deux signaux présentant une forte présence d’artefacts de mouvement. (B) De haut en bas : signal isossestique provenant du même canal fibre/imagerie représentant des artefacts de mouvement, signal RCaMP (à gauche) ou GCaMP (à droite) uniquement lissé (filtré à médiane, vert), même signal à fibre/canal d’imagerie corrigé par régression (cyan), signal multi-fibres corrigé par régression (bleu). (C) Gauche : Score moyen de bruit à haute fréquence pour chaque signal enregistré (n = 21 sessions). Le bruit à haute fréquence est calculé comme l’écart-type de la première dérivée de chaque signal. Les scores de bruit sont affichés pour 1) un signal corrigé par la ligne de base mais non corrigé (gris), 2) un signal corrigé par la ligne de base et filtré par la médiane (rose), 3) un signal corrigé par régression de base et par fibre unique avec lissage Savitsky-Golay (cyan), et 4) un signal corrigé par régression de base et multi-fibres avec lissage Savitsky-Golay (bleu). Centre : Pourcentage moyen de bruit éliminé en utilisant uniquement le lissage, la régression à fibre unique ou multi-fibres, agrégé sur les 23 signaux. Le pourcentage de bruit éliminé est calculé comme la fraction de bruit de haute fréquence résultant d’une correction par rapport au niveau global de bruit dans le signal brut corrigé de base. Les données sont rapportées comme moyenne ± SEM (test de Kruskal-Wallis, H = 501,9420, p < 0,001 ; tests post-hoc Mann-Whitney U avec correction de Bonferroni : régression lissée vs monosite ***p < 0,001, Régression lissée vs Multi-site ***p < 0,001, Single-site vs Multi-site **p = 0,0073). À droite : Pourcentage de signaux où le bruit à haute fréquence est statistiquement différent après une correction basée sur la régression entre une seule fibre et une multi-fibre, séparée par type de signal. Pour chacun des 23 signaux, les niveaux de bruit après correction à site unique ont été comparés à ceux après correction multi-sites en utilisant soit des tests t indépendants (s’ils sont normalement distribués), soit des tests Mann-Whitney U (alternative non paramétrique). Les segments bleus indiquent des signaux où la correction multi-site a entraîné un bruit significativement moindre que la correction à site unique (p < 0,05). Les segments gris indiquent des signaux sans différence significative entre les méthodes de correction. Les segments cyan indiquent des signaux où une correction sur un seul site a entraîné un bruit significativement moindre. E : Vgat- GCaMP6f. DA : ATTRAPE-DA. I : Vgat+ jRCaMP1b. (D) Identique à (C) mais représentant la variabilité de trame à image (jitter), calculée comme la dérivée du coefficient de variation de chaque signal (test de Kruskal-Wallis, H = 495,6685, p < 0,001 ; tests post-hoc de Mann-Whitney U avec correction de Bonferroni : régression lissée vs monosite ***p < 0,001, régression lissée vs multi-site ***p < 0,001, Single-site vs Multi-site ***p = 0,0006). Pour la liste complète des tests statistiques et des résultats, voir le Tableau 3. Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Figure 5. Présentation multisite de photométrie à fibre bicolore. Des signaux représentatifs issus de plusieurs types d’indicateurs enregistrés simultanément lors du comportement social (test résident-intrus, époque d’entrée de l’intrus). Les traces montrent des GCaMP6f non normalisés et soustraits de la ligne de base mais non corrigés par les artefacts de mouvement (neurones Cre-négatifs), jRCaMP1b (neurones GABAergiques Vgat-Cre+) et GRAB-DA1m (capteur de dopamine dans le noyau accumbens) sur les sites du réseau de comportement social chez le même animal. Veuillez cliquer ici pour voir une version agrandie de cette figurine.
Tableau 1. Coordonnées virales et longueurs de fibres pour chaque Réseau Comportemental Social ciblé Région du Cerveau. Veuillez cliquer ici pour télécharger ce tableau.
Tableau 2. Comparaison des méthodes de correction des artefacts de mouvement pour la photométrie à fibres multi-sites. Principes sous-jacents, cas d’usage optimaux, avantages et limites des méthodes suivantes de correction d’artefacts de mouvement : régression linéaire simple, régression Lasso avec régularisation L1 (utilisant une approche multi-prédictrice), modélisation générative bayésienne, scaling basé sur la fréquence, et moindres carrés itératifs reponderés. Veuillez cliquer ici pour télécharger ce tableau.
Tableau 3. Détails et résultats des tests statistiques utilisés dans la Figure 4. Veuillez cliquer ici pour télécharger ce tableau.
Dossier supplémentaire 1. Fichier d’impression 3D pour le connecteur de patchcord.Veuillez cliquer ici pour télécharger ce fichier.
Fichier supplémentaire 2. 3D fichier d’impression pour connecteur de câble sans trous traversants.Veuillez cliquer ici pour télécharger ce fichier.
Fichier supplémentaire 3. 3D fichier d’impression pour l’accessoire de perçage stéréotaxique.Veuillez cliquer ici pour télécharger ce fichier.
Dossier supplémentaire 4. Plan personnalisé du puck de polissage.Veuillez cliquer ici pour télécharger ce fichier.
Fichier supplémentaire 5. 3D fichier d’impression pour la matrice de photométrie multisite ciblant 12 sites SBN.Veuillez cliquer ici pour télécharger ce fichier.
Fichier supplémentaire 6. 3D fichier d’impression pour une matrice de photométrie multisite sans trous traversants ciblant les sites.Veuillez cliquer ici pour télécharger ce fichier.
Dossier supplémentaire 7. Fichier de configuration DAQ pour le logiciel DAQ.Veuillez cliquer ici pour télécharger ce fichier.
Dossier supplémentaire 8. Fichier de codage visuel pour le système de photométrie.Veuillez cliquer ici pour télécharger ce fichier.
Dossier supplémentaire 9. Le code Python utilisé dans le protocole Step 4 et pour les comparaisons d’artefacts de mouvement corrigé. Veuillez cliquer ici pour télécharger ce fichier.