Technika opisana w niniejszym artykule została zastosowana w szerszym badaniu mającym na celu wyjaśnienie mechanizmów leżących u podstaw inhibicji odwrotnej transkrypcji HIV-1 przez ludzkie białko antyretrowirusowe APOBEC3G (A3G)6. Rysunek 3 przedstawia reprezentatywne wyniki uzyskane po zastosowaniu protokołu w próbkach limfocytów T CEM-SS zainfekowanych wirusem HIV-1 z deficytem vif w obecności lub przy braku A3G. Całkowita liczba unikalnych odczytów uzyskanych z każdej próbki po odfiltrowaniu duplikatów PCR posiadających ten sam 6 nt kod kreskowy oraz tę samą długość (wykonane przez dostarczone oprogramowanie analityczne) została przedstawiona na Rysunku 3a. Wzrost poziomu A3G redukuje całkowitą liczbę odczytów, co odzwierciedla hamujący wpływ A3G na syntezę cDNA zapośredniczoną przez RT, opisany i zmierzony wcześniej za pomocą qPCR6,13,21,22. Na Rysunku 3b przedstawiono frakcję cząsteczek o każdej możliwej długości w obrębie pierwszych 182 nt. W przypadku infekcji HIV-1 przy braku A3G najliczniejszym gatunkiem jest sama główna cząsteczka strong-stop o długości 180 nt, z pewną akumulacją odczytów w krótszym zakresie (23 do 40 nt) (wykres górny, histogramy niebieskie). Dodanie A3G zmienia ten profil, gdyż wykrywany jest gwałtowny wzrost liczby krótszych, utrąconych cząsteczek cDNA w kilku bardzo specyficznych, powtarzalnych pozycjach (wykresy środkowy i dolny). Ponieważ A3G jest cytidyno-deaminazą, mutacje cytozyny na urydynę (identyfikowane jako C-to-T) w cDNA występują, gdy A3G jest obecne w zakaźnych wirionach21,23,24. Wykorzystując uzyskane informacje z sekwencjonowania, procent mutacji C-to-T naniesiono na ten sam wykres (czerwona linia przerywana). Należy zauważyć, że profil mutacyjny został wyprowadzony ze wszystkich połączonych unikalnych odczytów, a pokrycie każdego nukleotydu będzie się różnić. Jednakże, jeśli jest to wymagane, informacje o sekwencji mogą zostać powiązane z każdą cząsteczką i skorelowane z konkretnym końcem 3'. Przedstawione dane pochodzą z pracy Pollpeter et al.6, a wykazano, że korelacja między profilami mutacyjnymi a długością cDNA wynika z wykrywania i rozcinania deaminowanego cDNA przez komórkowy aparat naprawy DNA.
Kontrolę pozytywną dla metody mapowania końców 3' można łatwo przygotować, przetwarzając pulę syntetycznych oligonukleotydów o znanej sekwencji, długości i stężeniu. Kontrolę tę dodaje się podczas ligacji adapterów w kroku 3.3.2 i zaleca się jej uwzględnienie we wszystkich bibliotekach multipleksowanych. Dane uzyskane z próbki kontrolnej powinny wykazać wszystkie oligonukleotydy w oczekiwanych stosunkach wejściowych, przy bardzo niskim poziomie odczytów tła. Rycina 4 przedstawia wyniki zestawu kontroli pozytywnej składającego się z 17 chemicznie syntetyzowanych oligonukleotydów (sekwencje podano w Tabeli 2), które zostały zmieszane w stosunkach ekimolarnych. Zgodnie z oczekiwaniami, wszystkie cząsteczki występują w niemal jednakowej liczebności z jedynie niewielkimi odchyleniami (wykres górny). Podczas gdy większość pozycji w sekwencji DNA -sss, które nie były reprezentowane przez oligonukleotyd, zwraca zerową liczbę odczytów, zaobserwowaliśmy niewielką ilość produktów, które są o 1 lub 2 nt krótsze od rzeczywistych oligonukleotydów kontrolnych. Nie badaliśmy dalej tych produktów pobocznych, przyjmując, że reprezentują one produkty zdegradowane lub niekompletne, które potencjalnie znajdowały się w dostarczonych zapasach oligonukleotydów w momencie zakupu (oligonukleotydy zamówiono jako oczyszczone metodą HPLC, dla których producent wskazał czystość > 80%). Wykres dolny przedstawia próbkę kontrolną z innego przebiegu biblioteki, gdzie zmienność między 17 oligonukleotydami jest nieco większa i koreluje z całkowitą długością – dłuższe cząsteczki kontrolne są wykrywane wydajniej niż krótsze. Może to wynikać z niewielkiego błędu (bias) w reakcjach PCR lub podczas klastrowania w sekwencjonowaniu MiSeq, które posiada optymalną wielkość insertu i może wystąpić w bibliotekach o szczególnie szerokim zakresie wielkości insertów. Podstawowym sposobem na skorygowanie tego błędu jest zastosowanie czynnika normalizacji opartego na nachyleniu linii, która wskazuje błąd korelujący z długością cząsteczki (różowa linia). Wymagane obliczenia są zawarte w programie analitycznym (patrz krok 8.3 w protokole).

Rycina 1: Schemat przedstawiający pierwsze etapy odwrotnej transkrypcji HIV-1. Proces rozpoczyna się od hybrydyzacji tRNA(Lys,3) (pomarańczowy) z miejscem wiązania startera (PBS) w wirusowym RNA genomicznym (krok 1), co umożliwia inicjację i elongację wirusowego cDNA (niebieski, krok 2). Jednocześnie matrycowe RNA genomiczne jest degradowane przez aktywność RNaseH odwrotnej transkryptazy (RT) (krok 3). Pierwszym pełnym produktem pośrednim w procesie odwrotnej transkrypcji jest cDNA silnego stopu nici minus (-)sss, który powstaje, gdy polimeryzacja katalizowana przez RT osiąga koniec 5' regionu powtórzeń (R) gRNA (krok 3). Pośrednik (-)sss jest przenoszony na koniec 3' matrycy RNA genomicznego poprzez hybrydyzację z komplementarnym regionem R długiego powtórzenia końcowego (LTR) 3'. Stąd polimeryzacja jest kontynuowana (krok 4). W opisanym metodzie postęp odwrotnej transkrypcji jest określany poprzez mapowanie dokładnej długości nowo powstałego wirusowego cDNA (niebieski). PPT, trakt polipurynowy; U5, unikalna sekwencja 5'; U3, unikalna sekwencja 3'. Rycina została przedrukowana z wcześniejszej publikacji6. Kliknij tutaj, aby zobaczyć powiększoną wersję tej ryciny.

Rysunek 2: Schemat przebiegu prac oraz schemat strategii ligacji adapterów i amplifikacji PCR. (a) Schemat przebiegu prac przedstawiający główne etapy opisanej techniki służącej do wyznaczania końców 3' transkryptów odwrotnych HIV-1 w zainfekowanych komórkach. Rysunek został zaadaptowany z wcześniejszej publikacji6. (b) Schemat strategii ligacji adapterów i amplifikacji PCR. Nowo powstałe cząsteczki cDNA o różnej długości, które zostały oczyszczone w poprzednich etapach, są ligowane z jednoniciowym adapterem DNA przy użyciu ligazy T4 DNA. Projekt adaptera w formie spinki do włosów (nazwanego „full Kwok + MiSeq”, patrz Tabela 1) był inspirowany pracami Kwok et al.11. Adapter zawiera losową 6 nt sekwencję kodu kreskowego (barcode), która umożliwia parowanie zasad ułatwiające ligację i jednocześnie służy jako identyfikator dla unikalnych odczytów. Końce 3' adaptera posiadają łącznik (SpC3) w celu zapobiegania samoligacji. Produkty ligacji są oddzielane od nadmiaru adaptera za pomocą denaturującej elektroforezy w żelu poliakryloamidowym (PAGE). Kwasy nukleinowe w żelu są barwione i wycinane w trzech oddzielnych, równych fragmentach żelu w obszarze od poziomu powyżej adaptera do studzienki, zgodnie z metodą opisaną w25. Po elucji, wytrąceniu i resuspensji, produkty są amplifikowane metodą PCR z użyciem starterów przyłączających się do znanej sekwencji adaptera (starter 1, zestaw oligonukleotydów multipleksowych, patrz Tabela materiałów) oraz startera zawierającego pierwsze 22 nt sekwencji 5'-LTR HIV-1 bezpośrednio następującej po tRNA (starter 2, MP1.0 + 22HIV). Końce 5' wybranych starterów zawierają adaptery dla wybranej platformy sekwencjonowania (P5 i P7) oraz sekwencję indeksową w celu rozróżnienia poszczególnych próbek analizowanych w tej samej bibliotece. Oznaczono punkty startowe starterów odczytu sekwencjonowania. Niebieska ramka wskazuje region zainteresowania w celu wyznaczenia oryginalnych końców 3' uchwyconej cząsteczki. Rysunek ten został zaadaptowany z wcześniejszej publikacji6. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 3Reprezentatywne wyniki. (a) Całkowita liczba odczytów reprezentatywnych próbek przetworzonych zgodnie z opisanym protokołem. Obejmuje to wszystkie sekwencje zidentyfikowane jako unikalne odczyty cząsteczek HIV-1, których końce 3' znajdowały się w obrębie pierwszych 635 nt cDNA nici minus (do PPT, patrz Rycina 1). Zakażenie wirusem HIV-1 nieprzenoszącym A3G daje najwyższą liczbę odczytów, podczas gdy A3G hamuje syntezę cDNA, a tym samym zmniejsza całkowitą liczbę odczytów. Niezakażone komórki posłużyły jako kontrola negatywna, natomiast zestaw syntetycznych oligonukleotydów stanowi kontrolę pozytywną. b) Względna liczebność cDNA dla każdej długości między pozycjami nukleotydowymi 23 a 182 (pełnowymiarowe cDNA sss ma od 180 do 182 nt) HIV-1NL4.3 sekwencja (oś X) jest przedstawiona na niebieskich histogramach (skala na lewej osi Y). Względną obfitość cDNA obliczono, dzieląc bezwzględną liczbę sekwencji kończących się na danym nukleotydzie w obrębie sekwencji cDNA -sss przez sumę wszystkich odczytów o długości 182 nt lub mniejszej. Czerwone linie przerywane przedstawiają procent odczytów z mutacjami C-na-T/U w odpowiedniej pozycji (skala na prawej osi Y)-osie). Rycina 3b zostało przedrukowane z poprzedniej publikacji6. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 4: Reprezentatywne wyniki próbek kontrolnych. Przedstawiono dwa profile dla pul zawierających ekimolowe ilości 17 syntetycznych oligonukleotydów o różnych długościach. Oligonukleotydy te mają sekwencje pochodzące z HIV-1NL4.3 i zostały dobrane tak, aby obejmowały różne długości oraz prezentowały wszystkie 4 zasady jako nukleotyd 3'-końcowy (patrz Tabela 2). Górny wykres przedstawia próbkę kontrolną pozytywną z Rysunku 3a. Nie stwierdzono istotnego odchylenia w zależności od długości cząsteczek lub otwartych końców 3'. Dolny wykres przedstawia inny przebieg sekwencjonowania biblioteki, który spowodował niewielkie odchylenie długości w sekwencjonowaniu. W takim przypadku zaleca się zastosowanie współczynnika normalizacji, który jest wyprowadzony z nachylenia (zaznaczonego na różowo) reprezentującego odchylenie wielkości. Rysunek ten został przedrukowany z wcześniejszej publikacji6. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.
| Nazwa oligonukleotydu | Długość w nt | Sekwencja | Cel | Producent (oczyszczanie) |
| pełny Kwok + MiSeq | 61 | 5’-PHO-tgaagagcctagtcgctgttcannnnnnctgcccatagagagatcggaagagcacacgtct-SpC3-3’ | Adapter | Technologie DNA IDT (HPLC) |
| Przynęta 2xBiotin SS | 40 | 5’-biotyna-cagtgtggaaaatctctagcagtggcgcccgaacagggac-biotyna-3’ | Wzbogacanie hybrydyzacyjne | MWG Eurofins (HPLC) |
| Biotyna 1-16 ss | 22 | 5’-cagtgtggaaaatctctagcag-BiTEG-3’ | Hybrydyzacyjne wychwytywanie | MWG Eurofins (HPLC) |
| Biotynilowany tRNA + CTG | 16 | 5’-cagtggcgcccgaaca-BITEG-3’ | Hybrydowe wychwytywanie | MWG Eurofins (HPLC) |
| MP1.0 + 22HIV | 82 | 5’-aatgatacggcgaccaccgagatctacactctttccctacacgacgctcttccgatctcactgctagagattttccacactg-3’ | Amplifikacja PCR | MWG Eurofins (HPLC |
Tabela 1: Tabela oligonukleotydów zawierająca długość, sekwencje oraz modyfikacje wykorzystywane w opisanym protokole. Tabela została opracowana na podstawie wcześniejszej publikacji6. Kliknij tutaj, aby pobrać tę tabelę jako plik excel.
| Nazwa oligonu | Długość w nt | Sekwencja | Producent (Oczyszczanie) |
| HTP con long C | 120 | 5’-ctgctagagattttccacactgactaaaagggtctgagggatctctagttaccagagtcacacaacagacgggcacacactactttgagcactcaaggcaagctttattgaggcttaagc-3’ | MWG Eurofins (HPLC) |
| HTP con long G | 119 | 5’-ctgctagagattttccacactgactaaaagggtctgagggatctctagttaccagagtcacacaacagacgggcacacactactttgagcactcaaggcaagctttattgaggcttaag-3’ | MWG Eurofins (HPLC) |
| HTP con long T | 116 | 5’-ctgctagagattttccacactgactaaaagggtctgagggatctctagttaccagagtcacacaacagacgggcacacactactttgagcactcaaggcaagctttattgaggctt-3’ | MWG Eurofins (HPLC) |
| HTP con long A | 118 | 5’-ctgctagagattttccacactgactaaaagggtctgagggatctctagttaccagagtcacacaacagacgggcacacactactttgagcactcaaggcaagctttattgaggcttaa-3’ | MWG Eurofins (HPLC) |
| HTP con mid C | 76 | 5’-ctgctagagattttccacactgactaaaagggtctgagggatctctagttaccagagtcacacaacagacgggcac-3’ | MWG Eurofins (HPLC) |
| HTP con mid G (a) | 71 | 5’-ctgctagagattttccacactgactaaaagggtctgagggatctctagttaccagagtcacacaacagacg-3’ | MWG Eurofins (HPLC) |
| HTP con mid G (b) | 72 | 5’-ctgctagagattttccacactgactaaaagggtctgagggatctctagttaccagagtcacacaacagacgg-3’ | MWG Eurofins (HPLC) |
| HTP con mid A | 69 | 5’-ctgctagagattttccacactgactaaaagggtctgagggatctctagttaccagagtcacacaacaga-3’ | MWG Eurofins (HPLC) |
| HTP con mid T | 85 | 5’-ctgctagagattttccacactgactaaaagggtctgagggatctctagttaccagagtcacacaacagacgggcacacactactt-3’ | MWG Eurofins (HPLC) |
| HTP con short A | 40 | 5’-ctgctagagattttccacactgactaaaagggtctgaggga-3’ | MWG Eurofins (HPLC) |
| HTP con short T | 33 | 5’-ctgctagagattttccacactgactaaaagggt-3’ | MWG Eurofins (HPLC) |
| HTP con short G | 41 | 5’-ctgctagagattttccacactgactaaaagggtctgaggg-3’ | MWG Eurofins (HPLC) |
| HTP con short C | 34 | 5’-ctgctagagattttccacactgactaaaagggtc-3’ | MWG Eurofins (HPLC) |
| HTP Con 46 (T) | 46 | 5’-ctgctagagattttccacactg actaaaagggtctgagggatctct-3’ | MWG Eurofins (HPLC) |
| HTP Con 83 (C) | 83 | 5’-ctgctagagattttccacactg actaaaagggtctgagggatctctagttaccagagtcacacaacagacgggcacacactac-3’ | MWG Eurofins (HPLC) |
| HTP Con 103 (C) | 103 | 5’-ctgctagagattttccacactg actaaaagggtctgagggatctctagttaccagagtcacacaacagacgggcacacactactttgagcactcaaggcaagc-3’ | MWG Eurofins (HPLC) |
| HTP Con 107 (A) | 107 | 5’-ctgctagagattttccacactg actaaaagggtctgagggatctctagttaccagagtcacacaacagacgggcacacactactttgagcactcaaggcaagcttta-3’ | MWG Eurofins (HPLC) |
Tabela 2: Tabela 17 syntetycznych oligonukleotydów kontrolnych użytych jako próbka kontrolna dodatnia. Pierwsze 13 oligonukleotydów wybrano na podstawie wielkości [długie (116 do 120 nt), średnie (69 do 85 nt), krótkie (33 do 41 nt)] oraz ich końców 3'. Tabela została opracowana na podstawie wcześniejszej publikacji6. Kliknij tutaj, aby pobrać tę tabelę jako plik excel.