$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Para grabar simultáneamente desde hasta 12 regiones cerebrales de SBN, construimos una matriz multifibra personalizable siguiendo el protocolo anterior. Tras la recuperación de la cirugía, los ratones fueron registrados durante la interacción social (Figura 5). En todos los experimentos, se permitió que los ratones machos CD1 interactuaran con ratones Balb/c machos en su jaula de origen. Utilizamos ratones Vgat-Cre para expresar jRCaMP1b40 en neuronas GABAérgicas de manera dependiente de Cre en 11 regiones cerebrales de SBN, GCaMP6f39 en células no GABAérgicas de manera excluyentede Cre 45 en las mismas regiones cerebrales, y GRAB-DA1m33 en el núcleoaccumbens 42 (Tabla 1). La confirmación histológica de la segmentación de fibras puede encontrarse en nuestras publicaciones anteriores42,43. Aquí, buscamos probar el rendimiento de diferentes métodos para la corrección de artefactos de movimiento. Específicamente, comparamos la eficacia de usar: (1) un filtro mediano solo; (2) Regresión de lazo con regularización L1 usando únicamente señales isosmísicas de la misma fibra (enfoque de predictor único)11 y (3) Regresión de lazo usando señales isosmísicas de todas las fibras como predictores (enfoque multipredictor)41,42 (Figura 4A). De manera crítica, los métodos basados en lazo corrigen la señal en bruto usando la señal fluorescente independiente de la actividad generada en respuesta al LED de frecuencia isosmósica (415 nm), mientras que el suavizado del filtro mediano no considera la señal independiente de la actividad.
Evaluamos la corrección de movimiento utilizando dos medidas diferentes: ruido absoluto de alta frecuencia y jitter. El ruido de alta frecuencia se definió como la desviación estándar de la primera derivada de la señal, que refleja la magnitud de las rápidas fluctuaciones de la señal. El jitter se definió como el coeficiente de variación de la derivada de la señal, que refleja la magnitud de la variabilidad de la señal de trama a trama respecto a la media de la derivada de la señal. Descubrimos que ambos enfoques basados en Lazo (predictor único y multipredictor) superaron sustancialmente al filtro mediano en términos de reducción tanto del ruido de alta frecuencia como del jitter, con la corrección multipredictor superando ligeramente la corrección de predictor único (Figura 4C–D; Ruido de alta frecuencia, n=414 para 23 señales en 18 sesiones de grabación: prueba de Kruskal-Wallis, H = 501,9420, p < 0,001; Pruebas U de Mann-Whitney post-hoc con corrección de Bonferroni: Regresión suavizada vs de un solo sitio ***p < 0,001, Regresión suavizada vs multisitio***p < 0,001, Regresión monositio vs Multi-sitio **p = 0,0073; Jitter: Test de Kruskal-Wallis, H = 495.6685, p < 0.001; pruebas U de Mann-Whitney post-hoc con corrección de Bonferroni: Regresión suavizada vs de sitio único ***p < 0,001, Regresión suavizada vs multisitio***p < 0,001, Regresión monositio vs multisitio ***p = 0,0006). Específicamente, el enfoque multipredictor (utilizando todas las isosmísticas de fibra) mostró beneficios para las señales de canal rojo (RCaMP): 5 de 12 señales rojas mostraron un ruido residual significativamente menor en comparación con la corrección con predictor único (n = 18 grabaciones/región). Para el jitter, BNST (prueba t, t = 5,42, ***p < 0,001), POA (prueba t, t = 2,68, *p = 0,011), AH (prueba t, t = 2,66, *p = 0,012) y PAG (prueba U de Mann-Whitney, U = 279, ***p < 0,001) mostraron superioridad en múltiples predictores. Para ruido de alta frecuencia, BNST (prueba U de Mann-Whitney, U = 285, ***p < 0,001), POA (prueba U de Mann-Whitney, U = 292, ***p < 0,001), AH (prueba t, t = 7,10, ***p < 0,001), PAG (prueba t, t = 13,72, ***p < 0,001) y LHb (prueba t, t = 5,47, ***p < 0,001) mostraron superioridad multipredictora. Las señales de canal verde (GCaMP) no mostraron diferencias significativas entre los dos enfoques Lasso en las 12 regiones para métricas de jitter o ruido de alta frecuencia (todas p>0,05, para la lista completa de pruebas estadísticas véase Tabla 3). Esto sugiere que, para la fotometría de doble color, incorporar información isosmética cruzada mejora la corrección de movimiento para señales rojas, probablemente porque las isosmésticas de canal verde proporcionan información de movimiento más limpia. Estos datos demuestran que la corrección de artefactos de movimiento mediante métodos de regresión basados en isosbest supera ampliamente a los enfoques de filtrado simples, siendo la regresión multicanal una ventaja específica para grabaciones de doble color.

Figura 1. Grabación de la Red de Comportamiento Social utilizando fotometría de fibra multisitio y bicolor. (A) Representación esquemática de cómo construir un arreglo multifibra. (B) Representación esquemática del flujo de trabajo quirúrgico. (C) Representación esquemática de experimentos de imagen de calcio a dos colores. (D) Ejemplo de usar una señal isoscómica para calcular una señal corregida por movimiento. Por favor, haz clic aquí para ver una versión ampliada de esta figura.

Figura 2. Imágenes de procedimientos clave para fabricación y cirugía. (A) Cable de parche fabricado listo para grabación. (B) Matriz multifibra colocada en un soporte estereotáctico personalizado con la broca de 0,25 mm x 5,5 mm situada sobre un orificio pasante de fibra. (C) Matriz multifibra con todos los agujeros pasantes de fibra despejados y tuercas hexagonales pegadas. (D) Cortar una fibra usando un grabador de fibra. (E) Matriz multifibra con todas las fibras pegadas con superpegamento y colocadas sobre un disco de pulido personalizado. (F) Vista ampliada de cómo alinear y atornillar el conector del cable de parche a la matriz multifibra. (G) Matriz multifibra completamente construida con clavillas pegadas en su lugar. (H) Ratón fijado a un marco estereotáctico con incisión en el cuero cabelludo hecha para exponer el cráneo. (I) Matriz multifibra fabricada acoplada a un manipulador estereotáctico para maniobrar la matriz multifibra durante la cirugía. (J) Vista ampliada de la matriz multifibra siendo bajada al cerebro. (K) Matriz multifibra completamente bajada y fijada al cráneo con cemento dental. (L) Matriz multifibra en un animal que está conectada de forma segura al parche antes de la grabación. (M) Vista alejada del animal en jaula con cable conectado antes y/o durante la grabación. Por favor, haz clic aquí para ver una versión ampliada de esta figura.

Figura 3. Separación espectral. Ejemplos de trazas que demuestran una diafonía espectral mínima entre indicadores de calcio. Las señales GCaMP6f (azul) muestran una activación robusta durante la excitación de 470 nm pero una respuesta insignificante a la excitación de 560 nm. Las señales jRCaMP1b (rojas) presentan el patrón inverso con una fuerte activación desde 560 nm pero una respuesta mínima a una iluminación de 470 nm. (A) Simultáneamente grabó señales GCaMP6f y jRCaMP1b con iluminación de 470 nm. (B) Igual que (a), pero solo con iluminación de 560 nm. Por favor, haz clic aquí para ver una versión ampliada de esta figura.

Figura 4. Comparación de diferentes métodos para el procesamiento de señales multifibra. (A) Esquema que representa la corrección de la señal de fotometría basada en regresión multisitio. Este enfoque implica pasos secuenciales de procesamiento en los que todas las señales isosmísicas disponibles de 415nm se mapean linealmente a cada señal dependiente del calcio. Los componentes del movimiento predicho se restan a esta señal para producir la lectura final de la actividad. señales ΔF/F con puntuación z de dos señales con fuerte presencia de artefactos de movimiento. (B) De arriba hacia abajo: señal isosmódica del mismo canal de fibra/imagen que representa artefactos de movimiento, señal RCaMP solo suavizada (izquierda) o GCaMP (derecha) (filtrada por mediana, verde), misma señal corregida por regresión de fibra/canal de imagen (cian), señal corregida por regresión multifibra (azul). (C) Izquierda: Puntuación media de ruido de alta frecuencia para cada señal grabada (n = 21 sesiones). El ruido de alta frecuencia se calcula como la desviación estándar de la primera derivada de cada señal. Se muestran las puntuaciones de ruido para 1) señal corregida por línea base pero no corregida (gris), 2) señal corregida y filtrada por mediana (rosa), 3) señal corregida por regresión de línea base y de fibra única con suavizado Savitsky-Golay (cian), y 4) señal corregida por regresión basal y multifibra con suavizado Savitsky-Golay (azul). Centro: Porcentaje medio de ruido eliminado usando solo suavizado, regresión de fibra única o multifibra, agregado en las 23 señales. El porcentaje de ruido eliminado se calcula como la fracción de ruido de alta frecuencia resultante de la corrección respecto al nivel general de ruido en la señal bruta corregida por línea base. Los datos se informan como media ± SEM (prueba de Kruskal-Wallis, H = 501,9420, p < 0,001; pruebas post-hoc de Mann-Whitney U con corrección de Bonferroni: Regresión suavizada vs de un solo sitio ***p < 0,001, Regresión suavizada vs multisitio***p < 0,001, Regresión monositio vs multisitio **p = 0,0073). Derecha: porcentaje de señales en las que el ruido de alta frecuencia es estadísticamente diferente tras una corrección basada en regresión de una sola fibra frente a multifibra, separadas por tipo de señal. Para cada una de las 23 señales, los niveles de ruido tras la corrección en un solo sitio se compararon con los niveles de ruido tras la corrección multi-sitio utilizando pruebas t independientes (si se distribuyen normalmente) o pruebas Mann-Whitney U (alternativa no paramétrica). Los segmentos azules indican señales donde la corrección multisitio resultó en un ruido significativamente menor que la corrección en un solo sitio (p < 0,05). Los segmentos grises indican señales sin diferencias significativas entre los métodos de corrección. Los segmentos cian indican señales donde la corrección en un solo sitio resultó en un ruido significativamente menor. E: Vgat- GCaMP6f. DA: AGARRA-DA. I: Vgat+ jRCaMP1b. (D) Igual que (C) pero representando la variabilidad de cuadro a cuadro (jitter), calculada como la derivada del coeficiente de variación de cada señal (prueba de Kruskal-Wallis, H = 495,6685, p < 0,001; pruebas Mann-Whitney U post-hoc con corrección de Bonferroni: Regresión suavizada vs de un solo sitio ***p < 0,001, Regresión suavizada vs multisitio***p < 0,001, Un solo sitio vs multisitio ***p = 0,0006). Para la lista completa de pruebas estadísticas y resultados, consulte la Tabla 3. Por favor, haz clic aquí para ver una versión ampliada de esta figura.

Figura 5. Exposición multisitio de fotometría de fibra a doble color. Señales representativas de varios tipos de indicadores registradas simultáneamente durante el comportamiento social (prueba de residente-intruso, época de entrada del intruso). Las trazas muestran GCaMP6f no normalizada y restada de la línea basal pero no corregida por artefactos de movimiento (neuronas Cre-negativas), jRCaMP1b (neuronas GABAérgicas Vgat-Cre+) y GRAB-DA1m (sensor de dopamina en núcleo accumbens) en los sitios de la red de comportamiento social en el mismo animal. Por favor, haz clic aquí para ver una versión ampliada de esta figura.
Tabla 1. Coordenadas del virus y longitudes de fibra para cada región cerebral de la Red de Comportamiento Social Objetivo. Por favor, haz clic aquí para descargar esta tabla.
Tabla 2. Comparación de métodos de corrección de artefactos de movimiento para fotometría de fibra multisitio. Principios subyacentes, casos de uso óptimos, ventajas y limitaciones de los siguientes métodos de corrección de artefactos de movimiento: regresión lineal simple, regresión Lasso con regularización L1 (usando enfoque multipredictor), modelado generativo bayesiano, escalado basado en frecuencia y mínimos cuadrados reponderados iterativamente. Por favor, haz clic aquí para descargar esta tabla.
Tabla 3. Detalles y resultados de las pruebas estadísticas utilizadas en la Figura 4. Por favor, haz clic aquí para descargar esta tabla.
Expediente suplementario 1. Archivo de impresión 3D para conector de patchcable.Por favor, haga clic aquí para descargar este archivo.
Archivo suplementario 2. 3D archivo de impresión para conector de cable de parche sin agujeros pasantes.Por favor, haga clic aquí para descargar este archivo.
Archivo suplementario 3. 3D archivo de impresión para accesorio de perforación estereotáxica.Por favor, haga clic aquí para descargar este archivo.
Expediente suplementario 4. Plano personalizado del puck de pulido.Por favor, haga clic aquí para descargar este archivo.
Archivo suplementario 5. 3D archivo de impresión para matriz de fotometría multisitio dirigida a 12 sitios SBN.Por favor, haga clic aquí para descargar este archivo.
Archivo suplementario 6. 3D archivo de impresión para matriz de fotometría multisitio sin atravesar los sitios dirigidos.Por favor, haga clic aquí para descargar este archivo.
Expediente suplementario 7. Archivo de configuración DAQ para software DAQ.Por favor, haga clic aquí para descargar este archivo.
Expediente suplementario 8. Archivo de codificación visual para el sistema de fotometría.Por favor, haga clic aquí para descargar este archivo.
Expediente suplementario 9. Código en Python usado en el Protocolo Step 4 y para comparaciones corregidas de artefactos de movimiento. Por favor, haga clic aquí para descargar este archivo.