Artykuł metodologiczny

Badanie epitranskryptomów m6A i m5C po infekcji wirusowej: przykład z HIV

4.3K wyświetleń

⸱

DOI:

10.3791/62426

⸱

5 marca 2022

 ,  , 
1Institute of Microbiology, Lausanne University Hospital and University of Lausanne, 2Bioinformatics Core Facility, SIB Swiss Institute of Bioinformatics

W tym artykule

Podsumowanie

Rola modyfikacji RNA w infekcjach wirusowych dopiero zaczyna być badana i może uwypuklić nowe mechanizmy interakcji wirus-gospodarz. W tej pracy przedstawiamy proces badania modyfikacji RNA m6A i m5C w kontekście infekcji wirusowych.

Streszczenie

Rola modyfikacji RNA w procesach biologicznych była przedmiotem coraz większej liczby badań w ciągu ostatnich kilku lat i jest obecnie znana jako epitranskryptomika. Na cząsteczkach mRNA opisano m.in. modyfikacje RNA N6-metyloadenozyny (m6A) i 5-metylocytozyny (m5C), które mogą odgrywać rolę w modulowaniu procesów komórkowych. Epitranskryptomika jest zatem nową warstwą regulacji, którą należy wziąć pod uwagę oprócz analiz transkryptomicznych, ponieważ może być również zmieniana lub modulowana przez ekspozycję na dowolny czynnik chemiczny lub biologiczny, w tym infekcje wirusowe.

Tutaj prezentujemy przepływ pracy, który pozwala na analizę wspólnego komórkowego i wirusowego krajobrazu epitranskryptomicznego znaczników m6A i m5C jednocześnie, w komórkach zakażonych lub nie ludzkim wirusem niedoboru odporności (HIV). Po izolacji i fragmentacji mRNA z komórek zakażonych wirusem HIV i niezakażonych, zastosowaliśmy dwie różne procedury: MeRIP-Seq, technikę opartą na immunoprecypitacji RNA, w celu wzbogacenia o fragmenty RNA zawierające znacznik m6A oraz BS-Seq, technikę opartą na konwersji wodorosiarczynu, w celu identyfikacji znaku m5C w rozdzielczości pojedynczego nukleotydu. Po wychwyceniu specyficznym dla metylacji biblioteki RNA są przygotowywane do sekwencjonowania o wysokiej przepustowości. Opracowaliśmy również dedykowany potok bioinformatyczny w celu identyfikacji transkryptów metylowanych różnicowo (DM) niezależnie od ich podstawowego profilu ekspresji.

Ogólnie rzecz biorąc, metodologia pozwala na jednoczesną eksplorację wielu znaczników epitranskryptomicznych i dostarcza atlas transkryptów DM w przypadku infekcji wirusowej lub innych zaburzeń komórkowych. Podejście to stwarza nowe możliwości identyfikacji nowych graczy i nowych mechanizmów odpowiedzi komórek, takich jak czynniki komórkowe promujące lub ograniczające replikację wirusa.

Wprowadzenie

Od dawna wiadomo, że cząsteczki RNA mogą być modyfikowane, a do tej pory opisano ponad 150 modyfikacji potranskrypcyjnych1. Polegają one na dodawaniu grup chemicznych, głównie grup metylowych, do praktycznie dowolnej pozycji pierścieni pirymidynowych i purynowych cząsteczek RNA2. Wykazano już, że takie modyfikacje potranskrypcyjne są wysoce wzbogacone w transferowy RNA (tRNA) i rybosomalny RNA (rRNA), a ostatnio opisano je również w cząsteczkach mRNA.

Rozwój nowych technologii, takich jak sekwencjonowanie nowej generacji (NGS), oraz produkcja specyficznych przeciwciał rozpoznających określone modyfikacje chemiczne pozwoliły, po raz pierwszy, na zbadanie lokalizacji i częstotliwości specyficznych modyfikacji chemicznych na poziomie całego transkryptomu. Postępy te doprowadziły do lepszego zrozumienia modyfikacji RNA i mapowania kilku modyfikacji cząsteczek mRNA3,4.

Podczas gdy epigenetyka bada rolę modyfikacji DNA i histonów w regulacji transkryptomu, epitranskryptomika w podobny sposób skupia się na modyfikacjach RNA i ich roli. Badanie modyfikacji epitranskryptomicznych stwarza nowe możliwości zwrócenia uwagi na nowe mechanizmy regulacji, które mogą dostrajać różne procesy komórkowe (tj. splicing RNA, eksport, stabilność i translację)5. Nie było więc wielkim zaskoczeniem, że ostatnie badania ujawniły wiele modyfikacji epitranskryptomicznych po infekcji wirusowej zarówno w komórkowych, jak i wirusowych RNAs6. Wirusy przebadane do tej pory obejmują zarówno wirusy DNA, jak i RNA; wśród nich HIV można uznać za pionierski przykład. Podsumowując, odkrycie metylacji RNA w kontekście infekcji wirusowych może pozwolić na zbadanie jeszcze nieopisanych mechanizmów ekspresji lub replikacji wirusa, dostarczając w ten sposób nowych narzędzi i celów do ich kontrolowania7.

W dziedzinie epitranskryptomiki HIV, modyfikacje transkryptów wirusowych były szeroko badane i wykazały, że obecność tej modyfikacji była korzystna dla replikacji wirusa8,9,10,11,12,13. Do tej pory można stosować różne techniki wykrywania znaczników epitranskryptomicznych na poziomie całego transkryptomu. Najczęściej stosowane techniki identyfikacji m6A opierają się na technikach wytrącania immunologicznego, takich jak MeRIP-Seq i miCLIP. Podczas gdy MeRIP-Seq opiera się na fragmentacji RNA w celu wychwycenia fragmentów zawierających metylowane reszty, miCLIP opiera się na generowaniu mutacji sygnatury specyficznych dla przeciwciała α-m6A po sieciowaniu UV RNA-przeciwciało, umożliwiając w ten sposób bardziej precyzyjne mapowanie.

Wykrycie modyfikacji m5C może być osiągnięte albo za pomocą technologii opartych na przeciwciałach, podobnych do detekcji m6A (m5C RIP), albo przez konwersję bisulfitu, albo przez AZA-IP lub przez miCLIP. Zarówno Aza-IP, jak i m5C miCLIP wykorzystują specyficzną metylotransferazę jako przynętę do celowania w RNA podczas przechodzenia metylacji RNA. W Aza-IP komórki docelowe są wystawione na działanie 5-azacytydyny, co powoduje losowe wprowadzenie miejsc analogu 5-azacytydyny cytydyny do powstającego RNA. W miCLIP metylotransferaza NSun2 jest genetycznie modyfikowana tak, aby zawierała mutację C271A14,15.

W tej pracy skupiamy się na podwójnej charakterystyce modyfikacji m6A i m5C w zakażonych komórkach, używając HIV jako modelu. Po optymalizacji metodologicznej opracowaliśmy przepływ pracy, który łączy immunoprecypitację metylowanego RNA (MeRIP) i konwersję wodorosiarczynu RNA (BS), umożliwiając jednoczesną eksplorację znaczników epitranskryptomicznych m6A i m5C na poziomie całego transkryptomu, zarówno w kontekście komórkowym, jak i wirusowym. Ten przepływ pracy można zaimplementować na komórkowych ekstraktach RNA, a także na RNA wyizolowanym z cząstek wirusa.

Podejście do immunoprecypitacji metylowanego RNA (MeRIP)16 pozwalające na badanie m6A na poziomie transkryptomu jest dobrze ugruntowane, a tablica przeciwciał specyficznych dla m6A jest dostępna na rynku17. Metoda ta polega na selektywnym wychwytywaniu fragmentów RNA zawierających m6A przy użyciu przeciwciała swoistego dla m6A. Dwie główne wady tej techniki to (i) ograniczona rozdzielczość, która w dużym stopniu zależy od wielkości fragmentów RNA, a tym samym zapewnia przybliżoną lokalizację i region zawierający metylowaną resztę, oraz (ii) dużą ilość materiału potrzebnego do przeprowadzenia analizy. W poniższym zoptymalizowanym protokole ustandaryzowaliśmy rozmiar fragmentu do około 150 nt i zmniejszyliśmy ilość materiału wyjściowego z 10 μg RNA wyselekcjonowanego przez poli-A, co jest obecnie zalecaną ilością materiału wyjściowego, do zaledwie 1 μg RNA wyselekcjonowanego przez poli-A. Zmaksymalizowaliśmy również wydajność odzyskiwania fragmentów RNA m6A związanych ze specyficznymi przeciwciałami za pomocą elucji metodą konkurencji z peptydem m6A zamiast bardziej konwencjonalnych i mniej specyficznych metod elucji przy użyciu technik opartych na fenolu lub proteinazie K. Głównym ograniczeniem tego testu opartego na RIP pozostaje jednak nieoptymalna rozdzielczość, która nie pozwala na identyfikację precyzyjnego zmodyfikowanego nukleotydu A.

Analiza znaku C m5C może być obecnie przeprowadzona przy użyciu dwóch różnych podejść: metody opartej na RIP z przeciwciałami specyficznymi dla m5C i konwersją wodorosiarczynu RNA. Ponieważ RIP oferuje tylko ograniczoną rozdzielczość identyfikacji metylowanej pozostałości, zastosowaliśmy konwersję wodorosiarczynu, która może zapewnić rozdzielczość pojedynczego nukleotydu. Narażenie RNA na wodorosiarczyn (BS) prowadzi do deaminacji cytozyny, przekształcając w ten sposób resztę cytozyny w uracyl. Tak więc podczas reakcji konwersji wodorosiarczynu RNA każda niemetylowana cytozyna jest deaminowana i przekształcana w uracyl, podczas gdy obecność grupy metylowej w pozycji 5 cytozyny ma działanie ochronne, zapobiegając deaminacji wywołanej przez BS i zachowując resztę cytozyny. Podejście oparte na BS pozwala na wykrycie nukleotydu zmodyfikowanego m5C w rozdzielczości pojedynczej zasady oraz na ocenę częstotliwości metylacji każdego transkryptu, zapewniając wgląd w dynamikę modyfikacji m5C18. Główne ograniczenie tej techniki polega jednak na fałszywie dodatnim odsetku metylowanych pozostałości. Rzeczywiście, konwersja BS jest skuteczna w przypadku jednoniciowego RNA z dostępnymi resztami C. Jednak obecność ścisłej struktury drugorzędowej RNA może maskować pozycję N5C i utrudniać konwersję BS, co skutkuje niemetylowanymi resztami C, które nie są przekształcane w reszty U, a tym samym fałszywie dodatnimi. Aby obejść ten problem i zminimalizować odsetek wyników fałszywie dodatnich, zastosowaliśmy 3 rundy cykli denaturacji i konwersji wodorosiarczynu19. Wprowadziliśmy również 2 kontrole w próbkach, aby umożliwić oszacowanie wydajności konwersji wodorosiarczynu: zwiększyliśmy kontrolę sekwencjonowania ERCC (niemetylowane standaryzowane i dostępne na rynku sekwencje)20, a także RNA zubożone w poli-A, aby ocenić współczynnik konwersji wodorosiarczynu z jednej strony, a z drugiej strony zweryfikować za pomocą RT-PCR obecność znanego i dobrze zachowanego miejsca metylacji, C4447, z drugiej strony na rybosomalnym RNA 28S21.

W dziedzinie wirusologii, połączenie tych dwóch metod badań epitranskryptomicznych z sekwencjonowaniem nowej generacji i dokładną analizą bioinformatyczną pozwala na dogłębne badanie dynamiki m6A i m5C (tj. modyfikacji RNA, zmian czasowych, które mogą wystąpić po infekcji wirusowej i mogą odkryć szereg nowych terapeutycznie istotnych celów do użytku klinicznego).

Protokół

1. Przygotowanie komórek

UWAGA: W zależności od typu komórki i jej zawartości RNA, początkowa liczba komórek może się różnić.

  1. Mieć wystarczającą ilość komórek, aby uzyskać od 200 do 500 μg całkowitego RNA lub 5-7 μg RNA wyselekcjonowanego przez poli-A. Na przykład komórki 50 x 106 SupT1 powinny dać około 500 μg całkowitego RNA po ekstrakcji odczynnikami na bazie fenolu, a zatem jest to wymagane dla każdego badanego stanu.
  2. Przygotuj wymaganą liczbę komórek zgodnie z planem eksperymentalnym, a tym samym zgodnie z liczbą badanych warunków (infekcja, punkty czasowe, leczenie). Jeśli eksperyment ma na celu uzyskanie komórek niezakażonych i komórek zakażonych wirusem HIV w ciągu 24 godzin po zakażeniu, potrzeba w sumie 100 x 106 komórek, połowa w stanie niezakażonym i połowa w stanie zakażonym.

2. Ekstrakcja RNA

  1. Z komórek: ekstrakcja RNA za pomocą fenolu-chloroformu
    1. Dla każdego warunku zebrać komórki (np. 50 x 106) przez odwirowanie i odrzucić supernatant.
    2. Dodać 5 ml odczynnika na bazie fenolu do każdej6-komórkowej osadki o wymiarach 50 x 10 i wymieszać, kilkakrotnie pipetując w górę iw dół.
    3. Inkubować przez 5 minut w temperaturze pokojowej, aby umożliwić całkowitą lizę. Lizowane komórki mogą być przechowywane w temperaturze -80 °C lub przetwarzane bezpośrednio.
      UWAGA: W razie potrzeby komórki można również podzielić na porcje 10 x 106 komórek na probówkę w probówkach o pojemności 1,5 ml i poddać lizie w 1 ml odczynnika na bazie fenolu w celu wygodniejszego przechowywania.
    4. Dodać 1 ml chloroformu i wymieszać przez odwrócenie.
    5. Inkubować przez 3 minuty w temperaturze pokojowej.
    6. Wirować przez 15 minut przy 2.000 x g i 4 °C.
    7. Odpipetować fazę wodną (faza górna) i przenieść do nowej probówki. Zakończ przenoszenie fazy wodnej, ustawiając probówkę pod kątem 45° i ostrożnie odpipetując roztwór.
      UWAGA: Ilość fazy wodnej może się różnić w zależności od próbki, ale powinna być zbliżona do ilości chloroformu dodanego do próbki (tj. 1 ml). Nie przenoś żadnej warstwy interfazowej ani organicznej! Zastosowanie lamp z blokadą fazową lub roztrząsaczem faz może ułatwić ten proces.
    8. Dodać 0,5 ml izopropanolu o jakości cząsteczkowej 100% do fazy wodnej.
    9. Inkubować przez 1 godzinę w temperaturze -80 °C, aby umożliwić wytrącenie RNA.
    10. Wirować przez 10 minut w temperaturze 12 000 x g i temperaturze 4 °C w celu osadzenia wytrąconego RNA.
    11. Wyrzucić supernatant i ponownie zawiesić osad RNA w 1 ml 75% etanolu klasy biologii molekularnej. Wiruj krótko.
    12. Wirować przez 5 minut w temperaturze 7 500 x g i temperaturze 4 °C i odrzucić supernatant.
    13. Suszyć pelety na powietrzu przez 15 minut.
    14. Zawiesić osad w 20 μl wody wolnej od RNaz i przenieść do nowej probówki.
    15. Umyj pustą probówkę dodatkowymi 20 μl wody, aby zmaksymalizować odzysk RNA, i połącz z pierwszymi 20 μl objętości.
    16. Określ ilościowo całkowite RNA za pomocą spektrofotometru i oceń jakość RNA za pomocą analizatora fragmentów.
  2. Z cząstek wirusowych: ekstrakcja RNA za pomocą zestawu do ekstrakcji wirusowego RNA opartego na kolumnie
    UWAGA: Ekstrakcja RNA z cząstek wirusa za pomocą odczynnika na bazie fenolu skutkuje niską jakością wirusowego RNA i bibliotekami o niższej jakości. W związku z tym należy preferować ekstrakcję RNA opartą na kolumnach. Zestawy do ekstrakcji RNA wykorzystujące nośnik RNA do elucji i odzyskiwania RNA nie są odpowiednie do tej procedury i należy ich unikać. Ponieważ RNA HIV jest poli-adenylowane, bezpośrednia ekstrakcja RNA bez dalszej izolacji mRNA jest wystarczająca, aby wejść do rurociągów MeRIP-Seq i BS-Seq. Zwykle 1-2 ml supernatantu wirusowego z powszechnie zakażonych komórek powinno zapewnić wystarczającą ilość RNA do wykonania całego przepływu pracy.
    1. Przygotować bufor, dodając 150 μl beta-merkaptoetanolu do 30 ml buforu do lizy. Rozpuść Viral Wash Buffer, dodając 96 ml 100% etanolu.
    2. Zbierz supernatanty zawierające wirusy i odwiruj w celu osadzenia szczątków komórek, aby zminimalizować zanieczyszczenie komórkowym RNA.
    3. Przenieść 1 ml supernatantu wirusa do probówki o pojemności 15 ml.
    4. Dodaj 3 ml wirusowego buforu RNA do 1 ml próbki wirusa i wymieszaj przez wirowanie.
    5. Przenieść 700 μl próbki do kolumny, umieścić w probówce zbiorczej.
    6. Wirować przez 2 minuty przy 13 000 x g w temperaturze pokojowej.
    7. Odrzuć przepływ.
    8. Powtarzaj 3 poprzednie kroki, aż cała próbka zostanie przetworzona, a tym samym całe RNA zostanie wychwycone na kolumnie matrycowej na bazie krzemionki.
    9. Dodać 500 μl buforu do płukania wirusów do kolumny.
    10. Wirować przez 1 minutę przy 10 000 x g w temperaturze pokojowej. Odrzuć przepływ.
    11. Dodać 200 μl buforu do płukania wirusów do kolumny.
    12. Wirować przez 1 minutę przy 10 000 x g w temperaturze pokojowej. Odrzuć przepływ.
    13. Umieść kolumnę w pustej probówce zbiorczej.
    14. Wirować przez 1 minutę przy 10 000 x g w temperaturze pokojowej, aby dodatkowo usunąć wszelkie pozostałe zanieczyszczenia buforem do płukania.
    15. Ostrożnie przenieść kolumnę do probówki o pojemności 1,5 ml.
    16. Dodaj 20 μl wody wolnej od DNaz/RNaz bezpośrednio do środka matrycy kolumny i odwiruj przy 10 000 x g przez 30 s w temperaturze pokojowej.
    17. Dodaj dodatkowe 10 μl wody wolnej od DNaz/RNaz bezpośrednio do środka matrycy kolumny i ponownie odwiruj przez 30 sekund.
    18. Określ ilościowo całkowite RNA za pomocą spektrofotometru i oceń jakość RNA za pomocą analizatora fragmentów.
      UWAGA: Ekstrakcję RNA można przeprowadzić dowolną metodą, jeśli jakość pobranego RNA jest wysoka, z numerem integralności/jakości RNA > 9. Całkowite RNA może być przechowywane w temperaturze -80 °C do czasu dalszego przetwarzania.

3. Izolacja mRNA metodą selekcji poli-A za pomocą Oligo(dT)25

UWAGA: Ze względu na obecność wysoce metylowanego rybosomalnego RNA w ekstraktach komórkowych, zdecydowanie zaleca się izolowanie poli-A RNA albo przez zubożenie rRNA, albo preferencyjnie przez selekcję dodatnią poli-A. Ten krok jest opcjonalny i powinien być wykonywany tylko dla próbek komórkowego RNA, aby uzyskać wyniki sekwencjonowania w wyższej rozdzielczości. Jeśli analizujesz metylację niepoli-adenylowanych wirusowych RNA, faworyzuj raczej zubożenie rRNA niż selekcję poli-A lub ostatecznie przeprowadź analizę całkowitego RNA.

  1. Przygotowanie ściegu do wychwytywania poly-A
    1. Zawiesić fiolkę z kulkami magnetycznymi Oligo(dT)25, wirując przez >30 s.
    2. Przenieś 200 μl kulek magnetycznych do probówki o pojemności 1,5 ml. Przygotuj liczbę probówek z kulkami magnetycznymi zgodnie z całkowitą ilością próbek RNA do przetworzenia.
      UWAGA: Jedna probówka z 200 μl roztworu podstawowego Dynabead odpowiada 1 mg kulek i może pomieścić próbkę 75 μg całkowitego RNA.
    3. Umieść probówki na magnesie na 1 minutę i wyrzuć supernatant. Wyjmij rurki z magnesu.
    4. Dodać 1 ml buforu wiążącego (20 mM Tris-HCl, pH 7,5, 1,0 M LiCl, 2 mM EDTA) i zawiesić ponownie przez wirowanie. Umieść probówki na magnesie na 1 minutę i wyrzuć supernatant. Wyjmij rurki z magnesu. Powtórzyć.
    5. Ponownie zawiesić umyte kulki magnetyczne w 100 μl buforu wiążącego.
  2. Przygotowanie całkowitego RNA
    1. Rozcieńczyć całkowite RNA w końcowym stężeniu 0,75 μg/μl wodą wolną od RNaz, co odpowiada 75 μg/100 μL.
      UWAGA: Jeśli RNA ma niższe stężenie, postępuj zgodnie z opisem poniżej bez modyfikowania objętości.
    2. Podwielokrotnić całkowitą ilość RNA w wielu probówkach, dozując 100 μl próbki RNA na probówkę.
    3. Dodać 100 μl buforu wiążącego do każdej próbki RNA.
    4. Podgrzewać całkowite RNA do temperatury 65 °C przez 2 minuty, aby rozbić struktury drugorzędowe.
    5. Umieść natychmiast na lodzie, aż będziesz gotowy do przejścia do następnego kroku.
      UWAGA: Czas inkubacji może się różnić w zależności od liczby próbek do przetworzenia, ale nie powinien przekraczać 1 godziny, aby uniknąć degradacji RNA.
  3. Wybór Poly-A
    1. Do każdej probówki RNA (od kroku 3.2) dodaj 100 μl przemytych kulek magnetycznych (od kroku 3.1).
    2. Dokładnie wymieszać, pipetując w górę i w dół, a następnie pozostawić na obracającym się kole w temperaturze pokojowej na 15 minut.
    3. Otwórz wszystkie probówki, umieść je na magnesie na 1 minutę i ostrożnie usuń cały supernatant.
    4. Odzyskać supernatant w nowej probówce i odłożyć na bok do drugiej rundy wychwytywania RNA (krok 3.3.14), w celu poprawy ostatecznego odzysku poli-A.
    5. Wyjmij rurkę z magnesu i dodaj 200 μl buforu myjącego (10 mM Tris-HCl, pH 7,5, 0,15 M LiCl, 1 mM EDTA). Wymieszać ostrożnie pipetując 4 do 5 razy.
    6. Umieść probówkę na magnesie na 1 minutę i wyrzuć supernatant.
    7. Powtórz krok prania raz (powtórz kroki 3.3.5 i 3.3.6).
    8. Dodaj 20 μl lodowatego 10 mM Tris-HCl, aby wymyć poli-A RNA z kulek.
    9. Inkubować w temperaturze 80 °C przez 2 minuty.
    10. Umieść probówkę na magnesie i szybko przenieś supernatant zawierający poli-A RNA do nowej probówki wolnej od RNaz. Umieść rurkę na lodzie.
    11. Powtórzyć etap elucji (kroki 3.3.8-3.3.10) w celu zwiększenia wydajności.
    12. Umyj te same koraliki raz za pomocą 200 μl bufora do mycia. Wymieszać ostrożnie pipetując 4 do 5 razy.
    13. Umieść na magnesie na 1 minutę i wyrzuć bufor do mycia.
    14. Dodać przepływ z etapu 3.3.4 do kulek i powtórzyć procedurę od wiązania do elucji (kroki 3.3.2 do 3.3.10). Na razie trzymaj eluaty RNA w oddzielnych probówkach.
      UWAGA: Opcjonalnie, ponownie zachować supernatant odpowiadający krokowi 3.3.4 w nowej probówce, ponieważ może on być używany jako kontrola. Pod koniec procedury należy oczyścić i zagęścić RNA przez wytrącanie etanolu lub wybraną metodę opartą na kolumnie (tj. oczyszczanie i koncentrator RNA). Próbka ta odpowiada próbce RNA zubożonej w poli-A i może być stosowana jako kontrola konwersji wodorosiarczynu (etap 8.2.2).
    15. Oznaczyć ilościowo wymyty RNA za pomocą spektrofotometru i zachować podwielokrotność 2 μl w celu dalszej oceny jakości RNA za pomocą analizatora fragmentów.
      UWAGA: RNA Poly-A może być przechowywany w temperaturze -80 °C do czasu, gdy będzie potrzebny.

4. Przepływ pracy RNA

  1. Podziel próbki komórkowego poli-A RNA (mRNA) i wirusowego RNA na 2 podwielokrotności, przeznaczone dla odpowiedniego potoku analizy epitranskryptomicznej:
    (i) 5 μg komórkowego mRNA lub 1 μg wirusowego RNA dla MeRIP-Seq i kontroli wejściowej (przejdź do kroków od 5 do 7 i kroku 9).
    (ii) 1 μg komórkowego mRNA lub 500 ng wirusowego RNA dla BS-Seq (przejdź do kroków 8 i 9).

5. Fragmentacja RNA

UWAGA: Fragmentacja RNA jest przeprowadzana za pomocą odczynnika do fragmentacji RNA i jest przeznaczona dla MeRIP-Seq i kontrolnych próbek RNA. Jest to bardzo ważny krok, który wymaga starannej optymalizacji w celu uzyskania fragmentów o wielkości 100-200 nt.

  1. Podziel całkowitą objętość mRNA na 0,2 ml probówki PCR z 18 μl mRNA/probówkę.
    UWAGA: Pracuj szybko. Nie pracuj z więcej niż 8 próbkami jednocześnie, aby uzyskać powtarzalne wyniki. Skalowanie woluminu w górę nie gwarantuje powtarzalnego i jednolitego rozdrobnienia.
  2. Rozgrzej termocykler w temperaturze 70 °C.
  3. Dodać 2 μl odczynnika do rozdrabniania na brzeg każdej probówki do PCR.
  4. Zamknij probówkę i odwiruj w dół (tak, aby odczynnik wszedł w kontakt z RNA w tym samym czasie dla 8 probówek).
  5. Inkubować próbki przez 15 minut w temperaturze 70 °C w nagrzanym termocyklerze.
  6. Jak tylko inkubacja się zakończy, dodaj szybko 2 μl roztworu Stop do każdej probówki.
  7. Zakręć i pozostaw na lodzie, aż będziesz gotowy do przejścia do następnego kroku.
    UWAGA: Czas inkubacji może się różnić w zależności od liczby próbek do przetworzenia, ale nie powinien przekraczać 1 godziny, aby uniknąć degradacji RNA.
  8. Powtórzyć procedurę dla wszystkich próbek (jeżeli jest więcej niż 8 podwielokrotności).
  9. Połącz probówki razem i przystąp do oczyszczania RNA za pomocą zestawu do czyszczenia i koncentratora RNA (krok 6) lub dowolnego niestandardowego zestawu opartego na kolumnie, aby pozbyć się i odzyskać czyste rozdrobnione RNA w wodzie.

6. Oczyszczanie RNA

UWAGA: Ten krok może być przeprowadzony przez wytrącanie etanolu lub za pomocą dowolnej metody oczyszczania i zagęszczania RNA opartej na kolumnach (np. RNA Clean and Concentrator).

  1. Eluować lub ponownie zawiesić oczyszczony RNA w całkowitej objętości 50-75 μl wody wolnej od DNazy/RNaz.
    UWAGA: Jeśli stosowana jest metoda oparta na kolumnach, zdecydowanie zaleca się dwie rundy elucji, aby zapewnić maksymalne odzysk.
  2. Określ ilościowo oczyszczone pofragmentowane mRNA za pomocą spektrofotometru i oceń jakość RNA za pomocą analizatora fragmentów.
  3. Zachowaj 100 ng pofragmentowanego mRNA jako kontrolę wejściową do przygotowania biblioteki i sekwencjonowania (przejdź do kroku 9). Pozostałe pofragmentowane mRNA (minimum 2,5 μg) można wykorzystać do MeRIP (przejdź do kroku 7.2).

7. MeRIP

UWAGA: Minimum 2,5 μg rozdrobnionego mRNA jest wymagane dla każdej immunoprecypitacji (IP), albo przy użyciu specyficznego przeciwciała anty-m6A (warunek testowy), albo przy użyciu przeciwciała anty-IgG (kontrola negatywna).

  1. Przygotowanie kulek magnetycznych do immunoprecypitacji
    1. Dla każdej próbki przygotuj 4 ml 1x buforu IP w nowej stożkowej probówce, rozcieńczając 800 μl buforu mRNA IP 5x (50 mM Tris-HCl pH 7,4, 750 mM NaCl, 0,5% Igepal CA-630 i woda wolna od nukleaz) z 3,2 ml wody wolnej od nukleaz.
      UWAGA: Potrzebne są co najmniej 2 reakcje (jeden test i jedna kontrola IgG).
    2. Umieść rurkę na lodzie.
    3. Oznaczyć odpowiednią liczbę probówek mikrowirówkowych o pojemności 1,5 ml dla liczby pożądanych reakcji IP:
      n probówek (test) na przeciwciało anty-m6A.
      n probówek (kontrola ujemna) dla normalnej mysiej IgG.
    4. Ponownie zawieś kulki magnetyczne (np. Magna ChIP Protein A/G ) przez odwrócenie i wirowanie. Nie powinny być widoczne żadne kępy koralików.
    5. Dla każdej zaplanowanej reakcji przenieść 25 μl kulek magnetycznych do probówki mikrowirówkowej.
    6. Dodać dziesięć razy więcej bufora 1x IP (z kroku 7.1.1) w stosunku do pierwotnej objętości użytych kulek (tj. 250 μl 1x bufora IP na 25 μl kulek magnetycznych).
    7. Wymieszać kulki, delikatnie pipetując kilkakrotnie w górę i w dół, aby całkowicie się zawiesić.
    8. Umieść rurkę na separatorze magnetycznym na 1 minutę.
    9. Usuń i wyrzuć supernatant, uważając, aby nie zassać żadnych kulek magnetycznych. Wyjmij rurkę z magnesu.
    10. Powtórzyć etap prania (kroki od 7.1.6 do 7.1.9).
    11. Ponownie zawiesić koraliki w 100 μl 1x bufora IP na 25 μl pierwotnej objętości kulek magnetycznych.
    12. Dodać 5 μl przeciwciała (1 μg/μl) na 25 μl pierwotnej objętości kulek magnetycznych.
      n probówek (test) z przeciwciałem anty-m6A (klon 17-3-4-1) [1 μg/μL].
      n probówek (kontrola ujemna) z normalną mysią IgG (1 μg/μL).
    13. Inkubować na obracającym się kole przez 30 minut w temperaturze pokojowej, aby umożliwić koniugację przeciwciał z kulkami magnetycznymi.
    14. Umieścić probówkę na separatorze magnetycznym na 1 minutę. Wyrzucić supernatant. Wyjmij probówkę z magnesu i ponownie zawieś mieszaninę przeciwciała i kulki w 100 μl 1x bufor IP.
  2. Immunoprecypitacja RNA (RIP)
    1. Przygotować 500 μl mieszaniny reakcyjnej RIP na każdą próbkę mRNA o wielkości 2,5 μg w następujący sposób: 2,5 μg w 100 μl rozdrobnionego RNA (z kroku 6.12); 295 μl wody wolnej od nukleaz; 5 μl inhibitora RNazy 40 U/μL; i 100 μl buforu 5x IP.
    2. Dodać 500 μl mieszaniny reakcyjnej RIP do każdej mieszaniny przeciwciało-kulka (~100 μl z kroku 7.1.14). Wymieszać delikatnie pipetując kilka razy, aby całkowicie zawiesić kulki. Ułóż na lodzie.
    3. Inkubować wszystkie probówki RIP na obracającym się kole przez 2 godziny w temperaturze 4 °C.
    4. Odwiruj krótko reakcje MeRIP, aby odwirować kropelki cieczy z boków nasadki i probówki. Umieść probówki na separatorze magnetycznym na 1 minutę.
    5. Przenieść supernatant do nowej probówki wirówkowej, uważając, aby nie naruszyć kulek magnetycznych .
      UWAGA: Przepływ może być zachowany jako kontrola w celu sprawdzenia wydajności procesora RIP (przejdź do kroku 7.3.9).
    6. Wyjmij rurki z magnesu. Umyj kulki, dodając 500 μl zimnego bufora 1x IP. Wymieszaj kulki, delikatnie pipetując kilka razy, aby całkowicie zawiesić kulki.
    7. Umieść probówki na separatorze magnetycznym na 1 minutę i wyrzuć supernatant.
    8. Powtórzyć procedurę prania (kroki 7.2.6-7.2.7) dwukrotnie, co daje w sumie 3 prania.
    9. Umieść probówki na lodzie i natychmiast przystąp do elucji.
  3. Elucji
    1. Przygotować roztwór 20mM 6A, rozpuszczając 10 mg soli sodowej N6-metyloadenozyny, 5′-monofosforanu (m6A) w 1,3 ml wody wolnej od nukleaz. Przygotować 150 μl porcji i przechowywać w temperaturze -20 °C.
    2. Dla każdej próbki (badanie i kontrole): Przygotować 225 μl buforu elucyjnego, mieszając następujące składniki: 45 μl 5x buforu IP, 75 μl 20 mM6A, 3,5 μl 40U/μL inhibitora RNazy i 101,5 μl wody wolnej od nukleaz.
    3. Dodać 100 μl buforu elucyjnego (z kroku 7.3.2) do kulek (z kroku 7.2.9). Wymieszać, delikatnie pipetując kilka razy, aby całkowicie zawiesić kulki.
    4. Inkubować wszystkie probówki przez 1 godzinę, wytrząsając w sposób ciągły na bujaku w temperaturze 4 °C.
    5. Krótko odwirować reakcje RIP, aby odwirować kropelki cieczy z boków nasadki i probówki. Umieść probówki na separatorze magnetycznym na 1 minutę.
    6. Przenieść supernatant zawierający eluowane fragmenty RNA do nowej probówki do mikrowirówki o pojemności 1,5 ml. Uważaj, aby nie zasysać koralików, ponieważ zwiększy to hałas w tle.
    7. Powtórzyć etapy elucji (7.3.3-7.3.6), ponownie dodając 100 μl buforu elucyjnego, inkubując przez 1 godzinę w temperaturze 4 °C i zbierając eluat po rozdzieleniu magnetycznym.
    8. Połączyć wszystkie eluaty z tej samej próbki (całkowita objętość elucji powinna wynosić 200 μl).
    9. Oczyścić wymyty RNA i jego przepływ (opcjonalnie, od kroku 7.2.5) przez wytrącanie etanolu lub wybraną metodę opartą na kolumnie (tj. RNA Clean and Concentrator).
    10. Oceń ilość i jakość RNA w próbkach przepływowych i eluowanych za pomocą analizatora fragmentów przy użyciu zestawu do wykrywania o wysokiej czułości. Jeśli jakość RNA jest zadowalająca, przystąp do przygotowania biblioteki i sekwencjonowania o wysokiej przepustowości (krok 9).
      UWAGA: Ilość RNA pobranego przez MeRIP jest bardzo niska i bezwzględnie wymaga zestawów do wykrywania o wysokiej czułości, aby zapewnić kwantyfikację. Jeśli żaden bioanalizator nie jest dostępny, można przejść na ślepo do przygotowania biblioteki.

8. Konwersja wodorosiarczynu RNA

  1. Kontrola i przygotowanie odczynnika
    1. Kontrola skoków mieszanki ERCC: Dodaj mieszankę ERCC zgodnie z instrukcjami producenta, który zaleca dodanie 0,5 μl nierozcieńczonej mieszanki ERCC do 500 ng mRNA. Ta kontrola może pomóc w ocenie wydajności konwersji wodorosiarczynu.
    2. Spike Poly-A-zubożone RNA (z kroku 3.3.14) w stosunku 1/1000 (tj. 500 pg poli-A-zubożonego RNA na 500 ng mRNA). Próbka ta jest wzbogacona w rybosomalny RNA i dlatego powinna zawierać 28S rRNA, pozytywną kontrolę konwersji wodorosiarczynu.
      UWAGA: Całkowite RNA może być również używane jako kontrola pozytywna zamiast RNA zubożonego w poli-A.
    3. Przeprowadź konwersję wodorosiarczynu za pomocą zestawu do metylacji RNA (np. Zymo EZ).
    4. Bufor do płukania RNA: Przed użyciem dodaj 48 ml 100% etanolu (lub 52 ml 95% etanolu) do 12 ml koncentratu buforu do płukania RNA.
  2. Konwersja wodorosiarczynu
    UWAGA: Konwersję wodorosiarczynu przeprowadzono za pomocą dostępnego na rynku zestawu do konwersji wodorosiarczynu RNA zgodnie z procedurą producenta opisaną poniżej.
    1. W probówkach PCR o pojemności 0,2 ml dodaj 1000 ng mRNA (lub od 300 do 1000 ng). Dodaj kontrole skokowe: 1 μl mieszanki ERCC (krok 8.1.1) i 1000 pg RNA zubożonego w poli-A (krok 8.1.2). Pełna objętość do 20 μl z wodą wolną od DNaz/RNaz.
    2. Dodać 130 μl odczynnika do konwersji RNA do każdej próbki 20 μl RNA.
    3. Wymieszać próbkę, pipetując w górę i w dół.
    4. Odkręć na chwilę, aby upewnić się, że w nasadce lub bokach rurki nie ma kropel.
    5. Umieścić probówki PCR w termocyklerze i wykonać następujące czynności: denaturacja w temperaturze 70 °C przez 5 minut; konwersja w temperaturze 54 °C przez 45 minut; powtórzyć etapy denaturacji i konwersji, w sumie przez 3 cykle; a następnie utrzymać w temperaturze 4 °C przez czas nieokreślony.
      UWAGA: Trzy cykle denaturacji i konwersji wodorosiarczynu zapewniają całkowitą konwersję próbki w kierunku wodorosiarczynu. Próbki mogą być przechowywane w temperaturze -80 °C lub przetwarzane bezpośrednio.
    6. Kontynuować odsiarczanie w kolumnie. Umieść kolumnę w pustej probówce zbiorczej i dodaj do niej 250 μl buforu wiążącego RNA.
    7. Załaduj próbkę (~150 μl z kroku 8.2.5) do kolumny zawierającej bufor wiążący RNA i wymieszaj, pipetując w górę iw dół.
    8. Dodać 400 μl 95-100% etanolu do mieszaniny buforu wiążącego próbkę i RNA w kolumnie. Zamknij nakrętkę i natychmiast wymieszaj, kilkakrotnie odwracając kolumnę.
    9. Wirować z pełną prędkością (≥ 10 000 x g) przez 30 s. Wyrzucić przepływ.
    10. Dodać 200 μl buforu płuczącego RNA do kolumny i wirować z pełną prędkością przez 30 s.
    11. Dodać 200 μl buforu odsulfonacyjnego RNA do kolumny i inkubować w temperaturze pokojowej przez 30 minut. Po inkubacji wirować na pełnych obrotach przez 30 s. Odrzuć przepływ.
    12. Dodać 400 μl buforu płuczącego RNA do kolumny i wirować z pełną prędkością przez 30 s. Powtórzyć etap płukania z dodatkowymi 400 μl buforu płuczącego RNA. Odrzuć przepływ.
    13. Odwirować kolumnę w opróżnionej probówce zbiorczej z pełną prędkością przez 2 minuty. Przenieść kolumnę do probówki wolnej od RNaz.
    14. Dodać ≥ 10 μl wody wolnej od DNaz/RNaz bezpośrednio do matrycy kolumny i inkubować przez 1 minutę w temperaturze pokojowej. Wirować z pełną prędkością przez 30 s.
      UWAGA: Zazwyczaj eluujemy w objętości 20 μL. Uwolniony RNA może być zużyty natychmiast lub przechowywany w temperaturze -20 °C przez okres do 3 miesięcy. W przypadku długotrwałego przechowywania przechowywać w temperaturze -80 °C.
    15. Pobrać 2,5 μl do oceny jakości i ilości RNA za pomocą analizatora fragmentów, a następnie przejść do przygotowania biblioteki i sekwencjonowania o wysokiej przepustowości (krok 9).
    16. Weź 4 μl przekształconego RNA w celu kontroli wydajności konwersji wodorosiarczynu (krok 8.3).
  3. Kontrola konwersji wodorosiarczynu przez RT-PCR
    UWAGA: Ten krok zapewnia, że konwersja wodorosiarczynu zakończyła się pomyślnie przed przystąpieniem do sekwencjonowania. Rybosomalny RNA 28S z Homo sapiens zostanie wykorzystany jako kontrola pozytywna do analizy metylacji RNA, ponieważ reszta C w pozycji 4447 (akces GenBank # NR_003287) została opisana jako w 100% metylowana.
    Sekwencje podkładowe:
    Podkład H 28SF: 5'-GGGGTTTTAYGATTTTTTTTTTTTTGGG-3'
    Podkład H 28SR: 5'-CCAACTCACRTTCCCTATTAATAAATAAAC-3'
    1. Przygotuj mieszaninę reakcyjną odwrotnej transkrypcji (RT) przy użyciu zestawu do odwrotnej transkrypcji cDNA o dużej pojemności. Rozmrozić elementy zestawu na lodzie i przygotować główną mieszankę RT na lodzie w następujący sposób:
      4 μl RNA przekształconego w wodorosiarczyn (z kroku 8.2.14):
      2 μl buforu 10xRT
      0,8 μL 25x dNTP Mix [100 mM]
      2 μL 10 losowych starterów RT
      1 μl odwrotnej transkryptazy MultiScribe
      1 μl inhibitora RNazy
      9,2 μl H2O
      wolne od nukleaz UWAGA: Każda reakcja RT powinna zawierać 20 μl końcowej objętości w 0,2 ml probówkach PCR.
    2. Umieścić probówki w termocyklerze z następującym programem RT: 25 °C przez 10 min; 37 °C przez 120 min; 85 °C przez 5 min; następnie w 4 °C na czas nieokreślony.
    3. Przygotuj reakcję PCR, aby amplifikować specyficznie 28S rRNA za pomocą enzymu korygującego PCR. Rozmrozić elementy zestawu na lodzie, delikatnie odwirować i krótko odwirować. Przygotować mieszankę wzorcową PCR na lodzie lub na lodowatym metalowym uchwycie na płytkę w następujący sposób:
      0,6 μl 10 μM H 28SF startera
      0,6 μl 10 μM H 28SF startera
      6,5 μl matrycy cDNA
      22,5 μl głównej mieszanki polimerazy DNA
      UWAGA: Każda reakcja PCR powinna zawierać 20 μl końcowej objętości w 0,2 ml probówek PCR.
    4. Umieścić probówki w termocyklerze z następującym programem PCR: początkowa denaturacja w temperaturze 95 °C przez 5 minut; 45 cykli denaturacji (95 °C przez 15 s), wyżarzanie (57 °C przez 30 s) i wydłużenie (72 °C przez 15 s), końcowe wydłużenie w temperaturze 72°C przez 10 minut, a następnie utrzymywanie w temperaturze 4 °C przez czas nieokreślony.
    5. Przeprowadź 10 μl reakcji na 2% żelu agarozowym. Przewidywana wielkość pasma to 130 - 200 pz.
  4. Sekwencjonowanie produktów PCR
    1. Oczyść produkty PCR wybraną metodą kolumnową w celu usunięcia enzymów i reszt dNTP oraz eluuj amplifikowane DNA w co najmniej 20 μl wody wolnej od DNazy/RNaz.
    2. Oznacz ilościowo oczyszczone DNA za pomocą spektrofotometru.
    3. Reakcja sekwencjonowania
      1. Użyć 40 ng produktu PCR/reakcji sekwencjonowania.
      2. Sekwencjonuj w obu kierunkach za pomocą spłonek H 28SF i H 28SR.
      3. Dopasuj sekwencje do znanej nieprzekształconej sekwencji (28S rybosomalny N5 (RNA28SN5). Sprawdzić, czy w miejscu C4447 nie ma pozostałości C, a w innych miejscach nie ma pozostałości T.

9. Przygotowanie biblioteki i sekwencjonowanie o wysokiej przepustowości

  1. Przygotuj biblioteki do sekwencjonowania za pomocą zestawów mRNA (np. Illumina TruSeq Stranded), uruchamiając protokół na etapie Elute-Prime-Fragment i postępując zgodnie z instrukcjami producenta.
    1. Jednak w przypadku wejściowych próbek RNA-Seq i MeRIP-Seq należy inkubować próbki w temperaturze 80 °C przez 2 minuty, aby tylko je zagruntować, ale nie dalej rozdrobnić.
  2. Przeprowadź sekwencjonowanie za pomocą platform Illumina. Reakcje sekwencjonowania można przeprowadzać zgodnie z preferencjami i projektem eksperymentalnym, pojedynczymi lub sparowanymi końcami, o długości co najmniej 100 nt.

10. Analizy bioinformatyczne

  1. m6A Przetwarzanie danych
    1. Uruchom FASTQC24, aby ocenić jakość odczytu w m6A i wprowadzić pliki FASTQ z sekwencjonowania.
    2. Uruchom Atropos25, aby przyciąć sekwencje końcowe i adapterowe niskiej jakości z odczytów. Ustaw następujące parametry podczas uruchamiania Atropos.
      1. Usuń następujące sekwencje adapterów: AGATCGGAAGAG, CTCTTCCGATCT, AACACTCTTTCCCT, AGATCGGAAGAGCG, AGGGAAAGAGTGTT, CGCTCTCCCCGATCT.
      2. Użyj następującego odcięcia jakości Phred: 5, do przycinania końcówek niskiej jakości, zgodnie z zaleceniami producenta (https://support.illumina.com/downloads/illumina-adapter-sequences-document-1000000002694.html).
      3. Użyj następującej minimalnej długości odczytu po przycięciu: 25 par zasad.
    3. Scal ludzki genom GRh38 i odniesienie do HIV [Integrated linear pNL4-3Env-GFP] w formacie FASTA.
    4. Indeksuj scalone odniesienie za pomocą HISAT226.
    5. Uruchom HISAT2 na przyciętych odczytach, aby wyrównać do indeksowanego odwołania. Użyj domyślnych parametrów HISAT.
    6. Sortuj i indeksuj wyrównane odczyty za pomocą SAMtools27.
    7. Uruchom SAMtools stat i Qualimap 228, aby sprawdzić jakość sekwencjonowanych bibliotek po wyrównaniu.
    8. Opcjonalnie zbierz i podsumuj miary jakości z poprzedniego kroku za pomocą multiQC29.
    9. Genom HIV ma homologiczne sekwencje 634 pz w 5' LTR i 3' LTR: Dopasuj odczyty wielomapowania z 5' LTR do odpowiedniego regionu 3' LTR za pomocą SAMtools.
    10. Aby zidentyfikować piki m6A, uruchom oprogramowanie do wywoływania szczytów MACS230 (v 2.1.2). Starannie dobieraj parametry pracy MACS2, aby zapewnić prawidłowe działanie danych RNA-Seq, ponieważ poziom ekspresji genów może mieć wpływ na wywoływanie pików, a krótkie eksony mogą być błędnie nazywane pikami. W związku z tym sygnał wejściowy musi zostać odjęty od sygnału m6A, bez wygładzania rutynowo stosowanego przez MACS2 do danych opartych na DNA. Zastosuj następujące parametry do podkomendy "callpeak" z MACS2:
      -keep-dup auto (kontroluje zachowanie MACS2 w stosunku do zduplikowanych odczytów, 'auto' pozwala MACS obliczyć maksymalną liczbę odczytów w dokładnie tym samym miejscu na podstawie rozkładu dwumianowego przy użyciu 1e-5 jako odcięcia wartości p)
      -g 2.7e9 (wielkość ludzkiego genomu w BP)
      -q 0,01 (minimalna granica FDR do wywołania znaczących szczytów)
      -nomodel (aby pominąć budowanie modelu przesuwania, który jest dostosowany do eksperymentów ChIP-Seq)
      -slocal 0
      -llocal 0 (ustawienie tego i poprzedniego parametru na 0 pozwala MACS2 na bezpośrednie odejmowanie, bez wygładzania, odczyty wejściowe z m6A odczytuje)
      -extsize 100 (średnia długość fragmentów w bp)
      -B
    11. Uruchom podkomendę wywołania piku różnicowego MACS2, 'bdgdiff', aby porównać próbki zakażone i niezainfekowane. 'bdgdiff' przyjmuje jako dane wejściowe pliki bedGraph wygenerowane przez 'callpeak' w poprzednim kroku. Dla każdego punktu czasowego uruchom porównanie próbek zakażonych i niezakażonych za pomocą "bdgdiff", odejmując odpowiedni sygnał wejściowy od sygnału m6A i podając dodatkowe parametry: -g 60 -l 120.
  2. m5C Przetwarzanie danych
    1. Uruchom Cutadapt31, aby przyciąć sekwencje adapterów z nieprzetworzonych odczytów, z następującymi parametrami:
      adapter "AGATCGGAAGAGCACACGTCTGAAC"
      -minimalna-długość=25.
    2. Odwrotne uzupełnienie przyciętych odczytów przy użyciu seqkit32, ponieważ protokół sekwencjonowania generuje odczyty z odwrotnej nici.
    3. Uruchom FastQC, aby sprawdzić jakość odczytu.
    4. Scal ludzki genom GRh38 i odniesienie do HIV [Integrated linear pNL4-3Env-GFP] w formacie FISTA.
    5. Indeksuj scalone odniesienie z aplikacją meRanGh z pakietu meRanTK class33.
    6. Wyrównaj z meRanGh z następującymi parametrami:
      -UN umożliwiające zapisywanie niezmapowanych odczytów do plików wyjściowych
      -MM umożliwiające zapis odczytów z wieloma mapami do pliku wyjściowego
      -bg dla danych wyjściowych w bedGraph
      -mbgc 10 Filtruj raportowane regiony według zasięgu (co najmniej 10 odczytów pokrycia)
    7. Genom HIV ma homologiczne sekwencje 634 pz w 5' LTR i 3' LTR: dopasuj odczyty multimapowania z 5' LTR do odpowiedniego regionu 3' LTR za pomocą SAMtools.
    8. Uruchom wywołanie metylacji za pomocą narzędzia meRanCall, dostarczonego przez meRanTK, z następującymi parametrami:
      -rl = 126, długość odczytu
      -ei = 0,1, przedział błędu dla obliczenia wartości p szybkości metylacji
      -cr = 0,99, oczekiwana konwersja
    9. Uruchom estimateSizeFactors.pl narzędziowy MeRanTK, aby oszacować współczynniki wielkości każdej próbki. Współczynniki wielkości zostaną użyte jako parametry w następnym kroku.
    10. Uruchom MeRanCompare w celu przeprowadzenia różnicowej analizy metylacji osób niezakażonych i zakażonych w punktach czasowych 12, 24 i 36 godzin. Stosowane są następujące parametry: wartość istotności 0,01 jako minimalny próg raportowania i współczynniki wielkości z poprzedniego kroku.

Wyniki

Ten schemat postępowania okazał się użyteczny w badaniu roli metylacji m6A i m5C w kontekście zakażenia HIV. W tym celu wykorzystaliśmy model linii komórkowej limfocytów T CD4+ (SupT1), które poddano zakażeniu HIV lub pozostawiono nieleczone. Procedurę rozpoczęliśmy od 50 milionów komórek na każdą kondycję, uzyskując średnio 500 µg całkowitego RNA o wskaźniku jakości RNA wynoszącym 10 (Rycina 1A-B). Po selekcji poly-A odzyskaliśmy od 10 do 12 µg mRNA na kondycję (co stanowi około 2% całkowitego RNA) (Rycina 1B). Na tym etapie wykorzystaliśmy 5 µg RNA wyselekcjonowanego metodą poly-A do procedury MeRIP-Seq oraz 1 µg do procedury BS-Seq. Ponieważ RNA HIV jest poliadenylowane, nie są wymagane żadne dodatkowe działania i procedury MeRIP-Seq oraz BS-Seq mogą być zastosowane bezpośrednio.

Schemat procesu ekstrakcji RNA; kroki selekcji i sekwencjonowania mRNA; wykres danych przedstawia czystość RNA.
Rysunek 1: Przygotowanie RNA do dalszych zastosowań. A) Schemat blokowy przedstawiający przygotowanie i rozdział RNA dla równoległych procesów MeRIP-Seq i BS-Seq. Każdy wypełniony sześciokąt reprezentuje typ modyfikacji RNA, np. m6A (zielony) lub m5C (różowy). Wskazano ilości materiału RNA niezbędne do przeprowadzenia eksperymentu. B) Reprezentatywne wyniki przedstawiające oczekiwane profile rozkładu RNA (rozmiar i ilość) po ekstrakcji całkowitego RNA (górny panel) oraz selekcji poly-A (dolny panel). Próbki wprowadzono do analizatora fragmentów z zestawem o standardowej czułości w celu oceny jakości RNA przed przystąpieniem do specyficznych procedur MeRIP-Seq i BS-Seq. RQN: wskaźnik jakości RNA; nt: nukleotydy. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Procedura MeRIP-Seq to technika oparta na immunoprecypitacji RNA, która umożliwia badanie modyfikacji m6A w cząsteczkach RNA. W tym celu RNA jest najpierw fragmentowane, a następnie inkubowane z przeciwciałami specyficznymi dla m6A sprzężonymi z kulkami magnetycznymi w celu immunoprecypitacji i wychwytu. Fragmenty RNA wzbogacone w MeRIP oraz frakcja niepoddana immunoprecypitacji (input) są następnie sekwencjonowane i porównywane w celu zidentyfikowania regionów RNA zmodyfikowanych przez m6A, a tym samym transkryptów zmetylowanych w pozycji m6A (Rysunek 2A). Rozdzielczość tej techniki zależy od wydajności fragmentacji RNA. W rzeczywistości krótsze fragmenty pozwalają na dokładniejszą lokalizację reszty m6A. W tym przypadku komórkowe RNA wyselekcjonowane pod kątem ogona poly-A oraz RNA wirusowe poddano fragmentacji jonowej z użyciem buforu do fragmentacji RNA przez 15 min w końcowej objętości 20 µL, aby uzyskać fragmenty RNA o długości 100-150 nt. Zaczynając od 5 µg mRNA, odzyskano 4,5 µg sfragmentowanego RNA, co odpowiada wydajności odzysku na poziomie 90% (Rysunek 2B). Jako kontrolę input wykorzystano 100 ng sfragmentowanego, oczyszczonego RNA, który poddano bezpośrednio przygotowaniu biblioteki i sekwencjonowaniu. Pozostałe RNA (~4,4 µg) przetworzono zgodnie z procedurą MeRIP-Seq, która rozpoczyna się od inkubacji sfragmentowanego RNA z kulkami związanymi albo z przeciwciałami specyficznymi przeciwko m6A, albo z przeciwciałami anty-IgG jako kontrolą. Specyficzne RIP dla m6A (MeRIP) z 2,5 µg sfragmentowanego RNA pozwoliło na odzyskanie około 15 ng materiału wzbogaconego w m6A, który poddano przygotowaniu biblioteki i sekwencjonowaniu (Rysunek 2B). RIP z kontrolą anty-IgG, zgodnie z oczekiwaniami, nie dostarczyło wystarczającej ilości RNA, aby umożliwić dalszą analizę (Rysunek 2B).

Analiza mRNA metodą MeRIP-Seq; schemat m⁶A RIP przedstawia metodę sekwencjonowania, analizę danych i wyniki porównawcze.
Rycina 2: Schemat procedury MeRIP-Seq. A) Schematyczny przebieg procesu MeRIP-Seq i kontroli wejściowej (input control). Po selekcji poli-A próbki sfragmentowano na fragmenty o długości 120-150 nt, a następnie poddano je bezpośrednio sekwencjonowaniu (100 ng, kontrola wejściowa) lub wykorzystano do immunoprecypitacji RNA (2,5 µg, RIP) z użyciem przeciwciała specyficznego dla m6A lub przeciwciała anty-IgG jako kontroli negatywnej przed sekwencjonowaniem. B) Reprezentatywne wyniki pokazujące oczekiwane profile rozkładu RNA (rozmiar i ilość) po fragmentacji (panel górny) oraz RIP (panele dolne, MeRIP: lewo, kontrola IgG: prawo). Próbki naniesiono na analizator fragmentów w celu oceny jakości i stężenia RNA przed dalszym przetwarzaniem w ramach przygotowania bibliotek i sekwencjonowania. Analiza sfragmentowanego RNA została przeprowadzona przy użyciu zestawu RNA standard sensitivity, natomiast w przypadku RNA immunoprecypitowanego zastosowano zestaw high sensitivity. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Procedura BS-Seq umożliwia badanie modyfikacji m5C RNA z rozdzielczością pojedynczych nukleotydów i prowadzi do identyfikacji transkryptów zmetylowanych w pozycji m5C. W wyniku konwersji pirosiarczynowej niezmetylowane cytozyny są przekształcane w uracyl, podczas gdy zmetylowane cytozyny pozostają niezmienione (Rycyna 3A). Ze względu na drastyczne warunki procedury konwersji pirosiarczynowej (tzn. wysoką temperaturę i niskie pH), przekonwertowane mRNA ulegają znacznej degradacji (Rycina 3B), jednak nie koliduje to z przygotowaniem biblioteki i sekwencjonowaniem. Konwersja pirosiarczynowa jest wydajna wyłącznie w przypadku jednoniciowego RNA, zatem może być potencjalnie utrudniona przez wtórne dwuniciowe struktury RNA. Aby ocenić wydajność konwersji C-U, wprowadziliśmy dwie kontrole. Jako kontrolę pozytywną wykorzystaliśmy wcześniej opisaną obecność silnie zmetylowanej cytozyny w pozycji C4447 28S rRNA23. Po amplifikacji RT-PCR i sekwencjonowaniu fragmentu o długości 200 bp otaczającego miejsce metylacji zaobserwowaliśmy, że wszystkie cytozyny zostały pomyślnie przekonwertowane na uracyle, pojawiając się zatem jako tymidyny w sekwencji DNA, z wyjątkiem cytozyny w pozycji 4447, która pozostała niezmieniona. Jako kontrolę stopnia konwersji pirosiarczynowej użyliśmy dostępnych komercyjnie syntetycznych sekwencji RNA ERCC. Mieszanina ta składa się z puli znanych, niezmetylowanych i poliadenylowanych sekwencji RNA o zróżnicowanych długościach i strukturach wtórnych. Po przygotowaniu biblioteki i sekwencjonowaniu skupiliśmy się na tych sekwencjach ERCC, aby obliczyć stopień konwersji, co można zrobić poprzez zliczenie liczby przekonwertowanych C spośród wszystkich reszt C we wszystkich sekwencjach ERCC w każdej próbce. Uzyskaliśmy stopień konwersji na poziomie 99,5%, co potwierdza wydajność i powodzenie reakcji konwersji pirosiarczynowej (Rycina 3D).

Proces konwersji pirosiarczynowej, schemat sekwencjonowania mRNA, analiza DNA, wykres wydajności konwersji.
Rysunek 3: Procedura BS-Seq. A) Schematyczny przebieg procesu BS-Seq. Po selekcji poly-A próbki poddawane są działaniu pirosiarczynu, co prowadzi do konwersji C na U (w wyniku deaminacji) w przypadku niezmetylowanych reszt C. Z kolei zmetylowane reszty C (m5C) nie ulegają działaniu pirosiarczynu i pozostają niezmienione. B) Reprezentatywny wynik profilu rozkładu RNA poddanego konwersji pirosiarczynowej (rozmiar i ilość) w analizie na analizatorze fragmentów z zestawem o standardowej czułości. C) Elektroferogram przedstawiający reprezentatywny wynik sekwencjonowania amplikonu RT-PCR regionu otaczającego 100% zmetylowaną cytyzynę C w pozycji 4447 w 28S rRNA (zaznaczoną na niebiesko). W przeciwieństwie do niej, reszty C z sekwencji referencyjnej zostały zidentyfikowane jako reszty T w sekwencji amplikonu ze względu na poprawną konwersję pirosiarczynową. D) Ocena wydajności konwersji C-U poprzez analizę sekwencji spike-in ERCC w komórkach zakażonych i niezakażonych HIV. Średnia wydajność konwersji wynosi 99,5%. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Próbki wzbogacone w M6A, próbki poddane konwersji pirosiarczynowej oraz kontrole wejściowe są następnie przetwarzane w celu przygotowania bibliotek, sekwencjonowania i analizy bioinformatycznej (Rysunek 4). W zależności od projektu eksperymentalnego i stawianych pytań biologicznych można zastosować wiele analiz bioinformatycznych. Jako dowód koncepcji przedstawiamy tutaj reprezentatywne wyniki z jednej potencjalnej aplikacji (t.j. analizy różnicowego metylowania), która koncentruje się na identyfikacji transkryptów o różnym stopniu metylowania indukowanych po infekcji HIV. W skrócie zbadaliśmy poziom metylowania m6A lub m5C transkryptów, niezależnie od ich poziomu ekspresji genów, zarówno w komórkach niezinfekowanych, jak i zakażonych HIV, aby lepiej zrozumieć rolę metylacji RNA w cyklu życiowym wirusa. Po normalizacji ekspresji genów zidentyfikowano, że transkrypt ZNF469 był różnicowo metylowany pod kątem m6A w zależności od stanu zakażenia; w rzeczywistości transkrypt ten nie był metylowany w komórkach niezinfekowanych, podczas gdy po infekcji HIV wykazywał kilka pików metylowania (Rysunek 5A). Podobna analiza różnicowego metylowania m5C wykazała, że transkrypt PHLPP1 zawierał kilka reszt zmetylowanych, które wykazują tendencję do częstszego metylowania w warunkach zakażenia HIV (Rysunek 5B). W tym kontekście obie analizy sugerują, że infekcja HIV wpływa na komórkowy epitranskryptom.

Potoki sekwencjonowania RNA m6A i m5C; schemat blokowy, przetwarzanie plików FASTQ, dopasowanie, wyznaczanie pików.
Rycina 4: Schematyczna reprezentacja bioinformatycznego procesu analizy danych m6A i m5C. Kliknij tutaj, aby zobaczyć powiększoną wersję tej ryciny.

Mapowanie genomiczne ZNF469, PHLPP1; integracja HIV, wykres zmienności sekwencji, analiza chromosomowa.
Rysunek 5: Przykład transkryptów o różnym stopniu metylacji po infekcji. A) Reprezentatywny wynik pokazujący metylację m6A transkryptu ZNF459 w komórkach zainfekowanych HIV (zielone) i niezinfekowanych (szare). Intensywność piku (po odjęciu ekspresji wejściowej) jest przedstawiona na osi y, a pozycja w chromosomie na osi x. Analiza różnicowej metylacji wykazuje, że transkrypt ZFN469 jest hipermetylowany po infekcji HIV. B) Reprezentatywny wynik dla genu z metylacją m5C w komórkach zainfekowanych HIV (górna ścieżka) i niezinfekowanych (dolna ścieżka). Wysokość każdego słupka reprezentuje liczbę odczytów na nukleotyd i pozwala na ocenę pokrycia. Każdy reszt C jest zaznaczony na czerwono, a proporcja zmetylowanego C jest przedstawiona na niebiesko. Dokładny stopień metylacji (%) jest podany nad każdą resztą C. Strzałki wskazują statystycznie istotne różnice w metylacji C. Próbki wizualizowano za pomocą przeglądarki IGV. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Dyskusja

Rola modyfikacji RNA w infekcji wirusowej jest nadal w dużej mierze nieznana. Lepsze zrozumienie roli modyfikacji epitranskryptomicznych w kontekście infekcji wirusowej może przyczynić się do poszukiwania nowych celów leczenia przeciwwirusowego.

W tej pracy zapewniamy kompletny przepływ pracy, który umożliwia badanie epitranskryptomów m6A i m5C zakażonych komórek. W zależności od kwestii biologicznej zalecamy użycie RNA wyselekcjonowanego przez poli-A jako materiału wyjściowego. Chociaż jest to opcjonalne, ponieważ rurociąg może być używany z całkowitym RNA, ważne jest, aby pamiętać, że rRNA, jak również małe RNA, są wysoce zmodyfikowane i zawierają znaczną liczbę metylowanych reszt. Może to skutkować obniżeniem jakości i ilości istotnych danych sekwencjonowania.

Jeśli jednak badanie koncentruje się na niepoli-adenylowanym RNA, etap ekstrakcji RNA powinien być dostosowany, aby uniknąć odrzucenia małego RNA (w przypadku ekstrakcji RNA na podstawie kolumny) i uprzywilejować techniki dezubożenia rybosomów, a nie selekcję poli-A do wejścia do rurociągu.

W celu zapewnienia wysokiej jakości RNA, prawidłowej fragmentacji oraz odpowiedniej jakości RNA wzbogaconego w m6A i przetworzonego BS do przygotowania biblioteki, zdecydowanie zalecamy użycie analizatora fragmentów lub bioanalizatora. Jednak ten sprzęt nie zawsze jest dostępny. Alternatywnie, jakość RNA, mRNA i wielkość rozdrobnionego RNA można również ocenić za pomocą wizualizacji na żelu agarozowym. Alternatywnie, przygotowanie biblioteki można przeprowadzić bez wcześniejszej oceny ilości RNA.

Użyliśmy techniki MeRIP-Seq16 opartej na przeciwciałach, aby zbadać epitranskryptomiczny krajobraz m6A. Technika ta opiera się na immunoprecypitacji RNA i jest skuteczna; Jednak niektóre kroki wymagają starannej optymalizacji i mogą być krytyczne. Chociaż opisano metylacjęm6A zachodzącą głównie w sekwencji konsensusu RRA*CH, motyw ten jest bardzo częsty wzdłuż cząsteczek mRNA i nie pozwala na precyzyjną identyfikację miejsca metylacji. Dlatego tak ważne jest osiągnięcie powtarzalnej i spójnej fragmentacji RNA, generującej małe fragmenty RNA, w celu poprawy rozdzielczości opartej na RIP. W tym protokole zalecamy zoptymalizowaną procedurę, zapewniającą powtarzalne i spójne wyniki w naszych warunkach eksperymentalnych; Jednak ten etap fragmentacji może wymagać dalszej optymalizacji zgodnie z określonymi przykładowymi funkcjami.

Niedawno opisano nową technikę pozwalającą na bezpośrednie sekwencjonowanie m6A. Opiera się na wykorzystaniu specyficznych wariantów odwrotnej transkryptazy, które wykazują unikalne sygnatury RT w odpowiedzi na napotkanie modyfikacji RNA m6A24. Technologia ta, po starannej optymalizacji, może obejść główne ograniczenie stojące przed MeRIP-Seq (zmniejszenie ilości materiału wyjściowego i umożliwienie wyższej rozdzielczości). Aby zbadać modyfikację m5C, zdecydowaliśmy się użyć techniki konwersji wodorosiarczynu w celu wykrycia z rozdzielczością nukleotydową zmodyfikowanych reszt C. W celu zmniejszenia odsetka wyników fałszywie dodatnich ze względu na obecność struktur drugorzędowych RNA, przeprowadziliśmy 3 cykle denaturacji/konwersji wodorosiarczynu i dalej kontrolowaliśmy wydajność współczynnika konwersji wodorosiarczynu dzięki zastosowaniu kontroli ERCC z impulsami skokowymi. Jednym z ograniczeń związanych z tą techniką jest to, że konwersja wodorosiarczynu jest bardzo ostra, a trzy cykle denaturacji/konwersji wodorosiarczynu mogą zdegradować część RNA, a tym samym zmniejszyć rozdzielczość. Jednak w naszym otoczeniu zdecydowaliśmy się zadowolić potencjalnie nieco niższą rozdzielczością, aby zwiększyć jakość zestawu danych.

Dzięki tym optymalizacjom i kontrolom byliśmy w stanie zapewnić niezawodny i solidny przepływ pracy, który można wykorzystać do zbadania krajobrazu epitranskryptomicznego i jego zmian w kontekście infekcji wirusowych, interakcji gospodarz-patogen lub jakiejkolwiek ekspozycji na określone terapie.

Oświadczenia

Autorzy nie mają nic do ujawnienia.

Podziękowania

Ta praca była wspierana przez Szwajcarską Narodową Fundację Nauki (granty 31003A_166412 i 314730_188877).

Materiały

Lista materiałów użytych w tym artykule
NazwaFirmaNumer katalogowyKomentarze
AccuPrime Pfx SuperMixInvitrogen12344-040
Przeciwciało anty-m6A _Clone 17-3-4-1MilliporeMABE1006
ChloroformMerck67-66-3
ERCCInvitrogen4456740
EZZymo ResearchEZR5001
Zestaw RNA - Zestaw HS RNAAgilentDNF-472-0500
Zestaw analizatora fragmentów RNA - Zestaw RNAAgilentDNF-471-0500
Zestaw do odwrotnej transkrypcji cDNA o dużej pojemnościApplied Biosystem4368814
Illumina TruSeq Nici mRNAIllumina20020594
Kulki magnetyczne Mieszanka A/GMerck16-663
N6-metyloadenozyna, sól sodowa 5′-monofosforanu (m6A)Sigma AldrichM2780-10MG
Normalna mysia IgGMerk12371
Oligo(dT)25Life Technologies61005,
PCRapaceStratec 1020220300
Szybki zestaw wirusowy RNAZymo Research1034
RNA Czystość & KoncentratorZymo ResearchR1015
Odczynnik do fragmentacji RNAAmbionAM8740
Inhibitor RNazyAmbionAM2684
Trizol TRIzolOdczynnik 15596026
Zestaw do analizy fragmentów

Bibliografia

  1. Machnicka, M. A., et al. MODOMICS: a database of RNA modification pathways--2013 update. Nucleic Acids Research. 41, Database issue 262-267 (2013).
  2. Zaccara, S., Ries, R. J., Jaffrey, S. R. Reading, writing and erasing mRNA methylation. Nature Reviews Molecular Cell Biology. 20 (10), 608-624 (2019).
  3. Davalos, V., Blanco, S., Esteller, M. SnapShot: Messenger RNA Modifications. Cell. 174 (2), 498(2018).
  4. Saletore, Y., et al. The birth of the Epitranscriptome: deciphering the function of RNA modifications. Genome Biology. 13 (10), 175(2012).
  5. Zhao, B. S., Roundtree, I. A., He, C. Post-transcriptional gene regulation by mRNA modifications. Nature Reviews Molecular Cell Biology. 18 (1), 31-42 (2017).
  6. Netzband, R., Pager, C. T. Epitranscriptomic marks: Emerging modulators of RNA virus gene expression. Wiley Interdisciplinary Reviews: RNA. 11 (3), 1576(2020).
  7. Pereira-Montecinos, C., Valiente-Echeverria, F., Soto-Rifo, R. Epitranscriptomic regulation of viral replication. Biochimica et Biophysica Acta. 1860 (4), 460-471 (2017).
  8. Lichinchi, G., et al. Dynamics of the human and viral m(6)A RNA methylomes during HIV-1 infection of T cells. Nature Microbiology. 1, 16011(2016).
  9. Courtney, D. G., et al. Epitranscriptomic Addition of m(5)C to HIV-1 Transcripts Regulates Viral Gene Expression. Cell Host & Microbe. 26 (2), 217-227 (2019).
  10. Kennedy, E. M., et al. Posttranscriptional m(6)A Editing of HIV-1 mRNAs Enhances Viral Gene Expression. Cell Host & Microbe. 19 (5), 675-685 (2016).
  11. Tirumuru, N., Wu, L. HIV-1 envelope proteins up-regulate N (6)-methyladenosine levels of cellular RNA independently of viral replication. Journal of Biological Chemistry. 294 (9), 3249-3260 (2019).
  12. Tirumuru, N., et al. N(6)-methyladenosine of HIV-1 RNA regulates viral infection and HIV-1 Gag protein expression. Elife. 5, (2016).
  13. Cristinelli, S., Angelino, P., Janowczyk, A., Delorenzi, M., Ciuffi, A. HIV Modifies the m6A and m5C Epitranscriptomic Landscape of the Host Cell. Frontiers in Virology. 1 (11), (2021).
  14. Khoddami, V., Cairns, B. R. Transcriptome-wide target profiling of RNA cytosine methyltransferases using the mechanism-based enrichment procedure Aza-IP. Nature Protocols. 9 (2), 337-361 (2014).
  15. Hussain, S., Aleksic, J., Blanco, S., Dietmann, S., Frye, M. Characterizing 5-methylcytosine in the mammalian epitranscriptome. Genome Biology. 14 (11), 215(2013).
  16. Dominissini, D., Moshitch-Moshkovitz, S., Salmon-Divon, M., Amariglio, N., Rechavi, G. Transcriptome-wide mapping of N6-methyladenosine by m6A-seq based on immunocapturing and massively parallel sequencing. Nature Protocols. 8 (1), 176-189 (2013).
  17. Dominissini, D., et al. Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq. Nature. 485 (7397), 201-206 (2012).
  18. Shobbir Hussain, J. A., Blanco, S., Dietmann, S., Frye, M. Characterizing 5-methylcytosine in the mammalian epitranscriptome. Genome Biology. 14 (215), (2013).
  19. Amort, T., et al. Distinct 5-methylcytosine profiles in poly(A) RNA from mouse embryonic stem cells and brain. Genome Biology. 18 (1), 1(2017).
  20. Endrullat, C., Glökler, J., Franke, P., Frohme, M. Standardization and quality management in next-generation sequencing. Applied & Translational Genomics. 10, 2-9 (2016).
  21. Schaefer, M., Pollex, T., Hanna, K., Lyko, F. RNA cytosine methylation analysis by bisulfite sequencing. Nucleic Acids Research. 37 (2), 12(2009).
  22. Cristinelli, S., Angelino, P., Janowczyk, A., Delorenzi, M., Ciuffi, A. HIV Modifies the m6A and m5C Epitranscriptomic Landscape of the Host Cell. biorxiv. 1 (11), (2021).
  23. Squires, J. E., et al. Widespread occurrence of 5-methylcytosine in human coding and noncoding RNA. Nucleic Acids Research. 40 (11), 5023-5033 (2012).
  24. Aschenbrenner, J., et al. Engineering of a DNA Polymerase for Direct m(6) A Sequencing. Angewandte Chemie (International ed. in English). 57 (2), 417-421 (2018).
  25. Bioinformatics. , Available from: http://www.bioinformatics.babraham.ac.uk/projects/fastqc (2021).
  26. Didion, J. P., Martin, M., Collins, F. S. Atropos: specific, sensitive, and speedy trimming of sequencing reads. PeerJ. 5, 3720(2017).
  27. Kim, D., Langmead, B., Salzberg, S. L. HISAT: a fast spliced aligner with low memory requirements. Nature Methods. 12 (4), 357-360 (2015).
  28. Li, H., et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 25 (16), 2078-2079 (2009).
  29. Okonechnikov, K., Conesa, A., García-Alcalde, F. Qualimap 2: advanced multi-sample quality control for high-throughput sequencing data. Bioinformatics. 32 (2), Oxford, England. 292-294 (2016).
  30. Ewels, P., Magnusson, M., Lundin, S., Käller, M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 32 (19), Oxford, England. 3047-3048 (2016).
  31. Zhang, Y., et al. Model-based Analysis of ChIP-Seq (MACS). Genome Biology. 9 (9), 137(2008).
  32. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet Journal. 17 (1), (2011).
  33. Shen, W., Le, S., Li, Y., Hu, F. SeqKit: A Cross-Platform and Ultrafast Toolkit for FASTA/Q File Manipulation. PLOS ONE. 11 (10), 0163962(2016).
  34. Rieder, D., Amort, T., Kugler, E., Lusser, A., Trajanoski, Z. meRanTK: methylated RNA analysis ToolKit. Bioinformatics. 32 (5), 782-785 (2015).

Przedruki i uprawnienia

Tagi

Znaczniki epitranskryptomicznemodyfikacje RNAzakażenie HIVmetylacja m6Ametylacja m5CMeRIP-Seqsekwencjonowanie z bisulfitamimetylacja różnicowaizolacja mRNA
Badanie epitranskryptomów m6A i m5C po infekcji wirusowej: przykład z HIV
Skontaktuj się z nami
Poleć bibliotece

Prywatność

Warunki korzystania z usługi

Zasady

Prawa autorskie © 2026 MyJoVE Corporation. Wszelkie prawa zastrzeżone.