$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Per registrare simultaneamente da un massimo di 12 regioni cerebrali SBN, abbiamo costruito un array multifibra personalizzabile seguendo il protocollo sopra descritto. Dopo il recupero dall'intervento, i topi sono stati registrati durante l'interazione sociale (Figura 5). Per tutti gli esperimenti, ai maschi di soggetto CD1 è stato permesso di interagire con maschi Balb/c nella loro gabbia domestica. Abbiamo utilizzato topi Vgat-Cre per esprimere jRCaMP1b40 nei neuroni GABAergici in modo dipendente da Cre in 11 regioni cerebrali SBN, GCaMP6f39 in cellule non GABAergiche in modoesclusivo Cre 45 nelle stesse regioni cerebrali, e GRAB-DA1m33 nel Nucleus accumbens42 (Tabella 1). La conferma istologica del targeting delle fibre si trova nelle nostre pubblicazioniprecedenti 42,43. Qui, abbiamo cercato di testare le prestazioni di diversi metodi per la correzione degli artefatti di movimento. In particolare, abbiamo confrontato l'efficacia dell'uso: (1) di un filtro mediano da solo; (2) Regressione del lazo con regolarizzazione L1 usando solo segnali isosbestici della stessa fibra (approccio a singolo predittore)11 e (3) Regressione del lazo usando segnali isosbestici di tutte le fibre come predittori (approccio multi-predittore)41,42 (Figura 4A). In modo fondamentale, i metodi basati su Lasso correggono il segnale grezzo usando il segnale fluorescente indipendente dall'attività generato in risposta al LED a frequenza isosbest (415 nm), mentre il levigamento del filtro mediano non considera il segnale indipendente dall'attività.
Abbiamo valutato la correzione del movimento utilizzando due diverse misure: rumore ad alta frequenza assoluta e jitter. Il rumore ad alta frequenza è stato definito come la deviazione standard della prima derivata del segnale, che riflette l'entità delle fluttuazioni rapide del segnale. Il jitter è stato definito come il coefficiente di variazione della derivata del segnale, che riflette la grandezza della variabilità del segnale da frame a frame rispetto alla media della derivata del segnale. Abbiamo riscontrato che entrambi gli approcci basati su Lazo (singolo predittore e multi-predittore) hanno notevolmente superato il filtro mediano in termini di riduzione sia del rumore ad alta frequenza che del jitter, con la correzione multi-predittore che ha superato leggermente la correzione a singolo predittore (Figura 4C–D; Rumore ad alta frequenza, n=414 per 23 segnali in 18 sessioni di registrazione: test di Kruskal-Wallis, H = 501,9420, p < 0,001; Test post-hoc Mann-Whitney U con correzione di Bonferroni: Regressione smussata vs a sito singolo ***p < 0,001, regressione smussata vs multisito ***p < 0,001, A sito singolo vs multisito **p = 0,0073; Jitter: test di Kruskal-Wallis, H = 495,6685, p < 0,001; Test post-hoc di Mann-Whitney U con correzione di Bonferroni: Regressione smussata vs a sito singolo ***p < 0,001, Regressione smussata vs multisito ***p < 0,001, A sito singolo vs multisito ***p = 0,0006). In particolare, l'approccio multi-predittore (utilizzando tutte le isosbest in fibra) ha mostrato benefici per i segnali a canale rosso (RCaMP): 5 segnali rossi su 12 hanno mostrato un rumore residuo significativamente inferiore rispetto alla correzione con un singolo predittore (n = 18 registrazioni/regione). Per il jitter, BNST (test t, t = 5,42, ***p < 0,001), POA (test t, t = 2,68, *p = 0,011), AH (test t, t = 2,66, *p = 0,012) e PAG (test U di Mann-Whitney, U = 279, ***p < 0,001) hanno mostrato superiorità multi-predittori. Per il rumore ad alta frequenza, il BNST (test U di Mann-Whitney, U = 285, ***p < 0,001), POA (test U di 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) e LHb (test t, t = 5,47, ***p < 0,001) hanno mostrato superiorità multi-predittori. I segnali del canale verde (GCaMP) non hanno mostrato differenze significative tra i due approcci Lasso in tutte e 12 le regioni, sia per le metriche di jitter che per il rumore ad alta frequenza (tutte p>0,05, per l'elenco completo dei test statistici vedi Tabella 3). Questo suggerisce che per la fotometria a doppio colore, incorporare informazioni isosbest trasversali migliora la correzione del movimento per i segnali rossi, probabilmente perché le isosbest a canale verde forniscono informazioni di movimento più pulite. Questi dati dimostrano che la correzione degli artefatti di movimento utilizzando metodi di regressione basati su isosbest supera sostanzialmente gli approcci di filtraggio semplici, con la regressione multicanale che offre vantaggi specifici per le registrazioni a doppio colore.

Figura 1. Registrazione dalla Social Behavior Network utilizzando fotometria a fibra multisito e bicolore. (A) Rappresentazione schematica di come costruire un array multifibra. (B) Rappresentazione schematica del flusso di lavoro chirurgico. (C) Rappresentazione schematica di esperimenti di imaging del calcio a due colori. (D) Esempio di utilizzo di un segnale isosbest per calcolare un segnale corretto di movimento. Clicca qui per visualizzare una versione più grande di questa figura.

Figura 2. Immagini delle procedure chiave per la fabbricazione e la chirurgia. (A) Patch cord fabbricato pronto per la registrazione. (B) Array multifibra inserito in un supporto stereotattico personalizzato con la punta da 0,25 mm x 5,5 mm posizionata sopra un foro passante per fibra. (C) Array multifibra con tutti i fori passanti della fibra liberati e i dadi esagonali incollati. (D) Tagliare una fibra usando una scrittura a fibra. (E) Array multifibra con tutte le fibre incollate con supercolla e posizionate su un disco di lucidatura personalizzato. (F) Vista ingrandita di come allineare e avvitare il connettore patch al multifiber. (G) Array multifibra completamente costruito con perni a perno incollati in posizione. (H) Topo fissato al telaio stereotassico con incisione al cuoio capelluto per esporre il cranio. (I) Array multifibra fabbricato collegato a un manipolatore stereotattico per manovrare l'array multifibra durante l'intervento. (J) Vista ingrandita dell'array multifibra che viene abbassata nel cervello. (K) Array multifibra completamente abbassata e fissata al cranio con cemento dentale. (L) Matrice multifibra su un animale che è collegata saldamente al cavo patch prima della registrazione. (M) Vista allontanata dell'animale in gabbia con un patch collegato prima e/o durante la registrazione. Clicca qui per visualizzare una versione più grande di questa figura.

Figura 3. Separazione spettrale. Esempi di tracce che dimostrano una minima interferenza spettrale tra indicatori di calcio. I segnali GCaMP6f (blu) mostrano un'attivazione robusta durante l'eccitazione a 470 nm ma una risposta trascurabile all'eccitazione a 560 nm. I segnali jRCaMP1b (rossi) mostrano il pattern inverso con forte attivazione da 560 nm ma risposta minima a un'illuminazione a 470 nm. (A) Registrato simultaneamente i segnali GCaMP6f e jRCaMP1b con illuminazione solo a 470 nm. (B) Uguale a (a), ma con illuminazione solo a 560 nm. Clicca qui per visualizzare una versione più grande di questa figura.

Figura 4. Confronto tra diversi metodi per l'elaborazione del segnale multifibra. (A) Schema che raffigura la correzione basata su regressione multi-sito del segnale di fotometria. Questo approccio prevede passaggi di elaborazione sequenziali in cui tutti i segnali isosbestici disponibili a 415nm vengono mappati linearmente a ciascun segnale dipendente dal calcio. Le componenti del moto previsto vengono sottratte da questo segnale per produrre la lettura finale dell'attività. segnali ΔF/F con punteggio z provenienti da due segnali con forte presenza di artefatti di movimento. (B) Dall'alto verso il basso: segnale isosbest dallo stesso canale fibra/imaging che rappresenta artefatti di movimento, segnale RCaMP solo smussato, o GCaMP (destra) (filtrato mediano, verde), stesso segnale corretto/regressione a fibra/canale di imaging (ciano), segnale corretto a regressione multifibra (blu). (C) Sinistra: punteggio medio di rumore ad alta frequenza per ogni segnale registrato (n = 21 sessioni). Il rumore ad alta frequenza viene calcolato come deviazione standard della prima derivata di ciascun segnale. I punteggi di rumore sono mostrati per 1) segnale corretto alla base ma non corretto (grigio), 2) segnale corretto alla linea di base e filtrato alla mediana (rosa), 3) segnale corretto per regressione di base e a fibra singola con levigatura Savitsky-Golay (ciano), e 4) segnale corretto per regressione a base e multifibra con levigatura Savitsky-Golay (blu). Centro: percentuale media di rumore rimossa usando solo smoothing, regressione single-fiber o multi-fibra, aggregata su tutti e 23 i segnali. La percentuale di rumore rimosso viene calcolata come la frazione di rumore ad alta frequenza risultante dalla correzione rispetto al livello complessivo di rumore nel segnale grezzo corretto di base. I dati sono riportati come media ± SEM (test di Kruskal-Wallis, H = 501,9420, p < 0,001; test post-hoc di Mann-Whitney U con correzione di Bonferroni: Regressione smussata vs a sito singolo ***p < 0,001, regressione smussata vs multisito ***p < 0,001, Regressione a sito singolo vs multisito **p = 0,0073). Destra: Percentuale di segnali in cui il rumore ad alta frequenza è statisticamente diverso dopo una correzione basata su regressione a fibra singola rispetto a multi-fibra, separata per tipo di segnale. Per ciascuno dei 23 segnali, i livelli di rumore dopo la correzione a sito singolo sono stati confrontati con quelli dopo la correzione multi-sito utilizzando test t indipendenti (se normalmente distribuiti) o test Mann-Whitney U (alternativa non parametrica). I segmenti blu indicano segnali in cui la correzione multi-sito ha prodotto rumore significativamente inferiore rispetto alla correzione a sito singolo (p < 0,05). I segmenti grigi indicano segnali senza differenze significative tra i metodi di correzione. I segmenti cian indicano segnali in cui la correzione a sito singolo ha prodotto un rumore significativamente inferiore. E: Vgat- GCaMP6f. DA: AFFA-DA. I: Vgat+ jRCaMP1b. (D) Uguale a (C) ma rappresenta la variabilità frame a frame (jitter), calcolata come derivata del coefficiente di variazione di ciascun segnale (test di Kruskal-Wallis, H = 495.6685, p < 0.001; test post-hoc Mann-Whitney U con correzione di Bonferroni: Regressione smussata vs a sito singolo ***p < 0.001, regressione smussata vs multisito ***p < 0.001, Regressione a sito singolo vs multisito ***p = 0.0006). Per l'elenco completo dei test statistici e dei risultati si veda la Tabella 3. Clicca qui per visualizzare una versione più grande di questa figura.

Figura 5. Esposizione multi-sito di fotometria a fibra bicolore. Segnali rappresentativi da più tipi di indicatori registrati simultaneamente durante il comportamento sociale (test del residente-intruso, epoca di ingresso dell'intruso). Le tracce mostrano GCaMP6f (neuroni Cre-negativi), jRCaMP1b (neuroni Vgat-Cre+ GABAergici e sottratti alla base, non normalizzati e sottratti dalla linea di base) ma non corretti da artefatti di movimento, tra i siti della rete comportamentale sociale nello stesso animale. Clicca qui per visualizzare una versione più grande di questa figura.
Tabella 1. Coordinate del virus e lunghezze delle fibre per ogni regione cerebrale della rete comportamentale sociale target. Clicca qui per scaricare questa tabella.
Tabella 2. Confronto dei metodi di correzione degli artefatti di movimento per la fotometria a fibre multi-sito. Principi di base, casi d'uso ottimali, vantaggi e limitazioni dei seguenti metodi di correzione degli artefatti del movimento: regressione lineare semplice, regressione Lasso con regolarizzazione L1 (usando approccio multi-predittore), modellazione generativa bayesiana, scalatura basata sulla frequenza e minimi quadrati riponderati iterativamente. Clicca qui per scaricare questa tabella.
Tabella 3. Dettagli e risultati dei test statistici utilizzati nella Figura 4. Clicca qui per scaricare questa tabella.
Fascicolo Supplementare 1. File di stampa 3D per il connettore patchcord.Clicca qui per scaricare questo file.
File supplementare 2. 3D file di stampa per connettore patch senza fori passanti.Clicca qui per scaricare questo file.
File supplementare 3. 3D file di stampa per accessori di perforazione stereotatsica.Clicca qui per scaricare questo file.
Fascicolo Supplementare 4. Progetto personalizzato del disco per lucidatura.Clicca qui per scaricare questo file.
File supplementare 5. 3D file di stampa per array di fotometria multisito che mira a 12 siti SBN.Clicca qui per scaricare questo file.
File supplementare 6. 3D file di stampa per array di fotometria multisito senza fori passanti che puntano i siti.Clicca qui per scaricare questo file.
File Supplementare 7. File di configurazione DAQ per il software DAQ.Clicca qui per scaricare questo file.
Fascicolo Supplementare 8. File di codifica visiva per il sistema di fotometria.Clicca qui per scaricare questo file.
File Supplementare 9. Il codice Python utilizzato nel Protocollo Step 4 e per i confronti corretti con artefatti di movimento. Clicca qui per scaricare questo file.