Testowanie oprogramowania przeprowadzono z wykorzystaniem danych obrazowych chromosomów w metafazie uzyskanych z HC i CNL. Próbki krwi napromieniowano urządzeniem XRAD-320 (promienie X 250 kV, 12,5 mA, filtracja 2 mm Al, szybkość dawki: 0,92 lub 1,7 Gy/min) skalibrowanym za pomocą komory jonizacyjnej w HC i przetworzonym w obu laboratoriach. Próbki limfocytów krwi obwodowej hodowano, utrwalano i barwiono w każdej placówce zgodnie z uznanymi protokołami3,15. Obrazy metafaz z preparatów barwionych metodą Giemsy były niezależnie rejestrowane przez każde laboratorium przy użyciu zautomatyzowanego systemu mikroskopowego. Eksperci w każdym laboratorium ręcznie oceniali DC w kilku tych próbkach, konstruowali własne krzywe kalibracyjne i szacowali dawki próbek testowych o nieznanym stopniu ekspozycji. Szczegółowy opis tych zbiorów danych przedstawiono w Tabeli 1.
| Dawka fizyczna | Cel | Przygotowanie HC | Przygotowanie CNL |
| | Nazwa powiązana | liczba obrazów | Nazwa odniesienia | liczba obrazów |
| 0 Gy | Kalibracja | HC0Gy | 731 | CNL0Gy | 798 |
| 0,1 Gy | Kalibracja | HC01Gy | 2162 | ND | N/D |
| 0,25 Gy | Kalibracja | HC025Gy | 1826 | Brak danych | ND |
| 0,5 Gy | Kalibracja | HC05Gy | 1054 | CNL05Gy | 1532 |
| 0,75 Gy | Kalibracja | HC075Gy | 1233 | Brak danych | ND |
| 1 Gy | Kalibracja | HC1Gy | 1566 | CNL1Gy | 841 |
| 2 Gy | Kalibracja | HC2Gy | 1147 | CNL2Gy | 996 |
| 3 Gy | Kalibracja | HC3Gy | 1212 | CNL3Gy | 1188 |
| 4 Gy | Kalibracja | HC4Gy | 909 | CNL4Gy | 1635 |
| 5 Gy | Kalibracja | HC5Gy | 1019 | ND | Brak danych |
| 3,1 Gy | Test | HCS01 | 540 | CNLS01 | 500 |
| 2,3 Gy | Test | HCS08 | 637 | CNLS08 | 500 |
| 1,4 Gy | Test | HCS10 | 708 | Brak danych | NA |
| 1,8 Gy | Test | HCS04 | 600 | CNLS04 | 957 |
| 2,8 Gy | Test | HCS05 | 1136 | CNLS05 | 1527 |
| 3,4 Gy | Test | HCS07 | 477 | CNLS07 | 735 |
Tabela 1: Źródła danych obrazowych dostarczone przez HC i CNL do ewaluacji oprogramowania.
Przypis: Zmodyfikowano na podstawie Tabela 1 w Rogan i wsp., 20164Wcześniej z CNL dysponowaliśmy jedynie obrazami wybranymi ręcznie. Obecnie dostępne są obrazy niefiltrowane, a liczba obrazów została odpowiednio zaktualizowana. Dodatkowo zaprezentowano tutaj nowo pozyskane próbki HC (0,25 Gy, 0,75 Gy i 5 Gy).
Automatyczny wybór obrazów w próbkach
Jakość obrazu jest kluczowa dla prawidłowej detekcji DC w analizie DC. W konwencjonalnej analizie DC wybór obrazów przez specjalistów cytogenetycznych zazwyczaj odbywa się ręcznie. ADCI wykorzystuje ilościowe kryteria obrazu do automatycznego wyboru obrazów przed obliczeniem częstotliwości DC16. Użytkownicy mogą odfiltrować obrazy na podstawie specyficznej morfologii chromosomów i/lub sortować komórki zgodnie z względnymi proporcjami długości obiektów według znanych długości zdefiniowanych cytogenetycznie grup chromosomów w normalnym ludzkim kariotypie (metoda określana jako group-bin distance). Dostępne filtry morfologiczne wykorzystują progi niezależne od skali, aby odrzucić obrazy komórek z niepełnymi zestawami chromosomów, z wielokrotnymi metafazami, z chromosomami w prometafazie, z wyraźną dysocjacją chromatyd siostrzanych, z silnie zgiętymi i skręconymi chromosomami, z obiektami o gładkich konturach charakterystycznych dla nienaruszonych jąder oraz te, w których rozpoznano zbyt mało obiektów jako chromosomy. Rysunek 3 (a) i (b) przedstawia przykłady wybranych obrazów, natomiast Rysunek 3 (c) i (d) to przykłady obrazów odfiltrowanych przez oprogramowanie. Obrazy te pochodzą z próbki HCS05 (opisanej w Tabeli 1) i zostały wybrane przez predefiniowany model wyboru obrazów, który szereguje wszystkie obrazy według group bin distance, a następnie wybiera 250 najlepszych obrazów. Chromosomy na Rysunku 3 (a), (b) są dobrze odseparowane i wykazują zadowalającą morfologię. Rysunek 3 (c) zawiera nadmierną liczbę nakładających się skupisk chromosomów. Rysunek 3 (d) pokazuje silną separację chromatyd siostrzanych. Chromatydy siostrzane są całkowicie rozdzielone w przypadku co najmniej 8 chromosomów, a przewężenia centromeryczne są niejednoznaczne w większości pozostałych chromosomów.

Rycina 3: Przykłady obrazów metafazy w próbce HCS05 (powiększenie: 63X), zarówno niewybrane, jak i wybrane przez model „Group Bin Distance, Top 250 Images”. (A) i (B) wybrane obrazy. (C) i (D) są obrazami, które zostały odrzucone przez model. Obraz (C) wykluczono ze względu na zbyt dużą liczbę nakładających się chromosomów, a obraz (D) ze względu na nadmierną liczbę rozdzielonych chromatyd siostrzanych. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Efekty zastosowania tych modeli selekcji obrazów są widoczne po zbadaniu poziomu ufności detekcji DC w próbce. Występowanie DC w populacji komórek z napromieniowanej próbki odpowiada rozkładowi Poissona. Test chi-kwadrat dobroci dopasowania porównuje obserwowany rozkład częstości DC z oczekiwanym dopasowaniem do rozkładu Poissona. Modele, które prawidłowo filtrują dane z próbki, wykazują częstości DC, które nie różnią się znacząco od wartości wynikających z rozkładu Poissona (zazwyczaj poziom istotności >0,01). Rysunek 4 przedstawia występowanie DC oraz odpowiadające im dopasowania do rozkładów Poissona dla próbki HC4Gy obejmującej wszystkie obrazy w porównaniu z obrazami wybranymi przez model „group bin distance, top 250 images”. Rysunek 4 (b) pokazuje lepsze dopasowanie do rozkładu Poissona. Wartość p (0,36) dla przefiltrowanego zestawu obrazów znacząco przekracza wartość dla niefiltrowanego rozkładu DC na Rysunku 4 (a). Przy poziomie istotności 5% lub 1% niefiltrowana próbka na Rysunku 4 (a) jest mniej wiarygodna, ponieważ zawiera dane DC niższej jakości, a hipoteza zerowa o rozkładzie Poissona dla DC zostaje odrzucona.

Rycina 4: Zrzuty ekranu z oprogramowania przedstawiające dopasowanie proporcjonalnych częstotliwości DC do rozkładów Poissona dla próbki HC4Gy. (A) Wszystkie obrazy zostały dołączone, (B) Uwzględniono wyłącznie obrazy wybrane przez model (odległość grup binarnych, 250 najlepszych obrazów). Legenda (w prawym górnym rogu) przedstawia statystyki dopasowania do rozkładu Poissona (indeks dyspersji, test Mu oraz Lambda) oraz test chi-kwadrat dla dobroci dopasowania (wartość p). Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej figury.
Wykrywanie chromosomów dicentrycznych (DC)
Dokładna detekcja centromer dycentrycznych (DC) jest krytycznym warunkiem wstępnym ADCI. Prawidłowo wykryte DC oraz te, które zostały pominięte przez oprogramowanie, definiuje się odpowiednio jako prawdziwie dodatnie (TP) i fałszywie ujemne (FN). Obiekty, które nie są DC, ale zostały błędnie wykryte jako DC, określa się mianem fałszywie dodatnich (FP). Do FP należą chromosomy monocentryczne, fragmenty chromosomów, rozdzielone chromatydy siostrzane, nakładające się skupiska chromosomów oraz obiekty niechromosomowe. Rysunek 5 przedstawia wyniki detekcji DC na dwóch obrazach metafaz. Obiekty 1 i 3 są TP, natomiast obiekt 4 to FP składający się z dwóch odrębnych chromosomów monocentrycznych połączonych w obrębie ich krótkich ramion. Na Rysunku 5 (a) obiekt 2 był pierwotnie FP, ale został następnie skorygowany przez filtry FP w oprogramowaniu. Obiekt 5 i obiekt 6 na Rysunku 5 (b) to prawdopodobnie przykłady FN.

Rysunek 5: Zrzuty ekranu wskazują klasyfikację chromosomów w metafazie potencjalnych komórek DC. (A) Obraz w próbce CNL1Gy (powiększenie: 63X) przedstawiający 1 TP, obiekt „1”, oraz 1 skorygowany FP, obiekt „2”. (B) Obraz w próbce CNL4Gy (powiększenie: 63X) przedstawiający 1 TP, obiekt „3”; 1 FP, obiekt „4” oraz 2 potencjalne FN, obiekty „5” i „6”. TP, skorygowane FP, normalne chromosomy monocentryczne i chromosomy niezaklasyfikowane są odpowiednio obrysowane konturami w kolorze czerwonym, żółtym, zielonym i niebieskim. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej figury.
Szacowanie dawki próbek testowych
Końcowym wynikiem analiz ADCI są szacunki dawek dla próbek wyznaczone na podstawie krzywych kalibracyjnych. Szacunki dawek wykonane przez oprogramowanie dla próbek testowych wymienionych w Tabeli 1 przedstawiono w Tabelach 2 i 3. Dla porównania podano fizyczną dawkę promieniowania emitowaną oraz dawki określone manualnie przez ekspertów z HC dla próbek HCS01, HCS08 i HCS10. Analogicznie, dla próbek CNLS04, CNLS05 i CNLS07 przedstawiono dawki fizyczne i określone manualnie przez ekspertów z CNL.
Rysunek 6 przedstawia krzywe kalibracyjne z szacunkami dawek promieniowania dla próbek z laboratorium biodozymetrii Health Canada: HCS01, HCS08, HCS10, HCS04, HCS05 i HCS07. Krzywa kalibracyjna została wygenerowana przy użyciu próbek HC0Gy, HC1Gy, HC2Gy, HC3Gy i HC4Gy. Do wszystkich próbek zastosowano model wyboru obrazów zawierający 3 filtry oparte na wartości Z-score + „odległość kosza grupy, 250 najlepszych obrazów”. Szacunki dawek wraz z powiązanymi analizami statystycznymi przedstawiono w Tabeli 2.

Rysunek 6: Zrzut ekranu z szacowania dawki próbek testowych HC. Czarne kwadraty reprezentują próbki kalibracyjne. Obrazy w próbkach testowych i kalibracyjnych są wybierane przez model (3 filtry FP + odległość koszy grup, 250 najlepszych obrazów). Grube linie przerywane reprezentują mapowanie DC/metafaz poprzez krzywą kalibracyjną na szacowaną dawkę. Cienkie linie przerywane oznaczają górną i dolną 95% granicę ufności DC/metafaz. Kody kolorystyczne próbek testowych: jasnoczerwony, HC S01 (dawka fizyczna: 3,1 Gy, dawka wywnioskowana HC: 3,4 Gy, ADCI: 3 Gy); ciemnozielony, HC S04 (dawka fizyczna: 1,8 Gy, ADCI: 1,85 Gy); jasnoniebieski, HC S05 (dawka fizyczna: 2,8 Gy, ADCI: 2,95 Gy); ciemnoniebieski, HC S07 (dawka fizyczna: 3,4 Gy, ADCI: 2,35 Gy); ciemnoczerwony, HC S08 (dawka fizyczna: 2,3 Gy, dawka wywnioskowana HC: 2,5 Gy, ADCI: 2 Gy); jasnozielony, HC S10 (dawka fizyczna: 1,4 Gy, dawka wywnioskowana HC: 1,4 Gy, ADCI: 0,95 Gy). Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.
| Próbki | Dawka fizyczna | Dawka wnioskowana HC | Dawka szacowana ADCI | Dolna granica ufności dawki szacowanej (LCL) | Górna granica ufności dawki szacowanej (UCL) | Wartość p* |
| HCS01 | 3.1 | 3.4 | 3 | 2.3 | 3.8 | 0.17 |
| HCS08 | 2.3 | 2.5 | 2 | 1.4 | 2.7 | 0.815 |
| HCS10 | 1.4 | 1.4 | 0.95 | 0.5 | 1.5 | 0.21 |
| 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.002 |
Tabela 2: Wyniki szacowania dawki dla próbek testowych HC.
Przypis: Zmodyfikowano na podstawie Tabeli 3 w pracy Rogan et al., 20164. Wcześniej raportowane szacunki dawki ADCI opierały się na obrazach niefiltrowanych, a dopasowanie krzywej wykonano metodą najmniejszych kwadratów. Tutaj krzywą kalibracyjna dopasowano metodą największej wiarygodności, a przed szacowaniem dawki zastosowano model wyboru obrazów zawierający 3 filtry FP + „group bin distance, top 250 images”. Szacowane UCL i LCL dawki odnoszą się do górnej i dolnej 95% granicy ufności szacowania dawki na podstawie rozkładu Poissona dla wydajności DC. * Chi kwadrat dopasowania do teoretycznego rozkładu Poissona; NA: Nie podano wyników dawki wyznaczonej ręcznie.
Szacunkowe dawki promieniowania dla próbek z Canadian Nuclear Laboratories CNLS04, CNLS05, CNLS07, CNLS01 i CNLS08 przedstawiono w Rycina 7Krzywa kalibracyjna została wygenerowana przy użyciu próbek CNL0Gy, CNL0.5Gy, CNL1Gy, CNL2Gy, CNL3Gy oraz CNL4Gy. Do wszystkich próbek zastosowano model selekcji obrazów składający się z 6 filtrów FP. Wyniki wraz z analizami statystycznymi przedstawiono w Tabela 3.

Rycina 7: Zrzut ekranu z szacowania dawki dla próbek testowych CNL. Czarne kwadraty reprezentują próbki kalibracyjne. Obrazy w próbkach testowych i kalibracyjnych są wybierane przy użyciu 6 filtrów FP. Pogrubione linie przerywane reprezentują mapowanie liczby DC/metafaz poprzez krzywą kalibracyjną na szacowaną dawkę. Cienkie linie przerywane oznaczają górne i dolne 95% granice ufności dla liczby DC/metafaz. Kody kolorystyczne próbek testowych: jasnoczerwony, CNL S04 (dawka fizyczna: 1,8 Gy, dawka wnioskowana CNL: 1,7 Gy, ADCI: 1,95 Gy); ciemnoczerwony, CNL S05 (dawka fizyczna: 2,8 Gy, dawka wnioskowana CNL: 2,7 Gy, ADCI: 3,05 Gy); jasnozielony, CNL S07 (dawka fizyczna: 3,4 Gy, dawka wnioskowana CNL: 3,1 Gy, ADCI: 3,4 Gy); ciemnozielony, CNL S01 (dawka fizyczna: 3,1 Gy, ADCI: 3,75 Gy); niebieski, CNL S08 (dawka fizyczna: 2,3 Gy, ADCI: 2,8 Gy). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
| Próbki | Dawka fizyczna | Wywnioskowana dawka CNL | Szacowana dawka ADCI | Szacowana dawka LCL | Szacowana dawka UCL | wartość p* |
| 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 | ND | 3.75 | 3.35 | >4 | 7,63E-1 |
| CNLS08 | 2.3 | ND | 2.8 | 2.25 | 3.3 | 0.777 |
Tabela 3: Wyniki szacowania dawki dla próbek testowych CNL.
Przypisy: Zmodyfikowano na podstawie Tabela 3, Rogan i in., 20164Wcześniej raportowane szacunki dawki ADCI opierały się na obrazach niefiltrowanych (HC) lub selekcjonowanych ręcznie (CNL), a dopasowanie krzywej przeprowadzono metodą najmniejszych kwadratów. W niniejszej pracy krzywą kalibracyjną dopasowano metodą największej wiarygodności, a przed oszacowaniem dawki zastosowano model selekcji obrazów zawierający 3 filtry FP oraz parametr „group bin distance, top 250 images”. Szacowane UCL i LCL dawki odnoszą się odpowiednio do górnej i dolnej 95% granicy ufności dawki szacowanej, opartej na rozkładzie Poissona dla uzysku DC.
* Test chi-kwadrat zgodności z teoretycznym rozkładem Poissona; NA: Wyniki ręcznie wyznaczonej dawki były niedostępne.
Szacowanie dawki promieniowania w zakresie liniowym krzywej kalibracyjnej (<1 Gy) można przeprowadzić za pomocą oprogramowania, jednak zaleca się zastosowanie wartości Sigma wynoszącej 1,0, aby dodatkowo zmniejszyć częstotliwość błędnej klasyfikacji DC (Rysunek 8).

Rysunek 8: Zrzuty ekranu dwóch krzywych kalibracyjnych uzyskanych z próbek kalibracyjnych HC przy różnych wartościach sigma. (A) Próbki kalibracyjne HC: 0 Gy, 2 Gy, 3 Gy, 4 Gy oraz 5 Gy przy Sigma = 1,5. (B) Próbki kalibracyjne HC: 0 Gy, 0,25 Gy, 0,5 Gy, 0,75 Gy, 1 Gy, 2 Gy, 3 Gy, 4 Gy i 5 Gy przy użyciu SVM Sigma = 1,0. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Analizy te wskazują, że istnieją niewielkie, lecz dopuszczalne różnice między dawką fizyczną a dawką wnioskowaną biologicznie, interpretowaną przez ekspertów oraz przez oprogramowanie. Różnica między szacowaniem manualnym a programowym względem dawki fizycznej jest określana jako „błąd”. Błąd w dawkach wnioskowanych dla próbek ocenianych manualnie przez HC i CNL wynosi ≤0.3 Gy. Przetwarzanie zautomatyzowane przez oprogramowanie jest mniej dokładne niż ocena ekspertów, lecz zazwyczaj mieści się w granicach triażu określonych przez IAEA, wynoszących ± 0.5 Gy3. Dla większości próbek testowych w Tabelach 2 i 3 oprogramowanie uzyskało poprawne wyniki w obrębie tego progu. Jednakże HCS07 i CNLS01 wykazują słabe dopasowanie do rozkładu Poissona, co sugeruje potencjalne problemy z jakością obrazu i DC w tych próbkach, które nie zostały rozwiązane poprzez zastosowanie modeli selekcji obrazów i FP. Próg istotności wartości p wydaje się być zbyt rygorystyczny w przypadku HCS05, gdzie oprogramowanie dokładnie określiło poprawną dawkę.