Wygenerowaliśmy zestaw obrazów składający się z 25 pól/dołka, 54 dołków/płytki (3 populacje komórek x 6 stężeń leku x 3 powtórzenia), na trzech płytkach, co dało łącznie 4 050 pojedynczych obrazów. Zestawy obrazów wygenerowane w trakcie eksperymentu przeanalizowano przy użyciu autorskiego oprogramowania (patrz tabela materiałów) w celu wyodrębnienia różnych właściwości ilościowych komórek (t.j. morfologii, fluorescencji), które można było następnie wykorzystać do klasyfikacji subpopulacji komórek. Jednakże, ze względu na ograniczony dostęp do użytego oprogramowania komercyjnego, stworzono porównywalne potoki przetwarzania danych w programach CellProfiler i CellProfiler Analyst.
Klasyfikacja heterokomórkowa na subpopulacje
Jądra komórkowe zidentyfikowano i wysegmentowano na podstawie barwienia DNA (w tym przypadku Hoechst), a populacje komórek sklasyfikowano na podstawie fluorescencji lub morfologii (Rycina 1). W przypadku klasyfikacji opartej na fluorescencji, fibroblasty (CCD-19Lu) zostały wcześniej transdukowane lentiwirusem GFP. Zmierzono poziomy intensywności GFP dla każdego jądra; te, których wartości obliczone przekroczały przyjęty próg (oparty na sygnale tła), sklasyfikowano jako CCD-19Lu, natomiast te poniżej progu zidentyfikowano jako komórki nowotworowe (H3255). W przypadku klasyfikacji opartej na morfologii, komórki wcześniej zbarwiono nietoksycznym barwnikiem komórkowym (patrz tabela materiałów), co posłużyło do identyfikacji i segmentacji cytoplazmy. Algorytm uczenia maszynowego przeszkolono na ~50-100 komórkach z każdej populacji. Zidentyfikowano cechy morfologiczne, które istotnie różniły się między populacjami, a następnie wykorzystano je do zaprojektowania liniowego klasyfikatora pozwalającego rozróżnić komórki CCD-19Lu od H3255. Protokoły klasyfikacji fluorescencyjnej i morfologicznej wykazały zgodność na poziomie 97,4% (n = 1403) przy rozróżnianiu dwóch populacji komórek w warunkach nieleczonych oraz 92,5% (n = 916) w warunkach leczenia lekiem (1 µM erlotynibu) (Rycina 2).
Analizy fenotypowe subpopulacji
Oprócz rozróżniania typów komórek, dążyliśmy do scharakteryzowania właściwości fenotypowych każdej subpopulacji. Analizy multipleksowe oszczędzają czas i odczynniki, zwiększają spójność oraz dostarczają dodatkowych informacji na temat badanego układu. Istnieje wiele potencjalnych parametrów fenotypowych, które należy dobrać w zależności od stawianych pytań badawczych. W niniejszej pracy zbadano zmiany w morfologii komórek i ich żywotności w odpowiedzi na leczenie erlotynibem. Po trzech dniach podawania leku zaobserwowano zmniejszenie powierzchni jądra oraz zwiększenie powierzchni komórek H3255 (Rysunek 3A). Średnia różnica w powierzchni jądra między populacjami „bez leku” a „leczonymi” okazała się statystycznie istotna w dwustronnym teście t typu 2 (równa wariancja) (p = 7.92 x 10-16). Hipotezą jest, że obserwacja ta stanowi odpowiedź komórkową na stres wywołany leczeniem farmakologicznym.
Interesujące jest również zbadanie, czy lek wykazuje działanie cytotoksyczne (tj. wzrost liczby martwych komórek w czasie) lub cytostatyczne (tj. spadek liczby nowo powstających komórek w czasie), ponieważ ma to istotny wpływ kliniczny. Na przykład cytostatyczne działanie leku indukuje zatrzymanie wzrostu, ale nie eliminuje komórek z guza, co stwarza potencjał do ponownego zainicjowania proliferacji komórek nowotworowych po odstawieniu leku. Efekty działania leków są często zależne od kontekstu, stężenia oraz typu komórek. Wcześniej zaobserwowaliśmy, że erlotynib wywołuje odpowiedź cytotoksyczną w jednym typie komórek, podczas gdy w innym wykazuje odpowiedź cytostatyczną13.
Tradycyjne testy żywotności podają względną liczbę komórek i w związku z tym nie odróżniają zatrzymania wzrostu od śmierci komórkowej. W niniejszej pracy martwe komórki zidentyfikowano na podstawie barwienia jodkiem propidyny (Rycina 3B). Zaobserwowano zarówno efekty cytotoksyczne, jak i cytostatyczne erlotynibu w komórkach H3255, objawiające się wzrostem liczby zgonów i spadkiem liczby nowych komórek po zastosowaniu leku (Rycina 3C). Warto zauważyć, że liczba martwych komórek spada po pierwszym dniu, co prawdopodobnie wynika z usuwania szczątków komórkowych. Komórki CCD-19Lu nie zostały poddane wpływowi leku. Dodatkową zaletą tej platformy jest generowanie danych ilościowych. Na przykład w naszym eksperymencie z kokulturą początkowa subpopulacja 1,118 (75,8%) komórek H3255 zmieniła się do 2,817 (87,9%) lub 396 (57,2%) po trzech dniach odpowiednio bez lub z zastosowaniem erlotynibu (Rycina 4). Ponieważ możemy generować rzeczywiste zliczenia komórek zamiast względnych wartości procentowych (jak w metodach cytometrii przepływowej), wyciągamy wniosek, że zmiana składu podczas leczenia jest wynikiem spadku liczby komórek H3255, a nie wzrostu liczby komórek CCD-19Lu. Należy odnotować, że wskaźniki śmierci mogą być niedoszacowane ze względu na usuwanie szczątków komórkowych, co jest trudne do oceny eksperymentalnej i prawdopodobnie różni się w zależności od typu komórek.

Rycina 1: Przegląd protokołu analizy obrazu. Dwie potencjalne ścieżki analizy obrazu do klasyfikacji populacji heterokomórkowych z wykorzystaniem klasyfikacji opartej na morfologii lub klasyfikacji opartej na fluorescencji. Paski skali = 100 µm. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 2: Zgodność między klasyfikacją opartą na morfologii a klasyfikacją opartą na fluorescencji. (A) Wykres zgodności przedstawiający nakładanie się dwóch protokołów klasyfikacji. Te same komórki zostały sklasyfikowane jako H3255 przy użyciu klasyfikacji opartej na morfologii i fluorescencji. Oba protokoły były zgodne w przypadku 97,4% (n = 1403) komórek nieleczonych oraz 92,5% (n = 916) komórek traktowanych erlotynibem (Uwaga: biały obszar jest zbyt mały, aby był widoczny). (B) Obrazy 10X przedstawiające przykłady dobrej i słabej zgodności między klasyfikacją opartą na fluorescencji a klasyfikacją opartą na morfologii. Białe strzałki wskazują komórki, które zostały niespójnie sklasyfikowane między platformami. Obraz wejściowy: niebieski – jądra (Hoechst); zielony – CCD19Lu (GFP). Obrazy klasyfikacji: czerwony – H3255; zielony – CCD-19Lu. Pasek skali = 100 µm. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 3: Multipleksowe pomiary fenotypowe z jednego układu eksperymentalnego. (A) Cechy morfologiczne, takie jak jądra i powierzchnia komórek, zostały obliczone na poziomie pojedynczych komórek w obecności i przy braku leku. Uwaga: Powierzchnie komórek mniejsze niż 100 µm2 uznano za debris i wykluczono z analiz. Wykres pudełkowy przedstawia medianę z zakresami pierwszego i trzeciego kwartyla oraz słupkami błędów dla 95% przedziału ufności. (B) Komórki H3255 (niebieskie) i CCD-19Lu (zielone) były hodowane wspólnie, a martwe komórki zidentyfikowano na podstawie intensywności barwienia jodkiem propidium (czerwony) i obrazowano za pomocą obiektywu 10X. Pasek skali = 1 mm (panel górny); 100 µm (obrazy dolne). (C) Całkowitą liczbę żywych i martwych komórek obliczono w ciągu trzech dni z zastosowaniem lub bez leczenia lekiem, wykazując wyraźny spadek liczby żywych komórek i wzrost liczby martwych komórek po dodaniu erlotynibu. Słupki błędów reprezentują błąd standardowy średniej na podstawie trzech powtórzeń. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 4: Dynamika subpopulacji w czasie. Reprezentatywne obrazy 10X studni zawierających komórki H3255 (niebieskie) i CCD-19Lu (zielone) w dniu 0 lub dniu 3, z lekiem i bez leku. Policzono komórki należące do każdej subpopulacji, a proporcjonalne wykresy kołowe przedstawiają rzeczywistą zmianę składu populacji w poszczególnych próbkach. Paski skali = 1 mm (panele środkowe, pokazane na panelu skrajnie lewym), 1 mm (obraz górny), 100 µm (obrazy dolne). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.