W celu zbadania regulatorów epigenetycznych chemorezystencji w TNBC przedstawiono ogólny protokół z wykorzystaniem linii komórkowej SUM159PT na Rysunku 1. Najpierw określono wartości IC10 i IC50 Taxolu w komórkach SUM159PT poprzez traktowanie komórek rosnącymi stężeniami Taxolu (Rysunek 2A). Do wytworzenia opornych na Taxol komórek SUM159PT zastosowano metodę stopniowego zwiększania dawki. Po kilku miesiącach przesunięcie wartości IC50 obrazuje pomyślne uzyskanie fenotypu opornego przy użyciu tego protokołu w porównaniu z komórkami kontrolnymi, które były pasażowane równolegle bez ekspozycji na Taxol (Rysunek 2B).

Rysunek 1: Schemat screeningu EPIKOL w komórkach SUM159PT opornych na Taxol. Ogólny przebieg screeningu EPIKOL w wygenerowanych komórkach SUM159PT opornych na Taxol. Rysunek został stworzony za pomocą BioRender.com Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Po uzyskaniu linii komórek opornych przeprowadzono szereg eksperymentów w celu scharakteryzowania ich zmian fenotypowych i molekularnych. Zdolność do proliferacji komórek opornych podczas leczenia Taxolem oceniono za pomocą testu tworzenia kolonii w porównaniu z komórkami wyjściowymi. Jak pokazano na Rysunku 2C,D, komórki oporne zachowują zdolność do proliferacji podczas leczenia Taxolem, podczas gdy komórki wyjściowe nie są w stanie przetrwać. Test ten jest jednym z istotnych etapów charakteryzacji komórek opornych i demonstruje, w jaki sposób można ocenić zdolność do proliferacji pod wpływem leku w porównaniu między komórkami wyjściowymi a opornymi.
W ramach charakterystyki komórek opornych przeprowadzono sekwencjonowanie RNA komórek rodzicielskich i opornych, aby przeanalizować zmiany transkrypcyjne związane z chemoopornością (Rycina 2E). Wśród genów o najwyższym poziomie nadekspresji w komórkach opornych znalazł się ABCB1, a za nim kilka innych transporterów ABC, w tym ABCB4 i ABCA5. Nadekspresja ABCB1 została potwierdzona za pomocą ilościowej reakcji łańcuchowej polimerazy (qPCR), co obrazuje sposób walidacji zmian w ekspresji genów kandydackich z wykorzystaniem qPCR (Rycina 2F).

Rycina 2: Generowanie komórek opornych oraz ich charakterystyka. (A) Schemat przedstawiający generowanie komórek SUM159PT opornych na Taxol metodą stopniowego zwiększania dawki. (B) Krzywa zależności odpowiedzi od dawki wykazująca różnicę między wartościami IC50 dla komórek rodzicielskich a wygenerowanych komórek opornych. IC50 komórek rodzicielskich = 6 nM, a IC50 komórek opornych na Taxol = 355,4 nM. (C) Test tworzenia kolonii w obecności Taxolu. (D) Ilościowe oznaczenie powierzchni kolonii z punktu (C). Zastosowano dwuczynnikową analizę wariancji (ANOVA) z testem post hoc Tukeya przy użyciu programu Prism 8 (GraphPad Software). Poziomy istotności określono jako *P < 0,05, **P < 0,01, ***P < 0,001. (E) Wykresy wulkaniczne przedstawiające geny o zróżnicowanej ekspresji w komórkach Taxol-Res w porównaniu z komórkami rodzicielskimi, przy odcięciu LFC>2 i p<0,001. (F) Walidacja qPCR ekspresji mRNA ABCB1 w komórkach rodzicielskich i Taxol-Res. Wartości P określone za pomocą obustronnego testu t Studenta w porównaniu z grupą kontrolną; **P < 0,01. Panel A został stworzony za pomocą BioRender.com. Rysunek w panelu E został zmodyfikowany za zgodą autorów Yedier-Bayram i wsp.17. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Po wykazaniu nadekspresji ABCB1 oraz innych transporterów ABC w komórkach opornych na Taxol i biorąc pod uwagę szerokie spektrum substratów ABCB1, oceniono występowanie fenotypu MDR. W tym celu komórki oporne na Taxol poddano działaniu wielu leków, które są znanymi substratami transportera ABCB1 i posiadają różne mechanizmy działania, a następnie oceniono ich wrażliwość na te leki (Rysunek 3A,B). Zgodnie z zaobserwowanym profilem ekspresji transporterów ABC, komórki oporne na Taxol wykazały również oporność na doksorubicynę i winkrystynę, co potwierdziło obecność profilu MDR i zilustrowało sposób funkcjonalnej oceny fenotypu wielolekooporności.
Fenotyp MDR stanowi wyzwanie dla metody przesiewowej opartej na CRISPR, EPIKOL, ponieważ biblioteka wymaga selekcji puromycyną w celu wzbogacenia populacji komórek transdukowanych. Ze względu na to, że puromycyna może być usuwana z komórek poprzez ef luks leków zależny od ABCB1, selekcja puromycyną może być nieefektywna w komórkach opornych (Rysunek 3C). W konsekwencji komórki nietransdukowane mogą przetrwać selekcję, co zaburza wyniki przesiewu.
Aby wyeliminować to ograniczenie, przed przesiewaniem EPIKOL zoptymalizowano proces selekcji komórek chemoopornych o profilu MDR. Verapamil, inhibitor ABCB1 pierwszej generacji, zastosowano w celu zablokowania odpływu leków zapośredniczonego przez transportery ABC, co umożliwiło wewnątrzkomórkową akumulację środków chemioterapeutycznych lub puromycyny (Rycyna 3D). Choć sam verapamil nie wpływał na żywotność komórek, przywracał on wrażliwość na leki w połączeniu ze środkami chemioterapeutycznymi (Rycyna 3E,F). Podobnie krzywe przeżywalności po podaniu puromycyny w połączeniu z nietoksycznymi dawkami verapamilu wykazały, że inhibicja ABCB1 przywraca wrażliwość na puromycynę (Rycyna 3G,H). Zapewnia to, że komórki przeżywające selekcję są tymi, które zostały skutecznie transdukowane, a nie takimi, które uniknęły selekcji z powodu odpływu leków napędzanego przez ABCB1. Ustanowienie skutecznego schematu selekcji w komórkach chemoopornych jest kluczowe dla wiarygodnych wyników przesiewania EPIKOL.

Rysunek 3: Fenotyp wielolekowej oporności i wpływ werapamilu na komórki oporne. Krzywe zależności odpowiedzi od dawki pokazujące różnicę w wartościach IC50 między komórkami rodzicielskimi a wygenerowanymi komórkami opornymi w przypadku (A) doksorubicyny, (B) winkrystyny, (C) puromycyny. (D) Schemat przedstawiający selekcję komórek SUM159PT opornych na Taxol z dodatkiem i bez dodatku werapamilu. Rysunek został stworzony za pomocą BioRender.com. (E) Test żywotności komórek Taxol-Res w celu ustalenia optymalnego stężenia werapamilu z Taxolem. (F) Verapamil przywraca wrażliwość komórek opornych na Taxol na działanie Taxolu. (G) Verapamil nie zmienia odpowiedzi komórek rodzicielskich na puromycynę. (H) Verapamil przywraca wrażliwość na puromycynę w komórkach Taxol-Res. Panele A, B i F zostały wykorzystane ponownie za zgodą Yedier-Bayram i wsp.17. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Aby zidentyfikować nowe regulatory chemooporności w TNBC, przeprowadzono przesiew EPIKOL z wykorzystaniem komórek SUM159PT opornych na Taxol, zgodnie z ilustracją na Rysunku 4A. Ponieważ komórki oporne rosną wolniej niż komórki wyjściowe17, podczas całego przesiewu monitorowano poziom podwojenia populacji w grupach traktowanych Taxolem oraz DMSO, aby określić czas potrzebny do osiągnięcia 16 podwojeń populacji (Rysunek 4B). W związku z tym przesiew EPIKOL dla grupy traktowanej DMSO zakończono około 35 dnia, podczas gdy grupę traktowaną Taxolem hodowano przez dodatkowe 15 dni, aby osiągnąć ten sam poziom podwojenia populacji.
Po sekwencjonowaniu wykorzystano potok analizy MAGeCK w celu połączenia efektów sgRNA celujących w ten sam gen, aby ujawnić zmiany na poziomie genu po nokaucie24. Biblioteka EPIKOL zawiera sgRNA celujące w geny niezbędne jako kontrole pozytywne; pierwszym krokiem kontroli jakości była ocena, czy sgRNA te uległy wyczerpaniu zgodnie z oczekiwaniami. Wykresy kaskadowe (waterfall plots) ilustrują zmiany Log2 fold changes genów niezbędnych w różnych porównaniach (Rycina 4C–E). Kiedy próbki z końcowego punktu czasowego z warunków traktowania DMSO lub Taxolem porównano z punktem początkowym, prawie wszystkie geny niezbędne uległy wyczerpaniu po nokaucie, co wskazuje na efektywność badania (Rycina 4C,D). Jednakże to wyczerpanie nie jest wyraźnie obserwowane, gdy dwie próbki końcowe (traktowane DMSO i traktowane Taxolem) porównywano ze sobą, co było spodziewane, ponieważ geny niezbędne są wymagane do przeżycia komórek w obu warunkach (Rycina 4E).
Alternatywną wizualizacją wyników przesiewu jest wykres wulkaniczny (volcano plot), który przedstawia zmiany Log2-krotności (Log2 fold changes) wraz z istotnością statystyczną (wartości p lub FDR). W przesiewie EPIKOL przeprowadzonym w komórkach opornych na Taxol zidentyfikowano wiele genów kontroli pozytywnej jako istotnie uszczuplone, co dodatkowo potwierdza niezawodność przesiewu, ponieważ uszczuplenie sgRNA celujących w geny kontroli pozytywnej stanowi typową kontrolę jakości udanego przesiewu (Rysunek 4F). Oprócz kontroli genów niezbędnych, biblioteka EPIKOL zawiera specyficzne dla kontekstu sgRNA kontroli pozytywnej celujące w kilka genów transporterów ABC. Zgodnie z oczekiwaniami, znormalizowane liczby sgRNA dla ABCB1 w różnych grupach wykazały istotne uszczuplenie sgRNA celujących w ABCB1 specyficznie w grupie traktowanej Taxolem (Rysunek 4G). Ogólnie rzecz biorąc, pokazuje to, w jaki sposób zmiany na poziomie sgRNA mogą być interpretowane w celu identyfikacji genów kandydatów związanych z odpowiedzią na lek.

Rycina 4: Wyniki przesiewu EPIKOL w komórkach opornych na Taxol oraz przykładowe wykresy. (A) Schemat protokołu przesiewowego EPIKOL. (B) Poziom podwojenia populacji (PDL) komórek Taxol-Res. (C–E) Wykresy kaskadowe (waterfall plots) przedstawiające wyniki przesiewu EPIKOL w komórkach Taxol-Res, z zaznaczeniem genów niezbędnych; grupa traktowana DMSO w porównaniu do punktu początkowego w panelu C, grupa traktowana Taxolem w porównaniu do punktu początkowego w panelu D oraz grupa traktowana Taxolem w porównaniu do grupy traktowanej DMSO w panelu E. (F) Wykres wulkaniczny (volcano plot) przedstawiający Log2FoldChanges genów w próbkach traktowanych Taxolem w porównaniu do punktu początkowego. Geny z p<0.05 zostały wyróżnione kolorystycznie i opisane. Rycina została zmodyfikowana na podstawie wcześniejszej publikacji17. (G) Wykresy liczby odczytów na milion (rpm) dla 10 różnych sgRNA ABCB1 z EPIKOL. Panel A został przygotowany za pomocą BioRender.com. Panel E został wykorzystany ponownie za zgodą Yedier-Bayram i wsp.17. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Tabela uzupełniająca 1: Tabela rozwiązywania problemów.Kliknij tutaj, aby pobrać ten plik.
Tabela uzupełniająca 2: Obliczanie miana wirusa i MOI.Prosimy kliknąć tutaj, aby pobrać ten plik.
Tabela uzupełniająca 3: Lista sgRNA biblioteki EPIKOL.Kliknij tutaj, aby pobrać ten plik.
Tabela uzupełniająca 4: Primery użyte podczas przygotowania biblioteki do sekwencjonowania EPIKOL.Prosimy kliknąć tutaj, aby pobrać ten plik.