1. Raccolta dati di mappatura ottica Eseguire la mappatura ottica cardiaca utilizzando uno di una vasta gamma di modelli sperimentali tra cui cuori interi intatti e isolati6,18, atri isolati14,19, cunei ventricolari20, fette cardiache 21 Mieto , 22e monostrati cellulari23. Vedere i riferimenti associati per i progetti sperimentali per raccogliere dati di mappatura ottica non elaborati da questi preparati. A condizione che i dati ottenuti possano essere convertiti in una pila tiff o salvati in un file . MAT file, dovrebbe essere analizzabile utilizzando ElectroMap. Sono inclusi i dati di dimensioni variabili (quadrato/rettangolare) e le risoluzioni (massimo testato attualmente 2048 pixel x 2048 pixel). 2. Installazione e avvio del software NOTA: di seguito sono descritti in dettaglio i due metodi per l’installazione e l’esecuzione di ElectroMap, ovvero all’interno di MATLAB eseguito dal codice di origine (m) o come file eseguibile autonomo (.exe for windows). Il software finale e la sua funzionalità sono invarianti tra le due opzioni di installazione (diverse da alcune differenze nella navigazione nella directory). Pertanto, le considerazioni principali per la scelta della versione da installare sono l’accesso a MATLAB e alle caselle degli strumenti necessarie e se si desidera l’accesso al codice sorgente. Ove possibile, si consiglia di utilizzare la versione MATLAB per tempi di avvio più rapidi, tempi di elaborazione più brevi e una segnalazione degli errori più semplice. Configurazione 1: elettromappa in esecuzione all’interno di MATLAB Installare MATLAB. ElectroMap è stato progettato in MATLAB 2017a, tuttavia, il software è stato testato per l’uso in tutte le versioni successive di MATLAB (fino a 2018b al momento della scrittura). Sono necessarie le seguenti caselle degli strumenti: Elaborazione immagine, Elaborazione del segnale, Statistiche e Apprendimento automatico e Adattamento curva. Scaricare/clonare tutti i file dall’ultima versione ‘codice sorgente’ di ElectroMap dal repository GitHub (https://github.com/CXO531/ElectroMap). Decomprimere il contenuto scaricato nella posizione desiderata. Aprire MATLAB e passare al percorso della cartella che ospita il codice sorgente ElectroMap. Aprire quindi il file ElectroMap.m e premere run nell’editor oppure digitare ElectroMap nella finestra di comando e premere INVIO. Verrà avviata l’interfaccia utente ElectroMap, Figura 1A. Installazione 2: file .exe autonomo Scaricare il file del programma di installazione: https://drive.google.com/open?id=1nJyI07w9WIt5zWcit0aEyIbtg31tANxI. Seguire le istruzioni del programma di installazione, che scaricherà richiesto il runtime MATLAB dal web insieme al software ElectroMap. Eseguire ElectroMap.exe.NOTA: il tempo di avvio per la versione standalone può essere di diversi minuti. 3. Caricamento e pre-elaborazione delle immagini Premere Seleziona cartella e passare alla posizione dei file di dati da analizzare. In questo modo la casella di riepilogo a sinistra verrà popolata con tutti i file all’interno di tale directory che sono del tipo di file corretto (.tif o . MAT). . I file MAT devono contenere solo la variabile dello stack di immagini.NOTA: durante lo spostamento nel selettore di directory verranno visualizzate solo le cartelle e non i singoli file. Selezionare il file da caricare dall’interno dell’interfaccia e premere Carica immagini. Una volta caricato, apparirà il primo fotogramma e il contorno rosso indicherà la soglia automatica dell’immagine. Se necessario, ricaricare i ROI utilizzati in precedenza selezionando Salva/Carica ROI. In questo caso, saltare il passaggio 3.3. Per impostazione predefinita, la soglia si basa sulle intensità dei pixel nel primo fotogramma. Se lo si desidera, modificarla a una soglia in base all’ampiezza del percorso del tempo del segnale modificando l’opzione nel menu a discesa Immagine per soglia. Si prega di notare che una volta selezionata la soglia, viene quindi applicato per l’intero stack di immagini. Se lo si desidera, modificare l’opzione di soglia su manuale, che attiverà il dispositivo di scorrimento per regolare manualmente la soglia dell’immagine. Inoltre, ritagliare le immagini (Ritaglia immagine) e/o disegnare un’area di interesse personalizzata (ROI personalizzato) per l’analisi selezionando le caselle di controllo appropriate sotto le opzioni di soglia. Si noti che le opzioni avanzate per la selezione dell’area di interesse, ad esempio il numero di aree, sono disponibili da Selezione ROI nel menu in alto. Una volta applicata una soglia appropriata, premere Elabora immagini per applicare l’elaborazione. Le impostazioni per l’elaborazione sono descritte in dettaglio di seguito (passaggio 3.4.1-3.4.5). A questo punto, assicurarsi che siano state immesse le impostazioni corrette della fotocamera. Queste sono le dimensioni dei pixel in : (IMPORTANTE: questa è la dimensione in pixel dell’immagine, e non le dimensioni dei pixel che compongono il chip o l’hardware equivalente nel dispositivo di imaging) e Framerate in kHz. Per l’inversione del segnale, selezionare la casella di controllo Inverti dati per abilitare. Se segnalato segnale fluorescente è inversamente proporzionale al parametro di interesse (come con i coloranti potentiometrici comunemente utilizzati) il segnale può essere invertito. Per il filtro spaziale, selezionare Gaussiano o Media dal menu del kernel. La dimensione dell’area mediata spazialmente è controllata dall’ingresso Dimensione adiacente al menu a discesa Kernel (cioè 3 genera un kernel di filtro 3 x 3 pixel). Quando si applica un filtro gaussiano, la deviazione standard può essere impostata anche dall’ingresso Sigma. Per la correzione della linea di base, selezionare Correzione a cappello24 o polinomiale (4o 11 gradi)dal menu Linea di base. La correzione può essere applicata a ogni pixel singolarmente (tempo di elaborazione prolungato) o come media dell’intera immagine (più veloce ma presuppone alterazioni della linea di base omogenee). La correzione del cappello superiore può essere modificata anche impostando la lunghezza del cappello in millisecondi, adiacente al menu a discesa di selezione della linea di base. La lunghezza del kernel Top-Hat deve essere maggiore della scala cronologica dei singoli potenziali di azione/transitori di calcio. Per il filtro temporale, selezionare Filtro Savitzky-Goaly o IIR (Infinite Impulse) dal menu Filtro.NOTA: Oltre al segnale media del tessuto che appare in basso a sinistra, il filtraggio temporale viene applicato a ciascun pixel singolarmente al momento della quantificazione dei parametri dagli intervalli di immagini medi dell’ensemble. Questo è stato implementato per ridurre i tempi di elaborazione filtrando piccole sezioni di dati quando necessario, piuttosto che interi file. Per la rimozione dei fotogrammi, si noti che se è selezionata l’opzione Rimuovi fotogrammi, è possibile rimuovere grandi picchi con ampiezza maggiore del segnale di interesse dal set di immagini. Questo può essere utile in set di dati a ritmo ottico come il ritmo optogenetico in cui la depolarizzazione viene avviata dall’attivazione ottica di operazioni come la channelrhodopsin 211,12.NOTA: poiché la rimozione dei fotogrammi potrebbe introdurre modifiche non fisiologiche ai passaggi nei segnali dell’immagine, il filtraggio temporale può introdurre artefatti nei dati e quindi non è raccomandato qui. Si noti che il segnale verrà segmentato una volta che l’opzione Elabora immagini è stata selezionata in base alle opzioni di segmentazione,tuttavia questo può essere rapidamente modificato senza rielaborare l’intero set di dati (vedere la sezione 4). 4. Segmentazione dei dati e media dell’insieme NOTA: Una volta che il file è stato elaborato, i picchi nel segnale media dei tessuti (traccia in basso a destra, Figura 1A) saranno stati rilevati ed etichettati da cerchi rossi. Vengono conteggiati solo i picchi superiori a una soglia impostata (linea blu sulla traccia impostata dalla sogliadi picco). Inoltre, i picchi vengono conteggiati solo se sono sufficientemente ritardati rispetto ai picchi precedenti, impostati dall’ingresso Distanza di picco min. Il segnale viene quindi segmentato in base ai picchi rilevati. In primo luogo, la lunghezza effettiva del ciclo (CL) di ogni picco viene calcolata misurando il tempo tra esso e il picco successivo. Se un numero di picchi (impostato dall’input Numero minimo di picchi) ha CL simili (soglia per la quale è impostato da Input limite minimo), vengono raggruppati e il CL medio per tali picchi calcolati. Per un’ulteriore segmentazione dei dati, premere Segment Signal. Le opzioni di sottosegmentazione sono: Nessuno – tutti i picchi con lo stesso CL raggruppati insieme; Tutto: vengono identificati segmenti di npicchi entro i tempi CL costanti (npicchi sono impostati dall’input Dimensione segmento); Ultimo – Ipicchi finali n prima di una modifica CL vengono identificati e raggruppati e tutti gli altri non vengono analizzati; e Single Beat – Equivale ad applicare la segmentazione Tutti con npicchi , 1 e quindi non viene applicato alcun raggruppamento o media dell’insieme (vedere 4.5). Questo può essere applicato selezionando il pulsante Battuta singola. Applicare la segmentazione personalizzata del segnale ingrandendo un momento di interesse e selezionando Segnale di segmento. Verrà aggiunta un’opzione aggiuntiva intitolata Sezione ingrandita alla casella di riepilogo della sezione, corrispondente ai punti temporali selezionati. I risultati della segmentazione appariranno nella casella di riepilogo adiacente al segnale mediato del tessuto e mostreranno il numero di sezione e il CL stimato. Tutte le sezioni temporali segmentate sono indicate da colori diversi. Selezionare un segmento dalla casella di riepilogo per evidenziare la sezione in rosso. Questo attiverà automaticamente anche le analisi di questa sezione, come se fosse stato selezionato il pulsante Produce Maps (vedi sezione 5). Le analisi dei picchi raggruppati verranno eseguite sui dati “in media dell’insieme”. Ciò comporta la media dei picchi in un segmento insieme, con i tempi di riferimento i picchi identificati nel passaggio 4.2. Aggiornare l’intervallo di tempo per la media modificando gli ingressi prima e dopo e premendo Segment Signal. 5. Potenziale d’azione/durata transitoria del calcio e analisi della velocità di conduzione Una volta elaborate le immagini, il pulsante Produce Maps diventerà attivo. Premere Produce Maps per applicare la durata del potenziale di azione (APD), il tempo di attivazione, la velocità di conduzione e l’analisi SNR. Per impostazione predefinita, l’analisi verrà applicata al primo segmento di segnale. Seleziona altri segmenti dalla casella di riepilogo applicherà l’analisi al segmento selezionato.NOTA: i risultati dell’analisi vengono visualizzati nella tabella dei risultati, tra cui la media, la deviazione standard, l’errore standard, la varianza e l’analisi dal quinto al 95o percentile. Le mappe di durata sono tuttavia definite mappe “APD”, i segnali di calcio elaborati utilizzando le stesse impostazioni misureranno la durata transitoria del calcio. Selezionare Ottieni informazioni pixel per visualizzare una visualizzazione dettagliata del segnale da qualsiasi pixel all’interno dell’immagine e Confronta pixel per tracciare contemporaneamente i segnali da un massimo di 6 posizioni. Utilizzare il pannello Elaborazione segnale per regolare le impostazioni per l’analisi della durata. Questi sono: Durata – Tempo di percentuale di ripolarizzazione/decadimento per misurare dal picco; Base ‘APD’: periodo di tempo del segnale definito come linea di base di riferimento per le misurazioni di ampiezza; e ‘APD’ Ora di inizio– Ora di inizio per le misurazioni della durata. Queste sono le stesse opzioni per decidere il tempo di attivazione per le mappe isocroniche (discusse di seguito) e sono definite: Start ( d2F/dt2max), Upstroke (dF/dtmax), Depolarizzazione del punto medio (tempo di ampiezza del 50%), Picco (tempo di massima ampiezza). Queste definizioni applicate ai potenziali di azione di topi e cavie sono illustrate nella Figura 2A.NOTA: se si modifica una di queste opzioni, la mappa della durata e la tabella dei risultati verranno aggiornate automaticamente. Sono disponibili anche le opzioni di rimozione della scala della mappa e dell’outlier. Anche la velocità di conduzione viene misurata automaticamente all’interno dell’interfaccia software principale. Ciò si ottiene utilizzando il metodo multi-vettore di Bayly et al26 dalla mappa isoarta definita dalla misura di attivazione scelta (discussa nel passaggio 5.4). Premere Punti di attivazione per eseguire il rendering di una rappresentazione 3D della mappa di attivazione. Il metodo di misurazione della velocità di conduzione multivettore segmenta spazialmente la mappa isocronica in regioni di n x n pixel. Impostare il valore di n utilizzando l’input Dimensione finestra locale e impostare l’intervallo di tempi di attivazione per applicare l’analisi all’utilizzo degli input Tempi di attivazione adattamento.NOTA: per ogni regione locale è disponibile una superficie polinomiale, f, che meglio descrive la relazione tra tempo di attivazione e posizione spaziale (x,y). Il vettore di gradiente, CVlocal, di questa superficie viene quindi calcolato come:①dove indica l’operatore differenziale spaziale cartesiano bidimensionale26. Per ogni pixel nella mappa isocronica, viene calcolato un vettore locale che rappresenta la velocità e la direzione di conduzione. Selezionare Mappa isocronologia con vettori dal menu a discesa di visualizzazione per visualizzare questa analisi. SNR viene calcolato come rapporto tra l’ampiezza massima rispetto alla deviazione standard del segnale al basale. Questa analisi viene eseguita dopo tutti i passaggi di elaborazione. Premere il calcolo SNR nel menu in alto per modificare le impostazioni per il periodo del segnale definito come linea di base. 6. Modulo di analisi della conduzione Premere Conduzione per accedere ad un’analisi più dettagliata della velocità di conduzione. Si apre un modulo separato in cui la conduzione può essere quantificata utilizzando il metodo Multi-vector Bayly come nell’interfaccia principale, metodi a vettore singolo e come curva di attivazione. Premere Single Vector per analizzare la conduzione utilizzando il metodo a singolo vettore, in cui CV viene calcolato dal ritardo nel tempo di attivazione tra due punti. Questa operazione può essere eseguita utilizzando i metodi Automatic o Manual, selezionabili sotto il pulsante Vettore singolo. Per il metodo a vettore singolo automatico, selezionare una distanza e un punto iniziale da cui misurare la conduzione. Il software eseguirà quindi una sweep di 360 gradi dal punto selezionato, misurando il ritardo di tempo e calcolando la velocità di conduzione associata lungo tutte le direzioni in incrementi di 1 grado. I risultati di questa analisi vengono visualizzati nel grafico adiacente alla mappa e la direzione di conduzione più lenta viene visualizzata in rosso. Per il metodo a vettore singolo manuale, scegliere un punto iniziale e un punto finale dalla mappa isocronologica per calcolare la velocità di conduzione. Per selezionare un nuovo punto iniziale, premere Cancella punto iniziale. Premere Vettore locale per applicare il metodo multivettore, con le impostazioni corrispondenti a quelle dell’interfaccia principale. All’interno del modulo di conduzione, è possibile visualizzare la distribuzione delle velocità di conduzione, nonché la distribuzione angolare dei vettori calcolati e la dipendenza angolare della velocità di conduzione. Premere la curva di attivazione per tracciare la percentuale di tessuto attivata in funzione del tempo. Viene visualizzata automaticamente l’attivazione del tempo al 100%, mentre è possibile selezionare automaticamente i valori personalizzati per le percentuali di attivazione minime (blu) e massime (rosse) tra cui misurare. 7. Ulteriori analisi e moduli Oltre alle analisi della durata e della velocità di conduzione eseguite automaticamente, è possibile quantificare diversi altri parametri utilizzando ElectroMap. Queste analisi sono selezionabili dal menu a discesa sopra la mappa di visualizzazione. Selezionare una di queste opzioni per eseguire l’analisi e i risultati verranno visualizzati nella riga 4 della tabella dei risultati: 1) Intervallo diastolico – Tempo dal 90% di ripolarizzazione al tempo di attivazione del potenziale di azione successivo; 2) Frequenza dominante – Lo spettro di frequenza di ogni pixel viene calcolato utilizzando la trasformazione veloce di Fourier, e la frequenza con la maggior potenza è definita come la frequenza dominante. Le impostazioni avanzate di intervallo e finestra per l’analisi della frequenza dominante sono disponibili selezionando Mappatura frequenza; 3) Tempo di picco – Il tempo di aumento tra due percentuali selezionate dall’utente (predefinito da 10 a 90%) della fase di depolarizzazione del potenziale d’azione o del rilascio di calcio. I valori percentuali possono essere modificati selezionando Impostazioni TTP; e 4) Costante di rilassamento ()– La costante di rilassamento viene calcolata montando un decadimento mono-esponenziale della forma della forma:②dove il livello di fluorescenza al momento t dipende dalla fluorescenza di picco, F0e dal decadimento successivo (C è una costante)27. Il valore tra il quale adattare l’equazione 2 sono selezionabili all’interno delle interfacce utente principali di ElectroMap, nonché una bontà dei criteri di esclusione di adattamento in base al valore r 2. Premere Single File Analysis per aprire un modulo dedicato per la durata ad alta velocità effettiva e l’analisi della conduzione di ogni segmento identificato in un file. L’analisi può essere eseguita sull’intera immagine (durata, conduzione e tempo di attivazione) o su regioni selezionate o punti di interesse (attualmente solo durata). I risultati vengono restituiti in un file CSV.NOTA: per i valori APD dell’intera immagine, la prima colonna nel file CSV è la media, mentre la seconda colonna è la deviazione standard. Premere Alternans per avviare un modulo autonomo per l’analisi dedicata e la mappatura della variabilità beat-to-beat. Vedere O’Shea et al. 201913 per informazioni dettagliate sulle opzioni di elaborazione e analisi delle opzioni di elaborazione e analisi degli alternani. In particolare, questo modulo è progettato per identificare due oscillazioni del periodo, note come alternani. Sia gli alternari di durata che di ampiezza vengono calcolati e restituiti.NOTA: le alternanza di durata vengono misurate confrontando la misurazione della durata da un picco al successivo; cioè se il picco uno e due e APD1 e APD2 rispettivamente, allora la durata alternana③La misurazione della durata viene eseguita utilizzando le impostazioni nell’interfaccia principale. Nel frattempo, gli alternanzi di ampiezza possono essere quantificati e mappati tra preparazioni multicellulari come cambiamento assoluto (definito come percentuale in cui 0% – stessa ampiezza tra un battito e l’altro). Inoltre, gli effetti di fenomeni come il carico di calcio possono essere ulteriormente studiati misurando e confrontando le alternanarie di carico e rilascio, come è stato precedentemente riportato28. Se L è definito come l’ampiezza di picco dei grandi battiti (cioè dove l’ampiezza è maggiore del battitoprecedente), S l’ampiezza dei piccoli battiti eD il carico diastolico dei piccoli battiti, il rilascio alternans ( ) sono definiti come:④Al contrario, lealternanni di carico ( ) sono definite come:⑤Le misurazioni alternativi possono essere effettuate su tutto il tessuto e i risultati dell’analisi vengono visualizzati in basso a destra del modulo. Quando si utilizza il modulo per la prima volta, l’analisi viene eseguita su tutto il file sperimentale e i risultati visualizzati sono una differenza di beat-beat media nell’intero file. Tuttavia, l’analisi può essere limitata a orari specifici del file deselezionando Tieni premuto lo zoom, ingrandendo un periodo di tempo specifico e selezionando Analizza sezione ingrandita. Questo aggiornerà il pannello dei risultati per mostrare l’analisi del periodo di tempo selezionato. Selezionare play per visualizzare un video beat-to-beat dell’analisi alternano. Inoltre, selezionare Crea mappa media per esportare una mappa del comportamento alternano mediato dai punti temporali selezionati, impostati nel menu a comparsa quando si utilizza questa funzione. Premere Fase Mappa per avviare il modulo di mappatura di fase. Una trasformazione di Hilbert viene eseguita per calcolare la fase istantanea (tra -e ) dei segnali in ogni punto temporale. Premere il dispositivo di scorrimento riproduzione o trascinamento per visualizzare il comportamento della fase nel tempo e fare clic su un pixel per eseguire il rendering di un diagramma di fase. 8. Esportazione dei dati I dati vengono esportati da ElectroMap in una varietà di forme. Premere Esporta valori per salvare i valori della mappa attualmente visualizzata nell’interfaccia principale utilizzata. I valori misurati possono essere salvati come mappa (preservando le posizioni dei pixel) o condensati in un unico elenco e possono essere salvati come .csv,.txt o . MAT. Premere Esporta mappa per visualizzare un popup contenente la mappa attualmente visualizzata, che può quindi essere salvata in una varietà di formati di immagine. Le opzioni di visualizzazione per la mappa vengono controllate selezionando Impostazioni mappa, ma possono essere modificate anche dopo aver selezionato Esporta mappa. Ad esempio, è possibile aggiungere una barra dei colori selezionando questa icona dal menu in alto e la scala può essere impostata selezionando Modifica > Mappa colori. Premere Video di attivazione per eseguire il rendering di un’animazione della sequenza di attivazione, che può essere salvata come file .gif animato. Premere SegmentO Video per salvare un file video .avi del parametro attualmente visualizzato di ogni segmento identificato.