$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Para gravar simultaneamente de até 12 regiões cerebrais de SBN, construímos uma matriz multifibra personalizável seguindo o protocolo acima. Após a recuperação da cirurgia, os camundongos foram registrados durante a interação social (Figura 5). Para todos os experimentos, camundongos machos CD1 foram autorizados a interagir com ratos Balb/c machos em sua gaiola doméstica. Usamos camundongos Vgat-Cre para expressar jRCaMP1b40 em neurônios GABAérgicos de forma dependente de Cre em 11 regiões cerebrais de SBN, GCaMP6f39 em células não GABAérgicas de forma excluidade Cre 45 nas mesmas regiões cerebrais, e GRAB-DA1m33 no Núcleo accumbens42 (Tabela 1). A confirmação histológica do direcionamento das fibras pode ser encontrada em nossas publicações anteriores42,43. Aqui, buscamos testar o desempenho de diferentes métodos para correção de artefatos de movimento. Especificamente, comparamos a eficácia do uso: (1) apenas um filtro mediano; (2) Regressão de laço com regularização L1 usando apenas sinais isósméticos da mesma fibra (abordagem de preditor único)11 e (3) Regressão de laço usando sinais isosmíticos de todas as fibras como preditores (abordagem multi-preditor)41,42 (Figura 4A). De forma crítica, os métodos baseados em Laço corrigem o sinal bruto usando o sinal fluorescente independente da atividade gerado em resposta ao LED de frequência isosmítica (415 nm), enquanto o suavizamento do filtro mediano não considera o sinal independente da atividade.
Avaliamos a correção de movimento usando duas medidas diferentes: ruído absoluto de alta frequência e jitter. Ruído de alta frequência foi definido como o desvio padrão da primeira derivada do sinal, que reflete a magnitude das flutuações rápidas do sinal. Jitter foi definido como o coeficiente de variação da derivada do sinal, que reflete a magnitude da variabilidade do sinal quadro a quadro em relação à média da derivada do sinal. Descobrimos que ambas as abordagens baseadas em Laço (preditor único e multipreditor) superaram substancialmente o filtro mediano em termos de redução tanto do ruído de alta frequência quanto do jitter, com a correção multipreditor superando marginalmente a correção de preditor único (Figura 4C–D; Ruído de alta frequência, n=414 para 23 sinais em 18 sessões de gravação: teste de Kruskal-Wallis, H = 501,9420, p < 0,001; Testes pós-hoc de Mann-Whitney U com correção de Bonferroni: Regressão suavizada vs de sítio único ***p < 0,001, Regressão suavizada vs multi-sítio ***p < 0,001, Regressão mono-sítio vs Multi-sítio **p = 0,0073; Jitter: teste de Kruskal-Wallis, H = 495,6685, p < 0,001; testes pós-hoc de Mann-Whitney U com correção de Bonferroni: Regressão suavizada vs de sítio único ***p < 0,001, Regressão suavizada vs multi-sítio ***p < 0,001, Regressão mono-sítio vs Multi-sítio ***p = 0,0006). Especificamente, a abordagem multi-preditor (usando toda a isosbest de fibra) mostrou benefícios para sinais de canal vermelho (RCaMP): 5 de 12 sinais vermelhos apresentaram ruído residual significativamente menor em comparação com correção preditor único (n = 18 gravações/região). Para jitter, BNST (teste t, t = 5,42, ***p < 0,001), POA (teste t, t = 2,68, *p = 0,011), AH (teste t, t = 2,66, *p = 0,012) e PAG (teste U de Mann-Whitney, U = 279, ***p < 0,001) mostraram superioridade em múltiplos preditores. Para ruído de alta frequência, BNST (teste U de Mann-Whitney, U = 285, ***p < 0,001), POA (teste U de Mann-Whitney, U = 292, ***p < 0,001), AH (teste t, t = 7,10, ***p < 0,001), PAG (teste t, t = 13,72, ***p < 0,001) e LHb (teste t, t = 5,47, ***p < 0,001) mostraram superioridade em múltiplos preditores. Os sinais de canal verde (GCaMP) não mostraram diferenças significativas entre as duas abordagens do Lasso em todas as 12 regiões para métricas de jitter ou ruído de alta frequência (todas p>0,05, para a lista completa de testes estatísticos veja Tabela 3). Isso sugere que, para fotometria de dupla cor, incorporar informações isosbesticas cruzadas melhora a correção de movimento para sinais vermelhos, provavelmente porque isosmêsticas de canal verde fornecem informações de movimento mais limpas. Esses dados demonstram que a correção de artefatos de movimento usando métodos de regressão baseados em isosbest supera substancialmente as abordagens simples de filtragem, com a regressão multicanal oferecendo vantagens específicas para gravações em cores duplas.

Figura 1. Gravação da Rede de Comportamento Social usando fotometria de fibra multisite e duas cores. (A) Representação esquemática de como construir um array multifibra. (B) Representação esquemática do fluxo de trabalho cirúrgico. (C) Representação esquemática de experimentos de imagem de cálcio em duas cores. (D) Exemplo de usar sinal isósbico para calcular um sinal corrigido por movimento. Por favor, clique aqui para ver uma versão ampliada desta figura.

Figura 2. Imagens de procedimentos-chave para fabricação e cirurgia. (A) Cabo de conexão fabricado pronto para uso em gravação. (B) Matriz multifibra colocada em um suporte estereotáxico personalizado com a broca de 0,25 mm x 5,5 mm posicionada acima de um furo de fibra atravessante. (C) Matriz multifibra com todos os furos de passagem da fibra limpos e porcas hexagonales coladas. (D) Cortar uma fibra usando um gravador de fibra. (E) Matriz multifibra com todas as fibras coladas com supercola e colocadas em um disco de polimento personalizado. (F) Vista ampliada de como alinhar e rosquear o conector do cabo de conexão ao array multifibra. (G) Matriz multifibra totalmente construída com pinos de cavilha colados no lugar. (H) Rato fixado à estrutura estereotáxica com incisão no couro cabeludo feita para expor o crânio. (I) Matriz multifibra fabricada acoplada a um manipulador estereotáxico para manobrar a matriz multifibra durante a cirurgia. (J) Vista ampliada da matriz multifibra sendo baixada para dentro do cérebro. (K) Matriz multifibra totalmente abaixada e fixada ao crânio com cimento dentário. (L) Matriz multifibra em um animal que esteja conectada de forma segura ao cabo de conexão antes da gravação. (M) Vista ampliada do animal na gaiola com cabo conectado antes e/ou durante a gravação. Por favor, clique aqui para ver uma versão ampliada desta figura.

Figura 3. Separação espectral. Exemplos de traços demonstrando uma interferência espectral mínima entre indicadores de cálcio. Os sinais GCaMP6f (azul) apresentam ativação robusta durante excitação de 470 nm, mas resposta desprezível a excitação de 560 nm. Os sinais jRCaMP1b (vermelho) apresentam o padrão inverso com forte ativação a partir de 560 nm, mas resposta mínima a iluminação de 470 nm. (A) Registrava simultaneamente sinais GCaMP6f e jRCaMP1b com iluminação apenas de 470 nm. (B) Igual a (a), mas apenas com iluminação de 560 nm. Por favor, clique aqui para ver uma versão ampliada desta figura.

Figura 4. Comparação de diferentes métodos para processamento de sinais multifibra. (A) Esquema representando a correção baseada em regressão multi-site do sinal de fotometria. Essa abordagem envolve etapas sequenciais de processamento nas quais todos os sinais isosmíticos disponíveis de 415nm são linearmente mapeados para cada sinal dependente do cálcio. Componentes de movimento previstos são subtraídos desse sinal para produzir a leitura final da atividade. Sinais ΔF/F com pontuação z provenientes de dois sinais com forte presença de artefatos de movimento. (B) De cima a baixo: sinal isosmêstico do mesmo canal de fibra/imagem representando artefatos de movimento, sinal apenas suavizado RCaMP (esquerda) ou GCaMP (direita) (filtrado com mediana, verde), mesmo sinal de fibra/canal de imagem corrigido por regressão (ciano), sinal de regressão multi-fibra corrigido (azul). (C) Esquerda: Pontuação média de ruído de alta frequência para cada sinal gravado (n = 21 sessões). Ruído de alta frequência é calculado como o desvio padrão da primeira derivada de cada sinal. As pontuações de ruído são mostradas para 1) sinal corrigido pela linha de base mas não corrigido (cinza), 2) sinal corrigido pela linha de base e filtrado pela mediana (rosa), 3) sinal corrigido por regressão de linha de base e de fibra única com suavização Savitsky-Golay (ciano), e 4) sinal corrigido por regressão de linha de base e multifibra com suavização Savitsky-Golay (azul). Centro: Porcentagem média de ruído removida usando apenas suavização, regressão de fibra única ou multifibra, agregada em todos os 23 sinais. A porcentagem de ruído removido é calculada como a fração do ruído de alta frequência resultante da correção em relação ao nível geral de ruído no sinal bruto corrigido pela linha de base. Os dados são reportados como média ± SEM (teste de Kruskal-Wallis, H = 501,9420, p < 0,001; testes pós-hoc de Mann-Whitney U com correção de Bonferroni: Regressão suavizada vs de sítio único ***p < 0,001, regressão suavizada vs multi-sítio***p < 0,001, Regressão mono-local vs multi-sítio **p = 0,0073). Direita: Porcentagem de sinais em que o ruído de alta frequência é estatisticamente diferente após correção baseada em regressão de fibra única versus multifibra, separada por tipo de sinal. Para cada um dos 23 sinais, os níveis de ruído após correção em local único foram comparados aos níveis de ruído após correção multi-local, usando testes t independentes (se normalmente distribuídos) ou testes Mann-Whitney U (alternativa não paramétrica). Segmentos azuis indicam sinais onde a correção multi-site resultou em ruído significativamente menor do que a correção em local único (p < 0,05). Segmentos cinza indicam sinais sem diferença significativa entre os métodos de correção. Segmentos ciano indicam sinais onde a correção em um único local resultou em ruído significativamente menor. E: Vgat- GCaMP6f. DA: AGARRE-DA. I: Vgat+ jRCaMP1b. (D) Igual a (C), mas representando variabilidade quadro a quadro (jitter), calculada como a derivada do coeficiente de variação de cada sinal (teste de Kruskal-Wallis, H = 495,6685, p < 0,001; testes pós-hoc de Mann-Whitney U com correção de Bonferroni: Regressão suavizada vs de sítio único ***p < 0,001, Regressão suavizada vs multi-sítio***p < 0,001, Mono-sítio vs Multi-sítio ***p = 0,0006). Para a lista completa de testes e resultados estatísticos, consulte a Tabela 3. Por favor, clique aqui para ver uma versão ampliada desta figura.

Figura 5. Vitrine multi-site de fotometria de fibra dupla. Sinais representativos de múltiplos tipos de indicadores foram registrados simultaneamente durante o comportamento social (teste residente-intruso, época de entrada do intruso). Os traços mostram GCaMP6f (neurônios Cre-negativos), jRCaMP1b (neurônios Vgat-Cre+ GABAérgicos e subtraídos no ponto de partida), mas não corrigidos por artefatos de movimento, GCaMP1b (neurônios Vgat-Cre+ GABAérgicos) e GRAB-DA1m (sensor de dopamina no núcleo accumbens) em vários locais de rede de comportamento social no mesmo animal. Por favor, clique aqui para ver uma versão ampliada desta figura.
Tabela 1. Coordenadas do vírus e comprimentos das fibras para cada Rede de Comportamento Social Alvo da Região do Cérebro. Por favor, clique aqui para baixar esta Tabela.
Tabela 2. Comparação dos métodos de correção de artefatos de movimento para fotometria de fibra multi-site. Princípios subjacentes, casos de uso ótimos, vantagens e limitações dos seguintes métodos de correção de artefatos de movimento: regressão linear simples, regressão de Lasso com regularização L1 (usando abordagem multipreditor), modelagem generativa bayesiana, escalonamento baseado em frequência e mínimos quadrados reponderados iterativamente. Por favor, clique aqui para baixar esta Tabela.
Tabela 3. Detalhes e resultados dos testes estatísticos usados na Figura 4. Por favor, clique aqui para baixar esta Tabela.
Arquivo Suplementar 1. Arquivo de impressão 3D para conector de cabos de conexão.Por favor, clique aqui para baixar este arquivo.
Arquivo suplementar 2. 3D arquivo de impressão para conector de cabo de conexão sem furos passantes.Por favor, clique aqui para baixar este arquivo.
Arquivo suplementar 3. 3D arquivo de impressão para acessório de perfuração estereotáxica.Por favor, clique aqui para baixar este arquivo.
Arquivo Suplementar 4. Planta personalizada do puck de polimento.Por favor, clique aqui para baixar este arquivo.
Arquivo suplementar 5. 3D arquivo de impressão para matriz de fotometria multisite direcionada a 12 sites SBN.Por favor, clique aqui para baixar este arquivo.
Arquivo suplementar 6. 3D arquivo de impressão para matriz de fotometria multisite sem atravessar buracos direcionando os sítios.Por favor, clique aqui para baixar este arquivo.
Arquivo Suplementar 7. Arquivo de configuração DAQ para software DAQ.Por favor, clique aqui para baixar este arquivo.
Arquivo Suplementar 8. Arquivo de codificação visual para o sistema de fotometria.Por favor, clique aqui para baixar este arquivo.
Arquivo Suplementar 9. Código Python usado no Protocolo Step 4 e para comparações corrigidas de artefatos de movimento. Por favor, clique aqui para baixar este arquivo.