$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Aby zademonstrować skuteczność metody χCRAC, przeprowadzono eksperyment z udziałem szczepów drożdży wykazujących ekspresję białka Nrd1 znakowanego przez HTP. Szczegółowy schemat opisujący, jak działa ta metoda, znajduje się w Rysunek 1. Podobnie jak Nab3, Nrd1 bierze udział w rozpadzie jądrowego RNA różnych transkryptów RNA37. Wcześniejsze prace z laboratorium Corden sugerowały, że wiązanie Nrd1 z jego celami RNA zmienia się znacząco, gdy komórki są poddawane głodowi glukozy28,38. W związku z tym komórki rosnące wykładniczo w pożywce zawierającej glukozę (SD-TRP) zostały przesunięte na to samo podłoże bez glukozy (S-TRP) w ciągu czasu, aby monitorować dynamiczne zmiany w interakcjach Nrd1-RNA. Próbki pobrano i usieciowano w komorze Vari-X-linker (Rysunek 3A) przed zmianą, a następnie po 1, 2, 4, 8, 14 i 20 minutach. Pożywka używana do wzrostu komórek była celowo uboga w tryptofan, aby zmniejszyć absorpcję promieniowania UV przez ten aromatyczny aminokwas. Należy pamiętać, że najlepiej jest używać pożywki syntetycznej, która jest sterylizowana filtracyjnie, ponieważ autoklawowanie pożywki może prowadzić do karmelizacji cukrów. To z kolei zmniejsza 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 razem. Zamiast tego każdy z nich był uruchamiany indywidualnie na żelu. Jest to zalecane do wstępnych testów eksperymentalnych, aby wykazać, że białko skutecznie łączy się krzyżowo z RNA we wszystkich badanych punktach czasowych. Szczególnie intensywny sygnał zaobserwowano przy oczekiwanej masie cząsteczkowej RBP, reprezentującej białko związane z bardzo krótkimi, znakowanymi radioaktywnie RNA, które nie nadają się do sekwencjonowania. W związku z tym wyizolowano sygnał rozmazowy powyżej tego pasma, który jest białkiem usieciowanym z dłuższymi fragmentami RNA. Fragment został wycięty tuż nad pasmem białka plus około 30 kDa. Rysunek 4B pokazuje autoradiogram po wycięciu, z białkiem usieciowanym do krótkich RNA pozostawionych w żelu i wcześniej rozmazanym sygnałem teraz wyciętym.
Po odwrotnej transkrypcji, biblioteka cDNA musi zostać amplifikowana za pomocą PCR. Należy jednak unikać nadmiernej amplifikacji biblioteki, ponieważ może to wprowadzić stronniczość w kierunku sekwencji preferencyjnie amplifikowanych przez polimerazę i generować artefakty PCR. Nadmiernie wzmocnione biblioteki zawierają również dużą liczbę zduplikowanych sekwencji, które są marnotrawione podczas odczytywania sekwencera. W celu obliczenia idealnej liczby cykli PCR do amplifikacji biblioteki końcowej, podwielokrotność cDNA amplifikowano za pomocą qPCR przy użyciu oligonukleotydów P5 i BC. Pierwszy cykl, w którym biblioteka osiągnęła szczyt fluorescencji, został wybrany jako liczba cykli PCR. Rysunek 4C przedstawia przykład qPCR z typowej biblioteki cDNA, który dał szczytową liczbę cykli wynoszącą 16. Wartość ta została następnie wykorzystana do końcowego pcr χCRAC. W celu przetworzenia zsekwencjonowanych danych wykorzystaliśmy oprogramowanie opracowane wcześniej w naszym laboratorium (pyCRAC) i odpowiedni potok do analizy danych kinetycznych CRAC (Nues i in., 2017; https://git.ecdf.ed.ac.uk/sgrannem/pycrac, https://bitbucket.org/sgrann/kinetic_crac_pipeline/src/default/). Te narzędzia programowe typu open source umożliwiają demultipleksowanie i przycinanie danych, usuwanie duplikatów PCR, identyfikację statystycznie istotnych pików, odczyty klastrów w ciągłe sekwencje i identyfikowanie motywów wiązania39. Więcej informacji na temat działania tych narzędzi można znaleźć na ich stronach internetowych.
Zaczęliśmy również opracowywać protokół χCRAC dla komórek ssaków. Większość linii komórkowych ssaków jest hodowana jako monowarstwa, a taca w naszym sieciowniku z torbą przepuszczającą promieniowanie UV nie nadaje się do eksperymentów z przylegającymi komórkami. Aby przezwyciężyć ten problem, opracowaliśmy etap, w którym użytkownicy mogą naświetlać promieniami UV 1–2 szalki Petriego (o średnicy 150 mm i głębokości 25 mm) za pomocą przylegających komórek (Rysunek 3B). W pierwszym teście zmierzono skuteczność środka sieciującego dla komórek ssaków poprzez sieciowanie i wychwytywanie stabilnie znakowanego GFP-RBM7 przy użyciu przeciwciał anty-GFP i tradycyjnego oczyszczania opartego na CLIP. Jak pokazano w Rysunek 5A, sieciownik był w stanie odzyskać kompleksy białko-RNA z komórek ssaków hodowanych jako monowarstwa przy użyciu promieniowania UV o długości fali 254 nm z wydajnością porównywalną do powszechnie stosowanego urządzenia do naświetlania UV. Jednak standardowe tworzywa sztuczne do hodowli komórkowych zwykle używane do eksperymentów z sieciowaniem UV są nieprzepuszczalne dla promieniowania UV 254 nm. Dlatego w naszym sieciowniku ogniwa otrzymywałyby promieniowanie tylko z górnego banku lamp UV. Aby temu zaradzić, opracowaliśmy przepuszczalną dla promieniowania UV kwarcową szalkę Petriego do wzrostu komórek i sieciowania. Użycie kultur kwarcowych wykazało solidny odzysk kompleksów białko-RNA przy zaledwie 2 s naświetlania promieniowaniem UV (Rysunek 5B). W połączeniu z metodami wychwytywania RBP dla komórek ssaków, takimi jak technologie CLIP, te krótkie czasy sieciowania są podatne na przebieg czasowy w celu odzyskania czasoprzestrzennych profili wiążących RNA RBP w odpowiedzi na stres genotoksyczny lub szybkie wyczerpanie czynników białkowych, lub równolegle z synchronizacją transkrypcji lub cyklu komórkowego.
Rysunek 6 pokazuje kilka przykładów danych Nrd1 przetwarzanych przez potok χCRAC. Rysunek ten został przygotowany przy użyciu plików bedgraph wygenerowanych przez potok i pakietu GenomeBrowser Pythona (https://pypi.org/project/GenomeBrowser/1.6.3/), który zaprojektowaliśmy, aby uprościć tworzenie obrazów danych w przeglądarce genomu o jakości publikacji. Szare prostokąty reprezentują regiony genomowe, które wyrażają niekodujące RNA, takie jak tajemniczy niestabilny transkrypt (CUT), stabilne niescharakteryzowane transkrypty (SUT)40 i niestabilne transkrypty wrażliwe na Xrn1 (XUTs)41. Dane w Rysunek 6 pokazują, że Nrd1 wiąże się z wieloma z tych niekodujących transkryptów RNA, co jest zgodne z ideą, że białko to jest zaangażowane w degradację tej klasy transkryptów42. Rysunek 6A pokazuje region ~15 kb na chromosomie IV. W tym przypadku nastąpił znaczny wzrost wiązania Nrd1 z transkryptami kodującymi transportery glukozy o wysokim powinowactwie HXT6 i HXT7, z których oba są regulowane w górę podczas głodu glukozy. Jest prawdopodobne, że przerwanie transkrypcji przez kompleks NNS może wpływać na kinetykę indukcji tych genów podczas głodu glukozy. Rysunek 6B pokazuje przykład sieciowania Nrd1 do transkryptu Imd3, o którym wiadomo, że jest regulowany przez Nab343. W tym przypadku dane wykazały znaczące zmniejszenie wiązania przy głodzie glukozy. Poprzednie prace wykazały zmniejszone wiązanie Nab3 z transkryptem Tye7 podczas głodu glukozy44. Zgodnie z tą obserwacją, dane χCRAC sugerują, że wiązanie Nrd1 zmniejszyło się podczas głodu glukozy, a sieciowanie Nrd1 z Tye7 było najniższe po 8 minutach stresu (Figura 4C). Wydaje się jednak, że efekt ten był tylko przejściowy, ponieważ po 14 minutach głodu glukozy wiązanie Nrd1 wróciło do poziomów wyjściowych.

Rysunek 1: Schematyczne przedstawienie protokołu χCRAC. Oznaczone szczepy były uprawiane aż do pożądanej gęstości. RBP wskazuje na białko wiążące RNA. Następnie pobrano próbkę referencyjną i usieciowano ją światłem UV o długości fali 254 nm. Pozostałe komórki zebrano przez filtrację, a następnie szybko przeniesiono do pożywki wywołującej stres. W opisanym tutaj eksperymencie χCRAC próbki pobrano i usieciowano 1, 2, 4, 8, 14 i 20 minut po zmianie (1). Interesujący RBP został następnie oczyszczony przy użyciu bardzo rygorystycznego dwuetapowego oczyszczania powinowactwa (2). Następnie przechwycone usieciowane RNA zostały częściowo strawione RNazami, znakowane radioaktywnie na końcu 5', a adaptery zostały na nich podwiązane (3). Adaptery 5' zawierały unikalne sekwencje kodów kreskowych "w odczycie", dzięki czemu poszczególne próbki można było bioinformatyczną separacją po sekwencjonowaniu. Kompleksy RBP-RNA zostały następnie eluowane, połączone i wytrącone razem (4), rozwiązane przez SDS-PAGE i uwidocznione za pomocą autoradiografii (5). Następnie, pojedynczy plasterek żelu zawierający sygnał radioaktywny tuż nad głównym pasmem, zilustrowany przerywaną czerwoną ramką na obrazie autoradiograficznym, został wycięty z żelu (5). Wycinki żelu potraktowano proteazą K, a następnie wyekstrahowano RNA (6), przekształcono w cDNA i amplifikowano 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ć na jednym torze. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 2: Sieciowanie i filtracja próżniowa. (A) Środek sieciujący. Zawiesinę komórek wlewa się do lejka znajdującego się w prawym górnym rogu maszyny (zobacz również Rysunek 3A dla zbliżenia) i trzymany w przezroczystej dla UV torbie umieszczonej na środkowej tacy. Torba ta jest otoczona dwiema przesłonami, które pozostają zamknięte, dopóki użytkownik nie poinstruuje maszyny, aby rozpoczęła etap napromieniania. Komórki są naświetlane światłem UV z tacek zarówno powyżej, jak i poniżej. Maszyna jest dostarczana z lampami UV o długości fali 254 i 365 nm, przy czym ta ostatnia ma zastosowanie do eksperymentów PAR-CLIP. Maszyna jest obsługiwana za pomocą panelu dotykowego znajdującego się w prawym górnym rogu, który pozwala kontrolować dawkę UV lub czas ekspozycji. (B) Po usieciowaniu ogniwa są opróżniane z lewej strony maszyny. Zawiesiny komórek są odzyskiwane za pomocą próżni i odprowadzane do szklanej kolby, skąd można je następnie wlać 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, a między nim wkładany jest filtr. Cztery urządzenia filtracyjne były używane równolegle przez bardzo krótkie szeregi czasowe, aby nie tracić czasu w wyniku wymiany filtrów. (D) Po przefiltrowaniu supernatant pożywki został odsączony do kolb w celu późniejszego usunięcia. Zawory zostały zainstalowane poniżej urządzeń do filtracji próżniowej, aby utrzymać podciśnienie w systemie po wyjęciu filtra. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 3: Sieciowanie komórek zawieszonych i przylegających. (A) Środek sieciujący z komorą Vari-X-linker do komórek zawiesinowych. Kulturę komórkową wlewa się do wlotu próbki (lejka) znajdującego się w prawym górnym rogu tacy. (B) Taca, która może pomieścić plastikowe lub kwarcowe szalki Petriego do sieciowania przylegających komórek lub małych objętości komórek zawiesinowych. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 4: Przygotowanie biblioteki. (A) Przykład autoradiogramu z eksperymentu Nrd1-HTP χCRAC. Silny, skoncentrowany sygnał reprezentuje białko usieciowane z bardzo krótkimi RNA, podczas gdy powyższy rozmaz reprezentuje białko usieciowane z RNA o długości wystarczającej do sekwencjonowania. (B) Wymaz został wycięty, jak pokazano na autoradiogramie wykonanym po wycięciu żelu. (C) Reprezentatywny qPCR z biblioteki cDNA χCRAC. W tym przykładzie maksymalna amplifikacja cDNA została osiągnięta po 16 cyklach. W ten sposób do końcowego wzmocnienia użyto 16 cykli. Słupek błędu reprezentuje odchylenie standardowe trzech technicznych powtórzeń qPCR. (D) Przykład obrazu fosforowego z biblioteki cDNA na 6% żelu TBE. (E) analiza długości i jakości cDNA z elektroforezy kapilarnej opartej na chipie. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Ryc. 5: Test wysokiej RNazy w eksperymencie iCLIP w celu przetestowania sieciowania w komórkach ssaków. Pokazano autoradiogramy z eksperymentów GFP-RBM7 iCLIP, w których testowano skuteczność odzyskiwania RNP przy różnych energiach sieciowania. Immunoprecypitacje przeprowadzono przy użyciu przeciwciał anty-GFP sprzężonych z kulkami magnetycznymi na usieciowanych komórkach, które stabilnie wykazywały ekspresję GFP-RBM7. Immunoprecypaty inkubowano z wysokimi stężeniami RNazy I w celu przycięcia powiązanych RNA do krótkich, jednolitych długości. RNP wizualizowano za pomocą znakowania 32P i SDS-PAGE i migrują jako zdefiniowane pasmo, zbliżone do migracji białka nieusieciowanego. Kwantyfikacja wskazuje na wyniki analiz densytometrycznych znakowanego radioaktywnie sygnału RBM7-RNA znormalizowanego do sygnału western blot anty-GFP. (A) Przebieg czasowy sieciowania powszechnie stosowanego sieciownika UVP w porównaniu z naszym sieciownikiem (Vari-X-linker; VxL). (B) Sieciowanie przebiegu czasowego naszego sieciownika na kulturach kwarcowych (po lewej) i plastikowych (po prawej). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 6: Przykładowe wykresy przeglądarki genomu pokazujące moc χCRAC do pokazania różnicowego, czasowego wiązania Nrd1 z jego celami. Każde pole zawiera wykresy dla poszczególnych regionów genomu. Strzałki wskazują, na której nici są zakodowane geny (strzałka skierowana w lewo = nić ujemna; strzałka skierowana w prawo = nić plus). Punkty czasowe (min) są oznaczone symbolami t0, t1, t2 itd. na osiach y każdego podpowierzchni. Pokazane są cyfry rzymskie wskazujące chromosomy i współrzędne. (A) Po deprywacji glukozy, Nrd1 wiąże dwa transportery glukozy o wysokim powinowactwie, HXT6 i HXT7, które są regulowane w górę w tym stanie. (B) Zaobserwowano, że Nrd1 wiąże się z Imd3, już zatwierdzonym celem Nab344, ze zmniejszającą się intensywnością po głodzie glukozy. (C) Wiązanie Nrd1 Tye7 wykazuje dynamiczny i przejściowy charakter, zmniejszając się po głodzie glukozy do minimum po 8 minutach stresu. Jednak wiązanie powraca do poziomu podstawowego po 14 minutach. Odczyty zostały znormalizowane do "odczytów na milion" (RPM; oś y). Szare pola wskazują regiony kodujące niekodujące RNA. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.