Method Article
* Questi autori hanno contribuito in egual misura
L'analisi citogenetica del cromosoma dicentrico (DC) quantifica l'esposizione a radiazioni ionizzanti. Il software automatizzato di identificatore di cromosoma dicentrico e l'estimatore di Dose con precisione e rapidamente stima dose biologica da DCs in cellule di metafase. Esso distingue monocentric cromosomi e altri oggetti da DCs e le stime di dose di radiazione biologica dalla frequenza di DCs.
Dose di radiazione biologica può essere stimato da frequenze di cromosoma dicentrico in cellule di metafase. L'esecuzione di queste analisi citogenetica del cromosoma dicentrico è tradizionalmente un processo manuale e laborioso, non adatto per gestire il volume di campioni che possono richiedere l'esame a seguito di un evento di massa dell'incidente. Software automatizzato di identificatore di cromosoma dicentrico e Dose Estimator (ADCI) automatizza questo processo esaminando insiemi di immagini di metafase utilizzando tecniche di elaborazione immagine basata su apprendimento di macchina. Il Seleziona software immagini appropriate per analisi rimuovendo immagini inadatti, classifica ogni oggetto come un centromero-contenenti del cromosoma o non-cromosoma, più ulteriormente contraddistingue i cromosomi come monocentric cromosomi (MCs) o dicentrico cromosomi (DCs), determina la frequenza DC all'interno di un campione e le stime di dose di radiazione biologica confrontando la frequenza di campionamento DC con curve di calibrazione calcolate utilizzando campioni di calibrazione. Questo protocollo viene descritto l'utilizzo del software ADCI. In genere, sia calibrazione (dose nota) e insiemi di test (sconosciuto dose) di immagini di metafase sono importati per eseguire la stima accurata dose. Immagini ottimali per l'analisi possono essere trovati automaticamente utilizzando filtri di immagine preimpostati o possono anche essere filtrati attraverso ispezione manuale. Il software elabora le immagini all'interno di ogni campione e frequenze DC vengono calcolate a diversi livelli di rigore per la chiamata a controller di dominio, utilizzando un'approccio di apprendimento della macchina. Curve di taratura lineare-quadratica vengono generate in base sulle frequenze di DC in campioni di calibrazione esposti a dosi fisiche conosciute. Dosi di campioni esposti a livelli di radiazioni incerto sono stimate da loro frequenze DC utilizzando queste curve di calibrazione. Rapporti possono essere generati su richiesta e forniscono la sintesi dei risultati di uno o più campioni, di una o più curve di calibrazione o di stima di dose.
Radiazione biodosimetry utilizza marcatori biologici, principalmente le aberrazioni cromosomiche quali cromosomi dicentrici (DCs) e traslocazioni cromosomiche per misurare dosi di radiazioni che gli individui sono esposti a. Una dose biologicamente assorbita potrebbe essere diversa dalla dose fisica misurata dagli strumenti a causa della variabilità tra individui. Allo stesso modo, radiazione di una certa dose di fisica può produrre diverse esposizioni biologiche dovute a sottostanti condizioni fisiologiche o ambientali. Conoscenza della dose biologica è di particolare importanza per la diagnosi ed il trattamento.
Il dosaggio di DC è il gold standard della organizzazione mondiale della sanità (OMS) e Agenzia internazionale per energia atomica (AIEA) per valutare l'esposizione a radiazioni biologico nelle persone. E ' stato il primo test consigliato dall'AIEA e che per la valutazione della dose di radiazione. Frequenza DC è relativamente stabile per circa 4 settimane dopo l'esposizione di radiazioni1 e loro correlazione quantitativa con la dose di radiazione emessa è accurata, che rendono DCs il biomarcatore ideale. Il rapporto tra la dose di radiazione (a cui fa riferimento in unità Gray [Gy]) e la frequenza di DC (indicato come numero di DCs per cella) può essere espressa come una funzione lineare-quadratica.
L'analisi citogenetica di DC è stato lo standard industriale per circa 55 anni2. È stata eseguita manualmente, che richiedono 1-2 giorni per analizzare i dati di microscopio da un campione di sangue. Diverse centinaia a diverse migliaia di immagini sono necessarie per stimare con precisione l'esposizione a radiazioni a seconda della dose3. A dosi superiori a 1 Gy, AIEA raccomanda un minimo di 100 DCs essere rilevato. Esame di 250-500 immagini di metafase è pratica comune nei laboratori di citogenetica biodosimetry. Per campioni con esposizioni < 1 Gy, 3.000-5.000 immagini sono suggeriti dovuto la probabilità inferiore di formazione DC. In entrambi i casi, è un'attività intensa di lavoro.
Laboratori di citogenetica biodosimetry creano proprio in vitro le curve di calibrazione biodosimetry di radiazione prima di valutare le dosi biologiche in campioni di prova. I campioni di sangue da individui di controllo normali, sono esposti a radiazioni e linfociti sono quindi colta e preparati per l'analisi del cromosoma di metafase. Utilizzando questi campioni, hanno ricevute le dosi biologiche sono tarate per il note fisiche dosi emesse da una sorgente di radiazione standard. Dopo immagini delle cellule di metafase sono stati registrati, esperti esaminare immagini, contano DCs e calcolare le frequenze DC per ogni campione. Una curva di calibrazione è costruita montando una curva lineare-quadratica per le frequenze di DC a tutte le dosi. Quindi, l'esposizione in esempio di test da individui possono essere dedotte abbinando le frequenze DC per le dosi calibrate sulla curva o specificandoli nella formula quadratica lineare corrispondente.
Ci hanno automatizzato sia il rilevamento dei controller di dominio e determinazione per accelerare questa procedura utilizzando il software della dose. Identificatore di cromosoma dicentrico automatizzati e Dose Estimator (ADCI) utilizza immagine di macchina basata su apprendimento tecniche di lavorazione per rilevare e discriminare i cromosomi dicentrici (DCs) da monocentric cromosomi (MCs) e altri oggetti e automatizza radiazione stima della dose. Il software mira a significativamente ridurre o eliminare la necessità per la verifica manuale dei conteggi di DC e ad accelerare la stima della dose tramite l'automazione. Esso è stato sviluppato con il coinvolgimento dei laboratori di riferimento biodosimetry a salute Canada (HC) e laboratori nucleari canadesi (CNL). Loro feedback farà in modo che le prestazioni continuerà a soddisfare criteri di AIEA per questo test.
Il software svolge le seguenti funzioni: 1) filtro DCs e selezionando immagini cellulare ottimale metafase per analisi, riconoscimento 2) cromosoma, rilevamento di DC e DC frequenza determinazione e 3) stimare la dose di radiazioni da dose-calibrato, dati relativi alla radiazione citogenetica. Questo software elabora gruppi di immagini di metafase da stesso individuo (definiti un campione), conta il numero di controller di dominio in ogni utilizzo immagine tecniche di lavorazione e restituisce la dose di radiazione stimato ricevuta da ciascun campione in unità di grigi (Gy).
Il software è stato progettato per gestire una gamma di strutture cromosomiche, conteggi e densità. Tuttavia, l'algoritmo esegue in modo ottimale nelle immagini di metafase che contiene un complemento completo nei pressi di cromosomi ben separati, lineare4. Immagini che contengono altamente sovrapposti set di cromosomi, più cellule, cellule in metafase incompleta, sorella cromatidio della separazione, i nuclei, oggetti non-cromosomici e altri difetti possono ridurre la precisione dell'algoritmo. Modelli di selezione immagine e altri segmentazione oggetto soglie possono filtrare la maggior parte delle immagini sub-ottimale e falsi positivi DCs dedicati.
Cromosoma dicentrico rilevamento viene eseguito quando viene elaborata un'immagine. L'algoritmo tenta di determinare quali oggetti in un'immagine sono i cromosomi e quindi individua le due regioni più probabile di essere centromeri su ogni cromosoma. Quindi, una serie di supporto differenti Vector Machine (SVM) modelli di apprendimento distinguere cromosomi come DCs o normale, monocentric cromosomi. I modelli SVM differiscono nella sensibilità e specificità di rilevazione DC (Vedi punto 3.1.4 qui sotto), che può interessare le frequenze di DC che sono determinate in un campione.
ADCI elabora set di Giemsa - (o DAPI-) tinto immagini digitali di metafase (in formato TIFF o JPG) per uno o più campioni. Il software analizza DCs in entrambi esempi di taratura e prova. Le dosi di fisiche (in Gy) dei campioni di calibrazione sono noti e sono utilizzate nella generazione di una curva di calibrazione. Le dosi fisiche e biologiche degli individui con le esposizioni sconosciute vengono dedotti dal software dalla curva di calibrazione generato dal computer. Anche se laboratori utilizzano tecniche comparabili, le curve di calibrazione da diversi laboratori variano spesso3. Entrambi esempi di curva e test di calibrazione dal laboratorio stesso devono essere elaborati per la stima della dose esatta in campioni di prova.
Questo software offre velocità, precisione e scalabilità quali indirizzi la produttività necessaria per gestire un evento in cui molti individui contemporaneamente devono essere testati. E ' stato sviluppato dal 2008-20174,5,6,7,8,9,10,11,12 ,13. Utilizzando hardware computer recente, questo desktopSoftware per PC in grado di elaborare e dose di radiazioni di stima in un campione di 500 equivalenti di genoma di metafase in 10-20 min 4. Il codice è basato su un insieme di segmentazione di immagini proprietarie e machine learning algoritmi per l'analisi del cromosoma. Esperto di analisi di ogni cromosoma 3 Gy radiazioni ha dato esattezze paragonabili a ADCI. In una serie di 6 campioni di esposizioni sconosciute (precedentemente utilizzate in un esercizio di competenza internazionale), il software stimato dosi entro 0,5 Gy dei valori ottenuti da revisione manuale dei dati stessi di HC e CNL, che soddisfi i requisiti dell'AIEA per triage biodosimetry. Inoltre, standardizzazione interlaboratorio e, infine, la riproducibilità della dose stima beneficio dall'avere un controller di dominio comune, automatizzato, segnando algoritmo. Tuttavia, il software permette personalizzazione dell'immagine di filtraggio e criteri di selezione, che consente differenze nei metodi di preparazione del cromosoma e sorgenti di taratura di radiazione per essere presi in considerazione.
Questo software è un'interfaccia di utente grafica (GUI) - basati su sistema che analizza set di immagini di cromosoma contenente Giemsa (o DAPI) - macchiato cellule in metafase per le anomalie che derivano dall'esposizione a radiazioni ionizzanti. I set di immagine sono fotografati digitalmente con un sistema di microscopio luce (o epifluorescente) e ogni insieme corrisponde a un campione diverso. Il software utilizza tecniche per rilevare e discriminare le DCs da MCs e altri oggetti di elaborazione delle immagini. Empiricamente derivati segmentazione filtri eliminano quindi automaticamente i falsi positivi DCs senza intaccare il vero DCs. Infine, il software automaticamente filtra indesiderati immagini basati su varie strutture di immagine trovate immagini di metafase di scarsa qualità con modelli di selezione immagine pre-calcolate (o specificato dall'utente). Queste immagini includono quelli contenenti eccessiva o insufficiente numero di oggetti "rumorosi", più cromosomi sovrapposti, immagini mancanti cromosomi in metafase, un numero eccessivo di sorella cromatidi4. I dati di immagine automaticamente a cura vengono utilizzati per generare la curva di calibrazione di dose da campioni della dose di radiazione noti e sono utilizzati per stimare le esposizioni di prova campioni esposti a dosi sconosciute.
Output del software possono essere visualizzati e salvati come: 1) basato su testo output visualizzato nella console, 2) trame che possono essere salvati come immagini e 3) di report in formato HTML. Molti aspetti del software sono personalizzabili per soddisfare le esigenze specifiche dei diversi laboratori. Laboratori individuali forniscono solitamente i campioni di test e calibrazione preparato e raccolti basata sul protocollo citogenetico convalidato in quel laboratorio. Questo mantiene l'uniformità della preparazione del campione e permette curve di calibrazione generate dai campioni di calibrazione chiaramente da applicare per testare campioni derivati usando lo stesso protocollo. Curve di calibrazione possono essere create anche da coefficienti di curva o frequenze DC alle dosi definite. Le stime più accurate di dose sono ottenute filtrando immagini di qualità inferiore e falsi positivi DCs (FPs). Selezione dei sottoinsiemi di immagine ottimale all'interno di ogni campione viene eseguita utilizzando 'Modelli di immagine selezione' che eliminano immagini scadente che tendono a introdurre FPs. Una serie di modelli pre-convalidate sono inclusi con il software, tuttavia ulteriori modelli con soglie personalizzate e filtri possono essere creati e salvati dall'utente.
Una volta che il software viene caricato correttamente, la principale interfaccia grafica (GUI) è presentato (Vedi Figura 1). Da questa interfaccia, campioni, ciascuno composto da una cartella di file immagine di metafase cella, possono essere selezionate e lavorate per identificare controller di dominio, curve di calibrazione possono essere create e confrontate e dose di esposizione di radiazioni di campioni può essere determinata.
Figura 1: The settori principali dell'interfaccia di utente grafica includono: un elenco di campioni (1), un elenco di calibrazione delle curve (2), il processo di coda (3), che controlla lo stato del rilevamento di DC in ogni set di immagini di ogni campione, una trama display (4), che riassume le statistiche o altre proprietà quantitativa di un set di immagini nei campioni o curve di calibrazione e una console (5) che contiene un testo descrittivo come uscite di ogni operazione eseguita dal programma. Clicca qui per visualizzare una versione più grande di questa figura.
1. importazione e campioni di processo
2. visualizzazione e selezione di immagini (facoltativo, passaggio consigliato)
Nota: questo passaggio viene descritto l'utilizzo del Visualizzatore di immagine di metafase e la creazione di un modello di selezione immagine. Alcuni modelli di selezione immagine convalidati sono inclusi con il software che può essere utilizzato nella generazione della curva di calibrazione e stima della dose. Così, questo passaggio non è obbligatorio, tuttavia può essere utilizzato come una guida che descrive i passaggi necessari per eseguire questa operazione se si desidera.
3. Generazione della curva di
Figura 2: visualizzazione dell'effetto di modifica del valore di Sigma SVM dall'algoritmo sul vero positivo (TP) e falsi positivi (FP) DC conteggi, il Calue predittivo positivo (PPV) e vero positivo tassi (TPR). Clicca qui per visualizzare una versione più grande di questa figura.
4. Stima della dose
5. Segnalazione
Nota: il metodo utilizzato per denominare un report e selezionare una directory all'interno del quale viene salvato è comune a tutti i tipi di report. A ' nome Report ' devono essere fornite. Quando viene generato un report, verrà creata una directory contenente i file di report utilizzando questo nome automaticamente. Questa directory verrà inserita all'interno è il ' cartella Report '. Per impostazione predefinita, il ' cartella Report ' è una directory denominata ' report ' trovato nella directory dei dati specificata durante l'installazione.
6. Funzionalità di controllo
Nota: il software registra tutte le operazioni eseguite durante una sessione in un file di log. Il programma fornisce un'applicazione software accessorio che consente i file di registro essere guardato, cercato, utilizzata per valutare l'integrità di un'analisi e in alcuni casi, di recuperare i dati di esempio da incompleta o prematuramente terminata sessioni.
7. Curva e Dose stima statistiche opzioni
Test del software è stato effettuato con i dati immagine del cromosoma di metafase ottenuti da HC e CNL. I campioni di sangue sono stati irradiati da un'unità di XRAD-320 (250 kV i raggi x, 12,5 mA, 2mm Al filtraggio, rateo di dose: 0.92 o 1,7 Gy/min) calibrato con una camera a ioni a HC e lavorate ad entrambi i laboratori. Campioni di sangue periferico del linfocita erano coltivati, corretti e macchiati presso ogni impianto secondo protocolli stabiliti3,15. Immagini di metafase da vetrini colorati Giemsa furono catturati in modo indipendente da ogni laboratorio utilizzando un sistema di microscopia automatizzata. Esperti in ciascun laboratorio segnato DCs in molti di questi campioni manualmente, costruito le proprie curve di calibrazione e stimato dosi di campioni di prova di esposizioni sconosciute. Nella tabella 1 viene fornita una descrizione dettagliata di questi set di dati.
Dose fisica | Scopo | Preparazione di HC | Preparazione di CNL | ||
Cui nome | n # di immagini | Cui nome | n # di immagini | ||
0 Gy | Calibrazione | HC0Gy | 731 | CNL0Gy | 798 |
0.1 Gy | Calibrazione | HC01Gy | 2162 | NA | NA |
0,25 Gy | Calibrazione | HC025Gy | 1826 | NA | NA |
0,5 Gy | Calibrazione | HC05Gy | 1054 | CNL05Gy | 1532 |
0.75 Gy | Calibrazione | HC075Gy | 1233 | NA | NA |
1 Gy | Calibrazione | HC1Gy | 1566 | CNL1Gy | 841 |
2 Gy | Calibrazione | HC2Gy | 1147 | CNL2Gy | 996 |
3 Gy | Calibrazione | HC3Gy | 1212 | CNL3Gy | 1188 |
4 Gy | Calibrazione | HC4Gy | 909 | CNL4Gy | 1635 |
5 Gy | Calibrazione | HC5Gy | 1019 | NA | NA |
3.1 Gy | Prova | HCS01 | 540 | CNLS01 | 500 |
2.3 Gy | Prova | HCS08 | 637 | CNLS08 | 500 |
1.4 Gy | Prova | HCS10 | 708 | NA | NA |
1.8 Gy | Prova | HCS04 | 600 | CNLS04 | 957 |
2.8 Gy | Prova | HCS05 | 1136 | CNLS05 | 1527 |
3.4 Gy | Prova | HCS07 | 477 | CNLS07 | 735 |
Tabella 1: Fonti di dati di immagine forniti da HC e CNL per la valutazione del Software.
Nota: Modificata dalla tabella 1 in Rogan et al., 20164. Solo immagini manualmente preselezionati erano precedentemente disponibili a noi da CNL. Immagini non filtrate sono diventato disponibili e immagine conteggi vengono aggiornati di conseguenza. Inoltre, recentemente acquisito campioni HC (0.25Gy, 0.75Gy e 5Gy) sono presentati qui.
Selezione automatica dell'immagine nei campioni
Qualità dell'immagine è fondamentale per la corretta rilevazione di DC in analisi DC. Selezione delle immagini da specialisti citogenetici viene solitamente eseguita manualmente in analisi DC convenzionale. ADCI utilizza criteri quantitativi immagine per selezionare automaticamente le immagini prima di calcolo frequenza DC16. Gli utenti possono o filtro immagini basato su morfologie di cromosoma specifico e/o celle di ordinamento secondo proporzioni relative delle lunghezze degli oggetti secondo note lunghezze di citogenetica definiti gruppi di cromosomi in un cariotipo umano normale (definito il Metodo di distanza di gruppo-bin). I filtri morfologici disponibili utilizzano soglie invariante di scala per rifiutare immagini delle cellule con cromosoma incompleta di moda o con più metafasi, con i cromosomi di prometaphase, con prominente sorella dissociazione tra cromatidi fratelli, con altamente piegato e torto cromosomi, con oggetti che hanno contorni lisci caratteristici dei nuclei intatti e quelli in cui meno oggetti sono riconosciuti come cromosomi. Figura 3 (a) e (b) mostrare esempi di immagini selezionate, mentre nella figura 3 (c) e (d) sono esempi di immagini che vengono filtrati dal software. Queste immagini sono derivate dall'esempio HCS05 (descritto nella tabella 1) e vengono selezionate per il modello di selezione di immagini predefinito che truppa tutte le immagini di distanza del gruppo di bin, quindi seleziona le migliori 250 immagini. Cromosomi in Figura 3 (A), (b) sono ben separati e morfologia soddisfacente per mostre. Figura 3 (c) contiene un numero eccessivo di cluster di cromosoma sovrapposte. Figura 3 (d) spettacoli severa sorella separazione tra cromatidi fratelli. Sorella cromatidi sono completamente separati per almeno 8 dei cromosomi e le costrizioni centromeriche sono ambigue nella maggior parte dei altri cromosomi.
Figura 3: esempi di immagini di metafase nel campione HCS05 (ingrandimento: 63 X), Unselected sia selezionato dal modello 'Gruppo Bin distanza, Top 250 immagini'. Immagini selezionate (A) e (B) . (C) e (D) sono immagini che sono state eliminate dal modello. (C) è stato escluso perché conteneva troppi cromosomi sovrapposti e (D) aveva un numero eccessivo di separati sorella cromatidi. Clicca qui per visualizzare una versione più grande di questa figura.
Gli effetti dell'applicazione di questi modelli di selezione immagine è evidente esaminando il livello di probabilità di rilevazione di DC in un campione. Occorrenze di DCs in una popolazione delle cellule da un campione irradiato seguano una distribuzione di Poisson. Il test di bontà di adattamento di chi-quadrato confronta la distribuzione di frequenza osservata di DC per la misura prevista per la distribuzione di Poisson. Modelli che correttamente i dati di esempio di filtro presentano frequenze DC non significativamente differenti dalle Poisson previsto valori derivati (in genere il livello di significatività > 0,01). Figura 4 Visualizza le occorrenze DC e il corrispondente si adatta alla distribuzione di Poisson per il campione di HC4Gy di tutte le immagini vs solo immagini selezionate dal modello "distanza del gruppo di bin, top 250 immagini". Figura 4 (b) dimostra una migliore vestibilità per la distribuzione di Poisson. Il p-valore (0,36) di set filtrato di immagini supera notevolmente quello della distribuzione DC non filtrata in Figura 4 (a). Al 5% o 1% i livelli di significato, il campione non filtrato in Figura 4 (a) è meno affidabile, poiché contiene dati di DC qualità inferiori, come l'ipotesi di una distribuzione di Poisson di DCs viene rifiutata.
Figura 4: Screenshots delle frequenze proporzionate DC adatta a Poisson Dstributions di HC4Gy di esempio nel Software. (A) tutte le immagini sono inclusi, (B) solo immagini selezionate dal modello (distanza del gruppo di bin, top 250 immagini) sono inclusi. La leggenda (in alto a destra) indica le statistiche della fit per la distribuzione di Poisson (indice di dispersione, test di Mu e Lambda) e la bontà di chi-quadrato di prova adatta (valore p) Clicca qui per visualizzare una versione più grande di questa figura.
Rilevamento di cromosoma dicentrico (DC)
Rilevazione accurata di DC è il requisito fondamentale dei prerequisiti dell'ADCI. Correttamente rilevato DCs e quelli persi dal software sono rispettivamente definiti come veri positivi (TPs) e falsi negativi (FNs). Oggetti che non sono controller di dominio, ma erroneamente rilevato come controller di dominio, sono definiti come falsi positivi (FPs). FPs sono monocentric cromosomi, frammenti dei cromosomi, separata sorella cromatidi, cromosoma sovrapposte cluster e gli oggetti non-cromosomici. La figura 5 Mostra i risultati di rilevamenti di DC in due immagini di metafase. Oggetti 1 e 3 sono TPs, mentre oggetto 4 è un FP che comprende due cromosomi distinti monocentric congiunti lungo loro bracci corti. In Figura 5 (a), oggetto 2 era originariamente un FP, ma successivamente rettificato da FP filtri nel software. Oggetto 5 e 6 di oggetto in Figura 5 (b) sono esempi di probabili del FNs.
Figura 5: Screenshot per indicare la classificazione di cromosoma di metafase di potenziali DCs. (A) Un'immagine campione CNL1Gy (ingrandimento: 63 X) Risultati 1 TP, oggetto "1"; e 1 corretto FP, oggetto "2". (B) un'immagine campione CNL4Gy (ingrandimento: 63 X) Risultati 1 TP, oggetto "3"; 1 FP, oggetto "4"; e 2 potenziali FNs, oggetti "5" e "6". TPs, FPs corretto, normale monocentric e cromosomi non classificati sono descritte rispettivamente con contorni rossi, gialli, verdi e blu. Clicca qui per visualizzare una versione più grande di questa figura.
Stima della dose dei campioni di prova
Il risultato finale dell'analisi ADCI sono la dose stima dei campioni dedotti dalle curve di calibrazione. Le valutazioni di dose apportate dal software per i campioni di prova nella tabella 1 sono indicate nelle tabelle 2 e 3. Per confronto, la dose di radiazione fisica emessa e sono indicate le dosi manuale ha ottenute dagli esperti di HC per i campioni HCS01, HCS08 e HCS10. Allo stesso modo, la fisica e manuale segnato le dosi da esperti CNL sono mostrati per CNLS04, CNLS05 e CNLS07.
Figura 6 illustra le curve di calibrazione con le stime delle dosi di radiazione per campioni di laboratorio di Health Canada biodosimetry HCS01, HCS08, HCS10, HCS04, HCS05 e HCS07. La curva di taratura viene generata utilizzando campioni HC0Gy, HC1Gy, HC2Gy, HC3Gy e HC4Gy. Il modello di selezione immagine contenente 3 filtri basati su Z-punteggio + "gruppo bin distanza, immagini top 250" viene applicato a tutti i campioni. Le stime delle dosi insieme ad analisi statistiche associate sono riportati nella tabella 2.
Figura 6: Screenshot di stima di Dose di campioni di prova HC. Quadrati neri rappresentano campioni di calibrazione. Immagini in esempi di prova e taratura vengono selezionati dal modello (3 filtri FP + distanza del gruppo di bin, top 250 immagini). Linee spesse tratteggiate rappresentano il mapping di DCs/metafase tramite la curva di taratura per dose stimata. Sottili linee punteggiate indicano superiore e abbassare i limiti di confidenza del 95% di DCs/metafase. I codici dei campioni di prova colore: rosso brillante, HC S01 (dose fisica: 3.1Gy, dose HC dedotto: 3.4Gy, ADCI: 3Gy); verde scuro, HC S04 (dose fisica: 1.8Gy, ADCI: 1.85Gy); blu brillante, HC S05 (dose fisica: 2.8Gy, ADCI: 2.95Gy); blu scuro, HC S07 (dose fisica: 3.4Gy, ADCI: 2.35Gy); rosso scuro, HC S08 (dose fisica: 2.3Gy, dose HC dedotto: 2.5Gy, ADCI: 2Gy); verde brillante, HC S10 (dose fisica: 1.4Gy, dose HC dedotto: 1.4Gy, ADCI: 0.95Gy). Clicca qui per visualizzare una versione più grande di questa figura.
Campioni | Dose fisica | HC dedotto Dose | ADCI stimato Dose | Dose stimata LCL | Dose stimata UCL | P-valore * |
HCS01 | 3.1 | 3.4 | 3 | 2.3 | 3.8 | 0,117 |
HCS08 | 2.3 | 2.5 | 2 | 1.4 | 2.7 | 0,815 |
HCS10 | 1.4 | 1.4 | 0.95 | 0,5 | 1,55 | 0,211 |
HCS04 | 1.8 | NA | 1.85 | 1.25 | 2,55 | 0,0293 |
HCS05 | 2.8 | NA | 2.95 | 2.25 | 3.75 | 0.00354 |
HCS07 | 3.4 | NA | 2.35 | 1.7 | 3.1 | 0,0002 |
Tabella 2: Dose stima risultati di campioni di prova HC.
Nota: Modificata dalla tabella 3 in Rogan et al., 20164. Le stime delle dosi ADCI precedentemente segnalati sono stati basati su immagini non filtrate e adattamento della curva è stata effettuata usando il metodo dei minimi quadrati. Qui, la curva di taratura era adatta utilizzando il metodo di massima verosimiglianza e un modello di selezione immagine contenente 3 filtri FP + "gruppo bin distanza, immagini top 250" è stato applicato prima stima di dose. Dose stimata UCLe LCL si riferiscono all'alta stima della dose e abbassare i limiti di confidenza del 95% sulla base della natura di Poisson di rendimento DC. * Il bontà quadrato chi di adattamento per distribuzione di Poisson teorica; NA: Risultati della dose manualmente dedotto non erano forniti.
Stime di dose di radiazione per campioni da Canadian Nuclear CNLS04 di laboratori, CNLS05, CNLS07, CNLS01 e CNLS08 sono illustrati nella Figura 7. La curva di taratura viene generata utilizzando campioni CNL0Gy, CNL0.5Gy, CNL1Gy, CNL2Gy, CNL3Gy e CNL4Gy. Abbiamo applicato un modello di selezione immagine composto da 6 FP filtri per tutti i campioni. I risultati con analisi statistiche sono mostrati nella tabella 3.
Figura 7: Screenshot di stima di Dose di campioni CNL. Quadrati neri rappresentano campioni di calibrazione. Immagini in esempi di prova e taratura vengono selezionati con 6 filtri FP. Linee spesse tratteggiate rappresentano il mapping di DCs/metafase tramite la curva di taratura per dose stimata. Sottili linee punteggiate indicano superiore e abbassare i limiti di confidenza del 95% di DCs/metafase. I codici dei campioni di prova colore: rosso brillante, CNL S04 (dose fisica: 1.8Gy, dose CNL dedotto: 1.7Gy, ADCI: 1.95Gy); rosso scuro, S05 CNL (dose fisica: 2.8Gy, dose CNL dedotto: 2.7Gy, ADCI: 3.05Gy); verde brillante, CNL S07 (dose fisica: 3.4Gy, dose CNL dedotto: 3.1Gy, ADCI: 3.4Gy); verde scuro, CNL S01 (dose fisica: 3.1Gy, ADCI: 3.75Gy); blu, CNL S08 (dose fisica: 2.3Gy, ADCI: 2.8Gy). Clicca qui per visualizzare una versione più grande di questa figura.
Campioni | Dose fisica | CNL dedotto Dose | ADCI stimato Dose | Dose stimata LCL | Dose stimata UCL | P-valore * |
CNLS04 | 1.8 | 1.7 | 1.95 | 1.25 | 2.45 | 0,0545 |
CNLS05 | 2.8 | 2.7 | 3,05 | 2,75 | 3,35 | 0.325 |
CNLS07 | 3.4 | 3.1 | 3.4 | 3 | 3.75 | 0.473 |
CNLS01 | 3.1 | NA | 3.75 | 3,35 | > 4 | 7.63E-11 |
CNLS08 | 2.3 | NA | 2.8 | 2.25 | 3.3 | 0,777 |
Tabella 3: Risultati di stima di dose di CNL testare campioni.
Nota: Modificata dalla tabella 3, Rogan et al., 20164. Le stime delle dosi ADCI precedentemente segnalati sono stati basati su non filtrato (HC) o manualmente selezionato (CNL) immagini e adattamento della curva è stata effettuata usando il metodo dei minimi quadrati. Qui, la curva di taratura era adatta utilizzando il metodo di massima verosimiglianza e un modello di selezione immagine contenente 3 filtri FP + "gruppo bin distanza, immagini top 250" è stato applicato prima stima di dose. Dose stimata UCL e LCL, rispettivamente, si riferiscono per dose stimata superiore e abbassare i limiti di confidenza del 95% sulla base della natura di Poisson di rendimento DC.
* Il bontà quadrato chi di adattamento per distribuzione di Poisson teorica; NA: Risultati della dose manualmente dedotto non erano disponibili.
Stima della dose di radiazione all'interno della gamma lineare della curva di calibrazione (< 1 Gy) possono essere effettuati con il software, tuttavia un valore di 1.0 è consigliato di Sigma è quello di ridurre ulteriormente la frequenza di segnalarli DCs (Figura 8).
Figura 8: Screenshots delle due curve di taratura derivato da HC campioni di taratura presso Sigma diversi valori. (A) campioni di calibrazione HC: 0Gy, 2Gy, 3Gy, 4Gy e 5Gy presso Sigma = 1.5. Campioni di calibrazione di HC (B) : 0Gy, 0.25Gy, 0.5Gy, 0.75Gy, 1Gy, 2Gy, 3Gy, 4Gy, e 5Gy utilizzando SVM Sigma = 1.0. Clicca qui per visualizzare una versione più grande di questa figura.
Queste analisi indicano che ci sono piccole, ma accettabili differenze tra fisica e biologicamente inferita dose interpretate dagli esperti e dal software. La differenza tra manuale o software stima dalla dose fisica si riferisce a come "errore". L'errore nelle dosi dedotte di campioni manualmente segnati da HC e CNL è ≤ 0.3 Gy. Trattamento automatizzato dal software è meno preciso di quelli esperti, ma generalmente entro i limiti di triage AIEA-specificato di ± 0,5 Gy3. Per la maggior parte dei campioni nelle tabelle 2 e 3, il software prodotto risultati corretti all'interno di questa soglia. Tuttavia, HCS07 e CNLS01 presentano una scarsa bontà di adattamento per la distribuzione di Poisson, suggerendo che c'erano problemi potenziali in immagine e qualità DC in questi campioni che non sono stati risolti dall'applicazione dell'immagine e modelli di selezione FP. La soglia di rilevanza del valore p sembra essere eccessivamente severe in caso di HCS05, dove il software determinato con precisione la dose corretta.
Funzionalità e limitazioni del software
Il protocollo descritto in questo documento introduce la procedura graduale tipica usata in ADCI per importare ed elaborare immagini di metafase citogenetica, creare radiazioni curve di calibrazione e stimare la dose biologica nelle persone o nei campioni esposti a sconosciuto livelli di radiazione. Tuttavia, non è necessario eseguire queste istruzioni in sequenza. Ad esempio, molti campioni di prova della dose sconosciuta possono essere elaborati e analizzati utilizzando la stessa curva di taratura pre-calcolate. Inoltre, una volta completata l'elaborazione, la selezione delle immagini e DC modelli di filtraggio può essere iterata da parte dell'utente. Applicazione di un modello di selezione immagine appropriata dipende dalle caratteristiche e qualità dei dati dell'immagine di metafase, che a sua volta si basa sia sul protocollo di laboratorio utilizzato per preparare le cellule e i criteri di rigore utilizzati per selezionare le celle con automatizzato sistemi di cattura di metafase. Morfologie di cromosoma saranno diverso tra biodosimetry e laboratori di citogenetiche, e così, i modelli di selezione immagine dovrebbero essere valutati dall'utente per determinare se i modelli di selezione immagine predefiniti forniti con il software sarà sufficienti per produrre le stime delle dosi precise, o se personalizzato modelli definiti dall'utente soglie devono essere creati. Basandoci sulla nostra esperienza, l'efficacia di modelli di selezione immagine risente l'origine e la qualità delle immagini delle cellule. Gli utenti possono progettare i propri criteri di selezione immagine utilizzando diverse combinazioni di filtri per eliminare il falso positivo DCs e modelli di selezione immagine e i valori di soglia corrispondenti per selezionare le immagini desiderate. C'è flessibilità in ingresso di curve di calibrazione e stima della dose, come coefficienti della curva lineare-quadratica e frequenze DC possono essere modificati o inseriti manualmente.
Anche se il software è completamente automatizzato, immagini possono essere recensione e selezionati manualmente. Questa funzionalità è disponibile per includere o rimuovere immagini singolarmente elaborate attraverso la funzione di Visualizzatore di microscopio nella GUI principale. Tuttavia, grazie all'automazione, il software è molto più efficiente rispetto con il manuale notante di immagini di metafase e contando DCs. Un campione costituito da 1000 immagini possa essere elaborato in 20 (tiff) per 40 min (jpg) su una workstation di prestazioni multi-core. Questo software sarà particolarmente utile in situazioni critiche o laboriosi, come eventi in cui individui multipli sono stati esposti o sono stati ritenuti sospetto di essere stati esposti a radiazioni, o dove il tempo-sensibile diagnosi e trattamento le decisioni sono fondamentali.
Rilevamento preciso e accurato elevato throughput di DCs, come pure la stima di dose sono necessari per la valutazione automatica di radiazione. Altre alternative disponibili per il software non soddisfare entrambi questi requisiti. Un'analisi citogenetica utente-assistita, basata sulle immagini (DCScore, Metasystems17) sistema richiede la verifica manuale del candidato DCs, a causa di un errore elevato tasso attribuibile alla non corretta si sovrappone tra cromosomi e il sistema non determina dose di radiazioni. DCScore non sarebbero altrettanto efficaci come ADCI in un evento di radiazioni che coinvolgono un gran numero di soggetti potenzialmente esposti. Sistemi di microscopia di grande apertura in grado di raccogliere immagini di più cellule di metafase18, tuttavia, non analizzarli. "CABAS"19 e "Stima di Dose"20 software può generare calibrazione curve e stima della dose, ma non segnano DCs. Altri saggi di biodosimetry che non sono basati su analisi DC includono H2AX fluorescenza, sonde di ibridazione in situ di fluorescenza con DNA mirati a specifici cromosomi, espressione genica, micronucleo e l'urina e biomarcatori respiratorie. Questi metodi sono meno specifici e sensibili per le radiazioni ionizzanti, possono essere più costosi, in alcuni casi, sono più che richiede tempo e generalmente non sono stati standardizzati attraverso molteplici laboratori di riferimento. La maggior parte di queste tecniche non rilevano radiazioni stabile risposte, quindi non possono essere utilizzati per la valutazione a lungo termine (> post-esposizione di 7 giorni) della dose di radiazione. Al contrario, questo può valutare gli individui fino a 90 giorni di post-esposizione e in grado di elaborare dati da qualsiasi microscopio del laboratorio di citogenetica sistema di imaging. Tuttavia, se un campione è disegnato > 4 settimane post-esposizione, sensibilità è ridotta a causa del decadimento di aberrazioni dicentrico1,2,3 e il software non risolve attualmente frequenze DC per ritardi nel prelievo individui esposti.
Questo software ha alcune limitazioni. Modelli di selezione immagine esistente selezionare immagini di metafase accettabile per la maggior parte, ma in alcuni casi, non riescono a eliminare immagini insoddisfacente, che possono ridurre la precisione del rilevamento di DC. È ancora una questione aperta come progettare un modello di selezione immagine soddisfacente che elimina tutte le cellule in metafase inadatto. Il software fornisce una stima accurata per campioni esposti a dosi di radiazioni elevate (≥ 2 Gy). Nonostante i notevoli progressi nel ridurre il numero di falsi positivi DCs16, questi oggetti non sono stati eliminati. Abbassare le cellule di metafase di qualità a bassa dose di radiazioni (soprattutto < 1 Gy) sono più inclini a rilevamento del falso positivo DC. Di conseguenza, campioni di basso dosaggio non sono stati inclusi quando si genera la curva di taratura utilizzata per la stima della dose dei campioni di prova HC. Tuttavia, se una curva contenente campioni di basso dosaggio è desiderata, un valore inferiore di SVM Sigma riduce la conta falsi positivi nei campioni di basso dosaggio ma può provocare calo dei rendimenti DC in campioni di dose elevata. Figura 8 confronta la curva HC utilizzata per la stima della dose (Sigma = 1.5) con una curva di calibrazione misura con campioni di dose aggiuntiva di basso valore inferiore del sigma SVM (1.0). Nei campioni con un numero insufficiente di cellule in metafase e/o immagini di metafase di scarsa qualità, non è possibile stimare le esposizioni biologiche a basso dosaggio, potenzialmente con conseguente deviazioni da fisici dose superiore a 0,5 Gy.
Il software non può valutare con precisione tipi di radiazioni se le loro curve dose-risposta migliore adattano un modello lineare o quasi lineare. Finora, è stato testato solo con campioni esposti a X e raggi Gamma. Se un'altra fonte di radiazione è esaminata, gli utenti devono verificare i campioni di test e calibrazione sono esposti allo stesso tipo di radiazioni. Il software utilizza la massima verosimiglianza o minimi quadrati raccordo per costruire una curva dose-risposta utilizzando un modello lineare-quadratica. Attualmente non esiste alcuna opzione per imporre una rigorosa curva lineare, appropriato per le esposizioni di particelle ad alta energia, tuttavia tale funzionalità sarà disponibile in futuro.
Sviluppo futuro
Nostri sforzi sono concentrati sul miglioramento di modelli di selezione immagine e misura della dose precisa, in particolare dei campioni esposti a dosi di radiazione bassa. Versioni successive del software fornirà errore standard misure le stime delle dosi e intervalli di confidenza sulle curve di calibrazione. Inoltre, una versione high-performance computing del software per il supercomputer Blue Gene (BG/Q, IBM) è in fase di sviluppo per la tempestiva valutazione di individui esposti in un evento di massa-incidente radiazione. Alcuni componenti del software sono già stati testati e distribuiti su questa piattaformaLass = "xrif" > 11.
PKR e JHMK stato cofondatore di CytoGnomix, che è commercializzando ADCI e detiene brevetti correlati. YL e BCS sono dipendenti di CytoGnomix. ADCI è protetto da copyright, e il metodo di localizzazione del centromero in ADCI è brevettato (Pat. N. 8.605.981; Pat tedesca. N. 112011103687).
Siamo grati al Dr. Ruth Wilkins, radiobiologia e protezione Division presso Health Canada e Farrah Flegal, Canadian Nuclear Laboratories e il loro personale di laboratorio per l'accesso ai dati immagine metafase dai loro laboratori di citogenetica biodosimetry. Questa carta è stata sostenuta da un contratto da costruire in Canada Innovation Program a CytoGnomix (numero di serie. EN579-172270/001/SC). La versione iniziale di ADCI e sviluppo di algoritmi sono stati sostenuti dal fondo di innovazione occidentale; Scienze naturali e ingegneria Research Council of Canada (NSERC Discovery Grant 371758-2009); U.S. Public Health Service (DART-DOSE CMCR, 5U01AI091173-0); la Fondazione Canada per l'innovazione; Canada ricerca sedie e CytoGnomix Inc.
Name | Company | Catalog Number | Comments |
Automated Dicentric Chromosome Identifier and Dose Estimator (ADCI) | CytoGnomix | NA | ADCI software is released in a binary installation package file for Microsoft Windows 7, 8, 8.1 and 10; 235 Mb of disk storage are required for a typical installation. The software has been tested with Intel or AMD x86-64 processors; at least 1 Gb RAM is recommended. Analyses have been benchmarked on a computer configured with an Intel I7 processor and 16 Gb RAM. Operation of ADCI requires an active license and a USB-based hardware dongle, which must remain plugged in while the software is executing. The dongle encodes the software expiry date. Each time the software is started, this date is read. The software will allow access to the program if the current date and time precedes the expiration time-date stamp. Extending an expired software license can be accomplished by obtaining a new dongle or by renewing the license with an updated key at startup. |
Digital images of metaphase cell nuclei | Examples: Metasystems, Leica Microsystems | M-Search (Metasystems), Cytovision (Leica) software | High resolution TIFF format; typically >250 digital images generated with a microscope imaging capture system (minimum 63X magnification objective, 10X magnification ocular). |
MSI Leopard Pro (recommended, optional) | Micro-Star International | MSI GP62 6QF 480CA Leopard Pro | Multi-core performance workstation. |
Richiedi autorizzazione per utilizzare il testo o le figure di questo articolo JoVE
Richiedi AutorizzazioneThis article has been published
Video Coming Soon