$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Aby jednocześnie rejestrować z maksymalnie 12 obszarów mózgu SBN, zbudowaliśmy konfigurowalną matrycę wielowłóknową zgodnie z powyższym protokołem. Po powrocie do zdrowia po operacji myszy były rejestrowane podczas interakcji społecznych (Rysunek 5). We wszystkich eksperymentach samce myszy z grupy CD1 mogły wchodzić w interakcje z samcami myszy Balb/c w ich domowej klatce. Użyliśmy myszy Vgat-Cre do ekspresji jRCaMP1b40 w neuronach GABAergicznych w sposób zależny od Cre w 11 obszarach mózgu SBN, GCaMP6f39 w komórkach nie-GABAergicznych w sposób wykluczający Cre45 w tych samych obszarach mózgu oraz GRAB-DA1m33 w jądrze accumbens42 (Tabela 1). Histologiczne potwierdzenie targetowania światłowodów można znaleźć w naszych wcześniejszych publikacjach42,43. Tutaj chcieliśmy przetestować wydajność różnych metod korekcji artefaktów ruchu. Konkretnie porównaliśmy skuteczność stosowania: (1) samego filtra mediany; (2) regresja lasowa z regularizacją L1, wykorzystująca wyłącznie sygnały izosbesticzne tego samego światła (podejście pojedynczego predyktora)11 oraz (3) regresja lassowa z wykorzystaniem sygnałów izobesticznych ze wszystkich światłowodów jako predyktorów (podejście multi-predyktorów)41,42 (Rysunek 4A). Co istotne, metody oparte na Lasso korygują sygnał surowy, wykorzystując niezależny od aktywności sygnał fluorescencjonalny generowany w odpowiedzi na diodę LED o częstotliwości izobestowej (415 nm), podczas gdy wygładzanie filtra medianowego nie uwzględnia sygnału niezależnego od aktywności.
Ocenialiśmy korekcję ruchu za pomocą dwóch różnych miar: absolutnego szumu wysokiej częstotliwości oraz drgań. Szum wysokich częstotliwości został zdefiniowany jako odchylenie standardowe pierwszej pochodnej sygnału, które odzwierciedla wielkość szybkich fluktuacji sygnału. Jitter definiowano jako współczynnik zmienności pochodnej sygnału, który odzwierciedla wielkość zmienności sygnału klatka do klatki względem średniej pochodnej sygnału. Stwierdziliśmy, że oba podejścia oparte na Lasso (pojedynczy i wielopredyktorowy) znacznie przewyższały medianowy filtr pod względem redukcji zarówno szumu wysokich częstotliwości, jak i jitterów, przy czym korekcja wielopredyktorowa nieznacznie przewyższała korekcję pojedynczej predyktorowej (Rysunek 4C–D; Szum wysokich częstotliwości, n=414 dla 23 sygnałów na 18 sesjach nagraniowych: test Kruskala-Wallisa, H = 501,9420, s. < 0,001; Post-hoc testy Mann-Whitney U z korektą Bonferroniego: Regresja wygładzona vs jednomiejscowa ***p < 0,001, wygładzona vs wielomiejscowa regresja***p < 0,001, jednomiejscowa vs wielomiejscowa **p = 0,0073; Jitter: test Kruskala-Wallisa, H = 495,6685, s. < 0,001; post-hoc testy Mann-Whitney U z korektą Bonferroniego: Wygładzony vs regresja pojedynczego miejsca ***p < 0,001, Wygładzony vs regresja wielomiejscowa***p < 0,001, Jednomiejscowy vs wielomiejscowy ***p = 0,0006). Konkretnie, podejście multi-predyktorowe (wykorzystujące wszystkie izozbestiki światłowodowe) wykazało korzyści dla sygnałów kanału czerwonego (RCaMP): 5 z 12 sygnałów czerwonych wykazywało znacząco niższy szum resztkowy w porównaniu z korekcją pojedynczego predyktora (n = 18 nagrań/region). Dla jittera BNST (t-test, t = 5,42, ***p < 0,001), POA (t-test, t = 2,68, *p = 0,011), AH (t-test, t = 2,66, *p = 0,012) oraz PAG (test Mann-Whitney U, U = 279, ***p < 0,001) wykazały przewagę w multi-predyktorach. Dla szumu wysokich częstotliwości BNST (test Mann-Whitney U, U = 285, ***p < 0,001), POA (test Mann-Whitney U, U = 292, ***p < 0,001), AH (test t, t = 7,10, ***p < 0,001), PAG (test t, t = 13,72, ***p < 0,001) oraz LHb (test t, t = 5,47, ***p < 0,001) wykazały przewagę w multi-predyktorach. Sygnały Green Channel (GCaMP) nie wykazały istotnych różnic między dwoma podejściami Lasso we wszystkich 12 regionach pod względem zarówno metryk szumu jitter, jak i wysokiej częstotliwości (wszystkie na stronie s>0.05, pełna lista testów statystycznych zobaczy Tabelę 3). Sugeruje to, że w fotometrii dwukolorowej włączenie informacji izobesticznych międzykanałowych poprawia korekcję ruchu sygnałów czerwonych, prawdopodobnie dlatego, że izozbestia z zielonych kanałów dostarcza czystszych informacji o ruchu. Dane te pokazują, że korekcja artefaktów ruchu z wykorzystaniem metod regresji opartej na izobestie znacznie przewyższa proste metody filtrowania, a regresja wielokanałowa zapewnia szczególne zalety dla nagrań dwukolorowych.

Rysunek 1. Rejestrowanie z Social Behavior Network za pomocą wielomiejscowej, dwukolorowej fotometrii światłowodowej. (A) Schematyczne przedstawienie sposobu budowy wieloświatłowodowej matrycy. (B) Schematyczne przedstawienie przepływu pracy chirurgicznej. (C) Schematyczne przedstawienie dwukolorowych eksperymentów obrazowania wapnia. (D) Przykład wykorzystania sygnału izobesticznego do obliczenia sygnału korygowanego ruchem. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.

Rysunek 2. Zdjęcia kluczowych procedur w produkcji i operacji. (A) Wykonany patchcord gotowy do nagrywania. (B) Wielowłóknowa matryca umieszczona w niestandardowym stereotaktycznym uchwytie z wiertłem 0,25 mm x 5,5 mm umieszczonym nad otworem przez światłowód. (C) Wielowłóknowa matryca z oczyszczonymi otworami przez włókna i przyklejonymi nakrętkami sześciokątnymi. (D) Cięcie włókna za pomocą skrybatora do włókien. (E) Wielowłóknowa matryca, z wszystkimi włóknami sklejonymi superglue i umieszczonymi na niestandardowym płycie do polerowania. (F) Powiększony widok sposobu ustawienia i przykręcenia złącza patchcorda do wieloświatłowodowej matrycy. (G) W pełni zbudowana sieć wielowłóknowa z przyklejonymi kołkami do kołków. (H) Mysz przymocowana do stereotaktycznej ramy z nacięciem na skórze głowy odsłaniając czaszkę. (I) Wykonana matryca wielowłóknowa przymocowana do manipulatora stereotaktycznego do manewrowania matrycą wielowłóknową podczas operacji. (J) Powiększony widok na wielowłóknową matrycę opuszczaną do mózgu. (K) Wielowłóknowa matryca całkowicie opuszczona i przymocowana do czaszki cementem dentystycznym. (L) Wielowłóknowa matryca na zwierzęciu, która jest bezpiecznie podłączona do patchcorda przed nagraniem. (M) Widok z oddalonym widokiem zwierzęcia w klatce z podłączonym kablem patch przed i/lub podczas nagrywania. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.

Rysunek 3. Spektralne rozdzielenie. Przykładowe ślady wykazujące minimalne przesłuchy spektralne między wskaźnikami wapnia. Sygnały GCaMP6f (niebieskie) wykazują solidną aktywację podczas wzbudzenia 470 nm, ale pomijalną odpowiedź na wzbudzenie 560 nm. sygnały jRCaMP1b (czerwone) wykazują odwrotny wzór z silną aktywacją od 560 nm, ale minimalną odpowiedzią na oświetlenie 470 nm. (A) Jednocześnie nagrywał sygnały GCaMP6f i jRCaMP1b z oświetleniem tylko 470 nm. (B) To samo co (a), ale z oświetleniem tylko 560 nm. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.

Rysunek 4. Porównanie różnych metod przetwarzania sygnałów wielowłóknowych. (A) Schemat przedstawiający wielomiejscową korektę regresji sygnału fotometrycznego. Podejście to obejmuje kolejne etapy przetwarzania, w których wszystkie dostępne sygnały izobesticzne 415nm są liniowo mapowane na każdy sygnał zależny od wapnia. Przewidywane składowe ruchu są odejmowane od tego sygnału, aby uzyskać końcowy odczyt aktywności. sygnały ΔF/F z dwóch sygnałów o silnej obecności artefaktów ruchu. (B) Od góry do dołu: sygnał izobesticzny z tego samego kanału światłowodowy/obrazującego przedstawiający artefakty ruchu, sygnał RCaMP (tylko wygładzony) lub GCaMP (prawy) (filtrowany medianą, zielony), ten sam sygnał regresji światłowodowej/kanału obrazowego (cyjan), sygnał regresji wieloświatłowodowej (niebieski). (C) Left: Średni wynik szumu wysokich częstotliwości dla każdego nagranego sygnału (n = 21 sesji). Szum wysokich częstotliwości obliczany jest jako odchylenie standardowe pierwszej pochodnej każdego sygnału. Wyniki szumów są pokazane dla 1) sygnału skorygowanego linią bazową, ale nieskorygowanego (szary), 2) sygnału skorygowanego i filtrowanego przez medianę (różowy), 3) sygnału regresji bazowej i jednowłóknowej z wygładzeniem Savitsky-Golay (cyjan) oraz 4) sygnału regresji bazowej i wielowłóknowej z wygładzeniem Savitsky-Golay (niebieski). Center: Średni procent szumu usuniętego wyłącznie przy użyciu wygładzania, regresji pojedynczego światłowodu lub wieloświatłowowowy, zagregowany przez wszystkie 23 sygnały. Procent usuniętego szumu oblicza się jako ułamek szumu wysokich częstotliwości powstałego w wyniku korekcji względem całkowitego poziomu szumu w surowym, bazowo skorygowanym sygnale. Dane podawane są jako średnia ± SEM (test Kruskala-Wallisa, H = 501,9420, s. < 0,001; post-hoc testy Mann-Whitney U z korektą Bonferroniego: Wygładzony vs regresja pojedynczego miejsca ***p < 0,001, Wygładzony vs regresja wielomiejscowa***p < 0,001, Jednoomiejscowy vs wielomiejscowy **p = 0,0073). Po prawej: Procent sygnałów, w których szum wysokiej częstotliwości jest statystycznie różny po korekcji opartej na regresji pojedynczego światłowodu i wieloświatłowodowej, rozdzielony według typu sygnału. Dla każdego z 23 sygnałów poziomy szumu po korekcie pojedynczego miejsca porównywano z poziomami szumu po korekcie wielomiejscowej za pomocą niezależnych testów t (jeśli są normalnie rozłożone) lub testów Mann-Whitney U (alternatywa nieparametryczna). Segmenty niebieskie wskazują sygnały, w których korekcja wielomiejscowa skutkowała znacznie mniejszym szumem niż korekcja pojedynczej lokalizacji (p < 0,05). Segmenty szare wskazują sygnały bez istotnych różnic między metodami korekcji. Segmenty cyjanowe wskazują sygnały, w których korekcja pojedynczego miejsca skutkowała znacznie mniejszym szumem. E: Vgat- GCaMP6f. DA: GRAB-DA. I: Vgat+ jRCaMP1b. (D) To samo co (C), ale przedstawiające zmienność klatki do klatki (jitter), obliczane jako pochodna współczynnika zmienności każdego sygnału (test Kruskala-Wallisa, H = 495,6685, s. < 0,001; testy post-hoc Mann-Whitney U z korektą Bonferroniego: Smoothed vs Single-site regression ***p < 0,001, Smoothed vs Multi-site regresse***p < 0,001, Single-site vs Multi-site ***p = 0,0006). Pełną listę testów i wyników statystycznych można znaleźć w Tabeli 3. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.

Rysunek 5. Wielomiejscowa prezentacja dwukolorowej fotometrii światłowodowej. Reprezentatywne sygnały z wielu typów wskaźników rejestrowane jednocześnie podczas zachowań społecznych (test mieszkanie-intruz, epoka wejścia intruza). Ślady wykazują nieznormalizowane i odejmowane według wartości bazowej, ale nie korygowane przez artefakty ruchu, GCaMP6f (neurony Cre-ujemne), jRCaMP1b (neurony GABAergiczne Vgat-Cre+) oraz GRAB-DA1m (czujnik dopaminy w jądrze przyściennym) w miejscach sieci zachowań społecznych u tego samego zwierzęcia. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.
Tabela 1. Współrzędne wirusa i długości włókien dla każdego docelowego obszaru mózgu Sieci Zachowań Społecznych. Prosimy kliknąć tutaj, aby pobrać tę tabelę.
Tabela 2. Porównanie metod korekcji artefaktów ruchu dla fotometrii światłowodowej na wielu miejscach. Podstawowe zasady, optymalne zastosowania, zalety i ograniczenia następujących metod korekcji artefaktów ruchu: prosta regresja liniowa, regresja Lasso z regularizacją L1 (z użyciem podejścia wielopredyktorowego), bayesowskie modelowanie generatywne, skalowanie oparte na częstotliwości oraz iteracyjne ponownie ważone metody najmniejszych kwadratów. Prosimy kliknąć tutaj, aby pobrać tę tabelę.
Tabela 3. Szczegóły i wyniki testów statystycznych użytych na Rysunku 4. Prosimy kliknąć tutaj, aby pobrać tę tabelę.
Plik uzupełniający 1. Plik druku 3D do złącza patchcorda.Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 2. 3D plik druku dla złącza patchcord bez otworów przejściowych.Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 3. 3D plik druku dla akcesorium do wiercenia stereotaktycznego.Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 4. Niestandardowy plan do polerowania krążków.Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 5. 3D plik druku dla wielomiejscowej matrycy fotometrycznej skierowanej do 12 stacji SBN.Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 6. 3D plik druku dla wielomiejscowej tablicy fotometrycznej bez przejścia przez otwory celujące w lokalizacje.Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 7. Plik konfiguracyjny DAQ dla oprogramowania DAQ.Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 8. Plik kodowania wizualnego dla systemu fotometrycznego.Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 9. Kod Pythona używany w kroku 4 protokołu oraz do porównań z poprawą artefaktów ruchu. Kliknij tutaj, aby pobrać ten plik.