Artykuł metodologiczny

Dyskryminacja i charakterystyka populacji heterokomórkowych z wykorzystaniem technik obrazowania ilościowego

9.4K wyświetleń

⸱

DOI:

10.3791/55844

⸱

30 czerwca 2017

W tym artykule

Podsumowanie

Opracowaliśmy nowy protokół do badania dynamiki populacji heterokomórkowej w odpowiedzi na zakłócenia. Niniejszy manuskrypt opisuje platformę opartą na obrazowaniu, która tworzy ilościowe zestawy danych do jednoczesnej charakterystyki wielu fenotypów komórkowych populacji heterokomórkowych w solidny sposób.

Streszczenie

Procesy komórkowe są złożone i wynikają z wzajemnego oddziaływania między wieloma typami komórek i ich środowiskiem. Istniejące techniki biologii komórki często nie pozwalają na dokładną interpretację tej wzajemnej zależności. Stosując podejście oparte na obrazowaniu ilościowym, przedstawiamy protokół o wysokiej zawartości do charakteryzowania dynamicznych odpowiedzi fenotypowych (tj. zmian morfologii, proliferacji, apoptozy) heterogenicznych populacji komórek na zmiany bodźców środowiskowych. Podkreślamy naszą zdolność do rozróżniania typów komórek na podstawie intensywności fluorescencji lub nieodłącznych cech morfologii w zależności od zastosowania. Platforma ta pozwala na bardziej kompleksową charakterystykę reakcji subpopulacji na perturbacje przy jednoczesnym wykorzystaniu krótszego czasu, mniejszych ilości odczynników i mniejszego prawdopodobieństwa błędu niż tradycyjne testy biologii komórki. Jednak w niektórych przypadkach populacje komórek mogą być trudne do zidentyfikowania i ilościowego określenia na podstawie złożonych cech komórkowych i będą wymagały dodatkowego rozwiązywania problemów; Zwracamy uwagę na niektóre z tych okoliczności w protokole. Demonstrujemy tę aplikację za pomocą odpowiedzi na lek w modelu raka; Można go jednak z łatwością zastosować szerzej do innych procesów fizjologicznych. Protokół ten pozwala zidentyfikować subpopulacje w systemie kokultur i scharakteryzować szczególną reakcję każdej z nich na bodźce zewnętrzne.

Wprowadzenie

Testy oparte na komórkach były koniem pociągowym w podstawowych badaniach i opracowywaniu leków. Jednak ograniczenia tych standardowych testów stają się coraz bardziej widoczne wraz z rozbieżnością między danymi in vitro a danymi klinicznymi oraz nieuzyskaniem zatwierdzenia przez FDA większości leków. W tym miejscu przedstawiamy nowatorską metodę wykorzystania obrazowania ilościowego do jednoczesnej analizy fenotypów heterokomórkowych w odpowiedzi na odpowiednie, współwystępujące bodźce środowiskowe.

Tradycyjne testy komórkowe, które są używane do pomiaru żywotności komórek, obejmują: testy wykluczenia błękitu trypanowego, MTT/MTS i barwienie cytometrią przepływową Annexin V-FITC. Testy wykluczania błękitu trypanowego, choć proste i niedrogie, wymagają dużej liczby komórek, są czasochłonne i często wpływają na nie stronniczość użytkownika1. Testy MTT i MTS pośrednio mierzą żywotność komórek poprzez pomiary tempa metabolizmu mitochondriów. Jednak na aktywność metaboliczną komórek mogą wpływać różne warunki hodowli (takie jak pożywka lub stężenie tlenu), co prowadzi do niedokładnych wyników i uniemożliwia standaryzację między typami komórek i warunkami2,3. Inną poważną wadą tych technik jest ich niezdolność do rozróżnienia wielu typów komórek – większość systemów biologicznych jest heterokomórkowa. Chociaż metody cytometrii przepływowej mają zdolność rozróżniania wielu populacji komórek, wymagane są etykiety komórek, dynamiczne pobieranie próbek jest trudne, a w przypadku korzystania z przylegających komórek aplikacja ta staje się czasochłonna i podatna na błędy.

Inne ważne fenotypy komórkowe, w tym zmiany morfologiczne, występują w odpowiedzi na bodźce środowiskowe, ale nie są wychwytywane przez tradycyjne testy komórkowe. Profilowanie stanów komórek poprzez charakterystykę morfologiczną i mapowanie podobieństw między próbkami jest potężnym, bezstronnym narzędziem, które może dostarczyć nowatorskich spostrzeżeń w wielu aspektach badań podstawowych i translacyjnych, w tym podstawowej biologii komórki i odkrywania leków4. Ponadto wykazano, że morfologia komórek nowotworowych koreluje z podtypami nowotworów5 i agresywność6. W związku z tym bardzo interesujące jest badanie tych cech komórkowych i ich związku z określonymi zaburzeniami środowiskowymi. Dodatkowo można wykorzystać różnice w cechach morfologicznych do rozróżnienia subpopulacji w systemach kokulturowych. Znakowanie fluorescencyjne komórek ma wady (tj. zmiana wrodzonych właściwości komórek, jest czasochłonne) i dlatego korzystne są dodatkowe metody klasyfikacji typów komórek.

Obrazowanie oparte na mikroskopii jest alternatywną metodą profilowania fenotypów komórkowych w sposób multipleksowany, ilościowy i solidny. W tym manuskrypcie stosujemy nasz potok obrazowania ilościowego, aby podkreślić dynamikę ewolucyjną heterogenicznych populacji komórek w guzie. Koncentrujemy się na interakcji między komórkami niedrobnokomórkowego raka płuca (NSCLC) a fibroblastami związanymi z rakiem (CAF), najbardziej rozpowszechnionym typem komórek zrębu występującym w nowotworach. CAF są zaangażowane w inicjację, progresję i odpowiedź terapeutyczną guza; dlatego wykonywanie testów fenotypowych na komórkach nowotworowych przy braku CAF może być mylące7,8,9. W szczególności oceniliśmy wpływ CAF na komórki nowotworowe w odpowiedzi na erlotynib, małą cząsteczkę ukierunkowaną na receptor naskórkowego czynnika wzrostu (EGFR), który jest często stosowany w leczeniu klinicznym NSCLC. Do oceny wykorzystaliśmy platformę przesiewową o wysokiej zawartości i towarzyszące jej oprogramowanie do analizy obrazu; Próbując udostępnić tę metodologię innym badaczom, opracowaliśmy również porównywalny protokół downstream przy użyciu oprogramowania open source: CellProfiler10 i CellProfiler Analyst11. Większość opartych na obrazach testów przesiewowych o wysokiej zawartości jest analizowana za pomocą skomercjalizowanego oprogramowania specyficznego dla danego modelu urządzenia. Wyniki są trudne do powtórzenia w innych laboratoriach z innym oprogramowaniem, ponieważ podstawowe algorytmy są często zastrzeżone. Korzystając z tego potoku opartego na obrazie, zmierzono proliferację, śmierć i morfologię komórek każdej subpopulacji kultury heterokomórkowej w odpowiedzi na leczenie farmakologiczne przy użyciu klasyfikacji opartej zarówno na fluorescencji, jak i morfologii. Poniższy protokół zapewnia solidną metodologię sondowania złożonych procesów komórkowych.

Protokół

1. Hodowla komórkowa

UWAGA: Tutaj badano fenotypowe odpowiedzi komórek NSCLC (H3255) na czynnik celujący w EGFR, erlotynib, podczas wspólnej hodowli z fibroblastami płuc (CCD-19Lu GFP). Dla tego samego zestawu danych przeprowadzono klasyfikację opartą na fluorescencji i morfologii dwóch populacji komórek (H3255 i CCD-19Lu GFP) w celu zilustrowania ich zgodności. Należy jednak zastosować tylko jedną metodę klasyfikacji, którą należy wybrać na podstawie zastosowania.

  1. Przygotuj 500 ml RPMI-1640 uzupełnionego 10% inaktywowanym termicznie FBS i 1% penicyliny/streptomycyny.
    UWAGA: Pożywkę wzrostową i suplementy można w razie potrzeby wymienić na inne typy komórek.
  2. Hodować komórki w 10 cm3 płytki jako czyste populacje w 5% CO2 w temperaturze 37 °C i pasować w odpowiednim stosunku co 3-4 dni.

2. Przygotowanie komórek

  1. Aby przygotować zawiesinę komórkową, usuń pożywkę komórkową i umyj komórki 5 ml soli fizjologicznej buforowanej fosforanem 1X (PBS).
  2. Odessać z PBS, dodać 1 ml 0,05% trypsyny ogrzanej do 37 °C i inkubować przez 5 minut (lub do momentu, gdy komórki przestaną przylegać do płytki) w temperaturze 37 °C.
  3. Aby zneutralizować trypsynę, dodaj 4 ml RPMI do płytek GFP H3255 i CCD-19Lu i odpipetuj roztwory komórek do stożkowych probówek o pojemności 15 ml.
  4. Odwirować komórki przy 150 x g przez 5 minut w temperaturze pokojowej.
  5. Odessać supernatant i ponownie zawiesić granulki komórek w 10 ml pożywki RPMI w stożkowych probówkach, aby przygotować się do wysiewu.

3. Poszycie komórek

  1. Policz komórki, aby ustandaryzować wysiew.
    1. Odwróć probówki, aby wymieszać komórki w roztworze. Pobrać pipetą 10 μl zawiesiny każdej komórki do probówki mikrowirówkowej i dodać 10 μl błękitu trypanowego.
    2. Odpipetować 10 μl tego roztworu do szkiełka do liczenia komórek i włożyć do automatycznego licznika komórek.
      UWAGA: Można stosować inne metody zliczania komórek (np. hemocytometr).
    3. Powtórz liczbę komórek co najmniej trzy razy lub do momentu, gdy uzyskana liczba będzie spójna.
  2. Komórki nasienne na płytce wielodołkowej
    UWAGA: Liczba płytek do wysiewu zależy od liczby punktów czasowych, które zostaną zobrazowane. Alternatywnie, jedna płytka może być ponownie obrazowana w czasie, jeśli komórki są transdukowane histonem-2B-GFP (lub innymi stabilnymi barwnikami jądrowymi lentiwirusa, które zapewniają odpowiednią segmentację jądrową).
    1. Wysiewaj łącznie 1,500 komórek/studzienkę w 100 μl RPMI. Przygotuj trzy zawiesiny komórkowe do pokrycia trzech populacji komórek: H3255, CCD-19Lu GFP i 50% H3255 + 50% CCD-19Lu GFP. Wykonaj każdy z nich w trzech egzemplarzach z identycznymi ustawieniami płytek dla każdego punktu czasowego.
      UWAGA: Początkowe gęstości wysiewu powinny być zoptymalizowane dla każdej linii komórkowej i rozmiaru płytki.
    2. Komórki nasienne w 3 x 96-dołkowych płytkach za pomocą pipety wielokanałowej. Inkubować komórki O/N w temperaturze 37 °C, 5% CO2 .

4. Dawkowanie leków

  1. Dodać DMSO do proszku erlotynibu, aby uzyskać końcowy roztwór podstawowy erlotynibu o stężeniu 10 mM.
    UWAGA: Lek można zastąpić według uznania.
  2. Rozcieńczyć lek pożywką do hodowli komórkowych do 2x końcowego stężenia o najwyższym stężeniu końcowym przy 10 μM i seryjnie rozcieńczyć czterokrotnie w stosunku 1:10, co daje w sumie pięć stężeń leku i brak kontroli leku.
    UWAGA: Jeśli stężenie DMSO jest równe lub przekracza 0,1% v/v przy najwyższej dawce leku, kontrola bez leku powinna zawierać równoważną ilość DMSO, aby zapewnić, że obserwowane efekty nie są spowodowane toksycznością DMSO.
  3. Odpipetować 100 μl roztworu leku do odpowiedniego dołka, aby uzyskać końcowe stężenie leku 1x rozcieńczonego w pożywce.

5. Akwizycja obrazu

UWAGA: Obrazy zostały pozyskane w dniach 0, 2 i 3. W zależności od typów komórek i badanych procesów komórkowych można badać inne pożądane punkty czasowe.

  1. Wybarwiaj komórki, aby przygotować się do obrazowania.
    1. Uzupełnić roztwór barwnika w następujących końcowych stężeniach.
      1. Do klasyfikacji na podstawie fluorescencji należy przygotować 5 μg/ml barwienia jądrowego i 5 μM barwienia martwych komórek w PBS (patrz tabela materiałów).
      2. Do klasyfikacji opartej na morfologii należy przygotować 5 μg/ml barwienia jądrowego, 5 μM barwienia martwych komórek i 5 μM barwienia komórek w PBS (patrz tabela materiałów).
    2. Dodaj 20 μl roztworu barwnika do każdej studzienki. Inkubować przez 30 minut w temperaturze 37 °C, chroniąc przed światłem.
  2. Optymalizuj i pozyskuj obrazy.
    1. Wyjąć płytkę dołka z inkubatora, przetrzeć spód płytki 70% EtOH i umieścić płytkę w komorze obrazowania.
    2. W zakładce "Konfiguracja" kliknij przycisk "+" w sekcji "Wybór kanału", aby dodać odpowiednie kanały do obrazu (np. jasne pole, plama jądrowa, plama martwych komórek, sygnały GFP i RFP).
    3. Kliknij "Wybór układu" i wykonuj zdjęcia testowe ze stosem z co 2 μm, zaczynając od 0 μm i kończąc na 20 μm, aby zidentyfikować płaszczyznę ostrości. Wprowadź tę odległość dla każdego kanału w polu "Wysokość".
    4. Kliknij "Migawka" pod każdym kanałem, aby ocenić intensywność i zoptymalizować czas ekspozycji. W razie potrzeby wprowadź wyższe lub niższe wartości w polu "Czas".
    5. Po prawej stronie ekranu zaznacz odpowiednie dołki (do zobrazowania) na schemacie tablicy. Na poniższym schemacie studni zaznacz dwadzieścia pięć pól, które mają zostać zobrazowane w podobny sposób.
    6. W sekcji "Uruchom eksperyment" zidentyfikuj tabliczkę w "Nazwa tablicy" po lewej stronie i kliknij przycisk "Start" (poniżej), aby uruchomić protokół akwizycji obrazu z obiektywem 10X w celu wygenerowania obrazów TIFF w skali szarości w wyżej wymienionych kanałach.
      UWAGA: Instrukcje krok po kroku dotyczące protokołu obrazowania mogą się różnić w zależności od instrumentu, ale wartości parametrów powinny być nadal zoptymalizowane.
    7. Powtórz obrazowanie w drugim i trzecim dniu.

6. Analiza obrazu

UWAGA: Wszystkie analizy obrazów zostały przeprowadzone przy użyciu zastrzeżonego oprogramowania. Ponieważ jednak nie jest to publicznie dostępna, porównywalne analizy zaprojektowano również dla CellProfiler 2.2 i CellProfiler Analyst 2.0, z krótkim protokołem wymienionym poniżej i szczegółowym protokołem dostarczonym jako materiał uzupełniający (obrazy testowe są już załadowane do potoków w celu przetestowania przepływu pracy). Obrazy z tego eksperymentu zostały pogrupowane razem dla każdego dołka, tak aby każdy studzienek mógł być ładowany indywidualnie. Poniższe potoki CellProfiler zawierają wyrażenie regularne, które analizuje informacje o metadanych z każdej nazwy pliku obrazu i umożliwia dalsze grupowanie obrazów według kanałów.

  1. Pobierz i zainstaluj oprogramowanie open source: "CellProfiler" i "CellProfiler Analyst". Zaakceptuj wszystkie domyślne opcje i dodatki podczas instalacji.
  2. Klasyfikuj kultury heterokomórkowe do subpopulacji.
    1. Otwórz program 'CellProfiler', kliknij 'Plik'| 'Import Pipeline'| 'From File' i wybierz odpowiedni plik ("Fluorescence_Classification.cppipe" do klasyfikacji opartej na fluorescencji lub "Morphology_Classification.cppipe" do klasyfikacji opartej na morfologii).
      UWAGA: Pliki te zawierają protokoły podstawowego przetwarzania obrazu i segmentacji jąder/komórek. Ze względu na nierównomierne zabarwienie Hoechsta obecne w tych komórkach, jądra zostały podzielone na segmenty i powiększone, a następnie wykorzystane do zamaskowania obrazu, aby umożliwić segmentację jąder o mniejszej intensywności.
      1. Wybierz potok klasyfikacji oparty na fluorescencji, aby sklasyfikować komórki jako H3255 lub CCD-19Lu na podstawie intensywności eGFP i obliczyć cechy morfologii.
        UWAGA: Komórki CCD-19Lu transdukowano za pomocą GFP.
      2. Wybierz potok klasyfikacji oparty na morfologii, aby podzielić jądra i komórki na segmenty oraz obliczyć cechy morfologii.
        UWAGA: Do klasyfikacji potrzebna jest dalsza analiza w programie "CellProfiler Analyst". Martwe komórki są klasyfikowane za pomocą klasyfikacji opartej na fluorescencji.
    2. Kliknij "Wyświetl ustawienia wyjściowe" w lewym dolnym rogu ekranu. Wybierz lokalizację plików graficznych do analizy ("Domyślny folder wejściowy") i miejsce docelowe wyodrębnionych danych ("Domyślny folder wyjściowy").
    3. Kliknij "Analizuj obrazy", aby rozpocząć analizę. Po zakończeniu analizy kliknij "OK" i przejdź do "Domyślnego folderu wyjściowego", aby wyświetlić obliczone dane.
      UWAGA: Przykładowe wartości parametrów i obrazy są dostępne w pliku uzupełniającym 1-CellProfilerProtocol.pdf.
    4. W przypadku klasyfikacji opartej na morfologii (tylko) wykonaj następujące czynności.
      1. Otwórz oprogramowanie 'CellProfiler Analyst', wybierz plik właściwości bazy danych wygenerowany w CellProfiler i kliknij 'Otwórz'.
      2. Kliknij zakładkę "Klasyfikator" w lewym górnym rogu ekranu. Aby wywołać losowe obrazy komórek z eksperymentu, kliknij "Pobierz"; Obrazy pojawią się w oknie Niesklasyfikowane.
    5. Ręcznie sklasyfikuj komórki jako dodatnie (H3255) lub ujemne (CCD-19Lu), zaznaczając i przeciągając komórki do odpowiedniego kosza (patrz Plik uzupełniający 2-Pipeline.pdf: Rysunek B). Po sklasyfikowaniu co najmniej 50 komórek na subpopulację kliknij "Trenuj klasyfikator", a następnie kliknij "Sprawdź postęp".
      UWAGA: Komórki są klasyfikowane za pomocą nadzorowanego przez użytkownika algorytmu losowego uczenia maszynowego lasu12. Jeśli dokładność nie przekracza zatwierdzonego progu (>90%), może być konieczne cofnięcie się i optymalizacja segmentacji komórek w programie "CellProfiler" lub komórki mogą nie być idealne do klasyfikacji według morfologii.
    6. Kliknij "Oceń wszystko", aby wygenerować tabelę z liczbą komórek dla każdej subpopulacji.

Wyniki

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.

Schemat klasyfikacji komórek; mikroskopia, segmentacja, morfologia, analiza fluorescencji.
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.

Analiza efektu leku; schemat przedstawiający klasyfikacje morfologiczne i fluorescencyjne: wyniki zgodności.
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.

Analiza komórek pod wpływem leku: Wykresy pudełkowe jąder i powierzchni komórek (A); mikroskopia fluorescencyjna żywotności komórek z barwieniem martwych komórek (B); wykresy liniowe śledzące wzrost/śmierć komórek (C).
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.

Odpowiedź komórek w hodowli na lek w mikroskopii, obrazy z mikroskopii fluorescencyjnej, wykresy kołowe, H3255, CCD-19Lu.
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.

Dyskusja

Opisany powyżej protokół ulepsza obecne testy biologii komórki, zapewniając bardziej wszechstronny wgląd w dynamikę fenotypową wielu typów komórek w odpowiedzi na zakłócenia środowiskowe przy użyciu zredukowanych odczynników i czasu. Główną zaletą tego eksperymentalnego projektu jest możliwość analizowania wielu fenotypów za pomocą jednej konfiguracji i generowania danych ilościowych charakteryzujących te fenotypy na poziomie pojedynczej komórki. Jedną z technicznych zalet tej platformy jest łatwość początkowego rozwiązywania problemów w porównaniu z innymi testami. Ponieważ ta metoda jest oparta na obrazie, można wizualizować studnie pod kątem pozornego wysiewu nad/niedostatecznego. Wskazane jest, aby komórki znajdowały się w fazie wzrostu wykładniczego przez cały czas trwania eksperymentu i nie były ograniczone przez ograniczenia żywieniowe lub przestrzenne, ani nie były zakłócone przez starzenie się z powodu rzadkiego wysiewu. W przeciwnym razie wskaźniki urodzeń i zgonów mogą nie być powtarzalne między eksperymentami. Dla porównania, CellPD jest publicznie dostępnym programem do obliczania współczynników urodzeń i zgonów14. Dodatkowo można zwizualizować, czy do każdej studzienki dodano odpowiednie stężenia barwników. Problemy z pipetowaniem w pojedynczej studzience mogą powodować błędną segmentację i zniekształcone dane, ale można je łatwo wykryć za pomocą wyżej wymienionego protokołu.

Niestety, nie wszystkie typy komórek mogą być podatne na tę aplikację. Ważne jest, aby móc dokładnie segmentować jądra i komórki, dlatego analiza komórek, które organizują się w bardziej kuliste lub zbite struktury, może nie być odpowiednia. W przypadku niektórych komórek korzystne może być również użycie sitka komórkowego przed wysiewem, aby zapewnić wstępny wysiew pojedynczych komórek. Ponadto technika klasyfikatora liniowego ma zastosowanie tylko do typów komórek, które można łatwo odróżnić na podstawie cech morfologicznych.

Aby protokół był skuteczny, ważne jest, aby najpierw zoptymalizować warunki obrazowania, ponieważ ważność dalszych analiz zależy od jakości obrazów. Podczas gdy w tym przypadku przeprowadziliśmy eksperymenty przy użyciu platformy przesiewowej o wysokiej zawartości, akwizycję obrazu można również przeprowadzić za pomocą dowolnego mikroskopu fluorescencyjnego (chociaż zautomatyzowana platforma obrazowania jest idealna do podejść o wysokiej przepustowości). Obrazy testowe powinny być wykonywane przed każdym punktem czasowym obrazowania, aby upewnić się, że nie ma problemów z mikroskopem lub protokołem. Stosunek sygnału do szumu powinien być wysoki, szczególnie dla kanałów, które będą wykorzystywane do segmentacji (tj. jądra, barwienia komórek). Dodatkowo ważne jest, aby obrazować w optymalnej płaszczyźnie ostrości. Jeśli obrazy są nieostre, segmentacja staje się znacznie trudniejsza, a obliczone cechy morfologiczne będą prawdopodobnie niedokładne. Różnice w oświetleniu między polami mogą również powodować problemy z segmentacją obrazu. Duże różnice w jasności utrudniają automatyczny wybór wartości progowych. Dodatkowo, jeśli między komórkami występują niejednorodne intensywności fluorescencji, pojedyncza wartość progowa może nie wystarczyć do segmentacji wszystkich komórek na obrazie. W tej analizie problemy te zostały przezwyciężone poprzez stworzenie masek wokół jaśniejszych i ciemniejszych populacji komórek oraz segmentację każdej populacji osobno.

Podczas gdy fenotypy badane w tym protokole ograniczają się do charakterystyki żywej, martwej i morfologicznej, można je łatwo rozszerzyć w celu zbadania innych cech. Na przykład funkcjonalne badania genetyczne mogą być dodawane z RNAi, nadekspresją lub innymi zaburzeniami chemicznymi.

W niniejszej pracy zademonstrowano możliwości protokołu do pomiaru odpowiedzi komórek niedrobnokomórkowego raka płuca na erlotynib w obecności i bez CAF. Jest to jednak tylko jeden z wielu typów komórek i parametrów mikrośrodowiskowych, które można przetestować. Rozszerzyliśmy ten protokół, aby można go było stosować z innymi typami komórek i badaniami leków, w tym komórkami pierwotnymi wyizolowanymi z guzów pacjentów13,15.

Oświadczenia

Autorzy nie mają nic do ujawnienia.

Podziękowania

Ta praca została sfinansowana przez National Cancer Institute (NCI) Granty U54CA143798 i U54CA143907 na założenie Centrów Nauk Fizycznych i Onkologii (PS-OCs) odpowiednio w Instytucie Raka Dana-Farber i Uniwersytecie Południowej Kalifornii. S.M. Mumenthaler otrzymał nagrodę PS-OC transnetwork, która wspierała niektóre z tych prac.

Chcielibyśmy wyrazić naszą najgłębszą wdzięczność naszym filantropijnym sympatykom, szczególnie rodzinie Stephensonów, Emmetowi, Toni i Tessie, za ich darowiznę platformy Operetta HCS. Chcielibyśmy również podziękować J. Foo za wskazówki oraz członkom zespołu Centrum Stosowanej Medycyny Molekularnej: D. Agusowi za wskazówki kliniczne i mentoring, K. Patschowi za konstruktywne dyskusje na temat projektu eksperymentalnego, R. Rawatowi za pomoc w protokołach analizy obrazu, J. Katzowi za pomoc techniczną przy Operetce oraz P. Macklinowi i D. Rudermanowi za pomocne dyskusje i informacje zwrotne.

Materiały

Lista materiałów użytych w tym artykule
NazwaFirmaNumer katalogowyKomentarze
RPMI-1640Pożywka do hodowli komórkowychCorning10-040-CV
Płodowa surowica bydlęca (FBS)Gemini100-106średni suplement
Penicylina/streptomycynaGibco15-140-122średni suplement
TC20Licznik komórek1450102 Biorad
Szkiełka z błękitem trypanowymKomórka 1450003 Biorad szkiełka licznikowe
Płytki 96-dołkoweCorning3904przezroczyste dolne czarne płytki
ErlotinibLC LaboratoriesE-4007
DMSOVWR317275-100MLroztwór do ponownego zawieszenia leku
HoechstInvitrogenH21491barwnik jądrowy
Jodek propidyny (PI)InvitrogenP1304MPbarwienie martwych komórek
Cell Tracker Orange CMRALife TechnologiesC34551barwienie całych komórek
System obrazowania o wysokiej zawartości OperetkaPerkin Elmer
CellProfilerBroad Instituewersja 2.2.0 (rev 9969f42)http://cellprofiler.org/releases/
Inkubatordo hodowli komórkowychOdpowiedni będzie każdy inkubator do hodowli komórkowych - komórki hodowano pod 37 & ordm; C przy 5% CO2,
15 mL Probówki stożkowe FalconFalcon14-959-53A
10 cm2 płytki do hodowli komórkowychTPP93040

Bibliografia

  1. Kanemura, Y., et al. Evaluation of in vitro proliferative activity of human fetal neural stem/progenitor cells using indirect measurements of viable cells based on cellular metabolic activity. J Neurosci Res. 69 (6), 869-879 (2002).
  2. Plumb, J. A., Milroy, R., Kaye, S. B. Effects of the pH dependence of 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-tetrazolium bromide-formazan absorption on chemosensitivity determined by a novel tetrazolium-based assay. Cancer Res. 49 (16), 4435-4440 (1989).
  3. Hsu, S., et al. Green tea polyphenols induce differentiation and proliferation in epidermal keratinocytes. J Pharmacol Exp Ther. 306 (1), 29-34 (2003).
  4. Caicedo, J. C., Singh, S., Carpenter, A. E. Applications in image-based profiling of perturbations. Curr Opin Biotechnol. 39, 134-142 (2016).
  5. Sero, J. E., et al. Cell shape and the microenvironment regulate nuclear translocation of NF-kappaB in breast epithelial and tumor cells. Mol Syst Biol. 11 (3), 790(2015).
  6. Chen, J. F., et al. Subclassification of prostate cancer circulating tumor cells by nuclear size reveals very small nuclear circulating tumor cells in patients with visceral metastases. Cancer. 121 (18), 3240-3251 (2015).
  7. Ohlund, D., Elyada, E., Tuveson, D. Fibroblast heterogeneity in the cancer wound. J Exp Med. 211 (8), 1503-1523 (2014).
  8. Paraiso, K. H., Smalley, K. S. Fibroblast-mediated drug resistance in cancer. Biochem Pharmacol. 85 (8), 1033-1041 (2013).
  9. Kalluri, R. The biology and function of fibroblasts in cancer. Nat Rev Cancer. 16 (9), 582-598 (2016).
  10. Carpenter, A. E., et al. CellProfiler: image analysis software for identifying and quantifying cell phenotypes. Genome Biol. 7 (10), R100(2006).
  11. Jones, T. R., et al. CellProfiler Analyst: data exploration and analysis software for complex image-based screens. BMC Bioinformatics. 9, 482(2008).
  12. Dao, D., et al. CellProfiler Analyst: interactive data exploration, analysis and classification of large biological image sets. Bioinformatics. 32 (20), 3210-3212 (2016).
  13. Garvey, C. M., et al. A high-content image-based method for quantitatively studying context-dependent cell population dynamics. Sci Rep. 6, 29752(2016).
  14. Juarez, E. F., et al. Quantifying differences in cell line population dynamics using CellPD. BMC Syst Biol. 10 (1), 92(2016).
  15. Mumenthaler, S. M., et al. The Impact of Microenvironmental Heterogeneity on the Evolution of Drug Resistance in Cancer Cells. Cancer Inform. 14, Suppl 4. 19-31 (2015).

Przedruki i uprawnienia

Tagi

Klasyfikacja komórekIntensywność fluorescencjiCechy morfologiczneProtokół wysokoprzepustowej analizy obrazuBodźce środowiskoweAnaliza odpowiedzi na lekAnaliza w programie CellProfilerCharakterystyka subpopulacji