Hi-C 3.0 to ulepszony protokół Hi-C, który łączy sieciujące formaldehyd i glutaran disuccinimidylu z koktajlem enzymów restrykcyjnych DpnII i DdeI, aby zwiększyć stosunek sygnału do szumu i rozdzielczość wykrywania interakcji chromatyny.
Artykuł metodologiczny
Hi-C 3.0 to ulepszony protokół Hi-C, który łączy sieciujące formaldehyd i glutaran disuccinimidylu z koktajlem enzymów restrykcyjnych DpnII i DdeI, aby zwiększyć stosunek sygnału do szumu i rozdzielczość wykrywania interakcji chromatyny.
Przechwytywanie konformacji chromosomów (3C) służy do wykrywania trójwymiarowych oddziaływań chromatyny. Zazwyczaj chemiczne sieciowanie formaldehydem (FA) stosuje się do utrwalania oddziaływań chromatyny. Następnie trawienie chromatyny enzymem restrykcyjnym i późniejsza religacja końców fragmentów przekształca trójwymiarową (3D) bliskość w unikalne produkty ligacji. Wreszcie, po odwróceniu wiązań krzyżowych, usunięciu białka i izolacji DNA, DNA jest ścinane i przygotowywane do sekwencjonowania o wysokiej przepustowości. Częstość ligacji zbliżeniowej par loci jest miarą częstości ich kolokalizacji w przestrzeni trójwymiarowej w populacji komórek.
Sekwencjonowana biblioteka Hi-C dostarcza informacji o częstości interakcji między wszystkimi parami loci dla całego genomu. Rozdzielczość i precyzja Hi-C opiera się na wydajnym sieciowaniu, które utrzymuje kontakty chromatyny oraz częstą i równomierną fragmentację chromatyny. W artykule opisano ulepszony protokół Hi-C in situ, Hi-C 3.0, który zwiększa efektywność sieciowania poprzez połączenie dwóch środków sieciujących (formaldehydu [FA] i glutaranu disuccinimidylu [DSG]), a następnie dokładniejsze trawienie przy użyciu dwóch enzymów restrykcyjnych (DpnII i DdeI). Hi-C 3.0 to pojedynczy protokół do dokładnej kwantyfikacji cech fałdowania genomu w mniejszych skalach, takich jak pętle i domeny kojarzenia topologicznego (TAD), a także cech w większych skalach obejmujących jądro, takich jak przedziały.
Metoda uchwytywania konformacji chromosomów (Chromosome conformation capture) jest stosowana od 2002 roku1. Fundamentalnie każdy wariant uchwytywania konformacji opiera się na utrwaleniu oddziaływań DNA-białko oraz białko-białko w celu zachowania trójwymiarowej organizacji chromatyny. Następnie przeprowadza się fragmentację DNA, zazwyczaj poprzez trawienie restrykcyjne, a na koniec religację bliskich końców DNA, aby przekształcić loci bliskie przestrzennie w unikalne kowalencyjne sekwencje DNA. W początkowych protokołach 3C wykorzystywano PCR do badania specyficznych oddziaływań typu „jeden-do-jednego”. Kolejne testy 4C umożliwiły wykrywanie oddziaływań typu „jeden-do-wszystkich”2, natomiast 5C wykrywało oddziaływania typu „wiele-do-wielu”3. Metoda uchwytywania konformacji chromosomów w pełni rozwinęła się po wprowadzeniu wysokoprzepustowego sekwencjonowania nowej generacji (NGS), co pozwoliło na wykrywanie genomowych oddziaływań typu „wszystkie-do-wszystkich” przy użyciu ogennomowej metody Hi-C4 oraz technik porównywalnych, takich jak 3C-seq5, TCC6 i Micro-C7,8 (zob. również przegląd Denkera i De Laata9).
W metodzie Hi-C do znakowania nawisów 5′ po trawieniu i przed ligacją stosuje się nukleotydy biotynilowane (Rycina 1). Pozwala to na wyselekcjonowanie prawidłowo strawionych i religowanych fragmentów za pomocą kulek pokrytych streptawidyną, co odróżnia tę metodę od GCC10. Istotną aktualizację protokołu Hi-C wprowadzili Rao i wsp.11, którzy przeprowadzili trawienie i religację w nienaruszonych jądrach komórkowych (tzn. in situ), aby zredukować liczbę fałszywych produktów ligacji. Ponadto zastąpienie trawienia enzymem HindIII trawieniem MboI (lub DpnII) zmniejszyło rozmiar fragmentów i zwiększyło potencjalną rozdzielczość Hi-C. Wzrost ten umożliwił wykrywanie struktur o stosunkowo niewielkiej skali oraz bardziej precyzyjną lokalizację genomową punktów kontaktu, takich jak pętle DNA między małymi elementami cis, np. pętle między miejscami wiązania CTCF powstające w wyniku ekstruzji pętli11,12. Ten potencjał wiąże się jednak z pewnymi kosztami. Po pierwsze, dwukrotny wzrost rozdzielczości wymaga czterokrotnego (22) zwiększenia liczby odczytów sekwencjonowania13. Po drugie, niewielkie rozmiary fragmentów zwiększają prawdopodobieństwo pomylenia niestrawionych fragmentów sąsiednich z fragmentami strawionymi i religowanymi14. Jak wspomniano, w Hi-C fragmenty strawione i religowane różnią się od fragmentów niestrawionych obecnością biotyny w miejscu ligacji. Niemniej jednak, wymagane jest prawidłowe usunięcie biotyny z nieligowanych końców, aby zapewnić wyizolowanie wyłącznie miejsc ligacji14,15.
Wraz ze spadkiem kosztów NGS, możliwe staje się badanie fałdowania chromosomów w większych szczegółach. Aby zmniejszyć rozmiar fragmentów DNA, a tym samym zwiększyć rozdzielczość, protokół Hi-C można dostosować poprzez zastosowanie częściej tnących enzymów restrykcyjnych16 lub kombinacji enzymów restrykcyjnych17,18,19. Alternatywnie, w metodach Micro-C można miareczkować MNase7,8, a w DNase Hi-C miareczkować DNase20, aby uzyskać optymalny stopień trawienia.
Niedawna systematyczna ocena podstaw metod 3C wykazała, że detekcja cech fałdowania chromosomów we wszystkich skalach długości uległa znacznej poprawie dzięki sekwencyjnemu sieciowaniu 1% FA, a następnie 3 mM DSG17. Ponadto Hi-C z trawieniem za pomocą HindIII okazało się najlepszą opcją do wykrywania wielkoskalowych cech fałdowania, takich jak przedziały (compartments), natomiast Micro-C było skuteczniejsze w wykrywaniu małkoskalowych cech fałdowania, takich jak pętle DNA. Wyniki te doprowadziły do opracowania jednolitej strategii o wysokiej rozdzielczości „Hi-C 3.0”, która wykorzystuje kombinację sieciujących czynników FA i DSG, a następnie podwójne trawienie endonukleazami DpnII i DdeI21. Hi-C 3.0 stanowi skuteczną strategię do ogólnego zastosowania, ponieważ precyzyjnie wykrywa cechy fałdowania we wszystkich skalach długości17. Poniżej szczegółowo opisano część eksperymentalną protokołu Hi-C 3.0 oraz przedstawiono typowe wyniki, jakich można się spodziewać po sekwencjonowaniu.

Rysunek 1: Procedura Hi-C w sześciu krokach. Komórki są najpierw fiksowane za pomocą FA, a następnie DSG (1). Następnie po lizie następuje podwójne trawienie enzymami DdeI i DpnII (2). Biotyna jest wprowadzana poprzez wypełnianie nawisów, a proksymalne, stępione końce są ligowane (3) przed oczyszczaniem DNA (4). Biotyna jest usuwana z nieligowanych końców przed sonikacją i selekcją rozmiaru (5). Na koniec pull-down biotyny umożliwia ligację adapterów i amplifikację biblioteki metodą PCR (6). Skróty: FA = formaldehyd; DSG = glutaronian disukcynimidylu; B = biotyna. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.
1. Fiksacja poprzez sieciowanie
2. Przechwytywanie konformacji chromosomów
Tabela 1: Odczynniki do trawienia. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 2: Odczynniki do uzupełniania biotyny. *Należy pamiętać, że zmiana enzymów może wymagać zastosowania innych buforów i biotynylowanych dNTP. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 3: Odczynniki do mieszaniny ligacyjnej. Skróty: BSA = albuminy surowicy bydlęcej. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 4: Parametry nakładania na żel do oceny jakości i selekcji rozmiaru. Kliknij tutaj, aby pobrać tę tabelę.

Rycina 2: Żel agarozowy przedstawiający typowe wyniki kontroli jakości po oczyszczaniu DNA. (A) Kontrola CI powinna wykazać prążek DNA o wysokiej masie cząsteczkowej. (B) Próbki DC i Hi-C wykazują zakres wielkości DNA. Próbka Hi-C, w której fragmenty zostały połączone w większe jednostki, powinna mieć wyższą masę cząsteczkową niż DC. Zakres stężeń markerów umożliwia wyznaczenie krzywej wzorcowej. Należy zauważyć, że w tym przykładzie CI naniesiono na osobny żel, jednak zaleca się naniesienie i prowadzenie wszystkich próbek oraz kontroli jednocześnie. Skróty: CI = integralność chromatyny (chromatin integrity); DC = kontrola trawienia (digestion control); Hi-C = ligacja bliskościowa (proximity-ligated). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
3. Przygotowanie biblioteki do sekwencjonowania Hi-C
Tabela 5: Odczynniki i temperatury do usuwania biotyny Kliknij tutaj, aby pobrać tę tabelę.
Tabela 6: Parametry sonikacji. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 7: Odczynniki i temperatury do naprawy końców. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 8: Odczynniki i temperatury procesu A-tailingu. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 9: Startery PCR i oligosy adapterów typu paired-end wraz z odczynnikami do hybrydyzacji. Skrót: 5PHOS = fosforan 5'. Gwiazdki oznaczają fosforotioatowe zasady DNA. #Połącz indeksowany oligo z oligo uniwersalnym, aby uzyskać hybrydę w postaci indeksowanego adaptera. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 10: Odczynniki do ligacji adapterów. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 11: Odczynniki do PCR i parametry cyklu. Kliknij tutaj, aby pobrać tę tabelę.

Rycina 3: Żel agarozowy pokazujący typowe wyniki po selekcji wielkości fragmentów. Przedstawiono frakcje górną (Upper Fraction) i dolną (Lower Fraction) dla czterech próbek (ponumerowanych 1-4) Hi-C DpnII-DdeI. Pierwsza ścieżka dla każdej próbki zawiera frakcję górną, otrzymaną z mieszaniny kulek magnetycznych 0,8x, a druga i trzecia ścieżka zawierają rozcieńczenie frakcji dolnej otrzymanej z mieszaniny kulek magnetycznych 1,1x. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 4: Żel agarozowy z wynikami miareczkowania PCR. Począwszy od 5 cykli PCR, próbki pobierano co 2 cykle (5, 7, 9 i 11 cykli) dla każdej z czterech bibliotek. Na podstawie tego rysunku wybrano 6 cykli jako optymalną liczbę cykli dla każdej próbki. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 5: Końcowe produkty PCR. Po oczyszczaniu i selekcji wielkości, produkty PCR (Hi-C) nałożono obok frakcji tej samej biblioteki trawionej enzymem ClaI (ClaI). Fragmenty trawione ClaI wskazują na obecność poszukiwanych ligacji DpnII-DpnII. Należy zauważyć, że ClaI nie trawi złączy DpnII-DdeI, w związku z czym nie wszystkie ligacje przyczynią się do zmniejszenia wielkości fragmentów w wyniku tej restrykcji. Kliknij tutaj, aby zobaczyć powiększoną wersję tej ryciny.
Ryciny w niniejszym manuskrypcie zostały wygenerowane z osobnego, powtórzonego eksperymentu, innego niż ten opublikowany wcześniej przez Lafontaine i współpracowników21. Po uzyskaniu danych z sekwencjonowania wysokoprzepustowego, do przetworzenia danych Hi-C wykorzystano Open Chromatin Collective (Open2C: https://github.com/open2c). Podobny rurociąg przetwarzania można znaleźć w portalu danych projektu 4D Nucleome (https://data.4dnucleome.org/resources/data-analysis/hi_c-processing-pipeline). W skrócie, zaimplementowano rurociąg Nextflow distiller (https://github.com/open2c/distiller-nf), aby (1) dopasować sekwencje cząsteczek Hi-C do genomu referencyjnego, (2) sparsować pliki dopasowań .sam i utworzyć pliki z parami Hi-C, (3) odfiltrować duplikaty PCR oraz (4) zagregować pary w macierzach binowych interakcji Hi-C. Macierze te w formacie HDF5, nazywane coolerami, mogą być następnie (1) przeglądane na serwerze HiGlass (https://higlass.io/) oraz (2) analizowane przy użyciu szerokiego zestawu otwartych narzędzi obliczeniowych zawartych w kolekcji „cooltools” utrzymywanej przez Open Chromatin Collective (https://github.com/open2c/cooltools), w celu wyodrębnienia i ilościowego określenia cech fałdowania, takich jak przedziały (compartments), TADs i pętle.
Niektóre wskaźniki jakości bibliotek Hi-C3.0 można ocenić bezpośrednio po zmapowaniu par odczytów do genomu referencyjnego, korzystając z kilku prostych metryk/wskaźników. Po pierwsze, w przypadku komórek ludzkich zazwyczaj ok. 50% zsekwencjonowanych par odczytów można zmapować jednoznacznie. Ze względu na polimerową naturę chromosomów, większość tych zmapowanych odczytów (~60%-90%) reprezentuje oddziaływania w obrębie jednego chromosomu (cis), przy czym częstotliwość oddziaływań gwałtownie spada wraz ze wzrostem odległości genomicznej (spadek zależny od odległości). Spadek zależny od odległości najlepiej zobrazować na „wykresie skalowania” (scaling plot), który przedstawia prawdopodobieństwo kontaktu (na ramię chromosomu) jako funkcję odległości genomicznej. Stwierdzono, że zastosowanie różnych czynników sieciujących i enzymów może zmienić charakter spadku zależnego od odległości w zakresach długich i krótkich dystansów17. Dodanie sieciowania DSG zwiększa wykrywalność oddziaływań na krótkich dystansach w połączeniu z takimi enzymami jak Mnase oraz kombinacjami DpnII-DdeI, które generują mniejsze fragmenty (Rycina 6A).
Zależny od odległości spadek sygnału można również zaobserwować bezpośrednio na dwuwymiarowych macierzach oddziaływań: oddziaływania stają się rzadsze wraz ze zwiększaniem odległości od głównej przekątnej (Rysunek 6B). Ponadto cechy fałdowania genomu, takie jak przedziały (compartments), domeny TAD oraz pętle, można zidentyfikować na macierzach Hi-C i wykresach skalowania jako odchylenia od ogólnego średniego dla całego genomu spadku zależnego od odległości. Co istotne, sieciowanie za pomocą DSG w połączeniu z FA zmniejsza liczbę ligacji losowych, które ze względu na polimerową naturę chromosomów są nieograniczone i tym samym częściej występują między chromosomami (in trans) (Rysunek 6C). Ograniczenie ligacji losowych prowadzi do zwiększenia stosunku sygnału do szumu, szczególnie w przypadku oddziaływań międzychromosomalnych oraz bardzo długozasięgowych (>10-50 Mb) w obrębie jednego chromosomu.

Rysunek 6: Reprezentatywne wyniki zmapowanych i przefiltrowanych bibliotek Hi-C. (A) Wykresy skalowania z prawdopodobieństwem kontaktu i jego pochodną dla różnych enzymów, uporządkowane według długości fragmentu (góra) oraz sieciowania za pomocą FA lub FA + DSG (dół). Trawienie przy użyciu MNAse (microC) lub DpnII-DdeI (Hi-C 3.0) znacząco zwiększa liczbę kontaktów krótkozasięgowych (góra), podobnie jak dodanie DSG do FA (dół). (B) Kolumny przedstawiają mapy ciepła Hi-C dla trawienia DpnII po samym sieciowaniu FA oraz trawienia DpnII lub DdeI po sieciowaniu FA+DSG. Białe strzałki wskazują zwiększającą się intensywność „punktów” po sieciowaniu DSG i trawieniu DdeI, co sugeruje lepszą detekcję pętli DNA. Wiersze przedstawiają różne fragmenty chromosomu 3 przy zwiększającej się rozdzielczości, zgodnie z panelem C: górny wiersz: cały chromosom 3 (0-198 295 559 Mb); środkowy wiersz: 186-196 Mb; dolny wiersz: 191,0-191,5 Mb. (C) Wykresy pokrycia dla regionów przedstawionych w części A. Czarne strzałki wskazują niższe pokrycie (% odczytów cis) dla sieciowania samym FA. Skróty: FA = formaldehyd; DSG = glutarate disuccinimidyl; chr = chromosom. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Nie wszystkie zmapowane odczyty są użyteczne. Drugim wskaźnikiem jakości jest liczba duplikatów PCR. Jest wysoce nieprawdopodobne, aby identyczne odczyty wystąpiły przypadkowo po ligacji i sonikacji. Zatem takie odczyty prawdopodobnie są wynikiem amplifikacji PCR i muszą zostać odfiltrowane. Duplikaty często powstają, gdy do amplifikacji bibliotek o niskiej złożoności wymagana jest zbyt duża liczba cykli PCR. Zazwyczaj w przypadku Hi-C większość bibliotek wymaga jedynie 5-8 cykli końcowej amplifikacji PCR, co określa się poprzez PCR miareczkowy (patrz krok 3.9; Rycina 4). Niemniej jednak biblioteki o wystarczającej złożoności można uzyskać nawet po 14 cyklach amplifikacji PCR.
Kolejna kategoria duplikatów odczytów, tzw. duplikaty optyczne, może wynikać z procesu amplifikacji na platformach sekwencjonowania Illumina wykorzystujących wzorcowane ogniwa przepływowe (takich jak HiSeq4000). Duplikaty optyczne powstają albo w wyniku przeładowania ogniwa przepływowego, co powoduje, że (duże) klastry są rozpoznawane jako dwa oddzielne klastry, albo w wyniku lokalnej reklasteryzacji oryginalnej cząsteczki typu paired-end po pierwszej rundzie PCR. Ponieważ oba typy duplikatyw optycznych mają charakter lokalny, można je zidentyfikować i odróżnić od duplikatów PCR na podstawie ich położenia w ogniwie przepływowym. Podczas gdy biblioteki z zawartością >15% duplikatów PCR wymagałyby regeneracji, biblioteki z duplikatami optycznymi mogą zostać ponownie załadowane po zoptymalizowaniu procesu ładowania.
Tabela uzupełniająca S1: Bufory i roztwory. Proszę kliknąć tutaj, aby pobrać tę tabelę.
Krytyczne etapy obsługi komórek
Chociaż możliwe jest użycie mniejszej liczby komórek wejściowych, protokół ten został zoptymalizowany dla ~5 × 106 komórek na linię sekwencjonowania (~400 mln odczytów), aby zapewnić odpowiednią złożoność po głębokim sekwencjonowaniu. Komórki najlepiej liczyć przed utrwaleniem. W przypadku generowania bibliotek ultragłębokich zazwyczaj mnożymy liczbę linii (i komórek), aż do osiągnięcia żądanej głębokości odczytu. W celu optymalnego utrwalenia, pożywkę zawierającą surowicę należy zastąpić PBS przed utrwaleniem FA, a roztwory utrwalające należy dodać natychmiast i bez gradientów stężeń15,22. W przypadku zbierania komórek skrobanie jest preferowane zamiast trypsynizacji, ponieważ przejście od bardziej płaskiego do kulistego kształtu po trypsynizacji może wpływać na konformację jądrową. Po dodaniu DSG luźne i grudkowate granulki komórkowe łatwo się gubią. Na tym etapie należy zachować ostrożność podczas obchodzenia się z komórkami i dodać do 0,05% BSA, aby zmniejszyć zbrylanie.
Modyfikacje metody
Protokół ten został opracowany przy użyciu komórek ludzkich17. Jednak w oparciu o doświadczenie w wychwytywaniu konformacji chromosomów, protokół ten powinien działać w przypadku większości komórek eukariotycznych. W przypadku znacznie niższego wkładu (~1 × 106 komórek) zalecamy użycie połowy objętości do procedur wychwytywania lizy i konformacji [kroki 2.1-2.4]. Pozwoliłoby to również na przeprowadzenie izolacji DNA [krok 2.5] w wirówce stołowej z probówkami o pojemności 1,7 ml, co mogłoby poprawić granulowanie przy niskich stężeniach DNA. Kwantyfikacja DNA (krok 2.6.6) wskaże, jak należy postępować. W przypadku małych ilości wyizolowanego DNA (1-5 μg) sugerujemy pominięcie wyboru rozmiaru (krok 3.3) i przystąpienie do usuwania biotyny po zmniejszeniu objętości ze 130 μL do ~45 μL za pomocą CFU.
Protokół ten został opracowany specjalnie w celu zapewnienia wysokiej jakości danych po późniejszym sieciowaniu z FA i DSG oraz trawieniu z DpnII i DdeI. Jednak alternatywne strategie sieciowania, takie jak FA, a następnie EGS (bursztynian glikolu etylenowego)), który jest również stosowany w ChIP-seq23 i ChIA-PET24, mogą działać równie dobrze17. Podobnie do trawienia można stosować różne kombinacje enzymów, takie jak DpnII i HinfI18 lub MboI, MseI i NlaIII19 . Dostosowując kombinacje enzymów, należy pamiętać o stosowaniu biotynylowanych nukleotydów, które mogą wypełnić określone 5' zwisy i używać najbardziej optymalnych dla każdego koktajlu. DpnII jest dostarczany z własnym buforem, a producent enzymu zaleca określony bufor do trawienia DdeI. Jednak w przypadku podwójnego trawienia z DpnII i DdeI w tym protokole zalecany jest bufor restrykcyjny, ponieważ jest on oceniany na 100% aktywność dla obu enzymów.
Rozwiązywanie problemów z przechwytywaniem konformacji
Trzy kluczowe etapy wychwytywania konformacji chromosomów: sieciowanie, trawienie i religacja zostały wykonane, zanim wyniki będzie można zwizualizować na żelu. Aby określić jakość każdego z tych trzech etapów i określić, gdzie mogły wystąpić problemy, pobiera się podwielokrotności przed (CI) i po trawieniu (DC) i umieszcza je w żelu wraz z ligowaną próbką Hi-C (ryc. 2). Żel ten służy do określenia jakości próbki Hi-C i tego, czy warto będzie kontynuować protokół. Bez CI i DC trudno jest wskazać potencjalne nieoptymalne kroki. Warto zauważyć, że nieoptymalne podwiązanie może być spowodowane problemem w samym podwiązaniu, wypełnieniu lub problemem z sieciowaniem. Aby rozwiązać problemy z sieciowaniem, pamiętaj, aby nie używać więcej niż 1 × 107 komórek na bibliotekę i zacznij od świeżych odczynników do sieciowania i czystych komórek (tj. przepłukanych PBS). W przypadku podwiązania upewnij się, że komórki i mieszanina ligacji są trzymane na lodzie. Dodać ligazę DNA T4 tuż przed 4 godzinną inkubacją w temperaturze 16 °C i dobrze wymieszać.
Rozwiązywanie problemów z przygotowywaniem biblioteki
Jeśli potrzebnych jest więcej niż 10 cykli PCR lub po miareczkowaniu PCR na żelu nie widać produktu PCR (ryc. 4), istnieje kilka opcji zapisania próbki Hi-C. Wracając do miareczkowania PCR, pierwszą opcją jest ponowna próba PCR. Jeśli nadal nie ma wystarczającej ilości produktu, można spróbować kolejnej rundy ogona A i podwiązania adaptera (krok 3.6) po dwukrotnym umyciu kulek buforem 1x TLE. Po tym dodatkowym ogonie A i podwiązaniu adapterem można przejść do miareczkowania PCR jak poprzednio. Jeśli nadal nie ma produktu, ostatnią opcją jest ponowne rozpoznanie ułamka 0,8x z kroku 3.3 i kontynuowanie od tego miejsca.
Ograniczenia i zalety Hi-C3.0
Ważne jest, aby zdać sobie sprawę, że Hi-C jest metodą opartą na populacji, która rejestruje średnią częstotliwość interakcji między parami loci w populacji komórek. Niektóre analizy obliczeniowe mają na celu rozwikłanie kombinacji konformacji z populacji25, ale w zasadzie Hi-C jest ślepy na różnice między komórkami. Chociaż możliwe jest wykonanie pojedynczej komórki Hi-C26,27 i można dokonać wnioskowania obliczeniowego28, jednokomórkowa Hi-C nie nadaje się do uzyskiwania informacji 3C o ultrawysokiej rozdzielczości. Dodatkowym ograniczeniem Hi-C jest to, że wykrywa on tylko oddziaływania parami. Do wykrywania oddziaływań wielokontaktowych można użyć częstych obcinaczy w połączeniu z sekwencjonowaniem krótkiego odczytu (Illumina)16 lub wykonać wielokontaktowe 3C29 lub 4C30, używając sekwencjonowania długiego odczytu z platform PacBio lub Oxford Nanopore. Opracowano również pochodne Hi-C do specyficznego wykrywania kontaktów między chromatydami siostrzanymi i wzdłuż nich31,32.
Chociaż Hi-C19 i Micro-C33 mogą być używane do generowania map kontaktowych w rozdzielczościach subkilobazowych, oba wymagają dużej ilości odczytów sekwencjonowania, co może stać się kosztownym przedsięwzięciem. Aby uzyskać podobną, a nawet wyższą rozdzielczość bez ponoszenia kosztów, można zastosować wzbogacanie o określone regiony genomu (capture-C34) lub specyficzne interakcje białkowe (ChiA-PET35, PLAC-seq36, Hi-ChIP37). Zaletą i wadą tych aplikacji wzbogacających jest to, że próbkowana jest tylko ograniczona liczba interakcji. Przy takich wzbogaceniach zatraca się globalny aspekt Hi-C (i możliwość globalnej normalizacji).
Znaczenie i potencjalne zastosowania Hi-C3.0
Protokół ten został zaprojektowany, aby umożliwić uzyskanie wysokiej rozdzielczości, ultragłębokiego 3C, przy jednoczesnym wykrywaniu elementów składania na dużą skalę, takich jak TAD i komory17 (ilustracja 6). Protokół ten zaczyna się od 5 × 106 komórek na probówkę dla każdej biblioteki Hi-C, co powinno wystarczyć do sekwencjonowania jednej lub dwóch linii w komórce przepływowej, aby uzyskać do 1 miliarda odczytów sparowanych końców. Do ultragłębokiego sekwencjonowania należy przygotować wiele probówek po 5 × 106 komórek, w zależności od liczby zmapowanych odczytów i duplikatów PCR. Przy najwyższej rozdzielczości (<1 kb) interakcje pętlowe występują głównie między miejscami CTCF, ale można również wykryć interakcje promotor-wzmacniacz. Czytelnicy mogą zapoznać się z Akgol Oksuz et al.17 w celu uzyskania szczegółowego opisu analizy danych.
Autorzy nie mają do ujawnienia żadnych konfliktów interesów.
Chcielibyśmy podziękować Denisowi Lafontaine'owi za opracowanie protokołu i Sergeyowi Venevowi za pomoc bioinformatyczną. Praca ta została wsparta grantem z National Institutes of Health Common Fund 4D Nucleome Program dla J.D. (U54-DK107980, UM1-HG011536). J.D. jest badaczem w Instytucie Medycznym Howarda Hughesa.
Ten artykuł podlega polityce HHMI dotyczącej otwartego dostępu do publikacji. Szefowie laboratoriów HHMI wcześniej udzielali niewyłącznej licencji CC BY 4.0 dla publiczności i licencji podlegającej sublicencjonowaniu dla HHMI w swoich artykułach badawczych. Zgodnie z tymi licencjami, zaakceptowany przez autora manuskrypt tego artykułu może być bezpłatnie udostępniony na licencji CC BY 4.0 natychmiast po publikacji.
| Nazwa | Firma | Numer katalogowy | Komentarze |
|---|---|---|---|
| 1 kb Drabina | New England Biolabs | N3232L | |
| Agarose | Invitrogen | 16500100 | |
| Koraliki magnetyczne Agencourt AMPure XP, 60 ml | Beckman Coulter | A63881 | |
| Filtr odśrodkowy Amicon Ultra-0,5 (CFU) | EMD Milipore | UFC500396 | |
| bufor do wyżarzania (5x) | Zobacz przepis w materiałach | ||
| ATP 10 mM | ThermoFisher | R0441 | |
| Avanti J-25i Ultrawirówka z chłodzeniem o dużej prędkości Beckman | Coulter | ||
| Probówka ultrawirówkowa Beckman Coulter 35 ml | Beckman Coulter | 357002 | |
| Bufor wiążący (2x) | Zobacz przepis w materiałach | ||
| biotyna-14-dATP 0,4 mM | Invitrogen | 19524-016 | |
| BSA 10 mg/ml | New England Biolabs | B9000S | rozcieńczony od 20 mg/ml |
| Skrobak do komórek | Falcon | 353089 | |
| Skrobak do komórek | Corning | 3008 | |
| Stożkowe probówki polipropylenowe 50 ml | Denville | C1062-P | |
| Rurka stożkowa 15 ml | Denville | C1017-P | |
| Mikrorurka Covaris AFA z włóknem AFA Zatrzask 130 &MICRO; L | Covaris | 520045/520077 | |
| Covaris Sonicator | Covaris | E220/E220evolution/M220 | |
| Kolba hodowlana 175 cm2 | Falcon | 353112 | |
| Płytki hodowlane 150 mm x 25 mm | Corning | 430599 | |
| dATP 1 | mM Invitrogen | 56172 | |
| dCTP 10 mM | Invitrogen | 56173 | |
| DdeI | New England Biolabs | R0175L | |
| dGTP 10 mM | Invitrogen | 56174 | |
| DMSO | Sigma | D2650-5x10ML | |
| dNTP mix 25 mM | Homogenizator Invitrogen | 10297117 | |
| Dounce | DWK Life Sciences | 8853010002/8853030002 | |
| DPBS | Gibco | 14190-144 | |
| DPNII | New England Biolabs | R0543M | |
| DSG | ThermoScientific | 20593 | |
| dTTP 10 mM | Invitrogen | 56175 | |
| Etanol 70% | Fisher | A409-4 | Rozcieńczony ze 100% |
| bromku etydyny | Fisher | BP1302-10 | |
| Formaldehyd (37%) | Fisher | BP531-500 | |
| Barwnik ładujący w żelu (6x) | New England Biolabs | B7024S | |
| Glicyna w ultraczystej wodzie 2,5 M | Sigma | G8898-1KG | |
| HBSS | Gibco | 14025-092 | |
| Detergent Igepal CA-630 | MP Biomedicals | 198596 | |
| Klenow Polimeraza DNA 5 U/&mikro; L | New England Biolabs | M0210L | |
| Fragment Klenowa 3-->5' egzo-, 5 U/&mikro; L | New England Biolabs | M0212L | |
| (10x) | New England Biolabs | B7203S | |
| Probówka wirówkowa LoBind | |||
| Drabinka DNA o niskiej masie cząsteczkowej | England Biolabs | N3233L | |
| materiały | |||
| Magnetyczny separator cząstek | ThermoFisher | 12321D | |
| Probówki do mikrofuge 1,7 ml | Axygen | MCT-175-C | |
| MyOne Streptawidyna C1 kulki | Invitrogen | 65001 | |
| NEBuffer 2.1 (10x) | New England Biolabs | B7002S | |
| NEBuffer 3.1 (10x) | New England Biolabs | B7203S | |
| PBS | Gibco | 70013-032 | |
| Probówki do PCR (paskowe) Biorad | TBS0201/ TCS0803 | ||
| Termocykler do PCR | Biorad | T100 | |
| Pfu Ultra II Bufor (10x) | Agilent | W zestawie z | |
| polimerazą DNA Pfu Ultra PfuUltra II Fusion HS | Agilent | 600674 | |
| Rurka z blokadą fazową 15 ml | Qiagen | 129065 | |
| Rurki z blokadą fazową 2 ml | Qiagen | 129056 | |
| Fenol: chloroform: alkohol izoamylowy | Invitrogen | 15593-049 | |
| Koktajl inhibitorów proteazy | ThermoFisher | 78440 | |
| Proteinaza K w wodzie ultraczystej 10 mg/ml | Invitrogen | 25530-031 | |
| Wirówka z chłodzeniem | Eppendorf | 5810R | |
| RNaza A, bez DNaz i proteaz 10 mg/ml | Thermo Scientific | EN0531 | |
| Rotator | Argos technologies | EW-04397-40 lub platforma | |
| SDS 1% | Fisher | BP13111 | |
| Octan sodu pH = 5,2, 3 M | Sigma | ||
| Sub-Cell GT System Elektroforezy Poziomej | Biorad | 1704401 | |
| Ligaza DNA T4 1 U/&mikro; L | Invitrogen | 100004817 | |
| Polimeraza DNA | T4 New England Biolabs | M0203L | |
| Polimeraza DNA T4 | New England Biolabs | M0203L | |
| Polimeraza DNA T4 3 U/&mikro; L | Bufor ligacyjny New England Biolabs | M0203L | |
| T4 (5x) | Invitrogen | Y90001 | |
| Kinaza polinukleotydowa T4 10 U/&mikro; L | New England Biolabs | M0201L | |
| Wirówka stołowa | Eppendorf | 5425 | |
| Bufor TBE | Zobacz przepis w materiałach | ||
| Tris Low EDTA Buffer (TLE) | Zobacz przepis w materiałach | ||
| Triton X-100 (10%) | Sigma | 93443 | |
| Adapter Truseq oligos | Zintegrowane DNA Technologie (IDT) | )https://www.idtdna.com/site/order/oligoentry | 250 nmol i oczyszczony HPLC |
| detergent Tween 20 | Fisher | 9005-64-5 | |
| Bufor do płukania Tween | Zobacz przepis w materiałach | ||
| Vortex | Scientific Industries | (G560) SI-0236 |