Opracowano otwartoźródłową, samodzielną platformę do analizy morfologii mikrogleju, nazwaną Mapping Microglial Parameters Software (MMPS), aby zapewnić dostępny i ustandaryzowany proces analizy morfologii mikrogleju w oparciu o immunofluorescencję (IF). Platforma opiera się na wcześniejszych metodologiach analizy morfologicznej, nie wymaga wiedzy z zakresu programowania i posiada przystępny interfejs graficzny użytkownika. Program MMPS został zaprojektowany w celu poprawy wykrywania subtelnych rozgałęzień mikrogleju, przy jednoczesnym zastosowaniu wielu zabezpieczeń minimalizujących ryzyko przeszacowania masek mikrogleju.
MMPS zostało zwalidowane z wykorzystaniem modelu neurostanu zapalnego u gryzoni indukowanego LPS. Po obróbce tkanek i akwizycji obrazów (Rycina 1), swoistość IF potwierdzono przy użyciu przekrojów kontrolnych negatywnych, pozbawionych barwienia przeciwciałem pierwotnym (Rycina uzupełniająca 1). Przetwarzanie w MMPS rozpoczyna się od importu obrazów IF do interfejsu oprogramowania (Rycina 2A). Oprogramowanie oferuje wiele opcji przetwarzania obrazu w celu poprawy detekcji sygnału mikrogleju i separacji tła. Obowiązkowe odejmowanie tła jest wykonywane za pomocą algorytmu rolling-ball zaimplementowanego w pakiecie Scikit-image, natomiast opcjonalne filtry odszumiające i wyostrzające są dostępne poprzez filtrowanie gaussowskie oparte na SciPy. Po przetwarzaniu obrazu somy mikrogleju są zaznaczane ręcznie i identyfikowane za pomocą czerwonych znaczników wyboru (Rycina 2B). Granice som są następnie obrysowywane ręcznie w celu przeprowadzenia dalszych obliczeń morfologicznych (Rycina 2C).

Rysunek 1: Graficzny przegląd modelu neurostanu zapalnego i procedury przetwarzania. Szczurom wstrzyknięto LPS lub PBS, a po 24 h pobrano próbki poprzez perfuzję sercową z następującą po niej dekapitacją. Wyekstrahowano mózgi i wykonano przekroje o grubości 30 µm, które następnie wykorzystano w IF dla markera mikroglei Iba1, a następnie przetworzono za pomocą MMPS. Stworzone w Biorender. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 2: Przegląd zastosowania MMPS. (A) Przykład obrazu mikrogleju korowego (czerwony, Iba1) z jądrami komórkowymi (niebieski, DAPI) do analizy w MMPS (pasek skali = 75 µm, 40x). Po przetworzeniu obrazu za pomocą MMPS tło zostanie zredukowane, a morfologie zostaną dokładnie ocenione (B) Po przetworzeniu obrazu w MMPS somy są ręcznie zaznaczane czerwonym okręgiem do wykorzystania w kolejnych krokach szczegółowo opisanych w metodach (pasek skali = 75 µm, 40x). (C) Dwie reprezentatywne komórki mikrogleju zostały ręcznie obrysowane; poniżej przedstawiono zapisane obrazy obrysów somy. (D) Generowanie maski wybranych komórek automatycznie śledzi rozgałęzienia mikrogleju do 800 µm2, jednak wymaga to ręcznej kontroli jakości. Tutaj przedstawiono niedokładną maskę 800 µm2, a następnie dokładniejszą maskę 700 µm2. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.
Po obrysowaniu somy, program MMPS iteracyjnie generuje maski w zakresie wielkości 200–800 µm2 , stosując adaptacyjne progowanie ograniczone przez zdefiniowane przez użytkownika progi minimalnej intensywności. Mniejsze maski ściślej przestrzegają zdefiniowanego progu, natomiast w przypadku większych masek ograniczenia progowe są stopniowo rozluźniane, aby zmaksymalizować uwzględnienie wypustek. Wygenerowane maski poddawane są ręcznej kontroli jakości, aby zapobiec włączeniu sąsiednich komórek, nadmiernemu uwzględnieniu tła lub ucięciu przez krawędzie obrazu lub paski skali. Jak pokazano na Rysunku 2D, maska o wielkości 800 µm2 dla reprezentatywnej komórki mikrogleju rozciągała się na sąsiednią komórkę. Odrzucenie maski 800 µm2 i wybór maski 700 µm2 pozwoliło skorygować błąd segmentacji, zachowując jednocześnie dokładne wyznaczenie granic docelowej mikroglei.
Po przetworzeniu metodą MMPS oceniono aktywację mikrogleju, wykorzystując największą zaakceptowaną maskę dla każdej komórki, ponieważ maski te najdokładniej oddawały dystalne rozgałęzienia, minimalizując jednocześnie włączanie tła. Obliczone parametry morfologiczne podsumowano w Tabeli 1 i przedstawiono na Rysunku 3. Podanie LPS istotnie zmniejszyło średnią odległość centroidu oraz obwód komórki w porównaniu z kontrolą traktowaną nośnikiem (p < 0,05; Rysunek 3A,B), co jest zgodne z retrakcją wypustek mikrogleju podczas aktywacji. Nie zaobserwowano istotnych różnic w powierzchni maski, zaokrągleniu ani ekscentryczności pomiędzy grupami badawczymi (Rysunek 3C–E), pomimo wcześniejszych doniesień wiążących te parametry z fenotypami aktywowanego mikrogleju3,23. Niniejsza analiza nie posiadała mocy statystycznej pozwalającej na wykrycie subtelnych różnic w kształcie somy lub powierzchni maski, a subpopulacje mikrogleju mogą przyjmować heterogeniczne morfologie zapalne, które nie wykazują tych trendów w sposób jednolity3. Z kolei powierzchnia somy była istotnie zwiększona u zwierząt traktowanych LPS w porównaniu z grupą kontrolną z nośnikiem (p < 0,01; Rysunek 3F), co dodatkowo potwierdza przebudowę morfologiczną związaną z aktywacją22.

Rysunek 3: Porównanie parametrów morfologicznych mikrogleju. Dla wszystkich paneli przeprowadzono testy normalności (Shapiro-Wilka i Levene'a), w których we wszystkich przypadkach uzyskano wartości nieistotne (p>0,05), co doprowadziło do zastosowania testów t-Studenta. Kropki reprezentują poszczególne zwierzęta (n=5 na grupę leczoną) (A) Średnia odległość centroidu była istotnie większa w grupie kontrolnej (vehicle) w porównaniu z grupą leczoną LPS. (B) Obwód komórek mikrogleju był istotnie większy w grupie kontrolnej w porównaniu z grupą leczoną LPS. (C) Nie stwierdzono istotnej różnicy w powierzchni maski pomiędzy grupami leczonymi. (D) Nie stwierdzono istotnej różnicy w ekscentryczności somy pomiędzy grupami. (E) Nie ma istotnej różnicy w okrągłości somy pomiędzy grupami leczonymi. (F) Odnotowano istotny wzrost powierzchni somy w grupie LPS * p<0,05, ** p<0,01, *** p<0,001. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.
| Parametr | Definicja | Związek z aktywacją mikrogleju |
| Średnia odległość centroidów | Średnia odległość najdalszej gałęzi w pozycjach: górnej, prawej, dolnej i lewej | Zmniejszenia |
| Obwód | Pomiar obwodu komórki | Zmniejsza |
| Pole maskowania | Rzeczywisty obszar komórki, a nie obszar maski celu | Zmniejszenia |
| Ekscentryczność | Pomiar przybliżonego wydłużenia somy komórki (0: okrągła, 1: silnie wydłużona) | Zmniejszenia |
| Okrągłość | Pomiar przybliżonej kolistości somy komórki (0: wydłużona, 1: kolista) | Wzrosty |
| Wielkość somy | Pole powierzchni somy komórki | Zwiększenia |
Tabela 1: Parametry obliczone za pomocą programu MMPS i przewidywana zmiana w odpowiedzi na aktywację mikrogleju.Tabela podsumowuje parametry morfologiczne obliczone przez program Mapping Microglial Parameters Software (MMPS), w tym obwód, średnią odległość od centroidu, powierzchnię maski, ekscentryczność, okrągłość oraz powierzchnię somy.
Następnie oceniono powtarzalność międzyużytkownikową, przekazując ten sam zestaw obrazów do niezależnej analizy osobom zaślepionym i wcześniej nieszkolonym (Supplementary Figure 2). Nie stwierdzono statystycznie istotnych różnic między użytkownikami, a wszyscy oni odtworzyli trendy morfologiczne zgodne z aktywacją mikrogleju u zwierząt traktowanych LPS.
MMPS dodatkowo umożliwia zaawansowane analizy morfologiczne w dalszych etapach. Analiza korelacyjna parametrów wygenerowanych na Rysunku 3 wykazała związki między odrębnymi cechami morfologicznymi. U szczurów traktowanych placebo zaobserwowano jedynie słabą korelację między okrągłością a obwodem mikrogleju (R2 = 0,17), natomiast u szczurów traktowanych LPS stwierdzono silniejszą korelację (R2 = 0,64), co sugeruje, że zwiększona aktywacja mikrogleju wiązała się z jednoczesnym zmniejszeniem okrągłości i obwodu (Rysunek 4A). Ponadto pole maski wykazało silne dodatnie korelacje z przeciętną odległością centroidu w obu grupach – traktowanej placebo oraz LPS (odpowiednio R2 = 0,71 i 0,88), co jest zgodne ze zmianami w złożoności rozgałęzień związanych ze stanem aktywacji mikrogleju (Rysunek 4B).

Rysunek 4: Zaawansowana analiza parametrów mikrogleju. Wszystkie panele przedstawiają linie dopasowania na poziomie zwierząt, obliczone metodą zwykłych najmniejszych kwadratów, wraz z współczynnikiem determinacji (R2) oraz zacieniowanym 95% przedziałem ufności. (A) Okrągłość somy mikrogleju i całkowita długość obwodu wykazują słabą tendencję w mikrogleju z grupy nośnika (R2 = 0,17), natomiast korelacja jest silniejsza po ekspozycji na LPS (R2 = 0,64). (B) Średnia odległość od centroidu koreluje z powierzchnią mikrogleju w grupach z nośnikiem (R2 = 0,71) i LPS (R2 = 0,88). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Rycina uzupełniająca 1: Walidacja protokołu IF za pomocą kontroli negatywnej. Jeden skrawek mózgu z tego samego obszaru dla każdego szczura poddano protokołowi IF bez dodania przeciwciała pierwszorzędowego, aby potwierdzić brak autofluorescencji (pasek skali = 75 µm) (A) Reprezentatywny obraz tkanki barwionej nośnikiem dla DAPI (niebieski). (B) Reprezentatywny obraz tkanki barwionej nośnikiem dla Iba1 (czerwony). (C) Reprezentatywny obraz tkanki barwionej nośnikiem dla DAPI (niebieski) i Iba1 (czerwony). (D) Reprezentatywny obraz tkanki barwionej LPS dla DAPI (niebieski). (E) Reprezentatywny obraz tkanki barwionej LPS dla Iba1 (czerwony). (F) Reprezentatywny obraz tkanki barwionej LPS dla DAPI (niebieski) i Iba1 (czerwony). Proszę kliknąć tutaj, aby pobrać ten plik.
Rycina dodatkowa 2: Ocena zmienności międzyużytkownikowej MMPS. Kilku oślepionych użytkowników zastosowało MMPS w tych samych próbkach z Ryciny 3, a następnie porównano wyniki na poziomie zwierząt. Dla wszystkich próbek testy normalności (Shapiro-Wilka i Levene'a) okazały się nieistotne (p>0.05), co doprowadziło do zastosowania ANOVA dla pomiarów powtarzanych z post-hoc skorygowanym testem t dla par według metody Holma (p>0.05). (A) W każdej grupie nie stwierdzono istotnej różnicy między użytkownikami pod względem średniej odległości centroidu (p>0.05). (B) W każdej grupie nie stwierdzono istotnej różnicy między użytkownikami pod względem obwodu (p>0.05). (C) W każdej grupie nie stwierdzono istotnej różnicy między użytkownikami pod względem powierzchni maski (p>0.05). (D) W każdej grupie nie stwierdzono istotnej różnicy między użytkownikami pod względem ekscentryczności (p>0.05). (E) W każdej grupie nie stwierdzono istotnej różnicy między użytkownikami pod względem okrągłości (p>0.05). (F) W każdej grupie nie stwierdzono istotnej różnicy między użytkownikami pod względem powierzchni somy (p>0.05). Kliknij tutaj, aby pobrać ten plik.