Method Article
O estudo visa desenvolver tecnologia para medições de batimentos cardíacos sem anestesia em peixes-zebra em movimento. Nossa abordagem combina imagens infravermelhas de ondas curtas e rastreamento do coração baseado em aprendizado de máquina. É uma técnica não invasiva, sem rótulos e fácil de usar que se adapta a uma ampla gama de estudos sobre o modelo de peixe-zebra.
O peixe-zebra (Danio rerio) é um organismo modelo amplamente utilizado em pesquisas fisiológicas, farmacológicas e toxicológicas devido à sua semelhança genética com humanos e estágio embrionário transparente, o que facilita estudos cardiovasculares não invasivos. No entanto, os métodos atuais de avaliação da frequência cardíaca em peixe-zebra muitas vezes dependem da anestesia para imobilizar o sujeito, introduzindo alterações fisiológicas que comprometem a precisão e a reprodutibilidade dos dados. Este estudo apresenta uma nova técnica sem anestesia para medir os batimentos cardíacos em larvas de peixe-zebra em movimento livre, abordando uma limitação crítica na pesquisa cardiovascular. A abordagem proposta integra imagens infravermelhas de ondas curtas com rastreamento cardíaco baseado em aprendizado de máquina, permitindo o monitoramento preciso e contínuo da atividade cardíaca em espécimes não imobilizados. Uma rede neural convolucional foi treinada para detectar a região do coração e um sinal fotopletismográfico foi extraído das sequências de imagens para determinar a frequência cardíaca. A validação experimental demonstrou a confiabilidade e consistência do método em várias condições de teste. Um dos principais benefícios da metodologia é sua capacidade de preservar o estado fisiológico natural do peixe-zebra, minimizando os artefatos induzidos pelo estresse. Essa técnica não invasiva e sem rótulo oferece vantagens significativas para o estudo da fisiologia cardiovascular, cardiotoxicidade de drogas e toxicologia ambiental, expandindo as aplicações potenciais do peixe-zebra como modelo para pesquisa biomédica.
O peixe-zebra (Danio rerio), um pequeno peixe ciprinídeo, tornou-se um organismo modelo essencial devido ao seu pequeno tamanho, alta taxa reprodutiva e facilidade de manipulação genética 1,2,3. A avaliação da freqüência cardíaca em embriões transparentes de peixe-zebra é cada vez mais utilizada em fisiologia, embriologia, toxicologia e outros campos 4,5,6,7,8. Por um lado, essa utilidade se deve ao fato de que o genoma do peixe-zebra inclui genes associados a doenças cardiovasculares humanas9, e o coração Danio rerio compartilha estruturas e vias de sinalização semelhantes com os humanos10,11. Isso torna o peixe-zebra um modelo inestimável para estudar o desenvolvimento e as doenças do coração 11,12,13. Por outro lado, a frequência cardíaca do peixe-zebra é sensível a influências externas, tornando-se um excelente modelo para estudos fisiológicos e toxicológicos, comparando a função cardíaca em peixes tratados e não tratados 7,8,14.
Progressos significativos foram feitos no desenvolvimento de métodos ópticos não invasivos para avaliar a frequência cardíaca em embriões de peixes transparentes15,16. Essas técnicas oferecem a vantagem de uma coleta rápida de dados de grandes tamanhos de amostra. Consequentemente, abordagens totalmente automatizadas para avaliação da frequência cardíaca em embriões de peixes foram desenvolvidas 4,5,6,17.
No entanto, certas limitações atualmente restringem o uso dessas técnicas ao período de 3-4 dpf. A primeira limitação é a perda de transparência devido à pigmentação do corpo do peixe. A segunda é o aumento da mobilidade dos embriões ao longo do tempo. Estender o período de desenvolvimento inicial do peixe-zebra durante o qual as abordagens ópticas podem ser usadas aumentaria sua utilidade, permitindo projetos experimentais de longo prazo para estudar cardiomiopatia, defeitos cardíacos congênitos e vários impactos no sistema cardiovascular, incluindo o rastreamento da dinâmica dos efeitos ao longo do tempo. Nosso grupo recentemente abordou a questão da perda de transparência empregando imagens na faixa de infravermelho de ondas curtas de 900-1700 nm18. Este artigo se concentra em abordar a questão da mobilidade embrionária.
Normalmente, anestésicos como o metanossulfonato de tricaína (MS-222) são usados para imobilizar embriões e larvas de peixes nadando livremente antes da imagem 14,19,20. No entanto, o MS-222 reduz significativamente a freqüência cardíaca21,22, assim como outros anestésicos23. Torna-se um desafio discernir se as alterações observadas na função cardíaca são devidas ao tratamento experimental, ao anestésico ou a uma interação entre os dois. Outra maneira de estender o período de baixa mobilidade do embrião é reduzindo a temperatura durante o desenvolvimento inicial8. No entanto, essa abordagem nem sempre é compatível com os objetivos da pesquisa e apenas estende minimamente o período de inscrição.
Neste estudo, apresentamos um novo método para abordar a mobilidade embrionária durante o registro da frequência cardíaca. Treinamos uma rede neural convolucional para identificar a região de interesse do coração em registros de embriões de peixe-zebra nadando livremente. A variação periódica na intensidade do pixel dentro desta área é utilizada para derivar o fotopletismograma (PPG), que é posteriormente usado para calcular a frequência cardíaca. O aplicativo de desktop desenvolvido, AutoHR, utilizou treinamento de rede neural e processamento de pilha de imagens, garantindo facilidade de uso e reprodutibilidade do protocolo.
Os peixes-zebra foram criados e criados de acordo com os protocolos ZFIN estabelecidos24. Todos os procedimentos foram aprovados pelo Comitê de Bioética do Centro Científico e Tecnológico de Instrumentação Única da Academia Russa de Ciências (STC UI RAS), protocolo #3/24, datado de 21/08/2024, e seguem as diretrizes de cuidados com o peixe-zebra do STC UI RAS. Manuais para versões individuais estão disponíveis mediante solicitação.
1. Preparação do equipamento para medição
2. Aquisição de imagem
3. Treinando a rede neural para rotulagem
4. Treinando as redes neurais para detecção cardíaca
NOTA: Esta etapa é realizada uma vez para uma idade específica e condição de imagem. A GPU NVIDIA é fortemente recomendada para treinamento, pois acelera significativamente o processamento.
5. Quantificação da frequência cardíaca
6. Verificação dos resultados do algoritmo
A frequência cardíaca do peixe-zebra a 12 dpf foi determinada usando o protocolo descrito acima (Vídeo Suplementar 1). Os vídeos incluem uma sequência de imagens de larvas de peixe-zebra nadando livremente, um fotopletismograma derivado dessas sequências usando o protocolo proposto e a frequência cardíaca correspondente calculada a partir do fotopletismograma.
Os dados rotulados foram divididos aleatoriamente em conjuntos de treinamento, teste e validação em uma proporção de 3:2:1 durante o treinamento. A função de perda foi então implementada de acordo com a seguinte expressão:
onde x é logits de entrada, y é um alvo, tamanho do lote é o número de imagens processadas simultaneamente,
Para avaliar o desempenho do modelo durante o treinamento, foi empregada a função de precisão de predição, definida pela seguinte:
Aplicamos abordagens padrão para o desenvolvimento da função de perda, conforme detalhado, por exemplo, em25. A rede de atenção multiescala com arquitetura de modelo EfficientNet-B0 foi usada como codificador e camada de saída softmax. A rede neural foi treinada em uma amostra de teste por 50 épocas, utilizando uma taxa de aprendizado exponencialmente decrescente que começou em 0,02 e teve um fator de decaimento de 0,99.
A rede neural treinada foi usada para obter a máscara corporal e a máscara ocular (Figura 5С) do indivíduo em cada quadro de uma sequência arbitrária. O centro de massa da máscara ocular foi então calculado e a esqueletização foi aplicada à máscara corporal (Figura 5D). No cordão de peixe P resultante, foi identificado o ponto M, localizado mais próximo do centro de massa da máscara ocular (Figura 5F). A uma distância d do ponto M, igual à distância entre os olhos, o ponto N foi identificado ao longo da corda PPP em direção à cauda. Este ponto corresponde ao centro da área do coração. Essa área foi delimitada por um círculo com raio de r = d x 0,3, centralizado nesse ponto (Figura 5F).
O sinal PPG bruto foi calculado como o número médio de pixels dentro da área do coração em todos os quadros gravados (Figura 5G). Este sinal bruto foi então normalizado para seu valor médio. Posteriormente, o sinal foi filtrado usando um filtro Chebyshev Tipo II de segunda ordem com um design de4ª ordem, atenuação da banda de parada de 25 dB e uma frequência de borda da banda de parada de 12 Hz, consistente com os algoritmos de processamento padrão26,27 (Figura 6A). A frequência cardíaca foi determinada identificando-se a frequência dentro da faixa de 1-4 Hz que corresponde à potência de pico no módulo quadrado da transformada de Fourier do sinal filtrado (Figura 6A).
A equação (1) garante que a imagem do coração tenha aproximadamente 20 pixels x 20 pixels. Supõe-se que a dimensão longitudinal do coração seja responsável por aproximadamente 8% do comprimento total do corpo do peixe-zebra. Quando a imagem do coração é menor, o sinal resultante geralmente contém um componente de ruído significativo que supera o sinal útil, tornando-o inadequado para o cálculo da frequência cardíaca, conforme mostrado na Figura 6B.
Os resultados da aplicação do protocolo dependem significativamente da estabilidade da segmentação cardíaca para cada indivíduo. É crucial realizar experimentos de controle para demonstrar que as etapas 2 a 5 foram concluídas com sucesso e os dados obtidos refletem processos fisiológicos reais, em vez de resultados aleatórios. A validação do protocolo foi realizada seguindo o procedimento detalhado na etapa 6. No teste descrito na etapa 6.2, o sal foi adicionado à água a uma concentração de 5 mg/L como estímulo provocativo. Cinco larvas de peixe-zebra (12 pdf) foram incluídas em cada teste. A temperatura da água em um molde de agarose foi mantida em 23-25 °C. Os resultados indicam que a diferença entre a avaliação visual e o método proposto não ultrapassou 3% no primeiro teste. No segundo teste, a frequência cardíaca medida pelo método aumentou, como esperado (Figura 7).
Figura 1: Fluxo de trabalho esquemático. (A) Esquema e (B) aparência da configuração experimental Clique aqui para ver uma versão maior desta figura.
Figura 2: Imagens de larvas de peixe-zebra. Exemplos de (A) imagens de larvas de peixe-zebra de alta qualidade e (B-D) de baixa qualidade são (B) quadro subexposto, (С) quadro superexposto e (D) quadro desfocado. As imagens foram cortadas para destacar a larva. A barra de escala é de 1 mm. Clique aqui para ver uma versão maior desta figura.
Figura 3: Imagens representativas. Exemplos de (A) imagens iniciais e (B) rotuladas. A barra de escala é de 1 mm. Clique aqui para ver uma versão maior desta figura.
Figura 4: Contornos do corpo do peixe-zebra. Exemplos de (A) contornos de alta qualidade e (B-C) de baixa qualidade de todo o corpo do peixe-zebra: um peixe-zebra não se encaixa inteiramente no retângulo e o tamanho do retângulo é visivelmente maior que o do peixe-zebra. A barra de escala é de 1 mm. Clique aqui para ver uma versão maior desta figura.
Figura 5: Pipeline de processamento de imagem ilustrando os estágios de processamento primário e seus resultados intermediários correspondentes. Clique aqui para ver uma versão maior desta figura.
Figura 6: Quadros mostrando a frequência cardíaca larval. Exemplos de (A) quadros de alta qualidade e (B) de baixa qualidade com a área do coração larval com diâmetros de 20 pixels e 14 pixels, respectivamente, e espectros PPG e Fourier correspondentes de PPG. A barra de escala é de 1 mm. Abreviaturas: PPG = fotopletismografia. Clique aqui para ver uma versão maior desta figura.
Figura 7: Medições da frequência cardíaca. As medidas foram obtidas pelo método proposto em um teste com estímulo provocador. Número da amostra. A caixa abrange os percentis 25e 75, a linha horizontal é a mediana, os bigodes pretos indicam os pontos de dados mais extremos não considerados discrepantes Clique aqui para ver uma versão maior desta figura.
Vídeo Suplementar 1: Ilustração do algoritmo de rede neural proposto para detectar e calcular a frequência cardíaca. O vídeo apresenta larvas de peixe-zebra em movimento livre, um fotopletismograma derivado pelo protocolo proposto e calculado a partir da frequência cardíaca das larvas do fotopletismograma Clique aqui para baixar este arquivo.
Arquivo Suplementar 1: Função MATLAB que executa o aumento da imagem. Clique aqui para baixar este arquivo.
Neste estudo, apresentamos um protocolo experimental para medir os batimentos cardíacos de larvas de peixe-zebra nadando livremente. Avaliamos essa abordagem por meio de vários experimentos, demonstrando sua eficácia. Os principais componentes do método proposto incluem soluções de hardware e software. Primeiramente, utilizamos iluminação infravermelha para a imagem, que, como demonstrado anteriormente, evita problemas relacionados à pigmentação e aumenta a transparência do tecido, facilitando a determinação precisa dos batimentos cardíacos18. Em segundo lugar, propomos imagens dos peixes por baixo. Quando anestesiadas, as larvas do peixe-zebra geralmente ficam de lado devido à posição da bexiga natatória, resultando na diferença de orientação entre as configurações padrão e invertida do microscópio. No entanto, como o coração de um peixe nadador está sempre localizado ventralmente, o posicionamento do canal de imagem abaixo e do canal de iluminação acima aumenta a relação sinal-ruído. Finalmente, melhoramos significativamente o desempenho do processamento empregando uma rede neural para localizar o espécime dentro da imagem e identificar órgãos individuais. Usando conhecimento prévio e dados empíricos, desenvolvemos métodos para identificar o coração no corpo da larva do peixe-zebra, conforme detectado pela rede neural.
Outras vantagens da abordagem incluem o uso de uma fonte de luz LED, que, ao contrário das lâmpadas halógenas, não gera calor. A fixação do canal de iluminação no palco permite o alinhamento em relação ao tamanho do molde uma vez, eliminando a necessidade de ajustes de luz adicionais ao refocar. O algoritmo também pode ser aplicado para gerar grandes conjuntos de dados para o treinamento subsequente de uma rede neural para segmentação cardíaca. Ele acelera o processamento de dados e minimiza erros na determinação da posição do coração. Rotular diretamente o coração nas imagens apresenta um desafio significativo devido ao seu tamanho pequeno e às dificuldades em determinar com precisão seus limites e posição em um único quadro.
Vários fatores críticos devem ser abordados para alcançar os melhores resultados ao implementar este protocolo. Primeiro, o tempo de exposição não deve exceder 1,5 ms, pois tempos de exposição mais longos podem levar ao desfoque da imagem, principalmente durante movimentos rápidos de indivíduos. Esse desfoque dificultaria o cálculo preciso do fotopletismograma em pontos de tempo específicos (Figura 6B). Em segundo lugar, as forças de tensão superficial farão com que a superfície da água na cubeta se curve, levando à dispersão da luz e distorção nas imagens resultantes durante a captura de vídeo. Para mitigar esse efeito, uma camada fina e uniforme de água deve ser aplicada sobre a superfície da agarose para garantir um nível de água plano e consistente em toda a placa de Petri. Como a agarose absorve água gradualmente, recomendamos adicionar água aproximadamente a cada 10 minutos para manter um nível estável. Finalmente, a escolha da ampliação da lente deve ser adaptada ao estágio de desenvolvimento e aos requisitos específicos do experimento. Uma ampliação fixa dentro de uma faixa apropriada é suficiente para medições consistentes dentro da mesma faixa etária e sob condições uniformes. No entanto, uma lente de ampliação variável é mais adequada para estudos que envolvem diferentes estágios ou tarefas de desenvolvimento.
Se o teste revelar que o algoritmo falha, as etapas 2 e 3 devem ser revisitadas e reimplementadas. A etapa 3 geralmente apresenta desafios, principalmente durante a rotulagem de imagens. Garanta uma anotação cuidadosa dos olhos e do corpo, evitando qualquer sobreposição entre essas duas classes. Use imagens que capturem o indivíduo em várias posições, especialmente aquelas que exibem flexão corporal pronunciada, para aumentar a precisão da rotulagem. A etapa 2 deve ser repetida nos casos em que forem identificados problemas com intensidade de iluminação, tempo de exposição ou ampliação. Configurações incorretas podem reduzir a relação sinal-ruído, permitindo que o ruído sobrecarregue o sinal (Figura 6B). Além disso, pode surgir falha no algoritmo se os dados de treinamento e experimentais forem obtidos em diferentes configurações.
A configuração de hardware proposta utiliza uma câmera infravermelha com uma faixa de sensibilidade de 900-1700 nm. No entanto, as câmeras de alta resolução que operam nessa faixa espectral geralmente têm um custo proibitivo. Para resolver essa limitação, câmeras de vídeo industriais de espectro visível com sensibilidade do sensor que se estende até a faixa do infravermelho próximo podem ser empregadas após a remoção de seus filtros infravermelhos como uma alternativa econômica às câmeras infravermelhas. Além do uso de radiação infravermelha, os efeitos da pigmentação também podem ser minimizados com a aplicação de produtos químicos removedores de pigmento ou com o uso de peixes de linhagens genéticas específicas que carecem de pigmentação28,29. A versão atual do software foi projetada exclusivamente para medição da frequência cardíaca. No entanto, outros parâmetros de onda de pulso são igualmente essenciais para pesquisas abrangentes. Futuras atualizações de software se concentrarão na incorporação da análise da variabilidade da frequência cardíaca, um indicador-chave de várias doenças. Medições mais detalhadas, como o intervalo Q-T, são objeto de investigação adicional. Para desenvolver um modelo de rede neural universal capaz de operar com dados de várias idades e condições de aquisição de imagens, o conjunto de dados de treinamento deve incluir diversas amostras com pelo menos 128 imagens rotuladas de cada tipo.
Várias abordagens foram desenvolvidas anteriormente para detecção automatizada da área cardíaca e monitoramento da frequência cardíaca em embriões de peixe-zebra 4,6 e medaka 6,17. A marcação fluorescente do coração em peixe-zebra foi proposta para determinação da área cardíaca30. No entanto, todos os métodos publicados anteriormente limitam-se a trabalhar com embriões imóveis e transparentes durante o breve período pós-fertilização antes do início da motilidade embrionária. Esta é uma limitação significativa que reduz a aplicabilidade dessas técnicas. Outra questão descrita na literatura e em materiais complementares envolve o movimento brusco dos embriões durante o registro da frequência cardíaca 6,17. Esse movimento pode deslocar a área cardíaca alvo do software para o cálculo da frequência cardíaca. A abordagem descrita neste estudo aborda essas deficiências, permitindo o monitoramento de peixes-zebra móveis de até 30 dpf.
As vantagens dessa abordagem expandem significativamente as possibilidades de sua aplicação potencial. Nos últimos anos, o Danio rerio tornou-se um modelo amplamente difundido para o estudo de patologias cardiovasculares e cardiotoxicidade 11,12,31. A capacidade deste método de registrar de forma não invasiva a frequência cardíaca em diferentes estágios iniciais de desenvolvimento sem anestésicos oferece novas oportunidades para estudar a dinâmica das malformações cardíacas induzidas e os efeitos terapêuticos. Atualmente, o monitoramento da frequência cardíaca em embriões de peixe-zebra é usado para triagem de drogas em estudos pré-clínicos 32. As vantagens descritas do método proposto fornecem ferramentas adicionais para avaliar os efeitos retardados e a farmacodinâmica dos medicamentos. Finalmente, a utilização da freqüência cardíaca do peixe-zebra aumenta cada vez mais no monitoramento ecotoxicológico33,34. Neste campo, a abordagem proposta permite a avaliação dos efeitos tóxicos crônicos de poluentes em baixas concentrações durante períodos de exposição prolongados e os efeitos associados à bioacumulação.
Todos os autores revelaram quaisquer conflitos de interesse.
Este estudo foi apoiado pelo Programa de Tarefas do Estado Federal do STC UI RAS (FFNS-2025-0008). Este trabalho foi realizado utilizando o equipamento do Centro de Uso Coletivo do STC UI RAS [http:// https://ckp.ntcup.ru/en/].
Name | Company | Catalog Number | Comments |
Reagents | |||
Low melting agarose | Biozym | 850111 | |
Table salt | Pegasus | N/A | |
Tricaine (Ethyl 3-aminobenzoate methanesulfonate) | Sigma-Aldrich | E10505 | MS-222 |
Equipment | |||
Base with rod | Altami | SM-U1 | |
Collimator lens | JLLSCMGGX | Focal length 30 mm | |
Focusing mechanism | Altami | SM-12 | D=76 mm |
LED | Cree | TR-3535IR-3W | |
Lens | SFK Security | C-Mount, F1.6, 1/3”, | |
Near infrared camera | ToupTek | SWIR1300KMA | |
Pasteur pipette | PE-LD | 149293 | |
Petri Dish 35 x 15 mm | BD Falcon | 351008 | |
Plastic forms | N/A | N/A | Made by 3D printing |
Power supply | Unit-T | UTP3300TFL-II | |
Stage | N/A | N/A | Made by 3D printing |
Stationery knife | ErichKrause | 19145 | |
Test object | Wally Sky | MS-1-EB | |
Software | |||
EfficientDet | N/A | N/A | https://github.com/rwightman/efficientdet-pytorch |
EfficientNet-b0 model | N/A | N/A | https://arxiv.org/abs/1905.11946 |
Google API Client | N/A | Google API Client is a Python client library for Google's discovery-based APIs. https://github.com/googleapis/google-api-python-client | |
Hardware | |||
Multi-scale attention network | N/A | N/A | https://arxiv.org/abs/2209.14145 |
NVIDIA DIGITS | NVIDIA | N/A | NVIDIA DIGITS is a wrapper for Caffe that provides a graphical web interface. https://developer.nvidia.com/digits |
NVIDIA GPU | NVIDIA | N/A | An NVIDIA GPU is needed as some of the software frameworks below will not work otherwise. https://www.nvidia.com |
OpenCV | Intel | N/A | OpenCV is a library for computer vision. https://opencv.org |
Python | Python Software Foundation | N/A | Python is a programming language. https://www.python.org |
Solicitar permissão para reutilizar o texto ou figuras deste artigo JoVE
Solicitar PermissãoThis article has been published
Video Coming Soon
Copyright © 2025 MyJoVE Corporation. Todos os direitos reservados