Name: Oparte na morfologii rozróżnienie między zdrowymi i patologicznymi komórkami z wykorzystaniem transformaty Fouriera i samoorganizujących się map (Wideo) | JoVE Journal
Uploaded: 2018-10-28T14:00:00.000+00:00
Duration: 8 min 59 s
Description: Poznaj metodę automatycznej identyfikacji komórek zdrowych i patologicznych przy użyciu transformaty Fouriera i samoorganizujących się map. Technika ta wspiera diagnostykę nowotworów i chorób zapalnych poprzez analizę morfologii trójwymiarowej.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

Autorzy dziękują Benjaminowi Krause za owocną dyskusję i wsparcie. Autorzy dziękują również Robertowi Güntherowi za pomoc w mikroskopii żywych komórek.
Praca była wspierana przez DFG wsparcie finansowe NI1167/3-1 (JIMI) dla R.N. i Z.C., DFG wsparcie finansowe CRC 1278 PolyTarget Project Z01 dla Z.C., C01 w TRR130 do R.N. i SFB633, TRR130, Exc257 dla A.E.H. i J.B.S. BfR zapewnił wsparcie stacjonarne SFP1322-642 dla F.L.K i A.L.

1. Wymagania dotyczące protokołu
<ol>
	<li>Aby uzyskać obraz o wysokiej rozdzielczości, należy uzyskać dane z mikroskopii trójwymiarowej (3D) o wysokiej rozdzielczości, zdekonwoluowane zgodnie z kryterium Nyquista, z interwałem próbkowania co najmniej dwukrotnie większym niż najwyższa częstotliwość przestrzenna próbki.</li>
	<li>Użyj oprogramowania do renderowania 3D do rekonstrukcji i eksportu powierzchni.</li>
	<li>Użyj oprogramowania do animacji 3D, które może uruchamiać skrypty Pythona (skrypt Pythona można pobrać z repozytorium github: https://github.com/zcseresn/ShapeAnalysis), aby tworzyć projekcje 2D.</li>
	<li>Użyj Fiji13, aby przeanalizować rzuty 2D i wyodrębnić komponenty DFT.
	<ol>
		<li>Użyj bieżącej dystrybucji Fidżi. Jeśli istnieje już zainstalowana wersja Fidżi, upewnij się, że jest to najnowsza wersja. Można to łatwo osiągnąć, uruchamiając punkt menu Pomoc |Opcja aktualizacji.</li>
		<li>Użyj wtyczki Active Contour15, którą można pobrać z http://imagejdocu.tudor.lu/doku.php?id=plugin:segmentation:active_contour:start i skopiować do folderu wtyczek.</li>
		<li>Pobierz wtyczkę SHADE Fiji z repozytorium github i skopiuj ją do folderu plugins.</li>
	</ol>
	</li>
	<li>Korzystaj z oprogramowania do matematyki obliczeniowej, które jest w stanie obliczać samoorganizujące się mapy.</li>
</ol>
2. Zrekonstruuj obraz 3D.
UWAGA: Do celów testowych, przykładowy zestaw danych znajduje się w repozytorium github (patrz wyżej).
<ol>
	<li>Uruchom oprogramowanie do rekonstrukcji 3D i otwórz dane obrazu 3D.</li>
	<li>Utwórz powierzchnię 3D (wszystkich) obiektów.
	<ol>
		<li>Wybierz opcję Widok 3D i kliknij Powierzchnie. Kliknij przycisk Dalej (niebieskie kółko z białym trójkątem), aby kontynuować pracę z kreatorem tworzenia powierzchni.</li>
		<li>Wybierz kanał obrazu do rekonstrukcji powierzchni.</li>
		<li>Zastosuj funkcję wygładzania, aby uniknąć porowatych powierzchni.
		<ol>
			<li>Wybierz wartość wygładzenia, która nie ukrywa szczegółów powierzchni, ale unika powierzchni porowatych.</li>
		</ol>
		</li>
		<li>Wybierz metodę progowania, aby znaleźć powierzchnie.
		<ol>
			<li>Użyj bezwzględnego progu intensywności, gdy obiekty są dobrze oddzielone od tła i mają w przybliżeniu jednolity poziom jasności.</li>
			<li>Zastosuj próg kontrastu lokalnego, gdy obiekty różnią się intensywnością, ale nadal mogą być oddzielone od lokalnego tła i od innych otaczających je obiektów. Ustaw lokalny obszar wyszukiwania progów zgodnie z wartością oczekiwanej średnicy zrekonstruowanych obiektów.</li>
		</ol>
		</li>
		<li>Przefiltruj zrekonstruowane powierzchnie według interesujących ich parametrów morfologicznych, np. objętości, sferyczności, stosunku powierzchni do objętości itp., a następnie zakończ rekonstrukcję powierzchni.</li>
	</ol>
	</li>
	<li>Zapisz i wyeksportuj wygenerowane powierzchnie w formacie zgodnym z oprogramowaniem do animacji 3D, które będzie używane w następnym kroku.</li>
</ol>
3. Przekształć zrekonstruowane powierzchnie 3D w rzuty 2D
<ol>
	<li>Uruchom Blendera i przejdź do zakładki wyjścia w oknie po prawej stronie. Wybierz format TIFF z menu rozwijanego i ustaw głębię kolorów na 8 bitów RGBA.</li>
	<li>Przełącz się w tryb skryptowy i otwórz dostarczony plik skryptu "GUI_AutoRotate.py" z repozytorium dostarczonego z tą pracą (https://github.com/zcseresn/ShapeAnalysis).</li>
	<li>Kliknij Uruchom skrypt. Wybierz folder plików wrl po wyświetleniu monitu o wprowadzenie.</li>
	<li>W razie potrzeby utwórz więcej obrotów podczas pracy z bardziej złożonymi powierzchniami: przejdź do GUI i ustaw pole Rotacje na wartość powyżej 6, 
	UWAGA: Obrót o 6 różnych kątów może być wystarczający do rozróżnienia różnych populacji komórek. Nie zaleca się tworzenia mniej niż sześciu obrotów na powierzchnię ze względu na potencjalną utratę informacji.</li>
	<li>Uruchom skrypt, klikając przycisk Obróć w graficznym interfejsie użytkownika. Zapisz rzuty poszczególnych powierzchni w tym samym folderze, który został użyty jako folder wejściowy (krok 2.3). Domyślnie obrazy są zapisywane w 8-bitowym formacie Tiff (patrz krok 2.1), który jest formatem wymaganym przez wtyczkę Fiji SHADE.</li>
</ol>
4. Znajdź peryferia i oblicz składowe Fouriera za pomocą Fidżi.
<ol>
	<li>Otwórz Fidżi i wybierz SHADE w menu Wtyczki. Zacznij od wartości domyślnych i dostosuj parametry później. Kliknij przycisk OK, gdy program będzie gotowy do uruchomienia.
	<ol>
		<li>Wybierz wartość Próg gradientu dla progowania obrazu wejściowego.</li>
		<li>Wybierz liczbę iteracji. Im wyższa wartość Liczba iteracji, tym dokładniejsza jest rekonstrukcja peryferii. W przypadku prostszych kształtów zwykle wystarcza niższa liczba.</li>
		<li>Użyj parametru Liczba dylatacji, aby określić, o ile większa jest maska początkowa w porównaniu z rzeczywistą komórką. Zwykle bardziej złożone kształty wymagają większej liczby stopni dylatacji, aby prawidłowo znaleźć obrzeża.</li>
		<li>Zaznacz pole wyboru Ciemne tło, jeśli rzutowane kształty są jaśniejsze niż tło.</li>
		<li>Aktywuj pole wyboru Pokaż wyniki pośrednie tylko wtedy, gdy używasz małego zestawu danych testowych do określenia wydajności funkcji SHADE. Włączenie tej opcji w przypadku większych zestawów danych obniża wydajność obliczeniową i może spowodować zatrzymanie systemu z małą ilością pamięci wideo.</li>
		<li>Zaznacz pole wyboru Zapisz tabele wyników, aby użyć wyników funkcji SHADE jako danych wejściowych w kroku 5. Jeśli to pole jest zaznaczone, wszystkie wyniki są zapisywane w osobnych plikach csv. Podsumowanie danych wyjściowych jest zawsze generowane w pliku o nazwie "Result_collection_of_all_DFT_calculations.csv".</li>
	</ol>
	</li>
	<li>Wybierz folder danych wejściowych, który zawiera pliki TIFF utworzone w kroku 3.</li>
	<li>Podaj folder danych wyjściowych.</li>
	<li>Kliknij OK, aby uruchomić wtyczkę.</li>
</ol>
5. Mapy samoorganizujące się
UWAGA: Sieci SOM są w stanie klasyfikować dane tylko wtedy, gdy są trenowane na dużym zbiorze danych, który zawiera dane wejściowe ze wszystkich oczekiwanych typów komórek i warunków. W celach demonstracyjnych taki zestaw danych jest dostępny i można go znaleźć w naszym repozytorium ("AllCells_summary_normalised.csv" z https://github.com/zcseresn/ShapeAnalysis
<ol>
	<li>Postępuj zgodnie z tymi wskazówkami, jeśli nie ma jeszcze dostępnego przeszkolonego SOM dla danych wejściowych; w przeciwnym razie przejdź do kroku 5.2.
	<ol>
		<li>Uruchom oprogramowanie obliczeniowo-matematyczne zdolne do wykonywania klasyfikacji sieci neuronowych.</li>
		<li>Wybierz plik danych, który ma być używany do trenowania sieci SOM. Ten zestaw danych powinien zawierać wszystkie warunki eksperymentalne w celu wytrenowania SOM na określonych typach komórek i warunkach eksperymentalnych. 
		UWAGA: Możliwe jest również użycie dostarczonego AllCells_summary_normalised.csv do testowania systemu.</li>
		<li>Rozpocznij trening i poczekaj, aż trening zostanie zakończony, zanim przejdziesz dalej. Domyślnie skrypt jest ustawiony na uruchamianie 2000 iteracji ("Epoki"). 
		UWAGA: Liczba iteracji zależy od szybkości uczenia się SOM. W zależności od danych wejściowych wskazane jest przetestowanie zarówno wyższej, jak i dolnej liczby epok oraz obserwacja stabilności wzorca SOM. W przypadku korzystania z dostarczonego skryptu liczbę iteracji można zmienić w wierszu 32. Rozmiar sieci można zmienić w wierszu 34 (domyślnie jest ustawiony na 12 na 12).</li>
		<li>Po zakończeniu szkolenia należy sprawdzić topologię sieci (odległości sąsiadów, płaszczyzny wejściowe, trafienia próbki itp.). Sieć jest teraz przeszkolona i można ją zapisać do wykorzystania w przyszłości.</li>
	</ol>
	</li>
	<li>Załaduj SOM w przypadku korzystania z już wytrenowanej mapy (może to pochodzić z kroku 5.1 lub z innych źródeł) w celu klastrowania zestawu danych.
	<ol>
		<li>Zaimportuj plik csv, który ma być testowany ze wstępnie załadowanym wytrenowanym SOM. Wybierz wyjście csv wtyczki SHADE z kroku 4 w przypadku korzystania z danych przygotowanych przez wtyczkę SHADE. 
		UWAGA: Możliwe jest również użycie przykładowych plików danych "InteractingCells_summary_normalised.csv", "MobileCells_summary_normalised.csv" lub "PhagocytosingCells_summary_normalised.csv", które są dostarczane za pośrednictwem github.</li>
		<li>Po zakończeniu klasyfikacji oceń wyniki SOM jak w kroku 5.1.5.
		<ol>
			<li>Sprawdź mapę trafień wygenerowaną na podstawie pliku csv. Każda komórka mapy pokazuje, ile razy zestaw danych "trafia" w tę konkretną komórkę wytrenowanego SOM. Gdy grupa komórek jest zgrupowana na małym obszarze tej mapy, oznacza to, że zestaw danych jest dość jednorodny. Wiele klastrów będzie oznaczać, że w zestawie danych prawdopodobnie istnieją podgrupy.</li>
			<li>Zbadaj odległości wagowe w okolicy. Obszary tej mapy, które są dobrze oddzielone, odpowiadają grupom obiektów, które zachowują się zupełnie inaczej z punktu widzenia SOM. W przypadku komponentów DFT jako danych wejściowych oznacza to, że te grupy komórek mają bardzo różne kształty odpowiadających im powierzchni 3D.</li>
			<li>Sprawdź płaszczyzny wag, aby uzyskać informacje o udziale każdego elementu w wektorze cech. W przypadku użycia 20 komponentów DFT, jak opisano wcześniej, pojawi się tutaj 19 map. W przypadku korzystania z podanego przykładowego zestawu danych pierwsze 5 lub 6 płaszczyzn wagi będzie się różnić, ale reszta z nich będzie wyglądać dość podobnie. W tym przypadku można stwierdzić, że wystarczyłoby użycie około 7 komponentów DFT.</li>
		</ol>
		</li>
	</ol>
	</li>
</ol>

<ol><li>Masedunskas, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Intravital+microscopy%3A+a+practical+guide+on+imaging+intracellular+structures+in+live+animals%5BTitle%5D)+AND+%22Bioarchitecture%22%5BJournal%5D)">Intravital microscopy: a practical guide on imaging intracellular structures in live animals.</a> Bioarchitecture. 2 (5), 143-157 (2012).</li>
<li>Niesner, R. A., Hauser, A. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Recent+advances+in+dynamic+intravital+multi-photon+microscopy%5BTitle%5D)+AND+%22Cytometry+A%22%5BJournal%5D)">Recent advances in dynamic intravital multi-photon microscopy.</a> Cytometry A. 79 (10), 789-798 (2011).</li>
<li>Ho, S. Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((NeurphologyJ%3A+an+automatic+neuronal+morphology+quantification+method+and+its+application+in+pharmacological+discovery%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">NeurphologyJ: an automatic neuronal morphology quantification method and its application in pharmacological discovery.</a> BMC Bioinformatics. 12, 230(2011).</li>
<li>Yin, Z., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+screen+for+morphological+complexity+identifies+regulators+of+switch-like+transitions+between+discrete+cell+shapes%5BTitle%5D)+AND+%22Nature+Cell+Biology%22%5BJournal%5D)">A screen for morphological complexity identifies regulators of switch-like transitions between discrete cell shapes.</a> Nature Cell Biology. 15 (7), 860(2013).</li>
<li>Yu, H. Y., Lim, K. P., Xiong, S. J., Tan, L. P., Shim, W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Functional+Morphometric+Analysis+in+Cellular+Behaviors%3A+Shape+and+Size+Matter%5BTitle%5D)+AND+%22Advanced+Healthcare+Materials%22%5BJournal%5D)">Functional Morphometric Analysis in Cellular Behaviors: Shape and Size Matter.</a> Advanced Healthcare Materials. 2 (9), (2013).</li>
<li>Johnson, G. R., Buck, T. E., Sullivan, D. P., Rohde, G. K., Murphy, R. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Joint+modeling+of+cell+and+nuclear+shape+variation%5BTitle%5D)+AND+%22Molecular+Biology+of+the+Cell%22%5BJournal%5D)">Joint modeling of cell and nuclear shape variation.</a> Molecular Biology of the Cell. 26 (22), 4046-4056 (2015).</li>
<li>Wang, S. -H., Cheng, H., Phillips, P., Zhang, Y. -D. <a target="_blank" href="http://www.mdpi.com/1099-4300/20/4/254">Multiple Sclerosis Identification Based on Fractional Fourier Entropy and a Modified Jaya Algorithm.</a> Entropy. 20 (4), 254(2018).</li>
<li>Kriegel, F. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Cell+shape+characterization+and+classification+with+discrete+Fourier+transforms+and+self-organizing+maps%5BTitle%5D)+AND+%22Cytometry+Part+A%22%5BJournal%5D)">Cell shape characterization and classification with discrete Fourier transforms and self-organizing maps.</a> Cytometry Part A. 93 (3), 323-333 (2017).</li>
<li>Styner, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Framework+for+the+Statistical+Shape+Analysis+of+Brain+Structures+using+SPHARM-PDM%5BTitle%5D)+AND+%22Insight+Journal%22%5BJournal%5D)">Framework for the Statistical Shape Analysis of Brain Structures using SPHARM-PDM.</a> Insight Journal. (1071), 242-250 (2006).</li>
<li>El-Baz, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((3D+shape+analysis+for+early+diagnosis+of+malignant+lung+nodules%5BTitle%5D)+AND+%22Medical+Image+Computing+and+Computer+Assisted+Intervention%22%5BJournal%5D)">3D shape analysis for early diagnosis of malignant lung nodules.</a> Medical Image Computing and Computer Assisted Intervention. 14 (Pt 3), 175-182 (2011).</li>
<li>Williams, E. L., El-Baz, A., Nitzken, M., Switala, A. E., Casanova, M. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Spherical+harmonic+analysis+of+cortical+complexity+in+autism+and+dyslexia%5BTitle%5D)+AND+%22Translational+Neuroscience%22%5BJournal%5D)">Spherical harmonic analysis of cortical complexity in autism and dyslexia.</a> Translational Neuroscience. 3 (1), 36-40 (2012).</li>
<li>Kruggel, F. <a target="_blank" href="https://www.sciencedirect.com/science/article/pii/S1361841507001107">Robust parametrization of brain surface meshes.</a> Medical Image Analysis. 12 (3), 291-299 (2008).</li>
<li>Schindelin, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Fiji%3A+an+open-source+platform+for+biological-image+analysis%5BTitle%5D)+AND+%22Nature+Methods%22%5BJournal%5D)">Fiji: an open-source platform for biological-image analysis.</a> Nature Methods. 9 (7), 676-682 (2012).</li>
<li>Kohonen, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Essentials+of+the+self-organizing+map%5BTitle%5D)+AND+%22Neural+Networks%22%5BJournal%5D)">Essentials of the self-organizing map.</a> Neural Networks. 37, 52-65 (2013).</li>
<li>Andrey, P., Boudier, T. <a target="_blank" href="https://pdfs.semanticscholar.org/4b2d/57e828777e7f5578b1192fd804890f9a464f.pdf">Adaptive Active Contours.</a> ImageJ user and developer conference. , Luxembourg. (2006).</li>
<li>Burger, W., Burge, M. J. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aBurger%2C+W+%22Principles+of+Digital+Image+Processing%22">Principles of Digital Image Processing</a>. , Springer-Verlag. (2013).</li>
<li>Lestrel, P. E. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aLestrel%2C+PE+%22Fourier+Descriptors+and+their+Applications+in+Biology%22">Fourier Descriptors and their Applications in Biology</a>. , Cambridge University Press. (2008).</li>
<li>Bayerl, S. H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Time+lapse+in+vivo+microscopy+reveals+distinct+dynamics+of+microglia-tumor+environment+interactions-a+new+role+for+the+tumor+perivascular+space+as+highway+for+trafficking+microglia%5BTitle%5D)+AND+%22Glia%22%5BJournal%5D)">Time lapse in vivo microscopy reveals distinct dynamics of microglia-tumor environment interactions-a new role for the tumor perivascular space as highway for trafficking microglia.</a> Glia. 64 (7), 1210-1226 (2016).</li>
</ol>

Autorzy oświadczają, że nie mają konkurencyjnych interesów finansowych.

Identyfikacja potencjalnie patologicznych stanów przy użyciu małych, nienaruszonych próbek tkanek ma duże znaczenie. Takie techniki zapewnią szybką reakcję na choroby zakaźne i agresywne rodzaje raka. Odpowiedzi kinetyczne i morfologiczne różnych komórek odpornościowych, np. mikrogleju i makrofagów, są charakterystyczne dla odpowiedzi immunologicznej organizmu. Chociaż w większości przypadków monitorowanie kinetycznego zachowania tych komórek nie jest praktyczne ani nawet możliwe, dość łatwo jest uzyskać trójwymiarowe obrazy w celu odzyskania ich kształtu. Zazwyczaj komórki odpornościowe przyjmują złożony kształt w zdrowej tkance i znacznie prostszą formę w stanach zapalnych lub nowotworowych18. Podczas gdy zależne od czasu cechy takiej zmiany kształtu wzbogaciłyby naszą wiedzę na temat rozwoju odpowiedzi immunologicznej, użycie tylko kształtu 3D reprezentatywnej grupy komórek może być również wystarczające do określenia zdrowej lub patologicznej natury tkanki.
Scharakteryzowanie trójwymiarowej powierzchni komórki nie jest prostym zadaniem. Zastosowanie harmonicznych sferycznych jest sposobem na przedstawienie powierzchni 3D ze stosunkowo dużą liczbą (50-70) składników11,12. Ponadto wyznaczanie harmonicznych sferycznych jest kosztowne obliczeniowo; rzutowanie bardzo skomplikowanych kształtów na sferę jednostkową jest albo niemożliwe, albo bardzo trudne ze względu na konieczność nałożenia na sferę jednostkową wielu siatek o różnym stopniu rozdrobnienia; Wreszcie, sensowna interpretacja widm składowych harmonicznych sferycznych jest daleka od trywialnej.
W naszej pracy przedstawionej tutaj zastępujemy trudne zadanie bezpośredniej analizy powierzchni 3D znacznie prostszym podejściem polegającym na wykorzystaniu rzutów 2D oryginalnej powierzchni w celu uzyskania wystarczających informacji morfologicznych do identyfikacji stanów patologicznych. Zademonstrowaliśmy każdy etap tego procesu pracy, wykorzystując dane mikroskopowe 3D z komórek szpikowych, jednocześnie wyraźnie zaznaczając, że wszystkie kroki były proste do wykonania, a uzyskane dwuwymiarowe mapy były łatwe do interpretacji.
Oczywiście projekcja 3D do 2D doprowadzi do utraty informacji o strukturze powierzchni. W naszym przykładowym zestawie danych mikrogleju w mysim modelu guza korowego wystarczyło użyć sześciu kątów podczas tworzenia projekcji 2D. Jednak bardziej złożone kształty lub mniej widoczne zmiany morfologiczne mogą wymagać utworzenia większej liczby projekcji, aby móc wiarygodnie zidentyfikować podgrupy komórek za pomocą SOM. Z tego powodu nasze podejście zostało zaprojektowane tak, aby móc generować i analizować dowolną liczbę projekcji. Po prostu wybierając większą liczbę projekcji dla bardziej złożonych kształtów, możliwe jest skalowanie utraty informacji do tolerowanego minimum. Na przykład typ współpracującej komórki na rysunkach 4a i 4b wymagałby większej liczby rzutów, aby prawidłowo odwzorować złożoną powierzchnię.
Jak każda metoda przybliżona, proponowany w niniejszym sposób pracy musiał zostać przetestowany w porównaniu z wynikami ręcznego procesu klasyfikacji mikrogleju18. Zaprezentowane wcześniej wyniki potwierdziły niezawodność zautomatyzowanego przepływu pracy. Co więcej, przepływ pracy jest bardziej efektywny czasowo w porównaniu z konwencjonalną analizą. Ekspert medyczny, który ręcznie sklasyfikował komórki mikrogleju, potrzebował około 4 tygodni na analizę zestawu danych, podczas gdy nasz przepływ pracy potrzebował tylko około 1 dnia. Solidność naszego podejścia została również wyraźnie potwierdzona przez odtwarzalność wytrenowanego SOM do podzbioru danych, które należały do tego samego typu komórki, ale nie zostały wykorzystane do trenowania SOM, jak pokazano na rysunku 3c.
Mimo że nasze podejście nie uwzględniało informacji kinetycznych, zbadaliśmy wpływ czasu na analizę kształtu opartą na DFT. Najbardziej typowy przykład zachowania zależnego od czasu znaleziono w populacji komórek mobilnych, gdzie wkład wyżej indeksowanych komponentów DFT był wyraźnie widoczny, jak na rysunku 4a. Zwraca to uwagę na znaczenie wykorzystania wystarczająco dużej liczby komponentów DFT w przypadku typów komórek, które mogą zachowywać się w sposób bardzo zależny od czasu. Ze względu na zautomatyzowany charakter i wysoką szybkość wykonywania naszych narzędzi programowych, zwiększona liczba komponentów i prognoz DFT zwiększy precyzję i wiarygodność wyników, a jednocześnie nie ograniczy znacząco wydajności obliczeniowej.

Terminowe, proste i precyzyjne określenie stanu patologicznego tkanki biologicznej jest przedmiotem największego zainteresowania w badaniach biomedycznych. Modele mysie dostarczają środków do badania szeregu stanów patologicznych, takich jak reakcje immunologiczne lub rozwój raka, w połączeniu ze złożonymi technikami mikroskopii 3D i 4D (3 wymiary przestrzenne i czas). Badania mikroskopowe można wykonywać za pomocą mikroskopii 2-fotonowej przyżyciowej lub tkankowej z wyciętymi tkankami, mikroskopii arkuszowej światła i - do ograniczonej głębokości tkanki około 100 μm - za pomocą mikroskopii konfokalnej. Aby uzyskać informacje związane z czasem na temat zachowania komórek w warunkach fizjologicznych lub patologicznych, konieczne jest monitorowanie tkanki przez dłuższy czas, co zwykle wymaga obrazowania przyżyciowego1,2. Oczywiście zastosowanie tej techniki jest ograniczone do modeli zwierzęcych ze względu na jej inwazyjność. Dostępne są również techniki nieinwazyjne do zastosowań u ludzi, w tym różne metody tomografii (MSOT, CT itp.), ale wszystkie te metody nie mają niezbędnej rozdzielczości przestrzennej - i często czasowej - do badania zachowania na poziomie komórkowym.
Statyczne informacje dotyczące wyglądu komórek mogą być łatwiej dostępne dzięki różnym technikom obrazowania 3D wykonanym na wyciętych próbkach tkanek. W tym przypadku nie mierzy się zachowania kinetycznego komórek, dlatego konieczne jest zastosowanie nowatorskich technik analizy, które są w stanie określić stan patogenny badanych komórek wyłącznie na podstawie ich morfologii3. Takie podejście wykorzystano do powiązania kształtów komórek i tekstur tkanek z patologicznym zachowaniem4,5,6.
W nowej technice opisanej tutaj, komórki są rekonstruowane jako powierzchnie 3D, a ich kształty są charakteryzowane za pomocą projekcji 3D do 2D i kolejnej analizy kształtów peryferii opartej na Fourierze7,8. Zmniejszając wymiary z 3 do 2, problem jest uproszczony. Możliwe jest również scharakteryzowanie powierzchni komórek w 3D poprzez zastosowanie analizy harmonicznych sferycznych, tak jak to zostało zrobione w przypadku obrazów medycznych9. Jednak harmoniczne sferyczne nie radzą sobie dobrze z ostrymi i chropowatymi kształtami, co wymaga ustanowienia wieloskalowej siatki na sferze jednostkowej. Ponadto liczba niezbędnych składowych harmonicznych sferycznych może być duża (50-70), przy czym obliczenia leżące u ich podstaw są bardzo wymagające, a wyniki trudne do interpretacji10,11,12.
Dzięki naszej nowo zaproponowanej metodzie, zadanie jest zredukowane do serii opisów kształtów 2D, gdzie liczba projekcji 2D zależy od analityka i może być dostosowana do złożoności kształtu 3D. Projekcje są generowane automatycznie za pomocą skryptu Pythona, który działa w narzędziu do animacji 3D. Rzuty 2D są opisane przez dyskretne składowe transformaty Fouriera (DFT) ich peryferii, obliczone przez wtyczkę Fiji13, która jest tutaj dostarczana jako część naszego pakietu oprogramowania. DFT jest tutaj stosowany w celu rozłożenia złożonego konturu komórki na szereg funkcji grzechu i cos. W ten sposób możemy opisać konspekt za pomocą stosunkowo niewielkiej liczby składowych DFT, zmniejszając w ten sposób złożoność problemu (więcej szczegółów w sekcji Równania). Komponenty DFT są umieszczane w wytrenowanej mapie samoorganizującej się (SOM14), gdzie można obiektywnie przetestować istnienie klastrów kształtów8. SOM zapewniają konkurencyjne i nienadzorowane narzędzie uczenia się z dziedziny sztucznej inteligencji. Składają się one z połączonego układu sztucznych neuronów, które komunikują się ze sobą za pomocą ważonej funkcji odległości sąsiedztwa. Układ neuronalny reaguje na pierwszy element wejściowego zbioru danych, a neurony, których odpowiedź jest najsilniejsza, są "zgrupowane" bliżej siebie. W miarę jak układ neuronowy otrzymuje coraz więcej danych wejściowych, neurony danych, które wielokrotnie silnie reagują, zaczynają tworzyć dobrze zdefiniowane klastry w systemie. Po odpowiednim wytrenowaniu na dużym zbiorze danych, który zawiera informacje o kształcie 2D w postaci zestawu komponentów DFT, dowolne komponenty DFT pojedynczej komórki można umieścić w wytrenowanym SOM i ujawnić, czy komórka prawdopodobnie należy do zdrowej, czy patogennej grupy komórek. Oczekujemy, że takie narzędzie stanie się doskonałym uzupełnieniem metod diagnostyki naukowej i klinicznej.

<table><tbody><tr><td>Imaris 9.1.2, oprogramowanie</td><td>Bitplane, Z&uuml; bogata, Szwajcaria</td><td>v.9.1.2</td><td>Rekonstrukcja obrazu 3D i generowanie powierzchni; to było przez nas używane!</td></tr><tr><td>Blender 2.75a, oprogramowanie</td><td>https://www.blender.org/</td> <td>v.2.75a</td><td>Oprogramowanie do animacji 3D i 4D o otwartym kodzie źródłowym; 2.75a to wymagana wersja dla tego Pythona</td></tr><tr><td>Fiji / ImageJ, oprogramowanie</td><td>https://fiji.sc/</td> <td>ImageJ v.1.52b</td><td>Zestaw narzędzi do analizy obrazów multi-D o otwartym kodzie źródłowym</td></tr><tr><td>MATLAB</td><td>MathWorks, www.mathworks.com</td><td>R2017b</td><td>Ogólne oprogramowanie obliczeniowe matematyczne</td></tr><tr><td>Maszyna MATLAB Zestaw do nauki</td><td>MathWorks, www.mathworks.com</td><td>R2017b</td><td>Może być używany tylko razem z</td></tr><tr><td>wtyczkami MATLAB Fiji: SHADE</td><td></td><td>https://github.com/zcseresn/ShapeAnalysis v.1.0</td><td></td></tr><tr><td>Wtyczki Fidżi: ActiveContour</td><td>http://imagejdocu.tudor.lu/doku.php? id=plugin:segmentation:active_contour:start</td><td>absnake2</td><td></td></tr><tr><td></td><td>Komputer Dowolny</td><td>NA</td><td>Zapoznaj się z instrukcjami Imaris, aby uzyskać informacje o minimalnych wymaganiach dotyczących komputera</td></tr></tbody></table>

oparte na morfologii rozróżnienie między zdrowymi i patologicznymi komórkami z wykorzystaniem transformaty fouriera i samoorganizujących się map

Morphology-Based Distinction Between Healthy and Pathological Cells Utilizing Fourier Transforms and Self-Organizing Maps

Wygląd i ruchy komórek odpornościowych są zależne od ich środowiska. W reakcji na inwazję patogenu komórki odpornościowe są rekrutowane do miejsca zapalenia i są aktywowane, aby zapobiec dalszemu rozprzestrzenianiu się inwazji. Znajduje to również odzwierciedlenie w zmianach w zachowaniu i wyglądzie morfologicznym komórek odpornościowych. W tkance nowotworowej podobne zmiany morfokinetyczne zaobserwowano w zachowaniu komórek mikrogleju: mikroglej wewnątrznowotworowy ma mniej złożone 3-wymiarowe kształty, ma mniej rozgałęzione procesy komórkowe i porusza się szybciej niż te w zdrowej tkance. Badanie takich właściwości morfokinetycznych wymaga złożonych technik mikroskopii 3D, które mogą być niezwykle trudne, gdy są wykonywane wzdłużnie. W związku z tym rejestracja statycznego kształtu 3D komórki jest znacznie prostsza, ponieważ nie wymaga pomiarów przyżyciowych i może być wykonywana również na wyciętej tkance. Niezbędne jest jednak posiadanie narzędzi analitycznych, które pozwalają na szybki i precyzyjny opis kształtów 3D oraz umożliwiają diagnostyczną klasyfikację zdrowych i patogennych próbek tkanek wyłącznie na podstawie statycznych informacji związanych z kształtem. W tym miejscu przedstawiamy zestaw narzędzi, który analizuje dyskretne składowe Fouriera konturu zestawu projekcji 2D powierzchni komórek 3D za pomocą map samoorganizujących się. Zastosowanie metod sztucznej inteligencji pozwala naszej strukturze uczyć się o różnych kształtach komórek w miarę stosowania jej do coraz większej liczby próbek tkanek, podczas gdy przepływ pracy pozostaje prosty.

Zastosowaliśmy dyskretną transformatę Fouriera (DFT) w celu obliczenia głównych składowych kształtu odpowiadających projekcjom komórek. Deskryptory Fouriera uzyskano poprzez zastosowanie algorytmu DFT do par współrzędnych xy dopasowanego obrysu projekcji komórek, otrzymanych jako wynik działania modułu AbSnake w naszym schemacie analizy. Pary współrzędnych xy te można traktować jako dwuwymiarowy wektor g o wartościach zespolonych: 
<img alt="Ciąg równań, notacja matematyczna g=(g₀,g₁,,...,gₘ₋₁), analiza ciągów." src="/files/ftp_upload/58543/58543eq1.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 1">
Z wektora „g” wykorzystujemy DFT do obliczenia zespolonego widma Fouriera: 
<img alt="Wzór teorii grafów, G=(G0,G1,...,GM-1), reprezentacja sekwencyjna, równania matematyczne." src="/files/ftp_upload/58543/58543eq2.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 2"> 
W oparciu o powszechnie znane wzory na dyskretne widmo Fouriera oraz stosując zapis liczby zespolonej dla „g” jako: 
<img alt="Reprezentacja liczb zespolonych; wzór gₖ=xₖ+i·yₖ; równanie matematyczne." src="/files/ftp_upload/58543/58543eq3.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 3"> 
Otrzymujemy: 
<img alt="Wzór transformaty Fouriera, analiza zespolona, równanie matematyczne, koncepcja edukacyjna." src="/files/ftp_upload/58543/58543eq4.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 4"> (1) 
Możemy obliczyć część rzeczywistą („A”) i urojoną („B”) wartości <img alt="Równanie wzmocnienia zespolonego dla układów optycznych: Gm=(Am+i·Bm)." src="/files/ftp_upload/58543/58543eq5.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 5">: 
<img alt="Wzór analizy Fouriera, Σ dla dekompozycji sygnału, ilustracja matematyczna." src="/files/ftp_upload/58543/58543eq6.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 6"> (2) 
<img alt="Wzór równania równowagi statycznej, Bm równa się Im(Gm), reprezentacja szeregów Fouriera, schemat analizy matematycznej." src="/files/ftp_upload/58543/58543eq7.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 7"> (3) 
Tutaj pierwszy składnik DFT G0 odpowiada wartości m = 0, co daje: 
<img alt="Równanie równowagi statycznej, wzór A0=(1/M)Σ[xk*cos(θ)+yk*sin(θ)], koncepcja matematyczna." src="/files/ftp_upload/58543/58543eq8.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 8"> (4) 
<img alt="Równanie równowagi statycznej, Σ[y*cos(0)-z*sin(0)]=ȳ, analiza wzoru matematycznego." src="/files/ftp_upload/58543/58543eq9.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 9"> (5) 
W konsekwencji komponent ten opisuje środek geometryczny oryginalnego obiektu. 
Drugi element widma DFT w kierunku prostym, G1odpowiada m = 1: 
<img alt="Równanie transformaty Fouriera, schemat, G_k=G_1*e^(i2πk/M), koncepcja przetwarzania sygnałów." src="/files/ftp_upload/58543/58543eq10.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 10"> 
<img alt="Równanie szeregu Fouriera z komponentami trygonometrycznymi; wizualizacja wzoru." src="/files/ftp_upload/58543/58543eq11.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 1"> (6)
Z równania 6 wynika, że punkty te tworzą okrąg o promieniu <img alt="r₁=|G₁| równanie, wzór matematyczny, obliczanie wartości bezwzględnej, analiza, pojęcia algebraiczne" src="/files/ftp_upload/58543/58543eq12.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 12"> oraz kąt początkowy <img alt="Wzór trygonometryczny θ=tan⁻¹(B₁/A₁) do obliczania kąta w analizie matematycznej." src="/files/ftp_upload/58543/58543eq13.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 13">, gdzie okrąg opisuje jeden pełny obrót, podczas gdy kształt jest kreślony raz. Środek okręgu znajduje się w początku układu współrzędnych (0, 0), a promień wynosi |G1| a punktem wyjścia jest:
<img alt="Równania równowagi statycznej, symbol transformaty Fouriera, schemat matematyczny analizy widmowej." src="/files/ftp_upload/58543/58543eq14.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 14"> <img alt="równowaga statyczna; schemat ΣFx=0; bilans sił; analiza układu mechanicznego; narzędzie edukacyjne" src="/files/ftp_upload/58543/58543eq15.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 15"> (7)
Ogólnie, dla pojedynczego współczynnika Fouriera <img alt="Wzór matematyczny dla funkcji Gm jako pary Am, Bm w badaniu teoretycznym." src="/files/ftp_upload/58543/58543eq16.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 16">współrzędne opisano następująco:
<img alt="Równanie matematyczne, reprezentacja transformaty Fouriera, złożona funkcja wykładnicza, schemat." src="/files/ftp_upload/58543/58543eq17.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 17"> 
<img alt="Równanie transformaty Fouriera; zespolone wyrażenia wykładnicze w cosinusie i sinusie; analiza matematyczna." src="/files/ftp_upload/58543/58543eq18.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 18"> (8)
Podobnie jak równanie 6, równanie 8 również opisuje okrąg, lecz o promieniu Rm=|Gm|, kąt początkowy <img alt="Równanie funkcji arkus tangens, tan⁻¹(Bₘ/Aₘ), w wyrażeniu matematycznym." src="/files/ftp_upload/58543/58543eq19.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 19"> or punkt wyjścia w <img alt="Równanie: g₀^(m) = Gₘ; reprezentacja symboliczna; notacja matematyczna; zastosowanie w badaniach." src="/files/ftp_upload/58543/58543eq20.jpg" data-altupdatedat="1788250015000" data-alt="Równanie 20">, gdzie kontur jest obrysowany raz, podczas gdy okrąg wykonuje „m” pełnych obrotów16,17.
Parametry kształtu jako dane wejściowe SOM 
Przebieg procesu, opisany na Rysunku 1, zastosowano do zestawu danych z intravitalnej mikroskopii wielofotonowej komórek mikrogleju, poddanych dekonwolucji (z wykorzystaniem zmierzonej funkcji rozmycia punktu – Point Spread Function), aby scharakteryzować ich zmiany morfologiczne w zdrowej lub nowotworowej tkance korowej18. Dla każdej projekcji 2D zrekonstruowanych powierzchni 3D obliczono dwudziestą liczbę komponentów DFT, a wyniki posłużyły jako dane wejściowe do trenowania SOM. W warunkach fizjologicznych mikroglej wykazywał dość złożony kształt z licznymi, silnie rozgałęzionymi wypustkami (Rysunek 2a). W środowisku nowotworowym (model guza kory mózgowej) mikroglej zmienił się w prostszy kształt, bardziej przypominający wrzeciono (Rysunek 2b).
Przeszkolona SOM została przetestowana w celu oceny jej zdolności do rozróżniania komórek zdrowych od nowotworowych. Populacja komórek zdrowych została rzutowana na pojedynczy obszar SOM (Rycina 2c). SOM zareagowała na zestaw danych mikrogleju nowotworowego aktywnym obszarem w kształcie hantla (Rycina 2d). Ślepo zmieszany zestaw danych wejściowych, składający się z komponentów kształtu DFT zarówno z grupy zdrowej, jak i nowotworowej, został przez SOM rzutowany na dwie odrębne grupy, przy jednoczesnym zachowaniu kształtu ich poszczególnych konturów podobnego do tych z grup oddzielnych (Rycina 2e; porównać z 2c oraz 2d). Można wywnioskować, że zmieszany zestaw danych został pomyślnie zgrupowany przez SOM.
Wydajność SOM przetestowano, porównując jej projekcje z manualną analizą tych samych danych przeprowadzoną przez eksperta medycznego, który zaklasyfikował zbiór danych na podstawie ich zachowania czasoprzestrzennego. Ekspert zidentyfikował cztery odrębne grupy komórek (komórki spoczynkowe, komórki fagocytujące, komórki oddziałujące oraz komórki ruchome18), które zrekonstruowano i wykorzystano do wytrenowania sieci SOM o rozmiarze 12x12. Wytrenowana sieć (Rysunek 3a) wykazuje grupy sztucznych neuronów o wysokich wartościach trafień, szczególnie w dolnej lewej oraz środkowej części SOM. Odpowiedź wytrenowanej sieci przetestowano również przy użyciu czterech losowo wybranych podzbiorów obrazów (które nie stanowiły części zbioru treningowego), pochodzących z czterech różnych grup zidentyfikowanych przez eksperta18. Podzbiory te doprowadziły do uzyskania czterech dobrze zdefiniowanych odpowiedzi SOM, co przedstawiono na Rysunku 3b. Komórki spoczynkowe wykazują najbardziej złożony kształt i najwyższy poziom separacji w sieci neuronowej (Rysunek 3b, panel „resting”). Pozostałe trzy zidentyfikowane typy komórek dzieliły wspólny obszar SOM w dolnym lewym rogu, lecz w pozostałym zakresie zostały rozdzielone przez SOM. Zatem obszar SOM w dolnym lewym rogu odpowiada wartościom DFT o niższych indeksach.
Wytrzymałość podejścia SOM została przetestowana poprzez zastosowanie wytrenowanej mapy SOM do trzech losowych podzbiorów tego samego typu komórek w stanie spoczynku (które nie stanowiły części zbioru treningowego). Odpowiedź SOM na te dane wejściowe jest bardzo podobna (Rysunek 3c, podzbiory 1-3), co potwierdza odporność naszego podejścia.
Zależne od czasu zmiany kształtu komórek są precyzyjnie charakteryzowane za pomocą DFT 
W celu zbadania wpływu zależnych od czasu zmian kształtu komórek na komponenty DFT, od jednej do trzech komórek z podgrupy (patrz Rycina 3b) śledzono przez 13 do 28 punktów czasowych. Rycina 4 przedstawia pierwsze dziesięć komponentów DFT komórki mobilnej (Rycina 4a) oraz komórki wchodzącej w interakcje (Rycina 4b), wykreślone jako funkcja czasu. Komórka mobilna wykazuje stale zmieniający się kształt (patrz Supplementary Video 4 w 8), co odzwierciedla bardziej chropowata powierzchnia DFT. Skoki amplitudy DFT w pierwszej jednej trzeciej przebiegu czasowego dla komórki wchodzącej w interakcje pokrywają się z szybkimi i znacznymi zmianami kształtu komórki, co pokazano w Supplementary Video 5 w 8.
Przebieg czasowy wszystkich 19 komponentów DFT został również scharakteryzowany dla tych dwóch komórek w trzech oddzielnych punktach czasowych podczas śledzenia komórki mobilnej (Rysunek 5a) oraz komórki wchodzącej w interakcje (Rysunek 5b). Osie prostopadłe reprezentują tutaj sześć kątów obrotu i wskazują, że wszystkie projekcje są równie istotne dla charakterystyki kształtu obu typów komórek.
<img alt="Schemat przetwarzania danych mikroskopowych: rekonstrukcja 3D, projekcje 2D, DFT, analiza SOM, klasteryzacja." class="xfigimg" src="/files/ftp_upload/58543/58543fig1.jpg" data-altupdatedat="1788250015000" data-alt="Rycina 1"> 
Rysunek 1. Krok po kroku przedstawiony schemat przetwarzania danych w celu identyfikacji klastrów komórkowych na podstawie kształtu komórek. Powierzchnie zrekonstruowane w 3D posłużyły jako dane wejściowe do programu Blender w celu automatycznego wykonania projekcji 3D na 2D. Wyznaczono obrys każdej projekcji, a następnie obliczono komponenty DFT. Komponenty te stanowiły dane wejściowe dla wytrenowanej mapy samoorganizującej (SOM) w programie Matlab lub posłużyły do wytrenowania nowej mapy SOM. <a href="https://www.jove.com/files/ftp_upload/58543/58543fig1large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.</a>
<img alt="Zielone struktury molekularne w mikrografiach 3D wraz z diagramami do analizy kształtu w skali mikrometrycznej." class="xfigimg" src="/files/ftp_upload/58543/58543fig2.jpg" data-altupdatedat="1788250015000" data-alt="Rycina 2"> 
Rysunek 2. Typowy wygląd komórek mikrogleju korowego myszy w warunkach kontrolnych (a) oraz w tkance nowotworowej (b) Zrzuty ekranu zrekonstruowanych powierzchni mikrogleju. Projekcje SOM zostały stworzone z trzech grup próbek mikrogleju z kory mózgu myszy: komórek kontrolnych (nietumorowych) (c), komórek nowotworowych (d) oraz mieszanej populacji komórek (e). Rysunek został zmodyfikowany za zgodą autora.8. <a href="https://www.jove.com/files/ftp_upload/58543/58543fig2large.jpg" target="_blank">Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.</a>
<img alt="Stany aktywności komórek z mapami ciepła wag; analiza klastrów i podzbiorów; schemat do celów badawczych." class="xfigimg" src="/files/ftp_upload/58543/58543fig3.jpg" data-altupdatedat="1788250015000" data-alt="Rysunek 3"> 
Rycina 3. (a, lewa strona) Mapa Samoorganizująca (SOM) zbioru danych mikrogleju mysiego, składającego się z 768 wektorów cech wejściowych. Zbiór danych posłużył do wytrenowania sztucznej sieci neuronowej o rozmiarze 12x12, z zastosowaniem heksagonalnej geometrii sąsiedztwa, randomizacji początkowej i 200 epok. (a, prawa strona) Odpowiadające im płaszczyzny wejściowe SOM dla pierwszych 10 komponentów DFT (b) Reakcje SOM przedstawionej na (a) na jeden losowy podzbiór plików VRML dla każdego z czterech typów komórek: „mobilnych”, „oddziałujących”, „spoczynkowych” i „fagocytujących”, opisanych pierwotnie w Rycina 5 Bayerla i in.18(c) Reakcja tej samej SOM, co w (a, lewo), na trzy losowe podzbiory całego zbioru danych (które zatem nie były częścią zbioru treningowego) powierzchni 3D typu „komórki spoczywające”. Podobieństwo między tymi trzema reakcjami jest zauważalne. Rycina została zmodyfikowana za zgodą autora8. <a href="https://www.jove.com/files/ftp_upload/58543/58543fig3large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.</a>
<img alt="Wykres analizy Fouriera przedstawiający amplitudę w funkcji punktów czasowych; dane widmowe i analiza składowych DFT." class="xfigimg" src="/files/ftp_upload/58543/58543fig4.jpg" data-altupdatedat="1788250015000" data-alt="Rysunek 4"> 
Rysunek 4. (a) Zależność od czasu pierwszych 10 składowych DFT podczas eksperymentu obrazowania intravitalnego mikrogleju mysiego. Panel ten przedstawia dane dla komórki typu „Mobile Cells”. Oś x odpowiada punktom czasowym eksperymentu z rozdzielczością czasową 60 s, oś y pokazuje amplitudę składowych DFT w jednostkach dowolnych (a.u.), natomiast oś z odpowiada składowym DFT od 1 do 10. (b) Analogicznie jak w (a), ale dla komórki typu „Interacting Cells”. Rycina została zmodyfikowana za zgodą autora.8. <a href="https://www.jove.com/files/ftp_upload/58543/58543fig4large.jpg" target="_blank">Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.</a>
<img alt="Dyskretna analiza Fouriera, amplituda a identyfikator komponentu, wykres 3D, porównanie danych widmowych." class="xfigimg" src="/files/ftp_upload/58543/58543fig5.jpg" data-altupdatedat="1788250015000" data-alt="Rycina 5"> 
Rysunek 5. (a) Zachowanie wszystkich 19 komponentów DFT komórki typu „Mobile Cells” na początku, w połowie i na końcu eksperymentu. Liczby na osi x odpowiadają identyfikatorom komponentów DFT od 1 do 19. Oś y przedstawia amplitudę komponentów DFT w jednostkach dowolnych (a.u.), natomiast oś z oznacza sześć losowych kątów obrotu. (b) Analogicznie jak w (a), ale dla komórki typu „Interacting Cells”. Rycina została zmodyfikowana za zgodą autora8. <a href="https://www.jove.com/files/ftp_upload/58543/58543fig5large.jpg" target="_blank">Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.</a>

Tutaj udostępniamy przepływ pracy, który pozwala na identyfikację zdrowych i patologicznych komórek na podstawie ich 3-wymiarowego kształtu. Opisujemy proces wykorzystania konturów projekcji 2D opartych na powierzchniach 3D do trenowania Mapy Samoorganizującej się, która zapewni obiektywne grupowanie badanych populacji komórek.

Oparte na morfologii rozróżnienie między zdrowymi i patologicznymi komórkami z wykorzystaniem transformaty Fouriera i samoorganizujących się map

The appearance and the movements of immune cells are driven by their environment. As a reaction to a pathogen invasion, the immune cells are recruited to ...

Metoda ta może pomóc odpowiedzieć na kluczowe pytania w badaniach podstawowych i zastosowaniach klinicznych, takich jak badania nad rakiem i leczenie raka, ale także w dziedzinie immunologii. Główną zaletą tej metody jest to, że jest to metoda statyczna i łatwa do ułatwienia. Ta technika jest w stanie automatycznie zidentyfikować chore typy komórek na podstawie ich kształtu i charakterystyki ruchu.Implikacje tej techniki rozciągają się na diagnostykę i terapię raka lub chorób zapalnych. Metoda ta może zapewnić wgląd w działania między komórkami odpornościowymi a komórkami rakowymi, ale może być również stosowana w każdej dziedzinie, w której niezbędne są trójwymiarowe dane mikroskopowe. Może się zdarzyć, że osoby, które są nowicjuszami w tej metodzie, będą miały trudności z opisaniem i walidacją danych z mikroskopii trójwymiarowej z powodu problemów.Pomyśleliśmy, że wizualna demonstracja tej metody jest kluczowa, ponieważ narzędzia programowe, których tutaj używamy, nie są znane wielu naukowcom. Na początek uzyskaj zestaw danych z mikroskopii trójwymiarowej o wysokiej rozdzielczości, zgodnie z opisem w dołączonym protokole tekstowym. Załaduj dane obrazu 3D do oprogramowania do rekonstrukcji i rozpocznij tworzenie powierzchni 3D dla każdego obiektu.Aby to zrobić, wybierz opcję Widok 3D i kliknij Powierzchnie, a następnie kliknij przycisk Dalej, aby kontynuować pracę z Kreatorem tworzenia powierzchni. Teraz wybierz kanał obrazu do rekonstrukcji powierzchni. Wybierz wartość wygładzenia, która nie ukrywa szczegółów powierzchni, ale także unika porowatych powierzchni.Zastosuj funkcję wygładzania, klikając pole wyboru Opcja wygładzania i podając promień wygładzania. Podczas tworzenia trójwymiarowych konstrukcji obrazowych danych mikroskopowych szczególnie ważne jest zwrócenie uwagi na współczynnik wygładzania i odpowiednią metodę progową, aby nie stracić żadnych cech komórkowych ani kształtów. Następnie wybierz metodę progowania, aby znaleźć powierzchnie.Użyj bezwzględnego progu intensywności, gdy obiekty, takie jak te pokazane tutaj, są dobrze oddzielone od tła i mają w przybliżeniu jednolity poziom jasności. Jeśli obiekty różnią się intensywnością, ale nadal można je oddzielić od lokalnego tła i od innych otaczających je obiektów, zastosuj próg kontrastu lokalnego. Ustaw lokalny obszar wyszukiwania progu zgodnie z wartością oczekiwanej średnicy zrekonstruowanych obiektów.Następnie wybierz z listy opcji, aby przefiltrować zrekonstruowane powierzchnie według interesujących ich parametrów morfologicznych. Obejmuje to objętość, sferyczność, stosunek powierzchni do objętości i inne. Zapisz i wyeksportuj wygenerowane powierzchnie w formacie takim jak VRML, który jest zgodny z oprogramowaniem do animacji 3D, które będzie używane w następnym kroku.Uruchom Blendera i przejdź do zakładki Wyjście po prawej stronie okna. Wybierz format TIFF z menu rozwijanego i ustaw głębię kolorów na 8-bitowe RGBA. Następnie przełącz się w tryb skryptów i przejdź do dostarczonego pliku skryptu zatytułowanego:GUI_Autorotate.py, który można pobrać z repozytorium GitHub dla tej pracy. Wróć do okna głównego, kliknij Uruchom skrypt i wybierz folder z plikami wrl, gdy zostaniesz poproszony o wprowadzenie. Następnie przejdź do menu domyślnego.W tym miejscu ustaw obroty na wartość równą lub wyższą niż sześć. Uruchom skrypt, klikając przycisk Obróć. Obrót o sześć różnych kątów jest zwykle wystarczający do rozróżnienia różnych populacji komórkowych, jednak nie zalecamy używania mniej niż sześciu obrotów, aby nie stracić żadnych informacji o kształtach lub cechach.Rzuty poszczególnych powierzchni należy zapisać w tym samym folderze, który został użyty do wprowadzenia danych. Domyślnie obrazy są zapisywane w 8-bitowym formacie TIFF, który jest formatem wymaganym przez wtyczkę FIJI, Shade. Otwórz Fidżi i wybierz Cień w menu Wtyczki.Zacznij od wartości domyślnych, a później dostosuj parametry. Kliknij OK, gdy wszystko będzie gotowe do uruchomienia programu. Następnie wybierz folder źródłowy, który zawiera pliki TIFF, które zostały utworzone w poprzedniej sekcji, i kliknij Wybierz.Następnie podaj folder danych wyjściowych. Gdy pojawi się pierwszy obraz, narysuj prostokąt otaczający komórkę i zacznij od kliknięcia W porządku. Podczas działania wtyczki widzisz wstępne przetwarzanie obrazu na obrzeżach, znajdowanie komórek oznaczonych czerwonymi liniami.Współrzędne znalezionych krawędzi są teraz używane do obliczania dyskretnych składowych Fouriera. Aby po raz pierwszy wytrenować samoorganizujące się mapy, zacznij od załadowania programu MATLAB. Załaduj skrypt TrainSOM MATLAP, a następnie wybierz pozycję Uruchom, aby rozpocząć trenowanie.W razie potrzeby upewnij się, że ścieżka do pliku jest poprawna. Gdy zacznie działać, pojawi się dodatkowe okno z postępem. Poczekaj, aż szkolenie zostanie zakończone, zanim przejdziesz dalej.Domyślnie skrypt jest ustawiony na uruchamianie 2 000 iteracji lub epików. Po zakończeniu trenowania sprawdź wykresy topologii sieci. Oto przykład wykresu odległości sąsiada, wykresu example-hits i wykresu input-planes.Samoorganizujące się mapy są ważnym narzędziem do odnajdywania ukrytych relacji w zestawie danych. Grupują one dane w sposób obiektywny i mają tę zaletę, że nie potrzebują zestawu danych treningowych, ponieważ uczą się bez nadzoru. Sieć jest teraz przeszkolona.Kliknij prawym przyciskiem myszy plik w swoim obszarze roboczym i zapisz go do wykorzystania w przyszłości. Załaduj samoorganizujące się mapy, gdy używasz już wytrenowanej mapy, w celu klastrowania zestawu danych. Następnie zaimportuj plik CSV, który ma być testowany ze wstępnie załadowanymi wytrenowanymi mapami.Tutaj wybierzemy wyjście CSV wtyczki Shade. Po załadowaniu danych zmień typ danych wyjściowych na Macierz liczbowa, a następnie wybierz pozycję Importuj wybór. Po zakończeniu klasyfikacji użyj okna poleceń, aby ocenić różne wykresy.Pokazany tutaj obraz pochodzi z dekonwolucji przyżyciowej wielofotonowej miskroskopii danych komórek mikrogleju. W warunkach fizjologicznych mikroglej miał dość złożony kształt z licznymi, silnie rozgałęzionymi procesami. Po umieszczeniu w środowisku nowotworowym, takim jak ten model guza korowego, mikroglej zmienił się na prostszy, bardziej wrzecionowaty kształt.20 z tych komponentów opisowych w kształcie foriera zostało wykorzystanych jako dane wejściowe do trenowania samoorganizującej się mapy. Wytrenowana mapa została następnie przetestowana w celu oceny jej zdolności do odróżniania komórek zdrowych od nowotworowych. Zdrowa populacja komórek została rzutowana na pojedynczy obszar pokazany tutaj, podczas gdy zestaw danych dotyczących mikrogleju nowotworowego został przedstawiony jako obszar aktywny w kształcie hantli.Mapy mogą być również trenowane przez ekspertów medycznych w celu identyfikacji różnych grup komórek. W tym przypadku zidentyfikowano, zrekonstruowano i wykorzystano komórki spoczynkowe, komórki fagocytozy, komórki oddziałujące i komórki ruchome do trenowania mapy 12x12. Ta połączona mapa pokazuje grupy sztucznych neuronów o wysokiej wartości trafienia, szczególnie w lewym dolnym i środkowym obszarze mapy.Solidność podejścia do mapowania została przetestowana przy użyciu wytrenowanej samoorganizującej się mapy z 3 losowymi podzbiorami tego samego typu komórek spoczynkowych, które nie były częścią zestawu danych treningowych. Odpowiedź S.O.M na to dane wejściowe wykazuje bardzo podobną odpowiedź, która wskazuje prawidłowy typ komórki. Po opracowaniu tej techniki mamy teraz zestaw narzędzi, za pomocą którego możemy w prosty sposób kategoryzować zmiany kształtu komórki.Zmiany te zachodzą, na przykład podczas raka lub innych reakcji układu odpornościowego, i możemy je zmierzyć za pomocą mikroskopii, używając całych zwierząt, wyciętych tkanek, a nawet narządów na chipie.

morphology-based-distinction-between-healthy-pathological-cells

Research

JoVE Journal

Cancer Research

Badania nad rakiem

Prezentujemy metodę mikroskopii w ciemnym polu opartą na filtrowaniu typu Gabora do pomiaru dynamiki subkomórkowej w pojedynczych żywych komórkach. Technika ta jest czuła na zmiany w strukturze organelli, takie jak fragmentacja mitochondriów.

Mikroskopia rozproszeniowa optyczna oparta na dwuwymiarowych filtrach Gabora

Opracowaliśmy platformę programową wykorzystującą Imaris Neuroscience, ImarisXT oraz MATLAB do pomiaru zmian morfologii nieokreślonego kształtu uzyskanego z trójwymiarowej konfokalnej fluorescencji pojedynczych komórek. To nowatorskie podejście może być wykorzystane do ilościowego określania zmian w kształcie komórek po aktywacji receptorów, stanowiąc tym samym potencjalne dodatkowe narzędzie w procesie odkrywania leków.

Narzędzie analityczne do ilościowej oceny zmian morfologii komórkowej na podstawie trójwymiarowych obrazów fluorescencyjnych

Źródło: Huang, Y. et. al., Podłużne monitorowanie morfologiczne i fizjologiczne trójwymiarowych sferoid nowotworowych za pomocą optycznej tomografii koherentnej. J. Vis. Exp. (2019)Ten film demonstruje zastosowanie koherencji optycznej tomografia (OCT) do monitorowania sferoid guza w płytce wielodołkowej poprzez przechwytywanie obrazów przekroju poprzecznego o wysokiej rozdzielczości. Analizując odbicia światła od struktur komórkowych, OCT konstruuje modele 3D w celu oceny cech morfologicznych,...

Monitorowanie morfologiczne i fizjologiczne sferoid nowotworowych za pomocą optycznej koherentnej tomografii 

Demonstrujemy metody wykrywania zniekształceń architektonicznych w poprzednich mammografach. Zorientowane struktury są analizowane za pomocą filtrów Gabora i portretów fazowych w celu wykrycia miejsc promieniowania wzorców tkankowych. Każde miejsce jest scharakteryzowane i sklasyfikowane za pomocą miar reprezentujących wzory spikulacyjne. Metody te powinny pomóc w wykrywaniu raka piersi.

Wykrywanie zniekształceń architektonicznych we wcześniejszych mammografiach poprzez analizę zorientowanych wzorców

Ten artykuł opisuje szczegółową metodologię przygotowania matrycy Multiplexed Artificial Cellular MicroEnvironment (MACME) do wysokoprzepustowej manipulacji fizycznych i chemicznych sygnałów naśladujących mikrośrodowiska komórkowe in vivo oraz do identyfikacji optymalnego środowiska komórkowego dla ludzkich pluripotencjalnych komórek macierzystych (hPSC) z profilowaniem pojedynczych komórek.

Wytwarzanie multipleksowanej sztucznej macierzy mikrośrodowiska komórkowego

Mikroglej to komórki odpornościowe mózgu, które badają i reagują na zmienione fizjologie mózgu poprzez zmiany morfologiczne, które mogą być oceniane ilościowo. Protokół ten przedstawia protokół analizy oparty na ImageJ, aby przedstawić morfologię mikrogleju jako ciągłe dane zgodnie z takimi wskaźnikami, jak rozgałęzienie komórki, złożoność i kształt.

Kwantyfikacja morfologii mikrogleju na podstawie mikrofotografii immunohistochemicznej tkanki przygotowanej za pomocą ImageJ

Przedstawiona tutaj metoda wykorzystuje mikrowzorce wraz z obrazowaniem ilościowym, aby ujawnić organizację przestrzenną w kulturach ssaków. Technika ta jest łatwa do wdrożenia w standardowym laboratorium biologii komórki i oferuje elastyczny system do badania wzorców in vitro.

Mapowanie wyłaniającej się organizacji przestrzennej komórek ssaków za pomocą mikrowzorców i obrazowania ilościowego

Ten artykuł wyjaśnia, jak używać nadzorowanego symulacji uczenia maszynowego do analizy morfologii mitochondriów w obrazach z mikroskopii fluorescencyjnej stałych komórek.

Analiza morfologii mitochondriów za pomocą symulacji uczenia nadzorowanego

Protokół ten przedstawia multimodalne podejście łączące spektroskopię Ramana i tomograficzną mikroskopię fazową w celu bezznacznikowego, nieinwazyjnego fenotypowania żywych ludzkich komórek raka piersi. Zawiera przyjazny dla użytkownika, wolny od uprzedzeń człowieka potok analityczny do szybkiej, ilościowej analizy morfochemicznej, dostarczając szczegółowych informacji molekularnych i strukturalnych z zastosowaniami w badaniach nad rakiem i biologii komórki.

Ramy obrazowania multimodalnego w celu przyspieszenia fenotypowania żywych, bezznacznikowych komórek raka piersi

Oprogramowanie do mapowania parametrów mikrogleju (MMPS) to przyjazny dla użytkownika system analityczny, który nie wymaga znajomości kodowania i półautomatycznie analizuje morfologię pojedynczych komórek mikrogleju na obrazach immunofluorescencyjnych. MMPS, zweryfikowane w modelu neurozapalenia wywołanego lipopolisacharydem u gryzoni, w sposób powtarzalny wykrywa zmiany morfologiczne związane z aktywacją i umożliwia standaryzowane, skalowalne fenotypowanie mikrogleju.

Oprogramowanie do mapowania parametrów mikrogleju (MMPS): Otwarte, przyjazne dla użytkownika narzędzie do ilościowej analizy morfologii mikrogleju

Poznaj metodę automatycznej identyfikacji komórek zdrowych i patologicznych przy użyciu transformaty Fouriera i samoorganizujących się map. Technika ta wspiera diagnostykę nowotworów i chorób zapalnych poprzez analizę morfologii trójwymiarowej.

Oparte na morfologii rozróżnienie między zdrowymi i patologicznymi komórkami z wykorzystaniem transformaty Fouriera i samoorganizujących się map (Wideo) | JoVE Journal

Oparte na morfologii rozróżnienie między zdrowymi i patologicznymi komórkami z wykorzystaniem transformaty Fouriera i samoorganizujących się map

Przeglądaj więcej filmów