Wygenerowano obrazy 10-pleksowej IF oraz odpowiadające im wirtualne obrazy H&E dla 136 rdzeni macierzy tkankowej, obejmujące 34 typy tkanek ludzkich na dwóch slajdach przygotowanych różnymi metodologiami (Rycina 2). Przedstawione reprezentatywne obrazy demonstrują udane barwienie multipleksowe z zachowaną morfologią tkanki i wyraźną detekcją wielu biomarkerów w rdzeniach macierzy tkankowej. Jeden slajd został przygotowany przy użyciu automatycznego barwiarki do slajdów, a drugi w procesie ręcznego przygotowania slajdów. Większą częstość występowania artefaktów tkankowych, w tym odwarstwiania tkanki, a także utraty lub przesunięcia rdzeni, zaobserwowano przy zastosowaniu ręcznej metody opartej na komorze do odmaskowywania (Rycina 3). Artefakty te zmniejszają analizowalny obszar tkanki na próbkę lub rdzeń, co wpływa na późniejszą analizę obrazów multipleksowych. Obserwacje te potwierdzają, że protokół umożliwia niezawodne barwienie multipleksowe i obrazowanie w różnych typach tkanek.
Wskaźnikami pomyślnego obrazowania multipleksowego są zachowana morfologia tkanki, spójne dopasowanie pomiędzy poszczególnymi rundami barwienia oraz wyraźna i swoista detekcja sygnału biomarkerów. Z kolei nieoptymalne wyniki charakteryzują się odwarstwieniem tkanki, utratą sygnału lub artefaktami obrazowania, które ograniczają obszary poddane analizie i utrudniają dalszą obróbkę danych.
Następnie porównano profile barwienia markerów, w tym PANCK (tkanka nabłonkowa), CD31 (tkanka śródbłonkowa), aSMA (tkanka stropmalna), keratynę 15 (K15; nabłonek podstawny), vimentynę (komórki mezenchymalne), Ki67 (marker proliferacji), CD8 (cytotoksyczne limfocyty T) oraz FOXP3 (TREGs), dla obu technik przygotowania preparatów, przy czym w metodzie ręcznej zaobserwowano jakościowo (wizualnie) niewielki spadek sygnału (Rysunek 4A). Wykorzystując progi klasyfikacji biomarkerów w oprogramowaniu do patologii cyfrowej, próbki podzielono na populacje pozytywne i negatywne dla każdej z ok. 630 tys. komórek w zbiorze danych (Rysunek 4B). Ilościową ocenę sygnału markerów przeprowadzono na podstawie pomiarów intensywności uzyskanych z cyfrowego przepływu analizy patomorfologicznej. Szacunki gęstości populacji (gdzie pole pod krzywą reprezentuje 100% populacji) oraz wartości mediany intensywności wykazały niewielkie wzrosty dla 5 z 8 biomarkerów w populacjach pozytywnych przy zastosowaniu automatycznego barwienia preparatów w porównaniu do procesu ręcznego (Rysunek 4C).
Aby dopełnić panel, opisane powyżej markery zestawiono z dodatkowymi biomarkerami dla Lamin B1 (błona jądrowa) oraz p16 (inhibitor kinazy zależnej od cyklin), co pozwoliło na wizualizację względnej liczebności i rozkładu wszystkich 10 markerów w obrębie pojedynczej próbki (Rycina 5). Zdolność do profilowania dużej liczby biomarkerów jest szczególnie użyteczna w badaniach nad starzeniem się komórek (senescencją), które opierają się na ocenie wielu ortogonalnych markerów w celu zapewnienia prawidłowej klasyfikacji2. W tym przypadku udało nam się sprofilować zmiany intensywności Lamin B1 w populacjach z wysoką i niską ekspresją p16. Co ważne, ponieważ była to próbka multipleksowa, komórki te mogły zostać osadzone w kontekście przestrzennym poprzez porównanie ich z rozkładem innych typów komórek w tej samej tkance (Rycina 6A). Za pomocą progów klasyfikacji w oprogramowaniu Digital Pathology Software, barwienie to poddano kwantyfikacji, aby wykazać, że komórki p16-dodatnie wiązały się z niższymi poziomami Lamin B1 niż ich odpowiedniki p16-ujemne (Rycina 6B), co jest profilem zgodnym z wcześniej opisanymi wzorcami markerów związanymi z senescencją12. Zaobserwowaliśmy również, że komórki Lamin B1-ujemne wykazywały wyższe poziomy p16 (Rycina 6C). Wykorzystując funkcję fenotypów w oprogramowaniu Digital Pathology Software (która klasyfikuje komórki na podstawie segmentacji pojedynczych komórek i progowania wielu markerów), kombinację tych markerów wykorzystano do zdefiniowania populacji komórek jako senescentnych, tj. takich, które były jednocześnie p16-dodatnie i Lamin B1-ujemne. Następnie dokonano kwantyfikacji proporcji tych komórek we wszystkich rdzeniach w zależności od typu tkanki, a także każdego z 10 markerów z osobna (Rycina 6D). Wynik ten umożliwia ilościowe porównanie populacji komórek zdefiniowanych przez biomarkery w wielu typach tkanek w obrębie zbioru danych multipleksowych. Choć potok analizy działa na poziomie pojedynczych komórek, obecny zbiór danych jest przedstawiony głównie za pomocą rozkładów na poziomie populacyjnym. Stosunkowo niska liczebność komórek p16-dodatnich w tym zbiorze danych może ograniczać moc statystyczną dla rzetelnej analizy korelacji komórka po komórce między p16 a Lamin B1, dlatego takie zależności należy interpretować z ostrożnością.

Rycina 1: Schemat przedstawiający pełny przebieg pracy z immunofluorescencją multipleksową. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 2: Przegląd obrazowania multipleksowego: (A-B): Wirtualne obrazy H&E wygenerowane z barwienia DAPI i autofluorescencji Cy3 z rundy 1 w celu wizualizacji integralności tkanki na preparatach przygotowanych za pomocą automatycznego barwiarki do slajdów lub procesów manualnego przygotowania preparatów. (C-D): Multipleksowe barwienie IF wszystkich biomarkerów. Prosimy kliknąć tutaj, aby zobaczyć powiększoną wersję tej ryciny.

Rysunek 3: Przykłady typów artefaktów: Reprezentatywne przykłady różnych typów artefaktów, które mogą powstać głównie podczas przygotowywania preparatów. Częściowe odklejenie tkanki objawia się jako obszary poza ogniskową i może powodować problemy z dopasowaniem tkanki w kolejnych rundach. Brakujące rdzenie zmniejszają moc eksperymentalną i często występują nieproporcjonalnie częściej w delikatnych typach tkanek. Przesunięte lub nakładające się rdzenie powstają, gdy rdzeń odkleja się od swojej pozycji na TMA, ale przykleja się do innego obszaru, co potęguje stratę próbki, utrudniając analizę ukrytego rdzenia. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rysunek 4: Porównanie sygnału fluorescencji uzyskanego za pomocą automatycznego barwiarki BOND oraz ręcznego przygotowania preparatów: (A): Reprezentatywne barwienie IF dla markerów strukturalnych: pan-cytokeratyny (PANCK), CD31 oraz alfa-aktyny mięśni gładkich (aSMA). W pierwszej kolumnie przedstawiono nałożenie tych trzech obrazów. (B): Maski analizy klasyfikacji binarnej z oprogramowania do patologii cyfrowej, wskazujące na wykrycie pozytywnych markerów w różnych typach komórek: CD31 (czerwony), FOXP3 (turkusowy), Ki67 (żółty), CD8 (biały), PANCK (zielony). Biała ramka wskazuje obszar cyfrowego powiększenia w dolnych panelach. (C): Ilościowe oznaczenie intensywności biomarkerów w populacjach komórek pogrupowanych według metody przygotowania tkanek (automatyczna barwiarka do preparatów = czerwony, metoda ręczna = niebieski) oraz pozytywności markera (ciemny = negatywny, jasny = pozytywny). Pole pod każdym wykresem gęstości reprezentuje 100% analizowanej populacji komórek. AF = autofluorescencja. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 5: Pełny panel multipleksowy: Reprezentatywne obrazy z tkanki wyrostka robaczkowego pokazujące A-B: pełny 10-pleksowy panel IF C: PANCK (zielony), D: CD31 (czerwony), E:aSMA (żółty), F: K15 (różowy), G: Wimentyna (biała), H: Lamin B1 (turkusowy), I:p16 (złoty), J: Ki67 (jasnożółty), K: CD8 (jasnozielony) oraz L: FOXP3 (fioletowy). Biała ramka = powiększenie cyfrowe (Merge, CD31, aSMA, K15, Vim, LB1, p16), ramka przerywana = powiększenie cyfrowe (Ki67, CD8, FOXP3). Aby zobaczyć powiększoną wersję tego rysunku, kliknij tutaj.

Rysunek 6: Wykrywanie markerów starzenia: (A): Reprezentatywne obrazy IF z tkanki trzustki ukazujące kontrastujące wzorce barwienia Laminy B1 (turkusowy) i p16 (złoty) wraz z relacją przestrzenną do innych markerów: PANCK (zielony), CD31 (czerwony) i wimentyny (biały). (B): Ilościowe oznaczenie intensywności jąder Laminy B1 (LB1) w populacjach komórek p16-dodatnich i p16-ujemnych. (C): Ilościowe oznaczenie intensywności komórkowej p16 w populacjach komórek Laminy B1 (LB1)-dodatnich i Laminy B1 (LB1)-ujemnych. Pole pod krzywą na wykresie gęstości reprezentuje 100% populacji. (D): Z-score znormalizowana mapa ciepła średniego procentowego odsetka pozytywności dla każdego biomarkera lub fenotypu starzenia (p16+/LB1-ve) w różnych typach tkanek. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.