$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Profile trawienia chromatyny
Optymalizacja trawienia MNazy jest niezbędna dla powodzenia tego protokołu. Kluczowe znaczenie ma wygenerowanie profilu trawienia zdominowanego przez pojedyncze rozmiary fragmentów nukleosomów, a jednocześnie nie trawionego w nadmiarze, aby umożliwić odzyskanie fragmentów nukleosomów wyższego rzędu. Idealny profil trawienia składa się z większości pojedynczych fragmentów nukleosomów, przy czym niewielka frakcja reprezentuje fragmenty mniejsze i większe niż pojedyncze nukleosomy. Rysunek 1 pokazuje przykłady idealnych, nadmiernie trawionych i niedostatecznie trawionych profili rozkładu wielkości. Należy zauważyć, że nieoptymalne trawienie chromatyny będzie również widoczne w profilu biblioteki sekwencjonowania wygenerowanej z materiału IP (Rysunek 2).
Walidacja jakości biblioteki ndChIP-seq przez qPCR
qPCR to sprawdzona metoda oceny jakości ChIP18,19,20. Podczas wykonywania sekwencji ndChIP na 10 000 komórek wydajność kwasu nukleinowego po IP będzie poniżej 1 ng. Dlatego ważne jest, aby przeprowadzić qPCR po zbudowaniu biblioteki, aby ocenić względne wzbogacenie regionów docelowych w porównaniu z tłem. Aby zapewnić oszacowanie tła, generowane są biblioteki zbudowane z chromatyny trawionej przez MNazę (dane wejściowe). Dla każdej biblioteki IP potrzebne są dwa zestawy starterów (patrz Tabela uzupełniająca3, aby uzyskać listę starterów dla powszechnie stosowanych znaczników histonowych). Jeden zestaw starterów powinien być specyficzny dla regionu genomu, który jest konsekwentnie związany z modyfikacją histonów będącą przedmiotem zainteresowania (cel dodatni) i innego regionu, który nie jest oznaczony modyfikacją histonów będącą przedmiotem zainteresowania (cel ujemny). Jakość biblioteki ChIP-seq będzie oceniana jako wzbogacenie fałdy w odniesieniu do biblioteki wejściowej. Wzbogacenie krotności można obliczyć za pomocą następującego równania, które zakłada wykładniczą amplifikację docelowego regionu genomowego: 2Ctwejście - CtIP. Nasz niestandardowy pakiet oprogramowania statystycznego R, qcQpcr_v1.2, jest odpowiedni do analizy wzbogacania qPCR natywnych bibliotek ChIP-seq (Supplemental Code Files) o niskim nakładzie wejściowym. Rysunek 3 przedstawia wynik qPCR dla udanych i nieudanych bibliotek ChIP-seq. Minimalna oczekiwana wartość wzbogacenia składania dla dobrej jakości bibliotek ndChIP-seq wynosi 16 dla wąskich znaczników, takich jak H3K4me3, i 7 dla szerokich znaczników, na przykład H3K27me3.
Modelowanie dostępności MNase
Analiza obliczeniowa sekwencji ChIP jest złożona i unikalna dla każdego środowiska eksperymentalnego. Zestaw wytycznych opracowanych przez Międzynarodowe Konsorcjum Epigenomiczne Człowieka (IHEC) i Encyklopedię Elementów DNA (ENCODE) można wykorzystać do oceny jakości bibliotek ChIP-seq21. Ważne jest, aby pamiętać, że głębokość sekwencjonowania bibliotek ma wpływ na wykrywanie i rozpoznawanie wzbogaconych regionów20. Liczba wykrytych pików wzrasta i zbliża się do plateau wraz ze wzrostem głębokości odczytu. Zalecamy, aby biblioteki ndChIP-seq były sekwencjonowane zgodnie z zaleceniami IHEC dotyczącymi 50 milionów sparowanych odczytów (25 milionów fragmentów) dla wąskich znaczników (np. H3K4me3) i 100 milionów sparowanych odczytów (50 milionów fragmentów) dla szerokich znaczników (np. H3K27me3) i input22. Te głębokości sekwencjonowania zapewniają wystarczające wyrównanie sekwencji do wykrywania najważniejszych pików przy użyciu powszechnie używanych wywołujących piki ChIP-seq, takich jak MACS2 i HOMER, bez osiągania nasycenia23,24. Wysokiej jakości biblioteka ndChIP-seq ssaków ma wskaźnik duplikacji PCR <10% i wskaźnik wyrównania genomu referencyjnego > 90% (w tym zduplikowane odczyty). Udane biblioteki ndChIP-seq będą zawierały wysoce skorelowane repliki, przy czym znaczna część (> 40%) wyrównanych odczytów w MACS222 zidentyfikowanych wzbogaconych pików, a inspekcja wyrównanych odczytów w przeglądarce genomu powinna ujawnić wizualnie wykrywalne wzbogacenia w porównaniu z biblioteką wejściową (Rysunek 4). Ponadto ndChIP-seq może być używany do oceny gęstości nukleosomów przy użyciu algorytmu rozkładu mieszaniny Gaussa (w1 * n(x; μ1,σ 1) + w2 * n(x; μ2,σ 2) = 1) w MACS2 zidentyfikowano wzbogacone regiony do modelowania gęstości nukleosomów zdefiniowanej przez granice dostępne dla MNazy. W tym modelu w1 oznacza masę dystrybucji mononukleosomów, a w2 oznacza masę dystrybucji dinukleosomów. Gdzie w1 jest większe niż w2, występuje dominacja fragmentów mononukleosomów i odwrotnie. Ta analiza wymaga, aby biblioteki były sekwencjonowane w sposób sparowany, aby można było zdefiniować rozmiary fragmentów. Aby zastosować algorytm rozkładu mieszaniny Gaussa, najpierw identyfikuje się statystycznie istotne regiony wzbogacone. Sugerujemy wywołanie wartości szczytowej za pomocą MACS2 przy użyciu Input jako kontrolki i z domyślnymi ustawieniami dla wąskich znaków i odcięciem wartości q wynoszącym 0,01 dla szerokich znaków. Szereg pakietów statystycznych wykorzystujących algorytm rozkładu mieszaniny Gaussa jest dostępnych w powszechnie stosowanych pakietach oprogramowania statystycznego. Wykorzystując średni rozmiar fragmentu, określony przez granice odczytu sparowanych końców próbek IPed, rozkłady w MACS2 zidentyfikowały wzbogacone promotory, a algorytm rozkładu mieszaniny Gaussa można zastosować do każdego promotora przy użyciu pakietu R-statistical Mclust w wersji 3.025 w celu obliczenia rozkładu ważonego. W tym zastosowaniu zalecamy wyeliminowanie promotorów zawierających mniej niż 30 wyrównanych fragmentów, ponieważ poniżej tego progu wynikowe szacunki masy stają się niewiarygodne. Dobrej jakości biblioteka ndChIP-seq generuje rozkład mieszaniny Gaussa, który składa się z dwóch głównych składników o średnich wartościach odpowiadających długościom fragmentów mono-, dinukleosomów.

Rysunek 1: Ocena trawienia MNazy przed wygenerowaniem biblioteki. Profile analizatora elektroforezy kapilarnej opartej na chipach optymalnej chromatyny trawionej (A), niedostatecznie trawionej (B) i nadmiernie trawionej (C) chromatyny. Kontrpróby biologiczne są pokazane jako niebieskie, czerwone i zielone ślady. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 2: Ocena trawienia MNazy po wygenerowaniu biblioteki. (A) Profile konstrukcji po bibliotece optymalnie strawionych danych wejściowych (kontrpróby biologiczne; czerwony, zielony,) i IP (repliki biologiczne; cyjan, fioletowy, niebieski) oraz (B) nieoptymalnych danych wejściowych (repliki biologiczne; czerwony, zielony, niebieski) i IP (repliki biologiczne; cyjan, fioletowy, pomarańczowy) bibliotek. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 3: Ilościowy PCR po budowie biblioteki może być użyty do oceny jakości bibliotek ndChIP-seq. Wzbogacenie w krotność bibliotek IP H3K4me3 w odniesieniu do bibliotek wejściowych oblicza się jako 2 (Ct danych wejściowych - Ct adresu IP) dla celów dodatnich i ujemnych przy użyciu qcQpcr_v1.2. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 4: Reprezentatywna biblioteka ndChIP-seq zbudowana z 10 000 pierwotnych komórek krwi pępowinowej CD34+. Korelacja Pearsona sygnału H3K4me3 (odczyty na milion zmapowanych odczytów) obliczona w promotorach (TSS+/-2Kb) między 3 kontrpróbami biologicznymi, (A) replikatem 1 i 2, (B) powtórzeniem 1 i 3, (C) powtórzeniem 2 i 3. (D) Widok przeglądarki UCSC na klaster genów HOXA usieciowanych ChIP-seq wygenerowanych z 1 miliona komórek na IP, udany ndChIP-seq z 10 000 komórek na IP i nieudany ndChIP-seq z 10 000 komórek na IP. (czerwony: H3K27me3, zielony: H3K4me3 i: Wejście). (E) Frakcja zmapowanych odczytów w MACS2 zidentyfikowała wzbogacone regiony H3K4me3 () i H3K27me3 (szary). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
Skład bufora
A.1. Bufor immunoprecypitacyjny (IP)
20 mM Tris-HCl pH 7,5
2 mM EDTA
150 mM NaCl
0,1% Tryton X-100
0,1% Deoksycholan
10 mM maślanu sodu
A.2. Bufor o niskiej zawartości soli
20 mM Tris-HCl pH 8,0
2 mM EDTA
150 mM NaCl
1% Tryton X-100
0,1% SDS
A.3. Bufor do przemywania o wysokiej zawartości soli
20 mM Tris-HCl pH 8,0
2 mM EDTA
500 mM NaCl
1% Tryton X-100
0,1% SDS
A.4. Bufor elucyjny ChIP
100 mM NaHCO
3
1% SDS
A.5. 1x Bufor do lizy – 1 ml
0,1% Tryton
0,1% Deoksycholan
10 mM maślanu sodu
A.6. Bufor do rozcieńczania Ab
0,05% (w/v) Azydek
0,05% środki przeciwdrobnoustrojowe o szerokim spektrum działania (
np. ProClin 300)
w PBS
A.7. 30% roztwór kulek magnetycznych PEG/1M NaCl (odniesienie16)
30% PEG
1 M NaCl
10 mM Tris HCl pH 7,5
1 mM EDTA
1 ml umytych kulek superparamagnetycznych
A.8. 20% roztwór kulek magnetycznych PEG/1M NaCl (odniesienie16)
30% PEG
1 M NaCl
10 mM Tris HCl pH 7,5
1 mM EDTA
1 ml umytych kulek superparamagnetycznych
Tabela 1: kompozycja bufora ndChIP-seq.
Modyfikacja histonu
Stężenie (μg/μL)
H3K4me3 powiedział:
0,125 pkt.
H3K4me1 powiedział:
0,25 pkt.
H3K27me3 powiedział:
0,125 pkt.
H3K9me3 powiedział:
0,125 pkt.
H3K36me3 powiedział:
0,125 pkt.
H3K27ac powiedział:
0,125 pkt.
Tabela 2: Ilość przeciwciał wymagana dla ndChIP-seq.
powiedział:
szt.
szt.
| Reganet | Objętość (μL) |
| Ultra czysta woda | Z numerem 478 |
| 1 M Tris-HCl pH 7,5 | 10 |
| 0,5 mln EDTA | 10 |
| 5 M NaCl | cyfra arabska |
| glicerol | 500 |
| Całkowita objętość | 1 000 |
Tabela 3: Przepis na bufor rozcieńczający MNazę.
Odczynnik
Objętość (μL)
Ultra czysta woda
13
Naziemna telewizja cyfrowa 20 mM
1
10x Bufor MNazy
4
20 U/μl Mnase
cyfra arabska
Całkowita objętość
20
Tabela 4: Przepis na MNase Master Mix.
Odczynnik
Objętość (μL)
Bufor elucyjny
Rozdział 30
Bufor G2
8
Proteazy
cyfra arabska
Całkowita objętość
Rozdział 40
Tabela 5: Przepis na Mistrza Oczyszczania DNA.
Odczynnik
Objętość (μL)
Ultra czysta woda
3.3
10x represyjny bufor endonukleazy (
np. NEBuffer)
5
25 mM ATP
cyfra arabska
10 mM dNTP
cyfra arabska
Kinaza polinukleotydowa T4 (10 U/μl)
1
Polimeraza DNA T4 (3 U/μl)
1,5
Polimeraza DNA I, duży (Klenow) fragment (5 U/μl)
0,2 pkt.
Całkowita objętość
15
Tabela 6: Przepis na Mistrz Naprawy Końcowej.
Odczynnik
Objętość (μL)
Ultra czysta woda
8
10x represyjny bufor endonukleazy (
np. NEBuffer)
5
10 mM dATP
1
Klenow (3'-5' egzo-)
1
Całkowita objętość
15
Tabela 7: Przepis na A-Tailing Master Mix.
Odczynnik
Objętość (μL)
Ultra czysta woda
4.3
5x bufor do szybkiego podwiązywania
12
Ligaza DNA T4 (2000 U/μl)
6.7
Całkowita objętość
23
Tabela 8: Przepis na Mieszankę Główną Podwiązania Adaptera.
Odczynnik
Objętość (μL)
Ultra czysta woda
7
Starter do PCR 25 uM 1.0
cyfra arabska
5x bufor HF
12
Dmso
1,5
Polimeraza DNA
0,5
Całkowita objętość
23
Tabela 9: Przepis na PCR Master Mix.
Temperatura (°C)
Czas trwania (s)
Liczba cykli
Rozdział 98
60 Rozdział
1
Rozdział 98
Rozdział 30
65 Rozdział
15
10
72 Rozdział 72
15
72 Rozdział 72
300 szt.
1
4
trzymać
trzymać
Tabela 10: Metoda uruchamiania PCR.
Oligo (Oligo)
kolejność
Modyfikacji
PE_adapter 1
5'- /5Phos/GAT CGG AAG AGC GGT TCA GCA GGA ATG CCG AG -3'
Modyfikacja 5': fosforylacja
PE_adapter 2
5' - ACA CTC TTT CCC TAC ACG ACG CTC TTC CGA CGA TC*T - 3'
Modyfikacja 3': *T jest wiązaniem fosforotionowym
Tabela uzupełniająca 1: Sekwencje oligologiczne do generowania adaptera PE.
Nazwa podkładu
Sekwencja
indeks
IndexRevC (do użycia do sekwencjonowania)
Starter do odwrotnego indeksowania PCR 1
CAAGCAGAAGACGGCATACGAGATCGTGATCGGTCTCTCCTCTGCTGCTGAACCGCTCTTCCGATCT
CGTGAT (Organizacja Zdrowia)
ATCACG (Biblioteka ATCACG)
Starter do odwróconego indeksowania PCR 2
CAAGCAGAAGACGGCATACGAGATCTGATCCGGTCTCTCTCCTGCTGCTGAACCGCTCTTCCGATCT
Protokół CTGATC (CTGATC)
gatcag powiedział:
Starter do odwróconego indeksowania PCR 3
CAAGCAGAAGACGGCATACGAGATGGGGTCTCTCTCTCCTGCTCTGCTGCTGCTCTCTCTCTCTCTCTCTCTCTCTCT
GGGGTT (Język gggtt)
Biblioteka AACCCC
Starter do odwróconego indeksowania PCR 4
CAAGCAGAAGACGGCATACGAGATGGGTCGGTCTCTCTCCTGCTGCTGAACCGCTCTTCCGATCT
CTGGGT powiedział:
Acccag
Podkład do odwróconego indeksowania PCR 5
CAAGCAGAAGACGGCATACGAGATAGCGCTCGGTCTCTCCTGCTGCTGCTGCTGCTGCTGCTCTCTCTCTCTCTCTCTCTCT
AGCGCT powiedział:
agcgct powiedział:
Starter do odwrotnego indeksowania PCR 6
CAAGCAGAAGACGGCATACGAGATCTTGCGGTCTCTCTCCTGCTGCTGCTGAACCGCTCTTCCGATCT
Protokół CTTTTG (Mechanizm Monitorowania Terroryzmu
CAAAAG powiedział:
Starter do odwrotnego indeksowania PCR 7
CAAGCAGAAGACGGCATACGAGATTGTTGGTCTCTCTCCTGCTGCTGCTGAACCGCTCTTCCGATCT
TGTTGG powiedział:
CCAACA (Ccaaca)
Podkład do odwrotnego indeksowania PCR 8
CAAGCAGAAGACGGCATACGAGATAGCTAGCGGTCTCTCCTCTGCTGCTGAACCGCTCTTCCGATCT
Znacznik AGCTAG
ctagct powiedział:
Podkład do odwróconego indeksowania PCR 9
CAAGCAGAAGACGGCATACGAGATAGCATCCGGTCTCTCCTCTGCTGCTGAACCGCTCTTCCGATCT
AGCATC (Organizacja AGCATC)
Biblioteka GATGCT
Podkład do odwrotnego indeksowania PCR 10
CAAGCAGAAGACGGCATACGAGATCGATTCGGTCTCTCCTGCTGCTGAACCGCTCTTCCGATCT
CGATTA (język komputerowy)
Taatcg powiedział:
Starter do odwrotnego indeksowania PCR 11
CAAGCAGAAGACGGCATACGAGATTCACGGTCTCTCTCCTGCTGCTGAACCGCTCTTCCGATCT
CATTCA (Kotka)
tgaatg powiedział:
Podkład do odwróconego indeksowania PCR 12
CAAGCAGAAGACGGCATACGAGATGGAACTCGGTCTCTCCTCTGCTGCTGAACCGCTCTTCCGATCT
GGAAKT
agttcc powiedział:
Starter do odwrotnego indeksowania PCR 13
CAAGCAGAAGACGGCATACGAGATACATCGCGGTCTCTCCTGCTGCTGCTGAACCGCTCTTCCGATCT
ACATCG (Organizacja Współpracy Gospodarczej)
cgatgt powiedział:
Podkład do odwrotnego indeksowania PCR 14
CAAGCAGAAGACGGCATACGAGATAAGCTACGGTCTCTCCTGCTGCTGAACCGCTCTCTCTCTCTCTCT
AAGCTA (Agencja Kontroli Środowiska)
tagctt powiedział:
Podkład do odwrotnego indeksowania PCR 15
CAAGCAGAAGACGGCATACGAGATCAAGTTCGGTCTCTCCTGCTGCTGCTGCTGCTCTCTCTCTCTCTCTCTCTCTCTCTCTCT
CAAGTT (Biblioteka Główna)
Biblioteka AACTTG
Starter do odwróconego indeksowania PCR 16
CAAGCAGAAGACGGCATACGAGATGCCGGTCTCTCTCCTGCTGCTGCTGCTGCTCTTCCGATCT
GCCGGT (Zatoka Tysiąclecia)
Certyfikat
Podkład do odwrotnego indeksowania PCR 17
CAAGCAGAAGACGGCATACGAGATCGGTCTCTCTCTCTCCTGCTGCTGCTGCTGCTCTCTCTCTCT
CGGCCT powiedział:
aggccg powiedział:
Podkład do odwrotnego indeksowania PCR 18
CAAGCAGAAGACGGCATACGAGATTAGTTGCGGTCTCTCCTGCTGCTGCTGCTGCTCTCTCTCTCTCTCTCTCTCTCTCT
TAGTTG (Biblioteka TAGTTG)
Caacta (Grzybica)
Podkład do odwrotnego indeksowania PCR 19
CAAGCAGAAGACGGCATACGAGATGCGTGGTCTCTCTCCTGCTGCTGCTGAACCGCTCTTCCGATCT
Biblioteka GCGTGG
CCACGC powiedział:
Podkład do odwróconego indeksowania PCR 20
CAAGCAGAAGACGGCATACGAGATGTATAGCGGTCTCTCCTGCTGCTGCTGAACCGCTCTTCCGATCT
GTATAG (Znacznik GTATAG)
ctatac powiedział:
Podkład do odwrotnego indeksowania PCR 21
CAAGCAGAAGACGGCATACGAGATCCTTGCCGGTCTCTCTCTCCTGCTGCTGAACCGCTCTTCCGATCT
Konferencja CCTTGC (Międzynarodowa Organizacja Przemysł
gcaagg powiedział:
Starter do odwróconego indeksowania PCR 22
CAAGCAGAGACGGCATACGAGATGCTCTACGGTCTCTCTCCTGCTGCTGAACCGCTCTTCCGATCT
GCTGTA (Liga Mistrzowska)
tacagc powiedział:
Podkład do odwróconego indeksowania PCR 23
CAAGCAGAAGACGGCATACGAGATATGGCACGGTCTCTCCTGCTGCTGCTGAACCGCTCTTCCGATCT
ATGGCA (Międzynarodowa Sieć Operacyjna
tgccat powiedział:
Podkład do odwróconego indeksowania PCR 24
CAAGCAGAAGACGGCATACGAGATTGACATCGGTCTCTCCTCTGCTGCTGAACCGCTCTTCCGATCT
TGACAT (Dyrekcja Generalna)
ATGTCA
Tabela uzupełniająca 2: Sekwencje starterów z odwróconym indeksowaniem PCR.
Podkłady
Sekwencja
ZNF333_genic_H3K9me3_F
5'-AGCCTTCAATCAGCCATCATCCCT-3'
ZNF333_genic_H3K9me3_R
5'-TCTGGTATGGGTTCGCAGATGTGT-3'
HOXA9-10_F
5'-ACTGAAGTAATGAAGGCAGTGTCGT-3'
HOXA9-10_R
5'-GCAGCAYCAGAACTGGTCGGTG-3'
GAPDH_genic_H3K36me3_F
5'-AGGCAACTAGGATGGTGTGG-3'
GAPDH_genic_H3K36me3_R
5'-TTGATTTTGGAGGGATCTCG-3'
GAPDH-F
5'-TACTAGCGGTTTTACGGGCG-3'
GAPDH-R
5'-TCGAACAGGAGGAGCAGAGAGA-3'
Modyfikacja histonów
Pozytywny cel
Negatywny cel
H3K4me3 powiedział:
GAPDH
HOXA9-10 powiedział:
H3K4me1 powiedział:
GAPDH_genic
Zobacz materiał ZNF333
H3K27me3 powiedział:
HOXA9-10 powiedział:
Zobacz materiał ZNF333
H3K27ac powiedział:
GAPDH
Zobacz materiał ZNF333
H3K9me3 powiedział:
Zobacz materiał ZNF333
HOXA9-10 powiedział:
H3K36me3 powiedział:
GAPDH_genic
Zobacz materiał ZNF333
Tabela uzupełniająca 3: Lista starterów dla powszechnie używanych znaczników histonów (H3K4me3, H3K4me1, H3K27me3, H3K27ac, H3K9me3 i H3K36me3).
Plik uzupełniający 1: ndChIP-seq WorkSheet. Kliknij tutaj, aby pobrać ten plik.
Pliki kodu uzupełniającego: qcQpcr_v1.2. Pakiet oprogramowania statystycznego R do analizy wzbogacania qPCR natywnych bibliotek ChIP-seq o niskim nakładzie wejściowym. Kliknij tutaj, aby pobrać ten plik.