Modułowa konstrukcja Promoter Capture Hi-C
Promoter Capture Hi-C został zaprojektowany w celu wzbogacenia bibliotek Hi-C dla interakcji z udziałem promotorów. Interakcje te obejmują tylko podzbiór produktów ligacji obecnych w bibliotece Hi-C.
Capture Hi-C można łatwo zmodyfikować w celu wzbogacenia bibliotek Hi-C dla dowolnego regionu genomowego lub regionów zainteresowania poprzez zmianę systemu przechwytywania. Regiony wychwytywania mogą być ciągłymi segmentami genomu 44,45,46,48, wzmacniaczami, które zostały zidentyfikowane w PCHi-C ("Reverse Capture Hi-C'35") lub miejscach nadwrażliwych na DNazę I49. Rozmiar systemu przechwytywania można dostosować w zależności od zakresu eksperymentalnego. Na przykład Dryden i wsp. Celuj w 519 fragmentów przynęty w trzech pustyniach genowych związanych z rakiem piersi44. System wychwytywania Martina i wsp. jest ukierunkowany zarówno na ciągłe segmenty genomu ("przechwytywanie regionów": łącznie 211 regionów genomu; 2 131 fragmentów restrykcyjnych), jak i wybrane promotory (3 857 promotorów genów)45.
Biblioteki SureSelect są dostępne w różnych zakresach rozmiarów: od 1 kb do 499 kb (5 190–4 806), od 500 kb do 2,9 MB (5 190–4 816) i od 3 Mb do 5,9 Mb (5 190–4 831). Ponieważ każdy pojedynczy biotyna-RNA ma długość 120 nukleotydów, te systemy wychwytywania mieszczą maksymalnie 4 158, 24 166 i 49 166 indywidualnych sond wychwytujących. Odpowiada to odpowiednio 2 079, 12 083 i 24 583 ukierunkowanym fragmentom restrykcyjnym (należy zauważyć, że liczby dla fragmentów ograniczających są dolnymi granicami opartymi na założeniu, że dla każdego fragmentu restrykcji można zaprojektować dwie indywidualne sondy przechwytujące — w rzeczywistości ze względu na powtarzające się sekwencje nie będzie tak w przypadku każdego fragmentu restrykcji (zob. również rysunek 1B, C), co skutkuje większą liczbą możliwych do namierzenia fragmentów ograniczających dla stałej liczby dostępnych sond przechwytujących).
Opisany tutaj protokół opiera się na wykorzystaniu enzymu restrykcyjnego z miejscem rozpoznawania 6 pz w celu odkrycia interakcji dalekiego zasięgu. Możliwe jest również użycie enzymu restrykcyjnego z miejscem rozpoznawania 4 pz w celu uzyskania większej rozdzielczości bardziej proksymalnych interakcji40,49.
Ograniczenia PCHi-C
Jednym z nieodłącznych ograniczeń wszystkich testów wychwytywania konformacji chromosomów jest to, że ich rozdzielczość jest określana przez enzym restrykcyjny używany do generowania biblioteki. Interakcje zachodzące między elementami DNA znajdującymi się na tym samym fragmencie restrykcyjnym są niewidoczne dla testów typu C. Co więcej, w PCHi-C w niektórych przypadkach więcej niż jedno miejsce rozpoczęcia transkrypcji może znajdować się na tym samym fragmencie restrykcyjnym zawierającym promotor, a PIR w niektórych przypadkach zawierają zarówno aktywne, jak i represyjne znaczniki histonów, co utrudnia określenie, które elementy regulatorowe pośredniczą w interakcjach, oraz przewidywanie regulacyjnego wyniku interakcji promotora. Stosowanie enzymów restrykcyjnych z miejscami rozpoznawania 4 pz łagodzi ten problem, ale odbywa się kosztem znacznie zwiększonej złożoności biblioteki Hi-C (biblioteki Hi-C wygenerowane za pomocą enzymów ograniczających miejsce rozpoznawania 4 pz są co najmniej 100 razy bardziej złożone niż biblioteki Hi-C wygenerowane za pomocą enzymów restrykcyjnych miejsca rozpoznawania 6 pz) i związanych z tym kosztów sekwencjonowania nowej generacji.
Innym ograniczeniem jest to, że obecny protokół PCHi-C wymaga milionów komórek jako materiału wyjściowego, co uniemożliwia analizę interakcji promotorów w rzadkich typach komórek. Zmodyfikowana wersja PCHi-C, umożliwiająca badanie kontaktów promotora w populacjach komórek liczących od 10 000 do 100 000 komórek (na przykład komórek we wczesnym okresie rozwoju embrionalnego lub hematopoetycznych komórek macierzystych), byłaby zatem cennym uzupełnieniem zestawu narzędzi Capture Hi-C.
Wreszcie, podobnie jak wszystkie metody, które opierają się na utrwalaniu formaldehydu, PCHi-C rejestruje tylko te interakcje, które są "zamrożone" w momencie utrwalenia. W związku z tym, aby zbadać kinetykę i dynamikę oddziaływań promotorów, potrzebne są metody takie jak mikroskopia żywych komórek o wysokiej rozdzielczości wraz z PCHi-C.
Metody analizy przestrzennej organizacji chromosomów w wysokiej rozdzielczości
Ogromna złożoność bibliotek interakcji chromosomowych uniemożliwia wiarygodną identyfikację produktów interakcji między dwoma określonymi fragmentami restrykcyjnymi z istotnością statystyczną. Aby obejść ten problem, przechwytywanie sekwencji zostało wykorzystane do wzbogacenia bibliotek Hi-C 33,34,40,44 lub 3C50,51 dla określonych interakcji. Główną zaletą korzystania z bibliotek Hi-C w porównaniu z bibliotekami 3C na etapie wzbogacania jest to, że Hi-C, w przeciwieństwie do 3C, zawiera etap wzbogacania dla oryginalnych produktów ligacji. W konsekwencji odsetek prawidłowych odczytów w bibliotekach PCHi-C jest około 10-krotnie wyższy niż w bibliotekach Capture-C50, które zawierały około 5–8% prawidłowych odczytów po filtrowaniu HiCUP. Sahlen i in. bezpośrednio porównali Capture-C z HiCap, który podobnie jak PCHi-C używa bibliotek Hi-C do wzbogacania przechwytywania, w przeciwieństwie do Capture-C, który używa bibliotek 3C. Zgodnie z naszymi ustaleniami odkryli, że biblioteki Capture-C składają się głównie z nieligowanych fragmentów40. Ponadto biblioteki HiCap miały większą złożoność niż biblioteki Capture-C40.
Wariant Capture-C, zwany Capture-C52 nowej generacji (NG Capture-C) wykorzystuje jeden oligo na koniec fragmentu restrykcyjnego, jak wcześniej ustalono w PCHi-C 33,34, zamiast nakładających się sond używanych w oryginalnym protokole Capture-C50. Zwiększa to nieznacznie odsetek prawidłowych odczytów w porównaniu z Capture-C, ale NG Capture-C wykorzystuje dwie sekwencyjne rundy wzbogacania przechwytywania i stosunkowo dużą liczbę cykli PCR (łącznie od 20 do 24 cykli, w porównaniu do 11 cykli zwykle dla PCHi-C), co nieuchronnie skutkuje większą liczbą duplikatów sekwencji i mniejszą złożonością biblioteki. W eksperymentach próbnych podczas optymalizacji PCHi-C stwierdziliśmy, że odsetek unikalnych (tj. nie zduplikowanych) par odczytu wynosił tylko około 15%, gdy użyliśmy 19 cykli PCR (13 cykli przed wychwytywaniem + 6 cykli po przechwyceniu; dane nie pokazane), jednak optymalizacja do mniejszej liczby cykli PCR zwykle daje 75-90% unikalnych par odczytu. W związku z tym zmniejszenie liczby cykli PCR znacznie zwiększa ilość informacyjnych danych sekwencyjnych.
Najnowsza metoda łączy ChIP z Hi-C, aby skupić się na interakcjach chromosomowych, w których pośredniczy określone białko będące przedmiotem zainteresowania (HiChIP53). W porównaniu z ChIA-PET54, który opiera się na podobnym przesłaniu, dane HiChIP zawierają większą liczbę odczytów sekwencji informacyjnych, co pozwala na wywołanie interakcji o wyższym poziomie ufności53. Bardzo interesujące będzie bezpośrednie porównanie odpowiednich zestawów danych HiChIP i Capture Hi-C, gdy staną się one dostępne (na przykład HiChIP przy użyciu przeciwciała przeciwko jednostce kohezyny Smc1a53 z Capture Hi-C dla wszystkich fragmentów restrykcyjnych związanych z Smc1a) obok siebie. Jedną z nieodłącznych różnic między tymi dwoma podejściami jest to, że Capture Hi-C nie opiera się na immunoprecypitacji chromatyny, a zatem jest zdolny do badania interakcji chromosomowych niezależnie od zajętości białka. Umożliwia to porównanie organizacji genomu 3D w obecności lub braku wiązania określonego czynnika, co zostało wykorzystane do identyfikacji PRC1 jako kluczowego regulatora architektury przestrzennej genomu myszy ESC7.
PCHi-C i GWAS
Badania asocjacyjne całego genomu (GWAS) wykazały, że ponad 95% wariantów sekwencji związanych z chorobą znajduje się w niekodujących regionach genomu, często w dużej odległości od genów kodujących białka55. Warianty GWAS często znajdują się w bliskim sąsiedztwie miejsc nadwrażliwych na DNazę I, co jest cechą charakterystyczną sekwencji o potencjalnej aktywności regulatorowej. PCHi-C i Capture Hi-C były szeroko stosowane do łączenia promotorów z loci ryzyka GWAS związanymi z rakiem piersi44, rakiem jelita grubego48 i chorobą autoimmunologiczną 35,45,46. Badanie PCHi-C na 17 różnych typach ludzkich komórek krwiotwórczych wykazało, że SNP związane z chorobą autoimmunologiczną były wzbogacone w PIR w komórkach limfoidalnych, podczas gdy warianty sekwencji związane z cechami specyficznymi dla płytek krwi i czerwonych krwinek stwierdzono głównie odpowiednio w makrofagach i erytroblastach35,56. W związku z tym specyficzne dla typu tkanki interakcje promotorowe odkryte przez PCHi-C mogą pomóc w zrozumieniu funkcji niekodujących wariantów sekwencji związanych z chorobą i zidentyfikowaniu nowych potencjalnych genów chorobowych do interwencji terapeutycznej.
Charakterystyka regionów oddziałujących na siebie z promotorem
Istnieje kilka linii dowodowych łączących interakcje promotorów z kontrolą ekspresji genów. Po pierwsze, kilka badań PCHi-C wykazało, że regiony genomu oddziałujące z promotorami genów o (wysokiej) ekspresji są wzbogacone w znaczniki związane z aktywnością wzmacniającą, takie jak acetylacja H3K27 i wiązanie p300 33,34,37. Znaleźliśmy pozytywną korelację między poziomem ekspresji genów a liczbą oddziałujących wzmacniaczy, co sugeruje, że addytywne działanie wzmacniaczy powoduje zwiększenie poziomu ekspresji genów34,35. Po drugie, naturalnie występujące loci ilościowych cech ekspresji (eQTL) są wzbogacone w PIR, które są połączone z tymi samymi genami, na których ekspresję wpływają eQTL35. Po trzecie, integrując dane TRIP57 i PCHi-C, Cairns i in. odkryli, że geny reporterowe TRIP mapujące do PIR w mysich ESC wykazują silniejszą ekspresję genów reporterowych niż geny reporterowe w miejscach integracji w regionach nieoddziałujących z promotorem58, co wskazuje, że PIRR posiadają transkrypcyjną aktywność regulacyjną. Łącznie odkrycia te sugerują, że interaktomy promotorowe odkryte przez PCHi-C w różnych typach komórek mysich i ludzkich obejmują kluczowe moduły regulacyjne do kontroli ekspresji genów.
Warto zauważyć, że wzmacniacze stanowią tylko niewielki ułamek (~20%) wszystkich PIR-ów odkrytych przez PCHi-C33,34. Inne PIR mogą pełnić rolę strukturalną lub topologiczną, a nie bezpośrednie transkrypcyjne funkcje regulacyjne. Istnieją jednak również dowody na to, że PCHi-C może odkrywać elementy DNA o funkcji regulacyjnej, które nie zawierają klasycznych znamion wzmacniacza. W ludzkiej linii komórek limfoidalnych stwierdzono, że promotor BRD7 oddziałuje z regionem pozbawionym śladów wzmacniacza, który okazał się mieć aktywność wzmacniającą w testach genów reporterowych33. Elementy regulacyjne o podobnych cechach mogą być liczniejsze niż obecnie się uważa. Na przykład badanie przesiewowe regulatorowe DNA oparte na CRISPR zidentyfikowało nieoznaczone elementy regulatorowe (URE), które kontrolują ekspresję genów, ale są pozbawione znaczników wzmacniających59.
W innych przypadkach wykazano, że PIR zawierają znaczniki chromatyny związane z represją transkrypcji. PIR i oddziałujące promotory związane przez PRC1 w mysich ESC były zaangażowane w rozległą sieć przestrzenną stłumionych genów noszących represyjny znak H3K27me37. W ludzkich komórkach limfoblastoidalnych odległy element oddziałujący z promotorem BCL6 tłumił ekspresję genu reporterowegotransgenu 33, co sugeruje, że może on działać w celu tłumienia transkrypcji BCL6 w swoim natywnym kontekście.
Innąbrację PIR wzbogaconą o obecność białka izolatora chromatyny CTCF w ludzkich ESC i NEC37 mogą reprezentować jeszcze inną klasę PIR. Podsumowując, wyniki te sugerują, że PIR zawierają zbiór czynności regulacyjnych genów, które nie zostały jeszcze scharakteryzowane funkcjonalnie.