$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Ilościowa ilościowa objętość jądrowa w sferoidach nowotworowych w 3D
Zautomatyzowana mikroskopia 3D (z-stacks) umożliwia naukowcom wizualizację złożonych systemów wielokomórkowych, takich jak organoidy. Aby ilościowo określić właściwości sygnału morfologicznego i fluorescencyjnego na poziomie pojedynczych komórek, często wymagana jest segmentacja komórek w 3D. Poniższy przykład pokazuje, jak Cell-ACDC może segmentować jądra w organoidy zawierające tysiące komórek z kanału barwienia jądrowego i ilościowo określać ich objętość (dane z Ref.9). Segmentacja była przeprowadzana za pomocą niestandardowego modelu Cellpose, który został wytrenowany i opublikowany w Ref.9. Wytrenowany model Cellpose może być używany bezpośrednio w Cell-ACDC poprzez podanie ścieżki do pliku wag w interfejsie graficznym podczas wyboru parametrów modelu dla Cellpose v2. Segmentacja została wykonana za pomocą drugiego modułu Cell-ACDC do przetwarzania wszystkich obrazów w grupach (Rysunek 1A, moduł "Segment i śledzenie"). Następnie wynik został zwizualizowany w interfejsie GUI trzeciego modułu (Rysunek 1A, moduł "Visualize and correct"). Moduł ten został zoptymalizowany do wizualizacji tysięcy pojedynczych obiektów z wieloma opcjami adnotacji (np. kontury, nakładane maski segmentacyjne lub identyfikatory tekstu). Po obliczeniu pomiarów z masek 3D wszystkie dostępne pomiary były dokumentowane bezpośrednio w interfejsie graficznym. Dostępne są po przejściu do górnego paska menu (rysunek 6, "Menu bar"), wybraniu menu Pomiary , a następnie Ustaw pomiary.... Okno okienkowe pozwala użytkownikom uzyskać informacje specyficzne dla pomiarów z przycisków informacyjnych i wybrać, które pomiary chcą zapisać. Dla reprezentatywnych wyników na Rysunku 10 użyto "cell_vol_fl_3D", czyli objętości każdego obiektu obliczonej przez pomnożenie łącznej liczby wokseli w obiekcie przez kwadrat rozmiaru piksela i głębokość wokseli. Właściwości te są automatycznie wyodrębniane z pliku surowego mikroskopu lub dostarczane przez użytkownika. Obliczenie pomiarów z tego interfejsu wymaga załadowania obrazów RAW. Aby usprawnić proces i umożliwić przetwarzanie wsadowe, pomiary można również obliczyć z menu Utilities (rysunek 1B, "Utilities"), przechodząc do podmenu Pomiary , a następnie obliczaj pomiary dla jednego lub więcej eksperymentów. Na koniec rozkład objętości jądrowej został narysowany jako reprezentatywny wynik (Rysunek 10). Analiza ta ujawnia znaczącą część małych jąder, prawdopodobnie z powodu artefaktów segmentacji. Można je łatwo usunąć, filtrując małe obiekty z maski segmentacji. Jednocześnie bardzo duże jądra mogą być wynikiem połączenia jąder podczas segmentacji. Zawsze zaleca się wykreślanie rozkładu objętości obiektów (np. pojedynczych komórek) w celu identyfikacji artefaktów i uzyskania dodatkowych informacji biologicznych o rozmiarze komórek.
Ilościowa ilościowa analiza danych z mikroskopii poklatkowej
Dzięki mikroskopii poklatkowej dynamika komórkowa może być bezpośrednio obserwowana na poziomie pojedynczych komórek. Poza segmentacją komórek, wyodrębnianie dynamiki czasowej wymaga dodatkowej analizy, w tym śledzenia komórek i adnotacji rodowodu komórek. Ze względu na współzależność tych etapów analizy, błędy wprowadzone na początku potoku mogą przenosić się na kolejne etapy. W związku z tym konieczne jest ciągłe wizualizowanie i korekta błędów segmentacji, śledzenia i adnotacji. Poniżej pokazuje, że Cell-ACDC nadaje się do realizacji tych zadań. Wybrano dwa zestawy danych dwóch różnych organizmów modelowych: 1) drożdże pączkujące (szczep DCY001-1 z Ref.35, gdzie dwa białka histonów H2B, Htb1 i Htb2, są oznaczone mCytrynem) oraz 2) embrionalne komórki macierzyste myszy (mESC, dane z Ref.22).
Te dwa zestawy danych podkreślają dwa tryby adnotacji dostępne w Cell-ACDC: asymetryczny i symetryczny (czyli symetryczna cytokineza, czyli "normalna") podział komórek. Ponieważ sposób podziału się różni, oba organizmy wymagają innego systemu śledzenia i adnotacji.
W podziale asymetrycznym komórka matka tworzy pąk, który rośnie i ostatecznie rozdziela się, tworząc komórkę potomną. Po podziale komórka macierzysta zachowuje swój oryginalny identyfikator komórki, a jej numer generowania wzrasta o jeden, podczas gdy komórka potomna otrzymuje nowy identyfikator komórki, a jej numer generowania jest ustawiany na jeden. Faza pączkowania odpowiada fazom S/G2/M cyklu komórkowego i jest tak oznaczana w Cell-ACDC. Te opcje adnotacji pomagają odpowiadać na typowe pytania biologiczne związane z cyklem komórkowym u drożdży pączkujących.
W przypadku podziału "symetrycznego" (np. komórek ssaków) komórka macierzysta dzieli się na dwie komórki potomne. Komórka macierzysta, czyli jej identyfikacja, znika podczas podziału, a dwie komórki potomne otrzymują nowe identyfikatory. Dodatkowo liczba generacji komórek potomnych jest zwiększona o jeden względem komórki macierzystej. Cell-ACDC śledzi także identyfikator rodzica, identyfikator korzenia (oryginalną komórkę przodka na początku linii rodowej) oraz identyfikator siostrzany.
Dla obu trybów adnotacji opracowano innowacyjne ramy korygujące błędy adnotacji, gdzie poprawka jest automatycznie przekazywana do wszystkich przeszłych i przyszłych istotnych punktów czasowych. Wizualizacja, adnotacje i korekta zostały wykonane w interfejsie GUI trzeciego modułu (Rysunek 1A, moduł "Visualize and correct" oraz Rysunek 6). Aby segmentować i śledzić komórki w zbiorze danych 1, zastosowano model YeaZ_v226 do kanału kontrastu fazowego, natomiast dla kanału jądrowego (histonowego) zastosowano model StarDist25 . Następnie, korzystając z narzędzia Tracking and lineage > Track i/or count sub-cellular objects (Rysunek 1B, menu Utilities ), Cell-ACDC przypisał każdemu jądru identyfikator komórki odpowiadającej komórki, zapewniając spójność między tabelami generowanymi z masek cell a nucleus.
Dla zbioru danych 2 zastosowano model segmentacji DeepSea22 . Wszystkie trzy modele są już dostępne w Cell-ACDC, co pokazuje zaletę integracji kilku modeli segmentacji z oprogramowaniem.
Po poprawieniu błędów segmentacji i śledzenia anonotowano rodowody komórek, obliczano cechy numeryczne oraz przeprowadzano dalszą analizę. Dla zbioru danych 1 kolumna TaYFP_amount_autoBkgr oraz liczba segmentowanych jąder (rysunek 11A) zostały narysowane względem czasu. "TaYFP" to nazwa kanału nuklearnego. "amount_autoBkgr" to wskaźnik pośredni dla całkowitej ilości białka komórkowej wyekstrahowanej z obrazów epifluorescencyjnych36. Oblicza się ją jako różnicę między średnią intensywnością fluorescencji w każdej masce komórki a medianą tła, pomnożoną przez powierzchnię komórki (w pikselach). Tutaj mediana tła jest obliczana ze wszystkich pikseli, które nie są podzielone na komórki. Jak można się spodziewać, ilość H2B zaczyna rosnąć w momencie wyłaniania się pąków (rysunek 11A-ii) i osiąga stałą wartość przed podziałem jądrowym (rysunek 11A-iii). Jest to ważna kontrola jakości homeostazy białek histonowych, ponieważ ilość białek histonowych jest oczekiwana jako zależna od cyklu komórkowego. Dodatkowo, wykres liczby jąder w czasie potwierdza, że ilość białek histonowych osiąga maksimum mniej więcej w okolicach podziału jądrowego.
Dla zbioru danych 2 wykreślono powierzchnię komórki w czasie wybranej komórki podlegającej podziałowi komórki. Zgodnie z oczekiwaniami, powierzchnia komórki rośnie, aż osiągnie maksymalną wartość (rysunek 11B-i). Następnie zmniejsza się aż do podziału komórki (rysunek 11B-ii) wraz ze skurczem komórki. Na koniec cykl zaczyna się od nowa dla dwóch komórek potomnych. To kolejna zalecana analiza, ponieważ sprawdzenie zmian rozmiaru komórek w trakcie cyklu komórkowego jest niezbędne, aby potwierdzić, że komórki rosną i dzielą się zgodnie z oczekiwaniami (lub nie, w przypadku konkretnych mutantów).

Rysunek 1: Moduły Cell-ACDC. (A) Przegląd 4 głównych modułów, które mogą być wystrzeliwane z głównego wyrzutni Cell-ACDC. Po segmentacji, śledzeniu i adnotacji danych mikroskopijnych, cechy numeryczne można obliczyć zarówno z trzeciego modułu ("Wizualizuj i popraw"), jak i z (B) menu Narzędzia na górnym pasku menu. "Utilities" to procedury, które mogą działać automatycznie na wielu zbiorach danych bez udziału użytkownika. Oprócz obliczania pomiarów, inne narzędzia obejmują łączenie wielu tabel wyjściowych w jedną tabelę, śledzenie obiektów subkomórkowych oraz wstępne przetwarzanie obrazów. Cell-ACDC obsługuje dane 2D, 3D (z-stack lub time-lapse) oraz 4D (z-stacki w czasie), z dowolną liczbą dodatkowych kanałów. Tabela wyjściowa z cechami numerycznymi może być następnie wykorzystana do dalszej analizy i odkryć biologicznych (Rysunek 10 i Rysunek 11). W tym celu na stronie Cell-ACDC na GitHub dostępne są notatniki Jupyter, które zawierają przykłady wykresów dostępnych z tabeli wyjściowej. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.

Rysunek 2: Schemat decyzyjny cell-ACDC. Schemat blokowy określający moduł do wykorzystania w zależności od typu zbioru danych i wymagań analitycznych. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.

Rysunek 3: Pobierz przykładowe dane. (A) Otwarty przewodnik powitalny. (B) Pobieranie przykładowych danych potrzebnych do replikacji protokołu. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.

Rysunek 4: Ładowanie danych do przygotowania danych. (A) Załaduj dane do interfejsu GUI przygotowującego dane. (B) Wybierz kanał do załadowania. (C) Edycja i potwierdzenie metadanych obrazu. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.

Rysunek 5: Uruchom proces przygotowania danych. (A) Rozpocznij proces. (B) ROI pozycji (dla przycinania) oraz ROI tła. Prosimy kliknąć tutaj, aby zobaczyć większą wersję tej liczby.

Rysunek 6: Interfejs trzeciego modułu do wizualizacji i poprawy wyników. Zrzut ekranu interfejsu trzeciego modułu (GUI ("Visualise and correct" na Rysunku 1A) z podświetlonymi elementami. Zwróć uwagę, że większość przycisków na paskach ma podpowiedź (dostępną po nałożeniu kursora myszy na przycisk) wyjaśniającą, jak używać tej konkretnej funkcji. Selektor trybu może być używany do przełączania się między 5 trybami: "Viewer", "Segmentation and Tracking", "Cell cycle analysis" (dla komórek dzielących się asymetrycznie), "Normal division: Lineage tree" (dla komórek dzielących się symetrycznie, np. komórek ssaków) oraz "Niestandardowe adnotacje". Należy zauważyć, że pasek narzędzi lub pasek menu często znajduje się w innych interfejsach graficznych (np. moduł "Data pre-processing", Rysunek 1). Paska narzędzi Edycja zawiera wszystkie funkcje, które można użyć do edycji i korygowania błędów segmentacji i śledzenia (np. pędzel, gumka, identyfikator edycji itp.). Załadowany obraz jest wyświetlany w widoku dwóch paneli, co jest pomocne, gdy potrzebne są różne opcje adnotacji (np. informacje o cyklu komórki na lewym obrazie i ID na prawym obrazie). Prawy obraz można również wyłączyć (kliknij prawym przyciskiem myszy na obraz i odznacz Pokaż obraz lustrzany). Każdy panel obrazów ma suwak LUT po boku, który pozwala szybko regulować poziomy intensywności. Klikając prawym przyciskiem myszy na kontrolkę LUT, użytkownik może wybrać różne mapy kolorów dla obrazów intensywności. Dodatkowo, po prawej stronie, znajduje się selektor LUT do nałożenia koloru etykiet segmentacji na obrazy intensywności (opcja adnotacji zwana maskami Segm). Po lewej stronie opcji adnotacji dla lewego obrazu znajdują się dodatkowe przełączniki do sterowania ustawieniami, takimi jak automatyczne zapisywanie, rozmiar czcionki itp. Proszę kliknąć tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 7: Wizualizacja wstępnego przetwarzania w głównym interfejsie graficznym. (A) Otwórz dialog wstępnego przetwarzania. (B) Wstępne przetwarzanie dialogu z parametrami używanymi w inicjalizacji protokołu. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.

Rysunek 8: Wizualizuj wyjście segmentacji w głównym interfejsie graficznym. (A) Wybierz model segmentacji. (B) Ustalić parametry segmentacji. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.

Rysunek 9: Struktura teczek wymagana przez Cell-ACDC. Aby działać z Cell-ACDC, dane muszą być ułożone w określoną strukturę folderów. Chociaż Cell-ACDC udostępnia moduł do automatycznego generowania tej struktury, ważne jest, aby zrozumieć, jak powinna ona wyglądać. Po pierwsze, wszystkie pliki muszą znajdować się w folderze o nazwie Images. Następnie wszystkie muszą zacząć od tej samej nazwy, tzw. "basename_". Minimalnym wymaganym zestawem plików jest jednokanałowy plik TIFF (2D, 3D z-stack lub 3D+time) oraz plik CSV kończący się na "_metadata.csv". Ten plik musi być tabelą z dwiema kolumnami, pierwszą kolumną nazywaną Description , a drugą kolumną o wartościach, i powinien zawierać co najmniej wpisy dla SizeT oraz SizeZ dla liczby klatek i z-slice, odpowiednio. Jeśli plik obrazu nie zawiera wycinków z, rozmiar Z musi być ustawiony na 1. To samo dotyczy zdjęć bez poklatku czasowego, gdzie rozmiar T musi wynosić 1. Jako plik CSV, wpis to pojedyncza linia Description, value, oddzielona przecinkiem, np. SizeZ,1. Dla wielu kanałów musi być wygenerowany jeden plik TIFF na kanał. Folder Images należy następnie umieścić w folderze o nazwie "Position_1". Dozwolone są różne stanowiska, które muszą być oznaczone kolejnymi numerami. Podczas ładowania danych do dowolnego modułu Cell-ACDC użytkownik może wybrać konkretny folder Position lub cały folder eksperymentu. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.

Rysunek 10: 3D ilościowa ilościowa struktura organoidów nowotworowych. Zrzut ekranu reprezentatywnego organoidu nowotworowego (dane z9) załadowanego do interfejsu trzeciego modułu Cell-ACDC (po lewej), przykładowe z-wycinki z czerwonymi konturami podkreślającymi maski segmentacyjne (w środku) oraz histogram rozkładu objętości komórkowej obliczony na podstawie masek segmentacji 3D (po prawej). Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.

Rysunek 11: Ilościowa analiza danych z mikroskopii poklatkowej. (A) Zrzut ekranu z mikroskopii poklatkowej z komórek pączkujących drożdży załadowanych do interfejsu trzeciego modułu Cell-ACDC (po lewej) oraz ilościowa ilościowa ilość białka histonu H2B w czasie w reprezentatywnym cyklu komórkowym (po prawej). Powiększene obrazy pokazują przykładową komórkę i jej pąk (białe strzałki) na początku cyklu komórkowego (i), w momencie wyłaniania się pąków (ii) oraz przy podziale jądrowym (iii). Dane pochodzą z Chatzitheodoridou i in.35 (B) Zrzut ekranu z mikroskopii poklatkowej danych z mikroskopii poklatkowej załadowanych do trzeciego modułu GUI Cell-ACDC (po lewej) oraz powierzchni komórki (μm2) przedstawionej jako funkcja czasu reprezentatywnej komórki podlegającej podziałowi komórkowemu. Powiększone obrazy pokazują przykładową komórkę i jej potomki w maksymalnej powierzchni komórki przed podziałem (i), podziałem na dwie komórki potomne (ii) oraz ostatnią analizowaną klatkę po podziale (iii). Dane pochodzą z Zargari i in. 22,33. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.