Testy sferoidów w płytkach 96-dołkowych o ultra-niskim przyczepieniu zapewniają wysokoprzepustową ocenę fenotypową pod kątem potencjalnej onkogenności w kontekście, który lepiej odzwierciedla warunki fizjologiczne występujące w guzach in vivo. W rzeczywistości linie komórek nowotworowych MCF10DCIS.com i BT474 tworzą zwarte struktury sferoidalne (Rycina 1A), a badanie immunohistologiczne przekrojów sferoidów wykazało wyraźne zmiany przestrzenne w morfologii komórkowej i jądrowej. Z czasem niektóre sferoidy, takie jak sferoidy BT474, rozwijają obszary nekrotyczne, co jest częstą cechą agresywnych guzów litych (Rycina 1B). Niektóre sferoidy nie wykształcają rdzeni nekrotycznych, jak w przypadku linii komórkowej MDA-MB-231, ale wykazują znaczną zmienność w ekspresji markera proliferacji Ki67, który koreluje odwrotnie z ekspresją kaspazy-3 w formie rozszczepionej, będącej markerem apoptozy (Rycina 1C). Aby potwierdzić, że wiele linii komórkowych jest podatnych na wyciszanie genów za pomocą siRNA, komórki BT474, MCF10DCIS.com, MDA-MB-231 i JIMT1 poddano odwrotnej transfekcji siRNA przez siedem dni. Obecność odczynnika do transfekcji (mock) lub transfekcja kontrolnych siRNA nie wpłynęła na żywotność sferoidów, podczas gdy wyciszenie kluczowego genu Ubikwiny B (UBB) znacząco zredukowało żywotność sferoidów we wszystkich badanych liniach komórkowych (Rycina 1D).
Zaprojektowaliśmy ludzką bibliotekę siRNA obejmującą najczęściej mutowane geny w nieoselekcjonowanych, ER+, HER2+ oraz potrójnie ujemnych rakach piersi. Biblioteka ta składa się z genów o znanej funkcji, takich jak MYC, PIK3CA i TP53, oraz genów, których udział w kancerogenezie nie został jeszcze ustalony. Przesiew zawiera również kilka niecelujących kontrolnych siRNA (Kontrola nr 1, Kontrola nr 2) oraz siRNA celujących w geny niezbędne do przeżycia, takie jak PLK1 i UBB, które pełnią rolę kontroli śmierci komórkowej (Tabela 1). Wybraliśmy linie komórkowe raka piersi BT474, ponieważ łatwo tworzą one sferoidy bez konieczności dodawania zrekonstytuowanej błony podstawnej, są powszechnie stosowanymi liniami komórkowymi i mają znaną architekturę genomową. Na przykład komórki BT474 są dodatnie pod kątem receptora estrogenowego (ER+), wykazują nadekspresję ludzkiego receptora czynnika wzrostu naskórka 2 (HER2+) oraz posiadają mutacje w genach TP53 (E285K) i PIK3CA (K1N)13.
Stosując opisany powyżej protokół, monitorowaliśmy wpływ wyciszenia genów na wielkość i żywotność sferoidów po siedmiu dniach odwrotnej transfekcji siRNA (Ryc. 1E). Co ciekawe, większość genów nie miała istotnego wpływu na powierzchnię lub żywotność sferoidów (Ryc. 2A). Wyciszenie FOXO3, PIK3CA, ERBB2 i SF3B1 doprowadziło do najbardziej znaczącej, powtarzalnej redukcji wielkości sferoidów. Redukcję tę zaobserwowano również w żywotności sferoidów po wyciszeniu ERBB2 i SF3B1. Z satysfakcją potwierdziliśmy wpływ wyciszenia PIK3CA, ERBB2 i SF3B1 na wielkość sferoidów przy użyciu mikroskopii w polu jasnym (Ryc. 2B). Wcześniej zidentyfikowaliśmy SF3B1 jako gen niezbędny w wielu modelach linii komórkowych, zatem siSF3B1 stanowi dobrą kontrolę pozytywną dla śmierci komórek, obok UBB 14. Co ciekawe, spośród wszystkich 20 genów tylko wyciszenie E-Cadheriny doprowadziło do znaczącego wzrostu żywotności sferoidów (Ryc. 2A). Analiza morfologii sferoidów wykazała, że wyciszenie E-Cadheriny spowodowało całkowity rozpad architektury sferoidu, a żywe komórki osiadły na dnie studni o niskim stopniu adhezji (Ryc. 2B). Ręczna ponowna analiza danych z przesiewu objętości sferoidów wykazała, że zjawisko to również zaobserwowano, lecz zostało ono usunięte z kwantyfikacji powierzchni, ponieważ obiekt przekraczał ustawione ograniczenia wielkości. Jak podkreślono wcześniej, komórki BT474 nadmiernie eksponują receptorową tyrozynową kinazę HER2 i posiadają onkogeniczną mutację w PIK3CA (K1N). Potwierdziliśmy, że wyciszenie ERBB2 i PIK3CA doprowadziło do zmniejszenia żywotności sferoidów, podczas gdy transfekcja kontrolami nietargetującymi nie wywoła żadnego efektu (Ryc. 2C).
Następnie zbadano wpływ wyciszenia siRNA na histologię sferoidów. Sferoidy BT474 poddano odwrotnej transfekcji kontrolnym siRNA niecelującym oraz siRNA celującymi w PIK3CA, ERBB2 i UBB. Wyciszenie ERBB2 i UBB doprowadziło do zmniejszenia ekspresji markera proproliferacyjnego Ki67 w porównaniu z kontrolnym siRNA (Rysunek 2D). Aktywacja proapoptotycznego markera rozszczepionej kaspazy-3 została zaobserwowana wyłącznie po wyciszeniu UBB, co sugeruje, że usunięcie HER2 i PIK3CA nie spowodowało apoptozy, lecz działało cytostatycznie, a nie cytotoksycznie. Rzeczywiście, wyciszenie HER2 i PIK3CA doprowadziło do zwiększenia ekspresji białka p27, odpowiedzialnego za zatrzymanie cyklu komórkowego, w porównaniu ze sferoidami poddanymi transfekcji kontrolnej.
Podsumowując, wyniki te wskazują, że wzrost komórek BT474 w formie sferoidów 3D jest napędzany przez onkogenne szlaki sygnałowe HER2 oraz PIK3CA. Co ważniejsze, wyniki te dowodzą, że możliwe jest odporne i powtarzalne zaprojektowanie oraz wdrożenie spersonalizowanej biblioteki przesiewowej siRNA obejmującej setki genów w sferoidach linii komórkowych raka.

Rysunek 1: Optymalizacja linii komórkowych do wzrostu trójwymiarowego. A. Linie komórkowe raka piersi, MCF10DCIS.com i BT474, hodowano na płytkach o niskim stopniu adhezji przez 7 dni. Reprezentatywne obrazy w jasnym polu wykonano za pomocą mikroskopu odwróconego. Paski skali = 100 µm. B. Sferoidy MCF10DCIS.com i BT474 hodowano przez 28 dni. 100 μL uzupełniano świeże medium co 3–4 dni. Sferoidy utrwalono w 3,8% formaldehydzie, zakastowano, wykonano ich przekroje i wybarwiono hematoksyliną i eozyną (H&E). Przedstawiono reprezentatywne obrazy przy małym i dużym powiększeniu. Paski skali odpowiadają 100 µm i 33 µmodpowiednio. C. Sferoidy MDA-MB-231 hodowano przez 21 dni. 100 μL uzupełniano świeżą pożywką co 3-4 dni. Sferoidy utrwalono w 3,8% formaldehydzie, zatopiono w bloku, wycięto preparaty i zabarwiono na Ki67 oraz rozszczepioną kaspazę-3. Przedstawiono reprezentatywne obrazy. Paski skali = 100 µm. D. Linie komórkowe BT474, MCF10DCIS.com, MDA-MB-231 i JIMT1 poddano odwrotnej transfekcji za pomocą mock (sam odczyn do transfekcji) oraz kontrolnych siRNA oraz siRNA przeciwko UBB, a następnie płytki ultra-low attachment wirowano w celu utworzenia sferoidów. Świeżą pożywkę (100 µl) dodawano w dniach 1. i 4. Po 7 dniach określono żywotność komórek. Dane przedstawiają średnią ± SD z dwóch niezależnych powtórzeń biologicznych wykonanych w trzech kopiach, znormalizowanych względem kontroli nr 1. Istotność statystyczną obliczono za pomocą nieparzystego testu t Studenta (p< 0.05). E. Schemat blokowy podsumowujący protokół odwrotnej transfekcji wykorzystywany do funkcjonalnej analizy zależności genowej w sferoidach linii komórkowych nowotworowych. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 2: Funkcjonalne badania genomiczne sferoidów BT474 ujawniają zależności onkogenne. A. Komórki BT474 poddano odwrotnej transfekcji z wykorzystaniem biblioteki siRNA siGENOME dla 20 genów ludzkich w trzech powtórzeniach. Obserwowano wielkość i żywotność sferoidów. Surowe wartości danych znormalizowano względem mediany płytki, a następnie obliczono wartości z-score w celu zidentyfikowania istotnych wartości odstających, przekraczających 1,7-krotność odchylenia standardowego mediany płytki. 15Należy zwrócić uwagę na geny odstające: ERBB2, SF3B1, PLK1, UBB oraz E-cad. siRNA, które znacząco zwiększyły lub zmniejszyły powierzchnię i żywotność sferoidów, są zaznaczone kolorem niebieskim i czerwonym, przy czym kolor czerwony oznacza kontrolne siRNAs. B. Komórki poddano transfekcji odwrotnej za pomocą niespecyficznego siRNA kontrolnego lub siRNA przeciwko E-kadherynie (E-Cad), PIK3CA, ERBB2 lub UBB, a następnie wirowano w płytce o ultraniskim stopniu przylegania w celu utworzenia sferoidów. Po 7 dniach wykonano reprezentatywne obrazy sferoidów w jasnym polu za pomocą mikroskopu odwróconego. Paski skali = 100 µm. C. Komórki poddano odwrotnej transfekcji za pomocą niecelującego siRNA kontrolnego lub siRNA przeciwko PIK3CA, ERBB2 lub UBB, a następnie odwirowano w płytce o ultra-niskiej adhezji w celu utworzenia sferoidów. Po 7 dniach określono przeżywalność komórek. Dane przedstawiają średnią ± SD z dwóch niezależnych powtórzeń biologicznych wykonanych w trzech powtórzeniach technicznych, znormalizowanych względem kontroli nr 1. Istotność statystyczną obliczono przy użyciu nieparzystego testu t Studenta (p < 0.05). D. Po 7 dniach sferoidy utrwalono, zatopiono, pocięto na skrawki i zabarwiono na H&E.&E, Ki67, rozszczepiona kaspaza-3 i p27. Przedstawiono obrazy reprezentatywne. Paski skali = 100 µmZwróć uwagę na H&Sferoidy siControl wydają się mniejsze ze względu na artefakt przetwarzania oraz wybór poszczególnych nienaruszonych sferoidów do barwienia. Nie wpływa to jednak negatywnie na zmiany zaobserwowane na obrazach zdigitalizowanych oraz wyniki przeżywalności komórek. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej figury.
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 1 | 12 |
| | | | | | | | | | | |
| Kontrolny niecelujący 1 | TP53 | DST | KMT2C | Nietraktowane 1 | FCGBP | ARID1B | FBXM7 | TTC40 | Kontrolny niecelujący 2 | |
| GATA3 | TTN | MUC12 | MUC4 | AHNAK | HUWEI1 | DNAH1 | ITPR2 | ABCA13 | CREBBP | |
| MAP2K4 | PIK3CA | F5 | APOB | ANKRD30A | MUC17 | DNAH17 | LAMA2 | ACE | CSMD2 | |
| STARD9 | USH2A | FAT3 | LPR2 | CSMD3 | MYO18B | DNAH5 | MDN1 | ARHGAP5 | DNAH9 | |
| CXCR3 | MUC16 | RB1 | PKHD1L1 | DNAH2 | SYNE2 | DYNC1H1 | PCLO | CACNA1B | ERBB2 | |
| PLK1 | SYNE1 | LYST | PTEN | SPTA1 | Nietraktowane 2 | VHL | RYR1 | COL6A3 | UBB | |
| | | | | | | | | | | |
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 1 | 12 |
| | | | | | | | | | | |
| Kontrolny niecelujący 1 | ENAM | RYR2 | BRCA2 | Nietraktowane 1 | FMN2 | HECW1 | LAMB4 | SI | Kontrolny niecelujący 2 | |
| FHOD3 | MACF1 | RYR3 | C2ORF16 | DMD | FRG1 | HERC2 | MYH1 | STAB1 | ZDBF2 | |
| GOLGA6L2 | NEB | SMG1 | CACNA1E | DNA14 | GCC2 | HIVEP2 | NIPBL | TANC1 | ZNF536 | |
| HMCN1 | NF1 | UBR5 | CACNA1F | DYNC2H1 | GON4L | HYDIN | PKD1L1 | TF | ANK3 | |
| HRNR | OBSCN | USP34 | CYMA5 | FAM208B | GPR112 | ITSN2 | RNF213 | TPR | ASPM | |
| PLK1 | PCDH15 | XIRP2 | COL7A1 | FLG2 | Nietraktowane 2 | VHL | SAGE1 | UNC80 | UBB | |
| | | | | | | | | | | |
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 1 | 12 |
| | | | | | | | | | | |
| Kontrolny niecelujący 1 | DCHS2 | MUC5B | ZFHX4 | Nietraktowane 1 | MYO9A | SPHKAP | CXORF22 | NCOR1 | VPS13D | |
| ATR | DMXL2 | MXRA5 | ANK2 | KIAA1210 | PRUNE2 | TCHH | DANH6 | NOTCH2 | ANKRD12 | |
| BIRC6 | DNAH10 | TENM1 | ATM | LRP1 | SCN10A | VPS13C | DNAH7 | SPEN | C5ORF42 | |
| CDH1 | DNAH3 | PEG3 | DIDO1 | MAP1A | SCN2A | IPTR3 | ERBB3 | SRRM2 | CCDC8A | |
| CUBN | NOCK1 | RELN | DNAH8 | MED12 | SHROOM2 | CEP350 | FAT4 | SZT2 | CHD4 | |
| PLK1 | EYS | SACS | KIAA109 | MED13 | Nietraktowane 2 | VHL | KMT2A | VPS13A | UBB | |
| | | | | | | | | | | |
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 1 | 12 |
| | | | | | | | | | | |
| Kontrolny niecelujący | QSER1 | ARID1A | WDFY3 | Nietraktowane 1 | SDK1 | TEX15 | LAMA1 | | Kontrolny niecelujący 2 | |
| COL14A1 | SHROOM3 | ATRX | EFCAB5 | SF3B1 | CBFB | AHNAK2 | | | | |
| CSMD1 | TBX3 | KIAA0947 | FOXA1 | ITPR1 | DDX3X | KIF4A | | | | |
| MEFV | UBR4 | MYCBP2 | INPPL1 | FLG | HECTD4 | FAT2 | | | | |
| MGAM | VCAN | NBEAL1 | MAP3K1 | AKAP9 | GPR98 | FOXO3 | | | | |
| PLK1 | ZNF462 | SETX | NRP1 | HERC1 | Nietraktowane 2 | VHL | | | UBB | |
| | | | | | | | | | | |
Tabela 1: Układ płytki i dekonwolucja ludzkiej biblioteki siRNA.Tabela zawiera zawartość każdej z pul siRNA oraz układ płytek o niskim stopniu adhezji wykorzystanych w badaniu przesiewowym.