$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Um bis zu 12 SBN-Gehirnregionen gleichzeitig aufzuzeichnen, haben wir ein anpassbares Multifiber-Array nach dem oben genannten Protokoll gebaut. Nach der Genesung von der Operation wurden die Mäuse während der sozialen Interaktion erfasst (Abbildung 5). Bei allen Experimenten durften männliche CD1-Probenmäuse mit männlichen Balb/c-Mäusen in ihrem Heimkäfig interagieren. Wir nutzten Vgat-Cre-Mäuse, um jRCaMP1b40 in GABAergen Neuronen cre-abhängig über 11 SBN-Gehirnregionen zu exprimieren, GCaMP6f39 in nicht-GABAergen Zellen in Cre-ausschließender Weise45 über dieselben Gehirnregionen und GRAB-DA1m33 im Nucleus accumbens42 (Tabelle 1). Histologische Bestätigung der Fasertargetierung findet sich in unseren früheren Publikationen42,43. Hier wollten wir die Leistung verschiedener Methoden zur Korrektur von Bewegungsartefakten testen. Konkret verglichen wir die Wirksamkeit der Verwendung von: (1) einem Medianfilter allein; (2) Lasso-Regression mit L1-Regularisierung, die ausschließlich gleichfaserige isosbestische Signale verwendet (Ein-Prädiktor-Ansatz)11 und (3) Lasso-Regression unter Verwendung isosbestischer Signale aller Fasern als Prädiktoren (Multi-Prädiktor-Ansatz)41,42 (Abbildung 4A). Entscheidend ist, dass die Lasso-basierte Methoden das Rohsignal mit dem aktivitätsunabhängigen fluoreszierenden Signal korrigieren, das als Reaktion auf die isosbestische Frequenz (415 nm) LED erzeugt wird, während die mediane Filterglättung das aktivitätsunabhängige Signal nicht berücksichtigt.
Wir bewerteten die Bewegungskorrektur anhand von zwei verschiedenen Maßen: absolutem Hochfrequenzrauschen und Jitter. Hochfrequentes Rauschen wurde als Standardabweichung der ersten Ableitung des Signals definiert, die die Größe schneller Signalfluktuationen widerspiegelt. Jitter wurde als der Variationskoeffizient der Ableitung des Signals definiert, der die Größe der Frame-zu-Frame-Signalvariabilität relativ zum Mittelwert der Ableitung des Signals widerspiegelt. Wir stellten fest, dass beide Lasso-basierte Ansätze (Einzel-Prädiktor und Multi-Prädiktor) den Medianfilter deutlich übertroffen in Bezug auf die Reduzierung von Hochfrequenzrauschen und Jitter, wobei die Multi-Prädiktor-Korrektur die Korrektion mit einem einzelnen Prädiktor geringfügig übertraf (Abbildung 4C–D; Hochfrequentes Rauschen, n=414 für 23 Signale über 18 Aufzeichnungssitzungen: Kruskal-Wallis-Test, H = 501,9420, S. < 0,001; Post-hoc Mann-Whitney U-Tests mit Bonferroni-Korrektur: Glättete vs. Einzelstandortregression ***p < 0,001, Glättete vs. Multi-Site-Regression ***p < 0,001, Einzelstandort vs. Multi-Site **p = 0,0073; Jitter: Kruskal-Wallis-Test, H = 495,6685, S. < 0,001; Post-hoc Mann-Whitney-U-Tests mit Bonferroni-Korrektur: Glättete vs. Einzelstandortregression ***p < 0,001, Glättete vs. Multi-Site-Regression ***p < 0,001, Einzelstandort vs. Multi-Site-Regression ***p = 0,0006). Konkret zeigte der Multi-Prädiktor-Ansatz (unter Verwendung aller Glasfaserisosbestika) Vorteile für Rotkanal-(RCaMP)-Signale: 5 von 12 roten Signalen wiesen signifikant weniger Restrauschen auf als die Korrektur eines einzelnen Prädiktors (n = 18 Aufnahmen/Region). Für Jitter zeigten 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) und PAG (Mann-Whitney U-Test, U = 279, ***p < 0,001) eine Multi-Prädiktor-Überlegenheit. Für hochfrequentes Rauschen zeigten BNST (Mann-Whitney U-Test, U = 285, ***p < 0,001), POA (Mann-Whitney U-Test, U = 292, ***p < 0,001), AH (t-Test, t = 7,10, ***p < 0,001), PAG (t-Test, t = 13,72, ***p < 0,001) und LHb (t-Test, t = 5,47, ***p < 0,001) Multi-Prädiktor-Überlegenheit. Green-Channel-(GCaMP)-Signale zeigten keine signifikanten Unterschiede zwischen den beiden Lasso-Ansätzen in allen 12 Regionen sowohl für Jitter- als auch für Hochfrequenzrauschmetriken (alle p>0,05, für die vollständige Liste statistischer Tests siehe Tabelle 3). Dies deutet darauf hin, dass bei der zweifarbigen Photometrie die Einbeziehung von kanalübergreifenden isosbestischen Informationen die Bewegungskorrektur für rote Signale verbessert, wahrscheinlich weil grüne Kanalisosbestik sauberere Bewegungsinformationen liefern. Diese Daten zeigen, dass die Bewegungsartefaktkorrektur mit isosbestischen Regressionsmethoden einfache Filteransätze deutlich übertrifft, wobei die Mehrkanalregression spezifische Vorteile für Dual-Color-Aufnahmen bietet.

Abbildung 1. Aufzeichnung aus dem Social Behavior Network mittels Multisite-Zweifarb-Faserphotometrie. (A) Schematische Darstellung, wie man ein Multifiber-Array aufbaut. (B) Schematische Darstellung des Operationsablaufs. (C) Schematische Darstellung von Zweifarb-Kalziumbildungsexperimenten. (D) Beispiel für die Verwendung eines isosbestischen Signals zur Berechnung eines bewegungskorrigierten Signals. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 2. Bilder von wichtigen Verfahren für Herstellung und Operation. (A) Gefertigtes Patchcord, das zur Aufnahme bereit ist. (B) Multifiber-Array, das in einen individuellen stereotaktischen Halter mit dem 0,25 mm x 5,5 mm Bohrer über einem Glasfaserdurchgang platziert wird. (C) Multifiber-Array mit allen Glasfaser-Durchlöchern und eingeklebten Sechskantmuffen. (D) Eine Faser mit einem Faserschreiber spalten. (E) Multifiber-Array, bei dem alle Fasern mit Sekundenkleber befestigt und auf einen individuellen Polierpuck gelegt werden. (F) Vergrößerte Ansicht, wie man den Patchcord-Anschluss am Multifiber-Array ausrichtet und verschraubt. (G) Voll gebautes Multifiber-Array mit Dübelstiften, die an Ort und Stelle geklebt sind. (H) Maus am stereotaktischen Rahmen befestigt, mit einem Kopfhautschnitt, der den Schädel freilegt. (I) Hergestelltes Multifiber-Array, das an einen stereotaktischen Manipulator angeschlossen ist, um das Multifiber-Array während der Operation zu steuern. (J) Vergrößerte Ansicht des Multifiber-Arrays, das ins Gehirn abgesenkt wird. (K) Multifiber-Array vollständig abgesenkt und mit Zahnkleber am Schädel befestigt. (L) Multifiber-Array an einem Tier, das vor der Aufnahme sicher mit dem Patchcord verbunden ist. (M) Herausgezoomte Ansicht des Tieres im Käfig mit Patchcord, das vor und/oder während der Aufnahme angeschlossen ist. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 3. Spektrale Trennung. Beispielspuren, die minimale spektrale Übersprechen zwischen Kalziumindikatoren zeigen. GCaMP6f (blaue) Signale zeigen eine robuste Aktivierung während der 470-nm-Anregung, aber eine vernachlässigbare Antwort auf die 560-nm-Anregung. jRCaMP1b (rote) Signale zeigen das umgekehrte Muster mit starker Aktivierung ab 560 nm, aber minimaler Antwort auf 470 nm Beleuchtung. (A) Gleichzeitig wurden GCaMP6f- und jRCaMP1b-Signale mit nur 470 nm Beleuchtung aufgenommen. (B) Dasselbe wie (a), aber nur mit 560 nm Beleuchtung. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 4. Vergleich verschiedener Methoden für Multifiber-Signalverarbeitung. (A) Schaltplan, der die multistandortbasierte regressionsbasierte Korrektur des photometrischen Signals darstellt. Dieser Ansatz umfasst sequentielle Verarbeitungsschritte, bei denen alle verfügbaren 415-nm-isosbestischen Signale linear auf jedes calciumabhängige Signal abgebildet werden. Vorhergesagte Bewegungskomponenten werden von diesem Signal abgezogen, um die endgültige Auslesung der Aktivität zu erzeugen. z-bewertete ΔF/F-Signale von zwei Signalen mit starker Präsenz von Bewegungsartefakten. (B) Von oben nach unten: isosbestisches Signal vom selben Glasfaser-/Bildkanal, das Bewegungsartefakte darstellt, nur geglättetes RCaMP (links) oder GCaMP (rechts) Signal (mediangefiltert, grün), dasselbe Glasfaser-/Bildkanal-regressionskorrigiertes Signal (Cyan), Multifaser-regressionskorrigiertes Signal (blau). (C) Links: Mittlerer Hochfrequenzrauschwert für jedes aufgezeichnete Signal (n = 21 Sitzungen). Hochfrequentes Rauschen wird als Standardabweichung der ersten Ableitung jedes Signals berechnet. Rauschwerte werden angezeigt für 1) basislinienkorrigiertes, aber unkorrigiertes Signal (grau), 2) basislinienkorrigiertes und mediangefiltertes Signal (pink), 3) basislinien- und einzelfaserregressionskorrigiertes Signal mit Savitsky-Golay-Glättung (cyan) und 4) Basislinien- und Multifaser-regressionskorrigiertes Signal mit Savitsky-Golay-Glättung (blau). Zentrum: Mittlerer Prozentsatz des entfernten Rauschens, der nur durch Glättung, Einzelfaser- oder Mehrfaserregression, aggregiert über alle 23 Signale. Der Prozentsatz des entfernten Rauschens wird als der Anteil des hochfrequenten Rauschens berechnet, der durch Korrektur im Verhältnis zum Gesamtrauschpegel des rohen, baseline-korrigierten Signals entsteht. Die Daten werden als Mittelwert ± SEM angegeben (Kruskal-Wallis-Test, H = 501,9420, p < 0,001; post-hoc Mann-Whitney-U-Tests mit Bonferroni-Korrektur: Glättete vs. Einzelstandortregression ***p < 0,001, Glättete vs. Multi-Site-Regression ***p < 0,001, Einzelstandort vs. Multisite-Regression **p = 0,0073). Rechts: Prozentsatz der Signale, bei denen hochfrequentes Rauschen statistisch unterschiedlich ist, nach einer Korrektur von Einzelfaser- versus Mehrfaserregression, getrennt nach Signaltyp. Für jedes der 23 Signale wurden die Rauschpegel nach der Single-Site-Korrektur mit den Rauschpegeln nach der Multi-Site-Korrektur verglichen, entweder mit unabhängigen t-Tests (wenn normalverteilt) oder Mann-Whitney-U-Tests (nichtparametrische Alternative). Blaue Segmente zeigen Signale an, bei denen Multi-Site-Korrektur zu deutlich geringerem Rauschen führte als Single-Site-Korrektur (S. < 0,05). Graue Segmente zeigen Signale an, die keinen signifikanten Unterschied zwischen den Korrekturmethoden haben. Cyan-Segmente zeigen Signale an, bei denen die Korrektur eines einzelnen Standorts zu deutlich geringerem Rauschen führte. E: Vgat- GCaMP6f. DA: GRAB-DA. I: Vgat+ jRCaMP1b. (D) Dasselbe wie (C), aber mit Bild-zu-Bild-Variabilität (Jitter), berechnet als Ableitung des Variationskoeffizienten jedes Signals (Kruskal-Wallis-Test, H = 495,6685, p < 0,001; Post-hoc Mann-Whitney-U-Tests mit Bonferroni-Korrektur: Glättet vs. Einzelstandort-Regression ***p < 0,001, Glättete vs. Multi-Site-Regression ***p < 0,001, Einzelstandort vs. Multi-Site ***p = 0,0006). Eine vollständige Liste statistischer Tests und Ergebnisse finden Sie in Tabelle 3. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 5. Multi-Site-Präsentation der zweifarbigen Faserphotometrie. Repräsentative Signale mehrerer Indikatortypen, die während des sozialen Verhaltens gleichzeitig aufgezeichnet werden (Bewohner-Eindringlingstest, Eindringlings-Epoche). Spuren zeigen unnormalisierte und von Baseline subtratraherte, aber nicht bewegungsartefakt-korrigierte GCaMP6f (Cre-negative Neuronen), jRCaMP1b (GABAerge Vgat-Cre+-Neuronen) und GRAB-DA1m (Dopaminsensor im Nucleus accumbens) über soziale Verhaltensnetzwerke desselben Tieres. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
Tabelle 1. Viruskoordinaten und Faserlängen für jede gezielte soziale Verhaltensnetzwerk-Hirnregion. Bitte klicken Sie hier, um diese Tabelle herunterzuladen.
Tabelle 2. Vergleich von Methoden zur Korrektur von Bewegungsartefakten für Multi-Site-Faserphotometrie. Grundlegende Prinzipien, optimale Anwendungsfälle, Vorteile und Einschränkungen folgender Bewegungsartefaktkorrekturmethoden: einfache lineare Regression, Lasso-Regression mit L1-Regularisierung (unter Verwendung des Multi-Prädiktor-Ansatzes), Bayessche generative Modellierung, frequenzbasierte Skalierung und iterativ umgewichtete kleinste Quadrate. Bitte klicken Sie hier, um diese Tabelle herunterzuladen.
Tabelle 3. Details und Ergebnisse der in Abbildung 4 verwendeten statistischen Tests. Bitte klicken Sie hier, um diese Tabelle herunterzuladen.
Ergänzende Akte 1. 3D-Druckdatei für den Patchcord-Anschluss.Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Datei 2. 3D Druckdatei für den Patchcord-Anschluss ohne Durchgangslöcher.Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Akte 3. 3D Druckdatei für stereotaktisches Bohrzubehör.Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Akte 4. Individueller Polier-Puck-Bauplan.Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Datei 5. 3D Druckdatei für Multisite-Photometrie-Arrays, die 12 SBN-Standorte anspricht.Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Datei 6. 3D Druckdatei für Multisite-Photometrie-Array ohne Durchbohrungs-Zielstellen.Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Akte 7. DAQ-Konfigurationsdatei für DAQ-Software.Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzungsakte 8. Visuelle Codierungsdatei für das Photometriesystem.Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzungsakte 9. Python-Code, der in Protokoll Schritt 4 verwendet wird und für Bewegungsartefakte-korrigierte Vergleiche verwendet wird. Bitte klicken Sie hier, um diese Datei herunterzuladen.