Niniejszy raport opisuje pełny protokół SLIC-CAGE służący do otrzymywania bibliotek gotowych do sekwencjonowania z nanogramowych ilości wyjściowego całkowitego RNA (Rysunek 1). Aby uzyskać mieszaninę syntetycznego nośnika RNA, najpierw należy przygotować matryce nośnika PCR i oczyścić je za pomocą elektroforezy żelowej w celu usunięcia produktów ubocznych PCR (Rysunek 2A). Każda matryca do PCR (łącznie dziesięć) jest otrzymywana przy użyciu wspólnego prawego, lecz różnych lewych starterów (Tabela 2), co prowadzi do uzyskania matryc PCR o różnych długościach, aby umożliwić zróżnicowanie rozmiaru syntetycznych nośników RNA. Po oczyszczeniu matryce PCR są wykorzystywane do transkrypcji in vitro cząsteczek nośnika. Oczekiwanym rezultatem jest pojedynczy produkt nośnika RNA, jeżeli matryce zostały oczyszczone za pomocą elektroforezy żelowej (patrz reprezentacyjna analiza żelowa w Rycina 2B). Przygotowanie nośnika można przeskalować w zależności od potrzeb, a po przygotowaniu wymieszać i zamrozić w temperaturze -80 °C do późniejszego wykorzystania.
Stosując zalecaną minimalną ilość całkowitego RNA z próbki (10 ng) w połączeniu z 16-18 cyklami amplifikacji PCR, można uzyskać biblioteki SLIC-CAGE o wysokiej złożoności. Liczba cykli PCR wymaganych do amplifikacji końcowej biblioteki w dużym stopniu zależy od ilości użytego całkowitego RNA wejściowego (oczekiwana liczba cykli jest przedstawiona w Tabeli 4).
Po pierwszej rundzie degradacji, w wynikach qPCR (krok 17), oczekiwana różnica między wartościami Ct uzyskanymi przy użyciu startera adaptor_f1 a carrier_f1 wynosi 1-2, przy czym wartości Ct uzyskane dla adaptor_f1 są niższe niż dla carrier_f1.
Rozkład długości fragmentów w końcowej bibliotece wynosi od 200 do 2 000 bp, przy średniej wielkości fragmentu od 700 do 900 bp (na podstawie analizy regionów przy użyciu oprogramowania Bioanalyzer, Rysunek 4B,D). Krótsze fragmenty, przedstawione na Rysunku 4A,C, muszą zostać usunięte poprzez dodatkowe cykle wykluczania wielkością (kroki 20-21). Krótkie fragmenty te są artefaktami amplifikacji PCR, a nie docelową biblioteką. Należy pamiętać, że krótsze fragmenty lepiej skupiają się na przepływowych komórkach sekwencjonujących i mogą powodować problemy z sekwencjonowaniem.
Oczekiwana ilość materiału biblioteki uzyskanego z jednej próbki wynosi od 5-50 ng. Znacznie niższe ilości wskazują na stratę próbki podczas protokołu. Jeśli uzyskana niska ilość jest wystarczająca do sekwencjonowania (wymagane jest 2-3 ng połączonych bibliotek), biblioteki mogą charakteryzować się mniejszą złożonością (patrz poniżej).
W zależności od sekwenatora konieczna może być optymalizacja ilości biblioteki załadowanej na ogniwo przepływowe. W przypadku zastosowania systemu Illumina HiSeq 2500, załadowanie bibliotek SLIC-CAGE o stężeniu 8-12 pM pozwala uzyskać średnio 150-200 milionów odczytów, z czego >80% odczytów spełnia próg jakości Q30.
Uzyskane odczyty są następnie mapowane do genomu referencyjnego [dla odczytów 50 bp można zastosować program Bowtie212 z parametrami domyślnymi, które dopuszczają zero niedopasowań na sekwencję ziarna (22 bp)]. Oczekiwana wydajność mapowania zależy od całkowitej ilości wprowadzonego RNA i jest przedstawiona w Tabeli 5. Unikalnie zmapowane odczyty można następnie załadować do środowiska obliczeń graficznych i statystycznych R13 i przetworzyć przy użyciu pakietu CAGEr (pakiet Bioconductor14). Winieta pakietu jest łatwa do śledzenia i szczegółowo wyjaśnia przebieg prac oraz przetwarzanie zmapowanych danych. Prostym sposobem wizualnej kontroli złożoności biblioteki jest rozkład szerokości promotora, ponieważ biblioteki o niskiej złożoności będą miały sztucznie zwężone promotory (Rysunek 5A, biblioteka SLIC-CAGE pochodząca z 1 ng całkowitego RNA, szczegóły w poprzedniej publikacji10). Jednak nawet biblioteki SLIC-CAGE o niskiej złożoności pozwalają na identyfikację rzeczywistych CTSS z większą precyzją niż alternatywne metody mapowania TSS przy niskiej/średniej ilości materiału wejściowego (Rysunek 5B,C).

Rysunek 1: Etapy protokołu SLIC-CAGE. RNA z próbki miesza się z mieszaniną nośnika RNA, aby uzyskać 5 µg całkowitego materiału RNA. cDNA jest syntetyzowane poprzez transkrypcję odwrotną, a czapeczka zostaje utleniona przy użyciu nadjodynianu sodu. Utlenianie umożliwia przyłączenie biotyny do czapeczki za pomocą hydrazydy biotyny. Biotyna przyłącza się również do końca 3′ mRNA, ponieważ ten również jest utleniany nadjodynianem sodu. Aby usunąć biotynę z hybryd mRNA:cDNA z niepełnie zsyntetyzowanym cDNA oraz z końców 3′ mRNA, próbki traktuje się RNazą I. cDNA, które dotarło do końca 5′ mRNA, jest następnie wyodrębniane poprzez oczyszczanie powinowactwa na magnetycznych koralikach z streptawidyną (wychwytywanie czapeczki). Po uwolnieniu cDNA liguje się linkery 5′ oraz 3′. Cząsteczki biblioteki pochodzące z nośnika są degradowane za pomocą endonukleaz homingowych I-SceI i I-CeuI, a fragmenty są usuwane przy użyciu magnetycznych koralików SPRI. Następnie biblioteka jest amplifikowana metodą PCR. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 2: Reprezentatywna analiza żelowa matryc PCR nośnikowych i transkryptów nośnikowych in vitro. (A) Matryce PCR nośnikowe przed oczyszczaniem żelowym: pierwsza studnia zawiera marker 1 kbp, a następnie matryce PCR nośnikowe 1, 1-10. (B) Transkrypty nośnikowe in vitro: pierwsza studnia zawiera marker 1 kbp, a następnie transkrypty nośnikowe 1-10. Transkrypty nośnikowe zostały zdenaturowane przez ogrzewanie przez 5 min w 95 °C przed nałożeniem. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 3: Reprezentatywny zapis jakości DNA (chip DNA o wysokiej czułości) SLIC-CAGE przed pierwszą rundą degradacji nośnika. Aby zobaczyć powiększoną wersję tego rysunku, kliknij tutaj.

Rysunek 4: Reprezentatywne profile jakości DNA (chip DNA o wysokiej czułości) bibliotek SLIC-CAGE po amplifikacji PCR. (A) Biblioteka SLIC-CAGE wymagająca dodatkowej selekcji wielkości w celu usunięcia krótkich fragmentów. (B) Biblioteka SLIC-CAGE po selekcji wielkości z zastosowaniem kulek SPRI w stosunku 0,6x do próbki. (C) Biblioteka SLIC-CAGE o niższej wydajności, wymagająca selekcji wielkości w celu usunięcia krótkich fragmentów. (D) Biblioteka SLIC-CAGE o niższej wydajności po selekcji wielkości z zastosowaniem kulek SPRI w stosunku 0,6:1 do próbki. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 5: Walidacja bibliotek SLIC-CAGE. (A) Rozkład szerokości międzykwartylowej klastrów tagów w bibliotekach SLIC-CAGE przygotowanych z 1, 5 lub 10 ng całkowitego RNA S. cerevisiae oraz w bibliotece nAnT-iCAGE przygotowanej z 5 µg całkowitego RNA S. cerevisiae. Wysoka liczba wąskich klastrów tagów w bibliotece SLIC-CAGE 1 ng wskazuje na jej niską złożoność. (B) Krzywe ROC dla identyfikacji CTSS w bibliotekach SLIC-CAGE S. cerevisiae. Wszystkie CTSS nAnT-iCAGE S. cerevisiae zostały użyte jako zbiór prawdziwy. (C) Krzywe ROC dla identyfikacji CTSS w bibliotekach nanoCAGE S. cerevisiae. Wszystkie CTSS nAnT-iCAGE S. cerevisiae zostały użyte jako zbiór prawdziwy. Porównanie krzywych ROC pokazuje, że SLIC-CAGE znacznie przewyższa nanoCAGE w identyfikacji CTSS. Wykorzystano dane z ArrayExpress E-MTAB-6519. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Tabela 1: Sekwencja syntetycznego genu nośnika. Miejsca I-SceI są pogrubione, zapisane kursywą i zaznaczone kolorem fioletowym, a miejsca rozpoznawania I-CeuI są zielone. Kliknij tutaj, aby wyświetlić tę tabelę (kliknij prawym przyciskiem myszy, aby pobrać).
| nośnik | primer rewersyjny 5’-3’ | długość produktu PCR / bp |
| 1 | PCR_N6_r1: NNNNNNCTACGTGTCGCAGACGAATT | 1034 |
| 2 | PCR_N6_r2: NNNNNNTATCCAGATCGTTGAGCTGC | 966 |
| 3 | PCR_N6_r3: NNNNNNCACTGCGGGATCTCTTTACG | 889 |
| 4 | PCR_N6_r4: NNNNNNGCCGTCGATAACTTGTTCGT | 821 |
| 5 | PCR_N6_r5: NNNNNNAGTTGACCGCAGAAGTCTTC | 744 |
| 6 | PCR_N6_r6: NNNNNNGTGAAGAATTTCTGTTCCCA | 676 |
| 7 | PCR_N6_r7: NNNNNNCTCGCGGCTCCAGTCATAAC | 599 |
| 8 | PCR_N6_r8: NNNNNNTATACGCGATGTTGTCGTAC | 531 |
| 9 | PCR_N6_r9: NNNNNNACCGCCGCGCCTTCCGCAGG | 454 |
| 10 | PCR_N6_r10: NNNNNNCAGGACGTTTTTGCCCAGCA | 386 |
| * Primer foremowy jest taki sam dla wszystkich matryc nośnika. Podkreślono sekwencję promotora T7. PCR_GN5_f1: TAATACGACTCACTATAGNNNNNCAGCGTTCGCTA | |
Tabela 2: Primery do amplifikacji matrycy nośnikowej. Primer forward jest taki sam dla wszystkich matryc nośnikowych. Podkreślono sekwencję promotora T7. PCR_GN5_f1: TAATACGACTCACTATAGNNNNNCAGCGTTCGCTA. Dzięki zastosowaniu różnych primerów reverse otrzymuje się matryce PCR, a w konsekwencji RNA nośnikowe o różnej długości.
| nośnik | długość | bez zakrętki/µg | zamknięte/µg |
| 1 | 1034 | 3.96 | 0.45 |
| 2 | 966 | 8.36 | 0.95 |
| 3 | 889 | 4.4 | 0.5 |
| 4 | 821 | 6.6 | 0.75 |
| 5 | 744 | 4.4 | 0.5 |
| 6 | 676 | 3.08 | 0.35 |
| 7 | 599 | 4.4 | 0.5 |
| 8 | 531 | 3.96 | 0.45 |
| 9 | 454 | 2.64 | 0.3 |
| 10 | 386 | 2.2 | 0.25 |
Tabela 3: Mieszanina nośnika RNA. Łącznie 49 µg mieszaniny nośnika 0,3-1 kbp: niekapszowana = 44 µg, kapszowana = 5 µg.
| Całkowita ilość wprowadzonego RNA /ng | Cykle PCR |
| 1 ng | 18 |
| 2 ng | 17 |
| 5 ng | 16 |
| 10 ng | 15-16 |
| 25 ng | 14-15 |
| 50 ng | 13-15 |
| 100 ng | 12-14 |
Tabela 4: Przewidywana liczba cykli PCR w zależności od całkowitej ilości wprowadzonego RNA z próbki. Przybliżona liczba cykli opiera się na eksperymentach przeprowadzonych z użyciem całkowitego RNA z Saccharomyces cerevisiae, Drosophila melanogaster oraz Mus musculus.
| Całkowita ilość wprowadzonego RNA/ng | % całkowitego mapowania | % mapowania unikalnego | % nośnika |
| 1 ng | 30 | 20-30 | 30 |
| 2 ng | 60 | 20-50 | 10 |
| 5 ng | 60-70 | 40-60 | 5-10 |
| 10 ng | 60-70 | 40-60 | 5-10 |
| 25 ng | 65-80 | 40-70 | 0-5 |
| 50 ng | 65-80 | 40-70 | 0-3 |
| 100 ng | 70-85 | 40-70 | 0-2 |
Tabela 5: Przewidywana wydajność mapowania w zależności od całkowitej ilości wprowadzonego RNA. Przedstawiono wartości przybliżone, oparte na eksperymentach przeprowadzonych z wykorzystaniem całkowitego RNA z Saccharomyces cerevisiae i Mus musculus.
Tabela uzupełniająca 1: Sekwencje starterów. Kliknij tutaj, aby wyświetlić tę tabelę (kliknij prawym przyciskiem myszy, aby pobrać).
Tabela uzupełniająca 2: Hybrydyzacja linkerów 5’ i 3’. Aby wyświetlić tę tabelę, kliknij tutaj (kliknij prawym przyciskiem myszy, aby pobrać).
Tabela uzupełniająca 3: Wyżarzanie linkera 5’. Kliknij tutaj, aby wyświetlić tę tabelę (kliknij prawym przyciskiem myszy, aby pobrać).
Tabela uzupełniająca 4: Mieszanie linkerów 5’. Kliknij tutaj, aby wyświetlić tę tabelę (kliknij prawym przyciskiem myszy, aby pobrać).
Tabela uzupełniająca 5: Rozcieńczenie łącznika 5’. Kliknij tutaj, aby wyświetlić tę tabelę (kliknij prawym przyciskiem myszy, aby pobrać).
Tabela uzupełniająca 6: Wyżarzanie łącznika 3’. Kliknij tutaj, aby wyświetlić tę tabelę (Kliknij prawym przyciskiem myszy, aby pobrać).
Tabela uzupełniająca 7: Rozcieńczenie łącznika 3’. Kliknij tutaj, aby wyświetlić tę tabelę (kliknij prawym przyciskiem myszy, aby pobrać).
Tabela uzupełniająca 8: Przygotowanie standardowych rozcieńczeń seryjnych. Kliknij tutaj, aby wyświetlić tę tabelę (kliknij prawym przyciskiem myszy, aby pobrać).