$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Om gelijktijdig tot 12 SBN-hersengebieden op te nemen, hebben we een aanpasbare multifiberarray gebouwd volgens het bovenstaande protocol. Na herstel van de operatie werden de muizen geregistreerd tijdens sociale interactie (Figuur 5). Voor alle experimenten mochten mannelijke CD1-proefmuizen in hun thuiskooi interactie hebben met mannelijke Balb/c-muizen. We gebruikten Vgat-Cre-muizen om jRCaMP1b40 in GABAerge neuronen op een Cre-afhankelijke manier tot expressie te brengen over 11 SBN-hersengebieden, GCaMP6f39 in niet-GABAerge cellen op een Cre-exclusieve manier45 over dezelfde hersengebieden, en GRAB-DA1m33 in de nucleus accumbens42 (Tabel 1). Histologische bevestiging van vezeltargeting is te vinden in onze eerdere publicaties42,43. Hier probeerden we de prestaties van verschillende methoden voor motion artifact-correctie te testen. Specifiek vergeleken we de effectiviteit van het gebruik van: (1) alleen een mediaanfilter; (2) Lassoregressie met L1-regularisatie waarbij alleen isosbestische signalen van dezelfde vezel worden gebruikt (single-predictor benadering)11 en (3) Lassoregressie met isosbestische signalen van alle vezels als voorspellers (multi-predictor benadering)41,42 (Figuur 4A). Cruciaal is dat de Lasso-gebaseerde methoden het ruwe signaal corrigeren met behulp van het activiteitsonafhankelijke fluorescerende signaal dat wordt gegenereerd als reactie op de isosbestische frequentie (415 nm) LED, terwijl mediaanfiltersmoothing het activiteitsonafhankelijke signaal niet meeneemt.
We beoordeelden bewegingscorrectie met twee verschillende metingen: absolute hoogfrequente ruis en jitter. Hoogfrequente ruis werd gedefinieerd als de standaardafwijking van de eerste afgeleide van het signaal, die de omvang van snelle signaalfluctuaties weerspiegelt. Jitter werd gedefinieerd als de variatiecoëfficiënt van de afgeleide van het signaal, die de grootte van de frame-to-frame signaalvariabiliteit ten opzichte van het gemiddelde van de afgeleide van het signaal weerspiegelt. We ontdekten dat beide op lasso gebaseerde benaderingen (single-predictor en multi-predictor) het mediaanfilter aanzienlijk overtroffen in het verminderen van zowel hoogfrequente ruis als jitter, waarbij multi-predictor correctie marginaal beter presteerde dan single-predictor correctie (Figuur 4C–D; Hoogfrequente ruis, n=414 voor 23 signalen over 18 opnamesessies: Kruskal-Wallis-test, H = 501,9420, p < 0,001; post-hoc Mann-Whitney U-tests met Bonferroni-correctie: Gladde versus Single-site regressie ***p < 0,001, Smoothed vs Multi-site regressie ***p < 0,001, Single-site vs Multi-site **p = 0,0073; Jitter: Kruskal-Wallis test, H = 495,6685, p < 0,001; post-hoc Mann-Whitney U-tests met Bonferroni-correctie: Gladde vs Single-site regressie ***p < 0,001, Smoothed vs Multi-site regressie ***p < 0,001, Single-site vs Multi-site ***p = 0,0006). Specifiek toonde de multi-predictor benadering (met alle vezelisosbestics) voordelen voor rode kanaal (RCaMP) signalen: 5 van de 12 rode signalen vertoonden significant minder restruis vergeleken met correctie van enkele voorspeller (n = 18 opnames/regio). Voor jitter toonden 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) en PAG (Mann-Whitney U-test, U = 279, ***p < 0,001) superioriteit van meerdere voorspellers. Voor hoogfrequente ruis toonden 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) en LHb (t-test, t = 5,47, ***p < 0,001) superioriteit van multi-voorspellers. Green channel (GCaMP) signalen vertoonden geen significante verschillen tussen de twee Lasso-benaderingen in alle 12 regio's voor jitter- of hoogfrequente ruismetrieken (alle p>0,05, voor volledige lijst van statistische tests zie Tabel 3). Dit suggereert dat voor tweekleurige fotometrie het opnemen van cross-channel isosbestische informatie de bewegingscorrectie voor rode signalen verbetert, waarschijnlijk omdat groene kanaalisosbestica schonere bewegingsinformatie bieden. Deze gegevens tonen aan dat bewegingsartefactcorrectie met isosbestische regressiemethoden aanzienlijk beter presteert dan eenvoudige filtermethoden, waarbij multikanaalregressie specifieke voordelen biedt voor dual-color opnames.

Figuur 1. Opname vanaf het Social Behavior Network met behulp van multisite, tweekleurige vezelfotometrie. (A) Schematische weergave van hoe een multifiberarray gebouwd moet worden. (B) Schematische weergave van de operatieworkflow. (C) Schematische weergave van tweekleurige calciumbeeldvormingsexperimenten. (D) Voorbeeld van het gebruik van een isosbestisch signaal om een bewegingsgecorrigeerd signaal te berekenen. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 2. Afbeeldingen van belangrijke procedures voor fabricage en chirurgie. (A) Gefabriceerde patchkabel klaar om te worden gebruikt voor opname. (B) Multifiberarray geplaatst in een aangepaste stereotactische houder met de 0,25 mm x 5,5 mm boor boven een glasvezeldoorgang. (C) Multivezelarray met alle vezeldoorgangen vrijgemaakt en zeskantmoeren vastgelijmd. (D) Een vezel splijten met een vezelscribe. (E) Multifiberarray met alle vezels superlijm en geplaatst op een custom polijstpuck. (F) Vergrote weergave van hoe de patchcordconnector op de multifiberarray moet worden uitgelijnd en vastgeschroefd. (G) Volledig gebouwde multivezelarray met deuvelpennen die op hun plaats zijn gelijmd. (H) Muis bevestigd aan het stereotactische frame met een hoofdhuidincisie om de schedel bloot te leggen. (I) Gefabriceerde multifiberarray gekoppeld aan een stereotactische manipulator om de multifiberarray tijdens de operatie te manoeuvreren. (J) Vergroot beeld van de multifiberarray die in de hersenen wordt neergelaten. (K) Multivezelarray volledig neergelaten en met tandcement aan de schedel bevestigd. (L) Multifiberarray op een dier dat stevig is aangesloten op de patchkabel vóór opname. (M) Uitgezoomde weergave van dier in kooi met patchkabel aangesloten vóór en/of tijdens de opname. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 3. Spectrale scheiding. Voorbeeldsporen die minimale spectrale crosstalk tussen calciumindicatoren aantonen. GCaMP6f (blauwe) signalen tonen robuuste activatie tijdens 470 nm excitatie, maar een verwaarloosbare respons op 560 nm excitatie. jRCaMP1b (rood) signalen vertonen het omgekeerde patroon met sterke activatie vanaf 560 nm, maar minimale respons op 470 nm verlichting. (A) Gelijktijdig GCaMP6f- en jRCaMP1b-signalen opgenomen met alleen 470 nm verlichting. (B) Hetzelfde als (a), maar alleen met 560 nm verlichting. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 4. Vergelijking van verschillende methoden voor multivezel signaalverwerking. (A) Schema dat multi-site regressie-gebaseerde correctie van het fotometriesignaal weergeeft. Deze aanpak omvat sequentiële verwerkingsstappen waarbij alle beschikbare isosbestische signalen van 415nm lineair worden toegespand aan elk calciumafhankelijk signaal. Voorspelde bewegingscomponenten worden van dit signaal afgetrokken om de uiteindelijke uitlezing van activiteit te produceren. z-gescorede ΔF/F-signalen van twee signalen met sterke aanwezigheid van bewegingsartefacten. (B) Van boven naar beneden: isosbestisch signaal van hetzelfde vezel-/beeldvormingskanaal dat bewegingsartefacten weergeeft, alleen gladde RCaMP (links) of GCaMP (rechts) signaal (mediaangefilterd, groen), hetzelfde vezel-/beeldkanaal regressiegecorrigeerd signaal (cyaan), multivezel regressie-gecorrigeerd signaal (blauw). (C) Links: Gemiddelde hoogfrequente ruisscore voor elk opgenomen signaal (n = 21 sessies). Hoogfrequente ruis wordt berekend als de standaardafwijking van de eerste afgeleide van elk signaal. Ruisscores worden weergegeven voor 1) baseline-gecorrigeerd maar ongecorrigeerd signaal (grijs), 2) baseline-gecorrigeerd en mediaan-gefilterd signaal (roze), 3) baseline- en single-fiber regressie-gecorrigeerd signaal met Savitsky-Golay smoothing (cyaan), en 4) baseline- en multi-fiber regressie-gecorrigeerd signaal met Savitsky-Golay smoothing (blauw). Center: Gemiddeld percentage ruis verwijderd met alleen smoothing, single-fiber of multi-fiber regressie, geaggregeerd over alle 23 signalen. Het percentage verwijderde ruis wordt berekend als het aandeel hoogfrequente ruis dat voortkomt uit correctie ten opzichte van het totale ruisniveau in het ruwe, door de basislijn gecorrigeerde signaal. De gegevens worden gerapporteerd als gemiddelde ± SEM (Kruskal-Wallis test, H = 501,9420, p < 0,001; post-hoc Mann-Whitney U tests met Bonferroni-correctie: Gladde versus Single-site regressie ***p < 0,001, Smoothed vs Multi-site regressie ***p < 0,001, Single-site vs Multi-site **p = 0,0073). Rechts: Percentage signalen waarbij hoogfrequente ruis statistisch verschillend is volgens single-fiber versus multi-fiber regressie-gebaseerde correctie, gescheiden op signaaltype. Voor elk van de 23 signalen werden ruisniveaus na correctie op één locatie vergeleken met ruisniveaus na multi-site correctie, met behulp van onafhankelijke t-tests (indien normaal verdeeld) of Mann-Whitney U-tests (niet-parametrisch alternatief). Blauwe segmenten geven signalen aan waarbij multi-site correctie resulteerde in significant minder ruis dan single-site correctie (p < 0,05). Grijze segmenten geven signalen aan zonder significant verschil tussen correctiemethoden. Cyaansegmenten geven signalen aan waarbij correctie van één locatie tot aanzienlijk minder ruis leidde. E: Vgat- GCaMP6f. DA: GRAB-DA. I: Vgat+ jRCaMP1b. (D) Hetzelfde als (C) maar met frame-to-frame variabiliteit (jitter), berekend als de afgeleide van de variatiecoëfficiënt van elk signaal (Kruskal-Wallis test, H = 495,6685, p < 0,001; post-hoc Mann-Whitney U tests met Bonferroni-correctie: Gladde vs Single-site regressie ***p < 0,001, Smoothed vs Multi-site regressie ***p < 0,001, Single-site vs Multi-site ***p = 0,0006). Voor een volledige lijst van statistische tests en resultaten, zie Tabel 3. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 5. Multi-site presentatie van dual-color vezelfotometrie. Representatieve signalen van meerdere indicatortypes die gelijktijdig worden geregistreerd tijdens sociaal gedrag (resident-intruder test, indringer-intree-epoch). Sporen tonen niet-genormaliseerde en baseline-subtracted, maar niet bewegingartefact-gecorrigeerde GCaMP6f (Cre-negatieve neuronen), jRCaMP1b (GABAerge Vgat-Cre+ neuronen) en GRAB-DA1m (dopaminesensor in nucleus accumbens) over sociale gedragsnetwerksites in hetzelfde dier. Klik hier om een grotere versie van deze figuur te bekijken.
Tabel 1. Viruscoördinaten en vezellengtes voor elk doelwit sociaal gedragsnetwerk en hersengebied. Klik hier om deze tabel te downloaden.
Tabel 2. Vergelijking van bewegingsartefactcorrectiemethoden voor multi-site vezelfotometrie. Onderliggende principes, optimale gebruiksscenario's, voordelen en beperkingen van de volgende bewegingsartefactcorrectiemethoden: eenvoudige lineaire regressie, Lasso-regressie met L1-regularisatie (met multi-predictorbenadering), Bayesiaanse generatieve modellering, frequentiegebaseerde schaal en iteratief herwogen kleinste kwadraten. Klik hier om deze tabel te downloaden.
Tabel 3. Details en resultaten van de statistische tests die in Figuur 4 worden gebruikt. Klik hier om deze tabel te downloaden.
Aanvullend Dossier 1. 3D-printbestand voor patchcordconnector.Klik hier om dit bestand te downloaden.
Supplemental File 2. 3D printbestand voor patchcordconnector zonder doorgaande gaten.Klik hier om dit bestand te downloaden.
Supplemental File 3. 3D printbestand voor stereotactisch booraccessoire.Klik hier om dit bestand te downloaden.
Aanvullend Dossier 4. Op maat gemaakte polish-puckblauwdruk.Klik hier om dit bestand te downloaden.
Supplemental File 5. 3D printbestand voor multisite fotometrie-array gericht op 12 SBN-locaties.Klik hier om dit bestand te downloaden.
Supplemental File 6. 3D printbestand voor multisite fotometrie-array zonder doorgaande gaten die doelpunten targeten.Klik hier om dit bestand te downloaden.
Aanvullend dossier 7. DAQ-configuratiebestand voor DAQ-software.Klik hier om dit bestand te downloaden.
Aanvullend Dossier 8. Visueel codeerbestand voor fotometriesysteem.Klik hier om dit bestand te downloaden.
Aanvullend dossier 9. Python-code gebruikt in Protocol Step 4 en voor bewegingsartefacten corrigeerde vergelijkingen. Klik hier om dit bestand te downloaden.