$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Ten raport opisuje pełny protokół SLIC-CAGE do uzyskiwania bibliotek gotowych do sekwencjonowania z nanogramów początkowego całkowitego materiału RNA (Rysunek 1). Aby uzyskać mieszaninę nośników syntetycznego RNA, najpierw należy przygotować i oczyścić matryce nośników PCR w celu wyeliminowania produktów ubocznych PCR (Rysunek 2A). Każda matryca PCR (w sumie dziesięć) jest wytwarzana przy użyciu wspólnego startera do przodu, ale innego startera wstecznego (tabela 2), co prowadzi do różnych długości matrycy PCR, aby umożliwić zmienność wielkości nośników syntetycznego RNA. Po oczyszczeniu matryce PCR są używane do transkrypcji cząsteczek nośnikowych in vitro. Oczekuje się, że pojedynczy produkt nośnika RNA zostanie oczyszczony żelem (patrz reprezentatywna analiza żelu w Rysunek 2B). Przygotowanie nośnika można zwiększyć w zależności od potrzeb, a po przygotowaniu wymieszać i zamrozić w temperaturze -80 °C do wykorzystania w przyszłości.
Używając zalecanej minimalnej ilości próbki całkowitego RNA (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 wejściowego RNA (oczekiwaną liczbę cykli przedstawiono w tabeli 4).
Po pierwszej rundzie degradacji, w wynikach qPCR (krok 17), oczekiwana różnica między wartościami Ct uzyskanymi przy użyciu adaptor_f1 lub carrier_f1 startera wynosi 1-2, przy czym wartości Ct uzyskane przy adaptor_f1 niższych niż przy carrier_f1.
Rozkład długości fragmentów w końcowej bibliotece wynosi od 200 do 2 000 bp, a średni rozmiar fragmentu wynosi 700-900 pz (na podstawie analizy regionu za pomocą oprogramowania Bioanalyzer, Rysunek 4B,D). Krótsze fragmenty, jak pokazano w Rysunek 4A,C, muszą zostać usunięte przez dodatkowe rundy wykluczania rozmiaru (kroki 20-21). Te krótkie fragmenty są artefaktami amplifikacji PCR, a nie biblioteką docelową. Należy pamiętać, że krótsze fragmenty lepiej grupują się w komórkach przepływu sekwencjonowania i mogą powodować problemy z sekwencjonowaniem.
Oczekiwana ilość materiału bibliotecznego uzyskanego na próbkę wynosi od 5 do 50 ng. Znacznie niższe ilości wskazują na utratę próbki w trakcie trwania protokołu. Jeśli uzyskana niewielka ilość jest wystarczająca do sekwencjonowania (potrzebne są 2-3 ng bibliotek w puli), biblioteki mogą mieć mniejszą złożoność (patrz poniżej).
W zależności od maszyny sekwencjonującej, ilość biblioteki załadowanej do komórki przepływowej może wymagać optymalizacji. Korzystając z Illumina HiSeq 2500, ładowanie bibliotek 8-12 pM SLIC-CAGE daje średnio 150-200 milionów odczytów, przy czym >80% odczytów spełnia próg Q30.
Uzyskane odczyty są następnie mapowane do genomu referencyjnego [dla odczytów 50 pz, Bowtie212 może być używane z domyślnymi parametrami, które pozwalają na zero niezgodności na sekwencję nasion (22 pz)]. Oczekiwane wyniki mapowania zależą od całkowitej ilości wejściowego RNA i przedstawiono je w tabeli 5. Unikalnie zmapowane odczyty mogą być następnie załadowane do graficznego i statystycznego środowiska obliczeniowego R13 i przetworzone za pomocą CAGEr (Bioconductor package14). Winieta pakietu jest łatwa do naśladowania i szczegółowo wyjaśnia przepływ pracy i przetwarzanie zmapowanych danych. Łatwą wizualną kontrolą złożoności biblioteki jest rozkład szerokości promotora, ponieważ biblioteki o niskiej złożoności będą miały sztucznie zawęż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ę prawdziwych CTSS z większą precyzją niż alternatywne metody mapowania TSS o niskim/średnim wkładzie (Rysunek 5B, C).

Rysunek 1: Kroki w protokole SLIC-CAGE. Próbka RNA jest mieszana z mieszanką nośników RNA w celu uzyskania 5 μg całkowitego materiału RNA. cDNA jest syntetyzowane przez odwrotną transkrypcję, a czapeczka jest utleniana za pomocą periodynianu sodu. Utlenianie umożliwia przyłączenie biotyny do nakrętki za pomocą hydrazydu biotyny. Biotyna przyłącza się do końca 3′ mRNA, ponieważ jest również utleniana za pomocą periodynianu sodu. Aby wyeliminować biotynę z hybryd mRNA:cDNA z niecałkowicie zsyntetyzowanym cDNA i 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 wybierane przez oczyszczanie powinowactwa na kulkach magnetycznych streptawidyny (pułapkowanie czapek). Po uwolnieniu cDNA ligowane są łączniki 5′ i 3′. Cząsteczki biblioteczne, które pochodzą z nośnika, są degradowane za pomocą endonukleaz naprowadzających I-SceI i I-CeuI, a fragmenty są usuwane za pomocą kulek magnetycznych SPRI. Biblioteka jest następnie amplifikowana metodą PCR. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rycina 2: Reprezentatywna analiza żelowa matryc do PCR na nośniku i transkryptów nośnika in vitro. (A) Matryce do PCR na nośniku przed oczyszczeniem żelu: pierwsza studzienka zawiera marker 1 kbp, a następnie matryce 1 i 1-10 nośnika PCR. (B) Nośnik transkryptów in vitro: pierwsza studzienka zawiera marker 1 kbp, a następnie transkrypty nośne 1-10. Transkrypty nośników poddano denaturacji przez ogrzewanie przez 5 minut w temperaturze 95 °C przed załadowaniem. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rycina 3: Reprezentatywny ślad jakości DNA (chip DNA o wysokiej czułości) SLIC-CAGE przed pierwszą rundą degradacji nośnika. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rycina 4: Reprezentatywna jakość DNA (chip DNA o wysokiej czułości) śladów bibliotek SLIC-CAGE po amplifikacji PCR. (A) Biblioteka SLIC-CADE, która wymaga dodatkowego wyboru rozmiaru w celu usunięcia krótkich fragmentów. (B) Biblioteka SLIC-CAGE po wyborze rozmiaru przy użyciu 0,6x proporcji kulek SPRI do próbki. (C) Biblioteka SLIC-CAGE o mniejszej ilości wyjściowej, która wymaga wyboru rozmiaru w celu usunięcia krótkiego fragmentu. (D) Biblioteka SLIC-CAGE o niższej ilości wyjściowej po wybraniu rozmiaru przy użyciu proporcji kulek SPRI 0,6:1 do próbki. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 5: Walidacja bibliotek SLIC-CAGE. (A) Rozkład szerokości międzykwantylowych klastrów znacznikó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. Duża liczba wąskich klastrów znaczników w bibliotece 1 ng SLIC-CAGE wskazuje na jej niską złożoność. (B) Krzywe ROC do identyfikacji CTSS w bibliotekach S. cerevisiae SLIC-CADE. Wszystkie CTSS S. cerevisiae nAnT-iCAGE zostały użyte jako prawdziwy zestaw. (C) Krzywe ROC do identyfikacji CTSS w bibliotekach S. cerevisiae nanoCAGE. Wszystkie CTSS S. cerevisiae nAnT-iCAGE zostały użyte jako prawdziwy zestaw. 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 zobaczyć większą wersję tego rysunku.
Tabela 1: Sekwencja syntetycznego genu nosiciela. Witryny rozpoznawania I-SceI są pogrubione i oznaczone kursywą na fioletowo, a miejsca rozpoznawania I-CeuI są zielone. Kliknij tutaj, aby wyświetlić tę tabelę (kliknij prawym przyciskiem myszy, aby pobrać).
TGL
szt.
szt.
szt.
TGL
TGLI
pkt.
TGL
TGL
| przewoźnik | podkład odwrócony 5'-3' | Długość produktu PCR / bp |
| 1 | PCR_N6_r1: NNNNNNCTACGTGTCGCAGACGAATT | Numer katalogowy 1034 |
| cyfra arabska | 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 | Okręg wyborczy 531 |
| 9 | PCR_N6_r9: NNNNNNACCGCCGCGCCTTCCGCAGG | Rejon 454 |
| 10 | PCR_N6_r10: NNNNNNCAGGACGTTTTTGCCCAGCA | Okręg wyborczy 386 |
| * Podkład do przodu jest taki sam dla wszystkich szablonów nośników. Podkreślona jest sekwencja promotora T7. PCR_GN5_f1: TAATACGACTCACTATAGNNNNNCAGCGTTCGCTA | |
Tabela 2: Startery do wzmacniania szablonu nośnego. Podkład do przodu jest taki sam dla wszystkich szablonów nośników. Podkreślona jest sekwencja promotora T7. PCR_GN5_f1: TAATACGACTCACTATAGNNNNNCAGCGTTCGCTA. Przy użyciu różnych starterów odwrotnych wytwarza się matryce PCR, a tym samym nośne RNA o różnej długości.
TGL
pkt.
pkt.
szt.
pkt.
szt.
szt.
pkt.
TGL
TGLI
pkt.
pkt.
TGL
pkt.
pkt.
TGL
TGL
pkt.
| przewoźnik | długość | bez nakrętki/μg | Ograniczony/μg |
| 1 | Numer katalogowy 1034 | 3,96 | 0,45 |
| cyfra arabska | 966 | Godzina 8,36 | 0,95 |
| 3 | 889 | 4.4 | 0,5 |
| 4 | 821 | 6,6 | 0,75 |
| 5 | 744 | 4.4 | 0,5 |
| 6 | 676 | Pytanie 3,08 | 0,35 |
| 7 | 599 | 4.4 | 0,5 |
| 8 | Okręg wyborczy 531 | 3,96 | 0,45 |
| 9 | Rejon 454 | Z drugiej strony, 64 | 0,3 |
| 10 | Okręg wyborczy 386 | 2.2 | 0,25 |
Tabela 3: Mieszanka nośników RNA. Łącznie 49 μg mieszaniny nośników 0,3-1 kbp: bez korka = 44 μg, z nakrętką = 5 μg.
| Całkowity wkład RNA /ng | Cykle PCR |
| 1 ng | Rozdział 18 |
| 2 ng | 17 |
| 5 ngów | 16 |
| 10 ng | 15-16 |
| 25 ng | 14-15 |
| 50 ng | 13-15 |
| 100 ng | 12-14 |
Tabela 4: Oczekiwana liczba cykli PCR w zależności od całkowitej ilości próbki pobranej RNA. Przybliżona liczba cykli opiera się na eksperymentach przeprowadzonych przy użyciu Saccharomyces cerevisiae, Drosophila melanogaster i Mus musculus total RNA.
Rozdział
szt.
szt.
szt.
szt.
szt.
szt.
| Całkowity dopływ RNA/ng | % ogółem zmapowanych | % unikatowych odwzorowań | % przewoźnika |
| 1 ng | Rozdział 30 | 20-30 | Rozdział 30 |
| 2 ng | 60 | 20-50 | 10 |
| 5 ngów | 60-70 | 40-60 | 5-10 |
| 10 ng | 60-70 | 40-60 | 5-10 |
| 25 ng | 65-80 | Rozdział 40-70 | 0-5 |
| 50 ng | 65-80 | Rozdział 40-70 | 0-3 |
| 100 ng | 70-85 | Rozdział 40-70 | 0-2 |
Tabela 5: Oczekiwana skuteczność mapowania i zależność od całkowitej ilości wejściowego RNA. Przybliżone liczby są przedstawione i oparte na eksperymentach przeprowadzonych przy użyciu Saccharomyces cerevisiae i Mus musculus total RNA.
Tabela uzupełniająca 1: Sekwencje podkładowe. Kliknij tutaj, aby wyświetlić tę tabelę (kliknij prawym przyciskiem myszy, aby pobrać).
Tabela uzupełniająca 2: Wyżarzanie łączników 5' i 3'. Kliknij tutaj, aby wyświetlić tę tabelę (kliknij prawym przyciskiem myszy, aby pobrać).
Tabela uzupełniająca 3: 5' wyżarzanie łącznika. Kliknij tutaj, aby wyświetlić tę tabelę (kliknij prawym przyciskiem myszy, aby pobrać).
Tabela uzupełniająca 4: 5' mieszanie łączników. Kliknij tutaj, aby wyświetlić tę tabelę (kliknij prawym przyciskiem myszy, aby pobrać).
Tabela uzupełniająca 5: 5' rozcieńczenie łącznika. Kliknij tutaj, aby wyświetlić tę tabelę (kliknij prawym przyciskiem myszy, aby pobrać).
Tabela uzupełniająca 6: 3' wyżarzanie łącznika. Kliknij tutaj, aby wyświetlić tę tabelę (kliknij prawym przyciskiem myszy, aby pobrać).
Tabela uzupełniająca 7: 3' rozcieńczenie łącznika. Kliknij tutaj, aby wyświetlić tę tabelę (kliknij prawym przyciskiem myszy, aby pobrać).
Tabela uzupełniająca 8: Przygotowanie standardowych seryjnych rozcieńczeń. Kliknij tutaj, aby wyświetlić tę tabelę (kliknij prawym przyciskiem myszy, aby pobrać).