Aby zademonstrować skuteczność metody χCRAC, przeprowadzono eksperyment z analizą przebiegu czasowego z użyciem szczepów drożdży eksponujących białko Nrd1 z etykietą HTP. Szczegółowy schemat przedstawiający zasadę działania tej metody znajduje się na Rysunku 1. Podobnie jak Nab3, białko Nrd1 bierze udział w jądrowej degradacji różnorodnych transkryptów RNA37. Wcześniejsze prace zespołu Corden sugerowały, że wiązanie Nrd1 z jego celami RNA zmienia się znacząco, gdy komórki są poddane głodzeniu glukozowemu28,38. W związku z tym komórki rosnące wykładniczo w medium zawierającym glukozę (SD-TRP) przeniesiono do tego samego medium bez glukozy (S-TRP) w określonych odstępach czasu, aby monitorować dynamiczne zmiany w oddziaływaniach Nrd1-RNA. Próbki pobrano i poddano sieciowaniu w komorze Vari-X-linker (Rysunek 3A) przed przeniesieniem, a następnie po 1, 2, 4, 8, 14 i 20 min. Medium stosowane do wzrostu komórek było celowo pozbawione tryptofanu, aby ograniczyć absorpcję promieniowania UV przez ten aminokwas aromatyczny. Należy pamiętać, że najlepiej stosować syntetyczne medium sterylizowane przez filtrację, ponieważ autoklawowanie medium może prowadzić do karmelizacji cukrów, co następnie obniża wydajność sieciowania.
Rysunek 4A przedstawia reprezentatywny autoradiogram z eksperymentu χCRAC. Należy zauważyć, że w tym przykładzie próbki nie zostały połączone. Zamiast tego każda z nich została poddana osobnej elektroforezie w żelu. Jest to zalecane podczas wstępnych testów eksperymentalnych, aby wykazać, że białko skutecznie sieciuje się z RNA we wszystkich badanych punktach czasowych. Szczególnie intensywny sygnał zaobserwowano przy spodziewanej masie cząsteczkowej RBP, co odpowiada białku związanemu z bardzo krótkimi, znakowanymi promieniotwórczo cząsteczkami RNA, które nie nadają się do sekwencjonowania. Dlatego wyizolowano rozmyty sygnał powyżej tego prążka, który reprezentuje białko sieciowane z dłuższymi fragmentami RNA. Fragment wycięto z obszaru tuż powyżej prążka białka, obejmując dodatkowo około 30 kDa. Rysunek 4B przedstawia autoradiogram po wycięciu, na którym białko sieciowane z krótkimi RNA pozostało w żelu, a wcześniej rozmyty sygnał został usunięty.
Po odwrotnej transkrypcji bibliotekę cDNA należy powielić za pomocą PCR. Należy jednak unikać nadmiernej amplifikacji biblioteki, ponieważ może to wprowadzić błąd w postaci preferencyjnego powielania niektórych sekwencji przez polimerazę oraz generować artefakty PCR. Nadmiernie powielone biblioteki zawierają również dużą liczbę sekwencji duplikowanych, co prowadzi do marnowania odczytów w sekwenatorze. Aby obliczyć optymalną liczbę cykli PCR dla amplifikacji końcowej biblioteki, alikwot cDNA powielono metodą qPCR z użyciem oligonukleotydów P5 i BC. Jako liczbę cykli PCR wybrano pierwszy cykl, w którym biblioteka osiągnęła szczyt fluorescencji. Rycina 4C przedstawia przykład qPCR z typowej biblioteki cDNA, w której szczytowa liczba cykli wyniosła 16. Wartość ta została następnie zastosowana w końcowej reakcji PCR dla χCRAC. W celu przetworzenia danych z sekwencjonowania wykorzystano oprogramowanie opracowane wcześniej w naszym laboratorium (pyCRAC) oraz odpowiadający mu rurociąg analityczny do analizy kinetycznych danych CRAC (Nues et al., 2017; https://git.ecdf.ed.ac.uk/sgrannem/pycrac, https://bitbucket.org/sgrann/kinetic_crac_pipeline/src/default/). Te narzędzia programistyczne typu open source umożliwiają demultipleksowanie i przycinanie danych, usuwanie duplikatów PCR, identyfikację statystycznie istotnych pików, grupowanie odczytów w sekwencje ciągłe oraz identyfikację motywów wiązania39. Dalsze szczegóły dotyczące działania tych narzędzi znajdują się na ich odpowiednich stronach internetowych.
Rozpoczęliśmy również opracowywanie protokołu χCRAC dla komórek ssaków. Większość linii komórkowych ssaków hodowana jest w monowarstwie, a taca w naszym urządzeniu do sieciowania z workiem przepuszczającym promienie UV nie nadaje się do eksperymentów z komórkami adherentnymi. Aby rozwiązać ten problem, opracowaliśmy platformę, która umożliwia naświetlanie UV 1–2 szalek Petriego (o średnicy 150 mm i głębokości 25 mm) z komórkami adherentnymi (Rycina 3B). W ramach pierwszego testu wydajność sieciowania dla komórek ssaków zmierzono poprzez sieciowanie i wychwytywanie stabilnie znakowanego białka GFP-RBM7 przy użyciu przeciwciał anty-GFP i tradycyjnej puryfikacji opartej na metodzie CLIP. Jak pokazano na Rycynie 5A, urządzenie do sieciowania pozwoliło na odzyskanie kompleksów białko-RNA z komórek ssaków hodowanych w monowarstwie przy użyciu promieniowania UV o długości 254 nm z wydajnością porównywalną do powszechnie stosowanego urządzenia do naświetlania UV. Jednak standardowe plastikowe naczynia do hodowli komórkowej, zwykle używane w eksperymentach z sieciowaniem UV, są nieprzepuszczalne dla promieniowania UV o długości 254 nm. W związku z tym, w naszym urządzeniu komórki otrzymywałyby naświetlanie tylko z górnego rzędu lamp UV. Aby temu zaradzić, opracowaliśmy przepuszczalną dla UV kwarcową szalkę Petriego do hodowli komórek i sieciowania. Zastosowanie kwarcowych naczyń hodowlanych wykazało skuteczne odzyskiwanie kompleksów białko-RNA już przy 2 s naświetlania UV (Rycina 5B). W połączeniu z metodami wychwytywania RBP dla komórek ssaków, takimi jak technologie CLIP, te krótkie czasy sieciowania są odpowiednie do badań kinetycznych w celu odzyskania spatiotemporalnych profili wiązania RNA przez RBP w odpowiedzi na stresy genotoksyczne lub szybką deplecję czynników białkowych, a także równolegle z synchronizacją transkrypcyjną lub cyklu komórkowego.
Rycina 6 przedstawia kilka przykładów danych Nrd1 przetworzonych za pomocą potoku χCRAC. Rycina ta została przygotowana przy użyciu plików bedgraph wygenerowanych przez potok oraz pakietu python GenomeBrowser (https://pypi.org/project/GenomeBrowser/1.6.3/), który zaprojektowaliśmy w celu uproszczenia tworzenia obrazów z przeglądarki genomu w jakości publikacyjnej. Szare prostokąty reprezentują regiony genomowe, w których ekspresowane są niekodujące RNA, takie jak kryptyczne niestabilne transkrypty (CUTs), stabilne niecharakteryzowane transkrypty (SUTs)40 oraz niestabilne transkrypty wrażliwe na Xrn1 (XUTs)41. Dane na Rycini 6 pokazują, że Nrd1 wiąże się z wieloma z tych niekodujących transkryptów RNA, co jest zgodne z tezą, że białko to uczestniczy w degradacji tej klasy transkryptów42. Rycina 6A przedstawia region o długości ~15 kb na chromosomie IV. Odnotowano tutaj znaczący wzrost wiązania Nrd1 do transkryptów kodujących wysokowydajne transportery glukozy HXT6 i HXT7, z których oba ulegają nadekspresji podczas głodzenia glukozowego. Prawdopodobnie terminacja transkrypcji przez kompleks NNS może wpływać na kinetykę indukcji tych genów podczas głodzenia glukozowego. Rycina 6B pokazuje przykład sieciowania Nrd1 z transkryptem Imd3, o którym wiadomo, że jest regulowany przez Nab343. W tym przypadku dane wykazały znaczną redukcję wiązania podczas głodzenia glukozowego. Wcześniejsze prace wykazały zmniejszone wiązanie Nab3 do transkryptu Tye7 podczas głodzenia glukozowego44. Zgodnie z tą obserwacją dane χCRAC sugerują, że wiązanie Nrd1 zmniejszyło się podczas głodzenia glukozowego, a sieciowanie Nrd1 z Tye7 było najniższe po 8 min stresu (Rycina 4C). Wydaje się jednak, że efekt ten był jedynie przejściowy, ponieważ po 14 min głodzenia glukozowego wiązanie Nrd1 powróciło do poziomów wyjściowych.

Rycina 1: Schematyczny przebieg protokołu χCRAC. Oznakowane szczepy hodowano do osiągnięcia pożądanej gęstości. RBP oznacza białko wiążące RNA. Następnie pobrano próbkę referencyjną i poddano ją sieciowaniu światłem UV o długości 254 nm. Pozostałe komórki zebrano poprzez filtrację, a następnie szybko przeniesiono do medium indukującego stres. W opisanym tutaj eksperymencie χCRAC próbki pobierano i sieciowano 1, 2, 4, 8, 14 i 20 min po przeniesieniu (1). Następnie oczyszczono interesujące RBP za pomocą wysokostringentnego dwuetapowego oczyszczania powinowactwa (2). Dalej, wychwycone sieciowane RNA poddano częściowemu trawieniu RNazami, znakowaniu radioaktywnym na końcu 5’ i ligacji adapterów (3). Adaptery 5’ zawierały unikalne sekwencje kodów kreskowych „in-read”, aby poszczególne próbki można było oddzielić bioinformatycznie po sekwencjonowaniu. Kompleksy RBP-RNA zostały następnie eluowane, połączone i wspólnie wytrącone (4), rozdzielone za pomocą SDS-PAGE i zwizualizowane poprzez autoradiografię (5). Następnie z żelu wycięto pojedynczy pasek zawierający sygnał radioaktywny tuż nad głównym prążkiem, zaznaczony czerwoną przerywaną ramką na obrazie autoradiograficznym (5). Paski żelu potraktowano proteazą K, a następnie wyekstrahowano RNA (6), przekształcono w cDNA i namnożono za pomocą PCR (7). Etap PCR wprowadził dodatkowe kody kreskowe (żółty blok wprowadzony przez oligo P7), dzięki czemu wiele bibliotek można było zmultipleksować w jednej ścieżce. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rycina 2: Sieciowanie i filtracja próżniowa. (A) Układ sieciujący. Zawiesinę komórkową wlewa się do lejka znajdującego się w prawym górnym rogu urządzenia (zob. również Rycina 3A dla zbliżenia), a następnie utrzymuje się ją w przezroczystym dla UV worku umieszczonym w środkowej tacy. Worek jest flankowany przez dwie przesłony, które pozostają zamknięte do momentu wydania przez użytkownika polecenia rozpoczęcia etapu naświetlania. Komórki są naświetlane światłem UV z tac znajdujących się zarówno powyżej, jak i poniżej. Urządzenie jest wyposażone w lampy UV o długości fali 254 nm i 365 nm, przy czym ta druga jest stosowana w eksperymentach PAR-CLIP. Obsługa urządzenia odbywa się za pomocą panelu dotykowego w prawym górnym rogu, który umożliwia kontrolę dawki UV lub czasu ekspozycji. (B) Po sieciowaniu komórki są odprowadzane z lewej strony urządzenia. Zawiesiny komórkowe są odzyskiwane za pomocą próżni i odprowadzane do kolby szklanej, z której mogą zostać następnie przelane do urządzenia do filtracji próżniowej w celu zbioru. (C) Urządzenia do filtracji próżniowej. Są one otwierane i zamykane za pomocą klipsa, pomiędzy który wkłada się filtr. W przypadku bardzo krótkich serii czasowych równolegle użyto czterech urządzeń filtrujących, aby nie tracić czasu na wymianę filtrów. (D) Po filtracji nadsącz pożywki był odprowadzany do kolb w celu późniejszej utylizacji. Poniżej urządzeń do filtracji próżniowej zainstalowano zawory, aby utrzymać próżnię w układzie po wyjęciu filtra. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rycina 3: Sieciowanie komórek w zawiesinie w porównaniu do komórek adherentnych. (A) Urządzenie do sieciowania z komorą Vari-X-linker dla komórek w zawiesinie. Hodowlę komórkową wlewa się do wlotu próbki (lejka) znajdującego się w prawym górnym rogu tacki. (B) Tacka, w której można umieścić plastikowe lub kwarcowe szalki Petriego do sieciowania komórek adherentnych lub małych objętości komórek w zawiesinie. Kliknij tutaj, aby zobaczyć powiększoną wersję tej ryciny.

Rysunek 4: Przygotowanie biblioteki. (A) Przykład autoradiogramu z eksperymentu Nrd1-HTP χCRAC. Silny, skoncentrowany sygnał reprezentuje białko sieciowane z bardzo krótkimi RNA, natomiast smuga powyżej reprezentuje białko sieciowane z RNA o długości wystarczającej do sekwencjonowania. (B) Smugę wycięto, co pokazano na autoradiogramie wykonanym po wycięciu fragmentu żelu. (C) Reprezentatywny wynik qPCR z biblioteki cDNA χCRAC. W tym przykładzie maksymalną amplifikację cDNA osiągnięto w 16 cyklach. Zatem do końcowej amplifikacji wykorzystano 16 cykli. Słupek błędu reprezentuje odchylenie standardowe trzech technicznych powtórzeń qPCR. (D) Przykład obrazu fosforowego biblioteki cDNA na żelu TBE 6%. (E) Analiza długości i jakości cDNA za pomocą kapilarnej elektroforezy na chipie. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 5: Eksperyment iCLIP z wysoką aktywnością RNazy w celu przetestowania sieciowania w komórkach ssaków. Przedstawiono autoradiogramy z eksperymentów iCLIP GFP-RBM7, w których badano wydajność odzyskiwania RNP przy różnych energiach sieciowania. Immunoprecypitacje przeprowadzono z użyciem przeciwciał anty-GFP sprzężonych z kulkami magnetycznymi na sieciowanych komórkach stabilnie eksprimujących GFP-RBM7. Immunoprecypitaty inkubowano z wysokimi stężeniami RNazy I w celu przycięcia związanych RNA do krótkich, jednolitych długości. RNP wizualizowano poprzez znakowanie 32P i SDS-PAGE; migrują one jako wyraźny prążek, blisko poziomu migracji białka niesieciowanego. Kwantyfikacja przedstawia wyniki analiz densytometrycznych sygnału znakowanego promieniotwórczo RBM7-RNA znormalizowanego do sygnału western blot anty-GFP. (A) Kinetyka czasu sieciowania powszechnie stosowanego sieciowacza UVP w porównaniu z naszym sieciowaczem (Vari-X-linker; VxL). (B) Kinetyka czasu sieciowania naszego sieciowacza na naczyniach hodowlanych kwarcowych (lewo) i plastikowych (prawo). Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina 6: Przykładowe wykresy z przeglądarki genomu pokazujące zdolność metody χCRAC do wykazania różnicowego i czasowego wiązania Nrd1 z jego celami. Każda ramka przedstawia wykresy dla poszczególnych regionów genomowych. Strzałki wskazują nić, na której zakodowane są geny (strzałka skierowana w lewo = nić minus; strzałka skierowana w prawo = nić plus). Punkty czasowe (min) są oznaczone jako t0, t1, t2 itd. na osiach y każdego wykresu pomocniczego. Podano rzymskie cyfry oznaczające chromosomy oraz współrzędne. (A) Po pozbawieniu glukozy Nrd1 wiąże dwa wysokopowinowate transportery glukozy, HXT6 i HXT7, które oba wykazują zwiększoną ekspresję w tych warunkach. (B) Zaobserwowano wiązanie Nrd1 z Imd3, który jest już zwalidowanym celem Nab344, przy czym intensywność wiązania maleje po głodzeniu glukozowym. (C) Wiązanie Tye7 przez Nrd1 ma charakter dynamiczny i przejściowy, zmniejszając się po głodzeniu glukozowym do minimum po 8 min stresu. Następnie jednak wiązanie powraca do poziomów podstawowych po 14 min. Odczyty zostały znormalizowane do „liczby odczytów na milion” (RPM; oś y). Szare pola oznaczają regiony kodujące niekodujące RNA. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.