R-loops to trójniciowe struktury kwasów nukleinowych, które tworzą się głównie podczas transkrypcji po hybrydyzacji powstającego transkryptu RNA z matrycową nicią DNA. Powoduje to powstanie hybrydy RNA:DNA i powoduje przemieszczenie nici DNA bez matrycy w stanie jednoniciowym zapętlonym. Rekonstytucja biochemiczna1,2,3,4 i modelowanie matematyczne5, w połączeniu z innymi pomiarami biofizycznymi6,7, ustaliły, że prawdopodobieństwo wystąpienia pętli R w regionach, które wykazują określone korzystne cechy. Na przykład regiony, które wykazują asymetrię nici w dystrybucji guaniny (G) i cytozyny (C) tak, że RNA jest bogate w G, właściwość zwana dodatnim skosem GC, są faworyzowane do tworzenia pętli R podczas transkrypcji ze względu na wyższą stabilność termodynamiczną hybrydy DNA: RNA w porównaniu z klasą DNA duplex8. Regiony, które wyewoluowały dodatnie skośność GC, takie jak wczesne części wielu genów eukariotycznych4,9,10,11, są podatne na tworzenie pętli R in vitro i in vivo3,4,12. Ujemne naprężenie superhelical DNA również bardzo sprzyja tworzeniu struktury13,14, ponieważ pętle R skutecznie absorbują takie naprężenia topologiczne i przywracają otaczające włókno DNA do korzystnego stanu relaksu5,15.
Historycznie rzecz biorąc, uważano, że struktury pętli R są wynikiem rzadkich, spontanicznych splątań RNA z DNA podczas transkrypcji. Jednak rozwój immunoprecypitacji DNA:RNA (DRIP) w połączeniu z wysokoprzepustowym sekwencjonowaniem DNA (DRIP-seq) pozwolił na pierwsze mapowanie pętli R w całym genomie i ujawnił, że struktury te są znacznie bardziej rozpowszechnione niż oczekiwano w komórkach ludzkich4,16. Pętle R występują w dziesiątkach tysięcy konserwatywnych, transkrybowanych, genowych hotspotów w genomach ssaków, z predylekcją do wysp CpG skośnych GC, które nakładają się na pierwszy intron genów i końcowe regiony wielu genów17. Ogólnie rzecz biorąc, pętle R łącznie zajmują 3%-5% genomu w komórkach ludzkich, co jest zgodne z pomiarami w innych organizmach, w tym drożdżach, roślinach, muchach i myszach18,19,20,21,22.
Analiza gorących punktów tworzących pętlę R w ludzkich komórkach wykazała, że takie regiony kojarzą się ze specyficznymi sygnaturami chromatyny23. Ogólnie rzecz biorąc, pętle R znajdują się w regionach o niższym zajętości nukleosomów i wyższej gęstości polimerazy RNA. U promotorów pętle R wiążą się ze zwiększoną rekrutacją dwóch kotranskrypcyjnie zdeponowanych modyfikacji histonów, H3K4me1 i H3K36me317. Na końcach genów pętle R wiążą się z blisko ułożonymi genami, które przechodzą skuteczną terminację transkrypcji17, zgodnie z wcześniejszymi obserwacjami24. Wykazano również, że pętle R uczestniczą w inicjacji replikacji DNA w początkach replikacji genomów bakteriofagów, plazmidów, mitochondriów i drożdży25,26,27,28,29,30,31. Ponadto 76% podatnych na pętlę R ludzkich promotorów wysp CpG funkcjonuje jako wczesne, konstytutywne źródła replikacji32,33,34,35, co dodatkowo wzmacnia powiązania między pętlami R a początkami replikacji. Podsumowując, badania te sugerują, że pętle R reprezentują nowy typ sygnału biologicznego, który może wyzwalać określone wyjścia biologiczne w sposób zależny od kontekstu23.
Na początku wykazano, że pętle R tworzą się w sekwencjach przełączania klas podczas procesu rekombinacji przełączników klasy immunoglobulin3,36,37. Uważa się, że takie zaprogramowane pętle R inicjują rekombinację przełączników klas poprzez wprowadzenie dwuniciowych pęknięć DNA38. Od tego czasu szkodliwe powstawanie pętli R, ogólnie rozumiane jako wynikające z nadmiernego tworzenia pętli R, zostało powiązane z niestabilnością genomu i procesami, takimi jak hiperrekombinacja, kolizje transkrypcji-replikacji, replikacja i stres transkrypcyjny (do przeglądu39,40,41,42,43). W związku z tym ulepszone mapowanie struktur pętli R stanowi ekscytujące i zasadnicze wyzwanie dla lepszego rozszyfrowania rozmieszczenia i funkcji tych struktur w zdrowiu i chorobie.
DNA:RNA immunoprecypitacja (DRIP) opiera się na wysokim powinowactwie przeciwciała monoklonalnego S9.6 do hybryd DNA:RNA44. DRIP-seq umożliwia solidne profilowanie tworzenia pętli R w całym genomie4,45. Chociaż jest to przydatne, technika ta ma ograniczoną rozdzielczość ze względu na fakt, że enzymy restrykcyjne są używane do uzyskania delikatnej fragmentacji DNA. Ponadto sekwencja DRIP nie dostarcza informacji na temat kierunkowości powstawania pętli R. W tym miejscu przedstawiamy wariant sekwencji DRIP, który umożliwia mapowanie pętli R w wysokiej rozdzielczości w sposób specyficzny dla pasma. Metoda ta opiera się na sonikacji w celu fragmentacji genomu przed immunoprecypitacją, a zatem metoda ta nazywa się sDRIP-seq (sonikacja DNA: immunoprecypitacja RNA sprzężona z sekwencjonowaniem o wysokiej przepustowości) (Figura 1). Zastosowanie sonikacji pozwala na zwiększenie rozdzielczości i ogranicza ograniczenia związane z fragmentacją związaną z enzymami obserwowane w podejściach DRIP-seq46. sDRIP-seq tworzy mapy pętli R, które są silnie zgodne z wynikami zarówno DRIP-seq, jak i wcześniej opisanej metody DRIPc-seq o wysokiej rozdzielczości, w której biblioteki sekwencjonowania są budowane z nici RNA immunoprecypitowanych struktur pętli R45.
W obliczu mnóstwa metod do wyboru, użytkownicy mogą zastanawiać się, które konkretne podejście oparte na DRIP jest preferowane dla ich potrzeb. Oferujemy następujące porady. DRIP-seq, pomimo swoich ograniczeń, jest technicznie najłatwiejszy i jest najbardziej wytrzymały (najwyższe plony) ze wszystkich trzech omawianych tutaj metod; W związku z tym pozostaje on szeroko użyteczny. Opublikowano liczne zestawy danych DRIP-seq, które stanowią użyteczny punkt odniesienia dla nowych zestawów danych. Wreszcie, proces analizy bioinformatycznej jest prostszy, ponieważ dane nie są pozostawione na lodzie. Zaleca się, aby nowi użytkownicy zaczęli doskonalić swoje umiejętności mapowania pętli R za pomocą DRIP, a następnie ilościowej reakcji łańcuchowej polimerazy (qPCR) i sekwencji DRIP. sDRIP-seq reprezentuje nieco wyższy stopień trudności technicznej: wydajność jest nieznacznie zmniejszona z powodu sonikacji (omówionej poniżej), a proces biblioteki sekwencjonowania jest nieco bardziej złożony. Jednak zysk wynikający z utknięcia na mieliźnie i wyższej rozdzielczości jest nieoceniony. Należy zauważyć, że sDRIP-seq wychwyci zarówno dwuniciowe hybrydy RNA:DNA, jak i trójniciowe pętle R. Ze względu na etapy budowy biblioteki, DRIP-seq nie wychwyci dwuniciowych hybryd RNA:DNA. DRIPc-seq jest najbardziej wymagający technicznie i wymaga większej ilości materiałów wyjściowych. W zamian oferuje najwyższą rozdzielczość i osierocenia. Ponieważ biblioteki sekwencjonowania są zbudowane z części RNA pętli R lub hybryd, sekwencja DRIPc może ucierpieć z powodu możliwego zanieczyszczenia RNA, zwłaszcza, że S9.6 ma szczątkowe powinowactwo do dsRNA19,47,48. sDRIP-seq umożliwia mapowanie specyficzne dla nici w wysokiej rozdzielczości bez obaw o zanieczyszczenie RNA, ponieważ biblioteki sekwencjonowania pochodzą z nici DNA. Ogólnie rzecz biorąc, te trzy metody pozostają użyteczne i prezentują różne stopnie złożoności oraz nieco inne zastrzeżenia. Wszystkie trzy generują jednak wysoce spójne zestawy danych48 i są bardzo wrażliwe na wstępne traktowanie RNazy H, co stanowi istotną kontrolę zapewniającą specyficzność sygnału45,49. Należy zauważyć, że biorąc pod uwagę wybór wielkości narzucony bibliotekom sekwencjonowania, małe hybrydy (szacowane <75 pz), takie jak te, które tworzą się przejściowo wokół miejsc inicjacji replikacji replikacji DNA o opóźnionej nici (startery Okazaki) zostaną wykluczone. Podobnie, ponieważ wszystkie metody DRIP obejmują fragmentację DNA, niestabilne pętle R, które wymagają ujemnego superzwijania DNA dla swojej stabilności, zostaną utracone5. W związku z tym podejścia DRIP mogą zaniżać obciążenia pętli R, szczególnie w przypadku krótkich, niestabilnych pętli R, które najlepiej uchwycić za pomocą podejść in vivo45,48. Należy zauważyć, że pętle R mogą być również profilowane w sposób niezależny od S9.6 przy głębokim pokryciu, wysokiej rozdzielczości i w sposób specyficzny dla nici na pojedynczych cząsteczkach DNA po traktowaniu wodorosiarczynu sodu12. Dodatkowo, strategie wykorzystujące katalitycznie nieaktywny enzym RNaza H1 zostały wykorzystane do mapowania natywnych pętli R in vivo, podkreślając krótkie, niestabilne pętle R, które tworzą się głównie w wstrzymanych promotorach50,51,52.