Wygenerowano zestaw obrazów składający się z 25 pól/studzienkę, 54 studzienek/płytkę (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 za pomocą autorskiego oprogramowania (patrz tabela materiałów), aby wyodrębnić różne właściwości ilościowe komórek (t.j. morfologię, fluorescencję), które następnie posłużyły do klasyfikacji subpopulacji komórek. Jednakże ze względu na ograniczony dostęp do stosowanego oprogramowania komercyjnego, stworzono porównywalne potoki analizy w programach CellProfiler i CellProfiler Analyst.
Klasyfikacja heterokomórkowa na subpopulacje
Jądra komórkowe zidentyfikowano i poddano segmentacji na podstawie barwienia DNA (w tym przypadku Hoechst), a populacje komórek sklasyfikowano albo na podstawie fluorescencji, albo morfologii (Rycina 1). W przypadku klasyfikacji opartej na fluorescencji, fibroblasty (CCD-19Lu) zostały wcześniej transdukowane lentiwirusem GFP. Dla każdego jądra zmierzono poziomy intensywności GFP; te, których wartości obliczono powyżej przyjętego progu (opartego na sygnale tła), zaklasyfikowano jako CCD-19Lu, natomiast wartości poniżej progu zidentyfikowano jako komórki nowotworowe (H325). W klasyfikacji opartej na morfologii komórki zostały wcześniej zabarwione nietoksycznym barwnikiem komórkowym (patrz tabela materiałów), co posłużyło do identyfikacji i segmentacji cytoplazmy. Algorytm uczenia maszynowego wytrenowano na ~50-10 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 odróżnić komórki CCD-19Lu od H325. Protokoły klasyfikacji fluorescencyjnej i morfologicznej były zgodne w 97,4% (n = 1403) przy rozróżnianiu dwóch populacji komórek w warunkach nieleczonych oraz w 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. Multipleksowanie analiz pozwala zaoszczędzić czas i odczynniki, zwiększa spójność oraz dostarcza dodatkowych informacji na temat badanego układu. Istnieje wiele potencjalnych parametrów fenotypowych, które należy wybierać w zależności od stawianych pytań badawczych. W niniejszej pracy badano zmiany w morfologii komórek i ich żywotności w odpowiedzi na leczenie erlotinibem. Po trzech dniach podawania leku zaobserwowano zmniejszenie powierzchni jądra oraz zwiększenie powierzchni komórek H3255 (Ryc. 3A). Stwierdzono, że średnia różnica w powierzchni jądra między populacjami „bez leku” a „leczonymi” jest statystycznie istotna w dwustronnym teście t typu 2 (równe wariancje) (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 wywołuje efekt cytotoksyczny (t.j. wzrost liczby martwych komórek w czasie) czy cytostatyczny (t.j. spadek liczby nowych komórek w czasie), ponieważ ma to istotny wpływ kliniczny. Na przykład efekt cytostatyczny leku indukuje zatrzymanie wzrostu, ale nie eliminuje komórek z guza, w związku z czym istnieje możliwość ponownego zainicjowania proliferacji komórek nowotworowych po odstawieniu leku. Działanie leków często zależy od kontekstu, stężenia i 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, a zatem nie rozróżniają między zatrzymaniem wzrostu a śmiercią komórkową. W niniejszej pracy martwe komórki zidentyfikowano na podstawie barwienia jodkiem propidylożym (Rycina 3B). Zaobserwowano zarówno efekty cytotoksyczne, jak i cytostatyczne erlotynibu w komórkach H3255, objawiające się zwiększeniem liczby zgonów oraz zmniejszeniem liczby podziałów 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 były pod wpływem leku. Dodatkową zaletą tej platformy jest generowanie danych ilościowych. Na przykład w naszym eksperymencie z ko-kulturą, początkowa subpopulacja 1 18 (75,8%) komórek H325 wzrosła do 2 817 (87,9%) lub spadła do 396 (57,2%) po trzech dniach odpowiednio bez lub z zastosowaniem erlotynibu (Rycina 4). Ponieważ jesteśmy w stanie uzyskać rzeczywistą liczbę komórek zamiast względnego procentu (jak w metodach cytometrii przepływowej), wyciągamy wniosek, że zmiana składu podczas leczenia lekiem wynika ze spadku liczby komórek H325, a nie ze wzrostu liczby komórek CCD-19Lu. Należy zauważyć, że wskaźniki śmiertelności 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.

Rysunek 1: Przegląd protokołu analizy obrazu. Dwa potencjalne potoki analizy obrazu w dalszych etapach procesu, służące do klasyfikacji populacji heterokomórkowych przy użyciu klasyfikacji opartej na morfologii lub klasyfikacji opartej na fluorescencji. Skala = 100 µm. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 2: Zgodność między klasyfikacją opartą na morfologii a klasyfikacją opartą na fluorescencji. (A) Wykres zgodności obrazujący nakładanie się dwóch protokołów klasyfikacji. Te same komórki zostały sklasyfikowane jako H325 przy użyciu klasyfikacji opartej zarówno na morfologii, jak i na fluorescencji. Oba protokoły były zgodne w przypadku 97,4% (n = 1403) komórek nieleczonych oraz 92,5% (n = 916) komórek leczonych erlotinibem (Uwaga: biały obszar jest zbyt mały, aby był widoczny). (BObrazy 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 pomiędzy platformami. Obraz wejściowy: niebieski – jądra komórkowe (Hoechst); zielony – CCD19Lu (GFP). Obrazy klasyfikacji: czerwony – H325; zielony – CCD-19Lu. Pasek skali = 100 µm. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 3: Multiple pomiary fenotypowe z pojedynczego układu eksperymentalnego. (ACechy morfologiczne, takie jak jądra i powierzchnia komórki, obliczono na poziomie pojedynczych komórek w obecności i pod nieobecność leku. Uwaga: powierzchnie komórek mniejsze niż 100 µm2 uznano za zanieczyszczenia i wykluczono z analiz. Wykres pudełkowy przedstawia medianę wraz 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) poddano ko-kulturze, a martwe komórki zidentyfikowano na podstawie intensywności barwienia jodkiem propidium (czerwone) i obrazowano przy użyciu obiektywu 10X. Pasek skali = 1 mm (panel górny); 100 µm (dolne zdjęcia). (CObliczono całkowitą liczbę żywych i martwych komórek w ciągu trzech dni, z zastosowaniem leku lub bez niego, wykazując wyraźny spadek liczby żywych komórek oraz wzrost liczby martwych komórek po dodaniu erlotynibu. Słupki błędów reprezentują błąd standardowy średniej obliczony na podstawie trzech powtórzeń. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej figury.

Rysunek 4: Dynamika subpopulacji w czasie. Reprezentatywne obrazy 10X studni zawierających komórki H325 (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 w panelu skrajnie lewym), 1 mm (obraz górny), 100 µm (dolne obrazy). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.