Ten manuskrypt opisuje metodę badań przesiewowych na skalę genomu opartą na komórkach w celu identyfikacji interakcji receptor-ligand zewnątrzkomórkowy.
Artykuł metodologiczny
Ten manuskrypt opisuje metodę badań przesiewowych na skalę genomu opartą na komórkach w celu identyfikacji interakcji receptor-ligand zewnątrzkomórkowy.
Komunikacja międzykomórkowa, w której pośredniczą bezpośrednie interakcje między receptorami na powierzchni komórki osadzonej w błonie, jest kluczowa dla prawidłowego rozwoju i funkcjonowania organizmów wielokomórkowych. Wykrywanie tych interakcji pozostaje jednak wyzwaniem technicznym. Niniejszy manuskrypt opisuje systematyczne podejście do genetycznego nokautowania CRISPR/Cas9 w skali genomu, które ujawnia szlaki komórkowe wymagane dla określonych zdarzeń rozpoznawania powierzchni komórki. Test ten wykorzystuje rekombinowane białka wytwarzane w systemie ekspresji białek ssaków jako zapalone sondy wiążące do identyfikacji partnerów interakcji w genetycznym badaniu przesiewowym opartym na komórkach. Metoda ta może być wykorzystana do identyfikacji genów niezbędnych do interakcji na powierzchni komórki wykrywanych przez rekombinowane sondy wiążące odpowiadające ektodomenom receptorów osadzonych w błonie. Co ważne, biorąc pod uwagę charakter tego podejścia w skali genomu, ma ono również tę zaletę, że nie tylko identyfikuje bezpośredni receptor, ale także komponenty komórkowe, które są wymagane do prezentacji receptora na powierzchni komórki, dostarczając w ten sposób cennych informacji na temat biologii receptora.
Zewnętrzne interakcje białek receptorowych na powierzchni komórki kierują ważnymi procesami biologicznymi, takimi jak organizacja tkanek, rozpoznawanie gospodarza-patogenu i regulacja immunologiczna. Badanie tych interakcji jest interesujące dla szerszej społeczności biomedycznej, ponieważ receptory błonowe są celami systematycznie dostarczanych środków terapeutycznych, takich jak przeciwciała monoklonalne. Pomimo ich znaczenia, badanie tych interakcji pozostaje wyzwaniem technicznym. Dzieje się tak głównie dlatego, że receptory osadzone w błonie są amfipatyczne, co utrudnia ich wyizolowanie od błon biologicznych w celu manipulacji biochemicznej, a ich interakcje charakteryzują się słabym powinowactwem do interakcji (KDs w zakresie μM-mM)1. W związku z tym wiele powszechnie stosowanych metod nie nadaje się do wykrywania tej klasy interakcji białkowych1,2.
Opracowano szereg metod w celu szczegółowego badania interakcji receptor-ligand zewnątrzkomórkowy, które uwzględniają ich unikalne właściwości biochemiczne3. Wiele z tych podejść obejmuje ekspresję całej ektodomeny receptora jako rozpuszczalnego białka rekombinowanego w systemach opartych na komórkach ssaków lub owadów, aby zapewnić, że białka te zawierają modyfikacje potranslacyjne, które są ważne strukturalnie, takie jak glikany i wiązania dwusiarczkowe. Aby przezwyciężyć wiązanie o niskim powinowactwie, ektodomeny są często oligomeryzowane, aby zwiększyć ich awidność wiązania. Ektodomeny białkowe Avid zostały z powodzeniem wykorzystane jako sondy wiążące do identyfikacji partnerów interakcji w bezpośrednich rekombinowanych ekranach interakcji białko-białko4,5,6,7. Chociaż metody oparte na białkach rekombinowanych są szeroko skuteczne, wymagają, aby ektodomena receptora błonowego była wytwarzana jako białko rozpuszczalne. W związku z tym ma ogólne zastosowanie tylko do białek, które zawierają ciągły region zewnątrzkomórkowy (np. jednoprzebiegowe typu I, typu II lub zakotwiczone w GPI) i nie jest ogólnie odpowiedni dla kompleksów receptorowych i białek błonowych, które wielokrotnie obejmują błonę.
Techniki klonowania wyrażeń, w których biblioteka komplementarnych DNA (cDNA) jest transfekowana do komórek i testowana pod kątem fenotypu wzmocnienia wiązania, zostały również wykorzystane do identyfikacji zewnątrzkomórkowych interakcji białko-białko8. Dostępność dużych zbiorów sklonowanych i zsekwencjonowanych plazmidów ekspresji cDNA w ostatnich latach ułatwiła metody, w których linie komórkowe z nadekspresją cDNA kodujące receptory na powierzchni komórki są badane pod kątem wiązania białek rekombinowanych w celu identyfikacji interakcji9,10. Podejścia oparte na nadekspresji cDNA, w przeciwieństwie do metod opartych na rekombinacji białek, dają możliwość identyfikacji interakcji w kontekście błony komórkowej. Jednak sukces stosowania konstruktów ekspresji cDNA zależy od zdolności komórek do nadekspresji białka w prawidłowo sfałdowanej formie, ale często wymaga to komórkowych czynników dodatkowych, takich jak transportery, białka opiekuńcze i prawidłowe składanie oligomeryczne. Transfekcja pojedynczego cDNA może zatem nie wystarczyć do osiągnięcia ekspresji na powierzchni komórki.
Techniki przesiewowe wykorzystujące konstrukty cDNA lub rekombinowane sondy białkowe są zasobożerne i wymagają dużych zbiorów cDNA lub bibliotek białek rekombinowanych. Specjalnie zaprojektowane metody oparte na spektrometrii mas zostały ostatnio wykorzystane do identyfikacji oddziaływań zewnątrzkomórkowych, które nie wymagają składania dużych bibliotek. Techniki te wymagają jednak chemicznej manipulacji powierzchnią komórki, która może zmienić biochemiczny charakter cząsteczek obecnych na powierzchni komórek i obecnie ma zastosowanie tylko w przypadku interakcji, w których pośredniczą białka glikozylowane11,12. Większość obecnie dostępnych metod również w dużym stopniu koncentruje się na interakcjach między białkami, w dużej mierze ignorując wkład mikrośrodowiska błonowego, w tym cząsteczek takich jak glikany, lipidy i cholesterol.
Niedawny rozwój wysoce efektywnego celowania bialleleicznego przy użyciu metod opartych na CRISPR umożliwił biblioteki w skali genomu komórek pozbawionych zdefiniowanych genów w jednej puli, które mogą być badane w systematyczny i bezstronny sposób w celu identyfikacji komponentów komórkowych zaangażowanych w różne konteksty, w tym analizowanie procesów sygnalizacji komórkowej, identyfikacja zaburzeń nadających odporność na leki, toksyny, i patogenów oraz określanie swoistości przeciwciał13,14,15,16. W tym miejscu opisujemy test przesiewowy komórek nokautujących oparty na CRISPR w skali genomu, który stanowi alternatywę dla obecnych podejść biochemicznych do identyfikacji interakcji receptor-ligand zewnątrzkomórkowy. To podejście polegające na identyfikowaniu interakcji, w których pośredniczą receptory błonowe, za pomocą genetycznych badań przesiewowych jest szczególnie odpowiednie dla badaczy, którzy koncentrują się na pojedynczych ligandach, ponieważ pozwala uniknąć konieczności kompilowania dużych bibliotek cDNA lub rekombinowanych białek.
Ten test składa się z trzech głównych etapów: 1) Wysoce zapalone sondy wiążące białka rekombinowane, składające się z zewnątrzkomórkowych regionów receptora będącego przedmiotem zainteresowania, są produkowane i używane w testach wiązania opartych na cytometrii przepływowej opartej na fluorescencji; 2) Testy wiązania służą do identyfikacji linii komórkowej, która wyraża partnera interakcji rekombinowanej sondy białkowej; 3) Wytwarzana jest wersja linii komórkowej z ekspresją Cas9, która oddziałuje z białkiem będącym przedmiotem zainteresowania, a następnie przeprowadza się badanie przesiewowe nokautu oparte na CRISPR / Cas9 w skali genomu (Figura 1). W tym genetycznym badaniu przesiewowym wiązanie rekombinowanego białka z liniami komórkowymi jest wykorzystywane jako mierzalny fenotyp, w którym komórki w bibliotece nokautu, które utraciły zdolność wiązania sondy, są sortowane za pomocą opartego na fluorescencji aktywnego sortowania komórek (FACS) i genów, które spowodowały utratę fenotypu wiązania zidentyfikowanego przez sekwencjonowanie. W zasadzie identyfikowane są geny kodujące receptor odpowiedzialny za wiązanie sondy Avid oraz te, które są niezbędne do wyświetlania jej powierzchni komórki.
Pierwszy krok tego protokołu polega na produkcji rekombinowanych sond białkowych reprezentujących ektodomenę receptorów związanych z błoną. Wiadomo, że receptory te często zachowują swoje zewnątrzkomórkowe funkcje wiązania, gdy ich ektodomeny są wyrażane jako rekombinowane rozpuszczalne białko1. W przypadku białka będącego przedmiotem zainteresowania, rozpuszczalne białka rekombinowane mogą być wytwarzane w dowolnym odpowiednim systemie ekspresji białek eukariotycznych w dowolnym formacie, pod warunkiem, że mogą być oligomeryzowane w celu zwiększenia awidności wiązania i zawierają znaczniki, które można stosować w testach wiązania opartych na cytometrii przepływowej opartej na fluorescencji (np. znacznik FLAG, znacznik biotyny). Szczegółowe protokoły wytwarzania rozpuszczalnych ektodomen receptorów błonowych z wykorzystaniem systemu ekspresji białek HEK293, a także różne techniki multimeryzacji i konstrukty ekspresji białek do produkcji zarówno białek pentamerycznych, jak i białek monomerycznych zostały wcześniej opisane1,17. Protokół tutaj opisze etapy generowania fluorescencyjnych sond avid z monomerycznych biotynylowanych białek poprzez sprzężenie ich ze streptawidyną sprzężoną z fluorochromem (np. fikoerytryną lub PE), która może być stosowana bezpośrednio w testach wiązania opartych na komórkach i ma tę zaletę, że nie wymaga wtórnego przeciwciała do wykrywania. Ogólne protokoły wykonywania badań przesiewowych w skali genomu zostały już opisane20,21, dlatego protokół koncentruje się głównie na specyfice wykonywania badań przesiewowych wiązania białek rekombinowanych opartych na cytometrii przepływowej przy użyciu systemu badań przesiewowych CRISPR/Cas9 knockout przy użyciu 18.
1. Produkcja i oczyszczanie biotynylowanych białek znakowanych Hisem
2. Kwantyfikacja i oligomeryzacja monomerycznego białka biotynylowanego
UWAGA: Aby zwiększyć awidność wiązania, oligomeryzuj biotynylowane białka monomeryczne na tetramerycznej streptawidynie-PE przed użyciem ich w testach wiązania. Osiągnąć optymalne współczynniki koniugacji białek monomerycznych i tetramerycznej streptawidyny-PE poprzez badanie serii rozcieńczeń biotynylowanych monomerów w stosunku do ustalonego stężenia streptawidyny i empiryczne ustalenie minimalnego rozcieńczenia, przy którym nie można wykryć nadmiaru biotynylowanych monomerów.
3. Testy wiązania komórek oparte na cytometrii przepływowej
4. Określanie wkładu wiązania epitopów termolabilnych i łańcuchów bocznych siarczanu heparanu
UWAGA: Aktywność wielu białek jest niewymierna, więc utrata aktywności wiązania po obróbce cieplnej jest zachęcająca. Zaleca się określenie udziału ujemnie naładowanych glikozaminoglikanów, głównie siarczanu heparanu (HS), w pośredniczeniu w wiązaniu rekombinowanych białek. Dzieje się tak dlatego, że wiązanie przez HS w opisanym tutaj teście wiązania komórkowego może być addytywne, a nie współzależne od innych receptorów19. Oznacza to, że obserwowane wiązanie może być całkowicie pośredniczone przez łańcuchy boczne HS proteoglikanów na powierzchni komórki, a nie przez specyficzny receptor. Wiązanie z HS na powierzchni komórki niekoniecznie jest niespecyficzne, ale raczej jest właściwością białka, którą warto znać przed wykonaniem pełnego badania genetycznego.
5. Wybór linii komórkowych stabilnie wyrażających Cas9
UWAGA: Zanim linia komórkowa wiążąca sondę będzie mogła zostać wykorzystana w badaniach przesiewowych CRISPR, musi najpierw zostać zaprojektowana tak, aby wyrażała nukleazę Cas9 i wysoce aktywny klon selected19.
6. Wybieranie klonów o wysokiej aktywności Cas9
UWAGA: Poliklonalny Cas9 może być używany do skutecznego przeprowadzania badań genetycznych; jednak wybór klonu o wysokiej aktywności Cas9 poprawia wyniki badań przesiewowych18.
7. Generowanie biblioteki nokautów przesiewowych CRISPR-Cas9 dla całego genomu
8. Genetyczne badania przesiewowe pod kątem wiązania powierzchni komórek
9. Ekstrakcja genomowego DNA i pierwszy PCR do wzbogacenia gRNA
10. Druga runda PCR do indeksowania i sekwencjonowania
11. Analiza bioinformatyczna w celu identyfikacji receptora i powiązanych szlaków
Dane z sekwencjonowania pochodzące z dwóch reprezentatywnych przesiewów knockout w skali genomu, przeprowadzonych odpowiednio w komórkach NCI-SNU-1 i HEK293 w celu identyfikacji partnera wiązania ludzkiego TNFSF9 oraz P. falciparum RH5, zostały udostępnione (Tabela uzupełniająca 1). Zachowanie wiązania RH5 było zależne zarówno od siarczanu heparanu, jak i jego znanego receptora BSG24 (Rysunek 3C), natomiast TNFRSF9 wiązał się specyficznie ze swoim znanym receptorem TNFSF9 i nie utracił zdolności wiązania po preinkubacji z rozpuszczalną heparyną. Białko 3 na Rysunku 3B reprezentuje TNFRSF9.
Dla obu linii komórkowych przedstawiono również rozkład gRNA w kontrolnej bibliotece mutantów po 3 dniach (oraz 9, 14 i 16 dniach po transdukcji) (Tabela uzupełniająca 1). Rozkład gRNA wykazał, że złożoność biblioteki została zachowana przez cały czas trwania eksperymentu (Rysunek 5A). Przesiew genetyczny w celu identyfikacji liganda dla TNFSF9 przeprowadzono w 14. dniu po transdukcji, natomiast dla RH5 w 9. dniu po transdukcji. Jakość techniczną przesiewów oceniono poprzez analizę rozkładu obserwowanych zmian krotności (fold-change) gRNA celujących w referencyjny zestaw genów nieessentialnych w porównaniu z rozkładem dla referencyjnego zestawu genów essentialnych22 (Rysunek 5B). Ponadto analiza wzbogacenia na poziomie szlaków wykazała, że oczekiwane szlaki essentialne zostały zidentyfikowane i uległy znacznemu wzbogaceniu w populacji „drop-out” podczas porównywania próbki kontrolnej z oryginalną biblioteką plazmidową. Przykład dla próbki NCI-SNU-1 z 14. dnia przedstawiono na Rysunku 5C.
Do identyfikacji receptora z przesiewów fenotypowych wykorzystano rozkład gRNA w populacji kontrolnej w porównaniu z populacją posegregowaną, stosując funkcję -test programu MAGeCK (zob. tabelę uzupełniającą 1 w celu uzyskania podsumowania genów z MAGeCK).Zmodyfikowany wynik RRA raportowany przez MAGeCK w analizie na poziomie genów został naniesiony na wykres w odniesieniu do genów uszeregowanych według wartości p. Wynik RRA w MAGeCK stanowi miarę, w której gRNA są szeregowane konsekwentnie wyżej, niż oczekiwano. W przesiewie dla TNFRSF9 najlepszym wynikiem był TNFSF9, który jest znanym partnerem wiążącym dla TNFRSF9 (rycina 5D). Ponadto zidentyfikowano szereg genów związanych ze szlakiem TP53. W przypadku RH5, poza znanym receptorem (BSG) i genem niezbędnym do produkcji siarkowanych GAG (SLC35B2), zidentyfikowano również dodatkowy gen (SLC16A1) (rycina 5E). SLC16A1 jest chaperonem niezbędnym do transportu BSG na powierzchnię komórek25. Wspólnie wyniki te demonstrują zdolność przesiewu do identyfikacji bezpośrednio oddziałujących receptorów oraz komponentów komórkowych niezbędnych do ekspresji tego receptora na powierzchni komórek w formie funkcjonalnej.

Rycina 1: Przegląd podejścia do przesiewu genetycznego w celu identyfikacji receptorów powierzchniowych komórek. Test ten składa się z trzech głównych etapów: Po pierwsze, białka rekombinowane reprezentujące ektodomenę receptorów powierzchniowych komórek są ekspresowane w linii komórkowej zdolnej do wprowadzania krytycznych strukturalnie modyfikacji potranslacyjnych, takich jak komórki HEK293. Monomeryczne ektodomeny białek są oligomeryzowane poprzez koniugację ze streptawidyną-PE, aby zwiększyć ich awidność wiązania. Po drugie, te awidne sondy są wykorzystywane w testach wiązania komórkowego, w których silne barwienie linii komórkowych, wskazane przez wyraźne przesunięcie fluorescencji PE (na zielono) w porównaniu z białkiem kontrolnym negatywnym (na czarno), wykazuje obecność partnera wiążącego na powierzchni komórki. Po trzecie, wybierane są linie komórkowe wykazujące obecność receptora i ekspresujące Cas9, a następnie przeprowadzany jest przesiew w skali genomu z użyciem gRNA celujących w zdecydowaną większość genów kodujących białka. Podczas generowania bibliotek mutantów powszechnie stosuje się 30% wydajność transdukcji, co opiera się na prawdopodobieństwie rozkładu Poissona, które zapewnia, że każda komórka otrzyma pojedynczy gRNA, dzięki czemu wynikający fenotyp można przypisać do konkretnego nokautu. Marker BFP ekspresowany przez transdukowane komórki jest wykorzystywany do selekcji komórek zawierających gRNA za pomocą FACS. Przesiewy fenotypowe są przeprowadzane od 9 do 16 dni po transdukcji. W dniu przesiewu całkowita populacja komórek mutantów jest dzielona na dwie części. Jedna połowa jest zachowywana jako populacja kontrolna, a druga jest selekcjonowana pod kątem wiązania białka rekombinowanego. Komórki z biblioteki mutantów, które nie są już w stanie wiązać białka rekombinowanego, są sortowane za pomocą FACS, a wzbogacenie gRNA w populacji sortowanej w stosunku do populacji kontrolnej jest wykorzystywane do identyfikacji genów niezbędnych do wiązania znakowanej awidnej sondy na powierzchni komórki. Wskazano etapy protokołu wymagające znacznej ilości czasu. Rycina została zmodyfikowana na podstawie pracy Sharma i wsp.19. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 2: Ustalanie stosunków biotynilowanego białka do streptawidyny-PE przy użyciu metody opartej na teście ELISA. Przykład strategii koniugacji streptawidyny-PE wykorzystanej do wytworzenia sondy awidnej z biotynilowanego białka monomerycznego. Serię rozcieńczeń biotynilowanych monomerów inkubowano z ustalonym stężeniem streptawidyny. Minimalne rozcieńczenie, przy którym nie można było wykryć nadmiaru biotynilowanych monomerów, określono za pomocą testu ELISA. Test ELISA przeprowadzono z wstępną inkubacją różnych rozcieńczeń białka z 10 ng streptawidyny-PE lub bez niej. W obecności streptawidyny-PE wyznaczono minimalne rozcieńczenie, przy którym nie zidentyfikowano sygnału (określone czarnym kółkiem), a następnie obliczono ilość białka wymaganą do nasycenia w celu przygotowania 10-krotnego roztworu zapasowego zawierającego 4 µg/mL streptawidyny-PE. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 3: Reprezentatywne wiązanie białek z liniami komórkowymi. (A) Wiązanie białek z liniami komórkowymi wykazało wyraźny wzrost fluorescencji związanej z komórkami w porównaniu z próbką kontrolną. Obróbka termiczna (80 °C przez 10 min) białka rekombinowanego zniosła całe wiązanie, sprowadzając je do poziomu kontroli negatywnej, co dowodzi, że zdolność wiązania była zależna od prawidłowego sfałdowania białka. (B) Różne klasy zachowania wiązania białek z powierzchniami komórkowymi; zależność od GAG. Od lewej do prawej białka można podzielić na trzy typy: Białka typu 1 adsorbują wyłącznie do HS. Białka te tracą zdolność wiązania po preinkubacji z heparyną w stężeniach powyżej 0.2 mg/mL. Białka typu 2 wiążą się z HS oraz ze specyficznym receptorem. Białka te wykazują częściową utratę wiązania w eksperymentach z preblokowaniem. Białka typu 3 nie wiążą się z HS. Białka te nie tracą zdolności wiązania w porównaniu z liniami rodzicielskimi. (C) Przykład białka (tj. RH5), które wiąże się z HS i specyficznym receptorem w sposób addytywny. Celowanie w receptor (tj. BSG) lub enzymy niezbędne do syntezy HS (np. SLC35B2, EXTL3) powoduje jedynie częściową redukcję wiązania RH5 z komórkami w stosunku do kontroli. W takich eksperymentach w celu ustalenia zachowania wiązania można wykorzystać transdukowane linie poliklonalne. Rycina została zmodyfikowana na podstawie pracy Sharma et al.19. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 4: Wybór klonalnych linii komórkowych o wysokiej aktywności Cas9. Wydajność edycji genomu w liniach poliklonalnych i klonowanych linii komórkowych NCI-SNU-1 oceniono przy użyciu systemu reporterowego GFP-BFP, w którym linie komórkowe transdukowano wirusami z plazmidem targetującym gRNA kodującym GFP lub bez niego (tj. „pustym”). Przedstawiono schemat. Cytometrię przepływową wykorzystano do badania ekspresji BFP i GFP po transdukcji w porównaniu z nieinfekowaną kontrolą. Ekspresję GFP wykorzystano jako wskaźnik aktywności Cas9, natomiast ekspresję BFP do oznaczenia komórek transdukowanych. Profil dla komórek nieinfekowanych i zainfekowanych pustym plazmidem był podobny dla wszystkich klonów. Reprezentatywne profile przedstawiono na lewym panelu. Wszystkie pięć klonów linii komórkowej NCI-SNU-1 wykazało większą utratę GFP w porównaniu z linią poliklonalną (prawy panel), przy czym klon 4 wykazał najwyższą wydajność z najmniejszą populacją refrakterną. Rycina została zmodyfikowana na podstawie pracy Sharma et al.19. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 5: Reprezentatywne wyniki przesiewów genetycznych w celu identyfikacji partnerów wiążących się z powierzchnią komórki. (A) Wykresy funkcji rozkładu skumulowanego porównujące liczebność gRNA w bibliotece plazmidowej z bibliotekami mutantów komórek HEK-293-E i NCI-SNU-1 w 9., 14. i 16. dniu po transdukcji. Dla dowolnej liczby funkcja gęstości skumulowanej podaje procent punktów danych znajdujących się poniżej tego progu. Niewielkie przesunięcie populacji komórek mutantów w porównaniu z pierwotną populacją plazmidową reprezentuje ubytek w podzbiorze gRNA w porównaniu z biblioteką plazmidową. (B) Rozkład logarytmicznych zmian krotności (log-fold changes) w genach, które zostały wcześniej skategoryzowane jako niezbędne (czerwony) lub nieistotne (czarny) w liniach komórkowych HEK293 i NCI-SNU-1. Rozkład zmian krotności dla genów nieistotnych był skupiony wokół ~0, podczas gdy rozkład dla genów niezbędnych przesunął się w lewo w stronę ujemnych zmian krotności. (C) Istotnie wzbogacone szlaki w genach wyczerpanych w populacji kontrolnej mutantów NCI-SNU-1 14 dni po transdukcji. Zidentyfikowano przewidywane, znane szlaki niezbędne dla komórek. (D) Wynik Robust Rank Algorithm (RRA) dla genów wzbogaconych w posortowanych komórkach, które utraciły zdolność wiązania sondy TNFRSF9. Geny zostały uszeregowane zgodnie z wynikiem RRA. W przesiewie zidentyfikowano znany partner interakcji TNFSF9 oraz geny związane ze szlakiem TP53 (oznaczone na czerwono). (E) Uszeregowane wyniki RRA dla genów zidentyfikowanych w analizie wzbogacenia gRNA wymaganych do wiązania RH5 z komórkami HEK293 (lewy panel). SLC35B2 oraz SLC16A1 zostały zidentyfikowane w ramach progu fałszywych odkryć (FDR) wynoszącego 5%. Dwa dodatkowe geny w szlaku biosyntezy HS (tj. EXTL3 i NDST1) zostały zidentyfikowane w ramach FDR wynoszącego 25%. Schemat przedstawiający ogólny szlak biosyntezy GAG z naniesionymi odpowiednimi genami na odpowiadające im etapy (panel 2). W przesiewie nie zidentyfikowano genów wymaganych do przejścia w biogenezę siarczanu chondroityny (tj. CSGALNACT1/2). Rycina została zmodyfikowana na podstawie pracy Sharma i wsp.19. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
| Nazwa plazmidu | Plazmid nr | Zastosowanie |
| Konstrukt ekspresji białka: CD200RCD4d3+4-bio-linker-his | Addgene: 36153 | Produkcja rekombinowanego białka z sekwencjami CD4d3+4, biotyną i znacznikami 6-his. |
| pMD2.G | Addgene: 12259 | Plazmid wykazujący ekspresję białka otoczki VSV-G; produkcja lentiwirusa |
| psPAX2 | Addgene: 12260 | Plazmid pakujący lentiwirusowy, produkcja lentiwirusa |
| Konstrukt Cas9: pKLV2-EF1a-Cas9Bsd-W | Addgene: 68343 | Produkcja linii konstytutywnie eksponującej Cas9 |
| konstrukt ekspresyjny gRNA: pKLV2-U6gRNA5(BbsI)-PGKpuro2ABFP-W | Addgene: 67974 | Wektor ekspresyjny gRNA CRISPR z ulepszonym rusztowaniem oraz markerami puro/BFP |
| Udoskonalona ludzka biblioteka CRISPR do knockoutu w całym genomie | Addgene: 67989 | Biblioteka gRNA skierowana przeciwko 18 010 ludzkim genom, zaprojektowana do zastosowania w lentiwirusach. |
| Konstrukt GFP-BFP: pKLV2-U6gRNA5(gGFP)-PGKBFP2AGFP-W | Addgene: 67980 | Reporter aktywności Cas9 z BFP i GFP. |
| Konstrukt pusty: pKLV2-U6gRNA5(empty)-PGKBFP2AGFP-W | Addgene: 67979 | Reporter aktywności Cas9 (kontrola) z BFP i GFP. |
Tabela 1: Plazmidy wykorzystane w tym podejściu.
| Nazwa buforu | Składniki |
| HBS (10X) | 1.5 M NaCl oraz 200 mM HEPES w wodzie MiliQ, doprowadzić do pH 7.4 |
| PBS (10X) | 80 g NaCl, 2 g KCl, 14.4 g Na2HPO4 oraz 2.4 g KH2PO4 w wodzie MiliQ, doprowadzić do pH 7.4 |
| Bufor fosforanowy sodu (roztwór stokowy 80mM) | 7.1 g Na2HPO4.2H2O, 5.55 g NaH2PO4, doprowadzić do pH 7.4 |
| Bufor wiążący do oczyszczania białek His | 20 mM bufor fosforanowy sodu, 0.5 M NaCl oraz 20 mM imidazolu, doprowadzić do pH 7.4 |
| Bufor elucyjny do oczyszczania białek His | 20 mM bufor fosforanowy sodu, 0.5 M NaCl oraz 400 mM imidazolu, doprowadzić do pH 7.4 |
| Bufor dietanoloaminowy | 10% dietanoloaminy i 0.5 mM MgCl2 w wodzie MiliQ, doprowadzić do pH 9.2: |
| D10 | DMEM, 1% penicyliny-streptomycyny (100 units/mL) oraz 10% inaktywowanego ciepłem FBS |
Tabela 2: Bufory wymagane do niniejszego badania.
| Komponenty | szalka 10-cm | płytka 6-dołkowa |
| komórki 293FT | 70–80% konfluencji | 70–80% konfluencji |
| Pośrodki kompatybilne z transfekcją (Opti-MEM) (Krok 5.1.2) | 3 mL | 500 µL |
| Pośrodki kompatybilne z transfekcją (Opti-MEM) (Krok 5.1.4) | 5 mL | 2 ml |
| Lentiwirusowy wektor transferowy | 3 µg | 0.5 µg |
| psPax2 (patrz tabela 1) | 7.4 µg | 1.2 µg |
| pMD2.G (patrz tabela 1) | 1.6 µg | 0.25 µg |
| odczynnik PLUS | 12 µL | 2 µL |
| Lipofectamine LTX | 36 µL | 6 µL |
| D10 (Krok 7.1.7) | 5 mL | 1,5 mL |
| D10 (Krok 7.1.8 i 7.1.10) | 8 ml | 2 mL |
Tabela 3: Ilości i objętości odczynników dla mieszaniny do pakowania lentiwirusów.
Tabela 4: Sekwencje starterów do amplifikacji gRNA i NGS. Kliknij tutaj, aby wyświetlić ten plik (kliknij prawym przyciskiem myszy, aby pobrać).
| Odczynnik | Objętość na reakcję | Mieszanina Master mix (x38) |
| Q5 Hot Start High-Fidelity 2x | 25 μL | 950 μL |
| mieszanina starterów (L1/U1) (10 μM każdy) | 1 μL | 38 μL |
| DNA genomiczne (1 mg/ml) | 2 μL | 72 μL |
| H2O | 22 μL | 1100 μL |
| Suma | 50 μL | 1900 μL |
Tabela 5: PCR w celu amplifikacji gRNA z próbek o wysokiej złożoności.
| Numer cyklu | Denaturacja | Przyłączanie | Ekstensja |
| 1 | 98 °C, 30s | ||
| 2-24 | 98 °C, 10s | 61 °C, 15s | 72 °C, 20s |
| 25 | 72 °C, 2 min |
Tabela 6: Warunki PCR dla pierwszej reakcji PCR.
| Odczynnik | Objętość na reakcję |
| KAPA HiFi HotStart ReadyMix | 25 μL |
| Mieszanina starterów (starter PE1.0/starter indeksujący) (5 μM każdy) | 2μL |
| Pierwszy produkt PCR (40 pg/μL) | 5 μL |
| H2O | 18 μL |
| Suma | 50 μL |
Tabela 7: PCR do znakowania indeksowego sgRNA z przesiewów genetycznych.
| Numer cyklu | Denaturacja | Przyłączanie | Elongacja |
| 1 | 98 °C, 30s | ||
| 2-15 | 98 °C, 10s | 66 °C, 15s | 72 °C, 20s |
| 16 | 72 °C, 5 min |
Tabela 8: Warunki PCR dla drugiej reakcji PCR.

Rysunek uzupełniający S1: Przewodnik po wyznaczaniu bramek do sortowania populacji niewiążącej. (A) Idealny kandydat białkowy do screeningu powinien wykazywać wyraźne przesunięcie populacji wiążącej w porównaniu z populacją kontrolną, a wiązanie powinno zostać zachowane w komórkach pozbawionych aparatu do biosyntezy HS. Zamiast testowania na liniach komórkowych z nokautem SLC35B2 można przeprowadzić eksperyment z blokowaniem heparyną. (B) Zebrano komórki pozbawione barwienia powierzchniowego z ektodomeny białka, ale wykazujące fluorescencję BFP po transdukcji lentiwirusem. Przedstawione komórki pochodzą ze screeningu w celu identyfikacji receptora dla GABBR222. Rysunek został zmodyfikowany na podstawie pracy Sharma et al.19. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek uzupełniający S2: Wykres PCA oparty na transkryptomice glikoprotein powierzchni komórki z wykorzystaniem danych RNA-seq z ponad 1 000 linii komórkowych nowotworów. Linie komórkowe z bazy Cell Model Passport27 zostały zgrupowane metodą klastrowania K-średnich na podstawie wartości FPKM dla około 1 500 glikoprotein powierzchni komórki. Oznaczono reprezentatywne linie komórkowe z każdego klastra. Klaster 5 składał się w całości z linii komórkowych pochodzenia hematopoetycznego (patrz również Tabela uzupełniająca 2). Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina uzupełniająca S3: Wyniki istotności dla genów eksportu białek i glikozylacji N-łączeniowej z adnotacji KEGG uzyskane z wyników projektu. Skorygowane wyniki istotności Bayesa dla ~330 linii komórkowych (kolumny, nieoznaczone) przedstawiono dla genów szlaku eksportu białek i glikozylacji N-łączeniowej (oś X). Wyniki wyższe niż 0 oznaczają istotne wyczerpanie w populacji mutantów w porównaniu do oryginalnej biblioteki plazmidowej. Geny można podzielić na trzy odrębne klastry reprezentujące różne poziomy istotności w liniach komórkowych. Klastrowanie to może być wykorzystane do wyznaczenia dnia sortowania. Jeśli przesiew zostanie przeprowadzony w późniejszym punkcie czasowym (dzień 16), istnieje możliwość, że geny znane jako istotne dla komórek (klastry 1 i 3) nie zostaną zidentyfikowane. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.
Tabela uzupełniająca 1: Pliki z surową liczbą odczytów oraz pliki gene_summary wygenerowane przez oprogramowanie MAGeCK powiązane z reprezentatywnymi badaniami przesiewowymi genetycznymi. Kliknij tutaj, aby wyświetlić ten plik (kliknij prawym przyciskiem myszy, aby pobrać).
Tabela uzupełniająca 2: Grupowanie linii komórkowych według ekspresji receptorów powierzchniowych komórek. Kliknij tutaj, aby wyświetlić ten plik (kliknij prawym przyciskiem myszy, aby pobrać).
Opisano strategię badań przesiewowych opartą na CRISPR w celu identyfikacji genów kodujących komponenty komórkowe zaangażowane w rozpoznawanie komórek. Podobne podejście wykorzystujące aktywację CRISPR zapewnia również genetyczną alternatywę do identyfikacji bezpośrednio oddziałujących receptorów rekombinowanych białek bez konieczności generowania dużych bibliotek białek26. Jednak jedną z głównych zalet tego podejścia jest to, że ma ono zastosowanie do interakcji, w których pośredniczą cząsteczki powierzchniowe natywnie wyświetlane na komórce i nie zależy od nadekspresji receptorów, co może wpływać na awidność wiązania receptora. W przeciwieństwie do innych metod, technika ta nie przyjmuje żadnych założeń dotyczących biochemicznej natury lub biologii komórki receptorów i daje możliwość badania interakcji, w których pośredniczą białka, które zwykle są trudne do zbadania przy użyciu podejść biochemicznych, takich jak bardzo duże białka lub te, które wielokrotnie przechodzą przez błonę lub tworzą kompleksy z innymi białkami. oraz cząsteczki inne niż białka, takie jak glikany, glikolipidy i fosfolipidy. Biorąc pod uwagę charakter metody w skali genomu, podejście to ma również tę zaletę, że nie tylko identyfikuje receptor, ale także dodatkowe składniki komórkowe, które są wymagane do zdarzenia wiązania, zapewniając w ten sposób wgląd w biologię komórki receptora.
Jednym z głównych ograniczeń tej metody przy użyciu jej do identyfikacji receptora białka sierocego jest początkowy wymóg, aby najpierw zidentyfikować linię komórkową, która wiąże się z białkiem. Nie zawsze jest to łatwe, a identyfikacja linii komórkowej, która wykazuje fenotyp wiązania, który jest również dopuszczalny dla genetycznych badań przesiewowych, może być krokiem ograniczającym czas wdrożenia tego testu. Niektóre linie komórkowe mają tendencję do wiązania się z większą liczbą białek niż inne. Jest to szczególnie istotne w przypadku białek, które wiążą HS, ponieważ białka te mają tendencję do wiązania się z dowolną linią komórkową, która wykazuje łańcuchy boczne HS, niezależnie od natywnego kontekstu wiązania. Dodatkowo zaobserwowaliśmy, że regulacja w górę syndekanów (tj. proteoglikanów zawierających HS) w liniach komórkowych prowadzi do zwiększonego wiązania białek wiążących HS26. Może to być czynnik, który należy wziąć pod uwagę przy wyborze linii komórkowej do badań przesiewowych. Jednak ważne jest również, aby pamiętać, że addytywne wiązanie HS nie zakłóca wiązania z określonym receptorem. Oznacza to, że jeśli zaobserwuje się wiązanie, możliwe jest, że jest ono pośredniczone wyłącznie przez HS, ponieważ wiązanie, w którym pośredniczy HS w tym teście, jest addytywne, a nie współzależne19. W takim scenariuszu opisane podejście do blokowania heparyny może zidentyfikować takie zachowania bez konieczności przeprowadzania pełnego badania genetycznego.
Przydatnym źródłem do wyboru linii komórkowych jest Paszport Modelu Komórki, który zawiera informacje o genomice, transkryptomice i warunkach hodowli dla ~1,000 linii komórek rakowych27. W zależności od kontekstu biologicznego komórki można wybierać na podstawie ich profili ekspresji. Aby pomóc w selekcji linii komórkowych, zgrupowaliśmy ~ 1,000 linii komórkowych w paszporcie modelu komórki zgodnie z ekspresją ~ 1,500 wstępnie opisanych glikoprotein powierzchniowych komórek ludzkich28 (rysunek uzupełniający 2; informacje o klastrach dla każdej linii komórkowej wraz z warunkami wzrostu znajdują się w Tabeli Uzupełniającej 2). Podczas testowania wiązania białka o nieznanej funkcji przydatne jest wybranie panelu reprezentatywnych linii komórkowych z każdego klastra, aby zwiększyć szansę na pokrycie szerokiego zakresu receptorów. Mając wybór, zaleca się wybór linii komórkowych, które są łatwe do hodowli i łatwe do transdukcji. Ponieważ te linie komórkowe będą wykorzystywane w badaniach przesiewowych w skali genomu, preferowane jest, aby można je było łatwo hodować w dużych ilościach i były podatne na transdukcję lentiwirusową, ponieważ jest to najpowszechniej dostępna metoda dostarczania sgRNA do genetycznych badań przesiewowych opartych na CRISPR na późniejszych etapach.
Ogólnie rzecz biorąc, selekcja fenotypów jest przeprowadzana w jednym sortowaniu. Decyduje o tym jednak jasność wybarwionej populacji komórek w porównaniu z kontrolą. Iteracyjne rundy selekcji mogą być stosowane w scenariuszach, w których stosunek sygnału do szumu pożądanego fenotypu jest niski lub gdy celem badania przesiewowego jest identyfikacja mutantów o silnych fenotypach. W przypadku stosowania podejścia do selekcji iteracyjnej dla badań przesiewowych opartych na FACS ważne jest, aby wziąć pod uwagę, że proces sortowania może spowodować śmierć komórki, głównie ze względu na samą siłę sortera. W związku z tym nie wszystkie zebrane komórki będą reprezentowane w następnej rundzie sortowania.
Złożoność biblioteki jest bardzo ważnym czynnikiem w przeprowadzaniu udanych badań genetycznych, zwłaszcza w przypadku badań przesiewowych selekcji negatywnej, ponieważ zakres ich uszczuplenia można określić tylko poprzez porównanie wyników z tym, co było obecne w bibliotece wyjściowej. W przypadku ekranów selekcji negatywnej często utrzymuje się biblioteki o złożoności 500-1,000x. Ekrany selekcji pozytywnej są jednak bardziej odporne na rozmiary bibliotek, ponieważ w takich ekranach oczekuje się, że tylko niewielka liczba mutantów zostanie wybrana dla określonego fenotypu. W związku z tym, na opisanym tutaj ekranie wyboru pozytywnego, rozmiar biblioteki może zostać zmniejszony do złożoności 50-100x bez uszczerbku dla jakości ekranu. Na dodatek, w badaniach tych możliwe jest również użycie biblioteki kontrolnej dla danej linii komórkowej w danym dniu jako "kontroli ogólnej" dla wszystkich próbek posortowanych w danym dniu dla tej danej linii komórkowej. Zmniejszy to liczbę bibliotek kontrolnych, które muszą być tworzone i sekwencjonowane.
Inną ważną kwestią przy stosowaniu tego podejścia są ograniczenia badań przesiewowych utraty funkcji w identyfikacji genów, które są niezbędne do wzrostu komórek in vitro. Czas badań przesiewowych ma w tym względzie kluczowe znaczenie, ponieważ im dłużej zmutowane komórki są trzymane w hodowli, tym większe prawdopodobieństwo, że komórki z mutacjami w podstawowych genach staną się niezdolne do życia i nie będą już reprezentowane w bibliotece mutantów. Niedawne badania genetyczne przeprowadzone w ramach inicjatywy Project Score na ponad 300 liniach komórkowych pokazują, że wiele genów w szlaku wydzielania białek i N-glikozylacji z adnotacją KEGG jest często identyfikowanych jako niezbędne dla wielu linii komórkowych (rysunek uzupełniający 3)29. Można to wziąć pod uwagę przy projektowaniu badań przesiewowych, jeśli wpływ genów niezbędnych do proliferacji i żywotności ma być badany w kontekście procesu rozpoznawania komórkowego. W takim przypadku przeprowadzenie badań przesiewowych we wczesnym punkcie czasowym (np. 9 dnia po transdukcji) byłoby ogólnie właściwe. Jeśli jednak podejście to jest stosowane do identyfikacji kilku celów o silnym wpływie na wielkość, a nie ogólnych szlaków komórkowych, właściwe może być przeprowadzenie badań przesiewowych w późniejszym czasie (np. 15-16 dnia po transdukcji).
Wyniki badań przesiewowych są bardzo solidne; W ośmiu badaniach przesiewowych wiązania rekombinacji białek przeprowadzonych w przeszłości, receptor powierzchniowy komórki był największym hitem w każdym przypadku19. Stosując to podejście do identyfikacji partnera interakcji, należy się zatem spodziewać, że receptor i czynniki przyczyniające się do jego prezentacji na powierzchni będą identyfikowane z dużą pewnością statystyczną. Po przeprowadzeniu badania przesiewowego i walidacji trafienia za pomocą pojedynczego nokautu gRNA, można przeprowadzić dalsze obserwacje przy użyciu istniejących metod biochemicznych, takich jak AVEXIS4 i bezpośrednie saturywalne wiązanie oczyszczonych białek za pomocą powierzchniowego rezonansu plazmonowego. Opisane tutaj podejście ma zastosowanie do wszystkich białek, dla których możliwe jest wygenerowanie rozpuszczalnej rekombinowanej sondy wiążącej.
Podsumowując, jest to podejście oparte na nokaucie CRISPR w skali genomu w celu identyfikacji interakcji, w których pośredniczą białka powierzchniowe komórki. Metoda ta ma ogólne zastosowanie do identyfikacji szlaków komórkowych wymaganych do rozpoznawania powierzchni komórki w szerokim zakresie różnych kontekstów biologicznych, w tym między własnymi komórkami organizmu (np. rozpoznawanie neuronalne i immunologiczne), a także między komórkami gospodarza a białkami patogenu. Metoda ta stanowi genetyczną alternatywę dla podejść biochemicznych przeznaczonych do identyfikacji receptorów, a ponieważ nie wymaga żadnych wcześniejszych założeń dotyczących biochemicznej natury lub biologii komórki receptorów, ma ogromny potencjał do dokonywania zupełnie nieoczekiwanych odkryć.
Autorzy nie mają nic do ujawnienia.
Ta praca została wsparta przez grant Wellcome Trust numer 206194 przyznany GJW. Dziękujemy placówce Cytometry Core: Bee Ling Ng, Jennifer Graham, Sam Thompson i Christopher Hall za pomoc w FACS.
| Nazwa | Firma | Numer katalogowy | Komentarze |
|---|---|---|---|
| Przeciw-mysia fosfataza alkaliczna | Sigma | A4656 | |
| Blasticidin | Chem-Cruz | SC-204655 | |
| Krew i Cell Culture DNA Maxi Kit | Qiagen | 13362 | |
| BSA | Sigma | A9647-100G | |
| Dietanoloamina | Sigma | 398179 | |
| DMEM | Gibco | 31966-021 | |
| Zestaw do krwi i tkanek Dneasy | Qiagen | 69504 | |
| DynaMag-96 Magnes boczny | Invitrogen | 12331D | |
| HEK293T opakowanie komórki | ATCC | CRL-3216 | |
| Heparyna | Sigma | H4784-1G | |
| KAPA HiFi HotStart ReadyMix | Kapa | KK2602 | |
| Lipofectamina LTX z odczynnikiem PLUS | Invitrogen | 15338100 | |
| Sortownik komórek MoFlo XDP | BD | ||
| Ni2+-NTA koraliki agarozowe | Jena Bioscience | AC-501-25 | |
| OPTI-MEM | Life Technologies | 31985-070 | |
| OX-68 przeciwciało | AbD Serotec | MCA1022R | |
| fosforan p-nitrofenylu | Sigma | 1040-506 | |
| PD-10 kolumny odsalające | GE healthcare | 17085101 | |
| Polybrene | Millipore | TR-1003-G | |
| Probówki polipropylenowe o objętości łoża 5 ml | Qiagen | 34964 | |
| Proteinaza K, rekombinowana, klasa PCR | Roche | 3115879001 | |
| Puromycin | Gibco | A11138-03 | |
| Q5 Hot Start High-Fidelity 2 razy; | |||
| Zestaw do oczyszczania QIAquick PCR | Qiagen | 28104 | |
| Filtr SCFA | Nalgene | 190-2545 | |
| Sony Sortownik komórek | Koraliki Sony | ||
| SPRI (koraliki Agencourt AMPure XP) | Beckman | A63881 | |
| Płytki do mikromiareczkowania pokryte streptawidyną | Nalgene | 734-1284 | |
| Streptawidyna-PE | Biolegend | 405204 |